CN102853763A - 一种不规则表面文物体积的测量方法 - Google Patents
一种不规则表面文物体积的测量方法 Download PDFInfo
- Publication number
- CN102853763A CN102853763A CN2012102812705A CN201210281270A CN102853763A CN 102853763 A CN102853763 A CN 102853763A CN 2012102812705 A CN2012102812705 A CN 2012102812705A CN 201210281270 A CN201210281270 A CN 201210281270A CN 102853763 A CN102853763 A CN 102853763A
- Authority
- CN
- China
- Prior art keywords
- grid model
- triangular grid
- volume
- historical relic
- cloud data
- 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.)
- Pending
Links
Images
Landscapes
- Length Measuring Devices By Optical Means (AREA)
Abstract
本发明公开了一种不规则表面文物体积的测量方法,包括如下步骤:获取待测文物表面三维空间点云数据信息;根据上述获取的三维空间点云数据信息构建文物的三角格网模型;对上述三角格网模型进行仿射变换,以将该三角格网模型的中心平移至坐标原点;计算该仿射变换后的三角格网模型中所有三角形与坐标原点构成的四面体的体积的代数和,该代数和即为所述待测文物的体积。本发明所述的不规则表面文物体积的测量方法中,通过三维扫描仪获取文物表面的三维空间点云数据信息,并通过这些点云数据信息构建三角格网模型,进而基于该三角格网模型采用改进的有向体积法计算文物体积,能够提高文物体积的测量精度,且整个测量过程操作简便、速度快、效率高。
Description
技术领域
本发明涉及一种文物体积测量方法,尤其涉及一种不规则表面文物体积的测量方法。
背景技术
随着岁月的流逝和人类活动的影响,人类留下来的珍贵而丰富的自然、文化遗产在逐渐的损毁、破坏和消失。如何利用先进的技术手段来保护、修复这些宝贵文物,成为迫在眉睫的问题。
在文物修复中经常需要估算材料的用量,然而一些文物结构非常复杂,表面呈不规则形状,传统测量方法无法直接求得精确体积以估算材料用量,通常是采用把体积参数转换成质量和密度的关系而进行间接测量。但间接测量方式易受条件、经验的约束,且测量速度慢、精度低、自动化程度不高。
发明内容
本发明针对现有技术的弊端,提供一种不规则表面文物体积的测量方法。
本发明所述的不规则表面文物体积的测量方法,包括如下步骤:
步骤一,获取待测文物表面三维空间点云数据信息;
步骤二,根据上述获取的三维空间点云数据信息构建文物的三角格网模型;
步骤三,对上述三角格网模型进行仿射变换,以将该三角格网模型的中心平移至坐标原点;
步骤四,计算该仿射变换后的三角格网模型中所有三角形与坐标原点构成的四面体的体积的代数和,该代数和即为所述待测文物的体积。
本发明所述的不规则表面文物体积的测量方法的步骤一中,对获取的三维空间点云数据信息中的杂点信息和噪点信息进行剔除。
本发明所述的不规则表面文物体积的测量方法的步骤一中,对获取的三维空间点云数据信息进行拼合合并。
本发明所述的不规则表面文物体积的测量方法的步骤三中,确定三角格网模型中心的步骤如下:
获取三角格网模型顶点坐标在X轴、Y轴、Z轴方向上的最小值与最大值,并标记为xmin、xmax,ymin、ymax,zmin、zmax;
以三维坐标值为xmin、ymin、zmin的点和三维坐标值为xmax、ymax、z max的点为对角顶点,构建包围住三角格网模型的包围盒,所述包围盒的中心即为该三角格网模型的中心。
本发明所述的不规则表面文物体积的测量方法的步骤四中,将所述仿射变换后的三角格网模型中的任意三角形ΔABC的一个顶点A至坐标原点O方向的单位向量标记为N,将该任意三角形ΔABC所在平面一侧的单位法向量标记为NABC,所述任意三角形顶点A、B、C按逆时针排列时与N成右手系,则当向量N与向量NABC的内积大于0时,该任意三角形与坐标原点构成的四面体的体积为正数;若向量N与向量NABC的内积小于0时,则该任意三角形与坐标原点构成的四面体的体积为负数。
本发明所述的不规则表面文物体积的测量方法中,通过三维扫描仪获取文物表面的三维空间点云数据信息,并通过这些点云数据信息构建三角格网模型,进而基于该三角格网模型采用改进的有向体积法计算文物体积,能够提高文物体积的测量精度,且整个测量过程操作简便、速度快、效率高。
附图说明
图1为本发明所述不规则表面文物体积的测量方法的流程示意图;
图2所示为格网模型中任意一个三角形与坐标原点所构成的四面体的示意图。
具体实施方式
下面结合附图对本发明做进一步的详细说明,以令本领域技术人员参照说明书文字能够据以实施。
由于四面体是三维空间最简单的不规则体单元,是进行三维几何拓扑关系描述的最基本的元素之一,因此,本发明中利用计算四面体有向体积法来计算三角格网模型体积,其实质是利用三维激光扫描仪对文物进行扫描获取其表面几何信息点云,再经过对点云数据预处理并通过拓扑重建得到文物的三角格网模型,从而来计算表面复杂的文物体积。
如图1所示,本发明所述的不规则表面文物体积的测量方法,包括如下步骤:
步骤101,获取待测文物表面三维空间点云数据信息。
文物表面三维信息数据的获取也可称之为实物表面数字化,通过特定的测量设备和测量方法获取实物表面点的三维坐标,从而获取大量的数据点信息,这些数据点称之为点云。
本实施例中,是利用便携式关节臂测量机(CIM CORE Infinite 2.0),其三维空间坐标测量精度为0.045mm,获取文物表面高精度三维空间点云。
由于扫描实体模型的过程中难免会受到环境光和载物装置等干扰,因此会出现杂点信息与噪点信息等,故而需要进行剔除。同时,为减少信息处理量,还需要精减数据,例如,对于多视图扫描文件,需要对这些扫描文件中的三维空间点云数据信息进行拼合合并处理等。
步骤102,根据上述获取的三维空间点云数据信息构建文物的三角格网模型。
由于点云模型实质上是采用大量具有三维坐标的点来描述被测物体的表面形态,是无数个离散点的集合,点与点之间是没有拓扑关系的,这些点云之间并没有构成文物的实际表面,所以要得到文物有拓扑关系的真实表面,有限元算法就是确定两点之间空间关系的基本算法,常用的建模方法有采用基本格网构造法和三角格网构造法等。
本步骤中,采用三角格网构造法构建文物的三角格网模型。所述三角格网模型是一种重要的三维形体表示方法,从大量的点云数据中重建的三维形体转化为三角格网模型进行存储和管理,便于采用不同的方法进行后续处理和分析,便于提取不同实际应用所需要的重要信息。
步骤103,对上述三角格网模型进行仿射变换,以将该三角格网模型的中心平移至坐标原点。
按照步骤102构建得到的三角格网模型,其表面极其复杂且不规则,为了提高体积计算精度,本步骤中,在对文物体积进行计算之前,还需将模型数据进行预处理,其目的是对原始三角格网模型进行仿射变换以将该模型的中心平移到坐标原点。
本步骤中,首先获取三角格网模型顶点坐标在X轴、Y轴、Z轴方向上的最小值与最大值,并标记为xmin、xmax,ymin、ymax,zmin、zmax;
以三维坐标值为xmin、ymin、zmin的点和三维坐标值为xmax、ymax、z max的点为对角顶点,构建包围住三角格网模型的包围盒,所述包围盒的中心即为该三角格网模型的中心。
之后,将上述确定的三角格网模型的中心平移到坐标原点。在此基础上再进行体积计算,能够提高体积计算的精度。
步骤104,计算该仿射变换后的三角格网模型中所有三角形与坐标原点构成的四面体的体积的代数和,该代数和即为所述待测文物的体积。
对于任一给定的网格模型,只有在它的网格面片围成有限封闭空间的情况下,求其体积才有意义。在体积计算之前确保模型的每个面片都是三角形,三角形的法向量与三角形的三个顶点编号顺序成右手规则,即所有三角形面片的法向一致,均指向体外。本发明中涉及的网格模型就是指这类模型,并且满足以下条件:(1)网格上的每个面片都是平面凸多边形。(2)网格上的每条边由且只由2个面片共有并且共用网格上任一顶点的各面片在该顶点处形成的曲面与单位圆拓扑同构。
在前述文物三角格网模型构建好并仿射变换后,即可基于该格网模型计算文物体积。
如图2所示,设格网模型中的任意一个三角形的三个顶点为A(x1,y 1,z1),B(x2,y2,z2),C(x3,y3,z3),则以ΔABC为底(按逆时针顺序排列),以坐标原点O(0,0,0)为顶点的四面体体积可以表示为如下形式:
本发明中,将所述仿射变换后的三角格网模型中的三角形ΔABC的一个 顶点A至坐标原点O方向的单位向量标记为N,将该任意三角形ΔABC所在平面一侧的单位法向量标记为NABC,所述任意三角形顶点A、B、C按逆时针排列时与N成右手系,则当向量N与向量NABC的内积大于0时,该任意三角形与坐标原点构成的四面体的体积为正数;若向量N与向量NABC的内积小于0时,则该任意三角形与坐标原点构成的四面体的体积为负数。
设三角格网模型表面由n个分片三角形所组成,其外侧面第i个三角形的三个顶点为Ai1(xi1,yi1,zi1),Bi(xi2,yi2,zi2),Ci(x i3,yi3,zi3),三个顶点的顺序按曲面外侧逆时针方向排列,则三维模型总体积计算公式可以用如下公式表示:
其中,Vi为第i个三角形与坐标原点构成的带符号体积,其正负符号同样根据原点O和三角形中一顶点构成的向量与三角形的法向量的内积的正负号确定。
本发明所述的不规则表面文物体积的测量方法中,通过三维扫描仪获取文物表面的三维空间点云数据信息,并通过这些点云数据信息构建三角格网模型,进而基于该三角格网模型采用改进的有向体积法计算文物体积,能够提高文物体积的测量精度,且整个测量过程操作简便、速度快、效率高。
尽管本发明的实施方案已公开如上,但其并不仅仅限于说明书和实施方式中所列运用,它完全可以被适用于各种适合本发明的领域,对于熟悉本领域的人员而言,可容易地实现另外的修改,因此在不背离权利要求及等同范围所限定的一般概念下,本发明并不限于特定的细节和这里示出与描述的图例。
Claims (5)
1.一种不规则表面文物体积的测量方法,其特征在于,包括如下步骤:
步骤一,获取待测文物表面三维空间点云数据信息;
步骤二,根据上述获取的三维空间点云数据信息构建文物的三角格网模型;
步骤三,对上述三角格网模型进行仿射变换,以将该三角格网模型的中心平移至坐标原点;
步骤四,计算该仿射变换后的三角格网模型中所有三角形与坐标原点构成的四面体的体积的代数和,该代数和即为所述待测文物的体积。
2.如权利要求1所述的测量方法,其特征在于,所述步骤一中,对获取的三维空间点云数据信息中的杂点信息和噪点信息进行剔除。
3.如权利要求1所述的测量方法,其特征在于,所述步骤一中,对获取的三维空间点云数据信息进行拼合合并。
4.如权利要求1所述的测量方法,其特征在于,所述步骤三中,确定三角格网模型中心的步骤如下:
获取三角格网模型顶点坐标在X轴、Y轴、Z轴方向上的最小值与最大值,并标记为xmin、xmax,ymin、ymax,zmin、zmax;
以三维坐标值为xmin、ymin、zmin的点和三维坐标值为xmax、ymax、zmax的点为对角顶点,构建包围住三角格网模型的包围盒,所述包围盒的中心即为该三角格网模型的中心。
5.如权利要求1所述的测量方法,其特征在于,所述步骤四中,将所述仿射变换后的三角格网模型中的任意三角形△ABC的一个顶点A至坐标原点O方向的单位向量标记为N,将该任意三角形△ABC所在平面一侧的单位法向量标记为NABC,所述任意三角形顶点A、B、C按逆时针排列时与N成右手系,则当向量N与向量NABC的内积大于0时,该任意三角形与坐标原点构成的四面体的体积为正数;若向量N与向量NABC的内积小于0时,则该任意三角形与坐标原点构成的四面体的体积为负数。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2012102812705A CN102853763A (zh) | 2012-08-08 | 2012-08-08 | 一种不规则表面文物体积的测量方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2012102812705A CN102853763A (zh) | 2012-08-08 | 2012-08-08 | 一种不规则表面文物体积的测量方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN102853763A true CN102853763A (zh) | 2013-01-02 |
Family
ID=47400568
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN2012102812705A Pending CN102853763A (zh) | 2012-08-08 | 2012-08-08 | 一种不规则表面文物体积的测量方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102853763A (zh) |
Cited By (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103983188A (zh) * | 2014-05-06 | 2014-08-13 | 华东师范大学 | 一种运动物体的体积测量方法及其测量装置 |
CN104457574A (zh) * | 2014-12-11 | 2015-03-25 | 天津大学 | 一种非接触式测量不规则物体体积的装置和测量方法 |
CN106643494A (zh) * | 2016-12-22 | 2017-05-10 | 上海华测导航技术股份有限公司 | 一种矿山堆料体积测量方法及系统 |
CN109507723A (zh) * | 2017-09-15 | 2019-03-22 | 中国石油化工股份有限公司 | 微地震压裂裂缝模型的压裂体积计算方法及系统 |
CN109766661A (zh) * | 2019-01-28 | 2019-05-17 | 中国矿业大学 | 一种异形薄壁建筑结构加固方法 |
CN111156922A (zh) * | 2019-12-23 | 2020-05-15 | 苏州迭慧智能科技有限公司 | 一种利用外形轮廓测量的方法 |
CN112258650A (zh) * | 2020-09-21 | 2021-01-22 | 北京科技大学 | 一种膏体充填进度实时测量与可视化方法及系统 |
CN113048886A (zh) * | 2021-05-31 | 2021-06-29 | 山东捷瑞数字科技股份有限公司 | 计量工件不规则体尺寸的测量方法及设备 |
CN113251926A (zh) * | 2021-06-04 | 2021-08-13 | 山东捷瑞数字科技股份有限公司 | 一种不规则物体的尺寸测量方法及测量装置 |
CN113418445A (zh) * | 2021-05-28 | 2021-09-21 | 中北大学 | 基于激光点云的固体火箭发动机内腔体积测量方法及装置 |
CN113947630A (zh) * | 2021-10-20 | 2022-01-18 | 吉旗物联科技(上海)有限公司 | 一种物体的体积估算方法、装置及存储介质 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2001067495A (ja) * | 1999-08-25 | 2001-03-16 | Ricoh Co Ltd | 三角形メッシュの生成方法、装置及び記録媒体 |
JP2001351119A (ja) * | 2000-06-06 | 2001-12-21 | Ricoh Co Ltd | 三角形メッシュ生成方法及び該方法を実行するための記録媒体 |
-
2012
- 2012-08-08 CN CN2012102812705A patent/CN102853763A/zh active Pending
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2001067495A (ja) * | 1999-08-25 | 2001-03-16 | Ricoh Co Ltd | 三角形メッシュの生成方法、装置及び記録媒体 |
JP2001351119A (ja) * | 2000-06-06 | 2001-12-21 | Ricoh Co Ltd | 三角形メッシュ生成方法及び該方法を実行するための記録媒体 |
Non-Patent Citations (2)
Title |
---|
XIAOQING ZHANG ET AL.: "Calculation of Areas of Cultural Relics Surface Defects Based on the Triangular Mesh Model", 《ADVANCED MATERIALS RESEARCH》 * |
张小青 等: "基于四面体的不规则表面文物体积计算", 《测绘通报》 * |
Cited By (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103983188A (zh) * | 2014-05-06 | 2014-08-13 | 华东师范大学 | 一种运动物体的体积测量方法及其测量装置 |
CN103983188B (zh) * | 2014-05-06 | 2017-01-25 | 华东师范大学 | 一种运动物体的体积测量方法及其测量装置 |
CN104457574A (zh) * | 2014-12-11 | 2015-03-25 | 天津大学 | 一种非接触式测量不规则物体体积的装置和测量方法 |
CN106643494A (zh) * | 2016-12-22 | 2017-05-10 | 上海华测导航技术股份有限公司 | 一种矿山堆料体积测量方法及系统 |
CN109507723A (zh) * | 2017-09-15 | 2019-03-22 | 中国石油化工股份有限公司 | 微地震压裂裂缝模型的压裂体积计算方法及系统 |
CN109766661A (zh) * | 2019-01-28 | 2019-05-17 | 中国矿业大学 | 一种异形薄壁建筑结构加固方法 |
CN111156922A (zh) * | 2019-12-23 | 2020-05-15 | 苏州迭慧智能科技有限公司 | 一种利用外形轮廓测量的方法 |
CN112258650A (zh) * | 2020-09-21 | 2021-01-22 | 北京科技大学 | 一种膏体充填进度实时测量与可视化方法及系统 |
CN112258650B (zh) * | 2020-09-21 | 2024-03-19 | 北京科技大学 | 一种膏体充填进度实时测量与可视化方法及系统 |
CN113418445A (zh) * | 2021-05-28 | 2021-09-21 | 中北大学 | 基于激光点云的固体火箭发动机内腔体积测量方法及装置 |
CN113048886A (zh) * | 2021-05-31 | 2021-06-29 | 山东捷瑞数字科技股份有限公司 | 计量工件不规则体尺寸的测量方法及设备 |
CN113048886B (zh) * | 2021-05-31 | 2021-08-17 | 山东捷瑞数字科技股份有限公司 | 计量工件不规则体尺寸的测量方法及设备 |
CN113251926A (zh) * | 2021-06-04 | 2021-08-13 | 山东捷瑞数字科技股份有限公司 | 一种不规则物体的尺寸测量方法及测量装置 |
CN113947630A (zh) * | 2021-10-20 | 2022-01-18 | 吉旗物联科技(上海)有限公司 | 一种物体的体积估算方法、装置及存储介质 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102853763A (zh) | 一种不规则表面文物体积的测量方法 | |
CN110516388B (zh) | 基于调和映射的曲面离散点云模型环切刀轨生成方法 | |
CN102194253B (zh) | 一种面向三维地质层面结构的四面体网格生成方法 | |
CN106248035A (zh) | 一种基于点云模型的面轮廓度精确检测的方法及系统 | |
CN104282040A (zh) | 一种用于三维实体模型重建的有限元前处理方法 | |
CN103942837B (zh) | 基于序列线性规划的叶片点云模型截面曲线直接构造方法 | |
CN102853803B (zh) | 一种文物破损面积的测定方法 | |
CN103236043B (zh) | 一种植物器官点云修复方法 | |
CN107622530B (zh) | 一种高效鲁棒的三角网切割方法 | |
Han et al. | An adaptive GRP scheme for compressible fluid flows | |
CN103106632A (zh) | 一种基于均值漂移的不同精度三维点云数据的融合方法 | |
CN106482700B (zh) | 一种草图直接成图的数字化房产面积测量方法 | |
CN103759677A (zh) | 基于三角面积比法度量岩体结构面三维粗糙度的方法 | |
CN107545602B (zh) | 基于LiDAR点云的空间拓扑关系约束下的建筑物建模方法 | |
CN103914879A (zh) | 一种在抛物线方程中由三角面元数据生成立方网格数据的方法 | |
CN106097433A (zh) | 物体工业与影像模型的叠加方法及系统 | |
CN107886569A (zh) | 一种基于离散李导数的测度可控的曲面参数化方法及系统 | |
Song et al. | A comparison study of algorithms for surface normal determination based on point cloud data | |
CN103191958A (zh) | 一种板料成形回弹的光学检测方法 | |
CN104949628B (zh) | 基于二维正交曲率的柔性板状结构复杂形态重建方法 | |
Wang et al. | Feature-based solid model reconstruction | |
Stanko et al. | Surfacing curve networks with normal control | |
CN101546351B (zh) | 一种叶轮的变复杂度形状优化的几何参数化建模方法 | |
CN109767492A (zh) | 一种变电站三维模型的间距计算方法 | |
CN101877145A (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 | ||
C12 | Rejection of a patent application after its publication | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20130102 |