第4章计算机断层成像技术(2).ppt

PrinciplesofRadiationImaging4.2CT的基本概念一、基本概念1、像素：图像显示的小方块，CT图像的最小单元。3、CT值CTnumber4、ExponentialAttenuationofX-rayRay-SumofX-rayAttenuation5、ComputedTomography（CT）ReconstructionIdea反投影法AlgebraicReconstructionTechnique(ART)4.3中心切片定理1、图像（image）2、投影线（ray）3、投影变换（Shadow/RadonTransform）4、投影（projection）的性质5、投影与重建二、TheFourier(Central)SliceTheoremFromProjectionstoImage1、θ＝0°，在x轴投影2、θ＝90°，在y轴投影3、任意θ一般情况傅里叶变换重建图像三、物体、傅里叶、Radon空间4.4平行束反投影重建算法原始数据--反投影成像根据反投影算法x1=p5=5x6=p2+p3+p5=18 …平均化处理，除以投影线数目xi=xi/6从这个过程可以看出重建后的图像一方面确实在中心点给出了一个较高的密度值，但另一方面，在周围的象素上也出现了一些不希望有的值，这就是星状伪迹。产生星状伪迹的原因在于：反投影重建的本质是把取自有限物体空间的射线投影均匀地回抹(反投影)到射线所及的无限空间的各点之上，包括原先像素值为零的点这就是直接反投影方法的局限性。直接反投影法的图像校正二、滤波反投影法滤波反投影法FilteredBackprojectionExample:ProjectionExample:FilteringExample:FilteredBackprojection三、卷积反投影法卷积运算是信号处理中常用的一种运算。卷积反投影方法是指将测得的投影数据与一个“核函数”作卷积运算，然后再用所得的结果作反投影。卷积反投影法我们考虑孤立点源反投影重建，中心点A经n条投影线投影后，投影值均为1： p1=p2=...=pn=1因此重建后而其他点均为1/n这类伪迹成为星状伪迹选择“核函数”为序列[-1/3，1，-1/3]，用这个序列与测量到的投影数据[0，1，0]做卷积运算后可得到一个新的序列[-1/3，1，-1/3]，用这个序列数据作反投影可得如下图(a)～(d)所示的结果。把各个不同方向的投影结果相加，可得如图(e)的结果。与上图(f)相比可以看出，卷积反投影的结果消除了直接反投影图像中的伪象。Example:ProjectionExample:convolutionExample:ConvolutedBackprojection投影函数与反投影重建方法定义1、R-L滤波函数2、S-L滤波函数4.5扇形束反投影重建算法一、等角度扇形束扫描的图像重建公式推导1、投影修正二、等距扇形束扫描的图像重建公式推导1、投影修正三、数据重排算法4.6迭代法（代数重建算法）代数重建算法一种最简单的迭代算法是所谓的代数重建技术，右图给出了一幅由四个象素组成的图像，若四个象素的值分别为5、7、6、2。则可以分别获得6个投影数据，包括两个水平方向，两个垂直方向和两个对角线方向，分别是7、11、9、13、12和8。这就是所能得到的所有的已知数据。迭代开始，先令所有的重建单元的值为0，第一步计算出垂直方向的投影值，分别都是0，如右上图。把这个计算值与实测值11和9相比较后，将其差值除以2以后分别加到相应的单元上去，就可得到垂直方向的迭代结果，如右下图。在此基础上可以再做水平方向的迭代，此时有计算值10、10，实测值为12、8，将它们比较后求出差值，分别加到有关的象素上去，结果是：最后再作对角线方向的迭代，就得到了所要求的真实数据。实际上，要求重建的矩阵很大，因此选代算法是非常耗费时间的。2、卷积运算3、加权反投影f4f3f2f1X射线f1=5f2=7f3=6f4=21197812135.54.55.54.51010000000原像素值反投影后像素值f(x,y)f(x,y)P(??t)P’(??t)1)Filttheprojectionswithafilter