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
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
|
package require vtk
package require vtkinteraction
# Create the RenderWindow, Renderer and Interactor
#
vtkRenderer ren1
vtkRenderWindow renWin
renWin AddRenderer ren1
vtkRenderWindowInteractor iren
iren SetRenderWindow renWin
vtkTIFFReader imageIn
imageIn SetFileName "$VTK_DATA_ROOT/Data/beach.tif"
# "beach.tif" image contains ORIENTATION tag which is
# ORIENTATION_TOPLEFT (row 0 top, col 0 lhs) type. The TIFF
# reader parses this tag and sets the internal TIFF image
# orientation accordingly. To overwrite this orientation with a vtk
# convention of ORIENTATION_BOTLEFT (row 0 bottom, col 0 lhs ), invoke
# SetOrientationType method with parameter value of 4.
imageIn SetOrientationType 4
imageIn ReleaseDataFlagOff
imageIn Update
proc PowerOfTwo {amt} {
set pow 0
incr amt -1
while {1} {
set amt [expr $amt >> 1]
incr pow
if {$amt <= 0} {return [expr 1 << $pow];}
}
}
set orgX [expr [lindex [[imageIn GetOutput] GetWholeExtent] 1] - [lindex [[imageIn GetOutput] GetWholeExtent] 0] + 1]
set orgY [expr [lindex [[imageIn GetOutput] GetWholeExtent] 3] - [lindex [[imageIn GetOutput] GetWholeExtent] 2] + 1]
set padX [PowerOfTwo $orgX]
set padY [PowerOfTwo $orgY]
vtkImageConstantPad imagePowerOf2
imagePowerOf2 SetInputConnection [imageIn GetOutputPort]
imagePowerOf2 SetOutputWholeExtent 0 [expr $padX - 1] 0 [expr $padY - 1] 0 0
vtkImageRGBToHSV toHSV
toHSV SetInputConnection [imageIn GetOutputPort]
toHSV ReleaseDataFlagOff
vtkImageExtractComponents extractImage
extractImage SetInputConnection [toHSV GetOutputPort]
extractImage SetComponents 2
extractImage ReleaseDataFlagOff
vtkImageThreshold threshold1
threshold1 SetInputConnection [extractImage GetOutputPort]
threshold1 ThresholdByUpper 230
threshold1 SetInValue 255
threshold1 SetOutValue 0
threshold1 Update
set extent [[threshold1 GetOutput] GetWholeExtent]
vtkImageSeedConnectivity connect
connect SetInputConnection [threshold1 GetOutputPort]
connect SetInputConnectValue 255
connect SetOutputConnectedValue 255
connect SetOutputUnconnectedValue 0
connect AddSeed [lindex $extent 0] [lindex $extent 2]
connect AddSeed [lindex $extent 1] [lindex $extent 2]
connect AddSeed [lindex $extent 1] [lindex $extent 3]
connect AddSeed [lindex $extent 0] [lindex $extent 3]
vtkImageGaussianSmooth smooth
smooth SetDimensionality 2
smooth SetStandardDeviation 1 1
smooth SetInputConnection [connect GetOutputPort]
vtkImageShrink3D shrink
shrink SetInputConnection [smooth GetOutputPort]
shrink SetShrinkFactors 2 2 1
shrink AveragingOn
vtkImageDataGeometryFilter geometry
geometry SetInputConnection [shrink GetOutputPort]
vtkTextureMapToPlane geometryTexture
geometryTexture SetInputConnection [geometry GetOutputPort]
geometryTexture SetOrigin 0 0 0
geometryTexture SetPoint1 [expr $padX - 1] 0 0
geometryTexture SetPoint2 0 [expr $padY - 1] 0
vtkCastToConcrete geometryPD
geometryPD SetInputConnection [geometryTexture GetOutputPort]
vtkClipPolyData clip
clip SetInput [geometryPD GetPolyDataOutput]
clip SetValue 5.5
clip GenerateClipScalarsOff
clip InsideOutOff
clip InsideOutOn
[[clip GetOutput] GetPointData] CopyScalarsOff
clip Update
vtkTriangleFilter triangles
triangles SetInputConnection [clip GetOutputPort]
vtkDecimatePro decimate
decimate SetInputConnection [triangles GetOutputPort]
decimate BoundaryVertexDeletionOn
decimate SetDegree 25
decimate PreserveTopologyOn
vtkLinearExtrusionFilter extrude
extrude SetInputConnection [decimate GetOutputPort]
extrude SetExtrusionType 2
extrude SetScaleFactor -20
vtkPolyDataNormals normals
normals SetInputConnection [extrude GetOutputPort]
normals SetFeatureAngle 80
vtkStripper strip
strip SetInputConnection [extrude GetOutputPort]
vtkPolyDataMapper map
map SetInputConnection [strip GetOutputPort]
map SetInputConnection [normals GetOutputPort]
map ScalarVisibilityOff
vtkTexture imageTexture
imageTexture InterpolateOn
imageTexture SetInputConnection [imagePowerOf2 GetOutputPort]
vtkActor clipart
clipart SetMapper map
clipart SetTexture imageTexture
ren1 AddActor clipart
[clipart GetProperty] SetDiffuseColor 1 1 1
[clipart GetProperty] SetSpecular .5
[clipart GetProperty] SetSpecularPower 30
[clipart GetProperty] SetDiffuse .9
ren1 ResetCamera
set camera [ren1 GetActiveCamera]
$camera Azimuth 30
$camera Elevation -30
$camera Dolly 1.5
ren1 ResetCameraClippingRange
ren1 SetBackground 0.2 0.3 0.4
renWin SetSize 320 256
iren Initialize
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 .
|
