CN106875439B - 基于三维点云模型的单晶硅棒外形尺寸测量方法 - Google Patents

基于三维点云模型的单晶硅棒外形尺寸测量方法 Download PDF

Info

Publication number
CN106875439B
CN106875439B CN201710088274.4A CN201710088274A CN106875439B CN 106875439 B CN106875439 B CN 106875439B CN 201710088274 A CN201710088274 A CN 201710088274A CN 106875439 B CN106875439 B CN 106875439B
Authority
CN
China
Prior art keywords
single crystal
point
point cloud
silicon single
coordinate
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.)
Expired - Fee Related
Application number
CN201710088274.4A
Other languages
English (en)
Other versions
CN106875439A (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 CN201710088274.4A priority Critical patent/CN106875439B/zh
Publication of CN106875439A publication Critical patent/CN106875439A/zh
Application granted granted Critical
Publication of CN106875439B publication Critical patent/CN106875439B/zh
Expired - Fee Related 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/60Analysis of geometric attributes
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01BMEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
    • G01B21/00Measuring arrangements or details thereof, where the measuring technique is not covered by the other groups of this subclass, unspecified or not relevant
    • G01B21/10Measuring arrangements or details thereof, where the measuring technique is not covered by the other groups of this subclass, unspecified or not relevant for measuring diameters
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01BMEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
    • G01B21/00Measuring arrangements or details thereof, where the measuring technique is not covered by the other groups of this subclass, unspecified or not relevant
    • G01B21/20Measuring arrangements or details thereof, where the measuring technique is not covered by the other groups of this subclass, unspecified or not relevant for measuring contours or curvatures, e.g. determining profile
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/10Image acquisition modality
    • G06T2207/10028Range image; Depth image; 3D point clouds
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/30Subject of image; Context of image processing
    • G06T2207/30108Industrial image inspection
    • G06T2207/30148Semiconductor; IC; Wafer

Abstract

一种基于三维点云模型的单晶硅棒外形尺寸测量方法:垂直于单晶硅棒三维点云模型的轴线方向进行分层,共计N层;分别计算每一层横截面轮廓点云的几何中心点;对N个几何中心点进行直线拟合,得到单晶硅棒的几何中心轴线;计算单晶硅棒点云模型中任一点到拟合几何中心轴线的距离;所得单晶硅棒三维点云模型中最大y坐标和最小y坐标的差值即为该单晶硅棒的高度值。本发明可以有效避免由于单晶硅棒表面形状复杂多变以及工人通过游标卡尺等工具手工测量产生的误差。该专利算法简便易行,复杂度低,准确性高,数据易于保存和传输,可实现自动化测量,可有效提高单晶硅棒外形尺寸的测量效率和精度,为单晶硅片生产提供高效的质量保证。

Description

基于三维点云模型的单晶硅棒外形尺寸测量方法
技术领域
本发明涉及一种单晶硅棒外形尺寸测量方法。特别是涉及一种基于三维点云模型的单晶硅棒外形尺寸测量方法。
背景技术
从单晶硅棒加工成单晶硅抛光硅片的工艺流程主要包括切断、外径滚磨、平边或V型槽处理及切片的过程。经过切断处理后,单晶硅棒的两个断面平整且相互平行,但是通过单晶生长得到的单晶硅棒表面形状起伏变化大,直径不均匀,因此在外径滚磨前需要对单晶硅棒进行测量,得到硅棒最细直径值,从而预判该单晶硅棒是否能够加工成合格的单晶硅片,或者可加工成多大直径的单晶硅片。目前该行业一般采用工人通过游标卡尺等工具进行测量的方式完成。由于单晶硅棒表面的形状复杂性,这样测量容易偏离真实值,为后续加工工序带来不必要的麻烦。
发明内容
本发明所要解决的技术问题是,提供一种能够快速、准确地实现单晶硅棒外形尺寸,特别是最大、最小直径尺寸及相对位置三维测量的基于三维点云模型的单晶硅棒外形尺寸测量方法。
本发明所采用的技术方案是:一种基于三维点云模型的单晶硅棒外形尺寸测量方法,包括如下步骤:
1)垂直于单晶硅棒三维点云模型的轴线方向,即y坐标方向等间隔的对点云进行分层,得到平行于xoz坐标平面的横截面轮廓点云,共计N层;
2)分别计算每一层横截面轮廓点云的几何中心点,记为Ci(xcen,i,yi,zcen,i),i=1,…,N;
3)对步骤2)中得到的N个几何中心点进行直线拟合,得到的拟合直线即为单晶硅棒的几何中心轴线;
4)计算单晶硅棒点云模型中任一点到拟合几何中心轴线的距离,并记录最大距离值、最小距离值和对应点的y坐标,即测得单晶硅棒的最大直径、最小直径和位置;
5)所得单晶硅棒三维点云模型中最大y坐标和最小y坐标的差值即为该单晶硅棒的高度值。
步骤2)所述的计算每一层横截面轮廓点云的几何中心点坐标,是采用如下公式:
式中,xcen,i是几何中心点的x坐标;zcen,i是几何中心点的z坐标;Mi是第i层横截面轮廓点云的点的个数;xij是第i层横截面轮廓点云第j个点的x坐标;zij是第i层横截面轮廓点云第j个点的z坐标。
步骤3)拟合后的几何中心轴线方程如下:
式中m,l,n是几何中心轴线的方向向量,x0,y0,z0是拟合后几何中心轴线上任一点的x、y、z坐标分量。
步骤4)所述的计算单晶硅棒点云模型中任一点到拟合几何中心轴线的距离,是采用如下公式:
式中xk,yk,zk是单晶硅棒点云模型中任一点的x、y、z坐标分量。
本发明的基于三维点云模型的单晶硅棒外形尺寸测量方法,利用单晶硅棒的三维点云模型,采用几何建模的方法通过计算机程序自动实现,可以有效避免由于单晶硅棒表面形状复杂多变以及工人通过游标卡尺等工具手工测量产生的误差。该专利算法简便易行,复杂度低,准确性高,数据易于保存和传输,可实现自动化测量,可有效提高单晶硅棒外形尺寸的测量效率和精度,为单晶硅片生产提供高效的质量保证。
附图说明
图1是本发明的基于三维点云模型的单晶硅棒外形尺寸测量方法流程图;
图2是本发明中单晶硅棒三维点云模型、坐标系和水平横截面;
图3是第i层轮廓线点云及其几何中心坐标;
图4初始几何中心轴线和拟合后的几何中心轴线;
图5单晶硅棒原始点云(a)、最大直径对应圆柱(b)和最小直径对应圆柱(c)的示意图。
具体实施方式
下面结合实施例和附图对本发明的基于三维点云模型的单晶硅棒外形尺寸测量方法做出详细说明。
如图1所示,本发明的基于三维点云模型的单晶硅棒外形尺寸测量方法,包括如下步骤:
1)垂直于单晶硅棒三维点云模型的轴线方向,即y坐标方向等间隔的对点云进行分层,得到平行于xoz坐标平面的横截面轮廓点云,共计N层;
2)分别计算每一层横截面轮廓点云的几何中心点,记为Ci(xcen,i,yi,zcen,i),i=1,…,N;
所述的计算每一层横截面轮廓点云的几何中心点坐标,是采用如下公式:
式中,xcen,i是几何中心点的x坐标;zcen,i是几何中心点的z坐标;Mi是第i层横截面轮廓点云的点的个数;xij是第i层横截面轮廓点云第j个点的x坐标;zij是第i层横截面轮廓点云第j个点的z坐标。
3)对步骤2)中得到的N个几何中心点进行直线拟合,得到的拟合直线即为单晶硅棒的几何中心轴线;拟合后的几何中心轴线方程如下:
式中m,l,n是几何中心轴线的方向向量,x0,y0,z0是拟合后几何中心轴线上任一点的x、y、z坐标分量。
4)计算单晶硅棒点云模型中任一点到拟合几何中心轴线的距离,并记录最大距离值、最小距离值和对应点的y坐标,即测得单晶硅棒的最大直径、最小直径和位置;
5)所得单晶硅棒三维点云模型中最大y坐标和最小y坐标的差值即为该单晶硅棒的高度值。
所述的计算单晶硅棒点云模型中任一点到拟合几何中心轴线的距离,是采用如下公式:
式中xk,yk,zk是单晶硅棒点云模型中任一点的x、y、z坐标分量。
下面给了具体实例:
(1)如图2所示,垂直于单晶硅棒三维点云模型的高度方向(y坐标方向),以4mm间隔对点云进行分层。第i层的y坐标值yi为2mm+i×4mm,y坐标在yi±2mm范围内的点同属于一层轮廓点云,共得到N层轮廓线点云。这样处理后认为同层轮廓线点云具有相同的y坐标.
(2)图3是第i层轮廓点云,通过式(1)计算该层轮廓线点云的几何中心坐标,式中Mi是该层轮廓线中点的个数,(xj,zj)是该层轮廓线点云中的任意一点。计算得到的中心坐标(xcen,i,zcen,i)如图3中的中心圆点所示。
(3)对步骤(2)中得到的所有N个中心点进行直线拟合,得到拟合后的几何中心轴线,其直线方程可用式(2)表示。
(4)利用式(3)计算单晶硅棒点云模型中任一点pk(xk,yk,zk)到拟合后几何中心线的距离dk。对所有距离进行排序,找出最大离值、最小距离值和对应点的y坐标,即测得该单晶硅棒的最大、最小直径和位置.
(5)利用如下式(4)得到该单晶硅棒的高度值。
H=ymax-ymin (4)
图4是初始几何中心轴线和拟合后的几何中心轴线,离散点是步骤(2)得到的各层几何中心点,中间的细线是步骤(3)拟合得到的几何中心轴线。图5是测量结束后三维可视化显示结果,标记a指示的是单晶硅棒原始点云,标记b指示的是最大直径对应的外接圆柱,标记c指示的是最小直径对应的内接圆柱。表1给出了该单晶硅棒的具体测量结果和测量时间。
表1图2所示点云模型对应单晶硅棒外形关键尺寸测量结果

Claims (1)

1.一种基于三维点云模型的单晶硅棒外形尺寸测量方法,其特征在于,包括如下步骤:
1)垂直于单晶硅棒三维点云模型的轴线方向,即y坐标方向等间隔的对点云进行分层,得到平行于xoz坐标平面的横截面轮廓点云,共计N层;
2)分别计算每一层横截面轮廓点云的几何中心点,记为Ci(xcen,i,yi,zcen,i),i=1,…,N;
所述的计算每一层横截面轮廓点云的几何中心点坐标,是采用如下公式:
式中,xcen,i是几何中心点的x坐标;zcen,i是几何中心点的z坐标;Mi是第i层横截面轮廓点云的点的个数;xij是第i层横截面轮廓点云第j个点的x坐标;zij是第i层横截面轮廓点云第j个点的z坐标;
3)对步骤2)中得到的N个几何中心点进行直线拟合,得到的拟合直线即为单晶硅棒的几何中心轴线;拟合后的几何中心轴线方程如下:
式中m,l,n是几何中心轴线的方向向量,x0,y0,z0是拟合后几何中心轴线上任一点的x、y、z坐标分量;
4)计算单晶硅棒点云模型中任一点到拟合几何中心轴线的距离,并记录最大距离值、最小距离值和对应点的y坐标,即测得单晶硅棒的最大直径、最小直径和位置;
所述的计算单晶硅棒点云模型中任一点到拟合几何中心轴线的距离,是采用如下公式:
式中xk,yk,zk是单晶硅棒点云模型中任一点的x、y、z坐标分量,m,l,n是几何中心轴线的方向向量,x0,y0,z0是拟合后几何中心轴线上任一点的x、y、z坐标分量;
5)所得单晶硅棒三维点云模型中最大y坐标和最小y坐标的差值即为该单晶硅棒的高度值。
CN201710088274.4A 2017-02-20 2017-02-20 基于三维点云模型的单晶硅棒外形尺寸测量方法 Expired - Fee Related CN106875439B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710088274.4A CN106875439B (zh) 2017-02-20 2017-02-20 基于三维点云模型的单晶硅棒外形尺寸测量方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710088274.4A CN106875439B (zh) 2017-02-20 2017-02-20 基于三维点云模型的单晶硅棒外形尺寸测量方法

Publications (2)

Publication Number Publication Date
CN106875439A CN106875439A (zh) 2017-06-20
CN106875439B true CN106875439B (zh) 2019-06-07

Family

ID=59167104

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710088274.4A Expired - Fee Related CN106875439B (zh) 2017-02-20 2017-02-20 基于三维点云模型的单晶硅棒外形尺寸测量方法

Country Status (1)

Country Link
CN (1) CN106875439B (zh)

Families Citing this family (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108332713B (zh) * 2018-01-16 2020-11-13 安徽慧视金瞳科技有限公司 一种绝缘子引脚倾斜的检测方法
CN108846888B (zh) * 2018-04-23 2022-03-29 北京建筑大学 一种古木建筑构件精细尺寸信息自动化提取方法
CN111272091B (zh) * 2020-03-09 2021-03-23 东南大学 一种基于三维激光扫描的焊接工字钢大变形自动识别方法
CN111366113A (zh) * 2020-03-31 2020-07-03 陕西梵翌琨机电科技有限公司 单晶硅棒测量模块及其晶体生长系统
CN111750808B (zh) * 2020-06-28 2021-09-14 电子科技大学 基于激光雷达扫描仪的地表粗糙度获取方法
CN115795579B (zh) * 2022-12-23 2023-06-27 岭南师范学院 无特征复杂曲面测量误差分析的快速坐标对齐方法

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102855663A (zh) * 2012-05-04 2013-01-02 北京建筑工程学院 一种根据激光雷达栅格点云构建csg模型的方法
CN104359459A (zh) * 2014-12-04 2015-02-18 上海岩土工程勘察设计研究院有限公司 采用三维激光扫描反射率信息生成隧道内壁影像的方法
CN104457574A (zh) * 2014-12-11 2015-03-25 天津大学 一种非接触式测量不规则物体体积的装置和测量方法
CN105004274A (zh) * 2015-07-07 2015-10-28 深圳大学 一种基于三维视觉的圆柱面半径测量方法
GB2537681A (en) * 2015-04-24 2016-10-26 Isis Innovation A method of detecting objects within a 3D environment
CN106248035A (zh) * 2016-08-19 2016-12-21 苏州大学 一种基于点云模型的面轮廓度精确检测的方法及系统

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20140157579A1 (en) * 2012-12-08 2014-06-12 8 Tree, Llc Networked marketplace for custom 3D fabrication

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102855663A (zh) * 2012-05-04 2013-01-02 北京建筑工程学院 一种根据激光雷达栅格点云构建csg模型的方法
CN104359459A (zh) * 2014-12-04 2015-02-18 上海岩土工程勘察设计研究院有限公司 采用三维激光扫描反射率信息生成隧道内壁影像的方法
CN104457574A (zh) * 2014-12-11 2015-03-25 天津大学 一种非接触式测量不规则物体体积的装置和测量方法
GB2537681A (en) * 2015-04-24 2016-10-26 Isis Innovation A method of detecting objects within a 3D environment
CN105004274A (zh) * 2015-07-07 2015-10-28 深圳大学 一种基于三维视觉的圆柱面半径测量方法
CN106248035A (zh) * 2016-08-19 2016-12-21 苏州大学 一种基于点云模型的面轮廓度精确检测的方法及系统

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
3D Moving Object Reconstruction by Temporal Accumulation;Anas Abuzaina etc.;《2014 22nd International Conference on Pattern Recognition》;20141231;全文
圆柱面拟合方法研究;王崇倡 等;《测绘工程》;20140331;第23卷(第3期);全文
基于人体特征提取的模特体型尺寸自动测量方法;葛宝臻 等;《纺 织 学 报》;20120430;第33卷(第4期);全文

Also Published As

Publication number Publication date
CN106875439A (zh) 2017-06-20

Similar Documents

Publication Publication Date Title
CN106875439B (zh) 基于三维点云模型的单晶硅棒外形尺寸测量方法
CN106799745A (zh) 一种基于协同克里金的工业机械臂精度标定方法
CN100538261C (zh) 基于探路法的未知自由曲面自适应测量方法与测头装置
CN103777570B (zh) 基于nurbs曲面的加工误差快速检测补偿方法
CN104359404B (zh) 飞机零件上大量小尺寸导孔的快速视觉检测方法
CN103217688B (zh) 一种基于不规则三角网机载激光雷达点云平差计算方法
JP2002059340A (ja) 移動装置の位置誤差評価方法およびその評価結果に基づく移動精度向上方法
CN102622479B (zh) 一种基于三维草图的逆向工程cad建模方法
CN103302550A (zh) 透镜形状加工方法以及透镜形状加工装置
CN106406237B (zh) 一种具有自由曲面金属零件的加工方法
CN107063060A (zh) 一种确定表面平面度的方法及装置
CN108917600A (zh) 一种基于三维反求的复杂结构毛坯加工定位方法
CN102049731A (zh) 对盘形凸轮进行精确轮廓测量及精确加工的方法
CN110103071A (zh) 一种变形复杂零件的数字化寻位加工方法
CN102853754B (zh) 一种快速实现电感量仪检测叶片型面的方法
CN115026706B (zh) 一种航空发动机叶片抛光方法及系统
CN106568379A (zh) 一种基于空间定位的测量飞机部件对接面外形的方法
CN103591912B (zh) 一种环形零件异型特征位置度的测量方法
CN108645301A (zh) 一种直齿轮齿面偏差在机测量方法
CN109059821A (zh) 坐标测量机测量路径规划方法
CN111664801A (zh) 一种基于高斯曲率的光学靶标点投影规划方法
CN110470223A (zh) 一种基于1d位移传感器的三维空间坐标测量方法
CN108332784B (zh) 一种距离测量验证方法
CN106772332B (zh) 一种旋转式距离测量验证方法
US10132623B2 (en) Method for measuring slant wall thickness dimension of hub

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20190607

Termination date: 20210220

CF01 Termination of patent right due to non-payment of annual fee