Linear engineering settlement risk assessment method and device based on expandable sliding grid

Abstract

The invention provides a linear engineering settlement risk assessment method and device based on an expandable sliding grid. The method comprises the following steps: dividing geographic space data into grids according to a preset rule; s2, grids including linear engineering are selected as target grids; s3, taking the target grid as a center point, and expanding surrounding adjacent N grids as an integral grid; s4, determining a starting point grid, judging whether N grids adjacent to the starting point grid contain linear projects or not, if yes, sliding the whole grid to the position containing the linear projects, taking the starting point grid as a target grid, and returning to the step S3 until the sliding path of the whole grid completely passes through all the linear projects; s5, performing level division on the to-be-analyzed map; and S6, carrying out comprehensive risk assessment on the grid. The calculation resource consumption can be reduced, and the settlement condition around the linear engineering can be quickly judged, so that the subsequent settlement evaluation analysis is facilitated.

Description

Technical Field

[0001] The present invention relates to the technical field of risk assessment, and in particular to a linear engineering settlement risk assessment method and device based on an expandable sliding grid. Background Art

[0002] InSAR (Interferometric synthetic aperture radar, InSAR) combines synthetic aperture radar imaging technology and interferometric measurement technology, and uses the sensor's system parameters and imaging geometry to accurately measure the three-dimensional spatial position and tiny changes of a point on the surface. InSAR monitoring data includes the longitude and latitude of the monitoring point, as well as its sedimentation rate data (unit: mm / year) and cumulative sedimentation data (unit: mm). These results are recorded in the form of points (longitude and latitude). The monitoring point spacing of different InSAR satellite systems is different. For example, the monitoring data of the commonly used Sentinel satellite is spaced at 30m.

[0003] InSAR monitoring data can reflect the settlement conditions of the entire area on a large scale. However, for general linear projects such as subways, underground pipelines, power grids, etc., the focus is on the settlement conditions around their lines. Other external conditions are not within their consideration scope. On the contrary, a large amount of data not around the lines will become noise data, affecting the system's data analysis around the pipelines.

[0004] The traditional method is to select the surrounding area of ​​the line, adopt the line topology routing method, and perform one-by-one matching analysis with the settlement results map. This consumes a lot of computing resources, and a large number of linear projects do not actually have a complete line topology. If this method is to be used, a lot of manual calibration work is required, which greatly consumes manpower and time costs. Summary of the invention

[0005] The present invention aims to provide a linear engineering settlement risk assessment method and device based on an expandable sliding grid that overcomes the above-mentioned problems or at least partially solves the above-mentioned problems.

[0006] To achieve the above object, the technical solution of the present invention is specifically implemented as follows: One aspect of the present invention provides a linear engineering settlement risk assessment method based on an expandable sliding grid, comprising: S1, dividing the geographic spatial data into grids according to preset rules, each grid representing an area for storing and processing geographic information; S2, matching the size of the grid with the coordinate value of the linear project, and selecting a grid including the linear project as a target grid; S3, taking the target grid as the center point, expanding N adjacent grids around it as the overall grid, where N is a preset value; S4, determining a starting grid, judging whether N grids adjacent to the starting grid include the linear project, and if so, sliding the entire grid to a position including the linear project, taking the starting grid as the target grid and returning to execute S3, until the sliding path of the entire grid completely passes through all the linear projects; S5, dividing the map to be analyzed into levels; S6, obtaining the cumulative settlement values ​​of all the INSAR points included in the overall grid, evaluating the cumulative settlement values ​​of the INSAR points according to density, evaluating the cumulative settlement values ​​of the INSAR points according to the maximum and minimum differences, and evaluating the cumulative settlement values ​​of the INSAR points according to the maximum settlement; S7, performing a comprehensive risk assessment of the grid.

[0007] Optionally, dividing the geographic spatial data into grids according to a preset rule includes: The geographic spatial data is divided into grids whose sizes are consistent with the pixel sizes of electronic map tiles at different levels.

[0008] Optionally, matching the size of the grid with the coordinate value of the linear project and selecting the grid including the linear project as the target grid comprises: It is determined whether the centerline coordinate value of the linear project is within the range of the four-corner coordinate values ​​of the grid. If so, it is determined that the grid contains the linear project and the grid is used as the target grid.

[0009] Optionally, N=8.

[0010] Optionally, dividing the map to be analyzed into levels includes: The INSAR cumulative settlement values ​​in the map to be analyzed are divided into 15 levels according to the natural break point classification method to obtain the settlement range of each level.

[0011] Optionally, the evaluating according to the density based on the cumulative settlement value of the INSAR point includes: Calculate the first arithmetic mean A of the cumulative settlement value of the INSAR point within the maximum 5-level interval, and calculate the first proportion ; Calculate the second arithmetic mean B of the cumulative settlement value of the INSAR point in the middle 5-level interval, and calculate the second proportion ; Calculate the third arithmetic mean C of the cumulative settlement value of the INSAR point within the minimum 5-level interval, and calculate the third proportion ; Calculate D A+ B+ C, where if D is in the maximum 3-level interval, the density assessment value PD=1, if D is in the secondary 3 intervals, PD=0.5, and if D is in the remaining intervals, PD=0.

[0012] Optionally, the evaluating according to the maximum and minimum differences of the accumulated settlement values ​​of the INSAR points includes: Calculate the difference between the maximum and minimum values , if the difference The size that spans the middle 3-level range is recorded as , then the difference evaluation value P =1; if the difference The size that spans the middle 2-level interval is recorded as , then the difference evaluation value P =0.5; if the difference The size that spans the middle 1-level interval range is recorded as , then the difference evaluation value P =0.

[0013] Optionally, the evaluating according to the maximum settlement amount based on the cumulative settlement value of the INSAR point includes: Calculate the cumulative settlement of all INSAR settlement points If the maximum settlement value belongs to the maximum level 3 interval, it is recorded as , then the maximum evaluation value P =1, if the maximum settlement value belongs to the sub-3rd level interval, it is recorded as , then the maximum evaluation value P =0.5, if the maximum settlement value belongs to other intervals, it is recorded as , then the maximum evaluation value P =0.

[0014] Optionally, performing a comprehensive risk assessment of the grid includes: The comprehensive risk assessment value of the grid is calculated according to the following formula: calculate: +P

[0015] Wherein, if the comprehensive risk assessment value If it is greater than or equal to 2, it is a high risk, as described in the comprehensive risk assessment value If it is greater than or equal to 1, it is medium risk, and the rest is low risk.

[0016] Another aspect of the present invention provides a linear engineering settlement risk assessment device based on an expandable sliding grid, which uses the above-mentioned linear engineering settlement risk assessment method based on an expandable sliding grid to perform risk assessment.

[0017] It can be seen that the linear engineering settlement risk assessment method and device based on the expandable sliding grid provided by the present invention can reduce the consumption of computing resources and can quickly determine the settlement conditions around the linear engineering, so as to facilitate subsequent settlement assessment and analysis. BRIEF DESCRIPTION OF THE DRAWINGS

[0018] In order to more clearly illustrate the technical solutions of the embodiments of the present invention, the accompanying drawings required for use in the description of the embodiments will be briefly introduced below. Obviously, the accompanying drawings described below are only some embodiments of the present invention. For ordinary technicians in this field, other accompanying drawings can be obtained based on these accompanying drawings without paying creative work.

[0019] Figure 1 A flow chart of a linear engineering settlement risk assessment method based on an expandable sliding grid provided in an embodiment of the present invention; Figure 2 A schematic diagram of a grid sliding method provided in an embodiment of the present invention; Figure 3 A schematic diagram of a breakpoint classification method provided in an embodiment of the present invention. DETAILED DESCRIPTION

[0020] The exemplary embodiments of the present disclosure will be described in more detail below with reference to the accompanying drawings. Although the exemplary embodiments of the present disclosure are shown in the accompanying drawings, it should be understood that the present disclosure can be implemented in various forms and should not be limited by the embodiments set forth herein. On the contrary, these embodiments are provided to enable a more thorough understanding of the present disclosure and to fully convey the scope of the present disclosure to those skilled in the art.

[0021] Figure 1 A flow chart of a linear engineering settlement risk assessment method based on an expandable sliding grid provided by an embodiment of the present invention is shown. Figure 1 The linear engineering settlement risk assessment method based on an expandable sliding grid provided by an embodiment of the present invention includes: S1, divides the geospatial data into grids according to preset rules, each grid represents an area, and is used to store and process geographic information.

[0022] As an optional implementation of the embodiment of the present invention, dividing the geographic spatial data into grids according to a preset rule includes: dividing the geographic spatial data into grids whose sizes are consistent with the pixel sizes of electronic map tiles at different levels.

[0023] Specifically, the present invention first performs gridding, which means dividing geographic spatial data into grids according to certain rules, each grid representing an area, for storing and processing geographic information. The specific method of grid division is as follows: the grid size division in the present invention is consistent with the pixel size of electronic map tiles at different levels.

[0024] At present, electronic maps such as Google, Bing, AutoNavi, and Tiandi Map all use Google Maps’ tile pixel division method. The coordinate origin of Google Tile Maps’ standard is 180° east longitude and 85.05° north latitude, with the x-axis pointing right and the y-axis pointing downward. The maximum coverage range is from -180° to +180° longitude and -85.05° to +85.05° latitude.

[0025] The number of tiles contained in the electronic map varies according to the zoom level (from 0 to 21, a total of 22 levels). At the smallest map level 0, the world map consists of only one tile. The higher the tile level, the more tiles the world map consists of. At the same time, the number of tiles in the next level of the map is 4 times that of the previous level. Regardless of the level, the pixels in each tile are always 256*256.

[0026] S2, matching the size of the grid with the coordinate value of the linear project, and selecting the grid including the linear project as the target grid.

[0027] Specifically, the linear engineering of the present invention refers to linear engineering such as gas pipelines, subway lines, highways, etc. The coordinate value of the linear engineering refers to the center line coordinate of the linear engineering, including longitude and latitude, which is used to display the location of the pipeline in the geographic information system.

[0028] As an optional implementation of an embodiment of the present invention, matching the size of the grid with the coordinate value of the linear project, selecting the grid including the linear project as the target grid includes: determining whether the centerline coordinate value of the linear project is within the range of the four corner coordinate values ​​of the grid; if so, determining that the grid contains the linear project, and taking the grid as the target grid.

[0029] Specifically, each grid has the latitude and longitude coordinate values ​​of its center point, as well as the four corner coordinate values ​​of the square grid in which it is located, and also has the centerline coordinate value of the linear project. In this way, based on the size relationship between them, it can be determined whether the centerline coordinate value of the linear project is included in the four corner coordinate range of the grid.

[0030] S3, taking the target grid as the center point, expanding the surrounding N adjacent grids as the overall grid, where N is a preset value.

[0031] As an optional implementation of the embodiment of the present invention, N=8.

[0032] Specifically, the present invention takes the selected grid as the center point and expands the surrounding 8 adjacent grids as a whole. The expansion of the present invention is to link the central grid with the surrounding grids for comprehensive analysis. If it is not expanded, only a single grid is analyzed, especially in the case of a map with a very high zoom level, the number of points contained in a grid is too small, which is not conducive to comprehensive analysis.

[0033] S4, determine the starting grid, determine whether the N grids adjacent to the starting grid contain linear projects, if so, slide the entire grid to the position containing the linear projects, and return the starting grid as the target grid to execute S3 until the sliding path of the entire grid completely passes through all linear projects.

[0034] Specifically, in the above-mentioned matched grids, the present invention can use the first grid in the upper left corner of the entire map as the starting point, and judge whether its 8 adjacent grids contain the corresponding linear project. If so, slide to this point and repeat the above process (at this time, the remaining 8 grids involved in the previous point need to be eliminated). By analogy, the grid sliding coverage of the entire map can be completed to evaluate the settlement data in each sliding grid. If multiple grids adjacent to the starting grid all contain linear projects, slide to each of them respectively. Linear projects must be distributed within a certain map range. The sliding coverage of the grid of the entire map refers to the reference to the linear project. Figure 2 In the sliding method shown, the linear project starts from a certain starting point and the sliding path ends when it passes through all linear projects.

[0035] S5, dividing the map to be analyzed into levels.

[0036] As an optional implementation of the embodiment of the present invention, the level division of the map to be analyzed includes: dividing the INSAR cumulative settlement value in the map to be analyzed into 15 levels according to the natural break point classification method to obtain the settlement range of each level.

[0037] Specifically, the present invention first divides the INSAR cumulative settlement value in the entire image into 15 levels according to the natural break point classification method, and obtains the settlement range of each level, as shown in the example Figure 3 As shown, after classification, each row represents a level, and the range column is the specific range of values ​​contained in each level, which is the value of the above-mentioned InSAR cumulative precipitation in the present invention.

[0038] The map to be analyzed refers to, for example, if the object to be analyzed is Beijing subway, then Beijing subway is included in the map range of Beijing city, and the entire map refers to the map of Beijing city.

[0039] The natural break classification method is a common data classification method in the geographic information industry. This method was created by Professor Geogre Frederick Jenks. This method optimizes data grouping by using an algorithm to minimize the difference in data values ​​within the same category and maximize the difference in data values ​​between categories. This classification method seeks to minimize the average deviation from the class mean while maximizing the deviation from the mean of other groups. This method reduces the variance within the class and maximizes the variance between classes. Also known as the goodness of fit (GVF), it is equal to subtracting the SDCM (the sum of the squared deviations of the class means) from the SDAM (the sum of the squared deviations of the array means).

[0040] S6, obtaining the cumulative settlement values ​​of all INSAR points included in the overall grid, evaluating the INSAR points according to the density, evaluating the INSAR points according to the maximum and minimum differences, and evaluating the INSAR points according to the maximum settlement.

[0041] Specifically, the present invention requires confirmation of the evaluation scope. Each evaluation scope is the central grid and the eight adjacent grids outside, that is, a nine-square grid. The grids including the linear engineering participate in the data evaluation in the nine-square grid.

[0042] The present invention can be evaluated from three dimensions, including: settlement data based on density statistics, maximum and minimum difference evaluation and maximum settlement statistics.

[0043] When sedimentation data is calculated based on density: As an optional implementation of the embodiment of the present invention, evaluating according to the density based on the accumulated settlement value of the INSAR point includes: Calculate the first arithmetic mean A of the cumulative settlement values ​​of the INSAR points within the maximum 5-level interval, and calculate the first proportion ; Calculate the second arithmetic mean B of the cumulative settlement value of the INSAR point in the middle 5-level interval, and calculate the second proportion ; Calculate the third arithmetic mean C of the cumulative settlement value of the INSAR point within the minimum 5-level interval, and calculate the third proportion ; Calculate D A+ B+ C, where if D is in the maximum 3-level interval, the density assessment value PD=1, if D is in the secondary 3 intervals, PD=0.5, and if D is in the remaining intervals, PD=0.

[0044] Specifically, the cumulative settlement of all INSAR points included in the above-mentioned grids participating in the calculation is counted, and the value within the maximum 5-level interval is recorded as A, and A is the arithmetic mean of the points with values ​​within the maximum 5-level interval, and its number accounts for the statistical ; The values ​​in the middle 5-level interval are recorded as B, and B is the arithmetic mean of the points in the middle 5-level interval, and its number percentage statistics ; The value within the minimum 5-level interval is recorded as C, and C is the arithmetic mean of the points within the minimum 5-level interval, and its number percentage statistics Calculate D A+ B+ C. If D is in the maximum 3-level interval, the density assessment value PD=1; if D is in the secondary 3 intervals, PD=0.5; if D is in the remaining intervals, PD=0.

[0045] When evaluating the maximum and minimum differences: As an optional implementation of the embodiment of the present invention, evaluating according to the maximum and minimum differences of the accumulated settlement values ​​of the INSAR points includes: Calculate the difference between the maximum and minimum values , if the difference The size that spans the middle 3-level range is recorded as , then the difference evaluation value P =1; if the difference The size that spans the middle 2-level interval is recorded as , then the difference evaluation value P =0.5; if the difference The size that spans the middle 1-level interval range is recorded as , then the difference evaluation value P =0.

[0046] Specifically, the cumulative settlement of all INSAR settlement points included in the above-mentioned grids is counted, and the difference between the maximum and minimum values ​​is calculated. The difference value that spans the range of the middle 3 levels is recorded as , then the difference evaluation value P =1; the difference value that spans the range of the middle 2 levels is recorded as , then the difference evaluation value P =0.5; the difference value that spans the range of the middle level 1 interval is recorded as , then the difference evaluation value P =0.

[0047] When calculating the maximum settlement amount: As an optional implementation of the embodiment of the present invention, evaluating according to the maximum settlement amount based on the cumulative settlement value of the INSAR point includes: Calculate the cumulative settlement of all INSAR settlement points If the maximum settlement value belongs to the maximum level 3 interval, it is recorded as , then the maximum evaluation value P =1, if the maximum settlement value belongs to the sub-3rd level interval, it is recorded as , then the maximum evaluation value P =0.5, if the maximum settlement value belongs to other intervals, it is recorded as , then the maximum evaluation value P =0.

[0048] Specifically, the cumulative settlement of all INSAR settlement points included in the above-mentioned grids is counted. , the maximum settlement value belongs to the maximum level 3 interval is recorded as , then the maximum evaluation value P =1 is 1, and the maximum settlement value belongs to the sub-3rd level interval is recorded as , then the maximum evaluation value P =0.5, the maximum settlement value belongs to other intervals is recorded as , then the maximum evaluation value P =0.

[0049] S7, conduct a comprehensive risk assessment of the grid.

[0050] As an optional implementation of the embodiment of the present invention, performing a comprehensive risk assessment of a grid includes: The comprehensive risk assessment value of the grid is calculated according to the following formula: calculate: +P

[0051] Among them, if the comprehensive risk assessment value If it is greater than or equal to 2, it is a high risk, such as the comprehensive risk assessment value If it is greater than or equal to 1, it is medium risk, and the rest is low risk.

[0052] It can be seen that the linear engineering settlement risk assessment method based on the expandable sliding grid provided by the embodiment of the present invention can reduce the consumption of computing resources and can quickly determine the settlement situation around the linear engineering, so as to facilitate subsequent settlement assessment and analysis.

[0053] The present invention also provides a linear engineering settlement risk assessment device based on an expandable sliding grid, and adopts the above-mentioned linear engineering settlement risk assessment method based on an expandable sliding grid to perform risk assessment.

[0054] It can be seen that the linear engineering settlement risk assessment device based on the expandable sliding grid provided by the embodiment of the present invention can reduce the consumption of computing resources and can quickly determine the settlement conditions around the linear engineering, so as to facilitate subsequent settlement assessment and analysis.

[0055] The above are only embodiments of the present application and are not intended to limit the present application. For those skilled in the art, the present application may have various changes and variations. Any modification, equivalent replacement, improvement, etc. made within the spirit and principle of the present application should be included within the scope of the claims of the present application.

Claims

1. A linear engineering settlement risk assessment method based on an expandable sliding grid, characterized in that: include: S1, dividing the geographic spatial data into grids according to preset rules, each grid representing an area for storing and processing geographic information; S2, matching the size of the grid with the coordinate value of the linear project, and selecting a grid including the linear project as a target grid; S3, taking the target grid as the center point, expanding N adjacent grids around it as the overall grid, where N is a preset value; S4, determining a starting grid, judging whether N grids adjacent to the starting grid include the linear project, and if so, sliding the entire grid to a position including the linear project, taking the starting grid as the target grid and returning to execute S3, until the sliding path of the entire grid completely passes through all the linear projects; S5, dividing the map to be analyzed into levels; S6, obtaining the cumulative settlement values ​​of all the INSAR points included in the overall grid, evaluating the cumulative settlement values ​​of the INSAR points according to density, evaluating the cumulative settlement values ​​of the INSAR points according to the maximum and minimum differences, and evaluating the cumulative settlement values ​​of the INSAR points according to the maximum settlement; S7, performing a comprehensive risk assessment of the grid.

2. The method according to claim 1, characterized in that Dividing the geographic spatial data into grids according to preset rules includes: The geographic spatial data is divided into grids whose sizes are consistent with the pixel sizes of electronic map tiles at different levels.

3. The method according to claim 1, characterized in that The step of matching the size of the grid with the coordinate value of the linear project and selecting the grid including the linear project as the target grid comprises: It is determined whether the centerline coordinate value of the linear project is within the range of the four-corner coordinate values ​​of the grid. If so, it is determined that the grid contains the linear project and the grid is used as the target grid.

4. The method according to claim 1, characterized in that Said N=8.

5. The method according to any one of claims 1 to 4, characterized in that: The level division of the map to be analyzed includes: The INSAR cumulative settlement values ​​in the map to be analyzed are divided into 15 levels according to the natural break point classification method to obtain the settlement range of each level.

6. The method according to claim 5, characterized in that The evaluation according to the density based on the accumulated settlement value of the INSAR point includes: Calculate the first arithmetic mean A of the cumulative settlement value of the INSAR point within the maximum 5-level interval, and calculate the first proportion ; Calculate the second arithmetic mean B of the cumulative settlement value of the INSAR point in the middle 5-level interval, and calculate the second proportion ; Calculate the third arithmetic mean C of the cumulative settlement value of the INSAR point within the minimum 5-level interval, and calculate the third proportion ; Calculate D A+ B+ C, where if D is in the maximum 3-level interval, the density assessment value PD=1, if D is in the secondary 3 intervals, PD=0.5, and if D is in the remaining intervals, PD=0.

7. The method according to claim 6, characterized in that The evaluation according to the maximum and minimum differences of the accumulated settlement values ​​of the INSAR points includes: Calculate the difference between the maximum and minimum values , if the difference The size that spans the middle 3-level range is recorded as , then the difference evaluation value P =1; if the difference The size that spans the middle 2-level interval is recorded as , then the difference evaluation value P =0.5; if the difference The size that spans the middle 1-level interval range is recorded as , then the difference evaluation value P =0.

8. The method according to claim 7, characterized in that The evaluation according to the maximum settlement amount based on the cumulative settlement value of the INSAR point includes: Calculate the cumulative settlement of all INSAR settlement points If the maximum settlement value belongs to the maximum level 3 interval, it is recorded as , then the maximum evaluation value P =1, if the maximum settlement value belongs to the sub-3rd level interval, it is recorded as , then the maximum evaluation value P =0.5, if the maximum settlement value belongs to other intervals, it is recorded as , then the maximum evaluation value P =0.

9. The method according to claim 8, characterized in that The comprehensive risk assessment of the grid includes: The comprehensive risk assessment value of the grid is calculated according to the following formula: calculate: +P Wherein, if the comprehensive risk assessment value If it is greater than or equal to 2, it is a high risk, as described in the comprehensive risk assessment value If it is greater than or equal to 1, it is medium risk, and the rest is low risk.

10. A linear engineering settlement risk assessment device based on an expandable sliding grid, characterized in that: The risk assessment is performed using the linear engineering settlement risk assessment method based on an expandable sliding grid as described in any one of claims 1 to 9 above.

Images