］提出的改进相位曲率自聚焦方法（Ｗｅｉｇｈｔｅｄ　Ｐｈａｓｅ　Ｃｕｒｖａｔｕｒｅ　Ａｕｔｏｆｏｃｕｓ，ＷＰＣＡ），及钟雪莲［１２］提出的 ＭＴＰＴ（Ｍｕｌｔｉ－ｓｑｕｉｎｔ　Ｔｅｃｈｎｉｑｕｅｗｉｔｈ　Ｐｏｉｎｔ　Ｔａｒｇｅｔｓ）方法等。但前一类通过子孔径干涉相位估计 ＲＭＥ的方法不适用于非线性程度较高的 ＲＭＥ。后一类方法算法需要ＳＡＲ 影像包含高信噪比散射目标，大大缩 小了 该方法的适用范围。随着ＳＡＲ 系统的发展，多基线 ＳＡＲ 系统由于能提供多余观测，丰富成像信息 等优点，已成为主流发展趋势。但是，由于不同影像包含不同的 ＲＭＥ，在利用多基线ＳＡＲ 数据协同解算之前，必须对相位进行定标处理。由于上述两类方法只是用与单基线且存在不稳定性，难以满足多基线ＳＡＲ 数据对高精度相位定标的需求。为解决该问题，Ｂｉａｎｃｏ［１７］提出了一种基于永久散射体的多基线相位定标方法，该方法能有效校正多基线数据集的 ＲＭＥ，但由于在多基线定标的过程中均以主影像作为参照，因此定标后的 多 基 线 数 据 集 中 均 会 残 留 主 影 像 中 存 在的 ＲＭＥ。为解决上述问题，本文提出多基线干涉ＳＡＲ 联合定标两部法：首先，利 用基于永久散射体的多基线相位定标方法将数据集参照主影像进行相位校正。之后，引入多项式模型，从多基线相位定标后的干涉相位中估计出残留的主影像 ＲＭＥ。新方法被应用于 Ｅ－ＳＡＲ 系 统 的 Ｐ 波 段 机 载 重 轨 时 序 数据集。１　基于 ＰＳ的多基线相位定标两步法１．１　基于 ＰＳ技术的多基线相位定标方法传统 的 多 基 线 的 相 位 定 标 依 赖 于 地 面 控 制点［２０］，通过地面角反射器的高程信息，反解出相位误差，插值到整景影像后再进行校正。Ｂｉａｎｃｏ提出基于 ＰＳ多基线相位定标方法，利用 ＰＳ点代替地面的角反 射 器 作 为 控 制 点 进 行 数 据 集 的 相 位 定 标。ＰＳ点指的是在时序数据集中，相位幅度均保持相对稳定的像元点［１７］。它要求像元点在整个多基线数据集中包含的相位误差较小，通过相位误差较小的ＰＳ点，就能较为准确的估计出该点对应的地面 高程，从而较为准确的估计出其对应的相位误差［１６］。该方法首先要选出多基线数据集中的 ＰＳ候选点，可先计算配准后多基线数据集的每个像元对应的振幅离差因子 Ｄ（ｘ，ｒ）［１７］：Ｄ（ｘ，ｒ）＝ｕ（ｘ，ｒ）σ（ｘ，ｒ）． （１）式中：ｕ（ｘ，ｒ）为整个数据集坐标（ｘ ，ｒ）下幅度的平均值，σ（ｘ，ｒ）为对应幅度的标准差。通过选定合适的振幅离差因子阈值，对于振幅离差因子大于确定阈值的像元被选为 ＰＳ候选点。对于选定的 ＰＳ候选点而言，其干涉相位由以下成分构成：φ ＝φｍａｓ－φｓｌａ ＝φｆｌａｔ＋φｔｏｐｏ＋φｄｅｆ＋φｒｍｅ＋φｎｏｉｓｅ．（２）式中：φｆｌａｔ表示平地相位，可利用成像几何参数估计得出；φｄｅｆ表示形变相位，在时间基线很短的情况下，可忽略；φｒｍｅ表 示 ＲＭＥ 引起 的相位误 差，φｎｏｉｓｅ表示系统热噪 声 等 随 机 噪 声；φｔｏｐｏ表 示 地 形 相 位，可 表示为φｔｏｐｏ ＝ ＫｚＺｒｅａｌ．（３）式中：Ｚｒｅａｌ为雷达波在地表的相位中心高度，Ｋｚ表示干涉影像对应的相高转换参数：Ｋｚ ＝４πＢｃｏｓ（θ－α）λＲｓｉｎθ． （４）式中：Ｂ 为对应基线 长度，Ｒ 为 对 应 主 影 像 斜 距，λ为波长，θ 和α 分别为局部入射角和基线倾角。因此，对于时间基线短，去 除平地效 应后的多基线数据集中基线 Ｋ 对应的干涉相位为φｋ ＝ Ｋｚ，ｋＺｒｅａｌ＋φｒｍｅ，ｋ．（５）由上式可知，只要能求出地表相位中心Ｚｒｅａｌ，就可反解出 相 位 误 差φｒｍｅ，ｋ，完 成 相 位 定 标 工 作。而 根 据ＩｎＳＡＲ 测高原理可知，在多基线数据集中，即使轨道不同，干 涉 基 线 参 数 不 一 致，但 其 对 应 的Ｚｒｅａｌ一致，多基线数据给Ｚｒｅａｌ的求解提供多余观测，因此可采用 极 大 似 然 估 计 的 方 法 对 ＰＳ 候 选 点 的Ｚｒｅａｌ求解［１７］：γｋ（ｘ，ｒ）＝ ｍａｘｚ１Ｋ｜∑Ｋｋ＝１ｅｘｐ［φＫ（ｘ，ｒ）－Ｋｚ，ｋＺ］｜ ．（６）式中：φＫ为多基线影像的相位；Ｚ 为待求参数，搜索ｚ，使γ最大的ｚ 即为所求的Ｚｒｅａｌ；γ是该点在多基线数据集中相位稳定性的一个评价指标，γ 越接近于１，该点相位稳定性越高，求出来的Ｚｒｅａｌ越可靠，定标精度越高。通过设定γ 的阈值，最终选出高于γ 阈值的点为 ＰＳ点，利用 ＰＳ点解算的Ｚｒｅａｌ，通过式（４）反解出不同基线影像 ＰＳ点的相位误差，并将相位误差插值到整景影像 后 再进行相位 校正。该方法通过选取 ＰＳ点，通过在相位误差比较小的 ＰＳ点出估计出较为精准的地形高度，再利用地形高度反算出精准的相位误差，最后插值计算整景影像相位定标校正值。同时，利用 ＲＭＥ 在多基线数据集中呈现随机误差特性，探测多基线数据集中的 ＲＭＥ，但由于主影像中包含的 ＲＭＥ在各个干涉影像中是相·２５· 测　绘　工　程　　　　　　　　　 　 　　　　　　　　　第２７卷数据集中不呈现随机误差特性，而呈现系统误差特性，因此，该方法存在局限性，无法剔除主影像 ＲＭＥ的影响。１．２　基于多项式拟合的相位定标方法综上所述，定标后的干涉相位中仍然会保留主影像中的 ＲＭＥ。在条带成像模式下，ＳＡＲ 影像整个距离向 对 应 的 搭 载 平 台 运 动 参 数 是 一 致 的，而ＲＭＥ引起的相位波动可看成是搭载平台运动位置偏移量在雷达视线方向上的投影，这导致 ＲＭＥ 在距离向上会引起趋势变化的相位波动［１４］，而这种趋势变化，可采用一种高次线性多项式拟合的方法从干涉相位与外部 ＤＥＭ 的差分相位的距离向中估计并抑制消除该部分 ＲＭＥ。采用的多项式表达为［１９］φＲＭＥ（ｘ，ｙ）＝∑ｎｉ＝０ａｉｙｉ＋ａｈｈ． （７）其中，ｎ为多项式的阶数，ａｉ和ａｈ为模型的拟合参数，ｈ 为 外 部 ＤＥＭ 提 供 的 地 表 高 度。解 算 拟 合 参 数前，可考虑利用小波函数对定 标后相位滤 波，抑制地形残差等非线性相位成分对拟 合参数解算的影响，再利用最小二乘解出拟合 参数后，从 定标后相位中减去通过拟合参数计算的 ＲＭＥ，即可完成整个联合相位定标过程。多 项式拟合 的方法依赖于高精度的外部 ＤＥＭ，而通过基于 ＰＳ的多基线的相位定标方法 估 计 出 的Ｚｒｅａｌ能 有 效 代 替 高 精 度 的 外 部ＤＥＭ，因此，本文提出的联合定标方法具有不依赖高精度外部 ＤＥＭ 的优点。整个联合定标两步方法流程如图１所示。图１　基于 ＰＳ的多基线相位定标两步法流程２　实验分析２．１　实验数据本文采用的 Ｂｉｏｓａｒ２００８ 数据为位 于 瑞 典 北 部Ｋｒｙｃｋｌａｎ河流域的机载 Ｅ－ＳＡＲ 数据，海拔高度 范围为１１４～４０５ ｍ。其 中 本 文 选 用 了 获 取 于 ２００８－１０－１４的６景 ｐ波段时序 ＳＡＲ 数据，详细信息如表１所示。该ＳＡＲ 数据经多视处理后大小为６　４７２像元×１　５０１像元，距离向空间分辨率２．１２ｍ，方位向空间分 辨 率 １．６ ｍ。 图 ２（ａ）为 数 据 获 取 区 域 的Ｇｏｏｇｌｅ影像，示意实验区大概的地理位置及范围。图２（ｂ）为数据获取区域的 Ｐａｕｌｉ基分解影像。表１　影像数据集概况影像ＩＤ 波段 极化方式 获取日期 干涉基线０８ｂｉｏｓａｒ０１０３×１ Ｐ Ｆｕｌ　 ２００８－１０－１４ ＩｎＳＡＲ　Ｍａｓｔｅｒ　８ｍｌｅｆｔ０８ｂｉｏｓａｒ０１０５×１ Ｐ Ｆｕｌ　 ２００８－１０－１４ ＩｎＳＡＲ　Ｓｌａｖｅ　１６ｍｌｅｆｔ０８ｂｉｏｓａｒ０１０７×１ Ｐ Ｆｕｌ　 ２００８－１０－１４ ＩｎＳＡＲ　Ｓｌａｖｅ　２４ｍｌｅｆｔ０８ｂｉｏｓａｒ０１０９×１ Ｐ Ｆｕｌ　 ２００８－１０－１４ ＩｎＳＡＲ　Ｓｌａｖｅ　３２ｍｌｅｆｔ０８ｂｉｏｓａｒ０１１１×１ Ｐ Ｆｕｌ　 ２００８－１０－１４ ＩｎＳＡＲ　Ｓｌａｖｅ　４０ｍｌｅｆｔ０８ｂｉｏｓａｒ０１０１×１ Ｐ Ｆｕｌ　 ２００８－１０－１４ ＩｎＳＡＲ　Ｓｌａｖｅ　０ｍ第１２期　　　　　　　　　　　　　贺文杰，等：多基线机载 ＳＡＲ 联合相位定标两步法·３５·图２　实验区示意图２．２　实验结果及分析本文在对 Ｂｉｏｓａｒ２００８ 数 据集的相 位 定 标 处 理中选定振幅离差因子 Ｄ＞１５的点为 ＰＳ候选点，相位稳定指标γ高于０．８的点为最终的 ＰＳ点。联合定标过程中采用了三次多项式拟合多基线定标后数据中的 ＲＭＥ，最终结果如图３所示。图３是各个定标 过 程 中，各 干 涉 对 的 干 涉 相 位 与 ＬＩＤＡＲ 的ＤＥＭ 生成的 地形相位 的差分 相 位 图。因 为 ＲＭＥ在距离向上存在趋势变化，因此会在干涉相位中产生方位向的相位条纹，很容易在差分相位中辨识出来，图３　定标过程中的差分相位·４５· 测　绘　工　程　　　　　　　　　 　 　　　　　　　　　第２７卷多基线机载ＳＡＲ 联合相位定标两步法贺文杰，朱建军，付海强，王会强（中南大学 地球科学与信息物理学院，湖南 长沙 ４１００８３）摘　要：多基线机载 ＳＡＲ 数据相比单基线数据具有获取更丰富地表散射信息特点，已在地形测绘、植被监测领域广泛应用。在多基线数据应用前，需要对多基线数据集进行相位定标，校正数据集内由平台运动误差引起的相位偏差。本文提出一种多基线联合相位定标方法，解决传统永久散射体（Ｐｅｒｓｉｓｔｅｎｔ　Ｓｃａｔｔｅｒ，ＰＳ）相位定标方法无法剔除主影像残余运动误差（Ｒｅｓｉｄｕａｌ　Ｍｏｔｉｏｎ　Ｅｒｒｏｒｓ，ＲＭＥ）的问题。利用德国宇航局机载 Ｅ－ＳＡＲ 多基线数据对算法进行测试，试验结果表明，采用本文方法定标后，利用干涉相位估计的 ＤＥＭ 的精度为２．３０ｍ，优于传统相位定标方法解算的４．４７ｍ，有效解决传统定标方法无法兼顾主影像运动误差的问题。关键词：相位定标；多基线；永久散射体；机载 ＳＡＲ中图分类号：Ｐ２２８　　　文献标识码：Ａ　　　文章编号：１００６－７９４９（２０１８）１２－００５１－０６Ｍｕｌｔｉ－ｂａｓｅｌｉｎｅ　ａｉｒｂｏｒｎｅ　ＳＡＲ　ｐｈａｓｅ　ｃａｌｉｂｒａｔｉｏｎ　ｗｉｔｈ　ｔｗｏ－ｓｔｅｐ　ｍｅｔｈｏｄＨＥ　Ｗｅｎｊｉｅ，ＺＨＵ　Ｊｉａｎｊｕｎ，ＦＵ　Ｈａｉｑｉａｎｇ，ＷＡＮＧ　Ｈｕｉｑｉａｎｇ（Ｓｃｈｏｏｌ　ｏｆ　Ｇｅｏｓｃｉｅｎｃｅｓ　ａｎｄ　Ｉｎｆｏ－Ｐｈｙｓｉｃｓ，Ｃｅｎｔｒａｌ　Ｓｏｕｔｈ　Ｕｎｉｖｅｒｓｉｔｙ，Ｃｈａｎｇｓｈａ　４１００８３，Ｃｈｉｎａ）Ａｂｓｔｒａｃｔ：Ｔｈｅ　ｍｕｌｔｉ－ｂａｓｅｌｉｎｅ　ＳＡＲ　ｄａｔａｓｅｔｓ，ｄｕｅ　ｔｏ　ｔｈｅａｄｖａｎｔａｇｅｓ　ｏｆ　ｐｒｏｖｉｄｉｎｇ　ｒｅｄｕｎｄａｎｔ　ｏｂｓｅｒｖａｔｉｏｎ，ｈａｓｂｅｃｏｍｅ　ｔｈｅ　ｍａｉｎｓｔｒｅａｍ　ｔｒｅｎｄ　ｏｆ　ＳＡＲ　ｓｙｓｔｅｍ　ｄｅｖｅｌｏｐｍｅｎｔ．Ｐｒｉｏｒ　ｔｏ　ａｎｙ　ｄａｔａ　ｐｒｏｃｅｓｓｉｎｇ，ｉｔ　ｉｓ　ｎｅｃｅｓｓａｒｙ　ｔｏｃａｌｉｂｒａｔｅ　ｔｈｅ　ｐｈａｓｅ　ｉｎｆｏｒｍａｔｉｏｎ　ｉｎ　ｔｈｅ　ｄａｔａｓｅｔｓ　ｉｎ　ｏｒｄｅｒ　ｔｏ　ｒｅｍｏｖｅ　ｔｈｅ　ｐｈａｓｅ　ｅｒｒｏｒｓ　ｌｉｋｅ　ｒｅｓｉｄｕａｌ　ｍｏｔｉｏｎｅｒｒｏｒｓ，ｔｈｅｒｍａｌ　ｎｏｉｓｅ　ａｎｄ　ｓｏ　ｏｎ．Ｔｈｉｓ　ｐａｐｅｒ　ｐｒｏｐｏｓｅｓ　ａ　ｔｗｏ－ｓｔｅｐ　ｍｅｔｈｏｄ　ｔｏ　ｓｏｌｖｅ　ｔｈｅ　ｍａｓｔｅｒ　ｉｍａｇｅ　ＲＭＥｒｅｍａｉｎ　ｐｒｏｂｌｅｍ　ｃａｕｓｅｄ　ｂｙ　ｔｈｅ　ｔｒａｄｉｔｉｏｎａｌ　ｍｅｔｈｏｄ．Ｔｈｅ　ｎｅｗ　ｍｅｔｈｏｄ　ａｐｐｌｉｅｄ　ｔｏ　ｔｈｅ　ｐｈａｓｅ　ｃａｌｉｂｒａｔｉｏｎ　ｏｆ　Ｐｂａｎｄ　ＢｉｏＳＡＲ２００８ｄａｔａｓｅｔｓ　ｓｈｏｗｓ　ｔｈａｔ　ｔｈｅ　ＲＭＳＥ　ｏｆ　ＤＥＭ　ｅｒｒｏｒｓ　ｂｅｔｗｅｅｎ　ｔｈｅ　ＤＥＭ　ｍａｄｅ　ｂｙ　ＬＩＤＡＲ　ａｎｄｔｈｅ　ｏｔｈｅｒ　ｏｎｅ　ｉｎｖｅｒｓｅｄ　ｂｙ　ｉｎｔｅｒｆｅｒｏｍｅｔｒｉｃｐｈａｓｅ　ｃａｎ　ｂｅ　２．３０ ｍ ｗｉｔｈ　ｔｈｅ　ｐｒｏｐｏｓｅｄ　ｍｅｔｈｏｄ，ｗｉｔｈ　ｔｈｅｔｒａｄｉｔｉｏｎａｌ　ｍｅｔｈｏｄ　ａｓ　４．４７ｍ，ｂｅｉｎｇ　ｔｈｅ　ｐｒｏｐｏｓｅｄ　ｍｅｔｈｏｄ　ｔｈｕｓ　ｉｍｐｒｏｖｅｄ　ａ　ｌｏｔ．Ｋｅｙ　ｗｏｒｄｓ：ｐｈａｓｅ　ｃａｌｉｂｒａｔｉｏｎ，ｍｕｌｔｉ－ｂａｓｅｌｉｎｅ，ｐｅｒｓｉｓｔｅｎｔ　ｓｃａｔｔｅｒｅｒ，ａｉｒｂｏｒｎｅ　ＳＡＲ收稿日期：２０１８－０１－１７基金 项 目：国 家 自 然 科 学 基 金 资 助 项 目 （４１５３１０６８；４１３７１３３５；４１６７１３５６）；高分辨率对地观测系统重大专项基金资助项目（民用部分）（０３－Ｙ２０Ａ１１－９００１－１５／１６）第一作者简介：贺文杰（１９９３－），男，硕士研究生通信作者简介：朱建军（１９６２－），男，教授，博导．　　在机载ＳＡＲ 系统中，由于飞机平台在飞行过程中受气流摄动影响，会产生轨迹偏离和姿态变化等运动误差。这种运动误差的补偿通常依赖于飞机的导航定向定位系统（ＰＯＳ），通过记录实时成像的位置参数及姿态参数，在成像过程中对ＳＡＲ 影像进行运动补偿［１］。然而，现有导航设备精度指标无法满足运动补偿精度的需求，导致运动补偿后的 ＳＡＲ影像 仍 包 含 残 余 运 动 误 差 （Ｒｅｓｉｄｕａｌ　Ｍｏｔｉｏｎ　Ｅｒ－ｒｏｒｓ，ＲＭＥ）。它会造成 ＳＡＲ 影像出现散焦及相位波动［３］，其在雷达视线方向投影产生的相位误差信号，如果不加以定标校正，将直接掩 盖掉观测对象的细节相 位 信 号，从 而 影 响 后 续 数 据 处 理 及 相 关应用。传统的 ＲＭＥ 估计方法可分为 两 类［２］：一类利用子视ＳＡＲ 影像对应的 ＲＭＥ 不同，从子孔径的干涉相位中估计 ＲＭＥ。包括 Ｂｕｌｌｏｃｋ［４－５］等提出的双天线数据估计 ＲＭＥ 的方法，利用干涉对配准误差估计 ＲＭＥ的方法［７－８］，以及 Ｒｅｉｇｂｅｒ［６，９］提出的多斜视子孔径估计 ＲＭＥ方法等。另一类基于单景ＳＡＲ影像，利用自聚焦的方法估计 ＲＭＥ。包括 Ｃａｎｔａ－ｌｌｏｕｂｅ［１０］提出的利用孤立点目标进行相位跟踪估计标过程的差分相位图来评价其相位精度。图３（ａ）是定标前的差分相位，各干涉对均存在由 ＲＭＥ造成的相位条纹。图３（ｂ）是采用 ＰＳ多基线相位定标方法提取的相位误差，它满足各个干涉对互不相同且方位向变化的相位条纹的特征，可推断出该相位误差多由 ＲＭＥ产生。图３（ｃ）是多基线相位定标后数据集的差分相位，定标后数据集的相位精度较定标前有了较大的改善，但仍存在位置相同，变化特征一致相位条纹，可推断出其来源于主影像的 ＲＭＥ。传统多基线相位定标方法的实质是将各副影像参照主影像来校正各副影像中不同的相位误差，但主影像同样存在 ＲＭＥ 造成的相位误差，而主影像作为参考对象，其自身的相位误差利用该方法无法剔除，因此被保留在定标后的相位中。图３（ｄ）是本文提出的联合相位定标后的差分相位，上述相位条纹已被很好的抑制，可看到被掩盖的地形残差的相位细节。为了更直观 的评定出两种相位定标方法的相位定标精度，本文通过定标后干涉相位相高转换后的 ＤＥＭ 与实验区 ＬＩＤＡＲ 生成 ＤＥＭ 的地形残差的 ＲＭＳＥ来评定定标精度。各定标过程中的地形残差如图４所示，在不考虑相干性较低区域（取相干性＞０．８）的 情 况 下，各 干 涉 对 定 标 前 后 与 ＬＩＤＡＲ的 ＤＥＭ 的 ＲＭＳＥ如表２所示。由图４（ａ）可知，定标前相位精度 较 低，直接利用干涉 相位进行 ＤＥＭ反演，会存在数十米的误差，地表 起伏细节的相位信号被 ＲＭＥ 掩 盖。且 影 像 中 存 在 随 空 间 基 线 增大，残余运动误差导致的地形 残差变小的 现 象，这是由于残留运动误差引起的相位波 动量级一致的情况下，随空间基线 Ｂ 增大，相高转换参数Ｋｚ相应增大，因此，引起的误 差相位 转化成的地 形残差变小，基线较长的干涉对抗差能力强。图４（ｂ）是经过多基线相位定标后的地形残差图，在该数据集的影像上半部分校正较好，可以看到地形残差的细节信息，但仍存在量级较大且位置、量级均相同的 ＲＭＥ，参照多基线相位定标过程，可知其由主影像的 ＲＭＥ导致。图４（ｃ）是本文提出的联合相位定标方法处理后的地形残 差图，ＲＭＥ 导致的趋 势相位条纹被很好的抑制消除，地形残差信息保存相对完整。对照表２可知，提出的新方法在各个干涉对均有很好的效果，能极大提高多基线数据集的相位精度。图４　定标前后的 ＤＥＭ 误差第１２期　　　　　　　　　　　　　贺文杰，等：多基线机载 ＳＡＲ 联合相位定标两步法·５５·表２　各干涉对定标前后地形残差的 ＲＭＳＥ　 ｍ　干涉对编号 未定标 已有定标方法定标 联合相位定标０１０１－０１０３　 １５．４５　 ４．９９　 ３．２１０１０５－０１０３　 １４．２３　 ５．３７　 ３．８００１０７－０１０３　 ６．１０　 ４．６２　 ２．６００１０９－０１０３　 ６．２２　 ４．５２　 ２．４００１１１－０１０３　 ６．０８　 ４．４７　 ２．３０３　结　论针对多基线相位定标后主影像 ＲＭＥ 的残留问题，本文提出了一种联合相位定标方法。该方法通过引入高次线性多项式从干涉相位与地形相位的差分相位中来拟合存在距离向趋势变化的主影像ＲＭＥ，校正残留的主影像 ＲＭＥ，且新方法具有不依赖于外部 ＤＥＭ 的优点。实验表明，该方法能有效提高多基线数据集的相位 精度。但仍存在以下问题：在校正多基线数据集 ＲＭＥ 的同时会造成数据集中因 轨 道 位 置 变 化 而 产 生 的 相 位 差 异 信 号 的损失。参考文献：［１］　ＨＥＮＳＬＥＹ　Ｓ，ＭＩＣＨＥＬ　Ｔ，ＳＩＭＡＲＤ　Ｍ，ｅｔ　ａｌ．Ｒｅ－ｓｉｄｕａｌ　ｍｏｔｉｏｎ　ｅｓｔｉｍａｔｉｏｎ　ｆｏｒ　ＵＡＶＳＡＲ：Ｉｍｐｌｉｃａｔｉｏｎｓｏｆ　ａｎ　ｅｌｅｃｔｒｏｎｉｃａｌｌｙ　ｓｃａｎｎｅｄ　ａｒｒａｙ［Ｃ］／／ＩＥＥＥ　ＲａｄａｒＣｏｎｆｅｒｅｎｃｅ．ＩＥＥＥ，２００９：１－５．［２］　钟雪莲，向茂生，郭华东，等．机载重轨干涉合成孔径雷达的发展［Ｊ］．雷达学报，２０１３，２（３）：３６７－３８１．［３］　ＧＲＡＹ　Ａ　Ｌ，ＦＡＲＲＩＳ－ＭＡＮＮＩＮＧ　Ｐ　Ｊ．Ｒｅｐｅａｔ－ｐａｓｓｉｎｔｅｒｆｅｒｏｍｅｔｒｙ　ｗｉｔｈ　ａｉｒｂｏｒｎｅ　ｓｙｎｔｈｅｔｉｃ　ａｐｅｒｔｕｒｅ　ｒａｄａｒ［Ｊ］．Ｊｏｕｒｎａｌ　ｏｆ　Ｅｌｅｃｔｒｏｎｉｃｓ　＆ Ｉｎｆｏｒｍａｔｉｏｎ　Ｔｅｃｈｎｏｌｏ－ｇｙ，２００３，３１（１）：１８０－１９１．［４］　ＢＵＬＬＯＣＫ　Ｒ　Ｊ，ＶＯＬＥＳ　Ｒ，ＣＵＲＲＩＥ　Ａ，ｅｔ　ａｌ．Ｔｗｏ－ｌｏｏｋ　ｍｅｔｈｏｄ　ｆｏｒ　ｒｏｌｌ　ｃｏｒｒｅｃｔｉｏｎ　ｏｆ　ａｉｒｃｒａｆｔ－ｂｏｒｎｅ　ｉｎｔｅｒ－ｆｅｒｏｍｅｔｒｉｃ　ＳＡＲ　ｄａｔａ［Ｃ］／／ Ｒａｄａｒ　Ｉｎｔｅｒｆｅｒｏｍｅｔｒｙ．１９９７：２／１ －２／６．［５］　ＢＵＬＬＯＣＫ　Ｒ　Ｊ，ＢＲＥＮＮＡＮ　Ｐ　Ｖ，ＧＲＩＦＦＩＴＨＳ　Ｈ　Ｄ．Ｔｗｏ－ｌｏｏｋ　ｒｏｌｌ　ｃｏｍｐｅｎｓａｔｉｏｎ　ｆｏｒ　ａｉｒｃｒａｆｔ－ｂｏｒｎｅ　ｉｎｔｅｒｆｅｒ－ｏｍｅｔｒｉｃ　ＳＡＲ　ｗｉｔｈｏｕｔ　ｐｈａｓｅ　ｕｎｗｒａｐｐｉｎｇ［Ｊ］．Ｅｌｅｃｔｒｏｎ－ｉｃｓ　Ｌｅｔｔｅｒｓ，２００５，４１（６）：３６７－３６８．［６］　ＰＲＡＴＳ　Ｐ，ＭＡＬＬＯＲＱＵＩ　Ｊ　Ｊ．Ｅｓｔｉｍａｔｉｏｎ　ｏｆ　ａｚｉｍｕｔｈｐｈａｓｅ　ｕｎｄｕｌａｔｉｏｎｓ　ｗｉｔｈ　ｍｕｌｔｉｓｑｕｉｎｔ　ｐｒｏｃｅｓｓｉｎｇ　ｉｎ　ａｉｒ－ｂｏｒｎｅ　ｉｎｔｅｒｆｅｒｏｍｅｔｒｉｃ　ＳＡＲ　ｉｍａｇｅｓ［Ｊ］．ＩＥＥＥ　Ｔｒａｎｓａｃ－ｔｉｏｎｓ　ｏｎ　Ｇｅｏｓｃｉｅｎｃｅ　＆ Ｒｅｍｏｔｅ　Ｓｅｎｓｉｎｇ，２００３，４１（６）：１５３０－１５３３．［７］　ＨＥＮＳＬＥＹ　Ｓ，ＭＩＣＨＥＬ　Ｔ，ＳＩＭＡＲＤ　Ｍ，ｅｔ　ａｌ．Ｒｅ－ｓｉｄｕａｌ　ｍｏｔｉｏｎ　ｅｓｔｉｍａｔｉｏｎ　ｆｏｒ　ＵＡＶＳＡＲ：Ｉｍｐｌｉｃａｔｉｏｎｓｏｆ　ａｎ　ｅｌｅｃｔｒｏｎｉｃａｌｌｙ　ｓｃａｎｎｅｄ　ａｒｒａｙ［Ｃ］／／ＩＥＥＥ　ＲａｄａｒＣｏｎｆｅｒｅｎｃｅ．ＩＥＥＥ，２００９：１－５．［８］　ＳＣＨＥＩＢＥＲ　Ｒ，ＭＯＲＥＩＲＡ　Ａ．Ｃｏｒｅｇｉｓｔｒａｔｉｏｎ　ｏｆ　ｉｎｔｅｒ－ｆｅｒｏｍｅｔｒｉｃ　ＳＡＲ　ｉｍａｇｅｓ　ｕｓｉｎｇ　ｓｐｅｃｔｒａｌ　ｄｉｖｅｒｓｉｔｙ［Ｊ］．ＩＥＥＥ　Ｔｒａｎｓａｃｔｉｏｎｓ　ｏｎ　Ｇｅｏｓｃｉｅｎｃｅ　＆ Ｒｅｍｏｔｅ　Ｓｅｎｓｉｎｇ，２０００，３８（５）：２１７９－２１９１．［９］　ＲＥＩＧＢＥＲ　Ａ，ＰＲＡＴＳ　Ｐ，ＭＡＬＬＯＲＱＵＩ　Ｊ　Ｊ．Ｒｅｆｉｎｅｄｅｓｔｉｍａｔｉｏｎ　ｏｆ　ｔｉｍｅ－ｖａｒｙｉｎｇ　ｂａｓｅｌｉｎｅ　ｅｒｒｏｒｓ　ｉｎ　ａｉｒｂｏｒｎｅＳＡＲ　ｉｎｔｅｒｆｅｒｏｍｅｔｒｙ［Ｊ］．ＩＥＥＥ　Ｇｅｏｓｃｉｅｎｃｅ　＆ ＲｅｍｏｔｅＳｅｎｓｉｎｇ　Ｌｅｔｔｅｒｓ，２００６，３（１）：１４５－１４９．［１０］ＣＡＮＴＡＬＬＯＵＢＥ　Ｈ，ＤＵＢＯＩＳ－ＦＥＲＮＡＮＤＥＺ　Ｐ．Ａｉｒ－ｂｏｒｎｅ　Ｘ－ｂａｎｄ　ＳＡＲ　ｉｍａｇｉｎｇ　ｗｉｔｈ　１０ ｃｍ　ｒｅｓｏｌｕｔｉｏｎ：ｔｅｃｈｎｉｃａｌ　ｃｈａｌｌｅｎｇｅ　ａｎｄ　ｐｒｅｌｉｍｉｎａｒｙ　ｒｅｓｕｌｔｓ［Ｊ］．２００６，１５３（２）：１６３－１７６．［１１］ＭＡＣＥＤＯ　Ｋ　Ａ　Ｃ　Ｄ，ＳＣＨＥＩＢＥＲ　Ｒ，ＭＯＲＥＩＲＡ　Ａ．Ａｎ　Ａｕｔｏｆｏｃｕｓ　Ａｐｐｒｏａｃｈ　ｆｏｒ　Ｒｅｓｉｄｕａｌ　Ｍｏｔｉｏｎ　ＥｒｒｏｒｓＷｉｔｈ　Ａｐｐｌｉｃａｔｉｏｎ　ｔｏ　Ａｉｒｂｏｒｎｅ　Ｒｅｐｅａｔ－Ｐａｓｓ　ＳＡＲ　Ｉｎｔｅｒ－ｆｅｒｏｍｅｔｒｙ［Ｊ］．ＩＥＥＥ　Ｔｒａｎｓａｃｔｉｏｎｓ　ｏｎ　Ｇｅｏｓｃｉｅｎｃｅ　＆Ｒｅｍｏｔｅ　Ｓｅｎｓｉｎｇ，２００８，４６（１０）：３１５１－３１６２．［１２］ＺＨＯＮＧ　Ｘ，ＸＩＡＮＧ　Ｍ，ＹＵＥ　Ｈ，ｅｔ　ａｌ．Ａｌｇｏｒｉｔｈｍ　ｏｎｔｈｅ　Ｅｓｔｉｍａｔｉｏｎ　ｏｆ　Ｒｅｓｉｄｕａｌ　Ｍｏｔｉｏｎ　Ｅｒｒｏｒｓ　ｉｎＡｉｒｂｏｒｎｅＳＡＲ　Ｉｍａｇｅｓ［Ｊ］．ＩＥＥＥ　Ｔｒａｎｓａｃｔｉｏｎｓ　ｏｎ　Ｇｅｏｓｃｉｅｎｃｅ　＆Ｒｅｍｏｔｅ　Ｓｅｎｓｉｎｇ，２０１４，５２（２）：１３１１－１３２３．［１３］ＭＯＲＲＩＳＯＮ　Ｒ　Ｌ，ＤＯ　Ｍ　Ｎ，ＭＵＮＳＯＮ　Ｄ　Ｃ．ＳＡＲ　Ｉｍ－ａｇｅ　Ａｕｔｏｆｏｃｕｓ　Ｂｙ　Ｓｈａｒｐｎｅｓｓ　Ｏｐｔｉｍｉｚａｔｉｏｎ：ＡＴｈｅｏｒｅｔｉ－ｃａｌ　Ｓｔｕｄｙ［Ｊ］．ＩＥＥＥ　Ｔｒａｎｓａｃｔｉｏｎｓ　ｏｎ　Ｉｍａｇｅ　Ｐｒｏｃｅｓｓ－ｉｎｇ，２００７，１６（９）：２３０９－２３２１．［１４］ＦＵ　Ｈ　Ｑ，ＺＨＵ　Ｊ　Ｊ，ＷＡＮＧ　Ｃ　Ｃ，ｅｔ　ａｌ．Ａ　ＷａｖｅｌｅｔＤｅｃｏｍｐｏｓｉｔｉｏｎ　ａｎｄ　Ｐｏｌｙｎｏｍｉａｌ　Ｆｉｔｔｉｎｇ－Ｂａｓｅｄ　Ｍｅｔｈｏｄｆｏｒ　ｔｈｅ　Ｅｓｔｉｍａｔｉｏｎ　ｏｆ　Ｔｉｍｅ－Ｖａｒｙｉｎｇ　Ｒｅｓｉｄｕａｌ　ＭｏｔｉｏｎＥｒｒｏｒ　ｉｎ　Ａｉｒｂｏｒｎｅ　Ｉｎｔｅｒｆｅｒｏｍｅｔｒｉｃ　ＳＡＲ ［Ｊ］．ＩＥＥＥＴｒａｎｓａｃｔｉｏｎｓ　ｏｎ　Ｇｅｏｓｃｉｅｎｃｅ　＆ Ｒｅｍｏｔｅ　Ｓｅｎｓｉｎｇ，２０１７，ＰＰ（９９）：１－１１．［１５］ＰＡＲＤＩＮＩ　Ｍ．Ａｄｖａｎｃｅｓ　ａｎｄ　Ｅｘｐｅｒｉｍｅｎｔｓ　ｏｆ　Ｔｏｍｏ－ｇｒａｐｈｉｃ　ＳＡＲ　Ｉｍａｇｉｎｇ　ｆｏｒ　ｔｈｅ　Ａｎａｌｙｓｉｓ　ｏｆ　ＣｏｍｐｌｅｘＳｃｅ－ｎａｒｉｏｓ，Ｐｈ．Ｄ．Ｔｈｅｓｉｓ， Ｕｎｉｖｅｒｓｉｔｙ　ｏｆ　Ｐｉｓａ，Ｉｔａｌｙ，Ｊｕｎｅ　２０１０．［１６］ＴＥＢＡＬＤＩＮＩ　Ｓ，ＧＵＡＲＮＩＥＲＩ　Ａ　Ｍ．Ｏｎ　ｔｈｅ　Ｒｏｌｅ　ｏｆＰｈａｓｅ　Ｓｔａｂｉｌｉｔｙ　ｉｎ　ＳＡＲ　Ｍｕｌｔｉｂａｓｅｌｉｎｅ　Ａｐｐｌｉｃａｔｉｏｎｓ［Ｊ］．ＩＥＥＥ　Ｔｒａｎｓａｃｔｉｏｎｓ　Ｏｎ　Ｇｅｏｓｃｉｅｎｃｅ　Ａｎｄ　Ｒｅｍｏｔｅ　Ｓｅｎｓ－ｉｎｇ，２０１０，４８（７）：２９５３－２９６６．（下转第６３页）·６５· 测　绘　工　程　　　　　　　　　 　 　　　　　　　　　第２７卷３　结束语本文结合电子政务及地理信息的特点，提出一种更为灵活的权限管理方法：将部分资源封装成资源集，再 通 过 角 色 分 层 为 每 个 角 色 分 配 一 个 资 源集，当角色资源权限发生变化时，可以直接对其资源集进行服务的增删或对其中相关服务进行修改，而不需重新为角色分配成百上千个服务；研究一种固有权限与动态权限相结合的权限修改方法；同时将权限分配与部门所管辖区域关联起来，提出一种空间相 关 的 权 限 设 置 策 略。 本 文 以 襄 阳 市 “一 张图”空间信息云平台为例，详细说明本文权限管理策略的具体实现，依据该方法开发的权限管理软件模型的应用表明：本文提出的权限模型具有可实现性与可操作性，权限修改灵活，利于资源共享与协同办公。参考文献：［１］　田捷，侯海 波．基 于 ＷｅｂＧＩＳ的 一 站 式 电 子 政 务［Ｊ］．测绘科学，２００３，２８（１）：３３－３５．［２］　ＤＯＮＧ　Ｗ　Ｈ，ＬＩＵ　Ｊ　Ｐ，ＧＵＯ　Ｑ　Ｓ．Ｃｏｎｓｔｒｕｃｔｉｏｎ　ｏｆ　Ｅ－Ｇｏｖｅｒｎｍｅｎｔ　ＧＩＳ　ｂａｓｅｄ　ｏｎ．Ｎｅｔ　ｐｌａｔｆｏｒｍ　ａｎｄ　ＷｅｂＳｅｒｖｉｃｅ［Ｊ］．Ｊｏｕｒｎａｌ　ｏｆ　Ｌｉａｏｎｉｎｇ　Ｔｅｃｈｎｉｃａｌ　Ｕｎｉｖｅｒｓｉｔｙ，２００７，６７（２）：９２１－９２３．［３］　菅建华，董志．政务电子地图服务平台的设计与实现［Ｊ］．地理空间信息，２０１０，８（３）：１３４－１３７．［４］　曹先革，杨金 玲，黄 声 享，等．基 于 ＷｅｂＧＩＳ的 城 市 交通应急救 援 指 挥 系 统 设 计 与 实 现 ［Ｊ］．测 绘 工 程，２０１１，２０（２）：４６－４８．［５］　周海，李宏伟，杜泽欣，等．面向服务的城市管网综合管理系统设计［Ｊ］．测绘工程，２０１５，２４（３）：６９－７３．［６］　金兵兵，马桂云，刘洋，等．“多规合一”信息平台设计及关键技术研究［Ｊ］．测绘工程，２０１７，２６（７）：５２－５６．［７］　ＦＥＲＲＡＩＯＬＯ　Ｄ，ＣＵＧＩＮＩ　Ｊ，ＫＵＨＮ　Ｄ　Ｒ．Ｒｏｌｅ－ｂａｓｅｄａｃｃｅｓｓ　ｃｏｎｔｒｏｌ （ＲＢＡＣ）：ｆｅａｔｕｒｅｓ　ａｎｄ　ｍｏｔｉｖａｔｉｏｎｓ［Ｃ］／／Ｐｒｏｃｅｅｄｉｎｇｓ　ｏｆ　１１ｔｈ　Ａｎｎｕａｌ　Ｃｏｍｐｕｔｅｒ　ＳｅｃｕｒｉｔｙＡｐｐｌｉｃａｔｉｏｎ　Ｃｏｎｆｅｒｅｎｃｅ． Ｍｉｃｈｉｇａｎ：ＩＥＥＥ　ＣｏｍｐｕｔｅｒＳｏｃｉｅｔｙ　Ｐｒｅｓｓ，１９９５：２４１－２４８．［８］　蒋富伟，黄勤，刘益良，等．电子政务中基于角色的交叉集散授权策略［Ｊ］．信息安全与通信保密，２００６（６）：３９－４１．［９］　ＬＥＥ　Ｍ，ＣＨＯ　Ｎ　Ｄ，ＧＡＴＴＯＮ　Ｔ　Ｍ．Ａ　ｆｕｎｃｔｉｏｎ－ｂａｓｅｄｕｓｅｒ　ａｕｔｈｏｒｉｔｙ　ｄｅｌｅｇａｔｉｏｎ　ｍｏｄｅｌ［Ｊ］．Ｉｎｆｏｒｍａｔｉｏｎ　Ｓｃｉ－ｅｎｃｅｓ，２０１０，１８０（５）：７６５－７７５．［１０］赵再军．基于角色权限分配的协同电子政务访问控制模型研究［Ｊ］．计算机科学，２０１０，３７（８）：１４３－１４５．［１１］盘莉莉．电子政务中通用 ＲＢＡＣ模块的实现［Ｊ］．计算机工程与设计，２０１０，３１（１４）：３２１０－３２１３．［１２］胡领，毕利．改进的 ＲＢＡＣ 模型在电子政务系统中的应 用 研 究 ［Ｊ］．太 原 理 工 大 学 学 报，２０１４，４５（３）：３６３－３６８．［１３］王卫国，闾国年．电子政务信息系统的动态权限管理研究［Ｊ］．地球信息科学学报，２００５，７（４）：６７－７０．［１４］袁煜锋，吴健平．角 色 权 限 管 理 在 ＷｅｂＧＩＳ系 统 中 的应用［Ｊ］．计算机安全，２００８（１２）：３９－４２．［１５］胡腾波，叶建栲．适用于 ＷｅｂＧＩＳ的 ＲＢＡＣ 扩展模型研究［Ｊ］．计算机应用与软件，２００９，２６（６）：１８５－１８６．［责任编辑：李铭娜］（上接第５６页）［１７］ＦＥＲＲＥＴＴＩ　Ａ，ＰＲＡＴＩ　Ｃ，ＲＯＣＣＡ　Ｆ．ＰｅｒｍａｎｅｎｔＳｃａｔｔｅｒｅｒｓ　ｉｎ　ＳＡＲ　Ｉｎｔｅｒｆｅｒｏｍｅｔｒｙ［Ｊ］．ＩＥＥＥ　Ｔｒａｎｓａｃ－ｔｉｏｎｓ　ｏｎＧｅｏｓｃｉｅｎｃｅ　ａｎｄ　Ｒｅｍｏｔｅ　Ｓｅｎｓｉｎｇ，２００１，３９（１）：８－２０．［１８］ＢＩＡＮＣＯ　Ｖ，ＰＡＲＤＩＮＩ　Ｍ，ＰＡＰＡＴＨＡＮＡＳＳＩＯＵ　Ｋ，ｅｔ　ａｌ．Ｐｈａｓｅ　Ｃａｌｉｂｒａｔｉｏｎ　ｏｆ　Ｍｕｌｔｉｂａｓｅｌｉｎｅ　ＳＡＲ　ＤａｔａＳｔａｃｋｓ：Ａ　Ｍｉｎｉｍｕｍ　Ｅｎｔｒｏｐｙ　Ａｐｐｒｏａｃｈ［Ｃ］／／ Ｇｏｌｄ．ＤＬＲ，２０１２：１－４．［１９］ＦＵ　Ｈ，ＺＨＵ　Ｊ，ＷＡＮＧ　Ｃ，ｅｔ　ａｌ．Ｕｎｄｅｒｌｙｉｎｇ　Ｔｏｐｏｇ－ｒａｐｈｙ　Ｅｓｔｉｍａｔｉｏｎ　ｏｖｅｒ　Ｆｏｒｅｓｔ　Ａｒｅａｓ　Ｕｓｉｎｇ　Ｈｉｇｈ－Ｒｅｓｏ－ｌｕｔｉｏｎ　Ｐ－Ｂａｎｄ　Ｓｉｎｇｌｅ－Ｂａｓｅｌｉｎｅ　ＰｏｌＩｎＳＡＲ　Ｄａｔａ［Ｊ］．Ｒｅ－ｍｏｔｅ　Ｓｅｎｓｉｎｇ，２０１７，９（９）．［２０］ＰＡＲＤＩＮＩ　Ｍ．Ａｄｖａｎｃｅｓ　ａｎｄ　Ｅｘｐｅｒｉｍｅｎｔｓ　ｏｆ　Ｔｏｍｏ－ｇｒａｐｈｉｃ　ＳＡＲ　Ｉｍａｇｉｎｇ　ｆｏｒ　ｔｈｅ　Ａｎａｌｙｓｉｓ　ｏｆ　ＣｏｍｐｌｅｘＳｃｅ－ｎａｒｉｏｓ［Ｃ］．Ｐｈ．Ｄ．Ｔｈｅｓｉｓ，Ｕｎｉｖｅｒｓｉｔｙ　ｏｆ　Ｐｉｓａ，Ｉｔａｌｙ，Ｊｕｎｅ　２０１０．［责任编辑：李铭娜］第１２期　　　　　　　　　　　周思嘉，等：电子政务空间信息平台的权限管理方法·３６·