text: Remote Sensing http://www.mdpi.com/journal/remotesensing Latest open access articles published in Remote Sens. at http://www.mdpi.com/journal/remotesensing http://www.mdpi.com/2072-4292/7/7/8543 The rapidly developing urbanization since the last decade of the 20th century has led to extensive groundwater extraction, resulting in subsidence in Ho Chi Minh City, Vietnam. Recent advances in multi-temporal spaceborne SAR interferometry, especially with a persistent scatters interferometry (PSI) approach, has made this a robust remote sensing technique for measuring large-scale ground subsidence with millimetric accuracy. This work has presented an advanced PSI analysis, to provide an unprecedented spatial extent and continuous temporal coverage of the subsidence in Ho Chi Minh City from 2006 to 2010. The study shows that subsidence is most severe in the Holocene silt loam areas along the Sai Gon River and in the southwest of the city. The groundwater extraction resulting from urbanization and urban growth is mainly responsible for the subsidence. Subsidence in turn leads to more flooding and water nuisance. The correlation between the reference leveling velocity and the estimated PSI result is R2 = 0.88, and the root mean square error is 4.3 (mm/year), confirming their good agreement. From 2006 to 2010, the estimation of the average subsidence rate is -8.0 mm/year, with the maximum value up to -70 mm/year. After four years, in regions along Sai Gon River and in the southwest of the city, the land has sunk up to -12 cm. If not addressed, subsidence leads to the increase of inundation, both in frequency and spatial extent. Finally, regarding climate change, the effects of subsidence should be considered as appreciably greater than those resulting from rising sea level. It is essential to consider these two factors, because the city is inhabited by more than 7.5 million people, where subsidence directly impacts urban structures and infrastructure. Remote Sensing 2015-07-06 7 7 Article 10.3390/rs70708543 8543 8562 2072-4292 2015-07-06 doi: 10.3390/rs70708543 Dinh Minh Le Van Trung Thuy Toan http://www.mdpi.com/2072-4292/7/7/8516 River deltas belong to the most densely settled places on earth. Although they only account for 5% of the global land surface, over 550 million people live in deltas. These preferred livelihood locations, which feature flat terrain, fertile alluvial soils, access to fluvial and marine resources, a rich wetland biodiversity and other advantages are, however, threatened by numerous internal and external processes. Socio-economic development, urbanization, climate change induced sea level rise, as well as flood pulse changes due to upstream water diversion all lead to changes in these highly dynamic systems. A thorough understanding of a river delta's general setting and intra-annual as well as long-term dynamic is therefore crucial for an informed management of natural resources. Here, remote sensing can play a key role in analyzing and monitoring these vast areas at a global scale. The goal of this study is to demonstrate the potential of intra-annual time series analyses at dense temporal, but coarse spatial resolution for inundation characterization in five river deltas located in four different countries. Based on 250 m MODIS reflectance data we analyze inundation dynamics in four densely populated Asian river deltas—namely the Yellow River Delta (China), the Mekong Delta (Vietnam), the Irrawaddy Delta (Myanmar), and the Ganges-Brahmaputra (Bangladesh, India)—as well as one very contrasting delta: the nearly uninhabited polar Mackenzie Delta Region in northwestern Canada for the complete time span of one year (2013). A complex processing chain of water surface derivation on a daily basis allows the generation of intra-annual time series, which indicate inundation duration in each of the deltas. Our analyses depict distinct inundation patterns within each of the deltas, which can be attributed to processes such as overland flooding, irrigation agriculture, aquaculture, or snowmelt and thermokarst processes. Clear differences between mid-latitude, subtropical, and polar deltas are illustrated, and the advantages and limitations of the approach for inundation derivation are discussed. Remote Sensing 2015-07-06 7 7 Article 10.3390/rs70708516 8516 8542 2072-4292 2015-07-06 doi: 10.3390/rs70708516 Claudia Kuenzer Igor Klein Tobias Ullmann Efi Georgiou Roland Baumhauer Stefan Dech http://www.mdpi.com/2072-4292/7/7/8489 Random Forest (RF) is a widely used algorithm for classification of remotely sensed data. Through a case study in peatland classification using LiDAR derivatives, we present an analysis of the effects of input data characteristics on RF classifications (including RF out-of-bag error, independent classification accuracy and class proportion error). Training data selection and specific input variables (i.e., image channels) have a large impact on the overall accuracy of the image classification. High-dimension datasets should be reduced so that only uncorrelated important variables are used in classifications. Despite the fact that RF is an ensemble approach, independent error assessments should be used to evaluate RF results, and iterative classifications are recommended to assess the stability of predicted classes. Results are also shown to be highly sensitive to the size of the training data set. In addition to being as large as possible, the training data sets used in RF classification should also be (a) randomly distributed or created in a manner that allows for the class proportions of the training data to be representative of actual class proportions in the landscape; and (b) should have minimal spatial autocorrelation to improve classification results and to mitigate inflated estimates of RF out-of-bag classification accuracy. Remote Sensing 2015-07-06 7 7 Article 10.3390/rs70708489 8489 8515 2072-4292 2015-07-06 doi: 10.3390/rs70708489 Koreen Millard Murray Richardson http://www.mdpi.com/2072-4292/7/7/8469 Normally, polarimetric SAR classification is a high-dimensional nonlinear mapping problem. In the realm of pattern recognition, sparse representation is a very efficacious and powerful approach. As classical descriptors of polarimetric SAR, covariance and coherency matrices are Hermitian semidefinite and form a Riemannian manifold. Conventional Euclidean metrics are not suitable for a Riemannian manifold, and hence, normal sparse representation classification cannot be applied to polarimetric SAR directly. This paper proposes a new land cover classification approach for polarimetric SAR. There are two principal novelties in this paper. First, a Stein kernel on a Riemannian manifold instead of Euclidean metrics, combined with sparse representation, is employed for polarimetric SAR land cover classification. This approach is named Stein-sparse representation-based classification (SRC). Second, using simultaneous sparse representation and reasonable assumptions of the correlation of representation among different frequency bands, Stein-SRC is generalized to simultaneous Stein-SRC for multi-frequency polarimetric SAR classification. These classifiers are assessed using polarimetric SAR images from the Airborne Synthetic Aperture Radar (AIRSAR) sensor of the Jet Propulsion Laboratory (JPL) and the Electromagnetics Institute Synthetic Aperture Radar (EMISAR) sensor of the Technical University of Denmark (DTU). Experiments on single-band and multi-band data both show that these approaches acquire more accurate classification results in comparison to many conventional and advanced classifiers. Remote Sensing 2015-07-02 7 7 Article 10.3390/rs70708469 8469 8488 2072-4292 2015-07-02 doi: 10.3390/rs70708469 Fan Yang Wei Gao Bin Xu Jian Yang http://www.mdpi.com/2072-4292/7/7/8453 Airborne laser scanning (ALS) is increasingly being used to enhance the accuracy of biomass estimates in tropical forests. Although the technological development of ALS instruments has resulted in ever-greater pulse densities, studies in boreal and sub-boreal forests have shown consistent results even at relatively small pulse densities. The objective of the present study was to assess the effects of reduced pulse density on (1) the digital terrain model (DTM), and (2) canopy metrics derived from ALS data collected in a tropical rainforest in Tanzania. We used a total of 612 coordinates measured with a differential dual frequency Global Navigation Satellite System receiver to analyze the effects on DTMs at pulse densities of 8, 4, 2, 1, 0.5, and 0.025 pulses·m−2. Furthermore, canopy metrics derived for each pulse density and from four different field plot sizes (0.07, 0.14, 0.21, and 0.28 ha) were analyzed. Random variation in DTMs and canopy metrics increased with reduced pulse density. Similarly, increased plot size reduced variation in canopy metrics. A reliability ratio, quantifying replication effects in the canopy metrics, indicated that most of the common metrics assessed were reliable at pulse densities >0.5 pulses·m−2 at a plot size of 0.07 ha. Remote Sensing 2015-06-30 7 7 Article 10.3390/rs70708453 8453 8468 2072-4292 2015-06-30 doi: 10.3390/rs70708453 Endre Hansen Terje Gobakken Erik Næsset http://www.mdpi.com/2072-4292/7/7/8436 Spatially-detailed forest height data are useful to monitor local, regional and global carbon cycle. LiDAR remote sensing can measure three-dimensional forest features but generating spatially-contiguous forest height maps at a large scale (e.g., continental and global) is problematic because existing LiDAR instruments are still data-limited and expensive. This paper proposes a new approach based on an artificial neural network (ANN) for modeling of forest canopy heights over the China continent. Our model ingests spaceborne LiDAR metrics and multiple geospatial predictors including climatic variables (temperature and precipitation), forest type, tree cover percent and land surface reflectance. The spaceborne LiDAR instrument used in the study is the Geoscience Laser Altimeter System (GLAS), which can provide within-footprint forest canopy heights. The ANN was trained with pairs between spatially discrete LiDAR metrics and full gridded geo-predictors. This generates valid conjugations to predict heights over the China continent. The ANN modeled heights were evaluated with three different reference data. First, field measured tree heights from three experiment sites were used to validate the ANN model predictions. The observed tree heights at the site-scale agreed well with the modeled forest heights (R = 0.827, and RMSE = 4.15 m). Second, spatially discrete GLAS observations and a continuous map from the interpolation of GLAS-derived tree heights were separately used to evaluate the ANN model. We obtained R of 0.725 and RMSE of 7.86 m and R of 0.759 and RMSE of 8.85 m, respectively. Further, inter-comparisons were also performed with two existing forest height maps. Our model granted a moderate agreement with the existing satellite-based forest height maps (R = 0.738, and RMSE = 7.65 m (R2 = 0.52, and RMSE = 8.99 m). Our results showed that the ANN model developed in this paper is capable of estimating forest heights over the China continent with a satisfactory accuracy. Forth coming research on our model will focus on extending the model to the estimation of woody biomass. Remote Sensing 2015-06-30 7 7 Article 10.3390/rs70708436 8436 8452 2072-4292 2015-06-30 doi: 10.3390/rs70708436 Xiliang Ni Yuke Zhou Chunxiang Cao Xuejun Wang Yuli Shi Taejin Park Sungho Choi Ranga Myneni http://www.mdpi.com/2072-4292/7/7/8416 Soil lead content is an important parameter in environmental and industrial applications. Chemical analysis, the most commonly method for studying soil samples, are costly, however application of soil spectroscopy presents a more viable alternative. The first step in the method is usually to extract some appropriate spectral features and then regression models are applied to these extracted features. The aim of this paper was to design an accurate and robust regression technique to estimate soil lead contents from laboratory observed spectra. Three appropriate spectral features were selected according to information from other research as well as the spectrum interpretation of field collected soil samples containing lead. These features were then applied to common Multiple Linear Regression (MLR), Partial Least Square Regression (PLSR) and Neural Network (NN) regression models. Results showed that although NN had adequate accuracy, it produced unstable results (i.e., variation of response in different runs). This problem was addressed with application of a Fuzzy Neural Network (FNN) with a least square training strategy. In addition to the stabilized and unique response, the capability of the proposed FNN was proved in terms of regression accuracy where a Ratio of Performance to Deviation (RPD) of 8.76 was achieved for test samples. Remote Sensing 2015-06-29 7 7 Article 10.3390/rs70708416 8416 8435 2072-4292 2015-06-29 doi: 10.3390/rs70708416 Rohollah Goodarzi Mehdi Mokhtarzade M. Zoej http://www.mdpi.com/2072-4292/7/7/8391 Surface reflectance has a central role in the analysis of land surface for a broad variety of Earth System studies. An accurate atmospheric correction, obtained by an appropriate selection of aerosol model, is the first requirement for reliable surface reflectance estimation. In the aerosol model, the type is defined by the physical and chemical properties, while the loading is usually described by the optical thickness at 550 nm. The aim of this work is to evaluate the radiative impact of the aerosol model on the surface reflectance obtained from Compact High Resolution Imaging Spectrometer (CHRIS) hyperspectral data over land by using the specifically developed algorithm CHRIS Atmospherically Corrected Reflectance Imagery (CHRIS@CRI) based on the 6SV radiative transfer model. We employed five different aerosol models: one provided by the AERONET inversion products (used as reference), three standard aerosol models in 6SV, and one obtained from the output of the GEOS-Chem global chemistry-transport model (CTM). The results obtained for the two case studies selected over Benelux show that in the absence of AERONET data on the scene, the best performing aerosol model is the one derived from CTM output. Remote Sensing 2015-06-26 7 7 Article 10.3390/rs70708391 8391 8415 2072-4292 2015-06-26 doi: 10.3390/rs70708391 Cecilia Tirelli Gabriele Curci Ciro Manzo Paolo Tuccella Cristiana Bassani http://www.mdpi.com/2072-4292/7/7/8368 Land cover classification has been widely investigated in remote sensing for agricultural, ecological and hydrological applications. Landsat images with multispectral bands are commonly used to study the numerous classification methods in order to improve the classification accuracy. Thermal remote sensing provides valuable information to investigate the effectiveness of the thermal bands in extracting land cover patterns. k-NN and Random Forest algorithms were applied to both the single Landsat 8 image and the time series Landsat 4/5 images for the Attert catchment in the Grand Duchy of Luxembourg, trained and validated by the ground-truth reference data considering the three level classification scheme from COoRdination of INformation on the Environment (CORINE) using the 10-fold cross validation method. The accuracy assessment showed that compared to the visible and near infrared (VIS/NIR) bands, the time series of thermal images alone can produce comparatively reliable land cover maps with the best overall accuracy of 98.7% to 99.1% for Level 1 classification and 93.9% to 96.3% for the Level 2 classification. In addition, the combination with the thermal band improves the overall accuracy by 5% and 6% for the single Landsat 8 image in Level 2 and Level 3 category and provides the best classified results with all seven bands for the time series of Landsat TM images. Remote Sensing 2015-06-26 7 7 Article 10.3390/rs70708368 8368 8390 2072-4292 2015-06-26 doi: 10.3390/rs70708368 Liya Sun Karsten Schulz http://www.mdpi.com/2072-4292/7/7/8348 Identifying historical forest disturbances is difficult, especially in selectively logged areas. LiDAR is able to measure fine-scale variations in forest structure over multiple kilometers. We use LiDAR data from ca. 16 km2 of forest in Sierra Leone, West Africa, to discriminate areas of old-growth from areas recovering from selective logging for 23 years. We examined canopy height variation and gap size distributions. We found that though recovering blocks of forest differed little in height from old-growth forest (up to 3 m), they had a greater area of canopy gaps (average 10.2% gap fraction in logged areas, compared to 5.6% in unlogged area); and greater numbers of gaps penetrating to the forest floor (162 gaps at 2 m height in logged blocks, and 101 in an unlogged block). Comparison of LiDAR measurements with field data demonstrated that LiDAR delivered accurate results. We found that gap size distributions deviated from power-laws reported previously, with substantially fewer large gaps than predicted by power-law functions. Our analyses demonstrate that LiDAR is a useful tool for distinguishing structural differences between old-growth and old-secondary forests. That makes LiDAR a powerful tool for REDD+ (Reduction of Emissions from Deforestation and Forest Degradation) programs implementation and conservation planning. Remote Sensing 2015-06-26 7 7 Article 10.3390/rs70708348 8348 8367 2072-4292 2015-06-26 doi: 10.3390/rs70708348 Rafi Kent Jeremy Lindsell Gaia Laurin Riccardo Valentini David Coomes http://www.mdpi.com/2072-4292/7/7/8323 Interferometric Synthetic Aperture Radar (InSAR) capability to detect slow deformation over terrain areas is limited by temporal decorrelation, geometric decorrelation and atmospheric artefacts. Multitemporal InSAR methods such as Persistent Scatterer (PS-InSAR) and Small Baseline Subset (SBAS) have been developed to deal with various aspects of decorrelation and atmospheric problems affecting InSAR observations. Nevertheless, the applicability of both PS-InSAR and SBAS in mountainous regions is still challenging. Correct phase unwrapping in both methods is hampered due to geometric decorrelation in particular when using C-band SAR data for deformation analysis. In this paper, we build upon the SBAS method implemented in StaMPS software and improved the technique, here called ISBAS, to assess tectonic and volcanic deformation in the center of the Alborz Mountains in Iran using both Envisat and ALOS SAR data. We modify several aspects within the chain of the processing including: filtering prior to phase unwrapping, topographic correction within three-dimensional phase unwrapping, reducing the atmospheric noise with the help of additional GPS data, and removing the ramp caused by ionosphere turbulence and/or orbit errors to better estimate crustal deformation in this tectonically active region. Topographic correction is done within the three-dimensional unwrapping in order to improve the phase unwrapping process, which is in contrast to previous methods in which DEM error is estimated before/after phase unwrapping. Our experiments show that our improved SBAS approach is able to better characterize the tectonic and volcanic deformation in the center of the Alborz region than the classical SBAS. In particular, Damavand volcano shows an average uplift rate of about 3 mm/year in the year 2003–2010. The Mosha fault illustrates left-lateral motion that could be explained with a fault that is locked up to 17–18 km depths and slips with 2–4 mm/year below that depth. Remote Sensing 2015-06-25 7 7 Article 10.3390/rs70708323 8323 8347 2072-4292 2015-06-25 doi: 10.3390/rs70708323 Sanaz Vajedian Mahdi Motagh Faramarz Nilfouroushan http://www.mdpi.com/2072-4292/7/7/8300 Secondary forest succession on abandoned agricultural land has played a significant role in land cover changes in Europe over the past several decades. However, it is difficult to quantify over large areas. In this paper, we present a conceptual framework for mapping forest succession patterns using vegetation structure information derived from LiDAR data supported by national topographic vector data. This work was performed in the Szczawnica commune in the Polish Carpathians. Using object-based image analysis segments of no vegetation, and sparse/dense low/medium/high vegetation were distinguished and subsequently compared to the national topographic dataset to delineate agricultural land that is covered by vegetation, which indicates secondary succession on abandoned fields. The results showed that 18.7% of the arable land and 40.4% of grasslands, that is 31.0% of the agricultural land in the Szczawnica commune, may currently be experiencing secondary forest succession. The overall accuracy of the approach was assessed using georeferenced terrestrial photographs and was found to be 95.0%. The results of this study indicate that the proposed methodology can potentially be applied in large-scale mapping of secondary forest succession patterns on abandoned land in mountain areas. Remote Sensing 2015-06-25 7 7 Article 10.3390/rs70708300 8300 8322 2072-4292 2015-06-25 doi: 10.3390/rs70708300 Natalia Kolecka Jacek Kozak Dominik Kaim Monika Dobosz Christian Ginzler Achilleas Psomas http://www.mdpi.com/2072-4292/7/7/8271 The growing availability of high-resolution satellite imagery provides an opportunity for identifying road objects. Most studies associated with road detection are scene-related and also based on the digital number of each pixel. Because images can provide more details (including color, size, shape, and texture), object-based processing is more advantageous. Therefore, in this paper, to handle the existing uncertainty of satellite image pixel values, using type-2 fuzzy set theory in combination with object-based image analysis is proposed. Because the main challenges of the type-2 fuzzy set are parameter tuning and extensive computations, a hybrid genetic algorithm (GA) consisting of Pittsburgh and cooperative-competitive learning schemes is proposed to address these problems. The most prominent feature of our research in this work is to establish a comprehensive object-based type-2 fuzzy logic system that enables us to detect roads in high-resolution satellite images with no training data. The validation assessment of road detection results using the proposed framework for independent images demonstrates the capability and efficiency of our method in identifying road objects. For more evaluation, a type-1 fuzzy logic system with the same structure as type-2 is tuned. Evaluations show that type-1 fuzzy logic system quality in training is very similar to that of the proposed type-2 fuzzy framework. However, in general, its lower accuracy, as inferred by validation assessments, makes the type-1 fuzzy logic system significantly different from the proposed type-2. Remote Sensing 2015-06-25 7 7 Article 10.3390/rs70708271 8271 8299 2072-4292 2015-06-25 doi: 10.3390/rs70708271 Maryam Nikfar Mohammad Zoej Mehdi Mokhtarzade Mahdi Shoorehdeli http://www.mdpi.com/2072-4292/7/7/8250 Surface soil moisture (SM) plays a fundamental role in energy and water partitioning in the soil–plant–atmosphere continuum. A reliable and operational algorithm is much needed to retrieve regional surface SM at high spatial and temporal resolutions. Here, we provide an operational framework of estimating surface SM at fine spatial resolutions (using visible/thermal infrared images and concurrent meteorological data) based on a trapezoidal space defined by remotely sensed vegetation cover (Fc) and land surface temperature (LST). Theoretical solutions of the wet and dry edges were derived to achieve a more accurate and effective determination of the Fc/LST space. Subjectivity and uncertainty arising from visual examination of extreme boundaries can consequently be largely reduced. In addition, theoretical derivation of the extreme boundaries allows a per-pixel determination of the VI/LST space such that the assumption of uniform atmospheric forcing over the entire domain is no longer required. The developed approach was tested at the Tibetan Plateau Soil Moisture/Temperature Monitoring Network (SMTMN) site in central Tibet, China, from August 2010 to August 2011 using Moderate Resolution Imaging Spectroradiometer (MODIS) Terra images. Results indicate that the developed trapezoid model reproduced the spatial and temporal patterns of observed surface SM reasonably well, with showing a root-mean-square error of 0.06 m3·m−3 at the site level and 0.03 m3·m−3 at the regional scale. In addition, a case study on 2 September 2010 highlighted the importance of the theoretically calculated wet and dry edges, as they can effectively obviate subjectivity and uncertainties in determining the Fc/LST space arising from visual interpretation of satellite images. Compared with Land Surface Models (LSMs) in Global Land Data Assimilation System-1, the remote sensing-based trapezoid approach gave generally better surface SM estimates, whereas the LSMs showed systematic underestimation. Sensitivity analyses suggested that the trapezoid method is most sensitive to field capacity and temperature but less sensitive to other meteorological observations and parameters. Remote Sensing 2015-06-24 7 7 Article 10.3390/rs70708250 8250 8270 2072-4292 2015-06-24 doi: 10.3390/rs70708250 Yuting Yang Huade Guan Di Long Bing Liu Guanghua Qin Jun Qin Okke Batelaan http://www.mdpi.com/2072-4292/7/6/8224 Impervious surfaces (IS) are a key indicator of environmental quality, and mapping of urban IS is important for a wide range of applications including hydrological modelling, water management, urban and environmental planning and urban climate studies. This paper addresses the accuracy and applicability of vegetation indices (VI), from Landsat imagery, to estimate IS fractions for European cities. The accuracy of three different measures of vegetation cover is examined for eight urban areas at different locations in Europe. The Normalized Difference Vegetation Index (NDVI) and Soil Adjusted Vegetation Index (SAVI) are converted to IS fractions using a regression modelling approach. Also, NDVI is used to estimate fractional vegetation cover (FR), and consequently IS fractions. All three indices provide fairly accurate estimates (MAEs ≈ 10%, MBE's < 2%) of sub-pixel imperviousness, and are found to be applicable for cities with dissimilar climatic and vegetative conditions. The VI/IS relationship across cities is examined by quantifying the MAEs and MBEs between all combinations of models and urban areas. Also, regional regression models are developed by compiling data from multiple cities to examine the potential for developing and applying a single regression model to estimate IS fractions for numerous urban areas without reducing the accuracy considerably. Our findings indicate that the models can be applied broadly for multiple urban areas, and that the accuracy is reduced only marginally by applying the regional models. SAVI is identified as a superior index for the development of regional quantification models. The findings of this study highlight that IS fractions, and spatiotemporal changes herein, can be mapped by use of simple regression models based on VIs from remote sensors, and that the method presented enables simple, accurate and resource efficient quantification of IS. Remote Sensing 2015-06-19 7 6 Article 10.3390/rs70608224 8224 8249 2072-4292 2015-06-19 doi: 10.3390/rs70608224 Per Kaspersen Rasmus Fensholt Martin Drews http://www.mdpi.com/2072-4292/7/6/8202 Land subsidence is a worldwide problem that is typically caused by human activities, primarily the removal of groundwater. In Western Taiwan, groundwater has been pumped for industrial, residential, agricultural, and aquacultural uses for over 40 years. In this study, a multisensor monitoring system comprising GPS stations, leveling surveys, monitoring wells, and Persistent Scatterer Interferometric Synthetic Aperture Radar (PS-InSAR) was employed to monitor land subsidence in Western Taiwan. The results indicate that land subsidence in Yunlin County was mainly affected by the compaction of subsurface soils and over-pumping of groundwater from deep soils. The study area comprised western foothills, characterized by sediments containing predominantly gravel, and coastal areas, where clay was predominant. The subsidence in coastal areas was more severe than that in the western foothills, as a result of groundwater removal. An additional factor affecting subsidence was the compaction of deep layers caused by deep groundwater removal and the deep-layer compaction was difficult to recover. Based on multisensor monitoring results, severe subsidence is mainly affected by compaction of subsurface soils, over-pumping of groundwater from deep soils, and deep soil compaction. Remote Sensing 2015-06-19 7 6 Article 10.3390/rs70608202 8202 8223 2072-4292 2015-06-19 doi: 10.3390/rs70608202 Wei-Chen Hsu Hung-Cheng Chang Kuan-Tsung Chang En-Kai Lin Jin-King Liu Yuei-An Liou http://www.mdpi.com/2072-4292/7/6/8180 Quantifying post-fire effects in a forested landscape is important to ascertain burn severity, ecosystem recovery and post-fire hazard assessments and mitigation planning. Reporting of such post-fire effects assumes significance in fire-prone countries such as USA, Australia, Spain, Greece and Portugal where prescribed burns are routinely carried out. This paper describes the use of Terrestrial Laser Scanning (TLS) to estimate and map change in the forest understorey following a prescribed burn. Eighteen descriptive metrics are derived from bi-temporal TLS which are used to analyse and visualise change in a control and fire-altered plot. Metrics derived are Above Ground Height-based (AGH) percentiles and heights, point count and mean intensity. Metrics such as AGH50change, mean AGHchange and point countchange are sensitive enough to detect subtle fire-induced change (28%–52%) whilst observing little or no change in the control plot (0–4%). A qualitative examination with field measurements of the spatial distribution of burnt areas and percentage area burnt also show similar patterns. This study is novel in that it examines the behaviour of TLS metrics for estimating and mapping fire induced change in understorey structure in a single-scan mode with a minimal fixed reference system. Further, the TLS-derived metrics can be used to produce high resolution maps of change in the understorey landscape. Remote Sensing 2015-06-19 7 6 Article 10.3390/rs70608180 8180 8201 2072-4292 2015-06-19 doi: 10.3390/rs70608180 Vaibhav Gupta Karin Reinke Simon Jones Luke Wallace Lucas Holden http://www.mdpi.com/2072-4292/7/6/8154 The development of low-cost and relatively simple tools to identify emerging land degradation across complex regions is fundamental to plan monitoring and intervention strategies. We propose a procedure that integrates multi-spectral satellite observations and air temperature data to detect areas where the current status of local vegetation and climate shows evident departures from the mean conditions of the investigated region. Our procedure was tested in Basilicata (Italy), which is a typical bio-geographic example of vulnerable Mediterranean landscape. We grouped Landsat TM/ETM+ NDVI and air temperature (T) data by vegetation cover type to estimate the statistical distributions of the departures of NDVI and T from the respective land cover class means. The pixels characterized by contextual left tail NDVI values and right tail T values that persisted in time (2002–2006) were classified as critical to land degradation. According to our results, most of the critical areas (88.6%) corresponded to forests affected by erosion and to riparian buffers that are shaped by fragmentation, as confirmed by aerial and in-situ surveys. Our procedure enables cost-effective screenings of complex areas able to identify raising hotspots that require urgent and deeper investigations. Remote Sensing 2015-06-19 7 6 Article 10.3390/rs70608154 8154 8179 2072-4292 2015-06-19 doi: 10.3390/rs70608154 Maria Lanfredi Rosa Coppola Tiziana Simoniello Rosa Coluzzi Mariagrazia D\'Emilio Vito Imbrenda Maria Macchiato http://www.mdpi.com/2072-4292/7/6/8128 The current study presents an algorithm to retrieve surface Soil Moisture (SM) from multi-temporal Synthetic Aperture Radar (SAR) data. The developed algorithm is based on the Cumulative Density Function (CDF) transformation of multi-temporal RADARSAT-2 backscatter coefficient (BC) to obtain relative SM values, and then converts relative SM values into absolute SM values using soil information. The algorithm is tested in a semi-arid tropical region in South India using 30 satellite images of RADARSAT-2, SMOS L2 SM products, and 1262 SM field measurements in 50 plots spanning over 4 years. The validation with the field data showed the ability of the developed algorithm to retrieve SM with RMSE ranging from 0.02 to 0.06 m3/m3 for the majority of plots. Comparison with the SMOS SM showed a good temporal behaviour with RMSE of approximately 0.05 m3/m3 and a correlation coefficient of approximately 0.9. The developed model is compared and found to be better than the change detection and delta index model. The approach does not require calibration of any parameter to obtain relative SM and hence can easily be extended to any region having time series of SAR data available. Remote Sensing 2015-06-18 7 6 Article 10.3390/rs70608128 8128 8153 2072-4292 2015-06-18 doi: 10.3390/rs70608128 Sat Tomer Ahmad Al Bitar Muddu Sekhar Mehrez Zribi S. Bandyopadhyay K. Sreelash A.K. Sharma Samuel Corgne Yann Kerr http://www.mdpi.com/2072-4292/7/6/8107 Soil management practices that maintain crop residue cover and reduce tillage improve soil structure, increase organic matter content in the soil, positively influence water infiltration, evaporation and soil temperature, and play an important role in fixing CO2 in the soil. Consequently, good residue management practices on agricultural land have many positive impacts on soil quality, crop production quality and decrease the rate of soil erosion. Several studies have been undertaken to develop and test methods to derive information on crop residue cover and soil tillage using empirical and semi-empirical methods in combination with remote sensing data. However, these methods are generally not sufficiently rigorous and accurate for characterizing the spatial variability of crop residue cover in agricultural fields. The goal of this research is to investigate the potential of hyperspectral Hyperion (Earth Observing-1, EO-1) data and constrained linear spectral mixture analysis (CLSMA) for percent crop residue cover estimation and mapping. Hyperion data were acquired together with ground-reference measurements for validation purposes at the beginning of the agricultural season (prior to spring crop planting) in Saskatchewan (Canada). At this time, only bare soil and crop residue were present with no crop cover development. In order to extract the crop residue fraction, the images were preprocessed, and then unmixed considering the entire spectral range (427 nm–2355 nm) and the pure spectra (endmember). The results showed that the correlation between ground-reference measurements and extracted fractions from the Hyperion data using CLMSA showed that the model was overall a very good predictor for crop residue percent cover (index of agreement (D) of 0.94, coefficient of determination (R2) of 0.73 and root mean square error (RMSE) of 8.7%) and soil percent cover (D of 0.91, R2 of 0.68 and RMSE of 10.3%). This performance of Hyperion is mainly due to the spectral band characteristics, especially the availability of contiguous narrow bands in the short-wave infrared (SWIR) region, which is sensitive to the residue (lignin and cellulose absorption features). Remote Sensing 2015-06-18 7 6 Article 10.3390/rs70608107 8107 8127 2072-4292 2015-06-18 doi: 10.3390/rs70608107 Abderrazak Bannari Karl Staenz Catherine Champagne K. Khurshid http://www.mdpi.com/2072-4292/7/6/8102 This Editorial introduces the papers published in the special issue "Earth Observation for Ecosystems Monitoring in Space and Time" which includes the most important researchers in the field and the most challenging aspects of the application of remote sensing to study ecosystems. Remote Sensing 2015-06-18 7 6 Editorial 10.3390/rs70608102 8102 8106 2072-4292 2015-06-18 doi: 10.3390/rs70608102 Duccio Rocchini http://www.mdpi.com/2072-4292/7/6/8067 Solar surface radiation data of high quality is essential for the appropriate monitoring and analysis of the Earth's radiation budget and the climate system. Further, they are crucial for the efficient planning and operation of solar energy systems. However, well maintained surface measurements are rare in many regions of the world and over the oceans. There, satellite derived information is the exclusive observational source. This emphasizes the important role of satellite based surface radiation data. Within this scope, the new satellite based CM-SAF SARAH (Solar surfAce RAdiation Heliosat) data record is discussed as well as the retrieval method used. The SARAH data are retrieved with the sophisticated SPECMAGIC method, which is based on radiative transfer modeling. The resulting climate data of solar surface irradiance, direct irradiance (horizontal and direct normal) and clear sky irradiance are covering 3 decades. The SARAH data set is validated with surface measurements of the Baseline Surface Radiation Network (BSRN) and of the Global Energy and Balance Archive (GEBA). Comparison with BSRN data is performed in order to estimate the accuracy and precision of the monthly and daily means of solar surface irradiance. The SARAH solar surface irradiance shows a bias of 1.3 \(W/m^2\) and a mean absolute bias (MAB) of 5.5 \(W/m^2\) for monthly means. For direct irradiance the bias and MAB is 1 \(W/m^2\) and 8.2 \(W/m^2\) respectively. Thus, the uncertainty of the SARAH data is in the range of the uncertainty of ground based measurements. In order to evaluate the uncertainty of SARAH based trend analysis the time series of SARAH monthly means are compared to GEBA. It has been found that SARAH enables the analysis of trends with an uncertainty of 1 \(W/m^2/dec\); a remarkable good result for a satellite based climate data record. SARAH has been also compared to its legacy version, the satellite based CM-SAF MVIRI climate data record. Overall, SARAH shows a significant higher accuracy and homogeneity than its legacy version. With its high accuracy and temporal and spatial resolution SARAH is well suited for regional climate monitoring and analysis as well as for solar energy applications. Remote Sensing 2015-06-18 7 6 Article 10.3390/rs70608067 8067 8101 2072-4292 2015-06-18 doi: 10.3390/rs70608067 Richard Müller Uwe Pfeifroth Christine Träger-Chatterjee Jörg Trentmann Roswitha Cremer http://www.mdpi.com/2072-4292/7/6/8045 The main objective of the present study was to apply a slope-based spectral method to both dry and fresh pasture vegetation. Differences in eight spectral ranges were identified across the near infrared-shortwave infrared (NIR-SWIR) that were indicative of changes in chemical properties. Slopes across these ranges were calculated and a partial least squares (PLS) analytical model was constructed for the slopes vs. crude protein (CP) and neutral detergent fiber (NDF) contents. Different datasets with different numbers of fresh/dry samples were constructed to predict CP and NDF contents. When using a mixed-sample dataset with dry-to-fresh ratios of 85%:15% and 75%:25%, the correlations of CP (R2 = 0.95, in both) and NDF (R2 = 0.84 and 0.82, respectively) were almost as high as when using only dry samples (0.97 and 0.85, respectively). Furthermore, satisfactory correlations were obtained with a dry-to-fresh ratio of 50%:50% for CP (R2 = 0.92). The results of our study are especially encouraging because CP and NDF contents could be predicted even though some of the selected spectral regions were directly affected by atmospheric water vapor or water in the plants. Remote Sensing 2015-06-18 7 6 Article 10.3390/rs70608045 8045 8066 2072-4292 2015-06-18 doi: 10.3390/rs70608045 Rachel Lugassi Alexandra Chudnovsky Eli Zaady Levana Dvash Naftaly Goldshleger http://www.mdpi.com/2072-4292/7/6/8019 Topography affects the fraction of direct and diffuse radiation received on a pixel and changes the sun–target–sensor geometry, resulting in variations in the observed radiance. Retrieval of surface–atmosphere properties from top of atmosphere radiance may need to account for topographic effects. This study investigates how such effects can be taken into account for top of atmosphere radiance modeling. In this paper, a system for top of atmosphere radiance modeling over heterogeneous non-Lambertian rugged terrain through radiative transfer modeling is presented. The paper proposes an extension of "the four-stream radiative transfer theory" (Verhoef and Bach 2003, 2007 and 2012) mainly aimed at representing topography-induced contributions to the top of atmosphere radiance modeling. A detailed account for BRDF effects, adjacency effects and topography effects on the radiance modeling is given, in which sky-view factor and non-Lambertian reflected radiance from adjacent slopes are modeled precisely. The paper also provides a new formulation to derive the atmospheric coefficients from MODTRAN with only two model runs, to make it more computationally efficient and also avoiding the use of zero surface albedo as used in the four-stream radiative transfer theory. The modeling begins with four surface reflectance factors calculated by the Soil–Leaf–Canopy radiative transfer model SLC at the top of canopy and propagates them through the effects of the atmosphere, which is explained by six atmospheric coefficients, derived from MODTRAN radiative transfer code. The top of the atmosphere radiance is then convolved with the sensor characteristics to generate sensor-like radiance. Using a composite dataset, it has been shown that neglecting sky view factor and/or terrain reflected radiance can cause uncertainty in the forward TOA radiance modeling up to 5 (mW/m2·sr·nm). It has also been shown that this level of uncertainty can be translated into an over/underestimation of more than 0.5 in LAI (or 0.07 in fCover) in variable retrieval. Remote Sensing 2015-06-18 7 6 Article 10.3390/rs70608019 8019 8044 2072-4292 2015-06-18 doi: 10.3390/rs70608019 Alijafar Mousivand Wout Verhoef Massimo Menenti Ben Gorte http://www.mdpi.com/2072-4292/7/6/7995 Leaf biomass distribution is a key factor for modeling energy and carbon fluxes in forest canopies and for assessing fire behavior. We propose a new method to estimate 3D leaf bulk density distribution, based on a calibration of indices derived from T-LiDAR. We applied the method to four contrasted plots in a mature Quercus pubescens forest. Leaf bulk densities were measured inside 0.7 m-diameter spheres, referred to as Calibration Volumes. Indices were derived from LiDAR point clouds and calibrated over the Calibration Volume bulk densities. Several indices were proposed and tested to account for noise resulting from mixed pixels and other theoretical considerations. The best index and its calibration parameter were then used to estimate leaf bulk densities at the grid nodes of each plot. These LiDAR-derived bulk density distributions were used to estimate bulk density vertical profiles and loads and above four meters compared well with those assessed by the classical inventory-based approach. Below four meters, the LiDAR-based approach overestimated bulk densities since no distinction was made between wood and leaf returns. The results of our method are promising since they demonstrate the possibility to assess bulk density on small plots at a reasonable operational cost. Remote Sensing 2015-06-18 7 6 Article 10.3390/rs70607995 7995 8018 2072-4292 2015-06-18 doi: 10.3390/rs70607995 François Pimont Jean-Luc Dupuy Eric Rigolot Vincent Prat Alexandre Piboule http://www.mdpi.com/2072-4292/7/6/7987 Habitat quality is the ability of the environment to provide conditions appropriate for individual and species persistence. Measuring or monitoring habitat quality requires complex integration of many properties of the ecosystem, where traditional terrestrial data collection methods have proven extremely time-demanding. Remote sensing has known potential to map various ecosystem properties, also allowing rigorous checking of accuracy and supporting standardized processing. Our Special Issue presents examples where remote sensing has been successfully used for habitat mapping, quantification of habitat quality parameters, or multi-parameter modelling of habitat quality itself. New frontiers such as bathymetric scanning, grassland vegetation classification and operational use were explored, various new ecological verification methods were introduced and integration with ongoing habitat conservation schemes was demonstrated. These studies show that remote sensing and Geoinformation Science for habitat quality analysis have evolved from isolated experimental studies to an active field of research with a dedicated community. It is expected that these new methods will substantially contribute to biodiversity conservation worldwide. Remote Sensing 2015-06-17 7 6 Editorial 10.3390/rs70607987 7987 7994 2072-4292 2015-06-17 doi: 10.3390/rs70607987 András Zlinszky Hermann Heilmeier Heiko Balzter Bálint Czúcz Norbert Pfeifer http://www.mdpi.com/2072-4292/7/6/7959 Timely and accurate information on the global cropland extent is critical for applications in the fields of food security, agricultural monitoring, water management, land-use change modeling and Earth system modeling. On the one hand, it gives detailed location information on where to analyze satellite image time series to assess crop condition. On the other hand, it isolates the agriculture component to focus food security monitoring on agriculture and to assess the potential impacts of climate change on agricultural lands. The cropland class is often poorly captured in global land cover products due to its dynamic nature and the large variety of agro-systems. The overall objective was to evaluate the current availability of cropland datasets in order to propose a strategic planning and effort distribution for future cropland mapping activities and, therefore, to maximize their impact. Following a very comprehensive identification and collection of national to global land cover maps, a multi-criteria analysis was designed at the country level to identify the priority areas for cropland mapping. As a result, the analysis highlighted priority regions, such as Western Africa, Ethiopia, Madagascar and Southeast Asia, for the remote sensing community to focus its efforts. A Unified Cropland Layer at 250 m for the year 2014 was produced combining the fittest products. It was assessed using global validation datasets and yields an overall accuracy ranging from 82%–94%. Masking cropland areas with a global forest map reduced the commission errors from 46% down to 26%. Compared to the GLC-Share and the International Institute for Applied Systems Analysis-International Food Policy Research Institute (IIASA-IFPRI) cropland maps, significant spatial disagreements were found, which might be attributed to discrepancies in the cropland definition. This advocates for a shared definition of cropland, as well as global validation datasets relevant for the agriculture class in order to systematically assess existing and future cropland maps. Remote Sensing 2015-06-17 7 6 Article 10.3390/rs70607959 7959 7986 2072-4292 2015-06-17 doi: 10.3390/rs70607959 François Waldner Steffen Fritz Antonio Di Gregorio Pierre Defourny http://www.mdpi.com/2072-4292/7/6/7938 The hydrodynamics of low-lying riverine floodplains and wetlands play a critical role in hydrology and ecosystem processes. Because small topographic features affect floodplain storage and flow velocity, a hydrodynamic model setup of these regions imposes more stringent requirements on the input Digital Elevation Model (DEM) compared to upland regions with comparatively high slopes. This current study provides a systematic approach to evaluate the required relative vertical accuracy and spatial resolution of current and future satellite-based altimeters within the context of DEM requirements for 2-D floodplain hydrodynamic models. A case study is presented for the Atchafalaya Basin with a model domain of 1190 km2. The approach analyzes the sensitivity of modeled floodplain water elevation and velocity to typical satellite-based DEM grid-box scale and vertical error, using a previously calibrated version of the physically-based flood inundation model (LISFLOOD-ACC). Results indicate a trade-off relationship between DEM relative vertical error and grid-box size. Higher resolution models are the most sensitive to vertical accuracy, but the impact diminishes at coarser resolutions because of spatial averaging. The results provide guidance to engineers and scientists when defining the observation scales of future altimetry missions such as the Surface Water and Ocean Topography (SWOT) mission from the perspective of numerical modeling requirements for large floodplains of O[103] km2 and greater. Remote Sensing 2015-06-17 7 6 Article 10.3390/rs70607938 7938 7958 2072-4292 2015-06-17 doi: 10.3390/rs70607938 Hahn Jung Michael Jasinski http://www.mdpi.com/2072-4292/7/6/7914 Satellite remote sensing offers an effective approach to estimate indicators of air quality on a large scale. It is critically significant for air quality monitoring in areas experiencing rapid urbanization and consequently severe air pollution, like the Pearl River Delta (PRD) in China. This paper starts with examining ground observations of particulate matter (PM) and the relationship between PM10 (particles smaller than 10 μm) and aerosol optical thickness (AOT) by analyzing observations on the sampling sites in the PRD. A linear regression (R2 = 0.51) is carried out using MODIS-derived 500 m-resolution AOT and PM10 concentration from monitoring stations. Data of atmospheric boundary layer (ABL) height and relative humidity are used to make vertical and humidity corrections on AOT. Results after correction show higher correlations (R2 = 0.55) between extinction coefficient and PM10. However, coarse spatial resolution of meteorological data affects the smoothness of retrieved maps, which suggests high-resolution and accurate meteorological data are critical to increase retrieval accuracy of PM. Finally, the model provides the spatial distribution maps of instantaneous and yearly average PM10 over the PRD. It is proved that observed PM10 is more relevant to yearly mean AOT than instantaneous values. Remote Sensing 2015-06-17 7 6 Article 10.3390/rs70607914 7914 7937 2072-4292 2015-06-17 doi: 10.3390/rs70607914 Lili Li Jingxue Yang Yunpeng Wang http://www.mdpi.com/2072-4292/7/6/7892 The objective of this study is to develop new algorithms for automated urban forest inventory at the individual tree level using LiDAR point cloud data. LiDAR data contain three-dimensional structure information that can be used to estimate tree height, base height, crown depth, and crown diameter. This allows precision urban forest inventory down to individual trees. Unlike most of the published algorithms that detect individual trees from a LiDAR-derived raster surface, we worked directly with the LiDAR point cloud data to separate individual trees and estimate tree metrics. Testing results in typical urban forests are encouraging. Future works will be oriented to synergize LiDAR data and optical imagery for urban tree characterization through data fusion techniques. Remote Sensing 2015-06-16 7 6 Article 10.3390/rs70607892 7892 7913 2072-4292 2015-06-16 doi: 10.3390/rs70607892 Caiyun Zhang Yuhong Zhou Fang Qiu http://www.mdpi.com/2072-4292/7/6/7865 Due to technical limitations, it is impossible to have high resolution in both spatial and temporal dimensions for current NDVI datasets. Therefore, several methods are developed to produce high resolution (spatial and temporal) NDVI time-series datasets, which face some limitations including high computation loads and unreasonable assumptions. In this study, an unmixing-based method, NDVI Linear Mixing Growth Model (NDVI-LMGM), is proposed to achieve the goal of accurately and efficiently blending MODIS NDVI time-series data and multi-temporal Landsat TM/ETM+ images. This method firstly unmixes the NDVI temporal changes in MODIS time-series to different land cover types and then uses unmixed NDVI temporal changes to predict Landsat-like NDVI dataset. The test over a forest site shows high accuracy (average difference: −0.0070; average absolute difference: 0.0228; and average absolute relative difference: 4.02%) and computation efficiency of NDVI-LMGM (31 seconds using a personal computer). Experiments over more complex landscape and long-term time-series demonstrated that NDVI-LMGM performs well in each stage of vegetation growing season and is robust in regions with contrasting spatial and spatial variations. Comparisons between NDVI-LMGM and current methods (i.e., Spatial and Temporal Adaptive Reflectance Fusion Model (STARFM), Enhanced STARFM (ESTARFM) and Weighted Linear Model (WLM)) show that NDVI-LMGM is more accurate and efficient than current methods. The proposed method will benefit land surface process research, which requires a dense NDVI time-series dataset with high spatial resolution. Remote Sensing 2015-06-16 7 6 Article 10.3390/rs70607865 7865 7891 2072-4292 2015-06-16 doi: 10.3390/rs70607865 Yuhan Rao Xiaolin Zhu Jin Chen Jianmin Wang http://www.mdpi.com/2072-4292/7/6/7846 Validation of land cover products is a fundamental task prior to data applications. Current validation schemes and methods are, however, suited only for assessing classification accuracy and disregard the reliability of land cover products. The reliability evaluation of land cover products should be undertaken to provide reliable land cover information. In addition, the lack of high-quality reference data often constrains validation and affects the reliability results of land cover products. This study proposes a validation schema to evaluate the reliability of land cover products, including two methods, namely, result reliability evaluation and process reliability evaluation. Result reliability evaluation computes the reliability of land cover products using seven reliability indicators. Process reliability evaluation analyzes the reliability propagation in the data production process to obtain the reliability of land cover products. Fuzzy fault tree analysis is introduced and improved in the reliability analysis of a data production process. Research results show that the proposed reliability evaluation scheme is reasonable and can be applied to validate land cover products. Through the analysis of the seven indicators of result reliability evaluation, more information on land cover can be obtained for strategic decision-making and planning, compared with traditional accuracy assessment methods. Process reliability evaluation without the need for reference data can facilitate the validation and reflect the change trends of reliabilities to some extent. Remote Sensing 2015-06-16 7 6 Article 10.3390/rs70607846 7846 7864 2072-4292 2015-06-16 doi: 10.3390/rs70607846 Wenzhong Shi Xiaokang Zhang Ming Hao Pan Shao Liping Cai Xuzhe Lyu http://www.mdpi.com/2072-4292/7/6/7826 Bidirectional reflectance distribution function (BRDF) archetypes extracted from the Moderate Resolution Imaging Spectroradiometer (MODIS) BRDF/Albedo product over the global Earth Observing System Land Validation Core Sites can be used to simplify BRDF models. The present study attempts to evaluate the representativeness of BRDF archetypes for surface reflectance anisotropy. Five-year forward-modeled MODIS multi-angular reflectance (MCD-ref) and aditional actual MODIS multi-angular observations (MCD-obs) in four growing periods in 2008 over three tiles were taken as validation data. First, BRDF archetypes in the principal plane were qualitatively compared with the time-series MODIS BRDF product of randomly sampled pixels. Secondly, BRDF archetypes were used to fit MCD-ref, and the average root-mean-squared errors (RMSEs) over each tile were examined for these five years. Finally, both BRDF archetypes and the MODIS BRDF were used to fit MCD-obs, and the histograms of the fit-RMSEs were compared. The consistency of the directional reflectance between the BRDF archetypes and MODIS BRDFs in nadir-view, hotspot and entire viewing hemisphere at 30° and 50° solar geometries were also examined. The results confirm that BRDF archetypes are representative of surface reflectance anisotropy for available snow-free MODIS data. Remote Sensing 2015-06-15 7 6 Article 10.3390/rs70607826 7826 7845 2072-4292 2015-06-15 doi: 10.3390/rs70607826 Hu Zhang Ziti Jiao Yadong Dong Xiaowen Li http://www.mdpi.com/2072-4292/7/6/7809 Forest dynamics influence climate, biodiversity, and livelihoods at multiple scales, yet current resource policy addressing these dynamics is ineffective without reliable land use land cover change data. The collective impact of harvest decisions by many small forest owners can be substantial at the landscape scale, yet monitoring harvests and regrowth in these forests is challenging. Remote sensing is an obvious route to detect and monitor small-scale land use dynamics over large areas. Using an annual series of Landsat-5 Thematic Mapper (TM) images and a GIS shapefile of property boundaries, we identified units where harvests occurred from 2005 to 2011 using an Object-Based Change Detection (OBCD) approach. Percent of basal area harvested was verified using stand-level harvest data. Our method detected all harvests above 20% basal area removal in all forest types (northern hardwoods, mixed deciduous/coniferous, coniferous), on properties as small as 10 acres (0.4 ha; approximately four Landsat pixels). Our results had a resolution of about 10% basal area (that is, a selective harvest removal of 30% could be distinguished from one of 40%). Our method can be automated and used to measure annual harvest rates and intensities for large areas of the United States, providing critical information on land use transition. Remote Sensing 2015-06-15 7 6 Technical Note 10.3390/rs70607809 7809 7825 2072-4292 2015-06-15 doi: 10.3390/rs70607809 Riccardo Tortini Audrey Mayer Pieralberto Maianti http://www.mdpi.com/2072-4292/7/6/7785 This paper describes a novel strategy for the visualization of hyperspectral imagery based on the analysis of image pixel pairwise distances. The goal of this approach is to generate a final color image with excellent interpretability and high contrast at the cost of distorting a few pairwise distances. Specifically, the principle of equal variance is introduced to divide all hyperspectral bands into three subgroups and to ensure the energy is distributed uniformly between them, as in natural color images. Then, after detecting both normal and outlier pixels, these three subgroups are mapped into three color components of the output visualization using two different mapping (i.e., dimensionality reduction) schemes for the two types of pixels. The widely-used multidimensional scaling (MDS) is used for normal pixels and a new objective function, taking into account the weighting of pairwise distances, is presented for the outlier pixels. The pairwise distance weighting is designed such that small pairwise distances between the outliers and their respective neighbors are emphasized and large deviations are suppressed. This produces an image with high contrast and good interpretability while retaining the detailed information content. The proposed algorithm is compared with several state-of-the-art visualization techniques and evaluated on the well-known AVIRIS hyperspectral images. The effectiveness of the proposed strategy is substantiated both visually and quantitatively. Remote Sensing 2015-06-12 7 6 Article 10.3390/rs70607785 7785 7808 2072-4292 2015-06-12 doi: 10.3390/rs70607785 Yi Long Heng-Chao Li Turgay Celik Nathan Longbotham William Emery http://www.mdpi.com/2072-4292/7/6/7768 Cloud fraction (CF) is known as the dominant modulator of Earth's radiative fluxes. Ground-based CF observations are useful to characterize the cloudiness of a specific site and are valuable for comparison with satellite observations and numerical models. We present for the first time CF statistics (relative to liquid clouds only) for Bern, Switzerland, derived from the observations of a ground-based microwave radiometer. CF is derived with a new method involving the analysis of the integrated liquid water distribution measured by the radiometer. The 10-year analyzed period (2004–2013) allowed us to compute a CF climatology for Bern, showing a maximum CF of 60.9% in winter and a minimum CF of 42.0% in summer. The CF monthly anomalies are identified with respect to the climatological mean values, and they are confirmed through MeteoSwiss yearly climatological bulletins. The CF monthly mean variations are similar to the observations taken at another Swiss location, Payerne, suggesting a large-scale correlation between different sites on the Swiss Plateau. A CF diurnal cycle is also computed, and large intraseasonal variations are found. The overall mean CF diurnal cycle, however, shows a typical sinusoidal cycle, with higher values in the morning and lower values in the afternoon. Remote Sensing 2015-06-11 7 6 Article 10.3390/rs70607768 7768 7784 2072-4292 2015-06-11 doi: 10.3390/rs70607768 Federico Cossu Klemens Hocke Christian Mätzler http://www.mdpi.com/2072-4292/7/6/7753 In this paper, a novel current inversion algorithm from X-band marine radar images is proposed. The routine, for which deep water is assumed, begins with 3-D FFT of the radar image sequence, followed by the extraction of the dispersion shell from the 3-D image spectrum. Next, the dispersion shell is converted to a polar current shell (PCS) using a polar coordinate transformation. After removing outliers along each radial direction of the PCS, a robust sinusoidal curve fitting is applied to the data points along each circumferential direction of the PCS. The angle corresponding to the maximum of the estimated sinusoid function is determined to be the current direction, and the amplitude of this sinusoidal function is the current speed. For validation, the algorithm is tested against both simulated radar images and field data collected by a vertically-polarized X-band system and ground-truthed with measurements from an acoustic Doppler current profiler (ADCP). From the field data, it is observed that when the current speed is less than 0.5 m/s, the root mean square differences between the radar-derived and the ADCP-measured current speed and direction are 7.3 cm/s and 32.7°, respectively. The results indicate that the proposed procedure, unlike most existing current inversion schemes, is not susceptible to high current speeds and circumvents the need to consider aliasing. Meanwhile, the relatively low computational cost makes it an excellent choice in practical marine applications. Remote Sensing 2015-06-11 7 6 Article 10.3390/rs70607753 7753 7767 2072-4292 2015-06-11 doi: 10.3390/rs70607753 Chengxi Shen Weimin Huang Eric Gill Ruben Carrasco Jochen Horstmann http://www.mdpi.com/2072-4292/7/6/7732 In this study, an analysis of multi-temporal and multi-frequency Synthetic Aperture Radar data is performed to investigate the backscatter behavior of various semantic classes in the context of flood mapping in central Europe. The focus is mainly on partially submerged vegetation such as forests and agricultural fields. The test area is located at River Saale, Saxony-Anhalt, Germany, which is covered by a time series of 39 TerraSAR-X data acquired within the time interval December 2009 to June 2013. The data set is supplemented by ALOS PALSAR L-band and RADARSAT-2 C-band data. The time series covers two inundations in January 2011 and June 2013 which allows evaluating backscatter variations between flood periods and normal water level conditions using different radar wavelengths. According to the results, there is potential in detecting flooding beneath vegetation in all microwave wavelengths, even in X-band for sparse vegetation or leaf-off forests. Remote Sensing 2015-06-11 7 6 Article 10.3390/rs70607732 7732 7752 2072-4292 2015-06-11 doi: 10.3390/rs70607732 Sandro Martinis Christoph Rieke http://www.mdpi.com/2072-4292/7/6/7712 Management strategies and silvicultural treatments of fire-prone ecosystems often rely on knowledge of the regeneration potential and long-term recovery ability of vegetation types. Remote sensing and GIS applications are valuable tools providing cost-efficient information on vegetation recovery patterns and their associated environmental factors. In this study we used an ordinal classification scheme to describe the land cover changes induced by a wildfire that occurred in 1983 in Pinus brutia woodlands on Karpathos Aegean Island, south-eastern Greece. As a proxy variable that indicates ecosystem recovery, we also estimated the difference between the NDVI and NBR indices a few months (1984) and almost 30 years after the fire (2012). Environmental explanatory variables were selected using a digital elevation model and various thematic maps. To identify the most influential environmental factors contributing to woodland recovery, binary logistic regression and linear regression techniques were applied. The analyses showed that although a large proportion of the P. brutia woodland has recovered 26 years after the fire event, a considerable amount of woodland had turned into scrub vegetation. Altitude, slope inclination, solar radiation, and pre-fire woodland physiognomy were identified as dominant factors influencing the vegetation's recovery probability. Additionally, altitude and inclination are the variables that explain changes in the satellite remote sensing vegetation indices reflecting the recovery potential. Pinus brutia showed a good post-fire recovery potential, especially in parts of the study area with increased moisture availability. Remote Sensing 2015-06-10 7 6 Article 10.3390/rs70607712 7712 7731 2072-4292 2015-06-10 doi: 10.3390/rs70607712 Foula Nioti Fotios Xystrakis Nikos Koutsias Panayotis Dimopoulos http://www.mdpi.com/2072-4292/7/6/7695 The robust detection of ships is one of the key techniques in coastal and marine applications of synthetic aperture radar (SAR). Conventional SAR ship detectors involved multiple parameters, which need to be estimated or determined very carefully. In this paper, we propose a new ship detection approach based on multi-scale heterogeneities under the a contrario decision framework, with a few parameters that can be easily determined. First, multi-scale heterogeneity features are extracted and fused to build a heterogeneity map, in which ships are well highlighted from backgrounds. Second, a set of reference objects are automatically selected by analyzing the saliency of local regions in the heterogeneity map and then are used to construct a null hypothesis model for the final decision. Finally, the detection results are obtained by using an a contrario decision. Experimental results on real SAR images demonstrate that the proposed method not only works more stably for ships with different sizes, but also has better performance than conventional ship detectors. Remote Sensing 2015-06-10 7 6 Article 10.3390/rs70607695 7695 7711 2072-4292 2015-06-10 doi: 10.3390/rs70607695 Xiaojing Huang Wen Yang Haijian Zhang Gui-Song Xia http://www.mdpi.com/2072-4292/7/6/7671 Stable night-time light data from the Defense Meteorological Satellite Program (DMSP) Operational Line-scan System (OLS) provide a unique proxy for anthropogenic development. This paper presents a regional urban extent extraction method using a one-class classifier and combinations of DMSP/OLS stable night-time light (NTL) data, MODIS normalized difference vegetation index (NDVI) data, and land surface temperature (LST) data. We first analyzed how well MODIS NDVI and LST data quantify the properties of urban areas. Considering that urban area is the only class of interest, we applied the one-class support vector machine (OCSVM) to classify different combinations of the three datasets. We evaluated the effectiveness of the proposed method and compared with the locally optimized threshold method in regional urban extent mapping in China. The experimental results demonstrate that DMSP/OLS NTL data, MODIS NDVI and LST data provide different but complementary information sources to quantify the urban extent at a regional scale. The results also indicate that the OCSVM classification of the combination of all three datasets generally outperformed the locally optimized threshold method. The proposed method effectively and efficiently extracted the urban extent at a regional scale, and is applicable to other study areas. Remote Sensing 2015-06-09 7 6 Technical Note 10.3390/rs70607671 7671 7694 2072-4292 2015-06-09 doi: 10.3390/rs70607671 Xiya Zhang Peijun Li Cai Cai http://www.mdpi.com/2072-4292/7/6/7646 The provenance of observations from a Sensor Web enabled remote sensing application represents a great challenge. There are currently no representations or tracking methods. We propose a provenance method that represents and tracks remote sensing observations in the Sensor Web enabled environment. The representation can be divided into the description model, encoding method, and service implementation. The description model uses a tuple to define four objects (sensor, data, processing, and service) and their relationships at a time point or interval. The encoding method incorporates the description into the Observations & Measurements specification of the Sensor Web. The service implementation addresses the effects of the encoding method on the implementation of Sensor Web services. The tracking method abstracts a common provenance algorithm and four algorithms that track the four objects (sensor, data, processing, and service) in a remote sensing observation application based on the representation. We conducted an experiment on the representation and tracking of provenance information for vegetation condition products, such as the Normalized Difference Vegetation Index (NDVI) and the Vegetation Condition Index (VCI). Our experiments used raw Moderate Resolution Imaging Spectroradiometer (MODIS) data to produce daily NDVI, weekly NDVI, and weekly VCI for the 48 contiguous states of the United States, for May from 2000 to 2012. We also implemented inverse tracking. We evaluated the time and space requirements of the proposed method in this scenario. Our results show that this technique provides a solution for determining provenance information in remote sensing observations. Remote Sensing 2015-06-09 7 6 Article 10.3390/rs70607646 7646 7670 2072-4292 2015-06-09 doi: 10.3390/rs70607646 Zeqiang Chen Nengcheng Chen http://www.mdpi.com/2072-4292/7/6/7615 Wetlands are an important natural resource that requires monitoring. A key step in environmental monitoring is to map the locations and characteristics of the resource to better enable assessment of change over time. Synthetic Aperture Radar (SAR) systems are helpful in this way for wetland resources because their data can be used to map and monitor changes in surface water extent, saturated soils, flooded vegetation, and changes in wetland vegetation cover. We review a few techniques to demonstrate SAR capabilities for wetland monitoring, including the commonly used method of grey-level thresholding for mapping surface water and highlighting changes in extent, and approaches for polarimetric decompositions to map flooded vegetation and changes from one class of land cover to another. We use the Curvelet-based change detection and the Wishart-Chernoff Distance approaches to show how they substantially improve mapping of flooded vegetation and flagging areas of change, respectively. We recommend that the increasing availability SAR data and the proven ability of these data to map various components of wetlands mean SAR should be considered as a critical component of a wetland monitoring system. Remote Sensing 2015-06-09 7 6 Review 10.3390/rs70607615 7615 7645 2072-4292 2015-06-09 doi: 10.3390/rs70607615 Lori White Brian Brisco Mohammed Dabboor Andreas Schmitt Andrew Pratt http://www.mdpi.com/2072-4292/7/6/7597 Understanding the spatial and temporal dynamics of vegetation is essential in drylands. In this paper, we evaluated three vegetation indices, namely the Normalized Difference Vegetation Index (NDVI), the Soil-Adjusted Vegetation Index (SAVI) and the Enhanced Vegetation Index (EVI), derived from the Moderate Resolution Imaging Spectroradiometer (MODIS) Surface-Reflectance Product in the Xinjiang Uygur Autonomous Region, China (XUAR), to assess index time series' suitability for monitoring vegetation dynamics in a dryland environment. The mean annual VI and its variability were generated and analyzed from the three VI time series for the period 2001–2012 across XUAR. Two phenological metrics, start of the season (SOS) and end of the season (EOS), were detected and compared for each vegetation type. The mean annual VI images showed similar spatial patterns of vegetation conditions with varying magnitudes. The EVI exhibited high uncertainties in sparsely vegetated lands and forests. The phenological metrics derived from the three VIs are consistent for most vegetation types, with SOS and EOS generated from NDVI showing the largest deviation. Remote Sensing 2015-06-09 7 6 Article 10.3390/rs70607597 7597 7614 2072-4292 2015-06-09 doi: 10.3390/rs70607597 Linlin Lu Claudia Kuenzer Cuizhen Wang Huadong Guo Qingting Li http://www.mdpi.com/2072-4292/7/6/7571 Changes in soil moisture between two radar acquisitions can impact the observed coherence in differential interferometry: both coherence magnitude |Υ| and phase Φ are affected. The influence on the latter potentially biases the estimation of deformations. These effects have been found to be variable in magnitude and sign, as well as dependent on polarization, as opposed to predictions by existing models. Such diversity can be explained when the soil is modelled as a half-space with spatially varying dielectric properties and a rough interface. The first-order perturbative solution achieves–upon calibration with airborne L band data–median correlations ρ at HH polarization of 0.77 for the phase Φ, of 0.50 for |Υ|, and for the phase triplets ≡ of 0.56. The predictions are sensitive to the choice of dielectric mixing model, in particular the absorptive properties; the differences between the mixing models are found to be partially compensatable by varying the relative importance of surface and volume scattering. However, for half of the agricultural fields the Hallikainen mixing model cannot reproduce the observed sensitivities of the phase to soil moisture. In addition, the first-order expansion does not predict any impact on the HV coherence, which is however empirically found to display similar sensitivities to soil moisture as the co-pol channels HH and VV. These results indicate that the first-order solution, while not able to reproduce all observed phenomena, can capture some of the more salient patterns of the effect of soil moisture changes on the HH and VV DInSAR signals. Hence it may prove useful in separating the deformations from the moisture signals, thus yielding improved displacement estimates or new ways for inferring soil moisture. Remote Sensing 2015-06-08 7 6 Article 10.3390/rs70607571 7571 7596 2072-4292 2015-06-08 doi: 10.3390/rs70607571 Simon Zwieback Scott Hensley Irena Hajnsek http://www.mdpi.com/2072-4292/7/6/7545 Desert locusts (Schistocerca gregaria) represent a major threat for agro-pastoral resources and food security over almost 30 million km2 from northern Africa to the Arabian peninsula and India. Given the differential food preferences of this insect pest and the extent and remoteness of the their distribution area, near-real-time remotely-sensed information on potential habitats support control operations by narrowing down field surveys to areas favorable for their development and prone to gregarization and outbreaks. The development of dynamic greenness maps, which detect the onset of photosynthetic vegetation, allowed national control centers to identify potential habitats to survey, as locusts prefer green and fresh vegetation. Their successful integration into the daily control operations led to a new need: the near-real-time identification of the onset of dryness, a synonym for the loss of habitat attractiveness, likely to be abandoned by locusts. The timely availability of this information would enable control centers to focus their surveys on areas more prone to gregarization, leading to more efficiency in the allocation of resources and in decision making. In this context, this work developed an original method to detect in near-real-time the onset of vegetation senescence. The design of the detection relies on the temporal behavior of two indices: the Normalized Difference Vegetation Index, depending on the green vegetation, and the Normalized Difference Tillage Index, sensitive to both green and dry vegetation. The method is demonstrated in Mauritania, an ever-affected country, with 10-day MODIS mean composites for the years 2010 and 2011. The discrimination performance of three classes ("growth", "density reduction" and "drying") were analyzed for three classification methods: maximum likelihood (61.4% of overall accuracy), decision tree (71.5%) and support vector machine (72.3%). The classification accuracy is heterogeneous in both time and space and is affected by several factors, such as vegetation density, the north-south climatic gradient and the relief. Smoothing the vegetation time series resulted in an increase of the overall accuracy of about 5% at the expense of a loss in timeliness of ten days. To simulate near-real-time monitoring conditions, the decision tree was applied to the decade of 2010. Overall, the seasonal vegetation cycle appeared clear and consistent. The results obtained pave the way for an operational implementation of the senescence dynamic mapping and, consequently, to further strengthen the capacity of the locust control management. Remote Sensing 2015-06-08 7 6 Article 10.3390/rs70607545 7545 7570 2072-4292 2015-06-08 doi: 10.3390/rs70607545 Cécile Renier François Waldner Damien Jacques Mohamed Babah Ebbe Keith Cressman Pierre Defourny http://www.mdpi.com/2072-4292/7/6/7521 Especially in the remote sensing context, thematic classification is a desired product for coral reef surveys. This study presents a novel statistical-based image classification approach, namely Partial Least Square Discriminant Analysis (PLS-DA), capable of doing so. Three classification models were built and implemented for the images while the fourth was a combination of spectra from all three images together. The classification was optimised by using pre-processing transformations (PPTs) and post-classification low-pass filtering. Despite the fact that the images were acquired under different conditions and quality, the best classification model was achieved by combining spectral training samples from three images (accuracy 0.63 for all classes). PPTs improved the classification accuracy by 5%–15% and post-classification treatments further increased the final accuracy by 10%–20%. The fourth classification model was the most accurate one, suggesting that combining spectra from differ conditions improves thematic classification. Despite some limitations, available aerial sensors already provide an opportunity to implement the described classification and mark the next investigation step. Nonetheless, the findings of this study are relevant both to the field of remote sensing in general and to the niche of coral reef spectroscopy. Remote Sensing 2015-06-05 7 6 Article 10.3390/rs70607521 7521 7544 2072-4292 2015-06-05 doi: 10.3390/rs70607521 Tamir Caras Arnon Karnieli http://www.mdpi.com/2072-4292/7/6/7494 Dongting Lake, the second largest freshwater lake in China, is well known for its rapid seasonal fluctuations in inundation extents in the middle reach of the Yangtze River, and it is also the lake most affected by the Three Gorges Project. Significant inter-annual and seasonal variations in flood inundations were observed from Moderate Resolution Imaging Spectroradiometer (MODIS) time-series imagery between 2000 and 2012 in the Dongting Lake. Results demonstrated that temporal changes in inundation extents derived from MODIS data were accordant with variations in annual and monthly precipitation and runoff data. Spatial and temporal dynamics of some related parameters of flood regime were analyzed as well, which included flood inundation probability, duration and start/end date of the annual largest flood. Large areas with high flood inundation probability were identified in 2000 and 2002, but relatively small regions with great flood inundation probability occurred in 2001, 2006, and 2011. Long flood durations were observed in 2000, 2002, 2008, 2010, and 2012, whereas short flood durations occurred in 2001, 2006, and 2011. Correlation analysis techniques were applied to explore spatial-temporal relationships between parameters associated with flood regime and wetland landscape patterns from 2000 to 2012. In addition, this paper presented comprehensive discussions on development of related parameters of flood regime and their influences on wetland landscape pattern after impoundment of the Three Gorges Reservoir, changes in wetland landscape patterns after the flood period, and the role of flooding in wetland evolution and vegetation succession. These results can provide scientific guidance and baseline data for wetland management and long-term monitoring of wetland ecological environment in the Dongting Lake. Remote Sensing 2015-06-05 7 6 Article 10.3390/rs70607494 7494 7520 2072-4292 2015-06-05 doi: 10.3390/rs70607494 Yanxia Hu Jinliang Huang Yun Du Pengpeng Han Wei Huang http://www.mdpi.com/2072-4292/7/6/7471 The scattering of GNSS signals over a water surface is studied when the receiver is at a low height, as in GNSS-R coastal altimetry. The precise determination of the local sea level and wave state from the coast will provide useful altimetry and wave information as "dry" tide and wave gauges. An experiment has been conducted at the Canal d'Investigació i Experimentació Marítima (CIEM) wave channel for two simulated "sea" states. The GNSS-reflectometer used is the P(Y) and C/A ReflectOmeter (PYCARO) instrument, a closed-loop receiver with delay and Doppler tracking loops that uses the conventional GNSS-R technique for the GPS C/A code. After retracking of the scattered GPS signals, the coherent and incoherent components have been studied. To reproduce the transmitted GPS signals indoors, a Rohde and Schwarz signal generator is used. It is found that, despite the ratio of the coherent and incoherent components being ~1, the coherent component is strong enough that it can be tracked. The coherent component comes from clusters of points on the surface that approximately satisfy the specular reflection conditions ("roughed facet"). The Pearson's linear correlation coefficients of the derived "sea" surface height with the wave gauge data are: 0.78, 0.85 and 0.81 for a SWH = 36 cm and 0.34, 0.74, and 0.72 for a SWH = 64 cm, respectively, for transmitter elevation angles of = 60°, 75° and 86°, respectively. Finally, the rms phase of the received signal before the retracking processing is used to estimate the effective rms surface height of the 'facets', where the waves get scattered. It is found to be between 2.5- and 4.1-times smaller than the theoretical values corresponding to the half of the coherent reflectivity decaying factor. Remote Sensing 2015-06-05 7 6 Article 10.3390/rs70607471 7471 7493 2072-4292 2015-06-05 doi: 10.3390/rs70607471 Hugo Carreno-Luengo Adriano Camps http://www.mdpi.com/2072-4292/7/6/7447 Dual polarization is a typical operational mode of polarimetric synthetic aperture radar (SAR). However, few studies have considered the scattering mechanism extraction of dual-polarization SARs. A modified Cloude-Pottier decomposition is proposed to investigate the performance of the scattering mechanism extraction of dual-polarization SARs. It is theoretically demonstrated that only HH-VV SAR can discriminate the three canonical scattering mechanisms from an isotropic surface, horizontal dipole, and isotropic dihedral. Various experiments are conducted using 21 scenes from real datasets acquired by AIRSAR, Convair-580 SAR, EMISAR, E-SAR, Pi-SAR, and RADARSAT-2. Division of the dual-polarization H-α plane is experimentally obtained. The lack of cross-polarization induces the diffusion of scattering mechanisms and their overlap in the HH-VV H-α plane. However, the performance of HH-VV SAR for extracting scattering mechanisms is acceptable. Thus, HH-VV SAR is a suitable alternative to full-polarization SAR in certain cases. Meanwhile, the extraction performance of the other two dual-polarization SARs is badly degraded due to the lack of co-polarization. Therefore, HH-HV and HV-VV SARs cannot effectively extract the scattering mechanisms in the H-α plane. Remote Sensing 2015-06-05 7 6 Article 10.3390/rs70607447 7447 7470 2072-4292 2015-06-05 doi: 10.3390/rs70607447 Kefeng Ji Yonghui Wu http://www.mdpi.com/2072-4292/7/6/7425 The estimation of the Fraction of Absorbed Photosynthetically Active Radiation in forests (forest fAPAR) from multi-spectral Landsat-8 data is investigated in this paper using a physically based radiative transfer model (Invertible Forest Reflectance Model, INFORM) combined with an inversion strategy based on artificial neural nets (ANN). To derive the forest fAPAR for the Dabie mountain test site in China in 30 m spatial resolution (size approximately 3000 km2), a database of forest canopy spectral reflectances was simulated with INFORM taking into account structural variables such as leaf area index (LAI), crown coverage and stem density as well as leaf composition. To establish the relationship between forest fAPAR and the reflectance modeled by INFORM, a logarithmic relationship between LAI and fAPAR was used previously established using on-site field measurements. On this basis, predictive models between Landsat-8 reflectance and fAPAR were established using an artificial neural network. After calibrating INFORM for the test site, forty-two forest stands were used to validate the performance of the method. The results show that spectral signatures modeled by INFORM correspond reasonably well with the forest canopy reflectance spectra derived from Landsat data. Deviations increase with increasing angle between surface normal of the hilly terrain and sun incidence. The comparison of estimated and measured fAPAR (R2 = 0.47, RMSE = 0.11) demonstrates that INFORM can be inverted using neural nets to provide acceptable estimates of forest fAPAR. The accuracy of the predictions increased significantly when excluding pixels located in very steep terrain. This demonstrates that the applied topographic correction was not sufficiently accurate and should be improved for making optimum use of radiative transfer models such as INFORM. Remote Sensing 2015-06-05 7 6 Article 10.3390/rs70607425 7425 7446 2072-4292 2015-06-05 doi: 10.3390/rs70607425 Huili Yuan Ronghua Ma Clement Atzberger Fei Li Steven Loiselle Juhua Luo http://www.mdpi.com/2072-4292/7/6/7402 In environments that are hostile to Global Navigation Satellites Systems (GNSS), the precision achieved by a mobile light detection and ranging (LiDAR) system (MLS) can deteriorate into the sub-meter or even the meter range due to errors in the positioning and orientation system (POS). This paper proposes a novel least squares collocation (LSC)-based method to improve the accuracy of the MLS in these hostile environments. Through a thorough consideration of the characteristics of POS errors, the proposed LSC-based method effectively corrects these errors using LiDAR control points, thereby improving the accuracy of the MLS. This method is also applied to the calibration of misalignment between the laser scanner and the POS. Several datasets from different scenarios have been adopted in order to evaluate the effectiveness of the proposed method. The results from experiments indicate that this method would represent a significant improvement in terms of the accuracy of the MLS in environments that are essentially hostile to GNSS and is also effective regarding the calibration of misalignment. Remote Sensing 2015-06-03 7 6 Article 10.3390/rs70607402 7402 7424 2072-4292 2015-06-03 doi: 10.3390/rs70607402 Qingzhou Mao Liang Zhang Qingquan Li Qingwu Hu Jianwei Yu Shaojun Feng Washington Ochieng Hanlu Gong http://www.mdpi.com/2072-4292/7/6/7378 Greenhouse detection and mapping via remote sensing is a complex task, which has already been addressed in numerous studies. In this research, the innovative goal relies on the identification of greenhouse horticultural crops that were growing under plastic coverings on 30 September 2013. To this end, object-based image analysis (OBIA) and a decision tree classifier (DT) were applied to a set consisting of eight Landsat 8 OLI images collected from May to November 2013. Moreover, a single WorldView-2 satellite image acquired on 30 September 2013, was also used as a data source. In this approach, basic spectral information, textural features and several vegetation indices (VIs) derived from Landsat 8 and WorldView-2 multi-temporal satellite data were computed on previously segmented image objects in order to identify four of the most popular autumn crops cultivated under greenhouse in Almería, Spain (i.e., tomato, pepper, cucumber and aubergine). The best classification accuracy (81.3% overall accuracy) was achieved by using the full set of Landsat 8 time series. These results were considered good in the case of tomato and pepper crops, being significantly worse for cucumber and aubergine. These results were hardly improved by adding the information of the WorldView-2 image. The most important information for correct classification of different crops under greenhouses was related to the greenhouse management practices and not the spectral properties of the crops themselves. Remote Sensing 2015-06-03 7 6 Article 10.3390/rs70607378 7378 7401 2072-4292 2015-06-03 doi: 10.3390/rs70607378 Manuel Aguilar Andrea Vallario Fernando Aguilar Andrés Lorca Claudio Parente http://www.mdpi.com/2072-4292/7/6/7350 Search-centric, sample supervised image segmentation has been demonstrated as a viable general approach applicable within the context of remote sensing image analysis. Such an approach casts the controlling parameters of image processing—generating segments—as a multidimensional search problem resolvable via efficient search methods. In this work, this general approach is analyzed in the context of connected component segmentation. A specific formulation of connected component labeling, based on quasi-flat zones, allows for the addition of arbitrary segment attributes to contribute to the nature of the output. This is in addition to core tunable parameters controlling the basic nature of connected components. Additional tunable constituents may also be introduced into such a framework, allowing flexibility in the definition of connected component connectivity, either directly via defining connectivity differently or via additional processes such as data mapping functions. The relative merits of these two additional constituents, namely the addition of tunable attributes and data mapping functions, are contrasted in a general remote sensing image analysis setting. Interestingly, tunable attributes in such a context, conjectured to be safely useful in general settings, were found detrimental under cross-validated conditions. This is in addition to this constituent's requiring substantially greater computing time. Casting connectivity definitions as a searchable component, here via the utilization of data mapping functions, proved more beneficial and robust in this context. The results suggest that further investigations into such a general framework could benefit more from focusing on the aspects of data mapping and modifiable connectivity as opposed to the utility of thresholding various geometric and spectral attributes. Remote Sensing 2015-06-03 7 6 Article 10.3390/rs70607350 7350 7377 2072-4292 2015-06-03 doi: 10.3390/rs70607350 Christoff Fourie http://www.mdpi.com/2072-4292/7/6/7324 The mapping and forecasting of droughts and floods is an important potential field of application of global soil moisture and water storage products from satellites and models. Especially when extremes in near-surface soil moisture propagate into extremes in total water storage, agricultural production and water supply can be severely impacted. This study relates soil moisture from the WaterGAP Global Hydrology Model (WGHM) and the satellite sensors Advanced Microwave Scanning Radiometer—Earth Observing System (AMSR-E) and Advanced Scatterometer (ASCAT) to total water storage variations from the satellite gravity mission GRACE. A particular focus is on destructive hydrological extreme events, as listed in the International Disaster Database EM-DAT. Data sets are analyzed via correlation, time shift, and principal component analyses. The study area is the La Plata Basin in South America. The results indicate that most of the soil moisture anomalies are linked to periods of El Niño and La Niña and associated natural disasters. For the La Plata drought of 2008/2009 and the El Niño flooding of 2009/2010, soil moisture serves as an indicator for the later deficit or surplus in total water storage. These hydrological anomalies were strongest in the southern, central, and eastern parts of the basin, but more than one hundred thousand people were also affected in the northwestern part. Remote Sensing 2015-06-03 7 6 Article 10.3390/rs70607324 7324 7349 2072-4292 2015-06-03 doi: 10.3390/rs70607324 Sarah Abelen Florian Seitz Rodrigo Abarca-del-Rio Andreas Güntner http://www.mdpi.com/2072-4292/7/6/7298 Managers of forested water supply catchments require efficient and accurate methods to quantify changes in forest water use due to changes in forest structure and density after disturbance. Using Light Detection and Ranging (LiDAR) data with as few as 0.9 pulses m−2, we applied a local maximum filtering (LMF) method and normalised cut (NCut) algorithm to predict stocking density (SDen) of a 69-year-old Eucalyptus regnans forest comprising 251 plots with resolution of the order of 0.04 ha. Using the NCut method we predicted basal area (BAHa) per hectare and sapwood area (SAHa) per hectare, a well-established proxy for transpiration. Sapwood area was also indirectly estimated with allometric relationships dependent on LiDAR derived SDen and BAHa using a computationally efficient procedure. The individual tree detection (ITD) rates for the LMF and NCut methods respectively had 72% and 68% of stems correctly identified, 25% and 20% of stems missed, and 2% and 12% of stems over-segmented. The significantly higher computational requirement of the NCut algorithm makes the LMF method more suitable for predicting SDen across large forested areas. Using NCut derived ITD segments, observed versus predicted stand BAHa had R2 ranging from 0.70 to 0.98 across six catchments, whereas a generalised parsimonious model applied to all sites used the portion of hits greater than 37 m in height (PH37) to explain 68% of BAHa. For extrapolating one ha resolution SAHa estimates across large forested catchments, we found that directly relating SAHa to NCut derived LiDAR indices (R2 = 0.56) was slightly more accurate but computationally more demanding than indirect estimates of SAHa using allometric relationships consisting of BAHa (R2 = 0.50) or a sapwood perimeter index, defined as (BAHaSDen)½ (R2 = 0.48). Remote Sensing 2015-06-03 7 6 Article 10.3390/rs70607298 7298 7323 2072-4292 2015-06-03 doi: 10.3390/rs70607298 Dominik Jaskierniak George Kuczera Richard Benyon Luke Wallace http://www.mdpi.com/2072-4292/7/6/7272 As the largest natural source of methane, wetlands play an important role in the carbon cycle. High-resolution maps of wetland type and extent are required to quantify wetland responses to climate change. Mapping northern wetlands is particularly important because of a disproportionate increase in temperatures at higher latitudes. Synthetic aperture radar data from a spaceborne platform can be used to map wetland types and dynamics over large areas. Following from earlier work by Whitcomb et al. (2009) using Japanese Earth Resources Satellite (JERS-1) data, we applied the "random forests" classification algorithm to variables from L-band ALOS PALSAR data for 2007, topographic data (e.g., slope, elevation) and locational information (latitude, longitude) to derive a map of vegetated wetlands in Alaska, with a spatial resolution of 50 m. We used the National Wetlands Inventory and National Land Cover Database (for upland areas) to select training and validation data and further validated classification results with an independent dataset that we created. A number of improvements were made to the method of Whitcomb et al. (2009): (1) more consistent training data in upland areas; (2) better distribution of training data across all classes by taking a stratified random sample of all available training pixels; and (3) a more efficient implementation, which allowed classification of the entire state as a single entity (rather than in separate tiles), which eliminated discontinuities at tile boundaries. The overall accuracy for discriminating wetland from upland was 95%, and the accuracy at the level of wetland classes was 85%. The total area of wetlands mapped was 0.59 million km2, or 36% of the total land area of the state of Alaska. The map will be made available to download from NASA's wetland monitoring website. Remote Sensing 2015-06-03 7 6 Article 10.3390/rs70607272 7272 7297 2072-4292 2015-06-03 doi: 10.3390/rs70607272 Daniel Clewley Jane Whitcomb Mahta Moghaddam Kyle McDonald Bruce Chapman Peter Bunting http://www.mdpi.com/2072-4292/7/6/7257 This paper presents a technique to generate cirrus optical depth and particle effective size estimates from the cloud emissivities at 8.5, 11 and 12 μm contained in the Collection-6 (C6) MYD06 cloud product. This technique employs the latest scattering models and scattering radiative transfer approximations to estimate cloud optical depth and particle effective size using efficient analytical formulae. Two scattering models are tested. The first is the same scattering model as that used in the C6 MYD06 solar reflectance products. The second model is an empirical model derived from radiometric consistency. Both models are shown to generate optical depths that compare well to those from constrained CALIPSO retrievals and MYD06. In terms of effective radius retrievals, the results from the radiometric empirical model agree more closely with MYD06 than those from the C6 model. This analysis is applied to AQUA/MODIS data collocated with CALIPSO/CALIOP during January 2010. Remote Sensing 2015-06-03 7 6 Article 10.3390/rs70607257 7257 7271 2072-4292 2015-06-03 doi: 10.3390/rs70607257 Andrew Heidinger Yue Li Bryan Baum Robert Holz Steven Platnick Ping Yang http://www.mdpi.com/2072-4292/7/6/7231 Significant economic losses, large affected populations, and serious environmental damage caused by recurrent natural disaster events (NDE) worldwide indicate insufficiency in emergency preparedness and response. The barrier of full life cycle data preparation and information support is one of the main reasons. This paper adopts the method of integrated environmental modeling, incorporates information from existing event protocols, languages, and models, analyzes observation demands from different event stages, and forms the abstract full life cycle natural disaster event metamodel (FLCNDEM) based on meta-object facility. Then task library and knowledge base for floods are built to instantiate FLCNDEM, forming the FLCNDEM for floods (FLCNDEMF). FLCNDEMF is formalized according to Event Pattern Markup Language, and a prototype system, Natural Disaster Event Manager, is developed to assist in the template-based modeling and management. The flood in Liangzi (LZ) Lake of Hubei, China on 16 July 2010 is adopted to illustrate how to apply FLCNDEM in real scenarios. FLCNDEM-based modeling is realized, and the candidate remote sensing (RS) dataset for different observing missions are provided for LZ Lake flood. Taking the mission of flood area extraction as an example, the appropriate RS data are selected via the model of simplified general perturbation version 4, and the flood area in different phases are calculated and displayed on the map. The phase-based modeling and visualization intuitively display the spatial-temporal distribution and the evolution process of the LZ Lake flood, and it is of great significance for flood responding. In addition, through the extension mechanism, FLCNDEM can also be applied in other environmental applications, providing important support for full life cycle information sharing and rapid responding. Remote Sensing 2015-06-03 7 6 Article 10.3390/rs70607231 7231 7256 2072-4292 2015-06-03 doi: 10.3390/rs70607231 Nengcheng Chen Wenying Du Fan Song Zeqiang Chen http://www.mdpi.com/2072-4292/7/6/7212 Long-term snow depth/snow water equivalent (SWE) products derived from passive microwave remote sensing data are fundamental for climatological and hydrological studies. However, the temporal continuity of the products is affected by the updating or replacement of passive microwave sensors or satellite platforms. In this study, we inter-calibrated brightness temperature (Tb) data obtained from the Special Sensor Microwave Imager (SSM/I) and the Special Sensor Microwave Imager/Sounder (SSMI/S). Then, we evaluated the consistency of the snow cover area (SCA) and snow depth derived from the Scanning Multichannel Microwave Radiometer (SMMR), SSM/I and SSMI/S. The results indicated that (1) the spatial pattern of the SCA derived from the SMMR and SSM/I data was more consistent after calibration than before; (2) the relative biases in the SCA and snow depth in China between the SSM/I and SSMI/S data decreased from 42.42% to 1.65% and from 66.18% to −1.5%, respectively; and (3) the SCA and snow depth derived from the SSM/I data carried on F08, F11 and F13 were highly consistent. To obtain consistent snow depth and SCA products, inter-sensor calibrations between SMMR, SSM/I and SSMI/S are important. In consideration of the snow data product continuation, we suggest that the brightness temperature data from all sensors be calibrated based on SSMI/S. Remote Sensing 2015-06-02 7 6 Article 10.3390/rs70607212 7212 7230 2072-4292 2015-06-02 doi: 10.3390/rs70607212 Liyun Dai Tao Che Yongjian Ding http://www.mdpi.com/2072-4292/7/6/7181 This paper examines the spatial error structures of eight precipitation estimates derived from four different satellite retrieval algorithms including TRMM Multi-satellite Precipitation Analysis (TMPA), Climate Prediction Center morphing technique (CMORPH), Global Satellite Mapping of Precipitation (GSMaP) and Precipitation Estimation from Remotely Sensed Information using Artificial Neural Networks (PERSIANN). All the original satellite and bias-corrected products of each algorithm (3B42RTV7 and 3B42V7, CMORPH_RAW and CMORPH_CRT, GSMaP_MVK and GSMaP_Gauge, PERSIANN_RAW and PERSIANN_CDR) are evaluated against ground-based Asian Precipitation-Highly Resolved Observational Data Integration Towards Evaluation of Water Resources (APHRODITE) over Central Asia for the period of 2004 to 2006. The analyses show that all products except PERSIANN exhibit overestimation over Aral Sea and its surrounding areas. The bias-correction improves the quality of the original satellite TMPA products and GSMaP significantly but slightly in CMORPH and PERSIANN over Central Asia. 3B42RTV7 overestimates precipitation significantly with large Relative Bias (RB) (128.17%) while GSMaP_Gauge shows consistent high correlation coefficient (CC) (>0.8) but RB fluctuates between −57.95% and 112.63%. The PERSIANN_CDR outperforms other products in winter with the highest CC (0.67). Both the satellite-only and gauge adjusted products have particularly poor performance in detecting rainfall events in terms of lower POD (less than 65%), CSI (less than 45%) and relatively high FAR (more than 35%). Remote Sensing 2015-06-01 7 6 Article 10.3390/rs70607181 7181 7211 2072-4292 2015-06-01 doi: 10.3390/rs70607181 Hao Guo Sheng Chen Anming Bao Jujun Hu Abebe Gebregiorgis Xianwu Xue Xinhua Zhang http://www.mdpi.com/2072-4292/7/6/7157 New calculations of the relative optical air mass function are made over the 0°–87° range of apparent solar zenith angle θ, for various vertical profiles of background aerosol, diamond dust and thin cirrus cloud particle extinction coefficient in the Arctic and Antarctic atmospheres. The calculations were carried out by following the Tomasi and Petkov (2014) procedure, in which the above-mentioned vertical profiles derived from lidar observations were used as weighting functions. Different sets of lidar measurements were examined, recorded using: (i) the Koldewey-Aerosol-Raman Lidar (KARL) system (AWI, Germany) at Ny-Ålesund (Spitsbergen, Svalbard) in January, April, July and October 2013; (ii) the Cloud-Aerosol Lidar and Infrared Pathfinder Satellite Observations (CALIPSO) satellite-based sensor over Barrow (Alaska), Eureka (Nunavut, Canada) and Sodankylä (northern Finland), and Neumayer III, Mario Zucchelli and Mirny coastal stations in Antarctica in the local summer months of the last two years; (iii) the National Institute of Optics (INO), National Council of Research (CNR) Antarctic lidar at Dome C on the Antarctic Plateau for a typical "diamond dust" case; and (iv) the KARL lidar at Ny-Ålesund and the University of Rome/National Agency for New Technologies, Energy and Sustainable Economic Development (ENEA) lidar at Thule (northwestern Greenland) for some cirrus cloud layers in the middle and upper troposphere. The relative optical air mass calculations are compared with those obtained by Tomasi and Petkov (2014) to define the seasonal changes produced by aerosol particles, diamond dust and cirrus clouds. The results indicate that the corresponding air mass functions generally decrease as angle θ increases with rates that are proportional to the increase in the pure aerosol, diamond dust and cirrus cloud particle optical thickness. Remote Sensing 2015-06-01 7 6 Article 10.3390/rs70607157 7157 7180 2072-4292 2015-06-01 doi: 10.3390/rs70607157 Claudio Tomasi Boyan Petkov Mauro Mazzola Christoph Ritter Alcide di Sarra Tatiana di Iorio Massimo del Guasta http://www.mdpi.com/2072-4292/7/6/7126 Validation and performance evaluations are beneficial for developing methods that estimate the remotely sensed land surface temperature (LST). However, such evaluations for Advanced Spaceborne Thermal Emission and Reflection Radiometer (ASTER) data are rare. By selecting the middle reach of the Heihe River basin (HRB), China, as the study area, the atmospheric correction (AC), mono-window (MW), single-channel (SC), and split-window (SW) methods were evaluated based on in situ measured LSTs. Results demonstrate that the influences of surface heterogeneity on the validation are significant in the study area. For the AC, MW, and SC methods, the LSTs estimated from channel 13 are more accurate than those from channel 14 in general cases. When the in situ measured atmospheric profiles are available, the AC method has the highest accuracy, with a root-mean squared error (RMSE) of about 1.4–1.5 K at the homogenous oasis sites. In actual application without sufficient in situ measured inputs, the MW method is highly accurate; the RMSE is around 1.5–1.6 K. The SC method systematically overestimates LSTs and it is sensitive to error in the water vapor content. The two SW methods are simple to use but their performances are limited by accuracies, revealed by the simulation dataset. Therefore, when the in situ atmospheric profiles are available, the AC method is recommended to generate reliable ASTER LSTs for modeling the eco-hydrological processes in the middle reach of the HRB. When sufficient in situ measured inputs are not available, the MW method can be used instead. Remote Sensing 2015-05-29 7 6 Article 10.3390/rs70607126 7126 7156 2072-4292 2015-05-29 doi: 10.3390/rs70607126 Ji Zhou Mingsong Li Shaomin Liu Zhenzhen Jia Yanfei Ma http://www.mdpi.com/2072-4292/7/6/7105 In view of the fact that oil spill remote sensing could only generate the oil slick information at a specific time and that traditional oil spill simulation models were not designed to deal with dynamic conditions, a dynamic data-driven application system (DDDAS) was introduced. The DDDAS entails both the ability to incorporate additional data into an executing application and, in reverse, the ability of applications to dynamically steer the measurement process. Based on the DDDAS, combing a remote sensor system that detects oil spills with a numerical simulation, an integrated data processing, analysis, forecasting and emergency response system was established. Once an oil spill accident occurs, the DDDAS-based oil spill model receives information about the oil slick extracted from the dynamic remote sensor data in the simulation. Through comparison, information fusion and feedback updates, continuous and more precise oil spill simulation results can be obtained. Then, the simulation results can provide help for disaster control and clean-up. The Penglai, Xingang and Suizhong oil spill results showed our simulation model could increase the prediction accuracy and reduce the error caused by empirical parameters in existing simulation systems. Therefore, the DDDAS-based detection and simulation system can effectively improve oil spill simulation and diffusion forecasting, as well as provide decision-making information and technical support for emergency responses to oil spills. Remote Sensing 2015-05-29 7 6 Article 10.3390/rs70607105 7105 7125 2072-4292 2015-05-29 doi: 10.3390/rs70607105 Jining Yan Lizhe Wang Lajiao Chen Lingjun Zhao Bomin Huang http://www.mdpi.com/2072-4292/7/6/7080 This paper uses the refined Generalized Split-Window (GSW) algorithm to derive the land surface temperature (LST) from the data acquired by the Visible and Infrared Radiometer on FengYun 3B (FY-3B/VIRR). The coefficients in the GSW algorithm corresponding to a series of overlapping ranges for the mean emissivity, the atmospheric Water Vapor Content (WVC), and the LST are derived using a statistical regression method from the numerical values simulated with an accurate atmospheric radiative transfer model MODTRAN 4 over a wide range of atmospheric and surface conditions. The GSW algorithm is applied to retrieve LST from FY-3B/VIRR data in an arid area in northwestern China. Three emissivity databases are used to evaluate the accuracy of different emissivity databases for LST retrieval, including the ASTER Global Emissivity Database (ASTER_GED) at a 1-km spatial resolution (AG1km), an average of twelve ASTER emissivity data in the 2012 summer and emissivity spectra extracted from spectral libraries. The LSTs retrieved from the three emissivity databases are evaluated with ground-measured LST at four barren surface sites from June 2012 to December 2013 collected during the HiWATER field campaign. The results indicate that using emissivity extracted from ASTER_GED can achieve the highest accuracy with an average bias of 1.26 and −0.04 K and an average root mean square error (RMSE) of 2.69 and 1.38 K for the four sites during daytime and nighttime, respectively. This result indicates that ASTER_GED is a useful emissivity database for generating global LST products from different thermal infrared data and that using FY-3B/VIRR data can produce reliable LST products for other research areas. Remote Sensing 2015-05-29 7 6 Article 10.3390/rs70607080 7080 7104 2072-4292 2015-05-29 doi: 10.3390/rs70607080 Jinxiong Jiang Hua Li Qinhuo Liu Heshun Wang Yongming Du Biao Cao Bo Zhong Shanlong Wu http://www.mdpi.com/2072-4292/7/6/7062 The cumulative error at a point in a LiDAR-derived DEM consists of three components: propagated LiDAR-sensor error, propagated ground error, and interpolation error. To combine these error components so as to assess the vertical accuracy of a LiDAR-derived DEM, statistical methods based on the error propagation theory are often used. Due to the existence of systematic error, statistical methods are only effective if a large number of checkpoints are available, which may not be affordable in many practical applications. This paper presents approximation theory as an alternative methodology that departs from error propagation theory in fundamental ways. Using approximation theory, an error bound of the cumulative error at any point in the study site can be obtained, thus informing users conservatively of the spatial variation of DEM accuracy and pointing out the weakly determined areas. The new method is illustrated from DEM users' perspective by assessing whether a publicly available LiDAR-derived DEM meets FEMA's accuracy standard for flood risk mapping. The paper calls for a change in the existing methods of assessing and reporting the errors in a LiDAR-derived DEM, in particular those introduced during the ground filtering process. Remote Sensing 2015-05-29 7 6 Article 10.3390/rs70607062 7062 7079 2072-4292 2015-05-29 doi: 10.3390/rs70607062 XiaoHang Liu Hai Hu Peng Hu http://www.mdpi.com/2072-4292/7/6/7044 Imagery registration is a fundamental step, which greatly affects later processes in image mosaic, multi-spectral image fusion, digital surface modelling, etc., where the final solution needs blending of pixel information from more than one images. It is highly desired to find a way to identify registration regions among input stereo image pairs with high accuracy, particularly in remote sensing applications in which ground control points (GCPs) are not always available, such as in selecting a landing zone on an outer space planet. In this paper, a framework for localization in image registration is developed. It strengthened the local registration accuracy from two aspects: less reprojection error and better feature point distribution. Affine scale-invariant feature transform (ASIFT) was used for acquiring feature points and correspondences on the input images. Then, a homography matrix was estimated as the transformation model by an improved random sample consensus (IM-RANSAC) algorithm. In order to identify a registration region with a better spatial distribution of feature points, the Euclidean distance between the feature points is applied (named the S criterion). Finally, the parameters of the homography matrix were optimized by the Levenberg–Marquardt (LM) algorithm with selective feature points from the chosen registration region. In the experiment section, the Chang'E-2 satellite remote sensing imagery was used for evaluating the performance of the proposed method. The experiment result demonstrates that the proposed method can automatically locate a specific region with high registration accuracy between input images by achieving lower root mean square error (RMSE) and better distribution of feature points. Remote Sensing 2015-05-29 7 6 Article 10.3390/rs70607044 7044 7061 2072-4292 2015-05-29 doi: 10.3390/rs70607044 Xiangjun Wang Yang Li Hong Wei Feng Liu http://www.mdpi.com/2072-4292/7/6/7029 To meet growing food demand with limited land and reduced environmental impact, soil testing and formulated fertilization methods have been widely adopted around the world. However, conventional technology for investigating nitrogen fertilization rates (NFR) is time consuming and expensive. Here, we evaluated the use of visible near-infrared shortwave-infrared (VIS-NIR-SWIR: 400–2500 nm) spectroscopy for the assessment of NFR to provide necessary information for fast, cost-effective and precise fertilization rating. Over 2000 samples were collected from paddy-rice fields in 10 Chinese provinces; samples were added to the Chinese Soil Spectral Library (CSSL). Two kinds of modeling strategies for NFR, quantitative estimation of soil N prior to classification and qualitative by classification, were employed using partial least squares regression (PLSR), locally weighted regression (LWR), and support vector machine discriminant analogy (SVMDA). Overall, both LWR and SVMDA had moderate accuracies with Cohen's kappa coefficients of 0.47 and 0.48, respectively, while PLSR had fair accuracy (0.37). We conclude that VIS-NIR-SWIR spectroscopy coupled with the CSSL appears to be a viable, rapid means for the assessment of NFR in paddy-rice soil. Based on qualitative classification of soil spectral data only, it is recommended that the SVMDA be adopted for rapid implementation. Remote Sensing 2015-05-29 7 6 Article 10.3390/rs70607029 7029 7043 2072-4292 2015-05-29 doi: 10.3390/rs70607029 Shuo Li Wenjun Ji Songchao Chen Jie Peng Yin Zhou Zhou Shi http://www.mdpi.com/2072-4292/7/6/7007 Assessing nitrogen (N) status is important from economic and environmental standpoints. To date, many spectral indices to estimate cotton chlorophyll or N content have been purely developed using statistical analysis approach where they are often subject to site-specific problems. This study describes and tests a novel method of utilizing physical characteristics of N-fertilized cotton and combining field spectral measurements made at different spatial scales as an approach to estimate in-season chlorophyll or leaf N content of field-grown cotton. In this study, leaf greenness estimated from spectral measurements made at the individual leaf, canopy and scene levels was combined with percent ground cover to produce three different indices, named TCCLeaf, TCCCanopy, and TCCScene. These indices worked best for estimating leaf N at early flowering, but not for chlorophyll content. Of the three indices, TCCLeaf showed the best ability to estimate leaf N (R2 = 0.89). These results suggest that the use of green and red-edge wavelengths derived at the leaf scale is best for estimating leaf greenness. TCCCanopy had a slightly lower R2 value than TCCLeaf (0.76), suggesting that the utilization of yellow and red-edge wavelengths obtained at the canopy level could be used as an alternative to estimate leaf N in the absence of leaf spectral information. The relationship between TCCScene and leaf N was the lowest (R2 = 0.50), indicating that the estimation of canopy greenness from scene measurements needs improvement. Results from this study confirmed the potential of these indices as efficient methods for estimating in-season leaf N status of cotton. Remote Sensing 2015-05-29 7 6 Article 10.3390/rs70607007 7007 7028 2072-4292 2015-05-29 doi: 10.3390/rs70607007 Farrah Muharam Stephen Maas Kevin Bronson Tina Delahunty http://www.mdpi.com/2072-4292/7/6/6986 With the high resolution of optical data and the lack of weather effects of passive microwave data, we developed an algorithm to map daily cloud-free fractional snow cover (FSC) based on the Moderate Resolution Imaging Spectroradiometer (MODIS) standard daily FSC product, the Advanced Microwave Scanning Radiometer (AMSR2) snow water equivalent (SWE) product and digital elevation data. We then used the algorithm to produce a daily cloud-free FSC product with a resolution of 500 m for regions in China. In addition, we produced a high-resolution FSC map using a Landsat 8 Operational Land Imager (OLI) image as a true value to test the accuracy of the cloud-free FSC product developed in this study. The analysis results show that the daily cloud-free FSC product developed in this study can completely remove clouds and effectively improve the accuracy of snow area monitoring. Compared to the true value, the mean absolute error of our product is 0.20, and its root mean square error is 0.29. Thus, the synthesized product in this study can improve the accuracy of snow area monitoring, and the obtained snow area data can be used as reliable input parameters for hydrological and climate models. The land cover type and terrain factors are the main factors that limit the accuracy of the daily cloud-free FSC product developed in this study. These limitations can be further improved by improving the accuracy of the MODIS standard snow product for complicated underlying surfaces. Remote Sensing 2015-05-29 7 6 Article 10.3390/rs70606986 6986 7006 2072-4292 2015-05-29 doi: 10.3390/rs70606986 Jie Deng Xiaodong Huang Qisheng Feng Xiaofang Ma Tiangang Liang http://www.mdpi.com/2072-4292/7/6/6950 Typically, digital image processing for burned-areas detection combines the use of a spectral index and the seasonal differencing method. However, the seasonal differencing has many errors when applied to a long-term time series. This article aims to develop and test two methods as an alternative to the traditional seasonal difference. The study area is the Chapada dos Veadeiros National Park (Central Brazil) that comprises different vegetation of the Cerrado biome. We used the MODIS/Terra Surface Reflectance 8-Day composite data, considering a 12-year period. The normalized burn ratio was calculated from the band 2 (250-meter resolution) and the band 7 (500-meter resolution reasampled to 250-meter). In this context, the normalization methods aim to eliminate all possible sources of spectral variation and highlight the burned-area features. The proposed normalization methods were the standardized time-series and the interannual phenological deviation. The standardized time-series calculate for each pixel the z-scores of its temporal curve, obtaining a mean of 0 and a standard deviation of 1. The second method establishes a reference curve for each pixel from the average interannual phenology that is subtracted for every year of its respective time series. Optimal threshold value between burned and unburned area for each method was determined from accuracy assessment curves, which compare different threshold values and its accuracy indices with a reference classification using Landsat TM. The different methods have similar accuracy for the burning event, where the standardized method has slightly better results. However, the seasonal difference method has a very false positive error, especially in the period between the rainy and dry seasons. The interannual phenological deviation method minimizes false positive errors, but some remain. In contrast, the standardized time series shows excellent results not containing this type of error. This precision is due to the design method that does not perform a subtraction with a baseline (prior year or average phenological curve). Thus, this method allows a high stability and can be implemented for the automatic detection of burned areas using long-term time series. Remote Sensing 2015-05-29 7 6 Article 10.3390/rs70606950 6950 6985 2072-4292 2015-05-29 doi: 10.3390/rs70606950 Osmar de Carvalho Júnior Renato Guimarães Cristiano Silva Roberto Gomes http://www.mdpi.com/2072-4292/7/6/6932 Diverse scene stitching is a challenging task in aerial video surveillance. This paper presents a hybrid stitching method based on the observation that aerial videos captured in real surveillance settings are neither totally ordered nor completely unordered. Often, human operators apply continuous monitoring of the drone to revisit the same area of interest. This monitoring mechanism yields to multiple short, successive video clips that overlap in either time or space. We exploit this property and treat the aerial image stitching problem as temporal sequential grouping and spatial cross-group retrieval. We develop an effective graph-based framework that can robustly conduct the grouping, retrieval and stitching tasks. To evaluate the proposed approach, we experiment on the large-scale VIRATaerial surveillance dataset, which is challenging for its heterogeneity in image quality and diversity of the scene. Quantitative and qualitative comparisons with state-of-the-art algorithms show the efficiency and robustness of our technique. Remote Sensing 2015-05-28 7 6 Article 10.3390/rs70606932 6932 6949 2072-4292 2015-05-28 doi: 10.3390/rs70606932 Tao Yang Jing Li Jingyi Yu Sibing Wang Yanning Zhang http://www.mdpi.com/2072-4292/7/6/6908 In this study deformation processes in northern Zona Metropolitana del Valle de Mexico (ZMVM) are evaluated by means of advanced multi-temporal interferometry. ERS and ENVISAT time series, covering approximately an 11-year period (between 1999 and 2010), were produced showing mainly linear subsidence behaviour for almost the entire area under study, but increasing rates that reach up to 285 mm/yr. Important non-linear deformation was identified in certain areas, presumably suggesting interaction between subsidence and other processes. Thus, a methodology for identification of probable fracturing zones based on discrimination and modelling of the non-linear (quadratic function) component is presented. This component was mapped and temporal subsidence evolution profiles were constructed across areas where notable acceleration (maximum of 8 mm/yr2) or deceleration (maximum of −9 mm/yr2) is found. This methodology enables location of potential soil fractures that could impact relevant infrastructure such as the Tunel Emisor Oriente (TEO) (along the structure rates exceed 200 mm/yr). Additionally, subsidence behaviour during wet and dry seasons is tackled in partially urbanized areas. This paper provides useful information for geological risk assessment in the area. Remote Sensing 2015-05-28 7 6 Article 10.3390/rs70606908 6908 6931 2072-4292 2015-05-28 doi: 10.3390/rs70606908 Gabriela Siles Juan Alcérreca-Huerta Penélope López-Quiroz Wolfgang Niemeier http://www.mdpi.com/2072-4292/7/6/6886 The adjacency effect and non-uniform responses complicate the precise delimitation of the surface support of remote sensing data and their derived products. Thus, modeling spatial response characteristics (SRCs) prior to using remote sensing information has become important. A point spread function (PSF) is typically used to describe the SRCs of the observation cells from remote sensors and is always estimated in a laboratory before the sensor is launched. However, research on the SRCs of high-order remote sensing products derived from the observations remains insufficient, which is an obstacle to converting between multi-scale remote sensing products and validating coarse-resolution products. This study proposed a method that combines simulation and validation to establish SRC models of coarse-resolution albedo products. Two series of commonly used 500-m/1-km resolution albedo products, which are derived from Moderate Resolution Imaging Spectroradiometer (MODIS) reflectance data, were investigated using 30-m albedo products that provide the required sub-pixel information. The analysis proves that the size of the surface support of each albedo pixel is larger than the nominal resolution of the pixel and that the response weight is non-uniformly distributed, with an elliptical Gaussian shape. The proposed methodology is generic and applicable for analyzing the SRCs of other advanced remote sensing products. Remote Sensing 2015-05-28 7 6 Article 10.3390/rs70606886 6886 6907 2072-4292 2015-05-28 doi: 10.3390/rs70606886 Jingjing Peng Qiang Liu Lizhao Wang Qinhuo Liu Wenjie Fan Meng Lu Jianguang Wen http://www.mdpi.com/2072-4292/7/6/6862 The primary restriction on high resolution remote sensing data is the limit observation frequency. Using a network of multiple sensors is an efficient approach to increase the observations in a specific period. This study explores a leaf area index (LAI) inversion method based on a 30 m multi-sensor dataset generated from HJ1/CCD and Landsat8/OLI, from June to August 2013 in the middle reach of the Heihe River Basin, China. The characteristics of the multi-sensor dataset, including the percentage of valid observations, the distribution of observation angles and the variation between different sensor observations, were analyzed. To reduce the possible discrepancy between different satellite sensors on LAI inversion, a quality control system for the observations was designed. LAI is retrieved from the high quality of single-sensor observations based on a look-up table constructed by a unified model. The averaged LAI inversion over a 10-day period is set as the synthetic LAI value. The percentage of valid LAI inversions increases significantly from 6.4% to 49.7% for single-sensors to 75.9% for multi-sensors. LAI retrieved from the multi-sensor dataset show good agreement with the field measurements. The correlation coefficient (R2) is 0.90, and the average root mean square error (RMSE) is 0.42. The network of multiple sensors with 30 m spatial resolution can generate LAI products with reasonable accuracy and meaningful temporal resolution. Remote Sensing 2015-05-28 7 6 Article 10.3390/rs70606862 6862 6885 2072-4292 2015-05-28 doi: 10.3390/rs70606862 Jing Zhao Jing Li Qinhuo Liu Wenjie Fan Bo Zhong Shanlong Wu Le Yang Yelu Zeng Baodong Xu Gaofei Yin http://www.mdpi.com/2072-4292/7/6/6828 Component substitution (CS) and multi-resolution analysis (MRA) are the two basic categories in the extended general image fusion (EGIF) framework for fusing panchromatic (Pan) and multispectral (MS) images. Despite of the method diversity, there are some unaddressed questions and contradictory conclusions about fusion. For example, is the spatial enhancement of CS methods better than MRA methods? Is spatial enhancement and spectral preservation competitive? How to achieve spectral consistency defined by Wald et al. in 1997? In their definition any synthetic image should be as identical as possible to the original image once degraded to its original resolution. To answer these questions, this research first finds out that all the CS and MRA methods can be derived from the Bayesian fusion method by adjusting a weight parameter to balance contributions from the spatial injection and spectral preservation models. The spectral preservation model assumes a Gaussian distribution of the desired high-resolution MS images, with the up-sampled low-resolution MS images comprising the mean value. The spatial injection model assumes a linear correlation between Pan and MS images. Thus the spatial enhancement depends on the weight parameter but is irrelevant of which category (i.e., MRA or CS) the method belongs to. This paper then adds a spectral consistency model in the Bayesian fusion framework to guarantee Wald's spectral consistency with regard to arbitrary sensor point spread function. Although the spectral preservation in the EGIF methods is competitive to spatial enhancement, the Wald's spectral consistency property is complementary with spatial enhancement. We conducted experiments on satellite images acquired by the QuickBird and WorldView-2 satellites to confirm our analysis, and found that the performance of the traditional EGIF methods improved significantly after adding the spectral consistency model. Remote Sensing 2015-05-28 7 6 Article 10.3390/rs70606828 6828 6861 2072-4292 2015-05-28 doi: 10.3390/rs70606828 Hankui Zhang Bo Huang http://www.mdpi.com/2072-4292/7/6/6808 The 2008 Wenchuan earthquake destroyed large areas of vegetation in the Baisha River and Longxi River basins, in Dujiangyan County, China. There were several debris flow events in these mountainous river basins after 2008. Currently, these damaged vegetation areas are in various stages of recovery. This recovery vegetation improves the resistance of slopes to both surficial erosion and mass wasting. We introduce a probabilistic approach to determining the relationships between damaged vegetation and slope materials' stability, and model the sediment and flow (hydrological) connectivity index to detect the hydrological changes in a given river basin, using the multi-temporal (1994–2014) remote-sensing images to monitor the vegetation recovery processes. Our results demonstrated that the earthquake-damaged vegetation areas have coupling relationships with topographic environment and slope material properties, and can be used to assess the slope material stability. Further, our analysis results showed that the areas with horizontal distance to river streams <500 m are areas that actively contribute sediment to the stream channel network, and are main material sources for debris flow processes in one given mountainous basin. Remote Sensing 2015-05-28 7 6 Article 10.3390/rs70606808 6808 6827 2072-4292 2015-05-28 doi: 10.3390/rs70606808 Huaizhen Zhang Xiaomeng Wang Jianrong Fan Tianhe Chi Shun Yang Ling Peng http://www.mdpi.com/2072-4292/7/6/6784 A land-cover-based linear BRDF (bi-directional reflectance distribution function) unmixing (LLBU) algorithm based on the kernel-driven model is proposed to combine the compact airborne spectrographic imager (CASI) reflectance with the moderate resolution imaging spectroradiometer (MODIS) daily reflectance product to derive the BRDF/albedo of the two sensors simultaneously in the foci experimental area (FEA) of the Heihe Watershed Allied Telemetry Experimental Research (HiWATER), which was carried out in the Heihe River basin, China. For each land cover type, an archetypal BRDF, which characterizes the shape of its anisotropic reflectance, is extracted by linearly unmixing from the MODIS reflectance with the assistance of a high-resolution classification map. The isotropic coefficients accounting for the differences within a class are derived from the CASI reflectance. The BRDF is finally determined by the archetypal BRDF and the corresponding isotropic coefficients. Direct comparisons of the cropland archetypal BRDF and CASI albedo with in situ measurements show good agreement. An indirect validation which compares retrieved BRDF/albedo with that of the MCD43A1 standard product issued by NASA and aggregated CASI albedo also suggests reasonable reliability. LLBU has potential to retrieve the high spatial resolution BRDF/albedo product for airborne and spaceborne sensors which have inadequate angular samplings. In addition, it can shorten the timescale for coarse spatial resolution product like MODIS. Remote Sensing 2015-05-28 7 6 Article 10.3390/rs70606784 6784 6807 2072-4292 2015-05-28 doi: 10.3390/rs70606784 Dongqin You Jianguang Wen Qing Xiao Qiang Liu Qinhuo Liu Yong Tang Baocheng Dou Jingjing Peng http://www.mdpi.com/2072-4292/7/6/6763 Acquiring multi-temporal spatial information on vegetation condition at scales appropriate for site-specific agricultural management is often complicated by the need for meticulous field measurements. Understanding spatial/temporal crop cover heterogeneity within irrigated croplands may support sustainable land use, specifically in areas affected by land degradation due to secondary soil salinization. This study demonstrates the use of multi-temporal, high spatial resolution (10 m) SPOT-4/5 image data in an integrated change vector analysis and spectral mixture analysis (CVA-SMA) procedure. This procedure was implemented with the principal objective of mapping sub-field vegetation cover dynamics in irrigated lowland areas within the lowerlands of the Amu Darya River. CVA intensity and direction were calculated separately for the periods of 1998–2006 and 2006–2010. Cumulative change intensity and the overall directional trend were also derived for the entire observation period of 1998–2010. Results show that most of the vector changes were observed between 1998 and 2006; persistent conditions were seen within the study region during the 2006–2010 period. A decreasing vegetation cover trend was identified within 38% of arable land. Areas of decreasing vegetation cover were located principally in the irrigation system periphery where deficient water supply and low soil quality lead to substandard crop development. During the 2006–2010 timeframe, degraded crop cover conditions persisted in 37% of arable land. Vegetation cover increased in 25% of the arable land where irrigation water supply was adequate. This high sub-field crop performance spatial heterogeneity clearly indicates that current land management practices are inefficient. Such information can provide the basis for implementing and adapting irrigation applications and salt leaching techniques to site-specific conditions and thereby make a significant contribution to sustainable regional land management. Remote Sensing 2015-05-26 7 6 Article 10.3390/rs70606763 6763 6783 2072-4292 2015-05-26 doi: 10.3390/rs70606763 Olena Dubovyk Gunter Menz Alexander Lee Juergen Schellberg Frank Thonfeld Asia Khamzina http://www.mdpi.com/2072-4292/7/6/6741 Snow is an important water resource and greatly influences water availability in the downstream areas. In this study, snow cover variations of the Upper Heihe River Basin (UHRB) during hydrological years (HY) 2003–2013 (September through August) is examined using the flexible multiday-combined MODIS snow cover products. Spatial distribution and pattern of snow cover from year to year for the basin is found to be relatively stable, with maximum snow cover area (SCA) and snow cover days occurring in HY2004, HY2008 and HY2012. A method, based on correlation coefficients between SCA and climate factors (mainly air temperature and precipitation), is presented to identify the threshold altitude that determines contributions of climate factors to SCA. A threshold altitude of 3650 ± 150 m is found for the UHRB, where below this altitude, both air temperature (Tair) and precipitation are negative factors on SCA, except in the winter season when both are positive factors. Above the threshold altitude, precipitation acts as a positive factor except in summer, while Tair is a negative factor except in autumn. Overall, Tair is the primary controlling factor on SCA below the threshold altitude, while precipitation is the primary controlling factor on SCA above the threshold altitude. Remote Sensing 2015-05-26 7 6 Article 10.3390/rs70606741 6741 6762 2072-4292 2015-05-26 doi: 10.3390/rs70606741 Yunbo Bi Hongjie Xie Chunlin Huang Changqing Ke http://www.mdpi.com/2072-4292/7/6/6710 Open geospatial data sources provide opportunities for low cost 3D scene reconstruction. In this study, based on a sparse airborne laser scanning (ALS) point cloud (0.8 points/m2) obtained from open source databases, a building reconstruction pipeline for CAD building models was developed. The pipeline includes voxel-based roof patch segmentation, extraction of the key-points representing the roof patch outline, step edge identification and adjustment, and CAD building model generation. The advantages of our method lie in generating CAD building models without the step of enforcing the edges to be parallel or building regularization. Furthermore, although it has been challenging to use sparse datasets for 3D building reconstruction, our result demonstrates the great potential in such applications. In this paper, we also investigated the applicability of open geospatial datasets for 3D road detection and reconstruction. Road central lines were acquired from an open source 2D topographic database. ALS data were utilized to obtain the height and width of the road. A constrained search method (CSM) was developed for road width detection. The CSM method was conducted by splitting a given road into patches according to height and direction criteria. The road edges were detected patch by patch. The road width was determined by the average distance from the edge points to the central line. As a result, 3D roads were reconstructed from ALS and a topographic database. Remote Sensing 2015-05-26 7 6 Article 10.3390/rs70606710 6710 6740 2072-4292 2015-05-26 doi: 10.3390/rs70606710 Lingli Zhu Matti Lehtomäki Juha Hyyppä Eetu Puttonen Anssi Krooks Hannu Hyyppä http://www.mdpi.com/2072-4292/7/6/6683 Tracking land surface dynamics over cloud prone areas with complex mountainous terrain is an important challenge facing the Earth Science community. One such region is the Lake Kivu region in Central Africa. We developed a processing chain to systematically monitor the spatio-temporal land use/land cover dynamics of this region over the years 1988, 2001, and 2011 using Landsat data, complemented by ancillary data. Topographic compensation was performed on Landsat reflectances to avoid the strong illumination angle impacts and image compositing was used to compensate for frequent cloud cover and thus incomplete annual data availability in the archive. A systematic supervised classification was applied to the composite Landsat imagery to obtain land cover thematic maps with overall accuracies of 90% and higher. Subsequent change analysis between these years found extensive conversions of the natural environment as a result of human related activities. The gross forest cover loss for 1988–2001 and 2001–2011 period was 216.4 and 130.5 thousand hectares, respectively, signifying significant deforestation in the period of civil war and a relatively stable and lower deforestation rate later, possibly due to conservation and reforestation efforts in the region. The other dominant land cover changes in the region were aggressive subsistence farming and urban expansion displacing natural vegetation and arable lands. Despite limited data availability, this study fills the gap of much needed detailed and updated land cover change information for this biologically important region of Central Africa. These multi-temporal datasets will be a valuable baseline for land use managers in the region interested in developing ecologically sustainable land management strategies and measuring the impacts of biodiversity conservation efforts. Remote Sensing 2015-05-26 7 6 Article 10.3390/rs70606683 6683 6709 2072-4292 2015-05-26 doi: 10.3390/rs70606683 Bikash Basnet Anthony Vodacek http://www.mdpi.com/2072-4292/7/6/6663 Accurate determination of the amount of incoming solar radiation at Earth's surface is important for both climate studies and solar power applications. Satellite-based datasets of solar radiation offer wide spatial and temporal coverage, but careful validation of their quality is a necessary prerequisite for reliable utilization. Here we study the retrieval quality of one polar-orbiting satellite-based dataset (CLARA-A1) and one geostationary satellite-based dataset (SARAH), using in situ observations of solar radiation from the Finnish and Swedish meteorological measurement networks as reference. Our focus is on determining dataset quality over high latitudes as well as evaluating daily mean retrievals, both of which are aspects that have drawn little focus in previous studies. We find that both datasets are generally capable of retrieving the levels and seasonal cycles of solar radiation in Finland and Sweden well, with some limitations. SARAH exhibits a slight negative bias and increased retrieval uncertainty near the coverage edge, but in turn offers better precision (less scatter) in the daily mean retrievals owing to the high sampling rate of geostationary imaging. Remote Sensing 2015-05-26 7 6 Article 10.3390/rs70606663 6663 6682 2072-4292 2015-05-26 doi: 10.3390/rs70606663 Aku Riihelä Thomas Carlund Jörg Trentmann Richard Müller Anders Lindfors http://www.mdpi.com/2072-4292/7/6/6635 This paper presents a practical framework for the integration of unmanned aerial vehicle (UAV) based photogrammetry and terrestrial laser scanning (TLS) with application to open-pit mine areas, which includes UAV image and TLS point cloud acquisition, image and cloud point processing and integration, object-oriented classification and three-dimensional (3D) mapping and monitoring of open-pit mine areas. The proposed framework was tested in three open-pit mine areas in southwestern China. (1) With respect to extracting the conjugate points of the stereo pair of UAV images and those points between TLS point clouds and UAV images, some feature points were first extracted by the scale-invariant feature transform (SIFT) operator and the outliers were identified and therefore eliminated by the RANdom SAmple Consensus (RANSAC) approach; (2) With respect to improving the accuracy of geo-positioning based on UAV imagery, the ground control points (GCPs) surveyed from global positioning systems (GPS) and the feature points extracted from TLS were integrated in the bundle adjustment, and three scenarios were designed and compared; (3) With respect to monitoring and mapping the mine areas for land reclamation, an object-based image analysis approach was used for the classification of the accuracy improved UAV ortho-image. The experimental results show that by introduction of TLS derived point clouds as GCPs, the accuracy of geo-positioning based on UAV imagery can be improved. At the same time, the accuracy of geo-positioning based on GCPs form the TLS derived point clouds is close to that based on GCPs from the GPS survey. The results also show that the TLS derived point clouds can be used as GCPs in areas such as in mountainous or high-risk environments where it is difficult to conduct a GPS survey. The proposed framework achieved a decimeter-level accuracy for the generated digital surface model (DSM) and digital orthophoto map (DOM), and an overall accuracy of 90.67% for classification of the land covers in the open-pit mine. Remote Sensing 2015-05-26 7 6 Article 10.3390/rs70606635 6635 6662 2072-4292 2015-05-26 doi: 10.3390/rs70606635 Xiaohua Tong Xiangfeng Liu Peng Chen Shijie Liu Kuifeng Luan Lingyun Li Shuang Liu Xianglei Liu Huan Xie Yanmin Jin Zhonghua Hong http://www.mdpi.com/2072-4292/7/6/6611 Supervised target detection and anomaly detection are widely used in various applications, depending upon the availability of target spectral signature. Basically, they are based on a similar linear process, which makes them highly correlated. In this paper, we propose a novel adjusted spectral matched filter (ASMF) for hyperspectral target detection, which aims to effectively improve target detection performance with anomaly detection output. Specifically, a typical case is presented by using the Reed-Xiaoli (RX) anomaly detector to adjust the output of supervised constrained energy minimization (CEM) detector. The adjustment is appropriately controlled by a weighting parameter in different detection scenarios. Experiments were implemented by using both synthetic and real hyperspectral datasets. Compared to the traditional single detection method (e.g., CEM), the experimental results demonstrate that the proposed ASMF can effectively improve its performance by utilizing the result from an anomaly detector (e.g., RX), particularly in situations with a complex background or strong anomalies. Remote Sensing 2015-05-26 7 6 Article 10.3390/rs70606611 6611 6634 2072-4292 2015-05-26 doi: 10.3390/rs70606611 Lianru Gao Bin Yang Qian Du Bing Zhang http://www.mdpi.com/2072-4292/7/6/6576 Coal fires that are induced by natural spontaneous combustion or result from human activities occurring on the surface and in underground coal seams destroy coal resources and cause serious environmental degradation. Thermal infrared image data, which directly measure surface temperature, can be an important tool to map coal fires over large areas. As the first of two parts introducing our coal fire detection method, this paper proposes a self-adaptive threshold-based approach for coal fire detection using ASTER thermal infrared data: the self-adaptive gradient-based thresholding method (SAGBT). This method is based on an assumption that the attenuation of temperature along the coal fire's boundaries generates considerable numbers of spots with extremely high gradient values. The SAGBT method applied mathematical morphology thinning to skeletonize the potential high gradient buffers into the extremely high gradient lines, which provides a self-adaptive mechanism to generate thresholds according to the thermal spatial patterns of the images. The final threshold was defined as an average temperature value reading from the high temperature buffers (segmented by 1.0 σ from the mean) and along a sequence of extremely high gradient lines (thinned from the potential high gradient buffers and segmented within the lower bounds, ranging from 0.5 σ to 1.5 σ and with an upper bound of 3.2 σ, where σ is the standard deviation), marking the coal fire areas. The SAGBT method used the basic outer boundary of the coal-bearing strata to simply exclude false alarms. The intermediate thresholds reduced the coupling with the temperature and converged by changing the potential high gradient buffers. This simple approach can be economical and accurate in identifying coal fire areas. In addition, it allows for the identification of thresholds using multiple ASTER TIR scenes in a consistent and uniform manner, and supports long-term coal fire change analyses using historical images in local areas. This paper focuses on the introduction of the methodology. Furthermore, an improvement to SAGBT is proposed. In a subsequent paper, subtitled "Part 2, Validation and Sensitivity Analysis," we address satellite-field simultaneous observations and report comparisons between the retrieved thermal anomalies and field measurements in different aspects to prove that the coal fires are separable by the SAGBT method. These comparisons allowed us to estimate the accuracy and biases of the SAGBT method. As an application of the SAGBT, a relationship between coal fires' decadal variation and coal production was also examined. Our work documented a total area increase in the beginning of 2003, which correlates with increased mining activities and the rapid increase of energy consumption in China during the decade (2001–2011). Additionally, a decrease in the total coal fire area is consistent with the nationally sponsored fire suppression efforts during 2007–2008. It demonstrated the applicability of SAGBT method for long-term change detection with multi-temporal images. Remote Sensing 2015-05-26 7 6 Article 10.3390/rs70606576 6576 6610 2072-4292 2015-05-26 doi: 10.3390/rs70606576 Xiaomin Du Daiyong Cao Deepak Mishra Sergio Bernardes Thomas Jordan Marguerite Madden http://www.mdpi.com/2072-4292/7/6/6558 In rugged terrain, the accuracy of surface reflectance estimations is compromised by atmospheric and topographic effects. We propose a new method to simultaneously eliminate atmospheric and terrain effects in Landsat Thematic Mapper (TM) images based on a 30 m digital elevation model (DEM) and Moderate Resolution Imaging Spectroradiometer (MODIS) atmospheric products. Moreover, we define a normalized factor of a Bidirectional Reflectance Distribution Function (BRDF) to convert the sloping pixel reflectance into a flat pixel reflectance by using the Ross Thick-Li Sparse BRDF model (Ambrals algorithm) and MODIS BRDF/albedo kernel coefficient products. Sole atmospheric correction and topographic normalization were performed for TM images in the upper stream of the Heihe River Basin. The results show that using MODIS atmospheric products can effectively remove atmospheric effects compared with the Fast Line-of-Sight Atmospheric Analysis of Spectral Hypercubes (FLAASH) model and the Landsat Climate Data Record (CDR). Moreover, superior topographic effect removal can be achieved by considering the surface BRDF when compared with the surface Lambertian assumption of topographic normalization. Remote Sensing 2015-05-26 7 6 Article 10.3390/rs70606558 6558 6575 2072-4292 2015-05-26 doi: 10.3390/rs70606558 Yanli Zhang Xin Li Jianguang Wen Qinhuo Liu Guangjian Yan http://www.mdpi.com/2072-4292/7/6/6535 Asian countries strongly depend on rice production for food security. The major rice-growing season (June to October) is highly exposed to the risk of tropical storm related damage. Unbiased and transparent approaches to assess the risk of rice crop damage are essential to support mitigation and disaster response strategies in the region. This study describes and demonstrates a method for rapid, pre-event crop status assessment. The ex-post test case is Typhoon Haiyan and its impact on the rice crop in Leyte Province in the Philippines. A synthetic aperture radar (SAR) derived rice area map was used to delineate the area at risk while crop status at the moment of typhoon landfall was estimated from specific time series analysis of Moderate Resolution Imaging Spectroradiometer (MODIS) data. A spatially explicit indicator of risk of standing crop loss was calculated as the time between estimated heading date and typhoon occurrence. Results of the analysis of pre- and post-event SAR images showed that 6500 ha were flooded in northeastern Leyte. This area was also the region most at risk to storm related crop damage due to late establishment of rice. Estimates highlight that about 700 ha of rice (71% of which was in northeastern Leyte) had not reached maturity at the time of the typhoon event and a further 8400 ha (84% of which was in northeastern Leyte) were likely to be not yet harvested. We demonstrated that the proposed approach can provide pre-event, in-season information on the status of rice and other field crops and the risk of damage posed by tropical storms. Remote Sensing 2015-05-26 7 6 Article 10.3390/rs70606535 6535 6557 2072-4292 2015-05-26 doi: 10.3390/rs70606535 Mirco Boschetti Andrew Nelson Francesco Nutini Giacinto Manfron Lorenzo Busetto Massimo Barbieri Alice Laborte Jeny Raviz Francesco Holecz Mary Mabalay Alfie Bacong Eduardo Quilang http://www.mdpi.com/2072-4292/7/6/6510 Image time series of high temporal and spatial resolution capture land surface dynamics of heterogeneous landscapes. We applied the ESTARFM (Enhanced Spatial and Temporal Adaptive Reflectance Fusion Model) algorithm to multi-spectral images covering two semi-arid heterogeneous rangeland study sites located in South Africa. MODIS 250 m resolution and RapidEye 5 m resolution images were fused to produce synthetic RapidEye images, from June 2011 to July 2012. We evaluated the performance of the algorithm by comparing predicted surface reflectance values to real RapidEye images. Our results show that ESTARFM predictions are accurate, with a coefficient of determination for the red band 0.80 < R2 < 0.92, and for the near-infrared band 0.83 < R2 < 0.93, a mean relative bias between 6% and 12% for the red band and 4% to 9% in the near-infrared band. Heterogeneous vegetation at sub-MODIS resolution is captured adequately: A comparison of NDVI time series derived from RapidEye and ESTARFM data shows that the characteristic phenological dynamics of different vegetation types are reproduced well. We conclude that the ESTARFM algorithm allows us to produce synthetic remote sensing images at high spatial combined with high temporal resolution and so provides valuable information on vegetation dynamics in semi-arid, heterogeneous rangeland landscapes. Remote Sensing 2015-05-26 7 6 Article 10.3390/rs70606510 6510 6534 2072-4292 2015-05-26 doi: 10.3390/rs70606510 Andreas Tewes Frank Thonfeld Michael Schmidt Roelof Oomen Xiaolin Zhu Olena Dubovyk Gunter Menz Jürgen Schellberg http://www.mdpi.com/2072-4292/7/5/6489 This study analyzed the scaling problem of land surface temperature (LST) data retrieved with the Temperature Emissivity Separation (TES) algorithm. We compiled a remotely sensed dataset that included Thermal Airborne Hyperspectral Imager (TASI) and satellite-based Advanced Spaceborne Thermal Emission Reflection (ASTER) data, which were acquired simultaneously. This dataset provided the range of spatial heterogeneities of land surface necessary for the study, which was quantified by the dispersion variance. The LST scaling problem was studied by comparing the remotely sensed LST products in two ways. First, the LST products calculated in the distributed method and the lumped method were compared. Second, the airborne and satellite-based LST products derived from the TES algorithm were compared. Four upscaling methods of LST were used in the process. A scaling correction methodology was developed based on the comparisons. The results showed that the scaling effect could be as large as 0.8 when the spatial resolution of the TASI LST data was coarse. The scaling effect increases quickly with the spatial resolution until it reaches the characteristic scale of the landscape and is positively correlated with the spatial heterogeneity. The first two upscaling methods denoted as Methods 1–2 can upscale the LST more effectively when compared with the other two scaling methods (Methods 3–4). The scaling effect for the ASTER data is not notable. The comparison between the TASI and ASTER data showed that they were highly consistent, with a root mean square error (RMSE) of approximately 0.88 K, when the pixels were relatively homogeneous. When the spatial heterogeneity was significant, the RMSE was as large as 2.68 K The scaling correction methodology provided resolution-invariant results with scaling effects of less than 0.5 K. Remote Sensing 2015-05-22 7 5 Article 10.3390/rs70506489 6489 6509 2072-4292 2015-05-22 doi: 10.3390/rs70506489 Tian Hu Qinhuo Liu Yongming Du Hua Li Heshun Wang Biao Cao http://www.mdpi.com/2072-4292/7/5/6454 Mean Annual Precipitation is one of the most important variables used in water resource management. However, quantifying Mean Annual Precipitation at high spatial resolution, needed for advanced hydrological analysis, is challenging in developing countries which often present a sparse gauge network and a highly variable climate. In this work, we present a methodology to quantify Mean Annual Precipitation at 1 km spatial resolution using different precipitation products from satellite estimates and gauge observations at coarse spatial resolution (i.e., ranging from 4 km to 25 km). Examples of this methodology are given for South America and West Africa. We develop a downscaling method that exploits the relationship among satellite-derived rainfall, Digital Elevation Model and Enhanced Vegetation Index. Finally, we validate its performance using rain gauge measurements: comparable annual precipitation estimates for both South America and West Africa are retrieved. Validation indicates that high resolution Mean Annual Precipitation downscaled from CHIRP (Climate Hazards Group Infrared Precipitation) and GPCC (Global Precipitation Climatology Centre) datasets present the best ensemble of performance statistics for both South America and West Africa. Results also highlight the potential of the presented technique to downscale satellite-derived rainfall worldwide. Remote Sensing 2015-05-22 7 5 Article 10.3390/rs70506454 6454 6488 2072-4292 2015-05-22 doi: 10.3390/rs70506454 Guido Ceccherini Iban Ameztoy Claudia Hernández Cesar Moreno http://www.mdpi.com/2072-4292/7/5/6433 This study presents a revised temporal scaling method based on a detection algorithm for the temporal stability of the evaporative fraction (EF) to estimate total daytime evapotranspiration (ET) at a regional scale. The study area is located in the Heihe River Basin, which is the second largest inland river basin in China. The remote sensing data and field observations used in this study were obtained from the Heihe Watershed Allied Telemetry Experimental Research (HiWATER) project. The half-hourly EF values (EFEC) calculated using meteorological observations from an eddy covariance (EC) system and an automatic meteorological station (AMS) represented the diurnal pattern of the EF across the majority of the study area. The remotely sensed instantaneous midday EF (EFASTER), which indicates the spatial distribution of the midday EF over the entire study area, was calculated from an Advanced Spaceborne Thermal Emission and Reflection Radiometer (ASTER) image. The temporal stability of EFEC was examined using a detection algorithm. Intervals with inconsistent EFEC values were distinguished from those with consistent EFEC values; the total daytime ET (from 9:00 to 19:00) within these interval types was integrated separately. Validation of the total daytime ET at the satellite pixel scale was conducted using measurements from17 EC towers. Using the detection algorithm for the temporal stability of the EF and dynamic adjustment, the revised temporal scaling method resulted in a root-mean-square error (RMSE) of 0.54 (mm·d−1), a mean relative error (MRE) of 7.26% and a correlation coefficient (Corr.) of 0.81; all of these values were superior to those of the two other methods (i.e., the constant EF and variable EF methods). The revised method easily extends to other areas and exhibits a superior performance in flat and regularly-irrigated farmlands at the regional scale. Remote Sensing 2015-05-21 7 5 Article 10.3390/rs70506433 6433 6453 2072-4292 2015-05-21 doi: 10.3390/rs70506433 Yi Song Mingguo Ma Long Jin Xufeng Wang http://www.mdpi.com/2072-4292/7/5/6414 It remains a challenging issue to accurately estimate the fraction of absorbed photosynthetically-active radiation (FPAR) using remote sensing data, as the direct and diffuse radiation reaching the vegetation canopy have different effects for FPAR. In this research, a FPAR inversion model was developed that may distinguish direct and diffuse radiation (the DnD model) based on the energy budget balance principle. Taking different solar zenith angles and diffuse PAR proportions as inputs, the instantaneous FPAR could be calculated. As the leaf area index (LAI) and surface albedo do not vary in a short periods, the FPAR not only on a clear day, but also on a cloudy day may be calculated. This new method was used to produce the FPAR products in the Heihe River Basin with the Moderate-Resolution Imaging Spectroradiometer (MODIS) LAI and surface albedo products as the input data source. The instantaneous FPAR was validated by using field-measured data (RMSE is 0.03, R2 is 0.85). The daily average FPAR was compared with the MODIS FPAR product. The inversion results and the MODIS FPAR product are highly correlated, but the MODIS FPAR product is slightly high in forest areas, which is in agreement with other studies for MODIS FPAR products. Remote Sensing 2015-05-21 7 5 Article 10.3390/rs70506414 6414 6432 2072-4292 2015-05-21 doi: 10.3390/rs70506414 Li Li Yongming Du Yong Tang Xiaozhou Xin Hailong Zhang Jianguang Wen Qinhuo Liu http://www.mdpi.com/2072-4292/7/5/6380 The applications of object-based image analysis (OBIA) in remote sensing studies of wetlands have been growing over recent decades, addressing tasks from detection and delineation of wetland bodies to comprehensive analyses of within-wetland cover types and their change. Compared to pixel-based approaches, OBIA offers several important benefits to wetland analyses related to smoothing of the local noise, incorporating meaningful non-spectral features for class separation and accounting for landscape hierarchy of wetland ecosystem organization and structure. However, there has been little discussion on whether unique challenges of wetland environments can be uniformly addressed by OBIA across different types of data, spatial scales and research objectives, and to what extent technical and conceptual aspects of this framework may themselves present challenges in a complex wetland setting. This review presents a synthesis of 73 studies that applied OBIA to different types of remote sensing data, spatial scale and research objectives. It summarizes the progress and scope of OBIA uses in wetlands, key benefits of this approach, factors related to accuracy and uncertainty in its applications and the main research needs and directions to expand the OBIA capacity in the future wetland studies. Growing demands for higher-accuracy wetland characterization at both regional and local scales together with advances in very high resolution remote sensing and novel tasks in wetland restoration monitoring will likely continue active exploration of the OBIA potential in these diverse and complex environments. Remote Sensing 2015-05-21 7 5 Review 10.3390/rs70506380 6380 6413 2072-4292 2015-05-21 doi: 10.3390/rs70506380 Iryna Dronova http://www.mdpi.com/2072-4292/7/5/6358 The paper reports the recent progress in the radiative transfer model (RTM) development, which serves as the observation operator of a Land Data Assimilation System (LDAS), and its validation at two Planetary Boundary Layer (PBL) stations with different weather and land cover conditions: Wenjiang station of humid and cropped field and Gaize station of arid and bare soil field. In situ observed micrometeorological data were used as the driven data of LDAS, in which AMSR-E brightness temperatures (TB) were assimilated into a land surface model (LSM). Near surface soil moisture content output from LDAS, together with the one simulated by a LSM with default parameters, were compared to the in-situ soil moisture observation. The comparison results successfully validated the capability of LDAS with new RTM to simulate near surface soil moisture at various environments, supporting that LDAS can generally simulate soil moisture with a reasonable accuracy for both humid vegetated fields and arid bare soil fields while the LSM overestimates near surface soil moisture for humid vegetated fields and underestimates soil moisture for arid bare soil fields. Remote Sensing 2015-05-21 7 5 Article 10.3390/rs70506358 6358 6379 2072-4292 2015-05-21 doi: 10.3390/rs70506358 Hui Lu Kun Yang Toshio Koike Long Zhao Jun Qin http://www.mdpi.com/2072-4292/7/5/6336 Lidar data provide both geometric and radiometric information. Radiometric information is influenced by sensor and target factors and should be calibrated to obtain consistent energy responses. The radiometric correction of airborne lidar system (ALS) converts the amplitude into a backscatter cross-section with physical meaning value by applying a model-driven approach. The radiometric correction of terrestrial mobile lidar system (MLS) is a challenging task because it does not completely follow the inverse square range function at near-range. This study proposed a radiometric normalization workflow for MLS using a data-driven approach. The scope of this study is to normalize amplitude of road points for road surface classification, assuming that road points from different scanners or strips should have similar responses in overlapped areas. The normalization parameters for range effect were obtained from crossroads. The experiment showed that the amplitude difference between scanners and strips decreased after radiometric normalization and improved the accuracy of road surface classification. Remote Sensing 2015-05-21 7 5 Article 10.3390/rs70506336 6336 6357 2072-4292 2015-05-21 doi: 10.3390/rs70506336 Tee-Ann Teo Hui-Lin Yu http://www.mdpi.com/2072-4292/7/5/6320 In order to investigate the aerosol variability over the southwest region of Slovenia, an ultraviolet scanning elastic backscatter LiDAR was utilized to make the vertical scan for atmospheric probing. With the assumption of horizontal atmospheric homogeneity, aerosol optical variables were retrieved from the horizontal pixel data points of two-dimensional range-height-indicator (RHI) diagrams by using a multiangle retrieval method, in which optical depth is defined as the slope of the resulting linear function when height is kept constant. To make the data retrieval feasible and precise, a series of key procedures complemented the data processing, including construction of the RHI diagram, correction of Rayleigh scattering, assessment of horizontal atmospheric homogeneity and retrieval of aerosol optical variables. The measurement example demonstrated the feasibility of the ultraviolet scanning elastic backscatter LiDAR in the applications of the retrieval of aerosol extinction and determination of the atmospheric boundary layer height. Three months' data combined with the modeling of air flow trajectories using Hybrid Single Particle Lagrangian Integrated Trajectory Model were analyzed to investigate aerosol variability. The average value of aerosol extinction with the presence of land-based air masses from the European continent was found to be two-times larger than that influenced by marine aerosols from the Mediterranean or Adriatic Sea. Remote Sensing 2015-05-20 7 5 Article 10.3390/rs70506320 6320 6335 2072-4292 2015-05-20 doi: 10.3390/rs70506320 Fei Gao Samo Stanič Klemen Bergant Ying Li Songhui Li Dengxin Hua Longlong Wang http://www.mdpi.com/2072-4292/7/5/6296 Optical remotely sensed images in mountainous areas are subject to radiometric distortions induced by topographic effects, which need to be corrected before quantitative applications. Based on Li model and Sandmeier model, this paper proposed an improved physics-based model for the topographic correction of Landsat Thematic Mapper (TM) images. The model employed Normalized Difference Vegetation Index (NDVI) thresholds to approximately divide land targets into eleven groups, due to NDVI's lower sensitivity to topography and its significant role in indicating land cover type. Within each group of terrestrial targets, corresponding MODIS BRDF (Bidirectional Reflectance Distribution Function) products were used to account for land surface's BRDF effect, and topographic effects are corrected without Lambertian assumption. The methodology was tested with two TM scenes of severely rugged mountain areas acquired under different sun elevation angles. Results demonstrated that reflectance of sun-averted slopes was evidently enhanced, and the overall quality of images was improved with topographic effect being effectively suppressed. Correlation coefficients between Near Infra-Red band reflectance and illumination condition reduced almost to zero, and coefficients of variance also showed some reduction. By comparison with the other two physics-based models (Sandmeier model and Li model), the proposed model showed favorable results on two tested Landsat scenes. With the almost half-century accumulation of Landsat data and the successive launch and operation of Landsat 8, the improved model in this paper can be potentially helpful for the topographic correction of Landsat and Landsat-like data. Remote Sensing 2015-05-20 7 5 Article 10.3390/rs70506296 6296 6319 2072-4292 2015-05-20 doi: 10.3390/rs70506296 Ainong Li Qingfang Wang Jinhu Bian Guangbin Lei http://www.mdpi.com/2072-4292/7/5/6280 Endmember selection is the basis for sub-pixel land cover classifications using multiple endmember spectral mixture analysis (MESMA) that adopts variant endmember matrices for each pixel to mitigate errors caused by endmember variability in SMA. A spectral library covering a large number of endmembers can account for endmember variability, but it also lowers the computational efficiency. Therefore, an efficient endmember selection scheme to optimize the library is crucial to implement MESMA. In this study, we present an endmember selection method based on vector length. The spectra of a land cover class were divided into subsets using vector length intervals of the spectra, and the representative endmembers were derived from these subsets. Compared with the available endmember average RMSE (EAR) method, our approach improved the computational efficiency in endmember selection. The method accuracy was further evaluated using spectral libraries derived from the ground reference polygon and Moderate Resolution Imaging Spectroradiometer (MODIS) imagery respectively. Results using the different spectral libraries indicated that MESMA combined with the new approach performed slightly better than EAR method, with Kappa coefficient improved from 0.75 to 0.78. A MODIS image was used to test the mapping fraction, and the representative spectra based on vector length successfully modeled more than 90% spectra of the MODIS pixels by 2-endmember models. Remote Sensing 2015-05-20 7 5 Letter 10.3390/rs70506280 6280 6295 2072-4292 2015-05-20 doi: 10.3390/rs70506280 Yuanliu Xu Jiancheng Shi Jinyang Du http://www.mdpi.com/2072-4292/7/5/6257 In this paper, grid cell based spatial statistics were used to quantify the drivers of land-cover and land-use change (LCLUC) and habitat degradation in a tropical rainforest in Madagascar. First, a spectral database of various land-cover and land-use information was compiled using multi-year field campaign data and photointerpretation of satellite images. Next, residential areas were extracted from IKONOS-2 and GeoEye-1 images using object oriented feature extraction (OBIA). Then, Landsat Thematic Mapper (TM) and Enhanced Thematic Mapper Plus (ETM+) data were used to generate land-cover and land-use maps from 1990 to 2011, and LCLUC maps were developed with decadal intervals and converted to 100 m vector grid cells. Finally, the causal associations between LCLUC were quantified using ordinary least square regression analysis and Moran's I, and a forest disturbance index derived from the time series Landsat data were used to further confirm LCLUC drivers. The results showed that (1) local spatial statistical approaches were most effective at quantifying the drivers of LCLUC, and (2) the combined threats of habitat degradation in and around the reserve and increasing encroachment of invasive plant species lead to the expansion of shrubland and mixed forest within the former primary forest, which was echoed by the forest disturbance index derived from the Landsat data. Remote Sensing 2015-05-20 7 5 Article 10.3390/rs70506257 6257 6279 2072-4292 2015-05-20 doi: 10.3390/rs70506257 Abduwasit Ghulam Oghlan Ghulam Maitiniyazi Maimaitijiang Karen Freeman Ingrid Porton Matthew Maimaitiyiming