Based on the spatiotemporally weighted two-step method (STW-TSM), the spatiotemporal characteristics of the residual microwave brightness temperature (MBT) with the Mw7.9 Wenchuan earthquake on May 12, 2008 are revealed by satellite data from the Advanced Microwave Scanning Radiometer - Earth Observing System (AMSR-E) sensor. Two significant MBT positive anomalies are found to be exactly in spatial accordance with two geological Quaternary zones, and the detailed geometric information of the MBT positive anomaly is found to be correlated with the microwave frequencies. After eliminating other possible influential factors, including surface temperature, vegetation index, land-surface roughness, and surface soil moisture under the conditions of space, time, and magnitude, and according to the microwave radiative transfer model, the dielectric variation in the ground surface is suggested to be the primary contributor of the MBT positive anomaly. The positive-hole (P-hole) theory is applied to interpret the geological preference of the MBT positive anomalies through a chain process: crustal stress enhancing—P-hole producing and flowing down stress gradients—surface P-hole accumulating—dielectric constant decreasing—and microwave radiation increasing. The stress-resulting effect of the dielectric decrease on MBT the increase provides a novel mechanism for microwave remote-sensing monitoring of crustal stress field alteration, earthquake preparation, and upcoming shocks. This research has a particular significance for searching potential georelations between the tectonic earthquake preparation and the abnormal satellite MBT.
Building change detection (CD), important for its application in urban monitoring, can be performed in near real time by comparing prechange and postchange very-high-spatial-resolution (VHR) synthetic-aperture-radar (SAR) images. However, multitemporal VHR SAR images are complex as they show high spatial correlation, prone to shadows, and show an inhomogeneous signature. Spatial context needs to be taken into account to effectively detect a change in such images. Recently, convolutional-neural-network (CNN)-based transfer learning techniques have shown strong performance for CD in VHR multispectral images. However, its direct use for SAR CD is impeded by the absence of labeled SAR data and, thus, pretrained networks. To overcome this, we exploit the availability of paired unlabeled SAR and optical images to train for the suboptimal task of transcoding SAR images into optical images using a cycle-consistent generative adversarial network (CycleGAN). The CycleGAN consists of two generator networks: one for transcoding SAR images into the optical image domain and the other for projecting optical images into the SAR image domain. After unsupervised training, the generator transcoding SAR images into optical ones is used as a bitemporal deep feature extractor to extract optical-like features from bitemporal SAR images. Thus, deep change vector analysis (DCVA) and fuzzy rules can be applied to identify changed buildings (new/destroyed). We validate our method on two data sets made up of pairs of bitemporal VHR SAR images on the city of L’Aquila (Italy) and Trento (Italy).
Supervised classification algorithms require a sufficiently large set of representative training samples to generate accurate land-cover maps. Collecting reference data is difficult, expensive, and unfeasible at the large scale. To solve this problem, this article introduces a novel approach that aims to extract reliable labeled data from existing thematic products. Although these products represent a potentially useful information source, their use is not straightforward. They are not completely reliable since they may present classification errors. They are typically aggregated at polygon level, where polygons do not necessarily correspond to homogeneous areas. Finally, usually, there is a semantic gap between map legends and remote sensing (RS) data. In this context, we propose an approach that aims to: 1) perform a domain understanding to detect the discrepancies between the thematic map domain and the RS data domain; 2) use RS data contemporary to the map to decompose the thematic product from the semantic and spatial viewpoints; and 3) extract a database of informative and reliable training samples. The database of weak labeled units is used for training an ensemble of classifiers on recent data whose results are then combined in a majority voting rule. Two sets of experimental results obtained on MS images by extracting training samples from a crop type map and the 2018 Corine Land Cover (CLC) map, respectively, confirm the effectiveness of the proposed approach.
Full waveform inversion (FWI) of ground-penetrating radar (GPR) data is becoming a promising technique to facilitate the interpretation of surface-GPR data and the mapping of the subsurface. However, more general FWIs still require a sufficient amount of RAM memory, and it is difficult to produce an accurate and representative reconstruction result due to a large amount of the Hessian matrix calculations and singular value decomposition (SVD). In this article, we developed a novel full-waveform approach of the limited-memory Broyden–Fletcher–Goldfarb–Shanno (L-BFGS) algorithm for surface-GPR data that is based on a quasi-Newton framework and the total variation (TV) regularization. The proposed approach uses an L-BFGS algorithm and combines a scale-transformation regularization technique to mitigate the ill-posed problem of inversion, which can impose a biparameter preinformation constraint to ensure the stability of inversion, and adaptive regularization weights are applied to improve the convergence efficiency of inversion. To demonstrate the novelty and effectiveness of the proposed scheme, we tested our FWI algorithm using synthetic data and in-site field data. In the testing, we focus on analyzing the influence of different aspects of the FWI results, including different scale factors, regularization weights, inversion strategies, acquisition configurations, initial models, and the noisy data set. In particular, the FWI experiment is performed to demonstrate the applicability of the proposed algorithm. The results show that the proposed algorithm can effectively reconstruct the biparameter near the subsurface with high accuracy, which makes our approach very attractive for attribute analysis applications and makes the surface-GPR FWI commercially viable.
FILtering by COHerence (FILCOH) is a novel technique for mitigating ground clutter echoes of weather radar, particularly those transmitting in multiple pulse repetition time (PRT) schemes. FILCOH takes advantage of the differences in correlation time for separating rain and ground echoes. In short, the longest time-lag coefficients of the autocorrelation function are used to estimate the ground clutter contribution which is then subtracted from the shortest lag coefficients before the extraction of the meteorological parameters using the conventional estimators. The initial results are obtained for uniform PRT data provided by the Degreane Horizon wind profiler radar. The analysis illustrates the behavior of FILCOH filtering and its high performance in such cases. The simulation results are then presented for 2- and 3-PRT pulse schemes using standard PRT ratios. For the 2-PRT scheme, FILCOH filters out the ground clutter echoes up to a clutter to signal ratio (CSR) of 65 dB for the received power and up to 45 and 60 dB for the radial velocity and spectrum width, respectively. For the 3-PRT scheme, ground clutter filtering is effective up to $CSR = 45$ dB for the received power and up to $CSR = 35$ dB for radial velocity and spectrum width. These results are confirmed by 3-PRT real data issued from the French X- and C-band radars. The observed performances are, however, a little less than those of the simulation since we observed a maximum 40-dB attenuation of the ground clutter echoes of the studied rain events.
Human sensing using WiFi signal transmissions is attracting significant attention for future applications in e-healthcare, security, and the Internet of Things (IoT). The majority of WiFi sensing systems are based around processing of channel state information (CSI) data which originates from commodity WiFi access points (APs) that have been primed to transmit high data-rate signals with high repetition frequencies. However, in reality, WiFi APs do not transmit in such a continuous uninterrupted fashion, especially when there are no users on the communication network. To this end, we have developed a passive WiFi radar system for human sensing which exploits WiFi signals irrespective of whether the WiFi AP is transmitting continuous high data-rate Orthogonal Frequency-Division Multiplexing (OFDM) signals, or periodic WiFi beacon signals while in an idle status (no users on the WiFi network). In a data transmission phase, we employ the standard cross ambiguity function (CAF) processing to extract Doppler information relating to the target, while a modified version is used for lower data-rate signals. In addition, we investigate the utility of an external device that has been developed to stimulate idle WiFi APs to transmit usable signals without requiring any type of user authentication on the WiFi network. In this article, we present experimental data which verifies our proposed methods for using any type of signal transmission from a standalone WiFi device, and demonstrate the capability for human activity sensing.
In this article, we propose a phase-derived velocity measurement (PDVM) method for high-speed spatial targets based on the stepped-frequency chirp signal (SFCS). This method is capable of accurately measuring the velocity of high-speed targets and yields root-mean-squared error values at the level of centimeters per second; therefore, it has great potential for measuring the micromotion of targets and is of significant importance for target recognition. The traditional phase-derived measurement method is not applicable for high-speed targets. The main challenge that we have solved is how to extract the echo phase from the high-resolution range profile, which is corrupted by range migration, intrapulse motion, and range straddling under high-speed target conditions. To guide the implementation of the proposed method in radar systems, constraint conditions for the compensation accuracy are thoroughly derived and systematically justified under different radar parameter settings. The simulation results are presented to validate the high accuracy of the method under various circumstances. In addition, the small-amplitude micromotion measurement capability of the proposed method is verified, and reconstruction of the target micromotion trajectory is demonstrated.
Many modern phased-array radars (PARs) are multichannel systems that include multiple receivers for data acquisition. Each channel provides a signal from a group of Transmit/Receive modules comprising a section of the antenna. Channels typically consist of a full receive path, often with an independent local oscillator (LO) clock source. Such arrangement provides for beamforming flexibility on receive which can be applied in a digital domain. Consequently, the channel-to-channel phase and magnitude alignment is critical to maximizing the performance of the digital beamforming process and the accuracy of resulting detections and measurements. Herein, a novel method to improve such alignment using weather returns and achieve the improvement in the polarimetric variable estimates is described.
With the advances in radar technology, through-the-wall radar imaging (TWRI) has become a viable sensing modality that can allow fire-and-rescue personnel, police, and military forces to detect, localize, and identify targets behind opaque obstacles. Many of the existing TWRI approaches detect either stationary or moving targets but not both of them simultaneously. In this article, a method is proposed to detect both stationary and moving targets from a sequence of radar signals. The proposed method decomposes the 3-D radar data, i.e., frequency, space, and time data into a low-rank tensor and two sets of sparse images. One set of images comprises the stationary targets, and the other set of images contains the moving targets. Wall clutter removal and target detection are formulated into an optimization problem regularized by tensor low-rank, joint sparsity, and total variation constraints. Then, an alternating direction technique is developed to reconstruct the sets of stationary and moving target images. Experiments using simulated and real radar signals are conducted. The experimental results illustrate the effectiveness of the proposed method to detect and separate the stationary and moving targets into a pair of sparse images.
Model-based incoherent polarimetric decomposition is a frequently used technique to analyze multilook data of polarimetric synthetic aperture radars (POLSARs). The purpose of this study is to analyze and compare different model-based incoherent polarimetric decomposition algorithms from the polarimetric information change aspect. For the input of a model-based incoherent polarimetric decomposition algorithm, polarimetric entropy was used to represent the polarimetric information of a coherency matrix. For the output of a model-based incoherent polarimetric decomposition algorithm, there are usually several decomposed components. To quantitatively represent their total polarimetric information, a new concept, generalized polarimetric entropy, was proposed which generalized the concept of polarimetric entropy based on the information entropy additivity of information theory. Generalized polarimetric entropy consists of two parts named as polarimetric power entropy and polarimetric residual entropy, respectively. Polarimetric power entropy describes the distribution status of the Span values of all decomposed components. Polarimetric residual entropy represents the residual randomness of all decomposed components. With the three new concepts, eight model-based incoherent polarimetric decomposition algorithms were compared and analyzed. Two real POLSAR images, respectively, derived by the E-SAR airborne system of Germany and the GF-3 satellite of China were used for the experiments. Experimental results had illustrated several useful conclusions.
In high-frequency over-the-horizon radar (OTHR), it is a challenging work to detect targets in the nonhomogeneous range-Doppler (RD) map with multitarget interference and sharp/smooth clutter edges. The intensity transition of the clutter edge may be sharp or smooth due to the coexistence of atmospheric noise, sea clutter, and ionospheric clutter in OTHR. The analysis of the RD map shows the spatial correlation among neighboring cell-under-test (CUT) that varies from clutter to clutter. This article proposes an algorithm that uses the spatial relationship to estimate the statistical distribution parameters of every CUT by the adaptive tight frame and the weighted group-sparsity regularization. In the proposed algorithm, the spatial relationship is formulated mathematically by regularization terms and combined with the log-likelihood function of CUTs to construct the objective function. The proposed algorithm is verified by the simulated data and real RD maps collected from both trial sky-wave and surface-wave OTHRs in which it shows robust and improved detection.
Bistatic forward-looking spotlight synthetic aperture radar with stationary transmitter and maneuvering receiver (STMR-BFSSAR) is a promising sensor for various applications, such as the automatic navigation and landing of maneuvering vehicles. Because of the bistatic forward-looking configuration and the receiver’s maneuvers, conventional image formation algorithms suffer from high computational complexity or small size of a well-focused scene if applied to STMR-BFSSAR. In this article, we propose a wavenumber-domain algorithm for STMR-BFSSAR image formation, which is termed the bistatic-range-Doppler-aperture wavenumber algorithm (BDWA). First, a novel range model in bistatic-range and Doppler-aperture coordinate space instead of conventional Cartesian coordinate space is established by employing the elliptic polar coordinate system and the method of series reversion. The novel range model not only makes the echo’s samples to be regular along the direction of the bistatic-range wavenumber axis but also constructs a curved wavefront close to the true wavefront. Second, an operation termed wavenumber-domain gridding is conceived to regularize the echo’s samples along the Doppler-aperture wavenumber axis, which can be implemented by 1-D interpolation. The proposed algorithm significantly outperforms the conventional algorithms in terms of computational complexity and scene size limits. Both point and distributed targets are simulated for two STMR-BFSSAR systems with different parameters. The simulation results verify the validity and superiority of the proposed BDWA.
The microwave photonic (MWP) radar technique is capable of providing ultrawide frequency bandwidth waveforms to generate ultrahigh-resolution (UHR) inverse synthetic aperture radar (ISAR) imagery. Nevertheless, conventional ISAR imaging algorithms have limitations in focusing UHR MWP-ISAR imagery, where high-precision high-order range cell migration (RCM) and phase correction are crucially necessary. In this article, a UHR MWP-ISAR imaging algorithm integrating rotation estimation and high-order motion terms compensation is proposed. By establishing the relationship between parametric ISAR rotation model and high-order motion terms, an average range profile sharpness maximization (ARPSM) is developed to obtain rotation velocity by using nonuniform fast Fourier transform (NUFFT). Second-order range-dependent RCM is corrected with parametric compensation model by using the rotation velocity estimation. Furthermore, the spatial-variant high-order phase error is extracted to compensation by the entire image sharpness maximization (EISM). A new imaging framework is established with two one-dimensional (1-D) parameter estimations: ARPSM and EISM. Extensive experiments demonstrate that the proposed algorithm outperforms traditional ISAR imaging strategies in high-order RCM correction and azimuth focusing performance.
Sentinel-1 (S1) extra-wide (EW) swath data in cross-polarization (horizontal–vertical, HV or vertical–horizontal, VH) are strongly affected by the scalloping effect and thermal noise, particularly over areas with weak backscattered signals, such as sea surfaces. Although noise vectors in both the azimuth and range directions are provided in the standard S1 EW data for subtraction, the residual thermal noise still significantly affects sea ice detection by the EW data. In this article, we improve the denoising method developed in previous studies to remove the additive noise for the S1 EW data in cross-polarization. Furthermore, we propose a new method for eliminating the residual noise (i.e., multiplicative noise) at the subswath boundaries of the EW data, which cannot be well processed by simply subtracting the reconstructed 2-D noise field. The proposed method of removing both the additive and multiplicative noise was applied to EW HV-polarized images processed using different Instrument Processing Facility (IPF) versions. The results suggest that the proposed algorithm significantly improves the quality of EW HV-polarized images under various sea ice conditions and sea states in the marginal ice zone (MIZ) of the Arctic. This is of great support for the utilization of cross-polarization synthetic aperture radar (SAR) images in wide swaths for intensive sea ice monitoring in polar regions.
This article addresses the general problem of single-look multi-master SAR tomography. For this purpose, we establish the single-look multi-master data model, analyze its implications for the single and double scatterers, and propose a generic inversion framework. The core of this framework is the nonconvex sparse recovery, for which we develop two algorithms: one extends the conventional nonlinear least squares (NLS) to the single-look multi-master data model and the other is based on bi-convex relaxation and alternating minimization (BiCRAM). We provide two theorems for the objective function of the NLS subproblem, which lead to its analytic solution up to a constant phase angle in the 1-D case. We also report our findings from the experiments on different acceleration techniques for BiCRAM. The proposed algorithms are applied to a real TerraSAR-X data set and validated with the height ground truth made available by an SAR imaging geodesy and simulation framework. This shows empirically that the single-master approach, if applied to a single-look multi-master stack, can be insufficient for layover separation, and the multi-master approach can indeed perform slightly better (despite being computationally more expensive) even in the case of single scatterers. In addition, this article also sheds light on the special case of single-look bistatic SAR tomography, which is relevant for the current and future SAR missions such as TanDEM-X and Tandem-L.
The edges of a target provide essential geometric information and are extremely important for human visual perception and image recognition. However, due to the coherent superposition of received echoes, the continuous edges of targets are discretized in synthetic aperture radar (SAR) images, i.e., the edges become dispersed points, which seriously affects the extraction of visual and geometric information from SAR images. In this article, we focus on solving the problem of how to recover smooth linear edges (SLEs). By introducing multiangle observations, we propose an SAR parametric image reconstruction method (SPIRM) that establishes a parametric framework to recover SLEs from SAR echoes. At the core of the SPIRM is a novel physical characteristic parameter called the scattering-phase-mutation feature (SPMF), which reveals the most essential difference between the residual endpoints of a disappeared SLE and points. Numerical simulations and real-data experiments demonstrate the robustness and effectiveness of the proposed method.
The active recognition of interesting targets has been a vital issue for synthetic aperture radar (SAR) systems. The SAR recognition methods are mainly grouped as follows: extracting image features from the target amplitude image or matching the testing samples with the template ones according to the scattering centers extracted from the target complex data. For amplitude image-based methods, convolutional neural networks (CNNs) achieve nearly the highest accuracy for images acquired under standard operating conditions (SOCs), while scattering center feature-based methods achieve steady performance for images acquired under extended operating conditions (EOCs). To achieve target recognition with good performance under both SOCs and EOCs, a feature fusion framework (FEC) based on scattering center features and deep CNN features is proposed for the first time. For the scattering center features, we first extract the attributed scattering centers (ASCs) from the input SAR complex data, then we construct a bag of visual words from these scattering centers, and finally, we transform the extracted parameter sets into feature vectors with the $k$ -means. For the CNN, we propose a modified VGGNet, which can not only extract powerful features from amplitude images but also achieve state-of-the-art recognition accuracy. For the feature fusion, discrimination correlation analysis (DCA) is introduced to the FEC framework, which not only maximizes the correlation between the CNN and ASCs but also decorrelates the features belonging to different categories within each feature set. Experiments on Moving and Stationary Target Acquisition and Recognition (MSTAR) database demonstrate that the proposed FEC achieves superior effectiveness and robustness under both SOCs and EOCs.
Supervised change detection methods always face a big challenge that the current scene (target domain) is fully unlabeled. In remote sensing, it is common that we have sufficient labels in another scene (source domain) with a different but related data distribution. In this article, we try to detect changes in the target domain with the help of the prior knowledge learned from multiple source domains. To achieve this goal, we propose a change detection framework based on selective adversarial adaptation. The adaptation between multisource and target domains is fulfilled by two domain discriminators. First, the first domain discriminator regards each scene as an individual domain and is designed for identifying the domain to which each input sample belongs. According to the output of the first domain discriminator, a subset of important samples is selected from multisource domains to train a deep neural network (DNN)-based change detection model. As a result, not only the positive transfer is enhanced but also the negative transfer is alleviated. Second, as for the second domain discriminator, all the selected samples are thought from one domain. Adversarial learning is introduced to align the distributions of the selected source samples and the target ones. Consequently, it further adapts the knowledge of change from the source domain to the target one. At the fine-tuning stage, target samples with reliable labels and the selected source ones are used to jointly fine-tune the change detection model. As the target domain is fully unlabeled, homogeneity- and boundary-based strategies are exploited to make the pseudolabels from a preclassification map reliable. The proposed method is evaluated on three SAR and two optical data sets, and the experimental results have demonstrated its effectiveness and superiority.
Interferometric synthetic aperture radar (InSAR) allows for mapping of crustal deformation on land with high spatial resolution and precision in areas with high signal-to-noise ratios. Efforts to obtain precise displacement time series globally, however, are severely limited by radar path delays within the troposphere. The tropospheric delay is integrated along the full path length between the ground and the satellite, resulting in correlations between the interferometric phase and elevation that can vary dramatically in both space and time. We evaluate the performance of spatially variable, empirical removal of phase-elevation dependence within SAR interferograms through the use of the $K$ -means clustering algorithm. We apply this method to both synthetic test data, as well as to C-band Sentinel-1a/b time series acquired over a large area in south-central Mexico along the Pacific coast and inland—an area with a large elevation gradient that is of particular interest to researchers studying tectonic- and anthropogenic-related deformation. We show that the clustering algorithm is able to identify cases where tropospheric properties vary across topographic divides, reducing total root mean square (rms) by an average of 50%, as opposed to a spatially constant phase-elevation correction, which has insignificant error reduction. Our approach also reduces tropospheric noise while preserving test signals in synthetic examples. Finally, we show the average standard deviation of the residuals from the best-fit linear rate decreases from approximately 3 to 1.5 cm, which corresponds to a change in the error on the best-fit linear rate from 0.94 to 0.63 cm/yr.
In this study, a method to improve the reconstruction performance of antenna-pattern deconvolution based on the gradient iterative regularization scheme is proposed. The method exploits microwave measurements acquired by a multichannel radiometer to enhance their native spatial resolution. The proposed rationale consists of using the information carried on a high-frequency (finer spatial resolution) channel to ameliorate the spatial resolution of the lowest resolution radiometer channel. Experiments performed using both synthetic and real special sensor microwave/imager (SSM/I) radiometer data demonstrate that an enhanced spatial resolution 19.35-GHz channel can be obtained by ingesting in the algorithm information coming from 37.0-GHz channel. This multichannel spatial resolution method is also shown to outperform the conventional gradient-like regularization scheme in terms of both observation of smaller targets and reduction of ringings and fluctuations.
