CN111766580A - Method for successively recursion refinement of ground position of SAR remote sensing satellite big data - Google Patents

Method for successively recursion refinement of ground position of SAR remote sensing satellite big data Download PDF

Info

Publication number
CN111766580A
CN111766580A CN202010666204.4A CN202010666204A CN111766580A CN 111766580 A CN111766580 A CN 111766580A CN 202010666204 A CN202010666204 A CN 202010666204A CN 111766580 A CN111766580 A CN 111766580A
Authority
CN
China
Prior art keywords
sar
remote sensing
data
satellite
calculation
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
Application number
CN202010666204.4A
Other languages
Chinese (zh)
Other versions
CN111766580B (en
Inventor
尤红建
董磊
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Aerospace Information Research Institute of CAS
Original Assignee
Aerospace Information Research Institute of CAS
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Aerospace Information Research Institute of CAS filed Critical Aerospace Information Research Institute of CAS
Priority to CN202010666204.4A priority Critical patent/CN111766580B/en
Publication of CN111766580A publication Critical patent/CN111766580A/en
Application granted granted Critical
Publication of CN111766580B publication Critical patent/CN111766580B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO 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/00Systems 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/88Radar or analogous systems specially adapted for specific applications
    • G01S13/89Radar or analogous systems specially adapted for specific applications for mapping or imaging
    • G01S13/90Radar or analogous systems specially adapted for specific applications for mapping or imaging using synthetic aperture techniques, e.g. synthetic aperture radar [SAR] techniques
    • G01S13/9094Theoretical aspects
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO 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/00Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/02Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
    • G01S7/41Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00 using analysis of echo signal for target characterisation; Target signature; Target cross-section
    • G01S7/418Theoretical aspects

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)
  • Image Processing (AREA)
  • Radar Systems Or Details Thereof (AREA)
  • Position Fixing By Use Of Radio Waves (AREA)

Abstract

A method for successively recursion refinement of ground position of SAR remote sensing satellite big data comprises the following steps: calculating a first synergistic factor and a ground position by utilizing data of two SAR remote sensing satellites; extracting the ground point image coordinates of newly added SAR satellite remote sensing data, and calculating the gain value of the newly added data; calculating a three-dimensional position newly added with the SAR satellite remote sensing data for updating; calculating an updated co-factor after newly adding SAR satellite remote sensing data; and repeatedly executing recursive calculation until all the remote sensing data are processed. The method carries out successive recursive refined calculation according to the sequence of the SAR satellite remote sensing data, only one newly added SAR satellite data needs to be subjected to incremental calculation during recursive calculation, and processing from the beginning is not needed.

Description

Method for successively recursion refinement of ground position of SAR remote sensing satellite big data
Technical Field
The invention relates to the technical field of remote sensing image processing, in particular to a method for successively recurrently refining ground positions of SAR remote sensing satellite big data.
Background
With the development of the synthetic aperture radar (SAR for short, hereinafter, the SAR is generally expressed by SAR) technology, the number of SAR remote sensing satellites which run in orbit at home and abroad is greatly increased, and SAR remote sensing has the characteristic of big data, so that the data volume of the obtained SAR remote sensing data is larger and larger, and the types of the SAR remote sensing data are richer. The accuracy of the three-dimensional position of the ground point can be improved by utilizing SAR remote sensing big data, but the existing SAR satellite remote sensing big data system is obviously insufficient in ground position refinement processing.
The ground position calculation of the single SAR remote sensing satellite data depends on high-precision satellite attitude and speed measurement and radar system design, and the three-dimensional position precision of a ground point calculated by three-dimensional positioning based on the two SAR remote sensing satellite data is limited. The SAR remote sensing satellite big data can be used for carrying out multi-angle, multi-orbit, multi-height and multi-resolution observation on the same ground area by a plurality of SAR satellites, and the data has internal geometric consistency, so that the precision of the ground position can be improved by fully applying repeated observation processing of the SAR remote sensing satellite big data, and the precision of the three-dimensional position can be gradually improved along with the increase of the number of the SAR remote sensing satellites, therefore, the three-dimensional position of the ground point can be refined by utilizing the SAR remote sensing satellite big data.
When the SAR remote sensing satellite big data is used for carrying out three-dimensional refinement calculation, a traditional least square integral adjustment method can be adopted, once new SAR remote sensing satellite data is added, the SAR remote sensing satellite data needs to be integrated from the beginning, namely, all the existing SAR remote sensing satellite data and the newly added SAR remote sensing satellite data need to be processed again, along with the continuous increase of the SAR remote sensing satellite big data, the complexity of calculation is rapidly increased, and the SAR remote sensing satellite big data is not suitable for real-time calculation and large-scale engineering application.
Disclosure of Invention
In view of the above, the present invention provides a method for successively recurrently refining ground position of big data of a SAR remote sensing satellite, so as to partially solve at least one of the above technical problems.
In order to achieve the above object, as an aspect of the present invention, there is provided a method for successively recurrently refining a ground position of SAR remote sensing satellite big data, comprising the steps of:
calculating a first synergistic factor and a ground position by utilizing data of two SAR remote sensing satellites;
extracting the ground point image coordinates of newly added SAR satellite remote sensing data, and calculating the gain value of the newly added data;
calculating a three-dimensional position newly added with the SAR satellite remote sensing data for updating;
calculating an updated co-factor after newly adding SAR satellite remote sensing data;
and repeatedly executing recursive calculation until all the remote sensing data are processed.
Based on the technical scheme, compared with the prior art, the method for successively recurrently refining the ground position of the SAR remote sensing satellite big data has at least one of the following beneficial effects:
(1) successive recursive refined calculation is carried out according to the sequence of the SAR satellite remote sensing data, incremental calculation is carried out only on newly added SAR satellite data during recursive calculation, and processing from the beginning is not needed.
(2) The calculation of the recursion refinement process is simple, the matrix scale is within three orders, and the real-time calculation is facilitated.
(3) The three-dimensional position of the ground point is explicitly calculated in the refinement calculation, an initial approximate value is not needed, and iteration and an approximate process are not needed in the calculation.
Drawings
FIG. 1 is a flow chart of a method for successively recurrently refining ground positions of big data of SAR remote sensing satellites according to an embodiment of the invention.
Detailed Description
In order that the objects, technical solutions and advantages of the present invention will become more apparent, the present invention will be further described in detail with reference to the accompanying drawings in conjunction with the following specific embodiments.
As shown in fig. 1, a flowchart of a method for successively recurrently refining a ground position of SAR remote sensing satellite big data according to an embodiment of the present invention specifically includes the following steps:
step A: and performing first calculation of the co-factor and the ground position by using the data of the two SAR remote sensing satellites.
This step can be further divided into the following three substeps
Step A1: extracting two-dimensional image coordinates (x) of ground points on the two SAR remote sensing satellite data by using an image matching or artificial method1,y1)(x2,y2):
Step A2: calculating a first co-factor, wherein the formula for calculating the co-factor is as follows:
Figure BDA0002578439310000031
wherein Qk-1Representing the first-time co-factor of the calculation, is a 3-row, 3-column matrix.
Ak-1A coefficient term matrix with 3 rows and 3 columns, wherein each element of the matrix is specifically expressed as follows:
Figure BDA0002578439310000032
(XS1,YS1,ZS1) And (X)S2,YS2,ZS2) The three-dimensional positions of the corresponding satellites on the two SAR remote sensing data respectively (V)X1,VY1,VZ1) And (V)X2,VY2,VZ2) The satellite three-dimensional speeds corresponding to the two SAR remote sensing data can be extracted from the auxiliary data of the SAR remote sensing satellite data and are corresponding to the row coordinate y of the ground point1,y2The values are obtained by interpolation processing.
And A isT k-1Representation pair matrix Ak-1A transposition calculation of the matrix is performed.
()-1Indicating that the matrix is subjected to an inversion calculation.
The subscript with k in the letter represents the recursion times of the refinement calculation, k is valued as 2 in the first calculation, and the value is increased by 1 every recursion time in the subsequent calculation.
Step A3: calculating the first three-dimensional position of the ground point, wherein the calculation formula is as follows:
Figure BDA0002578439310000033
wherein Xk-1For the calculated three-dimensional position of the ground, a matrix of 3 rows and 1 column is used, representing three values of the three-dimensional position.
Qk-1Representing the calculated co-factor, is the 3 row, 3 column matrix calculated by step a 2.
Lk-1A constant matrix with 3 rows and 1 column, wherein each element of the matrix is specifically expressed as follows:
Figure BDA0002578439310000041
wherein x is1And x2The coordinate values of the image column directions of the ground point on the two pieces of SAR remote sensing data are respectively obtained in step A1. (X)S1,YS1,ZS1) And (X)S1,YS1,ZS1) The three-dimensional positions of the corresponding satellites on the 1 st and 2 nd SAR remote sensing data respectively, and (V)X1,VY1,VZ1) And (V)X2,Vt2,VZ2) The satellite three-dimensional speeds corresponding to the 1 st SAR remote sensing data and the 2 nd SAR remote sensing data can be extracted from auxiliary data of the SAR remote sensing satellite data and are according to the row coordinate y corresponding to the ground point1,y2The values are obtained by interpolation processing. r is1And r2Perigee slope, M, for 1 st and 2 nd SAR dataX1And MX2Pitch resolution, λ, for 1 st and 2 nd SAR data1And λ2Radar wavelength, f, for 1 st and 2 nd SAR dataD1And fD2The four parameters are known parameters of satellite-borne SAR imaging, namely Doppler center frequencies of 1 st SAR data and 2 nd SAR data, and can be directly read from auxiliary data of SAR satellites.
And B: after newly adding an SAR satellite remote sensing data, extracting the image coordinates of the ground point, and calculating the gain value of the newly added data, wherein the calculation formula is as follows:
Figure BDA0002578439310000042
wherein G iskFor the calculated gain value of the newly added SAR remote sensing data,
Akcoefficient term matrix of 3 rows and 3 columns
Qk-1Representing the co-factor of the last recursion calculation,
e denotes an identity matrix of 3 rows and 3 columns,
()-1indicating that the matrix is subjected to an inversion calculation.
And A iskEach element of the matrix is specifically represented as follows:
Figure BDA0002578439310000043
(XS1,YS1,ZS1) And (X)Sk,YSk,ZSk) The three-dimensional positions of the corresponding satellites on the 1 st SAR remote sensing data and the newly added (k th) SAR remote sensing data respectively, and (V)X1,VY1,VZ1) And (V)Xk,VYk,VZk) The satellite three-dimensional speeds corresponding to the No. 1 SAR remote sensing data and the newly added (No. k) SAR remote sensing data can be extracted from the auxiliary data of the SAR remote sensing satellite data and are according to the row coordinate y corresponding to the ground point1,y2The values are obtained by interpolation processing.
And C: calculating the updated three-dimensional position after newly adding SAR satellite remote sensing data, wherein the calculation formula is as follows:
Xk=Xk-1+GkLk
wherein XkFor the calculated three-dimensional position of the ground, Xk-1For the three-dimensional position, X, of the ground obtained from the last recursion calculationkAnd Xk-1Are all 3 rows and 1 column matrices representing three values for a three dimensional position.
GkAnd B, calculating a gain value of newly added remote sensing data according to the step B.
LkA constant term matrix with 3 rows and 1 column, wherein each element of the matrix is specifically expressed as follows:
Figure BDA0002578439310000051
wherein x is1The coordinate value of the ground point in the image column direction of the 1 st SAR data is obtained by the step A1, xkAdding the image column direction coordinate value of the (k th) SAR data for the ground point newly, and performing the step BExtracting to obtain the extract. (X)S1,YS1,ZS1) And (X)Sk,YSk,ZSk) The three-dimensional positions of the corresponding satellites on the SAR remote sensing data of the 1 st and newly added (k) frames respectively, and (V)X1,VY1,VZ1) And (V)Xk,VYk,VZk) The satellite three-dimensional speeds corresponding to the SAR remote sensing data of the 1 st and the newly added (k) frames can be extracted from the auxiliary data of the SAR remote sensing satellite data and are according to the row coordinate y corresponding to the ground point1,ykThe values are obtained by interpolation processing. r is1And rkShort-range slope, M, for 1 st and newly added (kth) SAR dataX1And MXkPitch resolution, λ, for 1 st and newly added (kth) SAR data1And λkRadar wavelength, f, for SAR data of 1 st and newly added (k-th) framesD1And fDkThe four parameters are known parameters of satellite-borne SAR imaging and can be directly read from auxiliary data of SAR satellites for the Doppler center frequency of the 1 st SAR data and the newly added (k th) SAR data.
Step D: calculating an updated co-factor after newly adding SAR satellite remote sensing data, wherein the calculation formula is as follows:
Qk=Qk-1-GkAkQk-1
Qkrepresenting the calculated updated co-factor matrix, Qk-1Representing the co-factor matrix of the last recursion calculation.
GkAnd B, calculating a gain value of newly added remote sensing data according to the step B.
AkAnd B, a coefficient term matrix of 3 rows and 3 columns is also obtained by calculation in the step B.
Step E: the calculation of step B to step D is repeatedly performed. And for SAR satellite remote sensing big data, repeating the calculation from the step B to the step D in sequence after adding one piece of remote sensing data from the 3 rd data, and recursively finishing the refinement of the three-dimensional position of the ground point once the calculation from the step B to the step D is finished.
In conclusion, the successive recursion refinement ground position method of SAR remote sensing satellite big data of the invention directly and explicitly calculates the three-dimensional position of the ground point by applying the platform parameter and the radar parameter of the satellite-borne SAR remote sensing image without any initial value setting and iteration and approximation. And after every satellite remote sensing data is added, the current refined position of the new SAR remote sensing satellite data added can be obtained by recursion calculation only by utilizing the refined value and the co-factor of the refined value which are already obtained by the previous SAR remote sensing satellite data. It can be seen that the recursive calculation is simple, and the matrix scale is within three orders. For newly-added remote sensing data, processing from the beginning is not needed, and only the estimation result of the prior data is utilized and the contribution of the newly-added data is considered.
The above-mentioned embodiments are intended to illustrate the objects, technical solutions and advantages of the present invention in further detail, and it should be understood that the above-mentioned embodiments are only exemplary embodiments of the present invention and are not intended to limit the present invention, and any modifications, equivalents, improvements and the like made within the spirit and principle of the present invention should be included in the protection scope of the present invention.

Claims (10)

1. A method for successively recursion refinement of ground position of SAR remote sensing satellite big data is characterized by comprising the following steps:
calculating a first synergistic factor and a ground position by utilizing data of two SAR remote sensing satellites;
extracting the ground point image coordinates of newly added SAR satellite remote sensing data, and calculating the gain value of the newly added data;
calculating a three-dimensional position newly added with the SAR satellite remote sensing data for updating;
calculating an updated co-factor after newly adding SAR satellite remote sensing data;
and repeatedly executing recursive calculation until all the remote sensing data are processed.
2. The method according to claim 1, characterized in that said step of calculating the first-time co-factor and the ground position using two SAR remote sensing satellite data comprises the following sub-steps:
extracting two-dimensional image coordinates (x) of ground points on two SAR remote sensing satellite data by using image matching or artificial method1,y1)(x2,y2);
Calculating a first co-factor, wherein the formula for calculating the co-factor is as follows:
Figure FDA0002578439300000011
wherein Q isk-1The first-time co-factor representing the calculation is a 3-row and 3-column matrix; a. thek-1A coefficient term matrix of 3 rows and 3 columns; a. theT k-1Representation pair matrix Ak-1Perform the transpose calculation of the matrix, ()-1Representing the inverse calculation of the matrix; the subscript with k in the letter represents the recursion times of the refinement calculation, k is valued as 2 in the first calculation, and the value is increased by 1 every recursion time in the subsequent calculation.
Calculating the first three-dimensional position of the ground point, wherein the calculation formula is as follows:
Figure FDA0002578439300000012
wherein, Xk-1The ground three-dimensional position obtained by calculation is a matrix with 3 rows and 1 column, and represents three numerical values of the three-dimensional position; l isk-1A constant matrix of 3 rows and 1 column.
3. The method of claim 2, wherein a isk-1Each element of the matrix is specifically represented as follows:
Figure FDA0002578439300000013
(XS1,YS1,ZS1) And (X)S2,YS2,ZS2) The three-dimensional positions of the corresponding satellites on the two SAR remote sensing data respectively (V)X1,VY1,VZ1) And (V)X2,VY2,VZ2) The satellite three-dimensional speeds corresponding to the two SAR remote sensing data can be extracted from the auxiliary data of the SAR remote sensing satellite data and are respectively corresponding to the row coordinate y of the ground point1,y2The values are obtained by interpolation processing.
4. The method of claim 2, wherein L isk-1Each element of the matrix is specifically represented as follows:
Figure FDA0002578439300000021
wherein x is1And x2The coordinate values of the image array directions of the ground points on the two SAR remote sensing data are respectively r1And r2Perigee slope, M, for 1 st and 2 nd SAR dataX1And MX2For the slant range resolution, λ, of the 1 st and 2 nd SAR data1And λ2Radar wavelength, f, for 1 st and 2 nd SAR dataD1And fD2The four parameters are known parameters of satellite-borne SAR imaging, namely Doppler center frequencies of 1 st SAR data and 2 nd SAR data, and can be directly read from auxiliary data of SAR satellites.
5. The method according to claim 1, wherein the formula for extracting the ground point image coordinates of the newly added SAR satellite remote sensing data and calculating the gain value of the newly added data is as follows:
Figure FDA0002578439300000022
wherein G iskFor calculated gain values of newly added SAR remote sensing data, AkA matrix of coefficient terms of 3 rows and 3 columns, Qk-1Represents the co-factor of the last recursion calculation, E represents the 3 rows and 3 columns identity matrix, ()-1Indicating that the matrix is subjected to an inversion calculation.
6. The method of claim 5, wherein the method is performed in a batch processIs characterized in that A iskEach element of the matrix is specifically represented as follows:
Figure FDA0002578439300000023
(XS1,YS1,ZS1) And (X)Sk,YSk,ZSk) The 1 st and newly added frames are respectively the corresponding satellite three-dimensional positions on the k frame SAR remote sensing data, and (V)X1,VY1,VZ1) And (V)Xk,VYk,VZk) The three-dimensional speeds of the satellites corresponding to the No. 1 SAR remote sensing data and the new SAR remote sensing data, namely the three-dimensional speeds of the satellites corresponding to the kth SAR remote sensing data, can be extracted from auxiliary data of the SAR remote sensing satellite data and are corresponding to the row-direction coordinate y of the ground point1,y2The values are obtained by interpolation processing.
7. The method according to claim 1, wherein in the step of calculating the updated three-dimensional position of the newly added SAR satellite remote sensing data, the calculation formula of the updated three-dimensional position is as follows:
Xk=Xk-1+GkLk
wherein, XkFor the calculated three-dimensional position of the ground, Xk-1For the three-dimensional position, X, of the ground obtained from the last recursion calculationkAnd Xk-1All are 3 rows and 1 column matrices, representing three values of three-dimensional position; l iskA matrix of constant terms in 3 rows and 1 column.
8. The method of claim 7, wherein L iskEach element of the matrix is specifically represented as follows:
Figure FDA0002578439300000031
wherein r is1And rkShort-range slope, M, for the 1 st and newly added frames, i.e., the kth SAR dataX1And MXkFor the 1 st and newly added frames, i.e. the slope resolution, lambda, of the k-th SAR data1And λkFor the 1 st and newly added frames, i.e. the radar wavelength, f, of the k-th SAR dataD1And fDkThe four parameters are known parameters of satellite-borne SAR imaging and can be directly read from auxiliary data of SAR satellites for the 1 st frame and a newly added frame, namely the k frame SAR data Doppler center frequency.
9. The method according to claim 1, wherein in the step of calculating the updated co-factor after adding new SAR satellite remote sensing data, the updated co-factor calculation formula is as follows:
Qk=Qk-1-GkAkQk-1
wherein Q iskRepresenting the calculated updated co-factor matrix, Qk-1Representing the co-factor matrix of the last recursion calculation.
10. The method according to claim 1, wherein in the step of repeatedly performing recursive calculation until all the remote sensing data are processed, the recursive calculation comprises extracting ground point image coordinates of newly added SAR satellite remote sensing data and calculating gain values of the newly added SAR satellite remote sensing data; calculating a three-dimensional position newly added with the SAR satellite remote sensing data for updating; and calculating an updated co-factor after newly adding the SAR satellite remote sensing data.
CN202010666204.4A 2020-07-10 2020-07-10 Method for successively recursion refinement of ground position of SAR remote sensing satellite big data Active CN111766580B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010666204.4A CN111766580B (en) 2020-07-10 2020-07-10 Method for successively recursion refinement of ground position of SAR remote sensing satellite big data

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010666204.4A CN111766580B (en) 2020-07-10 2020-07-10 Method for successively recursion refinement of ground position of SAR remote sensing satellite big data

Publications (2)

Publication Number Publication Date
CN111766580A true CN111766580A (en) 2020-10-13
CN111766580B CN111766580B (en) 2022-03-11

Family

ID=72726703

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010666204.4A Active CN111766580B (en) 2020-07-10 2020-07-10 Method for successively recursion refinement of ground position of SAR remote sensing satellite big data

Country Status (1)

Country Link
CN (1) CN111766580B (en)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112050735A (en) * 2020-09-09 2020-12-08 中国科学院空天信息创新研究院 Method for recursion refinement of ground position of optical remote sensing satellite big data and storage medium
CN115168788A (en) * 2022-09-07 2022-10-11 中国科学院空天信息创新研究院 Method, device, equipment and medium for determining satellite remote sensing big data

Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6911932B1 (en) * 2003-10-28 2005-06-28 Sandia Corporation Two antenna, two pass interferometric synthetic aperture radar
WO2007038068A2 (en) * 2005-09-26 2007-04-05 Honeywell International Inc. Methods and systems for measuring terrain height
CN106443676A (en) * 2016-09-29 2017-02-22 中国科学院电子学研究所 Scarce control point space-borne synthetic aperture radar image ground positioning method
CN107728145A (en) * 2017-10-24 2018-02-23 中国科学院电子学研究所 The method that ground point three-dimensional position is calculated based on sequence satellite-borne SAR image
EP3364212A1 (en) * 2017-02-15 2018-08-22 Deutsches Zentrum für Luft- und Raumfahrt e.V. A method and an apparatus for computer-assisted processing of sar raw data
CN108594224A (en) * 2018-03-30 2018-09-28 中国电力工程顾问集团中南电力设计院有限公司 Merge the three-D sequential deformation monitoring method of different platform and track SAR data
CN110412574A (en) * 2019-09-05 2019-11-05 河海大学 A kind of distributed object InSAR timing sequence process method and apparatus of temporal and spatial coherence enhancing
US20190339381A1 (en) * 2018-05-03 2019-11-07 University Of Electronic Science And Technology Of China Method, device and storage medium for extracting height and deformation information of high voltage transmission tower by using sar tomography

Patent Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6911932B1 (en) * 2003-10-28 2005-06-28 Sandia Corporation Two antenna, two pass interferometric synthetic aperture radar
WO2007038068A2 (en) * 2005-09-26 2007-04-05 Honeywell International Inc. Methods and systems for measuring terrain height
CN106443676A (en) * 2016-09-29 2017-02-22 中国科学院电子学研究所 Scarce control point space-borne synthetic aperture radar image ground positioning method
EP3364212A1 (en) * 2017-02-15 2018-08-22 Deutsches Zentrum für Luft- und Raumfahrt e.V. A method and an apparatus for computer-assisted processing of sar raw data
CN107728145A (en) * 2017-10-24 2018-02-23 中国科学院电子学研究所 The method that ground point three-dimensional position is calculated based on sequence satellite-borne SAR image
CN108594224A (en) * 2018-03-30 2018-09-28 中国电力工程顾问集团中南电力设计院有限公司 Merge the three-D sequential deformation monitoring method of different platform and track SAR data
US20190339381A1 (en) * 2018-05-03 2019-11-07 University Of Electronic Science And Technology Of China Method, device and storage medium for extracting height and deformation information of high voltage transmission tower by using sar tomography
CN110412574A (en) * 2019-09-05 2019-11-05 河海大学 A kind of distributed object InSAR timing sequence process method and apparatus of temporal and spatial coherence enhancing

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
NIANGANG JIAO: ""A New RD-RFM Stereo Geolocation Model for 3D Geo-Information Reconstruction of SAR-Optical Satellite Image Pairs"", 《IEEE ACCESS》 *
NIANGANG JIAO: ""Geolocation Accuracy Improvement of Multiobserved GF-3 Spaceborne SAR Imagery"", 《IEEE GEOSCIENCE AND REMOTE SENSING LETTERS 》 *
马艳鸽 等: ""函数模型与随机模型双约束的GPS-InSAR数据融合方法建立三维形变场"", 《大地测量与地球动力学》 *

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112050735A (en) * 2020-09-09 2020-12-08 中国科学院空天信息创新研究院 Method for recursion refinement of ground position of optical remote sensing satellite big data and storage medium
CN112050735B (en) * 2020-09-09 2022-02-22 中国科学院空天信息创新研究院 Method for recursion refinement of ground position of optical remote sensing satellite big data and storage medium
CN115168788A (en) * 2022-09-07 2022-10-11 中国科学院空天信息创新研究院 Method, device, equipment and medium for determining satellite remote sensing big data
CN115168788B (en) * 2022-09-07 2022-11-22 中国科学院空天信息创新研究院 Method, device, equipment and medium for determining satellite remote sensing big data

Also Published As

Publication number Publication date
CN111766580B (en) 2022-03-11

Similar Documents

Publication Publication Date Title
CN111766580B (en) Method for successively recursion refinement of ground position of SAR remote sensing satellite big data
CN111337893B (en) Off-grid DOA estimation method based on real-value sparse Bayesian learning
CN109284677B (en) Bayesian filtering target tracking algorithm
CN110728706B (en) SAR image fine registration method based on deep learning
CN114862932B (en) BIM global positioning-based pose correction method and motion distortion correction method
CN111768337B (en) Image processing method and device and electronic equipment
CN109298420B (en) Moving target iteration minimum entropy imaging method and device of synthetic aperture radar
CN110632559B (en) MIMO pulse radar waveform design and optimization method based on application characteristics
CN116500575B (en) Extended target tracking method and device based on variable decibel leaf theory
CN108562900B (en) SAR image geometric registration method based on elevation correction
CN114972078B (en) Method and system for improving uncontrolled geometric quality of domestic optical satellite image by SAR image
CN110954860A (en) DOA and polarization parameter estimation method
CN109696651B (en) M estimation-based direction-of-arrival estimation method under low snapshot number
CN117310682A (en) SAR equivalent radar speed estimation method based on dichotomy search
CN110361733B (en) Medium orbit SAR (synthetic aperture radar) large squint imaging method based on time-frequency joint resampling
CN110912535B (en) Novel non-pilot Kalman filtering method
CN116563354A (en) Laser point cloud registration method combining feature extraction and clustering algorithm
CN114357719B (en) Method and system for accurately correcting soil erosion theoretical value by using Kalman filtering
CN109459768A (en) A kind of quick satellite selection method based on Big Dipper satellite signal intensity weight optimization model
CN115221768A (en) Improved unscented particle filter direct tracking method based on multi-nested array
CN115453523A (en) Scanning radar sparse target batch processing super-resolution method
CN110261839B (en) Enhanced spectrum diversity azimuth offset estimation method based on double samples
CN114417917A (en) Direct positioning method under unknown mutual coupling condition
CN112050735B (en) Method for recursion refinement of ground position of optical remote sensing satellite big data and storage medium
CN115932907B (en) Single-target direct tracking method based on Kalman filter and particle swarm optimization

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