CN108445458A - A kind of synthetic aperture radar orbit error removing method and device - Google Patents
A kind of synthetic aperture radar orbit error removing method and device Download PDFInfo
- Publication number
- CN108445458A CN108445458A CN201810207673.2A CN201810207673A CN108445458A CN 108445458 A CN108445458 A CN 108445458A CN 201810207673 A CN201810207673 A CN 201810207673A CN 108445458 A CN108445458 A CN 108445458A
- Authority
- CN
- China
- Prior art keywords
- pixel
- synthetic aperture
- aperture radar
- data
- phase
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
- G01S13/88—Radar or analogous systems specially adapted for specific applications
- G01S13/89—Radar or analogous systems specially adapted for specific applications for mapping or imaging
- G01S13/90—Radar or analogous systems specially adapted for specific applications for mapping or imaging using synthetic aperture techniques, e.g. synthetic aperture radar [SAR] techniques
- G01S13/9004—SAR image acquisition techniques
- G01S13/9019—Auto-focussing of the SAR signals
-
- 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
- G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
- G01S13/88—Radar or analogous systems specially adapted for specific applications
- G01S13/89—Radar or analogous systems specially adapted for specific applications for mapping or imaging
- G01S13/90—Radar or analogous systems specially adapted for specific applications for mapping or imaging using synthetic aperture techniques, e.g. synthetic aperture radar [SAR] techniques
- G01S13/904—SAR modes
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
A kind of synthetic aperture radar orbit error removing method of offer of the embodiment of the present invention and device, are related to the microwave remote sensing technique field of remote sensing image.Wherein, the present invention is to carry out spectrum estimation based on data of synthetic aperture radar itself and main frequency is zeroed, theoretical tight, and does not need other data and be supported, and can effectively improve data processing precision.
Description
Technical field
The present invention relates to the microwave remote sensing technique fields of remote sensing image, in particular to a kind of synthetic aperture radar rail
Road error cancelling method and device.
Background technology
Synthetic aperture radar interferometry is that have the haplopia of small visual angle difference multiple using two scapes of covering the same area
Data form interference fringe picture and obtain the technology of target elevation or three-dimensional coordinate by being further processed, differential interferometry synthesis
Aperture radar technique is that the emulation phase for generating interference fringe picture and DEM carries out difference processing, can be obtained high-precisionly
Table deformation data.
But in the track of satellite, there are differential phase can be caused to there is remaining striped when error so that result of calculation is deposited
In systematic error, while the accuracy of observation of deformation data is greatly reduced.Currently, removing the main side of orbit error in the prior art
Method the precise orbit data to be released such as has, carries out calibrating these three modes using fitting of a polynomial and field survey control point,
In, wait the mode actual effect of precise orbit data to be released poor, and some satellites cannot provide precise orbit；Fitting of a polynomial
Mode be to select some stable points in stability region, using one trend surface of phase-fitting on stable point and deduct, the party
Formula is not rigorous and very strong to the dependence of operator's experience；The mode calibrated of actual measurement control point have it is strict theoretical and
Higher precision, but workload is very big.
Invention content
In view of this, the embodiment of the present invention be designed to provide a kind of synthetic aperture radar orbit error removing method and
Device, to improve the above problem.
Present pre-ferred embodiments provide a kind of synthetic aperture radar orbit error removing method, the method includes：
Data of synthetic aperture radar corresponding with target area is obtained, and is generated really according to the data of synthetic aperture radar
SAR images, and include according to the data of synthetic aperture radar the first orbit observation value, the second orbit observation value and radar in
Frequency of heart intercepts dem data corresponding with the target area from digital elevation model；
Over-sampling is carried out to the dem data, the dem data after over-sampling is mapped into SAR coordinate systems in image to establish
Mapping relations between dem data after over-sampling and SAR image coordinates, as initial look-up table；
Simulating image is generated according to the dem data, and the simulating image is registrated with the true SAR images,
Obtain the relativeness between each pixel of the simulating image and the true each pixel of SAR images；
The initial look-up table is verified to obtain look-up table of refining according to the relativeness；
Based on the look-up table of refining, the dem data after over-sampling is mapped into SAR coordinate systems in image to obtain
The corresponding height value of each pixel in true SAR images is stated, and according to the corresponding elevation of each pixel in the true SAR images
Value calculates the corresponding emulation interferometric phase of each pixel；
It is generated according to the first haplopia complex data for including in the data of synthetic aperture radar and the second haplopia complex data dry
Bar graph is related to, difference processing is carried out to obtain difference to the interference fringe picture and the corresponding emulation interferometric phase of each pixel
Divide interferometric phase image；
The differential interferometry phase diagram is expanded, and FFT transform, meter are carried out to the differential interferometry phase diagram after expansion
The intensity value of each pixel after FFT transform is calculated, and circumferential displacement is carried out centered on the maximum pixel position of intensity value；
FFT inverse transformations are carried out to obtain FFT phase diagrams, according to described to the differential interferometry phase diagram after completion circumferential displacement
The original size of interferometric phase image cuts the FFT phase diagrams dry to obtain the difference of the target area after error concealment
Relate to phase.
Further, the relativeness can be obtained by following steps：
The simulating image is registrated with the true SAR images based on correlation function algorithm；
Using quadratic polynomial to registration result carry out orientation and distance to fitting with calculate the simulating image with
Offset between the true SAR images；
The relativeness between the simulating image and the true SAR images is obtained based on the offset.
Further, the emulation interferometric phase is calculated by the following formula to obtain：
Wherein, Q is antenna model, S1For first antenna phase centre location coordinate, P is goal pels coordinate, B second
Relative position coordinates of the antenna relative to first antenna, λ are synthetic aperture radar wavelength.
Further, the step of expanding the differential interferometry phase diagram include：
The width and height of the differential interferometry phase diagram are expanded according to 2 integral number power, so as to complete after expanding
The width of differential interferometry phase diagram and high satisfaction：
Wherein, WinitFor original width, HinitFor original height, WFFTFor the width after expansion, HFFTTo expand further,
Before the step of carrying out circumferential displacement centered on pixel position maximum by intensity value, the method further includes：
Judge whether the position coordinates of the pixel position are integer, if it is not, then using s i nc functions and according to described
The position coordinates of pixel position carry out resampling to data so that the position coordinates are integer.
Present pre-ferred embodiments also provide a kind of synthetic aperture radar orbit error cancellation element, and described device includes：
Dem data generation module, for obtaining data of synthetic aperture radar corresponding with target area, and according to the synthesis
Aperture radar data generates true SAR images, and include according to the data of synthetic aperture radar the first orbit observation value, the
Two orbit observation values and radar center frequency intercept dem data corresponding with the target area from digital elevation model；
Initial look-up table establishes module, for carrying out over-sampling to the dem data, by the DEM numbers after over-sampling
According to mapping to SAR coordinate systems in image to establish the mapping relations between the dem data after over-sampling and SAR image coordinates, as
Initial look-up table；
Pixel registration module, for generating simulating image according to the dem data, and by the simulating image with it is described true
SAR images are registrated, and the opposite pass between each pixel of the simulating image and the true each pixel of SAR images is obtained
System；
Look-up table calibration module, for being verified the initial look-up table to be refined according to the relativeness
Look-up table；
Interferometric phase generation module is emulated, for being based on the look-up table of refining, the dem data after over-sampling is reflected
SAR coordinate systems in image is incident upon to obtain the corresponding height value of each pixel in the true SAR images, and according to described true
The corresponding emulation interferometric phase of each pixel of each corresponding height value calculating of pixel in SAR images；
Difference Calculation module, for according to the first haplopia complex data for including in the data of synthetic aperture radar and second
Haplopia complex data generates interference fringe picture, is carried out to the interference fringe picture and the corresponding emulation interferometric phase of each pixel
Difference processing is to obtain differential interferometry phase diagram；
FFT transform module, for expanding the differential interferometry phase diagram, and to the differential interferometry phase after expansion
Figure carries out FFT transform, calculates the intensity value of each pixel after FFT transform, and is with the maximum pixel position of intensity value
Center carries out circumferential displacement；
Phase calibration generation module, for complete circumferential displacement after differential interferometry phase diagram carry out FFT inverse transformations with
FFT phase diagrams are obtained, the FFT phase diagrams are cut according to the original size of the interferometric phase image and are disappeared with obtaining error
The differential interferometry phase of target area after removing.
Further, the pixel registration module includes：
Pixel registration unit matches the simulating image and the true SAR images for being based on correlation function algorithm
It is accurate；
Offset computation unit, for carrying out orientation and distance to being fitted to registration result using quadratic polynomial
To calculate the offset between the simulating image and the true SAR images；
Relation determination unit is obtained for being based on the offset between the simulating image and the true SAR images
Relativeness.
Further, the emulation interferometric phase is calculated by the following formula to obtain：
Wherein, Q is antenna model, S1For first antenna phase centre location coordinate, P is goal pels coordinate, B second
Relative position coordinates of the antenna relative to first antenna, λ are synthetic aperture radar wavelength.
Further, described device further includes：
Coordinate determination module, for judging whether the position coordinates of the pixel position are integer, if it is not, then using s i
Nc functions simultaneously carry out resampling so that the position coordinates are integer according to the position coordinates of the pixel position to data.
Present pre-ferred embodiments also provide a kind of computer readable storage medium, are stored thereon with computer program, institute
It states computer program and realizes step in above-mentioned synthetic aperture radar orbit error removing method when being executed.
Compared with prior art, a kind of synthetic aperture radar orbit error removing method and dress provided in an embodiment of the present invention
It sets, this method is supported without other data, is only to carry out spectrum estimation based on data of synthetic aperture radar itself and then disappear
Except orbit error, the demand in the prior art to smart rail data and control point, and data processing of the present invention essence can be effectively avoided
Degree is high, has a wide range of application.
To enable the above objects, features and advantages of the present invention to be clearer and more comprehensible, preferred embodiment cited below particularly, and coordinate
Appended attached drawing, is described in detail below.
Description of the drawings
In order to illustrate the technical solution of the embodiments of the present invention more clearly, below will be to needed in the embodiment attached
Figure is briefly described, it should be understood that the following drawings illustrates only certain embodiments of the present invention, therefore is not construed as pair
The restriction of range for those of ordinary skill in the art without creative efforts, can also be according to this
A little attached drawings obtain other relevant attached drawings.
Fig. 1 is the frame structure schematic diagram of electric terminal provided in an embodiment of the present invention.
Fig. 2 is the flow diagram of synthetic aperture radar orbit error removing method provided in an embodiment of the present invention.
Fig. 3 a- Fig. 3 c are respectively the schematic diagram of oblique distance DEM, interference fringe picture, differential interferometry figure.
Fig. 3 d- Fig. 3 f are respectively that differential interferometry figure after expanding, FFT transform result and picture element intensities query result show
It is intended to.
Fig. 3 g and Fig. 3 h are circular shifting result schematic diagram.
Fig. 3 i and Fig. 3 j are that FFT inverse transformed results and orbit error eliminate result schematic diagram.
Fig. 4 is the frame structure schematic diagram of synthetic aperture radar orbit error cancellation element provided in an embodiment of the present invention.
Fig. 5 is the frame structure schematic diagram of pixel registration module shown in Fig. 4.
Icon：10- electric terminals；100- synthetic aperture radar orbit error cancellation elements；110-DEM data generate mould
Block；120- initial look-up tables establish module；130- pixel registration modules；131- pixel registration units；132- offsets calculate single
Member；133- relation determination units；140- look-up table calibration modules；150- emulates interferometric phase generation module；160- Difference Calculations
Module；170-FFT conversion modules；180- phase calibration generation modules；200- memories；300- storage controls；400- processing
Device.
Specific implementation mode
Following will be combined with the drawings in the embodiments of the present invention, and technical solution in the embodiment of the present invention carries out clear, complete
Site preparation describes, it is clear that described embodiments are only a part of the embodiments of the present invention, instead of all the embodiments.Usually
The component for the embodiment of the present invention being described and illustrated herein in the accompanying drawings can be arranged and be designed with a variety of different configurations.Cause
This, the detailed description of the embodiment of the present invention to providing in the accompanying drawings is not intended to limit claimed invention below
Range, but it is merely representative of the selected embodiment of the present invention.Based on the embodiment of the present invention, those skilled in the art are not doing
The every other embodiment obtained under the premise of going out creative work, shall fall within the protection scope of the present invention.
It should be noted that：Similar label and letter indicate similar terms in following attached drawing, therefore, once a certain Xiang Yi
It is defined, then it further need not be defined and explained in subsequent attached drawing in a attached drawing.
As shown in Figure 1, for the frame structure schematic diagram of electric terminal 10 provided in an embodiment of the present invention, the electric terminal 10
Including synthetic aperture radar orbit error cancellation element 100, memory 200, storage control 300 and processor 400.Its
In, the memory 200, storage control 300,400 each element of processor are directly or indirectly electrically connected between each other, with
Realize the transmission or interaction of data.For example, being realized electrically by one or more communication bus or signal wire between these elements
Connection.The synthetic aperture radar orbit error cancellation element 100 can be stored including at least one in the form of software or firmware
In the memory 200 or the software function module that is solidificated in the operating system of the electric terminal 10.The processor
400 access the memory 200 under the control of the storage control 300, are stored for executing in the memory 200
Executable module, such as software function module and meter included by the synthetic aperture radar orbit error cancellation element 100
Calculation machine program etc..Optionally, the electric terminal 10 may be, but not limited to, smart mobile phone, PC (personal
Computer, PC), tablet computer, personal digital assistant (personal digital assistant, PDA), mobile Internet access set
Standby (mobile Internet device, MID) etc..
It is appreciated that structure shown in FIG. 1 is only to illustrate, the electric terminal 10 may also include more than shown in Fig. 1
Either less component or with the configuration different from shown in Fig. 1.Hardware, software may be used in each component shown in Fig. 1
Or combinations thereof realize.
Fig. 2 is please referred to, is a kind of synthetic aperture applied to the electric terminal 10 provided in an embodiment of the present invention
The flow diagram of radar orbit error removing method, method and step defined in the related flow of the method can be by described
Processor 400 is realized.Detailed process shown in Fig. 2 will be described in detail below.
Before synthetic aperture radar orbit error removing method is described, first to the synthetic aperture radar into
Row is briefly introduced：Synthetic aperture radar (Synthetic Aperture Radar, SAR) is using a miniature antenna along long line
The track constant speed movement of battle array simultaneously radiates coherent signal, and the echo that is received in different location is carried out Coherent processing, to obtain compared with
High-resolution imaging radar.It should be understood that the synthetic aperture radar is a single aerial system or dual-antenna system, wherein
A single aerial system observes target area along same track twice, and dual-antenna system is same to target area by two sets of antennas
When be observed.For a single aerial system, the first antenna being previously mentioned in the present embodiment is antenna when observing for the first time, institute
Antenna when the second antenna is second of observation is stated, for dual-antenna system, two sets of antennas are divided into primary antenna and slave antenna, institute
It is primary antenna to state first antenna, and second antenna is slave antenna；The main image, the first haplopia complex data, the first track are seen
Measured value is obtained based on the first antenna, and the subpictures, the second haplopia complex data, the second orbit observation value are i.e. based on described
Second antenna acquisition etc., details are not described herein for the present embodiment.Disappear to synthetic aperture radar orbit error with reference to specific steps
Except method is introduced.
Step S10 obtains data of synthetic aperture radar, and generates true SAR images according to the data of synthetic aperture radar,
And the first orbit observation value, the second orbit observation value and the radar center frequency for according to the data of synthetic aperture radar including from
Dem data corresponding with the target area is intercepted in digital elevation model.
In the present embodiment, the data of synthetic aperture radar is target area space flight synthetic aperture radar haplopia complex data
(such as the first haplopia complex data, the second haplopia complex data), the first orbit observation value, the first orbit observation value, radar center frequency,
Initial oblique distance, Doppler's geometric parameter etc..Wherein, the detailed process for intercepting the dem data includes：Track first is read to observe
Value, initial oblique distance and Doppler frequency calculate four angle point geographical coordinates of main image using Range-Doppler model, according to
Four angle point geographical coordinates of main image count coordinate range, and then are cut from DEM digital elevation models according to scope of statistics
Go out the dem data of corresponding target area.
Step S20 carries out over-sampling to the dem data, the dem data after over-sampling is mapped to SAR image coordinates
System is to establish the mapping relations between the dem data after over-sampling and SAR image coordinates, as initial look-up table.
In the present embodiment, when carrying out over-sampling to the dem data, adopting less than 2 times of SAR image resolutions can be used
Sample interval carries out data sampling.In addition, when generating initial look-up table, each in the dem data after over-sampling need to be calculated
The three-dimensional coordinate of grid points, and its image coordinate on the true SAR images is calculated by indirect location, and then obtain
The mapping relations between dem data and SAR image coordinates after to the over-sampling.
Step S30 generates simulating image according to the dem data, and by the simulating image and the true SAR images into
Row registration, obtains the relativeness between each pixel of the simulating image and the true each pixel of SAR images.
Step S40 verifies to obtain look-up table of refining the initial look-up table according to the relativeness.
In step S30 and step S40, when generating the simulating image, each grid points of DEM after over-sampling need to be calculated
Local incidence angle, i.e. the angle of antenna and grid points line and landform tangent plane residing for the grid points takes the cosine value of angle
As the simulating image.Correlation function algorithm is then based on to be registrated the simulating image with the true SAR images；Profit
Orientation and distance are carried out to being fitted to registration result with quadratic polynomial with calculate the simulating image with it is described true
Offset between SAR images；It is obtained based on the offset opposite between the simulating image and the true SAR images
Relationship.
The dem data after over-sampling is mapped to SAR coordinate systems in image by step S50 based on the look-up table of refining
To obtain the corresponding height value of each pixel in the true SAR images, and according to each pixel pair in the true SAR images
The height value answered calculates the corresponding emulation interferometric phase of each pixel.
Wherein, the corresponding height value of each pixel is dem data under coordinate systems in image in the true SAR images,
Alternatively referred to as oblique distance DEM (as shown in Figure 3a).In addition, in the present embodiment, it is corresponding each pixel can be calculated by following procedure
Interferometric phase is emulated, it is specific as follows.
(1) baseline estimations.Read the corresponding first orbit observation value of major and minor image and the second orbit observation value.And
In the time range for obtaining main image, 3 moment can be selected at equal intervals, go out each moment using Neville difference approach interpolations
Primary antenna positionTrack to the nearest position of primary antenna where calculating separately out slave antennaCalculate the opposite position of the two
It setsFinally use fitting of a polynomial relative position BiIt can be obtained baseline B.
(2) three-dimensional coordinate of a certain goal pels in the true SAR images is calculated.Read main image corresponding first
The height value of each pixel in orbit observation value, initial oblique distance, Doppler frequency and the true SAR images utilizes distance-
The three-dimensional coordinate of Doppler model and each pixel of ellipsoid equation calculation is shown below.
In formula, r1For the distance of first antenna to goal pels, can be obtained by initial oblique distance and sampling interval；S1It is
One antenna phase center, can be by obtaining the first orbit observation value interpolation；P=[Px Py Pz]TFor goal pels coordinate；fdop
For Doppler frequency；V velocity vectors can calculate the first derivative of time by the first orbit observation value and obtain；H is described true
The corresponding height value of pixel in real SAR images；RaReference ellipsoid long axis；RbFor reference ellipsoid short axle；λ is synthetic aperture radar
Wavelength.
(3) computer sim- ulation interferometric phase.The emulation interferometric phaseIt can be calculated by the following formula to obtain：
Wherein, Q is antenna model, and e.g., Q=1 indicates that the single-emission and double-receiving pattern of double antenna interference system, Q=2 indicate Dan Tian
The repeat track interference pattern of linear system system or " table tennis " pattern of double antenna interference system.S1For first antenna (primary antenna) phase
Position center position coordinates, P are goal pels coordinate, and B is relative position coordinates of second antenna relative to first antenna, and λ is to close
At aperture radar wavelength.
Step S60, according to the first haplopia complex data for including in the data of synthetic aperture radar and the second haplopia plural number
According to interference fringe picture is generated, difference processing is carried out to the interference fringe picture and the corresponding emulation interferometric phase of each pixel
To obtain differential interferometry phase diagram.
Wherein, the interference fringe picture to be generated in this implementation as shown in Figure 3b.When calculating differential interferometry phase diagram, can incite somebody to action
Emulation interferometric phase interference fringe corresponding with the interference fringe picture carries out complex conjugate multiplication, obtains differential interferometry phase
Position, can be as shown in Figure 3c.
Step S70 expands the differential interferometry phase diagram, and is carried out to the differential interferometry phase diagram after expansion
FFT transform, calculates the intensity value of each pixel after FFT transform, and centered on the maximum pixel position of intensity value into
Row circumferential displacement.
In the present embodiment, the step of expanding the differential interferometry phase diagram, includes：According to 2 integral number power to institute
The width and height for stating differential interferometry phase diagram are expanded, and are met so as to complete the wide and high of the differential interferometry phase diagram after expandingWherein, WinitFor original width, HinitFor original height, WFFTFor the width after expansion
Degree, HFFTIt is as shown in Figure 3d the differential interferometry figure after expansion for the height for expanding after adjusting.During actual implementation, work as WFFT>
WinitOr HFFT>HinitWhen, additional region be placed on initial data (the differential interferometry phase diagram before not expanding) right side and
Downside, and with 0 filling.
Further, it is as shown in Figure 3 e the transformation results in step S70 after FFT transform, and then is based on the calculating knot
Fruit can calculate the intensity value of each pixel after FFT transform, and be justified centered on the maximum pixel position of intensity value
Zhou Weiyi is as follows：
(1) assume that the intensity for calculating pixel I isIn formula, IpwrFor the intensity value of pixel, IreFor
The real part of pixel, IimFor the imaginary part of pixel, then, the maximum pixel of the intensity value of the pixel found can be remembered as illustrated in figure 3f
For (imax jmax)。
(2) with position (imax jmax) centered on, it takes 3 × 3 windows and the intensity value of each pixel is carried out with quadratic polynomial
Fitting obtains f=a0+a1x+a2y+a3x2+a4y2, wherein a0、a1、a2、a3And a4For undetermined coefficient, using least square method
Solve undetermined coefficient a0、a1、a2、a3And a4, and then obtain position (x when f is maximum valuemax ymax)。
(3) by position (xmax ymax) circular shifting to position the first row first row, such as Fig. 3 g and 3h show circumferential displacement
As a result.It should be noted here that before carrying out circular shifting, the position coordinates x for judging the pixel position is neededmaxOr ymaxWhether
For integer, if it is not, then using sinc functions and the position coordinates x according to the pixel positionmaxOr ymaxData are adopted again
Sample is so that position coordinates xmaxOr ymaxCircumferential displacement is carried out again for integer, wherein the sinc functions are two dimension normalization letter
Number.
Step S80, to the differential interferometry phase diagram progress FFT inverse transformations after completion circumferential displacement to obtain FFT phase diagrams,
The FFT phase diagrams are cut to obtain the target area after error concealment according to the original size of the interferometric phase image
Differential interferometry phase.
Wherein, FFT inverse transformations are carried out to the data after resampling, is as shown in fig. 2i the transformation knot after FFT inverse transformations
Fruit schematic diagram can be cut to obtain the synthetic aperture as shown in Fig. 3 j based on the transformation results according to original size and position
SAR differential interferometry phase schematic diagrames after radar error concealment.
In above-mentioned synthetic aperture radar orbit error removing method can clearly it can be seen from the embodiment of the present invention
It is that spectrum estimation is carried out based on data of synthetic aperture radar itself and main frequency is zeroed, it is theoretical tight, and do not need other
Data are supported, and can effectively improve data processing precision.In addition, the synthetic aperture radar orbit error that the present invention provides disappears
Except method can play an important role in differential interferometry synthetic aperture radar obtains Ground Deformation information processing.
Further, based on the description to above-mentioned synthetic aperture radar orbit error removing method, the present invention is preferably implemented
Example also provides a kind of computer readable storage medium, is stored thereon with computer program, and the computer program is real when being executed
Step in existing above-mentioned synthetic aperture radar orbit error removing method.
Further, as shown in figure 4, being synthetic aperture radar orbit error cancellation element provided in an embodiment of the present invention
100 frame structure schematic diagram, the synthetic aperture radar orbit error cancellation element 100 are applied to electric terminal 10, wherein
The synthetic aperture radar orbit error cancellation element 100 establishes module including dem data generation module 110, initial look-up table
120, pixel registration module 130, look-up table calibration module 140, emulation interferometric phase generation module 150, Difference Calculation module
160, FFT transform module 170 and phase calibration generation module 180.
The dem data generation module 110, for obtaining data of synthetic aperture radar corresponding with target area, and root
True SAR images are generated according to the data of synthetic aperture radar, and the first track for including according to the data of synthetic aperture radar
Observation, the second orbit observation value and radar center frequency intercept corresponding with the target area from digital elevation model
Dem data；In the present embodiment, the description as described in the dem data generation module 110 specifically refers to the detailed of above-mentioned steps S10
Thin description, that is, the step S10 can be executed by dem data generation module 110, thus does not illustrate more herein.
The initial look-up table establishes module 120, for carrying out over-sampling to the dem data, by the institute after over-sampling
It states dem data and maps to SAR coordinate systems in image to establish the pass of the mapping between the dem data after over-sampling and SAR image coordinates
System, as initial look-up table；In the present embodiment, the description as described in the initial look-up table establishes module 120 specifically refers to
The detailed description of step S20 is stated, is executed that is, the step S20 can be established module 120 by initial look-up table, thus herein
Do not illustrate more.
The pixel registration module 130, for according to the dem data generate simulating image, and by the simulating image with
The true SAR images are registrated, and are obtained between each pixel of the simulating image and the true each pixel of SAR images
Relativeness；In the present embodiment, the description as described in the pixel registration module 130 specifically refers to the detailed of above-mentioned steps S30
Thin description, that is, the step S30 can be executed by pixel registration module 130, thus does not illustrate more herein.Optionally,
It is determined as shown in figure 5, the pixel registration module 130 includes pixel registration unit 131, offset computation unit 132 and relationship
Unit 133.
The pixel registration unit 131, for being based on correlation function algorithm to the simulating image and the true SAR images
It is registrated；
The offset computation unit 132, for using quadratic polynomial to registration result progress orientation and distance to
It is fitted to calculate the offset between the simulating image and the true SAR images；
The relation determination unit 133 obtains the simulating image and the true SAR shadows for being based on the offset
Relativeness as between.
The look-up table calibration module 140, for according to the relativeness to the initial look-up table verified with
Obtain look-up table of refining；In the present embodiment, the description as described in the look-up table calibration module 140 specifically refers to above-mentioned steps
The detailed description of S40 that is, the step S40 can be executed by look-up table calibration module 140, thus is not spoken more more herein
It is bright.
The emulation interferometric phase generation module 150, for being based on the look-up table of refining, described in after over-sampling
Dem data maps to SAR coordinate systems in image to obtain the corresponding height value of each pixel in the true SAR images, and according to
The corresponding emulation interferometric phase of each pixel of each corresponding height value calculating of pixel in the true SAR images；The present embodiment
In, the description as described in the emulation interferometric phase generation module 150 specifically refers to the detailed description of above-mentioned steps S50, that is,
The step S50 can be executed by emulation interferometric phase generation module 150, thus not illustrated more herein.
The Difference Calculation module 160, for according to the first haplopia plural number for including in the data of synthetic aperture radar
Interference fringe picture is generated according to the second haplopia complex data, the interference fringe picture and the corresponding emulation of each pixel are interfered
Phase carries out difference processing to obtain differential interferometry phase diagram；In the present embodiment, the description as described in the Difference Calculation module 160
The detailed description of above-mentioned steps S60 is specifically referred to, that is, the step S60 can be executed by Difference Calculation module 160, because
And do not illustrate more herein.
The FFT transform module 170, for expanding the differential interferometry phase diagram, and to the difference after expansion
Interferometric phase image carries out FFT transform, calculates the intensity value of each pixel after FFT transform, and with the maximum picture of intensity value
Circumferential displacement is carried out centered on first position.In the present embodiment, the description as described in the FFT transform module 170 specifically refers to
The detailed description of step S70 is stated, that is, the step S70 can be executed by FFT transform module 170, thus is not made herein more
Explanation.
The phase calibration generation module 180, it is inverse for carrying out FFT to the differential interferometry phase diagram after completion circumferential displacement
Transformation cuts to obtain the FFT phase diagrams according to the original size of the interferometric phase image with obtaining FFT phase diagrams
The differential interferometry phase of target area after error concealment.In the present embodiment, about retouching for the phase calibration generation module 180
The detailed description for specifically referring to above-mentioned steps S80 is stated, that is, the step S80 can be held by phase calibration generation module 180
Row, thus do not illustrate more herein.
In conclusion a kind of synthetic aperture radar orbit error removing method provided in an embodiment of the present invention and device, it should
Method is supported without other data, is only to carry out spectrum estimation based on data of synthetic aperture radar itself and then eliminate track
Error can effectively avoid the demand in the prior art to smart rail data and control point, and data processing precision of the present invention is high, answers
It is wide with range.
In several embodiments that the embodiment of the present invention is provided, it should be understood that disclosed device and method also may be used
To realize by another way.Device and method embodiment described above is only schematical, for example, in attached drawing
Flow chart and block diagram show the device of multiple embodiments according to the present invention, the possibility of method and computer program product is realized
Architecture, function and operation.In this regard, each box in flowchart or block diagram can represent module, a program
A part for a part for section or code, the module, section or code includes that one or more is patrolled for realizing defined
Collect the executable instruction of function.It should also be noted that at some as the function of in the realization method replaced, being marked in box
It can occur in a different order than that indicated in the drawings.For example, two continuous boxes can essentially be held substantially in parallel
Row, they can also be executed in the opposite order sometimes, this is depended on the functions involved.It is also noted that block diagram and/or
The combination of each box in flow chart and the box in block diagram and or flow chart can use function or dynamic as defined in executing
The dedicated hardware based system made is realized, or can be realized using a combination of dedicated hardware and computer instructions.
In addition, each function module in each embodiment of the present invention can integrate to form an independent portion
Point, can also be modules individualism, can also two or more modules be integrated to form an independent part.
It, can be with if the function is realized and when sold or used as an independent product in the form of software function module
It is stored in a computer read/write memory medium.Based on this understanding, technical scheme of the present invention is substantially in other words
The part of the part that contributes to existing technology or the technical solution can be expressed in the form of software products, the meter
Calculation machine software product is stored in a storage medium, including some instructions are used so that a computer equipment (can be
People's computer, electronic equipment or network equipment etc.) execute all or part of step of each embodiment the method for the present invention
Suddenly.And storage medium above-mentioned includes：USB flash disk, read-only memory (ROM, Read-Only Memory), is deposited mobile hard disk at random
The various media that can store program code such as access to memory (RAM, Random Access Memory), magnetic disc or CD.
It should be noted that herein, the terms "include", "comprise" or its any other variant are intended to the packet of nonexcludability
Contain, so that the process, method, article or equipment including a series of elements includes not only those elements, but also includes
Other elements that are not explicitly listed, or further include for elements inherent to such a process, method, article, or device.
In the absence of more restrictions, the element limited by sentence "including a ...", it is not excluded that including the element
Process, method, article or equipment in there is also other identical elements.
The foregoing is merely the alternative embodiments of the present invention, are not intended to restrict the invention, for the skill of this field
For art personnel, the invention may be variously modified and varied.All within the spirits and principles of the present invention, any made by repair
Change, equivalent replacement, improvement etc., should all be included in the protection scope of the present invention.
Claims (10)
1. a kind of synthetic aperture radar orbit error removing method, which is characterized in that the method includes：
Data of synthetic aperture radar corresponding with target area is obtained, and true SAR is generated according to the data of synthetic aperture radar
Image, and the first orbit observation value, the second orbit observation value and the radar center that include according to the data of synthetic aperture radar
Frequency intercepts dem data corresponding with the target area from digital elevation model；
Over-sampling is carried out to the dem data, the dem data after over-sampling is mapped into SAR coordinate systems in image and was adopted with establishing
Mapping relations between dem data after sample and SAR image coordinates, as initial look-up table；
Simulating image is generated according to the dem data, and the simulating image is registrated with the true SAR images, is obtained
Relativeness between each pixel of simulating image and the true each pixel of SAR images；
The initial look-up table is verified to obtain look-up table of refining according to the relativeness；
Based on the look-up table of refining, it is described true to obtain that the dem data after over-sampling is mapped to SAR coordinate systems in image
The corresponding height value of each pixel in real SAR images, and according to the corresponding height value meter of each pixel in the true SAR images
Calculate the corresponding emulation interferometric phase of each pixel；
Interference item is generated according to the first haplopia complex data for including in the data of synthetic aperture radar and the second haplopia complex data
It is dry to obtain difference to carry out difference processing to the interference fringe picture and the corresponding emulation interferometric phase of each pixel for line figure
Relate to phase diagram；
The differential interferometry phase diagram is expanded, and FFT transform is carried out to the differential interferometry phase diagram after expansion, calculates warp
The intensity value of each pixel after FFT transform, and carry out circumferential displacement centered on the maximum pixel position of intensity value；
FFT inverse transformations are carried out to obtain FFT phase diagrams, according to the interference to the differential interferometry phase diagram after completion circumferential displacement
The original size of phase diagram cuts the FFT phase diagrams to obtain the differential interferometry phase of the target area after error concealment
Position.
2. synthetic aperture radar orbit error removing method according to claim 1, which is characterized in that the relativeness
It can be obtained by following steps：
The simulating image is registrated with the true SAR images based on correlation function algorithm；
Using quadratic polynomial to registration result carry out orientation and distance to fitting with calculate the simulating image with it is described
Offset between true SAR images；
The relativeness between the simulating image and the true SAR images is obtained based on the offset.
3. synthetic aperture radar orbit error removing method according to claim 1, which is characterized in that the emulation interference
Phase is calculated by the following formula to obtain：
Wherein, Q is antenna model, S1For first antenna phase centre location coordinate, P is goal pels coordinate, and B is the second antenna
Relative to the relative position coordinates of first antenna, λ is synthetic aperture radar wavelength.
4. synthetic aperture radar orbit error removing method according to claim 1, which is characterized in that dry to the difference
Relating to the step of phase diagram is expanded includes：
The width and height of the differential interferometry phase diagram are expanded according to 2 integral number power, so as to complete the difference after expanding
The width of interferometric phase image and high satisfaction：
Wherein, WinitFor original width, HinitFor original height, WFFTFor the width after expansion, HFFTFor the height for expanding after adjusting
Degree.
5. synthetic aperture radar orbit error removing method according to claim 1, which is characterized in that with intensity value most
Before the step of carrying out circumferential displacement centered on big pixel position, the method further includes：
Judge whether the position coordinates of the pixel position are integer, if it is not, then using sinc functions and according to the pixel position
The position coordinates set carry out resampling to data so that the position coordinates are integer.
6. a kind of synthetic aperture radar orbit error cancellation element, which is characterized in that described device includes：
Dem data generation module, for obtaining data of synthetic aperture radar corresponding with target area, and according to the synthetic aperture
Radar data generates true SAR images, and include according to the data of synthetic aperture radar the first orbit observation value, the second rail
Road observation and radar center frequency intercept dem data corresponding with the target area from digital elevation model；
Initial look-up table establishes module, and for carrying out over-sampling to the dem data, the dem data after over-sampling is mapped to
SAR coordinate systems in image is to establish the mapping relations between the dem data after over-sampling and SAR image coordinates, as initial lookup
Table；
Pixel registration module, for generating simulating image according to the dem data, and by the simulating image and the true SAR
Image is registrated, and the relativeness between each pixel of the simulating image and the true each pixel of SAR images is obtained；
Look-up table calibration module, for being verified to the initial look-up table with the lookup that obtains refining according to the relativeness
Table；
Interferometric phase generation module is emulated, for being based on the look-up table of refining, the dem data after over-sampling is mapped to
SAR coordinate systems in image is to obtain the corresponding height value of each pixel in the true SAR images, and according to the true SAR shadows
The corresponding emulation interferometric phase of each pixel of each corresponding height value calculating of pixel as in；
Difference Calculation module, for according to the first haplopia complex data for including in the data of synthetic aperture radar and the second haplopia
Complex data generates interference fringe picture, and difference is carried out to the interference fringe picture and the corresponding emulation interferometric phase of each pixel
Processing is to obtain differential interferometry phase diagram；
FFT transform module, for expanding the differential interferometry phase diagram, and to the differential interferometry phase diagram after expansion into
Row FFT transform calculates the intensity value of each pixel after FFT transform, and centered on the maximum pixel position of intensity value
Carry out circumferential displacement；
Phase calibration generation module, for carrying out FFT inverse transformations to the differential interferometry phase diagram after completion circumferential displacement to obtain
FFT phase diagrams, after being cut the FFT phase diagrams to obtain error concealment according to the original size of the interferometric phase image
Target area differential interferometry phase.
7. synthetic aperture radar orbit error cancellation element according to claim 6, which is characterized in that the pixel registration
Module includes：
Pixel registration unit is registrated the simulating image with the true SAR images for being based on correlation function algorithm；
Offset computation unit, for using quadratic polynomial to registration result carry out orientation and distance to fitting to calculate
Offset between the simulating image and the true SAR images；
Relation determination unit, for obtaining the phase between the simulating image and the true SAR images based on the offset
To relationship.
8. synthetic aperture radar orbit error cancellation element according to claim 6, which is characterized in that the emulation interference
Phase is calculated by the following formula to obtain：
Wherein, Q is antenna model, S1For first antenna phase centre location coordinate, P is goal pels coordinate, and B is the second antenna
Relative to the relative position coordinates of first antenna, λ is synthetic aperture radar wavelength.
9. synthetic aperture radar orbit error cancellation element according to claim 6, which is characterized in that described device is also wrapped
It includes：
Coordinate determination module, for judging whether the position coordinates of the pixel position are integer, if it is not, then using sinc functions
And resampling is carried out to data so that the position coordinates are integer according to the position coordinates of the pixel position.
10. a kind of computer readable storage medium, is stored thereon with computer program, which is characterized in that the computer program
The step in the synthetic aperture radar orbit error removing method described in any one of the claims 1-5 is realized when being executed
Suddenly.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810207673.2A CN108445458B (en) | 2018-03-14 | 2018-03-14 | Synthetic aperture radar track error elimination method and device |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810207673.2A CN108445458B (en) | 2018-03-14 | 2018-03-14 | Synthetic aperture radar track error elimination method and device |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN108445458A true CN108445458A (en) | 2018-08-24 |
| CN108445458B CN108445458B (en) | 2020-08-14 |
Family
ID=63194271
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201810207673.2A Active CN108445458B (en) | 2018-03-14 | 2018-03-14 | Synthetic aperture radar track error elimination method and device |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN108445458B (en) |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110501705A (en) * | 2019-07-29 | 2019-11-26 | 武汉大学 | The basin ground SAR library bank deformation monitoring device based on long-distance radio communication |
| CN112596055A (en) * | 2020-12-08 | 2021-04-02 | 中国地质大学（武汉） | Method for correcting residual system error of InSAR DEM |
| CN118011344A (en) * | 2024-04-02 | 2024-05-10 | 中国科学院空天信息创新研究院 | SAR interference calibration method assisted by external DEM |
Citations (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102520406A (en) * | 2011-12-20 | 2012-06-27 | 中国测绘科学研究院 | SAR (synthetic aperture radar) stereoscopic model construction and measurement method based on image coordinate mapping |
| CN103091676A (en) * | 2013-01-22 | 2013-05-08 | 中国矿业大学 | Mining area surface subsidence synthetic aperture radar interferometry monitoring and calculating method |
| US20140062764A1 (en) * | 2012-09-04 | 2014-03-06 | Fugro Earthdata, Inc. | Method and apparatus for mapping and characterizing sea ice from airborne simultaneous dual frequency interferometric synthetic aperture radar (ifsar) measurements |
| CN103675790A (en) * | 2013-12-23 | 2014-03-26 | 中国国土资源航空物探遥感中心 | Method for improving earth surface shape change monitoring precision of InSAR (Interferometric Synthetic Aperture Radar) technology based on high-precision DEM (Digital Elevation Model) |
| CN103728604A (en) * | 2013-11-19 | 2014-04-16 | 中国国土资源航空物探遥感中心 | Broadband synthetic aperture radar sub-band interferometric data processing method |
| CN105974411A (en) * | 2016-04-28 | 2016-09-28 | 电子科技大学 | DINSAR (differential interferometric synthetic aperture radar)-based high-voltage transmission tower top tilt displacement monitoring method |
| CN106908770A (en) * | 2017-01-24 | 2017-06-30 | 上海卫星工程研究所 | The ground integrated emulation mode of high-resolution microwave imaging satellite star |
| WO2018027332A1 (en) * | 2016-08-08 | 2018-02-15 | Comercial E Industrial Gesecology Limitada | Method and system for the analysis and generation of early or predictive alerts concerning the stability of slopes in open-pit mines |
-
2018
- 2018-03-14 CN CN201810207673.2A patent/CN108445458B/en active Active
Patent Citations (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102520406A (en) * | 2011-12-20 | 2012-06-27 | 中国测绘科学研究院 | SAR (synthetic aperture radar) stereoscopic model construction and measurement method based on image coordinate mapping |
| US20140062764A1 (en) * | 2012-09-04 | 2014-03-06 | Fugro Earthdata, Inc. | Method and apparatus for mapping and characterizing sea ice from airborne simultaneous dual frequency interferometric synthetic aperture radar (ifsar) measurements |
| CN103091676A (en) * | 2013-01-22 | 2013-05-08 | 中国矿业大学 | Mining area surface subsidence synthetic aperture radar interferometry monitoring and calculating method |
| CN103728604A (en) * | 2013-11-19 | 2014-04-16 | 中国国土资源航空物探遥感中心 | Broadband synthetic aperture radar sub-band interferometric data processing method |
| CN103675790A (en) * | 2013-12-23 | 2014-03-26 | 中国国土资源航空物探遥感中心 | Method for improving earth surface shape change monitoring precision of InSAR (Interferometric Synthetic Aperture Radar) technology based on high-precision DEM (Digital Elevation Model) |
| CN105974411A (en) * | 2016-04-28 | 2016-09-28 | 电子科技大学 | DINSAR (differential interferometric synthetic aperture radar)-based high-voltage transmission tower top tilt displacement monitoring method |
| WO2018027332A1 (en) * | 2016-08-08 | 2018-02-15 | Comercial E Industrial Gesecology Limitada | Method and system for the analysis and generation of early or predictive alerts concerning the stability of slopes in open-pit mines |
| CN106908770A (en) * | 2017-01-24 | 2017-06-30 | 上海卫星工程研究所 | The ground integrated emulation mode of high-resolution microwave imaging satellite star |
Non-Patent Citations (2)
| Title |
|---|
| ZHANG KUI等: ""Interferometric phase reconstruction using simplified coherence network"", 《ISPRS JOURNAL OF PHOTOGRAMMETRY AND REMOTE SENSING》 * |
| 王萌萌等: ""基于点相位改正的机载干涉SAR定标算法"", 《测绘学报》 * |
Cited By (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110501705A (en) * | 2019-07-29 | 2019-11-26 | 武汉大学 | The basin ground SAR library bank deformation monitoring device based on long-distance radio communication |
| CN112596055A (en) * | 2020-12-08 | 2021-04-02 | 中国地质大学（武汉） | Method for correcting residual system error of InSAR DEM |
| CN112596055B (en) * | 2020-12-08 | 2023-04-25 | 中国地质大学（武汉） | Method for correcting residual systematic error of InSAR DEM |
| CN118011344A (en) * | 2024-04-02 | 2024-05-10 | 中国科学院空天信息创新研究院 | SAR interference calibration method assisted by external DEM |
| CN118011344B (en) * | 2024-04-02 | 2024-06-11 | 中国科学院空天信息创新研究院 | SAR interference calibration method assisted by external DEM |
Also Published As
| Publication number | Publication date |
|---|---|
| CN108445458B (en) | 2020-08-14 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US8466834B2 (en) | Radar-imaging of a scene in the far-field of a one- or two-dimensional radar array | |
| CN111856461B (en) | Improved PFA-based bunching SAR imaging method and DSP implementation thereof | |
| CN108445458A (en) | A kind of synthetic aperture radar orbit error removing method and device | |
| CN104007440A (en) | Imaging method for acceleratedly factorized back-projection bunching synthetic aperture radar | |
| CN110865346B (en) | Satellite-borne SAR time parameter calibration method based on direct positioning algorithm | |
| AU2018206753A1 (en) | Method and device for filling invalid regions of terrain elevation model data | |
| CN107219527A (en) | A kind of single snap fast imaging method of the double array channel formula imaging systems of preiodic type | |
| CN112711021B (en) | Multi-resolution InSAR (interferometric synthetic Aperture Radar) interactive interference time sequence analysis method | |
| US20040004569A1 (en) | Method and system for mutual coherent synthetic aperture radiometry | |
| KR102151362B1 (en) | Image decoding apparatus based on airborn using polar coordinates transformation and method of decoding image using the same | |
| CN107390216A (en) | High speed super-resolution stationary point scan imaging method based on wave-number domain coherence factor | |
| RU2661941C1 (en) | Method for radar imaging using the radar station with the synthesized aperture | |
| CN111721241A (en) | GNSS-InBSAR and GB-InSAR cross-system fusion three-dimensional deformation measurement method | |
| US7402794B2 (en) | Radiometer imaging system and method thereof | |
| CN111127334B (en) | SAR image real-time geometric correction method and system based on RD plane pixel mapping | |
| CN109061645B (en) | Virtual correlation sparse imaging method and system for synthetic aperture radiometer | |
| CN115453530A (en) | Bistatic SAR (synthetic aperture radar) filtering back-projection two-dimensional self-focusing method based on parameterized model | |
| CN109541598B (en) | Side lobe ratio evaluation method and device for SAR imaging quality and storage medium | |
| CN108333562A (en) | A kind of dimensionality reduction method for registering images that landform altitude is adaptive | |
| Pipia et al. | A comparison of different techniques for atmospheric artefact compensation in GBSAR differential acquisitions | |
| Pepin | Accuracy of point scatterer localization across a range of 2D and 3D SAR imaging parameters | |
| US12025693B2 (en) | Method for correcting a synthetic aperture radar antenna beam image | |
| CN109581377B (en) | Resolution evaluation method and device for point echo imaging quality and storage medium | |
| Cheung | Digital beamforming array phase calibration techniques for multi-pass interferometric SAR | |
| CN114966679B (en) | Method and system for generating omnibearing DEM (digital elevation model) of low frequency InCSAR of airborne heavy rail |
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 |