CN106875439A - 基于三维点云模型的单晶硅棒外形尺寸测量方法 - Google Patents
基于三维点云模型的单晶硅棒外形尺寸测量方法 Download PDFInfo
- Publication number
- CN106875439A CN106875439A CN201710088274.4A CN201710088274A CN106875439A CN 106875439 A CN106875439 A CN 106875439A CN 201710088274 A CN201710088274 A CN 201710088274A CN 106875439 A CN106875439 A CN 106875439A
- Authority
- CN
- China
- Prior art keywords
- single crystal
- point cloud
- crystal silicon
- coordinate
- point
- 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
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
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Geometry (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Theoretical Computer Science (AREA)
- Crystals, And After-Treatments Of Crystals (AREA)
- Length Measuring Devices By Optical Means (AREA)
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 (4)
1.一种基于三维点云模型的单晶硅棒外形尺寸测量方法,其特征在于,包括如下步骤:
1)垂直于单晶硅棒三维点云模型的轴线方向,即y坐标方向等间隔的对点云进行分层,得到平行于xoz坐标平面的横截面轮廓点云,共计N层;
2)分别计算每一层横截面轮廓点云的几何中心点,记为Ci(xcen,i,yi,zcen,i),i=1,…,N;
3)对步骤2)中得到的N个几何中心点进行直线拟合,得到的拟合直线即为单晶硅棒的几何中心轴线;
4)计算单晶硅棒点云模型中任一点到拟合几何中心轴线的距离,并记录最大距离值、最小距离值和对应点的y坐标,即测得单晶硅棒的最大直径、最小直径和位置;
5)所得单晶硅棒三维点云模型中最大y坐标和最小y坐标的差值即为该单晶硅棒的高度值。
2.根据权利要求1所述的基于三维点云模型的单晶硅棒外形尺寸测量方法,其特征在于,步骤2)所述的计算每一层横截面轮廓点云的几何中心点坐标,是采用如下公式:
式中,xcen,i是几何中心点的x坐标;zcen,i是几何中心点的z坐标;Mi是第i层横截面轮廓点云的点的个数;xij是第i层横截面轮廓点云第j个点的x坐标;zij是第i层横截面轮廓点云第j个点的z坐标。
3.根据权利要求1所述的基于三维点云模型的单晶硅棒外形尺寸测量方法,其特征在于,步骤3)拟合后的几何中心轴线方程如下:
式中m,l,n是几何中心轴线的方向向量,x0,y0,z0是拟合后几何中心轴线上任一点的x、y、z坐标分量。
4.根据权利要求1所述的基于三维点云模型的单晶硅棒外形尺寸测量方法,其特征在于,步骤4)所述的计算单晶硅棒点云模型中任一点到拟合几何中心轴线的距离,是采用如下公式:
式中xk,yk,zk是单晶硅棒点云模型中任一点的x、y、z坐标分量。
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 true CN106875439A (zh) | 2017-06-20 |
CN106875439B 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) |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108332713A (zh) * | 2018-01-16 | 2018-07-27 | 安徽慧视金瞳科技有限公司 | 一种绝缘子引脚倾斜的检测方法 |
CN108846888A (zh) * | 2018-04-23 | 2018-11-20 | 北京建筑大学 | 一种古木建筑构件精细尺寸信息自动化提取方法 |
CN111272091A (zh) * | 2020-03-09 | 2020-06-12 | 东南大学 | 一种基于三维激光扫描的焊接工字钢大变形自动识别方法 |
CN111366113A (zh) * | 2020-03-31 | 2020-07-03 | 陕西梵翌琨机电科技有限公司 | 单晶硅棒测量模块及其晶体生长系统 |
CN111750808A (zh) * | 2020-06-28 | 2020-10-09 | 电子科技大学 | 基于激光雷达扫描仪的地表粗糙度获取方法 |
CN113628166A (zh) * | 2021-07-12 | 2021-11-09 | 广东电网有限责任公司广州供电局 | 一种开关柜断路器室的几何尺寸检测方法及装置 |
CN115795579A (zh) * | 2022-12-23 | 2023-03-14 | 岭南师范学院 | 一种无特征复杂曲面误差分析的快速坐标对齐方法 |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102855663A (zh) * | 2012-05-04 | 2013-01-02 | 北京建筑工程学院 | 一种根据激光雷达栅格点云构建csg模型的方法 |
US20140157579A1 (en) * | 2012-12-08 | 2014-06-12 | 8 Tree, Llc | Networked marketplace for custom 3D fabrication |
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 | 苏州大学 | 一种基于点云模型的面轮廓度精确检测的方法及系统 |
-
2017
- 2017-02-20 CN CN201710088274.4A patent/CN106875439B/zh not_active Expired - Fee Related
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102855663A (zh) * | 2012-05-04 | 2013-01-02 | 北京建筑工程学院 | 一种根据激光雷达栅格点云构建csg模型的方法 |
US20140157579A1 (en) * | 2012-12-08 | 2014-06-12 | 8 Tree, Llc | Networked marketplace for custom 3D fabrication |
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 |
---|
ANAS ABUZAINA ETC.: "3D Moving Object Reconstruction by Temporal Accumulation", 《2014 22ND INTERNATIONAL CONFERENCE ON PATTERN RECOGNITION》 * |
王崇倡 等: "圆柱面拟合方法研究", 《测绘工程》 * |
葛宝臻 等: "基于人体特征提取的模特体型尺寸自动测量方法", 《纺 织 学 报》 * |
Cited By (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108332713A (zh) * | 2018-01-16 | 2018-07-27 | 安徽慧视金瞳科技有限公司 | 一种绝缘子引脚倾斜的检测方法 |
CN108846888A (zh) * | 2018-04-23 | 2018-11-20 | 北京建筑大学 | 一种古木建筑构件精细尺寸信息自动化提取方法 |
CN108846888B (zh) * | 2018-04-23 | 2022-03-29 | 北京建筑大学 | 一种古木建筑构件精细尺寸信息自动化提取方法 |
CN111272091A (zh) * | 2020-03-09 | 2020-06-12 | 东南大学 | 一种基于三维激光扫描的焊接工字钢大变形自动识别方法 |
CN111272091B (zh) * | 2020-03-09 | 2021-03-23 | 东南大学 | 一种基于三维激光扫描的焊接工字钢大变形自动识别方法 |
CN111366113A (zh) * | 2020-03-31 | 2020-07-03 | 陕西梵翌琨机电科技有限公司 | 单晶硅棒测量模块及其晶体生长系统 |
CN111750808A (zh) * | 2020-06-28 | 2020-10-09 | 电子科技大学 | 基于激光雷达扫描仪的地表粗糙度获取方法 |
CN111750808B (zh) * | 2020-06-28 | 2021-09-14 | 电子科技大学 | 基于激光雷达扫描仪的地表粗糙度获取方法 |
CN113628166A (zh) * | 2021-07-12 | 2021-11-09 | 广东电网有限责任公司广州供电局 | 一种开关柜断路器室的几何尺寸检测方法及装置 |
CN115795579A (zh) * | 2022-12-23 | 2023-03-14 | 岭南师范学院 | 一种无特征复杂曲面误差分析的快速坐标对齐方法 |
Also Published As
Publication number | Publication date |
---|---|
CN106875439B (zh) | 2019-06-07 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106875439A (zh) | 基于三维点云模型的单晶硅棒外形尺寸测量方法 | |
CN102494657B (zh) | 一种曲面轮廓测量及检测的测头半径补偿方法 | |
CN103453849B (zh) | 多光学传感器协同的复杂曲面零件三维测量方法与系统 | |
CN103886593B (zh) | 一种基于三维点云曲面圆孔检测方法 | |
CN105354880A (zh) | 一种基于线激光扫描的喷砂机器人自动路径生成方法 | |
CN102001024A (zh) | 基于加工机床的自由曲面原位测量测量方法 | |
CN104392476A (zh) | 基于最小包围盒算法提取隧道三维轴线的方法 | |
CN103777570A (zh) | 基于nurbs曲面的加工误差快速检测补偿方法 | |
CN106248001A (zh) | 基于三维激光扫描仪的大型工艺管线对接平直度测量方法 | |
CN102853763A (zh) | 一种不规则表面文物体积的测量方法 | |
CN106769276B (zh) | 基于Dice相似度量的三维结构面代表性试样选取方法 | |
CN102506805B (zh) | 一种基于支持向量分类的多测点平面度评定方法 | |
CN103308021B (zh) | 一种测量工件圆度误差的方法 | |
CN104217083A (zh) | 一种基于多尺度分形函数的反射面天线面板建模方法 | |
CN103438828A (zh) | 一种激光检测螺杆转子端截形的方法 | |
CN106568365B (zh) | 一种球面孔系复合位置度误差的检测与评定方法 | |
CN103213068B (zh) | 超精密气囊抛光技术中工件边缘去除函数的测量方法 | |
CN101847262A (zh) | 一种快速三维点云搜索匹配方法 | |
CN103292654B (zh) | 一种计算圆柱体零件作用尺寸的方法 | |
CN105180886B (zh) | 一种测量冷轧薄钢板应变分布的方法 | |
CN106570338A (zh) | 结构面粗糙度尺寸效应研究中轮廓线采样精度确定方法 | |
CN114608461A (zh) | 一种非均匀壁厚零件的激光扫描测量方法 | |
CN112131682B (zh) | 一种叶片辊轧模具型面迭代优化的方法 | |
CN106643629A (zh) | 一种管状结构内表面粗糙度测量计算方法 | |
CN113418495A (zh) | 一种风洞喷管同轴度测量方法 |
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 | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20190607 Termination date: 20210220 |