1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94
|
package require vtk
package require vtkinteraction
package require vtktesting
# Create the RenderWindow, Renderer and both Actors
#
vtkRenderer ren1
vtkRenderWindow renWin
renWin AddRenderer ren1
vtkRenderWindowInteractor iren
iren SetRenderWindow renWin
# create pipeline
#
vtkVolume16Reader v16
v16 SetDataDimensions 64 64
[v16 GetOutput] SetOrigin 0.0 0.0 0.0
v16 SetDataByteOrderToLittleEndian
v16 SetFilePrefix "$VTK_DATA_ROOT/Data/headsq/quarter"
v16 SetImageRange 1 93
v16 SetDataSpacing 3.2 3.2 1.5
v16 Update
vtkImageMarchingCubes iso
iso SetInputConnection [v16 GetOutputPort]
iso SetValue 0 1150
iso SetInputMemoryLimit 1000
vtkPlane topPlane
topPlane SetNormal 0 0 1
topPlane SetOrigin 0 0 0.5
vtkPlane botPlane
botPlane SetNormal 0 0 -1
botPlane SetOrigin 0 0 137.0
vtkPlane sagPlane
sagPlane SetNormal 1 0 0
sagPlane SetOrigin 100.8 0 0
vtkPlaneCollection capPlanes
capPlanes AddItem topPlane
capPlanes AddItem botPlane
capPlanes AddItem sagPlane
vtkClipClosedSurface clip
clip SetClippingPlanes capPlanes
clip SetInputConnection [iso GetOutputPort]
clip SetBaseColor 0.9804 0.9216 0.8431
clip SetClipColor 1.0 1.0 1.0
clip SetActivePlaneColor 1.0 1.0 0.8
clip SetActivePlaneId 2
clip SetScalarModeToColors
clip GenerateOutlineOn
clip GenerateFacesOn
vtkPolyDataMapper isoMapper
isoMapper SetInputConnection [clip GetOutputPort]
isoMapper ScalarVisibilityOn
isoMapper ImmediateModeRenderingOn
vtkActor isoActor
isoActor SetMapper isoMapper
vtkOutlineFilter outline
outline SetInputConnection [v16 GetOutputPort]
vtkPolyDataMapper outlineMapper
outlineMapper SetInputConnection [outline GetOutputPort]
vtkActor outlineActor
outlineActor SetMapper outlineMapper
outlineActor VisibilityOff
# Add the actors to the renderer, set the background and size
#
ren1 AddActor outlineActor
ren1 AddActor isoActor
ren1 SetBackground 0.2 0.3 0.4
renWin SetSize 200 200
ren1 ResetCamera
[ren1 GetActiveCamera] Elevation 90
[ren1 GetActiveCamera] SetViewUp 0 0 -1
[ren1 GetActiveCamera] Azimuth 270
ren1 ResetCameraClippingRange
renWin Render
# render the image
#
iren AddObserver UserEvent {wm deiconify .vtkInteract}
# prevent the tk window from showing up then start the event loop
wm withdraw .
|