ParaView/Python/Tracking an object during animation

From KitwarePublic
< ParaView
Revision as of 02:59, 19 May 2011 by Raymondc (Talk | contribs)

(diff) ← Older revision | Latest revision (diff) | Newer revision → (diff)
Jump to: navigation, search

Here is a script that I wrote for someone on the mailing list. Add it as a macro from the Python shell for easy access. Essentially, this script sets up tracking for the active object in the active view. Look at the source file for instructions on how to use it. Note that the tracking is pretty basic: the camera moves with the center of the spatial bounds of the object. It would be easy to change the code to do fancier things. If you want to extract a block/cell/subset/whatever, first extract whatever you want tracked using the selection mechanism + extract selection filter and add tracking for that.

# This script demonstrates how to track a moving object.
# Each time it is executed, it marks the active view
# and object. It also adds a track to the animation
# that moves the camera with the center of the bounds
# of that object.
# To use, setup the view like you want the camera to
# look, run this script and then run the animation.
# To disable tracking, select the server object in the
# pipeline browser and rerun this script.
from paraview.simple import *
ans = None
    ans = GetAnimationScene()
except NameError:
    for i in servermanager.ProxyManager().GetProxiesInGroup("animation").values():
        if i.IsA("vtkSMAnimationSceneProxy"):
            ans = i
if ans:
    # Turn of caching so that we can get the bounds for the
    # data for each time step
    ans.Caching = False
    scene = ans.GetClientSideObject()
    # Remove the previous observer.
    except NameError:
    def callback(caller, *args):
        global rv, source, first_bounds, first_focal, first_pos, view_up
        if source:
            #The following line does not work in Paraview 3.10.1 -- the error printed is "AttributeError: GetViewUpdateTime". The script seems to work fine without it.
            bds = source.GetDataInformation().GetBounds()
            # Center of bounds
            bds = map(lambda x,y : (x+y)/2, bds[0:6:2], bds[1:6:2])
            # Difference between initial bounds and current bounds
            bds_diff = map(lambda x,y : x-y, bds, first_bounds)
            # Move the focal point as much as the center of the bounds
            # moved
            rv.CameraFocalPoint = map(lambda x,y: x+y, first_focal, bds_diff)
            # Move the camera as much as the center of the bounds
            # moved
            rv.CameraPosition = map(lambda x,y: x+y, first_pos, bds_diff)
            # Keep the original view up
            rv.CameraViewUp = view_up
    rv = GetActiveView()
    if rv:
        # Add the observer to the VTK scene with priority 1 so that it happens before
        # anything else
        callbackid = scene.AddObserver("AnimationCueTickEvent", callback, 1.0)
        source = GetActiveSource()
        bds = source.GetDataInformation().GetBounds()
        # Center of bounds
        first_bounds = map(lambda x,y : (x+y)/2, bds[0:6:2], bds[1:6:2])
        first_focal = rv.CameraFocalPoint[:]
        first_pos = rv.CameraPosition[:]
        view_up = rv.CameraViewUp[:]
def removeCameraObserver():
    global callbackid, scene