CN107678009A - Interference treatment ship carries radar array deformation error compensation and object detection method - Google Patents
Interference treatment ship carries radar array deformation error compensation and object detection method Download PDFInfo
- Publication number
- CN107678009A CN107678009A CN201711172469.3A CN201711172469A CN107678009A CN 107678009 A CN107678009 A CN 107678009A CN 201711172469 A CN201711172469 A CN 201711172469A CN 107678009 A CN107678009 A CN 107678009A
- Authority
- CN
- China
- Prior art keywords
- mrow
- array
- radar
- msub
- matrix
- 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
- 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/41—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00 using analysis of echo signal for target characterisation; Target signature; Target cross-section
Landscapes
- Engineering & Computer Science (AREA)
- Computer Networks & Wireless Communication (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
Present disclose provides a kind of interference treatment ship to carry radar array deformation error compensation and object detection method, including:Step S1, using two, object field scape or moving target carry out real aperture imaging to array echo data over the ground in the same time;Step S2, interference treatment is carried out to two width two dimension complex patterns, removes array deformation error phase;Step S3, the relational expression established between echo-signal fixed phase and complex pattern frequency coefficient vector product;Step S4, relational expression reconstruction image is solved using l1 norms optimization criterion;Step S5, contravariant shift to spatial domain and obtain the two dimensional image after array deformation error compensation;Under imaging pattern over the ground, frequency domain low-pass ripple processing can be also carried out to complex pattern 2-d spectrum after interference, contravariant shifts to the 3-D view that spatial domain obtains atural object scene.Method of disclosure, which obtains, focuses on good imaging results over the ground, can remove target and moves the influence that defocuses brought, improves radar target acquisition performance.
Description
Technical field
This disclosure relates to Radar Imaging Processing and moving target field of detecting, more particularly to one kind are based on interference treatment and pressure
The ship that contracting perceives carries radar conformal array deformation error compensating method and motive target imaging detection method.
Background technology
Stratospheric airship because its can be resident in high-altitude the characteristics of, in fields such as high altitude surveillance, regional early warning and communication services
Have broad application prospects.The huge volume of stratospheric airship is realizes that real aperture imaging provides bar using large-size antennae
(referring to Li Liechen, the ship of Li Daojing, the flat in Huangping based on the sparse compressed sensing of transform domain carries thinned array antenna radar reality hole to part
Footpath is imaged [J] radar journals, 2016,5 (1):109-117).Because dirigible volume is larger, the size of array antenna is also larger.
When platform moves, the deformation of array antenna is difficult to avoid that.Therefore, it is necessary to analyze shadow of the array deformation error to imaging
Ring, and compensation method corresponding to research.
Since CS theories are introduced into radar imagery field, related research work, which achieves, to be much in progress (referring to Zhu
X X and Bamler R.Superresolving SAR tomography for multidimensional imaging
of urban areas:compressive sensing-based TomoSAR inversion[J].SP Magazine,
2014,31(4):51-58).But the research contents of above-mentioned document is mostly the sparse scene in space, such as surface vessel, aerial
Aircraft etc..In practical application, spatial scene is non-sparse mostly.On the other hand, pair under the conditions of needing research transform domain sparse
Ground imaging problem (referring to Li Liechen, continuous scene thinned array SAR three-dimensional imaging [J] electricity of the Li Dao capital based on compressed sensing
Son and information journal, 2014,36 (9):2166-2172), but for ship the research work of conformal array deformation error compensation is carried not
More, it is necessary to further deeply.
Motive target imaging Detection Techniques are a study hotspots in radar imagery field, are all had in military and civilian
Play the role of important.Because moving target RCS is smaller, and movement locus is not known, and its echoing characteristics is compared with static atural object
For complexity, phenomena such as target glint during radar illumination also be present, along with background clutter, system noise, artificial disturbance etc.
Influence so that moving object detection turns into one of the difficult point in radar imagery field.Based on ship carry array radar, research radar into
As application problem of the technology in moving target detection, imaging detection technology is formed, is had to development modern radar technology important
Value.
The content of the invention
(1) technical problems to be solved
In view of above-mentioned technical problem, present disclose provides a kind of interference treatment ship to carry radar array deformation error compensation and mesh
Mark detection method.Under imaging pattern over the ground, methods described removes phase caused by array deformation error by interference treatment and missed
Difference, it is equivalent to realize that array deformation error accurately compensates by frequency-domain sparse image reconstruction, obtain and focus on good imaging knot over the ground
Fruit.Under moving target detection mode, when array echo data is obtained by aperture synthesis or burst mode, methods described
Target can be removed simultaneously and moves the influence that defocuses brought, improve radar target acquisition performance.
(2) technical scheme
According to an aspect of this disclosure, there is provided a kind of interference treatment ship carries radar array deformation error compensation and target
Detection method, including:Step S1, when array has deformation error, the array echo of radar two at different moments is carried using ship
Object field scape or moving target carry out real aperture imaging to data over the ground, obtain two width two dimension complex patterns;Step S2, mould is being imaged over the ground
Under formula or moving target detection mode, interference treatment is carried out to the two width two dimension complex pattern, removes the random first of resolution cell
Phase error caused by beginning phase and array deformation error;Wherein, under moving target detection mode, while target motion is removed
Caused phase error;Step S3, echo-signal-fixed phase as fixed phase, is established using a wherein width complex pattern phase
Relational expression between complex pattern frequency coefficient vector product;Step S4, utilizeNorm optimization criterion, returned to what is established
Relational expression between ripple signal-fixed phase and complex pattern frequency coefficient vector product is solved, to be restored after being interfered
The frequency coefficient vector of image;Step S5, obtained frequency coefficient vector contravariant is shifted into spatial domain, the atural object after being focused on
Scene or moving target two dimensional image, it is equivalent to realize array deformation error compensation.
According to another aspect of the disclosure, there is provided a kind of interference treatment ship carries radar array deformation error compensation and mesh
Detection method is marked, including:Step S1, under imaging pattern over the ground, when array has deformation error, radar two is carried using ship
Object field scape carries out real aperture imaging to array echo data at different moments over the ground, obtains two width two dimension complex patterns;Step S2, to two
Width two dimension complex pattern carry out interference treatment, remove resolution cell random initial phase and array deformation error caused by phase miss
Difference;Step S3, frequency domain low-pass ripple processing is carried out to the complex pattern 2-d spectrum after interference treatment;Step S4, will be low through frequency domain
2-d spectrum contravariant after pass filter processing shifts to spatial domain, the atural object scene two dimensional image after being focused on, equivalent realization battle array
Row deformation error compensation.
In certain embodiments, in the step S1, if dirigible is with speed vaAlong orientation linear uniform motion, ship carries
Radar two is respectively t at different moments1、t2, array echo data interval T=t1-t2;Under imaging pattern over the ground, t1, t2Moment
The oblique distance matrix of corresponding radar and each scattering unit of target is respectively r1、r2;Radar is operated under the conditions of positive side view, and r1For
The nearest oblique distance of radar and target, r1、r2Meet relational expression:
Wherein, T=t1-t2For array echo data interval (or radar pulse repetition period);The then multiple figure of two width two dimension
As α and αrefExpression formula is respectively:
In formula, A represents picture amplitude matrix,The phasing matrix of two width two dimension complex patterns is represented respectively,The Phase Moment determined under imaging pattern over the ground in two width two dimension complex patterns by signal trip path is represented respectively
Battle array,Random initial phase matrix corresponding to two width two dimension complex pattern resolution cells, Δ R, Δ R' difference table are represented respectively
Show that oblique distance caused by array deformation offsets moment matrix in two width two dimension complex patterns,Two two are represented respectively
Tie up array deformation phase error matrix in complex pattern.
In certain embodiments, in the step S2, due to the spy that ship platform has movement velocity slowly and approximation is hovered
Point, ship carry array deformation height correlation corresponding to radar former and later two array echo datas at different moments, and array deformation causes
Oblique distance offset approximately equal, i.e. Δ R ≈ Δ R',Then under imaging pattern over the ground, after interference treatment
The image expression formula of gained atural object scene is:
In certain embodiments, under moving target detection mode, if moving target is transported along orientation speed v uniform rectilinears
Dynamic (the orientation movement velocity typically results in image defocus), t1、t2When inscribe the distance matrix of radar and each scattering unit of target
Respectively R1、R2';Under the conditions of positive side view, R1For radar and the nearest oblique distance of target, then:
When target speed v is smaller, array echo data interval T is smaller, R1> > vT, then formula (4) Taylor's level
Number expands at v=0:
Wherein,Then two width two dimension the complex pattern α and αrefExpression formula is respectively:
In formula, A represents picture amplitude matrix,The phasing matrix of two images is represented respectively,
The phasing matrix determined under moving target detection mode in two width two dimension complex patterns by signal trip path is represented respectively,Random initial phase matrix corresponding to two width two dimension complex pattern resolution cells is represented respectively, and Δ R, Δ R' are represented respectively
Oblique distance caused by array deformation offsets moment matrix in two width two dimension complex patterns,Two width two dimension is represented respectively
Array deformation phase error matrix, v in complex patternaRepresent dirigible along orientation linear uniform motion speed.
In certain embodiments, in the step S2, the image expression formula of gained moving target is after interference treatment:
When target and the distant of radar, speed v are smaller and when time interval T is shorter,Smaller, it is present
Target imaging detection is not influenceed.
In certain embodiments, between the echo-signal-fixed phase and complex pattern frequency coefficient vector product established
Relational expression be:
Wherein, it is respectively t at different moments that ship, which carries radar two,1、t2, s is that ship carries radar in t1Corresponding to moment transmitting pulse
Echo-signal, Φ are calculation matrix, P t2The diagonal matrix that image reference phase corresponding to moment is formed, Ψ become for Fourier
Change base,Vectorial for the frequency coefficient of image to be restored after interference, image is α after interference treatmentnew。
In certain embodiments, the expression formula of the calculation matrix Φ is:
Φ={ pm(Rl-zn)}L×N (9)
In formula, pm(Rl-zn) signal form of m-th of submatrix transmitting in conformal array, R are carried for shiplFor l-th of sampled point
Position, l=1,2 ..., L, L be sampling number;znFor n-th of imaging unit position, n=1,2 ..., N, N be imaging unit
Number.
In certain embodiments, the diagonal matrix of the fixed phase composition is:
In formula,For the phase with reference to corresponding to complex pattern.
In certain embodiments, in the step S4, the optimization criterion is:
In formula, ε is measurement noise, and s is that ship carries radar in t1Echo-signal corresponding to moment transmitting pulse, Φ are measurement square
Battle array, P t2The diagonal matrix that image reference phase corresponding to moment is formed, Ψ is Fourier transformation base,It is extensive to be treated after interference
The frequency coefficient vector of complex pattern.
(3) beneficial effect
It can be seen from the above technical proposal that disclosure interference treatment ship carries radar array deformation error compensation and target is visited
Survey method at least has the advantages that one of them:
(1) radar conformal array deformation error is carried for ship, the array number of echoes of radar two at different moments is carried using ship
According to real aperture imaging is carried out, obtain two width complex patterns and carry out interference treatment, remove the random initial phase of scattering unit, simultaneously
Remove phase error caused by array deformation error;Picture signal Spectrum compression tremendously low frequency section after interference, passes through frequency-domain sparse
Image reconstruction, it equivalent can realize that array deformation error accurately compensates.
(2) existing method focuses on good imaging results over the ground to obtain, usually require that measurement accuracy better than ten/
One wavelength;And high-precision measurement capability, it will usually increase the volume weight and power consumption of array distortion measurement system equipment.This public affairs
The extraction of root, when array distortion measurement system accuracy is wavelength magnitude, you can obtain and focus on good imaging results over the ground, significantly
Reduce to array distortion measurement system accuracy requirement.
(3) ground or moving air target are directed to, when array has deformation error, to front and rear corresponding at different moments two
Width complex pattern carries out interference treatment, by frequency-domain sparse image reconstruction, equivalent can realize array deformation error compensation, remove simultaneously
Phase error caused by target motion, obtains and focuses on good moving target two dimensional image.
Brief description of the drawings
By the way that shown in accompanying drawing, above and other purpose, the feature and advantage of the disclosure will become apparent from.In whole accompanying drawings
Identical reference instruction identical part, does not deliberately draw accompanying drawing, it is preferred that emphasis is show by actual size equal proportion scaling
Go out the purport of the disclosure.
Fig. 1 is to carry radar according to disclosure ship to be imaged over the ground and moving target detection geometrical model schematic diagram.
Fig. 2 is to carry radar array deformation error compensation and target according to ship of the first embodiment of the present disclosure based on interference treatment
The flow chart of detection method.
Fig. 3 is using oblique caused by array deformation error in disclosure ship load radar conformal array deformation error compensating method
Away from error schematic diagram.
Fig. 4 is under imaging pattern over the ground, conventional imaging method imaging results is illustrated when the error of array deformation shown in Fig. 3 be present
Figure.
Fig. 5 is under imaging pattern over the ground, is based on using the first embodiment of the present disclosure when error of array deformation shown in Fig. 3 be present
The ship of interference treatment carries the imaging results schematic diagram of radar array deformation error compensation and object detection method.
Fig. 6 is under moving target pattern, and conventional imaging method imaging results are illustrated when the error of array deformation shown in Fig. 3 be present
Figure.
Fig. 7 is under moving target pattern, when the error of array deformation shown in Fig. 3 be present, using first embodiment of the present disclosure base
The imaging results schematic diagram of radar array deformation error compensation and object detection method is carried in the ship of interference treatment.
Fig. 8 is to carry radar array deformation error compensation and target according to ship of the first embodiment of the present disclosure based on interference treatment
The flow chart of detection method.
Fig. 9 is under imaging pattern over the ground, is based on using the second embodiment of the present disclosure when error of array deformation shown in Fig. 3 be present
The ship of interference treatment carries the imaging results schematic diagram of radar array deformation error compensation and object detection method.
Embodiment
For the purpose, technical scheme and advantage of the disclosure are more clearly understood, below in conjunction with specific embodiment, and reference
Accompanying drawing, the disclosure is further described.
It should be noted that in accompanying drawing or specification description, similar or identical part all uses identical figure number.It is attached
The implementation for not illustrating or describing in figure, it is form known to a person of ordinary skill in the art in art.In addition, though this
Text can provide the demonstration of the parameter comprising particular value, it is to be understood that parameter is worth accordingly without being definitely equal to, but be able to can connect
The error margin received is similar to be worth accordingly in design constraint.The direction term mentioned in embodiment, such as " on ", " under ",
"front", "rear", "left", "right" etc., only it is the direction of refer to the attached drawing.Therefore, the direction term used is for illustrating not to use
To limit the protection domain of the disclosure.
The present disclosure proposes a kind of interference treatment ship to carry radar array deformation error compensation and object detection method.Fig. 1 is
Disclosure ship carries radar and is imaged over the ground and moving target detection geometrical model schematic diagram.Imaging pattern over the ground or moving target into
As under detection mode, removing phase error caused by array deformation error by interference treatment, and pass through frequency-domain sparse image weight
Build, obtain and focus on good imaging over the ground or motive target imaging result, do not reduce imaging resolution, valency is applied with important
Value.
First, first embodiment
Fig. 2 is to carry radar array deformation error compensation and target acquisition side according to first embodiment of the present disclosure interference treatment ship
The flow chart of method.As shown in Fig. 2 ship of the present embodiment based on interference treatment carries radar array deformation error compensation and target acquisition
Method, including:
Step S1, when array has deformation error, the two neighboring array number of echoes at different moments of radar is carried using ship
Real aperture imaging is carried out according to object field scape over the ground or moving target, obtains two width two dimension complex patterns;
If dirigible is with speed v in step S1aAlong orientation linear uniform motion, ship load radar adjacent two moment is respectively t1、
t2;Under imaging pattern over the ground, t1、t2The oblique distance matrix of radar corresponding to moment and each scattering unit of target is respectively r1、r2, then
Wherein, T=t1-t2For array echo data interval (or radar pulse repetition period);It is sharp in the step S1
Carrying the array echo data of radar two at different moments with ship, object field scape carries out real aperture imaging over the ground, and obtained two width two dimension is multiple
Image α and αrefExpression formula is respectively:
Wherein, A represents picture amplitude matrix,The phasing matrix of two images is represented respectively,
The phasing matrix determined under imaging pattern over the ground in two images by signal trip path is represented respectively,Two are represented respectively
Random initial phase matrix corresponding to width image resolution unit, Δ R, Δ R' are represented in two images caused by array deformation respectively
Oblique distance offsets moment matrix,Array deformation phase error matrix in two images is represented respectively.
Under moving target detection mode simultaneously (target and dirigible move), if moving target along orientation speed v at the uniform velocity
Move along a straight line (the orientation movement velocity typically results in image defocus), t1、t2When inscribe radar and each scattering unit of target away from
It is respectively R from matrix1、R2', then:
When target speed v is smaller, array echo data interval T is smaller, R1> > vT, then above formula can use Taylor's level
Number expands into:
Wherein,(it is to represent R for the ease of simplifying herein2' and introduce variable R2);The then step
In rapid S1, moving target is entered using the array echo data of ship load radar two at different moments under moving target detection mode
The real aperture imaging of row, two obtained width two dimension complex pattern α and αrefExpression formula is:
Wherein, A represents picture amplitude matrix,The phasing matrix of two images is represented respectively,
The phasing matrix determined under moving target detection mode in two images by signal trip path is represented respectively,Respectively
Random initial phase matrix corresponding to two images resolution cell is represented, Δ R, Δ R' represent array deformation in two images respectively
Caused oblique distance offsets moment matrix,Array deformation phase error matrix in two images is represented respectively.
Step S2, under imaging pattern over the ground or moving target detection mode, two width two dimension complex patterns are carried out at interference
Reason, remove resolution cell random initial phase and array deformation error caused by phase error, in moving target detection mode
Under, while remove phase error caused by target motion.
Because ship platform has movement velocity slowly and the characteristics of approximate hovering, ship carries radar former and later two battle arrays at different moments
Array deformation height correlation corresponding to row echo data, oblique distance offset approximately equal caused by array deformation, i.e. Δ R ≈ Δs
R',
Then under imaging pattern over the ground, the image expression formula of gained atural object scene is after interference treatment:
Under moving target detection mode, the image expression formula of gained moving target is after interference treatment:
αnewRepresent to carry out complex conjugate multiplication to two width two dimension complex patterns, obtain the new complex pattern of a width, that is, gained after interfering
New complex pattern, it is the result that two width two dimension complex patterns obtained after interference treatment;
When target and the distant of radar, speed v are smaller and when time interval T is shorter,Smaller, it is present
Target imaging detection is not influenceed;Ship carries radar and is operated in positive side view condition in formula, when beam scanning be present in radar, above-mentioned formula
Derive similar, step S1 and S2 still have above-mentioned performance under the conditions of strabismus.
Step S3, echo-signal-fixed phase and complex pattern as fixed phase, are established using a wherein width complex pattern phase
Relational expression between frequency coefficient vector product;Echo-signal-the fixed phase and complex pattern frequency coefficient vector established multiply
Relational expression between product is:
Wherein, s is that ship carries echo-signal corresponding to the transmitting pulse of some moment of radar, and Φ is calculation matrix, and P is reference
The diagonal matrix that phase is formed, Ψ is Fourier transformation base,For the frequency coefficient vector of image to be restored after interference.
Step S4, utilizeNorm optimization criterion, to the echo-signal-fixed phase established and complex pattern frequency domain system
Relational expression between number vector product is solved, the frequency coefficient vector of image to be restored after being interfered;The measurement square
Battle array Φ expression formula be:
Φ={ pm(Rl-zn)}L×N;
Wherein, pm(Rl-zn) signal form of m-th of submatrix transmitting in conformal array, R are carried for shiplFor l-th of sampled point
Position, l=1,2 ..., L, L be sampling number;znFor n-th of imaging unit position, n=1,2 ..., N, N be imaging unit
Number.The diagonal matrix that the fixed phase is formed is:
Wherein,For the phase with reference to corresponding to complex pattern.The optimization criterion is:
Wherein, ε is measurement noise.
Step S5, gained frequency coefficient vector contravariant is shifted into spatial domain, obtains and focuses on good atural object scene or motion
Target two dimensional image, it is equivalent to realize array deformation error compensation.
Fig. 3 is using oblique caused by array deformation error in disclosure ship load radar conformal array deformation error compensating method
Away from error schematic diagram.As shown in figure 3, oblique distance error caused by array deformation is three rank multinomials in emulation, max value of error is
30mm.Simulating scenes size is 150m × 250m (orientation × distance to), includes ground and an elliptical cone, and elliptic cone is high
For 50m, major axis 100m, short axle about 50m.The backscattering coefficient of circular cone is 1, and the backscattering coefficient on ground is 0.3.Ship carries
Each submatrix transmitting linear FM signal in conformal array, displaced phase center number are 88 (by closely spaced array or thinned array production
It is raw) remaining simulation parameter is as shown in table 1.
The simulation parameter of table 1
Fig. 4 is conventional imaging method imaging results schematic diagram when the error of array deformation shown in Fig. 3 be present.As illustrated, work as
When conformal array has deformation error, if being directly imaged using measurement data, image can produce serious blooming effect,
Scene can not be rebuild.
Fig. 5 is under imaging pattern over the ground, is based on using the first embodiment of the present disclosure when error of array deformation shown in Fig. 3 be present
The ship of interference treatment carries the imaging results schematic diagram of radar array deformation error compensation and object detection method;As illustrated, this
Imaging results corresponding to open array deformation compensation method can reflect the feature of image scene.
Fig. 6 is under moving target pattern, and conventional imaging method imaging results are illustrated when the error of array deformation shown in Fig. 3 be present
Figure;If moving target is the cross target being made up of 9 point targets, each point is at intervals of 20m;Moving target is at the uniform velocity straight along orientation
Line moves, speed 50m/s, and remaining simulation parameter is as shown in table 1.As illustrated, due to phase caused by array deformation error
Error so that motive target imaging result has blooming effect.
Fig. 7 is under moving target pattern, when the error of array deformation shown in Fig. 3 be present, using first embodiment of the present disclosure base
The imaging results schematic diagram of radar array deformation error compensation and object detection method is carried in the ship of interference treatment;As shown in fig. 7,
Methods described defocuses phase by removing, and is rebuild by frequency-domain sparse and moving target is focused on again, obtained moving target
Imaging results focus on it is good, can correct display target positional information.
2nd, second embodiment
Fig. 8 is to carry radar array deformation error compensation and target acquisition side according to second embodiment of the present disclosure interference treatment ship
The flow chart of method.Including:
Step S1, under imaging pattern over the ground, when array has deformation error, the two neighboring difference of radar is carried using ship
Object field scape carries out real aperture imaging to the array echo data at moment over the ground, obtains two width two dimension complex patterns;
Step S2, interference treatment is carried out to two width two dimension complex patterns, removes the random initial phase and array of resolution cell
Phase error caused by deformation error;
Step S3, frequency domain low-pass ripple processing is carried out to the complex pattern 2-d spectrum after interference treatment;
Step S4, gained 2-d spectrum contravariant is shifted into spatial domain, obtains and focuses on good atural object scene two dimensional image, etc.
Effect realizes array deformation error compensation.
Fig. 9 is under imaging pattern over the ground, when the error of array deformation shown in Fig. 3 be present, using second embodiment of the present disclosure base
The imaging results schematic diagram of radar array deformation error compensation and object detection method is carried in the ship of interference treatment.
In addition, in second embodiment step S1, S2, S4 respectively with first embodiment over the ground imaging pattern when the step of
S1, S2, S5 detailed process are similar, and here is omitted.
So far, the embodiment of the present disclosure is described in detail combined accompanying drawing.According to above description, art technology
Personnel should have clear to ship of the disclosure based on interference treatment load radar array deformation error compensation and object detection method
Understanding.
It should be noted that in accompanying drawing or specification text, the implementation that does not illustrate or describe is affiliated technology
Form known to a person of ordinary skill in the art, is not described in detail in field.In addition, the above-mentioned definition to each element and method is simultaneously
Various concrete structures, shape or the mode mentioned in embodiment are not limited only to, those of ordinary skill in the art can carry out letter to it
Singly change or replace.
Certainly, according to being actually needed, ship of the disclosure based on interference treatment carries radar array deformation error compensation and target
The step of detection method also includes other, because the innovation of the same disclosure is unrelated, here is omitted.
Similarly, it will be appreciated that in order to simplify the disclosure and help to understand one or more of each open aspect,
Above in the description to the exemplary embodiment of the disclosure, each feature of the disclosure is grouped together into single implementation sometimes
In example, figure or descriptions thereof.However, the method for the disclosure should be construed to reflect following intention:I.e. required guarantor
The disclosure of shield requires features more more than the feature being expressly recited in each claim.It is more precisely, such as following
Claims reflect as, open aspect is all features less than single embodiment disclosed above.Therefore,
Thus the claims for following embodiment are expressly incorporated in the embodiment, wherein each claim is in itself
Separate embodiments all as the disclosure.
Particular embodiments described above, the purpose, technical scheme and beneficial effect of the disclosure are carried out further in detail
Describe in detail bright, should be understood that the specific embodiment that the foregoing is only the disclosure, be not limited to the disclosure, it is all
Within the spirit and principle of the disclosure, any modification, equivalent substitution and improvements done etc., the guarantor of the disclosure should be included in
Within the scope of shield.
Claims (10)
1. a kind of interference treatment ship carries radar array deformation error compensation and object detection method, including:
Step S1, carrying the array echo data of radar two at different moments using ship, object field scape or moving target carry out real hole over the ground
Footpath is imaged, and obtains two width two dimension complex patterns;
Step S2, under imaging pattern over the ground or moving target detection mode, the two width two dimension complex pattern is carried out at interference
Reason, remove resolution cell random initial phase and array deformation error caused by phase error;Wherein, detected in moving target
Under pattern, while remove phase error caused by target motion;
Step S3, echo-signal-fixed phase and complex pattern frequency domain as fixed phase, are established using a wherein width complex pattern phase
Relational expression between coefficient vector product;
Step S4, utilizeNorm optimization criterion, to echo-signal-fixed phase for being established and complex pattern frequency coefficient to
Relational expression between amount product is solved, the frequency coefficient vector of image to be restored after being interfered;
Step S5, obtained frequency coefficient vector contravariant is shifted into spatial domain, atural object scene or moving target after being focused on
Two dimensional image, it is equivalent to realize array deformation error compensation.
2. a kind of interference treatment ship carries radar array deformation error compensation and object detection method, including:
Step S1, under imaging pattern over the ground, the array echo data of radar two at different moments object field scape over the ground is carried using ship
Real aperture imaging is carried out, obtains two width two dimension complex patterns;
Step S2, interference treatment is carried out to two width two dimension complex patterns, removes random initial phase and the array deformation of resolution cell
Phase error caused by error;
Step S3, frequency domain low-pass ripple processing is carried out to the complex pattern 2-d spectrum after interference treatment;
Step S4, the 2-d spectrum contravariant after the processing of frequency domain low-pass ripple is shifted into spatial domain, the ground object field after being focused on
Scape two dimensional image, it is equivalent to realize array deformation error compensation.
3. interference treatment ship according to claim 1 or 2 carries radar array deformation error compensation and object detection method, its
In, in the step S1, if dirigible is with speed vaAlong orientation linear uniform motion, ship carries radar two and distinguished at different moments
For t1、t2, array echo data interval T=t1-t2;Under imaging pattern over the ground, t1、t2Radar corresponding to moment is each scattered with target
The oblique distance matrix for penetrating unit is respectively r1、r2;Radar is operated under the conditions of positive side view, and r1For radar and the nearest oblique distance of target,
r1、r2Meet relational expression:
<mrow>
<msub>
<mi>r</mi>
<mn>2</mn>
</msub>
<mo>=</mo>
<msqrt>
<mrow>
<msup>
<msub>
<mi>r</mi>
<mn>1</mn>
</msub>
<mn>2</mn>
</msup>
<mo>+</mo>
<msup>
<mrow>
<mo>(</mo>
<msub>
<mi>v</mi>
<mi>a</mi>
</msub>
<mi>T</mi>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
</mrow>
</msqrt>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</mrow>
Then two width two dimension the complex pattern α and αrefExpression formula is respectively:
In formula, A represents picture amplitude matrix,The phasing matrix of two width two dimension complex patterns is represented respectively,
The phasing matrix determined under imaging pattern over the ground in two width two dimension complex patterns by signal trip path is represented respectively,Point
Random initial phase matrix corresponding to two width two dimension complex pattern resolution cells is not represented, and Δ R, Δ R' represent that two width two dimension is multiple respectively
Oblique distance caused by array deformation offsets moment matrix in image,Battle array in two width two dimension complex patterns is represented respectively
Row deformation phase error matrix.
4. interference treatment ship according to claim 3 carries radar array deformation error compensation and object detection method, wherein,
In the step S2, Δ R ≈ Δ R',Under imaging pattern over the ground, gained ground object field after interference treatment
The image expression formula of scape is:
5. interference treatment ship according to claim 1 carries radar array deformation error compensation and object detection method, wherein,
Under moving target detection mode, if moving target is along orientation speed v linear uniform motion, t1、t2When inscribe radar and target
The distance matrix of each scattering unit is respectively R1、R2';Under the conditions of positive side view, R1For radar and the nearest oblique distance of target, then:
<mrow>
<msup>
<msub>
<mi>R</mi>
<mn>2</mn>
</msub>
<mo>&prime;</mo>
</msup>
<mo>=</mo>
<msqrt>
<mrow>
<msup>
<msub>
<mi>R</mi>
<mn>1</mn>
</msub>
<mn>2</mn>
</msup>
<mo>+</mo>
<msup>
<mrow>
<mo>(</mo>
<mi>v</mi>
<mi>T</mi>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>+</mo>
<msup>
<mrow>
<mo>(</mo>
<msub>
<mi>v</mi>
<mi>a</mi>
</msub>
<mi>T</mi>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
</mrow>
</msqrt>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>4</mn>
<mo>)</mo>
</mrow>
</mrow>
Work as R1During > > vT, then formula (4) is expanded into Taylor series at v=0:
<mrow>
<msup>
<msub>
<mi>R</mi>
<mn>2</mn>
</msub>
<mo>&prime;</mo>
</msup>
<mo>&ap;</mo>
<msub>
<mi>R</mi>
<mn>2</mn>
</msub>
<mo>+</mo>
<mfrac>
<msup>
<mrow>
<mo>(</mo>
<mi>v</mi>
<mi>T</mi>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mrow>
<mn>2</mn>
<msub>
<mi>R</mi>
<mn>2</mn>
</msub>
</mrow>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>5</mn>
<mo>)</mo>
</mrow>
</mrow>
Wherein,Then two width two dimension the complex pattern α and αrefExpression formula is respectively:
In formula, A represents picture amplitude matrix,The phasing matrix of two images is represented respectively,Point
The phasing matrix that Biao Shi do not determined in two width two dimension complex patterns by signal trip path under moving target detection mode,
Random initial phase matrix corresponding to two width two dimension complex pattern resolution cells is represented respectively, and Δ R, Δ R' represent two width two dimension respectively
Oblique distance caused by array deformation offsets moment matrix in complex pattern,Represent respectively in two width two dimension complex patterns
Array deformation phase error matrix, vaRepresent dirigible along orientation linear uniform motion speed.
6. interference treatment ship according to claim 5 carries radar array deformation error compensation and object detection method, wherein,
In the step S2, the image expression formula of gained moving target is after interference treatment:
7. interference treatment ship according to claim 1 carries radar array deformation error compensation and object detection method, wherein,
Relational expression between the echo-signal-fixed phase and complex pattern frequency coefficient vector product established is:
<mrow>
<mi>s</mi>
<mo>=</mo>
<msub>
<mi>&Phi;P&alpha;</mi>
<mrow>
<mi>n</mi>
<mi>e</mi>
<mi>w</mi>
</mrow>
</msub>
<mo>=</mo>
<mi>&Phi;</mi>
<mi>P</mi>
<mi>&Psi;</mi>
<mover>
<mi>&beta;</mi>
<mo>&RightArrow;</mo>
</mover>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>8</mn>
<mo>)</mo>
</mrow>
</mrow>
Wherein, it is respectively t at different moments that ship, which carries radar two,1、t2, s is that ship carries radar in t1Echo corresponding to moment transmitting pulse
Signal, Φ are calculation matrix, P t2The diagonal matrix that image reference phase corresponding to moment is formed, Ψ is Fourier transformation base,Vectorial for the frequency coefficient of image to be restored after interference, image is α after interference treatmentnew。
8. interference treatment ship according to claim 7 carries radar array deformation error compensation and object detection method, wherein,
The expression formula of the calculation matrix Φ is:
Φ={ pm(Rl-zn)}L×N (9)
In formula, pm(Rl-zn) signal form of m-th of submatrix transmitting in conformal array, R are carried for shiplFor l-th of sampling point position, l
=1,2 ..., L, L be sampling number;znFor n-th of imaging unit position, n=1,2 ..., N, N be imaging unit number.
9. interference treatment ship according to claim 7 carries radar array deformation error compensation and object detection method, wherein,
The diagonal matrix that the fixed phase is formed is:
In formula,For the phase with reference to corresponding to complex pattern.
10. interference treatment ship according to claim 3 carries radar array deformation error compensation and object detection method, its
In, in the step S4, the optimization criterion is:
<mrow>
<mi>m</mi>
<mi>i</mi>
<mi>n</mi>
<mo>|</mo>
<mo>|</mo>
<mover>
<mi>&beta;</mi>
<mo>&RightArrow;</mo>
</mover>
<mo>|</mo>
<msub>
<mo>|</mo>
<mn>1</mn>
</msub>
<mo>,</mo>
<mi>s</mi>
<mo>.</mo>
<mi>t</mi>
<mo>.</mo>
<mo>|</mo>
<mo>|</mo>
<mi>s</mi>
<mo>-</mo>
<mi>&Phi;</mi>
<mi>P</mi>
<mi>&Psi;</mi>
<mover>
<mi>&beta;</mi>
<mo>&RightArrow;</mo>
</mover>
<mo>|</mo>
<msubsup>
<mo>|</mo>
<mn>2</mn>
<mn>2</mn>
</msubsup>
<mo>&le;</mo>
<mi>&epsiv;</mi>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>11</mn>
<mo>)</mo>
</mrow>
</mrow>
In formula, ε is measurement noise, and s is that ship carries radar in t1Echo-signal corresponding to moment transmitting pulse, Φ is calculation matrix, P
For t2The diagonal matrix that image reference phase corresponding to moment is formed, Ψ is Fourier transformation base,For figure to be restored after interference
The frequency coefficient vector of picture.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711172469.3A CN107678009B (en) | 2017-11-22 | 2017-11-22 | Deformation error compensation and target detection method for interference processing of shipborne radar array |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711172469.3A CN107678009B (en) | 2017-11-22 | 2017-11-22 | Deformation error compensation and target detection method for interference processing of shipborne radar array |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107678009A true CN107678009A (en) | 2018-02-09 |
CN107678009B CN107678009B (en) | 2021-04-09 |
Family
ID=61149128
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201711172469.3A Active CN107678009B (en) | 2017-11-22 | 2017-11-22 | Deformation error compensation and target detection method for interference processing of shipborne radar array |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107678009B (en) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112098997A (en) * | 2020-09-18 | 2020-12-18 | 欧必翼太赫兹科技(北京)有限公司 | Three-dimensional holographic imaging security inspection radar image foreign matter detection method |
CN116643290A (en) * | 2023-06-16 | 2023-08-25 | 山西建筑工程集团有限公司 | Metering method and system for double-platform motion compensation of irregular contour |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20030020648A1 (en) * | 2001-07-27 | 2003-01-30 | Fienup James R. | System and method for coherent array aberration sensing |
US20060284762A1 (en) * | 2002-01-08 | 2006-12-21 | Garren David A | Process for mapping multiple-bounce ghosting artifacts from radar imaging data |
CN101369018A (en) * | 2007-08-17 | 2009-02-18 | 电子科技大学 | Satellite machine combined double-base synthetic aperture radar frequency domain imaging method |
US20130082870A1 (en) * | 2010-08-26 | 2013-04-04 | David H. Chambers | Synthetic aperture integration (sai) algorithm for sar imaging |
CN103682677A (en) * | 2013-11-14 | 2014-03-26 | 中国科学院电子学研究所 | Airship radar conformal thinned array antenna and its signal processing method |
CN103698764A (en) * | 2013-12-27 | 2014-04-02 | 中国科学院电子学研究所 | Interferometric synthetic aperture radar imaging method under sparse sampling condition |
WO2015168887A1 (en) * | 2014-05-07 | 2015-11-12 | 华为技术有限公司 | Phased array calibration method and phased array calibration circuit |
CN106680818A (en) * | 2016-12-30 | 2017-05-17 | 中国科学院电子学研究所 | Synthetic aperture radar 3D imaging method based on 2D encoding and frequency domain sparse matrix |
-
2017
- 2017-11-22 CN CN201711172469.3A patent/CN107678009B/en active Active
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20030020648A1 (en) * | 2001-07-27 | 2003-01-30 | Fienup James R. | System and method for coherent array aberration sensing |
US20060284762A1 (en) * | 2002-01-08 | 2006-12-21 | Garren David A | Process for mapping multiple-bounce ghosting artifacts from radar imaging data |
CN101369018A (en) * | 2007-08-17 | 2009-02-18 | 电子科技大学 | Satellite machine combined double-base synthetic aperture radar frequency domain imaging method |
US20130082870A1 (en) * | 2010-08-26 | 2013-04-04 | David H. Chambers | Synthetic aperture integration (sai) algorithm for sar imaging |
CN103682677A (en) * | 2013-11-14 | 2014-03-26 | 中国科学院电子学研究所 | Airship radar conformal thinned array antenna and its signal processing method |
CN103698764A (en) * | 2013-12-27 | 2014-04-02 | 中国科学院电子学研究所 | Interferometric synthetic aperture radar imaging method under sparse sampling condition |
WO2015168887A1 (en) * | 2014-05-07 | 2015-11-12 | 华为技术有限公司 | Phased array calibration method and phased array calibration circuit |
CN106680818A (en) * | 2016-12-30 | 2017-05-17 | 中国科学院电子学研究所 | Synthetic aperture radar 3D imaging method based on 2D encoding and frequency domain sparse matrix |
Non-Patent Citations (4)
Title |
---|
LIECHEN LI ET AL.: "Complex-valued interferometric inverse synthetic aperture radar image compression base on compressed sensing", 《THE JOURNAL OF ENGINEERING》 * |
LIECHEN LI, DAOJING LI: "AIRSHIP SPARSE ARRAY ANTENNA RADAR PERFORMANCE ANALYSIS", 《IGARSS 2013》 * |
周建卫 等: "基于共形稀疏阵列的艇载外辐射源雷达性能分析", 《电子与信息学报》 * |
李烈辰 等: "基于变换域稀疏压缩感知的艇载稀疏阵列天线雷达实孔径成像", 《雷达学报》 * |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112098997A (en) * | 2020-09-18 | 2020-12-18 | 欧必翼太赫兹科技(北京)有限公司 | Three-dimensional holographic imaging security inspection radar image foreign matter detection method |
CN116643290A (en) * | 2023-06-16 | 2023-08-25 | 山西建筑工程集团有限公司 | Metering method and system for double-platform motion compensation of irregular contour |
CN116643290B (en) * | 2023-06-16 | 2024-04-26 | 山西建筑工程集团有限公司 | Metering method and system for double-platform motion compensation of irregular contour |
Also Published As
Publication number | Publication date |
---|---|
CN107678009B (en) | 2021-04-09 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104931966B (en) | A kind of spaceborne video SAR image processing methods based on DCS algorithms | |
CN111797717B (en) | High-speed high-precision SAR image ship detection method | |
CN107462887B (en) | Compressed sensing based wide cut satellite-borne synthetic aperture radar imaging method | |
CN103682677B (en) | A kind of ship load radar conformal thinned array antenna and signal processing method thereof | |
CN102353945B (en) | Three-dimensional position reconstructing method based on ISAR (inverse synthetic aperture radar) image sequence for scattering point | |
CN107561533B (en) | A kind of C-band satellite-borne synthetic aperture radar motive target imaging method | |
CN103472450B (en) | Based on the nonuniform space configuration distributed SAR moving target three-D imaging method of compressed sensing | |
CN105388476B (en) | A kind of chromatography SAR imaging methods based on joint sparse model | |
CN104181531B (en) | A kind of three-dimensional relevance imaging method based on phased-array radar | |
CN105223572B (en) | A kind of positive forward sight Bistatic SAR image processing method based on PFA algorithms | |
CN107576961A (en) | A kind of relatively prime down-sampled sparse imaging method of interval synthetic aperture radar | |
CN107102324A (en) | A kind of close shot microwave imaging method and system | |
CN105242255B (en) | Dual-Channel SAR GMTI methods based on compressed sensing | |
CN107064933A (en) | The method that SAR based on circulation Power estimation chromatographs depth of building | |
CN103630903B (en) | The method of flow field, sea radial velocity is measured based on straight rail interference SAR | |
CN109188432A (en) | A kind of parallel biradical Spotlight SAR Imaging rapid bp imaging method | |
CN112835009A (en) | Data set expansion method based on FMCW radar RAI | |
CN105093225A (en) | Inverse synthetic aperture radar self-focusing imaging method based on double sparse constraints | |
CN105652271A (en) | Super-resolution processing method for augmented Lagrangian real-beam radar angle | |
CN103969628A (en) | Compressive-sensing-based synthetic aperture radar PFA imaging algorithm | |
CN106680818A (en) | Synthetic aperture radar 3D imaging method based on 2D encoding and frequency domain sparse matrix | |
CN107678009A (en) | Interference treatment ship carries radar array deformation error compensation and object detection method | |
CN107607945A (en) | A kind of scanning radar forword-looking imaging method based on spatial embedding mapping | |
CN109814074A (en) | Multiple targets tracking based on image aspects processing | |
CN109188384A (en) | The Electromagnetic Simulation method of extraterrestrial target echo dynamic observation |
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 |