核电子探索者

作为我们的原始图像，我们将使用 Shepp-Logan 幻影。在计算拉东变换时，我们需要确定要使用多少个投影角度。经验法则表明，投影的数量应该与对象横跨的像素数量大致相同（要了解原因，请考虑在重建过程中必须确定多少个未知像素值，并将此值与投影提供的测量值进行比较），我们在这里遵循此规则。

OSEM（Ordered Subset Expectation Maximization）是一种在医学成像中常用的迭代图像重建算法，特别是在PET（正电子发射断层扫描）成像中。函数内部，我们初始化了一个重建图像（recon），并计算了正弦图的归一化权重。然后，我们进行了多次迭代，在每次迭代中，我们更新了重建图像的估计。此外，还可以使用专门的库，如PyTomography ，它提供了更完整的实现和优化。然后，我们创建了一个示例“phantom”图像和一个对应的正弦图，最后使用OSEM算法进行图像重建。

【代码】OSEM算法实现——Jupyter笔记（2）

【代码】OSEM算法实现——Jupyter笔记（1）

【代码】PET图像重建MELM算法。

【代码】PET图像重建OSEM算法。