(19)国家知识产权局 (12)发明 专利 (10)授权公告 号 (45)授权公告日 (21)申请 号 202211194499.5 (22)申请日 2022.09.29 (65)同一申请的已公布的文献号 申请公布号 CN 115291215 A (43)申请公布日 2022.11.04 (73)专利权人 中国科学院空天信息创新研究院 地址 100190 北京市海淀区北四环西路19 号 (72)发明人 邓云凯　向卫　刘大成　王吉利　 张勇　 (74)专利代理 机构 北京科迪生专利代理有限责 任公司 1 1251 专利代理师 邓治平 (51)Int.Cl. G01S 13/90(2006.01) G01S 7/02(2006.01) G01B 21/32(2006.01) G01C 5/00(2006.01) G06T 7/13(2017.01)G06T 7/73(2017.01) (56)对比文件 JP 2020176852 A,2020.10.2 9 CN 108919 262 A,2018.1 1.30 CN 111366925 A,2020.07.0 3 EP 2392943 A1,201 1.12.07 CN 114612631 A,202 2.06.10 耿诗画.基 于多源数据的采空区地表裂缝活 动监测. 《万方学术 论文全文在线数据库》 .202 2, 第1-92页. Claudio De Luca 等.Com ments on “Study of Systematic Bias i n Measuri ng Surface Deformati on With SAR I nterferometry ”. 《IEEE Transacti ons on Geoscience and Remote Sensi ng》 .2021,第1- 5页. (续) 审查员 许兆山 (54)发明名称 基于升降轨SAR卫星的长时序二维形变快速 解算方法 (57)摘要 本发明涉及一种基于升降轨SAR卫星的长时 序二维形变快速解算方法， 包 括步骤： S1、 采用时 序InSAR方法分别解算升降轨卫星SAR数据， 得到 升降轨视线向时序形变图集； S2、 比较升降轨时 序形变图， 确定公共区域的坐标范围； S3、 根据公 共范围坐标分别提取范围约束后的升降轨时序 形变图集； S4、 选择升轨或 降轨任意一个时序形 变图集按照时间间隔为1天进行三次样条插值； S5、 由另一未选中进行插值的时序形变图集， 得 到其存在累积形变的时间点， 并联合插值后的时 序形变图集提取升降轨公共时间点及其对应的 时序形变图集； S6、 将时间跨度统一后的升降轨 时序形变图集分别与时序上第一幅形变图做差，得到时序形变图集； S7、 分别进行解算， 得到最终 具有统一时间跨度基准的时序二维形变 。 [转续页] 权利要求书2页 说明书6页 附图10页 CN 115291215 B 2022.12.20 CN 115291215 B (56)对比文件 黄兆欢.基 于真实二维I nSAR技术的北麓河 多年冻土区地表形变时空特 征分析与活动层厚度反演. 《中国优秀硕士学位 论文全文数据库基 础科学辑》 .2021,(第09期),第A0 08-64页.2/2 页 2[接上页] CN 115291215 B1.一种基于升降轨SAR卫星的长时序二维形变快速解算方法， 其特征在于， 包括如下步 骤： S1、 采用时序I nSAR方法解 算升降轨卫星SAR数据， 得到 升降轨视线向时序形变图集； S2、 比较升降轨时序形变图四个角点 坐标大小， 确定公共区域的坐标 范围； S3、 根据公共范围坐标分别提取 范围约束后的升降轨时序形变图集； S4、 选择升 轨或降轨任意 一个时序形变图集按照时间 间隔为1天进行三次样条插值； S5、 由另一未选中进行插值的时序形变图集， 得到其存在累积形变的时间点， 并联合插 值后的时序形变图集 提取升降轨公共时间点及其对应的时序形变图集； S6、 将时间跨度统一后的升降轨时序形变 图集分别与时序上第一幅形变图做差， 得到 以第一公共时间点 为起始数据的时序形变图集； S7、 分别采用经典二维形变解算模型进行解算， 得到最终具有统一时间跨度基准的时 序二维形变。 2.根据权利要求1所述的一种基于升降轨SAR卫星的长时序二维形变快速解算方法， 其 特征在于， 所述 步骤S1包括： S11、 在雷达影 像坐标系中分别对升 轨和降轨SAR影 像数据集进行差分干涉处 理； S12、 采用相位解缠技术对升降轨差分干涉图进行相位解缠， 获得解缠后的升降轨差分 干涉形变相位图集； S13、 然后分别 对升降轨差分干涉形变相位图集采用时序InSA R技术进行时序累积形变 量解算， 并通过地理编码得到地理坐标系中升降轨 视线向时序形变图集。 3.根据权利要求1所述的一种基于升降轨SAR卫星的长时序二维形变快速解算方法， 其 特征在于， 所述 步骤S2包括： S21、 分别 对升轨和降轨视线向时序形变图东西南北四个角点的经纬度坐标进行比较， 确定升降轨形变图四个角点 坐标之间的大小关系； S22、 通过截取西北角点经度最大值、 纬度最小值， 东北角点经度最小值、 纬度最小值， 东南角点经度最小值、 纬度最大值， 西南角点经度最大值、 纬度最大值， 进而确定升降轨时 序形变图集的公共区域范围。 4.根据权利要求1所述的一种基于升降轨SAR卫星的长时序二维形变快速解算方法， 其 特征在于， 所述 步骤S3包括： S31、 根据公共区域范围分别对所有升降轨 视线向时序形变图集进行裁 剪； S32、 分别存 储裁剪为相同区域的升降轨 视线向时序形变图集。 5.根据权利要求1所述的一种基于升降轨SAR卫星的长时序二维形变快速解算方法， 其 特征在于， 所述 步骤S4包括： S41、 在时间域分别提取升 轨或降轨时序形变图集中的每 个像素点的时序形变数据集； S42、 采用三次样条插值的方法对每 个像素点时序形变按照1 天时间间隔进行插值； S43、 然后按照以上方法依次对时序形变图集中的每个像素点进行插值计算， 进而获得 经三次样条插值后时间 间隔为1天的时序形变图集。 6.根据权利要求1所述的一种基于升降轨SAR卫星的长时序二维形变快速解算方法， 其 特征在于， 所述 步骤S5包括： S51、 提取另一未进行插值计算的轨道影 像时间点；权　利　要　求　书 1/2 页 2 CN 115291215 B 3