CN106875439B - 基于三维点云模型的单晶硅棒外形尺寸测量方法 - Google Patents
基于三维点云模型的单晶硅棒外形尺寸测量方法 Download PDFInfo
- 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
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/60—Analysis of geometric attributes
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01B—MEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
- G01B21/00—Measuring arrangements or details thereof, where the measuring technique is not covered by the other groups of this subclass, unspecified or not relevant
- G01B21/10—Measuring 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
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01B—MEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
- G01B21/00—Measuring arrangements or details thereof, where the measuring technique is not covered by the other groups of this subclass, unspecified or not relevant
- G01B21/20—Measuring 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
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/10—Image acquisition modality
- G06T2207/10028—Range image; Depth image; 3D point clouds
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/30—Subject of image; Context of image processing
- G06T2207/30108—Industrial image inspection
- G06T2207/30148—Semiconductor; 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坐标的差值即为该单晶硅棒的高度值。
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)
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)
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)
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 |
-
2017
- 2017-02-20 CN CN201710088274.4A patent/CN106875439B/zh not_active Expired - Fee Related
Patent Citations (6)
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)
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 |