直接解调成象方法在核医学图象重建中的应用.pdfVIP

直接解调成象方法在核医学图象重建中的应用 刘 力 （中国科学院高能物理研究所 北京 100039) 本文受国家攀登计划（B ）基金资助 ABSTRACT ： 本文将一种天文学成象的直接解调法用于核医学中的图象重建，就最简单的 一个切片时的蒙特卡罗样本进行计算分析，结果显示了该方法对小尺度放射源具有较高的灵 敏度和空间分辨率。同其他传统方法相比，有一定的优越性，可应用于核医学图象数据的分 析，以及对微小“点状”癌肿的早期诊断。 关键词 直接解调成象法 图象重建 非线性约束 一、引言 目前先进的核医学影象设备可以对人体进行断层成象。如 SPECT，PET 等，其空间分 辨率可分别达到 10mm,3mm 左右。影响空间分辨率的因素很多，如探测器尺寸，排列方式 等。从软件角度，传统的重建方法如滤波反投影法[1] 等有最大的截止频率ν ，从而无法分 max 辨比 1/ ν 更细小的结构。另一个重要因素就是统计涨落比较严重。病变小尺寸组织可能 max 有很高的放射性药物浓缩，但由于尺寸小，仍是一个低信噪比的信号。 在高能天文学中，也常遇到信噪比小的点源的探测与图象重建问题。我国天体物理工作 者近来提出一种带有非线性约束的直接解成象方程的方法[2 ] 。实践表明了该直接方法 （DM ） 比传统方法灵敏度高，空间分辨率好，尤其可用于统计性差的观测数据[3] 。 本文先简述DM 方法，并将之推广应用到核医学图象重建中。就一个切面的简单情形， 实现了由低维观测数据对高维图象的重建，并显示了该方法卓越的空间定位与分辨能力。 二、直接求解成象方程的成象方法 以高能天体物理观测为例，联系于观测数据与待求真实分布之间的成象方程为[2 ] ： ∑P f (i,j ) d (α,β) α,β=1,2,m x ,m y 简记为 Pf=d ( 1) αβ, ,i , j i , j 式中 f(i,j) , i,j=1, 2,n x ,n y 为第 （i ,j )个象素上的未知量 （比如放射性浓度，辐射强度等）， d (α,β) 为第( α,β)道上的观测数据（计数），P 为探测器的点扩展函数，也可理解为 αβ,ij 1 第（i ,j )个象素对第 ( α,β) 道上计数的贡献比，是已知的。 观测数据中存在统计噪声,方程 （1）是矛盾的和严重病态的，通常不能直接求解。在付 氏分析频率空间中看，完整地求解方程（1）就是设计一个倒数滤波器来进行高频恢复。但 噪声含有大量高频成份，故由此得到的 “严格解”是噪声占优势的解。对分立的点状结构的 信号源，它含有相当多的高频成份，这使得一切基于付氏分析的滤波器此时不能有效工作。 DM 方法的特点就是能对突变信号源进行有效的恢复，同时又能抑制噪声。 根据文[2]提出的直接方法，方程（1）作相关变换（拟合）而成为： ∑P ij ,i j f (i ,j ) d (i, j ) i,j , i,j =1,2 nx ,ny (2) ij 其中 d (i,j ) ∑P d (αβ, ) P ∑P P (3) αβ，ij ij ,i j αβ,ij αβ,i j αβ αβ 这样可改善迭代求解方程的收敛性，压抑统计涨落，控制虚假峰的出现。用迭代法求解方程 （2 ）。迭代过程中加入位形空间中的物理约束条件，这些条件具有物理上的直观性，准确 和简