Dry-hot valley basin vegetation engineering construction priority area construction method based on microclimate

Abstract

The invention relates to the technical field of region identification, and discloses a method for constructing a vegetation engineering construction priority region of a dry-hot valley basin based on microclimate. Comprising the following steps: drawing a construction region and a region boundary, dividing the construction region into sub-regions, screening out effective sub-regions from the sub-regions, obtaining comprehensive micro-meteorological parameters of the effective sub-regions, calculating a construction evaluation value, and constructing a construction priority region from the effective sub-regions; compared with the prior art, the method has the advantages that the sub-regions are drawn according to the region boundaries, the effect of breaking up a large-area region into parts in the dry-hot valley drainage basin is achieved, meanwhile, the comprehensive micro-meteorological parameters of the effective sub-regions are collected, refined collection and analysis of the micro-meteorological parameters of the small regions in the dry-hot valley drainage basin are achieved, and the analysis accuracy is improved. And on the basis of the construction priority value, the construction priority area meeting the vegetation engineering construction requirement is identified, an accurate and appropriate position is provided for subsequent vegetation engineering construction, and the aim of treating the debris flow source of the dry and hot valley basin is fulfilled.

Description

Technical Field

[0001] The present invention relates to the technical field of regional identification, and more specifically, to a method for constructing a priority area for vegetation engineering construction in a dry-hot river valley basin based on micrometeorology. Background Art

[0002] The dry and hot river valley basins have both low and medium altitude areas with extremely high temperatures and droughts that hinder the growth of vegetation, and high altitude areas with insufficient accumulated temperature. At the same time, extreme meteorological conditions allow gullies to fully develop, providing a rich source of material for the formation of debris flows. Vegetation engineering is an effective means of debris flow prevention and control, especially on slopes with fully developed gullies. Implementing vegetation engineering can meet the needs of debris flow source control. Therefore, it is necessary to accurately identify the vegetation engineering areas in the dry and hot river valley basins.

[0003] The patent application with reference publication number CN116680658A discloses a method and system for selecting a site for a heat wave monitoring station based on risk assessment. By giving priority to high-risk areas and taking into account the spatial characteristics of meteorological elements, it can quickly screen out candidate site locations and achieve the same effect as the traditional method with fewer heat wave monitoring stations, thereby significantly saving installation and maintenance costs. It has a wide range of adaptability and can take into account future scenarios. It is not only applicable to the whole country, but can also comprehensively consider future meteorological data and propose site selection plans under climate adaptation changes;

[0004] The prior art has the following deficiencies:

[0005] The existing construction priority areas are mainly evaluated for the suitability of vegetation engineering construction through comprehensive remote sensing technology, geographic information technology, combined with geographic statistical methods and ecological models, and then the priority level of vegetation engineering construction is determined according to the goals, and finally the construction priority areas are planned. However, in the dry and hot river valley basins, due to the lack of in-depth understanding of local hydrothermal processes, there is a large uncertainty in meteorological, hydrological and other spatial data, resulting in the inability to accurately collect micro-meteorological parameters in small areas, thereby reducing the calculation accuracy of regional priority values, resulting in low accuracy of the construction priority areas identified in the dry and hot river valley basins, which cannot meet the needs of debris flow source control in the dry and hot river valley basins.

[0006] In view of this, the present invention proposes a method for constructing priority areas for vegetation engineering construction in dry and hot river valley basins based on micrometeorology to solve the above problems. Summary of the invention

[0007] In order to overcome the above-mentioned defects of the prior art and to achieve the above-mentioned purpose, the present invention provides the following technical solution: a method for constructing a priority area for vegetation engineering construction in a dry and hot river valley basin based on micrometeorology, comprising:

[0008] S1: Mark the initial area of ​​the selected region on the electronic map, mark the regional critical points on the initial area, and draw the construction area;

[0009] S2: Mark the regional boundary lines of the construction area and divide the construction area into sub-areas based on the regional boundary lines. The regional boundary lines include slope boundary lines, ridge boundary lines and river bank boundary lines;

[0010] S3: Obtain basic geographic parameters of the sub-region, calculate the regional stability value, and screen out valid sub-regions from the sub-regions. The basic geographic parameters include altitude difference value and terrain flatness value;

[0011] S4: Obtain the comprehensive micro-meteorological parameters of the effective sub-area, and calculate the construction evaluation value based on the comprehensive micro-meteorological parameters, the comprehensive micro-meteorological parameters include the vertical temperature decrease rate, the unit precipitation difference value and the effective light duration value;

[0012] S5: Obtain the construction cost value of the effective sub-area, combine the construction assessment value and the construction cost value to generate a construction priority value, and construct a construction priority area from the effective sub-area.

[0013] Furthermore, the marking method of the initial area includes:

[0014] Query the maximum latitude, minimum latitude, maximum longitude and minimum longitude of the selected area through the geographic information system;

[0015] Mark the points corresponding to the maximum latitude, the minimum latitude, the maximum longitude, and the minimum longitude on the electronic map to obtain a first point, a second point, a third point, and a fourth point;

[0016] Draw a line horizontally through the positions of the first point and the second point to obtain two latitude lines, and draw a line vertically through the positions of the third point and the fourth point to obtain two longitude lines;

[0017] Extend the two latitude lines and the two longitude lines respectively until they intersect, and record the area enclosed by the two latitude lines and the two longitude lines as the initial area.

[0018] Furthermore, the method for drawing the construction area includes:

[0019] Based on the preset interval length, i vertical extension lines and i horizontal extension lines extending into the initial area are drawn on two latitude lines and two longitude lines respectively;

[0020] Mark p nodes equidistantly on i vertical extension lines and i horizontal extension lines, and query the altitude values ​​of the p nodes one by one through the geographic information system;

[0021] Along the drawing direction of the vertical extension line and the horizontal extension line, the first node whose elevation value is greater than the elevation values ​​of the two adjacent nodes is recorded as a valid node;

[0022] After summarizing the first point, the second point, the third point, the fourth point and all valid nodes, k regional critical points are obtained;

[0023] Taking the regional critical point corresponding to the first point as the starting point of the connection, k regional critical points are connected in a clockwise direction to draw the construction area.

[0024] Furthermore, the sub-region division method includes:

[0025] The elevation values ​​of all locations in the construction area are queried through the geographic information system, and a digital elevation model of the construction area is constructed in combination with virtual reality technology;

[0026] Based on the preset side length, the digital elevation model is divided into w cells. The center points of the cells are identified through computer vision technology, and the altitude values ​​of the center points are marked one by one.

[0027] Record the center point whose altitude value is greater than the preset first altitude threshold as the ridge center point, and record the cell where the ridge center point is located as the ridge cell;

[0028] The center point whose altitude value is less than the preset second altitude threshold is recorded as the river bank center point, and the cell where the river bank center point is located is recorded as the river bank cell;

[0029] Recording a center point whose altitude value is between a preset first altitude threshold and a preset second altitude threshold as a slope center point, and recording a cell where the slope center point is located as a slope cell;

[0030] The outer boundary of the cell is recorded as the target boundary, and lines are drawn along the target boundaries of the ridge cell, the river bank cell and the slope cell respectively to obtain the ridge dividing line, the river bank dividing line and the slope dividing line;

[0031] Based on the ridge boundary line, river bank boundary line and slope boundary line, the construction area is divided into m closed sub-areas.

[0032] Furthermore, the method for obtaining the altitude difference value includes:

[0033] Scan the m sub-regions one by one to obtain m regional scan images, and mark all the pixels in the m regional scan images to obtain n pixels;

[0034] The altitude values ​​of n pixels are queried one by one through the geographic information system, and the maximum altitude value of the pixel point and the minimum altitude value of the pixel point are screened out. The maximum altitude value of m pixels and the minimum altitude value of m pixels are subtracted to obtain m altitude difference values;

[0035] The expression for the altitude difference is:

[0036] HB lcm =HB zdm -HB zxm ;

[0037] In the formula, HB lcm is the elevation difference of the mth sub-area, HB zdm is the maximum altitude value of the pixel points in the mth sub-area, HB zxm is the minimum altitude value of the pixels in the mth sub-region.

[0038] Furthermore, the method for obtaining the terrain flatness value includes:

[0039] Randomly select v pixels as the center of the circle in the m area scans, draw a circle with a preset length as the radius, and draw V pixel circles;

[0040] Count the altitude values ​​of x pixels within the v pixel circle one by one, and add up the altitude values ​​of x pixels and calculate the average to obtain the altitude values ​​of v circles;

[0041] The expression for the circle altitude value is:

[0042]

[0043] In the formula, HB zymv is the circle height value of the vth pixel circle in the mth sub-region, HB xsmva is the altitude value of the a-th pixel point of the v-th pixel circle of the m-th sub-region;

[0044] After subtracting the vth circle elevation value from the v-1th circle elevation value, we get v-1 sub-flat values;

[0045] The expression of the sub-flat value is:

[0046] PT zzmv-1 =HB zymv -HB zymv-1 ;

[0047] Where, PT zzmv-1 is the v-1th sub-flat value of the mth sub-region, HB zymv-1 is the altitude value of the v-1th circle in the mth sub-area;

[0048] After removing the maximum and minimum values ​​of the sub-flat values, the remaining v-3 sub-flat values ​​are accumulated and averaged to obtain m terrain flatness values;

[0049] The expression for the flat terrain value is:

[0050]

[0051] Where, PT dsm is the terrain flatness value of the mth sub-area, PT zzmb is the bth sub-flat value of the mth sub-region.

[0052] Furthermore, the screening method of the effective sub-region includes:

[0053] Combine the m altitude drop values ​​and the m terrain flatness values ​​to obtain m regional stability values;

[0054] The expression of the regional stationary value is:

[0055]

[0056] In the formula, QY pwm is the regional stability value of the mth sub-region, σ1 and σ2 are proportional coefficients greater than 0;

[0057] The m regional stationary values ​​QY pwm Sequentially with the preset regional stability threshold QY yz Compare;

[0058] When QY pwm Less than QY yz When , the mth sub-region is not recognized as a valid sub-region;

[0059] When QY pwm Greater than or equal to QY yz , the mth sub-region is identified as a valid sub-region, and Z valid sub-regions are obtained.

[0060] Furthermore, the method for obtaining the vertical temperature decrease rate includes:

[0061] In the Z effective sub-areas, q temperature measurement locations with different altitudes are randomly marked to obtain q temperature measurement points;

[0062] Taking the same time as the starting time and the preset temperature measurement time as the interval standard, r temperature measurement moments are marked, and the temperature of the r temperature measurement moments is detected by the temperature sensor respectively to obtain r temperature values, and the r temperature values ​​are accumulated and averaged to obtain q sub-temperature measurement values;

[0063] The expression of the sub-temperature measurement value is:

[0064]

[0065] In the formula, CW zzq is the sub-temperature measurement value of the qth temperature measurement point in the zth effective sub-area, WD zzqc is the cth temperature value of the qth temperature measurement point in the zth effective sub-area;

[0066] The altitude values ​​of q temperature measurement points are queried one by one through the geographic information system to obtain q temperature measurement altitudes, and the qth sub-temperature measurement value is subtracted from the q-1th sub-temperature measurement value, and then compared with the difference between the qth temperature measurement altitude and the q-1th temperature measurement altitude to obtain the q-1 sub-decline rate;

[0067] The expression of the sub-decline rate is:

[0068]

[0069] In the formula, DJ lzq-1 is the q-1th sub-decline rate of the zth effective sub-region, CW zzq-1 is the sub-temperature measurement value of the q-1th temperature measurement point in the zth effective sub-area, HB cwzq is the qth temperature measurement altitude of the zth effective sub-area, HB cwzq-1 is the q-1th temperature measurement altitude of the zth effective sub-area;

[0070] The q-1 sub-decline rates are accumulated and averaged to obtain z vertical temperature lapse rates;

[0071] The expression of temperature vertical lapse rate is:

[0072]

[0073] In the formula, DJ wdz is the vertical temperature decrease rate of the zth effective sub-region, DJ lzd is the dth sub-decline rate of the zth effective sub-region.

[0074] Furthermore, the method for obtaining the unit precipitation difference value includes:

[0075] Taking the preset rainfall cycle as the standard, the precipitation data of q temperature measuring points are detected by the rainfall sensor to obtain q precipitation values, and the q precipitation values ​​are respectively compared with the duration corresponding to the preset rainfall cycle to obtain q unit precipitation;

[0076] The expression of unit precipitation is:

[0077]

[0078] In the formula, JS dwzq is the qth unit precipitation in the zth effective sub-area, JS zzqis the qth precipitation value of the zth effective sub-area, SC yl The duration corresponding to the preset rainfall cycle;

[0079] Subtract the q unit precipitation from the unit precipitation threshold in turn, and then add up the differences and average them to obtain the z unit precipitation difference values;

[0080] The expression of unit precipitation difference value is:

[0081]

[0082] In the formula, JS cez is the unit precipitation difference value of the zth effective sub-area, JS dwze The e-th unit precipitation in the z-th effective sub-area, JS yz is the unit precipitation threshold;

[0083] The expression to construct the evaluation value is:

[0084]

[0085] In the formula, GJ pgz is the construction evaluation value of the zth valid sub-region, GZ scz is the effective illumination duration value of the zth effective sub-area, and σ3, σ4, and σ5 are proportional coefficients greater than 0.

[0086] 11. Further, construct the expression of priority value:

[0087]

[0088] In the formula, JS yxz is the construction priority value of the zth effective sub-area, JS cbz is the construction cost value of the zth effective sub-area, σ6 and σ7 are proportional coefficients greater than 0;

[0089] Methods for constructing construction priority areas include:

[0090] Arrange the Z construction priority values ​​in descending order from largest to smallest, and mark the first, second, third, and fourth construction priority values;

[0091] The valid sub-regions corresponding to the first, second, third and fourth construction priority values ​​are recorded as the first region, the second region, the third region and the fourth region respectively;

[0092] The center positions of the first area, the second area, the third area, and the fourth area are identified one by one by computer vision technology to obtain a first center position, a second center position, a third center position, and a fourth center position;

[0093] Taking the first center position, the second center position, the third center position and the fourth center position as expansion base points respectively, expanding outwards on the first area, the second area, the third area and the fourth area to obtain a first expansion area, a second expansion area, a third expansion area and a fourth expansion area;

[0094] Stop expanding until the areas of the first expansion area, the second expansion area, the third expansion area and the fourth expansion area reach 70% of the area of ​​the first area, 40% of the area of ​​the second area, 20% of the area of ​​the third area and 10% of the area of ​​the fourth area respectively;

[0095] The first expansion area, the second expansion area, the third expansion area and the fourth expansion area after the expansion is stopped are recorded as construction priority areas.

[0096] The technical effects and advantages of the method for constructing priority areas for vegetation engineering construction in dry and hot river valley basins based on micrometeorology of the present invention are as follows:

[0097] The present invention marks an initial area of ​​a selected area on an electronic map, marks a regional critical point on the initial area, draws a construction area, draws a regional boundary line of the construction area, and divides the construction area into sub-areas based on the regional boundary line, obtains basic geographical parameters of the sub-areas, calculates regional stability values ​​of the sub-areas, and screens out effective sub-areas from the sub-areas, obtains comprehensive micro-meteorological parameters of the effective sub-areas, and calculates construction evaluation values ​​of the effective sub-areas based on the comprehensive micro-meteorological parameters, inquires about the construction cost value of the effective sub-area, combines the construction evaluation value and the construction cost value, generates a construction priority value, and constructs a construction priority area from the effective sub-area; compared with the prior art, by drawing a regional boundary line in the initial area, the boundaries of different geographical environment areas in a dry and hot river valley basin can be determined. Marking, and drawing sub-regions according to regional boundaries, to achieve the effect of breaking up larger areas in the dry and hot valley basin into smaller areas, which is convenient for the accurate identification of construction priority areas. At the same time, the regional stability value is calculated in combination with the basic geographical parameters, and the effective sub-regions that meet the requirements of vegetation engineering can be screened out from the sub-regions. By collecting the comprehensive micro-meteorological parameters of the effective sub-regions, the detailed collection and analysis of the micro-meteorological parameters of small areas in the dry and hot valley basin is realized, and the construction priority values ​​of the effective sub-regions can be calculated in different spatial dimensions, and then the construction priority areas that meet the needs of vegetation engineering construction can be identified, and accurate and appropriate regional locations can be provided for the subsequent construction of vegetation engineering in the dry and hot valley basin, ensuring that the vegetation in the construction priority areas can grow healthily and rapidly, meeting the ultimate goal of mudslide source control in the dry and hot valley basin. BRIEF DESCRIPTION OF THE DRAWINGS

[0098] Figure 1A schematic diagram of a process for constructing a priority area for vegetation engineering construction in a dry-hot river valley basin based on micrometeorology provided in the first embodiment of the present invention;

[0099] Figure 2 This is a schematic diagram of a system for constructing priority areas for vegetation engineering construction in dry-hot river valleys based on micrometeorology provided in Example 2 of the present invention. DETAILED DESCRIPTION

[0100] The following will be combined with the drawings in the embodiments of the present invention to clearly and completely describe the technical solutions in the embodiments of the present invention. Obviously, the described embodiments are only part of the embodiments of the present invention, not all of the embodiments. Based on the embodiments of the present invention, all other embodiments obtained by ordinary technicians in this field without creative work are within the scope of protection of the present invention.

[0101] Example 1: Please refer to Figure 1 As shown, the method for constructing a priority area for vegetation engineering construction in a dry-hot river valley basin based on micrometeorology in this embodiment includes:

[0102] S1: Mark the initial area of ​​the selected region on the electronic map, mark the regional critical points on the initial area, and draw the construction area;

[0103] The selected area refers to the area where the dry and hot river valley basin where vegetation engineering construction is required is located, which can be used as the basis for drawing the construction area on the electronic map, while the initial area is the area used to preliminarily represent the location of the selected area on the electronic map;

[0104] The initial area marking methods include:

[0105] Query the maximum latitude, minimum latitude, maximum longitude and minimum longitude of the selected area through the geographic information system;

[0106] The points corresponding to the maximum latitude, the minimum latitude, the maximum longitude and the minimum longitude are marked on the electronic map to obtain the first point, the second point, the third point and the fourth point; by selecting the points corresponding to the maximum latitude, the minimum latitude, the maximum longitude and the minimum longitude, the positions and directions of the four boundaries of the initial area can be preliminarily located, and point support can be provided for the subsequent determination of the boundaries of the initial area;

[0107] Draw a line horizontally through the positions of the first point and the second point to obtain two latitude lines, and draw a line vertically through the positions of the third point and the fourth point to obtain two longitude lines;

[0108] Extend the two latitude lines and the two longitude lines respectively until they intersect, and record the area enclosed by the two latitude lines and the two longitude lines as the initial area.

[0109] Since the area represented by the initial area is irregular, and the boundary of the initial area is also blurred due to interference from factors such as rivers and vegetation, in order to accurately identify the boundary of the initial area, it is necessary to mark the regional critical points on the boundary of the initial area, and draw the construction area based on the regional critical points;

[0110] Methods for drawing construction areas include:

[0111] Based on the preset interval length, i vertical extension lines and i horizontal extension lines extending into the initial area are drawn on two latitude lines and two longitude lines respectively; the preset interval length is used to represent the interval distance between two adjacent extension lines, which can ensure that the two adjacent extension lines are independent of each other, prevent the two adjacent extension lines from overlapping, and thus avoid the phenomenon of data duplication caused by the overlap of subsequent nodes;

[0112] Mark p nodes equidistantly on i vertical extension lines and i horizontal extension lines, and query the altitude values ​​of the p nodes one by one through the geographic information system;

[0113] Along the drawing direction of the vertical extension line and the horizontal extension line, the first node whose elevation value is greater than the elevation values ​​of the two adjacent nodes is recorded as a valid node; at this time, the elevation values ​​of the nodes before and after the valid node will be smaller than the elevation value of the valid node, so that the valid node can be used as the critical point between the construction area and the non-construction area;

[0114] After summarizing the first point, the second point, the third point, the fourth point and all valid nodes, k regional critical points are obtained;

[0115] Taking the regional critical point corresponding to the first point as the starting point of the connection, k regional critical points are connected in a clockwise direction to draw the construction area.

[0116] It should be noted that the construction area is a real and accurate location area representation of the selected area on the electronic map. At this time, the construction area corresponds to a larger area and is used to comprehensively represent the overall situation of the selected area.

[0117] S2: Mark the regional boundary of the construction area, and divide the construction area into sub-areas based on the regional boundary;

[0118] The regional boundary line refers to the boundary line that can divide the locations of different geographic attribute information in the construction area, so that the different geographic attribute information contained in the construction area can be effectively distinguished, and based on the regional boundary line, the overall and larger construction area is divided into smaller sub-areas, so that the type of geographic attribute information contained in each sub-area is relatively single and the amount is relatively small;

[0119] Regional boundaries include slope boundaries, ridge boundaries and river bank boundaries; slope boundaries refer to boundaries that can divide the mountain slopes within the construction area, ridge boundaries refer to boundaries that can divide the mountain ridges within the construction area, and river bank boundaries refer to boundaries that can divide the river banks within the construction area;

[0120] After the regional boundary line is drawn, the construction area can be divided into multiple relatively closed sub-areas based on the regional boundary line, so that each sub-area can provide a detailed representation of the selected area;

[0121] The sub-area division methods include:

[0122] The elevation values ​​of all locations in the construction area are queried through the geographic information system, and a digital elevation model of the construction area is constructed in combination with virtual reality technology; both virtual reality technology and digital elevation model are existing technologies that can construct specific regional location information into a model with precise coordinates, thereby facilitating the precise identification of subsequent demarcation lines;

[0123] Based on the preset side length, the digital elevation model is divided into w cells. The center point of each cell is identified by computer vision technology, and the altitude value of each center point is marked one by one. The preset side length is used to represent the side length of each cell, which ensures that the side length of each cell is the same, thus ensuring that the digital elevation model can be evenly divided.

[0124] The center point whose altitude value is greater than the preset first altitude threshold is recorded as the ridge center point, and the cell where the ridge center point is located is recorded as the ridge cell; the preset first altitude threshold is used to represent the minimum value of the altitude value identified as the ridge center point, so that the ridge center point can be accurately identified;

[0125] The center point whose altitude value is less than the preset second altitude threshold is recorded as the river bank center point, and the cell where the river bank center point is located is recorded as the river bank cell; the preset second altitude threshold is used to represent the maximum value of the altitude value identified as the river bank center point, so that the river bank center point can be accurately identified; the preset first altitude threshold and the preset second altitude threshold are obtained by respectively calculating the average value of the altitude values ​​corresponding to the minimum value of the ridge center point and the maximum value of the river bank center point in actual large numbers;

[0126] Recording a center point whose altitude value is between a preset first altitude threshold and a preset second altitude threshold as a slope center point, and recording a cell where the slope center point is located as a slope cell;

[0127] The outer boundary of the cell is recorded as the target boundary, and lines are drawn along the target boundaries of the ridge cell, the river bank cell and the slope cell to obtain the ridge boundary line, the river bank boundary line and the slope boundary line;

[0128] Based on the ridge boundary line, river bank boundary line and slope boundary line, the construction area is divided into m closed sub-areas.

[0129] It should be noted that when divided into sub-regions, the actual locations of the selected areas corresponding to two adjacent sub-regions are adjacent, so that the sub-regions can fully and continuously represent the location of the construction area, avoiding the omission of a certain area in the dry and hot valley basin.

[0130] S3: Obtain the basic geographic parameters of the sub-region, calculate the regional stability value, and screen out the effective sub-regions from the sub-regions;

[0131] Basic geographic parameters refer to data that can represent the stability of the geographical environment where the sub-region is located, and can serve as the data basis for the subsequent calculation of regional stability values, so that the calculated regional stability values ​​can preliminarily identify the stability of the geographical environment where each sub-region is located;

[0132] The basic geographic parameters include the altitude difference value and the terrain flatness value; the altitude difference value is used to numerically express the difference between the maximum altitude value and the minimum altitude value of the sub-region, that is, to express the height difference of the sub-region; the terrain flatness value is used to numerically express the flatness of the ground at the sub-region, that is, to express the flatness of the sub-region;

[0133] Methods for obtaining altitude difference values ​​include:

[0134] Scan the m sub-regions one by one to obtain m regional scan images, and mark all the pixels in the m regional scan images to obtain n pixels;

[0135] The altitude values ​​of n pixels are queried one by one through the geographic information system, and the maximum altitude value of the pixel point and the minimum altitude value of the pixel point are screened out. The maximum altitude value of m pixels and the minimum altitude value of m pixels are subtracted to obtain m altitude difference values;

[0136] The expression for the altitude difference is:

[0137] HB lcm =HB zdm -HB zxm ;

[0138] In the formula, HB lcm is the elevation difference of the mth sub-area, HB zdm is the maximum altitude value of the pixel points in the mth sub-area, HB zxm is the minimum altitude value of the pixels in the mth sub-region.

[0139] Methods for obtaining the terrain flatness value include:

[0140] Randomly select V pixel points as the center of the circle in the m area scan images, draw a circle with a preset length as the radius, and draw V pixel circles; the preset length is used to represent the radius of the pixel circle, so as to ensure that the drawn pixel circles can maintain the same size, and also ensure that each pixel circle can contain a sufficient number of pixels, so as to facilitate the rapid and accurate calculation of the circle altitude value;

[0141] Count the altitude values ​​of x pixels within the v pixel circle one by one, and add up the altitude values ​​of x pixels and calculate the average to obtain the altitude values ​​of v circles;

[0142] The expression for the circle altitude value is:

[0143]

[0144] In the formula, HB zvmv is the circle height value of the vth pixel circle in the mth sub-region, HB xsmva is the altitude value of the a-th pixel point of the v-th pixel circle of the m-th sub-region;

[0145] After subtracting the vth circle elevation value from the v-1th circle elevation value, we get v-1 sub-flat values;

[0146] The expression of the sub-flat value is:

[0147] PT zzmv-1 =HB zymv -HB zymv-1 ;

[0148] Where, PT zzmv-1 is the v-1th sub-flat value of the mth sub-region, HB zymv-1 is the altitude value of the v-1th circle in the mth sub-area;

[0149] After removing the maximum and minimum values ​​of the sub-flat values, the remaining v-3 sub-flat values ​​are accumulated and averaged to obtain m terrain flatness values;

[0150] The expression for the flat terrain value is:

[0151]

[0152] Where, PT dsm is the terrain flatness value of the mth sub-area, PT zzmb is the bth sub-flat value of the mth sub-region.

[0153] After obtaining the altitude drop value and the terrain flatness value, the regional stability value can be calculated according to the altitude drop value and the terrain flatness value, and by comparing the regional stability values, it can be judged whether the sub-region can be used as a valid sub-region, ensuring that the corresponding geographical area of ​​the selected valid sub-region has good stability and can be used as the direct screening range of the subsequent construction priority area;

[0154] The effective sub-area screening methods include:

[0155] Combine the m altitude drop values ​​and the m terrain flatness values ​​to obtain m regional stability values;

[0156] The expression of the regional stationary value is:

[0157]

[0158] In the formula, QY pwm is the regional stability value of the mth sub-region, σ1 and σ2 are proportional coefficients greater than 0;

[0159] Among them, σ1+σ2=1, for example, σ1 is 0.43, and σ2 is 57; σ1 and σ2 are to balance the proportion of the altitude difference value and the terrain flatness value in the regional stability value, so that the influence of the balanced altitude difference value and the terrain flatness value on the size of the regional stability value can be relatively balanced, thereby ensuring the accuracy of the regional stability value;

[0160] The m regional stationary values ​​QY pwm Sequentially with the preset regional stability threshold QY yz Comparison; The preset regional stability threshold is the critical value of the regional stability value when it is identified as a valid sub-region and when it is not identified as a valid sub-region, which can accurately distinguish whether the sub-region is identified as a valid sub-region, thereby ensuring that the identification and construction of the subsequent construction priority area remain accurate;

[0161] When QY pwm Less than QY yz When , it means that the stability value of the m-th region is less than the preset regional stability threshold. At this time, the stability of the geographical environment of the m-th sub-region is poor, and the m-th sub-region is not identified as a valid sub-region;

[0162] When QY pwm Greater than or equal to QYyz When , it means that the stability value of the m-th region is greater than or equal to the preset regional stability threshold. At this time, the stability of the geographical environment of the m-th sub-region is good, then the m-th sub-region is identified as a valid sub-region, and Z valid sub-regions are obtained.

[0163] S4: Obtain comprehensive micro-meteorological parameters of the effective sub-area, and calculate the construction evaluation value based on the comprehensive micro-meteorological parameters;

[0164] Comprehensive micro-meteorological parameters refer to parameters that can comprehensively represent the changes in micro-meteorology at the location of the effective sub-region, and serve as the basis for evaluating the effective sub-region as a construction priority area, thereby improving the calculation accuracy and comprehensiveness of the construction evaluation value;

[0165] Comprehensive micro-meteorological parameters include vertical temperature lapse rate, unit precipitation difference value and effective sunshine duration value;

[0166] The vertical temperature lapse rate refers to the decrease in the temperature at different heights of the effective sub-region as the height changes, which can be used to represent the temperature change in the effective sub-region. The larger the vertical temperature lapse rate, the larger the decrease in the temperature at different heights of the effective sub-region as the height changes, and the smaller the constructed evaluation value.

[0167] The method for obtaining the vertical temperature decrease rate includes:

[0168] In the Z effective sub-areas, q temperature measurement locations with different altitudes are randomly marked to obtain q temperature measurement points;

[0169] Taking the same moment as the starting moment and the preset temperature measurement time as the interval standard, r temperature measurement moments are marked, and the temperature of the r temperature measurement moments is detected by the temperature sensor respectively to obtain r temperature values, and the r temperature values ​​are accumulated and averaged to obtain q sub-temperature measurement values; the preset temperature measurement time is used to represent the numerical value of the temperature collection moment corresponding to the temperature measurement points at different heights in the effective sub-area, ensuring that the two collected temperature values ​​can be separated by a longer time, so that the two adjacent temperature values ​​of the temperature measurement points will change;

[0170] The expression of the sub-temperature measurement value is:

[0171]

[0172] In the formula, CW zzq is the sub-temperature measurement value of the qth temperature measurement point in the zth effective sub-area, WD zzqc is the cth temperature value of the qth temperature measurement point in the zth effective sub-area;

[0173] The altitude values ​​of q temperature measurement points are queried one by one through the geographic information system to obtain q temperature measurement altitudes, and the qth sub-temperature measurement value is subtracted from the q-1th sub-temperature measurement value, and then compared with the difference between the qth temperature measurement altitude and the q-1th temperature measurement altitude to obtain the q-1 sub-decline rate;

[0174] The expression of the sub-decline rate is:

[0175]

[0176] In the formula, DJ lzq-1 is the q-1th sub-decline rate of the zth effective sub-region, CW zzq-1 is the sub-temperature measurement value of the q-1th temperature measurement point in the zth effective sub-area, HB cwzq is the qth temperature measurement altitude of the zth effective sub-area, HB cwzq-1 is the q-1th temperature measurement altitude of the zth effective sub-area;

[0177] The q-1 sub-decline rates are accumulated and averaged to obtain z vertical temperature lapse rates;

[0178] The expression of temperature vertical lapse rate is:

[0179]

[0180] In the formula, DJ wdz is the vertical temperature decrease rate of the zth effective sub-region, DJ lzd is the dth sub-decline rate of the zth effective sub-region.

[0181] The unit precipitation difference value refers to the difference between the precipitation in the effective sub-area per unit time and the standard precipitation, which can represent the precipitation change in the effective sub-area. The larger the unit precipitation difference value, the larger the difference between the precipitation in the effective sub-area per unit time and the standard precipitation, and the smaller the constructed evaluation value.

[0182] The methods for obtaining the unit precipitation difference value include:

[0183] Taking the preset rainfall cycle as the standard, the precipitation data of q temperature measuring points are detected by the rainfall sensor to obtain q precipitation values, and the q precipitation values ​​are respectively compared with the duration corresponding to the preset rainfall cycle to obtain q unit precipitation. The preset rainfall cycle is used to numerically represent the interval between two adjacent precipitation values, so as to ensure that two adjacent precipitation values ​​can have sufficient precipitation when being collected, and at the same time reduce the amount of collected data of the precipitation value, thereby ensuring the calculation accuracy of the unit precipitation.

[0184] The expression of unit precipitation is:

[0185]

[0186] In the formula, JS dwzq is the qth unit precipitation in the zth effective sub-area, JS zzq is the qth precipitation value of the zth effective sub-area, SC yl The duration corresponding to the preset rainfall cycle;

[0187] Subtract the q unit precipitation from the unit precipitation threshold in turn, and then add up the differences and average them to obtain Z unit precipitation difference values; the unit precipitation threshold refers to the unit precipitation in the effective sub-area per unit time under normal circumstances, which can be used as a basis for judging whether the unit precipitation exceeds the limit;

[0188] The expression of unit precipitation difference value is:

[0189]

[0190] In the formula, JS cez is the unit precipitation difference value of the zth effective sub-area, JS dwze The e-th unit precipitation in the z-th effective sub-area, JS yz is the unit precipitation threshold.

[0191] The effective illumination duration value refers to the duration of time that the real-time illumination intensity of the effective sub-area is higher than the standard illumination intensity, which can represent the illumination change in the effective sub-area. The larger the effective illumination duration value, the longer the real-time illumination intensity of the effective sub-area is higher than the standard illumination intensity, and the larger the constructed evaluation value. The effective illumination duration value is obtained through detection by the illumination sensors deployed in the effective sub-area.

[0192] After obtaining the vertical temperature decrease rate, unit precipitation difference value and effective light duration value, the construction evaluation value can be calculated so that the construction evaluation value can numerically represent the order of construction priority areas for each effective sub-area;

[0193] The expression to construct the evaluation value is:

[0194]

[0195] In the formula, GJ pgz is the construction evaluation value of the zth valid sub-region, GZ scz is the effective illumination duration value of the zth effective sub-area, σ3, σ4, and σ5 are proportional coefficients greater than 0;

[0196] Among them, the setting logic of σ3, σ4, and σ5 is consistent with the setting logic of σ1 and σ2 mentioned above.

[0197] S5: Obtain the construction cost value of the effective sub-area, combine the construction assessment value and the construction cost value to generate a construction priority value, and construct a construction priority area from the effective sub-area;

[0198] The construction cost value refers to the data stored in the database that has the same vertical temperature decline rate, unit precipitation difference value and effective light duration value as the Z effective sub-areas, and is used to simulate and calculate the construction cost required for the construction of vegetation projects based on the vertical temperature decline rate, unit precipitation difference value and effective light duration value. This can provide a reference for the construction cost of the current Z effective sub-areas, and can also provide a data basis for the calculation of subsequent construction priority values;

[0199] The construction priority value is the final numerical representation of whether the Z effective sub-areas can be used as the priority areas for vegetation engineering construction, thus truly representing the suitability of the Z effective sub-areas as the priority areas for vegetation engineering construction;

[0200] The expression for constructing the priority value is:

[0201]

[0202] In the formula, JS yxz is the construction priority value of the zth effective sub-area, JS cbz is the construction cost value of the zth effective sub-area, σ6 and σ7 are proportional coefficients greater than 0;

[0203] Among them, the setting logic of σ6 and σ7 is consistent with the setting logic of σ1 and σ2 mentioned above.

[0204] The construction priority area refers to the sub-area in the dry and hot river valley basin that is suitable for vegetation engineering construction, so that the construction priority area can achieve an effective balance between the construction assessment value and the construction cost value, thereby meeting the needs of vegetation engineering construction in the dry and hot river valley basin;

[0205] Methods for constructing construction priority areas include:

[0206] Arrange the z construction priority values ​​in descending order from largest to smallest, and mark the first, second, third, and fourth construction priority values;

[0207] The valid sub-regions corresponding to the first, second, third and fourth construction priority values ​​are recorded as the first region, the second region, the third region and the fourth region respectively;

[0208] The center positions of the first area, the second area, the third area, and the fourth area are identified one by one by computer vision technology to obtain a first center position, a second center position, a third center position, and a fourth center position;

[0209] Taking the first center position, the second center position, the third center position and the fourth center position as expansion base points respectively, the first area, the second area, the third area and the fourth area are expanded outward to obtain the first expansion area, the second expansion area, the third expansion area and the fourth expansion area; the outward expansion method enables the expansion area to expand outward evenly and circumferentially on the basis of the center position, so as to ensure that the expansion area can always be within the corresponding area, thereby ensuring that the location of the subsequent construction priority area does not exceed the range;

[0210] The expansion is stopped until the areas of the first expansion area, the second expansion area, the third expansion area and the fourth expansion area reach 70% of the area of ​​the first area, 40% of the area of ​​the second area, 20% of the area of ​​the third area and 10% of the area of ​​the fourth area respectively; when the construction priority value is larger, it means that the area of ​​the effective sub-area that can be used as the construction priority area is larger, and the area of ​​the construction priority area selected in the effective sub-area is also larger, so that the effective sub-areas participating in the construction of the construction priority area can be reasonably and accurately screened;

[0211] The first expansion area, the second expansion area, the third expansion area and the fourth expansion area after the expansion is stopped are recorded as construction priority areas.

[0212] It should be noted that the construction priority area is not an area in a specific location in terms of geographical location, but an area composed of multiple different locations, so that the construction priority area may exist on the ridges, slopes or river banks of the dry and hot river valley basin, thus achieving the construction effect of multi-regional, multi-location and diversified vegetation engineering construction priority areas in the dry and hot river valley basin.

[0213] In this embodiment, the initial area of ​​the selected area is marked on the electronic map, and the regional critical points are marked on the initial area, the construction area is drawn, the regional boundary line of the construction area is drawn, and based on the regional boundary line, the construction area is divided into sub-areas, the basic geographical parameters of the sub-areas are obtained, the regional stability value of the sub-areas is calculated, and the effective sub-areas are screened out from the sub-areas, the comprehensive micro-meteorological parameters of the effective sub-areas are obtained, and the construction evaluation value of the effective sub-area is calculated based on the comprehensive micro-meteorological parameters, the construction cost value of the effective sub-area is queried, the construction evaluation value and the construction cost value are combined to generate a construction priority value, and a construction priority area is constructed from the effective sub-area; compared with the prior art, by drawing the regional boundary line in the initial area, the boundaries of different geographical environment areas in the dry and hot valley basin can be determined. Marking is carried out, and sub-regions are drawn according to regional boundaries, so as to achieve the effect of breaking up larger areas in the dry and hot river valley basin into smaller areas, which is convenient for the accurate identification of construction priority areas. At the same time, the regional stability value is calculated in combination with the basic geographical parameters, and the effective sub-regions that meet the requirements of vegetation engineering can be screened out from the sub-regions. By collecting the comprehensive micrometeorological parameters of the effective sub-regions, the detailed collection and analysis of the micrometeorological parameters of small areas in the dry and hot river valley basin is realized, and the construction priority values ​​of the effective sub-regions can be calculated in different spatial dimensions, and then the construction priority areas that meet the needs of vegetation engineering construction can be identified, and accurate and appropriate regional locations can be provided for the subsequent construction of vegetation engineering in the dry and hot river valley basin, ensuring that the vegetation in the construction priority areas can grow healthily and rapidly, meeting the ultimate goal of mudslide source control in the dry and hot river valley basin.

[0214] Example 2: Please refer to Figure 2 As shown, the part not described in detail in this embodiment is described in the first embodiment, and a system for constructing a priority area for vegetation engineering construction in a dry and hot river valley basin based on micrometeorology is provided, which is used to implement a method for constructing a priority area for vegetation engineering construction in a dry and hot river valley basin based on micrometeorology, including a region drawing module, a region division module, a region screening module, an evaluation calculation module and a priority area construction module, wherein each module is connected via a wired or wireless network;

[0215] The area drawing module is used to mark the initial area of ​​the selected area on the electronic map, mark the regional critical points on the initial area, and draw the construction area;

[0216] A region division module is used to mark the regional boundary lines of the construction area and divide the construction area into sub-areas based on the regional boundary lines. The regional boundary lines include slope boundary lines, ridge boundary lines and river bank boundary lines.

[0217] The regional screening module is used to obtain the basic geographical parameters of the sub-region, calculate the regional stability value, and screen out the effective sub-regions from the sub-regions. The basic geographical parameters include the altitude difference value and the terrain flatness value;

[0218] An evaluation calculation module is used to obtain the comprehensive micro-meteorological parameters of the effective sub-area and calculate the construction evaluation value based on the comprehensive micro-meteorological parameters. The comprehensive micro-meteorological parameters include the vertical temperature decrease rate, the unit precipitation difference value and the effective light duration value;

[0219] The priority area construction module is used to obtain the construction cost value of the effective sub-area, combine the construction assessment value and the construction cost value, generate the construction priority value, and construct the construction priority area from the effective sub-area.

[0220] The above description is only a specific implementation mode of the present invention, but the protection scope of the present invention is not limited thereto. Any technician familiar with the technical field can easily think of changes or substitutions within the technical scope disclosed by the present invention, which should be covered by the protection scope of the present invention.

Claims

1. A method for constructing priority areas for vegetation engineering construction in dry and hot river valley basins based on micrometeorology, characterized in that: include: S1: Mark the initial area of ​​the selected region on the electronic map, mark the regional critical points on the initial area, and draw the construction area; S2: Mark the regional boundary lines of the construction area and divide the construction area into sub-areas based on the regional boundary lines. The regional boundary lines include slope boundary lines, ridge boundary lines and river bank boundary lines; S3: Obtain basic geographic parameters of the sub-region, calculate the regional stability value, and screen out valid sub-regions from the sub-regions. The basic geographic parameters include altitude difference value and terrain flatness value; S4: Obtain the comprehensive micro-meteorological parameters of the effective sub-area, and calculate the construction evaluation value based on the comprehensive micro-meteorological parameters, the comprehensive micro-meteorological parameters include the vertical temperature decrease rate, the unit precipitation difference value and the effective light duration value; S5: Obtain the construction cost value of the effective sub-area, combine the construction assessment value and the construction cost value to generate a construction priority value, and construct a construction priority area from the effective sub-area.

2. The method for constructing priority areas for vegetation engineering construction in dry and hot river valley basins based on micrometeorology according to claim 1 is characterized in that: The marking method of the initial area comprises: Query the maximum latitude, minimum latitude, maximum longitude and minimum longitude of the selected area through the geographic information system; Mark the points corresponding to the maximum latitude, the minimum latitude, the maximum longitude, and the minimum longitude on the electronic map to obtain a first point, a second point, a third point, and a fourth point; Draw a line horizontally through the positions of the first point and the second point to obtain two latitude lines, and draw a line vertically through the positions of the third point and the fourth point to obtain two longitude lines; Extend the two latitude lines and the two longitude lines respectively until they intersect, and record the area enclosed by the two latitude lines and the two longitude lines as the initial area.

3. The method for constructing priority areas for vegetation engineering construction in dry and hot river valley basins based on micrometeorology according to claim 2 is characterized in that: The method for drawing the construction area includes: Based on the preset interval length, i vertical extension lines and i horizontal extension lines extending into the initial area are drawn on two latitude lines and two longitude lines respectively; Mark p nodes equidistantly on i vertical extension lines and i horizontal extension lines, and query the altitude values ​​of the p nodes one by one through the geographic information system; Along the drawing direction of the vertical extension line and the horizontal extension line, the first node whose elevation value is greater than the elevation values ​​of the two adjacent nodes is recorded as a valid node; After summarizing the first point, the second point, the third point, the fourth point and all valid nodes, k regional critical points are obtained; Taking the regional critical point corresponding to the first point as the starting point of the connection, k regional critical points are connected in a clockwise direction to draw the construction area.

4. The method for constructing priority areas for vegetation engineering construction in dry and hot river valley basins based on micrometeorology according to claim 3 is characterized in that: The sub-area division method includes: The elevation values ​​of all locations in the construction area are queried through the geographic information system, and a digital elevation model of the construction area is constructed in combination with virtual reality technology; Based on the preset side length, the digital elevation model is divided into w cells. The center points of the cells are identified through computer vision technology, and the altitude values ​​of the center points are marked one by one. Record the center point whose altitude value is greater than the preset first altitude threshold as the ridge center point, and record the cell where the ridge center point is located as the ridge cell; The center point whose altitude value is less than the preset second altitude threshold is recorded as the river bank center point, and the cell where the river bank center point is located is recorded as the river bank cell; Recording a center point whose altitude value is between a preset first altitude threshold and a preset second altitude threshold as a slope center point, and recording a cell where the slope center point is located as a slope cell; The outer boundary of the cell is recorded as the target boundary, and lines are drawn along the target boundaries of the ridge cell, the river bank cell and the slope cell respectively to obtain the ridge dividing line, the river bank dividing line and the slope dividing line; Based on the ridge boundary line, river bank boundary line and slope boundary line, the construction area is divided into m closed sub-areas.

5. The method for constructing priority areas for vegetation engineering construction in dry and hot river valley basins based on micrometeorology according to claim 4 is characterized in that: The method for obtaining the altitude difference value includes: Scan the m sub-regions one by one to obtain m regional scan images, and mark all the pixels in the m regional scan images to obtain n pixels; The altitude values ​​of n pixels are queried one by one through the geographic information system, and the maximum altitude value of the pixel point and the minimum altitude value of the pixel point are screened out. The maximum altitude value of m pixels and the minimum altitude value of m pixels are subtracted to obtain m altitude difference values; The expression for the altitude difference is: HB lcm =HB zdm -HB zxm ; In the formula, HB lcm is the elevation difference of the mth sub-area, HB zdm is the maximum altitude value of the pixel points in the mth sub-area, HB zxm is the minimum altitude value of the pixels in the mth sub-region.

6. The method for constructing priority areas for vegetation engineering construction in dry and hot river valley basins based on micrometeorology according to claim 5 is characterized in that: The method for obtaining the terrain flatness value comprises: Randomly select V pixels as the center of the circle in the m area scans, draw a circle with a preset length as the radius, and draw v pixel circles; Count the altitude values ​​of x pixels within the v pixel circle one by one, and add up the altitude values ​​of x pixels and calculate the average to obtain the altitude values ​​of V circles; The expression for the circle altitude value is: In the formula, HB zymv is the circle height value of the vth pixel circle in the mth sub-region, HB XSmva is the altitude value of the a-th pixel point of the V-th pixel circle of the m-th sub-region; After subtracting the vth circle elevation value from the v-1th circle elevation value, we get v-1 sub-flat values; The expression of the sub-flat value is: PT zzmv-1 =HB zvmv -HB zvmv-1 ; Where, PT zzmv-1 is the v-1th sub-flat value of the mth sub-region, HB zymv-1 is the altitude value of the v-1th circle in the mth sub-area; After removing the maximum and minimum values ​​of the sub-flat values, the remaining v-3 sub-flat values ​​are accumulated and averaged to obtain m terrain flatness values; The expression for the flat terrain value is: Where, PT dsm is the terrain flatness value of the mth sub-area, PT zzmb is the bth sub-flat value of the mth sub-region.

7. The method for constructing priority areas for vegetation engineering construction in dry-hot river valley basins based on micrometeorology according to claim 6 is characterized in that: The effective sub-region screening method comprises: Combine the m altitude drop values ​​and the m terrain flatness values ​​to obtain m regional stability values; The expression of the regional stationary value is: In the formula, QY pwm is the regional stability value of the mth sub-region, σ1 and σ2 are proportional coefficients greater than 0; The m regional stationary values ​​QY pwm Sequentially with the preset regional stability threshold QY yz Compare; When QY pwm Less than QY yz When , the mth sub-region is not recognized as a valid sub-region; When QY pwm Greater than or equal to QY yz , the mth sub-region is identified as a valid sub-region, and Z valid sub-regions are obtained.

8. The method for constructing priority areas for vegetation engineering construction in dry-hot river valley basins based on micrometeorology according to claim 7 is characterized in that: The method for obtaining the vertical temperature lapse rate includes: In the Z effective sub-areas, q temperature measurement locations with different altitudes are randomly marked to obtain q temperature measurement points; Taking the same time as the starting time and the preset temperature measurement time as the interval standard, r temperature measurement moments are marked, and the temperature of the r temperature measurement moments is detected by the temperature sensor respectively to obtain r temperature values, and the r temperature values ​​are accumulated and averaged to obtain q sub-temperature measurement values; The expression of the sub-temperature measurement value is: In the formula, CW zzq is the sub-temperature measurement value of the qth temperature measurement point in the zth effective sub-area, WD zzqc is the cth temperature value of the qth temperature measurement point in the zth effective sub-area; The altitude values ​​of q temperature measurement points are queried one by one through the geographic information system to obtain q temperature measurement altitudes, and the qth sub-temperature measurement value is subtracted from the q-1th sub-temperature measurement value, and then compared with the difference between the qth temperature measurement altitude and the q-1th temperature measurement altitude to obtain the q-1 sub-decline rate; The expression of the sub-decline rate is: In the formula, DJ lzq-1 is the q-1th sub-decline rate of the zth effective sub-region, CW zzq-1 is the sub-temperature measurement value of the q-1th temperature measurement point in the zth effective sub-area, HB cwzq is the qth temperature measurement altitude of the zth effective sub-area, HB cwzq-1 is the q-1th temperature measurement altitude of the zth effective sub-area; The q-1 sub-decline rates are accumulated and averaged to obtain z vertical temperature lapse rates; The expression of temperature vertical lapse rate is: In the formula, DJ wdz is the vertical temperature decrease rate of the zth effective sub-region, DJ lzd is the dth sub-decline rate of the zth effective sub-region.

9. The method for constructing priority areas for vegetation engineering construction in dry-hot river valley basins based on micrometeorology according to claim 8 is characterized in that: The method for obtaining the unit precipitation difference value includes: Taking the preset rainfall cycle as the standard, the precipitation data of q temperature measuring points are detected by the rainfall sensor to obtain q precipitation values, and the q precipitation values ​​are respectively compared with the duration corresponding to the preset rainfall cycle to obtain q unit precipitation; The expression of unit precipitation is: In the formula, JS dwzq is the qth unit precipitation in the zth effective sub-area, JS zzq is the qth precipitation value of the zth effective sub-area, SC yl The duration corresponding to the preset rainfall cycle; Subtract the q unit precipitation from the unit precipitation threshold in turn, and then add up the differences and average them to obtain the z unit precipitation difference values; The expression of unit precipitation difference value is: In the formula, JS cez is the unit precipitation difference value of the zth effective sub-area, JS dwze The e-th unit precipitation in the z-th effective sub-area, JS yz is the unit precipitation threshold; The expression to construct the evaluation value is: In the formula, GJ pgz is the construction evaluation value of the zth valid sub-region, GZ scz is the effective illumination duration value of the zth effective sub-area, and σ3, σ4, and σ5 are proportional coefficients greater than 0.

10. The method for constructing priority areas for vegetation engineering construction in dry-hot river valley basins based on micrometeorology according to claim 9 is characterized in that: The expression of the construction priority value is: In the formula, JS yxz is the construction priority value of the zth effective sub-area, JS cbz is the construction cost value of the zth effective sub-area, σ6 and σ7 are proportional coefficients greater than 0; Methods for constructing construction priority areas include: Arrange the z construction priority values ​​in descending order from largest to smallest, and mark the first, second, third, and fourth construction priority values; The valid sub-regions corresponding to the first, second, third and fourth construction priority values ​​are recorded as the first region, the second region, the third region and the fourth region respectively; The center positions of the first area, the second area, the third area, and the fourth area are identified one by one by computer vision technology to obtain a first center position, a second center position, a third center position, and a fourth center position; Taking the first center position, the second center position, the third center position and the fourth center position as expansion base points respectively, expanding outwards on the first area, the second area, the third area and the fourth area to obtain a first expansion area, a second expansion area, a third expansion area and a fourth expansion area; Stop expanding until the areas of the first expansion area, the second expansion area, the third expansion area and the fourth expansion area reach 70% of the area of ​​the first area, 40% of the area of ​​the second area, 20% of the area of ​​the third area and 10% of the area of ​​the fourth area respectively; The first expansion area, the second expansion area, the third expansion area and the fourth expansion area after the expansion is stopped are recorded as construction priority areas.

Images