[Paraview] Enabling vtkSplineFilter in ParaView

2010-06-01 Thread Jérôme
Hi,

I was able to enable the vtkSplineFilter in ParaView-3.7 with the
attached ServerManager XML. But with the current 3.9, I am having a
segfault after the filter updates. Here is the WinXP32 call stack
after the segfault :

ntdll.dll!7c91120e()
[Les frames ci-dessous sont peut-être incorrects et/ou manquants,
aucun symbole chargé pour ntdll.dll]
msvcr90d.dll!102d528b() 
msvcr90d.dll!102d5ef9() 
msvcr90d.dll!102c103e() 
msvcr90d.dll!1023d496() 
msvcr90d.dll!102cff17() 
msvcr90d.dll!102cff0e() 
msvcr90d.dll!102cfb2f() 
msvcr90d.dll!102cfadc() 
msvcr90d.dll!102db25b() 
QtCored4.dll!qMalloc(unsigned int size=48)  Ligne 55 + 0xa octets   
C++
QtCored4.dll!QVectorData::allocate(int size=48, int alignment=4)
Ligne 67 + 0x28 octets  C++
QtGuid4.dll!QVectorQScriptItem::malloc(int aalloc=1696132182)
Ligne 391 + 0x23 octets C++
00e084c8()  
msvcr90d.dll!1023d496() 
msvcr90d.dll!102d3261() 
msvcr90d.dll!102d3258() 
QtGuid4.dll!QSpanData::adjustSpanMethods()  Ligne 5143 + 0xf octets 
C++
QtGuid4.dll!QSpanData::setup(const QBrush  brush={...}, int
alpha=1697441976, QPainter::CompositionMode compositionMode=284019584)
 Ligne 5095 C++
10edcb80()  
QtGuid4.dll!blend_color_argb(int count=-1, const QT_FT_Span_ *
spans=0x00e08ae8, void * userData=0x10edcb88)  Ligne 3082   C++
QtGuid4.dll!qt_span_fill_clipRect(int count=14715600, const
QT_FT_Span_ * spans=0x00e10b34, void * userData=0x652d84e8)  Ligne
4714 + 0x14 octets  C++
00e08ac8()  

Note that I also get this segfault in Fedora 12.

I certainly missed something in the XML that was not fatal with the
previous PV version. Please, could someone take a look at my XML, try
it and tell if this is reproducible?

Thanks a lot!
Jerome
ServerManagerConfiguration
   ProxyGroup name=filters
  SourceProxy name=SplineFilter
   class=vtkSplineFilter
  label=Spline Filter
 Documentation
  long_help= 
  short_help= 
 /Documentation
 InputProperty
   name=Input
   command=SetInputConnection
ProxyGroupDomain name=groups
   Group name=sources/
   Group name=filters/
/ProxyGroupDomain
DataTypeDomain name=input_type
   DataType value=vtkPolyData/
/DataTypeDomain
 /InputProperty
 
 IntVectorProperty
   name=Subdivide
   command=SetSubdivide
   number_of_elements=1
   default_values=1

EnumerationDomain name=enum
   Entry value=0 text=Specified/
   Entry value=1 text=Length/
/EnumerationDomain
 /IntVectorProperty
 IntVectorProperty
   name=NumberOfSubdivisions
   command=SetNumberOfSubdivisions
   number_of_elements=1
   default_values=5
 /IntVectorProperty
 
 DoubleVectorProperty
  name=Length
  command=SetLength
  number_of_elements=1
  default_values=1
 /DoubleVectorProperty
  /SourceProxy
   /ProxyGroup 
/ServerManagerConfiguration
___
Powered by www.kitware.com

Visit other Kitware open-source projects at 
http://www.kitware.com/opensource/opensource.html

Please keep messages on-topic and check the ParaView Wiki at: 
http://paraview.org/Wiki/ParaView

Follow this link to subscribe/unsubscribe:
http://www.paraview.org/mailman/listinfo/paraview


Re: [Paraview] Multi-Threading issues

2010-06-01 Thread liuning
Hi Biddiscombe,

I do a similar thing. But I use a Qt thread instead of a timer. The thread
controls the entire progress. I define a iteration steps in the thread,
e.g.,1024. In each iteration, the thread triggers the slot accept() method
(in ParaView3/Qt/Components/pqObjectInspectorWidget), in which the pipeline
is forced to update, then the source will call its RequestData() method and
invoke the simulation program to get the data. Maybe it is more appropriate
to let the simulation the control the entire progress, that is, when one
step simulation is done, then it will notify the ParaView to do the
corresponding rendering stuff. I just do it in the reverse way, because it
is easier to implement.



//the added thread class

pqSimulationThread::pqSimulationThread(pqObjectInspectorWidget* _gui){
this-gui=_gui;
connect(this,SIGNAL(oneStepSimulationFinished()),this-gui,SLOT(accept()));
}

pqSimulationThread::~pqSimulationThread(){
}

void pqSimulationThread::run(){

for(int i=0;iNUM_ITERATION;i++){
emit this-oneStepSimulationFinished();

//for synchronization
mutex.lock();
renderringFinishedCondatin.wait(mutex);
mutex.unlock();
}
}

===

// in the method accept() of pqObjectInspectorWidget

//force the pipeline to update
// i add a property to the source
//and any filters used will need this property too

source-pqProxy::setModifiedState(pqProxy::MODIFIED);
vtkSMSourceProxy * vp  = vtkSMSourceProxy::SafeDownCast(source-getProxy());
vp-InvokeCommand(FakePropertyForUpdate);

===

Then do some other little changes it will work with the simulation program.
I hope this can be of some help.

By the way, I hate the added property FakePropertyForUpdate , but I do not
know how to remove it . I have tried many other methods to force the
pipeline to update but got nothing good. And I also want to change the logic
to be more natural, that is, let the simulation program to control the whole
progress. I appreciate any help from anybody.

Thanks.

-Ning

On Tue, Jun 1, 2010 at 3:06 AM, Biddiscombe, John A. biddi...@cscs.chwrote:

 Utkarsh

  It is possible that the timer times out in middle of a progress event
  which happens while filters are executing.


 OK. This we have detected, though what progress is going on I'm unsure
 about (!) - by which I mean that the only thing paraview should be doing is
 updating our stuff and we block more updates whilst one is happening.

 Question : Is there a mutex of any kind on the update vtk objects, or
 update pipeline stuff which we can make use of to ensure that only one pull
 from the server occurs at any given time?

 JB
 ___
 Powered by www.kitware.com

 Visit other Kitware open-source projects at
 http://www.kitware.com/opensource/opensource.html

 Please keep messages on-topic and check the ParaView Wiki at:
 http://paraview.org/Wiki/ParaView

 Follow this link to subscribe/unsubscribe:
 http://www.paraview.org/mailman/listinfo/paraview

___
Powered by www.kitware.com

Visit other Kitware open-source projects at 
http://www.kitware.com/opensource/opensource.html

Please keep messages on-topic and check the ParaView Wiki at: 
http://paraview.org/Wiki/ParaView

Follow this link to subscribe/unsubscribe:
http://www.paraview.org/mailman/listinfo/paraview


Re: [Paraview] Enabling vtkSplineFilter in ParaView

2010-06-01 Thread Jérôme
Hi again !

I tried to replace the SplineFilter with a ProgrammableFilter
containing the following script :

input = self.GetInput()
splineFilter = vtk.vtkSplineFilter()
splineFilter.SetInput( input )
splineFilter.SetSubdivideToLength()
splineFilter.SetLength( 1 )
splineFilter.Update( )
self.GetOutput().DeepCopy(splineFilter.GetOutput())

This way, I am able to get a spline-based resampling. It means that my
XML description contains an error, but I cannot find it.

Jerome


2010/6/1 Jérôme jerome.ve...@gmail.com:
 Hi,

 I was able to enable the vtkSplineFilter in ParaView-3.7 with the
 attached ServerManager XML. But with the current 3.9, I am having a
 segfault after the filter updates. Here is the WinXP32 call stack
 after the segfault :

        ntdll.dll!7c91120e()
        [Les frames ci-dessous sont peut-être incorrects et/ou manquants,
 aucun symbole chargé pour ntdll.dll]
        msvcr90d.dll!102d528b()
        msvcr90d.dll!102d5ef9()
        msvcr90d.dll!102c103e()
        msvcr90d.dll!1023d496()
        msvcr90d.dll!102cff17()
        msvcr90d.dll!102cff0e()
        msvcr90d.dll!102cfb2f()
        msvcr90d.dll!102cfadc()
        msvcr90d.dll!102db25b()
        QtCored4.dll!qMalloc(unsigned int size=48)  Ligne 55 + 0xa octets      
  C++
        QtCored4.dll!QVectorData::allocate(int size=48, int alignment=4)
 Ligne 67 + 0x28 octets  C++
        QtGuid4.dll!QVectorQScriptItem::malloc(int aalloc=1696132182)
 Ligne 391 + 0x23 octets C++
        00e084c8()
        msvcr90d.dll!1023d496()
        msvcr90d.dll!102d3261()
        msvcr90d.dll!102d3258()
        QtGuid4.dll!QSpanData::adjustSpanMethods()  Ligne 5143 + 0xf octets    
  C++
        QtGuid4.dll!QSpanData::setup(const QBrush  brush={...}, int
 alpha=1697441976, QPainter::CompositionMode compositionMode=284019584)
  Ligne 5095     C++
        10edcb80()
        QtGuid4.dll!blend_color_argb(int count=-1, const QT_FT_Span_ *
 spans=0x00e08ae8, void * userData=0x10edcb88)  Ligne 3082       C++
        QtGuid4.dll!qt_span_fill_clipRect(int count=14715600, const
 QT_FT_Span_ * spans=0x00e10b34, void * userData=0x652d84e8)  Ligne
 4714 + 0x14 octets      C++
        00e08ac8()

 Note that I also get this segfault in Fedora 12.

 I certainly missed something in the XML that was not fatal with the
 previous PV version. Please, could someone take a look at my XML, try
 it and tell if this is reproducible?

 Thanks a lot!
 Jerome

___
Powered by www.kitware.com

Visit other Kitware open-source projects at 
http://www.kitware.com/opensource/opensource.html

Please keep messages on-topic and check the ParaView Wiki at: 
http://paraview.org/Wiki/ParaView

Follow this link to subscribe/unsubscribe:
http://www.paraview.org/mailman/listinfo/paraview


Re: [Paraview] [*] Re: [*] Build problem

2010-06-01 Thread rhushabh
Hi,
The error related to command prompt was removed.
But now during building I am getting 80% of the projects as failed with
many types of errors.

Please could you give me the exact link (other than the wiki page for
build and install) to build it from scratch again.
The specifications of my computer are:
Windows XP.

For building and installation required softwares in my computer:
Qt version 4.6.2/Qt SDK by Nokia v2010.02.1 (open source)
Cmake 2.6-patch 4
MPICH2
Python version 2.6.5
Microsoft Visual Studio 2005 Version 8.0.50727.42

I had many problems when configuring the source in Cmake mainly with what
parameters to change and to what value should i change them to.

Please help me out as I have wasted many hours in just building the
software on my system without any success.

One doubt: Pleae could you tell me that can I add programmble filters in
the .exe version of paraview 3.8.0.

Thank you in advance

Regards
Rhushabh Bhandari

 Yes, it's a visual studio issue. You might want to google around to
 try to find what's missing in your setup.

 Utkarsh

 On Mon, May 31, 2010 at 10:00 AM,  rhush...@cse.iitk.ac.in wrote:
 Hi,
 Thanks.
 I tried the link which you had sent, but I didnt have the problem which
 the person was having.
 Is the problem related to visual studio and not paraview as then I will
 try other versions of it.
 Thanks
 Regards
 Rhushabh Bhandari

 It doesn't look like it has anything to do with enabling/disabling
 Icet, although why would you want to do that? IceT is what ParaView
 uses for parallel rendering, without which parallel rendering won't
 work. Besides, IceT is enabled only when you are compiling with MPI.

 As far as /cmd /. /exe issue goes, this may help:
 http://www.interact-sw.co.uk/iangblog/2005/09/12/cmdspawnerror

 Utkarsh

 On Mon, May 31, 2010 at 6:00 AM,  rhush...@cse.iitk.ac.in wrote:
  Hi,
  While building paraview 3.8.0  by ALL_BUILD in visual studio,
  I get many errors stating:
  'Error spawning ((cmd)(.)(exe) command prompt)'
  Please could you help me in this situation

  And i have disabled ice-t, what are the effects of this on the
 software.
  Thanks in advance

  Regards,
  Rhushabh
 TEST MAIL

 ___
 Powered by www.kitware.com

 Visit other Kitware open-source projects at
 http://www.kitware.com/opensource/opensource.html

 Please keep messages on-topic and check the ParaView Wiki at:
 http://paraview.org/Wiki/ParaView

 Follow this link to subscribe/unsubscribe:
 http://www.paraview.org/mailman/listinfo/paraview









 ___
 Powered by www.kitware.com

 Visit other Kitware open-source projects at
 http://www.kitware.com/opensource/opensource.html

 Please keep messages on-topic and check the ParaView Wiki at:
 http://paraview.org/Wiki/ParaView

 Follow this link to subscribe/unsubscribe:
 http://www.paraview.org/mailman/listinfo/paraview








___
Powered by www.kitware.com

Visit other Kitware open-source projects at 
http://www.kitware.com/opensource/opensource.html

Please keep messages on-topic and check the ParaView Wiki at: 
http://paraview.org/Wiki/ParaView

Follow this link to subscribe/unsubscribe:
http://www.paraview.org/mailman/listinfo/paraview


Re: [Paraview] [*] Re: [*] Build problem

2010-06-01 Thread Utkarsh Ayachit
The Wiki page pretty much covers all that's necessary to build the
software. You'll have to be a bit more elaborate than many types of
errors before we can offer some concrete help. Are you using Qt
binaries downloaded from Nokia's website? As mentioned on their
download page, Nokia's binaries are built using MinGW and cannot be
used with Visual Studio builds. You need to build Qt from source using
visual studio.

Yes, you can add python programmable filters in the .exe version of ParaView.

Utkarsh

On Tue, Jun 1, 2010 at 9:47 AM,  rhush...@cse.iitk.ac.in wrote:
 Hi,
 The error related to command prompt was removed.
 But now during building I am getting 80% of the projects as failed with
 many types of errors.

 Please could you give me the exact link (other than the wiki page for
 build and install) to build it from scratch again.
 The specifications of my computer are:
 Windows XP.

 For building and installation required softwares in my computer:
 Qt version 4.6.2/Qt SDK by Nokia v2010.02.1 (open source)
 Cmake 2.6-patch 4
 MPICH2
 Python version 2.6.5
 Microsoft Visual Studio 2005 Version 8.0.50727.42

 I had many problems when configuring the source in Cmake mainly with what
 parameters to change and to what value should i change them to.

 Please help me out as I have wasted many hours in just building the
 software on my system without any success.

 One doubt: Pleae could you tell me that can I add programmble filters in
 the .exe version of paraview 3.8.0.

 Thank you in advance

 Regards
 Rhushabh Bhandari

 Yes, it's a visual studio issue. You might want to google around to
 try to find what's missing in your setup.

 Utkarsh

 On Mon, May 31, 2010 at 10:00 AM,  rhush...@cse.iitk.ac.in wrote:
 Hi,
 Thanks.
 I tried the link which you had sent, but I didnt have the problem which
 the person was having.
 Is the problem related to visual studio and not paraview as then I will
 try other versions of it.
 Thanks
 Regards
 Rhushabh Bhandari

 It doesn't look like it has anything to do with enabling/disabling
 Icet, although why would you want to do that? IceT is what ParaView
 uses for parallel rendering, without which parallel rendering won't
 work. Besides, IceT is enabled only when you are compiling with MPI.

 As far as /cmd /. /exe issue goes, this may help:
 http://www.interact-sw.co.uk/iangblog/2005/09/12/cmdspawnerror

 Utkarsh

 On Mon, May 31, 2010 at 6:00 AM,  rhush...@cse.iitk.ac.in wrote:
  Hi,
  While building paraview 3.8.0  by ALL_BUILD in visual studio,
  I get many errors stating:
  'Error spawning ((cmd)(.)(exe) command prompt)'
  Please could you help me in this situation

  And i have disabled ice-t, what are the effects of this on the
 software.
  Thanks in advance

  Regards,
  Rhushabh
 TEST MAIL

 ___
 Powered by www.kitware.com

 Visit other Kitware open-source projects at
 http://www.kitware.com/opensource/opensource.html

 Please keep messages on-topic and check the ParaView Wiki at:
 http://paraview.org/Wiki/ParaView

 Follow this link to subscribe/unsubscribe:
 http://www.paraview.org/mailman/listinfo/paraview









 ___
 Powered by www.kitware.com

 Visit other Kitware open-source projects at
 http://www.kitware.com/opensource/opensource.html

 Please keep messages on-topic and check the ParaView Wiki at:
 http://paraview.org/Wiki/ParaView

 Follow this link to subscribe/unsubscribe:
 http://www.paraview.org/mailman/listinfo/paraview









___
Powered by www.kitware.com

Visit other Kitware open-source projects at 
http://www.kitware.com/opensource/opensource.html

Please keep messages on-topic and check the ParaView Wiki at: 
http://paraview.org/Wiki/ParaView

Follow this link to subscribe/unsubscribe:
http://www.paraview.org/mailman/listinfo/paraview


[Paraview] Finite-size voxels and Point vs. Cell Centered Image Data

2010-06-01 Thread Eric Nodwell
What is the standard way to represent image data in VTK/ParaView, where the 
image data is intended to represent finite sized voxels?  Here the 3D image 
data represent values at the centre of voxels (or presumably some sort of 
weighted value centred on the voxel centre).

This may seem like a trivial question, however I am unable to obtain the 
behaviour that seems most appropriate.  For example, suppose I want to display 
all voxels in a (3D) image that exceed a certain threshold value.  The obvious 
thing to do seems to be to save the data in a VTK .vti file and then use the 
Theshold filter in ParaView.  The Threshold filter will display Cells, or not 
display Cells, depending on the threshold.  But Cells of a vtkDataObject are 
hexagons connecting the points, not hexagons centred on the points.  (There is 
an option All Scalars for the Threshold filter which determines how the cell 
value is calculated from the neighbouring point values, but it has no effect on 
the location of the cell.)

One possibility may be that we are storing the image data incorrectly in the 
.vti file.  Right now we have are generating the data as vtkImageSource with 
Scalar Point Data.  Perhaps it would be more logical to generate and save this 
as vtkImageSource with Scalar Cell Data, and indeed the ParaView documentation 
seems to imply this (see 
http://paraview.org/OnlineHelpCurrent/XMLImageDataReader.html ).  However, 
vtkImageSource does not seem to be intended to be used in this way.  In the 
vtkImageData implementation, all the methods that access scalar data do not 
check for the existence of Scalars for Cell Data, but instead always access the 
Point Data Scalars, like this:

  vtkDataArray *scalars = this-GetPointData()-GetScalars();

Or to present this another way, consider the trivial example of a single voxel. 
 We would like to have this represented as a finite size box, but the data is 
inherently a single value.  Here is some example VTK python code that attempts 
to generate a single voxel:

  image = vtk.vtkImageData()
  image.SetScalarTypeToInt()
  image.SetSpacing (1, 1, 1)
  image.SetOrigin (0, 0, 0)
  image.SetDimensions (1, 1, 1)
  image.SetScalarComponentFromFloat (0,0,0,0,10)
  print image.GetBounds()

The output for the bounds is:

(0.0, 0.0, 0.0, 0.0, 0.0, 0.0)

This is not necessarily incorrect in the right context.  One example of such a 
context might be the case where the image data represent discrete samples taken 
at distinct points in space.  However it does seem inappropriate for data meant 
to represent voxels as small volume elements (e.g. most medical images), where 
we would like a single voxel to have finite real-space bounds and be displayed 
as such.  I have again the feeling that vtkImageData is not the appropriate 
object to use to represent voxel data, but I don't know what the correct 
alternative would be (something with scalar cell data presumably).

As an aside, if you save the above image object with vtkXMLImageDataWriter, the 
file can't even be opened in ParaView 3.18.0, as the following error occurs:

ERROR: In 
/Users/eric/source/ParaView-3.8.0/Servers/Filters/vtkTexturePainter.cxx, line 
286
vtkTexturePainter (0x11de86350): Incorrect dimensionality.

However, a vtkImageData object with 2x2x2 dimensions can be opened and will 
display (using for example the Threshold filter) as a 1x1x1 box in ParaView.

Any insight would be appreciated.

Thanks,
Eric

___
Powered by www.kitware.com

Visit other Kitware open-source projects at 
http://www.kitware.com/opensource/opensource.html

Please keep messages on-topic and check the ParaView Wiki at: 
http://paraview.org/Wiki/ParaView

Follow this link to subscribe/unsubscribe:
http://www.paraview.org/mailman/listinfo/paraview


Re: [Paraview] Finite-size voxels and Point vs. Cell Centered Image Data

2010-06-01 Thread Francois Bertel
Hello,

in vtkImageData, Dimension is the number of points along each axis,
not the number of cells along each axis. If you want one voxel, then
you need 2x2x2 points.
Note that using GetDimension() is OK but SetDimension() is only here
for backward compatibility, use SetExtent() instead:

image-SetExtent(0,1,0,1,0,1);

On Tue, Jun 1, 2010 at 1:18 PM, Eric Nodwell enodw...@ucalgary.ca wrote:
 What is the standard way to represent image data in VTK/ParaView, where the 
 image data is intended to represent finite sized voxels?  Here the 3D image 
 data represent values at the centre of voxels (or presumably some sort of 
 weighted value centred on the voxel centre).

 This may seem like a trivial question, however I am unable to obtain the 
 behaviour that seems most appropriate.  For example, suppose I want to 
 display all voxels in a (3D) image that exceed a certain threshold value.  
 The obvious thing to do seems to be to save the data in a VTK .vti file and 
 then use the Theshold filter in ParaView.  The Threshold filter will display 
 Cells, or not display Cells, depending on the threshold.  But Cells of a 
 vtkDataObject are hexagons connecting the points, not hexagons centred on the 
 points.  (There is an option All Scalars for the Threshold filter which 
 determines how the cell value is calculated from the neighbouring point 
 values, but it has no effect on the location of the cell.)

 One possibility may be that we are storing the image data incorrectly in the 
 .vti file.  Right now we have are generating the data as vtkImageSource with 
 Scalar Point Data.  Perhaps it would be more logical to generate and save 
 this as vtkImageSource with Scalar Cell Data, and indeed the ParaView 
 documentation seems to imply this (see 
 http://paraview.org/OnlineHelpCurrent/XMLImageDataReader.html ).  However, 
 vtkImageSource does not seem to be intended to be used in this way.  In the 
 vtkImageData implementation, all the methods that access scalar data do not 
 check for the existence of Scalars for Cell Data, but instead always access 
 the Point Data Scalars, like this:

  vtkDataArray *scalars = this-GetPointData()-GetScalars();

 Or to present this another way, consider the trivial example of a single 
 voxel.  We would like to have this represented as a finite size box, but the 
 data is inherently a single value.  Here is some example VTK python code that 
 attempts to generate a single voxel:

  image = vtk.vtkImageData()
  image.SetScalarTypeToInt()
  image.SetSpacing (1, 1, 1)
  image.SetOrigin (0, 0, 0)
  image.SetDimensions (1, 1, 1)
  image.SetScalarComponentFromFloat (0,0,0,0,10)
  print image.GetBounds()

 The output for the bounds is:

 (0.0, 0.0, 0.0, 0.0, 0.0, 0.0)

 This is not necessarily incorrect in the right context.  One example of such 
 a context might be the case where the image data represent discrete samples 
 taken at distinct points in space.  However it does seem inappropriate for 
 data meant to represent voxels as small volume elements (e.g. most medical 
 images), where we would like a single voxel to have finite real-space bounds 
 and be displayed as such.  I have again the feeling that vtkImageData is not 
 the appropriate object to use to represent voxel data, but I don't know what 
 the correct alternative would be (something with scalar cell data presumably).

 As an aside, if you save the above image object with vtkXMLImageDataWriter, 
 the file can't even be opened in ParaView 3.18.0, as the following error 
 occurs:

 ERROR: In 
 /Users/eric/source/ParaView-3.8.0/Servers/Filters/vtkTexturePainter.cxx, line 
 286
 vtkTexturePainter (0x11de86350): Incorrect dimensionality.

 However, a vtkImageData object with 2x2x2 dimensions can be opened and will 
 display (using for example the Threshold filter) as a 1x1x1 box in ParaView.

 Any insight would be appreciated.

 Thanks,
 Eric

 ___
 Powered by www.kitware.com

 Visit other Kitware open-source projects at 
 http://www.kitware.com/opensource/opensource.html

 Please keep messages on-topic and check the ParaView Wiki at: 
 http://paraview.org/Wiki/ParaView

 Follow this link to subscribe/unsubscribe:
 http://www.paraview.org/mailman/listinfo/paraview




-- 
François Bertel, PhD  | Kitware Inc. Suite 204
1 (518) 371 3971 x113 | 28 Corporate Drive
  | Clifton Park NY 12065, USA
___
Powered by www.kitware.com

Visit other Kitware open-source projects at 
http://www.kitware.com/opensource/opensource.html

Please keep messages on-topic and check the ParaView Wiki at: 
http://paraview.org/Wiki/ParaView

Follow this link to subscribe/unsubscribe:
http://www.paraview.org/mailman/listinfo/paraview


Re: [Paraview] paraview - client-server

2010-06-01 Thread burlen

Hi,

I have seen this error mostly relating to a plugin being loaded on the 
client side but not the server side. It looks like 
vtk1DTransferFunctionFilter, is part of a PointSprite plugin. You could 
check the tools-plugin manager dialog to see if the client and server 
have the same plugins loaded?


It would also be helpful to know what type of file have you tried to open?

good luck
Burlen



Pierre-Olivier Dallaire wrote:

Hi,

I'm facing an error under the client-server mode. As soon as I try to open a 
file, paraview crashes with this error on the
server side :

ERROR: In /home/podallaire/ParaView/Servers/Common/vtkProcessModule.cxx, line 
1065
vtkProcessModule (0x851da0): Cannot create object of type 
vtk1DTransferFunctionFilter.
while processing
Message 0 = New
  Argument 0 = string_value {vtk1DTransferFunctionFilter}
  Argument 1 = id_value {387}


ERROR: In /home/podallaire/ParaView/Servers/Common/vtkProcessModule.cxx, line 
1066
vtkProcessModule (0x851da0): Aborting execution for debugging purposes.




Any idea what is going on ? I tried different files - same result.

Regards,

PO

On 2010-04-30, at 2:49 PM, burlen wrote:

  

yeah, that's a typo in the bug report. Both last night and this morning I have 
been unable to access mantis or the the pv web site. Page never loads. I'll try 
later.

pat marion wrote:


Hey Burlen, on the bug report page for 10283, I think you need to fix the 
command line you are testing with :

$ ssh remote cmd1  cm2

will execute cmd1 on remote and cmd2 locally.  It should be:

$ ssh remote cmd1  cmd2

Pat

On Fri, Apr 30, 2010 at 9:12 AM, pat marion pat.mar...@kitware.com 
mailto:pat.mar...@kitware.com wrote:

   I have applied your patch.  I agree that paraview should explicity
   close the child process.  But... what I am pointing out is that
   calling QProcess::close() does not help in this situation.  What I
   am saying is that, even when paraview does kill the process, any
   commands run by ssh on the other side of the netpipe will be
   orphaned by sshd.  Are you sure you can't reproduce it?


   $ ssh localhost sleep 1d
   $  press control-c 
   $ pidof sleep
   $ # sleep is still running

   Pat


   On Fri, Apr 30, 2010 at 2:08 AM, burlen burlen.lor...@gmail.com
   mailto:burlen.lor...@gmail.com wrote:

   Hi Pat,

   From my point of view the issue is philosophical, because
   practically speaking I couldn't reproduce the orphans with out
   doing something a little odd namely, ssh ...   sleep 1d.
   Although the fact that a user reported suggests that it may
   occur in the real world as well. The question is this: should
   an application explicitly clean up resources it allocates? or
   should an application rely on the user not only knowing that
   there is the potential for a resource leak but also knowing
   enough to do the right thing to avoid it (eg ssh -tt ...)? In
   my opinion, as a matter of principle, if PV spawns a process
   it should explicitly clean it up and there should be no way it
   can become an orphan. In this case the fact that the orphan
   can hold ports open is particularly insidious, because further
   connection attempt on that port fails with no helpful error
   information. Also it is not very difficult to clean up a
   spawned process. What it comes down to is a little book
   keeping to hang on to the qprocess handle and a few lines of
   code called from pqCommandServerStartup destructor to make
   certain it's cleaned up. This is from the patch I submitted
   when I filed the bug report.

   +// close running process
   +if (this-Process-state()==QProcess::Running)
   +  {
   +  this-Process-close();
   +  }
   +// free the object
   +delete this-Process;
   +this-Process=NULL;

   I think if the cluster admins out there new which ssh options
   (GatewayPorts etc) are important for ParView to work
   seamlessly, then they might be willing to open them up. It's
   my impression that the folks that build clusters want tools
   like PV to be easy to use, but they don't necessarily know all
   the in's and out's of confinguring and running PV.

   Thanks for looking at this again! The -tt option to ssh is
   indeed a good find.

   Burlen

   pat marion wrote:

   Hi all!

   I'm bringing this thread back- I have learned a couple new
   things...

   ---
   No more orphans:

   Here is an easy way to create an orphan:

 $ ssh localhost sleep 1d
 $ press control c

   The ssh process is cleaned up, but sshd orphans the sleep
   process.  You can avoid this by adding '-t' to ssh:

$ ssh -t localhost sleep 1d

   Works like a charm!  But then there is another problem...
   try this command 

Re: [Paraview] extract cell volume with python

2010-06-01 Thread Fabian Braennstroem

Hello to all,

thanks for your help. Oh, I think I got it...

Regards!
Fabian

On 05/30/2010 02:36 PM, Andy Bauer wrote:

There are a couple of gradient filters in ParaView (Gradient and
Gradient of Unstructured DataSet in the filters menu).  The one to use
depends on the input grid you have.  The Mesh Quality filter can be used
to get the volume of tet and hex cells.

Andy

On Sun, May 30, 2010 at 5:49 AM, Fabian Braennstroem
f.braennstr...@gmx.de mailto:f.braennstr...@gmx.de wrote:

Hello,

sorry, but did I miss anything simple...
Would be nice, to get some hints!

Regards!
Fabian

On 05/13/2010 02:05 PM, Fabian Braennstroem wrote:

Hi,

I would like to extract the cell volume and gradients using python
programmable filter.
Does anyone have a hint, how to do this?

___
Powered by www.kitware.com http://www.kitware.com

Visit other Kitware open-source projects at
http://www.kitware.com/opensource/opensource.html

Please keep messages on-topic and check the ParaView Wiki at:
http://paraview.org/Wiki/ParaView


Follow this link to subscribe/unsubscribe:
http://www.paraview.org/mailman/listinfo/paraview




___
Powered by www.kitware.com

Visit other Kitware open-source projects at 
http://www.kitware.com/opensource/opensource.html

Please keep messages on-topic and check the ParaView Wiki at: 
http://paraview.org/Wiki/ParaView

Follow this link to subscribe/unsubscribe:
http://www.paraview.org/mailman/listinfo/paraview


Re: [Paraview] vtkCellLocator: Python

2010-06-01 Thread bastil2...@yahoo.de
 Thanks Eric, all,

but this does not seem to work at all. My loop looks something like this:

input = self.GetInputDataObject(0, 0)
output = self.GetOutputDataObject(0)
if input.IsA(vtkMultiBlockDataSet):
iter = input.NewIterator()
iter.UnRegister(None)
iter.InitTraversal()
i=1
while not iter.IsDoneWithTraversal():
curInput = iter.GetCurrentDataObject()
curInformation = curInput.GetInformation()
print i
print iter.GetCurrentMetaData()
# This is where I want to query for the current name, line does
not work
print iter.GetCurrentMetaData().Get(iter.GetCurrentMetaData.NAME)
i=i+1
iter.GoToNextItem();

The output of print iter.GetCurrentMetaData() looks like this:

vtkInformation (0x3808f00)
Debug: Off
Modified Time: 269417
Reference Count: 2
Registered Events: (none)
DONT_PRUNE: 1
NAME: my_given_name

I want to work wit the NAME attribute ('my_given_name in this case)
but I don't know how...?

bastil2...@yahoo.de wrote:


  Original-Nachricht 
 Betreff:  Re: [Paraview] vtkCellLocator: Python
 Datum:Mon, 31 May 2010 15:15:54 -0400
 Von:  Eric E. Monson emon...@cs.duke.edu
 An:   bastil2...@yahoo.de
 CC:   Utkarsh Ayachit utkarsh.ayac...@kitware.com, 'ParaView'
 paraview@paraview.org



 Hey Bastian,

 I think you should be able to just use the Get() method to pull out
 that information key. When I want to find current time step
 information on some pipeline data (with a programmable filter), which
 has information that looks like this:

 vtkInformation (0x11dd723d0)
 Debug: Off
 Modified Time: 70754
 Reference Count: 2
 Registered Events: (none)
 DATA_NUMBER_OF_GHOST_LEVELS: 0
 DATA_TIME_STEPS: 0
 DATA_PIECE_NUMBER: 0
 DATA_EXTENT_TYPE: 0
 DATA_NUMBER_OF_PIECES: 1

 I can do it with these commands:

 pdi = self.GetInputDataObject(0,0)
 current_time = pdi.GetInformation().Get(pdi.DATA_TIME_STEPS(),0)

 So you might be able to just put NAME() in where I put DATA_TIME_STEPS().

 Talk to you later,
 -Eric

 --
 Eric E Monson
 Duke Visualization Technology Group


 On May 31, 2010, at 1:23 PM, bastil2...@yahoo.de
 mailto:bastil2...@yahoo.de wrote:

 Thanks Utkarsh,

 the required information is in MetaData:

 print iter.GetCurrentMetaData()

 vtkInformation (0x3808f00)
 Debug: Off
 Modified Time: 269417
 Reference Count: 2
 Registered Events: (none)
 DONT_PRUNE: 1
 NAME: my_given_name

 Now I want to work with the NAME-Attribute (my_given_name) but I
 don't know how to get this into a variable? I can't find a
 GetName()-Method...

 print dir(iter.GetCurrentMetaData())

 ['AddObserver', 'Append', 'AppendUnique', 'BreakOnError', 'Clear',
 'Copy', 'CopyEntries', 'CopyEntry', 'DebugOff', 'DebugOn',
 'FastDelete', 'Get', 'GetAddressAsString', 'GetClassName',
 'GetDebug', 'GetGlobalWarningDisplay', 'GetKey', 'GetMTime',
 'GetNumberOfKeys', 'GetReferenceCount', 'GetRequest',
 'GlobalWarningDisplayOff', 'GlobalWarningDisplayOn', 'Has',
 'HasObserver', 'InvokeEvent', 'IsA', 'IsTypeOf', 'Length',
 'Modified', 'NewInstance', 'PrintRevisions', 'Register', 'Remove',
 'RemoveAllObservers', 'RemoveObserver', 'RemoveObservers',
 'SafeDownCast', 'Set', 'SetDebug', 'SetGlobalWarningDisplay',
 'SetReferenceCount', 'SetRequest', 'UnRegister']

 Thanks for help.

 Regards Bastian

 Am 30.05.2010 16:58, schrieb Utkarsh Ayachit:
 The following example shows how to create an iterator to iterate over
 all blocks in a composite dataset. The iterator
 (http://www.vtk.org/doc/nightly/html/classvtkCompositeDataIterator.html)
 has methods like HasCurrentMetaData(), GetCurrentMetaData() that you
 can use to get the meta-data for the current block. Note that by
 default the iterator skips over non-leaf nodes. If you want to access
 those, use VisitOnlyLeavesOff() on the iterator.

 http://www.itk.org/Wiki/Python_Programmable_Filter#Dealing_with_Composite_Datasets

 Utkarsh

 On Sat, May 29, 2010 at 1:15 PM, BastiL2001 bastil2...@yahoo.de wrote:
   
 Thanks Andi,

 I don't know since the ensight reader does this for me. I want to work with
 the input he produces.

 Is it in any way possible to link GetMetaData to the current block? Afaik I
 can not query each blocks meta-Data, can I?

 Are you naming the block in the meta-data?  In vtkCompositeDataSet there is
 a method called GetMetaData
 (http://www.vtk.org/doc/nightly/html/classvtkCompositeDataSet.html#abcf79d23fc3a67138da3336e4cda2afe)
 that may be what you're looking for.

 Andy

 On Fri, May 28, 2010 at 3:58 PM, BastiL2001 bastil2...@yahoo.de wrote:
 
 Thanks Andi,

 I will do this. What I am planning to do is to create a scalr field with
 some distance information as an output.

 Another issue: I want to loop oder a multi-block dataset with an iterator
 in a programmable filter (works) and find out the name of the currently
 visited block. Is that possible? Thanks.

 Regards Bastian

 
 

Re: [Paraview] Finite-size voxels and Point vs. Cell Centered Image Data

2010-06-01 Thread Eric Nodwell
Francois,

Thank you for the reply.

Using 2x2x2 points for one voxel is fine in principle.  The difficulty is, you 
can't really conveniently (as far as I can tell) use a vtkImageData object for 
this purpose.  Because now of course our scalar image data needs to be 
associated with cells, not to points.  But vtkImageData in many cases only 
deals with the scalar data attached to the PointData.  For example, this is 
from vtkImageData.cxx :

//
// This method returns a pointer to the origin of the vtkImageData.
void *vtkImageData::GetScalarPointer()
{
  if (this-PointData-GetScalars() == NULL)
{
vtkDebugMacro(Allocating scalars in ImageData);
this-AllocateScalars();
}
  return this-PointData-GetScalars()-GetVoidPointer(0);
}

Notice that is method just grabs the Scalars from PointData, and other methods 
of vtkImageData are similar, so if you want to store the scalars in CellData, 
which is necessary represent voxels, it is quite difficult.

Eric


 Date: Tue, 1 Jun 2010 13:26:15 -0400
 From: Francois Bertel francois.ber...@kitware.com
 Subject: Re: [Paraview] Finite-size voxels and Point vs. Cell Centered
   Image   Data
 To: paraview@paraview.org
 Message-ID:
   aanlktikz6uw4tofeudmdu0xocwc_1b4_0rs9n_z9s...@mail.gmail.com
 Content-Type: text/plain; charset=ISO-8859-1
 
 Hello,
 
 in vtkImageData, Dimension is the number of points along each axis,
 not the number of cells along each axis. If you want one voxel, then
 you need 2x2x2 points.
 Note that using GetDimension() is OK but SetDimension() is only here
 for backward compatibility, use SetExtent() instead:
 
 image-SetExtent(0,1,0,1,0,1);

___
Powered by www.kitware.com

Visit other Kitware open-source projects at 
http://www.kitware.com/opensource/opensource.html

Please keep messages on-topic and check the ParaView Wiki at: 
http://paraview.org/Wiki/ParaView

Follow this link to subscribe/unsubscribe:
http://www.paraview.org/mailman/listinfo/paraview


Re: [Paraview] Finite-size voxels and Point vs. Cell Centered Image Data

2010-06-01 Thread burlen

Hi Eric,

Because now of course our scalar image data needs to be associated with cells, 
not to points.  But vtkImageData in many cases only deals with the scalar data 
attached to the PointData.
Could you use the concept of a dual grid? In other words make your vis 
grid node centered with the nodes located on the cell centers of your 
computational grid? Your arrays end up in PointData that way with no 
interpolations required.


Burlen

Eric Nodwell wrote:

Francois,

Thank you for the reply.

Using 2x2x2 points for one voxel is fine in principle.  The difficulty is, you 
can't really conveniently (as far as I can tell) use a vtkImageData object for 
this purpose.  Because now of course our scalar image data needs to be 
associated with cells, not to points.  But vtkImageData in many cases only 
deals with the scalar data attached to the PointData.  For example, this is 
from vtkImageData.cxx :

//
// This method returns a pointer to the origin of the vtkImageData.
void *vtkImageData::GetScalarPointer()
{
  if (this-PointData-GetScalars() == NULL)
{
vtkDebugMacro(Allocating scalars in ImageData);
this-AllocateScalars();
}
  return this-PointData-GetScalars()-GetVoidPointer(0);
}

Notice that is method just grabs the Scalars from PointData, and other methods 
of vtkImageData are similar, so if you want to store the scalars in CellData, 
which is necessary represent voxels, it is quite difficult.

Eric


  

Date: Tue, 1 Jun 2010 13:26:15 -0400
From: Francois Bertel francois.ber...@kitware.com
Subject: Re: [Paraview] Finite-size voxels and Point vs. Cell Centered
Image   Data
To: paraview@paraview.org
Message-ID:
aanlktikz6uw4tofeudmdu0xocwc_1b4_0rs9n_z9s...@mail.gmail.com
Content-Type: text/plain; charset=ISO-8859-1

Hello,

in vtkImageData, Dimension is the number of points along each axis,
not the number of cells along each axis. If you want one voxel, then
you need 2x2x2 points.
Note that using GetDimension() is OK but SetDimension() is only here
for backward compatibility, use SetExtent() instead:

image-SetExtent(0,1,0,1,0,1);



___
Powered by www.kitware.com

Visit other Kitware open-source projects at 
http://www.kitware.com/opensource/opensource.html

Please keep messages on-topic and check the ParaView Wiki at: 
http://paraview.org/Wiki/ParaView

Follow this link to subscribe/unsubscribe:
http://www.paraview.org/mailman/listinfo/paraview
  


___
Powered by www.kitware.com

Visit other Kitware open-source projects at 
http://www.kitware.com/opensource/opensource.html

Please keep messages on-topic and check the ParaView Wiki at: 
http://paraview.org/Wiki/ParaView

Follow this link to subscribe/unsubscribe:
http://www.paraview.org/mailman/listinfo/paraview


Re: [Paraview] vtkCellLocator: Python

2010-06-01 Thread Utkarsh Ayachit
Try the following:

 print iter.GetCurrentMetaData().Get(vtk.vtkCompositeDataSet.NAME())

Utkarsh

On Tue, Jun 1, 2010 at 3:12 PM, bastil2...@yahoo.de bastil2...@yahoo.dewrote:

   Thanks Eric, all,

 but this does not seem to work at all. My loop looks something like this:

 input = self.GetInputDataObject(0, 0)
 output = self.GetOutputDataObject(0)
 if input.IsA(vtkMultiBlockDataSet):
 iter = input.NewIterator()
 iter.UnRegister(None)
 iter.InitTraversal()
 i=1
 while not iter.IsDoneWithTraversal():
 curInput = iter.GetCurrentDataObject()
 curInformation = curInput.GetInformation()
 print i
 print iter.GetCurrentMetaData()
 # This is where I want to query for the current name, line does not
 work
 print iter.GetCurrentMetaData().Get(iter.GetCurrentMetaData.NAME)
 i=i+1
 iter.GoToNextItem();

 The output of print iter.GetCurrentMetaData() looks like this:


  vtkInformation (0x3808f00)
  Debug: Off
  Modified Time: 269417
  Reference Count: 2
  Registered Events: (none)
  DONT_PRUNE: 1
  NAME: my_given_name

  I want to work wit the NAME attribute ('my_given_name in this case) but
 I don't know how...?

 bastil2...@yahoo.de wrote:



  Original-Nachricht   Betreff: Re: [Paraview]
 vtkCellLocator: Python  Datum: Mon, 31 May 2010 15:15:54 -0400  Von: Eric
 E. Monson emon...@cs.duke.edu emon...@cs.duke.edu  An:
 bastil2...@yahoo.de  CC: Utkarsh Ayachit 
 utkarsh.ayac...@kitware.comutkarsh.ayac...@kitware.com,
 'ParaView' paraview@paraview.org paraview@paraview.org


 Hey Bastian,

  I think you should be able to just use the Get() method to pull out that
 information key. When I want to find current time step information on some
 pipeline data (with a programmable filter), which has information that looks
 like this:

  vtkInformation (0x11dd723d0)
 Debug: Off
 Modified Time: 70754
 Reference Count: 2
 Registered Events: (none)
 DATA_NUMBER_OF_GHOST_LEVELS: 0
 DATA_TIME_STEPS: 0
 DATA_PIECE_NUMBER: 0
 DATA_EXTENT_TYPE: 0
 DATA_NUMBER_OF_PIECES: 1

  I can do it with these commands:

  pdi = self.GetInputDataObject(0,0)
 current_time = pdi.GetInformation().Get(pdi.DATA_TIME_STEPS(),0)

  So you might be able to just put NAME() in where I put DATA_TIME_STEPS().

  Talk to you later,
 -Eric

  --
 Eric E Monson
 Duke Visualization Technology Group


  On May 31, 2010, at 1:23 PM, bastil2...@yahoo.de wrote:

  Thanks Utkarsh,

 the required information is in MetaData:

 print iter.GetCurrentMetaData()

 vtkInformation (0x3808f00)
 Debug: Off
 Modified Time: 269417
 Reference Count: 2
 Registered Events: (none)
 DONT_PRUNE: 1
 NAME: my_given_name

 Now I want to work with the NAME-Attribute (my_given_name) but I don't know
 how to get this into a variable? I can't find a GetName()-Method...

 print dir(iter.GetCurrentMetaData())

 ['AddObserver', 'Append', 'AppendUnique', 'BreakOnError', 'Clear', 'Copy',
 'CopyEntries', 'CopyEntry', 'DebugOff', 'DebugOn', 'FastDelete', 'Get',
 'GetAddressAsString', 'GetClassName', 'GetDebug', 'GetGlobalWarningDisplay',
 'GetKey', 'GetMTime', 'GetNumberOfKeys', 'GetReferenceCount', 'GetRequest',
 'GlobalWarningDisplayOff', 'GlobalWarningDisplayOn', 'Has', 'HasObserver',
 'InvokeEvent', 'IsA', 'IsTypeOf', 'Length', 'Modified', 'NewInstance',
 'PrintRevisions', 'Register', 'Remove', 'RemoveAllObservers',
 'RemoveObserver', 'RemoveObservers', 'SafeDownCast', 'Set', 'SetDebug',
 'SetGlobalWarningDisplay', 'SetReferenceCount', 'SetRequest', 'UnRegister']

 Thanks for help.

 Regards Bastian

 Am 30.05.2010 16:58, schrieb Utkarsh Ayachit:

 The following example shows how to create an iterator to iterate over
 all blocks in a composite dataset. The iterator
 (http://www.vtk.org/doc/nightly/html/classvtkCompositeDataIterator.html)
 has methods like HasCurrentMetaData(), GetCurrentMetaData() that you
 can use to get the meta-data for the current block. Note that by
 default the iterator skips over non-leaf nodes. If you want to access
 those, use VisitOnlyLeavesOff() on the iterator.
 http://www.itk.org/Wiki/Python_Programmable_Filter#Dealing_with_Composite_Datasets

 Utkarsh

 On Sat, May 29, 2010 at 1:15 PM, BastiL2001 bastil2...@yahoo.de 
 bastil2...@yahoo.de wrote:


  Thanks Andi,

 I don't know since the ensight reader does this for me. I want to work with
 the input he produces.

 Is it in any way possible to link GetMetaData to the current block? Afaik I
 can not query each blocks meta-Data, can I?

 Are you naming the block in the meta-data?  In vtkCompositeDataSet there is
 a method called GetMetaData
 (http://www.vtk.org/doc/nightly/html/classvtkCompositeDataSet.html#abcf79d23fc3a67138da3336e4cda2afe)
 that may be what you're looking for.

 Andy

 On Fri, May 28, 2010 at 3:58 PM, BastiL2001 bastil2...@yahoo.de 
 bastil2...@yahoo.de wrote:


  Thanks Andi,

 I will do this. What I am planning to do is to create a scalr field 

Re: [Paraview] Finite-size voxels and Point vs. Cell Centered Image Data

2010-06-01 Thread Eric Nodwell
Burlen,

I assume by dual grid you mean creating two vtkImageData objects, which both 
contain pointers (properly registered and reference counted) to the same 
vtkDataArray object, but in one case, as the Scalars for the CellData, and in 
the other case as the Scalars for the PointsData.  This might work, although it 
does seem unnecessarily complicated for storing voxels as volume elements 
which would seem to be a common usage.  And I'm not sure it answers all the 
questions such as, how would you write the data to a file?  If you want it 
accurately rendered in ParaView, it needs to be stored with CellData Scalars, 
but that could be problematic if you (or someone else) wants to read it and 
process it with VTK again.

Also there seems to be no guarantee that vtkImageData won't clobber your 
CellData Scalars when you're not looking.  Consider the following method copied 
from vtkImageData.cxx .  I have no idea when this might be called, but notice 
that PointData Scalars will survive it, while CellData Scalars won't.

//
// Graphics filters reallocate every execute.  Image filters try to reuse
// the scalars.
void vtkImageData::PrepareForNewData()
{
  // free everything but the scalars
  vtkDataArray *scalars = this-GetPointData()-GetScalars();
  if (scalars)
{
scalars-Register(this);
}
  this-Initialize();
  if (scalars)
{
this-GetPointData()-SetScalars(scalars);
scalars-UnRegister(this);
}
}



On 2010-06-01, at 2:04 PM, burlen wrote:

 Hi Eric,
 Because now of course our scalar image data needs to be associated with 
 cells, not to points.  But vtkImageData in many cases only deals with the 
 scalar data attached to the PointData.
 Could you use the concept of a dual grid? In other words make your vis grid 
 node centered with the nodes located on the cell centers of your 
 computational grid? Your arrays end up in PointData that way with no 
 interpolations required.
 
 Burlen
 
 Eric Nodwell wrote:
 Francois,
 
 Thank you for the reply.
 
 Using 2x2x2 points for one voxel is fine in principle.  The difficulty is, 
 you can't really conveniently (as far as I can tell) use a vtkImageData 
 object for this purpose.  Because now of course our scalar image data needs 
 to be associated with cells, not to points.  But vtkImageData in many cases 
 only deals with the scalar data attached to the PointData.  For example, 
 this is from vtkImageData.cxx :
 
 //
 // This method returns a pointer to the origin of the vtkImageData.
 void *vtkImageData::GetScalarPointer()
 {
  if (this-PointData-GetScalars() == NULL)
{
vtkDebugMacro(Allocating scalars in ImageData);
this-AllocateScalars();
}
  return this-PointData-GetScalars()-GetVoidPointer(0);
 }
 
 Notice that is method just grabs the Scalars from PointData, and other 
 methods of vtkImageData are similar, so if you want to store the scalars in 
 CellData, which is necessary represent voxels, it is quite difficult.
 
 Eric
 
 
  
 Date: Tue, 1 Jun 2010 13:26:15 -0400
 From: Francois Bertel francois.ber...@kitware.com
 Subject: Re: [Paraview] Finite-size voxels and Point vs. Cell Centered
 Image   Data
 To: paraview@paraview.org
 Message-ID:
 aanlktikz6uw4tofeudmdu0xocwc_1b4_0rs9n_z9s...@mail.gmail.com
 Content-Type: text/plain; charset=ISO-8859-1
 
 Hello,
 
 in vtkImageData, Dimension is the number of points along each axis,
 not the number of cells along each axis. If you want one voxel, then
 you need 2x2x2 points.
 Note that using GetDimension() is OK but SetDimension() is only here
 for backward compatibility, use SetExtent() instead:
 
 image-SetExtent(0,1,0,1,0,1);

 
 ___
 Powered by www.kitware.com
 
 Visit other Kitware open-source projects at 
 http://www.kitware.com/opensource/opensource.html
 
 Please keep messages on-topic and check the ParaView Wiki at: 
 http://paraview.org/Wiki/ParaView
 
 Follow this link to subscribe/unsubscribe:
 http://www.paraview.org/mailman/listinfo/paraview
  
 
 

___
Powered by www.kitware.com

Visit other Kitware open-source projects at 
http://www.kitware.com/opensource/opensource.html

Please keep messages on-topic and check the ParaView Wiki at: 
http://paraview.org/Wiki/ParaView

Follow this link to subscribe/unsubscribe:
http://www.paraview.org/mailman/listinfo/paraview


[Paraview] Compiling error with 3.8.0

2010-06-01 Thread Fred Fred

Hello,
I have downloaded the 3.8 release source from the download page and I try to 
compile it on Mac OS X Leopard, using the classic cmake porcedure that has been 
successful with the 3.6.2 but I get a surprising error:

[ 34%] Building CXX object 
VTK/VolumeRendering/CMakeFiles/vtkVolumeRenderingPython.dir/vtkVolumeRenderingPythonInit.cxx.o
Linking CXX shared module ../../bin/libvtkVolumeRenderingPython.so
[ 34%] Built target vtkVolumeRenderingPython
[ 34%] Generating vtkCompositeZPassShader_fs.cxx, vtkCompositeZPassShader_fs.h
[ 34%] Generating vtkSocketCommunicatorHash.h
dyld: Library not loaded: libvtksys.pv3.8.dylib
  Referenced from: 
/usr/local/ParaView3_BUILD/VTK/Parallel/../../bin/vtkHashSource
  Reason: image not found
/bin/sh: line 1:   970 Trace/BPT trap  ../../bin/vtkHashSource 
/usr/local/ParaView-3.8.0/VTK/Parallel/vtkSocketCommunicator.cxx 
vtkSocketCommunicatorHash vtkSocketCommunicatorHash.h
make[2]: *** [VTK/Parallel/vtkSocketCommunicatorHash.h] Error 133
make[1]: *** [VTK/Parallel/CMakeFiles/vtkParallel.dir/all] Error 2
make: *** [all] Error 2
Macintosh-3:ParaView3_BUILD$ find . -name *libvtksys*
./bin/libvtksys.dylib
./bin/libvtksys.pv3.8.dylib

So the required library is available in ./bin, which should be included in the 
path to search for librairies, but it does not work.
What should I do to solve this issue?


  
_
Installez gratuitement les nouvelles Emoch'ticones !
http://www.ilovemessenger.fr/emoticones/telecharger-emoticones-emochticones.aspx___
Powered by www.kitware.com

Visit other Kitware open-source projects at 
http://www.kitware.com/opensource/opensource.html

Please keep messages on-topic and check the ParaView Wiki at: 
http://paraview.org/Wiki/ParaView

Follow this link to subscribe/unsubscribe:
http://www.paraview.org/mailman/listinfo/paraview


Re: [Paraview] Compiling error with 3.8.0

2010-06-01 Thread Dave Partyka
Turn VTK_USE_RPATH to ON or export DYLD_LIBRARY_PATH to your bin dir.

On Tue, Jun 1, 2010 at 5:07 PM, Fred Fred stan1...@hotmail.fr wrote:

  Hello,
 I have downloaded the 3.8 release source from the download page and I try
 to compile it on Mac OS X Leopard, using the classic cmake porcedure that
 has been successful with the 3.6.2 but I get a surprising error:

 [ 34%] Building CXX object
 VTK/VolumeRendering/CMakeFiles/vtkVolumeRenderingPython.dir/vtkVolumeRenderingPythonInit.cxx.o
 Linking CXX shared module ../../bin/libvtkVolumeRenderingPython.so
 [ 34%] Built target vtkVolumeRenderingPython
 [ 34%] Generating vtkCompositeZPassShader_fs.cxx,
 vtkCompositeZPassShader_fs.h
 [ 34%] Generating vtkSocketCommunicatorHash.h
 dyld: Library not loaded: libvtksys.pv3.8.dylib
   Referenced from:
 /usr/local/ParaView3_BUILD/VTK/Parallel/../../bin/vtkHashSource
   Reason: image not found
 /bin/sh: line 1:   970 Trace/BPT trap  ../../bin/vtkHashSource
 /usr/local/ParaView-3.8.0/VTK/Parallel/vtkSocketCommunicator.cxx
 vtkSocketCommunicatorHash vtkSocketCommunicatorHash.h
 make[2]: *** [VTK/Parallel/vtkSocketCommunicatorHash.h] Error 133
 make[1]: *** [VTK/Parallel/CMakeFiles/vtkParallel.dir/all] Error 2
 make: *** [all] Error 2
 Macintosh-3:ParaView3_BUILD$ find . -name *libvtksys*
 ./bin/libvtksys.dylib
 ./bin/libvtksys.pv3.8.dylib

 So the required library is available in ./bin, which should be included in
 the path to search for librairies, but it does not work.
 What should I do to solve this issue?



 --
 Vous voulez protéger votre vie privée ? La solution avec Internet Explorer
 8 http://clk.atdmt.com/FRM/go/232102477/direct/01/

 ___
 Powered by www.kitware.com

 Visit other Kitware open-source projects at
 http://www.kitware.com/opensource/opensource.html

 Please keep messages on-topic and check the ParaView Wiki at:
 http://paraview.org/Wiki/ParaView

 Follow this link to subscribe/unsubscribe:
 http://www.paraview.org/mailman/listinfo/paraview


___
Powered by www.kitware.com

Visit other Kitware open-source projects at 
http://www.kitware.com/opensource/opensource.html

Please keep messages on-topic and check the ParaView Wiki at: 
http://paraview.org/Wiki/ParaView

Follow this link to subscribe/unsubscribe:
http://www.paraview.org/mailman/listinfo/paraview


[Paraview] Gzipped image file load

2010-06-01 Thread Joohwi Lee
Hi,

I tried to load gzipped NRRD file in Paraview, but it can't recognize the
image saying that

  ERROR: In
/home/kitware/Dashboard/MyTests/ParaView-3-8/ParaView-3.8/ParaView/VTK/Parallel/vtkPNrrdReader.cxx,
line 365

vtkPNrrdReader (0x1d65a900): Unsupported encoding: gzip

 ERROR: In
/home/kitware/Dashboard/MyTests/ParaView-3-8/ParaView-3.8/ParaView/VTK/Filtering/vtkExecutive.cxx,
line 756

vtkCompositeDataPipeline (0x1d66b110): Algorithm vtkPNrrdReader(0x1d65a900)
returned failure for request: vtkInformation (0x1d675050)

Debug: Off

Modified Time: 71816

Reference Count: 1

Registered Events: (none)

Request: REQUEST_INFORMATION

ALGORITHM_AFTER_FORWARD: 1

FORWARD_DIRECTION: 0



Is there a way to load 'gzipped' image files? I need it to save my local
storage.
Thanks in advance,

Joohwi.
___
Powered by www.kitware.com

Visit other Kitware open-source projects at 
http://www.kitware.com/opensource/opensource.html

Please keep messages on-topic and check the ParaView Wiki at: 
http://paraview.org/Wiki/ParaView

Follow this link to subscribe/unsubscribe:
http://www.paraview.org/mailman/listinfo/paraview