CN113963050A - 一种基于点云的土方量计算方法及系统 - Google Patents
一种基于点云的土方量计算方法及系统 Download PDFInfo
- Publication number
- CN113963050A CN113963050A CN202110920858.XA CN202110920858A CN113963050A CN 113963050 A CN113963050 A CN 113963050A CN 202110920858 A CN202110920858 A CN 202110920858A CN 113963050 A CN113963050 A CN 113963050A
- Authority
- CN
- China
- Prior art keywords
- elevation
- point cloud
- areas
- average
- polygon
- 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
- 238000000034 method Methods 0.000 title claims abstract description 36
- 238000004364 calculation method Methods 0.000 claims abstract description 34
- 238000013461 design Methods 0.000 claims abstract description 29
- 238000009412 basement excavation Methods 0.000 claims abstract description 8
- 238000000605 extraction Methods 0.000 claims description 6
- 238000004590 computer program Methods 0.000 claims description 3
- 238000004891 communication Methods 0.000 claims description 2
- 238000010276 construction Methods 0.000 description 6
- 238000010586 diagram Methods 0.000 description 5
- 230000009286 beneficial effect Effects 0.000 description 2
- 239000000284 extract Substances 0.000 description 2
- 230000003287 optical effect Effects 0.000 description 2
- 239000004575 stone Substances 0.000 description 2
- 238000012935 Averaging Methods 0.000 description 1
- 238000000205 computational method Methods 0.000 description 1
- 230000007123 defense Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000001556 precipitation Methods 0.000 description 1
- 230000011218 segmentation Effects 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
Images
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
- G06T7/62—Analysis of geometric attributes of area, perimeter, diameter or volume
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/08—Construction
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T17/00—Three dimensional [3D] modelling, e.g. data description of 3D objects
- G06T17/20—Finite element generation, e.g. wire-frame surface description, tesselation
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Geometry (AREA)
- Business, Economics & Management (AREA)
- Human Resources & Organizations (AREA)
- Strategic Management (AREA)
- Tourism & Hospitality (AREA)
- Primary Health Care (AREA)
- General Business, Economics & Management (AREA)
- Marketing (AREA)
- General Health & Medical Sciences (AREA)
- Economics (AREA)
- Health & Medical Sciences (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Computer Graphics (AREA)
- Software Systems (AREA)
- Processing Or Creating Images (AREA)
Abstract
本发明属于土方工程领域,具体涉及一种基于点云的土方量计算方法,包括:获取地形点云数据和设计平面,将点云数据划分为若干高程区域;分别计算若干高程区域内点云高度坐标的平均值得到该高程区域的平均高程;将设计平面划分为若干多边形;判断所述若干多边形的顶点所位于的高程区域,获取所述若干多边形的顶点所位于高程区域的平均高程作为该多边形顶点的平均高程;并通过多边形顶点的平均高程得到该多边形的计算高程;通过所述若干多边形的面积和计算高程得到土方量,并通过土方量的正负判断是挖方量或填方量。本发明在保证了计算精度的前提下,简化了土方的计算量;并且本发明的计算边界和设计边界一致,避免了边界误差。
Description
技术领域
本发明属于土方工程领域,具体涉及一种基于点云的土方量计算方法。
背景技术
土方工程是建筑工程施工中主要工程之一,包括一切土方的挖梆、填筑、 运输以及排水、降水等方面。土木工程中,土石方工程有:场地平整、路基开 挖、人防工程开挖、地坪填土,路基填筑以及基坑回填。要合理安排施工计划, 尽量不要安排在雨季,同时为了降低土石方工程施工费用,贯彻不占或少占农 田和可耕地并有利于改地造田的原则,要作出土石方的合理调配方案,统筹安 排。土方量的计算作为施工结算和工作量结算的重要依据,与施工的各个周期 都息息相关。
现有技术下,计算基坑设计挖方量往往采用三角网法或者方格网法。三角 网法运算量大,尤其是对于大数据量的点云而言,必须借助高性能计算机进行 长时间运算才能实现。方格网法虽然运算量小,但是由于其确定的边界是近似 边界,导致计算结果不准确。如公开号为:CN111595403A的中国专利文献,其 公开了一种基于点云测量技术的工程土方计量方法,该方法采用三角形格网或 泰森多边形格网建对点云进行划分,从而对土方量进行计算;但是其计算量巨 大,需要很高的算力进行求解;并且该方法由于设计提取的点云和测量点云均 单独建立格网,会导致两者边界不统一,带来误差。
发明内容
本发明的目的在于克服现有技术中所存在的上述不足,提供一种基于点云 的土方量计算方法系统及可读存储介质,通该方法可以减少计算量,并且消除 了现有的计算方式带来的边界误差。
为实现上述目的,本发明采用的技术方案包括以下各方面:
一种基于点云的土方量计算方法,包括如下步骤:
S1.获取地形点云数据和设计平面,将点云数据划分为若干高程区域;
分别计算若干高程区域内点云高度坐标的平均值得到平均高程;
S2.将设计平面划分为若干多边形,获取所述若干多边形的顶点所在的高程 区域,将高程区域的平均高程作为所述若干多边形的顶点的平均高程;并通过 多边形顶点的平均高程得到该多边形的计算高程;
S3.通过所述若干多边形的面积和计算高程得到土方量,并通过土方量的正负 判断是挖方量或填方量。
进一步是,所述步骤S1包括如下步骤:
S11.将点云区域划分为若干边长相同的正方形高程区域,正方形高程区域边 长d的计算公式为:
其中,Xmax表示点云中所有点X坐标的最大值,Xmin表示点云中所有点X 坐标的最小值,Ymax表示点云中所有点Y坐标的最大值,Ymin表示点云中所有 点Y坐标的最小值;
S12.计算每一个正方形高程区域内点云的平均高程,
其中,Z(i,j)表示第i行、第j列正方形高程区域的平均高程,n表示方格内包 含点的数量,Z(k)表示正方形高程区域内第k个点云的高程。
进一步的是,步骤S2中,通过公式
判断多边形顶点所位于的高程区域;
进一步的是,步骤S2中,以该多边形各顶点至该多边形中心的距离为权重, 对该多边形各顶点的平均高程取平均值,得到该多边形的计算高程。
进一步的是,所述地形的点云数据通过三维模型提取、激光扫描装置提取和 无人机航测三者之中一种或多种方式相结合得到。
同时,本发明提供一种基于点云的土方量计算系统,包括点云数据获取模 块,至少一个处理器,以及与所述至少一个处理器通信连接的存储器;所述存 储器存储有可被所述至少一个处理器执行的指令,所述指令被所述至少一个处 理器执行,以使所述至少一个处理器能够执行上述计算方法。通过该系统将冗 杂的点云数据划分为不同的高程区域,简化了土方计算所需的高程计算;将设 计平面划分为若干计算平面,进一步的简化了土方计算所需的计算量;同时保 留了原有设计平面的设计边界,不会产生边界误差。
进一步的是,所述点云数据获取模块为激光扫描器或设置有激光扫描器的 飞行器。
同时,本发明还提供一种可读存储介质,其上存储有计算机程序,所述程 序被处理器执行实现上述计算方法。
与现有技术相比,本发明的有益效果是:通过划分高程区域在保证了计算精 度的前提下,可以简化土方的计算量;并且本发明采用统一的设计边界进行土 方量计算,计算边界和设计边界一致,避免了边界误差。
附图说明:
图1为本发明示例性实施例中一种基于点云的土方量计算方法逻辑框图;
图2为本发明示例性实施例中设计平面的分割图;
图3为本发明示例性实施例中的设计平面划分示意图;
图4为本发明示例性实施例中高程区域示意图;
图5为本发明示例性实施例中基坑区域、点云区域的示意图;
图6为本发明示例性实施例提供的一种基于点云的土方量计算系统示意图。
具体实施方式
下面结合试验例及具体实施方式对本发明作进一步的详细描述。但不应将 此理解为本发明上述主题的范围仅限于以下的实施例,凡基于本发明内容所实 现的技术均属于本发明的范围。
实施例1
图1示出了本发明示例性实施例中一种基于点云的土方量计算方法逻辑框图。 该实施例的计算方法主要包括:
S1.获取地形点云数据和设计平面,将点云数据划分为若干高程区域;
分别计算若干高程区域内点云高度坐标的平均值得到平均高程;
具体的,所述步骤S1可以包括:
S11.将点云区域划分为若干边长相同的正方形高程区域,正方形高程区 域边长d的计算公式为:
其中,Xmax表示点云中所有点X坐标的最大值,Xmin表示点云中所有点 X坐标的最小值,Ymax表示点云中所有点Y坐标的最大值,Ymin表示点云中所 有点Y坐标的最小值;
S12.计算每一个正方形高程区域内点云的平均高程,
其中,Z(i,j)表示第i行、第j列正方形高程区域的平均高程,n表示方格 内包含点的数量,Z(k)表示正方形高程区域内第k个点云的高程。通过上述方 法得到的高程区域如图2所示;为了保证数据计算的精准性,同时又可以节约 算力,本发明将数量庞大、复杂的点云划分为若干正方形的高程区域,并通过 高程区域内的点云高程平均值表征该高程区域的高程,简化了土方量的计算, 节省了计算时所需的算力。
S2.将设计平面划分为若干多边形,获取所述若干多边形的顶点所在的高程 区域,将高程区域的平均高程作为所述若干多边形的顶点的平均高程;并通过 多边形顶点的平均高程得到该多边形的计算高程;
进一步的是,步骤S2可以包括如下步骤:
S21.将设计平面中点与设计平面各顶点依次相连,将设计平面划分成若 干个三角形;
S22.将三角形与设计平面中点相连的两边等分,并连接两边对应的等分 点;并将连接三角形两边等分点的线段和三角形不与设计平面顶点的线段继续 等分;将新建立的等分点与设计平面顶点连接,得到若干计算平面。
S23.将计算平面进一步划分,划分为若干多边形,并获取若干多边形顶点 所在的高程区域,将改高程区域的平均高程作为所述若干多边形的顶点的平均 高程;同时以该多边形各顶点至该多边形中心的距离为权重,对该多边形各顶 点的平均高程取平均值,得到该多边形的计算高程。
作为示例性的实施方式,如图2所示,分别将设计平片中点与设计平面各 顶点连接,将设计平面划分为若干三角形。
S23.将各个三角形再次进行划分,如图3所示,O为设计平面中点,A、B 为设计平面的两个顶点。将OA和OB分别等分成u段,再将各等分点互相连接, 形成CD、EF等若干线段,再将CD、EF等线段再依次等分为w段,并相互连接 各等分点,这样,就将初始的一个大三角形,划分为1个小三角形和若干个小 四边形。u和w采用如下公式计算:
其中SOA表示OA的长度,SOB表示OB的长度,SCD表示CD的长度,SEF表示EF 的长度。
通过该步骤,将复杂设计平面简化、拆分;既可以保证计算的精度又可以 合理的减少所需计算量和算力。
S24.具体的,可以通过如下公式判断三角形或者四边形的顶点位于哪个正 方形方格:
在判断了三角形或者四边形的顶点位于哪个正方形方格后通过如下公式 判断三角形计算区域的计算高程:
S123=S1+S2+S3
S1、S2、S3分表表示三角形三个顶点到三角形中心点的距离,S123为三角形 三个顶点到三角形中心点距离和。
同样,基于上述公式的方法也可以计算四边形计算区域的计算高程。
S3.通过所述若干多边形的面积和计算高程得到土方量,并通过土方量的 正负判断是挖方量或填方量。
作为优选的,所述地形的点云数据可以通过三维模型提取、激光扫描装置 提取和无人机航测三者之中一种或多种方式相结合得到。
实施例2
某基坑工程如图5所示,该基坑所在区域的点云数据包含1000万个三维坐 标数据,在英特尔(R)核心(TM)I5 CPU,和8GB内存的PC环境下,分别采取本 发明方法和三角网格发进行挖方量计算,二者运算时间如表1所示,
表1
方法 | 电脑运算时长 |
本发明所述方法 | 1分45秒 |
三角网法 | 2小时11分 |
由此可见,本发明提出的方法很大程度上节约了电脑的算力,进而使得该 方法的应用较为灵活,并可以在低性能设备上使用。
实施例3
图6示出了根据本发明示例性实施例的一种基于点云的土方量计算系统, 即电子设备310(例如具备程序执行功能的计算机服务器),其包括至少一个处 理器311,电源314,点云数据获取模块,以及与所述至少一个处理器311通信 连接的存储器312和输入输出接口313;所述存储器312存储有可被所述至少一 个处理器311执行的指令,所述指令被所述至少一个处理器311执行,以使所 述至少一个处理器311能够执行前述任一实施例所公开的方法;所述输入输出 接口313可以包括显示器、键盘、鼠标、以及USB接口,用于输入输出数据; 电源314用于为电子设备310提供电能。
其中点云数据获取模块可以为激光扫描器或设置有激光扫描器的无人飞行 器、飞机等;无人飞行器可以为:无人固定翼飞机、无人垂直起降飞机、无人 飞艇、无人直升机、无人多旋翼飞行器、无人伞翼机等;飞机可以为直升飞机; 同时点云数据获取模块还可以通过已经建立好的三维模型如BIM模型中提取点 云数据。
本领域技术人员可以理解:实现上述方法实施例的全部或部分步骤可以通 过程序指令相关的硬件来完成,前述的程序可以存储于计算机可读取存储介质 中,该程序在执行时,执行包括上述方法实施例的步骤;而前述的存储介质包 括:移动存储设备、只读存储器(Read Only Memory,ROM)、磁碟或者光盘等 各种可以存储程序代码的介质。
当本发明上述集成的单元以软件功能单元的形式实现并作为独立的产品销 售或使用时,也可以存储在一个计算机可读取存储介质中。基于这样的理解, 本发明实施例的技术方案本质上或者说对现有技术做出贡献的部分可以以软件 产品的形式体现出来,该计算机软件产品存储在一个存储介质中,包括若干指 令用以使得一台计算机设备(可以是个人计算机、服务器、或者网络设备等) 执行本发明各个实施例所述方法的全部或部分。而前述的存储介质包括:移动 存储设备、ROM、磁碟或者光盘等各种可以存储程序代码的介质。
以上所述仅为本发明的较佳实施例而已,并不用以限制本发明,凡在本发明 的精神和原则之内所作的任何修改、等同替换和改进等,均应包含在本发明的 保护范围之内。
Claims (8)
1.一种基于点云的土方量计算方法,其特征在于,包括如下步骤:
S1.获取地形点云数据和设计平面,将点云数据划分为若干高程区域;分别计算若干高程区域内点云高度坐标的平均值得到该高程区域的平均高程;
S2.将设计平面划分为若干多边形,获取所述若干多边形的顶点所在的高程区域,将高程区域的平均高程作为所述若干多边形的顶点的平均高程;并通过多边形顶点的平均高程得到该多边形的计算高程;
S3.通过所述若干多边形的面积和计算高程得到土方量,并通过土方量的正负判断是挖方量或填方量。
4.如权利要求3所述的一种基于点云的土方量计算方法,其特征在于,步骤S2中,以该多边形各顶点至该多边形中心的距离为权重,对该多边形各顶点的平均高程取平均值,得到该多边形的计算高程。
5.如权利要求1至4任意一项所述的一种基于点云的土方量计算方法,其特征在于,所述地形的点云数据通过三维模型提取、激光扫描装置提取和无人机航测三者之中一种或多种方式相结合得到。
6.一种基于点云的土方量计算系统,其特征在于,包括点云数据获取模块,至少一个处理器,以及与所述至少一个处理器通信连接的存储器;所述存储器存储有可被所述至少一个处理器执行的指令,所述指令被所述至少一个处理器执行,以使所述至少一个处理器能够执行权利要求1至5中任一项所述的方法。
7.如权利要求6所述的一种计算系统,其特征在于,所述点云数据获取模块为激光扫描器或设置有激光扫描器的飞行器。
8.一种可读存储介质,其上存储有计算机程序,其特征在于,所述程序被处理器执行实现权利要求1至5中任一项所述的计算方法。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110920858.XA CN113963050A (zh) | 2021-08-11 | 2021-08-11 | 一种基于点云的土方量计算方法及系统 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110920858.XA CN113963050A (zh) | 2021-08-11 | 2021-08-11 | 一种基于点云的土方量计算方法及系统 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN113963050A true CN113963050A (zh) | 2022-01-21 |
Family
ID=79460501
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110920858.XA Pending CN113963050A (zh) | 2021-08-11 | 2021-08-11 | 一种基于点云的土方量计算方法及系统 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113963050A (zh) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115469328A (zh) * | 2022-08-29 | 2022-12-13 | 中国电建集团福建省电力勘测设计院有限公司 | 一种基于激光点云的数字化基坑验收方法 |
CN115482269A (zh) * | 2022-09-22 | 2022-12-16 | 佳都科技集团股份有限公司 | 一种土方量计算方法、装置、终端设备以及存储介质 |
CN115797288A (zh) * | 2022-11-30 | 2023-03-14 | 安徽开源路桥有限责任公司 | 基于地面点云数据计算填挖方量的方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN205262429U (zh) * | 2015-12-30 | 2016-05-25 | 中铁三局集团华东建设有限公司 | 一种施工现场土方量测定系统 |
CN105783810A (zh) * | 2016-04-15 | 2016-07-20 | 昆山数字城市信息技术有限公司 | 基于无人机摄影技术的工程土方量测量方法 |
US20200032490A1 (en) * | 2018-07-26 | 2020-01-30 | Built Robotics Inc. | Filling earth at a location within a dig site using an excavation vehicle |
CN112734929A (zh) * | 2020-12-28 | 2021-04-30 | 三峡大学 | 基于网格细分算法的土石坝复杂土料场开挖方量计算方法 |
-
2021
- 2021-08-11 CN CN202110920858.XA patent/CN113963050A/zh active Pending
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN205262429U (zh) * | 2015-12-30 | 2016-05-25 | 中铁三局集团华东建设有限公司 | 一种施工现场土方量测定系统 |
CN105783810A (zh) * | 2016-04-15 | 2016-07-20 | 昆山数字城市信息技术有限公司 | 基于无人机摄影技术的工程土方量测量方法 |
US20200032490A1 (en) * | 2018-07-26 | 2020-01-30 | Built Robotics Inc. | Filling earth at a location within a dig site using an excavation vehicle |
CN112734929A (zh) * | 2020-12-28 | 2021-04-30 | 三峡大学 | 基于网格细分算法的土石坝复杂土料场开挖方量计算方法 |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115469328A (zh) * | 2022-08-29 | 2022-12-13 | 中国电建集团福建省电力勘测设计院有限公司 | 一种基于激光点云的数字化基坑验收方法 |
CN115482269A (zh) * | 2022-09-22 | 2022-12-16 | 佳都科技集团股份有限公司 | 一种土方量计算方法、装置、终端设备以及存储介质 |
CN115797288A (zh) * | 2022-11-30 | 2023-03-14 | 安徽开源路桥有限责任公司 | 基于地面点云数据计算填挖方量的方法 |
CN115797288B (zh) * | 2022-11-30 | 2023-09-29 | 安徽开源路桥有限责任公司 | 基于地面点云数据计算填挖方量的方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN113963050A (zh) | 一种基于点云的土方量计算方法及系统 | |
CN108010092A (zh) | 一种基于低空摄影测量的城市高密度区太阳能利用潜力评估方法 | |
CN103020342B (zh) | 一种从地面LiDAR数据中提取建筑物轮廓和角点的方法 | |
CN109143257A (zh) | 无人机机载雷达矿山开采土地变化监测系统及方法 | |
CN103217688B (zh) | 一种基于不规则三角网机载激光雷达点云平差计算方法 | |
CN107767453B (zh) | 一种基于规则约束的建筑物lidar点云重构优化方法 | |
CN104835202A (zh) | 一种三维虚拟场景快速构建方法 | |
CN112100715A (zh) | 一种基于三维倾斜摄影技术的土方优化方法及系统 | |
CN102708587A (zh) | 一种快速获取三维建筑信息的方法及系统 | |
CN105761310B (zh) | 一种天空可视域数字地图的模拟分析及图像显示方法 | |
CN109508508A (zh) | 一种露天矿山治理勘查设计方法 | |
CN103871102A (zh) | 一种基于高程点和道路轮廓面的道路三维精细建模方法 | |
CN103256914B (zh) | 一种基于dem计算淤地坝淹没面积的方法及系统 | |
CN113566781B (zh) | 一种机场道面板角高程测量方法、装置、介质及电子设备 | |
KR100657870B1 (ko) | 항공 레이저 측량 데이터를 이용한 지반고 추출 방법 | |
CN114580310A (zh) | 一种基于palm实现wrf模拟风场降尺度处理的方法 | |
CN104809692A (zh) | 一种等高线生成dem的内插方法 | |
CN104931970B (zh) | 一种面向机载气象雷达仿真的三维云场生成方法 | |
CN113743027A (zh) | 一种基于cfd技术绘制风资源图谱的方法和装置 | |
CN106557600A (zh) | 一种输电塔长短腿及基础的仿真方法 | |
CN109635317A (zh) | 一种高海拔地区cfd仿真地形拓扑方法 | |
CN117523125A (zh) | 一种建筑物三维体块模型快速构建方法 | |
CN108663014B (zh) | 一种垄沟耕作坡耕地田块相对高程数据计算方法 | |
CN111898819B (zh) | 空间网格划分方法及装置 | |
CN114565735A (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 | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20220121 |