CN1912540A - 微纳结构3d轮廓测量中基于坐标变换的倾斜误差补偿方法 - Google Patents

微纳结构3d轮廓测量中基于坐标变换的倾斜误差补偿方法 Download PDF

Info

Publication number
CN1912540A
CN1912540A CNA2006100154965A CN200610015496A CN1912540A CN 1912540 A CN1912540 A CN 1912540A CN A2006100154965 A CNA2006100154965 A CN A2006100154965A CN 200610015496 A CN200610015496 A CN 200610015496A CN 1912540 A CN1912540 A CN 1912540A
Authority
CN
China
Prior art keywords
micro
nano structure
profile
tilt error
reference plane
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
CNA2006100154965A
Other languages
English (en)
Other versions
CN100402978C (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.)
Tianjin University
Original Assignee
Tianjin 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 Tianjin University filed Critical Tianjin University
Priority to CNB2006100154965A priority Critical patent/CN100402978C/zh
Publication of CN1912540A publication Critical patent/CN1912540A/zh
Application granted granted Critical
Publication of CN100402978C publication Critical patent/CN100402978C/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Length Measuring Devices By Optical Means (AREA)

Abstract

一种微纳结构3D轮廓测量中基于坐标变换的倾斜误差补偿方法,是由如下过程完成的:1.采集被测物表面干涉图,得到被测物表面初始轮廓图;2.确定被测物的基准平面;3.进行坐标系旋转,使新坐标系的XOY平面与基准平面平行,重构被测物表面3D轮廓。确定被测物的基准平面,是使用最小二乘法确定调平表面轮廓的基准平面,包括有最小二乘中线法和面内最小二乘法。本发明计算量小,使用简单,处理效率高;最大限度地避免了轮廓的变形;不仅可用于相移显微干涉***的测量应用中,也可以用于其他测量MEMS/NEMS表面3D轮廓的技术的图像处理中;可以根据不同应用的需求,选择不同的面作为基准,得到不同视角的表面三维轮廓信息,运用灵活。

Description

微纳结构3D轮廓测量中基于坐标变换的倾斜误差补偿方法
技术领域
本发明涉及一种微/纳测试技术。特别是涉及一种微纳结构3D轮廓测量中基于坐标变换的倾斜误差补偿方法。
背景技术
微/纳机电***(Micro/Nano ElectroMechanical System,MEMS/NEMS)是集传感、信息处理和执行于一体的集成微/纳***,近年来已经成为高新技术领域的研究热点。MEMS/NEMS的特征尺寸多为微米和纳米量级,具有体积小、成本低、可批量加工和可靠性高等优点,这使得MEMS/NEMS在传感器、执行器、射频、光学、军事、航空航天和生物医学等领域都具有广泛的应用前景。但是由于分析和测试MEMS/NEMS可靠性与结构特性参数技术的匮乏,其开发和商业化应用受到了很大制约。因此,对MEMS/NEMS的可靠性和结构特性参数进行分析和测试的研究越来越受到人们的重视。
在对MEMS/NEMS器件的可靠性和结构特性参数进行测量时,常常要求获取微纳结构的表面信息,如三维轮廓、表面粗糙度、微小的位移和变形等。测量微纳结构表面3D轮廓的方法有很多,相移显微干涉技术凭借其高精度、高垂直分辨率、测量简单快捷以及对器件无损伤等优点,成为该类精密测量中最常用的手段之一。这种方法通过获取被测微纳结构表面的干涉图像并对其进行相位提取和相位展开从而获得被测微纳结构表面的3D轮廓,垂直分辨率可以达到纳米量级。
然而,在相移显微干涉测量技术中,未经调平处理而直接得到的微纳结构的表面3D轮廓信息往往存在着一定的倾斜误差,这种倾斜误差主要是由于放置待测微纳结构的载物台通常具有一定的倾斜度造成的。虽然载物台的这种倾斜度可能比较微小,但是在微纳结构表面3D轮廓精密测量中,这种倾斜度对某些应用领域的影响可能非常大。例如,在某些微纳结构力学特性测量中,如果需要得到微纳结构表面相对于基底或某个平面的高度差,带有倾斜误差的结果与真实结果之间的相差就可能非常悬殊,从而导致测试结果的严重失真。因此,在微纳结构表面3D轮廓测量中,对载物台的倾斜误差进行补偿是非常必要的,然而传统的倾斜误差补偿方法往往是补偿精度不高,或容易造成表面3D轮廓的变形,或补偿算法非常复杂从而计算时间很长。
发明内容
本发明所要解决的技术问题是,提供一种在降低补偿算法复杂性的同时有效地提高了倾斜误差的补偿精度,最大限度地避免了表面3D轮廓的变形的微纳结构3D轮廓测量中基于坐标变换的倾斜误差补偿方法。
本发明所采用的技术方案是:一种微纳结构3D轮廓测量中基于坐标变换的倾斜误差补偿方法,是包括以下步骤:1.采集被测物微纳结构表面的干涉图,计算被测微纳结构表面的初始轮廓图;2.确定被测物的基准平面;3.进行坐标系旋转,使新坐标系的XOY平面与基准平面平行,重构被测微纳结构的表面3D轮廓。
所述的采集被测物微纳结构表面的干涉图,计算被测微纳结构表面的初始轮廓图,是通过相位提取和相位展开实现的。
所述的确定被测微纳结构的基准平面,是通过最小二乘中线法和面内最小二乘法实现的。
所述的进行坐标系旋转是基于欧式变换的坐标旋转法实现的,即先绕Y轴旋转初始表面轮廓在X方向的倾斜角度,补偿X方向的倾斜误差,然后绕X轴旋转初始表面轮廓在Y方向的倾斜角度,补偿Y方向的倾斜误差。
所述的重构被测微纳结构的表面3D轮廓是通过重构被测表面相对于其准平面的表面3D轮廓来实现的。
本发明的微纳结构3D轮廓测量中基于坐标变换的倾斜误差补偿方法,具有如下特点:
1、所基于的理论成熟,计算量小,使用简单,处理效率高;
2、使用旋转坐标系的方法重构三维信息,最大限度地避免了表面3D轮廓的变形,克服了传统调平方法的缺陷;
3、不仅可用于相移显微干涉***的测量应用中,也可以用于其他测量MEMS/NEMS表面3D轮廓的技术的图像处理中,嵌入简单;
4、可以根据不同应用的需求,选择不同的面作为基准,从而得到不同视角的表面三D轮廓信息,运用灵活;
5、该补偿方法能得到被测表面各象素点在旋转后的新坐标系中的X、Y、Z坐标值,在使用的显微物镜放大倍率很高或者被测表面面积很小且起伏较大时更能显示其高精度的特点。
附图说明
图1是本发明所使用的相移显微干涉法的原理图;
图2是本发明的旋转坐标系调平被测表面图像的示意图。
其中:
1:光源              2:干涉图
3:分光镜            4:参考镜
5:被测表面          6:被测表面
7:基准平面          8:原始坐标系
9:旋转后的新坐标系
具体实施方式
下面结合附图和实施例对本发明的微纳结构3D轮廓测量中基于坐标变换的倾斜误差补偿方法做出详细说明。
本发明的微纳结构3D轮廓测量中基于坐标变换的倾斜误差补偿方法,是在采用相移显微干涉测量微结构表面3D轮廓的图像处理过程中。通过发下步骤实现的:
1.采集被测物表面干涉图,计算被测微纳结构表面初始轮廓图。
首先得用如图1所示的相移显微干涉***,得到被测表面干涉图像2,采集一系列干涉图像并保存;然后对干涉图像2进行相位提取和相位展开得到被测微纳结构表面的初始轮廓图像,此时的表面3D轮廓数据包含有倾斜误差,如图2中的被测表面6。
2.确定被测物的基准平面。
使用最小二乘法确定3D表面轮廓的基准平面。所述的最小二乘法包括有最小二乘中线法和面内最小二乘法。
最小二乘中线法,要求在测量取样长度内,物体表面轮廓曲线上各点的轮廓偏距(轮廓曲线上各点到基准线的距离)的平方和最小。根据这个要求可以确定沿X或Y方向逐条基准线的斜率。将X方向逐条基准线的斜率求和取平均值作为X方向的倾斜斜率,将Y方向逐条基准线的斜率求和取平均作为Y方向的倾斜斜率。假设基准平面经过原点,综合X、Y两方向的斜率,即可得到基准平面方程ax+by+c+z=0;
面内最小二乘法:是在一个选定的区域内使用最小二乘法求基准平面,它要求所求得的基准平面各像素点与所选区域对应像素点的Z向高度差的平方和最小,从而确定基准平面方程ax+by+c+z=0中的系数a、b、c。
3.进行坐标系旋转,使新坐标系的XOY平面与基准平面平行。
在基准平面确定之后,需要得到被测表面相对于这个基准面的高度从而重构被测表面的三维轮廓。本发明采用基于欧式变换的坐标旋转法来实现这一过程。也就是:在确定基准平面之后,根据基准平面与原始坐标系XOY平面之间的夹角(即初始表面轮廓的倾斜角度)旋转直角坐标系使得新坐标系的XOY平面与基准平面平行,从而达到调平被测表面并补偿倾斜误差的目的。
根据求得的基准平面方程可以得到其与原始坐标系8的X轴和Y轴之间的夹角,依据欧式变换的原理,将原始坐标系8绕Y轴旋转基准面2与原始坐标系8的X轴间夹角的角度大小,便可以补偿原始坐标系8在X方向的倾斜;同理,再将得到的坐标系绕X轴旋转基准面与原始坐标系8的Y轴间夹角的角度大小,便可以补偿原始坐标系8在Y方向的倾斜。
设εX、εY分别为原始坐标系8与基准平面2在X和Y方向的夹角,X1Y1Z1是旋转之前的坐标值,X2Y2Z2为旋转之后的坐标值,则坐标旋转的公式为:
X 2 Y 2 Z 2 = R 1 ( ϵ X ) R 2 ( ϵ Y ) X 1 Y 1 Z 1 - - - ( 1 )
其中,
R 1 ( ϵ X ) = 1 0 0 0 cos ϵ X sin ϵ X 0 - sin ϵ X cos ϵ X - - - ( 2 )
R 2 ( ϵ Y ) = cos ϵ Y 0 - sin ϵ Y 0 1 0 sin ϵ Y 0 cos ϵ Y - - - ( 3 )
4、重构被测微纳结构的表面3D轮廓。
如图2所示,旋转之后的新坐标系9的XOY平面与基准平面2平行,相应地Z向的高度信息也就成为了相对于所选基准面的3D轮廓数据。即获得到被测微纳结构的表面6在新坐标系9中相对于基准面7的3D轮廓数据,重构了被测微纳结构的表面相对于基准平面的3D轮廓。根据应用的不同需求,应用此方法可以得到针对于不同基准平面的表面3D轮廓,相当于得到不同视角的三D轮廓信息而最大可能地避免表面3D轮廓信息的变形。
本发明公开和揭示的方法可借鉴本文公开内容,尽管本发明的方法已通过较佳实施例进行了描述,但是本领域技术人员明显能在不脱离本发明内容、精神和范围内对本文所述的方法进行拼接或改动,或增减某些环节的,更具体地说,所有相类似的替换和改动对本领域技术人员来说是显而易见的,它们都被视为包括在本发明精神、范围和内容中。

Claims (5)

1.一种微纳结构3D轮廓测量中基于坐标变换的倾斜误差补偿方法,其特征在于,是包括以下步骤:1.采集被测物微纳结构表面的干涉图,计算被测微纳结构表面的初始轮廓图;2.确定被测物的基准平面;3.进行坐标系旋转,使新坐标系的XOY平面与基准平面平行,重构被测微纳结构的表面3D轮廓。
2.根据权利要求1所述的微纳结构3D轮廓测量中基于坐标变换的倾斜误差补偿方法,其特征在于,所述的采集被测物微纳结构表面的干涉图,计算被测微纳结构表面的初始轮廓图,是通过相位提取和相位展开实现的。
3.根据权利要求1所述的微纳结构3D轮廓测量中基于坐标变换的倾斜误差补偿方法,其特征在于,所述的确定被测微纳结构的基准平面,是通过最小二乘中线法和面内最小二乘法实现的。
4.根据权利要求1所述的微纳结构3D轮廓测量中基于坐标变换的倾斜误差补偿方法,其特征在于,所述的进行坐标系旋转是基于欧式变换的坐标旋转法实现的,即先绕Y轴旋转初始表面轮廓在X方向的倾斜角度,补偿X方向的倾斜误差,然后绕X轴旋转初始表面轮廓在Y方向的倾斜角度,补偿Y方向的倾斜误差。
5.根据权利要求1所述的微纳结构3D轮廓测量中基于坐标变换的倾斜误差补偿方法,其特征在于,所述的重构被测微纳结构的表面3D轮廓是通过重构被测表面相对于其准平面的表面3D轮廓来实现的。
CNB2006100154965A 2006-08-30 2006-08-30 微纳结构3d轮廓测量中基于坐标变换的倾斜误差补偿方法 Expired - Fee Related CN100402978C (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CNB2006100154965A CN100402978C (zh) 2006-08-30 2006-08-30 微纳结构3d轮廓测量中基于坐标变换的倾斜误差补偿方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CNB2006100154965A CN100402978C (zh) 2006-08-30 2006-08-30 微纳结构3d轮廓测量中基于坐标变换的倾斜误差补偿方法

Publications (2)

Publication Number Publication Date
CN1912540A true CN1912540A (zh) 2007-02-14
CN100402978C CN100402978C (zh) 2008-07-16

Family

ID=37721541

Family Applications (1)

Application Number Title Priority Date Filing Date
CNB2006100154965A Expired - Fee Related CN100402978C (zh) 2006-08-30 2006-08-30 微纳结构3d轮廓测量中基于坐标变换的倾斜误差补偿方法

Country Status (1)

Country Link
CN (1) CN100402978C (zh)

Cited By (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101806585A (zh) * 2010-04-09 2010-08-18 中北大学 基于红外光干涉技术的mems器件形貌测量方法
CN102645202A (zh) * 2012-05-11 2012-08-22 厦门大学 大口径非球面工件轮廓的测量方法
CN102749041A (zh) * 2012-07-13 2012-10-24 合肥工业大学 一种螺旋桨型面轮廓误差测量仪及测量方法
CN102922066A (zh) * 2012-09-21 2013-02-13 大连理工大学 一种平板间微纳液膜厚度测试调平方法和装置
CN102999904A (zh) * 2012-11-14 2013-03-27 北京理工大学 基于平面的微纳物体图像倾斜校正方法
CN103077513A (zh) * 2012-11-16 2013-05-01 北京理工大学 基于二次曲面的微纳物体图像倾斜校正方法
CN106643602A (zh) * 2016-12-24 2017-05-10 上海与德信息技术有限公司 一种探针基准补偿***及方法
CN110006364A (zh) * 2019-03-18 2019-07-12 南京师范大学 基于圆条纹径向空间载波相移的三维实时显微测量方法
CN110455246A (zh) * 2019-08-27 2019-11-15 浙江大学 一种用于共形光学元件的面形测量装置及方法
CN113050660A (zh) * 2021-06-01 2021-06-29 新石器慧通(北京)科技有限公司 误差补偿方法、装置、计算机设备及存储介质

Family Cites Families (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP3532347B2 (ja) * 1996-05-24 2004-05-31 株式会社リコー 形状測定装置
JP3474448B2 (ja) * 1998-09-01 2003-12-08 株式会社リコー 座標軸直角度誤差の校正方法及び三次元形状測定装置
CN2784875Y (zh) * 2005-03-22 2006-05-31 范光照 一种基于光纤传像束的新型微三维形貌测量装置

Cited By (17)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101806585B (zh) * 2010-04-09 2011-06-22 中北大学 基于红外光干涉技术的mems器件形貌测量方法
CN101806585A (zh) * 2010-04-09 2010-08-18 中北大学 基于红外光干涉技术的mems器件形貌测量方法
CN102645202B (zh) * 2012-05-11 2014-09-10 厦门大学 大口径非球面工件轮廓的测量方法
CN102645202A (zh) * 2012-05-11 2012-08-22 厦门大学 大口径非球面工件轮廓的测量方法
CN102749041A (zh) * 2012-07-13 2012-10-24 合肥工业大学 一种螺旋桨型面轮廓误差测量仪及测量方法
CN102749041B (zh) * 2012-07-13 2014-10-08 合肥工业大学 一种螺旋桨型面轮廓误差测量仪及测量方法
CN102922066A (zh) * 2012-09-21 2013-02-13 大连理工大学 一种平板间微纳液膜厚度测试调平方法和装置
CN102922066B (zh) * 2012-09-21 2014-11-12 大连理工大学 一种平板间微纳液膜厚度测试调平方法和装置
CN102999904A (zh) * 2012-11-14 2013-03-27 北京理工大学 基于平面的微纳物体图像倾斜校正方法
CN102999904B (zh) * 2012-11-14 2015-06-10 北京理工大学 基于平面的微纳物体图像倾斜校正方法
CN103077513A (zh) * 2012-11-16 2013-05-01 北京理工大学 基于二次曲面的微纳物体图像倾斜校正方法
CN106643602A (zh) * 2016-12-24 2017-05-10 上海与德信息技术有限公司 一种探针基准补偿***及方法
CN106643602B (zh) * 2016-12-24 2019-03-12 重庆与德通讯技术有限公司 一种探针基准补偿***及方法
CN110006364A (zh) * 2019-03-18 2019-07-12 南京师范大学 基于圆条纹径向空间载波相移的三维实时显微测量方法
CN110455246A (zh) * 2019-08-27 2019-11-15 浙江大学 一种用于共形光学元件的面形测量装置及方法
CN110455246B (zh) * 2019-08-27 2020-11-03 浙江大学 一种用于共形光学元件的面形测量装置及方法
CN113050660A (zh) * 2021-06-01 2021-06-29 新石器慧通(北京)科技有限公司 误差补偿方法、装置、计算机设备及存储介质

Also Published As

Publication number Publication date
CN100402978C (zh) 2008-07-16

Similar Documents

Publication Publication Date Title
CN1912540A (zh) 微纳结构3d轮廓测量中基于坐标变换的倾斜误差补偿方法
US8160351B2 (en) Method and apparatus for mura detection and metrology
CN111207657B (zh) 一种大尺寸圆柱类零件素线直线度误差检测装置及方法
CN1180222C (zh) 双频共焦台阶高度显微测量装置
CN1851522A (zh) 激光扫描共聚焦显微镜扫描畸变现象的全场校正方法
Cai et al. An ultra-precision tool nanoindentation instrument for replication of single point diamond tool cutting edges
CN1236278C (zh) 光轴与安装基面平行度的检验方法
CN102507985A (zh) 用于测量具有垂直侧壁微结构的扫描隧道显微镜倾斜测量方法
CN211178284U (zh) 一种大尺寸圆柱类零件素线直线度误差检测装置
CN101046374A (zh) 基于原子力显微镜力曲线的纳米梁厚度的非破坏测量方法
KR101016229B1 (ko) 초정밀 리니어 스테이지의 운동오차 측정방법 및 측정시스템
CN1892178A (zh) 非接触式测量法
CN101819217B (zh) 一种微纳米尺度平面周期性结构的反演方法
Kang et al. New measurement method of Poisson’s ratio of thin films by applying digital image correlation technique
CN101033949A (zh) 一种基于错位相关原理的物体应变测量方法及装置
CN101046371A (zh) 基于显微干涉和有限差分的微/纳悬臂梁的曲率测量方法
Heikkinen et al. Step height standards based on self-assembly for 3D metrology of biological samples
Costa et al. Residual stress measurement in PVD optical coatings by microtopography
Xu et al. Surface form measurement and analysis of a cylindrical workpiece with microstructures
CN1932438A (zh) 对称变形场的单光束电子散斑干涉二维检测方法
CN109544639B (zh) 一种多镜面单相机三维振动测试装置及方法
CN2826363Y (zh) 双波长二维空间相移电子散斑干涉仪
Zheng et al. A non-contact swing-arm profilometer with the spectrally-resolved-interferometry distance sensor
Hu et al. An on-axis self-calibration approach for precision rotary metrology stages based on an angular artifact plate
Huntley et al. Large-scale full-field metrology using projected fringes: some challenges and solutions

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
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20080716

Termination date: 20210830