图像重建与可视化.docVIP

第七章 医学图像重建与可视化 数字图像处理在绝大多数情况下， 是对已获得的图像进行处理和转化， 即把处理技术作为系统， 其输入和输出均为数字图像。 因此在医学领域的应用中， 数字图像处理也习惯性的称作后处理。 但在某些情况下， 图像处理过程也涉及将数据进行计算并最终转化为图像的过程， 如图像重建，这种情况更多的出现在医学领域中，如 CT、MR 、超声成像和核医学成像。 这些医学成像系统往往通过获得尽可能多的目标的检测数据， 并将数据进行计算处理， 最终生成图像， 这就是图像重建 ( image reconstruction) ，基本过程如图 7-1 所示。这种情况发生在图像后处理之前。 图 7-1 图像重建过程 此外，上述的医学成像系统， 往往获得连续的人体断层数据和断层图像。 如 果将这些数据或图像进行整合处理， 仅通过图像处理技术获得新的切层位置和方 向的断层图像； 或将断层图像合成为体数据并显示成三维图像， 以加强人们对器 官解剖结构和病灶三维形态的观察和理解， 这种处理我们称为图像的三维可视化 ( three- dimensional visualization) 或三维重建 ( three- dimensional reconstruction) ， 其处理过程如图 7-2 所示。 图 7-2 图像三维可视化过程 图像重建与图像可视化都是较为复杂的计算与处理过程， 本章仅对医学图像的重建与可视化的基本理论与方法进行讨论。 第一节 医学图像重建概述 一、医学断层成像 在医学断层成像出现前， 传统的医学成像技术或手段是将成像区域内的三维 人体组织，投射于二维的成像范围内， 记录在胶片或显示屏上， 最终形成二维医学影像。这种方式势必造成人体组织信息在影像上的重叠和遮挡， 使病灶缺失某一维度（如深度）的信息，虽然可以通过多体位摄影进行适当补偿，但并不能根本解决问题。 如：在胸部影像的实际诊断过程中， 某些肺部病灶由于肋骨影像的重叠而无法确定位置，或者被肋骨或纵隔影像模糊。 这种影像重叠是三维影像在一维路径上产生了叠加，用积分形式可以表示 为： I d  x, y  I 0 exp  x, y, z dz  （ 7.1 ） I d  x, y  为影像记录装置上记录的  X 线强度分布，  x, y, z  为组织吸收  X 线系数 的分布函数， I 0 为入射 X 线强度。 为解决这一问题，曾出现了体层摄影术（ X - ray tomography），光源和影像记录装置（胶片或成像板） 沿相反方向运动， 但保持与成像区域的垂直距离不变，如图 7-3 所示。这样成像区域中的某层组织处于聚焦面上， 被记录于影像记录装置从而成像，其余层面组织影像由于受到光源和影像记录装置的运动而模糊， 无法清晰显示。 图 7-3 直线运动体层摄影 体层摄影根据光源的运动形式分为直线运动成像、 圆运动成像和摆线运动成像。这些运动方式要求光源和记录装置精确计算运动速度， 以保证计算准确。 同时严格限制运动时间，保证在被检者可接受的屏气时间内完成， 以限制呼吸伪影。 体层摄影成像的计算如下：探测器接受强度为 I d xd , yd t xd , yd I 0 xd , yd （ 7.2 ） M M L k m k m k m 直线运动形式的体层摄影射线源移动路径为： f x, y g x, y, z d vt rect x （ 7.3 ） y X 代入（ 7.1 ）式得最终的像强度为： I d xd , yd t xd , yd I 0 rect xd yd （ 7.4 ） M M X k m X k m 其中， V 是射线源在 X 方向上走过的速度， X 是射线源走过的距离， m 为放大率， k 为底片运动的换算因子， k+m 代表总放大率。 体层摄影曾一度作为一项重要的成像功能出现在 X 线成像设备中。但体层成像存在一些无法解决的问题，如辐射面积广、剂量大，要成一层组织的影像，则比成像区域大的多的组织体积均要接受 X 线照射，大大增加了被检者受照射量，增加了辐射损伤。另外，体层摄影的图像质量，尤其是对比度分辨率，相对 X 线平片并无明显提高。 影像信息中还掺杂了模糊区域的半影， 相当于影像中增加了更多的散射线，削弱了组织对比。体层摄影还需要被检者的屏气与制动配合，不适于床旁摄影和屏气困难的被检者。 真正的断层成像出现在 1971 年，始于世界上第一台可应用于临床的 CT，安装于 Atkinson- Morley 医院。这种技术也叫 CAT（computed aided tomography computed axial tomography），其图像重建的数学理论最早适于奥地利数学家 Radon 于