正电子发射断层扫描仪(Positron Emission Tomography, PET)是当前医学界公认的肿瘤、心脏、脑等疾病诊断与病理生理研究的重要方法。随着核医学影像设备的广泛应用和计算机技术的迅速发展，图像重建方法作为PET成像的一个关键环节，其研究工作也越发受到重视。 PET探测器检测注入人体的示踪剂在湮灭辐射过程中产生的射线，经过符合采集系统处理形成投影线，以SINO的方式存放于计算机硬盘中[1]。计算机调用图像重建模块，生成人体断层图像。目前，PET图像基础重建算法主要包括解析法和迭代法。 1. 解析法 解析法是以中心切片定理为基础的反投影方法，常用的是滤波反投影法(Filtered Back-Projection, FBP)。在FBP中，图像重建主要包含两个步骤：反投影和滤波。 我们在初中就已经学过投影与反投影的概念，从不同角度观察物体可以得到不同的信息，当我们从多种不同角度获取物体的投影，可以反向推出这个物体真实的形态。 图1 光线将物体的形状投射到一个平面称为投影 在成像原理上，PET和CT略有差异。CT是投射成像，X射线旋转360°，采集被扫描物体不

医学图像重建算法中OSEM算法的Matlab程序

有序子集期望最大化(OSEM)迭代算法是近年来发展较快的一种迭代类算法。但该算法在迭代过程中容易产生条纹状伪影、金属伪影或者散射伪影。构造了平滑约束矩阵作为先验信息引入到重建迭代过程, 建立了一种平滑约束OSEM(SC-OSEM)迭代重建算法。分别将中值滤波、全变差最小(TVM)方法作为平滑约束条件, 通过数值模拟, 针对不完备理想投影数据、含金属不完备投影数据、含噪声不完备投影数据三种情况, 重建出了与原始模型一致性较好的计算机层析成像技术(CT)图像, 比单独OSEM迭代算法重建质量高, 并且发现中值滤波约束重建图像的整体噪声较小, TVM算法使金属边界更清晰, 表明SC-OSEM迭代重建算法是一种精度高、适应性较强的CT重建算法。

基于优化子集顺序的三维OSEM图像重建算法

主要实现了基于正弦图的三种重建算法：FBP,MLEM,OSEM等，适合从事医学图像开发研究人员，包括但不限于CT, MRI，PET等。