CinemaIO¶

Pipeline description¶

This example first loads a cinema database of a simulation from disk, consisting of three dimensional image files, using the CinemaReader. This outputs a vtkTable.

The database is queried for a selection of images, using CinemaQuery, which supports SQL queries on vtkTables (bottom view shows query result in a spreadsheet view).

Each selected entry in the database is read by the CinemaProductReader, which outputs a vtkMultiBlock of the images.

The images are sliced with a plane, and each slice is visualized side-by-side using the GridLayout (top view in screenshot).

ForEach is used to loop through all slices, and then the ArrayEditor is used to add a FieldData value to each slice. In this case, we add the interval SampleInterval between each queried entry. Each slice is then written to a new cinema database with the CinemaWriter. Finally, the for-loop is terminated using EndFor.

ParaView¶

To reproduce the above screenshot, go to your ttk-data directory and enter the following command:

$ paraview --state=states/cinemaIO.pvsm

Python code¶

#!/usr/bin/env python

from paraview.simple import *

# create a new 'TTK CinemaReader'
viscousFingerscdb = TTKCinemaReader(DatabasePath='ViscousFingers.cdb')

# create a new 'TTK CinemaQuery'
tTKCinemaQuery1 = TTKCinemaQuery(InputTable=viscousFingerscdb)
tTKCinemaQuery1.SQLStatement = """SELECT * FROM InputTable0
WHERE Sim='run01' AND Time%10=0
ORDER BY Time
LIMIT 8 OFFSET 1"""

# create a new 'TTK CinemaProductReader'
tTKCinemaProductReader1 = TTKCinemaProductReader(Input=tTKCinemaQuery1)

# create a new 'Slice'
slice1 = Slice(Input=tTKCinemaProductReader1)
slice1.SliceType = 'Plane'
slice1.HyperTreeGridSlicer = 'Plane'
slice1.SliceOffsetValues = [0.0]

# init the 'Plane' selected for 'SliceType'
slice1.SliceType.Origin = [31.5, 31.5, 31.5]

# create a new 'TTK GridLayout'
tTKGridLayout1 = TTKGridLayout(Input=slice1)
tTKGridLayout1.ColumnAxis = 'Y'
tTKGridLayout1.ColumnGap = 8.0
tTKGridLayout1.RowAxis = 'Z'
tTKGridLayout1.NumberofRows = 1

# create a new 'TTK ForEach'
tTKForEach1 = TTKForEach(Input=tTKGridLayout1)
tTKForEach1.IterationMode = 'Block'
tTKForEach1.InputArray = ['POINTS', 'ImageFile']
tTKForEach1.OutputType = 'vtkPolyData'

# create a new 'TTK ArrayEditor'
tTKArrayEditor1 = TTKArrayEditor(Target=tTKForEach1,
    Source=None)
tTKArrayEditor1.TargetAttributeType = 'Field Data'
tTKArrayEditor1.DataString = 'SampleInterval, 10'
tTKArrayEditor1.TargetArray = ['POINTS', 'ImageFile']

# create a new 'TTK CinemaWriter'
tTKCinemaWriter1 = TTKCinemaWriter(Input=tTKArrayEditor1,
    DatabasePath='ViscousFingersSampled.cdb')
tTKCinemaWriter1.ScalarField = ['POINTS', 'ImageFile']

# create a new 'TTK EndFor'
tTKEndFor1 = TTKEndFor(Data=tTKCinemaWriter1,
    For=tTKForEach1)

UpdatePipeline()

Inputs¶

ViscousFingers.cdb: a cinema database of VTK Image files from a simulation.

Outputs¶

ViscousFingersSampled.cdb: a cinema database containing the sampled slices of the input cinema database.

C++/Python API¶