医学图像重建算法概述.pptVIP

CT图像重建算法概述 内容目录 图像重建原理及比尔定律 图像重建过程及相关概念 解析算法-滤波反投影FBP 迭代算法-ARTOSEM CT图像重建原理：相同射线强度，穿过不同物质，衰减不同，利用这一规律区别人体内的不同物质。 骨骼 软组织 空气 初始强度 比尔定律：描述X射线穿过物体时发生衰减的规律，如下： 均匀物质，衰减系数μ 初始强度 穿透强度 长度L 决定了是什么物质 u1 u2 u3 u4 un d d d 实际情况：物质并非均匀，组成复杂，因此，我们将物质分成许多小份，即有了体素（像素）的概念，如下d： 根据前面公式，即有： 初始强度 穿透强度 μ1+μ2+μ3+ μ4+μ5+μ6 …+μn μ11 μ12 μ13 μ14 … μ21 μ22 μ23 μ24 … μ31 μ32 μ33 μ34 … … … … … μ11+μ12+μ13+…+μ1n=… μ21+μ22+μ23+…+μ2n=… μ31+μ32+μ33+…+μ3n=… μ41+μ42+μ43+…+μ4n=… μ51+μ52+μ53+…+μ5n=… μ61+μ62+μ63+…+μ6n=… … … 被扫描物体 45 1 456 2 101 23 4 5 22 67 2 334 1 43 111 12 45 45 88 1 445 3.3 34 134 23 4 23 1 2 345 2.3 111 0 1 90 23 33 2 456 2 23 11 89 11 11 5 1 678 145 21 7 0 3 12 0 3 12 444 56 9 4 4 45 12 0-10 红色 10-50 蓝灰色 50-100 灰蓝色 100-200 浅蓝色 200 蓝色 CT值：实际情况，重建图像显示的并非衰减系数，而是用CT值描述，表示如下： 常用： 人体组织 HU 人体组织 HU 骨密质 1000 脂肪 -100 钙质 60 水 0 脑白质 36 血液 16 脑灰质 24 凝固的血 56-70 图像重建过程及相关概念： Central Park 投影 图像 反投影 我就是CT 扫描原始数据（即投影） 医学图像重建算法之——解析法。 解析法是以中心切片定理为基础的反投影方法，典型算法如滤波反投影法(Filtered Back-Projection，FBP)。解析法具有容易实现,速度较快,且能重建出高质量的图像的特点,但是对投影数据完备性要求高。 医学图像重建算法 3 2 4 0 5 4 7 2 解析算法之——反投影示例（数据完备） 假设，已知： 7 2 7 2 7 2 5 5 4 4 5 4 + + 3 6 6 0 3 2 4 3 + 3 6 18 15 21 9 6 3 = 9 6 12 0 减去基数9 3 2 4 0 图像还原 3 2 4 0 5 4 7 2 解析算法之——反投影示例（数据欠缺） 假设，已知： 7 2 7 2 7 2 5 5 4 4 5 4 + 12 7 11 6 = 3 -2 2 -3 减去基数9 差距较大 迭代重建算法简述 医学图像重建算法之——迭代法。 迭代法迭代法是利用求解线性方程组来重建图像，实际是从—幅假设的初始图像出发，采用逐步逼近的方法，将理论投影值同实际测量投影值进行比较，在某种最优化准则指导下寻找最优解。它能够在投影数据信噪较低条件下，获得高质量图像。 8 9 12 5 X1+X2=5 X3+X4=12 X2+X4=9 X1+X3=8 X2+X3=7 X1+X4=10 迭代算法简单示例 10 7 X1=3 X2=2 X3=5 X4=7 A）假定一初始图像； B）计算该图像投影； C ) 同测量投影值对比； D )计算校正系数并更新初始图像值； E ) 满足停步规则时，迭代中止，否则以新的重建图像作为初始图像从B步开始。 实际操作中迭代重建算法求解过程 典型迭代重建示例 5 10 15 20 15 35 25 20 25 30 这个图像不知道 知道这些 10 10 10 10 20 20 20 20 20 20 猜测图像 0.75 1.75 1.25 1.00 1.25 1.50 9.375 14.0625 21.875 32.8125 1次 迭代 反投影 20 20 20 20 20 20 0.75 1.75 1.25 1.00 1.25 1.50 5 10 15 20 15 35 25 20 25 30 知道这些 10 10 10 10 20 20 猜测图像 0.75 1.75 7.5 7.5 17.5 17.5 1次迭代 第一子集 0.75 1.75 这个图像不知道 有序子集迭代方法示例 7.5 7.5 17.5 17.5 6 9 14 21 二次迭代 第一子集 0.8 1.2 25 25 5 10 15 20