/[escript]/trunk/pyvisi/examples/isosurfacePlot.py
ViewVC logotype

Annotation of /trunk/pyvisi/examples/isosurfacePlot.py

Parent Directory Parent Directory | Revision Log Revision Log


Revision 637 - (hide annotations)
Thu Mar 23 10:55:31 2006 UTC (14 years, 3 months ago) by gross
File MIME type: text/x-python
File size: 1962 byte(s)
more copyright statements added
1 cochrane 337 # $Id: isosurfacePlot.py,v 1.4 2005/11/08 08:23:45 paultcochrane Exp $
2    
3     """
4     Example of plotting a set of isosurfaces with pyvisi
5 gross 637
6     @var __author__: name of author
7     @var __license__: licence agreement
8     @var __copyright__: copyrights
9     @var __url__: url entry point on documentation
10     @var __version__: version
11     @var __date__: date of the version
12 cochrane 337 """
13    
14 gross 637 __copyright__=""" Copyright (c) 2006 by ACcESS MNRF
15     http://www.access.edu.au
16     Primary Business: Queensland, Australia"""
17     __license__="""Licensed under the Open Software License version 3.0
18     http://www.opensource.org/licenses/osl-3.0.php"""
19     __author__="Paul Cochrane"
20     __url__="http://www.iservo.edu.au/esys"
21     __version__="$Revision$"
22     __date__="$Date$"
23    
24    
25 cochrane 337 import sys
26     numArgs = len(sys.argv)
27     if numArgs == 1:
28     ren_mod = "vtk"
29     else:
30     ren_mod = sys.argv[1]
31    
32     # set up some data to plot
33     from Numeric import *
34    
35     # example code for how a user would write a script in pyvisi
36 gross 605 from esys.pyvisi import * # base level visualisation stuff
37 cochrane 337 # import the objects to render the scene using the specific renderer
38     if ren_mod == "vtk":
39 gross 605 from esys.pyvisi.renderers.vtk import * # vtk
40 cochrane 337 elif ren_mod == "povray":
41 gross 605 from esys.pyvisi.renderers.povray import * # povray
42 cochrane 337 else:
43     raise ValueError, "Unknown renderer module"
44    
45     # define the scene object
46     # a Scene is a container for all of the kinds of things you want to put
47     # into your plot for instance, images, meshes, arrow/vector/quiver plots,
48     # contour plots, spheres etc.
49     scene = Scene()
50    
51     # create a IsosurfacePlot object
52     plot = IsosurfacePlot(scene)
53    
54     # add some helpful info to the plot
55     plot.title = 'Example isosurface plot'
56     plot.xlabel = 'x'
57     plot.ylabel = 'y'
58     plot.zlabel = 'z'
59    
60     # plot data defined in a vtk file
61     plot.setData(fname='temp-0500.vtk', format='vtk-xml')
62    
63     scene.render(pause=True, interactive=True)
64    
65     # save the plot
66     scene.save(fname="isosurfacePlot.png", format="png")
67    
68     # vim: expandtab shiftwidth=4:
69    

  ViewVC Help
Powered by ViewVC 1.1.26