In [1]:
import vtk
import pygslib
from vtkdatawidgets.vtk_binding import VtkJupyterBridge
In [2]:
def ReadImageData(path):
"""ReadImageData(object grid, str path)
Read a vtkImageData ('*.vti')
Parameters
----------
path : str
Extension (*.vti) will be added if not provided
Returns
-------
object : vtkImageData
"""
# add extension to path
if not path.lower().endswith('.vti'):
path = path + '.vti'
reader = vtk.vtkXMLImageDataReader()
reader.SetFileName(path)
reader.Update()
return reader.GetOutput()
In [9]:
reader = vtk.vtkXMLImageDataReader()
reader.SetFileName('vamap1.vti')
reader.Update()
vario3d = ReadImageData('vamap1.vti')
In [10]:
bridge = VtkJupyterBridge()
In [7]:
bridge.SetInputDataObject(vario3d)
bridge.Update()
In [8]:
from vtkdatawidgets import VtkRenderer
renderer = VtkRenderer(dataset=bridge.widget, background=(0.5, 0, 0), size=(600, 400))
renderer
In [ ]: