图像处理课件-chapter6汇编.pptVIP

(数字)图像处理 (Digital) Image Processing 第6章 图像重建 6.1 计算机断层扫描技术 6.2 投影定理 6.3 傅立叶投影定理 6.4 卷积逆投影重建 6.5 代数重建 6.6 三维图像重建的体绘制 发射断层成像系统 反射断层成像系统 透射式断层扫描成像系统 发射断层成像： 发射源在物体内部，将具有放射性的离子（放射元素）注入物体内部，在物体外部检测其经过物体吸收之后放射量。 发射投影成像 如，正电子发射成像（PET：Positron Emission Tomography） 采用在衰减时放出正电子的放射性离子，放出的正电子很快与负电子相撞湮灭而产生一对相背运动的光子。 相对放置的两个检测器接收到这两个光子就可以确定一条射线， 检测器围绕物体呈环形分布， 相对的两个检测器构成一组检测器对，检测由一对正负电子产生的光子。 射线投影成像的基本原理： 人体组织对X射线吸收和散射，造成衰减， 人体内的不同结构，比如脂肪、胰、骨骼对X射线吸收能力有所不同。 透射投影成像， 图6.2表示等强度的射线透过不同密度分布时的情况， 每块上的数字表示每块的密度或衰减，总的衰减是叠加的， 其中一条射线束通过均匀密度物质的厚块，另一射线通过不等密度的厚块组合，但检测器的记录相同， 因此，投影重建时需要一系列投影才能重建二维图像。 投影重建是利用人体(物体)对射线的能量吸收衰减作用,不同密度的组织具有不同的吸收能力。 6.1 计算机断层扫描技术 CT扫描成像示意图 第一代CT 第二代CT 第三代CT CT一次平移扫描所获得的输出信号 第四代CT 第三和第四代CT CT 扫描仪 数据获取系统 (DAS) 数据获取系统(DAS) 医学影像领域： Computed Tomography(CT)：获1979年诺贝尔奖（Nobel Price） 布尔赫、珀塞尔，获1952年诺贝尔奖，发现了核磁共振现象 劳特布尔（美）、P·曼斯菲尔德（英）获2003年年诺贝尔奖，核磁共振的研究 6.2 投影定理 一个N维函数f(x1,x2,x3,…,xN)在第N－1维上的映射称为函数f在第N-1维的投影。在简单的二维情况下，函数f(x,y)在x轴上（沿y方向）和在y轴上（沿x方向）的投影可分别表示为 : 设f(x,y)的傅立叶变换为F(u,v)，则根据傅立叶反变换式可知： 把f(x,y)傅立叶变换代入： 上式表明gy(x)是F(u,0)的傅立叶反变换。或者说gy(x)的傅立叶变换G(u)与F(u,0)相同。由此可知，函数f(x,y)在x轴上投影的傅立叶变换等于f(x,y)的傅立叶变换在(u,v)平面上沿u轴平面上的切片。 沿y轴的投影图示 假设函数f(x,y)投影到一条经过旋转的直线上t1，t是一条与t1平行经过原点的直线，与t垂直经过原点的直线为s，该直线s与x轴的夹角为θ，直线t1离开原点的距离为s1 函数f(x,y)沿着s方向(在t1轴上)的投影为 上式表明，f(x,y)在一条与x轴夹角为θ，离开原点距离为r的直线上的投影的傅立叶变换等于二维傅立叶变换在与u轴成θ方向上的切片，这就是投影定理，也称之为切片定理。 6.3 傅立叶投影重建 傅立叶投影重建的基础就是傅立叶投影定理。 根据投影定理，如果能将不同角度θ1,θ2,…, θn得到的投影值进行傅立叶变换，就可以得到F(u,v)分别在相应角度位置上的切片。当切片趋向无穷多，即取无穷多个投影时，就可获得在 (u,v)平面上的所有F(u,v)值，从而进行傅立叶反变换就可以重建图像f(x,y)。 由此可得，用傅立叶变换法重建图像的步骤如下： ① 根据式（6.12）或式（6.18）对N个不同θ方向上投影进行一维傅立叶变换。 ② 在傅立叶变换空间从极坐标向直角坐标插值。 ③ 利用式（6.15）或离散形式的傅立叶频谱进行反变换得到重建图像。 6.4 卷积逆投影重建 卷积逆投影重建法：以投影切片定理为基础； 傅立叶变换重建法：计算量比较小，但要二维插值，在射电天文学研究中得到应用广泛； 卷积逆投影法：能快速实现，噪声较小时可重建出准确清晰的图像，在X射线CT成像中应用广泛。 在计算投影的一维傅立叶变换F(R,θ)时，R为频域极轴变量。投影数据g(ρ,θ)总是被有限截断。当ρ的取样间隔为d时，在频率R的变化范围将是-d/2～d/2，于是投影切片定理可近似写成: 由(6.15)可知： 由式（6.21）可知，右边正是投影数据g(ρ,θ)与脉冲响应h(ρ)所表示的滤波器的卷积，h(·)为卷积函数。 求f ’(x,y,θ) 则是在θ角方向上卷积了的投影，因此从式（6.22）求f (x,y)可被认为是求逆投影过程，即卷积逆