CN115424136A - Forest canopy health evaluation method and system combining remote sensing and forest map - Google Patents
Forest canopy health evaluation method and system combining remote sensing and forest map Download PDFInfo
- Publication number
- CN115424136A CN115424136A CN202211051390.6A CN202211051390A CN115424136A CN 115424136 A CN115424136 A CN 115424136A CN 202211051390 A CN202211051390 A CN 202211051390A CN 115424136 A CN115424136 A CN 115424136A
- Authority
- CN
- China
- Prior art keywords
- forest
- index
- image
- spectral
- canopy
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
- 230000036541 health Effects 0.000 title claims abstract description 73
- 238000011156 evaluation Methods 0.000 title claims abstract description 54
- 230000003595 spectral effect Effects 0.000 claims abstract description 71
- 238000000034 method Methods 0.000 claims abstract description 35
- 238000002310 reflectometry Methods 0.000 claims abstract description 33
- 238000012545 processing Methods 0.000 claims abstract description 25
- 238000010606 normalization Methods 0.000 claims abstract description 19
- 238000012937 correction Methods 0.000 claims description 12
- 238000004458 analytical method Methods 0.000 claims description 7
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 6
- 239000003086 colorant Substances 0.000 claims description 4
- 238000007637 random forest analysis Methods 0.000 claims description 4
- 230000009466 transformation Effects 0.000 claims description 4
- 101100356682 Caenorhabditis elegans rho-1 gene Proteins 0.000 claims description 3
- 101100468652 Drosophila melanogaster rho-5 gene Proteins 0.000 claims description 3
- 238000010276 construction Methods 0.000 claims description 3
- 238000012847 principal component analysis method Methods 0.000 claims description 3
- 238000007619 statistical method Methods 0.000 claims description 3
- 238000012854 evaluation process Methods 0.000 abstract description 10
- 238000011158 quantitative evaluation Methods 0.000 abstract description 8
- 238000012544 monitoring process Methods 0.000 description 6
- 230000009525 mild injury Effects 0.000 description 4
- 230000009526 moderate injury Effects 0.000 description 4
- 230000009528 severe injury Effects 0.000 description 4
- 238000010586 diagram Methods 0.000 description 2
- 238000010587 phase diagram Methods 0.000 description 2
- 238000000513 principal component analysis Methods 0.000 description 2
- 239000000047 product Substances 0.000 description 2
- 238000005070 sampling Methods 0.000 description 2
- 238000012360 testing method Methods 0.000 description 2
- 208000035126 Facies Diseases 0.000 description 1
- 241000607479 Yersinia pestis Species 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 239000006227 byproduct Substances 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000008641 drought stress Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000012886 linear function Methods 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 230000000750 progressive effect Effects 0.000 description 1
- 238000004445 quantitative analysis Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000012216 screening Methods 0.000 description 1
- 230000000007 visual effect Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06V—IMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
- G06V20/00—Scenes; Scene-specific elements
- G06V20/10—Terrestrial scenes
- G06V20/188—Vegetation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06V—IMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
- G06V10/00—Arrangements for image or video recognition or understanding
- G06V10/20—Image preprocessing
- G06V10/24—Aligning, centring, orientation detection or correction of the image
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06V—IMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
- G06V20/00—Scenes; Scene-specific elements
- G06V20/10—Terrestrial scenes
- G06V20/13—Satellite images
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Multimedia (AREA)
- Theoretical Computer Science (AREA)
- Health & Medical Sciences (AREA)
- General Health & Medical Sciences (AREA)
- Astronomy & Astrophysics (AREA)
- Remote Sensing (AREA)
- Image Processing (AREA)
Abstract
The invention relates to a method and a system for evaluating canopy health by combining remote sensing and a forest map, wherein the method comprises the following steps: determining a multispectral reflectivity image of the area to be detected according to the high-resolution six-satellite wide image; obtaining a forest area image according to the forest phase image and the multispectral reflectivity image; determining a spectral index according to the forest area image; carrying out data normalization processing and weight determination on the spectral indexes to obtain weight values corresponding to the spectral indexes; constructing a canopy health evaluation index according to the weight value and the spectral index; and evaluating the canopy of the forest region according to the calculated canopy health evaluation index. The method integrates satellite images and spectral index information, carries out quantitative evaluation on the health condition of the canopy, improves the objectivity of the evaluation result, is easy to operate in the evaluation process, and improves the reliability of the evaluation process.
Description
Technical Field
The invention relates to the technical field of quantitative evaluation of canopy health, in particular to a canopy health evaluation method and system combining remote sensing and a forest facies map.
Background
The method for evaluating the health of the forest canopy at home and abroad is mainly based on the evaluation of ground sampling survey or by means of man-machine aerial survey and visual interpretation, the methods provide valuable information for forest managers, but the evaluation of the health of the forest canopy through manual survey not only needs to invest a large amount of manpower, material resources and financial resources and is difficult to ensure the time consistency in a large area, but also has complicated evaluation indexes and subjectivity; the manned and airborne survey has high taxi fee, is often limited by weather and airspace control and the like, and is difficult to acquire all-area data in time; the satellite remote sensing can effectively monitor the health state of the forest vegetation canopy by virtue of the advantages of rapidness and macroscopicity.
Compared with ground sampling survey and manned flight survey methods, the method for monitoring and evaluating the health condition of the canopy by using the satellite remote sensing technology has the advantages of wide monitoring range, capability of monitoring at intervals of time of year, month and day, and more detailed and more spatial results. However, the existing research mostly focuses on evaluating the overall health level of the regional ecological system or the forest, has multiple and complex indexes, the analysis on the forest growth condition is usually only aimed at single conditions such as drought stress or certain pest control, the selection of index factors is also different from person to person, and a simple, convenient and easy-to-operate quantitative method for monitoring and evaluating the health condition of the canopy is still lacking at present.
Disclosure of Invention
In order to overcome the defects of the prior art, the invention aims to provide a canopy health evaluation method and a canopy health evaluation system combining remote sensing and a forest phase diagram.
In order to achieve the purpose, the invention provides the following scheme:
a canopy health evaluation method combining remote sensing and a forest map comprises the following steps:
determining a multispectral reflectivity image of the area to be detected according to the high-resolution six-satellite wide image;
obtaining a forest area image according to the forest phase image and the multispectral reflectivity image;
determining a spectral index according to the forest area image;
carrying out data normalization processing and weight determination on the spectral indexes to obtain weight values corresponding to the spectral indexes;
constructing a canopy health evaluation index according to the weight value and the spectral index;
and evaluating the canopy of the forest region according to the calculated canopy health evaluation index.
Preferably, the determining the multispectral reflectivity image of the region to be detected according to the high-resolution six-number satellite wide image includes:
acquiring an L1A-level multispectral image of the wide image of the high-resolution six-number satellite in clear sky in the forest vegetation growing season in the area to be detected;
carrying out absolute radiometric calibration and atmospheric correction on the L1A-level multispectral image to obtain a first corrected image;
performing geometric correction on the first correction image based on the RPC of the product data and the reference image to obtain a second correction image;
and splicing, cutting and cloud mask processing are carried out on the second corrected image to obtain the multispectral reflectivity image.
Preferably, the obtaining a forest area image according to a forest map and the multispectral reflectivity image comprises:
merging the image spots of the forest map to generate a forest area and a non-forest area;
and performing mask processing on the multispectral reflectivity image, reserving original reflectivity values of all wave band values of the image covering the forest area, and setting the reflectivity values of all wave band values of the non-forest area to be 0 to obtain the forest area image.
Preferably, the determining a spectral index from the forest area image comprises:
carrying out statistical analysis on spectral indexes of the selected healthy and sub-healthy samples in the forest region image by adopting a random forest regression method to obtain an analysis result;
and selecting the spectral index which reflects that the correlation of the greenness, the moisture and the disturbed state of the tree crown leaves is greater than a preset threshold value according to the analysis result.
Preferably, the spectral indices include greenness normalized vegetation index, greenness ratio index, red-edge normalized vegetation index, yellow-edge normalized vegetation index, normalized water body index, and fire passing index.
Preferably, the formula for calculating the spectral index is:
GNDVI is the greenness normalized vegetation index, GVI is the greenness ratio index, RENDVI is the red-edge normalized vegetation index, YENDVI is the yellow-edge normalized vegetation index, NDWI is the normalized water body index, BAI is the fire passing index, and rho 1, rho 2, rho 3, rho 4, rho 5, rho 6 and rho 8 are the reflectivities of wave band 1, wave band 2, wave band 3, wave band 4, wave band 5, wave band 6 and wave band 8 of the high-resolution six-satellite wide-width image in sequence.
Preferably, the performing data normalization processing and weight determination on the spectral indexes to obtain weight values corresponding to the spectral indexes includes:
using a formulaNormalizing each spectral index to obtain normalized data; wherein, X ni The normalized value of a certain spectral index in a pixel i; x i The value of the spectral index in a pixel i is shown; x max And X min Respectively the maximum value and the minimum value of the statistics of the spectral index;
and performing linear transformation on the normalized data based on a principal component analysis method to obtain the weight value.
Preferably, the formula of the canopy health evaluation index is as follows:
FCHI=a×G+b×R+c×Y+d×V+e×N+f×B;
FCHI is the health evaluation index of the canopy, and a, B, c, d, e and f are weight values of G, R, Y, V, N and B in sequence; g, R, Y, V, N and B are values after GNDVI, RENDVI, GVI, NDWI and BAI normalization processing in sequence.
Preferably, evaluating the canopy of the forest region according to the calculated canopy health evaluation index includes:
using a formulaStandardizing the canopy health evaluation index to obtain a standardized index; wherein, FCHI · For the normalized index, FCHI max 、FCHI min Maximum and minimum health index values within the region, respectively;
and classifying the result obtained by the standardized index into healthy, mild damage, moderate damage and severe damage types according to a natural breakpoint method, and giving different colors to the result and then outputting the drawing to obtain an evaluation result.
A canopy health evaluation system combining remote sensing and a forest map comprises:
the first image determining unit is used for determining a multispectral reflectivity image of the area to be detected according to the high-resolution six-satellite wide image;
the second image determining unit is used for obtaining a forest area image according to the forest phase image and the multispectral reflectivity image;
the index determining unit is used for determining a spectral index according to the forest area image;
the weight determining unit is used for carrying out data normalization processing and weight determination on the spectral indexes to obtain weight values corresponding to the spectral indexes;
the index construction unit is used for constructing a canopy health evaluation index according to the weight value and the spectral index;
and the evaluation unit is used for evaluating the forest crowns in the forest area according to the calculated forest crown health evaluation index.
According to the specific embodiment provided by the invention, the invention discloses the following technical effects:
the invention provides a canopy health evaluation method and system combining remote sensing and a forest map, wherein the method comprises the following steps: determining a multispectral reflectivity image of the area to be detected according to the high-resolution six-satellite wide image; obtaining a forest area image according to the forest phase image and the multispectral reflectivity image; determining a spectral index according to the forest area image; carrying out data normalization processing and weight determination on the spectral indexes to obtain weight values corresponding to the spectral indexes; constructing a canopy health evaluation index according to the weight value and the spectral index; and evaluating the canopy of the forest region according to the calculated canopy health evaluation index. The method integrates satellite images and spectral index information, carries out quantitative evaluation on the health condition of the canopy, improves the objectivity of the evaluation result, is easy to operate in the evaluation process, and improves the reliability of the evaluation process.
Drawings
In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the drawings required in the embodiments will be briefly described below, and it is obvious that the drawings in the following description are only some embodiments of the present invention, and it is obvious for those skilled in the art that other drawings can be obtained according to these drawings without creative efforts.
FIG. 1 is a flow chart of a method in an embodiment provided by the present invention;
fig. 2 is a schematic diagram of an evaluation process in an embodiment provided by the present invention.
Detailed Description
The technical solutions in the embodiments of the present invention will be clearly and completely described below with reference to the drawings in the embodiments of the present invention, and it is obvious that the described embodiments are only a part of the embodiments of the present invention, and not all of the embodiments. All other embodiments, which can be obtained by a person skilled in the art without making any creative effort based on the embodiments in the present invention, belong to the protection scope of the present invention.
The invention aims to provide a method and a system for evaluating the health of a canopy by combining remote sensing and a forest phase diagram, which integrate satellite images and spectral index information, carry out quantitative evaluation on the health condition of the canopy, improve the objectivity of an evaluation result, facilitate the operation of an evaluation process and improve the reliability of the evaluation process.
The high-resolution six-satellite 16m wide data (GF-6 WFV) has the width of 800km, and a 4d repetition period obtains 8 spectral band image data of the same region, which can be used for distinguishing vegetation physiological characteristics, so that the method can be used for forest canopy health monitoring and evaluation in a large region. The invention provides a forest canopy health quantitative evaluation method integrating spectral index information of satellite images by combining GF-6WFV images and forest map data, and aims to provide an objective, reliable and easy-to-operate quantitative evaluation method for forest canopy health monitoring and managers.
In order to make the aforementioned objects, features and advantages of the present invention comprehensible, embodiments accompanied with figures are described in further detail below.
Fig. 1 is a flowchart of a method in an embodiment provided by the present invention, and as shown in fig. 1, the present invention provides a method for evaluating canopy health by combining remote sensing and a forest map, including:
step 100: determining a multispectral reflectivity image of the area to be detected according to the high-resolution six-satellite wide image;
step 200: obtaining a forest area image according to the forest map and the multispectral reflectivity image;
step 300: determining a spectral index according to the forest area image;
step 400: carrying out data normalization processing and weight determination on the spectral indexes to obtain weight values corresponding to the spectral indexes;
step 500: constructing a canopy health evaluation index according to the weight value and the spectral index;
step 600: and evaluating the canopy of the forest region according to the calculated canopy health evaluation index.
Preferably, the step 100 specifically includes:
acquiring an L1A-level multispectral image of the wide image of the high-resolution six-number satellite in clear sky in the forest vegetation growing season in the area to be detected;
carrying out absolute radiometric calibration and atmospheric correction on the L1A-level multispectral image to obtain a first corrected image;
performing geometric correction on the first corrected image based on the RPC of the product data and the reference image to obtain a second corrected image;
and splicing, cutting and cloud mask processing are carried out on the second corrected image to obtain the multispectral reflectivity image.
Fig. 2 is a schematic diagram of an evaluation process in an embodiment provided by the present invention, and as shown in fig. 2, the first step in this embodiment is image selection and processing, namely, selecting a clear sky high-resolution six-size satellite wide-width image (GF-6 WFV) L1A-level multispectral image in a forest vegetation growing season of a test area, then performing absolute radiometric calibration and atmospheric correction on the selected image, performing geometric correction on the multispectral image based on RPC and a reference image carried by product data, and finally performing splicing or clipping, cloud mask processing, and the like to obtain a cloud-free multispectral reflectance image with a spatial resolution of the test area of 16 m.
Preferably, the step 200 specifically includes:
merging the pattern spots of the forest map to generate a forest area and a non-forest area;
and performing mask processing on the multispectral reflectivity image, reserving original reflectivity values of all wave band values of the image covering the forest area, and setting the reflectivity values of all wave band values of the non-forest area to be 0 to obtain the forest area image.
Specifically, in this embodiment, the second step is to select a forest region image, merge forest map patches to generate a forest region and a non-forest region, perform mask processing on the processed GF-6WFV reflectance image, retain original reflectance values of each waveband of the image covering the forest region, and set reflectance values of each waveband of the rest of the non-forest region to 0.
Preferably, the step 300 specifically includes:
carrying out statistical analysis on spectral indexes of the selected healthy and sub-healthy samples in the forest region image by adopting a random forest regression method to obtain an analysis result;
and selecting the spectral index which reflects that the correlation of the greenness, the moisture and the disturbed state of the tree crown leaves is greater than a preset threshold value according to the analysis result.
Further, in the third step of this embodiment, spectral index screening is performed, a random forest regression method is adopted to statistically analyze spectral indexes of selected healthy and sub-healthy (mild damage, moderate damage, and severe damage) samples for GF-6WFV images of covered forest regions generated by processing, and then spectral indexes with high correlation of states such as greenness, moisture, disturbance and the like of forest canopy leaves are selected and used for constructing forest canopy health quantitative evaluation indexes.
Preferably, the spectral indices include greenness normalized vegetation index, greenness ratio index, red-edge normalized vegetation index, yellow-edge normalized vegetation index, normalized water body index, and fire passing index.
Preferably, the formula for calculating the spectral index is:
GNDVI is the greenness normalized vegetation index, GVI is the greenness ratio index, RENDVI is the red-edge normalized vegetation index, YENDVI is the yellow-edge normalized vegetation index, NDWI is the normalized water body index, BAI is the fire passing index, and rho 1, rho 2, rho 3, rho 4, rho 5, rho 6 and rho 8 are the reflectivities of wave band 1, wave band 2, wave band 3, wave band 4, wave band 5, wave band 6 and wave band 8 of the high-resolution six-satellite wide-width image in sequence.
Preferably, the step 400 specifically includes:
using a formulaNormalizing each spectral index to obtain normalized data; wherein, X ni The normalized value of a certain spectral index in the pixel i; x i The value of the spectral index in a pixel i is shown; x max And X min Respectively the maximum value and the minimum value of the statistics of the spectral index;
and performing linear transformation on the normalized data based on a principal component analysis method to obtain the weight value.
Optionally, the fourth step in this embodiment is data normalization processing, and since the selected and calculated exponential factor dimensions are not completely uniform, direct use for modeling may cause imbalance of each index. Therefore, the linear function normalization method is adopted to unify the dimensions of the exponential factors to [0,1], and the normalization processing is carried out on the exponential factors. The normalized formula is:
in the formula, X ni The value of a certain factor in the pixel i after normalization; x i The value of the factor in the pixel i; x max And X min Respectively the maximum and minimum of the factor statistic.
Further, in the fifth step of this embodiment, the weight of each factor is determined, and Principal Component Analysis (PCA) is used to convert the multidimensional remote sensing index information into a few variables that are easy to interpret and are not related to each other for 6 index factors after data normalization, where the first principal component (PC 1) of the variables contains most of the information in the original data set. The 6 variables can be transformed by linear transformation by utilizing the PC1 to construct the canopy health index, so that the result deviation caused by the conventional artificial subjective weight setting is avoided.
Preferably, the formula of the canopy health evaluation index is as follows:
FCHI=a×G+b×R+c×Y+d×V+e×N+f×B;
FCHI is the health evaluation index of the canopy, and a, B, c, d, e and f are weight values of G, R, Y, V, N and B in sequence; g, R, Y, V, N and B are values after GNDVI, RENDVI, GVI, NDWI and BAI normalization processing in sequence.
Specifically, in this embodiment, the sixth step is to construct a canopy health evaluation index, and calculate a canopy health evaluation index (FCHI) according to equation (8) by using weight values of 6 index factors in PC 1.
FCHI=a×G+b×R+c×Y+d×V+e×N+f×B
In the formula, a, B, c, d, e and f are weight values of G, R, Y, V, N and B in sequence; g, R, Y, V, N and B are values after GNDVI, RENDVI, GVI, NDWI and BAI normalization processing in sequence.
Preferably, the step 600 specifically includes:
using the formulaStandardizing the canopy health evaluation index to obtain a standardized index; wherein, FCHI · For the normalized index, FCHI max 、FCHI min Maximum and minimum health index values within the region, respectively;
and classifying the result obtained by the standardized index into healthy, mild damage, moderate damage and severe damage types according to a natural breakpoint method, and giving different colors to the result and then outputting the drawing to obtain an evaluation result.
Specifically, in this embodiment, the last step is the classification of canopy health, in order to facilitate the measurement and comparison of evaluation indexes, the FCHI result is normalized according to the following formula, the values are unified to [0, 100], and the types of health, mild damage, moderate damage and severe damage are classified according to the natural breakpoint method, and different colors are given for chart output.
In the formula: FCHI max 、FCHI min Are respectively asMaximum and minimum health index values within a region. The closer the FCHI value is to 100, the healthier the forest canopy leaves are; otherwise, the less healthy it is.
Corresponding to the method, the embodiment provides a canopy health evaluation system combining remote sensing and a forest map, which includes:
the first image determining unit is used for determining a multispectral reflectivity image of the area to be detected according to the high-resolution six-satellite wide image;
the second image determining unit is used for obtaining a forest area image according to the forest phase image and the multispectral reflectivity image;
the index determining unit is used for determining a spectral index according to the forest area image;
the weight determining unit is used for carrying out data normalization processing and weight determination on the spectral indexes to obtain weight values corresponding to the spectral indexes;
the index construction unit is used for constructing a canopy health evaluation index according to the weight value and the spectral index;
and the evaluation unit is used for evaluating the forest crowns in the forest area according to the calculated forest crown health evaluation index.
The invention has the following beneficial effects:
the method integrates satellite images and spectral index information, carries out quantitative evaluation on the health condition of the canopy, improves the objectivity of the evaluation result, is easy to operate in the evaluation process, and improves the reliability of the evaluation process.
The embodiments in the present description are described in a progressive manner, each embodiment focuses on differences from other embodiments, and the same and similar parts among the embodiments are referred to each other.
The principles and embodiments of the present invention have been described herein using specific examples, which are provided only to help understand the method and the core concept of the present invention; meanwhile, for a person skilled in the art, according to the idea of the present invention, the specific embodiments and the application range may be changed. In view of the foregoing, the description is not to be taken in a limiting sense.
Claims (10)
1. A canopy health evaluation method combining remote sensing and a forest map is characterized by comprising the following steps:
determining a multispectral reflectivity image of the area to be detected according to the high-resolution six-number satellite wide image;
obtaining a forest area image according to the forest map and the multispectral reflectivity image;
determining a spectral index according to the forest area image;
carrying out data normalization processing and weight determination on the spectral indexes to obtain weight values corresponding to the spectral indexes;
constructing a canopy health evaluation index according to the weight value and the spectral index;
and evaluating the canopy of the forest region according to the calculated canopy health evaluation index.
2. The method for canopy health assessment by remote sensing in combination with a forest map of claim 1, wherein determining a multispectral reflectance image of the area to be assessed from the high-resolution six-satellite wide-width image comprises:
acquiring an L1A-level multispectral image of the wide image of the high-resolution six-number satellite in clear sky in the forest vegetation growing season in the area to be detected;
carrying out absolute radiometric calibration and atmospheric correction on the L1A-level multispectral image to obtain a first corrected image;
performing geometric correction on the first correction image based on the RPC of the product data and the reference image to obtain a second correction image;
and splicing, cutting and cloud mask processing are carried out on the second corrected image to obtain the multispectral reflectivity image.
3. The method for canopy health assessment using remote sensing and forest map together of claim 1, wherein said obtaining a forest area image from a forest map and said multispectral reflectance image comprises:
merging the pattern spots of the forest map to generate a forest area and a non-forest area;
and performing mask processing on the multispectral reflectivity image, reserving original reflectivity values of all wave band values of the image covering the forest area, and setting the reflectivity values of all wave band values of the non-forest area to be 0 to obtain the forest area image.
4. The method for evaluating canopy health by combining remote sensing with a forest map of claim 1, wherein determining a spectral index from the forest area image comprises:
carrying out statistical analysis on spectral indexes of the health and sub-health samples in the selected forest region image by adopting a random forest regression method to obtain an analysis result;
and selecting the spectral index which reflects that the greenness, the moisture and the disturbed state of the tree leaves are more than a preset threshold value according to the analysis result.
5. The method of claim 4, wherein the spectral indices comprise greenness normalized vegetation index, greenness ratio index, red-edge normalized vegetation index, yellow-edge normalized vegetation index, normalized water body index, and fire index.
6. The method for forest crown health assessment by combining remote sensing and forest map according to claim 5, wherein the formula for calculating the spectral index is as follows:
GNDVI is the greenness normalized vegetation index, GVI is the greenness ratio index, RENDVI is the red-edge normalized vegetation index, YENDVI is the yellow-edge normalized vegetation index, NDWI is the normalized water body index, BAI is the fire passing index, and rho 1, rho 2, rho 3, rho 4, rho 5, rho 6 and rho 8 are the reflectivities of wave band 1, wave band 2, wave band 3, wave band 4, wave band 5, wave band 6 and wave band 8 of the high-resolution six-satellite wide-width image in sequence.
7. The method for canopy health assessment with combined remote sensing and forest map of claim 6, wherein the performing data normalization and weight determination on the spectral indices to obtain weight values corresponding to each of the spectral indices comprises:
using a formulaNormalizing each spectral index to obtain normalized data; wherein, X ni The normalized value of a certain spectral index in a pixel i; x i The value of the spectral index in a pixel i is obtained; x max And X min Are respectively provided withAre the statistical maximum and minimum values of the spectral index;
and performing linear transformation on the normalized data based on a principal component analysis method to obtain the weight value.
8. The method for forest crown health assessment by combining remote sensing and forest map according to claim 7, wherein the formula of the forest crown health assessment index is as follows:
FCHI=a×G+b×R+c×Y+d×V+e×N+f×B;
FCHI is the health evaluation index of the canopy, and a, B, c, d, e and f are weight values of G, R, Y, V, N and B in sequence; g, R, Y, V, N and B are values after GNDVI, RENDVI, GVI, NDWI and BAI normalization processing in sequence.
9. The method for evaluating the health of a forest crown by combining remote sensing and a forest map according to claim 8, wherein evaluating the forest crown of a forest area according to the calculated forest crown health evaluation index comprises:
using a formulaStandardizing the canopy health evaluation index to obtain a standardized index; wherein, FCHI · Is the normalized index, FCHI max 、FCHI min Maximum and minimum health index values within the region, respectively;
and classifying the result obtained by the standardized index into healthy, mild victim, moderate victim and severe victim types according to a natural breakpoint method, and giving different colors for drawing output to obtain an evaluation result.
10. A canopy health evaluation system combining remote sensing and a forest map is characterized by comprising:
the first image determining unit is used for determining a multispectral reflectivity image of the area to be detected according to the high-resolution six-satellite wide image;
the second image determining unit is used for obtaining a forest area image according to the forest phase image and the multispectral reflectivity image;
the index determining unit is used for determining a spectral index according to the forest area image;
the weight determining unit is used for carrying out data normalization processing and weight determination on the spectral indexes to obtain weight values corresponding to the spectral indexes;
the index construction unit is used for constructing a canopy health evaluation index according to the weight value and the spectral index;
and the evaluation unit is used for evaluating the canopy of the forest area according to the calculated canopy health evaluation index.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202211051390.6A CN115424136A (en) | 2022-08-31 | 2022-08-31 | Forest canopy health evaluation method and system combining remote sensing and forest map |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202211051390.6A CN115424136A (en) | 2022-08-31 | 2022-08-31 | Forest canopy health evaluation method and system combining remote sensing and forest map |
Publications (1)
Publication Number | Publication Date |
---|---|
CN115424136A true CN115424136A (en) | 2022-12-02 |
Family
ID=84200018
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202211051390.6A Pending CN115424136A (en) | 2022-08-31 | 2022-08-31 | Forest canopy health evaluation method and system combining remote sensing and forest map |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN115424136A (en) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN116935238A (en) * | 2023-07-07 | 2023-10-24 | 滁州学院 | Forest disturbance monitoring method, system, equipment and medium based on deep learning |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN203350719U (en) * | 2013-07-03 | 2013-12-18 | 广州地理研究所 | Multispectral remote sensing system of single rotor micro unmanned plane |
CN106324617A (en) * | 2016-10-31 | 2017-01-11 | 塔里木大学 | Satellite remote sensing monitoring method for content of jujube tree canopy carotenoid |
CN106504240A (en) * | 2016-10-31 | 2017-03-15 | 塔里木大学 | A kind of satellite remote-sensing monitoring method of jujube tree canopy Chlorophyll-a Content |
CN113689490A (en) * | 2021-09-01 | 2021-11-23 | 崔宁博 | Fruit tree leaf area index inversion method based on random forest algorithm |
CN113945527A (en) * | 2021-11-15 | 2022-01-18 | 江苏天汇空间信息研究院有限公司 | Method for obtaining water quality total phosphorus parameter inversion optimal model based on satellite data |
CN114694036A (en) * | 2022-03-18 | 2022-07-01 | 南京农业大学 | High-altitude area crop classification and identification method based on high-resolution images and machine learning |
-
2022
- 2022-08-31 CN CN202211051390.6A patent/CN115424136A/en active Pending
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN203350719U (en) * | 2013-07-03 | 2013-12-18 | 广州地理研究所 | Multispectral remote sensing system of single rotor micro unmanned plane |
CN106324617A (en) * | 2016-10-31 | 2017-01-11 | 塔里木大学 | Satellite remote sensing monitoring method for content of jujube tree canopy carotenoid |
CN106504240A (en) * | 2016-10-31 | 2017-03-15 | 塔里木大学 | A kind of satellite remote-sensing monitoring method of jujube tree canopy Chlorophyll-a Content |
CN113689490A (en) * | 2021-09-01 | 2021-11-23 | 崔宁博 | Fruit tree leaf area index inversion method based on random forest algorithm |
CN113945527A (en) * | 2021-11-15 | 2022-01-18 | 江苏天汇空间信息研究院有限公司 | Method for obtaining water quality total phosphorus parameter inversion optimal model based on satellite data |
CN114694036A (en) * | 2022-03-18 | 2022-07-01 | 南京农业大学 | High-altitude area crop classification and identification method based on high-resolution images and machine learning |
Non-Patent Citations (5)
Title |
---|
国务院第一次全国自然灾害综合风险普查领导小组办公室: "《气象灾害调查与风险评估》", 北京:应急管理出版社, pages: 76 * |
张悦琦 等: "基于多时相 GF-6 遥感影像的水稻种植面积提取", 《农业工程学报》, vol. 37, no. 17, pages 189 - 196 * |
曾庆伟 等: "GF-6WFV黄边波段在森林扰动监测中的应用潜力研究", vol. 34, no. 2, pages 211 * |
洪娇 等: "基于Hyperion数据的滇西北高寒山区高山松典型森林生态系统健康分级研究", vol. 36, no. 6, pages 80 - 84 * |
王仁军 等: "基于 GF⁃6 的可可西里典型湖泊水体面积提取方法研究", 《测绘通报》, vol. 2022, no. 5, pages 32 - 37 * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN116935238A (en) * | 2023-07-07 | 2023-10-24 | 滁州学院 | Forest disturbance monitoring method, system, equipment and medium based on deep learning |
CN116935238B (en) * | 2023-07-07 | 2024-02-23 | 滁州学院 | Forest disturbance monitoring method, system, equipment and medium based on deep learning |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Lu et al. | Decomposition of vegetation cover into woody and herbaceous components using AVHRR NDVI time series | |
CN107067334B (en) | Farmland fire and straw burning monitoring method and system based on multi-source remote sensing data | |
Ehrlich et al. | Applications of NOAA-AVHRR 1 km data for environmental monitoring | |
CN104089647B (en) | A kind of crop pest occurrence scope monitoring method and system | |
Gholizadeh et al. | Multi‐temporal assessment of grassland α‐and β‐diversity using hyperspectral imaging | |
CN110927120B (en) | Early warning method for vegetation coverage | |
King et al. | Modelling and mapping damage to forests from an ice storm using remote sensing and environmental data | |
Kuhnell et al. | Mapping woody vegetation cover over the state of Queensland using Landsat TM imagery | |
CN114821349A (en) | Forest biomass estimation method considering harmonic model coefficients and phenological parameters | |
CN115424136A (en) | Forest canopy health evaluation method and system combining remote sensing and forest map | |
Sabzchi-Dehkharghani et al. | Recognition of different yield potentials among rain-fed wheat fields before harvest using remote sensing | |
Trier et al. | Multi-sensor forest vegetation height mapping methods for Tanzania | |
CN113343808A (en) | Tropical forest resource measuring method based on satellite remote sensing technology | |
Pérez-Planells et al. | Retrieval of land surface emissivities over partially vegetated surfaces from satellite data using radiative transfer models | |
Sevillano Marco et al. | Improvement of existing and development of future Copernicus land monitoring products–the ECOLASS project | |
Kukkonen et al. | Volumes by tree species can be predicted using photogrammetric UAS data, Sentinel-2 images and prior field measurements | |
Hovi et al. | Synergistic use of multi-and hyperspectral remote sensing data and airborne LiDAR to retrieve forest floor reflectance | |
Dutta | Assessment of tea bush health and yield using geospatial techniques | |
CN117907248B (en) | Remote sensing monitoring method and system for root system soil water content in key growth period of winter wheat | |
Cihlar et al. | Microwave remote sensing of agricultural crops in Canada | |
Marshall et al. | Drought detection in semi-arid regions using remote sensing of vegetation indices and drought indices | |
Han et al. | Extraction Method of Stand Density Based on High-Resolution Remote Sensing Imagery | |
Zhuang et al. | Relaxation vegetation index in non-linear modelling of ground plant cover by satellite remote-sensing data | |
do Amaral et al. | Geoprocessing Applied to Crop Management | |
Hazaymeh | Development of a remote sensing-based agriculture monitoring drought index and its application over semi-arid region |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination |