目录
第一章 绪论
第二章 从投影重建图像算法I :反投影重建算法(累加法)
第三章 从投影重建图像算法lI :滤波反投影重建算法
第四章 扇束投影重建
第五章 从投影重建图像算法Ⅲ:迭代重建算法
第二章、从投影重建图像算法I :反投影重建算法(累加法)
1.反投影重建算法的物理概念
断层平面中某一点的密度值可看作这一平面内所有经过该点射线(反)投影之和(的平均值)。整幅重建图像可看作所有方向下的(反)投影累加而成。
例如,图2.2 为一与纸平面平行的断层,具有16个像素。像素值(代表密度)分别为
x
1
.
.
.
x
16
x_1...x_{16}
x1...x16。我们赋值如下:
x
6
=
5
,
x
7
=
2
,
x
10
=
1
x_6 = 5,x_7=2,x_{10}=1
x6=5,x7=2,x10=1 , 其余均为零。射线的标号见下图,根据投影的定义有:
p
1
=
x
3
+
x
7
+
x
11
+
x
15
=
2
,
p
2
=
x
2
+
x
0
+
x
10
+
x
14
=
5
+
1
=
6
,
p
3
=
x
5
+
x
8
+
x
7
+
x
8
=
5
+
2
=
7
,
p
4
=
x
8
+
x
10
+
x
11
+
x
12
=
1
,
p
5
=
x
1
+
x
3
+
x
11
+
x
18
=
5
,
p
6
=
x
4
+
x
7
+
x
10
+
x
13
=
1
+
2
=
3.
\begin{aligned} & p_1=x_3+x_7+x_{11}+x_{15}=2, \\ & p_2=x_2+x_0+x_{10}+x_{14}=5+1=6, \\ & p_3=x_5+x_8+x_7+x_8=5+2=7, \\ & p_4=x_8+x_{10}+x_{11}+x_{12}=1, \\ & p_5=x_1+x_3+x_{11}+x_{18}=5,\\ & p_6=x_4+x_7+x_{10}+x_{13}=1+2=3. \end{aligned}
p1=x3+x7+x11+x15=2,p2=x2+x0+x10+x14=5+1=6,p3=x5+x8+x7+x8=5+2=7,p4=x8+x10+x11+x12=1,p5=x1+x3+x11+x18=5,p6=x4+x7+x10+x13=1+2=3.
若将这些投影值用反投影法重建图像,得:
x
1
=
p
5
=
5
,
x
2
=
p
2
=
6
,
x
3
=
p
1
=
2
,
x
4
=
p
0
=
3
,
x
5
=
p
3
=
7
,
x
6
=
p
2
+
p
3
+
p
5
=
6
+
7
+
5
=
18
,
x
7
=
p
1
+
p
3
+
p
6
=
2
+
7
+
3
=
12
,
x
8
=
p
3
=
7
,
x
9
=
p
4
=
1
,
x
10
=
p
2
+
p
4
+
p
0
=
6
+
1
+
3
=
10
,
x
11
=
p
1
+
p
4
+
p
5
=
2
+
1
+
5
=
8
,
x
12
=
p
4
=
1
,
x
13
=
p
0
=
3
,
x
14
=
p
2
=
6
,
x
15
=
p
1
=
2
,
x
18
=
p
5
=
5
。
\begin{aligned} & x_1=p_5=5, \\ & x_2=p_2=6, \\ & x_3=p_1=2, \\ & x_4=p_0=3, \\ & x_5=p_3=7, \\ & x_6=p_2+p_3+p_5=6+7+5=18, \\ & x_7=p_1+p_3+p_6=2+7+3=12, \\ & x_8=p_3=7, \\ & x_9=p_4=1, \\ & x_{10}=p_2+p_4+p_0=6+1+3=10, \\ & x_{11}=p_1+p_4+p_5=2+1+5=8, \\ & x_{12}=p_4=1, \\ & x_{13}=p_0=3, \\ & x_{14}=p_2=6, \\ & x_{15}=p_1=2, \\ & x_{18}=p_5=5。 \end{aligned}
x1=p5=5,x2=p2=6,x3=p1=2,x4=p0=3,x5=p3=7,x6=p2+p3+p5=6+7+5=18,x7=p1+p3+p6=2+7+3=12,x8=p3=7,x9=p4=1,x10=p2+p4+p0=6+1+3=10,x11=p1+p4+p5=2+1+5=8,x12=p4=1,x13=p0=3,x14=p2=6,x15=p1=2,x18=p5=5。
重建后的图像如图2.3(b) 所示。可以看出原图中像素值不为零的点反投影重建后仍较突出。但原图中像素值为零的像素,经反投影重建后不再为零。即有伪迹。有时为了使重建后图像的像素值更接近于原图的像素值,在求反投影时,把数据(见图2.3(b) )除以投影午的数目n, (也就是可能经过像素的射线数。本例中n=4) 。
2.反投影重建算法的数学描述
3.正弦图(sinogram)