Environment assessment method for ecological restoration of mine

Abstract

The invention discloses an environment assessment method for ecological restoration of a mine, relates to the technical field of ecological assessment of the mine, solves the problem that an original environment assessment mode is not high in pertinence, and comprises the following steps: generating a mountain model of a corresponding mine through a high-definition image of the mine, and constructing a mountain model of the corresponding mine through a constructed base plane body. According to the method, related division is carried out on different areas in a mountain model, so that the whole mountain model is divided into a flat area or a slope area, vertical confirmation is carried out on the slope area again, and the slope area is divided into a high slope area or a low slope area based on a distance difference value between adjacent vertical lines, so that the whole mountain model is divided into different areas; according to the method, region division is carried out based on the overall characteristics of the mine mountain, follow-up dynamic restoration mountain evaluation is facilitated, abnormal related ratio values are recognized, and therefore abnormal regions are locked quickly and effectively, and follow-up related personnel can carry out related management in time conveniently.

Description

Technical Field

[0001] The invention relates to the technical field of mine ecological assessment, and in particular to an environmental assessment method for mine ecological restoration. Background Art

[0002] Mining usually causes serious damage to land, water, vegetation, etc. Land may collapse, soil impoverishment, etc. Water bodies may be polluted, affecting the quality of water resources in surrounding areas; vegetation may be cut down, and biodiversity may be reduced.

[0003] The application with publication number CN116030356A relates to the field of image data processing technology, and specifically to an environmental assessment method for mine ecological restoration, the method comprising: obtaining a histogram of the H component in an HSV image of the mine ecology to be restored; performing Gaussian mixture fitting on the target histogram; performing vegetation Gaussian model detection on a sub-Gaussian model set, and screening out sub-Gaussian models adjacent to the vegetation Gaussian model from the sub-Gaussian model set; performing attribution analysis and screening processing on each hue value between the target hue value and each reference hue value in the target histogram; performing boundary analysis processing on a set of overlapping hue values ​​corresponding to each reference hue value; determining a target vegetation coverage rate corresponding to the mine ecology to be restored; and evaluating and judging the environment of the mine ecology to be restored. The present invention improves the efficiency of evaluating the environment for mine ecological restoration by performing data processing on the target histogram, and is mainly used in environmental assessment of mine ecological restoration.

[0004] Generally, for the overall assessment of mine ecological restoration, the target vegetation coverage is analyzed by using the target histogram method to perform color screening and analysis, so as to conduct a comprehensive assessment. However, this assessment method is too comprehensive. First of all, there are several areas with different slopes in the mine. The target vegetation coverage rate of the flat area is definitely higher than that of the slope area. It is generally difficult for the target vegetation to exist in the related areas with larger slopes. Therefore, the original environmental assessment method is not very targeted, and the resulting assessment results can only be said to be relatively comprehensive, but not precise. Summary of the invention

[0005] In view of the deficiencies in the prior art, the present invention provides an environmental assessment method for mine ecological restoration, which solves the problem that the original environmental assessment method is not very targeted.

[0006] To achieve the above objectives, the present invention is implemented through the following technical solutions: an environmental assessment method for mine ecological restoration, comprising the following steps:

[0007] P1. Collect drone images of mines that need ecological restoration. After collecting the overall high-definition images of the mines, generate the original mountain model of the mines based on the image features inside the high-definition images;

[0008] P2. Perform a preliminary division on the relevant areas inside the constructed mountain model, give priority to determining the base surface of the mountain model, construct a set of base surface bodies with a specified height based on the base surface, and construct several sets of equidistant vertical lines in the base surface bodies. Based on the specific intersection results of the base surface bodies and the mountain model, divide the surface area of ​​the mountain model into slope areas or flat areas. The specific sub-steps are as follows:

[0009] P21. Based on the constructed mountain model, determine the base surface where the mountain model is located, and then determine a group of parallel planes parallel to the base surface, the distance between the parallel planes and the base surface is a specified value, and the specified value is a preset value. Based on the parallel planes determined by the base surface, connect the edge contours of the two planes with a vertical line to determine a group of base surface bodies, the upper surface of the base surface body is the parallel plane, and the lower surface is the base surface, and construct a number of groups of equidistant vertical lines in the base surface body, and the vertical lines are perpendicular to the upper and lower surfaces of the base surface body;

[0010] P22, the established base surface body starts from the bottom of the base surface of the mountain model and gradually moves upward, and records the intersection points between the inner vertical line of the base surface body and the outer surface of the mountain model during each group of movement, and marks the generated groups of intersection points on the mountain model;

[0011] P23. Randomly select a group of edge intersection points as the initial point, determine the horizontal angle between the initial point and the surrounding adjacent intersection points, construct a line segment between the initial point and the surrounding adjacent intersection points, and evaluate the angle J between this line segment and the horizontal plane. i , whose horizontal plane is parallel to the base plane, where i represents different line segments, if J i ≤Y1, mark this line segment as a flat segment, where Y1 is the preset value and the angle J i is the internal angle, which is the internal angle between the horizontal plane and the lower endpoint of the line segment. If J i >Y1, mark the corresponding line segment as a slope segment. After the horizontal angles between adjacent intersections around the initial point are confirmed, the horizontal angles between other adjacent intersections are confirmed one by one.

[0012] P24, mark the area between the consecutive flat sections as the flat area, mark the area between the consecutive slope sections as the slope area, if the slope area contains the flat area, identify the area ratio ZB of the flat area in the slope area, if ZB>20%, keep the original marked flat area unchanged, if ZB≤20%, divide the flat area contained in the slope area into the slope area, if the flat area contains the slope area, identify the area ratio ZB of the slope area in the entire flat area, if ZB>20%, keep the original marked slope area unchanged, if ZB≤20%, divide the slope area contained in the flat area into the flat area;

[0013] P3. For the determined slope area, a set of horizontal standard surfaces is determined above the mountain model. Several interference points are confirmed from the horizontal standard surfaces, and the interference points are vertically mapped to the slope area. Then, based on the change in the distance value between adjacent mapping lines, the determined slope area is divided into a high slope area or a low slope area. The specific method is as follows:

[0014] P31. Based on the slope area marked in the mountain model, a set of horizontal standard surfaces parallel to the base surface are constructed. The horizontal standard surfaces are located above the mountain model, and the area of ​​the horizontal standard surfaces is consistent with the area of ​​the base surface.

[0015] P32. Determine a number of interference points in the horizontal standard plane, and the distance between each interference point is consistent, so that the corresponding interference point is vertically mapped to the slope area to determine the mapping point and synchronously determine the mapping line between the interference point and the mapping point, and no mapping processing is performed on the flat area;

[0016] P33, determine the line length difference Cz between two adjacent groups of mapping lines, and Cz>0, if Cz>Y2, where Y2 is a preset value, the two groups of mapping points are marked as high slope points, otherwise, the two groups of mapping points are marked as low slope points;

[0017] P34, identify the specific area covered by several high slope points in the corresponding slope area and divide it into a high slope area, and divide the part of the slope area that does not belong to the high slope area into a low slope area;

[0018] P4. For the mines that have completed ecological restoration, drone image collection is performed again. Based on the different areas divided in the mountain model corresponding to the mine, it is synchronously divided again, and the images of each divided area are compared to determine the green proportion of the relevant area and identify the abnormal area. The specific method is as follows:

[0019] P41. Based on the different areas divided in the mountain model and the specific boundaries of the corresponding areas, the regional image is intercepted on the mine image after image acquisition to obtain high-definition images of each different area;

[0020] P42. Compare the high-definition image with the green spectrum, where the green spectrum is a preset color spectrum, divide the areas with the same color in the high-definition image into standard areas, identify the proportion of the standard area in the corresponding area and mark it as Zq, where q represents different areas, and Zq = the area of ​​the standard area ÷ the total area of ​​the corresponding area. The confirmed proportion value Zq of different areas is directly displayed through the display terminal;

[0021] P43. Analyze and process the different proportion values ​​Zq belonging to the same classification area to determine the abnormal area; the specific method for determining the abnormal area is:

[0022] P431, sorting several groups of different proportion values ​​Zq from small to large to determine the sorting sequence;

[0023] P432. Identify the difference CC between adjacent proportion values ​​Zq in the sorting sequence, and CC>0. If CC<Y3, where Y3 is a preset value, the corresponding proportion value Zq is divided into the same type of value. Otherwise, no division is performed. The proportion values ​​belonging to the same type of value in the sorting sequence are confirmed one by one, and the areas related to the proportion values ​​that do not belong to the same type are marked as abnormal areas and displayed directly.

[0024] The present invention provides an environmental assessment method for mine ecological restoration. Compared with the prior art, it has the following beneficial effects:

[0025] The present invention generates a mountain model of the corresponding mine from a high-definition image of the mine, and then divides different areas inside the mountain model through the constructed base surface, so that the entire mountain model is divided into a flat area or a slope area, and the slope area is vertically confirmed again, and the slope area is divided into a high slope area or a low slope area based on the distance difference between adjacent vertical lines, so that the entire mountain model is divided into different areas. Based on the overall characteristics of the mine mountain, the area division is carried out to facilitate the subsequent dynamic repair of the mountain assessment;

[0026] For the identified several different areas, analyze the proportion of each area after the repair is completed, and then based on the different proportions of different areas, perform numerical confirmation to identify the more abnormal related proportions, so as to quickly and effectively lock the abnormal areas, making it convenient for subsequent relevant personnel to carry out relevant management in a timely manner. BRIEF DESCRIPTION OF THE DRAWINGS

[0027] Figure 1 It is a schematic diagram of the process of the present invention;

[0028] Figure 2 It is a schematic diagram of the regional division of the mountain model of the present invention. DETAILED DESCRIPTION

[0029] 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.

[0030] See also Figure 1 , the present application provides an environmental assessment method for mine ecological restoration, comprising the following steps:

[0031] P1) Use drones to collect images of mines that need ecological restoration. After the overall high-definition image of the mine is collected, the original mountain model of the mine is generated based on the image features inside the high-definition image. Specifically, since the method of generating the mountain model based on the collected high-definition image is relatively common in the prior art, it will not be described in detail here. By analyzing the stereo image pair, the digital elevation model (DEM) and point cloud data can be extracted; DEM data is the basis for mountain modeling, which contains the elevation information of the terrain. To ensure the accuracy of the image data, internal processing is required, including radiation correction, geometric correction, etc., to eliminate the imaging process. Various distortions, generate shadow maps with terrain features through DEM data, enhance the three-dimensional visual effect of the terrain, set parameters including DEM stretch type, color table, terrain features, azimuth, elevation, etc. to adjust the details and overall effect of the shadow map, and use the relevant functions in the three-dimensional modeling software or GIS software on the basis of DEM data and shadow maps to build a three-dimensional model of the mountain according to the extracted point cloud data, apply the processed high-definition image as a texture map to the three-dimensional model, improve the realism and detail performance of the model, adjust and optimize the local details of the model, including the contour of the mountain, slope, bumps and other terrain features, to make the model more accurate and vivid;

[0032] P2) Perform an initial division on the relevant areas inside the constructed mountain model, give priority to determining the base surface of the mountain model, construct a set of base surface bodies with a specified height based on the base surface, and construct several sets of equidistant vertical lines in the base surface bodies, and divide the surface area of ​​the mountain model into slope areas or flat areas based on the specific intersection results of the base surface bodies and the mountain model, wherein the specific sub-steps for performing the specific division are:

[0033] P21, Combination Figure 2 , based on the constructed mountain model, determine the base surface where the mountain model is located (which can be understood as the bottom surface of the mountain model), and then determine a group of parallel planes parallel to the base surface, the distance between the parallel planes and the base surface is a specified value (that is, height), and the specified value is a preset value. Based on the parallel planes determined by the base surface, the edge contours of the two planes are connected by vertical lines to determine a group of base surface bodies, the upper surface of the base surface body is a parallel plane, and the lower surface is a base surface. Construct a number of groups of equidistant vertical lines in the base surface body, and the vertical lines are perpendicular to the upper and lower surfaces of the base surface body. The specific number of the vertical lines is set in advance by the operator, and is generally not less than 1500 groups or even more;

[0034] P22, the established base surface body starts from the bottom of the base surface of the mountain model and gradually moves upward, and records the intersection points between the inner vertical line of the base surface body and the outer surface of the mountain model during each group of movement, and marks the generated groups of intersection points on the mountain model;

[0035] P23. Randomly select a group of edge intersection points as the initial point, determine the horizontal angle between the initial point and the surrounding adjacent intersection points, construct a line segment between the initial point and the surrounding adjacent intersection points, and evaluate the angle J between this line segment and the horizontal plane. i , whose horizontal plane is parallel to the base plane, where i represents different line segments, if J i ≤Y1, mark this line segment as a flat segment, where Y1 is a preset value, and its specific value is determined by the operator based on experience, generally 25°, and the angle J i is the internal angle between the horizontal plane and the lower end point of the line segment. If it is a flat surface, the angle is generally very small. If it is a slope surface, it generally does not exceed 90°. If there is a severely inclined slope surface, it may exceed 90°. If J i >Y1, mark the corresponding line segment as a slope segment. After the horizontal angles between adjacent intersections around the initial point are confirmed, the horizontal angles between other adjacent intersections are confirmed one by one.

[0036] P24, mark the area between the consecutive flat sections as the flat area, mark the area between the consecutive slope sections as the slope area, if the slope area contains the flat area, identify the area ratio ZB of the flat area in the slope area, if ZB>20%, keep the original marked flat area unchanged, if ZB≤20%, divide the flat area contained in the slope area into the slope area (that is, the flat area is changed and divided into the slope area);

[0037] Similarly, if the flat area contains a slope area, identify the area ratio ZB of the slope area in the entire flat area. If ZB>20%, keep the originally calibrated slope area unchanged. If ZB≤20%, divide the slope area contained in the flat area into the flat area.

[0038] Example: Assume that in the first base body (0-10 meters) at the bottom, we recorded multiple intersections and found that the angles between these points were all less than 25°, so this part was marked as a flat area. As the base body moves upward to 10-20 meters, the angles between the intersections begin to increase, and some exceed 25°, and these parts are marked as slope sections. If we find a small flat area (angle less than 25°) in a slope section, but the area of ​​this flat area accounts for only 15% of ZB, then according to the rules, this small flat area will be included in the slope area.

[0039] P3) For the determined slope area, a set of horizontal standard surfaces is determined above the mountain model, and a number of interference points are confirmed within the horizontal standard surfaces, so that the interference points are vertically mapped onto the slope area, and then based on the change in the distance value between adjacent mapping lines, the determined slope area is divided into a high slope area or a low slope area, wherein the specific method for determination is:

[0040] P31. Based on the slope area marked in the mountain model, a set of horizontal standard surfaces parallel to the base surface are constructed. The horizontal standard surfaces are located above the mountain model, and the area of ​​the horizontal standard surfaces is consistent with the area of ​​the base surface.

[0041] P32. Determine several interference points from the horizontal standard plane, and the distance between each interference point is consistent, so that the corresponding interference point is vertically mapped to the slope area to determine the mapping point (that is, the vertical line) and simultaneously determine the mapping line between the interference point and the mapping point. No mapping processing is performed for the flat area, and the number of interference points is prepared in advance by the operator, generally not less than 1000;

[0042] P33, determine the line length difference Cz between two adjacent groups of mapping lines, and Cz>0, if Cz>Y2, where Y2 is a preset value, and its specific value is determined by the operator based on experience, the two groups of mapping points are marked as high slope points, otherwise, the two groups of mapping points are marked as low slope points;

[0043] P34, identify the specific area covered by several high-slope points in the corresponding slope area (several high-slope points can be connected to determine the specific area covered, and the specific area is the area included between the connecting lines) and divide it into a high-slope area, and divide the part of the slope area that does not belong to the high-slope area into a low-slope area;

[0044] This processing method can also confirm the leveling area by adjusting specific related parameters. Similarly, the high slope area and the low slope area can also be confirmed by determining the corresponding angles in step P2.

[0045] For example, a slope area is marked in the mountain model, and the base area of ​​this area is 100 square meters. We will construct a horizontal standard surface with an area of ​​100 square meters above this area.

[0046] Determine the interference and mapping points:

[0047] On the horizontal standard plane, we evenly selected 25 interference points (5x5 grid), with equal distance between each interference point; then, we vertically mapped down from each interference point to the slope area surface and determined the corresponding mapping points. These mapping points are connected to the interference points by vertical lines to form mapping lines.

[0048] Calibrate high slope points and low slope points:

[0049] Next, we analyze the line length difference Cz between two adjacent groups of mapping lines. Assume that the preset value Y2 is 0.5 meters. If Cz is greater than 0.5 meters, we mark the two groups of mapping points as high-slope points; if Cz is less than or equal to 0.5 meters, the mapping points are marked as low-slope points.

[0050] Divide high slope area and low slope area:

[0051] Finally, we identify the area covered by all the high slope points in the slope area. By connecting these high slope points, we can determine the specific area they cover, which is classified as the high slope area. The remaining part, which does not belong to the high slope area, is classified as the low slope area.

[0052] P4) UAV image collection is conducted again for the mines that have completed ecological restoration. Based on the different areas divided in the mountain model corresponding to the mine, it is divided again synchronously, and the images of each divided area are compared to determine the green proportion of the relevant area and identify the abnormal area:

[0053] The specific method for determination is as follows:

[0054] P41. Based on the several different areas (flat area, low slope area and high slope area) divided in the mountain model, based on the specific edge lines of the corresponding areas, the regional image is intercepted in the mine image after the image acquisition is completed (the initial model of the mine image can be generated based on the overall edge contour of the image acquired again, and then based on the specific regional division process and characteristics, the picture of the corresponding area in the mine image is divided and confirmed, so as to confirm the regional image of each area) to obtain a high-definition image of each different area;

[0055] P42. Compare the high-definition image with the green spectrum, where the green spectrum is a preset chromatogram, which is prepared in advance by the operator and has several different green related colors set inside. Divide the areas with the same color in the high-definition image into standard areas, identify the proportion of the standard area in the corresponding area and mark it as Zq, where q represents different areas, and its Zq = standard area ÷ total area of ​​the corresponding area. The confirmed proportion value Zq of different areas is directly displayed through the display terminal;

[0056] P43. Analyze and process the different proportion values ​​Zq belonging to the same classification area (all belonging to the flat area, high slope area or low slope area) to determine the abnormal area. The specific method of determination is:

[0057] P431, sorting several groups of different proportion values ​​Zq from small to large to determine the sorting sequence;

[0058] P432, identify the difference CC of adjacent proportion values ​​Zq in the sorting sequence, and CC>0, if CC<Y3, where Y3 is a preset value, the specific value of which is determined by the operator based on experience, generally 10%, and classify the corresponding proportion value Zq into the same type of value, otherwise, no classification is performed, and the proportion values ​​belonging to the same type of value in the sorting sequence are confirmed one by one, and the area related to the proportion values ​​that do not belong to the same type of value is marked as an abnormal area, and directly displayed;

[0059] Example: Based on the preset green map, which contains a variety of green samples, it is used to evaluate the vegetation coverage or environmental compliance of the area; Image comparison: Compare the high-definition image of each area with the green map to identify the area in the image with the same color as the green map;

[0060] Calculate the proportion value Zq: For each area, calculate the proportion of the qualified area (the area whose color matches the green map) to the total area of ​​the area Zq. For example, if the total area of ​​the flattened area is 200 square meters, of which 150 square meters match the green map, then Zq is 75%.

[0061] Display percentage value: The Zq value of each area is directly displayed to the operator through the display terminal (such as the large screen or computer interface in the monitoring room).

[0062] P43: Analyze and process the proportion of different regions Zq

[0063] Assume that we have three flat areas with Zq values ​​of 60%, 70% and 80%.

[0064] Sorting: Sort these Zq values ​​from small to large and get the sequence: 60%, 70%, 80%.

[0065] Calculate the difference CC: Calculate the difference CC between adjacent Zq values. In this example, CC1 = 70% - 60% = 10%, CC2 = 80% - 70% = 10%.

[0066] Identifying similar values: If the preset value Y3 is 10%, then both CC values ​​are less than Y3, so these Zq values ​​are classified as similar values.

[0067] Identify abnormal areas: In this example, since all Zq values ​​are classified as the same value, there are no abnormal areas. If a CC value is greater than Y3, the corresponding Zq value will be considered abnormal, and the relevant area will be marked as an abnormal area and displayed.

[0068] Through this process, we can effectively identify and display the compliance status of each area and quickly identify areas where there may be problems. This helps the environmental monitoring and management of the mine.

[0069] Specifically, the reason for calibrating abnormal areas is that when the same area is undergoing the same degree of ecological restoration, some areas perform better, that is, the proportion is higher, but some areas perform poorly, that is, the proportion is too low. Then, areas with such low proportions are abnormal areas and need to be displayed to facilitate external management to understand.

[0070] Some of the data in the above formulas are dimensionless and numerically calculated. Meanwhile, the contents not described in detail in this specification belong to the prior art known to those skilled in the art.

[0071] The above embodiments are only used to illustrate the technical method of the present invention rather than to limit it. Although the present invention has been described in detail with reference to the preferred embodiments, those skilled in the art should understand that the technical method of the present invention may be modified or replaced by equivalents without departing from the spirit and scope of the technical method of the present invention.

Claims

1. An environmental assessment method for mine ecological restoration, characterized in that: The following steps are involved: P1. Collect drone images of mines that need ecological restoration. After collecting the overall high-definition images of the mines, generate the original mountain model of the mines based on the image features inside the high-definition images; P2. Perform a preliminary division of the relevant areas inside the constructed mountain model, give priority to determining the base surface of the mountain model, construct a set of base surface bodies with a specified height based on the base surface, and construct several sets of equidistant vertical lines in the base surface bodies. Based on the specific intersection results of the base surface bodies and the mountain model, divide the surface area of ​​the mountain model into a slope area or a flat area; P3. For the determined slope area, a set of horizontal standard surfaces is determined above the mountain model, and a number of interference points are confirmed from the horizontal standard surfaces, so that the interference points are vertically mapped to the slope area, and then the determined slope area is divided into a high slope area or a low slope area based on the change in the distance value between adjacent mapping lines; P4. UAV image collection is conducted again on the mines that have completed ecological restoration. Based on the several different areas divided in the mountain model corresponding to this mine, it is divided again synchronously, and the images of each divided area are compared to determine the green proportion of the relevant area and identify abnormal areas.

2. The environmental assessment method for mine ecological restoration according to claim 1 is characterized in that: In step P2, the specific sub-steps of dividing the surface area of ​​the mountain model into slope areas or flat areas are: P21. Based on the constructed mountain model, determine the base surface where the mountain model is located, and then determine a group of parallel planes parallel to the base surface, the distance between the parallel planes and the base surface is a specified value, and the specified value is a preset value. Based on the parallel planes determined by the base surface, connect the edge contours of the two planes with a vertical line to determine a group of base surface bodies, the upper surface of the base surface body is the parallel plane, and the lower surface is the base surface, and construct a number of groups of equidistant vertical lines in the base surface body, and the vertical lines are perpendicular to the upper and lower surfaces of the base surface body; P22, the established base surface body starts from the bottom of the base surface of the mountain model and gradually moves upward, and records the intersection points between the inner vertical line of the base surface body and the outer surface of the mountain model during each group of movement, and marks the generated groups of intersection points on the mountain model; P23. Randomly select a group of edge intersection points as the initial point, determine the horizontal angle between the initial point and the surrounding adjacent intersection points, construct a line segment between the initial point and the surrounding adjacent intersection points, and evaluate the angle J between this line segment and the horizontal plane. i , whose horizontal plane is parallel to the base plane, where i represents different line segments, if J i ≤Y1, mark this line segment as a flat segment, where Y1 is the preset value and the angle J i is the internal angle, which is the internal angle between the horizontal plane and the lower endpoint of the line segment. If J i >Y1, mark the corresponding line segment as a slope segment. After the horizontal angles between adjacent intersections around the initial point are confirmed, the horizontal angles between other adjacent intersections are confirmed one by one. P24. Mark the area between consecutive flat sections as a flat area, and mark the area between consecutive slope sections as a slope area. If the slope area contains a flat area, identify the area ratio ZB of the flat area in the slope area. If ZB>20%, keep the original marked flat area unchanged. If ZB≤20%, divide the flat area contained in this slope area into this slope area.

3. The environmental assessment method for mine ecological restoration according to claim 2 is characterized in that: In P24: If the flat area contains a slope area, identify the area ratio ZB of the slope area in the entire flat area. If ZB>20%, keep the originally calibrated slope area unchanged. If ZB≤20%, divide the slope area contained in the flat area into this flat area.

4. The environmental assessment method for mine ecological restoration according to claim 1 is characterized in that: In P3, the specific method of dividing the slope area into a high slope area or a low slope area is: P31. Based on the slope area marked in the mountain model, a set of horizontal standard surfaces parallel to the base surface are constructed. The horizontal standard surfaces are located above the mountain model, and the area of ​​the horizontal standard surfaces is consistent with the area of ​​the base surface. P32. Determine a number of interference points in the horizontal standard plane, and the distance between each interference point is consistent, so that the corresponding interference point is vertically mapped to the slope area to determine the mapping point and synchronously determine the mapping line between the interference point and the mapping point, and no mapping processing is performed on the flat area; P33, determine the line length difference Cz between two adjacent groups of mapping lines, and Cz>0, if Cz>Y2, where Y2 is a preset value, the two groups of mapping points are marked as high slope points, otherwise, the two groups of mapping points are marked as low slope points; P34. Identify the specific area covered by several high slope points in the corresponding slope area and divide it into the high slope area, and divide the part of the slope area that does not belong to the high slope area into the low slope area.

5. The environmental assessment method for mine ecological restoration according to claim 1 is characterized in that: In P4, the specific method for determining the green proportion values ​​of different areas is: P41. Based on the different areas divided in the mountain model and the specific boundaries of the corresponding areas, the regional image is intercepted on the mine image after image acquisition to obtain high-definition images of each different area; P42. Compare the high-definition image with the green spectrum, where the green spectrum is a preset color spectrum, divide the areas with the same color in the high-definition image into standard areas, identify the proportion of the standard area in the corresponding area and mark it as Zq, where q represents different areas, and Zq = the area of ​​the standard area ÷ the total area of ​​the corresponding area. The confirmed proportion value Zq of different areas is directly displayed through the display terminal; P43. Analyze and process the different proportion values ​​Zq belonging to the same classification area to determine the abnormal area.

6. The environmental assessment method for mine ecological restoration according to claim 5 is characterized in that: In P43, the specific method of determining the abnormal area is: P431, sorting several groups of different proportion values ​​Zq from small to large to determine the sorting sequence; P432. Identify the difference CC between adjacent proportion values ​​Zq in the sorting sequence, and CC>0. If CC<Y3, where Y3 is a preset value, the corresponding proportion value Zq is divided into the same type of value. Otherwise, no division is performed. The proportion values ​​belonging to the same type of value in the sorting sequence are confirmed one by one, and the areas related to the proportion values ​​that do not belong to the same type are marked as abnormal areas and displayed directly.

Images