2-D and 3-D processing and interpretation of multi-fold ground penetrating radar data:a case history from an archaeological site

A 2.5-D and 3-D multi-fold GPR survey was carried out in the Archaeological Park of Aquileia (northern Italy). The primary objective of the study was the identification of targets of potential archaeological interest in an area designated by local archaeological authorities. The second geophysical objective was to test 2-D and 3-D multi-fold methods and to study localised targets of unknown shape and dimensions in hostile soil conditions. Several portions of the acquisition grid were processed in common offset (CO), common shot (CSG) and common mid point (CMP) geometry. An 8×8 m area was studied with orthogonal CMPs thus achieving a 3-D subsurface coverage with azimuthal range limited to two normal components. Coherent noise components were identified in the pre-stack domain and removed by means of FK filtering of CMP records. Stack velocities were obtained from conventional velocity analysis and azimuthal velocity analysis of 3-D pre-stack gathers. Two major discontinuities were identified in the area of study. The deeper one most probably coincides with the paleosol at the base of the layer associated with activities of man in the area in the last 2500 years. This interpretation is in agreement with the results obtained from nearby cores and excavations. The shallow discontinuity is observed in a part of the investigated area and it shows local interruptions with a linear distribution on the grid. Such interruptions may correspond to buried targets of archaeological interest. The prominent enhancement of the subsurface images obtained by means of multi-fold techniques, compared with the relatively poor quality of the conventional single-fold georadar sections, indicates that multi-fold methods are well suited for the application to high resolution studies in archaeology.

Remote sensing of ocean waves by along-track interferometric synthetic aperture radar

A new wave retrieval method for the Along-Track Interferometric Synthetic Aperture Radar (AT-InSAR) phase image is presented. The new algorithm, named parametric retrieval algorithm (PRA), uses the full nonlinear mapping relations. It differs from previous retrieval algorithms in that it does not require a priori information about the sea state or the wind vector from scatterometer data. Instead, it combines the observed AT-InSAR phase spectrum and assumed wind vector to estimate the wind sea spectrum. The method has been validated using several C-band and X-band HH-polarized AT-InSAR observations collocated with spectral buoy measurements. In this paper, X-band and C-band HH-polarized AT-InSAR phase images of ocean waves are first used to study AT-InSAR wave imaging fidelity. The resulting phase spectra are quantitatively compared with forward-mapped in situ directional wave spectra collocated with the AT-InSAR observations. Subsequently, we combine the parametric retrieval algorithm (PRA) with X-band and C-band HH-polarized AT-InSAR phase images to retrieve ocean wave spectra. The results show that the ocean wavelengths, wave directions, and significant wave heights estimated from the retrieved ocean wave spectra are in agreement with the buoy measurements.

Near-Space Wide-Swath Radar Imaging With Multiaperture Antenna

Near-space, defined as the altitude region between 20 and 100 km, offers many capabilities that are not accessible for low Earth-orbit (LEO) satellites or airplanes because it is above storm and not constrained by orbital mechanics and high fuel consumption. Hence, a high flying speed can be obtained for the maneuvering vehicles operating in near-space. This offers a promising solution to simultaneous high-resolution and wide-swath synthetic aperture radar (SAR) imaging. As such, one near-space wide-swath SAR imaging technique is presented in this letter. The system configuration, signal model, and imaging scheme are described. An example near-space SAR system is designed, and its imaging performance is analyzed. Simulation results show that near-space maneuvering vehicle SAR indeed seems to be a promising solution to wide-swath SAR imaging.

On the nonlinear integral transform of an ocean wave spectrum into an along-track interferometric synthetic aperture radar image spectrum

We present a new nonlinear integral transform relating the ocean wave spectrum to the along-track interferometric synthetic aperture radar (AT-INSAR) image spectrum. The AT-INSAR, which is a synthetic aperture radar (SAR) employing two antennas displaced along the platform's flight direction, is considered to be a better instrument for imaging ocean waves than the SAR. This is because the AT-INSAR yields the phase spectrum and not only the amplitude spectrum as with the conventional SAR. While the SAR and AT-INSAR amplitude spectra depend strongly on the modulation of the normalized radar cross section (NRCS) by the long ocean waves, which is poorly known, the phase spectrum depends only weakly on this modulation. By measuring the phase difference between the signals received by both antennas, AT-INSAR measures the radial component of the orbital velocity associated with the ocean waves, which is related to the ocean wave height field by a well-known transfer function. The nonlinear integral transform derived in this paper differs from the one previously derived by Bao et al. [1999] by an additional term containing the derivative of the radial component of the orbital velocity associated with the long ocean waves. By carrying out numerical simulations, we show that, in general, this additional term cannot be neglected. Furthermore, we present two new quasi-linear approximations to the nonlinear integral transform relating the ocean wave spectrum to the AT-INSAR phase spectrum.

Ocean wave imaging mechanism by imaging radar

Analytical representations of the high frequency spectra of ocean wave and its variation due to the variation of ocean surface current are derived from the wave-number spectrum balance equation. The ocean surface imaging formulation of real aperture radar (RAR) is given using electromagnetic wave backscattering theory of ocean surface and the modulations of ocean surface winds, currents and their variations to RAR are described. A general representation of the phase modulation induced by the ocean surface motion is derived according to standard synthetic aperture radar (SAR) imaging theory. The detectability of ocean current and sea bottom topography by imaging radar is discussed. The results constitute the theoretical basis for detecting ocean wave fields, ocean surface winds, ocean surface current fields, sea bottom topography, internal wave and so on.

A parametric method of retrieving ocean wave spectra from synthetic aperture radar images

A parametric method that extracts the ocean wave directional spectra from synthetic aperture radar (SAR) image is presented. The 180 degrees ambiguity of SAR image and the loss of information beyond the azimuthal cutoff can be overcome with this method. The ocean wave spectra can be obtained from SAR image directly by using iteration inversion mapping method with forward nonlinear mapping. Some numerical experiments have been made by using ERS-1 satellite SAR imagette data. The ocean wave direction retrieved from SAR imagette data is in agreement with the wind direction from the scatterometer data.

Numerical simulation and validation of ocean waves measured by an Along-Track Interferometric Synthetic Aperture Radar

A new nonlinear integral transform of ocean wave spectra into Along-Track Interferometric Synthetic Aperture Radar (ATI-SAR) image spectra is described. ATI-SAR phase image spectra are calculated for various sea states and radar configurations based on the nonlinear integral transform. The numerical simulations show that the slant range to velocity ratio (R/V), significant wave height to ocean wavelength ratio (H-s/lambda), the baseline (2B) and incident angle (theta) affect ATI-SAR imaging. The ATI-SAR imaging theory is validated by means of Two X-band, HH-polarized ATI-SAR phase images of ocean waves and eight C-band, HH-polarized ATI-SAR phase image spectra of ocean waves. It is shown that ATI-SAR phase image spectra are in agreement with those calculated by forward mapping in situ directional wave spectra collected simultaneously with available ATI-SAR observations. ATI-SAR spectral correlation coefficients between observed and simulated are greater than 0.6 and are not sensitive to the degree of nonlinearity. However, the ATI-SAR phase image spectral turns towards the range direction, even if the real ocean wave direction is 30 degrees. It is also shown that the ATI-SAR imaging mechanism is significantly affected by the degree of velocity bunching nonlinearity, especially for high values of R/V and H-s/lambda.

A new method for the estimation of oceanic mixed-layer depth using shipboard X-band radar images

Shipboard X-band radar images acquired on 24 June 2009 are used to study nonlinear internal wave characteristics in the northeastern South China Sea. The studied images show three nonlinear internal waves in a packet. A method based on the Radon Transform technique is introduced to calculate internal wave parameters such as the direction of propagation and internal wave velocity from backscatter images. Assuming that the ocean is a two-layer finite depth system, we can derive the mixed-layer depth by applying the internal wave velocity to the mixed-layer depth formula. Results show reasonably good agreement with in-situ thermistor chain and conductivity-temperature-depth data sets.

Measurements of ocean wave and current field using dual polarized X-band radar

A new ocean wave and sea surface current monitoring system with horizontally-(HH) and vertically-(VV) polarized X-band radar was developed. Two experiments into the use of the radar system were carried out at two sites, respectively, for calibration process in Zhangzi Island of the Yellow Sea, and for validation in the Yellow Sea and South China Sea. Ocean wave parameters and sea surface current velocities were retrieved from the dual polarized radar image sequences based on an inverse method. The results obtained from dual-polarized radar data sets acquired in Zhangzi Island are compared with those from an ocean directional buoy. The results show that ocean wave parameters and sea surface current velocities retrieved from radar image sets are in a good agreement with those observed by the buoy. In particular, it has been found that the vertically-polarized radar is better than the horizontally-polarized radar in retrieving ocean wave parameters, especially in detecting the significant wave height below 1.0 m.

Processing of Synthetic Aperture Radar Images by a Multiscale Boundary Contour System and Feature Contour System

An improved Boundary Contour System (BCS) and Feature Contour System (FCS) neural network model of preattentive vision is applied to large images containing range data gathered by a synthetic aperture radar (SAR) sensor. The goal of processing is to make structures such as motor vehicles, roads, or buildings more salient and more interpretable to human observers than they are in the original imagery. Early processing by shunting center-surround networks compresses signal dynamic range and performs local contrast enhancement. Subsequent processing by filters sensitive to oriented contrast, including short-range competition and long-range cooperation, segments the image into regions. The segmentation is performed by three "copies" of the BCS and FCS, of small, medium, and large scales, wherein the "short-range" and "long-range" interactions within each scale occur over smaller or larger distances, corresponding to the size of the early filters of each scale. A diffusive filling-in operation within the segmented regions at each scale produces coherent surface representations. The combination of BCS and FCS helps to locate and enhance structure over regions of many pixels, without the resulting blur characteristic of approaches based on low spatial frequency filtering alone.

Processing of Synthetic Aperture Radar Images by the Boundary Contour System and Feature Contour System

An improved Boundary Contour System (BCS) and Feature Contour System (FCS) neural network model of preattentive vision is applied to two large images containing range data gathered by a synthetic aperture radar (SAR) sensor. The goal of processing is to make structures such as motor vehicles, roads, or buildings more salient and more interpretable to human observers than they are in the original imagery. Early processing by shunting center-surround networks compresses signal dynamic range and performs local contrast enhancement. Subsequent processing by filters sensitive to oriented contrast, including short-range competition and long-range cooperation, segments the image into regions. Finally, a diffusive filling-in operation within the segmented regions produces coherent visible structures. The combination of BCS and FCS helps to locate and enhance structure over regions of many pixels, without the resulting blur characteristic of approaches based on low spatial frequency filtering alone.

Human respiration rate estimation using ultra-wideband distributed cognitive radar system

It has been shown that remote monitoring of pulmonary activity can be achieved using ultra-wideband (UWB) systems, which shows promise in home healthcare, rescue, and security applications. In this paper, we first present a multi-ray propagation model for UWB signal, which is traveling through the human thorax and is reflected on the air/dry-skin/fat/muscle interfaces. A geometry-based statistical channel model is then developed for simulating the reception of UWB signals in the indoor propagation environment. This model enables replication of time-varying multipath profiles due to the displacement of a human chest. Subsequently, a UWB distributed cognitive radar system (UWB-DCRS) is developed for the robust detection of chest cavity motion and the accurate estimation of respiration rate. The analytical framework can serve as a basis in the planning and evaluation of future measurement programs. We also provide a case study on how the antenna beamwidth affects the estimation of respiration rate based on the proposed propagation models and system architecture