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