CN101726366B - 一种基于多针孔板的波前测量方法和装置 - Google Patents

一种基于多针孔板的波前测量方法和装置 Download PDF

Info

Publication number
CN101726366B
CN101726366B CN2009102313226A CN200910231322A CN101726366B CN 101726366 B CN101726366 B CN 101726366B CN 2009102313226 A CN2009102313226 A CN 2009102313226A CN 200910231322 A CN200910231322 A CN 200910231322A CN 101726366 B CN101726366 B CN 101726366B
Authority
CN
China
Prior art keywords
spininess
orifice plate
pin hole
fourier transform
measurement
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.)
Expired - Fee Related
Application number
CN2009102313226A
Other languages
English (en)
Other versions
CN101726366A (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.)
Shandong Normal University
Original Assignee
Shandong Normal University
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 Shandong Normal University filed Critical Shandong Normal University
Priority to CN2009102313226A priority Critical patent/CN101726366B/zh
Publication of CN101726366A publication Critical patent/CN101726366A/zh
Application granted granted Critical
Publication of CN101726366B publication Critical patent/CN101726366B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Length Measuring Devices By Optical Means (AREA)
  • Investigating Materials By The Use Of Optical Means Adapted For Particular Applications (AREA)
  • Investigating Or Analysing Materials By Optical Means (AREA)

Abstract

本发明公开了一种基于多针孔板的波前测量方法和装置。它将被测量物波照射到含有一个参考针孔和多个测量针孔的多针孔板上;用图像传感器记录透过多针孔板的物波的夫琅和费衍射光场的强度分布图样;该夫琅和费衍射光场正比于透过多针孔板的物波的傅立叶变换;对所记录的强度分布图样做逆傅里叶变换,将该逆傅里叶变换图样中与多针孔板上各测量针孔的中心位置所对应的点的函数值提取出来,该函数值正比于被测物波在各测量针孔处的复振幅值;利用该复振幅值可在计算机中重现被成像物体。本发明提高了衍射成像速度,结构简单,调节方便,成本较低,适用于多种不同光源,可实现纯相位物体、复数物体或三维物体的实时波前传感和衍射成像。

Description

一种基于多针孔板的波前测量方法和装置
技术领域
本发明涉及一种测量入射光波波前或复振幅分布的装置,特别是一种利用多针孔板实现对入射光波波前或复振幅分布进行记录和重现的基于多针孔板的波前测量方法和装置。
背景技术
在自适应光学、光学精密测量、显微成像等应用领域,经常需要测量光波的波前或复振幅分布。尤其是在自适应光学系统中,为了实现对波前的动态控制,需要实时测量波前信息。目前,已经发展了多种测量光波的波前或复振幅分布的技术,如全息干涉测量技术、哈特曼波前传感技术、相位反演波前测量技术等。
全息干涉测量技术是利用被测光波和与之相干的参考波的干涉实现光波波前的记录,通过对干涉图样的分析实现光波波前的测量和重现。该技术需要使用分光装置产生所需要的相干参考波。虽然同轴全息干涉测量技术可以避免使用分光装置,但需要克服背景光和共轭项的影响。
哈特曼波前传感技术的基本原理是利用阵列微透镜把入射光波分割为若干子光波,然后利用二维阵列图像传感器(如CCD)测量各个子光波经微透镜的聚焦光斑的质心漂移量得到波前斜率,然后用各种波前复原算法得到被测光波的波前相位。哈特曼波前传感器结构简单,并且可以实现波前的实时动态测量。因此,该技术已经得到广泛应用。但是微透镜阵列的使用使哈特曼传感器的分辨率和应用领域受到限制。因为在有些波段中,如电子波、X光和极紫外光波段,制备所需要的微透镜阵列是非常困难的。
相位反演波前测量技术则是利用物波的衍射光场或像的强度分布实现光波波前或复振幅分布的测量。由于可以避免成像透镜孔径和像差对分辨率的限制,这种技术特别适用于像X射线、电子束这类缺乏或很难制备出高质量成像透镜的领域。相位反演波前测量技术的关键是怎样从一幅或多幅衍射强度图样中准确快速地将被测物波的振幅和相位信息恢复出来。目前,解决这一问题的一个主要途径是采用迭代算法。传统的迭代方法一般都需要较长的迭代时间,迭代结果存在不确定性;对被测物体也存在一些苛刻的限制条件,如物体要特殊固定、要求被测物体是纯振幅物体或纯相位物体,等等。
申请人在先申请的中国专利申请号“200910014336.2”中提供了一种基于多针孔板的相干衍射成像方法及其处理装置,它是让被测量波前先通过一个可旋转的多针孔板,然后再用图像传感器记录物波经过多针孔板后的夫琅和费衍射强度分布,最后借助计算机图像处理技术提取出被测物波波前的振幅和相位分布信息。该技术的一个不足之处是,为了实现二维波前的测量需要旋转多针孔板进行多次记录,影响了波前测量速度,不利于在自适应光学等需要对波前进行高速动态测量领域中的应用。
发明内容
本发明的目的在于克服现有技术的不足,在申请号为“200910014336.2”的中国专利的基础上,提供一种不需要旋转针孔板进行多次记录就能实现二维波前测量的基于多针孔板的波前测量方法和装置。
本发明的目的可通过以下技术措施实现:
一种基于多针孔板的波前测量方法,包括以下步骤:
第一步,将被测量物波照射到含有一个参考针孔和多个测量针孔的多针孔板上;
第二步,用图像传感器记录透过多针孔板的物波的夫琅和费衍射光场的强度分布图样;记录平面上的衍射光场正比于多针孔板平面上透过多针孔板的物波的傅立叶变换;
第三步,对该衍射强度图样做逆傅里叶变换,将该逆傅里叶变换图样中与多针孔板上各测量针孔的中心位置所对应的点的函数值提取出来,该函数值正比于被测物波在各测量针孔处的复振幅值;利用该复振幅值可在计算机中重现物波波前和成像。
一种基于多针孔板的波前测量方法的处理装置,它包括依次设置的多针孔板、图像传感器和处理装置,图像传感器与多孔板之间的距离Z满足夫琅和费衍射条件,即图像传感器记录平面上的光场分布正比于透过多针孔板的物波的傅立叶变换。
一种基于多针孔板的波前测量方法的处理装置,它包括依次设置的多针孔板、傅立叶变换透镜、图像传感器和处理装置,所述图像传感器位于傅里叶变换透镜的后焦面上,即图像传感器记录平面上的光场分布正比于透过多针孔板的物波的傅立叶变换。
所述多针孔板是一块不透明薄板上设置一个参考针孔和若干测量针孔组成,所有测量针孔按照正三角形或称六方形二维周期点阵格式均匀排列,相邻测量针孔中心的间隔L等于或大于测量针孔直径D的2.5倍;参考针孔位于由相邻的三个测量针孔构成的一个正三角形的几何中心处。
所述多针孔板上的测量针孔排列成两个相互错位半个周期的正方形阵列,阵列周期T等于或大于针孔直径D的2倍;该两个正方形阵列分别位于参考针孔两侧,参考针孔中心与最近邻的测量针孔中心之间的距离等于阵列周期T的一半。
所述多针孔板上的参考针孔和测量针孔形状为圆形或正方形或其它形状。
本发明的原理是:被测量光波先经过多针孔板抽样,然后由图像传感器记录下透过多针孔板的光波的夫琅和费衍射的强度图样IP。设透过多针孔板上第m个测量针孔的物波的复振幅值为其中Am
Figure G2009102313226D00032
分别为透过第m个测量针孔的物波的振幅和相位,
Figure G2009102313226D00033
为针孔板所在平面上的位置矢量(以参考针孔的中心为坐标原点),
Figure G2009102313226D00034
为第m个测量针孔中心的位置矢量,j为虚数,
Figure G2009102313226D00035
为第m个测量针孔的孔径函数;在记录平面上得到的夫琅和费衍射光场的强度分布可表示为:
其中,I0为积分常数,
Figure G2009102313226D00037
为傅立叶变换算符,M为针孔板上测量针孔的数目,m=0对应参考针孔,
Figure G2009102313226D00038
为记录平面上的空间频率坐标。对该衍射强度图样做逆傅立叶变换得到被多针孔板抽样的光波波前的自相关函数
Figure G2009102313226D00039
该自相关函数可以表示成以下形式(见Gong-Xiang Wei,Cheng-Shan Guo,et al.,Appl.Opt.48,5099,2009),即
g ( r → ) = C 0 A 0 2 δ ( r → )
Figure G2009102313226D000312
Figure G2009102313226D000313
Figure G2009102313226D000314
其中,A0
Figure G2009102313226D000315
分别为透过参考针孔的物波的振幅和相位值,C0为积分常数。由(2)式可见,该自相关函数可分成四项,其中第二项就是待测量的物波波前的抽样阵列。一般情况下,(2)式中的第二项和其他各项在空间上是相互重叠的,需要采用复杂的迭代算法才有可能将其提取出来。由于本发明中所采用的多针孔板上的抽样针孔的位置矢量满足以下条件
r → m - r → n ≠ r → k and r → m + r → n ≠ 0 , ( m , n , k = 1,2 , . . . , M ) , - - - ( 3 )
使得(2)式中的第二项与其它各项在空间上是相互分离的。这样,只需一个简单的抽样函数 P ( r → ) = Σ m = 1 M δ ( r → - r → m ) 就可以将(2)式中的第二项直接提取出来,即
Figure G2009102313226D000318
利用(4)式得到的被测物波波前的抽样阵列就可以在计算机中实现被测物体的数字成像或重现。
本发明相对于在先申请的中国专利200910014336.2避免了旋转和同步扫描装置,仅从一幅强度图样中就可以将待测波前恢复出来,其意义是可以实现复振幅物体的动态实时测量,在自适应光学领域和超快成像领域有重要应用价值。
本发明与现有技术相比有以下优点:利用本发明所公开的基于多针孔板的波前测量装置进行二维波前测量时不需要旋转多针孔板和进行多次测量,从而真正实现了二维波前的动态和实时测量。本发明所公开的基于多针孔板的波前测量装置完全避免了精密旋转和同步控制元件的使用,简化了记录过程,大大降低了系统成本。本发明所采用的多针孔板容易制备,可适用于像X射线、电子波这类难以使用传统的哈特曼波前传感器进行动态波前测量的领域。
附图说明
图1是本发明的波前测量方法的过程方框图;
图2是本发明的波前测量方法中的多针孔板的第一种设计实例的针孔分布示意图;
图3是本发明的多针孔板的第二种设计实例的针孔分布示意图;
图4是本发明的波前测量处理装置的第一种实施方式的结构示意图;
图5是本发明的波前测量处理装置的第二种实施方式的结构示意图;
图6a是采用本发明的波前测量方法和装置进行相干衍射成像的一个实验结果实例;
图6b是采用本发明的波前测量方法和装置进行相干衍射成像的一个实验结果实例;
图6c是采用本发明的波前测量方法和装置进行相干衍射成像的一个实验结果实例。
其中,1.多针孔板、2.参考针孔、3.测量针孔、4.图像传感器、5.计算机、6.傅里叶变换透镜。
具体实施方式
下面结合附图与实施例对本发明做进一步的描述。
如图1所示,本发明方法包括如下步骤:(1)首先用一个多针孔板1对被测量物波的波前进行抽样,该多针孔板1上的针孔由一个参考针孔2和若干个测量针孔3构成;(2)透过多针孔板1的物波继续传播,在远场(或通过一个傅里叶变换透镜后在透镜后焦面)形成夫琅和费衍射光场;(3)用一个二维图像传感器4(如CCD)记录透过多针孔板1的物波的夫琅和费衍射强度图样;(4)利用图像处理技术对该夫琅和费衍射强度图样作逆傅里叶变换得到透过多针孔板1的物波复振幅的相关函数图样;将该相关函数图样中与多针孔板1上各测量针孔3的中心位置所对应的点的函数值提取出来,就得到透过多针孔板1的物波复振幅在各测量针孔3处的相对振幅和位相值;利用该测量数据就可以在计算机中通过程序模拟光的衍射过程,重现被成像物体。
图2是本发明的多针孔板1的一种典型实施方式。该多针孔板1上的针孔由直径均为D的一个参考针孔2和多个测量针孔3构成,测量针孔3按照正三角形(或六方形)二维周期点阵格式均匀排列,相邻测量针孔3中心的间隔L等于或大于测量针孔3直径D的2.5倍;参考针孔2位于由相邻的三个测量针孔3构成的一个正三角形的几何中心处。参考针孔2和测量针孔3的形状可以是圆形或正方形或其它形状。
图3是本发明的方法中所述的多针孔板1的第二种实施方式。该多针孔板1上的测量针孔3排列成两个相互错位半个周期的正方形阵列,阵列周期T等于或大于针孔直径D的2倍;该两个正方形阵列分别位于参考针孔2两侧,参考针孔2中心与最近邻的测量针孔3中心之间的距离等于阵列周期T的一半。参考针孔2和测量针孔3的形状可以是圆形或正方形或其它形状。
图4是本发明的波前测量处理装置的一种典型实施方式。该装置包括多针孔板1、图像传感器4、计算机5。多针孔板1到图像传感器4的感光平面之间的距离Z保证满足夫琅和费衍射条件,使透过多针孔板1的物波在图像传感器感光平面上的衍射光场正比于该物波的傅里叶变换。图像传感器4由计算机5控制实现对衍射强度图像的数字记录。从所记录的衍射强度图样中提取物波的复振幅信息并成像由根据本发明方法设计的计算机程序完成。
图5是本发明的装置的第二种实施方式,比图4多了一个傅里叶变换透镜6,即在多针孔板1和图像传感器4之间放置傅里叶变换透镜6,图像传感器记录平面位于傅里叶变换透镜6的后焦面上。这样可以使装置紧凑,并且可以通过改变透镜6的焦距方便地放大或缩小传播到传感器记录平面上的夫琅和费衍射强度图样的大小。
图6a-图6c是采用本发明的波前测量方法和装置进行波前测量和相干衍射成像的一个实验结果实例。实验中,光源采用He-Ne激光,输出光波的波长为0.6328微米;被成像物体为字符“山”的缩微图片,如图6a所示;多针孔板采用图2所示的六方形阵列模式。针孔直径为50微米,相邻测量针孔的相距为400微米。实验装置采用图5所示的第二种实施方式,其中傅里叶变换透镜的焦距为240毫米,图像传感器为CCD数码相机,像素数为1300×1030,像素大小为6.7微米。图6b为用CCD记录的透过多针孔板的被测物波的夫琅和费衍射强度图样。图6c为利用图6b所示的强度图样在计算机中通过计算机程序得到的被测物体的波前重现像。
上述方法及实施例都是通过记录被测物波通过一个多针孔板后的夫琅和费衍射强度图样和通过对该图样的图像处理来达到重现物波的振幅和相位分布的目的。本发明的实施不局限于上述具体实施方案。只要是通过记录被测物波经过上述方法所述的多针孔板后的夫琅和费衍射强度图样来实现重现物波的振幅和相位分布目的的方法、装置与系统,均属于本发明的保护范围。

Claims (6)

1.一种基于多针孔板的波前测量方法,其特征在于,包括以下步骤:
第一步,将被测量物波照射到含有一个参考针孔和多个测量针孔的多针孔板上;
第二步,用图像传感器记录透过多针孔板的物波的夫琅和费衍射光场的强度分布图样;记录平面上的衍射光场正比于多针孔板平面上透过多针孔板的物波的傅立叶变换;
第三步,对该衍射强度分布图样做逆傅里叶变换,将该逆傅里叶变换图样中与多针孔板上各测量针孔的中心位置所对应的点的复振幅值提取出来,该复振幅值正比于被测物波在各测量针孔处的复振幅值;利用该复振幅值可在计算机中重现被成像物体。
2.一种权利要求1所述基于多针孔板的波前测量方法的处理装置,其特征在于,它包括依次设置的多针孔板、图像传感器和计算机,图像传感器与多孔板之间的距离Z满足夫琅和费衍射条件,即图像传感器记录平面上的光场分布正比于透过多针孔板的物波的傅立叶变换。
3.一种权利要求1所述基于多针孔板的波前测量方法的处理装置,其特征在于,它包括依次设置的多针孔板、傅立叶变换透镜、图像传感器和计算机,所述图像传感器位于傅里叶变换透镜的后焦面上,即图像传感器记录平面上的光场分布正比于透过多针孔板的物波的傅立叶变换。
4.如权利要求2或3所述的基于多针孔板的波前测量方法的处理装置,其特征在于,所述多针孔板是一块不透明薄板上设置一个参考针孔和若干测量针孔组成,所有测量针孔按照正三角形或六方形二维周期点阵格式均匀排列,相邻测量针孔中心的间隔L等于或大于测量针孔直径D的2.5倍;参考针孔位于由相邻的三个测量针孔构成的一个正三角形的几何中心处。
5.如权利要求2或3所述的基于多针孔板的波前测量方法的处理装置,其特征在于,所述多针孔板上的测量针孔排列成两个相互错位半个周期的正方形阵列,阵列周期T等于或大于针孔直径D的2倍;该两个正方形阵列分别位于参考针孔两侧,参考针孔中心与最近邻的测量针孔中心之间的距离等于阵列周期T的一半。
6.如权利要求2或3所述的基于多针孔板的波前测量方法的处理装置,其特征在于,所述多针孔板上的参考针孔和测量针孔形状为圆形或正方形或其它形状。
CN2009102313226A 2009-12-02 2009-12-02 一种基于多针孔板的波前测量方法和装置 Expired - Fee Related CN101726366B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN2009102313226A CN101726366B (zh) 2009-12-02 2009-12-02 一种基于多针孔板的波前测量方法和装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN2009102313226A CN101726366B (zh) 2009-12-02 2009-12-02 一种基于多针孔板的波前测量方法和装置

Publications (2)

Publication Number Publication Date
CN101726366A CN101726366A (zh) 2010-06-09
CN101726366B true CN101726366B (zh) 2011-11-23

Family

ID=42447589

Family Applications (1)

Application Number Title Priority Date Filing Date
CN2009102313226A Expired - Fee Related CN101726366B (zh) 2009-12-02 2009-12-02 一种基于多针孔板的波前测量方法和装置

Country Status (1)

Country Link
CN (1) CN101726366B (zh)

Families Citing this family (14)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101968383B (zh) * 2010-09-02 2012-01-25 北京理工大学 一种抗扰动的时频域波前检测方法
CN102169020B (zh) * 2010-12-08 2012-03-21 中国科学院光电技术研究所 一种哈特曼波前传感器质心运算的流水线装置
CN103292910B (zh) * 2012-02-29 2017-05-24 深圳光启高等理工研究院 基于超材料的哈特曼波前传感器
US11468557B2 (en) 2014-03-13 2022-10-11 California Institute Of Technology Free orientation fourier camera
US9829695B2 (en) * 2015-01-26 2017-11-28 California Institute Of Technology Array level Fourier ptychographic imaging
US11092795B2 (en) 2016-06-10 2021-08-17 California Institute Of Technology Systems and methods for coded-aperture-based correction of aberration obtained from Fourier ptychography
EP3327413B1 (en) * 2016-11-24 2022-04-27 IMEC vzw A method, an apparatus and a system for holographic wavefront sensing
CN108375421A (zh) * 2018-03-20 2018-08-07 深圳大学 一种激光波前测量仪及激光波前测量方法
CN110187627B (zh) * 2019-06-28 2021-04-02 赣南师范大学 一种复振幅传感成像装置和方法
CN112504479B (zh) * 2020-12-02 2022-05-03 中国人民解放军国防科技大学 一种新型波前测量装置及方法、设备及介质
CN112504480B (zh) * 2020-12-02 2021-12-17 中国人民解放军国防科技大学 一种高精度大视场波前测量方法、电子设备及介质
CN112414565B (zh) * 2020-12-02 2021-09-21 中国人民解放军国防科技大学 一种大视场波前测量装置及方法、设备及介质
CN112629679B (zh) * 2020-12-02 2021-11-30 中国人民解放军国防科技大学 一种适应于背景纹影的高精度测量方法、电子设备及介质
CN113325596B (zh) * 2021-06-04 2023-04-04 江西欧迈斯微电子有限公司 一种衍射光学元件及光学设备

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5294971A (en) * 1990-02-07 1994-03-15 Leica Heerbrugg Ag Wave front sensor
CN1245904A (zh) * 1998-08-26 2000-03-01 中国科学院光电技术研究所 光学波前传感器
CN1385675A (zh) * 2001-05-16 2002-12-18 中国科学院光电技术研究所 波前传感器

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5294971A (en) * 1990-02-07 1994-03-15 Leica Heerbrugg Ag Wave front sensor
CN1245904A (zh) * 1998-08-26 2000-03-01 中国科学院光电技术研究所 光学波前传感器
CN1385675A (zh) * 2001-05-16 2002-12-18 中国科学院光电技术研究所 波前传感器

Also Published As

Publication number Publication date
CN101726366A (zh) 2010-06-09

Similar Documents

Publication Publication Date Title
CN101726366B (zh) 一种基于多针孔板的波前测量方法和装置
CN102589414B (zh) 可实时测量的同步相移斐索干涉装置
EP3677969B1 (en) Holographic imaging device and data processing method therefor
EP3065001B1 (en) Holographic microscope and data processing method for high-resolution hologram image
CN105423911A (zh) 一种基于光栅离焦的共路数字全息显微装置与方法
CN106441085A (zh) 一种双载频共路数字全息显微装置及显微方法
CN102865811B (zh) 基于正交双光栅的同步移相共光路干涉显微检测装置及检测方法
CN103592108A (zh) Ccd芯片调制传递函数测试装置及方法
CN104833977A (zh) 一种基于微波片阵列的瞬时遥感偏振成像装置及其实现方法
CN109000781B (zh) 一种结构微振动线域测量装置及方法
WO2009113068A1 (en) Intraoral imaging system and method based on conoscopic holography
CN103913127A (zh) 一种基于子孔径相位拼接的数字全息球面面型检测装置
CN103425035A (zh) 一种采用异型掩模板的全息体视图打印系统
CN105242512A (zh) 基于远心光学结构的透射式数字全息显微成像装置
CN107356195B (zh) 基于二维周期光栅和点衍射的三视场数字全息检测装置与方法
CN101957171B (zh) 一种可有效抑制零级和共轭像的同轴数字全息方法
CN201885805U (zh) 环形共光路点衍射干涉波前传感装置
CN105823563B (zh) 动态高分辨率波前测量装置及方法
CN102865810B (zh) 基于正交双光栅的同步相移共光路干涉检测装置及检测方法
CN102878930B (zh) 一种位相物体位相分布的定量测量方法和装置及其应用
CN107421437B (zh) 基于二维相位光栅和点衍射的三视场数字全息检测装置与方法
CN101532908A (zh) 离散成像器件统计调制传递函数激光散斑测量方法和装置
CN103411687A (zh) 基于三角数字全息测量空间相干性的系统与方法
CN107356194B (zh) 基于二维周期光栅和点衍射的四视场数字全息检测装置与方法
CN104634253A (zh) 基于形貌相关计算的三维位移测量方法

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20111123

Termination date: 20141202

EXPY Termination of patent right or utility model