Bistatic System Research Articles

Range-Domain Subspace Detector in the Presence of Direct Blast for Forward Scattering Detection in Shallow-Water Environments

This paper aims to detect a target that crosses the baseline connecting the source and the receiver in shallow-water environments, which is a special scenario for a bistatic sonar system. In such a detection scenario, an intense sound wave, known as the direct blast, propagates directly from the source to the receiver without target scattering. This direct blast usually arrives at the receiver simultaneously with the forward scattering signal and exhibits a larger intensity than the signal, posing a significant challenge for target detection. In this paper, a range-domain subspace is constructed by the horizontal distance between the source/target and each element of a horizontal linear array (HLA) when the ranges of environmental parameters are known a priori. Meanwhile, a range-domain subspace detector based on direct blast suppression (RSD-DS) is proposed for forward scattering detection. The source and the target are located at different positions, so the direct blast and the scattered signal are in different range-domain subspaces. By projecting the received data onto the orthogonal complement subspace of the direct blast subspace, the direct blast can be suppressed and the signal that lies outside the direct blast subspace is used for target detection. The simulation results indicate that the proposed RSD-DS exhibits a performance close to the generalized likelihood ratio detector (GLRD) while requiring less prior knowledge of environments (only known are the ranges of the sediment sound speed and the bottom sound speed), and its robustness to environmental uncertainties is better than that of the latter. Moreover, the proposed RSD-DS exhibits better immunity against the direct blast than the GLRD, since it can still work effectively at a signal-to-direct blast ratio (SDR) of −30 dB, while the GLRD stops working in this case.

The Design of a Long-Distance Signal Transmission System in a Bistatic RCS Testing System

In wireless communication and radar systems, long-distance signal transmission poses significant challenges that affect overall system performance. In this paper, we propose an innovative bistatic radar cross section (RCS) testing system designed to address these challenges, with a particular focus on its long-distance signal transmission capabilities. This system is capable of accurately measuring the RCS of a target and improving multipath channel modeling accuracy by using precise RCS values. The integrated upper computer software extracts amplitude and phase information from received echo signals, processes these data, and provides detailed outputs including the target’s RCS, one-dimensional image, and two-dimensional image. The experimental results obtained prove that this system can not only achieve effective long-distance signal transmission but also substantially enhance the accuracy of RCS measurements, offering reliable support for multipath channel modeling. However, the conclusions drawn are preliminary and require further experimental validation to fully substantiate this system’s performance. Future work will focus on improving system accuracy, minimizing the impact of environmental noise, and optimizing data-processing methods to enhance the efficiency of wireless communication and radar applications.

A Novel SAR Imaging Method for GEO Satellite–Ground Bistatic SAR System with Severe Azimuth Spectrum Aliasing and 2-D Spatial Variability

The satellite–ground bistatic configuration, which uses geosynchronous synthetic aperture radar (GEO SAR) for illumination and ground equipment for reception, can achieve wide coverage, high revisit, and continuous illumination of interest areas. Based on the analysis of the signal characteristics of GEO satellite–ground bistatic SAR (GEO SG-BiSAR), it is found that the bistatic echo signal has problems of azimuth spectrum aliasing and 2-D spatial variability. Therefore, to overcome those problems, a novel SAR imaging method for a GEO SG-BiSAR system with severe azimuth spectrum aliasing and 2-D spatial variability is proposed. Firstly, based on the geometric configuration of the GEO SG-BiSAR system, the time-domain and frequency-domain expressions of the signal are derived in detail. Secondly, in order to avoid the increasing cost caused by traditional multi-channel reception technology and the processing burden caused by inter-channel errors, the azimuth deramping is executed to solve the azimuth spectrum aliasing of the signal under the special geometric structure of GEO SG-BiSAR. Thirdly, based on the investigation of azimuth and range spatial variability characteristics of GEO SG-BiSAR in the Range Doppler (RD) domain, the azimuth spatial variability correction strategy is proposed. The signal corrected by the correction strategy has the same migration characteristics as monostatic radar. Therefore, the traditional chirp scaling function (CSF) is also modified to solve the range spatial variability of the signal. Finally, the two-dimensional spectrum of GEO SG-BiSAR with modified chirp scaling processing is derived, followed by the SPECAN operation to obtain the focused SAR image. Furthermore, the completed flowchart is also given to display the main composed parts for GEO SG-BiSAR imaging. Both azimuth spectrum aliasing and 2-D spatial variability are taken into account in the imaging method. The simulated data and the real data obtained by the Beidou navigation satellite are used to verify the effectiveness of the proposed method.

Monitoring Mosquito Abundance: Comparing an Optical Sensor with a Trapping Method.

Optical sensors have shown significant promise in offering additional data to track insect populations. This article presents a comparative study between abundance measurements obtained from a novel near-infrared optical sensor and physical traps. The optical instrument, named an Entomological Bistatic Optical Sensor System, or eBoss, is a non-destructive sensor operating in the near-infrared spectral range and designed to continuously monitor the population of flying insects. The research compares the mosquito aerial density (#/m3) obtained through the eBoss with trap counts from eight physical traps during an eight-month field study. The eBoss recorded over 302,000 insect sightings and assessed the aerial density of all airborne insects as well as male and female mosquitoes specifically with a resolution of one minute. This capability allows for monitoring population trends throughout the season as well as daily activity peaks. The results affirmed the correlation between the two methods. While optical instruments do not match traps in terms of taxonomic accuracy, the eBoss offered greater temporal resolution (one minute versus roughly three days) and statistical significance owing to its much larger sample size. These outcomes further indicate that entomological optical sensors can provide valuable complementary data to more common methods to monitor flying insect populations, such as mosquitoes or pollinators.

Imaging Seafloor Features Using Multipath Arrival Structures

In this paper, we propose an imaging method for seafloor features based on multipath arrival structures. The bistatic sonar system employed consists of a vertical transmitting array and a horizontal towed array. The conventional back projection (BP) method, which considers the direct path from the source to the seafloor scatterer and then to the receiver, is used in this system. However, discrepancies between the calculated delay values and the actual propagation delay result in projection deviations and offsets in the seafloor features within sound intensity images. To address this issue, we analyze the multipath structures from the source to the scatterer and then to the receiver based on ray theory. The delay at each grid is calculated using different multipaths, considering the distances from the seafloor grids to the source and the receiver. In the direct zone, the delay is determined using the direct ray and the surface reflection ray, while in the bottom bounce area, the delay is calculated using the bottom–surface reflection ray and the surface–bottom–surface reflection ray. Numerical simulations and experimental results demonstrate that the proposed method rectifies the delay calculation issues inherent in the conventional method. This adjustment enhances the accuracy of the projection, thereby improving the imaging quality of seafloor features.

Spreading Sea Clutter Suppression for High-Frequency Hybrid Sky-Surface Wave Radar Using Orthogonal Projection in Spatial–Temporal Domain

In recent years, the high-frequency hybrid sky-surface wave radar (HSSWR) has been increasingly used for target detection applications. Nevertheless, the specific bistatic system layout and the phase path disturbances induced by the ionospheric propagation channel may severely spread the sea clutter spectrum, thereby deteriorating the detection ability of the HSSWR for slow-moving targets. In this work, a novel subspace method based on the hybrid use of the amplitude and phase estimator (APES) and the orthogonal projection (OP) in the spatial–temporal domain, denoted as the APES-OP method, is proposed to suppress the spreading first-order sea clutter of the HSSWR. The distribution characteristics of targets and first-order sea clutter in the spatial–temporal domain were investigated, and a time-domain subspace signal model was adopted to describe targets perturbed by ionospheric phase path modulation. An APES filter was adopted to filter out the potential targets with a traversal approach to avoid attenuating desired signals while suppressing sea clutter. After that, sampling data from multi-channels and slow-time domains at the cell under test were employed to construct a spatial–temporal matrix, which was then utilized to obtain the sea clutter subspace by singular value decomposition. Simulation results indicate that the proposed algorithm can suppress sea clutter while retaining the target, even if the target is buried by sea clutter. The processing results of measured data further demonstrate the efficiency of the proposed algorithm. After sea clutter suppression, the target obscured by clutter can be revealed, and the signal-to-clutter ratio of the target is greatly improved.

Interferometric Calibration Model for the LuTan-1 Mission: Enhancing Digital Elevation Model Accuracy

The LuTan-1 (LT-1) mission, China’s first civilian bistatic spaceborne Synthetic Aperture Radar (SAR) mission, comprises two L-band SAR satellites. These satellites operate in bistatic InSAR strip map mode, maintaining a formation flight with an adjustable baseline to generate global digital elevation models (DEMs) with high accuracy and spatial resolution. This research introduces a dedicated interferometric calibration model for LT-1, tackling the unique challenges of the bistatic system, such as interferometric parameter coupling and the π-ambiguity problem caused by synchronization phase errors. This study validates the model using SAR images from LT-1 and Xinjiang corner reflector data, achieving interferometric phase accuracy better than 0.1 rad and baseline accuracy better than 2 mm, thereby producing high-precision DEMs with a height accuracy meeting the 5 m requirement.

Long-Range Imaging LiDAR with Multiple Denoising Technologies

The ability to capture and record high-resolution images over long distances is essential for a wide range of applications, including connected and autonomous vehicles, defense and security operations, as well as agriculture and mining industries. Here, we demonstrate a self-assembled bistatic long-range imaging LiDAR system. Importantly, to achieve high signal-to-noise ratio (SNR) data, we employed a comprehensive suite of denoising methods including temporal, spatial, spectral, and polarization filtering. With the aid of these denoising technologies, our system has been validated to possess the capability of imaging under various complex usage conditions. In terms of distance performance, the test results achieved ranges of over 4000 m during daylight with clear weather, 19,200 m at night, 6700 m during daylight with haze, and 2000 m during daylight with rain. Additionally, it offers an angular resolution of 0.01 mrad. These findings demonstrate the potential to offer comprehensive construction strategies and operational methodologies to individuals seeking long-range LiDAR data.

Highlight Model of Underwater Target Acoustic Scattering in the Bistatic System

<p>A highlight model for the target signal in the bistatic system is proposed based on the calculation results of the acoustic field, which means that the target can be equivalent to several highlights according to certain characteristics of the target when calculating the scattering acoustic field. The parameters of the highlight model were given based on the calculation results. Two sets of experiments on the highlight model of the bistatic system were conducted with a fixed transmitter and a target, one being the receiver rotating around the target, and the other being the target rotating around itself at different receiver angles. The comparison of multi angle pseudo color images through experiments and theoretical simulations shows that the scattering highlight model can describe the scattering characteristics of the target to a certain extent. Only geometric highlights were simulated in this article, and it can be foreseen that bistatic system highlight model of scattering waves can be accurately described with a full understanding of the target.</p> <p> </p>

Automated Estimation of Sub-Canopy Topography Combined with Single-Baseline Single-Polarization TanDEM-X InSAR and ICESat-2 Data

TanDEM-X bistatic interferometric system successfully generated a high-precision, high-resolution global digital elevation model (DEM). However, in forested areas, two core problems make it difficult to obtain sub-canopy topography: (1) the penetrability of short-wave signals is limited, and the DEM obtained in dense forest areas contains a significant forest signal, that is, the scattering phase center (SPC) height; and (2) the single-baseline and single-polarization TanDEM-X interferometric synthetic aperture radar (InSAR) data cannot provide sufficient observations to make the existing physical model reversible for estimating the real surface phase, whereas the introduction of optical data makes it difficult to ensure data synchronization and availability of cloud-free data. To overcome these problems in accurately estimating sub-canopy topography from TanDEM-X InSAR data, this study proposes a practical method of sub-canopy topography estimation based on the following innovations: (1) An orthogonal polynomial model was established using TanDEM-X interferometric coherence and slope to estimate the SPC height. Interferometric coherence records forest height and dielectric property information from an InSAR perspective and has spatiotemporal consistency with the InSAR-derived DEM. (2) Introduce Ice, Cloud, and Land Elevation Satellite-2 (ICESat-2) data to provide more observational information and automatically screen ICESat-2 control points with similar forest and slope conditions in the local area to suppress forest spatial heterogeneity. (3) A weighted least squares criterion was used to solve this model to estimate the SPC height. The results were validated at four test sites using high-precision airborne light detection and ranging (LiDAR) data as a reference. Compared to the InSAR-derived DEM, the accuracy of the sub-canopy topography was improved by nearly 60%, on average. Furthermore, we investigated the necessity of local modeling, confirming the potential of the proposed method for estimating sub-canopy topography by relying only on TanDEM-X and ICESat-2 data.

Analysis of bistatic multiphoton quantum radar cross section for the cylindrical surface.

A closed-form model of bistatic multiphoton quantum radar cross section (QRCS) for the cylindrical surface, the main structure of typical aircraft, especially missiles, is established to analyze the system and scattering characteristics. The influence of curvature of the three-dimensional target on QRCS is analyzed. By comparing and analyzing the bistatic multiphoton QRCS for a cylinder and a rectangular plate, we find that the QRCS for the convex surface target is the extension of the QRCS for the planar target with inhomogeneous atomic arrangement intervals and patterns. The characteristics of cylindrical QRCS are discussed by combining the transceiver system and the photon number of the transmitted signal, and the influences of the cylindrical radius, cylindrical length, and incident photon number on QRCS are analyzed. The bistatic results provide guidance on potential strong scattering directions for the target under various directions of photon incidence. Compared with the plane target, the cylindrical target amplifies scattering intensity near the target surface at the scattering angle side in the bistatic system. A bistatic multiphoton quantum radar system can achieve sharpening and amplification of the main lobe of the QRCS for a cylinder in an extensive scattering angle range. Bistatic multiphoton quantum radar has better visibility for the cylinder with a smaller length. These characteristics will provide prior information for research in many fields, such as photonic technology, radar technology, and precision metrology.

Terahertz bistatic three-dimensional computational imaging of hidden objects through random media

Random media pose limitations on the imaging capability of photoelectric detection devices. Currently, imaging techniques employed through random media primarily operate within the laser wavelength range, leaving the imaging potential of terahertz waves unexplored. In this study, we present an approach for terahertz bistatic three-dimensional imaging (TBTCI) of hidden objects through random media. By deducing the field distribution of bistatic terahertz time-domain spectroscopy system, and proposing an explicit point spread function of the random media, we conducted three-dimensional imaging of hidden objects obscured by the random media. Our proposed method exhibits promising applications in imaging scenarios with millimeter-wave radar, including non-invasive testing and biological imaging.

Метод оценки отражательных характеристик объектов с использованием мульти-фокусных зондирующих полей

The possibility of using special (multi-focus) fields for solving problems related to the experimental evaluation of the reflective characteristics of objects has been considered. Regularities in the formation of multi-focus fields in three-dimensional cases have been specified. The advantages of the proposed approach have been demonstrated in terms of flexible adjustment of signal processing algorithms based on digital signal filtering methods. This ensures obtaining radar characteristics of objects, considering the peculiarities of constructing single-position, bi-static, and multi-static location systems.

Active Control of the Reflection Coefficient of an Underwater Surface

From a strategic point of view, it is essential to protect underwater vehicles from being detected by opponents. Modern mono- or bistatic sonar systems are capable of precisely determining the position of a watercraft. In order to triangulate the positions of watercrafts, the sonar sends out acoustic signals that are reflected by the vehicles’ surfaces. These deflected sound waves are subsequently detected and evaluated. How well an object can be detected using a sonar can be measured according to the target strength. Through their shape, construction and choice of materials, modern underwater vehicles are optimized in relation to minimizing their radiated and reflected sound waves; in this way, their target strength is minimized. These passive measures are particularly effective in the medium- and high-frequency range down to 1500 Hz. To effectively reduce reflections at lower frequencies, an active system is developed in this study and evaluated in a laboratory test with a water-filled impedance tube. The incident sound waves were measured in front of an active surface and then processed using an adaptive control system based on an FPGA platform. The system operates with a very thin piezoceramic applied to the surface in order to minimize the reflections of the sound waves. The laboratory tests showed the high effectiveness of the system under the influence of sonar-like signals.

First Assessment of Bistatic Geometric Calibration and Geolocation Accuracy of Innovative Spaceborne Synthetic Aperture Radar LuTan-1

LuTan-1 (LT-1) is a bistatic synthetic aperture radar (BiSAR) system consisting of two identical L-band SAR satellites. The bistatic mode of LT-1 plays a critical role in generating high-precision digital elevation models (DEMs), which requires precise geometric calibration of initial range and azimuth times for both SARs to ensure the reliability and quality of geolocation. However, existing geometric calibration methods predominantly focus on monostatic SAR systems, with limited literature on slave SAR calibration in bistatic systems. This research addresses this gap by establishing geometric calibration models for both SARs based on signal echo history and the range–Doppler model. The geometric errors are effectively resolved using corner reflector data from Xinjiang, China. Through statistical analysis of LT-1 SAR images acquired between July and November in bistatic mode, this paper has demonstrated range delay accuracy of better than 5 ns and azimuth time accuracy of better than 0.1 ms. This level of precision translates into a positional accuracy better than 0.8 m. The proposed models have been successfully applied to geometric calibration, providing precise geolocation for LT-1, thus enhancing its utility for a wide range of Earth observation applications. This paper is the first endeavor to present the assessment of the geometric calibration and geolocation accuracy of LT-1 and discuss the results of the bistatic geometric calibration of the master and slave SARs in a BiSAR formation.

Improving the Pulse-Limited Footprint Resolution of GNSS-R Based on the Novel Joint Bandwidth Method

The bistatic global navigation satellite system’s (GNSS) signal reflection technology has become an effective means of space-based sea surface wind field retrieval and height retrieval. By adopting a wider signal bandwidth, a higher pulse-limited footprint resolution can be achieved. However, for the GNSS-Reflectometry (GNSS-R) system, its signal bandwidth is affected by the signal bandwidth of the GNSS satellite, which limits the further improvement of the pulse-limited footprint resolution. This article proposes a method based on the novel signal bandwidth joint principle to improve the resolution of GNSS-R pulse-limited footprints. Firstly, currently in-orbit GNSS-R satellites use the traditional single frequency band (TSFB) method, which is limited by the GNSS satellite’s signals and has a theoretical upper limit on its signal bandwidth. In response to this issue, this article proposes the novel joint bandwidth (NJBW) method (Galileo E5a and E5b signals) based on the auto-correlation function (ACF) signal ambiguity theory. The NJBW method reduces the main lobe width of the ACF of the GNSS-R signal by jointly processing the signals of E5a and E5b frequency bands, thus improving the pulse-limit footprint resolution of GNSS-R. Secondly, in order to verify the improvement effect of the novel joint bandwidth method on the pulse-limited footprint resolution of GNSS-R, this paper designs and fabricates an NJBW antenna verification prototype for the joint Galileo E5a and E5b frequency band and tests it in a microwave anechoic chamber. The test results indicate that the radio frequency (RF) bandwidth of the NJBW antenna validation prototype can cover both the frequency bands of E5a and E5b, making it suitable for use as the NJBW method for the GNSS-R receiving antenna. The bandwidth test values of the NJBW antenna validation prototype are consistent with the design values, which verifies the correctness of the NJBW antenna design model and further proves the feasibility of the NJBW method. Thirdly, based on the joint Galileo E5a and E5b frequency band signals, the NJBW method was applied to analyze the improvement effect of the pulse-limited footprint resolution. Compared to the TSFB method, the application of the NJBW method can increase the resolution of the GNSS-R pulse-limiting footprint by 1.73 times, which effectively improves the performance of the GNSS-R system. The NJBW method proposed in this article provides the theoretical method foundation and key technical support for sea surface wind field retrieval and height retrieval and the antenna design for the future high-precision and high pulse-limited footprint resolution GNSS-R sea surface wind field retrieval and height retrieval verification satellite.

Moon Imaging Performance of FAST Radio Telescope in Bistatic Configuration with Other Radars

Ground-based radar has been used for Moon imaging for more than 60 years. Five-hundred-meter Aperture Spherical radio Telescope (FAST), as the largest radio telescope on Earth, holds significant potential for celestial imaging missions with its exceptional sensitivity. A bistatic Synthetic Aperture Radar (SAR) Moon imaging model that incorporates FAST and other transmitting radars is presented. The objective of this paper is to design the imaging parameters of this bistatic configuration based on the required resolution, and to estimate the resolution performance based on a given bistatic system capability. Considering the ultra-far range and the ultra-long observation time between the radars and the Moon, the geometric relationship involved in this bistatic configuration is significantly distinct from the bistatic configuration of airborne and spaceborne radars. Therefore, this paper accurately derives the two-dimensional resolution on the Moon’s surface. First of all, the models of the Earth’s surface and the Moon’s surface, and the celestial motion of the Earth and Moon are established using WGS-84 and JPL-DE421, given by STK. Secondly, the bistatic range history within the observation time is calculated in terms of continuous celestial motion instead of the popular ‘stop-and-go’ assumption. Thirdly, no approximation is used in the resolution derivation process, and, in addition to the two-dimensional resolutions, the incident angle and the included angle are also given to describe the imaging performance. This method can also be extended to other bistatic-station and single-station celestial imaging, providing support for radar location and parameters design, for observation time span selection, for observation area selection, and for imaging performance estimation. The echo generation and imaging for point targets set on the Moon are shown. The simulation results prove the validity and accuracy of the proposed method in the paper.

Algorithm for the Joint Flight of Two Uncrewed Aerial Vehicles Constituting a Bistatic Radar System for the Soil Remote Sensing

The study of soil agrophysical and agrochemical properties is based on ground-based point measurements and measurements conducted using radiometric remote sensing systems (satellite or airborne). A disadvantage of the existing remote sensing systems using normal surface irradiation is the insignificant depth of penetration of the probing radiation into the soil layer. It is proposed to use a radar system for remote sensing agricultural lands to eliminate this drawback. The system uses a method for assessing the soil’s physical and chemical properties based on the interference measurements of direct and reflected electromagnetic waves at incidence angles that provide a total refraction effect, i.e., close to Brewster’s angle. The possibility of using this method for remote assessment of soil’s physical and chemical properties, including the subsurface layer moisture, was established. A feature of the bistatic system is that it is necessary to coordinate the mutual arrangement of the transmitting and receiving positions, which imposes special requirements on the UAVs’ flight algorithm. The UAVs’ relative position makes it possible to form the conditions for the manifestation of the total refraction effect, to determine the current value of Brewster’s angle, and to fix these conditions for the subsequent flight, making it possible to measure the soil’s physical and chemical parameters. The research results can be used to implement precision farming technology in hard-to-reach places, large agricultural areas, and digital agriculture.

MVFRnet: A Novel High-Accuracy Network for ISAR Air-Target Recognition via Multi-View Fusion

Inverse Synthetic Aperture Radar (ISAR) is a promising technique for air target imaging and recognition. However, the traditional monostatic ISAR only can provide partial features of the observed target, which is a challenge for high-accuracy recognition. In this paper, to improve the recognition accuracy of air targets, we propose a novel recognition network based on multi-view ISAR imaging and fusion, called Multi-View Fusion Recognition network (MVFRnet). The main structure of MVFRnet consists of two components, the image fusion module and the target recognition module. The fusion module is used for multi-view ISAR data and image preprocessing and mainly performs imaging spatial match, image registration, and weighted fusion. The recognition network consists of the Skip Connect Unit and the Gated Channel Transformation (GCT) attention module, where the Skip Connect Unit ensures the extraction of global depth features of the image and the attention module enhances the perception of shallow contour features of the image. In addition, MVFRnet has a strong perception of image details and suppresses the effect of noise. Finally, simulated and real data are used to verify the effectiveness of the proposed scheme. Multi-view ISAR echoes of six types of aircraft are produced by electromagnetic simulation software. In addition, we also build a millimeter wave ground-based bistatic ISAR experiment system and collect multi-view data from an aircraft model. The simulation and experiment results demonstrate that the proposed scheme can obtain a higher recognition accuracy compared to other state-of-the-art methods. The recognition accuracy can be improved by approximately 30% compared with traditional monostatic recognition.

Baseline Calibration of L-Band Spaceborne Bistatic SAR TwinSAR-L for DEM Generation

The Terrain Wide-swath Interferometric L-band Synthetic Aperture Radar (TwinSAR-L) mission is a spaceborne bistatic synthetic aperture radar (SAR) mission to derive a high-quality global digital elevation model (DEM). The prerequisite of the high-accuracy DEM is knowing the interferometric baseline with high precision. The challenging problem is that the baseline of the bistatic system is highly dynamic due to the fast relative motion between the two satellites. In this paper, a pixel-related baseline model based on the geometrical shift is proposed to accurately reflect the position change of satellites. The baseline error is then calibrated using height gradient information and a small number of point targets with a slight incidence angle difference, eliminating the need for low-frequency corner reflectors and avoiding the difficulty of selecting a calibration site. The proposed method has been successfully exploited during the initial Commissioning Phase of TwinSAR-L, demonstrating its effectiveness in evaluating the precise baseline and supporting the generation of high-precision DEM.