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
|
package require vtk
# This script shows the result of an ideal highpass filter in spatial domain
# Image pipeline
vtkImageReader2Factory createReader
set reader [createReader CreateImageReader2 "$VTK_DATA_ROOT/Data/fullhead15.png"]
$reader SetFileName "$VTK_DATA_ROOT/Data/fullhead15.png"
vtkImageFFT fft
fft SetInputConnection [$reader GetOutputPort]
vtkImageIdealHighPass highPass
highPass SetInputConnection [fft GetOutputPort]
highPass SetXCutOff 0.1
highPass SetYCutOff 0.1
highPass ReleaseDataFlagOff
vtkImageRFFT rfft
rfft SetInputConnection [highPass GetOutputPort]
vtkImageExtractComponents real
real SetInputConnection [rfft GetOutputPort]
real SetComponents 0
vtkImageViewer viewer
viewer SetInputConnection [real GetOutputPort]
viewer SetColorWindow 500
viewer SetColorLevel 0
viewer Render
$reader UnRegister viewer
|
