Downscaling MODIS LST in the East African mountains using elevation gradient and land-cover information

Land-surface temperature (LST) is strongly affected by altitude and surface albedo. In mountain regions where steep slopes and heterogeneous land cover are predominant, LST can vary significantly within short distances. Although remote sensing currently provides opportunities for monitoring LST in inaccessible regions, the coarse resolution of some sensors may result in large uncertainties at sub-pixel scales. This study aimed to develop a simple methodology for downscaling 1 km Moderate Resolution Spectroradiometer (MODIS) LST pixels, by accounting for sub-pixel LST variation associated with altitude and land-cover spatial changes. The approach was tested in Mount Kilimanjaro, Tanzania, where changes in altitude and vegetation can take place over short distances. Daytime and night-time MODIS LST estimates were considered separately. A digital elevation model (DEM) and normalized difference vegetation index (NDVI), both at 250 m spatial resolution, were used to assess altitude and land-cover changes, respectively. Simple linear regressions and multivariate regressions were used to quantify the relationship between LST and the independent variables, altitude and NDVI. The results show that, in Kilimanjaro, altitude variation within the area covered by a 1 km MODIS LST pixel can be up to ±300 m. These altitude changes can cause sub-pixel variation of up to ±2.13°C for night-time and ±2.88°C for daytime LST. NDVI variation within 1 km pixels ranged between –0.2 and 0.2. For night-time measurements, altitude explained up to 97% of LST variation, while daytime LST was strongly affected by land cover. Using multivariate regressions, the combination of altitude and NDVI explained up to 94% of daytime LST variation in Kilimanjaro. Finally, the downscaling approach proposed in this study allowed an improved representation of the influence of landscape features on local-scale LST patterns.     

Super-resolution land cover mapping with indicator geostatistics

Many satellite images have a coarser spatial resolution than the extent of land cover patterns on the ground, leading to mixed pixels whose composite spectral response consists of responses from multiple land cover classes. Spectral unmixing procedures only determine the fractions of such classes within a coarse pixel without locating them in space. Super-resolution or sub-pixel mapping aims at providing a fine resolution map of class labels, one that displays realistic spatial structure (without artifact discontinuities) and reproduces the coarse resolution fractions. In this paper, existing approaches for super-resolution mapping are placed within an inverse problem framework, and a geostatistical method is proposed for generating alternative synthetic land cover maps at the fine (target) spatial resolution; these super-resolution realizations are consistent with all the information available.More precisely, indicator coKriging is used to approximate the probability that a pixel at the fine spatial resolution belongs to a particular class, given the coarse resolution fractions and (if available) a sparse set of class labels at some informed fine pixels. Such Kriging-derived probabilities are used in sequential indicator simulation to generate synthetic maps of class labels at the fine resolution pixels. This non-iterative and fast simulation procedure yields alternative super-resolution land cover maps that reproduce: (i) the observed coarse fractions, (ii) the fine resolution class labels that might be available, and (iii) the prior structural information encapsulated in a set of indicator variogram models at the fine resolution. A case study is provided to illustrate the proposed methodology using Landsat TM data from SE China.     

基于局部连续性与全局相似性的光谱保持型亚像元映射

遥感图像的像元级分类精度受混合像元的影响. 亚像元映射以像元分解获得的丰度值为基础,在地物分布规律的约束下,细化估计各类地物的亚像元级分布模式. 本文同时考虑了地物分布的空间与光谱信息,提出了一种基于局部连续性与全局相似性的光谱保持型亚像元映射算法. 针对地物的空间分布特性,提出了利用类内离散度对局部连续性进行建模,并通过相似分布像元表示误差引入全局相似性约束项. 针对地物的光谱特性,采用最小化光谱误差约束了亚像元映射过程中的光谱无失真性. 模拟数据与真实数据上的实验结果表明,本文算法比其他同类算法具有更高的估计精度,且更适合于实际应用.     

Soft-then-hard super-resolution mapping based on a spatial attraction model with multiscale sub-pixel shifted images

Super-resolution mapping (SRM) is a technique for exploring spatial distribution information of the land-cover classes at finer spatial resolution. The soft-then-hard super-resolution mapping (STHSRM) algorithm is a type of SRM algorithm that first estimates the soft class values for sub-pixels at the target fine spatial resolution and then predicts the hard class labels for sub-pixels. The sub-pixel shifted images from the same area can be incorporated to improve the accuracy of STHSRM algorithm. In this article, multiscale sub-pixel shifted images (MSSI) based on the fine-scale model and the coarse-scale model are utilized to increase the accuracy of STHSRM. First, class fraction images are derived from multiple sub-pixel shifted coarse spatial resolution images by soft classification. Then using the sub-pixel/sub-pixel spatial attraction model as fine-scale and the sub-pixel/pixel spatial attraction model as coarse scale, all MSSI can be derived from fraction images. The MSSI for each class are then integrated to obtain the desired fine spatial resolution images. Finally, the integrated fine spatial resolution images are used to allocate classes for sub-pixel. Experiments on two synthetic remote sensing images and a real hyperspectral remote sensing imagery show that the proposed method produces higher mapping accuracy result.     

Spatially adaptive smoothing parameter selection for Markov random field based sub-pixel mapping of remotely sensed images

Sub-pixel mapping is a process to provide the spatial distributions of land cover classes with finer spatial resolution than the size of a remotely sensed image pixel. Traditional Markov random field-based sub-pixel mapping (MRF_SPM) adopts a fixed smoothing parameter estimated based on the entire image to balance the spatial and spectral energies. However, the spectra of the remotely sensed pixels are always spatially variable. Adopting a fixed smoothing parameter disregards the local properties provided by each pixel spectrum, and may probably lead to insufficient smoothing in the homogeneous region and over-smoothing between class boundaries simultaneously. This article proposes a spatially adaptive parameter selection method for the MRF_SPM model to overcome the limitation of the fixed parameter. As pixel class proportions are indicators of the type and proportion of land cover classes within each coarse pixel, in the proposed method, fraction images providing pixel class proportions as local properties of each pixel spectrum are employed to constrain the smoothing parameter. Consequently, the smoothing parameter is spatially adaptive to each pixel spectrum of the remotely sensed image. Synthetic images and IKONOS multi-spectral images were employed. Results showed that compared with the hard classification method and the non-spatially adaptive MRF_SPM adopting a fixed smoothing parameter, the spatially adaptive MRF_SPM with the smoothing parameter constrained to each pixel spectrum yielded sub-pixel maps not only with higher accuracy but also with shapes and boundaries visually reconstructed more closely to the reference map.     

Evaluating MODIS soil fractional cover for arid regions,using albedo from high-spatial resolution satellite imagery

Broad-scale high-temporal frequency satellite imagery is increasingly used for environmental monitoring. While the normalized difference vegetation index (NDVI) is the most commonly used index to track changes in vegetation cover, newer spectral mixture approaches aim to quantify sub-pixel fractions of photosynthesizing vegetation, non-photosynthesizing vegetation, and exposed soil. Validation of the unmixing products is essential to enable confident use of the products for management and decision-making. The most frequently used validation method is by field data collection, but this is very time consuming and costly, in particular in remote regions where access is difficult.

This study developed and demonstrates an alternative method for quantifying land-cover fractions using high-spatial resolution satellite imagery. The research aimed to evaluate the bare soil fraction in a sub-pixel product, MODIS Fract-G, for the natural arid landscapes of the far west of South Australia. Twenty-two sample regions, of 3400 sampling points each, were investigated across several arid land types in the study area. Albedo thresholds were carefully determined in Advanced Land Observing Satellite Panchromatic Remote-sensing Instrument Stereo Mapping (ALOS PRISM) images (2.5 m spatial resolution), which separated predominantly bare soil from predominantly vegetated or covered soil, and created classified images. Correlation analysis was carried out between MODIS Fract-G bare soil fractional cover and ALOS PRISM bare soil proportions for the same areas. Results showed much lower correlations than expected, though limited agreement was found in some specific areas. It is posited that the Moderate Resolution Imaging Spectroradiometer (MODIS) fractional cover product, which is based on unmixing using the NDVI and a cellulose absorption index (CAI) proxy, may be generally unable to separate soil from vegetation in situations where both indices are low. In addition, separation is hampered by the lack of ‘pure pixels’ in this heterogeneous landscape. This suggests that the MODIS fractional cover product, at least in its present form, is unsuited to monitor sparsely vegetated arid landscapes.     

Vectorial boundary-based sub-pixel mapping method for remote-sensing imagery

This article presents a vectorial boundary-based sub-pixel mapping (VBSPM) method to obtain the land-cover distribution with finer spatial resolution in mixed pixels. With inheritance from the geometric SPM (GSPM), VBSPM first geometrically partitions a mixed pixel using polygons, and then utilizes a vectorial boundary extraction model (VBEM), rather than the rasterization method in GSPM, to determine the location and length of each edge in the polygon, while these edges are located at the boundary of and within the interior of the mixed pixel. Furthermore, VBSPM uses a decay function to manage the mixed pixels along the image boundary region due to the missing parts of their neighbours. Finally, a ray-crossing algorithm is employed to determine the land-cover class of each sub-pixel in terms of vectorial boundaries. The experiments with artificial and remotely sensed images have demonstrated that VBSPM can reduce the inconsistency between the boundaries of different land-cover classes, approximately calculating errors with an odd zoom factor, and achieve more accurate sub-pixel mapping results than the hard classification methods and GSPM.     

Land cover change detection at coarse spatial scales based on iterative estimation and previous state information

This study focuses on the use of coarse spatial resolution (CR, pixel size about 1 km²) remote sensing data for land cover change detection and qualification. Assuming the linear mixing model for CR pixels, the problem is that both the multitemporal class features and the pixel composition in terms of classes are unknown. The proposed algorithm is then based on the iterative alternate estimation of each unknown variable. At each iteration, the class features are estimated, thanks to the knowledge of the composition of some pixels, and then the pixel composition is re-estimated knowing the class features. The subset of known composition pixels is the subset of pixels where no change has occurred, i.e. the previous land cover map is still valid. It is derived automatically by removing at each iteration the pixels where the new composition estimation disagrees with the former one. Finally, for the final estimation of the pixel composition, a Markovian chain model is used to guide the solution, i.e. the previous land cover map is used as a ‘reminder’ or ‘memory’ term.This approach has been first validated using simulated data with different spatial resolution ratios. Then, the detection of forest change with SPOT/VGT-S10 has been considered as an actual application case. Finally, the method has been applied to change detection on the Val de Saône watershed between the 1980s and 2000. The results obtained from three coarse resolution series, NOAA/AVHRR, SPOT/VGT-S10 and SPOT/VGT-P, have been compared.     

Variation in NDVI values with change in spatial resolution for semi-arid savanna vegetation: a case study in northwestern South Africa

Natural vegetation and crop-greening patterns in semi-arid savannas are commonly monitored using normalized difference vegetation index (NDVI) values from low spatial resolution sensors such as the Advanced Very High Resolution Radiometer (AVHRR) (1 km, 4 km) and Moderate Resolution Imaging Spectroradiometer (MODIS) (250 m, 500 m). However, because semi-arid savannas characteristically have scattered tree cover, the NDVI values at low spatial resolution suffer from the effect of aggregation of near-infrared and red energy from adjacent vegetated and non-vegetated cover types. This effect is seldom taken into consideration or quantified in NDVI analyses of the vegetation of semi-arid lands. This study examined the effect of pixel size on NDVI values of land-cover features for a semi-arid area, using the 1000 m, 250 m and 10 m pixel sizes. A rainy season Système Pour l'Observation de la Terre 5 (SPOT 5) High Resolution Geometric (HRG) image at 10 m spatial resolution was utilized. Following radiometric and geometric preprocessing, the 10 m pixel size of the image was aggregated to 250 m and 1000 m to simulate imagery at these pixel sizes, and then NDVI images at the spatial resolution scales of 10 m (NDVI_{10 m}), 250 m (NDVI_{250 m}), and 1000 m (NDVI_{1000 m}) derived from the respective images. The simulation of the NDVI_{250 m} image was validated against a concurrent 16 day MODIS NDVI composite (MOD13Q1) image, and the accuracy derived from the validation was generalized to the NDVI_{1000 m} image. With change from low to high spatial resolution, extreme magnitude NDVI values shifted towards the centre (mode) of the resulting approximately Gaussian NDVI distributions. There was a statistically significant difference in NDVI values at the three pixel sizes. Low spatial magnitude vegetation sites (woodland, cropland) had reductions of up to 28% in NDVI value between the NDVI_{10 m} and NDVI_{1000 m} scales. The results indicate that vegetation monitoring using low spatial resolution imagery in semi-arid savannas may only be indicative and needs to be supplemented by higher spatial resolution imagery.     

Super-resolution land-cover mapping using multiple sub-pixel shifted remotely sensed images

Super-resolution land-cover mapping is a promising technology for prediction of the spatial distribution of each land-cover class at the sub-pixel scale. This distribution is often determined based on the principle of spatial dependence and from land-cover fraction images derived with soft classification technology. However, the resulting super-resolution land-cover maps often have uncertainty as no information about sub-pixel land-cover patterns within the low-resolution pixels is used in the model. Accuracy can be improved by incorporating supplemental datasets to provide more land-cover information at the sub-pixel scale; but the effectiveness of this is limited by the availability and quality of these additional datasets. In this paper, a novel super-resolution land-cover mapping technology is proposed, which uses multiple sub-pixel shifted remotely sensed images taken by observation satellites. These satellites take images over the same area once every several days, but the images are not identical because of slight orbit translations. Low-resolution pixels in these remotely sensed images therefore contain different land-cover fractions that can provide useful information for super-resolution land-cover mapping. We have constructed a Hopfield Neural Network (HNN) model to solve it. Maximum spatial dependence is the goal of the proposed model, and the fraction maps of all images are constraints added to the energy function of HNN. The model was applied to synthetic artificial images as well as to a real degraded QuickBird image. The output maps derived from different numbers of images at different zoom factors were compared visually and quantitatively to the super-resolution map generated from a single image. The resulting land-cover maps with multiple remotely sensed images were more accurate than was the single image map. The use of multiple remotely sensed images is therefore a promising method for decreasing the uncertainty of super-resolution land-cover mapping. Moreover, remotely sensed images with similar spatial resolution from different satellite platforms can be used together, allowing a fusion of information obtained from remotely sensed imagery.     

Use of MODIS NDVI data to improve forest-area estimation

The objective of this article is to develop and test a methodology capable of using medium spatial resolution satellite imagery to improve forest-area statistics derived from ground sampling. The methodology builds on the evidence that multitemporal Normalized Difference Vegetation Index (NDVI) images bring significant information on the spatial distribution of forest surfaces. Consequently, Moderate Resolution Imaging Spectroradiometer (MODIS) NDVI images are potentially useful to improve forest-area assessment based on ground data. This expectation is verified in Tuscany (central Italy) using forest-area references extracted from the Coordination of Information on Environment (CORINE) land-cover map. The accuracy of forest-area statistics obtained at province level by different reference samplings is first assessed. Next, locally calibrated regression analyses are applied to multitemporal MODIS NDVI images in order to obtain per-pixel forest-area estimates. Two statistical methods (the direct expansion and the regression estimator) are finally used to combine these estimates with the ground data and produce corrected per-province statistics. The experimental results confirm that MODIS NDVI data contain relevant information on forest distribution, which can be efficiently extended over the land surface by locally calibrated regressions. The obtained estimates can be combined with the ground data for enhancing forest-area assessment at province level. To this aim, the regression estimator gives the best performance for all sampling densities of the reference data.     

Sub-pixel mapping of remotely sensed imagery with hybrid intra- and inter-pixel dependence

Sub-pixel mapping of remotely sensed imagery is often performed by assuming that land cover is spatially dependent both within and between image pixels. Intra- and inter-pixel dependencies are two widely used approaches to represent different land-cover spatial dependencies at present. However, merely using intra- or inter-pixel dependence alone often fails to fully describe land-cover spatial dependence, making current sub-pixel mapping models defective. A more reasonable object for sub-pixel mapping is maximizing both intra- and inter-pixel dependencies simultaneously instead of using only one of them. In this article, the differences between intra- and inter-pixel dependencies are discussed theoretically, and a novel sub-pixel mapping model aiming to maximize hybrid intra- and inter-pixel dependence is proposed. In the proposed model, spatial dependence is formulated as a weighted sum of intra-pixel dependence and inter-pixel dependence to satisfy both intra- and inter-pixel dependencies. By application to artificial and synthetic images, the proposed model was evaluated both visually and quantitatively by comparing with three representative sub-pixel mapping algorithms: the pixel swapping algorithm, the sub-pixel/pixel attraction algorithm, and the pixel swapping initialized with sub-pixel/pixel attraction algorithm. The results showed increased accuracy of the proposed algorithm when compared with these traditional sub-pixel mapping algorithms.     

Hierarchical Land Cover Classification based on MODIS NDVI Time-series

In this paper,we mainly used MODIS NDVI time-series dataset at 16-days temporal resolution and 250-meters spatial resolution to analyze land cover mapping of northeastern China.We used two different filter methods to fit NDVI time-series dataset,and compared their average classes’ separability based on Jeffries-Matusita distance index.In addition,we made use of hierarchical classification method to complete classification,combined with short-wave infrared spectral reflectance data and DEM.We conformed to the principle that separate area hierarchically into several parts first and then classify each part further,and use a single characteristic band first and then multiple feature bands.In the process of classification,we adopted threshold value method,support vector machine,artificial net neural and C5.0 decision tree classification to distinguish each land-cover type hierarchically.Finally,we evaluated the accuracy of the final classification of study area using known land-cover classification data and high-resolution remote sensing imagery,overall accuracy is 84.61%,Kappa coefficient is 0.8262.     

遥感影像亚像元制图方法研究进展综述

遥感影像混合像元的普遍存在给遥感影像解译造成困扰。有效处理混合像元问题,细化分类结果,获得更为精细的地物细节信息就需要进行亚像元绘图。目前亚像元制图方法主要包括3个步骤:① 混合像元分解;② 提取软信息;③ 亚像元制图。总结归纳了近年来遥感影像亚像元绘图领域的研究进展和成果,详细阐述了亚像元制图的步骤及涉及的研究方法。依据辅助信息的类型将亚像元绘图方法大致划分为:基于空间相关性、基于空间结构信息、基于神经网络、基于像元交换途径的4类亚像元分类方法,并分别对各种方法的优缺点进行了分析对比。最后,评述了亚像元制图的发展趋势。     

Retrieval of sub-pixel snow cover information in the Himalayan region using medium and coarse resolution remote sensing data

Numerous constrained and unconstrained algorithms have been used to retrieve sub-pixel snow-cover information quantitatively using medium and coarse spatial resolution multispectral images from the Advanced Wide Field Sensor (AWiFS) and Moderate Resolution Imaging Spectrometer (MODIS) sensors over the Himalayan region. Both the methods give slow convergence rates and inaccurate estimation of sub-pixel components analysed using root mean square (RMS) error and image deviation. Multiplicative iterative algorithms such as the Expectation Maximization Maximum Likelihood Method (EMML) and the Image Space Reconstruction Algorithm (ISRA) based on the minimization of least squares and Kullback–Leibler distances have been attempted to compute the endmembers' abundances in unmixing of satellite data. In this paper we discuss the eigenvalues of minimum noise fraction (MNF) transformation bands, data noise removal using MNF transformation and selection of pure endmembers using satellite images. The normalized difference snow index (NDSI) is also estimated using field spectral reflectance results and satellite images in green and shortwave infrared (SWIR) wavelength regions in order to carry out a comparative analysis for its variations with sub-pixel snow cover fractions. The present analysis shows the advantage of iterative over direct (constrained and unconstrained) methods; constraints are easily handled and allow better regularization of the solution for the ill-conditioned cases. Iterative methods are found to be faster compared to those of direct methods and can be used operationally for all resolution data for accurate estimation of sub-pixel snow cover.     

基于元胞自动机模型的遥感图像亚像元定位

由于遥感图像中普遍存在混合像元,因此传统分类方法得到的结果通常会存在较大误差,应用混合像元分解技术,虽然可以得到混合像元中各端元组分的丰度,但是却不能得到各端元组分的空间分布状态,而亚像元定位则是在混合像元分解的基础上,将混合像元剖分为亚像元,再利用端元组分的丰度及像元空间分布的特点,将亚像元赋予不同端元组分来得到各端元组分的空间分布情况,以提高遥感图像分类的精度。为了更好地解决亚像元定位问题,结合亚像元定位的理论模型,提出了一种新的元胞自动机模型,并通过模拟数据和实际数据对该模型进行了检验,结果表明,该模型是一种简单有效的解决亚像元定位问题的方法。     

Obtaining crop-specific time profiles of NDVI: the use of unmixing approaches for serving the continuity between SPOT-VGT and PROBA-V time series

The study examined the potential of two unmixing approaches for deriving crop-specific normalized difference vegetation index (NDVI) profiles so that upon availability of Project for On-Board Autonomy – Vegetation (PROBA-V) imagery in winter 2013, this new data set can be combined with existing Satellite Pour l’Observation de la Terre – VEGETATION (SPOT-VGT) data despite the differences in spatial resolution (300 m of PROBA-V versus 1 km of SPOT-VGT). To study the problem, two data sets were analysed: (1) a set of 10 temporal NDVI images, with 300 and 1000 m spatial resolution, from the state of São Paulo (Brazil) synthesized from 30 m Landsat Thematic Mapper (TM) images, and (2) a corresponding set of 10 observed Moderate Resolution Imaging Spectroradiometer (MODIS) images (250 m spatial resolution). To mimic the influence of noise on the retrieval accuracy, different sensor/atmospheric noise levels were applied to the first data set. For the unmixing analysis, a high-resolution land-cover (LC) map was used. The LC map was derived beforehand using a different set of Landsat TM images. The map distinguishes nine classes, with four different sugarcane stages, two agricultural sub-classes, plus forest, pasture, and urban/water. Unmixing aiming at the retrieval of crop-specific NDVI profiles was done at administrative level. For the synthesized data set it was demonstrated that the ‘true’ NDVI temporal profiles of different land-cover classes (from 30 m TM data) can generally be retrieved with high accuracy. The two simulated sensors (PROBA-V and SPOT-VGT) and the two unmixing algorithms gave similar results. Analysing the MODIS data set, we also found a good correspondence between the modelled NDVI profiles (both approaches) and the (true) Landsat temporal endmembers.     

A multiple dataset approach for 30-m resolution land cover mapping: a case study of continental Africa

Recent developments in global land-cover mapping have focused on spatial resolution improvement with more heterogeneous features to integrate spatial, spectral and temporal information. In this study, hundreds of features derived from four seasonal Landsat 8 OLI (Operational Land Imager) spectral bands, Moderate Resolution Imaging Spectroradiometer (MODIS) time series vegetation index (VI) data, night-time light (NTL), digital elevation models (DEM) and climatic variables were used for land cover mapping with a target 30-m resolution for the whole African continent. In total, 49,007 training samples (from 11,231 locations) and 23,803 validation samples (from 5,414 locations) interpreted from seasonal Landsat, MODIS Normalized Difference Vegetation Index (NDVI) time series and high-resolution images in Google Earth were used for classifier training (Random Forest) and map validation. Overall accuracy was 76% at 30-m spatial resolution, which is better than previous land cover mapping for the African continent. Besides, accuracies for cropland were improved dramatically by more than 10%. Our method also addressed many remaining issues for 30-m mapping (e.g. boundary effects and declines in resolution). This framework is promising for automatic and efficient global land cover mapping resulting in better visual effects and classification accuracy.     

遥感影像亚像元定位研究综述

遥感影像亚像元定位是在混合像元分解基础上,利用地物空间分布特征确定不同地物类型在混合像元中的具体位置,得到亚像元尺度的地物分类图,是一种有效解决混合像元空间不确定性的方法。首先介绍遥感影像亚像元定位的基本概念,分析亚像元定位的理论模型和求解算法;然后总结亚像元定位模型的误差来源、精度评价方法以及结果不确定性的表达手段,同时讨论利用辅助数据源提高亚像元定位精度的主要方法;最后对亚像元定位的研究趋势做了进一步展望。     

Analysis of multi-temporal SPOT NDVI images for small-scale land-use mapping

Land-use information is required for a number of purposes such as to address food security issues, to ensure the sustainable use of natural resources and to support decisions regarding food trade and crop insurance. Suitable land-use maps often either do not exist or are not readily available. This article presents a novel method to compile spatial and temporal land-use data sets using multi-temporal remote sensing in combination with existing data sources. Satellite Pour l'Observation de la Terre (SPOT)-Vegetation 10-day composite normalized difference vegetation index (NDVI) images (1998–2002) at 1km² resolution for a part of the Nizamabad district, Andhra Pradesh, India, were linked with available crop calendars and information about cropping patterns. The NDVI images were used to stratify the study area into map units represented by 11 distinct NDVI classes. These were then related to an existing land-cover map compiled from high resolution Indian Remote Sensing (IRS)-images (Liss-III on IRS-1C), reported crop areas by sub-district and practised crop calendar information. This resulted in an improved map containing baseline information on both land cover and land use. It is concluded that each defined NDVI class represents a varying but distinct mix of land-cover classes and that the existing land-cover map consists of too many detailed ‘year-specific’ features. Four groups of the NDVI classes present in agricultural areas match well with four categories of practised crop calendars. Differences within a group of NDVI classes reveal area specific variations in cropping intensities. The remaining groups of NDVI classes represent other land-cover complexes. The method illustrated in this article has the potential to be incorporated into remote sensing and Geographical Information System (GIS)-based drought monitoring systems.