提出新算法,让全息成像在高分辨率下也能高效去斑。
Monte Carlo Maximum Likelihood Reconstruction for Digital Holography with Speckle
- 用随机线性代数避开高维矩阵求逆,加速最大似然估计
- 在真实孔径模型下重建质量显著提升,速度比旧方法快数倍
- 适合需要高精度与高效率的数字全息成像应用
相干成像中的散斑被统计建模为乘性噪声,严重阻碍图像重建。最大似然估计(MLE)虽具理论优势,但因有限孔径数字全息系统中高维矩阵求逆成本过高,尤其在高分辨率下难以应用。本文提出基于随机线性代数的可扩展MLE优化方法,无需显式矩阵求逆即可计算梯度。通过利用传感矩阵结构并结合共轭梯度法评估似然梯度,所提算法支持精确的孔径建模,无需简化假设。该方法称为投影梯度下降结合蒙特卡洛估计(PGD-MC)。实验表明,PGD-MC在多种物理真实孔径模型下均表现稳健,显著提升重建质量与计算效率,并可有效扩展至高分辨率数字全息。使用三种典型去噪器作为正则化项的对比实验显示,其在精度和速度上均优于现有插件式模型迭代重建方法。代码已开源:https://github.com/Computational-Imaging-RU/MC_Maximum_Likelihood_Digital_Holography_Speckle。
原文摘要 · Abstract (English)
In coherent imaging, speckle is statistically modeled as multiplicative noise, posing a fundamental challenge for image reconstruction. While maximum likelihood estimation (MLE) provides a principled framework for speckle mitigation, its application to coherent imaging system such as digital holography with finite apertures is hindered by the prohibitive cost of high-dimensional matrix inversion, especially at high resolutions. This computational burden has prevented the use of MLE-based reconstruction with physically accurate aperture modeling. In this work, we propose a randomized linear algebra approach that enables scalable MLE optimization without explicit matrix inversions in gradient computation. By exploiting the structural properties of sensing matrix and using conjugate gradient for likelihood gradient evaluation, the proposed algorithm supports accurate aperture modeling without the simplifying assumptions commonly imposed for tractability. We term the resulting method projected gradient descent with Monte Carlo estimation (PGD-MC). The proposed PGD-MC framework (i) demonstrates robustness to diverse and physically accurate aperture models, (ii) achieves substantial improvements in reconstruction quality and computational efficiency, and (iii) scales effectively to high-resolution digital holography. Extensive experiments incorporating three representative denoisers as regularization show that PGD-MC provides a flexible and effective MLE-based reconstruction framework for digital holography with finite apertures, consistently outperforming prior Plug-and-Play model-based iterative reconstruction methods in both accuracy and speed. Our code is available at: https://github.com/Computational-Imaging-RU/MC_Maximum_Likelihood_Digital_Holography_Speckle.
Thank you to arXiv for use of its open access interoperability. PaperDance 不是 arXiv 官方产品;中文卡片由大模型生成,请以原文为准。