ParaView/Automation with Python

1. Extract surface using paraview, the code will look like

from paraview.simple import *

#========== read in data file
print 'Read phasta file'
dataObj = PhastaReader( FileName= 'restart_one.pht' )  # you have to change the name of this pht file
SRF = ExtractSurface(Input=dataObj)
MGBLK1 = MergeBlocks(Input=SRF)
vtkWriter = DataSetWriter(FileName = 'Surface.vtk' , Input = MGBLK1, FileType = 1) # you have to change the name of this vtk file
vtkWriter.UpdatePipeline()

Get the python file from

/users/ychen/NGC/geometry2012/pvExtractSuface.py

Generating Python Scripts for Data Extraction

Generating python scripts can be done in the Paraview GUI and is actually pretty easy and straightforward. This process is documented here in an example. Once the the script you want is generated, that script can be edited to automate data extraction. Again, a general example of this is outlined here. Finally, once the script is run, data will be exported and look something like this.