CN114019532A - 工程进度核查方法及设备 - Google Patents

Abstract

本发明提供一种工程进度核查方法及设备，方法包括：获取被测物在第一时刻对应的激光点云数据，并根据激光点云数据确定被测物在第一时刻的第一体积；获取被测物在第二时刻的第二体积，其中第二时刻在第一时刻之前；根据第一体积和第二体积确定被测物在目标时间段内的实际体积变化量，其中目标时间段的起始时刻是第二时刻，终止时刻是第一时刻；获取目标时间段对应的上报体积变化量，并将实际体积变化量与上报体积变化量进行比较，生成工程进度核查结果，可以确定上报的体积变化量，即确定工程进度是否准确，实现工程进度的核查，提高上报的体积变化量的可靠性，即提高上报的工程进度的可靠性。

Description

工程进度核查方法及设备

技术领域

本发明实施例涉及工程技术领域，尤其涉及一种工程进度核查方法及设备。

背景技术

土石方量是各项土石方工程量(挖方工程量、填方工程量)之总和。目前，在工程(例如，高速公路建设工程)施工过程中，经常需要人工测量某段时期内的土石方量的变化量，然后将该段时期内的土石方的变化量输入到终端，即上报该段时期的土石方的用量。

现有技术中，在工程(例如，高速公路建设工程)施工过程中，经常需要人工测量某段时期内的土石方量的变化量，然后将该段时期内的土石方的变化量输入到终端，即上报该段时期的土石方的用量。

然而，发明人发现：由于是人工上报土石方的用量，容易出现误报(即上报的土石方与用量与实际用量不符)等问题，造成上报的土石方的用量可靠性较低。

发明内容

本发明实施例提供一种工程进度核查方法及设备，以解决现有技术中可靠性较低的问题。

第一方面，本发明实施例提供一种工程进度核查方法，包括：

获取被测物在第一时刻对应的激光点云数据，并根据所述激光点云数据确定所述被测物在所述第一时刻的第一体积；

获取所述被测物在第二时刻的第二体积，其中所述第二时刻在所述第一时刻之前；

根据所述第一体积和所述第二体积确定所述被测物在目标时间段内的实际体积变化量，其中所述目标时间段的起始时刻是所述第二时刻，终止时刻是所述第一时刻；

获取所述目标时间段对应的上报体积变化量，并将所述实际体积变化量与上报体积变化量进行比较，生成工程进度核查结果。

第二方面，本发明实施例提供一种工程进度核查设备，包括：

数据获取模块，用于获取被测物在第一时刻对应的激光点云数据，并根据所述激光点云数据确定所述被测物在所述第一时刻的第一体积；

处理模块，用于获取所述被测物在第二时刻的第二体积，其中所述第二时刻在所述第一时刻之前；

所述处理模块，还用于根据所述第一体积和所述第二体积确定所述被测物在目标时间段内的实际体积变化量，其中所述目标时间段的起始时刻是所述第二时刻，终止时刻是所述第一时刻；

所述处理模块，还用于获取所述目标时间段对应的上报体积变化量，并将所述实际体积变化量与上报体积变化量进行比较，生成工程进度核查结果。

第三方面，本发明实施例提供一种电子设备，包括：至少一个处理器和存储器；

所述存储器存储计算机执行指令；

所述至少一个处理器执行所述存储器存储的计算机执行指令，使得所述至少一个处理器执行如第一方面任一项所述的工程进度核查方法。

第四方面，本发明实施例提供一种计算机可读存储介质，所述计算机可读存储介质中存储有计算机执行指令，当处理器执行所述计算机执行指令时，实现如上第一方面以及第一方面各种可能的设计所述的工程进度核查方法。

本发明实施例提供的工程进度核查方法及设备，该方法通过获取被测物在第一时刻对应的激光点云数据，并根据所述激光点云数据确定所述被测物在所述第一时刻的第一体积，获取所述被测物在第二时刻的第二体积，其中所述第二时刻在所述第一时刻之前，根据所述第一体积和所述第二体积确定所述被测物在目标时间段内的实际体积变化量，其中所述目标时间段的起始时刻是所述第二时刻，终止时刻是所述第一时刻，获取所述目标时间段对应的上报体积变化量，并将所述实际体积变化量与上报体积变化量进行比较，生成工程进度核查结果。本发明实施例根据被测物在第一时刻对应的激光点云数据计算被测物在第一时刻的体积，然后根据被测物在第一时刻的体积和被测物在第二时刻的体积确定该被测物确定被测物在目标时间段内的实际用量，即确定被测物在目标时间段内的实际体积变化量，实现体积变化量的自动检测，且将计算得到的实际体积变化量与人工上报的上报体积变化量进行比较以确定上报体积变化量是否与实际体积变化量相符，从而确定上报体积变化量是否准确，实现工程进度的核查，即实现体积变化量的核查，提高上报体积变化量的可靠性，从而提高上报的工程进度的可靠性。

附图说明

为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作一简单地介绍，显而易见地，下面描述中的附图是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动性的前提下，还可以根据这些附图获得其他的附图。

图1为本发明实施例提供的工程进度核查系统的示意图；

图2为本发明实施例提供的工程进度核查方法的流程示意图一；

图3为本发明实施例提供的工程进度核查方法的流程示意图二；

图4为本发明实施例提供的被测物对应的底部区域的示意图；

图5为本发明实施例提供的工程进度核查方法的流程示意图三；

图6为本发明实施例提供的工程进度核查设备的结构示意图；

图7为本发明实施例提供的电子设备的硬件结构示意图。

具体实施方式

为使本发明实施例的目的、技术方案和优点更加清楚，下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。

现有技术中，在工程施工过程中，经常需要人工测量某段时期内的土石方量的变化量，即测量某段时期内土石方的体积变化量，然后上报该段时期的土石方的体积变化量，即土石方的用量。但由于是人工上报土石方的用量，容易出现误报等问题，造成上报的土石方的用量可靠性较低。

因此，针对上述问题，本发明的技术构思是获取土石方在目标时间段中的起始时刻对应的激光点云数据，并获取土石方在目标时间段中的终止时刻对应的激光点云数据，然后按照预设自动化算法，对土石方在起始时刻对应的点云数据进行体积计算，得到土石方在起始时刻的体积，同理，确定土石方在终止时刻的体积，计算土石方在终止时刻的体积与土石方在起始时刻的体积的差值，得到土石方在目标时间段内的体积的实际变化量，即得到土石方在目标时间段内的土石方量的实际变化量，实现土石方量的变化量的自动计算，然后自动将土石方在目标时间段内的土石方量的实际变化量与上报的土石方在目标时间段内的土石方量的变化量在进行比较，生成核查报告，以表明上报的土石方在目标时间段内的土石方量的变化量是否存在误报、错报等问题，可以有效的验证上报的土石方量的变化情况，提高上报的土石方量的可靠性，从而起到土石方量进度核查的作用。

下面以具体地实施例对本发明的技术方案进行详细说明。下面这几个具体的实施例可以相互结合，对于相同或相似的概念或过程可能在某些实施例不再赘述。

图1为本发明实施例提供的工程进度核查系统的架构示意图，如图1所示，该系统包括激光雷达装置101、电子设备102和被测物103。激光雷达装置101可以周期性地对被测物进行激光扫描，得到该被测物在当前时刻的激光点云数据，并将该激光数据发送给电子设备102，以供电子设备103进行相关体积计算。

其中，电子设备102可以为服务器或其它具备处理数据能力的设备。被测物可以为土石方，也可以为煤炭堆等呈堆积放置的物体。

需要说明的是，上述被测物103的形状仅为一种示例，被测物可以呈任意形状，只需呈堆积放置即可，本发明实施例并不以此为限。

图2为本发明实施例提供的工程进度核查方法的流程示意图一，本实施例中的执行主体可以电脑、服务器等电子设备，如图2所示，该方法包括：

S201、获取被测物在第一时刻对应的激光点云数据，并根据激光点云数据确定被测物在第一时刻的第一体积。

在本实施例中，电子设备周期性的获取被测物的激光点云数据。电子设备可以通过多种方式获取被测物的激光点云数据，例如，由激光扫描装置周期性地对该被测物进行激光扫描，得到相应的激光点云数据，并将该激光点云数据发送给该电子设备，以使该电子设备进行相关体积计算。

可选的，激光雷达设备可以集成在该电子设备上，也可以为独立的外部设备，用户可以根据实际情况进行设置。激光雷达设备可以每隔预设时间向电子设备发送被测物的激光点云数据，例如，每隔30天发送一次。该预设间隔可以根据实际情况进行合理设置，在此，不该预设时间进行限定。

在本实施例中，第一时刻可以为当前时刻。假设激光扫描装置每隔30天向电子设备发送一次激光点云数据，电子设备每次接收到的激光点云数据均为被测物当前的激光点云数据，即为被测物在当前时刻对应的激光点云数据，均按照本实施例的后续流程进行执行。

其中，激光点云数据包括激光点的坐标值。

在本实施例中，在得到被测物在第一时刻对应的激光点云数据后，根据该被测物在第一时刻对应的激光点云数据确定该被测物在第一时刻的体积，即确定被测物当前的体积，并将该体积作为第一体积，实现被测物在某一时刻的体积的自动检测计算。

S202、获取被测物在第二时刻的第二体积，其中第二时刻在第一时刻之前。

S203、根据第一体积和第二体积确定被测物在目标时间段内的实际体积变化量，其中目标时间段的起始时刻是第二时刻，终止时刻是第一时刻。

在本实施例中，第二时刻在第一时刻之前，其可以为第一时刻的上一时刻。电子设备在获取到被测物在当前时刻的体积后，还需确定被测物在上一时刻的体积，以利用被测物在上一时刻的体积确定从第二时刻开始截止到第一时刻这段时期内的体积变化量。

其中，被测物在第二时刻的第二体积为预先存储的，即预先存储在指定位置中。该第二体积是电子设备根据被测物在第二时刻对应的激光点云数据确定的。

可选的，指定位置可以为数据库，本地硬盘或其它用户指定存储数据的位置，在此，不对指定位置进行限定。

具体的，在确定体积变化量时，根据被测物在第一时刻的体积以及被测物在第二时刻的体积确定该被测物在目标时间段内的实际体积变化量，该目标时间段便为从第二时刻开始截止到第一时刻的这段时期，该实际体积变化量表明被测物在目标时间段内的实际用量，实现被测物在目标时间段内的用量的自动检测。

S204、获取目标时间段对应的上报体积变化量，并将实际体积变化量与上报体积变化量进行比较，生成工程进度核查结果。

在本实施例中，获取被测物在目标时间段对应的上报土石方变化量，该上报体积变化量为相关人员人工测量并上报的被测物的体积变化量，即被测物的用量。将实际体积变化量与上报体积变化量进行比较，生成工程进度核查结果，从而可以根据该工程进度核查结果确定在目标时间段内的被测物的用量是否符合实际用量。

其中，工程进度核查结果为上报的目标时间段内的被测物的用量符合实际，或者为上报的目标时间段内的被测物的用量不符合实际。

具体地，在将该实际体积变化量与上报体积变化量进行比较时，计算该实际变化量与上报体积变化量的差值，若该差值在第一预设差值范围内，表明上报的目标时间段内的体积的体积变化量符合实际，则确定工程进度核查结果为上报的目标时间段内的被测物的用量符合实际，否则，则确定工程进度核查结果为上报的目标时间段内的被测物的用量不符合实际。

其中，第一预设差值范围可以根据实际需求进行设置，在此，不对其进行限定。

另外，在生成工程进度核查结果后，可以将其从电子设备中导出，电子设备也可以发送给需要该核查报告的终端。

在本实施例中，将实际体积变化量与上报体积变化量进行比较，自动生成相应的工程进度核查结果，该工程进度核查结果表明上报的在目标时间段内的土石方的用量是否符合实际用量，可以避免出现弄虚作假的现象，提高上报体积变化量的可靠性，即提高上报的工程进度的可靠性，无需人工进行核查，提高被测物的用量的核查效率。

在本实施例中，根据被测物在第一时刻对应的激光点云数据计算被测物在第一时刻的体积，然后根据被测物在第一时刻的体积和被测物在第二时刻的体积确定该被测物确定被测物在目标时间段内的实际用量，即确定被测物在目标时间段内的实际体积变化量，实现体积变化量的自动检测，无需人工进行计算，从而也不会出现由于人工导致体积变化量的计算效率低以及计算结果可靠性较低的问题。且将计算得到的实际体积变化量与人工上报的上报体积变化量进行比较以确定上报体积变化量是否与实际体积变化量相符，从而确定上报体积变化量是否准确，实现工程进度的核查，即提高上报体积变化量的可靠性，从而保证上报的工程进度的可靠性。

另外，在确定被测物在第一时刻的体积时，可以根据该被测物对应的激光点云数据将该被测物拆分成若干个长方体，分别计算各长方体的体积，然后对其进行累加求和，实现被测物在第一时刻的体积的快速准确获取，下面采用详细的实施例对确定被测物在第一时刻的体积的具体过程进行详细地说明。

图3为本发明实施例提供的工程进度核查方法的流程示意图二，本实施例在图2实施例的基础上，对确定被测物在第一时刻的体积的具体实现过程进行了详细说明，如图3所示，该方法包括：

S301、获取被测物在第一时刻对应的激光点云数据。

在本实施例中，激光扫描装置在对被测物进行激光扫描时，实际是对被测区域进行激光扫描，该被测区域为被测物所在的区域。因此，扫描得到的激光点云数据可能会包括很多噪点数据，例如植被、建筑、车辆、以及被测区域的点云数据等等。电子设备在接收到激光扫描装置发送的激光点云数据后，对该激光点云数据进行分类，然后将不属于被测物的激光点云数据剔除，实现噪点数据的剔除。

可选的，被测区域可以是人为设定的，也可是由电子设备进行划定和裁切得到的，以尽量减少噪点，即使被测区域尽量仅包括被测物。

S302、获取被测物对应的底部区域，其中底部区域包括多个网格。

在本实施例中，获取被测物对应的底部区域，该底部区域可以覆盖在当前时刻及之前时刻的被测物的投影区域，即包括被测物在各被测时刻对应的激光点云数据中的所有激光点的投影。

其中，被测时刻为激光扫描该被测物，得到该被测物对应的激光点云数据的时刻，即激光扫描装置将激光点云数据上报给电子设备的时刻，由于激光扫描装置是周期性上报，因此，被测时刻存在多个。该被测时刻包括上述的第一时刻和第二时刻。

另外，还可以对被测物对应的激光点云数据进行可视化以清楚地展示被测物变化，用户可以通过选择演示方式，然后按照用户选择的可视化方式进行演示。

其中，可视化方式包括单批次点云演示方式、多批次点云多窗口联动展示方式、多批次点云的时间轴演示方式等等。

可选的，预设三维坐标系由x轴，y轴和z轴组成。被测物对应的底部区域可以根据电子设备已经获取到的所有激光点云进行确定，即获取被测物对应的所有激光点云数据，其中所有激光点云数据包括被测物在各被测时刻对应的激光点云数据。查找所有激光点云数据中的第一最大坐标值和第一最小坐标值，其中第一最大坐标值为所有激光点云数据在预设三维坐标系中的x轴上的最大值，第一最小坐标值为所有激光点云数据在x轴上的最小值。查找所有激光点云数据中的第二最大坐标值和第二最小坐标值，其中第二最大坐标值为所有激光点云数据在预设三维坐标系中的y轴上的最大值，第二最小坐标值为所有激光点云数据在y轴上的最小值。根据第一最大坐标值、第一最小坐标值、第二最大坐标值和第二最小坐标值确定被测物对应的底部区域。

在本实施例中，电子设备获取其接收到的被测物对应的所有激光点云数据，被测物对应的所有激光点云数据包括多个激光点的坐标值，该坐标值包括激光点在z轴上的坐标值，激光点在x轴上的坐标值和激光点在y轴上的坐标值。查找所有激光点云数据在x轴上的最大值，即查找该多个激光点中在x轴上的坐标值最大的激光点，并将该激光点在x轴上的坐标值作为第一最大坐标值，以及查找所有激光点云数据中在x轴上的最小值，即查找该多个激光点中在x轴上的坐标值最小的激光点，并将该激光点在x轴上的坐标值作为第一最小坐标值。

同理，按照查找第一最大坐标值和第一最小坐标值的方式查找第二最大坐标值及第二最小坐标值。

以一个具体应用场景为例，所有激光点云数据包括5个激光点的坐标值，分别为激光点1的坐标值，激光点2的坐标值，激光点3的坐标值，激光点4的坐标值和激光点5的坐标值，激光点3的坐标值中的x值最大，则将该x值作为第一最大坐标值，激光点1的坐标值中的x值最小，则将该x值作为第一最小坐标值，激光点4的坐标值中的y值最大，则将该y值作为第二最大坐标值，激光点5的坐标值中的y值最小，则将该y值作为第二最小坐标值。

可选的，根据第一最大坐标值、第一最小坐标值、第二最大坐标值和第二最小坐标值确定被测物对应的底部区域，包括：根据第一最大坐标值、第一最小坐标值、第二最大坐标值和第二最小坐标值绘制矩形区域，按照预设网格边长对矩形区域进行划分，得到被测物对应的底部区域。

在本实施例中，根据第一最大坐标值、第一最小坐标值、第二最大坐标值和第二最小坐标值绘制矩形区域，即分别将第一最大坐标值对应的激光点、第一最小坐标值对应的激光点、第二最大坐标值对应的激光点和第二最小坐标值对应激光点作为顶点，然后依次将各个顶点进行连接，得到该待测物对应的矩形区域。由于是将被测物拆分计算体积，因此，还需按照预设网格边长对该矩形区域进行划分，以将矩形区域划分为多个网格，该网格为正方形，边长为该预设网格边长，例如，对图4中矩形区域进行划分得到相应的底部区域，该底部区域包括网格G1、网格G2、网格G3和网格G4。

另外，在按照预设网格边长对该矩形区域进行划分时，在得到第一最大坐标值、第一最小坐标值、第二最大坐标值和第二最小坐标值后，还可以按照以下公式处理得到的坐标值：

x0＝floor(x_min/step)*step，y0＝floor(y_min/step)*step；

nx＝ceil((x_max-x0)/step)；ny＝ceil((y_max-y0)/step)。

其中，x0为处理后的第一最小坐标值，x_min为第一最小坐标值，step为预设网格边长，y0为处理后的第二最小坐标值，y_min为第二最小坐标值，x_max为第一最大坐标值，y_max为第二最大坐标值。

可选的，预设网格边长可以按照实际情况进行设定，预设网格边长值越小，计算得到的体积越准确。

另外，为了方便被测物的体积计算以及提高体积计算的准确性，可以将当前接收的被测物在某一时刻对应的激光点云数据与之前已经接收的激光点云数据统一基准和配准，即进行三维激光点云数据拼接，以将当前接收的被测物在某一时刻对应的激光点云数据与之前已经接收的激光点云数据统一到同一坐标系下。

在将激光点云数据统一到同一坐标系下时，可以利用人为设定的固定特征标记，并通过旋转、平移等操作，把在不同坐标系下测得的激光点云数据进行坐标转换，以将其统一到统一坐标系下。

现有在采集被测物的点云数据时，一般是通过摄像机获取被测物的图像以根据该图像重建三维点云结构；或者是通过无人机获取被测物的无人机点云数据以及通过三维激光点云获取被测物的三维激光点云数据，将三维激光点云数据和无人机点云数据转化为为相同格式，然后以无人机点云数据为参考系，三维激光点云数据为移动系，基于三维点云数据的RGB颜色类型选取多个同名点，利用全局匹配方式，将两者的点云数据建立相关性匹配。

而本申请是直接采用激光雷达设备来采集被测物在不同时刻对应的激光点云数据，将然后利用固定特征标记将被测物在不同时刻对应的激光点云数据统一到同一坐标系下，无需利用RGB颜色值选同名点，更无需将地面点云与无人机点云匹配，实现点云数据的快速准确获取。

S303、对于每个网格，获取第一激光点及其对应的坐标值，并根据第一激光点对应的坐标值以及预设网格面积确定网格在第一时刻对应的体积，其中第一激光点的投影在网格中。

在本实施例中，由于被测物对应的底部区域包括多个网格，在计算该被测物的体积时，相当于将该被测物拆分为若干个长方体，每个长方体的底部为该底部区域中的一个网格，然后计算各长方体的体积，即在计算被测物在某一时刻的体积时，分别计算该被测物对应的底部区域中每个网格对应的体积。

具体地，在计算底部区域中的网格对应体积时，获取投影落在该网格上的激光点，该激光点的x值在该网格对应的x值范围内，该激光点的y值在该网格对应的y值范围内，且该激光点为被测物在第一时刻对应的激光点云数据中的激光点，将获取的激光点作为第一激光点。然后利用该第一激光点的坐标值中的z值，即高程值以及该网格对应的网格面积确定该网格在第一时刻的体积。

可选的，根据第一激光点对应的坐标值以及预设网格面积确定网格在第一时刻对应的体积，包括：

获取各第一激光点在预设三维坐标系中z轴上的坐标值。

通过各第一激光点在z轴上的坐标值计算得到处理值，并将处理值作为网格在第一时刻对应的第一高程值。

获取网格在初始时刻对应的初始高程值。

获取第一高程值与初始高程值的差值，并根据差值与网格对应的预设网格面积，得到网格在第一时刻对应的体积。

具体地，根据网格对应的各第一激光点在z轴上的坐标值计算得到处理值，并将该处理值作为该网格在第一时刻对应的第一高程值，即该网格在第一时刻对应的高度。

其中，第一激光点在z轴上的坐标值实际为第一激光点对应的高程值，即第一激光点的高度值。

由于在初始时刻，被测物的体积可能并不是零，即在被测区域已经存在被测物，则在计算被测物在某一时刻的体积时，还需减去被测物在初始时刻的体积，因此，在计算网格在第一时刻对应的体积时，还需减去该网格在初始时刻对应的体积，则获取该网格在初始时刻对应的初始高程值，将该第一高程值减去该初始高程值，得到相应的差值，然后将该差值与该网格对应的网格面积的乘积作为该网格在第一时刻对应的体积。

其中，处理值可以为平均值，方差，标准差等常用统计指标。

以一个具体应用场景为例，当处理值为平均值时，如图4所示，在计算网格G1对应的体积时，获取被测物在第一时刻对应的激光点云数据中的投影落在网格G1上的激光点，其分别为激光点P1和激光点P2，即激光点P1和激光点P2在网格G1的垂直上方，其为第一激光点，获取激光点P1在z轴上的坐标值，即获取激光点P1对应的坐标值中的z值，获取激光点P2在z轴上的坐标值，即获取激光点P2对应的坐标值中的z值，计算激光点P1对应的z值和激光点P2对应的z值的平均值，并将该平均值作为该网格在第一时刻对应的第一高程值h1，获取网格G1在初始时刻对应的高度值，即网格G1的初始高程值h2，计算h1与h2的差值，并将该差值与网格G1对应的网格面积的乘积作为网格G1在第一时刻的对应的体积，然后按照此过程，分别计算网格G2对应的体积、网格G3对应的体积和网格G4对应的体积。

另外，网格在初始时刻对应的初始高程值可以是相关人员预先测量得到的，也可以是按照确定网格在第一时刻对应的第一高程值的过程确定的。

S304、根据各个网格在第一时刻对应的体积，得到被测物在第一时刻的第一体积。

在本实施例中，在得到底部区域中的各个网格在第一时刻对应的体积后，利用各网格在第一时刻对应的体积确定该被测物在第一时刻的体积，即确定第一体积。

可选的，在利用各网格在第一时刻对应的体积确定该被测物在第一时刻的体积时，获取各个网格在第一时刻对应的体积的和，并将其确定为所述被测物在所述第一时刻的第一体积，即计算底部区域中的所有网格在第一时刻对应的体积的和，并将该和作为该被测物整体在第一时刻的体积。

承接上述应用场景，在得到网格G1对应的体积网格，G2对应的体积，网格G3对应的体积和网格G4对应的体积后，将网格G1对应的体积网格，G2对应的体积，网格G3对应的体积和网格G4对应的体积进行累加，得到被测物在第一时刻的体积。

另外，在得到被测物在第一时刻的体积后，也可以对其进行上报或导出，例如，将该体积发送给目标终端，以使目标终端对应的相关人员可以获知被测物在第一时刻的体积。

S305、获取被测物在第二时刻的第二体积，其中第二时刻在第一时刻之前。

S306、根据第一体积和第二体积确定被测物在目标时间段内的实际体积变化量，其中目标时间段的起始时刻是第二时刻，终止时刻是第一时刻。

S307、获取目标时间段对应的上报体积变化量，并将实际体积变化量与上报体积变化量进行比较，生成工程进度核查结果。

本实施例的S305-S307的具体实施方式，与上述实施例中的S202至S204类似，此处不再赘述。

现有在计算被测物的体积时，一般是通过建模来进行计算，例如对获取的图像进行特征点提取与匹配、三维重建，在三维点云中加入特征点间的实际距离与在相对坐标系与真实世界的距离比例作为尺度信息，利用特征点三维坐标，构建Delaunay三角网模型，利用Delaunay三角网模型计算被测物的体积。又例如，从目标区域中划分出RTK(Real-timekinematic)待测区域和PPK(Post processed kinematic)待测区域，在目标区域开展飞行设备航摄作业，得到倾斜摄影图像，根据倾斜摄影图像进行点云处理，得到点云数据，根据点云数据建模，例如，构建dem数字高程模型，得到初步挖填土石方工程量，分别在待测区域扫描，得到植被覆盖区，并获取植被覆盖区的植被总方量，根据初步挖填土石方工程量和植被总方量的差值得到最终挖填土石方工程量；又例如在施工期间不同阶段对观测对象所在测区进行扫描，得到测区对应的点云数据，将施工期间不同阶段的点云数据构造不规则三角网模型TIN后，选中已构造的不规则三角网模型TIN投影到相同基准面上进行土方计算，以得到土方变化量。

而本申请在计算被测物的体积时，直接对被测物对应的激光点云数据所包括的激光点的坐标值以及预设网格边长将被测物拆分成若干个小长方体，分别计算每个长方体的体积，然后对将每个长方体的体积进行相加，得到该被测物的体积，检测方法简单，无需提取特征点以及计算特征点之间的实际距离，更无需构建TIN，DEM，Delaunay等模型，由于在进行建模时，需要相关人员使用相应的建模工具进行建模，而本申请无需进行建模，因此可以大幅度提高自动化程度，并且由于无需涉及到大量计算，可以提高工程进度核查的效率以及准确率。

另外，本申请可以通过设定预设网格边长的值来调整被测物的体积的计算精度，当对被测物的体积的计算精度要求较高时，设定预设网格变长的值较小，而当对被测物的体积的计算精度要求较低时，设定预设网格变长的值较大以提高体积的计算效率。

具体地，当被测物为土石方时，在得到土石方在目标时间段内的实际体积变化量后，还可以将其与标准土石方变化量进行比较，以生成相应的施工工程土石方进度核查结果，下面采用详细的实施例对生成施工工程土石方进度核查结果的具体过程进行详细地说明。

图5为本发明实施例提供的工程进度核查方法的流程示意图三，本实施例在图2实施例的基础上，对生成施工工程土石方进度核查结果的的具体实现过程进行了详细说明，如图5所示，该方法包括：

S501、获取被测物在第一时刻对应的激光点云数据，并根据激光点云数据确定被测物在第一时刻的第一体积。

S502、获取被测物在第二时刻的第二体积，其中第二时刻在第一时刻之前。

S503、根据第一体积和第二体积确定被测物在目标时间段内的实际体积变化量，其中目标时间段的起始时刻是第二时刻，终止时刻是第一时刻。

本实施例的S501至S503的具体实施方式，与上述实施例中的S201至S203类似，此处不再赘述。

S504、获取目标时间段对应的标准土石方变化量。

S505、将实际体积变化量与标准土石方变化量进行比较，生成施工工程土石方进度核查结果。

在本实施例中，获取土石方在目标时间段对应的标准土石方变化量，该标准土石方变化量为相关施工计划中制定的土石方量的计划变化量。将实际体积变化量与标准土石方变化量进行比较，生成施工工程土石方进度核查结果，从而可以根据该土石方核查报告确定在目标时间段内的土石方的用量是否符合计划。

其中，施工工程土石方进度核查结果为目标时间段内的土石方的用量符合计划，或者为目标时间段内的土石方的用量不符合计划。

另外，在将该实际体积变化量与标准土石方变化量进行比较时，计算该实际变化量与标准土石方变化量的差值，若该差值在第二预设差值范围内，表明目标时间段内的土石方的实际体积变化量符合计划，则确定施工工程土石方进度核查结果为目标时间段内的土石方的用量符合计划，否则，则确定施工工程土石方进度核查结果为目标时间段内的土石方的用量不符合计划。

其中，第二预设差值范围和第一预设差值范围可以相同也可以不同，可以根据实际需求进行设置，在此，不对其进行限定。

另外，可选的，根据施工工程土石方进度核查结果和/或工程进度核查结果生成相应的土石方核查报告，以使相关人员可以根据土石方核查报告确定上报的土石方量的变化量是否存在误报、错报以及不符合计划等问题，可以有效的验证上报的土石方量的变化情况，提高上报的土石方量的可靠性，从而起到土石方量进度核查的作用。

在本实施例中，本申请应用于公路施工领域，在公路施工过程中，需要上报人员上报土石方的用量，即土石方量的变化量，因此，需要验证上报人员上报的土石方的用量是否存在误报、错报等问题。在此基础上，本申请根据土石方在不同时刻的激光点云数据计算土石方量的变化量，然后根据土石方量的变化量生成相应的土石方核查报告，以使相关人员根据土石方核查报告确定上报人员上报的土石方量变化量是否存在误报、错报等问题，起到土石方量进度核查的作用。

在本实施例中，将实际体积变化量与标准土石方变化量进行比较，自动生成相应的施工工程土石方进度核查结果，该施工工程土石方进度核查结果表明在目标时间段内的土石方的用量是否符合计划，以使相关人员根据施工工程土石方进度核查结果及时调整土石方的用量，避免出现土石方用量过多或过少的问题，无需人工进行核查，提高土石方的用量的核查效率。

图6为本发明实施例提供的工程进度核查设备的结构示意图。如图6所示，该工程进度核查设备600包括：数据获取模块601和处理模块602。

其中，数据获取模块601，用于获取被测物在第一时刻对应的激光点云数据，并根据激光点云数据确定被测物在第一时刻的第一体积。

处理模块602，还用于获取被测物在第二时刻的第二体积，其中第二时刻在第一时刻之前。

处理模块602，还用于根据第一体积和第二体积确定被测物在目标时间段内的实际体积变化量，其中目标时间段的起始时刻是第二时刻，终止时刻是第一时刻。

处理模块602，还用于获取目标时间段对应的上报体积变化量，并将实际体积变化量与上报体积变化量进行比较，生成工程进度核查结果。

在一种可能的设计中，激光点云数据包括激光点的坐标值。

数据获取模块601，还用于获取被测物对应的底部区域，其中底部区域包括多个网格。

数据获取模块601，还用于对于每个网格，获取第一激光点及其对应的坐标值，并根据第一激光点对应的坐标值以及预设网格面积确定网格在第一时刻对应的体积，其中第一激光点的投影在网格中。

数据获取模块601，还用于根据各个网格在第一时刻对应的体积，得到被测物在第一时刻的第一体积。

在一种可能的设计中，数据获取模块601，还用于获取各第一激光点在预设三维坐标系中z轴上的坐标值。

数据获取模块601，还用于通过各第一激光点在z轴上的坐标值计算得到处理值，并将处理值作为网格在第一时刻对应的第一高程值。

数据获取模块601，还用于获取网格在初始时刻对应的初始高程值。

数据获取模块601，还用于获取第一高程值与初始高程值的差值，并根据差值与网格对应的预设网格面积，得到网格在第一时刻对应的体积。

在一种可能的设计中，数据获取模块601，还用于在获取被测物对应的底部区域之前，获取被测物对应的所有激光点云数据，其中所有激光点云数据包括被测物在各被测时刻对应的激光点云数据，被测时刻为激光扫描被测物的时刻。

数据获取模块601，还用于获取所有激光点云数据中的第一最大坐标值和第一最小坐标值，其中第一最大坐标值为所有激光点云在预设三维坐标系中的x轴上的最大值，第一最小坐标值为所有激光点云在x轴上的最小值。

数据获取模块601，还用于获取所有激光点云数据中的第二最大坐标值和第二最小坐标值，其中第二最大坐标值为所有激光点云在预设三维坐标系中的y轴上的最大值，第二最小坐标值为所有激光点云在y轴上的最小值。

数据获取模块601，还用于根据第一最大坐标值、第一最小坐标值、第二最大坐标值和第二最小坐标值确定被测物对应的底部区域。

在一种可能的设计中，数据获取模块601，还用于根据第一最大坐标值、第一最小坐标值、第二最大坐标值和第二最小坐标值绘制矩形区域。

数据获取模块601，还用于按照预设网格边长对矩形区域进行划分，得到被测物对应的底部区域。

在一种可能的设计中，数据获取模块601，还用于获取各个网格在第一时刻对应的体积的和，并将其确定为被测物在第一时刻的第一体积。

在一种可能的设计中，被测物为土石方。

处理模块602，还用于获取目标时间段对应的标准土石方变化量。

处理模块602，还用于将实际体积变化量与标准土石方变化量进行比较，生成施工工程土石方进度核查结果。

在一种可能的设计中，被测物在第二时刻的第二体积为预先存储的。

本实施例提供的工程进度核查设备，可用于执行上述方法实施例的技术方案，其实现原理和技术效果类似，本实施例此处不再赘述。

图7为本发明实施例提供的电子设备的硬件结构示意图。如图7所示，本实施例的电子设备700包括：至少一个处理器701和存储器702。其中，处理器701、存储器702通过总线703连接。

在具体实现过程中，至少一个处理器701执行存储器702存储的计算机执行指令，使得至少一个处理器701执行上述方法实施例中的工程进度核查方法。

处理器701的具体实现过程可参见上述方法实施例，其实现原理和技术效果类似，本实施例此处不再赘述。

在上述的图7所示的实施例中，应理解，处理器可以是中央处理单元(英文：Central Processing Unit，简称：CPU)，还可以是其他通用处理器、数字信号处理器(英文：Digital Signal Processor，简称：DSP)、专用集成电路(英文：Application SpecificIntegrated Circuit，简称：ASIC)等。通用处理器可以是微处理器或者该处理器也可以是任何常规的处理器等。结合发明所公开的方法的步骤可以直接体现为硬件处理器执行完成，或者用处理器中的硬件及软件模块组合执行完成。

存储器可能包含高速RAM存储器，也可能还包括非易失性存储NVM，例如至少一个磁盘存储器。

总线可以是工业标准体系结构(Industry Standard Architecture，ISA)总线、外部设备互连(Peripheral Component，PCI)总线或扩展工业标准体系结构(ExtendedIndustry Standard Architecture，EISA)总线等。总线可以分为地址总线、数据总线、控制总线等。为便于表示，本申请附图中的总线并不限定仅有一根总线或一种类型的总线。

本发明实施例还提供一种计算机可读存储介质，所述计算机可读存储介质中存储有计算机执行指令，当处理器执行所述计算机执行指令时，实现如上所述的工程进度核查方法。

在本发明所提供的几个实施例中，应该理解到，所揭露的设备和方法，可以通过其它的方式实现。例如，以上所描述的设备实施例仅仅是示意性的，例如，所述模块的划分，仅仅为一种逻辑功能划分，实际实现时可以有另外的划分方式，例如多个模块可以结合或者可以集成到另一个系统，或一些特征可以忽略，或不执行。另一点，所显示或讨论的相互之间的耦合或直接耦合或通信连接可以是通过一些接口，设备或模块的间接耦合或通信连接，可以是电性，机械或其它的形式。

所述作为分离部件说明的模块可以是或者也可以不是物理上分开的，作为模块显示的部件可以是或者也可以不是物理单元，即可以位于一个地方，或者也可以分布到多个网络单元上。可以根据实际的需要选择其中的部分或者全部模块来实现本实施例方案的目的。

另外，在本发明各个实施例中的各功能模块可以集成在一个处理单元中，也可以是各个模块单独物理存在，也可以两个或两个以上模块集成在一个单元中。上述模块成的单元既可以采用硬件的形式实现，也可以采用硬件加软件功能单元的形式实现。

上述以软件功能模块的形式实现的集成的模块，可以存储在一个计算机可读取存储介质中。上述软件功能模块存储在一个存储介质中，包括若干指令用以使得一台计算机设备(可以是个人计算机，服务器，或者网络设备等)或处理器(英文：processor)执行本申请各个实施例所述方法的部分步骤。

应理解，上述处理器可以是中央处理单元(英文：Central Processing Unit，简称：CPU)，还可以是其他通用处理器、数字信号处理器(英文：Digital Signal Processor，简称：DSP)、专用集成电路(英文：Application Specific Integrated Circuit，简称：ASIC)等。通用处理器可以是微处理器或者该处理器也可以是任何常规的处理器等。结合发明所公开的方法的步骤可以直接体现为硬件处理器执行完成，或者用处理器中的硬件及软件模块组合执行完成。

存储器可能包含高速RAM存储器，也可能还包括非易失性存储NVM，例如至少一个磁盘存储器，还可以为U盘、移动硬盘、只读存储器、磁盘或光盘等。

总线可以是工业标准体系结构(Industry Standard Architecture，ISA)总线、外部设备互连(Peripheral Component，PCI)总线或扩展工业标准体系结构(ExtendedIndustry Standard Architecture，EISA)总线等。总线可以分为地址总线、数据总线、控制总线等。为便于表示，本申请附图中的总线并不限定仅有一根总线或一种类型的总线。

上述存储介质可以是由任何类型的易失性或非易失性存储设备或者它们的组合实现，如静态随机存取存储器(SRAM)，电可擦除可编程只读存储器(EEPROM)，可擦除可编程只读存储器(EPROM)，可编程只读存储器(PROM)，只读存储器(ROM)，磁存储器，快闪存储器，磁盘或光盘。存储介质可以是通用或专用计算机能够存取的任何可用介质。

一种示例性的存储介质耦合至处理器，从而使处理器能够从该存储介质读取信息，且可向该存储介质写入信息。当然，存储介质也可以是处理器的组成部分。处理器和存储介质可以位于专用集成电路(Application Specific Integrated Circuits，简称：ASIC)中。当然，处理器和存储介质也可以作为分立组件存在于电子设备或主控设备中。

本领域普通技术人员可以理解：实现上述各方法实施例的全部或部分步骤可以通过程序指令相关的硬件来完成。前述的程序可以存储于一计算机可读取存储介质中。该程序在执行时，执行包括上述各方法实施例的步骤；而前述的存储介质包括：ROM、RAM、磁碟或者光盘等各种可以存储程序代码的介质。

最后应说明的是：以上各实施例仅用以说明本发明的技术方案，而非对其限制；尽管参照前述各实施例对本发明进行了详细的说明，本领域的普通技术人员应当理解：其依然可以对前述各实施例所记载的技术方案进行修改，或者对其中部分或者全部技术特征进行等同替换；而这些修改或者替换，并不使相应技术方案的本质脱离本发明各实施例技术方案的范围。

Claims (10)

1.一种工程进度核查方法，其特征在于，包括：

获取被测物在第一时刻对应的激光点云数据，并根据所述激光点云数据确定所述被测物在所述第一时刻的第一体积；

获取所述被测物在第二时刻的第二体积，其中所述第二时刻在所述第一时刻之前；

根据所述第一体积和所述第二体积确定所述被测物在目标时间段内的实际体积变化量，其中所述目标时间段的起始时刻是所述第二时刻，终止时刻是所述第一时刻；

获取所述目标时间段对应的上报体积变化量，并将所述实际体积变化量与上报体积变化量进行比较，生成工程进度核查结果。

2.根据权利要求1所述的方法，其特征在于，所述激光点云数据包括激光点的坐标值；

所述根据所述激光点云数据确定所述被测物在所述第一时刻的第一体积，包括：

获取所述被测物对应的底部区域，其中所述底部区域包括多个网格；

对于每个网格，获取第一激光点及其对应的坐标值，并根据所述第一激光点对应的坐标值以及预设网格面积确定所述网格在第一时刻对应的体积，其中所述第一激光点的投影在所述网格中；

根据各个所述网格在第一时刻对应的体积，得到所述被测物在所述第一时刻的第一体积。

3.根据权利要求2所述的方法，其特征在于，所述根据所述第一激光点对应的坐标值以及预设网格面积确定所述网格在第一时刻对应的体积，包括：

获取各第一激光点在预设三维坐标系中z轴上的坐标值；

通过各第一激光点在z轴上的坐标值计算得到处理值，并将所述处理值作为所述网格在第一时刻对应的第一高程值；

获取所述网格在初始时刻对应的初始高程值；

获取所述第一高程值与所述初始高程值的差值，并根据所述差值与所述网格对应的预设网格面积，得到所述网格在第一时刻对应的体积。

4.根据权利要求2所述的方法，其特征在于，在所述获取所述被测物对应的底部区域之前，还包括：

获取被测物对应的所有激光点云数据，其中所述所有激光点云数据包括所述被测物在各被测时刻对应的激光点云数据，所述被测时刻为激光扫描所述被测物的时刻；

查找所述所有激光点云数据中的第一最大坐标值和第一最小坐标值，其中所述第一最大坐标值为所述所有激光点云在预设三维坐标系中的x轴上的最大值，所述第一最小坐标值为所述所有激光点云在所述x轴上的最小值；

查找所述所有激光点云数据中的第二最大坐标值和第二最小坐标值，其中所述第二最大坐标值为所述所有激光点云在所述预设三维坐标系中的y轴上的最大值，所述第二最小坐标值为所述所有激光点云在所述y轴上的最小值；

根据所述第一最大坐标值、所述第一最小坐标值、所述第二最大坐标值和所述第二最小坐标值确定所述被测物对应的底部区域。

5.根据权利要求4所述的方法，其特征在于，所述根据所述第一最大坐标值、所述第一最小坐标值、所述第二最大坐标值和所述第二最小坐标值确定所述被测物对应的底部区域，包括：

根据所述第一最大坐标值、所述第一最小坐标值、所述第二最大坐标值和所述第二最小坐标值绘制矩形区域；

按照预设网格边长对所述矩形区域进行划分，得到所述被测物对应的底部区域。

6.根据权利要求2所述的方法，其特征在于，所述根据各个网格在第一时刻对应的体积，得到所述被测物在所述第一时刻的第一体积，包括：

获取各个网格在第一时刻对应的体积的和，并将其确定为所述被测物在所述第一时刻的第一体积。

7.根据权利要求1至6任一项所述的方法，其特征在于，所述被测物为土石方，所述方法还包括：

获取所述目标时间段对应的标准土石方变化量；

将所述实际体积变化量与所述标准土石方变化量进行比较，生成施工工程土石方进度核查结果。

8.根据权利要求1至6任一项所述的方法，其特征在于，所述被测物在第二时刻的第二体积为预先存储的。

9.一种工程进度核查设备，其特征在于，包括：

数据获取模块，用于获取被测物在第一时刻对应的激光点云数据，并根据所述激光点云数据确定所述被测物在所述第一时刻的第一体积；

处理模块，用于获取所述被测物在第二时刻的第二体积，其中所述第二时刻在所述第一时刻之前；

所述处理模块，还用于根据所述第一体积和所述第二体积确定所述被测物在目标时间段内的实际体积变化量，其中所述目标时间段的起始时刻是所述第二时刻，终止时刻是所述第一时刻；

所述处理模块，还用于获取所述目标时间段对应的上报体积变化量，并将所述实际体积变化量与上报体积变化量进行比较，生成工程进度核查结果。

10.一种计算机可读存储介质，其特征在于，所述计算机可读存储介质中存储有计算机执行指令，当处理器执行所述计算机执行指令时，实现如权利要求1至8任一项所述的工程进度核查方法。

Images