Search code examples
pythonvtkopenfoam

vtk filter for OpenFOAM


I want to filter a pressure by thresholding using python-vtk

    import vtk
    filename = "a.foam"

    reader = vtk.vtkOpenFOAMReader()
    reader.SetFileName(filename)
    reader.CreateCellToPointOn()
    reader.DecomposePolyhedraOn()
    reader.EnableAllCellArrays()
    reader.Update()

    tArray = vtk_to_numpy(reader.GetTimeValues())
    reader.UpdateTimeStep(tArray[-1]) 
    reader.Update()

    filter_threshold = vtk.vtkThreshold()
    filter_threshold.SetInputConnection(reader.GetOutputPort()) 

how should I select "pressure" on filter class ?


Solution

  • Use SetInputArrayToProcess method.

    https://vtk.org/doc/nightly/html/classvtkAlgorithm.html#a6bea16e1329609dbccce0ff8d2367484