Installation detection method for building special-shaped multi-curved outer envelope structure of unmanned aerial vehicle three-dimensional reconstruction

Abstract

The application discloses an installation detection method for a building special-shaped multi-curved outer envelope structure of a three-dimensional reconstruction of a UAV, belongs to the technical field of construction methods, and performs three-dimensional reconstruction on a target building image sequence photographed by the UAV, performs precision rechecking on corresponding positions of a structure installation position of a risk-sensitive area of a real scene three-dimensional model and a BIM design model, generates deviation data between the real scene three-dimensional model and the BIM design model through registration deviation analysis, marks the positions with deviations on the real scene three-dimensional model, can identify the structure installation positions affected by galloping vibration and vortex-induced vibration, and each outer envelope structure affecting long-term stability, performs precision rechecking on the real scene three-dimensional model and the BIM design model, and guides construction personnel to rectify the quality of an installation project through registration deviation analysis data, overcomes the limitations of traditional technologies, significantly improves construction quality management efficiency and precision, and thus reduces construction risks.

Description

Technical Field

[0001] This disclosure belongs to the fields of computer-aided design technology and construction management technology, specifically involving a method for the installation and inspection of irregularly shaped multi-curved building envelope structures reconstructed by UAV 3D reconstruction. Background Technology

[0002] In the contemporary architecture industry, buildings with irregular and curved facades, with their unique artistic expression and spatial creativity, are gradually becoming an important choice for urban landmark buildings. These buildings break away from the regular forms of traditional architecture, characterized by complex shapes such as curved surfaces, broken lines, and irregular blocks. Their key components, such as roofs and facades, often exhibit non-linear spatial variations, placing extremely high demands on the accuracy of construction surveying and the adaptability of detailed design. When hyperbolic facade envelope structures and roof metal panel systems are applied to high-rise buildings, large public buildings (such as airports and stadiums), and special buildings (large-scale, functionally complex, and densely populated buildings), the building will experience vertical compressive deformation between the structural layers and the outer envelope under the influence of external risk factors such as gravity loads or wind loads. To guide airflow and reduce wind pressure and vibration on the external envelope (or improve the roof's resistance to overturning, etc.), hyperboloid facade designs are often based on aerodynamic optimization by architects. Insufficient layout accuracy can lead to deviations between the actual installed curved surface and the design model (unevenness and discontinuous curvature in the actual formed area). These deviations disrupt the streamlined design, generating turbulence and other adverse factors, increasing local wind loads, and even inducing harmful vibrations. Furthermore, curved envelope systems are prone to accumulating errors during installation, and errors in the frame installation will further affect the subsequent processing and production of curved panels. Traditional measurement methods, limited by the complex spatial angles and irregular shapes of irregular structures, rely on ground-based measurement methods such as total stations or 3D scanning. These methods often suffer from blind spots, making it difficult to cover all key points, resulting in low measurement efficiency and an inability to accurately capture the continuous changing trends of the structure, leading to significant discrepancies between the acquired data and the actual structure. Summary of the Invention

[0003] The purpose of this invention is to propose a method for the installation and inspection of irregularly shaped multi-curved building envelope structures using UAV 3D reconstruction, in order to solve one or more computer-aided design problems existing in the prior art, and at least provide a beneficial option or create conditions.

[0004] To achieve the above objectives, according to one aspect of this disclosure, a method for detecting the installation of irregularly shaped, multi-curved building envelope structures reconstructed by unmanned aerial vehicles (UAVs) is provided, the method comprising the following steps:

[0005] S100: Import the BIM design model of the target building and mark the location of the external envelope structure on the BIM design model as the structural installation location.

[0006] S200 uses computer-aided design simulation software to simulate and obtain the surface wind pressure distribution of the BIM design model;

[0007] S300, under different wind speed conditions, loads the surface wind pressure distribution onto the BIM design model to obtain the vibration displacement of each structural installation location under different average wind speeds, and identifies and marks the structural installation locations with risk-sensitive areas on the BIM design model.

[0008] The S400 uses drones to capture a sequence of images of a target building, and then performs 3D reconstruction on the image sequence to obtain a real-world 3D model of the target building.

[0009] S500 performs accuracy verification on the corresponding positions of the structural installation locations in risk-sensitive areas of the real-world 3D model and the BIM design model during the main structure installation stage, the enclosure structure frame installation stage, and the enclosure structure panel installation stage. It also generates deviation data between the real-world model and the BIM design model through registration deviation analysis.

[0010] The S600 marks the locations with deviations on the real-world 3D model, outputs and displays the real-world 3D model.

[0011] Furthermore, in S100, the BIM design model of the target building is a static BIM 3D model of the target building constructed using BIM technology based on the design drawings of the target building.

[0012] The external envelope includes a curtain wall system and a metal roofing system. The curtain wall system includes a metal curtain wall or glass curtain wall, a supporting structure, and a connection and sealing system. The supporting structure includes connectors or purlins or purlin brackets. The connection and sealing system includes silicone sealant, sealing strips or hinges, and bolts.

[0013] Furthermore, in S200, the method for simulating and obtaining the surface wind pressure distribution of the BIM design model using computer-aided design simulation software is as follows: import the BIM design model into the computer-aided design simulation software, and perform flow field analysis under given wind speed boundary conditions to obtain the surface wind pressure distribution of the BIM design model. The specific steps of the flow field analysis include: importing the BIM design model, surface segmentation, mesh generation, and mesh attribute definition; setting material parameters, geometric properties, wind speed boundary conditions, load conditions, and defining the connectivity of elements.

[0014] Preferably, the flow field analysis specifically involves: simulating transient wind pressure distribution under different wind angles (0°~360°, interval 22.5°) in the surrounding flow field domain of the BIM design model of the target building; coupling the wind pressure data to the structural finite element model; applying dynamic wind pressure loads; performing time history analysis and identifying resonance areas; and exporting the time history data (.csv / .dat format) of the surface wind pressure coefficient of the BIM design model, as well as exporting the wind pressure distribution data and displacement cloud map of the BIM design model.

[0015] Furthermore, the BIM design model is in RVT or 3DM format.

[0016] Furthermore, S100 also includes segmenting and extracting the 3D mesh of the wall surface of the BIM design model, calculating the 3D support plane fitted by the least squares method, calculating the Euclidean distance between the 3D mesh and the 3D support plane, and normalizing it to generate a depth map of the wall surface; generating a texture image of the wall surface through 3D spatial texture mapping based on the calibrated camera parameters; and ensuring that the length and width dimensions of the depth map and texture map of the wall surface are consistent.

[0017] Preferably, the method for mesh generation is the sweep method.

[0018] The given wind speed boundary conditions range from 5 m / s to 36.9 m / s.

[0019] Preferably, different wind speeds are set at different heights in the BIM design model of the target building for simulation. The given wind speed boundary conditions are simulated as follows: the wind speed in the area corresponding to the location of the target building's BIM design model at 0-10 meters is set to 5 m / s, the wind speed at 10-30 meters is set to 8.7 m / s, the wind speed at 30-90 meters is set to 15 m / s, and the wind speed above 90 meters is set to 36.9 m / s.

[0020] Note: The "Code for Design of Building Structures" GB 50009-2012 stipulates that high-rise buildings require a wind speed of level 12 (32.7~36.9m / s) as the design benchmark; while the wind speed at a height of 10 meters is 5m / s, at 30 meters it reaches 8.7m / s, and at 90 meters it rises to 15m / s. See reference: Fang Xilin, Sun Pengou. Disaster Prevention Reader for Building Engineering [M]. China Environmental Science Press, 2014.

[0021] Preferably, the wind speed boundary conditions are set based on local meteorological data (probability distribution of wind speed and wind direction, refer to the "Code for Design of Building Structures" (GB 50009-2012)).

[0022] Preferably, the computer-aided design simulation software is the open-source OpenFOAM software or Elmer FEM software.

[0023] Furthermore, in S300, different wind speed conditions are wind speeds ranging from 5 m / s to 36.9 m / s.

[0024] Furthermore, in S300, the method for identifying and marking the structural installation locations of risk-sensitive areas on the BIM design model is as follows:

[0025] Mark the wind pressure orientation area in the adjacent structural installation locations of the current structural installation location;

[0026] Based on linear wind pressure conditions, structural installation locations in risk-sensitive areas are selected within the wind pressure orientation region.

[0027] Furthermore, the wind pressure orientation area is marked in the adjacent structural installation locations of the current structural installation location, specifically as follows:

[0028] Based on the surface wind pressure distribution, obtain the wind pressure values ​​corresponding to all points at each structural installation location in the BIM design model. Record the vibration displacement of the geometric center of gravity at each structural installation location from the minimum wind pressure value to the maximum wind pressure value as the wind pressure limit displacement. (Due to differences in material parameters, geometric properties, wind speed boundary conditions, and load conditions, the wind pressure limit displacement will vary at different locations and heights in the BIM design model when subjected to wind pressure.)

[0029] Other structural installation positions adjacent to the current structural installation position are recorded as the adjacent structural installation positions of the current structural installation position (regardless of whether the external envelope is a component-type building curtain wall, unitized curtain wall, glass curtain wall, stone curtain wall, metal panel curtain wall, or all-glass curtain wall, they are all interconnected between individual structural installation positions. Therefore, each structural installation position is surrounded by other curtain walls. The meaning of adjacent is: other structural installation positions that are the shortest distance from the current structural installation position).

[0030] Let MaxCUP be the geometric centroid of the structural installation position with the largest wind pressure limit displacement among all adjacent structural installation positions of the current structural installation position; let MinCUP be the geometric centroid of the structural installation position with the smallest wind pressure limit displacement among all adjacent structural installation positions of the current structural installation position; then, a ray CUPL is formed from MinCUP to MaxCUP; all structural installation positions that CUPL passes through in sequence are recorded as the wind pressure orientation region of the current structural installation position; (all wind pressure orientation regions that CUPL passes through in sequence are areas with the largest limit displacement and continuous wind vibration or large wind load risk at the current structural installation position. Since the unit panels of the glass curtain wall are usually installed in an interconnected manner to ensure the airtightness, watertightness, and wind pressure resistance of the overall structure, the connection method, force transmission, and deformation adaptation between the various glass curtain wall units are usually affected by the adjacent structural installation positions. Due to the continuous wind vibration, the dynamic deformation of the structural displacement under wind vibration will cause stress to gradually diffuse and transfer from the area with large wind pressure limit displacement to the area with low wind pressure limit displacement. The range of diffusion and transfer is the wind pressure orientation region marked in the above steps).

[0031] Furthermore, based on linear wind pressure conditions, structural installation locations in risk-sensitive areas are selected within the wind pressure orientation region, specifically as follows:

[0032] Let WV(j) be the j-th wind pressure orientation region at the current structural installation location; j is the sequence number, and the order of the wind pressure orientation regions is the order in which CUPL passes through the wind pressure orientation regions.

[0033] Determine each WV(j). If WV(j) satisfies the linear wind pressure condition, then mark WV(j) as the structural installation location of the risk-sensitive area.

[0034] The linear wind pressure condition is: MV(WV(j-1))>MV(WV(j)) and MV(WV(j))<MV(WV(j+1)), where MV(WV(j)) is the wind pressure limit displacement corresponding to WV(j);

[0035] The above methods quickly locate risk-sensitive structural installation locations by determining the wind pressure limit displacement of discontinuous installation locations. These risk-sensitive areas cannot rely on adjacent structural installation locations for stress transfer and are therefore likely to be the areas most affected by continuous wind vibration. However, if the structural installation location can only be detected as a discrete point, the above methods have low accuracy when identifying a continuous area affected by wind vibration. Therefore, identification is required using the following linear wind pressure conditions:

[0036] Preferably, the linear wind pressure condition is:

[0037] DE(WV(j-1))>DE(WV(j)) and DE(WV(j))<DE(WV(j+1)),

[0038] Where DE(WV(j)) is the inertial wind pressure index of the j-th wind pressure orientation region WV(j) at the current structural installation location;

[0039] Preferably, the inertial wind pressure index DE(WV(j)) is calculated as follows:

[0040] DE(WV(j))=|MinA1-MinA2|÷|MaxA1-MaxA2|;

[0041] Where MaxA1 is the maximum value of the wind pressure limit displacement corresponding to the wind pressure orientation region from the 1st to the (j-1th)th region; MinA1 is the minimum value of the wind pressure limit displacement corresponding to the wind pressure orientation region from the 1st to the (j-1th)th region; MaxA2 is the maximum value of the wind pressure limit displacement corresponding to the wind pressure orientation region from the (j+1th)th to the SuMth region; MinA2 is the minimum value of the wind pressure limit displacement corresponding to the wind pressure orientation region from the (j+1th)th to the SuMth region; and SuM is the number of wind pressure orientation regions.

[0042] The inertial wind pressure index DE(WV(j)) calculated by the above method is used to determine whether the wind pressure influence has continuous inertia and whether it affects the vertical compression deformation, expansion deformation, and stress concentration of the surrounding adjacent outer envelope due to wind vibration in the form of the ratio of the convergent wind pressure difference to the transmitted wind pressure difference. However, due to the characteristics of wind pressure, when the wind bypasses the edge of the building, the Karman vortex street effect will occur periodically, resulting in galloping and vortex-induced resonance caused by alternating wind pressure. The outer envelope continuously absorbs energy in the wind due to aerodynamic negative damping, resulting in an increase in amplitude. The above method cannot identify the risk-sensitive areas of the structural installation location affected by galloping and vortex-induced resonance, but can only identify the structural installation location affected by continuous wind vibration. In order to accurately identify the structural installation location affected by galloping and vortex-induced resonance, it is necessary to identify it through the following linear wind pressure conditions, specifically:

[0043] Preferably, the inertial wind pressure index DE(WV(j)) is calculated as follows:

[0044] The specific method for calculating the convergent and transmitted air pressure differences at all structural installation locations is as follows:

[0045] The convergent wind pressure difference is defined as the difference between the average value of the wind pressure limit displacement corresponding to the wind pressure orientation region from the 1st to the (j-1th)th structural installation position and the corresponding wind pressure limit displacement at the structural installation position; the transmitted wind pressure difference is defined as the difference between the average value of the wind pressure limit displacement corresponding to the wind pressure orientation region from the (j+1th)th to the Nth structural installation position and the corresponding wind pressure limit displacement at the structural installation position. (The transmitted wind pressure difference represents the degree to which the wind pressure limit displacement increases along the wind pressure orientation region starting from the structural installation position, while the convergent wind pressure difference represents the degree to which the wind pressure limit displacement decreases along the wind pressure orientation region starting from the structural installation position.)

[0046] Traverse the wind pressure orientation regions corresponding to each structural installation location. Let MaxS represent the number of wind pressure orientation regions where the convergent wind pressure difference is greater than the transmitted wind pressure difference; let MinS represent the number of wind pressure orientation regions where the convergent wind pressure difference is less than the transmitted wind pressure difference; then the inertial wind pressure index DE(WV(j)) is calculated as follows: Inertial wind pressure index DE(WV(j)) = MaxS ÷ (MaxS + MinS).

[0047] The inertial wind pressure index DE(WV(j)) calculated by the above method is determined by the ratio of the convergent wind pressure difference to the transmitted wind pressure difference at the structural installation location. It judges the degree of inertial diffusion of wind pressure limit displacement caused by galloping and vortex-induced resonance due to alternating wind pressure generated by the Karman vortex street effect in certain continuous areas. This accurately quantifies and marks the structural installation location with a high wind pressure limit displacement, which will be affected by wind pressure inertia and thus fail to effectively adapt to dynamic deformation under wind vibration. It identifies the structural installation location affected by galloping and vortex-induced resonance, and various external envelope structures that affect long-term stability.

[0048] Furthermore, in the S400, the method for capturing image sequences of target buildings using drones is as follows: using drones as carriers, generating acquisition flight paths, and the image sequences of drones flying along the flight paths ultimately form image sequences of the target buildings.

[0049] Furthermore, in S400, the real-world 3D model of the target building is obtained by performing 3D reconstruction on the image sequence: the real-world 3D model is obtained by performing 3D reconstruction on the image sequence using 3D Reconstruction or MeshRoomCL 3D reconstruction tools.

[0050] Furthermore, in S500, during the main structure installation stage, the enclosure structure frame installation stage, and the enclosure structure panel installation stage, the accuracy of the corresponding positions of the structural installation locations in risk-sensitive areas between the real-world 3D model and the BIM design model is verified. The method for generating deviation data between the real-world model and the BIM design model through registration deviation analysis is as follows:

[0051] The ICP algorithm is used to verify the accuracy of the corresponding positions of the structural installation locations in risk-sensitive areas between the real-scene 3D model and the BIM design model. The deviation data between the real-scene model and the BIM design model is generated through registration deviation analysis.

[0052] The ICP algorithm-based fusion analysis involves registering the real-world 3D model with the BIM design model and calculating the deviation data for the corresponding structural installation locations in risk-sensitive areas. The deviation data is displayed using color coding to show the dimensional deviation ranges for different areas.

[0053] The deviation data includes overall deviation statistics (mean, standard deviation, maximum deviation), deviation heatmap (color mapping distance values, visually displaying error distribution), and structural installation location error analysis data for structures in risk-sensitive areas: deviations are grouped and statistically analyzed according to structures (external envelope structures) in risk-sensitive areas.

[0054] The ICP algorithm is the ICP registration algorithm in the open-source 3D data processing library Open3D.

[0055] During the main structure installation phase, the corresponding positions of the structural installation locations in the real-world 3D model and the BIM design model include: the steel structure / concrete main body (the load-bearing core of irregular-shaped buildings), connectors or purlins or purlin supports, steel columns, steel beams or space trusses (which need to be pre-bent into the design curve), and concrete shear walls or irregular-shaped core tubes (which achieve curved shapes through template customization) in the real-world 3D model and the BIM design model.

[0056] Among them, the structural installation positions in the above main structure installation stage provide stable support for the irregular-shaped outer enclosure structure, ensuring the positioning benchmark for the subsequent frame and panels.

[0057] During the enclosure structure frame installation phase, the corresponding positions of the structural installation locations in the real-world 3D model and the BIM design model include: main keel: usually an arc-shaped steel beam or aluminum profile, arranged along the normal direction of the curved surface; secondary keel: perpendicular to the main keel, forming a grid-like support (which needs to be pre-bent into a double curvature shape); and connection nodes: three-dimensional adjustable hinged supports (adapting to multi-directional degree of freedom adjustment).

[0058] In the main structure installation phase, the structural installation location is a secondary support system adapted to irregular curved surfaces, which provides an installation carrier for the panels.

[0059] The corresponding positions of the structural installation locations in the real-world 3D model and the BIM design model during the enclosure panel installation phase include: metal panels (aluminum panels, titanium-zinc panels): CNC stamping forming of hyperbolic panels; glass curtain walls: cold bending technology or triangular plane glass splicing fitting curved surfaces; composite material panels (GFRP, ETFE): factory prefabricated integral curved surface units; and the corresponding positions of the connection and sealing system, including silicone sealant, sealing strips or hinges, and bolts, in the real-world 3D model and the BIM design model.

[0060] Among them, the structural installation position in the enclosure panel installation stage is to complete the precise installation of irregularly shaped multi-curved panels to ensure the appearance and airtightness.

[0061] Furthermore, in the S500, the method for outputting and displaying the real-world 3D model is as follows: adaptive histogram equalization is performed on the locations of deviations in the real-world 3D model to enhance local contrast.

[0062] Preferably, it also includes transmitting the real-scene 3D model to a client or mobile device for display output.

[0063] Preferably, the method for transmitting a real-scene 3D model to a client or mobile device for display output includes the following steps: processing the real-scene 3D model using an edge-folding algorithm to obtain a simplified model;

[0064] The simplified model is divided into 1m×1m blocks and compressed separately.

[0065] Dynamically push visible blocks based on the view of the client or mobile device (using HTTP / 2 ServerPush).

[0066] When loading on the client or mobile device, check GPU performance and select the corresponding LOD level:

[0067] High-end equipment: LOD0 (original precision);

[0068] Mid-range equipment: LOD1 (50% surface area);

[0069] Low-end equipment: LOD2 (20% surface area + solid color heatmap);

[0070] The installation locations of structures in risk-sensitive areas will be highlighted in real time.

[0071] Among them, high-end equipment: supports hardware-accelerated ray tracing, floating-point computing power ≥3 TFLOPS, and memory bandwidth ≥100GB / s;

[0072] Mid-range devices: Support Vulkan 1.1 and above API, with floating-point computing power of 1~3 TFLOPS;

[0073] Low-end devices: Only support OpenGL ES 3.0, floating-point computing power <1 TFLOPS.

[0074] The beneficial effects of this invention are as follows: This invention provides a method for detecting the installation of irregularly shaped, multi-curved building envelope structures using UAV 3D reconstruction. It accurately quantifies and marks the structural installation locations with high wind pressure limit displacement that are affected by wind pressure inertia and thus fail to effectively adapt to dynamic deformation under wind vibration. It can identify the structural installation locations affected by galloping and vortex-induced resonance, as well as various envelope structures that affect long-term stability. It performs accuracy verification between the real-world 3D model and the BIM design model, accurately quantifies the differences between construction and design, and guides construction personnel to rectify the quality of the installation project through registration deviation analysis data. It overcomes the limitations of traditional technologies, significantly improves the efficiency and accuracy of construction quality management, and thus reduces construction quality risks. Attached Figure Description

[0075] Figure 1 The diagram shows a flowchart of the installation and inspection method for irregularly shaped, multi-curved building envelope structures reconstructed by UAV in 3D. Detailed Implementation

[0076] The following will provide a clear and complete description of the concept, specific structure, and resulting technical effects of this disclosure in conjunction with the embodiments and accompanying drawings, so as to fully understand the purpose, solution, and effects of this disclosure. It should be noted that, unless otherwise specified, the embodiments and features described in this application can be combined with each other.

[0077] Example 1

[0078] Figure 1 The diagram shows a flowchart of the installation and inspection method for irregularly shaped, multi-curved building envelope structures reconstructed by UAV 3D reconstruction. The following section will combine... Figure 1 This invention describes a method for detecting the installation of irregularly shaped, multi-curved building envelope structures using UAV 3D reconstruction according to an embodiment of the present invention. The method includes the following steps:

[0079] S100: Import the BIM design model of the target building and mark the location of the external envelope structure on the BIM design model as the structural installation location.

[0080] S200 uses computer-aided design simulation software to simulate and obtain the surface wind pressure distribution of the BIM design model;

[0081] S300, under different wind speed conditions, loads the surface wind pressure distribution onto the BIM design model to obtain the vibration displacement of each structural installation location under different average wind speeds, and identifies and marks the structural installation locations with risk-sensitive areas on the BIM design model.

[0082] The S400 uses drones to capture a sequence of images of a target building, and then performs 3D reconstruction on the image sequence to obtain a real-world 3D model of the target building.

[0083] S500 performs accuracy verification of the corresponding positions of structural installation locations in risk-sensitive areas between the real-world 3D model and the BIM design model, and generates deviation data between the real-world model and the BIM design model through registration deviation analysis.

[0084] The S600 marks the locations with deviations on the real-world 3D model, outputs and displays the real-world 3D model.

[0085] Furthermore, in S100, the BIM design model of the target building is a static BIM 3D model of the target building constructed using BIM technology based on the design drawings of the target building.

[0086] The external envelope includes a curtain wall system and a metal roofing system. The curtain wall system includes a metal curtain wall or glass curtain wall, a supporting structure, and a connection and sealing system. The supporting structure includes connectors or purlins or purlin brackets. The connection and sealing system includes silicone sealant, sealing strips or hinges, and bolts.

[0087] Furthermore, in S200, the method for simulating and obtaining the surface wind pressure distribution of the BIM design model using computer-aided design simulation software is as follows: import the BIM design model into the computer-aided design simulation software, and perform flow field analysis under given wind speed boundary conditions to obtain the surface wind pressure distribution of the BIM design model. The specific steps of the flow field analysis include: importing the BIM design model, surface segmentation, mesh generation, and mesh attribute definition; setting material parameters, geometric properties, wind speed boundary conditions, load conditions, and defining the connectivity of elements.

[0088] Preferably, the flow field analysis specifically involves: simulating transient wind pressure distribution under different wind angles (0°~360°, interval 22.5°) in the surrounding flow field domain of the BIM design model of the target building; coupling the wind pressure data to the structural finite element model; applying dynamic wind pressure loads; performing time history analysis and identifying resonance areas; and exporting the time history data (.csv format) of the wind pressure coefficient on the surface of the BIM design model, as well as exporting the wind pressure distribution data and displacement cloud map of the BIM design model.

[0089] Furthermore, the BIM design model is in 3dm format.

[0090] Furthermore, S100 also includes segmenting and extracting the 3D mesh of the wall surface of the BIM design model, calculating the 3D support plane fitted by the least squares method, calculating the Euclidean distance between the 3D mesh and the 3D support plane, and normalizing it to generate a depth map of the wall surface; generating a texture image of the wall surface through 3D spatial texture mapping based on the calibrated camera parameters; and ensuring that the length and width dimensions of the depth map and texture map of the wall surface are consistent.

[0091] Preferably, the method for mesh generation is the sweep method.

[0092] The given wind speed boundary conditions range from 5 m / s to 36.9 m / s.

[0093] Preferably, different wind speeds are set at different heights in the BIM design model for simulation. The given wind speed boundary conditions are simulated as follows: the wind speed in the area corresponding to the location of the BIM design model at 0-10 meters is set to 5 m / s, the wind speed at 10-30 meters is set to 8.7 m / s, the wind speed at 30-90 meters is set to 15 m / s, and the wind speed above 90 meters is set to 36.9 m / s.

[0094] Preferably, the computer-aided design simulation software is the open-source OpenFOAM software.

[0095] Furthermore, in S300, different wind speed conditions are wind speeds ranging from 5 m / s to 36.9 m / s.

[0096] Furthermore, in S300, the method for identifying and marking the structural installation locations of risk-sensitive areas on the BIM design model is as follows:

[0097] Based on the surface wind pressure distribution, obtain the wind pressure values ​​corresponding to all points at each structural installation location in the BIM design model. Record the vibration displacement of the geometric center of gravity at each structural installation location from the time of minimum wind pressure to the time of maximum wind pressure as the wind pressure limit displacement.

[0098] Record other structural installation positions adjacent to the current structural installation position as the adjacent structural installation positions of the current structural installation position;

[0099] Let MaxCUP be the geometric centroid of the structural installation position with the largest wind pressure limit displacement among all adjacent structural installation positions of the current structural installation position; let MinCUP be the geometric centroid of the structural installation position with the smallest wind pressure limit displacement among all adjacent structural installation positions of the current structural installation position; then CUPL is formed from MinCUP to MaxCUP; all structural installation positions that CUPL passes through in sequence are recorded as the wind pressure orientation region of the current structural installation position.

[0100] Let WV(j) be the j-th wind pressure orientation region at the current structural installation location; j is the sequence number, and the order of the wind pressure orientation regions is the order in which CUPL passes through the wind pressure orientation regions.

[0101] Determine each WV(j). If WV(j) satisfies the linear wind pressure condition, then mark WV(j) as the structural installation location of the risk-sensitive area.

[0102] The linear wind pressure condition is: MV(WV(j-1))>MV(WV(j)) and MV(WV(j))<MV(WV(j+1)), where MV(WV(j)) is the wind pressure limit displacement corresponding to WV(j).

[0103] The key source code in C# language describing the specific implementation of the method for identifying and marking the structural installation locations of risk-sensitive areas on the BIM design model in Embodiment 1 is as follows:

[0104] using System;

[0105] using System.Collections.Generic;

[0106] using System.Linq;

[0107] using UnityEngine; / / Requires the Unity engine's math library. Adjust accordingly if using another 3D math library.

[0108] public class CurtainWallRiskAnalyzer

[0109] { / / Structure installation location data structure

[0110] public class CurtainWallRegion

[0111] {public int RegionId { get; set;}

[0112] public Vector3 GeometricCenter { get; set;} / / Geometric center of gravity

[0113] public List <vector3>SurfacePoints { get; set;} / / Set of surface points

[0114] public float[] WindPressureValues ​​{ get; set;} / / Wind pressure values ​​at various points

[0115] public float MaxWindPressureDisplacement { get; set;} / / Wind pressure limit displacement

[0116] public List <int>AdjacentRegions { get; set;} / / List of adjacent region IDs

[0117] public bool IsAtRisk { get; set;} / / Risk marker}

[0118] / / / Identify and mark the installation locations of structures in risk-sensitive areas

[0119] / / / <param name="allRegions"> All structural installation location data

[0120] public static void IdentifyRiskRegions(List <curtainwallregion>(allRegions)

[0121] Step 1: Calculate the wind pressure limit displacement for all regions.

[0122] CalculateMaxDisplacements(allRegions);

[0123] / / Step 2: Determine the wind pressure direction for each area

[0124] Dictionary <int, List <curtainwallregion>>directionRegionsMap = new();

[0125] foreach (var region in allRegions)

[0126] {var directionRegions = GetWindDirectionRegions(region, allRegions);

[0127] directionRegionsMap.Add(region.RegionId, directionRegions);}

[0128] / / Step 3: Check linear wind pressure conditions and mark risks

[0129] foreach (var region in allRegions)

[0130] {if (directionRegionsMap.TryGetValue(region.RegionId, out varwvList))

[0131] {for (int j = 1; j <wvList.Count - 1; j++)

[0132] {if (CheckLinearCondition(wvList, j))

[0133] {wvList[j].IsAtRisk = true;Debug.Log($"Risk marker: Region {wvList[j].RegionId}");

[0134] }}}}}

[0135] / / / Calculate the ultimate displacement under wind pressure (the displacement of the center of gravity from the minimum wind pressure to the maximum wind pressure).

[0136] private static void CalculateMaxDisplacements(List <curtainwallregion>regions)

[0137] {foreach (var region in regions)

[0138] {if(region.WindPressureValues==null||region.WindPressureValues.Length==0)

[0139] continue?

[0140] / / Simulate displacement changes under different wind pressures

[0141] float minPressure = region.WindPressureValues.Min();

[0142] float maxPressure = region.WindPressureValues.Max();

[0143] region.MaxWindPressureDisplacement = (maxPressure - minPressure) *0.1f;

[0144] }}

[0145] / / / Get the wind pressure direction of the current region (the region traversed by the MinCUP->MaxCUP ray).

[0146] private static List <curtainwallregion>GetWindDirectionRegions(

[0147] CurtainWallRegion currentRegion,

[0148] List <curtainwallregion>(allRegions)

[0149] / / Find the extrema centroid in the adjacent region

[0150] var adjacentRegions = allRegions.Where(r =>

[0151] currentRegion.AdjacentRegions.Contains(r.RegionId)).ToList();

[0152] if (adjacentRegions.Count == 0)

[0153] return new List <curtainwallregion>();

[0154] var maxCUP = adjacentRegions.OrderByDescending(r =>r.MaxWindPressureDisplacement)

[0155] .First().GeometricCenter;

[0156] var minCUP = adjacentRegions.OrderBy(r =>r.MaxWindPressureDisplacement)

[0157] .First().GeometricCenter;

[0158] / / Construction Ray Direction

[0159] Vector3 rayDirection = (maxCUP - minCUP).normalized;

[0160] Ray cupRay = new Ray(minCUP, rayDirection);

[0161] / / Ray detection of all areas traversed

[0162] var piercedRegions = new List <curtainwallregion>();

[0163] foreach (var region in allRegions)

[0164] {if (RayIntersectsRegion(cupRay, region))

[0165] piercedRegions.Add(region);}

[0166] / / Sort by the order in which rays pass through (simplified to sort by centroid distance)

[0167] return piercedRegions.OrderBy(r =>

[0168] Vector3.Distance(minCUP, r.GeometricCenter)).ToList();}

[0169] / / / Detection of intersection between ray and structural installation location (using bounding boxes to simplify calculations)

[0170] private static bool RayIntersectsRegion(Ray ray, CurtainWallRegionregion)

[0171] {Bounds bounds = new Bounds(region.GeometricCenter, Vector3.zero);

[0172] foreach (var point in region.SurfacePoints)

[0173] bounds.Encapsulate(point);

[0174] return bounds.IntersectRay(ray);}

[0175] / / / Check the linear wind pressure condition MV(WV(j-1))>MV(WV(j)) <MV(WV(j+1))

[0176] private static bool CheckLinearCondition(List <curtainwallregion>wvRegions, int j)

[0177] {if (wvRegions.Count<3 || j<= 0 || j>= wvRegions.Count - 1)

[0178] return false;

[0179] float prevMV = wvRegions[j - 1].MaxWindPressureDisplacement;

[0180] float currentMV = wvRegions[j].MaxWindPressureDisplacement;

[0181] float nextMV = wvRegions[j + 1].MaxWindPressureDisplacement;

[0182] return prevMV>currentMV&¤tMV <nextMV;}}。

[0183] Furthermore, in S500, during the main structure installation stage, the enclosure structure frame installation stage, and the enclosure structure panel installation stage, the accuracy of the corresponding positions of the structural installation locations in risk-sensitive areas between the real-world 3D model and the BIM design model is verified. The method for generating deviation data between the real-world model and the BIM design model through registration deviation analysis is as follows:

[0184] The ICP algorithm is used to verify the accuracy of the corresponding positions of the structural installation locations in risk-sensitive areas between the real-scene 3D model and the BIM design model. The deviation data between the real-scene model and the BIM design model is generated through registration deviation analysis.

[0185] The ICP algorithm-based fusion analysis involves registering the real-world 3D model with the BIM design model and calculating the deviation data for the corresponding structural installation locations in risk-sensitive areas. The deviation data is displayed using color coding to show the dimensional deviation ranges for different areas.

[0186] The deviation data includes overall deviation statistics (mean, standard deviation, maximum deviation), deviation heatmap (color mapping distance values, visually displaying error distribution), and structural installation location error analysis data for structures in risk-sensitive areas: deviations are grouped and statistically analyzed according to structures (external envelope structures) in risk-sensitive areas.

[0187] The ICP algorithm is the ICP registration algorithm in the open-source 3D data processing library Open3D.

[0188] Furthermore, in the S500, the method for outputting and displaying the real-world 3D model is as follows: adaptive histogram equalization is performed on the locations of deviations in the real-world 3D model to enhance local contrast.

[0189] Preferably, it also includes transmitting the real-scene 3D model to a client or mobile device for display output.

[0190] Preferably, the method for transmitting a real-scene 3D model to a client or mobile device for display output includes the following steps: processing the real-scene 3D model using an edge-folding algorithm to obtain a simplified model;

[0191] The simplified model is divided into 1m×1m blocks and compressed separately.

[0192] Dynamically push visible blocks based on the view of the client or mobile device (using HTTP / 2 ServerPush).

[0193] When loading on the client or mobile device, check GPU performance and select the corresponding LOD level:

[0194] High-end equipment: LOD0 (original precision);

[0195] Mid-range equipment: LOD1 (50% surface area);

[0196] Low-end equipment: LOD2 (20% surface area + solid color heatmap);

[0197] The installation locations of structures in risk-sensitive areas will be highlighted in real time.

[0198] Among them, high-end equipment: supports hardware-accelerated ray tracing, floating-point computing power ≥3 TFLOPS, and memory bandwidth ≥100GB / s;

[0199] Mid-range devices: Support Vulkan 1.1 and above API, with floating-point computing power of 1~3 TFLOPS;

[0200] Low-end devices: Only support OpenGL ES 3.0, floating-point computing power <1 TFLOPS.

[0201] Example 2

[0202] This second embodiment replaces the linear wind pressure condition with that of the first embodiment, specifically:

[0203] Preferably, the linear wind pressure condition is:

[0204] DE(WV(j-1))>DE(WV(j)) and DE(WV(j))<DE(WV(j+1)),

[0205] Where DE(WV(j)) is the inertial wind pressure index of the j-th wind pressure orientation region WV(j) at the current structural installation location;

[0206] Preferably, the inertial wind pressure index DE(WV(j)) is calculated as follows:

[0207] DE(WV(j))=|MinA1-MinA2|÷|MaxA1-MaxA2|;

[0208] Where MaxA1 is the maximum value of the wind pressure limit displacement corresponding to the wind pressure orientation region from the 1st to the (j-1th)th region; MinA1 is the minimum value of the wind pressure limit displacement corresponding to the wind pressure orientation region from the 1st to the (j-1th)th region; MaxA2 is the maximum value of the wind pressure limit displacement corresponding to the wind pressure orientation region from the (j+1th)th to the SuMth region; MinA2 is the minimum value of the wind pressure limit displacement corresponding to the wind pressure orientation region from the (j+1th)th to the SuMth region; and SuM is the number of wind pressure orientation regions.

[0209] The key source code in C# language describing the specific implementation of the linear wind pressure condition in Example 2 is as follows:

[0210] public class WindPressureOptimizer

[0211] Check if the linear wind pressure condition is met.

[0212] / / / DE(WV(j-1))>DE(WV(j)) and DE(WV(j))<DE(WV(j+1))

[0213] / / / <param name="deValues"> Array of inertial wind pressure indices for all regions (index = region number - 1)

[0214] / / / <param name="currentRegion"> Current region j (starting from 1)

[0215] public static bool IsLinearWindPressureCondition(double[] deValues,int currentRegion)

[0216] {if (deValues ​​== null || currentRegion<2 || currentRegion>=deValues.Length)

[0217] return false;

[0218] return deValues[currentRegion - 2]>deValues[currentRegion - 1]&&

[0219] deValues[currentRegion - 1] <deValues[currentRegion];}

[0220] / / / Calculate the improved inertial wind pressure index DE(WV(j)) = |MinA1-MinA2| / |MaxA1-MaxA2|

[0221] / / / <param name="windLoadLimits"> Wind pressure limit displacement array (index = region number - 1)

[0222] / / / <param name="windDirectionRegions"> Wind pressure orientation area dictionary (Key=area number, Value=list of corresponding orientation areas)

[0223] / / / <param name="currentRegion"> Current region j (starting from 1)

[0224] public static double CalculateEnhancedInertialIndex(

[0225] double[] windLoadLimits,

[0226] Dictionary <int, List <int>windDirectionRegions,

[0227] int currentRegion)

[0228] {if (windLoadLimits == null || windLoadLimits.Length == 0)

[0229] throw new ArgumentException("Invalid wind pressure limit displacement data");

[0230] if (currentRegion<1 || currentRegion>windLoadLimits.Length)

[0231] throw new ArgumentException("Region ID is out of range");

[0232] / / Get all wind pressure directions in the current area

[0233] var currentDirectionRegions = windDirectionRegions[currentRegion];

[0234] if (currentDirectionRegions == null || currentDirectionRegions.Count== 0)

[0235] return 0;

[0236] / / Step 1: Extract the displacement values ​​belonging to the current orientation region from the first j-1 regions.

[0237] var previousRegions = Enumerable.Range(0, currentRegion - 1)

[0238] .Where(i =>windDirectionRegions[i + 1].Intersect(currentDirectionRegions).Any()).Select(i =>windLoadLimits[i])

[0239] .ToList();

[0240] / / Step 2: Extract the displacement values ​​belonging to the current orientation region from j+1 to SuM regions.

[0241] var nextRegions = Enumerable.Range(currentRegion,windLoadLimits.Length - currentRegion).Where(i=>windDirectionRegions[i+ ].Intersect(currentDirectionRegions).Any()).Select(i =>windLoadLimits[i]).ToList();

[0242] / / Calculate extreme values

[0243] double minA1 = previousRegions.Any() ? previousRegions.Min() : 0;

[0244] double maxA1 = previousRegions.Any() ? previousRegions.Max() : 0;

[0245] double minA2 = nextRegions.Any() ? nextRegions.Min() : 0;

[0246] double maxA2 = nextRegions.Any() ? nextRegions.Max() : 0;

[0247] / / Prevent the denominator from being zero (returns 0 or a specific default value when MaxA1 == MaxA2)

[0248] double denominator = Math.Abs(maxA1 - maxA2);

[0249] if (denominator < 1e-10) / / Floating-point precision tolerance

[0250] return 0;

[0251] return Math.Abs(minA1 - minA2) / denominator;}

[0252] public static List <int>FindOptimalRegions(double[] windLoadLimits, Dictionary<int, List <int>> directionRegions)

[0253] {var deValues = new double[windLoadLimits.Length];

[0254] for (int j = 1; j<= windLoadLimits.Length; j++)

[0255] {deValues[j-1]=CalculateEnhancedInertialIndex(windLoadLimits,directionRegions, j);}

[0256] var optimalRegions = new List <int>();

[0257] for (int j = 2; j<windLoadLimits.Length; j++)

[0258] {if (IsLinearWindPressureCondition(deValues, j))

[0259] optimalRegions.Add(j);}

[0260] return optimalRegions;}}.

[0261] Example 3

[0262] This example 3 is based on the calculation method of the inertial wind pressure index DE(WV(j)) in example 1 is replaced, specifically:

[0263] Preferably, the calculation method of the inertial wind pressure index DE(WV(j)) is:

[0264] Calculate the convergence wind pressure difference and the transfer wind pressure difference of all structure installation positions, the specific method is:

[0265] The difference between the average value of the wind pressure limit displacement corresponding to the wind pressure orientation area from the first to the j-1th structure installation position and the corresponding wind pressure limit displacement of the structure installation position is the convergence wind pressure difference; the difference between the average value of the wind pressure limit displacement corresponding to the wind pressure orientation area from the j+1th to the Nth structure installation position and the corresponding wind pressure limit displacement of the structure installation position is the transfer wind pressure difference; (the transfer wind pressure difference represents the degree of increase of the wind pressure limit displacement along the wind pressure orientation area from the structure installation position, and the convergence wind pressure difference represents the degree of reduction of the wind pressure limit displacement along the wind pressure orientation area from the structure installation position);

[0266] Traverse the wind pressure orientation area corresponding to each structure installation position, MaxS represents the number of wind pressure orientation areas whose convergence wind pressure difference is greater than the transfer wind pressure difference in each wind pressure orientation area corresponding to the structure installation position; MinS represents the number of wind pressure orientation areas whose convergence wind pressure difference is less than the transfer wind pressure difference in each wind pressure orientation area corresponding to the structure installation position; then the calculation method of the inertial wind pressure index DE(WV(j)) is: inertial wind pressure index DE(WV(j))= MaxS÷(MaxS+MinS).

[0267] Wherein, the key source code of the specific implementation of the calculation method of the inertial wind pressure index DE(WV(j)) of this example 3 in C# language is:

[0268] public class WindPressureCalculator

[0269] Calculate the inertial wind pressure index DE(WV(j)) at the specified structural installation location j.

[0270] / / / <param name="windLoadLimits"> An array of wind pressure limit displacement values ​​for all structural installation locations (subscript = area number - 1).

[0271] / / / <param name="windDirectionRegions"> Wind pressure orientation zone dictionary for all structural installation locations (Key=zone number, Value=list of corresponding orientation zones)

[0272] / / / <param name="j"> The current calculated structural installation location number (starting from 1).

[0273] public static double CalculateInertialIndex(

[0274] double[] windLoadLimits,

[0275] Dictionary <int, List <int>>windDirectionRegions,int j)

[0276] { / / Parameter validation

[0277] if (windLoadLimits == null || windLoadLimits.Length == 0)

[0278] throw new ArgumentException("Invalid wind pressure limit displacement array");

[0279] if (j<1 || j>windLoadLimits.Length)

[0280] throw new ArgumentException("Structure installation location number is out of range");

[0281] int maxS = 0; / / Number of oriented regions where the convergent wind pressure difference is greater than the propagating wind pressure difference

[0282] int minS = 0; / / Converging wind pressure difference < Number of oriented regions transmitting wind pressure difference

[0283] / / Iterate through all wind pressure orientation areas in the current region j

[0284] foreach (int directionRegion in windDirectionRegions[j])

[0285] / / Calculate the convergent wind pressure difference (the difference between the average of the previous j-1 regions and the current region)

[0286] double convergenceDiff = CalculateConvergenceDifference(windLoadLimits, j, directionRegion);

[0287] / / Calculate the transmitted wind pressure difference (the average difference between the current region and the region from j+1 to N).

[0288] double transmissionDiff = CalculateTransmissionDifference(windLoadLimits, j, directionRegion);

[0289] / / Statistics on MaxS and MinS

[0290] if (convergenceDiff>transmissionDiff)maxS++;

[0291] else if (convergenceDiff<transmissionDiff)minS++;}

[0292] / / Calculate the inertial wind pressure index (avoid division by zero)

[0293] return (maxS + minS) == 0? 0 : (double)maxS / (maxS + minS);}

[0294] / / / Calculate the convergence wind pressure difference (the average of the wind pressure limit displacement of the previous j-1 regions and the difference with the current region j)

[0295] private static double CalculateConvergenceDifference(

[0296] double[] windLoadLimits,

[0297] int currentRegion,

[0298] int directionRegion)

[0299] {double sum = 0;int count = 0;

[0300] / / Add the displacement values of the previous 1 to j-1 regions that belong to the current direction region

[0301] for (int i = 0; i<currentRegion - 1; i++) / / Note that the array index starts from 0

[0302] {sum += windLoadLimits[i];count++;}

[0303] double average = count>0? sum / count : 0;

[0304] return average - windLoadLimits[currentRegion - 1]; / / Note that the array index is converted

[0305] / / / Calculate transmission difference (average of wind pressure limit displacement of regions from j+1 to N minus current region j)

[0306] private static double CalculateTransmissionDifference(

[0307] double[] windLoadLimits,

[0308] int currentRegion, int directionRegion)

[0309] {double sum = 0;int count = 0;

[0310] / / Accumulate displacement values of regions from j+1 to end that belong to current direction region

[0311] for (int i = currentRegion; i<windLoadLimits.Length; i++) / / i corresponds to region number i+1

[0312] { / / Also need to judge whether region i+1 belongs to the current direction region

[0313] sum += windLoadLimits[i];count++;}

[0314] double average = count>0? sum / count : 0;

[0315] return average - windLoadLimits[currentRegion - 1];}}.

[0316] Further, the building engineer can mark the corresponding position of the position with deviation on the real three-dimensional model to perform any one of the following constructions: setting a telescopic joint node, using a deflectable connecting piece, and selecting a high-elasticity sealing material as a reinforced sealant.

[0317] Although the description of the present disclosure has been quite detailed and particularly described with respect to several described embodiments, it is not intended to be limited to any of these details or embodiments or any special embodiment, so as to effectively cover the intended scope of the present disclosure. In addition, the present disclosure is described above in embodiments that the inventor can foresee, and the purpose is to provide a useful description, and non-essential modifications to the present disclosure that have not yet been foreseen can still represent equivalent modifications of the present disclosure.< / int> < / int> < / int> < / int> < / int> < / curtainwallregion> < / curtainwallregion> < / curtainwallregion> < / curtainwallregion> < / curtainwallregion> < / curtainwallregion> < / curtainwallregion> < / curtainwallregion> < / int>

Claims

1. A method for detecting the installation of irregularly shaped, multi-curved building envelope structures using UAV 3D reconstruction, characterized in that... The method includes the following steps: S100: Import the BIM design model of the target building and mark the location of the external envelope structure on the BIM design model as the structural installation location. S200 uses computer-aided design simulation software to simulate and obtain the surface wind pressure distribution of the BIM design model; S300, under different wind speed conditions, loads the surface wind pressure distribution onto the BIM design model to obtain the vibration displacement of each structural installation location under different average wind speeds, and identifies and marks the structural installation locations with risk-sensitive areas on the BIM design model. The S400 uses a drone to capture a sequence of images of a target building, and then performs 3D reconstruction on the image sequence to obtain a realistic 3D model of the main structure of the target building. S500 performs accuracy verification on the corresponding positions of the structural installation locations in risk-sensitive areas of the real-world 3D model and the BIM design model during the main structure installation stage, the enclosure structure frame installation stage, and the enclosure structure panel installation stage. It also generates deviation data between the real-world model and the BIM design model through registration deviation analysis. S600 marks the locations with deviations on the real-world 3D model, outputs and displays the real-world 3D model; In S300, the specific method for identifying and marking the structural installation locations of risk-sensitive areas on the BIM design model is as follows: Mark the wind pressure orientation area in the adjacent structural installation locations of the current structural installation location; Based on linear wind pressure conditions, structural installation locations in risk-sensitive areas are selected within the wind pressure orientation region; The linear wind pressure condition is as follows: Let WV(j) be the j-th wind pressure orientation region at the current structural installation position; j is the sequence number, MV(WV(j-1))>MV(WV(j)) and MV(WV(j))<MV(WV(j+1)), where MV(WV(j)) is the wind pressure limit displacement corresponding to WV(j).

2. The method for installation and inspection of irregularly shaped multi-curved building envelope structures reconstructed by UAV according to claim 1, characterized in that, The method for marking the wind pressure orientation area in the adjacent structural installation locations of the current structural installation location is as follows: Based on the surface wind pressure distribution, obtain the wind pressure values ​​corresponding to all points at each structural installation location in the BIM design model. Record the vibration displacement of the geometric center of gravity at each structural installation location from the time of minimum wind pressure to the time of maximum wind pressure as the wind pressure limit displacement. Record other structural installation positions adjacent to the current structural installation position as the adjacent structural installation positions of the current structural installation position; Let MaxCUP be the geometric centroid of the structural installation position with the largest wind pressure limit displacement among all adjacent structural installation positions of the current structural installation position; let MinCUP be the geometric centroid of the structural installation position with the smallest wind pressure limit displacement among all adjacent structural installation positions of the current structural installation position; then a ray CUPL is formed from MinCUP to MaxCUP; all structural installation positions that CUPL passes through in sequence are recorded as the wind pressure orientation region of the current structural installation position.

3. The method for installation and inspection of irregularly shaped multi-curved building envelope structures reconstructed by UAV according to claim 2, characterized in that, Replace the linear wind pressure condition with: DE(WV(j-1))>DE(WV(j)) and DE(WV(j))<DE(WV(j+1)), Where DE(WV(j)) is the inertial wind pressure index of the j-th wind pressure orientation region WV(j) at the current structural installation location.

4. The method for installation and inspection of irregularly shaped multi-curved building envelope structures reconstructed by UAV according to claim 3, characterized in that, The inertial wind pressure index DE(WV(j)) is calculated as follows: DE(WV(j))=|MinA1-MinA2|÷|MaxA1-MaxA2|; Where MaxA1 is the maximum value of the wind pressure limit displacement corresponding to the wind pressure orientation region from the 1st to the (j-1th)th region; MinA1 is the minimum value of the wind pressure limit displacement corresponding to the wind pressure orientation region from the 1st to the (j-1th)th region; MaxA2 is the maximum value of the wind pressure limit displacement corresponding to the wind pressure orientation region from the (j+1th)th to the SuMth region; MinA2 is the minimum value of the wind pressure limit displacement corresponding to the wind pressure orientation region from the (j+1th)th to the SuMth region; and SuM is the number of wind pressure orientation regions.

5. The method for installation and inspection of irregularly shaped multi-curved building envelope structures reconstructed by UAV according to claim 4, characterized in that, The calculation method for the inertial wind pressure index DE(WV(j)) is replaced with: The specific method for calculating the convergent and transmitted air pressure differences at all structural installation locations is as follows: The difference between the average value of the wind pressure limit displacement corresponding to the wind pressure orientation area from the 1st to the (j-1th)th structural installation position and the corresponding wind pressure limit displacement of the structural installation position is the convergent wind pressure difference. The difference between the average value of the wind pressure limit displacement corresponding to the wind pressure orientation area from the (j+1)th to the Nth structural installation position and the wind pressure limit displacement corresponding to the structural installation position is used as the transmitted wind pressure difference. Traverse the wind pressure orientation regions corresponding to each structural installation location. Let MaxS represent the number of wind pressure orientation regions where the convergent wind pressure difference is greater than the transmitted wind pressure difference; let MinS represent the number of wind pressure orientation regions where the convergent wind pressure difference is less than the transmitted wind pressure difference; then the inertial wind pressure index DE(WV(j)) is calculated as follows: Inertial wind pressure index DE(WV(j)) = MaxS ÷ (MaxS + MinS).

6. The method for installation and inspection of irregularly shaped multi-curved building envelope structures reconstructed by UAV according to claim 1, characterized in that, In the S500, the method for outputting and displaying the real-world 3D model is to perform adaptive histogram equalization on the locations of deviations in the real-world 3D model to enhance local contrast.

7. The method for installation and inspection of irregularly shaped multi-curved building envelope structures reconstructed by UAV according to claim 1, characterized in that, It also includes transmitting the real-world 3D model to the client or mobile device for display output.

8. The method for installation and inspection of irregularly shaped multi-curved building envelope structures reconstructed by UAV according to claim 7, characterized in that, The method for transmitting a real-scene 3D model to a client or mobile device for display output includes the following steps: processing the real-scene 3D model through an edge-folding algorithm to obtain a simplified model; The simplified model is divided into spatial blocks, and each block is compressed separately. Dynamically push visible blocks based on the field of view of the client or mobile device; When loading on a client or mobile device, the GPU performance is detected and the corresponding LOD level is selected.

Images