CN102103077A - MODIS data-based agricultural drought monitoring method - Google Patents
MODIS data-based agricultural drought monitoring method Download PDFInfo
- Publication number
- CN102103077A CN102103077A CN 200910248457 CN200910248457A CN102103077A CN 102103077 A CN102103077 A CN 102103077A CN 200910248457 CN200910248457 CN 200910248457 CN 200910248457 A CN200910248457 A CN 200910248457A CN 102103077 A CN102103077 A CN 102103077A
- Authority
- CN
- China
- Prior art keywords
- index
- pairing
- pixel point
- vswi
- monitored
- 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
Landscapes
- Image Processing (AREA)
Abstract
The invention discloses a moderate-resolution imaging spectroradiometer (MODIS) data-based agricultural drought monitoring method. In the farmland drought monitoring process, an agricultural drought index is determined by a crop water supply index and a rainfall anomaly index; in the drought monitoring process, a vegetation index and the surface temperature are inverted by MODIS data, and the crop water supply index is calculated by the vegetation index and the surface temperature; the rainfall anomaly index is calculated by rainfall data; and drought severity is determined by dividing the level of the drought index.
Description
Technical field
The present invention relates to a kind of agricultural drought monitoring method.
Background technology
The generating process of arid is potential, is not easy to find; The occurrence characteristics of agricultural drought is that coverage is big, brings serious catastrophic effect and economic loss; Research, estimate that arid takes place and the process of development, can take the corresponding drought resisting measure of preventing and reducing natural disasters, reduce the agricultural disaster loss.Along with the development of remote sensing technology, remote sensing is with its dynamic, real-time, multispectral, cheap advantage, for new approach has been opened up in the damage caused by a drought monitoring.Vegetation index that remote sensing is obtained and surface temperature are to describe two ten minutes important parameters of earth surface feature, when therefore arid takes place, can disclose the plant physiology off-note by the variation of vegetation index or surface temperature, reflect farmland hydro-thermal stress state indirectly.Can be divided into according to the employed wave band of remotely-sensed data: visible light, near infrared, thermal infrared, microwave etc.Because employed wave band difference has produced numerous models and method.As water deficit exponential model, temperature vegetation drought index model etc.Though it is a lot of to be used for the model of agricultural drought disaster monitoring, great majority are experimental research.Therefore it is very necessary to explore a kind of accurate relatively agricultural drought monitoring method.
Summary of the invention
The purpose of this invention is to provide a kind of new method of utilizing remote sensing technology to carry out face of land draught monitor, this agricultural drought monitoring utilizes crop water supply exponential sum rainfall anomaly index to determine agriculture damage caused by a drought index.When carrying out the damage caused by a drought monitoring, utilize MODIS data inversion vegetation index and surface temperature, utilize vegetation index and surface temperature to calculate crop water supply index.Utilize precipitation data computation precipitation anomaly index.At last the damage caused by a drought index is carried out rank and divide to determine the arid order of severity.。
The present invention adopts and monitors based on the agricultural arid monitoring model of vegetation index and surface temperature, in the process of monitoring the model Ecological Parameter is improved, the relevant parameter of the face of land damage caused by a drought of explaining is determined as follows: based on the spectral reflectance data of each wave band pixel of the pairing MODIS of monitored farmland massif;
A kind of agricultural arid monitoring method based on the MODIS data mainly is to inquire into a kind of objective, dynamic, real-time, accurate, monitoring method of being easy to realize by the remote sensing technology means on big regional scale.Adopt in the method based on the agricultural arid monitoring model of vegetation index and surface temperature and monitor, in the process of monitoring the model Ecological Parameter is improved, the relevant parameter of the face of land damage caused by a drought of explaining is determined as follows:
Spectral reflectance data based on each wave band pixel of the pairing MODIS of monitored farmland massif;
Obtain EVI, T
s, VSWI and SDI parameter, remotely-sensed data to selected wave band is carried out unified pre-service, mainly comprise according to a conventional method and carry out (in remote sensing softwares such as ERDAS or ENVI, handling) radiant quantity calculating, geometry correction, cloud detection, projection conversion, so that further mate calculating; When calculating rainfall anomaly index, at first the rainfall data should be carried out interpolation processing, obtain precipitation anomaly figure; The unified raster data that changes into after interpolation is finished, so as with the remotely-sensed data comparative analysis.
1) calculating of vegetation index:
EVI is for strengthening vegetation index, ρ
NIR, ρ
RedAnd ρ
BlueBe respectively the spectral reflectivity of near-infrared band, red spectral band and the blue wave band pixel of the pairing MODIS of monitored farmland massif; L is that background is adjusted item; C
1And C
2Be fitting coefficient; G is a gain factor; When calculating MODIS-EVI, L=1, C
1=6, C
2=7.5, G=2.5; Calculate the EVI data of pairing each the pixel point of remote sensing image of monitored farmland massif respectively;
2) Surface Temperature Retrieval
Surface temperature is a basic parameter of agricultural drought disaster monitoring, its computing method adopt (Qin Zhihao such as Qin Zhi person of outstanding talent, Gao Maofang, Qin Xiaomin etc., surface temperature remote sensing inversion method in the agricultural drought disaster monitoring-with MODIS data instance [J]. the disaster journal, 2005,14 (4): the 64-71.) window algorithm of splitting based on middle large scale MODIS data of Ti Chuing, its computing formula is:
T
s=A
0+A
1T
31-A
2T
32
A
0=E
1a
31-E
2a
32
A
1=1+A+E
1b
31
A
2=A+E
2b
32
A=D
31/E
0
E
1=D
32(1-C
31-D
31)/E
0
E
2=D
31(1-C
32-D
32)/E
0
E
0=D
32C
31-D
31C
32
C
t=ε
iτ
i
D
t=[1+(1-ε
i)τ
i]
T in the formula
sBe surface temperature (K), T
31And T
32Be respectively the brightness temperature of MODIS the 31st and 32 wave bands, A
0, A
1And A
2Be the parameter of splitting window algorithm, a
31, b
31, a
31And b
32Be constant, the desirable a of difference in surface temperature 0-50 ℃ scope
31=-64.60363, b
31=0.440817, a
32=-68.72575, b
32=0.473453; The calculating of other intermediate parameters sees the division window algorithm (Qin Zhihao at the MODIS data such as Qin Zhihao for details, Gao Maofang, Qin Xiaomin etc., the surface temperature remote sensing inversion method in the agricultural drought disaster monitoring-with MODIS data instance [J]. the disaster journal, 2005,14 (4): 64-71.); Utilize formula to calculate the T of pairing each the pixel point of remote sensing image of monitored farmland massif respectively then
sData;
Wherein i is meant the 31st and 32 wave bands of the pairing MODIS image of monitored farmland massif, is respectively i=31 or 32; τ
iBe the atmospheric transmittance of the wave band i of pairing each the pixel point of remote sensing image of monitored farmland massif, ε
iIt is the face of land emissivity of the wave band i of pairing each the pixel point of remote sensing image of monitored farmland massif.
ε
t=ε
iw+P
vR
vε
iv+(1-P
v)R
sε
is
ε
Iw, ε
IvAnd ε
IsBe respectively water body, vegetation and the exposed soil of pairing each the pixel point of remote sensing image of monitored farmland massif face of land emissivity, get ε respectively at wave band i wave band
31w=0.99683, ε
32w=0.99254, ε
31v=0.98672, ε
32v=0.98990, ε
31s=0.96767, ε
31s=0.97790; R
vAnd R
sBe respectively the vegetation of the pairing remote sensing image of monitored farmland massif and the radiation ratio of exposed soil, (Qin Z such as Qin Zhihao, Karnieli A.Progress in the Remote Sensing of Land SurfaceTemperature and Ground Emissivity Using NOAA-AVHRR[J] .InternationalJournal of Remote Sensing, 1999,20 (12): 2367-2397.), the R of calculating
v=0.99240, R
s=1.00744.P
vBe the vegetation coverage of pairing each the pixel point of monitored farmland massif, estimate by vegetation index:
NDVI is the vegetation index of pairing each the pixel point of monitored farmland massif in the formula, NDVI
vAnd NDVI
sBe respectively the NDVI value of dense vegetation covering and complete exposed soil pixel, get NDVI usually
v=0.674, NDVI
s=0.039 (Sun Rui, Zhu Qijiang. the analysis [J] of the clean primary productivity of clean primary productivity model of vegetation and China. Beijing Normal University's journal (natural science edition), 1998,34:132-137.);
B in the formula
1And B
2It is respectively the reflectivity of MODIS image the 1st and 2 wave bands.
Present near linear relationship between atmospheric transmittance and the atmosphere vapour content, according to Qin Zhihao etc. (
Qin?Z,Karnieli?A,Berliner?P.A?Mono- window?Algorithm?for?RetrievingLand?Surface?Temperature?from?Landsat?TM?Data?and?Its?Application
When moisture content at 0.4~2.0g/cm
2The time:
τ
31=0.99513-0.0808w
τ
32=0.99376-0.11369w
When moisture content at 2~4.0g/cm
2The time:
τ
31=1.08692-0.12759w
τ
32=1.07900-0.15925w
When moisture content at 4.0~6.0g/cm
2The time:
τ
31=1.07268-0.12571w
τ
32=0.93821-0.12613w
The formula that calculates atmosphere vapour content is:
w=[(α-lnT
w)/β]
2
Wherein w is the atmosphere vapour content of pairing each the pixel point of monitored farmland massif; T
wBe the ratio of the ground surface reflectance of the reflectivity of pairing each the pixel point of monitored farmland massif MODIS image the 18th wave band and pairing each the pixel point of the 2nd wave band, α, β are parameters, get α=0.02 respectively, β=0.651;
3) calculating of crop water supply index:
VSWI is the crop water supply index of pairing each the pixel point of monitored farmland massif in the formula; EVI is the enhancing vegetation index of each pixel point of the pairing MODIS image of monitored farmland massif; T
sSurface temperature for pairing each the pixel point of remote sensing image of monitored farmland massif;
4) crop water supply index is carried out standardization:
SDI=(VSWI-VSWI
d)/(VSWI
w-VSWI
d)×100%
SDI is the water supply index of each pixel point crop of the pairing remote sensing image of monitored farmland massif after the standardization in the formula, gets 0~100%, and wherein SDI=0 represents severe drought, and SDI=100% represents very moistening; VSWI
dAnd VSWI
wWhen being respectively the non-irrigateiest and the crop water supply index when the most moistening; The classification step-length of EVI can be made as d, when d=0.05, and when the temperature space of suitable plant growth is 20 ℃~45 ℃, VSWI
d=(n * d)/45, VSWI
w=(n * d)/20, n is for strengthening the number of vegetation index step-length, the positive integer of n 〉=1.For farmland ecosystem, VAWI
dAnd VSWI
wValue under the temperature of suitable growth is respectively:
5) calculating of precipitation anomaly index
SRI is the precipitation anomaly index of monitored each pixel point of farmland massif in the formula, and SRI is big more moistening more; R is for working as the ten days rainfall amount; R
wWhen ten days rainfall amount mean value, get nearest 10 years for for many years in the method when ten days rainfall amount mean value.As R>2R
wThe time for extremely moistening, get SRI=100%, R>R
wThe time be normal, get SRI=50%, R is less than an average for many years half, is arid extremely, this moment, the SRI value was 0~25%.
Arid is a rainless meteorologic phenomena of long period, and a ten days rainlessly arid might not occur with interior, and the rainless of above time in a ten days just can produce arid problem, therefore the rainfall amount that needs to consider simultaneously early stage.This method has been considered the precipitation affects in nearest 8 ten days, thereby calculates comprehensive precipitation anomaly index, and its formula is:
SMRI=A
0×SRI
0+A
1×SRI
1+A
2×SRI
2+A
3×SRI
3+…+A
8×SRI
8
6) determine agriculture damage caused by a drought index
DI=B
1×SDI+B
2×SMRI
DI is the agriculture damage caused by a drought index of pairing each the pixel point of remote sensing image of monitored farmland massif in the formula, and it is the coupling of crop water supply exponential sum rainfall anomaly index, gets 0~100%, and DI=0 represents very arid, and DI=100% represents very moistening; SDI is a standardization water supply index, B
1Be its weight, get 0.6; SMRI is a drought index of considering rain factor of many ten days, B
2Be its weight, get 0.4; B
1, B
2The weight value be that value according to SDI conforms to the most with actual conditions.According to the rank criteria for classifying, judge the degree of arid then;
Herein, 1%≤DI≤15% drought of attaching most importance to, 15%<DI≤30% is middle drought, and 30%<DI≤50% is light drought, and 50%<DI≤70% be that normally 70%<DI≤100% is moistening.
The invention has the advantages that:
Utilize the MODIS remotely-sensed data that draught monitor is carried out in the farmland in big zone.In monitoring, model parameter is improved, adopted enhancing vegetation index (EVI) to replace normalized differential vegetation index (NDVI) that crop water supply index is calculated.NDVI and EVI are the vegetation indexs that Moderate Imaging Spectroradiomete data (MODIS) are selected for use, and EVI improves as the succession of NDVI with to some defective of NDVI, overall treatment Soil Background, atmospheric noise and ruddiness saturation problem.And the calculating of NDVI is the non-linear stretching near infrared and red spectral band, consequently strengthened low value part, suppressed high value part, and it is saturated that signal is appearred in high vegetation-covered area easily.It is easily saturated that EVI has then overcome the high areal coverage of vegetation, and the low areal coverage of vegetation is subjected to soil vegetative cover to influence bigger shortcoming.The present invention utilizes remote sensing technology to carry out the agricultural drought dynamic monitoring, and through practical application check, this method is easy, efficient, easy operating, result are accurate, can be widely used among the agricultural drought monitoring of China;
Embodiment
With the monitoring of North China's agricultural drought is example.Study area mainly comprises Hebei, Beijing, Tianjin, Shandong, Henan, Jiangsu and Anhui Province.This district is the main farming region of China, and major weather disaster is an arid, its occurrence frequency height, and longer duration seriously influences national economy.What local area often occurred is spring drought; Therefore the damage caused by a drought in North China April has been carried out study on monitoring.
At first select corresponding MODISL1B data according to the study area scope, (in remotely-sensed data softwares such as ERDAS or ENVI, handle) according to a conventional method the remotely-sensed data of selected wave band is carried out unified pre-service, utilize the remote sensing software data to carry out geometry correction, cloud detection and projection conversion, so that further mate calculating.Utilize each corresponding wave band data to calculate EVI, the T of each pixel point of monitored area respectively then
sWith crop water supply index.And select its pairing VSWI according to the size of the EVI value of monitored district each pixel point of remote sensing image
dAnd VSWI
wValue is carried out standardization to crop water supply index.Utilize the comprehensive rainfall anomaly of the rainfall data computation precipitation anomaly exponential sum index in the upper, middle and lower ten days in April of each meteorological site of study area; In ArcGIS software,, obtain and to delete trrellis diagram with the precipitation anomaly index of Remote Sensing Images Matching comprehensive precipitation anomaly exponential interpolation.At last according to the standardized crop water supply of precipitation anomaly exponential sum Index for Calculation agricultural damage caused by a drought index.
A kind of agricultural arid monitoring method based on the MODIS data mainly is to inquire into a kind of objective, dynamic, real-time, accurate, monitoring method of being easy to realize by the remote sensing technology means on big regional scale.Adopt in the method based on the agricultural arid monitoring model of vegetation index and surface temperature and monitor, in the process of monitoring the model Ecological Parameter is improved, the relevant parameter of the face of land damage caused by a drought of explaining is determined as follows:
Spectral reflectance data based on each wave band pixel of the pairing MODIS of monitored farmland massif;
1) calculating of vegetation index:
EVI is for strengthening vegetation index, ρ
NIR, ρ
RedAnd ρ
BlueBe respectively the spectral reflectivity of near-infrared band, red spectral band and the blue wave band pixel of the pairing MODIS of monitored farmland massif; L is that background is adjusted item; C
1And C
2Be fitting coefficient; G is a gain factor; When calculating MODIS-EVI, L=1, C
1=6, C
2=7.5, G=2.5; Calculate the EVI data of pairing each the pixel point of remote sensing image of monitored farmland massif respectively;
2) Surface Temperature Retrieval
Surface temperature is a basic parameter of agricultural drought disaster monitoring, its computing method adopt (Qin Zhihao such as Qin Zhi person of outstanding talent, Gao Maofang, Qin Xiaomin etc., surface temperature remote sensing inversion method in the agricultural drought disaster monitoring-with MODIS data instance [J]. the disaster journal, 2005,14 (4): the 64-71.) window algorithm of splitting based on middle large scale MODIS data of Ti Chuing, its computing formula is:
T
s=A
0+A
1T
31-A
2T
32
A
0=E
1a
31-E
2a
32
A
1=1+A+E
1b
31
A
2=A+E
2b
32
A=D
31/E
0
E
1=D
32(1-C
31-D
31)/E
0
E
2=D
31(1-C
32-D
32)/E
0
E
0=D
32C
31-D
31C
32
C
i=ε
iτ
i
D
t=[1+(1-ε
i)τ
i]
T in the formula
sBe surface temperature (K), T
31And T
32Be respectively the brightness temperature of MODIS the 31st and 32 wave bands, A
0, A
1And A
2Be the parameter of splitting window algorithm, a
31, b
31, a
31And b
32Be constant, the desirable a of difference in surface temperature 0-50 ℃ scope
31=-64.60363, b
31=0.440817, a
32=-68.72575, b
32=0.473453; The calculating of other intermediate parameters sees the division window algorithm (Qin Zhihao at the MODIS data such as Qin Zhihao for details, Gao Maofang, Qin Xiaomin etc., the surface temperature remote sensing inversion method in the agricultural drought disaster monitoring-with MODIS data instance [J]. the disaster journal, 2005,14 (4): 64-71.); Utilize formula to calculate the T of pairing each the pixel point of remote sensing image of monitored farmland massif respectively then
sData;
Wherein i is meant the 31st and 32 wave bands of the pairing MODIS image of monitored farmland massif, is respectively i=31 or 32; τ
iBe the atmospheric transmittance of the wave band i of pairing each the pixel point of remote sensing image of monitored farmland massif, ε
iIt is the face of land emissivity of the wave band i of pairing each the pixel point of remote sensing image of monitored farmland massif.
ε
i=ε
iw+P
vR
vε
iv+(1-P
v)R
sε
is
ε
Iw, ε
IvAnd ε
IsBe respectively water body, vegetation and the exposed soil of pairing each the pixel point of remote sensing image of monitored farmland massif face of land emissivity, get ε respectively at wave band i wave band
31w=0.99683, ε
32w=0.99254, ε
31v=0.98672, ε
32v=0.98990, ε
31s=0.96767, ε
31s=0.97790; R
vAnd R
sBe respectively the vegetation of the pairing remote sensing image of monitored farmland massif and the radiation ratio of exposed soil, (Qin Z such as Qin Zhihao, Karnieli A.Progress in the Remote Sensing of Land SurfaceTemperature and Ground Emissivity Using NOAA-AVHRR[J] .InternationalJournal of Remote Sensing, 1999,20 (12): 2367-2397.), the R of calculating
v=0.99240, R
s1.00744.P
vBe the vegetation coverage of pairing each the pixel point of monitored farmland massif, estimate by vegetation index:
NDVI is the vegetation index of pairing each the pixel point of monitored farmland massif in the formula, NDVI
vAnd NDVI
sBe respectively the NDVI value of dense vegetation covering and complete exposed soil pixel, get NDVI usually
v=0.674, NDVI
s=0.039 (Sun Rui, Zhu Qijiang. the analysis [J] of the clean primary productivity of clean primary productivity model of vegetation and China. Beijing Normal University's journal (natural science edition), 1998,34:132-137.);
B in the formula
1And B
2It is respectively the reflectivity of MODIS image the 1st and 2 wave bands.
Present near linear relationship between atmospheric transmittance and the atmosphere vapour content, according to Qin Zhihao etc. (
Qin?Z,Karnieli?A,Berliner?P.A?Mono- window?Algorithm?for?RetrievingLand?Surface?Temperature?from?Landsat?TM?Data?and?Its?Application?to
When moisture content at 0.4~2.0g/cm
2The time:
τ
31=0.99513-0.0808w
τ
32=0.99376-0.11369w
When moisture content at 2~4.0g/cm
2The time:
τ
31=1.08692-0.12759w
τ
32=1.07900-0.15925w
When moisture content at 4.0~6.0g/cm
2The time:
τ
31=1.07268-0.12571w
τ
32=0.93821-0.12613w
The formula that calculates atmosphere vapour content is:
w=[(α-lnT
w)/β]
2
Wherein w is the atmosphere vapour content of pairing each the pixel point of monitored farmland massif; T
wBe the ratio of the ground surface reflectance of the reflectivity of pairing each the pixel point of monitored farmland massif MODIS image the 18th wave band and pairing each the pixel point of the 2nd wave band, α, β are parameters, get α=0.02 respectively, β=0.651;
3) calculating of crop water supply index:
VSWI is the crop water supply index of pairing each the pixel point of monitored farmland massif in the formula; EVI is the enhancing vegetation index of each pixel point of the pairing MODIS image of monitored farmland massif; T
sSurface temperature for pairing each the pixel point of remote sensing image of monitored farmland massif;
4) crop water supply index is carried out standardization:
SDI=(VSWI-VSWI
d)/(VSWI
w-VSWI
d)×100%
SDI is the water supply index of each pixel point crop of the pairing remote sensing image of monitored farmland massif after the standardization in the formula, gets 0~100%, and wherein SDI=0 represents severe drought, and SDI=100% represents very moistening; VSWI
dAnd VSWI
wWhen being respectively the non-irrigateiest and the crop water supply index when the most moistening; The classification step-length of EVI can be made as d, when d=0.05, and when the temperature space of suitable plant growth is 20 ℃~45 ℃, VSWI
d=(n * d)/45, VSWI
w=(n * d)/20, n be for strengthening the number of vegetation index step-length, the positive integer of n 〉=1 (in this example n determine that by every square kilometre of monitoring farmland massif area be 1 pixel, n is the pixel number).For farmland ecosystem, VAWI
dAnd VSWI
wValue under the temperature of suitable growth is respectively:
5) calculating of precipitation anomaly index
SRI is the precipitation anomaly index of monitored each pixel point of farmland massif in the formula, and SRI is big more moistening more; R is for working as the ten days rainfall amount; R
wWhen ten days rainfall amount mean value, get nearest 10 years for for many years in the method when ten days rainfall amount mean value.As R>2R
wThe time for extremely moistening, get SRI=100%, R>R
wThe time be normal, get SRI=50%, R is less than an average for many years half, is arid extremely, this moment, the SRI value was 0~25%.
Arid is a rainless meteorologic phenomena of long period, and a ten days rainlessly arid might not occur with interior, and the rainless of above time in a ten days just can produce arid problem, therefore the rainfall amount that needs to consider simultaneously early stage.This method has been considered the precipitation affects in nearest 8 ten days, thereby calculates comprehensive precipitation anomaly index, and its formula is:
SMRI=A
0×SRI
0+A
1×SRI
1+A
2×SRI
2+A
3×SRI
3+…+A
8×SRI
8
6) determine agriculture damage caused by a drought index
DI=B
1×SDI+B
2×SMRI
DI is the agriculture damage caused by a drought index of pairing each the pixel point of remote sensing image of monitored farmland massif in the formula, and it is the coupling of crop water supply exponential sum rainfall anomaly index, gets 0~100%, and DI=0 represents very arid, and DI=100% represents very moistening; SDI is a standardization water supply index, B
1Be its weight, get 0.6; SMRI is a drought index of considering rain factor of many ten days, B
2Be its weight, get 0.4; B
1, B
2The weight value be that value according to SDI conforms to the most with actual conditions.According to the rank criteria for classifying, judge the degree of arid then;
Herein, 1%≤DI≤15% drought of attaching most importance to, 15%<DI≤30% is middle drought, and 30%<DI≤50% is light drought, and 50%<DI≤70% be that normally 70%<DI≤100% is moistening.
In actual applications,, need add up area suffered from drought by calculating each province, the different damage caused by a drought grades in district and the ratio of the suffering from drought evaluation degree (table 1) of suffering from drought to the damage caused by a drought result in order to provide the damage caused by a drought result to relevant department more intuitively.
Table 1: each province in April, North China agricultural damage caused by a drought monitoring result
From this month precipitation event, western part, this month Hebei and big the precipitation in Henan are more on the low side, because cold and heat air is all very active, so strong convective weather is very many, the middle ten days and the last ten days, continuous thunderstorm weather has appearred in some places in Jiangsu, Anhui.But find out roughly that by the website observation data the general not enough 20mm of most of precipitation region adds that antecedent precipitation is on the low side, strong wind weather is many, and temperature significantly gos up, and causes the soil loss of moisture very fast, the high temperature and little rainfall that continue make the arid sustainable development in these places, to the end of month drought-hit area enlarge to some extent.From the rainfall data, this monitoring method can be good at reflecting the basic trend of agriculture damage caused by a drought.
Claims (2)
1. agricultural arid monitoring method based on the MODIS data is based on the spectral reflectance data of each wave band pixel of the pairing MODIS of monitored farmland massif; It is characterized in that:
1) calculating of vegetation index:
EVI is for strengthening vegetation index, ρ
NIR, ρ
RedAnd ρ
BlueBe respectively the spectral reflectivity of near-infrared band, red spectral band and the blue wave band pixel of the pairing MODIS of monitored farmland massif; L is that background is adjusted item; C
1And C
2Be fitting coefficient; G is a gain factor; When calculating MODIS-EVI, L=1, C
1=6, C
2=7.5, G=2.5; Calculate the EVI data of pairing each the pixel point of remote sensing image of monitored farmland massif respectively;
2) Surface Temperature Retrieval, its computing formula is:
T
s=A
0+A
1T
31-A
2T
32
A
0=E
1a
31-E
2a
32
A
1=1+A+E
1b
31
A
2=A+E
2b
32
A=D
31/E
0
E
1=D
32(1-C
31-D
31)/E
0
E
2=D
31(1-C
32-D
32)/E
0
E
0=D
32C
31-D
31C
32
C
i=ε
iτ
i
D
i=[1+(1-ε
i)τ
i]
T in the formula
sBe surface temperature (K), T
31And T
32Be respectively the brightness temperature of MODIS the 31st and 32 wave bands, A
0, A
1And A
2Be the parameter of splitting window algorithm, a
31, b
31, a
31And b
32Be constant, the desirable a of difference in surface temperature 0-50 ℃ scope
31=-64.60363, b
31=0.440817, a
32=-68.72575, b
32=0.473453; Utilize formula to calculate the T of pairing each the pixel point of remote sensing image of monitored farmland massif respectively then
sData;
Wherein i is meant the 31st and 32 wave bands of the pairing MODIS image of monitored farmland massif, is respectively i=31 or 32; τ
iBe the atmospheric transmittance of the wave band i of pairing each the pixel point of remote sensing image of monitored farmland massif, ε
iIt is the face of land emissivity of the wave band i of pairing each the pixel point of remote sensing image of monitored farmland massif.
ε
i=ε
iw+P
vR
vε
iv+(1-P
v)R
sε
is
ε
Iw, ε
IvAnd ε
IsThe water body, vegetation and the exposed soil that are pairing each the pixel point of remote sensing image of monitored farmland massif are respectively got ε respectively at the face of land of i wave band emissivity
31w=0.99683, ε
32w=0.99254, ε
31v=0.98672, ε
32v=0.98990, ε
31s=0.96767, ε
31s=0.97790; R
vAnd R
sBe respectively the vegetation of the pairing remote sensing image of monitored farmland massif and the radiation ratio of exposed soil, the R of calculating
v=0.99240, R
s=1.00744.P
vBe the vegetation coverage of pairing each the pixel point of monitored farmland massif, estimate by vegetation index:
NDVI is the vegetation index of pairing each the pixel point of monitored farmland massif in the formula, NDVI
vAnd NDVI
sBe respectively the NDVI value of dense vegetation covering and complete exposed soil pixel, get NDVI usually
v=0.674, NDVI
s=0.039;
B1 and B2 are respectively the reflectivity of MODIS image the 1st and 2 wave bands in the formula.
Present near linear relationship between atmospheric transmittance and the atmosphere vapour content, relational expression is
When moisture content at 0.4~2.0g/cm
2The time:
τ
31=0.99513-0.0808w
τ
32=0.99376-0.11369w
When moisture content at 2~4.0g/cm
2The time:
τ
31=1.08692-0.12759w
τ
32=1.07900-0.15925w
When moisture content at 4.0~6.0g/cm
2The time:
τ
31=1.07268-0.12571w
τ
32=0.93821-0.12613w
The formula that calculates atmosphere vapour content is:
w=[(α-lnT
w)/β]
2
Wherein w is the atmosphere vapour content of pairing each the pixel point of monitored farmland massif; T
wBe the ratio of the ground surface reflectance of the reflectivity of pairing each the pixel point of monitored farmland massif MODIS image the 18th wave band and pairing each the pixel point of the 2nd wave band, α, β are parameters, get α=0.02 respectively, β=0.651;
3) calculating of crop water supply index:
VSWI is the crop water supply index of pairing each the pixel point of monitored farmland massif in the formula; EVI is the enhancing vegetation index of each pixel point of the pairing MODIS image of monitored farmland massif; T
sSurface temperature for pairing each the pixel point of remote sensing image of monitored farmland massif;
4) crop water supply index is carried out standardization:
SDI=(VSWI-VSWI
d)/(VSWI
w-VSWI
d)×100%
SDI is the water supply index of each pixel point crop of the pairing remote sensing image of monitored farmland massif after the standardization in the formula, gets 0~100%, and wherein SDI=0 represents severe drought, and SDI=100% represents very moistening; VSWI
dAnd VSWI
wWhen being respectively the non-irrigateiest and the crop water supply index when the most moistening; The classification step-length of EVI can be made as d, when d=0.05, and when the temperature space of suitable plant growth is 20 ℃~45 ℃, VSWI
d=(n * d)/45, VSWI
w=(n * d)/20, n is for strengthening the number of vegetation index step-length, the positive integer of n 〉=1.For farmland ecosystem, VAWI
dAnd VSWI
wValue under the temperature of suitable growth is respectively:
5) calculating of precipitation anomaly index
SRI is the precipitation anomaly index of monitored each pixel point of farmland massif in the formula, and SRI is big more moistening more; R is for working as the ten days rainfall amount; R
wWhen ten days rainfall amount mean value, get nearest 10 years for for many years in the method when ten days rainfall amount mean value.As R>2R
wThe time for extremely moistening, get SRI=100%, R>R
wThe time be normal, get SRI=50%, R is less than an average for many years half, is arid extremely, this moment, the SRI value was 0~25%.
Arid is a rainless meteorologic phenomena of long period, and a ten days rainlessly arid might not occur with interior, and the rainless of above time in a ten days just can produce arid problem, therefore the rainfall amount that needs to consider simultaneously early stage.This method has been considered the precipitation affects in nearest 8 ten days, thereby calculates comprehensive precipitation anomaly index, and its formula is:
SMRI=A
0×SRI
0+A
1×SRI
1+A
2×SRI
2+A
3×SRI
3+…+A
8×SRI
8
6) determine agriculture damage caused by a drought index
DI=B
1×SDI+B
2×SMRI
DI is the agriculture damage caused by a drought index of pairing each the pixel point of remote sensing image of monitored farmland massif in the formula, and it is the coupling of crop water supply exponential sum rainfall anomaly index, gets 0~100%, and DI=0 represents very arid, and DI=100% represents very moistening; SDI is a standardization water supply index, B
1Be its weight, get 0.6; SMRI is a drought index of considering rain factor of many ten days, B
2Be its weight, get 0.4; B
1, B
2The weight value be that value according to SDI conforms to the most with actual conditions.According to the rank criteria for classifying, judge the degree of arid then;
Herein, 1%≤DI≤15% drought of attaching most importance to, 15%<DI≤30% is middle drought, and 30%<DI≤50% is light drought, and 50%<DI≤70% be that normally 70%<DI≤100% is moistening.
2. according to the described method of claim 1, it is characterized in that: obtain EVI, T
s, VSWI and SDI parameter, the remotely-sensed data of selected wave band is carried out unified pre-service, mainly comprise and carry out radiant quantity calculating, geometry correction, cloud detection, projection conversion according to a conventional method, so that further mate calculating; When calculating rainfall anomaly index, at first the rainfall data should be carried out interpolation processing, obtain precipitation anomaly figure; The unified raster data that changes into after interpolation is finished, so as with the remotely-sensed data comparative analysis.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN 200910248457 CN102103077A (en) | 2009-12-16 | 2009-12-16 | MODIS data-based agricultural drought monitoring method |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN 200910248457 CN102103077A (en) | 2009-12-16 | 2009-12-16 | MODIS data-based agricultural drought monitoring method |
Publications (1)
Publication Number | Publication Date |
---|---|
CN102103077A true CN102103077A (en) | 2011-06-22 |
Family
ID=44156009
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN 200910248457 Pending CN102103077A (en) | 2009-12-16 | 2009-12-16 | MODIS data-based agricultural drought monitoring method |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102103077A (en) |
Cited By (32)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103238493A (en) * | 2013-05-20 | 2013-08-14 | 江苏省农业科学院 | Responses of different peach stock types for persistent drought and evaluation method of drought resistance of different peach stock types |
CN103424405A (en) * | 2013-08-27 | 2013-12-04 | 华南农业大学 | Drought monitoring method based on HJ-1A/1B CCD data |
CN103439297A (en) * | 2013-09-05 | 2013-12-11 | 太原理工大学 | Remote sensing estimation method for fresh weights of green plants in desert grassland |
CN103678914A (en) * | 2013-12-16 | 2014-03-26 | 中国科学院遥感与数字地球研究所 | Alpine grassland soil respiration estimation method based on satellite remote sensing data |
CN103761674A (en) * | 2014-01-27 | 2014-04-30 | 林兴志 | Crop growing period alarming and intervening method based on remote sensing and mass climate information |
CN103902839A (en) * | 2014-04-18 | 2014-07-02 | 山东科技大学 | Single channel method surface temperature retrieval method supported by satellite remote sensing product |
CN103994976A (en) * | 2013-11-28 | 2014-08-20 | 江苏省水利科学研究院 | MODIS data-based agricultural drought remote sensing monitoring method |
CN104484576A (en) * | 2014-12-31 | 2015-04-01 | 中国科学院深圳先进技术研究院 | Method for monitoring meteorological drought |
CN105023072A (en) * | 2015-08-19 | 2015-11-04 | 苏州奥诺遥感科技有限公司 | Multi-drought index fusion method based on structure inference |
CN105628618A (en) * | 2015-12-21 | 2016-06-01 | 中国矿业大学(北京) | Method for detecting earth surface evapotranspiration based on remote sensing earth surface temperature and vegetation coverage two-stage space |
CN105760814A (en) * | 2016-01-25 | 2016-07-13 | 中国水利水电科学研究院 | Data mining-based drought monitoring method |
CN105929406A (en) * | 2016-04-25 | 2016-09-07 | 珠江水利委员会珠江水利科学研究院 | Agricultural drought remote sensing monitoring method |
CN106226260A (en) * | 2016-08-10 | 2016-12-14 | 武汉大学 | A kind of combination microwave and the Soil Moisture Inversion method of infrared remote sensing image |
CN106525753A (en) * | 2016-10-12 | 2017-03-22 | 天津大学 | Convenient and simple remote-sensing soil moisture monitoring method |
CN106650238A (en) * | 2016-11-18 | 2017-05-10 | 北京师范大学 | Remote-sensing-based watershed vegetation ecological system dry season water consumption estimation method |
CN107077650A (en) * | 2014-09-12 | 2017-08-18 | 克莱米特公司 | Method and system for managing rural activity |
CN107085712A (en) * | 2017-04-28 | 2017-08-22 | 山东省农业可持续发展研究所 | A kind of agricultural arid monitoring method based on MODIS data |
CN107219226A (en) * | 2017-07-24 | 2017-09-29 | 中国科学院遥感与数字地球研究所 | Image collecting device and enhancing vegetation index monitoring system |
CN108154193A (en) * | 2018-01-16 | 2018-06-12 | 黄河水利委员会黄河水利科学研究院 | A kind of long-term sequence precipitation data NO emissions reduction method |
CN108169161A (en) * | 2017-12-12 | 2018-06-15 | 武汉大学 | A kind of corn planting regional soil humidity appraisal procedure based on modified MODIS indexes |
CN108896188A (en) * | 2018-06-06 | 2018-11-27 | 北京空间飞行器总体设计部 | A kind of high rail remote sensing satellite infrared channel atmospheric correction method of Shared aperture |
CN108982548A (en) * | 2018-07-20 | 2018-12-11 | 浙江大学 | A kind of soil moisture inversion method based on passive microwave remote sensing data |
CN108985959A (en) * | 2018-08-09 | 2018-12-11 | 安徽大学 | Wheat powdery mildew remote sensing monitoring method based on surface temperature inversion technology |
CN109063553A (en) * | 2018-06-22 | 2018-12-21 | 中国矿业大学 | Field-crop growth defect area's remote sensing fast diagnosis method after a kind of land control |
CN109186774A (en) * | 2018-08-30 | 2019-01-11 | 清华大学 | Surface temperature information acquisition method, device, computer equipment and storage medium |
CN110111379A (en) * | 2019-05-06 | 2019-08-09 | 新疆农业大学 | A kind of modified surface temperature daily difference NO emissions reduction evaluation method based on MODIS Yu Landsat data |
CN110163472A (en) * | 2019-04-11 | 2019-08-23 | 中国水利水电科学研究院 | A wide range of extreme drought emergency monitoring and impact evaluation method and system |
CN110175932A (en) * | 2019-05-16 | 2019-08-27 | 中国农业科学院农业资源与农业区划研究所 | The Agriculture Drought monitoring method coupled with more ten days precipitation is extracted based on remote sensing information |
CN110243409A (en) * | 2019-06-18 | 2019-09-17 | 中国农业科学院农业资源与农业区划研究所 | A kind of eco-drought monitoring and forecasting system and method based on earth's surface water-heat process |
CN110263409A (en) * | 2019-06-14 | 2019-09-20 | 中国农业科学院农业资源与农业区划研究所 | A kind of staple crops Low temperature disaster remote sensing monitoring forecast system and method |
CN110287944A (en) * | 2019-07-04 | 2019-09-27 | 吉林大学 | The crop pests monitoring method of multi-spectrum remote sensing image based on deep learning |
CN117093813A (en) * | 2023-09-06 | 2023-11-21 | 中国科学院地理科学与资源研究所 | Composite drought index calculation method |
-
2009
- 2009-12-16 CN CN 200910248457 patent/CN102103077A/en active Pending
Cited By (49)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103238493A (en) * | 2013-05-20 | 2013-08-14 | 江苏省农业科学院 | Responses of different peach stock types for persistent drought and evaluation method of drought resistance of different peach stock types |
CN103238493B (en) * | 2013-05-20 | 2015-02-25 | 江苏省农业科学院 | Responses of different peach stock types for persistent drought and evaluation method of drought resistance of different peach stock types |
CN103424405B (en) * | 2013-08-27 | 2015-04-29 | 华南农业大学 | Drought monitoring method based on HJ-1A/1B CCD data |
CN103424405A (en) * | 2013-08-27 | 2013-12-04 | 华南农业大学 | Drought monitoring method based on HJ-1A/1B CCD data |
CN103439297A (en) * | 2013-09-05 | 2013-12-11 | 太原理工大学 | Remote sensing estimation method for fresh weights of green plants in desert grassland |
CN103439297B (en) * | 2013-09-05 | 2016-03-23 | 太原理工大学 | A kind of desert steppe green plants fresh weight remote sensing estimation method |
CN103994976B (en) * | 2013-11-28 | 2016-07-06 | 江苏省水利科学研究院 | Agriculture Drought remote-sensing monitoring method based on MODIS data |
CN103994976A (en) * | 2013-11-28 | 2014-08-20 | 江苏省水利科学研究院 | MODIS data-based agricultural drought remote sensing monitoring method |
CN103678914A (en) * | 2013-12-16 | 2014-03-26 | 中国科学院遥感与数字地球研究所 | Alpine grassland soil respiration estimation method based on satellite remote sensing data |
CN103761674A (en) * | 2014-01-27 | 2014-04-30 | 林兴志 | Crop growing period alarming and intervening method based on remote sensing and mass climate information |
CN103761674B (en) * | 2014-01-27 | 2016-11-23 | 林兴志 | The crop growth phase based on remote sensing with magnanimity climatic information alerts and interference method |
CN103902839A (en) * | 2014-04-18 | 2014-07-02 | 山东科技大学 | Single channel method surface temperature retrieval method supported by satellite remote sensing product |
CN107077650A (en) * | 2014-09-12 | 2017-08-18 | 克莱米特公司 | Method and system for managing rural activity |
CN104484576A (en) * | 2014-12-31 | 2015-04-01 | 中国科学院深圳先进技术研究院 | Method for monitoring meteorological drought |
CN104484576B (en) * | 2014-12-31 | 2018-04-13 | 中国科学院深圳先进技术研究院 | A kind of method for monitoring meteorological drought |
CN105023072A (en) * | 2015-08-19 | 2015-11-04 | 苏州奥诺遥感科技有限公司 | Multi-drought index fusion method based on structure inference |
CN105023072B (en) * | 2015-08-19 | 2019-01-18 | 嘉兴市南湖区翊轩塑料五金厂(普通合伙) | A kind of more drought index fusion methods based on structure inferring |
CN105628618A (en) * | 2015-12-21 | 2016-06-01 | 中国矿业大学(北京) | Method for detecting earth surface evapotranspiration based on remote sensing earth surface temperature and vegetation coverage two-stage space |
CN105628618B (en) * | 2015-12-21 | 2018-03-09 | 中国矿业大学(北京) | Method based on Remote Sensing temperature Yu vegetation cover degree two benches space exploration Remote sensing |
CN105760814A (en) * | 2016-01-25 | 2016-07-13 | 中国水利水电科学研究院 | Data mining-based drought monitoring method |
CN105760814B (en) * | 2016-01-25 | 2019-07-02 | 中国水利水电科学研究院 | A kind of drought monitoring method based on data mining |
CN105929406B (en) * | 2016-04-25 | 2018-11-02 | 珠江水利委员会珠江水利科学研究院 | A kind of Agriculture Drought remote-sensing monitoring method |
CN105929406A (en) * | 2016-04-25 | 2016-09-07 | 珠江水利委员会珠江水利科学研究院 | Agricultural drought remote sensing monitoring method |
CN106226260A (en) * | 2016-08-10 | 2016-12-14 | 武汉大学 | A kind of combination microwave and the Soil Moisture Inversion method of infrared remote sensing image |
CN106226260B (en) * | 2016-08-10 | 2019-04-12 | 武汉大学 | A kind of Soil Moisture Inversion method of combination microwave and infrared remote sensing image |
CN106525753A (en) * | 2016-10-12 | 2017-03-22 | 天津大学 | Convenient and simple remote-sensing soil moisture monitoring method |
CN106650238A (en) * | 2016-11-18 | 2017-05-10 | 北京师范大学 | Remote-sensing-based watershed vegetation ecological system dry season water consumption estimation method |
CN107085712A (en) * | 2017-04-28 | 2017-08-22 | 山东省农业可持续发展研究所 | A kind of agricultural arid monitoring method based on MODIS data |
CN107219226A (en) * | 2017-07-24 | 2017-09-29 | 中国科学院遥感与数字地球研究所 | Image collecting device and enhancing vegetation index monitoring system |
CN108169161A (en) * | 2017-12-12 | 2018-06-15 | 武汉大学 | A kind of corn planting regional soil humidity appraisal procedure based on modified MODIS indexes |
CN108169161B (en) * | 2017-12-12 | 2019-12-24 | 武汉大学 | Corn planting area soil humidity assessment method based on improved MODIS index |
CN108154193A (en) * | 2018-01-16 | 2018-06-12 | 黄河水利委员会黄河水利科学研究院 | A kind of long-term sequence precipitation data NO emissions reduction method |
CN108154193B (en) * | 2018-01-16 | 2021-10-08 | 黄河水利委员会黄河水利科学研究院 | Long-time sequence precipitation data downscaling method |
CN108896188A (en) * | 2018-06-06 | 2018-11-27 | 北京空间飞行器总体设计部 | A kind of high rail remote sensing satellite infrared channel atmospheric correction method of Shared aperture |
CN109063553A (en) * | 2018-06-22 | 2018-12-21 | 中国矿业大学 | Field-crop growth defect area's remote sensing fast diagnosis method after a kind of land control |
CN108982548A (en) * | 2018-07-20 | 2018-12-11 | 浙江大学 | A kind of soil moisture inversion method based on passive microwave remote sensing data |
CN108982548B (en) * | 2018-07-20 | 2020-01-17 | 浙江大学 | Surface soil moisture retrieval method based on passive microwave remote sensing data |
CN108985959A (en) * | 2018-08-09 | 2018-12-11 | 安徽大学 | Wheat powdery mildew remote sensing monitoring method based on surface temperature inversion technology |
CN108985959B (en) * | 2018-08-09 | 2021-05-28 | 安徽大学 | Wheat powdery mildew remote sensing monitoring method based on surface temperature inversion technology |
CN109186774A (en) * | 2018-08-30 | 2019-01-11 | 清华大学 | Surface temperature information acquisition method, device, computer equipment and storage medium |
CN110163472A (en) * | 2019-04-11 | 2019-08-23 | 中国水利水电科学研究院 | A wide range of extreme drought emergency monitoring and impact evaluation method and system |
CN110111379A (en) * | 2019-05-06 | 2019-08-09 | 新疆农业大学 | A kind of modified surface temperature daily difference NO emissions reduction evaluation method based on MODIS Yu Landsat data |
CN110111379B (en) * | 2019-05-06 | 2023-01-31 | 新疆农业大学 | Improved earth surface temperature daily difference downscaling estimation method based on MODIS and Landsat data |
CN110175932A (en) * | 2019-05-16 | 2019-08-27 | 中国农业科学院农业资源与农业区划研究所 | The Agriculture Drought monitoring method coupled with more ten days precipitation is extracted based on remote sensing information |
CN110263409A (en) * | 2019-06-14 | 2019-09-20 | 中国农业科学院农业资源与农业区划研究所 | A kind of staple crops Low temperature disaster remote sensing monitoring forecast system and method |
CN110243409A (en) * | 2019-06-18 | 2019-09-17 | 中国农业科学院农业资源与农业区划研究所 | A kind of eco-drought monitoring and forecasting system and method based on earth's surface water-heat process |
CN110287944A (en) * | 2019-07-04 | 2019-09-27 | 吉林大学 | The crop pests monitoring method of multi-spectrum remote sensing image based on deep learning |
CN110287944B (en) * | 2019-07-04 | 2022-07-01 | 吉林大学 | Crop pest monitoring method based on multispectral remote sensing image of deep learning |
CN117093813A (en) * | 2023-09-06 | 2023-11-21 | 中国科学院地理科学与资源研究所 | Composite drought index calculation method |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102103077A (en) | MODIS data-based agricultural drought monitoring method | |
Becker-Reshef et al. | A generalized regression-based model for forecasting winter wheat yields in Kansas and Ukraine using MODIS data | |
CN109829234B (en) | A kind of across scale Dynamic High-accuracy crop condition monitoring and yield estimation method based on high-definition remote sensing data and crop modeling | |
Xiao et al. | A continuous measure of gross primary production for the conterminous United States derived from MODIS and AmeriFlux data | |
Xin et al. | Mapping crop key phenological stages in the North China Plain using NOAA time series images | |
CN107782701A (en) | A kind of agricultural arid monitoring method of multi- source Remote Sensing Data data | |
Le-Xiang et al. | Impacts of land use and cover change on land surface temperature in the Zhujiang Delta | |
Meng et al. | Corn yield forecasting in northeast china using remotely sensed spectral indices and crop phenology metrics | |
CN110927120B (en) | Early warning method for vegetation coverage | |
Wang et al. | GIS and RS based assessment of cultivated land quality of Shandong Province | |
Tian et al. | Investigating the urban-induced microclimate effects on winter wheat spring phenology using Sentinel-2 time series | |
Kimura | Effect of the strong wind and land cover in dust source regions on the Asian dust event over Japan from 2000 to 2011 | |
CN113887024A (en) | Surface soil moisture inversion method based on normalized temperature construction and drought index | |
Wang et al. | Research on cropping intensity mapping of the Huai River Basin (China) based on multi-source remote sensing data fusion | |
Shahzad et al. | Spatio-temporal variations in trends of vegetation and drought changes in relation to climate variability from 1982 to 2019 based on remote sensing data from East Asia | |
Sahu et al. | Land use land cover (LU/LC) change analysis using earth observation data sets over Jharsuguda districts of Odisha | |
Liu et al. | Vegetation cover change and its response to climate extremes in the Yellow River Basin | |
Chai et al. | Estimating fractional cover of non-photosynthetic vegetation in a typical grassland area of northern China based on Moderate Resolution Imaging Spectroradiometer (MODIS) image data | |
Darmawan et al. | Seasonal analysis of precipitation, drought and vegetation index in Indonesian paddy field based on remote sensing data | |
Liu et al. | Crop drought area extraction based on remote sensing time series spatial-temporal fusion vegetation index | |
Kanemasu et al. | Estimating grassland biomass using remotely sensed data | |
Matongera et al. | Characterizing bracken fern phenological cycle using time series data derived from Sentinel-2 satellite sensor | |
Käfer et al. | Modeling Pinus elliottii growth with multitemporal Landsat data: a study case in southern Brazil | |
Shah et al. | Evaluation of land surface temperature and normalized difference vegetation index relationship using landsat 8 Satellite images in Mehar Taluka, Dadu | |
Xie et al. | Spatial-temporal variation of satellite-based gross primary production estimation in wheat-maize rotation area during 2000–2015 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C02 | Deemed withdrawal of patent application after publication (patent law 2001) | ||
WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20110622 |