(19)国家知识产权局 (12)发明 专利申请 (10)申请公布号 (43)申请公布日 (21)申请 号 202211010660.9 (22)申请日 2022.08.23 (71)申请人 自然资源部国土卫星遥感应用中心 地址 100048 北京市海淀区紫竹院百胜 村1 号 (72)发明人 郭莉　王霞　唐新明　付冬暇　 刘书含　岳明宇　刘佳星　董立停　 黄鹏飞　赵小满　 (74)专利代理 机构 北京智乾知识产权代理事务 所(普通合伙) 11552 专利代理师 王晋 (51)Int.Cl. G06T 5/00(2006.01) G06T 7/10(2017.01) G06T 11/00(2006.01)G06V 10/46(2022.01) G06V 10/74(2022.01) (54)发明名称 制作大区域星载光学和SAR影像DOM的方法 和系统 (57)摘要 本发明提出制作大区域星 载光学和SAR影像 DOM的方法和系统。 其中， 方法包括： 基于资源三 号数字正射影像和外部已有DEM数据， 分别利用 星载光学和SAR影像与数字正射影像间特征和相 位一致性匹配的原理， 得到光学和SAR影像与数 字正射影像库间对应的像点坐标， 随后通过提取 的像控点和外部已有DEM数据参与光学和SAR联 合区域网平差和正射纠正处理， 最后通过匀光匀 色和镶嵌裁切得到星载光学和SAR影像DOM产品。 本发明减少人工干预带来的效率低下或人工选 点带来的刺点误差， 可实现大区域星载光学和 SAR影像自动和高精度纠 正处理； 有效保障了大 区域光学和SAR影 像联合DOM的产品质量。 权利要求书3页 说明书9页 附图16页 CN 115423702 A 2022.12.02 CN 115423702 A 1.一种制作大区域星载光学和SAR影 像DOM的方法， 其特 征在于， 所述方法包括： 步骤S1、 选取测区内原始星载光学影像及RPC模型、 原始星载SAR影像及RD模型、 数字正 射影像库和外 部参考DEM数据； 步骤S2、 对所述原始星载SA R影像进行多视和滤波 处理， 得到多视滤波 处理后SA R影像； 对所述RD模型开展RPC方程 替代得到拟合后RPC模型； 步骤S3、 基于光学影像特征一致性原则， 开展所述原始星载光学影像和数字正射影像 间像控点的匹配， 得到原 始星载光学影 像控制点； 步骤S4、 利用异源影像相位一致性原则， 开展所述多视滤波处理后SAR影像和数字正射 影像间像控点的匹配， 得到多视滤波处 理后SAR影 像控制点； 步骤S5、 分别利用所述原始星载光学影像控制点和多视滤波处理后SAR影像控制点， 基 于所述原始星载光学影像及RPC模型， 所述多视滤波处理后SAR影像和拟合后RPC模型开展 联合区域网平差和正 射纠正处 理， 得到正 射纠正后光学和SAR影 像； 步骤S6、 对所述正射纠正后光学和SAR影像开展自动匀光和匀色处理， 得到匀光匀色处 理后光学和SAR影 像； 步骤S7、 对所述匀光匀色处理后光学和SAR影像进行镶嵌处理， 得到镶嵌后光学和SAR 正射纠正产品； 步骤S8、 对所述镶嵌后光学和SA R正射纠正产品按标准分 幅进行裁切处理， 得到标准分 幅光学和SAR影 像DOM产品。 2.根据权利 要求1所述的一种制作大区域星载光学和SAR影像DOM的方法， 其特征在于， 在所述步骤S1中， 所述方法还 包括： 检查筛选所述原始星载光学影像和原始星载SAR影像的数据覆盖度和影像质量， 确定 数字正射影像库和 外部参考DEM数据与原始星载光学影像和原始星载SAR影像的匹配套合 程度。 3.根据权利 要求1所述的一种制作大区域星载光学和SAR影像DOM的方法， 其特征在于， 在所述步骤S2中， 所述对 所述原始星 载SAR影像进 行多视和滤波处理， 得到多视滤波处理后 SAR影像的方法包括： 步骤S21、 对原始星载SAR影像单视斜距复数据方位向或距离向做平均， 得到多视处理 后SAR影像； 然后通过自适应滤波、 中值滤波和Lee滤波 方法消除所述多视处理后SAR影像的 干涉图中系统性的热噪声、 由叠掩、 遮挡和时相 变化引起的去相关噪声以及局部匹配不准 引起的噪声， 得到所述多视滤波处 理后SAR影 像； 步骤S22、 对所述RD模型开展RPC方程替代得到拟合后RPC模型的方法包括： 结合所述 RPC模型的特性， 利用SAR严格成像几何模 型建立地面点的立体空间格网和影像面之间的对 应关系作为拟合的控制点， 求解RPC参数； RPC参数求解选用与地形无关和与地形相关两种 解算方案， 得到所述拟合后RPC模型。 4.根据权利 要求1所述的一种制作大区域星载光学和SAR影像DOM的方法， 其特征在于， 在所述步骤S3中， 所述基于光学影像特征一致性原则， 开展所述原始星载光学影像和数字 正射影像间像控点的匹配， 得到原 始星载光学影 像控制点的方法包括： 利用所述原始星载光学影像和数字正射影像库， 基于光学影像间特征一致性原则， 通 过相关系数法、 SIFT算法和SGM影像匹配算法获取预设数量的同名点， 得到所述原始星载光权　利　要　求　书 1/3 页 2 CN 115423702 A 2学影像控制点。 5.根据权利 要求1所述的一种制作大区域星载光学和SAR影像DOM的方法， 其特征在于， 在所述步骤S5中， 所述分别利用所述原始星载光学影像控制点和多视滤波处理后SAR影像 控制点， 基于所述原始星 载光学影像及RP C模型， 所述多视滤波处理后SAR影像和拟合后RP C 模型开展联合区域网平差和正 射纠正处 理， 得到正 射纠正后光学和SAR影 像的方法包括： 步骤S51、 使用所述原始星载光学影像控制点和多视滤波处理后SAR影像控制点， 通过 影像间控制点的约束关系补偿RPC模型的系统误差； 式中， (x， y)是原始星载光学影像和多视滤波处理后SAR影像上的量测坐标,(sample， line)为地面控制点利用RPC参数投影到影 像面的投影值； 对每个控制点列如下线性方程， 进行光束法平差: V＝At+BX‑l t＝(Δpx0 Δpx1 Δpx2 Δpy0 Δpy1 Δpy2)T， x＝(Δlat  Δlon Δheitht)T， 其中， V是残差， A是关于仿射参数的系数矩阵， B是物方坐标的系数矩阵， t为仿射参数 改正数， X为物方坐标改正数， l 为常数项部分； 其中， RPC模型为原 始光学影 像RPC模型和SAR影 像拟合而成； 步骤S52、 基于所述外部参考DEM数据、 原始星载光学和SAR影像区域网平差后的RPC参 数， 对原始星载光学和多视滤波处 理后SAR影 像进行正 射纠正。 6.根据权利 要求1所述的一种制作大区域星载光学和SAR影像DOM的方法， 其特征在于， 在所述步骤S7中， 所述对 所述匀光匀色处理后光学和SAR影像进 行镶嵌处理， 得到镶嵌后光 学和SAR正 射纠正产品的方法包括： 在所述匀光匀色处理后光学和SAR影像的预设定区域内， 增加闪避点， 调整镶嵌线两侧 局部区域影像的色调； 通过增加闪避点和编辑镶嵌线， 使预设定区域内的影像间的色调均 匀过渡。 7.根据权利 要求1所述的一种制作大区域星载光学和SAR影像DOM的方法， 其特征在于， 在所述步骤S8中， 所述标准分幅内图廓的四个角点像元中心点坐标的最小外接矩形， 以所 述最小外接矩形向外扩展预设置数值P个像素； 所述标准分幅内图廓的四个角点的像元中 心点坐标计算公式如下： Xmin＝int[min(X1， X2， X3， X4)/d] ×d‑P×d Ymin＝int[min(Y1， Y2， Y 3， Y4)/d] ×d‑P×d权　利　要　求　书 2/3 页 3 CN 115423702 A 3