CN114813643A - 一种基于预置散射介质的穿散射介质关联成像方法及系统 - Google Patents

一种基于预置散射介质的穿散射介质关联成像方法及系统 Download PDF

Info

Publication number
CN114813643A
CN114813643A CN202210258421.9A CN202210258421A CN114813643A CN 114813643 A CN114813643 A CN 114813643A CN 202210258421 A CN202210258421 A CN 202210258421A CN 114813643 A CN114813643 A CN 114813643A
Authority
CN
China
Prior art keywords
scattering medium
light field
detected
preset
obtaining
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.)
Pending
Application number
CN202210258421.9A
Other languages
English (en)
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.)
Beijing Institute of Technology BIT
Original Assignee
Beijing Institute of Technology BIT
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 Beijing Institute of Technology BIT filed Critical Beijing Institute of Technology BIT
Priority to CN202210258421.9A priority Critical patent/CN114813643A/zh
Publication of CN114813643A publication Critical patent/CN114813643A/zh
Pending legal-status Critical Current

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N21/00Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
    • G01N21/17Systems in which incident light is modified in accordance with the properties of the material investigated
    • G01N21/47Scattering, i.e. diffuse reflection
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N21/00Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
    • G01N21/17Systems in which incident light is modified in accordance with the properties of the material investigated
    • G01N21/59Transmissivity
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04NPICTORIAL COMMUNICATION, e.g. TELEVISION
    • H04N25/00Circuitry of solid-state image sensors [SSIS]; Control thereof
    • H04N25/70SSIS architectures; Circuits associated therewith
    • H04N25/71Charge-coupled device [CCD] sensors; Charge-transfer registers specially adapted for CCD sensors

Landscapes

  • Chemical & Material Sciences (AREA)
  • Biochemistry (AREA)
  • Pathology (AREA)
  • Physics & Mathematics (AREA)
  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Immunology (AREA)
  • Analytical Chemistry (AREA)
  • General Physics & Mathematics (AREA)
  • General Health & Medical Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Multimedia (AREA)
  • Signal Processing (AREA)
  • Length Measuring Devices By Optical Means (AREA)

Abstract

本发明公开一种基于预置散射介质的穿散射介质关联成像方法及系统,方法包括以下步骤:获取激光光束,基于所述激光光束,获得赝热光场;基于所述赝热光场,获得探测光场和参考光场;基于环境中的散射介质和所述探测光场,获得待测物体光强;基于预置散射介质和所述参考光场,获得散斑光强;基于所述待测物体光强和所述散斑光强,获得重构的待测物体图像。本发明为了实现穿散射介质关联成像,提出了一种基于预置散射介质的穿散射介质关联成像方法,具有成本低,系统复杂度低的特点。

Description

一种基于预置散射介质的穿散射介质关联成像方法及系统
技术领域
本发明属于计算成像技术领域,尤其涉及一种基于预置散射介质的穿散射介质关联成像方法及系统
背景技术
生活中常见的雾、霾、沙尘等散射介质会改变光的传播路径,导致基于“点对点成像机制”的传统成像方式成像质量下降。穿透散射介质成像,在遥感、医学成像、浓烟环境中的成像等场景中扮演重要的角色,具有巨大的应用潜力。
现有的穿散射介质成像技术中,反馈优化波前法时效性不高,系统复杂;基于传输矩阵的重聚焦技术系统复杂且稳定性要求较高;基于光学记忆效应的散斑自相关成像技术的成像视场较小,受到散射介质记忆效应范围的限制。
关联成像,是一种非局域性的成像技术。赝热照明光场通过分束器分成探测光场和参考光场,探测光通过探测光路与物体相互作用后的总光强值被桶探测器探测,参考光场通过参考光路,其光场分布被CCD记录,进行多次测量,最后对测得的信号进行二阶关联运算即可重构出目标的图像。这种成像机制具有“离物成像”的优势,为实现穿散射介质成像提供了新思路。
相关研究人员在理论、实验方面进行了许多探索,研究结果表明,与传统之间成像相比,关联成像对位于目标物体与桶探测器之间的散射介质具有抵抗能力。然而,位于分束器与物体之间的散射介质会对使探测光场的光场分布与参考光路中CCD相机采集的参考光场分布相关性降低甚至失去相关性,从而对关联成像的成像结果造成干扰。针对这个问题,国内外也有很多相关报道,但是都存在一些劣势,比如有的仍需知道散射介质的传输矩阵、有的需要知道待成像目标的形状等先验信息等。目前报道的方法都未较好的利用关联成像的非局域特性,本文提出“将参考臂中关于分束器对称的位置引入与探测臂中存在的相同特性的散射介质”的方案,实现了预置散射介质对参考光场进行调制后进行二阶关联计算,进而实现穿散射介质关联成像。
发明内容
本发明公开一种基于预置散射介质的穿散射介质关联成像方法及系统。基于关联成像非局域性的成像机制,在参考光路中,关于分束器对称的位置预置一种与存在于探测光路中的散射介质散射特性相同的散射介质,再进行二阶关联重构,实现穿散射介质关联成像。
为实现上述目的,一方面,本发明提供了一种基于预置散射介质的穿散射介质关联成像方法,包括以下步骤:
获取激光光束,基于所述激光光束,获得赝热光场;
基于所述赝热光场,获得探测光场和参考光场;
基于环境中的散射介质和所述探测光场,模拟探测光场通过环境中散射介质的场景,并获得透过待测物体的总光强;
基于预置散射介质和所述参考光场,获得参考散斑光场;
基于所述透过待测物体的总光强和所述参考散斑光场,通过二阶关联算法重构出待测物体图像。
可选地,获得透过待测物体的总光强的过程为:所述探测光场通过所述环境中的散射介质后,再照明待测物体,获得所述透过待测物体的总光强。
可选地,所述透过待测物体总光强的分布计算公式为:
Figure BDA0003549297530000031
可选地,获得参考散斑光场的过程为:所述参考光场通过所述预置散射介质后,获得所述参考散斑光场。
可选地,所述参考散斑光场的计算公式为:
Figure BDA0003549297530000032
可选地,重构待测物体图像的过程为:对所述透过待测物体的总光强和所述参考散斑光场进行二阶关联计算,重构出待测物体图像。
可选地,所述二阶关联计算公式为:
H=<a′-<A′》<b′-<B′》。
另一方面,本发明提供了一种基于预置散射介质的穿散射介质关联成像系统,包括:激光器、旋转毛玻璃、分束器、环境中的散射介质、预置散射介质、桶探测器、CCD相机和成像模块;
所述激光器用于产生激光光束;
所述旋转毛玻璃用于基于所述激光光束,获得赝热光场;
所述分束器用于基于所述赝热光场,获得探测光场和参考光场;
所述环境中的散射介质用于模拟所述探测光场通过散射介质时的场景,通过散射介质后的探测光场照明待测物体后,获得透过待测物体总光强;
所述预置散射介质用于基于所述预置散射介质和所述参考光场,获得参考散斑光场;
所述桶探测器用于探测所述透过待测物体的总光强;
所述CCD相机用于采集所述参考散斑光场;
所述成像模块用于基于所述透过待测物体的总光强和所述参考散斑光场,重构出待测物体图像。
与现有技术相比,本发明技术效果:
本发明为了实现穿散射介质关联成像,提出了一种基于预置散射介质的穿散射介质关联成像方法,具有成本低,系统复杂度低的特点。
附图说明
构成本申请的一部分的附图用来提供对本申请的进一步理解,本申请的示意性实施例及其说明用于解释本申请,并不构成对本申请的不当限定。在附图中:
图1是本发明所述的关联成像的装置示意图。
图2是本发明所述的当分束器与待测物体之间存在散射介质时的装置示意图。
图3是本发明所述的在参考光路中预置与探测光路中的散射介质具有相同散射特性的散射介质时的装置示意图。
图4是成像结果对比图。其中:(a)是光路中无散射介质时的成像结果,可以重构出待测物体图像;(b)是分束器与物体之间存在散射介质时的成像结果,无法重构出待测物体图像;(c)是本发明所述的方法下的成像结果,可以重构出待测物体图像。
具体实施方式
需要说明的是,在不冲突的情况下,本申请中的实施例及实施例中的特征可以相互组合。下面将参考附图并结合实施例来详细说明本申请。
实施例1
现有技术中,如图1所示,当光路中不存在散射介质时,激光透过旋转毛玻璃后生成赝热光场,然后通过分束器分为探测光场和参考光场。探测光场在探测光路中照明待测物体后,透射的总光强值被桶探测器探测,参考光场在参考光路中被CCD相机采集。设待测物体H的尺寸为m×n,每一个像素位置记为O(x,y),以矩阵形式可作一列进行表示:
H=[O(1,1) O(2,1) ... O(m,n)]T
参考光路中CCD相机采集的散斑图像尺寸也为m×n,每一个像素位置的光强值记为IR(x,y),第k次探测到的参考光场分布以矩阵形式可作一行表示:
Figure BDA0003549297530000051
探测光场照明待测物体后,透过的总光强值被桶探测器探测,第k次测量结果满足如下数学关系式:
b(k)=a(k)H
进行K次测量后,采集到的所有参考光场的光场分布可表示为:
Figure BDA0003549297530000061
采集到的对应桶探测值可表示为:
B=[b(1) b(2) ... b(K)]T
满足数学关系式:
B=AH
通过以下数学运算可解算出待测物体图像矩阵:
H=A-1B
如图2所示,当分束器与物体之间存在散射介质时,散射介质对光场的调制过程可用该散射介质的传输矩阵F进行表示。激光透过旋转毛玻璃后生成赝热光场,然后通过分束器分为探测光场和参考光场。探测光场在探测光路中通过散射介质后,第k次测量时,每个像素位置光强I′T(x,y)=IT(x,y)·F,K次测量后可表示为:
Figure BDA0003549297530000062
桶探测器采集到的透过待测物体的总共强,即桶探测值为:
B′=(A·F)H=GH
而利用CCD相机K次测量后,采集到的所有参考光场的光场分布与该桶探测值进行关联计算后,重构出的结果为:
H′=A-1B′=A-1·GH=FH≠H
即无法重构出待测目标物体。
如图3所示,本发明提供一种基于预置散射介质的穿散射介质关联成像方法,方法流程为:获取激光光束,基于激光光束,获得赝热光场;基于赝热光场,获得探测光场和参考光场;基于环境中的散射介质,模拟所述探测光场通过散射介质时的场景,通过散射介质后的探测光场照明待测物体后,获得透过待测物体总光强;基于预置散射介质,参考光场通过该预置的散射介质,获得参考散斑光场。基于透过待测物体的总光强和参考散斑光场,通过二阶关联算法重构出待测物体图像。
具体的,在参考光路中预置与探测光路中散射介质散射特性相同的散射介质,该散射介质对参考光场的调制也可用散射介质传输矩阵F表示,参考光场通过该预置散射介质后,第k次测量时每个像素位置光强I′R(x,y)=IR(x,y)·F,K次测量后可表示为A′=AF,即:
Figure BDA0003549297530000071
利用CCD相机K次测量后,采集到的所有参考散斑光场与该桶探测值进行关联计算后,重构出的结果为:
H″=A′-1B′=A′-1·(A·F)H=(AF)-1·(A·F)H=H
即通过对测得的两部分信号进行二阶关联计算:H=<a′-<A′>><b′-<B′>>可以重构出待测物体图像,实现穿透散射介质关联成像。
实施例2
现有技术中,如图1所示,关联成像的装置示意图,主要的组成部分包括:532nm激光器、旋转毛玻璃、50:50分束器、CCD相机、待测物体镂空双缝以及桶探测器。
激光器发出激光打在旋转毛玻璃上产生赝热光场,通过分束器分为两部分,一部分为探测光场,另一部分为参考光场。待测物体与分束器距离为Z1,CCD相机与分束器距离也为Z1,桶探测器探测面紧贴待测物体。探测光场在探测光路中传播距离Z1后与待测物体物体相互作用,作用后光场的总光强值被桶探测器探测。参考光场在参考光路中传播距离Z1后被CCD相机采集。经过780次采样,将采集到的信号进行二阶关联计算,重构结果如图4(a)。
如图2所示,当分束器与待测物体之间存在散射介质时的装置示意图。盛放散射介质的容器(散射介质1)与分束器的距离为Z2,探测光场与散射介质1作用后,再与待测物体物体相互作用,作用后的总光强值被桶探测器探测。参考光场传播距离Z1后被CCD相机采集。经过780次采样,将采集到的信号进行二阶关联计算,重构结果如图4(b)。
如图3所示,本发明提供一种基于预置散射介质的穿散射介质关联成像系统,包括:532nm激光器、旋转毛玻璃、50:50分束器、环境中的散射介质、预置散射介质模块、待测物体镂空双缝、桶探测器和CCD相机;
当分束器与待测物体之间存在散射介质时,本发明所述的在参考光路中预置与探测光路中的散射介质具有相同散射特性的散射介质时的装置示意图。盛放预置散射介质2的容器,即预置散射介质模块与分束器的距离也为Z2,参考光场与预置的散射2相互作用后,传播到CCD像面所在位置时的光场分布被CCD相机采集。经过780次采样,将采集到的信号进行二阶关联计算H=<a′-<A′>><b′-<B′>>,重构结果如图4(c),尺寸为128×128。
进一步的,计算重构结果的衬噪比(CNR)和可见度(V),作为评价指标。衬噪比的计算公式为:
Figure BDA0003549297530000091
可见度的计算公式为:
Figure BDA0003549297530000092
其中,H(xin)和H(xout)分别代表镂空双缝透射区域和反射区域;
如图4所示为成像结果对比图。重构结果的衬噪比和可见度以(CNR,V)的形式标记在每副图像的下方。其中:(a)是光路中无散射介质时的成像结果,可以重构出待测物体图像;(b)是分束器与物体之间存在散射介质时的成像结果,无法重构出待测物体图像,衬噪比和可见度明显减小;(c)是本发明所述的方法下的成像结果,可以重构出待测物体图像,衬噪比和可见度与(b)相比有所提高。
本发明为了实现穿散射介质关联成像,提出了一种基于预置散射介质的穿散射介质关联成像方法,具有成本低,系统复杂度低的特点。
以上所述,仅为本申请较佳的具体实施方式,但本申请的保护范围并不局限于此,任何熟悉本技术领域的技术人员在本申请揭露的技术范围内,可轻易想到的变化或替换,都应涵盖在本申请的保护范围之内。因此,本申请的保护范围应该以权利要求的保护范围为准。

Claims (8)

1.一种基于预置散射介质的穿散射介质关联成像方法,其特征在于,包括以下步骤:
获取激光光束,基于所述激光光束,获得赝热光场;
基于所述赝热光场,获得探测光场和参考光场;
基于环境中的散射介质和所述探测光场,模拟探测光场通过环境中散射介质的场景,并获得透过待测物体的总光强;
基于预置散射介质和所述参考光场,获得参考散斑光场;
基于所述透过待测物体的总光强和所述参考散斑光场,通过二阶关联算法重构出待测物体图像。
2.根据权利要求1所述的基于预置散射介质的穿散射介质关联成像方法,其特征在于,获得透过待测物体总光强的过程为:所述探测光场通过所述探测光路中存在的散射介质后,再照明待测物体获得透过所述待测物体的总光强。
3.根据权利要求2所述的基于预置散射介质的穿散射介质关联成像方法,其特征在于,所述透过待测物体总光强的分布计算公式为:
Figure FDA0003549297520000011
4.根据权利要求3所述的基于预置散射介质的穿散射介质关联成像方法,其特征在于,获得参考散斑光场的过程为:所述参考光场通过所述预置散射介质后,获得所述参考散斑光场。
5.根据权利要求4所述的基于预置散射介质的穿散射介质关联成像方法,其特征在于,所述参考散斑光场的计算公式为:
Figure FDA0003549297520000021
6.根据权利要求5所述的基于预置散射介质的穿散射介质关联成像方法,其特征在于,重构待测物体图像的过程为:对所述透过待测物体总光强和所述参考散斑光场进行二阶关联计算,重构出待测物体图像。
7.根据权利要求6所述的基于预置散射介质的穿散射介质关联成像方法,其特征在于,所述二阶关联计算公式为:
H=<a′-<A′>><b′-<B′>>。
8.一种基于预置散射介质的穿散射介质关联成像系统,其特征在于,包括:激光器、旋转毛玻璃、分束器、环境中的散射介质、预置散射介质、桶探测器、CCD相机和成像模块;
所述激光器用于产生激光光束;
所述旋转毛玻璃用于基于所述激光光束,获得赝热光场;
所述分束器用于基于所述赝热光场,获得探测光场和参考光场;
所述环境中的散射介质用于模拟所述探测光场通过散射介质时的场景,通过散射介质后的探测光场照明待测物体后,获得透过待测物体总光强;
所述预置散射介质用于基于所述预置散射介质和所述参考光场,获得参考散斑光场;
所述桶探测器用于探测所述透过待测物体的总光强;
所述CCD相机用于采集所述参考散斑光场;
所述成像模块用于基于所述透过待测物体的总光强和所述参考散斑光场,通过二阶关联算法重构出待测物体图像。
CN202210258421.9A 2022-03-16 2022-03-16 一种基于预置散射介质的穿散射介质关联成像方法及系统 Pending CN114813643A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202210258421.9A CN114813643A (zh) 2022-03-16 2022-03-16 一种基于预置散射介质的穿散射介质关联成像方法及系统

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202210258421.9A CN114813643A (zh) 2022-03-16 2022-03-16 一种基于预置散射介质的穿散射介质关联成像方法及系统

Publications (1)

Publication Number Publication Date
CN114813643A true CN114813643A (zh) 2022-07-29

Family

ID=82528920

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202210258421.9A Pending CN114813643A (zh) 2022-03-16 2022-03-16 一种基于预置散射介质的穿散射介质关联成像方法及系统

Country Status (1)

Country Link
CN (1) CN114813643A (zh)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115553723A (zh) * 2022-09-20 2023-01-03 湖南大学 用于血液中异常细胞筛查的基于高速调制随机介质掺杂光纤的关联成像方法
CN117347324A (zh) * 2023-10-27 2024-01-05 北京菲舍波特科技发展有限公司 一种蒸汽干度在线检测装置

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101699312A (zh) * 2009-09-08 2010-04-28 中国科学院上海光学精密机械研究所 强度关联提高散射介质中物体成像质量的装置
CN103364794A (zh) * 2013-07-12 2013-10-23 华东师范大学 基于关联光学的单次成像及存取方法
CN203365677U (zh) * 2013-07-12 2013-12-25 华东师范大学 基于关联光学的单次成像及存取装置
CN110333205A (zh) * 2019-06-26 2019-10-15 西安电子科技大学 一种透过动态液体介质的散射成像装置和方法

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101699312A (zh) * 2009-09-08 2010-04-28 中国科学院上海光学精密机械研究所 强度关联提高散射介质中物体成像质量的装置
CN103364794A (zh) * 2013-07-12 2013-10-23 华东师范大学 基于关联光学的单次成像及存取方法
CN203365677U (zh) * 2013-07-12 2013-12-25 华东师范大学 基于关联光学的单次成像及存取装置
CN110333205A (zh) * 2019-06-26 2019-10-15 西安电子科技大学 一种透过动态液体介质的散射成像装置和方法

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115553723A (zh) * 2022-09-20 2023-01-03 湖南大学 用于血液中异常细胞筛查的基于高速调制随机介质掺杂光纤的关联成像方法
CN115553723B (zh) * 2022-09-20 2024-05-10 湖南大学 用于血液中异常细胞筛查的基于高速调制随机介质掺杂光纤的关联成像方法
CN117347324A (zh) * 2023-10-27 2024-01-05 北京菲舍波特科技发展有限公司 一种蒸汽干度在线检测装置
CN117347324B (zh) * 2023-10-27 2024-05-31 北京菲舍波特科技发展有限公司 一种蒸汽干度在线检测装置

Similar Documents

Publication Publication Date Title
Kirmani et al. Looking around the corner using transient imaging
CN114813643A (zh) 一种基于预置散射介质的穿散射介质关联成像方法及系统
CN109124615B (zh) 一种可选区高动态激光散斑血流成像装置及方法
US6819436B2 (en) Image capturing apparatus and distance measuring method
CN106772310B (zh) 一种基于计算鬼成像的目标鉴定装置及其方法
CN108780035A (zh) 成像装置以及成像方法
CN105988123B (zh) 一种线扫描成像装置
US20120330162A1 (en) Modulated aperture imaging for automatic moving target detection
CN110455834B (zh) 基于光强传输方程的x射线单次曝光成像装置及方法
CN111596310B (zh) 一种基于点探测的运动目标鬼成像系统与方法
CN210036591U (zh) 一种基于频域oct技术的三维彩色动态成像装置
JP2007510963A (ja) デジタル画像化組立品、及びその方法
CN112585521B (zh) 用于捕获具有湍流衰减的显微全光图像的设备和方法
CN113302541A (zh) 用于捕获任意平面之间的全光图像的过程和装置
CN113790676B (zh) 一种基于编码孔径和光场分布的三维空间光谱成像方法及装置
CN114966739A (zh) 基于光纤激光相控阵的目标探测成像系统及方法
Shi et al. Non-invasive imaging through scattering medium beyond the memory effect via polarization-modulation
CN110703276A (zh) 在强散射条件下的傅里叶成像装置及方法
CN109557070A (zh) 一种基于空间编码光的拉曼成像系统
CN110187498B (zh) 一种真热光关联成像系统
CN111369627B (zh) 一种非侵入式散斑定向成像方法
CN209264563U (zh) 一种折射率显微测量系统
CN102878930A (zh) 一种位相物体位相分布的定量测量方法和装置及其应用
CN109142273A (zh) 一种折射率显微测量系统
CN115494005A (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