[Paraview] Writing IntegrateAttributes to CSV using Python

Fri Apr 10 04:30:28 EDT 2009

Good morning Utkarsh,

Here's the data replied by the import sys and print sys.path
It shows no old versions. I also included the filters with
dir(servermanager.filters)

Python 2.5.1 (r251:54863, Dec  6 2008, 10:49:39) 
[GCC 4.2.1 (SUSE Linux)] on linux2
>>> from paraview.simple import *
>>> import sys
>>> print sys.path
['/home/bart/bin/lib/paraview-3.5', '/usr/lib/python25.zip',
'/usr/lib64/python2.5', '/usr/lib64/python2.5/plat-linux2',
'/usr/lib64/python2.5/lib-tk', '/usr/lib64/python2.5/lib-dynload',
'/usr/lib64/python2.5/site-packages']
>>> dir(servermanager.filters)
['AlltoN', 'AnnotateTime', 'AppendAttributes', 'AppendDatasets',
'AppendGeometry', 'Balance', 'BlockScalars', 'CacheKeeper', 'Calculator',
'CellCenters', 'CellDatatoPointData', 'Clean', 'CleantoGrid',
'ClientServerMoveData', 'Clip', 'Clip(GenericDatasets)',
'ComputeDerivatives', 'Connectivity', 'Contour',
'Contour(GenericDatasets)', 'ConvertSelection', 'Crop', 'Curvature', 'D3',
'Decimate', 'Delaunay2D', 'Delaunay3D', 'Elevation', 'ExtractAttributes',
'ExtractBlock', 'ExtractCTHFragments', 'ExtractCTHParts',
'ExtractCellsbyRegion', 'ExtractDatasets', 'ExtractEdges', 'ExtractLevel',
'ExtractSelection', 'ExtractSelection(internal)', 'ExtractSubset',
'ExtractSurface', 'ExtractSurface(GenericDatasets)',
'FFTofSelectionOverTime', 'FeatureEdges', 'GenerateIds',
'GenerateQuadraturePoints', 'GenerateQuadratureSchemeDictionary',
'GenerateSurfaceNormals', 'GeometryFilter', 'Glyph', 'Glyph(CustomSource)',
'Gradient', 'Gradient(Unstructured)', 'GradientMagnitude', 'GroupDatasets',
'Histogram', 'ImageShrink', 'IntegrateVariables',
'InterpolatetoQuadraturePoints', 'IntersectCTHFragments', 'LevelScalars',
'LinearExtrusion', 'LoopSubdivision', 'MPIMoveData', 'MaskPoints',
'Median', 'MergeBlocks', 'MeshQuality', 'MinMax', 'NormalGlyphs',
'OctreeDepthLimit', 'OctreeDepthScalars', 'OrderedCompositeDistributor',
'Outline', 'Outline(GenericDataSets)', 'Outline(curvilinear)',
'OutlineCorners', 'ParticlePathlines', 'ParticleTracer',
'PlotGlobalVariablesoverTime', 'PlotOverLine', 'PlotSelectionOverTime',
'PointDatatoCellData', 'PolylinetoRectilinearGrid', 'ProbeLocation',
'ProcessIdScalars', 'ProgrammableFilter', 'PythonCalculator',
'QuadraturePointStatistics', 'QuadricClustering', 'RandomVectors',
'RectilinearGridGeometryFilter', 'ReductionFilter', 'Reflect',
'ResampleWithDataset', 'Ribbon', 'RotationalExtrusion', 'ScatterPlot',
'Shrink', 'Slice', 'Slice(GenericDatasets)', 'Smooth', 'StreamTracer',
'StreamTracer(CustomSource)', 'StreamTracer(GenericDatasets)', 'Subdivide',
'SurfaceFlow', 'SurfaceVectors', 'TableFFT', 'TableToPoints',
'TableToStructuredGrid', 'TemporalCache', 'TemporalInterpolator',
'TemporalShiftScale', 'TemporalSnaptoTimeStep', 'TemporalStatistics',
'Tessellate', 'Tessellator(GenericDatasets)', 'Tetrahedralize',
'TextureMaptoCylinder', 'TextureMaptoPlane', 'TextureMaptoSphere',
'Threshold', 'Transform', 'TriangleStrips', 'Triangulate', 'Tube',
'UpdateSuppressor2', 'Warp(scalar)', 'Warp(vector)', '__class__',
'__delattr__', '__dict__', '__doc__', '__getattribute__', '__hash__',
'__init__', '__module__', '__new__', '__reduce__', '__reduce_ex__',
'__repr__', '__setattr__', '__str__', 

Still, a lot of thanks for helping, would love to see it working in a while
:)

Bart

On Thu, 9 Apr 2009 11:04:49 -0400, Utkarsh Ayachit
<utkarsh.ayachit at kitware.com> wrote:
> Bart,
> 
> Make sure your paraview is built as correctly. Then
> import sys
> print sys.path
> 
> Ensure that the path does not include any old version of paraview
> libraries. If it does then you are ending up import old versions.
> 
> Utkarsh
> 
> On Thu, Apr 9, 2009 at 10:27 AM,  <bart.boonacker at actiflow.com> wrote:
>> hey Utkarsh,
>>
>> I'm sorry that my writing was a bit unclear, but I said that I can find
>> the
>> .cxx and .h files, but I can't find them in the dir(servermanager) and
>> dir(servermanager.filters).
>>
>> The dataStamp.cmake contains this, (today indeed):
>>
>> # Do not edit!  Generated by kwsysDateStamp.py
>>
>> # KWSys version date year component.  Format is CCYY.
>> SET(KWSYS_DATE_STAMP_YEAR  2009)
>>
>> # KWSys version date month component.  Format is MM.
>> SET(KWSYS_DATE_STAMP_MONTH 04)
>>
>> # KWSys version date day component.  Format is DD.
>> SET(KWSYS_DATE_STAMP_DAY   09)
>>
>>
>> Bart
>>
>>
>> On Thu, 9 Apr 2009 10:16:19 -0400, Utkarsh Ayachit
>> <utkarsh.ayachit at kitware.com> wrote:
>>> If there's no Servers/Filters/vtkAttributeDataToTableFilter.h and
>>> Servers/Filters/vtkAttributeDataToTableFilter.cxx then I am guessing
>>> you are not updating correctly. There may be some sticky tags set with
>>> your CVS checkout. Can you look at
>>> ParaView3/VTK/Utilities/kwsys/kwsysDataStamp.cmake and verify that the
>>> date mentioned is indeed today's?
>>>
>>> Utkarsh
>>>
>>> On Thu, Apr 9, 2009 at 9:39 AM,  <bart.boonacker at actiflow.com> wrote:
>>>> Hello again,
>>>>
>>>> I did ran it from the python shell, so that problem is fixed.
>>>> I just updated the CVS and built the update, but still not
>>>> AttributeDataToTableFilter :(
>>>> not in dir(servermanager.filters) nor in dir(servermanager)
>>>>
>>>> When looking for files I can find the
vtkAttributeDataToTableFilter.cxx
>>>> and
>>>> .h but this is not available as well.
>>>>
>>>> Bart
>>>>
>>>>
>>>> On Thu, 9 Apr 2009 09:21:14 -0400, Utkarsh Ayachit
>>>> <utkarsh.ayachit at kitware.com> wrote:
>>>>> Oh, are you running the script through the python shell in the
>>>>> ParaView client? In that case, sorry, but you don't need the
Connect()
>>>>> :). I assumed you were using pvpython.
>>>>>
>>>>> Are you sure you built your paraview after updating? I just tried the
>>>>> CVS and I indeed have the AttributeDataToTableFilter.
>>>>>
>>>>> Utkarsh
>>>>>
>>>>> On Thu, Apr 9, 2009 at 9:18 AM,  <bart.boonacker at actiflow.com> wrote:
>>>>>> Ok, than that should be the problem I guess;
>>>>>>
>>>>>> Here's my "full" script, almost the same as provided before:
>>>>>>
>>>>>> from paraview import servermanager as sm
>>>>>> sm.Connect()
>>>>>>
>>>>>>
>>>>
>>
###############################################################################################
>>>>>> # create implicit device to populate Slice
>>>>>> impPlane =
>>>>>> sm.implicit_functions.Plane(registrationGroup='implicit_fuctions')
>>>>>>
>>>>
>>
###############################################################################################
>>>>>>
>>>>>> caseName = 'Case'
>>>>>> caseLocation = '/home/bart/Desktop/Test_Python/'
>>>>>> showScalarBar = 2 #0=none, 1=vertical, 2=horizontal
>>>>>>
>>>>>> #read data
>>>>>> VTKreader01 =
>>>>>>
>>>>
>>
sm.sources.LegacyVTKreader(FileNames='/home/bart/Desktop/Test_Python/1/VTK/1_1000.vtk',registrationGroup='sources',registrationName='RR')
>>>>>>
>>>>>> Slice_in =
>>>>>>
>>>>
>>
sm.filters.Slice(registrationGroup='sources',registrationName='slice02',Input=VTKreader01)
>>>>>> Slice_inprop = Slice_in.GetProperty('CutFunction')
>>>>>> domainSlice_in = Slice_inprop.GetDomain('proxy_list')
>>>>>> domainSlice_in.AddProxy(impPlane.SMProxy)
>>>>>>
>>>>>> cutplane_in = impPlane
>>>>>> cutplane_in.Normal = [0, 1, 0]
>>>>>> cutplane_in.Offset = 0
>>>>>> cutplane_in.Origin = [0, 0, 0]
>>>>>> Slice_in.SliceType= cutplane_in
>>>>>>
>>>>>> Calculator01 =
>>>>>>
>>>>
>>
sm.filters.Calculator(registrationGroup='sources',registrationName='calculator',Input=Slice_in)
>>>>>> Calculator01.ResultArrayName = 'Pt'
>>>>>> Calculator01.Function = '1.225*p+0.5*1.225*mag(U)^2'
>>>>>> Calculator01.AttributeMode = 1 # 1=pointdata,2=celldata
>>>>>> Calculator01.VectorVariable = ['U','U','0','1','2']
>>>>>> Calculator01.ScalarVariable =
>>>>>> ['p','p','0','U_0','U','0','U_1','U','1','U_2','U','2']
>>>>>>
>>>>>> filter1 =
>>>>>>
>>>>
>>
sm.filters.IntegrateVariables(registrationGroup='sources',registrationName='Integrator',Input=Calculator01)
>>>>>>
>>>>>> ## AttributeDataToTableFilter converts any dataset to a vtkTable
>>>>>> which
>>>>>> can
>>>>>> be
>>>>>> ## written as a CSV file. One needs to choose which attribute to
save
>>>>>> as
>>>>>> well.
>>>>>> ## First save points.
>>>>>> convertor =
>>>>>>
>>>>
>>
sm.filters.AttributeDataToTableFilter(Input=filter1,FieldAssociation="Points",AddMetaData=1)
>>>>>> csvWriter = sm.writers.CSVWriter(Input=convertor)
>>>>>> csvWriter.FileName =
>> "/home/bart/Desktop/Test_Python/output.points.csv"
>>>>>>
>>>>>> ## this writes the output.
>>>>>> csvWriter.UpdatePipeline()
>>>>>>
>>>>>> ### Now save cells.
>>>>>> convertor.FieldAssociation="Cells"
>>>>>> csvWriter = sm.writers.CSVWriter(Input=convertor)
>>>>>> csvWriter.FileName =
>>>>>> "/home/bart/Desktop/Test_Python/output.cells.csv"
>>>>>>
>>>>>> ### this writes the output.
>>>>>> csvWriter.UpdatePipeline()
>>>>>>
>>>>>>
>>>>>> when trying to Connect it gives:
>>>>>> Traceback (most recent call last):
>>>>>>  File "<string>", line 2, in <module>
>>>>>>  File
>>>>>>
>>>>
>>
"/home/bart/ParaView3-bin/Utilities/VTKPythonWrapping/paraview/servermanager.py",
>>>>>> line 1861, in Connect
>>>>>>    raise RuntimeError, "Cannot create a connection through python.
>>>>>> Use
>>>>>> the
>>>>>> GUI to setup the connection."
>>>>>> RuntimeError: Cannot create a connection through python. Use the GUI
>> to
>>>>>> setup the connection.
>>>>>>
>>>>>> Bart
>>>>>>
>>>>>>
>>>>>>
>>>>>> On Thu, 9 Apr 2009 09:13:58 -0400, Utkarsh Ayachit
>>>>>> <utkarsh.ayachit at kitware.com> wrote:
>>>>>>> Can you give your full script? If you are not calling Connect() at
>>>>>>> all, it won't work since Connect() is required.
>>>>>>>
>>>>>>> Utkarsh
>>>>>>>
>>>>>>> On Thu, Apr 9, 2009 at 9:08 AM,  <bart.boonacker at actiflow.com>
>>>>>>> wrote:
>>>>>>>> Utkarsh,
>>>>>>>>
>>>>>>>> I used the CVS from yesterday but it returns :
>>>>>>>>      Traceback (most recent call last):
>>>>>>>>      File "<string>", line 41, in <module>
>>>>>>>>      AttributeError: 'PVModule' object has no attribute
>>>>>>>> 'AttributeDataToTableFilter'
>>>>>>>>
>>>>>>>> I have not used the following as the Connect part did not work,
>>>> perhaps
>>>>>>>> thats the cause?
>>>>>>>> I replaced the sm by servermanager.
>>>>>>>>       from paraview import servermanager as sm
>>>>>>>>       sm.Connect()
>>>>>>>>
>>>>>>>> Bart
>>>>>>>>
>>>>>>>>
>>>>>>>> On Thu, 9 Apr 2009 08:55:25 -0400, Utkarsh Ayachit
>>>>>>>> <utkarsh.ayachit at kitware.com> wrote:
>>>>>>>>> Bart,
>>>>>>>>>
>>>>>>>>> Please cc the mailing list as well so other can use this
>> information
>>>>>>>>> as well as respond.
>>>>>>>>>
>>>>>>>>> I think one dubious statement I see in your code is:
>>>>>>>>> ## First save points.
>>>>>>>>> convertor =
>>>> sm.filters.AttributeDataToTableFilter(Input=filter1Output,
>>>>>>>>>   FieldAssociation="Points",
>>>>>>>>>   AddMetaData=1)
>>>>>>>>>
>>>>>>>>> You are setting the input to this filter be the data object
>> returned
>>>>>>>>> by using Fetch(). In ParaView python scripts the input to filters
>>>>>>>>> cannot be other data objects, but filters or sources. So, the
>> proper
>>>>>>>>> usage would be:
>>>>>>>>>
>>>>>>>>> convertor = sm.filters.AttributeDataToTableFilter(Input=filter1,
>>>>>>>>>   FieldAssociation="Points",
>>>>>>>>>   AddMetaData=1)
>>>>>>>>>
>>>>>>>>> Does that solve the issue? Also creation of the representation is
>>>>>>>>> not
>>>>>>>>> necessary unless you are rendering into a view, which doesn't
seem
>>>>>>>>> to
>>>>>>>>> be the case.
>>>>>>>>>
>>>>>>>>> Utkarsh
>>>>>>>>>
>>>>>>>>> On Thu, Apr 9, 2009 at 8:32 AM,  <bart.boonacker at actiflow.com>
>>>>>>>>> wrote:
>>>>>>>>>> Dear Utkarsh,
>>>>>>>>>>
>>>>>>>>>> Thank you for your fast replies. I got the CVS version to work,
>> but
>>>> I
>>>>>>>>>> am
>>>>>>>>>> still having some trouble to rewrite all the filters and
>> functions.
>>>>>>>>>> When running my code, see below, Paraview crashes and gives a
>>>>>>>>>> Segmentation
>>>>>>>>>> Fault.
>>>>>>>>>> I am not a real programming expert, so I tried to neglect the
>>>>>>>>>> representations of all filters. This removes the crash, but
still
>>>>>>>>>> posses
>>>>>>>>>> some problem with the code. The code was 1-to-1 copied from the
>> 3.4
>>>>>>>>>> version, which probably causes the errors.
>>>>>>>>>> Is there a way or manual to see the new options and maybe
>>>> conversions
>>>>>>>>>> between 3.4 and the CVS version ?
>>>>>>>>>>
>>>>>>>>>> Thanks in advance,
>>>>>>>>>>
>>>>>>>>>> Bart
>>>>>>>>>>
>>>>>>>>>> Attached code ( I removed all representation code and view
>>>> settings):
>>>>>>>>>>
>>>>>>>>>>
>>>>>>>>>>
>>>>>>>>
>>>>>>
>>>>
>>
###############################################################################################
>>>>>>>>>> # create implicit device to populate Slice
>>>>>>>>>> impPlane =
>>>>>>>>>>
>>>>>>>>
>>>>>>
>>>>
>>
servermanager.implicit_functions.Plane(registrationGroup='implicit_fuctions')
>>>>>>>>>> impBox =
>>>>>>>>>>
>>>>>>>>
>>>>>>
>>>>
>>
servermanager.implicit_functions.Box(registrationGroup='implicit_fuctions')
>>>>>>>>>> impSphere =
>>>>>>>>>>
>>>>>>>>
>>>>>>
>>>>
>>
servermanager.implicit_functions.Sphere(registrationGroup='implicit_fuctions')
>>>>>>>>>>
>>>>>>>>
>>>>>>
>>>>
>>
###############################################################################################
>>>>>>>>>>
>>>>>>>>>> caseName = 'Case'
>>>>>>>>>> caseLocation = '/home/bart/Desktop/Test_Python/'
>>>>>>>>>> showScalarBar = 2 #0=none, 1=vertical, 2=horizontal
>>>>>>>>>>
>>>>>>>>>> #read data
>>>>>>>>>> VTKreader01 =
>>>>>>>>>>
>>>>>>>>
>>>>>>
>>>>
>>
servermanager.sources.LegacyVTKreader(FileNames='/home/bart/Desktop/Test_Python/1/VTK/1_1000.vtk',registrationGroup='sources',registrationName='RR')
>>>>>>>>>> #repVTKreader01 =
>>>>>>>>>>
>>>>>>>>
>>>>>>
>>>>
>>
servermanager.CreateRepresentation(VTKreader01,view,registrationGroup='representations')
>>>>>>>>>> #repVTKreader01.Visibility = 0
>>>>>>>>>> #repVTKreader01.Representation = 3
>>>>>>>>>>
>>>>>>>>>> Slice02 =
>>>>>>>>>>
>>>>>>>>
>>>>>>
>>>>
>>
servermanager.filters.Slice(registrationGroup='sources',registrationName='slice02',Input=VTKreader01)
>>>>>>>>>> cfSlice02 = Slice02.GetProperty('CutFunction')
>>>>>>>>>> domainSlice02 = cfSlice02.GetDomain('proxy_list')
>>>>>>>>>> domainSlice02.AddProxy(impPlane.SMProxy)
>>>>>>>>>> domainSlice02.AddProxy(impBox.SMProxy)
>>>>>>>>>> domainSlice02.AddProxy(impSphere.SMProxy)
>>>>>>>>>>
>>>>>>>>>> cutplane02 = impPlane
>>>>>>>>>> cutplane02.Normal = [0, 1, 0]
>>>>>>>>>> cutplane02.Offset = 0
>>>>>>>>>> cutplane02.Origin = [0, 0, 0]
>>>>>>>>>>
>>>>>>>>>> #set cutplane
>>>>>>>>>> Slice02.CutFunction= cutplane02
>>>>>>>>>>
>>>>>>>>>> Calculator01 =
>>>>>>>>>>
>>>>>>>>
>>>>>>
>>>>
>>
servermanager.filters.Calculator(registrationGroup='sources',registrationName='calculator',Input=Slice02)
>>>>>>>>>> Calculator01.ResultArrayName = 'Pt'
>>>>>>>>>> Calculator01.Function = '1.225*p+0.5*1.225*mag(U)^2'
>>>>>>>>>> Calculator01.AttributeMode = 1 # 1=pointdata,2=celldata
>>>>>>>>>> Calculator01.AddVectorVariable = ['U','U','0','1','2']
>>>>>>>>>> Calculator01.AddScalarVariable =
>>>>>>>>>> ['p','p','0','U_0','U','0','U_1','U','1','U_2','U','2']
>>>>>>>>>>
>>>>>>>>>> filter1 =
>>>>>>>>>>
>>>>>>>>
>>>>>>
>>>>
>>
servermanager.filters.IntegrateAttributes(registrationGroup='sources',registrationName='Integrator',Input=Calculator01)
>>>>>>>>>> filter1Output = servermanager.Fetch(filter1)
>>>>>>>>>> #repfilter1 =
>>>>>>>>>>
>>>>>>>>
>>>>>>
>>>>
>>
servermanager.CreateRepresentation(filter1,view,registrationGroup='representations')
>>>>>>>>>>
>>>>>>>>>> # From Utkarsh Ayachit:
>>>>>>>>>> from paraview import servermanager as sm
>>>>>>>>>> sm.Connect()
>>>>>>>>>>
>>>>>>>>>> ## AttributeDataToTableFilter converts any dataset to a vtkTable
>>>>>>>>>> which
>>>>>>>>>> can
>>>>>>>>>> be
>>>>>>>>>> ## written as a CSV file. One needs to choose which attribute to
>>>> save
>>>>>>>>>> as
>>>>>>>>>> well.
>>>>>>>>>> ## First save points.
>>>>>>>>>> convertor =
>>>>>>>>>> sm.filters.AttributeDataToTableFilter(Input=filter1Output,
>>>>>>>>>> FieldAssociation="Points",
>>>>>>>>>> AddMetaData=1)
>>>>>>>>>> csvWriter = sm.writers.CSVWriter(Input=convertor)
>>>>>>>>>> csvWriter.FileName =
>>>>>>>>>> "/home/bart/Afstuderen/Straight_Auto/output.points.csv"
>>>>>>>>>>
>>>>>>>>>> ## this writes the output.
>>>>>>>>>> csvWriter.UpdatePipeline()
>>>>>>>>>>
>>>>>>>>>> ## Now save cells.
>>>>>>>>>> convertor.FieldAssociation="Cells"
>>>>>>>>>> csvWriter = sm.writers.CSVWriter(Input=convertor)
>>>>>>>>>> csvWriter.FileName =
>>>>>>>>>> "/home/bart/Afstuderen/Straight_Auto/output.cells.csv"
>>>>>>>>>>
>>>>>>>>>> ## this writes the output.
>>>>>>>>>> csvWriter.UpdatePipeline()
>>>>>>>>>>
>>>>>>>>>>
>>>>>>>>>>
>>>>>>>>
>>>>>>
>>>>
>>