基于柱谐函数分解法约束EM迭代初值的中子源三维图像重建方法

本发明涉及一种脉冲中子源图像诊断方法，具体涉及一种基于柱谐函数分解法约束EM迭代初值的中子源三维图像重建方法，解决现有柱谐函数分解法和迭代类方法均难以获得高精度、高准确性中子源三维图像的技术问题。该一种基于柱谐函数分解法约束EM迭代初值的中子源三维图像重建方法，包括以下步骤：1】获取N个中子源二维投影图像I0，并确定其对应的投影角，N＞1；2】获得N个中子源二维投影图像I；3】获得柱谐函数重建结果；4】利用EM算法迭代更新，直至第k次迭代结果与第k‑1次迭代结果的差异小于10‑4时停止迭代，并将最后一次迭代结果作为输出结果，k≥1；5】将输出结果重组为三维矩阵，并将重组后三维矩阵作为中子源重建结果。