9CT图像重建原理PPT.pptVIP

第三节 CT图像重建 一、从投影重建图像的原理 ——中心切片定理 二、从投影重建图像的算法 一、图像构成概念 （1）体层、像素、体素 （2）图像矩阵 （3）投影 投照受检体后出射的X线束强度I称为投影，投影的数值称为投影值，投影值的分布称为投影函数。 二、图像重建的基本要求 CT的图像重建是根据采集的数据，由计算机求解出图像矩阵中各个像素所对应的吸收系数，在求解计算过程中应满足图像重建的基本要求。 应不失真地反映被测人体层面上的图像信息； 要在尽可能短的时间内完成计算； 要在理论上和技术上可行。 0 0 θ R x y f(x,y) gθ(R) 投影线L 三、从投影重建图像的原理 ——中心切片定理 0 θ v u f(u,v) 1DFT Gθ(ρ) 物体空间 f(x,y) Radon空间 gθ(R) 傅立叶空间 F(ρ,θ) R R-1 F1 F2 F2-1 中心切片定理指出：密度函数f(x,y)在某一方向上的投影函数gθ(R)的一维傅立叶变换函数Gθ(ρ)是原密度函数f(x,y)的二维傅立叶变换函数F(ρ, θ)在(ρ, θ)平面上沿同一方向且过原点的直线上的值。 2D IFT 空间域 频域 1D FT 插值 四、图像重建算法 ——平行束反投影重建算法 （1）直接反投影法 （2）滤波反投影法 （3）卷积反投影法 （1）直接反投影法 基本原理是将所测得的投影值按其原路径反投影到投影线的各个像素上（即指定投影线上所有点的值等于所测得的投影值），再将所有反投影数据进行叠加，从而推断出原图像。 缺点是：图像四周会出现云晕状阴影，造成影像边缘的不清晰。 注：运算中的基数等于所有体素的特征参数的总和。 重建的物体图像不是正方形，变成了“星”状物，中心处吸收系数μ值越大，离中心越远，这就是图像的边缘失锐。 消除云晕状阴影的数理基础 直接反投影得到的密度函数： 反投影密度函数与实际密度函数之间的关系： 时域中两个函数的卷积的傅立叶变换是这两个函数的傅立叶变换的乘积。 2D IFT 空间域 频域 2D FT ×ρ 直接反投影 （2）滤波反投影法 原理：对一维投影函数作滤波处理，得到修正后的投影函数。然后再将此投影后的函数作反投影运算，得出原始的密度函数。 关键问题：如何修正投影函数使之在作反投影后能重建得出原始密度函数。 密度函数的二维傅立叶反变换： 1D FT 空间域 频域 1D IFT （3）卷积反投影法 卷积反投影函数改写成卷积的形式： 算法特点：将投影函数gθ(R)直接在空域中进行修正，即将gθ(R)与一个事先设计好地卷积函数C(R)进行卷积运算，然后将卷积后的结果作反投影。 卷积函数C(R) 卷积函数 因 不可积，故对 重新定义如下： 最后得 卷积函数的选择是卷积反投影方法中的关键问题。在实际的系统中选择卷积函数时还要考虑到许多其他的因素，包括系统的带宽、信噪比与分辨率等。