CN107390195A - Two-dimensional encoded synthetic aperture radar three-dimensional imaging array deformation error compensating method - Google Patents
Two-dimensional encoded synthetic aperture radar three-dimensional imaging array deformation error compensating method Download PDFInfo
- Publication number
- CN107390195A CN107390195A CN201710602673.8A CN201710602673A CN107390195A CN 107390195 A CN107390195 A CN 107390195A CN 201710602673 A CN201710602673 A CN 201710602673A CN 107390195 A CN107390195 A CN 107390195A
- Authority
- CN
- China
- Prior art keywords
- dimensional
- signal
- echo
- deformation error
- array deformation
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
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
-
- 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
-
- 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
Abstract
Present disclose provides a kind of two-dimensional encoded synthetic aperture radar three-dimensional imaging array deformation error compensating method, including:Straight rail is carried out to echo-signal to cross rail to spatial modulation with two-dimentional MURA codings;To array deformation error be present in cross rail and along straight rail to slow time-varying under the conditions of, two-way echo-signal after modulation is carried out down to regard three-dimensional imaging;Interference treatment is carried out to two width three-dimensional complex patterns, removes phase error caused by random initial phase and array deformation error;The SAR complex pattern frequency coefficient relational expressions established after the interference of time domain echo-signal;Using l1Norm optimization criterion, is solved to relational expression;Gained frequency coefficient vector contravariant is shifted into spatial domain, obtains target scene 3-D view, it is equivalent to realize array deformation error compensation.The three-dimensional imaging array deformation error compensating method, still obtained when array distortion measurement system accuracy is wavelength magnitude and focus on good 3-D view, reduce array distortion measurement system accuracy requirement.
Description
Technical field
This disclosure relates to Radar Imaging Processing field, more particularly to it is a kind of based on MURA codings, sparse spectrum and compressed sensing
Synthetic aperture radar three-dimensional imaging array deformation error compensating method.
Background technology
Airborne array SAR obtains distance to high-resolution by transmitting/receiving wide-band signal, utilizes antenna platform and target
Relative motion obtain orientation high-resolution, obtain cross rail to height using the array antenna structure for being laid in carrier aircraft cross rail direction
Resolution ratio.Lower under mode of operation, airborne array SAR can avoid shadow effect, obtain the three of complicated landform and urban target
Tie up information.SAR system complexity and carrier aircraft load can be reduced using sparse array structure, according to displaced phase center principle and dilute
Image reconstruction algorithm is dredged, graing lobe and high secondary lobe problem that sparse sampling is brought can be avoided.
At present, the correlative study of array deformation error compensation is predominantly estimated on receiving wave beam to the direction of arrival of Scattering Targets
Error correction is counted, array SAR imaging performances are studied less.Airborne D S AR imaging arrays deformation error compensation is studied at present
Article mainly for one dimensional linear array deformation problems by a small margin, and utilize distributed location and attitude measurement system (Position
And Orientation System, POS) measurement deformation when, its position and array distortion measurement required precision reach ten/
One wavelength magnitude.When array sizes are larger, array deformation error is smart along straight rail to slow time-varying, array distortion measurement systematic survey
, it is necessary to further optimization method when spending inadequate.
The high precision position information of each submatrix can be obtained using distributed POS, and obtains array deformation control information, its
Positional accuracy measurement about 5cm~6cm, ensure the required precision of array deformation error compensation (referring to Li Daojing, Teng Xiumin, Pan Zhou
Great distributed locations and the concept of attitude measurement system and application direction [J] radar journals, 2 (4):400-405).Therefore, may be used
Consider the submatrix positional informations that provide of Combined Treatment distribution POS and carry out curve fitting (referring to Zhou Jianwei, Li Daojing, Tian He,
Ship load external illuminators-based radar performance evaluation [J] electronics and information journal of Pan Jie, the Hu Xuan based on conformal thinned array, 2017,
39(5):1058-1063.).On this basis, the array shape based on the optimization of SAR three-dimensional imagings performance can be carried out to echo data
Become error compensation, good image is focused on to obtain.
The content of the invention
(1) technical problems to be solved
In view of above-mentioned technical problem, present disclose provides a kind of two-dimensional encoded synthetic aperture radar three-dimensional imaging array deformation
Error compensating method, SAR system complexity and carrier aircraft load are significantly reduced using sparse array structure.Gone by interference treatment
It is equivalent to realize that array deformation error is accurately mended by frequency-domain sparse image reconstruction except phase is defocused caused by array deformation error
Repay.Disclosure three-dimensional imaging array deformation error compensating method is when array distortion measurement system accuracy is wavelength magnitude, still
It can obtain and focus on good three-dimensional imaging result, array distortion measurement system accuracy requirement is greatly decreased.
(2) technical scheme
According to an aspect of this disclosure, there is provided a kind of two-dimensional encoded synthetic aperture radar three-dimensional imaging array deformation misses
Poor compensation method, including:
Step S1, SAR is imaged to array three-dimensional to airborne cross rail, under frequency division orthogonal signalling MIMO mode of operation,
Carries out straight rail to the echo-signal in the synthetic aperture time with MURA codings respectively to modulate to-cross rail to two-dimensional space, formation two
Road echo-signal;
Step S2, to array deformation error be present in cross rail and along straight rail to slow time-varying under conditions of, MURA is encoded
Two-way echo-signal after modulation carries out down regarding three-dimensional imaging respectively, obtains two width three-dimensional complex patterns;
Step S3, interference treatment is carried out to the two width three-dimensional complex pattern, removes the random initial phase of resolution cell, together
When remove phase error caused by array deformation error, picture signal three-dimensional frequency spectrum is compressed into low-frequency range after interference treatment;
Step S4, using a wherein width three-dimensional complex pattern phase as fixed phase, introduce compressive sensing theory, distance to
The relational expression established between the SAR complex pattern frequency coefficients after echo-signal-interference;
Step S5, utilizes l1Norm optimization criterion, to the SAR complex pattern frequency domains after echo-signal-interference for being established
Relational expression between coefficient is solved, frequency coefficient of the SAR image to be restored under Fourier transformation base Ψ after being interfered
Vector
Step S6, by gained frequency coefficient vectorContravariant shifts to spatial domain, obtains and focuses on good target scene three-dimensional
Image, it is equivalent to realize array deformation error compensation.
In the disclosure some embodiments, the MURA is encoded to two-dimensional pseudo-random code, including the positive and negative codings of MURA,
Its radix-minus-one complement is that code is negated by turn, and property is identical with code.
In the disclosure some embodiments, in step sl, with the positive and negative codings of two-dimentional MURA in the synthetic aperture time
Echo-signal carries out straight rail and modulated to-cross rail to two-dimensional space, including:By MURA it is positive and negative coding respectively with each distance to
Straight rail corresponding to sampled point carries out dot product to-cross rail to two-dimentional echo-signal.
In the disclosure some embodiments, the distance is to straight rail corresponding to sampled point to-cross rail to two-dimentional echo-signal
Be distance to frequency-domain sparse sampling or it is fully sampled under the conditions of obtain.
In the disclosure some embodiments, in step s 2, array deformation error in the two width three-dimensional complex pattern phase
Caused phase error expresses formula:
In formula,Phase error caused by array deformation error in two images is represented respectively, and Δ R, Δ R ' are respectively
Represent oblique distance offset caused by array deformation in two images.
In the disclosure some embodiments, array deformation occurs during echo wave signal acquisition, cross rail to-straight rail to
Two-dimensional space, the corresponding same signal acquisition process of the positive and negative codings of MURA staggeredly, array deformation corresponding to two width three-dimensional complex patterns
It is highly correlated, its corresponding oblique distance offset approximately equal, i.e. Δ R ≈ Δ R ',
In the disclosure some embodiments, in step s 4, the relational expression is:
Wherein, s be two-dimensional space modulation echo-signal by straight rail to-cross rail to the signal matrix after two-dimensional imaging, H
For sampling matrix, Φ is calculation matrix, and P is the diagonal matrix that fixed phase is formed, and Ψ is Fourier transformation base,After interference
The frequency coefficient vector of SAR image to be restored;
The calculation matrix is:
Wherein, KrFor frequency modulation rate, RlBe l-th of distance to sampling point position, l=1,2 ..., L, L be distance to sampled point
Number;znBe distance to n-th of imaging unit position, n=1,2 ..., N, N be distance to imaging unit number.
In the disclosure some embodiments, the diagonal matrix that the fixed phase is formed is:
In formula,For a wherein width three-dimensional complex pattern phase.
In the disclosure some embodiments, in the step S5, the optimization criterion is:
Wherein, ε is measurement noise, and s is the echo-signal of two-dimensional space modulation after straight rail is to-cross rail to two-dimensional imaging
Signal matrix, H is sampling matrix, and Φ is calculation matrix, and P is the diagonal matrix that fixed phase is formed, and Ψ is Fourier transformation
Base,For the frequency coefficient vector of SAR image to be restored after interference.
In the disclosure some embodiments, the code length of the MURA codings is returning closest within the synthetic aperture time
Prime number from ripple signal cross rail to position number, and to MURA codings or echo-signal in cross rail to, straight rail to carrying out both ends respectively
Zero padding, make MURA coding code lengths identical with echo-signal position number.
(3) beneficial effect
It can be seen from the above technical proposal that the disclosure is based on two-dimensional encoded synthetic aperture radar three-dimensional imaging array shape
Become error compensating method and at least have the advantages that one of them:
(1) for the cross rail along straight rail to time-varying to array deformation error, interference treatment is carried out to two width three-dimensional complex patterns,
The random initial phase of scattering unit is removed, while removes phase error caused by array deformation error;Image letter after interference
Number three-dimensional Spectrum compression tremendously low frequency section, by frequency-domain sparse image reconstruction, equivalent can realize that array deformation error accurately compensates.
(2) disclosure is based on two-dimensional encoded synthetic aperture radar three-dimensional imaging array deformation error compensating method, in battle array
When row distortion measurement system accuracy is wavelength magnitude, obtains and focus on good three-dimensional imaging result, position and array is greatly decreased
The requirement of distortion measurement system accuracy, it is 1/10th wavelength to usually require that measurement accuracy.
(3) disclosure is introduced based on two-dimensional encoded synthetic aperture radar three-dimensional imaging array deformation error compensating method
Frequency domain compressed sensing, effectively removing image secondary lobe influences.
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 based on two-dimensional encoded synthetic aperture radar three-dimensional imaging array deformation error compensating method according to the disclosure
Flow chart;
Fig. 2 is that the disclosure is used for airborne cross rail to regarding three-dimensional imaging model schematic under thinned array antenna SAR;
Fig. 3 is that the MURA of 127 × 127 sizes encodes code pattern schematic diagram;
Fig. 4 is that the MURA of 127 × 127 sizes encodes auto-correlation function schematic diagram;
Fig. 5 is to emulate in data verification to observe schematic diagram of a scenario using disclosure imaging method;
Fig. 6 is using oblique distance caused by cross rail to array deformation error of the disclosure imaging method along orientation slow time-varying
Error Two dimensional Distribution schematic diagram;
Fig. 7 is conventional imaging method imaging results schematic diagram when the error of array deformation shown in Fig. 6 be present;
Fig. 8 uses the disclosure three-dimensional imaging array deformation error compensating method when being and the error of array deformation shown in Fig. 6 be present
Imaging results schematic diagram.
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 based on two-dimensional encoded synthetic aperture radar three-dimensional imaging array deformation error compensating method.Adopt
SAR system complexity and carrier aircraft load can be greatly reduced with sparse array structure.Methods described removes array by interference treatment
Phase is defocused caused by deformation error, error is accurately compensated.Method, which can obtain, focuses on good three-dimensional imaging result, reason
Imaging resolution is not reduced by upper, there is important actual application value.
Flows of the Fig. 1 according to the two-dimensional encoded synthetic aperture radar three-dimensional imaging array deformation error compensating method of the disclosure
Figure.As shown in figure 1, this method includes:
Step S1, to cross rail to thinned array three-dimensional imaging SAR, under frequency division orthogonal signalling MIMO mode of operation,
Carries out straight rail to the echo-signal in the synthetic aperture time to-cross rail to two-dimensional space tune with the positive and negative codings of two-dimentional MURA respectively
System, form two-way echo-signal.
Step S2, in cross rail to array deformation error be present and along straight rail to slow time-varying under conditions of, to MURA just,
Two-way echo-signal after Gray code modulation carries out down regarding respectively three-dimensional imaging, two width three-dimensional complex patterns of acquisition, be denoted as respectively α with
αref, expression formula is:
Wherein A represents the magnitude matrix of three-dimensional complex pattern,The Phase Moment of two width three-dimensional complex patterns is represented respectively
Battle array,The phase determined in two images by trip path is represented respectively,Represent to divide in two images respectively
Random initial phase in unit is distinguished,For phase error caused by array deformation error.
Step S3, interference treatment is carried out to two width three-dimensional complex patterns, removes the random initial phase of resolution cell, go simultaneously
Except phase error caused by array deformation error.Picture signal three-dimensional Spectrum compression tremendously low frequency section, expression formula are after interference:
Step S4, using a wherein width three-dimensional complex pattern phase as fixed phase, introduce compressive sensing theory, distance to
The relational expression established between the SAR complex pattern frequency coefficients after echo-signal-interference, the relational expression are:
Wherein, s be two-dimensional space modulation echo-signal by straight rail to-cross rail to the signal matrix after two-dimensional imaging, H
For sampling matrix, Φ is calculation matrix, and P is the diagonal matrix that fixed phase is formed, and Ψ is Fourier transformation base,After interference
The frequency coefficient vector of SAR image to be restored;
Step S5, utilizeNorm optimization criterion, to the SAR complex pattern frequency domains after echo-signal-interference for being established
Relational expression between coefficient is solved, frequency coefficient of the SAR image to be restored under Fourier transformation base Ψ after being interfered
VectorThe optimization criterion is:
Wherein, ε is measurement noise.
Step S6, by gained frequency coefficient vectorContravariant shifts to spatial domain, obtains and focuses on good target scene three-dimensional
Image, it is equivalent to realize that array deformation error accurately compensates.
In step sl, straight rail is carried out to echo-signal using the positive and negative codings of MURA to modulate to-cross rail to two-dimensional space,
Specially:By the positive and negative codings of MURA respectively with each distance to the corresponding straight rail of sampled point to-cross rail to two-dimentional echo-signal
Carry out dot product.Straight rail corresponding to the distance to sampled point can sample to-cross rail to two-dimentional echo-signal in distance to frequency-domain sparse
Or it is fully sampled under the conditions of obtain.The MURA codings are a kind of two-dimensional pseudo-random codes, and its radix-minus-one complement is that code is negated by turn, property
Matter is identical with code.
In step s 2, phase error expression formula is caused by array deformation error in two width three-dimensional complex pattern phases:
Wherein Δ R, Δ R ' represent oblique distance offset caused by array deformation in two images respectively.Array deformation occurs
During echo wave signal acquisition, adopted in cross rail to-straight rail to two-dimensional space, the corresponding same signal of the positive and negative codings of MURA staggeredly
Collection process, array deformation corresponding to two images is height correlation, corresponding oblique distance offset approximately equal, i.e. Δ R ≈ Δs
R ',Though phase error caused by array deformation error can defocus image α and α ref, pass through interference treatment and frequency
Sparse rebuild in domain can remove its influence, equivalent to realize that array deformation error accurately compensates.
In step s 4, the calculation matrix is:
Wherein, KrFor frequency modulation rate, RlBe l-th of distance to sampling point position, l=1,2 ..., L, L be distance to sampled point
Number;znBe distance to n-th of imaging unit position, n=1,2 ..., N, N be distance to imaging unit number.It is described to refer to phase
Position form diagonal matrix be:
Wherein,For a wherein width three-dimensional complex pattern phase.
The code length of MURA coding is the prime number closest to echo-signal cross rail to position number, and to MURA codings or
Echo-signal, to, straight rail to both ends zero padding is carried out respectively, makes MURA coding code lengths identical with echo-signal position number in cross rail.
Fig. 2 show in the disclosure under airborne cross rail thinned array SAR depending on three-dimensional imaging geometrical system schematic diagram and adopt
The thinned array antenna formed with M submatrix, wherein X, Y, Z axis represent straight rail direction, cross rail direction and elevation direction respectively, carry
Machine flying height is H, flying speed v.Cross rail is to thinned array antenna active redundancy(-ce) linear array, wherein between antenna submatrix minimum
Away from being d/2, array antenna length L=(M-1) × d for d, displaced phase center minimum spacing.As M=13, cross rail is to equivalent
Phase centre location number is 65.
Fig. 3 is that the MURA of 127 × 127 sizes encodes code pattern schematic diagram, and Fig. 4 is that the MURA of 127 × 127 sizes is encoded
Auto-correlation function schematic diagram.It can be seen that MURA encodes dutycycle about 50%, and autocorrelation performance is good.
Fig. 5 is to emulate in data verification to observe schematic diagram of a scenario using disclosure imaging method.To the simulation observation in Fig. 5
Scene imaging, parameter are as shown in table 1.
The simulation parameter of table 1
Fig. 6 is using oblique distance caused by cross rail to array deformation error of the disclosure imaging method along orientation slow time-varying
Error Two dimensional Distribution schematic diagram.Wherein cross rail is distributed to array deformation error along three rank multinomials;Straight rail is to kinematic error along just
Chord curve is distributed.Oblique distance max value of error caused by two-dimentional sample plane array deformation error is 0.05m.
Fig. 7 is conventional imaging method imaging results schematic diagram when the error of array deformation shown in Fig. 6 be present, and Fig. 8 is Fig. 6 be present
Using the imaging results schematic diagram of the disclosure three-dimensional imaging array deformation error compensating method during shown array deformation error.As a result
Show that method of disclosure has carried out accurate compensation to array deformation error, correct three-dimensional reconstruction is carried out to target scene.
So far, the embodiment of the present invention is described in detail combined accompanying drawing.According to above description, art technology
Personnel should to the present invention based on two-dimensional encoded synthetic aperture radar three-dimensional imaging array deformation error compensating method have it is clear
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, the present invention is missed based on two-dimensional encoded synthetic aperture radar three-dimensional imaging array deformation
The step of poor compensation method also includes other, because the innovation with the present invention is unrelated, here is omitted.
In addition, unless specifically described or the step of must sequentially occur, the order of above-mentioned steps, which has no, is limited to above institute
Row, and can change or rearrange according to required design.And above-described embodiment can based on design and reliability consideration, that
This mix and match uses using or with other embodiment mix and match, i.e., the technical characteristic in different embodiments can be with independent assortment
Form more embodiments.
Algorithm and display be not inherently related to any certain computer, virtual system or miscellaneous equipment provided herein.
Various general-purpose systems can also be used together with enlightenment based on this.As described above, required by constructing this kind of system
Structure be obvious.In addition, the disclosure is not also directed to any certain programmed language.It should be understood that it can utilize various
Programming language realizes content of this disclosure described here, and the description done above to language-specific is to disclose this public affairs
The preferred forms opened.
The disclosure can be by means of including the hardware of some different elements and by means of properly programmed computer
Realize.The all parts embodiment of the disclosure can realize with hardware, or to be run on one or more processor
Software module is realized, or is realized with combinations thereof.It will be understood by those of skill in the art that can be in practice using micro-
Processor or digital signal processor (DSP) are some or all in the relevant device according to the embodiment of the present disclosure to realize
The some or all functions of part.The disclosure be also implemented as a part for performing method as described herein or
Whole equipment or program of device (for example, computer program and computer program product).Such journey for realizing the disclosure
Sequence can store on a computer-readable medium, or can have the form of one or more signal.Such signal can
Obtained with being downloaded from internet website, either provide on carrier signal or provided in the form of any other.
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 two-dimensional encoded synthetic aperture radar three-dimensional imaging array deformation error compensating method, including:
Step S1, SAR is imaged to array three-dimensional to airborne cross rail, under frequency division orthogonal signalling MIMO mode of operation, respectively
Carries out straight rail to the echo-signal in the synthetic aperture time with MURA codings to modulate to-cross rail to two-dimensional space, form two-way time
Ripple signal;
Step S2, to array deformation error be present in cross rail and along straight rail to slow time-varying under conditions of, to MURA coded modulations
Two-way echo-signal afterwards carries out down regarding three-dimensional imaging respectively, obtains two width three-dimensional complex patterns;
Step S3, interference treatment is carried out to the two width three-dimensional complex pattern, the random initial phase of resolution cell is removed, goes simultaneously
Except phase error caused by array deformation error, picture signal three-dimensional frequency spectrum is compressed into low-frequency range after interference treatment;
Step S4, using a wherein width three-dimensional complex pattern phase as fixed phase, compressive sensing theory is introduced, in distance to foundation
The relational expression between SAR complex pattern frequency coefficients after echo-signal-interference;
Step S5, utilizes l1Norm optimization criterion, to the SAR complex patterns frequency coefficient after echo-signal-interference for being established it
Between relational expression solved, after being interfered SAR image to be restored under Fourier transformation base Ψ frequency coefficient vector
Step S6, by gained frequency coefficient vectorContravariant shifts to spatial domain, obtains and focuses on good target scene 3-D view,
It is equivalent to realize array deformation error compensation.
2. two-dimensional encoded synthetic aperture radar three-dimensional imaging array deformation error compensating method according to claim 1, its
In, the MURA is encoded to two-dimensional pseudo-random code, including the positive and negative codings of MURA, and its radix-minus-one complement is that code is negated by turn, property
It is identical with code.
3. two-dimensional encoded synthetic aperture radar three-dimensional imaging array deformation error compensating method according to claim 2, its
In, in step sl, with two-dimentional MURA positive and negative codings the echo-signal in the synthetic aperture time is carried out straight rail to-cross rail to
Two-dimensional space is modulated, including:By the positive and negative codings of MURA respectively with each distance to the corresponding straight rail of sampled point to-cross rail to
Two-dimentional echo-signal carries out dot product.
4. two-dimensional encoded synthetic aperture radar three-dimensional imaging array deformation error compensating method according to claim 3, its
In, the distance is to be sampled in distance to frequency-domain sparse or full to-cross rail to two-dimentional echo-signal to straight rail corresponding to sampled point
Obtained under sampling condition.
5. two-dimensional encoded synthetic aperture radar three-dimensional imaging array deformation error compensating method according to claim 1, its
In, in step s 2, phase error expression formula is caused by array deformation error in the two width three-dimensional complex pattern phase:
In formula,Phase error caused by array deformation error in two images is represented respectively, and Δ R, Δ R ' are represented respectively
Oblique distance offset caused by array deformation in two images.
6. two-dimensional encoded synthetic aperture radar three-dimensional imaging array deformation error compensating method according to claim 5, its
In, array deformation occurs during echo wave signal acquisition, in cross rail to-straight rail to two-dimensional space, the positive and negative volumes of MURA staggeredly
The corresponding same signal acquisition process of code, array deformation corresponding to two width three-dimensional complex patterns are highly correlated, and its corresponding oblique distance is inclined
Shifting amount approximately equal, i.e. Δ R ≈ Δ R ',
7. two-dimensional encoded synthetic aperture radar three-dimensional imaging array deformation error compensating method according to claim 1, its
In, in step s 4, the relational expression is:
<mrow>
<mi>s</mi>
<mo>=</mo>
<mi>H</mi>
<mi>&Phi;</mi>
<mi>P</mi>
<mi>&Psi;</mi>
<mover>
<mi>&beta;</mi>
<mo>&RightArrow;</mo>
</mover>
</mrow>
Wherein, s be two-dimensional space modulation echo-signal by straight rail to-cross rail to the signal matrix after two-dimensional imaging, H is adopts
Sample matrix, Φ are calculation matrix, and P is the diagonal matrix that fixed phase is formed, and Ψ is Fourier transformation base,It is extensive to be treated after interference
The frequency coefficient vector of multiple SAR image;
The calculation matrix is:
<mrow>
<mi>&Phi;</mi>
<mo>=</mo>
<msub>
<mrow>
<mo>{</mo>
<msup>
<mi>e</mi>
<mrow>
<msub>
<mi>j&pi;K</mi>
<mi>r</mi>
</msub>
<msup>
<mrow>
<mo>(</mo>
<msub>
<mi>R</mi>
<mi>l</mi>
</msub>
<mo>-</mo>
<msub>
<mi>z</mi>
<mi>n</mi>
</msub>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
</mrow>
</msup>
<mo>}</mo>
</mrow>
<mrow>
<mi>L</mi>
<mo>&times;</mo>
<mi>N</mi>
</mrow>
</msub>
</mrow>
Wherein, KrFor frequency modulation rate, RlBe l-th of distance to sampling point position, l=1,2 ..., L, L be distance to sampling number;
znBe distance to n-th of imaging unit position, n=1,2 ..., N, N be distance to imaging unit number.
8. two-dimensional encoded synthetic aperture radar three-dimensional imaging array deformation error compensating method according to claim 7, its
In, the diagonal matrix that the fixed phase is formed is:
In formula,For a wherein width three-dimensional complex pattern phase.
9. two-dimensional encoded synthetic aperture radar three-dimensional imaging array deformation error compensating method according to claim 1, its
In, in the step S5, the optimization criterion is:
<mfenced open = "" close = "">
<mtable>
<mtr>
<mtd>
<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>
</mrow>
</mtd>
<mtd>
<mrow>
<mi>s</mi>
<mo>.</mo>
<mi>t</mi>
<mo>.</mo>
<mo>|</mo>
<mo>|</mo>
<mi>s</mi>
<mo>-</mo>
<mi>H</mi>
<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>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
Wherein, ε is measurement noise, s be two-dimensional space modulation echo-signal by straight rail to-cross rail to the letter after two-dimensional imaging
Number matrix, H are sampling matrix, and Φ is calculation matrix, and P is the diagonal matrix that fixed phase is formed, and Ψ is Fourier transformation base,
For the frequency coefficient vector of SAR image to be restored after interference.
10. two-dimensional encoded synthetic aperture radar three-dimensional imaging array deformation error compensating method according to claim 1, its
In, the code length of the MURA codings is matter of the echo-signal cross rail closest within the synthetic aperture time to position number
Number, and to MURA codings or echo-signal in cross rail to, straight rail to both ends zero padding is carried out respectively, make MURA coding code lengths and echo
Signal location number is identical.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710602673.8A CN107390195A (en) | 2017-07-21 | 2017-07-21 | Two-dimensional encoded synthetic aperture radar three-dimensional imaging array deformation error compensating method |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710602673.8A CN107390195A (en) | 2017-07-21 | 2017-07-21 | Two-dimensional encoded synthetic aperture radar three-dimensional imaging array deformation error compensating method |
Publications (1)
Publication Number | Publication Date |
---|---|
CN107390195A true CN107390195A (en) | 2017-11-24 |
Family
ID=60337085
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710602673.8A Pending CN107390195A (en) | 2017-07-21 | 2017-07-21 | Two-dimensional encoded synthetic aperture radar three-dimensional imaging array deformation error compensating method |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107390195A (en) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109031300A (en) * | 2018-09-03 | 2018-12-18 | 中科卫星应用德清研究院 | Synthetic aperture radar monitors Dangerous Rock Body deformation method and system |
CN109031222A (en) * | 2018-07-09 | 2018-12-18 | 中国科学院电子学研究所 | It navigated again array synthetic aperture radar three-dimensional imaging kinematic error compensation method |
CN110632615A (en) * | 2019-10-30 | 2019-12-31 | 中国科学院电子学研究所 | Synthetic aperture laser radar three-dimensional imaging method based on sparse aperture |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103682677A (en) * | 2013-11-14 | 2014-03-26 | 中国科学院电子学研究所 | Airship radar conformal thinned array antenna and its signal processing method |
CN106508081B (en) * | 2012-12-27 | 2015-06-10 | 中国科学院电子学研究所 | A kind of spaceborne radar large-scale antenna thinned array is realized and signal processing method |
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-07-21 CN CN201710602673.8A patent/CN107390195A/en active Pending
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106508081B (en) * | 2012-12-27 | 2015-06-10 | 中国科学院电子学研究所 | A kind of spaceborne radar large-scale antenna thinned array is realized and signal processing method |
CN103682677A (en) * | 2013-11-14 | 2014-03-26 | 中国科学院电子学研究所 | Airship radar conformal thinned array antenna and its signal processing method |
CN106680818A (en) * | 2016-12-30 | 2017-05-17 | 中国科学院电子学研究所 | Synthetic aperture radar 3D imaging method based on 2D encoding and frequency domain sparse matrix |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109031222A (en) * | 2018-07-09 | 2018-12-18 | 中国科学院电子学研究所 | It navigated again array synthetic aperture radar three-dimensional imaging kinematic error compensation method |
CN109031300A (en) * | 2018-09-03 | 2018-12-18 | 中科卫星应用德清研究院 | Synthetic aperture radar monitors Dangerous Rock Body deformation method and system |
CN110632615A (en) * | 2019-10-30 | 2019-12-31 | 中国科学院电子学研究所 | Synthetic aperture laser radar three-dimensional imaging method based on sparse aperture |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103901428B (en) | Look side ways high-order nonlinear frequency modulation before Missile-borne SAR sub-aperture and become mark formation method | |
CN102662171B (en) | Synthetic aperture radar (SAR) tomography three-dimensional imaging method | |
CN105759263B (en) | A kind of spaceborne Squint SAR radar imaging method under high-resolution large scene | |
JP6660857B2 (en) | Method of generating image of region of interest using radar system | |
CN102053247B (en) | Phase correction method for three-dimensional imaging of multi-base line synthetic aperture radar | |
CN107390195A (en) | Two-dimensional encoded synthetic aperture radar three-dimensional imaging array deformation error compensating method | |
CN105388476B (en) | A kind of chromatography SAR imaging methods based on joint sparse model | |
WO2008073011A1 (en) | A sar radar system and a method relating thereto | |
CN104316924A (en) | Autofocus motion compensation method of airborne ultra-high resolution SAR (Synthetic Aperture Radar) back projection image | |
JP5679904B2 (en) | Image radar device | |
CN104181531B (en) | A kind of three-dimensional relevance imaging method based on phased-array radar | |
CN106680818B (en) | Based on two-dimensional encoded and frequency-domain sparse array synthetic aperture radar three-dimensional imaging method | |
CN103869316A (en) | Method for super-resolution imaging of foresight array SAR based on sparse representation | |
CN107561533B (en) | A kind of C-band satellite-borne synthetic aperture radar motive target imaging method | |
CN110487241B (en) | Method for extracting elevation control points of building area by satellite laser height measurement | |
CN104122549B (en) | Radar angle super-resolution imaging method based on deconvolution | |
CN104122554B (en) | The attribute scattering center feature extracting method of High Resolution SAR Images target | |
CN103308913A (en) | Foresight SAR ambiguity resolving algorithm of double antennas carried by high-speed aircraft | |
CN105954750A (en) | Strip-map synthetic aperture radar non-sparse scene imaging method based on compressed sensing | |
CN104701637B (en) | Anamorphic array antenna electric performance compensation method based on mechanical-electric coupling Yu least square method | |
CN108008389A (en) | A kind of fast frequency-domain rear orientation projection three-D imaging method based on GPU | |
CN107607948A (en) | The motion compensation of medium earth orbit satellite SAR a kind of and imaging method | |
CN109031222A (en) | It navigated again array synthetic aperture radar three-dimensional imaging kinematic error compensation method | |
CN103605121A (en) | Broadband radar data fusion method based on rapid sparse Bayesian learning algorithm | |
JP2014013180A (en) | Radar processor |
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 | ||
RJ01 | Rejection of invention patent application after publication | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20171124 |