CN113989661B - Normalized vegetation index conversion method of MERSI-2 and MODIS - Google Patents

Normalized vegetation index conversion method of MERSI-2 and MODIS Download PDF

Info

Publication number
CN113989661B
CN113989661B CN202111205935.XA CN202111205935A CN113989661B CN 113989661 B CN113989661 B CN 113989661B CN 202111205935 A CN202111205935 A CN 202111205935A CN 113989661 B CN113989661 B CN 113989661B
Authority
CN
China
Prior art keywords
modis
resolution imaging
imaging spectrometer
vegetation index
mersi
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.)
Active
Application number
CN202111205935.XA
Other languages
Chinese (zh)
Other versions
CN113989661A (en
Inventor
周广胜
何奇瑾
任鸿瑞
刘二华
汲玉河
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Taiyuan University of Technology
China Agricultural University
Chinese Academy of Meteorological Sciences CAMS
Original Assignee
Taiyuan University of Technology
China Agricultural University
Chinese Academy of Meteorological Sciences CAMS
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Taiyuan University of Technology, China Agricultural University, Chinese Academy of Meteorological Sciences CAMS filed Critical Taiyuan University of Technology
Priority to CN202111205935.XA priority Critical patent/CN113989661B/en
Publication of CN113989661A publication Critical patent/CN113989661A/en
Application granted granted Critical
Publication of CN113989661B publication Critical patent/CN113989661B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2218/00—Aspects of pattern recognition specially adapted for signal processing
    • G06F2218/02—Preprocessing
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10—Complex mathematical operations
    • G06F17/15—Correlation function computation including computation of convolution operations
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10—Complex mathematical operations
    • G06F17/18—Complex mathematical operations for evaluating statistical data, e.g. average values, frequency distributions, probability functions, regression analysis
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2218/00—Aspects of pattern recognition specially adapted for signal processing
    • G06F2218/08—Feature extraction
    • Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
    • Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
    • Y02A90/10—Information and communication technologies [ICT] supporting adaptation to climate change, e.g. for weather forecasting or climate simulation

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Data Mining & Analysis (AREA)
  • Pure & Applied Mathematics (AREA)
  • Computational Mathematics (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Mathematical Physics (AREA)
  • General Engineering & Computer Science (AREA)
  • Software Systems (AREA)
  • Databases & Information Systems (AREA)
  • Algebra (AREA)
  • Artificial Intelligence (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Computing Systems (AREA)
  • Signal Processing (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Bioinformatics & Cheminformatics (AREA)
  • Bioinformatics & Computational Biology (AREA)
  • Evolutionary Biology (AREA)
  • Operations Research (AREA)
  • Probability & Statistics with Applications (AREA)
  • Investigating Or Analysing Materials By Optical Means (AREA)

Abstract

The invention provides a normalized vegetation index conversion method of FY-3E MERSI-2 and MODIS, relating to the technical field of remote sensing science, the reflection spectrum information of MODIS and FY-3E MERSI-2 in corresponding wavelength ranges is extracted from a green plant reflection spectrum characteristic curve graph to obtain a normalized vegetation index expression of a resolution imaging spectrometer in the MODIS and a normalized vegetation index expression of the resolution imaging spectrometer in the FY-3E MERSI-2, and then the normalized vegetation index of the MODIS is converted into a conversion expression of the normalized vegetation index of the FY-3E MERSI-2 by combining a linear function expression of the resolution imaging spectrometer in the MODIS on wavelength and reflectivity, the normalized vegetation index conversion of FY-3E MERSI-2 and MODIS is realized, and the consistency of the two on earth vegetation and ecological environment monitoring is effectively ensured.

Description

Normalized vegetation index conversion method of MERSI-2 and MODIS
Technical Field
The invention relates to the technical field of remote sensing science, in particular to a method, a system, electronic equipment, a non-transient computer readable storage medium and a computer program product for converting normalized vegetation indexes of FY-3E MERSI-2 and MODIS.
Background
The normalized vegetation index (NDVI) is the ratio of the difference between the reflectivity of the near-infrared band and the reflectivity of the red-light band to the sum of the reflectivity of the near-infrared band and the reflectivity of the red-light band, and is one of the most commonly used remote sensing vegetation indexes.
The first advanced remote sensing satellite Terra of the polar earth orbit environment, which successfully launched the Earth Observation System (EOS) in the united states, 18 months 1999; in 2002, 5, 4 months and successfully emits the remote sensing satellite aquaa. A moderate resolution imaging spectrometer (MODIS) carried on two satellites of Terra and Aqua is an important instrument for observing global biological and physical processes in the United states earth observation system plan, the orbit height of the MODIS is 705km, the scanning swath is 2330km, the bandwidth is 10km, 36 observation bands are provided, three spatial resolutions of 250m, 500m and 1000m are provided, and full spectral coverage from 0.4 mu m (visible light) to 14.4 mu m (thermal infrared) can be realized. The MODIS are simultaneously mounted on two satellites, Terra (morning star) and Aqua (afternoon star), and observe the earth's surface once every one to two days. Both Terra and Aqua are polar orbiting satellites, Terra crosses the equator north to south in the morning when local (10:30), Aqua crosses the equator north to south in the afternoon (1:30), Terra and Aqua are observed for most of the central countries up to 4 times per day. In the NDVI calculation based on MODIS, the spatial resolution of the red light band and the near infrared band is 250m, and the corresponding wavelength ranges are 620-670 nm and 841-876 nm respectively.
At present, domestic wind and cloud satellites in China are developed rapidly, and the wind and cloud satellites are applied widely. The wind cloud No. three E star (FY-3E) is the fifth satellite of a wind cloud No. three polar orbit meteorological satellite series, is emitted in 7 months and 5 days in 2021, is the first morning and evening orbit earth observation satellite in the wind cloud satellite series and even world business meteorological satellites, a plurality of sets of advanced remote sensing instruments are carried on the satellite, and the medium-resolution imaging spectrometer MERSI-2 is one of core instruments thereof, and plays an important role in improving the precision and the timeliness of vegetation and ecological environment monitoring. The morning and evening orbit is a satellite observation orbit about 6:00 of local descending intersection point time (ETC), the morning and evening orbit observation can effectively make up for the defect that the assimilation time window of every 6 hours has a blank satellite observation data, and the morning and evening orbit satellite three-star networking observation can enable 100% of satellite data in the assimilation window of 6 hours to cover the whole world, and has an important effect on improving the accuracy and the timeliness of global numerical weather forecast. MERSI-2 is one of core instruments of the FY-3E satellite, integrates the functions of a resolution ratio spectral imager (MERSI-1) and a visible light infrared scanning radiometer (VIRR) in the original wind cloud three-number satellite, has 25 channels, can seamlessly obtain a 250m resolution ratio true color image of a global daily scale, can carry out continuous comprehensive observation on the atmosphere, the land and the ocean, and realizes high-precision quantitative remote sensing inversion of geophysical elements such as cloud characteristics, aerosol, land surface characteristics, ocean water color and the like. In the NDVI index calculation based on FY-3E MERSI-2, the spatial resolution of the red light band and the near infrared band is 250m, and the corresponding wavelength ranges are 600-700 nm and 815-915 nm respectively.
Although the spatial resolution of the near infrared band and the red light band of the MERSI-2 and the MODIS are the same, the corresponding wavelength ranges are different, and the transit time in the Chinese region is different, so that the NDVI obtained by the FY-3E MERSI-2 and the NDVI obtained by the MODIS in the same region on the same date are different, and the monitoring results of the two on the earth vegetation and the ecological environment are different.
Disclosure of Invention
The invention provides a conversion method of normalized vegetation indexes of FY-3E MERSI-2 and MODIS, which realizes the conversion of the normalized vegetation indexes of FY-3E MERSI-2 and MODIS and effectively ensures the consistency of the two indexes in monitoring earth vegetation and ecological environment.
The invention provides a normalized vegetation index conversion method of FY-3E MERSI-2 and MODIS, which comprises the following steps:
acquiring a green plant reflection spectrum characteristic curve graph of FY-3E MERSI-2 and MODIS;
extracting reflection spectrum information of a resolution imaging spectrometer in MODIS, which corresponds to a wavelength of 620-670 nm, and reflection spectrum information of a corresponding wavelength of 841-876 nm from a green plant reflection spectrum characteristic curve graph, and obtaining a normalized vegetation index expression of the resolution imaging spectrometer in MODIS;
extracting reflection spectrum information with the wavelength of 600-700 nm corresponding to the resolution imaging spectrometer in the FY-3E MERSI-2 and reflection spectrum information with the wavelength of 815-915 nm corresponding to the resolution imaging spectrometer in the green plant reflection spectrum characteristic curve graph, and obtaining a normalized vegetation index expression of the resolution imaging spectrometer in the FY-3E MERSI-2;
obtaining a linear function expression of the resolution imaging spectrometer in the MODIS about wavelength and reflectivity according to a green plant reflection spectrum characteristic curve graph;
obtaining a conversion expression for converting the normalized vegetation index of MODIS into the normalized vegetation index of FY-3E MERSI-2 according to the normalized vegetation index expression of the resolution imaging spectrometer in MODIS, the normalized vegetation index expression of the resolution imaging spectrometer in FY-3E MERSI-2 and a linear function expression of the resolution imaging spectrometer in MODIS on wavelength and reflectivity, and realizing the normalized vegetation index conversion of FY-3E MERSI-2 and MODIS through the conversion expression.
According to the normalized vegetation index conversion method of FY-3E MERSI-2 and MODIS provided by the invention, the normalized vegetation index expression of the resolution imaging spectrometer in the MODIS is as follows:
NDVImod＝(B2-A2)/(B2+A2)，
wherein NDVImodThe method comprises the steps of representing a normalized vegetation index of a resolution imaging spectrometer in the MODIS, representing reflection spectrum information of the resolution imaging spectrometer in the MODIS, corresponding to the wavelength of 620-670 nm by A2, and representing reflection spectrum information of the resolution imaging spectrometer in the MODIS, corresponding to the wavelength of 841-876 nm by B2;
the normalized vegetation index expression of the FY-3E MERSI-2 medium-resolution imaging spectrometer is as follows:
NDVIFY＝[(B1+B2+B3)-(A1+A2+A3)]/[(B1+B2+B3)+(A1+A2+A3)]，
wherein NDVIFYThe method comprises the steps of representing a normalized vegetation index of a resolution imaging spectrometer in the FY-3E MERSI-2, (A1+ A2+ A3) representing reflection spectrum information of the resolution imaging spectrometer in the FY-3E MERSI-2 corresponding to the wavelength of 600-700 nm, and (B1+ B2+ B3) representing reflection spectrum information of the resolution imaging spectrometer in the FY-3E MERSI-2 corresponding to the wavelength of 815-915 nm.
According to the normalized vegetation index conversion method of FY-3E MERSI-2 and MODIS provided by the invention, a linear function expression of a resolution imaging spectrometer in MODIS about wavelength and reflectivity is obtained according to a green plant reflection spectrum characteristic curve graph, and specifically two linear function expressions are obtained, wherein the two linear function expressions are respectively as follows:
f＝(B1+B3)/(B2+A2)＝C×NDVImod+D，
g＝(A1+A3)/(B2+A2)＝F×NDVImod+G，
f represents a first linear function expression, g represents a second linear function expression, and (B1+ B3) represents the reflection spectrum information of the resolution imaging spectrometer in the MODIS corresponding to the wavelengths of 815-841 nm and 876-915 nm; (A1+ A3) represents the reflection spectrum information of the resolution imaging spectrometer in the MODIS corresponding to the wavelengths of 600-620 nm and 670-700 nm, (B2+ A2) represents the reflection spectrum information of the resolution imaging spectrometer in the MODIS corresponding to the wavelengths of 841-876 nm and 620-670 nm, and NDVImodThe normalized vegetation index of the resolution imaging spectrometer in MODIS is shown, and C, D, F, G respectively show different relation parameters.
According to the normalized vegetation index conversion method of FY-3E MERSI-2 and MODIS provided by the invention, the conversion expression for converting the normalized vegetation index of MODIS into the normalized vegetation index of FY-3E MERSI-2 is as follows:
NDVIFY＝(a×NDVImod+b)/(c×NDVImod+d)，
wherein NDVIFYIndicating normalized vegetation index, NDVI, of a resolution imaging spectrometer in FY-3E MERSI-2modThe normalized vegetation index of the resolution imaging spectrometer in MODIS is expressed, a, b, C, D represent different regression parameters, respectively, and a ═ C-F +1, b ═ D-G, C ═ C + F, D ═ 1+ D + G.
The method for converting the normalized vegetation index of FY-3E MERSI-2 and MODIS provided by the invention further comprises the following steps: actually measuring the observation data of the resolution imaging spectrometer in the MODIS and the resolution imaging spectrometer in the FY-3E MERSI-2 on the green plants, and obtaining a functional relation diagram describing the relation between the normalized vegetation index of the MODIS and the normalized vegetation index of the FY-3E MERSI-2 according to the observation data so as to verify the conversion effect of converting the normalized vegetation index of the MODIS into the normalized vegetation index of the FY-3E MERSI-2 through a conversion expression.
According to the normalized vegetation index conversion method of FY-3E MERSI-2 and MODIS provided by the invention, the information contained in the functional relation graph describing the relation between the normalized vegetation index of MODIS and the normalized vegetation index of FY-3E MERSI-2 comprises the correlation coefficient information of the fitting relation of the conversion expression and the significance level information of the conversion expression.
The invention also provides a system for converting the normalized vegetation index of FY-3E MERSI-2 and MODIS, which comprises:
the green plant reflection spectrum characteristic curve graph obtaining module is used for obtaining green plant reflection spectrum characteristic curve graphs of FY-3E MERSI-2 and MODIS;
an MODIS normalized vegetation index expression obtaining module, which is used for extracting the reflection spectrum information of the resolution imaging spectrometer in the MODIS, which corresponds to the wavelength of 620-670 nm, and the reflection spectrum information of the resolution imaging spectrometer, which corresponds to the wavelength of 841-876 nm, from the green plant reflection spectrum characteristic curve graph, and obtaining a normalized vegetation index expression of the resolution imaging spectrometer in the MODIS;
the normalized vegetation index expression obtaining module is used for extracting reflection spectrum information of a resolution imaging spectrometer with the wavelength of 600-700 nm and reflection spectrum information of a resolution imaging spectrometer with the wavelength of 815-915 nm in the FY-3E MERSI-2 from a green plant reflection spectrum characteristic curve graph and obtaining a normalized vegetation index expression of the resolution imaging spectrometer in the FY-3E MERSI-2;
the MODIS linear function expression obtaining module is used for obtaining a linear function expression of the resolution imaging spectrometer in the MODIS about wavelength and reflectivity according to the green plant reflection spectrum characteristic curve graph;
the normalized vegetation index conversion module of the FY-3E MERSI-2 and the MODIS is used for obtaining a conversion expression for converting the normalized vegetation index of the MODIS into the normalized vegetation index of the FY-3E MERSI-2 according to a normalized vegetation index expression of a resolution imaging spectrometer in the MODIS, a normalized vegetation index expression of the resolution imaging spectrometer in the FY-3EMERSI-2 and a linear function expression of the resolution imaging spectrometer in the MODIS on wavelength and reflectivity, and realizing the normalized vegetation index conversion of the FY-3E MERSI-2 and the MODIS through the conversion expression.
The invention also provides an electronic device, which comprises a memory, a processor and a computer program stored on the memory and capable of running on the processor, wherein the processor executes the program to realize the steps of the dynamic knowledge graph-based trusted attack detection method.
The present invention also provides a non-transitory computer readable storage medium having stored thereon a computer program which, when executed by a processor, implements the steps of any of the methods for normalized vegetation index conversion of FY-3E mers-2 and MODIS as described above.
The invention also provides a computer program product comprising a computer program which, when executed by a processor, implements the steps of the method for converting normalized vegetation index of any of the FY-3E mers-2 and MODIS described above.
The invention provides a method for converting normalized vegetation indexes of FY-3E MERSI-2 and MODIS, by extracting the reflection spectrum information of different wavelength ranges corresponding to the resolution imaging spectrometer in MODIS and the resolution imaging spectrometer in FY-3E MERSI-2 under the condition that the spatial resolution of the red wave band and the near infrared wave band is 250m from the green plant reflection spectrum characteristic curve graph, respectively obtaining a normalized vegetation index expression of the resolution imaging spectrometer in the MODIS and a normalized vegetation index expression of the resolution imaging spectrometer in the FY-3E MERSI-2, then obtaining a conversion expression for converting the normalized vegetation index of the MODIS into the normalized vegetation index of the FY-3E MERSI-2 by combining a linear function expression of the resolution imaging spectrometer in the MODIS on wavelength and reflectivity, to realize the normalized vegetation index conversion of FY-3E MERSI-2 and MODIS.
Although the transit time of the resolution imaging spectrometer in the MODIS and the resolution imaging spectrometer in the FY-3E MERSI-2 in the Chinese area is different when the earth vegetation and the ecological environment are monitored, the conversion of the normalized vegetation index of the FY-3E MERSI-2 and the normalized vegetation index of the MODIS can be realized through the invention, and the consistency of the two on the earth vegetation and the ecological environment monitoring is effectively ensured.
The MODIS data has incomparable advantages of large space, high timeliness, low economy and other sensors, and has become the first choice data for developing regional, national, intercontinental and global vegetation and ecological environment monitoring at present. The FY-3E satellite is emitted at 7 months and 5 days in 2021, and related research results or technical reports are not seen at present. In order to fully excavate the application potential of FY-3E satellite data in vegetation and ecological environment monitoring research work, the data of vegetation and ecological environment monitoring engineering which is used for a long time and is verified on the basis of MODIS NDVI in the global scope is fully utilized, the ground verification work of vegetation and ecological environment monitoring engineering which is established on the basis of FY-3E MERSI-2NDVI is reduced to the maximum extent, the NDVI of the FY-3E MERSI-2 and the NDVI of the MODIS are converted by the method, the FY-3E MERSI-2NDVI is used for replacing the MODIS NDVI, and the wide application of the FY-3E MERSI-2 remote sensing data can be effectively promoted.
Drawings
In order to more clearly illustrate the technical solutions of the present invention or the prior art, the drawings needed for the description of the embodiments or the prior art will be briefly described below, and it is obvious that the drawings in the following description are some embodiments of the present invention, and those skilled in the art can also obtain other drawings according to the drawings without creative efforts.
FIG. 1 is a schematic flow chart of the normalized vegetation index conversion method of FY-3E MERSI-2 and MODIS provided by the present invention.
FIG. 2 shows a graph of the reflection spectrum characteristics of green plants of FY-3E MERSI-2 and MODIS according to the present invention.
FIG. 3 illustrates a functional diagram depicting MODIS normalized vegetation index versus FY-3E MERSI-2 normalized vegetation index.
FIG. 4 shows a block diagram of a normalized vegetation index conversion system for implementing FY-3E MERSI-2 and MODIS provided by the present invention.
Fig. 5 is a schematic structural diagram of an electronic device provided in the present invention.
Detailed Description
In order to make the objects, technical solutions and advantages of the present invention clearer, the technical solutions of the present invention will be clearly and completely described below with reference to the accompanying drawings, and it is obvious that the described embodiments are some, but not all embodiments of the present invention. All other embodiments, which can be derived by a person skilled in the art from the embodiments given herein without making any creative effort, shall fall within the protection scope of the present invention.
The normalized vegetation index conversion method of FY-3E MERSI-2 and MODIS of the present invention is described below with reference to FIGS. 1-5.
A method for converting normalized vegetation index of FY-3E MERSI-2 and MODIS, as shown in fig. 1, includes:
s1: obtaining a green plant reflection spectrum characteristic curve graph of FY-3E MERSI-2 and MODIS.
Specifically, FIG. 2 shows a green plant reflection spectrum characteristic graph of FY-3E MERSI-2 and MODIS, wherein B2 represents the reflection spectrum information corresponding to the wavelength range 841-876 nm of the resolution imaging spectrometer in MODIS, (B1+ B2+ B3) represents the reflection spectrum information corresponding to the wavelength range 815-915 nm of the resolution imaging spectrometer in FY-3E MERSI-2, A2 represents the reflection spectrum information corresponding to the wavelength range 620-670 nm of the resolution imaging spectrometer in MODIS, and (A1+ A2+ A3) represents the reflection spectrum information corresponding to the wavelength range 600-700 nm of the resolution imaging spectrometer in FY-3E MERSI-2; b1 and B3 correspond to reflection spectrum information in the wavelength ranges of 815-841 nm and 876-915 nm respectively; a1 and A3 correspond to reflection spectrum information in the wavelength ranges of 600-620 nm and 670-700 nm, respectively.
S2: and extracting the reflection spectrum information of the resolution imaging spectrometer in the MODIS, which corresponds to the wavelength of 620-670 nm, and the reflection spectrum information of the resolution imaging spectrometer which corresponds to the wavelength of 841-876 nm from the green plant reflection spectrum characteristic curve graph, and obtaining a normalized vegetation index expression of the resolution imaging spectrometer in the MODIS.
Specifically, the normalized vegetation index expression of the resolution imaging spectrometer in the MODIS is as follows:
NDVImod＝(B2-A2)/(B2+A2) (1)，
wherein NDVImodThe normalized vegetation index of the resolution imaging spectrometer in MODIS is shown, A2 in MODISThe resolution imaging spectrometer corresponds to the reflection spectrum information with the wavelength of 620-670 nm, and B2 represents the reflection spectrum information with the wavelength of 841-876 nm corresponding to the resolution imaging spectrometer in the MODIS.
S3: and extracting reflection spectrum information of the resolution imaging spectrometer with the corresponding wavelength of 600-700 nm and reflection spectrum information of the corresponding wavelength of 815-915 nm in the FY-3E MERSI-2 from a green plant reflection spectrum characteristic curve graph, and obtaining a normalized vegetation index expression of the resolution imaging spectrometer with the FY-3E MERSI-2.
Specifically, the normalized vegetation index expression of the FY-3E MERSI-2 medium-resolution imaging spectrometer is as follows:
NDVIFY＝[(B1+B2+B3)-(A1+A2+A3)]/[(B1+B2+B3)+(A1+A2+A3)](2)，
wherein NDVIFYThe method comprises the steps of representing a normalized vegetation index of a resolution imaging spectrometer in the FY-3E MERSI-2, (A1+ A2+ A3) representing reflection spectrum information of the resolution imaging spectrometer in the FY-3E MERSI-2 corresponding to the wavelength of 600-700 nm, and (B1+ B2+ B3) representing reflection spectrum information of the resolution imaging spectrometer in the FY-3E MERSI-2 corresponding to the wavelength of 815-915 nm.
Further, to simplify equation (2), one may let:
f＝(B1+B3)/(B2+A2) (3)，
g＝(A1+A3)/(B2+A2) (4)，
by substituting formulae (1), (3), and (4) for formula (2), formula (2) can be simplified to:
NDVIFY＝(f+NDVImod–g)/(f+1+g) (5)。
according to the green plant reflection spectrum characteristic graphs of the formulas (1), (3) and (4) and the FY-3E MERSI-2 and MODIS, the denominators of the formulas (1), (3) and (4) are (B2+ A2) but the numerators are different, and in the green plant reflection spectrum characteristic graphs of the FY-3E MERSI-2 and MODIS, it can be seen that the reflection information of the (B1+ B3), (A1+ A3) and (B2-A2) on the same underlying surface is relatively fixed, namely the values are proportional, namely linear relation, so that the formulas (3) and (4) can be established.
S4: and obtaining a linear function expression of the resolution imaging spectrometer in the MODIS about wavelength and reflectivity according to the green plant reflection spectrum characteristic curve graph.
Specifically, according to the green plant reflection spectrum characteristic curve diagram, a linear function expression of the resolution imaging spectrometer in the MODIS about the wavelength and the reflectivity is obtained, and specifically two linear function expressions are obtained, which are respectively:
f＝(B1+B3)/(B2+A2)＝C×NDVImod+D (6)，
g＝(A1+A3)/(B2+A2)＝F×NDVImod+G (7)，
f represents a first linear function expression, g represents a second linear function expression, and (B1+ B3) represents the reflection spectrum information of the resolution imaging spectrometer in the MODIS corresponding to the wavelengths of 815-841 nm and 876-915 nm; (A1+ A3) represents the reflection spectrum information of the resolution imaging spectrometer in the MODIS corresponding to the wavelengths of 600-620 nm and 670-700 nm, (B2+ A2) represents the reflection spectrum information of the resolution imaging spectrometer in the MODIS corresponding to the wavelengths of 841-876 nm and 620-670 nm, and NDVImodThe normalized vegetation index of the resolution imaging spectrometer in MODIS is shown, and C, D, F, G respectively show different relation parameters.
Further, by substituting formula (6) and formula (7) for formula (5), it is possible to obtain:
NDVIFY＝[(1+C-F)×NDVImod+(D-G)]/[(C+F)×NDVImod+(1+D+G)](8)。
s5: obtaining a conversion expression for converting the normalized vegetation index of MODIS into the normalized vegetation index of FY-3E MERSI-2 according to the normalized vegetation index expression of the resolution imaging spectrometer in MODIS, the normalized vegetation index expression of the resolution imaging spectrometer in FY-3E MERSI-2 and a linear function expression of the resolution imaging spectrometer in MODIS on wavelength and reflectivity, and realizing the normalized vegetation index conversion of FY-3EMERSI-2 and MODIS through the conversion expression.
Specifically, the conversion expression of converting the normalized vegetation index of MODIS into the normalized vegetation index of FY-3EMERSI-2 is as follows:
NDVIFY＝(a×NDVImod+b)/(c×NDVImod+d) (9)，
wherein NDVIFYIndicating normalized vegetation index, NDVI, of a resolution imaging spectrometer in FY-3E MERSI-2modThe normalized vegetation index of the resolution imaging spectrometer in MODIS is expressed, a, b, C, D represent different regression parameters, respectively, and a ═ C-F +1, b ═ D-G, C ═ C + F, D ═ 1+ D + G.
S6: actually measuring the observation data of the resolution imaging spectrometer in the MODIS and the resolution imaging spectrometer in the FY-3E MERSI-2 on the green plants, and obtaining a functional relation diagram describing the relation between the normalized vegetation index of the MODIS and the normalized vegetation index of the FY-3E MERSI-2 according to the observation data so as to verify the conversion effect of converting the normalized vegetation index of the MODIS into the normalized vegetation index of the FY-3E MERSI-2 through a conversion expression.
Specifically, the information contained in the functional relationship diagram describing the relationship between the normalized vegetation index of MODIS and the normalized vegetation index of FY-3E MERSI-2 includes correlation coefficient information of the fitting relationship of the conversion expression and significance level information of the conversion expression.
For example, the spectral observation of the same underlying surface is carried out by adopting a FieldSpec Pro portable spectral radiometer produced by American company, the spectral observation of 6:00 is used for simulating the observation data of FY-3E MERSI-2, the spectral observation of 10:30 is used for simulating the observation data of Terra MODIS, and the observation data of 49 corn sample prescriptions and 33 wheat sample prescriptions of a certain agricultural meteorological national field scientific observation research station in 25-27 days in 2021 year and 7-31 months in 2021 and 10 grassland sample prescriptions of a certain grassland in 28-31 days are obtained; meanwhile, for forest vegetation, performing spectrum observation on the same underlying surface by adopting an unmanned aerial vehicle hyperspectral imaging system produced by American company, simulating FY-3E MERSI-2 observation data by spectrum observation of 6:00, simulating Terra MODIS observation data by spectrum observation of 10:30, obtaining 103 pixel data of the underlying surface of the forest at 28-31 days in 7 months in 2021, calculating parameters a, b, c and d in the formula (9), and calculating the result as follows:
a＝1.98，b＝-0.07，c＝0.96，d＝1.04(R2＝0.92，P<0.001，n＝195)。
meanwhile, the relationship between the normalized vegetation index describing MODIS and the normalized vegetation index of FY-3E MERSI-2 is drawnA functional relationship diagram of the system, a functional relationship diagram describing the relationship between the normalized vegetation index of MODIS and the normalized vegetation index of FY-3E MERSI-2 is shown in FIG. 3, wherein R2The correlation coefficient values representing the fitting relationship, P representing the significance level value, n representing the observation data quantity, Grassland representing Grassland, Wheat representing Wheat, Maize representing Maize, Forest representing Forest.
It can be seen that the functional relation graph describing the relation between the normalized vegetation index of MODIS and the normalized vegetation index of FY-3E MERSI-2 visually shows the distribution of observed data (observation points/observation values), the value of the relevant coefficient of the fitting relation reaches 0.92, that is, the interpretation rate reaches 92%, it is quantitatively shown that the interpretation of the observation values (points) by the conversion expression for converting the normalized vegetation index of MODIS into the normalized vegetation index of FY-3E MERSI-2 in the normalized vegetation index conversion method of FY-3E MERSI-2 of the present invention reaches 92%, indicating that the fitting relation is very good; and p represents a significance level, and p <0.001 represents that the difference between experimental data is extremely significant, and the test is meaningful, namely the conversion effect of converting the normalized vegetation index of MODIS into the normalized vegetation index of FY-3E MERSI-2 is good through the conversion expression.
The invention provides a method for converting normalized vegetation indexes of FY-3E MERSI-2 and MODIS, by extracting the reflection spectrum information of different wavelength ranges corresponding to the resolution imaging spectrometer in MODIS and the resolution imaging spectrometer in FY-3E MERSI-2 under the condition that the spatial resolution of the red wave band and the near infrared wave band is 250m from the green plant reflection spectrum characteristic curve graph, respectively obtaining a normalized vegetation index expression of the resolution imaging spectrometer in the MODIS and a normalized vegetation index expression of the resolution imaging spectrometer in the FY-3E MERSI-2, then obtaining a conversion expression for converting the normalized vegetation index of the MODIS into the normalized vegetation index of the FY-3E MERSI-2 by combining a linear function expression of the resolution imaging spectrometer in the MODIS on wavelength and reflectivity, to realize the normalized vegetation index conversion of FY-3E MERSI-2 and MODIS.
Although the transit time of the resolution imaging spectrometer in the MODIS and the resolution imaging spectrometer in the FY-3E MERSI-2 in the Chinese area is different when the earth vegetation and the ecological environment are monitored, the conversion of the normalized vegetation index of the FY-3E MERSI-2 and the normalized vegetation index of the MODIS can be realized through the invention, and the consistency of the two on the earth vegetation and the ecological environment monitoring is effectively ensured.
The MODIS data has incomparable advantages of large space, high timeliness, low economy and other sensors, and has become the first choice data for developing regional, national, intercontinental and global vegetation and ecological environment monitoring at present. The FY-3E satellite is emitted at 7 months and 5 days in 2021, and related research results or technical reports are not seen at present. In order to fully excavate the application potential of FY-3E satellite data in vegetation and ecological environment monitoring research work, the data of vegetation and ecological environment monitoring engineering which is used for a long time and is verified on the basis of MODIS NDVI in the global scope is fully utilized, the ground verification work of vegetation and ecological environment monitoring engineering which is established on the basis of FY-3E MERSI-2NDVI is reduced to the maximum extent, the NDVI of the FY-3E MERSI-2 and the NDVI of the MODIS are converted by the method, the FY-3E MERSI-2NDVI is used for replacing the MODIS NDVI, and the wide application of the FY-3E MERSI-2 remote sensing data can be effectively promoted.
In addition, it is to be noted that the method of the present invention is also applicable to the conversion between FY-3E MERSI-2NDVI and Aqua MODIS NDVI, the conversion between FY-3D MERSI-2NDVI and Terra MODIS NDVI, the conversion between FY-3D MERSI-2NDVI and Aqua MODIS NDVI, and the like.
The following describes a system, an apparatus, a non-transitory computer-readable storage medium, and a computer program product for converting the normalized vegetation index of FY-3E MERSI-2 and MODIS provided by the present invention, and the system, the apparatus, the non-transitory computer-readable storage medium, and the computer program product for converting the normalized vegetation index of FY-3E MERSI-2 and MODIS described below and the above-described method for converting the normalized vegetation index of FY-3E MERSI-2 and MODIS may be referred to correspondingly.
The invention also discloses a normalized vegetation index conversion system of FY-3E MERSI-2 and MODIS, as shown in FIG. 4, comprising:
and a green plant reflection spectrum characteristic curve graph obtaining module 410, configured to obtain green plant reflection spectrum characteristic curves of FY-3E mers-2 and MODIS.
The MODIS normalized vegetation index expression obtaining module 420 is used for extracting reflection spectrum information of the resolution imaging spectrometer in the MODIS, corresponding to the wavelength of 620-670 nm, and reflection spectrum information of the resolution imaging spectrometer in the MODIS, corresponding to the wavelength of 841-876 nm, in a green plant reflection spectrum characteristic curve graph, and obtaining a normalized vegetation index expression of the resolution imaging spectrometer in the MODIS.
Specifically, the normalized vegetation index expression of the resolution imaging spectrometer in the MODIS is as follows:
NDVImod＝(B2-A2)/(B2+A2)，
wherein NDVImodThe normalized vegetation index of the resolution imaging spectrometer in the MODIS is represented, A2 represents the reflection spectrum information of the resolution imaging spectrometer in the MODIS, corresponding to the wavelengths of 620-670 nm, and B2 represents the reflection spectrum information of the resolution imaging spectrometer in the MODIS, corresponding to the wavelengths of 841-876 nm.
The FY-3E MERSI-2 normalized vegetation index expression obtaining module 430 is used for extracting reflection spectrum information of 600-700 nm corresponding to the wavelength of the resolution imaging spectrometer and reflection spectrum information of 815-915 nm corresponding to the wavelength in the FY-3E MERSI-2 from a green plant reflection spectrum characteristic curve graph and obtaining a normalized vegetation index expression of the resolution imaging spectrometer in the FY-3E MERSI-2.
Specifically, the normalized vegetation index expression of the FY-3E MERSI-2 medium-resolution imaging spectrometer is as follows:
NDVIFY＝[(B1+B2+B3)-(A1+A2+A3)]/[(B1+B2+B3)+(A1+A2+A3)]，
wherein NDVIFYThe method comprises the steps of representing a normalized vegetation index of a resolution imaging spectrometer in the FY-3E MERSI-2, (A1+ A2+ A3) representing reflection spectrum information of the resolution imaging spectrometer in the FY-3E MERSI-2 corresponding to the wavelength of 600-700 nm, and (B1+ B2+ B3) representing reflection spectrum information of the resolution imaging spectrometer in the FY-3E MERSI-2 corresponding to the wavelength of 815-915 nm.
The MODIS linear function expression obtaining module 440 is configured to obtain a linear function expression of the resolution imaging spectrometer in MODIS with respect to wavelength and reflectivity according to the green plant reflection spectrum characteristic curve.
Specifically, according to the green plant reflection spectrum characteristic curve diagram, a linear function expression of the resolution imaging spectrometer in the MODIS about the wavelength and the reflectivity is obtained, and specifically two linear function expressions are obtained, which are respectively:
f＝(B1+B3)/(B2+A2)＝C×NDVImod+D，
g＝(A1+A3)/(B2+A2)＝F×NDVImod+G，
f represents a first linear function expression, g represents a second linear function expression, and (B1+ B3) represents reflection spectrum information of the resolution imaging spectrometer in the MODIS corresponding to the wavelengths of 815-841 nm and 876-915 nm; (A1+ A3) represents the reflection spectrum information of the resolution imaging spectrometer in the MODIS corresponding to the wavelengths of 600-620 nm and 670-700 nm, (B2+ A2) represents the reflection spectrum information of the resolution imaging spectrometer in the MODIS corresponding to the wavelengths of 841-876 nm and 620-670 nm, and NDVImodThe normalized vegetation index of the resolution imaging spectrometer in MODIS is shown, and C, D, F, G respectively show different relation parameters.
The conversion module 450 for the normalized vegetation index of FY-3E MERSI-2 and MODIS is configured to obtain a conversion expression for converting the normalized vegetation index of MODIS into the normalized vegetation index of FY-3E MERSI-2 according to the normalized vegetation index expression of the resolution imaging spectrometer in MODIS, the normalized vegetation index expression of the resolution imaging spectrometer in FY-3E MERSI-2, and a linear function expression of the resolution imaging spectrometer in MODIS with respect to wavelength and reflectivity, and implement the conversion of the normalized vegetation index of FY-3E MERSI-2 and MODIS through the conversion expression.
Specifically, the conversion expression for converting the normalized vegetation index of MODIS into the normalized vegetation index of FY-3E MERSI-2 is as follows:
NDVIFY＝(a×NDVImod+b)/(c×NDVImod+d)，
wherein NDVIFYIndicating normalized vegetation index, NDVI, of a resolution imaging spectrometer in FY-3E MERSI-2modNormalized vegetation finger for representing resolution imaging spectrometer in MODISThe numbers a, b, C, D represent different regression parameters, and a ═ C-F +1, b ═ D-G, C ═ C + F, D ═ 1+ D + G, respectively.
And the conversion expression inspection module is used for actually measuring the observation data of the resolution imaging spectrometer in the MODIS and the resolution imaging spectrometer in the FY-3E MERSI-2 on the green plants, and obtaining a functional relation diagram describing the relation between the normalized vegetation index of the MODIS and the normalized vegetation index of the FY-3E MERSI-2 according to the observation data so as to inspect the conversion effect of converting the normalized vegetation index of the MODIS into the normalized vegetation index of the FY-3E MERSI-2 through a conversion expression.
Fig. 5 illustrates a physical structure diagram of an electronic device, which may include, as shown in fig. 5: a processor (processor)810, a communication Interface 820, a memory 830 and a communication bus 840, wherein the processor 810, the communication Interface 820 and the memory 830 communicate with each other via the communication bus 840. The processor 810 may invoke logic instructions in the memory 830 to perform a method of normalized vegetation index conversion of FY-3E MERSI-2 and MODIS, the method comprising:
acquiring a green plant reflection spectrum characteristic curve graph of FY-3E MERSI-2 and MODIS;
extracting reflection spectrum information of a resolution imaging spectrometer in MODIS, which corresponds to a wavelength of 620-670 nm, and reflection spectrum information of a corresponding wavelength of 841-876 nm from a green plant reflection spectrum characteristic curve graph, and obtaining a normalized vegetation index expression of the resolution imaging spectrometer in MODIS;
extracting reflection spectrum information of a resolution imaging spectrometer with the corresponding wavelength of 600-700 nm and reflection spectrum information of the corresponding wavelength of 815-915 nm in the FY-3E MERSI-2 from a green plant reflection spectrum characteristic curve graph, and obtaining a normalized vegetation index expression of the resolution imaging spectrometer in the FY-3E MERSI-2;
obtaining a linear function expression of the resolution imaging spectrometer in the MODIS about wavelength and reflectivity according to a green plant reflection spectrum characteristic curve graph;
obtaining a conversion expression for converting the normalized vegetation index of MODIS into the normalized vegetation index of FY-3E MERSI-2 according to the normalized vegetation index expression of the resolution imaging spectrometer in MODIS, the normalized vegetation index expression of the resolution imaging spectrometer in FY-3E MERSI-2 and a linear function expression of the resolution imaging spectrometer in MODIS on wavelength and reflectivity, and realizing the normalized vegetation index conversion of FY-3E MERSI-2 and MODIS through the conversion expression.
In addition, the logic instructions in the memory 830 may be implemented in software functional units and stored in a computer readable storage medium when the logic instructions are sold or used as independent products. Based on such understanding, the technical solution of the present invention may be embodied in the form of a software product, which is stored in a storage medium and includes instructions for causing a computer device (which may be a personal computer, a server, or a network device) to execute all or part of the steps of the method according to the embodiments of the present invention. And the aforementioned storage medium includes: a U-disk, a removable hard disk, a Read-Only Memory (ROM), a Random Access Memory (RAM), a magnetic disk, or an optical disk, and various media capable of storing program codes.
In another aspect, the present invention further provides a computer program product, the computer program product including a computer program, the computer program being stored on a non-transitory computer readable storage medium, and when the computer program is executed by a processor, the computer being capable of executing the normalized vegetation index conversion method of FY-3E mers-2 and MODIS provided by the above methods, the method including:
acquiring a green plant reflection spectrum characteristic curve graph of FY-3E MERSI-2 and MODIS;
extracting reflection spectrum information of a resolution imaging spectrometer in MODIS, which corresponds to a wavelength of 620-670 nm, and reflection spectrum information of a corresponding wavelength of 841-876 nm from a green plant reflection spectrum characteristic curve graph, and obtaining a normalized vegetation index expression of the resolution imaging spectrometer in MODIS;
extracting reflection spectrum information of a resolution imaging spectrometer with the corresponding wavelength of 600-700 nm and reflection spectrum information of the corresponding wavelength of 815-915 nm in the FY-3E MERSI-2 from a green plant reflection spectrum characteristic curve graph, and obtaining a normalized vegetation index expression of the resolution imaging spectrometer in the FY-3E MERSI-2;
obtaining a linear function expression of the resolution imaging spectrometer in the MODIS about wavelength and reflectivity according to a green plant reflection spectrum characteristic curve graph;
obtaining a conversion expression for converting the normalized vegetation index of MODIS into the normalized vegetation index of FY-3E MERSI-2 according to the normalized vegetation index expression of the resolution imaging spectrometer in MODIS, the normalized vegetation index expression of the resolution imaging spectrometer in FY-3E MERSI-2 and a linear function expression of the resolution imaging spectrometer in MODIS on wavelength and reflectivity, and realizing the normalized vegetation index conversion of FY-3E MERSI-2 and MODIS through the conversion expression.
In another aspect, the present invention further provides a non-transitory computer-readable storage medium, on which a computer program is stored, the computer program being implemented by a processor to perform the method for converting normalized vegetation index of FY-3E measure-2 and MODIS provided by the above methods, the method comprising:
acquiring a green plant reflection spectrum characteristic curve graph of FY-3E MERSI-2 and MODIS;
extracting reflection spectrum information of a resolution imaging spectrometer in MODIS, which corresponds to a wavelength of 620-670 nm, and reflection spectrum information of a corresponding wavelength of 841-876 nm from a green plant reflection spectrum characteristic curve graph, and obtaining a normalized vegetation index expression of the resolution imaging spectrometer in MODIS;
extracting reflection spectrum information of a resolution imaging spectrometer with the corresponding wavelength of 600-700 nm and reflection spectrum information of the corresponding wavelength of 815-915 nm in the FY-3E MERSI-2 from a green plant reflection spectrum characteristic curve graph, and obtaining a normalized vegetation index expression of the resolution imaging spectrometer in the FY-3E MERSI-2;
obtaining a linear function expression of the resolution imaging spectrometer in the MODIS about wavelength and reflectivity according to a green plant reflection spectrum characteristic curve graph;
obtaining a conversion expression for converting the normalized vegetation index of MODIS into the normalized vegetation index of FY-3E MERSI-2 according to the normalized vegetation index expression of the resolution imaging spectrometer in MODIS, the normalized vegetation index expression of the resolution imaging spectrometer in FY-3E MERSI-2 and a linear function expression of the resolution imaging spectrometer in MODIS on wavelength and reflectivity, and realizing the normalized vegetation index conversion of FY-3E MERSI-2 and MODIS through the conversion expression.
The above-described embodiments of the apparatus are merely illustrative, and the units described as separate parts may or may not be physically separate, and parts displayed as units may or may not be physical units, may be located in one place, or may be distributed on a plurality of network units. Some or all of the modules may be selected according to actual needs to achieve the purpose of the solution of the present embodiment. One of ordinary skill in the art can understand and implement it without inventive effort.
Through the above description of the embodiments, those skilled in the art will clearly understand that each embodiment can be implemented by software plus a necessary general hardware platform, and certainly can also be implemented by hardware. With this understanding in mind, the above-described technical solutions may be embodied in the form of a software product, which can be stored in a computer-readable storage medium such as ROM/RAM, magnetic disk, optical disk, etc., and includes instructions for causing a computer device (which may be a personal computer, a server, or a network device, etc.) to execute the methods described in the embodiments or some parts of the embodiments.
Finally, it should be noted that: the above examples are only intended to illustrate the technical solution of the present invention, but not to limit it; although the present invention has been described in detail with reference to the foregoing embodiments, it will be understood by those of ordinary skill in the art that: the technical solutions described in the foregoing embodiments may still be modified, or some technical features may be equivalently replaced; and such modifications or substitutions do not depart from the spirit and scope of the corresponding technical solutions of the embodiments of the present invention.

Claims (7)

1. A normalized vegetation index conversion method of FY-3E MERSI-2 and MODIS is characterized by comprising the following steps:
acquiring a green plant reflection spectrum characteristic curve graph of FY-3E MERSI-2 and MODIS;
extracting reflection spectrum information of a resolution imaging spectrometer in MODIS, which corresponds to a wavelength of 620-670 nm, and reflection spectrum information of a corresponding wavelength of 841-876 nm from a green plant reflection spectrum characteristic curve graph, and obtaining a normalized vegetation index expression of the resolution imaging spectrometer in MODIS;
extracting reflection spectrum information of a resolution imaging spectrometer with the corresponding wavelength of 600-700 nm and reflection spectrum information of the corresponding wavelength of 815-915 nm in the FY-3E MERSI-2 from a green plant reflection spectrum characteristic curve graph, and obtaining a normalized vegetation index expression of the resolution imaging spectrometer in the FY-3E MERSI-2;
obtaining a linear function expression of the resolution imaging spectrometer in the MODIS about wavelength and reflectivity according to a green plant reflection spectrum characteristic curve graph;
obtaining a conversion expression for converting the normalized vegetation index of MODIS into the normalized vegetation index of FY-3E MERSI-2 according to a normalized vegetation index expression of a resolution imaging spectrometer in MODIS, a normalized vegetation index expression of the resolution imaging spectrometer in FY-3E MERSI-2 and a linear function expression of the resolution imaging spectrometer in MODIS on wavelength and reflectivity, and realizing the normalized vegetation index conversion of FY-3E MERSI-2 and MODIS through the conversion expression;
the normalized vegetation index expression of the resolution imaging spectrometer in the MODIS is as follows:
NDVImod＝(B2-A2)/(B2+A2)，
wherein NDVImodNormalized vegetation index for resolution imaging spectrometer in MODISA2 represents the reflection spectrum information of the resolution imaging spectrometer in the MODIS with the wavelength of 620-670 nm, and B2 represents the reflection spectrum information of the resolution imaging spectrometer in the MODIS with the wavelength of 841-876 nm;
the normalized vegetation index expression of the FY-3E MERSI-2 medium-resolution imaging spectrometer is as follows:
NDVIFY＝[(B1+B2+B3)-(A1+A2+A3)]/[(B1+B2+B3)+(A1+A2+A3)]，
wherein NDVIFYThe method comprises the steps of (A) representing a normalized vegetation index of a resolution imaging spectrometer in FY-3E MERSI-2, (A1+ A2+ A3) representing reflection spectrum information of the resolution imaging spectrometer in the FY-3E MERSI-2 corresponding to the wavelength of 600-700 nm, and (B1+ B2+ B3) representing reflection spectrum information of the resolution imaging spectrometer in the FY-3E MERSI-2 corresponding to the wavelength of 815-915 nm;
according to the green plant reflection spectrum characteristic curve graph, obtaining a linear function expression of the resolution imaging spectrometer in the MODIS about wavelength and reflectivity, and specifically obtaining two linear function expressions, wherein the two linear function expressions are respectively as follows:
f＝(B1+B3)/(B2+A2)＝C×NDVImod+D，
g＝(A1+A3)/(B2+A2)＝F×NDVImod+G，
f represents a first linear function expression, g represents a second linear function expression, and (B1+ B3) represents the reflection spectrum information of the resolution imaging spectrometer in the MODIS corresponding to the wavelengths of 815-841 nm and 876-915 nm; (A1+ A3) represents the reflection spectrum information of the resolution imaging spectrometer in the MODIS corresponding to the wavelengths of 600-620 nm and 670-700 nm, (B2+ A2) represents the reflection spectrum information of the resolution imaging spectrometer in the MODIS corresponding to the wavelengths of 841-876 nm and 620-670 nm, and NDVImodThe normalized vegetation index of the resolution imaging spectrometer in MODIS is shown, and C, D, F, G respectively show different relation parameters.
2. The method of claim 1, wherein the conversion expression for converting the normalized vegetation index of MODIS into the normalized vegetation index of FY-3E MERSI-2 is as follows:
NDVIFY＝(a×NDVImod+b)/(c×NDVImod+d)，
wherein NDVIFYIndicating normalized vegetation index, NDVI, of a resolution imaging spectrometer in FY-3E MERSI-2modThe normalized vegetation index of the resolution imaging spectrometer in MODIS is expressed, a, b, C, D represent different regression parameters, respectively, and a ═ C-F +1, b ═ D-G, C ═ C + F, D ═ 1+ D + G.
3. The method of claim 1, further comprising the steps of: actually measuring the observation data of the resolution imaging spectrometer in the MODIS and the resolution imaging spectrometer in the FY-3E MERSI-2 on the green plants, and obtaining a functional relation diagram describing the relation between the normalized vegetation index of the MODIS and the normalized vegetation index of the FY-3E MERSI-2 according to the observation data so as to verify the conversion effect of converting the normalized vegetation index of the MODIS into the normalized vegetation index of the FY-3E MERSI-2 through a conversion expression.
4. The method of claim 1, wherein the functional relationship diagram describing the relationship between the normalized vegetation index of MODIS and the normalized vegetation index of FY-3E MERSI-2 includes information including correlation coefficient information of the fitting relationship of the conversion expression and significance level information of the conversion expression.
5. A normalized vegetation index conversion system of FY-3E MERSI-2 and MODIS, comprising:
the green plant reflection spectrum characteristic curve graph obtaining module is used for obtaining green plant reflection spectrum characteristic curve graphs of FY-3E MERSI-2 and MODIS;
an MODIS normalized vegetation index expression obtaining module, which is used for extracting reflection spectrum information of a resolution imaging spectrometer in MODIS, corresponding to the wavelength of 620-670 nm and reflection spectrum information of the corresponding wavelength of 841-876 nm, from a green plant reflection spectrum characteristic curve graph and obtaining a normalized vegetation index expression of the resolution imaging spectrometer in MODIS;
the normalized vegetation index expression obtaining module is used for extracting reflection spectrum information of a resolution imaging spectrometer with the wavelength of 600-700 nm and reflection spectrum information of a resolution imaging spectrometer with the wavelength of 815-915 nm in the FY-3E MERSI-2 from a green plant reflection spectrum characteristic curve graph and obtaining a normalized vegetation index expression of the resolution imaging spectrometer in the FY-3E MERSI-2;
the MODIS linear function expression obtaining module is used for obtaining a linear function expression of the resolution imaging spectrometer in the MODIS about wavelength and reflectivity according to the green plant reflection spectrum characteristic curve graph;
the normalized vegetation index conversion module of FY-3E MERSI-2 and MODIS is used for obtaining a conversion expression for converting the normalized vegetation index of MODIS into the normalized vegetation index of FY-3E MERSI-2 according to the normalized vegetation index expression of the resolution imaging spectrometer in MODIS, the normalized vegetation index expression of the resolution imaging spectrometer in FY-3E MERSI-2 and the linear function expression of the resolution imaging spectrometer in MODIS on wavelength and reflectivity, and realizing the normalized vegetation index conversion of FY-3E MERSI-2 and MODIS through the conversion expression;
the normalized vegetation index expression of the resolution imaging spectrometer in the MODIS is as follows:
NDVImod＝(B2-A2)/(B2+A2)，
wherein NDVImodThe method comprises the steps of representing a normalized vegetation index of a resolution imaging spectrometer in the MODIS, representing reflection spectrum information of the resolution imaging spectrometer in the MODIS, corresponding to the wavelength of 620-670 nm by A2, and representing reflection spectrum information of the resolution imaging spectrometer in the MODIS, corresponding to the wavelength of 841-876 nm by B2;
the normalized vegetation index expression of the FY-3E MERSI-2 medium-resolution imaging spectrometer is as follows:
NDVIFY＝[(B1+B2+B3)-(A1+A2+A3)]/[(B1+B2+B3)+(A1+A2+A3)]，
wherein NDVIFYExpressing the normalized vegetation index of the resolution imaging spectrometer in the FY-3E MERSI-2, (A1+ A2+ A3) expressing that the resolution imaging spectrometer in the FY-3E MERSI-2 corresponds toReflection spectrum information with the wavelength of 600-700 nm, (B1+ B2+ B3) represents the reflection spectrum information with the wavelength of 815-915 nm corresponding to a resolution imaging spectrometer in FY-3E MERSI-2;
according to the green plant reflection spectrum characteristic curve graph, obtaining a linear function expression of the resolution imaging spectrometer in the MODIS about wavelength and reflectivity, and specifically obtaining two linear function expressions, wherein the two linear function expressions are respectively as follows:
f＝(B1+B3)/(B2+A2)＝C×NDVImod+D，
g＝(A1+A3)/(B2+A2)＝F×NDVImod+G，
f represents a first linear function expression, g represents a second linear function expression, and (B1+ B3) represents the reflection spectrum information of the resolution imaging spectrometer in the MODIS corresponding to the wavelengths of 815-841 nm and 876-915 nm; (A1+ A3) represents the reflection spectrum information of the resolution imaging spectrometer in the MODIS corresponding to the wavelengths of 600-620 nm and 670-700 nm, (B2+ A2) represents the reflection spectrum information of the resolution imaging spectrometer in the MODIS corresponding to the wavelengths of 841-876 nm and 620-670 nm, and NDVImodThe normalized vegetation index of the resolution imaging spectrometer in MODIS is shown, and C, D, F, G respectively show different relation parameters.
6. An electronic device comprising a memory, a processor and a computer program stored on the memory and executable on the processor, wherein the processor when executing the program implements the steps of the method of normalized vegetation index conversion of FY-3E measure si-2 and MODIS as claimed in any one of claims 1 to 4.
7. A non-transitory computer readable storage medium having stored thereon a computer program, wherein the computer program when executed by a processor implements the steps of the method for normalized vegetation index conversion of FY-3E measure si-2 and MODIS as claimed in any one of claims 1 to 4.
CN202111205935.XA 2021-10-15 2021-10-15 Normalized vegetation index conversion method of MERSI-2 and MODIS Active CN113989661B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202111205935.XA CN113989661B (en) 2021-10-15 2021-10-15 Normalized vegetation index conversion method of MERSI-2 and MODIS

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202111205935.XA CN113989661B (en) 2021-10-15 2021-10-15 Normalized vegetation index conversion method of MERSI-2 and MODIS

Publications (2)

Publication Number Publication Date
CN113989661A CN113989661A (en) 2022-01-28
CN113989661B true CN113989661B (en) 2022-05-17

Family

ID=79738909

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202111205935.XA Active CN113989661B (en) 2021-10-15 2021-10-15 Normalized vegetation index conversion method of MERSI-2 and MODIS

Country Status (1)

Country Link
CN (1) CN113989661B (en)

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108955883A (en) * 2018-07-19 2018-12-07 国家卫星气象中心 Remote sensor Calibration Method and scaling system based on high-spectral data
CN109753916A (en) * 2018-12-28 2019-05-14 厦门理工学院 A kind of vegetation index spatial scaling model building method and device
CN111175231A (en) * 2019-12-31 2020-05-19 国家卫星气象中心(国家空间天气监测预警中心) Inversion method and device of canopy vegetation index and server
CN112881327A (en) * 2021-01-25 2021-06-01 安徽皖南烟叶有限责任公司 Tobacco leaf SPAD value estimation method based on novel vegetation index
AU2021105579A4 (en) * 2021-08-16 2021-10-14 Institute Of Agricultural Resources And Regional Planning, Chinese Academy Of Agricultural Sciences Method for Retrieving Land Surface Temperature from FY-3D/MERSI-2 Data

Family Cites Families (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103345707A (en) * 2013-06-04 2013-10-09 中国科学院遥感与数字地球研究所 Crop maturation stage remote sensing prediction method based on multi-source remote sensing data
CN106547840B (en) * 2016-10-13 2019-09-06 国家卫星气象中心 A kind of parsing and management method of whole world three-dimensional atmospheric data
CN111402162B (en) * 2020-03-13 2023-11-07 北京华云星地通科技有限公司 Method for processing clear sky data set of satellite remote sensing image
CN111982837A (en) * 2020-08-27 2020-11-24 中国气象科学研究院 Conversion method of vegetation ecological parameter remote sensing estimation model

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108955883A (en) * 2018-07-19 2018-12-07 国家卫星气象中心 Remote sensor Calibration Method and scaling system based on high-spectral data
CN109753916A (en) * 2018-12-28 2019-05-14 厦门理工学院 A kind of vegetation index spatial scaling model building method and device
CN111175231A (en) * 2019-12-31 2020-05-19 国家卫星气象中心(国家空间天气监测预警中心) Inversion method and device of canopy vegetation index and server
CN112881327A (en) * 2021-01-25 2021-06-01 安徽皖南烟叶有限责任公司 Tobacco leaf SPAD value estimation method based on novel vegetation index
AU2021105579A4 (en) * 2021-08-16 2021-10-14 Institute Of Agricultural Resources And Regional Planning, Chinese Academy Of Agricultural Sciences Method for Retrieving Land Surface Temperature from FY-3D/MERSI-2 Data

Also Published As

Publication number Publication date
CN113989661A (en) 2022-01-28

Similar Documents

Publication Publication Date Title
Yang et al. Modeling grassland above-ground biomass based on artificial neural network and remote sensing in the Three-River Headwaters Region
Rodríguez-Veiga et al. Quantifying forest biomass carbon stocks from space
Sinha et al. Estimation of leaf area index using PROSAIL based LUT inversion, MLRA-GPR and empirical models: Case study of tropical deciduous forest plantation, North India
CN111678866B (en) Soil water content inversion method for multi-model ensemble learning
CN109974665B (en) Aerosol remote sensing inversion method and system for short-wave infrared data lack
Sánchez et al. On the synergy of airborne GNSS-R and Landsat 8 for soil moisture estimation
Schweiger et al. Plant beta-diversity across biomes captured by imaging spectroscopy
Carrer et al. Land surface albedo derived on a ten daily basis from Meteosat Second Generation Observations: The NRT and climate data record collections from the EUMETSAT LSA SAF
Gara et al. Leaf to canopy upscaling approach affects the estimation of canopy traits
Fu et al. Spatial downscaling of MODIS Chlorophyll-a using Landsat 8 images for complex coastal water monitoring
CN114740180A (en) Soil organic carbon estimation method and device based on multi-source remote sensing data
Zhu et al. UAV-based hyperspectral analysis and spectral indices constructing for quantitatively monitoring leaf nitrogen content of winter wheat
CN115494007A (en) Random forest based high-precision rapid detection method and device for soil organic matters
Wu et al. Winter wheat LAI inversion considering morphological characteristics at different growth stages coupled with microwave scattering model and canopy simulation model
dos Santos et al. Vegetation cover monitoring in tropical regions using SAR-C dual-polarization index: seasonal and spatial influences
CN110836870A (en) GEE-based large-area lake transparency rapid drawing method
Weiß et al. Evaluation of different radiative transfer models for microwave backscatter estimation of wheat fields
Huang et al. Constructing a finer-resolution forest height in China using icesat/glas, landsat and Alos Palsar data and height patterns of natural forests and plantations
Park et al. Decadal measurements of the first Geostationary Ocean Color Satellite (GOCI) compared with MODIS and VIIRS data
CN111079835B (en) Himapari-8 atmospheric aerosol inversion method based on deep full-connection network
Lou et al. An effective method for canopy chlorophyll content estimation of marsh vegetation based on multiscale remote sensing data
CN114241331B (en) Remote sensing modeling method for ground biomass of reed in wetland by taking UAV as ground and Septinel-2 medium
Kaiser et al. Potential of Sentinel-1 data for spatially and temporally high-resolution detection of drought affected forest stands
Zahran et al. Remote sensing based water resources and agriculture spatial indicators system
CN113989661B (en) Normalized vegetation index conversion method of MERSI-2 and MODIS

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
GR01 Patent grant
GR01 Patent grant