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
|
package require vtk
package require vtkinteraction
# Quadric definition
vtkQuadric quadric
quadric SetCoefficients .5 1 .2 0 .1 0 0 .2 0 0
vtkSampleFunction sample
sample SetSampleDimensions 30 30 30
sample SetImplicitFunction quadric
sample Update
#sample Print
sample ComputeNormalsOff
# Extract cells that contains isosurface of interest
vtkConnectivityFilter conn
conn SetInputConnection [sample GetOutputPort]
conn ScalarConnectivityOn
conn SetScalarRange 0.6 0.6
conn SetExtractionModeToCellSeededRegions
conn AddSeed 105
# Create a surface
vtkContourFilter contours
contours SetInputConnection [conn GetOutputPort]
# contours SetInputConnection [sample GetOutputPort]
contours GenerateValues 5 0.0 1.2
vtkDataSetMapper contMapper
# contMapper SetInputConnection [contours GetOutputPort]
contMapper SetInputConnection [conn GetOutputPort]
contMapper SetScalarRange 0.0 1.2
vtkActor contActor
contActor SetMapper contMapper
# Create outline
vtkOutlineFilter outline
outline SetInputConnection [sample GetOutputPort]
vtkPolyDataMapper outlineMapper
outlineMapper SetInputConnection [outline GetOutputPort]
vtkActor outlineActor
outlineActor SetMapper outlineMapper
eval [outlineActor GetProperty] SetColor 0 0 0
# Graphics
# create a window to render into
vtkRenderer ren1
vtkRenderWindow renWin
renWin AddRenderer ren1
# create a renderer
# interactiver renderer catches mouse events (optional)
vtkRenderWindowInteractor iren
iren SetRenderWindow renWin
ren1 SetBackground 1 1 1
ren1 AddActor contActor
ren1 AddActor outlineActor
ren1 ResetCamera
[ren1 GetActiveCamera] Zoom 1.4
iren AddObserver UserEvent {wm deiconify .vtkInteract}
iren Initialize
wm withdraw .
|
