Summary of the invention
In view of the above-mentioned deficiencies of the prior art, the technical problem to be solved by the present invention is to provide one kind to be based on multi-functional thunder
The intertidal zone monitoring method reached, should intertidal zone monitoring method based on multifunction radar by draw out in radar return land with
Beach boundary separates map, greatly reduces calculation amount and given benchmark is easy to monitor, by coherent accumulation enhance echo and
Preferably inhibit sea clutter by two threshold value thresholding T1, T2 parameters in optimization SIC interframe accumulation algorithm, after then handling
Echo data generate picture in real time using OpenGL;It is real-time using every three seconds reading current images of Canny function in OpenCV
Show the region of current intertidal zone;Finally using interpretation, in-site measurement is manually visualized to coastline data correction, to reach pair
The monitoring of intertidal zone.
In order to solve the above technical problems, the technical solution adopted by the present invention is that:
A kind of intertidal zone monitoring method based on multifunction radar, includes the following steps.
Step 1, in radar return land and beach border map drafting: using TinyXML technology and OpenGL it is more
Side shape lattice method completes the drafting on land and beach border map in radar return.
Step 2, enhance radar return image: radar return image being enhanced using FFT coherent accumulation technology.
Step 3, the inhibition of sea clutter: first two threshold value threshold parameters T1, T2 in SIC interframe accumulation algorithm are carried out excellent
Change, obtains optimization threshold value thresholding;Then, after using the SIC interframe accumulation algorithm with optimization threshold value thresholding to enhancing in step 2
Radar return image data handled, to inhibit sea clutter;The calculation formula of two optimizations threshold value threshold parameter T1, T2
Are as follows:
T1=μ1+σ1*K
T2=μ1+Fα(x)*σ1*K
Wherein, FαIt (x) is the probability density function of envelope, calculation formula are as follows:
In formula, K is by false-alarm probability PfaAnd the constant that the distribution function of video data is derived by.
μ1It is the mean value of noise profile in original video data, σ1For the standard deviation of noise profile in original video data,
In, μ1And σ1Calculation formula are as follows:
In formula, 2 σ2For the mean-square value of clutter amplitude,For the mean value of clutter amplitude, x is sea clutter amplitude instantaneous value.
Step 4, the radar return image data of ocean clutter cancellation processing, benefit coastline real-time rendering: will be passed through in step 3
Current radar return image data are generated into picture with OpenGL, were read using Canny function in OpenCV every three seconds current
Picture simultaneously draws coastline, the region of the current intertidal zone of real-time display.
Further include step 5, coastline data refine just: for, to the position of coastline real-time rendering, utilizing people in step 4
The data in work visual interpretation and the practical coastline of in-site measurement are modified, so as to Tideland resources real-time monitoring.
The calculation method of clutter amplitude mean-square value are as follows: choose arbitrary region, it is assumed that sea clutter meets uniformly in selection area
Characteristic, to be currently detected the range statistics sea clutter characteristic for opening up MxN size centered on unit (i, j);Pass through multi-functional thunder
Up to the instantaneous echo amplitude value of selection area is obtained, the mean-square value of current clutter amplitude is found out.
In step 1, firstly, the initial data to mapping is extracted in initial data using TinyXML technology
LineString element and Polygon element;Then, it using longitude and latitude sample point interpolation technique, obtains in practical map most
Different Road in outer layer closed boundary region and region;Finally, completing side using OpenGL polygon lattice technology
It draws on boundary.
In step 2, the specific method that radar return image is enhanced using FFT coherent accumulation technology, including it is as follows
Step:
Step 21, calculate accumulation umber of pulse N: in azimuth scan radar half-power beam width, received umber of pulse N by
Following formula is calculated:
In formula, θα, 0.5For half-power antenna bearingt beam angle;ΩαFor antenna bearingt scanning speed;θeFor target elevation;
frFor the pulse recurrence frequency of radar.
Step 22, FFT coherent accumulation: coherent accumulation is done to accumulation umber of pulse using FFT Fast Fourier Transform (FFT), to target
Each of radar return image shifts to addition:
Wherein e-jwtEffect be offset by frequency w caused f (0) within the t short time phase shift, j is imaginary number,
That is:
The invention has the following beneficial effects: completed using TinyXML technology and OpenGL polygon lattice method
The drafting on land and beach boundary, considerably reduces calculation amount and given benchmark is easy to monitor;Pass through optimization SIC interframe accumulation
Two threshold value thresholding T1, T2 parameters more preferably inhibit sea clutter in algorithm, improve the recognizable degree of intertidal zone;Pass through two
Optimize threshold value thresholding and processed echo is generated into picture by OpenGL, was read and worked as using Canny function in OpenCV every three seconds
Preceding picture draws coastline, and the region of the current intertidal zone of real-time display improves the accuracy of identification of intertidal zone.
Specific embodiment
Xia Mianjiehefutuhejuti compare Jia Shishifangshiduibenfamingzuojinyibuxiangxishuoming.
As shown in Figure 1, a kind of intertidal zone monitoring method based on multifunction radar, includes the following steps.
Step 1, in radar return land and beach border map drafting: using TinyXML technology and OpenGL it is more
Side shape lattice method completes the drafting on land and beach border map in radar return.
In step 1, firstly, the initial data to mapping is extracted in initial data using TinyXML technology
LineString element (route) and Polygon element (polygon);Then, it using longitude and latitude sample point interpolation technique, obtains
Take Road different in outermost layer closed boundary region in practical map and region;Finally, utilizing OpenGL polygon point
Technology of formatting completes boundary mapping.
Step 2, enhance radar return image: radar return image being enhanced using FFT coherent accumulation technology.
Wherein, specific method radar return image enhanced using FFT coherent accumulation technology, including walk as follows
It is rapid:
Step 21, calculate accumulation umber of pulse N: the antenna bearingt scanning speed for setting radar is preferentially set as 12r/min, thunder
The pulse recurrence frequency reached is fr, distance resolution is θl(being determined by radar pulsewidth), in azimuth scan radar Half Power Beamwidth
In degree, received umber of pulse N is calculated by following formula:
In formula, θα, 0.5For half-power antenna bearingt beam angle;ΩαFor antenna bearingt scanning speed;θeFor target elevation;
frFor the pulse recurrence frequency of radar.
Step 22, FFT coherent accumulation: coherent accumulation is done to accumulation umber of pulse using FFT Fast Fourier Transform (FFT), to target
Each of radar return image shifts to addition:
Wherein e-jwtEffect be offset by frequency w caused f (O) within the t short time phase shift, j is imaginary number,
That is:
If had if when not containing frequency w in signal f (t), and containing frequency w1:
Step 3, the inhibition of sea clutter.
First two threshold value threshold parameters T1, T2 in SIC interframe accumulation algorithm are optimized, obtain optimization threshold value door
Limit;Then, using the SIC interframe accumulation algorithm with optimization threshold value thresholding to radar return image number enhanced in step 2
According to being handled, to inhibit sea clutter.
The algorithm of two threshold value thresholdings T1, T2 are as follows in traditional SIC interframe correlation accumulation:
Having incoherent condition and the slow target echo interframe of movement completely in interframe using sea clutter is strong phase
The characteristic of pass.Interframe data is accumulated by following formula:
Aij_next=aAij+(1-a)Aij_last
A in formulaijIndicate the echo amplitude of (i, j) a processing unit of present frame, Aij_lastIt is (i, j) of present frame
The history echo amplitude of a processing unit, the i.e. output amplitude of former frame, Aij_nextRefer to that (i, j) of present frame a processing is single
First output amplitude.Effective echo will be therefrom extracted by SI treated data by SC processing.
SIC interframe correlation accumulation needs two threshold Ts 1 and T2, thresholding T1 to act on the video data after SI processing,
Thresholding T2 acts on original video data.If original video data is more than threshold T 2, and SI treated amplitude is more than thresholding
T1 then exports current value, otherwise exports zero, and the calculation formula of threshold value is as follows:
K is the constant shifted by the distribution function of false-alarm probability and video data in formula, and derivation process is existing
Technology,It is the mean value of the noise profile of original video data, δ is standard deviation.
Two threshold value thresholdings are determined by false-alarm probability and noise profile function in traditional SIC interframe correlation accumulation,
And the value of two threshold values is only applicable to the environment of uniform sea clutter in traditional SIC interframe correlation accumulation, and same
Threshold value application condition obtained in region is big, or even can curb moving target algorithm is made to fail.
In practical application scene, sea clutter is often non-homogeneous clutter, and the distribution character of different zones sea clutter is different,
Therefore sea clutter is unpredictable, it is also not possible to identical threshold value be arranged to different regions, so to the specificity analysis of sea clutter
Carry out real-time statistics.
The calculation method of clutter amplitude mean-square value are as follows: first choose any smaller area, it is assumed that sea clutter accords in selection area
Uniform properties are closed, to be currently detected the range statistics sea clutter characteristic for opening up MxN size centered on unit (i, j), such as following table
It is shown.
MN sea clutter statistics of features region
Then, the instantaneous echo amplitude value of selection area is obtained by multifunction radar, finds out the square of current clutter amplitude
Value.
The probability density function of sea clutter correlation generation method are as follows:
In formula, x is sea clutter amplitude instantaneous value, 2 σ2For the mean-square value of clutter amplitude.Assuming that extra large in original video data
Clutter probability submits to normal distribution, so the probability density function of its envelope is rayleigh distributed, it is a kind of one-parameter distribution,
The distribution function of accumulation are as follows:
So the calculation formula of two optimizations threshold value threshold parameter T1, T2 are as follows:
T1=μ1+σ1*K
T2=μ1+f(α)*σ1*K
In formula, K is by false-alarm probability PfaAnd the constant that the distribution function of video data is derived by, false-alarm probability Pfa、
The distribution function of video data and the derivation process of K are the prior art, and which is not described herein again for specific derivation process.
μ1It is the mean value of noise profile in original video data, σ1For the standard deviation of noise profile in original video data,
In, μ1And σ1Calculation formula are as follows:
In formula, 2 σ2For the mean-square value of clutter amplitude,For the mean value of clutter amplitude, x is sea clutter amplitude instantaneous value.
Step 4, the radar return image data of ocean clutter cancellation processing, benefit coastline real-time rendering: will be passed through in step 3
Current radar return image data are generated into picture with OpenGL, were read using Canny function in OpenCV every three seconds current
Picture simultaneously draws coastline, the region of the current intertidal zone of real-time display.
Step 5, the refine of coastline data is being just: for, to the position of coastline real-time rendering, utilizing artificial visual in step 4
The data in interpretation and the practical coastline of in-site measurement are modified, so as to Tideland resources real-time monitoring.
The present invention completes drawing for land and beach boundary using TinyXML technology and OpenGL polygon lattice method
System, considerably reduces calculation amount and given benchmark is easy to monitor;Pass through two threshold value thresholdings in optimization SIC interframe accumulation algorithm
T1, T2 parameter more preferably inhibit sea clutter, improve the recognizable degree of intertidal zone;It will be handled by two optimization threshold value thresholdings
It crosses echo and picture is generated by OpenGL, draw coastline using every three seconds reading current images of Canny function in OpenCV, it is real
When show the region of current intertidal zone, improve the accuracy of identification of intertidal zone.
The preferred embodiment of the present invention has been described above in detail, still, during present invention is not limited to the embodiments described above
Detail a variety of equivalents can be carried out to technical solution of the present invention within the scope of the technical concept of the present invention, this
A little equivalents all belong to the scope of protection of the present invention.