examples/pyvisi/seriesofcuts.py

# Import the necessary modules.
from esys.pyvisi import Scene, DataCollector, MapOnPlaneCut, Camera
from esys.pyvisi.constant import *
import os

PYVISI_EXAMPLE_MESHES_PATH = "data_meshes"
PYVISI_EXAMPLE_IMAGES_PATH = "data_sample_images"
X_SIZE = 400
Y_SIZE = 400

SCALAR_FIELD_POINT_DATA = "temperature"
FILE_3D = "interior_3D.xml"
IMAGE_NAME = "seriesofcuts"
JPG_RENDERER = Renderer.ONLINE_JPG

# Create a Scene.
s = Scene(renderer = JPG_RENDERER, num_viewport = 1, x_size = X_SIZE, 
        y_size = Y_SIZE)

# Create a DataCollector reading from a XML file.
dc1 = DataCollector(source = Source.XML)
dc1.setFileName(file_name = os.path.join(PYVISI_EXAMPLE_MESHES_PATH, FILE_3D))
dc1.setActiveScalar(scalar = SCALAR_FIELD_POINT_DATA)

# Create a MapOnPlaneCut.
mopc1 = MapOnPlaneCut(scene = s, data_collector = dc1, 
        viewport = Viewport.SOUTH_WEST, lut = Lut.COLOR, cell_to_point = False,
        outline = True)
mopc1.setPlaneToYZ(offset = 0.1)

# Create a Camera.
c1 = Camera(scene = s, viewport = Viewport.SOUTH_WEST)
c1.isometricView()

# Render the object with multiple cuts from a series of translation.
for i in range(0, 5):
    s.render(image_name = os.path.join(PYVISI_EXAMPLE_IMAGES_PATH, IMAGE_NAME +
            "%02d.jpg") % i)
    mopc1.translate(0.6,0,0)
1	# Import the necessary modules.
2	from esys.pyvisi import Scene, DataCollector, MapOnPlaneCut, Camera
3	from esys.pyvisi.constant import *
4	import os
5
6	PYVISI_EXAMPLE_MESHES_PATH = "data_meshes"
7	PYVISI_EXAMPLE_IMAGES_PATH = "data_sample_images"
8	X_SIZE = 400
9	Y_SIZE = 400
10
11	SCALAR_FIELD_POINT_DATA = "temperature"
12	FILE_3D = "interior_3D.xml"
13	IMAGE_NAME = "seriesofcuts"
14	JPG_RENDERER = Renderer.ONLINE_JPG
15
16	# Create a Scene.
17	s = Scene(renderer = JPG_RENDERER, num_viewport = 1, x_size = X_SIZE,
18	y_size = Y_SIZE)
19
20	# Create a DataCollector reading from a XML file.
21	dc1 = DataCollector(source = Source.XML)
22	dc1.setFileName(file_name = os.path.join(PYVISI_EXAMPLE_MESHES_PATH, FILE_3D))
23	dc1.setActiveScalar(scalar = SCALAR_FIELD_POINT_DATA)
24
25	# Create a MapOnPlaneCut.
26	mopc1 = MapOnPlaneCut(scene = s, data_collector = dc1,
27	viewport = Viewport.SOUTH_WEST, lut = Lut.COLOR, cell_to_point = False,
28	outline = True)
29	mopc1.setPlaneToYZ(offset = 0.1)
30
31	# Create a Camera.
32	c1 = Camera(scene = s, viewport = Viewport.SOUTH_WEST)
33	c1.isometricView()
34
35	# Render the object with multiple cuts from a series of translation.
36	for i in range(0, 5):
37	s.render(image_name = os.path.join(PYVISI_EXAMPLE_IMAGES_PATH, IMAGE_NAME +
38	"%02d.jpg") % i)
39	mopc1.translate(0.6,0,0)