医学成像技术课件（下）.ppt

课程要求 了解X线CT的投影重建的基本原理 重点掌握CT图像重建的中心切片定理 重点掌握反投影重建算法 重点掌握滤波反投影重建算法 CT重建图像 CT重建图像 CT重建图像 几种重建算法： 1. 线性代数方法 2.傅立叶重建 3.反向投影重建 4.滤波反向投影重建 CT重建图像 CT重建图像 CT重建图像 CT重建图像 傅立叶重建图像的原理 傅立叶重建图像的原理 傅立叶重建图像的原理 傅立叶重建图像的原理-中心切片定理 傅立叶重建图像的原理-中心切片定理 傅立叶重建图像的原理 傅立叶重建图像的原理 傅立叶重建图像的原理 傅立叶重建图像的原理 傅立叶重建图像的原理 傅立叶重建图像的原理 1. 中心切片定理指出了从投影重建图像的可能性。 2.对某物质断层采集不同角度下的投影函数，理论上对采集0-180度范围内连续取值的无穷多个投影。 3. 一个物体在任何角度下投影的傅立叶变换等于同一物体二维傅立叶变换在同一方向得到的切片直线（过0的切片） 4.将上述无穷多个切片直线组成图像的二维傅立叶变换 5. 对上述二维傅立叶变换求逆变换得到待重建图像 傅立叶重建图像的原理 的投影函数 可以看成是以R和 为参数的二维空间的函数。 傅立叶重建图像的原理 傅立叶重建图像的原理 Radon变换例子matlab程序 I = zeros(100,100); I(25:75, 25:75) = 1; imshow(I) [R,xp] = radon(I,[0 15 30 45 60 75 90 105 120 135 150 165 180]); figure(1); plot(xp,R(:,1)); title(R_{0^o} (x\prime)) figure(2); plot(xp,R(:,4)); title(R_{45^o} (x\prime)) 傅立叶重建图像的原理 傅立叶重建图像的原理 傅立叶重建图像的原理 傅立叶变换重建图像的流程图 傅立叶重建图像的原理 反投影法重建 反投影法重建 反投影法重建 反投影法重建 反投影法的基本思想： 将每次“投影”重新放回矩阵中，并继续取它们的射线和。 反投影法重建 反投影法重建 反投影法重建 反投影法重建 反投影法重建 最后计算是：从每个图像单元减去背景值（等于某投射角情况下个投影之和），本例中背景强度等于6，再将各吸收系数除以最大公约数 反投影法重建 反投影法重建 分析下面四体素矩阵反投影的重建过程 反投影法重建 反投影法重建 反投影法重建 反投影法重建 反投影法重建 反投影重建matlab程序： P = phantom(128); Angle=linspace(0,179,180) R = radon(P,Angle); I1 = iradon(R,Angle); I2 = iradon(R,Angle,linear,none); figure(1), imshow(P), title(Original) figure(2), imshow(I2,[]), title(Unfiltered backprojection) 反投影法重建 反投影法重建 反投影法重建 反投影法重建 1. 直接反投影法也称累加法，是最简单，最老式的方法。 2. 存在星形伪影或云晕状伪影，由于有这种失真存在被扫描的物体边缘不清晰。 3. 由于反投影本身固有的缺点，得出的图像模糊，目前已无实用价值。 滤波反投影 直接反投影所产生的图像边缘的云晕状伪影在数学上称为对原图像的一次褶积，要去除伪影就需要再做褶积解除，这一数学修正，它也叫滤波，就是在每一个投影波形上加上一个修正用的函数波形，这有二种方法，一种是褶积处理在空城中滤波，另一种是傅立叶转换在频域中滤波，经滤波处理后，每个投影波形不仅包含了代表X线强度的正向脉冲，同时其相邻二边又加上了反向的修正脉冲。将这些滤波函数与投影波形相加，云晕状阴影就被抵消掉了。抵消得越彻底，反投影后的重建图像就越接近原来物体。 反投影重建存在的问题 反投影重建存在的问题 反投影重建存在的问题 反投影重建存在的问题 反投影重建存在的问题 反投影重建存在的问题 反投影重建存在的问题 反投影重建存在的问题 直接反投影修正法 直接反投影修正法 滤波反投影重建 滤波反投影重建 滤波反投影重建 滤波反投影重建 滤波反投影重建 滤波反投影重建 （1）对某一角度下的投影函数作一维傅立叶变换 （2）对(1)的变换结果乘上一维权重因子 （3）对（2）的加权结果作一维逆傅立叶变换 （4）用（3）中得出的修正过的投影函数作直接反投影 （5）改变投影角度，重复（1）-（4）的过程，直到完成全部180度下的反投影 滤波反投影重建 滤波反投影重建 滤波反投影重建