Data analysis-based stony desertification monitoring data analysis method and system

Abstract

The invention discloses a stony desertification monitoring data analysis method and system based on data analysis, relates to the field of environmental governance, and solves the problem of poor application effect of an existing stony desertification monitoring data analysis method.The method comprises the steps that S1, a stony desertification monitoring area is divided into a plurality of stony desertification monitoring sub-areas, the method comprises the following steps: S1, carrying out historical rock coverage monitoring on each stony desertification monitoring sub-region to obtain surface rock monitoring data, S2, carrying out historical vegetation coverage monitoring on each stony desertification monitoring sub-region according to the surface vegetation monitoring data to obtain surface vegetation monitoring data, and carrying out historical rock coverage monitoring on each stony desertification monitoring sub-region according to the surface vegetation monitoring data. And S3, performing treatment qualification evaluation on the rocky desertification monitoring area according to the surface rock monitoring data and the surface vegetation monitoring data. The rocky desertification monitoring data analysis method is high in accuracy and comprehensiveness.

Description

Technical Field

[0001] The present invention belongs to the field of environmental governance and relates to data analysis technology, in particular to a rocky desertification monitoring data analysis method and system based on data analysis. Background Art

[0002] The existing rocky desertification monitoring data analysis methods have the following defects when evaluating the effects of rocky desertification control: 1. Existing desertification monitoring data analysis methods often only analyze environmental data for the current period, without effectively combining historical data, which leads to inaccurate evaluation of desertification control effects; 2. The existing rocky desertification monitoring data analysis method does not analyze the rock coverage and vegetation coverage, resulting in a lack of comprehensiveness in the analysis results.

[0003] To this end, we propose a rocky desertification monitoring data analysis method and system based on data analysis. Summary of the invention

[0004] In view of the shortcomings of the prior art, the purpose of the present invention is to provide a rocky desertification monitoring data analysis method and system based on data analysis. The present invention aims to improve the accuracy and comprehensiveness of the rocky desertification monitoring data analysis method.

[0005] In order to achieve the above-mentioned purpose, the present invention adopts the following technical scheme: a rocky desertification monitoring data analysis method based on data analysis, comprising the following specific steps: Step S1: Obtain a rocky desertification monitoring area and divide the rocky desertification monitoring area into several rocky desertification monitoring sub-areas, conduct historical rock cover monitoring on each rocky desertification monitoring sub-area, obtain the rock cover evolution coefficient corresponding to the rocky desertification monitoring sub-area according to the monitoring results, and obtain surface rock monitoring data; Step S2: Conduct historical vegetation coverage monitoring for each rocky desertification monitoring sub-area according to the surface vegetation monitoring data, obtain the vegetation coverage evolution coefficient corresponding to each rocky desertification monitoring sub-area according to the monitoring results, and obtain the surface vegetation monitoring data; Step S3: Conducting rocky desertification monitoring analysis on the rocky desertification monitoring area based on the surface rock monitoring data and the surface vegetation monitoring data, and conducting a rocky desertification control qualification assessment on the rocky desertification monitoring area based on the analysis results.

[0006] Furthermore, the step S1 further includes the following specific steps: Step S11: acquiring a rocky desertification monitoring area, dividing the acquired rocky desertification monitoring area into a plurality of rocky desertification monitoring sub-areas of equal area, and selecting a sample rocky desertification sub-area from the acquired plurality of rocky desertification monitoring sub-areas; Step S12: Performing historical monitoring of rock cover in the sample rocky desertification sub-region to obtain a rock cover evolution coefficient corresponding to the sample rocky desertification sub-region; Step S13: acquiring the rock cover evolution coefficient corresponding to each rocky desertification monitoring sub-area respectively to obtain a plurality of rock cover evolution coefficients; Step S14: defining a plurality of rocky desertification monitoring sub-areas and the rock cover evolution coefficient corresponding to each rocky desertification monitoring sub-area as surface rock monitoring data.

[0007] Furthermore, the step S12 further includes the following specific steps: Step S121: in the process of performing historical monitoring of rock cover on the sample rocky desertification sub-area, a number of historical monitoring time points with equal intervals are marked, and the marked rock monitoring historical time points are marked in chronological order as Y1 historical monitoring time point to Ya historical monitoring time point; Step S122: monitoring the rock coverage of the sample rocky desertification sub-area at the Ya historical monitoring time point to obtain the rock coverage index value of the Ya area; Step S123: rock coverage monitoring is performed on the sample rocky desertification sub-areas at the historical monitoring time point Y1 to the historical monitoring time point Ya-1, and the rock coverage index value of the Y1 area to the rock coverage index value of the Ya-1 area is obtained; Step S124: In the existing plane rectangular coordinate system, the historical monitoring time point is marked as the horizontal coordinate, and the regional rock coverage index value is marked as the vertical coordinate to obtain the rock coverage coordinate system. In the rock coverage coordinate system, the Y1 historical monitoring time point and the Y1 regional rock coverage index value are obtained as the vertical coordinate coordinate point to obtain the Y1 rock coverage coordinate point. Similarly, the Ya historical monitoring time point and the Ya regional rock coverage index value are obtained as the vertical coordinate coordinate point to obtain the Ya rock coverage coordinate point. Step S125: in the rock cover coordinate system, obtain the slope of the line connecting the Y1 rock cover coordinate point and the Y2 rock cover coordinate point to obtain the slope value of the Y1 line, obtain the slope of the line connecting the Y2 rock cover coordinate point and the Y3 rock cover coordinate point to obtain the slope value of the Y2 line, and so on, obtain the slope of the line connecting the Ya-1 rock cover coordinate point and the Ya rock cover coordinate point to obtain the slope value of the Ya-1 line; Step S126: sum the slope values ​​of the line connecting Y1 to the line connecting Ya-1 to obtain the rock cover evolution coefficient corresponding to the rocky desertification sub-area of ​​the sample.

[0008] Furthermore, the step S122 further includes the following specific steps: Obtain the satellite remote sensing map of the sample rocky desertification sub-area at the Ya historical monitoring time point to obtain the Ya satellite remote sensing map; Using image recognition algorithms, several different types of cover rocks are extracted from Ya satellite remote sensing maps and named as the first type of cover rock to the mth type of cover rock; The coverage area values ​​corresponding to the first type of covered rock to the mth type of covered rock are obtained through the Ya satellite remote sensing map, and the coverage area values ​​of the first rock to the mth rock are obtained; The area value of the sample rocky desertification sub-area in the Ya satellite remote sensing map is obtained to obtain the area value of the rocky desertification area; A plurality of first type covered rock samples are selected, and the soil material content corresponding to each first type covered rock sample is respectively obtained to obtain a plurality of soil material contents, and the average of the obtained plurality of soil material contents is calculated to obtain the first type rock soil content.

[0009] Furthermore, the step S122 further includes the following specific steps: Obtaining the rock soil content of the second type of covering rock to the mth type of covering rock respectively, and obtaining the rock soil content of the second type of covering rock to the mth type of covering rock; The rock cover index value of the Ya region is obtained by calculating the rock desertification area value, the first rock cover area value to the mth rock cover area value, and the first type of rock soil content to the mth type of rock soil content; The rock cover index value of the Ya area is calculated, and the specific formula is as follows: ; Among them, Yzba is the rock coverage index value of the Ya area, Sfgi is the value of the i-th rock coverage area, Smz is the value of the rocky desertification area area, Nthi is the soil content of the i-th type of rock, and m is the number of types corresponding to the covering rock.

[0010] Furthermore, the step S2 further includes the following specific steps: Step S21: acquiring surface vegetation monitoring data, acquiring multiple rocky desertification monitoring sub-areas respectively according to the surface vegetation monitoring data, and selecting a characteristic rocky desertification sub-area from the acquired multiple rocky desertification monitoring sub-areas; Step S22: Performing historical monitoring of vegetation coverage in the characteristic rocky desertification sub-region to obtain a vegetation coverage evolution coefficient corresponding to the characteristic rocky desertification sub-region; Step S23: respectively obtaining the vegetation coverage evolution coefficient corresponding to each rocky desertification monitoring sub-area to obtain surface vegetation monitoring data.

[0011] Furthermore, the step S22 further includes the following specific steps: Step S221: in the process of historical monitoring of vegetation coverage in the characteristic rocky desertification sub-area, a number of historical monitoring time points with equal intervals are marked, and the marked historical monitoring time points of vegetation are marked in chronological order as historical monitoring time point Z1 to historical monitoring time point Zb; Step S222: monitoring vegetation coverage of the characteristic rocky desertification sub-region at the historical monitoring time point Zb, and obtaining the vegetation coverage index value of the Zb region; Step S223: performing vegetation coverage monitoring on the characteristic rocky desertification sub-regions at the historical monitoring time point Z1 to the historical monitoring time point Zb-1, respectively, to obtain vegetation coverage index values ​​of the Z1 region to the Zb-1 region; Step S224: in the existing plane rectangular coordinate system, the historical monitoring time point is marked as the horizontal coordinate, and the regional vegetation coverage index value is marked as the vertical coordinate to obtain the vegetation coverage coordinate system. In the vegetation coverage coordinate system, the coordinate point with the vertical coordinate of the Z1 historical monitoring time point and the Z1 regional vegetation coverage index value is obtained to obtain the Z1 vegetation coverage coordinate. Similarly, the coordinate point with the vertical coordinate of the Zb historical monitoring time point and the Zb regional vegetation coverage index value is obtained to obtain the Zb vegetation coverage coordinate point. Step S225: in the vegetation coverage coordinate system, obtain the slope of the line connecting the Z1 vegetation coverage coordinate point and the Z2 vegetation coverage coordinate point to obtain the Z1 line slope value, obtain the slope of the line connecting the Z2 vegetation coverage coordinate point and the Z3 vegetation coverage coordinate point to obtain the Z2 line slope value, and so on, obtain the slope of the line connecting the Zb-1 vegetation coverage coordinate point and the Zb vegetation coverage coordinate point to obtain the Zb-1 line slope value; Step S226: summing the slope values ​​of the lines connecting Z1 to Zb-1 to obtain the vegetation cover evolution coefficient corresponding to the characteristic rocky desertification sub-area.

[0012] Furthermore, the step S222 further includes the following specific steps: Obtain the satellite remote sensing map of the characteristic rocky desertification sub-area at the Zb historical monitoring time point to obtain the Zb satellite remote sensing map; Using image recognition algorithms, several different types of ground cover vegetation were extracted from the Zb satellite remote sensing map, and they were named the first type of cover vegetation to the nth type of cover vegetation. Obtain the coverage area values ​​corresponding to the first type of vegetation to the nth type of vegetation through the Zb satellite remote sensing map, and obtain the first vegetation coverage area value to the nth vegetation coverage area value; The area value of the characteristic rocky desertification sub-area in the Zb satellite remote sensing map is obtained to obtain the area value of the rocky desertification area; Selecting a plurality of first-type covered vegetation samples, respectively obtaining the root moisture content corresponding to each first-type covered vegetation sample, obtaining a plurality of root moisture contents, and averaging the obtained plurality of root moisture contents to obtain the root moisture content of the first-type vegetation; Repeat the process of obtaining the root moisture content of the first type of vegetation, and obtain the root moisture content of the vegetation from the second type of covering vegetation to the nth type of covering vegetation, respectively, to obtain the root moisture content of the second type of vegetation to the nth type of vegetation; The rocky desertification area value, the first vegetation coverage area value to the nth vegetation coverage area value and the first type of vegetation root moisture content to the nth type of vegetation root moisture content are calculated to obtain the Zb regional vegetation coverage index value; The vegetation coverage index value of the Zb area is calculated, and the specific formula is as follows: ;

[0013] Among them, Zhsb is the vegetation coverage index value of Zb area, Zfgi is the value of the i-th vegetation coverage area, Smz is the value of the rocky desertification area area, Zhsi is the root moisture content of the i-th type of vegetation, and n is the number of types corresponding to the covered vegetation.

[0014] Furthermore, the step S3 further includes the following specific steps: Obtaining surface rock monitoring data, and obtaining the rock cover evolution coefficient corresponding to each rocky desertification monitoring sub-area based on the surface rock monitoring data; If the rock cover evolution coefficient is greater than or equal to 0, the corresponding rocky desertification monitoring sub-area is judged to be a rock cover growth area; If the rock cover evolution coefficient is less than 0, the corresponding rocky desertification monitoring sub-area is judged to be a rock cover decline area; Obtaining surface vegetation monitoring data, and obtaining the vegetation coverage evolution coefficient corresponding to each rocky desertification monitoring sub-area based on the surface vegetation monitoring data; If the vegetation cover evolution coefficient is greater than 0, the corresponding rocky desertification monitoring sub-area is judged to be a vegetation cover growth area; If the vegetation cover evolution coefficient is less than or equal to 0, the corresponding rocky desertification monitoring sub-area is judged to be a vegetation cover decline area; Mark the rocky desertification monitoring sub-areas that are both in the vegetation cover growth area and the rock cover decline area as the rocky desertification control qualified area, and count the number of rocky desertification control qualified areas, and the number of control qualified areas; Obtain the number of rocky desertification monitoring sub-areas into which the rocky desertification monitoring area is divided, and obtain the number of monitoring sub-areas; Calculate the ratio of the number of qualified areas to the number of monitored sub-areas to obtain the qualified rate of treatment corresponding to the rocky desertification monitoring area; Obtain the preset governance pass rate. If the governance pass rate is greater than or equal to the preset governance pass rate, the rocky desertification governance in the rocky desertification monitoring area is evaluated to be qualified. If the governance pass rate is less than the preset governance pass rate, the rocky desertification governance in the rocky desertification monitoring area is evaluated to be unqualified.

[0015] The rocky desertification monitoring data analysis system based on data analysis includes: Rock monitoring module: Obtain the rocky desertification monitoring area and divide it into several rocky desertification monitoring sub-areas, conduct historical rock cover monitoring on each rocky desertification monitoring sub-area, obtain the rock cover evolution coefficient corresponding to the rocky desertification monitoring sub-area according to the monitoring results, and obtain the surface rock monitoring data; Vegetation monitoring module: Conduct historical vegetation coverage monitoring for each desertification monitoring sub-area based on surface vegetation monitoring data, obtain the vegetation coverage evolution coefficient corresponding to each desertification monitoring sub-area based on the monitoring results, and obtain surface vegetation monitoring data; Monitoring and analysis module: Conduct desertification monitoring and analysis on the desertification monitoring area based on the surface rock monitoring data and surface vegetation monitoring data, and conduct desertification control qualification assessment on the desertification monitoring area based on the analysis results.

[0016] In summary, due to the adoption of the above technical solution, the beneficial effects of the present invention are: 1. The present invention divides the rocky desertification monitoring area into several rocky desertification monitoring sub-areas, and analyzes the rocky desertification monitoring data of each rocky desertification monitoring sub-area at different time points to evaluate the effect of rocky desertification control, thereby ensuring the accuracy of the evaluation results; 2. The present invention simultaneously conducts historical rock cover monitoring and historical vegetation cover monitoring on each rocky desertification monitoring sub-area, and conducts a rocky desertification control qualification assessment on the rocky desertification monitoring area based on the monitoring results. BRIEF DESCRIPTION OF THE DRAWINGS

[0017] In order to facilitate understanding by those skilled in the art, the present invention is further described below with reference to the accompanying drawings.

[0018] Figure 1 It is a diagram of the implementation steps of the present invention; Figure 2 It is the overall system block diagram of the present invention. DETAILED DESCRIPTION

[0019] The technical solution of the present invention will be clearly and completely described below in conjunction with the embodiments. 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.

[0020] Example 1: Please refer to Figure 1 The present invention provides a technical solution: a rocky desertification monitoring data analysis method and system based on data analysis, comprising the following specific steps: Step S1: Acquire a rocky desertification monitoring area and divide the rocky desertification monitoring area into a number of rocky desertification monitoring sub-areas, and perform historical rock cover monitoring on each rocky desertification monitoring sub-area to obtain surface rock monitoring data; The step S1 further includes the following specific steps: Step S11: acquiring a rocky desertification monitoring area, dividing the acquired rocky desertification monitoring area into a plurality of rocky desertification monitoring sub-areas of equal area, and selecting a sample rocky desertification sub-area from the acquired plurality of rocky desertification monitoring sub-areas; Step S12: Performing historical monitoring of rock cover in the sample rocky desertification sub-region to obtain a rock cover evolution coefficient corresponding to the sample rocky desertification sub-region; The step S12 further includes the following specific steps: Step S121: in the process of performing historical monitoring of rock cover on the sample rocky desertification sub-area, a number of historical monitoring time points with equal intervals are marked, and the marked rock monitoring historical time points are marked in chronological order as Y1 historical monitoring time point to Ya historical monitoring time point; Step S122: monitoring the rock coverage of the sample rocky desertification sub-area at the Ya historical monitoring time point to obtain the rock coverage index value of the Ya area; The step S122 further includes the following specific steps: Obtain the satellite remote sensing map of the sample rocky desertification sub-area at the Ya historical monitoring time point to obtain the Ya satellite remote sensing map; Using image recognition algorithms, several different types of cover rocks are extracted from Ya satellite remote sensing maps and named as the first type of cover rock to the mth type of cover rock; The coverage area values ​​corresponding to the first type of covered rock to the mth type of covered rock are obtained through the Ya satellite remote sensing map, and the coverage area values ​​of the first rock to the mth rock are obtained; The area value of the sample rocky desertification sub-area in the Ya satellite remote sensing map is obtained to obtain the area value of the rocky desertification area; Selecting a plurality of first-type covered rock samples, respectively obtaining the soil content corresponding to each of the first-type covered rock samples, obtaining a plurality of soil contents, and averaging the obtained plurality of soil contents to obtain the first-type rock soil content; Obtaining the rock soil content of the second type of covering rock to the mth type of covering rock respectively, and obtaining the rock soil content of the second type of covering rock to the mth type of covering rock; The rock cover index value of the Ya region is obtained by calculating the rock desertification area value, the first rock cover area value to the mth rock cover area value, and the first type of rock soil content to the mth type of rock soil content; The rock cover index value of the Ya area is calculated, and the specific formula is as follows: ; Among them, Yzba is the rock coverage index value of the Ya region, Sfgi is the value of the i-th rock coverage area, Smz is the value of the rocky desertification area, Nthi is the soil content of the i-th type of rock, and m is the number of types corresponding to the covering rock; Step S123: rock coverage monitoring is performed on the sample rocky desertification sub-areas at the historical monitoring time point Y1 to the historical monitoring time point Ya-1, and the rock coverage index value of the Y1 area to the rock coverage index value of the Ya-1 area is obtained; Step S124: In the existing plane rectangular coordinate system, the historical monitoring time point is marked as the horizontal coordinate, and the regional rock coverage index value is marked as the vertical coordinate to obtain the rock coverage coordinate system. In the rock coverage coordinate system, the Y1 historical monitoring time point and the Y1 regional rock coverage index value are obtained as the vertical coordinate coordinate point to obtain the Y1 rock coverage coordinate point. Similarly, the Ya historical monitoring time point and the Ya regional rock coverage index value are obtained as the vertical coordinate coordinate point to obtain the Ya rock coverage coordinate point. Step S125: in the rock cover coordinate system, obtain the slope of the line connecting the Y1 rock cover coordinate point and the Y2 rock cover coordinate point to obtain the slope value of the Y1 line, obtain the slope of the line connecting the Y2 rock cover coordinate point and the Y3 rock cover coordinate point to obtain the slope value of the Y2 line, and so on, obtain the slope of the line connecting the Ya-1 rock cover coordinate point and the Ya rock cover coordinate point to obtain the slope value of the Ya-1 line; Step S126: summing the slope values ​​of the line connecting Y1 to the slope values ​​of the line connecting Ya-1 to obtain the rock cover evolution coefficient corresponding to the rocky desertification sub-region of the sample; Step S13: acquiring the rock cover evolution coefficient corresponding to each rocky desertification monitoring sub-area respectively to obtain a plurality of rock cover evolution coefficients; Step S14: defining a plurality of rocky desertification monitoring sub-areas and a rock cover evolution coefficient corresponding to each rocky desertification monitoring sub-area as surface rock monitoring data; Step S2: Conduct historical vegetation coverage monitoring for each rocky desertification monitoring sub-area according to the surface vegetation monitoring data, obtain the vegetation coverage evolution coefficient corresponding to each rocky desertification monitoring sub-area according to the monitoring results, and obtain the surface vegetation monitoring data; The step S2 further includes the following specific steps: Step S21: acquiring surface vegetation monitoring data, acquiring multiple rocky desertification monitoring sub-areas respectively according to the surface vegetation monitoring data, and selecting a characteristic rocky desertification sub-area from the acquired multiple rocky desertification monitoring sub-areas; Step S22: Performing historical monitoring of vegetation coverage in the characteristic rocky desertification sub-region to obtain a vegetation coverage evolution coefficient corresponding to the characteristic rocky desertification sub-region; The step S22 further includes the following specific steps: Step S221: in the process of historical monitoring of vegetation coverage in the characteristic rocky desertification sub-area, a number of historical monitoring time points with equal intervals are marked, and the marked historical monitoring time points of vegetation are marked in chronological order as historical monitoring time point Z1 to historical monitoring time point Zb; Step S222: monitoring vegetation coverage of the characteristic rocky desertification sub-region at the historical monitoring time point Zb, and obtaining the vegetation coverage index value of the Zb region; The step S222 further includes the following specific steps: Obtain the satellite remote sensing map of the characteristic rocky desertification sub-area at the Zb historical monitoring time point to obtain the Zb satellite remote sensing map; Using image recognition algorithms, several different types of ground cover vegetation were extracted from the Zb satellite remote sensing map, and they were named the first type of cover vegetation to the nth type of cover vegetation. Obtain the coverage area values ​​corresponding to the first type of vegetation to the nth type of vegetation through the Zb satellite remote sensing map, and obtain the first vegetation coverage area value to the nth vegetation coverage area value; The area value of the characteristic rocky desertification sub-area in the Zb satellite remote sensing map is obtained to obtain the area value of the rocky desertification area; Selecting a plurality of first-type covered vegetation samples, respectively obtaining the root moisture content corresponding to each first-type covered vegetation sample, obtaining a plurality of root moisture contents, and averaging the obtained plurality of root moisture contents to obtain the root moisture content of the first-type vegetation; Repeat the process of obtaining the root moisture content of the first type of vegetation, and obtain the root moisture content of the vegetation from the second type of covering vegetation to the nth type of covering vegetation, respectively, to obtain the root moisture content of the second type of vegetation to the nth type of vegetation; The rocky desertification area value, the first vegetation coverage area value to the nth vegetation coverage area value and the first type of vegetation root moisture content to the nth type of vegetation root moisture content are calculated to obtain the Zb regional vegetation coverage index value; The vegetation coverage index value of the Zb area is calculated, and the specific formula is as follows: ; Among them, Zhsb is the vegetation coverage index value of Zb area, Zfgi is the value of the i-th vegetation coverage area, Smz is the value of the rocky desertification area, Zhsi is the root moisture content of the i-th type of vegetation, and n is the number of types corresponding to the covered vegetation; Step S223: performing vegetation coverage monitoring on the characteristic rocky desertification sub-regions at the historical monitoring time point Z1 to the historical monitoring time point Zb-1, respectively, to obtain vegetation coverage index values ​​of the Z1 region to the Zb-1 region; Step S224: In the existing plane rectangular coordinate system, the historical monitoring time point is marked as the horizontal coordinate, and the regional vegetation coverage index value is marked as the vertical coordinate to obtain the vegetation coverage coordinate system. In the vegetation coverage coordinate system, the coordinate point with the vertical coordinate of the Z1 historical monitoring time point and the Z1 regional vegetation coverage index value is obtained to obtain the Z1 vegetation coverage coordinate. Similarly, the coordinate point with the vertical coordinate of the Zb historical monitoring time point and the Zb regional vegetation coverage index value is obtained to obtain the Zb vegetation coverage coordinate point. Step S225: in the vegetation coverage coordinate system, obtain the slope of the line connecting the Z1 vegetation coverage coordinate point and the Z2 vegetation coverage coordinate point to obtain the Z1 line slope value, obtain the slope of the line connecting the Z2 vegetation coverage coordinate point and the Z3 vegetation coverage coordinate point to obtain the Z2 line slope value, and so on, obtain the slope of the line connecting the Zb-1 vegetation coverage coordinate point and the Zb vegetation coverage coordinate point to obtain the Zb-1 line slope value; Step S226: summing the slope values ​​of the lines from Z1 to Zb-1 to obtain the vegetation cover evolution coefficient corresponding to the characteristic rocky desertification sub-area; Step S23: respectively obtaining the vegetation coverage evolution coefficient corresponding to each rocky desertification monitoring sub-area to obtain surface vegetation monitoring data; Step S3: performing rocky desertification monitoring analysis on the rocky desertification monitoring area according to the surface rock monitoring data and the surface vegetation monitoring data, and conducting a rocky desertification control qualification assessment on the rocky desertification monitoring area according to the analysis results; The step S3 further includes the following specific steps: Obtaining surface rock monitoring data, and obtaining the rock cover evolution coefficient corresponding to each rocky desertification monitoring sub-area based on the surface rock monitoring data; If the rock cover evolution coefficient is greater than or equal to 0, the corresponding rocky desertification monitoring sub-area is judged to be a rock cover growth area; If the rock cover evolution coefficient is less than 0, the corresponding rocky desertification monitoring sub-area is judged to be a rock cover decline area; Obtaining surface vegetation monitoring data, and obtaining the vegetation coverage evolution coefficient corresponding to each rocky desertification monitoring sub-area based on the surface vegetation monitoring data; If the vegetation cover evolution coefficient is greater than 0, the corresponding rocky desertification monitoring sub-area is judged to be a vegetation cover growth area; If the vegetation cover evolution coefficient is less than or equal to 0, the corresponding rocky desertification monitoring sub-area is judged to be a vegetation cover decline area; Mark the rocky desertification monitoring sub-areas that are both in the vegetation cover growth area and the rock cover decline area as the rocky desertification control qualified area, and count the number of rocky desertification control qualified areas, and the number of control qualified areas; Obtain the number of rocky desertification monitoring sub-areas into which the rocky desertification monitoring area is divided, and obtain the number of monitoring sub-areas; Calculate the ratio of the number of qualified areas to the number of monitored sub-areas to obtain the qualified rate of treatment corresponding to the rocky desertification monitoring area; Obtain the preset governance pass rate. If the governance pass rate is greater than or equal to the preset governance pass rate, the rocky desertification governance in the rocky desertification monitoring area is evaluated to be qualified. If the governance pass rate is less than the preset governance pass rate, the rocky desertification governance in the rocky desertification monitoring area is evaluated to be unqualified.

[0021] In this application, if corresponding calculation formulas appear, the above calculation formulas are all dimensionless and take their numerical calculations. The weight coefficients, proportional coefficients and other coefficients in the formulas are set to a result value obtained by quantifying each parameter. The size of the weight coefficient and the proportional coefficient can be determined as long as it does not affect the proportional relationship between the parameter and the result value.

[0022] Example 2: Please refer to Figure 2 , based on another concept of the same invention, a rocky desertification monitoring data analysis system based on data analysis is now proposed, the data analysis system includes a rock monitoring module, a vegetation monitoring module, a monitoring and analysis module and a server, the rock monitoring module, the vegetation monitoring module and the monitoring and analysis module are respectively connected to the server, and the server controls the rock monitoring module, the vegetation monitoring module and the monitoring and analysis module respectively; The rock monitoring module acquires the rocky desertification monitoring area and divides the rocky desertification monitoring area into several rocky desertification monitoring sub-areas, and conducts historical rock coverage monitoring on each rocky desertification monitoring sub-area to obtain surface rock monitoring data; The details are as follows: Acquire a rocky desertification monitoring area, divide the acquired rocky desertification monitoring area into a plurality of rocky desertification monitoring sub-areas of equal area, and select a sample rocky desertification sub-area from the acquired plurality of rocky desertification monitoring sub-areas; Conduct historical monitoring of rock cover in the sample rocky desertification sub-area and obtain the rock cover evolution coefficient corresponding to the sample rocky desertification sub-area; The details are as follows: In the process of historical monitoring of rock cover in the sample rocky desertification sub-area, several historical monitoring time points with equal intervals are marked, and the marked rock monitoring historical time points are marked in chronological order as Y1 historical monitoring time point to Ya historical monitoring time point; It should be noted here that: In this application, Y referred to herein is an identifier corresponding to a historical time point of rock monitoring, a is a quantity value corresponding to a historical time point of rock monitoring, and a is an integer greater than 0; In this application, the time interval between two consecutive historical monitoring time points is one year; In the present application, the Y1 historical monitoring time point is the historical monitoring time point with the longest time interval from the current moment, and the Ya historical monitoring time point involved here is the historical monitoring time point corresponding to the current moment.

[0023] The rock coverage of the sample rocky desertification sub-area at the Ya historical monitoring time point was monitored to obtain the rock coverage index value of the Ya area; The details are as follows: Obtain the satellite remote sensing map of the sample rocky desertification sub-area at the Ya historical monitoring time point to obtain the Ya satellite remote sensing map; Using image recognition algorithms, several different types of cover rocks are extracted from Ya satellite remote sensing maps and named as the first type of cover rock to the mth type of cover rock; It should be noted here that: In the present application, m referred to here is the type quantity value corresponding to the covering rock, and m is an integer greater than 0.

[0024] The first type of covering rock involved here may be granite, the second type of covering rock may be sandstone, the third type of covering rock may be limestone, and the fourth type of covering rock may be basalt.

[0025] The coverage area values ​​corresponding to the first type of covered rock to the mth type of covered rock are obtained through the Ya satellite remote sensing map, and the coverage area values ​​of the first rock to the mth rock are obtained; The area value of the sample rocky desertification sub-area in the Ya satellite remote sensing map is obtained to obtain the area value of the rocky desertification area; Selecting a plurality of first-type covered rock samples, respectively obtaining the soil content corresponding to each of the first-type covered rock samples, obtaining a plurality of soil contents, and averaging the obtained plurality of soil contents to obtain the first-type rock soil content; Repeat the process of obtaining the soil content of the first type of rock, and obtain the soil content of the rock from the second type of covering rock to the mth type of covering rock, to obtain the soil content of the second type of rock to the mth type of rock; The rock cover index value of the Ya region is obtained by calculating the rock desertification area value, the first rock cover area value to the mth rock cover area value, and the first type of rock soil content to the mth type of rock soil content; The rock cover index value of the Ya area is calculated, and the specific formula is as follows: ; Among them, Yzba is the rock coverage index value of the Ya region, Sfgi is the value of the i-th rock coverage area, Smz is the value of the rocky desertification area, Nthi is the soil content of the i-th type of rock, and m is the number of types corresponding to the covering rock; It should be noted here that: In the present application, the i-th rock coverage area value involved here may be any rock coverage area value from the first rock coverage area value to the m-th rock coverage area value, and the i-th type of rock soil content involved here may be any type of rock soil content from the first type of rock soil content to the m-th type of rock soil content; In the specific implementation, there are the following experimental test data: The area value of rocky desertification is Smz = 10 square kilometers, and the number of rock types covered is m = 3; The value of rock cover area Sfgi (unit: square kilometers): Sfg1=3, Sfg2=4, Sfg3=2; Soil content of rock of type i Nthi (unit: percentage): Nth1=5%, Nth2=7.5%, Nth3=8%; The regional rock cover index value can be calculated to be 0.839.

[0026] Repeat the process of obtaining the rock coverage index value of the Ya region, and conduct rock coverage monitoring on the sample rocky desertification sub-regions at the historical monitoring time point of Y1 to the historical monitoring time point of Ya-1, and obtain the rock coverage index value of the Y1 region to the rock coverage index value of the Ya-1 region; In the existing plane rectangular coordinate system, the historical monitoring time point is marked as the horizontal coordinate, and the regional rock coverage index value is marked as the vertical coordinate to obtain the rock coverage coordinate system. In the rock coverage coordinate system, the Y1 historical monitoring time point and the Y1 regional rock coverage index value are obtained as the vertical coordinate coordinate point to obtain the Y1 rock coverage coordinate. Similarly, the Ya historical monitoring time point and the Ya regional rock coverage index value are obtained as the vertical coordinate coordinate point to obtain the Ya rock coverage coordinate point. In the rock cover coordinate system, obtain the slope of the line connecting the Y1 rock cover coordinate point and the Y2 rock cover coordinate point to obtain the slope value of the Y1 line, obtain the slope of the line connecting the Y2 rock cover coordinate point and the Y3 rock cover coordinate point to obtain the slope value of the Y2 line, and so on, obtain the slope of the line connecting the Ya-1 rock cover coordinate point and the Ya rock cover coordinate point to obtain the slope value of the Ya-1 line; The slope values ​​of the line connecting Y1 and Ya-1 are summed to obtain the rock cover evolution coefficient corresponding to the rocky desertification sub-area of ​​the sample; Repeat the process of obtaining the rock cover evolution coefficient corresponding to the sample rocky desertification sub-area, respectively obtain the rock cover evolution coefficient corresponding to each rocky desertification monitoring sub-area, and obtain multiple rock cover evolution coefficients; The multiple rocky desertification monitoring sub-areas and the rock cover evolution coefficient corresponding to each rocky desertification monitoring sub-area are defined as surface rock monitoring data; The rock monitoring module acquires the surface rock monitoring data and transmits it to the vegetation monitoring module and the monitoring and analysis module; The vegetation monitoring module conducts historical vegetation coverage monitoring for each rocky desertification monitoring sub-area based on the surface vegetation monitoring data, obtains the vegetation coverage evolution coefficient corresponding to each rocky desertification monitoring sub-area based on the monitoring results, and obtains the surface vegetation monitoring data; The details are as follows: Acquire surface vegetation monitoring data, acquire multiple rocky desertification monitoring sub-areas according to the surface vegetation monitoring data, and select a characteristic rocky desertification sub-area from the acquired multiple rocky desertification monitoring sub-areas; Conduct historical monitoring of vegetation coverage in characteristic rocky desertification sub-regions and obtain vegetation coverage evolution coefficients corresponding to characteristic rocky desertification sub-regions; The details are as follows: In the process of historical monitoring of vegetation coverage in the characteristic rocky desertification sub-area, several historical monitoring time points with equal intervals are marked, and the marked historical vegetation monitoring time points are marked in chronological order as Z1 historical monitoring time point to Zb historical monitoring time point; It should be noted here that: In this application, Z referred to herein is an identifier corresponding to a historical time point of vegetation monitoring, b is a quantity value corresponding to a historical time point of vegetation monitoring, and b is an integer greater than 0; In this application, the time interval between two consecutive historical monitoring time points is one year; In the present application, the Z1 historical monitoring time point is the historical monitoring time point with the longest time interval from the current moment, and the Zb historical monitoring time point involved here is the historical monitoring time point corresponding to the current moment.

[0027] Monitor the vegetation coverage of the characteristic rocky desertification sub-area at the historical monitoring time point Zb, and obtain the vegetation coverage index value of the Zb area; The details are as follows: Obtain the satellite remote sensing map of the characteristic rocky desertification sub-area at the Zb historical monitoring time point to obtain the Zb satellite remote sensing map; Using image recognition algorithms, several different types of ground cover vegetation were extracted from the Zb satellite remote sensing map, and they were named the first type of cover vegetation to the nth type of cover vegetation. It should be noted here that: In the present application, n referred to here is the number of types corresponding to the covering vegetation, and n is an integer greater than 0.

[0028] The first type of covering vegetation involved here may be meadow plants, the second type of covering vegetation may be small and medium-sized trees, and the third type of covering vegetation may be mosses and lichens.

[0029] Obtain the coverage area values ​​corresponding to the first type of vegetation to the nth type of vegetation through the Zb satellite remote sensing map, and obtain the first vegetation coverage area value to the nth vegetation coverage area value; The area value of the characteristic rocky desertification sub-area in the Zb satellite remote sensing map is obtained to obtain the area value of the rocky desertification area; Selecting a plurality of first-type covered vegetation samples, respectively obtaining the root moisture content corresponding to each first-type covered vegetation sample, obtaining a plurality of root moisture contents, and averaging the obtained plurality of root moisture contents to obtain the root moisture content of the first-type vegetation; Repeat the process of obtaining the root moisture content of the first type of vegetation, and obtain the root moisture content of the vegetation from the second type of covering vegetation to the nth type of covering vegetation, respectively, to obtain the root moisture content of the second type of vegetation to the nth type of vegetation; The rocky desertification area value, the first vegetation coverage area value to the nth vegetation coverage area value and the first type of vegetation root moisture content to the nth type of vegetation root moisture content are calculated to obtain the Zb regional vegetation coverage index value; The vegetation coverage index value of the Zb area is calculated, and the specific formula is as follows: ; Among them, Zhsb is the vegetation coverage index value of Zb area, Zfgi is the value of the i-th vegetation coverage area, Smz is the value of the rocky desertification area, Zhsi is the root moisture content of the i-th type of vegetation, and n is the number of types corresponding to the covered vegetation; It should be noted here that: In the present application, the i-th vegetation coverage area value involved here may be any vegetation coverage area value from the first vegetation coverage area value to the n-th vegetation coverage area value, and the i-th type of vegetation root moisture content involved here may be any type of vegetation root moisture content from the first type of vegetation root moisture content to the n-th type of vegetation root moisture content; In the specific implementation, there are the following experimental test data: The area value of rocky desertification is Snz=10 square kilometers, and the number of covered vegetation types is n=3; The value of the i-th vegetation coverage area Sfgi (unit: square kilometers): Zfg1=2, Zfg2=1, Zfg3=3; Root moisture content of vegetation type i Zhsi (unit: percentage): Zhs1=10%, Zhs2=15%, Zhs3=30%; The regional vegetation coverage index value can be calculated to be 0.725.

[0030] Repeat the process of obtaining the vegetation coverage index value of the Zb area, and monitor the vegetation coverage of the characteristic rocky desertification sub-areas from the historical monitoring time point of Z1 to the historical monitoring time point of Zb-1, and obtain the vegetation coverage index value of the Z1 area to the vegetation coverage index value of the Zb-1 area; In the existing plane rectangular coordinate system, the historical monitoring time point is marked as the horizontal coordinate, and the regional vegetation coverage index value is marked as the vertical coordinate to obtain the vegetation coverage coordinate system. In the vegetation coverage coordinate system, the coordinate point with the vertical coordinate of the Z1 historical monitoring time point and the Z1 regional vegetation coverage index value is obtained to obtain the Z1 vegetation coverage coordinate. Similarly, the coordinate point with the vertical coordinate of the Zb historical monitoring time point and the Zb regional vegetation coverage index value is obtained to obtain the Zb vegetation coverage coordinate point. In the vegetation coverage coordinate system, obtain the slope of the line connecting the Z1 vegetation coverage coordinate point and the Z2 vegetation coverage coordinate point to obtain the Z1 line slope value, obtain the slope of the line connecting the Z2 vegetation coverage coordinate point and the Z3 vegetation coverage coordinate point to obtain the Z2 line slope value, and so on, obtain the slope of the line connecting the Zb-1 vegetation coverage coordinate point and the Zb vegetation coverage coordinate point to obtain the Zb-1 line slope value; The slope values ​​of the lines connecting Z1 and Zb-1 are summed to obtain the vegetation cover evolution coefficient corresponding to the characteristic rocky desertification sub-area. Repeat the process of obtaining the vegetation cover evolution coefficient corresponding to the characteristic rocky desertification sub-area, obtain the vegetation cover evolution coefficient corresponding to each rocky desertification monitoring sub-area, and obtain the surface vegetation monitoring data; The vegetation monitoring module acquires the surface vegetation monitoring data and transmits it to the monitoring and analysis module; The monitoring and analysis module conducts rocky desertification monitoring and analysis on the rocky desertification monitoring area based on the surface rock monitoring data and the surface vegetation monitoring data, and conducts a rocky desertification control qualification assessment on the rocky desertification monitoring area based on the analysis results; The details are as follows: Obtaining surface rock monitoring data, and obtaining the rock cover evolution coefficient corresponding to each rocky desertification monitoring sub-area based on the surface rock monitoring data; If the rock cover evolution coefficient is greater than or equal to 0, the corresponding rocky desertification monitoring sub-area is judged to be a rock cover growth area; If the rock cover evolution coefficient is less than 0, the corresponding rocky desertification monitoring sub-area is judged to be a rock cover decline area; Obtaining surface vegetation monitoring data, and obtaining the vegetation coverage evolution coefficient corresponding to each rocky desertification monitoring sub-area based on the surface vegetation monitoring data; If the vegetation cover evolution coefficient is greater than 0, the corresponding rocky desertification monitoring sub-area is judged to be a vegetation cover growth area; If the vegetation cover evolution coefficient is less than or equal to 0, the corresponding rocky desertification monitoring sub-area is judged to be a vegetation cover decline area; Mark the rocky desertification monitoring sub-areas that are both in the vegetation cover growth area and the rock cover decline area as the rocky desertification control qualified area, and count the number of rocky desertification control qualified areas, and the number of control qualified areas; Obtain the number of rocky desertification monitoring sub-areas into which the rocky desertification monitoring area is divided, and obtain the number of monitoring sub-areas; Calculate the ratio of the number of qualified areas to the number of monitored sub-areas to obtain the qualified rate of treatment corresponding to the rocky desertification monitoring area; Obtain the preset governance pass rate. If the governance pass rate is greater than or equal to the preset governance pass rate, the rocky desertification governance in the rocky desertification monitoring area is evaluated to be qualified. If the governance pass rate is less than the preset governance pass rate, the rocky desertification governance in the rocky desertification monitoring area is evaluated to be unqualified.

[0031] It should be noted here that: In the present application, the preset control qualification rate involved here is 0.6. If the control qualification rate corresponding to the rocky desertification monitoring area is 0.56, then the corresponding rocky desertification monitoring area is evaluated as unqualified for rocky desertification control.

[0032] The preferred embodiments of the present invention disclosed above are only used to help explain the present invention. The preferred embodiments do not describe all the details in detail, nor do they limit the invention to only specific implementation methods. Obviously, many modifications and changes can be made according to the content of this specification. This specification selects and specifically describes these embodiments in order to better explain the principles and practical applications of the present invention, so that those skilled in the art can understand and use the present invention well. The present invention is limited only by the claims and their full scope and equivalents.

Claims

1. A rocky desertification monitoring data analysis method based on data analysis, characterized in that: include: Step S1: Acquire the rocky desertification monitoring area, divide the rocky desertification monitoring area into several rocky desertification monitoring sub-areas, conduct historical rock cover monitoring on each rocky desertification monitoring sub-area, and obtain the rock cover evolution coefficient corresponding to the rocky desertification monitoring sub-area to obtain surface rock monitoring data; Step S2: Perform historical vegetation coverage monitoring on each rocky desertification monitoring sub-area according to the surface vegetation monitoring data, and obtain the vegetation coverage evolution coefficient corresponding to each rocky desertification monitoring sub-area to obtain the surface vegetation monitoring data; Step S3: Conduct rocky desertification monitoring analysis on the rocky desertification monitoring area based on the surface rock monitoring data and the surface vegetation monitoring data, and conduct governance assessment on the rocky desertification monitoring area based on the analysis results.

2. The rocky desertification monitoring data analysis method based on data analysis according to claim 1 is characterized in that: The step S1 further includes the following specific steps: Step S11: acquiring a rocky desertification monitoring area, dividing the acquired rocky desertification monitoring area into a plurality of rocky desertification monitoring sub-areas of equal area, and selecting a sample rocky desertification sub-area from the acquired plurality of rocky desertification monitoring sub-areas; Step S12: Performing historical monitoring of rock cover in the sample rocky desertification sub-region to obtain a rock cover evolution coefficient corresponding to the sample rocky desertification sub-region; Step S13: acquiring the rock cover evolution coefficient corresponding to each rocky desertification monitoring sub-area to obtain a plurality of rock cover evolution coefficients; Step S14: defining a plurality of rocky desertification monitoring sub-areas and the rock cover evolution coefficient corresponding to each rocky desertification monitoring sub-area as surface rock monitoring data.

3. The rocky desertification monitoring data analysis method based on data analysis according to claim 2 is characterized in that: The step S12 further includes the following specific steps: step S121: In the process of historical monitoring of rock cover in the sample rocky desertification sub-area, mark the Y1 historical monitoring time point to the Ya historical monitoring time point; Step S122: monitoring the rock coverage of the sample rocky desertification sub-area at the Ya historical monitoring time point to obtain the rock coverage index value of the Ya area; Step S123: rock coverage monitoring is performed on the sample rocky desertification sub-area at the historical monitoring time point Y1 to the historical monitoring time point Ya-1, and the rock coverage index value of the Y1 area to the rock coverage index value of the Ya-1 area is obtained; Step S124: Mark the historical monitoring time point as the horizontal coordinate, and mark the regional rock coverage index value as the vertical coordinate to obtain the rock coverage coordinate system. In the rock coverage coordinate system, the Y1 historical monitoring time point and the Y1 regional rock coverage index value as the vertical coordinate coordinate point are obtained to obtain the Y1 rock coverage coordinate point. Similarly, the Ya historical monitoring time point and the Ya regional rock coverage index value as the vertical coordinate coordinate point are obtained to obtain the Ya rock coverage coordinate point. Step S125: in the rock cover coordinate system, obtain the slope of the line connecting the Y1 rock cover coordinate point and the Y2 rock cover coordinate point, obtain the slope value of the Y1 line, and obtain the slope value of the Ya-1 line; Step S126: sum the slope values ​​of the line connecting Y1 to the line connecting Ya-1 to obtain the rock cover evolution coefficient corresponding to the rocky desertification sub-area of ​​the sample.

4. The rocky desertification monitoring data analysis method based on data analysis according to claim 3 is characterized in that: The step S122 includes the following steps: Obtain the satellite remote sensing map of the sample rocky desertification sub-area at the Ya historical monitoring time point to obtain the Ya satellite remote sensing map; Several different types of cover rocks are extracted from the Ya satellite remote sensing map and named as the first type of cover rock to the mth type of cover rock; The coverage area values ​​corresponding to the first type of covered rock to the mth type of covered rock are obtained through the Ya satellite remote sensing map, and the coverage area values ​​of the first rock to the mth rock are obtained; The area value of the sample rocky desertification sub-area in the Ya satellite remote sensing map is obtained to obtain the area value of the rocky desertification area; Select multiple first-type covered rock samples, obtain the soil material content corresponding to each first-type covered rock sample, obtain multiple soil material contents, and average the obtained multiple soil material contents to obtain the first-type rock soil content.

5. The rocky desertification monitoring data analysis method based on data analysis according to claim 4 is characterized in that: The step S122 further includes the following specific steps: Obtaining the rock soil content of the second type of covering rock to the mth type of covering rock respectively, and obtaining the rock soil content of the second type of covering rock to the mth type of covering rock; The rock cover index value of the Ya region is obtained by calculating the rock desertification area value, the first rock cover area value to the mth rock cover area value, and the first type of rock soil content to the mth type of rock soil content; The rock cover index value of the Ya area is calculated, and the specific formula is as follows: ； Among them, Yzba is the rock coverage index value of the Ya area, Sfgi is the value of the i-th rock coverage area, Smz is the value of the rocky desertification area area, Nthi is the soil content of the i-th type of rock, and m is the number of types corresponding to the covering rock.

6. The rocky desertification monitoring data analysis method based on data analysis according to claim 1 is characterized in that: The step S2 further includes the following specific steps: Step S21: acquiring surface vegetation monitoring data, acquiring multiple rocky desertification monitoring sub-areas respectively according to the surface vegetation monitoring data, and selecting a characteristic rocky desertification sub-area from the acquired multiple rocky desertification monitoring sub-areas; Step S22: Performing historical monitoring of vegetation coverage in the characteristic rocky desertification sub-region to obtain a vegetation coverage evolution coefficient corresponding to the characteristic rocky desertification sub-region; Step S23: respectively obtaining the vegetation coverage evolution coefficient corresponding to each rocky desertification monitoring sub-area to obtain surface vegetation monitoring data.

7. The rocky desertification monitoring data analysis method based on data analysis according to claim 6 is characterized in that: The step S22 further includes the following specific steps: Step S221: in the process of historical monitoring of vegetation coverage in the characteristic rocky desertification sub-area, marking the historical monitoring time point Z1 to the historical monitoring time point Zb; Step S222: monitoring vegetation coverage of the characteristic rocky desertification sub-region at the historical monitoring time point Zb, and obtaining the vegetation coverage index value of the Zb region; Step S223: performing vegetation coverage monitoring on the characteristic rocky desertification sub-regions at the historical monitoring time point Z1 to the historical monitoring time point Zb-1, respectively, to obtain vegetation coverage index values ​​of the Z1 region to the Zb-1 region; Step S224: in the existing plane rectangular coordinate system, the historical monitoring time point is marked as the horizontal coordinate, and the regional vegetation coverage index value is marked as the vertical coordinate to obtain the vegetation coverage coordinate system. In the vegetation coverage coordinate system, the coordinate point with the vertical coordinate of the Z1 historical monitoring time point and the Z1 regional vegetation coverage index value is obtained to obtain the Z1 vegetation coverage coordinate. Similarly, the coordinate point with the vertical coordinate of the Zb historical monitoring time point and the Zb regional vegetation coverage index value is obtained to obtain the Zb vegetation coverage coordinate point. Step S225: in the vegetation coverage coordinate system, obtain the slope of the line connecting the Z1 vegetation coverage coordinate point and the Z2 vegetation coverage coordinate point to obtain the Z1 line slope value, and similarly obtain the slope of the line connecting the Zb-1 vegetation coverage coordinate point and the Zb vegetation coverage coordinate point to obtain the Zb-1 line slope value; Step S226: summing the slope values ​​of the lines connecting Z1 to Zb-1 to obtain the vegetation cover evolution coefficient corresponding to the characteristic rocky desertification sub-area.

8. The rocky desertification monitoring data analysis method based on data analysis according to claim 7 is characterized in that: The step S222 further includes the following specific steps: Obtain the satellite remote sensing map of the characteristic rocky desertification sub-area at the Zb historical monitoring time point to obtain the Zb satellite remote sensing map; Using image recognition algorithms, several different types of ground cover vegetation were extracted from the Zb satellite remote sensing map, and they were named the first type of cover vegetation to the nth type of cover vegetation. Obtain the coverage area values ​​corresponding to the first type of vegetation to the nth type of vegetation through the Zb satellite remote sensing map, and obtain the first vegetation coverage area value to the nth vegetation coverage area value; The area value of the characteristic rocky desertification sub-area in the Zb satellite remote sensing map is obtained to obtain the area value of the rocky desertification area; Selecting multiple first-type covered vegetation samples, obtaining the root moisture content corresponding to each first-type covered vegetation sample, obtaining multiple root moisture contents, and averaging the obtained multiple root moisture contents to obtain the root moisture content of the first-type vegetation; The root moisture content of vegetation is obtained for the second type of cover vegetation to the nth type of cover vegetation respectively, and the root moisture content of the second type of vegetation to the nth type of vegetation is obtained; The rocky desertification area value, the first vegetation coverage area value to the nth vegetation coverage area value and the first type of vegetation root moisture content to the nth type of vegetation root moisture content are calculated to obtain the Zb regional vegetation coverage index value; The vegetation coverage index value of the Zb area is calculated, and the specific formula is as follows: ； Among them, Zhsb is the vegetation coverage index value of Zb area, Zfgi is the value of the i-th vegetation coverage area, Smz is the value of the rocky desertification area area, Zhsi is the root moisture content of the i-th type of vegetation, and n is the number of types corresponding to the covered vegetation.

9. The rocky desertification monitoring data analysis method based on data analysis according to claim 1, characterized in that: The step S3 further includes the following specific steps: Obtaining surface rock monitoring data, and obtaining the rock cover evolution coefficient corresponding to each rocky desertification monitoring sub-area based on the surface rock monitoring data; If the rock cover evolution coefficient is greater than or equal to 0, the corresponding rocky desertification monitoring sub-area is judged to be a rock cover growth area; If the rock cover evolution coefficient is less than 0, the corresponding rocky desertification monitoring sub-area is judged to be a rock cover decline area; Obtaining surface vegetation monitoring data, and obtaining the vegetation coverage evolution coefficient corresponding to each rocky desertification monitoring sub-area based on the surface vegetation monitoring data; If the vegetation cover evolution coefficient is greater than 0, the corresponding rocky desertification monitoring sub-area is judged to be a vegetation cover growth area; If the vegetation cover evolution coefficient is less than or equal to 0, the corresponding rocky desertification monitoring sub-area is judged to be a vegetation cover decline area; Mark the rocky desertification monitoring sub-areas that are both in the vegetation cover growth area and the rock cover decline area as the rocky desertification control qualified area, and count the number of rocky desertification control qualified areas, and the number of control qualified areas; Obtain the number of rocky desertification monitoring sub-areas into which the rocky desertification monitoring area is divided, and obtain the number of monitoring sub-areas; Calculate the ratio of the number of qualified areas to the number of monitored sub-areas to obtain the qualified rate of treatment corresponding to the rocky desertification monitoring area; Obtain the preset governance pass rate. If the governance pass rate is greater than or equal to the preset governance pass rate, the rocky desertification governance in the rocky desertification monitoring area is evaluated to be qualified. If the governance pass rate is less than the preset governance pass rate, the rocky desertification governance in the rocky desertification monitoring area is evaluated to be unqualified.

10. A rocky desertification monitoring data analysis system based on data analysis, applicable to the rocky desertification monitoring data analysis method based on data analysis according to any one of claims 1 to 9, characterized in that: The data analysis system comprises: Rock monitoring module: Obtain the rocky desertification monitoring area and divide it into several rocky desertification monitoring sub-areas, conduct historical rock cover monitoring on each rocky desertification monitoring sub-area, obtain the rock cover evolution coefficient corresponding to the rocky desertification monitoring sub-area, and obtain the surface rock monitoring data; Vegetation monitoring module: Conduct historical vegetation coverage monitoring for each desertification monitoring sub-area based on surface vegetation monitoring data, obtain the vegetation coverage evolution coefficient corresponding to each desertification monitoring sub-area, and obtain surface vegetation monitoring data; Monitoring and analysis module: Conduct rocky desertification monitoring and analysis on the rocky desertification monitoring area based on the surface rock monitoring data and surface vegetation monitoring data, and conduct a governance qualification assessment on the rocky desertification monitoring area based on the analysis results.

Images