医学影像系统原理4CT重建教程方案.pptVIP

思考题 能否采用比体元数少的投影方程数重建CT图像?为什么? 如果可能,请给出你的实现方案与步骤. * * 谢谢! Thank you! 即傅里叶切片定理。 角度为? 的线投影的一维傅立叶变换G? (?)等于物体的二维傅立叶变换域中的一条直线，该直线过原点且角度为?。 傅里叶变换 从投影到图像（中心切片定理） * 二维傅立叶变换法 假设图像矩阵f(x,y)的二维傅立叶变换为: 若空间频率坐标以极坐标来表示,则有: 则二维傅立叶变换可写为: (1) X射线投影路径 * X射线投影路径 由于在X射线扫描中, X射线的投影P总是与X射线路径L有关,为此引进一个新的坐标系(极坐标t- )来描述X射线路径L的位置. 设X射线路径L到极坐标中心O的距离为t,与y轴夹角为 ,则X射线路径L用直线方程表示为: 当在某一 角度时,则沿X射线路径L的投影为: (2) * 投影数据 Radon变换 * 投影与线积分 沿t方向的一维傅立叶变换 从投影数据的一维傅立叶变换恢复二维数据的傅立叶变换 推而广之 不同方向投影数据的傅立叶变换就可以得到二维傅立叶变换的近似估计值 投影变换定理图解 吸收系数图象的恢复 * 二维傅立叶变换法 将(1)式的投影函数进行一维的傅立叶变换,是在 角度时投影的傅立叶变换形式: 由上面分析可见,图像在某一 角度上投影的傅立叶变换正好等于该图像吸收系数相对应角度( ) 剖面上的二维傅立叶变换.若在整个 平面上添满了全部角度投影的傅立叶变换值,那么就可以用图像的二维傅立叶变换来表示原图像f(x,y)的变换值,再由二维傅立叶反变换即可得到重建后的原图像. * 二维傅立叶变换法 * * 傅立叶变换反投影法 二维傅立叶变换法被认为是最理想的图像重建算法之一.但因该方法需要进行正反两次傅立叶变换,计算量比较大,在实际应用中不易实现. 在反投影法中,反投影重建图像的吸收系数与实际图像之间存在一个模糊因子,如果利用二维傅立叶变换表示,则: 上式说明校正模糊失真的步骤可先将fb(x,y)做二维傅立叶变换,然后将变换结果用 加权,即可得到真正图像的二维傅立叶变换.通过反变换就可重建真实图像. 滤波反投影法 滤波反投影也是解析法中的一种.这种方法为了消除模糊因子的影响,并将二维傅立叶变换改变为只进行一维傅立叶变换,既可校正失真,又可简化计算,提高重建速度. 采用卷积计算的滤波反投影法是当前CT成像装置中应用最为广泛,故也称之卷积反投影法(Convolution Back Projection, CBP). * 滤波反投影的基本原理 根据反投影法中关系式: 式中的投影 用一维傅立叶反变换代替,同时用 乘除变换上式得到: (3) * 滤波反投影的基本原理 (3)式与二维傅立叶变换法公式比较,会发现二维傅立叶变换重建原图像吸收系数f(x,y)与反投影法重建图像吸收系数fb(x,y)之间正好相差一个 .这里 正是 的傅立叶变换形式,即反投影重建图像fb(x,y)由于模糊因子 造成图像模糊,使fb(x,y)与f(x,y)有误差. 消除模糊因子的影响,可以采取对每一投影的傅立叶变换值用 加权,以产生失真的重建图像. * 卷积计算的实现 根据傅立叶变换的卷积定理: 上式说明在频域内,投影的傅立叶变换 用 进行变换或滤波,等效于投影 与滤波函数 的傅立叶反变换进行卷积计算. 的傅立叶反变换在时域中可以构造成滤波函数h(t).通过选取不同的滤波函数,对于投影 进行有效地滤波,达到满意的重建图像效果. 卷积计算中存在的主要问题是如何确定空间滤波函数h(t) ,它是CT成像软件中的核心机密. * * * 滤波反投影法（Filter Back projection） 示例 投影 反投影 反投影重建 滤波投影数据 滤波 滤波反投影重建 滤波反投影 * * 医学影像系统原理:CT图像重建算法 丁明跃 华中科技大学生物医学工程系 “图像信息处理与智能控制”教育部重点实验室 * * 一、反投影法 二、傅立叶变换法 三、滤波反投影法 目录 * X射线衰减的计算 根据朗勃-比尔(Lamber-Beer)吸收定理,X射线穿过长度为L,吸收系数为 的体元,则穿过该体元后的X射线强度为: 穿过第二个体元时的X射线强度为: 类似地,穿过第n个体元时的X射线强度为: (1) * X射线投影与CT扫描(CT scan) 在CT重建过程中,为了确定不同点的吸收系数,我们需要首先计算出每个体元的吸收系数.