CN102538706B - 绘制表面轮廓的方法和装置 - Google Patents

绘制表面轮廓的方法和装置 Download PDF

Info

Publication number
CN102538706B
CN102538706B CN201110376531.7A CN201110376531A CN102538706B CN 102538706 B CN102538706 B CN 102538706B CN 201110376531 A CN201110376531 A CN 201110376531A CN 102538706 B CN102538706 B CN 102538706B
Authority
CN
China
Prior art keywords
index
height
vector field
logical value
target object
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.)
Active
Application number
CN201110376531.7A
Other languages
English (en)
Other versions
CN102538706A (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.)
Gaoying Technology Co ltd
Original Assignee
Koh Young Technology Inc
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 Koh Young Technology Inc filed Critical Koh Young Technology Inc
Publication of CN102538706A publication Critical patent/CN102538706A/zh
Application granted granted Critical
Publication of CN102538706B publication Critical patent/CN102538706B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/50Depth or shape recovery
    • G06T7/521Depth or shape recovery from laser ranging, e.g. using interferometry; from the projection of structured light
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01BMEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
    • G01B11/00Measuring arrangements characterised by the use of optical techniques
    • G01B11/24Measuring arrangements characterised by the use of optical techniques for measuring contours or curvatures
    • G01B11/25Measuring arrangements characterised by the use of optical techniques for measuring contours or curvatures by projecting a pattern, e.g. one or more lines, moiré fringes on the object
    • G01B11/2518Projection by scanning of the object

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • General Physics & Mathematics (AREA)
  • Optics & Photonics (AREA)
  • Theoretical Computer Science (AREA)
  • Length Measuring Devices By Optical Means (AREA)

Abstract

本发明公开一种绘制表面轮廓的方法和装置。该方法包括:沿至少两个方向依次地向目标物体投射狭缝图案光,以获得由目标物体反射的图案图像,根据方向使用图案图像获得高度,获得示出高度的最大变化方向的矢量场,获得与至少两个方向对应的高度的可靠性指数,使用可靠性指数和矢量场获得整合的矢量场,以及,使用整合的矢量场计算目标物体的每个位置的高度。因此,提高了精确度。

Description

绘制表面轮廓的方法和装置
相关申请的交叉引用
本申请要求2010年11月19日提交的韩国专利申请No.10-2010-0115351的优先权,为了在此充分列举的目的,该专利以引用的方式并入本文。
技术领域
本发明的示例性实施例涉及用于绘制表面轮廓的方法和装置,或更具体地说,涉及用于绘制表面轮廓的非接触式方法和非接触式装置。
背景技术
已经以各种方式研发出了各种绘制三维目标物体的表面轮廓的非接触式方法。在非接触式方法中,经常使用利用激光的方法和利用莫尔图案(moirépattern)的方法。
根据使用激光的方法或使用莫尔图案的方法,激光或狭缝图案光沿倾斜方向照射目标物体。然而,当目标物体以大于激光或狭缝图案光照射角度的角度突出时,可能会产生激光或狭缝图案光无法到达的阴影区域。
为了解决上述问题,激光或狭缝图案光沿各种方向照射到基板上,以通过获取基板的每个位置(x,y)的高度来获取目标物体的轮廓。然而,通过照射激光或狭缝图案光所获得的高度通常不一致,并且由于高度包含缺陷,导致难以获得精确的高度。
发明内容
本发明的示例性实施例提供了一种绘制表面轮廓的方法,其能够提高高度测量的精确度。
本发明的示例性实施例还提供了能够实现上述方法的装置。
本发明的附加特征将在下述说明中阐明,并且部分将从说明中明显地得出,或者可以通过实施本发明而获知。
本发明的示例性实施例披露了一种绘制表面轮廓的方法,包括:沿至少两个方向依次地向目标物体投射狭缝图案光,以获得由目标物体反射的图案图像,根据方向使用图案图像获得高度,获得示出高度的最大变化方向的矢量场,获得与至少两个方向对应的高度的可靠性指数,使用可靠性指数和矢量场获得整合的矢量场,以及,使用整合的矢量场计算目标物体的每个位置的高度。
例如,可靠性指数可以包括可见度指数、反射指数和阴影指数中的至少一者。
例如,当可见度大于可见度阈值时,可见度指数可以为逻辑值真(1);而当可见度等于或小于可见度阈值时,可见度指数可以为逻辑值假(0)。当通过累加大于反射阈值的强度的逻辑值所生成的累积强度大于或等于整合的反射阈值时,反射指数可以为逻辑值真(1);而当累积强度小于整合的反射阈值时,反射指数可以为逻辑值假(0)。当通过累加小于阴影阈值的强度的逻辑值所生成的累积强度大于或等于整合的阴影阈值时,阴影指数可以为逻辑值真(1);而当累积强度小于整合的阴影阈值时,阴影指数可以为逻辑值假(0)。与至少两个方向对应的高度的可靠性指数可以是使用可见度指数、反射指数和阴影指数获得的。
例如,与至少两个方向对应的高度的可靠性指数可以是通过将可见度指数、反射指数和阴影指数逻辑相加而获得的。
例如,与至少两个方向对应的高度的可靠性指数是通过将可见度指数与可选的逻辑值真或假相乘,将反射指数与可选的逻辑值相乘,将阴影指数与可选的逻辑值相乘,再将上述三个乘积逻辑相加而获得的。
例如,整合的矢量场可以是通过选择可靠性指数的逻辑值为真的矢量场以及对所选择的矢量场进行取中位数操作而获得的。
例如,计算目标物体的每个位置的高度可以是这样执行的:列泊松方程,其中方程右侧通过对整合的矢量场应用散度算子而获得,方程的左侧通过对高度应用拉普拉斯算子而获得;以及,求解泊松方程,以得到与高度对应的解。
例如,求解泊松方程可以是这样执行的:将具有连续算子 的泊松方程转换为具有离散算子的离散泊松方程;以及,利用傅里叶级数表示离散泊松方程的解,以得到傅里叶级数的系数。
本发明的示例性实施例披露了一种绘制表面轮廓的装置,包括:支撑部分,其支撑设有目标物体的基板;图案图像投射部分,其沿至少两个方向向目标物体投射狭缝图案光;图像捕获部分,其捕获依次被目标物体反射的与至少两个方向对应的图案图像;存储单元,其存储由图像捕获部分捕获的图案图像;以及,计算部分,其使用图案图像获得与至少两个方向对应的高度,通过对高度进行梯度操作获得与至少两个方向对应的矢量场,获得与至少两个方向对应的高度的可靠性指数,通过使用可靠性指数和矢量场获得整合的矢量场,并且使用整合的矢量场计算目标物体的每个位置的高度。
例如,计算部分可以包括:高度计算单元,其依次使用图案图像计算与至少两个方向对应的高度;可靠性指数计算单元,其计算与至少两个方向对应的高度的可靠性指数;存储单元,其存储计算出的高度和计算出的可靠性指数;矢量场计算单元,其计算矢量场;整合的矢量场计算单元,其使用与至少两个方向对应的可靠性指数和矢量场计算整合的矢量场;边界设定单元,其设定目标物体的边界;以及泊松方程求解单元,其通过对由边界设定单元设定的边界应用边界条件来求解泊松方程,泊松方程的右侧通过对整合的矢量场应用散度算子而获得,泊松方程的左侧通过对高度应用拉普拉斯算子而获得。
应当理解,前述概括说明和以下的详细说明是示例性和说明性的,其目的在于为本发明提供进一步的说明。
附图说明
为提供对本发明的进一步的理解,本发明包含并入说明书并且构成说明书的一部分的附图,附图示出本发明的实施例并且连同说明书一起用于解释本发明的原理。
图1是示出根据本发明的示例性实施例的绘制表面轮廓的方法的流程图。
图2是示出根据本发明的示例性实施例的绘制表面轮廓的装置的示意图。
图3是示出根据示例性实施例的图2中的存储单元和计算部分的框图。
具体实施方式
在下文中,将参考示出了本发明的实施例的附图更详细描述本发明。然而,本发明可以体现为多种不同的形式,并且不应解释为仅局限于在此列举的实施例。相反,提供这些示例性实施例,使得本公开完整,并且将本发明的范围充分地传达给本领域的技术人员。在附图中,为了清晰可见,可能夸大层和区域的尺寸和相对尺寸。图中相似的附图标记代表相似的元件。
下面,将参考附图详细地描述本发明。
图1是示出根据本发明的示例性实施例的绘制表面轮廓的方法的流程图。
参考图1,根据本发明的绘制表面轮廓的方法,狭缝图案光依次沿至少两个方向照射到目标物体上,并且依次获得被目标物体反射的目标物体的图案图像(步骤S110)。
例如,狭缝图案光可以沿斜角照射到目标物体上,并且可以通过使光传送光栅图案来生成狭缝图案光,该光栅图案上选择性地形成有透明部分和不透明部分。当狭缝图案光沿一个方向照射到目标物体上以测量该目标的高度时,狭缝图案光可能不会到达目标物体的相反部分,从而产生了阴影区域。在这种情况下,狭缝图案光可以沿其他方向照射到目标物体上,以便补偿目标物体的相反部分。因此,可以测量得到更精确的高度。
在照射狭缝图案光时,例如沿每个方向移动光栅图案而获得图案图像。例如,在3桶算法中,通过将光栅图案移动与光栅图案的栅距的1/3相等的距离,获得3个图案图像;并且在4桶算法中,通过将光栅图案移动与光栅图案的栅距的1/4相等的距离,获得4个图案图像。
然后,使用每个方向上的图案图像计算目标物体的高度(步骤S120)。
根据穿过光栅图案的狭缝图案光的图案图像可以近似为下述表达式1所示的正弦函数。
表达式1
Ibi=A+I0sin(φibi)
在上述表达式中,下标‘b’表示多个方向中的一个方向(或通道),下标‘i’与光栅图案被移动(i-1)次时的第i个图像对应(例如,在N桶算法中,‘i’可以是从1到N的一个自然数),‘I0’是图案光的最大强度(或调幅),‘A’是N桶算法中图案图像的平均强度,而因移动光栅图案而形成的相位αbi可以由下述表达式2表示。
表达式2
α bi = 2 π ( i - 1 ) N
使用上述表达式1,沿每个方向b的相位可以由下述表达式3表示。
表达式3
φ b = tan - 1 [ Σ i = 1 N I bi sin ( α bi ) Σ i = 1 N I bi cos ( α bi ) ]
使用上述表达式3,在每个方向b上的高度Hb可以近似为下述表达式4。
表达式4
H b ( x , y ) = l 0 2 π f 0 d 0 φ b ( x , y )
在上述表达式中,l0是目标物体与用于捕捉目标物体的图像的图像捕捉部分之间的距离,d0是图像捕捉部分与投射狭缝图案光的图案图像投射部分之间的距离,而f0是投射到基准面上的正弦图案(sin fridge pattern)的频率。
使用上述表达式4,可以获得每个方向(b)上的目标物体的每个点(x,y)的高度Hb
另一方面,上述获得每个方向(b)上的目标物体的每个点(x,y)的高度Hb的方法是一个实例,并且可以使用其他的方法得到高度Hb。
下面,将详细描述从每个方向(b)上的目标物体的每个点(x,y)的高度Hb中去除噪声以及将每个方向上(b)的高度Hb整合为一个高度的方法。
在每个方向(b)上获得的高度[Hb(x,y)]可以包括例如各种类型的缺陷,诸如由阴影引起的缺陷,由全反射引起的缺陷(或反射缺陷)等。因此,在每个方向(b)上获得的高度[Hb(x,y)]可以彼此不同。因此,可以通过整合各个方向上的高度[H(x,y)]获得精确的高度。
为了该目的,获取示出高度的最大变化方向的矢量场(步骤S130)。
可以使用下述表达式5获取在每个方向(b)上的矢量场(Gb(x,y))。
表达式5
G b ( x , y ) = ▿ H b = ( ∂ H b ( x , y ) ∂ x , ∂ H b ( x , y ) ∂ y )
也就是说,通过对每个方向(b)上的高度Hb进行梯度操作,获得每个方向(b)上的矢量场(Gb(x,y))。
然后,获得关于每个方向(b)上的高度Hb的可靠性指数(Cb)(步骤S140)。可靠性指数(Cb)可以包括可见度指数、反射指数和阴影指数中的至少一者。
为了获得可见度指数,在4桶算法的情况下,表达式1中的最大强度(或调幅)I0可以由下述表达式6表示。
表达式6
I 0 = ( I b 1 - I b 3 ) 2 + ( I b 2 - I b 4 ) 2 2
此外,表达式1中的平均强度A可以由下述表达式7表示。
表达式7
A b = 1 N Σ i = 1 N I bi
然后,可见度VF可以由使用表达式6和表达式7获得的下述表达式8表示。
表达式8
VF b = I 0 A b
根据本示例性实施例的可见度指数VIb可以由使用可见度VF获得的下述表达式9表示。
表达式9
VIb=(VFb>vthr)
在上述表达式中,vthr是可见度的阈值,并且该可见度阈值vthr可以通过实验来确定。表达式9的含义如下。即,例如,当多个方向(b)中的第一方向的可见度指数(VI1)大于阈值(vthr)时,将逻辑值真(1)赋予第一方向的可见度指数(VI1),而当第一方向的可见度指数(VI1)小于或等于阈值(vthr)时,将逻辑值假(0)赋予可见度指数(VI1)。
此外,与全反射引起的反射缺陷有关的多个方向(b)上的整合的反射指数(SpIb)可以由下述表达式10表示。
表达式10
SpI b = ( Σ i = 1 N ( I bi > sp thr ) ≥ Sp thr )
在上述表达式中,反射阈值(spthr)和整合的反射阈值(Spthr)可以通过实验来确定。表达式10的含义如下。即,例如,当多个方向中的第一方向的第i个强度(I1i)(N桶算法中的第i个强度)大于反射阈值(spthr)时,将逻辑值真(1)赋予第i个强度,并且当第i个强度(I1i)小于或等于反射阈值(spthr)时,将逻辑值假(0)赋予第i个强度(I1i)。然后,将有关于N桶算法的N个图像的逻辑值相加。当逻辑值的总和大于或等于整合的反射阈值(Spthr)时,将逻辑值真(1)赋予第一方向的反射指数(SpI1),而当逻辑值的总和小于整合的反射阈值(Spthr)时,将逻辑值假(0)赋予第一方向的反射指数(SpI1)。
此外,与阴影引起的阴影缺陷有关的多个方向(b)上的整合的阴影指数(ShIb)可以由下述表达式11表示。
表达式11
ShI b = ( &Sigma; i = 1 N ( I bi < sh thr ) &GreaterEqual; Sh thr )
在上述表达式中,阴影阈值(shthr)和整合的阴影阈值(Shthr)可以通过实验来确定。表达式11的含义如下。即,例如,当第一方向的N桶算法的第i个强度(I1i)小于阴影阈值(shthr)时,将逻辑值真(1)赋予第i个强度(I1i),并且当第i个强度(I1i)大于或等于阴影阈值(shthr)时,将逻辑值假(0)赋予第i个强度(I1i)。然后,将有关于N桶算法的N个图像的逻辑值相加。当逻辑值的总和大于或等于整合的阴影阈值(Shthr)时,将逻辑值真(1)赋予第一方向的阴影指数(ShI1),而当逻辑值的总和小于整合的阴影阈值(Shthr)时,将逻辑值假(0)赋予第一方向的阴影阈值(ShI1)。
利用分别从上述表达式9、表达式10和表达式11获得的可见度指数(VIb)、反射指数(SpIb)和阴影指数(ShIb),使用下述表达式12计算每个方向(b)的可靠性指数(Cb)。
表达式12
Cb=VM·VIb|SpM·SpIb|ShM·Shb
在上述表达式中,VM、SpM和ShM具有逻辑值真(1)或假(0),并且‘|’运算符表示“或”操作。如上文所述,通过控制VM、SpM和ShM,可以选择可靠性指数(Cb)所需的项。例如,当将逻辑值真(1)赋予全部VM、SpM和ShM时,考虑全部可见度、反射缺陷和阴影缺陷获得每个方向的可靠性指数(Cb)。当将逻辑值假(0)赋予VM和SpM,并且将逻辑值真(1)赋予ShM时,仅考虑阴影缺陷获得每个方向的可靠性质数(Cb)。
在本实施例中,例如,首先获得示出高度的最大变化和最大变化方向的矢量场(步骤S130),然后获得每个方向的可靠性指数(步骤S140)。然而,上述执行顺序是示例性的。因此,可以首先获得每个方向的可靠性指数,然后再获得矢量场。
然后,利用使用表达式12获得的可靠性指数(Cb)和使用表达式5获得的矢量场(Gb(x,y)),获得整合的矢量场(G(x,y))(步骤S150)。例如,可以使用下述表达式13获得整合的矢量场(G(x,y))。
表达式13
G(x,y)=Median(G1·C1,G2·C2,...,GM·CM)
即,在全部方向的矢量场中,选择逻辑值为真(1)的矢量场,然后对逻辑值为真(1)的矢量场进行取中位数操作以获得整合的矢量场(G(x,y))。
然后,利用使用表达式13获得的整合的矢量场(G(x,y))获得目标物体的每个位置(x,y)的高度(H(x,y))(步骤S160),高度(H(x,y))合并了全部方向(b)的全部高度(Hb(x,y))。为了该目的,将下述表达式14最小化以获得高度(H(x,y))。
表达式14
| &dtri; H - G |
为了将表达式14最小化,可以使用表示为下述表达式15的泊松方程。
表达式15
&dtri; 2 H = &dtri; &CenterDot; G
如上述表达式15所述,将散度算子应用于整合的矢量场以设定方程的右侧,将拉普拉斯算子应用于高度(H(x,y))以设定方程的左侧。然后,满足上述泊松方程的高度(H(x,y))是最终的目标。通过沿方向(b)投射狭缝图案所获得的高度(Hb(x,y))包括诸如因阴影、全反射等引起的各种缺陷。然而,从整合的矢量场(G(x,y))中获得了没有上述缺陷的上述表达式15的高度(H(x,y))。因此,高度(H(x,y))不包括干扰。
例如,为了解出上述表达式15以获得高度(H(x,y))的解,例如将具有拉普拉斯算子中的连续算子和散度算子的连续泊松方程转换为(或近似为)具有离散算子的离散泊松方程。离散泊松方程可以由下述表达式16表示。
表达式16
h m + 1 , n - 2 h m , n + h m - 1 , n &Delta;x 2 + h m , n + 1 - 2 h m , n + h m , n - 1 &Delta;y 2 = g m , n
为了解上述偏微分方程,将边界的高度设定为地。因此,应用h0n=hMn=hm0=hmN=0的边界条件。
在上述表达式16中,当假设沿X轴的单位长度和沿y轴的单位长度等于1(Δx=Δy=1)时,表达式16可以表示为下述表达式17。
表达式17
hm+1,n+hm-1,n+hm,n+1+hm,n-1-4hm,n=gm,n
另一方面,可以使用下述表达式18所示的傅里叶级数将任意两个参数函数展开。
表达式18
h m , n = &Sigma; j = 1 M - 1 &Sigma; k = 1 N - 1 U j , k sin ( &pi;mj M ) sin ( &pi;nk N )
当将上述表达式18所表示的傅里叶级数假设为由表达式17所表示的方程的解并应用于表达式17时,可以得到下述表达式19的方程。
表达式19
&Sigma; j = 1 M - 1 &Sigma; k = 1 N - 1 U j , k [ sin ( &pi; ( m + 1 ) j M ) + sin ( &pi; ( m - 1 ) j M ] sin ( &pi;nk N ) + [ sin ( &pi; ( n + 1 ) k N ) + sin ( &pi; ( n - 1 ) k N ] sin ( &pi;mj M ) - 4 sin ( &pi;mj M ) sin ( &pi;nk N )
= &Sigma; j = 1 M - 1 &Sigma; k = 1 N - 1 F jk sin ( &pi;mj M ) sin ( &pi;nk N )
在上述表达式中,方程右侧的系数Fjk可以由下述表达式20表示。
表达式20
F jk = 1 MN &Sigma; m = 1 M - 1 &Sigma; n = 1 N - 1 g mn sin ( &pi;mj M ) sin ( &pi;nk N )
在上述表达式20中,可以通过将散度算子应用于表达式15的整合的矢量场并进行傅里叶展开来获得系数gmn
可以使用正弦函数的可加性公式将表达式19表示为下述表达式21。
表达式21
&Sigma; j = 1 M - 1 &Sigma; k = 1 N - 1 [ ( 2 cos ( &pi; j M ) - 2 cos ( &pi; k N ) - 4 ) U jk - F jk ] sin ( &pi;mj M ) sin ( &pi;nk N ) = 0
因此,表达式18的系数Ujk可以由下述表达式22表示,从而可以获得表达式15的高度的解。
表达式22
U jk = F jk 2 cos ( &pi;j M ) + 2 cos ( &pi;k N ) - 4
根据本示例性实施例,在通过沿各个方向应用狭缝图案所获得的高度中,选择出没有诸如阴影缺陷、反射缺陷等缺陷的高度。并且,将没有缺陷的高度整合以提高精确度。
图2是示出根据本发明的示例性实施例的绘制表面轮廓的装置的示意图,而图3是示出根据示例性实施例的图2中的存储单元和计算部分的框图。
参考图2,根据本发明的示例性实施例的绘制表面100的轮廓的装置包括:支撑部分140,其支撑并承载设有目标物体的基板S;至少两个图案图像投射部分110,其向基板S投射图案图像;图像捕获部分130,其捕获基板S的图像;存储单元120,其存储捕获的图像;以及计算部分150,其执行各种计算。
图案图像投射部分110将图案图像投射到基板S上以获取基板S上的目标物体的轮廓。例如,图案图像投射部分110包括:光源112,其产生光;光栅114,其将光源112产生的光转换为图案图像;光栅移动器116,其移动光栅114;以及透镜118,其向目标物体投射由光栅114转换的图案图像。
可以利用诸如压电致动器(PZT)等光栅移动器116使光栅114移动2π/N以改变图案图像,其中‘N’是大于1的自然数。绘制表面100的轮廓的装置可以包括相对于图像捕获部分130沿着圆周设置的多个图案图像投射部分110,以提高精确度。多个图案图像投射部分110相对于基板S倾斜地设置,以便沿多个方向将图案图像投射到基板S上。
图像捕获部分130捕获其上由图案图像投射部分110投射有图案图像的基板S的图像。例如,图像捕获部分130可以直接设置在基板S的上方。图像捕获部分130可以包括用于捕获图像的CCD照相机或CMOS照相机。
存储单元120存储由图像捕获部分130捕获的图像。存储单元120包括多个部分,并且在多个方向中的一个方向上根据N桶算法所捕获的N个图像可以存储在其中一个部分里。
计算部分150执行各种操作,以使用储存在存储单元120中的图案图像、示出高度的最大变化方向的矢量场、与方向对应的高度的可靠性指数获得与各个方向对应的高度;使用可靠性指数和矢量场获得整合的矢量场;然后使用整合的矢量场获得目标物体的每个位置的高度。
参考图3,执行上述操作的计算部分150可以包括高度计算单元151、可靠性指数计算单元152、存储单元153、矢量场计算单元154、整合的矢量场计算单元155、边界设定单元157和泊松方程求解单元156。
高度计算单元151使用图案图像分别计算与各种方向对应的高度。具体地说,例如使用通过在与第一方向对应的第一通道(Ch1)上投射狭缝图案图像而与N桶算法对应的N个图像(I11、I12,...,I1N),通过表达式4计算与第一方向对应的高度H1,以及,例如使用通过在与第m个方向对应的第m通道(Chm)上投射狭缝图案图像而与N桶算法对应的N个图像(Im1、Im2,...,ImN),通过表达式4计算与第m方向对应的高度Hm
可靠性指数计算单元152通过表达式12计算与各个方向(1至m)的高度(H1,H2,...Hm)相关的可靠性指数(C1,C2,...,Cm)。
存储单元153存储计算出的高度(H1,H2,...Hm)和计算出的可靠性指数(C1,C2,...,Cm)。
矢量场计算单元154使用存储在存储单元153里的高度(H1,H2,...Hm)通过表达式5计算示出最大变化方向的矢量场(G1,G2,...Gm)。
整合的矢量场计算单元155使用与方向(1至m)对应的可靠性指数(C1,C2,...,Cm)和矢量场(G1,G2,...Gm),通过表达式13计算整合的矢量场(G)。
边界设定单元157设定目标物体的边界。边界设定单元157可以使用由下述表达式23表示的加权平均方法设定目标物体的边界。
表达式23
h i = &Sigma; j = 0 m h j &CenterDot; w j &Sigma; j = 0 m w j
在上述表达式23中,wj是加权因子。在这种情况下,边界不限于目标物体的边缘,而是可以向设置有目标物体的基板扩展。
泊松方程求解单元156通过应用从边界设定单元157获得的边界条件来求解表达式15的泊松方程以得到目标物体的高度,该方程的右侧通过将散度算子应用于整合的矢量场(G)而设定,该方程的左侧通过将拉普拉斯算子应用于高度而设定。
已参考图1详细描述了每个单元的操作方法。因此,将省略对其进一步描述。
根据本发明的示例性实施例的绘制表面100的轮廓的装置利用图案图像投射部分110向基板S上的目标物体投射狭缝图案,并捕获被目标物体反射的图像以获得目标物体的轮廓。然而,图2的绘制表面100的轮廓的装置是示例性的。根据本发明的绘制表面轮廓的装置可以具有各种硬件的变化。例如,绘制表面轮廓的装置可以具有两个以上的图案图像投射部分110并且图案捕获部分具有不同的元件。
对本领域的技术人员显而易见的是,在不脱离本发明的精神和范围的情况下,可以对本发明做出各种修改和变化。因此,本发明旨在覆盖落入所附权利要求及其等同内容的范围内的本发明的修改和变化。

Claims (10)

1.一种绘制表面轮廓的方法,包括:
沿至少两个方向依次地向目标物体投射狭缝图案光,以获得由所述目标物体反射的图案图像;
根据所述方向使用图案图像获得高度;
获得示出高度的最大变化方向的矢量场;
获得与至少两个方向对应的高度的可靠性指数;
使用所述可靠性指数和所述矢量场获得整合的矢量场;以及
使用所述整合的矢量场计算所述目标物体的每个位置的高度。
2.根据权利要求1所述的绘制表面轮廓的方法,其中,所述可靠性指数包括可见度指数、反射指数和阴影指数中的至少一者。
3.根据权利要求2所述的绘制表面轮廓的方法,其中,当可见度大于可见度阈值时,所述可见度指数为逻辑值真(1);而当所述可见度等于或小于所述可见度阈值时,所述可见度指数为逻辑值假(0),
当通过将大于反射阈值的强度的逻辑值累加所生成的累积强度大于或等于整合的反射阈值时,所述反射指数为逻辑值真(1);而当通过将大于反射阈值的强度的逻辑值累加所生成的累积强度小于所述整合的反射阈值时,所述反射指数为逻辑值假(0),
当通过将小于阴影阈值的强度的逻辑值累加所生成的累积强度大于或等于整合的阴影阈值时,所述阴影指数为逻辑值真(1);而当通过将小于阴影阈值的强度的逻辑值累加所生成的累积强度小于所述整合的阴影阈值时,所述阴影指数为逻辑值假(0),并且
与至少两个方向对应的高度的可靠性指数是使用所述可见度指数、所述反射指数和所述阴影指数获得的。
4.根据权利要求3所述的绘制表面轮廓的方法,其中与至少两个方向对应的高度的可靠性指数是通过将所述可见度指数、所述反射指数和所述阴影指数逻辑相加而获得的。
5.根据权利要求3所述的绘制表面轮廓的方法,其中,与至少两个方向对应的高度的可靠性指数是通过将所述可见度指数与可选的逻辑值真或假相乘,将所述反射指数与可选的逻辑值相乘,将所述阴影指数与可选的逻辑值相乘,再将上述三个乘积逻辑相加而获得的。
6.根据权利要求5所述的绘制表面轮廓的方法,其中这样获得所述整合的矢量场:
选择可靠性指数的逻辑值为真的矢量场;以及
通过对所选择的矢量场进行取中位数操作获得所述整合的矢量场。
7.根据权利要求6所述的绘制表面轮廓的方法,其中,计算所述目标物体的每个位置的高度是这样执行的:
列泊松方程,其中方程右侧通过对所述整合的矢量场应用散度算子而获得,以及方程的左侧通过对所述高度应用拉普拉斯算子而获得;以及
求解所述泊松方程,以得到与所述高度对应的解。
8.根据权利要求7所述的绘制表面轮廓的方法,其中,求解所述泊松方程是这样执行的:
将具有连续算子的所述泊松方程转换为具有离散算子的离散泊松方程;以及
利用傅里叶级数表示所述离散泊松方程的解,以得到所述傅里叶级数的系数。
9.一种绘制表面轮廓的装置,包括:
支撑部分,其支撑设有目标物体的基板;
图案图像投射部分,其沿至少两个方向朝所述目标物体投射狭缝图案光;
图像捕获部分,其捕获依次被所述目标物体反射的与所述至少两个方向对应的图案图像;
存储单元,其存储由所述图像捕获部分捕获的所述图案图像;以及
计算部分,其使用所述图案图像获得与至少两个方向对应的高度,通过对所述高度进行梯度操作获得与所述至少两个方向对应的矢量场,获得与至少两个方向对应的高度的可靠性指数,通过使用所述可靠性指数和所述矢量场获得整合的矢量场,并且使用所述整合的矢量场计算所述目标物体的每个位置的高度。
10.根据权利要求9所述的绘制表面轮廓的装置,其中,所述计算部分包括:
高度计算单元,其依次使用所述图案图像计算与所述至少两个方向对应的所述高度;
可靠性指数计算单元,其计算与所述至少两个方向对应的高度的可靠性指数;
存储单元,其存储计算出的高度和计算出的可靠性指数;
矢量场计算单元,其计算所述矢量场;
整合的矢量场计算单元,其使用与所述至少两个方向对应的所述可靠性指数和所述矢量场计算所述整合的矢量场;
边界设定单元,其设定所述目标物体的边界;以及
泊松方程求解单元,其通过对由所述边界设定单元设定的所述边界应用边界条件来求解所述泊松方程,所述泊松方程的右侧通过对所述整合的矢量场应用散度算子而获得,以及所述泊松方程的左侧通过对所述高度应用拉普拉斯算子而获得。
CN201110376531.7A 2010-11-19 2011-11-21 绘制表面轮廓的方法和装置 Active CN102538706B (zh)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
KR1020100115351A KR101174676B1 (ko) 2010-11-19 2010-11-19 표면형상 측정방법 및 측정장치
KR10-2010-0115351 2010-11-19

Publications (2)

Publication Number Publication Date
CN102538706A CN102538706A (zh) 2012-07-04
CN102538706B true CN102538706B (zh) 2014-09-24

Family

ID=46064019

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201110376531.7A Active CN102538706B (zh) 2010-11-19 2011-11-21 绘制表面轮廓的方法和装置

Country Status (4)

Country Link
US (1) US8982330B2 (zh)
JP (1) JP5546522B2 (zh)
KR (1) KR101174676B1 (zh)
CN (1) CN102538706B (zh)

Families Citing this family (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103105285B (zh) * 2013-01-22 2015-03-18 中国计量学院 基于暗影剔除的光纤面板极限分辨率测试方法及仪器
KR101490359B1 (ko) * 2013-05-20 2015-02-06 (주)펨트론 표면 형상 측정 장치
JP6507653B2 (ja) * 2015-01-13 2019-05-08 オムロン株式会社 検査装置及び検査装置の制御方法
CN104713497B (zh) * 2015-03-13 2017-07-04 香港应用科技研究院有限公司 相位偏移校准方法、3d形状检测的方法、系统及投影系统
WO2016145582A1 (zh) * 2015-03-13 2016-09-22 香港应用科技研究院有限公司 相位偏移校准方法、3d形状检测的方法、系统及投影系统
CN105043300A (zh) * 2015-05-27 2015-11-11 东莞市盟拓光电科技有限公司 可从多方位对被测物同时投影的三维测量系统
JP7040325B2 (ja) * 2018-07-04 2022-03-23 日本製鉄株式会社 凹凸部検出方法及び凹凸部検出装置
US11610350B2 (en) 2021-08-05 2023-03-21 Adobe Inc. Automatically filling a geometric shape with graphical cells aligned with the geometric shape
US11651536B2 (en) * 2021-09-16 2023-05-16 Adobe Inc. Directional pattern generation for digital images
US11704848B2 (en) 2021-10-19 2023-07-18 Adobe Inc. Systems for generating stroked paths
CN118742784A (zh) 2022-03-09 2024-10-01 雅马哈发动机株式会社 三维计测用运算装置、三维计测用程序、记录介质、三维计测装置及三维计测用运算方法

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1437000A (zh) * 2002-02-09 2003-08-20 沈阳同联集团高新技术有限公司 投影栅线测量物体三维表面形状的方法和装置
JP2007322162A (ja) * 2006-05-30 2007-12-13 3D Media Co Ltd 3次元形状測定装置及び3次元形状測定方法

Family Cites Families (18)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4087788A (en) * 1977-01-14 1978-05-02 Ncr Canada Ltd - Ncr Canada Ltee Data compression system
US5300347A (en) * 1991-03-01 1994-04-05 Kimberly-Clark Corporation Embossed facial tissue
US7440590B1 (en) * 2002-05-21 2008-10-21 University Of Kentucky Research Foundation System and technique for retrieving depth information about a surface by projecting a composite image of modulated light patterns
US7286246B2 (en) * 2003-03-31 2007-10-23 Mitutoyo Corporation Method and apparatus for non-contact three-dimensional surface measurement
JP4077754B2 (ja) 2003-04-04 2008-04-23 オリンパス株式会社 3次元形状測定装置
US8675176B2 (en) * 2005-02-25 2014-03-18 Asml Netherlands B.V. Parameter control in a lithographic apparatus using polarization
JP4611782B2 (ja) 2005-03-28 2011-01-12 シチズンホールディングス株式会社 3次元形状測定方法及び測定装置
JP4230525B2 (ja) * 2005-05-12 2009-02-25 有限会社テクノドリーム二十一 3次元形状計測方法およびその装置
CA2528791A1 (en) * 2005-12-01 2007-06-01 Peirong Jia Full-field three-dimensional measurement method
US7545512B2 (en) * 2006-01-26 2009-06-09 Koh Young Technology Inc. Method for automated measurement of three-dimensional shape of circuit boards
FR2902197B1 (fr) * 2006-06-07 2008-07-25 Univ Bretagne Occidentale Etab Procede de localisation d'un objet magnetique, produit programme d'ordinateur, moyen de stockage et dispositif de localisation correspondants
US20080117438A1 (en) * 2006-11-16 2008-05-22 Solvision Inc. System and method for object inspection using relief determination
JP4930834B2 (ja) 2007-01-10 2012-05-16 株式会社ニコン 形状測定方法
JP5016520B2 (ja) 2008-02-26 2012-09-05 パナソニック株式会社 3次元形状計測方法および装置
JP5055191B2 (ja) 2008-04-24 2012-10-24 パナソニック株式会社 3次元形状計測方法および装置
DE102010064593A1 (de) 2009-05-21 2015-07-30 Koh Young Technology Inc. Formmessgerät und -verfahren
US20110218428A1 (en) * 2010-03-04 2011-09-08 Medical Scan Technologies, Inc. System and Method for Three Dimensional Medical Imaging with Structured Light
US8391590B2 (en) * 2010-03-04 2013-03-05 Flashscan3D, Llc System and method for three-dimensional biometric data feature detection and recognition

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1437000A (zh) * 2002-02-09 2003-08-20 沈阳同联集团高新技术有限公司 投影栅线测量物体三维表面形状的方法和装置
JP2007322162A (ja) * 2006-05-30 2007-12-13 3D Media Co Ltd 3次元形状測定装置及び3次元形状測定方法

Also Published As

Publication number Publication date
KR20120054130A (ko) 2012-05-30
JP5546522B2 (ja) 2014-07-09
CN102538706A (zh) 2012-07-04
US20120127305A1 (en) 2012-05-24
US8982330B2 (en) 2015-03-17
KR101174676B1 (ko) 2012-08-17
JP2012112952A (ja) 2012-06-14

Similar Documents

Publication Publication Date Title
CN102538706B (zh) 绘制表面轮廓的方法和装置
Salvi et al. A state of the art in structured light patterns for surface profilometry
CN100520285C (zh) 投射多频光栅的物体表面三维轮廓的视觉测量方法
US8462207B2 (en) Depth ranging with Moiré patterns
EP1643210B1 (en) Method and apparatus for measuring shape of an object
TWI436052B (zh) 缺陷檢查裝置及缺陷檢查方法
US9568304B2 (en) Image sequence and evaluation method and system for structured illumination microscopy
JP2011064482A (ja) 高速三次元計測装置及び高速三次元計測方法
JPWO2006013635A1 (ja) 3次元形状計測方法及びその装置
CN103697815A (zh) 基于相位编码的混频结构光三维信息获取方法
JPWO2007088789A1 (ja) 表面形状の測定方法およびこれを用いた装置
US6876458B2 (en) Method and device for determining the absolute coordinates of an object
CN103940370A (zh) 基于周期互质混合编码的目标物体三维信息获取方法
CN103942802A (zh) 基于随机模板的结构光动态场景深度获取方法
CN205718873U (zh) 一种双频相移三维投影测量仪
JP6682251B2 (ja) 検査システム、及び検査方法
US20180112975A1 (en) Fringe Projection Method, Fringe Projection Device and Computer Program Product (as translated)
JP2010060420A (ja) 表面形状および/または膜厚測定方法およびその装置
KR101549310B1 (ko) 투명체의 결함검출을 위한 위상천이 영사장치 및 그 방법
US20150324657A1 (en) Method and apparatus for measuring 3d shape by using derivative moire
JP2010151781A (ja) 周波数推定方法、周波数推定装置、表面形状測定方法及び表面形状測定装置
JP6130805B2 (ja) 距離測定装置および方法
KR101538557B1 (ko) 위상차 이미지를 이용한 결함검출 장치 및 그 방법
Tan et al. Feature extraction from Moiré pattern images for tilt sensing
KR102677574B1 (ko) 측정시편의 이동상황에서 편향 측정법을 이용한 자유곡면 분석시스템 및 분석방법

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
CP01 Change in the name or title of a patent holder

Address after: Han Guoshouer

Patentee after: Gaoying Technology Co.,Ltd.

Address before: Han Guoshouer

Patentee before: KOH YOUNG TECHNOLOGY Inc.

CP01 Change in the name or title of a patent holder