ParaView and Python

Introduction

ParaView offers a rich and powerful Python interface. This allows users to automate processing of their data, and gives access to powerful tools in the Visualization Tool Kit (VTK). This tutorial will describe ParaView and Python. It shows a user how to drive ParaView using Python commands, and how to automate the creation and use of these commands.

Overview

ParaView is a client/ server architecture. The client includes the ParaView GUI and display. The server reads the user's data, processes the data, and passes these images to the client. We can use Python to control ParaView either in the GUI, at the client level, or directly on the server.

A simple Python toy example within ParaView

• Start ParaView.
• Start the Python Interpreter Tools → Python Shell

Notes

You can copy commands from elsewhere and paste them into the Python Shell.

Python is case sensitive. Be sure to use correct capitalization as shown below.

Python is indent sensitive. Be sure to not indent, as shown below.

Lets create and display a sphere.

(Type the following into the Python Shell)

sphere=Sphere()

Show()

Render()

We have now created a sphere in the pipeline, turned on it's visibility, and re-rendered.

Next, lets add a shrink filter. We hide the sphere, add the shrink filter, and re-render.

Hide()

shrink=Shrink()

Show()

Render()

ParaView will allow us to use either the GUI controls or Python. For instance:

Select the Sphere in the pipeline browser.

In the Python Shell, type the following:

clip=Clip()

Show()

Render()

Help! (How do we find out what commands are available?)

To see all commands available in ParaView:

dir()

To see all of the options for the Clip command:

dir(Clip)

To see all of the options for the clip we created above:

dir(clip)

To see lots of detail on a command, create the object and ask for help on that object..

help(clip)

Change! (Lets look at, and change, something)

Print the Theta Resolution

print(sphere.ThetaResolution)

Change it to 64

sphere.ThetaResolution=64

Show()

Render()

Control input

Lets delete the clip

Delete(clip)

Lets add a filter to the sphere, without selecting it first

wireframe=ExtractEdges(Input=sphere)

Show()

Render()

A simple Python example reading a datafile and writing a screenshot

• Within the ParaView GUI, Edit → Reset Session
• Start the Python Interpreter Tools → Python Shell

Read in data, use a filter and save a screenshot

Lets read in can.exo, clip can.exo, paint can.exo and save a screenshot.

canex2=OpenDataFile('D:\\directoryName\\can.ex2')

clip=Clip()

Show()

ResetCamera()

Render()

Hide(clip)

Render()

Show(canex2)

Render()

Hide(canex2)

Show(clip)

Render()

Cut and paste friendly version (replica of above):

canex2=OpenDataFile('D:\\directoryName\\can.ex2')

clip=Clip(); Show(); ResetCamera(); Render()

Hide(clip); Render()

Show(canex2); Render()

Hide(canex2); Show(clip); Render()

SaveScreenshot('D:\\directoryName\\picture.jpg')

Move forward one timestep

We want to move forward one timestep, so min and max are set correctly for a variable

animationScene1 = GetAnimationScene()

animationScene1.GoToNext()

Playing through all time is done with the following command.

animationScene1.Play()

First timestep is found as follows:

animationScene1.GoToFirst()

Playing through all time is done with the following command.

animationScene1.GoToLast()

Available commands are found using:

dir(animationScene1) (after you have created the animationScene1 variable)

dir(GetAnimationScene())

Control the camera

We want to move the camera.

First, get the camera and reset the camera to a known good position.

camera=GetActiveCamera()

camera.SetFocalPoint(0,0,0)

camera.SetPosition(0,0,-10)

camera.SetViewUp(0,1,0)

How to move the camera closer or further away

camera.Dolly(10)

Render()

camera.Dolly(.1)

Render()

How to rotate the camera around the view direction 45 degrees, centered on the dataset. After the reset above, rotate around the X axis.

camera.Roll(45)

Render()

How to rotate the camera around the vector up, centered on the dataset. After the reset above, rotate around the Y axis.

camera.Azimuth(45)

Render()

How to rotate the camera around the vector up, centered on the Y axis. After the reset above, rotate around the Y axis.

camera.Yaw(45)

Render()

Available commands are found using:

dir(camera) (after you have created the animationScene1 variable)

dir(GetActiveCamera())

Paint by a variable

We want to color by the variable.

Be sure to Show the Clip, and not the Can.

canex2.ElementVariables = ['EQPS']

canex2.PointVariables = ['DISPL', 'VEL', 'ACCL']

canex2.GlobalVariables = ['KE', 'XMOM', 'YMOM', 'ZMOM', 'NSTEPS', 'TMSTEP']

renderView1 = GetActiveViewOrCreate('RenderView')

canex2Display = Show(clip, renderView1)

ColorBy(canex2Display, ('POINTS', 'DISPL'))

Render() (Actually not needed)

Cut and paste friendly version (replica of above):

canex2.ElementVariables = ['EQPS']; canex2.PointVariables = ['DISPL', 'VEL', 'ACCL']

canex2.GlobalVariables = ['KE', 'XMOM', 'YMOM', 'ZMOM', 'NSTEPS', 'TMSTEP']

renderView1 = GetActiveViewOrCreate('RenderView'); canex2Display = Show(clip, renderView1)

ColorBy(canex2Display, ('POINTS', 'DISPL'))

Render() (Actually not needed)

Trace Recorder

ParaView includes a tool to automatically generate Python scripts for us. It is called the Trace Recorder. An example is as follows.

Read in can.exo, clip can, paint by EQPS, change the camera to +Y, write out a screenshot and write out a movie

• Tools → Start Trace Select Show Incremental Trace.
• File → Open. Open can.exo. OK.
• Turn all variables on.
• Apply.
• +Y
• Clip. Y Normal. Unselect Show Plane. Apply.
• Color by EQPS.
• Last timestep.
• Rescale to Data Range
• First timestep.
• File → Save Screenshot. Save as .png.
• File → Save Animation. Save as .avi.
• Tools → Stop Trace
• File → Save. Save to a known location.

Another way to find Python for ParaView is through Save State. This should be a last resort, but it may include commands that the Trace Recorder missed. File → Save State → Python State File.

Running Scripts

ParaView allows a user to run a script. This is done as follows:

• Tools → Python Shell
• Run Script

Now, browse to your script, and select OK.

Macros

ParaView can save and use Python scripts that have been placed in a known location. When you create a trace, you have the option to File → Save As Macro. You also have the option on the Macros menu to Add new macro. Macros will be added to the Macro toolbar at the top of the ParaView GUI. You can edit and delete these Macros through the Macro menu.

As an example, lets add the python script that we created above.

• Macros → Add new macro, find your macro, and click OK.
• Click on your Macro on the toolbar.

Example Scripts

First

Second

Third

Forth

Python Help

Python documentation (out of date) http://www.paraview.org/Wiki/ParaView/Python_Scripting

The ParaView Guide (Python scattered throughout the guide) http://www.paraview.org/paraview-guide/

Where do you go next?

• ParaView Batch.

Acknowledgements

Sandia is a multiprogram laboratory operated by Sandia Corporation, a Lockheed Martin Company, for the United States Department of Energy’s National Nuclear Security Administration under contract DE-AC04-94AL85000.