Search

Craters are among the most abundant features on the surface of many planets with great importance for planetary scientists. They reveal chronology information about planets and may be used for autonomous spacecraft navigation and landing. Although numerous research efforts have been carried out in the field of crater detection, existing crater detection algorithms (CDAs) are only helpful in a limited number of applications. A promising crater detection approach involves two main steps: 1) hypothesis generation (HG) and 2) hypothesis verification (HV). During HG, potential crater locations are detected. The validity of the hypothesized crater locations is then tested in a HV step. In this context, we discuss some commonly used algorithms for HG such as highlight-shadow region detection and Hough transform as well as our novel and enhanced algorithms based on interest point detection and convex grouping. A key objective of this paper is to analyze their performance while paying special attention to how they affect the accuracy of the verification step. To deal with different size craters, we focus on multiscale HG. For HV, we have chosen convolutional neural networks which have recently achieved state-of-the-art performance in many computer vision applications. Due to the variation of test sets in the literature, it is often challenging to compare the performance of different CDAs in a fair way. In this paper, we present a comprehensive performance evaluation and comparison of CDAs. Each algorithm has been trained/tested using common data sets generated by a systematic approach. [ABSTRACT FROM AUTHOR]

In this paper, we consider the development of algorithms for the automatic detection of buried threats using ground penetrating radar (GPR) measurements. GPR is one of the most studied and successful modalities for automatic buried threat detection (BTD), and a large variety of BTD algorithms have been proposed for it. Despite this, large-scale comparisons of GPR-based BTD algorithms are rare in the literature. In this paper, we report the results of a multi-institutional effort to develop advanced BTD algorithms for a real-world GPR BTD system. The effort involved five institutions with substantial experience with the development of GPR-based BTD algorithms. In this paper, we report the technical details of the advanced algorithms submitted by each institution, representing their latest technical advances, and many state-of-the-art GPR-based BTD algorithms. We also report the results of evaluating the algorithms from each institution on the large experimental data set used for development. The experimental data set comprised 120 000 m2 of GPR data using surface area, from 13 different lanes across two U.S. test sites. The data were collected using a vehicle-mounted GPR system, the variants of which have supplied data for numerous publications. Using these results, we identify the most successful and common processing strategies among the submitted algorithms, and make recommendations for GPR-based BTD algorithm design. [ABSTRACT FROM AUTHOR]

Linear spectral unmixing is the practice of decomposing the mixed pixel into a linear combination of the constituent endmembers and the estimated abundances. This paper focuses on unsupervised spectral unmixing where the endmembers are unknown a priori. Conventional approaches use either geometrical- or statistical-based approaches. In this paper, we address the challenges of spectral unmixing with unsupervised deep learning models, in specific, the autoencoder models, where the decoder serves as the endmembers and the hidden layer output serves as the abundances. In several recent attempts, part-based autoencoders have been designed to solve the unsupervised spectral unmixing problem. However, the performance has not been satisfactory. In this paper, we first discuss some important findings we make on issues with part-based autoencoders. By proof of counterexample, we show that all existing part-based autoencoder networks with nonnegative and tied encoder and decoder are inherently defective by making these inappropriate assumptions on the network structure. As a result, they are not suitable for solving the spectral unmixing problem. We propose a so-called untied denoising autoencoder with sparsity, in which the encoder and decoder of the network are independent, and only the decoder of the network is enforced to be nonnegative. Furthermore, we make two critical additions to the network design. First, since denoising is an essential step for spectral unmixing, we propose to incorporate the denoising capacity into the network optimization in the format of a denoising constraint rather than cascading another denoising preprocessor in order to avoid the introduction of additional reconstruction error. Second, to be more robust to the inaccurate estimation of a number of endmembers, we adopt an $l_{21}$ -norm on the encoder of the network to reduce the redundant endmembers while decreasing the reconstruction error simultaneously. The experimental results demonstrate that the proposed approach outperforms several state-of-the-art methods, especially for highly noisy data. [ABSTRACT FROM AUTHOR]

To advance the utility of polarimetric synthetic aperture radar observations for extracting biophysical information about vegetation canopies, it is important to develop good understanding of the scattering mechanisms that give rise to the polarimetric scattering response and to apply it toward the development of effective decomposition models of the polarimetric covariance matrices. In this paper, we introduce a more rigorous approach to characterizing the volume scattering component of the three-component scattering model developed by Freeman and Durden. The improved rigor has two aspects: 1) the use of a T-matrix model for computing the polarimetric scattering by an inclined dielectric cylinder and 2) the incorporation of the change in local incidence angle as a function of the orientation of the vegetation cylinders. The improvements address some of the current limitations of the Freeman–Durden model. This paper also provides a sensitivity analysis of the various components of the covariance matrix as a function of several physical parameters including cylinder size, its dielectric constant, and the orientation distribution of cylinders. Such an analysis is a precursor to the development of improved inversion algorithms. [ABSTRACT FROM AUTHOR]

Capon’s minimum-variance method (MVM) and amplitude and phase estimation (APES) spectral estimation algorithms can be applied to synthetic aperture radar (SAR) processing to improve the resolution and suppress sidelobe levels. In this paper, we use Capon-/APES-based SAR reprocessing algorithms to increase the persistent scatterer (PS) density in PS interferometry (PSI). We propose a PS candidate (PSC) selection algorithm applicable to the superresolution reprocessed images and the corresponding processing chain. The performance of the proposed algorithm is evaluated by a number of simulations and a stack of TerraSAR-X data. The results show that the Capon algorithm outperforms others in PSC selection. We present a full PSI time-series analysis on the PSCs extracted from the Capon-reprocessed stacks. The results show that the PS density is increased between 50% and 60%, while their interferometric quality is maintained. [ABSTRACT FROM AUTHOR]

In most of the existing ground penetrating radar (GPR) imaging algorithms, using either full or sparse data collection, the ground is modeled as a single half-space layer. In this paper, a generalized sparse imaging approach with total variation minimization (TVM) for multiple-input multiple-output (MIMO) GPR imaging through multilayered subsurface is proposed. The multilayered media Green’s function is incorporated in the imaging algorithm to take into account the complex wave propagation effects under multilayered subsurface. An analytical expression of the layered subsurface Green’s function is derived using the stationary-phase method, which significantly reduces the computation time and complexity. On the other hand, as TVM minimizes the gradient of the image, its incorporation in the imaging algorithm results in a reconstruction that preserves the geometry and edges of the targets better than the standard L1-minimization-based sparsity-driven imaging. The number of antenna elements and frequency measurements in MIMO GPR system can be significantly reduced using the proposed technique without degradation of the image quality. Although MIMO configuration is investigated in this paper, the presented approach can be simply applied to monostatic synthetic aperture radar. [ABSTRACT FROM AUTHOR]

Sea ice concentration (SIC) is an important sea ice parameter for sea ice navigation, environmental research, and weather and ice forecasting. We have developed and tested a method for estimation of the Baltic Sea SIC using SENTINEL-1 synthetic aperture radar (SAR) and Advanced Microwave Scanning Radiometer 2 passive microwave radiometer (MWR) data. Here, we present the method and results for January 2016. Ice concentration grids of Finnish Meteorological Institute daily ice charts have been used as reference data in this paper. We present a comparison of four SIC estimation methods with our reference data. In addition to the combined SAR/MWR SIC estimation method, we also compare SIC estimates produced using SAR alone and two MWR-based methods. The main target of this paper was to develop and test a high-resolution SIC estimation method suitable for operational use. [ABSTRACT FROM PUBLISHER]

High-resolution thermal infrared (TIR) remote sensing has a wide range of applications. In this paper, we describe the different applications and requirements identified in a literature review and during a consultation meeting with researcher experts in different fields. As a result, more than 30 applications were identified within three different fields: 1) land and solid Earth; 2) health and hazards; and 3) security and surveillance. A complete set of requirements (spatial, temporal, and radiometric resolution, algorithms used, and supporting data, among others) for each application is also provided. The results presented in this paper provide useful information to enhance the importance of high-resolution TIR data for civil applications and may serve as a reference document for future TIR mission concepts. [ABSTRACT FROM AUTHOR]

Surface roughness is a land-surface parameter that is widely used in terrain analysis. Some typical roughness details, which have important effects on surface analysis, fail to be characterized on previous roughness maps. The objective of this paper is to provide a more accurate small-to-large scale roughness overview. The new roughness method is designed based on a complete orthogonal system called the V-system. The V-system roughness utilizes the special functions to detect and extract the roughness characteristics from high-resolution digital elevation models (DEMs). In this paper, Lunar Orbiter Laser Altimeterderived DEMs are used as the source data for the roughness calculation. Compared with the global root-mean-square slope and Fourier-based roughness maps, the V-system roughness maps show that more typical roughness details have been added to clearly indicate the small roughness variations on the large map. Furthermore, the reliability and practicability of V-system roughness are demonstrated based on the multiresolution DEMs. As an example, the statistical parameters of the roughness characteristics in the lunar Maria and highlands identify the fact that the highlands are rougher at all scales than the Maria. And this difference corresponds to the basic roughness property. [ABSTRACT FROM AUTHOR]

The automatic detection of buried threats in ground penetrating radar (GPR) data is an active area of research due to GPR's ability to detect both metal and nonmetal subsurface objects. Recent work on algorithms designed to distinguish between threats and nonthreats in GPR data has utilized computer vision methods to advance the state-of-the-art detection and discrimination performance. Feature extractors, or descriptors, from the computer vision literature have exhibited excellent performance in representing 2-D GPR image patches and allow for robust classification of threats from nonthreats. This paper aims to perform a broad study of feature extraction methods in order to identify characteristics that lead to improved classification performance under controlled conditions. The results presented in this paper show that gradient-based features, such as the edge histogram descriptor and the scale invariant feature transform, provide the most robust performance across a large and varied data set. These results indicate that various techniques from the computer vision literature can be successfully applied to target detection in GPR data and that more advanced techniques from the computer vision literature may provide further performance improvements. [ABSTRACT FROM PUBLISHER]

Soil Moisture and Ocean Salinity (SMOS) mission has successfully contributed to global soil moisture products since 2009. Validation and calibration activities were conducted worldwide, yet some of the validation results do not fulfill the targeted accuracy of ±0.04 $\text{m}^{3}\text{m}^{-3}$. This paper presented the site-specific calibration of the V620 retrieval algorithm with in situ data collected at selected agricultural sites in the humid tropical regions, Malaysia. This set of data has been validated where low accuracy of SMOS soil moisture products was found. To improve the SMOS soil moisture retrieval, calibration of SMOS soil moisture retrieval algorithm based on the L-band Microwave Emission and Biosphere model and SMOS Level 1C $\text{T}_{\mathrm {B}}$ products, considering the local parameters was conducted. The calibration proves that these site-specific parameters improve the product’s accuracy. Validation of SMOS Level 2 product with in situ data showed bias, root-mean-square error (RMSE), and unbiased RMSE (ubRMSE) ranging from 0.050 to 0.118 $\text{m}^{3}\text{m}^{-3}$ , 0.068 to 0.142 $\text{m}^{3}\text{m}^{-3}$ , and 0.069 to 0.103 $\text{m}^{3}\text{m}^{-3}$ , respectively. The soil moisture retrieval based on the calibrated model showed an improved bias of 0.020–0.056 $\text{m}^{3}\text {m}^{-3}$ and RMSE of 0.026–0.065 $\text{m}^{3}\text{m}^{-3}$. The ubRMSE ranges from 0.017 to 0.034 $\text{m}^{3}\text{m}^{-3}$. Recently released SMOS-IC V105 product was also validated, where small improvements were noticed when compared to the accuracy of SMOS Level 2. This paper shows the importance of local parameters in retrieving soil moisture with higher accuracy compared to the use of global generalized parameters that are used in the original SMOS soil moisture retrieval algorithm. [ABSTRACT FROM AUTHOR]

For the quantitative applications of the Suomi National Polar-orbiting Partnership (SNPP) Advanced Technology Microwave Sounder (ATMS), the geolocation accuracy of its sensor data records must be quantified during its on-orbit operation. In this paper, a refined coastline inflection point method is used to evaluate the on-orbit geolocation accuracy of SNPP ATMS. It is disclosed that for SNPP ATMS, the static error term with scan-angle-dependent feature is a dominant part among all the geolocation error sources. A mathematical model is then developed to convert the in-track and cross-track geolocation errors to the beam pointing Euler angles defined in the spacecraft coordinate system, which can be further used to construct the correction matrix for on-orbit geolocation process. By using the correction matrix built in this paper, the geolocation error is obviously reduced both at nadir and at the edge of the scan. The total geolocation error at nadir before/after correction is 3.8/0.8 km at K-band, 5.6/0.8 km at Ka-band, 3.3/0.4 km at V-band, and 1.5/0.1 km at W-band. The geolocation bias at the edge of the scan line before/after correction is 4.6/1.3 km at K-band, 9.4/1.8 km at Ka-band, 4.4/2.4 km at V-band, and 3.2/0.8 km at W-band. After correction, the scan-angle-dependent feature in geolocation error is also largely reduced. [ABSTRACT FROM AUTHOR]

Band selection is an effective approach to mitigate the “Hughes phenomenon” of hyperspectral image (HSI) classification. Recently, sparse representation (SR) theory has been successfully introduced to HSI band selection, and many SR-based methods have been developed and shown great potential and superiority. However, due to the inherent limitations of the SR scheme, i.e., individually representing each band with only a few other bands from the same subspace, the SR-based methods cannot effectively capture the global structures of the data, which limit the band selection performance. In this paper, to overcome this obstacle, the novel Laplacian-regularized low-rank subspace clustering (LLRSC) algorithm is proposed for HSI band selection. On the one hand, the low-rank subspace clustering model is introduced to capture the global structure information for the learned representation coefficient matrix and deal with the HSI band selection task in the clustering framework. On the other hand, considering the high correlation between adjacent bands, 1-D Laplacian regularization is utilized to incorporate the neighboring band information and further reduce the representation bias. Lastly, an eigenvalue analysis algorithm based on band mutation information is utilized to estimate the appropriate size of the band subset. The experimental results indicate that the proposed LLRSC algorithm outperforms the other state-of-the-art methods and achieves a very competitive band selection performance for HSIs. [ABSTRACT FROM AUTHOR]

Accurate and high throughput extraction of crop phenotypic traits, as a crucial step of molecular breeding, is of great importance for yield increasing. However, automatic stem–leaf segmentation as a prerequisite of many precise phenotypic trait extractions is still a big challenge. Current works focus on the study of the 2-D image-based segmentation, which are sensitive to illumination and occlusion. Light detection and ranging (LiDAR) can obtain accurate 3-D information with its active laser scanning and strong penetration ability, which breaks through phenotyping from 2-D to 3-D. However, few researches have addressed the problem of the LiDAR-based stem–leaf segmentation. In this paper, we proposed a median normalized-vector growth (MNVG) algorithm, which can segment stem and leaf with four steps, i.e., preprocessing, stem growth, leaf growth, and postprocessing. The MNVG method was tested by 30 maize samples with different heights, compactness, leaf numbers, and densities from three growing stages. Moreover, phenotypic traits at leaf, stem, and individual levels were extracted with the truly segmented instances. The mean accuracy of segmentation at point level in terms of the recall, precision, F-score, and overall accuracy were 0.92, 0.93, 0.92, and 0.93, respectively. The accuracy of phenotypic trait extraction in leaf, stem, and individual levels ranged from 0.81 to 0.95, 0.64 to 0.97, and 0.96 to 1, respectively. To our knowledge, this paper proposed the first LiDAR-based stem–leaf segmentation and phenotypic trait extraction method in agriculture field, which may contribute to the study of LiDAR-based plant phonemics and precise agriculture. [ABSTRACT FROM AUTHOR]

This paper introduces a bistatic fast backprojection synthetic aperture radar (SAR) algorithm that is available for different bistatic geometries. The proposed algorithm is tested with the bistatic CARABAS-like data, and the results indicate that the algorithm is a good candidate for bistatic SAR data processing. The image quality measurements are quite similar to the referenced values obtained with the bistatic global backprojection algorithm. That is, the peak sidelobe ratio (PSLR) is −13.7 dB in comparison to the referenced PSLR of −13.8 dB. The half-power beamwidths (HPBWs) measured on the cuts in $x$ and $y$ and the direction where the peak sidelobes locate are 2.8, 3.7, and 3.5 m, respectively, which are approximate to the referenced HPBWs. The small differences in the measured results mainly come from the interpolation step of the algorithm. [ABSTRACT FROM AUTHOR]

This paper deals with signal reconstruction for multichannel high-resolution wide-swath synthetic aperture radar (HRWS-SAR) systems. One of the commonly used algorithms for processing HRWS-SAR data is Krieger's digital beamforming (DBF) algorithm. This algorithm works quite well for cases with uniform space–time sampling as well as moderately nonuniform sampling. However, for cases with highly nonuniform sampling, Krieger's DBF algorithm in its original form suffers from a severe degradation in the signal-to-noise ratio, and for cases with coinciding sampling, this algorithm even fails to function. In this paper, we will solve these problems such that an improved DBF (IDBF) algorithm is developed. The key of the IDBF algorithm consists of two procedures: The first one implements an operation called ambiguity-index screening, which is based on two novel quantities, termed the equivalent sampling spacing and the effective reconstructed Doppler bandwidth; and the other key procedure implements a Doppler-spectrum weighting based on a newly defined parameter that is termed the maximum-to-minimum ratio of the noise spectrum. Simulation results show that the newly proposed IDBF algorithm shares a pronounced improvement in the imaging performance over Krieger's DBF algorithm. Moreover, simulation results also demonstrate that the IDBF algorithm exhibits a more robust and consistent overall performance than several other existing reconstruction algorithms. [ABSTRACT FROM PUBLISHER]

In this paper, we discuss the structure and performance of the maximum entropy temperature-emissivity separation (MaxEnTES) algorithm for assessing land temperature and emissivity from thermal infrared hyperspectral images. This procedure derives the emissivity spectrum and the temperature of the target adopting the maximum entropy (MaxEnt) estimation approach. The main advantage of the MaxEnt statistical inference is the absence of any external hypothesis, which is, instead, the main critical point characterizing any other temperature-emissivity separation (TES) algorithm. The MaxEnTES algorithm carries out the TES task adopting a modified version of the subgradient Shor's r-algorithm adopted for numerical optimization of a MaxEnt objective function. For this purpose, we have utilized the C/C++ Solvopt code from the University of Gratz to develop a practical data processing implementation. In this paper, we discuss the mathematical structure of the MaxEnTES algorithm and analyze its performance in depth using numerical simulations and remote sensing Multispectral Infrared/Visible Imaging Spectrometer images. We show that the MaxEnTES algorithm provides improved accuracy for temperature and emissivity estimation, lowering the standard estimation error to a fraction of degree Kelvin. In agreement with previous investigations, we find that the estimation accuracy grows when increasing the number of available spectral channels. The systematic errors affecting the temperature estimates (e.g., bias) are thoroughly evaluated. We prove that the MaxEnTES algorithm retrieves the correct shape of the target emissivity spectrum even in presence of a significant temperature estimation error. [ABSTRACT FROM PUBLISHER]

Estimation of the surface elevation change of the Greenland Ice Sheet (GrIS) is essential for understanding its response to recent and future climate change. Laser measurements from the NASA’s Ice, Cloud, and land Elevation Satellite (ICESat) created altimetric surveys of GrIS surface elevations over the 2003–2009 operational period of the mission. This paper compares four change detection methods using Release 634 ICESat laser altimetry data: repeat tracks (RTs), crossovers (XOs), overlapping footprints (OFPs), and triangulated irregular networks (TINs). All four methods begin with a consistently edited data set and yield estimates of volumetric loss of ice from the GrIS ranging from −193 to −269 km3/yr. Using a uniform approach for quantifying uncertainties, we find that volume change rates at the drainage system scale from the four methods can be reconciled within 1- $\sigma $ uncertainties in just 5 of 19 drainage systems. Ice-sheet-wide volume change estimates from the four methods cannot be reconciled within 1- $\sigma $ uncertainties. Our volume change estimates lie within the range of previously published estimates, highlighting that the choice of method plays a dominant role in the scatter of volume change estimates. We find that for much of the GrIS, the OFP and TIN methods yield the lowest volume change uncertainties because of their superior spatial distribution of elevation change rate estimates. However, the RT and XO methods offer inherent advantages, and the future work to combine the elevation change detection methods to produce better estimates is warranted. [ABSTRACT FROM PUBLISHER]

The availability of multiangular information from the NASA Multi-angle Imaging SpectroRadiometer (MISR) instrument provides an excellent opportunity for the characterization of surface anisotropy, which can be used for improving surface albedo estimation. However, the MISR data have been reported with large uncertainties and data gaps due to inaccurate aerosol estimation and/or cloud masking limiting its otherwise broader applications. To mitigate these issues, two approaches were proposed to estimate land surface albedo directly from surface reflectance (LSA_sfc) and Top-of-Atmosphere reflectance (LSA_toa), respectively. As a further development of the traditional albedo algorithms, this is the first attempt to simultaneously utilize multispectral and multiangular information in surface albedo estimation without any prior constraining information. Validations at AmeriFlux sites show that the proposed algorithms can achieve accuracies similar to that of the MISR product with respective bias and RMSE of 0.004 and 0.032 for LSA_sfc and 0.005 and 0.032 for LSA_toa algorithms. We found that the LSA_toa algorithm can significantly reduce data gaps and provide accurate surface albedo retrievals with two to three times more valid data than the current MISR product. In addition, these approaches can be easily applied to other optical sensors to produce accurate and gap-free clear-sky surface albedo estimations. The results of this paper also highlight the importance of having two to three simultaneous observations with sufficient angular sampling, which can improve albedo accuracy and reduce data gaps. [ABSTRACT FROM PUBLISHER]

In this paper, we propose a long-time coherent integration method for a maneuvering target with the first-, second-, and third-order range migrations, and the complex Doppler frequency migration. In this method, after range compression, the echo signal is first transformed into the range-frequency and Doppler domain based on series reversion, and then an azimuth matched filtering procedure is implemented in the 2-D frequency domain. It can eliminate the coupling effects between range and azimuth jointly caused by the radial velocity, radial acceleration, and radial acceleration rate of a moving target. Due to the linear transform property, the proposed method can work well under low signal-to-clutter and noise ratio. In addition, the Doppler ambiguities, when target azimuth spectrum either is within a pulse repetition frequency (PRF) or spans over neighboring PRF bands, can be well solved. Both simulated and real synthetic aperture radar data processing results are provided to validate the effectiveness of the proposed algorithm. [ABSTRACT FROM AUTHOR]

Hyperspectral unmixing has attracted much attention in recent years. Single sparse unmixing assumes that a pixel in a hyperspectral image consists of a relatively small number of spectral signatures from large, ever-growing, and available spectral libraries. Joint-sparsity (or row-sparsity) model typically enforces all pixels in a neighborhood to share the same set of spectral signatures. The two sparse models are widely used in the literature. In this paper, we propose a joint-sparsity-blocks model for abundance estimation problem. Namely, the abundance matrix of size $m\times n$ is partitioned to have one row block and $s$ column blocks and each column block itself is joint-sparse. It generalizes both the single (i.e., $s=n$) and the joint (i.e., $s=1$) sparsities. Moreover, concatenating the proposed joint-sparsity-blocks structure and low rankness assumption on the abundance coefficients, we develop a new algorithm called joint-sparse-blocks and low-rank unmixing. In particular, for the joint-sparse-blocks regression problem, we develop a two-level reweighting strategy to enhance the sparsity along the rows within each block. Simulated and real-data experiments demonstrate the effectiveness of the proposed algorithm. [ABSTRACT FROM AUTHOR]

A newly designed grounded electrical-source airborne transient electromagnetics (GREATEM) system was introduced recently. Detailed data preprocessing techniques to acquire the high-precision measured magnetic field are discussed here. Different from the previous work in which the reconstruction of the underground structure is performed in 1-D, we interpret the GREATEM data in 3-D by the volume integral equation (VIE) method in the frequency domain. Therefore, the VIE in the forward electromagnetic scattering model is formulated in the low-frequency regime. It is solved by using the stabilized biconjugate gradient fast Fourier transform (BCGS-FFT) method. In the nonlinear inversion, the Born iterative method (BIM) and the conjugate gradient method are adopted to minimize the cost function. A synthetic model of GREATEM survey is used to validate the proposed 3-D forward and inversion algorithms. Then, the field data from two GREATEM surveys are used to test the effectiveness and accuracy of the proposed inversion algorithm. The reconstructed conductivity structures are consistent with geological drilling results, confirming the potential of our method for solving the 3-D GREATEM inversion problems in geophysical engineering applications. This paper represents the first application of the BCGS-FFT and BIM algorithms to a GREATEM system. [ABSTRACT FROM AUTHOR]

With the rapid increase of remote sensing images in temporal, spectral, and spatial resolutions, it is urgent to develop effective techniques for joint interpretation of spatial-temporal images. Multitype change detection (CD) is a significant research topic in multitemporal remote sensing image analysis, and its core is to effectively measure the difference degree and represent the difference among the multitemporal images. In this paper, we propose a novel difference representation learning (DRL) network and present an unsupervised learning framework for multitype CD task. Deep neural networks work well in representation learning but rely too much on labeled data, while clustering is a widely used classification technique free from supervision. However, the distribution of real remote sensing data is often not very friendly for clustering. To better highlight the changes and distinguish different types of changes, we combine difference measurement, DRL, and unsupervised clustering into a unified model, which can be driven to learn Gaussian-distributed and discriminative difference representations for nonchange and different types of changes. Furthermore, the proposed model is extended into an iterative framework to imitate the bottom-up aggregative clustering procedure, in which similar change types are gradually merged into the same classes. At the same time, the training samples are updated and reused to ensure that it converges to a stable solution. The experimental studies on four pairs of multispectral data sets demonstrate the effectiveness and superiority of the proposed model on multitype CD. [ABSTRACT FROM AUTHOR]

This paper presents a new method based on recent optimization technique to detect slow-moving landslides (<150m/year) in time series of displacement field generated by satellite images. Sparse optimization is applied simultaneously on the 3-D data set in space as well as in time. The proposed method takes into account the distinctive signal physical properties in space and time, by enforcing a sparse representation by blocks in space, but a continuing and monotonous representation in time of the landslides. As a result, we show that a mixed $\ell _{1,2}$ -norm is the most suitable norm for this detection problem, compared to pure $\ell _{1}$ -norm or $\ell _{2}$ -norm. Moreover, an outlier estimation step is included that sets apart the Gaussian noise from locally sparse processing errors in the data. The performance of this approach is tested by applying it both on synthetic data and on a time series of displacements fields over 16 dates in the Colca Valley, Peru. This detection presents commission and omission errors for landslides of 29% and 14%, respectively, using a medium resolution (10 m) data set of optical satellite images. It detects all important landslides, already known from field investigations. Moreover, it also points out other smaller or unknown landslides, increasing the existing slow-moving landslide inventory by +50%. [ABSTRACT FROM AUTHOR]

A rain-gauge-adjusted algorithm for global satellite mapping of precipitation (GSMaP) that estimates the surface precipitation rate with resolutions of 0.1° and 1 h over the globe is described herein. Precipitation is one of the most important parameters of the earth’s system, and its global distribution and changes are essential data for modeling the water cycle, maintaining ecosystems, increasing agricultural production, improving weather forecasting precision, and implementing flood warning systems. In the Global Precipitation Measurement project, integrated products of high-resolution mapping of precipitation, obtained from microwave measurements made by a constellation satellite and infrared radiometers in geostationary orbit, are developed and supplied to the public. However, these high-resolution products, such as GSMaP_MVK, sometimes underestimate surface precipitation, introducing large errors into hydrological modeling. This paper combines the global gauge data set with GSMaP_MVK, using a new algorithm [gauge-adjusted GSMaP (GSMaP_Gauge)], described and evaluated herein using local radar and rain-gauge data sets. This algorithm outperforms other GSMaP products in all validation tests. [ABSTRACT FROM AUTHOR]

Despeckling is a fundamental preprocessing step for applications using polarimetric synthetic aperture radar data in most cases. In this paper, a guided filter with nonlinear weight kernels and adaptive filtering windows is presented for PolSAR image despeckling, in which the guidance image is constructed by a weighted average using the statistical traits of the speckled image. The output result is then estimated by another weighted average, with the aid of the fully polarimetric information from both the guidance image and the speckled image. In the experimental part, the filtering results obtained with both simulated and real PolSAR images reveal the positive performance of the proposed method in both reducing speckle and retaining details, when compared with some of the state-of-the-art algorithms. Furthermore, the relatively low computational complexity is another strength of the proposed method. [ABSTRACT FROM AUTHOR]

In this paper, we investigate the temperature and emissivity separation (TES) problem from hyperspectral data acquired in the long-wave infrared region (LWIR) of the electromagnetic spectrum. We derive a general class of TES algorithms [subspace-based TES (SBTES)] relying on the assumption that the emissivity spectra of natural and man-made materials can be well represented in a given subspace of the original data space. Specifically, by exploiting the subspace representation and the Gaussian model for the noise affecting LWIR hyperspectral data, we approach TES under a statistical perspective by obtaining the maximum likelihood estimates of both the temperature and the spectral emissivity. The proposed approach originates several algorithms whose specific form depends on the particular basis matrix adopted to address the emissivity subspace. We study the performance of the presented class of algorithms and derive theoretical bounds on the accuracy of the temperature and emissivity estimators. Furthermore, by specifying two basis matrices for the emissivity subspace, we propose two different algorithms within the SBTES class. Finally, we present the results of an extensive experimental analysis carried out over simulated data to assess and compare the performance of the two presented algorithms. [ABSTRACT FROM AUTHOR]

Establishing causal relations between random variables from observational data is perhaps the most important challenge in today’s science. In remote sensing and geosciences, this is of special relevance to better understand the earth’s system and the complex interactions between the governing processes. In this paper, we focus on an observational causal inference, and thus, we try to estimate the correct direction of causation using a finite set of empirical data. In addition, we focus on the more complex bivariate scenario that requires strong assumptions and no conditional independence tests can be used. In particular, we explore the framework of (nondeterministic) additive noise models, which relies on the principle of independence between the cause and the generating mechanism. A practical algorithmic instantiation of such principle only requires: 1) two regression models in the forward and backward directions and 2) the estimation of statistical independence between the obtained residuals and the observations. The direction leading to more independent residuals is decided to be the cause. We instead propose a criterion that uses the sensitivity (derivative) of the dependence estimator, the sensitivity criterion allows to identify samples most affecting the dependence measure, and hence, the criterion is robust to spurious detections. We illustrate the performance in a collection of 28 geoscience causal inference problems, a database of radiative transfer models simulations and machine learning emulators in vegetation parameter modeling involving 182 problems, and assessing the impact of different regression models in a carbon cycle problem. The criterion achieves the state-of-the-art detection rates in all cases, and it is generally robust to noise sources and distortions. The presented approach confirms the validity in observational bivariate problems in the earth sciences. [ABSTRACT FROM AUTHOR]

Relationships between vegetation indices (VIs) and leaf area index (LAI) tend to saturate in the nadir direction, and vary with crop canopy structure and view zenith angles (VZAs). The objective of this paper was to improve the monitoring accuracy and angular stability of VIs for estimating LAI using multiangular remote sensing data. The relationship between LAI and ground-based hyperspectral spectral reflectance was quantified in winter wheat (Triticum aestivum L.) exhibiting erectophile and planophile growth habits. To further reduce the saturation, species specificity, and angular sensitivity, we developed a saturation factor (SF), based on near-infrared and green bands. Multiplying all VIs by SF greatly improved the association with LAI across all VZAs ($R^{2} = 0.73{-}0.82$). Most VI $\times$ SF values, particularly optimized soil-adjusted VI $\times$ SF, were able to construct universal algorithms across VZAs for accurate estimation of LAI due to the sensitivity of SF to LAI in a dense canopy, and the insensitivity of SF to view effects with larger VZAs. This approach is also promising for the exploitation of multiangular satellite data for the design and calibration of nonview-angle-corrected spectral reflectance, for which the sensor is only deployed at fixed observation direction. [ABSTRACT FROM AUTHOR]

The advanced Himawari imager (AHI) aboard the Himawari-8 geostationary satellite provides high-frequency observations with broad coverage, multiple spectral channels, and high spatial resolution. In this paper, AHI data were used to develop an algorithm for joint retrieval of aerosol optical depth (AOD) over land and land surface bidirectional reflectance. Instead of performing surface reflectance estimation before calculating AOD, the AOD and surface bidirectional reflectance were retrieved simultaneously using an optimal estimation method. The algorithm uses an atmospheric radiative transfer model coupled with a surface bidirectional reflectance factor (BRF) model. Based on the assumption that the surface bidirectional reflective properties are invariant during a short time period (i.e., a day), multiple temporal AHI observations were combined to calculate the AOD and surface BRF. The algorithm was tested over East Asia for year 2016, and the AOD retrieval results were validated against the aerosol robotic network (AERONET) sites observation and compared with the Moderate Resolution Imaging Spectroradiometer Collection 6.0 AOD product. The validation of the retrieved AOD with AERONET measurements using 14 713 colocation points in 2016 over East Asia shows a high correlation coefficient: $R = 0.88$ , root-mean-square error = 0.17, and approximately 69.9% AOD retrieval results within the expected error of $\pm 0.2\cdot {\mathrm {AOD}}_{\mathrm {AERONET}}\pm 0.05$. A brief comparison between our retrieval and AOD product provided by Japan Meteorological Agency is also presented. The comparison and validation demonstrates that the algorithm has the ability to estimate AOD with considerable accuracy over land. [ABSTRACT FROM AUTHOR]

This paper presents a processing framework to separate moving targets from the clutter, under multichannel synthetic aperture radar (SAR) scenarios, and addresses the moving target imaging and velocity estimation problems for ground moving target indication (GMTI) applications. A practical implementation is introduced to break the SAR/GMTI problem into two processing stages, and the sparsity of the moving targets in the observed scene is exploited throughout the stages. The two-stage process extracts the moving targets from the monitored region via a sparsity-based iterative decomposition algorithm and subsequently estimates the complete velocity vectors of moving targets by enforcing sparsity constraints. The model is sufficiently versatile to incorporate digital elevation map information, which further improves the moving target relocation accuracy. The effectiveness of the presented framework is demonstrated using the Air Force Research Laboratory Gotcha GMTI challenge data. [ABSTRACT FROM AUTHOR]

The airborne Operation IceBridge and spaceborne CryoSat-2 missions observe polar sea ice at different spatial and temporal scales as well as with different sensor suites. Comparison of data products from IceBridge and CryoSat-2 is complicated by the fact that they use different geophysical corrections: reference ellipsoid, geoid model, tide model, and atmospheric corrections to derive surface elevation and sea-ice freeboard. In this paper, we compare sea-ice surface elevation and freeboard using eight coincident CryoSat-2, Airborne Topographic Mapper (ATM), and Land, Vegetation, and Ice Sensor (LVIS) observations with direct IceBridge underflights of CryoSat-2 ground tracks. We apply identical geophysical corrections to CryoSat-2 and IceBridge data to eliminate elevation biases due to these differences and focus on differences due to retracker performance. The IceBridge ATM and LVIS elevation and freeboard and Snow Radar snow depth data sets are averaged to each CryoSat-2 footprint for comparison. With snow depth measurements, we are able to compare elevations and freeboards at the snow/ice interface for five different CryoSat-2 retrackers (ESA, GSFCv1, AWI, JPL, and GSFCv2) and IceBridge. The overall mean of freeboard differences between GSFCv2, ESA, AWI, JPL retrackers, and ATM are in agreement within 0.05 m. However, the five different CryoSat-2 retrackers show distinct differences in mean elevation over leads and over floes. This suggests that the physical interpretation of the different retrackers needs to be considered depending on usage, for example, elevations from CryoSat-2 retrackers need to be carefully calibrated before comparing with elevation from other satellites for long-term surface elevation trends. [ABSTRACT FROM AUTHOR]

The simplex growing algorithm (SGA) has been widely used for finding endmembers. It can be considered as a sequential version of the well-known endmember finding algorithm, N-finder algorithm (N-FINDR), which finds endmembers one at a time by growing simplexes. However, one of the major hurdles for N-FINDR and SGA is the calculation of simplex volume (SV) which poses a great challenge in designing any algorithm using SV as a criterion for finding endmembers. This paper develops an orthogonal projection (OP)-based SGA (OP-SGA) which essentially resolves this computational issue. It converts the issue of calculating SV to calculating the OP on previously found simplexes without computing matrix determinants. Most importantly, a recursive Kalman filter-like OP-SGA, to be called recursive OP-SGA (ROP-SGA), can be also derived to ease computation. By virtue of ROP-SGA, several advantages and benefits in computational savings and hardware implementation can be gained for which N-FINDR and SGA do not have. [ABSTRACT FROM AUTHOR]

This paper develops and characterizes the algorithms used to generate the Level 1 (L1) science data products of the Cyclone Global Navigation Satellite System (CYGNSS) mission. The L1 calibration consists of two parts: the Level 1a (L1a) calibration converts the raw Level 0 delay–Doppler maps (DDMs) of processed counts into received power in units of watts. The L1a DDMs are then converted to Level 1b DDMs of bistatic radar cross section values by unwrapping the forward scattering model and generating two additional DDMs: one of unnormalized bistatic radar cross section values (in units of square meters) and a second of bin-by-bin effective scattering areas. The L1 data products are generated in such a way as to allow for flexible processing of variable areas of the DDM (which correspond to different regions on the surface). The application of the L1 data products to the generation of input observables for the CYGNSS Level 2 (L2) wind retrievals is also presented. This includes a demonstration of using only near-specular DDM bins to calculate a normalized bistatic radar cross section (unitless, i.e., m2/m2) over a subset of DDM pixels, or DDM area. Additionally, an extensive term-by-term error analysis has been performed using this example extent of the DDM to help quantify the sensitivity of the L1 calibration as a function of key internal instrument and external parameters in the near-specular region. [ABSTRACT FROM AUTHOR]

This paper introduces robust algorithms for extracting the ground points in laser scanning 3-D point cloud data. Global polynomial functions have been used for filtering algorithms for point cloud data; however, it is not suitable as it may lead to bias for the filtering algorithms and can cause misclassification errors when many different objects are present. In this paper, robust statistical approaches are coupled with locally weighted 2-D regression that fits without any predefined global function for the variables of interest. Algorithms are performed iteratively on 2-D profiles: $x-z$ and $y-z$. The $z$ (elevation) values are robustly down weighted based on the residuals for the fitted points. The new set of down-weighted $z$ values, along with the corresponding $x$ (or $y$) values, is used to get a new fit for the lower surface level. The process of fitting and down weighting continues until the difference between two consecutive fits is insignificant. The final fit is the required ground level, and the ground surface points are those that fall within the ground level and the level after adding some threshold value with the ground level for $z$ values. Experimental results are compared with the recently proposed segmentation method through simulated and real mobile laser scanning point clouds from urban areas that include many objects that appear in road scenes such as short walls, large buildings, electric poles, signposts, and cars. Results show that the proposed robust methods efficiently extract ground surface points with better than 97% accuracy. [ABSTRACT FROM PUBLISHER]

This paper proposes a synchronization algorithm for bistatic synthetic aperture radar (BiSAR) imaging in a spaceborne/stationary configuration. In real bistatic systems, synchronization errors are generally introduced into the received data. Additionally, the lack of precise imaging parameters, such as the position of the transmitter and the accurate sampling time, could affect the imaging quality greatly. Fortunately, the image could be well focused by the proposed algorithm in the case of lack of the accurate position of a transmitter and the sampling time. First, a preprocessing step is employed to remove synchronization errors through matching an echo signal with a direct signal. Then, a modified chirp scaling factor containing an error phase term is constructed, and the accurate position of the transmitter and the sampling time can be acquired by the phase extraction of the direct signal and the searching method based on contrast optimization. After that, the corresponding imaging process can be implemented. Finally, the proposed algorithm is validated by the simulation and experimental results, where TerraSAR-X is used as the illuminator. [ABSTRACT FROM PUBLISHER]

In blind hyperspectral unmixing (HU), the pure-pixel assumption is well known to be powerful in enabling simple and effective blind HU solutions. However, the pure-pixel assumption is not always satisfied in an exact sense, especially for scenarios where pixels are heavily mixed. In the no-pure-pixel case, a good blind HU approach to consider is the minimum volume enclosing simplex (MVES). Empirical experience has suggested that MVES algorithms can perform well without pure pixels, although it was not totally clear why this is true from a theoretical viewpoint. This paper aims to address the latter issue. We develop an analysis framework wherein the perfect endmember identifiability of MVES is studied under the noiseless case. We prove that MVES is indeed robust against lack of pure pixels, as long as the pixels do not get too heavily mixed and too asymmetrically spread. The theoretical results are supported by numerical simulation results. [ABSTRACT FROM PUBLISHER]

This paper describes a clutter suppression approach and the corresponding moving target imaging algorithm for a multichannel in azimuth high-resolution and wide-swath (MC-HRWS) synthetic aperture radar (SAR) system. Incorporated with digital beamforming processing, MC-HRWS SAR systems are able to suppress the Doppler ambiguities to allow for HRWS SAR imaging and null the clutter directions to suppress clutter for ground moving target indication. In this paper, the degrees of freedom in azimuth for the multichannel SAR systems are employed to implement clutter suppression. First, the clutter and moving target echoes are transformed into the range compression and azimuth chirp Fourier transform frequency domain, i.e., coarse-focused images formation, when the clutter echoes are with azimuth Doppler ambiguity. Considering that moving targets are sparse in the imaging scene and that there is a difference between clutter and a moving target in the spatial domain, a series of spatial domain filters are constructed to extract moving target echoes. Then, using an extracted moving target echo, two groups of signals are formed, and slant-range velocity of a moving target can be estimated based on baseband Doppler centroid estimation algorithm and multilook cross-correlation Doppler centroid ambiguity number resolving approach. After the linear range cell migration correction and azimuth focus processing, a well-focused moving target image can be obtained. In addition, the proposed clutter suppression and imaging approach is not only adapted for uniformly displaced phase center sampling but also for the nonuniform sampling cases. Some simulation experiments are taken to demonstrate our proposed algorithms. Finally, some real measured data results are presented to validate the theoretical investigations and the proposed approaches. [ABSTRACT FROM PUBLISHER]

A time or spatial series of drop counts is but one realization of a multiple stochastic process. In this paper, a method is presented that extracts more of the information contained in the time series of 1-min Joss-Waldvogel disdrometer counts in rain than a simple analysis of the magnitudes of the counts would provide. This is done by greatly increasing the size of a data set using a Bayesian analysis of drop count measurements in 17 size bins. Using the empirical copula statistical technique of probability density function transformations, a 1391-min time series of drop counts was expanded to the equivalent of 40 000 min. This dramatic increase in sample size permits a deeper characterization of the rain. Using this single disdrometer, it also allows one to translate these counts into a 200 × 200 grid filled at each point with drop size distributions of mean drop concentrations consistent with the observed statistical properties of the rain. Such a field can be used for remote sensing studies of the effect of partial beam filling and for algorithm development. Moreover, since there is nothing unique to this set of drop counts, this approach can be applied to any other set of count data, including snow and clouds. [ABSTRACT FROM PUBLISHER]

The aim of this paper was to test the capabilities of the Sentinel-1 radar data in downscaling Soil Moisture Active Passive (SMAP) radiometer data for high-resolution soil moisture estimation. Three different active–passive downscaling algorithms, including the brightness temperature-based downscaling algorithm (BTBDA), the soil moisture-based downscaling algorithm (SMBDA), and a change detection method (CDM), were analyzed using pairs of Sentinel-1 active and SMAP passive observations collected over a semiarid landscape in southeastern Australia from May 2015 to May 2016. While these algorithms have been tested previously, this is the first study to evaluate the three algorithms using real Sentinel-1 radar and SMAP radiometer data. The SMAP passive observations were disaggregated to 9-, 3-, and 1-km scales and then compared with ground soil moisture measurements. The results suggest that the root-mean-square error (RMSE) in downscaled soil moisture at 9-km resolution was 0.057, 0.056, and 0.067 cm3/cm3 for the BTBDA, SMBDA, and CDM, respectively. The accuracy of downscaling methods was generally decreased when applied at the finer spatial resolution. The SMBDA had overall better performance in terms of correctly detecting the soil moisture pattern and relatively lower RMSE values, and is, therefore, recommended for the combined Sentinel-1 radar and SMAP radiometer setup for soil moisture monitoring. The influence of incidence angle normalization of Sentinel-1 SAR data on downscaled soil moisture was also investigated and found to be minimal. [ABSTRACT FROM AUTHOR]

In this paper, a simple algorithm is proposed to perform the nowcasting of rainfall in the tropical region. The algorithm applies global positioning system-derived precipitable water vapor (PWV) values and its second derivative for the short-term prediction of rainfall. The proposed algorithm incorporates the seasonal dependency of PWV values for the prediction of a rain event in the coming 5 min based on the past 30 min of PWV data. This proposed algorithm is based on the statistical study of four-year PWV and rainfall data from a station in Singapore and is validated using two-year independent data for the same station. The results show that the algorithm can achieve an average true detection rate and a false alarm rate of 87.7% and 38.6%, respectively. To analyze the applicability of the proposed algorithm, further validations are done using one-year data from one independent station from Singapore and two-year data from one station from Brazil. It is shown that the proposed algorithm performs well for both the independent stations. For the station from Brazil, the average true detection and false alarm rates are around 84.7% and 37%, respectively. All these observations suggest that the proposed algorithm is reliable and works well with a good detection rate. [ABSTRACT FROM AUTHOR]

Accurate information about ocean surface currents is of high importance for many purposes. Traditional ways of measuring surface currents using instruments immersed in water are associated with various challenges, as well as high installation and maintenance costs. Algorithms and hardware used to measure surface currents based on marine X-band radar images have been steadily improved during the recent decades. Such systems can now provide the user with current measurements with high accuracy in local areas of interest and with no need for underwater equipment, both from fixed sites and from moving installations, such as vessels. The primary focus of the work presented in this paper is to examine the measurement accuracy of an X-band radar-based system. For this purpose, field trials have been arranged and data acquired from a fixed platform in the North Sea in November 2015 and December 2015 and from a research vessel in the Norwegian Sea and the Barents Sea in November 2016. Data are analyzed, measured time series are compared to reference measurements and tidal current models, and statistical estimates of the measurement accuracy are calculated. Root-mean-square measurement errors of 0.032 m/s and 9.1° for magnitude and direction, respectively, are estimated for the fixed installation. Correlation coefficients are 0.93 and 0.94 for east–west and north–south components, respectively. A pooled standard deviation of 0.059 m/s is estimated for the moving installation. [ABSTRACT FROM AUTHOR]

A novel algorithm is put forward to separate radar target and moving clutter based on a combination of the low-rank matrix optimization (LRMO) and the decision tree. Similar to the moving object detection in the field of automated video analysis, the proposed separation method, which is carried out in the range-Doppler domain, makes use of different motion variations of radar target and clutter in the spectrogram sequence. The technique is very general, but the focus of this paper is on narrowband moving clutter suppression in a weather radar. The first step in implementing this method is the generation of a range-Doppler spectrogram sequence. Then, the LRMO is applied to the obtained sequence to divide target and moving clutter into foreground and background. From the foreground sequence which is obtained by solving the LRMO, foreground frequency and spectral width are combined in a decision tree to obtain a filtering mask to mitigate the narrowband moving clutter and noise. Data collected by a polarimetric Doppler weather radar known as the IRCTR Drizzle Radar are used to validate the performance of the proposed algorithm. Moreover, its effectiveness in removing narrowband moving clutter is quantitatively assessed through comparisons with another clutter mitigation method. [ABSTRACT FROM AUTHOR]

This paper concerns two polarization orientation parameters of built-up areas derived from the polarimetric synthetic aperture radar (PolSAR) data considering two modeling orthogonal dihedral structures. The first orientation parameter is derived from the well-known circular polarization algorithm with the enrichment of arc distance median filtering using an adaptive neighborhood. The derivation of the second orientation parameter is realized by combining the slope-induced changes in polarimetric orientation angle with the shape-from-shading technique. The combination provides the possibility to measure the incidence angle and the azimuth component of the terrain slopes from the cross-pol SAR intensity image. With reference to the cross scattering model, a doubled cross scattering model (DCSM) is introduced by incorporating the orientation parameters, thus serving to guide the model-based decomposition. Using the DCSM refines the estimation of the cross-pol component by enabling us to further reveal the scattering characteristics of built-up areas. Following a novel criterion, the decomposition is implemented at two layers: one for urban areas and one for nonurban areas. The performance of parameter derivation is demonstrated and evaluated with airborne synthetic aperture radar, uninhabited aerial vehicle synthetic aperture radar, and GF-3 fully PolSAR data over different test sites. The decomposed results are consistent with the reference information provided by the National Land Cover Database 2011 about the land cover classification of test sites and encourage the use of the proposed decomposition scheme for different applications. [ABSTRACT FROM AUTHOR]

Anomaly detection has been known to be a challenging problem due to the uncertainty of anomaly and the interference of noise. In this paper, we focus on anomaly detection in hyperspectral images (HSIs) and propose a novel detection algorithm based on spectral unmixing and dictionary-based low-rank decomposition. The innovation is threefold. First, due to the highly mixed nature of pixels in HSI data, instead of using the raw pixel directly for anomaly detection, the proposed algorithm applies spectral unmixing to obtain the abundance vectors and uses these vectors for anomaly detection. We show that the abundance vectors possess more distinctive features to identify anomaly from background. Second, to better represent the highly correlated background and the sparse anomaly, we construct a dictionary based on the mean shift clustering of the abundance vectors to improve both the discriminative and representative powers of the algorithm. Finally, a low-rank matrix decomposition method based on the constructed dictionary is proposed to encourage the coefficients of the dictionary, instead of the background itself, to be low rank, and the residual matrix to be sparse. Anomalies can then be extracted by summing up the columns of the residual matrix. The proposed algorithm is evaluated on both synthetic and real data sets. Experimental results show that the proposed approach constantly achieves high detection rate, while maintaining low false alarm rate regardless of the type of images tested. [ABSTRACT FROM AUTHOR]

This paper focuses on the endmember extraction (EE) technique for analyzing hyperspectral images. We first prove that the reconstruction errors (REs) and abundance anomalies (AAs) (abundances that fail to satisfy the abundance constraints) are effective in extracting undetected endmembers. Then, according to the spatial continuity of the endmember objects and differing from noise or outliers with a sparse distribution, the endmembers are assumed to be located at some salient areas in the RE and AA maps. A novel EE algorithm termed saliency-based endmember detection (SED) is proposed, where the visual saliency model is introduced to explore and analyze the spatial information that is contained in the AA and RE maps. Specifically, the AA and RE maps are regarded as the visual inputs, whereas the endmembers are treated as the visual stimuli. In SED, we assume that the pure pixel assumption holds. Based on the characteristics of the human visual system, the proposed method can not only extract endmembers in homogenous areas, but it can also highlight the small targets whose abundances may be spatially varied. In addition, since the spatial information is exploited in the reconstruction, the capability of the endmembers to represent the hyperspectral scene is automatically considered in the process of EE, and the detected endmembers are both accurate and reliable. The experimental results obtained on both simulated and real hyperspectral data confirm the merits and viability of the proposed algorithm. [ABSTRACT FROM AUTHOR]

Many hyperspectral anomaly detectors are designed based on the traditional Mahalanobis distance-based RX algorithm, which is usually considered as an inverse operation of the principal components analysis. Such detectors include the uniform target detector (UTD) algorithm, RX-UTD algorithm, and so on. However, the possibility of background statistical contamination caused by anomalies still exists. In order to alleviate this problem, in this paper, we propose a spectral matched filter (background-anomaly component projection and separation optimized filter) to minimize the average output energy of separate image components and the output values of the weighted background regular term for hyperspectral image anomaly detection, which could strengthen the separation between anomalies and backgrounds. By calculating the optimal solution to the background-anomaly component projection and separation function, we obtain the optimal projection, where we can effectively suppress the background while highlighting the anomalies. Proposed algorithm has the following research advantages: 1) it creates a novel collaborative component projection and robust background optimization function to separate the background and anomalies and 2) it analyzes the intrinsic statistical distribution of pixels and applies appropriate iterative shrinkage-thresholding algorithm to solve the $\ell _{1}$ -min problem. Experiments were conducted on three real hyperspectral data sets. The detection results demonstrate that the proposed algorithm is superior to other state-of-the-art anomaly detection algorithms. [ABSTRACT FROM AUTHOR]

The Soil Moisture Active Passive (SMAP) satellite mission launched in January 2015 provides worldwide soil moisture (SM) monitoring based on L-band brightness temperature ($T_{B}^{\mathrm{ p}}$) measurements at horizontal ($\text{T}_{B}^{\mathrm{ H}}$) and vertical ($T_{B}^{\mathrm{ V}}$) polarizations. This paper presents a performance assessment of SMAP soil emission model and SM retrieval algorithms for a Tibetan desert ecosystem. It is found that the SMAP emission model largely underestimates the SMAP measured $T_{B}^{\mathrm{ H}}$ (≈ 15 K), and the $T_{B}^{\mathrm{ V}}$ is underestimated during dry-down episodes. A cold bias is noted for the SMAP effective temperature due to underestimation of soil temperature, leading to the $T_{B}^{\mathrm{ p}}$ underestimation (>5 K). The remaining $T_{B}^{\mathrm{ H}}$ underestimation is found to be related to the surface roughness parameterization that underestimates its effect on modulating the $T_{B}^{\mathrm{ p}}$ measurements. Further, the topography and uncertainty of soil information are found to have minor impacts on the $T_{B}^{\mathrm{ p}}$ simulations. The SMAP baseline SM products produced by single-channel algorithm (SCA) using the $T_{B}^{\mathrm{ V}}$ measurements capture the measured SM dynamics well, while an underestimation is noted for the dry-down periods because of $T_{B}^{\mathrm{ V}}$ underestimation. The products based on the SCA with $T_{B}^{\mathrm{ H}}$ measurements underestimate the SM due to underestimation of $T_{B}^{\mathrm{ H}}$ , and the dual-channel algorithm overestimates the SM. After implementing a new surface roughness parameterization and improving the soil temperature and texture information, the deficiencies noted above in $T_{B}^{\mathrm{ p}}$ simulation and SM retrieval are greatly resolved. This indicates that the SMAP SM retrievals can be enhanced by improving both surface roughness and adopted soil temperature and texture information for Tibetan desert ecosystem. [ABSTRACT FROM AUTHOR]

The spatial modulation coincidence imaging (SMCI), as a novel kind of microwave coincidence imaging method, is proposed in this paper. The SMCI system provides a new way to produce the time-space independent signal instead of multitransmitter architecture with wideband randomly modulated signal in radar coincidence imaging. Due to some special features, metamaterial plate is utilized as the reflective surface to modulate the incident signal to construct random radiation field. The resolution of SMCI system is analyzed under large viewing angle with two different transmitting signals. Reflective surface is nonuniformly divided to derive the expression of resolution. The analysis results show that the resolution of SMCI system is mainly determined by the size of reflective surface and center frequency, which is similar to the traditional aperture. The SMCI system is low cost and flexible in design. Simultaneously, it can avoid the synchronization problem between subsources. Moreover, the SMCI system can achieve the resolution of space target through single transmitter–single receiver radar system. High-resolution image can be reconstructed since the tests are nonlinear. Finally, a series of simulation experiments is presented based on the nondirect-viewing scene we proposed. Using the algorithm based on a compressed sensing theory, we reconstructed the target image with high resolution. [ABSTRACT FROM AUTHOR]

This paper describes the adaptation of a Bayesian sea ice detection algorithm for the scatterometer on-board the European Remote Sensing (ERS) satellites (ERS-1 and ERS-2). The algorithm is based on statistics of distances to ocean wind and sea ice geophysical model functions (GMFs) and its performance is validated against coincident active and passive microwave data. We furthermore propose a new model for sea ice backscatter at the C-band in vertical polarization based on the sea ice GMFs derived from ERS and advanced scatterometer data. The model characterizes the dependence of sea ice backscatter on the incidence angle and the sea ice type, allowing a more precise incidence angle correction than afforded by the usual linear transformation. The resulting agreement between the ERS, QuikSCAT, and special sensor microwave imager sea ice extents during the year 2000 is high during the fall and winter seasons, with an estimated ice edge accuracy of about 20 km, but shows persistent biases between scatterometer and radiometer extents during the melting period, with scatterometers being more sensitive to summer (lower concentration and rotten) sea ice types. [ABSTRACT FROM AUTHOR]