CN107678029A - A kind of rear orientation projection's imaging method based on the average cross-correlation information of random reference - Google Patents
A kind of rear orientation projection's imaging method based on the average cross-correlation information of random reference Download PDFInfo
- Publication number
- CN107678029A CN107678029A CN201710764779.8A CN201710764779A CN107678029A CN 107678029 A CN107678029 A CN 107678029A CN 201710764779 A CN201710764779 A CN 201710764779A CN 107678029 A CN107678029 A CN 107678029A
- Authority
- CN
- China
- Prior art keywords
- mrow
- msub
- cross
- echo
- imaged
- 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.)
- Granted
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO 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/00—Systems 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/88—Radar or analogous systems specially adapted for specific applications
- G01S13/89—Radar or analogous systems specially adapted for specific applications for mapping or imaging
- G01S13/90—Radar or analogous systems specially adapted for specific applications for mapping or imaging using synthetic aperture techniques, e.g. synthetic aperture radar [SAR] techniques
- G01S13/9004—SAR image acquisition techniques
- G01S13/9017—SAR image acquisition techniques with time domain processing of the SAR signals in azimuth
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO 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/00—Systems 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/88—Radar or analogous systems specially adapted for specific applications
- G01S13/89—Radar or analogous systems specially adapted for specific applications for mapping or imaging
- G01S13/90—Radar or analogous systems specially adapted for specific applications for mapping or imaging using synthetic aperture techniques, e.g. synthetic aperture radar [SAR] techniques
- G01S13/904—SAR modes
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO 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/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/02—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
- G01S7/28—Details of pulse systems
- G01S7/2813—Means providing a modification of the radiation pattern for cancelling noise, clutter or interfering signals, e.g. side lobe suppression, side lobe blanking, null-steering arrays
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Radar, Positioning & Navigation (AREA)
- Physics & Mathematics (AREA)
- Computer Networks & Wireless Communication (AREA)
- General Physics & Mathematics (AREA)
- Electromagnetism (AREA)
- Signal Processing (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
A kind of rear orientation projection's imaging method based on the average cross-correlation information of random reference, it belongs to radar image processing technology field, solves the problems, such as that traditional rear orientation projection's imaging method secondary lobe and interference are higher, imaging performance is poor, detailed process is as follows：One, by radar gathered data, by calculating, target to be imaged is involutory into the time delay of aperture location and transmitting pulse signal centre frequency obtains the echo data matrix of point to be imaged；Two, three groups of cross correlation number vectors are calculated using the random reference echo vector of middle reference echo vector and the left and right sides, and calculate the average value of the cross-correlation coefficient of three groups of cross correlation number vectors respectively；Three, the average value of optimal cross-correlation coefficient is selected by setting two-stage threshold value；Four, calculate the range value of point to be imaged using the average value and echo vector paragraph of obtained optimal cross-correlation coefficient；Repeat the above steps one to four, scan all points to be imaged, obtain imaging results.The present invention can be used for radar image process field.
Description
Technical field
The present invention relates to a kind of radar image treatment technology, and in particular to one kind based on the average cross-correlation of random reference
Rear orientation projection's imaging method of information, belongs to radar imaging technology field.
Background technology
Radar imaging technology is that the target in scene carries out two dimension or three-dimensional display directly perceived, is easy to follow-up
Object detection and recognition, to obtain the geometry and physical message of object in scene.
In terms of traditional rear orientation projection's imaging, document《TRM-SAR imaging techniques are studied in ultra-wideband ground-penetrating radar (uw-gpr)》(electronics
University of Science and Technology's journal, 2011-05-30) mainly by the way that time reversal mirror imaging technique and SAR imaging techniques are combined, with
By means of its statistics from average characteristics and space-time matched filtering characteristic, contrasted with time domain rear orientation projection imaging algorithm, Ke Yiti
For the imaging results of higher resolution；Document《The parallel fast back projection algorithm of pixel-oriented》(Institutes Of Technology Of Nanjing's journal,
2014-10-30) mainly solves back-projection algorithm operand greatly and poor real by parallel fast back projection algorithm
Problem；Document《Improvement back-projection algorithm based on dechirp Missile-borne SARs》(Acta Physica Sinica, 2013-11-08) is mainly pin
The problem of and original back-projection algorithm difficult to the compensation of missile borne SAR range migration computationally intensive poor real,
A kind of improvement BP imaging algorithms based on solution line frequency modulation Missile-borne SAR are proposed, target area can be realized while parallel processing
Domain accurately image；Patent《GPR rear orientation projection imaging method based on amplitude weighting》(CN 105974405A) is mainly
Response amplitude construction window function corresponding to the point on hyperbola echo is obtained using edge extracting, the then scattering to each passage rings
Amplitude is answered to be weighted imaging；Patent《A kind of SAR rear orientation projection imaging method of feature based apart from subspace》(CN
104155653A) mainly by the angle that based on apart from history vectors, can be focused on from traditional BP algorithm propose away from
From history vectors matching degree, from the optimal amendment measurement distance history vectors in characteristic distance subspace, the essence to compensating phase is realized
Really estimation, obtains more preferable focusing effect；Patent《A kind of scattering strength for GPR Imaging weights processing method》(CN
102253371A) weight coefficient is sought using the scattering data in imaging region at corresponding time lag curve.
Due to being limited by hardware advances, treatment technology is relatively fixed traditional rear orientation projection's imaging method, is not present
Phase information, it can not go to improve imaging performance using phase information, while traditional rear orientation projection's imaging technique is also by near field
The limitation of environment, transmission signal environment are mostly Gaussian pulse signal, therefore traditional rear orientation projection's imaging method secondary lobe and interference
It is higher.
The content of the invention
The present invention is in order to solve in traditional rear orientation projection's imaging method secondary lobe and disturb higher, imaging performance is poor to ask
A kind of topic, there is provided rear orientation projection's imaging method based on the average cross-correlation information of random reference.
Detailed process of the present invention is as follows：
Step 1, radar gathered data is involutory into the time delay and transmitting arteries and veins of aperture location by calculating target to be imaged
Rush signal center frequency and obtain the echo data matrix of point to be imaged；
Step 2, three groups are calculated using middle reference echo vector paragraph and the random reference echo vector paragraph of the left and right sides
Cross correlation number vector, and the average value of the cross-correlation coefficient of three groups of cross correlation number vectors is calculated respectively；
Step 3, the average value of optimal cross-correlation coefficient is selected by setting two-stage threshold value；
Step 4, point to be imaged is calculated using the average value and echo vector paragraph of obtained optimal cross-correlation coefficient
Range value；
Repeat the above steps 1 to 4, scan all points to be imaged, obtain imaging results.
Further, the step 1 is specially：
K-th of synthetic aperture transmitting antenna T of step 11)kPosition be (xk, Δ y/2,0), k-th of synthetic aperture receives day
Line RkPosition be (xk,-Δ y/2,0), target A to be imaged position is (xA,0,zA)；
Step 12) target A to be imaged relative to k-th of synthetic aperture position delay, τA,kFor
The number of synthetic aperture is Np, target A to be imaged to each synthetic aperture position time delay vector：
The centre frequency of step 13) target A transmission signals to be imaged is f0, equivalent sampling frequency is fs, at each time delay
Centered on position, in the upper and lower echo data for taking S length respectively of the place-centric as pending echo data, this segment signal
Length be：
Wherein L is odd number, and echo data and the pending echo data collectively constitute the dimension echo of L × 1 at delay positions
Vector paragraph, symbolExpression rounds symbol；
Echo vector paragraph is step 14) accordingly
Corresponding echo section is target A to be imaged echo section, s at k-th of synthetic aperture position of transmitting antennaA,k
Echo vectors of the target A to be imaged at k-th of synthetic aperture position of transmitting antenna is represented, if sample range is beyond sampled point
Number scope M, or less than 0, then with 0 come polishing, one L × N on imaging point of generationpThe echo data matrix of dimension
Further, the step 2 is specially：
Step 21)For middle reference echo vector paragraph, then one is randomly choosed on the left of middle reference echo section
Road left side reference echo vector paragraphOne of right side reference echo vector paragraph is randomly choosed on the right side of middle reference echo section
Step 22) is according to the echo vector paragraph that each synthetic aperture position obtains and middle reference echo vector paragraph, left side
The correlation of Control echo vector paragraph and right side reference echo vector paragraph, three groups of cross correlation number vectors are calculatedWithWherein cross-correlation coefficient
ρA,1k、ρA,2k、ρA,3kCalculating use Pearson correlation coefficient：
Wherein Cov (i, j) function representation vector i and vectorial j covariance；
Step 23) calculates the average value of the cross-correlation coefficient in three groups of cross correlation number vectors respectivelyI.e.
Further, the step 3 is specially：
Step 31) judges to select the average value of optimal cross-correlation coefficient using dual threshold, if threshold value is respectively th1With
th2, the average value of described cross-correlation coefficientIn maximum be
Step 32) selects the average value of optimal cross-correlation coefficientFor
Further, the step 4 is specially：
The range value of point to be imaged is calculated using the average value and echo vector paragraph of obtained optimal cross-correlation coefficient
EAMethod be：
The beneficial effects of the invention are as follows：Rear orientation projection's imaging method of the present invention is by calculating echo vector paragraph with referring back to
The cross-correlation coefficient of wave band vector, using cross-correlation coefficient weighted imaging range value, while the present invention is selected using dual threshold thresholding
Optimal average cross correlation coefficient is selected, adds the robustness of average cross correlation coefficient, can be avoided because strong jamming is referring back to
Performance difference caused by wave band vector position, amplitude of rear orientation projection's imaging method in target location in the present invention are protected
Stay, be inhibited in the amplitude of secondary lobe and interference sections, the performance of effective suppressed sidelobes and interference can be played.
The present invention can be widely applied to radar image processing technology field, by suppressed sidelobes and interference, obtain more preferable
Imaging performance.
Brief description of the drawings
Fig. 1 is a kind of stream of rear orientation projection's imaging method based on the average cross-correlation information of random reference of the present invention
Cheng Tu.
Fig. 2 is a kind of rear orientation projection's imaging method based on the average cross-correlation information of random reference of the present invention double
Model of place figure under mode of operation of standing.
Fig. 3 is a kind for the treatment of for rear orientation projection's imaging method based on the average cross-correlation information of random reference of the present invention
Echo data figure of the imageable target in each synthetic aperture opening position.
Fig. 4 is in a kind of rear orientation projection's imaging method based on the average cross-correlation information of random reference of the present invention
Between refer back to the location drawing of wave band, left side reference echo section and right side reference echo section.
Fig. 5 is the imaging results figure of traditional rear orientation projection's imaging method.
Fig. 6 is the imaging results figure using rear orientation projection's imaging method of random reference cross-correlation information.
Fig. 7 be random reference cross-correlation information rear orientation projection's imaging method with traditional rear orientation projection's imaging method vertical
The Profile Correlation figure of opening position.
Fig. 8 is Fig. 7 partial enlarged drawing.
Fig. 9 is rear orientation projection's imaging method and the traditional rear orientation projection's imaging method of random reference cross-correlation information in level
The Profile Correlation figure of opening position.
Figure 10 is Fig. 9 partial enlarged drawing.
Figure 11 is rear orientation projection's imaging method based on the average cross-correlation information of random reference of the present invention in strong jamming ring
The echo data figure of target location under border.
Figure 12 is Figure 11 partial enlarged drawing.
Figure 13 is rear orientation projection's imaging results figure that cross-correlation coefficient average value is sought merely with middle reference echo vector paragraph.
Figure 14 is Figure 13 partial enlarged drawing.
Figure 15 is rear orientation projection's imaging method based on the average cross-correlation information of random reference of the present invention in strong jamming ring
Imaging results figure under border.
Figure 16 is Figure 15 partial enlarged drawing.
In figure：1 represents k-th of synthetic aperture position of transmitting antenna；2 represent k-th of synthetic aperture reception antenna position；3
Represent intermediate solid round dot.
Embodiment
The technical scheme in the embodiment of the present invention is carried out below in conjunction with the accompanying drawing in the embodiment of the present invention clear, complete
Ground describes, it is clear that described embodiment is only part of the embodiment of the present invention, rather than whole embodiments.Based on this
Embodiment in invention, the every other reality that those of ordinary skill in the art are obtained under the premise of creative work is not made
Example is applied, belongs to the scope of protection of the invention.
Embodiment one, with reference to Fig. 1 illustrate present embodiment.It is a kind of based on the average cross-correlation information of random reference
Rear orientation projection's imaging method, detailed process are as follows：
Step 1, radar gathered data is involutory into the time delay and transmitting arteries and veins of aperture location by calculating target to be imaged
Rush signal center frequency and obtain the echo data matrix of point to be imaged；
Step 2, three groups are calculated using middle reference echo vector paragraph and the random reference echo vector paragraph of the left and right sides
Cross correlation number vector, and the average value of the cross-correlation coefficient of three groups of cross correlation number vectors is calculated respectively；
Step 3, the average value of optimal cross-correlation coefficient is selected by setting two-stage threshold value；
Step 4, point to be imaged is calculated using the average value and echo vector paragraph of obtained optimal cross-correlation coefficient
Range value；
Repeat the above steps 1 to 4, scan all points to be imaged, obtain imaging results.
Embodiment two, referring to Fig. 2 and 3 explanation present embodiment.Present embodiment is to embodiment one
Further explanation, step 1 is specially：
K-th of synthetic aperture transmitting antenna T of step 11)kPosition be (xk, Δ y/2,0), k-th of synthetic aperture receives day
Line RkPosition be (xk,-Δ y/2,0), target A to be imaged position is (xA,0,zA)；
Step 12) target A to be imaged relative to k-th of synthetic aperture position delay, τA,kFor
The number of synthetic aperture is Np, target A to be imaged to each synthetic aperture position time delay vector：
The centre frequency of step 13) target A transmission signals to be imaged is f0, equivalent sampling frequency is fs, at each time delay
Centered on position, in the upper and lower echo data for taking S length respectively of the place-centric as pending echo data, this segment signal
Length be：
Wherein L is odd number, and echo data and the pending echo data collectively constitute the dimension echo of L × 1 at delay positions
Vector paragraph, symbolExpression rounds symbol；
Echo vector paragraph is step 14) accordingly
Corresponding echo section is target A to be imaged echo section, s at k-th of synthetic aperture position of transmitting antennaA,k
Echo vectors of the target A to be imaged at k-th of synthetic aperture position of transmitting antenna is represented, if sample range is beyond sampled point
Number scope M, or less than 0, then with 0 come polishing, one L × N on imaging point of generationpThe echo data matrix of dimension
Under dual station mode of operation, rear orientation projection's imaging method of the invention based on the average cross-correlation information of random reference
Model of place as shown in Fig. 2 in x, y and z three dimensions, the triangle in y-axis positive axis direction represents synthetic aperture and connect
Aerial position is received, the triangle of minus half direction of principal axis of y-axis represents synthetic aperture position of transmitting antenna, and adjacent synthetic aperture
The distance between reception antenna position is Δ l, the distance between adjacent synthetic aperture position of transmitting antenna is Δ l.Such as Fig. 3 institutes
Show, for using the target A to be imaged that the method for step 1 obtains each synthetic aperture opening position echo data figure.
Embodiment three, present embodiment are the further explanations to embodiment one, and step 2 is specially：
Step 21)For middle reference echo vector paragraph, then one is randomly choosed on the left of middle reference echo section
Road left side reference echo vector paragraphOne of right side reference echo vector paragraph is randomly choosed on the right side of middle reference echo section
Step 22) is according to the echo vector paragraph that each synthetic aperture position obtains and middle reference echo vector paragraph, left side
The correlation of Control echo vector paragraph and right side reference echo vector paragraph, three groups of cross correlation number vectors are calculatedWithWherein cross-correlation coefficient
ρA,1k、ρA,2k、ρA,3kCalculating use Pearson correlation coefficient：
Wherein Cov (i, j) function representation vector i and vectorial j covariance；
Step 23) calculates the average value of the cross-correlation coefficient in three groups of cross correlation number vectors respectivelyI.e.
In the present embodiment, as shown in figure 4, black circle represents sound of the point target on channel time delay position in figure
Should, and intermediate solid round dot 3 represents the response in middle reference delay positions, in middle reference delay positions arranged on left and right sides
A delay positions, the delay positions randomly selected in arranged on left and right sides passage and middle reference time delay are respectively randomly selected in passage
Position is collectively as with reference to delay positions.
Embodiment four, present embodiment are the further explanations to embodiment one, and step 3 is specially：
Step 31) judges to select the average value of optimal cross-correlation coefficient using dual threshold, if threshold value is respectively th1With
th2, the average value of described cross-correlation coefficientIn maximum be
Step 32) selects the average value of optimal cross-correlation coefficientFor
Embodiment five, present embodiment are the further explanations to embodiment one, and step 4 is specially：
The range value of point to be imaged is calculated using the average value and echo vector paragraph of obtained optimal cross-correlation coefficient
EAMethod be：
With reference to integration secondary lobe than come the quantitative assessment present invention based on the backward of the average cross-correlation information of random reference
Projection imaging method is to secondary lobe and the rejection ability of interference.
Integration secondary lobe ratio is to be used for weighing imaging method to secondary lobe and the performance of AF panel degree in Radar Signal Processing
Index, integration definition of the secondary lobe than ISLR are：
E in formulatotalRepresent the gross energy in entire image imaging results, EmainRepresent target in entire image imaging results
Energy in main lobe.
For the ease of analyzing the integration secondary lobe ratio of imaging method of the present invention and traditional rear orientation projection's imaging method, the present invention adopts
With simple image scene as shown in Figure 4, a point target is only existed in the scene, and signal is preferable carrierfree Gaussian pulse
Signal, the neighbouring N of point targetpThere is DELAY RESPONSE of the target in synthetic aperture position in track data, i.e. echo shaping is identical, only
It is difference in delay positions, so hyp shape occurs.
Assuming that the amplitude of delay positions where point target is B, represented in Fig. 4 with black circle, other positions are with hollow
Round dot represents.In radar detection environment, always with the presence of various noises or interference, it is assumed herein that the amplitude η of other positions
Very little, and obey N (μ, σ2) Gaussian Profile, wherein μ absolute value is much smaller than target amplitude value B.In theory deduction below
It is middle to be expressed as with T as the range value of rear target, it is expressed as the range value as rear secondary lobe and interference with C.The rear orientation projection of the present invention
When imaging method is imaged for point to be imaged, if point A to be imaged is target point, in obtained cross correlation number vector
Each cross-correlation coefficient approximation is all 1, after each cross-correlation coefficient is averagedIf point to be imaged is non-targeted point, obtain
Each cross-correlation coefficient in cross correlation number vector is 1 except Control echo vector paragraph position, the cross-correlation coefficient of other positions
For the value ρ of very littlemin, so just cause target and the contrast of interference to greatly increase.Traditional rear orientation projection's imaging method be by
All N of point to be imagedpIndividual amplitude is directly cumulative, and rear orientation projection's imaging method of the invention is by all N of point to be imagedpIndividual amplitude
Accumulated value multiply optimal average cross correlation coefficient, the amplitude T of target point is respectively in two kinds of algorithms
T in formulaBPRepresent the amplitude of target point after traditional rear orientation projection's imaging method imaging, TSRACBPRepresent the present invention's
The amplitude of target point after rear orientation projection's imaging method imaging based on the average cross-correlation information of random reference,Represent cross-correlation
The average of the cross-correlation coefficient of coefficient vector, at aiming spot
For interference magnitude C, when now first considering to be not present the echo amplitude of target in echo vector paragraph, traditional is backward
The interference magnitude of projection imaging method is approximately Npμ, the rear orientation projection of the invention based on the average cross-correlation information of random reference
Cross-correlation coefficient can be close to 1/N by sum-average arithmetic, in theory average cross correlation coefficient in imaging methodp, can then draw：
CBP=Npμ,CSRACBP=μ
In formula, CBPRepresent the interference magnitude after the non-targeted point imaging of traditional rear orientation projection's imaging method, CSRACBPRepresent this
Interference magnitude after the non-targeted point imaging of rear orientation projection's imaging method based on the average cross-correlation information of random reference of invention, η
It is one and obeys N (μ, σ2) Gaussian Profile very little amplitude.
If consider a target point echo be present in non-targeted echo vector paragraph, then traditional rear orientation projection's imaging method
Interference magnitude C 'BPIt is approximately Npμ+B, rear orientation projection's imaging method of the invention based on the average cross-correlation information of random reference
The average value of cross-correlation coefficient be about 1/Np, therefore interference magnitude CS′RACBPIt is approximately μ+B/Np, can then draw：
C′BP=Npμ+B,C′SRACBP=μ+B/Np
For above-mentioned two situations, because the amplitude of target point is the same order of magnitude, and it is approximately the same, therefore target point
Amplitude is with interference magnitude than being mainly reflected in the difference of interference magnitude.T and C ratio reflects target under different situations and done
The contrast after imaging is disturbed, calculates traditional rear orientation projection's imaging method and the present invention respectively based on the average cross-correlation letter of random reference
The T and C of rear orientation projection's imaging method of breath ratio
For above-mentioned two situations it can be seen that the amplitude after the imaging of target location and the amplitude after the imaging of non-targeted position
There is nearly NpPerformance boost again.If imaging region is divided into M × N number of pixel, and picture shared by interference magnitude in above-mentioned two situations
Prime number mesh is respectively S1、S2, wherein S2＜ S1, number of pixels shared by target is 1, then the ISLR values of two kinds of imaging method results
Respectively：
Wherein, ISLRBPIt is the ISLR, ISLR after traditional rear orientation projection's imaging method imagingSRACBPBe the present invention based on
ISLR after rear orientation projection's imaging method imaging of the average cross-correlation information of machine reference, as can be seen from the above equation, base of the invention
ISLR after rear orientation projection's imaging method imaging of the average cross-correlation information of random reference is than traditional rear orientation projection's imaging method
ISLR after imaging reduces 20lg (Np), that is, illustrate the rear orientation projection based on the average cross-correlation information of random reference of the invention
The performance of imaging method has stronger secondary lobe and AF panel performance compared to more traditional rear orientation projection's imaging method.
The effect of the present invention is described further with reference to emulation experiment.
(1) simulation parameter
Transmission signal centre frequency f0=400MHz, equivalent sampling frequency fs=25GHz, transmission signal are believed for Gaussian pulse
Number, noise is white Gaussian noise, and signal to noise ratio 10dB, the x-axis direction sampling interval is 2cm, and the z-axis direction sampling interval is 6mm,
Sending and receiving antenna spacing is 20cm, and the positions of two targets is (0.8m, 0m, 2m) and (1.4m, 0m, 6m), the letter of target location
Miscellaneous ratio is -10dB.
(1) simulation result
Fig. 5 is the imaging results figure of traditional rear orientation projection's imaging method；Backward throwing based on random reference cross-correlation information
The imaging results of shadow imaging method are as shown in Figure 6；Fig. 7 is the rear orientation projection's imaging method and tradition of random reference cross-correlation information
Profile Correlation figure of rear orientation projection's imaging method in upright position；Accordingly, Fig. 8 is Fig. 7 partial enlarged drawing；Fig. 9 be with
Machine with reference to cross-correlation information rear orientation projection's imaging method and traditional rear orientation projection's imaging method horizontal position section pair
Than figure；Accordingly, Figure 10 is Fig. 9 partial enlarged drawing, from profile it can be seen that the present invention's is averagely mutual based on random reference
For rear orientation projection's imaging method of relevant information compared with traditional back-projection algorithm, secondary lobe and interfering energy are relatively low；Figure 13 is only
Rear orientation projection's imaging results figure of cross-correlation coefficient average value is sought using middle reference echo vector paragraph；Accordingly, Figure 14 is figure
The partial enlarged drawing of 13 imaging results；Figure 15 be the present invention the rear orientation projection based on the average cross-correlation information of random reference into
Imaging results figure of the image space method under strong interference environment；Accordingly, Figure 16 is the partial enlarged drawing of Figure 15 imaging results；From two
Imaging results Figure 14 and Figure 16 of kind method can be seen that to be judged to select optimal cross-correlation coefficient using the method for dual threshold
Average value, strong jamming can be avoided to cause degradation in target location, that is, there is stronger antijamming capability.
Above to a kind of rear orientation projection's imaging method based on the average cross-correlation information of random reference provided by the present invention
It is described in detail, while applies specific embodiment and the principle and embodiment of the present invention are set forth, the above is real
The explanation for applying example is only intended to help the method and its core concept for understanding the present invention；For those of ordinary skill in the art,
According to the thought of the present invention, there will be changes in specific embodiments and applications, in summary, in this specification
Appearance should not be construed as limiting the invention.
Claims (5)
1. a kind of rear orientation projection's imaging method based on the average cross-correlation information of random reference, it is characterised in that including following step
Suddenly：
Step 1, by radar gathered data, by calculating, target to be imaged is involutory into the time delay of aperture location and transmitting pulse is believed
Number centre frequency obtains the echo data matrix of point to be imaged；
Step 2, it is mutual using three groups of middle reference echo vector paragraph and the random reference echo vector paragraph of left and right sides calculating
Relation number vector, and the average value of the cross-correlation coefficient of three groups of cross correlation number vectors is calculated respectively；
Step 3, the average value of optimal cross-correlation coefficient is selected by setting two-stage threshold value；
Step 4, the amplitude of point to be imaged is calculated using the average value and echo vector paragraph of obtained optimal cross-correlation coefficient
Value；
Repeat the above steps 1 to 4, scan all points to be imaged, obtain imaging results.
2. rear orientation projection's imaging method according to claim 1 based on the average cross-correlation information of random reference, its feature
It is, the step 1 is specially：
K-th of synthetic aperture transmitting antenna T of step 11)kPosition be (xk, Δ y/2,0), k-th of synthetic aperture reception antenna Rk
Position be (xk,-Δ y/2,0), target A to be imaged position is (xA,0,zA)；
Step 12) target A to be imaged relative to k-th of synthetic aperture position delay, τA,kFor
<mrow>
<msub>
<mi>&tau;</mi>
<mrow>
<mi>A</mi>
<mo>,</mo>
<mi>k</mi>
</mrow>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<mn>2</mn>
<msqrt>
<mrow>
<msup>
<mrow>
<mo>(</mo>
<msub>
<mi>x</mi>
<mi>A</mi>
</msub>
<mo>-</mo>
<msub>
<mi>x</mi>
<mi>k</mi>
</msub>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>+</mo>
<msup>
<mrow>
<mo>(</mo>
<mi>&Delta;</mi>
<mi>y</mi>
<mo>/</mo>
<mn>2</mn>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>+</mo>
<msubsup>
<mi>z</mi>
<mi>A</mi>
<mn>2</mn>
</msubsup>
</mrow>
</msqrt>
</mrow>
<mi>c</mi>
</mfrac>
</mrow>
The number of synthetic aperture is Np, target A to be imaged to each synthetic aperture position time delay vector：
The centre frequency of step 13) target A transmission signals to be imaged is f0, equivalent sampling frequency is fs, with position at each time delay
Centered on, in the upper and lower echo data for taking S length respectively of the place-centric as pending echo data, the length of this segment signal
Spend and be：
Wherein L is odd number, at delay positions echo data and the pending echo data collectively constitute the dimension echo section of L × 1 to
Amount, symbolExpression rounds symbol；
The echo vector paragraph of the corresponding echo section of step 14) is
Corresponding echo section is target A to be imaged in k-th of synthetic aperture antenna opening position echo section, sA,kExpression is treated into
As target A is in the echo vector of k-th synthetic aperture antenna opening position, if sample range beyond sampling number scope M, or
Less than 0, then with 0 come polishing, one L × N on imaging point of generationpThe echo data matrix of dimension
3. rear orientation projection's imaging method according to claim 1 based on the average cross-correlation information of random reference, its feature
It is, the step 2 is specially：
Step 21)For middle reference echo vector paragraph, then random selection is left together on the left of middle reference echo section
Side Control echo vector paragraphOne of right side reference echo vector paragraph is randomly choosed on the right side of middle reference echo section
Step 22) is according to the echo vector paragraph that each synthetic aperture antenna position obtains and middle reference echo vector paragraph, left side
The correlation of Control echo vector paragraph and right side reference echo vector paragraph, three groups of cross correlation number vectors are calculatedWithWherein cross-correlation coefficient
ρA,1k、ρA,2k、ρA,3kCalculating use Pearson correlation coefficient：
<mrow>
<msub>
<mi>&rho;</mi>
<mrow>
<mi>A</mi>
<mo>,</mo>
<mn>2</mn>
<mi>k</mi>
</mrow>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<mi>C</mi>
<mi>o</mi>
<mi>v</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>x</mi>
<mrow>
<mi>A</mi>
<mo>,</mo>
<mi>k</mi>
</mrow>
</msub>
<mo>,</mo>
<msub>
<mi>x</mi>
<mrow>
<mi>A</mi>
<mo>,</mo>
<msub>
<mi>k</mi>
<mi>L</mi>
</msub>
</mrow>
</msub>
<mo>)</mo>
</mrow>
</mrow>
<msqrt>
<mrow>
<mi>C</mi>
<mi>o</mi>
<mi>v</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>x</mi>
<mrow>
<mi>A</mi>
<mo>,</mo>
<mi>k</mi>
</mrow>
</msub>
<mo>,</mo>
<msub>
<mi>x</mi>
<mrow>
<mi>A</mi>
<mo>,</mo>
<mi>k</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mo>&CenterDot;</mo>
<mi>C</mi>
<mi>o</mi>
<mi>v</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>x</mi>
<mrow>
<mi>A</mi>
<mo>,</mo>
<msub>
<mi>k</mi>
<mi>L</mi>
</msub>
</mrow>
</msub>
<mo>,</mo>
<msub>
<mi>x</mi>
<mrow>
<mi>A</mi>
<mo>,</mo>
<msub>
<mi>k</mi>
<mi>L</mi>
</msub>
</mrow>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</msqrt>
</mfrac>
</mrow>
<mrow>
<msub>
<mi>&rho;</mi>
<mrow>
<mi>A</mi>
<mo>,</mo>
<mn>3</mn>
<mi>k</mi>
</mrow>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<mi>C</mi>
<mi>o</mi>
<mi>v</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>x</mi>
<mrow>
<mi>A</mi>
<mo>,</mo>
<mi>k</mi>
</mrow>
</msub>
<mo>,</mo>
<msub>
<mi>x</mi>
<mrow>
<mi>A</mi>
<mo>,</mo>
<msub>
<mi>k</mi>
<mi>R</mi>
</msub>
</mrow>
</msub>
<mo>)</mo>
</mrow>
</mrow>
<msqrt>
<mrow>
<mi>C</mi>
<mi>o</mi>
<mi>v</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>x</mi>
<mrow>
<mi>A</mi>
<mo>,</mo>
<mi>k</mi>
</mrow>
</msub>
<mo>,</mo>
<msub>
<mi>x</mi>
<mrow>
<mi>A</mi>
<mo>,</mo>
<mi>k</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mo>&CenterDot;</mo>
<mi>C</mi>
<mi>o</mi>
<mi>v</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>x</mi>
<mrow>
<mi>A</mi>
<mo>,</mo>
<msub>
<mi>k</mi>
<mi>R</mi>
</msub>
</mrow>
</msub>
<mo>,</mo>
<msub>
<mi>x</mi>
<mrow>
<mi>A</mi>
<mo>,</mo>
<msub>
<mi>k</mi>
<mi>R</mi>
</msub>
</mrow>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</msqrt>
</mfrac>
</mrow>
Wherein Cov (i, j) function representation vector i and vectorial j covariance；
Step 23) calculates the average value of the cross-correlation coefficient in three groups of cross correlation number vectors respectivelyI.e.
<mrow>
<msub>
<mover>
<mi>&rho;</mi>
<mo>&OverBar;</mo>
</mover>
<mrow>
<mi>A</mi>
<mo>,</mo>
<mi>i</mi>
</mrow>
</msub>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<msub>
<mi>N</mi>
<mi>p</mi>
</msub>
</mfrac>
<mo>&CenterDot;</mo>
<mo>|</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>k</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<msub>
<mi>N</mi>
<mi>p</mi>
</msub>
</munderover>
<msub>
<mi>&rho;</mi>
<mrow>
<mi>A</mi>
<mo>,</mo>
<mi>i</mi>
<mi>k</mi>
</mrow>
</msub>
<mo>|</mo>
<mo>,</mo>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
<mo>,</mo>
<mn>2</mn>
<mo>,</mo>
<mn>3.</mn>
</mrow>
4. rear orientation projection's imaging method according to claim 1 based on the average cross-correlation information of random reference, its feature
It is, the step 3 is specially：
Step 31) judges to select the average value of optimal cross-correlation coefficient using dual threshold, if threshold value is respectively th1And th2, institute
The average value for the cross-correlation coefficient statedIn maximum be
<mrow>
<msub>
<mover>
<mi>P</mi>
<mo>&OverBar;</mo>
</mover>
<mrow>
<mi>A</mi>
<mo>,</mo>
<mi>m</mi>
<mi>a</mi>
<mi>x</mi>
</mrow>
</msub>
<mo>=</mo>
<mi>m</mi>
<mi>a</mi>
<mi>x</mi>
<mo>{</mo>
<msub>
<mover>
<mi>P</mi>
<mo>&OverBar;</mo>
</mover>
<mrow>
<mi>A</mi>
<mo>,</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>,</mo>
<msub>
<mover>
<mi>P</mi>
<mo>&OverBar;</mo>
</mover>
<mrow>
<mi>A</mi>
<mo>,</mo>
<mn>2</mn>
</mrow>
</msub>
<mo>,</mo>
<msub>
<mover>
<mi>P</mi>
<mo>&OverBar;</mo>
</mover>
<mrow>
<mi>A</mi>
<mo>,</mo>
<mn>3</mn>
</mrow>
</msub>
<mo>}</mo>
</mrow>
Step 32) selects the average value of optimal cross-correlation coefficientFor
5. rear orientation projection's imaging method according to claim 1 based on the average cross-correlation information of random reference, its feature
It is, the average value and echo vector paragraph of the optimal cross-correlation coefficient that the specific described utilization of the step 4 obtains calculate
The range value E of point to be imagedAMethod be：
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710764779.8A CN107678029B (en) | 2017-08-30 | 2017-08-30 | Backward projection imaging method based on random reference average cross-correlation information |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710764779.8A CN107678029B (en) | 2017-08-30 | 2017-08-30 | Backward projection imaging method based on random reference average cross-correlation information |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN107678029A true CN107678029A (en) | 2018-02-09 |
| CN107678029B CN107678029B (en) | 2020-08-07 |
Family
ID=61134164
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201710764779.8A Active CN107678029B (en) | 2017-08-30 | 2017-08-30 | Backward projection imaging method based on random reference average cross-correlation information |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN107678029B (en) |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109901164A (en) * | 2019-03-21 | 2019-06-18 | 桂林电子科技大学 | A kind of distributed rear orientation projection's imaging method of synthetic aperture radar |
| CN112904334A (en) * | 2021-01-26 | 2021-06-04 | 中国人民解放军空军工程大学 | Ground penetrating radar back projection fast imaging method based on cross correlation |
| CN114966560A (en) * | 2022-07-29 | 2022-08-30 | 中南大学 | Ground penetrating radar backward projection imaging method and system |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104391295A (en) * | 2014-09-02 | 2015-03-04 | 电子科技大学 | Compressive sensing SAR sparse self-focusing imaging method with optimum image entropy |
| CN106646466A (en) * | 2016-11-04 | 2017-05-10 | 深圳市航天华拓科技有限公司 | Imaging method of weighted back projection algorithm based on principal component analysis |
-
2017
- 2017-08-30 CN CN201710764779.8A patent/CN107678029B/en active Active
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104391295A (en) * | 2014-09-02 | 2015-03-04 | 电子科技大学 | Compressive sensing SAR sparse self-focusing imaging method with optimum image entropy |
| CN106646466A (en) * | 2016-11-04 | 2017-05-10 | 深圳市航天华拓科技有限公司 | Imaging method of weighted back projection algorithm based on principal component analysis |
Non-Patent Citations (4)
| Title |
|---|
| GUANGMIN ZHANG ET AL.: "Novel Imaging Method Based on Cross-correlation Function for Suppressing the Interference of Noise", 《2016 IEEE INTERNATIONAL CONFERENCE ON SIGNAL AND IMAGE PROCESSING》 * |
| LIN ZHOU ET AL.: "A Fast Back-Projection Algorithm Based on Cross Correlation for GPR Imaging", 《IEEE GEOSCIENCE AND REMOTE SENSING LETTERS》 * |
| R. ZETIK ET AL.: "Modified cross-correlation back projection for UWB imaging: numerical examples", 《2005 IEEE INTERNATIONAL CONFERENCE ON ULTRA-WIDEBAND》 * |
| 周琳 等: "基于互相关的探地雷达反向投影成像算法", 《电子与信息学报》 * |
Cited By (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109901164A (en) * | 2019-03-21 | 2019-06-18 | 桂林电子科技大学 | A kind of distributed rear orientation projection's imaging method of synthetic aperture radar |
| CN112904334A (en) * | 2021-01-26 | 2021-06-04 | 中国人民解放军空军工程大学 | Ground penetrating radar back projection fast imaging method based on cross correlation |
| CN112904334B (en) * | 2021-01-26 | 2023-04-25 | 中国人民解放军空军工程大学 | Ground penetrating radar backward projection rapid imaging method based on cross correlation |
| CN114966560A (en) * | 2022-07-29 | 2022-08-30 | 中南大学 | Ground penetrating radar backward projection imaging method and system |
| CN114966560B (en) * | 2022-07-29 | 2022-10-28 | 中南大学 | Ground penetrating radar backward projection imaging method and system |
Also Published As
| Publication number | Publication date |
|---|---|
| CN107678029B (en) | 2020-08-07 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN104851097B (en) | The multichannel SAR GMTI methods aided in based on target shape and shade | |
| CN107561507B (en) | Clutter cancellation method for external radiation source radar | |
| CN107678029A (en) | A kind of rear orientation projection's imaging method based on the average cross-correlation information of random reference | |
| CN109298418B (en) | Radar detection false alarm suppression method and device based on building internal structure characteristics | |
| CN109407055A (en) | The Beamforming Method utilized based on multipath | |
| CN105954745A (en) | Imaging method suitable for through-wall radar multipath phantom inhibition | |
| CN109298417B (en) | Building internal structure detection method and device based on radar signal processing | |
| KR20160012284A (en) | Method and Apparatus for suppressing jammer signals and estimating Angle Of Arrival of original signal using orthogonal of transmitting signal waveform | |
| CN112612005A (en) | Radar main lobe interference resisting method based on deep learning | |
| CN111580099A (en) | Wall clutter suppression method of through-wall imaging radar based on joint entropy | |
| Rasool et al. | Biologically inspired processing of radar waveforms for enhanced delay-Doppler resolution | |
| CN110146881A (en) | A kind of scanning radar super-resolution imaging method based on improvement total variation | |
| CN109471097A (en) | A kind of through-wall radar Signal optimum processing method and device | |
| CN112882016A (en) | Multi-person vital sign detection method based on improved robust adaptive beam forming | |
| CN108416105B (en) | Steady adaptive beam-forming algorithm under pulse and Gaussian noise | |
| Gao et al. | TWR-MCAE: A data augmentation method for through-the-wall radar human motion recognition | |
| CN107783111A (en) | A kind of radar foresight super-resolution imaging method based on maximum entropy criterion | |
| CN115575921B (en) | Pitching-direction-based multichannel multi-interference-base suppression interference suppression method | |
| CN108564962A (en) | Unmanned plane voice signal Enhancement Method based on tetrahedron microphone array | |
| CN116125421B (en) | Array radar multi-echo signal target detection method based on deep learning | |
| CN107132532A (en) | Small target detecting method based on extra large peak restrained and multi-frame joint | |
| CN116930963A (en) | Through-wall imaging method based on wireless communication system | |
| CN110133641A (en) | A kind of through-wall imaging radar target tracking method of dimension self-adaption | |
| Lim et al. | Time-and frequency-domain MIMO FLGPR imaging | |
| CN106371095A (en) | Pulse compression technique-based range imaging method and range imaging system |
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 | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |