CN102930595B - 基于单幅扫描电子显微镜图像的样品表面三维重建方法 - Google Patents
基于单幅扫描电子显微镜图像的样品表面三维重建方法 Download PDFInfo
- Publication number
- CN102930595B CN102930595B CN201210362683.6A CN201210362683A CN102930595B CN 102930595 B CN102930595 B CN 102930595B CN 201210362683 A CN201210362683 A CN 201210362683A CN 102930595 B CN102930595 B CN 102930595B
- Authority
- CN
- China
- Prior art keywords
- image
- electron microscope
- scanning electron
- pixel
- sample
- 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
Landscapes
- Analysing Materials By The Use Of Radiation (AREA)
Abstract
本发明提供了一种基于单幅扫描电子显微镜图像的样品表面三维重建方法。该方法的特点是仅使用一幅由扫描电子显微镜从由单一元素组成的样品表面上方垂直拍摄的图像,通过分析图像中各个像素点的灰度特征就可以对由单一元素组成的样品的表面形貌进行三维重建。该方法包括:接受来自扫描电子显微镜的图像数据,分析图像灰度变化规律进行归一化处理,根据辐照度方程对归一化处理后的数据进行运算,得到各个点的空间位置并对三维重建结果进行三维可视化。
Description
技术领域
本发明属于微纳加工制造技术领域,具体涉及一种用于对使用扫描电子显微镜观察的具有三维结构的微纳器件进行三维重建的方法。
背景技术
目前的集成电路技术对硅表面的加工已经推进到22纳米,这不但接近光刻的极限,同时也付出了高昂的代价,在集成电路代工厂的设备中光刻机及其相关设备的价格都是千万美元量级的天文数字,投入与产出的矛盾十分突出,大大制约了适用于微纳机电系统的制造技术的发展。因此,业界研究的重点开始转向三维微纳结构(主要是硅结构)。如英特尔公司在2011年推出的三栅极(Tri-Gate)三维晶体管设计,不仅集成度提高而且可以减少50%以上的漏电流。对三维微纳结构的研究都需要一种能够获得待研究样品表面三维结构的方法和技术。目前基于扫描电子显微镜图像对微小结构进行三维重建的方法是使用同一位置多角度的照片对目标样品表面进行重建,典型的是奥地利Alicona公司研制的MEX软件。但对于人造的微纳结构来讲,微纳结构所在的载体通常很小(如单片CPU)难以对同一位置进行多角度拍摄,因此利用公知的方法难以完成对三维微纳结构三维表面的获取。
发明内容
本发明目的是解决现有方法难以完成对三维微纳结构三维表面获取的问题,提供一种基于单幅扫描电子显微镜图像的样品表面三维重建方法,作为三维微纳加工工艺研发定量的观察验证手段,同时也为三维微纳结构的应用提供研究材料。
本发明提供了一种基于单幅垂直俯视视角扫描电子显微镜图像的三维重建方法,所述方法仅使用一幅垂直俯视视角的扫描电子显微镜图像,通过分析图像中各个像素点的灰度特征就可以对由单一元素组成的样品的表面形貌进行三维重建。
所述方法具体包括以下步骤:
第1、接受来自扫描电子显微镜的图像数据;该图像是由扫描电子显微镜从由单一元素组成的样品表面的上方垂直拍摄的图像。
具体步骤是:将由单一元素组成的样品水平放置在扫描电子显微镜的载物台上,使用扫描电子显微镜从样品的正上方对目标区域进行拍摄;通过与扫描电子显微镜相连的计算机获得拍摄照片;
第2、分析图像灰度变化特征,寻找整幅图像灰度的最小值作为图像特征值k;
第3、使用图像特征值k对扫描电子显微镜图像进行归一化处理,归一化时使用的公式如下所示:
A(x,y)=E(x,y)/k
其中,(x,y)是指扫描电子显微镜图像中正在进行归一化处理的像素点的横纵坐标,E(x,y)是指该位置像素点在扫描电子显微镜图像中的灰度值,k是图像特征值,A(x,y)是指图像中该像素对应的实际样品表面位置被扫描电子显微镜激发的电子数量描述值,对整幅图像依据这一公式进行运算,就可以得到实际样品表面被扫描电子显微镜激发的电子数量描述矩阵A;
第4、根据辐照度方程对电子数量描述矩阵A进行运算,进而得到各个点以像素为单位的空间位置坐标,随后根据扫描电子显微镜拍摄时所使用的比例尺将像素单位换算成实际长度单位,所使用的辐照度方程如下所示:
∫∫(A(x,y)-R(p(x,y),q(x,y)))2dxdy
其中,p(x,y)是指图像中(x,y)像素点对应的实际样品表面位置的沿x方向的梯度,q(x,y)是指图像中(x,y)像素点对应的实际样品表面位置的沿y方向的梯度。R(p(x,y),q(x,y))是指根据电子激发的理论公式计算出来的图像(x,y)像素点对应实际样品表面位置被扫描电子显微镜激发的电子数量描述值。对于不同的p(x,y)和q(x,y)值可以根据公式得到不同R(p(x,y),q(x,y))值,所有点的R(p(x,y),q(x,y))与A(x,y)间的差距总和越小就说明此时的p(x,y)和q(x,y)值越与实际相符合。当R(p(x,y),q(x,y))与A(x,y)间的差距总和最小时,即辐照度方程取得最小值时,即可得到各个点的空间位置。因此通过求取辐照度方程的最小值就可以获得扫描电子显微镜图像的三维重建结果。这一过程中所使用的电子激发的公式如下所示:
其中n=(p(x,y),q(x,y),-1)是指图像中(x,y)像素点对应的实际样品表面位置的法向量,ns是指扫描电子显微镜进行拍摄时在所使用的扫描电子束的入射方向。由于扫描电子显微镜从样品上方垂直向下进行拍摄,所以ns=(0,0,-1)。对于图像中各个像素点都使用电子激发的理论公式进行运算并带入辐照度方程,最终通过求取辐照度方程最小值来获得图像中各个像素点的高度,即得到各个点的空间位置;
第5、对三维重建结果进行三维可视化显示。
本发明的优点和积极效果:
本发明实现了一种基于单幅垂直俯视视角扫描电子显微镜图像对单一元素组成的样品表面进行三维重建的新方法。能够方便地对载体较小的三维微纳结构(如黑硅表面微纳结构)进行三位重建,为三维微纳加工工艺研发定量的观察验证手段、同时也为三维微纳结构的应用提供研究材料。
附图说明
图1为实施例1所使用的扫描电子显微镜图像。
图2为本发明使用扫描电子显微镜设备拍摄样本的示意图.
图3为本发明基于单幅垂直俯视视角扫描电子显微镜图像的三维重建方法的流程图。
图4为最终得到的三维重建结果示意图。
具体实施方式
实施例1
本实施例所使用的样品是表面刻蚀有微纳结构的黑硅片,本发明提供的基于单幅扫描电子显微镜图像的样品表面三维重建方法可以为微纳结构进行三维重建。具体实施情况如下:
1、将由单一的硅元素组成的黑硅样品水平放置在扫描电子显微镜的载物台上,使用扫描电子显微镜从样品的正上方对目标区域进行拍摄;通过与扫描电子显微镜相连的计算机获得拍摄照片如图1,是一幅216*222个像素的图像;
图2所示为使用扫描电子显微镜对待重建的样品进行拍摄以获取待重建区域图的操作示意图,其中待重建的样品必须由单一元素组成。拍摄时需要保证样品底部水平,黑硅片的底部通常十分平整,这时只需保证载物台A水平,这时拍摄的扫描电子显微镜照片为垂直俯视视角。与此同时,可以在与扫描电子显微镜相连的计算机中获得待重建的图像在显示器B中显示。
随后对待重建图像操作以获得重建结果,具体的步骤如图3所示。
2、进行分析图像灰度变化特征,寻找整幅图像灰度的最小值(本例最小值为18),作为图像特征值k;
3、使用图像特征值k对扫描电子显微镜图像进行归一化处理,归一化时使用的公式如下所示:
A(x,y)=E(x,y)/k
其中,(x,y)是指扫描电子显微镜图像中正在进行归一化处理的像素点的横纵坐标,E(x,y)是指该位置像素点在扫描电子显微镜图像中的灰度值,k是图像特征值,A(x,y)是指图像中该像素对应的实际样品表面位置被扫描电子显微镜激发的电子数量描述值,对整幅图像依据这一公式进行运算,就可以得到实际样品表面被扫描电子显微镜激发的电子数量描述矩阵A;
在本例中,在(x,y)=(1,1)时,A(x,y)=3,(x,y)=(2,2)时,A(x,y)=3以此类推,最终获得描述矩阵 共216行,222列。
4、根据辐照度方程对电子数量描述矩阵A进行运算,进而得到各个点以像素为单位的空间位置坐标,随后根据扫描电子显微镜拍摄时所使用的比例尺将像素单位换算成实际长度单位,本例中所使用的比例尺是,200像素比5微米,最终得到空间位置矩阵为:
在计算所使用的辐照度方程如下所示:
∫∫(A(x,y)-R(p(x,y),q(x,y)),1)2dxdy
其中,(x,y)是指扫描电子显微镜图像中正在进行归一化处理的像素点的横纵坐标(扫描电子显微镜图像的位置与实际位置存在一一对应的关系),A(x,y)是指图像中(x,y)像素点对应的实际样品表面位置被扫描电子显微镜激发的电子数量描述值,(p(x,y),q(x,y),-1)是指图像中(x,y)像素点对应的实际样品表面位置的坡度向量,R(p(x,y),q(x,y))是指根据图像中(x,y)像素点对应的实际样品表面位置的坡度和电子激发的理论公式计算出来的被扫描电子显微镜激发的电子数量描述值,所使用的电子激发的理论公式如下所示:
其中n=(p(x,y),q(x,y),-1)是指图像中(x,y)像素点对应的实际样品表面位置的坡度向量,ns是指扫描电子显微镜进行拍摄时在所使用的扫描电子束的入射方向,由于扫描电子显微镜从样品上方垂直向下进行拍摄,所以ns=(0,0,-1),对于图像中各个像素点都使用电子激发的理论公式进行运算并带入辐照度方程,最终通过求取辐照度方程最小值来获得图像中各个像素点的高度,即得到各个点的空间位置。如果需要可根据扫描电子显微镜拍照时所使用的比例尺对重建结果进行像素和实际距离之间进行换算;
5、得到的三维重建结果将被三维可视化并显示在计算机的显示器上,本例的三维重建结果如图4所示。
以上实施方式仅用于说明本发明,而并非对本发明的限制,有关技术领域的普通技术人员,在不脱离本发明的精神和范围的情况,还可以做出各种变化,因此所有的等同的技术方案也属于本发明的范畴,本发明的专利保护范围应由其权利要求限定。
Claims (3)
1.一种基于单幅扫描电子显微镜图像的样品表面三维重建方法,其特征在于该方法包括:
第1、扫描电子显微镜从样品上方垂直向下进行拍摄,计算机接受来自扫描电子显微镜的图像数据;
第2、分析图像灰度变化特征,寻找整幅图像灰度的最小值作为图像特征值 ;
第3、使用图像特征值对扫描电子显微镜图像进行归一化处理,归一化时使用的公式如下所示:
其中,是指扫描电子显微镜图像中正在进行归一化处理的像素点的横纵坐标,是指点在扫描电子显微镜图像中的灰度值,是指图像中点对应的样品表面位置被扫描电子显微镜激发的电子数量的归一化值,对整幅图像依据这一公式进行运算,就能够得到样品表面被扫描电子显微镜激发的电子数量的归一化矩阵;
第4、根据辐照度方程对电子数量的归一化矩阵A进行运算,进而得到各个点以像素为单位的空间位置坐标,随后根据扫描电子显微镜拍摄时所使用的比例尺将像素单位换算成实际长度单位,所使用的辐照度方程如下所示:
其中,是指图像中像素点对应的样品表面位置的沿方向的梯度,是指图像中像素点对应的样品表面位置的沿方向的梯度,是指根据电子激发的理论公式计算出来的图像像素点对应样品表面位置被扫描电子显微镜激发的电子数量描述值;对于不同的和值可以根据公式得到不同值,所有点的与间的差距总和越小就说明此时的和值越与实际相符合;当与间的差距总和最小时,即辐照度方程取得最小值时,即可得到各个点的空间位置;因此通过求取辐照度方程的最小值就可以获得扫描电子显微镜图像的三维重建结果,这一过程中所使用的电子激发的理论公式如下所示:
其中是指图像中像素点对应的样品表面位置的法向量,是指扫描电子显微镜进行拍摄时在所使用的扫描电子束的入射方向,由于扫描电子显微镜从样品上方垂直向下进行拍摄,所以);对于图像中各个像素点都使用该电子激发的理论公式进行运算并带入辐照度方程,最终通过求取辐照度方程最小值来获得图像中各个像素点的高度,即得到各个点的空间位置;
第5、对三维重建结果进行三维可视化显示。
2.根据权利要求1所述的方法,其特征在于,第1步所述样品由单一元素组成。
3.根据权利要求2所述的方法,其特征在于,对黑硅表面进行三维重建时待重建的样品黑硅均由硅元素组成。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210362683.6A CN102930595B (zh) | 2012-09-25 | 2012-09-25 | 基于单幅扫描电子显微镜图像的样品表面三维重建方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210362683.6A CN102930595B (zh) | 2012-09-25 | 2012-09-25 | 基于单幅扫描电子显微镜图像的样品表面三维重建方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102930595A CN102930595A (zh) | 2013-02-13 |
CN102930595B true CN102930595B (zh) | 2015-04-08 |
Family
ID=47645385
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201210362683.6A Active CN102930595B (zh) | 2012-09-25 | 2012-09-25 | 基于单幅扫描电子显微镜图像的样品表面三维重建方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102930595B (zh) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103295203A (zh) * | 2013-06-12 | 2013-09-11 | 中国科学院光电技术研究所 | 一种基于b样条的硅片表面形貌重构的方法 |
CN105590338B (zh) * | 2015-12-07 | 2018-08-10 | 中国科学院微电子研究所 | 一种扫描电子显微图像的三维重构方法 |
CN111191624B (zh) * | 2020-01-04 | 2021-02-19 | 圣点世纪科技股份有限公司 | 一种基于法向量信息的指静脉特征提取方法 |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101236659A (zh) * | 2008-02-26 | 2008-08-06 | 天津商业大学 | 基于相关识别的三维表面重建方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7885483B2 (en) * | 2006-11-24 | 2011-02-08 | National Synchrotron Radiation Research Center | Image alignment method |
-
2012
- 2012-09-25 CN CN201210362683.6A patent/CN102930595B/zh active Active
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101236659A (zh) * | 2008-02-26 | 2008-08-06 | 天津商业大学 | 基于相关识别的三维表面重建方法 |
Non-Patent Citations (2)
Title |
---|
Environmental scanning electron microscopy in colour;G.D.Danilatos等;《Journal of Microscopy》;19860630;第142卷(第3期);第317-325页 * |
基于显微图像处理的MEMS运动测量方法;赵新等;《纳米技术与精密工程》;20120115;第10卷(第1期);第1-5页 * |
Also Published As
Publication number | Publication date |
---|---|
CN102930595A (zh) | 2013-02-13 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP5393550B2 (ja) | 走査荷電粒子顕微鏡を用いた画像生成方法及び装置、並びに試料の観察方法及び観察装置 | |
Zhu et al. | 3D nanostructure reconstruction based on the SEM imaging principle, and applications | |
CN106062824B (zh) | 边缘检测装置和边缘检测方法 | |
Guery et al. | Characterization of SEM speckle pattern marking and imaging distortion by digital image correlation | |
US9823458B2 (en) | Imaging system and method for multi-scale three-dimensional deformation or profile output | |
CN106709958A (zh) | 一种基于灰度梯度和颜色直方图的图像质量评价方法 | |
CN103076347B (zh) | 基于原位x射线断层照相的脆性材料力学损伤的测量方法 | |
Rossi et al. | Assessment of the metrological performance of an in situ storage image sensor ultra-high speed camera for full-field deformation measurements | |
CN102930595B (zh) | 基于单幅扫描电子显微镜图像的样品表面三维重建方法 | |
CN103778616A (zh) | 基于区域的对比度金字塔图像融合方法 | |
CN109285137A (zh) | 页岩孔隙贡献度的获取方法及装置 | |
CN106062815A (zh) | 使用准周期栅格的图像校正的方法及系统 | |
Gontard et al. | Assessment of engineered surfaces roughness by high-resolution 3D SEM photogrammetry | |
KR20150003783A (ko) | X선 검사 방법 및 x선 검사 장치 | |
Toeneboehn et al. | Stereovision combined with particle tracking velocimetry reveals advection and uplift within a restraining bend simulating the Denali fault | |
Dai et al. | The geometric phase analysis method based on the local high resolution discrete Fourier transform for deformation measurement | |
Chen et al. | Intelligent adaptive sampling guided by Gaussian process inference | |
Zhu et al. | A computer vision approach to study surface deformation of materials | |
US10324021B2 (en) | Magnetophorisis measuring system for determining motion status of object and quantifying amount of magnetic particles contained therein | |
WO2021048104A1 (en) | Wafer inspection methods and systems | |
Wang et al. | Resolution enhancement for topography measurement of high-dynamic-range surfaces via image fusion | |
Senin et al. | The use of areal surface topography analysis for the inspection of micro-fabricated thin foil laser targets for ion acceleration | |
CN107578001B (zh) | 测试指纹采集设备的分辨率的方法及装置 | |
CN106122430A (zh) | 一种基于特征图像的小模数齿轮边缘检测精度计算方法 | |
CN107590787A (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 |