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 95 96 97 98 99 100 101 102 103 104 105 106 107 108 109 110 111 112 113 114 115 116 117 118 119 120 121 122 123 124 125 126 127 128 129 130 131 132
|
package require vtk
package require vtkinteraction
# create pipeline - structured grid
#
vtkMultiBlockPLOT3DReader pl3d
pl3d SetXYZFileName "$VTK_DATA_ROOT/Data/combxyz.bin"
pl3d SetQFileName "$VTK_DATA_ROOT/Data/combq.bin"
pl3d SetScalarFunctionNumber 100
pl3d SetVectorFunctionNumber 202
pl3d Update
set output [[pl3d GetOutput] GetBlock 0]
vtkDataSetSurfaceFilter gf
gf SetInputData $output
vtkPolyDataMapper gMapper
gMapper SetInputConnection [gf GetOutputPort]
eval gMapper SetScalarRange [$output GetScalarRange]
vtkActor gActor
gActor SetMapper gMapper
vtkDataSetSurfaceFilter gf2
gf2 SetInputData $output
vtkPolyDataMapper g2Mapper
g2Mapper SetInputConnection [gf2 GetOutputPort]
eval g2Mapper SetScalarRange [$output GetScalarRange]
vtkActor g2Actor
g2Actor SetMapper g2Mapper
g2Actor AddPosition 0 15 0
# create pipeline - poly data
#
vtkDataSetSurfaceFilter gf3
gf3 SetInputConnection [gf GetOutputPort]
vtkPolyDataMapper g3Mapper
g3Mapper SetInputConnection [gf3 GetOutputPort]
eval g3Mapper SetScalarRange [$output GetScalarRange]
vtkActor g3Actor
g3Actor SetMapper g3Mapper
g3Actor AddPosition 0 0 15
vtkDataSetSurfaceFilter gf4
gf4 SetInputConnection [gf2 GetOutputPort]
gf4 UseStripsOn
vtkPolyDataMapper g4Mapper
g4Mapper SetInputConnection [gf4 GetOutputPort]
eval g4Mapper SetScalarRange [$output GetScalarRange]
vtkActor g4Actor
g4Actor SetMapper g4Mapper
g4Actor AddPosition 0 15 15
# create pipeline - unstructured grid
#
vtkSphere s
eval s SetCenter [$output GetCenter]
s SetRadius 100.0; #everything
vtkExtractGeometry eg
eg SetInputData $output
eg SetImplicitFunction s
vtkDataSetSurfaceFilter gf5
gf5 SetInputConnection [eg GetOutputPort]
vtkPolyDataMapper g5Mapper
g5Mapper SetInputConnection [gf5 GetOutputPort]
eval g5Mapper SetScalarRange [$output GetScalarRange]
vtkActor g5Actor
g5Actor SetMapper g5Mapper
g5Actor AddPosition 0 0 30
vtkDataSetSurfaceFilter gf6
gf6 SetInputConnection [eg GetOutputPort]
gf6 UseStripsOn
vtkPolyDataMapper g6Mapper
g6Mapper SetInputConnection [gf6 GetOutputPort]
eval g6Mapper SetScalarRange [$output GetScalarRange]
vtkActor g6Actor
g6Actor SetMapper g6Mapper
g6Actor AddPosition 0 15 30
# create pipeline - rectilinear grid
#
vtkRectilinearGridReader rgridReader
rgridReader SetFileName "$VTK_DATA_ROOT/Data/RectGrid2.vtk"
rgridReader Update
vtkDataSetSurfaceFilter gf7
gf7 SetInputConnection [rgridReader GetOutputPort]
vtkPolyDataMapper g7Mapper
g7Mapper SetInputConnection [gf7 GetOutputPort]
eval g7Mapper SetScalarRange [[rgridReader GetOutput] GetScalarRange]
vtkActor g7Actor
g7Actor SetMapper g7Mapper
g7Actor SetScale 3 3 3
vtkDataSetSurfaceFilter gf8
gf8 SetInputConnection [rgridReader GetOutputPort]
gf8 UseStripsOn
vtkPolyDataMapper g8Mapper
g8Mapper SetInputConnection [gf8 GetOutputPort]
eval g8Mapper SetScalarRange [[rgridReader GetOutput] GetScalarRange]
vtkActor g8Actor
g8Actor SetMapper g8Mapper
g8Actor SetScale 3 3 3
g8Actor AddPosition 0 15 0
# Create the RenderWindow, Renderer and both Actors
vtkRenderer ren1
vtkRenderWindow renWin
renWin AddRenderer ren1
vtkRenderWindowInteractor iren
iren SetRenderWindow renWin
ren1 AddActor gActor
ren1 AddActor g2Actor
ren1 AddActor g3Actor
ren1 AddActor g4Actor
ren1 AddActor g5Actor
ren1 AddActor g6Actor
ren1 AddActor g7Actor
ren1 AddActor g8Actor
renWin SetSize 340 550
set cam1 [ren1 GetActiveCamera]
$cam1 SetClippingRange 84 174
$cam1 SetFocalPoint 5.22824 6.09412 35.9813
$cam1 SetPosition 100.052 62.875 102.818
$cam1 SetViewUp -0.307455 -0.464269 0.830617
iren Initialize
# render the image
#
iren AddObserver UserEvent {wm deiconify .vtkInteract}
# prevent the tk window from showing up then start the event loop
wm withdraw .
|