def __init__(self):
self.fig = mlab.figure()
self.scene = self.fig.scene
src = ArraySource(transpose_input_array=False)
src.scalar_data = np.random.randn(10,10,10)
src.origin = [-50,-50,-50]
src.spacing = [10,10,10]
mlab.pipeline.image_plane_widget(src, figure=self.fig)
def setUp(self):
"""Initial setting up of test fixture, automatically called by TestCase before any other test method is invoked"""
e = NullEngine()
# Uncomment to see visualization for debugging etc.
#e = Engine()
e.start()
s = e.new_scene()
self.e = e
self.s = s
############################################################
# Create a new scene and set up the visualization.
d = ArraySource()
sc, vec = self.make_data()
d.origin = (-5, -5, -5)
d.scalar_data = sc
d.vector_data = vec
e.add_source(d)
# Create an outline for the data.
o = Outline()
e.add_module(o)
# View the data.
st = Streamline()
e.add_module(st)
widget = st.seed.widget
widget.set(radius=1.0,
center=(-4.0, -4.0, -4.0),
theta_resolution=4,
phi_resolution=4)
st = Streamline(streamline_type='ribbon')
seed = st.seed
seed.widget = seed.widget_list[1]
e.add_module(st)
seed.widget.set(point1=(-5.0, -4.5, -4.0), point2=(-5.0, -4.5, 4.0))
st.ribbon_filter.width = 0.25
st = Streamline(streamline_type='tube')
seed = st.seed
seed.widget = seed.widget_list[2]
e.add_module(st)
seed.widget.set(center=(-5.0, 1.5, -2.5))
st.tube_filter.radius = 0.15
st = Streamline(streamline_type='tube')
seed = st.seed
seed.widget = seed.widget_list[3]
e.add_module(st)
seed.widget.position = (-5.0, 3.75, 3.75)
st.tube_filter.radius = 0.2
self.st = st
self.scene = e.current_scene
return
def setUp(self):
"""Initial setting up of test fixture, automatically called by TestCase before any other test method is invoked"""
e = NullEngine()
# Uncomment to see visualization for debugging etc.
#e = Engine()
e.start()
s=e.new_scene()
self.e=e
self.s=s
############################################################
# Create a new scene and set up the visualization.
d = ArraySource()
sc, vec = self.make_data()
d.origin = (-5, -5, -5)
d.scalar_data = sc
d.vector_data = vec
e.add_source(d)
# Create an outline for the data.
o = Outline()
e.add_module(o)
# View the data.
st = Streamline()
e.add_module(st)
widget = st.seed.widget
widget.set(radius=1.0, center=(-4.0, -4.0, -4.0),
theta_resolution=4, phi_resolution=4)
st = Streamline(streamline_type='ribbon')
seed = st.seed
seed.widget = seed.widget_list[1]
e.add_module(st)
seed.widget.set(point1=(-5.0, -4.5, -4.0), point2=(-5.0, -4.5, 4.0))
st.ribbon_filter.width = 0.25
st = Streamline(streamline_type='tube')
seed = st.seed
seed.widget = seed.widget_list[2]
e.add_module(st)
seed.widget.set(center=(-5.0, 1.5, -2.5))
st.tube_filter.radius = 0.15
st = Streamline(streamline_type='tube')
seed = st.seed
seed.widget = seed.widget_list[3]
e.add_module(st)
seed.widget.position=(-5.0, 3.75, 3.75)
st.tube_filter.radius = 0.2
self.st = st
self.scene = e.current_scene
return
def setUp(self):
"""Initial setting up of test fixture, automatically called by TestCase before any other test method is invoked"""
e = NullEngine()
# Uncomment to see visualization for debugging etc.
#e = Engine()
e.start()
s=e.new_scene()
self.e=e
self.s=s
############################################################
# Create a new scene and set up the visualization.
d = ArraySource()
sc, vec = self.make_data()
d.origin = (-5, -5, -5)
d.scalar_data = sc
d.vector_data = vec
e.add_source(d)
# Create an outline for the data.
o = Outline()
e.add_module(o)
# Glyphs for the scalars
g = Glyph()
e.add_module(g)
g.glyph.glyph_source.glyph_position = 'center'
g.glyph.glyph.vector_mode = 'use_normal'
g.glyph.glyph.scale_factor = 0.5
g.actor.property.line_width = 1.0
v = VectorCutPlane()
glyph = v.glyph
gs = glyph.glyph_source
gs.glyph_position = 'tail'
gs.glyph_source = gs.glyph_list[1]
e.add_module(v)
v.implicit_plane.set(normal=(0, 1, 0), origin=(0, 3, 0))
v = VectorCutPlane()
glyph = v.glyph
gs = glyph.glyph_source
gs.glyph_source = gs.glyph_list[2]
gs.glyph_position = 'head'
e.add_module(v)
v.implicit_plane.set(normal=(0, 1, 0), origin=(0, -2, 0))
self.g=g
self.v=v
self.scene = e.current_scene
return
def add_cortical_surf(self):
# lifted from Gael Varoquax
# this is fairly brittle-- don't know what will result if
# the brain is not skull-stripped
# brain_image is (currently) a copy of the integer indices
# (not the "real valued" scalars)
brain_image = self.master_src.blender.main.image_arr.copy()
np.putmask(brain_image, brain_image>255, 0)
# since this is skull stripped from BET, the boundary will be
# wherever the image drops off to 0
arr = ndimage.gaussian_filter(
(brain_image > 0).astype('d'), 6
)
mask = ndimage.binary_fill_holes(arr > .5)
# iterations x voxel size ~= erosion depth??
mask = ndimage.binary_erosion(mask, iterations=5)
arr_blurred = ndimage.gaussian_filter(
mask.astype('d'), 2
)
# This AA filter will rather dangle off the pipeline..
# it will be the "source" of the ProbeFilter, but we want
# to keep it pipeline-enabled so that we can dynamically
# change colors
point_scalars = surf_to_component[self.surface_component]
surf_colors = mlab.pipeline.set_active_attribute(
self.master_src, point_scalars=point_scalars
)
# Now, make a new ArraySource with the attributes copied
# from master_src
anat_blurred = ArraySource(transpose_input_array=False)
anat_blurred.scalar_data = arr_blurred
anat_blurred.scalar_name = 'blurred'
anat_blurred.origin = self.master_src.data.origin
anat_blurred.spacing = self.master_src.data.spacing
anat_blurred = mlab.pipeline.add_dataset(anat_blurred)
## anat_blurred.update_pipeline()
contour = mlab.pipeline.contour(anat_blurred)
decimated = mlab.pipeline.decimate_pro(contour)
extracted = mlab.pipeline.user_defined(
decimated, filter=self.poly_extractor
)
sampler = tvtk.ProbeFilter()
sampler.source = surf_colors.outputs[0]
surf_points = mlab.pipeline.user_defined(extracted, filter=sampler)
self.cortical_surf = mlab.pipeline.surface(
surf_points,
opacity=.95
)
self.cortical_surf.actor.property.backface_culling = True
self.bcontour = contour
self.surf_colors = surf_colors
# 'mayavi' is always defined on the interpreter. e = OffScreenEngine() #e = Engine() e.start() win = e.new_scene(magnification=1) win.scene.isometric_view() # Make the data and add it to the pipeline. src = ArraySource(transpose_input_array=True) src.scalar_data = data.buf src.spacing=(data.dx, data.dy, -data.dz) src.origin=(data.x_orig, data.y_orig, -data.z_orig) #print src.all_trait_names() # e.add_source(src) # Visualize the data. # o = Outline() # e.add_module(o) # lut=e.scenes[0].children[0].children[0].scalar_lut_manager # lut.data_range=[-1,600] # lut.show_legend = True # lut.data_name = 'Stack' # Create one ContourGridPlane normal to the 'x' axis. cgp = ContourGridPlane()
def do(self):
############################################################
# Imports.
script = self.script
from enthought.mayavi.sources.array_source import ArraySource
from enthought.mayavi.modules.outline import Outline
from enthought.mayavi.modules.glyph import Glyph
from enthought.mayavi.modules.vector_cut_plane import VectorCutPlane
############################################################
# Create a new scene and set up the visualization.
s = self.new_scene()
d = ArraySource()
sc, vec = self.make_data()
d.origin = (-5, -5, -5)
d.scalar_data = sc
d.vector_data = vec
script.add_source(d)
# Create an outline for the data.
o = Outline()
script.add_module(o)
# Glyphs for the scalars
g = Glyph()
script.add_module(g)
g.glyph.glyph_source.glyph_position = 'center'
g.glyph.glyph.vector_mode = 'use_normal'
g.glyph.glyph.scale_factor = 0.5
g.actor.property.line_width = 1.0
v = VectorCutPlane()
glyph = v.glyph
gs = glyph.glyph_source
gs.glyph_position = 'tail'
gs.glyph_source = gs.glyph_list[1]
script.add_module(v)
v.implicit_plane.set(normal=(0, 1, 0), origin=(0, 3, 0))
v = VectorCutPlane()
glyph = v.glyph
gs = glyph.glyph_source
gs.glyph_source = gs.glyph_list[2]
gs.glyph_position = 'head'
script.add_module(v)
v.implicit_plane.set(normal=(0, 1, 0), origin=(0, -2, 0))
# Set the scene to a suitable view.
self.set_view(s)
self.check()
############################################################
# Test if the modules respond correctly when the components
# are changed.
g.actor = g.actor.__class__()
glyph = g.glyph
g.glyph = glyph.__class__()
g.glyph = glyph
glyph = v.glyph
v.glyph = glyph.__class__()
v.glyph = glyph
v.actor = v.actor.__class__()
v.cutter = v.cutter.__class__()
ip = v.implicit_plane
v.implicit_plane = ip.__class__()
v.implicit_plane = ip
s.render()
self.check()
############################################################
# Test if saving a visualization and restoring it works.
# Save visualization.
f = StringIO()
f.name = abspath('test.mv2') # We simulate a file.
script.save_visualization(f)
f.seek(0) # So we can read this saved data.
# Remove existing scene.
engine = script.engine
engine.close_scene(s)
# Load visualization
script.load_visualization(f)
s = engine.current_scene
# Set the scene to a suitable view.
self.set_view(s)
self.check()
############################################################
# Test if the MayaVi2 visualization can be deepcopied.
# Pop the source object.
sources = s.children
s.children = []
# Add it back to see if that works without error.
s.children.extend(sources)
self.set_view(s)
self.check()
# Now deepcopy the source and replace the existing one with
# the copy. This basically simulates cutting/copying the
# object from the UI via the right-click menu on the tree
# view, and pasting the copy back.
sources1 = copy.deepcopy(sources)
s.children[:] = sources
self.set_view(s)
self.check()
def do(self):
############################################################
# Imports.
script = self.script
from enthought.mayavi.sources.array_source import ArraySource
from enthought.mayavi.modules.outline import Outline
from enthought.mayavi.modules.glyph import Glyph
from enthought.mayavi.modules.vector_cut_plane import VectorCutPlane
############################################################
# Create a new scene and set up the visualization.
s = self.new_scene()
d = ArraySource()
sc, vec = self.make_data()
d.origin = (-5, -5, -5)
d.scalar_data = sc
d.vector_data = vec
script.add_source(d)
# Create an outline for the data.
o = Outline()
script.add_module(o)
# Glyphs for the scalars
g = Glyph()
script.add_module(g)
g.glyph.glyph_source.glyph_position = 'center'
g.glyph.glyph.vector_mode = 'use_normal'
g.glyph.glyph.scale_factor = 0.5
g.actor.property.line_width = 1.0
v = VectorCutPlane()
glyph = v.glyph
gs = glyph.glyph_source
gs.glyph_position = 'tail'
gs.glyph_source = gs.glyph_list[1]
script.add_module(v)
v.implicit_plane.set(normal=(0, 1, 0), origin=(0, 3, 0))
v = VectorCutPlane()
glyph = v.glyph
gs = glyph.glyph_source
gs.glyph_source = gs.glyph_list[2]
gs.glyph_position = 'head'
script.add_module(v)
v.implicit_plane.set(normal=(0, 1, 0), origin=(0, -2, 0))
# Set the scene to a suitable view.
self.set_view(s)
self.check()
############################################################
# Test if the modules respond correctly when the components
# are changed.
g.actor = g.actor.__class__()
glyph = g.glyph
g.glyph = glyph.__class__()
g.glyph = glyph
glyph = v.glyph
v.glyph = glyph.__class__()
v.glyph = glyph
v.actor = v.actor.__class__()
v.cutter = v.cutter.__class__()
ip = v.implicit_plane
v.implicit_plane = ip.__class__()
v.implicit_plane = ip
s.render()
self.check()
############################################################
# Test if saving a visualization and restoring it works.
# Save visualization.
f = StringIO()
f.name = abspath('test.mv2') # We simulate a file.
script.save_visualization(f)
f.seek(0) # So we can read this saved data.
# Remove existing scene.
engine = script.engine
engine.close_scene(s)
# Load visualization
script.load_visualization(f)
s = engine.current_scene
# Set the scene to a suitable view.
self.set_view(s)
self.check()
############################################################
# Test if the MayaVi2 visualization can be deepcopied.
# Pop the source object.
sources = s.children
s.children = []
# Add it back to see if that works without error.
s.children.extend(sources)
self.set_view(s)
self.check()
# Now deepcopy the source and replace the existing one with
# the copy. This basically simulates cutting/copying the
# object from the UI via the right-click menu on the tree
# view, and pasting the copy back.
sources1 = copy.deepcopy(sources)
s.children[:] = sources
self.set_view(s)
self.check()
data.buf.shape = (data.nx, data.ny, data.nz)
max_ib = numpy.argmax(data.buf)
ix, iy, iz = data.get_ix_iy_iz(max_ib)
# 'mayavi' is always defined on the interpreter.
e = OffScreenEngine()
#e = Engine()
e.start()
win = e.new_scene(magnification=1)
win.scene.isometric_view()
# Make the data and add it to the pipeline.
src = ArraySource(transpose_input_array=True)
src.scalar_data = data.buf
src.spacing = (data.dx, data.dy, -data.dz)
src.origin = (data.x_orig, data.y_orig, -data.z_orig)
#print src.all_trait_names()
# e.add_source(src)
# Visualize the data.
# o = Outline()
# e.add_module(o)
# lut=e.scenes[0].children[0].children[0].scalar_lut_manager
# lut.data_range=[-1,600]
# lut.show_legend = True
# lut.data_name = 'Stack'
# Create one ContourGridPlane normal to the 'x' axis.
cgp = ContourGridPlane()
# e.add_module(cgp)
def do(self):
############################################################
# Imports.
script = self.script
from enthought.mayavi.sources.array_source import ArraySource
from enthought.mayavi.modules.outline import Outline
from enthought.mayavi.modules.streamline import Streamline
############################################################
# Create a new scene and set up the visualization.
s = self.new_scene()
d = ArraySource()
sc, vec = self.make_data()
d.origin = (-5, -5, -5)
d.scalar_data = sc
d.vector_data = vec
script.add_source(d)
# Create an outline for the data.
o = Outline()
script.add_module(o)
# View the data.
st = Streamline()
script.add_module(st)
widget = st.seed.widget
widget.set(radius=1.0,
center=(-4.0, -4.0, -4.0),
theta_resolution=4,
phi_resolution=4)
st = Streamline(streamline_type='ribbon')
seed = st.seed
seed.widget = seed.widget_list[1]
script.add_module(st)
seed.widget.set(point1=(-5.0, -4.5, -4.0), point2=(-5.0, -4.5, 4.0))
st.ribbon_filter.width = 0.25
st = Streamline(streamline_type='tube')
seed = st.seed
seed.widget = seed.widget_list[2]
script.add_module(st)
seed.widget.set(center=(-5.0, 1.5, -2.5))
st.tube_filter.radius = 0.15
st = Streamline(streamline_type='tube')
seed = st.seed
seed.widget = seed.widget_list[3]
script.add_module(st)
seed.widget.position = (-5.0, 3.75, 3.75)
st.tube_filter.radius = 0.2
# Set the scene to a suitable view.
s.scene.z_plus_view()
c = s.scene.camera
c.azimuth(-30)
c.elevation(30)
s.render()
# Now compare the image.
self.compare_image(s, 'images/test_streamline.png')
############################################################
# Test if the modules respond correctly when the components
# are changed.
tf = st.tube_filter
st.tube_filter = tf.__class__()
st.tube_filter = tf
st.ribbon_filter = st.ribbon_filter.__class__()
seed = st.seed
st.seed = seed.__class__()
st.seed = seed
st.actor = st.actor.__class__()
tracer = st.stream_tracer
st.stream_tracer = tracer.__class__()
st.stream_tracer = tracer
s.render()
# Now compare the image.
self.compare_image(s, 'images/test_streamline.png')
s.render()
############################################################
# Test if saving a visualization and restoring it works.
bg = s.scene.background
# Save visualization.
f = StringIO()
f.name = abspath('test.mv2') # We simulate a file.
script.save_visualization(f)
f.seek(0) # So we can read this saved data.
# Remove existing scene.
engine = script.engine
engine.close_scene(s)
# Load visualization
script.load_visualization(f)
s = engine.current_scene
# Set the scene to a suitable view.
s.scene.z_plus_view()
c = s.scene.camera
c.azimuth(-30)
c.elevation(30)
s.render()
s.scene.background = bg
# Now compare the image.
self.compare_image(s, 'images/test_streamline.png')
############################################################
# Test if the MayaVi2 visualization can be deepcopied.
# Pop the source object.
sources = s.children
s.children = []
# Add it back to see if that works without error.
s.children.extend(sources)
s.scene.reset_zoom()
# Now compare the image.
self.compare_image(s, 'images/test_streamline.png')
# Now deepcopy the source and replace the existing one with
# the copy. This basically simulates cutting/copying the
# object from the UI via the right-click menu on the tree
# view, and pasting the copy back.
sources1 = copy.deepcopy(sources)
s.children[:] = sources
s.scene.reset_zoom()
self.compare_image(s, 'images/test_streamline.png')
def plot_slice_mayavi(dat_filename,output_file,hyp_x,hyp_y,hyp_z,search_grid_file_name,max_stack_value):
base_path=os.getenv('WAVELOC_PATH')
lib_path="%s/lib"%base_path
# grid geometry
hdr_file=lib_path + os.sep + search_grid_file_name
# detection
detection=50
# stations
stations_file="%s/coord_stations_piton"%lib_path
sta=StationList()
sta.read_from_file(stations_file)
# create the object to contain the stations
pd = tvtk.PolyData()
pd.points = [[s.x/1000.0, s.y/1000.0, -s.elev/1000.0] for s in sta.stations.values()]
# create the object to contain the stations
try:
pd_hyp = tvtk.PolyData()
pd_hyp.points=[[hyp_x,hyp_y,hyp_z]]
except TypeError:
pass
# read the dat file
print dat_filename
data=QDGrid()
data.read_NLL_hdr_file(hdr_file)
data.buf=numpy.fromfile(dat_filename, dtype=numpy.int16)
max_ib=numpy.argmax(data.buf)
print max_ib
max_val=data.buf[max_ib]
ix,iy,iz=data.get_ix_iy_iz(max_ib)
#data.buf=numpy.array(data.buf, dtype=numpy.float)
data.buf.shape = (data.nx, data.ny, data.nz)
# 'mayavi' is always defined on the interpreter.
e = OffScreenEngine()
#e = Engine()
e.start()
win = e.new_scene(magnification=1)
e.current_scene.scene.off_screen_rendering = True
win.scene.isometric_view()
# Make the data and add it to the pipeline.
src = ArraySource(transpose_input_array=True)
src.scalar_data = data.buf
src.spacing=(data.dx, data.dy, -data.dz)
src.origin=(data.x_orig, data.y_orig, -data.z_orig)
e.add_source(src)
# Visualize the data.
o = Outline()
e.add_module(o)
lut=e.scenes[0].children[0].children[0].scalar_lut_manager
lut.data_range=[-1,max_stack_value]
lut.show_legend = True
lut.data_name = 'Stack'
# Create one ContourGridPlane normal to the 'x' axis.
cgp = ContourGridPlane()
e.add_module(cgp)
# Set the position to the middle of the data.
if max_val > detection:
cgp.grid_plane.position = ix
else:
cgp.grid_plane.position = data.nx/2
cgp.contour.filled_contours = True
cgp.actor.property.opacity = 0.6
output=cgp.grid_plane.outputs[0]
x_data=numpy.array(output.point_data.scalars.to_array())
# Another with filled contours normal to 'y' axis.
cgp = ContourGridPlane()
e.add_module(cgp)
# Set the axis and position to the middle of the data.
cgp.grid_plane.axis = 'y'
if max_val > detection:
cgp.grid_plane.position = iy
else:
cgp.grid_plane.position = data.ny/2
cgp.contour.filled_contours = True
cgp.actor.property.opacity = 0.6
output=cgp.grid_plane.outputs[0]
y_data=numpy.array(output.point_data.scalars.to_array())
# Another with filled contours normal to 'z' axis.
cgp = ContourGridPlane()
e.add_module(cgp)
# Set the axis and position to the middle of the data.
cgp.grid_plane.axis = 'z'
if max_val > detection:
cgp.grid_plane.position = iz
else:
cgp.grid_plane.position = data.nz/2
cgp.contour.filled_contours = True
cgp.actor.property.opacity = 0.6
output=cgp.grid_plane.outputs[0]
z_data=numpy.array(output.point_data.scalars.to_array())
a=Axes()
e.add_module(a)
d=VTKDataSource()
d.data=pd
e.add_source(d)
g=Glyph()
e.add_module(g)
g.glyph.glyph_source.glyph_source=g.glyph.glyph_source.glyph_list[4]
g.glyph.glyph_source.glyph_source.radius=0.1
d=VTKDataSource()
d.data=pd_hyp
e.add_source(d)
g=Glyph()
e.add_module(g)
g.glyph.glyph_source.glyph_source=g.glyph.glyph_source.glyph_list[4]
g.glyph.glyph_source.glyph_source.radius=0.5
g.actor.property.color=(0.0,0.0,0.0)
#view(azimuth=-60,elevation=60,distance=120)
win.scene.save(output_file,size=(800,800))
e.stop()
del win
del e
return (x_data, y_data, z_data)
def do(self):
############################################################
# Imports.
script = self.script
from enthought.mayavi.sources.array_source import ArraySource
from enthought.mayavi.modules.outline import Outline
from enthought.mayavi.modules.surface import Surface
from enthought.mayavi.modules.vectors import Vectors
############################################################
# Create a new scene and set up the visualization.
s = self.new_scene()
d = ArraySource()
self.check_input_validation(d)
sc, vec = self.make_2d_data()
d.origin = (-1, -1, 0)
d.scalar_data = sc
d.vector_data = vec
script.add_source(d)
# Create an outline for the data.
o = Outline()
script.add_module(o)
# View the data.
s = Surface()
script.add_module(s)
v = Vectors()
script.add_module(v)
# Add a 3D data source
d = ArraySource()
sc, vec = self.make_3d_data()
d.scalar_data = sc
d.vector_data = vec
script.add_source(d)
# Create an outline for the data.
o = Outline()
script.add_module(o)
# View a slice.
s = Surface()
script.add_module(s)
v = Vectors()
script.add_module(v)
# Set the scene to a suitable view.
s.scene.z_plus_view()
c = s.scene.camera
c.azimuth(-30)
c.elevation(30)
self.check()
############################################################
# Test if saving a visualization and restoring it works.
bg = s.scene.background
# Save visualization.
f = StringIO()
f.name = abspath('test.mv2') # We simulate a file.
script.save_visualization(f)
f.seek(0) # So we can read this saved data.
# Remove existing scene.
engine = script.engine
engine.close_scene(s)
# Load visualization
script.load_visualization(f)
s = engine.current_scene
# Set the scene to a suitable view.
s.scene.z_plus_view()
c = s.scene.camera
c.azimuth(-30)
c.elevation(30)
s.scene.background = bg
self.check()
############################################################
# Test if the MayaVi2 visualization can be deepcopied.
# Pop the source object.
sources = s.children
s.children = []
# Add it back to see if that works without error.
s.children.extend(sources)
s.scene.reset_zoom()
self.check()
# Now deepcopy the source and replace the existing one with
# the copy. This basically simulates cutting/copying the
# object from the UI via the right-click menu on the tree
# view, and pasting the copy back.
sources1 = copy.deepcopy(sources)
s.children[:] = sources
s.scene.reset_zoom()
self.check()
def plot_slice_mayavi(dat_filename, output_file, hyp_x, hyp_y, hyp_z,
search_grid_file_name, max_stack_value):
base_path = os.getenv('WAVELOC_PATH')
lib_path = "%s/lib" % base_path
# grid geometry
hdr_file = lib_path + os.sep + search_grid_file_name
# detection
detection = 50
# stations
stations_file = "%s/coord_stations_piton" % lib_path
sta = StationList()
sta.read_from_file(stations_file)
# create the object to contain the stations
pd = tvtk.PolyData()
pd.points = [[s.x / 1000.0, s.y / 1000.0, -s.elev / 1000.0]
for s in sta.stations.values()]
# create the object to contain the stations
try:
pd_hyp = tvtk.PolyData()
pd_hyp.points = [[hyp_x, hyp_y, hyp_z]]
except TypeError:
pass
# read the dat file
print dat_filename
data = QDGrid()
data.read_NLL_hdr_file(hdr_file)
data.buf = numpy.fromfile(dat_filename, dtype=numpy.int16)
max_ib = numpy.argmax(data.buf)
print max_ib
max_val = data.buf[max_ib]
ix, iy, iz = data.get_ix_iy_iz(max_ib)
#data.buf=numpy.array(data.buf, dtype=numpy.float)
data.buf.shape = (data.nx, data.ny, data.nz)
# 'mayavi' is always defined on the interpreter.
e = OffScreenEngine()
#e = Engine()
e.start()
win = e.new_scene(magnification=1)
e.current_scene.scene.off_screen_rendering = True
win.scene.isometric_view()
# Make the data and add it to the pipeline.
src = ArraySource(transpose_input_array=True)
src.scalar_data = data.buf
src.spacing = (data.dx, data.dy, -data.dz)
src.origin = (data.x_orig, data.y_orig, -data.z_orig)
e.add_source(src)
# Visualize the data.
o = Outline()
e.add_module(o)
lut = e.scenes[0].children[0].children[0].scalar_lut_manager
lut.data_range = [-1, max_stack_value]
lut.show_legend = True
lut.data_name = 'Stack'
# Create one ContourGridPlane normal to the 'x' axis.
cgp = ContourGridPlane()
e.add_module(cgp)
# Set the position to the middle of the data.
if max_val > detection:
cgp.grid_plane.position = ix
else:
cgp.grid_plane.position = data.nx / 2
cgp.contour.filled_contours = True
cgp.actor.property.opacity = 0.6
output = cgp.grid_plane.outputs[0]
x_data = numpy.array(output.point_data.scalars.to_array())
# Another with filled contours normal to 'y' axis.
cgp = ContourGridPlane()
e.add_module(cgp)
# Set the axis and position to the middle of the data.
cgp.grid_plane.axis = 'y'
if max_val > detection:
cgp.grid_plane.position = iy
else:
cgp.grid_plane.position = data.ny / 2
cgp.contour.filled_contours = True
cgp.actor.property.opacity = 0.6
output = cgp.grid_plane.outputs[0]
y_data = numpy.array(output.point_data.scalars.to_array())
# Another with filled contours normal to 'z' axis.
cgp = ContourGridPlane()
e.add_module(cgp)
# Set the axis and position to the middle of the data.
cgp.grid_plane.axis = 'z'
if max_val > detection:
cgp.grid_plane.position = iz
else:
cgp.grid_plane.position = data.nz / 2
cgp.contour.filled_contours = True
cgp.actor.property.opacity = 0.6
output = cgp.grid_plane.outputs[0]
z_data = numpy.array(output.point_data.scalars.to_array())
a = Axes()
e.add_module(a)
d = VTKDataSource()
d.data = pd
e.add_source(d)
g = Glyph()
e.add_module(g)
g.glyph.glyph_source.glyph_source = g.glyph.glyph_source.glyph_list[4]
g.glyph.glyph_source.glyph_source.radius = 0.1
d = VTKDataSource()
d.data = pd_hyp
e.add_source(d)
g = Glyph()
e.add_module(g)
g.glyph.glyph_source.glyph_source = g.glyph.glyph_source.glyph_list[4]
g.glyph.glyph_source.glyph_source.radius = 0.5
g.actor.property.color = (0.0, 0.0, 0.0)
#view(azimuth=-60,elevation=60,distance=120)
win.scene.save(output_file, size=(800, 800))
e.stop()
del win
del e
return (x_data, y_data, z_data)
def do(self):
############################################################
# Imports.
script = self.script
from enthought.mayavi.sources.array_source import ArraySource
from enthought.mayavi.modules.outline import Outline
from enthought.mayavi.modules.streamline import Streamline
############################################################
# Create a new scene and set up the visualization.
s = self.new_scene()
d = ArraySource()
sc, vec = self.make_data()
d.origin = (-5, -5, -5)
d.scalar_data = sc
d.vector_data = vec
script.add_source(d)
# Create an outline for the data.
o = Outline()
script.add_module(o)
# View the data.
st = Streamline()
script.add_module(st)
widget = st.seed.widget
widget.set(radius=1.0, center=(-4.0, -4.0, -4.0),
theta_resolution=4, phi_resolution=4)
st = Streamline(streamline_type='ribbon')
seed = st.seed
seed.widget = seed.widget_list[1]
script.add_module(st)
seed.widget.set(point1=(-5.0, -4.5, -4.0), point2=(-5.0, -4.5, 4.0))
st.ribbon_filter.width = 0.25
st = Streamline(streamline_type='tube')
seed = st.seed
seed.widget = seed.widget_list[2]
script.add_module(st)
seed.widget.set(center=(-5.0, 1.5, -2.5))
st.tube_filter.radius = 0.15
st = Streamline(streamline_type='tube')
seed = st.seed
seed.widget = seed.widget_list[3]
script.add_module(st)
seed.widget.position=(-5.0, 3.75, 3.75)
st.tube_filter.radius = 0.2
# Set the scene to a suitable view.
s.scene.z_plus_view()
c = s.scene.camera
c.azimuth(-30)
c.elevation(30)
s.render()
# Now compare the image.
self.compare_image(s, 'images/test_streamline.png')
############################################################
# Test if the modules respond correctly when the components
# are changed.
tf = st.tube_filter
st.tube_filter = tf.__class__()
st.tube_filter = tf
st.ribbon_filter = st.ribbon_filter.__class__()
seed = st.seed
st.seed = seed.__class__()
st.seed = seed
st.actor = st.actor.__class__()
tracer = st.stream_tracer
st.stream_tracer = tracer.__class__()
st.stream_tracer = tracer
s.render()
# Now compare the image.
self.compare_image(s, 'images/test_streamline.png')
s.render()
############################################################
# Test if saving a visualization and restoring it works.
bg = s.scene.background
# Save visualization.
f = StringIO()
f.name = abspath('test.mv2') # We simulate a file.
script.save_visualization(f)
f.seek(0) # So we can read this saved data.
# Remove existing scene.
engine = script.engine
engine.close_scene(s)
# Load visualization
script.load_visualization(f)
s = engine.current_scene
# Set the scene to a suitable view.
s.scene.z_plus_view()
c = s.scene.camera
c.azimuth(-30)
c.elevation(30)
s.render()
s.scene.background = bg
# Now compare the image.
self.compare_image(s, 'images/test_streamline.png')
############################################################
# Test if the MayaVi2 visualization can be deepcopied.
# Pop the source object.
sources = s.children
s.children = []
# Add it back to see if that works without error.
s.children.extend(sources)
s.scene.reset_zoom()
# Now compare the image.
self.compare_image(s, 'images/test_streamline.png')
# Now deepcopy the source and replace the existing one with
# the copy. This basically simulates cutting/copying the
# object from the UI via the right-click menu on the tree
# view, and pasting the copy back.
sources1 = copy.deepcopy(sources)
s.children[:] = sources
s.scene.reset_zoom()
self.compare_image(s, 'images/test_streamline.png')
