CN113189101A - 一种带有负反馈调节的无透镜成像方法 - Google Patents
一种带有负反馈调节的无透镜成像方法 Download PDFInfo
- Publication number
- CN113189101A CN113189101A CN202110459324.1A CN202110459324A CN113189101A CN 113189101 A CN113189101 A CN 113189101A CN 202110459324 A CN202110459324 A CN 202110459324A CN 113189101 A CN113189101 A CN 113189101A
- Authority
- CN
- China
- Prior art keywords
- flat plate
- amplitude
- complex amplitude
- diffraction
- detected
- 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
Links
- 238000003384 imaging method Methods 0.000 title claims abstract description 41
- 230000009125 negative feedback regulation Effects 0.000 title claims abstract description 10
- 238000005457 optimization Methods 0.000 claims abstract description 19
- 238000000034 method Methods 0.000 claims abstract description 13
- 238000001514 detection method Methods 0.000 claims abstract description 4
- 238000004364 calculation method Methods 0.000 claims description 15
- 238000001228 spectrum Methods 0.000 claims description 5
- 230000003287 optical effect Effects 0.000 claims description 4
- 238000011478 gradient descent method Methods 0.000 claims description 3
- 238000012545 processing Methods 0.000 claims description 3
- 238000012360 testing method Methods 0.000 claims description 3
- 238000011084 recovery Methods 0.000 description 5
- 230000000694 effects Effects 0.000 description 4
- 230000007547 defect Effects 0.000 description 3
- 238000010586 diagram Methods 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 238000005259 measurement Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000007796 conventional method Methods 0.000 description 1
- 230000003111 delayed effect Effects 0.000 description 1
- 238000005286 illumination Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012805 post-processing Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
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/84—Systems specially adapted for particular applications
Landscapes
- Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Chemical & Material Sciences (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- General Physics & Mathematics (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Testing Of Optical Devices Or Fibers (AREA)
- Optical Modulation, Optical Deflection, Nonlinear Optics, Optical Demodulation, Optical Logic Elements (AREA)
Abstract
本发明公开了一种带有负反馈调节的无透镜成像方法,该方法在目标域设置一个支持域,并每次迭代过程中支持域之外区域估计的复振幅与迭代输入的复振幅带权重叠加作为下一次迭代的输入;支持域之内待测波前使用非线性优化方法实现待测波前的更新。本发明实现了使用一张强度图实现目标物振幅和相位的精确重建,所需实验系统简单,且无需移动图像传感器,适用于实时波前检测和动态目标物成像。
Description
技术领域
本发明涉及光学测量和成像技术领域,尤其涉及一种带有负反馈调节的无透镜成像方法。
背景技术
计算成像极大地改变了传统成像所见即所得的模式,通过计算机的后期处理,克服了传统成像系统成像不完善的弊病,而且降低了成像系统的尺寸和重量。无透镜成像是计算成像的一种实现方式。无透镜成像方法无需透镜可以直接对待测目标成像,因其高分辨率、大视场及无相差等特性成为计算成像领域的一个研究热点。传统的无透镜成像技术是通过不同衍射距离下的多幅强度图样迭代重建出样品的完整波前信息,即多距离相位恢复技术。多距离相位恢复成像系统存在倾斜照明、收敛迟滞、初始距离无法直接测量、分辨率受限等问题。
鉴于传统多距离无透镜成像方法的弊端,本发明提出了一种带有负反馈调节的无透镜成像方法及系统,使用一幅衍射光斑,无需移动图像采集装置,即可实现对待测目标的直接成像。
发明内容
本发明的目的在于针对现有技术的不足,提供一种带有负反馈调节的无透镜成像方法,使用一幅衍射光斑,无需移动图像采集装置,即可实现对待测目标的直接成像,整套装置的成本低,且可以用于实时动态成像。
本发明的目的是通过以下技术方案来实现的:
一种带有负反馈调节的无透镜成像方法,该方法在激光器的出射光路的光轴上依次布置有扩束器、待测平板、图像传感器,该方法包括如下步骤:
S1:使用所述图像传感器在某一位置上采集一幅含有所述的待测平板信息的衍射光斑;
S2:分别设置所述待测平板的支持域D、所述待测平板与所述图像传感器之间的距离d、负反馈权重系数β、相位优化步长stepphase、振幅优化步长stepamp、待测平板重建的迭代总数N、波前检测的初始迭代计数i=1、待测平板重建复振幅初始估计值g1(x,y),其中,(x,y)为所述待测平板的坐标;
S3:构建优化目标函数E
其中,I是在d位置处采集的衍射光斑,U是待测平板在位置d处的衍射积分;W(u,v)是权重分布,用来移除图像传感器上信噪比低的像素点和坏点;
S4:对所述待测平板的支持域D内的部分,首先采用公式(2)~(5)计算波前优化的相位梯度和振幅梯度Δa,然后更新待测波前的振幅和相位,并通过公式(6)对估计的波前进行梯度优化,得到所述支持域D内的复振幅
其中,gi(x,y)表示第i次迭代的待测平板的复振幅输入估计值,Gi(x,y)代表从待测面衍射计算到图像采集面的复振幅,代表从图像采集面逆衍射计算到待测面的复振幅,j是虚数,*代表取复共轭,Im()代表取虚部,Re()代表取实部,和分别代表衍射计算算子和逆衍射计算算子;
其中,β是权重常数,一般取0.75;
S5:计算下一次迭代的输入复振幅
S6:如果i<N,i=i+1,并返回S4,否则结束迭代,获得满足衍射光斑约束的复振幅估计值;
S7:对S6获得的复振幅估计值的振幅和相位分别使用去噪算法处理,移除噪声影响,并截取支持域之内的复振幅估计值得到真实的待测平板复振幅。
进一步地,所述S4中的衍射计算算子优先使用角谱衍射理论计算。
进一步地,所述S7中的去噪算法优先使用BM3D去噪算法。
进一步地,S4中的梯度优化的方法优选梯度下降法。
本发明的有益效果如下:
本发明的带有负反馈调节的无透镜成像方法利用一幅衍射光斑图实现无透镜待测平板的相位和振幅的重建,通过在支持域之外叠加负反馈实现了算法的快速收敛,并克服了传统相位恢复算法中孪生像的问题;支持域内使用梯度优化算法提高了算法灵活度。另外,本发明的方法使用单图成像,克服了传统无透镜成像方法中多图成像中轴向位置测不准以及横向位置不匹配造成的误差,且本发明的一幅衍射光斑图实现波前重建能够适用于实时动态成像。
附图说明
图1是本发明的带有负反馈调节的无透镜成像的装置示意图;
图2是本发明的带有负反馈调节的无透镜成像方法的流程图;
图3是本发明的带有负反馈调节的无透镜成像方法与传统方法的结果对比图,其中,图(a1)和(a2)分别为本发明所提方法的恢复的振幅和相位,图(b1)和(b2)分别为多图无透镜成像方法恢复的振幅和相位,图(c1)和(c2)分别为真实图的振幅和相位。
具体实施方式
下面根据附图和优选实施例详细描述本发明,本发明的目的和效果将变得更加明白,应当理解,此处所描述的具体实施例仅仅用以解释本发明,并不用于限定本发明。
如图1所示,本发明的带有负反馈调节的无透镜成像方法的装置,在激光器1的出射光路的光轴上依次布置有扩束器2、待测平板3、图像传感器4,该方法包括如下步骤:
S1:使用图像传感器4,在某一位置上采集一幅含有待测平板信息的衍射光斑;
S2:分别设置待测平板3的支持域D、待测平板3与图像传感器4之间的距离d、负反馈权重系数β、相位优化步长stepphase、振幅优化步长stepamp、待测平板重建的迭代总数N、波前检测的初始迭代计数i=1、待测平板重建复振幅初始估计值g1(x,y),其中,(x,y)为待测平板3的坐标;
S3:构建优化目标函数E
其中,I是在d位置处采集的衍射光斑,U是待测平板在位置d处的衍射积分;W(u,v)是权重分布,用来移除图像传感器4上信噪比低的像素点和坏点;
其中,gi(x,y)表示第i次迭代的待测平板的复振幅输入估计值,Gi(x,y)代表从待测面衍射计算到图像采集面的复振幅,代表从图像采集面逆衍射计算到待测面的复振幅,j是虚数,*代表取复共轭,Im()代表取虚部,Re()代表取实部,和分别代表衍射计算算子和逆衍射计算算子;
其中,β是权重常数,一般取0.75;
S5:计算下一次迭代的输入复振幅
S6:如果i<N,i=i+1,并返回S4,否则结束迭代,获得满足衍射光斑约束的复振幅估计值;
S7:对S6获得的复振幅估计值的振幅和相位分别使用去噪算法处理,移除噪声影响,并截取支持域之内的复振幅估计值得到真实的待测平板复振幅。
考虑到无透镜成像模型及角谱衍射理论计算精度最高,所述S4和S6中的衍射计算算子优先使用角谱衍射理论计算。
考虑到去噪效果好及精度高,所述S8中的去噪算法优先使用BM3D去噪算法。
为了保证梯度优化的效果最好,因此S4中的梯度优化优选梯度下降法。
下面给出本发明的方法的一个具体实施例,对该方法的技术效果进行说明
这里选择,d=200mm,口径D=5mm,相位优化步长stepphase=-0.4、振幅优化步长stepamp=0.01、相位成像的迭代总数N=5000,β=0.75。
在该实施例中,采集一幅图像进行波前的相位和振幅重建,所述的图像传感器与待测平板之间的距离为200mm,选用的衍射计算模型为角谱衍射模型,图3(a1)和(a2)分别为本发明所提方法的恢复的振幅和相位,图3(b1)和(b2)分别为多图无透镜成像方法恢复的振幅和相位。从图中可以看出,本发明提出的方法用一幅图即可精准恢复波前的振幅和相位,恢复结果与真实值形貌轮廓一致,传统的多图无透镜成像方法恢复的振幅和相位精度较差。
本领域普通技术人员可以理解,以上所述仅为发明的优选实例而已,并不用于限制发明,尽管参照前述实例对发明进行了详细的说明,对于本领域的技术人员来说,其依然可以对前述各实例记载的技术方案进行修改,或者对其中部分技术特征进行等同替换。凡在发明的精神和原则之内,所做的修改、等同替换等均应包含在发明的保护范围之内。
Claims (4)
1.一种带有负反馈调节的无透镜成像方法,该方法在激光器的出射光路的光轴上依次布置有扩束器、待测平板、图像传感器,其特征在于,该方法包括如下步骤:
S1:使用所述图像传感器在某一位置上采集一幅含有所述的待测平板信息的衍射光斑;
S2:分别设置所述待测平板的支持域D、所述待测平板与所述图像传感器之间的距离d、负反馈权重系数β、相位优化步长stepphase、振幅优化步长stepamp、待测平板重建的迭代总数N、波前检测的初始迭代计数i=1、待测平板重建复振幅初始估计值g1(x,y),其中,(x,y)为所述待测平板的坐标;
S3:构建优化目标函数E
其中,I是在d位置处采集的衍射光斑,U是待测平板在位置d处的衍射积分;W(u,v)是权重分布,用来移除图像传感器上信噪比低的像素点和坏点。
S4:对所述待测平板的支持域D内的部分,首先采用公式(2)~(5)计算波前优化的相位梯度和振幅梯度Δa,然后更新待测波前的振幅和相位,并通过公式(6)对估计的波前进行梯度优化,得到所述支持域D内的复振幅
其中,gi(x,y)表示第i次迭代的待测平板的复振幅输入估计值,Gi(x,y)代表从待测面衍射计算到图像采集面的复振幅,代表从图像采集面逆衍射计算到待测面的复振幅,j是虚数,*代表取复共轭,Im()代表取虚部,Re()代表取实部,和分别代表衍射计算算子和逆衍射计算算子;
其中,β是权重常数,一般取0.75;
S5:计算下一次迭代的输入复振幅
S6:如果i<N,i=i+1,并返回S4,否则结束迭代,获得满足衍射光斑约束的复振幅估计值;
S7:对S6获得的复振幅估计值的振幅和相位分别使用去噪算法处理,移除噪声影响,并截取支持域之内的复振幅估计值得到真实的待测平板复振幅。
2.根据权利要求1所述的带有负反馈调节的无透镜成像方法,其特征在于,所述S4中的衍射计算算子优先使用角谱衍射理论计算。
3.根据权利要求1所述的带有负反馈调节的无透镜成像方法,其特征在于,所述S7中的去噪算法优先使用BM3D去噪算法。
4.根据权利要求1所述的带有负反馈调节的无透镜成像方法,其特征在于,S4中的梯度优化的方法优选梯度下降法。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110459324.1A CN113189101B (zh) | 2021-04-27 | 2021-04-27 | 一种带有负反馈调节的无透镜成像方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110459324.1A CN113189101B (zh) | 2021-04-27 | 2021-04-27 | 一种带有负反馈调节的无透镜成像方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113189101A true CN113189101A (zh) | 2021-07-30 |
CN113189101B CN113189101B (zh) | 2024-01-30 |
Family
ID=76979375
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110459324.1A Active CN113189101B (zh) | 2021-04-27 | 2021-04-27 | 一种带有负反馈调节的无透镜成像方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113189101B (zh) |
Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102721477A (zh) * | 2012-06-04 | 2012-10-10 | 中国科学院光电技术研究所 | 一种用于相位差波前传感器的离焦量测量方法 |
CN105954994A (zh) * | 2016-06-30 | 2016-09-21 | 深圳先进技术研究院 | 无透镜数字全息显微成像的图像增强方法 |
CN106289543A (zh) * | 2016-07-25 | 2017-01-04 | 南京理工大学 | 大动态离焦范围单帧双幅图像光强采集装置及其采集方法 |
CN107101974A (zh) * | 2017-07-03 | 2017-08-29 | 西京学院 | 一种新型三步无透镜相干衍射成像方法 |
CN108204949A (zh) * | 2017-12-28 | 2018-06-26 | 南京理工大学 | 基于反射式交叠衍射成像的相位恢复算法 |
CN110160751A (zh) * | 2019-05-16 | 2019-08-23 | 浙江大学 | 一种基于相位恢复的宽频段波前误差检测装置及检测方法 |
CN110470245A (zh) * | 2019-08-21 | 2019-11-19 | 浙江大学 | 一种基于菲涅尔波带片衍射信息融合的相位恢复检测装置及相位恢复方法 |
US20200003709A1 (en) * | 2017-10-09 | 2020-01-02 | Nanjing University | Imaging device, imaging method, and imaging system |
US20200012234A1 (en) * | 2017-03-27 | 2020-01-09 | Ke Li | Method for generating a holographic image, processor, holographic image display device, and apparatus |
CN112697751A (zh) * | 2020-12-07 | 2021-04-23 | 杭州电子科技大学 | 多角度照明无透镜成像方法、系统及装置 |
-
2021
- 2021-04-27 CN CN202110459324.1A patent/CN113189101B/zh active Active
Patent Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102721477A (zh) * | 2012-06-04 | 2012-10-10 | 中国科学院光电技术研究所 | 一种用于相位差波前传感器的离焦量测量方法 |
CN105954994A (zh) * | 2016-06-30 | 2016-09-21 | 深圳先进技术研究院 | 无透镜数字全息显微成像的图像增强方法 |
CN106289543A (zh) * | 2016-07-25 | 2017-01-04 | 南京理工大学 | 大动态离焦范围单帧双幅图像光强采集装置及其采集方法 |
US20200012234A1 (en) * | 2017-03-27 | 2020-01-09 | Ke Li | Method for generating a holographic image, processor, holographic image display device, and apparatus |
CN107101974A (zh) * | 2017-07-03 | 2017-08-29 | 西京学院 | 一种新型三步无透镜相干衍射成像方法 |
US20200003709A1 (en) * | 2017-10-09 | 2020-01-02 | Nanjing University | Imaging device, imaging method, and imaging system |
CN108204949A (zh) * | 2017-12-28 | 2018-06-26 | 南京理工大学 | 基于反射式交叠衍射成像的相位恢复算法 |
CN110160751A (zh) * | 2019-05-16 | 2019-08-23 | 浙江大学 | 一种基于相位恢复的宽频段波前误差检测装置及检测方法 |
CN110470245A (zh) * | 2019-08-21 | 2019-11-19 | 浙江大学 | 一种基于菲涅尔波带片衍射信息融合的相位恢复检测装置及相位恢复方法 |
CN112697751A (zh) * | 2020-12-07 | 2021-04-23 | 杭州电子科技大学 | 多角度照明无透镜成像方法、系统及装置 |
Non-Patent Citations (1)
Title |
---|
王智博;张骏;李拓;史诗;: "基于已知振幅信息约束的叠层光电成像", 中国科学院大学学报, no. 05 * |
Also Published As
Publication number | Publication date |
---|---|
CN113189101B (zh) | 2024-01-30 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107533759B (zh) | 对于关注图案图像群体的异常检测 | |
CN107917676B (zh) | 一种基于条纹图像频谱分析的干涉测量方法 | |
US20060048089A1 (en) | System and method for simulating an aerial image | |
TWI646303B (zh) | 決定裝置之儀器轉移函數之方法以及系統 | |
CN107564091A (zh) | 一种基于快速对应点搜索的三维重建方法及装置 | |
CN106324995B (zh) | 光刻机原位快速高空间分辨率波像差检测装置及方法 | |
CN109870129A (zh) | 一种基于相位偏折原理的晶圆表面粗糙度检测装置 | |
CN114460093B (zh) | 一种航空发动机缺陷检测方法及系统 | |
CN102252763A (zh) | 光学成像系统光学波前的单焦面高精度测试方法 | |
JP5595463B2 (ja) | 波面光学測定装置 | |
CN111561877B (zh) | 一种基于点衍射干涉仪的可变分辨率相位解包裹方法 | |
CN110895792B (zh) | 一种图像拼接方法及装置 | |
CN112903121B (zh) | 一种基于交叉迭代自动对焦的波前检测方法 | |
CN113188671B (zh) | 一种基于交叉迭代自动位置矫正的波前检测方法 | |
WO2016201788A1 (zh) | 光刻机原位多通道成像质量检测装置及方法 | |
CN112504240A (zh) | 一种激光投线仪校准系统及校准方法 | |
Alici | Extraction of modulation transfer function by using simulated satellite images | |
Rodríguez | Online self-calibration for mobile vision based on laser imaging and computer algorithms | |
CN113189101A (zh) | 一种带有负反馈调节的无透镜成像方法 | |
Bai et al. | Recent Progress of Full-Field Three-Dimensional Shape Measurement Based on Phase Information | |
JP2021051038A (ja) | 収差推定方法、収差推定装置、プログラムおよび記録媒体 | |
CN113432731B (zh) | 一种光栅横向剪切干涉波前重建过程中的补偿方法 | |
CN111369610B (zh) | 基于可信度信息的点云数据粗差定位和剔除方法 | |
JP6539812B1 (ja) | レンズの検査方法 | |
CN113532473B (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 |