A method and device for repairing vegetation damage in high-cold and arid coal mine areas

Abstract

The application provides a kind of method and device for high-cold arid coal mine area vegetation damage repair, it is related to mine vegetation technical field, the present application is divided into sub-region to coal mine area, the water content monitoring is carried out to each sub-region, the meteorological data of meteorological station is combined, to analyze the evapotranspiration of each sub-region, combined with evapotranspiration, rainfall, vegetation water absorption, to analyze whether the water content of sub-region meets the requirement of plant growth, and output irrigation index reflecting the current demand irrigation degree of sub-region, while also image shooting and multispectral data acquisition are carried out to sub-region, generate greening evaluation index for reflecting the growth condition of sub-region vegetation, combined with greening evaluation index and irrigation index analysis, obtain sub-region repair index, finally output sub-region repair priority area, priority is carried out targeted vegetation damage repair to repair priority area, greatly improve repair efficiency.

Description

TECHNICAL FIELD

[0001] The present application relates to the technical field of mine vegetation, in particular to a method and device for repairing damaged vegetation in an alpine and arid coal mining area. BACKGROUND

[0002] The damaged vegetation in an alpine and arid coal mining area refers to the phenomenon of decreased vegetation coverage, reduced plant species and degraded ecosystem function caused by human activities and natural environment during coal mining. Due to the special climate conditions, the annual average precipitation in these areas is limited, and the soil moisture evaporates quickly. In addition, phenomena such as surface subsidence, soil erosion and pollution caused by coal mining. During the mining process, the soil is stripped and the vegetation is damaged, which leads to changes in microclimate and further affects the physical and chemical properties of the soil and biodiversity. In addition, the dust, noise and wastewater discharge around the mining area also cause serious damage to the habitat, further threatening the growth of local vegetation and ecological balance.

[0003] The damaged vegetation will lead to a decrease in greening rate, which may trigger various disaster risks. Currently, there is no specific way to repair damaged vegetation in an alpine and arid coal mining area. Paving vegetation damage repair not only consumes a lot of manpower and resources, but also wastes a lot of time and is low in efficiency.

[0004] The above information disclosed in the background section is only used to enhance the understanding of the background of the present disclosure, and therefore it can include information that does not constitute prior art known to those of ordinary skill in the art. SUMMARY

[0005] The purpose of the present application is to provide a method and device for repairing damaged vegetation in an alpine and arid coal mining area to solve the problems raised in the background.

[0006] To achieve the above-mentioned purpose, the present application provides the following technical solutions:

[0007] A method for repairing damaged vegetation in an alpine and arid coal mining area, comprising the following specific steps:

[0008] S1, query the maximum latitude, minimum latitude, maximum longitude and minimum longitude of the coal mining area through the geographic information system; mark the points corresponding to the maximum latitude, minimum latitude, maximum longitude and minimum longitude on the electronic map, respectively, and the four points form a quadrilateral area, which is taken as the coal mining area;

[0009] S2, the coal mine area is divided into N sub-regions, and the sub-regions are numbered, the three-dimensional point cloud data of the sub-regions is collected and stored in the database, the three-dimensional model of the sub-region is drawn through Surpac three-dimensional software, the soil moisture sensor is set in each sub-region to collect the water content, the meteorological data is collected, the meteorological data includes net radiation, soil heat flux density, dry air specific heat capacity, wind speed, temperature, relative humidity data, the three-dimensional model of the sub-region is imported into ArcGIS Pro, the water content and meteorological data are imported into ArcGIS Pro, and the corresponding spatial points are associated, and the water content and meteorological data of the three-dimensional model of the sub-region are visualized through the 3D Analyst module;

[0010] S3, the temperature data is analyzed, the saturated vapor pressure is generated, and the actual vapor pressure and the slope of the saturated vapor pressure are further generated, the actual vapor pressure, the saturated vapor pressure, the slope of the saturated vapor pressure and the meteorological data are analyzed, the mean evapotranspiration is generated, and the irrigation index is generated by analyzing the mean evapotranspiration;

[0011] S4, the image of the sub-region is shot and the multispectral data is obtained, the vegetation image data of the sub-region is obtained through the edge algorithm, the vegetation coverage is generated by analyzing the correlation of the image data, the vegetation health index is generated by analyzing the correlation of the multispectral data, and the greening evaluation index is generated by analyzing the correlation of the vegetation coverage and the vegetation health index;

[0012] S5, the greening evaluation index LP and the irrigation index GZ are analyzed, the sub-region repair index XF is generated, the sub-region repair index data is imported into ArcGIS Pro, is associated with the corresponding spatial point, a repair model is constructed, and a sub-region repair priority area is generated.

[0013] Further, the meteorological data is the average value of the day, the temperature T is analyzed to generate the saturated vapor pressure e s , and the formula is:

[0014]

[0015] Wherein, T is the average value of the temperature monitored by the weather station on the day, the saturated vapor pressure e s is used to reflect the maximum pressure that water vapor can exist in the air at the current temperature;

[0016] The saturated vapor pressure e s and the relative humidity H are analyzed to generate the actual vapor pressure e a , and the formula is:

[0017]

[0018] Actual vapor pressure e a Actual vapor pressure e refers to the pressure of water vapor actually existing in the air under the current environment, which is used to reflect the actual content of water vapor in the air.

[0019] Correlation analysis is performed on the temperature T and the saturated vapor pressure e s , and the slope Δ of the saturated vapor pressure is generated, and the formula is as follows:

[0020]

[0021] The slope Δ of the saturated vapor pressure is used to reflect the relationship between the saturated vapor pressure and the temperature.

[0022] Further, correlation analysis is performed on the actual vapor pressure, the saturated vapor pressure, the slope of the saturated vapor pressure, and the meteorological data, and the mean evapotranspiration ET c is generated, and the formula is as follows:

[0023]

[0024] The mean evapotranspiration ET c is used to reflect the total amount of water lost due to evaporation in the sub-region on the same day, wherein R n is the net radiation, G is the soil heat flux density, γ is the specific heat capacity of dry air, and u is the wind speed.

[0025] Correlation analysis is performed on the mean evapotranspiration ET c , and the irrigation index GZ is generated, and the formula is as follows:

[0026]

[0027] Wherein, W is the water required in the growth cycle of the vegetation in the sub-region, C is the water content in the sub-region, R e is the rainfall in the sub-region on the same day, and the irrigation index GZ is used to reflect the degree of irrigation required for the growth of vegetation in the sub-region on the same day.

[0028] Further, correlation analysis is performed on the sub-region, and vegetation data i m in the sub-region is collected, wherein the vegetation data i m is the number of the mth vegetation in the ith sub-region, and correlation analysis is performed on the vegetation data i m , and the water required in the growth cycle of the vegetation in the sub-region W is generated, and the formula is as follows:

[0029]

[0030] Wherein, C m is the water required for the growth of the mth vegetation per day.

[0031] Further, the unmanned aerial vehicle is equipped with a multi-spectral sensor, multi-spectral data of a sub-region is captured through the multi-spectral sensor, the multi-spectral data contains information of different wave bands, including near-infrared reflectance NIR and red reflectance RED, correlation analysis is performed on the multi-spectral data, and vegetation health index JP is generated according to the following formula:

[0032]

[0033] The vegetation health index JP is used for reflecting the health degree of the vegetation in the sub-region.

[0034] Correlation analysis is performed on the image data, and vegetation coverage Y is generated according to the following formula:

[0035]

[0036] Wherein, k is the number of pixel points occupied by the vegetation in the image data, and K is the number of pixel points occupied by the sub-region in the image.

[0037] Correlation analysis is performed on the vegetation coverage Y and the vegetation health index JP, and a greening evaluation index LP is generated according to the following formula:

[0038] LP=Y*alpha*JP

[0039] Wherein, alpha is a coverage weight factor, and the value range is (0, 1], the coverage weight factor is the weight degree of the sub-region vegetation coverage, and the greening evaluation index LP is used for reflecting the growth degree of the vegetation in the sub-region.

[0040] Further, correlation analysis is performed on the greening evaluation index LP and the irrigation index GZ, and a sub-region repair index XF is generated according to the following formula:

[0041]

[0042] The sub-region repair index XF is used for reflecting the growth of the vegetation in the sub-region.

[0043] Further, a repair model is constructed, the sub-region repair index XF is arranged in descending order, and the corresponding sub-region is arranged, finally, the sub-region arrangement order is output, and the sub-region with the sub-region arrangement serial number 1 is the sub-region repair priority area.

[0044] The application also provides a vegetation damage repair device for high-cold and arid coal mine areas, which is used for executing a vegetation damage repair method for high-cold and arid coal mine areas.

[0045] The area construction module inquires latitude maximum value, latitude minimum value, longitude maximum value and longitude minimum value of the coal mining area through a geographic information system; points corresponding to the latitude maximum value, the latitude minimum value, the longitude maximum value and the longitude minimum value are marked on an electronic map respectively, four points form a quadrilateral area, and the quadrilateral area is taken as the coal mine area;

[0046] The weather data acquisition module divides the coal mine area into N sub-areas, numbers the sub-areas, acquires three-dimensional point cloud data of the sub-areas and stores the three-dimensional point cloud data in a database, draws a three-dimensional model of the sub-areas through Surpac three-dimensional software, sets a soil moisture sensor in each sub-area to acquire water content, acquires weather data, the weather data includes net radiation, soil heat flux density, dry air specific heat capacity, wind speed, temperature and relative humidity data, imports the three-dimensional model of the sub-areas into ArcGIS Pro, imports the water content and the weather data into ArcGIS Pro, correlates the water content and the weather data to corresponding spatial points, and visualizes the water content and the weather data of the three-dimensional model of the sub-areas through a 3D Analyst module;

[0047] The data analysis module is used for correlation analysis on temperature data, generation of saturated vapor pressure, further generation of actual vapor pressure and slope of the saturated vapor pressure, correlation analysis on the actual vapor pressure, the saturated vapor pressure, the slope of the saturated vapor pressure and weather data, generation of mean evapotranspiration, correlation analysis on the mean evapotranspiration, and generation of an irrigation index;

[0048] The greening evaluation analysis module is used for image shooting of the sub-areas and acquisition of multispectral data, acquisition of vegetation image data of the sub-areas through an edge algorithm, correlation analysis on the image data, generation of vegetation coverage, correlation analysis on the multispectral data, generation of a vegetation health index, correlation analysis on the vegetation coverage and the vegetation health index, and generation of a greening evaluation index;

[0049] The to-be-repaired output module performs correlation analysis on the greening evaluation index LP and the irrigation index GZ, generates a sub-area to-be-repaired index XF, imports the sub-area to-be-repaired index data into ArcGIS Pro, correlates the sub-area to-be-repaired index data to corresponding spatial points, constructs a repair model, and generates a sub-area repair priority area.

[0050] Compared with the prior art, the present application has the following beneficial effects:

[0051] The present application divides the coal mine area into sub-regions, monitors the water content of each sub-region, analyzes the evapotranspiration of each sub-region in combination with meteorological data, analyzes whether the water content of the sub-region meets the requirements of plant growth in combination with evapotranspiration, rainfall and vegetation water absorption, outputs an irrigation index reflecting the current irrigation requirement degree of the sub-region, simultaneously takes images of the sub-region and collects multispectral data to generate a greening evaluation index reflecting the growth status of the vegetation in the sub-region, analyzes the greening evaluation index and the irrigation index to obtain a sub-region repair index, and finally outputs a sub-region repair priority area, so that the repair efficiency is greatly improved. BRIEF DESCRIPTION OF DRAWINGS

[0052] Figure 1 It is a whole method flowchart of the present application;

[0053] Figure 2 It is a whole system flowchart of the present application;

[0054] Figure 3 It is a fitting curve diagram of the greening evaluation index-repair index of the present application;

[0055] Figure 4 It is a fitting curve diagram of the irrigation index-repair index of the present application. DETAILED DESCRIPTION

[0056] In order to make the purpose, technical scheme and advantages of the present application more clear, the present application is further described in detail below in combination with specific embodiments.

[0057] It should be noted that, unless otherwise defined, the technical terms or scientific terms used in the present application should be understood as the general meaning understood by those skilled in the art to which the present application belongs. The "first", "second" and similar words used in the present application do not represent any order, quantity or importance, but are only used to distinguish different components. "Include" or "contain" and similar words mean that the elements or objects before the word cover the elements or objects listed after the word and their equivalents, without excluding other elements or objects. "Connected" or "connected" and similar words are not limited to physical or mechanical connection, but can include electrical connection, whether direct or indirect. "Up", "down", "left", "right" and the like are only used to represent relative positional relationship, when the absolute position of the described object changes, the relative positional relationship may also change accordingly.

[0058] Embodiment:

[0059] Please refer to Figure 1 The present application provides a technical scheme:

[0060] A kind of method for repairing vegetation damage in alpine arid coal mine area, specific steps include:

[0061] Step 1, query the latitude maximum value, latitude minimum value, longitude maximum value and longitude minimum value of the coal mining area through geographic information system;Mark the point corresponding to the latitude maximum value, latitude minimum value, longitude maximum value and longitude minimum value on the electronic map respectively, and the four points form a quadrilateral region, which is taken as the coal mine area;

[0062] Step 2, divide the coal mine area into N sub-regions, and number the sub-regions, collect the three-dimensional point cloud data of the sub-regions and store them in the database, draw the three-dimensional model of the sub-region by Surpac three-dimensional software, set up soil moisture sensor in each sub-region to collect water content, collect meteorological data, including net radiation, soil heat flux density, dry air specific heat capacity, wind speed, temperature, relative humidity data, import the three-dimensional model of the sub-region into ArcGIS Pro, ArcGIS Pro software can combine three-dimensional model and meteorological data, etc., more intuitive display of the relevant data of the point, import the water content and meteorological data into ArcGIS Pro, and correlate to the corresponding space point, visualize the water content and meteorological data of the three-dimensional model of the sub-region by 3D Analyst module;

[0063] At least three groups of sensors are arranged in each sub-region of the coal mine area, each group containing 3 sensors. Each sensor is buried at least 15-30 cm deep in the soil to obtain accurate soil moisture data. The sensor is a soil moisture sensor based on resistivity, and the obtained water content is averaged.

[0064] Step 3, correlation analysis is performed on the temperature data to generate saturated vapor pressure, and further generate actual vapor pressure and slope of saturated vapor pressure, correlation analysis is performed on the actual vapor pressure, saturated vapor pressure, slope of saturated vapor pressure and meteorological data to generate mean evapotranspiration, correlation analysis is performed on the mean evapotranspiration to generate irrigation index;

[0065] The meteorological data are all average values of the day, and correlation analysis is performed on the temperature T to generate saturated vapor pressure e s , and the formula is:

[0066]

[0067] Wherein, T is the average value of the temperature monitored by the weather station on the day, and saturated vapor pressure e s is used to reflect the maximum pressure that water vapor can exist in the air at the current temperature;

[0068] Correlation analysis is performed on the saturated vapor pressure e sCorrelation analysis is performed on the actual vapor pressure e and the relative humidity H to generate the actual vapor pressure e a , and the formula is:

[0069]

[0070] The actual vapor pressure e a refers to the actual pressure of water vapor existing in the air under the current environment, which is used to reflect the actual content of water vapor in the air.

[0071] Correlation analysis is performed on the temperature T and the saturated vapor pressure e s to generate the slope Δ of the saturated vapor pressure, and the formula is:

[0072]

[0073] The slope Δ of the saturated vapor pressure is used to reflect the relationship between the saturated vapor pressure and the temperature.

[0074] Correlation analysis is performed on the actual vapor pressure, the saturated vapor pressure, the slope of the saturated vapor pressure, and the meteorological data to generate the mean evapotranspiration ET c , and the formula is:

[0075]

[0076] The mean evapotranspiration ET c is used to reflect the total amount of water loss in the sub-region due to evaporation on the same day, where R n is the net radiation, G is the soil heat flux density, γ is the specific heat capacity of dry air, and u is the wind speed.

[0077] Evapotranspiration mainly includes the conversion of water in the soil into water vapor, and the release of water vapor from the stomata of plant leaves into the atmosphere after the plant roots absorb water. The evapotranspiration of plants is affected by temperature, humidity, wind speed, and soil moisture conditions, while the evaporation of soil depends on soil moisture, temperature, and environmental conditions. Overall, the greater the slope Δ of the saturated vapor pressure and the specific heat capacity γ of dry air, the more difficult it is for water to evaporate. The wind speed and the water vapor pressure difference affect the water transfer rate. Further comprehensive consideration of factors such as solar radiation, temperature, humidity, and wind speed determines the mean evapotranspiration ET c of the sub-region on the same day.

[0078] Correlation analysis is performed on the mean evapotranspiration ET c to generate the irrigation index GZ, and the formula is:

[0079]

[0080] where W is the water required for the growth cycle of the vegetation in the sub-region, C is the water content of the sub-region, and R eFor the sub-area of the day rainfall, irrigation index GZ is used to reflect the sub-area of the day the degree of vegetation growth needs irrigation. The existing measured water content of the sub-area is added to the rainfall and then subtracted by the evapotranspiration, which is the water content of the sub-area on the day, which is compared with the water required in the growth cycle of the sub-area vegetation, and the irrigation index GZ is output. When the irrigation index GZ is greater than 1, it means that the water content of the sub-area is sufficient, and when the irrigation index GZ is less than 1, it means that the water content of the sub-area is insufficient. The greater the irrigation index GZ, the better the water content of the region, and the better the vegetation growth condition, and the less the need for vegetation damage repair.

[0081] Among them, the correlation analysis is carried out on the sub-area, and the vegetation data i m of the sub-area is collected m The number of the mth vegetation in the ith sub-area is analyzed, and the correlation analysis is carried out on the vegetation data i m The water required in the growth cycle of the sub-area vegetation W is generated, and the formula is as follows:

[0082]

[0083] Among them, C m is the water required by the mth vegetation in the current growth cycle per day.

[0084] Step 4, image shooting and multi-spectral data acquisition of the sub-area, vegetation image data of the sub-area is obtained through edge algorithm, correlation analysis is carried out on the image data, vegetation coverage rate is generated, correlation analysis is carried out on the multi-spectral data, vegetation health index is generated, correlation analysis is carried out on the vegetation coverage rate and the vegetation health index, and the green evaluation index is generated;

[0085] The unmanned aerial vehicle with high stability, long endurance and strong load capacity is selected, which is suitable for carrying multi-spectral sensor and high-resolution camera for data collection. The multi-spectral sensor can obtain five waveband information of red light, blue light, green light, red edge and near infrared of the sub-area. The brightness and color of the multi-spectral sensor are calibrated to ensure the consistency of the data under different time and environmental conditions. The image data collected by each flight should include high-resolution image data and multi-spectral data, and the information such as time, weather condition and GPS coordinate of each flight should be recorded. Agisoft Metashape software is used to process the multi-spectral data and generate the vegetation health index.

[0086] The unmanned aerial vehicle carries multi-spectral sensor, and the multi-spectral data of the sub-area is captured by the multi-spectral sensor. The multi-spectral data contains information of different wavebands, including near-infrared reflectance NIR and red light reflectance RED. The correlation analysis is carried out on the multi-spectral data, and the vegetation health index JP is generated, and the formula is as follows:

[0087]

[0088] The vegetation health index JP is used to reflect the health degree of the vegetation in the sub-region;

[0089] Healthy plants absorb light energy in the red light band for photosynthesis, so the reflectivity is low. In the near-infrared band, the reflectivity of the plants is relatively high. By combining the reflectivity of the two bands, JP can effectively reflect the health status of the vegetation. Similarly, the leaves of healthy plants contain a large amount of chlorophyll, which makes them absorb significantly in the red light band and reflect strongly in the near-infrared band. Therefore, the NIR value of healthy vegetation is higher than the RED value, resulting in the value of JP tends to 1. The value of JP is between -1 and 1, the negative value corresponds to no vegetation in the sub-region, and the value close to 0 usually indicates sparse vegetation in the sub-region, and the JP value close to 1 indicates dense or healthy vegetation.

[0090] Correlation analysis is performed on the image data to generate the vegetation coverage Y, and the formula is:

[0091]

[0092] wherein k is the number of pixel points occupied by the vegetation in the image data, and K is the number of pixel points occupied by the sub-region in the image;

[0093] Correlation analysis is performed on the vegetation coverage Y and the vegetation health index JP to generate the greening evaluation index LP, and the formula is:

[0094] LP=Y*α*JP

[0095] wherein α is a coverage weight factor, and the value range is (0, 1], the coverage weight factor is the weight degree of the vegetation coverage of the sub-region. The higher the vegetation coverage Y and the vegetation health index JP, the higher the greening evaluation index LP, and the higher the greening evaluation index LP, the better the growth of the vegetation in the sub-region, that is, the less repair is needed.

[0096] Step 5, correlation analysis is performed on the greening evaluation index LP and the irrigation index GZ to generate the sub-region repair index XF, and the sub-region repair index data is imported into ArcGIS Pro and associated with the corresponding spatial points to construct a repair model and generate a sub-region repair priority area.

[0097] Correlation analysis is performed on the greening evaluation index LP and the irrigation index GZ to generate the sub-region repair index XF, and the formula is:

[0098]

[0099] The sub-region to-be-repaired index XF is used to reflect the vegetation growth of the sub-region, wherein, the beta is a weight factor of the irrigation index, used to represent the influence degree of the irrigation index on the sub-region to-be-repaired index, the greater the value is, the higher the influence degree is, in the sub-region to-be-repaired index XF, the higher the irrigation index is, the better the soil moisture content of the sub-region is, the more smoothly the vegetation grows, and the less the damage repair is needed, the constant c is used to adjust the sensitivity of the greening evaluation index LP, and can prevent the formula from being invalid when the LP value is 0, the greater the greening evaluation index LP is, the better the vegetation growth of the sub-region is, and the less the damage repair is needed, the greater the irrigation index GZ is, the better the water content of the region is, and the better the vegetation growth condition is, and the less the vegetation damage repair is needed, therefore, the smaller the value of the sub-region to-be-repaired index XF is, the less the damage repair of the sub-region is needed.

[0100] The repair model is constructed, the sub-region to-be-repaired index XF is arranged in descending order, and the corresponding sub-region is arranged, and finally the sub-region arrangement order is output, and the sub-region with the sub-region arrangement serial number 1 is the sub-region repair priority area.

[0101] Reference Figure 2 The application also provides a vegetation damage repair device for high-cold and arid coal mine areas, which is used to execute a vegetation damage repair method for high-cold and arid coal mine areas.

[0102] The region construction module queries the latitude maximum value, the latitude minimum value, the longitude maximum value and the longitude minimum value of the coal mining region through a geographic information system; marks the points corresponding to the latitude maximum value, the latitude minimum value, the longitude maximum value and the longitude minimum value on an electronic map respectively, and the four points form a quadrilateral region, and the quadrilateral region is taken as the coal mine area.

[0103] The weather data acquisition module divides the coal mine area into N sub-regions, numbers the sub-regions, acquires the three-dimensional point cloud data of the sub-regions and stores the data in a database, draws a three-dimensional model of the sub-region through Surpac three-dimensional software, sets a soil moisture sensor in each sub-region to acquire the water content, acquires weather data, the weather data includes net radiation, soil heat flux density, dry air specific heat capacity, wind speed, temperature and relative humidity data, imports the three-dimensional model of the sub-region into ArcGIS Pro, imports the water content and weather data into ArcGIS Pro, associates the data to corresponding spatial points, and visualizes the water content and weather data of the three-dimensional model of the sub-region through a 3D Analyst module.

[0104] The data analysis module is used for correlation analysis of temperature data, generation of saturated vapor pressure, and further generation of actual vapor pressure and slope of saturated vapor pressure, correlation analysis of actual vapor pressure, saturated vapor pressure, slope of saturated vapor pressure, and meteorological data, generation of mean evapotranspiration, correlation analysis of mean evapotranspiration, and generation of irrigation index.

[0105] The green evaluation analysis module is used for shooting images of sub-regions and obtaining multispectral data, obtaining vegetation image data of the sub-regions through an edge algorithm, performing correlation analysis on the image data to generate vegetation coverage, performing correlation analysis on the multispectral data to generate vegetation health index, and performing correlation analysis on the vegetation coverage and the vegetation health index to generate a green evaluation index.

[0106] The to-be-repaired output module performs correlation analysis on the green evaluation index LP and the irrigation index GZ to generate a sub-region to-be-repaired index XF, and imports the sub-region to-be-repaired index data into ArcGIS Pro, associates it to the corresponding spatial point, constructs a repair model, and generates a sub-region repair priority area.

[0107] In this embodiment, Table 1 is a statistical situation of 20 groups of green evaluation indexes, irrigation indexes, and sub-region to-be-repaired indexes, as shown in the following table:

[0108] Table 1

[0109]

[0110]

[0111] According to Table 1, a fitting curve graph of the green evaluation index-to-be-repaired index and a fitting curve graph of the irrigation index-to-be-repaired index are fitted.

[0112] Referring to Figure 3 The graph shows the relationship between the "to-be-repaired index" and the "green evaluation index". By comparing these two indicators, it can be seen that they show a clear negative correlation trend. The horizontal axis represents the green evaluation index, ranging from 0 to 3, and the vertical axis represents the to-be-repaired index, ranging from 0.5 to 2.5. The black squares represent the to-be-repaired index data points, which gradually decrease as the green evaluation index increases.

[0113] Referring to Figure 4 The graph presents the relationship between the "to-be-repaired index" and the "irrigation index", showing a clear negative correlation trend between the two. The horizontal axis represents the irrigation index, ranging from 1 to 2, and the vertical axis represents the to-be-repaired index, ranging from 0.5 to 2.5. The black squares in the graph represent the actual measurements of the to-be-repaired index, which significantly decrease as the irrigation index increases.

[0114] The above formulas are all dimensionless values calculated, the formula is obtained by collecting a large amount of data to simulate the most recent real situation, and the preset parameters in the formula are set by a person skilled in the art according to the actual situation.

[0115] The above embodiments can be implemented wholly or partially by software, hardware, firmware or any other combination. When implemented by software, the above embodiments can be implemented wholly or partially in the form of a computer program product. Those skilled in the art can realize that the units and algorithm steps of the examples described in connection with the embodiments disclosed herein can be realized by electronic hardware or a combination of computer software and electronic hardware. Whether the functions are executed by hardware or software methods depends on the specific application and design constraints of the technical solutions.

[0116] The units described as separate components can or can not be physically separated, and the components shown as units can or can not be physical units, which can be located in one place or distributed on multiple network units. Part or all of the units can be selected to achieve the purpose of the embodiments according to actual needs.

[0117] The above is only a specific implementation of the present application, but the protection scope of the present application is not limited thereto, any person skilled in the art can easily think of changes or replacements within the technical range disclosed in the present application, which should be covered within the protection scope of the present application.

Claims

1. A method for repairing vegetation damage in high-cold and arid coal mine areas, characterized in that, The specific steps include: S1, query the boundary point coordinates of the coal mining area, based on the boundary point coordinates to form a quadrilateral space area, define the scope of the coal mine area, and divide the coal mine area in the digital space; S2, divide the entire area into N sub-regions using the equal division algorithm, assign a unique number to each sub-region, collect three-dimensional point cloud data for the sub-regions, draw a three-dimensional model of the sub-regions using Surpac three-dimensional software, collect water content and weather data for the sub-regions, import the three-dimensional model of the sub-regions, water content and weather data into ArcGIS Pro, use spatial coordinates to achieve geographic registration, and visualize the water content and weather data of the three-dimensional model of the sub-regions using the 3D Analyst module; S3, perform temperature correlation analysis based on the weather data of the sub-regions, derive the air saturation vapor pressure, and then calculate the actual vapor pressure and the slope of the saturation vapor pressure with respect to temperature, perform multivariate correlation analysis on the actual vapor pressure, saturation vapor pressure, slope and weather data, calculate the mean evapotranspiration of the sub-regions, and further calculate the irrigation index based on the mean evapotranspiration results; S4, collect images and multispectral data of the vegetation in the sub-regions; use edge detection algorithm to preprocess the vegetation images, extract effective vegetation boundaries and feature information, perform correlation analysis based on image data, calculate the vegetation coverage of each sub-region, and calculate the vegetation health index based on multispectral data; integrate the vegetation coverage and vegetation health index, and use statistical correlation model to calculate and generate the greening evaluation index; S5, construct a sub-region repair index model, correlate the greening evaluation index with the irrigation index, calculate and generate the repair index of each sub-region; import the repair index data into ArcGIS Pro, correlate it to the corresponding spatial point, construct a repair model, and generate a sub-region repair priority area.

2. The method for repairing vegetation damage in high-cold and arid coal mine areas according to claim 1, characterized in that: The meteorological data are all average values of the day, and the correlation analysis of temperature T generates saturated vapor pressure e s , and the formula is: wherein T is the average value of the temperature monitored by the weather station on the day, the saturated vapor pressure e s for reflecting the maximum pressure at which water vapor can exist in air at the current temperature; Correlation analysis was performed on the saturated vapor pressure e s and the relative humidity H, to generate the actual vapor pressure e a , based on the formula: actual vapor pressure e a actual vapor pressure e refers to the pressure of water vapor actually existing in the air under the current environment, and is used to reflect the actual content of water vapor in the air; The temperature T and the saturated vapor pressure e s A correlation analysis was performed to generate the slope Δ of the saturated vapor pressure, according to the formula: The slope of the saturation vapor pressure is used to reflect the relationship between the saturation vapor pressure and the temperature.

3. The method for repairing vegetation damage in high-cold and arid coal mine areas according to claim 2, characterized in that: Correlation analysis is performed on actual vapor pressure, saturated vapor pressure, slope of saturated vapor pressure, and meteorological data to generate the mean evapotranspiration ET c , and the formula used is: Mean evapotranspiration ET c to reflect the total amount of water lost due to evaporation in the sub-region on the day, wherein R n is the net radiation, G is the soil heat flux density, γ is the specific heat capacity of dry air, and u is the wind speed; The correlation analysis is performed on the mean evapotranspiration ET c An irrigation index GZ is generated according to the formula: Wherein, W is the water required in the vegetation growth cycle of the sub-region, C is the water content of the sub-region, R e is the rainfall of the sub-region on the day, and the irrigation index GZ is used to reflect the degree of irrigation required by the vegetation growth of the sub-region on the day.

4. The method for repairing vegetation damage in high-cold and arid coal mine areas according to claim 3, characterized in that: Correlation analysis is performed on the sub-regions, and vegetation data i in the sub-regions is collected m , the vegetation data i m is the number of the mth vegetation in the ith sub-region, and correlation analysis is performed on the vegetation data i m to generate the required water amount W in the growth cycle of the vegetation in the sub-region, and the formula is as follows: where C m is the water requirement for the mth plant for each day in the current growth cycle.

5. The method for repairing vegetation damage in high-cold and arid coal mine areas according to claim 1, characterized in that: The unmanned aerial vehicle is equipped with a multispectral sensor, which captures multispectral data of the sub-region through the multispectral sensor. The multispectral data contains information of different wavebands, including near-infrared reflectance NIR and red reflectance RED. Correlation analysis is performed on the multispectral data to generate the vegetation health index JP, which is based on the formula: The vegetation health index JP is used to reflect the health degree of the vegetation in the sub-region; Correlation analysis is performed on the image data to generate the vegetation coverage Y, which is based on the formula: Where k is the number of pixel points occupied by vegetation in the image data, and K is the number of pixel points occupied by the sub-region in the image; Correlation analysis is performed on the vegetation coverage Y and the vegetation health index JP to generate the greening evaluation index LP, which is based on the formula: LP=Y*α*JP Where α is the coverage weight factor, which is in the range of (0, 1], and is the weight degree of the vegetation coverage in the sub-region. The greening evaluation index LP is used to reflect the growth quality of the vegetation in the sub-region.

6. The method for repairing vegetation damage in high-cold and arid coal mine areas according to claim 1, characterized in that: Correlation analysis is performed on the greening evaluation index LP and the irrigation index GZ to generate the sub-region repair index XF, which is based on the formula: The sub-region repair index XF is used to reflect the growth of the vegetation in the sub-region.

7. The method for repairing vegetation damage in high-cold and arid coal mine areas according to claim 6, characterized in that: The repair model is constructed, the sub-area to-be-repaired index XF is arranged in descending order, and the corresponding sub-area is arranged, and finally the sub-area arrangement order is output, and the sub-area arrangement serial number of the top 20% of the sub-area is the sub-area repair priority area.

8. A device for repairing vegetation damage in alpine arid coal mine areas, used to perform the method for repairing vegetation damage in alpine arid coal mine areas according to claim 1, characterized in that, Comprise: The area construction module queries the latitude maximum value, latitude minimum value, longitude maximum value and longitude minimum value of the coal mining area through the geographic information system; Mark the points corresponding to the latitude maximum value, latitude minimum value, longitude maximum value and longitude minimum value on the electronic map, and the four points form a quadrilateral area, which is the coal mine area; The weather data acquisition module divides the coal mine area into N sub-areas, numbers the sub-areas, acquires and stores the three-dimensional point cloud data of the sub-areas in the database, draws the three-dimensional model of the sub-areas by Surpac three-dimensional software, sets soil moisture sensors in each sub-area to collect water content, acquires weather data, the weather data includes net radiation, soil heat flux density, dry air specific heat capacity, wind speed, temperature, relative humidity data, imports the three-dimensional model of the sub-areas into ArcGIS Pro, imports the water content and weather data into ArcGIS Pro, and associates them with the corresponding spatial points. The water content and weather data of the three-dimensional model of the sub-area are visualized by the 3D Analyst module; The data analysis module is used for correlation analysis of temperature data, generation of saturated vapor pressure, and further generation of actual vapor pressure and slope of saturated vapor pressure, correlation analysis of actual vapor pressure, saturated vapor pressure, slope of saturated vapor pressure and weather data, generation of mean evapotranspiration, correlation analysis of mean evapotranspiration, and generation of irrigation index; The green evaluation analysis module is used for taking images of the sub-area and acquiring multispectral data, obtaining vegetation image data of the sub-area by edge algorithm, performing correlation analysis on the image data to generate vegetation coverage, performing correlation analysis on the multispectral data to generate vegetation health index, and performing correlation analysis on the vegetation coverage and vegetation health index to generate green evaluation index; The to-be-repaired output module performs correlation analysis on the green evaluation index LP and the irrigation index GZ to generate the sub-area to-be-repaired index XF, imports the sub-area to-be-repaired index data into ArcGIS Pro, associates it with the corresponding spatial points, constructs a repair model, and generates a sub-area repair priority area.

Images