CN111965646A - Satellite-borne radar data processing method, device, equipment and storage medium - Google Patents
Satellite-borne radar data processing method, device, equipment and storage medium Download PDFInfo
- Publication number
- CN111965646A CN111965646A CN202011121051.1A CN202011121051A CN111965646A CN 111965646 A CN111965646 A CN 111965646A CN 202011121051 A CN202011121051 A CN 202011121051A CN 111965646 A CN111965646 A CN 111965646A
- Authority
- CN
- China
- Prior art keywords
- interferogram
- phase
- baseline
- calculating
- track
- 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
- 238000003672 processing method Methods 0.000 title claims abstract description 11
- 230000008859 change Effects 0.000 claims abstract description 54
- 238000000034 method Methods 0.000 claims abstract description 30
- 238000012545 processing Methods 0.000 claims abstract description 21
- 238000004364 calculation method Methods 0.000 claims description 23
- 238000012544 monitoring process Methods 0.000 claims description 13
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Chemical compound O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 8
- 238000004804 winding Methods 0.000 claims description 6
- 238000004590 computer program Methods 0.000 claims description 4
- 238000012937 correction Methods 0.000 claims description 4
- 239000000126 substance Substances 0.000 description 10
- 238000010586 diagram Methods 0.000 description 4
- 230000009471 action Effects 0.000 description 3
- 238000004422 calculation algorithm Methods 0.000 description 3
- 238000004891 communication Methods 0.000 description 3
- 230000008878 coupling Effects 0.000 description 3
- 238000010168 coupling process Methods 0.000 description 3
- 238000005859 coupling reaction Methods 0.000 description 3
- 230000000694 effects Effects 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000009897 systematic effect Effects 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 1
- 230000006870 function Effects 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000005305 interferometry Methods 0.000 description 1
- 238000002844 melting Methods 0.000 description 1
- 230000008018 melting Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 239000002689 soil Substances 0.000 description 1
Images
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
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Radar, Positioning & Navigation (AREA)
- Physics & Mathematics (AREA)
- Electromagnetism (AREA)
- Computer Networks & Wireless Communication (AREA)
- General Physics & Mathematics (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
The application provides a method, a device and equipment for processing satellite-borne radar data and a storage medium. The satellite-borne radar data processing method comprises the following steps: acquiring SLC images of the earth surface area to be monitored, which are generated by the satellite-borne radar at different time points, and generating an interference pattern of the earth surface area to be monitored according to the interference pair images; calculating to obtain an initial baseline of the interferogram; removing the first plano-phase and the first topographic phase in the interferogram according to the initial baseline, and obtaining a first differential interferogram of the interferogram; calculating the change frequency of the track residual fringe; correcting an initial baseline of the interference pattern according to the change frequency of the track residual fringe and obtaining a compensated baseline of the interference pattern; removing the orbit error in the interferogram according to the compensated baseline of the interferogram. The method and the device can remove the orbit error in the satellite-borne radar data, and further can improve the inversion precision of the deformation of the earth surface area to be monitored.
Description
Technical Field
The application relates to the field of interferometry, in particular to a method, a device and equipment for processing satellite-borne radar data and a storage medium.
Background
At present, the time sequence InSAR technology based on the long time sequence SAR image is widely applied to the aspects of surface subsidence, landslide, earthquake, fault activity, volcanic activity, frozen soil melting, deformation monitoring of human infrastructure and the like.
The accuracy of the time series InSAR technique is susceptible to a variety of error factors, including atmospheric delay errors, terrain residual errors, orbit errors, incoherent noise, and the like, where orbit errors are particularly significant when monitoring large areas.
Disclosure of Invention
An object of the embodiment of the application is to provide a method, a device, equipment and a storage medium for processing satellite-borne radar data, so as to remove orbit errors in the satellite-borne radar data, and further improve the inversion accuracy of deformation of a to-be-monitored ground surface area.
To this end, the first aspect of the present application discloses a method for processing satellite-borne radar data, the method comprising:
acquiring SLC images generated by the satellite-borne radar at different time points and aiming at an earth surface area to be monitored, wherein the SLC images at every two time points form interference pair images;
generating an interference image of the earth surface area to be monitored according to the interference pair image;
calculating to obtain an initial baseline of the interferogram;
removing a first plano phase and a first topographic phase in the interferogram according to the initial baseline, and obtaining a first differential interferogram of the interferogram;
when the concentration or the significance of the track residual fringes in the first differential interference pattern meets a preset threshold value, calculating the change frequency of the track residual fringes;
correcting the initial baseline of the interference pattern according to the change frequency of the track residual fringe and obtaining a compensated baseline of the interference pattern;
and removing the track error in the interferogram according to the compensated base line of the interferogram.
In the first aspect of the application, a first plano phase and a first topographic phase in an interferogram of an earth surface area to be monitored are removed through an initial baseline, so that a first differential interferogram of the earth surface area to be monitored can be obtained, further, when the density or the significance of track residual fringes in the first differential interferogram meets a preset threshold, the change frequency of the track residual fringes is calculated, further, the initial baseline of the interferogram is corrected according to the change frequency of the track residual fringes, a compensated baseline of the interferogram is obtained, and finally, a track error in the interferogram can be removed according to the compensated baseline of the interferogram, so that the accuracy of a deformation monitoring result of the earth surface area to be monitored is improved.
In the first aspect of the present application, as an optional implementation manner, after the computing obtains the initial baseline of the interferogram, before the removing the first geostationary phase and the second geostationary phase in the interferogram according to the initial baseline and obtaining the first differential interferogram of the interferogram, the method further includes:
calculating the first earth phase and the first terrain phase in the interferogram from the external SRTM DEM data and the initial baseline.
In this alternative embodiment, the first geostationary phase and the first topographical phase in the banded interferogram can be calculated from the external SRTM DEM data and the initial baseline.
In the first aspect of the present application, as an optional implementation manner, the calculating a change frequency of the track residual fringe includes:
and calculating the change frequency of the track residual fringe according to the fast discrete Fourier transform.
In this optional embodiment, the change frequency of the track residual fringe can be calculated by a fast discrete fourier transform algorithm.
In the first aspect of the present application, as an optional implementation manner, the correcting the initial baseline of the interferogram according to the variation frequency of the track residual fringes and obtaining the compensated baseline of the interferogram includes:
calculating to obtain a vertical baseline deviation according to the change frequency of the track residual fringe;
and compensating the initial baseline of the interferogram according to the vertical baseline deviation, and obtaining a compensated baseline of the interferogram.
In this optional embodiment, the vertical baseline deviation can be calculated by the change frequency of the track residual fringes, and then the initial baseline of the interferogram can be compensated according to the vertical baseline deviation, so as to obtain a compensated baseline of the interferogram.
In the first aspect of the present application, as an optional implementation manner, the removing the orbit error in the interferogram according to the compensated baseline of the interferogram includes:
calculating a second earth phase and a second topographic phase of the interferogram according to the compensated baseline of the interferogram and the external SRTM DEM data;
removing the second earth phase and the second topographic phase from the interferogram and obtaining a second differential interferogram;
phase unwrapping is carried out on the second differential interference pattern, and a wrapped interference pair is obtained;
fitting and calculating to obtain the orbit error in the winding interference pair and obtain the orbit error in the interference pattern;
removing the orbit error from the interferogram.
In this optional embodiment, the second terrestrial phase and the second terrestrial phase of the interferogram can be calculated according to the compensated baseline of the interferogram and the external SRTM DEM data, and then the second terrestrial phase and the second terrestrial phase can be removed from the interferogram, and a second differential interferogram is obtained, and then the second differential interferogram is subjected to phase unwrapping to obtain a wrapped interference pair, and then a track residual in the wrapped interference pair is obtained through fitting calculation, so that a track error can be obtained, and finally the track error can be removed from the interferogram.
In the first aspect of the present application, as an optional implementation manner, after the removing the second flat phase and the second ground phase from the interferogram and obtaining the second differential interferogram, the method further includes, before the phase unwrapping the second differential interferogram and obtaining a wrapped interferometric pair:
acquiring external water vapor data;
and removing the atmospheric delay phase in the second differential interference pattern according to the external water vapor data pair.
In this application embodiment, through obtaining outside steam data, and then can be according to outside steam data to the atmospheric delay phase place that gets rid of in the second difference interferogram, so, can avoid the atmospheric phase place to the estimation of ancient tunnel error.
In the first aspect of the present application, as an optional implementation manner, the calculating to obtain the initial baseline of the interferogram includes:
and calculating to obtain an initial baseline of the interferogram according to the orbit state vector of the satellite-borne radar.
In this alternative embodiment, the initial baseline of the interferogram can be calculated from the orbit state vector of the satellite-borne radar.
The second aspect of the present application discloses a satellite-borne radar data processing apparatus, the apparatus includes:
the system comprises a first acquisition module, a second acquisition module and a monitoring module, wherein the first acquisition module is used for acquiring SLC images generated by a satellite-borne radar at different time points and aiming at a to-be-monitored ground surface area, and the SLC images at every two time points form interference pair images;
the generation module is used for generating an interferogram of the earth surface area to be monitored according to the interference pair image;
the first calculation module is used for calculating to obtain an initial baseline of the interferogram;
the first removing module is used for removing a first plano phase and a first topographic phase in the interferogram according to the initial baseline and obtaining a first differential interferogram of the interferogram;
the second calculation module is used for calculating the change frequency of the track residual fringes when the concentration or the significance of the track residual fringes in the first differential interference pattern meets a preset threshold;
the correction module is used for correcting the initial baseline of the interference pattern according to the change frequency of the track residual fringe and obtaining a compensated baseline of the interference pattern;
and the second removing module is used for removing the track error in the interferogram according to the compensated base line of the interferogram.
According to the device of the second aspect of the application, by executing the satellite-borne radar data processing method, the first plano phase and the first topographic phase in the interferogram of the earth surface area to be monitored can be removed through the initial baseline, so that the first differential interferogram of the earth surface area to be monitored is obtained, further, when the density or the significance of the track residual fringes in the first differential interferogram meets a preset threshold value, the change frequency of the track residual fringes is calculated, the initial baseline of the interferogram is corrected according to the change frequency of the track residual fringes, the compensated baseline of the interferogram is obtained, finally, the track error in the interferogram can be removed according to the compensated baseline of the interferogram, and the accuracy of the deformation monitoring result of the earth surface area to be monitored is improved.
A third aspect of the present application discloses a satellite-borne radar data processing apparatus, the apparatus comprising:
a processor; and
a memory configured to store machine readable instructions, which when executed by the processor, cause the processor to perform the method of the first aspect of the present application for on-board radar data processing.
According to the device of the third aspect of the application, by executing the satellite-borne radar data processing method, the first plano phase and the first topographic phase in the interferogram of the earth surface area to be monitored can be removed through the initial baseline, so that the first differential interferogram of the earth surface area to be monitored is obtained, further, when the density or the significance of the track residual fringes in the first differential interferogram meets a preset threshold, the change frequency of the track residual fringes is calculated, further, the initial baseline of the interferogram is corrected according to the change frequency of the track residual fringes, the compensated baseline of the interferogram is obtained, finally, the track error in the interferogram can be removed according to the compensated baseline of the interferogram, and the accuracy of the deformation monitoring result of the earth surface area to be monitored is improved.
A fourth aspect of the present application discloses a storage medium storing a computer program, which is executed by a processor to perform the method for processing the satellite-borne radar data of the first aspect of the present application.
According to the storage medium of the fourth aspect of the application, by executing the satellite-borne radar data processing method, the first plano phase and the first topographic phase in the interferogram of the to-be-monitored surface area can be removed through the initial baseline, so that the first differential interferogram of the to-be-monitored surface area is obtained, further, when the density or the significance of the track residual fringe in the first differential interferogram meets a preset threshold value, the change frequency of the track residual fringe is calculated, further, the initial baseline of the interferogram is corrected according to the change frequency of the track residual fringe, the compensated baseline of the interferogram is obtained, finally, the track error in the interferogram can be removed according to the compensated baseline of the interferogram, and the accuracy of the deformation monitoring result of the to-be-monitored surface area is improved.
Drawings
In order to more clearly illustrate the technical solutions of the embodiments of the present application, the drawings that are required to be used in the embodiments of the present application will be briefly described below, it should be understood that the following drawings only illustrate some embodiments of the present application and therefore should not be considered as limiting the scope, and that those skilled in the art can also obtain other related drawings based on the drawings without inventive efforts.
Fig. 1 is a schematic flow chart illustrating a method for processing satellite-borne radar data according to an embodiment of the present application;
FIG. 2 is a schematic structural diagram of a satellite-borne radar data processing device disclosed in an embodiment of the present application;
fig. 3 is a schematic structural diagram of a satellite-borne radar data processing device disclosed in the embodiment of the present application.
Detailed Description
The technical solutions in the embodiments of the present application will be described below with reference to the drawings in the embodiments of the present application.
Referring to fig. 1, fig. 1 is a schematic flowchart illustrating a method for processing satellite-borne radar data according to an embodiment of the present disclosure. As shown in fig. 1, the method of the embodiment of the present application includes the steps of:
101. acquiring SLC images generated by the satellite-borne radar at different time points and aiming at an earth surface area to be monitored, wherein the SLC images at every two time points form interference pair images;
102. generating an interference image of the earth surface area to be monitored according to the interference pair image;
103. calculating to obtain an initial baseline of the interferogram;
104. removing the first plano-phase and the first topographic phase in the interferogram according to the initial baseline, and obtaining the interferogram;
105. when the concentration or the significance of the track residual fringes in the first differential interference pattern meets a preset threshold value, calculating the change frequency of the track residual fringes;
106. correcting an initial baseline of the interference pattern according to the change frequency of the track residual fringe and obtaining a compensated baseline of the interference pattern;
107. removing the orbit error in the interferogram according to the compensated baseline of the interferogram.
In the embodiment of the present application, for example, assuming that the Jinshajiang white lattice landslide is a ground surface area to be monitored, on the other hand, by acquiring a multi-view SLC image of the Jinshajiang white lattice landslide within a period of 7 months to 2019 months and 4 months, and then performing multi-view processing on the multi-view SLC image of the Jinshajiang white lattice landslide based on a short space-time baseline, a plurality of pairs of interference pair images can be obtained. It should be noted that, as to how to perform multi-view processing on a multi-view SLC image of the Jinsha river white lattice landslide based on a short space-time baseline, please refer to the prior art, which is not described in detail in the embodiments of the present application.
In the embodiment of the application, the interferogram of the Jinsha Jiangbangge landslide is generated by complex conjugate multiplication of two SLC images.
In the embodiment of the application, the interferogram of the Jinsha Jiangbai landslide may include an atmospheric phase, a first geostationary phase and a first topographic phase, wherein the most significant phase components in the interferogram are the first geostationary phase and the second geostationary phase respectively.
In the embodiment of the application, a first plano phase and a first topographic phase in an interferogram of an earth surface area to be monitored are removed through an initial baseline, so that a first differential interferogram of the earth surface area to be monitored can be obtained, further, when the density or the significance of track residual fringes in the first differential interferogram meets a preset threshold, the change frequency of the track residual fringes is calculated, further, the initial baseline of the interferogram is corrected according to the change frequency of the track residual fringes, a compensated baseline of the interferogram is obtained, finally, a track error in the interferogram can be removed according to the compensated baseline of the interferogram, and the accuracy of a deformation monitoring result of the earth surface area to be monitored is improved.
In this embodiment of the present application, optionally, step 102: the specific way of calculating the initial baseline of the interferogram is as follows:
and calculating to obtain an initial baseline of the interferogram according to the orbit state vector of the satellite-borne radar.
In the embodiment of the application, the orbit state vector of the satellite-borne radar comprises the instantaneous position and the speed of the satellite in a Cartesian coordinate system with the earth as the center.
In the embodiment of the present application, further optionally, the step of: calculating an initial baseline of an interferogram according to an orbit state vector of the satellite-borne radar, and comprising the following sub-steps of:
calculating according to the orbit state vector of the satellite-borne radar to obtain an orbit azimuth angle and a radar incident angle;
calculating to obtain a parallel baseline and a vertical baseline of the interferogram according to the track azimuth angle and the radar incident angle;
an initial baseline of the interferogram is determined from the parallel baseline and the perpendicular baseline.
Specifically, the parallel baseline is calculated by the following formula:
where b1 denotes a parallel base line,which is indicative of the angle of incidence of the radar,representing the track azimuth and b representing the spatial baseline.
Specifically, the vertical baseline is calculated by the following formula:
where, b2 denotes the vertical base line,which is indicative of the angle of incidence of the radar,representing the track azimuth and b representing the spatial baseline.
It can be seen that in this alternative embodiment, the initial baseline of the interferogram can be calculated from the orbit state vector of the satellite-borne radar.
In the embodiment of the present application, as an optional implementation manner, in step 102: after calculating the initial baseline of the interferogram, step 103: before removing the first geosynchronous phase and the second geosynchronous phase in the interferogram according to the initial baseline and obtaining the first differential interferogram of the interferogram, the method of the embodiment of the application further comprises the steps of:
a first plano-phase and a first topographical phase in the interferogram are calculated from the external SRTM DEM data and the initial baseline.
Specifically, the first average phase is calculated by the following formula:
wherein the content of the first and second substances,represents a first earth phase, anddenotes the wavelength, b1Parallel baselines are indicated.
Specifically, the first terrain phase is calculated by the following formula:
wherein the content of the first and second substances,a phase of a first terrain is represented and,denotes the wavelength, b2Representing a vertical baseline, R represents the distance between the satellite and the terrain, and z represents elevation.
It should be noted that the external SRTM DEM data includes the distance R, elevation z between the satellite and the ground feature.
It can thus be seen that in this alternative embodiment, the first plateau phase and the first terrain phase in the banded interferogram can be calculated from the external SRTM DEM data and the initial baseline.
In the embodiment of the present application, the differential interferogram obtained by removing the first terrestrial phase and the first terrestrial phase in the interferogram according to the initial baseline has the track residual, so that the track error needs to be further removed according to the intensity or the significance of the track residual fringes in the interferogram, wherein the track error comprises the residual terrestrial phase.
In the embodiment of the present application, the residual flat phase is calculated by the following formula:
wherein the content of the first and second substances,it indicates that the phase of the remaining flat ground,denotes the wavelength, b2A vertical baseline is shown and is,which is indicative of the angle of incidence of the radar,representing the orbital azimuth, R represents the distance between the satellite and the earth,indicating that the radar angle of incidence is spatially poor,indicating that the elevation is spatially different.
In the embodiment of the present application, the track error mainly appears as a systematically periodically varying fringe in the first differential interference pattern, wherein the fringes of the track error have different frequencies in space, i.e., the fringes of a part of the track error in the first differential interference pattern are very dense, and the whole scene image has a plurality of phase periods (2 pi); the fringe of the partial orbit error in the first differential interference image changes slowly, and the variation of the whole scene image is less than a period; still other interferograms exhibit only very insignificant systematic variation and no fringes are formed. For this situation, the embodiment of the present application removes the track error corresponding to the fringe (track residual fringe) when the intensity or significance in the first differential interference pattern satisfies the preset threshold. It should be noted that the preset threshold is that the stripes are greater than or equal to about 1 cycle.
In the embodiment of the present application, as an optional implementation manner, step 105: calculating the change frequency of the track residual stripes, comprising the following steps:
and calculating the change frequency of the track residual fringe according to the fast discrete Fourier transform.
In this optional embodiment, the change frequency of the track residual fringe can be calculated by a fast discrete fourier transform algorithm.
In the embodiment of the present application, the window size of the fast discrete fourier transform (FFT) may be 256, 512, or 2048 or 4096, wherein the specific value of the window size of the fast discrete fourier transform (FFT) may be determined according to the overall fringe variation frequency in the first differential interferogram.
In the embodiment of the application, the formulaIt can be seen that the orbit error is related to the vertical baseline, and therefore, as an alternative embodiment, step 106: correcting the initial baseline of the interferogram according to the frequency of the change of the track residual fringes and obtaining a compensated baseline of the interferogram, comprising the sub-steps of:
calculating to obtain a vertical baseline deviation according to the change frequency of the track residual fringe;
and compensating the initial baseline of the interferogram according to the vertical baseline deviation, and obtaining a compensated baseline of the interferogram.
In the embodiment of the present application, the specific calculation formula of the vertical baseline deviation is as follows:
wherein the content of the first and second substances,denotes the vertical baseline deviation, f denotes the frequency of change of the track residual fringes,which represents the wavelength of the light emitted by the light source,representing the radar angle of incidence and R representing the distance between the satellite and the ground object.
In the embodiment of the present application, the compensated baseline may be expressed as:
it can be seen that, in the optional embodiment, the vertical baseline deviation can be calculated through the change frequency of the track residual fringes, and then the initial baseline of the interferogram can be compensated according to the vertical baseline deviation, so as to obtain the compensated baseline of the interferogram.
In the embodiment of the present application, as an optional implementation manner, step 106: removing the orbit error in the interferogram from the compensated baseline of the interferogram, comprising the sub-steps of:
calculating a second earth phase and a second topographic phase of the interferogram according to the compensated baseline of the interferogram and the external SRTM DEM data;
removing the second earth phase and the second topographic phase from the interferogram and obtaining a second differential interferogram;
phase unwrapping is carried out on the second differential interference pattern, and a wrapping interference pair is obtained;
fitting calculation is carried out to obtain a track error in the winding interference pair and obtain a track error in the interference image;
the orbit error is removed from the interferogram.
In the embodiment of the present application, a specific polynomial fitting formula for obtaining the orbit error in the winding interference pair by fitting calculation is as follows:
wherein the content of the first and second substances,andcoordinates representing the azimuth direction and the distance direction,andrepresenting the weight parameter to be found.
In the embodiment of the present application, since the orbit error may affect the phase unwrapping result, the orbit error in the wrapping interference pair is obtained through the fitting calculation, and the orbit error in the interferogram can be further calculated, thereby further reducing the orbit error in the interferogram.
It should be noted that, the phase unwrapping and fitting calculation can be performed on the differential interferogram for multiple times to continuously screen the orbit error, so as to further reduce the orbit error in the interferogram.
It can be seen that, in the alternative embodiment, the second terrestrial phase and the second terrestrial phase of the interferogram can be calculated according to the compensated baseline of the interferogram and the external SRTM DEM data, so that the second terrestrial phase and the second terrestrial phase can be removed from the interferogram, the second differential interferogram is obtained, the second differential interferogram is subjected to phase unwrapping, a wrapped interference pair can be obtained, a track residual error in the wrapped interference pair is obtained through fitting calculation, a track error can be obtained, and finally the track error can be removed from the interferogram.
In this embodiment of the present application, optionally, in step 102: after generating an interferogram of the earth surface region to be monitored according to the interference pair image, step 103: before calculating the initial baseline of the interferogram, the method of the embodiment of the application further comprises the steps of:
and directly eliminating the interference pattern with the orbit error only existing in one corner of the image.
In the embodiment of the application, by eliminating the interferograms of which the orbit errors only exist in one corner of the image, the orbit errors in the interferograms can be calculated in a fitting mode through a polynomial fitting formula.
In the embodiment of the present application, as an optional implementation manner, in the step: after removing the second earth phase and the second topographic phase from the interferogram and obtaining a second differential interferogram, the steps of: before phase unwrapping the second differential interferogram and obtaining a wrapped interference pair, the method of an embodiment of the present application further includes:
acquiring external water vapor data;
and removing the atmospheric delay phase in the second differential interference pattern according to the external water vapor data pair.
In this application embodiment, through obtaining outside steam data, and then can be according to outside steam data to the atmospheric delay phase place that gets rid of in the second difference interferogram, so, can avoid the atmospheric phase place to the estimation of ancient tunnel error.
Example two
Referring to fig. 2, fig. 2 is a schematic structural diagram of a satellite-borne radar data processing apparatus according to an embodiment of the present disclosure. As shown in fig. 2, the apparatus of the embodiment of the present application includes:
the first obtaining module 201 is configured to obtain SLC images generated by the satellite-borne radar at different time points and specific to an earth surface area to be monitored, where the SLC images at every two time points form an interference pair image;
the generation module 202 is configured to generate an interferogram of a ground surface area to be monitored according to the interference pair image;
a first calculation module 203, configured to calculate an initial baseline of the interferogram;
a first removing module 204, configured to remove the first plano phase and the first topographic phase in the interferogram according to the initial baseline, and obtain a first differential interferogram of the interferogram;
a second calculating module 205, configured to calculate a change frequency of the track residual fringes when the density or the significance of the track residual fringes in the first differential interference pattern meets a preset threshold;
a correction module 206, configured to correct an initial baseline of the interferogram according to a change frequency of the track residual fringes and obtain a compensated baseline of the interferogram;
a second removal module 207 for removing the orbit error in the interferogram according to the compensated baseline of the interferogram.
According to the device, by executing the satellite-borne radar data processing method, the first plano phase and the first topographic phase in the interferogram of the earth surface area to be monitored can be removed through the initial baseline, the first differential interferogram of the earth surface area to be monitored is further obtained, further, when the density or the significance of the track residual fringes in the first differential interferogram meets a preset threshold value, the change frequency of the track residual fringes is calculated, the initial baseline of the interferogram is corrected according to the change frequency of the track residual fringes, the compensated baseline of the interferogram is obtained, finally, the track error in the interferogram can be removed according to the compensated baseline of the interferogram, and the accuracy of the deformation monitoring result of the earth surface area to be monitored is improved.
In the embodiment of the application, the interferogram of the Jinsha Jiangbangge landslide is generated by complex conjugate multiplication of two SLC images.
In the embodiment of the application, the interferogram of the Jinsha Jiangbai landslide may include an atmospheric phase, a first geostationary phase and a first topographic phase, wherein the most significant phase components in the interferogram are the first geostationary phase and the second geostationary phase respectively.
In the embodiment of the application, a first plano phase and a first topographic phase in an interferogram of an earth surface area to be monitored are removed through an initial baseline, so that a first differential interferogram of the earth surface area to be monitored can be obtained, further, when the density or the significance of track residual fringes in the first differential interferogram meets a preset threshold, the change frequency of the track residual fringes is calculated, further, the initial baseline of the interferogram is corrected according to the change frequency of the track residual fringes, a compensated baseline of the interferogram is obtained, finally, a track error in the interferogram can be removed according to the compensated baseline of the interferogram, and the accuracy of a deformation monitoring result of the earth surface area to be monitored is improved.
In this embodiment of the present application, optionally, the specific way for the first calculation module 203 to perform the calculation to obtain the initial baseline of the interferogram is as follows:
and calculating to obtain an initial baseline of the interferogram according to the orbit state vector of the satellite-borne radar.
In the embodiment of the application, the orbit state vector of the satellite-borne radar comprises the instantaneous position and the speed of the satellite in a Cartesian coordinate system with the earth as the center.
In this embodiment of the present application, further optionally, the specific way for the first calculation module 203 to calculate the initial baseline of the interferogram according to the orbit state vector of the satellite-borne radar is as follows:
calculating according to the orbit state vector of the satellite-borne radar to obtain an orbit azimuth angle and a radar incident angle;
calculating to obtain a parallel baseline and a vertical baseline of the interferogram according to the track azimuth angle and the radar incident angle;
an initial baseline of the interferogram is determined from the parallel baseline and the perpendicular baseline.
Specifically, the parallel baseline is calculated by the following formula:
where b1 denotes a parallel base line,which is indicative of the angle of incidence of the radar,representing the track azimuth and b representing the spatial baseline.
Specifically, the vertical baseline is calculated by the following formula:
where, b2 denotes the vertical base line,which is indicative of the angle of incidence of the radar,representing the track azimuth and b representing the spatial baseline.
It can be seen that in this alternative embodiment, the initial baseline of the interferogram can be calculated from the orbit state vector of the satellite-borne radar.
In this embodiment of the present application, as an optional implementation manner, the apparatus of this embodiment of the present application further includes:
a third calculation module to calculate a first plano phase and a first topographical phase in the interferogram from the external SRTM DEM data and the initial baseline.
Specifically, the first average phase is calculated by the following formula:
wherein the content of the first and second substances,represents a first earth phase, anddenotes the wavelength, b1Parallel baselines are indicated.
Specifically, the first terrain phase is calculated by the following formula:
wherein the content of the first and second substances,a phase of a first terrain is represented and,denotes the wavelength, b2Means of being perpendicularBaseline, R represents the distance between the satellite and the terrain, and z represents elevation.
It should be noted that the external SRTM DEM data includes the distance R, elevation z between the satellite and the ground feature.
It can thus be seen that in this alternative embodiment, the first plateau phase and the first terrain phase in the banded interferogram can be calculated from the external SRTM DEM data and the initial baseline.
In the embodiment of the present application, the differential interferogram obtained by removing the first terrestrial phase and the first terrestrial phase in the interferogram according to the initial baseline has the track residual, so that the track error needs to be further removed according to the intensity or the significance of the track residual fringes in the interferogram, wherein the track error comprises the residual terrestrial phase.
In the embodiment of the present application, the residual flat phase is calculated by the following formula:
wherein the content of the first and second substances,it indicates that the phase of the remaining flat ground,denotes the wavelength, b2A vertical baseline is shown and is,which is indicative of the angle of incidence of the radar,representing the orbital azimuth, R represents the distance between the satellite and the earth,indicating that the radar angle of incidence is spatially poor,indicating that the elevation is spatially different.
In the embodiment of the present application, the track error mainly appears as a systematically periodically varying fringe in the first differential interference pattern, wherein the fringes of the track error have different frequencies in space, i.e., the fringes of a part of the track error in the first differential interference pattern are very dense, and the whole scene image has a plurality of phase periods (2 pi); the fringe of the partial orbit error in the first differential interference image changes slowly, and the variation of the whole scene image is less than a period; still other interferograms exhibit only very insignificant systematic variation and no fringes are formed. For this situation, the embodiment of the present application removes the track error corresponding to the fringe (track residual fringe) when the intensity or significance in the first differential interference pattern satisfies the preset threshold. It should be noted that the preset threshold is that the stripes are greater than or equal to about 1 cycle.
In this embodiment, as an optional implementation manner, the specific way for the second calculation module 205 to calculate the change frequency of the track residual fringe is as follows:
and calculating the change frequency of the track residual fringe according to the fast discrete Fourier transform.
In this optional embodiment, the change frequency of the track residual fringe can be calculated by a fast discrete fourier transform algorithm.
In the embodiment of the present application, the window size of the fast discrete fourier transform (FFT) may be 256, 512, or 2048 or 4096, wherein the specific value of the window size of the fast discrete fourier transform (FFT) may be determined according to the overall fringe variation frequency in the first differential interferogram.
In the embodiment of the application, theFormula (II)It can be seen that the track error is related to the vertical baseline, and therefore, as an alternative embodiment, the specific way for the correction module 206 to correct the initial baseline of the interferogram according to the variation frequency of the track residual fringes and obtain the compensated baseline of the interferogram is as follows:
calculating to obtain a vertical baseline deviation according to the change frequency of the track residual fringe;
and compensating the initial baseline of the interferogram according to the vertical baseline deviation, and obtaining a compensated baseline of the interferogram.
In the embodiment of the present application, the specific calculation formula of the vertical baseline deviation is as follows:
wherein the content of the first and second substances,denotes the vertical baseline deviation, f denotes the frequency of change of the track residual fringes,which represents the wavelength of the light emitted by the light source,representing the radar angle of incidence and R representing the distance between the satellite and the ground object.
In the embodiment of the present application, the compensated baseline may be expressed as:
it can be seen that, in the optional embodiment, the vertical baseline deviation can be calculated through the change frequency of the track residual fringes, and then the initial baseline of the interferogram can be compensated according to the vertical baseline deviation, so as to obtain the compensated baseline of the interferogram.
In the embodiment of the present application, as an optional implementation manner, the specific manner in which the second removing module 207 performs removing the orbit error in the interferogram according to the compensated baseline of the interferogram is as follows:
calculating a second earth phase and a second topographic phase of the interferogram according to the compensated baseline of the interferogram and the external SRTM DEM data;
removing the second earth phase and the second topographic phase from the interferogram and obtaining a second differential interferogram;
phase unwrapping is carried out on the second differential interference pattern, and a wrapping interference pair is obtained;
fitting calculation is carried out to obtain a track error in the winding interference pair and obtain a track error in the interference image;
the orbit error is removed from the interferogram.
In the embodiment of the present application, a specific polynomial fitting formula for obtaining the orbit error in the winding interference pair by fitting calculation is as follows:
wherein the content of the first and second substances,andcoordinates representing the azimuth direction and the distance direction,andrepresenting the weight parameter to be found.
In the embodiment of the present application, since the orbit error may affect the phase unwrapping result, the orbit error in the wrapping interference pair is obtained through the fitting calculation, and the orbit error in the interferogram can be further calculated, thereby further reducing the orbit error in the interferogram.
It should be noted that, the phase unwrapping and fitting calculation can be performed on the differential interferogram for multiple times to continuously screen the orbit error, so as to further reduce the orbit error in the interferogram.
It can be seen that, in the alternative embodiment, the second terrestrial phase and the second terrestrial phase of the interferogram can be calculated according to the compensated baseline of the interferogram and the external SRTM DEM data, so that the second terrestrial phase and the second terrestrial phase can be removed from the interferogram, the second differential interferogram is obtained, the second differential interferogram is subjected to phase unwrapping, a wrapped interference pair can be obtained, a track residual error in the wrapped interference pair is obtained through fitting calculation, a track error can be obtained, and finally the track error can be removed from the interferogram.
In this embodiment of the present application, optionally, the apparatus in this embodiment of the present application further includes:
and the rejecting module is used for rejecting the interference pattern of which the track error only exists in one corner of the image.
In the embodiment of the application, by eliminating the interferograms of which the orbit errors only exist in one corner of the image, the orbit errors in the interferograms can be calculated in a fitting mode through a polynomial fitting formula.
In this embodiment of the present application, as an optional implementation manner, the apparatus of this embodiment of the present application further includes:
the second acquisition module is used for acquiring external water vapor data;
and the third removing module is used for removing the atmospheric delay phase in the second differential interference pattern according to the external water vapor data pair.
In this application embodiment, through obtaining outside steam data, and then can be according to outside steam data to the atmospheric delay phase place that gets rid of in the second difference interferogram, so, can avoid the atmospheric phase place to the estimation of ancient tunnel error.
EXAMPLE III
Referring to fig. 3, fig. 3 is a schematic structural diagram of a satellite-borne radar data processing device according to an embodiment of the present application. As shown in fig. 3, the apparatus of the embodiment of the present application includes:
a processor 301; and
the memory 302 is configured to store machine-readable instructions, and when the instructions are executed by the processor 301, the processor 301 executes the method for processing the radar data on the satellite according to the first embodiment of the present application.
The device provided by the embodiment of the application can remove the first plano phase and the first topographic phase in the interferogram of the to-be-monitored surface area through the initial baseline by executing the satellite-borne radar data processing method, further obtain the first differential interferogram of the to-be-monitored surface area, further calculate the change frequency of the track residual fringe when the density or the significance of the track residual fringe in the first differential interferogram meets a preset threshold, further correct the initial baseline of the interferogram according to the change frequency of the track residual fringe and obtain a compensated baseline of the interferogram, finally remove the track error in the interferogram according to the compensated baseline of the interferogram, and improve the accuracy of the deformation monitoring result of the to-be-monitored surface area.
Example four
The embodiment of the application discloses a storage medium, wherein a computer program is stored in the storage medium, and the computer program is executed by a processor to execute the satellite-borne radar data processing method in the embodiment of the application.
The storage medium of the embodiment of the application can remove the first plateau phase and the first terrain phase in the interferogram of the ground surface area to be monitored through the initial baseline by executing the satellite-borne radar data processing method, so as to obtain the first differential interferogram of the ground surface area to be monitored, further calculate the change frequency of the track residual fringe when the density or the significance of the track residual fringe in the first differential interferogram meets a preset threshold, further correct the initial baseline of the interferogram according to the change frequency of the track residual fringe and obtain a compensated baseline of the interferogram, and finally remove the track error in the interferogram according to the compensated baseline of the interferogram, so as to improve the accuracy of the deformation monitoring result of the ground surface area to be monitored.
In the embodiments provided in the present application, it should be understood that the disclosed apparatus and method may be implemented in other ways. The above-described embodiments of the apparatus are merely illustrative, and for example, the division of the units is only one logical division, and there may be other divisions when actually implemented, and for example, a plurality of units or components may be combined or integrated into another system, or some features may be omitted, or not executed. In addition, the shown or discussed mutual coupling or direct coupling or communication connection may be an indirect coupling or communication connection of devices or units through some communication interfaces, and may be in an electrical, mechanical or other form.
In addition, units described as separate parts may or may not be physically separate, and parts displayed as units may or may not be physical units, may be located in one place, or may be distributed on a plurality of network units. Some or all of the units can be selected according to actual needs to achieve the purpose of the solution of the embodiment.
Furthermore, the functional modules in the embodiments of the present application may be integrated together to form an independent part, or each module may exist separately, or two or more modules may be integrated to form an independent part.
It should be noted that the functions, if implemented in the form of software functional modules and sold or used as independent products, may be stored in a computer readable storage medium. Based on such understanding, the technical solution of the present application or portions thereof that substantially contribute to the prior art may be embodied in the form of a software product stored in a storage medium and including instructions for causing a computer device (which may be a personal computer, a server, or a network device) to execute all or part of the steps of the method according to the embodiments of the present application. And the aforementioned storage medium includes: various media capable of storing program codes, such as a usb disk, a removable hard disk, a Read-Only Memory (ROM), a Random Access Memory (RAM), a magnetic disk, or an optical disk.
In this document, relational terms such as first and second, and the like may be used solely to distinguish one entity or action from another entity or action without necessarily requiring or implying any actual such relationship or order between such entities or actions.
The above description is only an example of the present application and is not intended to limit the scope of the present application, and various modifications and changes may be made by those skilled in the art. Any modification, equivalent replacement, improvement and the like made within the spirit and principle of the present application shall be included in the protection scope of the present application.
Claims (10)
1. A method for processing spaceborne radar data, which is characterized by comprising the following steps:
acquiring SLC images generated by the satellite-borne radar at different time points and aiming at an earth surface area to be monitored, wherein the SLC images at every two time points form interference pair images;
generating an interference image of the earth surface area to be monitored according to the interference pair image;
calculating to obtain an initial baseline of the interferogram;
removing a first plano phase and a first topographic phase in the interferogram according to the initial baseline, and obtaining a first differential interferogram of the interferogram;
when the concentration or the significance of the track residual fringes in the first differential interference pattern meets a preset threshold value, calculating the change frequency of the track residual fringes;
correcting the initial baseline of the interference pattern according to the change frequency of the track residual fringe and obtaining a compensated baseline of the interference pattern;
and removing the track error in the interferogram according to the compensated base line of the interferogram.
2. The method of claim 1, wherein after said computing an initial baseline of said interferogram, said method further comprises, before said removing a first geostationary phase and a second geostationary phase in said interferogram from said initial baseline and obtaining a first differential interferogram of said interferogram:
calculating the first earth phase and the first terrain phase in the interferogram from the external SRTM DEM data and the initial baseline.
3. The method of claim 1, wherein said calculating a frequency of change of said track residual fringes comprises:
and calculating the change frequency of the track residual fringe according to the fast discrete Fourier transform.
4. The method of claim 3, wherein said correcting the initial baseline of the interferogram according to the frequency of change of the track residual fringes and obtaining a compensated baseline of the interferogram comprises:
calculating to obtain a vertical baseline deviation according to the change frequency of the track residual fringe;
and compensating the initial baseline of the interferogram according to the vertical baseline deviation, and obtaining a compensated baseline of the interferogram.
5. The method of claim 4, wherein removing the orbital error in the interferogram from the compensated baseline of the interferogram comprises:
calculating a second earth phase and a second topographic phase of the interferogram according to the compensated baseline of the interferogram and the external SRTM DEM data;
removing the second earth phase and the second topographic phase from the interferogram and obtaining a second differential interferogram;
phase unwrapping is carried out on the second differential interference pattern, and a wrapped interference pair is obtained;
fitting and calculating to obtain the orbit error in the winding interference pair and obtain the orbit error in the interference pattern;
removing the orbit error from the interferogram.
6. The method of claim 5, wherein after said removing the second flat earth phase and the second ground phase from the interferogram and obtaining the second differential interferogram, the method further comprises, before said phase unwrapping the second differential interferogram and obtaining a wrapped interferometric pair:
acquiring external water vapor data;
and removing the atmospheric delay phase in the second differential interference pattern according to the external water vapor data pair.
7. The method of claim 1, wherein said calculating an initial baseline for said interferogram comprises:
and calculating to obtain an initial baseline of the interferogram according to the orbit state vector of the satellite-borne radar.
8. An on-board radar data processing apparatus, the apparatus comprising:
the system comprises a first acquisition module, a second acquisition module and a monitoring module, wherein the first acquisition module is used for acquiring SLC images generated by a satellite-borne radar at different time points and aiming at a to-be-monitored ground surface area, and the SLC images at every two time points form interference pair images;
the generation module is used for generating an interferogram of the earth surface area to be monitored according to the interference pair image;
the first calculation module is used for calculating to obtain an initial baseline of the interferogram;
the first removing module is used for removing a first plano phase and a first topographic phase in the interferogram according to the initial baseline and obtaining a first differential interferogram of the interferogram;
the second calculation module is used for calculating the change frequency of the track residual fringes when the concentration or the significance of the track residual fringes in the first differential interference pattern meets a preset threshold;
the correction module is used for correcting the initial baseline of the interference pattern according to the change frequency of the track residual fringe and obtaining a compensated baseline of the interference pattern;
and the second removing module is used for removing the track error in the interferogram according to the compensated base line of the interferogram.
9. An on-board radar data processing apparatus, characterized in that the apparatus comprises:
a processor; and
a memory configured to store machine readable instructions that, when executed by the processor, cause the processor to perform the method of on-board radar data processing according to any of claims 1-7.
10. A storage medium, characterized in that the storage medium stores a computer program which is executed by a processor to perform the on-board radar data processing method according to any one of claims 1 to 7.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011121051.1A CN111965646B (en) | 2020-10-20 | 2020-10-20 | Satellite-borne radar data processing method, device, equipment and storage medium |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011121051.1A CN111965646B (en) | 2020-10-20 | 2020-10-20 | Satellite-borne radar data processing method, device, equipment and storage medium |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111965646A true CN111965646A (en) | 2020-11-20 |
CN111965646B CN111965646B (en) | 2021-01-26 |
Family
ID=73387068
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202011121051.1A Active CN111965646B (en) | 2020-10-20 | 2020-10-20 | Satellite-borne radar data processing method, device, equipment and storage medium |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111965646B (en) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112485790A (en) * | 2020-11-23 | 2021-03-12 | 湖南中大检测技术集团有限公司 | K-waveband radar-based track non-contact deformation high-precision measurement method |
CN116736306A (en) * | 2023-08-15 | 2023-09-12 | 成都理工大学 | Time sequence radar interference monitoring method based on third high-resolution |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2005085900A1 (en) * | 2004-03-01 | 2005-09-15 | Gamma Remote Sensing Research And Consulting Ag | Method for identifying radar point targets |
CN103728604A (en) * | 2013-11-19 | 2014-04-16 | 中国国土资源航空物探遥感中心 | Broadband synthetic aperture radar sub-band interferometric data processing method |
CN104111456A (en) * | 2014-07-23 | 2014-10-22 | 中国国土资源航空物探遥感中心 | Ground surface deformation high-resolution interferometric synthetic aperture radar (InSAR) monitoring method along high speed railway |
CN106526590A (en) * | 2016-11-04 | 2017-03-22 | 山东科技大学 | Method for monitoring and resolving three-dimensional ground surface deformation of industrial and mining area by means of multi-source SAR image |
CN108663017A (en) * | 2018-08-13 | 2018-10-16 | 伟志股份公司 | A kind of method of ground settlement along monitoring city underground |
-
2020
- 2020-10-20 CN CN202011121051.1A patent/CN111965646B/en active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2005085900A1 (en) * | 2004-03-01 | 2005-09-15 | Gamma Remote Sensing Research And Consulting Ag | Method for identifying radar point targets |
CN103728604A (en) * | 2013-11-19 | 2014-04-16 | 中国国土资源航空物探遥感中心 | Broadband synthetic aperture radar sub-band interferometric data processing method |
CN104111456A (en) * | 2014-07-23 | 2014-10-22 | 中国国土资源航空物探遥感中心 | Ground surface deformation high-resolution interferometric synthetic aperture radar (InSAR) monitoring method along high speed railway |
CN106526590A (en) * | 2016-11-04 | 2017-03-22 | 山东科技大学 | Method for monitoring and resolving three-dimensional ground surface deformation of industrial and mining area by means of multi-source SAR image |
CN108663017A (en) * | 2018-08-13 | 2018-10-16 | 伟志股份公司 | A kind of method of ground settlement along monitoring city underground |
Non-Patent Citations (1)
Title |
---|
殷硕文: "基于相干目标DInSAR方法的南方公路沉降监测应用研究", 《中国博士学位论文全文数据库 基础科学辑》 * |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112485790A (en) * | 2020-11-23 | 2021-03-12 | 湖南中大检测技术集团有限公司 | K-waveband radar-based track non-contact deformation high-precision measurement method |
CN112485790B (en) * | 2020-11-23 | 2023-11-24 | 中大智能科技股份有限公司 | Track non-contact deformation high-precision measurement method based on K-band radar |
CN116736306A (en) * | 2023-08-15 | 2023-09-12 | 成都理工大学 | Time sequence radar interference monitoring method based on third high-resolution |
CN116736306B (en) * | 2023-08-15 | 2023-10-24 | 成都理工大学 | Time sequence radar interference monitoring method based on third high-resolution |
Also Published As
Publication number | Publication date |
---|---|
CN111965646B (en) | 2021-01-26 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111965646B (en) | Satellite-borne radar data processing method, device, equipment and storage medium | |
Zhang et al. | A novel multitemporal InSAR model for joint estimation of deformation rates and orbital errors | |
CN107102333B (en) | Satellite-borne InSAR long and short baseline fusion unwrapping method | |
CN105487065B (en) | A kind of sequential spaceborne radar data processing method and device | |
CN104459696A (en) | SAR interference baseline precise estimating method based on flat-earth phase | |
CN108007476B (en) | Interference calibration method and system for space-based interference imaging radar altimeter | |
CN114114181B (en) | Satellite-borne SAR interference baseline correction method based on orbit error phase basis | |
Pinheiro et al. | Generation of highly accurate DEMs over flat areas by means of dual-frequency and dual-baseline airborne SAR interferometry | |
Xia et al. | Comparative analysis for combination of unwrapping and de-noising of phase data with high speckle decorrelation noise | |
Li et al. | An adaptive phase optimization algorithm for distributed scatterer phase history retrieval | |
CN103235477A (en) | Pure phase holographic projection method for inclined plane | |
Zhang et al. | Real-time computer-generated frequency-carrier Moiré profilometry with three-frequency heterodyne temporal phase unwrapping | |
Cai et al. | Gamma error correction algorithm for phase shift profilometry based on polar angle average | |
Ghulam et al. | A filtering approach to improve deformation accuracy using large baseline, low coherence DInSAR phase images | |
US20170123375A1 (en) | Apparatus and method for generating hologram | |
McMillan et al. | Tide model accuracy in the Amundsen Sea, Antarctica, from radar interferometry observations of ice shelf motion | |
CN110297242A (en) | Compressed sensing based synthetic aperture radar chromatography three-D imaging method and device | |
CN110658521A (en) | Winding phase-based GBInSAR atmospheric correction method and system | |
CN115494500A (en) | Goaf rapid detection method and system based on remote sensing interferometry and application | |
CN113341410A (en) | Large-range under-forest terrain estimation method, device, equipment and medium | |
Zheng et al. | Nonlinear error self-correction for fringe projection profilometry with minimum phase probability variance | |
CN116736306B (en) | Time sequence radar interference monitoring method based on third high-resolution | |
CN115950359B (en) | Three-dimensional reconstruction method and device and electronic equipment | |
Lee et al. | Topographic phase corrected coherence estimation using multi-pass differential SAR interferometry: differential coherence | |
Hasni et al. | Correcting ionospheric and orbital errors in spacebome SAR differential interferograms |
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 |