Method For Land Cover Mapping Research Articles

Land Cover Mapping with Convolutional Neural Networks Using Sentinel-2 Images: Case Study of Rome

Land cover monitoring is crucial to understand land transformations at a global, regional and local level, and the development of innovative methodologies is necessary in order to define appropriate policies and land management practices. Deep learning techniques have recently been demonstrated as a useful method for land cover mapping through the classification of remote sensing imagery. This research aims to test and compare the predictive models created using the convolutional neural networks (CNNs) VGG16, DenseNet121 and ResNet50 on multitemporal and single-date Sentinel-2 satellite data. The most promising model was the VGG16 both with single-date and multi-temporal images, which reach an overall accuracy of 71% and which was used to produce an automatically generated EAGLE-compliant land cover map of Rome for 2019. The methodology is part of the land mapping activities of ISPRA and exploits its main products as input and support data. In this sense, it is a first attempt to develop a high-update-frequency land cover classification tool for dynamic areas to be integrated in the framework of the ISPRA monitoring activities for the Italian territory.

Learning spectral-spatial representations from VHR images for fine-scale crop type mapping: A case study of rice-crayfish field extraction in South China

Accurate information on crop type spatial distributions is crucial for yield estimation, agroecological modeling, and agrarian policy development. In comparison with moderate resolution (10–30 m) satellite images, very high-resolution (VHR, ≤ 1 m) images are superior at describing the spatial details and complex texture features of agricultural land. However, VHR images generally suffer from insufficient spatiotemporal coverage and high spectral variation in crop type, making fine-scale crop type mapping at regional or larger scales challenging. In this study, we developed a novel data-driven residual attention U-shape network (RAUNet) to identify crop types at the regional scale based on VHR images. RAUNet was adapted from UNet by introducing a residual module and a gated attention mechanism to strengthen the multilevel representation of crop features from VHR images. We selected rice-crayfish fields (RCFs) in Hubei Province, China, as the case crop type and area to test the performance of RAUNet. Seven models, including RAUNetpan, ResUNet, AttnUNet, UNet, DeepLabV3 + and random forest at the pixel and object levels, were adopted for comparison with RAUNet. Moreover, the model transferability was evaluated over regions with different agricultural landscapes. Our results showed that RAUNet-derived RCF maps obtained an average F1-score and MCC of 0.90 and 0.85 at different cultivation stages, which significantly outperformed other methods. Additionally, RAUNet had good spatiotemporal transferability, with an average F1-score and MCC of 0.93 and 0.88 in four transfer areas. The visualization results of the attention weights indicated RAUNet’s prominent capability in combining useful and relevant information from different feature levels to produce refinement predictions for the target crop type. Furthermore, the high-level features derived from RAUNet can well separate RCFs and non-RCFs over different landscapes, which were also consistent across different cultivation periods. These encouraging results suggest that RAUNet works well for fine-scale crop type mapping in smallholder farming systems, and it is an efficient method for land cover mapping based on VHR images at the regional scale.

A simple, fast, and accurate method for land cover mapping in Mongolia

Low accuracy of global land cover (LC) products at local and regional scales is concerning, because of its huge impact on downstream applications. Therefore, developing a method for high accuracy LC maps at local and regional scale is of utmost importance. Taking Mongolia as a case study, we proposed a simple, fast, and accurate method to produce annual LC maps with 250 m spatial resolution from 2001 to 2020 using MODIS data. Our products have higher spatial resolution and higher accuracy (Overall Accuracy ∼90%) compared to MCD12Q1. These new maps are critical e.g., for land degradation research, and desertification/forest cover monitoring. Especially, information on grassland ecosystems is of utmost importance for Mongolia since more than half of the country economically depends on grassland resources. Therefore, Mongolia will benefit from the new dataset. Furthermore, due to the simplicity of the method, it can be easily applied and transferred to other regions.

Integration of Optical and SAR Imagery for Dual PolSAR Features Optimization and Land Cover Mapping

Integrating optical and synthetic aperture radar (SAR) data at a feature level has proven to be a powerful method for land cover mapping. However, blurring effects and resolution loss induced by the despeckling procedure can degrade the estimation accuracy of SAR features and therefore introduce artifacts into data fusion. This article presents a method to solve this problem using SAR features estimation driven by both optical and dual polarimetric SAR images. The incentive behind this method is that compared to a low signal-to-noise ratio (SNR) SAR image, labeling sample from the optical image with higher SNR and spectral resolution can enhance sample homogeneity. The homogeneous sample following the same population can ensure SAR feature estimation accuracy and simultaneously preserve the spatial resolution. The optimum SAR features, together with optical images, are then integrated using the machine learning algorithm to obtain a classification map. We further fully evaluate the influence of SAR features estimation on fusion classification. The effectiveness of the proposed method is validated by moderate resolution Sentinel-1/2 datasets over the Saint Louis area, where abundant land cover types and rich texture present a relative challenge for classification. The quantitative analysis indicates that the proposed method can increase around 5%–10% in overall accuracy than the state-of-the-art methods. The results confirm the necessity of SAR features with high quality for land cover classification, especially over the textural region.

Assessing the effects of irrigated agricultural expansions on Lake Urmia using multi-decadal Landsat imagery and a sample migration technique within Google Earth Engine

Irrigated agricultural expansion is one of the main reasons for water scarcity in the Lake Urmia basin. Although previous studies have analyzed the impact of cropland expansion on the Lake Urmia Shrinkage, there is a lack of comprehensive annual assessment of historical irrigation expansion in the Lake Urmia basin and its impact on water resources of this region. In this study, we developed an automatic and efficient workflow using Landsat and Gravity Recovery and Climate Experiment (GRACE) data, GRACE Follow-On (GRACE-FO) data, and a sample migration technique within the Google Earth Engine cloud computing platform to comprehensively investigate the impact of irrigated agricultural expansion on the shrinkage of Lake Urmia, as one of the most severe environmental crisis in the world. Additionally, using the global surface water data, we proposed a fully automatic procedure to obtain reference samples from water bodies. The Lake Urmia basin was first classified into the water, irrigated, and Non-Water/Irrigated classes using the random forest algorithm. The average overall accuracy of the produced annual land cover maps during 1987–2020 was 92.2%, representing the great potential of the developed method for land cover mapping. We found that the irrigated lands expanded by nearly 890 km2 during the study period. Coincident with this change, although the area of water bodies in Lake Urmia partially recovered after 2015 (reached from 1,050 km2 in 2015 to 3,370 km2 in 2020), it is currently far beyond its original condition (i.e., ∼5,400 km2, average record during 1987–2000). Moreover, the information of the Terrestrial Water Storage (TWS) from the GRACE and GRACE-FO data between 2003 and 2020 showed a dramatic decrease in TWS level (∼−11.5 cm). The findings of this research will assist the local stakeholders and authorities to better understanding the environmental costs of irrigation expansion in the Lake Urmia basin.

Phenology based classification index method for land cover mapping from hyperspectral imagery

Remote sensing imagery classification contributes assistance to real-time applications for comfort and secures the society. The imagery of satellites entirely depends on the sensor type in satellites. Phenology reflection varies based on the land cover type, which absorbs external energy. Multispectral high-resolution imagery has the maximum details about the earth’s surface. This research work defines phenology based classification approach, which can produce precise high precision land cover classification. The need to develop a phenology based methodology reflects on the vegetation development classification and produces a much more suitable land cover map based on reflection values. The RGB channel values of the image do not influence this technique of reflection phenology classification. Phenology Based Classification Index (PBCI) supervised method is used to classify the high-resolution multispectral imagery with improved phenology classification methods. PBCI works on the passive sensor satellite images, without clouds and shadow in classification. The proposed method has compared with existing phenology classification methods using more than seven quality metrics.

Land cover classification in Thailand's Eastern Economic Corridor (EEC) using convolutional neural network on satellite images

Land cover plays an integral role in urban management as a source of information to support authorities’ decision making. Recently, computer vision methods, machine learning algorithms in particular, are increasingly used for land cover and land use mapping, which can help make these processes more efficient and more affordably. To this end, this paper focuses on leveraging artificial intelligence using convolutional neural network (CNN) to propose a new method for land cover mapping. A first CNN model is trained with a large number of 4 image classes to obtain a land cover model. We then apply the model directly on satellite images that are cropped into images of the same size as the training images. The results of the direct application of the model is not satisfactory; in particular, it confuses water and forested areas. The paper propose a different approach where land cover model is obtained by combining 2 models in series. The first model is a binary-class CNN model which contains classification of two classes (water, land) while the second model is a three-class CNN model containing classification of land, excluding non-water areas. Global accuracy obtained is 98% and 91% for the binary- and three-class CNN model respectively. This approach is used successfully on large area satellite images.

Unseen Land Cover Classification from High-Resolution Orthophotos Using Integration of Zero-Shot Learning and Convolutional Neural Networks

Zero-shot learning (ZSL) is an approach to classify objects unseen during the training phase and shown to be useful for real-world applications, especially when there is a lack of sufficient training data. Only a limited amount of works has been carried out on ZSL, especially in the field of remote sensing. This research investigates the use of a convolutional neural network (CNN) as a feature extraction and classification method for land cover mapping using high-resolution orthophotos. In the feature extraction phase, we used a CNN model with a single convolutional layer to extract discriminative features. In the second phase, we used class attributes learned from the Word2Vec model (pre-trained by Google News) to train a second CNN model that performed class signature prediction by using both the features extracted by the first CNN and class attributes during training and only the features during prediction. We trained and tested our models on datasets collected over two subareas in the Cameron Highlands (training dataset, first test dataset) and Ipoh (second test dataset) in Malaysia. Several experiments have been conducted on the feature extraction and classification models regarding the main parameters, such as the network’s layers and depth, number of filters, and the impact of Gaussian noise. As a result, the best models were selected using various accuracy metrics such as top-k categorical accuracy for k = [1,2,3], Recall, Precision, and F1-score. The best model for feature extraction achieved 0.953 F1-score, 0.941 precision, 0.882 recall for the training dataset and 0.904 F1-score, 0.869 precision, 0.949 recall for the first test dataset, and 0.898 F1-score, 0.870 precision, 0.838 recall for the second test dataset. The best model for classification achieved an average of 0.778 top-one, 0.890 top-two and 0.942 top-three accuracy, 0.798 F1-score, 0.766 recall and 0.838 precision for the first test dataset and 0.737 top-one, 0.906 top-two, 0.924 top-three, 0.729 F1-score, 0.676 recall and 0.790 precision for the second test dataset. The results demonstrated that the proposed ZSL is a promising tool for land cover mapping based on high-resolution photos.

Analysis of Land Cover Changing and Vegetation Index at Kuranji Watershed in Padang, West Sumatera, Indonesia

Various activities around Kuranji watershed included the land conversioncan be impacted to topographic condition and also contributed to altering the vegetation density. Remote sensing technology is an effective methodfor land cover mapping. The objectives of the present study were to analyze the changing of land cover and classifying the vegetation density index in the upstream Kuranji Watershed. This study was conducted at Kuranji Watershed in Padang, West Sumatera Province. Two Landsat images representing the changing of the watershed area during 2017 and 2018 as well as obtaining the classification of vegetation density during corresponding years.Landsat 8 OLI images were classified using a supervised classification technique, then computed the vegetation index using the Normalized Difference Vegetation Index (NDVI). The result showed that the extension of forest area, settlement area and paddy field (283.92; 35.06; and 27 Ha, respectively) and decline of mix dryland agriculture, shrub and garden area (93.68; 277.43; and 190.95 Ha respectively). Decreasing of dense vegetation found at lower dense class (6.47 Ha) and highest dense class (5535.35 Ha). Therefore, the increasing area found at the cloud, dense and higher dense class (93.17; 5525.1; and 109.94 Ha, respectively). So, it is highlighted that changing land cover and vegetation index happen during the only one-year period.

Category-Sensitive Domain Adaptation for Land Cover Mapping in Aerial Scenes

Since manually labeling aerial images for pixel-level classification is expensive and time-consuming, developing strategies for land cover mapping without reference labels is essential and meaningful. As an efficient solution for this issue, domain adaptation has been widely utilized in numerous semantic labeling-based applications. However, current approaches generally pursue the marginal distribution alignment between the source and target features and ignore the category-level alignment. Therefore, directly applying them to land cover mapping leads to unsatisfactory performance in the target domain. In our research, to address this problem, we embed a geometry-consistent generative adversarial network (GcGAN) into a co-training adversarial learning network (CtALN), and then develop a category-sensitive domain adaptation (CsDA) method for land cover mapping using very-high-resolution (VHR) optical aerial images. The GcGAN aims to eliminate the domain discrepancies between labeled and unlabeled images while retaining their intrinsic land cover information by translating the features of the labeled images from the source domain to the target domain. Meanwhile, the CtALN aims to learn a semantic labeling model in the target domain with the translated features and corresponding reference labels. By training this hybrid framework, our method learns to distill knowledge from the source domain and transfers it to the target domain, while preserving not only global domain consistency, but also category-level consistency between labeled and unlabeled images in the feature space. The experimental results between two airborne benchmark datasets and the comparison with other state-of-the-art methods verify the robustness and superiority of our proposed CsDA.

DtwSat: Time-Weighted Dynamic Time Warping for Satellite Image Time Series Analysis in R

The opening of large archives of satellite data such as LANDSAT, MODIS and the SENTINELs has given researchers unprecedented access to data, allowing them to better quantify and understand local and global land change. The need to analyze such large data sets has led to the development of automated and semi-automated methods for satellite image time series analysis. However, few of the proposed methods for remote sensing time series analysis are available as open source software. In this paper we present the R package dtwSat. This package provides an implementation of the time-weighted dynamic time warping method for land cover mapping using sequence of multi-band satellite images. Methods based on dynamic time warping are flexible to handle irregular sampling and out-of-phase time series, and they have achieved significant results in time series analysis. Package dtwSat is available from the Comprehensive R Archive Network (CRAN) and contributes to making methods for satellite time series analysis available to a larger audience. The package supports the full cycle of land cover classification using image time series, ranging from selecting temporal patterns to visualizing and assessing the results.

Mapping Inter-Annual Land Cover Variations Automatically Based on a Novel Sample Transfer Method

Most land cover mapping methods require the collection of ground reference data at the time when the remotely sensed data are acquired. Due to the high cost of repetitive collection of reference data, however, it limits the production of annual land cover maps to a short time span. In order to reduce the mapping cost and to improve the timeliness, an object-based sample transfer (OBST) method was presented in this study. The object-based analysis with strict constrains in area, shape and index values is expected to reduce the accident errors in selecting and transferring samples. The presented method was tested and compared with same-year mapping (SY), cross-year mapping (CY) and multi-index automatic classification (MI). For the study years of 2001–2016, both the overall accuracies (above 90%) and detailed accuracy indicators of the presented method were very close to the SY accuracy and higher than accuracies of CY and MI. With the presented method, the times-series land cover map of Guangzhou, China were derived and analyzed. The results reveal that the city has undergone rapid urban expansion and the pressure on natural resources and environment has increased. These results indicate the proposed method could save considerable cost and time for mapping the spatial-temporal changes of urban development. This suggests great potential for future applications as more satellite observations have become available all over the globe.

Multi-Feature Classification of Multi-Sensor Satellite Imagery Based on Dual-Polarimetric Sentinel-1A, Landsat-8 OLI, and Hyperion Images for Urban Land-Cover Classification.

This paper focuses on evaluating the ability and contribution of using backscatter intensity, texture, coherence, and color features extracted from Sentinel-1A data for urban land cover classification and comparing different multi-sensor land cover mapping methods to improve classification accuracy. Both Landsat-8 OLI and Hyperion images were also acquired, in combination with Sentinel-1A data, to explore the potential of different multi-sensor urban land cover mapping methods to improve classification accuracy. The classification was performed using a random forest (RF) method. The results showed that the optimal window size of the combination of all texture features was 9 × 9, and the optimal window size was different for each individual texture feature. For the four different feature types, the texture features contributed the most to the classification, followed by the coherence and backscatter intensity features; and the color features had the least impact on the urban land cover classification. Satisfactory classification results can be obtained using only the combination of texture and coherence features, with an overall accuracy up to 91.55% and a kappa coefficient up to 0.8935, respectively. Among all combinations of Sentinel-1A-derived features, the combination of the four features had the best classification result. Multi-sensor urban land cover mapping obtained higher classification accuracy. The combination of Sentinel-1A and Hyperion data achieved higher classification accuracy compared to the combination of Sentinel-1A and Landsat-8 OLI images, with an overall accuracy of up to 99.12% and a kappa coefficient up to 0.9889. When Sentinel-1A data was added to Hyperion images, the overall accuracy and kappa coefficient were increased by 4.01% and 0.0519, respectively.

Exploring image fusion of ALOS/PALSAR data and LANDSAT data to differentiate forest area

Remote sensing data utilize valuable information via various satellite sensors that have different specifications. Image fusion allows the user to combine different spatial and spectral resolutions to improve the information for purposes such as forest monitoring and land cover mapping. In this study, I assessed the contribution of dual-polarized Advanced Land Observing Satellite/Phased Array type L-band Synthetic Aperture Radar data to multispectral Landsat imagery. The research investigated the separability of forested areas using different image fusion techniques. Quality analysis of the fused images was conducted using qualitative and quantitative analyses. I applied the support vector machine image classification method for land cover mapping. Among all methods examined, the à trous wavelet transform method best differentiated the forested area with an overall accuracy (OA) of 94.316%, while Landsat had an OA of 92.626%. The findings of this study indicated that optical-SAR-fused images improve land cover classification, which results in higher quality forest inventory data and mapping.

Using phenological metrics and the multiple classifier fusion method to map land cover types

Feature selection and multiple classifier fusion (MCF) are effective approaches to improve land cover classification accuracy. In this study, we combined phenological metrics and the MCF method to map land cover types in Jiangsu province of China during the second crop growing season using moderate resolution imaging spectroradiometer time-series data. Eight phenological metrics were developed and calculated, and a MCF scheme was proposed by combining a simple majority vote and the measurement of posterior probabilities. The four base classifiers (i.e., the maximum likelihood classifier, the Mahalanobis distance classifier, the support vector machine classifier, and the neural networks classifier) and the MCF method were used in classifications using two spectral indices from the original satellite data (direct classification) and the computed metric data (metrics-based classification). Accuracy assessments indicated that the overall accuracies and kappa coefficients of the metrics-based classifications were all higher than those of direct classifications. The average overall accuracy and kappa coefficient of metrics-based classifications were 8.36% and 0.1 higher than that of direct classifications, respectively. Similarly, the overall accuracy and kappa coefficient of MCF generally were close to or exceeded the highest accuracy among all the base classifiers. The highest overall accuracy and kappa coefficient was achieved by classification with the MCF method based on phenological metrics (m-MCF), which were 88% and 0.85, respectively. Our results suggested that combining phenological metrics and MCF in classification is a promising method for land cover mapping in regions where strong phenological signals can be detected.

Mapping land cover gradients through analysis of hyper-temporal NDVI imagery

The green cover of the earth exhibits various spatial gradients that represent gradual changes in space of vegetation density and/or in species composition. To date, land cover mapping methods differentiate at best, mapping units with different cover densities and/or species compositions, but typically fail to express such differences as gradients. Present interpretation techniques still make insufficient use of freely available spatial-temporal Earth Observation (EO) data that allow detection of existing land cover gradients. This study explores the use of hyper-temporal NDVI imagery to detect and delineate land cover gradients analyzing the temporal behavior of NDVI values. MODIS-Terra MVC-images (250m, 16-day) of Crete, Greece, from February 2000 to July 2009 are used. The analysis approach uses an ISODATA unsupervised classification in combination with a Hierarchical Clustering Analysis (HCA). Clustering of class-specific temporal NDVI profiles through HCA resulted in the identification of gradients in landcover vegetation growth patterns. The detected gradients were arranged in a relational diagram, and mapped. Three groups of NDVI-classes were evaluated by correlating their class-specific annual average NDVI values with the field data (tree, shrub, grass, bare soil, stone, litter fraction covers). Multiple regression analysis showed that within each NDVI group, the fraction cover data were linearly related with the NDVI data, while NDVI groups were significantly different with respect to tree cover (adj. R2=0.96), shrub cover (adj. R2=0.83), grass cover (adj. R2=0.71), bare soil (adj. R2=0.88), stone cover (adj. R2=0.83) and litter cover (adj. R2=0.69) fractions. Similarly, the mean Sorenson dissimilarity values were found high and significant at confidence interval of 95% in all pairs of three NDVI groups. The study demonstrates that hyper-temporal NDVI imagery can successfully detect and map land cover gradients. The results may improve land cover assessment and aid in agricultural and ecological studies.

COMBINATION OF GENETIC ALGORITHM AND DEMPSTER-SHAFER THEORY OF EVIDENCE FOR LAND COVER CLASSIFICATION USING INTEGRATION OF SAR AND OPTICAL SATELLITE IMAGERY

Abstract. The integration of different kinds of remotely sensed data, in particular Synthetic Aperture Radar (SAR) and optical satellite imagery, is considered a promising approach for land cover classification because of the complimentary properties of each data source. However, the challenges are: how to fully exploit the capabilities of these multiple data sources, which combined datasets should be used and which data processing and classification techniques are most appropriate in order to achieve the best results. In this paper an approach, in which synergistic use of a feature selection (FS) methods with Genetic Algorithm (GA) and multiple classifiers combination based on Dempster-Shafer Theory of Evidence, is proposed and evaluated for classifying land cover features in New South Wales, Australia. Multi-date SAR data, including ALOS/PALSAR, ENVISAT/ASAR and optical (Landsat 5 TM+) images, were used for this study. Textural information were also derived and integrated with the original images. Various combined datasets were generated for classification. Three classifiers, namely Artificial Neural Network (ANN), Support Vector Machines (SVMs) and Self-Organizing Map (SOM) were employed. Firstly, feature selection using GA was applied for each classifier and dataset to determine the optimal input features and parameters. Then the results of three classifiers on particular datasets were combined using the Dempster-Shafer theory of Evidence. Results of this study demonstrate the advantages of the proposed method for land cover mapping using complex datasets. It is revealed that the use of GA in conjunction with the Dempster-Shafer Theory of Evidence can significantly improve the classification accuracy. Furthermore, integration of SAR and optical data often outperform single-type datasets.

Using a sub-pixel mapping model to improve the accuracy of landscape pattern indices

The assessment of landscape spatial patterns is a key issue in landscape management. Landscape pattern indices (LPIs) are tools appropriate for analyzing landscape spatial patterns. LPIs are often derived from raster land cover maps that are extracted from remotely sensed data through hard classification. However, pixel-based hard classification methods suffer from the mixed pixel problem (in which pixels contain more than one land cover class), making for inaccurate classification maps and LPIs. In addition, LPIs generated by hard classification methods are characterized by grain sizes (the sampling unit sizes) that limit the derived landscape pattern to a certain scale. Sub-pixel mapping (SPM) models can enable fine-scale estimation of the spatial patterns of land cover classes without requiring additional data; hence, this is an appropriate downscaling method for land cover mapping. The fraction images generated by soft classification estimate the area proportion of each land cover class within each pixel, and using these images as input enables SPM models to alleviate the mixed pixel problem. At the same time, by transforming fraction images into a finer-scaled hard classification map, SPM models can minimize the influence of grain size on LPIs calculation. In this research, simulated landscape thematic patterns that can provide different landscape spatial patterns, eight commonly used LPIs and a SPM model that maximizes the spatial dependence between neighbouring sub-pixels were applied to assess the efficiency of deriving LPIs from sub-pixel model maps. Results showed that the SPM model can more precisely characterize landscape patterns than hard classification methods can. Landscape fragmentation, class abundance, the uncertainty in SPM, and the spatial resolution of the remotely sensed data influenced LPIs derived from sub-pixel maps. The largest patch index, landscape division, and patch cohesion derived from remotely sensed data with different spatial resolutions through the SPM model were suitable for inter-comparison, whereas the patch density, mean patch area, edge density, landscape shape index, and area-weighted mean shape index derived from the sub-pixel maps were sensitive to the spatial resolution of the remotely sensed data.

Using earth observation data from multiple sources to map rare habitats in a coastal conservation area

Kenfig NNR (National Nature Reserve) is a coastal sand dune system in south Wales, UK. The site is an important location for the conservation of the fen orchid Liparis loeselii, a significant proportion of the UK population is found solely on the site. Approaches to the mapping and monitoring of the habitats at Kenfig NNR using EO (Earth Observation) meth- ods are investigated. Typical airborne EO missions over such sites produce more than a single source of EO data; these may include various optical imaging sensors with different spectral ranges, film cameras and ranging devices to measure topography. Conservation managers are thus presented with the problem of which sources of data to use when producing a land cover map of the site of interest. Using a data set gathered over the Kenfig NNR site, we investigate land cover mapping methods for conservation. The land cover types of interest typically cover small areas within a much larger site so they present a hard problem for the EO data and associated classification methods to solve. Land cover classifications produced from the data sets provide a set of competing hypotheses of land cover type for the site. Methods we use to resolve this competition between the data sets include voting methods, data fusion methods and a method utilising fuzzy logic to aggregate information. This paper is intended to act as an introduction to some of the issues involved in using EO data for habitat mapping in highly heterogeneous coastal dune environments and to present some preliminary results of the performance of each method.