基于微元法的高精度系统响应矩阵建模.pdfVIP

第31卷 第 10期 核 电子学与探测技术 Vo1．31 No．10 2011年 10月 NuclearElectronics＆ DetectionTechnology Oct． 2Oll 基于微元法的高精度系统响应矩阵建模 邓贞宙 ，段志文 ，谢庆国 ， (1．武汉光电国家实验室(筹)，湖北武汉430074；2．华中科技大学生命科学与技术学院， 湖北武汉430074；3．华中科技大学数学与统计学院，湖北武汉430074) 摘要：系统响应矩阵是以矩阵形式表征的从激励到响应的数值化关系。由于绝大多数正电子发射 成像的重建算法都包含正投影和反投影的计算步骤 ，系统响应矩阵精度是研究者普遍关心的重要因素 之一。通过对系统响应矩阵进行微元法建模，提高了系统响应矩阵的精度。数值结果显示，随着微元的 尺寸变小，64点微元正投影的信噪比提高了100倍以上，16线微元反投影的信噪比提高了7个数量级。 关键词：微元法；系统响应矩阵建模；PET图像重建 中图分类号：TB867 文献标识码：A 文章编号： 0258-0934(2011)10．1125-06 系统响应矩阵(SystemResponseMatrix)是 子工艺，使得计算机的存储器读写速度和存储 正 电子发射 成像 (PositronEmissionTomo— 器体积大幅提高。于是在计算系统响应矩阵 graphy，以下简称 PET)中表征活度分布 (激励) 时，实现手段分化成了两类：一类是非外存储型 到符合计数(响应)之间的量化关系，可以用 的计算手段l3一’J，即实时计算系统响应矩阵，计 =g表示，其中向量，是像素个数为 N的图像， 算出系统响应矩阵部分后，不进行硬盘存储而 向量g是 条响应线上探测到的符合计数值。 直接进行正＼反投影操作，计算完成后亦不予以 h则是系统响应矩阵，用行数为M，列数为N的 保留；一类是预存储式 的计算手段 ，即事 矩阵表示。系统响应矩阵的精度直接关系到正 先计算好系统响应矩阵，存人硬盘设备中。在 投影和反投影的精度，从而影响图像重建的精 使用系统响应矩阵进行正反投影操作时直接读 度。在与PET成像数学模型相似的其他成像 取硬盘中的系统响应矩阵。预存储的系统响应 中(例如CT、SPECT、MRI)，系统响应矩阵同样 矩阵在第一次计算时可以采用解析 的计算 是普遍关心的问题 】J。 方式，也可以采用仿真软件 ¨或者迭代 ¨ 现有的系统响应矩阵计算方法有很多种。 的计算方式。由于预存储式的系统响应矩阵只 应用在这些方法时，计算精度和计算速度是相 决定于存储器的读写时间，应用预存储系统响 互掣肘的两个因素-】J。精确度高的建模方法 应矩阵进行正反投影的耗时对建模的复杂度不 计算复杂，耗时较长；计算简便，耗时较短的建 敏感。对于不同的系统响应矩阵建模方式，可 模方法其精确度较低。近年来高速发展的微电 以根据具体应用环境决定是否进行事先计算系 统响应矩阵。当计算复杂、图像像素较多、响应 收稿 日期 ：2011—03—29 线条数多时，选择后者 ；当实现简单、图像像素 基金项目：国家自然科学基金，科技部国 较少、响应线条数少时可以选择前者。尤其是 际合作重点项 目(2009DFR30580)。 当读取全部系统响应矩阵的时间比计算系统响 作者简介：邓贞宙(1987一)，男，江西南昌人，华中科 应矩阵的时间少时，更应该选择预存储的计算 技大学在读博士研究生，从事不