Multiple Synthetic Aperture Radar Research Articles

ROBUST REGISTRATION FOR OPTICAL AND SAR IMAGES BASED ON SPATIAL GEOMETRIC CONSTRAINT AND STRUCTURE FEATURES

Abstract. Automatic registration of optical and synthetic aperture radar (SAR) images is a challenging task due to significant geometric deformation and radiometric differences between two images. To address this issue, this paper proposes an automatic registration method for optical and SAR images based on spatial geometric constraint and structure features. Firstly, the Harris detector with a block strategy is used to extract evenly distributed feature points in the images. Subsequently, a local geometric correction is performed by using the Rational Function Model, which eliminates the rotation and scale differences between optical and SAR images. Secondly, orientated gradient information of images is used to construct a geometric structural feature descriptor. Then, the feature descriptor is transformed into the frequency domain, and the three-dimensional (3-D) phase correlation is used as the similarity metric to achieve correspondences by employing a template matching scheme. Finally, mismatches are eliminated based on spatial geometric constraint relationship between images, followed by a process of geometric correction to achieve the image registration. Experimental results with multiple high-resolution optical and SAR images show that the proposed method can achieve reliable registration accuracy, and outperforms the state of the art methods.

Signal Processing Options for High Resolution SAR Tomography of Natural Scenarios

Synthetic Aperture Radar (SAR) Tomography is a technique to provide direct three-dimensional (3D) imaging of the illuminated targets by processing SAR data acquired from different trajectories. In a large part of the literature, 3D imaging is achieved by assuming mono-dimensional (1D) approaches derived from SAR Interferometry, where a vector of pixels from multiple SAR images is transformed into a new vector of pixels representing the vertical profile of scene reflectivity at a given range, azimuth location. However, mono-dimensional approaches are only suited for data acquired from very closely-spaced trajectories, resulting in coarse vertical resolution. In the case of continuous media, such as forests, snow, ice sheets and glaciers, achieving fine vertical resolution is only possible in the presence of largely-spaced trajectories, which involves significant complications concerning the formation of 3D images. The situation gets even more complicated in the presence of irregular trajectories with variable headings, for which the one theoretically exact approach consists of going back to raw SAR data to resolve the targets by 3D back-projection, resulting in a computational burden beyond the capabilities of standard computers. The first aim of this paper is to provide an exhaustive discussion of the conditions under which high-quality tomographic processing can be carried out by assuming a 1D, 2D, or 3D approach to image formation. The case of 3D processing is then further analyzed, and a new processing method is proposed to produce high-quality imaging while largely reducing the computational burden, and without having to process the original raw data. Furthermore, the new method is shown to be easily parallelized and implemented using GPU processing. The analysis is supported by results from numerical simulations as well as from real airborne data from the ESA campaign AlpTomoSAR.

A New Concept of Full Ocean Current Vector Retrieval With Spaceborne SAR Based on Intrapulse Beam-Switching Technique

Synthetic aperture radar (SAR) has proven to be an important tool for the retrieval of ocean surface currents. However, existing spaceborne SAR systems can usually measure only one component of an ocean current vector. In this article, we propose a new concept of dual-look spaceborne SAR system for ocean surface motion remote sensing. The key of the proposed SAR system lies in the use of a technique called intrapulse beam switching (IPBS) implemented by a phased array antenna. One example for the application of the proposed IPBS-based SAR system is its combination with the Doppler centroid anomaly (DCA) method. Through the combination, the DCA method, which in its conventional version can measure only the radial component of a current vector, is generalized to measure full ocean current vectors. A system design example and its performance evaluation demonstrate that the accuracies for ocean surface velocity vector measurement better than 5 cm/s in both azimuth and range, and better than 10° in velocity direction for a mild current, and a swath width of about 100 km, can be achieved. A comparison with other beam-generating schemes indicates that a relatively large swath width, without compromising the current-measuring accuracy, can be achieved by the proposed SAR system with one single antenna array and one single physical output channel—a feature that means a high level of integration which in turn allows for the flexibility to operate this system in multiple SAR modes.

The Multiple Aperture SAR Interferometry (MAI) Technique for the Detection of Large Ground Displacement Dynamics: An Overview

This work presents an overview of the multiple aperture synthetic aperture radar interferometric (MAI) technique, which is primarily used to measure the along-track components of the Earth’s surface deformation, by investigating its capabilities and potential applications. Such a method is widely used to monitor the time evolution of ground surface changes in areas with large deformations (e.g., due to glaciers movements or seismic episodes), permitting one to discriminate the three-dimensional (up–down, east–west, north–south) components of the Earth’s surface displacements. The MAI technique relies on the spectral diversity (SD) method, which consists of splitting the azimuth (range) Synthetic Aperture RADAR (SAR) signal spectrum into separate sub-bands to get an estimate of the surface displacement along the azimuth (sensor line-of-sight (LOS)) direction. Moreover, the SD techniques are also used to correct the atmospheric phase screen (APS) artefacts (e.g., the ionospheric and water vapor phase distortion effects) that corrupt surface displacement time-series obtained by currently available multi-temporal InSAR (MT-InSAR) tools. More recently, the SD methods have also been exploited for the fine co-registration of SAR data acquired with the Terrain Observation with Progressive Scans (TOPS) mode. This work is primarily devoted to illustrating the underlying rationale and effectiveness of the MAI and SD techniques as well as their applications. In addition, we present an innovative method to combine complementary information of the ground deformation collected from multi-orbit/multi-track satellite observations. In particular, the presented technique complements the recently developed Minimum Acceleration combination (MinA) method with MAI-driven azimuthal ground deformation measurements to obtain the time-series of the 3-D components of the deformation in areas affected by large deformation episodes. Experimental results encompass several case studies. The validity and relevance of the presented approaches are clearly demonstrated in the context of geospatial analyses.

Azimuth Cutoff Compensation Method for SAR Wave Observation Based on Multiview Wave Spectrum Data Fusion

Azimuth cutoff is an inherent disadvantage of synthetic aperture radar (SAR) waves observation. The waves shorter than some certain wavelength (azimuth cutoff) cannot be imaged in original form (structures) by SAR. For a single SAR observation, the problem of the azimuth cutoff for ocean waves can be resolved to some extent by cooperative observations of SAR satellites. Multiple SAR satellites are required to achieve simultaneous observation of an ocean area to obtain multiview SAR ocean wave synchronization data. Currently, these data cannot be acquired from in-orbit SARs. In this study, imaging simulations of multiview SAR ocean wave synchronization data based on small SAR satellites were carried out for the first time with X-band, 4-m resolution, stripmap mode and Single Look Complex (SLC) product. The Max Planck Institute (MPI) method was used to obtain the optimum wave spectrum of the synchronous data. The influencing factors of the azimuth cutoff wavelength were analyzed by using measured and simulated SAR data. The analysis results were used to develop a novel multiview wave spectrum data fusion method for azimuth cutoff compensation. The azimuth cutoff compensation effect was evaluated by comparing the inversion results before and after data fusion: the azimuth cutoff decreased by 9.76% on average, the root mean square error (RMSE) of the significant wave height is 0.06m, and RMSE of the mean wave period is 0.58s. The azimuth cutoff compensation method can be applied to SAR data for medium and low sea states (that is, wind speeds of 5-15 m/s). These results show that the proposed method of multiview wave spectrum data fusion effectively compensates for azimuth cutoff.

Oil Tank Extraction Based on Joint-Spatial Saliency Analysis for Multiple SAR Images

The lack of true color and the presence of background clutter reduce the accuracy rate of the saliency analysis for oil tank extraction in the synthetic aperture radar (SAR) images. This letter proposes a specially designed unsupervised method to extract oil tanks using the joint-spatial saliency analysis (JSSA) for multiple SAR images. First, the intrasaliency analysis is established on a saliency driven iterative clustering. This considers the spatial intensity and texture feature within a single image and suppresses most backgrounds. Second, the cospatial residual and the local grayscale statistics are considered independently in the intersaliency analysis. The common salient parts among the input series are extracted and used to overcome the problem of the lack of true color. Third, to make the fusion of the two kinds of saliency maps, the low-rank matrix is introduced. The weights of different maps are calculated and the saliency cues are integrated efficiently. Finally, after the statistics of the highlight points within the candidates, the location of the oil tanks is refined. The experiments show the superiority of the proposed method in both the pixel level and the geometric segmentation. The result of the JSSA model appears to improve the accuracy with fewer missing objects compared with the competing algorithms.

Random cross-observation intensity consistency method for large-scale SAR images mosaics: An example of Gaofen-3 SAR images covering China

Large-scale microwave remote sensing research requires multiple synthetic aperture radar (SAR) image mosaics. However, owing to multiple factors, different SAR images frequently exhibit intensity differences. This leads to poor intensity continuity after mosaicking, making SAR image interpretation extremely difficult. The existing method assumes that SAR images have been subject to strict radiometric calibration and mainly focuses on intensity balancing to minimize the intensity difference between adjacent images. This method is not effective when the intensity difference between SAR images is large. In this paper, a novel random cross-observation intensity consistency (RCOIC) algorithm is proposed to eliminate the intensity differences for SAR image mosaics. The proposed algorithm is based on global intensity constraints. First, source images are downsampled. Then, the gain correction coefficients of each image in a study area are accurately estimated through random cross-observation. Further, multi-view local intensity correction of each image in the study area is performed using the overlapping areas of adjacent images, and a low-resolution intensity reference map covering the study area is generated. Finally, the intensities of high-resolution source images are corrected by constructing local independent linear models of images based on the low-resolution intensity reference map. Gaofen-3 satellite images with large ranges and large time differences, covering most of China's land area, were selected to test the effect of the newly developed algorithm. After processing, the overall visual effect was good after mosaicking, the intensity transitions at the edges of the images were smooth, the intensity distributions corresponded better to actual ground objects, and the statistical result of adjacent images overlapping regions was good. This validates the reliability of the algorithm.

Spaceborne tomography of multi-species Indian tropical forests

Synthetic Aperture Radar (SAR) backscatter is the coherent combination of the scattering from multiple individual scatterers within the radar resolution cell, which results into a 2-D radar image. Tomographic SAR (TomoSAR) takes advantage of multiple SAR acquisitions to provide 3-D vertical structure of the imaged target. Over forests, parameters such as canopy structure, canopy density, leaf area and phenology contribute to different SAR scattering mechanisms along the vertical dimension of the trees. Tomography has been demonstrated in various earlier works for forest bio-physical parameter estimation using airborne data. In this research work, multi-polarimetric space-borne TanDEM-X tomograms are examined for the first time over a multi-species Indian tropical forest. Multiple TanDEM-X acquisitions are focused using Capon beamforming at five polarimetric channels — HH, HV, VV, HH − VV and HH + VV. Four distinct forest species compartments are selected representing different canopy structure and density. The tomograms obtained in different polarizations for these species are analyzed in detail to understand the scattering patterns across different forest species. Field surveys carried out in several forest locations provide in situ observations of forest height and vertical structure. It was observed that canopy gaps and leaf density play a crucial role for X-band SAR signal penetration through ground. For dense canopy species the backscatter contributions are spread through the canopy, while in case of sparse canopy species, the ground scattering is dominant. Vertical profiles obtained at surveyed plot locations are plotted in all polarizations, and provide a good agreement with field observations. Further, the obtained TomoSAR backscatter layers have a good correlation with field measured above-ground biomass (AGB). The AGB is modeled from TomoSAR with the HH-pol TomoSAR backscatter layer at 27 m canopy height leading to an AGB estimation with correlation r = 0.76 and RMSE of 50.4 t/ha.

Estimation of rice yield from a C-band radar remote sensing image by integrating a physical scattering model and an optimization algorithm

Accurate estimation of rice yield with remotely sensed data plays a role in ensuring food security at local or national scales. In this study, an estimation scheme integrating a rice canopy scattering model (RCSM) and a genetic algorithm optimization tool (GAOT) was proposed on a basis of radar remote sensing technology. A C-band Radarsat-2 synthetic aperture radar (SAR) image acquired at rice maturity stage in Northeast, China was tested to simulate three yield-related rice parameters via an integrated RCSM-GAOT scheme, with a parallel computing environment. Rice yield was then estimated via a regression analysis by linking the simulated parameters to rice ear weights. Results showed that three parameters, ear length, ear diameter and ear density, can be simulated by the RCSM-GAOT scheme with a simulation error of 1.2 cm, 0.11 cm and 29 ears/m2, respectively. Rice yield is estimated with an average error of 0.28 kg/m2 and 0.22 kg/m2 for fresh ear weight and dry ear weight, respectively. In comparison with the need of multiple SAR data acquisitions in past studies, this study demonstrates the capability of one time Radarsat-2 quad-polarized SAR image in regional rice yield estimation. At the same time, it also indicates that the proposed RCSM-GAOT scheme has a potential in future precision agriculture applications.

Imbalanced Learning-Based Automatic SAR Images Change Detection by Morphologically Supervised PCA-Net

Change detection is a quite challenging task due to the imbalance between unchanged and changed class. In addition, the traditional difference map generated by log-ratio is subject to the speckle, which will reduce the accuracy. In this letter, an imbalanced learning-based change detection is proposed based on PCA network (PCA-Net), where a supervised PCA-Net is designed to obtain the robust features directly from given multitemporal synthetic aperture radar (SAR) images instead of a difference map. Furthermore, to tackle with the imbalance between changed and unchanged classes, we propose a morphologically supervised learning method, where the knowledge in the pixels near the boundary between two classes is exploited to guide network training. Finally, our proposed PCA-Net can be trained by the data sets with available reference maps and applied to a new data set, which is quite practical in change detection projects. Our proposed method is verified on five sets of multiple temporal SAR images. It is demonstrated from the experiment results that with the knowledge in training samples from the boundary, the learned features benefit change detection and make the proposed method outperform than supervised methods trained by randomly drawing samples.

Multi-View Tensor Sparse Representation Model for SAR Target Recognition

The local structure feature of the target in synthetic aperture radar (SAR) image and the inner correlation among multiple SAR images of the same target can effectively improve the recognition performance. In order to utilize these two kinds of information, we propose the multi-view tensor sparse representation (MTSR) model for SAR target recognition. In the proposed model, the SAR image is treated as tensor and represented by the dictionaries on each dimension of the image and corresponding sparse representation coefficient tensor (SRCT). For multi-view SAR images of one target, we design the proposed model by letting the SRCTs of those images have the same structure, which means that the non-zero elements of those SRCTs have the same coordinates. In order to solve the proposed model, we propose the joint tensor orthogonal matching pursuit (JT-OMP) algorithm to calculate the SRCTs of multiple views. JT-OMP ensures that all the SRCTs have the same structure by looking for the atoms of dictionaries that can contribute the total maximum energy for all multi-view tensors in every iteration. To achieve recognition, we construct two dictionaries for each class and compare the total sparse representation error of multi-view SAR images in each class. The experiments conducted on the moving and stationary target acquisition recognition database (MSTAR) verify the performance of the proposed algorithm.

Impact of Dielectric Changes on L-Band 3-D SAR Reflectivity Profiles of Forest Volumes

Synthetic aperture radar (SAR) tomography (TomoSAR) allows the reconstruction of the vertical distribution of the power backscattered by natural volumes by combining multiple SAR images acquired with slightly different incidence angles. Being a “radar” quantity, the profile depends on the radar frequency, polarization, and acquisition geometry, the 3-D distribution of the scattering elements and their dielectric properties. The characterization of each one of these factors is crucial to enable the extraction of physical 3-D structure attributes from TomoSAR profiles. The objective of this paper is to investigate how the vertical distribution of the backscattered power at L-band is affected by seasonal- and weather-induced changes. Radiometric (e.g., ground and volume powers) and geometric (e.g., center of mass of the volume-only profiles and phase centers of the volume scattering layers) parameters have been estimated under different weather and season conditions and compared. Then, TomoSAR data sets affected by dielectric nonstationarity (i.e., variability) have been considered in order to assess the invariance degree of each radiometric and geometric parameters. This analysis has been carried out by processing four L-band airborne TomoSAR data sets acquired before and after a rainfall and in spring and autumn over the Traunstein forest (south of Germany).

Multiple Mode SAR Raw Data Simulation and Parallel Acceleration for Gaofen-3 Mission

Gaofen-3 is China's first meter-level multipolarization synthetic aperture radar (SAR) satellite with 12 imaging modes for the scientific and commercial applications. In order to evaluate the imaging performance of these modes, the multiple mode SAR raw data simulation is highly demanded. In the paper, the multiple mode SAR simulation framework will be briefly introduced to expose how the raw data simulation guarantees the development of Gaofen-3 and its ground processing system. As an engineering simulation, the complex working modes and practical evaluation requirements of Gaofen-3 mission put forward to the higher demand for simulation simplification and data input/output (I/O) efficiency. To meet the requirements, two improvements have been proposed. First, the stripmap mode based multiple mode decomposition method is introduced to make a solid and simplified system simulation structure. Second, the cloud computing and graphics processing unit (GPU) are integrated to simulate the practical huge volume raw data, resulting in improved calculation and data I/O efficiency. The experimental results of sliding spotlight imaging prove the effectiveness of the Gaofen-3 mission simulation framework and the decomposition idea. The results for efficiency assessment show that the GPU cloud method greatly improves the computing power of a 16-core CPU parallel method about $40\times$ speedup and the data throughput with the Hadoop distributed file system. These results prove that the simulation system has the merits of coping with multiple modes and huge volume raw data simulation and can be extended to the future space-borne SAR simulation.

Raw Signal Simulator for SAR With Trajectory Deviation Based on Spatial Spectrum Analysis

This paper presents a whole new method for simulating a raw signal of synthetic aperture radar (SAR) with trajectory deviation. In this paper, a SAR raw signal is viewed as a spatial function in terms of radar location, at a fixed range frequency component of the radar signal. It is disclosed and verified that the SAR raw signals is band limited from the perspective of spatial spectrum analysis. Based on this finding, SAR raw signals under condition of trajectory deviation are obtained by interpolating the SAR raw signal along a set of ideally linear trajectories. The unprecedent performance achieved by the new method mainly includes the following. First, it provides a way to update existing simulators designed for ideally linear trajectory, giving them the capability to take trajectory deviation into consideration. Second, it is an efficient SAR raw signal simulator whose validity limit is greatly expanded compared with existing efficient simulators, as many constraints are eliminated without sacrificing computational efficiency or phase accuracy. Third, it is a uniform framework applicable for multiple SAR acquisition modes. Finally, the computational complexity of Monte Carlo experiment under various trajectory deviations is dramatically reduced, at the expense of storage memory requirement. Theoretical analyses are demonstrated by experimental results.

A Unified Multimode SAR Raw Signal Simulation Method Based on Acquisition Mode Mutation

Raw signal simulators for synthetic aperture radar (SAR) in different acquisition modes have been studied individually. This letter is dedicated to presenting a method for simulating multimode SAR raw signal in a unified framework. To this end, we first simulate stripmap SAR raw signal, and then mutate the raw signal from stripmap mode into the desired acquisition mode. The acquisition mode mutation is implemented by azimuth-time-variant and range-frequency-dependent bandpass filtering. While processing in range-frequency domain brings forth high accuracy, processing in azimuth-time domain makes the method applicable for flexible antenna steering laws specified by multiple SAR acquisition modes, such as staring spotlight mode, sliding spotlight mode, and terrain observation by progressive scans mode. The proposed method is validated by the simulation results.

An Improved Multi-Image Matching Method in Stereo-Radargrammetry

The new generation of synthetic aperture radar (SAR) sensors provides us with an opportunity to match multiple high-resolution SAR images. Moreover, the multiple SAR image matching methods have recently gained a lot of attention due to the fact that they can obtain more accurate, better distributed, and more reliable matches than the stereo matching methods. In this letter, we present an improved multi-image matching method to simultaneously identify matches from multiple SAR amplitude images. The proposed method makes better use of the relationships between the pixels in the deformed correlation window and integrates geometric and radiometric information from multiple SAR images. Experiments on Chinese Academy of Surveying and Mapping Synthetic Aperture Radar (CASMSAR) data sets demonstrate that the improved multi-image matching method is capable of providing more accurate and better distributed matches, as well as offering a better multi-image matching solution in stereo-radargrammetry under the conditions of geometric and radiometric distortions, especially in low-texture areas.

Multi-sensor InSAR Deformation Monitoring over Urban Area of Bratislava (Slovakia)

The integrated use of multiple Synthetic Aperture Radar (SAR) platforms for the deformation monitoring via satellite radar interferometry offers several perspectives for investigation of the behaviour of new and ageing structures, such as buildings and infrastructures, under varying or hazardous environment. Spanning almost 24 years of space-borne radar observations, this study aims to perform classical PSInSAR (Persistent Scatterer Interferometric Synthetic Aperture Radar) analysis incorporating measurements of ERS, Envisat, TerraSAR-X, Sentinel-1A and Radarsat-2 satellites. The results from the processing of different sensing geometries over Bratislava (Slovakia) urban area are presented, focusing on the description of characteristics associated with the specifics of every satellite platform in use. The discussion over technical feasibility of infrastructure monitoring is accompanied by the outline of possible future needs for the utilisation of the wealth source of information provided by the satellite radar imagery.

Environmental monitoring of natural disasters using synthetic aperture radar image multi-directional characteristics

Synthetic aperture radar (SAR) imaging can penetrate rain, snow, fog, and mist, providing an effective precipitation detector under a wide range of atmospheric conditions. This powerful technique operates under any weather conditions and can detect environmental changes or help in the evaluation of natural disasters. Multi-temporal change detection is important for monitoring disasters, but commonly applied wavelet transforms are not ideal for capturing change information. This paper presents a new multi-directional change detection (MDCD) method designed to improve the change detection accuracy using multiple SAR images. This method employs double-density dual-tree complex wavelet transforms (DDDT-CWT) theory, which allows the capture of multi-directional information. The MDCD method can provide information for 16 directions at any decomposed scale, which allows for change detection in multiple SAR images collected over time. We used the MDCD algorithm method to analyse SAR images from actual natural disasters and successfully identified environmental changes over time using SAR images.

Three-dimensional positions of scattering centers reconstruction from multiple SAR images based on radargrammetry

A method and procedure is presented to reconstruct three-dimensional (3D) positions of scattering centers from multiple synthetic aperture radar (SAR) images. Firstly, two-dimensional (2D) attribute scattering centers of targets are extracted from 2D SAR images. Secondly, similarity measure is developed based on 2D attributed scatter centers’ location, type, and radargrammetry principle between multiple SAR images. By this similarity, we can associate 2D scatter centers and then obtain candidate 3D scattering centers. Thirdly, these candidate scattering centers are clustered in 3D space to reconstruct final 3D positions. Compared with presented methods, the proposed method has a capability of describing distributed scattering center, reduces false and missing 3D scattering centers, and has fewer restrictionson modeling data. Finally, results of experiments have demonstrated the effectiveness of the proposed method.

Spaceborne MIMO Synthetic Aperture Radar for Multimodal Operation

In this paper, we introduce a novel multiple-input multiple-output (MIMO) synthetic aperture radar (SAR) concept for multimodal operation. The proposed system employs waveforms based on the orthogonal frequency division multiplexing (OFDM) technique and digital beamforming (DBF) on receive. Thereby, it becomes feasible to maximize spatial degrees of freedom, which are necessary for the multimodal operation. The proposed MIMO SAR system produces multiple high-resolution wide-swath SAR imageries that are used for coherent postprocessing. Through this paper, we aim to open up a new perspective of using MIMO concept for a wide-swath SAR imaging with high resolution in interferometric and polarimetric modes, based on OFDM and DBF techniques. Therefore, this paper encompasses broad theoretical backgrounds and general system design issues as well as specific signal processing techniques and aspects.