0

I have two vtk files obtained after the simulation that needs to be merged for visualization. Initially, I had written the code with vtkPolyData, but since it is unstructured, I changed it to vtkAppendFiler. In both case, I am unable to see any visualization although the files contains the output contents.

Below is the code structure:

from vtk import vtkAppendFilter, vtkUnstructuredGridWriter, vtkPolyData, vtkUnstructuredGridReader

reader = vtkUnstructuredGridReader()
append = vtkAppendFilter()

filenames = ['file1.vtk', 'file2.vtk']
for file in filenames:
    reader.SetFileName(file)
    reader.Update()
    polydata = vtkPolyData()
    polydata.ShallowCopy(reader.GetOutput())
    append.AddInputData(polydata)

append.Update()    

writer = vtkUnstructuredGridWriter()
#writer = vtkPolyDataWriter()
writer.SetFileName('output.vtk')
writer.SetInputData(append.GetOutput())
writer.Write()
rasunag27
  • 55
  • 4
  • 11

1 Answers1

1

The output type of the reader is vtkUnstructuredGrid. As you put it into a vtkpolydata, you try to append empty datasets.

Edit

A vtk file contains only one dataset. Its type is defined at the top of the file and does not depends on the reader. If you know that in your case it is always an Unstructured Grid, you can use the vtkUnstructuredGridReader. And the output will always be a vtkUnstructuredGrid.

You can change your script:

- polydata = vtkPolyData()
- polydata.ShallowCopy(reader.GetOutput())
+ unstructured = vtkUnstructuredGrid()
+ unstructured.ShallowCopy(reader.GetOutput())
Nico Vuaille
  • 2,310
  • 1
  • 6
  • 14
  • You mean the extraction of polydata from vtk files..? As written with `polydata=vtkPolyData()` ..? How do I change it so that it reads only UnstructuredGrid? – rasunag27 Sep 27 '21 at 10:02
  • answer updated. Also for ref, [this cross post](https://discourse.vtk.org/t/how-to-merge-vtk-files/6678) – Nico Vuaille Sep 27 '21 at 11:51