[ginop]

I would like to use foamToVTK generated files for simulation post-process in python.

At the moment, I have difficulties in getting cell centers. I read the vtk file, and get the data arrays for each cell. BUT how do I get the cell centers/points list?

I am a bit lost getting cell IDs, points, connectivities really...

So I do:

reader = vtk.vtkUnstructuredGridReader()
reader.SetFileName(filepath)
reader.ReadAllFieldsOn()
reader.ReadAllVectorsOn()
reader.ReadAllScalarsOn()
reader.Update()

vtkdata = reader.GetOutput()
celldata = vtkdata.GetCellData()

And that's okay for the actual data (U, p, etc.).

Then I tried

cells = vtkdata.GetCells()

but still Icould not manage to get point IDs/coords for each cell?

From the VTK user manual, it's not completely clear to me (haven't gone through the whole thing though...).
I hope someone can give me some light on this.

Cheers,
Gino

[ginop]

Found a solution thanks to this post

The code extension should look like:

cells = vtkdata.GetCells()
cellpoints = vtk_to_numpy(cells.GetData())
i=0
cell_point_id = list()
while i < cells.GetSize():
    pts = cellpoints[i]
    cell_point_id.append(cellpoints[i+1:i+pts+1])
    i += pts + 1