Search

Multiple-input multiple-output (MIMO) technology is widely used in the field of radar imaging. Array sparse optimization reduces the hardware cost of MIMO radar, while virtual aperture and the equivalent phase center (EPC) principle simplify the radar signal model and reduce the computation and complexity of imaging algorithms. However, the application of sparse array structure and the EPC principle produces a non-negligible phase error, which affects the imaging quality. This paper simplifies the MIMO radar signal model based on the phase center approximation, analyzes the phase error generated by this method, and proposes an improved phase error correction method to solve the problem that the target cannot be well-focused at non-reference distance during imaging. In addition, this paper designs a sparse linear MIMO array with a periodic structure, which reduces the number of transmitting and receiving units, system complexity, and hardware costs. The proposed phase correction method was combined with the wavenumber domain algorithm to simulate and experiment on the designed antenna array, and good experimental results were obtained to verify the effectiveness of the proposed method.

Ground-penetrating radar (GPR) is a classic geophysical exploration method that utilizes the emission and reception of electromagnetic waves to non-destructively detect target objects in the target medium. It has been widely applied in various fields such as pipeline detection, cavity detection, and rebar detection. However, GPR systems are susceptible to environmental clutter interference, which poses challenges for data interpretation and subsequent processing. In this paper, the separability of clutter and target signal in the frequency-wavenumber (F-K) domain is validated through modeling, leading to the proposal of a comprehensive clutter removal method based on the F-K domain for complex scenarios. The direct coupling wave is initially eliminated by applying a peak matching mean subtraction filter, which avoids the artifacts. Subsequently, the F-K domain transformation is performed and surface clutter undulations are effectively removed using a method based on singular value decomposition and k-means clustering. Finally, an angle filter with Gaussian tapering at the edges is designed based on physical models to efficiently eliminate linear interference without undesired ringing interference. The commonly used clutter removal algorithms, including mean subtraction (MS), singular value decomposition (SVD), robust principal component analysis (RPCA), and traditional F-K filtering methods, are compared with the proposed algorithm on both the numerical simulated data and actual GPR data. The results from visual and quantitative analysis confirm that our proposed method is more effective than current commonly used clutter suppression algorithms. We have successfully enhanced the Signal-to-Clutter Ratio (SCR) of the GPR data, resulting in an Improvement Factor (IF) of 30.63 dB, 23.59 dB, and 30.60 dB for simulated data, experimental data, and TU1208 public data, respectively. The detection capability of buried targets is enhanced, thereby establishing a solid foundation for subsequent data interpretation and target identification.

Ultra-wideband (UWB) vital detection radar plays an important role in post-disaster search and rescue, but the vital signal acquired in practice is often submerged in noise. In this paper, an advanced signal processing algorithm based on low-rank block-sparse representation is proposed to enhance the vital signal in life detection radar applications. The preprocessed echo signal can be decomposed into low-rank and block-sparse parts. The alternate direction method (ADM) is employed to obtain the block-sparse part containing the desired vital signal. We solve the subproblems involved in the ADM method using the Douglas/Peaceman–Rachford (DR) monotone operator splitting method. The projection method is applied to accelerate the calculation. Simulation and experimental results show that the proposed algorithm outperforms existing methods in terms of output signal-to-noise ratio (SNR).

With the development of millimeter-wave radar systems, substrate-integrated waveguide (SIW) slot antennas have been widely used. This paper presents a computational approach for calculating the equivalent normalized admittance of slot in low sidelobe SIW longitudinal slot antennas. In contrast to the existing methods, this method takes into account the impact of the N-slots coupling on the admittance. Through this method, fitting functions for resonant length and conductance are obtained, which describe the relationship between the resonant length, resonant conductance, and offset of slots. Moreover, based on the two fitting functions, a low sidelobe SIW longitudinal slot linear antenna is designed for millimeter-wave radar applications at 77 GHz. This antenna consists of a transition from the standard WR-12 waveguide to the SIW and a set of radiating slots. The simulation results indicate an impedance bandwidth of 4.8% and a gain of 13.2 dBi. The sidelobe level in the H-plane is below −28 dB, which meets the desired specifications and confirms the accuracy of the presented method.

In this study, repeat-pass synthetic aperture radar interferometry (repeat-pass THz InSAR) is first extended to the terahertz band, and it has tremendous potential in the application of high-resolution three-dimensional (3D) imaging due to its shorter wavelength, larger bandwidth, and greater sensitivity to elevation variation. The super-resolution and high sensitivity of THz InSAR pose greater demands on the baseline calibration for high-precision digital elevation model (DEM) generation. To meet the elevation accuracy requirement of THz InSAR, we propose a baseline calibration method relying on the estimation of the azimuth fringe frequency (EAFF) of the interferometric phase. Initially, a model for non-parallel sampling path errors within the squint SAR repeat-pass interferometry was established, and then, we conducted the theoretical analysis of the phase errors induced by the non-parallel errors. Following this, using a reference DEM, the relationship between the fringe frequency of the error phase and the bias in the repeat-path positioning was established. This allowed the estimation of the position errors to be transformed into the frequency spectrum estimation based on the FFT, which would mitigate the impact of unknown SAR sampling positions. Ultimately, we investigated the accuracy of the proposed EAFF calibration method, and the simulation showed that it can achieve the theoretical accuracy when the correlation coefficient exceeds 0.3. Furthermore, we configured the repeat-pass THz InSAR system with the 0.3 THz stepped-frequency radar. Compared to the conventional calibration based on ground control points (GCPs), the 3D reconstruction of both a knife and a terrain model, calibrated using the proposed EAFF algorithm, demonstrated that the elevation accuracy can achieve millimeter-level precision across the entire image swath. The above results also proved the great potential of THz InSAR in high-precision 3D imaging and remote sensing.

Ground-penetrating radar (GPR) has been widely used in estimating the permittivity of mediums. The radar echo amplitude method is an important method used by GPR in this estimation, the basic step of which is to deduce the magnitude of the permittivity according to the relationship between the reflection coefficient and the permittivity. Based on the basic principle of the radar echo amplitude method, this paper proposes a full-wave inversion real-time permittivity estimation method that can be used for stepped-frequency GPR (SFGPR), which offers high efficiency, accuracy, and generalization ability. The characteristics of this method are mainly reflected in the following four aspects: Using the SFGPR system and introducing a layered media detection model, we can complete waveform compensation optimization with high precision. The distance between the antenna and the surface of the reflective medium is extracted from the time domain waveform without manual measurement, avoiding human measurement errors. The inversion of the total reflection waveform at the required height works under the principle of an electromagnetic field, eliminating the need for repeated metal plate calibration experiments and improving work efficiency and waveform accuracy. In a continuous measurement line, the total reflection waveform inversion on each measurement point can be efficiently completed, and the change of permittivity on the measurement line can be obtained. To evaluate the feasibility of the proposed method, we performed experiments on a wall of known thickness, and the permittivity estimation was basically consistent with that of the dielectric probe, physical model calculation, and wall penetration. We also successfully applied this method to the dielectric property analysis of adobe samples. The results indicate that the proposed method can help grasp the condition of a measured medium, which can ensure the accuracy of detection and improve subsequent data processing efficiency.

In this article, a novel three-dimensional (3-D) imaging method based on the range decomposing algorithm (RDA) is proposed for millimeter wave imaging. We combined it with binomial theory and we derive the theoretical formulation of RDA applied to single-input–single-output (SISO)/multiple-input–multiple-output (MIMO) array; meanwhile, its computational complexity and computational error are analyzed. Compared to the classical Fourier algorithm, such as the range migration algorithm (RMA) and the phase shift migration (PSM), the proposed algorithm can replace the time-consuming interpolation and accumulation operations with reasonable approximations and transformations offering a more efficient approach, while maintaining the image quality. In addition, a method based on RDA which is applicable to the transformation between MIMO and SISO, is proposed to further enhance the processing efficiency. Proof-of-principle simulation using echo data collected from a large number of antennas, verifies that the proposed algorithm has higher efficiency. In order to better verify the feasibility of the proposed algorithm, a scanning prototype located in the millimeter wave band is designed. The experimental results of different targets demonstrate that the proposed algorithm achieves significantly higher reconstruction efficiency when compared to the traditional algorithms.

An efficient spectrum reconstruction algorithm based on the Tikhonov regularization for terahertz (THz) synthetic aperture radar (SAR) azimuth non-uniform sampling is proposed in this article. The high bandwidth, high azimuth resolution, and high frame rate characteristics of THz SAR contribute to its wide application prospects in both military and civilian remote sensing fields. However, the higher azimuth sampling rate also leads to the more severe non-uniform sampling issues of THz SAR. Traditional methods based on the hardware adjustment of pulse repetition frequency (PRF) and simple interpolation for azimuth resampling struggle to meet the higher imaging quality requirements. The back projection algorithm (BPA) can accurately focus non-uniformly sampled data but requires significant computational resources. The algorithm proposed in this paper, which can reconstruct the wavenumber spectrum of SAR azimuth non-uniformly sampled signals, transforms the spectrum reconstruction problem into a linear equation system and solves it using Tikhonov regularization, thereby exhibiting higher computational efficiency compared to BPA. Furthermore, the proposed algorithm is derived from precise theoretical formulations and controls the solution error by utilizing a regularization parameter, leading to a superior imaging quality compared to the azimuth resampling algorithm. In this paper, an accurate spectrum reconstruction formula of non-uniform sampling signals with a finite length is derived, the influence of noise error on the solution is analyzed, and the THz SAR azimuth non-uniform sampling signals are processed from the wavenumber domain. Finally, simulation and experimental results verify the effectiveness of the proposed algorithm.

Ultra-wideband (UWB) multiple-input multiple-output (MIMO) through-wall radar is widely used in through-wall human target detection for its good penetration characteristics and resolution. However, in actual detection scenarios, weak target masking and adjacent target unresolving will occur in through-wall imaging due to factors such as resolution limitations and differences in human reflectance, which will reduce the probability of target detection. An improved U-Net model is proposed in this paper to improve the detection probability of through-wall targets. In the proposed detection method, a ResNet module and a squeeze-and-excitation (SE) module are integrated in the traditional U-Net model. The ResNet module can reduce the difficulty of feature learning and improve the accuracy of detection. The SE module allows the network to perform feature recalibration and learn to use global information to emphasize useful features selectively and suppress less useful features. The effectiveness of the proposed method is verified via simulations and experiments. Compared with the order statistics constant false alarm rate (OS-CFAR), the fully convolutional networks (FCN) and the traditional U-Net, the proposed method can detect through-wall weak targets and adjacent unresolving targets effectively. The detection precision of the through-wall target is improved, and the missed detection rate is minimized.

The group velocity of lightning electromagnetic signals plays an important role in lightning location systems using the time difference of arrival (TDOA) method. Accurate estimation of group velocity is difficult due to the space- and time-varying properties of the Earth-ionosphere waveguide. Besides, the analytical solution of the group velocity is difficult to obtain from the classic mode theory, especially when higher-order modes, anisotropic geomagnetic background, diffuse ionosphere profile, and propagation path segmentation are all taken into consideration. To overcome these challenges, a novel numerical method is proposed in this paper to estimate the group velocity of the lightning signal during ionospheric quiet periods. The well-known Long Wavelength Propagation Capability (LWPC) code is used to model the propagation of VLF/LF radio waves. Since LWPC uses a simplified ionospheric model which is unable to describe the subtle variations of ionospheric parameters over time and space, the IRI-2016 model is incorporated into the numerical modeling process to provide more accurate ionosphere parameters. Experimental results of a VLF/LF lightning location network are demonstrated and analyzed to show the effectiveness of our method. The proposed method is also applicable when there is a sudden ionospheric disturbance as long as the parameters of the ionosphere are obtained in real time by remote sensing methods.

On 1 December 2020, China’s Chang’E-5 (CE-5) probe successfully landed in the northeastern Oceanus Procellarum. This work mainly presents the results of Lunar Regolith Penetrating Radar (LRPR) equipped on the CE-5 Lander. The lunar regolith structure of the landing site from the surface to 3-m depth is unveiled by LRPR, which found that abundant rock fragments are distributed in uniform lunar regolith. The imaging result proved that the drilling and sampling process was prevented by big rocks at about 100 cm depth. On the basis of the response of lunar soil to electromagnetic (EM) wave, the EM properties of the landing site estimate that the relative dielectric constant and the loss tangent are 2.520 ± 0.186 and 0.0133 ± 0.0020, respectively.

Drift detection is one of the important detection modes in a digital ionosonde system. In this paper, a new data processing method is presented for boosting the automatic and high-quality drift measurement, which is helpful for long-term ionospheric observation, and has been successfully applied to the Chinese Academy of Sciences, Digital Ionosonde (CAS-DIS). Based on Doppler interferometry principle, this method can be successively divided into four constraint steps: extracting the stable echo data; restricting the ionospheric detection region; extracting the reliable reflection cluster, including Doppler filtering and coarse clustering analysis; and calculating the drift velocity. Ordinary wave (O-wave) data extraction, complementary code pulse compression and other data preprocessing techniques are used to improve the signal-to-noise ratio (SNR) of echo data. For the purpose of eliminating multiple echoes, the ionospheric region is determined by combining the optimal height range and detection frequencies obtained from the ionogram. Successively, Doppler filtering and coarse clustering analysis extract reliable reflection clusters. Finally, the weighting factor is brought in, and then weighted least-squares (WLS) is used to fit the drift velocity. The entire data processing process can be implemented automatically without constantly changing parameter settings due to changes in external conditions. This is the first time coarse clustering analysis has been used to extract the paracentral reflection cluster to eliminate scattered reflection points and outer reflection clusters, which further reduces the impacts of external conditions on parameter settings and improves the ability of automatic drift measurement. Compared with the previous method possessed by Digisonde Protable Sounder 4D (DPS4D), the new method can achieve comparable drift detection precision and results even with fewer reflection points. In 2021–2022, several experiments on F region drift detection were carried out in Hainan, China. Results indicate that drift velocities fitted by the new method have diurnal variation and change more gently; the trends of drift velocities fitted by the new method and the previous method are semblable; and this new method can be widely applied to digital ionosondes.

In multiple-input-multiple-output synthetic aperture radar (MIMO-SAR) systems, sparse arrays are usually applied, resulting in increased sidelobes of the point spread function. In this paper, a phase shift migration (PSM) imaging algorithm based on the explosion reflection model with modified coherent factor was proposed for sidelobe suppression in MIMO-SAR three-dimensional (3D) imaging application. By defining the virtual difference wavenumber, reconstructing the raw echo by data rearrangement in wavenumber domain, the original coherent factor algorithm operating in spatial domain can be achieved by the PSM algorithm frame in the wavenumber domain, which means two orders of magnitude increase in computational efficiency. The correctness of the theory is verified by simulation. Finally, a bistatic prototype imaging system in the 0.3 THz band was designed for the proof-of-principle experiments. The experimental results show that the proposed algorithm has a 0.948 structural similarity value to the original coherent factor back-projection algorithm (CF-BPA) which means comparable image quality with much superior efficiency.

Ultra-wideband (UWB) impulse radar is widely used for through-wall human respiration detection due to its high range resolution and high penetration capability. UWB impulse radar emits very narrow time pulses, which can directly obtain the impulse response of the target. However, the time interval between successive pulses emitted is not ideally fixed because of timing jitter. This results in the impulse response position of the same target not being fixed, but it is related to slow-time. The clutter scattered by the stationary target becomes non-stationary clutter, which affects the accurate extraction of the human respiration signal. In this paper, we propose a method for reducing timing jitter’s impact in through-wall human detection by UWB impulse radar. After the received signal is processed by the Fast Fourier transform (FFT) in slow-time, we model the range-frequency matrix in the frequency domain as a superposition of the low-rank representation of jitter-induced clutter data and the sparse representation of human respiratory data. By only extracting the sparse component, the impact of timing jitter in human respiration detection can be reduced. Both numerical simulated data and experimental data demonstrate that our proposed method can effectively remove non-stationary clutter induced by timing jitter and improve the accuracy of the human target signal extraction.

Terahertz (THz) imaging technology has received increased attention in recent years and has been widely applied, whereas the three-dimensional (3D) imaging for moving targets remains to be solved. In this paper, an adaptive 3D imaging scheme is proposed based on a single input and multi-output (SIMO) interferometric inverse synthetic aperture radar (InISAR) imaging system to achieve 3D images of moving targets in THz band. With a specially designed SIMO antenna array, the angular information of the targets can be determined using the phase response difference in different receiving channels, which then enables accurate tracking by adaptively adjusting the antenna beam direction. On the basis of stable tracking, the high-resolution imaging can be achieved. A combined motion compensation method is proposed to produce well-focused and coherent inverse synthetic aperture radar (ISAR) images from different channels, based on which the interferometric imaging is performed, thus forming the 3D imaging results. Lastly, proof-of-principle experiments were performed with a 0.2 THz SIMO imaging system, verifying the effectiveness of the proposed scheme. Non-cooperative moving targets were accurately tracked and the 3D images obtained clearly identify the targets. Moreover, the dynamic imaging results of the moving targets were achieved. The promising results demonstrate the superiority of the proposed scheme over the existing THz imaging systems in realizing 3D imaging for moving targets. The proposed scheme shows great potential in detecting and monitoring moving targets with non-cooperative movement, including unmanned military vehicles and space debris.