CN111650203A - 一种微球内表面缺陷测量方法 - Google Patents
一种微球内表面缺陷测量方法 Download PDFInfo
- Publication number
- CN111650203A CN111650203A CN202010335374.4A CN202010335374A CN111650203A CN 111650203 A CN111650203 A CN 111650203A CN 202010335374 A CN202010335374 A CN 202010335374A CN 111650203 A CN111650203 A CN 111650203A
- Authority
- CN
- China
- Prior art keywords
- light
- complex amplitude
- microsphere
- diffraction
- defects
- 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
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
- G01N21/88—Investigating the presence of flaws or contamination
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)
- Investigating Materials By The Use Of Optical Means Adapted For Particular Applications (AREA)
Abstract
本发明公开了一种微球内表面缺陷测量方法,获取微球外表面测试光与参考光的干涉强度图,微球内表面测试光和参考光的干涉强度图以及参考光的光强分布图;获得外表面干涉光场和内表面干涉光场的复振幅分布;引入虚拟透镜进行复振幅重构获得重构后的外表面干涉光场和内表面干涉光场的复振幅分布;将外表面测试光的复振幅转化成透射光重构复振幅;获得消去调制后的内表面重构复振幅;根据衍射距离,采用逆向衍射计算获取衍射面复振幅分布,确定内表面缺陷相位信息。本发明引入了虚拟透镜,将曲面间的衍射计算转换成平面与平面间的衍射计算,避免了球面相位的引入所带来的采样问题。
Description
技术领域
本发明属于光学检测领域,具体为一种微球内表面缺陷测量方法。
背景技术
在ICF实验中,多束高能激光直接或者间接轰击微型靶丸产生聚变反应。直径介于几百微米至几毫米的微型靶丸是激光聚变的核心器件之一,其靶丸的内表面任何微小缺陷都有可能引发瑞利—泰勒不稳定性,进而引起不对称压缩,导致实验效率降低,甚至打靶失败的严重后果。
在现有显微系统的内表面缺陷测量中,根据表面检测方法可以分为:点探测法、面探测法。在点探测法中在国外以里弗莫尔国家实验室为代表,通过在传统的AFM探头的位置加入光纤干涉壁厚测量传感器。这使得能够通过沿相同轮廓截面测量外表面轮廓和壳层厚度,从而间接获得靶丸的内表面轮廓,然而在实际测量和对准期间,难以确保AFM探针和光纤干涉壁厚测量传感器探头具有完全重叠的测量轨迹,并且所使用光纤白光干涉仪的光斑直径大约为100um,因此,它的横向分辨率很低,不能满足靶丸内表面轮廓高精度测量要求。在国内北京理工大学提出使用激光差动共聚焦的方法,通过显微物镜径向扫面分别位于ICF靶丸的内表面和外表面,并使用内外表面测量位置参数和ICF靶丸的折射率通过光线追迹获得内表面的轮廓信息,可以实现0.166um的横向分辨率和10nm的垂直分辨率。但是其逐点扫描时间长、效率低。测量单个靶丸需要超过2000个子孔径,测量时间超过两天,如此低的检测效率无法获得大规模运用,并且差动共聚焦敏感特性曲线易受表面倾角的影响。
在面探测法中以南京理工大学研制的表面缺陷检测装置为代表,该装置中具有短相干光源模块和显微干涉成像模块。针对靶丸的中空结构,通过控制腔长选择正确的干涉条纹同时实现光源内部实现参考光与测试光的分离以及移相,获取微球表面缺陷信息。但是,获取的内表面缺陷信息存在离焦并且受到外表面缺陷调制的影响,不能满足微球内表面缺陷高精度测量的要求。
发明内容
本发明的目的在于提供一种微球内表面缺陷测量方法,针对微球内表面缺陷测量存在外表面缺陷调制的问题,消除了外表面缺陷的影响,解决了成像在外表面时,内表面缺陷离焦的问题。
实现本发明目的的技术方案为:一种微球内表面缺陷测量方法,具体步骤为:
获取微球外表面测试光与参考光的干涉强度图,微球内表面测试光和参考光的干涉强度图以及参考光的光强分布图;
采用四部移相算法确定参考光与测试光之间的相位差、来自微球其他面的反射光强以及测试光光强,获得外表面干涉光场和内表面干涉光场的复振幅分布;
引入虚拟透镜进行复振幅重构获得重构后的外表面干涉光场和内表面干涉光场的复振幅分布;
将外表面测试光的复振幅转化成透射光重构复振幅;消除内表面测试光中光从介质到空气中外表面缺陷对内表面调制的影响,获得消去调制后的内表面重构复振幅;
采用自聚焦算法获取衍射距离;
根据衍射距离,采用逆向衍射计算获取衍射面复振幅分布,确定内表面缺陷相位信息。
优选地,采用四部移相算法确定参考光与测试光之间的相位差、来自微球其他面的反射光强以及测试光光强的具体步骤为:
分别将外表面干涉强度图和内表面干涉强度图分离出四张相位不同的干涉图,且每一张干涉图中的相位差在参考光与测试光相位差上分别增加了0、π/2、π、3π/2相位;
根据采集的参考光光强以及四张相位不同的干涉图确定参考光与测试光之间的相位差、来自微球其他面的反射光强以及测试光光强。
优选地,获得的外表面干涉光场和内表面干涉光场的复振幅分布分别为:
优选地,在偏振相机所处平面建立xy直角坐标系,并在此平面引入的虚拟透镜,虚拟透镜具体数学表达式为:
式中,j为虚数单位,k为波数,df表示参考光和测试光的汇聚点到偏振相机的距离。
优选地,重构后的外表面干涉光场和内表面干涉光场的复振幅分布分别为:
优选地,外表面测试光的复振幅转化成的透射光重构复振幅具体为:
消去调制后的内表面重构复振幅具体为:
式中,T为微球透射率,n微球折射率。
优选地,采用自聚焦算法获取衍射距离的具体步骤为:
根据透射光重构复振幅E″out(x,y)以及消去调制后的内表面重构复振幅E″inner(x,y),选取ROI区域;
采用光线追迹的方法预估对焦区间[L1,L2],对该区间以Δx为间隔生成一系列衍射距离d1,d2,...,dn;
提取ROI区域的一系列复振幅分布,计算ROI区域的聚焦评价函数在各个衍射距离下的值,选取最大值所对应的衍射距离为ROI区域的衍射距离。
优选地,所述聚焦评价函数具体为:
式中,N表示ROI区域中的像素总个数;Ids表示逆向衍射距离为d的时候,第s个元素的光强值。
优选地,逆向衍射计算具体公式为:
式中,E″(x,y)表示透射光重构复振幅或消去调制后的内表面重构复振幅,E″′(x,y)表示衍射面复振幅分布,H(fx,fy)为正向衍射过程的传递函数。
优选地,内表面缺陷相位具体为:
本发明与现有技术相比,其显著优点为:1)本发明引入了虚拟透镜,将曲面间的衍射计算转换成平面与平面间的衍射计算,避免了球面相位的引入所带来的采样问题;2)本发明改进了传统的自聚焦算法,消除衍射带来的能量涨落对聚焦距离判定的影响;3)本发明不会增加系统的复杂度,不需要额外引入新元件,从而避免了进一步的装调误差。
附图说明
图1为本发明微球内表面缺陷测量方法流程图。
图2为短相干光源模块示意图。
图3为显微干涉成像测量模块示意图。
图4为ROI区域内聚焦评价函数Md随距离d变化的曲线示意图。
图5为未消除外表面缺陷调制和内表面缺陷离焦的内表面缺陷测量结果示意图。
图6为消除外表面调制和内表面缺陷离焦的结果示意图。
具体实施方式
如图1所示,一种微球内表面缺陷测量方法,具体步骤如下:
步骤1:利用微球内表面缺陷测量装置获取微球外表面测试光与参考光的干涉强度图,微球内表面测试光和参考光的干涉强度图以及参考光的光强分布图。
所述微球内表面缺陷测量装置包括短相干光源模块和显微干涉成像测量模块两部分。如图2所示,短相干光源模块中短相干激光器出射光经过衰片和1/2波片,这里的1/2波片用于调节后续的S光和P光的振幅比,再经过PBS偏振分光分为S光和P光,经角锥棱镜反射通过PBS合束,最后经光纤耦合器耦合进光纤。在这里光源模块的主要作用用于产生一对正交具有一定延迟短相干线偏振光。如图3所示,耦合进光纤的光从光纤端口出射,经PBS偏振分光分为S光和P光。其中P光经显微物镜1、1/4波片照射到微球上,再经显微物镜1收集携带微球信息的外表面测试光和内表面测试光,再次经过1/4波片,最后经过PBS和1/4波片到达偏振相机。S光经显微物镜2,1/4波片到达参考镜,再次经过显微物镜和1/4波片,最后经过PBS和1/4波片到达偏振相机。在成像光路中,微球外表面顶点和偏振相机中心共轭,通过调节短相干光源模块中的平移台,当参考光与外表面测试光的光程差接近零时,此时参考光与外表面测试光干涉,通过偏振相机获取微球外表面测试光和参考光的干涉强度图。再次推动短相干光源模块中平移台,当参考光与内表面测试光光程差接近零时,此时参考光与内表面测试光干涉,通过偏振相机获取微球内表面测试光和参考光的干涉强度图。遮挡来自微球的反射光,通过偏振相机获取参考光的光强分布图。
步骤2:分别将外表面干涉强度图和内表面干涉强度图都分离出四张相位不同的干涉图,并且满足每一张干涉图中的相位差在参考光与测试光相位差基础上额外增加了0、π/2、π、3π/2的相位。假设上述四张干涉图的强度分布为I01、I02、I03、I04,有如下表达式:
根据上述采用四步移相算法获取的结果,在偏振相机所处平面建立xy直角坐标系,外表面干涉光场和内表面干涉光场的复振幅分布Eout(x,y)、Einner(x,y)表达式如下:
式中为参考光的相位因子;Iout、Iinner分别为外表面反射光强、内表面反射光强,其中Iinner=T2Iout,T为微球透射率;Ir为参考光强;为外表面反射光与参考光的相位差;为内表面反射光与参考光的相位差,该相位因子被外表面缺陷调制了两次。
步骤3:在步骤2建立的坐标系前提下,引入虚拟透镜进行复振幅重构,其中虚拟透镜表达式如下:
其中,j为虚数单位,k为波数,df表示参考光和测试光的汇聚点到偏振相机的距离。干涉场复振幅经过虚拟透镜之后,消除球面相位因子之后外表面干涉光场和内表面干涉光场的复振幅分布在步骤2建立的坐标系下的表达式为:
步骤4:在步骤2建立的坐标系前提下,根据微球在自准直情况下的透射系数因子和反射系数因子的关系,将外表面测试光的复振幅转化成透射光复振幅,其透射光重构复振幅表达式为:
根据透射系数和微球折射率的关系,消除内表面测试光中光从介质到空气中外表面缺陷对内表面调制的影响,消去调制后的内表面重构复振幅分布为:
其中,T为微球透射率,n微球折射率。
步骤5:采用改进的自聚焦算法获取衍射距离。
根据步骤4重构后的复振幅分布E″out(x,y)、E″inner(x,y),选取ROI(region ofinterest)区域。
采用光线追迹的方法预估一个对焦区间[L1,L2],对该区间以Δx为间隔生成一系列衍射距离d1,d2,...,dn。
提取ROI区域的一系列复振幅分布,将ROI区域的复振幅分布带入聚焦评价函数,聚焦评价函数表达式如下所示:
其中,N表示ROI区域中的像素总个数;Ids表示逆向衍射距离为d的时候,第s个元素的光强值。
计算ROI区域的聚焦评价函数在各个衍射距离下的值,选取最大值所对应的衍射距离即为ROI区域的衍射距离,结果如图4所示。
步骤6:采用逆向衍射计算消除内表面缺陷离焦影响和光从空气进入微球外表面缺陷调制的影响,其逆向衍射计算表达式如下:
其中E″(x,y)表示步骤4重构的复振幅分布,E″′(x,y)表示衍射面复振幅分布,H(fx,fy)为正向衍射过程的传递函数,其表达式是关于衍射距离的函数:
其中,j为虚数单位,d为衍射距离,k为波数,fx为x方向的频谱坐标,fy为y方向的频谱坐标。
对步骤4重构的复振幅分布E″out(x,y)、E″inner(x,y),结合逆向衍射计算以及步骤5获取的衍射距离,获取衍射面复振幅分布E″′out(x,y)、E″′inner(x,y)。结合衍射面的复振幅分布消除光从空气进入微球外表面缺陷调制的影响,获得内表面缺陷相位信息,具体表达式如下:
Claims (10)
1.一种微球内表面缺陷测量方法,其特征在于,具体步骤为:
获取微球外表面测试光与参考光的干涉强度图,微球内表面测试光和参考光的干涉强度图以及参考光的光强分布图;
采用四部移相算法确定参考光与测试光之间的相位差、来自微球其他面的反射光强以及测试光光强,获得外表面干涉光场和内表面干涉光场的复振幅分布;
引入虚拟透镜进行复振幅重构获得重构后的外表面干涉光场和内表面干涉光场的复振幅分布;
将外表面测试光的复振幅转化成透射光重构复振幅;消除内表面测试光中光从介质到空气中外表面缺陷对内表面调制的影响,获得消去调制后的内表面重构复振幅;
采用自聚焦算法获取衍射距离;
根据衍射距离,采用逆向衍射计算获取衍射面复振幅分布,确定内表面缺陷相位信息。
2.根据权利要求1所述的微球内表面缺陷测量方法,其特征在于,采用四部移相算法确定参考光与测试光之间的相位差、来自微球其他面的反射光强以及测试光光强的具体步骤为:
分别将外表面干涉强度图和内表面干涉强度图分离出四张相位不同的干涉图,且每一张干涉图中的相位差在参考光与测试光相位差上分别增加了0、π/2、π、3π/2相位;
根据采集的参考光光强以及四张相位不同的干涉图确定参考光与测试光之间的相位差、来自微球其他面的反射光强以及测试光光强。
7.根据权利要求1所述的微球内表面缺陷测量方法,其特征在于,采用自聚焦算法获取衍射距离的具体步骤为:
根据透射光重构复振幅E″out(x,y)以及消去调制后的内表面重构复振幅E″inner(x,y),选取ROI区域;
采用光线追迹的方法预估对焦区间[L1,L2],对该区间以Δx为间隔生成一系列衍射距离d1,d2,...,dn;
提取ROI区域的一系列复振幅分布,计算ROI区域的聚焦评价函数在各个衍射距离下的值,选取最大值所对应的衍射距离为ROI区域的衍射距离。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010335374.4A CN111650203B (zh) | 2020-04-24 | 2020-04-24 | 一种微球内表面缺陷测量方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010335374.4A CN111650203B (zh) | 2020-04-24 | 2020-04-24 | 一种微球内表面缺陷测量方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111650203A true CN111650203A (zh) | 2020-09-11 |
CN111650203B CN111650203B (zh) | 2023-03-28 |
Family
ID=72346504
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010335374.4A Active CN111650203B (zh) | 2020-04-24 | 2020-04-24 | 一种微球内表面缺陷测量方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111650203B (zh) |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112345551A (zh) * | 2020-11-05 | 2021-02-09 | 菲特(天津)检测技术有限公司 | 一种飞机导管喇叭口内外壁表面缺陷检测方法及系统 |
CN112630232A (zh) * | 2020-12-01 | 2021-04-09 | 北京理工大学 | 差动共焦定面干涉靶丸内、外表面缺陷检测方法与装置 |
CN112666172A (zh) * | 2020-12-01 | 2021-04-16 | 北京理工大学 | 差动共焦定面干涉靶丸外表面缺陷检测方法与装置 |
CN112683918A (zh) * | 2020-12-01 | 2021-04-20 | 北京理工大学 | 差动共焦定面干涉靶丸内表面缺陷检测方法与装置 |
CN113533366A (zh) * | 2021-07-06 | 2021-10-22 | 常州联影智融医疗科技有限公司 | 导光束端面的缺陷检测装置及其检测方法 |
CN115598147A (zh) * | 2022-09-30 | 2023-01-13 | 南京理工大学(Cn) | 基于白光显微干涉的微球内外表面缺陷检测装置及方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2005062113A (ja) * | 2003-08-20 | 2005-03-10 | Olympus Corp | 屈折率分布の測定方法及び測定装置 |
CN103196361A (zh) * | 2013-02-28 | 2013-07-10 | 哈尔滨工业大学 | 用于微球表面形貌快速检测的短相干瞬时移相干涉测量仪及测量方法 |
CN109990709A (zh) * | 2019-04-19 | 2019-07-09 | 北京理工大学 | 双边错位差动共焦干涉靶丸形貌轮廓参数测量方法与装置 |
CN110030941A (zh) * | 2019-03-08 | 2019-07-19 | 北京理工大学 | 激光共焦干涉核聚变靶丸形貌轮廓参数测量方法与装置 |
CN110030942A (zh) * | 2019-03-08 | 2019-07-19 | 北京理工大学 | 激光差动共焦干涉核聚变靶丸形貌参数测量方法与装置 |
-
2020
- 2020-04-24 CN CN202010335374.4A patent/CN111650203B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2005062113A (ja) * | 2003-08-20 | 2005-03-10 | Olympus Corp | 屈折率分布の測定方法及び測定装置 |
CN103196361A (zh) * | 2013-02-28 | 2013-07-10 | 哈尔滨工业大学 | 用于微球表面形貌快速检测的短相干瞬时移相干涉测量仪及测量方法 |
CN110030941A (zh) * | 2019-03-08 | 2019-07-19 | 北京理工大学 | 激光共焦干涉核聚变靶丸形貌轮廓参数测量方法与装置 |
CN110030942A (zh) * | 2019-03-08 | 2019-07-19 | 北京理工大学 | 激光差动共焦干涉核聚变靶丸形貌参数测量方法与装置 |
CN109990709A (zh) * | 2019-04-19 | 2019-07-09 | 北京理工大学 | 双边错位差动共焦干涉靶丸形貌轮廓参数测量方法与装置 |
Cited By (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112345551A (zh) * | 2020-11-05 | 2021-02-09 | 菲特(天津)检测技术有限公司 | 一种飞机导管喇叭口内外壁表面缺陷检测方法及系统 |
CN112630232A (zh) * | 2020-12-01 | 2021-04-09 | 北京理工大学 | 差动共焦定面干涉靶丸内、外表面缺陷检测方法与装置 |
CN112666172A (zh) * | 2020-12-01 | 2021-04-16 | 北京理工大学 | 差动共焦定面干涉靶丸外表面缺陷检测方法与装置 |
CN112683918A (zh) * | 2020-12-01 | 2021-04-20 | 北京理工大学 | 差动共焦定面干涉靶丸内表面缺陷检测方法与装置 |
CN112630232B (zh) * | 2020-12-01 | 2021-12-03 | 北京理工大学 | 差动共焦定面干涉靶丸内、外表面缺陷检测方法与装置 |
CN112666172B (zh) * | 2020-12-01 | 2022-02-11 | 北京理工大学 | 差动共焦定面干涉靶丸外表面缺陷检测方法与装置 |
CN112683918B (zh) * | 2020-12-01 | 2022-02-11 | 北京理工大学 | 差动共焦定面干涉靶丸内表面缺陷检测方法与装置 |
CN113533366A (zh) * | 2021-07-06 | 2021-10-22 | 常州联影智融医疗科技有限公司 | 导光束端面的缺陷检测装置及其检测方法 |
CN113533366B (zh) * | 2021-07-06 | 2024-05-24 | 常州联影智融医疗科技有限公司 | 导光束端面的缺陷检测装置及其检测方法 |
CN115598147A (zh) * | 2022-09-30 | 2023-01-13 | 南京理工大学(Cn) | 基于白光显微干涉的微球内外表面缺陷检测装置及方法 |
CN115598147B (zh) * | 2022-09-30 | 2024-05-17 | 南京理工大学 | 基于白光显微干涉的微球内外表面缺陷检测装置及方法 |
Also Published As
Publication number | Publication date |
---|---|
CN111650203B (zh) | 2023-03-28 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111650203B (zh) | 一种微球内表面缺陷测量方法 | |
CN102494623B (zh) | 镜头中光学表面中心间距的非接触式测量装置的测量方法 | |
CN109211934B (zh) | 基于干涉显微的微球面缺陷检测装置及其检测方法 | |
CN109556531B (zh) | 一种基于图像信息的点衍射干涉仪光路精确校准系统及方法 | |
CN111121675B (zh) | 一种用于微球表面显微干涉测量的视场扩展方法 | |
CN111257227B (zh) | 基于偏振自相关的暗场共焦显微测量装置和方法 | |
US11635289B2 (en) | Surface shape measurement device and surface shape measurement method | |
CN108562241B (zh) | 基于光纤束的数字全息柔性测量的装置与方法 | |
CN115930773A (zh) | 轻离轴数字全息检测装置 | |
CN113125387A (zh) | 一种增强表面等离子体共振相位成像横向分辨率的方法 | |
CN111289479B (zh) | 基于非线性热像反演的相位缺陷检测装置和方法 | |
CN115598147B (zh) | 基于白光显微干涉的微球内外表面缺陷检测装置及方法 | |
CN112630232B (zh) | 差动共焦定面干涉靶丸内、外表面缺陷检测方法与装置 | |
CN110631510A (zh) | 一种基于迈克尔逊结构的高精度测角装置及测角方法 | |
CN112666172B (zh) | 差动共焦定面干涉靶丸外表面缺陷检测方法与装置 | |
CN112683918B (zh) | 差动共焦定面干涉靶丸内表面缺陷检测方法与装置 | |
CN114371549B (zh) | 一种基于多焦复用透镜的定量相位成像方法及系统 | |
TWI632361B (zh) | 數位全像顯微斷層之方法及裝置 | |
CN113418470B (zh) | 光谱扫描共焦单次曝光数字全息测量系统及测量方法 | |
CN113009801B (zh) | 高速多方向线共焦数字全息三维显微成像方法及装置 | |
CN104534980A (zh) | 一种反射型无透镜数字全息测量装置 | |
CN205538737U (zh) | 反射型双波长全息术的光学元件高景深表面疵病检测装置 | |
CN112711030A (zh) | 一种显微镜面阵扫频测量装置和方法 | |
He et al. | On-line beam diagnostics based on single-shot beam splitting phase retrieval | |
JP3871183B2 (ja) | 光学素子の3次元形状測定方法及び測定装置 |
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 |