Greetings,<br><br>The story: I started using python/numpy to perform data manipulation (3D connected components, etc..) and processing of 3D arrays. Next I managed to use vtk to perform volume rendering etc. to visualize the data. Now, I want to save the vtk data to a format that can be used in paraview so that a side by side comparison with CFD simulations can be performed. I am having trouble with saving the vtk data to a file to be used in paraview.<br>
<br>System: Windows 7, Python 2.7, vtkVersion:04859120.<br><br>1) Numpy array to vtk (example<a href="http://www.vtk.org/Wiki/VTK/Examples/Python/vtkWithNumpy"> here</a>) works fine:<br><br>data=np.array() #sum numpy array dim=20x20x40<br>
dataImporter = vtk.vtkImageImport()<br>data_string = data_matrix.tostring()<br>dataImporter.CopyImportVoidPointer(data_string, len(data_string))<br>dataImporter.SetDataScalarTypeToUnsignedChar()<br>dataImporter.SetNumberOfScalarComponents(1)<br>
dataImporter.SetDataExtent(0, 20, 0, 20, 0, 40)<br>dataImporter.SetWholeExtent(0, 20, 0, 20, 0, 40)<br><br>2) Displaying the volume (example<a href="http://www.vtk.org/Wiki/VTK/Examples/Python/vtkWithNumpy"> here</a>) works fine:<br>
# This class describes how the volume is rendered (through ray tracing).<br>compositeFunction = vtk.vtkVolumeRayCastCompositeFunction()<br># We can finally create our volume. We also have to specify the data for it, as well as how the data will be rendered.<br>
volumeMapper = vtk.vtkVolumeRayCastMapper()<br>volumeMapper.SetVolumeRayCastFunction(compositeFunction)<br>volumeMapper.SetInputConnection(dataImporter.GetOutputPort())<br> <br># The class vtkVolume is used to pair the preaviusly declared volume as well as the properties to be used when rendering that volume.<br>
volume = vtk.vtkVolume()<br>volume.SetMapper(volumeMapper)<br>volume.SetProperty(volumeProperty)<br><br>etc....<br><br>3) Writing data to some vtk file (can't figure out how to do this):<br>writer=vtk.vtkImageWriter()<br>
writer.SetInputConnection(dataImporter.GetOutputPort())<br>writer.SetFileName('./test.vtk')<br>writer.Update()<br>writer.Write()<br><br>When i try opening the test.vtk file with paraview, there is no information in it.<br>
<br>Any help would be much appreciated!<br><br>Thanks<br>Justin<br>