CN111796265A - Low-speed target detection method based on improved area bilateral smoothing - Google Patents

Low-speed target detection method based on improved area bilateral smoothing Download PDF

Info

Publication number
CN111796265A
CN111796265A CN202010536205.7A CN202010536205A CN111796265A CN 111796265 A CN111796265 A CN 111796265A CN 202010536205 A CN202010536205 A CN 202010536205A CN 111796265 A CN111796265 A CN 111796265A
Authority
CN
China
Prior art keywords
target
area
doppler velocity
improved
low
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Pending
Application number
CN202010536205.7A
Other languages
Chinese (zh)
Inventor
孟凡
曹阳
贾倩茜
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
724th Research Institute of CSIC
Original Assignee
724th Research Institute of CSIC
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by 724th Research Institute of CSIC filed Critical 724th Research Institute of CSIC
Priority to CN202010536205.7A priority Critical patent/CN111796265A/en
Publication of CN111796265A publication Critical patent/CN111796265A/en
Pending legal-status Critical Current

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S13/00Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
    • G01S13/66Radar-tracking systems; Analogous systems
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S7/00Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/02Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00

Abstract

The invention discloses a low-speed target detection method based on improved area bilateral smoothing. The method is suitable for designing a detection area by taking a target extrapolation position as a center aiming at a key focus target under a noise environment, then preprocessing radar data and converting the radar data into image data, adopting an improved area bilateral smoothing algorithm, proposing a Gaussian function value corresponding to Doppler velocity as a weight term of a value domain, combining the weight term with a weight coefficient kernel of a spatial distance, enhancing the processing of the target area and improving the contrast, then extracting the contour of the target area from the background by using a self-adaptive threshold segmentation method, and carrying out target condensation by combining the Doppler velocity.

Description

Low-speed target detection method based on improved area bilateral smoothing
Technical Field
The invention belongs to the field of radar data processing.
Background
Due to the radar system and the influence of the target background environment, noise and clutter are often accompanied in the echo signal receiving process. The method is used for detecting the low-speed target under the clutter background, on one hand, radar clutter suppression performance needs to be improved, and on the other hand, echo energy needs to be accumulated.
In the traditional detection of moving targets, adaptive matched filtering methods are often adopted, such as a Generalized Likelihood Ratio detector (GLRT) and an adaptive matched filtering detector (AMF), and coherent accumulation can be realized through Doppler filtering, so that the signal-to-noise Ratio of target echoes is effectively improved, and the capability of detecting moving targets is enhanced. However, for low-speed targets, especially in clutter environment, the conventional signal processing method such as MTI or MTD is difficult to detect such targets because the MTI filter has zero-frequency notch depth and width, which cause the snr loss of the low-speed targets after filtering, and the constant false alarm processing after detecting the MTD generates false alarm due to clutter. The method is used for detecting the low-speed target in the clutter from the angle of image processing, firstly, region enhancement is carried out, the contrast between the target and the clutter is improved, the target is separated from the background to the greatest extent (the target is a foreground, and the clutter is the background), the target is completely extracted, and meanwhile, noise is removed as far as possible. When the target is in the clutter area, a traditional radar target detection method is adopted, a strict detection threshold is set, the clutter removal effect is good, but the target is also removed, and the detection fails. The invention introduces the image processing correlation theory into the radar target detection field aiming at the echo data after signal processing, adopts the improved bilateral region smoothing algorithm to detect the target in the detection region, simultaneously considers the airspace information and the speed similarity, and enhances the contrast ratio of the target region. Not only can extract the target, but also can achieve the effect of denoising.
Disclosure of Invention
The invention mainly solves the problem of detecting the low-speed moving target in the complex environment.
Aiming at a key focus target, a target extrapolation position is taken as a center, a detection area is designed according to the change of Doppler velocity, an improved bilateral smoothing technology is provided for carrying out target detection in the area, the Doppler velocity corresponding to each point is used for replacing the amplitude as a weight item of a Gaussian kernel function, weights are calculated from two angles of space proximity and velocity similarity respectively, a convolution template is obtained, the detection area is smoothed, and then a self-adaptive threshold value segmentation method is adopted for extracting a target contour to carry out target condensation.
Bilateral smoothing is based on gaussian filtering, which has the same smoothing degree in each direction, and does not change the edge trend of the original image, and simultaneously, the gaussian kernel functions are single-valued functions, which monotonically decrease in all directions (consistent with target characteristics), and the position coordinates of each position, i.e. the weight of each amplitude point, corresponds to the function value on the two-dimensional gaussian function, when the value of σ is small, as shown in fig. 1: when x is less than-1 and x is greater than 1, the value of f (x) is already close to 0, which means that the amplitude value of the more distant position has no meaning for calculating the amplitude of the current position, namely, the purpose of highlighting the target area is achieved.
In order to realize the purpose of extracting a target and removing noise, a double-Gaussian filter is adopted for bilateral smoothing, in an image flat area, the pixel value change is very small, the corresponding pixel range domain weight is close to 1, and at the moment, the spatial domain weight plays a main role and is equivalent to Gaussian blur; in the edge area of the image, the pixel value is greatly changed, the pixel range area weight is increased, so that the information of the edge is kept, and the Euclidean distance of the target amplitude is considered, and the difference between the amplitude in the detection area and the central amplitude is also considered.
However, in actual data, due to the influence of fluctuation of the amplitude, the amplitude similarity cannot be accurately used as a basis for target detection. Therefore, the invention uses the corresponding Doppler velocity to replace the amplitude as the weight item of the Gaussian kernel function, and calculates the weight from two angles of space proximity and velocity similarity respectively, thereby obtaining the convolution template and realizing the smoothness of the detection area.
When radial velocity exists between the target and the radar, Doppler frequency shift is generated on the receiver, so that the target velocity can be obtained by using the Doppler change rate, and the moving target is detected. By definition, the Doppler frequency fdThe relationship to radial velocity is:
Figure BDA0002537074680000021
wherein v isrIn the radial velocity, λ is the wavelength, f is the frequency of the target, and c is the propagation velocity of the electromagnetic wave, as shown in fig. 2.
Therefore, the target in the clutter region is enhanced from the space domain and the value domain, the target extrapolation position is the center of the Gaussian function, as shown in formula (3), namely, in the neighborhood, the point with the speed value closer to the speed value of the central point has larger weight, and the point with the larger speed value difference has smaller weight. The weight is determined by a value domain Gaussian function, and the weight coefficients of the two are multiplied to obtain the final convolution template.
(1) Spatial distance: the Euclidean distance between the current point and the central point is referred, and the spatial domain Gaussian function is as follows:
Figure BDA0002537074680000022
wherein (x)i,yi) Is the current point position, (x)c,yc) σ is the spatial domain standard deviation for the center point position.
(2) Speed distance: the absolute value of the difference between the current point speed and the central point speed is indicated, and the value domain Gaussian function is as follows:
Figure BDA0002537074680000023
wherein, doppler (x)i,yi) Current point speed value, doppler (x)c,yc) For center point velocity values, σ is the range standard deviation.
The width of the gaussian function is characterized by a parameter σ, the value of which determines the amplitude of the variation of the gaussian function, the larger σ the wider the band of the gaussian filter and the better the smoothness, preferably a 3 × 3 template σ value of 0.85.
According to the result of the enhancement of the target area, the target contour is extracted by adopting an adaptive threshold segmentation method, the method considers that an image consists of a foreground and a background or two groups of gray values, the gray values are distributed into peak states on a histogram, and the low valley between the two peaks is the threshold value of the image segmentation.
And finally, performing target condensation according to a threshold segmentation result, wherein the traditional radar target detection method is to calculate the target centroid by condensation by utilizing amplitude information. Because the amplitude information is easily interfered by static noise and the detection probability is often reduced, the method introduces Doppler velocity to carry out target condensation aiming at the detected target contour area.
Figure BDA0002537074680000031
Figure BDA0002537074680000032
Wherein D is1nFor successive range cell values, A, of the orientation within the target profilenAmplitude value after azimuthal agglomeration within the target contour, fnIs the maximum Doppler velocity at a certain position in the target area in the target contour, and the Doppler velocity of the target after f' condensation.
The method of the invention is adopted to detect the focus target, can realize the enhancement of the target area under the complex environment, completely reserve the outline information of the target, remove the clutter and provide accurate measurement data for the target tracking.
Drawings
Fig. 1 is a diagram of doppler velocity.
Fig. 2 is a gaussian kernel function curve.
FIG. 3 is a flowchart of a low-speed target detection method based on improved bilateral smoothing of a region.
FIG. 4 is a diagram illustrating the effect of region object enhancement and threshold segmentation using the method of the present invention.
Fig. 5 is a histogram of the original detection area.
Fig. 6 is a histogram of the detected region after processing by the method of the present invention.
Detailed Description
The technical scheme of the invention is further explained by combining the drawings and the embodiment.
The invention provides a low-speed target detection method based on improved area bilateral smoothing, which comprises the following steps (see fig. 3):
(1) aiming at an important target, a detection area is designed according to the change of Doppler velocity by taking a target extrapolation position as a center, the width and the height are respectively equal to dsin (theta) + kh ', and equal to dcos (theta) + kw', and an alternative unit is formed according to the sub-units. Wherein d is the distance between the filtered value and the extrapolated value of the current period, h 'and w' are the height and width of the target initial detection region, k is a range coefficient, wherein,
Figure BDA0002537074680000033
Δ v is the difference between the radial velocity of the current period and the mean value of the radial velocities of the previous n periods, f is the maximum value of the Doppler velocity in the detection area of the current period,
Figure BDA0002537074680000034
is the mean value of the maximum values of the Doppler velocities of the first n periods, theta=arctan(vy/vx),vxAnd vyRadial velocities in the x and y directions, respectively.
(2) Preprocessing the data in the region and converting the data into image data;
(3) the target extrapolation position is used as the center of a Gaussian function, improved bilateral smoothing is adopted to enhance the target area, the test effect is shown in figure 4, the middle figure is the effect after enhancement treatment, noise and clutter are obviously smoothed, and the target area is completely stored;
(4) and adaptively selecting a threshold value. Fig. 5 is a histogram of an original target region, where only one peak value is present in the histogram, and the target and noise cannot be accurately distinguished, and fig. 6 is a histogram of a target region after an improved bilateral region smoothing process, where two peaks obviously appear in the histogram, a peak with an amplitude value between 30 and 40 represents smoothed background noise, a peak with an amplitude value around 50 is a target, and a trough between two peaks is a selection threshold.
(5) And extracting the target contour according to the result of threshold segmentation. If there are still multiple targets after the segmentation, as shown in the third diagram in fig. 4, it is described that the histogram distribution of the target coincides with the key target, and the target cannot be distinguished by the threshold segmentation, so that the target contour closest to the center of the region is extracted by using the euclidean distance, which is the key target of interest. Then, the original data is mapped according to the position of the target contour, and target agglomeration is performed by using equations (4) and (5).

Claims (2)

1. The low-speed target detection method based on improved area bilateral smoothing is characterized by comprising the following steps: aiming at a key focus target, a target detection area is designed according to the change of Doppler velocity by taking a target extrapolation position as a center, an improved bilateral smoothing technology is provided for carrying out target detection in the area, the Doppler velocity corresponding to each point is used for replacing the amplitude as a weight item of a Gaussian kernel function, weights are calculated from two angles of space proximity and velocity similarity respectively, so that a convolution template is obtained, the detection area is smoothed, and finally a self-adaptive threshold value segmentation method is adopted for extracting a target contour to carry out target condensation.
2. The improved region bilateral smoothing-based low-speed target detection method according to claim 1, characterized in that: target agglomeration is performed using the doppler velocity for the detected target contour region, and is represented by the following equation:
Figure FDA0002537074670000011
Figure FDA0002537074670000012
wherein D is1nFor successive range cell values, A, of the orientation within the target profilenAs the amplitude value after azimuthal agglomeration within the target profile, fnIs the maximum Doppler velocity at a certain position in the target profile, and f' is the Doppler velocity of the target after condensation.
CN202010536205.7A 2020-06-12 2020-06-12 Low-speed target detection method based on improved area bilateral smoothing Pending CN111796265A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010536205.7A CN111796265A (en) 2020-06-12 2020-06-12 Low-speed target detection method based on improved area bilateral smoothing

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010536205.7A CN111796265A (en) 2020-06-12 2020-06-12 Low-speed target detection method based on improved area bilateral smoothing

Publications (1)

Publication Number Publication Date
CN111796265A true CN111796265A (en) 2020-10-20

Family

ID=72803434

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010536205.7A Pending CN111796265A (en) 2020-06-12 2020-06-12 Low-speed target detection method based on improved area bilateral smoothing

Country Status (1)

Country Link
CN (1) CN111796265A (en)

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20090317015A1 (en) * 2008-06-20 2009-12-24 Porikli Fatih M Method for Filtering of Images with Bilateral Filters and Power Images
CN104680553A (en) * 2014-11-26 2015-06-03 沈阳工业大学 Bilateral filtration-based variable kernel function target tracing method
CN107064902A (en) * 2017-05-12 2017-08-18 安徽四创电子股份有限公司 A kind of target condensing method for airport surface detection radar system
CN109100697A (en) * 2018-07-31 2018-12-28 安徽四创电子股份有限公司 A kind of target condensing method based on ground surveillance radar system

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20090317015A1 (en) * 2008-06-20 2009-12-24 Porikli Fatih M Method for Filtering of Images with Bilateral Filters and Power Images
CN104680553A (en) * 2014-11-26 2015-06-03 沈阳工业大学 Bilateral filtration-based variable kernel function target tracing method
CN107064902A (en) * 2017-05-12 2017-08-18 安徽四创电子股份有限公司 A kind of target condensing method for airport surface detection radar system
CN109100697A (en) * 2018-07-31 2018-12-28 安徽四创电子股份有限公司 A kind of target condensing method based on ground surveillance radar system

Non-Patent Citations (5)

* Cited by examiner, † Cited by third party
Title
PETER SAND: "Particle Video: Long-Range Motion Estimation using Point Trajectories", COMPUTER SOCIETY, pages 3 *
刘天时;惠霄霄;: "基于暗通道的图像去雾改进算法研究", 智能计算机与应用, no. 01 *
张业宏;陈恩平;么跃轩;刘宝华;: "基于双边滤波与离散余弦变换的NLM去噪算法", 燕山大学学报, no. 03 *
曹瑛;李志永;卢晓鹏;邹谋炎;: "基于自适应邻域双边滤波的点目标检测预处理算法", 电子与信息学报, no. 08, pages 1910 *
田红彬;: "基于双边滤波和巴氏距离的红外小目标检测", 核电子学与探测技术, no. 10 *

Similar Documents

Publication Publication Date Title
CN109767439B (en) Target detection method for multi-scale difference and bilateral filtering of self-adaptive window
CN106918807B (en) A kind of Targets Dots condensing method of radar return data
CN102222322A (en) Multiscale non-local mean-based method for inhibiting infrared image backgrounds
CN112684428B (en) Multi-target constant false alarm rate detection method based on signal agent
CN109597045B (en) Static target steady recognition method based on two times of clutter suppression
CN111161308A (en) Dual-band fusion target extraction method based on key point matching
CN108647693B (en) Sea surface infrared target detection method based on binary significance characteristics
CN107479037A (en) A kind of PD radar clutters area judging method
CN112669332B (en) Method for judging sea-sky conditions and detecting infrared targets based on bidirectional local maxima and peak value local singularities
CN116953674A (en) Rapid target detection algorithm in sonar imaging
CN111796265A (en) Low-speed target detection method based on improved area bilateral smoothing
CN115797374B (en) Airport runway extraction method based on image processing
KR101813357B1 (en) Clutter elimination Method and Apparatus for applying variable data map and variable statistic characteristics information for identifying underwater fixed target and clutter
CN115828552A (en) Line spectrum purification method based on shore-based passive sonar and optimal path tracking algorithm
CN114429593A (en) Infrared small target detection method based on rapid guided filtering and application thereof
Meng et al. A Novel Method for Low-Speed Dim Small Target Detection
CN114283322A (en) SAR image ship detection method
CN109613489B (en) Radar signal sorting method based on fuzzy function landform characteristics
Ning et al. Ship detection of infrared image in complex scene based on bilateral filter enhancement
Bae et al. Recursive multi-SEs NWTH method for small target detection in infrared images
He et al. Marine multi-target detection based on improved wavelet transform
Faramarzi et al. Infrared dim and small targets detection based on multi-scale local contrast measure utilizing efficient spatial filters
Zhen et al. A Sonar Image Segmentation Algorithm based on Two-Dimensional Spatio-Temporal Fuzzy Entropy
Kuiying et al. Sar automatic target recognition based on shadow contour
Zhang et al. Oil tank detection based on linear clustering saliency analysis for synthetic aperture radar images

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
CB02 Change of applicant information

Address after: 210003 No. 346, Zhongshan North Road, Jiangsu, Nanjing

Applicant after: 724 Research Institute of China Shipbuilding Corp.

Address before: 210003 No. 346, Zhongshan North Road, Jiangsu, Nanjing

Applicant before: 724TH RESEARCH INSTITUTE OF CHINA SHIPBUILDING INDUSTRY Corp.

CB02 Change of applicant information