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
|
from InsightToolkit import *
from sys import argv
scale = itkScaleTransform2_New()
reader = itkImageFileReaderUS2_New()
writer = itkImageFileWriterUS2_New()
reader.SetFileName( argv[1] )
writer.SetFileName( argv[2] )
parameters = scale.GetParameters()
parameters.SetElement( 0, eval( argv[3] ) )
parameters.SetElement( 1, eval( argv[3] ) )
reader.Update()
inputImage = reader.GetOutput()
size = inputImage.GetLargestPossibleRegion().GetSize()
centralPixel = itkIndex3()
centralPixel.SetElement( 0, size.GetElement(0) / 2 )
centralPixel.SetElement( 1, size.GetElement(1) / 2 )
centralPoint = itkPointD2()
spacing = inputImage.GetSpacing()
interpolator = itkLinearInterpolateImageFunctionUS2D_New()
centralPoint.SetElement(0, centralPixel.GetElement(0) )
centralPoint.SetElement(1, centralPixel.GetElement(1) )
scale.SetCenter( centralPoint )
scale.SetParameters( parameters )
resampler = itkResampleImageFilterUS2US2_New()
resampler.SetInput( reader.GetOutput() )
resampler.SetTransform( scale.GetPointer() )
resampler.SetInterpolator( interpolator.GetPointer() )
resampler.SetSize( size )
resampler.SetOutputSpacing( spacing )
resampler.Update()
writer.SetInput( resampler.GetOutput() )
writer.Update()
|
