<br>Hello Doug,<br><br>Your points are being rendered white, which happens to be your background color. If you change your background color to something dark, you will see the points of your mesh. You could also extract a slice along each dimension of the mesh, which might give you something more along the lines of what you're looking for:
<br><br>#!/usr/bin/env python<br>import vtk<br>import math<br>from vtk.util.colors import *<br><br>rMin=0.5<br>rMax=1.0<br>dims = [13,11,11]<br>size = dims[0]*dims[1]*dims[2]<br><br>points = vtk.vtkPoints()<br>points.Allocate
(size,size)<br><br>deltaZ = 2.0 / (dims[2]-1)<br>deltaRad = (rMax-rMin) / (dims[1]-1)<br>x = [0.0, 0.0, 0.0]<br>for k in range(dims[2]):<br> x[2] = -1.0 + k*deltaZ<br> kOffset = k * dims[0] * dims[1]<br> for j in range(dims[1]):
<br> radius = rMin + j*deltaRad<br> jOffset = j * dims[0]<br> for i in range(dims[0]):<br> theta = i * 15.0 * vtk.vtkMath.DegreesToRadians()<br> x[0] = radius * math.cos(theta)<br> x[1] = radius *
math.sin(theta)<br> offset = i + jOffset + kOffset<br> points.InsertPoint(offset,x)<br><br> # Create the structured grid.<br>sgrid = vtk.vtkStructuredGrid()<br>sgrid.SetDimensions(dims) # Set dimensions of grid
<br>sgrid.SetPoints(points) # Load data for grid point locations<br><br><br>outline = vtk.vtkOutlineFilter()<br>outline.SetInput(sgrid)<br><br>outline_mapper = vtk.vtkPolyDataMapper()<br>outline_mapper.SetInputConnection(
outline.GetOutputPort())<br><br>outline_actor = vtk.vtkActor()<br>outline_actor.SetMapper(outline_mapper)<br><br>def createMeshActor(imin, imax, jmin, jmax, kmin, kmax):<br><br> mesh = vtk.vtkStructuredGridGeometryFilter
()<br> mesh.SetInput(sgrid)<br> mesh.SetExtent(imin, imax, jmin, jmax, kmin, kmax)<br><br> meshMapper = vtk.vtkPolyDataMapper()<br> meshMapper.SetInputConnection(mesh.GetOutputPort())<br> meshMapper.ScalarVisibilityOff
()<br><br> meshActor = vtk.vtkActor()<br> meshActor.SetMapper(meshMapper)<br> meshActor.GetProperty().SetRepresentationToWireframe()<br> <br> return meshActor<br> <br> # Create the usual rendering stuff<br>
n = 0<br>renderer = vtk.vtkRenderer()<br>renderer.AddActor(createMeshActor(n,n,0,dims[1],0,dims[2]))<br>renderer.AddActor(createMeshActor(0,dims[0],n,n,0,dims[2]))<br>renderer.AddActor(createMeshActor(0,dims[0],0,dims[1],n,n))
<br><br>renWin = vtk.vtkRenderWindow()<br>renWin.AddRenderer(renderer)<br><br>iren = vtk.vtkRenderWindowInteractor()<br>iren.SetRenderWindow(renWin)<br><br>#renderer.AddActor(outline_actor)<br><br>renderer.SetBackground(0.25
,0.25,0.5)<br>renderer.ResetCamera()<br>renderer.GetActiveCamera().Elevation(60.0)<br>renderer.GetActiveCamera().Azimuth(30.0)<br>renderer.GetActiveCamera().Zoom(1.25)<br>renWin.SetSize(600,600)<br><br> # interact with data
<br>renWin.Render()<br>iren.Start()<br><br><br clear="all"><br>-- <br>Ed Bachta<br>Visualization and Interactive Spaces<br>Pervasive Tech Labs at Indiana University