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: 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: In this paper, the classification via sprepresentation and multitask learning is presented for target recognition in SAR image. To capture the characteristics of SAR image, a multidimensional generalization of the analytic signal, namely the monogenic signal, is employed. The original signal can be then orthogonally decomposed into three components: 1) local amplitude; 2) local phase; and 3) local orientation. Since the components represent the different kinds of information, it is beneficial by jointly considering them in a unifying framework. However, these components are infeasible to be directly utilized due to the high dimension and redundancy. To solve the problem, an intuitive idea is to define an augmented feature vector by concatenating the components. This strategy usually produces some information loss. To cover the shortage, this paper considers three components into different learning tasks, in which some common information can be shared. Specifically, the component-specific feature descriptor for each monogenic component is produced first. Inspired by the recent success of multitask learning, the resulting features are then fed into a joint sparse representation model to exploit the intercorrelation among multiple tasks. The inference is reached in terms of the total reconstruction error accumulated from all tasks. The novelty of this paper includes 1) the development of three component-specific feature descriptors; 2) the introduction of multitask learning into sparse representation model; 3) the numerical implementation of proposed method; and 4) extensive comparative experimental studies on MSTAR SAR dataset, including target recognition under standard operating conditions, as well as extended operating conditions, and the capability of outliers rejection.

Description: Offsets of synthetic aperture radar (SAR) images have played an important role in deriving complete three-dimensional (3-D) surface displacement fields in geoscientific applications. However, offset maps often suffer from multiple outliers and patch-like artifacts, because the standard offset-measurement method is a regular moving-window operation that does not consider the scattering characteristics of the ground. Here, we show that by focusing the offset measurements on predetected strong reflectors, the reliability and accuracy of SAR offsets can be significantly improved. Application to the 2011 Van (Turkey) earthquake reveals a clear deformation signal from an otherwise decorrelated interferogram, making derivation of the 3-D coseismic displacement field possible. Our proposed method can improve mapping of coseismic deformation and other ground displacements, such as glacier flow and landslide movement when strong reflectors exist.

Description: We present an accelerated probabilistic learning concept and its prototype implementation for mining heterogeneous Earth observation images, e.g., multispectral images, synthetic aperture radar (SAR) images, image time series, or geographical information systems (GIS) maps. The system prototype combines, at pixel level, the unsupervised clustering results of different features, extracted from heterogeneous satellite images and geographical information resources, with user-defined semantic annotations in order to calculate the posterior probabilities that allow the final probabilistic searches. The system is able to learn different semantic labels based on a newly developed Bayesian networks algorithm and allows different probabilistic retrieval methods of all semantically related images with only a few user interactions. The new algorithm reduces the computational cost, overperforming existing conventional systems, under certain conditions, by several orders of magnitude. The achieved speed-up allows the introduction of new feature models improving the learning capabilities of knowledge-driven image information mining systems and opening them to Big Data environments.

Description: The mineralogy and environmental history of Mars are being extensively studied through remote sensing observations paired with laboratory and in situ experiments. A significant portion of these experiments is being devoted to the identification and mapping of different iron oxides present in the Martian terrains. Among these compounds, goethite has been an object of great interest since its occurrence can be interpreted as mineralogical evidence of past aqueous activity on those landscapes. Although such experiments can provide valuable information regarding the presence of these minerals, the scope of the resulting observations may be hindered by logistics and cost-related constraints. We believe that predictive computer simulations can be employed to mitigate some of these constraints and contribute to the generation and validation of hypotheses in this area. Accordingly, we propose the use of SPLITS ( Sp ectral Li ght T ransport Model for S and) in investigations involving the spectral signatures of iron-bearing regions of Mars. In this paper, we initially demonstrate the predictive capabilities of the SPLITS model in this context through qualitative comparisons of modeled results with actual observations and measured data. Using the resulting modeled reflectance curves as our baseline data, we then perform a series of controlled computational experiments to investigate how variations on goethite and hematite content affect the spectral responses of Martian sand-textured soils.

Description: Thermal infrared (TIR) time series are exploited by many methods based on Earth observation (EO), for such applications as agriculture, forest management, and meteorology. However, due to physical limitations, data acquired by a single sensor are often unsatisfactory in terms of spatial or temporal resolution. This issue can be tackled by using remotely sensed data acquired by multiple sensors with complementary features. When nonreal-time functioning or at least near real-time functioning is admitted, the measurements can be profitably fed to a sequential Bayesian algorithm, which allows to account for the correlation embedded in the successive acquisitions. In this work, we focus on applications that allow the batch processing of the whole data sequences acquired in a fixed time interval. In this case, multiple options for improving the final product are offered by the Bayesian framework, based on both sequential and smoothing techniques. We consider several such Bayesian strategies and comparatively assess their performances in practical applications and through real thermal data acquired by the SEVIRI and MODIS sensors, encompassing the presence of multiple disturbance source, e.g., the cloud coverage of the illuminated scene.

Description: Land surface albedo, qualifying the ratio of the radiant flux reflected from the land surface to the incident flux, is a key forcing parameter controlling the Earth’s energy budget. Previously, several BRDF archetypes were distilled from high-quality MODIS BRDF/Albedo products. In this study, we propose a method that largely relies on matching observed multiangular reflectances with the most appropriate of these prior BRDF archetypes to determine the amplitude and shape of the actual surface BRDFs, when directional signatures are insufficient. This method is first evaluated using an assortment of multisource BRDF data sets to demonstrate its viability for surface albedo estimates, and then is applied to airborne wide-angle infrared dual-mode line/area array scanner (WIDAS) from the Watershed Allied Telemetry Experimental Research (WATER) campaign in the Heihe River Basin of China in 2008. This algorithm makes use of the linear MODIS BRDF model to determine the BRDF archetypes needed as prior knowledge for intrinsic spectral albedo estimates. The intrinsic spectral albedos are then used to estimate actual spectral albedos by considering the proportion of direct and diffuse solar radiation. A spectral-to-broadband conversion is performed to generate the broadband albedo at shortwave regimes through the use of conversion coefficients derived from extensive radiative transfer simulations. A further validation confirms that the estimated albedos are consistent with in situ field measured albedos over available corn crop sites. This method provides a major advantage on utilizing generalized BRDF information derived from MODIS in conjunction with other instrument data that are acquired with less angular variation.

Description: Ocean current is highly related to the interaction between ocean and atmosphere. By measuring the speed and direction of the ocean current from space, we can investigate the ocean–atmosphere interaction on a global scale. The ocean–atmosphere interaction helps to maintain the balance that is essential for planet habitability. However, the conventional scatterometer is unable to measure the ocean current vector. To achieve this, a potentially feasible approach is to use a bigger antenna, a higher PRF, and measure the interferometric phase of two successive echoes. This paper derives four decorrelation factors, and provides the phase error model first. Then, an end-to-end simulation model is established, and it is used to analyze the feasibility of ocean surface current measurement from space. Based on the simulation model, the system parameters are optimized. The simulation results show that the current speed standard deviation (Std), which means the measurement accuracy, in along-track and cross-track direction is smaller than 0.1 m/s when the wind speed is larger than 4 m/s. The swath can be used for current vector inversion that is greater than 70% when the wind speed is larger than 7 m/s. Meanwhile, ${{bf K}_{{bf pc}}}$ of the modified scatterometer is computed and the results show that ${{bf K}_{{bf pc}}}$ is better than the traditional pencil-beam rotating scatterometer when the wind speed is larger than 6 m/s.

Description: High-resolution three-dimensional (3-D) radar imaging of space targets in micromotion plays a significant role in space target recognition and space situation awareness; thus, it has attracted extensive attention in recent years. Because of the fast rotation, some scattering centers are occluded by others, i.e., the scattering centers cannot be continuously illuminated by radar in the imaging interval, and their radar echoes are discontinuous. In this paper, a nonparametric 3-D imaging method based on scattering center trajectory association is proposed. It deals with target occlusion using the Riemannian manifold optimization and obtains focused imaging of targets in complex micromotion. The effectiveness of the proposed method is validated using simulated data.

Description: Recently, satellite-based systems have been introduced that utilize angle-of-arrival (AOA) measurements to geo-locate objects of interest. In the previous work, we considered the application of nonlinear optimization to AoA-based geolocation to these systems. This previous work, however, assumed that all noise sources were independent. In the case of fast-moving objects, however, there is a significant source of error due to the propagation time inherent in satellite-based observation of objects due to the difference between the location of the object when it is observed by a satellite, and the location of the object when it emitted the signal that is being measured. This introduces a systematic error into the system that cannot be resolved by the system proposed by Burchett et al. In this paper, we extend our prior work to account for the time-delay inherent in satellite-based geolocation systems, making this system accurate for fast-movers as well as fixed or slow-moving objects. Results demonstrating significant improvement in geolocation performance both in terms of accuracy and estimated error bounds are presented.

Description: Most of existing change detection methods could be classified into three groups, the traditional pixel-based change detection (PBCD), the object-based change detection (OBCD), and the hybrid change detection (HCD). Nevertheless, both PBCD and OBCD have disadvantages, and classical HCD methods belong to intuitive decision-level fusion schemes of PBCD and OBCD. There is no optimum HCD method as of yet. Analyzing the complementarities of PBCD and OBCD method, we propose a new unsupervised algorithm-level fusion scheme (UAFS-HCD) in this paper to improve the accuracy of PBCD using spatial context information through: 1) getting the preliminary change mask with PBCD at first to estimate some parameters for OBCD; 2) deriving the unchanged area mask to eliminate the areas without changes, reducing error amplification phenomenon of OBCD; and 3) obtaining the final change mask by means of OBCD method. Taking flood detection with multitemporal SAR data as an example, we compared the new scheme with some classical methods, including PBCD, OBCD, and HCD method and supervised manual trial-and-error procedure (MTEP). The experimental results of flood detection showed that the new scheme was efficient and robust, and its accuracy sometimes can even exceed MTEP.

Description: This paper describes a methodology to extract a consistent human settlement extent layer using Landsat data and its implementation in the Google Earth Engine platform. The approach allows the extraction of human settlement extents by means of the existing Landsat 5 and 7 data sets, allowing to check their evolution at 30-m spatial resolution. Since human settlements are the main proxy to people geographical distribution and to building locations, this layer may serve as a mean to disaggregate people/building counts at the regional/national level. The approach is tested in several parts of the world against existing ground truth data at the same spatial resolution in Brazil and China, as well as against extents manually extracted from VHR data in three different geographical areas: 1) Brazil; 2) South East China; and 3) Indonesia.

Description: Remote sensing of urban areas is mainly conducted at high frequencies, to obtain highly resolved images for classification, target detection, or urban areas monitoring for instance. We propose on the contrary to investigate the use of VHF-band for observation, between 120 and 360 MHz. In particular, our concern is to allow target detection, when the object is not in the line of sight (LOS) of the radar, which is typically an issue in urban areas. The benefits of low frequencies are highligted using measurements and simulations over a simple metallic model of two buildings surrounding a street.

Description: The atmospheric condition parameters used in the radiative transfer-based atmospheric correction (AC) are often uncertain. This uncertainty propagates to the estimated reflectance. The reflectance, is, however, not equally sensitive to all the parameters. A sensitivity analysis (SA) helps in prioritizing the parameters. The objective of this study was to perform an SA of reflectance to water vapor concentration ( $wv$ ) and aerosol optical thickness ( $AOT$ ). SA was performed using the Fourier amplitude sensitivity test (FAST) method, which computes sensitivity indices ( $text{SI}$ ) of these parameters. Besides variation in the two parameters, we also studied the effect of surface albedo on the $text{SI}$ by quantifying $text{SI}$ for three target surfaces (in the spectral range $text{0.44{-}0.96};upmu$ m): 1) a dark target (water); 2) a bright target (bare soil); and 3) a target having low albedo in the visible and high albedo in near-infrared range (forest). For $AOT$ , high ( $approx!text{0.9}$ ) $text{SI}$ values were observed at the nonwater absorption wavelengths. For $wv$ , high $text{SI}$ values were observed at wavelengths, where strong absorption features are loca- ed and when the surface albedo was high. For the dark target, the effect of $AOT$ was prominent throughout the spectral range. We found that the sensitivity of reflectance to $wv$ and $AOT$ is a function of wavelength, strength of the absorption features, and surface albedo. We conclude that $AOT$ is a more important parameter for dark targets than $wv$ even at the principal absorption feature. For bright targets, the importance of $wv$ and $AOT$ depends on the strength of the absorption feature.

Description: The feasibility study of the HALESIS (High-Altitude Luminous Events Studied by Infrared Spectro-imagery) project is presented. The purpose of this experiment is to measure the atmospheric perturbation in the minutes following the occurrence of transient luminous events (TLEs) from a stratospheric balloon in the altitude range of 20–40 km. The instrumentation will include a spectro-imager embedded in a pointing gondola. Infrared signatures of a single blue jet were simulated under the assumption of local thermodynamic equilibrium (LTE), and were then compared with a panel of commercially available instrument specifications. The sensitivity of the signatures with a local perturbation of the main vibrational energy level populations of ${mathbf{rm CO}_{2}}$ , CO, NO, ${mathbf{rm O}_{3}}$ , and ${mathbf{rm H}_{2}}{rm O}$ was measured and the infrared signatures of a single blue jet taking into account non-LTE hypotheses were compared with the same panel of commercially available instrument specifications. Lastly, the feasibility of the study is discussed.

Description: To improve accuracy and efficiency of object detection and classification with hyperspectral imagery (HSI), we propose a novel smoothing algorithm by coupling of a Laplacian-based reaction term to a classical divergence-based anisotropic diffusion partial differential equation (PDE). In addition, an adaptive parameter is introduced to regularize this nonlinear reaction-diffusion PDE by explicitly integrating the interband correlations with the noise level of each band in HSI. It is also well-known that the interband correlations can be implicitly embedded into the diffusion coefficient of the divergence-based PDE, to allow a selective smoothing that reduces the local homogeneous area variability while preventing smoothing across edges. Therefore, the interband correlations in HSI are exploited in the proposed method in both explicit and implicit ways. As a result, our algorithm is more effective at controlling the behavior of the diffusion evolution when compared to previous multi/hyperspectral diffusion algorithms. The simulations based on both synthetic data and real hyperspectral remote sensing data show that our algorithm can improve the hyperspectral data quality in terms of both visual inspection and image quality indices.

Description: This paper proposes a novel approach to classify hyperspectral (HS) images using both spectral and spatial information. It first consists of a supervised spectral dimension reduction step that transforms the HS image into a score image that has fewer channels. These channels are chosen so as to enhance distances between classes to be discriminated and to reduce background variability, thus leading to edges that correspond to actual class borders. In the second step, applying an edge-preserving spatial regularization on this score image leads to a lowered background variability. Therefore, in the third step, the pixel-wise classification of the regularized score image is greatly improved. We implement this approach using the partial least squares (PLS) method for spectral dimension reduction and the anisotropic diffusion for spatial regularization. We then compare linear discriminant analysis (LDA), K-nearest neighbors (KNN), and support vector machine (SVM) for the class decision. The effectiveness of our method was evaluated with three remotely sensed HS images. Its robustness was also assessed for different training sets, since the latter has a crucial influence on classification performance. On average, our method gave better results in terms of classification accuracy and was more robust than other classification methods tested with the same images.

Description: Due to its simple, fast, and good generalization ability, extreme learning machine (ELM) has recently drawn increasing attention in the pattern recognition and machine learning fields. To investigate the performance of ELM on the hyperspectral images (HSIs), this paper proposes two spatial–spectral composite kernel (CK) ELM classification methods. In the proposed CK framework, the single spatial or spectral kernel consists of activation–function-based kernel and general Gaussian kernel, respectively. The proposed methods inherit the advantages of ELM and have an analytic solution to directly implement the multiclass classification. Experimental results on three benchmark hyperspectral datasets demonstrate that the proposed ELM with CK methods outperform the general ELM, SVM, and SVM with CK methods.

Description: This paper proposes a framework for hyperspectral images (HSIs) classification with composite kernels discriminant analysis (CKDA). The CKDA uses the spectral and spatial information extracted by Gaussian weighted local mean operator (GWLM) and is suitable to solve few labeled samples classification problem of HSI, which has very important practical significance for the case that training samples are insufficient due to high cost. Experimental results show that the spatial information extracted by GWLM can greatly improve the performance, and demonstrate the superiority of CKDA for HSI classification in the case of few labeled samples. Compared with other state-of-the-art spectral–spatial kernel methods, the proposed methods also show very good advantages, especially the parallel kernel method.

Description: Accurate generation of a land cover map using hyperspectral data is an important application of remote sensing. Multiple classifier system (MCS) is an effective tool for hyperspectral image classification. However, most of the research in MCS addressed the problem of classifier combination, while the potential of selecting classifiers dynamically is least explored for hyperspectral image classification. The goal of this paper is to assess the potential of dynamic classifier selection/dynamic ensemble selection (DCS/DES) for classification of hyperspectral images, which consists in selecting the best (subset of) optimal classifier(s) relative to each input pixel by exploiting the local information content of the image pixel. In order to have an accurate as well as computationally fast DCS/DES, we proposed a new DCS/DES framework based on extreme learning machine (ELM) regression and a new spectral–spatial classification model, which incorporates the spatial contextual information by using the Markov random field (MRF) with the proposed DES method. The proposed classification framework can be considered as a unified model to exploit the full spectral and spatial information. Classification experiments carried out on two different airborne hyperspectral images demonstrate that the proposed method yields a significant increase in the accuracy when compared to the state-of-the-art approaches.

Description: Advanced classifiers, e.g., partial least squares discriminant analysis (PLS-DA) and random forests (RF), have been recently used to model reflectance spectral data in general, and of soil properties in particular, since their spectra are multivariate and highly collinear. Preprocessing transformations (PPTs) can improve the classification accuracy by increasing the variability between classes while decreasing the variability within classes. Such PPTs are common practice prior to a PLS-DA, but are rarely used for RF. The objectives of this paper are twofold: to compare the performances of PLS-DA and RF for modeling the spectral reflectance of soil in changed land-uses with different treatments and to compare the effects of nine different PPTs on the prediction accuracy of each of these classification methods. Differences in six physical, biological, and chemical soil properties of changed land-uses from the northern Negev Desert in Israel were evaluated. Significant differences were found between soil properties, which are used to classify land-uses and treatments. Depending on the dataset, different PPTs improved the classification accuracy by 11%–24% and 32%–42% for PLS-DA and RF, respectively, in comparison to the spectra without PPT. Out of the PPTs tested, the generalized least squares weighting (GLSW)-based transformations were found to be the most effective for most classifications using both PLS-DA and RF. Our results show that both PLS-DA and RF are suitable classifiers for spectral data, provided that an appropriate PPT is applied.

Description: This paper contributes the concept of spectral–spatial kernel-based multivariate analysis (KMVSSA) based on the statistical principle of multivariate statistics. The essence of proposed framework is to expose the inherent structure and meaning revealed within spectral and spatial features through various statistical methods in hyperspectral remotely sensed data. This kernel-based framework is investigated to incorporate the spectral and spatial information simultaneously for dimension reduction and classification of hyperdimensional datasets. The method uses multivariate analysis to choose and apply a transform matrix that the transformed components are as orthogonal as possible. This nonlinear framework is derived by means of the theory of complete orthonormal systems. KMVSSA exhibits great flexibility by the combination of spectral and spatial features. We investigate the possibility of using KMVSSA for the classification of hyperspectral images and dimension reduction. The proposed framework is examined and compared in different merits with several hyperspectral images in different conditions (urban/agricultural area and size of the training set). Experimental results show that the proposed framework can meaningfully enhance the dimensionality reduction and also it greatly improves the overall as well as per class classification accuracies. We demonstrate a comprehensive comparison of some state of the art hyperspectral image classification methods.

Description: Development of intelligent decision-making systems for complex problems, such as land covers classification of hyperspectral remote sensing (HSRS) images, requires efficient interpretation of available information through conceptual rather than numerical level. Granular neural network (GNN) in combination with the granular representation of information using linguistic terms is one such system. GNN takes the fuzzified input information and processes them with neural network (NN) architecture, where the network structure is transparent enough to interpret the processing steps. Further, knowledge encoding has been considered as one of the principal elements of intelligent decision-making systems. This paper proposes a new model of knowledge-encoded GNNs for land cover classification of HSRS images. Knowledge encoding is done using neighborhood rough sets (NRSs) that explore the local/contextual information. The encoded knowledge using NRS is obtained in the form of dependency rules with respect to the output class labels of land covers and these rules determine appropriate number of hidden nodes of GNNs. The dependency factors obtained during rule generation are used for initializing the connecting weights of GNNs. NRS is also used here in the selection of a subset of features for reducing the burden of high-dimensional fuzzy-granulated feature space of HSRS image. The proposed model thus explores jointly the advantages of fuzzy granulation, GNNs, and feature selection and knowledge encoding using NRS. Superiority of the model to similar other methods are justified in land covers classification of two HSRS images acquired by different remotely placed sensors.

Description: Several popular endmember extraction and unmixing algorithms are based on the geometrical interpretation of the linear mixing model, and assume the presence of pure pixels in the data. These endmembers can be identified by maximizing a simplex volume, or finding maximal distances in subsequent subspace projections, while unmixing can be considered a simplex projection problem. Since many of these algorithms can be written in terms of distance geometry, where mutual distances are the properties of interest instead of Euclidean coordinates, one can design an unmixing chain where other distance metrics are used. Many preprocessing steps such as (nonlinear) dimensionality reduction or data whitening, and several nonlinear unmixing models such as the Hapke and bilinear models, can be considered as transformations to a different data space, with a corresponding metric. In this paper, we show how one can use different metrics in geometry-based endmember extraction and unmixing algorithms, and demonstrate the results for some well-known metrics, such as the Mahalanobis distance, the Hapke model for intimate mixing, the polynomial post-nonlinear model, and graph-geodesic distances. This offers a flexible processing chain, where many models and preprocessing steps can be transparently incorporated through the use of the proper distance function.

Description: Progress in mapping plant species remotely with imaging spectroscopy data is limited by the traditional classification framework, which carries the requirement of exhaustively defining all classes (species) encountered in a landscape. As the research objective may be to map only one or a few species of interest, we need to explore alternative classification methods that may be used to more efficiently detect a single species. We compared the performance of three support vector machine (SVM) methods designed for single-class detection—binary (one-against-all) SVM, one-class SVM, and biased SVM—in detecting five focal tree and shrub species using data collected by the Carnegie Airborne Observatory over an African savanna. Prior to this comparison, we investigated the effects of training data amount and balance on binary SVM and evaluated alternative methods for tuning one-class and biased SVMs. A key finding was that biased SVM was generally best parameterized by crown-level cross validation paired with the tuning criterion proposed by Lee and Liu [1] . Among the different single-class methods, binary SVM showed the best overall performance (average F-scores 0.43–0.78 among species), whereas one-class SVM showed very poor performance (F-scores 0.09–0.46). However, biased SVM produced results similar to those obtained with binary SVM (F-scores 0.40–0.72), despite using labeled training data from only the focal class. Our results indicate that both binary and biased SVMs can work well for remote single-species detection, while both methods, particularly biased SVM, greatly reduce the amount of training data required compared with traditional multispecies classification.

Description: Hundreds of narrow contiguous spectral bands collected by a hyperspectral sensor have provided the opportunity to identify the various materials present on the surface. Moreover, spatial information, enforcing the assumption that the adjacent pixels belong to the same class with a high probability, is a valuable complement to the spectral information. In this paper, two predominant approaches have been developed to exploit the spatial information. First, by decomposing each pixel and the spatial neighborhood into a low-rank form, the spatial information can be efficiently integrated into the spectral signatures. Meanwhile, in order to describe the low-rank structure of the decomposed data more precisely, an $ell_{1/2}$ norm regularization is introduced and a discrete algorithm is proposed to solve the combined optimization problem by the augmented Lagrange multiplier (ALM) and a half-threshold operator. Second, a graph cuts segmentation algorithm has been applied on the sparse-representation-based probability estimates of the hyperspectral data to further improve the spatial homogeneity of the material distribution. Experimental results on four real hyperspectral data with different spectral and spatial resolutions have demonstrated the effectiveness and versatility of the proposed spatial information-fused approaches for hyperspectral image classification.

Description: Classification of hyperspectral imagery using few labeled samples is a challenging problem, considering the high dimensionality of hyperspectral imagery. Classifiers trained on limited samples with abundant spectral bands tend to overfit, leading to weak generalization capability. To address this problem, we have developed an enhanced ensemble method called multiclass boosted rotation forest (MBRF), which combines the rotation forest algorithm and a multiclass AdaBoost algorithm. The benefit of this combination can be explained by bias-variance analysis, especially in the situation of inadequate training samples and high dimensionality. Furthermore, MBRF innately produces posterior probabilities inherited from AdaBoost, which are served as the unary potentials of the conditional random field (CRF) model to incorporate spatial context information. Experimental results show that the classification accuracy by MBRF as well as its integration with CRF consistently outperforms the other referenced state-of-the-art classification methods when limited labeled samples are available for training.

Description: Curse of dimensionality is a major disadvantage for classification of hyperspectral imagery since a large number of bands need to be dealt with. Band selection is a task to reduce the number of bands. An unsupervised band selection method is proposed in this article. It is a three-step procedure. In the first step, characteristics (attributes) of the bands are found out. Next, redundancy among the bands is removed by executing clustering operation. At last, the remaining bands, which are nonredundant among themselves, are ranked according to their discriminating capability. Discriminating capability is calculated by measuring the capacitory discrimination of the bands. Results are compared with four state-of-the-art methods: a band elimination method, a ranking-based, and two clustering-based band selection methods to demonstrate the effectiveness of the proposed method. Four evaluation measures, namely: 1) classification accuracy; 2)  Kappa coefficient; 3) class separability, and 4) entropy, are calculated over the selected bands to assess the efficiency of the selected bands. The proposed method shows promising results compared to them.

Description: A fast forward feature selection algorithm is presented in this paper. It is based on a Gaussian mixture model (GMM) classifier. GMM are used for classifying hyperspectral images. The algorithm selects iteratively spectral features that maximizes an estimation of the classification rate. The estimation is done using the k-fold cross validation (k-CV). In order to perform fast in terms of computing time, an efficient implementation is proposed. First, the GMM can be updated when the estimation of the classification rate is computed, rather than re-estimate the full model. Secondly, using marginalization of the GMM, submodels can be directly obtained from the full model learned with all the spectral features. Experimental results for two real hyperspectral data sets show that the method performs very well in terms of classification accuracy and processing time. Furthermore, the extracted model contains very few spectral channels.

Description: Band selection is an essential step toward effective and efficient hyperspectral image classification. Traditional supervised band selection methods are often hindered by the problem of lacking enough training samples. To address this problem, we propose a semisupervised band selection method that allows contribution from both labeled and unlabeled hyperspectral pixels. This method first builds a hypergraph model from all hyperspectral samples to measure the similarity among pixels. We show that hypergraph can capture relationship among pixels in both spectral and spatial domain. In the second step, a semisupervised learning method is introduced to propagate class labels to unlabeled samples. Then a linear regression model with group sparsity constraint is used for band selection. Finally, hyperspectral pixels with selected bands are used to train a support vector machine (SVM) classifier. The proposed method is tested on three benchmark datasets. Experimental results demonstrate its advantages over several other band selection methods.

Description: Spectral unmixing is an important task for remotely sensed hyperspectral data exploitation. It amounts the identification of pure spectral signatures ( endmembers ) in the data, and the estimation of the abundance of each endmember in each (possibly mixed) pixel. A challenging problem in spectral unmixing is how to determine the number of endmembers in a given scene. One of the most popular and widely used techniques for this purpose is the HySime algorithm but, due to the complexity and high dimensionality of hyperspectral scenes, this technique is computational expensive. Reconfigurable field-programmable gate arrays (FPGAs) are promising platforms that allow hardware/software codesign and the potential to provide powerful onboard computing capabilities and flexibility at the same time. In this paper, we present the first FPGA design for the HySime algorithm. Our system includes a direct memory access (DMA) module and implements a prefetching technique to hide the latency of the input/output communications. The proposed method has been implemented on a Virtex-7 XC7VX690T FPGA and tested using real hyperspectral data collected by NASAs airborne visible infra-red imaging spectrometer (AVIRIS) over the Cuprite mining district in Nevada and the World trade center (WTC) in New York. Experimental results demonstrate that our hardware version of the HySime algorithm can significantly outperform a software version, which makes our reconfigurable system appealing for onboard hyperspectral data processing.

Description: Spatial–spectral classification is a very important topic in the field of remotely sensed hyperspectral imaging. In this work, we develop a parallel implementation of a novel supervised spectral–spatial classifier, which models the likelihood probability via ${bm{l}_{mathbf{1}}} - {bm{l}_{mathbf{2}}}$ sparse representation and the spatial prior as a Gibbs distribution. This classifier takes advantage of the spatial piecewise smoothness and correlation of neighboring pixels in the spatial domain, but its computational complexity is very high which makes its application to time-critical scenarios quite limited. In order to improve the computational efficiency of the algorithm, we optimized its serial version and developed a parallel implementation for commodity graphics processing units (GPUs). Our parallel spatial–spectral classifier with sparse representation and Markov random fields (SSC-SRMRF-P) exploits the low-level architecture of GPUs. The parallel optimization of the proposed method has been carried out using the compute unified device architecture (CUDA). The performance of the parallel implementation is evaluated and compared with the serial and multicore implementations on central processing units (CPUs). In fact, the proposed method has been designed to adequately exploit the massive data parallel capacities of GPUs together with the control and logic capacities of CPUs, thus resorting to a heterogeneous CPU–GPU framework in the design of the parallel algorithm. Experimental results using real hyperpsectral images demonstrate very high performance for the proposed CPU–GPU parallel method, both in terms of classification accuracy and computational performance.

Description: Classification is one of the most important analysis techniques for hyperspectral image analysis. Sparse representation is an extremely powerful tool for this purpose, but the high computational complexity of sparse representation-based classification techniques limits their application in time-critical scenarios. To improve the efficiency and performance of sparse representation classification techniques for hyperspectral image analysis, this paper develops a new parallel implementation on graphics processing units (GPUs). First, an optimized sparse representation model based on spatial correlation regularization and a spectral fidelity term is introduced. Then, we use this approach as a case study to illustrate the advantages and potential challenges of applying GPU parallel optimization principles to the considered problem. The first GPU optimization algorithm for sparse representation classification (SRCSC_P) of hyperspectral images is proposed in this paper, and a parallel implementation of the proposed method is developed using compute unified device architecture (CUDA) on GPUs. The GPU parallel implementation is compared with the serial and multicore implementations on CPUs. Experimental results based on real hyperspectral datasets show that the average speedup of SRCSC_P is more than $mathbf{130} times$ , and the proposed approach is able to provide results accurately and fast, which is appealing for computationally efficient hyperspectral data processing.

Description: This paper reports the outcomes of the 2014 Data Fusion Contest organized by the Image Analysis and Data Fusion Technical Committee (IADF TC) of the IEEE Geoscience and Remote Sensing Society (IEEE GRSS). As for previous years, the IADF TC organized a data fusion contest aiming at fostering new ideas and solutions for multisource remote sensing studies. In the 2014 edition, participants considered multiresolution and multisensor fusion between optical data acquired at 20-cm resolution and long-wave (thermal) infrared hyperspectral data at 1-m resolution. The Contest was proposed as a double-track competition: one aiming at accurate landcover classification and the other seeking innovation in the fusion of thermal hyperspectral and color data. In this paper, the results obtained by the winners of both tracks are presented and discussed.

Description: As current and future satellite systems provide both hyperspectral and multispectral images, a need has arisen for image fusion using hyperspectral and multispectral images to improve the fusion quality. This study introduces a hyperspectral image fusion algorithm using multispectral images with a higher spatial resolution and partially different wavelength range compared with the corresponding hyperspectral images. This study focuses on an image fusion technique that enhances the spatial quality and preserves the spectral information of hyperspectral images. The proposed algorithm generates a simulated multispectral band via a spectral unmixing technique and extracts high-frequency information based on blocks of associated bands. The algorithm was applied to Compact Airborne Spectrographic Imager (CASI) datasets acquired in two modes and was compared with two existing methods. Although the wavelength range of the multispectral image did not coincide with that of the hyperspectral image, the proposed algorithm efficiently improved the spatial details and preserved the spectral information of the fused results.

Description: Restoration is important in preprocessing hyperspectral images (HSI) to improve their visual quality and the accuracy in target detection or classification. In this paper, we propose a new low-rank spectral nonlocal approach (LRSNL) to the simultaneous removal of a mixture of different types of noises, such as Gaussian noises, salt and pepper impulse noises, and fixed-pattern noises including stripes and dead pixel lines. The low-rank (LR) property is exploited to obtain precleaned patches, which can then be better clustered in our spectral nonlocal method (SNL). The SNL method takes both spectral and spatial information into consideration to remove mixed noises as well as preserve the fine structures of images. Experiments on both synthetic and real data demonstrate that LRSNL, although simple, is an effective approach to the restoration of HSI.

Description: Extreme learning machine (ELM) is an efficient learning algorithm that has been recently applied to hyperspectral image classification. In this paper, the first implementation of the ELM algorithm fully developed for graphical processing unit (GPU) is presented. ELM can be expressed in terms of matrix operations so as to take advantage of the single instruction multiple data (SIMD) computing paradigm of the GPU architecture. Additionally, several techniques like the use of ensembles, a spatial regularization algorithm, and a spectral–spatial classification scheme are applied and projected to GPU in order to improve the accuracy results of the ELM classifier. In the last case, the spatial processing is based on the segmentation of the hyperspectral image through a watershed transform. The experiments are performed on remote sensing data for land cover applications achieving competitive accuracy results compared to analogous support vector machine (SVM) strategies with significantly lower execution times. The best accuracy results are obtained with the spectral–spatial scheme based on applying watershed and a spatially regularized ELM.

Description: Leaf area index (LAI) is a basic quantity indicating crop growth situation and plays a significant role in ecological model and interaction model between earth surface and atmosphere. However, nonlinear estimation processes of LAI from heterogeneous remote sensing data would induce a scaling bias. The purpose of this study is to provide a method to evaluate and correct the scaling bias. For the effectiveness of the method, first both statistical and physical models were built to estimate LAI directly from modified soil-adjusted vegetation index (MSAVI) as a function with univariate and also from red and near infrared reflectances as a bivariate model. The analysis of wavelet transform and fractal theory revealed that the scaling bias and the high-frequency coefficient from LAI at fine resolution decomposed by Haar wavelet were fractal relation. Based on the wavelet–fractal method, scaling bias could be well denoted by high-frequency coefficient in log–log coordinate for both univariate model and bivariate model, and the root-mean-square error (RMSE) and relative error (RE) of estimated LAI caused by the scaling bias could be greatly reduced after scaling correction. Additionally, to analyze the influence of spatial heterogeneity and nonlinearity of the retrieval model, the scaling bias was investigated on horizontal comparison of LAI retrieval models with univariate and bivariate at a certain resolution and was longitudinally discussed in a retrieval model at different aggregation scales. This study suggests that it is feasible to successfully correct and analyze the scaling bias using the wavelet–fractal method.

Description: Hyperspectral unmixing is an important technique for remotely sensed hyperspectral data exploitation. Linear spectral unmixing is frequently used to characterize mixed pixels in hyperspectral data. Over the last few years, many techniques have been proposed for identifying pure spectral signatures (endmembers) in hyperspectral images. The iterated constrained endmembers (ICE) algorithm is an iterative method that uses the linear model to extract endmembers and abundances simultaneously from the data set. This approach does not necessarily require the presence of pixels in the hyperspectral image as it can automatically derive the signatures of endmembers even if these signatures are not present in the data. As it is the case with other endmember identification algorithms, ICE suffers from high computational complexity. In this paper, a complete and scalable adaptation of the ICE algorithm is implemented using the parallel nature of commodity graphics processing units (GPUs). This gives significant speed increase over the traditional ICE method and allows for processing of larger data set with an increased number of endmembers.

Description: Detecting and mapping plant invaders using hyperspectral remote sensing is necessary in mitigating the extensive ecologic and economic damage these alien plants induce on our forest ecosystems. Using AISA Eagle image data, this study investigated the capability of two unsupervised classification methods for the detection and mapping of Solanum mauritianum located within commercial forestry ecosystems. The existing random forest (RF) outlier detection method when used in conjunction with Anselins Moran’s I produced a detection rate (DR) of 89% with a false positive rate (FPR) of 9.26%. In comparison, the newly developed methodology which is based on the decomposition of the RF proximity matrix using principal component analysis (PCA) resulted in a DR of 95% with a lower FPR (6.39%). Overall, this research has demonstrated the potential of utilizing an unsupervised and accurate RF framework for the detection and mapping of alien invasive plants.

Description: Hyperspectral image classification has been an active topic of research. In recent years, it has been found that light detection and ranging (LiDAR) data provide a source of complementary information that can greatly assist in the classification of hyperspectral data, in particular when it is difficult to separate complex classes. This is because, in addition to the spatial and the spectral information provided by hyperspectral data, LiDAR can provide very valuable information about the height of the surveyed area that can help with the discrimination of classes and their separability. In the past, several efforts have been investigated for fusion of hyperspectral and LiDAR data, with some efforts driven by the morphological information that can be derived from both data sources. However, a main challenge for the learning approaches is how to exploit the information coming from multiple features. Specifically, it has been found that simple concatenation or stacking of features such as morphological attribute profiles (APs) may contain redundant information. In addition, a significant increase in the number of features may lead to very high-dimensional input features. This is in contrast with the limited number of training samples often available in remote-sensing applications, which may lead to the Hughes effect. In this work, we develop a new efficient strategy for fusion and classification of hyperspectral and LiDAR data. Our approach has been designed to integrate multiple types of features extracted from these data. An important characteristic of the presented approach is that it does not require any regularization parameters, so that different types of features can be efficiently exploited and integrated in a collaborative and flexible way. Our experimental results, conducted using a hyperspectral image and a LiDAR-derived digital surface model (DSM) collected over the University of Houston campus and the neighboring urban area, indicate that the proposed fram- work for multiple feature learning provides state-of-the-art classification results.

Description: The availability of graphics processing units (GPUs) provides a low-cost solution to real-time processing, which may benefit many remote sensing applications. In this paper, a spectral–spatial classification scheme for hyperspectral images is specifically adapted for computing on GPUs. It is based on wavelets, extended morphological profiles (EMPs), and support vector machine (SVM). Additionally, a preprocessing stage is used to remove noise in the original hyperspectral image. The local computation of the techniques used in the proposed scheme makes them particularly suitable for parallel processing by blocks of threads in the GPU. Moreover, a block-asynchronous updating process is applied to the EMP to speedup the morphological reconstruction. The results over different hyperspectral images show that the execution can be speeded up to $8.2times$ compared to an efficient OpenMP parallel implementation, achieving real-time hyperspectral image classification while maintaining the high classification accuracy values of the original classification scheme.

Description: This study is an application of the hybrid regression method to improve the accuracy of retrieved humidity profile from infrared hyperspectral sounding observations. In hybrid regression method, a weighted average of two regression products that are derived using two different forms of predictand is computed. Regression coefficients for each form of predictand are computed using principal components of MetOp-IASI radiance spectra. First regression product uses logarithm of specific humidity as predictand, whereas second regression product uses only specific humidity as predictand. The weights used in hybrid regression are computed at different pressure levels based on error statistics of humidity retrieval from different predictands. The hybrid regression-based method shows improvement over the state-of-the-art regression method. Humidity profiles retrieved from different regression methods are validated with collocated ECMWF humidity profiles and radiosonde observations for dry, wet, and combined atmospheric conditions. For all cases, humidity retrieved from hybrid regression method is found to be the most accurate at all pressure levels.

Description: Target and anomaly detection are important techniques for remotely sensed hyperspectral data interpretation. Due to the high dimensionality of hyperspectral data and the large computational complexity associated to processing algorithms, developing fast techniques for target and anomaly detection has received considerable attention in recent years. Although several high-performance architectures have been evaluated for this purpose, field programmable gate arrays (FPGAs) offer the possibility of onboard hyperspectral data processing with low-power consumption, reconfigurability and radiation tolerance, which make FPGAs a relevant platform for hyperspectral processing. In this paper, we develop a novel FPGA-based technique for efficient target detection in hyperspectral images. The proposed method uses a streaming background statistics (SBS) approach for optimizing the constrained energy minimization (CEM) and Reed-Xiaoli (RX) algorithms, which are widely used techniques for target and anomaly detection, respectively. Specifically, these two algorithms are implemented in streaming fashion on FPGAs. Most importantly, we present a dual mode that implements a flexible datapath to decide in real time which one among these two algorithms should be used, thus allowing for the dynamic adaptation of the hardware to either target detection or anomaly detection scenarios. Our experiments, conducted with several well-known hyperspectral scenes, indicate the effectiveness of the proposed implementations.

Description: A sensor aboard an unmanned aerial vehicle (UAV) is vulnerable to vibration and natural conditions such as erratic winds. Considering the difference between laboratory and vicarious environments of calibration, a vicarious calibration is closer to the real environment and is a complement to laboratory calibrations for remote sensors. The existing vicarious calibration of UAVs only uses a reflectance-based method, rather than irradiance-based method. Therefore, the error caused by aerosol-type assumptions, which is the largest uncertainty for reflectance-based method, is not considered sufficiently during vicarious calibration of UAVs. Considering the difference in the upward radiative transfer path between satellites and UAVs, we propose an improved irradiance-based method. A simulation experiment was designed to compare the relative differences between two aerosol types under different aerosol optical thicknesses (AOTs) and heights for both the reflectance-based and the improved methods. Additionally, two field-calibration campaigns, under different weather conditions, were performed to calibrate a Headwall hyperspectral imager payload on a UAV, with the help of calibration tarps and MODTRAN5 radiative transfer code. When weather conditions were unsatisfactory, the total uncertainties of the original and improved methods were c. 5.9%–6.7% and c. 2.3%–3.5%, respectively, and the uncertainties caused by aerosol-type assumption were c. 15.8%–18.7% and c. 3.5%–8.0%, respectively. The results of the simulation and field experiments verified that the improved method has higher accuracy and lower uncertainty and is more suitable for the vicarious calibration of UAV hyperspectral remote sensors.

Description: Nature Climate Change 5 790 doi: 10.1038/nclimate2762

Description: A much-anticipated 'monster' El Niño failed to materialize in 2014, whereas an unforeseen strong El Niño is developing in 2015. El Niño continues to surprise us, despite decades of research into its causes. Natural variations most probably account for recent events, but climate change may also have played a role. Nature Climate Change 5 791 doi: 10.1038/nclimate2775

Description: Nature Climate Change 5 803 doi: 10.1038/nclimate2785

Description: Nature Climate Change 5 802 doi: 10.1038/nclimate2778

Description: Connected and automated vehicles enable new business models, such as self-driving taxis, that could transform transportation. These models have the potential to reduce energy consumption and greenhouse-gas emissions, but only if they are developed with energy use in mind. Nature Climate Change 5 804 doi: 10.1038/nclimate2700

Description: Some countries have pledged to become carbon neutral, while others' emissions continue to rise. Differences in their political attributes could explain the discrepancy in ambitions. Nature Climate Change 5 806 doi: 10.1038/nclimate2764

Description: Here it is argued that air pollution over West African cities needs greater consideration. The effects of aerosol pollution on clouds and solar and thermal radiation can be expected to alter regional climate and impact human health and food security. Nature Climate Change 5 815 doi: 10.1038/nclimate2727

Description: Mixed crop and livestock farms are the backbone of African agriculture, yet there is little information on how these systems may be affected by changes in climate. Addressing this knowledge gap could help smallholders adapt to climate change. Nature Climate Change 5 830 doi: 10.1038/nclimate2754

Description: Nature Climate Change 5 811 doi: 10.1038/nclimate2789

Description: Nature Climate Change 5 788 doi: 10.1038/nclimate2761

Description: Nature Climate Change 5 803 doi: 10.1038/nclimate2788

Description: The UN's climate negotiation process is no longer the 'only show in town', but there is little agreement among particpants on alternatives to replace it. Nature Climate Change 5 805 doi: 10.1038/nclimate2767

Description: A key element of the West African monsoon is changing faster than in the surrounding areas but the reason is unknown. Now research assesses the specific behaviour of the temperature over the Saharan desert. Nature Climate Change 5 807 doi: 10.1038/nclimate2773

Description: Projections of African ecological responses to climate change diverge widely. This Perspective unpicks some of the reasons for this uncertainty and reveals the importance of accounting for the influences of disturbancesand climate on vegetation. Nature Climate Change 5 823 doi: 10.1038/nclimate2753

Description: Water, energy and food security in southern Africa are interdependent and exposed to the climate. This Review considers the extent to which spatial and sectoral interdependencies can be, and are being, considered. Nature Climate Change 5 837 doi: 10.1038/nclimate2735

Description: Autonomous vehicles move passengers without human intervention. Modelling suggests that autonomous taxis could reduce transport emissions by 87–94% per mile in 2030 and save approximately 7 billion barrels of oil. Nature Climate Change 5 860 doi: 10.1038/nclimate2685

Description: Consumers in Germany are much more likely to purchase expensive ‘green’ energy produced from renewable resources if they have to actively opt out if they do not want it. In absence of such a ‘nudge’, behaviour depends more on political allegiance. Nature Climate Change 5 868 doi: 10.1038/nclimate2681

Description: Bioclimatic modelling suggests that as species distributions shift in response to climate change, few currently isolated but closely related species are likely to come into contact, implying that hybridization and competition risks will remain small. Nature Climate Change 5 883 doi: 10.1038/nclimate2699

Description: A combination of retreating sea ice and different rates of warming in the Greenland and Iceland seas is reducing winter air–sea heat fluxes. These fluxes drive ocean convection and are projected to decrease further. Nature Climate Change 5 877 doi: 10.1038/nclimate2688

Description: Shifts in the growth rate of a model green alga cultured in the presence of one or a combination of up to eight environmental drivers can generally be explained by the response to a single dominant driver, such as pH or temperature. Nature Climate Change 5 892 doi: 10.1038/nclimate2682

Description: Rapid climate warming has been linked to increasing shrub dominance in the Arctic tundra. Research now shows that climate–shrub growth relationships vary spatially and according to site characteristics such as soil moisture and shrub height. Nature Climate Change 5 887 doi: 10.1038/nclimate2697

Description: This Review looks at the state of knowledge on the El Niño/Southern Oscillation (ENSO), a natural climate phenomenon. It discusses recent advances and insights into how climate change will affect this natural climate varibility cycle. Nature Climate Change 5 849 doi: 10.1038/nclimate2743

Description: Survey data shows that policymakers are starting to seriously consider alternative climate governance forums to the United Nations Framework Convention on Climate Change. Nature Climate Change 5 864 doi: 10.1038/nclimate2684

Description: Nature Climate Change 5 803 doi: 10.1038/nclimate2787

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.