CN110297242A - Compressed sensing based synthetic aperture radar chromatography three-D imaging method and device - Google Patents
Compressed sensing based synthetic aperture radar chromatography three-D imaging method and device Download PDFInfo
- Publication number
- CN110297242A CN110297242A CN201910664899.XA CN201910664899A CN110297242A CN 110297242 A CN110297242 A CN 110297242A CN 201910664899 A CN201910664899 A CN 201910664899A CN 110297242 A CN110297242 A CN 110297242A
- Authority
- CN
- China
- Prior art keywords
- complex pattern
- image
- phase
- level land
- interference
- 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/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 a kind of compressed sensing based synthetic aperture radar chromatography three-D imaging method and device.Wherein, method includes: for each interference complex pattern, and the level land phase in the interferometric phase of complex pattern is interfered in removal, obtains the interferometric phase gone after level land；According to compression sensing method, respectively interfere the corresponding assisted image of complex pattern and remove the interferometric phase after level land, obtains the three-dimensional imaging result of ground object target；Wherein, the haplopia complex pattern for the ground object target that image satellite acquires supplemented by assisted image；It interferes complex pattern, is obtained after assisted image is registrated with reference image；It is the haplopia complex pattern for the ground object target that reference satellite acquires with reference to image.Compressed sensing based synthetic aperture radar chromatography three-D imaging method and device provided in an embodiment of the present invention can more accurately carry out phase error compensation by the phase error compensation method removed based on elliptical earth phase, can obtain the better imaging results of effect.
Description
Technical field
The present invention relates to electronic signal process technical fields, more particularly, to a kind of compressed sensing based synthesis hole
Diameter radar chromatographs three-D imaging method and device.
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, so as to cause synthetic aperture thunder
It is poor up to the low efficiency of chromatography three-dimensional imaging, imaging effect.
Summary of the invention
The embodiment of the present invention provides a kind of compressed sensing based synthetic aperture radar chromatography three-D imaging method and device,
To solve or at least be partially solved the defect of prior art imaging effect difference.
In a first aspect, the embodiment of the present invention provides a kind of compressed sensing based synthetic aperture radar chromatography three-dimensional imaging side
Method, comprising:
For each interference complex pattern, the level land phase in the interferometric phase of the interference complex pattern is removed, acquisition goes to put down
Interferometric phase after ground；
According to compression sensing method, the corresponding assisted image of each interference complex pattern and the interference gone after level land
Phase obtains the three-dimensional imaging result of ground object target；
Wherein, the haplopia complex pattern for the ground object target that image satellite acquires supplemented by the assisted image；The interference
Complex pattern, the assisted image obtain after being registrated with reference image；It is described with reference to image be reference satellite acquire
The haplopia complex pattern of the ground object target.
Preferably, the specific steps of the level land phase in the interferometric phase of the removal interference complex pattern include:
Interpolation is carried out to the interference complex pattern, Fourier transformation is carried out to the interference complex pattern after interpolation, is obtained
The maximum spectrum of the interferometric phase of the interference complex pattern；
The maximum spectrum is moved at zero-frequency spectrum, the interferometric phase after level land is removed described in acquisition.
Preferably, the specific steps for carrying out interpolation to the interference complex pattern include:
The interpolation that multiple is not less than 50 times is carried out to the interference complex pattern.
Preferably, according to compression sensing method, the corresponding assisted image of each interference complex pattern and it is described go level land it
Interferometric phase afterwards, the specific steps for obtaining the three-dimensional imaging result of ground object target include:
Each interference complex pattern corresponding assisted image progress first phase error is gone with reference to image according to described
It removes；
According to compression sensing method, carry out the corresponding auxiliary of each interference complex pattern after the removal of first phase error
Image and the interferometric phase gone after level land, obtain the three-dimensional imaging result of ground object target.
Preferably, it is described according to described with reference to image to each interference complex pattern corresponding assisted image progress first
Phase error removal specific steps include:
Image is referred to according to described, obtains compensation factor；
First phase error is carried out to the corresponding assisted image of each interference complex pattern according to the compensation factor to go
It removes.
Preferably, the compensation factor is
Wherein, j indicates imaginary unit；R indicate the distance of the ground object target to coordinate；λ indicates radar wavelength.
Preferably for each interference complex pattern, the level land phase in the interferometric phase of the interference complex pattern is removed, is obtained
Before removing the interferometric phase after level land further include:
Each assisted image is registrated with described with reference to image respectively, obtains each interference complex pattern.
Second aspect, the embodiment of the present invention provide a kind of compressed sensing based synthetic aperture radar chromatography three-dimensional imaging dress
It sets, comprising:
Level land removes module, for removing in the interferometric phase for interfering complex pattern for each interference complex pattern
Level land phase obtains the interferometric phase gone after level land；
Image-forming module is compressed, for according to compression sensing method, the corresponding assisted image of each interference complex pattern and institute
The interferometric phase after level land is stated, the three-dimensional imaging result of ground object target is obtained；
Wherein, the haplopia complex pattern for the ground object target that image satellite acquires supplemented by the assisted image；The interference
Complex pattern, the assisted image obtain after being registrated with reference image；It is described with reference to image be reference satellite acquire
The haplopia complex pattern of the ground object target.
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
Compressed sensing based synthetic aperture radar provided by any possible implementation chromatographs three-D imaging method in mode
The step of.
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
Provided by the implementation of energy the step of compressed sensing based synthetic aperture radar chromatography three-D imaging method.
Compressed sensing based synthetic aperture radar chromatography three-D imaging method and device provided in an embodiment of the present invention, lead to
The phase error compensation method based on elliptical earth phase removal is crossed, realizes compressed sensing based synthetic aperture radar three-dimensional imaging,
Energy is more acurrate, simply, easily carries out phase error compensation, can obtain the better synthetic aperture radar chromatography three-dimensional imaging of effect
As a 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 to chromatograph three-D imaging method according to compressed sensing based synthetic aperture radar provided in an embodiment of the present invention
Flow diagram；
Fig. 2 is the geometrical model schematic diagram of TomoSAR system provided in an embodiment of the present invention；
Fig. 3 is to chromatograph three-dimensional image forming apparatus according to compressed sensing based synthetic aperture radar provided in an embodiment of the present invention
Structural schematic diagram；
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 provides a kind of compressed sensing based synthetic aperture
Radar chromatography three-D imaging method and device, inventive concept are to be removed using the method for flat earth because elliptical earth phase draws
The phase error risen carries out third dimension imaging using the interferometric phase gone after level land, more preferable so as to obtain imaging effect
Imaging results.
Fig. 1 is to chromatograph three-D imaging method according to compressed sensing based synthetic aperture radar provided in an embodiment of the present invention
Flow diagram.As shown in Figure 1, this method comprises: step S101, for each interference complex pattern, removal interference complex pattern
Interferometric phase in level land phase, obtain and remove interferometric phase after level land.
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 compressed sensing based synthetic aperture radar that bright embodiment provides chromatographs three-D imaging method, 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 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).
Practical geometrical relationship shown according to fig. 2, s is far smaller than r, therefore has
To above formula in biTaylor expansion is carried out at=0, is had
Similarly, since s is far smaller than r, then have
Therefore, the three-dimensional mathematical model of SAR image can be further simplified as
It enablesThen have
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.
Compare si(x, r, s) and giThe formula of (x, r, s), it can be deduced that include in phase error
siTwo of them phase factor in the formula of (x, r, s)With exp (j2 π fiS) it constitutes
Interferometric phase between the assisted image and reference image of i-th auxiliary image satellite acquisition.Wherein,
Corresponding to the elliptical earth phase phase in interferometric phase, exp (j2 π fiS) correspond to the interference as caused by oblique distance vertical variation
Phase.
Interferometric phase can be obtained by interference complex pattern.
Image satellite s auxiliary for i-thiThe assisted image of (i=1,2 ..., N-1) acquisition, by the assisted image and reference
After image is registrated, corresponding interference complex pattern is obtained.
According to the interference complex pattern, the interferometric phase of the interference complex pattern can be extracted, interferometric phase is in striated, and
The level land phase in the interferometric phase is obtained according to the interferometric phase of the interference complex pattern, and removes the level land phase, realization pair
The removal of elliptical earth phase.
Removal to elliptical earth phase, referred to as flat earth, referred to as " going to level land ".
Level land phase can be obtained using any method for going to level land, can accordingly remove the level land phase in interferometric phase
Phase error compensation is realized in position, so as to obtain the interferometric phase after level land, i.e. acquisition phase factor exp (j2 π
fis)。
For example, level land phase can be obtained using the level land phase frequency spectrum estimation technique, and carry out corresponding removal level land phase
Processing.
Step S102, it according to compression sensing method, the corresponding assisted image of each interference complex pattern and goes dry after level land
Phase is related to, the three-dimensional imaging result of ground object target is obtained.
Wherein, the haplopia complex pattern for the ground object target that image satellite acquires supplemented by assisted image；Interfere complex pattern, is auxiliary
What image obtained after being registrated with reference image；It 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
After having executed phase error compensation step, the three-dimensional mathematical model of SAR image becomes
As can be seen from the above equation, i-th auxiliary image satellite si(i=1,2 ..., N-1) obtained according to compression sensing method
SAR image corresponds to spectrum component fiDistance it is vertical on ground object target discrete Fourier transform.
It therefore, can be directly according to going the interferometric phase after level land to carry out third dimension imaging, it is hereby achieved that atural object
The three-dimensional imaging result of target.
According to each auxiliary image satellite (i.e. s1To sN-1) SAR image that obtains, so as to obtain ground object target it is three-dimensional at
As a result, obtaining the three-dimensional point cloud of ground object target.
The embodiment of the present invention is realized compressed sensing based by the phase error compensation method removed based on elliptical earth phase
Synthetic aperture radar three-dimensional imaging, energy is more acurrate, simply, easily carries out phase error compensation, can obtain effect and preferably close
Three-dimensional imaging result is chromatographed at aperture radar.
The specific steps of the level land phase in the interferometric phase of complex pattern are interfered in content based on the various embodiments described above, removal
Include: that interpolation is carried out to interference complex pattern, Fourier transformation is carried out to the interference complex pattern after interpolation, obtains interference complex pattern
The maximum spectrum of interferometric phase；Maximum spectrum is moved at zero-frequency spectrum, the interferometric phase gone after level land is obtained.
Specifically, the level land in the interferometric phase of each interference complex pattern can be obtained using the level land phase frequency spectrum estimation technique
Phase.
In order to obtain more accurate level land phase, before it complex pattern will be interfered to transform to frequency domain by Fourier transformation,
Interpolation processing first is carried out to interference data (i.e. interference complex pattern).
By Fourier transformation, the interference complex pattern after interpolation is transformed into frequency domain, so as to obtain interference complex pattern
Interferometric phase maximum spectrum.
The maximum spectrum is exactly frequency spectrum corresponding to the phase of level land, therefore, which can be moved to zero-frequency spectrum
Place retains the interferometric phase gone after level land so as to complete level land phase.
The embodiment of the present invention obtains level land phase by the level land phase frequency spectrum estimation technique and removes level land phase, can obtain more
Accurate level land phase, so as to it is more acurrate, simply, easily carry out phase error compensation, effect can be obtained and preferably synthesized
Aperture radar chromatographs three-dimensional imaging result.
Content based on the various embodiments described above includes: to the multiple figure of interference to the specific steps that interference complex pattern carries out interpolation
The interpolation for being not less than 50 times as carrying out multiple.
Specifically, when carrying out interpolation to interference complex pattern, the interpolation processing of 50 times of interpolation or more high power is specifically carried out.
By the interpolation processing of high magnification numbe, more accurate level land phase can be obtained using the level land phase frequency spectrum estimation technique.
The interpolation that the embodiment of the present invention is not less than 50 times by carrying out multiple can be obtained based on the level land phase frequency spectrum estimation technique
More accurate level land phase can obtain the better synthetic aperture radar of effect so as to more accurately carry out phase error compensation
Chromatograph three-dimensional imaging result.
Content based on the various embodiments described above, according to compression sensing method, the corresponding assisted image of each interference complex pattern and
The interferometric phase after level land is removed, the specific steps for obtaining the three-dimensional imaging result of ground object target include: according to reference image pair
The corresponding assisted image of each interference complex pattern carries out the removal of first phase error.
Specifically, first phase error can be eliminated after removing elliptical earth phase
First phase error can be eliminated using the common phase compensating method based on reference image
According to the corresponding assisted image of each interference complex pattern after compression sensing method, progress first phase error removal
With remove the interferometric phase after level land, obtain the three-dimensional imaging result of ground object target.
Specifically, phase error is removedWithIt later, can be according to compressed sensing
The corresponding TomoSAR imaging algorithm such as method, utilizes the three-dimensional mathematical model of SAR imageObtain the SAR figure of each auxiliary image satellite
Picture, so as to obtain the three-dimensional imaging of ground object target according to the SAR image of each auxiliary image satellite as a result, obtaining atural object mesh
Target three-dimensional point cloud.
The embodiment of the present invention can remove phase error, more accurately by carrying out the removal of first phase error so as to obtain
It obtains the better synthetic aperture radar of effect and chromatographs three-dimensional imaging result.
Content based on the various embodiments described above carries out the corresponding assisted image of each interference complex pattern according to reference image
The specific steps of first phase error removal include: to obtain compensation factor according to reference image；According to compensation factor to each dry
It relates to the corresponding assisted image of complex pattern and carries out the removal of first phase error.
Specifically, available compensation factor corresponding with first phase error, it is corresponding for each interference complex pattern
Assisted image removes first phase error by the assisted image multiplied by the compensation factor.
The embodiment of the present invention is schemed by obtaining compensation factor according to reference image according to compensation factor is multiple to each interference
As the progress first phase error removal of corresponding assisted image, phase error can be accurately removed, it is more preferable so as to obtain effect
Synthetic aperture radar chromatograph three-dimensional imaging result.
Content based on the various embodiments described above, compensation factor are
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 compensation factor
The distance of ground object target to coordinate, between the scattering point range-azimuth unit fallen in and reference satellite away from
From.
Content based on the various embodiments described above, for each interference complex pattern, removal is interfered in the interferometric phase of complex pattern
Level land phase, obtain before removing the interferometric phase after level land further include: by each assisted image respectively with reference image into
Row registration, obtains each interference complex pattern.
Specifically, for si(i=1,2 ..., N-1) acquisition each assisted image, by each assisted image respectively with reference
Image is registrated, and the corresponding interference complex pattern of the assisted image is obtained.
The embodiment of the present invention obtains the multiple figure of each interference by the way that each assisted image to be registrated with reference image respectively
Picture, so as to extract interferometric phase and level land phase according to interference complex pattern, and then by going to level land to realize more accurate phase
Position error compensation can obtain the better synthetic aperture radar chromatography three-dimensional imaging result of effect.
It should be noted that after carrying out phase error compensation distance can also be obtained according to phase error compensation result
Bandwidth on vertical, apart from the vertical sparse number of maximum that can be reconstructed, be able to carry out the most short baseline total length effectively reconstructed.
Phase error compensation is carried out, phase error is removedWithLater, SAR schemes
The three-dimensional mathematical model of picture is
SAR image is the bandwidth B on vertical in the third dimensionWsFor
It follows that in the resolution ratio ρ on verticalsFor
Therefore, the vertical range of the distance that can be restored is
Wherein, r indicate ground object target distance to coordinate；λ indicates radar wavelength.
According to the resolution ratio apart from vertical range and on vertical, can obtain dilute apart from the vertical maximum that can be reconstructed
Dredge number KmaxFor
Apart from the vertical sparse number of maximum that can be reconstructed, refer to the distance effectively reconstructed it is vertical on maximum number of points.?
That is meeting RIP (limited equidistant) property when obtaining the point target position on vertical using compressed sensing technology
In the case where, being at best able to the sparse number effectively reconstructed is Kmax。
It chromatographs sparse several Ks of the SAR on oblique distance is vertical and is usually quite small, often take 3 or 4, K < < Kmax, therefore,
The most short baseline total length effectively reconstructed can be determined to according to above-mentioned condition.
The value (by taking K=4 as an example) of K is substituted intoIt is availableThat is, when K=4, it can
The most short baseline total length effectively reconstructed isIf baseline total length be less thanCompression sense can not then be carried out
Know.
The matrix right and wrong that RIP property (limited equidistant property) requires the every 2K column vector extracted from perception matrix to constitute
Unusual, this is to guarantee to perceive matrix and two different sparse signals will not be mapped in the same set (guarantee original
One by one mapping relations of the space to evacuated space).When this condition stub carries out compressed sensing imaging in SAR image, even if
Oblique distance is vertical, and there are sparse scattering points, but these sparse scattering points cannot be too near to, it is ensured that extract from perception matrix
The matrix that column vector is constituted is nonsingular it is necessary to guarantee that scattering point cannot be fallen in the vertical resolution cell of the same oblique distance, otherwise, nothing
Method distinguishes two points.Intrinsic resolution capability possessed by this and system itself are vertical in oblique distance is explained consistent.
According to it is above-mentioned apart from the vertical sparse number of maximum that can be reconstructed, be able to carry out the most short baseline overall length effectively reconstructed
Degree, can more effectively dispose TomoSAR system, avoid the case where can not carrying out compressed sensing, improve the effect of SAR tomography
Rate and success rate.
Fig. 3 is to chromatograph three-dimensional image forming apparatus according to compressed sensing based synthetic aperture radar provided in an embodiment of the present invention
Structural schematic diagram.Content based on the various embodiments described above, as shown in figure 3, the device includes level land removal module 301 and compression
Image-forming module 302, in which:
Level land removes module 301, for for each interference complex pattern, removal to be interfered flat in the interferometric phase of complex pattern
Ground phase obtains the interferometric phase gone after level land；
Image-forming module 302 is compressed, for according to compression sensing method, the corresponding assisted image of each interference complex pattern and going to put down
Interferometric phase after ground obtains the three-dimensional imaging result of ground object target；
Wherein, the haplopia complex pattern for the ground object target that image satellite acquires supplemented by assisted image；Interfere complex pattern, is auxiliary
What image obtained after being registrated with reference image；It is the haplopia complex pattern for the ground object target that reference satellite acquires with reference to image.
Specifically, removal module 301 in level land can extract the interferometric phase of each interference complex pattern, multiple according to the interference
The interferometric phase of image obtains the level land phase in the interferometric phase, and removes the level land phase, and elliptical earth phase is removed in realization
It removes, retains the interferometric phase gone after level land.
Compression image-forming module 302 directly basis can go the interferometric phase after level land to carry out third dimension imaging, so as to
To obtain the SAR image that each auxiliary image satellite obtains, according to the SAR image that each auxiliary image satellite obtains, so as to obtain
The three-dimensional imaging result of ground object target.
Compressed sensing based synthetic aperture radar provided in an embodiment of the present invention chromatographs three-dimensional image forming apparatus, for executing
The compressed sensing based synthetic aperture radar that the various embodiments described above of the present invention provide chromatographs three-D imaging method, should be based on compression
Each module that the synthetic aperture radar chromatography three-dimensional image forming apparatus of perception includes realizes that the specific method of corresponding function and process are detailed
See the embodiment of above-mentioned compressed sensing based synthetic aperture radar chromatography three-D imaging method, details are not described herein again.
The compressed sensing based synthetic aperture radar chromatography three-dimensional image forming apparatus is for foregoing embodiments based on pressure
The synthetic aperture radar of contracting perception chromatographs three-D imaging method.Therefore, the compressed sensing based conjunction in foregoing embodiments
At the description and definition in aperture radar chromatography three-D imaging method, it can be used for the reason of each execution module in the embodiment of the present invention
Solution.
The embodiment of the present invention is realized compressed sensing based by the phase error compensation method removed based on elliptical earth phase
Synthetic aperture radar three-dimensional imaging, energy is more acurrate, simply, easily carries out phase error compensation, can obtain effect and preferably close
Three-dimensional imaging result is chromatographed at aperture radar.
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
Compressed sensing based synthetic aperture radar chromatographs three-D imaging method, for example, for each interference complex pattern, removal is dry
The level land phase in the interferometric phase of complex pattern is related to, the interferometric phase gone after level land is obtained；According to compression sensing method, each dry
It relates to the corresponding assisted image of complex pattern and removes the interferometric phase after level land, obtain the three-dimensional imaging result of ground object target；Wherein,
The haplopia complex pattern for the ground object target that image satellite acquires supplemented by assisted image；Interfere complex pattern, be assisted image with reference to shadow
As being obtained after being registrated；It is the haplopia complex pattern for the ground object target that reference satellite acquires with reference to 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, it is three-dimensional that computer is able to carry out compressed sensing based synthetic aperture radar chromatography provided by above-mentioned each method embodiment
Imaging method, for example, for each interference complex pattern, removal is interfered the level land phase in the interferometric phase of complex pattern, obtained
Remove the interferometric phase after level land；According to compression sensing method, the corresponding assisted image of each interference complex pattern and go level land it
Interferometric phase afterwards obtains the three-dimensional imaging result of ground object target；Wherein, the atural object mesh that image satellite acquires supplemented by assisted image
Target haplopia complex pattern；It interferes complex pattern, is obtained after assisted image is registrated with reference image；It is reference with reference to image
The haplopia complex pattern of the ground object target of satellite acquisition.
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 execute computer provided by above-mentioned each method embodiment based on compressed sensing
Synthetic aperture radar chromatograph three-D imaging method, for example, for each interference complex pattern, removal interference complex pattern it is dry
The level land phase in phase is related to, the interferometric phase gone after level land is obtained；It is corresponding according to compression sensing method, each interference complex pattern
Assisted image and remove the interferometric phase after level land, obtain the three-dimensional imaging result of ground object target；Wherein, supplemented by assisted image
The haplopia complex pattern of the ground object target of image satellite acquisition；Interfere complex pattern, is after assisted image is registrated with reference image
It obtains；It is the haplopia complex pattern for the ground object target that reference satellite acquires with reference to 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 compressed sensing based synthetic aperture radar chromatographs three-D imaging method characterized by comprising
For each interference complex pattern, remove the level land phase in the interferometric phase of the interference complex pattern, acquisition go level land it
Interferometric phase afterwards；
According to compression sensing method, the corresponding assisted image of each interference complex pattern and the interference phase gone after level land
Position, obtains the three-dimensional imaging result of ground object target；
Wherein, the haplopia complex pattern for the ground object target that image satellite acquires supplemented by the assisted image；The multiple figure of interference
Picture, the assisted image obtain after being registrated with reference image；It is described with reference to image be reference satellite acquisition described in
The haplopia complex pattern of ground object target.
2. compressed sensing based synthetic aperture radar according to claim 1 chromatographs three-D imaging method, feature exists
In the specific steps of the level land phase in the interferometric phase of the removal interference complex pattern include:
Interpolation is carried out to the interference complex pattern, Fourier transformation is carried out to the interference complex pattern after interpolation, described in acquisition
Interfere the maximum spectrum of the interferometric phase of complex pattern；
The maximum spectrum is moved at zero-frequency spectrum, the interferometric phase after level land is removed described in acquisition.
3. compressed sensing based synthetic aperture radar according to claim 2 chromatographs three-D imaging method, feature exists
In the specific steps for carrying out interpolation to the interference complex pattern include:
The interpolation that multiple is not less than 50 times is carried out to the interference complex pattern.
4. compressed sensing based synthetic aperture radar according to claim 1 chromatographs three-D imaging method, feature exists
In, according to compression sensing method, the corresponding assisted image of each interference complex pattern and the interferometric phase gone after level land,
The specific steps of three-dimensional imaging result for obtaining ground object target include:
According to it is described with reference to image to the corresponding assisted image progress first phase error removal of each interference complex pattern；
According to compression sensing method, carry out the corresponding assisted image of each interference complex pattern after the removal of first phase error
With the interferometric phase gone after level land, the three-dimensional imaging result of ground object target is obtained.
5. compressed sensing based synthetic aperture radar according to claim 4 chromatographs three-D imaging method, feature exists
In, it is described according to described with reference to image to the corresponding assisted image progress first phase error removal of each interference complex pattern
Specific steps include:
Image is referred to according to described, obtains compensation factor；
The removal of first phase error is carried out to the corresponding assisted image of each interference complex pattern according to the compensation factor.
6. compressed sensing based synthetic aperture radar according to claim 5 chromatographs three-D imaging method, feature exists
In the compensation factor is
Wherein, j indicates imaginary unit；R indicate the distance of the ground object target to coordinate；λ indicates radar wavelength.
7. compressed sensing based synthetic aperture radar according to any one of claims 1 to 6 chromatographs three-D imaging method,
It is characterized in that, for each interference complex pattern, removes the level land phase in the interferometric phase of the interference complex pattern, acquisition goes to put down
Before interferometric phase after ground further include:
Each assisted image is registrated with described with reference to image respectively, obtains each interference complex pattern.
8. a kind of compressed sensing based synthetic aperture radar chromatographs three-dimensional image forming apparatus characterized by comprising
Level land removes module, for removing the level land in the interferometric phase of the interference complex pattern for each interference complex pattern
Phase obtains the interferometric phase gone after level land；
Image-forming module is compressed, for according to compression sensing method, the corresponding assisted image of each interference complex pattern and described going
Interferometric phase after level land obtains the three-dimensional imaging result of ground object target；
Wherein, the haplopia complex pattern for the ground object target that image satellite acquires supplemented by the assisted image；The multiple figure of interference
Picture, the assisted image obtain after being registrated with reference image；It is described with reference to image be reference satellite acquisition described in
The haplopia complex pattern of ground object target.
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 as described in any one of claim 1 to 7 be based on when executing described program
The synthetic aperture radar of compressed sensing chromatographs the step of three-D imaging method.
10. a kind of non-transient computer readable storage medium, is stored thereon with computer program, which is characterized in that the computer
Compressed sensing based synthetic aperture radar layer as described in any one of claim 1 to 7 is realized when program is executed by processor
The step of analysing three-D imaging method.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910664899.XA CN110297242A (en) | 2019-07-23 | 2019-07-23 | Compressed sensing based synthetic aperture radar chromatography three-D imaging method and device |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910664899.XA CN110297242A (en) | 2019-07-23 | 2019-07-23 | Compressed sensing based synthetic aperture radar chromatography three-D imaging method and device |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN110297242A true CN110297242A (en) | 2019-10-01 |
Family
ID=68031609
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910664899.XA Pending CN110297242A (en) | 2019-07-23 | 2019-07-23 | Compressed sensing based synthetic aperture radar chromatography three-D imaging method and device |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN110297242A (en) |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111640069A (en) * | 2020-04-17 | 2020-09-08 | 上海交通大学 | Compressive imaging method, system and device based on light sensing network and phase compensation |
| CN114545407A (en) * | 2021-09-24 | 2022-05-27 | 中国科学院精密测量科学与技术创新研究院 | Satellite-borne differential tomography SAR imaging method based on distributed compressed sensing |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101581779A (en) * | 2008-05-14 | 2009-11-18 | 中国科学院电子学研究所 | Method for generating three-dimensional imaging original echoed signals of chromatography synthetic aperture radars |
| CN102183766A (en) * | 2011-01-25 | 2011-09-14 | 中国船舶重工集团公司第七一五研究所 | Stratum chromatography method by synthetic aperture |
-
2019
- 2019-07-23 CN CN201910664899.XA patent/CN110297242A/en active Pending
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101581779A (en) * | 2008-05-14 | 2009-11-18 | 中国科学院电子学研究所 | Method for generating three-dimensional imaging original echoed signals of chromatography synthetic aperture radars |
| CN102183766A (en) * | 2011-01-25 | 2011-09-14 | 中国船舶重工集团公司第七一五研究所 | Stratum chromatography method by synthetic aperture |
Non-Patent Citations (3)
| Title |
|---|
| 孙希龙 等: "基于模拟干涉相位去斜的SAR层析处理方法", 《国防科技大学学报》 * |
| 彭星 等: "基于干涉信号的SAR层析技术重建建筑物三维模型", 《中国地球科学联合学术年会》 * |
| 郭立文: "多基线InSAR数据处理技术及应用研究", 《中国优秀硕士学位论文全文数据库 信息科技辑》 * |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111640069A (en) * | 2020-04-17 | 2020-09-08 | 上海交通大学 | Compressive imaging method, system and device based on light sensing network and phase compensation |
| CN111640069B (en) * | 2020-04-17 | 2022-11-11 | 上海交通大学 | Compressive imaging method, system and device based on light sensing network and phase compensation |
| CN114545407A (en) * | 2021-09-24 | 2022-05-27 | 中国科学院精密测量科学与技术创新研究院 | Satellite-borne differential tomography SAR imaging method based on distributed compressed sensing |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Bergamasco et al. | WASS: An open-source pipeline for 3D stereo reconstruction of ocean waves | |
| WO2018127007A1 (en) | Depth image acquisition method and system | |
| CN104730520B (en) | Circumference SAR back projection self-focusing method based on subaperture synthesis | |
| CN108038902A (en) | A kind of high-precision three-dimensional method for reconstructing and system towards depth camera | |
| CN110349195A (en) | A kind of target object 3D measurement parameter acquisition methods, system and storage medium based on depth image | |
| CN103713287B (en) | A kind of height reconstruction method based on relatively prime many baselines and device | |
| CN112415515B (en) | Method for separating targets with different heights by airborne circular track SAR | |
| CN111239736B (en) | Single-baseline-based surface elevation correction method, device, equipment and storage medium | |
| CN107390216A (en) | High speed super-resolution stationary point scan imaging method based on wave-number domain coherence factor | |
| CN110297242A (en) | Compressed sensing based synthetic aperture radar chromatography three-D imaging method and device | |
| CN113311433B (en) | InSAR interferometric phase two-step unwrapping method combining quality map and minimum cost flow | |
| CN109239710B (en) | Method and device for acquiring radar elevation information and computer-readable storage medium | |
| CN114663597B (en) | Real-time structured light reconstruction method and device based on normalized expanded polar line geometry | |
| CN112034463A (en) | High-precision settlement monitoring method and device for SAR data | |
| Gao et al. | A novel two-step noise reduction approach for interferometric phase images | |
| CN110297243A (en) | Synthetic aperture radar chromatographs the phase error compensation method and device in three-dimensional imaging | |
| CN103630901B (en) | The method of array 3-D SAR imaging is looked under airborne | |
| CN103871030B (en) | The filtering method and equipment of a kind of interference image | |
| CN103076608A (en) | Contour-enhanced beaming-type synthetic aperture radar imaging method | |
| CN103942775B (en) | Phase place based on maximum kernel density Estimation is correlated with sub-pix matching process | |
| CN107610161B (en) | A kind of InSAR fast image registration method based on Quadtree Partition | |
| CN106093902B (en) | The scattering signatures transform method of similar purpose | |
| CN114002666A (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 | |
| CN113671497A (en) | Single-channel SAR target three-dimensional coordinate extraction method based on cylindrical symmetric model |
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: 20191001 |