《CT原理与算法》庄天戈 -学习笔记 之 二

目录

第一章 绪论
第二章 从投影重建图像算法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 x65,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)

  

4.反投影重建算法的实现


总结

  • 0
    点赞
  • 1
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值