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
|
package require vtk
vtkImageReader reader
reader ReleaseDataFlagOff
reader SetDataByteOrderToLittleEndian
reader SetDataExtent 0 63 0 63 1 93
reader SetFilePrefix "$VTK_DATA_ROOT/Data/headsq/quarter"
reader SetDataMask 0x7fff
set rangeStart 0.0
set rangeEnd 0.2
vtkLookupTable LUT
LUT SetTableRange 0 1800
LUT SetSaturationRange 1 1
LUT SetHueRange $rangeStart $rangeEnd
LUT SetValueRange 1 1
LUT SetAlphaRange 1 1
LUT Build
# added these unused default arguments so that the prototype
# matches as required in python.
proc changeLUT { { a 0 } { b 0 } } {
global rangeStart rangeEnd
set rangeStart [expr $rangeStart + 0.1]
set rangeEnd [expr $rangeEnd + 0.1]
if { $rangeEnd > 1.0 } {
set rangeStart 0.0
set rangeEnd 0.2
}
LUT SetHueRange $rangeStart $rangeEnd
LUT Build
}
vtkImageMapToColors mapToRGBA
mapToRGBA SetInputConnection [reader GetOutputPort]
mapToRGBA SetOutputFormatToRGBA
mapToRGBA SetLookupTable LUT
mapToRGBA AddObserver "EndEvent" changeLUT
vtkMemoryLimitImageDataStreamer streamer
streamer SetInputConnection [mapToRGBA GetOutputPort]
streamer SetMemoryLimit 100
streamer UpdateWholeExtent
# set the window/level to 255.0/127.5 to view full range
vtkImageViewer viewer
viewer SetInputConnection [streamer GetOutputPort]
viewer SetColorWindow 255.0
viewer SetColorLevel 127.5
viewer SetZSlice 50
wm withdraw .
viewer Render
|