The sidelobe canceler in radar systems is a highly computational demanding problem. It can be efficiently tackled by resorting to the QR decomposition mapped onto a systolic array processor. The paper reports several mapping strategies by using massive parallel computers available on the market. MIMD as well as SIMD machines have been used, specifically MEIKO Computing Surface, nCUBE2, Connection Machine CM-200, and MasPar MP-1. The achieved data throughput values have been measured for a number of operational situations of practical interest.

Radar signals fluctuate because of the incoherent scattering of raindrops. Dual-polarization radar estimates rainfall rates from differential reflectivity (ZDR) and horizontal reflectivity (ZH). Here, ZDR and ZH are extracted from fluctuating radar signals by averaging. Therefore, instrumentally measured ZDR and ZH always have errors, so that estimated rainfall rates also have errors. This paper evaluates rainfall rate errors caused by signal fluctuation. Computer simulation based on a physical raindrop model is used to investigate the standard deviation of rainfall rate. The simulation considers acquisition time, and uses both simultaneous and alternate sampling of horizontal and vertical polarizations for square law and logarithmic estimators at various rainfall rates and elevation angles. When measuring rainfall rates that range from 1.0 to 10.0mm/h with the alternate sampling method, using a logarithmic estimator at a relatively large elevation angle, the estimated rainfall rates have significant errors. The simultaneous sampling method is effective in reducing these errors.

A rigorous radar cross section (RCS) analysis of a two-dimensional parallel-plate waveguide cavity with three-layer material loading is carried out for the E- and H-polarized planc wave incidence using the Wiener-Hopf technique. Introducing the Fourier transform for the scattered field and applying boundary conditions in the transform domain, the problem is formulated in terms of the simultaneous Wiener-Hopf equations satisfied by the unknown spectral functions. The Wiener-Hopf equations are solved via the factorization and decomposition procedure together with rigorous asymptotics, leading to the efficient approximate solution. The scattered field in the real space is evaluated by taking the inverse Fourier transform and applying the saddle point method. Representative numerical examples on the RCS are given for various physical parameters. It is shown that the three-layer lossy material loading inside the cavity results in significant RCS reduction over broad frequency range.

The techniques for imaging optically opaque region using an electromagnetic wave radar are being developed. One important application of these techniques is the detection of buried pipes and cables. The image quality of subsurface radar often becomes low because the electromagnetic waves are affected by the attenuation and inhomogeneity of soil. Hence, a method which improves the quality of the radar images has been required. The migration method is utilized in reflective seismic processing and is derived based on the solution of the wave equation represented in spatial frequency domain. It is classified into the F-K and the phase-shift (P-S) migration method. The former is derived on the assumption that propagation velocity of the wave is uniform in the soil while the latter is assumed that the propagation velocity is varying depending on the depth from the ground surface. The P-S method gives relatively good quality images but it requires very long computation time. In this paper, we propose the block migration method in which the F-K method is applied to the divided image blocks with local propagation velocity. In order to solve a problem concerning the connection between the contiguous blocks we present two approaches which are the processings using the overlapped regions and the Lapped Orthogonal Transform (LOT). Some experimental results point out that the block migration method has a good capability of improving the image quality and the processing time using LOT becomes one tenth in comparison with the P-S method.

The plane wave diffraction by a two-dimensional parallel-plate waveguide cavity with partial material loading is rigorously analyzed for both the E and the H polarization using the Wiener-Hopf technique. Introducing the Fourier transform for the scattered field and applying boundary conditions in the transform domain, the problem is formulated in terms of the simultaneous Wiener-Hopf equations satisfied by the unknown spectral functions. The Wiener-Hopf equations are solved exactly via the factorization and decomposition procedure leading to the formal solution, which involves branch-cut integrals with unknown integrands as well as infinite series with unknown coefficients. Applying rigorous asymptotics with the aid of the edge condition, the approximate solution to the Wiener-Hopf equations is derived in the form suitable for numerical computations. The scattered field inside and outside the cavity is evaluated by taking the inverse Fourier transform together with the use of the saddle point method. Numerical examples of the radar cross section are presented for various physical parameters, and the far field backscattering characteristics of the cavity are discussed in detail. Some comparisons with a high-frequency technique are also given to validate the present method.

We developed an automatic data processing algorithm for a ground-probing radar which is essential in analyzing a large amount of data by a non-expert. Its aim is to obtain an optimum result that the conventional technique can give, without the assistance of an experienced operator. The algorithm is general except that it postulates the existence of at least one isolated target in the radar image. The raw images of underground objects are compressed in the vertical and the horizontal directions by using a pulse-compression filter and the aperture synthesis technique, respectively. The test function needed to configure the compression filter is automatically selected from the given image. The sensitivity of the compression filter is adjusted to minimize the magnitude of spurious responses. The propagation velocity needed to perform the aperture synthesis is determined by fitting a hyperbola to the selected echo trace. We verified the algorithm by applying it to the data obtained at two test sites with different magnitude of clutter echoes.

Many previous works state that a multiple Sidelobe canceller (MSLC) with two auxiliary antennas is successful in suppressing two interference signals received simultaneously by sidelobes of a main antenna. In this paper, we show that the MSLC does not always guarantee such capability in three dimensional applications where the incident direction of interference signals is defined by two angles (elevation and azimuth). We show the singularity of the autocorrelation matrix for the auxiliary channel signals induces the degradation of the capability by analyzing characteristics of MSLC's in three dimensional applications from the view point of the eigenvalue problem. To overcome this singularity, we propose a novel MSLC controlling the placement of auxiliary antennas by means of switching over three antennas arranged triangularly. Some simulations are conducted to show the effectiveness of the proposed MSLC.

This paper applies the principle of radar polarimetry to the synthetic aperture frequency modulated continuous wave radar. First, the principle of monochromatic wave radar polarimetry using scattering matrix and polarization ratio necessary for introducing polarimetric imaging is given. In order to accommodate this principle to a wideband radar, a scattering matrix must be introduced, because FM-CW radar utilizes a wideband signal. This paper points out that the polarimetric target reflection coefficient obtained by the synthetic aperture FM-CW radar works as the scattering matrix element. This replacement, i.e., polarimetric reflection coefficient = the scattering matrix element, was verified by an experiment based on the polarization ratio which maximizes and minimizes a target. A radar system operative in the microwave X-band was successfully applied to the polarimetric detection of a metallic pipe of different orientations, demonstrating the validity of FM-CW radar polarimetry, and indicating an establishment of full polarimetric radar system.

A two-dimensional quasi-exact active imaging method for detecting the conducting objects buried in a dielectric half-space is proposed. In this imaging method, an image function which is a projection of buried object to an arbitrary direction, is introduced exactly by taking account of the presence of the planar boundary. The image function is synthesized from the scattering fields which are measured by moving a transmitting antenna (a current source) and a receiving antenna (an observation point) simultaneously along the ground surface. The scattering field is generated by the physical optics current assumed on the surface of buried object. Because the effectiveness of physical optics approximation has been confirmed for this problem, this is a quasi-exact active imaging method. The validity of this imaging method is confirmed by some numerical simulations and an experiment.

The size and weight of marine pulse radar systems must be limited in order to mount them on board boats. However, the azimuthal resolution of a marine radar with a small antenna is degraded by the antenna beam width. It is desirable to use signal processing techniques to increase both the azimuthal resolution and the range resolution of such systems without changing their external configuration. This paper introduces a resolution enhancement method based on deconvolution, which is a kind of inversion. The frequency domain deconvolution method is described first. The effectiveness of the proposed method is shown by simulation. Then, an example of resolution enhancement processing is applied to a pulse radar. The results of practical experiments show that this method is a promising way of upgrading radars by simply processing the received signals.

This paper presents a collision warning system that uses laser radar to measure the distance to a preceding vehicle and issues an audible warning to alert the driver if a safe headway is not maintained. The laser radar system is of the cooperative type in that it detects light reflected from a reflex reflector attached at the rear of other vehicles. With a 10-watt pulsed laser, a maximum detection range of over 100m is obtained. The construction and operation of the collision warning system are described along with the configuration of the optical system used in the laser radar head and the results of detection performance evaluations.

This paper proposes that the statistical property of the wave form obtained by a pulse type subsurface radar follows the Weibull probability density distribution. The shape parameter of this distribution is related to the underground condition. By using the shape parameter, we calculated the statistical variance. The ratio of the variance of target area to that of non-target area in invisible medium is evaluated for the effect of the radar signal processing. Over 20dB improvement, for example, can be obtained by means of Log/CFAR processing. It made clear that the cell size of processing should be selected the length corresponding to self-correlation.

An FM-CW radar system for the detection of objects buried in sandy ground is explored and is applied to a field measurement. The key factors for underground FM-CW radar performance are the center frequency and bandwidth determining the depth at which the radar can detect targets and the resolution in the range direction. In order for FM-CW radar sounding, two ridged horn antennas are employed in the system, which are operative in the frequency range of 250-1000MHz. The impedance matching to the ground is optimized by measuring the echo strength from a fixed target as a function of the spacing interval between the antenna aperture and the ground surface. It is shown that the radar with an output power of 18dBm could detect a metallic plate (30100cm) and a pipe (10cmφ) buried at the depth of 1.2m. Also the synthetic aperture technique together with an averaging and subtracting method produced fine image in shallow region up to 100cm in the sandy ground.

In the ISAR (Inverse Synthetic Aperture Radar), when a target is to be recognized by use of the radar image produced from the radar echoes, it is important first to estimate the scale of the target. To estimate the scale, the rotating motion of the target must be estimated. This paper describes a method for estimating the scale of the target from the information on the radar image by converting the target figure into a simple model and estimating the rotating motion of the target.

This paper discusses methods and numerical simulations of one and two dimensional profilings for an arbitrary convex conducting target using the electromagnetic backscattering. The inversions for profile reconstructions are based upon the modified extended physical optics method (EPO). The modified EPO method assumes the modified physical optics current properly over the entire surface of conducting scatterers. First, the cross sectional area along a line of sight is reconstructed by performing iteratively the Fourier transform of the backscattering field in the frequency domain. Second, the two dimensional profile is reconstructed by synthesizing the above one dimensional results for several incident angles. Numerical simulation results of the target profiling are shown for spheroids and cone-spheroid.

This paper describes a new single-unit underground radar for detecting underground buried pipes. The pipe depth can be calculated from the hyperbolic shape in the cross-sectional image of radar echoes. The edge contour of the image is extracted, and the buried depth is judged from the similarity between the extracted hyperbolic curve and the theoretical curve. A suitable amplification rate is estimated by choosing the best image from numerous cross-sectional images formed during one antenna movement repeated at different amplification rates. The best image has few pixels corresponding to weak and saturated signals. The new radar is very compact, so it can be operated by one person. Objects buried up to 2.0m deep can be detected.

A fully focused Synthetic Aperture Radar (SAR) image can be obtained only if the raw data processing procedure takes into account the space-variance of the SAR system transfer function. This paper presents a nonconventional Fast Fourier Transform (FFT) algorithm which allows an efficient compensation of the space-variant effect. It is specially designed for the SAR data of the Japanese Earth Resources Satellite (JERS-1) but can be extended to different cases.

Radar imaging technique is one of the most powerful tool for underground detection. However, performance of conventional methods is not sufficiently high when the observational direction or the aperture size is restricted. In the present paper, an image reconstruction method based on a model fitting with nonlinear least-squares has been developed, which is applicable to arbitrarily arranged arrays. Reconstruction is executed on the assumption that targets consist of discrete point scatterers embedded in a homogeneous medium. Model fitting is iterated as the number of point target in the assumed model is increased, until the residual in fitting becomes unchanged or small enough. A penalty function is used in nonlinear least-squares to make the algorithm stable. Fundamental characteristics of the method revealed with computer simulation are described. This method focuses a much sharper image than that obtained by the conventional aperture synthesis technique.

A multi-beam airborne pulsed-Doppler radar (MBR) system is presented and its clutter rejection performance compared with conventional phased array radar (PAR)'s by PRF tuning is discussed. The pulsed-Doppler radar equations taking account of the multi-beam operation are introduced and some kinds of computer simulations for seeking the conditions to get maximum signal to clutter ratio are carried out. As a results of this, it is cleared that same order of signal to clutter ratio improvement gotten in high PRF operation by conventional PAR can be realized at lower PRF operation by MBR on clutter free area, and higher clutter rejection effect, which is proportional to beam numbers, is obtained under affection of both of mainlobe and sidelobe clutters with order of beam numbers. This also means observable numbers of range bin are increased in MBR operation.

The antennas for subsurface radar are usually covered with a conducting cavity to prevent the radiation field from affecting the electromagnetic environment and to protect the received field from external noises. Furthermore, radiowave absorber is attached to the interior wall of the cavity in order to suppress the multiple reflections in the cavity. In this paper, the characteristics of the two-dimensional cavity-backed antenna having the absorber and the over-all properties of this subsurface radar due to buried objects are numerically analyzed by the Finite-Difference Time-Domain method. It is shown that the pulse propagation in the ground is confined to the narrow region due to the cavity. It is also shown that the multiple reflections in the cavity are effectively suppressed by choosing the suitable absorber, and so that the distinctive pulse echo can be obtained.