Generic Workflow

This shows an example workflow with vapor. We begin by creating a session and opening a dataset. You can have multiple sessions open at the same time.

import example_utils
from vapor import session, renderer, dataset, camera

Warning: sysroot "/Applications/Xcode_12.4.app/Contents/Developer/Platforms/MacOSX.platform/Developer/SDKs/MacOSX.sdk" not found (ignoring for now).

Vapor 3.9.3
Python 3.9.19 (/opt/anaconda3/envs/vapor)
OpenGL 4.1 Metal - 83.1

ses = session.Session()
data = example_utils.OpenExampleDataset(ses)

Dump the dataset metadata

print("Time Coordinate Variable Name:", data.GetTimeCoordVarName())
print("Coordinate Variable Names:", data.GetCoordVarNames())

print("Dimensions:")
for dim in data.GetDimensionNames():
    print(f"  {dim}:", data.GetDimensionLength(dim, 0))

print("Data Variables:")
for var in data.GetDataVarNames():
    print(f"  {var}")
    print(f"    Time Varying:", bool(data.IsTimeVarying(var)))
    print(f"    Dimensionality:", data.GetVarGeometryDim(var))
    print(f"    Coordinates:", data.GetVarCoordVars(var, True))
    print("     Data Range:", data.GetDataRange(var))

Time Coordinate Variable Name: 
Coordinate Variable Names: ['__regCoord_64_x', '__regCoord_64_y', '__regCoord_64_z']
Dimensions:
  __regDim_64: 64
Data Variables:
  U10
    Time Varying: False
    Dimensionality: 2
    Coordinates: ['__regCoord_64_x', '__regCoord_64_y']
     Data Range: [-1.9993805885314941, 1.9994730949401855]
  V10
    Time Varying: False
    Dimensionality: 2
    Coordinates: ['__regCoord_64_x', '__regCoord_64_y']
     Data Range: [-0.2711198925971985, 0.2711198925971985]
  V
    Time Varying: False
    Dimensionality: 3
    Coordinates: ['__regCoord_64_x', '__regCoord_64_y', '__regCoord_64_z']
     Data Range: [-0.3319709300994873, 0.3319709300994873]

Render the first 2D variable as a pseudocolor

first_2d_var = data.GetDataVarNames(2)[0]
print(f"Rendering 2D variable {first_2d_var}")

ren = data.NewRenderer(renderer.TwoDDataRenderer)
ren.SetVariableName(first_2d_var)
ren.GetPrimaryTransferFunction().SetMinMapValue(-1)
ren.GetPrimaryTransferFunction().SetMaxMapValue(1)

ses.GetCamera().ViewAll()
ses.Show()

ses.DeleteRenderer(ren)

Rendering 2D variable U10
../_images/bd262a4f26896cfc0a4b0b74cd17da8db9684e4b5a0a8aa279aa5464a3b1fa84.png

Render U10 and V10 as barbs over a map

# If your dataset is geo-referenced, this will automatically render a geographically correct map.
# map_ren = data.NewRenderer(renderer.ImageRenderer)

barbs = data.NewRenderer(renderer.BarbRenderer)
barbs.SetDimensions(2)
barbs.SetFieldVariableNames(['U10', 'V10'])
barbs.SetLineThickness(2)

ses.Show()
ses.DeleteRenderer(barbs)
../_images/d473047c12f5e4996855e951f18a8e3618ff3cbbb36e0939f45411ab25a92b43.png

Simulate and render a flow advection

flow = data.NewRenderer(renderer.FlowRenderer)
flow.SetFieldVariableNames(['U10', 'V10'])

ses.Show()

ses.DeleteRenderer(flow)
../_images/cb04e6e8c0107ebf850daa4f60bdee8db80d3935c8c2679de388641c9285bb38.png

Volume render a 3D variable

volume = data.NewRenderer(renderer.VolumeRenderer)
volume.SetVariableName("V")

tf = volume.GetPrimaryTransferFunction()
tf.SetOpacityList([1, 0, 0, 1])

ses.GetCamera().ViewAll()
ses.Show()

# Show a colorbar for the volume rendering
tf.ShowMatPlotLibColorbar(label="V")

UNSUPPORTED (log once): POSSIBLE ISSUE: unit 3 GLD_TEXTURE_INDEX_3D is unloadable and bound to sampler type (Float) - using zero texture because texture unloadable
../_images/8a4b3fbaf0b08ec9f21b4a178a1a1c06302ee503fc8da0b576f3190e1e263ece.png ../_images/377c6c3b4d0ea73be775ffe4a365913836a60cfa8978cb9b40ea4d93ff3e719b.png

Scale the dataset Z axis

Vapor will automatically scale the Z axis of a dataset to produce reasonable results. This can be manually adjusted as shown below.

scales = data.GetTransform().GetScales()
print("Default dataset scaling =", scales);
scales[2] *= 0.3
data.GetTransform().SetScales(scales)
print("New dataset scaling =", data.GetTransform().GetScales());
ses.Show()

Default dataset scaling = [1.0, 1.0, 1.0]
New dataset scaling = [1.0, 1.0, 0.3]
../_images/16380d51fa5dd906e5d78dcb55d8aa78ba5dbd7d4f6a64238344849f103669a7.png

Add axis annotations to the rendering

annotations = ses.GetAxisAnnotations()
annotations.SetAxisAnnotationEnabled(True)
ses.Show()
../_images/38438cb34d58cdf7f3f9a26c05dfe812c986c7ec5583e3dba138597099b22646.png

Export the session for use in the Vapor GUI application

Sessions created in Python can be saved as a .vs3 file. These files can then be opened in the Vapor GUI application and explored interactively. Conversely, sessions created in the Vapor GUI can be loaded into Python with Session.Load(path)

Since this example uses a dynamically generated dataset, the session cannot be saved as it would point to a dataset that does not exist on disk. If you were using a physical dataset, this would work.

ses.Save("tutorial.vs3")

[Session.cpp:121] Cannot save session that contains data dynamically loaded from python (dataset.PYTHON)

-1
On this page