examples/pyvisi/seriesofcuts.py


########################################################
#
# Copyright (c) 2003-2008 by University of Queensland
# Earth Systems Science Computational Center (ESSCC)
# http://www.uq.edu.au/esscc
#
# Primary Business: Queensland, Australia
# Licensed under the Open Software License version 3.0
# http://www.opensource.org/licenses/osl-3.0.php
#
########################################################

__copyright__="""Copyright (c) 2003-2008 by University of Queensland
Earth Systems Science Computational Center (ESSCC)
http://www.uq.edu.au/esscc
Primary Business: Queensland, Australia"""
__license__="""Licensed under the Open Software License version 3.0
http://www.opensource.org/licenses/osl-3.0.php"""
__url__="http://www.uq.edu.au/esscc/escript-finley"

"""
Author: John Ngui, john.ngui@uq.edu.au
"""

# 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 using 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	ksteube	1811
2			########################################################
3			#
4			# Copyright (c) 2003-2008 by University of Queensland
5			# Earth Systems Science Computational Center (ESSCC)
6			# http://www.uq.edu.au/esscc
7			#
8			# Primary Business: Queensland, Australia
9			# Licensed under the Open Software License version 3.0
10			# http://www.opensource.org/licenses/osl-3.0.php
11			#
12			########################################################
13
14			__copyright__="""Copyright (c) 2003-2008 by University of Queensland
15			Earth Systems Science Computational Center (ESSCC)
16			http://www.uq.edu.au/esscc
17			Primary Business: Queensland, Australia"""
18			__license__="""Licensed under the Open Software License version 3.0
19			http://www.opensource.org/licenses/osl-3.0.php"""
20			__url__="http://www.uq.edu.au/esscc/escript-finley"
21
22	jongui	1203	"""
23			Author: John Ngui, john.ngui@uq.edu.au
24			"""
25
26	ksteube	1147	# Import the necessary modules.
27			from esys.pyvisi import Scene, DataCollector, MapOnPlaneCut, Camera
28			from esys.pyvisi.constant import *
29	jongui	1148	import os
30	ksteube	1147
31	jongui	1148	PYVISI_EXAMPLE_MESHES_PATH = "data_meshes"
32			PYVISI_EXAMPLE_IMAGES_PATH = "data_sample_images"
33	ksteube	1147	X_SIZE = 400
34			Y_SIZE = 400
35
36			SCALAR_FIELD_POINT_DATA = "temperature"
37			FILE_3D = "interior_3D.xml"
38			IMAGE_NAME = "seriesofcuts"
39			JPG_RENDERER = Renderer.ONLINE_JPG
40
41			# Create a Scene.
42			s = Scene(renderer = JPG_RENDERER, num_viewport = 1, x_size = X_SIZE,
43			y_size = Y_SIZE)
44
45			# Create a DataCollector reading from a XML file.
46			dc1 = DataCollector(source = Source.XML)
47	jongui	1148	dc1.setFileName(file_name = os.path.join(PYVISI_EXAMPLE_MESHES_PATH, FILE_3D))
48	ksteube	1147	dc1.setActiveScalar(scalar = SCALAR_FIELD_POINT_DATA)
49
50			# Create a MapOnPlaneCut.
51			mopc1 = MapOnPlaneCut(scene = s, data_collector = dc1,
52			viewport = Viewport.SOUTH_WEST, lut = Lut.COLOR, cell_to_point = False,
53			outline = True)
54			mopc1.setPlaneToYZ(offset = 0.1)
55
56			# Create a Camera.
57	jongui	1148	c1 = Camera(scene = s, viewport = Viewport.SOUTH_WEST)
58	ksteube	1147	c1.isometricView()
59
60	jongui	1210	# Render the object with multiple cuts using a series of translation.
61	ksteube	1147	for i in range(0, 5):
62	jongui	1148	s.render(image_name = os.path.join(PYVISI_EXAMPLE_IMAGES_PATH, IMAGE_NAME +
63			"%02d.jpg") % i)
64	ksteube	1147	mopc1.translate(0.6,0,0)