《星载INSAR系统的外定标方法.pdf》由会员分享，可在线阅读，更多相关《星载INSAR系统的外定标方法.pdf（14页珍藏版）》请在专利查询网上搜索。
1、(10)申请公布号 CN 103364766 A (43)申请公布日 2013.10.23 CN 103364766 A *CN103364766A* (21)申请号 201210096161.6 (22)申请日 2012.04.01 G01S 7/40(2006.01) (71)申请人 中国科学院电子学研究所 地址 100080 北京市海淀区北四环西路 19 号 (72)发明人 王宇 丁赤飚 梁兴东 (74)专利代理机构 中科专利商标代理有限责任 公司 11021 代理人 宋焰琴 (54) 发明名称 星载 InSAR 系统的外定标方法 (57) 摘要 本发明公开了一种星载 InSAR 系统的。
2、外定标 方法。该外定标方法包括 ： 由卫星高度计测得海 面高度数据， 将其去除时变影响因素后获得平均 海面高度数据 ； 由星载 InSAR 系统获取的预设时 间和空间分布的海洋图像对和未定标的干涉参 数， 得到星载 InSAR 系统测量的海面高度数据 ； 仿 真计算在预设时间和空间分布条件下除传播影响 外的时变影响因素带来的海面高度变化 ； 由平均 海面高度和除传播影响外的时变影响因素带来 的海面高度变化， 获取时变合成海面高度数据 ； 将时变合成海面高度数据作为标准数据， 对星载 InSAR 系统测量的海面高度数据进行校正， 实现 星载 InSAR 系统的定标。本发明采用基于海洋场 的外定标。
3、方法， 可对星载 InSAR 系统误差中的慢 变化成分进行长期动态监测。 (51)Int.Cl. 权利要求书 2 页 说明书 7 页 附图 4 页 (19)中华人民共和国国家知识产权局 (12)发明专利申请 权利要求书2页 说明书7页 附图4页 (10)申请公布号 CN 103364766 A CN 103364766 A *CN103364766A* 1/2 页 2 1. 一种星载 InSAR 系统的外定标方法， 包括 ： 由卫星高度计测得海面高度数据， 将其去除观测时刻和观测点经纬度对应的时变影响 因素后获得平均海面高度数据， 该时变影响因素包括 ： 传播影响、 海况影响、 潮汐影响及大 。
4、气压影响 ； 由星载 InSAR 系统获取的预设时间和空间分布的海洋图像对和未定标的干涉参数， 得 到星载 InSAR 系统测量的海面高度数据 ； 仿真计算在所述预设时间和空间分布条件下除传播影响外的时变影响因素带来的海 面高度变化 ； 由所述平均海面高度和所述除传播影响外的时变影响因素带来的海面高度变化， 获取 时变合成海面高度数据 ； 及 将所述时变合成海面高度数据作为标准数据， 对所述星载 InSAR 系统测量的海面高度 数据进行校正， 实现星载 InSAR 系统的定标。 2.根据权利要求1所述的星载InSAR系统的外定标方法， 其中， 所述将由卫星高度计测 得海面高度数据去除观测时刻和。
5、观测点经纬度对应的时变影响因素后获得平均海面高度 数据的步骤包括 ： 获取由海洋卫星高度计测得的预设时间内的数据记录 ； 根据上述数据记录中的每个点数据记录的校正项的标志位状态判别其是否为无效数 据， 从数据记录中剔除无效点数据记录 ； 从有效点数据记录中去除观测时刻和观测点经纬度对应的时变影响因素对海面高度 数据的影响， 计算有效点数据记录对应的海面高度数据 ； 将预设时间内的去除观测时刻和观测点经纬度对应的时变影响因素的海面高度数据 进行平均， 获取平均海面高度数据。 3.根据权利要求2所述的星载InSAR系统的外定标方法， 其中， 所述将预设时间内的去 除观测时刻和观测点经纬度对应的时变。
6、影响因素影响的海面高度数据进行平均， 获取平均 海面高度的步骤之后还包括 ： 将所述由单颗卫星高度计获取的平均海面高度数据作为初始海面高度数据 ； 将由多个海洋卫星高度计获得的所述初始平均海面高度数据再次平均， 获取平均海面 高度数据。 4. 根据权利要求 2 所述的星载 InSAR 系统的外定标方法， 其中， 所述卫星高度计为 ： 美 国 Jason-1 卫星高度计、 美国 Jason-2 卫星高度计或法国研究机构 CNES 的 CLS2011 模型中 各卫星高度计。 5.根据权利要求1所述的星载InSAR系统的外定标方法， 其中， 所述将时变合成海面高 度数据作为标准数据， 对星载 InS。
7、AR 系统测量的海面高度数据进行校正的步骤包括 ： 步骤 S110a， 由时变合成海面高度数据和星载 InSAR 系统测量的海面高度数据相比较， 得到高程误差 ； 步骤 S110b， 判断该高程误差是否满足预设的精度， 如果是， 则采用该高程误差对星载 InSAR 系统进行定标， 流程结束 ； 否则， 执行下一步骤 ； 步骤 S110c， 通过高程误差， 进行基于敏感度方程的干涉参数误差解算， 得到干涉参数 误差向量， 其中， 敏感度方程中的参数为基线长度、 基线倾角和干涉相位 ； 权 利 要 求 书 CN 103364766 A 2 2/2 页 3 步骤 S110d， 利用该干涉参数误差向量。
8、对上一次的干涉参数进行修正 ； 步骤 S110e， 由星载 InSAR 系统获取的预设时间和空间分布的海洋图像对和修正后的 干涉参数进行地形三维重建， 得到星载 InSAR 系统观测的海面高度， 执行步骤 S110a。 6.根据权利要求5所述的星载InSAR系统的外定标方法， 其中， 所述步骤S110c中敏感 度方程为 ： 其中， Pz为 L1 高程数据误差向量 ； F 为 L3 敏感度矩阵 ； X 为 31 待估计的干 涉参数误差向量， l 表示基线长度误差， 表示基线倾角误差， 表示干涉相位误差 ； h 为高程值 ； L 为定标数据点的数目。 7.根据权利要求1所述的星载InSAR系统的外。
9、定标方法， 其中， 所述由平均海面高度数 据和所述除传播影响外的时变影响因素带来的海面高度变化， 获取时变合成海面高度数据 的步骤是通过以下公式实现的 ： 其中 ： SSSH(T， lat， fong)为时变合成海面高度数据， MSSH(lat， long)为平均海面高度数据， 分别表示潮汐、 海况、 逆大气压时变影响因素带来的海面高度变 化。 8.根据权利要求1至7中任一项所述的星载InSAR系统的外定标方法， 其中， 所述时变 影响因素对应的模型为 ： 海况影响采用 Chambers 经验参数模型 ； 大气压影响中， 逆大气压影响采用 ECMWF 气象预报模型 ； 潮汐影响中， 洋潮、 负。
10、荷潮、 地球潮、 极潮影响均采用 GOT00.2 全球潮汐模型 ； 或 传播影响中， 干对流层和湿对流层影响采用 ECMWF 气象预报模型 ； 电离影像采用双频 电离校正模型。 权 利 要 求 书 CN 103364766 A 3 1/7 页 4 星载 InSAR 系统的外定标方法 技术领域 0001 本发明涉及电子行业雷达技术领域， 尤其涉及一种星载干涉合成孔径雷达 (Interferometric Synthetic Aperture Radar， 简称 InSAR) 系统的外定标方法。 背景技术 0002 星载 InSAR 系统是利用合成孔径雷达的相位信息提取地表的三维信息和高程变 化信。
11、息的一种技术， 在卫星或航天飞机上通过两部天线同时观测或两次单部天线的两次观 测获得地面场景的高程信息。 随着数字高程模型(Digital Elevation Model， 简称DEM)测 量精度的要求日益提高， 星载 InSAR 系统的定标方法的研究也愈发重要。 0003 InSAR 中， 影响 DEM 测量精度的主要因素包括 ： 基线长度和基线倾角测量误差、 系 统延时误差、 干涉相位偏置等， 这些误差最终在 InSAR 干涉处理中引入 DEM 误差， 上述误差 按特性可分为固定、 慢变和随机三种。因此， 这些参数误差的定标是有必要的。 0004 外定标目的就是通过目标的已知位置和高程信息。
12、， 标定基线长度、 基线倾角， 干涉 相位偏执， 系统延时等干涉参数的误差， 并去除这些参数误差， 利用准确的干涉参数进行高 程反演， 最终提高高程测量精度。 0005 传统的 InSAR 系统外定标方法利用人工点目标定标场 ( 地面布设角反射器 ) 或自 然分布目标定标场 ( 如热带雨林或其他已知 DEM 的区域等 ) 来对进行定标。利用人工定标 场或自然定标场的已知位置和高程信息， 求解 InSAR 系统的基线长度、 基线倾角、 干涉相位 偏置、 系统延时等参数。 0006 人工点目标定标场在地形复杂地区布设困难， 成本高， 不能经常对定标场成像， 只 能静态标定系统的固定误差， 难以准确。
13、反映误差的变化特性 ； 自然分布目标定标场同样存 在不能经常重访定标场， 不能长期连续的监测星载 InSAR 系统， 造成采样不足， 不能反映 InSAR 系统误差的慢变部分。 发明内容 0007 ( 一 ) 要解决的技术问题 0008 为解决上述的一个或多个问题， 本发明提供了一种星载 InSAR 系统的外定标方 法， 以解决了传统外定标方法难以长期动态监测的缺点。 0009 ( 二 ) 技术方案 0010 根据本发明的一个方面， 提供而来一种星载 InSAR 系统的外定标方法， 包括 ： 由卫 星高度计测得海面高度数据， 将其去除观测时刻和观测点经纬度对应的时变影响因素后获 得平均海面高度。
14、数据， 该时变影响因素包括 ： 传播影响、 海况影响、 潮汐影响及大气压影响 ； 由星载 InSAR 系统获取的预设时间和空间分布的海洋图像对和未定标的干涉参数， 得到星 载 InSAR 系统测量的海面高度数据 ； 仿真计算在预设时间和空间分布条件下除传播影响外 的时变影响因素带来的海面高度变化 ； 由平均海面高度和除传播影响外的时变影响因素带 来的海面高度变化， 获取时变合成海面高度数据 ； 及将时变合成海面高度数据作为标准数 说 明 书 CN 103364766 A 4 2/7 页 5 据， 对星载 InSAR 系统测量的海面高度数据进行校正， 实现星载 InSAR 系统的定标。 0011。
15、 ( 三 ) 有益效果 0012 从上述技术方案可以看出， 相比与传统的外定标方法， 本发明星载 InSAR 系统的 定标方法具有以下有益效果 ： 0013 (1) 本发明中， 采用基于海洋场的外定标方法， 可对星载 InSAR 系统误差中的慢变 化成分进行长期动态监测， 解决了传统外定标方法难以长期动态监测的缺点 ； 0014 (2) 本发明中， 利用卫星高度计数据， 不需要布设人工定标场， 避免了传统外定标 方法设计实施困难和采样不足的缺点， 从而可以精确获取星载 InSAR 系统误差特性， 有效 提高了星载 InSAR 系统的高程测量精度和可靠性。 附图说明 0015 图 1 为本发明实。
16、施例星载 InSAR 系统定标方法的原理框图 ； 0016 图 2 为本发明实施例星载 InSAR 系统定标方法中采用 Jason-1 数据计算的 MSSH 图 ； 0017 图 3 为本发明实施例星载 InSAR 系统定标方法中采用 CLS2011 模型的 MSSH 图 ； 0018 图4为本发明实施例星载InSAR系统定标方法中洋潮影响因素的高度计算结果的 仿真图 ； 0019 图5为本发明实施例星载InSAR系统定标方法中负荷潮影响因素的高度计算结果 的仿真图 ； 0020 图6为本发明实施例星载InSAR系统定标方法中地球潮影响因素的高度计算结果 的仿真图 ； 0021 图7为本发明实。
17、施例星载InSAR系统定标方法中极潮影响因素的高度计算结果的 仿真图 ； 0022 图8为本发明实施例星载InSAR系统定标方法中海况影响因素的高度计算结果的 仿真图 ； 0023 图9为本发明实施例星载InSAR系统定标方法中逆大气压影响因素的高度计算结 果的仿真图 ； 0024 图 10 为本发明实施例星载 InSAR 系统定标方法中干涉定标处理流程图。 具体实施方式 0025 为使本发明的目的、 技术方案和优点更加清楚明白， 以下结合具体实施例， 并参照 附图， 对本发明进一步详细说明。 0026 需要说明的是， 在附图或说明书描述中， 相似或相同的部分都使用相同的图号。 附 图中未绘示。
18、或描述的步骤或实现方式， 为所属技术领域中普通技术人员所知的形式。 另外， 虽然本文可提供包含特定值的参数的示范， 但应了解， 参数无需确切等于相应的值， 而是可 在可接受的误差容限或设计约束内近似于相应的值。 0027 本发明提出一种利用覆盖全球超过三分之二的海平面作为定标场， 借助海洋卫星 雷达高度计数据源可以实现星载 InSAR 系统基线长度、 基线倾角、 干涉相位偏执、 系统延时 等的标定， 从而可以对SAR系统设备进行长期连续的动态监测， 有助于提高星载InSAR系统 说 明 书 CN 103364766 A 5 3/7 页 6 的高程测量精度。 0028 在本发明的一个示例性实施例。
19、中， 提出了一种星载InSAR系统外定标方法。 图1为 本发明实施例星载 InSAR 系统外定标方法原理框图。如图 1 所示， 本实施例星载 InSAR 系 统外定标方法包括 ： 0029 步骤 S102， 由卫星高度计测得海面高度数据， 将其去除观测时刻、 观测点的经纬度 对应的时变影响因素后获得平均海面高度(Mean Sea Surface Height， 简称MSSH)数据， 该 时变影响因素包括 ： 传播影响、 海况影响、 潮汐影响及大气压影响 ； 0030 其中， 平均海面高度数据指海面相对于地球参考椭球的一定时间内的平均高度， 用来去除海面高度的时变影响因素以及无效海面高度数据的模。
20、型， 平均海面高度模型优选 地利用了美国 Jason-1 卫星高度计数据或法国研究机构 CNES 的 CLS2011 模型的数据。对 于定标精度要求不高的情况可采用单颗卫星高度计数据短时间平均获得 MSSH 模型。单颗 卫星高度计的 MSSH 模型的计算可细分三个子步骤 ： 0031 S102a， 获取由海洋卫星高度计测得的一定时间内的数据记录 ； 0032 以 Jason-1 为例， 其轨道周期 9.915 天， 每个周期有 127 条轨道， 其数据产品一般 按半条轨道 ( 称为 pass) 文件存放， 每个轨道周期 254 个 pass， 每个 pass 中包含 1 个文件 头和N个点数据。
21、记录， 每个点数据记录是由440个字节高度计测量信息构成， 包括该点的经 度、 纬度、 测量时间、 轨道高度、 高度计测距数据、 各种校正项、 各种标志位等信息。 0033 S102b， 根据上述数据记录中的每个点数据记录的轨道、 高度计和校正项的标志位 状态判别是否为无效数据， 并剔除无效数据记录 ； 0034 海面高度数据 (Sea Surface Height SSH) 可从子步骤 S102a 所述的每个 440 字 节的点数据记录的对应字节提取， 并通过公式 (1) 得到。 0035 0036 其中， t， lat， long 代表点数据记录的观测时刻和经纬度为卫星轨道高 度， 即轨道。
22、相对于地球参考椭球的高度 ；为高度计测距数据， 即卫星到海面的距离， 是对时变影响因素误差的校正。卫星轨道测量装置和高度计仪器设备的故障都会 带来和的无效， 从而导致 SSH 的无效，校正的无效也会导致 SSH 无 效。因此， 计算实施时， 首先需要根据每个点数据记录的轨道、 高度计和校正项的标志位状 态判别是否为无效数据， 并剔除无效数据。 0037 S102c， 从有效点数据记录中去除观测时刻和观测点经纬度对应的时变影响因素， 即公式 (1) 中的对海面高度数据的影响， 并计算有效点数据记录对应的海面高度 数据 ； 0038 卫星高度计直接观测的海面高度受多种时变影响因素影响， 计算 MS。
23、SH 之前还需 要对这些因素按观测时刻和观测位置进行校正。校正过程可表示为 ： 0039 (2) 0040 0041 其中， 上标 t， lat， long 代表观测时刻和观测点的经纬度， dry， wet， iono 分别代 表干、 湿对流层和电离层影响， 并统称为传播影响 ； ssb 代表海况影响 ； otide， ltide， etide 和 ptide 分别代表洋潮， 负荷潮， 固体地球潮和极潮影响， 并统称为潮汐影响 ； ibc 代表逆大 气压影响， 各影响因素的具体模型在步骤 S104 中描述。 说 明 书 CN 103364766 A 6 4/7 页 7 0042 至此， 可通过。
24、公式 (1) 计算每个有效点数据记录的 SSH。 0043 S102d， 将一定时间内的海面高度数据进行平均， 获取平均海面高度。 0044 MSSH 是海面相对于地球参考椭球的一定时间内的平均高度， 即对子步骤 iii) 中 的 SSH 进行时间平均， 可表示为 ： 0045 0046 附图 2 利用美国 Jason-1 卫星的数个完整轨道周期 ( 轨道周期约为 10 天 ) 进行 计算给出了利用第 340-353( 对应 2011 年 3 月 30 日至 8 月 16 日 ) 轨道计算的 MSSH。 0047 对于精度要求高的应用， 可以做更长时间跨度的平均计算， 以去初长周期变化的 影响。
25、因素， 计算步骤类似于单颗卫星高度计步骤。还可以采用多颗卫星高度计长时间平均 获得 MSSH 数据， 例如采用法国 CNES 的 CLS2011 模型的结果， 该模型利用 TOPEX/POSEIDON， ERS-2， GFO， JASON-1， ENVISAT 多个卫星高度计共计 15 年的数据， 得到的 1/30 度间距的精 细 MSSH， 其计算结果如附图 3 所示。 0048 步骤 S104， 由星载 InSAR 系统获取的预设时间和空间分布的海洋图像对和未定标 的干涉参数进行地形三维重建， 得到星载 InSAR 系统测量的海面高度数据 ； 0049 本步骤中， 由海洋图像获取海面高度的。
26、技术属于现有技术， 此处不再详细描述。 该 步骤指由星载 InSAR 系统海洋图像对和未定标的干涉参数， 进行地形三维重建， 地形三维 重建中包括干涉处理即由海洋图像对的复数据共轭相乘得到干涉相位， 干涉相位的展开， 在得到展开相位后， 利用未定标的干涉参数进行三维重建， 得到星载 InSAR 系统测量的海 面高度 (DEM)， 也就是获得了海洋上图像点的平面位置和高程信息。 0050 步骤 S106， 仿真计算除传播影响之外的时变影响因素在上述时间和空间条件下带 来的海面高度变化 ； 0051 在步骤 S102 的子步骤 S102c 中已提到海面高度的时变影响因素主要分为 4 类， 即 传播。
27、影响， 海况影响， 潮汐影响和大气压影响。这些时变影响因素的造成高度计测量误差， 并且使不同观测时刻的同一区域的海面高度不同， 其计算分别描述如下 ： 0052 传播影响也称为地球物理影响或环境影响， 主要指卫星高度计的微波测量信号传 播速度受到的地球大气折射的影响， 从而影响高度测量， 包括干、 湿对流层和电离层影响， 干对流层影响可以通过公式 (4) 计算。 0053 0054 其中，为以毫巴为单位的海面压力， t， lat， long 分别代表所处的观测时刻 和位置，可以通过欧洲 ECMWF 气象预报模型获得。湿对流层采用类似方法借助欧洲 ECMWF 气象预报模型计算 ； 电离层影响通过。
28、 Jason-1 卫星高度计的双频电离校正模型计算。 0055 潮汐是直接观测海面高度的重要组成部分， 比其他时变海洋信号的影响更大， 并 且可以预测。潮汐由四主要部分组成 ： 洋潮， 负荷潮， 地球潮， 极潮。前三种与地球的月日作 用力有关， 极潮是由于地球转动带来的与月日作用力无关。潮汐影响采用 GOT00.2 全球潮 汐模型， 模型利用球谐波的方法求出各潮汐分量， 如公式 (5) ： 0056 说 明 书 CN 103364766 A 7 5/7 页 8 0057 其中 h(t， lat， long) 是 t 时刻相对于参考椭球上的点 (lat， long) 的海面高度， 整数 N 代表。
29、潮汐成分的总数， k是潮汐成分 k 的角速度， Ak(lat， long) 和分别 是潮汐成分 k 的幅度和相位。对各类潮汐因素仿真计算的结果如图 4- 图 7 所示。 0058 海况影响指卫星高度计观测海面区域的海浪变化的影响， 对 Jason-1 卫星采用 Chambers 经验参数模型计算， 仿真计算的结果如图 8 所示。 0059 本实施例中， 大气压影响主要是逆大气压影响， 对其仿真计算的结果如图 9 所示。 0060 综上所述， 各时变影响因素及量级如表 1 所示。 0061 表 1 各时变影响因素模型 0062 0063 步骤 S108， 由平均海面高度 (MSSH) 和除传播影。
30、响之外的时变影响因素带来的海 面高度变化， 得到时变合成海面高度数据 (Synthetic Sea surface Height， 简称 SSSH) ； 0064 由于星载 InSAR 系统观测区域在时间和空间上与卫星高度计往往不重合， 为了利 用卫星高度计的数据进行星载 InSAR 系统的定标， 需要利用卫星高度计数据合成出与星载 InSAR 系统观测时间和空间一致的合成海面高度数据， 合成海面高度数据由两部分组成， 一 部分是静态的平均海面高度 (MSSH)， 一部分是时变的海面高度变化。合成时需要恢复的时 变部分主要指潮汐， 海况和逆大气压影响， 可用公式 (6) 表示。 0065 (6。
31、) 0066 0067 步骤 S110， 将时变合成海面高度 (SSSH) 数据作为标准数据， 对星载 InSAR 系统测 量的海面高度数据进行校正， 实现星载 InSAR 系统的定标。 0068 本步骤是指干涉定标利用步骤 S106 得到的 SSSH 作为已知位置和高程的定标数 据， 并采用基于敏感度方程的定标方法， 对星载 InSAR 系统测量的海面高度进行标定。 0069 干涉定标处理如图 10 所示， 通过卫星高度计间接的 SSSH 作为星载 InSAR 系统的 说 明 书 CN 103364766 A 8 6/7 页 9 定标数据 ( 参考高程数据 )， 星载 InSAR 系统三维重。
32、建模型、 敏感度方程构成了星载 InSAR 系统定标处理的基本组成部分。该步骤又可以分为以下子步骤 ： 0070 步骤S110a， 由SSSH和星载InSAR系统测量的海面高度(DEM)相比较， 得到高程误 差， 也就是比较相同空间位置的对应的 SSSH 和星载 InSAR 系统测量的海面高度， 得到高程 误差 ； 0071 步骤 S110b， 判断该高程误差是否满足预设的精度， 如果是， 执行步骤 S110c， 否 则， 执行步骤 S110c ； 0072 步骤 S110c， 星载 InSAR 系统的定标过程结束， 采用该高程误差对星载 InSAR 系 统进行定标 ； 0073 步骤 S11。
33、0c， 通过高程误差和敏感度方程， 进行基于敏感度方程的干涉参数误差解 算， 得到干涉参数误差向量 ； 0074 该步骤中， 当 S110a 得到的高程误差不满足星载 InSAR 系统高程精度测量要求时 执行， 当精度满足要求则定标结束或不需要定标， 一般采用未定标的干涉参数获得的高程 误差是不满足要求的。 0075 星载 InSAR 获得高程的处理中， 两部天线的物理间隔长度为基线长度， 基线相对 于大地水平面的夹角为基线倾角， 两部天线获得的干涉图像的复数据共轭相乘的相位差为 干涉相位， 这三个参数是星载 InSAR 获得高程不可缺少的参数， 也就是需要定标的参数。 0076 通过对基线长。
34、度 l， 基线倾角 ， 干涉相位 参数求偏导数， 可以得到下面的敏感 度方程， 干涉定标处理也就是求解敏感度方程的过程 ： 0077 0078 公式中 Pz为 L1 高程数据误差向量 ； F 为 L3 敏感度矩阵 ； X 为 31 待估 计的干涉参数误差向量， l 表示基线长度误差， 表示基线倾角误差， 表示干涉相位 误差 ； h 为高程值 ； L 为定标数据点的数目。 0079 这里定标的干涉参数主要为 3 个即基线长度， 基线倾角， 干涉相位， 理论上在 InSAR 测绘带的距离向有 3 个已知的参考高程就可以解算， 为提高精度可取 6 个， 例如对于 美国航天飞机 X 波段 InSAR，。
35、 测绘带宽 50km， 以 10km 里为间距， 沿距离向通过 SSSH 模型合 成出 6 个 SSSH 高度数据视为准确的参考高程数据向量， 并和 InSAR 观测得到的 DEM 高程求 差得到高程数据误差向量， 带入敏感度方程， 求解关于3个未知数有6个方程组成的超定方 程组， 得到干涉参数误差向量 ； 0080 步骤 S110d， 利用该干涉参数误差向量对上一次的干涉参数进行修正 ； 0081 步骤 S110e， 重新由星载 InSAR 系统海洋图像对和修正后的干涉参数进行地形三 维重建， 得到星载 InSAR 系统观测的海面高度， 执行步骤 S110a ； 0082 重复步骤 S110。
36、a 至 S110c， 直至高程精度达到定标要求的精度等级， 完成星载 InSAR 系统的干涉定标处理。 0083 需要说明的是， 上述对各步骤的定义并不仅限于实施方式中提到的各种具体结构 说 明 书 CN 103364766 A 9 7/7 页 10 或形状， 本领域的普通技术人员可对其进行简单地熟知地替换， 例如 ： 0084 (1)MSSH 的计算还可以利用美国 Jason-2 卫星高度计数据计算， 甚至利用我国神 州四号或海洋 2 号卫星数据形式计算 ； 所选用时变影响因素也可以是其他模型， 0085 (2) 干涉定标处理中的敏感度方程中不仅限于基线长度、 基线倾角和干涉相位三 个参数，。
37、 还可以引入更多参数如卫星轨道高度、 卫星姿态参数等， 从而实现对更多干涉参数 的计算和求解。 0086 以上所述的具体实施例， 对本发明的目的、 技术方案和有益效果进行了进一步详 细说明， 所应理解的是， 以上所述仅为本发明的具体实施例而已， 并不用于限制本发明， 凡 在本发明的精神和原则之内， 所做的任何修改、 等同替换、 改进等， 均应包含在本发明的保 护范围之内。 说 明 书 CN 103364766 A 10 1/4 页 11 图 1 图 2 说 明 书 附 图 CN 103364766 A 11 2/4 页 12 图 3 图 4 图 5 说 明 书 附 图 CN 103364766 A 12 3/4 页 13 图 6 图 7 图 8 图 9 说 明 书 附 图 CN 103364766 A 13 4/4 页 14 图 10 说 明 书 附 图 CN 103364766 A 14 。