CT原理(下).ppt
《CT原理(下).ppt》由会员分享,可在线阅读,更多相关《CT原理(下).ppt(67页珍藏版)》请在淘文阁 - 分享文档赚钱的网站上搜索。
1、CTCT 定义定义 奥地利数学家雷登奥地利数学家雷登(J Radon)(J Radon)提出了提出了图像重建图像重建理论。理论。 这一理论指出由物体的一组横断面的投影来重建物体的内这一理论指出由物体的一组横断面的投影来重建物体的内部结构的图像。部结构的图像。 已被广泛应用于放射学、非破坏性工业测试和数据压缩等已被广泛应用于放射学、非破坏性工业测试和数据压缩等许多领域。许多领域。 CT CT是图像重建在医学上获得的最重要的应用之一。是图像重建在医学上获得的最重要的应用之一。CTCT图像重建方法图像重建方法 用此数学方法重建的图像应能正确反映被扫描用此数学方法重建的图像应能正确反映被扫描部位的物质
2、特性和确切位置,同时能被计算机所接部位的物质特性和确切位置,同时能被计算机所接受,受, 用计算机去快速处理,实现这种算法。用计算机去快速处理,实现这种算法。 图像重建方法的分类图像重建方法的分类: : 直接方法直接方法: (1) (1) 矩阵法矩阵法 (2)(2)迭代法迭代法 直接计算线形方程系数直接计算线形方程系数的方法。的方法。 间接方法间接方法: (1)(1)反投影法反投影法 (2)(2)卷积法卷积法 (3) (3) 傅利叶重建法傅利叶重建法 先计算傅利叶变换系数,再导出吸收系数先计算傅利叶变换系数,再导出吸收系数的方法。的方法。 CTCT有关的图象重建方法有有关的图象重建方法有: :(
3、1)(1)直接反投影法直接反投影法(2)(2)迭代法迭代法 A. A. 同时迭代重建法同时迭代重建法 B. B. 代数重建法代数重建法 C. C. 迭代最小二乘法迭代最小二乘法(3)(3)解析法解析法 A. A. 二维傅利叶重建法二维傅利叶重建法 B B 空间滤波反投影法空间滤波反投影法 C. C. 卷积法反投影法卷积法反投影法图像重建方法:图像重建方法:矩阵法矩阵法 又叫解联立方程法又叫解联立方程法 A射线方向:射线方向:X1 + X2 = 2 式式B射线方向:射线方向:X3 + X4 = 4 式式C射线方向:射线方向:X1 + X3 = 1 式式D射线方向:射线方向:X2 + X4 = 5
4、 式式式式 + 式式 式式 + 式式所以有一个是派生出的方程。所以有一个是派生出的方程。图像重建方法:图像重建方法:矩阵法矩阵法 又叫解联立方程法又叫解联立方程法 E射线方向:射线方向:X1 + X4 = 3 F射线方向:射线方向:X2 + X3 = 3 图像重建方法:图像重建方法:矩阵法矩阵法 又叫解联立方程法又叫解联立方程法 解联立方程,解联立方程,得:得:X1 = 0X2 = 3X3 = 1X4 = 3主要缺点:主要缺点:矩阵大时,计算矩阵大时,计算方法时间太长。方法时间太长。 图像重建方法:图像重建方法: 反投影法,又叫反投影法,又叫“投影反投影反馈法馈法”,或,或“总和法总和法”,是
5、是CTCT重建的基本方法。重建的基本方法。 反投影法是一种基本的反投影法是一种基本的重建方法,基本思想是重建方法,基本思想是将每次将每次“投影投影”重新放重新放回矩阵之中,并继续取回矩阵之中,并继续取它们的射线和。它们的射线和。 对矩阵作对矩阵作0 0o o、4545o o、9090o o、135135o o投影投影( (即扫描即扫描) ),再将,再将投影值反投回原投影值反投回原矩阵的对应位置矩阵的对应位置( (即扫描过的各即扫描过的各个体素个体素) )上,即上,即可解出原矩阵中可解出原矩阵中元素的元素的值。值。运算中的基数等于所有体素的特征参数的总和,运算中的基数等于所有体素的特征参数的总和
6、,这个总和也等于任一方向上的投影值的总和。这个总和也等于任一方向上的投影值的总和。此算法由计算机执行此算法由计算机执行, ,实现起来简单而不需要实现起来简单而不需要很复杂的数学,因此计算速度比较快。很复杂的数学,因此计算速度比较快。 反投影法的缺点:反投影法的缺点: 会出现图像会出现图像边缘失锐现边缘失锐现象象( (也是一也是一 种伪影种伪影) )。 图像重建方法:图像重建方法:卷积反投影法卷积反投影法 为了消除反投影法产生的图像的边缘失锐,为了消除反投影法产生的图像的边缘失锐, 在实际中采用的算法是在实际中采用的算法是卷卷积积反投影法反投影法 ( filtered back projecti
7、on)。 此方法是先把获得的投影函数作卷积处理,此方法是先把获得的投影函数作卷积处理, 所谓卷积即人为设计一种滤波函数,用它所谓卷积即人为设计一种滤波函数,用它 对得到的投影函数进行改造。对得到的投影函数进行改造。 然后把这些改造过的投影函数进行相加等然后把这些改造过的投影函数进行相加等 处理,就可以达到消除星状伪影的目的处理,就可以达到消除星状伪影的目的 。CTCT机现在普遍采用机现在普遍采用的数学算法为的数学算法为卷积卷积反投影法反投影法。右图说明有无卷积右图说明有无卷积对图像重建的影响。对图像重建的影响。 不同的滤波函数在不同的滤波函数在增加空间分辨率或增加空间分辨率或边缘增强的同时,边
8、缘增强的同时,也增加了图像噪声。也增加了图像噪声。 图像重建方法:图像重建方法:傅立叶变换法傅立叶变换法对数也对数也是一种是一种变换变换 是数字信号处理领域一种很重要的算法。是数字信号处理领域一种很重要的算法。傅立叶原理表明:任何连续测量的时序或信号,都可以傅立叶原理表明:任何连续测量的时序或信号,都可以表示为不同频率的正弦波信号的无限叠加。表示为不同频率的正弦波信号的无限叠加。傅立叶变换将原来难以处理的时域信号转换成了易于分傅立叶变换将原来难以处理的时域信号转换成了易于分析的频域信号(信号的频谱),可以利用一些工具对这析的频域信号(信号的频谱),可以利用一些工具对这些频域信号进行处理、加工。
9、些频域信号进行处理、加工。 图像重建方法:图像重建方法:傅立叶变换法傅立叶变换法 图像的图像的频率频率是表征图像中灰度变化剧烈程度的指标,是表征图像中灰度变化剧烈程度的指标,是灰度在平面空间上的梯度。傅立叶变换是将图像是灰度在平面空间上的梯度。傅立叶变换是将图像从空间从空间域转换到频率域域转换到频率域,换句话说,傅立叶变换的物理意义是将,换句话说,傅立叶变换的物理意义是将图像的灰度分布函数变换为图像的频率分布函数。图像的灰度分布函数变换为图像的频率分布函数。 由于空间是三维的,图像是二维的,因此空间中物体由于空间是三维的,图像是二维的,因此空间中物体在另一个维度上的关系就由梯度来表示,这样我们
10、可以通在另一个维度上的关系就由梯度来表示,这样我们可以通过观察图像得知物体在三维空间中的对应关系。过观察图像得知物体在三维空间中的对应关系。 图像重建方法:图像重建方法:傅立叶变换法傅立叶变换法 是在是在解矩阵方程解矩阵方程时常用的方法。时常用的方法。应用迭代法时,开始时可任意设出矩阵中的应用迭代法时,开始时可任意设出矩阵中的值值( (一般一般都假设图像是均匀的都假设图像是均匀的) ),然后将计算值与投影实测值,然后将计算值与投影实测值比较,并以计算值与实测值之间的差加以修正,然后比较,并以计算值与实测值之间的差加以修正,然后一遍遍地重复,直到假设值与测量值一样或在允许的一遍遍地重复,直到假设
11、值与测量值一样或在允许的误差范围内为止误差范围内为止图像重建方法:图像重建方法:迭代法迭代法, ,又叫又叫“逐步近似法逐步近似法” 第一步:给出初始值第一步:给出初始值 第二步:在垂直行中的新的数第二步:在垂直行中的新的数相加,产生新射线和,并与初相加,产生新射线和,并与初始测量的垂直射线和比较始测量的垂直射线和比较, 将将结果均分后放入新单元中。结果均分后放入新单元中。 第三步:在新建单元中取水平第三步:在新建单元中取水平方向投影,并把它与待重建单方向投影,并把它与待重建单元水平方向投影相比较,并将元水平方向投影相比较,并将结果均分后放入新单元中。结果均分后放入新单元中。 第四步:在第四步:
12、在45o及及135o方向投方向投影,对角线上的射线和也重复影,对角线上的射线和也重复上述过程并完成第一次迭代。上述过程并完成第一次迭代。 1.速度速度在速度方面,迭代法是最慢的。卷积反投影法和傅里叶变在速度方面,迭代法是最慢的。卷积反投影法和傅里叶变换重建法对于每投影在记录后便能立即处理,因此全部重换重建法对于每投影在记录后便能立即处理,因此全部重建在最后一次投影完成后就能在极短时间内显示出来。建在最后一次投影完成后就能在极短时间内显示出来。2.准确性准确性卷积反投影法和傅里叶变换重建法对具有良好的圆对称的卷积反投影法和傅里叶变换重建法对具有良好的圆对称的物体很适用,而迭代法更有利于处理对称性
13、较差的物体。物体很适用,而迭代法更有利于处理对称性较差的物体。 数据足够的条件下,可得精确得结果。数据足够的条件下,可得精确得结果。 CTCT CT图像从数据采集到图像从数据采集到图像重建要经过图像重建要经过X射线射线产生产生,信号采集信号采集, 预处预处理理, 卷积和反投影卷积和反投影 多个环节的运算处理,多个环节的运算处理,以确保得到的数据和以确保得到的数据和图像能正确反映被扫图像能正确反映被扫描物体的组织结构状描物体的组织结构状况。况。 今天今天CT能够成为一种可行的医疗设备的真正秘能够成为一种可行的医疗设备的真正秘密,既不是前面讨论的重建算法,也不是后面将要密,既不是前面讨论的重建算法
14、,也不是后面将要介绍的图像显示方法,而是数据的处理和校正方法。介绍的图像显示方法,而是数据的处理和校正方法。 事实上,这是所有事实上,这是所有CT制造商很少公开讨论的技制造商很少公开讨论的技术。这些校正步骤包括预处理和后处理的组成部分。术。这些校正步骤包括预处理和后处理的组成部分。CTCT成像分为以下三个基本过程:成像分为以下三个基本过程:数据扫描数据扫描图像重建图像重建图像显示图像显示CTCT用特殊的算法来进行数学计算,必须保证采集到用特殊的算法来进行数学计算,必须保证采集到的数据不受机器或其它外界因来的干扰,的数据不受机器或其它外界因来的干扰, 所接收到所接收到的衰减数据必须真实。的衰减数
15、据必须真实。 在进行图像重建之前,如要得到准确的重在进行图像重建之前,如要得到准确的重 建图像数据,必须对这些数据进行处理,也叫建图像数据,必须对这些数据进行处理,也叫 CTCT数据预处理。数据预处理。在图像处理机中,数据预处理是产生优质图像在图像处理机中,数据预处理是产生优质图像的必不可少的一部分工作,而相应表格的修正工作的必不可少的一部分工作,而相应表格的修正工作则定期由一定经验的工程师来负责进行。则定期由一定经验的工程师来负责进行。 预处理的内容预处理的内容n 补偿误差补偿误差 n KV脉动脉动 n DAS灵敏度漂移灵敏度漂移 n 射线硬化射线硬化 n 空间几何形状偏差空间几何形状偏差
16、n 由于扫描机架几何形状引起各通道灵敏度的差异由于扫描机架几何形状引起各通道灵敏度的差异 n 由于扇形束硬化效应及通道的非直线性而引起灵敏度差异由于扇形束硬化效应及通道的非直线性而引起灵敏度差异n 水的定标水的定标 这些误差的产生是不可避免,然而我们在建立图像前可这些误差的产生是不可避免,然而我们在建立图像前可以设法校正这些不足之处。在开始扫描时,记下机器以设法校正这些不足之处。在开始扫描时,记下机器的实际值,这就是说测出机器的重要特征。的实际值,这就是说测出机器的重要特征。 如放大器的补偿灵敏度等,然后制表储存于磁盘中。在如放大器的补偿灵敏度等,然后制表储存于磁盘中。在图像重建前的扫描处理过
17、程中使用这些表中的数值进图像重建前的扫描处理过程中使用这些表中的数值进行修正行修正, ,从而得从而得“理想化理想化”的数据。这些表格很有用的数据。这些表格很有用, ,被称为被称为“数据处理数据处理”表。表。 事实上,没有任何两个元件的参数是一致的,也没有事实上,没有任何两个元件的参数是一致的,也没有任何两个放大器的灵敏度是一致的。任何两个放大器的灵敏度是一致的。 分类分类在分类预处理步骤中在分类预处理步骤中,若探测器有若探测器有512个单元通个单元通道道,则按则按1512的顺序分类。并对所有数字信的顺序分类。并对所有数字信号进行校对和译码。号进行校对和译码。 预处理要做以下几方面工作预处理要做
18、以下几方面工作 补偿校正补偿校正 因为作为电子通道的众多放大器部件,存在着因为作为电子通道的众多放大器部件,存在着不同的零点漂移现象不同的零点漂移现象,即在没有输入信号的情况下即在没有输入信号的情况下,产生不同的输出信号,而系统所需求的精确性不允产生不同的输出信号,而系统所需求的精确性不允许忽视这种影响。所以补偿校正是有必需的。许忽视这种影响。所以补偿校正是有必需的。 补偿校正完成后存放在补偿校正完成后存放在OFFSET表中,如果表表中,如果表中的数据不好,重建的图像将出现中的数据不好,重建的图像将出现 伪影。伪影。 对数运算对数运算 因因X X线束是按指数规律衰减的,故需对所测得的线束是按指
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- CT 原理
限制150内