CT图像重建资料讲解

C T图像重建

昆明理工大学信息工程与自动化学院学生实验报告

（ 2009—2010学年第一学期）

一、实验目的与意义

医学成像技术是生物医学工程专业的一门重要的专业课程,课程主要涉及X光仪器，CT 仪器，MRI仪器和核医学仪器的工作原理及成像方法。其中CT算法的出现又为后来数字化医学成像技术的发展提供了基础。该门课程为生物医学工程专业的专业基础课。

CT技术是医学成像系统中的一种重要手段。它通过特定的算法，利用计算机的高速运算功能，可以在短时间内快速呈现人体断层图像。让学生练习CT图像的重建有助于学生理解CT算法的内容，熟悉数字图像重建的过程。同时也能培养学生的团队精神和解决实际问题的能力。

二、实验算法原理

1、MATLAB处理数字图像的基本函数；

2、X-CT三维图像重建的基本算法。

CT图象重建有四种基本的算法:矩阵法,迭代法,傅立叶算法,反投影算法.我们采用的方法为卷积反投影. 卷积反投影有:平行光束投影的卷积反投影算法, 等角扇形光来投影的重建算法.

1）.平行光束投影的卷积反投影算法

仅供学习与交流，如有侵权请联系网站删除谢谢2

仅供学习与交流，如有侵权请联系网站删除 谢谢3

从投影重建三维物体的图像，就是重建一个个横断面。这样三堆图像的重建就归结为二维图象的重建。二维图像的重建问题可以从数学上描述如下。

假定),(y x g 表示一个二维的未知函数，通过),(y x g 的直线称为光钱(见图2．1)。沿光线),(y x g 的积分称作光线积分。沿相同方向的一组光线积分，就构成一个投影。图2．1中垂直于直线'CC (与X 轴夹角为θ)的光线所形成。

图2.1 ),(y x g 在θ方向的投影)(t P θ 的投影)(t P θ，称之为),(y x g 在θ方向的投影。光线积分和投影在数学上可以定义如下：

在图2.1中直线AB 的方程为：

1sin cos t Y X =+θθ (2.1) 其中1t 是AB 到原点的距离，),(y x g 沿AB 的积分为：

dxdy t y x y x g ds y x g t P AB )sin cos (),(),()(11-+==??+∞

∞-θθδθ (2.2) 对于给定的θ，),(y x g 在θ方向的投影)(t P θ是t 的函数。如果),(y x g 在各个方向的投影已知，),(y x g 就可以唯一确定。下面就讨论卷积反投影重建算法。

仅供学习与交流，如有侵权请联系网站删除 谢谢4 假定投影方向θ，如图2.2，将坐标),(y x 旋转θ角（逆时针方向）形成坐标系),(s t 。),(y x g 在),(s t 坐标系中为),(s t g 。

图2.2 傅立叶切片定理示意图

坐标系),(s t 与),(y x 之间的关系为：

???

? ?????? ??-=???? ??y x s t θθθθcos sin sin cos

（2.3）

显然

()ds s t g t P ?+∞

∞-=),(θ （2.4）

令)(w S θ为)(t P θ的傅立叶变换则

dt wt j t P w S )2ex p()()(πθθ-=?+∞

∞- dsdt wt j s t g )2ex p(),(π-=?+∞

∞- （2.5）

将上式变换到),(y x 坐标系中，注意到变换的可比行列式