CN102661956A - 超光滑表面缺陷检测系统及其畸变校正方法 - Google Patents

超光滑表面缺陷检测系统及其畸变校正方法 Download PDF

Info

Publication number
CN102661956A
CN102661956A CN2012101229019A CN201210122901A CN102661956A CN 102661956 A CN102661956 A CN 102661956A CN 2012101229019 A CN2012101229019 A CN 2012101229019A CN 201210122901 A CN201210122901 A CN 201210122901A CN 102661956 A CN102661956 A CN 102661956A
Authority
CN
China
Prior art keywords
gauge plate
image
distortion
subimage
pixel
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
CN2012101229019A
Other languages
English (en)
Other versions
CN102661956B (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 University ZJU
Original Assignee
Zhejiang University ZJU
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 University ZJU filed Critical Zhejiang University ZJU
Priority to CN 201210122901 priority Critical patent/CN102661956B/zh
Publication of CN102661956A publication Critical patent/CN102661956A/zh
Application granted granted Critical
Publication of CN102661956B publication Critical patent/CN102661956B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Image Processing (AREA)

Abstract

本发明提出了一种超光滑表面缺陷检测系统及其畸变校正方法。本发明解决了超光滑表面缺陷检测系统中因为存在光学畸变而造成的子图像拼接时的缺陷断裂问题。本发明的技术特点在于,设计了一种超光滑表面缺陷检测系统和畸变校正标准板及标准板的夹持装置,利用检测系统对标准板采集得到暗场畸变图像,通过该畸变图像和计算机按照物面尺寸和像面像素关系重构的标准板理想图像的匹配建立畸变退化模型,并提出一种基于二次极坐标正反变换和二次灰度线性插值的畸变校正方法。本方法能够校正超光滑表面缺陷检测系统中因畸变造成的相邻子图像的拼接错位,同时也适用于校正其他基于图像拼接的大视场光学系统中存在的畸变。

Description

超光滑表面缺陷检测系统及其畸变校正方法
技术领域
本发明涉及一种超光滑表面缺陷检测系统及其畸变校正方法。
背景技术
大尺寸精密光学元件表面缺陷的定量检测一直是国际上公认的光学检测难题之一。被检表面是一个宏观量级,而表面上待分辨的缺陷有时需要微观到微米量级。发明人研制的“超光滑表面缺陷检测系统”利用光学显微散射暗场成像、子图像扫描及基于特征匹配的拼接等技术能够实现对精密光学元件表面缺陷的自动化检测和数字化评价。超光滑表面缺陷检测系统可以完成最大尺寸为430mm×430mm的光学元件全口径表面缺陷自动化定量检测,横向分辨率可以达到0.5μm。超光滑表面缺陷检测系统中存在枕形畸变。如果不在子图像进行图像处理之前修正这种畸变,那么相邻子孔径图像进行拼接时,可能会造成位于重叠区域的缺陷(尤其是划痕)出现断裂的情况;枕形畸变的存在也会使对缺陷的特征识别过程的结果(如长度、曲率)出现错误。在保证超光滑表面缺陷检测系统所要求的参数指标的前提下,改变光学设计成本较高而且难于实现。因此,本发明提出的畸变校正方法对于超光滑表面缺陷检测系统工作时子图像快速的拼接和缺陷的准确检测就显得极为关键。
发明内容
本发明的目的是克服现有技术的不足,提供一种超光滑表面缺陷检测系统及其畸变校正方法。
超光滑表面缺陷检测系统包括二维移导装置、计算机、光学显微成像装置、标准板和标准板的夹持装置,光学显微成像装置包括LED环形照明光源、变倍显微镜、CCD探测器;CCD探测器、变倍显微镜、LED环形照明光源顺次连接组成检测系统的光学显微成像装置,CCD探测器通过通信线缆与计算机相连,光学显微成像装置安放在二维移导装置上,二维移导装置带动光学显微成像装置一起运动;LED环形照明光源产生奇数束光束经过准直变成平行光,以角度α为27°~33°入射到待检元件或标准板表面,其中一束入射光1,照射到表面缺陷,变成散射光3进入显微镜,在CCD探测器形成待检表面的暗场子图像,传输并保存到计算机中,二维移导装置带动光学显微成像装置相对于待检元件,从待检元件左上边缘开始,进行路线为S形的扫描运动,并以子图像行列坐标位置A11,A21,A31,...,AM1,AM2,...,AMN命名图像矩阵,完成对整个待检元件的全通光口径图像采集,相邻子图像之间具有1/4~1/6的重叠区域,然后计算机拼接子图像,得到大尺寸元件全口径的疵病图像。
所述的标准板采用电子束曝光、反应离子束刻蚀的方法在石英玻璃板上制作图案,标准板分为若干区域,每个区域分别刻有相同的网格阵列,标准板从侧面插入标准板的夹持装置,通过螺钉固紧;标准板的夹持装置设有通光孔,刻蚀的网格线与表面缺陷等效,网格线诱发入射光散射后进入显微镜,形成标准板的暗场图像。
畸变校正方法包括如下步骤:
1)将标准板和标准板的夹持装置放置在系统显微镜的工作距离处,调整标准板所在平面与显微镜物面重合,显微镜高倍下,使用网格较密的区域;低倍下,使用网格较稀疏的区域,使获得的标准板畸变图像有足够多的网格交点,标准板经过显微镜在CCD探测器形成暗场畸变图像,通过标准板尺寸大小、显微镜的放大倍率以及CCD像素大小之间的关系重构标准板不存在畸变时的理想网格图像;
2)在标准板理想网格图像上选取m个网格交点,m个网格交点应等间距分布在经过图像中心的一半对角线上,同时选取m个网格交点在标准板畸变图像上的m个对应点,获得以图像左上角为原点的像素直角坐标值,将坐标值转换为以图像中心为原点的极坐标值;
3)建立一个基于极坐标变换的畸变退化多项式模型,
ρ ′ = T ρ ( ρ ) = Σ i = 0 n a i ρ i = a 0 + a 1 ρ + a 2 ρ 2 + . . . + a n ρ n - - - ( 1 )
θ′=Tθ(θ)=θ
式(1)的物理意义在于标准板理想图像上一点的极坐标(ρ,θ)经过畸变退化变换Tρ,Tθ,变成标准板畸变图像上的对应点(ρ′,θ′),将式(1)中各次项系数写成向量形式
Figure BDA0000156818780000022
a → t = [ a 0 , a 1 , a 2 , . . . , a n ]
将标准板理想图像上取的m个像素的极坐标第一项ρj,j=1,2,...,m,及其各次项写成矩阵A,
A = 1 ρ 1 ρ 1 2 . . . ρ 1 n 1 ρ 2 ρ 2 2 . . . ρ 2 n . . . . . . . . . . 1 ρ m ρ m 2 . . . ρ m n
将理想图像上取的m个点极坐标第一项ρ′j,j=1,2,...,m,写成向量形式
Figure BDA0000156818780000032
ρ ′ → t = [ ρ 1 ′ , ρ 2 ′ , . . . , ρ m ′ ]
则,畸变退化多项式(1)中各次项系数
Figure BDA0000156818780000034
为:
a → t = [ A - 1 ρ ′ → ] t = [ a 0 , a 1 , . . . , a n ] - - - ( 2 )
4)对系统检测时的每一幅子图像,将子图像不存在畸变时的理想子图像上每个像素P的以图像左上角为原点的直角坐标值(x,y)转换为以图像中心为原点的极坐标值(ρ,θ),代入式(1),获得其对应的畸变子图像像素坐标值(ρ′,θ′),再将此极坐标值转换为以图像左上角为原点的直角坐标值(u,v),若u,v为整数,则选取畸变子图像上对应像素P′(u,v)的灰度值作为理想子图像像素P的灰度值G(P);若u,v为非整数,则P为亚像素,采用畸变子图像上与(u,v)最为邻近的四个像素A(i,j),B(i+1,j),C(i,j+1),D(i+1,j+1),用四个像素的灰度值进行双线性插值后的灰度值作为理想子图像上像素P′的灰度值G(P′);
设α=u-i,β=v-j
E为线段
Figure BDA0000156818780000036
上一点,F为线段
Figure BDA0000156818780000037
上一点,E,F也为亚像素,
Figure BDA0000156818780000038
平行于
Figure BDA0000156818780000039
Figure BDA00001568187800000310
首先用第一次线性插值计算出E,F两处的灰度值G(E)和G(F)
G(E)=β[G(C)-G(A)]+G(A)            (3)
G(F)=β[G(D)-G(B)]+G(B)
式(3)中G(A),G(B),G(C),G(D)分别代表A,B,C,D四个像素的灰度值,同理,再用第二次线性插值计算出畸变子图像上像素P′的灰度值G(P′),即理想子图像上对应像素P灰度值G(P);
G(P)=G(P′)=α[G(F)-G(E)]+G(E)    (4)
本发明提出的超光滑表面缺陷检测系统中,待检元件采用立式放置方式,可以减小元件自重产生的应力对于待检元件外形改变的影响;本发明提出的畸变校正方法,能够很好地校正子图像拼接时缺陷的断裂和错位现象,为超光滑表面缺陷检测系统后续正确的特征识别和标定提供了基础。同时,这种方法同样可应用于校正其他基于图像拼接的大视场光学检测系统中的畸变。
附图说明
图1是超光滑表面缺陷检测系统的原理及结构示意图;
图2是超光滑表面缺陷检测系统工作时子图像扫描路径;
图3(a)是畸变校正中使用的二元光学网格标准板1设计图;
图3(b)是畸变校正中使用的二元光学网格标准板2设计图;
图4(a)是畸变校正网格标准板夹持装置正视图;
图4(b)是畸变校正网格标准板夹持装置侧视图;
图5是标准板的理想图像上选取像素的位置分布图
图6是灰度二次线性插值原理图;
图7(a)是检测系统对网格标准板采样获得的暗场畸变图像;
图7(b)是通过计算机利用物像关系重构的网格标准板的理想图像;
图7(c)是网格标准板畸变图像通过本发明畸变校正方法得到的校正后图像;
图8是超光滑表面缺陷检测系统及其畸变校正方法流程图。
具体实施方式
如图1所示,超光滑表面缺陷检测系统包括二维移导装置S5、计算机S6、光学显微成像装置S7、标准板S8和标准板的夹持装置S9,光学显微成像装置S7包括LED环形照明光源S2、变倍显微镜S3、CCD探测器S4;CCD探测器S4、变倍显微镜S3、LED环形照明光源S2顺次连接组成检测系统的光学显微成像装置S7,CCD探测器S4通过通信线缆与计算机S6相连,光学显微成像装置S7安放在二维移导装置S5上,二维移导装置S5带动光学显微成像装置一起运动;LED环形照明光源S2产生奇数束光束经过准直变成平行光,以角度α为27°~33°入射到待检元件S1或标准板S8表面,其中一束入射光1,照射到表面缺陷,变成散射光3进入显微镜S3,在CCD探测器S4形成待检表面的暗场子图像,传输并保存到计算机S6中;若表面光滑无缺陷,入射光1变成反射光2,不会进入显微镜S3;二维移导装置S5带动光学显微成像装置S7相对于待检元件S1,从待检元件S1左上边缘开始,进行路线为S形的扫描运动,如图2所示,并以子图像行列坐标位置A11,A21,A31,...,AM1,AM2,...,AMN命名图像矩阵,完成对整个待检元件S1的全通光口径图像采集,相邻子图像之间具有1/4~1/6的重叠区域,然后计算机S6拼接子图像,得到大尺寸元件全口径的疵病图像。
图3是所述的标准板的设计图,标准板S8采用电子束曝光、反应离子束刻蚀的方法在石英玻璃板上制作图案,标准板分为若干区域,每个区域分别刻有相同的网格阵列,不同区域的网格边长各不相同;图4是所述的标准板的夹持装置,使用时标准板S8从侧面插入标准板的夹持装置S9,通过螺钉固紧;标准板的夹持装置S9设有通光孔,在无刻蚀的光滑区域,入射光反射或透射偏出,不会进入显微镜;而刻蚀的网格线与表面缺陷等效,网格线诱发入射光散射后进入显微镜S3,形成标准板的暗场图像。
图8为使用上述装置的基于二次极坐标正反变换畸变校正方法流程图,畸变校正方法包括如下步骤:
1)将标准板S8和标准板的夹持装置S9放置在系统显微镜S3的工作距离处,调整标准板所在平面与显微镜物面重合,显微镜高倍下,使用网格较密的区域;低倍下,使用网格较稀疏的区域,使获得的标准板畸变图像有足够多的网格交点,标准板S8经过显微镜S3在CCD探测器S4形成暗场畸变图像,如图7(a)所示;通过标准板S8尺寸大小、显微镜S3的放大倍率以及CCD探测器S4像素大小之间的关系重构标准板不存在畸变时的理想网格图像,如图7(b)所示;
2)在标准板理想网格图像上选取m个网格交点,m个网格交点应等间距分布在经过图像中心的一半对角线上,如图5所示,同时选取m个网格交点在标准板畸变图像上的m个对应点,获得以图像左上角为原点的像素直角坐标值,将坐标值转换为以图像中心为原点的极坐标值;
3)建立一个基于极坐标变换的畸变退化多项式模型,
ρ ′ = T ρ ( ρ ) = Σ i = 0 n a i ρ i = a 0 + a 1 ρ + a 2 ρ 2 + . . . + a n ρ n - - - ( 1 )
θ′=Tθ(θ)=θ
式(1)的物理意义在于标准板理想图像上一点的极坐标(ρ,θ)经过畸变退化变换Tρ,Tθ,变成标准板畸变图像上的对应点(ρ′,θ′),将式(1)中各次项系数写成向量形式
a → t = [ a 0 , a 1 , a 2 , . . . , a n ]
将标准板理想图像上取的m个像素的极坐标第一项ρj,j=1,2,...,m,及其各次项写成矩阵形式,
A = 1 ρ 1 ρ 1 2 . . . ρ 1 n 1 ρ 2 ρ 2 2 . . . ρ 2 n . . . . . . . . . . 1 ρ m ρ m 2 . . . ρ m n
将理想图像上取的m个点极坐标第一项ρ′j,j=1,2,...,m,写成向量形式
Figure BDA0000156818780000062
ρ ′ → t = [ ρ 1 ′ , ρ 2 ′ , . . . , ρ m ′ ]
A,
Figure BDA0000156818780000065
代入式(1),得到,
A a → = ρ ′ →
左右两边求逆,则畸变退化多项式(1)中各次项系数
Figure BDA0000156818780000067
为:
a → t = [ A - 1 ρ ′ → ] t = [ a 0 , a 1 , . . . , a n ] - - - ( 2 )
4)对系统检测时的每一幅子图像,将子图像不存在畸变时的理想子图像上每个像素P的以图像左上角为原点的直角坐标值(x,y)转换为以图像中心为原点的极坐标值(ρ,θ),代入式(1),获得其对应的畸变子图像像素坐标值(ρ′,θ′),再将此极坐标值转换为以图像左上角为原点的直角坐标值(u,v),若u,v为整数,则选取畸变子图像上对应像素P′(u,v)的灰度值作为理想子图像像素P的灰度值G(P);若u,v为非整数,则P为亚像素,采用畸变子图像上与(u,v)最为邻近的四个像素A(i,j),B(i+1,j),C(i,j+1),D(i+1,j+1),如图6所示,用四个像素的灰度值进行双线性插值后的灰度值作为理想子图像上像素P′的灰度值G(P′);
设α=u-i,β=v-j
E为线段
Figure BDA0000156818780000069
上一点,F为线段
Figure BDA00001568187800000610
上一点,E,F也为亚像素,
Figure BDA00001568187800000611
平行于
Figure BDA00001568187800000612
Figure BDA00001568187800000613
首先用第一次线性插值计算出E,F两处的灰度值G(E)和G(F)
G(E)=β[G(C)-G(A)]+G(A)          (3)
G(F)=β[G(D)-G(B)]+G(B)
式(3)中G(A),G(B),G(C),G(D)分别代表A,B,C,D四个像素的灰度值,同理,再用第二次线性插值计算出畸变子图像上像素P′的灰度值G(P′),即理想子图像上对应像素P灰度值G(P);
G(P)=G(P′)=α[G(F)-G(E)]+G(E)  (4)
在逐一对每幅子图像进行畸变校正后,再将校正后的子图像按照重叠区域拼接,组成待检元件的全通光口径的检测图片。
实施例
本发明应用一种超光滑表面缺陷检测系统及其畸变校正方法实例如下。
如图1所示,超光滑表面缺陷检测系统包括二维移导装置S5、计算机S6、光学显微成像装置S7、标准板S8和标准板的夹持装置S9,光学显微成像装置S7包括LED环形照明光源S2、变倍显微镜S3、CCD探测器S4;CCD探测器S4、变倍显微镜S3、LED环形照明光源S2顺次连接组成检测系统的光学显微成像装置S7,CCD探测器S4通过通信线缆与计算机S6相连,光学显微成像装置S7安放在二维移导装置S5上,二维移导装置S5带动光学显微成像装置一起运动;LED环形照明光源S2产生奇数束光束经过准直变成平行光,以角度α为27°~33°入射到待检元件S1或标准板S8表面,其中一束入射光1,照射到表面缺陷,变成散射光3进入显微镜S3,在CCD探测器S4形成待检表面的暗场子图像,传输并保存到计算机S6中;若表面光滑无缺陷,入射光1变成反射光2,不会进入显微镜S3;二维移导装置S5带动光学显微成像装置S7相对于待检元件S1,从待检元件S1左上边缘开始,进行路线为S形的扫描运动,如图2所示,并以子图像行列坐标位置A11,A21,A31,...,AM1,AM2,...,AMN命名图像矩阵,完成对整个待检元件S1的全通光口径图像采集,相邻子图像之间具有1/4~1/6的重叠区域,然后计算机S6拼接子图像,得到大尺寸元件全口径的疵病图像。
图3是所述的标准板的设计图,标准板S8采用电子束曝光、反应离子束刻蚀的方法在石英玻璃板上制作图案,标准板分为若干区域,每个区域分别刻有相同的网格阵列,不同区域的网格边长各不相同,显微镜高倍率下使用网格边长较短的区域,显微镜低倍下使用网格边长较长的区域。本例中,显微镜放大倍率为1x时,选取标准板上网格阵列宽度为30mm×30mm,每个网格的边长0.3mm。
图8为使用上述装置的基于二次极坐标正反变换畸变校正方法流程图,畸变校正方法包括如下步骤:
1)图4是所述的标准板的夹持装置,使用时标准板S8从侧面插入标准板的夹持装置S9,通过螺钉固紧。将标准板S8和标准板的夹持装置S9放置在系统显微镜S3的工作距离处,调整标准板所在平面与显微镜物面重合。标准板S8经过显微镜S3在CCD探测器S4形成暗场畸变图像,如图7(a)所示;通过标准板S8尺寸大小、显微镜S3的放大倍率以及CCD探测器S4像素大小之间的关系重构标准板不存在畸变时的理想网格图像,如图7(b)所示;
2)在标准板理想网格图像上选取m个网格交点,m个网格交点应等间距分布在经过图像中心的一半对角线上,如图5所示,同时选取m个网格交点在标准板畸变图像上的m个对应点,获得以图像左上角为原点的像素直角坐标值,将坐标值转换为以图像中心为原点的极坐标值;
3)建立一个基于极坐标变换的畸变退化多项式模型,
ρ ′ = T ρ ( ρ ) = Σ i = 0 n a i ρ i = a 0 + a 1 ρ + a 2 ρ 2 + . . . + a n ρ n - - - ( 1 )
θ′=Tθ(θ)=θ
式(1)的物理意义在于标准板理想图像上一点的极坐标(ρ,θ)经过畸变退化变换Tρ,Tθ,变成标准板畸变图像上的对应点(ρ′,θ′),
对畸变模型式进行简化,假设只有ρ′被展开到四次项,即n=4,而畸变前后θ没有变化,
ρ′=a0+a1ρ+a2ρ2+a3ρ3+a4ρ4        (1’)
θ′=θ
将式(1)中各次项系数写成向量形式
Figure BDA0000156818780000082
a → t = [ a 0 , a 1 , a 2 , a 3 , a 4 ]
将标准板理想图像上取的m个像素的极坐标第一项ρj,j=1,2,...,m,及其各次项写成矩阵形式,
A = 1 ρ 1 ρ 1 2 . . . ρ 1 n 1 ρ 2 ρ 2 2 . . . ρ 2 n . . . . . . . . . . 1 ρ m ρ m 2 . . . ρ m n
将理想图像上取的m个点极坐标第一项ρ′j,j=1,2,...,m,写成向量形式
Figure BDA0000156818780000092
ρ ′ → t = [ ρ 1 ′ , ρ 2 ′ , . . . , ρ m ′ ]
A,代入式(1),得到,
A a → = ρ ′ →
左右两边求逆,则畸变退化多项式(1)中各次项系数
Figure BDA0000156818780000097
为:
a → t = [ A - 1 ρ ′ → ] t = [ a 0 , a 1 , a 2 , , a 3 , a 4 ] - - - ( 2 )
各系数值见表1.
表1.畸变退化多项式模型中各次项系数
  系数   a0   a1   a2   a3   a4
  值   0.135077806   0.9962497   2.01E-05   -3.32E-08   2.43E-11
标准板校正后的结果如图7(c)所示。
4)对系统检测时的每一幅子图像,将子图像不存在畸变时的理想子图像上每个像素P的以图像左上角为原点的直角坐标值(x,y)转换为以图像中心为原点的极坐标值(ρ,θ),代入式(1’),获得其对应的畸变子图像像素坐标值(ρ′,θ′),再将此极坐标值转换为以图像左上角为原点的直角坐标值(u,v),若u,v为整数,则选取畸变子图像上对应像素P′(u,v)的灰度值作为理想子图像像素P的灰度值G(P);若u,v为非整数,则P为亚像素,采用畸变子图像上与(u,v)最为邻近的四个像素A(i,j),B(i+1,j),C(i,j+1),D(i+1,j+1),如图6所示,用四个像素的灰度值进行双线性插值后的灰度值作为理想子图像上像素P′的灰度值G(P′);
设α=u-i,β=v-j
E为线段
Figure BDA0000156818780000101
上一点,F为线段
Figure BDA0000156818780000102
上一点,E,F也为亚像素,
Figure BDA0000156818780000103
平行于
Figure BDA0000156818780000104
Figure BDA0000156818780000105
首先用第一次线性插值计算出E,F两处的灰度值G(E)和G(F)
G(E)=β[G(C)-G(A)]+G(A)         (3)
G(F)=β[G(D)-G(B)]+G(B)
式(3)中G(A),G(B),G(C),G(D)分别代表A,B,C,D四个像素的灰度值,同理,再用第二次线性插值计算出畸变子图像上像素P′的灰度值G(P′),即理想子图像上对应像素P灰度值G(P);
G(P)=G(P′)=α[G(F)-G(E)]+G(E) (4)
在逐一对每幅子图像进行畸变校正后,再将校正后的子图像按照重叠区域拼接,组成待检元件的全通光口径的检测图片。

Claims (3)

1.一种超光滑表面缺陷检测系统,其特征在于包括二维移导装置(S5)、计算机(S6)、光学显微成像装置(S7)、标准板(S8)和标准板的夹持装置(S9),光学显微成像装置(S7)包括LED环形照明光源(S2)、变倍显微镜(S3)、CCD探测器(S4);CCD探测器(S4)、变倍显微镜(S3)、LED环形照明光源(S2)顺次连接组成检测系统的光学显微成像装置(S7),CCD探测器(S4)通过通信线缆与计算机(S6)相连,光学显微成像装置(S7)安放在二维移导装置(S5)上,二维移导装置(S5)带动光学显微成像装置一起运动;LED环形照明光源(S2)产生奇数束光束经过准直变成平行光,以角度α为27°~33°入射到待检元件(S1)或标准板(S8)表面,其中一束入射光1,照射到表面缺陷,变成散射光3进入显微镜(S3),在CCD探测器(S4)形成待检表面的暗场子图像,传输并保存到计算机(S6)中,二维移导装置(S5)带动光学显微成像装置(S7)相对于待检元件(S1),从待检元件(S1)左上边缘开始,进行路线为S形的扫描运动,并以子图像行列坐标位置A11,A21,A31,...,AM1,AM2,...,AMN命名图像矩阵,完成对整个待检元件(S1)的全通光口径图像采集,相邻子图像之间具有1/4~1/6的重叠区域,然后计算机(S6)拼接子图像,得到大尺寸元件全口径的疵病图像。
2.根据权利要求1所述的一种超光滑表面缺陷检测系统,其特征在于所述的标准板(S8)采用电子束曝光、反应离子束刻蚀的方法在石英玻璃板上制作图案,标准板分为若干区域,每个区域分别刻有相同的网格阵列,标准板(S8)从侧面插入标准板的夹持装置(S9),通过螺钉固紧;标准板的夹持装置(S9)设有通光孔,刻蚀的网格线与表面缺陷等效,网格线诱发入射光散射后进入显微镜(S3),形成标准板的暗场图像。
3.一种使用权利要求1所述系统的畸变校正方法,其特征在于包括如下步骤:
1)将标准板(S8)和标准板的夹持装置(S9)放置在系统显微镜(S3)的工作距离处,调整标准板所在平面与显微镜物面重合,显微镜高倍下,使用网格较密的区域;低倍下,使用网格较稀疏的区域,使获得的标准板畸变图像有足够多的网格交点,标准板(S8)经过显微镜(S3)在CCD探测器(S4)形成暗场畸变图像,通过标准板(S8)尺寸大小、显微镜(S3)的放大倍率以及CCD(S4)像素大小之间的关系重构标准板不存在畸变时的理想网格图像;
2)在标准板理想网格图像上选取m个网格交点,m个网格交点应等间距分布在经过图像中心的一半对角线上,同时选取m个网格交点在标准板畸变图像 上的m个对应点,获得以图像左上角为原点的像素直角坐标值,将坐标值转换为以图像中心为原点的极坐标值;
3)建立一个基于极坐标变换的畸变退化多项式模型,
Figure FDA0000156818770000021
θ′=Tθ(θ)=θ
式(1)的物理意义在于标准板理想图像上一点的极坐标(ρ,θ)经过畸变退化变换Tρ,Tθ,变成标准板畸变图像上的对应点(ρ′,θ′),n为畸变退化多项式中ρ的最高次项的次数,将式(1)中各次项系数写成向量形式 
Figure FDA0000156818770000022
Figure FDA0000156818770000023
在标准板的理想图像上取的m个像素,m>n,将m个像素的极坐标第一项ρj,j=1,2,...,m,及其各次项写成矩阵形式,矩阵A,
将理想图像上取的m个点极坐标第一项ρ′j,j=1,2,...,m,写成向量形式 
Figure FDA0000156818770000025
Figure FDA0000156818770000026
则,畸变退化多项式(1)中各次项系数 
Figure FDA0000156818770000027
为:
Figure FDA0000156818770000028
4)对系统检测时的每一幅子图像,将子图像不存在畸变时的理想子图像上每个像素P的以图像左上角为原点的直角坐标值(x,y)转换为以图像中心为原点的极坐标值(ρ,θ),代入式(1),获得其对应的畸变子图像像素坐标值(ρ′,θ′),再将此极坐标值转换为以图像左上角为原点的直角坐标值(u,v),若u,v为整数,则选取畸变子图像上对应像素P′(u,v)的灰度值作为理想子图像像素P的灰度值G(P);若u,v为非整数,则P为亚像素,采用畸变子图像上与(u,v)最为邻近的四个像素A(i,j),B(i+1,j),C(i,j+1),D(i+1,j+1),用四个像素的灰度值进行 双线性插值后的灰度值作为理想子图像上像素P′的灰度值G(P′);
设α=u-i,β=v-j
E为线段 
Figure FDA0000156818770000031
上一点,F为线段 
Figure FDA0000156818770000032
上一点,E,F也为亚像素, 
Figure FDA0000156818770000033
平行于 
Figure FDA0000156818770000034
和 
Figure FDA0000156818770000035
首先用第一次线性插值计算出E,F两处的灰度值G(E)和G(F)
G(E)=β[G(C)-G(A)]+G(A)          (3)
G(F)=β[G(D)-G(B)]+G(B)
式(3)中G(A),G(B),G(C),G(D)分别代表A,B,C,D四个像素的灰度值,同理,再用第二次线性插值计算出畸变子图像上像素P′的灰度值G(P′),即理想子图像上对应像素P灰度值G(P);
G(P)=G(P′)=α[G(F)-G(E)]+G(E)  (4) 。
CN 201210122901 2012-04-24 2012-04-24 超光滑表面缺陷检测系统的畸变校正方法 Active CN102661956B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN 201210122901 CN102661956B (zh) 2012-04-24 2012-04-24 超光滑表面缺陷检测系统的畸变校正方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN 201210122901 CN102661956B (zh) 2012-04-24 2012-04-24 超光滑表面缺陷检测系统的畸变校正方法

Publications (2)

Publication Number Publication Date
CN102661956A true CN102661956A (zh) 2012-09-12
CN102661956B CN102661956B (zh) 2013-06-12

Family

ID=46771488

Family Applications (1)

Application Number Title Priority Date Filing Date
CN 201210122901 Active CN102661956B (zh) 2012-04-24 2012-04-24 超光滑表面缺陷检测系统的畸变校正方法

Country Status (1)

Country Link
CN (1) CN102661956B (zh)

Cited By (14)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103438832A (zh) * 2013-08-30 2013-12-11 解则晓 基于线结构光的三维影像测量方法
CN104567590A (zh) * 2013-10-12 2015-04-29 上海和辉光电有限公司 一种缺陷定位部件
CN105115987A (zh) * 2015-09-30 2015-12-02 河南科技大学 基于数字滤波的圆锥滚子倒装缺陷检测方法
CN106447637A (zh) * 2016-09-30 2017-02-22 北京纳米维景科技有限公司 一种拼接探测器几何校正体模及校正方法
CN104655654B (zh) * 2015-02-04 2017-03-15 天津大学 基于尺寸畸变率的微细管道管内缺陷特征参数的校正方法
CN106645197A (zh) * 2016-12-29 2017-05-10 中国工程物理研究院激光聚变研究中心 检测精密光学元件表面颗粒物的在线检测系统及应用方法
CN106709885A (zh) * 2016-12-28 2017-05-24 大连鉴影光学科技有限公司 一种亚像素级畸变校正方法及装置
CN108152294A (zh) * 2017-12-26 2018-06-12 华中光电技术研究所(中国船舶重工集团公司第七七研究所) 一种超光滑镜片疵病检测装置及其方法
CN108508011A (zh) * 2017-02-24 2018-09-07 通用汽车环球科技运作有限责任公司 用于表征发动机缸体镗孔表面的方法及其装置
CN108645869A (zh) * 2018-08-20 2018-10-12 中国印刷科学技术研究院有限公司 凹印版辊表面缺陷智能检测的非缺陷排除方法及其装置
CN108760766A (zh) * 2018-05-25 2018-11-06 哈尔滨工业大学 一种大口径光学晶体表面微缺陷检测用的图像拼接方法
CN112304968A (zh) * 2020-09-11 2021-02-02 嘉兴驭光光电科技有限公司 用于微纳光学元件的检测系统及检测方法
CN112649985A (zh) * 2020-12-25 2021-04-13 武汉华星光电技术有限公司 一种显示面板及显示装置
CN114878583A (zh) * 2022-07-08 2022-08-09 四川大学 用于畸变光斑照明疵病暗场成像的图像处理方法及系统

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20040253616A1 (en) * 2003-06-12 2004-12-16 Cytyc Corporation Method and system of determining the stain quality of slides using a scatter plot distributions
CN1563957A (zh) * 2004-04-09 2005-01-12 浙江大学 光滑表面疵病的自动化检测方法及其系统
CN101086482A (zh) * 2007-07-17 2007-12-12 浙江大学 表面疵病检测图像拼接时物像坐标误差调整装置及方法
US20090224174A1 (en) * 2008-01-30 2009-09-10 Ralf Netz Calibration device and laser scanning microscope with such a calibration device
CN101639450A (zh) * 2009-08-28 2010-02-03 上海派尼科技实业有限公司 印刷质量的测试装置
CN102033308A (zh) * 2010-10-22 2011-04-27 浙江大学 一种超高分辨率的光学显微成像方法及装置

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20040253616A1 (en) * 2003-06-12 2004-12-16 Cytyc Corporation Method and system of determining the stain quality of slides using a scatter plot distributions
CN1563957A (zh) * 2004-04-09 2005-01-12 浙江大学 光滑表面疵病的自动化检测方法及其系统
CN101086482A (zh) * 2007-07-17 2007-12-12 浙江大学 表面疵病检测图像拼接时物像坐标误差调整装置及方法
US20090224174A1 (en) * 2008-01-30 2009-09-10 Ralf Netz Calibration device and laser scanning microscope with such a calibration device
CN101639450A (zh) * 2009-08-28 2010-02-03 上海派尼科技实业有限公司 印刷质量的测试装置
CN102033308A (zh) * 2010-10-22 2011-04-27 浙江大学 一种超高分辨率的光学显微成像方法及装置

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
杨甬英 等: "超光滑表面瑕疵的光学显微成像和数字化评价系统", 《红外与激光工程》, vol. 39, no. 2, 30 April 2010 (2010-04-30), pages 325 - 329 *

Cited By (21)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103438832B (zh) * 2013-08-30 2016-01-20 解则晓 基于线结构光的三维影像测量方法
CN103438832A (zh) * 2013-08-30 2013-12-11 解则晓 基于线结构光的三维影像测量方法
CN104567590A (zh) * 2013-10-12 2015-04-29 上海和辉光电有限公司 一种缺陷定位部件
CN104655654B (zh) * 2015-02-04 2017-03-15 天津大学 基于尺寸畸变率的微细管道管内缺陷特征参数的校正方法
CN105115987A (zh) * 2015-09-30 2015-12-02 河南科技大学 基于数字滤波的圆锥滚子倒装缺陷检测方法
CN105115987B (zh) * 2015-09-30 2017-11-14 河南科技大学 基于数字滤波的圆锥滚子倒装缺陷检测方法
CN106447637B (zh) * 2016-09-30 2019-03-15 北京纳米维景科技有限公司 一种拼接探测器几何校正体模及校正方法
CN106447637A (zh) * 2016-09-30 2017-02-22 北京纳米维景科技有限公司 一种拼接探测器几何校正体模及校正方法
CN106709885A (zh) * 2016-12-28 2017-05-24 大连鉴影光学科技有限公司 一种亚像素级畸变校正方法及装置
CN106645197A (zh) * 2016-12-29 2017-05-10 中国工程物理研究院激光聚变研究中心 检测精密光学元件表面颗粒物的在线检测系统及应用方法
CN106645197B (zh) * 2016-12-29 2024-01-30 中国工程物理研究院激光聚变研究中心 检测精密光学元件表面颗粒物的在线检测系统及应用方法
CN108508011A (zh) * 2017-02-24 2018-09-07 通用汽车环球科技运作有限责任公司 用于表征发动机缸体镗孔表面的方法及其装置
CN108152294A (zh) * 2017-12-26 2018-06-12 华中光电技术研究所(中国船舶重工集团公司第七七研究所) 一种超光滑镜片疵病检测装置及其方法
CN108760766A (zh) * 2018-05-25 2018-11-06 哈尔滨工业大学 一种大口径光学晶体表面微缺陷检测用的图像拼接方法
CN108760766B (zh) * 2018-05-25 2020-12-01 哈尔滨工业大学 一种大口径光学晶体表面微缺陷检测用的图像拼接方法
CN108645869A (zh) * 2018-08-20 2018-10-12 中国印刷科学技术研究院有限公司 凹印版辊表面缺陷智能检测的非缺陷排除方法及其装置
CN112304968A (zh) * 2020-09-11 2021-02-02 嘉兴驭光光电科技有限公司 用于微纳光学元件的检测系统及检测方法
CN112649985A (zh) * 2020-12-25 2021-04-13 武汉华星光电技术有限公司 一种显示面板及显示装置
CN112649985B (zh) * 2020-12-25 2022-02-22 武汉华星光电技术有限公司 一种显示面板及显示装置
CN114878583A (zh) * 2022-07-08 2022-08-09 四川大学 用于畸变光斑照明疵病暗场成像的图像处理方法及系统
CN114878583B (zh) * 2022-07-08 2022-09-20 四川大学 用于畸变光斑照明疵病暗场成像的图像处理方法及系统

Also Published As

Publication number Publication date
CN102661956B (zh) 2013-06-12

Similar Documents

Publication Publication Date Title
CN102661956B (zh) 超光滑表面缺陷检测系统的畸变校正方法
CN106856003B (zh) 轴类工件侧表面缺陷检测图像的展开校正方法
Huang et al. Accurate projector calibration based on a new point-to-point mapping relationship between the camera and projector images
CN101996398B (zh) 用于晶圆对准的图像匹配方法及设备
CN103615980B (zh) 一种板件上圆孔参数的测量方法及系统
CN109751964B (zh) 一种高精度非接触式管径测量方法及装置
CN111707187B (zh) 一种大型零件的测量方法及系统
EP3123501A1 (en) Variable image field curvature for object inspection
CN102322813A (zh) 三维网格应变测量方法
Zou et al. High-accuracy calibration of line-structured light vision sensors using a plane mirror
CN103531500A (zh) 晶圆缺陷检测设备的校准方法
CN114331924B (zh) 大工件多相机视觉测量方法
Zhou et al. Omnidirectional stereo vision sensor based on single camera and catoptric system
CN109544642B (zh) 一种基于n型靶标的tdi-ccd相机参数标定方法
CN102680477A (zh) 大尺寸光学元件高精度调平方法与装置
CN106780391A (zh) 一种用于全视角三维测量仪光学系统的畸变矫正算法
CN107228864A (zh) 显示面板表面缺陷的检测系统
CN104168414A (zh) 拍摄并拼接物件影像的方法
Liu et al. Large-curvature specular surface phase measuring deflectometry with a curved screen
CN102003941B (zh) 基于视觉的大型系泊链五环长测量方法及其测量装置
CN110864878A (zh) 一种高效大视场平显系统显示畸变检测的方法
CN103575742A (zh) 一键式光学测量仪
US20130162816A1 (en) Device for measuring the shape of a mirror or of a specular surface
KR101833245B1 (ko) 위상측정 광선편향법을 이용한 표면 검사 장치의 교정 방법 및 시스템
CN115046498A (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
EE01 Entry into force of recordation of patent licensing contract

Application publication date: 20120912

Assignee: Hangzhou Zheda Three Color Instrument Co., Ltd.

Assignor: Zhejiang University

Contract record no.: 2014330000345

Denomination of invention: Super-smooth surface defect detection system and distortion correction method thereof

Granted publication date: 20130612

License type: Common License

Record date: 20140901

LICC Enforcement, change and cancellation of record of contracts on the licence for exploitation of a patent or utility model