Mapping the spatial distribution of winter crops at sub-pixel level using AVHRR NDVI time series and neural nets

78Citations
Citations of this article
108Readers
Mendeley users who have this article in their library.

Abstract

For large areas, it is difficult to assess the spatial distribution and inter-annualvariation of crop acreages through field surveys. Such information, however, is of greatvalue for governments, land managers, planning authorities, commodity traders andenvironmental scientists. Time series of coarse resolution imagery offer the advantage ofglobal coverage at low costs, and are therefore suitable for large-scale crop type mapping.Due to their coarse spatial resolution, however, the problem of mixed pixels has to beaddressed. Traditional hard classification approaches cannot be applied because ofsub-pixel heterogeneity. We evaluate neural networks as a modeling tool for sub-pixel cropacreage estimation. The proposed methodology is based on the assumption that differentcover type proportions within coarse pixels prompt changes in time profiles of remotelysensed vegetation indices like the Normalized Difference Vegetation Index (NDVI).Neural networks can learn the relation between temporal NDVI signatures and the soughtcrop acreage information. This learning step permits a non-linear unmixing of the temporalinformation provided by coarse resolution satellite sensors. For assessing the feasibilityand accuracy of the approach, a study region in central Italy (Tuscany) was selected. Thetask consisted of mapping the spatial distribution of winter crops abundances within 1 kmAVHRR pixels between 1988 and 2001. Reference crop acreage information for networktraining and validation was derived from high resolution Thematic Mapper/EnhancedThematic Mapper (TM/ETM+) images and official agricultural statistics. Encouragingresults were obtained demonstrating the potential of the proposed approach. For example, thespatial distribution of winter crop acreage at sub-pixel level was mapped with a cross-validatedcoefficient of determination of 0.8 with respect to the reference information from highresolution imagery. For the eight years for which reference information was available, theroot mean squared error (RMSE) of winter crop acreage at sub-pixel level was 10%. Whencombined with current and future sensors, such as MODIS and Sentinel-3, the unmixing ofAVHRR data can help in the building of an extended time series of crop distributions andcropping patterns dating back to the 80s. © 2013 by the authors; licensee MDPI, Basel, Switzerland.

Figures

  • Figure 1. Study area. (Left) AVHRR-NDVI image of the study area Tuscany. The location of Tuscany’s capital city, Florence, is indicated with a red dot (coordinates: 43°46′N/ 11°15′E). The small inlet shows the location of Tuscany within Italy. (Center) Average climatic conditions for Florence. (Right) Generalized dominant slope class map of Tuscany (source: European Soil Database, www.eusoils.jrc.ec.europa.eu).
  • Figure 2. Maps displaying the spatial distribution of area fraction images (AFI) of major CORINE land cover classes resampled to AVHRR pixel size: forest, arable land, grass- and shrubland, plantation and urban/bare land. The two pie charts in the central bottom part of the figure indicate the general land cover of Tuscany (large) and the distribution of five AFIw classes within the arable land (small). All data derived from official AGRIT statistics [30] and CORINE. The geographical location of the area is indicated in Figure 1.
  • Figure 3. NDVI profiles of winter crops (in black) and the remaining crop land (labeled “summer” crops) (in red) in Tuscany derived from pixels with a high proportion of arable land (CORINE). The white areas highlight the dekads (ten-day periods) used for modeling (7–13 and 17–25).
  • Table 1. Average land occupation (in percent) within five arable land (AFIa) classes (fractional coverages of 0–20%, 20%–40%, 80%–100%) of Tuscany. The values indicate the number of pixels in relation to the total number of pixels within the AFIa class (in percent). The relative number of AVHRR pixels for each of the five arable land AFIa classes is indicated in the first column (in parentheses) in relation to the total pixel number of the study region. The information is derived from official AGRIT statistics [30] and CORINE data.
  • Figure 4. Reference winter wheat surface data for the region of Tuscany, Italy (source: AGRIT [30]). Large squares indicate the eight years for which high resolution TM/ETM+ imagery is available.
  • Figure 5. Principle functioning of a three-layer neural net as used for the current study. In the illustration, eight input nodes are shown representing four fractional land cover classes (AFI1 to AFI4) and four vegetation index values (VI7 to VI10). The hidden layer consists of three nodes receiving weighted inputs from the input nodes. The output node (target) outputs the fractional winter crop coverage for the provided input. After calculation of the weighted sums, the resulting values are passed through a (tan-sigmoid) transfer function. Weights and biases are optimized during the training by presenting the net with samples of input and target values and minimizing the error between reference and simulated target value. The current study used a very similar net with, however, twenty-one input nodes (five broad land cover classes plus 16 vegetation index values).
  • Figure 6. Inter-annual variability of winter crop acreages for AVHRR pixels with arable land. The map is derived from official AGRIT statistics and multi-annual (eight yrs) Landsat TM/ETM+ data. The inter-annual variability is expressed in standard deviations and grouped into three classes: low variability (STD < 5%), medium (5% < STD < 10%) and high variability (STD > 10%). Pixels without arable land (<2%) are shown in gray.
  • Table 2. Overview describing the experimental approach for mapping sub-pixel winter crop acreages from NOAA-AVHRR data. For training and evaluation, a dataset covering eight years was used. Independent estimates were obtained by means of cross-validation.

References Powered by Scopus

A simple method for reconstructing a high-quality NDVI time-series data set based on the Savitzky-Golay filter

1956Citations
N/AReaders
Get full text

A perfect smoother

1065Citations
N/AReaders
Get full text

Introduction neural networks in remote sensing

873Citations
N/AReaders
Get full text

Cited by Powered by Scopus

Advances in remote sensing of agriculture: Context description, existing operational monitoring systems and major information needs

792Citations
N/AReaders
Get full text

Remote sensing of ecosystem services: A systematic review

268Citations
N/AReaders
Get full text

Using low resolution satellite imagery for yield prediction and yield anomaly detection

253Citations
N/AReaders
Get full text

Register to see more suggestions

Mendeley helps you to discover research relevant for your work.

Already have an account?

Cite

CITATION STYLE

APA

Atzberger, C., & Rembold, F. (2013). Mapping the spatial distribution of winter crops at sub-pixel level using AVHRR NDVI time series and neural nets. Remote Sensing, 5(3), 1335–1354. https://doi.org/10.3390/rs5031335

Readers over time

‘12‘13‘14‘15‘16‘17‘18‘19‘20‘21‘22‘23‘24‘2506121824

Readers' Seniority

Tooltip

PhD / Post grad / Masters / Doc 49

63%

Researcher 20

26%

Professor / Associate Prof. 5

6%

Lecturer / Post doc 4

5%

Readers' Discipline

Tooltip

Earth and Planetary Sciences 25

35%

Environmental Science 25

35%

Agricultural and Biological Sciences 16

22%

Engineering 6

8%

Save time finding and organizing research with Mendeley

Sign up for free
0