### SPECT图像的最大似然断层重建 #### 一、引言 SPECT（单光子发射计算机断层成像）是一种重要的医学成像技术，它通过测量体内放射性同位素发射的γ射线来生成人体内部组织的图像。传统SPECT反投影断层重建技术往往无法提供足够的细节清晰度，特别是对于那些需要高分辨率图像的应用场景。为此，研究人员开发了一种基于统计模型的最大似然断层重建技术，该技术能够显著提高图像质量，尤其是能够有效补偿随机干扰、衰减、散射等因素导致的图像退化。 #### 二、SPECT成像原理与挑战 ##### 2.1 成像机理 SPECT成像的基本过程包括：患者体内注射带有放射性核素的示踪剂，这些核素会在特定的组织或器官中积累，并以一定的概率发射γ射线。通过围绕患者旋转探测器，可以获得多个角度下的γ射线投影数据。根据这些数据，可以使用不同的算法重构出组织或器官的横截面图像。 ##### 2.2 挑战 尽管SPECT成像技术已经取得了很大的进展，但它仍然面临着几个关键的挑战： - **随机性**：探测器上接收到的γ射线数量遵循泊松分布，这增加了图像的不确定性。 - **衰减和散射**：在组织内部传播的过程中，γ射线会发生衰减和散射，这会降低图像的质量。 - **低剂量限制**：为了减少患者接受的辐射剂量，通常使用较低的放射性示踪剂剂量，这导致采集到的数据较少。 #### 三、基于统计模型的最大似然断层重建 ##### 3.1 统计模型 为了克服上述挑战，基于统计模型的最大似然估计方法被引入到SPECT断层重建中。这种方法的核心在于建立一个统计模型来描述γ射线的分布情况，并以此为基础进行图像重建。 - **泊松分布**：探测器上每个像素点接收到的γ射线数遵循泊松分布，参数λ表示该像素对应的γ射线平均数，λ与该像素处的放射性核素浓度成正比。 - **最大似然估计**：通过寻找使观测数据最有可能发生的参数值，即最大化观测数据的似然函数，来进行图像重建。 ##### 3.2 算法实现 - **重建算法**：最大似然估计的断层重建通常采用迭代算法实现，如EM（期望最大化）算法。EM算法通过不断优化似然函数来逐步逼近最优解。 - **修正的EM算法**：为了解决原始EM算法存在的问题（例如收敛速度慢、容易陷入局部最优解），研究者们提出了一些改进的方法，比如最大后验概率（MAP）和有代价的最大似然（PML）准则，以及各种修正的EM算法。 #### 四、实验结果与分析 通过对实际数据进行模拟实验，结果显示最大似然断层重建技术相比于传统的反投影法，在提高图像清晰度方面具有明显优势。这种优势尤其体现在对微小结构的检测能力上，这对于早期疾病诊断至关重要。 #### 五、结论与展望 最大似然断层重建技术为提高SPECT图像质量提供了一种新的途径。尽管这种方法在计算效率和噪声控制方面还存在一些挑战，但随着算法优化和硬件性能的提升，未来有望在临床上得到更广泛的应用。 通过综合考虑统计模型和迭代算法，最大似然断层重建不仅能够显著提高图像质量，还能有效地补偿随机干扰、衰减和散射等因素的影响，为医学成像领域带来了革命性的进步。

Derenzo.m ：提取构成 Derenzo 幻像的圆的质心和半径（需要 derenzo.png 才能工作） derenzo_main.m ：应用示例； 给出一个 .txt 文件作为输出，其中包含定义几何体和幻影来源所需的 Gate 命令。 还生成一个绘图作为最终幻影的预览。 主要参数（活性、尺寸、体模材料..）可根据需要更改。 欢迎提出建议和错误报告！

CASToR是一款核医学成像领域，基于PET的图像重建过程的C++代码程序软件。该程序代码全部开源，有相应的开发维护团队。该压缩包是改程序的安装源代码

2021-2027全球与中国PET和SPECT扫描仪市场现状及未来发展趋势.docx

SPECT图像重建与结果分析，利用matlab对SPECT生成的原始数据进行处理，重建图像，进行滤波去噪处理，并对比分析前后的图像质量。

利用典型的加权平均融合算法进行灰度或彩色多模态医学图像融合，程序具体很好的通用性，并且提供几种图像融合客观评价指标，还给出3组宝贵的已配准的待融合图像。