CN115290601A - 一种宽谱非相干光散斑自相关成像探测的低冗余模拟方法 - Google Patents

一种宽谱非相干光散斑自相关成像探测的低冗余模拟方法 Download PDF

Info

Publication number
CN115290601A
CN115290601A CN202211226822.2A CN202211226822A CN115290601A CN 115290601 A CN115290601 A CN 115290601A CN 202211226822 A CN202211226822 A CN 202211226822A CN 115290601 A CN115290601 A CN 115290601A
Authority
CN
China
Prior art keywords
wavelength
light
wide
spectrum
detection
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN202211226822.2A
Other languages
English (en)
Other versions
CN115290601B (zh
Inventor
王歆
施钧辉
梁其传
王钰琪
陈睿黾
李驰野
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Zhejiang Lab
Original Assignee
Zhejiang Lab
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Zhejiang Lab filed Critical Zhejiang Lab
Priority to CN202211226822.2A priority Critical patent/CN115290601B/zh
Publication of CN115290601A publication Critical patent/CN115290601A/zh
Application granted granted Critical
Publication of CN115290601B publication Critical patent/CN115290601B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N21/00—Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
    • G01N21/17—Systems in which incident light is modified in accordance with the properties of the material investigated
    • G01N21/47—Scattering, i.e. diffuse reflection
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01J—MEASUREMENT OF INTENSITY, VELOCITY, SPECTRAL CONTENT, POLARISATION, PHASE OR PULSE CHARACTERISTICS OF INFRARED, VISIBLE OR ULTRAVIOLET LIGHT; COLORIMETRY; RADIATION PYROMETRY
    • G01J1/00—Photometry, e.g. photographic exposure meter

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Spectroscopy & Molecular Physics (AREA)
  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Chemical & Material Sciences (AREA)
  • Analytical Chemistry (AREA)
  • Biochemistry (AREA)
  • General Health & Medical Sciences (AREA)
  • Immunology (AREA)
  • Pathology (AREA)
  • Investigating Or Analysing Materials By Optical Means (AREA)

Abstract

本发明公开了一种宽谱非相干光散斑自相关成像探测的低冗余模拟方法。包括以下步骤：设置系统参数；选择物面中心点作为标定目标点；生成单一波长下的点扩展函数；将宽谱光源离散成多个单色波长；选出点扩展函数间相关性较小的波长作为模拟中的有效波长；生成各个有效波长对应的物体光场非相干传播到探测面所产生的光强分布；将各个光强分布按照一定比例加权平均得到物体在宽谱照明下在探测面产生的光强分布。本发明可以极大地降低模拟中的冗余，减小计算量，提高模拟效率；可以应用在非相干光的散射成像模拟实验中，为实际应用提供参考，在水下探测，透雾成像等方面都有很大应用前景。

Description

一种宽谱非相干光散斑自相关成像探测的低冗余模拟方法
技术领域
本申请涉及非相干光成像模拟领域，尤其涉及一种宽谱非相干光散斑自相关成像探测的低冗余模拟方法。
背景技术
透过散射介质的新型成像技术可以有效解决光透过雾、霾、云或生物组织等散射介质时，出射光场紊乱所导致的传统成像图模糊甚至无效的问题，具有很多优势。其中利用非相干光的散射成像技术更接近自然界中的被动照明场景，应用潜力巨大。仅散斑相关技术领域，就涌现出许多重要的研究成果。自Bertolotti团队在论文《J. Bertolotti, etal., “Non-invasive imaging through opaque scattering layers,” Nature 491,232–234 (2012).》展示了其基于散射介质的光学记忆效应，利用强度相关性重建了隐藏于散射介质后的荧光目标开始，先后有多个团队从成像速度，成像功能等方面对其做了进一步提高。比如，Katz团队在论文《O. Katz, et al., “Non-invasive single-shot imagingthrough scattering layers and around corners via speckle correlations,” Nat.Photon. 8, 784–790 (2014).》中率先提出使用单帧散斑来实现隐藏目标的重建，极大缩减了采样时间，Cua等人所在团队在论文《M. Cua, et al., “Imaging moving targetsthrough scattering Media,” Opt. Express 25, 3935 (2017).》中展示了一种利用背景不变性完成对明亮场景中移动目标的图案重建方法，西安电子科技大学邵晓鹏课题组则在论文《C. Guo, et al., “Tracking moving targets behind a scattering medium viaspeckle correlation,” Appl. Optics 57, 905-913 (2018).》中展示了散斑相关对目标旋转等运动状态的追踪能力。尽管这一成像技术在逐步完善，但是大部分针对宽谱非相干光的成像研究都以实验为主，虽然部分研究有针对单色的非相干赝热光照明进行成像模拟，但这在散斑对比度和有效采样面积等方面都与实际的宽谱非相干光照明有很大区别。所以，亟需提出一种全新的宽谱非相干光散射成像模拟方法，以从根本上做好理论与实验之间的衔接，更好的分析各个过程对应的参数阈值，了解应用极限，获得更详实的研究成果。
发明内容
针对现有技术的不足，本申请实施例的目的是提供一种宽谱非相干光散斑自相关成像探测的低冗余模拟方法，以波长之间的相关性作为筛选条件，只对相关性较低的各个波长分步执行成像仿真操作，以在较少重复步骤下实现有效的宽谱光学模拟。
根据本申请实施例的第一方面，提供一种宽谱非相干光散斑自相关成像探测的低冗余模拟方法，应用于宽谱非相干光散射成像探测模拟系统，包括：
步骤（1）：设置宽谱非相干光散射成像探测模拟系统的系统参数；
步骤（2）：选择物面中心点作为标定目标点；
步骤（3）：将波长范围为
Figure 998623DEST_PATH_IMAGE001
的宽谱光源按照预定波长间距对波长进行均匀的 离散性操作，并生成所述标定目标点对应的各个波长的点扩展函数；
步骤（4）：从
Figure 322288DEST_PATH_IMAGE002
开始，选出点扩展函数间相关性低于相关性阈值的波长作为采样模 拟中的有效波长，并始终保留最大波长
Figure 633184DEST_PATH_IMAGE003
作为有效波长；
步骤（5）：生成各个有效波长对应的物体光场非相干传播到所述探测面所产生的光强分布并按照各自的预定权重进行加权平均，得到物体在宽谱照明下在探测面产生的光强分布。
进一步地，步骤（1）中的所述系统参数包括物体形状，物面阵列尺寸、物面最小单元尺寸，物体到散射介质的距离u，散射介质面阵列尺寸、最小单元尺寸，散射介质特征函数，光阑直径，散射介质到探测面的距离v，探测面阵列大小、探测面像素大小，波长范围参数。
进一步地，步骤（3）包括：
对波长范围为
Figure 434918DEST_PATH_IMAGE001
的宽谱光源，按照波长间距
Figure 780448DEST_PATH_IMAGE004
对波长进行均匀的离散性操 作，得到波长为
Figure 743856DEST_PATH_IMAGE002
+
Figure 542048DEST_PATH_IMAGE005
，...，
Figure 147473DEST_PATH_IMAGE003
的一系列单色波长，其中，
Figure 347510DEST_PATH_IMAGE004
在0.2nm- 2nm之间；
生成波长为
Figure 744469DEST_PATH_IMAGE002
+
Figure 29957DEST_PATH_IMAGE005
，...，
Figure 32548DEST_PATH_IMAGE003
的一系列单色波长中各个波长 对应的点扩展函数。
进一步地，生成点扩展函数的过程包括：
将标定目标点按照菲涅尔衍射的方式传播距离u得到散射介质前表面的入射光场
Figure 962458DEST_PATH_IMAGE006
；
将入射光场
Figure 923460DEST_PATH_IMAGE006
先后和散射介质特征函数以及光阑特征函数相乘，得到散射介质 后表面的出射光场
Figure 571611DEST_PATH_IMAGE007
；
对光场
Figure 112313DEST_PATH_IMAGE007
按照菲涅尔衍射的方式光场传输距离v，得到探测面的光场分布
Figure 162309DEST_PATH_IMAGE008
，将 光场
Figure 28634DEST_PATH_IMAGE008
与其共轭光场
Figure 164080DEST_PATH_IMAGE009
相乘，得到单一波长下位于所述标定目标点的点光源在探测面处 所形成的光强度分布，即该波长下的点扩展函数。
进一步地，步骤（4）包括：
将
Figure 774053DEST_PATH_IMAGE002
作为第一个有效波长，将波长
Figure 537610DEST_PATH_IMAGE002
+
Figure 981360DEST_PATH_IMAGE004
，
Figure 197578DEST_PATH_IMAGE002
+
Figure 486608DEST_PATH_IMAGE010
...
Figure 370250DEST_PATH_IMAGE011
对应的点扩展函数依次 与
Figure 453744DEST_PATH_IMAGE002
对应的点扩展函数做相关性计算，选择相关性初次低于相关性阈值的波长作为第二个 有效波长
Figure 422837DEST_PATH_IMAGE012
；将波长
Figure 247049DEST_PATH_IMAGE013
+
Figure 250777DEST_PATH_IMAGE004
，
Figure 770751DEST_PATH_IMAGE013
+
Figure 227140DEST_PATH_IMAGE010
...，
Figure 717028DEST_PATH_IMAGE003
对应的点扩展函数依次与
Figure 185049DEST_PATH_IMAGE013
对应的点扩展函 数做相关性计算，选择相关性初次低于相关性阈值的波长作为第三个有效波长
Figure 558DEST_PATH_IMAGE014
；以此类 推，得到有效波长
Figure 819610DEST_PATH_IMAGE002
，
Figure 113188DEST_PATH_IMAGE013
，
Figure 701295DEST_PATH_IMAGE014
，...，
Figure 422127DEST_PATH_IMAGE015
；最大波长
Figure 587529DEST_PATH_IMAGE003
即使不满足所述相关性低于相关性阈值 的条件也始终保留作为有效波长；
其中，所述相关性阈值不大于1。
进一步地，步骤（5）包括：
生成物体矩阵，针对有效波长
Figure 560164DEST_PATH_IMAGE016
，对物体矩阵的每一点都生成对应的单一波长
Figure 392991DEST_PATH_IMAGE016
下探测面的光强分布，并将所得各个光强分布叠加，得到物体在该波长下对应的光强分布
Figure 160090DEST_PATH_IMAGE017
，从而得到有效波长
Figure 78367DEST_PATH_IMAGE002
，
Figure 589114DEST_PATH_IMAGE013
，
Figure 10868DEST_PATH_IMAGE014
，...，
Figure 211518DEST_PATH_IMAGE003
各自对应的物体光场非相干传播到探测面所产生 的光强分布
Figure 351512DEST_PATH_IMAGE018
，
Figure 931529DEST_PATH_IMAGE019
，
Figure 207790DEST_PATH_IMAGE020
，...，
Figure 706904DEST_PATH_IMAGE021
，其中k=1，2，…，n；
将各个有效波长下物体光场非相干传播到探测面所产生的光强分布
Figure 209561DEST_PATH_IMAGE018
，
Figure 452323DEST_PATH_IMAGE019
，
Figure 724036DEST_PATH_IMAGE020
，...，
Figure 128472DEST_PATH_IMAGE021
按照预定比例加权平均，得到物体在宽谱照明下在探测面产生的光强分布。
本申请的实施例提供的技术方案可以包括以下有益效果：
由上述实施例可知，本申请利用波长之间的相关性作为筛选条件，只对相关性较低的各个波长分步执行成像仿真操作，可以在较少重复步骤下实现有效的宽谱光学模拟，极大地降低模拟中的冗余，减小计算量，提高模拟效率；可以应用在散斑自相关等一系列使用非相干光的散射成像模拟实验中，将其与相位恢复算法等图像重建方式相结合，能够为实际应用提供参考，在水下探测，透雾成像等方面都有很大应用前景。
应当理解的是，以上的一般描述和后文的细节描述仅是示例性和解释性的，并不能限制本申请。
附图说明
此处的附图被并入说明书中并构成本说明书的一部分，示出了符合本申请的实施例，并与说明书一起用于解释本申请的原理。
图1是根据一示例性实施例示出的一种宽谱非相干光散斑自相关成像探测的低冗余模拟方法的流程图。
图2是根据一示例性实施例示出的不同波长对应的点扩展函数之间的相关性及变化规律示意图。
图3是根据一示例性实施例示出的物体在宽谱照明下在探测面产生的光强分布示意图。
图4是根据一示例性实施例示出的宽谱非相干光散射成像探测模拟系统对应的实际光路传输及探测装置的结构示意图。
附图标记：
1、光源；2、散射介质；3、光阑；4、面阵探测器；5、计算机；a、物面。
具体实施方式
这里将详细地对示例性实施例进行说明，其示例表示在附图中。下面的描述涉及附图时，除非另有表示，不同附图中的相同数字表示相同或相似的要素。以下示例性实施例中所描述的实施方式并不代表与本申请相一致的所有实施方式。
在本申请使用的术语是仅仅出于描述特定实施例的目的，而非旨在限制本申请。在本申请和所附权利要求书中所使用的单数形式的“一种”、“所述”和“该”也旨在包括多数形式，除非上下文清楚地表示其他含义。还应当理解，本文中使用的术语“和/或”是指并包含一个或多个相关联的列出项目的任何或所有可能组合。
图1是根据一示例性实施例示出的一种宽谱非相干光散斑自相关成像探测的低冗余模拟方法的流程图，如图1所示，该方法应用于宽谱非相干光散射成像探测模拟系统中，可以包括以下步骤：
步骤（1）：设置宽谱非相干光散射成像探测模拟系统的系统参数；
步骤（2）：选择物面中心点作为标定目标点；
步骤（3）：将波长范围为
Figure 118425DEST_PATH_IMAGE001
的宽谱光源按照预定波长间距对波长进行均匀的 离散性操作，并生成所述标定目标点对应的各个波长的点扩展函数；
步骤（4）：从
Figure 899299DEST_PATH_IMAGE002
开始，选出点扩展函数间相关性低于相关性阈值的波长作为采样模 拟中的有效波长，并始终保留最大波长
Figure 150152DEST_PATH_IMAGE003
作为有效波长；
步骤（5）：生成各个有效波长对应的物体光场非相干传播到所述探测面所产生的光强分布并按照各自的预定权重进行加权平均，得到物体在宽谱照明下在探测面产生的光强分布。
由上述实施例可知，本申请利用波长之间的相关性作为筛选条件，只对相关性较低的各个波长分步执行成像仿真操作，可以在较少重复步骤下实现有效的宽谱光学模拟，极大地降低模拟中的冗余，减小计算量，提高模拟效率；可以应用在散斑自相关等一系列使用非相干光的散射成像模拟实验中，将其与相位恢复算法等图像重建方式相结合，能够为实际应用提供参考，在水下探测，透雾成像等方面都有很大应用前景。
在步骤（1）的具体实施中，所述系统参数包括物体形状，物面阵列尺寸、物面最小单元尺寸，物体到散射介质的距离u，散射介质面阵列尺寸、最小单元尺寸，散射介质特征函数，光阑直径，散射介质到探测面的距离v，探测面阵列大小、探测面像素大小，波长范围参数；
在步骤（2）的具体实施中，模拟中一个有强度分布的物体用矩阵表示，若物体的矩阵表征为n*n，那么物面中心点就是这个n*n的矩阵的中心。
不同于相干光，非相干光的相位是无规则的，传输模型也更加复杂，对于具有一定带宽的非相干光而言，模拟需要涵盖其各个波长。非相干光的光场之间不发生干涉，光强具有叠加特性，一物体被宽谱非相干光照明后在探测器上所形成的光强分布相当于该物体上的各个点被各个波长的光照明后在探测器上所形成的光强分布的叠加。虽然不同波长对应的光强分布不同，但是波长相近的光形成的光强分布是高度相关的，相关性随着波长差的增大逐渐降低，甚至完全不相关。因此，在步骤（3）的具体实施中，可以包括以下子步骤：
（3.1）对波长范围为
Figure 132014DEST_PATH_IMAGE001
的宽谱光源，按照波长间距
Figure 468318DEST_PATH_IMAGE004
对波长进行均匀的离散 性操作，得到波长为
Figure 928249DEST_PATH_IMAGE002
+
Figure 299188DEST_PATH_IMAGE005
，...，
Figure 920793DEST_PATH_IMAGE003
的一系列单色波长，其中，
Figure 744392DEST_PATH_IMAGE004
在 0.2nm-2nm之间；
具体地，对宽谱光源进行均匀的离散性操作从而得到单色波长为本领域的常规技 术手段，此处不作赘述，为避免太过多余的计算，同时保留选择的合理与严谨性，建议
Figure 5085DEST_PATH_IMAGE004
的 选择范围是0.2nm-2nm，在该范围内，相邻波长依然还保持比较高的相关性，相关性变化对 应一个相对比较平滑的曲线。
（3.2）生成波长为
Figure 230530DEST_PATH_IMAGE002
+
Figure 147670DEST_PATH_IMAGE005
，...，
Figure 599511DEST_PATH_IMAGE003
的一系列单色波长中各个 波长对应的点扩展函数。
生成点扩展函数的过程包括：
（3.2.1）将标定目标点按照菲涅尔衍射的方式传播距离u得到散射介质前表面的 入射光场
Figure 525879DEST_PATH_IMAGE006
；
（3.2.2）将入射光场
Figure 481197DEST_PATH_IMAGE006
先后和散射介质特征函数以及光阑特征函数相乘，得到散 射介质后表面的出射光场
Figure 834818DEST_PATH_IMAGE007
；
（3.2.3）对光场
Figure 508376DEST_PATH_IMAGE007
按照菲涅尔衍射的方式光场传输距离v，得到探测面的光场分 布
Figure 238434DEST_PATH_IMAGE008
，将光场
Figure 48258DEST_PATH_IMAGE008
与其共轭光场
Figure 838360DEST_PATH_IMAGE009
相乘，得到单一波长下位于所述标定目标点的点光源在探 测面处所形成的光强度分布，即该波长下的点扩展函数。
具体地，步骤（3.2.1）-步骤（3.2.3）所描述的生成点扩展函数的过程为本领域的常规技术手段，此处不作赘述。需要说明的是，散射介质特征函数是本领域的常规技术手段，本领域技术人员可根据实际情况自行设置对应的散射介质特征函数以进行不同散射介质的模拟。
在步骤（4）的具体实施中，将
Figure 858269DEST_PATH_IMAGE002
作为第一个有效波长，将波长
Figure 1805DEST_PATH_IMAGE002
+
Figure 56349DEST_PATH_IMAGE004
，
Figure 892718DEST_PATH_IMAGE002
+
Figure 399922DEST_PATH_IMAGE010
...
Figure 81571DEST_PATH_IMAGE011
对应的点扩展函数依次与
Figure 256200DEST_PATH_IMAGE002
对应的点扩展函数做相关性计算，选择相关性初次低于相 关性阈值的波长作为第二个有效波长
Figure 994961DEST_PATH_IMAGE013
；将波长
Figure 255041DEST_PATH_IMAGE013
+
Figure 740381DEST_PATH_IMAGE004
，
Figure 503937DEST_PATH_IMAGE013
+
Figure 541163DEST_PATH_IMAGE010
...，
Figure 163906DEST_PATH_IMAGE003
对应的点扩展函数 依次与
Figure 311990DEST_PATH_IMAGE013
对应的点扩展函数做相关性计算，选择相关性初次低于相关性阈值的波长作为第 三个有效波长
Figure 336578DEST_PATH_IMAGE014
；
以此类推，得到有效波长
Figure 544706DEST_PATH_IMAGE002
，
Figure 389165DEST_PATH_IMAGE013
，
Figure 340940DEST_PATH_IMAGE014
，...，
Figure 220035DEST_PATH_IMAGE015
；最大波长
Figure 599063DEST_PATH_IMAGE003
即使不满足所述相关性 低于相关性阈值的条件也始终保留作为有效波长；
其中，所述相关性阈值不大于1，一般选择0.9即可。图2展示了参考波长
Figure 930819DEST_PATH_IMAGE022
分别450 nm、550 nm、650 nm，
Figure 420706DEST_PATH_IMAGE004
为2 nm时，波长
Figure 278940DEST_PATH_IMAGE022
对应的点扩展函数与
Figure 969816DEST_PATH_IMAGE023
对应的点扩展函数之间的 相关系数变化曲线，来展示不同波长之间的相关性变化规律。图中横坐标为波长差
Figure 913501DEST_PATH_IMAGE024
， 单位nm，纵坐标为相关系数。可见，不同波长对应的点扩展函数之间的相关性随波长差的增 大逐渐降低，并且波长越长相关性下降越慢，这也从侧面说明同样的相关性在不同波段对 应不同的波长间距，相比于固定波长间距下将宽谱光离散为多个单一波长的操作，统一相 关性阈值筛选有效单一波长的做法可以减少数据冗余，帮助执行更高效的模拟。
在步骤（5）的具体实施中，可以包括以下子步骤：
（5.1）生成物体矩阵，针对有效波长
Figure 102953DEST_PATH_IMAGE016
，对物体矩阵的每一点都生成对应的单一波 长
Figure 550115DEST_PATH_IMAGE016
下探测面的光强分布，并将所得各个光强分布叠加，得到物体在该波长下对应的光强 分布
Figure 411892DEST_PATH_IMAGE017
，从而得到有效波长
Figure 842873DEST_PATH_IMAGE002
，
Figure 674563DEST_PATH_IMAGE013
，
Figure 117177DEST_PATH_IMAGE014
，...，
Figure 274489DEST_PATH_IMAGE003
各自对应的物体光场非相干传播到探测面所 产生的光强分布
Figure 68132DEST_PATH_IMAGE018
，
Figure 703513DEST_PATH_IMAGE019
，
Figure 266212DEST_PATH_IMAGE020
，...，
Figure 594426DEST_PATH_IMAGE021
，其中k=1，2，…，n；
（5.2）将各个有效波长下物体光场非相干传播到探测面所产生的光强分布
Figure 344207DEST_PATH_IMAGE018
，
Figure 48858DEST_PATH_IMAGE019
，
Figure 466064DEST_PATH_IMAGE020
，...，
Figure 965178DEST_PATH_IMAGE021
按照预定比例加权平均，得到物体在宽谱照明下在探测面产生的光强分布。
实施例1：
本实施例中，采用Matlab软件构建宽谱非相干光散射成像探测模拟系统，配合64位普通台式机执行模拟操作。该宽谱非相干光散射成像探测模拟系统对应的实际光路传输及探测装置的结构示意图如图4所示，该系统主要包括光源1，散射介质2，光阑3，面阵探测器4，计算机5。
物面a上的物体被光源1照射，从物面a出射的光经过所述的散射介质2和光阑3后入射到所述的面阵探测器4，所述计算机5与所述面阵探测器4相连，用于存储所探测的光强信号即上述步骤（5.2）中得到的物体在宽谱照明下在探测面产生的光强分布并进行有关计算。
步骤（1）：设置宽谱非相干光散射成像探测模拟系统的系统参数，其中，选择数字“5”作为物体形状，设置物面阵列尺寸为300*300、物面最小单元尺寸为1μm，设置物体到散射介质的距离u为6.6cm，设置散射介质面阵列尺寸为600*600、最小单元尺寸为2μm，设置散射介质特征函数，设置光阑直径为1.2mm，设置散射介质到探测面的距离v为4cm，设置探测面阵列大小为2000*2000、探测面像素大小为3.45μm，设置波长范围参数为501nm-600nm；
步骤（2）：选择物面中心点作为标定目标点；
步骤（3）：对波长范围为
Figure 467835DEST_PATH_IMAGE001
，即501nm-600nm，的宽谱光源，按照波长间距
Figure 179439DEST_PATH_IMAGE004
= 1nm对波长进行均匀的离散性操作，得到波长为
Figure 841364DEST_PATH_IMAGE002
+
Figure 383817DEST_PATH_IMAGE025
（504nm），...，
Figure 498403DEST_PATH_IMAGE003
（600nm）的一系列单色波长；通过以下过程生成各个波长对应的点扩展函 数：
将标定目标点按照菲涅尔衍射的方式传播距离u（6.6cm）得到散射介质前表面的 入射光场
Figure 154644DEST_PATH_IMAGE006
；将入射光场
Figure 405496DEST_PATH_IMAGE006
先后和散射介质特征函数以及光阑特征函数相乘，得到散射介 质后表面的出射光场
Figure 121780DEST_PATH_IMAGE007
；对光场
Figure 723662DEST_PATH_IMAGE007
按照菲涅尔衍射的方式光场传输距离v（4cm），得到探 测面的光场分布
Figure 183594DEST_PATH_IMAGE008
，将光场
Figure 554532DEST_PATH_IMAGE008
与其共轭光场
Figure 441717DEST_PATH_IMAGE009
相乘，得到单一波长下位于所述标定目标点 的点光源在探测面处所形成的光强度分布，即该波长下的点扩展函数。
步骤（4）：将
Figure 265316DEST_PATH_IMAGE002
作为第一个有效波长，将波长
Figure 528938DEST_PATH_IMAGE002
+
Figure 754383DEST_PATH_IMAGE004
，
Figure 812469DEST_PATH_IMAGE002
+
Figure 388944DEST_PATH_IMAGE010
...
Figure 49732DEST_PATH_IMAGE011
对应的点扩展 函数依次与
Figure 5050DEST_PATH_IMAGE002
对应的点扩展函数做相关性计算，选择相关性初次低于相关性阈值的波长作 为第二个有效波长
Figure 624250DEST_PATH_IMAGE013
（506nm）；将波长
Figure 294879DEST_PATH_IMAGE013
+
Figure 759358DEST_PATH_IMAGE004
（507nm），
Figure 569182DEST_PATH_IMAGE013
+
Figure 359284DEST_PATH_IMAGE026
...，
Figure 254558DEST_PATH_IMAGE003
对应的点 扩展函数依次与
Figure 788308DEST_PATH_IMAGE013
对应的点扩展函数做相关性计算，选择相关性初次低于相关性阈值的波 长作为第三个有效波长
Figure 718218DEST_PATH_IMAGE014
（513nm）；以此类推，得到有效波长
Figure 413641DEST_PATH_IMAGE002
（501nm），
Figure 920846DEST_PATH_IMAGE013
（506nm），
Figure 868073DEST_PATH_IMAGE014
（513nm），...，
Figure 777124DEST_PATH_IMAGE015
（592nm）；最大波长
Figure 518815DEST_PATH_IMAGE003
（600nm）即使不满足所述相关性低于相关性阈值的 条件也始终保留作为有效波长；其中，所述相关性阈值取0.9。
步骤（5）：生成物体矩阵，针对有效波长
Figure 513315DEST_PATH_IMAGE016
，对物体矩阵的每一点都生成对应的单 一波长
Figure 264234DEST_PATH_IMAGE016
下探测面的光强分布，并将所得各个光强分布叠加，得到物体在该波长下对应的 光强分布
Figure 293370DEST_PATH_IMAGE017
，从而得到有效波长
Figure 205962DEST_PATH_IMAGE002
，
Figure 687759DEST_PATH_IMAGE013
，
Figure 101423DEST_PATH_IMAGE014
，...，
Figure 857502DEST_PATH_IMAGE003
各自对应的物体光场非相干传播到探测 面所产生的光强分布
Figure 65629DEST_PATH_IMAGE018
，
Figure 644509DEST_PATH_IMAGE019
，
Figure 861864DEST_PATH_IMAGE020
，...，
Figure 475379DEST_PATH_IMAGE021
，其中k=1，2，…，n；将各个有效波长下物体光场非相 干传播到探测面所产生的光强分布
Figure 119987DEST_PATH_IMAGE018
，
Figure 451742DEST_PATH_IMAGE019
，
Figure 941629DEST_PATH_IMAGE020
，...，
Figure 799864DEST_PATH_IMAGE021
按照预定比例加权平均，得到物体在宽 谱照明下在探测面产生的光强分布，如图3所示。本实施例不针对某一个固定光源，所以按 照1:1:1…:1比例平均加权。
从图3可以看出，探测图案从中心向周围呈放射状发散特性，满足宽谱光成像的主要特征，说明了本方法的有效性，即利用波长之间的相关性作为筛选条件，只对相关性较低的各个波长分步执行成像仿真操作，可以在较少重复步骤下实现有效的宽谱光学模拟，极大地降低模拟中的冗余，减小计算量，提高模拟效率；可以应用在散斑自相关等一系列使用非相干光的散射成像模拟实验中，将其与相位恢复算法等图像重建方式相结合，能够为实际应用提供参考，在水下探测，透雾成像等方面都有很大应用前景。
本领域技术人员在考虑说明书及实践这里公开的内容后，将容易想到本申请的其它实施方案。本申请旨在涵盖本申请的任何变型、用途或者适应性变化，这些变型、用途或者适应性变化遵循本申请的一般性原理并包括本申请未公开的本技术领域中的公知常识或惯用技术手段。
应当理解的是，本申请并不局限于上面已经描述并在附图中示出的精确结构，并且可以在不脱离其范围的基础上进行各种修改和改变。

Claims (6)

1.一种宽谱非相干光散斑自相关成像探测的低冗余模拟方法，其特征在于，应用于宽谱非相干光散射成像探测模拟系统，包括：
步骤（1）：设置宽谱非相干光散射成像探测模拟系统的系统参数；
步骤（2）：选择物面中心点作为标定目标点；
步骤（3）：将波长范围为
Figure 865434DEST_PATH_IMAGE001
的宽谱光源按照预定波长间距对波长进行均匀的离散 性操作，并生成所述标定目标点对应的各个波长的点扩展函数；
步骤（4）：从
Figure 543540DEST_PATH_IMAGE002
开始，选出点扩展函数间相关性低于相关性阈值的波长作为采样模拟中 的有效波长，并始终保留最大波长
Figure 978064DEST_PATH_IMAGE003
作为有效波长；
步骤（5）：生成各个有效波长对应的物体光场非相干传播到所述探测面所产生的光强分布并按照各自的预定权重进行加权平均，得到物体在宽谱照明下在探测面产生的光强分布。
2.根据权利要求1所述的方法，其特征在于，步骤（1）中的所述系统参数包括物体形状，物面阵列尺寸、物面最小单元尺寸，物体到散射介质的距离u，散射介质面阵列尺寸、最小单元尺寸，散射介质特征函数，光阑直径，散射介质到探测面的距离v，探测面阵列大小、探测面像素大小，波长范围参数。
3.根据权利要求1所述的方法，其特征在于，步骤（3）包括：
对波长范围为
Figure 690805DEST_PATH_IMAGE001
的宽谱光源，按照波长间距
Figure 552582DEST_PATH_IMAGE004
对波长进行均匀的离散性操作，得 到波长为
Figure 983563DEST_PATH_IMAGE002
+
Figure 953269DEST_PATH_IMAGE005
，...，
Figure 520516DEST_PATH_IMAGE003
的一系列单色波长，其中，
Figure 553194DEST_PATH_IMAGE004
在0.2nm-2nm之 间；
生成波长为
Figure 940313DEST_PATH_IMAGE002
+
Figure 982219DEST_PATH_IMAGE005
，...，
Figure 403973DEST_PATH_IMAGE003
的一系列单色波长中各个波长对应的 点扩展函数。
4.根据权利要求1所述的方法，其特征在于，生成点扩展函数的过程包括：
将标定目标点按照菲涅尔衍射的方式传播距离u得到散射介质前表面的入射光场
Figure 873131DEST_PATH_IMAGE006
；
将入射光场
Figure 154071DEST_PATH_IMAGE006
先后和散射介质特征函数以及光阑特征函数相乘，得到散射介质后表面 的出射光场
Figure 593143DEST_PATH_IMAGE007
；
对光场
Figure 10348DEST_PATH_IMAGE007
按照菲涅尔衍射的方式光场传输距离v，得到探测面的光场分布
Figure 509463DEST_PATH_IMAGE008
，将光场
Figure 12120DEST_PATH_IMAGE008
与其共轭光场
Figure 989303DEST_PATH_IMAGE009
相乘，得到单一波长下位于所述标定目标点的点光源在探测面处所形 成的光强度分布，即该波长下的点扩展函数。
5.根据权利要求1所述的方法，其特征在于，步骤（4）包括：
将
Figure 651228DEST_PATH_IMAGE002
作为第一个有效波长，将波长
Figure 196610DEST_PATH_IMAGE002
+
Figure 45618DEST_PATH_IMAGE004
，
Figure 987945DEST_PATH_IMAGE002
+
Figure 238798DEST_PATH_IMAGE010
...
Figure 955081DEST_PATH_IMAGE011
对应的点扩展函数依次与
Figure 291385DEST_PATH_IMAGE002
对应的点扩展函数做相关性计算，选择相关性初次低于相关性阈值的波长作为第二个有效 波长
Figure 751316DEST_PATH_IMAGE012
；将波长
Figure 122254DEST_PATH_IMAGE012
+
Figure 9439DEST_PATH_IMAGE004
，
Figure 833038DEST_PATH_IMAGE012
+
Figure 955715DEST_PATH_IMAGE010
...，
Figure 322106DEST_PATH_IMAGE003
对应的点扩展函数依次与
Figure 239246DEST_PATH_IMAGE012
对应的点扩展函数做相 关性计算，选择相关性初次低于相关性阈值的波长作为第三个有效波长
Figure 425508DEST_PATH_IMAGE013
；以此类推，得到 有效波长
Figure 351876DEST_PATH_IMAGE002
，
Figure 572772DEST_PATH_IMAGE012
，
Figure 926393DEST_PATH_IMAGE013
，...，
Figure 459006DEST_PATH_IMAGE014
；最大波长
Figure 64431DEST_PATH_IMAGE003
即使不满足所述相关性低于相关性阈值的条件 也始终保留作为有效波长；
其中，所述相关性阈值不大于1。
6.根据权利要求1所述的方法，其特征在于，步骤（5）包括：
生成物体矩阵，针对有效波长
Figure 264468DEST_PATH_IMAGE015
，对物体矩阵的每一点都生成对应的单一波长
Figure 661427DEST_PATH_IMAGE015
下探 测面的光强分布，并将所得各个光强分布叠加，得到物体在该波长下对应的光强分布
Figure 946914DEST_PATH_IMAGE016
，从 而得到有效波长
Figure 90451DEST_PATH_IMAGE002
，
Figure 144995DEST_PATH_IMAGE012
，
Figure 715784DEST_PATH_IMAGE013
，...，
Figure 488568DEST_PATH_IMAGE003
各自对应的物体光场非相干传播到探测面所产生的光强 分布
Figure 435796DEST_PATH_IMAGE017
，
Figure 344846DEST_PATH_IMAGE018
，
Figure 211171DEST_PATH_IMAGE019
，...，
Figure 81038DEST_PATH_IMAGE020
，其中k=1，2，…，n；
将各个有效波长下物体光场非相干传播到探测面所产生的光强分布
Figure 956590DEST_PATH_IMAGE017
，
Figure 595513DEST_PATH_IMAGE018
，
Figure 898318DEST_PATH_IMAGE019
，...，
Figure 989902DEST_PATH_IMAGE020
按 照预定比例加权平均，得到物体在宽谱照明下在探测面产生的光强分布。
CN202211226822.2A 2022-10-09 2022-10-09 一种宽谱非相干光散斑自相关成像探测的低冗余模拟方法 Active CN115290601B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202211226822.2A CN115290601B (zh) 2022-10-09 2022-10-09 一种宽谱非相干光散斑自相关成像探测的低冗余模拟方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202211226822.2A CN115290601B (zh) 2022-10-09 2022-10-09 一种宽谱非相干光散斑自相关成像探测的低冗余模拟方法

Publications (2)

Publication Number Publication Date
CN115290601A true CN115290601A (zh) 2022-11-04
CN115290601B CN115290601B (zh) 2023-01-24

Family

ID=83833492

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202211226822.2A Active CN115290601B (zh) 2022-10-09 2022-10-09 一种宽谱非相干光散斑自相关成像探测的低冗余模拟方法

Country Status (1)

Country Link
CN (1) CN115290601B (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116242252A (zh) * 2023-05-11 2023-06-09 之江实验室 一种具有定位和尺寸测量功能的散射成像方法

Citations (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20130100525A1 (en) * 2011-10-19 2013-04-25 Su Yu CHIANG Optical imaging system using structured illumination
CN107589542A (zh) * 2017-09-21 2018-01-16 中国科学院长春光学精密机械与物理研究所 宽波段相位差图像重建中的中心波长的选择方法
CN107907483A (zh) * 2017-08-14 2018-04-13 西安电子科技大学 一种基于散射介质的超分辨光谱成像系统及方法
CN110285755A (zh) * 2019-07-10 2019-09-27 吉林大学 一种用于筛选关联成像采样数据的装置和方法
US20200348229A1 (en) * 2019-04-19 2020-11-05 Tsinghua Shenzhen International Graduate School Non-invasive scattering imaging method beyond memory effect range based on connected component optimization
CN112051239A (zh) * 2019-06-05 2020-12-08 中国科学院上海光学精密机械研究所 基于动态散射系统的在探测面积受限情况下的成像方法
CN112161953A (zh) * 2020-08-25 2021-01-01 西安电子科技大学 一种基于散射介质的宽光谱单帧散射成像方法
WO2021144795A1 (en) * 2020-01-14 2021-07-22 Pxe Computation Imaging Ltd. System and method for optical imaging and measurement of objects
CN113218914A (zh) * 2021-03-24 2021-08-06 杭州电子科技大学 一种非侵入式散射介质点扩展函数获取装置及方法
CN113376993A (zh) * 2021-06-18 2021-09-10 哈尔滨理工大学 编码孔径相关全息术的虚拟点扩散函数记录方法

Patent Citations (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20130100525A1 (en) * 2011-10-19 2013-04-25 Su Yu CHIANG Optical imaging system using structured illumination
CN107907483A (zh) * 2017-08-14 2018-04-13 西安电子科技大学 一种基于散射介质的超分辨光谱成像系统及方法
CN107589542A (zh) * 2017-09-21 2018-01-16 中国科学院长春光学精密机械与物理研究所 宽波段相位差图像重建中的中心波长的选择方法
US20200348229A1 (en) * 2019-04-19 2020-11-05 Tsinghua Shenzhen International Graduate School Non-invasive scattering imaging method beyond memory effect range based on connected component optimization
CN112051239A (zh) * 2019-06-05 2020-12-08 中国科学院上海光学精密机械研究所 基于动态散射系统的在探测面积受限情况下的成像方法
CN110285755A (zh) * 2019-07-10 2019-09-27 吉林大学 一种用于筛选关联成像采样数据的装置和方法
WO2021144795A1 (en) * 2020-01-14 2021-07-22 Pxe Computation Imaging Ltd. System and method for optical imaging and measurement of objects
CN112161953A (zh) * 2020-08-25 2021-01-01 西安电子科技大学 一种基于散射介质的宽光谱单帧散射成像方法
CN113218914A (zh) * 2021-03-24 2021-08-06 杭州电子科技大学 一种非侵入式散射介质点扩展函数获取装置及方法
CN113376993A (zh) * 2021-06-18 2021-09-10 哈尔滨理工大学 编码孔径相关全息术的虚拟点扩散函数记录方法

Non-Patent Citations (6)

* Cited by examiner, † Cited by third party
Title
AARON M. KHO ET AL: "Incoherent excess noise spectrally encodes broadband light sources", 《LIGHT: SCIENCE & APPLICATIONS》 *
RUI MA ET AL: "Wavelength-dependent speckle multiplexing for imaging through opacity", 《OPTICS AND LASERS IN ENGINEERING》 *
SUJIT KUMAR SAHOO ET AL.: "Single-shot multispectral imaging with a monochromatic camera", 《OPTICA》 *
XIAOQING XU ET AL.: "Imaging of objects through a thin scattering layer using a spectrally and spatially separated reference", 《OPTICS EXPRESS》 *
王歆等: "基于多波长照明的透过散射介质成像", 《光学学报》 *
谢向生等: "散斑相关成像:从点扩展函数到光场全要素", 《光学学报》 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116242252A (zh) * 2023-05-11 2023-06-09 之江实验室 一种具有定位和尺寸测量功能的散射成像方法
CN116242252B (zh) * 2023-05-11 2023-08-15 之江实验室 一种具有定位和尺寸测量功能的散射成像方法

Also Published As

Publication number Publication date
CN115290601B (zh) 2023-01-24

Similar Documents

Publication Publication Date Title
Huang et al. Continuous-wave terahertz multi-plane in-line digital holography
CN108169173A (zh) 一种大视场高分辨三维衍射层析显微成像方法
CN115290601B (zh) 一种宽谱非相干光散斑自相关成像探测的低冗余模拟方法
CN106990694A (zh) 一种部分相干光照明下的非迭代相位恢复装置及方法
CN110243398A (zh) 一种基于相干检测的相位物体的关联成像方法
CN109521580A (zh) 基于涡旋光子筛望远镜系统的多层成像装置和成像方法
CN112051239B (zh) 基于动态散射系统的在探测面积受限情况下的成像方法
Zhang et al. Application of short-coherence lensless Fourier-transform digital holography in imaging through diffusive medium
CN111141706A (zh) 一种透过散射介质宽场的三维成像方法
CN112229822B (zh) 对流动液体内多目标的反射式单帧散射成像装置及方法
Gröger et al. Advantages of holographic imaging through fog
Xin et al. Research on triggering properties enhancement of polarization detection geiger-mode APD LIDAR
Hajjarian et al. Spatially multiplexed multi-input-multi-output optical imaging system in a turbid, turbulent atmosphere
CN110488397A (zh) 复合透镜、可编码单帧成像装置及成像方法
Li Target infrared radiation calculation model and method based on finite element analysis method in infrared photoelectric detection system
Yuan et al. Multiple-object ghost imaging with a single-pixel detector
Vollmerhausen et al. Modeling the target acquisition performance of laser-range-gated imagers
CN114966739A (zh) 基于光纤激光相控阵的目标探测成像系统及方法
CN111596465A (zh) 一种漂移光束计量超分辨鬼成像质量的装置
CN114518658A (zh) 基于全光场调控的散射介质光场聚焦方法与装置
CN114518659A (zh) 基于复振幅光场调控的散射介质光场聚焦方法与装置
CN113252613A (zh) 一种基于波前整形的散射介质吸收增强装置及方法
CN114764220B (zh) 一种基于离轴数字全息改善散斑自相关重建效果的方法
Meijiao et al. Study on second order correlation imaging technique based on He-Ne laser
CN116337225B (zh) 基于涡旋光的提高光电信号探测信噪比的方法及实验装置

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant