Minimum-Entropy-Based Autofocus Algorithm for SAR Data Using Chebyshev Approximation and Method of Series Reversion, and Its Implementation in a Data Processor

Focusing bistatic synthetic aperture radar (SAR) data in frequency domain requires two-dimensional (2D) point target reference spectrum (PTRS). Lofield's bistatic formula (LBF) and the Method of Series Reversion (MSR) have been introduced recently to compute PTRS of bistatic SAR. In this paper, flrstly we generalize the original LBF (OLBF) by introducing the Doppler contribution functions of transmitter and receiver. Thus, OLBF and its derivatives (e.g., extended LBF) can be viewed as special forms of the generalized LBF with constant Doppler contributions. Based on this, secondly the ideal LBF (ILBF) with no computing error, except the error resulting from the principle of stationary phase, is also presented. The ILBF reveals that the theoretical PTRS of bistatic SAR consists of only two monostatic terms, but it does not include bistatic deformation term in comparison with OLBF. It supplies us with a target when we deduce the PTRS for bistatic SAR. Finally, to get the precise analytical PTRS for general bistatic SAR, an approximated ILBF (AILBF) is proposed. It expresses the Doppler contributions of the transmitter and receiver as power series and can approach the ILBF very well. AILBF can keep the precision as MSR and inherit a simple form from LBF. In addition, error limit for the validity of bistatic PTRS is also given. The results in this paper can be used to develop imaging algorithms for extreme bistatic (e.g., spaceborne/airborne) and high squint (e.g., bistatic forward-looking) cases.

The speed and direction of a missile shifts sharply in the dive phase, making the azimuth frequency modulation (FM) rate change with the azimuthal position, leading to azimuth ambiguities and image distortion. To solve this problem, a modified nonlinear chirp scaling (NLCS) algorithm was adopted to compensate for the azimuth FM rate. First, the geometric configuration and echo signal model of the spaceborne missile bistatic synthetic aperture radar (SAR) were built, and then the Doppler frequency correction was performed, and the 2-D spectrum of the signal was derived by the method of series reversion. Next, range migration correction and range compression were finished in the 2-D frequency domain. Following this, a modified NLCS algorithm was proposed to solve the space variance of Doppler phase problem. After compensating for the azimuth FM rate, the azimuth compression focusing was completed and the imaging result was obtained. Finally, by comparing the calculation amount, imaging effect, and performance index with the traditional NLCS algorithm, it can be concluded that the algorithm reduced the calculation amount by 1.0128 × 108 floating point operations per second (FLOPs) compared with the traditional NLCS algorithm, and the azimuth focusing effect of the edge point was greatly improved. Its resolution, peak sidelobe ratio (PSLR), and integrated sidelobe ratio (ISLR) were improved by 0.87 m, 3.32 dB, and 1.79 dB, respectively, which proved the effectiveness and feasibility of this method.

One of the problems associated with synthetic aperture radar (SAR) data analysis is the presence of random noise or speckle SAR data, being achromatic in nature, which offers very limited scope for the detection and delineation of various terrain features. ERS-1 SAR data for the coastal region of West Bengal, India were processed (a) to suppress the random noise using various filters, (b) to generate the intensity, hue and saturation (IHS) transform from temporal SAR data, and (c) to study the synergism of SAR data with optical sensor data. The results indicate that the Gamma MSP filter with a 5 5 pixel kernel size has been the most efficient in suppressing the noise and concurrently improving the image contrast. The IHS transform of temporal SAR data made it easier to discriminate between various wetland categories. This was also the case with hybrid image generated by the Indian Remote Sensing Satellite (IRS-1B) Linear Imaging and Self-scanning Sensor (LISS-II) data when compared to SAR data alone.

For a synthetic aperture radar (SAR) onboard a platform with a rectilinear track, the range history of a point target can be accurately expressed hyperbolically. The track can be curvilinear for a maneuverable SAR platform. The hyperbolic equation becomes inadequate, and an expression with high-order terms is needed. Using the method of series reversion, we derived the 2-D spectrum for the return signal of the curvilinear SAR. There were five independent variables in the spectrum, but available imaging algorithms could only handle three in the focusing using the spectrum. Thus, a variable-decoupling method was developed to reparameterize the initial spectrum so that only three variables remained. After the incorporation of the decoupling method into the chirp-scaling algorithm, simulations of the SAR with a curvilinear track were studied. Promising results were obtained.

Despite recent research on the potential of dual- (DP) and full-polarimetry (FP) Synthetic Aperture Radar (SAR) data for crop mapping, the capability of compact polarimetry (CP) SAR data has not yet been thoroughly investigated. This is of particular concern, given the availability of such data from RADARSAT Constellation Mission (RCM) shortly. Previous studies have illustrated potential for accurate crop mapping using DP and FP SAR features, yet what contribution each feature makes to the model accuracy is not well investigated. Accordingly, this study examined the potential of the early- to mid-season (i.e., May to July) RADARSAT-2 SAR images for crop mapping in an agricultural region in Manitoba, Canada. Various classification scenarios were defined based on the extracted features from FP SAR data, as well as simulated DP and CP SAR data at two different noise floors. Both overall and individual class accuracies were compared for multi-temporal, multi-polarization SAR data using the pixel- and object-based random forest (RF) classification schemes. The late July C-band SAR observation was the most useful data for crop mapping, but the accuracy of single-date image classification was insufficient. Polarimetric decomposition features extracted from CP and FP SAR data produced relatively equal or slightly better classification accuracies compared to the SAR backscattering intensity features. The RF variable importance analysis revealed features that were sensitive to depolarization due to the volume scattering are the most important FP and CP SAR data. Synergistic use of all features resulted in a marginal improvement in overall classification accuracies, given that several extracted features were highly correlated. A reduction of highly correlated features based on integrating the Spearman correlation coefficient and the RF variable importance analyses boosted the accuracy of crop classification. In particular, overall accuracies of 88.23%, 82.12%, and 77.35% were achieved using the optimized features of FP, CP, and DP SAR data, respectively, using the object-based RF algorithm.

Method for image classification considering Probability Density Function (PDF) based on simplified beta distributions is proposed. In this paper, image classification for Synthetic Aperture Radar (SAR) data is concerned. In particular, Probability Density Function (PDF) of SAR data is followed by not multivariate normal distribution but Chi-Square like distribution. It, however, is not always true that the PDF of SAR data is followed by Chi-Square distribution. Due to the mismatch between Chi-Square distribution and actual distribution, classification performance gets worth. In this paper, simplified beta distribution is assumed for the PDF of the SAR data. Furthermore, it is used to add texture information to the SAR data when the Maximum Likelihood classification is applied. In the paper, “Contrast” of texture feature is added to the SAR data. Through the experiments with real SAR data, it is found that matching error between real PDF and the proposed simplified beta distribution is smaller than the normal distribution. It is also found that applying the proposed distribution-adaptive maximum likelihood method using the simplified beta-distribution could achieve a classification accuracy improvement of 94.7% and 12.1%.

In geosynchronous (GEO) synthetic aperture radar (SAR), the high orbit height and the long integration time would lead to the significant round-trip delay time and the curved flight path, thus the traditional range model with ‘stop-go’ assumption will be ineffective. In this Paper, based on the analysis of instantaneous Doppler frequency of GEO SAR, an equivalent midpoint mathematical model is proposed to accurately describe the range history, which takes the effect of ‘stop-go’ assumption into account and is more applicable to SAR processing. Based on the equivalent model, the two dimensional (2-D) spectrum of the GEO SAR can be obtained easily by the method of series reversion (MSR). The simulations show the validity of the presented method.

Since the launch of the European Space Agency’s (ESA’s) Sentinel-1 mission, the synthetic aperture radar (SAR) user community has grown exponentially, adding a large number of users to the field who have only limited experience with the peculiarities of SAR and interferometric SAR (InSAR) data. This growth has been fueled by the increasing availability of globally acquired, free and open SAR data and by a growing number of freely accessible, open source processing tools, such as the Sentinel Application Platform, InSAR Scientific Computing Environment, GMTSAR, and the Alaska Satellite Facility’s cloud-based Hybrid Pluggable Processing Pipeline and OpenScienceLab platforms. <p xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink">Anticipated to launch in early 2025, the NASA–Indian Space Research Organisation SAR (NISAR) mission will amplify this trend by bringing an unprecedented amount of new, free and open SAR data to the international SAR science community. As the first mission to provide L-band data freely, on a global scale, and in a regularly sampled manner, NISAR will further diversify the SAR application portfolio, likely accelerating community growth. Regional access to simultaneous L- and S-band data facilitated by NISAR will furthermore enable new SAR technology development that may lead to exciting new scientific discoveries. <p xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink">Focusing on efforts made by NASA and the broader U.S. SAR community, this article summarizes innovative approaches to making the massive datasets of Sentinel-1 and NISAR more accessible to and usable by the SAR community. We summarize novel concepts for data discovery and distribution and highlight tools and services that enable working with SAR data directly at the archive. We also discuss recent education and training efforts developed across the U.S. SAR community to help new users of SAR build capacity to use this exciting dataset. Here, we introduce a set of independent, yet interrelated, projects aimed at developing educational materials and open source science workflows for SAR users with different experience levels.

The paper provides a review into the steps involved in acquiring and processing synthetic aperture radar (SAR) data transmitted by the European remote sensing (ERS) satellites. The paper reports on a simple MATLAB-based SAR processing system, that reads the image out of the complex SAR data files and that is suitable for use in the classroom to demonstrate one of the procedures used in SAR data processing. The paper can also help a beginner in the field of SAR signal processing to get information and understand the basics that are necessary to acquire and process a SAR image. The data, provided by the Alaskan Satellite Facility (ASF), is categorized into different levels and the paper describes the process of obtaining the level-1 basic image from the level-0 raw data file provided by ASF.

The systematic evaluation of synthetic aperture radar (SAR) data analysis tools, such as segmentation and classification algorithms for geographic information systems, is difficult given the unavailability of ground-truth data in most cases. Therefore, testing is typically limited to small sets of pseudoground-truth data collected manually by trained experts, or primitive synthetic sets composed of simple geometries. To address this issue, we investigate the potential of employing an alternative approach, which involves the synthesis of SAR data and corresponding label fields from real SAR data for use as a reliable evaluation testbed. Given the scale-dependent nonstationary nature of SAR data, a new modeling approach that combines a resolution-oriented hierarchical method with a region-oriented binary tree structure is introduced to synthesize such complex data in a realistic manner. Experimental results using operational RADARSAT SAR sea-ice data and SIR-C/X-SAR land-mass data show that the proposed hierarchical approach can better model complex nonstationary scale structures than local MRF approaches and existing nonparametric methods, thus making it well suited for synthesizing SAR data and the corresponding label fields for potential use in the systematic evaluation of SAR data analysis tools.

The real-time processing of synthetic aperture radar (SAR) data has a high requirement for the processor, which is a difficult problem in SAR real-time processing. With the rapid development of optoelectronic devices, traditional electrical SAR data processing can be converted into optoelectronic processing to improve the processing speed. In this paper, a new type of optical device is proposed to improve the processing speed of SAR data. With the help of a spatial light modulator (SLM), the initial SAR signal and matched filter function are loaded on the input plane and spectrum plane of the 4f system, respectively. Using an optical lens with the function of the Fourier transform, the Fourier transform and inverse Fourier transform of the SAR signal are carried out to realize the fast imaging of SAR. In theory, the processing speed of SAR data is the speed of light. Compared with traditional methods such as the range-Doppler (RD) algorithm, it is no longer necessary to carry out a one-dimensional Fourier transform but to carry out matching filtering for the azimuth and range of the spectrum plane of 4f system at the same time. In this way, it is not necessary to introduce a cylindrical lens, only a spherical lens is needed to realize the Fourier transform imaging of SAR. Finally, a two-dimensional SAR processing optical system is built to obtain the SAR image in real time.

The objectives of this research are to develop robust methods for segmentation of multitemporal synthetic aperture radar (SAR) and optical data and to investigate the fusion of multitemporal ENVISAT advanced synthetic aperture radar (ASAR) and Chinese HJ-1B multispectral data for detailed urban land-cover mapping. Eight-date multiangle ENVISAT ASAR images and one-date HJ-1B charge-coupled device image acquired over Beijing in 2009 are selected for this research. The edge-aware region growing and merging (EARGM) algorithm is developed for segmentation of SAR and optical data. Edge detection using a Sobel filter is applied on SAR and optical data individually, and a majority voting approach is used to integrate all edge images. The edges are then used in a segmentation process to ensure that segments do not grow over edges. The segmentation is influenced by minimum and maximum segment sizes as well as the two homogeneity criteria, namely, a measure of color and a measure of texture. The classification is performed using support vector machines. The results show that our EARGM algorithm produces better segmentation than eCognition, particularly for built-up classes and linear features. The best classification result (80%) is achieved using the fusion of eight-date ENVISAT ASAR and HJ-1B data. This represents 5%, 11%, and 14% improvements over eCognition, HJ-1B, and ASAR classifications, respectively. The second best classification is achieved using fusion of four-date ENVISAT ASAR and HJ-1B data (78%). The result indicates that fewer multitemporal SAR images can achieve similar classification accuracy if multitemporal multiangle dual-look-direction SAR data are carefully selected.

Fresh water is a scarce resource in the West-African Sahel region, seasonally influenced by droughts and floods. Particularly in terms of climate change, the importance of wetlands increases for flora, fauna, human population, agriculture, livestock and fishery. Hence, access to open water is a key factor. Long-term monitoring of water dynamics is of great importance, especially with regard to the spatio-temporal extend of wetlands and drylands. It can predict future trends and facilitate the development of adequate management strategies. Lake Tabalak, a Ramsar wetland of international importance, is one of the most significant ponds in Niger and a refuge for waterbirds. Nevertheless, human population growth increased the pressure on this ecosystem, which is now degrading for all uses. The main objective of the study is a long-term monitoring of the Lake Tabalak’s water dynamics to delineate permanent and seasonal water bodies, using weather- and daytime-independent multi-sensor and multi-temporal Synthetic Aperture Radar (SAR) data available for the study area. Data of the following sensors from 1993 until 2016 are used: Sentinel-1A, TerraSARX, ALOS PALSAR-1/2, Envisat ASAR, RADARSAT-1/2, and ERS-1/2. All SAR data are processed with the Multi-SAR-System, unifying the different characteristics of all above mentioned sensors in terms of geometric, radiometric and polarimetric resolution to a consistent format. The polarimetric representation in Kennaugh elements allows fusing single-polarized data acquired by older sensors with multi-polarized data acquired by current sensors. The TANH-normalization guarantees a consistent and therefore comparable description in a closed data range in terms of radiometry. The geometric aspect is solved by projecting all images to an earth-fixed coordinate system correcting the brightness by the help of the incidence angle. The elevation model used in the geocoding step is the novel global model produced by the TanDEM-X satellite mission. The advantage of the Multi-SAR-System is that it comprises ortho-rectification, radiometric enhancement, normalization and Kennaugh decomposition, independent from sensors, modes, polarizations or acquisition date of SAR data. In addition, optical satellite data can be included as well, to fill gaps where SAR data are missing due to the special normalization scheme. This kind of pre-processing is exclusively implemented at the Earth Observation Center of the German Aerospace Center in Oberpfaffenhofen, Germany. Therefore, the dynamic change of the open water of the Lake Tabalak could be classified over dry and rainy seasons and years, using different SAR data. The study provides a unique database and contributes to a better understanding of wetland systems in the Sahel region influenced by human pressure and climate change.

Recent drastic reductions in the Arctic sea-ice cover have raised an interest in understanding the role of sea ice in the global system as well as pointed out a need to understand the physical processes that lead to such changes. Satellite remote-sensing data provide important information about remote ice areas, and Synthetic Aperture Radar (SAR) data have the advantages of penetration of the omnipresent cloud cover and of high spatial resolution. A challenge addressed in this paper is how to extract information on sea-ice types and sea-ice processes from SAR data. We introduce, validate and apply geostatistical and statistical approaches to automated classification of sea ice from SAR data, to be used as individual tools for mapping sea-ice properties and provinces or in combination. A key concept of the geostatistical classification method is the analysis of spatial surface structures and their anisotropies, more generally, of spatial surface roughness, at variable, intermediate-sized scales. The geostatistical approach utilizes vario parameters extracted from directional vario functions, the parameters can be mapped or combined into feature vectors for classification. The method is flexible with respect to window sizes and parameter types and detects anisotropies. In two applications to RADARSAT and ERS-2 SAR data from the area near Point Barrow, Alaska, it is demonstrated that vario-parameter maps may be utilized to distinguish regions of different sea-ice characteristics in the Beaufort Sea, the Chukchi Sea and in Elson Lagoon. In a third and a fourth case study the analysis is taken further by utilizing multi-parameter feature vectors as inputs for unsupervised and supervised statistical classification. Field measurements and high-resolution aerial observations serve as basis for validation of the geostatistical-statistical classification methods. A combination of supervised classification and vario-parameter mapping yields best results, correctly identifying several sea-ice provinces in the shore-fast ice and the pack ice. Notably, sea ice does not have to be static to be classifiable with respect to spatial structures. In consequence, the geostatistical-statistical classification may be applied to detect changes in ice dynamics, kinematics or environmental changes, such as increased melt ponding, increased snowfall or changes in the equilibrium line.

Detailed information on spatial distribution of wetland classes is crucial for monitoring this important productive ecosystem using advanced remote sensing tools and data. Although the potential of full- and dual-polarimetric (FP and DP) Synthetic Aperture Radar (SAR) data for wetland classification has been well examined, the capability of compact polarimetric (CP) SAR data has not yet been thoroughly investigated. This is of great significance, since the upcoming RADARSAT Constellation Mission (RCM), which will soon be the main source of SAR observations in Canada, will have CP mode as one of its main SAR configurations. This also highlights the necessity to fully exploit such important Earth Observation (EO) data by examining the similarities and dissimilarities between FP and CP SAR data for wetland mapping. Accordingly, this study examines and compares the discrimination capability of extracted features from FP and simulated CP SAR data between pairs of wetland classes. In particular, 13 FP and 22 simulated CP SAR features are extracted from RADARSAT-2 data to determine their discrimination capabilities both qualitatively and quantitatively in three wetland sites, located in Newfoundland and Labrador, Canada. Seven of 13 FP and 15 of 22 CP SAR features are found to be the most discriminant, as they indicate an excellent separability for at least one pair of wetland classes. The overall accuracies of 87.89%, 80.67%, and 84.07% are achieved using the CP SAR data for the three wetland sites (Avalon, Deer Lake, and Gros Morne, respectively) in this study. Although these accuracies are lower than those of FP SAR data, they confirm the potential of CP SAR data for wetland mapping as accuracies exceed 80% in all three sites. The CP SAR data collected by RCM will significantly contribute to the efforts ongoing of conservation strategies for wetlands and monitoring changes, especially on large scales, as they have both wider swath coverage and improved temporal resolution compared to those of RADARSAT-2.