基于分数匹配生成模型的无透镜成像方法

基于分数匹配生成模型的无透镜成像方法 Lens-less imaging via score-based generative model first-author 伍春花（1997-），女，江西南昌人，硕士研究生。现就读于南昌大学电子信息工程系，主要从事计算光学成像方面的研究。E-mail： wuchunhua@email.ncu.edu.cn http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34991825&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34991831&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34991828&type= 22.01332855 32.00399780 伍春花 （1997-），女，江西南昌人，硕士研究生。现就读于南昌大学电子信息工程系，主要从事计算光学成像方面的研究。E-mail： wuchunhua@email.ncu.edu.cn wuchunhua@email.ncu.edu.cn corresp E-mail： wanwenbo@ncu.edu.cn E-mail： wanwenbo@ncu.edu.cn 万文博（1990-），男，甘肃平凉人，博士，讲师。现就职于南昌大学信息工程学院电子信息工程系，主要从事生物医学光子学及光学计算成像等方面的研究。Email： wanwenbo@ncu.edu.cn http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34991834&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34991840&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34991837&type= 22.01332855 32.00400543 万文博 （1990-），男，甘肃平凉人，博士，讲师。现就职于南昌大学信息工程学院电子信息工程系，主要从事生物医学光子学及光学计算成像等方面的研究。Email： wanwenbo@ncu.edu.cn wanwenbo@ncu.edu.cn 南昌大学 信息工程学院，江西 南昌 330031 南昌大学 信息工程学院，江西 南昌 330031 School of Information Engineering， Nanchang University， Nanchang 330031， China School of Information Engineering， Nanchang University， Nanchang 330031， China 无透镜成像受到同轴全息图中孪生像噪声的影响，一直面临着重建信噪比差和成像分辨率低的问题。针对该问题，本文提出一种基于分数匹配生成模型的无透镜成像方法。在训练阶段，通过连续随机微分方程（Stochastic Differential Equation， SDE）缓慢添加高斯噪声扰动数据分布，然后训练具有去噪分数匹配的连续时间相关的分数函数，用于求解反向SDE生成目标样本数据。在测试阶段，使用单张菲涅尔波带片作为掩膜，在非相干光照明下实现无透镜编码调制，然后使用预测-校正的方法在数值求解器SDE和数据保真项步骤之间轮换更新进行图像重建。在LSUN-bedroom和LSUN-church数据集上的验证结果表明，提出的算法能够有效消除孪生像噪声，峰值信噪比和结构相似性分别可达25.23 dB和0.65。与传统的基于反向传播和基于压缩感知的无透镜成像结果相比，峰值信噪比分别提高17.49 dB、7.16 dB，结构相似度分别提高0.42、0.35，从而实现图像重建质量的有效提升。 Lens-less imaging is affected by twinning noise occurring in in-line holograms， and the reconstructed results continuously face poor reconstruction signal-to-noise ratio and low imaging resolution. This study proposes a lens-less imaging via a score-based generation model. In the training phase， the proposed model perturbs data distribution by gradually adding Gaussian noise by using a continuous stochastic differential equation （SDE）. A continuous time-dependent score-based function with denoising score matching is then trained and used to solve the inverse SDE required to generate object sample data. In the testing phase， a single Fresnel zone aperture is used as a mask to achieve lens-less encoding modulation under incoherent illumination. The prediction-correction method is then used to alternate iteration steps between the numerical SDE solver and data-fidelity term to achieve lens-less imaging reconstruction. Validation results on LSUN-bedroom and LSUN-church datasets show that the proposed algorithm can effectively eliminate twin image noise， and the peak signal-to-noise ratio （PSNR） and structural similarity （SSIM） of the reconstruction results can reach 25.23 dB and 0.65， respectively. The PSNR values of the reconstruction results are 17.49 dB and 7.16 dB， which is higher than that of lens-less imaging algorithms based on traditional back propagation or compressed sensing， respectively. In addition， the corresponding SSIM values were 0.42 and 0.35 higher， respectively. Therefore， the reconstruction quality of the lens-less imaging is effectively improved. 无透镜成像 编码成像 生成模型 图像重建 lens-less imaging encoding imaging generative model image reconstruction 1 引 言 1 基于掩膜的无透镜成像系统通过在图像传感器前几毫米处放置一个薄掩膜来替代镜头组件，直接应用图像传感器（如CCD和CMOS）记录场景信息，从而简化光学硬件 ［ 1 1 - 3 3 1-3 ］ 。该成像方式打破了场景到图像一一对应的采样方式，将成像的重心由硬件转移到计算上，使三维荧光显微镜、热成像和重聚焦摄影等成像系统实现结构紧凑化、轻量化、低能耗和低成本成为可能。在过去几年里，研究人员提出了各种各样的基于无透镜相干系统的成像技术，如，相干衍射成像 ［ 4 4 - 5 5 4-5 ］ 和片上显微镜 ［ 6 6 - 7 7 6-7 ］ 等。然而，这些成像技术对光的相干性要求很高，极大地限制了成像系统应用。 菲涅尔非相干全息术是新发展的一种非相干光成像技术。非相干光成像技术最早起源于20世纪90年代，Mertz和Young ［ 8 8 ］ 提出了用菲涅尔波带片（Fresnel Zone Aperture，FZA）编码全息成像的建议，将全息成像的概念和应用范围拓展到了非相干光领域。与相干成像技术相比，非相干成像技术的照明光源获取容易，且记录过程不会产生类似于激光全息术中固有的散斑噪声，可以有效提高再现图像的质量。因此，非相干技术在荧光三维成像、非相干全息三维显示、自适应光学、天文成像等领域具有独特的应用优势 ［ 9 9 ］ 。在无透镜成像系统中，使用非相干光照明将场景中无数个点光源线透过FZA，这些点源投影在传感器上非相干叠加后最后形成了编码图像。理论上，编码图像既包含物体振幅/透射率的信息，又包含相位延迟的信息。由于图像传感器只对入射光强度敏感，在图像采集过程中只能采集图像的强度信息而丢失了相位信息，这导致重建图像中出现物体的孪生像 ［ 10 10 - 12 12 10-12 ］ 。孪生像是同轴全息中的固有问题，表现为重建图像中会有物体的共轭像叠加在原始图像上，从而导致重建图像质量下降。 为克服无透镜成像中的孪生像效应，通常采用四步相移法来去除孪生像。该方法需要采集经过四组不同相位波带片的编码图像 ［ 13 13 - 14 14 13-14 ］ ，在实际操作中，不仅需要更换波带片，而且需要图像间的配准，增加了实验的复杂性。为了弥补这个缺点，可以采用单张FZA实现无透镜成像系统的编码，这样既不需要波带片和传感器之间的校准，又能简化实际操作。然而，使用单张FZA编码方式并不能消除孪生像，因此需要通过重建算法消除孪生像的影响。 针对无透镜成像过程因丢失相位信息导致重建图像不唯一的问题，如何建立场景与图像之间的联系，并通过求解逆问题反演图像是解决编码成像问题的关键。传统方法是直接通过使用反向传播（Back Propagation，BP）算法对同轴全息图进行重建，然而孪生像干扰严重影响了重建图像。正则化技术，特别是压缩感知（Compressive Sensing，CS）理论中提出的基于稀疏的全变分（Total Variation，TV）正则化方法是处理这种不适定问题的有效方式。2009年，Brady等 ［ 15 15 ］ 首次将压缩感知理论应用于全息成像中，实现了单幅全息图的深度重建和孪生像的去除。Wu等 ［ 16 16 ］ 根据孪生像和原始图像在梯度域稀疏程度的不同提出一种基于CS的无透镜成像方法，它通过在目标函数中加入TV正则化约束以改善问题的不适定性。在图像处理和计算成像技术中，通过外加约束先验的稀疏性在逆问题中可以增强图像统计的先验信息。虽然正则化先验技术常被用于不适定问题的求解。然而这类人为选择的正则化器存在校准误差和手动调整参数等不定因素，并不能表示真正的数据。 近年来，基于分数匹配生成模型 ［ 17 17 - 19 19 17-19 ］ 和去噪扩散概率模型（DDPMs） ［ 20 20 - 21 21 20-21 ］ 的新的生成模型 ［ 17 17 - 21 21 17-21 ］ 引起了人们的广泛兴趣，它在不需要对抗训练的情况下能够实现高质量的样本。生成模型本质是用一个已知的概率模型来拟合所给的数据样本。虽然神经网络（Neural Networks）也能拟合任意函数，但是概率分布有“非负”和“归一化”的要求，因此不能随意拟合一个概率分布。近来，Song等人 ［ 22 22 ］ 提出了一个新的基于广义离散分数匹配方法模型。该模型通过概率密度函数对观察到的样本进行数据建模，在利用已知输入不变性建模数据分布方面表现出了出色的性能。生成模型在拟合参数学习的过程会尽可能地增加数据的相似度，因而可以学习良好的先验信息。近年来，有不少该领域的研究工作聚焦于如何更好地利用算法提升来改善成像质量，尤其是针对深度重构的真实度。基于此，本文利用生成模型强大的先验信息建模能力和成像的拟合约束相结合的方法，将生成模型引入到无透镜成像逆问题求解。本文提出一种基于分数匹配生成模型的无透镜方法（Lens-less imaging via Score-based Generative Model，LSGM）。具体而言，无透镜成像系统在非相干照明下使用单张FZA实现点源编码，无需精确校准和严格标定，具有更高的信噪比。然而，由于孪生像的伪影很容易导致重建质量差以及颜色重建不准确等问题。针对这个问题，本文提出一种新的LSGM算法对编码图像进行解码重建。该模型通过利用梯度数据分布对图像样本进行建模，并从给定的观测数据的条件分布中采样数据，然后将采样的数据作为先验信息和成像的拟合约束相结合来求解无透镜成像的逆问题。该方法模型训练稳定、结构灵活，在数据分布建模方面表现出了强大的性能，并且能获得丰富的先验信息，可以有效地用于无透镜成像中逆问题的求解，并有望扩展至其他计算成像领域。 2 基本原理及算法过程 1 2.1　无透镜成像原理 2 如 图1 图1 所示，无透镜成像的主要流程主要包含两大部分：光学编码和计算解码。不同于传统光学成像——“所见即所得”。无透镜成像直接应用图像传感器记录由编码掩膜调制的场景信息，然后通过算法重建场景图像。 若将FZA的透射率函数近似为理想波带片，其振幅传递函数为： <math display="block" id="M1"><mi>t</mi><mo stretchy="false">(</mo><mi>r</mi><mo stretchy="false">)</mo><mo>=</mo><mfenced open="{" close="" separators="|"><mrow><mtable><mtr><mtd><mn mathvariant="normal">1</mn><mo>,</mo></mtd><mtd><mi mathvariant="normal">s</mi><mi mathvariant="normal">i</mi><mi mathvariant="normal">n</mi><mfenced separators="|"><mrow><mfrac><mrow><mi mathvariant="normal">π</mi><msup><mrow><mi>r</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup></mrow><mrow><msubsup><mrow><mi>r</mi></mrow><mrow><mn mathvariant="normal">1</mn></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msubsup></mrow></mfrac></mrow></mfenced><mo>≥</mo><mn mathvariant="normal">0</mn></mtd></mtr><mtr><mtd><mn mathvariant="normal">0</mn><mo>,</mo></mtd><mtd><mi mathvariant="normal">s</mi><mi mathvariant="normal">i</mi><mi mathvariant="normal">n</mi><mfenced separators="|"><mrow><mfrac><mrow><mi mathvariant="normal">π</mi><msup><mrow><mi>r</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup></mrow><mrow><msubsup><mrow><mi>r</mi></mrow><mrow><mn mathvariant="normal">1</mn></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msubsup></mrow></mfrac></mrow></mfenced><mo><</mo><mn mathvariant="normal">0</mn></mtd></mtr></mtable></mrow></mfenced></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34991855&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34991852&type= 40.80933380 22.18266678 . 对 式（1） 式（1） 作傅里叶展开为 <math id="M2"><mi>T</mi><mo stretchy="false">(</mo><mi>r</mi><mo stretchy="false">)</mo><mo>=</mo><mfenced open="[" close="]" separators="|"><mrow><mn mathvariant="normal">1</mn><mo>+</mo><mi mathvariant="normal">c</mi><mi mathvariant="normal">o</mi><mi mathvariant="normal">s</mi><mo stretchy="false">(</mo><mi mathvariant="normal">π</mi><msup><mrow><mi>r</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup><mo>/</mo><msubsup><mrow><mi>r</mi></mrow><mrow><mn mathvariant="normal">1</mn></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msubsup><mo stretchy="false">)</mo></mrow></mfenced><mo>/</mo><mn mathvariant="normal">2</mn></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34991860&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34991858&type= 45.04266739 4.57200003 。其中： <math id="M3"><mi>r</mi><mo>=</mo><mo stretchy="false">(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo stretchy="false">)</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34991867&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34991864&type= 15.15533447 3.97933316 ； <math id="M4"><msub><mrow><mi>r</mi></mrow><mrow><mn mathvariant="normal">1</mn></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34991874&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34991871&type= 2.62466669 3.72533321 为常数，表示FZA最内侧区域的半径。然后根据欧拉公式转化为： <math display="block" id="M5"><mi>T</mi><mo stretchy="false">(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo stretchy="false">)</mo><mo>=</mo><mfrac><mrow><mn mathvariant="normal">1</mn></mrow><mrow><mn mathvariant="normal">2</mn></mrow></mfrac><mo>+</mo><mfrac><mrow><mn mathvariant="normal">1</mn></mrow><mrow><mn mathvariant="normal">4</mn></mrow></mfrac><mfenced open="[" close="]" separators="|"><mrow><mi>h</mi><mo stretchy="false">(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo stretchy="false">)</mo><mo>+</mo><msup><mrow><mi>h</mi></mrow><mrow><mi mathvariant="normal">*</mi></mrow></msup><mo stretchy="false">(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo stretchy="false">)</mo></mrow></mfenced></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34991880&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34991877&type= 59.60533142 7.87400007 ， 其中： <math id="M6"><mi>h</mi><mo stretchy="false">(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo stretchy="false">)</mo><mo>=</mo><mi mathvariant="normal">e</mi><mi mathvariant="normal">x</mi><mi mathvariant="normal">p</mi><mo stretchy="false">(</mo><mi>i</mi><mo stretchy="false">(</mo><mi>π</mi><mo>/</mo><msubsup><mrow><mi>r</mi></mrow><mrow><mn mathvariant="normal">1</mn></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msubsup><mo stretchy="false">)</mo><mo stretchy="false">(</mo><msup><mrow><mi>x</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup><mo>+</mo><msup><mrow><mi>y</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup><mo stretchy="false">)</mo><mo stretchy="false">)</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34991885&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34991883&type= 51.56200027 5.16466665 ， <math id="M7"><msup><mrow><mi>h</mi></mrow><mrow><mi mathvariant="normal">*</mi></mrow></msup><mo stretchy="false">(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo stretchy="false">)</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34991891&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34991888&type= 12.78466606 4.23333359 是 <math id="M8"><mi>h</mi><mo stretchy="false">(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo stretchy="false">)</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34991904&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34991901&type= 11.51466751 3.97933316 的共轭。当 <math id="M9"><msubsup><mrow><mi>r</mi></mrow><mrow><mn mathvariant="normal">1</mn></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msubsup><mo>=</mo><mi>λ</mi><mi>d</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34991898&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34991896&type= 11.76866722 3.97933316 时， <math id="M10"><mi>h</mi><mo stretchy="false">(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo stretchy="false">)</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34991904&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34991901&type= 11.51466751 3.97933316 与菲涅尔传播核函数 <math id="M11"><mi mathvariant="normal">e</mi><mi mathvariant="normal">x</mi><mi mathvariant="normal">p</mi><mo stretchy="false">(</mo><mi>i</mi><mo stretchy="false">(</mo><mi mathvariant="normal">π</mi><mo>/</mo><mi>λ</mi><mi>d</mi><mo stretchy="false">)</mo><mo stretchy="false">(</mo><msup><mrow><mi>x</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup><mo>+</mo><msup><mrow><mi>y</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup><mo stretchy="false">)</mo><mo stretchy="false">)</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34991907&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34991905&type= 36.57600021 4.23333359 相同。 目标物体被放置在离FZA一定距离处，在非相干光源照明下，经过FZA调制的所有点源叠加的集合被传感器捕获。其成像过程的表达式即为目标图像与FZA函数的卷积： <math display="block" id="M12"><mtable><mtr><mtd><mi>I</mi><mo stretchy="false">(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo stretchy="false">)</mo><mo>=</mo><mi>O</mi><mo stretchy="false">(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo stretchy="false">)</mo><mi mathvariant="normal">*</mi><mi>T</mi><mo stretchy="false">(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo stretchy="false">)</mo><mo>=</mo><mfrac><mrow><mn mathvariant="normal">1</mn></mrow><mrow><mn mathvariant="normal">2</mn></mrow></mfrac><mi>O</mi><mo stretchy="false">(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo stretchy="false">)</mo><mo>+</mo><mfrac><mrow><mn mathvariant="normal">1</mn></mrow><mrow><mn mathvariant="normal">4</mn></mrow></mfrac><mi>O</mi><mo stretchy="false">(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo stretchy="false">)</mo><mi mathvariant="normal">*</mi><mfenced open="[" close="]" separators="|"><mrow><mi>h</mi><mo stretchy="false">(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo stretchy="false">)</mo><mo>+</mo><msup><mrow><mi>h</mi></mrow><mrow><mi mathvariant="normal">*</mi></mrow></msup><mo stretchy="false">(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo stretchy="false">)</mo></mrow></mfenced><mo>=</mo></mtd></mtr><mtr><mtd><mi>C</mi><mo>+</mo><mfrac><mrow><mn mathvariant="normal">1</mn></mrow><mrow><mn mathvariant="normal">4</mn></mrow></mfrac><mfenced open="[" close="]" separators="|"><mrow><mi>U</mi><mo stretchy="false">(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo stretchy="false">)</mo><mo>+</mo><msup><mrow><mi>U</mi></mrow><mrow><mi mathvariant="normal">*</mi></mrow></msup><mo stretchy="false">(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo stretchy="false">)</mo></mrow></mfenced><mo>=</mo><mi>C</mi><mo>+</mo><mi mathvariant="normal">R</mi><mi mathvariant="normal">e</mi><mfenced open="{" close="}" separators="|"><mrow><mi>U</mi><mo stretchy="false">(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo stretchy="false">)</mo></mrow></mfenced><mo>.</mo></mtd></mtr></mtable></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34991912&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34991909&type= 122.08934021 17.52599907 其中：*表示卷积； <math id="M13"><mi>O</mi><mo stretchy="false">(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo stretchy="false">)</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34991919&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34991915&type= 12.61533356 3.97933316 为目标图像； <math id="M14"><mi>U</mi><mo stretchy="false">(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo stretchy="false">)</mo><mo>=</mo><mi>O</mi><mo stretchy="false">(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo stretchy="false">)</mo><mi mathvariant="normal">*</mi><mi>h</mi><mo stretchy="false">(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo stretchy="false">)</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34991924&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34991920&type= 43.68799973 3.97933316 可以看作是波长为 <math id="M15"><mi>λ</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34991928&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34991926&type= 1.77800000 2.79399991 的光波在距离 <math id="M16"><mi>d</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34991933&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34991930&type= 2.28600001 2.79399991 处传播的衍射波前图像； <math id="M17"><msup><mrow><mi>U</mi></mrow><mrow><mi mathvariant="normal">*</mi></mrow></msup><mo stretchy="false">(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo stretchy="false">)</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34991944&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34991942&type= 14.22399998 4.23333359 与 <math id="M18"><mi>U</mi><mo stretchy="false">(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo stretchy="false">)</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34992028&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34992027&type= 12.86933327 3.97933316 相互共轭， 式（3） 式（3） 中 <math id="M19"><msup><mrow><mi>U</mi></mrow><mrow><mi mathvariant="normal">*</mi></mrow></msup><mo stretchy="false">(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo stretchy="false">)</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34991944&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34991942&type= 14.22399998 4.23333359 相关的项即为孪生像； <math id="M20"><mi>C</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34991948&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34991946&type= 2.62466669 2.79399991 为常数，可直接通过直流滤波去除。将上式卷积计算转换到傅里叶频率域，计算公式表示为： <math display="block" id="M21"><mi>I</mi><mo>=</mo><mfrac><mrow><mn mathvariant="normal">1</mn></mrow><mrow><mn mathvariant="normal">2</mn></mrow></mfrac><mi mathvariant="normal">R</mi><mi mathvariant="normal">e</mi><mfenced open="{" close="}" separators="|"><mrow><msup><mrow><mi>F</mi></mrow><mrow><mo>-</mo><mn mathvariant="normal">1</mn></mrow></msup><mi>H</mi><mi>F</mi><mi>O</mi></mrow></mfenced></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34991951&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34991949&type= 31.91933441 7.87400007 ， 其中， <math id="M22"><mi>I</mi><mo>∈</mo><msup><mrow><mi mathvariant="bold">R</mi></mrow><mrow><msub><mrow><mi>N</mi></mrow><mrow><mi>x</mi><mi>y</mi></mrow></msub></mrow></msup></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34991955&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34991953&type= 10.75266743 3.55599999 为传感器测的图像数据； <math id="M23"><mi>O</mi><mo>∈</mo><msup><mrow><mi mathvariant="bold">R</mi></mrow><mrow><msub><mrow><mi>N</mi></mrow><mrow><mi>x</mi><mi>y</mi></mrow></msub></mrow></msup></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34991959&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34991957&type= 12.36133385 3.55599999 为目标图像； <math id="M24"><mi>F</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34991963&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34991961&type= 2.62466669 2.79399991 和 <math id="M25"><msup><mrow><mi>F</mi></mrow><mrow><mo>-</mo><mn mathvariant="normal">1</mn></mrow></msup></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34991967&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34991965&type= 5.58799982 3.04800010 分别为傅里叶变换算子和逆傅里叶变换算子； <math id="M26"><mi>H</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34991971&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34991969&type= 3.30200005 2.79399991 表示乘以传递函数 <math id="M27"><mi>H</mi><mo>=</mo><mi>i</mi><mi mathvariant="normal">e</mi><mi mathvariant="normal">x</mi><mi mathvariant="normal">p</mi><mo stretchy="false">(</mo><mo>-</mo><mi>i</mi><mi>π</mi><mi>λ</mi><mi>z</mi><mo stretchy="false">(</mo><msup><mrow><mi>u</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup><mo>+</mo><msup><mrow><mi>v</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup><mo stretchy="false">)</mo><mo stretchy="false">)</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34991975&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34991974&type= 42.50266647 4.14866638 的操作算子。由于 <math id="M28"><mi>O</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34992081&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34992079&type= 2.96333337 2.87866688 为实函数、 <math id="M29"><mi>H</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34991984&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34991982&type= 3.30200005 2.79399991 为中心对称函数，则 <math id="M30"><mi mathvariant="normal">R</mi><mi mathvariant="normal">e</mi><mfenced open="{" close="}" separators="|"><mrow><mi>H</mi></mrow></mfenced><mo>=</mo><mi mathvariant="normal">s</mi><mi mathvariant="normal">i</mi><mi mathvariant="normal">n</mi><mo stretchy="false">(</mo><mi>π</mi><msubsup><mrow><mi>r</mi></mrow><mrow><mn mathvariant="normal">1</mn></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msubsup><mo stretchy="false">(</mo><msup><mrow><mi>u</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup><mo>+</mo><msup><mrow><mi>v</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup><mo stretchy="false">)</mo><mo stretchy="false">)</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34991988&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34991986&type= 44.02666473 3.97933316 是 <math id="M31"><mi mathvariant="normal">c</mi><mi mathvariant="normal">o</mi><mi mathvariant="normal">s</mi><mo stretchy="false">(</mo><mo stretchy="false">(</mo><mi>π</mi><mo>/</mo><mi>λ</mi><mi>d</mi><mo stretchy="false">)</mo><mo stretchy="false">(</mo><msup><mrow><mi>x</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup><mo>+</mo><msup><mrow><mi>y</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup><mo stretchy="false">)</mo><mo stretchy="false">)</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34991992&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34991990&type= 34.29000092 4.23333359 的归一化傅里叶变换。假设系统转换矩阵 <math id="M32"><mi mathvariant="bold-italic">Η</mi><mo>=</mo><msup><mrow><mi mathvariant="bold-italic">F</mi></mrow><mrow><mi mathvariant="normal">*</mi></mrow></msup><mi mathvariant="bold-italic">Λ</mi><mi mathvariant="bold-italic">F</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34991995&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34991994&type= 18.11866570 3.04800010 ，其中 <math id="M33"><mi mathvariant="bold-italic">F</mi><mo>∈</mo><msup><mrow><mi mathvariant="bold">R</mi></mrow><mrow><msub><mrow><mi>N</mi></mrow><mrow><mi>x</mi><mi>y</mi></mrow></msub><mo>×</mo><msub><mrow><mi>N</mi></mrow><mrow><mi>x</mi><mi>y</mi></mrow></msub></mrow></msup></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34991999&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34991997&type= 17.44133186 3.55599999 为二维离散傅里叶变换矩阵， <math id="M34"><mi mathvariant="bold-italic">Λ</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34992003&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34992001&type= 2.96333337 2.79399991 为对角矩阵，其非零项为 <math id="M35"><mi mathvariant="normal">R</mi><mi mathvariant="normal">e</mi><mfenced open="{" close="}" separators="|"><mrow><mi>H</mi></mrow></mfenced></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34992008&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34992006&type= 11.43000031 3.97933316 的离散值。则用 <math id="M36"><mi mathvariant="bold-italic">O</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34992012&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34992010&type= 2.96333337 2.79399991 表示 <math id="M37"><mi mathvariant="bold-italic">I</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34992016&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34992014&type= 1.86266661 2.79399991 的相关函数即为正向转换模型： <math display="block" id="M38"><mi mathvariant="bold-italic">I</mi><mo>=</mo><mi mathvariant="bold-italic">H</mi><mi mathvariant="bold-italic">O</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34992020&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34992018&type= 12.36133385 2.79399991 ， 其中， <math id="M39"><mi mathvariant="bold-italic">H</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34992025&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34992022&type= 3.47133350 2.79399991 为系统测量矩阵或者传递函数矩阵。上述正向模型是一个典型的逆问题。在 式（3） 式（3） 中， <math id="M40"><mi>U</mi><mo stretchy="false">(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo stretchy="false">)</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34992028&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34992027&type= 12.86933327 3.97933316 的虚部可以被任意赋值，使得求解不唯一，从而产生孪生像。孪生像的伪影很容易导致重建质量差以及颜色重建不准确等问题，引入合适的图像先验是解决该问题的关键。 传统的方法包括变换域中图像的稀疏性、图像梯度的稀疏性、正则化先验、低秩等。随着生成模型的研究发展，图像的生成式建模的进展已经远远超出传统的基于先验的算法。新的一类生成模型 ［ 22 22 - 24 24 22-24 ］ 被提出，它们提供了一种使用梯度数据分布对图像进行建模的强大方法，利用学习的分数函数作为先验，可以很容易地解决成像中的逆问题。近来，Song等人 ［ 22 22 ］ 提出使用随机微分方程（Stochastic Differential Equation，SDE）缓慢地添加噪声将复杂数据分布转换为已知的先验分布，并使用相应的反向SDE缓慢地去除噪声将先验分布转换回数据分布。至关重要的是，反向SDE仅取决于扰动数据分布的时间相关梯度场 ［ 25 25 ］ ，也称为分数函数。为方便表述，本文将基于分数的SDE称为基于分数匹配的生成模型。 2.2　基于分数匹配生成模型的无透镜成像方法 2 如 图2 图2 所示，生成模型可以简单地概括为用概率分布方式描述图像的生成，通过对概率分布采样产生数据。 图2 图2 右边是一个训练数据集，里面所有的数据都是从某个数据 <math id="M41"><msub><mrow><mi>p</mi></mrow><mrow><mi mathvariant="normal">d</mi><mi mathvariant="normal">a</mi><mi mathvariant="normal">t</mi><mi mathvariant="normal">a</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34992033&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34992031&type= 5.75733376 4.57200003 中独立同分布取出的随机样本。左边是其生成模型（即概率分布），在这种概率分布中，找出一个分布 <math id="M42"><msub><mrow><mi>p</mi></mrow><mrow><mi>θ</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34992044&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34992042&type= 3.38666677 4.57200003 使得它离 <math id="M43"><msub><mrow><mi>p</mi></mrow><mrow><mi mathvariant="normal">d</mi><mi mathvariant="normal">a</mi><mi mathvariant="normal">t</mi><mi mathvariant="normal">a</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34992040&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34992038&type= 5.75733376 4.57200003 的距离最近。接着在 <math id="M44"><msub><mrow><mi>p</mi></mrow><mrow><mi>θ</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34992044&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34992042&type= 3.38666677 4.57200003 上采新的样本，可以获得源源不断的新数据。通常，使用分数函数（概率密度函数对数的导数）来表示概率分布，称为分数匹配的生成模型。相比概率密度函数模型，它没有归一化的要求，因而模型选择更加灵活。 基于分数匹配的生成模型是一种通过优化以 <math id="M45"><mi>θ</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34992178&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34992177&type= 1.86266661 2.79399991 参数化的分数网络 <math id="M46"><msub><mrow><mi>S</mi></mrow><mrow><mi>θ</mi></mrow></msub><mo stretchy="false">(</mo><mi>x</mi><mo stretchy="false">)</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34992058&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34992056&type= 9.39799976 4.48733330 来估计数据分布的分数函数的方法。用多个噪声尺度对数据进行扰动是分数匹配生成模型方法成功的关键。它定义了一个正向扩散过程（即SDE），用于将数据转换为噪声并通过其逆向过程从噪声中生成数据。如 图3 图3 所示，分数模型使用一个正向SDE通过添加高斯噪声扰动数据分布，为了从数据分布中进行采样，可以训练一个神经网络来估计对数数据分布的梯度（即 <math id="M47"><msub><mrow><mo>∇</mo></mrow><mrow><mi>x</mi></mrow></msub><mi mathvariant="normal">l</mi><mi mathvariant="normal">o</mi><mi mathvariant="normal">g</mi><msub><mrow><mi>p</mi></mrow><mrow><mi>t</mi></mrow></msub><mo stretchy="false">(</mo><mi>x</mi><mo stretchy="false">)</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34992062&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34992060&type= 18.96533203 4.82600021 ），并用它来数值求解逆时间SDE从噪声中生成数据。 本文通过使用分数模型对目标图像数据集进行训练，将学习到的图像先验信息引入无透镜成像逆问题，然后通过预测-校正采样算法以实现高质量的无透镜成像重建。 具体来说，训练阶段，通过将目标图像数据集输入基于分数匹配生成模型进行训练。该模型使用SDE通过缓慢添加高斯噪声扰动数据分布，然后使用方差爆炸（Variance Exploding，VE）-SDE的去噪分数匹配方法训练了一个连续时间相关的分数函数 <math id="M48"><msub><mrow><mi>S</mi></mrow><mrow><mi>θ</mi></mrow></msub><mo stretchy="false">(</mo><mi>O</mi><mo stretchy="false">)</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34992073&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34992070&type= 10.07533360 4.48733330 估计对数数据分布的梯度，即模型学习的梯度先验 <math id="M49"><msub><mrow><mo>∇</mo></mrow><mrow><mi>O</mi></mrow></msub><mi mathvariant="normal">l</mi><mi mathvariant="normal">o</mi><mi mathvariant="normal">g</mi><msub><mrow><mi>p</mi></mrow><mrow><mi>t</mi></mrow></msub><mo stretchy="false">(</mo><mi>O</mi><mo stretchy="false">)</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34992077&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34992075&type= 20.06599998 4.82600021 ，其中 <math id="M50"><mi>O</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34992081&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34992079&type= 2.96333337 2.87866688 为目标图像数据集。测试阶段，根据训练的分数函数（即目标图像先验）和逆时间SDE数值求解器迭代更新重建目标图像。在非相干光照明下由单张FZA调制的编码图像，通过训练好的分数网络先经过数值求解器SDE一次预测，再通过退火朗之万方法在数值求解器SDE和数据保真项步骤之间轮换更新以实现编码图像的重建。基于分数匹配生成模型的无透镜成像方法的训练和重建算法将分别在2.3和2.4节展开介绍。 2.3　基于分数匹配生成模型的训练 2 图4 图4 为LSGM算法的训练与重建流程图。本论文提出的LSGM方法的训练过程如 图4 图4 上半部分所示。首先，模型构造了一个连续扩散过程 <math id="M51"><msubsup><mrow><mfenced open="{" close="}" separators="|"><mrow><mi>O</mi><mo stretchy="false">(</mo><mi>t</mi><mo stretchy="false">)</mo></mrow></mfenced></mrow><mrow><mi>t</mi><mo>=</mo><mn mathvariant="normal">0</mn></mrow><mrow><mi mathvariant="normal">T</mi></mrow></msubsup></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34992085&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34992083&type= 15.57866764 5.84200001 ，其中， <math id="M52"><mi>O</mi><mo stretchy="false">(</mo><mi>t</mi><mo stretchy="false">)</mo><mo>∈</mo><msup><mrow><mi mathvariant="bold">R</mi></mrow><mrow><msub><mrow><mi>N</mi></mrow><mrow><mi>x</mi><mi>y</mi></mrow></msub></mrow></msup></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34992089&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34992086&type= 17.10266685 3.64066648 ， <math id="M53"><mi>t</mi><mo>∈</mo><mo stretchy="false">[</mo><mn mathvariant="normal">0</mn><mo>,</mo><mi>T</mi><mo stretchy="false">]</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34992093&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34992091&type= 14.98600006 3.55599999 是时间变量。正向过程是将数据分布 <math id="M54"><mi>O</mi><mo stretchy="false">(</mo><mn mathvariant="normal">0</mn><mo stretchy="false">)</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34992190&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34992189&type= 8.38199997 3.55599999 通过随机过程扰动得到含噪图像 <math id="M55"><mi>O</mi><mo stretchy="false">(</mo><mi>T</mi><mo stretchy="false">)</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34992100&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34992098&type= 9.31333351 3.55599999 。其中， <math id="M56"><mi>O</mi><mo stretchy="false">(</mo><mn mathvariant="normal">0</mn><mo stretchy="false">)</mo><mo>∼</mo><msub><mrow><mi>p</mi></mrow><mrow><mi mathvariant="normal">d</mi><mi mathvariant="normal">a</mi><mi mathvariant="normal">t</mi><mi mathvariant="normal">a</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34992146&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34992144&type= 17.78000069 4.82600021 ， <math id="M57"><msub><mrow><mi>p</mi></mrow><mrow><mi mathvariant="normal">d</mi><mi mathvariant="normal">a</mi><mi mathvariant="normal">t</mi><mi mathvariant="normal">a</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34992152&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34992151&type= 5.75733376 4.57200003 表示目标图像数据分布； <math id="M58"><mi>O</mi><mo stretchy="false">(</mo><mi>T</mi><mo stretchy="false">)</mo><mo>∼</mo><msub><mrow><mi>p</mi></mrow><mrow><mi>T</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34992142&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34992140&type= 16.76399994 4.82600021 ， <math id="M59"><msub><mrow><mi>p</mi></mrow><mrow><mi>T</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34992150&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34992148&type= 3.80999994 4.57200003 表示先验数据分布。那么，该随机过程可以建模为求解以下SDE过程： <math display="block" id="M60"><mi mathvariant="normal">d</mi><mi>O</mi><mo>=</mo><mi>f</mi><mo stretchy="false">(</mo><mi>O</mi><mo>,</mo><mi>t</mi><mo stretchy="false">)</mo><mi mathvariant="normal">d</mi><mi>t</mi><mo>+</mo><mi>g</mi><mo stretchy="false">(</mo><mi>t</mi><mo stretchy="false">)</mo><mi mathvariant="normal">d</mi><mi>w</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34992120&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34992118&type= 41.65599823 3.97933316 ， 其中： <math id="M61"><mi>f</mi><mo>∈</mo><msup><mrow><mi mathvariant="bold">R</mi></mrow><mrow><msub><mrow><mi>N</mi></mrow><mrow><mi>x</mi><mi>y</mi></mrow></msub></mrow></msup></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34992124&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34992122&type= 11.26066685 4.06400013 和 <math id="M62"><mi>g</mi><mo>∈</mo><mi mathvariant="bold">R</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34992128&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34992126&type= 8.46666718 3.64066648 分别是漂移系数和扩散系数， <math id="M63"><mi>w</mi><mo>∈</mo><msup><mrow><mi mathvariant="bold">R</mi></mrow><mrow><msub><mrow><mi>N</mi></mrow><mrow><mi>x</mi><mi>y</mi></mrow></msub></mrow></msup></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34992132&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34992130&type= 12.36133385 3.55599999 是标准维纳过程或称布朗运动。 然后，通过求解逆时间SDE来重建编码图像。具体而言，从正向SDE训练过程学习的目标相关的先验信息 <math id="M64"><mi>O</mi><mo stretchy="false">(</mo><mi>T</mi><mo stretchy="false">)</mo><mo>∼</mo><msub><mrow><mi>p</mi></mrow><mrow><mi>T</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34992142&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34992140&type= 16.76399994 4.82600021 开始求解 式（6） 式（6） 的逆向问题可以得到 <math id="M65"><mi>O</mi><mo stretchy="false">(</mo><mn mathvariant="normal">0</mn><mo stretchy="false">)</mo><mo>∼</mo><msub><mrow><mi>p</mi></mrow><mrow><mi mathvariant="normal">d</mi><mi mathvariant="normal">a</mi><mi mathvariant="normal">t</mi><mi mathvariant="normal">a</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34992146&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34992144&type= 17.78000069 4.82600021 目标图像，其中 <math id="M66"><msub><mrow><mi>p</mi></mrow><mrow><mi>T</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34992150&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34992148&type= 3.80999994 4.57200003 是一个包含 <math id="M67"><msub><mrow><mi>p</mi></mrow><mrow><mi mathvariant="normal">d</mi><mi mathvariant="normal">a</mi><mi mathvariant="normal">t</mi><mi mathvariant="normal">a</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34992152&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34992151&type= 5.75733376 4.57200003 信息的非结构化先验分布。那么，逆SDE过程为： <math display="block" id="M68"><mtable columnalign="left"><mtr><mtd><mi mathvariant="normal">d</mi><mi>O</mi><mo>=</mo><mfenced open="[" close="]" separators="|"><mrow><mi>f</mi><mo stretchy="false">(</mo><mi>O</mi><mo>,</mo><mi>t</mi><mo stretchy="false">)</mo><mo>-</mo><mi>g</mi><msup><mrow><mo stretchy="false">(</mo><mi>t</mi><mo stretchy="false">)</mo></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup><msub><mrow><mo>∇</mo></mrow><mrow><mi>O</mi></mrow></msub><mi mathvariant="normal">l</mi><mi mathvariant="normal">o</mi><mi mathvariant="normal">g</mi><msub><mrow><mi>p</mi></mrow><mrow><mi>t</mi></mrow></msub><mo stretchy="false">(</mo><mi>O</mi><mo stretchy="false">)</mo></mrow></mfenced><mi mathvariant="normal">d</mi><mi>t</mi><mo>+</mo></mtd></mtr><mtr><mtd><mi>g</mi><mo stretchy="false">(</mo><mi>t</mi><mo stretchy="false">)</mo><mi mathvariant="normal">d</mi><mover accent="true"><mi>w</mi><mo>¯</mo></mover><mo>,</mo></mtd></mtr></mtable></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34992154&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34992153&type= 65.10866547 10.58333302 其中： <math id="M69"><mi mathvariant="normal">d</mi><mi>t</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34992156&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34992155&type= 3.13266683 2.79399991 是无穷小的负时间步长， <math id="M70"><mover accent="true"><mi>w</mi><mo>¯</mo></mover></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34992158&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34992157&type= 2.96333337 2.79399991 是反向的布朗运动。本文为了产生更高质量的样本，通过选取函数 <math id="M71"><mi>f</mi><mo>=</mo><mn mathvariant="normal">0</mn></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34992161&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34992159&type= 8.46666718 3.64066648 和 <math id="M72"><mi>g</mi><mo>=</mo><mroot><mrow><mi mathvariant="normal">d</mi><mo stretchy="false">[</mo><msup><mrow><mi>σ</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup><mo stretchy="false">(</mo><mi>t</mi><mo stretchy="false">)</mo><mo stretchy="false">]</mo><mo>/</mo><mi mathvariant="normal">d</mi><mi>t</mi></mrow><mrow/></mroot></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34992163&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34992162&type= 30.31066895 5.07999992 构造了一个VE-SDE ［ 22 22 ］ 。 为了求解 公式（7） 公式（7） ，需要已知所有 <math id="M73"><mi>t</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34992200&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34992199&type= 1.26999998 2.79399991 的分数函数 <math id="M74"><msub><mrow><mo>∇</mo></mrow><mrow><mi>O</mi></mrow></msub><mi mathvariant="normal">l</mi><mi mathvariant="normal">o</mi><mi mathvariant="normal">g</mi><msub><mrow><mi>p</mi></mrow><mrow><mi>t</mi></mrow></msub><mo stretchy="false">(</mo><mi>O</mi><mo stretchy="false">(</mo><mi>t</mi><mo stretchy="false">)</mo><mo stretchy="false">)</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34992168&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34992166&type= 24.97666740 4.82600021 。通过训练一个基于时间的评分网络 <math id="M75"><msub><mrow><mi>S</mi></mrow><mrow><mi>θ</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34992172&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34992169&type= 3.30200005 3.72533321 来匹配 <math id="M76"><msub><mrow><mi>S</mi></mrow><mrow><mi>θ</mi></mrow></msub><mo stretchy="false">(</mo><mi>O</mi><mo>,</mo><mi>t</mi><mo stretchy="false">)</mo><mo>≃</mo><msub><mrow><mo>∇</mo></mrow><mrow><mi mathvariant="normal">O</mi></mrow></msub><mi mathvariant="normal">l</mi><mi mathvariant="normal">o</mi><mi mathvariant="normal">g</mi><msub><mrow><mi>p</mi></mrow><mrow><mi>t</mi></mrow></msub><mo stretchy="false">(</mo><mi>O</mi><mo stretchy="false">(</mo><mi>t</mi><mo stretchy="false">)</mo><mo stretchy="false">)</mo><mtext> </mtext></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34992175&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34992173&type= 43.43399811 4.82600021 ［ 26 26 ］ ，然后将这个网络函数替换 式（7） 式（7） 中相应项。因此，通过对分数网络的参数 <math id="M77"><mi>θ</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34992178&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34992177&type= 1.86266661 2.79399991 进行优化： <math display="block" id="M78"><mtable columnalign="left"><mtr><mtd><msup><mrow><mi>θ</mi></mrow><mrow><mi mathvariant="normal">*</mi></mrow></msup><mo>=</mo><munder><mrow><mi mathvariant="normal">a</mi><mi mathvariant="normal">r</mi><mi mathvariant="normal">g</mi><mi mathvariant="normal">m</mi><mi mathvariant="normal">i</mi><mi mathvariant="normal">n</mi></mrow><mrow><mi>θ</mi></mrow></munder><msub><mrow><mi>Ε</mi></mrow><mrow><mi>t</mi><mo>∼</mo><mi>U</mi><mo stretchy="false">(</mo><mn mathvariant="normal">0,1</mn><mo stretchy="false">)</mo></mrow></msub><mfenced open="{" close="" separators="|"><mrow><mi>λ</mi><mo stretchy="false">(</mo><mi>t</mi><mo stretchy="false">)</mo><msub><mrow><mi>Ε</mi></mrow><mrow><mi>O</mi><mo stretchy="false">(</mo><mn mathvariant="normal">0</mn><mo stretchy="false">)</mo></mrow></msub><msub><mrow><mi>Ε</mi></mrow><mrow><mi>O</mi><mo stretchy="false">(</mo><mi>t</mi><mo stretchy="false">)</mo><mo stretchy="false">|</mo><mi>O</mi><mo stretchy="false">(</mo><mn mathvariant="normal">0</mn><mo stretchy="false">)</mo></mrow></msub></mrow></mfenced><mo>·</mo></mtd></mtr><mtr><mtd><mfenced open="" close="}" separators="|"><mrow><mfenced open="[" close="]" separators="|"><mrow><msubsup><mrow><mfenced open="‖" close="‖" separators="|"><mrow><msub><mrow><mi>S</mi></mrow><mrow><mi>θ</mi></mrow></msub><mo stretchy="false">(</mo><mi>O</mi><mo stretchy="false">(</mo><mi>t</mi><mo stretchy="false">)</mo><mo>,</mo><mi>t</mi><mo stretchy="false">)</mo><mo>-</mo><msub><mrow><mo>∇</mo></mrow><mrow><mi>O</mi></mrow></msub><mi mathvariant="normal">l</mi><mi mathvariant="normal">o</mi><mi mathvariant="normal">g</mi><msub><mrow><mi>p</mi></mrow><mrow><mn mathvariant="normal">0</mn><mi>t</mi></mrow></msub><mo stretchy="false">(</mo><mi>O</mi><mo stretchy="false">(</mo><mi>t</mi><mo stretchy="false">)</mo><mo stretchy="false">|</mo><mi>O</mi><mo stretchy="false">(</mo><mn mathvariant="normal">0</mn><mo stretchy="false">)</mo><mo stretchy="false">)</mo></mrow></mfenced></mrow><mrow><mn mathvariant="normal">2</mn></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msubsup></mrow></mfenced></mrow></mfenced><mo>,</mo></mtd></mtr></mtable></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34992180&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34992179&type= 70.10399628 16.17133331 其中： E 表示期望， <math id="M79"><mi>λ</mi><mo>:</mo><mo stretchy="false">[</mo><mn mathvariant="normal">0</mn><mo>,</mo><mi>T</mi><mo stretchy="false">]</mo><mo>→</mo><msub><mrow><mi mathvariant="bold">R</mi></mrow><mrow><mo>></mo><mn mathvariant="normal">0</mn></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34992182&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34992181&type= 24.55333328 4.48733330 为正权重函数， <math id="M80"><mi>t</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34992200&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34992199&type= 1.26999998 2.79399991 在 <math id="M81"><mo stretchy="false">[</mo><mn mathvariant="normal">0</mn><mo>,</mo><mi>T</mi><mo stretchy="false">]</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34992186&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34992185&type= 10.15999985 3.55599999 上均匀采样， <math id="M82"><msub><mrow><mi>p</mi></mrow><mrow><mn mathvariant="normal">0</mn><mi>t</mi></mrow></msub><mo stretchy="false">(</mo><mi>O</mi><mo stretchy="false">(</mo><mi>t</mi><mo stretchy="false">)</mo><mo stretchy="false">|</mo><mi>O</mi><mo stretchy="false">(</mo><mn mathvariant="normal">0</mn><mo stretchy="false">)</mo><mo stretchy="false">)</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34992188&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34992187&type= 25.06133461 4.82600021 是以 <math id="M83"><mi>O</mi><mo stretchy="false">(</mo><mn mathvariant="normal">0</mn><mo stretchy="false">)</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34992190&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34992189&type= 8.38199997 3.55599999 为中心的高斯扰动核， <math id="M84"><msub><mrow><mo>∇</mo></mrow><mrow><mi>O</mi></mrow></msub><mi mathvariant="normal">l</mi><mi mathvariant="normal">o</mi><mi mathvariant="normal">g</mi><msub><mrow><mi>p</mi></mrow><mrow><mn mathvariant="normal">0</mn><mi>t</mi></mrow></msub><mo stretchy="false">(</mo><mi>O</mi><mo stretchy="false">(</mo><mi>t</mi><mo stretchy="false">)</mo><mo stretchy="false">|</mo><mi>O</mi><mo stretchy="false">(</mo><mn mathvariant="normal">0</mn><mo stretchy="false">)</mo><mo stretchy="false">)</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34992192&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34992191&type= 35.30599976 4.82600021 表示扰动核的梯度， <math id="M85"><mi>σ</mi><mo stretchy="false">(</mo><mi>t</mi><mo stretchy="false">)</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34992194&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34992193&type= 6.77333355 3.55599999 是噪声尺度的单调递增函数。通过 公式（8） 公式（8） 训练的网络模型应当满足 <math id="M86"><msub><mrow><mi>S</mi></mrow><mrow><mi>θ</mi></mrow></msub><mo stretchy="false">(</mo><mi>O</mi><mo>,</mo><mi>t</mi><mo stretchy="false">)</mo><mo>≃</mo><msub><mrow><mo>∇</mo></mrow><mrow><mi>O</mi></mrow></msub><mi mathvariant="normal">l</mi><mi mathvariant="normal">o</mi><mi mathvariant="normal">g</mi><msub><mrow><mi>p</mi></mrow><mrow><mi>t</mi></mrow></msub><mo stretchy="false">(</mo><mi>O</mi><mo stretchy="false">(</mo><mi>t</mi><mo stretchy="false">)</mo><mo stretchy="false">)</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34992196&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34992195&type= 42.67200089 4.82600021 ，这意味着 <math id="M87"><msub><mrow><mo>∇</mo></mrow><mrow><mi>O</mi></mrow></msub><mi mathvariant="normal">l</mi><mi mathvariant="normal">o</mi><mi mathvariant="normal">g</mi><msub><mrow><mi>p</mi></mrow><mrow><mi>t</mi></mrow></msub><mo stretchy="false">(</mo><mi>O</mi><mo stretchy="false">)</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34992198&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34992197&type= 20.06599998 4.82600021 对所有 <math id="M88"><mi>t</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34992200&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34992199&type= 1.26999998 2.79399991 都被训练的 <math id="M89"><msub><mrow><mi>S</mi></mrow><mrow><mi>θ</mi></mrow></msub><mo stretchy="false">(</mo><mi>O</mi><mo>,</mo><mi>t</mi><mo stretchy="false">)</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34992214&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34992213&type= 13.03866673 4.48733330 所知。因此，可以将其带入逆时间SDE求解 公式（7） 公式（7） 并对其进行模拟以进行采样。 2.4　基于分数匹配生成模型的重建 2 无透镜成像问题归结为求解方程（5），求解该逆问题的经典方法是转化为优化以下约束问题： <math display="block" id="M90"><mi mathvariant="normal">a</mi><mi mathvariant="normal">r</mi><mi mathvariant="normal">g</mi><munder><mrow><mi mathvariant="normal">m</mi><mi mathvariant="normal">i</mi><mi mathvariant="normal">n</mi></mrow><mrow><mi>O</mi></mrow></munder><msubsup><mrow><mfenced open="‖" close="‖" separators="|"><mrow><mi>I</mi><mo>-</mo><mi>H</mi><mi>O</mi></mrow></mfenced></mrow><mrow><mn mathvariant="normal">2</mn></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msubsup><mo>+</mo><mi>τ</mi><mi>Ψ</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34992204&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34992203&type= 39.20066452 6.01133299 ， 其中： <math id="M91"><mi>Ψ</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34992206&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34992205&type= 3.47133350 2.79399991 是正则化项， <math id="M92"><mi>τ</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34992212&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34992211&type= 1.69333339 2.79399991 是正则化系数。在求解上式时，通常使用变量分裂 ［ 27 27 ］ 或投影到凸集 ［ 28 28 ］ 等近似算法将先验信息项和数据保项优化解耦，然后根据轮换求解思路在两个子问题之间交替更新实现最优解： <math display="block" id="M93"><mfenced open="{" close="" separators="|"><mrow><mtable columnalign="left"><mtr><mtd><msub><mrow><mi>O</mi></mrow><mrow><mi>i</mi></mrow></msub><mo>=</mo><mi>G</mi><mo stretchy="false">(</mo><msub><mrow><mi>O</mi></mrow><mrow><mi>i</mi><mo>+</mo><mn mathvariant="normal">1</mn></mrow></msub><mo stretchy="false">)</mo></mtd></mtr><mtr><mtd><msub><mrow><mi>O</mi></mrow><mrow><mi>i</mi><mo>-</mo><mn mathvariant="normal">1</mn></mrow></msub><mo>=</mo><mi mathvariant="normal">a</mi><mi mathvariant="normal">r</mi><mi mathvariant="normal">g</mi><munder><mrow><mi mathvariant="normal">m</mi><mi mathvariant="normal">i</mi><mi mathvariant="normal">n</mi></mrow><mrow><mi>O</mi></mrow></munder><msubsup><mrow><mfenced open="‖" close="‖" separators="|"><mrow><mi>I</mi><mo>-</mo><mi>H</mi><mi>O</mi></mrow></mfenced></mrow><mrow><mn mathvariant="normal">2</mn></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msubsup><mo>+</mo><mi>τ</mi><mfenced open="‖" close="‖" separators="|"><mrow><mi>O</mi><mo>-</mo><msub><mrow><mi>O</mi></mrow><mrow><mi>i</mi></mrow></msub></mrow></mfenced></mtd></mtr></mtable></mrow></mfenced></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34992210&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34992209&type= 65.27799988 12.19200039 ， 其中，第一个公式为生成；第二个公式的第一项是数据保真项，第二项是正则化项， <math id="M94"><mi>τ</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34992212&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34992211&type= 1.69333339 2.79399991 是正则化系数。 根据贝叶斯最大后验估计，如果生成模型能够准确地从目标图像数据集中估计复杂先验数据分布，可以实现高质量的无透镜成像重建。在2.3节中，详细阐述了基于分数匹配生成模型的先验学习过程。在上述训练过程中，一旦训练好分数网络 <math id="M95"><msub><mrow><mi>S</mi></mrow><mrow><mi>θ</mi></mrow></msub><mo stretchy="false">(</mo><mi>O</mi><mo>,</mo><mi>t</mi><mo stretchy="false">)</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34992214&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34992213&type= 13.03866673 4.48733330 ，并将其带入 公式（7） 公式（7） 中，即可以求解逆时间VE-SDE，实现无透镜编码图像重建： <math display="block" id="M96"><mi mathvariant="normal">d</mi><mi>O</mi><mo>=</mo><mo>-</mo><mi mathvariant="normal">d</mi><msup><mrow><mi>σ</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup><mo stretchy="false">(</mo><mi>t</mi><mo stretchy="false">)</mo><mo>/</mo><msub><mrow><mi>S</mi></mrow><mrow><mi mathvariant="normal">θ</mi></mrow></msub><mo stretchy="false">(</mo><mi>O</mi><mo>,</mo><mi>t</mi><mo stretchy="false">)</mo><mo>+</mo><mroot><mrow><mi mathvariant="normal">d</mi><msup><mrow><mi>σ</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup><mo stretchy="false">(</mo><mi>t</mi><mo stretchy="false">)</mo><mo>/</mo><mi mathvariant="normal">d</mi><mi>t</mi></mrow><mrow/></mroot><mi mathvariant="normal">d</mi><mover accent="true"><mi>w</mi><mo>¯</mo></mover></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34992216&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34992215&type= 66.88666534 4.82600021 . 对于该过程的求解，可以使用欧拉离散化 ［ 22 22 ］ 求解VE-SDE的数值。为了校正离散化方向SDE演化中的错误，根据预测-校正的思想引入了一个预测-校正采样器（Predictor-Corrector，PC）。基于分数匹配生成模型先通过数值求解逆时间VE-SDE从得到的目标图像的先验分布中获得初步预测的重建图像，然后通过使用退火朗之万 ［ 29 29 ］ 作为校正器来校正这个初步预测结果。与此同时，在求解逆时间VE-SDE和退火朗之万采样时，每个无条件采样更新步骤之后都需要进行一次保真。基于无透镜成像的迭代重建算法如下： 预测：该预测器可以是逆时间VE-SDE的任意数值解算器，并具有固定的离散化策略。从目标图像先验分布中生成目标图像可以通过求解由 式（11） 式（11） 得到： <math display="block" id="M97"><msub><mrow><mi>O</mi></mrow><mrow><mi>i</mi></mrow></msub><mo>=</mo><mo stretchy="false">(</mo><msubsup><mrow><mi>σ</mi></mrow><mrow><mi>i</mi><mo>+</mo><mn mathvariant="normal">1</mn></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msubsup><mo>-</mo><msubsup><mrow><mi>σ</mi></mrow><mrow><mi>i</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msubsup><mo stretchy="false">)</mo><msub><mrow><mi>s</mi></mrow><mrow><mi>θ</mi></mrow></msub><mo stretchy="false">(</mo><msub><mrow><mi>O</mi></mrow><mrow><mi>i</mi><mo>+</mo><mn mathvariant="normal">1</mn></mrow></msub><mo>,</mo><msub><mrow><mi>σ</mi></mrow><mrow><mi>i</mi><mo>+</mo><mn mathvariant="normal">1</mn></mrow></msub><mo stretchy="false">)</mo><mo>+</mo><mroot><mrow><msubsup><mrow><mi>σ</mi></mrow><mrow><mi>i</mi><mo>+</mo><mn mathvariant="normal">1</mn></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msubsup><mo>-</mo><msubsup><mrow><mi>σ</mi></mrow><mrow><mi>i</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msubsup></mrow><mrow/></mroot><mi>z</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34992218&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34992217&type= 71.71266174 4.82600021 ， <math display="block" id="M98"><msub><mrow><mi>O</mi></mrow><mrow><mi>i</mi></mrow></msub><mo>=</mo><msup><mrow><mi>H</mi></mrow><mrow><mi mathvariant="normal">T</mi></mrow></msup><mo stretchy="false">(</mo><mi>I</mi><mo>-</mo><mi>Η</mi><msub><mrow><mi>O</mi></mrow><mrow><mi>i</mi></mrow></msub><mo stretchy="false">)</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34992220&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34992219&type= 29.88733101 4.65666676 ， 其中： <math id="M99"><msub><mrow><mi>σ</mi></mrow><mrow><mi>i</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34992222&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34992221&type= 2.45533323 3.72533321 噪声尺度， <math id="M100"><mi>i</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34992234&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34992233&type= 1.10066664 2.79399991 是总迭代次数， <math id="M101"><mi>z</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34992226&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34992225&type= 1.77800000 2.79399991 是零均值的高斯白噪声。 式（13） 式（13） 称为数据保真步骤，可以有效增强数据的一致性。 式（12） 式（12） 和（13）即为该算法的预测过程。 校正：本文使用朗之万马可夫链 ［ 29 29 ］ 校正算法校正梯度上升方向，得到校正器步骤： <math display="block" id="M102"><msub><mrow><mi>O</mi></mrow><mrow><mi>i</mi></mrow></msub><mo>=</mo><msub><mrow><mi>O</mi></mrow><mrow><mi>i</mi><mo>+</mo><mn mathvariant="normal">1</mn></mrow></msub><mo>+</mo><msub><mrow><mi>ε</mi></mrow><mrow><mi>i</mi></mrow></msub><msub><mrow><mi>S</mi></mrow><mrow><mi>θ</mi></mrow></msub><mo stretchy="false">(</mo><msub><mrow><mi>O</mi></mrow><mrow><mi>i</mi><mo>+</mo><mn mathvariant="normal">1</mn></mrow></msub><mo>,</mo><msub><mrow><mi>σ</mi></mrow><mrow><mi>i</mi><mo>+</mo><mn mathvariant="normal">1</mn></mrow></msub><mo stretchy="false">)</mo><mo>+</mo><mroot><mrow><mn mathvariant="normal">2</mn><msub><mrow><mi>ε</mi></mrow><mrow><mi>i</mi></mrow></msub></mrow><mrow/></mroot><mi>z</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34992228&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34992227&type= 59.60533142 4.57200003 ， <math display="block" id="M103"><msub><mrow><mi>O</mi></mrow><mrow><mi>i</mi></mrow></msub><mo>=</mo><msup><mrow><mi>H</mi></mrow><mrow><mi mathvariant="normal">T</mi></mrow></msup><mo stretchy="false">(</mo><mi>I</mi><mo>-</mo><mi>H</mi><msub><mrow><mi>O</mi></mrow><mrow><mi>i</mi></mrow></msub><mo stretchy="false">)</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34992230&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34992229&type= 30.05666542 4.65666676 ， 其中， <math id="M104"><msub><mrow><mi>ε</mi></mrow><mrow><mi>i</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34992232&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34992231&type= 2.20133328 3.72533321 是第 <math id="M105"><mi>i</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34992234&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34992233&type= 1.10066664 2.79399991 次迭代的步长。 式（14） 式（14） 和（15）称该算法的校正过程。 本文最终产生的LSGM的重建算法的伪代码如 算法1 算法1 所示。LSGM重建过程包含两层循环：（1）外层循环次数为逆时间VE-SDE的离散化步骤数，无透镜编码图像输入网络先进行一次预测，再进入内层循环校正；（2）内层循环的校正过程是通过退火朗之万迭代1 000次。整个循环迭代过程中，每一次循环均对预测和校正过程的数据先验项和数据保真项进行更新，经过1 000次的迭代实现无透镜成像重建。 图4 图4 下半部分展示了本文提出的LSGM算法的重建流程图。在训练过程，利用生成模型对随机样本训练估计数据分布的梯度先验；在重建过程，先利用训练的梯度先验进行一次预测，然后再利用退火朗之万方法迭代进行校正，并且在每一次迭代完成，都会在预测和校正步骤之后对数据保真项和先验项进行更新，以保证实现最优求解。 Algorithm 1 LSGM Training Dataset： Object image dataset Output： <math id="M106"><msub><mrow><mi>S</mi></mrow><mrow><mi>θ</mi></mrow></msub><mo stretchy="false">(</mo><mi>O</mi><mo>,</mo><mi>t</mi><mo stretchy="false">)</mo><mo>≃</mo><msub><mrow><mo>∇</mo></mrow><mrow><mi>O</mi></mrow></msub><mi mathvariant="normal">l</mi><mi mathvariant="normal">o</mi><mi mathvariant="normal">g</mi><msub><mrow><mi>p</mi></mrow><mrow><mi>t</mi></mrow></msub><mo stretchy="false">(</mo><mi>O</mi><mo stretchy="false">)</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34992236&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34992235&type= 32.51200104 4.14866638 Reconstruction Setting： Score network <math id="M107"><msub><mrow><mi>S</mi></mrow><mrow><mi>θ</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34992238&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34992237&type= 2.87866688 3.21733332 ， Noise iteration step <math id="M108"><mfenced open="{" close="}" separators="|"><mrow><msub><mrow><mi>ε</mi></mrow><mrow><mi>i</mi></mrow></msub></mrow></mfenced></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34992240&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34992239&type= 4.65666676 3.55599999 ， Noise scale <math id="M109"><mfenced open="{" close="}" separators="|"><mrow><msub><mrow><mi>σ</mi></mrow><mrow><mi>i</mi></mrow></msub></mrow></mfenced></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34992242&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34992241&type= 4.82600021 3.55599999 ，Number of discretization steps for the reverse-time SDE <math id="M110"><mi>N</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34992244&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34992243&type= 2.79399991 2.45533323 ， Number of corrector steps <math id="M111"><mi>M</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34992256&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34992255&type= 3.30200005 2.45533323 1： <math id="M112"><msub><mrow><mi>O</mi></mrow><mrow><mi>N</mi></mrow></msub><mo>∼</mo><mi>N</mi><mo stretchy="false">(</mo><mn>0</mn><mo>,</mo><msubsup><mrow><mi>σ</mi></mrow><mrow><mi mathvariant="normal">m</mi><mi mathvariant="normal">a</mi><mi mathvariant="normal">x</mi></mrow><mrow><mn>2</mn></mrow></msubsup><mi>I</mi><mo stretchy="false">)</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34992248&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34992247&type= 22.86000061 4.06400013 2： For <math id="M113"><mi>i</mi><mo>=</mo><mi>N</mi><mo>-</mo><mn>1</mn></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34992250&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34992249&type= 13.29266548 2.45533323 to 0 do （Outer loop） 3：  Update <math id="M114"><msub><mrow><mi>O</mi></mrow><mrow><mi>i</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34992252&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34992251&type= 3.21733332 3.21733332 via Eq. （ 12 12 ） and （ 13 13 ）   【Predictor】 4：  For <math id="M115"><mi>j</mi><mo>=</mo><mn>1</mn></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34992254&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34992253&type= 7.11199999 3.13266683 to <math id="M116"><mi>M</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34992256&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34992255&type= 3.30200005 2.45533323 do （Inner loop） 5：   Update <math id="M117"><msub><mrow><mi>O</mi></mrow><mrow><mi>i</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34992258&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34992257&type= 3.21733332 3.21733332 via Eq. （ 14 14 ） and （ 15 15 ） 【Corrector】 6：  End for 7： End for 8： Return <math id="M118"><msub><mrow><mi>O</mi></mrow><mrow><mn>0</mn></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34992260&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34992259&type= 3.55599999 3.21733332 3 实验结果及分析 1 在本章节中，接下来的实验将从定性和定量两方面对本文提出的方法进行评估。 3.1　数据集和模型训练 2 由于目前还没有用于无透镜编码成像的公共数据集，本实验所做的训练和测试均在LSUN数据集上进行。LSUN数据集是一个大型的彩色图像数据集，包含10个场景类别和20个对象类别，每个类别包含约100万幅标记图像。主要选取其中两个子数据集：室内场景LSUN-bedroom和教堂LSUN-church数据集，分别观察网络模型在两不同训练集上的训练和测试效果。 具体地说，使用LSUN官方网站给出的用于训练的bedroom和church数据集分别训练，图像设置大小为 <math id="M119"><mn mathvariant="normal">256</mn><mo>×</mo><mn mathvariant="normal">256</mn></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34992388&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34992387&type= 16.00200081 2.79399991 ，其像素值均缩放至 <math id="M120"><mo stretchy="false">[</mo><mn mathvariant="normal">0,1</mn><mo stretchy="false">]</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34992266&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34992265&type= 8.97466660 3.55599999 之间。使用噪声数值范围为 <math id="M121"><mn mathvariant="normal">0.01</mn><mo>~</mo><mn mathvariant="normal">380</mn></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34992268&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34992267&type= 16.42533302 2.79399991 的高斯噪声扰动数据分布。同时，为规避数值问题，参数设置 <math id="M122"><mi>ε</mi><mo>=</mo><msup><mrow><mn mathvariant="normal">10</mn></mrow><mrow><mo>-</mo><mn mathvariant="normal">5</mn></mrow></msup></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34992270&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34992269&type= 12.86933327 3.04800010 ，信噪比为 <math id="M123"><mn mathvariant="normal">0.16</mn></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34992272&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34992271&type= 7.61999989 2.79399991 。将梯度裁剪的最大值设为1.0，参数采用指数移动平均，流速为0.999，每种数据集训练50万次。此外，为了进一步验证LSGM网络的性能表现，在每个数据集中选取100张图做测试。在PC采样算法重建中的内外层循环迭代次数分别为 <math id="M124"><mi>N</mi><mo>=</mo><mn mathvariant="normal">1</mn><mtext> </mtext><mn mathvariant="normal">000</mn></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34992356&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34992355&type= 16.00200081 2.87866688 、 <math id="M125"><mi>M</mi><mo>=</mo><mn mathvariant="normal">1</mn></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34992276&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34992275&type= 10.32933331 2.87866688 ，每张图迭代1 000次，取最大的峰值信噪比及其相应的结构相似性值，然后对这100对数据取平均。基于仿真实验使用固定学习率0.000 01，优化采用自适应估计方法（Adaptive Moment Estimation，Adam），其中 <math id="M126"><msub><mrow><mi>β</mi></mrow><mrow><mn mathvariant="normal">1</mn></mrow></msub><mo>=</mo><mn mathvariant="normal">0.9</mn></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34992278&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34992277&type= 13.97000027 4.57200003 ， <math id="M127"><msub><mrow><mi>β</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msub><mo>=</mo><mn mathvariant="normal">0.999</mn></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34992280&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34992279&type= 17.61066628 4.57200003 。 3.2　评价方法 2 由于人眼的视觉特性，经常出现评价结果与人的主观感觉不一致的情况。为有效定量评价本文的实验，主要使用以下几个指标来评估重建质量： 峰值信噪比（Peak Signal to Noise Ratio，PSNR， I PSNR ）是基于对应像素点间的误差，单位为dB。该指标是基于误差敏感的图像质量评价，数值越大表示失真越小。它是使用最为广泛的一种图像客观评价指标，定义如下： <math display="block" id="M128"><msub><mrow><mi>I</mi></mrow><mrow><mi mathvariant="normal">P</mi><mi mathvariant="normal">S</mi><mi mathvariant="normal">N</mi><mi mathvariant="normal">R</mi></mrow></msub><mo stretchy="false">(</mo><mover accent="true"><mi>x</mi><mo>̂</mo></mover><mo>,</mo><mi>x</mi><mo stretchy="false">)</mo><mo>=</mo><mfrac><mrow><mi mathvariant="normal">m</mi><mi mathvariant="normal">a</mi><msub><mrow><mi mathvariant="normal">x</mi></mrow><mrow><mo stretchy="false">(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo stretchy="false">)</mo></mrow></msub><mfenced open="{" close="}" separators="|"><mrow><msup><mrow><mover accent="true"><mi>x</mi><mo>̂</mo></mover></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup><mo stretchy="false">(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo stretchy="false">)</mo></mrow></mfenced></mrow><mrow><msub><mrow><mi>I</mi></mrow><mrow><mi mathvariant="normal">M</mi><mi mathvariant="normal">S</mi><mi mathvariant="normal">E</mi></mrow></msub><mo stretchy="false">(</mo><mover accent="true"><mi>x</mi><mo>̂</mo></mover><mo>,</mo><mi>x</mi><mo stretchy="false">)</mo></mrow></mfrac></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34992282&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34992281&type= 47.15933609 9.82133293 ， <math display="block" id="M129"><msub><mrow><mi>I</mi></mrow><mrow><mi mathvariant="normal">M</mi><mi mathvariant="normal">S</mi><mi mathvariant="normal">E</mi></mrow></msub><mo>=</mo><mfrac><mrow><mn mathvariant="normal">1</mn></mrow><mrow><mi>L</mi><mo>×</mo><mi>W</mi></mrow></mfrac><mstyle displaystyle="true"><munderover><mo largeop="true">∑</mo><mrow><mi>i</mi><mo>=</mo><mn mathvariant="normal">1</mn></mrow><mrow><mi>H</mi></mrow></munderover></mstyle><mrow><mstyle displaystyle="true"><munderover><mo largeop="true">∑</mo><mrow><mi>j</mi><mo>=</mo><mn mathvariant="normal">1</mn></mrow><mrow><mi>W</mi></mrow></munderover></mstyle><mrow><msup><mrow><mfenced separators="|"><mrow><mover accent="true"><mi>x</mi><mo>̂</mo></mover><mo stretchy="false">(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo stretchy="false">)</mo><mo>-</mo><mi>x</mi><mo stretchy="false">(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo stretchy="false">)</mo></mrow></mfenced></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup></mrow></mrow></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34992284&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34992283&type= 61.12933731 9.05933285 ， 其中： <math id="M130"><mover accent="true"><mi>x</mi><mo>̂</mo></mover></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34992300&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34992299&type= 2.20133328 2.87866688 和 <math id="M131"><mi>x</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34992302&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34992301&type= 2.20133328 2.79399991 分别表示目标图像和重建图像， <math id="M132"><mi>L</mi><mo>×</mo><mi>W</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34992290&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34992289&type= 11.34533310 2.87866688 为图像 <math id="M133"><mover accent="true"><mi>x</mi><mo>̂</mo></mover></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34992300&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34992299&type= 2.20133328 2.87866688 和 <math id="M134"><mi>x</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34992302&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34992301&type= 2.20133328 2.79399991 的尺寸， <math id="M135"><mover accent="true"><mi>x</mi><mo>̂</mo></mover><mo stretchy="false">(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo stretchy="false">)</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34992296&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34992295&type= 9.82133293 3.97933316 和 <math id="M136"><mi>x</mi><mo stretchy="false">(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo stretchy="false">)</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34992298&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34992297&type= 9.82133293 3.97933316 分别为 <math id="M137"><mover accent="true"><mi>x</mi><mo>̂</mo></mover></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34992300&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34992299&type= 2.20133328 2.87866688 和 <math id="M138"><mi>x</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34992302&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34992301&type= 2.20133328 2.79399991 在 <math id="M139"><mo stretchy="false">(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo stretchy="false">)</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34992304&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34992303&type= 7.11199999 3.97933316 位置的像素。峰值信噪比数值越大表明重建图像质量越好，反之越差。 结构相似性（Structural SIMilarity，SSIM， I SSIM ）：是一种衡量两幅图像结构上相似程度的指标，它分别从宽度、对比度、结构三方面度量图像之间的相似性： <math display="block" id="M140"><msub><mrow><mi>I</mi></mrow><mrow><mi mathvariant="normal">S</mi><mi mathvariant="normal">S</mi><mi mathvariant="normal">I</mi><mi mathvariant="normal">M</mi></mrow></msub><mo stretchy="false">(</mo><mover accent="true"><mi>x</mi><mo>̂</mo></mover><mo>,</mo><mi>x</mi><mo stretchy="false">)</mo><mo>=</mo><mfrac><mrow><mo stretchy="false">(</mo><mn mathvariant="normal">2</mn><msub><mrow><mi>μ</mi></mrow><mrow><mover accent="true"><mi>x</mi><mo>̂</mo></mover></mrow></msub><msub><mrow><mi>μ</mi></mrow><mrow><mi>x</mi></mrow></msub><mo>+</mo><msub><mrow><mi>c</mi></mrow><mrow><mn mathvariant="normal">1</mn></mrow></msub><mo stretchy="false">)</mo><mo stretchy="false">(</mo><mn mathvariant="normal">2</mn><msub><mrow><mi>σ</mi></mrow><mrow><mover accent="true"><mi>x</mi><mo>̂</mo></mover><mi>x</mi></mrow></msub><mo>+</mo><msub><mrow><mi>c</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msub><mo stretchy="false">)</mo></mrow><mrow><mo stretchy="false">(</mo><msubsup><mrow><mi>μ</mi></mrow><mrow><mover accent="true"><mi>x</mi><mo>̂</mo></mover></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msubsup><mo>+</mo><msubsup><mrow><mi>μ</mi></mrow><mrow><mi>x</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msubsup><mo>+</mo><msub><mrow><mi>c</mi></mrow><mrow><mn mathvariant="normal">1</mn></mrow></msub><mo stretchy="false">)</mo><mo stretchy="false">(</mo><msubsup><mrow><mi>σ</mi></mrow><mrow><mover accent="true"><mi>x</mi><mo>̂</mo></mover></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msubsup><mo>+</mo><msubsup><mrow><mi>σ</mi></mrow><mrow><mi>x</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msubsup><mo>+</mo><msub><mrow><mi>c</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msub><mo stretchy="false">)</mo></mrow></mfrac></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34992306&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34992305&type= 66.04000092 9.39799976 ， 其中： <math id="M141"><msub><mrow><mi>μ</mi></mrow><mrow><mover accent="true"><mi>x</mi><mo>̂</mo></mover></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34992308&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34992307&type= 3.38666677 4.57200003 和 <math id="M142"><msub><mrow><mi>μ</mi></mrow><mrow><mi>x</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34992310&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34992309&type= 3.38666677 4.57200003 分别表示目标图像和重建图像的均值， <math id="M143"><msub><mrow><mi>σ</mi></mrow><mrow><mover accent="true"><mi>x</mi><mo>̂</mo></mover></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34992312&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34992311&type= 2.96333337 3.72533321 和 <math id="M144"><msub><mrow><mi>σ</mi></mrow><mrow><mi>x</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34992314&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34992313&type= 2.96333337 3.72533321 分别表示目标图像和重建图像的方差， <math id="M145"><msub><mrow><mi>σ</mi></mrow><mrow><mover accent="true"><mi>x</mi><mo>̂</mo></mover><mi>x</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34992316&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34992315&type= 3.97933316 3.72533321 表示图像的协方差， <math id="M146"><msub><mrow><mi>c</mi></mrow><mrow><mn mathvariant="normal">1</mn></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34992318&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34992317&type= 2.70933342 3.72533321 和 <math id="M147"><msub><mrow><mi>c</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34992320&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34992319&type= 2.70933342 3.72533321 是用于保持稳定的常数。结构相似性的取值范围是 <math id="M148"><mo stretchy="false">[</mo><mn mathvariant="normal">0,1</mn><mo stretchy="false">]</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34992322&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34992321&type= 8.97466660 3.55599999 ，其值越大，表示图像失真越小，越能被人眼所接受。 3.3　实验设置及环境 2 本实验的模拟仿真实验参数设置具体如下：FZA作为编码掩膜放置在图像传感器 <math id="M149"><mn mathvariant="normal">3</mn><mtext> </mtext><mi mathvariant="normal">m</mi><mi mathvariant="normal">m</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34992324&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34992323&type= 8.55133343 2.87866688 处，目标大小为15 mm×15 mm，置于距离编码掩膜20 mm处，掩膜板最内侧区域半径 <math id="M150"><msub><mrow><mi>r</mi></mrow><mrow><mn mathvariant="normal">1</mn></mrow></msub><mo>=</mo><mn mathvariant="normal">0.23</mn></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34992326&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34992325&type= 15.07066536 3.72533321 mm。模拟传感器的像素数为 <math id="M151"><mn mathvariant="normal">256</mn><mo>×</mo><mn mathvariant="normal">256</mn></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34992388&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34992387&type= 16.00200081 2.79399991 ，像素间距为 <math id="M152"><mn mathvariant="normal">0.014</mn><mtext> </mtext><mi mathvariant="normal">μ</mi><mi mathvariant="normal">m</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34992330&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34992329&type= 15.15533447 3.64066648 。传感器上每个像素探测器记录的是编码掩膜上多个孔径所成图像的一个叠加，即来自场景中多个位置入射光场的叠加。 软件仿真环境为：实验使用基于PyTorch的深度学习框架，在一张2080 TI的GPU上进行LSGM算法的训练和测试。为了更好地体现算法的可重复性，我们将代码公开在网址： https：//github.com/yqx7150/LSGM. https://github.com/yqx7150/LSGM. 3.4　实验结果 2 3 . 4 . 1　LSUN-bedroom数据集上的实验结果及分析 3 近年来，随着数据量的急剧增长和计算能力的不断提高，基于分数匹配生成模型的成像方法已经被应用于磁共振成像 ［ 30 30 ］ 和CT成像重建 ［ 31 31 ］ ，然而尚未应用于无透镜成像。为证明本实验的有效性，本实验挑选了传统的BP算法和目前实验结果较好的CS算法 ［ 16 16 ］ 作为对比算法。 为从视觉上比较， 图5 图5 从LSUN-bedroom数据集中挑选出四组不同实验结果进行展示，其中每组实验结果包含原图、编码图和BP、CS、LSGM算法重建的图。可以很直观地看出，BP算法实验结果最差，可以观察到很明显的孪生像效应。其次，CS算法的重建结果与原图对比，存在明显的颜色差异，而且其图片清晰度存在很大差距。例如，在第二行和第三行中CS重建结果与原始图像有很明显的颜色差距，中间有很明显的颜色扩散以及图片边缘黑色区域变色现象。显然， 图5 图5 中第五列的图像质量明显优于另外两组结果，三种算法中LSGM的重建结果最优。 为进一步进行定量分析，三个算法的平均PSNR和SSIM值如 表1 表1 所示。从表中分析可知，传统的BP算法的平均PSNR值仅为7.78 dB。相比之下，其余两种算法通过迭代的方式可以获得更好的结果。对比表中数据，CS和LSGM的PSNR值分别比BP算法高10.69 dB、18.32 dB。本文提出的LSGM算法具有最优的成像质量，在迭代校正之后，其PSNR值为26.10 dB、SSIM值能达到0.64，比CS算法的PSNR和SSIM分别高出7.63 dB和0.33。 注： 加粗为最优结果 Method PSNR/dB SSIM BP 7.78 0.23 CS 18.47 0.31 LSGM 26.10 0.64 3 . 4 . 2　LSUN-church数据集上的实验结果及分析 3 为了更好地验证LSGM方法的鲁棒性，本实验还选取了另一个图像灰度、亮度和结构变化更丰富的LSUN-church（ <math id="M153"><mn mathvariant="normal">256</mn><mo>×</mo><mn mathvariant="normal">256</mn></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34992337&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34992336&type= 16.17133331 2.79399991 ）数据集进行实验。在使用LSUN-church数据集对LSGM模型训练完成之后，选取其中100张图像作为测试输入。 图6 图6 展示了BP、CS和本文所提出的LSGM三种算法在LSUN-church数据集上重建图像的视觉比较。和上小节实验相同，本实验展示了四组不同的重建图像，每组包含原图、编码图和BP、CS、LSGM算法重建结果。与在LSUN-bedroom数据集上的实验结果相似，BP算法重建结果最差；CS算法重建的图像上总是会有与图片不相关的颜色出现。例如， 图6 图6 中第一、二和四行的CS重建结果中能观察到明显的紫色和绿色孪生像，尤其是图片中间区域和边缘黑色区域更为明显。虽然CS算法重建结果并不影响图像的辨识，但图像质量仍然较模糊。对比 图6 图6 中第三行的CS和LSGM算法重建结果，可发现LSGM算法具有更高的清晰度。根据 图6 图6 中给出对应的PSNR和SSIM，LSGM算法重建的定量评价指标也是最高的，LSGM的PSNR和SSIM比CS算法分别提高6.68 dB和0.36。整体对比 图6 图6 中各组重建结果，最后一列图像不论在图像颜色、亮度和清晰度上都优于第三列和第四列。因此，LSGM算法是表现最好的，它在保证颜色的同时，还能清晰地重建出图像细节，更加接近原始图像。 定量对比 表2 表2 中PSNR和SSIM两项指标数据，在三个算法中，LSGM的PSNR和SSIM值最高，在该实验数据集上的平均PSNR值接近于LSUN-bedroom数据集上的结果。在更换成一个灰度、亮度和结构变化更为丰富的数据集后，BP算法重建记录的评价指标几乎不变。CS算法的PSNR值下降了0.68 dB，而SSIM值下降0.02。整体分析，本文提出的LSGM基算法在两项指标上均表现得相对稳定且最好，PSNR值最高达到24.47 dB、SSIM值达到0.65。相对BP和CS算法，LSGM算法的平均PSNR值依旧分别提高16.66 dB和6.68 dB，平均SSIM值分别提高0.43和0.36。因此，本文提出的方法重建图像最好，不仅能够在保持颜色、亮度和清晰度的同时，还能显著去除孪生像的伪影噪声。 注： 加粗为最优结果 Method PSNR/dB SSIM BP 7.81 0.22 CS 17.79 0.29 LSGM 24.47 0.65 综合上述两组实验结果可以得出结论，在不同数据集上，与BP算法、CS算法相比，LSGM方法都能达到更高的重建质量，模型性能更加稳定。 3 . 4 . 3　TV去网格效应 3 由于无透镜成像本身特性，重建图像存在很明显的“网格”现象，重建图像和原始图像存在一定差距，在图像纹理细节部分的恢复结果仍然不够理想。为了解决这个问题，本文将基于变分理论的TV方法应用于重建过程。该方法是依靠梯度下降流对图像进行平滑，希望在图像内部尽可能对图像进行平滑（相邻像素的差值较小），而在图像边缘（图像轮廓）尽可能不去平滑。前两个实验证明了本文提出的LSGM模型生成的先验信息作为正则化项求解无透镜成像逆问题，并取得了较好的数值结果。 图7 图7 展示了使用及不使用TV方法的LSGM算法的重建结果以及原始图像。从 图7 图7 （c）可以看出，结合了TV方法的重建图像中“网格”效应更为平滑，视觉效果更令人舒适，但其相应代价是失去了图像中的部分高频信息。由两种方法重建结果的定量指标可以得出，提出的LSGM算法比不使用TV方法的LSGM算法的PSNR和SSIM值分别平均提高0.80 dB和0.10。 4 讨 论 1 4.1　消融实验 2 上述实验表明，本文提出的LSGM方法能够实现高质量的无透镜成像结果。然而，使用基于分数匹配生成模型进行图像重建时，迭代式重建方法在耗时方面有着明显的局限性。为了更好地验证LSGM模型的性能，本实验探究了在不同 <math id="M154"><msub><mrow><mi>σ</mi></mrow><mrow><mi mathvariant="normal">m</mi><mi mathvariant="normal">a</mi><mi mathvariant="normal">x</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34992348&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34992347&type= 5.33400011 3.72533321 下的图像质量和离散化步骤数之间的关系。 如 图8 图8 所示，本实验选取了LSUN-church数据集对测试需要的离散化步骤数与图像质量之间的关系进行实验验证。当训练阶段使用 <math id="M155"><mi>N</mi><mo>=</mo><mn mathvariant="normal">1</mn><mtext> </mtext><mn mathvariant="normal">000</mn></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34992350&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34992349&type= 16.00200081 2.87866688 个离散化步骤来进行正向SDE，重建阶段使用 <math id="M156"><mi>N</mi><mo>=</mo><mn mathvariant="normal">1</mn><mtext> </mtext><mn mathvariant="normal">000</mn></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34992352&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34992351&type= 15.83266735 2.87866688 个离散化步骤时，需要大约10 min的推理时间。从 <math id="M157"><mi>N</mi><mo>=</mo><mn mathvariant="normal">1</mn></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34992354&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34992353&type= 9.65200043 2.87866688 到 <math id="M158"><mi>N</mi><mo>=</mo><mn mathvariant="normal">1</mn><mtext> </mtext><mn mathvariant="normal">000</mn></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34992356&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34992355&type= 16.00200081 2.87866688 每隔1个离散化步骤记录重建图像的PSNR以及SSIM。可以观察到随着离散化步骤的增加，重建图像PSNR以及SSIM呈现一定的规律性。在 <math id="M159"><mi>N</mi><mo>=</mo><mn mathvariant="normal">700</mn></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34992358&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34992357&type= 13.71599960 2.87866688 的情况下，可以实现高质量的重建图像，PSNR以及SSIM达到最高值。当离散化步骤数逐渐增加，PSNR值趋于稳定，SSIM会下降。由于在不同 <math id="M160"><msub><mrow><mi>σ</mi></mrow><mrow><mi mathvariant="normal">m</mi><mi mathvariant="normal">a</mi><mi mathvariant="normal">x</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34992360&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34992359&type= 5.33400011 3.72533321 下，PSNR随着离散步骤数增加，将趋向于稳定，并且稳定值在24 dB左右，而SSIM波动较大，所以本文选取当PSNR达到稳定值下的离散步骤数来进行所有实验。在之后的工作中，我们将尝试使用去噪扩散模型（DDPM） ［ 32 32 ］ 来进一步减少离散化步骤数。 4.2　成像视场分析 2 如 图9 图9 所示，为测试视场角对本文提出算法的影响，进一步分析了不同视场角下重建图像的PSNR和SSIM指标。在本组实验中，FZA的参数与前文一致，在保持目标尺寸不变的情况下，通过改变目标物与FZA掩膜间的距离来实现不同视场角。从 图9 图9 所示结果可以看出，随着视场角的增大，PSNR的值先增大后减小，随着目标图像的增大，边界处部分衍射信息丢失，造成图像重建质量的恶化。当视场角在 <math id="M163"><mn mathvariant="normal">32</mn><mo>°</mo><mo>∼</mo><mn mathvariant="normal">46</mn><mo>°</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34992369&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34992368&type= 13.80066681 3.04800010 范围内，相应的PSNR的值可达到 <math id="M164"><mn mathvariant="normal">26</mn><mo>∼</mo><mn mathvariant="normal">28</mn><mtext> </mtext></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34992371&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34992370&type= 12.53066730 2.87866688 dB、SSIM的值为 <math id="M165"><mn mathvariant="normal">0.85</mn><mo>∼</mo><mn mathvariant="normal">0.68</mn></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34992373&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34992372&type= 19.38866615 2.79399991 ，可以实现最好的成像效果。 4.3　成像分辨率分析 2 受制于工艺制造水平，传感器和FZA掩膜的像元尺寸不能无限小，系统重建的带宽有限。为了验证提出算法能实现的最小分辨率，本文采用一张尺寸为 <math id="M166"><mn mathvariant="normal">20</mn><mtext> </mtext><mi mathvariant="normal">c</mi><mi mathvariant="normal">m</mi><mo>×</mo><mn mathvariant="normal">20</mn><mtext> </mtext><mi mathvariant="normal">c</mi><mi mathvariant="normal">m</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34992378&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34992377&type= 22.86000061 2.87866688 、像素数为 <math id="M167"><mn mathvariant="normal">1</mn><mtext> </mtext><mn mathvariant="normal">280</mn><mo>×</mo><mn mathvariant="normal">1</mn><mtext> </mtext><mn mathvariant="normal">280</mn></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34992380&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34992379&type= 21.25133324 2.87866688 的USAF-1951测试图像进行测试，该成像目标与FZA距离为 <math id="M168"><mn mathvariant="normal">267</mn><mtext> </mtext><mi mathvariant="normal">m</mi><mi mathvariant="normal">m</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34992382&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34992381&type= 12.27666664 2.87866688 ，成像视场角为 <math id="M169"><mn mathvariant="normal">41.11</mn><mo>°</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34992384&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34992383&type= 10.49866676 2.96333337 。由于本文训练模型中传感器的像素数为 <math id="M170"><mn mathvariant="normal">256</mn><mo>×</mo><mn mathvariant="normal">256</mn></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34992388&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34992387&type= 16.00200081 2.79399991 ，为了适应本文算法，将USAF-1951原图经过FZA调制后的结果进行像素归并至像素数为 <math id="M171"><mn mathvariant="normal">256</mn><mo>×</mo><mn mathvariant="normal">256</mn></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34992388&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34992387&type= 16.00200081 2.79399991 ，以得到测量结果。 图像重建结果如 图10 图10 所示， 图10 图10 （a）为USAF-1951测试卡原图、 - 1组线对的细节以及相应的剖线。为了测试不同FZA参数 r 1 对结果的影响，本文在 <math id="M172"><msub><mrow><mi>r</mi></mrow><mrow><mn mathvariant="normal">1</mn></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34992396&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34992395&type= 2.62466669 3.72533321 分别为0.35、0.25、0.23、0.20 mm时进行了重建。可以很明显看出， <math id="M173"><msub><mrow><mi>r</mi></mrow><mrow><mn mathvariant="normal">1</mn></mrow></msub><mo>=</mo><mn mathvariant="normal">0.35</mn><mtext> </mtext><mi mathvariant="normal">m</mi><mi mathvariant="normal">m</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34992392&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34992391&type= 21.92866707 3.80999994 时，难以分辨 - 1组线对。随着 <math id="M174"><msub><mrow><mi>r</mi></mrow><mrow><mn mathvariant="normal">1</mn></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34992396&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34992395&type= 2.62466669 3.72533321 参数的减少，图像分辨率得到有效提高，当 <math id="M175"><msub><mrow><mi>r</mi></mrow><mrow><mn mathvariant="normal">1</mn></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34992396&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34992395&type= 2.62466669 3.72533321 减小到 <math id="M176"><mn mathvariant="normal">0.23</mn><mtext> </mtext><mi mathvariant="normal">m</mi><mi mathvariant="normal">m</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34992398&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34992397&type= 14.22399998 2.87866688 及 <math id="M177"><mn mathvariant="normal">0.20</mn><mtext> </mtext><mi mathvariant="normal">m</mi><mi mathvariant="normal">m</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34992400&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=34992399&type= 14.22399998 2.87866688 时，能够有效分辨出 - 1组4号线对。 5 结 论 1 本文提出了一种基于分数匹配生成模型框架的无透镜成像方法，即LSGM。该模型利用数据分布梯度对图像进行建模，并从给定的观测数据的条件分布中采样数据，然后将采样的数据作为图像先验数据用于图像重建。所提出的LSGM算法在计算重建过程中能够有效消除无透镜成像系统固有的病态性而带来的孪生像伪影噪声及其他影响成像的不良影响，使得成像质量大幅提高。在LSUN-bedroom和LSUN-church数据集上测试结果表明，本文提出的方法和其他传统算法相比，LSGM算法图像重建的平均PSNR值和平均SSIM值分别达到25.23 dB和0.65，量化指标以及视觉质量均获得了令人满意的结果。本文提出的模型框架结构灵活，训练稳定，在建模数据分布方面的性能表现出色，能够获得高质量样本。本文通过大量的实验验证了该方法在质量和实用性方面的优越性，有望应用于其他计算成像领域。 参考文献 SAHOO S K ， TANG D L ， DANG C . Single-shot multispectral imaging with a monochromatic camera ［J］. Optica ， 2017 ， 4 （ 10 ）： 1209 - 1213 . doi: 10.1364/optica.4.001209 http://dx.doi.org/10.1364/optica.4.001209 Single-shot multispectral imaging with a monochromatic camera BOOMINATHAN V ， ADAMS J K ， ROBINSON J T ， et al . PhlatCam： designed phase-mask based thin lensless camera ［J］. IEEE Transactions on Pattern Analysis and Machine Intelligence ， 2020 ， 42 （ 7 ）： 1618 - 1629 . doi: 10.1109/tpami.2020.2987489 http://dx.doi.org/10.1109/tpami.2020.2987489 PhlatCam： designed phase-mask based thin lensless camera 王雪 ， 刘虹遥 ， 路鑫超 ， 等 . 无透镜全息显微细胞成像 ［J］. 光学 精密工程 ， 2020 ， 28 （ 8 ）： 1644 - 1650 . 无透镜全息显微细胞成像 WANG X ， LIU H Y ， LU X C ， et al . Cell imaging by holographic lens-free microscopy ［J］. Opt. Precision Eng. ， 2020 ， 28 （ 8 ）： 1644 - 1650 . （in Chinese） Cell imaging by holographic lens-free microscopy CHAPMAN H N ， NUGENT K A . Coherent lensless X-ray imaging ［J］. Nature Photonics ， 2010 ， 4 （ 12 ）： 833 - 839 . doi: 10.1038/nphoton.2010.240 http://dx.doi.org/10.1038/nphoton.2010.240 Coherent lensless X-ray imaging WITTE S ， TENNER V T ， NOOM D W ， et al . Lensless diffractive imaging with ultra-broadband table-top sources： from infrared to extreme-ultraviolet wavelengths ［J］. Light： Science & Applications ， 2014 ， 3 （ 3 ）： e163 . doi: 10.1038/lsa.2014.44 http://dx.doi.org/10.1038/lsa.2014.44 Lensless diffractive imaging with ultra-broadband table-top sources： from infrared to extreme-ultraviolet wavelengths BISHARA W ， SU T W ， COSKUN A F ， et al . Lensfree on-chip microscopy over a wide field-of-view using pixel super-resolution ［J］. Optics Express ， 2010 ， 18 （ 11 ）： 11181 - 11191 . doi: 10.1364/oe.18.011181 http://dx.doi.org/10.1364/oe.18.011181 Lensfree on-chip microscopy over a wide field-of-view using pixel super-resolution GREENBAUM A ， LUO W ， SU T W ， et al . Imaging without lenses： achievements and remaining challenges of wide-field on-chip microscopy ［J］. Nature Methods ， 2012 ， 9 （ 9 ）： 889 - 895 . doi: 10.1038/nmeth.2114 http://dx.doi.org/10.1038/nmeth.2114 Imaging without lenses： achievements and remaining challenges of wide-field on-chip microscopy MERTZ L ， YOUNG N O . Fresnel transformations of images ［J］. SPIE milestone series ms ， 1996 ， 128 ： 44 - 49 . Fresnel transformations of images 万玉红 ， 刘超 ， 满天龙 ， 等 . 非相干全息术成像特性及研究进展 ［J］. 激光与光电子学进展 ， 2021 ， 58 （ 18 ）： 1 - 24 . 非相干全息术成像特性及研究进展 WAN Y H ， LIU C ， MAN T L ， et al . Incoherent correlation digital holography： principle， development， and applications ［J］. Laster & Optoelectronics Progress ， 2021 ， 58 （ 18 ）： 1 - 24 . （in Chinese） Incoherent correlation digital holography： principle， development， and applications LATYCHEVSKAIA T ， FINK H W . Solution to the twin image problem in holography ［J］. Physical Review Letters ， 2007 ， 98 （ 23 ）： 233901 . doi: 10.1103/physrevlett.98.233901 http://dx.doi.org/10.1103/physrevlett.98.233901 Solution to the twin image problem in holography HENNELLY B ， KELLY D ， PANDEY N ， et al . Review of twin reduction and twin removal techniques in holography ［J］. in Procoodings of CIICT ， 2009 ： 241 - 245 . Review of twin reduction and twin removal techniques in holography DENIS L ， FOURNIER C ， FOURNEL T ， et al . Twin-image noise reduction by phase retrieval in in-line digital holography ［C］. Optics and Photonics 2005 . Proc SPIE 5914， Wavelets XI， San Diego， California， USA . 2005， 5914 ： 148 - 161 . doi: 10.1117/12.617405 http://dx.doi.org/10.1117/12.617405 Twin-image noise reduction by phase retrieval in in-line digital holography SHIMANO T ， NAKAMURA Y ， TAJIMA K ， et al . Lensless light-field imaging with Fresnel zone aperture： quasi-coherent coding ［J］. Applied Optics ， 2018 ， 57 （ 11 ）： 2841 - 2850 . doi: 10.1364/ao.57.002841 http://dx.doi.org/10.1364/ao.57.002841 Lensless light-field imaging with Fresnel zone aperture： quasi-coherent coding SAO M ， NAKAMURA Y ， TAJIMA K ， et al . Lensless close-up imaging with Fresnel zone aperture ［J］. Japanese Journal of Applied Physics ， 2018 ， 57 （ 9S1 ）： 09 SB05. doi: 10.7567/jjap.57.09sb05 http://dx.doi.org/10.7567/jjap.57.09sb05 Lensless close-up imaging with Fresnel zone aperture BRADY D J ， CHOI K ， MARKS D L ， et al . Compressive holography ［J］. Optics Express ， 2009 ， 17 （ 15 ）： 13040 - 13049 . doi: 10.1364/oe.17.013040 http://dx.doi.org/10.1364/oe.17.013040 Compressive holography WU J ， ZHANG H ， ZHANG W H ， et al . Single-shot lensless imaging with Fresnel zone aperture and incoherent illumination ［J］. Light： Science & Applications ， 2020 ， 9 ： 53 . doi: 10.1038/s41377-020-0289-9 http://dx.doi.org/10.1038/s41377-020-0289-9 Single-shot lensless imaging with Fresnel zone aperture and incoherent illumination HYVARINEN A ， DAYAN P . Estimation of non-normalized statistical models by score matching ［J］. Journal of Machine Learning Research ， 2005 ， 6 （ 4 ）. Estimation of non-normalized statistical models by score matching SONG Y ， ERMON S . Generative modeling by estimating gradients of the data distribution ［J］. Advances in Neural Information Processing Systems ， 2019 ， 32 . Generative modeling by estimating gradients of the data distribution SONG Y ， GARG S ， SHI J ， et al . Sliced score matching： a scalable approach to density and score estimation ［C］. Uncertainty in Artificial Intelligence. PMLR ， 2020 ： 574 - 584 . Sliced score matching： a scalable approach to density and score estimation SOHL-DICKSTEIN J ， WWISS E ， MAHESWARANATHAN N ， et al . Deep unsupervised learning using nonequilibrium thermodynamics ［C］. International Conference on Machine Learning. PMLR ， 2015 ： 2256 - 2265 . Deep unsupervised learning using nonequilibrium thermodynamics HO J ， JAIN A ， ABBEEL P . Denoising diffusion probabilistic models ［J］. Advances in Neural Information Processing Systems ， 2020 ， 33 ： 6840 - 6851 . Denoising diffusion probabilistic models SONG Y ， SOHL-DICKSTEIN J ， KINGMA D P ， et al . Score-based generative modeling through stochastic differential equations ［EB/OL］. 2020： arXiv ： 2011 .13456［cs.LG］. https：//arxiv.org/abs/2011.13456 https://arxiv.org/abs/2011.13456 Score-based generative modeling through stochastic differential equations NICHOL A Q ， DHARIWAL P . Improved denoising diffusion probabilistic models ［C］. International Conference on Machine Learning. PMLR ， 2021 ： 8162 - 8171 . Improved denoising diffusion probabilistic models DHARIWAL P ， NICHOL A . Diffusion models beat GANs on image synthesis ［J］. Advances in Neural Information Processing Systems ， 2021 ， 34 ： 8780 - 8794 . Diffusion models beat GANs on image synthesis ANDERSON B D O . Reverse-time diffusion equation models ［J］. Stochastic Processes and Their Applications ， 1982 ， 12 （ 3 ）： 313 - 326 . doi: 10.1016/0304-4149(82)90051-5 http://dx.doi.org/10.1016/0304-4149(82)90051-5 Reverse-time diffusion equation models VINCENT P . A connection between score matching and denoising autoencoders ［J］. Neural Computation ， 2011 ， 23 （ 7 ）： 1661 - 1674 . doi: 10.1162/neco_a_00142 http://dx.doi.org/10.1162/neco_a_00142 A connection between score matching and denoising autoencoders BOYD S ， PARIKH N ， CHU E ， et al .. Distributed optimization and statistical learning via the alternating direction method of multipliers ［J］. Foundations and Trends ® in Machine learning ， 2011 ， 3 （ 1 ）： 1 - 122 . Distributed optimization and statistical learning via the alternating direction method of multipliers SAMSONOV A A ， KHOLMOVSKI E G ， PARKER D L ， et al . POCSENSE： POCS-based reconstruction for sensitivity encoded magnetic resonance imaging ［J］. Magnetic Resonance in Medicine ， 2004 ， 52 （ 6 ）： 1397 - 1406 . doi: 10.1002/mrm.20285 http://dx.doi.org/10.1002/mrm.20285 POCSENSE： POCS-based reconstruction for sensitivity encoded magnetic resonance imaging PARISI G . Correlation functions and computer simulations ［J］. Nuclear Physics B ， 1981 ， 180 （ 3 ）： 378 - 384 . doi: 10.1016/0550-3213(81)90056-0 http://dx.doi.org/10.1016/0550-3213(81)90056-0 Correlation functions and computer simulations QUAN C ， ZHOU J J ， ZHU Y Z ， et al . Homotopic gradients of generative density priors for MR image reconstruction ［J］. IEEE Transactions on Medical Imaging ， 2021 ， 40 （ 12 ）： 3265 - 3278 . doi: 10.1109/tmi.2021.3081677 http://dx.doi.org/10.1109/tmi.2021.3081677 Homotopic gradients of generative density priors for MR image reconstruction HE Z N ， ZHANG Y K ， GUAN Y ， et al . Iterative reconstruction for low-dose CT using deep gradient priors of generative model ［J］. IEEE Transactions on Radiation and Plasma Medical Sciences ， 2022 ， 6 （ 7 ）： 741 - 754 . doi: 10.1109/trpms.2022.3148373 http://dx.doi.org/10.1109/trpms.2022.3148373 Iterative reconstruction for low-dose CT using deep gradient priors of generative model SONG J ， MENG C ， ERMON S . Denoising diffusion implicit models ［J］. arXiv preprint arXiv ： 2010.02502 ， 2020 . Denoising diffusion implicit models 10.37188/OPE.20223018.2280 TP751 A 1004-924X（2022）18-2280-15 伍春花,彭鸿,刘且根等.基于分数匹配生成模型的无透镜成像方法[J].光学精密工程,2022,30(18):2280-2294. 伍春花,彭鸿,刘且根等.基于分数匹配生成模型的无透镜成像方法[J].光学精密工程,2022,30(18):2280-2294. WU Chunhua,PENG Hong,LIU Qiegen,et al.Lens-less imaging via score-based generative model[J].Optics and Precision Engineering,2022,30(18):2280-2294. WU Chunhua,PENG Hong,LIU Qiegen,et al.Lens-less imaging via score-based generative model[J].Optics and Precision Engineering,2022,30(18):2280-2294. 国家优秀青年科学基金项目（No. 62122033）；国家自然科学基金青年科学基金项目（No.62105138）；国家自然科学基金面上项目（No. 61871206） 2022-03-25 2022-06-07 2022-09-25 2022-09-26 光学精密工程 18 30