CN112906124A - 一种土方量计算方法、装置、设备和介质 - Google Patents

一种土方量计算方法、装置、设备和介质 Download PDF

Info

Publication number
CN112906124A
CN112906124A CN202110360582.4A CN202110360582A CN112906124A CN 112906124 A CN112906124 A CN 112906124A CN 202110360582 A CN202110360582 A CN 202110360582A CN 112906124 A CN112906124 A CN 112906124A
Authority
CN
China
Prior art keywords
vertex
dimensional model
live
action
points
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
Application number
CN202110360582.4A
Other languages
English (en)
Other versions
CN112906124B (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.)
South GNSS Navigation Co Ltd
Original Assignee
South GNSS Navigation Co Ltd
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 South GNSS Navigation Co Ltd filed Critical South GNSS Navigation Co Ltd
Priority to CN202110360582.4A priority Critical patent/CN112906124B/zh
Publication of CN112906124A publication Critical patent/CN112906124A/zh
Application granted granted Critical
Publication of CN112906124B publication Critical patent/CN112906124B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/10—Geometric CAD
    • G06F30/13—Architectural design, e.g. computer-aided architectural design [CAAD] related to design of buildings, bridges, landscapes, production plants or roads
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/20—Design optimisation, verification or simulation
    • 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/10—Constructive solid geometry [CSG] using solid primitives, e.g. cylinders, cubes
    • 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
    • Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
    • Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
    • Y02A90/10—Information and communication technologies [ICT] supporting adaptation to climate change, e.g. for weather forecasting or climate simulation

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Geometry (AREA)
  • Theoretical Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Software Systems (AREA)
  • Computer Graphics (AREA)
  • General Engineering & Computer Science (AREA)
  • Architecture (AREA)
  • Civil Engineering (AREA)
  • Structural Engineering (AREA)
  • Computational Mathematics (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Pure & Applied Mathematics (AREA)
  • Image Processing (AREA)
  • Processing Or Creating Images (AREA)

Abstract

本申请公开了一种土方量计算方法、装置、设备和介质，获取待测量区域的实景三维模型，并对实景三维模型进行三角网化，获取实景三维模型的顶点；对实景三维模型的顶点中的非地面点进行过滤，得到过滤后的实景三维模型；当过滤后的实景三维模型存在空洞区域时，对空洞区域进行插值，并通过最小二乘法对插值后的空洞区域进行曲面拟合，得到拟合曲面；提取过滤后的实景三维模型和拟合曲面的高程点，并基于高程点计算待测量区域的土方量。本申请改善了现有的土方量计算方法需要对原始地形实时测量大量的野外测量点的真实平面坐标数据和高程数据，存在测量过程复杂，计算效率低的技术问题。

Description

一种土方量计算方法、装置、设备和介质
技术领域
本申请涉及土方量计算技术领域，尤其涉及一种土方量计算方法、装置、设备和介质。
背景技术
土方量的计算是工程施工的一个重要步骤，工程设计阶段必须对土石方量进行预算，它直接关系到工程的费用及方案选优。在进行土石方量计算时，必须对原始地形进行测绘，根据土方量计算精度和地形条件确定测图比例尺，并且必须标明原始地形高程、原设计高程和变更设计高程，该方法需要进行实时测量，以获取大量的野外测量点的真实平面坐标数据和高程数据来计算土方量，测量过程复杂，计算效率低。
发明内容
本申请提供了一种土方量计算方法、装置、设备和介质，用于改善现有的土方量计算方法需要对原始地形实时测量大量的野外测量点的真实平面坐标数据和高程数据，存在测量过程复杂，计算效率低的技术问题。
有鉴于此，本申请第一方面提供了一种土方量计算方法，包括：
获取待测量区域的实景三维模型，并对所述实景三维模型进行三角网化，获取三角形的顶点，得到所述实景三维模型的顶点；
对所述实景三维模型的顶点中的非地面点进行过滤，得到过滤后的所述实景三维模型；
当过滤后的所述实景三维模型存在空洞区域时，对所述空洞区域进行插值，并通过最小二乘法对插值后的所述空洞区域进行曲面拟合，得到拟合曲面；
提取过滤后的所述实景三维模型和所述拟合曲面的高程点，并基于所述高程点计算所述待测量区域的土方量。
可选的，所述对所述实景三维模型的顶点中的非地面点进行过滤，得到过滤后的所述实景三维模型，包括：
从所述实景三维模型的顶点中提取扫描线；
获取各所述扫描线上的最低顶点和最高顶点，所述最低顶点为地面点，所述最高顶点为非地面点；
基于各所述扫描线上的所述最低顶点对各所述扫描线上的顶点进行粗分类，得到非地面点集；
以各所述扫描线上的所述最高顶点为中心点设置预置大小的等高线，将所述等高线内的顶点加入到所述非地面点集，所述等高线的形状为矩形或不规则多边形；
基于所述非地面点集对所述实景三维模型的顶点中的非地面点进行过滤，得到过滤后的所述实景三维模型。
可选的，所述基于各所述扫描线上的所述最低顶点对各所述扫描线上的顶点进行粗分类，得到非地面点集，包括：
将各所述扫描线上的所述最低顶点作为目标顶点；
计算各所述扫描线上的所述目标顶点和所述目标顶点的相邻顶点的坡度值，并基于所述坡度值判断所述目标顶点的相邻顶点是否为非地面点；
将所述目标顶点的相邻顶点作为新的所述目标顶点，返回所述计算各所述扫描线上的所述目标顶点和所述目标顶点的相邻顶点的坡度值，并基于所述坡度值判断所述目标顶点的相邻顶点是否为非地面点的步骤，直至遍历完各所述扫描线上的所述顶点，得到非地面点集。
可选的，所述将各所述扫描线上的所述最低顶点作为目标顶点，之前还包括：
检测各所述扫描线上的顶点中是否存在离群点，若是，则剔除所述离群点。
可选的，所述当过滤后的所述实景三维模型存在空洞区域时，对所述空洞区域进行插值，包括：
当过滤后的所述实景三维模型存在空洞区域时，将所述空洞区域的边界点与所述空洞区域周围预置大小的区域内的顶点构建三角网；
通过三角网插值法在所述空洞区域***点，当所述空洞区域的边界特征以直线为主时，***的该点的高程值为所述空洞区域的边界点的高程均值，当所述空洞区域的边界特征以曲线为主时，***的该点的高程值通过距离该点最近的三角网插值得到。
可选的，所述基于所述高程点计算所述待测量区域的土方量，包括：
基于所述高程点生成三角网，得到若干三角锥；或，将所述待测量区域进行方格划分，通过插值获取各方格顶点的高程值，基于各方格顶点和所述高程点生成若干三角锥；
计算各所述三角锥的体积，并统计所有所述三角锥的体积，得到所述待测量区域的土方量。
本申请第二方面提供了一种土方量计算装置，包括：
获取单元，用于获取待测量区域的实景三维模型，并对所述实景三维模型进行三角网化，获取三角形的顶点，得到所述实景三维模型的顶点；
过滤单元，用于对所述实景三维模型的顶点中的非地面点进行过滤，得到过滤后的所述实景三维模型；
插值和拟合单元，用于当过滤后的所述实景三维模型存在空洞区域时，对所述空洞区域进行插值，并通过最小二乘法对插值后的所述空洞区域进行曲面拟合，得到拟合曲面；
计算单元，提取过滤后的所述实景三维模型和所述拟合曲面的高程点，并基于所述高程点计算所述待测量区域的土方量。
可选的，所述过滤单元具体包括：
提取子单元，用于从所述实景三维模型的顶点中提取扫描线；
获取子单元，用于获取各所述扫描线上的最低顶点和最高顶点，所述最低顶点为地面点，所述最高顶点为非地面点；
分类子单元，用于基于各所述扫描线上的所述最低顶点对各所述扫描线上的顶点进行粗分类，得到非地面点集；
设置子单元，用于以各所述扫描线上的所述最高顶点为中心点设置预置大小的等高线，将所述等高线内的顶点加入到所述非地面点集，所述等高线的形状为矩形或不规则多边形；
过滤子单元，用于基于所述非地面点集对所述实景三维模型的顶点中的非地面点进行过滤，得到过滤后的所述实景三维模型。
本申请第三方面提供了一种土方量计算设备，所述设备包括处理器以及存储器；
所述存储器用于存储程序代码，并将所述程序代码传输给所述处理器；
所述处理器用于根据所述程序代码中的指令执行第一方面任一种所述的土方量计算方法。
本申请第四方面提供了一种计算机可读存储介质，所述计算机可读存储介质用于存储程序代码，所述程序代码用于执行第一方面任一种所述的土方量计算方法。
从以上技术方案可以看出，本申请具有以下优点：
本申请提供了一种土方量计算方法，包括：获取待测量区域的实景三维模型，并对实景三维模型进行三角网化，获取三角形的顶点，得到实景三维模型的顶点；对实景三维模型的顶点中的非地面点进行过滤，得到过滤后的实景三维模型；当过滤后的实景三维模型存在空洞区域时，对空洞区域进行插值，并通过最小二乘法对插值后的空洞区域进行曲面拟合，得到拟合曲面；提取过滤后的实景三维模型和拟合曲面的高程点，并基于高程点计算待测量区域的土方量。
本申请中，获取待测量区域的实景三维模型，对实景三维模型进行一系列处理，以过滤掉实景三维模型中的地表建筑和数木等非地面点，然后获取高程数据来计算土方量，不需要进行大量的复杂的测量工作，计算效率高；进一步，当过滤后的实景三维模型存在空洞区域时，对空洞区域进行插值，并通过最小二乘法对插值后的空洞区域进行曲面拟合，形成符合真实地表形状的数字高程模型，保证了土方量计算结果的准确性，从而在保证土方量计算结果的准确性的基础上，改善了现有的土方量计算方法需要对原始地形实时测量大量的野外测量点的真实平面坐标数据和高程数据，存在测量过程复杂，计算效率低的技术问题。
附图说明
为了更清楚地说明本申请实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本申请的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动性的前提下，还可以根据这些附图获得其它的附图。
图1为本申请实施例提供的一种土方量计算方法的一个流程示意图；
图2为本申请实施例提供的一种土方量计算装置的一个结构示意图。
具体实施方式
本申请提供了一种土方量计算方法、装置、设备和介质，用于改善现有的土方量计算方法需要对原始地形实时测量大量的野外测量点的真实平面坐标数据和高程数据，存在测量过程复杂，计算效率低的技术问题。
为了使本技术领域的人员更好地理解本申请方案，下面将结合本申请实施例中的附图，对本申请实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅是本申请一部分实施例，而不是全部的实施例。基于本申请中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本申请保护的范围。
为了便于理解，请参阅图1，本申请提供的一种土方量计算方法的一个实施例，包括：
步骤101、获取待测量区域的实景三维模型，并对实景三维模型进行三角网化，获取三角形的顶点，得到实景三维模型的顶点。
可以通过摄像头等采集待测量区域的图像数据，然后基于图像数据构建该待测量区域的实景三维模型。对该实景三维模型进行三角网化，提取三角形顶点，得到该实景三维模型的顶点。所有的顶点构成的三角形形成了该实景三维模型的最小的纹理贴图单位，可以根据实际需要进行三角形顶点的平滑插值或抽稀，以平衡计算精度和计算时间。其中，在进行三角网插值时，设定一个相对水平面的斜率阈值a，根据该阈值a划分若干个等级，在插值计算过程中，不同等级对应三角形顶点V1、V2、V3有不同的权重，以此来保证***的点接近真实模型。
本申请实施例通过采集待测量区域的图像数据构建实景三维模型，来获取高程数据，不需要进行复杂的测量工作。
步骤102、对实景三维模型的顶点中的非地面点进行过滤，得到过滤后的实景三维模型。
在计算土方量之前，需过滤掉顶点中的非地面点。具体过程如下：
S1、从实景三维模型的顶点中提取扫描线。
可以通过扫描线方法从实景三维模型的顶点中提取扫描线。
S2、获取各扫描线上的最低顶点和最高顶点，最低顶点为地面点，最高顶点为非地面点。
S3、基于各扫描线上的最低顶点对各扫描线上的顶点进行粗分类，得到非地面点集。
基于各扫描线上的最低顶点对各扫描线上的顶点进行粗分类，得到非地面点集的具体过程为：
S30、将各扫描线上的最低顶点作为目标顶点；
在获取最低顶点之前，还可以根据离群点检测算法检测各扫描线上的顶点中是否存在离群点，若是，则剔除离群点，若否，则获取各扫描线上的最低顶点。
S31、计算各扫描线上的目标顶点和目标顶点的相邻顶点的坡度值，并基于坡度值判断目标顶点的相邻顶点是否为非地面点；
计算各扫描线上的目标顶点和目标顶点的相邻顶点的坡度值，通过将坡度值与坡度阈值进行比较判断目标顶点的相邻顶点是地面点还是非地面点，其中，坡度值的计算属于现有技术，在此不再进行赘述。
S32、将目标顶点的相邻顶点作为新的目标顶点，返回步骤S31，直至遍历完各扫描线上的顶点，得到非地面点集。
S4、以各扫描线上的最高顶点为中心点设置预置大小的等高线，将等高线内的顶点加入到非地面点集，等高线的形状为矩形或不规则多边形；
在进行上述的粗分类后，还可以进一步再进行精细的过滤，将各扫描线上的最高顶点作为标记点，该点通常为建筑物或树冠等的中心点位置，定义一个矩形或不规则多边形进行向外扩展，按等高线扩展，定义线间阈值L，以及圈内点集密度u，对非地面点进行再次过滤，过滤时调整线间阈值L及密度u，通过把往外扩张的阈值把满足条件的点(即在等高线内的顶点)加入非地面点集中。
S5、基于非地面点集对实景三维模型的顶点中的非地面点进行过滤，得到过滤后的实景三维模型。
根据获取的非地面点集，对于实景三维模型中为非地面点的顶点进行过滤，从而去除地表建筑物和数木等非地面点，保留地面点。
步骤103、当过滤后的实景三维模型存在空洞区域时，对空洞区域进行插值，并通过最小二乘法对插值后的空洞区域进行曲面拟合，得到拟合曲面。
去除非地面点后的实景三维模型会存在空洞区域，需要对空洞区域进行插值。当过滤后的实景三维模型存在空洞区域时，将空洞区域的边界点与空洞区域周围预置大小的区域内的顶点构建三角网；通过三角网插值法在空洞区域***点，当空洞区域的边界特征以直线为主时，***的该点的高程值为空洞区域的边界点的高程均值，当空洞区域的边界特征以曲线为主时，***的该点的高程值通过距离该点最近的三角网插值得到。空洞区域的插值方法有两种。根据空洞区域的边界特征来选择空洞区域内部高程的变化类型，若边界特征以直线为主，则认为空洞区域内部的高程值为均值；若边界特征以曲线为主，则***空洞区域的点的高程值通过该点最近的三角网插值得到。最后采用最小二乘法进行曲面拟合，形成与去除非地面点后的真实地形接近的曲面。
步骤104、提取过滤后的实景三维模型和拟合曲面的高程点，并基于高程点计算待测量区域的土方量。
提取过滤后的实景三维模型和拟合曲面的高程点时，根据三维模型的存储模式，根据指定网格宽度来提取高程点。为了提高提取效率，可以采用多线程的方式分区域提取高程点，最后合并高程点。
基于高程点计算待测量区域的土方量有两种方法，在一种实施例中，基于高程点生成三角网，得到若干三角锥；计算各三角锥的体积，并统计所有三角锥的体积，得到待测量区域的土方量。根据上述获取的高程点生成三角网，使整个待测量区域的地形形成由三角锥组合的集合。根据用户给定的设计高程确定零平面，以“+”表示填方进行体积计算，以“-”表示挖方进行体积计算，最后统计所有的三角锥的挖填方量，得到待测量区域的土方量。
在另一种实施例中，将待测量区域进行方格划分，通过插值获取各方格顶点的高程值，基于各方格顶点和高程点生成若干三角锥；计算各三角锥的体积，并统计所有三角锥的体积，得到待测量区域的土方量。将待测量区域按20m或40m等预置大小进行划分，形成方格，方格顶点的高程值通过插值产生，基于方格顶点与前述步骤提取的高程点生成三角锥，同样的，以“+”表示填方进行体积计算，以“-”表示挖方进行体积计算，最后统计所有的三角锥的挖填方量，得到待测量区域的土方量。
本申请实施例中，获取待测量区域的实景三维模型，对实景三维模型进行一系列处理，以过滤掉实景三维模型中的地表建筑和数木等非地面点，然后获取高程数据来计算土方量，不需要进行大量的复杂的测量工作，计算效率高；进一步，当过滤后的实景三维模型存在空洞区域时，对空洞区域进行插值，并通过最小二乘法对插值后的空洞区域进行曲面拟合，形成符合真实地表形状的数字高程模型，保证了土方量计算结果的准确性，从而在保证土方量计算结果的准确性的基础上，改善了现有的土方量计算方法需要对原始地形实时测量大量的野外测量点的真实平面坐标数据和高程数据，存在测量过程复杂，计算效率低的技术问题。
以上为本申请提供的一种土方量计算方法的一个实施例，以下为本申请提供的一种土方量计算装置的一个实施例。
请参考图2，本申请实施例提供的一种土方量计算装置，包括：
获取单元，用于获取待测量区域的实景三维模型，并对实景三维模型进行三角网化，获取三角形的顶点，得到实景三维模型的顶点；
过滤单元，用于对实景三维模型的顶点中的非地面点进行过滤，得到过滤后的实景三维模型；
插值和拟合单元，用于当过滤后的实景三维模型存在空洞区域时，对空洞区域进行插值，并通过最小二乘法对插值后的空洞区域进行曲面拟合，得到拟合曲面；
计算单元，提取过滤后的实景三维模型和拟合曲面的高程点，并基于高程点计算待测量区域的土方量。
作为进一步地改进，过滤单元具体包括：
提取子单元，用于从实景三维模型的顶点中提取扫描线；
获取子单元，用于获取各扫描线上的最低顶点和最高顶点，最低顶点为地面点，最高顶点为非地面点；
分类子单元，用于基于各扫描线上的最低顶点对各扫描线上的顶点进行粗分类，得到非地面点集；
设置子单元，用于以各扫描线上的最高顶点为中心点设置预置大小的等高线，将等高线内的顶点加入到非地面点集，等高线的形状为矩形或不规则多边形；
过滤子单元，用于基于非地面点集对实景三维模型的顶点中的非地面点进行过滤，得到过滤后的实景三维模型。
作为进一步地改进，分类子单元，具体用于：
将各扫描线上的最低顶点作为目标顶点；
计算各扫描线上的目标顶点和目标顶点的相邻顶点的坡度值，并基于坡度值判断目标顶点的相邻顶点是否为非地面点；
将目标顶点的相邻顶点作为新的目标顶点，返回计算各扫描线上的目标顶点和目标顶点的相邻顶点的坡度值，并基于坡度值判断目标顶点的相邻顶点是否为非地面点的步骤，直至遍历完各扫描线上的顶点，得到非地面点集。
作为进一步地改进，过滤单元还包括：
检测子单元，用于检测各扫描线上的顶点中是否存在离群点，若是，则剔除离群点。
本申请实施例中，获取待测量区域的实景三维模型，对实景三维模型进行一系列处理，以过滤掉实景三维模型中的地表建筑和数木等非地面点，然后获取高程数据来计算土方量，不需要进行大量的复杂的测量工作，计算效率高；进一步，当过滤后的实景三维模型存在空洞区域时，对空洞区域进行插值，并通过最小二乘法对插值后的空洞区域进行曲面拟合，形成符合真实地表形状的数字高程模型，保证了土方量计算结果的准确性，从而在保证土方量计算结果的准确性的基础上，改善了现有的土方量计算方法需要对原始地形实时测量大量的野外测量点的真实平面坐标数据和高程数据，存在测量过程复杂，计算效率低的技术问题。
本申请实施例还提供了一种土方量计算设备，设备包括处理器以及存储器；
存储器用于存储程序代码，并将程序代码传输给处理器；
处理器用于根据程序代码中的指令执行前述方法实施例中的土方量计算方法。
本申请实施例还提供了一种计算机可读存储介质，计算机可读存储介质用于存储程序代码，程序代码用于执行前述方法实施例中的土方量计算方法。
所属领域的技术人员可以清楚地了解到，为描述的方便和简洁，上述描述的装置和单元的具体工作过程，可以参考前述方法实施例中的对应过程，在此不再赘述。
在本申请所提供的几个实施例中，应该理解到，所揭露的装置和方法，可以通过其它的方式实现。例如，以上所描述的装置实施例仅仅是示意性的，例如，所述单元的划分，仅仅为一种逻辑功能划分，实际实现时可以有另外的划分方式，例如多个单元或组件可以结合或者可以集成到另一个***，或一些特征可以忽略，或不执行。另一点，所显示或讨论的相互之间的耦合或直接耦合或通信连接可以是通过一些接口，装置或单元的间接耦合或通信连接，可以是电性，机械或其它的形式。
所述作为分离部件说明的单元可以是或者也可以不是物理上分开的，作为单元显示的部件可以是或者也可以不是物理单元，即可以位于一个地方，或者也可以分布到多个网络单元上。可以根据实际的需要选择其中的部分或者全部单元来实现本实施例方案的目的。
另外，在本申请各个实施例中的各功能单元可以集成在一个处理单元中，也可以是各个单元单独物理存在，也可以两个或两个以上单元集成在一个单元中。上述集成的单元既可以采用硬件的形式实现，也可以采用软件功能单元的形式实现。
所述集成的单元如果以软件功能单元的形式实现并作为独立的产品销售或使用时，可以存储在一个计算机可读取存储介质中。基于这样的理解，本申请的技术方案本质上或者说对现有技术做出贡献的部分或者该技术方案的全部或部分可以以软件产品的形式体现出来，该计算机软件产品存储在一个存储介质中，包括若干指令用以通过一台计算机设备(可以是个人计算机，服务器，或者网络设备等)执行本申请各个实施例所述方法的全部或部分步骤。而前述的存储介质包括：U盘、移动硬盘、只读存储器(英文全称：Read-OnlyMemory，英文缩写：ROM)、随机存取存储器(英文全称：RandomAccess Memory，英文缩写：RAM)、磁碟或者光盘等各种可以存储程序代码的介质。
以上所述，以上实施例仅用以说明本申请的技术方案，而非对其限制；尽管参照前述实施例对本申请进行了详细的说明，本领域的普通技术人员应当理解：其依然可以对前述各实施例所记载的技术方案进行修改，或者对其中部分技术特征进行等同替换；而这些修改或者替换，并不使相应技术方案的本质脱离本申请各实施例技术方案的精神和范围。

Claims (10)

1.一种土方量计算方法，其特征在于，包括：
获取待测量区域的实景三维模型，并对所述实景三维模型进行三角网化，获取三角形的顶点，得到所述实景三维模型的顶点；
对所述实景三维模型的顶点中的非地面点进行过滤，得到过滤后的所述实景三维模型；
当过滤后的所述实景三维模型存在空洞区域时，对所述空洞区域进行插值，并通过最小二乘法对插值后的所述空洞区域进行曲面拟合，得到拟合曲面；
提取过滤后的所述实景三维模型和所述拟合曲面的高程点，并基于所述高程点计算所述待测量区域的土方量。
2.根据权利要求1所述的土方量计算方法，其特征在于，所述对所述实景三维模型的顶点中的非地面点进行过滤，得到过滤后的所述实景三维模型，包括：
从所述实景三维模型的顶点中提取扫描线；
获取各所述扫描线上的最低顶点和最高顶点，所述最低顶点为地面点，所述最高顶点为非地面点；
基于各所述扫描线上的所述最低顶点对各所述扫描线上的顶点进行粗分类，得到非地面点集；
以各所述扫描线上的所述最高顶点为中心点设置预置大小的等高线，将所述等高线内的顶点加入到所述非地面点集，所述等高线的形状为矩形或不规则多边形；
基于所述非地面点集对所述实景三维模型的顶点中的非地面点进行过滤，得到过滤后的所述实景三维模型。
3.根据根据权利要求2所述的土方量计算方法，其特征在于，所述基于各所述扫描线上的所述最低顶点对各所述扫描线上的顶点进行粗分类，得到非地面点集，包括：
将各所述扫描线上的所述最低顶点作为目标顶点；
计算各所述扫描线上的所述目标顶点和所述目标顶点的相邻顶点的坡度值，并基于所述坡度值判断所述目标顶点的相邻顶点是否为非地面点；
将所述目标顶点的相邻顶点作为新的所述目标顶点，返回所述计算各所述扫描线上的所述目标顶点和所述目标顶点的相邻顶点的坡度值，并基于所述坡度值判断所述目标顶点的相邻顶点是否为非地面点的步骤，直至遍历完各所述扫描线上的所述顶点，得到非地面点集。
4.根据权利要求3所述的土方量计算方法，其特征在于，所述将各所述扫描线上的所述最低顶点作为目标顶点，之前还包括：
检测各所述扫描线上的顶点中是否存在离群点，若是，则剔除所述离群点。
5.根据权利要求1所述的土方量计算方法，其特征在于，所述当过滤后的所述实景三维模型存在空洞区域时，对所述空洞区域进行插值，包括：
当过滤后的所述实景三维模型存在空洞区域时，将所述空洞区域的边界点与所述空洞区域周围预置大小的区域内的顶点构建三角网；
通过三角网插值法在所述空洞区域***点，当所述空洞区域的边界特征以直线为主时，***的该点的高程值为所述空洞区域的边界点的高程均值，当所述空洞区域的边界特征以曲线为主时，***的该点的高程值通过距离该点最近的三角网插值得到。
6.根据权利要求1所述的土方量计算方法，其特征在于，所述基于所述高程点计算所述待测量区域的土方量，包括：
基于所述高程点生成三角网，得到若干三角锥；或，将所述待测量区域进行方格划分，通过插值获取各方格顶点的高程值，基于各方格顶点和所述高程点生成若干三角锥；
计算各所述三角锥的体积，并统计所有所述三角锥的体积，得到所述待测量区域的土方量。
7.一种土方量计算装置，其特征在于，包括：
获取单元，用于获取待测量区域的实景三维模型，并对所述实景三维模型进行三角网化，获取三角形的顶点，得到所述实景三维模型的顶点；
过滤单元，用于对所述实景三维模型的顶点中的非地面点进行过滤，得到过滤后的所述实景三维模型；
插值和拟合单元，用于当过滤后的所述实景三维模型存在空洞区域时，对所述空洞区域进行插值，并通过最小二乘法对插值后的所述空洞区域进行曲面拟合，得到拟合曲面；
计算单元，提取过滤后的所述实景三维模型和所述拟合曲面的高程点，并基于所述高程点计算所述待测量区域的土方量。
8.根据权利要求6所述的土方量计算装置，其特征在于，所述过滤单元具体包括：
提取子单元，用于从所述实景三维模型的顶点中提取扫描线；
获取子单元，用于获取各所述扫描线上的最低顶点和最高顶点，所述最低顶点为地面点，所述最高顶点为非地面点；
分类子单元，用于基于各所述扫描线上的所述最低顶点对各所述扫描线上的顶点进行粗分类，得到非地面点集；
设置子单元，用于以各所述扫描线上的所述最高顶点为中心点设置预置大小的等高线，将所述等高线内的顶点加入到所述非地面点集，所述等高线的形状为矩形或不规则多边形；
过滤子单元，用于基于所述非地面点集对所述实景三维模型的顶点中的非地面点进行过滤，得到过滤后的所述实景三维模型。
9.一种土方量计算设备，其特征在于，所述设备包括处理器以及存储器；
所述存储器用于存储程序代码，并将所述程序代码传输给所述处理器；
所述处理器用于根据所述程序代码中的指令执行权利要求1-6任一项所述的土方量计算方法。
10.一种计算机可读存储介质，其特征在于，所述计算机可读存储介质用于存储程序代码，所述程序代码用于执行权利要求1-6任一项所述的土方量计算方法。
CN202110360582.4A 2021-04-02 2021-04-02 一种土方量计算方法、装置、设备和介质 Active CN112906124B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110360582.4A CN112906124B (zh) 2021-04-02 2021-04-02 一种土方量计算方法、装置、设备和介质

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110360582.4A CN112906124B (zh) 2021-04-02 2021-04-02 一种土方量计算方法、装置、设备和介质

Publications (2)

Publication Number Publication Date
CN112906124A true CN112906124A (zh) 2021-06-04
CN112906124B CN112906124B (zh) 2023-06-13

Family

ID=76109880

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110360582.4A Active CN112906124B (zh) 2021-04-02 2021-04-02 一种土方量计算方法、装置、设备和介质

Country Status (1)

Country Link
CN (1) CN112906124B (zh)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114527461A (zh) * 2021-12-22 2022-05-24 上海振懋机电科技有限公司 一种空间测量装置及其测量方法
CN115482269A (zh) * 2022-09-22 2022-12-16 佳都科技集团股份有限公司 一种土方量计算方法、装置、终端设备以及存储介质
CN118052947A (zh) * 2024-03-28 2024-05-17 四川省公路规划勘察设计研究院有限公司 基于大数据的三维地理模型建立方法及设备

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2002081008A (ja) * 2000-09-11 2002-03-22 Kokusai Kogyo Co Ltd デジタル地図データを用いた土工設計支援システム
JP2002203234A (ja) * 2000-12-28 2002-07-19 Inst Of Physical & Chemical Res 等高線データの作成方法と装置
CN109191581A (zh) * 2018-08-29 2019-01-11 李东帅 一种基于bim技术的道路工程土方量高精度计算方法
CN109658477A (zh) * 2017-10-12 2019-04-19 西南科技大学 一种基于lidar数据的dem生成算法
CN110285792A (zh) * 2019-07-02 2019-09-27 山东省交通规划设计院 一种无人机倾斜摄影精细格网土方计量方法
CN110657855A (zh) * 2019-09-23 2020-01-07 广州市中海达测绘仪器有限公司 土方量测量方法、装置、计算机设备和存储介质
CN111325837A (zh) * 2020-01-23 2020-06-23 江西理工大学 一种基于地面三维激光点云的边坡dem生成方法

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2002081008A (ja) * 2000-09-11 2002-03-22 Kokusai Kogyo Co Ltd デジタル地図データを用いた土工設計支援システム
JP2002203234A (ja) * 2000-12-28 2002-07-19 Inst Of Physical & Chemical Res 等高線データの作成方法と装置
CN109658477A (zh) * 2017-10-12 2019-04-19 西南科技大学 一种基于lidar数据的dem生成算法
CN109191581A (zh) * 2018-08-29 2019-01-11 李东帅 一种基于bim技术的道路工程土方量高精度计算方法
CN110285792A (zh) * 2019-07-02 2019-09-27 山东省交通规划设计院 一种无人机倾斜摄影精细格网土方计量方法
CN110657855A (zh) * 2019-09-23 2020-01-07 广州市中海达测绘仪器有限公司 土方量测量方法、装置、计算机设备和存储介质
CN111325837A (zh) * 2020-01-23 2020-06-23 江西理工大学 一种基于地面三维激光点云的边坡dem生成方法

Non-Patent Citations (6)

* Cited by examiner, † Cited by third party
Title
刘国光: "基于离散地形点土方量计算研究", 《江西建材》 *
孙成科: "基于三维模型的土石方量核算方法", 《中国高新技术企业》 *
席靖智: "基于DEM的工程土方计算方法优化分析及应用研究", 《中国优秀硕士学位论文全文数据库工程科技II辑》 *
张蒙等: "等高线约束的Delaunay三角网在土石方量计算中的应用", 《地理与地理信息科学》 *
强辉强: "基于DTM土方量计算方法优化分析及应用研究", 《西部大开发(土地开发工程研究)》 *
柯晓山 等: "采用不规则三角网插值进行土地整理项目前期平整土方量的计算", 《农业工程学报》 *

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114527461A (zh) * 2021-12-22 2022-05-24 上海振懋机电科技有限公司 一种空间测量装置及其测量方法
CN115482269A (zh) * 2022-09-22 2022-12-16 佳都科技集团股份有限公司 一种土方量计算方法、装置、终端设备以及存储介质
CN118052947A (zh) * 2024-03-28 2024-05-17 四川省公路规划勘察设计研究院有限公司 基于大数据的三维地理模型建立方法及设备

Also Published As

Publication number Publication date
CN112906124B (zh) 2023-06-13

Similar Documents

Publication Publication Date Title
CN112906124A (zh) 一种土方量计算方法、装置、设备和介质
CN103871102B (zh) 一种基于高程点和道路轮廓面的道路三维精细建模方法
CN111784840B (zh) 基于矢量数据自动分割lod层级三维数据单体化方法及***
Zhu et al. Seed point selection method for triangle constrained image matching propagation
CN111323788B (zh) 建筑物变化的监测方法、装置及计算机设备
CN114332291A (zh) 一种倾斜摄影模型建筑物外轮廓规则提取方法
Tseng et al. Extraction of building boundary lines from airborne LiDAR point clouds
EP2150940A2 (en) Geospatial modeling system providing data thinning of geospatial data points and related methods
CN111667569A (zh) 一种基于Rhino及Grasshopper的三维实景土方可视化精准测算方法
JP4619504B2 (ja) ３次元デジタル地図作成装置
CN115082699A (zh) 一种轮廓形状提取方法、装置、电子设备及存储介质
CN114332134A (zh) 一种基于密集点云的建筑物立面提取方法和装置
CN114241155A (zh) 一种基于车载激光点云数据的城市树木三维可视化方法
KR101079475B1 (ko) 포인트 클라우드 필터링을 이용한 3차원 도시공간정보 구축 시스템
CN112687007B (zh) 一种基于lod技术的立体网格图生成方法
CN109064482A (zh) 一种在三维倾斜摄影数据场景中自动获取房屋轮廓的方法及装置
CN111583406A (zh) 杆塔脚基点坐标计算方法、装置及终端设备
CN116246069B (zh) 自适应地形点云滤波的方法、装置、智能终端及存储介质
CN112884884A (zh) 一种候选区域生成方法及***
CN112687005A (zh) 一种基于三维重建的煤堆体积测量方法
CN116012613B (zh) 一种基于激光点云的露天矿的土方量变化测算方法和***
CN108320324B (zh) 一种城市变化监测方法及***
CN112700538B (zh) Lod生成方法及***
CN113160419B (zh) 一种建筑物立面模型建立方法和装置
CN115713548A (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