Digital construction method based on three-dimensional surveying and mapping of unmanned aerial vehicle

Abstract

The invention relates to the technical field of remote sensing surveying and mapping, in particular to a digital construction method based on three-dimensional surveying and mapping of an unmanned aerial vehicle, which comprises the following steps: importing a digital earth surface model and identifying an obstacle, calling identification information to calculate a path avoidance cost to plan an unmanned aerial vehicle path, and adjusting a flight height in combination with an image definition demand. And analyzing the measurement data to reconstruct a three-dimensional model, analyzing the actual project progress and the size deviation, and obtaining a construction project measurement record. According to the invention, through regional obstacle identification, bypass cost score calculation combined with multiple measurement paths, safety and rationality of route planning are optimized, and a flight height and image definition linkage adjustment mode is utilized, so that dynamic control of imaging quality and construction of construction engineering model components are realized; the construction progress identification and the dimensional deviation detection are realized, the progress judgment accuracy and the space monitoring capability of digital construction are effectively improved, a data basis is provided for project management decision, and the construction quality is improved.

Description

Technical Field

[0001] The present invention relates to the technical field of remote sensing mapping, and particularly to a digital construction method based on three-dimensional mapping by unmanned aerial vehicles. Background Art

[0002] The technical field of remote sensing mapping includes obtaining surface information using aerospace platforms and performing data processing, analysis, and modeling. The core of this technical field lies in obtaining and expressing the spatial positions and morphological characteristics of targets such as terrain, landforms, and buildings through means such as remote sensing images, lidar, and aerial photography. Remote sensing mapping involves multiple sub-links, including sensor platform mounting, flight path planning, image acquisition and correction, three-dimensional modeling, and result expression. Its systematicness is reflected in aspects such as platform diversification, data acquisition automation, processing process standardization, and result expression multi-dimensionalization, and it is applied in multiple aspects such as geographic information systems, urban and rural construction, resource and environmental monitoring, and engineering surveying.

[0003] Among them, a digital construction method based on three-dimensional mapping by unmanned aerial vehicles refers to using an unmanned aerial vehicle platform as a data acquisition carrier, combining with three-dimensional real-scene modeling of the construction area to achieve data visualization and model expression throughout the construction process, and is used to assist construction planning, progress management, and on-site control. This patent theme mainly covers aspects such as airspace approval at the construction site, landing point layout, obstacle identification, terrain-following flight path design, setting of flight height and overlap parameters, etc. By importing DWG drawings or DSM models, terrain analysis and task zoning are completed, images or point cloud data with spatial reference attributes are collected, and later model splicing is achieved through multi-region overlap zone settings, providing three-dimensional basic data support for the digital construction process. The method is based on three-dimensional modeling by unmanned aerial vehicles and runs through stages such as pre-construction preparation, mid-term monitoring, and post-construction re-survey, and is an important component for realizing digital management of construction information.

[0004] In terms of operation logic, traditional unmanned aerial vehicle three-dimensional mapping and digital construction technologies mainly focus on static modeling or phased data splicing, lacking a fine division of the dynamic structure of regional elevation and an abnormal identification process, and it is difficult to actively construct a continuous terrain risk identification model. In terms of path planning, a preset path or a flight path construction method based on simple obstacle avoidance rules is adopted, ignoring the trade-off analysis between avoidance cost and path burden, resulting in unreasonable path selection. In the data acquisition stage, fixed-height acquisition is mainly used, and the imaging quality and focal length changes of images cannot be sensed in real time, resulting in blurred images or unstable measurement accuracy. In terms of construction progress judgment, it relies on phased model comparison or manual evaluation methods, lacking an automatic identification mechanism for component-level spatial errors and direction changes, resulting in delayed progress judgment or missed detection of component installation errors, affecting on-site control and post-construction re-survey efficiency, and leading to repeated measurements, data offset accumulation, and path execution conflicts in actual construction, reducing the data transfer efficiency of digital construction. Summary of the Invention

[0005] In order to solve the technical problems existing in the prior art, an embodiment of the present invention provides a digital construction method based on UAV three-dimensional surveying and mapping. The technical solution is as follows: In order to achieve the above object, the present invention adopts the following technical solution. A digital construction method based on UAV three-dimensional surveying and mapping includes the following steps: S1: Import the regional digital surface model, divide the construction area into multiple grids and identify the elevation peak points, calculate the elevation offset of the peak points by combining the elevation values of adjacent positions of the peak points, identify multiple obstacles in the area, and generate obstacle identification information; S2: Call the obstacle identification information, extract the spatial position and geometric shape of the obstacles, calculate the avoidance cost score by analyzing the increased path length, flight height adjustment amount and path turning angle in each UAV path for avoidance flight, and obtain the UAV measurement path; S3: Obtain the UAV measurement path, dynamically adjust the flight height control parameters of the UAV according to the acquisition accuracy requirements and in combination with the real-time image clarity, and collect the image information and lidar point cloud data of the construction area in real time to generate a regional imaging data set; S4: According to the regional imaging data set, use the flight height value to analyze the image and lidar point cloud data, reconstruct the construction engineering model of the construction area, compare it with the building information model of the target project, calculate the completion degree of the project in real time and identify the project progress, and generate construction progress information.

[0006] As a further solution of the present invention, the obstacle identification information is specifically an obstacle spatial distribution layer, obstacle geometric attribute data, and an obstacle unique identification code. The UAV measurement path includes a path segment start and end coordinate set, a path heading angle sequence, and a path corresponding obstacle avoidance strategy parameter set. The regional imaging data set includes an image resolution label, a point cloud coverage density map, and an image positioning coordinate set. The construction progress information specifically refers to the block completion status, the corresponding time identifier of the process node, and the construction completion ratio.

[0007] As a further solution of the present invention, the step of importing the regional digital surface model, dividing the construction area into multiple grids and identifying the elevation peak points, calculating the elevation offset of the peak points by combining the elevation values of adjacent positions of the peak points, identifying multiple obstacles in the area, and generating obstacle identification information is specifically as follows: S101: Import the regional digital surface model, divide the construction area into multiple grids, extract the elevation data in each grid, identify the elevation peak points of each grid and record the position coordinates, and obtain the peak point coordinate set; S102: Based on the peak point coordinate set, combine with the elevation values of adjacent positions, calculate the elevation offsets of the peak points in multiple directions, and obtain elevation offset data; S103: According to the elevation offset data, use the elevation offsets to detect elevation anomaly positions in the area and identify obstacles, including tower cranes and power towers. By extracting the size and position information of the target obstacles, generate obstacle identification information.

[0008] As a further solution of the present invention, the steps of calling the obstacle identification information, extracting the spatial position and geometric shape of the obstacle, and calculating the avoidance cost score by analyzing the increased path length, flight height adjustment amount, and path turning angle during avoidance flight in each UAV path to obtain the UAV measurement path are specifically as follows: S201: Call the obstacle identification information, extract the spatial position coordinates, vertical dimensions, and boundary contour points of each obstacle, and identify the positions with spatial overlap in each UAV path according to the spatial envelope range of each obstacle to obtain a path intersection segment data set; S202: According to the path intersection segment data set, analyze the path length, flight height adjustment amount, and path turning angle required for avoidance flight at multiple overlapping positions, calculate the avoidance cost score of each path, and obtain a path avoidance cost score set; S203: According to the path avoidance cost score set, calculate the priorities of multiple UAV paths according to the avoidance cost scores, and combine with the predicted flight duration to obtain the UAV measurement path.

[0009] As a further solution of the present invention, the specific formula for calculating the avoidance cost score of each path is: ; Calculate the comprehensive avoidance cost score value of the path to obtain a path avoidance cost score set; Wherein, is the comprehensive avoidance cost score value of the path numbered , is the unique identification number of the path, is the index number of the intersection segment in the path, is the total number of intersection segments in the path , is the avoidance path length value of the th intersection segment in the path , is the flight height adjustment amount required for the th intersection segment in the path , is the actual turning angle value on the th intersection segment in the path . is the average value of the turning angles corresponding to all paths at the th intersection segment, is the weight coefficient of the path length factor, is the weight coefficient of the flight altitude adjustment factor, is the weight coefficient of the turning angle deviation factor.

[0010] As a further solution of the present invention, the steps of obtaining the measurement path of the UAV, dynamically adjusting the flight altitude control parameters of the UAV according to the acquisition accuracy requirements and combining the real-time image clarity, and collecting the image information and lidar point cloud data of the construction area in real time to generate a regional imaging data set are specifically as follows: S301: Obtain the measurement path of the UAV, calculate the focal length parameters required by the image acquisition device according to the regional image acquisition accuracy requirements, and generate focal length-flight altitude comparison data by combining the flight altitude reference values required by the corresponding focal lengths; S302: According to the focal length-flight altitude comparison data, evaluate the clarity of the image by extracting the contrast value, edge sharpening degree value and target edge recognition rate of the image in real time during the measurement process, and obtain the image imaging clarity; S303: Dynamically update the flight altitude parameters according to the image imaging clarity, record the image acquisition time point and image sequence number, and obtain a regional imaging data set by combining the lidar point cloud data.

[0011] As a further solution of the present invention, the specific formula for evaluating the clarity of the image is: ; Calculate the image imaging clarity; where, represents the image imaging clarity, represents the maximum gray value in the image, represents the minimum gray value in the image, is a constant to prevent division by zero, is the number of edge pixel points, is the th gradient amplitude of the edge pixel, is the average gradient amplitude of all edge pixels, is the number of regions in the image, represents the th target boundary recognition rate of the region, is the average target boundary recognition rate of all regions, represents the index of the edge pixels in the image, represents the index of the regions in the image.

[0012] As a further solution of the present invention, according to the regional imaging dataset, using the flight altitude value, parsing the image and lidar point cloud data, reconstructing the construction engineering model of the construction area, and comparing it with the building information model of the target project, the steps of calculating the completion degree of the project in real time and identifying the project progress to generate the construction progress information are specifically as follows: S401: According to the regional imaging dataset, extract the flight altitude value and time stamp corresponding to each frame in the image information, and combine the spatial coordinate information of each measurement point in the lidar point cloud data to perform spatial coordinate registration on the image data and the point cloud data to generate a spatial data alignment result; S402: Based on the spatial data alignment result, reconstruct the building model in the target construction engineering area, and extract the geometric dimension information of the model to obtain the construction engineering model; S403: According to the construction engineering model, call the building information model of the target project and compare it with the actual construction model, calculate the completion degree of the project in multiple grids, identify the project progress and calculate the progress offset difference to obtain the construction progress information.

[0013] As a further solution of the present invention, the method further includes: S5: Call the construction progress information, extract the geometric center point coordinates, spatial extension direction and boundary dimensions of multiple components in the construction engineering model, combine the design position parameters and direction vectors of the corresponding components in the building information model, and calculate the spatial position offset degree of each component by comparing the spatial position differences between geometric center points, the deviation angles between direction vectors, and the coordinate distribution differences between boundary dimensions to obtain the construction engineering measurement record; The construction engineering measurement record includes the component actual spatial position set, the component offset angle matrix, and the component size deviation value distribution map.

[0014] As a further solution of the present invention, the steps of calling the construction progress information, extracting the geometric center point coordinates, spatial extension direction and boundary dimensions of multiple components in the construction engineering model, combining the design position parameters and direction vectors of the corresponding components in the building information model, and calculating the spatial position offset degree of each component by comparing the spatial position differences between geometric center points, the deviation angles between direction vectors, and the coordinate distribution differences between boundary dimensions to obtain the construction engineering measurement record are specifically as follows: S501: Call the construction progress information, extract the three-dimensional boundary data of the components in the construction engineering model in each grid, and generate a component spatial attribute dataset by extracting the geometric center point coordinates, spatial extension direction vector and boundary dimension value of each component; S502: Extract the design coordinates, design direction vectors, and standard boundary dimension data of the corresponding components in the building information model according to the component spatial attribute dataset, and compare them to obtain component spatial difference data; S503: Analyze the spatial position offset degree of each component according to the component spatial difference data, calculate the offset level, and map it to the construction engineering model to obtain the construction engineering measurement record.

[0015] The beneficial effects brought by the technical solutions provided in the embodiments of the present invention at least include: Through regional obstacle recognition, combined with the calculation of the avoidance cost scores of multiple measurement paths, the safety and rationality of the flight path planning are optimized. By using the linkage adjustment method of flight altitude and image clarity, the dynamic control of the imaging quality is realized. For the components of the construction engineering model, the construction progress is recognized and the dimensional deviation is detected, effectively improving the accuracy of the construction progress judgment and the spatial monitoring ability of the digital construction, providing a data basis for project management decisions and improving the construction quality. BRIEF DESCRIPTION OF THE DRAWINGS

[0016] In order to more clearly illustrate the technical solutions in the embodiments of the present invention, the following will briefly introduce the drawings required for the description of the embodiments. Obviously, the drawings in the following description are only some embodiments of the present invention. For those of ordinary skill in the art, without creative efforts, other drawings can be obtained based on these drawings.

[0017] Figure 1 It is a schematic diagram of the working process of the present invention. DETAILED DESCRIPTION OF THE EMBODIMENTS

[0018] The following will describe the technical solutions in the present invention with reference to the drawings.

[0019] In the embodiments of the present invention, words such as "exemplarily" and "for example" are used to represent examples, illustrations, or explanations. Any embodiment or design solution described as "exemplary" in the present invention should not be construed as being more preferred or having more advantages than other embodiments or design solutions. Exactly speaking, the use of the word "example" is intended to present concepts in a specific way. In addition, in the embodiments of the present invention, the meaning expressed by "and / or" can be both, or either one of the two.

[0020] In the embodiments of the present invention, "image" and "picture" can sometimes be used interchangeably. It should be noted that when the difference is not emphasized, the meanings they express are the same. "(of)", "corresponding", and "corresponding" can sometimes be used interchangeably. It should be noted that when the difference is not emphasized, the meanings they express are the same.

[0021] In the embodiments of the present invention, sometimes subscripts such as W 1 may be written in a non-subscript form such as W1. When the difference is not emphasized, the meanings they express are the same.

[0022] To make the technical problems, technical solutions, and advantages to be solved by the present invention clearer, the following will be described in detail with reference to the accompanying drawings and specific embodiments.

[0023] Please refer to Figure 1 , the present invention provides a technical solution, a digital construction method based on UAV three-dimensional mapping, including the following steps: S1: Import the regional digital surface model, divide the construction area into multiple grids and identify the elevation peak points, calculate the elevation offset of the peak points by combining the elevation values of adjacent positions of the peak points, identify multiple obstacles in the area, and generate obstacle identification information; S2: Call the obstacle identification information, extract the spatial position and geometric shape of the obstacles, calculate the avoidance cost score by analyzing the increased path length, flight height adjustment amount, and path turning angle in each UAV path for avoidance flight, and obtain the UAV measurement path; S3: Obtain the UAV measurement path, dynamically adjust the flight height control parameters of the UAV according to the acquisition accuracy requirements and in combination with the real-time image clarity, and collect the image information and lidar point cloud data of the construction area in real time to generate a regional imaging data set; S4: According to the regional imaging data set, use the flight height value to analyze the image and lidar point cloud data, reconstruct the construction engineering model of the construction area, compare it with the building information model of the target project, calculate the completion degree of the project in real time and identify the project progress, and generate construction progress information; Call the construction progress information, extract the geometric center point coordinates, spatial extension direction, and boundary dimensions of multiple components in the construction engineering model, combine the design position parameters and direction vectors of the corresponding components in the building information model, and calculate the spatial position deviation degree of each component by comparing the spatial position differences between geometric center points, the deviation angles between direction vectors, and the coordinate distribution differences between boundary dimensions, and obtain the construction engineering measurement record.

[0024] The obstacle identification information specifically includes the obstacle spatial distribution layer, obstacle geometric attribute data, and obstacle unique identification code. The UAV measurement path includes the start and end coordinate sets of the path segments, the path heading angle sequence, and the path corresponding obstacle avoidance strategy parameter set. The regional imaging data set includes the image resolution label, the point cloud coverage density map, and the image positioning coordinate set. The construction progress information specifically refers to the block completion status, the time identification corresponding to the process node, and the construction completion ratio. The construction engineering measurement record includes the actual spatial position set of the components, the component offset angle matrix, and the component size deviation value distribution map.

[0025] Import the regional digital surface model, divide the construction area into multiple grids and identify the elevation peak points, combine the elevation values ​​of the adjacent positions of the peak points, calculate the elevation offset of the peak points, identify multiple obstacles in the area, and generate obstacle identification information. The specific steps are as follows: S101: importing a regional digital surface model, dividing the construction area into multiple grids, extracting elevation data in each grid, identifying the elevation peak point of each grid and recording the position coordinates, and obtaining a peak point coordinate set; After importing the regional digital surface model, the two-dimensional range of the construction area is divided into 100×100 standard grids according to the rule of equal side length. The side length of each grid unit is set to 5 meters. After the grid division is completed, it is numbered in sequence, and the elevation data extraction operation is performed on each numbered grid. The extraction method is to screen all the elevation point data falling within the range of the grid boundary in the digital surface model, and record the spatial three-dimensional coordinate value of each elevation point. Then, all the elevation points in each grid are sorted by elevation value, and the coordinate point corresponding to the maximum elevation value is selected as the elevation peak point of the current grid. The spatial coordinates of this point are recorded in the peak point coordinate set. The judgment rule for the maximum value in the sorting is the largest of all elevation values ​​in the current grid, and the difference between it and the second largest elevation value in the same grid exceeds 1.2 The threshold is set by the minimum identifiable height of artificial structures in the project site. To avoid incorrect identification of low-undulation disturbances, if the difference is less than 1.2 meters, the current grid does not record the peak point. Through simulation experiments, if the maximum elevation is 47.2 meters and the second largest value is 46.3 meters, the difference is 0.9 meters, and the grid is not recorded as a peak point. If the maximum value is 52.5 meters and the second highest value is 50.8 meters, the difference is 1.7 meters, which exceeds the threshold, then the point is a valid peak point. Continue to traverse all grids to complete the peak extraction operation and obtain a peak point coordinate set with continuous spatial distribution. If the total number of grids is 20,000, 1,275 peak points are finally obtained. The constructed peak point coordinate set consists of the X coordinate, Y coordinate and Z value of each point to form a three-dimensional coordinate array in the form of ,in Indicates the peak point number, is the longitude coordinate, is the latitude coordinate, The array is the elevation value. It is used as the input basis for the next step of elevation offset calculation to obtain the peak point coordinate set.

[0026] S102: Based on the peak point coordinate set and in combination with the elevation values ​​of adjacent positions, the elevation offsets of the peak point in multiple directions are calculated to obtain elevation offset data; Based on the extracted peak point coordinate set, the elevation offset calculation operation is performed for each peak point. In this operation, each peak point is taken as the center point, and the center elevation points of the adjacent grids in the eight directions around it are retrieved. The elevation values ​​of the adjacent grids in each direction are differenced with the elevation value of the central peak point. The difference calculation adopts the formula ,in is the elevation value of the current peak point, For the The center point elevation of the adjacent grid in the direction. , the directions are arranged according to the main azimuth east, southeast, south, southwest, west, northwest, north, and northeast. After completing the difference calculation of the eight directions, the elevation offset value set of the peak point is recorded. This set reflects the degree of fluctuation of the peak point in each direction and can be used for subsequent abnormal area identification and judgment of the protruding structure of the structure. If the difference in a certain direction is greater than 2.5 meters, it is recorded as a strong protrusion direction. This threshold is derived from the average vertical elevation characteristics of artificial structures. According to the on-site investigation, the average height of the tower crane main structure from the base to the rotating platform is 3.2 meters. Considering the allowable error of DSM accuracy, the lower limit of the tolerance is set to 2.5 meters. If the peak point elevation is 51.6 meters and the east neighboring grid elevation is 48.5 meters, then , this direction meets the strong protrusion standard. If the other directions are all lower than 2.5 meters, it is recorded that this point has the protrusion feature only in the east direction. Finally, the 8-direction offset values ​​of all peak points form a matrix. The rows of the matrix correspond to the peak point numbers, and the columns correspond to the direction numbers, which are recorded as ,in Indicates the peak point, Indicates the direction number. After completing all calculations, the elevation offset data is generated.

[0027] S103: According to the elevation offset data, using the elevation offset, detecting the elevation abnormality position in the area and identifying obstacles, including tower cranes and power towers, and generating obstacle identification information by extracting the size and position information of the target obstacle; According to the elevation offset data, sequentially traverse the direction offset matrices of all peak points, and extract the points where the elevation difference exceeds 2.5 meters in any three consecutive directions. Mark these points as candidate points for the elevation anomaly area. Then, extract the number of other peak points within the 5×5 grid range around the candidate points. If the number exceeds 5, it is determined as a concentrated height difference fluctuation area, and enter the obstacle determination process. Further analyze the distribution density of the elevation peak points and the local maximum elevation value in this area. If the local maximum elevation exceeds 3 times the standard deviation of the average elevation in all adjacent grids around, and the number of discontinuous fluctuation directions around is less than 3, then identify this area as a candidate area for stable structure obstacles. Perform the structure recognition operation within this area, extract the vertical span value and horizontal projection area of all abnormal points, compare the structure morphology of the area where the span value exceeds 10 meters and the area is between 4 square meters and 36 square meters, match the standard outline dimensions of target structures such as tower cranes and electric towers, evaluate the similar structure type according to the dimension coincidence rate and boundary fitting degree, construct a structure label record, and establish an obstacle number in combination with the spatial position coordinates, and output the obstacle type and the corresponding coordinate and dimension information to complete the recording of the spatial obstacle characteristics and generate the obstacle identification information.

[0028] Call the obstacle identification information, extract the spatial position and geometric shape of the obstacle, and calculate the avoidance cost score by analyzing the increased path length, flight height adjustment amount, and path turning angle of the avoidance flight in each UAV path. The specific steps to obtain the UAV measurement path are as follows: S201: Call the obstacle identification information, extract the spatial position coordinates, vertical dimension, and boundary contour points of each obstacle, identify the positions with spatial overlap in each UAV path according to the spatial envelope range of each obstacle, and obtain the path intersection segment data set; Call the obstacle identification information. First, extract the spatial position coordinate information of each obstacle, including the longitude and latitude coordinates of the center point and the corresponding elevation value. At the same time, extract the vertical dimension data of the obstacle in three-dimensional space, record this dimension as the continuous elevation difference of the obstacle from the base to the top, and perform distributed sampling on each obstacle boundary contour point, and form a closed boundary set of its projection boundary in the horizontal coordinate system. On this basis, construct a three-dimensional envelope body according to the spatial envelope range of each obstacle and calibrate its spatial occupancy volume. Obtain the planned UAV measurement path data set, extract the sequence of flight segment coordinate points of each path, linearly connect all consecutive points in the path segment and construct a set of flight path line segments. Perform a spatial overlap judgment operation between the obstacle envelope body and the path line segment set. The judgment criterion is that there is an intersection between the path line segment and the obstacle envelope surface in any direction, and the intersection is within the area more than 1 meter below the top of the obstacle, which is regarded as a spatial intersection segment. By traversing each obstacle and all path segments in turn and judging whether the overlap condition is satisfied, record the starting point and ending point numbers of the path segments that meet the conditions and their respective path numbers, and generate a data set with the structure of path number, starting point index, ending point index, and intersection coordinates. If there are 36 paths in the survey area and 12 obstacles, and 54 intersection segments are obtained after detection, form a path intersection segment data set and obtain the path intersection segment data set.

[0029] S202: According to the path intersection segment data set, analyze the path length, flight height adjustment amount, and path turning angle required for avoiding flight around multiple overlapping positions, calculate the avoidance cost score for each path, and obtain the path avoidance cost score set; The specific formula for calculating the avoidance cost score for each path is: ; Calculate the comprehensive score value of the path avoidance cost and obtain the path avoidance cost score set; Among them, is the comprehensive avoidance cost score value of the path with path number , is the unique identification number of the path, is the index number of the intersection segment in the path, is the path the total number of intersection segments in, is the path in the th is the avoidance path length value of the th intersection segment in the path, is the flight height adjustment amount required for the th intersection segment in the path, is the actual turning angle value on the th intersection segment in the path, is at the The average of the turning angles corresponding to all paths at a crossing segment is the weight coefficient of the path length factor is the weight coefficient of the flight altitude adjustment factor is the weight coefficient of the turning angle deviation factor

[0030] Formula: Detailed explanation of the formula and the derivation process of formula calculation: The formula is used to calculate the comprehensive avoidance cost score of the path with path number on all crossing segments, and is used as the basis for the evaluation and sorting of the optimal scheduling of the UAV path

[0031] Meaning and setting values of parameters: is the current path number, set to path 05; is the number of crossing segments detected in path 05, determined by the overlapping relationship between the path and the spatial envelope of the obstacle, and the result of spatial registration analysis is 3 segments is the crossing segment index number, from 1 to 3; is the avoidance path length of the th segment in path 05, set to 28.3 meters is 32.7 meters is 25.1 meters; is the flight altitude adjustment amount of path 05 at the th segment, set to 6.2 meters is 4.7 meters is 5.9 meters; is the turning angle of path 05 at the th segment, set to 34°, is 48°, is 27°; is the average of the turning angles of all paths at the crossing segment , set , , ; , , are the weight coefficients of length, flight altitude, and angle deviation, set to , , ; Substitute the parameters into the formula for calculation: Cost score of the first segment: ; Cost score of the second segment: ; Cost score of the third segment: ; Sum: ; The result 50.29 represents the avoidance cost score of path 05 among all intersection segments, considering the detour length, flight altitude adjustment amount, and angle deviation. This value is a key indicator used for sorting in the subsequent path optimization scheduling.

[0032] S203: According to the set of path avoidance cost scores, based on the avoidance cost scores, combined with the predicted flight duration, calculate the priorities of multiple UAV paths to obtain the UAV measurement paths; According to the set of path avoidance cost scores, extract the cost scoring results of all paths, and calculate the flight duration by combining the flight distance of each path with the set flight speed. The flight duration calculation method is to divide the total path length by the flight speed, with the unit unified as seconds. Set the UAV flight speed to 8 meters per second. If the path length is 620 meters, the flight duration is 77.5 seconds. Normalize the cost score and flight duration of each path. The normalization value uses the interval mapping formula: , where is the original value, are the minimum and maximum values among all paths respectively. Set the two normalized values of the avoidance cost and flight duration to and respectively. Finally, combine them with weights to obtain the path priority index. The weights are 0.6 and 0.4 respectively. If the of a certain path, , then the priority value is . Arrange all paths in ascending order according to the priority value, extract the numbers of the top several paths and record their waypoint index sequences, and combine them to form the final UAV measurement path to obtain the UAV measurement path.

[0033] To obtain the UAV measurement path, according to the acquisition accuracy requirements, combined with the real-time image clarity, dynamically adjust the UAV flight altitude control parameters, and collect the image information and lidar point cloud data of the construction area in real time. The steps to generate the regional imaging dataset are as follows: S301: Obtain the UAV measurement path. According to the accuracy requirements of regional image acquisition, calculate the focal length parameters required by the image acquisition device, and combine with the flight height reference value required for the corresponding focal length to generate focal length - flight height comparison data; Obtain the UAV measurement path. According to the accuracy requirements of regional image acquisition, first extract the flight segment coordinate sequence and its corresponding number of waypoints in each UAV path, calculate the horizontal distance between waypoints for each flight segment, and combine with the sensor size and pixel size of the imaging device, and the vertical distance between the flight platform and the ground , and reverse - calculate the focal length parameters according to the standard spatial resolution formula: , where is the ground resolution, with the unit of centimeters per pixel, is the lens focal length, with the unit of millimeters. Set the target resolution to 3 cm / pixel, the pixel size to 4.2 μm, and the initial flight height value to 90 meters for substitution and calculation. Then the focal length is millimeters. To ensure imaging clarity, take multiple groups of flight height values at different heights and substitute them into the above formula respectively to obtain the corresponding required focal length values, construct a height - focal length mapping table, and organize each group of flight heights and their corresponding focal length values recorded in the table into data pairs , and sort them in ascending order of flight height to support subsequent dynamic adjustment operations. For a typical operation area, there are 5 commonly used heights of 70 meters, 80 meters, 90 meters, 100 meters, and 110 meters respectively. After calculation, the required focal lengths are 1.05 mm, 1.2 mm, 1.35 mm, 1.5 mm, and 1.65 mm respectively, and the corresponding records are taken as a set of comparison items to obtain the focal length - flight height comparison data.

[0034] S302: According to the focal length - flight height comparison data, evaluate the clarity of the image by extracting the contrast value, edge sharpness value, and target edge recognition rate of the image in real - time during the measurement process to obtain the image imaging clarity; The specific formula for evaluating the clarity of the image is: ; Calculate the image imaging clarity; where, represents the image imaging clarity, represents the maximum gray value in the image, represents the minimum gray value in the image, is a constant to prevent division by zero, is the number of edge pixel points, is the th gradient amplitude of the edge pixel, is the average gradient magnitude of all edge pixels, is the number of regions in the image, represents the target boundary recognition rate of the th region, is the average target boundary recognition rate of all regions, represents the index of the edge pixels in the image,

[0035] Formula: ; Detailed explanation of the formula and the derivation process of the formula calculation: This formula is used to calculate the imaging clarity of the image , and its result is used to evaluate the quality of the image. The formula includes three parts, namely the gray contrast part, the edge sharpening part, and the edge recognizability part.

[0036] Meaning of parameters and set values: is the maximum gray value of the image, set to 200, which reflects the maximum value of the image brightness; is the minimum gray value of the image, set to 50, which reflects the minimum value of the image brightness; is a constant to prevent division by zero, set to 0.01, which is used to prevent division by zero problems in calculations; is the number of edge pixels, set to 5000, is the gradient magnitude of the th edge pixel, representing the intensity of the edge part in the image, is the average gradient magnitude of all edge pixels, set to 110, representing the average intensity of all edge pixels, set ; is the number of regions in the image, set to 10, indicating that the image is divided into 10 regions, is the target boundary recognition rate of the th region, representing the proportion of the boundary of this region that is correctly recognized, is the average boundary recognition rate of all regions, set to 0.80, representing the average boundary recognition rate of all regions, set .

[0037] Substitute the parameters into the formula for calculation: Substitute the above parameters into the formula to calculate the clarity of the image: ; ; Result It shows the comprehensive clarity of the image. The image has a high contrast, good edge sharpening effect, and a high target boundary recognition rate, indicating that the image quality is good and suitable for further analysis and modeling.

[0038] S303: Dynamically update the flight altitude parameter according to the image imaging clarity, record the image acquisition time point and the image sequence number, and combine with the lidar point cloud data to obtain the regional imaging dataset; According to the image imaging clarity data, dynamically compare the evaluation results of each image. Record the flight altitude corresponding to the images that fail the clarity evaluation into the regulation queue, and correct the altitude of the flight segments to which these images belong. The correction value is set according to the missing items of the clarity index, and the altitude adjustment ratio is set to 5% step of the adjacent altitude groups in the focal length comparison table. If the shooting altitude of the current image is 90 meters and the recognition rate index is lower than the set threshold, then the flight altitude is corrected to meters, and mark the corresponding flight time point and image number information. At the same time, record the point cloud data frame number of the lidar system at the corresponding time. Synchronize the image acquisition and the laser point cloud with the same time stamp as the main key, and combine and construct a regional information record including the image frame number, shooting time, flight altitude, and point cloud frame number. Summarize the synchronization sequences of the images and point clouds in all operation segments to generate a full-time series dataset of the survey area and obtain the regional imaging dataset.

[0039] According to the regional imaging dataset, use the flight altitude value to analyze the image and lidar point cloud data, reconstruct the construction engineering model of the construction area, and compare it with the building information model of the target project. The steps to calculate the completion degree of the project in real time and identify the project progress to generate the construction progress information are as follows: S401: According to the regional imaging dataset, extract the flight altitude value and time mark corresponding to each frame in the image information, and combine with the spatial coordinate information of each measurement point in the lidar point cloud data to perform spatial coordinate registration on the image data and the point cloud data to generate the spatial data alignment result; According to the regional imaging dataset, first extract the flight altitude data corresponding to each frame of the image, record the acquisition time mark information of each frame of the image, and convert the pixel coordinates of each frame of the image into ground coordinates. The conversion process is to call the flight altitude value and the focal length , and use the interior orientation formula , , where is the pixel coordinate value of a point in the image coordinate system. In an actual case, if the pixel coordinates of a point in the image are (1200, 1500), the focal length mm, and the flight altitude meters, then the corresponding ground coordinates are mm, mm, which is (80 meters, 100 meters). Extract the spatial coordinates of each measurement point from the lidar point cloud data, including the longitude, latitude, and elevation coordinates of each measurement point in the point cloud data, and record them separately as . Use the ground coordinates obtained by image coordinate transformation to perform spatial alignment operations with the lidar point cloud measurement points. Specifically, compare the point coordinates after image coordinate transformation with the point cloud data measurement point coordinates, and adopt the minimum distance matching principle, that is, taking the calculated Euclidean distance value as an index. If the distance is less than the set spatial registration threshold of 2 meters, it is considered that the point is successfully registered. Traverse all image frame pixel coordinates and point cloud positions in sequence for registration matching operations. In an actual example, the coordinates of a certain point cloud point are (81 meters, 99 meters, 90.5 meters), and the calculated distance is meters, which is less than the threshold of 2 meters and is recorded as a successfully registered point. After completing the registration of all images and point clouds, record the image serial number of each group of successfully registered points and the corresponding point cloud measurement point index to obtain the spatial data alignment result.

[0040] S402: Based on the spatial data alignment result, reconstruct the building model in the target construction project area, and extract the geometric dimension information of the model to obtain the construction project model; Based on the spatial data alignment result, unify the point cloud data point sets of successfully registered points into the same spatial coordinate system, call all point cloud data, calculate the surface coordinate set of the area where each building target is located, and reconstruct the external surface form of the target building through the three-dimensional point cloud data. By extracting the three-dimensional boundary points of each part of the building, use the spatial coordinate solution method to obtain the geometric dimensions of the building components, including length, width, and height. Specifically, perform a search and matching operation on the farthest point pairs along the three axes of the boundary point set of each component. The search process performs a difference operation between the maximum and minimum values of the X, Y, and Z coordinates of the component boundary points. For example, if the maximum value of the X coordinate of the component point cloud in a certain construction area is 115 meters and the minimum value is 100 meters, then the length dimension of the component in the X direction is 15 meters. Determine the dimensions of the component in the Y and Z directions in the same way. Assuming that the coordinate difference in the Y direction is 12 meters and the coordinate difference in the Z direction is 8 meters, then the three-dimensional geometric dimensions of the component are obtained as length 15 meters, width 12 meters, and height 8 meters. Measure the dimensions of all components in the area in sequence, summarize and record the spatial positions and geometric dimension parameters of all components in the area, and establish a complete regional building structure information set to obtain the construction project model.

[0041] S403: According to the construction project model, call the building information model of the target project and compare it with the actual construction model, calculate the project completion degree in multiple grids, identify the project progress, and calculate the progress offset difference to obtain the construction progress information; According to the construction engineering model, extract the spatial positions and corresponding geometric dimension data of each component in the model, call the target engineering building information model that has completed the design, compare the actual construction model with the design model at the component level according to the spatial coordinates. Set the grid cell size to 5 meters × 5 meters, and use the actual construction volume of the construction model within each grid cell and the target volume of the design model for comparison, and calculate the project completion degree as . For example, if the volume of the design model of a certain grid cell is 200 cubic meters and the measured volume of the construction model is 160 cubic meters, then the completion degree of this grid cell is . If the completion degree is less than the set completion degree benchmark value of 90%, it is determined as a unit with slow construction progress. Separate records and marks are made for all slow grid cells. By calculating the difference between the actual completion degree and the design completion degree within all grid cells, that is, calculating the progress offset difference as . The difference of the above grid cell is cubic meters, that is, the offset difference is negative, indicating that the construction progress lags behind by 40 cubic meters. Summarize the completion degrees and their offset differences of all grid cells to establish a complete construction area progress information set and obtain the construction progress information.

[0042] Call the construction progress information, extract the geometric center point coordinates, spatial extension directions, and boundary dimensions of multiple components in the construction engineering model, combine the design position parameters and direction vectors of the corresponding components in the building information model, and calculate the spatial position offset degree of each component by comparing the spatial position differences between geometric center points, the deviation angles between direction vectors, and the coordinate distribution differences between boundary dimensions. The specific steps for obtaining the construction engineering measurement record are as follows: S501: Call the construction progress information, extract the three-dimensional boundary data of the components in the construction engineering model within each grid, and generate a component spatial attribute data set by extracting the geometric center point coordinates, spatial extension direction vectors, and boundary dimension values of each component; Call the construction progress information, sequentially read the completed component units within each spatial grid, extract the point cloud boundary point sets formed by them in the three-dimensional model, extract the farthest point pairs along the X, Y, and Z directions respectively for each component boundary point set, and perform difference operations on their coordinates to determine the boundary dimension values of the component in the three-axis directions. Suppose the maximum value of a certain component in the X-axis direction is 102.4 meters and the minimum value is 98.6 meters, then the X-direction dimension is meters. At the same time, calculate the arithmetic mean of all boundary point coordinates from the boundary point set to obtain the geometric center point coordinates. For example, if the X coordinates of the boundary points of a certain component are 98.6, 99.1, 100.4, and 102.4 in sequence, then the center point X coordinate is In meters, further based on the connection relationships of each boundary point in space, the principal component analysis method is used to obtain the direction of the maximum direction projection vector as the spatial extension direction vector. Let the extension vector be (0.87, 0.49, 0), indicating that the component extends in the northeast direction along the XY plane. After calculating the geometric center point coordinates, spatial extension direction vectors, and boundary dimensions of all components, they are structurally summarized according to the component numbers to form a structure data set with the component numbers as the primary key, and a component spatial attribute data set is generated.

[0043] S502: According to the component spatial attribute data set, extract the design coordinates, design direction vectors, and standard boundary dimension data of the corresponding components in the building information model and compare them to obtain the component spatial difference data; According to the component spatial attribute data set, read the corresponding design component parameters in the building information model one by one, compare the design coordinates and actual center point coordinates of each component with the same number, and perform three-dimensional space difference calculation. The position difference calculation formula is: , where are the center point coordinates of the component in the actual model respectively, is the center point coordinate of the component in the design model. If the actual component center point is (100.1, 36.7, 4.8) and the design center point is (100.0, 36.6, 4.9), then the position difference is meters, and the judgment of spatial direction difference is carried out. If the included angle is greater than 15 degrees, it is regarded as a direction deviation. In addition, compare the length, width, and height differences of the boundary dimensions. If the single-axis dimension difference exceeds 10% of the design dimension, it is marked as a component with abnormal dimensions. After comparing the positions, directions, and boundary dimensions of all components, output the difference value records of each component in the three types of dimensions to obtain the component spatial difference data.

[0044] S503: According to the component spatial difference data, analyze the spatial position offset degree of each component, calculate the offset level and map it to the construction engineering model to obtain the construction engineering measurement record; According to the component spatial difference data, extract the difference value ranges of each component in three dimensions: position, direction, and boundary dimensions. Set the threshold for dividing the offset levels. The position offset level is set as follows: level 1 for 0 - 0.1 meters, level 2 for 0.1 - 0.3 meters, and level 3 for greater than 0.3 meters. The direction offset level is divided as level 1 for 0 - 5 degrees, level 2 for 5 - 15 degrees, and level 3 for greater than 15 degrees. The dimension offset level is set as level 1 for a relative error of 0 - 5%, level 2 for 5% - 10%, and level 3 for greater than 10%. Assign level labels according to the interval where the difference value is located. Make a combined judgment on the level values of each component for the three types of indicators, and use the maximum offset level as the overall offset level of the component. For example, if the position offset of a component is 0.28 meters (level 2), the direction offset is 17 degrees (level 3), and the dimension offset is 4% (level 1), then the offset level of this component is level 3. Correlate this offset level with the component number, and mark its corresponding spatial coordinate position and level information in the construction model to generate an offset level distribution layer for all components in the area, and obtain the construction project measurement records.

[0045] The above embodiments can be implemented in whole or in part by software, hardware (such as circuits), firmware, or any other combination. When implemented using software, the above embodiments can be implemented in whole or in part in the form of a computer program product. A computer program product includes one or more computer instructions or computer programs. When the computer instructions or computer programs are loaded or executed on a computer, the processes or functions according to the embodiments of the present invention are generated in whole or in part. The computer can be a general-purpose computer, a special-purpose computer, a computer network, or other programmable devices. The computer instructions can be stored in a computer-readable storage medium, or transmitted from one computer-readable storage medium to another. For example, the computer instructions can be transmitted from one website, computer, server, or data center to another website, computer, server, or data center by wired (such as infrared, wireless, microwave, etc.) means. The computer-readable storage medium can be any available medium that the computer can access, or a data storage device such as a server or data center that contains one or more collections of available media. The available media can be magnetic media (such as floppy disks, hard disks, magnetic tapes), optical media (such as DVDs), or semiconductor media. The semiconductor media can be a solid-state drive.

[0046] It should be understood that the term "and / or" in this article is merely a description of the association relationship between associated objects, indicating that three relationships can exist. For example, A and / or B can represent: A exists alone, A and B exist simultaneously, and B exists alone. Here, A and B can be singular or plural. In addition, the character " / " in this article generally represents an "or" relationship between the associated objects before and after, but it may also represent an "and / or" relationship, which can be specifically understood with reference to the context before and after.

[0047] In the present invention, "at least one" means one or more, and "a plurality" means two or more. "At least one of the following" or similar expressions refer to any combination of these items, including any combination of single item(s) or plural item(s). For example, at least one of a, b, or c can represent: a, b, c, a - b, a - c, b - c, or a - b - c, where a, b, and c can be single or plural.

[0048] It should be understood that in various embodiments of the present invention, the magnitude of the serial numbers of the above - mentioned processes does not mean the sequence of execution. The execution sequence of each process should be determined by its function and internal logic, and should not constitute any limitation to the implementation process of the embodiments of the present invention.

[0049] Those of ordinary skill in the art can realize that the units and algorithm steps of each example described in combination with the embodiments disclosed herein can be implemented by electronic hardware, or a combination of computer software and electronic hardware. Whether these functions are executed in a hardware or software manner depends on the specific application and design constraints of the technical solution. Professionals can use different methods for each specific application to implement the described functions, but such implementation should not be considered to exceed the scope of the present invention.

[0050] Those skilled in the art can clearly understand that for the convenience and brevity of description, the specific working processes of the devices, apparatuses, and units described above can refer to the corresponding processes in the foregoing method embodiments, and will not be elaborated herein.

[0051] In several embodiments provided by the present invention, it should be understood that the disclosed devices, apparatuses, and methods can be implemented in other ways. For example, the device embodiments described above are merely illustrative. For example, the division of the units is only a logical function division. In actual implementation, there can be other division methods. For example, multiple units or components can be combined or integrated into another device, or some features can be ignored or not executed. Another point is that the displayed or discussed couplings or direct couplings or communication connections to each other can be through some interfaces. The indirect couplings or communication connections of the devices or units can be in electrical, mechanical, or other forms.

[0052] The units described as separate components may or may not be physically separated, and the components displayed as units may or may not be physical units, that is, they can be located in one place, or can be distributed to multiple network units. Some or all of the units can be selected according to actual needs to achieve the purpose of the solution of this embodiment.

[0053] In addition, in each embodiment of the present invention, each functional unit may be integrated in a processing unit, may exist physically alone for each unit, or two or more units may be integrated in one unit.

[0054] If the above-mentioned functions are implemented in the form of software functional units and sold or used as independent products, they can be stored in a computer-readable storage medium. Based on such an understanding, the technical solution of the present invention, in essence, or the part that contributes to the prior art, or a part of the technical solution, can be embodied in the form of a software product. The computer software product is stored in a storage medium and includes several instructions for causing a computer device (which may be a personal computer, a server, or a network device, etc.) to execute all or part of the steps of the methods described in the various embodiments of the present invention. The aforementioned storage medium includes: various media such as USB flash drives, mobile hard disks, read-only memory (ROM), random access memory (RAM), magnetic disks, or optical discs that can store program codes.

[0055] The above is only the specific implementation manner of the present invention, but the protection scope of the present invention is not limited thereto. Any person skilled in the art within the technical scope disclosed by the present invention can easily think of changes or substitutions, which should all be covered by the protection scope of the present invention. Therefore, the protection scope of the present invention should be subject to the protection scope of the claims.

Claims

1. A digital construction method based on drone three-dimensional mapping, characterized in that: The method comprises: S1: Import the regional digital surface model, divide the construction area into multiple grids and identify the elevation peak points, combine the elevation values ​​of the adjacent positions of the peak points, calculate the elevation offset of the peak points, identify multiple obstacles in the area, and generate obstacle identification information; S2: calling the obstacle identification information, extracting the spatial position and geometric shape of the obstacle, calculating the avoidance cost score by analyzing the path length, flight altitude adjustment and path turning angle added by the avoidance flight in each UAV path, and obtaining the UAV measurement path; S3: Obtain the measurement path of the UAV, dynamically adjust the flight altitude control parameters of the UAV according to the acquisition accuracy requirements and the real-time image clarity, collect image information and lidar point cloud data of the construction area in real time, and generate a regional imaging data set; S4: Based on the regional imaging data set, the flight altitude value is used to parse the image and the LiDAR point cloud data, reconstruct the construction project model of the construction area, and compare it with the building information model of the target project, calculate the completion of the project in real time, identify the project progress, and generate construction progress information.

2. The digital construction method based on drone three-dimensional mapping according to claim 1 is characterized in that: The obstacle identification information specifically includes the obstacle spatial distribution layer, obstacle geometric attribute data, and obstacle unique identification code. The UAV measurement path includes the path segment start and end coordinate sets, the path heading angle sequence, and the path corresponding obstacle avoidance strategy parameter set. The regional imaging data set includes image resolution labels, point cloud coverage density maps, and image positioning coordinate sets. The construction progress information specifically refers to the block completion status, the corresponding time mark of the process node, and the construction completion ratio.

3. The digital construction method based on drone three-dimensional mapping according to claim 1 is characterized in that: Import the regional digital surface model, divide the construction area into multiple grids and identify the elevation peak points, combine the elevation values ​​of the adjacent positions of the peak points, calculate the elevation offset of the peak points, identify multiple obstacles in the area, and generate obstacle identification information. The specific steps are as follows: S101: importing a regional digital surface model, dividing the construction area into multiple grids, extracting elevation data in each grid, identifying the elevation peak point of each grid and recording the position coordinates, and obtaining a peak point coordinate set; S102: Based on the peak point coordinate set and in combination with the elevation values ​​of adjacent positions, the elevation offsets of the peak points in multiple directions are calculated to obtain elevation offset data; S103: According to the elevation offset data, using the elevation offset, detecting elevation abnormality locations in the area and identifying obstacles, including tower cranes and power towers, and generating obstacle identification information by extracting size and location information of target obstacles.

4. The digital construction method based on drone three-dimensional mapping according to claim 3 is characterized in that: The obstacle identification information is called, the spatial position and geometric shape of the obstacle are extracted, and the avoidance cost score is calculated by analyzing the path length, flight altitude adjustment and path turning angle added by the avoidance flight in each UAV path. The specific steps of obtaining the UAV measurement path are as follows: S201: calling the obstacle identification information, extracting the spatial position coordinates, vertical dimensions and boundary contour points of each obstacle, identifying the spatially overlapping positions in each UAV path according to the spatial envelope range of each obstacle, and obtaining a path intersection segment dataset; S202: Analyzing the path lengths, flight altitude adjustment amounts, and path turning angles required for avoidance flights at multiple overlapping positions according to the path intersection segment data set, calculating the avoidance cost score of each path, and obtaining a path avoidance cost score set; S203: According to the path avoidance cost score set, according to the avoidance cost score, and in combination with the predicted flight duration, the priorities of multiple UAV paths are calculated to obtain the UAV measurement path.

5. The digital construction method based on drone three-dimensional mapping according to claim 4 is characterized in that: The specific formula for calculating the avoidance cost score of each path is: ； Calculate the comprehensive score of the path avoidance cost and obtain the path avoidance cost score set; in, The path number is The comprehensive score of the path avoidance cost, is the unique identification number of the path. is the index number of the intersection segment in the path, For path The total number of crossover segments, For path Middle The length of the detour path of the intersection segment, For path Middle The required flight altitude adjustment for each intersection segment is: For path Middle The actual steering angle value on the intersection segment, For the The average value of the steering angles corresponding to all paths at the intersection segment, is the weight coefficient of the path length factor, is the weight coefficient of the flight altitude adjustment factor, is the weight coefficient of the steering angle deviation factor.

6. The digital construction method based on drone three-dimensional mapping according to claim 4 is characterized in that: The steps of obtaining the UAV measurement path, dynamically adjusting the flight altitude control parameters of the UAV according to the acquisition accuracy requirements and the real-time image clarity, and collecting the image information and lidar point cloud data of the construction area in real time to generate the regional imaging data set are as follows: S301: Obtain the measurement path of the drone, calculate the focal length parameters required by the image acquisition device according to the regional image acquisition accuracy requirements, and generate focal length flight height comparison data in combination with the flight height reference value required for the corresponding focal length; S302: According to the focal length flying height comparison data, by extracting the contrast value, edge sharpness value and target edge recognizability of the image in real time during the measurement process, evaluating the clarity of the image to obtain the image imaging clarity; S303: dynamically update the flight altitude parameters according to the image clarity, record the image acquisition time point and image sequence number, and obtain the regional imaging data set in combination with the laser radar point cloud data.

7. The digital construction method based on drone three-dimensional mapping according to claim 6 is characterized in that: The specific formula for evaluating the clarity of the image is: ； Calculate image clarity; in, Represents the image clarity, Represents the maximum gray value in the image. Represents the minimum gray value in the image. To prevent division by zero, is the number of edge pixels, For the The gradient magnitude of edge pixels, is the average gradient magnitude of all edge pixels, is the number of regions in the image, Representative The target boundary recognition rate of the region is is the average target boundary recognition rate of all regions, represents the index of edge pixels in the image, Represents the index of the region in the image.

8. The digital construction method based on drone three-dimensional mapping according to claim 6 is characterized in that: According to the regional imaging data set, the flight altitude value is used to analyze the image and the LiDAR point cloud data, reconstruct the construction project model of the construction area, and compare it with the building information model of the target project, calculate the completion degree of the project in real time and identify the project progress. The specific steps of generating the construction progress information are as follows: S401: extracting the flight altitude value and time mark corresponding to each frame in the image information according to the regional imaging data set, combining the spatial coordinate information of each measurement point in the laser radar point cloud data, performing spatial coordinate registration on the image data and the point cloud data, and generating a spatial data alignment result; S402: reconstructing a building model in the target construction project area based on the spatial data alignment result, extracting geometric dimension information of the model, and obtaining a construction project model; S403: According to the construction project model, the building information model of the target project is called and compared with the actual construction model, the project completion degree in multiple grids is calculated, the project progress is identified and the progress offset difference is calculated, and the construction progress information is obtained.

9. The digital construction method based on drone three-dimensional mapping according to claim 1 is characterized in that: The method further comprises: S5: calling the construction progress information, extracting the geometric center point coordinates, spatial extension directions and boundary dimensions of multiple components in the construction project model, combining the design position parameters and direction vectors of the corresponding components in the building information model, and calculating the spatial position deviation degree of each component by comparing the spatial position difference between the geometric center points, the deviation angle between the direction vectors, and the coordinate distribution difference between the boundary dimensions, so as to obtain the construction project measurement record; The construction engineering measurement record includes a component actual spatial position set, a component offset angle matrix, and a component size deviation value distribution diagram.

10. The digital construction method based on drone three-dimensional mapping according to claim 9 is characterized in that: The construction progress information is called to extract the geometric center coordinates, spatial extension directions and boundary dimensions of multiple components in the construction project model. Combined with the design position parameters and direction vectors of the corresponding components in the building information model, the spatial position difference between the geometric center points, the deviation angle between the direction vectors, and the coordinate distribution difference between the boundary dimensions are compared to calculate the spatial position deviation degree of each component. The steps of obtaining the construction project measurement record are specifically as follows: S501: calling the construction progress information, extracting the three-dimensional boundary data of the components in the construction engineering model in each grid, and generating a component space attribute data set by extracting the geometric center point coordinates, space extension direction vector and boundary size value of each component; S502: extracting design coordinates, design direction vectors and standard boundary dimension data of corresponding components in the building information model according to the component space attribute data set and performing comparison to obtain component space difference data; S503: Analyze the spatial position offset degree of each component according to the component spatial difference data, calculate the offset grade and map it to the construction project model, and obtain the construction project measurement record.

Images