Medical Imaging
Medical Imaging
Medical Imaging
3D
Radon Transform :
~
P (,z) =
G(x)e
ikx
Back Projection :
P (,z)d
0 ~
Map :
The PSF associated with the simple Bach projection is: PSF Or (x, y ) = O(x, y )
BF
1 r
1 x2 + y2
Or (x, y ) 123 4 4
O(x, y ) 123
1 x2 + y2 1 24 4 3 1 k
Or (k x ,k y ) so
O(k x ,k y )
O(k x ,k y ) = k Or (k x ,k y )
1. 2. P ( , z) =
Equate the z axis with a tilted reference frame x ' || z, y ' z and x ' = x cos( ) + y sin( ) x = x ' cos( ) y ' sin( ) y = x ' sin( ) + y ' cos( )
3.
Substitute # 2 into #1 and change integral to dx ' dy ' (still over all space) P ( , z) =
'
4.
Integrate along x ' and note that z is only a point along the x ' axis. P , x ' =
5.
ix ' k
x'
dx ' dy '
) O(x, y )e
x'
dxdy
Define the tilted k space coordinate system. k x = k x ' cos( ) ky ' sin( ) k y = k x ' sin( ) k y ' cos( )
8.
) O(x, y )e
dxdy
k ' =0
y
O(x, y
)e ik xeik y dxdy
x y y
k ' =0
y
Consider a 2-dimensional example of an emission imaging system. O(x,y) is the object function, describing the source distribution. The projection data, is the line integral along the projection direction.
P 0o, y =
( ) O(x, y )dx
The Central Slice Theorem can be seen as a consequence of the separability of a 2-D Fourier Transform.
o(kx ,k y )=
~
O(x, y)e
ikx x iky y
dxdy
iky y
dy
iky y
dxdy e dxdy
iky y i 0x
= o(0, ky )
Radon transform
P , y '
( )
2D FT
1D FT
o(kx ,k y )
p ,k y '
( )