Per-pixel vs. object-based classification of urban land cover extraction using high spatial resolution imagery

In using traditional digital classification algorithms, a researcher typically encounters serious issues in identifying urban land cover classes employing high resolution data. A normal approach is to use spectral information alone and ignore spatial information and a group of pixels that need to be considered together as an object. We used QuickBird image data over a central region in the city of Phoenix, Arizona to examine if an object-based classifier can accurately identify urban classes. To demonstrate if spectral information alone is practical in urban classification, we used spectra of the selected classes from randomly selected points to examine if they can be effectively discriminated. The overall accuracy based on spectral information alone reached only about 63.33%. We employed five different classification procedures with the object-based paradigm that separates spatially and spectrally similar pixels at different scales. The classifiers to assign land covers to segmented objects used in the study include membership functions and the nearest neighbor classifier. The object-based classifier achieved a high overall accuracy (90.40%), whereas the most commonly used decision rule, namely maximum likelihood classifier, produced a lower overall accuracy (67.60%). This study demonstrates that the object-based classifier is a significantly better approach than the classical per-pixel classifiers. Further, this study reviews application of different parameters for segmentation and classification, combined use of composite and original bands, selection of different scale levels, and choice of classifiers. Strengths and weaknesses of the object-based prototype are presented and we provide suggestions to avoid or minimize uncertainties and limitations associated with the approach.     

Object-based urban detailed land cover classification with high spatial resolution IKONOS imagery

Improvement in remote sensing techniques in spatial/spectral resolution strengthens their applicability for urban environmental study. Unfortunately, high spatial resolution imagery also increases internal variability in land cover units and can cause a ‘salt-and-pepper’ effect, resulting in decreased accuracy using pixel-based classification results. Region-based classification techniques, using an image object (IO) rather than a pixel as a classification unit, appear to hold promise as a method for overcoming this problem. Using IKONOS high spatial resolution imagery, we examined whether the IO technique could significantly improve classification accuracy compared to the pixel-based method when applied to urban land cover mapping in Tampa Bay, FL, USA. We further compared the performance of an artificial neural network (ANN) and a minimum distance classifier (MDC) in urban detailed land cover classification and evaluated whether the classification accuracy was affected by the number of extracted IO features. Our analysis methods included IKONOS image data calibration, data fusion with the pansharpening (PS) process, Hue–Intensity–Saturation (HIS) transferred indices and textural feature extraction, and feature selection using a stepwise discriminant analysis (SDA). The classification results were evaluated with visually interpreted data from high-resolution (0.3 m) digital aerial photographs. Our results indicate a statistically significant difference in classification accuracy between pixel- and object-based techniques; ANN outperforms MDC as an object-based classifier; and the use of more features (27 vs. 9 features) increases the IO classification accuracy, although the increase is statistically significant for the MDC but not for the ANN.     

Object-based image classification for burned area mapping of Creus Cape, Spain, using NOAA-AVHRR imagery

Due to the ability of the NOAA-AVHRR sensor to cover a wide area and its high temporal frequency, it is possible to quickly obtain a general overview of the prevailing situation over a large area of terrain and, more specifically, quickly assess the damage caused by a recent large forest fire by mapping the extent of the burned area. The aim of this work was to map a large forest fire that recently took place on the Spanish Mediterranean coast using innovative image classification techniques and low spatial resolution imagery. The methodology involved developing an object-based classification model using spectral as well as contextual object information. The burned area map resulting from the image classification was compared with the fire perimeter provided by the Catalan Environmental Department in terms of spatial overlap and size in order to determine to what extent they were compatible. Results of the comparison indicated a high degree (≈90%) of spatial agreement. The total burned area of the classified image was found to be 6900 ha, compared to a fire perimeter of 6000 ha produced by the Catalan Environmental Department. It was concluded that, although the object-oriented classification approach was capable of affording very promising results when mapping a recent burn on the Spanish Mediterranean coast, the method in question required further assessment to ascertain its ability to map other burned areas in the Mediterranean.     

Fully fuzzy supervised classification of land cover from remotely sensed imagery with an artificial neural network

Mixed pixels are a major problem in mapping land cover from remotely sensed imagery. Unfortunately, such imagery may be dominated by mixed pixels, and the conventional hard image classification techniques used in mapping applications are unable to appropriately represent the land cover of mixed pixels. Fuzzy classification techniques can, however, accommodate the partial and multiple class membership of mixed pixels, and be used to derive an appropriate land cover representation. This is, however, only a partial solution to the mixed pixel problem in supervised image classification. It must be reognised that the land cover on the ground is fuzzy, at the scale of the pixel, and so it is inappropriate to use procedures designed for hard data in the training and testing stages of the classification. Here an approach for land cover classification in which fuzziness is accommodated in all three stages of a supervised classification is presented. Attention focuses on the classification of airborne thematic mapper data with an artificial neural network. Mixed pixels could be accommodated in training the artificial neural network, since the desired output for each training pixel can be specified. A fuzzy land cover representation was derived by outputting the activation level of the network's output units. The activation level of each output unit was significantly correlated with the proportion of the area represented by a pixel which was covered with the class associated with the unit (r>0.88, significant at the 99% level of confidence). Finally, the distance between the fuzzy land cover classification derived from the artificial neural network and the fuzzy ground data was used to illustrate the accuracy of the land cover representation derived. The dangers of hardening the classification output and ground data sets to enable a conventional assessment of classification accuracy are also illustrated; the hardened data sets were over three times more distant from each other than the fuzzy data sets.     

Contribution of multispectral and multitemporal information from MODIS images to land cover classification

The goal of this study is to evaluate the relative usefulness of high spectral and temporal resolutions of MODIS imagery data for land cover classification. In particular, we highlight the individual and combinatorial influence of spectral and temporal components of MODIS reflectance data in land cover classification. Our study relies on an annual time series of twelve MODIS 8-days composited images (MOD09A1) monthly acquired during the year 2000, at a 500 m nominal resolution. As our aim is not to propose an operational classifier directed at thematic mapping based on the most efficient combination of reflectance inputs — which will probably change across geographical regions and with different land cover nomenclatures — we intentionally restrict our experimental framework to continental Portugal. Because our observation data stream contains highly correlated components, we need to rank the temporal and the spectral features according not only to their individual ability at separating the land cover classes, but also to their differential contribution to the existing information. To proceed, we resort to the median Mahalanobis distance as a statistical separability criterion. Once achieved this arrangement, we strive to evaluate, in a classification perspective, the gain obtained when the dimensionality of the input feature space grows. We then successively embedded the prior ranked measures into the multitemporal and multispectral training data set of a Support Vector Machines (SVM) classifier. In this way, we show that, only the inclusion of the approximately first three dates substantially increases the classification accuracy. Moreover, this multitemporal factor has a significant effect when coupled with combinations of few spectral bands, but it turns negligible as soon as the full spectral information is exploited. Regarding the multispectral factor, its beneficence on classification accuracy remains more constant, regardless of the number of dates being used.     

A linearised pixel-swapping method for mapping rural linear land cover features from fine spatial resolution remotely sensed imagery

Accurate maps of rural linear land cover features, such as paths and hedgerows, would be useful to ecologists, conservation managers and land planning agencies. Such information might be used in a variety of applications (e.g., ecological, conservation and land management applications). Based on the phenomenon of spatial dependence, sub-pixel mapping techniques can be used to increase the spatial resolution of land cover maps produced from satellite sensor imagery and map such features with increased accuracy. Aerial photography with a spatial resolution of 0.25 m was acquired of the Christchurch area of Dorset, UK. The imagery was hard classified using a simple Mahalanobis distance classifier and the classification degraded to simulate land cover proportion images with spatial resolutions of 2.5 and 5 m. A simple pixel-swapping algorithm was then applied to each of the proportion images. Sub-pixels within pixels were swapped iteratively until the spatial correlation between neighbouring sub-pixels for the entire image was maximised. Visual inspection of the super-resolved output showed that prediction of the position and dimensions of hedgerows was comparable with the original imagery. The maps displayed an accuracy of 87%. To enhance the prediction of linear features within the super-resolved output, an anisotropic modelling component was added. The direction of the largest sums of proportions was calculated within a moving window at the pixel level. The orthogonal sum of proportions was used in estimating the anisotropy ratio. The direction and anisotropy ratio were then used to modify the pixel-swapping algorithm so as to increase the likelihood of creating linear features in the output map. The new linear pixel-swapping method led to an increase in the accuracy of mapping fine linear features of approximately 5% compared with the conventional pixel-swapping method.     

Updating the 2001 National Land Cover Database land cover classification to 2006 by using Landsat imagery change detection methods

The recent release of the U.S. Geological Survey (USGS) National Land Cover Database (NLCD) 2001, which represents the nation's land cover status based on a nominal date of 2001, is widely used as a baseline for national land cover conditions. To enable the updating of this land cover information in a consistent and continuous manner, a prototype method was developed to update land cover by an individual Landsat path and row. This method updates NLCD 2001 to a nominal date of 2006 by using both Landsat imagery and data from NLCD 2001 as the baseline. Pairs of Landsat scenes in the same season in 2001 and 2006 were acquired according to satellite paths and rows and normalized to allow calculation of change vectors between the two dates. Conservative thresholds based on Anderson Level I land cover classes were used to segregate the change vectors and determine areas of change and no-change. Once change areas had been identified, land cover classifications at the full NLCD resolution for 2006 areas of change were completed by sampling from NLCD 2001 in unchanged areas. Methods were developed and tested across five Landsat path/row study sites that contain several metropolitan areas including Seattle, Washington; San Diego, California; Sioux Falls, South Dakota; Jackson, Mississippi; and Manchester, New Hampshire. Results from the five study areas show that the vast majority of land cover change was captured and updated with overall land cover classification accuracies of 78.32%, 87.5%, 88.57%, 78.36%, and 83.33% for these areas. The method optimizes mapping efficiency and has the potential to provide users a flexible method to generate updated land cover at national and regional scales by using NLCD 2001 as the baseline.     

Object-based analysis and change detection of major wetland cover types and their classification uncertainty during the low water period at Poyang Lake, China

Productive wetland systems at land-water interfaces that provide unique ecosystem services are challenging to study because of water dynamics, complex surface cover and constrained field access. We applied object-based image analysis and supervised classification to four 32-m Beijing-1 microsatellite images to examine broad-scale surface cover composition and its change during November 2007-March 2008 low water season at Poyang Lake, the largest freshwater lake-wetland system in China (> 4000 km²). We proposed a novel method for semi-automated selection of training objects in this heterogeneous landscape using extreme values of spectral indices (SIs) estimated from satellite data. Dynamics of the major wetland cover types (Water, Mudflat, Vegetation and Sand) were investigated both as transitions among primary classes based on maximum membership value, and as changes in memberships to all classes even under no change in a primary class. Fuzzy classification accuracy was evaluated as match frequencies between classification outcome and a) the best reference candidate class (MAX function) and b) any acceptable reference class (RIGHT function). MAX-based accuracy was relatively high for Vegetation (≥ 90%), Water (≥ 82%), Mudflat (≥ 76%) and the smallest-area Sand (≥ 75%) in all scenes; these scores improved with the RIGHT function to 87-100%. Classification uncertainty assessed as the proportion of fuzzy object area within a class at a given fuzzy threshold value was the highest for all classes in November 2007, and consistently higher for Mudflat than for other classes in all scenes. Vegetation was the dominant class in all scenes, occupying 41.2-49.3% of the study area. Object memberships to Vegetation mostly declined from November 2007 to February 2008 and increased substantially only in February-March 2008, possibly reflecting growing season conditions and grazing. Spatial extent of Water both declined and increased during the study period, reflecting precipitation and hydrological events. The “fuzziest” Mudflat class was involved in major detected transitions among classes and declined in classification accuracy by March 2008, representing a key target for finer-scale research. Future work should introduce Vegetation sub-classes reflecting differences in phenology and alternative methods to discriminate Mudflat from other classes. Results can be used to guide field sampling and top-down landscape analyses in this wetland.     

Spatial-temporal dynamics of land surface temperature in relation to fractional vegetation cover and land use/cover in the Tabriz urban area, Iran

Rapid changes of land use and land cover (LULC) in urban areas have become a major environmental concern due to environmental impacts, such as the reduction of green spaces and development of urban heat islands (UHI). Monitoring and management plans are required to solve this problem effectively. The Tabriz metropolitan area in Iran, selected as a case study for this research, is an example of a fast growing city. Multi-temporal images acquired by Landsat 4, 5 TM and Landsat 7 ETM+ sensors on 30 June 1989, 18 August 1998, and 2 August 2001 respectively, were corrected for radiometric and geometric errors, and processed to extract LULC classes and land surface temperature (LST). The relationship between temporal dynamics of LST and LULC was then examined. The temperature vegetation index (TVX) space was constructed in order to study the temporal variability of thermal data and vegetation cover. Temporal trajectory of pixels in the TVX space showed that most changes due to urbanization were observable as the pixels migrated from the low temperature-dense vegetation condition to the high temperature-sparse vegetation condition in the TVX space. The uncertainty analysis revealed that the trajectory analysis in the TVX space involved a class-dependant noise component. This emphasized the need for multiple LULC control points in the TVX space. In addition, this research suggests that the use of multi-temporal satellite data together with the examination of changes in the TVX space is effective and useful in urban LULC change monitoring and analysis of urban surface temperature conditions as long as the uncertainty is addressed.     

Identification and analysis of urban surface temperature patterns in Greater Athens, Greece, using MODIS imagery

Thermal infrared images are being acquired by satellites for more than two decades enabling studies of the human-induced Urban Heat Island (UHI) phenomenon. As a result, the requirement of the scientific community for fast and efficient methods for extracting and analyzing the thermal patterns from a vast volume of acquired data has emerged. The present paper proposes an innovative object-based image analysis procedure to extract thermal patterns for the quantitative analysis of satellite-derived Land Surface Temperature (LST) maps. The spatial and thermal attributes associated with these objects are then calculated and used for the analyses of the intensity, the position and the spatial extent of UHIs. A case study was conducted in the Greater Athens Area, Greece. More than 3000 LST images of the area acquired by MODIS sensor over a decade were analyzed. Three daytime hot-spots were identified and studied (Megara, Elefsina-Aspropyrgos and Mesogeia). They were all found to exhibit similar behavior, gradually increasing their maximum temperature during the summer season and reaching their maxima in mid-July. The hot-spots' thermal intensities compared to a suburban area were of 9-10 °C and were found to be highly correlated to their areal extent. During the night-time, Athens center developed a typical UHI spatially coinciding with the dense urban fabric. The nighttime maximum LST peaked (on average) at the end of July, two weeks later than the daytime surface patterns. The mean spatial extent of UHI in Athens was 55.2 km², whilst its mean intensity was 5.6 °C. The proposed automatic extraction process can be customized for other cities and potentially used for comparison of LST patterns and UHI behavior between different cities.     

Land cover mapping of large areas using chain classification of neighboring Landsat satellite images

Satellite imagery is the major data source for regional to global land cover maps. However, land cover mapping of large areas with medium-resolution imagery is costly and often constrained by the lack of good training and validation data. Our goal was to overcome these limitations, and to test chain classifications, i.e., the classification of Landsat images based on the information in the overlapping areas of neighboring scenes. The basic idea was to classify one Landsat scene first where good ground truth data is available, and then to classify the neighboring Landsat scene using the land cover classification of the first scene in the overlap area as training data. We tested chain classification for a forest/non-forest classification in the Carpathian Mountains on one horizontal chain of six Landsat scenes, and two vertical chains of two Landsat scenes each. We collected extensive training data from Quickbird imagery for classifying radiometrically uncorrected data with Support Vector Machines (SVMs). The SVMs classified 8 scenes with overall accuracies between 92.1% and 98.9% (average of 96.3%). Accuracy loss when automatically classifying neighboring scenes with chain classification was 1.9% on average. Even a chain of six images resulted only in an accuracy loss of 5.1% for the last image compared to a reference classification from independent training data for the last image. Chain classification thus performed well, but we note that chain classification can only be applied when land cover classes are well represented in the overlap area of neighboring Landsat scenes. As long as this constraint is met though, chain classification is a powerful approach for large area land cover classifications, especially in areas of varying training data availability.     

Cross-border comparison of land cover and landscape pattern in Eastern Europe using a hybrid classification technique

Eastern Europe has experienced drastic changes in political and economic conditions following the breakdown of the Soviet Union. Furthermore, these changes often differ among neighboring countries. This offers unique possibilities to assess the relative importance of broad-scale political and socioeconomic factors on land cover and landscape pattern. Our question was how much land cover differed in the Polish, the Slovak, and the Ukrainian Carpathian Mountains and to what extent these differences can be related to dissimilarities in societal, economic, and political conditions. We used a hybrid classification technique, combining advantages from supervised and unsupervised methods, to derive a land cover map from three Landsat Thematic Mapper (TM) and Enhanced Thematic Mapper Plus (ETM+) images from 2000. Results showed marked differences in land cover between the three countries. Forest cover and composition was different for the three countries, for example Slovakia and Poland had about 20% more forest cover at higher elevations than Ukraine. Broadleaved forest dominated in Slovakia while high percentages of conifers were found in Poland and Ukraine. Agriculture was most abundant in Slovakia where the lowest level of agricultural fragmentation was found (22% core area compared to less than 5% in Poland and Ukraine). Post-socialist land change was greatest in Ukraine, were we found high agricultural fragmentation and widespread early-successional shrublands indicating extensive land abandonment. Concerning forests, differences can largely be explained by socialist forest management. The abundance and pattern of arable land and grassland can be explained by two factors: land tenure in socialist times and economic transition since 1990. These results suggest that broad-scale socioeconomic and political factors are of major significance for land cover patterns in Eastern Europe, and possibly elsewhere.     

Sub-pixel mapping of urban land cover using multiple endmember spectral mixture analysis: Manaus, Brazil

The spatial and spectral variability of urban environments present fundamental challenges to deriving accurate remote sensing products for urban areas. Multiple endmember spectral mixture analysis (MESMA) is a technique that potentially addresses both challenges. MESMA models spectra as the linear sum of spectrally pure endmembers that vary on a per-pixel basis. Spatial variability is addressed by mapping sub-pixel components of land cover as a combination of endmembers. Spectral variability is addressed by allowing the number and type of endmembers to vary from pixel to pixel. This paper presents an application of MESMA to map the physical components of urban land cover for the city of Manaus, Brazil, using Landsat Enhanced Thematic Mapper (ETM+) imagery.We present a methodology to build a regionally specific spectral library of urban materials based on generalized categories of urban land-cover components: vegetation, impervious surfaces, soil, and water. Using this library, we applied MESMA to generate a total of 1137 two-, three-, and four-endmember models for each pixel; the model with the lowest root-mean-squared (RMS) error and lowest complexity was selected on a per-pixel basis. Almost 97% of the pixels within the image were modeled within the 2.5% RMS error constraint. The modeled fractions were used to generate continuous maps of the per-pixel abundance of each generalized land-cover component. We provide an example to demonstrate that land-cover components have the potential to characterize trajectories of physical landscape change as urban neighborhoods develop through time. Accuracy of land-cover fractions was assessed using high-resolution, geocoded images mosaicked from digital aerial videography. Modeled vegetation and impervious fractions corresponded well with the reference fractions. Modeled soil fractions did not correspond as closely with the reference fractions, in part due to limitations of the reference data. This work demonstrates the potential of moderate-resolution, multispectral imagery to map and monitor the evolution of the physical urban environment.     

基于高分辨率遥感影像的面向对象城市土地覆被分类比较研究

针对高分辨率遥感影像的城市土地覆被信息提取,根据分类目的与精度要求的不同,分别引入了优化与广义两种面向对象分类方案,并对分类的结果进行分析比较。结果表明:①优化方案的分类结果总体上要比广义方案好,前者的总体精度为86.50%,相比后者的80.50%提高了6.0%,而总体Kappa系数提高了0.0851,但是该方案效率低,可移植性差;②广义方案的分类结果虽然精度略低,但是该方案具有很强的适用性与可移植性,能够在精度可控范围内,很大程度提高分类效率,实现系统而有效的自动分类;③广义方案得到的分类结果具有一致的精度,在利用其建立城市生态模型中能够保证数据之间的系统性与鲁棒性。因此,利用优化方案能够提高分类结果的绝对精度,而广义方案对于实时精确获取城市土地覆被信息、小尺度上定量监测与评价城市化的生态后果以及有效开展城市土地规划与管理具有更重要的意义。     

ECOCLIMAP-II: An ecosystem classification and land surface parameters database of Western Africa at 1 km resolution for the African Monsoon Multidisciplinary Analysis (AMMA) project

This work is devoted to a presentation of the ECOCLIMAP-II database for Western Africa, which is an upgrade for this region of the former initiative, ECOCLIMAP-I, implemented at global scale. ECOCLIMAP-II is a dual database at 1-km resolution that comprises an ecosystem classification and a coherent set of land surface parameters. This new physiographic information (e.g. leaf area index, fractional vegetation cover, albedo and land cover classification), was especially developed in the framework of the African Monsoon Multidisciplinary Analysis (AMMA) programme in order to support the modelling of land-atmosphere interactions, which stresses the importance of the present study. Criteria for coherence between prevalent land cover classifications and the analysis of time series of the satellite leaf area index (LAI) between 2000 and 2007 constitute the analysis tools for setting up ECOCLIMAP-II. The LAI and inferred fraction of vegetation cover are spatially distributed per land cover unit. The fraction of vegetation cover is handled to split the land surface albedo into vegetation and bare soil albedo components, as is required for a large number of applications. The new ECOCLIMAP-II land cover product is improved with regard to the spatial coherence compared to former version. The reliability of the physiographic details is also confirmed through verification with land cover products at higher resolution.     

Integrated modelling of urban spatial development under uncertain climate futures: A case study in Hungary

To provide fundamental decision support information for climate risk assessment in Hungary, an urban spatial development model of land cover change and population age structure dynamics was developed and applied to local integrated scenarios of climate change and stakeholder-derived socio-economic change. The four integrated scenarios for Hungary produced contrasting projections for urban patterns to 2100, but peri-urbanisation around Budapest was estimated to occur under all scenarios, together with a decline in working age population in the centres of the capital and major towns. This suggests that future urban planning needs to take into consideration the potential for underutilised urban infrastructure in the centre of the capital and pressures for social service provisioning in its outskirt. The integrated scenarios and model developed can be used in future studies to test the effectiveness of inter-sectoral policy responses in adapting urban planning to multiple climate and socio-economic challenges.     

Comparisons of land cover and LAI estimates derived from ETM+ and MODIS for four sites in North America: a quality assessment of 2000/2001 provisional MODIS products

The MODIS land science team produces a number of standard products, including land cover and leaf area index (LAI). Critical to the success of MODIS and other sensor products is an independent evaluation of product quality. In that context, we describe a study using field data and Landsat ETM+ to map land cover and LAI at four 49-km² sites in North America containing agricultural cropland (AGRO), prairie grassland (KONZ), boreal needleleaf forest, and temperate mixed forest. The purpose was to: (1) develop accurate maps of land cover, based on the MODIS IGBP (International Geosphere-Biosphere Programme) land cover classification scheme; (2) derive continuous surfaces of LAI that capture the mean and variability of the LAI field measurements; and (3) conduct initial MODIS validation exercises to assess the quality of early (i.e., provisional) MODIS products. ETM+ land cover maps varied in overall accuracy from 81% to 95%. The boreal forest was the most spatially complex, had the greatest number of classes, and the lowest accuracy. The intensive agricultural cropland had the simplest spatial structure, the least number of classes, and the highest overall accuracy. At each site, mapped LAI patterns generally followed patterns of land cover across the site. Predicted versus observed LAI indicated a high degree of correspondence between field-based measures and ETM+ predictions of LAI. Direct comparisons of ETM+ land cover maps with Collection 3 MODIS cover maps revealed several important distinctions and similarities. One obvious difference was associated with image/map resolution. ETM+ captured much of the spatial complexity of land cover at the sites. In contrast, the relatively coarse resolution of MODIS did not allow for that level of spatial detail. Over the extent of all sites, the greatest difference was an overprediction by MODIS of evergreen needleleaf forest cover at the boreal forest site, which consisted largely of open shrubland, woody savanna, and savanna. At the agricultural, temperate mixed forest, and prairie grassland sites, ETM+ and MODIS cover estimates were similar. Collection 3 MODIS-based LAI estimates were considerably higher (up to 4 m² m⁻²) than those based on ETM+ LAI at each site. There are numerous probable reasons for this, the most important being the algorithms' sensitivity to MODIS reflectance calibration, its use of a prelaunch AVHRR-based land cover map, and its apparent reliance on mainly red and near-IR reflectance. Samples of Collection 4 LAI products were examined and found to consist of significantly improved LAI predictions for KONZ, and to some extent for AGRO, but not for the other two sites. In this study, we demonstrate that MODIS reflectance data are highly correlated with LAI across three study sites, with relationships increasing in strength from 500 to 1000 m spatial resolution, when shortwave-infrared bands are included.     

A scalable approach to mapping annual land cover at 250 m using MODIS time series data: A case study in the Dry Chaco ecoregion of South America

Land use and land cover (LULC) maps from remote sensing are vital for monitoring, understanding and predicting the effects of complex human-nature interactions that span local, regional and global scales. We present a method to map annual LULC at a regional spatial scale with source data and processing techniques that permit scaling to broader spatial and temporal scales, while maintaining a consistent classification scheme and accuracy. Using the Dry Chaco ecoregion in Argentina, Bolivia and Paraguay as a test site, we derived a suite of predictor variables from 2001 to 2007 from the MODIS 250 m vegetation index product (MOD13Q1). These variables included: annual statistics of red, near infrared, and enhanced vegetation index (EVI), phenological metrics derived from EVI time series data, and slope and elevation. For reference data, we visually interpreted percent cover of eight classes at locations with high-resolution QuickBird imagery in Google Earth. An adjustable majority cover threshold was used to assign samples to a dominant class. When compared to field data, we found this imagery to have georeferencing error < 5% the length of a MODIS pixel, while most class interpretation error was related to confusion between agriculture and herbaceous vegetation. We used the Random Forests classifier to identify the best sets of predictor variables and percent cover thresholds for discriminating our LULC classes. The best variable set included all predictor variables and a cover threshold of 80%. This optimal Random Forests was used to map LULC for each year between 2001 and 2007, followed by a per-pixel, 3-year temporal filter to remove disallowed LULC transitions. Our sequence of maps had an overall accuracy of 79.3%, producer accuracy from 51.4% (plantation) to 95.8% (woody vegetation), and user accuracy from 58.9% (herbaceous vegetation) to 100.0% (water). We attributed map class confusion to limited spectral information, sub-pixel spectral mixing, georeferencing error and human error in interpreting reference samples. We used our maps to assess woody vegetation change in the Dry Chaco from 2002 to 2006, which was characterized by rapid deforestation related to soybean and planted pasture expansion. This method can be easily applied to other regions or continents to produce spatially and temporally consistent information on annual LULC.     

Dissimilarity-based classification in the absence of local ground truth: Application to the diagnostic interpretation of chest radiographs

In this paper classification on dissimilarity representations is applied to medical imaging data with the task of discrimination between normal images and images with signs of disease. We show that dissimilarity-based classification is a beneficial approach in dealing with weakly labeled data, i.e. when the location of disease in an image is unknown and therefore local feature-based classifiers cannot be trained. A modification to the standard dissimilarity-based approach is proposed that makes a dissimilarity measure multi-valued, hence, able to retain more information. A multi-valued dissimilarity between an image and a prototype becomes an image representation vector in classification. Several classification outputs with respect to different prototypes are further integrated into a final image decision. Both standard and proposed methods are evaluated on data sets of chest radiographs with textural abnormalities and compared to several feature-based region classification approaches applied to the same data. On a tuberculosis data set the multi-valued dissimilarity-based classification performs as well as the best region classification method applied to the fully labeled data, with an area under the receiver operating characteristic (ROC) curve (A_z) of 0.82. The standard dissimilarity-based classification yields A_z=0.80. On a data set with interstitial abnormalities both dissimilarity-based approaches achieve A_z=0.98 which is closely behind the best region classification method.