CT图像重建

《CT图像重建》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、..-XX理工大学信息工程与自动化学院学生实验报告〔2009—2010学年第一学期〕课程名称：医学成像系统与放射治疗装置开课实验室：32082008年12月24日年级、专业、班学号成绩实验工程名称CT图像重建指导教师X利军教师评语教师签名：年月日一、实验目的与意义医学成像技术是生物医学工程专业的一门重要的专业课程,课程主要涉及X光仪器，CT仪器，MRI仪器和核医学仪器的工作原理及成像方法。其中CT算法的出现又为后来数字化医学成像技术的开展提供了根底。该门课程为生物医学工程专业的专业根底课。CT技术是医学成像系统中的一种重要手段。它通过特定的算法，利用计算机的高速运算

2、功能，可以在短时间内快速呈现人体断层图像。让学生练习CT图像的重建有助于学生理解CT算法的内容，熟悉数字图像重建的过程。同时也能培养学生的团队精神和解决实际问题的能力。二、实验算法原理1、MATLAB处理数字图像的根本函数；2、X-CT三维图像重建的根本算法。CT图象重建有四种根本的算法:矩阵法,迭代法,傅立叶算法,反投影算法.我们采用的方法为卷积反投影.卷积反投影有:平行光束投影的卷积反投影算法,等角扇形光来投影的重建算法.-.word.zl-..-1〕.平行光束投影的卷积反投影算法从投影重建三维物体的图像，就是重建一个个横断面。这样三堆图像的重建就归结为二维图象

3、的重建。二维图像的重建问题可以从数学上描述如下。假定表示一个二维的未知函数，通过的直线称为光钱(见图2．1)。沿光线的积分称作光线积分。沿一样方向的一组光线积分，就构成一个投影。图2．1中垂直于直线(与X轴夹角为)的光线所形成。图2.1在方向的投影的投影，称之为在方向的投影。光线积分和投影在数学上可以定义如下：在图2.1中直线AB的方程为：(2.1)其中是AB到原点的距离，沿AB的积分为：(2.2)对于给定的，在方向的投影是t的函数。如果在各个方向的投影，-.word.zl-..-就可以唯一确定。下面就讨论卷积反投影重建算法。假定投影方向，如图2.2，将坐标旋转角〔

4、逆时针方向〕形成坐标系。在坐标系中为。图2.2傅立叶切片定理示意图坐标系与之间的关系为：〔2.3〕显然〔2.4〕令为的傅立叶变换那么〔2.5〕将上式变换到坐标系中，注意到变换的可比行列式〔2.6〕-.word.zl-..-从而得到：〔2.7〕其中〔2.8〕假设令的傅立叶变换为，由〔2.8〕可知〔2.9〕假设的傅立叶变换为的极坐标表示。这说明在方向的投影傅立叶变换等于在与u轴成角的直线上的值。这就是著名的傅立叶投影切片定理。可见在整个平面可以利用各个方向的投影来得到，从而也可以通过求的傅立叶反交换的方法求得：〔2.10〕变换到极坐标中，得到〔2.11〕经推导得〔2.1

5、2〕其中假设令-.word.zl-..-〔2.13〕那么(2.14)(2.13)式右端是两频谱函数和的乘积的傅立叶反变换。是投影傅立叶变换。假设的傅立叶反变换为，那么根据卷积定理有：〔2.15〕或其中〔2.16〕当图像的频谱是有限带宽时，那么上式变为〔2.17〕由于图象及其频谱都是离散采样的，假定图象采样间隔为，那么根据采样定理。为了进展数学处理，只需知道h(t)在有限带宽上的离散采样点的值．这样我们有〔2.18〕其中n为正负整数。〔2.18〕的离散形式为〔2.19〕假定在之外的值为0，那么上式变为〔2.20〕-.word.zl-..-或〔2.21〕其中从而可见为确

6、定的N个采样点上的的值，需要使用的2N—1个点上的值，从n=一(N—1)到(N—1)。为求得，利用傅立叶变换计算卷积是比拟快的方法，为去除循环卷积的周期交叠效应，实际上取2N个点，补0，使之有(2N—1)个元素，那么在N个采样点上就防止了交叠，如果使用以2为基的FFT(快速傅立叶变换)算法，和都必需朴0至(2N一1)个元素，(2N一1)为大于等于2N—l的最小的2的整数幂。计算的过程可以写为〔2.22〕其中FFT和IFFT分别表示快速傅立叶变换和反变换，光滑窗是在滤波过程中参加的光滑因子，例如引用汉明窗，有时可以改良重建效果。对于各个方向的投影，得到之后就可以由(2

7、.22)来计算。重建步骤可以归纳为：第一步：卷积，也称滤波，由(2．22)对每个方向计算。第二步。反投影，由(2．14)的近似形式〔2.23〕来计算的近似值。M为投影个数为投影方向角，他们均匀的分布在0~的范围内。当计算时，，不一定在的整离散点上，这就需要插值求得，预先将插值加密，即最靠近的点，可以提高计算速度。2).等角扇形光来投影的重建算法几乎所有的快遗CT设备都是用的扇形光束。这里表达的是等角度光束投影，如图2.3，测量投影数据的探测器等间距地分布在弧上，弧的半径为2D，D为光源到图像中心的距离。在下文中，-.word.zl-..-图象在极坐标中的表示。表