Next Article in Journal
Ecosystem Resistance and Resilience after Dry and Wet Events across Central Asia Based on Remote Sensing Data
Next Article in Special Issue
High-Resolution Azimuth Missing Data SAR Imaging Based on Sparse Representation Autofocusing
Previous Article in Journal
Distribution of Grazing Paths and Their Influence on Mountain Vegetation in the Traditional Grazing Area of the Tien-Shan Mountains
Previous Article in Special Issue
The Calibration Method of Multi-Channel Spatially Varying Amplitude-Phase Inconsistency Errors in Airborne Array TomoSAR
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Static High Target-Induced False Alarm Suppression in Circular Synthetic Aperture Radar Moving Target Detection Based on Trajectory Features

1
Radar Monitoring Technology Laboratory, School of Information Science and Technology, North China University of Technology, Beijing 100144, China
2
School of Electronic and Information Engineering, Beihang University, Beijing 100191, China
*
Author to whom correspondence should be addressed.
Remote Sens. 2023, 15(12), 3164; https://doi.org/10.3390/rs15123164
Submission received: 25 May 2023 / Revised: 12 June 2023 / Accepted: 16 June 2023 / Published: 18 June 2023

Abstract

:
The new mode of Circular Synthetic Aperture Radar (CSAR) has several advantages including multi-aspect and long-time observation, which can generate high-frame-rate image sequences to detect moving targets with a single-channel system. Nonetheless, due to CSAR being sensitive to 3D structures, static high targets are observed in scene display rotational motion within CSAR subaperture image sequences. Such motion can cause false alarms rising when utilizing image sequence-based moving target detection methods like logarithm background subtraction (LBS). To address this issue, this paper first thoroughly analyzes the moving target and static high target’s difference for the trajectory in an image sequence. Two new trajectory features of the rotation angle and moving distance are proposed to differentiate them. Based on the features, a new false alarm suppression method is proposed. The method first utilizes LBS to obtain coarse binary detection results comprising both moving and static high targets, then employs morphological filtering to eliminate noise. Next, DBSCAN and target tracking steps are employed to extract the trajectory features of the target and false alarm. Finally, false alarms are suppressed with trajectory-based feature discriminators to output detection results. The W-band CSAR open dataset is used to validate the proposed method’s effectiveness.

1. Introduction

Circular Synthetic Aperture Radar (CSAR) is a new SAR imaging mode capable of continuously observing areas of interest with a circular track [1,2]. Its long-time and multi-aspect capabilities are conducive to target detection and tracking, as well as estimating the velocity of moving targets [3,4]. Recently, such advantages for moving target detection have become a hot topic due to their potential to improve the detection performance [5,6,7].
Although the multi-channel system and corresponding methods are a common scheme in detecting moving targets (such as space–time adaptive processing) [8,9,10], it is important to develop single channel-based detection methods. The single-channel SAR still plays an important role in current SAR application. Several recent studies have shown that utilizing the advantages of CSAR can achieve a good performance in detecting moving targets, even with a single-channel system [11,12,13,14]. One topic of interest is the use of image sequences to detect moving targets in single-channel CSAR [15,16]. This kind of research can be categorized into two classes: detection with the moving target’s shadow [17,18,19] and detection with the target’s smear/displaced signature in the image sequence [20,21,22]. The former method relies on detecting the moving shadow to locate the target. However, this approach requires a high signal-to-noise ratio and only works in a higher frequency range like W-band SAR [23]. On the other hand, the latter method detects the target’s smear/displaced signature in image sequences, which has fewer constraints and is therefore more applicable. One of the typical image sequence-based methods is logarithm background subtraction (LBS) [24,25], which uses a median filter to obtain the static background image, and then extracts the moving target signature by subtracting it from the original images. Another typical example is low-rank-based methods like the Robust Principal Component Analysis (RPCA) and Go Decomposition (GoDec) [26,27,28]. These methods can separate the background and foreground images by utilizing their low-rank (background) and sparsity (foreground) nature in data.
Further investigation reveals that CSAR’s curved geometry makes static high targets such as buildings and towers have a rotational motion in subaperture image sequences. This happens due to the projections of these targets onto different positions relative to the corresponding azimuth viewing angle, ultimately leading to false alarms raising in the detection. Figure 1 shows the rotation of a highly stationary target in the image sequence. The four images are images of different azimuths in the image sequence at different times.
At present, there is no research on the suppression of a high-target false alarm in circular SAR (to the best of our knowledge), but according to the three-dimensional terrain extraction of circular SAR, there have been many studies on the use of circular SAR to conduct three-dimensional terrain [29,30]. In 2017, a team at Nanjing University of Posts and Telecommunications conducted a circular SAR DEM extraction of the scene and studied the extraction method [31]. In summary, the above research is only used to extract terrain or the target height. Although it is not combined with moving target research, its signal model has laid a foundation for moving target false alarm suppression. Since the existing moving target detection algorithm cannot solve the false alarm problem caused by the defocusing trajectory of the stationary elevation target, it is necessary to distinguish the moving target from the highly stationary target. The false alarm suppression algorithm for the false alarm problem caused by the stationary elevation target is studied to improve the detection accuracy and detection efficiency of the moving target in the observation scene using the CSAR subaperture image sequence. The high-precision CSAR moving target detection results are helpful for further research on target recognition and tracking applications.
A viable way to address this issue is to utilize the trajectory shape difference of the moving target and false alarm (static high target). Hence, this paper first thoroughly analyzes the moving target and static high target’s difference for the trajectory in an image sequence. Two new trajectory features of the rotation angle and moving distance are proposed to differentiate them. Based on the features, a new false alarm suppression method is proposed. The method first utilizes LBS to obtain coarse binary detection results comprising both moving and static high targets, then employs morphological filtering to eliminate noise. Next, DBSCAN and target tracking steps are employed to extract the trajectory features of the target and false alarm. Finally, false alarms are suppressed with trajectory-based feature discriminators to output detection results.
The rest of this paper is organized as follows. Section 2 is for the dataset. Section 3 is for the signal model. It analyzes the rotation motion of the static high target and makes a comparison with the moving target’s signature, which provide the basis for using the trajectory feature to suppress the false alarm. Section 4 is for the proposed false alarm suppression algorithm. Section 5 is for the real data experiment; a comparison experiment with other image sequence-based moving target detection methods is also described. Section 6 is the discussion and the last section is the conclusion.

2. Materials and Dataset

The proposed false alarm suppression algorithm was evaluated using open W-band CSAR data [23], and the parameters are provided in Table 1. A specific area of interest in the data was selected for processing purposes, measuring 300 × 200 pixels and labeled with a red box, as illustrated in Figure 2. This region includes both a moving target, represented by a green circle, and a static high target, designated in a yellow circle. The outcomes of each step throughout the processing procedure are elaborated in detail in this section.
Figure 3 illustrates the rotation phenomenon of billboards; Figure 3a displays its optical image while Figure 3b presents its rotations in 4 sub-images. For instance, the detection results obtained using the logarithm background subtraction (LBS) method [24,25] are shown in Figure 4, where the green box indicates the moving target, and the red circle represents the billboard. Specifically, Figure 4a shows the original images and Figure 4b shows the corresponding detection result, which identifies that the billboard remains after applying the LBS method. Therefore, it is crucial to suppress the false alarms caused by static high targets to achieve a good detection performance.

3. Signal Model of Static High Target’s Rotation Motion

In this section, the signal model of the static high target’s rotation motion is described and analyzed with point target simulation. Specifically, the first subsection involves derivation of a formula for the rotation radius of an ideal static point target at a fixed height. In the second subsection, we present a theoretical analysis of the rotation radius formula, then investigate the differences between the trajectories of static high targets and moving targets. The comparison provides the basis for using a trajectory feature to suppress a false alarm.

3.1. Defocus Radius of Static High Target

Figure 5 details the geometric configuration of CSAR, where (a) represents the three-dimensional model and (b) provides a frontal view of the imaging plane A P P B B . An approximate circular path is traversed by the radar, which is used to fit the real application; A denotes the position of the radar in this setup, while P (x0, y0, Δ h ) signifies the coordinate location of the ideal point target.
As depicted in Figure 3b, the distance between the radar and the target is represented as R, while H signifies the vertical height distance to the target. In instances where the imaging plane is established at the same elevation as that of the target Δ h , the object under observation is imaged at its actual position. The horizontal distance from the imaging plane to the target is considered as the ground range
R g = R 2 H 2
When the imaging plane selects the ground plane, the height difference is Δ h . Then, the target is projected on a ground plane based on the iso-range principle. Its ground range is
In situations where the imaging plane is aligned with the ground plane; a height difference (represented as “ Δ h ”) between the target and the imaging plane exists. In such instances, the target is projected onto the ground plane using the iso-range principle, representing the horizontal distance from the radar to the target on the earth’s surface. The corresponding quantity is effectively the ground range of the target
R g = R 2 ( H + Δ h ) 2
The ground range difference P′C represents the defocus radius, which is given as
Δ R g = R g R g = Δ h ( 2 H + Δ h ) R 2 H 2 + R 2 ( H + Δ h ) 2
Equation (3) is the exact formula to calculate the range difference. To give the reader a simple physical meaning, here, we further use an approximation in deduction. In the case of a CSAR system, it is customary for the height of the radar to exceed the elevation of the target. As such, Formula (3) may be simplified and expressed as follows:
Δ R g H Δ h R 2 H 2 = Δ h tan ( φ )
In Equation (4), it is found out that the distance between the target and its projection is only dependent on target elevation and the instantaneous incidence angle. It does not require the information of the radar’s three-dimensional position. This provides the basis for the analysis on designing and extracting the features for discrimination. Under an ideal circular case, for targets located at the scene center, the incidence angle remains constant, leading to a defocused circular signal within the full aperture image. In contrast, when the target is situated off-center, the resulting defocused signal is in an approximated circular shape. For both cases, we can see the rotation. In the next subsection, we delve into simulation-based investigations of this phenomenon and compare it with the trajectory observed in moving targets’ scenarios.

3.2. Trajectory Features of Moving Target and Static High Target

For this section, simulation was conducted to confirm the efficacy of the rotation radius formula. Since the simulation aimed to validate the signal model and phenomenon, here, we used X-band SAR parameters for simulation for convenience. Specifically, Table 2 illustrates the parameters, with the utilization of two target points located at 0 m and 5 m from the scene center. This study aimed to validate the accuracy of the aforementioned formula in an ideal circular track and constant-speed moving target environment.
The simulation results are presented in Figure 6, where Figure 6a depicts the subaperture image of a 45° azimuth angle, while Figure 6b presents a combination of all 24 subaperture images. In Figure 6a, it can be observed that target A is not offset, whereas target B is projected. By extracting the distance between targets A and B in Figure 6a, the rotation radius was estimated to be 77.85 m, which is illustrated as the red dashed line. The theoretical distance calculated with Equation (4) was 77.6 m. The minor discrepancy between these two values could be attributed to maximum selection during extraction, thereby validating the effectiveness of Equation (4). Figure 6b further corroborates the simulation analysis, depicting a focused target A at the scene center and a circular trajectory, indicating the phenomenon of rotation exhibited by target B.
Subsequently, a comparison was drawn between the trajectory features of a moving target and those of a stationary high target. As discussed in [32], the shape of a constantly moving target differs notably from the near-circle shape analyzed in Equation (4). To facilitate this comparison, we present the simulation results in Figure 6.
The simulation parameters mentioned in Table 1 were implemented, with the moving target placed within the range and azimuth dimensions and set to move at a constant speed of 0.2 m/s in simulation. Figure 7 shows the resulting trajectory feature comparison between the moving target and a stationary high target, with subaperture images covering an azimuth aperture interval of 15° from 0° to 360° presented in Figure 7a, which is the height target of non-origin and a complete image combination for the static high target shown in Figure 7b. The trajectory of the static high target, as evidenced by its circular path of rotation around the target location, contrasts significantly with the ‘W’- or ‘V’-shaped signatures generated by targets moving along the X-label and Y-label, which are depicted in Figure 7c,d, respectively. These observations provide the theoretical foundations necessary for distinguishing false alarms.

4. False Alarm Suppression Algorithm

The logarithm background subtraction was adopted to obtain the coarse detection result. As illustrated in Figure 2, the detection outcomes comprise of moving targets, false alarms, and residual clutters. To discriminate the moving target and suppress a false alarm, we developed a three-stage algorithm to address these challenges. Firstly, the use of morphological processing facilitated the removal of residual clutter. The second stage involved discrimination between false alarms and moving targets via the employment of a trajectory feature, which utilized clustering and tracking methods for feature extraction. Lastly, false alarm removal was achieved through the application of trajectory feature-based discrimination, with the resulting detection outcomes being reported as output.
A flowchart of the proposed approach is depicted in Figure 8, featuring the following steps: (1) inputting an image sequence of CSAR data; (2) processing with the logarithm background subtraction method; (3) morphological processing; (4) DBSCAN clustering implementation; (5) target tracking; (6) trajectory feature extraction; (7) discrimination based on trajectory features; and (8) false alarm removal and output processing results.
Next, the key step is introduced in detail.

4.1. Logarithm Background Subtraction

The logarithm background subtraction (LBS) method represents a single-channel CSAR moving target detection algorithm [24,25]. This algorithm extracts background (i.e., static scene) images from the sequence, followed by subsequent background subtraction procedures to facilitate detection of the moving target-containing foreground image. Due to this being the preprocessing step in this paper, details are not given here. For other image sequence-based detection methods like the low-rank method, they can be used as replacement. However, these methods also face the false alarm problem due to their sensitivity to the fast changes in an image sequence. During this first-stage detection, coarse detection comprises the signatures of residual clutter, static high targets, and moving targets.
The input CSAR image sequence has N images, which are denoted as I 0 ( n ) ,   n = 1 N . The coarse detection result I 1 ( n ) is given as
I 1 ( n ) = L B S [ I 0 ( n ) ]

4.2. Morphological Processing

In the second step, morphological processing [30] was employed to eliminate residual clutter and alleviate target signal breakage issues due to subtraction. The resultant detection outcomes are binary images in Figure 9, rendering the utilization of dilation and erosion-based morphological processing approaches ideal for enhancing image quality. Furthermore, false alarms are labeled with a red circle, and moving targets are labeled using green boxes while yellow circles indicate residual clutter.
It is observed that residual clutter occupies a relatively small number of pixels, rendering them susceptible to elimination through erosion techniques. For those small sizes of the static high target such as antennas, they may be suppressed in this step. For a larger-size false alarm like a tower, the rest of the steps are necessary. The erosion process involves employing a template, typically a box, where the size of the template serves as the primary parameter. The success of the algorithm is heavily dependent on achieving an appropriate template size. If the template size is too small, residual clutter remnants may remain while an excessive size could result in loss of the target signal in I 1 ( n ) . In practice, B e is the erosion template that is applied to binary images to obtain satisfactory results, whereas the outcome of this process is presented as the erosion results I 2 ( n ) :
I 2 ( n ) = erosion { I 1 ( n ) , B e }
Subsequent to the erosion process, a technique known as dilation is utilized in order to selectively repair broken sections of both target and false alarm signals, which are labelled with green boxes and red circles, respectively, in Figure 8. The effectiveness of the dilation method is highly dependent on selecting an appropriate template size as was previously observed with the erosion process. Here, the dilation template utilized in this procedure is referred to. The resulting output from this method is given, serving as a vital component for further processing and analyses, as I 3 ( n ) :
I 3 ( n ) = dilation { I 2 ( n ) , B d }

4.3. DBSCAN Clustering

Once the majority of residual clutter is successfully eliminated, the subsequent procedure involves clustering the moving target and false alarm signals. This process serves as a vital prerequisite for feature extraction. Through cluster analysis techniques [25], individual pixels can be grouped together, thereby rendering the shape, angle, center-of-mass, and other associated attributes easier to discern. Furthermore, the use of clustering methods can also offer potential solutions to issues relating to signal breakage, as broken segments may potentially be assigned membership to specific clusters.
The DBSCAN requires the consideration of two essential parameters, namely Eps and MinPts. Specifically, Eps is the radius of a given pixel’s neighborhood, whereas MinPts refers to the density threshold, the minimum required number of pixels that make up a cluster. By determining and setting these parameters at appropriate values, one can effectively cluster the pixels within the image. Ultimately, the output of DBSCAN results in the set of clustered data points I 4 ( n ) :
I 4 ( n ) = DBSCAN   [ I 3 ( n ) ]   ,   n = 1 , 2 , N
Although the targets and false alarm are clustered, the task of associating the same target in consecutive frames is not addressed. As such, an additional tracking step is necessary to ensure that this task is completed.

4.4. Target Tracking

The current section involves the implementation of a target tracking process for clusters with the aim of connecting the same target in the image sequence. It is important to note that the input image sequence makes use of overlapped subaperture images, meaning that the same target will inevitably overlap across adjacent images. The signal-to-noise ratio is greater in airborne conditions as compared to spaceborne scenarios. Hence, a straightforward approach such as utilizing Intersection Over Union (IoU) of the target in adjacent images can be employed for the purpose of tracking. The IoU ratio of the same target in two adjacent images is:
IoU = B i ( n ) B j ( n + 1 ) B i ( n ) B j ( n + 1 ) ,   n = 1 , 2 , , N
It can be expressed that B i ( n ) and B j ( n + 1 ) represent the ‘i-th’ and ‘j-th’ clusters in the ‘n-th’ and ‘(n+1)th’ images, respectively. The numerator in the given formula represents the total number of pixels that overlap between these two clusters, while the denominator refers to the sum of all the pixels in both clusters. By applying this equation to all the clusters in adjacent images, their corresponding IoU value can be obtained. Additionally, if two clusters belong to the same target and have a relatively high IoU value, they can be categorized as being assigned to an index. In practice, a suitable threshold for connecting the same target across an entire image sequence can be established using the aforementioned equation, and by setting the appropriate associated IoU-value C 0 threshold via the formula previously mentioned.
{ I area C 0 ,   same   target I area C 0 ,   not   same   target
The presented approach treats all clusters within the initial frame as potential targets for tracking purposes. The IoU tracker is subsequently applied to associate clusters in adjacent frames and successfully link such targets across multiple frames. If a target’s signature enters the image from the boundary, it is treated as a new potential target for tracking purposes; conversely, if a target moves out of the image boundary, the tracking will stop.
The next step was to extract the trajectory features and distinguish the moving target from the false alarm.

4.5. Trajectory Feature Extraction

Following the target tracking procedure, identification and labeling of the false alarm and moving targets were performed. This stage relies on the utilization of rotation angles and moving distances for differentiation. The above-mentioned two trajectory features were established through the incorporation of four fundamental features, which necessitated their preliminary extraction.

4.5.1. Basic Features

Table 3 presents the four fundamental features, which are visually depicted in Figure 10. Specifically, Figure 10a is the equivalent ellipse denoted by a red ellipse. The centroid of the target is represented as a red star in Figure 10b, whereas Figure 10c showcases the major axis length through a red dotted line. Orientation feature β is presented in Figure 10d.

4.5.2. Trajectory Features

(1)
Rotation angle
Figure 11 displays the orientation feature of false alarms. The red arrow indicates the direction of rotation, while the yellow point corresponds to the centroid and the green point denotes the endpoint of the ellipse. The figure reveals that the false alarm exhibits near circular motion in the image sequence, where the circle represents the theoretical trajectory. Notably, the orientation angle around the endpoint has significant changes, whereas the endpoint remains stationary. Conversely, when the moving target is observed within the limited frames, its trajectory can be approximated as nearly linear motion. For instance, as presented in Figure 11b, the moving target moves toward the endpoint in a linear direction without rotation, resulting in minimal change in the red orientation angle.
In practical scenarios, the discrepancies in the target shape and size within image sequences can result in centroid errors that lead to inaccurate orientation angle estimation. Consequently, the utilization of the aforementioned orientation feature is not a viable approach. As a solution, the recalibration of the angle via basic features is necessary. Hence, the accumulation of the orientation angle was adopted for calculating the rotation angle. Given the nth and first orientation angles of the tracked ith cluster in an image sequence, labeled as θ i ( n ) and θ i ( 1 ) , respectively, the corresponding rotation angle, denoted as ϕ i , can be obtained through the following equation:
ϕ i ( n ) = | θ i ( n ) θ i ( 1 ) |   ,   i = 1 , 2 , 3 ,   n = 1 , 2 , 3 N
The ϕ i ( n ) of the false alarm evidently changes along with the frame number, while the ϕ i ( n ) of the moving target is around zero.
(2)
Moving distance
Another crucial feature for discerning between two signal types is the distance traveled during the image sequence. Due to the restricted rotational motion within a relatively confined area, the false alarm exhibits limited displacement. In contrast, the moving target displays much longer travel distances, as demonstrated in the preceding section.
The moving distance feature is depicted in Figure 12, comprising illustrations of both false alarm and moving target cases in Figure 12a,b, respectively. This attribute is characterized by determining the standard deviation of the endpoint positions of the target across the image sequence. Specifically, as illustrated in Figure 12a, the endpoint of the ellipse follows a circular trajectory around the centroid with a radius of r. Consequently, the resulting moving distance, which corresponds to the standard deviation, is small. In contrast, as presented in Figure 12b, the moving target endpoint undergoes a considerably longer displacement due to the distinctive ‘W’-shaped pattern identified in the previous section.
The ith target’s position in the nth image is represented by p i ( n ) . The two-dimensional mean value of the trajectory’s length is denoted as p m , and it reflects the center coordinate of the trajectory. Furthermore, the moving distance or standard deviation, defined as σ , can be expressed as:
{ σ = [ p i ( n ) - p m ] 2 N 1   ,   i = 1 , 2 , 3 ,   n = 1 , 2 , 3 N p m = 1 N 1 p i ( n )   ,   i = 1 , 2 , 3 ,   n = 1 , 2 , 3 N
After extracting trajectory features, the next step was to use these features to separate the moving target from the false alarm, so as to remove the false alarm.

4.6. Apply Trajectory Feature for Discrimination

The last step was using the rotation angle and moving distance to distinguish the false alarm and moving target, and suppress the false alarm.
In real data, two distinct signature types may have changes in the shape or size within the image sequence. These changes can be attributed to two causes. Firstly, azimuth viewing angle gradual changes may change the projection of the building. Secondly, the LBS and associated processing may either introduce or reduce the target occupied area, leading to morphological alterations. Consequently, it is essential to establish appropriate thresholds for discrimination.
Two threshold values were established for both features. The threshold for the rotation angle ϕ is denoted as ϕ 0 , and the moving distance σ is denoted as σ 0 . The corresponding formulas are given below.
{ ϕ ( n ) > ϕ 0 ,   σ < σ 0   false   alarm ϕ ( n ) < ϕ 0 ,   σ > σ 0   moving   target
The thresholds applied in this study were established via experimental testing. The quantities analysis on how to determine them is one of the planned future works. Once the moving target and false alarm signals are distinguished between based on these criteria, the removal process can be conducted. As the signals are clustered and tracked, they are assigned a unique index to facilitate subsequent processing. Hence, static high target corresponding signals can be easily excluded from the final output detection results.

5. Experiment

5.1. Processing Results

5.1.1. LBS Result

The initial step involved obtaining an image sequence comprising 607 frames, which corresponded to the overlapped subaperture images. For illustrative purposes, frame number 171 was chosen as a representative example. Figure 13 displays the results of the LBS process on this frame.
Figure 13a displays the 171st frame. The extracted background image is displayed in Figure 13b and was obtained with a median filter in the image index dimension. As both the target and shadow signals manifested as significantly deviating from the pixel values of the static scene, the median filter proved effective in extracting the background image. It becomes apparent in Figure 13b that the target and shadow signals have disappeared, while the rotated billboard only covers a relatively small proportion of its base compared to Figure 13a. Subtraction of the original intensity image from the background extracted via LBS retains the signature of both the target and shadow along with the billboard, as illustrated in Figure 13c. Since the original LBS algorithm focused on detecting high-value targets rather than low-value shadows, the last step involved performing CFAR processing to generate coarse detection results, presented in Figure 13d. Upon examination of Figure 13d, one can observe that the static high target signal and moving target signals remain in the image, thus they were inputted into the subsequent steps for further processing.

5.1.2. Morphological Processing

Following CFAR detection, morphological processing was performed on the binary image. Firstly, erosion was employed to eliminate any remaining clutter present in the image. Subsequently, dilation was executed to repair fractured regions within the image. Figure 14a displays the 171st frame image in the absence of erosion/dilation procedures. Figure 14b illustrates the outcome of implementing both procedures on the same image. Evidently, application of these morphological operators effectively removes the clutter situated within the right-hand side image and repairs the minor abnormalities surrounding both the moving target and false alarm signals.

5.1.3. DBSCAN Clustering

The subsequent step undertook the application of a DBSCAN; it has two key parameters: Eps and MinPts. In this study, we set Eps to 0.07 and MinPts to 15 pixels based on prior testing experiences. The effectiveness of the DBSCAN was validated with the clustering result obtained from analyzing the 171st frame, represented in Figure 15. Specifically, the moving target and billboard are denoted by green boxes and red circles, respectively, as depicted in (a). Utilizing the DBSCAN, the extracted data are further presented in (b), with each signature accommodating allocation within unique clusters possessing distinct indices. From the results, the size and shape characteristics of false alarms remain unchanged before and after undergoing the clustering process, which corroborates the effectiveness of the clustering step.

5.1.4. Target Tracking

Subsequent to the clustering step, we established associations between identical targets situated within consecutive frames in the image sequence via tracking algorithms. Figure 16 displays the tracking results of a false alarm situated in frames 171 and 305. The false alarm is designated by a red rectangle, with its trajectory exhibited by the red line present in (b) that was drawn based upon the center of mass calculation. As demonstrated through these results, the false alarm can be tracked across multiple frames, utilizing a unique index for proper identification.

5.1.5. Trajectory Feature Extraction and Discrimination

In this section, we present the results related to trajectory features. Figure 17a showcases the rotation angle pertaining to the false alarm signal, superimposed with two adjacent frames (172 and 173). The equivalent ellipse extracted from these frames is represented by a red dotted line with the centroid indicated by a red star and the endpoint signified by a green dot. Subfigure (b) shows the feature observed in the rotation angle associated with the moving target. The extraction outcomes demonstrate the feature has detectable variance within the key trajectory features distinctive between the two signals, even across adjacent frames such as 172 and 173.
After extracting the aforementioned features on a frame-by-frame basis, we present the curves relating the rotation angle as a function of frame number for both the false alarm and moving target in Figure 18a,b, respectively. The results indicate a noticeable increase in the rotation angle exhibited by the false alarm with an increasing frame count. Notably, some frames demonstrate minor declines in the rotation angle attributable to centroid errors arising from irregular shapes during angle extraction, notwithstanding the underlying trend remaining consistent. Cumulatively, the false alarm experienced an approximate 60° rotation across a span of 134 frames.
Due to the data having a limited scene size coupled with the rapid motion of the moving target, the data presence within the scene was confined to a mere 10 frames. Hence, we utilized all 10 frames for further analysis purposes. In Figure 18b, the rotation angle for the moving target is changed a little. We established an initial threshold level for a cumulative rotation angle at 60°.
The following analysis pertains to the moving distance feature. As delineated in Figure 18, the extraction of endpoint frames for moving targets served as information for our subsequent calculations. Utilizing Equation (13), we calculated the standard deviation values associated with the moving distance features of both false alarm and moving target signals to be 3.8093 m and 23.2384 m, respectively, as depicted within Table 4.
The threshold of the motion distance was set as 5 m; if the moving distance is lower than 5 m, it is judged as a false alarm. Therefore, the moving target and false alarm can be distinguished.

5.1.6. False Alarm Suppression

After distinguishing the moving target and false alarm signals based on their trajectory features, we proceeded with suppression of the false alarm via its unique tracking index. The effectiveness of this technique is exemplified within Figure 19, where (a) depicts the original CSAR image captured at the 171st frame, (b) showcases the results after background subtraction and morphological processing, comprising both moving target and static high target signals, and (c) illustrates the subsequent outcomes upon suppressing the false alarm signal. Of note, signature signals associated with the billboard are no longer present in the results, further validating the efficacy of our proposed methodology.

5.1.7. Method Comparison Experiment

For this section, we conducted the comparison experiment between the proposed method, original logarithm background subtraction (LBS) method, and GoDec method. As described in the previous section, current research has not focused on the static high target’s rotation motion-induced false alarm rate rising problem. As will be shown in this subsection, for such image sequence-based methods, the false alarm should be dealt with, which demonstrates the necessity of the proposed method. Here, we focus on the processing results of the GoDec method, then give the comparison of the final detection results of the proposed method, original LBS method, and GoDec algorithm.
The details of the GoDec method can be found in references [26,33]. We use the same parameters as in 5.2.1 for processing with the GoDec method, i.e., the same image sequence to filter the background and foreground. However, GoDec and LBS have one unique difference. For a given image sequence, the LBS can only have one background and all frames have to subtract it to obtain the foreground. The GoDec method uses all frame data to estimate every image’s low-rank part and sparse part with optimization. The low-rank part represents the stable part, i.e., the background, and the sparse part represents the fast-change part, i.e., the moving target and rotated static high target.
The GoDec result of the 171st, 221st, and 271st frames is given in Figure 20. The three columns from the left to right are the original image (the corresponding background and foreground). The three rows from the top to bottom are the 171st, 221st, and 271st frames, respectively. From the original images, we can see the evident rotation phenomenon and in the 171st frame, we can see the moving target. In the corresponding background in the second column, the image is quite like the original images while the noise is reduced and the moving target in the 171st frame is removed. And the billboard target has a slight change compared with the original image. In foreground images, we can see that most of the stable backgrounds are eliminated in three frames and the moving target is well persevered in the 171st frame. This is the basis for moving target detection. On the other hand, the billboard is also preserved, which proves the necessity of the false alarm removal still existing in the GoDec method. Then, we determined the final detection results of the three methods. The GoDec method and original LBS method conduct the CFAR detection step to output the binary detection result. The proposed method further distinguishes between the moving target and false alarm. The comparison results are given in Figure 21.
In Figure 21, the three columns from the left to right are the binary detection results of the GoDec method, the original LBS method, and the proposed method. We can clearly see that the moving target in the three frames is always detected. The GoDec method and proposed method have less residual clutter. This is because the GoDec method can reduce background and the proposed method has the morphological filter step to deal with it. The static high target, i.e., the billboard, is removed only in the proposed method, as labeled with a red circle. From the comparison experiments illustrated in the figure, it can be seen that the first two methods cannot remove the false alarm, and the proposed false alarm suppression algorithm in this paper can effectively remove it. It shows the better performance and necessity of the proposed method.

6. Conclusions

This paper proposed a new trajectory feature-based method to suppress the false alarm generated by a static high target in CSAR moving target detection. It first thoroughly analyzed the moving target and static high target’s difference for the trajectory in the image sequence. Two new trajectory features of the rotation angle and moving distance were proposed to differentiate them. Based on the features, a new false alarm suppression method was proposed. The method first utilizes LBS to obtain coarse binary detection results comprising both moving and static high targets, then employs morphological filtering to eliminate noise. Next, DBSCAN and target tracking steps are employed to extract the trajectory features of the target and false alarm. Finally, false alarms are suppressed with trajectory-based feature discriminators to output detection results. The W-band CSAR open dataset was used to validate the proposed method’s effectiveness.
Future works can include developing the moving target’s shadow detection method, conducting a W-band CSAR experiment for a threshold selection principle analysis, and thoroughly analyzing the performance of the proposed method.

Author Contributions

Under the supervision of W.S., F.D. performed the experiments and methods. W.S. and F.D. wrote the manuscript. Y.W., J.S., Y.L. (Yang Li), Y.L. (Yun Lin), W.J. and S.W. gave valuable advice on manuscript writing. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the National Natural Science Foundation of China under grants 62201011 (youth program), 61971456 (general program), and 62131001 (key program), and the R&D Program of the Beijing Municipal Education Commission (grant: KM202210009004).

Data Availability Statement

Not applicable.

Acknowledgments

We thank the good advice from anonymous reviewers.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Frolind, P.O.; Ulander, L.M.H.; Gustavsson, A.; Stenstrom, G. VHF/UHF band SAR imaging using circular tracks. In Proceedings of the IEEE International Geoscience and Remote Sensing Symposium, Munich, Germany, 22–27 July 2012; pp. 7409–7411. [Google Scholar]
  2. Kim, S.; Yu, J.; Jeon, S.Y.; Dewantari, A.; Ka, M.H. Signal Processing for a Multiple-Input, Multiple-Output (MIMO) Video Synthetic Aperture Radar (SAR) with Beat Frequency Division Frequency-Modulated Continuous Wave (FMCW). Remote Sens. 2017, 9, 491. [Google Scholar] [CrossRef] [Green Version]
  3. Zhang, J.; Suo, Z.; Li, Z.; Zhang, Q. DEM Generation Using Circular SAR Data Based on Low Rank and Sparse Matrix Decomposition. IEEE Geosci. Remote Sens. Lett. 2018, 15, 724–728. [Google Scholar] [CrossRef]
  4. Cerutti-Maori, D.; Klare, J.; Brenner, A.R.; Ender, J.H.G. Wide-area traffic monitoring with the SAR/GMTI system PAMIR. IEEE Trans. Geosci. Remote Sens. 2008, 46, 3019–3030. [Google Scholar] [CrossRef]
  5. Damini, A.; Balaji, B.; Parry, C.; Mantle, V. A video SAR mode for the X-band wide band experimental air borne radar. In Proceedings of the Algorithms for Synthetic Aperture Radar Imagery XVII, Orlando, FL, USA, 18 April 2010; p. 76990E. [Google Scholar]
  6. Bao, J.; Zhang, X.; Zhang, T.; Shi, J.; Wei, S. A Novel Guided Anchor Siamese Network for Arbitrary Target-of-Interest Tracking in Video-SAR. Remote Sens. 2021, 13, 4504. [Google Scholar] [CrossRef]
  7. Mapelli, D.; Guarnieri, A.M.; Giudici, D. Generation and Calibration of High-Resolution DEM From Single-Baseline Spaceborne Interferometry: The “Split-Swath” Approach. IEEE Trans. Geosci. Remote Sens. 2013, 52, 4858–4867. [Google Scholar] [CrossRef]
  8. Jao, J.K. Theory of synthetic aperture radar imaging of a moving target. IEEE Trans. Geosci. Remote Sens. 2001, 39, 1984–1992. [Google Scholar]
  9. Lightstone, L.; Faubert, D.; Rempel, G. Multiple phase centre DPCA for airborne radar. In Proceedings of the 1991 IEEE National Radar Conference, Los Angeles, CA, USA, 12–13 March 1991; pp. 36–40. [Google Scholar]
  10. Makhoul, E.; Baumgartner, S.V.; Jager, M.; Broquetas, A. Multichannel SAR-GMTI in Mari time Scenarios with F-SAR and Terra SAR-X Sensors. IEEE J. Sel. Top. Appl. Earth Obs. Remote Sens. 2015, 8, 5052–5067. [Google Scholar] [CrossRef]
  11. Kirscht, M. Detection and imaging of arbitrarily moving targets with single-channel SAR. IEEE Proc. Radar Sonar Navig. 2003, 150, 7–11. [Google Scholar] [CrossRef]
  12. Ouchi, K. On the multilook images of moving targets by synthetic aperture radars. IEEE Trans. Antennas Propag. 1985, 33, 823–827. [Google Scholar] [CrossRef] [Green Version]
  13. Gierull, C.H. Statistical analysis of multilook SAR interferograms for CFAR detection of ground moving targets. IEEE Trans. Geosci. Remote Sens. 2004, 42, 691–701. [Google Scholar] [CrossRef]
  14. Henke, D.; Dominguez, E.M.; Small, D.; Schaepman, M.E.; Meier, E. Moving target tracking in single and multichannel SAR. IEEE Trans. Geosci. Remote Sens. 2015, 53, 3146–3159. [Google Scholar] [CrossRef]
  15. Perlovsky, L.; Ilin, R.; Deming, R.; Linnehan, R.; Lin, F. Moving target detection and characterization with circular SAR. In Proceedings of the IEEE Radar Conference, Washington, DC, USA, 10–14 May 2010; pp. 661–666. [Google Scholar]
  16. Henke, D.; Magnard, C.; Frioud, M.; Small, D.; Meier, E.; Schaepman, M.E. Moving-target tracking in single-channel wide-beam SAR. IEEE Trans. Geosci. Remote Sens. 2012, 50, 4735–4747. [Google Scholar] [CrossRef] [Green Version]
  17. Tian, X.; Liu, J.; Mallick, M. Simultaneous Detection and Tracking of Moving-Target Shadows in ViSAR Imagery. IEEE Trans. Geosci. Remote Sens. 2020, 59, 1182–1199. [Google Scholar] [CrossRef]
  18. Wang, H.; Chen, Z.; Zheng, S. Preliminary research of low RCS moving target detection based on Ka-band Video-SAR. IEEE Geosci. Remote Sens. Lett. 2017, 14, 811–815. [Google Scholar] [CrossRef]
  19. Ying, Z.; Daiyin, Z.; Xiang, Y.; Mao, X. Approach to moving targets shadow detection for Video SAR. J. Electron. Inf. Technol. 2017, 39, 2197–2202. [Google Scholar]
  20. Wang, K.; Qian, J.; Huang, S.; Xing, M.; Yang, X.; Yang, L. Single channel sar slow moving target detection in time-frequency domain. In Proceedings of the 2019 6th Asia-Pacific Conference on Synthetic Aperture Radar (APSAR), Xiamen, China, 26–29 November 2019. [Google Scholar]
  21. Yi, H.; Fan, C.Y.; Yang, J.G.; Huang, X.T. Imaging and locating multiple ground moving targets based on keystone transform and FrFT for single channel SAR system. In Proceedings of the 2009 2nd Asian-Pacific Conference on Synthetic Aperture Radar, Xian, China, 26–30 October 2009; pp. 771–774. [Google Scholar] [CrossRef]
  22. Huang, Y.; Liao, G.; Li, J.; Xu, J. Ground moving target indication and parameter estimation with single channel for SAR system. In Proceedings of the 2016 IEEE International Geoscience and Remote Sensing Symposium (IGARSS), Beijing, China, 10–15 July 2016; pp. 6417–6420. [Google Scholar] [CrossRef]
  23. Palm, S.; Sommer, R.; Janssen, D.; Tessmann, A.; Stilla, U. Airborne Circular W-Band SAR for Multiple Aspect Urban Site Monitoring. IEEE Trans. Geosci. Remote Sens. 2019, 57, 6996–7016. [Google Scholar] [CrossRef]
  24. Shen, W.; Lin, Y.; Yu, L.; Xue, F.; Hong, W. Single Channel Circular SAR Moving Target Detection Based on Logarithm Background Subtraction Algorithm. Remote Sens. 2018, 10, 742. [Google Scholar] [CrossRef] [Green Version]
  25. Shen, W.; Hong, W.; Han, B.; Wang, Y.; Lin, Y. Moving Target Detection with Modified Logarithm Background Subtraction and Its Application to the GF-3 Spotlight Mode. Remote Sens. 2019, 11, 1190. [Google Scholar] [CrossRef] [Green Version]
  26. Zhang, Y.; Zhu, D.; Wang, P.; Zhang, G.; Leung, H. Vision-Based Vehicle Detection for VideoSAR Surveillance Using Low-Rank Plus Sparse Three-Term Decomposition. IEEE Trans. Veh. Technol. 2020, 69, 4711–4726. [Google Scholar] [CrossRef]
  27. Li, Y.; Li, C.; Yang, W.; Wang, Y.; Li, Y. Moving Vehicle Detection Based on RPCA Using Multisquint Spaceborne SAR Images. IEEE Geosci. Remote Sens. Lett. 2022, 19, 4005005. [Google Scholar] [CrossRef]
  28. Bouwmans, T.; Sobral, A.; Javed, S.; Jung, S.K.; Zahzah, E.H. Decomposition into low-rank plus additive matrices for background/foreground separation: A review for a comparative evaluation with a large-scale dataset. Comput. Sci. Rev. 2017, 23, 1–71. [Google Scholar] [CrossRef] [Green Version]
  29. Feng, S.; Lin, Y.; Wang, Y.; Yang, Y.; Shen, W.; Teng, F.; Hong, W. DEM Generation with a Scale Factor Using Multi-Aspect SAR Imagery Applying Radar Grammetry. Remote Sens. 2020, 12, 556–573. [Google Scholar] [CrossRef] [Green Version]
  30. An, D.; Wang, W.; Zhou, Z. Refocusing of Ground Moving Target in Circular Synthetic Aperture Radar. IEEE Sens. J. 2019, 19, 8668–8674. [Google Scholar] [CrossRef]
  31. Li, J.; An, D.; Wang, W.; Zhou, Z.; Chen, M. A Novel Method for Single-Channel CSAR Ground Moving Target Imaging. IEEE Sens. J. 2019, 19, 8642–8649. [Google Scholar] [CrossRef]
  32. Shen, W.; Lin, Y.; Chen, S.; Xue, F.; Yang, Y.; Hong, W. Apparent Trace Analysis of Moving Target with Linear Motion in Circular SAR Imagery. In Proceedings of the EUSAR 2018: 12th European Conference on Synthetic Aperture Radar, Aachen, Germany, 4–7 June 2018; pp. 1–4. [Google Scholar]
  33. Zhou, T.; Tao, D. GoDec: Randomized Lowrank & Sparse Matrix Decomposition in Noisy Case. In Proceedings of the 28th International Conference on Machine Learning, ICML 2011, Bellevue, WA, USA, 28 June–2 July 2011; pp. 33–40. [Google Scholar]
Figure 1. Example image of W-band CSAR subaperture data, and the highly stationary target is in the red circle.
Figure 1. Example image of W-band CSAR subaperture data, and the highly stationary target is in the red circle.
Remotesensing 15 03164 g001
Figure 2. W-band CSAR data.
Figure 2. W-band CSAR data.
Remotesensing 15 03164 g002
Figure 3. Static high target and moving targets in CSAR image sequence: (a) optical image of billboard; (b) rotating motion of billboard at different times, and its in the red circle (i.e., different azimuth viewing angle).
Figure 3. Static high target and moving targets in CSAR image sequence: (a) optical image of billboard; (b) rotating motion of billboard at different times, and its in the red circle (i.e., different azimuth viewing angle).
Remotesensing 15 03164 g003
Figure 4. Moving target detection results after logarithm background subtraction. Green box is moving target, red circle is billboard. (a) Original subaperture image; (b) detection result.
Figure 4. Moving target detection results after logarithm background subtraction. Green box is moving target, red circle is billboard. (a) Original subaperture image; (b) detection result.
Remotesensing 15 03164 g004
Figure 5. CSAR geometry: (a) 3D geometry; (b) front view of A P P B B .
Figure 5. CSAR geometry: (a) 3D geometry; (b) front view of A P P B B .
Remotesensing 15 03164 g005
Figure 6. Simulation results: (a) subaperture image at azimuth viewing angle of 45°; (b) the combination of all 24 subaperture images.
Figure 6. Simulation results: (a) subaperture image at azimuth viewing angle of 45°; (b) the combination of all 24 subaperture images.
Remotesensing 15 03164 g006
Figure 7. Comparison of trajectories for static high target and moving target: (a) combination of subaperture images of static high target; (b) combination of all subaperture images of static high target; (c) simulation result of full aperture moving target along azimuth direction; (d) full aperture simulation of moving target along range direction.
Figure 7. Comparison of trajectories for static high target and moving target: (a) combination of subaperture images of static high target; (b) combination of all subaperture images of static high target; (c) simulation result of full aperture moving target along azimuth direction; (d) full aperture simulation of moving target along range direction.
Remotesensing 15 03164 g007
Figure 8. Flow chart of false alarm suppression algorithm.
Figure 8. Flow chart of false alarm suppression algorithm.
Remotesensing 15 03164 g008
Figure 9. Coarse detection results: (a) false alarm (red) and clutter (yellow); (b) moving target (green) and clutter (yellow).
Figure 9. Coarse detection results: (a) false alarm (red) and clutter (yellow); (b) moving target (green) and clutter (yellow).
Remotesensing 15 03164 g009
Figure 10. Four basic features: (a) equivalent ellipse; (b) centroid; (c) major axis length; (d) orientation.
Figure 10. Four basic features: (a) equivalent ellipse; (b) centroid; (c) major axis length; (d) orientation.
Remotesensing 15 03164 g010
Figure 11. Illustration of rotation angle feature: (a) false alarm; (b) moving target.
Figure 11. Illustration of rotation angle feature: (a) false alarm; (b) moving target.
Remotesensing 15 03164 g011
Figure 12. Illustration of moving distance feature: (a) false alarm; (b) moving target.
Figure 12. Illustration of moving distance feature: (a) false alarm; (b) moving target.
Remotesensing 15 03164 g012
Figure 13. Background subtraction detection results: (a) original SAR image; (b) background image; (c) forehead image; (d) After LBS image.
Figure 13. Background subtraction detection results: (a) original SAR image; (b) background image; (c) forehead image; (d) After LBS image.
Remotesensing 15 03164 g013
Figure 14. Morphological processing results: (a) original image; (b) after dilation and erosion.
Figure 14. Morphological processing results: (a) original image; (b) after dilation and erosion.
Remotesensing 15 03164 g014
Figure 15. Clustering result of 171st frame: (a) before DBSCAN; (b) after DBSCAN.
Figure 15. Clustering result of 171st frame: (a) before DBSCAN; (b) after DBSCAN.
Remotesensing 15 03164 g015
Figure 16. Target tracking results: (a) tracking results of 171st frame; (b) tracking results of 305th frame.
Figure 16. Target tracking results: (a) tracking results of 171st frame; (b) tracking results of 305th frame.
Remotesensing 15 03164 g016
Figure 17. Rotation angle extraction results. (a) False alarm results of adjacent frames 172 and 173; (b) moving target results of adjacent frames 172 and 173.
Figure 17. Rotation angle extraction results. (a) False alarm results of adjacent frames 172 and 173; (b) moving target results of adjacent frames 172 and 173.
Remotesensing 15 03164 g017
Figure 18. The curve of target rotation angle changing with frame number: (a) false alarm; (b) moving target.
Figure 18. The curve of target rotation angle changing with frame number: (a) false alarm; (b) moving target.
Remotesensing 15 03164 g018
Figure 19. False alarm suppression results, a moving target, represented by a green circle, and a static high target, designated in a red circle. 171st frame: (a) original CSAR image; (b) original detection result; (c) false alarm suppression result.
Figure 19. False alarm suppression results, a moving target, represented by a green circle, and a static high target, designated in a red circle. 171st frame: (a) original CSAR image; (b) original detection result; (c) false alarm suppression result.
Remotesensing 15 03164 g019
Figure 20. The GoDec result on frame 171, 221, and 271.
Figure 20. The GoDec result on frame 171, 221, and 271.
Remotesensing 15 03164 g020
Figure 21. Final detection results of three methods. Columns from left to right represent GoDec, original logarithm background subtraction (LBS), and the proposed method on frames 171, 221, and 271, respectively. The red circle labels the static high target false alarm.
Figure 21. Final detection results of three methods. Columns from left to right represent GoDec, original logarithm background subtraction (LBS), and the proposed method on frames 171, 221, and 271, respectively. The red circle labels the static high target false alarm.
Remotesensing 15 03164 g021
Table 1. Parameters of W-band CSAR data.
Table 1. Parameters of W-band CSAR data.
ParametersDescription
Center Frequency94 GHz
Flight Height300 m
Resolution15 cm × 12 cm
Frame Rate20 Hz
Total Frames607
Subaperture Overlap Rate80%
Table 2. Parameters for simulation.
Table 2. Parameters for simulation.
ParametersValue
Aircraft height (m)7251
Aircraft radius (m)467
Target A height (m)0
Target B height (m) 5
Incidence angle (°)15
Subaperture degree (°)2
Table 3. Four basic features’ description.
Table 3. Four basic features’ description.
Basic Feature NameDescription
Equivalent ellipseMinimum bounding ellipse with equal pixels
CentroidCenter of mass
Major axis lengthThe major axis of the ellipse
OrientationThe angle between ellipse’s major axis and X-axis
Table 4. Moving distance (standard deviation) of false alarm and moving target.
Table 4. Moving distance (standard deviation) of false alarm and moving target.
False AlarmMoving Target
Moving distance (m)3.809323.2384
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Shen, W.; Ding, F.; Wang, Y.; Li, Y.; Sun, J.; Lin, Y.; Jiang, W.; Wang, S. Static High Target-Induced False Alarm Suppression in Circular Synthetic Aperture Radar Moving Target Detection Based on Trajectory Features. Remote Sens. 2023, 15, 3164. https://doi.org/10.3390/rs15123164

AMA Style

Shen W, Ding F, Wang Y, Li Y, Sun J, Lin Y, Jiang W, Wang S. Static High Target-Induced False Alarm Suppression in Circular Synthetic Aperture Radar Moving Target Detection Based on Trajectory Features. Remote Sensing. 2023; 15(12):3164. https://doi.org/10.3390/rs15123164

Chicago/Turabian Style

Shen, Wenjie, Fan Ding, Yanping Wang, Yang Li, Jinping Sun, Yun Lin, Wen Jiang, and Shuo Wang. 2023. "Static High Target-Induced False Alarm Suppression in Circular Synthetic Aperture Radar Moving Target Detection Based on Trajectory Features" Remote Sensing 15, no. 12: 3164. https://doi.org/10.3390/rs15123164

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop