Difference between revisions of "ParaView/Python/Extracting Multiple Blocks"

From KitwarePublic
Jump to navigationJump to search
Line 40: Line 40:
Back to [[ParaView/PythonRecipes]].

Latest revision as of 13:39, 16 October 2018

I had some fun doing this in Python. Here is a script that creates a Python filter for each block and names it with the block name. To run it, you need to select the reader (or any other filter) and then run this script using the Python Shell under the Tools menu. This requires ParaView version >= 3.7 (June 23, 2009).

Another hint: remember that you can turn on/off the visibility of multiple object by selecting them (using shift+click or control+click) in the pipeline browser and clicking on the visibility icon (the eye) of one of them.

from paraview.simple import *

def extract_block(source, outputType, name, di, idx):
   if name:
       pf = ProgrammableFilter(source, registrationName=name)
       pf = ProgrammableFilter(source)
   pf.OutputDataSetType = outputType
   pf.Script = """
input = self.GetInputDataObject(0, 0)
self.GetOutputDataObject(0).ShallowCopy(input.GetBlock(%d))""" % idx
   if outputType == 'vtkImageData' or outputType ==
'vtkStructuredGrid' or outputType == 'vtkRectilinearGrid':
       pf.RequestInformationScript = """
from paraview import util
util.SetOutputWholeExtent(self, %r)""" % list(di.GetExtent())

source = GetActiveSource()
if not source:
   raise RuntimeError, "This macro needs an active pipeline object"
cdi = source.GetDataInformation().GetCompositeDataInformation()
blocks = cdi.GetNumberOfChildren()
for i in range(blocks):
   di = cdi.GetDataInformation(i)
   eb = extract_block(source, di.GetDataClassName(), cdi.GetName(i), di, i)

Back to ParaView/PythonRecipes.

ParaView: [Welcome | Site Map]