一种时序星载雷达数据处理方法和装置的制造方法
一种时序星载雷达数据处理方法和装置的制造方法
【技术领域】
[0001] 本发明设及干设测量技术领域，尤其设及一种时序星载雷达数据处理方法和装 置。
【背景技术】
[0002] 自上世纪八十年代W来，基于星载雷达数据的干设测量技术已被广泛用于全天 候、全天时、大范围和高精度地监测地表及地面设施的形变，包括桥梁、大巧、机场及建筑综 合体等。目前，InSAR(Synthetic Aperture Radar Interferomehy，合成孔径雷达干设）测 量技术已经成为监测及评估地震、火山及滑坡等自然灾害的重要手段之一，在环境恶劣区 域，甚至是唯一手段。其干设技术原理可简单概括为通过比较在同一地区相同成像几何条 件下，不同时刻获取的复数雷达数据的相位信息来推算地表的形变。雷达卫星一般位于高 度为800千米左右的轨道上，由于当前定轨技术精度的限制，对同一地区卫星的重访位置难 W-致，从而造成在对地物目标两次成像时存在入射角度的差异，该差异可W用两次成像 时卫星之间的距离来描述，运个距离通常被称之为基线。基线长度的精确测定对利用干设 测量技术获取地面高程及形变具有至关重要的作用。直观上，如果知道两次卫星的坐标即 可算出基线长度，但由于定轨技术精度所限，所确定的卫星坐标通常会存在数厘米或数十 厘米的误差。运种轨道误差会传播到基线上，形成基线误差，最终体现到干设相位中，从而 成为InSAR技术的误差源之一。
[0003] 目前已有多种方法来估计基线误差。由于基线误差会通过平地效应和高程传播到 干设相位中，现有方法均通过对相位的分析来计算轨道误差。首先假设干设相位成功解缠， 即不存在相位模糊度，并且观测区的形变、地形残差和大气延迟误差可W忽略不计，运时解 缠相位的贡献主要与基线的平行雷达视线和垂直雷达视线运两个方向的分量有关，通过在 成像区域均匀选取高程已知的控制点，建立观测方程，利用最小二乘来解算影像中屯、点的 基线和基线变化率。但是运种方法需要干设相位的成功解缠，但事实上解缠的成功率在实 际中无法保证，特别是在低相干区域，例如植被区、地形复杂区及高楼林立的城区，从而影 响到基线估计的精确度，另外，该方法需要对每个干设基线进行改正，在利用InSAR技术分 析海量雷达干设数据时效率低下。
【发明内容】
[0004] 本发明的目的在于提供一种时序星载雷达数据处理方法和装置，旨在解决现有技 术中因相位误差，如解缠误差、长波形变信号、大气延迟误差W及地形残差等干扰而导致的 基线估计的精确度不高，W及在利用InSAR技术分析海量雷达干设数据时需对每一个干设 图进行基线精化导致的效率低下的问题。
[0005] 本发明第一方面，提供一种时序星载雷达数据处理方法，包括：
[0006] 根据获取的时序星载雷达数据，构建时序干设测量模型，所述时序星载雷达数据 包括雷达影像，所述时序干设测量模型用于计算目标卫星的轨道误差和高程误差；
[0007] 根据所述时序干设测量模型，计算所述轨道误差和所述高程误差；
[0008] 基于所述轨道误差和所述高程误差重构基线误差相位，所述基线误差相位用于确 定基线及基线的变化率。
[0009] 本发明第二方面，提供一种时序星载雷达数据处理装置，包括：。
[0010] 模型构建模块，用于根据获取的时序星载雷达数据，构建时序干设测量模型，所述 时序星载雷达数据包括雷达影像，所述时序干设测量模型用于计算目标卫星的轨道误差和 高程误差；
[0011] 误差计算模块，用于根据所述时序干设测量模型，计算所述轨道误差和所述高程 误差；
[0012] 重构模块，用于基于所述轨道误差和所述高程误差重构基线误差相位，所述基线 误差相位用于确定基线及基线的变化率。
[0013] 本发明与现有技术相比存在的有益效果是:本发明提供的技术方案根据时序星载 雷达数据，构建能够同时计算轨道误差和高程误差的时序干设测量模型，从而实现从缠绕 的雷达干设相位中直接精确分离出轨道误差和高程误差，并基于轨道误差和高程误差重构 基线误差相位，从而提高基线估计的精确度和效率。
【附图说明】
[0014] 图1是本发明实施例一提供的时序星载雷达数据处理方法的流程图；
[0015] 图2是本发明实施例一提供的时序星载雷达数据处理方法中轨道误差对基线估计 的影响的示意图；
[0016] 图3是本发明实施例一提供的时序星载雷达数据处理方法中雷达成像的几何示意 图；
[0017] 图4是本发明实施例一提供的时序星载雷达数据处理方法中典型的干设图中的基 线误差示意图；
[0018] 图5是本发明实施例二提供的时序星载雷达数据处理方法的流程图；
[0019] 图6是本发明实施例Ξ提供的时序星载雷达数据处理装置的组成示意图；
[0020] 图7是本发明实施例四提供的时序星载雷达数据处理装置的组成示意图。
【具体实施方式】
[0021] 为了使本发明的目的、技术方案及优点更加清楚明白，W下结合附图及实施例，对 本发明进行进一步详细说明。应当理解，此处所描述的具体实施例仅仅用W解释本发明，并 不用于限定本发明。
[0022] W下结合具体附图对本发明的实现进行详细的描述。
[0023] 实施例一：
[0024] 图1是本发明实施例一提供的时序星载雷达数据处理方法的流程图，具体包括步 骤S101至S103,详述如下：
[0025] S101、根据获取的时序星载雷达数据，构建时序干设测量模型，该时序星载雷达数 据包括雷达影像，该时序干设测量模型用于计算目标卫星的轨道误差和高程误差。
[0026] 载有合成孔径雷达的卫星可在800千米的轨道上对地球进行重复观测，地面站通 过星载GPS(Global化sitioning System,全球定位系统）等设备对卫星进行追踪及定位， 但由于定位精度的限制，星载GPS所确定的卫星位置不可避免存地在误差，两次卫星观测位 置间的距离即为基线。显然地，卫星的轨道误差会传播至基线上产生基线误差。图2示出了 轨道误差对基线估计的影响，其中S1和S2表示卫星对地观测时的位置，Bo为S1和S2之间的 矢量距离即基线，α为观测基线与水平方向夹角，β为误差基线与水平方向夹角，Θ为雷达波 入射角。真实的基线矢量Β等于Bo和由轨道误差引起的基线误差矢量η之和。
[0027] 具体地，根据研究区域内所获取的所有时序星载雷达数据，构建一个时序干设测 量模型，该时序干设测量模型能够将轨道误差和高程误差同时解算出来。
[0028] 考虑到海量的雷达数据，对于同一研究区域，每一期影像通常有数百万个有效监 测点，若数十期影像的数据同时处理，数据量将极为庞大。为了解决相位模糊度问题W及抑 制大气延迟误差，需要对有效监测点进行采样，并对采样点进行空间差分处理，将采样点的 差分相位作为基本观测参照值，可W采用顾及监测密度和相位梯度的四叉树采样来确定基 本观测参照值。基本观测参照值确定之后，需要建立待估参数与基本观测参照值之间的关 系，即时序干设测量模型，其中待估参数包括采样点的形变序列、地形残差和轨道误差，其 中地形残差包括了高程误差。为了提高计算效率，可W根据干设图及雷达影像的线性组合 关系，将干设图轨道误差转为雷达影像的虚拟轨道误差，即将干设图中的轨道误差分量分 解到每一个雷达影像上，通过获取雷达影像的虚拟轨道误差参数来重构干设图的轨道误 差，建立待估参数与基本观测参照值之间的关系，作为时序干设测量模型。
[0029] S102、根据时序干设测量模型，计算轨道误差和高程误差。
[0030] 具体地，步骤S101得到的时序干设测量模型可W是一个超大型线性稀疏系统，采 用共辆梯度法，通过迭代方式完成对模型中待估参数的解算，从而分离出轨道误差和高程 误差。
[0031] S103、基于轨道误差和高程误差重构基线误差相位，该基线误差相位用于确定基 线及基线的变化率。
[0032] 由于InSAR测量技术是一种利用相位测距的技术，是通过测量两次卫星到地面目 标的距离，即斜距的变化来推算目标高程或形变。图3示出了卫星雷达成像的几何示意图， 其中S1和S2表示卫星对地观测时的位置，B为S1和S2之间的距离即基线，P为地表目标即相 干点，P'为地表目标在楠球面的投影，Z为目标点的高程，Η为卫星高，0代表楠球原点，α为基 线方向角，Θ为目标点对应的雷达波入射角，θ〇为忽略目标点高程的雷达波入射角，Ρ为卫星 至地面目标的距离，阴影区域表示地形。可W看到地面目标Ρ的高程及形变都会产生斜距 差，同样地，基线Β的基线误差也会产生斜距差，运些斜距差反映到卫星雷达的干设相位上 则显示为干设条纹，如图4所示。在图4显示的干设图中可W明显看到因基线误差引起的干 设条纹，条纹越密说明基线误差越大，该示例数据来自于C波段的化Vi sat/ASAR星载雷达， 成像区为罗布泊。因此，基线误差相位可W有两部分组成，包括平行基线相位和高程作用下 的垂直基线相位。
[0033] 具体地，在重构基线误差相位时，根据时序干设测量模型计算得到的轨道误差和 高程误差，分别重构水平基线相位和垂直基线相位，从而完成基线误差相位的重构。对重构 后的基线误差相位，可W利用最小二乘原理来确定基线W及基线的变化率。
[0034] 本实施例中，根据时序星载雷达数据，构建能够同时计算轨道误差和高程误差的 时序干设测量模型，从而实现从缠绕的雷达干设相位中直接精确分离轨道误差和高程误 差，并基于轨道误差和高程误差重构基线误差相位，从而提高基线估计的精确度和效率。 [003引实施例二：
[0036] 图5 是本发明实施例二提供的时序星载雷达数据处理方法的流程图，具体包括步 骤S201至S204,详述如下：
[0037] S201、根据获取的时序星载雷达数据，确定时序干设测量模型的基本观测参照值， 该时序星载雷达数据包括雷达影像，该时序干设测量模型用于计算目标卫星的轨道误差和 高程误差。
[0038] 具体地，确定时序干设测量模型的基本观测参照值可W通过步骤S2011至S2016实 现，详细说明如下：
[0039] S2011、获取研究区域的所有时序星载雷达数据并进行干设组合得到干设对。
[0040] 具体地，将研究区域内所获取的所有时序星载雷达数据进行干设组合，得到干设 对，并计算干设对的时间基线和空间基线，根据基线信息对干设对进行初选。
[0041] 干设对表示两幅雷达影像的相位差，假设有N+1幅雷达影像，最多可形成的干设组 合数量为〇.5XNX(N+l)。时间基线为构成干设对的两幅雷达影像的时间差。空间基线为卫 星在获取构成干设对的两幅影像时的位置差异。时间基线和空间基线可通过雷达影像的参 数文件计算得到。
[0042] 在进行干设对初选时，将所有干设对按基线信息在二维空间进行排列，用具备等 步长增长的捜索窗从零开始捜索，当捜索窗内所包含的干设对包含所有影像并且干设图数 量最小时停止捜索。
[0043] S2012、对干设对进行采样和相干性计算，并统计相干性高于预设阔值的采样点和 采样点的数量。
[0044] 具体地，针对步骤S2011中初选得到的干设对，将包含于运些干设对中的雷达影像 进行基于灰度的四叉树采样，并对采样数据进行相干性计算，然后统计各干设对的相干性 高于预设阔值的采样点和采样点的数量。
[0045] 初选的干设对数量可能很大，为了提高计算效率，需要进一步筛选高质量的干设 对。在雷达干设测量技术中一般采用相干性来衡量干设对的质量，但逐个像素计算相干性 比较费时，因此可W先通过四叉树方法对雷达强度影像图根据灰度信息进行四叉树采样， 然后计算采样点的相干性，用有限采样点的相干性来代表整幅干设图的质量。
[0046] 基于灰度的四叉树采样方法首先把雷达强度影像图等分成四个分区，如果检查到 某个分区所有格网的灰度值之间的差值小于或等于预置数值，则该某个分区就不再往下分 害d;否则，把该某个分区再分割成四个子区域，运样递归地分割，直至每个分区所有格网的 灰度值之间的差值小于或等于预置数值为止。
[0047] 相干性是一个介于0和1之间的值，可由如下公式计算：
[004引
[0049]其中，Cl和C2为同一地面目标在两次卫星获取影像时所反射的雷达波信号，该信号 为复数值，包含了地面目标对应的强度和相位信息，(Λ'为C2的复数共辆，Φ<1β*为地形和形变 贡献的相位，为取模算子。
[0050] S2013、根据采样点和采样点的数量，构建最小生成树。
[0051] 具体地，根据步骤S2012得到的采样点的相干性信息和采样点的数量，构建最小生 成树(Minimum Spanning化ee，MST)。该生成树可W在N景雷达影像中筛选出质量最优的N-1个干设对，其中N表示雷达影像数量，时序干设测量模型的基本观测参照值将从运些干设 对中选取。
[0052] 生成树是连通图的极小连通子图。极小是指若在树中任意增加一条边，则将出现 一个回路，若去掉一条边，将会使之变成非连通图。生成树各边的权值总和称为生成树的 权。权最小的生成树称为最小生成树。
[0053] 进一步地，将步骤S2012得到的采样点的相干性均值的倒数作为权进行MST捜索， 按照预定的构建最小生成树的方法构建最小生成树。
[0054] S2014、对包含于最小生成树中的干设对进行相干性评估，选取时域相干点。
[0055] 具体地，为了减少相干性估计的偏差，在对包含于最小生成树中的干设对进行相 干性评估时，首先要对采样窗的像素点进行同质点选取，即选取与中屯、像素点散射特征的 偏差在预定范围内的周边像素点，再计算所选取的像素点的相干性，并根据所计算的相干 性选取相干性小于经验性阔值的像素点，作为时域相干点。
[0056] 假设采样窗内有Μ个同质像素点，则相干性可W由W下公式计算：
[0化7]
[0058] 其中，Cl和C2为同一地面目标在两次卫星获取影像时所反射的雷达波信号，该信号 为复数值，包含了地面目标对应的强度和相位信息，C为C2的复数共辆，Φ det为地形和形变 贡献的相位，1[1为取模算子。
[0059] 利用上述公式所计算的相干性，通过设定一个经验性阔值可将时域相干点筛选出 来，经验性阔值可W设置为0.1。
[0060] S2015、对时域相干点根据空间分布和密度进行采样，得到采样相干点。
[0061] 具体地，对步骤S2014选出的时域相干点根据空间分布和密度进行四叉树采样，得 到采样相干点。通常选取一万个采样相干点已经足够恢复出轨道相位，即当四叉树划分的 窗口数量达到一万时可W停止操作。
[0062] S2016、对采样相干点进行空间高密度构网，并将网的每条弧段的缠绕相位差作为 时序干设测量模型的基本观测参照值。
[0063] 具体地，对步骤S2015得到的采样相干点进行空间高密度构网，网的每条边，即每 条弧段的缠绕相位差作为时序干设测量模型的基本观测参照值。
[0064] 高密度构网是把采样点密集地连接起来。考虑到大气延迟的空间相关性，将研究 区域划成1km X 1 km的格网，W1 km X 1 km格网的顶点为中屯、、1 km为半径的范围内捜索采样 点，将捜索到的采样点构建为Delaunays角网，实现采样点的连接。
[0065] 缠绕相位差即采样点的相位差。采样点反射的雷达波的原始干设相位Φ永远介于 (-η，π]之间，总是与实际的相位差如的整数倍，所W称为缠绕相位。假设两个采样点上的缠 绕相位分别为Φ谢Φ2,则缠绕相位差
，其中mod( ·)为余数算子。
[0066] S202、建立待估参数与基本观测参照值之间的联系，作为时序干设测量模型。
[0067] 具体地，待估参数主要包括采样相干点的形变序列、地形残差W及轨道误差，其中 地形残差包括了高程误差，通过建立待估参数与基本观测参照值即缠绕相位差之间的观测 方程，体现待估参数与基本观测参照值之间的联系，作为时序干设测量模型。
[0068] 观测方程的具体建立过程如下：
[0069] -方面，轨道误差贡献的差分相位与轨道参数之间的关系为： 剛
（3)
[0071] 其中，Δ为差分算子，N为N+1景影像形成的N个干设对，G为采样相干点构建的弧段 的个数，Dorb为系数阵，Porb为参数阵；
[0072] 系数阵Dnrb的计算公式为： W73]
(4)
[0074]其中，A表示获取的干设图与雷达影像的构件关系，D表示弧段上轨道多项式参数 对应的系数矩阵；
[0077] 其中，dX、dY、dXY分别为弧段上轨道多项式参数对应的系数;j表示弧段序号；
[0078] 参数阵Pnrb的计算公式为：
[0079]
所
[0080] 其中，w、bj、Cj为雷达影像的虚拟轨道误差参数。
[0081] 另一方面，在差分相位中除了轨道误差的贡献，还有高程和地表形变的贡献。具体 地，在第i个干设图中，第P个采样点的地形残差相位与该点的高程误差A hp存在确定 性关系，可表示如下：
[0082]
巧）
[0083] 其中，β?为一个与基线、入射角及斜距有关的系数。
[0084] 类似地，可W建立地表形变序列贡献的相位祐6，,Ρ与形变Vp的关系，可表示如下：
[0085]
(9)
[0086] 其中，钱是一个与时间有关的系数。
[0087] 则根据公式(8)和公式(9)，可W得到地形残差及形变序列共同贡献的相位为：
[008引 (10)
[0089] 其中B为系数阵，N为干设图数量。
[0090] 根据公式（10)可进一步得到采样点间的相位差矢量与地形残差及形变序列的关 系为： 幽]
(11)
[0092] 其4
C为反映弧段和采样相干点之间关系 的邻接矩阵，Ppar为包含了采样点上的地形残差及形变序列参数的参数阵。
[0093] 综合公式(3)和公式（11)，缠绕相位差与轨道误差，地形残差及形变序列之间的观 测方程可W表示为：
[0094]
(12)
[0095] 其中，D=[D0rb Dpar]T，P=[P0rb Ppar]T，W为缠绕相位差中除了轨道误差、地形残差 及形变序列外的噪声相位矢量。
[0096] 该观测方程（12)作为时序干设测量模型，体现了待估参数与基本观测参照值之间 的联系。在进行轨道误差参数化时，为了减少参数个数，考虑到干设图和雷达影像的线性组 合关系，将干设图中的轨道误差分量分解到每个雷达影像上，一旦获得雷达影像的虚拟轨 道误差参数，即可重构任何一个干设图的轨道误差。干设图轨道误差可W是相邻弧段的雷 达影像虚拟轨道误差的代数差。
[0097] 具体地，通过雷达影像的虚拟轨道误差参数重构干设图的轨道误差的过程如下：
[0098] 将雷达影像1和雷达影像2表示为相邻弧段的雷达影像；
[0099] 雷达影像1的虚拟轨道误差多项式为：
[0100]
(巧) [0101] 其中，ai、bi、ci为雷达影像1的虚拟轨道误差参数，X和y分别是W雷达影像左上角 点为坐标原点的像素横坐标和像素纵坐标；
[0102] 雷达影像2的虚拟轨道误差多项式为：
[01。引
（14)
[0104] 其中，a2、b2、C2为雷达影像2的虚拟轨道误差参数，X和y分别是W雷达影像左上角 点为坐标原点的像素横坐标和像素纵坐标。
[0105] 则根据多项式(13)和多项式(14)，可W得到干设图的轨道误差多项式为：
[0106]
(15)
[0107] 其中，a、b、c为干设图轨道误差参数，X和y分别是W雷达影像左上角点为坐标原点 的像素横坐标和像素纵坐标。
[010引通过观测方程解算出多项式（13)的参数ai、bi、C1和多项式（14)的参数日2、b2、C2，可 W得到多项式(15)的参数a、b、c，进而得到干设图的轨道误差分量。
[0109 ] S203、根据时序干设测量模型，计算轨道误差和高程误差。
[0110] 具体地，步骤S202中构建的观测方程是一个超大型线性稀疏系统，可W采用共辆 梯度法，通过迭代方式实现对参数的解算，进而从干设相位中分离出轨道误差分量及高程 误差分量。
[0111] S204、基于轨道误差和高程误差重构基线误差相位，该基线误差相位用于确定基 线及基线的变化率。
[0112] 由于InSAR测量技术是一种利用相位测距的技术，是通过测量两次卫星到地面目 标的距离，即斜距的变化来推算目标高程或形变。图3示出了卫星雷达成像的几何示意图， 其中S1和S2表示卫星对地观测时的位置，B为S1和S2之间的距离即基线，P为地表目标即相 干点，P'为地表目标在楠球面的投影，Z为目标点的高程，Η为卫星高，0代表楠球原点，α为基 线方向角，Θ为目标点对应的雷达波入射角，θ〇为忽略目标点高程的雷达波入射角，Ρ为卫星 至地面目标的距离，阴影区域表示地形。可W看到地面目标Ρ的高程及形变都会产生斜距 差，同样地，基线Β的基线误差也会产生斜距差，运些斜距差反映到卫星雷达的干设相位上 则显示为干设条纹，如图4所示。在图4显示的干设图中可W明显看到因基线误差引起的干 设条纹，条纹越密说明基线误差越大，该示例数据来自于C波段的化Vi sat/ASAR星载雷达， 成像区为罗布泊。因此，基线误差相位可W有两部分组成，包括平行基线相位和高程作用下 的垂直基线相位。
[0113] 具体地，在重构基线误差相位时，首先将时序干设测量模型计算的高程误差改正 到外部高程数据上，W改善外部高程数据的精度，然后根据采样点的垂直基线信息重构垂 直基线相位。
[0114] 定义垂直基线相位重构公式为：
[0115]
(1食)
[0116] 其中，λ为雷达波长，P为卫星到给定地面点的距离，Θ为该给定地面点所对应的雷 达波束入射角，Β丄为该给定地面点对应的垂直基线，h为将高程误差改正到外部高程数据上 后得到的该给定地面点精化后的高程。
[0117] 另一方面，平行基线贡献的相位仅与基线有关，且在空线上任意一点可用一个多 项式代替，即平行基线相位重构公式可定义为：
[0118]
(17；
[0119] 其中，λ为雷达波长，B=为所述地面点对应的平行基线，X和y分别是W雷达影像左 上角点为坐标原点的像素横坐标和像素纵坐标，a、b、c为干设图轨道误差参数，可由虚拟轨 道误差多项式线性组合得到。
[0120] 通过垂直基线相位重构公式（16)和平行基线相位重构公式（17)即可完成基线误 差相位的重构。对重构后的基线误差相位，可W利用最小二乘原理来确定基线W及基线的 变化率。
[0121 ]本实施例中，首先对获取的时序星载雷达数据进行干设组合，将得到的干设对进 行采样和相干性计算后，对所构建的最小生成树进行相干性评估，选取时域相干点，并对时 域相干点进行空间分布和密度的四叉树采样得到采样相干点，对采样相干点进行空间高密 度构网得到时序干设测量模型的基本观测参照值，并建立待估参数和基本观测参照值之间 的联系，作为时序干设测量模型，通过解算该模型实现从缠绕的雷达干设相位中直接精确 分离出轨道误差和高程误差，并基于轨道误差和高程误差重构基线误差相位，从而提高基 线估计的精确度和效率。
[0122] 实施例
[0123] 图6是本发明实施例Ξ提供的时序星载雷达数据处理装置的组成示意图，为了便 于说明，仅示出了与本发明实施例相关的部分。图8示例的用于时序星载雷达数据处理的装 置可W是前述实施例一提供的用于时序星载雷达数据处理的方法的执行主体，其可W是应 用设备或者应用设备中的一个功能模块。图8示例的用于时序星载雷达数据处理的装置主 要包括:模型构建模块81、误差计算模块82、重构模块83。各功能模块详细说明如下：
[0124] 模型构建模块81，用于根据获取的时序星载雷达数据，构建时序干设测量模型，所 述时序星载雷达数据包括雷达影像，所述时序干设测量模型用于计算目标卫星的轨道误差 和高程误差；
[0125] 误差计算模块82，用于根据所述时序干设测量模型，计算所述轨道误差和所述高 程误差；
[0126] 重构模块83,用于基于所述轨道误差和所述高程误差重构基线误差相位，所述基 线误差相位用于确定基线及基线的变化率。
[0127] 本实施例中的时序星载雷达数据处理的装置中，各模块实现各自功能的过程，参 见前述图1所示实施例的描述，此处不再寶述。
[0128] 从上述图8示例的用于时序星载雷达数据处理的装置可知，本实施例中，根据时序 星载雷达数据，构建能够同时计算轨道误差和高程误差的时序干设测量模型，从而实现从 缠绕的雷达干设相位中直接精确分离轨道误差和高程误差，并基于轨道误差和高程误差重 构基线误差相位，从而提高基线估计的精确度和效率。
[0129] 实施例四：
[0130] 图7是本发明实施例四提供的时序星载雷达数据处理装置的组成示意图，为了便 于说明，仅示出了与本发明实施例相关的部分。图9示例的用于时序星载雷达数据处理的装 置可W是前述实施例二提供的用于时序星载雷达数据处理的方法的执行主体，其可W是应 用设备或者应用设备中的一个功能模块。图9示例的用于时序星载雷达数据处理的装置主 要包括:模型构建模块91、误差计算模块92、重构模块93。各功能模块详细说明如下：
[0131] 模型构建模块91，用于根据获取的时序星载雷达数据，构建时序干设测量模型，所 述时序星载雷达数据包括雷达影像，所述时序干设测量模型用于计算目标卫星的轨道误差 和高程误差；
[0132] 误差计算模块92,用于根据所述时序干设测量模型，计算所述轨道误差和所述高 程误差；
[0133] 重构模块93,用于基于所述轨道误差和所述高程误差重构基线误差相位，所述基 线误差相位用于确定基线及基线的变化率。
[0134] 具体地，模型构建模块91，还包括：
[0135] 观测确定子模块911，用于根据获取的时序星载雷达数据，确定时序干设测量模型 的基本观测参照值，所述时序星载雷达数据包括雷达影像，所述时序干设测量模型用于计 算目标卫星的轨道误差和高程误差；
[0136] 联系建立子模块912,用于建立待估参数与所述基本观测参照值之间的联系，作为 时序干设测量模型。
[0137] 具体地，观测确定子模块911还用于：
[0138] 获取研究区域的所有时序星载雷达数据并进行干设组合得到干设对；
[0139] 对所述干设对进行采样和相干性计算，并统计相干性高于预设阔值的采样点和所 述采样点的数量；
[0140] 根据所述采样点和所述采样点的数量，构建最小生成树；
[0141 ]对包含于所述最小生成树中的干设对进行相干性评估，选取时域相干点；
[0142] 对所述时域相干点根据空间分布和密度进行采样，得到采样相干点；
[0143] 对所述采样相干点进行空间高密度构网，并将网的每条弧段的缠绕相位差作为时 序干设测量模型的基本观测参照值。
[0144] 具体地，联系建立子模块912还用于：
[0145] 所述待估参数包括形变序列、地形残差W及轨道误差，所述地形残差包括高程误 差；
[0146] 定义所述缠绕相位差与所述待估参数之间的观测方程为：
[0150] 其中，Δ为差分算子，N为N+1景所述雷达影像形成的N个干设对，G为所述采样相干 点构建的所述弧段的个数，W为所述缠绕相位差中除了所述形变序列、所述地形残差W及所 述轨道误差外的噪声相位矢量；
[0151]
C为反映所述弧段和所述采样相干点之间 关系的邻接矩阵，B为系数阵；
[0152] Ppar为包含了采样点上的所述形变序列和所述地形残差参数的参数阵；
[0153]
\表示获取的干设图与所述雷达影像的构件关系，D表示所述 弧段上轨道多项式参数对应的系数矩阵；
[0154]
j、bj、cj为所述雷达影像的虚拟轨道误差参数，j表示弧段序 号，相邻弧段的所述雷达影像的虚拟轨道误差参数的代数差为所述干设图的轨道误差参 数。具体地，重构模块93还用于：
[01W] 定义垂直基线相位重构公式为：
[0156]
[0157] 其中，λ为雷达波长，p为卫星到给定地面点的距离，θ为所述地面点所对应的雷达 波束入射角，Β丄为所述地面点对应的垂直基线，h为将所述高程误差改正到外部高程数据上 后得到的所述地面点精化后的高程；< br>[0158] 定义平行基线相位重构公式为：
[0159]
[0160] 其中，λ为雷达波长，B=为所述地面点对应的平行基线，X和y分别是W所述雷达影 像左上角点为坐标原点的像素横坐标和像素纵坐标，a、b、c为所述干设图轨道误差参数；
[0161] 通过所述垂直基线相位重构公式和所述平行基线相位重构公式完成基线误差相 位的重构。
[0162] 本实施例中的时序星载雷达数据处理的装置中，各模块实现各自功能的过程，参 见前述图5所示实施例的描述，此处不再寶述。
[0163] 从上述图9示例的用于时序星载雷达数据处理的装置可知，本实施例中，首先对获 取的时序星载雷达数据进行干设组合，将得到的干设对进行采样和相干性计算后，对所构 建的最小生成树进行相干性评估，选取时域相干点，并对时域相干点进行空间分布和密度 的四叉树采样得到采样相干点，对采样相干点进行空间高密度构网得到时序干设测量模型 的基本观测参照值，并建立待估参数和基本观测参照值之间的联系，作为时序干设测量模 型，通过解算该模型实现从缠绕的雷达干设相位中直接精确分离出轨道误差和高程误差， 并基于轨道误差和高程误差重构基线误差相位，从而提高基线估计的精确度和效率。
[0164] 需要说明的是，本说明书中的各个实施例均采用递进的方式描述，每一个实施例 重点说明的都是与其他实施例的不同之处，各个实施例之间相同或者相似的部分互相参见 即可。对于装置类实施例而言，由于其与方法实施例基本相似，所W描述的比较简单，相关 之处参见方法实施例的部分说明即可。
[0165] 值得注意的是，上述装置实施例中，所包括的各个模块只是按照功能逻辑进行划 分的，但并不局限于上述的划分，只要能够实现相应的功能即可；另外，各功能模块的具体 名称也只是为了便于相互区分，并不用于限制本发明的保护范围。
[0166] 本领域普通技术人员可W理解，实现上述各实施例方法中的全部或部分步骤是可 W通过程序来指令相关的硬件来完成，相应的程序可W存储于一计算机可读取存储介质 中，所述的存储介质，如R0M/RAM、磁盘或光盘等。
[0167] W上所述仅为本发明的较佳实施例而已，并不用W限制本发明，凡在本发明的精 神和原则之内所作的任何修改、等同替换和改进等，均应包含在本发明的保护范围之内。
【主权项】
1. 一种时序星载雷达数据处理方法，其特征在于，包括： 根据获取的时序星载雷达数据，构建时序干涉测量模型，所述时序星载雷达数据包括 雷达影像，所述时序干涉测量模型用于计算目标卫星的轨道误差和高程误差； 根据所述时序干涉测量模型，计算所述轨道误差和所述高程误差； 基于所述轨道误差和所述高程误差重构基线误差相位，所述基线误差相位用于确定基 线及基线的变化率。2. 根据权利要求1所述的时序星载雷达数据处理方法，其特征在于，所述根据获取的时 序星载雷达数据，构建时序干涉测量模型包括： 根据获取的时序星载雷达数据，确定时序干涉测量模型的基本观测参照值； 建立待估参数与所述基本观测参照值之间的联系，作为时序干涉测量模型。3. 根据权利要求2所述的时序星载雷达数据处理方法，其特征在于，所述根据获取的时 序星载雷达数据，确定时序干涉测量模型的基本观测参照值包括： 获取研究区域的所有时序星载雷达数据并进行干涉组合得到干涉对； 对所述干涉对进行采样和相干性计算，并统计相干性高于预设阈值的采样点和所述采 样点的数量； 根据所述采样点和所述采样点的数量，构建最小生成树； 对包含于所述最小生成树中的干涉对进行相干性评估，选取时域相干点； 对所述时域相干点根据空间分布和密度进行采样，得到采样相干点； 对所述采样相干点进行空间高密度构网，并将网的每条弧段的缠绕相位差作为时序干 涉测量模型的基本观测参照值。4. 根据权利要求3所述的时序星载雷达数据处理方法，其特征在于，所述建立待估参数 与所述基本观测参照值之间的联系包括： 所述待估参数包括形变序列、地形残差以及轨道误差，所述地形残差包括高程误差； 定义所述缠绕相位差与所述待估参数之间的观测方程为：其中，A为差分算子，N为N+1景所述雷达影像形成的N个干涉对，G为所述采样相干点构 建的所述弧段的个数，W为所述缠绕相位差中除了所述形变序列、所述地形残差以及所述轨 道误差外的噪声相位矢量；C为反映所述弧段和所述采样相干点之间关系 的邻接矩阵，B为系数阵； Ppar为包含了采样点上的所述形变序列和所述地形残差参数的参数阵；A表示获取的干涉图与所述雷达影像的构件关系，D表示所述弧段 上轨道多项式参数对应的系数矩阵；，W、匕、W为所述雷达影像的虚拟轨道误差参数，j表示弧段序号，相 邻弧段的所述雷达影像的虚拟轨道误差参数的代数差为所述干涉图的轨道误差参数。5. 根据权利要求4所述的时序星载雷达数据处理方法，其特征在于，所述基于所述轨道 误差和所述高程误差重构基线误差相位包括： 定义垂直基线相位重构公式为：其中，λ为雷达波长，P为卫星到给定地面点的距离，Θ为所述地面点所对应的雷达波束 入射角，B丄为所述地面点对应的垂直基线，h为将所述高程误差改正到外部高程数据上后得 到的所述地面点精化后的高程； 定义平行基线相位重构公式为：其中，λ为雷达波长，B=为所述地面点对应的平行基线，X和y分别是以所述雷达影像左 上角点为坐标原点的像素横坐标和像素纵坐标，a、b、c为所述干涉图轨道误差参数； 通过所述垂直基线相位重构公式和所述平行基线相位重构公式完成基线误差相位的 重构。6. -种时序星载雷达数据处理装置，其特征在于，包括： 模型构建模块，用于根据获取的时序星载雷达数据，构建时序干涉测量模型，所述时序 星载雷达数据包括雷达影像，所述时序干涉测量模型用于计算目标卫星的轨道误差和高程 误差； 误差计算模块，用于根据所述时序干涉测量模型计算所述轨道误差和所述高程误差； 重构模块，用于基于所述轨道误差和所述高程误差重构基线误差相位，所述基线误差 相位用于确定基线及基线的变化率。7. 根据权利要求6所述的时序星载雷达数据处理装置，其特征在于，所述模型构建模块 包括： 观测确定子模块，用于根据获取的时序星载雷达数据，确定时序干涉测量模型的基本 观测参照值，所述时序星载雷达数据包括雷达影像，所述时序干涉测量模型用于计算目标 卫星的轨道误差和高程误差； 联系建立子模块，用于建立待估参数与所述基本观测参照值之间的联系，作为时序干 涉测量模型。8. 根据权利要求7所述的时序星载雷达数据处理装置，其特征在于，所述观测确定子模 块还用于： 获取研究区域的所有时序星载雷达数据并进行干涉组合得到干涉对； 对所述干涉对进行采样和相干性计算，并统计相干性高于预设阈值的采样点和所述采 样点的数量； 根据所述采样点和所述采样点的数量，构建最小生成树； 对包含于所述最小生成树中的干涉对进行相干性评估，选取时域相干点； 对所述时域相干点根据空间分布和密度进行采样，得到采样相干点； 对所述采样相干点进行空间高密度构网，并将网的每条弧段的缠绕相位差作为时序干 涉测量模型的基本观测参照值。9. 根据权利要求8所述的时序星载雷达数据处理装置，其特征在于，所述联系建立子模 块还用于： 所述待估参数包括形变序列、地形残差以及轨道误差，所述地形残差包括高程误差； 定义所述缠绕相位差与所述待估参数之间的观测方程为：其中，A为差分算子，N为N+1景所述雷达影像形成的N个干涉对，G为所述采样相干点构 建的所述弧段的个数，W为所述缠绕相位差中除了所述形变序列、所述地形残差以及所述轨 道误差外的噪声相位矢量；C为反映所述弧段和所述采样相干点之间关系 的邻接矩阵，B为系数阵； Ppar为包含了采样点上的所述形变序列和所述地形残差参数的参数阵；A表示获取的干涉图与所述雷达影像的构件关系，D表示所述弧段 上轨道多项式参数对应的系数矩阵；W、匕、^为所述雷达影像的虚拟轨道误差参数，j表示弧段序号，相 邻弧段的所述雷达影像的虚拟轨道误差参数的代数差为所述干涉图的轨道误差参数。10. 根据权利要求9所述的时序星载雷达数据处理装置，其特征在于，所述重构模块还 用于： 定义垂直基线相位重构公式为：其中，λ为雷达波长，P为卫星到给定地面点的距离，Θ为所述地面点所对应的雷达波束 入射角，B丄为所述地面点对应的垂直基线，h为将所述高程误差改正到外部高程数据上后得 到的所述地面点精化后的高程； 定义平行基线相位重构公式为：其中，λ为雷达波长，B=为所述地面点对应的平行基线，X和y分别是以所述雷达影像左 上角点为坐标原点的像素横坐标和像素纵坐标，a、b、c为所述干涉图轨道误差参数； 通过所述垂直基线相位重构公式和所述平行基线相位重构公式完成基线误差相位的 重构。
【专利摘要】本发明涉及干涉测量技术领域，提供了一种时序星载雷达数据处理方法和装置，以解决现有技术中因相位误差干扰而导致的基线估计的精确度不高，以及在利用InSAR技术分析海量雷达干涉数据时需对每一个干涉图进行基线精化导致的效率低下的问题。所述方法包括：根据获取的时序星载雷达数据，构建时序干涉测量模型；根据时序干涉测量模型计算轨道误差和高程误差；基于轨道误差和高程误差重构基线误差相位，用于确定基线及基线的变化率。本发明的技术方案通过构建能够同时计算轨道误差和高程误差的时序干涉测量模型，实现从缠绕的雷达干涉相位中直接精确分离轨道误差相位和高程误差，并重构基线误差相位，从而提高基线估计的精确度和效率。
【IPC分类】G01S7/41
【公开号】CN105487065
【申请号】CN201610016363
【发明人】张磊, 丁晓利
【申请人】香港理工大学深圳研究院
【公开日】2016年4月13日
【申请日】2016年1月8日
【技术领域】
[0001] 本发明设及干设测量技术领域，尤其设及一种时序星载雷达数据处理方法和装 置。
【背景技术】
[0002] 自上世纪八十年代W来，基于星载雷达数据的干设测量技术已被广泛用于全天 候、全天时、大范围和高精度地监测地表及地面设施的形变，包括桥梁、大巧、机场及建筑综 合体等。目前，InSAR(Synthetic Aperture Radar Interferomehy，合成孔径雷达干设）测 量技术已经成为监测及评估地震、火山及滑坡等自然灾害的重要手段之一，在环境恶劣区 域，甚至是唯一手段。其干设技术原理可简单概括为通过比较在同一地区相同成像几何条 件下，不同时刻获取的复数雷达数据的相位信息来推算地表的形变。雷达卫星一般位于高 度为800千米左右的轨道上，由于当前定轨技术精度的限制，对同一地区卫星的重访位置难 W-致，从而造成在对地物目标两次成像时存在入射角度的差异，该差异可W用两次成像 时卫星之间的距离来描述，运个距离通常被称之为基线。基线长度的精确测定对利用干设 测量技术获取地面高程及形变具有至关重要的作用。直观上，如果知道两次卫星的坐标即 可算出基线长度，但由于定轨技术精度所限，所确定的卫星坐标通常会存在数厘米或数十 厘米的误差。运种轨道误差会传播到基线上，形成基线误差，最终体现到干设相位中，从而 成为InSAR技术的误差源之一。
[0003] 目前已有多种方法来估计基线误差。由于基线误差会通过平地效应和高程传播到 干设相位中，现有方法均通过对相位的分析来计算轨道误差。首先假设干设相位成功解缠， 即不存在相位模糊度，并且观测区的形变、地形残差和大气延迟误差可W忽略不计，运时解 缠相位的贡献主要与基线的平行雷达视线和垂直雷达视线运两个方向的分量有关，通过在 成像区域均匀选取高程已知的控制点，建立观测方程，利用最小二乘来解算影像中屯、点的 基线和基线变化率。但是运种方法需要干设相位的成功解缠，但事实上解缠的成功率在实 际中无法保证，特别是在低相干区域，例如植被区、地形复杂区及高楼林立的城区，从而影 响到基线估计的精确度，另外，该方法需要对每个干设基线进行改正，在利用InSAR技术分 析海量雷达干设数据时效率低下。
【发明内容】
[0004] 本发明的目的在于提供一种时序星载雷达数据处理方法和装置，旨在解决现有技 术中因相位误差，如解缠误差、长波形变信号、大气延迟误差W及地形残差等干扰而导致的 基线估计的精确度不高，W及在利用InSAR技术分析海量雷达干设数据时需对每一个干设 图进行基线精化导致的效率低下的问题。
[0005] 本发明第一方面，提供一种时序星载雷达数据处理方法，包括：
[0006] 根据获取的时序星载雷达数据，构建时序干设测量模型，所述时序星载雷达数据 包括雷达影像，所述时序干设测量模型用于计算目标卫星的轨道误差和高程误差；
[0007] 根据所述时序干设测量模型，计算所述轨道误差和所述高程误差；
[0008] 基于所述轨道误差和所述高程误差重构基线误差相位，所述基线误差相位用于确 定基线及基线的变化率。
[0009] 本发明第二方面，提供一种时序星载雷达数据处理装置，包括：。
[0010] 模型构建模块，用于根据获取的时序星载雷达数据，构建时序干设测量模型，所述 时序星载雷达数据包括雷达影像，所述时序干设测量模型用于计算目标卫星的轨道误差和 高程误差；
[0011] 误差计算模块，用于根据所述时序干设测量模型，计算所述轨道误差和所述高程 误差；
[0012] 重构模块，用于基于所述轨道误差和所述高程误差重构基线误差相位，所述基线 误差相位用于确定基线及基线的变化率。
[0013] 本发明与现有技术相比存在的有益效果是:本发明提供的技术方案根据时序星载 雷达数据，构建能够同时计算轨道误差和高程误差的时序干设测量模型，从而实现从缠绕 的雷达干设相位中直接精确分离出轨道误差和高程误差，并基于轨道误差和高程误差重构 基线误差相位，从而提高基线估计的精确度和效率。
【附图说明】
[0014] 图1是本发明实施例一提供的时序星载雷达数据处理方法的流程图；
[0015] 图2是本发明实施例一提供的时序星载雷达数据处理方法中轨道误差对基线估计 的影响的示意图；
[0016] 图3是本发明实施例一提供的时序星载雷达数据处理方法中雷达成像的几何示意 图；
[0017] 图4是本发明实施例一提供的时序星载雷达数据处理方法中典型的干设图中的基 线误差示意图；
[0018] 图5是本发明实施例二提供的时序星载雷达数据处理方法的流程图；
[0019] 图6是本发明实施例Ξ提供的时序星载雷达数据处理装置的组成示意图；
[0020] 图7是本发明实施例四提供的时序星载雷达数据处理装置的组成示意图。
【具体实施方式】
[0021] 为了使本发明的目的、技术方案及优点更加清楚明白，W下结合附图及实施例，对 本发明进行进一步详细说明。应当理解，此处所描述的具体实施例仅仅用W解释本发明，并 不用于限定本发明。
[0022] W下结合具体附图对本发明的实现进行详细的描述。
[0023] 实施例一：
[0024] 图1是本发明实施例一提供的时序星载雷达数据处理方法的流程图，具体包括步 骤S101至S103,详述如下：
[0025] S101、根据获取的时序星载雷达数据，构建时序干设测量模型，该时序星载雷达数 据包括雷达影像，该时序干设测量模型用于计算目标卫星的轨道误差和高程误差。
[0026] 载有合成孔径雷达的卫星可在800千米的轨道上对地球进行重复观测，地面站通 过星载GPS(Global化sitioning System,全球定位系统）等设备对卫星进行追踪及定位， 但由于定位精度的限制，星载GPS所确定的卫星位置不可避免存地在误差，两次卫星观测位 置间的距离即为基线。显然地，卫星的轨道误差会传播至基线上产生基线误差。图2示出了 轨道误差对基线估计的影响，其中S1和S2表示卫星对地观测时的位置，Bo为S1和S2之间的 矢量距离即基线，α为观测基线与水平方向夹角，β为误差基线与水平方向夹角，Θ为雷达波 入射角。真实的基线矢量Β等于Bo和由轨道误差引起的基线误差矢量η之和。
[0027] 具体地，根据研究区域内所获取的所有时序星载雷达数据，构建一个时序干设测 量模型，该时序干设测量模型能够将轨道误差和高程误差同时解算出来。
[0028] 考虑到海量的雷达数据，对于同一研究区域，每一期影像通常有数百万个有效监 测点，若数十期影像的数据同时处理，数据量将极为庞大。为了解决相位模糊度问题W及抑 制大气延迟误差，需要对有效监测点进行采样，并对采样点进行空间差分处理，将采样点的 差分相位作为基本观测参照值，可W采用顾及监测密度和相位梯度的四叉树采样来确定基 本观测参照值。基本观测参照值确定之后，需要建立待估参数与基本观测参照值之间的关 系，即时序干设测量模型，其中待估参数包括采样点的形变序列、地形残差和轨道误差，其 中地形残差包括了高程误差。为了提高计算效率，可W根据干设图及雷达影像的线性组合 关系，将干设图轨道误差转为雷达影像的虚拟轨道误差，即将干设图中的轨道误差分量分 解到每一个雷达影像上，通过获取雷达影像的虚拟轨道误差参数来重构干设图的轨道误 差，建立待估参数与基本观测参照值之间的关系，作为时序干设测量模型。
[0029] S102、根据时序干设测量模型，计算轨道误差和高程误差。
[0030] 具体地，步骤S101得到的时序干设测量模型可W是一个超大型线性稀疏系统，采 用共辆梯度法，通过迭代方式完成对模型中待估参数的解算，从而分离出轨道误差和高程 误差。
[0031] S103、基于轨道误差和高程误差重构基线误差相位，该基线误差相位用于确定基 线及基线的变化率。
[0032] 由于InSAR测量技术是一种利用相位测距的技术，是通过测量两次卫星到地面目 标的距离，即斜距的变化来推算目标高程或形变。图3示出了卫星雷达成像的几何示意图， 其中S1和S2表示卫星对地观测时的位置，B为S1和S2之间的距离即基线，P为地表目标即相 干点，P'为地表目标在楠球面的投影，Z为目标点的高程，Η为卫星高，0代表楠球原点，α为基 线方向角，Θ为目标点对应的雷达波入射角，θ〇为忽略目标点高程的雷达波入射角，Ρ为卫星 至地面目标的距离，阴影区域表示地形。可W看到地面目标Ρ的高程及形变都会产生斜距 差，同样地，基线Β的基线误差也会产生斜距差，运些斜距差反映到卫星雷达的干设相位上 则显示为干设条纹，如图4所示。在图4显示的干设图中可W明显看到因基线误差引起的干 设条纹，条纹越密说明基线误差越大，该示例数据来自于C波段的化Vi sat/ASAR星载雷达， 成像区为罗布泊。因此，基线误差相位可W有两部分组成，包括平行基线相位和高程作用下 的垂直基线相位。
[0033] 具体地，在重构基线误差相位时，根据时序干设测量模型计算得到的轨道误差和 高程误差，分别重构水平基线相位和垂直基线相位，从而完成基线误差相位的重构。对重构 后的基线误差相位，可W利用最小二乘原理来确定基线W及基线的变化率。
[0034] 本实施例中，根据时序星载雷达数据，构建能够同时计算轨道误差和高程误差的 时序干设测量模型，从而实现从缠绕的雷达干设相位中直接精确分离轨道误差和高程误 差，并基于轨道误差和高程误差重构基线误差相位，从而提高基线估计的精确度和效率。 [003引实施例二：
[0036] 图5 是本发明实施例二提供的时序星载雷达数据处理方法的流程图，具体包括步 骤S201至S204,详述如下：
[0037] S201、根据获取的时序星载雷达数据，确定时序干设测量模型的基本观测参照值， 该时序星载雷达数据包括雷达影像，该时序干设测量模型用于计算目标卫星的轨道误差和 高程误差。
[0038] 具体地，确定时序干设测量模型的基本观测参照值可W通过步骤S2011至S2016实 现，详细说明如下：
[0039] S2011、获取研究区域的所有时序星载雷达数据并进行干设组合得到干设对。
[0040] 具体地，将研究区域内所获取的所有时序星载雷达数据进行干设组合，得到干设 对，并计算干设对的时间基线和空间基线，根据基线信息对干设对进行初选。
[0041] 干设对表示两幅雷达影像的相位差，假设有N+1幅雷达影像，最多可形成的干设组 合数量为〇.5XNX(N+l)。时间基线为构成干设对的两幅雷达影像的时间差。空间基线为卫 星在获取构成干设对的两幅影像时的位置差异。时间基线和空间基线可通过雷达影像的参 数文件计算得到。
[0042] 在进行干设对初选时，将所有干设对按基线信息在二维空间进行排列，用具备等 步长增长的捜索窗从零开始捜索，当捜索窗内所包含的干设对包含所有影像并且干设图数 量最小时停止捜索。
[0043] S2012、对干设对进行采样和相干性计算，并统计相干性高于预设阔值的采样点和 采样点的数量。
[0044] 具体地，针对步骤S2011中初选得到的干设对，将包含于运些干设对中的雷达影像 进行基于灰度的四叉树采样，并对采样数据进行相干性计算，然后统计各干设对的相干性 高于预设阔值的采样点和采样点的数量。
[0045] 初选的干设对数量可能很大，为了提高计算效率，需要进一步筛选高质量的干设 对。在雷达干设测量技术中一般采用相干性来衡量干设对的质量，但逐个像素计算相干性 比较费时，因此可W先通过四叉树方法对雷达强度影像图根据灰度信息进行四叉树采样， 然后计算采样点的相干性，用有限采样点的相干性来代表整幅干设图的质量。
[0046] 基于灰度的四叉树采样方法首先把雷达强度影像图等分成四个分区，如果检查到 某个分区所有格网的灰度值之间的差值小于或等于预置数值，则该某个分区就不再往下分 害d;否则，把该某个分区再分割成四个子区域，运样递归地分割，直至每个分区所有格网的 灰度值之间的差值小于或等于预置数值为止。
[0047] 相干性是一个介于0和1之间的值，可由如下公式计算：
[004引
[0049]其中，Cl和C2为同一地面目标在两次卫星获取影像时所反射的雷达波信号，该信号 为复数值，包含了地面目标对应的强度和相位信息，(Λ'为C2的复数共辆，Φ<1β*为地形和形变 贡献的相位，为取模算子。
[0050] S2013、根据采样点和采样点的数量，构建最小生成树。
[0051] 具体地，根据步骤S2012得到的采样点的相干性信息和采样点的数量，构建最小生 成树(Minimum Spanning化ee，MST)。该生成树可W在N景雷达影像中筛选出质量最优的N-1个干设对，其中N表示雷达影像数量，时序干设测量模型的基本观测参照值将从运些干设 对中选取。
[0052] 生成树是连通图的极小连通子图。极小是指若在树中任意增加一条边，则将出现 一个回路，若去掉一条边，将会使之变成非连通图。生成树各边的权值总和称为生成树的 权。权最小的生成树称为最小生成树。
[0053] 进一步地，将步骤S2012得到的采样点的相干性均值的倒数作为权进行MST捜索， 按照预定的构建最小生成树的方法构建最小生成树。
[0054] S2014、对包含于最小生成树中的干设对进行相干性评估，选取时域相干点。
[0055] 具体地，为了减少相干性估计的偏差，在对包含于最小生成树中的干设对进行相 干性评估时，首先要对采样窗的像素点进行同质点选取，即选取与中屯、像素点散射特征的 偏差在预定范围内的周边像素点，再计算所选取的像素点的相干性，并根据所计算的相干 性选取相干性小于经验性阔值的像素点，作为时域相干点。
[0056] 假设采样窗内有Μ个同质像素点，则相干性可W由W下公式计算：
[0化7]
[0058] 其中，Cl和C2为同一地面目标在两次卫星获取影像时所反射的雷达波信号，该信号 为复数值，包含了地面目标对应的强度和相位信息，C为C2的复数共辆，Φ det为地形和形变 贡献的相位，1[1为取模算子。
[0059] 利用上述公式所计算的相干性，通过设定一个经验性阔值可将时域相干点筛选出 来，经验性阔值可W设置为0.1。
[0060] S2015、对时域相干点根据空间分布和密度进行采样，得到采样相干点。
[0061] 具体地，对步骤S2014选出的时域相干点根据空间分布和密度进行四叉树采样，得 到采样相干点。通常选取一万个采样相干点已经足够恢复出轨道相位，即当四叉树划分的 窗口数量达到一万时可W停止操作。
[0062] S2016、对采样相干点进行空间高密度构网，并将网的每条弧段的缠绕相位差作为 时序干设测量模型的基本观测参照值。
[0063] 具体地，对步骤S2015得到的采样相干点进行空间高密度构网，网的每条边，即每 条弧段的缠绕相位差作为时序干设测量模型的基本观测参照值。
[0064] 高密度构网是把采样点密集地连接起来。考虑到大气延迟的空间相关性，将研究 区域划成1km X 1 km的格网，W1 km X 1 km格网的顶点为中屯、、1 km为半径的范围内捜索采样 点，将捜索到的采样点构建为Delaunays角网，实现采样点的连接。
[0065] 缠绕相位差即采样点的相位差。采样点反射的雷达波的原始干设相位Φ永远介于 (-η，π]之间，总是与实际的相位差如的整数倍，所W称为缠绕相位。假设两个采样点上的缠 绕相位分别为Φ谢Φ2,则缠绕相位差
，其中mod( ·)为余数算子。
[0066] S202、建立待估参数与基本观测参照值之间的联系，作为时序干设测量模型。
[0067] 具体地，待估参数主要包括采样相干点的形变序列、地形残差W及轨道误差，其中 地形残差包括了高程误差，通过建立待估参数与基本观测参照值即缠绕相位差之间的观测 方程，体现待估参数与基本观测参照值之间的联系，作为时序干设测量模型。
[0068] 观测方程的具体建立过程如下：
[0069] -方面，轨道误差贡献的差分相位与轨道参数之间的关系为： 剛
（3)
[0071] 其中，Δ为差分算子，N为N+1景影像形成的N个干设对，G为采样相干点构建的弧段 的个数，Dorb为系数阵，Porb为参数阵；
[0072] 系数阵Dnrb的计算公式为： W73]
(4)
[0074]其中，A表示获取的干设图与雷达影像的构件关系，D表示弧段上轨道多项式参数 对应的系数矩阵；
[0077] 其中，dX、dY、dXY分别为弧段上轨道多项式参数对应的系数;j表示弧段序号；
[0078] 参数阵Pnrb的计算公式为：
[0079]
所
[0080] 其中，w、bj、Cj为雷达影像的虚拟轨道误差参数。
[0081] 另一方面，在差分相位中除了轨道误差的贡献，还有高程和地表形变的贡献。具体 地，在第i个干设图中，第P个采样点的地形残差相位与该点的高程误差A hp存在确定 性关系，可表示如下：
[0082]
巧）
[0083] 其中，β?为一个与基线、入射角及斜距有关的系数。
[0084] 类似地，可W建立地表形变序列贡献的相位祐6，,Ρ与形变Vp的关系，可表示如下：
[0085]
(9)
[0086] 其中，钱是一个与时间有关的系数。
[0087] 则根据公式(8)和公式(9)，可W得到地形残差及形变序列共同贡献的相位为：
[008引 (10)
[0089] 其中B为系数阵，N为干设图数量。
[0090] 根据公式（10)可进一步得到采样点间的相位差矢量与地形残差及形变序列的关 系为： 幽]
(11)
[0092] 其4
C为反映弧段和采样相干点之间关系 的邻接矩阵，Ppar为包含了采样点上的地形残差及形变序列参数的参数阵。
[0093] 综合公式(3)和公式（11)，缠绕相位差与轨道误差，地形残差及形变序列之间的观 测方程可W表示为：
[0094]
(12)
[0095] 其中，D=[D0rb Dpar]T，P=[P0rb Ppar]T，W为缠绕相位差中除了轨道误差、地形残差 及形变序列外的噪声相位矢量。
[0096] 该观测方程（12)作为时序干设测量模型，体现了待估参数与基本观测参照值之间 的联系。在进行轨道误差参数化时，为了减少参数个数，考虑到干设图和雷达影像的线性组 合关系，将干设图中的轨道误差分量分解到每个雷达影像上，一旦获得雷达影像的虚拟轨 道误差参数，即可重构任何一个干设图的轨道误差。干设图轨道误差可W是相邻弧段的雷 达影像虚拟轨道误差的代数差。
[0097] 具体地，通过雷达影像的虚拟轨道误差参数重构干设图的轨道误差的过程如下：
[0098] 将雷达影像1和雷达影像2表示为相邻弧段的雷达影像；
[0099] 雷达影像1的虚拟轨道误差多项式为：
[0100]
(巧) [0101] 其中，ai、bi、ci为雷达影像1的虚拟轨道误差参数，X和y分别是W雷达影像左上角 点为坐标原点的像素横坐标和像素纵坐标；
[0102] 雷达影像2的虚拟轨道误差多项式为：
[01。引
（14)
[0104] 其中，a2、b2、C2为雷达影像2的虚拟轨道误差参数，X和y分别是W雷达影像左上角 点为坐标原点的像素横坐标和像素纵坐标。
[0105] 则根据多项式(13)和多项式(14)，可W得到干设图的轨道误差多项式为：
[0106]
(15)
[0107] 其中，a、b、c为干设图轨道误差参数，X和y分别是W雷达影像左上角点为坐标原点 的像素横坐标和像素纵坐标。
[010引通过观测方程解算出多项式（13)的参数ai、bi、C1和多项式（14)的参数日2、b2、C2，可 W得到多项式(15)的参数a、b、c，进而得到干设图的轨道误差分量。
[0109 ] S203、根据时序干设测量模型，计算轨道误差和高程误差。
[0110] 具体地，步骤S202中构建的观测方程是一个超大型线性稀疏系统，可W采用共辆 梯度法，通过迭代方式实现对参数的解算，进而从干设相位中分离出轨道误差分量及高程 误差分量。
[0111] S204、基于轨道误差和高程误差重构基线误差相位，该基线误差相位用于确定基 线及基线的变化率。
[0112] 由于InSAR测量技术是一种利用相位测距的技术，是通过测量两次卫星到地面目 标的距离，即斜距的变化来推算目标高程或形变。图3示出了卫星雷达成像的几何示意图， 其中S1和S2表示卫星对地观测时的位置，B为S1和S2之间的距离即基线，P为地表目标即相 干点，P'为地表目标在楠球面的投影，Z为目标点的高程，Η为卫星高，0代表楠球原点，α为基 线方向角，Θ为目标点对应的雷达波入射角，θ〇为忽略目标点高程的雷达波入射角，Ρ为卫星 至地面目标的距离，阴影区域表示地形。可W看到地面目标Ρ的高程及形变都会产生斜距 差，同样地，基线Β的基线误差也会产生斜距差，运些斜距差反映到卫星雷达的干设相位上 则显示为干设条纹，如图4所示。在图4显示的干设图中可W明显看到因基线误差引起的干 设条纹，条纹越密说明基线误差越大，该示例数据来自于C波段的化Vi sat/ASAR星载雷达， 成像区为罗布泊。因此，基线误差相位可W有两部分组成，包括平行基线相位和高程作用下 的垂直基线相位。
[0113] 具体地，在重构基线误差相位时，首先将时序干设测量模型计算的高程误差改正 到外部高程数据上，W改善外部高程数据的精度，然后根据采样点的垂直基线信息重构垂 直基线相位。
[0114] 定义垂直基线相位重构公式为：
[0115]
(1食)
[0116] 其中，λ为雷达波长，P为卫星到给定地面点的距离，Θ为该给定地面点所对应的雷 达波束入射角，Β丄为该给定地面点对应的垂直基线，h为将高程误差改正到外部高程数据上 后得到的该给定地面点精化后的高程。
[0117] 另一方面，平行基线贡献的相位仅与基线有关，且在空线上任意一点可用一个多 项式代替，即平行基线相位重构公式可定义为：
[0118]
(17；
[0119] 其中，λ为雷达波长，B=为所述地面点对应的平行基线，X和y分别是W雷达影像左 上角点为坐标原点的像素横坐标和像素纵坐标，a、b、c为干设图轨道误差参数，可由虚拟轨 道误差多项式线性组合得到。
[0120] 通过垂直基线相位重构公式（16)和平行基线相位重构公式（17)即可完成基线误 差相位的重构。对重构后的基线误差相位，可W利用最小二乘原理来确定基线W及基线的 变化率。
[0121 ]本实施例中，首先对获取的时序星载雷达数据进行干设组合，将得到的干设对进 行采样和相干性计算后，对所构建的最小生成树进行相干性评估，选取时域相干点，并对时 域相干点进行空间分布和密度的四叉树采样得到采样相干点，对采样相干点进行空间高密 度构网得到时序干设测量模型的基本观测参照值，并建立待估参数和基本观测参照值之间 的联系，作为时序干设测量模型，通过解算该模型实现从缠绕的雷达干设相位中直接精确 分离出轨道误差和高程误差，并基于轨道误差和高程误差重构基线误差相位，从而提高基 线估计的精确度和效率。
[0122] 实施例
[0123] 图6是本发明实施例Ξ提供的时序星载雷达数据处理装置的组成示意图，为了便 于说明，仅示出了与本发明实施例相关的部分。图8示例的用于时序星载雷达数据处理的装 置可W是前述实施例一提供的用于时序星载雷达数据处理的方法的执行主体，其可W是应 用设备或者应用设备中的一个功能模块。图8示例的用于时序星载雷达数据处理的装置主 要包括:模型构建模块81、误差计算模块82、重构模块83。各功能模块详细说明如下：
[0124] 模型构建模块81，用于根据获取的时序星载雷达数据，构建时序干设测量模型，所 述时序星载雷达数据包括雷达影像，所述时序干设测量模型用于计算目标卫星的轨道误差 和高程误差；
[0125] 误差计算模块82，用于根据所述时序干设测量模型，计算所述轨道误差和所述高 程误差；
[0126] 重构模块83,用于基于所述轨道误差和所述高程误差重构基线误差相位，所述基 线误差相位用于确定基线及基线的变化率。
[0127] 本实施例中的时序星载雷达数据处理的装置中，各模块实现各自功能的过程，参 见前述图1所示实施例的描述，此处不再寶述。
[0128] 从上述图8示例的用于时序星载雷达数据处理的装置可知，本实施例中，根据时序 星载雷达数据，构建能够同时计算轨道误差和高程误差的时序干设测量模型，从而实现从 缠绕的雷达干设相位中直接精确分离轨道误差和高程误差，并基于轨道误差和高程误差重 构基线误差相位，从而提高基线估计的精确度和效率。
[0129] 实施例四：
[0130] 图7是本发明实施例四提供的时序星载雷达数据处理装置的组成示意图，为了便 于说明，仅示出了与本发明实施例相关的部分。图9示例的用于时序星载雷达数据处理的装 置可W是前述实施例二提供的用于时序星载雷达数据处理的方法的执行主体，其可W是应 用设备或者应用设备中的一个功能模块。图9示例的用于时序星载雷达数据处理的装置主 要包括:模型构建模块91、误差计算模块92、重构模块93。各功能模块详细说明如下：
[0131] 模型构建模块91，用于根据获取的时序星载雷达数据，构建时序干设测量模型，所 述时序星载雷达数据包括雷达影像，所述时序干设测量模型用于计算目标卫星的轨道误差 和高程误差；
[0132] 误差计算模块92,用于根据所述时序干设测量模型，计算所述轨道误差和所述高 程误差；
[0133] 重构模块93,用于基于所述轨道误差和所述高程误差重构基线误差相位，所述基 线误差相位用于确定基线及基线的变化率。
[0134] 具体地，模型构建模块91，还包括：
[0135] 观测确定子模块911，用于根据获取的时序星载雷达数据，确定时序干设测量模型 的基本观测参照值，所述时序星载雷达数据包括雷达影像，所述时序干设测量模型用于计 算目标卫星的轨道误差和高程误差；
[0136] 联系建立子模块912,用于建立待估参数与所述基本观测参照值之间的联系，作为 时序干设测量模型。
[0137] 具体地，观测确定子模块911还用于：
[0138] 获取研究区域的所有时序星载雷达数据并进行干设组合得到干设对；
[0139] 对所述干设对进行采样和相干性计算，并统计相干性高于预设阔值的采样点和所 述采样点的数量；
[0140] 根据所述采样点和所述采样点的数量，构建最小生成树；
[0141 ]对包含于所述最小生成树中的干设对进行相干性评估，选取时域相干点；
[0142] 对所述时域相干点根据空间分布和密度进行采样，得到采样相干点；
[0143] 对所述采样相干点进行空间高密度构网，并将网的每条弧段的缠绕相位差作为时 序干设测量模型的基本观测参照值。
[0144] 具体地，联系建立子模块912还用于：
[0145] 所述待估参数包括形变序列、地形残差W及轨道误差，所述地形残差包括高程误 差；
[0146] 定义所述缠绕相位差与所述待估参数之间的观测方程为：
[0150] 其中，Δ为差分算子，N为N+1景所述雷达影像形成的N个干设对，G为所述采样相干 点构建的所述弧段的个数，W为所述缠绕相位差中除了所述形变序列、所述地形残差W及所 述轨道误差外的噪声相位矢量；
[0151]
C为反映所述弧段和所述采样相干点之间 关系的邻接矩阵，B为系数阵；
[0152] Ppar为包含了采样点上的所述形变序列和所述地形残差参数的参数阵；
[0153]
\表示获取的干设图与所述雷达影像的构件关系，D表示所述 弧段上轨道多项式参数对应的系数矩阵；
[0154]
j、bj、cj为所述雷达影像的虚拟轨道误差参数，j表示弧段序 号，相邻弧段的所述雷达影像的虚拟轨道误差参数的代数差为所述干设图的轨道误差参 数。具体地，重构模块93还用于：
[01W] 定义垂直基线相位重构公式为：
[0156]
[0157] 其中，λ为雷达波长，p为卫星到给定地面点的距离，θ为所述地面点所对应的雷达 波束入射角，Β丄为所述地面点对应的垂直基线，h为将所述高程误差改正到外部高程数据上 后得到的所述地面点精化后的高程；< br>[0158] 定义平行基线相位重构公式为：
[0159]
[0160] 其中，λ为雷达波长，B=为所述地面点对应的平行基线，X和y分别是W所述雷达影 像左上角点为坐标原点的像素横坐标和像素纵坐标，a、b、c为所述干设图轨道误差参数；
[0161] 通过所述垂直基线相位重构公式和所述平行基线相位重构公式完成基线误差相 位的重构。
[0162] 本实施例中的时序星载雷达数据处理的装置中，各模块实现各自功能的过程，参 见前述图5所示实施例的描述，此处不再寶述。
[0163] 从上述图9示例的用于时序星载雷达数据处理的装置可知，本实施例中，首先对获 取的时序星载雷达数据进行干设组合，将得到的干设对进行采样和相干性计算后，对所构 建的最小生成树进行相干性评估，选取时域相干点，并对时域相干点进行空间分布和密度 的四叉树采样得到采样相干点，对采样相干点进行空间高密度构网得到时序干设测量模型 的基本观测参照值，并建立待估参数和基本观测参照值之间的联系，作为时序干设测量模 型，通过解算该模型实现从缠绕的雷达干设相位中直接精确分离出轨道误差和高程误差， 并基于轨道误差和高程误差重构基线误差相位，从而提高基线估计的精确度和效率。
[0164] 需要说明的是，本说明书中的各个实施例均采用递进的方式描述，每一个实施例 重点说明的都是与其他实施例的不同之处，各个实施例之间相同或者相似的部分互相参见 即可。对于装置类实施例而言，由于其与方法实施例基本相似，所W描述的比较简单，相关 之处参见方法实施例的部分说明即可。
[0165] 值得注意的是，上述装置实施例中，所包括的各个模块只是按照功能逻辑进行划 分的，但并不局限于上述的划分，只要能够实现相应的功能即可；另外，各功能模块的具体 名称也只是为了便于相互区分，并不用于限制本发明的保护范围。
[0166] 本领域普通技术人员可W理解，实现上述各实施例方法中的全部或部分步骤是可 W通过程序来指令相关的硬件来完成，相应的程序可W存储于一计算机可读取存储介质 中，所述的存储介质，如R0M/RAM、磁盘或光盘等。
[0167] W上所述仅为本发明的较佳实施例而已，并不用W限制本发明，凡在本发明的精 神和原则之内所作的任何修改、等同替换和改进等，均应包含在本发明的保护范围之内。
【主权项】
1. 一种时序星载雷达数据处理方法，其特征在于，包括： 根据获取的时序星载雷达数据，构建时序干涉测量模型，所述时序星载雷达数据包括 雷达影像，所述时序干涉测量模型用于计算目标卫星的轨道误差和高程误差； 根据所述时序干涉测量模型，计算所述轨道误差和所述高程误差； 基于所述轨道误差和所述高程误差重构基线误差相位，所述基线误差相位用于确定基 线及基线的变化率。2. 根据权利要求1所述的时序星载雷达数据处理方法，其特征在于，所述根据获取的时 序星载雷达数据，构建时序干涉测量模型包括： 根据获取的时序星载雷达数据，确定时序干涉测量模型的基本观测参照值； 建立待估参数与所述基本观测参照值之间的联系，作为时序干涉测量模型。3. 根据权利要求2所述的时序星载雷达数据处理方法，其特征在于，所述根据获取的时 序星载雷达数据，确定时序干涉测量模型的基本观测参照值包括： 获取研究区域的所有时序星载雷达数据并进行干涉组合得到干涉对； 对所述干涉对进行采样和相干性计算，并统计相干性高于预设阈值的采样点和所述采 样点的数量； 根据所述采样点和所述采样点的数量，构建最小生成树； 对包含于所述最小生成树中的干涉对进行相干性评估，选取时域相干点； 对所述时域相干点根据空间分布和密度进行采样，得到采样相干点； 对所述采样相干点进行空间高密度构网，并将网的每条弧段的缠绕相位差作为时序干 涉测量模型的基本观测参照值。4. 根据权利要求3所述的时序星载雷达数据处理方法，其特征在于，所述建立待估参数 与所述基本观测参照值之间的联系包括： 所述待估参数包括形变序列、地形残差以及轨道误差，所述地形残差包括高程误差； 定义所述缠绕相位差与所述待估参数之间的观测方程为：其中，A为差分算子，N为N+1景所述雷达影像形成的N个干涉对，G为所述采样相干点构 建的所述弧段的个数，W为所述缠绕相位差中除了所述形变序列、所述地形残差以及所述轨 道误差外的噪声相位矢量；C为反映所述弧段和所述采样相干点之间关系 的邻接矩阵，B为系数阵； Ppar为包含了采样点上的所述形变序列和所述地形残差参数的参数阵；A表示获取的干涉图与所述雷达影像的构件关系，D表示所述弧段 上轨道多项式参数对应的系数矩阵；，W、匕、W为所述雷达影像的虚拟轨道误差参数，j表示弧段序号，相 邻弧段的所述雷达影像的虚拟轨道误差参数的代数差为所述干涉图的轨道误差参数。5. 根据权利要求4所述的时序星载雷达数据处理方法，其特征在于，所述基于所述轨道 误差和所述高程误差重构基线误差相位包括： 定义垂直基线相位重构公式为：其中，λ为雷达波长，P为卫星到给定地面点的距离，Θ为所述地面点所对应的雷达波束 入射角，B丄为所述地面点对应的垂直基线，h为将所述高程误差改正到外部高程数据上后得 到的所述地面点精化后的高程； 定义平行基线相位重构公式为：其中，λ为雷达波长，B=为所述地面点对应的平行基线，X和y分别是以所述雷达影像左 上角点为坐标原点的像素横坐标和像素纵坐标，a、b、c为所述干涉图轨道误差参数； 通过所述垂直基线相位重构公式和所述平行基线相位重构公式完成基线误差相位的 重构。6. -种时序星载雷达数据处理装置，其特征在于，包括： 模型构建模块，用于根据获取的时序星载雷达数据，构建时序干涉测量模型，所述时序 星载雷达数据包括雷达影像，所述时序干涉测量模型用于计算目标卫星的轨道误差和高程 误差； 误差计算模块，用于根据所述时序干涉测量模型计算所述轨道误差和所述高程误差； 重构模块，用于基于所述轨道误差和所述高程误差重构基线误差相位，所述基线误差 相位用于确定基线及基线的变化率。7. 根据权利要求6所述的时序星载雷达数据处理装置，其特征在于，所述模型构建模块 包括： 观测确定子模块，用于根据获取的时序星载雷达数据，确定时序干涉测量模型的基本 观测参照值，所述时序星载雷达数据包括雷达影像，所述时序干涉测量模型用于计算目标 卫星的轨道误差和高程误差； 联系建立子模块，用于建立待估参数与所述基本观测参照值之间的联系，作为时序干 涉测量模型。8. 根据权利要求7所述的时序星载雷达数据处理装置，其特征在于，所述观测确定子模 块还用于： 获取研究区域的所有时序星载雷达数据并进行干涉组合得到干涉对； 对所述干涉对进行采样和相干性计算，并统计相干性高于预设阈值的采样点和所述采 样点的数量； 根据所述采样点和所述采样点的数量，构建最小生成树； 对包含于所述最小生成树中的干涉对进行相干性评估，选取时域相干点； 对所述时域相干点根据空间分布和密度进行采样，得到采样相干点； 对所述采样相干点进行空间高密度构网，并将网的每条弧段的缠绕相位差作为时序干 涉测量模型的基本观测参照值。9. 根据权利要求8所述的时序星载雷达数据处理装置，其特征在于，所述联系建立子模 块还用于： 所述待估参数包括形变序列、地形残差以及轨道误差，所述地形残差包括高程误差； 定义所述缠绕相位差与所述待估参数之间的观测方程为：其中，A为差分算子，N为N+1景所述雷达影像形成的N个干涉对，G为所述采样相干点构 建的所述弧段的个数，W为所述缠绕相位差中除了所述形变序列、所述地形残差以及所述轨 道误差外的噪声相位矢量；C为反映所述弧段和所述采样相干点之间关系 的邻接矩阵，B为系数阵； Ppar为包含了采样点上的所述形变序列和所述地形残差参数的参数阵；A表示获取的干涉图与所述雷达影像的构件关系，D表示所述弧段 上轨道多项式参数对应的系数矩阵；W、匕、^为所述雷达影像的虚拟轨道误差参数，j表示弧段序号，相 邻弧段的所述雷达影像的虚拟轨道误差参数的代数差为所述干涉图的轨道误差参数。10. 根据权利要求9所述的时序星载雷达数据处理装置，其特征在于，所述重构模块还 用于： 定义垂直基线相位重构公式为：其中，λ为雷达波长，P为卫星到给定地面点的距离，Θ为所述地面点所对应的雷达波束 入射角，B丄为所述地面点对应的垂直基线，h为将所述高程误差改正到外部高程数据上后得 到的所述地面点精化后的高程； 定义平行基线相位重构公式为：其中，λ为雷达波长，B=为所述地面点对应的平行基线，X和y分别是以所述雷达影像左 上角点为坐标原点的像素横坐标和像素纵坐标，a、b、c为所述干涉图轨道误差参数； 通过所述垂直基线相位重构公式和所述平行基线相位重构公式完成基线误差相位的 重构。
【专利摘要】本发明涉及干涉测量技术领域，提供了一种时序星载雷达数据处理方法和装置，以解决现有技术中因相位误差干扰而导致的基线估计的精确度不高，以及在利用InSAR技术分析海量雷达干涉数据时需对每一个干涉图进行基线精化导致的效率低下的问题。所述方法包括：根据获取的时序星载雷达数据，构建时序干涉测量模型；根据时序干涉测量模型计算轨道误差和高程误差；基于轨道误差和高程误差重构基线误差相位，用于确定基线及基线的变化率。本发明的技术方案通过构建能够同时计算轨道误差和高程误差的时序干涉测量模型，实现从缠绕的雷达干涉相位中直接精确分离轨道误差相位和高程误差，并重构基线误差相位，从而提高基线估计的精确度和效率。
【IPC分类】G01S7/41
【公开号】CN105487065
【申请号】CN201610016363
【发明人】张磊, 丁晓利
【申请人】香港理工大学深圳研究院
【公开日】2016年4月13日
【申请日】2016年1月8日
最新回复(0)