CN102538706B - 绘制表面轮廓的方法和装置 - Google Patents
绘制表面轮廓的方法和装置 Download PDFInfo
- 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
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/50—Depth or shape recovery
- G06T7/521—Depth or shape recovery from laser ranging, e.g. using interferometry; from the projection of structured light
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01B—MEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
- G01B11/00—Measuring arrangements characterised by the use of optical techniques
- G01B11/24—Measuring arrangements characterised by the use of optical techniques for measuring contours or curvatures
- G01B11/25—Measuring 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/2518—Projection 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(φi+αbi)
在上述表达式中,下标‘b’表示多个方向中的一个方向(或通道),下标‘i’与光栅图案被移动(i-1)次时的第i个图像对应(例如,在N桶算法中,‘i’可以是从1到N的一个自然数),‘I0’是图案光的最大强度(或调幅),‘A’是N桶算法中图案图像的平均强度,而因移动光栅图案而形成的相位αbi可以由下述表达式2表示。
表达式2
使用上述表达式1,沿每个方向b的相位可以由下述表达式3表示。
表达式3
使用上述表达式3,在每个方向b上的高度Hb可以近似为下述表达式4。
表达式4
在上述表达式中,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
也就是说,通过对每个方向(b)上的高度Hb进行梯度操作,获得每个方向(b)上的矢量场(Gb(x,y))。
然后,获得关于每个方向(b)上的高度Hb的可靠性指数(Cb)(步骤S140)。可靠性指数(Cb)可以包括可见度指数、反射指数和阴影指数中的至少一者。
为了获得可见度指数,在4桶算法的情况下,表达式1中的最大强度(或调幅)I0可以由下述表达式6表示。
表达式6
此外,表达式1中的平均强度A可以由下述表达式7表示。
表达式7
然后,可见度VF可以由使用表达式6和表达式7获得的下述表达式8表示。
表达式8
根据本示例性实施例的可见度指数VIb可以由使用可见度VF获得的下述表达式9表示。
表达式9
VIb=(VFb>vthr)
在上述表达式中,vthr是可见度的阈值,并且该可见度阈值vthr可以通过实验来确定。表达式9的含义如下。即,例如,当多个方向(b)中的第一方向的可见度指数(VI1)大于阈值(vthr)时,将逻辑值真(1)赋予第一方向的可见度指数(VI1),而当第一方向的可见度指数(VI1)小于或等于阈值(vthr)时,将逻辑值假(0)赋予可见度指数(VI1)。
此外,与全反射引起的反射缺陷有关的多个方向(b)上的整合的反射指数(SpIb)可以由下述表达式10表示。
表达式10
在上述表达式中,反射阈值(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
在上述表达式中,阴影阈值(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
为了将表达式14最小化,可以使用表示为下述表达式15的泊松方程。
表达式15
如上述表达式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
为了解上述偏微分方程,将边界的高度设定为地。因此,应用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
当将上述表达式18所表示的傅里叶级数假设为由表达式17所表示的方程的解并应用于表达式17时,可以得到下述表达式19的方程。
表达式19
在上述表达式中,方程右侧的系数Fjk可以由下述表达式20表示。
表达式20
在上述表达式20中,可以通过将散度算子应用于表达式15的整合的矢量场并进行傅里叶展开来获得系数gmn。
可以使用正弦函数的可加性公式将表达式19表示为下述表达式21。
表达式21
因此,表达式18的系数Ujk可以由下述表达式22表示,从而可以获得表达式15的高度的解。
表达式22
根据本示例性实施例,在通过沿各个方向应用狭缝图案所获得的高度中,选择出没有诸如阴影缺陷、反射缺陷等缺陷的高度。并且,将没有缺陷的高度整合以提高精确度。
图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
在上述表达式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所述的绘制表面轮廓的装置,其中,所述计算部分包括:
高度计算单元,其依次使用所述图案图像计算与所述至少两个方向对应的所述高度;
可靠性指数计算单元,其计算与所述至少两个方向对应的高度的可靠性指数;
存储单元,其存储计算出的高度和计算出的可靠性指数;
矢量场计算单元,其计算所述矢量场;
整合的矢量场计算单元,其使用与所述至少两个方向对应的所述可靠性指数和所述矢量场计算所述整合的矢量场;
边界设定单元,其设定所述目标物体的边界;以及
泊松方程求解单元,其通过对由所述边界设定单元设定的所述边界应用边界条件来求解所述泊松方程,所述泊松方程的右侧通过对所述整合的矢量场应用散度算子而获得,以及所述泊松方程的左侧通过对所述高度应用拉普拉斯算子而获得。
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)
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)
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)
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 |
-
2010
- 2010-11-19 KR KR1020100115351A patent/KR101174676B1/ko active IP Right Grant
-
2011
- 2011-11-17 US US13/298,969 patent/US8982330B2/en active Active
- 2011-11-21 CN CN201110376531.7A patent/CN102538706B/zh active Active
- 2011-11-21 JP JP2011253957A patent/JP5546522B2/ja active Active
Patent Citations (2)
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 |