光学系统与图像处理端到端协同设计及其应用

光学系统与图像处理端到端协同设计及其应用 End-to-end co-design of optics and image processing and its applications first-author 顿 雄（1986-），男，博士，特聘研究员，博士生导师，2007年于哈尔滨工业大学获得学士学位，2016年于北京理工大学获得博士学位，主要从事光学系统设计、计算成像和智能感知的研究。E-mail： dunx@tongji.edu.cn http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=37126027&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=37126036&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=37126032&type= 22.01333237 32.00399780 顿 雄 （1986-），男，博士，特聘研究员，博士生导师，2007年于哈尔滨工业大学获得学士学位，2016年于北京理工大学获得博士学位，主要从事光学系统设计、计算成像和智能感知的研究。E-mail： dunx@tongji.edu.cn dunx@tongji.edu.cn corresp E-mail： chengxb@tongji.edu.cn E-mail： chengxb@tongji.edu.cn 程鑫彬（1980-），男，博士，教授，博士生导师，2004年、2008年于同济大学分别获得学士、博士学位，主要从事微纳光学、智能感知研究。E-mail： chengxb@tongji.edu.cn http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=37126041&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=37126052&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=37126045&type= 22.01332855 32.00400162 程鑫彬 （1980-），男，博士，教授，博士生导师，2004年、2008年于同济大学分别获得学士、博士学位，主要从事微纳光学、智能感知研究。E-mail： chengxb@tongji.edu.cn chengxb@tongji.edu.cn 同济大学 物理科学与工程学院 精密光学工程技术研究所，先进微结构材料教育部重点实验室，上海市数字光学前沿科学研究基地，上海市全光谱高性能光学薄膜器件与应用专业技术服务平台，上海 200092 同济大学 物理科学与工程学院 精密光学工程技术研究所，先进微结构材料教育部重点实验室，上海市数字光学前沿科学研究基地，上海市全光谱高性能光学薄膜器件与应用专业技术服务平台，上海 200092 Institute of Precision Optical Engineering， MOE Key Laboratory of Advanced Micro-Structured Materials， Shanghai Frontiers Science Center of Digital Optics， Shanghai Professional Technical Service Platform for Full-Spectrum and High-Performance Optical Thin Film Devices and Applications， School of Physics Science and Engineering， Tongji University， Shanghai 200092， China Institute of Precision Optical Engineering， MOE Key Laboratory of Advanced Micro-Structured Materials， Shanghai Frontiers Science Center of Digital Optics， Shanghai Professional Technical Service Platform for Full-Spectrum and High-Performance Optical Thin Film Devices and Applications， School of Physics Science and Engineering， Tongji University， Shanghai 200092， China 计算成像是一种通过联合光学系统和图像处理来实现特定成像功能的新兴研究领域，长期以来，计算成像中的光学与算法联合都采取的是顺序设计模式，即光学系统与图像处理各自分开设计，但这样难以全面发挥二者协同的优势。近年来，随着深度学习技术的飞速发展，基于深度学习架构的光学系统与图像处理算法端到端协同设计方法开启了解决这一问题的大门。一方面，端到端协同设计通过全面探索整个解空间，可以实现光学与图像处理的自动最佳协同；另一方面，端到端协同设计更使得研制基于任务的最优成像系统成为可能。本文首先介绍了基于光学系统与图像处理端到端协同设计框架的进展，然后介绍了我们基于这一框架在平面透镜的宽谱成像、平面透镜的大视场成像、大景深成像、超分辨成像和快照式光谱成像方面的研究进展。 Computational imaging is an emerging research field that combines optical systems and image processing to achieve specific imaging features. For a long time， this kind of combination in computational imaging has adopted a sequential design mode， that is， the optical system and image processing are designed separately， making it difficult to exploit the synergistic advantages of both entirely. With the rapid development of deep learning， the end-to-end co-design method of optics and image processing algorithms based on deep learning architecture has opened the door to solving this problem. On the one hand， end-to-end co-design can realise the automatic optimal collaboration of optics and image processing by comprehensively exploring the entire solution space. On the other hand， end-to-end design makes it possible to develop specific optical imaging systems. In this paper， we review recent advances in end-to-end co-design of optics and image processing， including full-colour imaging of flat lenses， large field of view imaging， extended depth of field imaging， and super-resolution imaging， and snapshot spectral imaging. 计算成像 端到端设计 基于任务的成像系统 computational imaging end-to-end design task domain imaging system 1 引 言 1 光学成像系统作为人类视觉功能延伸扩展的工具，广泛应用于军事和民用领域。传统成像系统是“所见即所得”的成像模式 ［ 1 1 ］ ，受光学系统设计、强度成像机理、成像衍射极限等因素制约，在像差矫正、空间分辨、光谱分辨、信息维度与探测灵敏度等方面存在很大局限性。近年来，随着光学、电子学、信息处理、计算方法和高性能计算机的发展，计算成像技术应运而生 ［ 2 2 - 4 4 2-4 ］ ，它打破了传统“所见即所得”的模式，采取了更加灵活的非直接成像模式，通过联合光学系统和图像处理来实现最终需要的成像功能。这种灵活的设计模式可以改变光学成像采集信息的模式，并平衡光学硬件空间和图像处理算法软件空间所依赖的资源，从而赋予其诸多传统成像系统难以获得的革命性优势，如简化系统（单镜片成像 ［ 5 5 - 7 7 5-7 ］ ），获得额外功能（光谱成像 ［ 8 8 ］ ），突破设备原有分辨率限制（超分辨成像 ［ 9 9 ］ ）等。但目前的计算成像系统多采取顺序设计模式，即光学系统和图像处理算法依据经验规则独立设计，这种人为设定的经验规则很多时候并非最优，特别是从计算成像自身的角度来说，人为制定的规则可能并不能与后续图像处理算法相匹配，而且这种顺序设计模式难以考虑到诸如加工误差等不确定因素的影响。总的来说，基于顺序设计方法的计算成像系统对整体优化空间探索不足，性能和鲁棒性都有很大的提升空间。 因此，人们也一直在探索光学系统与图像处理算法协同设计的方法，但由于传统光学设计过程是不可自动微分的，传统优化设计理念下难以建立起它们之间协同优化架构，进展有限。2018年，在深度学习优化模型的启发下，我们建立了针对单个光学元件的光学系统与图像处理算法协同优化框架 ［ 5 5 ］ ，从而可以真正意义上开展光学系统与图形处理算法的协同设计。一方面，端到端协同设计通过突破光学设计和图像后处理之间的壁垒，全面探索整个解空间，找到光学和算法部分在硬件成本、加工可行性、体积重量、成像质量、算法复杂度以及特殊功能间的最佳折中，从而实现在设计要求下的最优方案，极大地提升成像系统的性能。另一方面，端到端协同设计使得研制基于特定任务的最优成像系统成为可能，有望实现成像系统的多样性和多功能化。自端到端光学联合设计（End-to-end Co-design）策略提出以来，迅速吸引了学界的广泛关注，经过近几年的发展已成为计算光学研究的热点分支，在消色差成像 ［ 5 5 - 6 6 5-6 ， 10 10 ］ 、大景深成像 ［ 5 5 ， 11 11 - 14 14 11-14 ］ 、光谱成像 ［ 8 8 ， 15 15 - 16 16 15-16 ］ 、单相机三维成像 ［ 17 17 - 19 19 17-19 ］ 、高动态范围成像 ［ 20 20 - 21 21 20-21 ］ 、超分辨成像 ［ 5 5 ， 7 7 ， 9 9 ］ 等方面都取得了广泛应用。 本文主要介绍基于端到端协同设计技术的基本架构，包括可微分成像模型和计算降维方法。此外，介绍了我们基于这一框架在平面透镜的宽谱成像、平面透镜的大视场成像、大景深成像，超分辨成像和快照式光谱成像方面的研究进展，其中有部分内容是作者在沙特阿卜杜拉国王科技大学的工作。 2 光学系统与图像处理算法端到端设计基本架构 1 光学系统与图像处理端到端协同设计架构如 图1 图1 所示，其包括成像模块、重建模块以及误差模块。长期以来，端到端设计难以实现的根本原因在于缺乏可微分的成像模块。2018年，我们基于波动光学推导针对单透镜的可微分成像模块 ［ 5 5 ］ ，并于2021年又发展了基于可微分射线追踪的成像模型 ［ 22 22 ］ 。在重建模块和误差模块上我们也发展了众多方法。 2.1　可微分成像模块 2 如 图1 图1 所示，光学成像模块主要包括光学系统点扩散函数计算过程 <math id="M1"><msub><mrow><mi>f</mi></mrow><mrow><mi mathvariant="normal">o</mi><mi mathvariant="normal">p</mi><mi mathvariant="normal">t</mi><mi mathvariant="normal">i</mi><mi mathvariant="normal">c</mi><mi mathvariant="normal">s</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=37126078&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=37126074&type= 6.68866634 4.57200003 ，点扩散函数（Point Spread Function，PSF）与目标图像卷积，以及探测器积分采样过程 <math id="M2"><msub><mrow><mi>f</mi></mrow><mrow><mi mathvariant="normal">s</mi><mi mathvariant="normal">e</mi><mi mathvariant="normal">n</mi><mi mathvariant="normal">s</mi><mi mathvariant="normal">o</mi><mi mathvariant="normal">r</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=37126200&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=37126198&type= 7.02733374 4.57200003 三部分。PSF计算过程 <math id="M3"><msub><mrow><mi>f</mi></mrow><mrow><mi mathvariant="normal">o</mi><mi mathvariant="normal">p</mi><mi mathvariant="normal">t</mi><mi mathvariant="normal">i</mi><mi mathvariant="normal">c</mi><mi mathvariant="normal">s</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=37126099&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=37126095&type= 6.68866634 4.57200003 采取波动光学理论推导得到，具体如下： <math display="block" id="M4"><msub><mrow><mi>p</mi></mrow><mrow><mi>λ</mi></mrow></msub><mo>=</mo><msup><mrow><mfenced open="|" close="|" separators="|"><mrow><mi mathvariant="normal">e</mi><mi mathvariant="normal">x</mi><mi mathvariant="normal">p</mi><mfenced separators="|"><mrow><mi>j</mi><mi>k</mi><mfenced separators="|"><mrow><mroot><mrow><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>+</mo><msup><mrow><mi>z</mi></mrow><mrow><mi>'</mi></mrow></msup></mrow><mrow/></mroot><mo>+</mo><mo stretchy="false">(</mo><mi>n</mi><mo>-</mo><mn mathvariant="normal">1</mn><mo stretchy="false">)</mo><mi>ϕ</mi><mo stretchy="false">(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo stretchy="false">)</mo></mrow></mfenced></mrow></mfenced><mi mathvariant="normal">*</mi><mi mathvariant="normal">e</mi><mi mathvariant="normal">x</mi><mi mathvariant="normal">p</mi><mfenced separators="|"><mrow><mi>j</mi><mfrac><mrow><mi>k</mi></mrow><mrow><mn mathvariant="normal">2</mn><mi>z</mi></mrow></mfrac><mfenced separators="|"><mrow><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></mrow></mfenced></mrow></mfenced></mrow></mfenced></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=37126108&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=37126103&type= 111.42134094 10.66800022 ， 其中： <math id="M5"><msub><mrow><mi>p</mi></mrow><mrow><mi>λ</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=37126119&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=37126114&type= 3.38666677 4.57200003 为单透镜光学系统的PSF， <math id="M6"><mi>j</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=37126129&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=37126124&type= 1.60866666 3.64066648 为虚数单位， <math id="M7"><msup><mrow><mi>z</mi></mrow><mrow><mi>'</mi></mrow></msup></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=37126142&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=37126136&type= 2.28600001 3.04800010 表示待模拟目标场景距镜头的距离， <math id="M8"><mi>n</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=37126152&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=37126147&type= 1.86266661 2.79399991 表示透镜材料的折射率， <math id="M9"><mi>ϕ</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=37126162&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=37126158&type= 11.68400002 3.97933316 表示透镜的面型（含可能的加工误差）， <math id="M10"><mi>z</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=37126171&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=37126166&type= 1.77800000 2.79399991 表示透镜距探测器的距离， <math id="M11"><mi>λ</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=37126181&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=37126176&type= 1.77800000 2.79399991 表示目标场景的波长， <math id="M12"><mi>k</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=37126193&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=37126187&type= 1.69333339 2.79399991 表示波数。有了透镜的点扩散函数后，完整的探测器输出计算过程 <math id="M13"><msub><mrow><mi>f</mi></mrow><mrow><mi mathvariant="normal">s</mi><mi mathvariant="normal">e</mi><mi mathvariant="normal">n</mi><mi mathvariant="normal">s</mi><mi mathvariant="normal">o</mi><mi mathvariant="normal">r</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=37126200&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=37126198&type= 7.02733374 4.57200003 可以描述为： <math display="block" id="M14"><msub><mrow><mi>I</mi></mrow><mrow><mi mathvariant="normal">c</mi></mrow></msub><mo stretchy="false">(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo stretchy="false">)</mo><mo>=</mo><mi>S</mi><mo stretchy="false">[</mo><mo>∫</mo><mo stretchy="false">(</mo><msub><mrow><mi>I</mi></mrow><mrow><mi>λ</mi></mrow></msub><mi>*</mi><msub><mrow><mi>p</mi></mrow><mrow><mi>λ</mi></mrow></msub><mo stretchy="false">)</mo><mo stretchy="false">(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo stretchy="false">)</mo><msub><mrow><mi>R</mi></mrow><mrow><mi mathvariant="normal">c</mi></mrow></msub><mo stretchy="false">(</mo><mi>λ</mi><mo stretchy="false">)</mo><mi mathvariant="normal">d</mi><mi>λ</mi><mo stretchy="false">]</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=37126209&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=37126205&type= 60.87533569 6.60400009 ， 其中： <math id="M15"><msub><mrow><mi>R</mi></mrow><mrow><mi mathvariant="normal">c</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=37126220&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=37126216&type= 3.80999994 3.72533321 表示探测器光谱响应， <math id="M16"><msub><mrow><mi>I</mi></mrow><mrow><mi>λ</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=37126228&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=37126224&type= 2.45533323 3.72533321 表示目标场景图像， <math id="M17"><mi>S</mi><mo stretchy="false">[</mo><mo>•</mo><mo stretchy="false">]</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=37126241&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=37126233&type= 7.95866632 3.55599999 表示离散采样以及探测器噪声模拟操作。从 公式（1） 公式（1） 可以看出，点扩散函数的模拟过程为二维运算，其计算量随着光学元件口径以平方关系增加。比如设计一个8 mm口径的单透镜，采样间隔为2 μm，数据量高达16 000 000，这给端到端误差方向传播过程中梯度存储造成了很大的内存需求挑战。我们在2020年提出了基于同心圆环分解的点扩散函数计算降维方法 ［ 6 6 ］ ，如 图2 图2 所示。考虑到光学系统一般都是旋转对称性结构，我们将原本的二维面型表述转化成一系列同心圆孔的线性组合，并且利用圆孔的点扩散函数（已知的一阶贝塞尔函数）成功地将原先二维PSF模拟转化为一维线性组合，使得整个原始的内存需求降低了1个数量级，可以进一步设计更大口径光学系统。 前述基于波动光学的可微分成像模块只能够用来设计单透镜，虽然已经实现了很多功能，但单个元件的调控能力难免不足。2021年，我们发展了针对透镜组的可微分光线追踪成像模型 ［ 22 22 ］ ，如 图3 图3 所示。可微分射线光线追踪广泛应用于图形学中的动画渲染领域，我们将其引入到镜头成像中。具体来说，我们采取了反向追踪的方法，首先将探测器按照像素大小划分成对应数量的网格，然后对每个网格发出一组光线逆向追踪到目标场景空间，根据每根光线与目标场景的交点位置计算每个光线对应的强度，最后把它们相加作为所对应目标像素的图像信号。为了减小每个像素点所需追踪的光线数量，我们采取了多次蒙特卡洛追迹后平均的方法，这使得我们每次仅需追踪64根光线，10次平均即可得到接近百万根光线的效果。这种可微分光线追踪成像模型没有使用诸如徬轴近似的假设，理论上其可以用于任何复杂光学系统和图像处理算法的端到端协同设计。 2.2　图像重建模块 2 图像重建模块用于将探测器输出的原始图像进行重建，常用的重建方法包括基于物理模型的重建方法、基于深度学习的重建方法以及物理模型结合深度学习的重建方法。基于物理模型的重建方法最典型的就是修正维纳滤波 ［ 5 5 ］ ，与传统维纳滤波算法不同，为了更便于在设计初期让深度学习框架给与光学部分更大的权重，我们对原始图像额外进行了加权： <math display="block" id="M18"><msub><mrow><mover><mrow><mi>I</mi></mrow><mrow><mi>˜</mi></mrow></mover></mrow><mrow><mi mathvariant="normal">c</mi></mrow></msub><mo>=</mo><msup><mrow><mi>F</mi></mrow><mrow><mo>-</mo><mn mathvariant="normal">1</mn></mrow></msup><mfenced separators="|"><mrow><mfrac><mrow><mi>F</mi><mo stretchy="false">(</mo><msub><mrow><mi>p</mi></mrow><mrow><mi mathvariant="normal">c</mi></mrow></msub><msup><mrow><mo stretchy="false">)</mo></mrow><mrow><mi mathvariant="normal">*</mi></mrow></msup></mrow><mrow><msup><mrow><mfenced open="|" close="|" separators="|"><mrow><mi>F</mi><mo stretchy="false">(</mo><msub><mrow><mi>p</mi></mrow><mrow><mi mathvariant="normal">c</mi></mrow></msub><mo stretchy="false">)</mo></mrow></mfenced></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup><mo>+</mo><mi>γ</mi></mrow></mfrac><mi>F</mi><mo stretchy="false">(</mo><mo stretchy="false">(</mo><mn mathvariant="normal">1</mn><mo>+</mo><mi>γ</mi><mo stretchy="false">)</mo><msub><mrow><mi>I</mi></mrow><mrow><mi>c</mi></mrow></msub><mo stretchy="false">)</mo></mrow></mfenced></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=37126273&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=37126269&type= 58.50466156 13.54666710 ， 其中： <math id="M19"><mi>γ</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=37126290&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=37126284&type= 1.94733346 3.64066648 为修正系数，属于可优化的变量，随着优化的进行， <math id="M20"><mi>γ</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=37126290&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=37126284&type= 1.94733346 3.64066648 会逐渐减小； <math id="M21"><mi>F</mi><mo stretchy="false">(</mo><msub><mrow><mi>p</mi></mrow><mrow><mi mathvariant="normal">c</mi></mrow></msub><mo stretchy="false">)</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=37126300&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=37126296&type= 9.56733322 4.82600021 为光学传递函数； <math id="M22"><msup><mrow/><mrow><mi mathvariant="normal">*</mi></mrow></msup></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=37126308&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=37126305&type= 1.18533325 3.21733332 表示复共轭； <math id="M23"><msub><mrow><mi>I</mi></mrow><mrow><mi mathvariant="normal">c</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=37126319&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=37126313&type= 2.37066650 3.72533321 为探测器采集的图像； <math id="M24"><mi>F</mi><mo stretchy="false">(</mo><mo>·</mo><mo stretchy="false">)</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=37126327&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=37126323&type= 7.19666624 3.55599999 和 <math id="M25"><msup><mrow><mi>F</mi></mrow><mrow><mo>-</mo><mn mathvariant="normal">1</mn></mrow></msup><mo stretchy="false">(</mo><mo>·</mo><mo stretchy="false">)</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=37126332&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=37126329&type= 10.32933331 3.80999994 分别表示傅里叶变换和傅里叶逆变换。 基于深度学习的重建方法也称为数据驱动的重建方法，典型的网络包括Unet和Resnet，我们在2020年提出了一种Unet和Resnet的结合体——Res-Unet ［ 6 6 ］ ，见 图4 图4 。它既保留了Unet对图像在多个尺度上操作的特点，又可以像Resnet那样通过不断扩展网络来增加重建网络的性能，适合模糊程度一般的图像。近年来人们也发展了物理模型结合深度学习的重建方法。如Monakhova等使用了一个循环的交替ADMM（Alternating Direction Method of Multipliers）网络（Le-ADMM） ［ 23 23 ］ ，他们将ADMM循环中的每一个超参数作为优化变量，通过网络训练学习；Ethan Tseng等在2021年提出了在深度学习特征空间进行维纳滤波的重建方法 ［ 10 10 ］ ，其分为3个过程，首先是一个多尺度的特征提取网络，然后是维纳滤波层，最后是特征重建网络，非常适合中等模糊的图像。Kyrollos Yanny等在2022年提出MultiWienerNet网络，采取了维纳滤波和Unet级联的方案 ［ 24 24 ］ ，由于维纳滤波在最前面，这种方式适合处理非常模糊的图像。 1;2; 2.3　误差模块 2 误差模块用于计算重建图像与目标图像之间偏差，并反向传播用来优化重建网络和光学系统结构参数。一般多用图像质量评价指标作为误差函数。如重建图像与目标图像的L2范数（均方差MSE），或者L1范数。也可以直接使用结构相似性（Structural Similarity， SSIM）或者峰值估噪比（Peak Signal to Noise Ratio， PSNR）作为误差函数。近年来，感知图像评价指标（Perceptual Loss）也被用来作为误差函数，其通常是取VGG网络中某几层输出的特征来计算重建图像与目标图像的接近程度。对抗网络也常用来评估重建图像与目标图像的偏差，例如我们在2021年在复杂镜头端到端设计中就使用了全局和局部对抗网络 ［ 22 22 ］ ，如 图5 图5 。通常选择前述多个的加权组合做为误差函数。 3 光学系统与图形处理算法端到端设计应用实例 1 3.1　大景深成像 2 传统成像镜头的景深与收集光的能力相互制约，计算成像可以打破这种制约，通过设计一个具备随深度近似不变的PSF的光学系统，然后应用深度平移不变的反卷积重建全深度清晰的图像。虽然研究者们对于这个PSF应该是什么样的做了很多研究，但脱离了后端图像处理算法的任何设计总是不尽如人意。2018年，我们使用端到端设计架构分别实现了基于单个自由曲面透镜和衍射光学元件的大景深成像系统，在该方法下，不同深度的PSF到底要设计成什么样可由图像处理算法以不同深度的成像质量为优化目标直接优化得到。仿真性能表明使用端到端协同设计的光学系统与最经典的三次相位板的景深延拓方法提升了6 dB ［ 5 5 ］ 。使用端到端协同设计，不仅成功避免了传统反卷积离面极容易出现的振铃现象，更显著提升了噪声鲁棒性。2021年，我们使用新发展的基于可微分射线追踪的成像模型框架，在3片透镜的系统中也实现了类似的效果 ［ 22 22 ］ ，如 图6 图6 。 3.2　超分辨成像 2 通常成像系统的分辨率由探测器像元大小和光学系统焦距共同决定。为了提升分辨率，要么增加光学系统焦距，要么更换像素单元更小、像元规模更高的探测器，但两种方法都会带来巨大硬件成本。以当前超快光学里面用的单光子雪崩二极管（SPAD）相机为例，目前其像素规模只有数十万量级，远低于当前CMOS相机，迫切需要研究提升其分辨率的方式。而且SPAD相机通常填充率很低，这通常会带来显著的空间频率混叠，从而使得该问题比传统意义上的图像超分辨更加复杂。2019年，我们采取了端到端设计理念，将光学系统、欠采样探测器、超分辨算法作为一个整体进行考虑，通过设计一个光学相位板对图像进行滤波，在抑制欠采样探测器导致混频的同时尽可能地为后续超分辨网络提供更多的信息，成功地实现了4倍超分辨成像 ［ 9 9 ］ 。 图7 图7 为端到端设计的相位板实现SPAD相机4倍超分辨率成像仿真实验表明，采用端到端设计的光学相位板提供的滤波器相比于之前通常用的高斯滤波器在空间分辨率上有显著提升，见 图7 图7 （a）中的图像。而实际成像实验也验证了基于端到端优化的相位板获得了明显的分辨率提升，见 图7 图7 （b）。 3.3　快照式光谱成像 2 快照式光谱因为具备快速响应的优点，是当前光谱成像系统研究的热点方向，但传统的快照式光谱成像系统通常需要很多光学元件，如色散光学元件、光强编码元件、中继镜头和成像镜头，导致快照式光谱成像系统体积都很大。衍射光学元件等新型微结构器件可以实现光谱编码与成像一体化，而且具备轻薄化的特点，可极大的缩减系统体积，但是如何在考虑衍射光学元件实际光谱调控能力的情况下设计光谱编码是难点。我们基于端到端设计理念，协同设计衍射光学元件的光学编码与后续深度学习图像重建算法，实现了单片衍射光学元件的快照式光谱成像系统 ［ 8 8 ］ ，如 图8 图8 。我们将衍射光学元件设计成可对光谱有各向异性响应，这样可以使得在PSF的大小保持不变时，其形状随着波长的变化快速旋转。由于没有了色散元件和光学编码元件，光谱重建的不确定度也显著增加，通过使用基于模型的重建算法专门用于反卷积的光谱图像重建，仿真和实验均表明这种基于衍射元件的快照式光谱成像可以获得与多数快照式光谱成像系统相当的光谱分辨率和空间分辨率。 3.4　基于平面微纳透镜的宽谱成像 2 基于人工微结构材料的微纳光学器件因其平面化、可调能力强等优势，有望颠覆传统复杂光学系统，实现轻薄化、小型化成像系统，目前已成为现在学术界和产业界关注的热点方向之一。但是微纳光学元件，如衍射光学元件（Diffractive Optical Element， DOE） ［ 26 26 - 28 28 26-28 ］ 、超表面 ［ 29 29 - 33 33 29-33 ］ 等都存在极强的内在色差，这严重制约了其在实际宽谱段如彩色成像中的应用。近年来，研究者发现可以通过逆向设计实现消色差微纳光学元件，其原理是将微纳结构的色差转移到波长无关的残留像差，然后再通过后续图像处理算法对波长无关的残留像差进行校正。长期以来，研究者们采取的也都是顺序设计模式，如前面所述，这种模式从计算成像的角度并非最优。我们通过应用端到端设计框架来解决这一问题 ［ 6 6 ］ 。如 图9 图9 所示，采用端到端设计的方式，其图像质量相比顺序设计模式在峰值信噪比（Peak Signal-to-noise Ratio， PSNR）、结构相似性（Structural Similarity， SSIM）和光谱角制图（Spectual Angular Mapper， SAM）分别比提升了1.3 dB、0.015和0.01，这表明端到端设计的消色差平面微纳透镜具有更好的宽谱段成像性能。而且通过端到端设计，我们研究发现，对于计算成像的消色差成像，其设计目标是使得所有波长在高频处的光学传递函数（OTF）都尽可能一致且尽可能高。 3.5　基于平面微纳透镜的大视场成像 2 对于一个成像系统，其核心指标包括F数、视场、焦距等，从目前的发展趋势来看，成像系统的F数越来越小，视场越来越大。为了实现这一需求，势必需要增加系统的复杂度和长度。但人们一直在思考是否有方法能够打破这一限制。结合计算成像的成像系统设计为这一方案提供了可能，但受制于设计理念，以前单镜片的成像视场通常只有10°到20°。我们针对大视场的需求，参考大景深成像采用的深度不变点扩散函数设计目标，设计了具有视场近似不变的点扩散函数，并针对性的设计了图像重构网络，最终实现了视场53°，F数1.8，焦距43 mm的单镜片成像系统 ［ 7 7 ］ 。2021年，我们针对原始算法在暗区表现不佳的问题，又发展了场景相关的图像恢复算法，使得全视场的成像质量又有所提升 ［ 34 34 ］ 。早期我们的端到端设计模型只用于徬轴近似条件，在2021年我们发展的基于可微分射线追踪成像模型打破了徬轴近似的限制，结合端到端设计架构实现了2个镜片的大视场成像。该系统与我们前期没有采用端到端设计的镜头相比，获得了更好的全视场成像质量，特别是在虚假目标抑制方面，大视场成像效果与对比见 图10 图10 ［ 22 22 ］ 。 4 总 结 1 计算成像本质上是联合光学硬件与图像处理算法来协同实现新功能、高性能的成像技术。光学系统与图像处理算法端到端协同设计为各环节深入协同提供了基础，它打破了传统的成像系统中光学、传感器、图像后处理算法以及显示处理等环节之间的壁垒，降低每个环节对经验的依赖，目标是为诸多场景提供傻瓜化的全新解决方案。目前光学成像系统的发展趋势正朝着轻小型化和高成像质量方向发展，例如消费电子和消费光子学要求成像系统的F数小、视场大，又要求成像系统重量轻、体积小。光学系统与图像处理算法端到端协同设计提高了这类光学系统设计的可能性，目前在平面透镜的宽谱成像、平面透镜的大视场成像、大景深成像，超分辨成像，快照式光谱成像方面取得了不错的进展。我们正在推进该技术在红外光电对抗领域的应用。 端到端设计架构在单镜片以及复杂镜片设计中都进行了有益的尝试并取得了不错的结果，如可微复合透镜光学模型，不仅克服了单一光学表面和单一材料的限制，而且支持优化离轴区域。但这种基于可微分射线追踪成像目前还是过于复杂，在优化复杂镜头时难以收敛，如何在保证精度的情况下对其进行降维是未来的研究方向之一，比如是否可以使用深度学习来进行透镜PSF的模拟，即只要给定结构参数，可快速通过一个网络计算出其所对应的各视场的点扩散函数。此外计算成像系统应用推进所面临的的计算资源需求过大问题也亟需解决，如何将前端光学系统、探测器采样、后端算法甚至后端算法可用硬件资源一起协同也是重要的研究方向。 另外，未来的光学成像系统要想真正实现轻薄化，不仅需要在设计方法上有突破，还需要在新材料上有突破，目前超表面这种小于波长的超薄人工结构是最有可能突破现有光学系统设计瓶颈的材料。虽然超表面有很大潜力，但其内在更大的自由度和自身固有的结构色散导致目前的成像质量和传统折射型光学系统还存在一定的差距。美国普林斯顿大学Heide等人使用端到端架构设计超表面实现拥有全色覆盖、宽视场成像的超表面相机 ［ 10 10 ］ 。在未来，端到端架构和超表面这种具有强大调控能力的人工材料相结合可能会碰撞出新的火花，实现更加强大的光学成像功能。 参考文献 左超 ， 陈钱 . 计算光学成像：何来，何处，何去，何从？ ［J］. 红外与激光工程 ， 2022 ， 51 （ 2 ）： 158 - 338 . doi: 10.3788/IRLA20220110 http://dx.doi.org/10.3788/IRLA20220110 计算光学成像：何来，何处，何去，何从？ ZUO CH ， CHEN Q . Computational optical imaging： an overview ［J］. Infrared and Laser Engineering ， 2022 ， 51 （ 2 ）： 158 - 338 . （in Chinese） . doi: 10.3788/IRLA20220110 http://dx.doi.org/10.3788/IRLA20220110 Computational optical imaging： an overview LEVOY M . Light fields and computational imaging ［J］. Computer ， 2006 ， 39 （ 8 ）： 46 - 55 . doi: 10.1109/mc.2006.270 http://dx.doi.org/10.1109/mc.2006.270 Light fields and computational imaging SUNG H Y ， CHANG C W . Method for designing computational optical imaging system ： US8331712 ［P］. 2012-12-11 . Method for designing computational optical imaging system MAIT J N ， EULISS G W ， ATHALE R A . Computational imaging ［J］. Advances in Optics and Photonics ， 2018 ， 10 （ 2 ）： 409 - 483 . doi: 10.1364/aop.10.000409 http://dx.doi.org/10.1364/aop.10.000409 Computational imaging SITZMANN V ， DIAMOND S ， PENG Y F ， et al . End-to-end optimization of optics and image processing for achromatic extended depth of field and super-resolution imaging ［J］. ACM Transactions on Graphics ， 2018 ， 37 （ 4 ）： 114 . doi: 10.1145/3197517.3201333 http://dx.doi.org/10.1145/3197517.3201333 End-to-end optimization of optics and image processing for achromatic extended depth of field and super-resolution imaging DUN X ， IKOMA H ， WETZSTEIN G ， et al . Learned rotationally symmetric diffractive achromat for full-spectrum computational imaging ［J］. Optica ， 2020 ， 7 （ 8 ）： 913 - 922 . doi: 10.1364/optica.394413 http://dx.doi.org/10.1364/optica.394413 Learned rotationally symmetric diffractive achromat for full-spectrum computational imaging PENG Y F ， SUN Q L ， DUN X ， et al . Learned large field-of-view imaging with thin-plate optics ［J］. ACM Transactions on Graphics ， 2019 ， 38 （ 6 ）： 219 . doi: 10.1145/3355089.3356526 http://dx.doi.org/10.1145/3355089.3356526 Learned large field-of-view imaging with thin-plate optics JEON D S ， BAEK S H ， YI S ， et al . Compact snapshot hyperspectral imaging with diffracted rotation ［J］. ACM Transactions on Graphics ， 2019 ， 38 （ 4 ）： 117 . doi: 10.1145/3306346.3322946 http://dx.doi.org/10.1145/3306346.3322946 Compact snapshot hyperspectral imaging with diffracted rotation SUN Q L ， ZHANG J ， DUN X ， et al . End-to-end learned， optically coded super-resolution SPAD camera ［J］. ACM Transactions on Graphics ， 2020 ， 39 （ 2 ）： 9 . doi: 10.1145/3372261 http://dx.doi.org/10.1145/3372261 End-to-end learned， optically coded super-resolution SPAD camera TSENG E ， COLBURN S ， WHITEHEAD J ， et al . Neural nano-optics for high-quality thin lens imaging ［J］. Nature Communications ， 2021 ， 12 ： 6493 . doi: 10.1038/s41467-021-26443-0 http://dx.doi.org/10.1038/s41467-021-26443-0 Neural nano-optics for high-quality thin lens imaging LÉVÊQUE O ， KULCSÁR C ， LEE A ， et al . Co-designed annular binary phase masks for depth-of-field extension in single-molecule localization microscopy ［J］. Optics Express ， 2020 ， 28 （ 22 ）： 32426 - 32446 . doi: 10.1364/oe.402752 http://dx.doi.org/10.1364/oe.402752 Co-designed annular binary phase masks for depth-of-field extension in single-molecule localization microscopy JIN L B ， TANG Y B ， WU Y C ， et al . Deep learning extended depth-of-field microscope for fast and slide-free histology ［J］. Proceedings of the National Academy of Sciences of the United States of America ， 2020 ， 117 （ 52 ）： 33051 - 33060 . doi: 10.1073/pnas.2013571117 http://dx.doi.org/10.1073/pnas.2013571117 Deep learning extended depth-of-field microscope for fast and slide-free histology LIU Y K ， ZHANG C Y ， KOU T D ， et al . End-to-end computational optics with a singlet lens for large depth-of-field imaging ［J］. Optics Express ， 2021 ， 29 （ 18 ）： 28530 - 28548 . doi: 10.1364/oe.433067 http://dx.doi.org/10.1364/oe.433067 End-to-end computational optics with a singlet lens for large depth-of-field imaging BANERJI S ， MEEM M ， MAJUMDER A ， et al . Extreme-depth-of-focus imaging with a flat lens ［J］. Optica ， 2020 ， 7 （ 3 ）： 214 - 217 . doi: 10.1364/optica.384164 http://dx.doi.org/10.1364/optica.384164 Extreme-depth-of-focus imaging with a flat lens LUO Y ， MENGU D ， YARDIMCI N T ， et al . Design of task-specific optical systems using broadband diffractive neural networks ［J］. Light： Science & Applications ， 2019 ， 8 ： 112 . doi: 10.1038/s41377-019-0223-1 http://dx.doi.org/10.1038/s41377-019-0223-1 Design of task-specific optical systems using broadband diffractive neural networks SONG H Y ， MA Y G ， HAN Y B ， et al . Deep-learned broadband encoding stochastic filters for computational spectroscopic instruments ［J］. Advanced Theory and Simulations ， 2021 ， 4 （ 3 ）： 2000299 . doi: 10.1002/adts.202000299 http://dx.doi.org/10.1002/adts.202000299 Deep-learned broadband encoding stochastic filters for computational spectroscopic instruments CHANG J L ， WETZSTEIN G . Deep optics for monocular depth estimation and 3D object detection ［C］. 2019 IEEE/CVF International Conference on Computer Vision （ICCV）. Seoul ， Korea （South） . IEEE ， 2019 ： 10192 - 10201 . doi: 10.1109/iccv.2019.01029 http://dx.doi.org/10.1109/iccv.2019.01029 Deep optics for monocular depth estimation and 3D object detection NEHME E ， FREEDMAN D ， GORDON R ， et al . DeepSTORM3D： dense 3D localization microscopy and PSF design by deep learning ［J］. Nature Methods ， 2020 ， 17 （ 7 ）： 734 - 740 . doi: 10.1038/s41592-020-0853-5 http://dx.doi.org/10.1038/s41592-020-0853-5 DeepSTORM3D： dense 3D localization microscopy and PSF design by deep learning NEHME E ， FERDMAN B ， WEISS L E ， et al . Learning an optimal PSF-pair for ultra-dense 3D localization microscopy ［EB/OL］. 2020： arXiv ： 2009 . 14303 . https：//arxiv.org/abs/2009.14303 https://arxiv.org/abs/2009.14303 . doi: 10.1038/s41592-020-0853-5 http://dx.doi.org/10.1038/s41592-020-0853-5 Learning an optimal PSF-pair for ultra-dense 3D localization microscopy SUN Q L ， TSENG E ， FU Q ， et al . Learning rank-1 diffractive optics for single-shot high dynamic range imaging ［C］. 2020 IEEE/CVF Conference on Computer Vision and Pattern Recognition （CVPR）. Seattle ， WA ， USA . IEEE ， 2020 ： 1383 - 1393 . doi: 10.1109/cvpr42600.2020.00146 http://dx.doi.org/10.1109/cvpr42600.2020.00146 Learning rank-1 diffractive optics for single-shot high dynamic range imaging MARTEL J N P ， MÜLLER L K ， CAREY S J ， et al . Neural sensors： learning pixel exposures for HDR imaging and video compressive sensing with programmable sensors ［J］. IEEE Transactions on Pattern Analysis and Machine Intelligence ， 2020 ， 42 （ 7 ）： 1642 - 1653 . doi: 10.1109/tpami.2020.2986944 http://dx.doi.org/10.1109/tpami.2020.2986944 Neural sensors： learning pixel exposures for HDR imaging and video compressive sensing with programmable sensors SUN Q L ， WANG C L ， FU Q ， et al . End-to-end complex lens design with differentiate ray tracing ［J］. ACM Transactions on Graphics ， 2021 ， 40 （ 4 ）： 71 . doi: 10.1145/3450626.3459674 http://dx.doi.org/10.1145/3450626.3459674 End-to-end complex lens design with differentiate ray tracing VOLLMER M ， MÖLLMANN K P . Infrared Thermal Imaging： Fundamentals， Research and Applications ［M］. Weinheim， Germany ： Wiley-VCH Verlag GmbH & Co. KGaA ， 2017 . doi: 10.1002/9783527693306 http://dx.doi.org/10.1002/9783527693306 Infrared Thermal Imaging： Fundamentals， Research and Applications YANNY K ， MONAKHOVA K ， SHUAI R W ， et al . Deep learning for fast spatially varying deconvolution ［J］. Optica ， 2022 ， 9 （ 1 ）： 96 - 99 . doi: 10.1364/optica.442438 http://dx.doi.org/10.1364/optica.442438 Deep learning for fast spatially varying deconvolution MONAKHOVA K ， YURTSEVER J ， KUO G ， et al . Learned reconstructions for practical mask-based lensless imaging ［J］. Optics Express ， 2019 ， 27 （ 20 ）： 28075 - 28090 . doi: 10.1364/oe.27.028075 http://dx.doi.org/10.1364/oe.27.028075 Learned reconstructions for practical mask-based lensless imaging BANERJI S ， MEEM M ， MAJUMDER A ， et al . Imaging with flat optics： metalenses or diffractive lenses？ ［J］. Optica ， 2019 ， 6 （ 6 ）： 805 - 810 . doi: 10.1364/optica.6.000805 http://dx.doi.org/10.1364/optica.6.000805 Imaging with flat optics： metalenses or diffractive lenses？ MEEM M ， BANERJI S ， PIES C ， et al . Large-area， high-numerical-aperture multi-level diffractive lens via inverse design： errata ［J］. Optica ， 2020 ， 7 （ 3 ）： 252 - 253 . doi: 10.1364/optica.388697 http://dx.doi.org/10.1364/optica.388697 Large-area， high-numerical-aperture multi-level diffractive lens via inverse design： errata WANG P ， MOHAMMAD N ， MENON R . Chromatic-aberration-corrected diffractive lenses for ultra-broadband focusing ［J］. Scientific Reports ， 2016 ， 6 ： 21545 . doi: 10.1038/srep21545 http://dx.doi.org/10.1038/srep21545 Chromatic-aberration-corrected diffractive lenses for ultra-broadband focusing TSENG M L ， HSIAO H ， CHU C H ， et al . Metalenses： advances and applications ［J］. Advanced Optical Materials ， 2018 ， 6 （ 18 ）： 1800554 . doi: 10.1002/adom.201800554 http://dx.doi.org/10.1002/adom.201800554 Metalenses： advances and applications AIETA F ， KATS M A ， GENEVET P ， et al . Applied optics. Multiwavelength achromatic metasurfaces by dispersive phase compensation ［J］. Science ， 2015 ， 347 （ 6228 ）： 1342 - 1345 . doi: 10.1126/science.aaa2494 http://dx.doi.org/10.1126/science.aaa2494 Applied optics. Multiwavelength achromatic metasurfaces by dispersive phase compensation KHORASANINEJAD M ， CHEN W T ， DEVLIN R C ， et al . Metalenses at visible wavelengths： Diffraction-limited focusing and subwavelength resolution imaging ［J］. Science ， 2016 ， 352 （ 6290 ）： 1190 - 1194 . doi: 10.1126/science.aaf6644 http://dx.doi.org/10.1126/science.aaf6644 Metalenses at visible wavelengths： Diffraction-limited focusing and subwavelength resolution imaging KHORASANINEJAD M ， ZHU A Y ， ROQUES-CARMES C ， et al . Polarization-insensitive metalenses at visible wavelengths ［J］. Nano Letters ， 2016 ， 16 （ 11 ）： 7229 - 7234 . doi: 10.1021/acs.nanolett.6b03626 http://dx.doi.org/10.1021/acs.nanolett.6b03626 Polarization-insensitive metalenses at visible wavelengths PHAN T ， SELL D ， WANG E W ， et al . High-efficiency， large-area， topology-optimized metasurfaces ［J］. Light： Science & Applications ， 2019 ， 8 ： 48 . doi: 10.1038/s41377-019-0159-5 http://dx.doi.org/10.1038/s41377-019-0159-5 High-efficiency， large-area， topology-optimized metasurfaces QI B Y ， CHEN W ， DUN X ， et al . All-day thin-lens computational imaging with scene-specific learning recovery ［J］. Applied Optics ， 2022 ， 61 （ 4 ）： 1097 - 1105 . doi: 10.1364/ao.448155 http://dx.doi.org/10.1364/ao.448155 All-day thin-lens computational imaging with scene-specific learning recovery 10.37188/OPE.20223021.2827 O439 TP751 A 1004-924X（2022）21-2827-12 顿雄,张健,冯诗淇等.光学系统与图像处理端到端协同设计及其应用[J].光学精密工程,2022,30(21):2827-2838. 顿雄,张健,冯诗淇等.光学系统与图像处理端到端协同设计及其应用[J].光学精密工程,2022,30(21):2827-2838. DUN Xiong,ZHANG Jian,FENG Shiqi,et al.End-to-end co-design of optics and image processing and its applications[J].Optics and Precision Engineering,2022,30(21):2827-2838. DUN Xiong,ZHANG Jian,FENG Shiqi,et al.End-to-end co-design of optics and image processing and its applications[J].Optics and Precision Engineering,2022,30(21):2827-2838. 上海市浦江人才计划资助项目（No.20PJ1414200）；国家自然科学基金资助项目（No.62105243） 2022-07-16 2022-08-20 2022-11-10 2022-11-22 光学精密工程 21 30