蒙特卡罗方法在针孔SPECT图像重建中应用.pdfVIP

第三篇蒙特卡罗方法应用 蒙特卡罗方法在针孑LSPECT图像重建中的应用 代秋声，漆玉金 (中国科学院上海应用物理研究所，上海201800) 【擅薹】 本文探讨用蒙特卡罗方法来模拟针孔SPECT成像，提供较精确的投影矩阵用于三维图像霞建，获得 高分辨的重建图像。我们采用蒙特卡罗方法模拟计算每个体素与像素之间的物理映射关系，建立各角度的投影矩 阵，用于基于统计迭代的OS．EM图像重建算法里，实现高分辨的图像重建；然后，用标准的Jaszczak模具来进行 实验验证与评估。初步的实验结果表明：跟理想的针孔几何映射计算提供的投影矩阵相比，南蒙特卡罗方法提供 的投影矩阵能够显著地提高重建图像的整体质量，并改善其空间分辨率。因此，利用蒙特卡罗方法提供的投影矩 阵来实现高分辨的针孔SPECT图像重建是可行的。 [关键词】 蒙特卡罗方法；针孔SPECT：图像霞建 【中圈分类号】0411．3；0571．33【文献标识硼 A 0引言 由于针孔成像技术具有放大功能，可以有效克服探测器内在分辨率对系统窄问分辨率的制约，已经被 (maximumlikelihood subsets expectation expectation 修正，显著改善图像质最，而且收敛速度快，因此得到广泛应用p一。 矩阵之间的映射关系，其数学描述被称为投影矩阵。然1fi『，在针孔成像的几何映射模型中，针孔通常是被 作为理想的点来处理的，而实际的针孔都具有一定的几何大小，点源经过针孔后会在探测器I：形成一个放 大的光斑图像(见图1)，其数学描述被称为点源响应函数(PSF)。因此，当用简单的几何映射来代替复杂 的物理映射时，不可避免的在投影矩阵中引入了大量的误差，从而影响了重建图像的质量。 为了继续提高针孔SPECT的空间分辨，我们需要准确地知道体素矩阵与像素矩阵间的物理映射关系， 即各个体素的PSF矩阵。PSF矩阵的获取通常有两种方法，一是实验测型7】，二是理论计算。前者虽然准 Carlo 确性较高，但是操作复杂，较难实现；后者虽然简瞽易得，但是准确性不高。此处，我们采用Monte 模拟方法来获取PSF矩阵，组建投影矩阵，应用于图像重建。 1针孔成像的物理模型 下面介绍针孔成像的物理模型。图1是针孔成像的示意图，其主要f}l成像物体、针孔准直器和探测器 三部分组成。伽玛光子针孔成像的基本过程是：注入生物体内的放射性核素衰变放出的伽玛光子向各个方 向等概率发射，满足针孔准直器几何条件的伽玛光子，通过针孔到达探测器，在探测器表面投影出伽玛源 的强度分布图。在此过程中，准直器是通过吸收来滤掉不满足几何条件的伽玛光子。根据射线与物质的相 陲金璜目I国家自然科学基金(批准号：上海市“漓江人才计划”资助项目 I作者简介l代秋声(1977一)。男，湖北宜城人，博士研究生，从事针孔SPECT成像研究 332 ∞I目全目荣特p罗^法A其＆用学术含论立 互作用规律．不可避免的会有一部分伽玛光子通过透射或散射的方式到达探测器，从而影响了引孔成像的 卒问分辨率。 目I孔^*i自日 l FigD…ofpinholeimaging 2蒙特卡罗模拟 为了获取各个体素的PSF矩阵，我们采用蒙特}罗方法米模拟针孔成像的物理过程，且体方法如下： 首先根据预期的体素尺、j_将成像物体离散成j维体素矩阵，然后根据成像的物理模型依次对母个体素进行 投影模拟计算。在对体素进行投影模拟时．首先在体系所处的几何卒问中依照均匀分布原则抽样产生伽玛 光f的初始位霄，再根据辑向同性原mⅡ抽样产生韧始发射方|f町，然后对能够选到引扎准直嚣空问范围内的 伽玛光干跟踪记录。 在使用蒙特#罗方法进行模拟计算过科中，由于体素数最巨人．选I旷嚣级，共至皿高，冈此需要埘 模扭过程进行简化