CN102592284B - 零件表面形貌三维高密度点云数据转化为灰度图像的方法 - Google Patents

零件表面形貌三维高密度点云数据转化为灰度图像的方法 Download PDF

Info

Publication number
CN102592284B
CN102592284B CN201210046615.9A CN201210046615A CN102592284B CN 102592284 B CN102592284 B CN 102592284B CN 201210046615 A CN201210046615 A CN 201210046615A CN 102592284 B CN102592284 B CN 102592284B
Authority
CN
China
Prior art keywords
matrix
cloud data
piece surface
dimensional high
dimensional
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
Application number
CN201210046615.9A
Other languages
English (en)
Other versions
CN102592284A (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.)
Shanghai Jiaotong University Chi bang technology Co., Ltd.
Original Assignee
Shanghai Jiaotong 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 Shanghai Jiaotong University filed Critical Shanghai Jiaotong University
Priority to CN201210046615.9A priority Critical patent/CN102592284B/zh
Publication of CN102592284A publication Critical patent/CN102592284A/zh
Application granted granted Critical
Publication of CN102592284B publication Critical patent/CN102592284B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Image Processing (AREA)
  • Image Analysis (AREA)

Abstract

零件表面形貌三维高密度点云数据转化为灰度图像的方法。属于零件表面形貌检测和点云数据处理技术领域。本发明方法包括:首先,通过三维高分辨率表面形貌测量技术测量零件表面得到三维高密度点云数据。然后,利用MATLAB软件读入三维高密度点云数据,通过插值命令将点云数据的Z坐标值转换为栅格各点的值,存储为二维矩阵I,并剔除矩阵I中不属于零件表面的元素。最后将矩阵I转换为灰度图像。本发明生成的灰度图像,能够覆盖整个零件表面,有效的反映零件表面整体形貌和加工纹理,改进了传统图像采集系统需通过局部图像匹配、拼接获取零件表面全貌的缺点和三维高分辨率表面形貌测量技术反映加工纹理特征方面的不足。

Description

零件表面形貌三维高密度点云数据转化为灰度图像的方法
技术领域:
本发明涉及一种零件表面形貌三维高密度点云数据转化为灰度图像的方法,通过将三维高密度点云数据中反映零件表面形貌的高度值映射为像素的灰度值,生成能够反映零件整体表面形貌特别是纹理特征的灰度图像。属于零件表面形貌检测和点云数据处理技术领域。
背景技术:
随着对精密零件表面加工要求的不断提高,作为一种重要的产品质量控制手段,机械加工零件表面形貌检测技术被广泛应用于机械表面的加工过程。铣削、磨削、刨削等工艺加工出的零件表面,由于加工过程刀具行程的变化、材料本身特性以及振动等原因,会在机械加工零件表面形成不同的纹理。对表面加工纹理的分析研究,一方面可以对零件的加工质量及机床本身的精度进行综合评价,另一方面可以用于表面质量缺陷分析及内在根源追溯。
传统的零件加工表面纹理图像大多由图像采集系统直接获得。已有技术中,黎明在论文“机械加工零件表面纹理缺陷检测”(《中国图象图形学报》2004年第9卷第3期,318-322页)中提到的零件表面纹理图像采集系统由荧光光源、显微镜和CCD摄像机组成。该零件表面纹理图像采集系统,由于测量精度的要求以及受视觉范围的限制,只能取得零件表面的局部形貌,无法获取零件表面的整体形貌,必须通过局部图像的匹配、拼接来获取零件表面全貌。另外,由于生成的零件表面形貌图像清晰度较低,损失了零件表面形貌特征的部分细节。
专利号为WO2008070746名称为“SYSTEM AND METHOD FOR SHIFTING PHASE IN AMULTI-WAVELENGTH INTERFEROMETRIC IMAGING SYSTEM”的美国专利,提出的三维高分辨率表面形貌测量技术能够对零件表面形貌进行整体检测,能够生成反映零件整体表面形貌的三维高密度点云数据,该数据为X、Y、Z三维坐标格式,以及彩色编码的三维表面形貌图像。该三维表面形貌图像虽然能够很好的反映零件表面各点的高度,但是对于反映零件加工表面留下痕迹的深浅、疏密、形状,特别是纹理特征的表现则不太理想,不能有效地对零件表面纹理进行描述和识别。
发明内容:
为了克服三维高分辨率表面形貌测量技术生成的彩色编码三维表面形貌图像,在显示零件表面加工纹理方面的不足,本发明提出了一种将零件表面形貌三维高密度点云数据转换为灰度图像的方法。由于三维高密度点云数据测点采样密度高,可达到每平方毫米40个测点,故可以将三维高密度点云数据中的点转化为图像中的像素,而不损失零件表面形貌的细节信息。该方法生成的灰度图像既能反映零件表面的整体形貌,又可以有效的体现零件表面的加工纹理特征。
本发明是按照下述技术方案实现的。本发明方法包括下述步骤:
第一步:通过采用三维高分辨率表面形貌测量技术对零件表面进行测量,得到三维高密度点云数据,该数据以X、Y、Z三维坐标格式显示,并得到彩色编码的三维表面形貌图像。
第二步:读入三维高密度点云数据。在计算机上利用软件分别将三维高密度点云数据各点的X、Y、Z三维坐标读入,存储为N行3列的矩阵A。其中行数N为点云数量;第一、二列存储各点的X、Y坐标值,表示各点在零件表面的位置;第三列存储各点的Z坐标值,表示各点相对于基准面的高度。
第三步:对读入的三维高密度点云数据进行插值。利用矩阵A第一、二列存储的各点X、Y的坐标值,生成覆盖零件表面的X、Y方向的矩形域栅格,通过软件插值命令将矩阵A中第三列存储的Z坐标值转换为栅格各点的值,并存储为二维矩阵I。
第四步:进行边界判断。利用矩阵A中前两列存储的X、Y坐标值生成零件表面轮廓,以判断矩阵I中各元素是否属于零件表面上的测点,通过矩阵I中元素位置计算得到该矩阵I中的测点在零件表面的实际位置,与零件表面轮廓比较,如果比较结果小于设定的边界阈值,则判定该矩阵I中的测点属于零件表面测点,并保留该矩阵I中的测点数值,如果大于设定的阈值,则判定该矩阵I中的测点不属于零件表面测点,将该矩阵I中的测点的值设为空。边界阈值大于等于栅格的步长。
第五步:将矩阵I转换为灰度图像。将二维矩阵I中每个元素的值转化为灰度值,并将不属于零件表面的点设置为背景色,生成能够有效反映零件表面整体形貌和加工纹理的灰度图像,最后对灰度图像进行显示和处理。
本发明的有益效果:
本发明涉及一种零件表面形貌三维高密度点云数据转化为灰度图像的方法。本发明生成的灰度图像,能够覆盖整个零件表面,可以有效的反映零件表面的整体形貌和加工纹理,改进了传统图像采集系统需通过局部图像的匹配、拼接来获取零件表面全貌的缺点,以及三维高分辨率表面形貌测量技术生成的三维表面形貌图像不足以反映零件加工表面纹理特征的缺陷。同时,由于三维高密度点云数据密集,保证了零件表面灰度图像的清晰度,从而为后续基于图像处理的零件形貌评价、模式匹配及质量控制方法提供了坚实的基础,具有重要的工程实用价值。
附图说明:
图1是三维高分辨率表面形貌测量技术测量发动机缸体顶面生成的三维表面形貌图像
图2是本发明三维高密度点云数据各点的X、Y坐标值生成的发动机缸体顶面轮廓图像
图3是本发明三维高密度点云数据转化为的发动机缸体顶面形貌灰度图像
具体实施方式:
下面以某发动机厂生产的发动机缸体顶面的测量为例,结合附图对本发明的具体实施作进一步描述。
如图1、图2、图3所示,本发明方法包括下述步骤:
第一步:通过采用三维高分辨率表面形貌测量技术对发动机缸体顶面进行测量,得到三维高密度点云数据,该数据以X、Y、Z三维坐标格式显示,并得到彩色编码的三维表面形貌图,如图1所示。其X、Y方向分辨率均为150μm,Z方向测量精度为1μm,点云数据密度可达到每平方毫米40个测点,测点总数可达4百万点。
第二步:读入三维高密度点云数据。在计算机上利用MATLAB软件分别将三维高密度点云数据各点的X、Y、Z三维坐标读入,存储为N行3列的矩阵A。其中行数N为点云数量;第一、二列存储各点的X、Y坐标值,表示各点在发动机缸体顶面的位置;第三列存储各点的Z坐标值,表示各点相对于基准面的高度。
第三步:对读入的三维高密度点云数据进行插值。利用矩阵A第一、二列存储的X、Y坐标值,生成覆盖发动机缸体顶面的X、Y方向的矩形域栅格,其中栅格的极值为三维高密度点云数据中X、Y方向的极值,栅格步长根据X、Y方向分辨率设置,应大于等于150μm。通过插值命令将矩阵A中第三列存储的Z坐标值转换为栅格各点的值,并存储为二维矩阵I。
第四步:进行边界判断,利用矩阵A中前两列存储的X、Y坐标值生成发动机缸体顶面轮廓,如图2所示,并选取二维矩阵I中属于发动机缸体顶面的测点数据。通过矩阵I中元素位置计算得到该矩阵I中的测点在发动机缸体顶面的实际位置,与发动机缸体顶面轮廓比较,如果比较结果小于设定的边界阈值,则判定该矩阵I中的测点属于发动机缸体顶面测点,并保留该矩阵I中的测点数值,如果大于设定的边界阈值,则判定该矩阵I中的测点不属于发动机缸体顶面测点,将该矩阵I中的测点的值设为空。边界阈值大于等于栅格的步长。
第五步:将二维矩阵I转换为灰度图像。利用MATLAB软件将二维矩阵I中每个元素的值转化为像素的灰度,灰度的取值范围为[01]或[0255],并将不属于发动机缸体顶面的点设置为背景色,生成能够有效反映发动机缸体顶面整体形貌和加工纹理的灰度图像,最后对灰度图像进行显示和处理,结果如图3所示。

Claims (5)

1.一种零件表面形貌三维高密度点云数据转化为灰度图像的方法,其特征在于,该方法包括下述步骤:
第一步:通过采用三维高分辨率表面形貌测量技术对零件表面进行测量,得到三维高密度点云数据,该数据以X、Y、Z三维坐标格式显示,并得到彩色编码的三维表面形貌图像;
第二步:读入三维高密度点云数据,在计算机上利用软件分别将三维高密度点云数据各点的X、Y、Z三维坐标读入,存储为N行3列的矩阵A,其中行数N为点云数量;第一、二列存储各点的X、Y坐标值,表示各点在零件表面的位置;第三列存储各点的Z坐标值,表示各点相对于基准面的高度;
第三步:对读入的三维高密度点云数据进行插值,利用矩阵A第一、二列存储的各点X、Y的坐标值,生成覆盖零件表面的X、Y方向的矩形域栅格,通过软件插值命令将矩阵A中第三列存储的Z坐标值转换为栅格各点的值,并存储为二维矩阵I;
第四步:进行边界判断,利用矩阵A中前两列存储的X、Y坐标值生成零件表面轮廓,以判断矩阵I中各元素是否属于零件表面上的测点,通过矩阵I中元素位置计算得到该矩阵I中的测点在零件表面的实际位置,与零件表面轮廓比较,如果比较结果小于设定的边界阈值,则判定该矩阵I中的测点属于零件表面测点,并保留该矩阵I中的测点数值,如果大于设定的阈值,则判定该矩阵I中的测点不属于零件表面测点,将该矩阵I中的测点的值设为空,边界阈值大于等于栅格的步长;
第五步:将矩阵I转换为灰度图像,将二维矩阵I中每个元素的值转化为灰度值,并将不属于零件表面的点设置为背景色,生成能够有效反映零件表面整体形貌和加工纹理的灰度图像,最后对灰度图像进行显示和处理。
2.根据权利要求1所述的零件表面形貌三维高密度点云数据转化为灰度图像的方法,其特征是所说的步骤一中,三维高分辨率表面形貌测量技术的X、Y方向分辨率均为150μm,Z方向测量精度为1μm,点云数据密度可达到每平方毫米40个测点,测点总数可达4百万点。
3.根据权利要求1所述的零件表面形貌三维高密度点云数据转化为灰度图像的方法,其特征是所说的步骤二中,读入三维高密度点云数据所用的软件为MATLAB。
4.根据权利要求1所述的零件表面形貌三维高密度点云数据转化为灰度图像的方法,其特征是所说的步骤三中,栅格的极值为三维高密度点云数据中X、Y方向的极值,栅格步长根据X、Y方向分辨率设置,应大于等于150μm,插值命令所用的软件为MATLAB。
5.根据权利要求1所述的零件表面形貌三维高密度点云数据转化为灰度图像的方法,其特征是所说的步骤五中,是用MATLAB软件将二维矩阵I中每个元素的值转化为像素的灰度,灰度的取值范围为[01]或[0255]。
CN201210046615.9A 2012-02-27 2012-02-27 零件表面形貌三维高密度点云数据转化为灰度图像的方法 Active CN102592284B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201210046615.9A CN102592284B (zh) 2012-02-27 2012-02-27 零件表面形貌三维高密度点云数据转化为灰度图像的方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201210046615.9A CN102592284B (zh) 2012-02-27 2012-02-27 零件表面形貌三维高密度点云数据转化为灰度图像的方法

Publications (2)

Publication Number Publication Date
CN102592284A CN102592284A (zh) 2012-07-18
CN102592284B true CN102592284B (zh) 2014-08-06

Family

ID=46480874

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201210046615.9A Active CN102592284B (zh) 2012-02-27 2012-02-27 零件表面形貌三维高密度点云数据转化为灰度图像的方法

Country Status (1)

Country Link
CN (1) CN102592284B (zh)

Families Citing this family (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103544694A (zh) * 2013-09-22 2014-01-29 上海交通大学 基于高清晰测量灰度图像的零件表面三维形貌评价方法
CN104316541A (zh) * 2014-11-12 2015-01-28 京东方科技集团股份有限公司 缺陷检测装置及偏光片贴附设备
CN105550682B (zh) * 2015-11-26 2019-10-22 北京市计算中心 钟鼎碑刻拓印方法
CN105737802B (zh) * 2016-01-26 2017-11-17 中国科学院水利部成都山地灾害与环境研究所 基于体感摄影技术的堆积剖面空间结构信息分析方法
CN106023303B (zh) * 2016-05-06 2018-10-26 西安电子科技大学 一种基于轮廓有效性提高三维重建点云稠密程度的方法
CN106949848B (zh) * 2017-03-15 2019-11-29 深圳市深视智能科技有限公司 一种高精度激光3d轮廓手机结构件检测方法
CN108168495A (zh) * 2018-01-10 2018-06-15 天津大学 一种材料表面采样参数的确定方法
CN108682012B (zh) * 2018-05-15 2021-07-27 佛山市南海区广工大数控装备协同创新研究院 一种基于线扫激光的3d曲面玻璃表面平整性缺陷检测方法
US11440119B2 (en) * 2018-10-12 2022-09-13 Teradyne, Inc. System and method for weld path generation
CN111168136B (zh) * 2020-01-16 2021-08-06 上海交通大学 基于高清晰测量的面铣削表面刀纹分割方法及系统
CN112419360B (zh) * 2020-11-16 2023-02-21 北京理工大学 一种基于立体成像的背景去除和目标图像分割方法
CN113095195B (zh) * 2021-04-03 2023-04-07 西北工业大学 一种基于表面形貌自特征的零件唯一识别方法

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102289784A (zh) * 2011-07-29 2011-12-21 中山大学 一种基于点云模型的数字几何图像简化方法

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2008070746A2 (en) * 2006-12-05 2008-06-12 Coherix, Inc. System and method for shifting phase in a multi-wavelength interferometric imaging system
US20110222757A1 (en) * 2010-03-10 2011-09-15 Gbo 3D Technology Pte. Ltd. Systems and methods for 2D image and spatial data capture for 3D stereo imaging

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102289784A (zh) * 2011-07-29 2011-12-21 中山大学 一种基于点云模型的数字几何图像简化方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
张明礼 等.点云曲面的复合几何图像表示及其应用.《浙江大学学报(工学版)》.2004, *

Also Published As

Publication number Publication date
CN102592284A (zh) 2012-07-18

Similar Documents

Publication Publication Date Title
CN102592284B (zh) 零件表面形貌三维高密度点云数据转化为灰度图像的方法
CN102003938B (zh) 大型高温锻件热态在位检测方法
CN111932671A (zh) 基于密集点云数据的三维实体模型重建方法
CN107514984A (zh) 一种基于光学显微的三维表面形貌测量方法和系统
CN104331699A (zh) 一种三维点云平面化快速搜索比对的方法
CN205352391U (zh) 基于地面三维激光扫描的实体测量系统
CN101986350A (zh) 基于单目结构光的三维建模方法
CN110070488B (zh) 一种基于卷积神经网络的多角度遥感影像森林高度提取方法
CN114419130A (zh) 一种基于图像特征和三维点云技术的散料体积测量方法
CN201858967U (zh) 一种工业激光三角测量与高速图像存储系统
CN108230442A (zh) 一种盾构隧道三维仿真方法
Cheng et al. Study on reverse engineering of historical architecture based on 3D laser scanner
CN110006374B (zh) 一种通过多维插值获取图像基准的cfrp粗糙度测量方法
TWI731604B (zh) 三維點雲資料處理方法
Grifoni et al. 3D multi-modal point clouds data fusion for metrological analysis and restoration assessment of a panel painting
Polo et al. Accurate 3D models in both geometry and texture: An archaeological application
CN112665528B (zh) 一种激光扫描三维成像的矫正方法
CN114170377A (zh) 一种基于bim的历史建筑特色构件族库构建方法及系统
Zhang et al. On the research of cultural relic restoration under reverse design
Jiang et al. Underwater 3D reconstruction based on laser line scanning
Li et al. Computer-assisted archaeological line drawing
Wan et al. Application of Structured-light 3D Reconstruction in Virtual Assembly
Faresin et al. Micrometer multiresolution laser scanning of a renaissance medallion
Patrucco et al. Digital Replicas of British Museum Artefacts
Russo Integrated Reverse Modeling techniques for the survey of complex shapes in Industrial Design

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
TR01 Transfer of patent right
TR01 Transfer of patent right

Effective date of registration: 20170831

Address after: 200030 Huashan Road, Shanghai, No. 1954, No.

Co-patentee after: Ling Weiguo

Patentee after: Jin Sun

Co-patentee after: Chen Ming

Co-patentee after: Chen Zhixiang

Co-patentee after: Luo Lei

Co-patentee after: Shi Lun

Address before: 200240 Dongchuan Road, Shanghai, No. 800, No.

Patentee before: Shanghai Jiao Tong University

TR01 Transfer of patent right
TR01 Transfer of patent right

Effective date of registration: 20180208

Address after: The new town of Pudong New Area Nanhui lake west two road 201306 Shanghai City No. 888 building C

Patentee after: Shanghai Zhi Bangyuan Intelligent Technology Co Ltd

Address before: 200030 Huashan Road, Shanghai, No. 1954, No.

Co-patentee before: Ling Weiguo

Patentee before: Jin Sun

Co-patentee before: Chen Ming

Co-patentee before: Chen Zhixiang

Co-patentee before: Luo Lei

Co-patentee before: Shi Lun

TR01 Transfer of patent right
TR01 Transfer of patent right

Effective date of registration: 20200116

Address after: 201306 Building 1, No.66, Feidu Road, Nanhui new town, Pudong New Area, Shanghai

Patentee after: Shanghai Jiaotong University Chi bang technology Co., Ltd.

Address before: The new town of Pudong New Area Nanhui lake west two road 201306 Shanghai City No. 888 building C

Patentee before: Shanghai Zhi Bangyuan Intelligent Technology Co Ltd