apply a CleanToGrid filter before the python programmable filter and
this is a corrected version of your script:
pdi = self.GetInput()
pdo = self.GetOutput()
newPoints = vtk.vtkPoints()
numPoints = pdi.GetNumberOfPoints()
for i in range(0, numPoints):
coord = pdi.GetPoint(i)
x, y, z = coord[:3]
r = x * sin(y)
t = x * cos(y)
newPoints.InsertPoint(i, r, t, z)
pdo.SetPoints(newPoints)
Felipe
Le 03/04/2013 04:55, Nikolaos Beratlis a écrit :
Hi,
I am reading the following xmf file (attached as var3d.xmf) into
Paraview:
<?xml version="1.0" ?>
<!DOCTYPE Xdmf SYSTEM "Xdmf.dtd" []>
<Xdmf Version="2.2">
<Domain>
<Grid GridType="Uniform">
<Topology TopologyType="3DRectMesh" Dimensions=" 2 4 3"/>
<Geometry GeometryType="VXVYVZ">
<DataItem Dimensions=" 3" NumberType="Float" Precision="4"
Format="XML">
0.0000 1.0000 2.0000
</DataItem>
<DataItem Dimensions=" 4" NumberType="Float" Precision="4"
Format="XML">
0.0000 2.0944 4.1888 6.2832
</DataItem>
<DataItem Dimensions=" 2" NumberType="Float" Precision="4"
Format="XML">
1.0000 2.0000
</DataItem>
</Geometry>
<Attribute Name="var" AttributeType="Scalar" Center="Node">
<DataItem Dimensions="24" NumberType="Float" Precision="4" Format="XML">
0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23
</DataItem>
</Attribute>
</Grid>
</Domain>
</Xdmf>
The grid is a 3D orthogonal grid read in cartesian coordinates (see
attached image showing grid and contours of variable). I then want to
transform the grid from cartesian to cylindrical coordinates. I tried
using the following Python script:
pdi = self.GetPolyDataInput()
pdo = self.GetPolyDataOutput()
newPoints = vtk.vtkPoints()
numPoints = pdi.GetNumberOfPoints()
for i in range(0, numPoints):
coord = pdi.GetPoint(i)
x, y, z = coord[:3]
r = x * sin(y)
t = x * cos(y)
newPoints.InsertPoint(i, r, t, z)
pdo.SetPoints(newPoints)
with Output Data Set Type chosen to be the Same as Input, but I get
this error:
Traceback (most recent call last):
File "<string>", line 26, in <module>
File "<string>", line 12, in RequestData
AttributeError: 'NoneType' object has no attribute 'SetPoints'
I have no experience with Python scripting. What am I doing wrong?
Thank you,
Nikos
_______________________________________________
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
--
Felipe Bordeu Weldt
Ingénieur de Recherche
-------------------------------------
Tél. : 33 (0)2 40 37 16 57
Fax. : 33 (0)2 40 74 74 06
felipe.bor...@ec-nantes.fr
Institut GeM - UMR CNRS 6183
École Centrale Nantes
1 Rue de La Noë, 44321 Nantes, FRANCE
-------------------------------------
_______________________________________________
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