CN110297243A - Synthetic aperture radar chromatographs the phase error compensation method and device in three-dimensional imaging - Google Patents
Synthetic aperture radar chromatographs the phase error compensation method and device in three-dimensional imaging Download PDFInfo
- Publication number
- CN110297243A CN110297243A CN201910665019.0A CN201910665019A CN110297243A CN 110297243 A CN110297243 A CN 110297243A CN 201910665019 A CN201910665019 A CN 201910665019A CN 110297243 A CN110297243 A CN 110297243A
- Authority
- CN
- China
- Prior art keywords
- image
- assisted image
- phase error
- assisted
- error compensation
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
- G01S13/88—Radar or analogous systems specially adapted for specific applications
- G01S13/89—Radar or analogous systems specially adapted for specific applications for mapping or imaging
- G01S13/90—Radar or analogous systems specially adapted for specific applications for mapping or imaging using synthetic aperture techniques, e.g. synthetic aperture radar [SAR] techniques
- G01S13/9021—SAR image post-processing techniques
- G01S13/9023—SAR image post-processing techniques combined with interferometric 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/40—Means for monitoring or calibrating
-
- 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)
- Remote Sensing (AREA)
- Radar, Positioning & Navigation (AREA)
- Physics & Mathematics (AREA)
- Computer Networks & Wireless Communication (AREA)
- General Physics & Mathematics (AREA)
- Electromagnetism (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
The embodiment of the present invention provides the phase error compensation method and device in a kind of chromatography three-dimensional imaging of synthetic aperture radar.Wherein, method includes: to obtain phase compensating factor according to the corresponding baseline parameter of assisted image for each assisted image；Phase error compensation is carried out to assisted image according to phase compensating factor；Wherein, the haplopia complex pattern for the ground object target that image satellite acquires supplemented by assisted image.Phase error compensation method and device in synthetic aperture radar chromatography three-dimensional imaging provided in an embodiment of the present invention, phase compensating factor is obtained according to the corresponding baseline parameter of assisted image, phase error compensation is carried out to assisted image according to phase compensating factor, energy is more acurrate, quickly and easily carries out phase error compensation, compensation effect can be improved, so as to obtain the better synthetic aperture radar chromatography three-dimensional imaging result of effect.
Description
Technical field
The present invention relates to electronic signal process technical fields, chromatograph more particularly, to a kind of synthetic aperture radar three-dimensional
Phase error compensation method and device in imaging.
Background technique
Synthetic aperture radar (Synthetic Aperture Radar, SAR) is a kind of active earth observation technology,
Compared with traditional optical sensor, round-the-clock, round-the-clock observation real-time over the ground is may be implemented in SAR, and has certain earth's surface
Penetration capacity.The development of SAR technology brings the sharply expansion of data volume and calculation amount, and with the increasing of real time imagery demand
By force, more important to the Real-time processing of extensive SAR echo data.
Satellite-borne synthetic aperture radar chromatography (SAR chromatography or chromatography SAR, SAR Tomograghy, TomoSAR) technology is more
One of dimension SAR imaging, is a kind of Three Dimensional Ground information inversion technique to grow up on the basis of interference SAR (InSAR), poor
Layering analysis can also be finally inversed by four-dimensional information.
Compressed sensing is a kind of main stream approach of chromatography SAR imaging.In compressed sensing based chromatography SAR imaging process
In, the randomness of phase error has seriously affected the precision of chromatography SAR three-dimensional imaging.Before carrying out chromatography SAR imaging, it is necessary to
Phase compensation is carried out, third dimension tomography is just can be carried out, is otherwise difficult to verify the correctness of its three-dimensional point cloud.Existing phase is missed
Poor compensation method fails to carry out phase error compensation specifically for SAR tomography, usually uses for reference and is imaged for interference SAR
In multidate differential interferometry technology (such as Permanent scatterers interferometry technology, Small Baseline Subset technology, enhancing space difference skill
Art etc.) carry out phase error compensation.But above-mentioned multidate differential interferometry technology exists, and step is complicated, computing cost is big, efficiency
It is low, and there is preferable compensation effect only for partial picture, without the deficiency of universality.Therefore, existing synthetic aperture
Phase error compensation method in radar chromatography three-dimensional imaging is difficult to obtain correct three-dimensional point cloud, it is difficult to effectively be mended
It repays, there are the deficiencies of compensation effect difference.
Summary of the invention
The embodiment of the present invention provides phase error compensation method and dress in a kind of chromatography three-dimensional imaging of synthetic aperture radar
It sets, to solve or at least be partially solved the defect of prior art compensation effect difference.
In a first aspect, the embodiment of the present invention provides the phase error compensation in a kind of chromatography three-dimensional imaging of synthetic aperture radar
Method, comprising:
For each assisted image, phase compensating factor is obtained according to the corresponding baseline parameter of the assisted image；
Phase error compensation is carried out to the assisted image according to the phase compensating factor；
Wherein, the haplopia complex pattern for the ground object target that image satellite acquires supplemented by the assisted image.
Preferably, described for each assisted image, phase is obtained according to the corresponding baseline parameter of the assisted image and is mended
The specific steps for repaying the factor include:
For each assisted image, the inclination angle, described auxiliary of the corresponding baseline of the assisted image and horizontal plane is obtained
Help the distance of the length of the corresponding baseline of image, the reference satellite observation visual angle on ground, radar wavelength and ground object target to seat
Mark；
According to the length of the corresponding baseline of assisted image baseline corresponding with the inclination angle of horizontal plane, the assisted image
Degree, the reference satellite observation visual angle on ground, radar wavelength and ground object target distance to coordinate, obtain the comprehensive compensation factor.
Preferably, described for each assisted image, phase is obtained according to the corresponding baseline parameter of the assisted image and is mended
The specific steps for repaying the factor include:
For each assisted image, the inclination angle, described auxiliary of the corresponding baseline of the assisted image and horizontal plane is obtained
The length of the corresponding baseline of image, the visual angle on reference satellite observation ground and radar wavelength are helped, it is corresponding according to the assisted image
Baseline baseline corresponding with the inclination angle of horizontal plane, the assisted image length, reference satellite observation ground visual angle and thunder
Up to wavelength, the first compensation factor is obtained；Obtain ground object target distance to coordinate, according to the corresponding baseline of the assisted image
The length of baseline corresponding with the inclination angle of horizontal plane, the assisted image, the visual angle on reference satellite observation ground, radar wavelength and
The distance of ground object target to coordinate, obtain the second compensation factor；
Alternatively, obtaining inclination angle, the institute of the corresponding baseline of the assisted image and horizontal plane for each assisted image
State the distance of the length of the corresponding baseline of assisted image, the reference satellite observation visual angle on ground, radar wavelength and ground object target to
Coordinate；According to the length of the corresponding baseline of assisted image baseline corresponding with the inclination angle of horizontal plane, the assisted image,
Reference satellite observe the distance of the visual angle on ground, radar wavelength and ground object target to coordinate, obtain second compensation factor.
Preferably, the comprehensive compensation factor is
Wherein, j indicates imaginary unit；biIndicate the corresponding baseline of assisted image；θ indicates the view on reference satellite observation ground
Angle；αiIndicate the inclination angle of assisted image corresponding baseline and horizontal plane；R indicate ground object target distance to coordinate；λ indicates thunder
Up to wavelength；I indicates the number of auxiliary image satellite, and 1≤i≤N-1, N-1 indicate the sum of auxiliary image satellite.
Preferably, first compensation factor isSecond compensation factor is
Wherein, j indicates imaginary unit；biIndicate the corresponding baseline of assisted image；θ indicates the view on reference satellite observation ground
Angle；αiIndicate the inclination angle of assisted image corresponding baseline and horizontal plane；R indicate ground object target distance to coordinate；λ indicates thunder
Up to wavelength；I indicates the number of auxiliary image satellite, and 1≤i≤N-1, N-1 indicate the sum of auxiliary image satellite.
Preferably, described that the assisted image also wrap before phase error compensation according to the phase compensating factor
It includes:
The assisted image is pre-compensated for according to reference image；
Wherein, the haplopia complex pattern with reference to the ground object target that image is reference satellite acquisition.
Preferably, the specific steps pre-compensated for according to reference image to the assisted image include:
The factor is pre-compensated for reference to image capturing according to describedIt is assisted according to the precompensation factor pair
Image is pre-compensated for；
Wherein, j indicates imaginary unit；R indicate ground object target distance to coordinate；λ indicates radar wavelength.
Second aspect, the embodiment of the present invention provide the phase error compensation in a kind of chromatography three-dimensional imaging of synthetic aperture radar
Device, comprising:
Because of sub-acquisition module, for obtaining phase according to the corresponding baseline parameter of assisted image for each assisted image
Compensation factor；
Phase compensation block, for carrying out phase error compensation to the assisted image according to the phase compensating factor；
Wherein, the haplopia complex pattern for the ground object target that image satellite acquires supplemented by the assisted image.
The third aspect, the embodiment of the present invention provides a kind of electronic equipment, including memory, processor and is stored in memory
Computer program that is upper and can running on a processor, realizes the various possible realizations such as first aspect when executing described program
Phase error compensation side in mode in the three-dimensional imaging of the chromatography of synthetic aperture radar provided by any possible implementation
The step of method.
Fourth aspect, the embodiment of the present invention provide a kind of non-transient computer readable storage medium, are stored thereon with calculating
Machine program, when which is executed by processor realize as first aspect various possible implementations in it is any can
Synthetic aperture radar provided by the implementation of energy chromatographs the step of phase error compensation method in three-dimensional imaging.
Phase error compensation method and device in synthetic aperture radar chromatography three-dimensional imaging provided in an embodiment of the present invention,
Phase compensating factor is obtained according to the corresponding baseline parameter of assisted image, phase is carried out to assisted image according to phase compensating factor
Error compensation, energy is more acurrate, quickly and easily carries out phase error compensation, compensation effect can be improved, so as to obtain effect more
Good synthetic aperture radar chromatographs three-dimensional imaging result.
Detailed description of the invention
In order to more clearly explain the embodiment of the invention or the technical proposal in the existing technology, to embodiment or will show below
There is attached drawing needed in technical description to be briefly described, it should be apparent that, the accompanying drawings in the following description is this hair
Bright some embodiments for those of ordinary skill in the art without creative efforts, can be with root
Other attached drawings are obtained according to these attached drawings.
Fig. 1 is the phase error compensation side chromatographed in three-dimensional imaging according to synthetic aperture radar provided in an embodiment of the present invention
The flow diagram of method；
Fig. 2 is the geometrical model schematic diagram of TomoSAR system provided in an embodiment of the present invention；
Fig. 3 is that the phase error compensation chromatographed in three-dimensional imaging according to synthetic aperture radar provided in an embodiment of the present invention fills
The structural schematic diagram set；
Fig. 4 is the entity structure schematic diagram according to electronic equipment provided in an embodiment of the present invention.
Specific embodiment
In order to make the object, technical scheme and advantages of the embodiment of the invention clearer, below in conjunction with the embodiment of the present invention
In attached drawing, technical scheme in the embodiment of the invention is clearly and completely described, it is clear that described embodiment is
A part of the embodiment of the present invention, instead of all the embodiments.Based on the embodiments of the present invention, those of ordinary skill in the art
Every other embodiment obtained without creative efforts, shall fall within the protection scope of the present invention.
In order to overcome the above problem of the prior art, the embodiment of the present invention provide a kind of synthetic aperture radar chromatography it is three-dimensional at
As in phase error compensation method and device, inventive concept be do not need according to interference complex pattern extract interferometric phase, and
Phase error compensation is carried out according to interferometric phase, it is only necessary to according to the corresponding baseline parameter of each assisted image, standard can be realized
Really, fast, easily phase error compensation.
Fig. 1 is the phase error compensation side chromatographed in three-dimensional imaging according to synthetic aperture radar provided in an embodiment of the present invention
The flow diagram of method.As shown in Figure 1, this method comprises: step S101, for each assisted image, according to assisted image pair
The baseline parameter answered obtains phase compensating factor.
Wherein, the haplopia complex pattern for the ground object target that image satellite acquires supplemented by assisted image.
Fig. 2 is the geometrical model schematic diagram of TomoSAR system provided in an embodiment of the present invention.It should be noted that this hair
The phase error compensation method in synthetic aperture radar chromatography three-dimensional imaging that bright embodiment provides, can be used for as shown in Figure 2
TomoSAR system.
TomoSAR is a kind of SAR 3-D imaging system, as shown in Fig. 2, s in the system0And si(i=1,2 ..., N-1) be
Along the N satellite of base direction b distribution (N is positive integer), wherein s0For reference satellite, siShadow supplemented by (i=1,2 ..., N-1)
As satellite.When carrying out TomoSAR imaging, pass through s0And siThe haplopia that (i=1,2 ..., N-1) acquires ground object target respectively is schemed again
Picture, siThe haplopia complex pattern of (i=1,2 ..., N-1) acquisition respectively with s0The haplopia complex pattern of acquisition is registrated, and N-1 width is obtained
Interfere complex pattern.siThe haplopia complex pattern of (i=1,2 ..., N-1) acquisition is known as assisted image, s0The haplopia complex pattern of acquisition claims
For with reference to image.Each width assisted image respectively corresponds an auxiliary image satellite, and a corresponding width interferes complex pattern.
B is the baseline total length of TomoSAR system, is equal to auxiliary image satellite si(i=1,2 ..., N-1) along base direction b
The difference of the maxima and minima of coordinate.θ is reference satellite s0Observe the visual angle on ground；αiIndicate baseline biWith inclining for horizontal plane
Angle, baseline biFor i-th auxiliary image satellite si(i=1,2 ..., N-1) and reference satellite s0Between the baseline that is formed, length is
bi；S indicates the vertical coordinate of oblique distance, and oblique distance is vertical parallel with base direction；R be distance to coordinate, indicate reference satellite s0To certain
The distance of one range-azimuth unit；rkIndicate vertical upper k-th of scattering point σ of oblique distancekTo reference satellite s0Distance；rikIndicate oblique
Away from vertical upper k-th of scattering point σkTo i-th auxiliary image satellite siThe distance between (i=1,2 ..., N-1).
The distance of ground object target to coordinate, between the scattering point range-azimuth unit fallen in and reference satellite away from
From.
The mathematical model of the haplopia complex pattern of satellite-borne SAR two dimensional image is
Wherein, s (x, r) indicates two-dimentional haplopia complex pattern；The coordinate of x expression orientation；R indicate distance to coordinate；σ
(x, r) indicates the backscattering coefficient of range-azimuth unit (x, r)；RaAnd RrRespectively indicate orientation and distance to resolution
Rate；λ indicates radar wavelength.
All scattering points of the Target scalar on oblique distance is vertical can all be fallen in same range-azimuth unit, by above-mentioned two
The mathematical model for tieing up the haplopia complex pattern of image is generalized to three-dimensional, then i-th auxiliary image satellite si(i=1,2 ..., N-1)
The three-dimensional mathematical model of SAR image can be expressed as
Wherein, si(x, r, s) indicates the assisted image of i-th auxiliary image satellite acquisition；Indicate convolution；J indicates imaginary number list
Position；σk(x, r, s) indicates the backscattering coefficient of k-th of the scattering point fallen in same range-azimuth unit；rik(r, s) i.e.
rik, indicate vertical upper k-th of scattering point σ of oblique distancekTo i-th auxiliary image satellite siThe distance between (i=1,2 ..., N-1).
For reduced model, ignore conventional part, then the three-dimensional mathematical model of above-mentioned SAR image can be abbreviated as
Practical geometrical relationship, has shown according to fig. 2
Due to s and biFar smaller than r, therefore, in two-dimensional imaging, above formula is in (s=0, bi=0) first order Taylor is used at
Expansion.It, need to be to above formula (s=0, b but when carrying out three-dimensional imagingi=0) the second Taylor series are carried out at, and above formula is indicated
At using oblique distance is vertical and the length of baseline is the function of independent variable, have
Further, it is possible to which it is as follows to acquire each local derviation
Each local derviation is substituted into, the mathematical equivalent model approximate expression of available distance is
It further arranges, has
According to the mathematical equivalent model approximate expression of distance, the three-dimensional mathematical model of SAR can be further simplified as
It enablesThen have
It can further be write as when the third dimension is imaged with the related part third dimension space variable s
It further arranges and simplifies, it is available
Wherein,
And the chromatography SAR three-dimensional imaging of various methods, it is all based on the theoretical model being shown below and is imaged:
Wherein, giThe corresponding SAR 3-dimensional image of the auxiliary image satellite of (x, r, s) expression i-th；σk(x, r, s) expression falls in same
The backscattering coefficient of k-th of scattering point in one range-azimuth unit；fiFor spectrum component, baseline b is indicatedi(i.e. i-th of base
Line) it is formed by frequency.
It follows that on vertical, different baseline biCorresponding SAR image can be regarded as on vertical
Atural objectIn frequency fiThe discrete Fourier transform result at place.
But apart from vertical s, different baseline biThere is also spatial offsets for corresponding SAR imageIt wants
When carrying out the compressed sensing imaging of the third dimension using various algorithms, it is necessary to first eliminate the phase as brought by this spatial offset
Difference.
Therefore, phase error includes at least level land phaseAnd spatial offset
Bring error
Due toTherefore
In general, above-mentioned phase error can be removed according to the method for flat earth, but need the frequency by interferometric phase
Spectrum carries out the very high level land phase estimation of precision, it is desirable that quite accurately can just eliminate space to the spectrum estimation of interferometric phase
Phase difference brought by deviating, this is in real data processing, since data volume is big, level land frequency spectrum brought by spatial offset
Deviation is again very small, accurately to estimate this level land frequency spectrum deviation, need just obtain in the sizable multiple of time domain interpolation, this
Data volume can be greatly increased, processing speed is influenced.Meanwhile this method needs to carry out by each assisted image and with reference to image
Interference processing, process is cumbersome, takes a long time, not easy to operate, inefficiency, it is difficult to realize.Also, usual SAR chromatography it is three-dimensional at
Image space method is also without discovery spatial offset bring phase deviation, more without carrying out phase for spatial offset bring phase deviation
Position error compensation, thus compensation effect is poor.
Therefore, the embodiment of the present invention directlys adopt the method for carrying out phase compensation based on baseline parameter, for each auxiliary
Accurately, fast and easily phase error compensation can be realized according to the corresponding baseline parameter of the assisted image in image.
Assisted image is by auxiliary image satellite si(i=1,2 ..., N-1) is acquired, then the corresponding baseline parameter of the assisted image
Including baseline biWith the inclination alpha of horizontal planei, baseline biLength, reference satellite s0Observe the view angle theta on ground, radar wavelength λ and
Ground object target distance to coordinate r.
The assisted image of image satellite acquisition auxiliary for i-th, level land phaseThe space and
OffsetBring errorIt is related with above-mentioned baseline parameter, therefore, according to upper
Baseline parameter is stated, the phase compensating factor for carrying out phase error compensation to the assisted image can be obtained.
Step S102, phase error compensation is carried out to assisted image according to phase compensating factor.
Specifically, after obtaining phase compensating factor, phase compensating factor can be multiplied with assisted image, so as to
Remove level land phaseAnd spatial offsetBring error As spatial offset bring phase error.
It should be noted that only effectively eliminating spatial offset bring phase error, effective phase just can be carried out
Error compensation can just obtain corresponding compressed sensing result.
The embodiment of the present invention obtains phase compensating factor according to the corresponding baseline parameter of assisted image, according to phase compensation because
Son carries out phase error compensation to assisted image, and energy is more acurrate, quickly and easily carries out phase error compensation, can improve compensation
Effect, so as to obtain the better synthetic aperture radar chromatography three-dimensional imaging result of effect.
Content based on the various embodiments described above obtains each assisted image according to the corresponding baseline parameter of assisted image
The specific steps for taking phase compensating factor include: to obtain the corresponding baseline of assisted image and horizontal plane for each assisted image
Inclination angle, the length of the corresponding baseline of assisted image, the reference satellite observation visual angle on ground, radar wavelength and ground object target away from
The coordinate of descriscent.
Specifically, directly once phase error compensation can be finished according to baseline parameter.
According to assisted image and image can be referred to, obtain the corresponding baseline parameter of the assisted image, including the auxiliary shadow
As the length of corresponding baseline baseline corresponding with the inclination angle of horizontal plane, the assisted image, reference satellite observe ground visual angle,
The distance of radar wavelength and ground object target to coordinate.
According to the length of the corresponding baseline of assisted image baseline corresponding with the inclination angle of horizontal plane, assisted image, with reference to defending
Star observe the distance of the visual angle on ground, radar wavelength and ground object target to coordinate, obtain the comprehensive compensation factor.
Specifically, according to the corresponding baseline of assisted image baseline corresponding with the inclination angle of horizontal plane, the assisted image
Length, the reference satellite observation visual angle on ground, radar wavelength and ground object target distance to coordinate, available comprehensive compensation
The factor.
The comprehensive compensation factor, for being compensated simultaneously to level land phase and spatial offset bring phase error.
It is specifically included correspondingly, carrying out phase error compensation to assisted image according to phase compensating factor: by comprehensive compensation
The factor is multiplied with assisted image, can finish to level land phase and spatial offset bring phase error compensation.
The embodiment of the present invention is by directly, once missing level land phase and spatial offset bring phase according to baseline parameter
Difference compensates, and compensation process is extremely simple, can greatly improve compensation efficiency, can greatly save opening for the hardware resources such as computer
Pin.
Content based on the various embodiments described above obtains each assisted image according to the corresponding baseline parameter of assisted image
The specific steps for taking phase compensating factor include: to obtain the corresponding baseline of assisted image and horizontal plane for each assisted image
Inclination angle, the corresponding baseline of assisted image length, reference satellite observation ground visual angle and radar wavelength, according to assisted image
The length of corresponding baseline baseline corresponding with the inclination angle of horizontal plane, assisted image, the visual angle on reference satellite observation ground and thunder
Up to wavelength, the first compensation factor is obtained；Obtain ground object target distance to coordinate, according to the corresponding baseline of assisted image and water
The inclination angle of plane, the length of the corresponding baseline of assisted image, reference satellite observation ground visual angle, radar wavelength and ground object target
Distance to coordinate, obtain the second compensation factor.
Specifically, the first compensation factor and second can be obtained respectively and is mended according to the corresponding baseline parameter of the assisted image
Repay the factor.
First compensation factor, for being compensated to level land phase.
Second compensation factor, for being compensated to spatial offset bring phase error.
Correspondingly, phase error compensation is carried out to assisted image according to phase compensating factor to specifically include: first mending first
It repays the factor to be multiplied with assisted image, level land phase compensation is finished, then multiplied result is multiplied with the second compensation factor, to space
Offset bring phase error compensation finishes.
For each assisted image, the specific steps of phase compensating factor are obtained according to the corresponding baseline parameter of assisted image
It include: that the corresponding baseline of assisted image baseline corresponding with the inclination angle of horizontal plane, assisted image is obtained for each assisted image
Length, the reference satellite observation visual angle on ground, radar wavelength and ground object target distance to coordinate；According to assisted image pair
The length of the baseline answered baseline corresponding with the inclination angle of horizontal plane, assisted image, the visual angle on reference satellite observation ground, radar wave
Long and ground object target distance to coordinate, obtain the second compensation factor.
It, can also be only according to the inclination angle of the corresponding baseline of assisted image and horizontal plane, auxiliary shadow for each assisted image
As the length of corresponding baseline, reference satellite observe the distance of the visual angle on ground, radar wavelength and ground object target to coordinate, obtain
Take the second compensation factor.
Correspondingly, it carries out phase error compensation to assisted image according to phase compensating factor to specifically include: according to general essence
The method of the removal elliptical earth phase of degree first compensates level land phase；It is carried out according to result of second compensation factor to removal elliptical earth phase
Phase error compensation.
The embodiment of the present invention is carried out according to baseline parameter to spatial offset bring by first carrying out the removal of level land phase
The compensation of phase error, without estimating level land frequency spectrum in time domain interpolation, compensation process is simple, can improve compensation efficiency, can save
The expense of the hardware resources such as computer.
Content based on the various embodiments described above, the comprehensive compensation factor are
Wherein, j indicates imaginary unit；biIndicate the corresponding baseline of assisted image；θ indicates the view on reference satellite observation ground
Angle；αiIndicate the inclination angle of assisted image corresponding baseline and horizontal plane；R indicate ground object target distance to coordinate；λ indicates thunder
Up to wavelength；I indicates the number of auxiliary image satellite, and 1≤i≤N-1, N-1 indicate the sum of auxiliary image satellite.
Content based on the various embodiments described above, the first compensation factor areSecond compensation factor is
Wherein, j indicates imaginary unit；biIndicate the corresponding baseline of assisted image；θ indicates the view on reference satellite observation ground
Angle；αiIndicate the inclination angle of assisted image corresponding baseline and horizontal plane；R indicate ground object target distance to coordinate；λ indicates thunder
Up to wavelength；I indicates the number of auxiliary image satellite, and 1≤i≤N-1, N-1 indicate the sum of auxiliary image satellite.
Content based on the various embodiments described above, according to phase compensating factor to assisted image carry out phase error compensation before
Further include: assisted image is pre-compensated for according to reference image.
It wherein, is the haplopia complex pattern for the ground object target that reference satellite acquires with reference to image.
It should be noted that phase error further includesThe common phase based on reference image can be used
Phase factor is eliminated in position compensation method
Inventive embodiments can accurately remove phase error by pre-compensating for according to reference image to assisted image.
Content based on the various embodiments described above, the specific steps packet that assisted image is pre-compensated for according to reference image
It includes: the factor is pre-compensated for according to reference image capturingIt is pre-compensated for according to precompensation factor pair assisted image.
Wherein, j indicates imaginary unit；R indicate ground object target distance to coordinate；λ indicates radar wavelength.
Specifically, according to reference image, the distance of available ground object target to coordinate and radar wavelength, so as to
Obtain the precompensation factor
The factor will be pre-compensated forIt is multiplied with assisted image, phase error can be eliminatedIt carries out pre-
Compensation deals and then according to phase compensating factor to precompensation processing after assisted image carry out phase error compensation.
Fig. 3 is that the phase error compensation chromatographed in three-dimensional imaging according to synthetic aperture radar provided in an embodiment of the present invention fills
The structural schematic diagram set.Content based on the various embodiments described above, as shown in figure 3, the device includes because of sub-acquisition module 301 and phase
Position compensating module 302, in which:
Because of sub-acquisition module 301, for obtaining phase according to the corresponding baseline parameter of assisted image for each assisted image
Position compensation factor；
Phase compensation block 302, for carrying out phase error compensation to assisted image according to phase compensating factor；
Wherein, the haplopia complex pattern for the ground object target that image satellite acquires supplemented by assisted image.
Specifically, because sub-acquisition module 301 can obtain phase compensation according to the corresponding baseline parameter of each assisted image
The factor.
Phase compensation block 302 carries out phase error compensation according to the phase compensating factor, to the assisted image, eliminates flat
Ground phaseAnd spatial offsetBring error
Phase error compensation device in synthetic aperture radar chromatography three-dimensional imaging provided in an embodiment of the present invention, for holding
The phase error compensation method in synthetic aperture radar chromatography three-dimensional imaging that row the various embodiments described above of the present invention provide, the synthesis
Aperture radar chromatography three-dimensional imaging in phase error compensation device include each module realize corresponding function specific method and
Process is detailed in the embodiment of the phase error compensation method in above-mentioned synthetic aperture radar chromatography three-dimensional imaging, no longer superfluous herein
It states.
The synthetic aperture radar chromatographs the synthesis that the phase error compensation device in three-dimensional imaging is used for foregoing embodiments
Aperture radar chromatographs the phase error compensation method in three-dimensional imaging.Therefore, the synthetic aperture radar in foregoing embodiments
The description and definition in the phase error compensation method in three-dimensional imaging are chromatographed, can be used for respectively executing mould in the embodiment of the present invention
The understanding of block.
The embodiment of the present invention obtains phase compensating factor according to the corresponding baseline parameter of assisted image, according to phase compensation because
Son carries out phase error compensation to assisted image, and energy is more acurrate, quickly and easily carries out phase error compensation, can improve compensation
Effect, so as to obtain the better synthetic aperture radar chromatography three-dimensional imaging result of effect.
Fig. 4 is the structural block diagram according to electronic equipment provided in an embodiment of the present invention.Content based on the above embodiment, such as
Shown in Fig. 4, which may include: processor (processor) 401, memory (memory) 402 and bus 403；Its
In, processor 401 and memory 402 pass through bus 403 and complete mutual communication；Processor 401 is stored in for calling
In reservoir 402 and the computer program instructions that can be run on processor 401, to execute provided by above-mentioned each method embodiment
Synthetic aperture radar chromatographs the phase error compensation method in three-dimensional imaging, for example, for each assisted image, according to auxiliary
The corresponding baseline parameter of image is helped to obtain phase compensating factor；Phase error benefit is carried out to assisted image according to phase compensating factor
It repays；Wherein, the haplopia complex pattern for the ground object target that image satellite acquires supplemented by assisted image.
Another embodiment of the present invention discloses a kind of computer program product, and computer program product is non-transient including being stored in
Computer program on computer readable storage medium, computer program include program instruction, when program instruction is held by computer
When row, the phase that computer is able to carry out in the chromatography three-dimensional imaging of synthetic aperture radar provided by above-mentioned each method embodiment is missed
Poor compensation method, for example, for each assisted image, according to the corresponding baseline parameter of assisted image obtain phase compensation because
Son；Phase error compensation is carried out to assisted image according to phase compensating factor；Wherein, image satellite acquisition supplemented by assisted image
The haplopia complex pattern of ground object target.
In addition, the logical order in above-mentioned memory 402 can be realized by way of SFU software functional unit and conduct
Independent product when selling or using, can store in a computer readable storage medium.Based on this understanding, originally
The technical solution of the inventive embodiments substantially part of the part that contributes to existing technology or the technical solution in other words
It can be embodied in the form of software products, which is stored in a storage medium, including several fingers
It enables and using so that a computer equipment (can be personal computer, server or the network equipment etc.) executes the present invention respectively
The all or part of the steps of a embodiment method.And storage medium above-mentioned includes: USB flash disk, mobile hard disk, read-only memory
(ROM, Read-Only Memory), random access memory (RAM, Random Access Memory), magnetic or disk
Etc. the various media that can store program code.
Another embodiment of the present invention provides a kind of non-transient computer readable storage medium, non-transient computer readable storages
Medium storing computer instruction, computer instruction make computer execute synthetic aperture radar provided by above-mentioned each method embodiment
Chromatograph the phase error compensation method in three-dimensional imaging, for example, corresponding according to assisted image for each assisted image
Baseline parameter obtains phase compensating factor；Phase error compensation is carried out to assisted image according to phase compensating factor；Wherein, it assists
The haplopia complex pattern for the ground object target that image satellite acquires supplemented by image.
The apparatus embodiments described above are merely exemplary, wherein unit can be as illustrated by the separation member
Or may not be and be physically separated, component shown as a unit may or may not be physical unit, i.e.,
It can be located in one place, or may be distributed over multiple network units.It can select according to the actual needs therein
Some or all of the modules achieves the purpose of the solution of this embodiment.Those of ordinary skill in the art are not paying creative labor
In the case where dynamic, it can understand and implement.
Through the above description of the embodiments, those skilled in the art can be understood that each embodiment can
It realizes by means of software and necessary general hardware platform, naturally it is also possible to pass through hardware.Such understanding, above-mentioned skill
Substantially the part that contributes to existing technology can be embodied in the form of software products art scheme in other words, the calculating
Machine software product may be stored in a computer readable storage medium, such as ROM/RAM, magnetic disk, CD, including some instructions are used
So that a computer equipment (can be personal computer, server or the network equipment etc.) executes above-mentioned each implementation
The method of certain parts of example or embodiment.
Finally, it should be noted that the above embodiments are merely illustrative of the technical solutions of the present invention, rather than its limitations；Although
Present invention has been described in detail with reference to the aforementioned embodiments, those skilled in the art should understand that: it still may be used
To modify the technical solutions described in the foregoing embodiments or equivalent replacement of some of the technical features；
And these are modified or replaceed, technical solution of various embodiments of the present invention that it does not separate the essence of the corresponding technical solution spirit and
Range.
Claims (10)
1. a kind of phase error compensation method in synthetic aperture radar chromatography three-dimensional imaging characterized by comprising
For each assisted image, phase compensating factor is obtained according to the corresponding baseline parameter of the assisted image；
Phase error compensation is carried out to the assisted image according to the phase compensating factor；
Wherein, the haplopia complex pattern for the ground object target that image satellite acquires supplemented by the assisted image.
2. the phase error compensation method in synthetic aperture radar chromatography three-dimensional imaging according to claim 1, feature
It is, it is described for each assisted image, the tool of phase compensating factor is obtained according to the corresponding baseline parameter of the assisted image
Body step includes:
For each assisted image, the inclination angle of the corresponding baseline of the assisted image and horizontal plane, the auxiliary shadow are obtained
As the length of corresponding baseline, reference satellite observe the distance of the visual angle on ground, radar wavelength and ground object target to coordinate；
According to the length of the corresponding baseline of assisted image baseline corresponding with the inclination angle of horizontal plane, the assisted image, ginseng
Examine the distance of the visual angle on moonscope ground, radar wavelength and ground object target to coordinate, obtain the comprehensive compensation factor.
3. the phase error compensation method in synthetic aperture radar chromatography three-dimensional imaging according to claim 1, feature
It is, it is described for each assisted image, the tool of phase compensating factor is obtained according to the corresponding baseline parameter of the assisted image
Body step includes:
For each assisted image, the inclination angle of the corresponding baseline of the assisted image and horizontal plane, the auxiliary shadow are obtained
As the length of corresponding baseline, the visual angle on reference satellite observation ground and radar wavelength, according to the corresponding base of the assisted image
The length of line baseline corresponding with the inclination angle of horizontal plane, the assisted image, the visual angle on reference satellite observation ground and radar wave
It is long, obtain the first compensation factor；Obtain ground object target distance to coordinate, according to the corresponding baseline of the assisted image and water
The inclination angle of plane, the length of the corresponding baseline of the assisted image, reference satellite observation ground visual angle, radar wavelength and atural object
The distance of target to coordinate, obtain the second compensation factor；
Alternatively, for each assisted image, the inclination angle, described auxiliary of the corresponding baseline of the assisted image and horizontal plane is obtained
Help the distance of the length of the corresponding baseline of image, the reference satellite observation visual angle on ground, radar wavelength and ground object target to seat
Mark；According to the length of the corresponding baseline of assisted image baseline corresponding with the inclination angle of horizontal plane, the assisted image, reference
The distance at the visual angle on moonscope ground, radar wavelength and ground object target to coordinate, obtain second compensation factor.
4. the phase error compensation method in synthetic aperture radar chromatography three-dimensional imaging according to claim 2, feature
It is, the comprehensive compensation factor is
Wherein, j indicates imaginary unit；biIndicate the corresponding baseline of assisted image；θ indicates the visual angle on reference satellite observation ground；αi
Indicate the inclination angle of assisted image corresponding baseline and horizontal plane；R indicate ground object target distance to coordinate；λ indicates radar wave
It is long；I indicates the number of auxiliary image satellite, and 1≤i≤N-1, N-1 indicate the sum of auxiliary image satellite.
5. the phase error compensation method in synthetic aperture radar chromatography three-dimensional imaging according to claim 3, feature
It is, first compensation factor isSecond compensation factor is
Wherein, j indicates imaginary unit；biIndicate the corresponding baseline of assisted image；θ indicates the visual angle on reference satellite observation ground；αi
Indicate the inclination angle of assisted image corresponding baseline and horizontal plane；R indicate ground object target distance to coordinate；λ indicates radar wave
It is long；I indicates the number of auxiliary image satellite, and 1≤i≤N-1, N-1 indicate the sum of auxiliary image satellite.
6. the phase error compensation method in synthetic aperture radar chromatography three-dimensional imaging according to any one of claims 1 to 5,
It is characterized in that, it is described according to the phase compensating factor to the assisted image carry out phase error compensation before further include:
The assisted image is pre-compensated for according to reference image；
Wherein, the haplopia complex pattern with reference to the ground object target that image is reference satellite acquisition.
7. the phase error compensation method in synthetic aperture radar chromatography three-dimensional imaging according to claim 6, feature
It is, the specific steps pre-compensated for according to reference image to the assisted image include:
The factor is pre-compensated for reference to image capturing according to describedAccording to assisted image described in the precompensation factor pair into
Row precompensation；
Wherein, j indicates imaginary unit；R indicate ground object target distance to coordinate；λ indicates radar wavelength.
8. the phase error compensation device in a kind of synthetic aperture radar chromatography three-dimensional imaging characterized by comprising
Because of sub-acquisition module, for obtaining phase compensation according to the corresponding baseline parameter of assisted image for each assisted image
The factor；
Phase compensation block, for carrying out phase error compensation to the assisted image according to the phase compensating factor；
Wherein, the haplopia complex pattern for the ground object target that image satellite acquires supplemented by the assisted image.
9. a kind of electronic equipment including memory, processor and stores the calculating that can be run on a memory and on a processor
Machine program, which is characterized in that the processor realizes synthesis as described in any one of claim 1 to 7 when executing described program
Aperture radar chromatographs the step of phase error compensation method in three-dimensional imaging.
10. a kind of non-transient computer readable storage medium, is stored thereon with computer program, which is characterized in that the computer
It is realized when program is executed by processor in synthetic aperture radar chromatography three-dimensional imaging as described in any one of claim 1 to 7
The step of phase error compensation method.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910665019.0A CN110297243B (en) | 2019-07-23 | 2019-07-23 | Phase error compensation method and device in synthetic aperture radar tomography three-dimensional imaging |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910665019.0A CN110297243B (en) | 2019-07-23 | 2019-07-23 | Phase error compensation method and device in synthetic aperture radar tomography three-dimensional imaging |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN110297243A true CN110297243A (en) | 2019-10-01 |
| CN110297243B CN110297243B (en) | 2021-03-30 |
Family
ID=68031636
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910665019.0A Active CN110297243B (en) | 2019-07-23 | 2019-07-23 | Phase error compensation method and device in synthetic aperture radar tomography three-dimensional imaging |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN110297243B (en) |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113702972A (en) * | 2021-08-30 | 2021-11-26 | 中国科学院空天信息创新研究院 | Airborne multi-channel radar amplitude-phase error estimation method based on terrain prior |
| CN116577784A (en) * | 2023-06-16 | 2023-08-11 | 北京建筑大学 | Microwave vision three-dimensional imaging method and device for synthetic aperture radar and electronic equipment |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102221697A (en) * | 2011-03-25 | 2011-10-19 | 电子科技大学 | Airborne multi-antenna SAR chromatography three dimensional imaging system and imaging method thereof |
| EP2802896B1 (en) * | 2012-01-10 | 2017-04-19 | Raytheon Company | Sar autofocus for ground penetration radar |
| CN109828272A (en) * | 2019-03-29 | 2019-05-31 | 云南电网有限责任公司电力科学研究院 | A kind of SAR difference chromatography method and device |
| CN110018474A (en) * | 2019-01-25 | 2019-07-16 | 北京理工大学 | Three-D imaging method based on geostationary orbit synthetic aperture radar chromatographic technique |
-
2019
- 2019-07-23 CN CN201910665019.0A patent/CN110297243B/en active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102221697A (en) * | 2011-03-25 | 2011-10-19 | 电子科技大学 | Airborne multi-antenna SAR chromatography three dimensional imaging system and imaging method thereof |
| EP2802896B1 (en) * | 2012-01-10 | 2017-04-19 | Raytheon Company | Sar autofocus for ground penetration radar |
| CN110018474A (en) * | 2019-01-25 | 2019-07-16 | 北京理工大学 | Three-D imaging method based on geostationary orbit synthetic aperture radar chromatographic technique |
| CN109828272A (en) * | 2019-03-29 | 2019-05-31 | 云南电网有限责任公司电力科学研究院 | A kind of SAR difference chromatography method and device |
Non-Patent Citations (2)
| Title |
|---|
| 孙希龙: "SAR层析与差分层析成像技术研究", 《中国博士学位论文全文数据库 信息科技辑》 * |
| 徐西桂 等: "多基线层析SAR技术的研究现状分析", 《测绘通报》 * |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113702972A (en) * | 2021-08-30 | 2021-11-26 | 中国科学院空天信息创新研究院 | Airborne multi-channel radar amplitude-phase error estimation method based on terrain prior |
| CN116577784A (en) * | 2023-06-16 | 2023-08-11 | 北京建筑大学 | Microwave vision three-dimensional imaging method and device for synthetic aperture radar and electronic equipment |
| CN116577784B (en) * | 2023-06-16 | 2024-04-26 | 北京建筑大学 | Synthetic aperture radar three-dimensional imaging method and device based on Fresnel diffraction zone |
Also Published As
| Publication number | Publication date |
|---|---|
| CN110297243B (en) | 2021-03-30 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN106780590B (en) | Method and system for acquiring depth map | |
| US8355564B2 (en) | Corresponding point searching method and three-dimensional position measuring method | |
| Li et al. | Image coregistration in SAR interferometry | |
| CN106249236B (en) | A kind of spaceborne InSAR long-short baselines image joint method for registering | |
| CN104730520B (en) | Circumference SAR back projection self-focusing method based on subaperture synthesis | |
| CN110349195A (en) | A kind of target object 3D measurement parameter acquisition methods, system and storage medium based on depth image | |
| CN111190182B (en) | Terahertz radar ultrahigh-resolution imaging method | |
| CN103544710A (en) | Image registration method | |
| CN107390216A (en) | High speed super-resolution stationary point scan imaging method based on wave-number domain coherence factor | |
| CN110297243A (en) | Synthetic aperture radar chromatographs the phase error compensation method and device in three-dimensional imaging | |
| Eichhardt et al. | Affine correspondences between central cameras for rapid relative pose estimation | |
| CN103824287A (en) | Robust plane fitting-based phase correlation sub-pixel matching method | |
| CN111239736A (en) | Single-baseline-based surface elevation correction method, device, equipment and storage medium | |
| CN112415515A (en) | Method for separating targets with different heights by airborne circular track SAR | |
| CN113311433A (en) | InSAR interferometric phase two-step unwrapping method combining quality map and minimum cost flow | |
| Gao et al. | A novel two-step noise reduction approach for interferometric phase images | |
| CN110297242A (en) | Compressed sensing based synthetic aperture radar chromatography three-D imaging method and device | |
| CN105353374B (en) | A kind of single-frequency radar imaging method for the target that spins | |
| CN103871030B (en) | The filtering method and equipment of a kind of interference image | |
| CN103942775A (en) | Phase related-sub-pixel matching method based on maximum-kernel-density estimation | |
| CN107610161B (en) | A kind of InSAR fast image registration method based on Quadtree Partition | |
| CN114002666B (en) | Method and equipment for extracting satellite-borne ATI-SAR ocean current flow velocity under any antenna configuration | |
| CN113640798A (en) | Radar target multi-angle reconstruction method and device and storage medium | |
| CN115409735A (en) | Filtering method and device for InSAR interference image | |
| CN114241022A (en) | Unmanned aerial vehicle image automatic registration method and system |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant | ||
| EE01 | Entry into force of recordation of patent licensing contract | ||
| EE01 | Entry into force of recordation of patent licensing contract |
Application publication date: 20191001 Assignee: Beijing Xingtong Weishi Technology Co.,Ltd. Assignor: Beijing University of Civil Engineering and Architecture Contract record no.: X2023980035223 Denomination of invention: Phase Error Compensation Method and Device in Synthetic Aperture Radar Tomography 3D Imaging Granted publication date: 20210330 License type: Common License Record date: 20230504 |