Method, device and apparatus for quantitatively characterizing the location of geological hazards in slopes

Abstract

This application belongs to the field of geological disasters and environment, and particularly relates to a method, device, and equipment for quantitatively characterizing the development location of geological disasters in slopes. The method of this application includes: obtaining the slope boundary of the target slope and constructing a triangle that matches the slope boundary; obtaining the slope terrain point data within the triangle based on the triangle; for any slope terrain point D i , respectively adding the preset attribute information corresponding to D i ; wherein, the preset attribute information at least includes an angle attribute and a distance ratio attribute; screening the slope terrain point data to obtain disaster point data, and based on the angle attribute and distance ratio attribute of the disaster point data, plotting and displaying the disaster points within a preset area in the polar coordinate system. This solution can be used to solve the problem that the lateral position distribution law of disasters cannot be analyzed in the current susceptibility analysis of geological disasters, so as to carry out engineering construction such as houses and roads on slopes more scientifically.

Description

Technical Field

[0001] The present application belongs to the field of geological disasters and environment, and in particular, relates to a method, device and electronic equipment for quantitatively characterizing the development location of geological disasters on a slope in a high mountain canyon area. Background Art

[0002] The occurrence of geological disasters on slopes in high mountain canyon areas is closely related to the spatial position of the geological disasters in the slopes. By finding out the distribution of disasters on the slopes, revealing the laws of disaster development on the slopes, and finding the disaster-prone parts of the slopes, the threat of geological disasters on slopes to structures can be effectively reduced by avoiding or reinforcing them in advance when building houses or roads. This is crucial to ensuring the safe construction of projects and the safety of people's lives and property. Current research mainly focuses on the longitudinal position of the disaster in the slope, such as the elevation of the disaster, the tug-of-war height, and the relative height in the slope, but does not focus on the lateral position of the disaster in the slope.

[0003] Based on this, the present application proposes a method, device and electronic equipment for quantitatively characterizing the development location of geological hazards in a slope. Summary of the invention

[0004] The embodiments of the present application provide a method, device and electronic device for quantitatively characterizing the development location of geological disasters in slopes, which can be used to solve the problem that the lateral distribution pattern of disasters in slopes cannot be analyzed in the current geological disaster susceptibility analysis, so as to carry out construction of houses, roads and other projects on slopes in a more scientific way.

[0005] The present application embodiment adopts the following technical solutions:

[0006] The present application provides a method for quantitatively describing the location of geological hazards in a slope, the method comprising:

[0007] Acquire a slope boundary of a target slope, and construct a triangle matching the slope boundary, wherein three vertices of the triangle are A, B, and C respectively;

[0008] Based on the triangle, obtaining slope terrain point data within the triangle range;

[0009] For any slope topographic point D i , respectively add and D i The corresponding preset attribute information includes at least an angle attribute and a distance ratio attribute. The angle attribute is based on D i The target angle value calculated from the angle formed by point A and point B, the distance ratio attribute is based on each D i The angle properties and D i The distances to point A are calculated as the distance ratio;

[0010] The slope terrain point data is screened to obtain disaster point data, and based on the angle attribute and distance ratio attribute of the disaster point data, the disaster point is drawn and displayed in a preset area in a polar coordinate system.

[0011] In a specific embodiment, the angle attribute is based on D i The target angle value calculated from the angle formed by point A and point B includes:

[0012] For any slope terrain point D i , calculate D i ∠BAD formed by A and B i The initial angle value of the target slope is A, and point B is located on the left side of point C.

[0013] The initial angle values ​​are converted into target angle values ​​α according to the pre-acquired conversion coefficients to obtain the slope terrain points D i The corresponding angle property.

[0014] In a specific implementation, the transformation coefficients are obtained as follows:

[0015] The angle value of ∠BAC of the triangle is calculated, and the angle value of ∠BAC is multiplied by a coefficient φ. If the result of multiplication by the coefficient φ is 90°, the coefficient φ is used as a pre-acquired transformation coefficient.

[0016] In a specific implementation, each initial angle value is converted into each target angle value according to the pre-acquired transformation coefficient, specifically including:

[0017] For any slope topographic point D i , multiply its initial angle value by the previously acquired transformation coefficient, and then reduce it by 45° to obtain D i The target angle value α.

[0018] In a specific implementation, the distance ratio attribute is based on each D i The distance ratios calculated from the distance data to point A include:

[0019] Based on each slope topographic point D i angle attribute, dividing the slope terrain point data into a plurality of angle classes according to a preset angle value interval;

[0020] For any slope terrain point D in each angle class i , calculate D i The distance L from A i , and determine the D in each angle class iThe maximum distance L from A max ;

[0021] For any slope terrain point D in each angle class i , the distance between it and A is L i The maximum distance L in its angle class max The ratio β is taken as the value of D i The corresponding distance ratio attribute.

[0022] In a specific implementation manner, the preset angle value interval is 1°, based on each slope terrain point D i According to the angle attribute of the slope terrain point data, the slope terrain point data can be divided into 90 angle classes according to the angle value interval of 1°.

[0023] In a specific implementation, the slope terrain point data is screened to obtain disaster point data, specifically including:

[0024] Obtaining the disaster range of the target slope; wherein the disaster range is the disaster range pre-circled within the triangular range;

[0025] Extract the disaster information of each slope terrain point data within the disaster range, and add the disaster information to the corresponding slope terrain point D i Disaster attributes;

[0026] Disaster point data are screened out from the slope terrain point data according to disaster attributes.

[0027] In a specific implementation, based on the angle attribute and the distance ratio attribute of the disaster point data, the disaster point is drawn and displayed in a preset area in a polar coordinate system, specifically including:

[0028] Pre-construct the sector in polar coordinates;

[0029] Using the target angle value in the angle attribute and the distance ratio value in the distance ratio attribute of the disaster point data as the polar angle and polar diameter of the disaster point in the polar coordinate system;

[0030] The disaster point is drawn and displayed in the constructed fan-shaped area according to the polar angle and polar diameter of the disaster point.

[0031] The present application also provides a device for quantitatively describing the location of geological hazards in a slope, the device comprising:

[0032] A slope boundary acquisition module, used for acquiring the slope boundary of the target slope and constructing a triangle matching the slope boundary, wherein the three vertices of the triangle are A, B and C respectively;

[0033] A terrain point data acquisition module, used for acquiring slope terrain point data within the triangle range based on the triangle;

[0034] The attribute acquisition module is used to obtain the attribute of any slope terrain point D i , respectively add and D i The corresponding preset attribute information includes at least an angle attribute and a distance ratio attribute; wherein the angle attribute is based on D i The target angle value calculated from the angle formed by point A and point B, the distance ratio attribute is based on each D i The angle properties and D i The distances to point A are calculated as the distance ratio;

[0035] A disaster point screening module is used to screen the slope terrain point data to obtain disaster point data;

[0036] The display module is used to draw and display the disaster point in a preset area in a polar coordinate system based on the angle attribute and distance ratio attribute of the disaster point data.

[0037] The present application also provides an electronic device, including:

[0038] processor;

[0039] The memory is used to store a computer program, and the computer program is executed by a processor to perform any of the methods described above.

[0040] At least one of the above technical solutions adopted in the embodiments of the present application can achieve the following beneficial effects: This solution generalizes the target slope in the three-dimensional space into a triangle in the plane, thereby obtaining the slope terrain point data within the triangle range, and then adds at least an angle attribute and a distance ratio attribute to each slope terrain point data, and then screens out the disaster point data, and uses the angle attribute and distance ratio attribute of the disaster point data to draw and display the disaster point in a preset area in the polar coordinate system. In other words, the present application uses the target angle value α of the angle formed by the geological disaster point and the vertex of the triangle to characterize the lateral position of the geological disaster point in the target slope, and uses the ratio β (i.e., L ) of the distance between the geological disaster point and the vertex of the triangle and the distance from the vertex to the bottom boundary in that direction. i / L max ), to describe the longitudinal position of the geological hazard point in the target slope, and finally to describe the shape of the target slope based on the transformation coefficient φ. Thus, the spatial position of the geological hazard in the target slope can be described by (φ, α, β). This scheme can be used to solve the problem that the horizontal position distribution law of the hazard in the slope cannot be analyzed in the current geological hazard susceptibility analysis, so as to carry out the construction of houses, roads and other projects on the slope more scientifically. BRIEF DESCRIPTION OF THE DRAWINGS

[0041] The drawings described herein are used to provide a further understanding of the present application and constitute a part of the present application. The illustrative embodiments of the present application and their descriptions are used to explain the present application and do not constitute an improper limitation on the present application. In the drawings:

[0042] Figure 1 It is a flow chart of a method for quantitatively describing the location of geological hazards in a slope in an embodiment of the present application;

[0043] Figure 2 is a schematic diagram of obtaining a target slope in a triangular range by generalizing the target slope based on a remote sensing image in an embodiment of the present application;

[0044] Figure 3 It is the display of the distribution of disaster points in the embodiment of this application;

[0045] Figure 4 It is a framework diagram of a specific embodiment of a method for quantitatively describing the location of geological hazards in a slope of the present application;

[0046] Figure 5 It is a schematic diagram of a device for quantitatively depicting the development location of geological hazards on a slope in an embodiment of the present application. DETAILED DESCRIPTION

[0047] The present application proposes a method, device and electronic device for quantitatively characterizing the location of geological hazards in a slope. In order to make the purpose, technical solution and advantages of the present application clearer, the technical solution of the present application will be clearly and completely described below in combination with the specific embodiments of the present application and the corresponding drawings. Obviously, the described embodiments are only part of the embodiments of the present application, not all of the embodiments. Based on the embodiments in the present application, all other embodiments obtained by ordinary technicians in this field without making creative work are within the scope of protection of this application.

[0048] First refer to Figure 1 , Figure 1 1 is a flow chart of a method for quantitatively describing the location of geological hazards in a slope in an embodiment of the present application. Figure 1 As shown, the method of the present application includes:

[0049] S110 , obtaining a slope boundary of the target slope, and constructing a triangle matching the slope boundary, wherein three vertices of the triangle are A, B, and C respectively.

[0050] For the target slope, the corresponding target slope image can be obtained based on the remote sensing image, and then the slope range can be manually circled in the target slope image. In this way, the slope boundary of the target slope can be directly obtained, that is, the shape of the target slope can be obtained. The target slope is further generalized into a triangle. Specifically, a triangle is constructed based on the shape of the slope so that the three sides of the triangle coincide with the slope boundary as much as possible, that is, a triangle matching the slope boundary is constructed, and the three vertices of the triangle are A, B, and C respectively. Figure 2 , Figure 2 Schematic diagram of obtaining a target slope in a triangular range by generalizing the target slope based on a remote sensing image in an embodiment of the present application. Figure 2 As shown in Figure 1, based on the remote sensing image, the boundary of the target slope is first circled ( Figure 2 The target slope is then divided into two groups, and then a triangle 2 is drawn based on the boundary range 1. The three sides of the triangle basically coincide with the boundary range 1, and a triangle matching the slope boundary is obtained. Among them, the vertex A of the triangle is the vertex of the target slope, and the vertex B of the triangle is located to the left of the vertex C. After obtaining the generalized triangle, the target slope within the triangle range can be retained.

[0051] S120: Acquire slope terrain point data within the triangle range based on the triangle.

[0052] In this step, the constructed triangle is used as the boundary of the target slope to obtain the "triangular state" of the target slope. As an example, based on the "triangular state" of the target slope, the digital elevation model of the slope can be clipped using ArcGIS software to obtain a "triangular state" slope terrain file named "slope terrain.tif", and then the raster conversion function of ArcGIS is used to convert the "slope terrain.tif" file into slope terrain point data named "slope point.shp". Those skilled in the art can also obtain the slope terrain point data of the "triangular state" target slope by other means, and this application does not specifically limit this.

[0053] S130, for any slope terrain point D i , respectively add and D i Corresponding preset attribute information; wherein the preset attribute information at least includes angle attribute and distance ratio attribute.

[0054] As an example, for any slope terrain point D i, you can add multiple attributes to it, these attributes can include X attribute, Y attribute, Z attribute, JD attribute, JD2 attribute, JDL attribute, JL attribute and JLB attribute. The X attribute is used to store the X-axis coordinate value, the Y attribute is used to store the Y-axis coordinate value, the Z attribute is used to store the Z-axis coordinate value (i.e., elevation), the JD attribute is used to store the initial angle value, the JD2 attribute is used to store the target angle value, the JDL attribute is used to store the angle class, the JL attribute is used to store the vertex distance, and the JLB attribute is used to store the distance ratio. Among them, for any slope terrain point D i The X-axis coordinate value, Y-axis coordinate value and Z-axis coordinate value of the target slope can be calculated by a suitable method, for example, the computational geometry function of ArcGIS can be used to directly obtain the X-axis coordinate value, Y-axis coordinate value and Z-axis coordinate value of the target slope. i The X-axis coordinate value and Y-axis coordinate value of the slope terrain point D are extracted by using the function of extracting multiple values ​​to points, that is, the elevation information contained in "slope terrain.tif", that is, the Z-axis coordinate value, and added to the corresponding slope terrain point D i The present embodiment does not specifically limit how to obtain the X attribute, the Y attribute, and the Z attribute.

[0055] In step S130, the angle attribute refers to the aforementioned JD2 attribute, which is based on D i The target angle value is calculated from the angle formed by point A and point B. As an example, when calculating the angle attribute, for any slope terrain point D i , first calculate D i ∠BAD formed by A and B i Then, each initial angle value α′ is converted into each target angle value α according to the pre-acquired transformation coefficient to obtain each slope terrain point D i Corresponding angle attribute. The initial angle value α′ is stored in the JD attribute, and the target angle value α is stored in the angle attribute (ie, JD2 attribute). Among them, the pre-acquired transformation coefficient is obtained as follows: calculate the angle value of the triangle ∠BAC, multiply the angle value of ∠BAC by the coefficient φ, and if the result of multiplying by the coefficient φ is 90°, then use the coefficient φ as the pre-acquired transformation coefficient. That is to say, before calculating the target angle value, first calculate the angle value of ∠BAC, and then multiply the angle value of ∠BAC by the transformation coefficient φ to make the angle of ∠BAC become 90°. Then calculate ∠BAD i The initial angle value α′, the target angle value α=φ*α′. Add α to D i The corresponding angle attribute (ie JD2 attribute).

[0056] In an optional embodiment, each initial angle value is converted into each target angle value according to the pre-acquired conversion coefficient, specifically including: for any slope terrain point Di , multiply its initial angle value by the previously acquired transformation coefficient, and then reduce it by 45° to obtain D i The target angle value α. That is, in this embodiment, the target angle value α = φ * α′ - 45°. Then add this α to the value corresponding to D i The corresponding angle attribute (i.e., JD2 attribute) can be designed in this way to improve the display effect of the disaster point in the final result, which will be explained later.

[0057] In step S130, the distance ratio attribute refers to the aforementioned JLB attribute, which is based on the D i The distance ratio is calculated from the distance data to point A. As an example, when calculating the distance ratio, the following steps are included:

[0058] S31, based on each slope topographic point D i The slope terrain point data is divided into multiple angle classes according to the preset angle value interval.

[0059] In S31, the preset angle value interval can be set to 1°, that is, each slope terrain point D i In the angle attribute of , the points with the target angle value of [-45, -44) are classified as the first category, the points with the target angle value of [-44, -43) are classified as the second category, the points with the target angle value of [-43, -42) are classified as the third category, and so on. A total of 90 angle categories can be obtained, and each slope terrain point D i The corresponding classification values ​​are stored in the aforementioned JDL attributes.

[0060] In this example, the classification is carried out according to [-45,-44), [-44,-43), [-43,-42)..., and the corresponding target angle value α is calculated according to the formula α=φ*α′-45°.

[0061] S32, for any slope terrain point D in each angle class i , calculate D i The distance L from A i , and determine the D in each angle class i The maximum distance L from A max .

[0062] In S32, as an example, according to point D i The X-axis coordinate value and Y-axis coordinate value of each slope terrain point D in each angle class are calculated using the Euclidean distance formula. i The distance L to the triangle vertex A i , the distance L can be i Save in the above and D i In the corresponding JL attribute. For each point D in the angle classi , we can determine the angle class D i The maximum distance L from A max This embodiment provides a method for calculating D i The distance to the triangle vertex A is not specifically limited, and a suitable algorithm can be selected according to actual needs. The Euclidean distance formula will not be described in detail.

[0063] S33, for any slope terrain point D in each angle class i , the distance between it and A is L i The maximum distance L in its angle class max The ratio β is taken as the value of D i The corresponding distance ratio attribute.

[0064] In S33, as an example, after calculating any slope terrain point D i The distance L to A i Afterwards, according to any slope topographic point D i The angle class in which it is located, calculate the distance ratio β=L i / L max Among them, L i is the current slope terrain point D i The distance to point A, L max This is the current slope topographic point D i Each point in the angle class D i The maximum distance value among the distances from A to A. The calculated distance ratio β is saved in the distance ratio attribute (ie, JLB attribute).

[0065] In summary, this embodiment is mainly for any slope terrain point D i , respectively add and D i The corresponding angle attribute and distance ratio attribute. In the above example, for any slope terrain point D i Adding multiple attributes makes it easier to calculate D i The angle attribute and distance ratio attribute of .

[0066] S140, screening the slope terrain point data to obtain disaster point data, and based on the angle attribute and distance ratio attribute of the disaster point data, drawing and displaying the disaster point in a preset area in the polar coordinate system.

[0067] In this step, the slope terrain point data is screened to obtain disaster point data, which may include:

[0068] S41. Obtain the disaster range of the target slope; wherein the disaster range is the disaster range pre-circled within the triangle range.

[0069] As an example, after obtaining the corresponding target slope image based on the remote sensing image, the disaster range in the target slope image can be circled. Figure 2 , Figure 2 The middle circle depicts the range of the target slope. Similarly, the disaster range within the target slope can be circled. Based on the remote sensing image, the disaster range of the target slope can be clearly seen. One way is to circle it manually, which will not be described in detail here.

[0070] S42, extracting disaster information from the data of each slope terrain point within the disaster range, and adding the disaster information to the corresponding slope terrain point D i disaster attributes.

[0071] In step S42, for any slope terrain point D i In addition to the multiple attributes added above, you can also add i Add the ZH attribute (i.e., disaster attribute). As an example, you can use the ArcGIS surface-to-raster function to convert the "disaster range" obtained by S41 into a raster file named "disaster range.tif", and set the raster value to 1, and then extract the information of "disaster range.tif" to the corresponding slope terrain point D i It can be understood that for all slope terrain point data, the slope terrain point D with a ZH attribute value of 1 i The slope terrain point D is a disaster point, and the ZH attribute value is not 1 i It is a non-disaster point.

[0072] S43, filtering out disaster point data from the slope terrain point data according to disaster attributes.

[0073] In this way, the disaster point data can be screened out from the slope terrain point data through steps S41-S43.

[0074] Further, in step S140, after the disaster point is screened out, the disaster point is drawn and displayed in a preset area in the polar coordinate system based on the angle attribute and distance ratio attribute of the disaster point data, including:

[0075] A fan-shaped area is pre-constructed in the polar coordinate system; the target angle value in the angle attribute of the disaster point data and the distance ratio value in the distance ratio attribute are used as the polar angle and polar diameter of the disaster point in the polar coordinate system; the disaster point is drawn and displayed in the constructed fan-shaped area according to the polar angle and polar diameter of the disaster point.

[0076] As an example, you can use drawing software, select the polar coordinate system, and draw a sector with a vertex of 90° and a radius of 1 in two-dimensional space to represent the slope in space. Export the screened disaster point data and save it as a "disaster point.dbf" file. Use the α value under the angle attribute and the β value under the distance ratio attribute in the "disaster point.dbf" file as the polar angle and polar diameter of the disaster point in the polar coordinate system, input them into the drawing software in sequence, draw the disaster point in the sector area and display it. It should be noted here that when calculating the angle attribute, the initial angle value is multiplied by the transformation coefficient and then reduced by 45°. The final display effect is as follows: Figure 3 If the angle is not reduced by 45°, the displayed sector area may appear as Figure 3 The middle sector area is rotated 45° counterclockwise or clockwise.

[0077] Reference Figure 4 , Figure 4 This is a framework diagram of a specific embodiment of a method for quantitatively describing the location of geological hazards in a slope. Figure 4 As shown, the inventive concept of the present application mainly includes: slope data collection, including terrain data, remote sensing images, and disaster data. After collecting the corresponding data, the slope terrain is converted into points, including circling the slope range, generalizing the slope triangle, and obtaining the coordinates of the slope points. Then the lateral position of any point in the slope is calculated, including angle value calculation, angle value transformation, and slope point classification. Then the longitudinal position of any point in the slope is calculated, including vertex distance calculation, maximum distance calculation, and distance ratio determination. Finally, the position of the disaster point in the slope is obtained. This embodiment will not be described in detail here, and please refer to the detailed description above. Those skilled in the art will understand that the method of the present application does not limit the execution order of the specific steps, and the optimal execution method can be selected according to actual calculation needs.

[0078] The present application also provides a device for quantitatively characterizing the location of geological hazards in a slope, such as Figure 5 As shown, the device comprises:

[0079] A slope boundary acquisition module 10 is used to acquire the slope boundary of the target slope and construct a triangle matching the slope boundary, wherein the three vertices of the triangle are A, B and C respectively;

[0080] A terrain point data acquisition module 20, for acquiring slope terrain point data within the triangle range based on the triangle;

[0081] The attribute acquisition module 30 is used to obtain the attribute of any slope terrain point D i , respectively add and D i Corresponding preset attribute information; wherein the preset attribute information at least includes an angle attribute and a distance ratio attribute; wherein the angle attribute is based on D iThe target angle value calculated from the angle formed by point A and point B; the distance ratio attribute is based on each D i The angle properties and D i The distances to point A are calculated as the distance ratio;

[0082] A disaster point screening module 40 is used to screen the slope terrain point data to obtain disaster point data;

[0083] The display module 50 is used to draw and display the disaster point in a preset area in the polar coordinate system based on the angle attribute and distance ratio attribute of the disaster point data.

[0084] The specific implementation of the device is described in the above method, which will not be repeated here. It can be understood by those skilled in the art that the device of the present application does not limit the execution order of the specific steps, and the optimal execution mode can be selected according to the actual calculation needs.

[0085] The present application also provides an electronic device, including:

[0086] processor;

[0087] The memory is used to store a computer program, and the computer program is executed by the processor to perform the above method.

[0088] In summary, the present application proposes a method, device and apparatus for quantitatively characterizing the location of geological hazards in a slope. First, the target slope in three-dimensional space is generalized into a triangle in a plane, and then the target angle value α formed by the geological hazard point and the triangle vertex is used to characterize the lateral position of the geological hazard point in the target slope, and the ratio β (i.e., L) of the distance between the geological hazard point and the triangle vertex and the distance from the vertex to the bottom boundary in that direction is used to characterize the lateral position of the geological hazard point in the target slope. i / L max ), to characterize the longitudinal position of the geological hazard point in the target slope. Finally, the shape of the target slope is characterized based on the transformation coefficient φ. Thus, the spatial position of the geological hazard in the target slope can be characterized by (φ, α, β). This scheme can be used to solve the problem that the horizontal position distribution law of the hazard in the slope cannot be analyzed in the current geological hazard susceptibility analysis, so as to carry out the construction of houses, roads and other projects on the slope more scientifically.

[0089] The above describes specific embodiments of this specification, and other embodiments are within the scope of the appended claims. In some cases, the actions or steps described in the claims can be performed in an order different from that in the embodiments and still achieve the desired results. The same and similar parts between the various embodiments in this specification can be referenced to each other, and each embodiment focuses on the differences from other embodiments.

[0090] The above is only an embodiment of this specification and is 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 in the scope of the claims of the present application.

Claims

1. A method for quantitatively characterizing the location of geological hazards in a slope, characterized in that: The method includes: Acquire a slope boundary of a target slope, and construct a triangle matching the slope boundary, wherein three vertices of the triangle are A, B, and C respectively; Based on the triangle, obtaining slope terrain point data within the triangle range; For any slope topographic point D i , respectively add and D i The corresponding preset attribute information includes at least an angle attribute and a distance ratio attribute; wherein the angle attribute is based on ∠BAD i The target angle value is calculated from the initial angle value of the target slope, A is the vertex of the target slope; the distance ratio attribute is based on each D i The angle properties and D i The distances to point A are calculated as the distance ratio; Screening the slope terrain point data to obtain disaster point data, and drawing and displaying the disaster point in a preset area in a polar coordinate system based on the angle attribute and distance ratio attribute of the disaster point data; The distance ratio attribute is based on each D i The angle properties and D i The distances to point A are calculated as follows: Based on each slope topographic point D i angle attribute, dividing the slope terrain point data into a plurality of angle classes according to a preset angle value interval; For any slope terrain point D in each angle class i , calculate D i The distance Li between A and D is determined in each angle category. i The maximum distance L from A max ; For any slope terrain point D in each angle class i , the distance between it and A is L i The maximum distance L in its angle class max The ratio β is taken as the value of D i The corresponding distance ratio attribute.

2. The method according to claim 1, characterized in that The angle property is based on ∠BAD i The target angle value calculated from the initial angle value includes: For any slope terrain point D i , calculate D i ∠BAD formed by A and B i The initial angle value of the target slope is A, and point B is located on the left side of point C. The initial angle values ​​are converted into target angle values ​​α according to the pre-acquired conversion coefficients to obtain the slope terrain points D i The corresponding angle property.

3. The method according to claim 2, characterized in that The transform coefficients are obtained as follows: The angle value of ∠BAC of the triangle is calculated, and the angle value of ∠BAC is multiplied by a coefficient φ. If the result of multiplication by the coefficient φ is 90°, the coefficient φ is used as a pre-acquired transformation coefficient.

4. The method according to claim 3, characterized in that The initial angle values ​​are converted into target angle values ​​according to the pre-acquired transformation coefficients, specifically including: For any slope topographic point D i , multiply its initial angle value by the previously acquired transformation coefficient, and then reduce it by 45° to obtain D i The target angle value α.

5. The method according to claim 2, characterized in that: The preset angle value interval is 1°, based on each slope terrain point D i According to the angle attribute of the slope terrain point data, the slope terrain point data can be divided into 90 angle classes according to the angle value interval of 1°.

6. The method according to any one of claims 2 to 5, characterized in that The slope terrain point data is screened to obtain disaster point data, specifically including: Obtaining the disaster range of the target slope; wherein the disaster range is the disaster range pre-circled within the triangular range; Extract the disaster information of each slope terrain point data within the disaster range, and add the disaster information to the corresponding slope terrain point D i Disaster attributes; Disaster point data are screened out from the slope terrain point data according to disaster attributes.

7. The method according to claim 6, characterized in that Based on the angle attribute and the distance ratio attribute of the disaster point data, the disaster point is drawn and displayed in a preset area in the polar coordinate system, specifically including: Pre-construct the sector in polar coordinates; Using the target angle value in the angle attribute and the distance ratio value in the distance ratio attribute of the disaster point data as the polar angle and polar diameter of the disaster point in the polar coordinate system; The disaster point is drawn and displayed in the constructed fan-shaped area according to the polar angle and polar diameter of the disaster point.

8. A device for quantitatively describing the location of geological hazards in a slope, characterized in that: The device includes: A slope boundary acquisition module, used for acquiring the slope boundary of the target slope and constructing a triangle matching the slope boundary, wherein the three vertices of the triangle are A, B and C respectively; A terrain point data acquisition module, used for acquiring slope terrain point data within the triangle range based on the triangle; The attribute acquisition module is used to obtain the attribute of any slope terrain point D i , respectively add and D i The corresponding preset attribute information includes at least an angle attribute and a distance ratio attribute; wherein the angle attribute is based on ∠BAD i The target angle value is calculated from the initial angle value of the target slope, A is the vertex of the target slope; the distance ratio attribute is based on each D i The angle properties and D i The distances to point A are calculated as the distance ratio; A disaster point screening module is used to screen the slope terrain point data to obtain disaster point data; A display module, used for drawing and displaying the disaster point in a preset area in a polar coordinate system based on the angle attribute and distance ratio attribute of the disaster point data; The distance ratio attribute is based on each D i The angle properties and D i The distances to point A are calculated as follows: Based on each slope topographic point D i angle attribute, dividing the slope terrain point data into a plurality of angle classes according to a preset angle value interval; For any slope terrain point D in each angle class i , calculate D i The distance L from A i , and determine the D in each angle class i The maximum distance L from A max ; For any slope terrain point D in each angle class i , the distance between it and A is L i The maximum distance L in its angle class max The ratio β is taken as the value of D i The corresponding distance ratio attribute.

9. An electronic device, characterized in that: include: processor; A memory for storing a computer program, wherein the computer program is executed by a processor to perform the method according to any one of claims 1 to 7.

Images