农作物长势综合监测——以印度为例

提出农作物长势综合监测方法,利用卫星遥感得到的NDVI时间序列数据,综合采用实时监测、过程监测和时间序列聚类监测方法,明确不同方法适用的监测尺度及监测目的,对不同范围农作物长势进行监测。改进了Crop Watch全球农情遥感速报系统运行化作物长势监测方法,克服了原有作物长势监测中实时监测方法无法反映相同区域苗情在整个生长过程中的连续变化情况的缺点。实现对相同区域作物长势连续变化的定量描述,可对作物长势进行更准确的判断。利用官方发布的作物单产变幅数据,对单产变幅较大的12个作物主产省区作物长势监测结果的准确性进行判断,结果表明:6个邦的实时监测和聚类监测方法所得结果一致,都符合作物单产变化的实际状况;4个邦的聚类监测方法所得结果对作物长势监测更为准确,更符合该区域作物单产的实际变化;1个邦实时监测结果对作物长势监测比聚类监测方法更为准确;只有1个邦采用两种方法对作物长势的监测存在误差,聚类监测方法在对农作物生长过程的连续监测及空间分布的定量化表述方面,比实时监测更为准确。3种方法可以综合使用,实现业务化运行的农作物长势监测。     

Development and validation of spectrometeorological yield models of mustard oil seed crop in Rajasthan state

Attempt has been made to develop spectro meteorological yield models using normalized difference vegetation index (NDVI) derived from NOAA AVHRR data over the crop growth period and monthly rainfall data for predicting yield of mustard crop. The AVHRR data spanning seven crop growing seasons, the rain gauze station-level rainfall data and crop yield data determined from crop cutting experiments (CCE) conducted by state Directorate of Economics and Statistics (DES) are the basic input data. A methodology has been developed to normalize the multi-temporal NDVIs for the minimisation of atmospheric effects, which is found to reduce the noise in NDVI due to varying atmospheric conditions from season to season and improve the predictability of statistical multiple linear regression yield models developed for nine geographically large districts of Rajasthan state. The spectro meteorological yield models had been validated by comparing the predicted district level yields with those estimated from the crop cutting experiments.     

Relation of wheat yield with parameters derived from a spectral growth profile

An attempt has been made to generate crop growth profiles using multi-date NOAA AVHRR data of wheat-growing season of 1987–88 for the districts of Punjab and Haryana states of India. A profile model proposed by Badhwar was fitted to the multi-date Normalised Difference Vegetation Index (NDVI) values obtained from geographically referenced samples in each district. A novel approach of deriving a set of physiologically meaningful profile parameters has been outlined and the relation of these parameters with district wheat yields has been studied in order to examine the potential of growth profiles for crop-yield modelling. The parameter ‘area under the profile’ is found to be the best estimator of yield. However, with such a parameter time available for prediction gets reduced. Combination of different profile parameters shows improvement in correlation but lacks the consistency for individual state data.     

Interpretation of crop growth patterns by means of NDVI‐time series in Zambia

A method to correlate crop production in Zambia to the yearly evolution of the Normalized Difference Vegetation Index (NDVI) is proposed. The method consists of the analysis of remote sensing data together with meteorological data and simulated crop production to obtain indicators of crop production. The accuracy of these indicators is assessed with statistical data.

The main objective was to assess whether the NDVI‐time series extracted from NOAA‐AVHRR‐images , having a pixel resolution of 73 km may give reliable information on crop production in Zambia where agricultural areas cover just 1% of the land area.

The mean NDVI‐value of several pixels, e.g. for one province or other administrative units, relates to the dominant type of vegetation in the area under consideration.

It is shown that the 7.3 km NDVI‐data give reliable indications on crop production in Zambia, when small areas (200–450 km² large ) are considered where agricultural land use is intensive. This implies that preliminary analysis is required to localize the agricultural areas. This has been done by means of high resolution satellite images i.e. LANDSAT‐MultiSpectral Scanner.

Consequently, the NDVI‐time series of the ‘agricultural ‘ pixels are used to calculate crop growth indicators which can be applied to assess the crop production.     

农作物长势综合遥感监测方法

作物收获之前进行大范围作物生长状况评价 ,可以尽早的获得有关作物产量信息。介绍了中国农情遥感监测系统的综合作物长势监测方法。以遥感数据标准化处理、云标识、云污染去除和非耕地去除为基础 ,生成质量一致的遥感数据产品集 ,提取区域作物生长过程。作物长势监测分为实时作物长势监测和作物生长趋势分析。实时的作物长势监测可以定性和定量地在空间上分析作物生长状况 ,分级显示作物生长状况 ,分区域统计水田和旱地中不同长势占的比重。作物生长趋势分析可以进行年际间的生长过程对比 ,从时间轴上反映作物持续生长的差异性 ,统计全国、主产区、省和区划单元 4个尺度的耕地、水田、旱地作物生长过程曲线年际间差异 ,从而为早期的产量预测提供信息。通过处理流程的系统化 ,建设了运行化的作物长势遥感监测分析系统 ,为用户构建了综合的作物实时生长状况 ,苗情的生长趋势分析环境。同时可以依据野外地面实测信息对遥感监测结果进行标定和检验。 1998年以来 ,系统在满足日常运行的前提下 ,技术方法逐渐改进和完善 ,监测结果的精度和可靠性不断得到提高。     

中国农情遥感速报系统

介绍了中国农情遥感速报系统的建设情况 ,系统内容包括农作物长势监测、农作物种植面积监测、农作物单产预测与粮食产量估算、作物时空结构监测和粮食供需平衡预警等。简要介绍了 1998年以来中国农情遥感速报系统在监测内容与监测范围、监测频率、技术发展以及质量控制与过程检验体系建立等方面的进展 ,并就中国农情遥感速报系统的发展方向提出了展望。     

陆面温度反演算法——劈窗算法的敏感度分析

首先给出ＣＯ２ 倍增下遥感光合作物产量的概念模型,之后分析未受ＣＯ２ 倍增的遥感光合作物产量估测模型;在考虑ＣＯ２ 倍增对作物产量的影响后,对影响干物质累积的作物光合速率的模型进行修正,进而修正遥感光合作物产量估测模型。建立ＣＯ２ 倍增下作物产量响应模型,求取各参数,并在ＣＯ２ 倍增下对我国华北地区冬小麦产量响应进行填图,表明模型的估测结果有良好的可比性。     

Analysis of sorghum yield dependence on crop growth characteristics using noaa-avhrr data

Growth profiles of 1987-88 rabi sorghum crop cultivated in spatially extensive sites in six tehsils of Solapur and Ahmadnagar districts in Maharashtra have been generated using multidate NOAA AVHRR data based on crop growth equation suggested by Badhwar (1980). The sensitive parameters for sorghum yield modelling have been identified. The correlation of final grain yield with growth parameters shows that yield relationship is stronger when logarithmic senescence rate and timeintegrated logarithmic senescence rate are considered as the parameters instead of its value on any day during 30 days senescence period after attaining maximum vegetative cover.     

Extracting Regional Pattern of Wheat Sowing Dates Using Multispectral and High Temporal Observations from Indian Geostationary Satellite

Monitoring of Agricultural crops using remote sensing data is an emerging tool in recent years. Spatial determination of sowing date is an important input of any crop model. Geostationary satellite has the capability to provide data at high temporal interval to monitor vegetation throughout the entire growth period. A study was conducted to estimate the sowing date of wheat crop in major wheat growing states viz. Punjab, Haryana, Uttar Pradesh (UP), Madhya Pradesh (MP), Rajasthan and Bihar. Data acquired by Charged Couple Detector (CCD) onboard Indian geostationary satellite INSAT 3A have continental (Asia) coverage at 1 km?×?1 km spatial resolution in optical spectral bands with high temporal frequency. Daily operational Normalized Difference Vegetation Index (NDVI) product from INSAT 3A CCD available through Meteorological and Oceanographic Satellite Data Archival Centre (MOSDAC) was used to estimate sowing date of wheat crop in selected six states. Daily NDVI data acquired from September 1, 2010 to December 31, 2010 were used in this study. A composite of 7 days was prepared for further analysis of temporal profile of NDVI. Spatial wheat crop map derived from AWiFS (56 m) were re-sampled at INSAT 3A CCD parent resolution and applied over each 7 day composite. The characteristic temporal profiles of 7 day NDVI composite was used to determine sowing date. NDVI profile showed decreasing trend during maturity of kharif crop, minimum value after harvest and increasing trend after emergence of wheat crop. A mathematical model was made to capture the persistent positive slope of NDVI profile after an inflection point. The change in behavior of NDVI profile was detected on the basis of change in NDVI threshold of 0.3 and sowing date was estimated for wheat crop in six states. Seven days has been deducted after it reached to threshold value with persistent positive slope to get sowing date. The clear distinction between early sowing and late sowing regions was observed in study area. Variation of sowing date was observed ranging from November 1 to December 20. The estimated sowing date was validated with the reported sowing date for the known wheat crop regions. The RMSD of 3.2 (n?=?45) has been observed for wheat sowing date. This methodology can also be applied over different crops with the availability of crop maps.     

基于MODIS数据的小麦生物量估算模型研究

地表生物量对农作物估产、植被长势评估具有很重要的意义。随着遥感技术的发展与应用,遥感为生物量估算提供了一种新的手段。本文以唐山市为例,利用小麦种植区的MODIS遥感影像数据和同期野外调查获得的16组32个生物量数据,对比分析了归一化植被指数(NDVI)、增强型植被指数(EVI)与小麦生物量多个回归方程的相关系数,进而建立了NDVI、EVI与小麦生物量的线性回归模型。结果显示,使用MODIS数据的植被指数能够很好地对研究区地上生物量进行估算,其中使用EVI的三次函数模型拟合精度最高,并且对每组数据进行平均处理会使模型精度进一步提高。     

区域作物生长过程的遥感提取方法

提出利用时序NDVI数据提取作物生长过程方法。遥感数据在采集过程中受云、大气因子的影响 ,以及混合像元问题 ,造成时序植被指数值变得没有规律 ,对比性不强。采用基于最小二次方拟合的谐函数分析方法 ,依据作物轮作规律和生长周期性特征 ,用主要频率的正弦、余弦谐函数重建时序图像 ,去除了影像中云污染的影响。以中国的旱地为例 ,考虑到像元内旱地对NDVI值的贡献率 ,计算区域内旱地像元加权平均值来反映其作物生长过程。同时与区域所有像元的平均值、旱地平均值等统计方法的结果进行对比分析 ,表明区域内旱地的加权平均值能够削弱旱地比例和地域间的差异 ,突出耕地上作物的生长过程特征。通过与地面实测数据分析 ,平滑前后的作物生长过程与叶面积指数相关性增加 5 %— 11% ,采用区域加权平均的方法得到的作物生长过程 ,比旱地平均和NDVI平均的结果与叶面积指数的相关性增加 14 %— 17%。     

用NOAA气象卫星AVHRR的定量资料计算冬小麦种植面积的两种方法

本文提出了用NOAA气象卫星AVHRR的定量资料计算冬小麦种植面积的方法,即绿度指数法和单通道法。绿度指数法是根据冬小麦在特定生育期内绿度值基本上保持为常数的特点,用几块巳知样地的种植面积,推算整体种植面积;单通道法是只使用AVHRR第二通道的反照率来计算冬小麦种植面积。这两种方法的优点在于不需考虑大气削弱的订正处理问题,从而使计算得以简化。     

Wheat crop production estimation using satellite data

This paper reports acreage, yield and production forecasting of wheat crop using remote sensing and agrometeorological data for the 1998–99 rabi season. Wheat crop identification and discrimination using Indian Remote Sensing (IRS) ID LISS III satellite data was carried out by supervised maximum likelihood classification. Three types of wheat crop viz. wheat-1 (high vigour-normal sown), wheat-2 (moderate vigour-late sown) and wheat-3 (low vigour-very late sown) have been identified and discriminated from each other. Before final classification of satellite data spectral separability between classes were evaluated. For yield prediction of wheat crop spectral vegetation indices (RVI and NDVI), agrometeorological parameters (ET_max and TD) and historical crop yield (actual yield) trend analysis based linear and multiple linear regression models were developed. The estimated wheat crop area was 75928.0 ha. for the year 1998–99, which sowed ?2.59% underestimation with land record commissioners estimates. The yield prediction through vegetation index based and vegetation index with agrometeorological indices based models were 1753 kg/ha and 1754 kg/ha, respectively and have shown relative deviation of 0.17% and 0.22%, the production estimates from above models when compared with observed production show relative deviation of ?2.4% and ?2.3% underestimations, respectively.     

多时相MODIS影像水田信息提取研究

水稻种植及其分布信息是土地覆被变化、作物估产、甲烷排放、粮食安全和水资源管理分析的重要数据源。基于遥感的水田利用监测中,通常采用时序NDVI植被指数法和影像分类法分别进行AVHRR和TM影像的水田信息获取。针对8天合成MODIS陆地表面反射比数据的特点和水稻生长特征,选取水稻种植前的休耕期、秧苗移植期、秧苗生长期和成熟期等多时相MODIS地表反射率影像数据,通过归一化植被指数、增强植被指数及利用对土壤湿度和植被水分含量较敏感的短波红外波段计算得到的陆表水指数进行水田信息获取。将提取结果与基于ETM＋影像的国土资源调查水田数据,通过网格化计算处理并进行对比分析,结果表明,利用MODIS影像的8天合成地表反射率数据,进行区域甚至全国的水田利用监测是可行的。     

Monitoring winter wheat growth in North China by combining a crop model and remote sensing data

Both of crop growth simulation models and remote sensing method have a high potential in crop growth monitoring and yield prediction. However, crop models have limitations in regional application and remote sensing in describing the growth process. Therefore, many researchers try to combine those two approaches for estimating the regional crop yields. In this paper, the WOFOST model was adjusted and regionalized for winter wheat in North China and coupled through the LAI to the SAIL–PROSPECT model in order to simulate soil adjusted vegetation index (SAVI). Using the optimization software (FSEOPT), the crop model was then re-initialized by minimizing the differences between simulated and synthesized SAVI from remote sensing data to monitor winter wheat growth at the potential production level. Initial conditions, which strongly impact phenological development and growth, and which are hardly known at the regional scale (such as emergence date or biomass at turn-green stage), were chosen to be re-initialized. It was shown that re-initializing emergence date by using remote sensing data brought simulated anthesis and maturity date closer to measured values than without remote sensing data. Also the re-initialization of regional biomass weight at turn-green stage led that the spatial distribution of simulated weight of storage organ was more consistent to official yields. This approach has some potential to aid in scaling local simulation of crop phenological development and growth to the regional scale but requires further validation.     

Comparison of eMODIS and MOD/MYD13A2 NDVI products during 2012–2014 spring green-up periods in Alaska and northwest Canada

Accurate monitoring of vegetation dynamics is required to understand the inter-annual variability and long term trends in terrestrial carbon exchange in tundra and boreal ecoregions. In western North America, two Normalized Vegetation Index (NDVI) products based on spectral reflectance data from the Moderate Resolution Imaging Spectroradiometer (MODIS) are available. The MOD/MYD13A2 NDVI product is available as a 16-day composite product in a sinusoidal projection as global hdf tiles. The eMODIS Alaska NDVI product is available as a 7-day composite geotif product in a regional equal area conic projection covering Alaska and the entire Yukon River Basin. These two NDVI products were compared for the 2012–2014 late May–late June spring green-up periods in Alaska and the Yukon Territory. Relative to the MOD/MYD13A2 NDVI product, it is likely that the eMODIS NDVI product contained more cloud-contaminated NDVI values. For example, the MOD/MYD13A2 product flagged substantially fewer pixels as “good quality” in each 16-day composite period compared to the corresponding MODIS Alaska NDVI product from a 7-day composite period. During the spring green-up period, when field-based NDVI increases, the eMODIS NDVI product averaged 43 percent of pixels that declined by at least 0.05 NDVI between 2 composite periods, consistent with cloud-contamination problems, while the MOD/MYD13A2 NDVI averaged only 6 percent of pixels. Based on a cloudy Landsat-8 scene, the eMODIS compositing process selected 23 percent pixels, while the MOD/MYD13A2 compositing process selected less than 0.003 percent pixels. Based on the results, it appears that the MOD/MYD13A2 NDVI product is superior for scientific applications based on NDVI phenology in the tundra and boreal regions of northwestern North America.     

Wheat Monitoring by Using Kernel Based Possibilistic <Emphasis Type="Italic">c</Emphasis>-Means Classifier: Bi-sensor Temporal Multi-spectral Data

The existence of mixed pixels in the satellite images has always been an area of concern. Adding to the challenge is an occurrence of non-linearity between the classes, which is generally overlooked. The study makes an attempt to solve the two frequently occurring problems by kernel based fuzzy approach. This research work deals with Possibilistic c-Means (PCM) classifier with local, global, spectral angle and hyper tangent kernels for wheat crop (Triticum aestivum) identification in Haridwar, Uttarakhand, India. The multi-temporal vegetation index data of Formosat-2 have been used which covers the whole phenology of wheat crop. The additional sensor Landsat-8 OLI imagery has been filled the crucial gap of Formosat-2 temporal datasets. Nine types of proposed kernels based PCM classifier have been applied on three temporal datasets (four, five and six date combinations) to classify two classes early sown and late sown wheat crop. These test results have been concluded that at optimized weighted constant KMOD and polynomial kernel was found effective to separate wheat crop. The five and six date combination were sufficient to discriminate early sown and late sown wheat crop.     

Assessment of crop growth parameters of wheat under stress condition through ground based spectral data

A field experiment was conducted on wheat crop during rabi seasons of 1995–96, 1996–97 and 1997–98 to study the spectral response of wheat crop (between 490 to 1080 nm) under water and nutrient stress condition. An indigenously developed ground truth radiometer having narrow band in visible and near infrared region (490 – 1080 nm) was used. Vegetation indices derived using different band combinations and related to crop growth parameters. The near infrared spectral region of 710 – 1025 nm was found most important for monitoring stress condition. Relationship has been developed between crop growth parameters and vegetation indices. Leaf Area Index (LAI) and chlorophyll could be predicted by knowing different reflectance ratios at milking stage of crop with R² value of 0.78 and 0.89, respectively. Dry biomass (DBM), Plant Water Content (PWC) and grain yield are also significantly related with reflectance ratios at flowering stage of crop with R² value of 0.90, 0.98 and 0.74, respectively.     

Relationship of wheat grain yield with spectral indices

An experiment was conducted during 1996–97 and 1997–98 to study spectral indices and their relationships with grain yield of wheat. Variations of ratio vegetation index (RVI), normalized differences vegetation index (NDVI). difference vegetation index (DVI), transformed vegetation index (TVI), perpendicular vegetation index (PVI) and greenness vegetation index (GVI) have been studied at anthesis stage under different moisture and nitrogen levels. Spectral indices were correlated with crop parameters and it was found that GVI was the best index for yield estimation (r = 0.91 ).     

Crop Evapotranspiration Estimation for Wheat (<Emphasis Type="Italic">Triticum aestivum</Emphasis> L.) Using Remote Sensing Data in Semi-Arid Region of Maharashtra

Penman–Monteith method adapted to satellite data was used for the estimation of wheat crop evapotranspiration during the entire growth period using satellite data together with ground meteorological measurements. The IRS-1D/IRS-P6 LISS-III sensor data at 23.5 m spatial resolution for path 096 and row 059 covering the study area were used to derive, albedo, normalized difference vegetation index, leaf area index and crop height and then to estimate wheat crop evapotranspiration referred to as actual evapotranspiration (ETact). The ETact varied from 0.86 to 3.41 mm/day during the crop growth period. These values are on an average 16.40 % lower than wheat crop potential evapotranspiration (ETc) estimated as product of reference crop evapotranspiration estimated by Penman–Monteith method and lysimetric crop coefficient (Kc). The deviation of ETact from ETc is significant, when both the values were compared with t test for paired two sample means. Though the observations on ETact were taken from well maintained unstressed experimental plot of 120 × 120 m size, there was significant deviation. This deviation could be attributed to, the satellite images representing the actual crop evapotranspiration as function crop canopy biophysical parameters, condition of the crop stand, climatic and soil conditions and the microclimate variation over area of one hectare. However, Penman–Monteith method represents a flat rate of specific growth stage of the crop.