ALBERT

Description: Postdisaster search and rescue is an important application of ultrawideband (UWB) radar systems, which mainly detect trapped victims by their respiratory-motion response. The development of a respiratory-motion detection (RMD) algorithm that can eliminate nonstationary clutter and noise is a challenging task for the application. A new algorithm is proposed to deal with the task in this letter. It uses the multichannel singular spectrum analysis (MSSA) technique to reconstruct the respiratory-motion response detected by a UWB radar. During the reconstruction, the periodicity and range interrelation characteristics of the response are exploited to adaptively identify signal subspaces. The performance of the algorithm is verified both by simulated and real data. The results show its improved performance over the reference algorithms, e.g., a singular-value-decomposition-based algorithm. The adaptive-MSSA-based RMD algorithm has great promise not only in practical use but also for future research of UWB-radar-based human being remote sensing.

Description: Land–sea contamination observed in Soil Moisture and Ocean Salinity (SMOS) brightness temperature images is found to have two main contributions: the floor error inherent of image reconstruction and a multiplicative error either in the antenna temperature or in the visibility samples measured by the correlator. The origin of this last one is traced down to SMOS calibration parameters to yield a simple correction scheme, which is validated against several geophysical scenarios. Autoconsistency rules in interferometric synthesis together with redundant and complementary calibration procedures provide a robust SMOS calibration scheme.

Description: This letter proposes a novel algorithm, which is based on the generalized method of moments (GMM), for the estimation and correction of phase errors induced in synthetic aperture radar (SAR) imagery. The GMM algorithm is used to replace the original phase-estimation kernel in the basic structure of the phase-gradient-autofocus algorithm. Since this novel algorithm does not require the observed signal to be a certain distribution model, it is able to estimate arbitrary phase errors. The GMM algorithm has the ability of estimating range-dependent phase errors, which makes it an efficient estimator. As a result, higher accuracy of the estimated phase errors and a better focused image can be achieved. Excellent results have been obtained in autofocusing and imaging experiments on real SAR data.

Description: Strong clutter reflections of terrain and marine surfaces obscure the contrast between the target-of-interest and clutter (terrain and marine surface reflections) in synthetic aperture radar (SAR) images and consequently hinder the efficiency of image interpretation and analysis. To overcome this problem, this letter proposes an efficient clutter suppression method in SAR images, which is named shedding irrelevant patterns (SIP). The essence is to construct a regression function that can suppress clutter and preserve the target patterns concurrently. We assume that the clutter is irrelevant to the target-of-interest and distinguishable in patterns in terms of image-pixel distribution and intensity (spatial information). Experimental results show the efficiency of the proposed method in both clutter suppression and target pattern preservation.

Description: A concern in hyperspectral image classification is the high number of required training samples. When traditional classifiers are applied, feature reduction (FR) techniques are the most common approaches to deal with this problem. Subspace-based classifiers, which are developed based on high-dimensional space characteristics, are another way to handle the high dimension of hyperspectral images. In this letter, a novel subspace-based classification approach is proposed and compared with basic and improved subspace-based classifiers. The proposed classifier is also compared with traditional classifiers that are accompanied by an FR technique and the well-known support vector machine classifier. Experimental results prove the efficiency of the proposed method, especially when a limited number of training samples are available. Furthermore, the proposed method has a very high level of automation and simplicity, as it has no parameters to be set.

Description: To improve the spatial density of measurement points of persistent-scatterer interferometry, distributed scatterer (DS) should be considered and processed. An important procedure in DS interferometry is the phase triangulation (PT). This letter introduces two modified PT algorithms (i.e., equal-weighted PT and coherence-weighted PT) and analyzes the mathematical relations between different published PT methods (i.e., the maximum-likelihood phase estimator, least squares estimator, and eigendecomposition-based phase estimators). The analysis shows that the above five PT methods share very similar mathematical forms with different weight values in the estimation procedure.

Description: Sparse representation-based classifier and its variants have been widely adopted for hyperspectral image (HSI) classification recently. However, sparse representation is unstable so that similar features might obtain significantly different sparse codes. Despite the instability, we find that the sparse codes follow a class-dependent distribution under the structured dictionary consisting of training samples from all classes. Based on this observation, a novel discriminative feature, sparse code histogram (SCH), is developed for HSI classification. By counting the SCH of each sample from the sparse codes of its spatial neighbors, we can statistically obtain the distribution pattern of sparse codes of the class to which the sample belongs, and then treat the SCH as a new feature for classification. To reduce the possible outliers among the neighbors, a shape-adaptive neighborhood extractor is also employed to enhance the stability of the histogram feature. Experimental results demonstrate that SCH enjoys a strong discriminative power, which can achieve notably better performance than several state-of-the-art methods for HSI classification with limited training samples.

Description: Several detection statistics have been proposed for detecting fine ground disturbances between two synthetic aperture radar (SAR) images, such as vehicle tracks. The standard method involves estimating a local correlation coefficient between images. Other methods have been proposed using various statistical hypothesis tests. One of these alternative methods is a generalized likelihood ratio test (GLRT), which compares a full-correlation image model to a no-correlation image model. In this letter, we expand the GLRT to polarimetric SAR data and derive the appropriate GLRT detection statistics. Additionally, we explore relaxing the equal variance/equal polarimetric covariance assumptions used in previous results and find improved performance on macroscopic scene changes.

Description: For downward-looking linear array 3-D synthetic aperture radar, the resolution in cross-track direction is much lower than the ones in range and azimuth. Hence, superresolution reconstruction algorithms are desired. Since the cross-track signal to be reconstructed is sparse in the object domain, compressive sensing algorithm has been used. However, the imaging processing on the 3-D scene brings large computational loads, which renders challenges in both data acquisition and processing. To cover this shortage, truncated singular value decomposition is utilized to reconstruct a reduced-redundancy spatial measurement matrix. The proposed algorithm provides advantages in terms of computational time while maintaining the quality of the scene reconstructions. Moreover, our results on uniform linear array are generally applicable to sparse nonuniform linear array. Superresolution properties and reconstruction accuracies are demonstrated using simulations under the noise and clutter scenarios.

Description: This letter proposes a signal processing method of passive bistatic radar (PBR) exploiting an uncooperative radar as an illuminator. Compared with other opportunity illuminators, the transmitting signal of a radar usually has a better ambiguity function, which leads to a higher range resolution. Two channels are needed in PBR system. The reference channel is used to estimate radar signal parameters and reconstruct directly propagated signal. The surveillance channel is used to receive scattered wave. An array antenna and a simultaneous multibeam algorithm are necessary in the surveillance channel due to the flexible beam scanning of the uncooperative radar. The procedure of the proposed method is explained in detail, which is then followed by a field experiment. Preliminary results from the field experiment show that the proposed method can be applied to target angle and bistatic range measurement successfully.

Description: In order to achieve 3-D imaging with an airborne down-looking linear-array synthetic aperture radar (LASAR), a uniform virtual antenna array may be obtained by aperture synthesis of the cross-track sparse multiple-input–multiple-output array. However, the actual 3-D imaging quality is unavoidably degraded by errors in the virtual element position. In this letter, we investigate the effects of these errors on the forms and the degrees of image quality degradation by decomposing the error-related stochastic processes via an orthogonal transform based on discrete Legendre polynomials. It should be noted that these analyses are helpful for designing a LASAR system and providing a reference for specifying the requisite precision of measurement devices and calibration methods. Finally, we briefly consider the use of calibration methods to eliminate the effects of errors.

Description: In problems where labeled data are scarce, semisupervised learning (SSL) techniques are an attractive framework that can exploit both labeled and unlabeled data. These approaches typically rely on a smoothness assumption such that examples that are similar in input space should also be similar in label space. In many domains, such as remotely sensed hyperspectral image (HSI) classification, the data violate this assumption. In response, we propose a general method by which a neighborhood graph used in SSL is learned using binary classifiers that are trained to predict whether a pair of pixels shares the same label. Working within the framework of semisupervised neural networks (SSNNs), we show that our approach improves on the performance of the SSNN on two HSI data sets.

Description: In this letter, new models for the spatial correlation of sea clutter texture and intensity are proposed as improved versions of current power law models or exponential decay model. The models for texture have three unknown parameters, and thus can be called triparametric models. The structure of the models is a weighted sum of two components, which can describe the decaying process of the correlation coefficient with spatial lags, as well as the periodic behavior due to the existence of transient coherent structures in sea clutter. Unknown parameters are optimized by the nonlinear least square fit method. Models for sea clutter intensity can be obtained through a linear transform for uncorrelated speckle based on the compound-Gaussian representation of sea clutter. The proposed models are validated and compared with current models using S- and C-band measured sea clutter data. Analysis results indicate the effectiveness of the proposed models in that they can describe the behavior of spatial correlation coefficients with higher accuracy.

Description: Terrestrial laser scanning (TLS) has become a popular tool for acquiring source data points which can be used to construct digital elevation models (DEMs) for a wide number of applications. A TLS point cloud often has a very fine spatial resolution, which can represent well the spatial variation of a terrain surface. However, the uncertainty in DEMs created from this relatively new type of source data is not well understood, which forms the focus of this letter. TLS survey data representing four terrain surfaces of different characteristics were used to explore the effects of surface complexity and typical TLS data density (in terms of data point spacing) on DEM accuracy. The spatial variation in TLS data can be decomposed into parts corresponding to the signal of spatial variation (of terrain surfaces) and noise due to measurement error. We found a linear relation between the DEM error and the typical TLS data spacings considered (30–100 mm) which arises as a function of the interpolation error, and a constant contribution from the propagated data noise. This letter quantifies these components for each of the four surfaces considered and shows that, for the interpolation method considered here, higher density sampling would not be beneficial.

Description: A novel way to estimate the live fuel moisture content (LFMC) was explored from the ratio of canopy water content (CWC) and foliage dry biomass (FDB). The CWC was estimated using the PROSAIL (PROSPECT + SAIL) radiative transfer model from the Landsat 8 product. A weak constraint 4-D variational data assimilation method was employed to assimilate the temporally estimated leaf area index into a soil-water-atmosphere-plant (SWAP) model for optimizing the model control variables. Then, the SWAP model was reinitialized with this optimum set of control variables, and better prediction of FDB was obtained. Results showed that a high accuracy level was achieved for the estimated CWC ( $R^{2}=0.91$ , $mbox{RMSE}=84.74 mbox{g/m}^2$ ) and FDB ( $R^2=0.88$ , $mbox{RMSE}=48.54 mbox{g/m} ^2$ ) when compared with in situ measured values. However, the accuracy level of estimated LFMC was poor ( $R^2=0.59$ , $mbox{RMSE} =30.85%$ ) . Further analyses find that the estimated LFMC is reliable for low LFMC but challenged for high LFMC, which indicates that the presented method still makes sense to the assessment of wildfire risk since the wildfire generally occurs when the vegetation is in low LFMC condition.

Description: In this letter, we present the use of experimental human micro-Doppler signature data gathered by a multistatic radar system to discriminate between unarmed and potentially armed personnel walking along different trajectories. Different ways of extracting suitable features from the spectrograms of the micro-Doppler signatures are discussed, particularly empirical features such as Doppler bandwidth, periodicity, and others, and features extracted from singular value decomposition (SVD) vectors. High classification accuracy of armed versus unarmed personnel (between 90% and 97% depending on the walking trajectory of the people) can be achieved with a single SVD-based feature, in comparison with using four empirical features. The impact on classification performance of different aspect angles and the benefit of combining multistatic information is also evaluated in this letter.

Description: Automatic urban area detection in remote sensing images is an important application in the field of earth observation. Most of the existing methods employ feature classifiers and thereby contain a data training process. Moreover, some methods cannot detect urban areas in complex scenes accurately. This letter proposes an automatic urban area detection method that uses multiple features that have different resolutions. First, a downsampled low-resolution image is used to segment the candidate area. After the corner points of the urban area are extracted, a weighted Gaussian voting matrix technique is employed to integrate the corner points into the candidate area. Then, the edge features and homogeneous region are extracted by using the original high-resolution image. Using these results as the input, the processes of guided filtering and contrast enhancement can finally detect accurately the urban areas. This method combines multiple features, such as corner, edge, and regional characteristics, to detect the urban areas. The experimental results show that the proposed method has better detection accuracy for urban areas than the existing algorithms.

Description: In marine sciences, time series are often nonlinear and nonstationary. Adequate and specific methods are needed to analyze such series. In this letter, an application of the empirical mode decomposition method (EMD) associated to the Hilbert spectral analysis (HSA) is presented. Furthermore, EMD-based time-dependent intrinsic correlation (TDIC) analysis is applied to consider the correlation between two nonstationary time series. Four temperature time series obtained from automatic measurements in nearshore waters of the Réunion island are considered, recorded every 10 min from July 2011 to January 2012. The application of the EMD on these series and the estimation of their power spectra using the HSA are illustrated. The authors identify low-frequency tidal waves and display the pattern of correlations at different scales and different locations. By TDIC analysis, it was concluded that the high-frequency modes have small correlation, whereas the trends are perfectly correlated.

Description: Detecting vehicles in aerial images provides important information for traffic management and urban planning. Detecting the cars in the images is challenging due to the relatively small size of the target objects and the complex background in man-made areas. It is particularly challenging if the goal is near-real-time detection, i.e., within few seconds, on large images without any additional information, e.g., road database and accurate target size. We present a method that can detect the vehicles on a 21-MPixel original frame image without accurate scale information within seconds on a laptop single threaded. In addition to the bounding box of the vehicles, we extract also orientation and type (car/truck) information. First, we apply a fast binary detector using integral channel features in a soft-cascade structure. In the next step, we apply a multiclass classifier on the output of the binary detector, which gives the orientation and type of the vehicles. We evaluate our method on a challenging data set of original aerial images over Munich and a data set captured from an unmanned aerial vehicle (UAV).

Description: We compare five slope correction methods developed by Walter et al. , Montes et al. , Schleppi et al. , España et al. , and Gonsamo et al. (referred to as WAL, MON, SCH, ESP, and GON, respectively) using artificial fisheye pictures simulated by graphics software and a lookup table (LUT) retrieval method. The LUT is built by simulating the directional gap fraction as a function of leaf area index (LAI) and average leaf inclination angle (ALIA) using the Poisson law. LAI and ALIA estimates correspond to the case of the LUT that provides the lowest root-mean-square error between the observed gap fractions after slope correction and the simulated ones. Three LAI values (1.5, 3.5, and 5.5), four ALIA values (26.8°, 45°, 57.5°, and 63.2°), and three slope angles (0°, 20°, and 50°) constituted 36 samples of random scenes. ESP is recommended because its results are accurate and independent on the leaf angle distribution (LAD), while GON only performs well for spherical LAD. The three other methods present less good performances with underestimation or overestimation of LAI and/or ALIA depending on the LAD, and the recommended order for them is MON, SCH, and WAL.

Description: In this letter, an improved phase correlation (PC) method based on 2-D plane fitting and the maximum kernel density estimator (MKDE) is proposed, which combines the idea of Stone's method and robust estimator MKDE. The proposed PC method first utilizes a vector filter to minimize the noise errors of the phase angle matrix and then unwraps the filtered phase angle matrix by the use of the minimum cost network flow unwrapping algorithm. Afterward, the unwrapped phase angle matrix is robustly fitted via MKDE, and the slope coefficients of the 2-D plane indicate the subpixel shifts between images. The experiments revealed that the improved method can effectively avoid the impact of outliers on the phase angle matrix during the plane fitting and is robust to aliasing and noise. The matching accuracy can reach 1/50th of a pixel using simulated data. The real image sequence tracking experiment was also undertaken to demonstrate the effectiveness of the proposed PC method with a registration accuracy of root-mean-square error better than 0.1 pixels.

Description: Accurately mapping forest carbon density by combining sample plots and remotely sensed images has become popular because this method provides spatially explicit estimates. However, mixed pixels often impede the improvement of the estimation. In this letter, regression modeling and spectral unmixing analysis were integrated to improve the estimation of forest carbon density for the You County of Hunan, China, using Landsat Thematic Mapper images. Linear spectral unmixing with and without a constraint (LSUWC and LSUWOC) and nonlinear spectral unmixing (NSU) were compared to derive the fractions of five endmembers, particularly forests. Stepwise regression, logistic regression, and polynomial regression (PR) with and without the forest fraction used as an independent variable and the product of the forest fraction image and the map from the best model without the forest fraction were compared. The models were developed using 56 sample plots, and their results were validated using 26 test plots. The decomposition of mixed pixels was assessed using higher spatial resolution SPOT images and a corresponding land cover map. The results showed that 1) LSUWC more accurately estimated the endmember fractions than LSUWOC and NSU, 2) PR had the greatest estimation accuracy of forest carbon, and 3) combining regression modeling and spectral unmixing increased the estimation accuracy by 31%–39%, and introducing the forest fraction into the regressions performed better than the product of forest fraction image and the results from PR without the fraction. This implied that the integrations provided great potential in reducing the impacts of mixed pixels in mapping forest carbon.

Description: This letter proposes a multiresolution technique to address the high computational cost in remote sensing image registration. The scale-invariant feature transform is applied to detect keypoints and descriptors, and then, global information combined with descriptors is utilized to establish keypoint mappings. Keypoints are first classified according to their octaves. Then, in the lowest resolution, the keypoints of the largest octave are mapped with descriptors and the global information, giving an initial affine transformation $T_0$ . In the next octave, the keypoints of the second largest octave are mapped by employing $T_0$ to narrow the space of matching keypoints. By this means, the process of establishing keypoint correspondences is conducted from one resolution (octave) to the next as the obtained transformation gets finer until we get to the highest resolution. Due to the high computational expense of computing global information, the proposed technique is important for aligning large-size remote sensing imagery. Experimental results show that the proposed method can achieve a comparable registration accuracy but with a less computational cost than directly building keypoint mappings on images of large size.

Description: Presents corrections to the paper, "Estimation of forest biomass from two-level model inversion of single-pass InSAR data" (Soja, M.J., et al.,Trans. Geosci. Remote Sens., vol. 53, no. 9, pp. 5083???5099, Sep. 2015).

Description: In this letter, we present an efficient parallel implementation of composite kernels in support vector machines (SVMs) for hyperspectral image (HSI) classification. Our implementation makes effective use of commodity graphics processing units (GPUs). Specifically, we port the calculation of composite kernels to GPUs, perform intensive computations based on NVidia's compute unified device architecture, and execute the rest of the operations related with control and small data calculations in the CPU. Our experimental results, conducted using real hyperspectral data sets and NVidia GPU platforms, indicate significant improvements in terms of computational effectiveness, achieving near-real-time performance of spatial–spectral HSI classification for the first time in the literature.

Description: Compared with airborne laser scanning, terrestrial laser scanning (TLS) offers ground-based point cloud data of trees and provides greater potential to accurately estimate tree and stand parameters. However, there is a lack of effective methods to accurately identify locations of individual trees from TLS point cloud data. It is also unknown whether the estimation accuracy of the parameters, including tree height (H), diameter at breast height (DBH), and so on, using TLS can meet the requirement of forest management and planning. In this letter, a novel method to effectively process point cloud data and further determine the locations of individual trees in a stand based on the central coordinates of point cloud data on a defined grid according to the largest DBH was developed. Moreover, a point-cloud-data-based convex hull algorithm and the cylinder method were, respectively, used to estimate DBH and H of individual trees. This study was conducted in a pure Chinese fir plantation of 45 trees located in Huang-Feng-Qiao forest farm, You County of Hunan, China. The comparison of the estimated and observed values showed that the obtained tree locations had errors of less than 20 cm, and the relative root mean square errors for the estimates of both DBH and H were less than 5%. This implies that TLS is very promising for the retrieval of tree and stand parameters in forest stands. For the applications of these methods to mixed forests with a structure of multilayer canopies, further examination is needed.

Description: In this letter, a novel algorithm for attitude measurement based on a 3-D electromagnetic model (3-D em-model) is proposed. The 3-D em-model is established offline based on the geometric structure of the target, and it can be used to predict the scattering features at different target attitudes. In order to measure the attitude of the air target, we design a bistatic step frequency radar system. The directions of the two radars' lines of sight (LOSs) relative to the target are acquired by matching the high-resolution range profiles (HRRPs) from the target echoes to the HRRPs generated from the 3-D em-model. Since the directions of two radars' LOSs relative to the Earth are already known, the absolute attitude of the target can be acquired. The innovative contributions of this letter are as follows: 1) A comprehensive theoretical analysis of air target attitude measurement based on its own 3-D em-model is proposed; 2) the method can be applied to different kinds of air targets such as aircraft, satellite, missile, etc.; 3) the proposed attitude measurement method does not require target motion model in advance; and 4) the proposed algorithm can be applied to any kind of step frequency waveforms. Experiments using both data predicted by a high-frequency electromagnetic code and data measured in the chamber verify the validity of the method.

Description: The Soil Moisture and Ocean Salinity (SMOS) and Soil Moisture Active Passive (SMAP) missions provide Level-1 brightness temperature (Tb) observations that are used for global soil moisture estimation. However, the nature of these Tb data differs: the SMOS Tb observations contain atmospheric and select reflected extraterrestrial (“Sky”) radiation, whereas the SMAP Tb data are corrected for these contributions, using auxiliary near-surface information. Furthermore, the SMOS Tb observations are multiangular, whereas the SMAP Tb is measured at 40° incidence angle only. This letter discusses how SMOS Tb, SMAP Tb, and radiative transfer modeling components can be aligned in order to enable a seamless exchange of SMOS and SMAP Tb data in soil moisture retrieval and assimilation systems. The aggregated contribution of the atmospheric and reflected Sky radiation is, on average, about 1 K for horizontally polarized Tb and 0.5 K for vertically polarized Tb at 40° incidence angle, but local and short-term values regularly exceed 5 K.

Description: Urban areas are one of the most important components of human society. Their extents have been continuously growing during the last few decades. Accurate and timely measurements of the extents of urban areas can help in analyzing population densities and urban sprawls and in studying environmental issues related to urbanization. Urban extents detected from remotely sensed data are usually a by-product of land use classification results, and their interpretation requires a full understanding of land cover types. In this study, for the first time, we mapped urban extents in the continental United States using a novel one-class classification method, i.e., positive and unlabeled learning (PUL), with multi-temporal Moderate Resolution Imaging Spectroradiometer (MODIS) data for the year 2010. The Defense Meteorological Satellite Program Operational Linescan System (DMSP-OLS) night stable light data were used to calibrate the urban extents obtained from the one-class classification scheme. Our results demonstrated the effectiveness of the use of the PUL algorithm in mapping large-scale urban areas from coarse remote-sensing images, for the first time. The total accuracy of mapped urban areas was 92.9% and the kappa coefficient was 0.85. The use of DMSP-OLS night stable light data can significantly reduce false detection rates from bare land and cropland far from cities. Compared with traditional supervised classification methods, the one-class classification scheme can greatly reduce the effort involved in collecting training datasets, without losing predictive accuracy.

Description: The knowledge about impacts of changes in precipitation regimes on terrestrial ecosystems is fundamental to improve our understanding of global environment change, particularly in the context that heavy precipitation is expected to increase according to the 5th Intergovernmental Panel on Climate Change (IPCC) assessment. Based on observed climate data and the Advanced Very High Resolution Radiometer (AVHRR) Global Inventory Modeling and Mapping Studies (GIMMS) satellite-derived normalized difference vegetation index (NDVI), here we analyzed the spatio-temporal changes in grassland NDVI, covering 1.64 × 106 km2, in northern China and their linkages to changes in precipitation and temperature during the period 1982–2011. We found that mean growing season (April–October) grass NDVI is more sensitive to heavy precipitation than to moderate or light precipitation in both relatively arid areas (RAA) and relatively humid areas (RHA), whereas the sensitivities of grass NDVI to temperature are comparable to total precipitation in RHA. Heavy precipitation showed the strongest impacts in more than half of northern China (56%), whereas impacts of light precipitation on grass NDVI were stronger in some areas (21%), mainly distributed in northwestern China, a typical arid and semi-arid area. Our findings suggest that responses of grasslands are divergent with respect to changes in precipitation intensities.

Description: Multiple stable states are established in coastal tidal wetlands (marshes, mangroves, deltas, seagrasses) by ecological, hydrological, and geomorphological feedbacks. Catastrophic shifts between states can be induced by gradual environmental change or by disturbance events. These feedbacks and outcomes are key to the sustainability and resilience of vegetated coastlines, especially as modulated by human activity, sea level rise, and climate change. Whereas multiple stable state theory has been invoked to model salt marsh responses to sediment supply and sea level change, there has been comparatively little empirical verification of the theory for salt marshes or other coastal wetlands. Especially lacking is long-term evidence documenting if or how stable states are established and maintained at ecosystem scales. Laboratory and field-plot studies are informative, but of necessarily limited spatial and temporal scope. For the purposes of long-term, coastal-scale monitoring, remote sensing is the best viable option. This review summarizes the above topics and highlights the emerging promise and challenges of using remote sensing-based analyses to validate coastal wetland dynamic state theories. This significant opportunity is further framed by a proposed list of scientific advances needed to more thoroughly develop the field.

Description: Aquatic vegetation serves many important ecological and socioeconomic functions in lake ecosystems. The presence of floating algae poses difficulties for accurately estimating the distribution of aquatic vegetation in eutrophic lakes. We present an approach to map the distribution of aquatic vegetation in Lake Taihu (a large, shallow eutrophic lake in China) and reduce the influence of floating algae on aquatic vegetation mapping. Our approach involved a frequency analysis over a 2003–2013 time series of the floating algal index (FAI) based on moderate-resolution imaging spectroradiometer (MODIS) data. Three phenological periods were defined based on the vegetation presence frequency (VPF) and the growth of algae and aquatic vegetation: December and January composed the period of wintering aquatic vegetation; February and March composed the period of prolonged coexistence of algal blooms and wintering aquatic vegetation; and June to October was the peak period of the coexistence of algal blooms and aquatic vegetation. By comparing and analyzing the satellite-derived aquatic vegetation distribution and 244 in situ measurements made in 2013, we established a FAI threshold of −0.025 and VPF thresholds of 0.55, 0.45 and 0.85 for the three phenological periods. We validated the accuracy of our approach by comparing the results between the satellite-derived maps and the in situ results obtained from 2008–2012. The overall classification accuracy was 87%, 81%, 77%, 88% and 73% in the five years from 2008–2012, respectively. We then applied the approach to the MODIS images from 2003–2013 and obtained the total area of the aquatic vegetation, which varied from 265.94 km2 in 2007 to 503.38 km2 in 2008, with an average area of 359.62 ± 69.20 km2 over the 11 years. Our findings suggest that (1) the proposed approach can be used to map the distribution of aquatic vegetation in eutrophic algae-rich waters and (2) dramatic changes occurred in the distribution of aquatic vegetation in Lake Taihu during the 11-year study.

Description: The temporal variations (diurnal and annual) in arboreal $(varepsilon_mathrm{Tree})$ and bare soil $(varepsilon_mathrm{Soil})$ dielectric constants and their correlation with precipitation were examined for several trees in Japan. A significant (1 $sigma$ (standard deviation) and 2 $sigma$ ) $varepsilon_mathrm{Tree}$ increase is observed after rainfall at 89.8% and 90.5% probability. However, rainfall does not always induce significant $varepsilon_mathrm{Tree}$ increases. Rainfall of more than 5 mm/day can induce 1 $sigma$ $varepsilon_mathrm{Tree}$ Tree increase at a 59.6% probability. In order to examine whether the increase in $varepsilon_mathrm{Tree}$ affects the L-band $sigma^{0}$ variation in a forest, the four-year temporal variation of the L-band backscattering coefficient $(sigma^{0})$ was estimated from observations by the Advanced Land Observing Satellite Phased Array type L-band Synthetic Aperture Radar. Observed maximum absolute deviations from the mean over the forest area were 1.0 and 1.2 dB for $sigma_{mathrm{HH}}^{0}$ and $sigma_mathrm{HV}^{0}$ , respectively, and 4.0 and 3.0 dB over open land. $sigma^{0}$ and rainfall correlations show that $varepsilon_mathrm{Tree}$ and $sigma_mathrm{Forest}^{0}$ are proportional to precipitation integrated over seven or eight days; $varepsilon_mathrm{Soil}$ and $sigma_mathrm{Open land}^{0}$ are proportional to precipitation integrated over three days. This finding indicates that $varepsilon_mathrm{Tree}$ variations influence $sigma_{mathrm{Forest areas}}^{0}$ . A stronger correlation between $sigma_mathrm{HV}^{0}$ and precipitation is observed in several sites with low $sigma_mathrm{HV}^{0}$ , where less biomass is expected, and several sites with high $sigma_mathrm{HV}^{0}$ , where more biomass is expected. A weaker correlation between $sigma_mathrm{HV}^{0}$ and precipitation is observed for several sites with high $sigma_mathrm{HV}^{0}$ . These differences may be explained by the different contributions of double bounce scattering and potential transpiration, which is a measure of the ability of the atmosphere to remove water from the surface th

Description: The spectral and radiometric quality of airborne imaging spectrometer data is affected by the anisotropic reflectance behavior of the imaged surface. Illumination and observation angle-dependent patterns of surface reflected radiation propagate into products, hinder quantitative assessment of biophysical/biochemical parameters, and decrease the comparability of data from multiple flight lines. The Ross–Li model, originally developed for multiangular observations, can be inverted to estimate and correct for surface anisotropy effects. This requires land cover be stratified into distinct types of scattering behavior. When the observations subsumed in these classes cover a range of view angles, a pseudo multiangular view on the surface can be employed to invert the Ross–Li model. A discrete land cover classification, however, bears the risk of inappropriate scattering correction resulting in spatial artifacts in the corrected data, predominantly in transition regions of two land cover types (e.g., soil and sparse vegetation with varying fractions). We invert the Ross–Li model on continuous land cover fraction layers. We decompose land cover in dominating structural types using linear spectral unmixing. Ross–Li kernel weights and formulations are estimated for each type independently; the correction is then applied pixel-wise according to the fractional distribution. The corrected Airborne Prism EXperiment imaging spectrometer data show significant reduction of anisotropic reflectance effects of up to 90% (average 60% to 75%, $p=0.05$ ), measured in the overlapping regions of adjacent flight lines. No spatial artifacts or spectral irregularities are observed after correction.

Description: Most spaceborne sensors have a tradeoff between high spatial and high temporal resolutions. This tradeoff limits the use of remote sensing data in various applications that require images in both the high spatial and high temporal resolutions. In this paper, we propose a novel technique to create a fine spatial and high temporal resolution images at a ground-based data processing system. Resourcesat-2 is one of the Indian Space Research Organization missions, and it carries the Linear Imaging Self-Scanning sensors (LISS III and LISS IV) and an Advanced Wide-Field Sensor (AWiFS) onboard. The spatial resolution of LISS III is 23.5 m, and that of AWiFS is 56 m. The temporal resolution of LISS III is 24 days, and that of AWiFS is five days. The proposed method creates a synthetic LISS-III image at 23.5-m spatial and five-day temporal resolutions. It is based on the subpixel relationship between a single AWiFS–LISS-III image pair, which is acquired before or after the prediction date. In temporal data composition, spurious spatial discontinuities are inevitable for land-cover type changes. These discontinuities were identified with temporal edge primitives and were smoothed with a spatial-profile-averaging method. A synthetic LISS-III image for time $t_{k}$ is predicted from an AWiFS image at time $t_{k}$ and a single AWiFS–LISS-III image pair at time $t_{0}$ , where $t_{0}ne t_{k}$ . Experimental results demonstrated that the proposed method is superior in terms of the computational efficiency and prediction accuracy with the other existing methods.

Description: Glacier avalanches are natural hazards that could damage infrastructures and threaten lives in high-altitude mountainous terrains. On April 7, 2012, a massive ice avalanche struck a Pakistani base at Gayari sector, Saltoro Valley, and buried/killed 148 soldiers and civilians. Keeping in view the catastrophe, a study was designed with the objectives to: 1) model and simulate the Gayari sector glacier avalanche incident in terms of height, extent, velocity, pressure, and momentum; 2) generate hazard risk assessment of possible other glacier avalanches in the Saltoro Valley through modeling and simulation; and 3) suitability analysis of current camp sites and recommendation of new safe camps sites locations in the Saltoro valley. To simulate the Gayari sector glacier event and other Glacier possible avalanches, a physical process-based rapid mass movements (RAMMS) was used. The RAMMS has two main components, i.e., Voellmy–Salm (VS) model and random kinetic energy, which deals with variables such as avalanche height and the mean avalanche velocity during the course of simulation. The suitability analysis of current camp sites were achieved using weighted overlay analysis with different constraints in ArcGIS Spatial Analyst. The RAMMS model simulation of the Gayari avalanche event predicted a maximum velocity of 74 ms -1 , generating a pressure of 5074 kPa and attaining a height of 45 m, whereas the predicted debris volume on the ground was 3.8145 million m 3 . A good agreement was found between actual debris height and extent, as compared with the RAMMS model output. The RAMMS model simulated all the potential tributary glacier avalanches of Saltoro valley very well. The weighted overlay analyses in ArcGIS Spatial Analyst revealed that the existing camp sites are safe and were not threatened by the glacier avalanche hazard. However, it was recommended that the Gayari camp should not be constructed at the same location and should be relocated- to the proposed safe camp sites identified in this research study. The proposed methodology developed in the current study could be applied in the Siachen conflict zone for avalanche hazard/risk analysis of all the camp sites located in the valley.

Description: Object-oriented remote sensing image classification is becoming more and more popular because it can integrate spatial information from neighboring regions of different shapes and sizes into the classification procedure to improve the mapping accuracy. However, object identification itself is difficult and challenging. Superpixels, which are groups of spatially connected similar pixels, have the scale between the pixel level and the object level and can be generated from oversegmentation. In this paper, we establish a new classification framework using a superpixel-based graphical model. Superpixels instead of pixels are applied as the basic unit to the graphical model to capture the contextual information and the spatial dependence between the superpixels. The advantage of this treatment is that it makes the classification less sensitive to noise and segmentation scale. The contribution of this paper is the application of a graphical model to remote sensing image semantic segmentation. It is threefold. 1) Gradient fusion is applied to multispectral images before the watershed segmentation algorithm is used for superpixel generation. 2) A probabilistic fusion method is designed to derive node potential in the superpixel-based graphical model to address the problem of insufficient training samples at the superpixel level. 3) A boundary penalty between the superpixels is introduced in the edge potential evaluation. Experiments on three real data sets were conducted. The results show that the proposed method performs better than the related state-of-the-art methods tested.

Description: We introduce a practical and accurate model, referred to as “GO4,” to describe near-nadir microwave scattering from the sea surface, and at the same time, we address the issue of the filtered mean square slope (mss) conventionally used in the geometrical optics model. GO4 is a simple correction of this last model, taking into account the diffraction correction induced by the rough surface through what we call an effective mean square curvature (msc). We evaluate the effective msc as a function of the surface wavenumber spectrum and the radar frequency and show that GO4 reaches the same accuracy as the physical optics model in a wide range of incidence and frequency bands with the sole knowledge of the mss and msc parameters. The key point is that the mss entering in GO4 is not the filtered but the total slope. We provide estimation of the effective msc on the basis of classical sea spectrum models. We also evaluate the effective msc from near-nadir satellite data in various bands and show that it is consistent with model predictions. Non-Gaussian effects are discussed and shown to be incorporated in the effective msc. We give some applications of the method, namely, the estimation of the total sea surface mss and the recalibration of relative radar cross sections.

Description: This paper evaluates the radiometric accuracy of observations from the Advanced Technology Microwave Sounder (ATMS) onboard Suomi National Polar-orbiting Partnership and Sondeur Atmospherique du Profil d' Humidité Intropicale par Radiométrie (SAPHIR) onboard Megha-Tropiques through intercalibration and validation versus in situ radiosonde and Global Positioning System Radio Occultation (GPS-RO) observations. SAPHIR and ATMS water vapor channels operate at slightly different frequencies. We calculated the bias due to radiometric errors as the difference between the observed and simulated differences between the two instruments. This difference, which is often referred to as double difference, ranges between 0.3 and 0.7 K, which shows good consistency between the instruments. We used a radiative transfer model to simulate the satellite brightness temperatures (Tbs) using radiosonde and GPS-RO profiles and then compared simulated and observed Tbs. The difference between radiosonde and ATMS Tbs for the middle and upper tropospheric temperature sounding channels was less than 0.5 K at most stations, but the difference between radiosonde and ATMS/SAPHIR Tbs for water vapor channels was between 0.5 and 2.0 K. The larger bias for the water vapor channels is mainly due to several errors in radiosonde humidity observations. The mean differences between the ATMS observations and the Tbs simulated using GPS-RO profiles were 0.2, 0.3, 0.4, 0.2, and −0.2 K for channels 10–14, respectively; and the uncertainty increases from 0.02 K for channel 10 to 0.07 K for channel 14.

Description: This paper presents a completely automatic processing chain for orthorectification of optical pushbroom sensors. The procedure is robust and works without manual intervention from raw satellite image to orthoimage. It is modularly divided in four main steps: metadata extraction, automatic ground control point (GCP) extraction, geometric modeling, and orthorectification. The GCP extraction step uses georeferenced vector roads as a reference and produces a file with a list of points and their accuracy estimation. The physical geometric model is based on collinearity equations and works with sensor-corrected (level 1) optical satellite images. It models the sensor position and attitude with second-order piecewise polynomials depending on the acquisition time. The exterior orientation parameters are estimated in a least squares adjustment, employing random sample consensus and robust estimation algorithms for the removal of erroneous points and fine-tuning of the results. The images are finally orthorectified using a digital elevation model and positioned in a national coordinate system. The usability of the method is presented by testing three RapidEye images of regions with different terrain configurations. Several tests were carried out to verify the efficiency of the procedure and to make it more robust. Using the geometric model, subpixel accuracy on independent check points was achieved, and positional accuracy of orthoimages was around one pixel. The proposed procedure is general and can be easily adapted to various sensors.

Description: This paper presents a new ground moving target imaging (GMTIm) algorithm for airborne synthetic aperture radar (SAR) based on a novel time-frequency representation (TFR), Lv's distribution (LVD). We first analyze generic moving target signatures for a multichannel SAR and then derive the analytical spectrum of a point target moving at a constant velocity by a polar format algorithm for SAR image formation. SAR motion deviation from a predetermined flight track is considered to facilitate airborne SAR applications. LVD, as a recently developed TFR for the analysis of multicomponent linear-frequency-modulated signal, is adopted to represent the target kinematic spectrum in the Doppler centroid frequency and chirp rate domain. As a result, the proposed SAR-GMTIm algorithm is capable of imaging multiple moving targets even when they are located at the same range resolution cell. Some practical issues such as imaging maneuvering targets and small/weak targets are discussed to enhance the applicability of the proposed algorithm. Simulation results with isotropic point moving targets are presented to validate the effectiveness and superiority of the proposed algorithm. Raw data collected by an airborne multichannel SAR are also used to verify the performance improvement made by the proposed algorithm.

Description: Satellite soil moisture estimates have received increasing attention over the past decade. This paper examines the applicability of estimating soil moisture states and soil hydraulic parameters through two particle filter (PF) methods: The PF with commonly used sampling importance resampling (PF-SIR) and the PF with recently developed Markov chain Monte Carlo sampling (PF-MCMC) methods. In a synthetic experiment, the potential of assimilating remotely sensed near-surface soil moisture measurements into a 1-D mechanistic soil water model (HYDRUS-1D) using both the PF-SIR and PF-MCMC algorithms is analyzed. The effects of satellite temporal resolution and accuracy, soil type, and ensemble size on the assimilation of soil moisture are analyzed. In a real data experiment, we first validate the Advanced Microwave Scanning Radiometer--Earth Observing System (AMSR-E) soil moisture products in the Oklahoma Little Washita Watershed. Aside from rescaling the remotely sensed soil moisture, a bias correction algorithm is implemented to correct the deep soil moisture estimate. Both the ascending and descending AMSR-E soil moisture data are assimilated into the HYDRUS-1D model. The synthetic assimilation results indicated that, whereas both updating schemes showed the ability to correct the soil moisture state and estimate hydraulic parameters, the PF-MCMC scheme is consistently more accurate than PR-SIR. For real data case, the quality of remotely sensed soil moisture impacts the benefits of their assimilation into the model. The PF-MCMC scheme brought marginal gains than the open-loop simulation in RMSE at both surface and root-zone soil layer, whereas the PF-SIR scheme degraded the open-loop simulation.

Description: We propose a new deterministic approach for remote sensing retrieval, called modified total least squares (MTLS), built upon the total least squares (TLS) technique. MTLS implicitly determines the optimal regularization strength to be applied to the normal equation first-order Newtonian retrieval using all of the noise terms embedded in the residual vector. The TLS technique does not include any constraint to prevent noise enhancement in the state space parameters from the existing noise in measurement space for an inversion with an ill-conditioned Jacobian. To stabilize the noise propagation into parameter space, we introduce an additional empirically derived regularization proportional to the logarithm of the condition number of the Jacobian and inversely proportional to the L2-norm of the residual vector. The derivation, operational advantages and use of the MTLS method are demonstrated by retrieving sea surface temperature from GOES-13 satellite measurements. An analytic equation is derived for the total retrieval error, and is shown to agree well with the observed error. This can also serve as a quality indicator for pixel-level retrievals. We also introduce additional tests from the MTLS solutions to identify contaminated pixels due to residual clouds, error in the water vapor profile and aerosols. Comparison of the performances of our new and other methods, namely, optimal estimation and regression-based retrieval, is performed to understand the relative prospects and problems associated with these methods. This was done using operational match-ups for 42 months of data, and demonstrates a relatively superior temporally consistent performance of the MTLS technique.

Description: Radar-based remote sensing for measurement of ocean surface waves presents advantages over conventional point sensors such as wave buoys. As its use becomes more widespread, it is important to understand the sensitivity of the extracted wave parameters to the characteristics of the radar and the scatterers. To examine such issues, experiments were performed offshore of the Scripps Institution of Oceanography pier in July 2010. Radar measurements in low wind speeds were performed with a dual-polarized high-resolution X-band pulse-Doppler radar at low grazing angles along with two independent measurements of the surface waves using conventional sensors, a GPS-based buoy, and an ultrasonic array. Comparison between radar cross section (RCS) and Doppler modulations shows peak values occurring nearly in-phase, in contrast with tilt modulation theory. Spectral comparisons between Doppler-based and RCS-based spectra show that Doppler-based spectra demonstrate greater sensitivity to swell-induced modulations, whereas RCS-based spectra show greater sensitivity to small-scale modulations (or generally have more noise at high frequency), and they equally capture energy at the wind wave peak. Doppler estimates of peak period were consistent with the conventional sensors, whereas the RCS differed in assignment of peak period to wind seas rather than swell in a couple of cases. Higher order period statistics of both RCS and Doppler were consistent with the conventional sensors. Radar-based significant wave heights are lower than buoy-based values and contain nontrivial variability of ∼33%. Comparisons between HH and VV polarization data show that VV data more accurately represent the wave field, particularly as the wind speeds decrease.

Description: Polarimetric synthetic aperture radar satellite and ground-based Ku- and X-band scatterometer measurements are used to explore the scattering mechanism for ice in shallow Arctic lakes, wherein strong radiometric responses are seen for floating ice, and low returns are evident where the ice has grounded. Scatterometer measurements confirm that high backscatter is from the ice/water interface, whereas polarimetric decomposition suggests that the dominant scattering mechanism from that interface is single bounce. Using Fresnel equations, a simple model for surface bounce from the ice/water interface is proposed, and its predictions are supported by experimental parameters such as co-pol phase difference, co-pol ratio, and the results of rigorous numerical modeling. Despite early research suggesting double-bounce scattering from columnar air bubbles and the ice/water interface as the dominant scattering mechanism in shallow lakes, this paper strongly supports a single-bounce model.

Description: The on-orbit radiometric calibration of the reflective solar bands (RSBs) of the Visible Infrared Imaging Radiometer Suite (VIIRS) aboard the Suomi National Polar-orbiting Partnership satellite is carried out primarily through observations of a fully illuminated solar diffuser (SD) panel. Accurate knowledge of the solar spectral radiance scattered from the SD is available. The sensor aperture spectral radiance is assumed to be a quadratic polynomial function of a VIIRS detector's background-subtracted response in digital number. The coefficients of the polynomial were initially determined prelaunch. Once on orbit, we assume that these coefficients change uniformly by a common calibration factor, which is referred to as the $F$ -factor. The known solar spectral radiance scattered from a fully illuminated SD allows for the determination of these $F$ -factors. We describe the methodology and the associated algorithms used in the calculation of the RSB $F$ -factors. Our results show that the $F$ -factors change over time, with the largest change occurring at a wavelength of 862 nm (with a value of about 1.55 on day 950 after the satellite launch, relative to its value at the beginning of the launch) . In addition, we estimate the relative error standard deviations of the computed top-of-the-atmosphere reflectance at the detector pixel level. On day 950 of the mission, the relative error standard deviations are all less or equal to 0.016, except for the M11 band (band central wavelength of 2257 nm) , which has a relative error standard deviation of about 0.049 due to a very low signal-to-noise ratio.

Description: Yield and quality estimations provide vital information to fruit growers, yet require accurate monitoring throughout the growing season. To this end, the temporal dependency of fruit yield and quality estimations through spectral vegetation indices was investigated in irrigated and rainfed pear orchards. Both orchards were monitored throughout three consecutive growing seasons, including spectral measurements (i.e., hyperspectral canopy reflectance measurements) as well as yield determination (i.e., total yield and number of fruits per tree) and quality assessment (i.e., fruit firmness, total soluble solids and fruit color). The results illustrated a clear association between spectral vegetation indices and both fruit yield and fruit quality (|r| > 0.75; p 〈 0.001). However, the correlations between vegetation indices and production variables varied throughout the growing season, depending on the phenological stage of fruit development. In the irrigated orchard, index values showed a strong association with production variables near time of harvest (|r| > 0.6; p 〈 0.001), while in the rainfed orchard, index values acquired during vegetative growth periods presented stronger correlations with fruit parameters (|r| > 0.6; p 〈 0.001). The improved planning of remote sensing missions during (rainfed orchards) and after (irrigated orchards) vegetative growth periods could enable growers to more accurately predict production outcomes and improve the production process.

Description: Since the late 1970s, the Chinese government has initiated ecological restoration programs in the Three North Shelter Forest System Project (TNSFSP) area. Whether accelerated climate change will help or hinder these efforts is still poorly understood. Using the updated and extended AVHRR NDVI3g dataset from 1982 to 2011 and corresponding climatic data, we investigated vegetation variations in response to climate change. The results showed that the overall state of vegetation in the study region has improved over the past three decades. Vegetation cover significantly decreased in 23.1% and significantly increased in 21.8% of the study area. An increase in all three main vegetation types (forest, grassland, and cropland) was observed, but the trend was only statistically significant in cropland. In addition, bare and sparsely vegetated areas, mainly located in the western part of the study area, have significantly expanded since the early 2000s. A moisture condition analysis indicated that the study area experienced significant climate variations, with warm-wet conditions in the western region and warm-dry conditions in the eastern region. Correlation analysis showed that variations in the Normalized Difference Vegetation Index (NDVI) were positively correlated with precipitation and negatively correlated with temperature. Ultimately, climate change influenced vegetation growth by controlling the availability of soil moisture. Further investigation suggested that the positive impacts of precipitation on NDVI have weakened in the study region, whereas the negative impacts from temperature have been enhanced in the eastern study area. However, over recent years, the negative temperature impacts have been converted to positive impacts in the western region. Considering the variations in the relationship between NDVI and climatic variables, the warm–dry climate in the eastern region is likely harmful to vegetation growth, whereas the warm–wet conditions in the western region may promote vegetation growth.

Description: Multispectral, as well as multi-temporal, satellite images, coupled with measurements, in situ, have been widely applied to the water quality monitoring of reservoirs. However, the spatial resolutions of the current multispectral satellite imageries are inadequate for trophic state mapping of small reservoirs which merely cover several hectares. Moreover, the temporal gap between effective satellite imaging and measurements, in situ, is usually a few days or weeks; this time lag hampers the establishment of regression models between band ratios and water quality parameters. In this research, the RGB and NIR sensors carried on an unmanned aerial vehicle (UAV) were applied to the trophic state mapping of Tain-Pu reservoir, which is one of the small reservoirs in Kinmen, Taiwan. Due to the limited sampling points and the uncertainty of water fluidity, the average method and the matching pixel-by-pixel (MPP) method were employed to search for the optimal regression models. The experimental results indicate that the MPP method can lead to better regression models than the average method, and the trophic state maps show that the averages of Chl-a, TP, and SD are 179.7 μg·L−1, 108.4 μg·L−1, and 1.4 m, respectively.

Description: In this study, a methodology for glacier elevation reconstruction from Digital Elevation Model (DEM) time series (tDEM) is described for modeling the evolution of glacier elevation and estimating related volume change, with focus on medium-resolution and noisy satellite DEMs. The method is robust with respect to outliers in individual DEM products. Fox Glacier and Franz Josef Glacier in New Zealand are used as test cases based on 31 Advanced Spaceborne Thermal Emission and Reflection Radiometer (ASTER) DEMs and the Shuttle Radar Topography Mission (SRTM) DEM. We obtained a mean surface elevation lowering rate of −0.51 ± 0.02 m·a−1 and −0.09 ± 0.02 m·a−1 between 2000 and 2014 for Fox and Franz Josef Glacier, respectively. The specific volume difference between 2000 and 2014 was estimated as −0.77 ± 0.13 m·a−1 and −0.33 ± 0.06 m·a−1 by our tDEM method. The comparably moderate thinning rates are mainly due to volume gains after 2013 that compensate larger thinning rates earlier in the series. Terminus thickening prevailed between 2002 and 2007.

Description: The aim of this paper is to propose a methodology combing multi-temporal X-band SAR images (TerraSAR-X) with continuous ground thetaprobe measurements, for the retrieval of surface soil moisture and texture at a high spatial resolution. Our analysis is based on seven radar images acquired at a 36° incidence angle in the HH polarization, over a semi-arid site in Tunisia (North Africa). The soil moisture estimations are based on an empirical change detection approach using TerraSAR-X data and ground auxiliary thetaprobe network measurements. Two assumptions were tested: (1) roughness variations during the three-month radar acquisition campaigns were not accounted for; (2) a simple correction for temporal variations in roughness was included. The results reveal a small improvement in the estimation of soil moisture when a correction for temporal variations in roughness is introduced. By considering the estimated temporal dynamics of soil moisture, a methodology is proposed for the retrieval of clay and sand content (expressed as percentages) in soil. Two empirical relationships were established between the mean moisture values retrieved from the seven acquired radar images and the two soil texture components over 36 test fields. Validation of the proposed approach was carried out over a second set of 34 fields, showing that highly accurate clay estimations can be achieved. Maps of soil moisture, clay and sand percentages at the studied site are derived.

Description: In close-range photogrammetry, images are difficult to acquire and organize primarily because of the limited field of view (FOV) of digital cameras when long focal lenses are used to measure large targets. To overcome this problem, we apply a scanning photography method that acquires images by rotating the camera in both horizontal and vertical directions at one station. This approach not only enlarges the FOV of each station but also ensures that all stations are distributed in order without coverage gap. We also conduct a modified triangulation according to the traits of the data overlapping among images from the same station to avoid matching all images with one another. This algorithm synthesizes the images acquired from the same station into synthetic images, which are then used to generate a free network. Consequently, we solve the exterior orientation elements of each original camera image in the free network and perform image matching among original images to obtain tie points. Finally, all original images are combined in self-calibration bundle adjustment with control points. The feasibility and precision of the proposed method are validated by testing it on two fields using 300 and 600 mm lenses. The results confirm that even with a small amount of control points, the developed scanning photogrammetry can steadily achieve millimeter scale accuracy at distances ranging from 40 m to 250 m.

Description: Rice farming in Northeast China is crucially important for China’s food security and sustainable development. A key challenge is how to optimize nitrogen (N) management to ensure high yield production while improving N use efficiency and protecting the environment. Handheld chlorophyll meter (CM) and active crop canopy sensors have been used to improve rice N management in this region. However, these technologies are still time consuming for large-scale applications. Satellite remote sensing provides a promising technology for large-scale crop growth monitoring and precision management. The objective of this study was to evaluate the potential of using FORMOSAT-2 satellite images to diagnose rice N status for guiding topdressing N application at the stem elongation stage in Northeast China. Five farmers’ fields (three in 2011 and two in 2012) were selected from the Qixing Farm in Heilongjiang Province of Northeast China. FORMOSAT-2 satellite images were collected in late June. Simultaneously, 92 field samples were collected and six agronomic variables, including aboveground biomass, leaf area index (LAI), plant N concentration (PNC), plant N uptake (PNU), CM readings and N nutrition index (NNI) defined as the ratio of actual PNC and critical PNC, were determined. Based on the FORMOSAT-2 imagery, a total of 50 vegetation indices (VIs) were computed and correlated with the field-based agronomic variables. Results indicated that 45% of NNI variability could be explained using Ratio Vegetation Index 3 (RVI3) directly across years. A more practical and promising approach was proposed by using satellite remote sensing to estimate aboveground biomass and PNU at the panicle initiation stage and then using these two variables to estimate NNI indirectly (R2 = 0.52 across years). Further, the difference between the estimated PNU and the critical PNU can be used to guide the topdressing N application rate adjustments.

Description: Light Detection and Ranging (Lidar) is a state of the art technology to assess forest aboveground biomass (AGB). To date, methods developed to relate Lidar metrics with forest parameters were built upon the vertical component of the data. In multi-layered tropical forests, signal penetration might be restricted, limiting the efficiency of these methods. A potential way for improving AGB models in such forests would be to combine traditional approaches by descriptors of the horizontal canopy structure. We assessed the capability and complementarity of three recently proposed methods for assessing AGB at the plot level using point distributional approach (DM), canopy volume profile approach (CVP), 2D canopy grain approach (FOTO), and further evaluated the potential of a topographical complexity index (TCI) to explain part of the variability of AGB with slope. This research has been conducted in a mountainous wet evergreen tropical forest of Western Ghats in India. AGB biomass models were developed using a best subset regression approach, and model performance was assessed through cross-validation. Results demonstrated that the variability in AGB could be efficiently captured when variables describing both the vertical (DM or CVP) and horizontal (FOTO) structure were combined. Integrating FOTO metrics with those of either DM or CVP decreased the root mean squared error of the models by 4.42% and 6.01%, respectively. These results are of high interest for AGB mapping in the tropics and could significantly contribute to the REDD+ program. Model quality could be further enhanced by improving the robustness of field-based biomass models and influence of topography on area-based Lidar descriptors of the forest structure.

Description: This paper presents a rigorous validation of five widely used global land cover products, i.e., GLCC (Global Land Cover Characterization), UMd (University of Maryland land cover product), GLC2000 (Global Land Cover 2000 project data), MODIS LC (Moderate Resolution Imaging Spectro-radiometer Land Cover product) and GlobCover (GLOBCOVER land cover product), and a national land cover map GLCD-2005 (Geodata Land Cover Dataset for year 2005) against an independent reference data set over China. The land cover reference data sets in three epochs (1990, 2000, and 2005) were collected on a web-based prototype system using a sampling-based labeling approach. Results show that, in China, the highest overall accuracy is observed in GLCD-2005 (72.3%), followed by MODIS LC (68.9%), GLC2000 (65.2%), GlobCover (57.7%) and GLCC (57.2%), while UMd has the lowest accuracy (48.6%); all of the products performed best in representing “Trees” and “Others”, well with “Grassland” and “Cropland”, but problematic with “Water” and “Urban” across China in general. Moreover, in respect of GLCD-2005, there are significant accuracy differences across seven geographical locations of China, ranging from 46.3% in the Southwest, 77.5% in the South, 79.2% in the Northwest, 80.8% in the North, 81.8% in the Northeast, 82.6% in the Central, to 89.0% in the East. This study indicates that a regionally focused land cover map would in fact be more accurate than extracting the same region from a globally produced map.

Description: This study proposes a method to recognize façades from large-scale urban scenes based on multi-level image features utilizing a recently developed oblique aerial photogrammetry technique. The method involves the use of multi-level image features, a bottom-up feature extraction procedure to produce regions of interest through monoscopic analysis, and then a coarse-to-fine feature matching strategy to characterise and match the regions in a stereoscopic model. Feature extraction from typical urban Manhattan scenes is based on line segments. Windows are re-organised based on the spatial constraints of line segments and the homogeneous structure of the spectrum. Façades as regions of interest are successfully constructed with a remarkable single edge and evidence from windows to get rid of occlusion. Feature matching is hierarchically performed beginning from distinctive facades and regularly distributed windows to the sub-pixel point primitives. The proposed strategy can effectively solve ambiguity and multi-solution problems in the complex urban scene matching process, particularly repetitive and poor-texture façades in oblique view.

Description: In this study, an algorithm to retrieve the gravimetric vegetation water content (GVWC, %) of corn was developed. First, the method for obtaining the optical depth from L-band (1.4 GHz) bi-angular, dual-polarized brightness temperatures (TB) for short vegetation was investigated. Then, the quantitative relationship between the corn optical depth, corn GVWC and corn leaf area index (LAI) was constructed. Finally, using the Polarimetric L-band Microwave Radiometer (PLMR) airborne data in the 2012 Heihe Watershed Allied Telemetry Experimental Research (HiWATER) project, the Global Land Surface Satellite (GLASS) LAI product, the height and areal density of the corn stalks, the corn GVWC was estimated (corn GLASS-GVWC). Both the in situ measured corn GVWC and the corn GVWC retrieved based on the in situ measured corn LAI (corn LAINET-GVWC) were used to validate the accuracy of the corn GLASS-GVWC. The results show that the GVWC retrieval method proposed in this study is feasible for monitoring the corn GVWC. However, the accuracy of the retrieval results is highly sensitive to the accuracy of the LAI input parameters.

Description: As the fifth largest freshwater lake in China, Lake Chaohu has drawn increasing attention due to the decline in water quality and the occurrence of massive algal blooms. We applied an algae pixel-growing algorithm to MODIS Terra or Aqua data (2100 images) to characterize surface floating algae bloom dynamics from 2000 to 2013 with respect to meteorological and lake nutrient conditions. The results show an increase in surface algal bloom coverage, frequency, and duration with a trend toward earlier bloom formation. Importantly, spatial and temporal patterns in the historically less compromised eastern and middle lake areas show that water quality conditions are deteriorating. This has occurred at the same time as lake management has made a catchment scale effort to reduce impact. Our results show that nutrient concentrations were not the main driver of inter-annual bloom variations. Local meteorological conditions, in particular wind speed and temperature, played an important role in the dynamics of floating algal bloom. This highlights the important challenges for lake management.

Description: Burn severity metrics and classification have yet to be tested for many eastern U.S. deciduous vegetation types, but, if suitable, would be valuable for documenting and monitoring landscape-scale restoration projects that employ prescribed fire treatments. Here we present a performance analysis of the Composite Burn Index (CBI) and its relationship to spectral data (differenced Normalized Burn Ratio (dNBR) and its relative form (RdNBR)) across an oak woodland - grassland landscape in southwestern Oklahoma, USA. Correlation and regression analyses were used to compare CBI strata, assess models describing burn severity, and determine thresholds for burn severity classes. Confusion matrices were used to assess burn severity classification accuracy. Our findings suggest that dNBR and RdNBR, thresholded using total CBI, can produce an accurate burn severity map in oak woodlands, particularly from an initial assessment period. Lower accuracies occurred for burn severity classifications of grasslands and raises questions related to definitions and detection of burn severity for grasslands, particularly in transition to more densely treed structures such as savannas and woodlands.

Description: We present simplified expressions for the cross-polarized backscatter of a randomly rough surface predicted by the second-order small-slope approximation (SSA2). The simplification is based on appropriate polynomial approximations of the SSA2 kernel function. We obtain numerically efficient expressions for the cross-polarized backscattering amplitude of a deterministic surface in the form of a single space integral involving only the surface elevation and the second (mixed) derivative of the surface elevation. The ensemble average normalized radar cross section is then derived under a Gaussian random process assumption for the surface. The resulting expression has the form of a Kirchhoff integral involving the roughness correlation function and its second- and fourth-order cross-derivatives. Further simplification is achieved for off-nadir observations using a high-frequency approximation; the result is an analytical formula involving only the resonant curvature and the radar-filtered mean square slope in the out-of-plane direction. A numerical validation of the simplified expressions is provided by comparison with exact SSA2 predictions in representative test cases. The dependence of cross-polarized backscattering on the incidence angle as well as wind speed and direction is then investigated for the case of a directional sea surface model. At near nadir incidence, a clear maximum in azimuth of the cross-polarized backscatter is observed for radar look directions 45 $^{circ}$ from the wind direction.

Description: The wide field of view (WFV) is an optical imaging sensor on-board the Gao Fen 1 (GF-1). The WFV lacks an on-board calibrator, so on-orbit radiometric calibration is required. Zhong et al. proposed a method for cross-calibrating the charge-coupled device on-board the Chinese Huan Jing 1 (HJ-1/CCD) that can be applied to the GF-1/WFV. However, the accuracy is limited because of the wider radiometric dynamic range and the higher spatial resolution of the GF-1/WFV. Therefore, Landsat-8 Operational Land Imager (OLI) imagery with a radiometric resolution similar to that of the GF-1/WFV and DEM extracted from ZY-3 three-line array panchromatic camera (TLC) with a higher spatial resolution were used. A calibration site with uniform surface material and a natural topographic variation was selected, and a model of this site’s bidirectional reflectance distribution function (BRDF) was developed. The model has excellent agreement with the real situation, as shown by the comparison of the simulations to the actual OLI surface reflectance. Then, the model was used to calibrate the WFV. Compared with the TOA reflectance from synchronized Landsat-8/OLI images, all errors calculated with the calibration coefficients retrieved in this paper are less than 5%, much less than the errors calculated with the calibration coefficients given by the China Centre for Resource Satellite Data and Application (CRESDA).

Description: Taking photographs with a commercially available digital camera is an efficient and objective method for determining the green fractional vegetation cover (FVC) for field validation of satellite products. However, classifying leaves under shadows in processing digital images remains challenging and results in classification errors. To address this problem, an automatic shadow-resistant algorithm in the Commission Internationale d’Eclairage L*a*b* color space (SHAR-LABFVC) based on a documented FVC estimation algorithm (LABFVC) is proposed in this paper. The hue saturation intensity (HSI) is introduced in SHAR-LABFVC to enhance the brightness of shaded parts of the image. The lognormal distribution is used to fit the frequency of vegetation greenness and to classify vegetation and the background. Real and synthesized images are used for evaluation, and the results are in good agreement with the visual interpretation, particularly when the FVC is high and the shadows are deep, indicating that SHAR-LABFVC is shadow resistant. Without specific improvements to reduce the shadow effect, the underestimation of FVC can be up to 0.2 in the flourishing period of vegetation at a scale of 10 m. Therefore, the proposed algorithm is expected to improve the validation accuracy of remote sensing products.

Description: Due to their image-based nature, “contextual” approaches are very attractive to estimate evapotranspiration (ET) from remotely-sensed land surface temperature (LST) data. Their application is however limited to highly heterogeneous areas where the soil and vegetation temperature endmembers (Tends) can be observed at the thermal sensor resolution. This paper aims to develop a simple theoretical approach to estimate Tends independently from LST images. Soil Tends are simulated by a soil energy balance model forced by meteorological data. Vegetation Tends are obtained from soil Tends and air temperature. Model-derived soil Tends are first evaluated with in situ measurements made over an irrigated area in Morocco. The root mean square difference (RMSD) between modeled and ground-based soil Tends is estimated as 2.4 ºC. Model-derived soil Tends are next compared with the soil Tends retrieved from 90-m resolution ASTER (Advanced Spaceborne Thermal Emission and Reflection Radiometer) data collected over two irrigated areas in Mexico and Spain. Such a comparison reveals a strong consistency between model-derived and high-resolution image-based soil Tends. A recent contextual ET model (SEB-1S) is then applied to 90-m resolution and to 1-km resolution (aggregated) ASTER data using the model-derived or image-based Tends as the input. The RMSD between 90-m resolution SEB-1S and in situ ET is estimated as 65 and 82 W-m2, and the RMSD between 1-km resolution SEB-1S and aggregated SEB-1S ET is estimated as 78 and 56 W-m2 for the image-based and model-derived Tends, respectively. In light of the above results, Tends should be estimated a priori when contextual models are applied to low resolution images. Moreover, the consistency over highly heterogeneous areas between model-derived and high-resolution image-based Ten

Description: The Earth’s surface net radiation controls the energy and water exchanges between the Earth’s surface and the atmosphere, and can be derived from satellite observations. The ability to monitor the net surface radiation over large areas at high spatial and temporal resolution is essential for many applications, such as weather forecasting, short-term climate prediction or water resources management. The objective of this paper is to derive the net surface radiation in the shortwave domain at high temporal (half-hourly) and spatial resolution (~1 km) using visible imagery from Geostationary Operational Environmental Satellite (GOES). The retrieval algorithm represents an adaptation to GOES data of a standard algorithm initially developed for the NASA-operated Clouds and Earth’s Radiant Energy System (CERES) scanner. The methodology relies on: (1) the estimation of top of atmosphere shortwave radiation from GOES spectral measurements; and (2) the calculation of net surface shortwave (SW) radiation accounting for atmospheric effects. Comparison of GOES-retrieved net surface shortwave radiation with ground-measurements at the National Oceanic and Atmospheric Administration’s (NOAA) Surface Radiation (SURFRAD) stations yields very good agreement with average bias lower than 5 W·m−2 and root mean square difference around 70 W·m−2. The algorithm performance is usually higher over areas characterized by low spatial variability in term of land cover type and surface biophysical properties. The technique does not involve retrieval and assessment of cloud properties and can be easily adapted to other meteorological satellites around the globe.

Description: The Soil Moisture and Ocean Salinity (SMOS) mission was initiated in 2009 with the goal of acquiring global soil moisture data over land using multi-angular L-band radiometric measurements. Specifically, surface soil moisture was estimated using the L-band Microwave Emission of the Biosphere (L-MEB) radiative transfer model. This study evaluated the applicability of this model to the Heihe River Basin in Northern China for specific underlying surfaces by simulating brightness temperature (BT) with the L-MEB model. To analyze the influence of a ground sampling strategy on the simulations, two resampling methods based on ground observations were compared. In the first method, the simulated BT of each point observation was initially acquired. The simulations were then resampled at a 1 km resolution. The other method was based on gridded data with a resolution of 1 km averaged from point observations, such as soil moisture, soil temperature, and soil texture. The simulated BTs at a 1 km resolution were then obtained using the L-MEB model. Because of the large variability in soil moisture, the resampling method based on gridded data was used in the simulation. The simulated BTs based on the calibrated parameters were validated using airborne L-band data from the Polarimetric L-band Multibeam Radiometer (PLMR) acquired during the HiWATER project. The root mean square errors (RMSEs) between the simulated results and the PLMR data were 6 to 7 K for V-polarization and 3 to 5 K for H-polarization at different angles. These results demonstrate that the model effectively represents agricultural land surfaces, and this study will serve as a reference for applying the L-MEB model in arid regions and for selecting a ground sampling strategy.

Description: How to produce the difference data of the two temporal images is a crucial factor in image change detection. In this letter, we propose multicontextual mutual information data (MMID) based on the bivariate Gaussian distribution (BGD) for synthetic aperture radar (SAR) image change detection and illustrate their superiorities over the classical difference data. MMID, which are an improved form of image spatial mutual information, are constructed based on the quadrilateral Markov random field (QMRF) and can be factored into the linear combination of the entropies. Then to adapt MMID to the change detection, we construct the 2-D entropies based on the BGD. In this way, MMID are able to capture the intertemporal statistical dependence of the two temporal images and thus can be taken as the feature-level difference data rather than the pixel-level data. The maximum-likelihood method, the automatic threshold method, and the Markov random field method are performed on the MMID of the real two temporal SAR images for the change detection. Experimental results demonstrate the superiorities of MMID over the traditional difference data.

Description: Spectral unmixing has been a popular technique for analyzing remotely sensed hyperspectral images. The goal of unmixing is to find a collection of pure spectral constituents (called endmembers ) that can explain each (possibly mixed) pixel of the scene as a combination of endmembers, weighted by their coverage fractions in the pixel or abundances . Over the last years, many algorithms have been presented to address the three main parts of the spectral unmixing chain: 1) estimation of the number of endmembers; 2) identification of the endmember signatures; and 3) estimation of the per-pixel fractional abundances. However, to date, there is no standardized tool that integrates these algorithms in a unified framework. In this letter, we present HyperMix, an open-source tool for spectral unmixing that integrates different approaches for spectral unmixing and allows building unmixing chains in graphical fashion, so that the end-user can define one or several spectral unmixing chains in fully configurable mode. HyperMix provides efficient implementations of most of the algorithms used for spectral unmixing, so that the tool automatically recognizes if the computer has a graphics processing unit (GPU) available and optimizes the execution of these algorithms in the GPU. This allows for the execution of spectral unmixing chains on large hyperspectral scenes in computationally efficient fashion. The tool is available online from http://hypercomphypermix.blogspot.com.es and has been validated with real hyperspectral scenes, providing state-of-the-art unmixing results.

Description: Seismic signals are nonlinear, and the seismic state-space model can be described as a nonlinear system. The particle filter (PF) method, as an effective method for estimating the state of a nonlinear system, can be applied to deal with seismic random noise attenuation. However, PF suffers from sample impoverishment caused by resampling, which results in serious loss of valid seismic information and leads to inaccurate representation of the reflected signal. To address the impoverishment issue and to further improve the particle quality, we propose a novel method to suppress seismic random noise—the adaptive fission particle filter (AFPF). In AFPF, all the particles undergo a fission process and produce “offspring” particles to maintain particle diversity. To implement the adaptation and to monitor the degree of fission, we apply a fission factor, which takes into account weights that indicate the quality of the particles. This leads to significant improvements in the particle quality, i.e., the proportion of highly weighted particles is increased. The effective seismic information provided by the resulting particles reproduces the true signal more reliably, reducing the bias of PF. In addition, we establish a dynamic state-space model suitable for seismic signals. Experimental results on synthetic records and field data illustrate the superior performance of AFPF in noise attenuation and reflected signal preservation compared with the PF.

Description: Due to an oversight by the authors there is a correction necessary in this publication [1]. We would like to change the word “phase” to “step” throughout the manuscript. We apologize for any inconvenience caused to readers.[...]

Description: Spartina alterniflora is one of the most hazardous invasive plant species in China. Monitoring the changes in dominant plant species can help identify the invasion mechanisms of S. alterniflora, thereby providing scientific guidelines on managing or controlling the spreading of this invasive species at Jiuduansha Wetland in Shanghai, China. However, because of the complex terrain and the inaccessibility of tidal wetlands, it is very difficult to conduct field experiments on a large scale in this wetland. Hence, remote sensing plays an important role in monitoring the dynamics of plant species and its distribution on both spatial and temporal scales. In this study, based on multi-spectral and high resolution (〈10 m) remote sensing images and field observational data, we analyzed spectral characteristics of four dominant plant species at different green-up phenophases. Based on the difference in spectral characteristics, a decision tree classification was built for identifying the distribution of these plant species. The results indicated that the overall classification accuracy for plant species was 87.17%, and the Kappa Coefficient was 0.81, implying that our classification method could effectively identify the four plant species. We found that the area of Phragmites australi showed an increasing trend from 1997 to 2004 and from 2004 to 2012, with an annual spreading rate of 33.77% and 31.92%, respectively. The area of Scirpus mariqueter displayed an increasing trend from 1997 to 2004 (12.16% per year) and a decreasing trend from 2004 to 2012 (−7.05% per year). S. alterniflora has the biggest area (3302.20 ha) as compared to other species, accounting for 51% of total vegetated area at the study region in 2012. It showed an increasing trend from 1997 to 2004 and from 2004 to 2012, with an annual spreading rate of 130.63% and 28.11%, respectively. As a result, the native species P. australi was surrounded and the habitats of S. mariqueter were occupied by S. alterniflora. The high proliferation ability and competitive advantage for S. alterniflora inhibited the growth of other plant species and we anticipate a continuous expansion of this invasive species at Jiuduansha Wetland. Effective measures should be taken to control the invasion of S. alterniflora.

Description: The accuracy of different workflows using Structure-from-Motion and Multi-View-Stereo techniques (SfM-MVS) is tested. Twelve point clouds of the Corral del Veleta rock glacier, in Spain, were produced with two different software packages (123D Catch and Agisoft Photoscan), using Low Dynamic Range images and High Dynamic Range compositions (HDR) for three different years (2011, 2012 and 2014). The accuracy of the resulting point clouds was assessed using benchmark models acquired every year with a Terrestrial Laser Scanner. Three parameters were used to estimate the accuracy of each point cloud: the RMSE, the Cloud-to-Cloud distance (C2C) and the Multiscale-Model-to-Model comparison (M3C2). The M3C2 mean error ranged from 0.084 m (standard deviation of 0.403 m) to 1.451 m (standard deviation of 1.625 m). Agisoft Photoscan overcome 123D Catch, producing more accurate and denser point clouds in 11 out 12 cases, being this work, the first available comparison between both software packages in the literature. No significant improvement was observed using HDR pre-processing. To our knowledge, this is the first time that the geometrical accuracy of 3D models obtained using LDR and HDR compositions are compared. These findings may be of interest for researchers who wish to estimate geomorphic changes using SfM-MVS approaches.

Description: This study investigated the spatiotemporal variation of vegetation growth and the influence of climatic drivers from 1982 to 2011 across China using datasets from the Normalized Difference Vegetation Index (NDVI) and climatic drivers. Long term trends, significance and abrupt change points of interannual NDVI time series were analyzed. We applied both simple regression and multi-regression models to quantify the effects of climatic drivers on vegetation growth and compare their relative contributions. Results show that on average, the growing season NDVI significantly increased by 0.0007 year

Description: The upcoming satellite mission EnMAP offers the opportunity to retrieve information on the seasonal development of vegetation parameters on a regional scale based on hyperspectral data. This study aims to investigate whether an analysis method for the retrieval of leaf area index (LAI), developed and validated on the 4 m resolution scale of six airborne datasets covering the 2012 growing period, is transferable to the spaceborne 30 m resolution scale of the future EnMAP mission. The widely used PROSAIL model is applied to generate look-up-table (LUT) libraries, by which the model is inverted to derive LAI information. With the goal of defining the impact of different selection criteria in the inversion process, different techniques for the LUT based inversion are tested, such as several cost functions, type and amount of artificial noise, number of considered solutions and type of averaging method. The optimal inversion procedure (Laplace, median, 4% inverse multiplicative noise, 350 out of 100,000 averages) is identified by validating the results against corresponding in-situ measurements (n = 330) of LAI. Finally, the best performing LUT inversion (R2 = 0.65, RMSE = 0.64) is adapted to simulated EnMAP data, generated from the airborne acquisitions. The comparison of the retrieval results to upscaled maps of LAI, previously validated on the 4 m scale, shows that the optimized retrieval method can successfully be transferred to spaceborne EnMAP data.

Description: The synoptic determination of ocean circulation using the data acquired from space, with a coherent depiction of its turbulent characteristics, remains a fundamental challenge in oceanography. This determination has the potential of revealing all aspects of the ocean dynamic variability on a wide range of spatiotemporal scales and will enhance our understanding of ocean–atmosphere exchanges at superresolution, as required in the present context of climate change. Here, we show a four-year time series of spatial superresolution (4 km) turbulent ocean dynamics generated from satellite data using emerging ideas in signal processing coming from nonlinear physics, low-resolution dynamics, and superresolution oceanic sea surface temperature data acquired from optical sensors. The method at its core consists in propagating across the scales the low-resolution dynamics in a multiresolution analysis computed on adimensional critical transition information.

Description: This paper presents the theory, algorithm, and results of a new bistatic interferometry synthetic aperture radar (InSAR) method. It employs the data acquired in an innovative bistatic configuration, which uses the orbital sensors as transmitters of opportunity and the stationary receivers on the ground, to generate a digital elevation model (DEM). In the bistatic spaceborne/stationary InSAR configuration, the interferometric phase only depends on the target-receiver range, which could not be obtained directly from the measured bistatic range. Therefore, the conventional transforming relationship between the interferometric phase and the topographic height is no longer practical. In order to solve the problem, we introduce a new conversion relationship between the interferometric phase and the topographic height, which is derived by the model of the ellipsoidal projection in the bistatic configuration. Meanwhile, the error analysis of the new conversion is carried out through a simulation. Both the simulated and measured data are used to test and verify the feasibility of the new bistatic InSAR method. In the spaceborne/stationary InSAR experiment, YaoGan-3 (an L-band spaceborne SAR system launched by China) was selected as the transmitter and two stationary receivers were mounted on the top of a tall building. The generated DEM of high quality shows that the presented method performs very well in the bistatic InSAR data process.

Description: We present a nonparametric Bayesian hierarchical model (HDP_IBPs) to classify very high resolution panchromatic satellite images in an unsupervised way, in which the hierarchical Dirichlet process (HDP) and Indian buffet process (IBP) are combined on multiple scenes. The main contribution of this paper is a novel application framework to solve the problems of traditional probabilistic topic models and achieve the effective unsupervised classification of very high resolution (VHR) panchromatic satellite images. In this framework, a VHR satellite image is first oversegmented into basic processing units and divided into a set of subimages. We use the Chinese restaurant franchise process as a construct method of the HDP to capture the latent semantic structures (i.e., classes) and the class proportion (i.e., co-occurrence of topics) for each subimage. Meanwhile, the subimages are grouped into different scenes based on the class proportion. Finally, the IBP is employed to select the most appropriate classes for each subimage from all of the classes based on different scenes in turn. The hierarchical structure of our model transmits the spatial information from the original image to the scene layer implicitly and provides useful cues of classification. The experimental results show that HDP_IBPs outperforms state-of-the-art models in terms of both qualitative and quantitative evaluations.

Description: This paper describes a linear-image-transform-based algorithm for reducing stripe noise, track line artifacts, and motion-induced errors in remote sensing data. Developed for multibeam bathymetry (MB), the method has also been used for removing scalloping in synthetic aperture radar images. The proposed image transform is the composition of an invertible edge detection operator and a fast discrete Radon transform (DRT) due to Götz, Druckmüller, and Brady. The inverse DRT is computed by using an iterative method and exploiting an approximate inverse algorithm due to Press. The edge operator is implemented by circular convolution with a Laplacian point spread function modified to render the operator invertible. In the transformed image, linear discontinuities appear as high-intensity spots, which may be reset to zero. In MB data, a second noise signature is linked to motion-induced errors. A Chebyshev approximation of the original image is subtracted before applying the transform, and added back to the denoised image; this is necessary to avoid boundary effects. It is possible to process data faster and suppress motion-induced noise further by filtering images in nonoverlapping blocks using a matrix representation for the inverse DRT. Processed test images from several MB data sets had less noise and distortion compared with those obtained with standard low-pass filters. Denoising also improved the accuracy in statistical classification of geomorphological type by 10–28% for two sets of invariant terrain features.

Description: Sparse graph embedding (SGE) is a promising technique useful for the nonlinear feature extraction (FE) of hyperspectral images (HSIs). However, such images exhibit spatial variability and spectral multimodality, presenting challenges to existing FE methods, including SGE. To address this issue, this paper presents two novel SGE methods for HSI classification. One method, which is termed simultaneous SGE (SSGE), is designed to consider the spatial variability of spectral signatures by using a simultaneous sparse representation (SSR) model integrated with a shape-adaptive neighborhood building approach. In addition, a sparse graph is constructed via matrix computation based on sparse codes. Then, low-dimensional features are produced by employing linear graph embedding (LGE) based on the constructed sparse graph. The other method, which is termed simultaneous sparse multimanifold learning (SSMML), is proposed to handle the multimodality of an HSI. In SSMML, multiple views are generated to represent different modalities. Then, multiview-oriented submanifolds are produced by adopting SSGE, and they are further integrated via coregularization. SSGE is capable of modeling both local and global data structures. Furthermore, SSMML serves as a prototype that can model multimodal data structures. The proposed methods are evaluated by using sparse multinomial logistic regression for HSI classification. Experimental results with two popular hyperspectral data sets validate the good performance of the two methods in producing more representative low-dimensional features and yielding superior classification results compared with other related approaches.

Description: We present a new algorithm that simultaneously retrieves aerosol properties and land surface bidirectional reflectance distribution function (BRDF) over Australian from Advanced Along-Track Scanning Radiometer images. Three key components are addressed: 1) an analytical radiative transfer formulation, based on Green's function, linking top-of-atmosphere (TOA) reflectance to the surface BRDF; 2) a novel approach to modeling BRDF using an extended compositing period; and 3) a set of representative aerosol models based on a published typology of Australian aerosols. Due to the generally low aerosol loadings and widespread bright surfaces over Australia, BRDF modeling is crucial. By using a 9-month compositing period, 90% of the Australian continent can be modeled with an error in the forward-to-nadir reflectance ratio of 2.5% or less. Comparison with suitably processed Moderate Resolution Imaging Spectroradiometer BRDF/albedo products demonstrates satisfactory agreement. For the studying period from 2002 to 2008, validation of aerosol optical depth (AOD) against eight sun photometers across Australia encompassing widely different atmospheric and surface regimes shows high accuracy, with a mean absolute error in AOD at 550 nm of 0.03 and a bias of 0.007. About 60% of the matchup points are within an absolute error of 0.03, 80% are within 0.05, and 96% are within 0.1. The algorithm selects for each cell an optimal aerosol type from a set of four predefined candidate models. Continental aerosol maps derived from the new method indicate broadscale agreement with known seasonal aerosol sources, while providing new insights into the spatial and temporal distributions of aerosol over the Australian continent.

Description: This investigation examines the phenomenology effects of the squint angle on the morphology of moving target smears in spotlight synthetic aperture radar (SAR). This analysis includes both the smears resulting from standard image formation applied to simulated radar measurements as well as the theoretical predictions for the central contours of the signatures. In particular, this paper generates the down-range and cross-range components of the predicted central 2-D contours of mover signatures, including the locations of the cross-range offsets. The analytics for squinted geometry include additional contributions in the signature contour equations that do not arise for the case of broadside imaging. These terms can affect the overall contour morphology, particularly in terms of shape and extent. Numerous examples are presented to demonstrate that the signature prediction equations yield excellent agreement with standard image formation with simulated radar data. Therefore, this analysis can provide an effective tool in predicting the shape, extent, and location of smears due to arbitrarily moving surface targets for squinted spotlight SAR.

Description: In this paper, we consider signal subspace estimation based on low-rank representation for hyperspectral imagery. It is often assumed that major signal sources occupy a low-rank subspace. Due to the mixed nature of hyperspectral remote sensing data, the underlying data structure may include multiple subspaces instead of a single subspace. Therefore, in this paper, we propose the use of low-rank subspace representation to estimate the number of subspaces in hyperspectral imagery. In particular, we develop simple estimation approaches without user-defined parameters because these parameters can be fixed as constants. Both real data experiments and computer simulations demonstrate excellent performance of the proposed approaches over those currently in the literature.

Description: The conventional spaced-receiver approach uses amplitude scintillations to estimate equatorial ionospheric irregularity drift velocities. This approach is less applicable at high latitudes where there is a lack of substantial amplitude scintillations. This paper presents a method to estimate ionosphere irregularity horizontal drift velocities based on GPS signal carrier phase measurements. Joint time–frequency analysis of the carrier phase measurements using an adaptive periodogram technique generates time-varying spectrograms of ionospheric irregularity-induced phase fluctuations. Cross correlation of the spectrograms between antenna pairs provides time lag information on propagating radio signals through the same ionospheric structure. The time lag information is combined with known positions of the receiver array, satellite orbits, and assumed irregularity altitude to infer ionospheric irregularity horizontal drift velocity. This paper presents the methodology and demonstrates its feasibility using data collected by a GPS receiver array at Gakona, Alaska. The potential error sources of this method are also analyzed.

Description: Forest inventories based on field sample surveys, supported by auxiliary remotely sensed data, have the potential to provide transparent and confident estimates of forest carbon stocks required in climate change mitigation schemes such as the REDD+ mechanism. The field plot size is of importance for the precision of carbon stock estimates, and better information of the relationship between plot size and precision can be useful in designing future inventories. Precision estimates of forest biomass estimates developed from 30 concentric field plots with sizes of 700, 900, …, 1900 m2, sampled in a Tanzanian rainforest, were assessed in a model-based inference framework. Remotely sensed data from airborne laser scanning (ALS) and interferometric synthetic aperture radio detection and ranging (InSAR) were used as auxiliary information. The findings indicate that larger field plots are relatively more efficient for inventories supported by remotely sensed ALS and InSAR data. A simulation showed that a pure field-based inventory would have to comprise 3.5–6.0 times as many observations for plot sizes of 700–1900 m2 to achieve the same precision as an inventory supported by ALS data.

Description: This paper presents a rapid multi-return ALS-based (Airborne Laser Scanning) tree trunk detection approach. The multi-core Divide & Conquer algorithm uses a CBH (Crown Base Height) estimation and 3D-clustering approach to isolate points associated with single trunks. For each trunk, a principal-component-based linear model is fitted, while a deterministic modification of LO-RANSAC is used to identify an optimal model. The algorithm returns a vector-based model for each identified trunk while parameters like the ground position, zenith orientation, azimuth orientation and length of the trunk are provided. The algorithm performed well for a study area of 109 trees (about 2/3 Norway Spruce and 1/3 European Beech), with a point density of 7.6 points per m2, while a detection rate of about 75% and an overall accuracy of 84% were reached. Compared to crown-based tree detection methods, the aTrunk approach has the advantages of a high reliability (5% commission error) and its high tree positioning accuracy (0:59m average difference and 0:78m RMSE). The usage of overlapping segments with parametrizable size allows a seamless detection of the tree trunks.

Description: 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).

Description: 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.

Description: 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.

Description: 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.

Description: 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.

Description: 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.

Description: 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.

Description: 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.

Description: 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.

Description: 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.

Description: 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.

Description: 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.

Description: 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.

Description: Fractional vegetation cover (FVC) plays an important role in earth surface process simulations, climate modeling, and global change studies. Several global FVC products have been generated using medium spatial resolution satellite data. However, the validation results indicate inconsistencies, as well as spatial and temporal discontinuities of the current FVC products. The objective of this paper is to develop a reliable estimation algorithm to operationally produce a high-quality global FVC product from the Moderate Resolution Imaging Spectroradiometer (MODIS) surface reflectance. The high-spatial-resolution FVC data were first generated using Landsat TM/ETM+ data at the global sampling locations, and then, the general regression neural networks (GRNNs) were trained using the high-spatial-resolution FVC data and the reprocessed MODIS surface reflectance data. The direct validation using ground reference data from validation of land European Remote Sensing instruments sites indicated that the performance of the proposed method ( $mathbf{R}^{2}=0.809$ , $mathbf{RMSE} =0.157$ ) was comparable with that of the GEOV1 FVC product ( $mathbf{R}^{2}=0.775$ , $mathbf{RMSE} =0.166$ ) , which is currently considered to be the best global FVC product from SPOT VEGETATION data. Further comparison indicated that the spatial and temporal continuity of the estimates from the proposed method was superior to that of the GEOV1 FVC product.

Description: This paper proposes a region kernel to measure the region-to-region distance similarity for hyperspectral image (HSI) classification. The region kernel is designed to be a linear combination of multiscale box kernels, which can handle the HSI regions with arbitrary shape and size. Integrating labeled pixels and labeled regions, we further propose a region-kernel-based support vector machine (RKSVM) classification framework. In RKSVM, three different composite kernels are constructed to describe the joint spatial–spectral similarity. Particularly, we design a desirable stack composite kernel that consists of the point-based kernel, the region-based kernel, and the cross point-to-region kernel. The effectiveness of the proposed RKSVM is validated on three benchmark hyperspectral data sets. Experimental results show the superiority of our region kernel method over the classical point kernel methods.

Description: Atmospheric refraction in the troposphere causes the propagation speed of electromagnetic signals to be less than the light speed. This creates a difference between the actual propagation path delay and the distance of the geometrical straight-line path, i.e., a quantity known as the tropospheric delay. As classical imaging algorithms for spaceborne synthetic aperture radar (SAR) do not take the tropospheric delay into account, imaging filters are designed based on the assumption of rectilinear propagation with the light speed. Therefore, a residual phase exists in imaging results, which affects focusing quality under the condition of high resolution. In order to compensate for the impact of tropospheric delay on focusing performance, this paper modifies the spaceborne SAR echo model and then proposes an imaging compensation algorithm. The key to this algorithm is to fit a range delay coefficient based on the European Geostationary Navigation Overlay Service model of zenith delay and Niell mapping function, which projects the zenith delay onto the looking direction. After range compensation, classical imaging, and azimuth compensation, which compose the proposed algorithm, the processed results are well focused.