CN105403170A - 一种显微3d形貌测量方法及装置 - Google Patents
一种显微3d形貌测量方法及装置 Download PDFInfo
- Publication number
- CN105403170A CN105403170A CN201510922156.XA CN201510922156A CN105403170A CN 105403170 A CN105403170 A CN 105403170A CN 201510922156 A CN201510922156 A CN 201510922156A CN 105403170 A CN105403170 A CN 105403170A
- Authority
- CN
- China
- Prior art keywords
- image
- micro
- objective table
- measurement method
- topography measurement
- 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
Abstract
本发明涉及一种显微3D形貌测量方法,步骤如下:1)分别在焦前与焦后拍摄载物台上的待测物体的图像If=If(i,j)、Ib=Ib(i,j),0≤i<I,0≤j<J;其中,I是图像总行数,J是图像总列数;2)根据如下公式计算得到待测物体在每一个位置(i,j)的待测图像I(i,j),公式为:
Description
技术领域
本发明涉及数字显微领域,更具体地说,涉及一种显微3D形貌测量方法,以及一种显微3D形貌测量装置。
背景技术
在先进数字显微镜以及现今追求高精度分辨率观察分析仪器中,有的具有高达微米甚至纳米的横向分辨率,然而纵向分辨率常常局限在几个微米,纵向分辨率难以提高。
专门为了实现光学切片、提高纵向分辨率的激光扫描共聚焦显微镜,在纵向上的分辨率也难以达到亚微米,而且3D成像速度非常缓慢。需要在纵向上不同的Z位置获取多层光学切片图像,之后来重建,获取每一个水平位置XY的表面高度Z值。由于需要在X方向、Y方向,Z方向3个方向激光扫描,因此该方法非常耗时。
另外在电动数字显微镜成像过程中,实时判断目标表面距离物镜成像面的距离也是难以解决的问题。
发明内容
本发明的目的在于克服现有技术的不足,提供一种用于数字显微镜或者放大倍数大于1的数字成像系统中快速获得3D重建或者3D形貌的显微3D形貌测量方法,以及一种显微3D形貌测量装置。
本发明的技术方案如下:
一种显微3D形貌测量方法,步骤如下:
1)分别在焦前与焦后拍摄载物台上的待测物体的图像If=If(i,j)、Ib=Ib(i,j),0≤i<I,0≤j<J;其中,I是图像总行数,J是图像总列数;
2)根据如下公式计算得到待测物体在每一个位置(i,j)的待测图像I(i,j),公式为:
3)将待测图像I(i,j)与预先刻度的I(h)曲线对比,确定待测物体在每一个位置(i,j)表面形貌相对高度h(i,j),0≤i<I,0≤j<J。
作为优选,步骤1)中,图像If=If(i,j)的获取步骤为:通过移动显微镜的物镜或者移动显微镜的载物台,使得载物台表面在显微镜物镜焦面前方υf;将摄像头置于显微成像系统的成像面上,拍摄一副放置在载物台上的待测物体的图像If=If(i,j),0≤i<I,0≤j<J;其中,I是图像总行数,J是图像总列数。
作为优选,步骤1)中,图像Ib=Ib(i,j)的获取步骤为:通过移动显微镜的物镜或者移动显微镜的载物台,使得载物台表面在显微镜物镜焦面后方υb;将摄像头置于显微成像系统的成像面上,拍摄一副放置在载物台上的待测物体的图像Ib=Ib(i,j),0≤i<I,0≤j<J;其中,I是图像总行数,J是图像总列数。
作为优选,步骤1)中,图像If=If(i,j)与Ib=Ib(i,j)的失焦距离相同。
作为优选,步骤3)中,I(h)曲线的获取步骤为:依次在载物台上放置若干不同高度的全反射平面的物体,对每个物体分别在焦前与焦后拍摄图像IF(i,j)、IB(i,j),对每一点(i,j)根据如下公式计算出参考图像I(h,i,j),绘制得到相应的I(h)曲线,公式为:
其中,0≤i<I,0≤j<J。
作为优选,图像IF(i,j)的获取步骤为:将摄像头放置在显微成像系统的成像面的位置,通过移动显微镜的物镜或者移动显微镜的载物台,使得载物台表面在显微镜物镜焦面前方υF;依次放置N≥2个不同高度hi=k*δ的具有全反射平面的物体于载物台,对每一个物体拍摄一张图像IF(i,j),其中;i≤N,
作为优选,图像IB(i,j)的获取步骤为:将摄像头放置在显微成像系统的成像面的位置,通过移动显微镜的物镜或者移动显微镜的载物台,使得载物台表面在显微镜物镜焦面后方υB;依次放置N≥2个不同高度hi=k*δ的具有全反射平面的物体于载物台,对每一个物体拍摄一张图像IB(i,j),其中;i≤N,
作为优选,图像IF(i,j)、IB(i,j)的失焦距离相同。
一种显微3D形貌测量装置,运行所述的显微3D形貌测量方法,载物台通过电机移动至特设的位置,通过所述的显微3D形貌测量方法获取3D形貌或者进行3D重建。
本发明的有益效果如下:
本发明所述的方法实现的纵向形貌测量效率高,仅仅获取两幅图像就能够以非接触的方式实现多点测量待测物的表面高度,不需要采用传统方法进行层层扫描来获取多幅图像。虽然获得刻度曲线时候需要多个图像,但是一台仪器只需要刻度一次,即可一直用于后续使用,使用不同的标本也不需要重新刻度。
本发明所述的方法实现的纵向形貌测量精度高,高于显微镜物镜的纵向光学分辨率。
本发明所述的设备结构简单:对于传统的数字显微镜,仅仅需要添加一个电机如步进电机使得载物台或者物镜可以纵向移动到指定距离,即可通过所述的显微3D形貌测量方法获取3D形貌或者进行3D重建。
具体实施方式
本发明为了解决现有技术存在的不足,提供一种显微3D形貌测量方法,以及一种显微3D形貌测量装置。本发明的核心技术方案是:通过实时拍摄少量图像,再通过与预先标定的曲线对比,计算出目标平面上每一点离开物镜的Z距离。
所述的显微3D形貌测量方法,步骤如下:
步骤1)分别在焦前与焦后拍摄载物台上的待测物体的图像If=If(i,j)、Ib=Ib(i,j),0≤i<I,0≤j<J;其中,I是图像总行数,J是图像总列数。
具体地,图像If=If(i,j)的获取步骤为:通过移动显微镜的物镜或者移动显微镜的载物台,使得载物台表面在显微镜物镜焦面前方υf;将摄像头置于显微成像系统的成像面上,拍摄一幅放置在载物台上的待测物体的图像If=If(i,j),0≤i<I,0≤j<J;其中,I是图像总行数,J是图像总列数。
图像Ib=Ib(i,j)的获取步骤为:通过移动显微镜的物镜或者移动显微镜的载物台,使得载物台表面在显微镜物镜焦面后方υb;将摄像头置于显微成像系统的成像面上,拍摄一幅放置在载物台上的待测物体的图像Ib=Ib(i,j),0≤i<I,0≤j<J;其中,I是图像总行数,J是图像总列数。
图像If=If(i,j)与Ib=Ib(i,j)的失焦距离相同,即载物台表面在显微镜物镜焦面前方的距离υf与载物台表面在显微镜物镜焦面后方的距离υb大小相等,方向相反。
步骤2)根据如下公式计算得到待测物体在每一个位置(i,j)的待测图像I(i,j),公式为:
步骤3)将待测图像I(i,j)与预先刻度的I(h)曲线对比,确定待测物体在每一个位置(i,j)表面形貌相对高度h(i,j),0≤i<I,0≤j<J。
I(h)曲线的获取步骤为:依次在载物台上放置若干不同高度的全反射平面的物体,对每个物体分别在焦前与焦后拍摄图像IF(i,j)、IB(i,j),对每一点(i,j)根据如下公式计算出参考图像I(h,i,j),绘制得到相应的I(h)曲线,公式为:
其中,0≤i<I,0≤j<J。
具体地,图像IF(i,j)的获取步骤为:将摄像头放置在显微成像系统的成像面的位置,通过移动显微镜的物镜或者移动显微镜的载物台,使得载物台表面在显微镜物镜焦面前方υF;依次放置N≥2个不同高度hi=k*δ的具有全反射平面的物体于载物台,对每一个物体拍摄一张图像IF(i,j),其中;i≤N,
图像IB(i,j)的获取步骤为:将摄像头放置在显微成像系统的成像面的位置,通过移动显微镜的物镜或者移动显微镜的载物台,使得载物台表面在显微镜物镜焦面后方υB;依次放置N≥2个不同高度hi=k*δ的具有全反射平面的物体于载物台,对每一个物体拍摄一张图像IB(i,j),其中;i≤N,
图像IF(i,j)、IB(i,j)的失焦距离相同,即载物台表面在显微镜物镜焦面前方的距离υF与载物台表面在显微镜物镜焦面后方的距离υB大小相等,方向相反。
基于所述的显微3D形貌测量方法,本发明还提供一种显微3D形貌测量装置,用于实现上述的显微3D形貌测量方法,载物台通过电机移动至特设的位置,通过所述的显微3D形貌测量方法获取3D形貌或者进行3D重建。
上述实施例仅是用来说明本发明,而并非用作对本发明的限定。只要是依据本发明的技术实质,对上述实施例进行变化、变型等都将落在本发明的权利要求的范围内。
Claims (9)
1.一种显微3D形貌测量方法,其特征在于,步骤如下:
1)分别在焦前与焦后拍摄载物台上的待测物体的图像If=If(i,j)、Ib=Ib(i,j),0≤i<1,0≤j<J;其中,I是图像总行数,J是图像总列数;
2)根据如下公式计算得到待测物体在每一个位置(i,j)的待测图像I(i,j),公式为:
3)将待测图像I(i,j)与预先刻度的I(h)曲线对比,确定待测物体在每一个位置(i,j)表面形貌相对高度h(i,j),0≤i<I,0≤j<J。
2.根据权利要求1所述的显微3D形貌测量方法,其特征在于,步骤1)中,图像If=If(i,j)的获取步骤为:通过移动显微镜的物镜或者移动显微镜的载物台,使得载物台表面在显微镜物镜焦面前方υf;将摄像头置于显微成像系统的成像面上,拍摄一幅放置在载物台上的待测物体的图像If=If(i,j),0≤i<I,0≤j<J;其中,I是图像总行数,J是图像总列数。
3.根据权利要求1所述的显微3D形貌测量方法,其特征在于,步骤1)中,图像Ib=Ib(i,j)的获取步骤为:通过移动显微镜的物镜或者移动显微镜的载物台,使得载物台表面在显微镜物镜焦面后方υb;将摄像头置于显微成像系统的成像面上,拍摄一幅放置在载物台上的待测物体的图像Ib=Ib(i,j),0≤i<I,0≤j<J;其中,I是图像总行数,J是图像总列数。
4.根据权利要求1所述的显微3D形貌测量方法,其特征在于,步骤1)中,图像If=If(i,j)与Ib=Ib(i,j)的失焦距离相同。
5.根据权利要求1所述的显微3D形貌测量方法,其特征在于,步骤3)中,I(h)曲线的获取步骤为:依次在载物台上放置若干不同高度的全反射平面的物体,对每个物体分别在焦前与焦后拍摄图像IF(i,j)、IB(i,j),对每一点(i,j)根据如下公式计算出参考图像I(h,i,j),绘制得到相应的I(h)曲线,公式为:
其中,0≤i<I,0≤j<J。
6.根据权利要求5所述的显微3D形貌测量方法,其特征在于,图像IF(i,j)的获取步骤为:将摄像头放置在显微成像系统的成像面的位置,通过移动显微镜的物镜或者移动显微镜的载物台,使得载物台表面在显微镜物镜焦面前方υF;依次放置N≥2个不同高度hi=k*δ的具有全反射平面的物体于载物台,对每一个物体拍摄一张图像IF(i,j),其中;i≤N,
7.根据权利要求5所述的显微3D形貌测量方法,其特征在于,图像IB(i,j)的获取步骤为:将摄像头放置在显微成像系统的成像面的位置,通过移动显微镜的物镜或者移动显微镜的载物台,使得载物台表面在显微镜物镜焦面后方υB;依次放置N≥2个不同高度hi=k*δ的具有全反射平面的物体于载物台,对每一个物体拍摄一张图像IB(i,j),其中;i≤N,
8.根据权利要求5所述的显微3D形貌测量方法,其特征在于,图像IF(i,j)、IB(i,j)的失焦距离相同。
9.一种显微3D形貌测量装置,其特征在于,运行权利要求1至8任一项所述的显微3D形貌测量方法,载物台通过电机移动至特设的位置,通过所述的显微3D形貌测量方法获取3D形貌或者进行3D重建。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510922156.XA CN105403170B (zh) | 2015-12-11 | 2015-12-11 | 一种显微3d形貌测量方法及装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510922156.XA CN105403170B (zh) | 2015-12-11 | 2015-12-11 | 一种显微3d形貌测量方法及装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105403170A true CN105403170A (zh) | 2016-03-16 |
CN105403170B CN105403170B (zh) | 2018-05-25 |
Family
ID=55468809
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510922156.XA Active CN105403170B (zh) | 2015-12-11 | 2015-12-11 | 一种显微3d形貌测量方法及装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105403170B (zh) |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2005090905A1 (en) * | 2004-03-16 | 2005-09-29 | Spiral Scratch Limited | Optical profilometer apparatus and method |
DE102006013316A1 (de) * | 2005-03-21 | 2006-09-28 | Daimlerchrysler Ag | 3D-Rekonstruktion von statischen Szenen durch vereinheitlichte Kombination der Methode des Depth from Defokus und der des Shape from Motion |
CN102188290A (zh) * | 2010-03-19 | 2011-09-21 | 卡尔斯特里姆保健公司 | 利用远心散焦的3d成像 |
JP2011215120A (ja) * | 2010-03-31 | 2011-10-27 | Nexco-Engineering Hokkaido Co Ltd | 浮遊する粒子の位置測定法 |
WO2012039636A2 (en) * | 2010-09-24 | 2012-03-29 | Christian Soeller | "3d localisation microscopy and 4d localisation microscopy and tracking methods and systems" |
CN102692364A (zh) * | 2012-06-25 | 2012-09-26 | 上海理工大学 | 一种基于模糊图像处理的动态颗粒测量装置及方法 |
-
2015
- 2015-12-11 CN CN201510922156.XA patent/CN105403170B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2005090905A1 (en) * | 2004-03-16 | 2005-09-29 | Spiral Scratch Limited | Optical profilometer apparatus and method |
DE102006013316A1 (de) * | 2005-03-21 | 2006-09-28 | Daimlerchrysler Ag | 3D-Rekonstruktion von statischen Szenen durch vereinheitlichte Kombination der Methode des Depth from Defokus und der des Shape from Motion |
CN102188290A (zh) * | 2010-03-19 | 2011-09-21 | 卡尔斯特里姆保健公司 | 利用远心散焦的3d成像 |
JP2011215120A (ja) * | 2010-03-31 | 2011-10-27 | Nexco-Engineering Hokkaido Co Ltd | 浮遊する粒子の位置測定法 |
WO2012039636A2 (en) * | 2010-09-24 | 2012-03-29 | Christian Soeller | "3d localisation microscopy and 4d localisation microscopy and tracking methods and systems" |
CN102692364A (zh) * | 2012-06-25 | 2012-09-26 | 上海理工大学 | 一种基于模糊图像处理的动态颗粒测量装置及方法 |
Non-Patent Citations (3)
Title |
---|
J RAYALA ET AL: "Estimation of depth from defocus as polynomial system identification", 《IEE PROCEEDINGS - VISION, IMAGE AND SIGNAL PROCESSING》 * |
孙玲 等: "一种基于矩不变原理的散焦测距方法", 《弹箭与制导学报》 * |
王鑫 等: "与纹理无关的散焦测距方法的研究", 《应用科技》 * |
Also Published As
Publication number | Publication date |
---|---|
CN105403170B (zh) | 2018-05-25 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102636496B (zh) | 光学表面疵病暗场检测中疵病宽度标定标准化方法 | |
CN102156037B (zh) | 光学检测装置及检测方法 | |
CN109751964B (zh) | 一种高精度非接触式管径测量方法及装置 | |
CN109061229B (zh) | 一种光场Micro-PIV系统的标定方法 | |
CN101526336A (zh) | 基于量块的线结构光三维视觉传感器标定方法 | |
CN111598931B (zh) | 一种单目视觉系统成像参数标定装置及方法 | |
CN113568153B (zh) | 一种显微成像设备以及纳米尺度三维形貌测量系统 | |
CN101852677A (zh) | 一种提高长焦距平行光管焦距检测精度的方法 | |
WO2022033391A1 (zh) | 一种超景深显微快速测量装置及测量方法 | |
CN108088381B (zh) | 一种基于图像处理的非接触式微小间隙宽度测量方法 | |
CN204479272U (zh) | 一种焦距测量仪 | |
CN106969719A (zh) | 一种光纤阵列纤芯间距的检测方法及装置 | |
CN103438803B (zh) | 计算机视觉技术跨视场精确测量矩形零件尺寸的方法 | |
CN104880913B (zh) | 一种提高工艺适应性的调焦调平系统 | |
CN210639338U (zh) | 多光谱显微自动聚焦装置 | |
CN105928464B (zh) | 基于图像拼接的自动测量系统及测量方法 | |
CN102003941B (zh) | 基于视觉的大型系泊链五环长测量方法及其测量装置 | |
CN105403170A (zh) | 一种显微3d形貌测量方法及装置 | |
Abdelsalam et al. | PIV camera response to high frequency signal: comparison of CCD and CMOS cameras using particle image simulation | |
CN106556350B (zh) | 显微镜载玻片曲面高度值的测算方法以及一种显微镜 | |
CN112414316B (zh) | 一种应变片敏感栅尺寸参数测量方法 | |
CN102636181A (zh) | 检测高精度空间航拍相机焦距的方法 | |
CN114113115A (zh) | 一种大口径元件表面微缺陷高精度自动定位方法 | |
CN106502074A (zh) | 一种用于像面数字全息显微测量的自动聚焦方法 | |
KR101349447B1 (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 | ||
GR01 | Patent grant | ||
GR01 | Patent grant |