CN113777606B - 分布式geo sar三维形变反演多角度选取方法及装置 - Google Patents
分布式geo sar三维形变反演多角度选取方法及装置 Download PDFInfo
- Publication number
- CN113777606B CN113777606B CN202110921920.7A CN202110921920A CN113777606B CN 113777606 B CN113777606 B CN 113777606B CN 202110921920 A CN202110921920 A CN 202110921920A CN 113777606 B CN113777606 B CN 113777606B
- Authority
- CN
- China
- Prior art keywords
- observation
- dimensional deformation
- coefficient
- angle
- distributed
- 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.)
- Active
Links
- 238000010187 selection method Methods 0.000 title claims abstract description 19
- 239000011159 matrix material Substances 0.000 claims abstract description 55
- 238000005259 measurement Methods 0.000 claims abstract description 38
- 238000000034 method Methods 0.000 claims abstract description 25
- 239000013598 vector Substances 0.000 claims description 21
- 238000005305 interferometry Methods 0.000 claims description 16
- 238000005286 illumination Methods 0.000 claims description 2
- 230000000875 corresponding effect Effects 0.000 description 14
- 238000010586 diagram Methods 0.000 description 7
- 238000004088 simulation Methods 0.000 description 5
- 238000004458 analytical method 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
- 239000000243 solution Substances 0.000 description 3
- 238000012795 verification Methods 0.000 description 3
- 238000004891 communication Methods 0.000 description 2
- 238000011156 evaluation Methods 0.000 description 2
- 238000003384 imaging method Methods 0.000 description 2
- 230000001678 irradiating effect Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000012545 processing Methods 0.000 description 2
- 238000010845 search algorithm Methods 0.000 description 2
- 230000001174 ascending effect Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 230000002596 correlated effect Effects 0.000 description 1
- 230000007423 decrease Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 239000012895 dilution Substances 0.000 description 1
- 238000010790 dilution Methods 0.000 description 1
- 230000006870 function Effects 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 230000002265 prevention Effects 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 238000007670 refining Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 230000035945 sensitivity Effects 0.000 description 1
- 239000000126 substance Substances 0.000 description 1
- 238000012876 topography Methods 0.000 description 1
- 230000000007 visual effect Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
- G01S13/88—Radar or analogous systems specially adapted for specific applications
- G01S13/89—Radar or analogous systems specially adapted for specific applications for mapping or imaging
- G01S13/90—Radar or analogous systems specially adapted for specific applications for mapping or imaging using synthetic aperture techniques, e.g. synthetic aperture radar [SAR] techniques
- G01S13/9021—SAR image post-processing techniques
- G01S13/9023—SAR image post-processing techniques combined with interferometric techniques
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
- G01S13/88—Radar or analogous systems specially adapted for specific applications
- G01S13/89—Radar or analogous systems specially adapted for specific applications for mapping or imaging
- G01S13/90—Radar or analogous systems specially adapted for specific applications for mapping or imaging using synthetic aperture techniques, e.g. synthetic aperture radar [SAR] techniques
- G01S13/904—SAR modes
- G01S13/9058—Bistatic or multistatic SAR
-
- 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/9094—Theoretical aspects
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
本发明提供一种分布式GEO SAR三维形变反演多角度选取方法及装置,所述方法包括获取分布式GEO SAR***对观测目标全部照射的空间范围内的速度、坐标、场景坐标,得到三维形变模型系数矩阵;通过各观测角度下的干涉图相位方差获得相位误差的协方差矩阵;根据所述三维形变模型系数矩阵和所述相位误差的协方差矩阵,计算所述分布式GEO SAR***对观测目标在全部照射时间段内的测量精度系数最小定位精度系数;基于所述分布式GEO SAR***对观测目标在全部照射时间段内的测量精度系数最小定位精度系数,进行最优组合选择。根据本发明的方案,可有效提高地表三维形变反演精度。
Description
技术领域
本发明涉及雷达领域,尤其涉及一种分布式GEO SAR三维形变反演多角度选取方法及装置。
背景技术
地球同步轨道合成孔径雷达(Geosynchronous Synthetic Aperture Radar,GEOSAR)运行在约36000km轨道高度上,具有快速重访能力和连续覆盖性能。但是,对于中等倾角的星下点为“小8”轨迹的GEO SAR而言,其在实现形变反演应用中面临目标散射时间去相关、大气扰动去相关,且其东西方向角度差异小,三维形变反演性能较差。在单星GEO SAR基础上增加被动接收的从星构成的分布式GEO SAR,可克服上述不足,对观测地区实现厘米级至毫米级的应急响应三维形变测量,在防灾减灾和灾害监测方面具有重要意义。
然而,与单基地SAR不同,分布式GEO SAR在三维形变反演方面面临以下问题:(1)收发平台分置,单基地SAR的三维形变反演模型及精度分析方法失效;(2)***观测角度多,无法选择最优的观测角度进行三维形变反演。因此,需建立分布式GEO SAR三维形变反演精度分析方法,并选择反演性能最佳的观测角度组合。但现有的形变反演方法并不能解决上述问题。
发明内容
为解决上述技术问题,本发明提出了一种分布式GEO SAR三维形变反演多角度选取方法及装置,所述方法及装置,用以解决现有技术中分布式GEO SAR在三维形变反演方面面临的技术问题。
根据本发明的第一方面,提供一种分布式GEO SAR三维形变反演多角度选取方法,所述方法包括以下步骤:
步骤S101:获取分布式GEO SAR***对观测目标全部照射的空间范围内的速度、坐标、场景坐标;其中场景中心记为P0,根据场景中心P0计算正北、正东和垂直三个方向的单位矢量选取任意三个观测角度,获取各观测角度对应的轨道孔径中心时刻的差分干涉测量的等效孔径中心位置SAn,计算轨道孔径中心时刻的双基地角βn和GEO卫星到目标斜距的单位矢量/>得到三维形变模型系数矩阵Θ,其中n=1,2,3;
步骤S102:由干涉图的多视视数N和相关系数γ计算各观测角度下的干涉图相位方差,通过各观测角度下的干涉图相位方差获得相位误差的协方差矩阵CΦ;
步骤S103:根据所述三维形变模型系数矩阵Θ和所述相位误差的协方差矩阵CΦ,计算所述分布式GEO SAR***对观测目标在全部照射时间段内的测量精度系数最小定位精度系数;
步骤S104:基于所述分布式GEO SAR***对观测目标在全部照射时间段内的测量精度系数最小定位精度系数,进行最优组合选择:
在所述分布式GEO SAR***对所述观测目标的全部照射轨道范围Sa上,获得最优的三个观测角度组合,其中,SA1,SA2,SA3分别为各观测角度对应的轨道孔径中心时刻的差分干涉测量的等效孔径中心位置,P0为观测场景中心,PDOPd(SA1,SA2,SA3,P0)为分布式SAR***在SA1,SA2,SA3位置对P0进行三维形变反演的测量精度系数。
根据本发明第二方面,提供一种分布式GEO SAR三维形变反演多角度选取装置,所述装置包括:
系数矩阵获取模块:配置为获取分布式GEO SAR***对观测目标全部照射的空间范围内的速度、坐标、场景坐标;其中场景中心记为P0,根据场景中心P0计算正北、正东和垂直三个方向的单位矢量选取任意三个观测角度,获取各观测角度对应的轨道孔径中心时刻的差分干涉测量的等效孔径中心位置SAn,计算轨道孔径中心时刻的双基地角βn和GEO卫星到目标斜距的单位矢量/>得到三维形变模型系数矩阵Θ,其中n=1,2,3;
协方差矩阵获取模块:配置为由干涉图的多视视数N和相关系数γ计算各观测角度下的干涉图相位方差,通过各观测角度下的干涉图相位方差获得相位误差的协方差矩阵CΦ;
精度系数计算模块:配置为根据所述三维形变模型系数矩阵Θ和所述相位误差的协方差矩阵CΦ,计算所述分布式GEO SAR***对观测目标在全部照射时间段内的测量精度系数最小定位精度系数;
确定观测角度模块:配置为基于所述分布式GEO SAR***对观测目标在全部照射时间段内的测量精度系数最小定位精度系数,进行最优组合选择:
在所述分布式GEO SAR***对所述观测目标的全部照射轨道范围Sa上,获得最优的三个观测角度组合,其中,SA1,SA2,SA3分别为各观测角度对应的轨道孔径中心时刻的差分干涉测量的等效孔径中心位置,P0为观测场景中心,PDOPd(SA1,SA2,SA3,P0)为分布式SAR***在SA1,SA2,SA3位置对P0进行三维形变反演的测量精度系数。
根据本发明第三方面,提供一种分布式GEO SAR三维形变反演多角度选取***,包括:
处理器,用于执行多条指令;
存储器,用于存储多条指令;
其中,所述多条指令,用于由所述存储器存储,并由所述处理器加载并执行如前所述的分布式GEO SAR三维形变反演多角度选取方法。
根据本发明第四方面,提供一种计算机可读存储介质,所述存储介质中存储有多条指令;所述多条指令,用于由处理器加载并执行如前所述的分布式GEO SAR三维形变反演多角度选取方法。
根据本发明的上述方案,本发明的方法,建立了分布式GEO SAR三维形变反演模型及其精度分析方法;提供了如何选择最优的观测角度进行三维形变反演的方法和在观测角度选取的过程中的三维形变反演精度评估的准则,实现对三维形变反演精度综合直观的评估。通过选择最优的角度组合进行形变反演,可有效提高地表三维形变反演精度。
上述说明仅是本发明技术方案的概述,为了能够更清楚了解本发明的技术手段,并可依照说明书的内容予以实施,以下以本发明的较佳实施例并配合附图详细说明如后。
附图说明
构成本发明的一部分的附图用来提供对本发明的进一步理解,本发明提供如下附图进行说明。在附图中:
图1为本发明一个实施方式的分布式GEO SAR三维形变反演多角度选取方法流程示意图;
图2为本发明一个实施方式的三维形变反演多角度选取方法示意图;
图3为本发明一个实施方式的SAR数据获取原理示意图;
图4(A)为本发明一个实施方式的仿真场景高程图;
图4(B)为本发明一个实施方式的高度方向形变量示意图;
图5(A)为本发明一个实施方式的南北方向的形变反演仿真验证结果示意图;
图5(B)为本发明一个实施方式的东西方向的形变反演仿真验证结果示意图;
图5(C)为本发明一个实施方式的高度方向的形变反演仿真验证结果示意图;
图6为本发明一个实施方式的分布式GEO SAR三维形变反演多角度选取装置结构框图。
具体实施方式
首先结合图1说明本发明一个实施方式的分布式GEO SAR三维形变反演多角度选取方法流程。如图1-2所示,所述方法包括以下步骤:
步骤S101:获取分布式GEO SAR***对观测目标全部照射的空间范围内的速度、坐标、场景坐标;其中场景中心记为P0,根据场景中心P0计算正北、正东和垂直三个方向的单位矢量选取任意三个观测角度,获取各观测角度对应的轨道孔径中心时刻的差分干涉测量的等效孔径中心位置SAn,计算轨道孔径中心时刻的双基地角βn和GEO卫星到目标斜距的单位矢量/>得到三维形变模型系数矩阵Θ,其中n=1,2,3;
步骤S102:由干涉图的多视视数N和相关系数γ计算各观测角度下的干涉图相位方差,通过各观测角度下的干涉图相位方差获得相位误差的协方差矩阵CΦ;
步骤S103:根据所述三维形变模型系数矩阵Θ和所述相位误差的协方差矩阵CΦ,计算所述分布式GEO SAR***对观测目标在全部照射时间段内的测量精度系数最小定位精度系数;
步骤S104:基于所述分布式GEO SAR***对观测目标在全部照射时间段内的测量精度系数最小定位精度系数,进行最优组合选择:
在所述分布式GEO SAR***对所述观测目标的全部照射轨道范围Sa上,获得最优的三个观测角度组合,其中,SA1,SA2,SA3分别为各观测角度对应的轨道孔径中心时刻的差分干涉测量的等效孔径中心位置,P0为观测场景中心,PDOPd(SA1,SA2,SA3,P0)为分布式SAR***在SA1,SA2,SA3位置对P0进行三维形变反演的测量精度系数。
如图3所示,所述步骤S101:获取分布式GEO SAR***对观测目标全部照射的空间范围内的速度、坐标、场景坐标;其中场景中心记为P0,根据场景中心P0计算正北、正东和垂直三个方向的单位矢量选取任意三个观测角度,获取各观测角度对应的轨道孔径中心时刻的差分干涉测量的等效孔径中心位置SAn,计算轨道孔径中心时刻的双基地角βn和GEO卫星到目标斜距的单位矢量/>得到三维形变模型系数矩阵Θ,其中n=1,2,3,其中:
所述三维形变模型系数矩阵Θ,计算方式为:
其中λ为雷达波长,<·,·>表示向量内积。
具体的三维形变反演模型建立原理如下所示。
假设SA1、SA2和SA3分别代表分布式GEO SAR在三个不同观测角度下获得的等效孔径中心位置,对应测得的形变干涉相位分别为Φ1、Φ2和Φ3,则建立的基于多角度观测的三维形变反演模型表示为
Φ=Θd (3)
其中,Φ=(Φ1,Φ2,Φ3)T为相位测量矢量,Θ为所述三维形变模型系数矩阵,d=(dN,dE,dU)T表示直角坐标系下的地表三维形变信息,dN为南北方向形变量,dE为东西方向形变量,dU为垂直方向形变量,T为转置。
D-InSAR(Differential Synthetic Aperture Radar Interferometry,D-InSAR)测量的形变场并不是地表真实的垂直向、东西向、南北向形变,而是地表形变场在卫星视线向上的投影,因此需要确定分布式GEO SAR在不同观测视角对目标场景区域的形变测量敏感方向和观测形变量,这一点是现有的形变反演研究所没有考虑的。
定义ST和SR分别为发射机和接收机位置,SA为等效的孔径中心位置。P为场景观测目标的位置,下标1和2分别代表前后两次测量。干涉相位的表示为
其中
对式(5)进行Taylor一阶展开得到
将式(6)代入式(4)中,可分离出形变相位表示为
形变相位对形变求导得到
其中β为和/>的夹角,即双基地角。
由式(8)可知对形变最敏感的方向也就是测量方向
结合公式(3)和(8),即可得到三维形变模型系数矩阵Θ表达式(2)。
所述步骤S102:由干涉图的多视视数N和相关系数γ计算各观测角度下的干涉图相位方差,通过各观测角度下的干涉图相位方差获得相位误差的协方差矩阵CΦ,其中:
所述相关系数γ可以基于形变测量的时间和RFI去相关等影响进行假设。
一般地,GEO SAR***对目标的覆盖角度范围很大,不同观测角度的GEO SAR成像子孔径之间没有重叠。因此,这些子孔径数据可以被认为是独立的,不同观测角度下获得的干涉图之间也是不相关的。这时,所述协方差矩阵CΦ中非对角线上的元素均为0,CΦ只由对角线上的元素确定表示为
其中,各观测角度下的干涉图相位方差由对应干涉图的多视视数N和相关系数γ计算,
所述步骤S103:根据所述三维形变模型系数矩阵Θ和所述相位误差的协方差矩阵CΦ,计算所述分布式GEO SAR***对观测目标在全部照射时间段内的测量精度系数最小定位精度系数PDOPd,其中:
本实施例中,区分与不区分各个观测角度下生成干涉图的相位误差大小的情况,计算所述分布式GEO SAR***对观测目标在全部照射时间段内的测量精度系数最小定位精度系数PDOPd的方式不同,
区分各个观测角度下生成干涉图的相位误差大小,
不区分各个观测角度下生成干涉图的相位误差大小,
其中,tr(·)为矩阵的迹,
本实施例中,引入最小定位精度系数(Position Dilution of Precision,PDOP)的概念来评估不同子孔径数据组合下场景目标三维形变测量精度的综合表现。根据经典的PDOP定义,定义三维形变的测量精度系数PDOPd表示为
其中,Λ2表示三个维度形变方差和,表征了三维形变反演的综合精度;表示三张干涉图的相位方差和。
对式(3)使用最小二乘方法得到三维形变量的估计表示为
其中,CΦ为相位误差的协方差矩阵。
根据式(14),估计获得的形变量的协方差矩阵表示为
利用多角度观测方法获得三维形变量的方差矢量表示为
将式(9)和式(16)代入式(13)得到PDOPd表达式为
若不区分各个子孔径下生成干涉图的相位误差大小,则上式可以简化为
所述步骤S104:基于所述分布式GEO SAR***对观测目标在全部照射时间段内的测量精度系数最小定位精度系数,进行最优组合选择:
在所述分布式GEO SAR***对所述观测目标的全部照射轨道范围Sa上,获得最优的三个观测角度组合,其中,SA1,SA2,SA3分别为各观测角度对应的轨道孔径中心时刻的差分干涉测量的等效孔径中心位置,P0为观测场景中心,PDOPd(SA1,SA2,SA3,P0)为分布式SAR***在SA1,SA2,SA3位置对P0进行三维形变反演的测量精度系数,其中:
分布式GEO SAR***对观测目标在全部照射时间段内选取的子孔径生成的干涉图的相位噪声是稳定的,三维形变反演的精度主要取决于用于三维形变反演的子孔径的空间几何关系,即观测角度位置的选择。因此,需要在分布式GEO SAR***对观测目标全部照射的空间范围内寻找最优的观测角度位置数据进行三维形变反演,使反演的三维形变综合精度达到最佳。考虑基于所述分布式GEO SAR***对观测目标在全部照射时间段内的测量精度系数最小定位精度系数,进行最优组合选择,
由于GEO SAR的轨迹非常弯曲和复杂,很难直接获得式(1)的解析解。因此,可以在分布式GEO SAR***对目标的全部照射轨道范围Sa上,设置步进量或根据应用的精度需求进一步插值细化步进间隔,利用搜索算法获得最优的三个观测角度组合。
以下结合具体参数说明一种分布式GEO SAR三维形变反演多角度选取方法。
在本实例中,我们考虑的分布式GEO SAR***,由一颗具备收发能力的主星,和一颗仅被动接收信号从星伴随飞行组成。卫星轨道根数如表1所示,平台参数如表2所示。实现SAR三维形变反演至少需要三对SAR图像对,其中两对图像对由主卫星多次照射同一个地区获取,一对图像对由主星发射从星接收照射同一个地区获取,SAR数据获取原理如图3所示。
SAR成像的场景中心位于东经104.4°、北纬36.9°。场景中间放置了一个底座为1km×1km、高度为1m的金字塔。塔顶具有5cm的高度向形变量,从金字塔顶点到塔底存在沿着金字塔地形变化梯度方向线性下降的高度向形变量。东西向和南北向无形变。假设形变测量的时间基线为1天,卫星对该场景成像,每幅图像大小为120像素×120像素,像素间隔为10m。仿真场景的三维图像和高度形变量见图4(A)-4(B)。
表1
轨道根数 | 主星 | 从星 |
半长轴(km) | 42164 | 42164 |
离心率 | 0 | 0 |
倾角(°) | 16 | 16 |
近地点幅角(°) | 0 | 0 |
升交点赤经经度(°) | 88 | 127.8 |
平近点角(°) | 0 | 0 |
表2
平台参数 | 数值 |
波长(m) | 0.24 |
带宽(MHz) | 60 |
脉宽(μs) | 20 |
脉冲重复频率(Hz) | 50 |
合成孔径时间(s) | 1000 |
首先按照步骤1,获取分布式GEOSAR***对观测目标全部照射的空间范围内的速度和坐标、场景坐标等信息(这里通过软件Satellite Tool Kit,设置卫星和目标位置后获取)。根据场景中心P0计算正北、正东和垂直三个方向的单位矢量计算对应每个观测角度的轨道孔径中心时刻的差分干涉测量的等效孔径中心位置SAn,n=1,2,3,计算孔径中心时刻的双基地角βn,n=1,2,3和GEO卫星到目标斜距的单位矢量/>n=1,2,3,得到,得到三维形变模型系数矩阵Θ。
步骤2,由等效多视视数N和干涉图的相关系数γ各观测角度下的干涉图相位方差n=1,2,3,获得相位误差的协方差矩阵CΦ,此处N=1,相关系数γ为0.8。
步骤3,根据所得到的三维形变模型系数矩阵Θ和相位误差的协方差矩阵CΦ,不区分各个子孔径下生成干涉图的相位误差大小,使用式(12)计算测量精度系数PDOPd值。
步骤4,利用PDOP准则下的最优角度组合选择方法,
在分布式GEO SAR***对目标全部照射的空间范围Sa上,设置步进量为600s,利用搜索算法获得最优的三个观测角度组合,孔径中心位置对应的真近点角分别为9.9°,89.4°和124.1°。
依据选取最佳的三个观测角度,然后在这三个角度对应的轨道位置对场景进行差分干涉处理,图5(A)-5(C)展示了最后的三维形变反演结果。可以看出高度向具有明显的形变,这与设置的高度向形变是一致的。
同时,比较使用选取方法与不使用选取方法,即任意选取三个观测角度组合,孔径中心位置对应的真近点角分别为39.7°,121.7°和86.9°的两组形变反演精度对比结果如表3,证明了本方法的有效性。
表3三维形变反演精度评估
使用选取方法 | 不使用选取方法 | |
南北方向(cm) | 0.7 | 6.8 |
东西方向(cm) | 1.5 | 1.0 |
高度方向(cm) | 3.9 | 2.5 |
PDOPd | 6.2 | 21.6 |
本发明实施例进一步给出一种分布式GEO SAR三维形变反演多角度选取装置,如图6所示,所述装置包括:
系数矩阵获取模块:配置为获取分布式GEO SAR***对观测目标全部照射的空间范围内的速度、坐标、场景坐标;其中场景中心记为P0,根据场景中心P0计算正北、正东和垂直三个方向的单位矢量选取任意三个观测角度,获取各观测角度对应的轨道孔径中心时刻的差分干涉测量的等效孔径中心位置SAn,计算轨道孔径中心时刻的双基地角βn和GEO卫星到目标斜距的单位矢量/>得到三维形变模型系数矩阵Θ,其中n=1,2,3;
协方差矩阵获取模块:配置为由干涉图的多视视数N和相关系数γ计算各观测角度下的干涉图相位方差,通过各观测角度下的干涉图相位方差获得相位误差的协方差矩阵CΦ;
精度系数计算模块:配置为根据所述三维形变模型系数矩阵Θ和所述相位误差的协方差矩阵CΦ,计算所述分布式GEO SAR***对观测目标在全部照射时间段内的测量精度系数最小定位精度系数;
确定观测角度模块:配置为基于所述分布式GEO SAR***对观测目标在全部照射时间段内的测量精度系数最小定位精度系数,进行最优组合选择:
在所述分布式GEO SAR***对所述观测目标的全部照射轨道范围Sa上,获得最优的三个观测角度组合,其中,SA1,SA2,SA3分别为各观测角度对应的轨道孔径中心时刻的差分干涉测量的等效孔径中心位置,P0为观测场景中心,PDOPd(SA1,SA2,SA3,P0)为分布式SAR***在SA1,SA2,SA3位置对P0进行三维形变反演的测量精度系数。
本发明实施例进一步给出一种分布式GEO SAR三维形变反演多角度选取***,包括:
处理器,用于执行多条指令;
存储器,用于存储多条指令;
其中,所述多条指令,用于由所述存储器存储,并由所述处理器加载并执行如前所述的分布式GEO SAR三维形变反演多角度选取方法。
本发明实施例进一步给出一种计算机可读存储介质,所述存储介质中存储有多条指令;所述多条指令,用于由处理器加载并执行如前所述的分布式GEO SAR三维形变反演多角度选取方法。
需要说明的是,在不冲突的情况下,本发明中的实施例及实施例中的特征可以相互组合。
在本发明所提供的几个实施例中,应该理解到,所揭露的***,装置和方法,可以通过其它的方式实现。例如,以上所描述的装置实施例仅仅是示意性的,例如,所述单元的划分,仅仅为一种逻辑功能划分,实际实现时可以有另外的划分方式,例如,多个单元或组件可以结合或者可以集成到另一个***,或一些特征可以忽略,或不执行。另一点,所显示或讨论的相互之间的耦合或直接耦合或通信连接可以是通过一些接口,装置或单元的间接耦合或通信连接,可以是电性,机械或其它的形式。
所述作为分离部件说明的单元可以是或者也可以不是物理上分开的,作为单元显示的部件可以是或者也可以不是物理单元,即可以位于一个地方,或者也可以分布到多个网络单元上。可以根据实际的需要选择其中的部分或者全部单元来实现本实施例方案的目的。
另外,在本发明各个实施例中的各功能单元可以集成在一个处理单元中,也可以是各个单元单独物理存在,也可以两个或两个以上单元集成在一个单元中。上述集成的单元既可以采用硬件的形式实现,也可以采用硬件加软件功能单元的形式实现。
上述以软件功能单元的形式实现的集成的单元,可以存储在一个计算机可读取存储介质中。上述软件功能单元存储在一个存储介质中,包括若干指令用以使得一台计算机装置(可以是个人计算机,实体机服务器,或者网络云服务器等,需安装Ubuntu操作***)执行本发明各个实施例所述方法的部分步骤。而前述的存储介质包括:U盘、移动硬盘、只读存储器(Read-Only Memory,ROM)、随机存取存储器(Random Access Memory,RAM)、磁碟或者光盘等各种可以存储程序代码的介质。
以上所述,仅是本发明的较佳实施例而已,并非对本发明作任何形式上的限制,依据本发明的技术实质对以上实施例所作的任何简单修改、等同变化与修饰,均仍属于本发明技术方案的范围内。
Claims (10)
1.一种分布式GEO SAR三维形变反演多角度选取方法,其特征在于,所述方法包括以下步骤:
步骤S101:获取分布式GEO SAR***对观测目标全部照射的空间范围内的速度、坐标、场景坐标;其中场景中心记为P0,根据场景中心P0计算正北、正东和垂直三个方向的单位矢量选取任意三个观测角度,获取各观测角度对应的轨道孔径中心时刻的差分干涉测量的等效孔径中心位置SAn,计算轨道孔径中心时刻的双基地角βn和GEO卫星到目标斜距的单位矢量/>得到三维形变模型系数矩阵Θ,其中n=1,2,3;
步骤S102:由干涉图的多视视数N和相关系数γ计算各观测角度下的干涉图相位方差,通过各观测角度下的干涉图相位方差获得相位误差的协方差矩阵CΦ;
步骤S103:根据所述三维形变模型系数矩阵Θ和所述相位误差的协方差矩阵CΦ,计算所述分布式GEO SAR***对观测目标在全部照射时间段内的测量精度系数最小定位精度系数;
步骤S104:基于所述分布式GEO SAR***对观测目标在全部照射时间段内的测量精度系数最小定位精度系数,进行最优组合选择:
在所述分布式GEO SAR***对所述观测目标的全部照射轨道范围Sa上,获得最优的三个观测角度组合,其中,SA1,SA2,SA3分别为各观测角度对应的轨道孔径中心时刻的差分干涉测量的等效孔径中心位置,P0为观测场景中心,PDOPd(SA1,SA2,SA3,P0)为分布式SAR***在SA1,SA2,SA3位置对P0进行三维形变反演的测量精度系数。
2.如权利要求1所述的分布式GEO SAR三维形变反演多角度选取方法,其特征在于,所述三维形变模型系数矩阵Θ,计算方式为:
其中λ为雷达波长,<·,·>表示向量内积。
3.如权利要求2所述的分布式GEO SAR三维形变反演多角度选取方法,其特征在于,所述步骤S102,所述相位误差的协方差矩阵CΦ表示为
其中,各观测角度下的干涉图相位方差由对应干涉图的多视视数N和相关系数γ计算,
4.如权利要求3所述的分布式GEO SAR三维形变反演多角度选取方法,其特征在于,所述步骤S103计算所述分布式GEO SAR***对观测目标在全部照射时间段内的测量精度系数最小定位精度系数PDOPd的方式:
在区分各个观测角度下生成干涉图的相位误差大小时,
在不区分各个观测角度下生成干涉图的相位误差大小时,
其中,tr(·)为矩阵的迹。
5.一种分布式GEO SAR三维形变反演多角度选取装置,其特征在于,所述装置包括:
系数矩阵获取模块:配置为获取分布式GEO SAR***对观测目标全部照射的空间范围内的速度、坐标、场景坐标;其中场景中心记为P0,根据场景中心P0计算正北、正东和垂直三个方向的单位矢量选取任意三个观测角度,获取各观测角度对应的轨道孔径中心时刻的差分干涉测量的等效孔径中心位置SAn,计算轨道孔径中心时刻的双基地角βn和GEO卫星到目标斜距的单位矢量/>得到三维形变模型系数矩阵Θ,其中n=1,2,3;
协方差矩阵获取模块:配置为由干涉图的多视视数N和相关系数γ计算各观测角度下的干涉图相位方差,通过各观测角度下的干涉图相位方差获得相位误差的协方差矩阵CΦ;
精度系数计算模块:配置为根据所述三维形变模型系数矩阵Θ和所述相位误差的协方差矩阵CΦ,计算所述分布式GEO SAR***对观测目标在全部照射时间段内的测量精度系数最小定位精度系数;
确定观测角度模块:配置为基于所述分布式GEO SAR***对观测目标在全部照射时间段内的测量精度系数最小定位精度系数,进行最优组合选择:
在所述分布式GEO SAR***对所述观测目标的全部照射轨道范围Sa上,获得最优的三个观测角度组合,其中,SA1,SA2,SA3分别为各观测角度对应的轨道孔径中心时刻的差分干涉测量的等效孔径中心位置,P0为观测场景中心,PDOPd(SA1,SA2,SA3,P0)为分布式SAR***在SA1,SA2,SA3位置对P0进行三维形变反演的测量精度系数。
6.如权利要求5所述的分布式GEO SAR三维形变反演多角度选取装置,其特征在于,所述三维形变模型系数矩阵Θ,计算方式为:
其中λ为雷达波长,<·,·>表示向量内积。
7.如权利要求6所述的分布式GEO SAR三维形变反演多角度选取装置,其特征在于,所述相位误差的协方差矩阵CΦ表示为
其中,各观测角度下的干涉图相位方差由对应干涉图的多视视数N和相关系数γ计算,
8.如权利要求7所述的分布式GEO SAR三维形变反演多角度选取装置,其特征在于,所述计算所述分布式GEO SAR***对观测目标在全部照射时间段内的测量精度系数最小定位精度系数PDOPd的方式:
在区分各个观测角度下生成干涉图的相位误差大小时,
在不区分各个观测角度下生成干涉图的相位误差大小时,
其中,tr(·)为矩阵的迹。
9.一种分布式GEO SAR三维形变反演多角度选取***,包括:
处理器,用于执行多条指令;
存储器,用于存储多条指令;
其中,所述多条指令,用于由所述存储器存储,并由所述处理器加载并执行如权利要求1-4中任一项所述的分布式GEO SAR三维形变反演多角度选取方法。
10.一种计算机可读存储介质,所述存储介质中存储有多条指令;所述多条指令,用于由处理器加载并执行如权利要求1-4中任一项所述的分布式GEO SAR三维形变反演多角度选取方法。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110921920.7A CN113777606B (zh) | 2021-08-12 | 2021-08-12 | 分布式geo sar三维形变反演多角度选取方法及装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110921920.7A CN113777606B (zh) | 2021-08-12 | 2021-08-12 | 分布式geo sar三维形变反演多角度选取方法及装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113777606A CN113777606A (zh) | 2021-12-10 |
CN113777606B true CN113777606B (zh) | 2023-12-26 |
Family
ID=78837502
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110921920.7A Active CN113777606B (zh) | 2021-08-12 | 2021-08-12 | 分布式geo sar三维形变反演多角度选取方法及装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113777606B (zh) |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105607058A (zh) * | 2015-12-24 | 2016-05-25 | 中国科学院电子学研究所 | 利用geosar相位定标信息进行短弧段轨道精化的方法 |
CN105866777A (zh) * | 2016-03-29 | 2016-08-17 | 北京理工大学 | 多角度多时段导航卫星双基地PS-InSAR三维形变反演方法 |
CN112444188A (zh) * | 2021-02-01 | 2021-03-05 | 深圳大学 | 一种多视角InSAR海堤高精度三维形变测量方法 |
CN112862924A (zh) * | 2020-12-31 | 2021-05-28 | 哈尔滨工业大学(威海) | 多模态成像中图像重建方法、装置和多模态成像技术*** |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110058236B (zh) * | 2019-05-21 | 2023-04-07 | 中南大学 | 一种面向三维地表形变估计的InSAR和GNSS定权方法 |
-
2021
- 2021-08-12 CN CN202110921920.7A patent/CN113777606B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105607058A (zh) * | 2015-12-24 | 2016-05-25 | 中国科学院电子学研究所 | 利用geosar相位定标信息进行短弧段轨道精化的方法 |
CN105866777A (zh) * | 2016-03-29 | 2016-08-17 | 北京理工大学 | 多角度多时段导航卫星双基地PS-InSAR三维形变反演方法 |
CN112862924A (zh) * | 2020-12-31 | 2021-05-28 | 哈尔滨工业大学(威海) | 多模态成像中图像重建方法、装置和多模态成像技术*** |
CN112444188A (zh) * | 2021-02-01 | 2021-03-05 | 深圳大学 | 一种多视角InSAR海堤高精度三维形变测量方法 |
Also Published As
Publication number | Publication date |
---|---|
CN113777606A (zh) | 2021-12-10 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Zhang et al. | Modeling PSInSAR time series without phase unwrapping | |
Pu et al. | Motion errors and compensation for bistatic forward-looking SAR with cubic-order processing | |
CN111538006A (zh) | 基于动态基线的InSAR数字高程模型构建方法及*** | |
CN106526590A (zh) | 一种融合多源sar影像工矿区三维地表形变监测及解算方法 | |
US7511655B2 (en) | Method and apparatus for 3-D sub-voxel position imaging with synthetic aperture radar | |
US7446705B1 (en) | Method and apparatus for determining parameters for a parametric expression characterizing the phase of an acquired signal | |
CN112986949B (zh) | 针对角反射器的sar高精度时序形变监测方法和装置 | |
Chatterjee et al. | Parallax and kinematics of PSR B0919+ 06 from VLBA astrometry and interstellar scintillometry | |
Mrstik et al. | Terrain height measurement accuracy of interferometric synthetic aperture radars | |
Baumgartner et al. | Dual-platform large along-track baseline GMTI | |
MacGougan et al. | Accuracy and reliability of tightly coupled GPS/ultra-wideband positioning for surveying in urban environments | |
CN105487094A (zh) | 一种数据链与卫星导航协同定位方法及定位*** | |
Hellwich et al. | Geocoding SAR interferograms by least squares adjustment | |
Nico | Exact closed-form geolocation for SAR interferometry | |
Ma et al. | Time synchronization requirement of global navigation satellite system augmentation system based on pseudolite | |
CN113777606B (zh) | 分布式geo sar三维形变反演多角度选取方法及装置 | |
GB2446615A (en) | Interleaved bi-static interferometric synthetic aperture radar technique for determining height information for an imaged area | |
RU2670976C1 (ru) | Способ определения местоположения источника радиоизлучения с периодической структурой сигнала и вращающейся направленной антенной | |
Burns et al. | IFSAR for the rapid terrain visualization demonstration | |
Zhang et al. | Spatio-temporal linking of multiple SAR satellite data from medium and high resolution Radarsat-2 images | |
RU2556699C1 (ru) | Способ пеленгации источников радиоизлучения на одной частоте | |
Baumann-Ouyang et al. | Estimating 3D displacement vectors from line-of-sight observations with application to MIMO-SAR | |
CN112346028A (zh) | 分布式InSAR卫星方位向空间同步在轨测试方法及*** | |
CN112505693A (zh) | 干涉逆合成孔径雷达成像配准方法、***和存储介质 | |
Huang | Statistical theory for the detection of persistent scatterers in InSAR imagery |
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 |