Method and system for automatically generating three-dimensional coordinates of surface measurement points of large-scale urban building complexes

Abstract

The invention discloses a method and system for automatically generating three-dimensional coordinates of surface measuring points of a large-scale urban building complex, which solves the problems of low efficiency, inability to realize regular and uniform arrangement, and poor arrangement adaptability of existing measuring point layout methods, and specifically comprises the following steps: obtaining effective building GIS data and projection coordinate system code; parsing building attribute information and outline data; transforming the projection coordinate system; creating building outline data and binding building attribute information; calculating all outline side lengths, orientations and midpoint coordinates of any building one by one, presetting the horizontal measuring point spacing, and calculating the plane coordinates of all measuring points on the building outline of the selected building by using a linear interpolation function until all outlines of the building are traversed; presetting the vertical measuring point spacing, and calculating the coordinates of all measuring points in the vertical direction from bottom to top until all floors of the building are traversed to obtain the three-dimensional coordinates of all measuring points; and finally obtaining the three-dimensional coordinates of all building surface measuring points.

Description

Technical Field

[0001] The invention relates to a system and method for automatically generating coordinates, and in particular to a method and system for automatically generating three-dimensional coordinates of surface measurement points of a large-scale urban building complex. Background Art

[0002] Cities are spaces where various elements are highly concentrated. With the gathering of a large population, resources and wealth, the development level of modern cities has become an important indicator of a country's comprehensive strength. Within urban space, large-scale building complexes are the key material carriers to ensure the normal operation of the city, and are also the main damage targets of air shock waves generated by medium- and low-altitude / ground-touching large-yield explosions. Conducting research on the damage assessment technology of explosion shock waves on large-scale urban building complexes and building a scientific, efficient and intelligent urban building complex damage effect calculation and assessment technology system are of great significance to urban disaster prevention and mitigation and emergency plan formulation.

[0003] However, in the above research process, there are two problems that need to be solved first: one is to develop a suitable method to calculate the shock wave transmission environment inside the complex urban space, and the other is to establish a method that can accurately, efficiently and automatically extract the explosion loads borne by thousands or even hundreds of thousands of buildings on this basis. The former will give the calculation results of the shock wave flow field environment, and the latter will complete the batch extraction of the time history data of the explosion load on the surface of the building group from the calculation results. If the building group load data cannot be fully and accurately extracted, the accuracy of subsequent series of work such as building load distribution, structural response calculation and damage assessment, explosion resistance and reinforcement design will be seriously restricted.

[0004] For the first problem mentioned above, considering that modern cities usually contain a large number of building targets with different shapes and sizes and different arrangements, the distribution of shock wave flow field and its evolution over time will be extremely complex. Traditional theoretical calculation methods or empirical formulas, charts, etc. of various shock wave environmental parameters can no longer accurately describe the above-mentioned global and all-time flow field information. At present, computational fluid dynamics (CFD) is the main tool to solve the above problems, that is, using various CFD programs to directly model and calculate the three-dimensional complex urban space containing a large number of building complexes, and then reproduce the transmission evolution process of shock waves in urban building complexes. In recent years, thanks to the rapid development of high-performance computing software and hardware, the scale and complexity of urban space explosion transmission problems that can be handled by CFD programs have also increased dramatically. It is now possible to carry out refined flow field calculations for urban spaces with a scale of tens or even hundreds of square kilometers and covering tens of thousands or even more buildings.

[0005] However, the above research only stays at the stage of calculating the flow field transmission environment, that is, the first problem mentioned above, and does not involve the second problem. Although the existing research can calculate the explosion shock wave in the combined space of buildings with different arrangements, different shapes and different sizes, the calculation results are mainly roughly qualitatively displayed in the form of cloud maps, and lack the means to batch and automatically extract the explosion loads borne by the surfaces of various real buildings within the calculation range. In fact, the prerequisite for extracting the explosion load history around the building complex is that the generation of the three-dimensional coordinates of all measuring points must be completed in the CFD pre-processing stage. The current method is to manually select the position one by one in the graphical user interface (GUI) of the pre-processing modeling software, and then obtain the three-dimensional coordinate value at the position. This method is completely dependent on manual operation and is only applicable to obtaining the position coordinates of a very small number of measuring points in a single building or several buildings. Obviously, it is impossible to carry out the automatic layout of large-scale measuring points. Especially when facing the scene of a super-large scale of hundreds of square kilometers, including hundreds of thousands of buildings with different geometric shapes and sizes, the above method has three obvious disadvantages: ① The layout of measuring points is inefficient and the number is very limited. If we want to fully and accurately reflect the temporal and spatial distribution characteristics of the load of urban buildings and support subsequent response calculations, thousands of measurement points may be required for only one conventional high-rise building. For a city, there are usually hundreds of thousands of buildings, so the number of measurement points required is extremely large. The existing method of obtaining the coordinates of the measurement points will be unbearable in terms of time alone. ② The measurement points cannot be arranged regularly and evenly, and the success rate is difficult to guarantee. When manually selecting the location, except for the relatively accurate positioning at the vertex of the building outline, the arrangement of measurement points on the building facade depends entirely on human subjective judgment. The intervals between the measurement points are uneven, and the regular and uniform effect cannot be achieved, which is not conducive to the subsequent load distribution and response research. In addition, the coordinates of the picked measurement points are often not strictly attached to the outer surface of the building. It is common for them to be far away from or invade the interior of the building model, resulting in the failure of the measurement point load extraction. ③ The adaptability of the measurement point layout is poor, and it cannot change with the changes in the calculation scene. When the target city or the calculation range changes, the existing measurement point coordinates will become invalid and need to be re-acquired one by one. Summary of the invention

[0006] The purpose of the present invention is to provide a method and system for automatically generating three-dimensional coordinates of surface measurement points of large-scale urban building complexes, so as to solve the technical problems of low efficiency and very limited number of measurement point layout methods, inability to achieve regular and uniform layout, and poor layout adaptability, so as to support the automatic and batch extraction of surface load history data of building complexes when calculating CFD flow field environments in complex urban spaces.

[0007] In order to achieve the above object, the present invention adopts the following technical solution:

[0008] A method for automatically generating three-dimensional coordinates of surface measurement points of a large-scale urban building complex is special in that it includes the following steps:

[0009] Step 1: Obtain original building GIS data, and select valid building GIS data from the original building GIS data;

[0010] Step 2: Get the projection coordinate system code;

[0011] Step 3: For each building in the valid building GIS data selected in step 1, the building attribute information and building outline data are analyzed in turn; the building attribute information includes the number of floors and floor height information of each building;

[0012] Step 4: Based on the building outline data parsed in step 3 and the projection coordinate system code obtained in step 2, the projection coordinate system is transformed to obtain the position coordinate data of the building outline data in the Cartesian coordinate system;

[0013] Step 5: Select the minimum value of all vertex coordinates in the projection coordinate system as the coordinate origin of the Cartesian coordinate system, or directly reset the origin coordinates of the Cartesian coordinate system in step 4 to zero, and synchronously move the position coordinate data of all building outlines obtained in step 4;

[0014] Step 6: for each building in the position coordinate data after the move in step 5, create building outline data and bind building attribute information in turn;

[0015] Step 7: Select one of the buildings, calculate the side length, orientation and midpoint coordinates of all the contours of the building one by one, preset the horizontal measurement point spacing, and use the linear interpolation function to calculate the plane coordinates of all the measurement points on the building contour of the selected building until all the contours of the building are traversed;

[0016] Step 8: Preset the vertical distance between measuring points, combine the number of floors and floor height information of the building selected in step 7, calculate the vertical coordinates of all measuring points from bottom to top, until all floors of the building are traversed, and obtain the three-dimensional coordinates of all measuring points;

[0017] Step 9, return to step 7, and select another building until all buildings are traversed to obtain the three-dimensional coordinates of all building surface measurement points, and complete the automatic generation of three-dimensional coordinates of surface measurement points of large-scale urban building complexes.

[0018] Furthermore, in step 1, original building GIS data is obtained from the geophysical information system data of the target city, and the building GIS data is two-dimensional vector data in shapefile format.

[0019] Furthermore, in step 1, the effective building GIS data refers to the original building GIS data including the two-dimensional plane geometric outline, location coordinates, number of building floors and building height of the target city building.

[0020] Furthermore, in step 3, the building outline data of each building is analyzed in turn as follows:

[0021] The latitude and longitude coordinate values ​​of the plane outline vertices, the area of ​​the closed polygon and the geometric center of each building in the valid building GIS data selected in step 1 are analyzed in turn and numbered; step 7 is specifically as follows:

[0022] Select the first building in the order of numbers, calculate the side length, orientation and midpoint coordinates of all the contours of the building one by one, preset the horizontal measurement point spacing, and use the linear interpolation function to calculate the plane coordinates of all the measurement points on the building contour of the selected building until all the contours of the building are traversed;

[0023] Step 9 is as follows:

[0024] Return to step 7, select another building in the order of the numbers, until all buildings are traversed, and the three-dimensional coordinates of all building surface measurement points are obtained, thus completing the automatic generation of three-dimensional coordinates of surface measurement points of large-scale urban building complexes.

[0025] Furthermore, in step 6, the building outline data is created specifically as follows:

[0026] Regenerate the building GIS data in the plane rectangular coordinate system in the order of point-line-surface-attribute information.

[0027] Furthermore, the method further comprises step 10:

[0028] Bind the three-dimensional coordinates of all measurement points obtained in step 9 with the attribute information of the corresponding buildings to obtain more complete data.

[0029] A system for automatically generating three-dimensional coordinates of surface measurement points of large-scale urban building complexes is used to execute the above-mentioned method for automatically generating three-dimensional coordinates of surface measurement points of large-scale urban building complexes. The system is special in that it includes an original building GIS data module, a user master control parameter module, a building GIS data parsing module, a projection transformation module, a building GIS data creation module, and a measurement point batch generation module.

[0030] The original building GIS data module is used to obtain the original building GIS data and select the valid building GIS data in the original building GIS data;

[0031] The user master control parameter module is used to obtain the projection coordinate system code, the preset horizontal measurement point spacing, and the preset vertical measurement point spacing;

[0032] The building GIS data analysis module is used to analyze the building attribute information and building outline data of each building in the acquired valid building GIS data in turn;

[0033] The projection transformation module transforms the projection coordinate system based on the parsed building outline data and the projection coordinate system code to obtain the position coordinate data of the building outline data in the Cartesian coordinate system; and is used to select the minimum value of all vertex coordinates as the coordinate origin of the Cartesian coordinate system or directly return the origin coordinate of the Cartesian coordinate system to zero, and simultaneously move the position coordinates of all building outlines synchronously;

[0034] The building GIS data creation module is used to create building outline data and bind building attribute information for each building in the moved position coordinate data in sequence; the building attribute information includes the number of floors and floor height information of each building;

[0035] The measuring point batch generation module is used to calculate the side length, orientation and midpoint coordinates of all the contours of each building one by one, and use a linear interpolation function to calculate the plane coordinates of all the measuring points on the building contour of the selected building until all the contours of the building are traversed; it is also used to combine the number of floors and floor height information of the building, calculate the coordinates of all the measuring points in the vertical direction from bottom to top, until all floors of the building are traversed, and the three-dimensional coordinates of all the measuring points are obtained; and then the three-dimensional coordinates of the measuring points on the surface of each building are obtained.

[0036] Furthermore, it also includes a data automatic output module;

[0037] The automatic data output module is used to automatically save and output the results generated by the building GIS data creation module and the measurement point batch generation module.

[0038] Beneficial effects of the present invention:

[0039] 1. The method and system provided by the present invention only need to specify the distance between measuring points in the horizontal direction and the height direction. According to these two control parameters, the regular, uniform and automatic layout of a large number of measuring points (more than 100,000) in a large-scale building complex (more than 100 buildings) can be achieved, and the three-dimensional coordinates of all measuring points can be obtained, and the layout quality of measuring points is high. The method is not limited by factors such as the scale of urban building complexes, the geometric shape of building complexes, and the arrangement method. In principle, there is no upper limit on the number of measuring points, which mainly depends on the size of the computer hard disk and the rounding error of the computer floating point operation.

[0040] 2. The present invention effectively solves the current technical problem of being unable to automatically and batch extract relevant load history data on the surface of building complexes in the calculation of shock wave flow field transmission environment in complex urban spaces, and builds a bridge between the detailed calculation of CFD flow field in complex urban spaces and the analysis of building structure response, providing important tool support for a series of subsequent studies on building load distribution, structural response calculation and damage assessment, explosion resistance / wind resistance and reinforcement design.

[0041] 3. The method and system provided by the present invention are highly efficient in arranging measuring points, and have low requirements on computer performance. Using an ordinary office computer, the arrangement of 100,000 measuring points can be completed in minutes. It has been verified that all measuring points are strictly attached to the outer surface of the building, and the success rate of measuring point load extraction can reach 100%.

[0042] 4. The method and system provided by the present invention are not only suitable for extracting explosion loads of complex urban space buildings, but also can be used for extracting wind loads of buildings in urban wind field CFD calculations. It can be used not only as a pressure measurement point, but also as a measurement point for various structural dynamic response parameters in related structural dynamics calculation programs, such as displacement, velocity, acceleration and strain.

[0043] 5. The method and system provided by the present invention can also generate measurement points and provide them to the structural dynamics calculation program, that is, they can be used as structural response measurement points such as displacement / velocity / acceleration / strain, etc. BRIEF DESCRIPTION OF THE DRAWINGS

[0044] Figure 1 It is a structural schematic diagram of a first embodiment of a system for automatically generating three-dimensional coordinates of surface measurement points of large-scale urban building complexes according to the present invention;

[0045] Figure 2 It is a flow chart of Embodiment 1 of the method for automatically generating three-dimensional coordinates of surface measurement points of large-scale urban building complexes of the present invention;

[0046] Figure 3 This is a comparison diagram of the building GIS data before and after the projection transformation in the first embodiment of the present invention, wherein the left side is the outline diagram of 8 multi-story buildings in the spherical coordinate system, and the right side is the plan diagram of 8 multi-story buildings in the Cartesian coordinate system;

[0047] Figure 4 It is a display diagram of the data file automatically generated in the first embodiment of the present invention, wherein (a) is a display diagram of the regenerated Shapefile data format file; (b) is a display diagram of the data file containing all the measurement point layout information;

[0048] Figure 5 The following are effect diagrams of measuring point layouts with different spacings in Example 1 of the present invention; (a) is an effect diagram of measuring point spacing of 1.0 m (a total of 2937 points); (b) is an effect diagram of measuring point spacing of 2.5 m (a total of 372 points);

[0049] Figure 6 It is 2.5km in a certain city in the second embodiment of the present invention. 2 The overall layout of buildings in the area;

[0050] Figure 7 This is a display diagram of the measurement point layout information data file in the second embodiment of the present invention;

[0051] Figure 8 This is a rendering of the arrangement of surface measurement points around some buildings in a large-scale urban building complex in the second embodiment of the present invention;

[0052] Fig. 9 is a location diagram of four typical measuring points in the building complex in the second embodiment of the present invention;

[0053] Fig.10 is a graph of the time history of the explosion load at the measuring point A in the large equivalent explosion scene in the second embodiment of the present invention;

[0054] Fig.11 It is a time history curve diagram of explosion load at measuring point B, measuring point C and measuring point D in a large equivalent explosion scene in the second embodiment of the present invention. DETAILED DESCRIPTION

[0055] The following will be combined with the drawings in the embodiments of the present invention to clearly and completely describe the technical solutions in the embodiments of the present invention. Obviously, the described embodiments are only part of the embodiments of the present invention, not all of the embodiments. Based on the embodiments of the present invention, all other embodiments obtained by ordinary technicians in this field without creative work are within the scope of protection of the present invention.

[0056] Example 1: Layout of measurement points in 8 multi-storey buildings

[0057] This embodiment aims at the problem that the surface load history data of building groups is difficult to extract automatically and in batches in the calculation of CFD flow field transmission environment in complex urban space. It innovatively proposes a system for automatically generating three-dimensional coordinates of surface measurement points of large-scale urban building groups. Figure 1 As shown in the figure, the system mainly consists of three parts: input layer, calculation layer and output layer, which contains 7 functional modules in total. The 7 functional modules are closely linked and interconnected, and finally achieve the purpose of batch layout of surface measurement points of large-scale urban buildings and automatic generation of three-dimensional coordinates. The functions of each module and the relationship between modules are as follows:

[0058] (1) The input layer contains all the input data required by the present invention and is composed of two functional modules: the original building GIS data module ① and the user master control parameter module ②.

[0059] The original building GIS data module ① reads the geophysical information system (GIS) data of the target urban building obtained through public or special channels. In this embodiment, two-dimensional vector data in shapefile format is used. This data type not only contains the two-dimensional plane geometric outline of the urban building and its position coordinates, but also contains other building attribute information such as building height in the form of an attribute list. For the original building GIS data, its building outline is represented by a closed vector polygon, and the vertex coordinates of the polygon are expressed in longitude and latitude.

[0060] The user master parameter module ② contains specific parameters such as the projection transformation coordinate system code, the preset horizontal measurement point spacing, and the vertical measurement point spacing, without other input. The projection transformation coordinate system is used to determine the method of transforming the target city building GIS data in the spherical coordinate system into the plane rectangular coordinate system, which varies from city to city. For a specific city, there may be multiple options for the projection transformation coordinate system, but the most suitable projection coordinate system can usually be found based on the longitude and latitude coordinates of the city combined with surveying and mapping knowledge, thereby reducing the distortion or distortion of the building outline before and after the projection transformation. The other two control parameters (horizontal and vertical measurement point spacing) are parameters in the Cartesian coordinate system after the urban building GIS projection transformation. They determine the spacing between measurement points and are used to control the sparseness or density of the measurement point layout, which is determined by the user according to needs.

[0061] (2) The calculation layer is the core calculation link involved in batch layout of measurement points and generation of three-dimensional coordinates of measurement points. It consists of four functional modules: building GIS data analysis module ③, projection transformation module ④, building GIS data creation module ⑤ and measurement point batch generation module ⑥.

[0062] Building GIS data parsing module ③ is a parsing module for building GIS data. It is mainly used to parse the latitude and longitude coordinate values ​​of the building plane outline vertices, the area of ​​the closed polygon, the geometric center and other information in the original building GIS data in sequence, and number and classify them for storage. The parsed data stream is used as the input condition for the next stage.

[0063] The projection transformation module ④ will receive the projection transformation coordinate system code specified by the user from the input layer, and transform all the analyzed data from the spherical coordinate system into the plane rectangular coordinate system according to the existing projection transformation theory and method, and further take the minimum value of all vertex coordinates as the coordinate origin or directly return the origin coordinate to zero to complete the movement of all the building plane outline position coordinates.

[0064] Building GIS data creation module ⑤ is mainly used to create building GIS data in a plane rectangular coordinate system. This module uses the building polygon vertices after projection transformation as input conditions, adopts the bottom-up construction principle, and regenerates building GIS data in the order of point-line-surface-attribute information. Different from the original building GIS data, the newly generated building GIS data can be directly used to construct a three-dimensional CFD calculation model of a complex urban space.

[0065] The measurement point batch generation module ⑥ is a measurement point batch layout module. This module will receive the measurement point spacing size preset by the user from the input layer, and determine the number of measurement points required for each polygonal contour by calculating the length and orientation of the building polygon outline and combining the measurement point spacing parameters. The linear interpolation function is used to determine the coordinate value of each measurement point, and further generate the three-dimensional coordinates of the measurement points on the surface of the entire building based on the number of floors and building height information. Repeat the above process for each building, and the three-dimensional coordinates of the measurement points of all buildings included in the urban building GIS file can be generated.

[0066] (3) The output layer includes the automatic data output module ⑦, which is mainly used to automatically save and output the results generated by the building GIS data creation module ⑤ and the measurement point batch generation module ⑥ to the local hard disk. Among them, the results generated in the building GIS data creation module ⑤ will be saved in the Shapefile format, and the results generated by the measurement point batch generation module ⑥ will be saved in a commonly used data format file.

[0067] like Figure 2 As shown, the basic working process of the above system is as follows:

[0068] Step 1: Read the original building GIS data prepared in advance and determine whether it is valid. If not, terminate the operation; if valid, proceed to the next step.

[0069] Step 2: Read the projection coordinate system code from the user master module and save it, then proceed to the next step.

[0070] Step 3: Start executing the building GIS data analysis module. For each building, perform building attribute information analysis, building outline data analysis, etc. in turn until all buildings are traversed and enter the next step.

[0071] Step 4: Enter the projection transformation module. For the architectural GIS data parsed in step 3, combined with the projection coordinate system code entered in step 2, complete the transformation of the projection coordinate system, and after generating the coordinates in metric units under the Cartesian coordinate system, proceed to the next step.

[0072] Step 5: Determine the appropriate coordinate origin and complete the movement and transformation of all building outline vertex coordinates to proceed to the next step.

[0073] Step 6: Create the building GIS data module. For each building, create the building outline data, bind the building attribute information, and so on. After all the buildings are traversed, proceed to the next step. The created building GIS data is automatically output and saved locally.

[0074] Step 7: Enter the measurement point batch generation module. For one of the buildings, calculate the side length, orientation, etc. of all polygonal contours of the building one by one, and then call different interpolation functions to obtain the plane coordinates of other measurement points according to the horizontal measurement point spacing. After all contours of the building are traversed, proceed to the next step.

[0075] Step 8: Arrange other measuring points from bottom to top based on the vertical measuring point spacing and the number of building floors and floor height information. After all floors of the building are traversed, generate the three-dimensional coordinates of all measuring points and proceed to the next step.

[0076] Step 9: Bind the generated three-dimensional coordinates of the measuring point with the attribute information of the building.

[0077] Step 10: Determine whether the traversal of all buildings has been completed. If not, continue to execute steps 7 to 9; if completed, proceed to the next step.

[0078] Step 11: Automatically output the three-dimensional coordinates of all measuring points and save them locally.

[0079] At this point, all steps for batch layout of surface measurement points of large-scale urban building complexes and automatic generation of coordinates are completed.

[0080] To further verify the effectiveness of the method of this embodiment, a typical application case is used as an example for explanation below.

[0081] Combination Figure 3 As shown in Table 1 and Table 2, this case includes 8 multi-story buildings, and the building plan outlines are quadrilateral and pentagonal, and the overall layout is an inverted triangle. Figure 3 A comparison diagram of the building GIS data before and after the projection coordinate transformation is given. It can be seen that after the original building GIS data passes through the functional modules such as the building GIS data parsing module, the projection transformation module and the building GIS data creation module established in this embodiment, the projection transformation of all building polygon outlines is automatically completed, and the coordinate origin is determined at the same time, laying the foundation for the subsequent measurement point layout.

[0082] Table 1 Vertex coordinates of building outlines before projection transformation

[0083]

[0084] Table 2 Vertex coordinates of building outlines after projection transformation

[0085]

[0086]

[0087] like Figure 4 As shown, Figure 4 (a) is the regenerated Shapefile data format file. Figure 4 (b) is a data file containing all the measurement point layout information, corresponding to the measurement point spacing of 1m and 2.5m. The file contains the measurement point number, the corresponding building number, the coordinate value in the three-dimensional rectangular coordinate system, and other auxiliary information.

[0088] To further verify whether the measurement points generated in this embodiment are strictly attached to the surrounding surfaces of the building complex, ParaView software is used to simultaneously open the newly created building GIS data file and the measurement point layout information data file. The measurement point layout effects of the two working conditions are shown in the following figure. Figure 5 As shown in the figure, whether the distance between measuring points is 1.0m or the distance between measuring points is increased to 2.5m, the method of this embodiment can automatically and batch generate regularly and evenly arranged measuring points on the surface around the building complex. All measuring points are attached to the outer surface of the building and will not be separated from or invaded into the interior of the building model.

[0089] Example 2: 2.5km in a city 2 Measurement point layout of buildings in the area

[0090] This embodiment is basically the same as the first embodiment, except that the example provided is a city with a radius of 2.5 km. 2 Measurement point layout of buildings in the area. Figure 6 As shown in the figure, a 2.5km×2.5km area in a city is selected as the measurement point layout area. The area contains a total of 1759 buildings with the number of floors ranging from 1 to 70.

[0091] Using the method of this embodiment, batch layout of measuring points is performed at intervals of 6m in the horizontal direction and 3m in the vertical direction. The obtained measuring point layout information data file is as follows: Figure 7 As shown, it contains more than 193,000 measurement points. It is particularly emphasized that the embodiment of this case was completed on a Lenovo computer with an 8-core Intel Core i7-9700 processor and 16GB RAM, and the total time consumed was less than 2 minutes.

[0092] Use ParaView software to view the layout of measurement points in the 3D urban building model, such as Figure 8As shown. Although the geometric contours of the buildings are complex and changeable, the method of this embodiment successfully generates regularly and evenly distributed measuring points on the surface around the building complex. The three-dimensional coordinates of the measuring points are further provided to the CFD program to calculate the explosion flow field transmission environment. After the calculation is completed, the explosion load data of 4 measuring points are randomly selected from more than 193,000 measuring points, as shown in Figures 9-11 As shown in the figure, the explosion load data shows that the measuring points are strictly attached to the surface of the building. After inspection, other measuring points all have pressure data, and the success rate of measuring point layout is 100%.

[0093] The above is only a specific embodiment of the present invention, but the protection scope of the present invention is not limited thereto. Any changes or substitutions within the technical scope disclosed by the present invention should be included in the protection scope of the present invention. Therefore, the protection scope of the present invention should be based on the protection scope of the claims.

Claims

1. A method for automatically generating three-dimensional coordinates of surface measurement points of large-scale urban buildings, characterized in that: The following steps are involved: Step 1: Obtain original building GIS data, and select valid building GIS data from the original building GIS data; The effective building GIS data refers to the original building GIS data including the two-dimensional plane geometric outline, location coordinates, number of building floors and building height of the target city building; Step 2: Get the projection coordinate system code; Step 3: For each building in the valid building GIS data selected in step 1, the building attribute information and building outline data are analyzed in turn; the building attribute information includes the number of floors and floor height information of each building; Step 4: Based on the building outline data parsed in step 3 and the projection coordinate system code obtained in step 2, the projection coordinate system is transformed to obtain the position coordinate data of the building outline data in the Cartesian coordinate system; Step 5: Select the minimum value of all vertex coordinates in the projection coordinate system as the coordinate origin of the Cartesian coordinate system, or directly reset the origin coordinates of the Cartesian coordinate system in step 4 to zero, and synchronously move the position coordinate data of all building outlines obtained in step 4; Step 6: for each building in the position coordinate data after the move in step 5, create building outline data and bind building attribute information in turn; Step 7: Select one of the buildings, calculate the side length, orientation and midpoint coordinates of all the contours of the building one by one, preset the horizontal measurement point spacing, and use the linear interpolation function to calculate the plane coordinates of all the measurement points on the building contour of the selected building until all the contours of the building are traversed; Step 8: Preset the vertical distance between measuring points, combine the number of floors and floor height information of the building selected in step 7, calculate the vertical coordinates of all measuring points from bottom to top, until all floors of the building are traversed, and obtain the three-dimensional coordinates of all measuring points; Step 9, return to step 7, and select another building until all buildings are traversed to obtain the three-dimensional coordinates of all building surface measurement points, and complete the automatic generation of three-dimensional coordinates of surface measurement points of large-scale urban building complexes.

2. The method for automatically generating three-dimensional coordinates of surface measurement points of large-scale urban buildings according to claim 1 is characterized in that: In step 1, specifically, original building GIS data is obtained from the geophysical information system data of the target city, and the building GIS data is two-dimensional vector data in shapefile format.

3. The method for automatically generating three-dimensional coordinates of surface measurement points of large-scale urban building complexes according to claim 1 or 2, characterized in that: In step 3, the building outline data of each building is analyzed in turn as follows: The latitude and longitude coordinate values ​​of the plane outline vertices, the area of ​​the closed polygon and the geometric center of each building in the valid building GIS data selected in step 1 are analyzed in turn and numbered; step 7 is specifically as follows: Select the first building in the order of numbers, calculate the side length, orientation and midpoint coordinates of all the contours of the building one by one, preset the horizontal measurement point spacing, and use the linear interpolation function to calculate the plane coordinates of all the measurement points on the building contour of the selected building until all the contours of the building are traversed; Step 9 is as follows: Return to step 7, select another building in the order of the numbers, until all buildings are traversed, and the three-dimensional coordinates of all building surface measurement points are obtained, thus completing the automatic generation of three-dimensional coordinates of surface measurement points of large-scale urban building complexes.

4. The method for automatically generating three-dimensional coordinates of surface measurement points of large-scale urban buildings according to claim 3 is characterized in that: In step 6, the building outline data is created specifically as follows: Regenerate the building GIS data in the plane rectangular coordinate system in the order of point-line-surface-attribute information.

5. The method for automatically generating three-dimensional coordinates of surface measurement points of large-scale urban buildings according to claim 4 is characterized in that: Also includes step 10: Bind the three-dimensional coordinates of all measurement points obtained in step 9 with the attribute information of the corresponding buildings to obtain more complete data.

6. A system for automatically generating three-dimensional coordinates of surface measurement points of large-scale urban buildings, used to execute the method for automatically generating three-dimensional coordinates of surface measurement points of large-scale urban buildings according to any one of claims 1 to 5, characterized in that: It includes original building GIS data module, user master parameter module, building GIS data analysis module, projection transformation module, building GIS data creation module and measurement point batch generation module; The original building GIS data module is used to obtain the original building GIS data and select the valid building GIS data in the original building GIS data; the valid building GIS data refers to the original building GIS data including the two-dimensional plane geometric outline, location coordinates, number of building floors and building height of the target city building; The user master control parameter module is used to obtain the projection coordinate system code, the preset horizontal measurement point spacing, and the preset vertical measurement point spacing; The building GIS data analysis module is used to analyze the building attribute information and building outline data of each building in the acquired valid building GIS data in turn; The projection transformation module transforms the projection coordinate system based on the parsed building outline data and the projection coordinate system code to obtain the position coordinate data of the building outline data in the Cartesian coordinate system; and is used to select the minimum value of all vertex coordinates as the coordinate origin of the Cartesian coordinate system or directly return the origin coordinate of the Cartesian coordinate system to zero, and simultaneously move the position coordinates of all building outlines synchronously; The building GIS data creation module is used to sequentially create building outline data and bind building attribute information for each building in the moved position coordinate data; The building attribute information includes the number of floors and floor height information of each building; The measuring point batch generation module is used to calculate the side length, orientation and midpoint coordinates of all the contours of each building one by one, and use a linear interpolation function to calculate the plane coordinates of all the measuring points on the building contour of the selected building until all the contours of the building are traversed; it is also used to combine the number of floors and floor height information of the building, calculate the coordinates of all the measuring points in the vertical direction from bottom to top, until all floors of the building are traversed, and the three-dimensional coordinates of all the measuring points are obtained; and then the three-dimensional coordinates of the measuring points on the surface of each building are obtained.

7. The system for automatically generating three-dimensional coordinates of surface measurement points of large-scale urban buildings according to claim 6 is characterized in that: It also includes a data automatic output module; The automatic data output module is used to automatically save and output the results generated by the building GIS data creation module and the measurement point batch generation module.

Images