融合SAR与光学偏移量技术的地表维形变计算方法及系统.pdf
《融合SAR与光学偏移量技术的地表维形变计算方法及系统.pdf》由会员分享,可在线阅读,更多相关《融合SAR与光学偏移量技术的地表维形变计算方法及系统.pdf(15页完成版)》请在专利查询网上搜索。
1、(19)中华人民共和国国家知识产权局 (12)发明专利申请 (10)申请公布号 (43)申请公布日 (21)申请号 201911324493.3 (22)申请日 2019.12.20 (71)申请人 长安大学 地址 710064 陕西省西安市雁塔区二环南 路中段126号 (72)发明人 赵超英陈立权刘晓杰陈恒祎 王宝行丁明涛 (74)专利代理机构 西安恒泰知识产权代理事务 所 61216 代理人 李婷金艳婷 (51)Int.Cl. G01S 13/90(2006.01) G01S 13/86(2006.01) G01B 11/16(2006.01) G01B 21/32(2006.01) (5。
2、4)发明名称 融合SAR与光学偏移量技术的地表维形变计 算方法及系统 (57)摘要 本发明公开了融合SAR与光学偏移量技术的 地表三维形变计算方法及系统, 包括: 获取待监 测区域形变前后的两景SAR数据, 利用SAR偏移量 跟踪方法根据两景单视复数影像计算待监测区 域地表的方位向形变DAZI和视线向形变DLOS; 获取 待监测区域形变前后的两景光学, 根据两景光学 数据计算地表的南北向形变DNS和东西向形变 DEW; 最后通过公式(1), 获得地表的真实垂直向 形变Dup、 真实南北向形变Dns和真实东西向形变 Dew。 本发明方法和系统能够获得地表三维真实 形变场, 为诸如地震、 滑坡等地。
3、质灾害的研究提 供重要的技术支撑。 权利要求书2页 说明书8页 附图4页 CN 111077525 A 2020.04.28 CN 111077525 A 1.融合SAR与光学偏移量技术的地表三维形变计算方法, 其特征在于, 包括以下步骤: 步骤1, 获取待监测区域形变前后的两景SAR数据, 将两景SAR数据处理为两景单视复数 影像, 根据两景单视复数影像计算待监测区域地表的方位向形变DAZI和视线向形变DLOS; 步骤2, 获取待监测区域形变前后的两景光学数据, 两景光学数据与两景SAR数据覆盖 区域相同, 根据两景光学数据计算地表的南北向形变DNS和东西向形变DEW; 步骤3, 对步骤1获。
4、得的方位向形变DAZI、 视线向形变DLOS和步骤2获得的南北向形变DNS和 东西向形变DEW进行重采样, 使其具有相同的空间分辨率; 步骤4, 通过公式(1), 获得地表的真实垂直向形变Dup、 真实南北向形变Dns和真实东西向 形变Dew; 式中, 表示卫星入射角, 表示卫星飞行方位角。 2.如权利要求1所述的融合SAR与光学偏移量技术的地表三维形变计算方法, 其特征在 于, 所述的步骤1中利用SAR偏移量跟踪方法计算待监测区域地表的方位向形变DAZI和视线 向形变DLOS。 3.如权利要求1所述的融合SAR与光学偏移量技术的地表三维形变计算方法, 其特征在 于, 所述的步骤2包括: 步骤。
5、2.1, 对于获取的两景光学数据, 判断两景光学数据的覆盖范围大小是否相同, 若 相同, 进行步骤2.2; 若不同, 裁剪出两景数据相同的公共区域, 对两景公共区域进行步骤 2.2; 步骤2.2, 使用基于ENVI软件的COSI-Corr光学偏移量计算工具获得地表的南北向形变 量DNS和东西向形变量DEW。 4.如权利要求1所述的融合SAR与光学偏移量技术的地表三维形变计算方法, 其特征在 于, 所述的步骤2还包括对南北向形变DNS和东西向形变DEW的去除误差过程: 根据得到的南 北向形变DNS和东西向形变DEW掩模形变区域, 在远离形变区域的位置选取若干像素点, 根据 这些像素点的形变值及地。
6、理位置利用一次拟合多项式拟合整个监测区域的轨道误差趋势 面, 然后将原始形变图与此误差趋势面相减, 得到去除误差后的南北向形变DNS和东西向形 变DEW。 5.如权利要求1所述的融合SAR与光学偏移量技术的地表三维形变计算方法, 其特征在 于, 所述的步骤4中, 在使用公式(1)计算Dup、 Dns和Dew时, 采用赫尔默特方差分量估计的方法 迭代来最终确定SAR与光学两类形变观测值的权, 进而求得地表的真实垂直向形变量Dup、 真 实南北向形变量Dns和真实东西向形变量Dew。 6.融合SAR与光学偏移量技术的地表三维形变计算系统, 其特征在于, 该系统包括: SAR偏移量计算模块, 用于获。
7、取待监测区域形变前后的两景SAR数据, 将两景SAR数据处 理为两景单视复数影像, 根据两景单视复数影像, 利用SAR偏移量跟踪方法计算待监测区域 地表的方位向形变DAZI和视线向形变DLOS; 权利要求书 1/2 页 2 CN 111077525 A 2 光学偏移量计算模块, 用于获取待监测区域形变前后的两景光学数据, 两景光学数据 与SAR偏移量计算模块中两景SAR数据覆盖区域相同, 根据两景光学数据计算地表的南北向 形变DNS和东西向形变DEW; 重采样模块, 用于对SAR偏移量计算模块获得的方位向形变DAZI、 视线向形变DLOS和光学 偏移量计算模块获得的南北向形变DNS和东西向形变。
8、DEW进行重采样, 使其具有相同的空间 分辨率; 地表三维形变计算模块, 用于通过公式(1), 获得地表的真实垂直向形变Dup、 真实南北 向形变Dns和真实东西向形变Dew; 式中, 表示卫星入射角, 表示卫星飞行方位角。 7.如权利要求6所述的融合SAR与光学偏移量技术的地表三维形变计算系统, 其特征在 于, 所述的SAR偏移量计算模块中利用SAR偏移量跟踪方法计算待监测区域地表的方位向形 变DAZI和视线向形变DLOS。 8.如权利要求6所述的融合SAR与光学偏移量技术的地表三维形变计算系统, 其特征在 于, 所述的光学偏移量计算模块包括: 对于获取的两景光学数据, 判断两景光学数据的覆。
9、盖范围大小是否相同, 若相同, 使用 基于ENVI软件的COSI-Corr光学偏移量计算工具获得地表的南北向形变量DNS和东西向形变 量DEW; 若不同, 裁剪出两景数据相同的公共区域, 再使用基于ENVI软件的COSI-Corr光学偏 移量计算工具生成地表的南北向形变量DNS和东西向形变量DEW。 9.如权利要求6所述的融合SAR与光学偏移量技术的地表三维形变计算系统, 其特征在 于, 所述的光学偏移量计算模块还包括误差去除模块, 用于去除南北向形变DNS和东西向形 变DEW中的误差, 具体为: 根据得到的南北向形变DNS和东西向形变DEW掩模形变区域, 在非形 变区域的位置选取若干像素点,。
10、 根据这些像素点的形变值及地理位置利用一次拟合多项式 模拟整个待监测区域的轨道误差趋势面, 然后将原始形变图与此误差趋势面相减, 得到的 去除误差后的南北向形变DNS和东西向形变DEW。 10.如权利要求6所述的融合SAR与光学偏移量技术的地表三维形变计算系统, 其特征 在于, 所述的地表三维形变计算模块中, 在使用公式(1)计算Dup、 Dns和Dew时, 采用赫尔默特 方差分量估计的方法迭代来最终确定SAR与光学两类形变观测值的权, 进而求得地表的真 实垂直向形变量Dup、 真实南北向形变量Dns和真实东西向形变量Dew。 权利要求书 2/2 页 3 CN 111077525 A 3 融合。
11、SAR与光学偏移量技术的地表维形变计算方法及系统 技术领域 0001 本发明属于地表三维形变监测技术领域, 具体涉及一种融合SAR与光学偏移量技 术的地表三维形变计算方法及系统, 可应用于地震、 滑坡、 火山等地质灾害监测领域的研究 中。 背景技术 0002 目前, 遥感技术能够有效应用于针对诸如地震、 滑坡、 火山、 地裂缝等地质灾害的 调查、 监测及预警工作中。 其中基于SAR数据和光学数据的两类技术作为遥感领域重要的两 个分支, 在地质灾害研究中均充当了重要的角色。 0003 针对光学的应用主要通过目视解译或多个波段数据的组合变换进而提取表征地 质灾害体特征的因素, 此类方法皆是对灾害体。
12、展开定性的研究, 利用光学偏移量技术可以 实现地质灾害所引起地表形变的定量分析, 然而传统的方法往往只能获取地表二维形变。 0004 基于SAR数据的地表形变监测技术主要包括合成孔径雷达干涉测量(InSAR)技术 及SAR偏移量技术等, 目前利用SAR数据的合成孔径雷达干涉测量(InSAR)及偏移量技术能 够有效的应用于地表三维形变场的恢复, 但为了获取研究区域地表三维形变, 传统的方法 至少需要三个来自于不同雷达成像几何结构的SAR数据。 当可用SAR卫星数据源较少, 造成 该方法在实际应用中不足以满足三维形变场计算条件时, 这些技术便无法开展。 发明内容 0005 针对现有技术的缺陷和不足。
13、, 本发明给出了一种融合SAR与光学偏移量技术的地 表三维形变计算方法及系统, 解决现有的地表三维形变解算在SAR数据源不足时使用受限 的问题。 0006 为了实现上述目的, 本发明采用如下技术方案予以实现: 0007 融合SAR与光学偏移量技术的地表三维形变计算方法, 包括以下步骤: 0008 步骤1, 获取待监测区域形变前后的两景SAR数据, 将两景SAR数据处理为两景单视 复数影像, 根据两景单视复数影像计算待监测区域地表的方位向形变DAZI和视线向形变 DLOS; 0009 步骤2, 获取待监测区域形变前后的两景光学数据, 两景光学数据与两景SAR数据 覆盖区域相同, 根据两景光学数据。
14、计算地表的南北向形变DNS和东西向形变DEW; 0010 步骤3, 对步骤1获得的方位向形变DAZI、 视线向形变DLOS和步骤2获得的南北向形变 DNS和东西向形变DEW进行重采样, 使其具有相同的空间分辨率; 0011 步骤4, 通过公式(1), 获得地表的真实垂直向形变Dup、 真实南北向形变Dns和真实 东西向形变Dew; 说明书 1/8 页 4 CN 111077525 A 4 0012 0013式中, 表示卫星入射角, 表示卫星飞行方位角。 0014 具体的, 所述的步骤1中利用SAR偏移量跟踪方法计算待监测区域地表的方位向形 变DAZI和视线向形变DLOS。 0015 具体的, 。
15、所述的步骤2包括: 0016 步骤2.1, 对于获取的两景光学数据, 判断两景光学数据的覆盖范围大小是否相 同, 若相同, 进行步骤2.2; 若不同, 裁剪出两景数据相同的公共区域, 对两景公共区域进行 步骤2.2; 0017 步骤2.2, 使用基于ENVI软件的COSI-Corr光学偏移量计算工具获得地表的南北向 形变量DNS和东西向形变量DEW。 0018 优选的, 所述的步骤2还包括对南北向形变DNS和东西向形变DEW的去除误差过程: 根据得到的南北向形变DNS和东西向形变DEW掩模形变区域, 在远离形变区域的位置选取若 干像素点, 根据这些像素点的形变值及地理位置利用一次拟合多项式拟合。
16、整个监测区域的 轨道误差趋势面, 然后将原始形变图与此误差趋势面相减, 得到去除误差后的南北向形变 DNS和东西向形变DEW。 0019 优选的, 所述的步骤4中, 在使用公式(1)计算Dup、 Dns和Dew时, 采用赫尔默特方差分 量估计的方法迭代来最终确定SAR与光学两类形变观测值的权, 进而求得地表的真实垂直 向形变量Dup、 真实南北向形变量Dns和真实东西向形变量Dew。 0020 本发明还公开了融合SAR与光学偏移量技术的地表三维形变计算系统, 该系统包 括: 0021 SAR偏移量计算模块, 用于获取待监测区域形变前后的两景SAR数据, 将两景SAR数 据处理为两景单视复数影像。
17、, 根据两景单视复数影像, 利用SAR偏移量跟踪方法计算待监测 区域地表的方位向形变DAZI和视线向形变DLOS; 0022 光学偏移量计算模块, 用于获取待监测区域形变前后的两景光学数据, 两景光学 数据与SAR偏移量计算模块中两景SAR数据覆盖区域相同, 根据两景光学数据计算地表的南 北向形变DNS和东西向形变DEW; 0023 重采样模块, 用于对SAR偏移量计算模块获得的方位向形变DAZI、 视线向形变DLOS和 光学偏移量计算模块获得的南北向形变DNS和东西向形变DEW进行重采样, 使其具有相同的 空间分辨率; 0024 地表三维形变计算模块, 用于通过公式(1), 获得地表的真实垂。
18、直向形变Dup、 真实 南北向形变Dns和真实东西向形变Dew; 0025 说明书 2/8 页 5 CN 111077525 A 5 0026式中, 表示卫星入射角, 表示卫星飞行方位角。 0027 具体的, 所述的SAR偏移量计算模块中利用SAR偏移量跟踪方法计算待监测区域地 表的方位向形变DAZI和视线向形变DLOS。 0028 具体的, 所述的光学偏移量计算模块包括: 0029 对于获取的两景光学数据, 判断两景光学数据的覆盖范围大小是否相同, 若相同, 使用基于ENVI软件的COSI-Corr光学偏移量计算工具获得地表的南北向形变量DNS和东西向 形变量DEW; 若不同, 裁剪出两景数。
19、据相同的公共区域, 再使用基于ENVI软件的COSI-Corr光 学偏移量计算工具生成地表的南北向形变量DNS和东西向形变量DEW。 0030 进一步的, 所述的光学偏移量计算模块还包括误差去除模块, 用于去除南北向形 变DNS和东西向形变DEW中的误差, 具体为: 根据得到的南北向形变DNS和东西向形变DEW掩模形 变区域, 在非形变区域的位置选取若干像素点, 根据这些像素点的形变值及地理位置利用 一次拟合多项式模拟整个待监测区域的轨道误差趋势面, 然后将原始形变图与此误差趋势 面相减, 得到的去除误差后的南北向形变DNS和东西向形变DEW。 0031 优选的, 所述的地表三维形变计算模块中。
20、, 在使用公式(1)计算Dup、 Dns和Dew时, 采 用赫尔默特方差分量估计的方法迭代来最终确定SAR与光学两类形变观测值的权, 进而求 得地表的真实垂直向形变量Dup、 真实南北向形变量Dns和真实东西向形变量Dew。 0032 与现有技术相比, 本发明的有益效果是: 0033 本发明通过融合SAR与光学偏移量技术得到一种用于计算地表三维形变的方法, 该方法的优点在于利用SAR与光学数据, 通过偏移量技术获得地表形变, 然后通过建立地表 真实三维形变求解模型, 且通过赫尔默特方差分量估计的方法最终确定SAR与光学两类形 变观测值的权重, 最终能够有效地获得地表三维真实形变, 解决了现有仅。
21、依赖SAR数据的地 表三维形变解算在SAR数据源不足时无法恢复地表三维形变的问题, 为地表三维形变计算 提供了一种新的方法, 此工作对于研究地表形变具有重要意义。 附图说明 0034 图1为本发明实施例记载的地表三维形变计算方法流程图。 0035 图2为实施例记载的凯库拉地震沿方位向形变监测结果。 0036 图3为实施例记载的凯库拉地震沿视线向形变监测结果。 0037 图4为实施例记载的凯库拉地震在南北向上的形变监测结果, 图中沿北方向形变 为正, 沿南方向形变为负。 0038 图5为实施例记载的凯库拉地震在东西向上的形变监测结果, 图中沿东方向形变 为正, 沿西方向形变为负。 0039 图6。
22、为利用本发明所得凯库拉地震地表三维形变场; 图中黑色箭头表示地表水平 向形变, 箭头方向为水平形变方向, 长度表示形变大小, 颜色变化表示垂直向形变。 0040 以下结合说明书附图和具体实施方式对本发明做具体说明。 具体实施方式 0041 以下给出本发明的具体实施方式, 需要说明的是本发明并不局限于以下具体实施 例中, 凡在本申请技术方案基础上做的等同变换均落入本发明的保护范围。 说明书 3/8 页 6 CN 111077525 A 6 0042 如图1所示, 本发明实施例提供的一种融合SAR与光学偏移量技术的地表三维形变 计算方法包括以下步骤: 0043 步骤1, 获取待监测区域的两景SAR。
23、数据, 如地震前后同一区域的两景数据; 将两景 SAR数据处理为两景单视复数影像, 利用SAR偏移量跟踪方法根据两景单视复数影像计算待 监测区域地表的方位向形变DAZI和视线向形变DLOS。 0044 SAR偏移量跟踪方法的具体思路是: 选定其中一景作为主影像, 将另一景从影像配 准重采样到该主影像, 使主影像和从影像具有相同的空间几何结构; 利用SAR幅度信息对两 景SAR影像进行粗配准, 计算主、 从影像之间的初始偏移量, 然后在粗配准的基础上, 选取一 定大小的搜索窗口并计算主从影像之间的归一化互相关系数 (x,y), 如式(2)所示, 0045 0046 式中, x和y为主影像坐标, 。
24、xs和ys为两个方向上的偏移量, m和n为窗口大小, fm和fs 分别为主影像幅度和从影像幅度。 0047 当互相关系数 (x,y)达到最大时便可获得精确的子像素级的偏移量, 最终获得地 表的方位向形变DAZI和视线向形变DLOS。 0048 步骤2, 获取待监测区域形变前后的两景光学数据, 两景光学数据与步骤1中两景 SAR数据获取时间点大致相同、 覆盖区域相同。 根据两景光学数据计算地表的光学偏移量, 即南北向形变DNS和东西向形变DEW; 0049 需要注意的是, 在利用两景光学数据计算地表的光学偏移量时, 首先判断两景光 学数据的覆盖范围大小是否相同, 若相同, 则使用基于ENVI软件。
25、的COSI-Corr光学偏移量计 算工具生成地表的南北向形变量DNS和东西向形变量DEW; 若不同, 需要先裁剪出两景数据相 同的公共区域, 对两景公共区域再使用基于ENVI软件的COSI-Corr光学偏移量计算工具生 成地表的南北向形变量DNS和东西向形变量DEW。 0050 对于得到的南北向形变量DNS和东西向形变量DEW, 由于一些光学影像数据没有经 过严格的正射校正和几何校正, 所以在计算所得形变场中存在明显的轨道误差, 因此需要 去除该误差, 具体为: 根据得到的南北向形变DNS和东西向形变DEW掩模形变区域, 在非形变 区域的位置选取若干像素点, 根据这些像素点的形变值及地理位置利。
26、用一次拟合多项式拟 合整个待监测区域的轨道误差趋势面, 然后将原始形变图与此误差趋势面相减, 得到去除 误差后的南北向形变DNS和东西向形变DEW。 0051 步骤3, 对步骤1获得的方位向形变DAZI、 视线向形变DLOS和步骤2获得的南北向形变 DNS和东西向形变DEW进行重采样, 使其具有相同的空间分辨率, 并裁剪为相同的范围。 0052 另外, 若光学偏移量结果和SAR偏移量结果的坐标系不一致时, 需要将两者转化为 相同的坐标系, 如光学偏移量结果通常为UTM坐标系, 而经过地理编码后的SAR偏移量结果 往往为WGS 84地理坐标系, 因此首先需要将光学偏移量结果由UTM坐标系转为与S。
27、AR结果相 同的地理坐标系。 0053 步骤4, 根据SAR卫星成像几何关系建立方位向形变和视线向形变与地表真实三维 形变之间的关系式(3)和(4), 0054 说明书 4/8 页 7 CN 111077525 A 7 0055 0056 式中, Dup表示地表的真实垂直向形变, Dns表示真实南北向形变, Dew表示真实东西 向形变; 0057 再结合光学所得南北向和东西向形变, 可建立SAR与光学偏移量联合求解得到地 表三维形变的观测方程, 如式(1)所示。 0058 0059式中, 表示卫星入射角, 表示卫星飞行方位角。 0060 最后通过公式(1), 获得地表的真实垂直向形变Dup、 。
28、真实南北向形变Dns和真实东 西向形变Dew。 0061 在使用公式(1)计算Dup、 Dns和Dew时, 采用赫尔默特方差分量估计的方法通过迭代 定权, 最终确定SAR与光学形变观测值在计算过程中的权重, 并最终求得地表三维形变场 Dup、 Dns和Dew。 具体步骤如下: 0062 步骤1), 首先确定初始方差, 通过给定一个遍历窗口, 分别计算方位向、 视线向、 东 西向及南北向形变量文件中每个中心像元在该窗口内的方差作为初始方差, 利用式(5)权 计算公式计算初始权, 式中c为任意给定常数, i2为第i类观测方差, pi即为权值。 令计算所 得SAR与光学两类形变观测值的权分别为p1和。
29、p2; 0063 0064步骤2), 令式(4)中L1DLOS DAZIT, L2DNS DEWT,A 为系数矩阵, AA1 A2T, 其中, 0065 0066 0067 将式(4)改写为误差方程,如式(8)所示, 0068 0069 由式(8)可得法方程为 0070 0071式中, 0072步骤3), 最终建立残差平方和和之间的关系式如式(10)所示, 其中, S矩阵如式(11)所示, 说明书 5/8 页 8 CN 111077525 A 8 0073 0074 0075 根据公式(10)求得两类观测值迭代后单位权方差, 计算公式如式(12)所示, 0076 0077 步骤4), 根据所求。
30、得两类观测值迭代后单位权方差, 利用公式(13)再次定权, 式 中, k为迭代次数, 0078 0079 然后重复步骤2)和步骤3), 进行迭代, 直至两类观测值迭代后的单位权方差相等 迭代停止, 最终求得即求得地表三维形变。 0080 本发明的实施例还提供了一种融合SAR与光学偏移量技术的地表三维形变计算系 统, 该系统包括: 0081 SAR偏移量计算模块, 用于获取待监测区域形变发生前后的两景SAR数据, 将两景 SAR数据处理为两景单视复数影像, 利用SAR偏移量跟踪方法根据两景单视复数影像中计算 待监测区域地表的方位向形变DAZI和视线向形变DLOS。 SAR偏移量跟踪方法的具体思路。
31、同上 文记载。 0082 光学偏移量计算模块, 用于获取待监测区域形变发生前后的两景光学数据, 两景 光学数据与SAR偏移量计算模块中两景SAR数据覆盖区域相同, 根据两景光学数据, 利用光 学偏移量技术计算地表的南北向形变DNS和东西向形变DEW; 0083 需要注意的是, 在利用两景光学数据计算地表的光学偏移量时, 首先判断两景光 学数据的覆盖范围大小是否相同, 若相同, 则使用基于ENVI软件的COSI-Corr光学偏移量计 算工具生成地表的南北向形变量DNS和东西向形变量DEW; 若不同, 需要先裁剪出两景数据的 公共区域, 对两景公共区域再使用基于ENVI软件的COSI-Corr光学。
32、偏移量计算工具生成地 表的南北向形变量DNS和东西向形变量DEW。 0084 对于得到的南北向形变量DNS和东西向形变量DEW, 由于一些光学影像数据没有经 过严格的正射校正和几何校正, 所以在计算所得形变场中存在明显的轨道误差, 因此还设 置了误差去除模块, 误差去除模块用于通过掩模得到的南北向形变DNS和东西向形变DEW形 变区域, 在非形变区域的位置选取若干像素点, 根据这些像素点的形变值及地理位置利用 一次拟合多项式拟合整个待监测区域的轨道误差趋势面, 然后将原始形变图与此误差趋势 面相减, 得到去除误差后的南北向形变DNS和东西向形变DEW。 0085 重采样模块, 用于对SAR偏移。
33、量计算模块获得的方位向形变DAZI、 视线向形变DLOS和 光学偏移量计算模块获得的南北向形变DNS和东西向形变DEW进行重采样, 使其具有相同的 空间分辨率; 0086 另外, 若光学偏移量结果和SAR偏移量结果的坐标系不一致时, 需要设置一个坐标 系转化模块, 该模块用于将两个不同的坐标系转化为相同的坐标系, 如光学偏移量结果通 常为UTM坐标系, 而经过地理编码后的SAR偏移量结果往往为WGS 84地理坐标系, 因此首先 需要将光学偏移量结果由UTM坐标系转为与SAR结果相同的地理坐标系。 说明书 6/8 页 9 CN 111077525 A 9 0087 地表三维形变计算模块, 用于通。
34、过公式(1), 获得地表的真实垂直向形变Dup、 真实 南北向形变Dns和真实东西向形变Dew; 0088 0089式中, 表示卫星入射角, 表示卫星飞行方位角。 0090 在地表三维形变计算模块中, 使用公式(1)计算Dup、 Dns和Dew时, 采用赫尔默特方差 分量估计的方法通过迭代定权, 最终确定SAR与光学形变观测值在计算过程中的权重, 并最 终求得地表三维形变场Dup、 Dns和Dew。 具体步骤同上文中的步骤1)至步骤4)。 0091 实施例 0092 本实施例选择位于新西兰东北部的凯库拉县地区作为实验区域, 该地区于2016年 11月14日发生了MW7.8级大地震, 地震发生之。
35、后, 地表破裂明显, 且形变量级较大。 实验选用 的SAR数据为2016年2月20日和2016年12月10日覆盖该地区的两景地震发生前后Sentinel- 1数据, 光学数据为2016年2月19日和2017年2月23日覆盖该地区的两景地震发生前后 Sentinel-2数据, 由于本实施例仅研究该地震同震三维形变的计算, 且同震形变量级远远 大于震间形变, 因此震后的SAR与光学数据获取时间节点相差时间内地表发生的震间形变 可忽略不计, 即可视为数据获取时间大致相同。 0093 实验过程 0094 步骤1, 首先对于获取的两景Sentinel-1数据进行预处理生成单视复数影像 (SLC), 然后。
36、选定2016年2月20日影像作为主影像, 将另一景从影像配准重采样到该主影像, 使主从影像具有相同的空间几何结构。 然后利用SAR偏移量技术, 设置搜索窗口为64*64, 搜 索步长为8*8, 计算方位向和视线向二维形变场, 分别如图2和图3所示。 图中Azimuth为卫星 飞行方向, 即为方位向, LOS为沿雷达视线方向, 即为视线向, 远离卫星方向形变量为负, 靠 近卫星方向形变量为正。 0095 步骤2, 获取与SAR数据覆盖相同地区的2016年2月19日和2017年2月23日两景 Sentinel-2数据, 首先进行数据预处理, 包括辐射定标和大气校正, 最终提取第8波段进行 光学偏移。
37、量的计算。 由于数据覆盖范围大小不同, 因此需裁剪出两景数据的公共区域。 然后 利用基于ENVI软件的COSI-Corr光学偏移量计算工具计算光学偏移量, 设置搜索窗口大小 为32*32, 移动步长为8*8, 掩膜阈值为0.9, 迭代次数为2次, 计算获得地表在南北向和东西 向上的形变。 由于一些光学影像数据没有经过严格的正射校正和几何校正, 所以在计算所 得形变场中存在明显的轨道误差, 因此在非形变区域的位置选取若干像素点, 根据这些像 素点的形变值及地理位置利用一次拟合多项式模拟整个研究区域的轨道误差趋势面, 然后 将原始形变图与此误差趋势面相减, 达到去除轨道误差的目的, 最终获得地表在。
38、南北向和 东西向上的形变如图4和图5所示。 0096 步骤3, 对SAR与光学偏移量结果进行重采样, 使二者具有相同的空间分辨率和坐 标系统。 由于Sentinel-2光学偏移量的结果为UTM坐标系, 因此首先需要将光学偏移量结果 由UTM坐标系转为WGS 84下的地理坐标系, 然后将SAR与光学结果重采样为相同的空间分辨 率, 并且裁剪为相同的范围, 最终导出为TIFF格式文件。 说明书 7/8 页 10 CN 111077525 A 10 0097 步骤4, 根据式(1), 通过步骤1)至步骤4)获得地表的真实垂直向形变Dup、 真实南北 向形变Dns和真实东西向形变Dew, 如图6所示为地表的三维形变场。 由图6中矢量箭头的方向 及垂直方向的形变可以看出地震造成各板块之间相互作用, 引起各板块的地表发生不同的 形变, 从而能够有效地解释地震所造成板块运动的机理。 说明书 8/8 页 11 CN 111077525 A 11 图1 说明书附图 1/4 页 12 CN 111077525 A 12 图2 图3 说明书附图 2/4 页 13 CN 111077525 A 13 图4 图5 说明书附图 3/4 页 14 CN 111077525 A 14 图6 说明书附图 4/4 页 15 CN 111077525 A 15 。
- 内容关键字: 融合 SAR 光学 偏移 技术 地表 形变 计算方法 系统
变速箱.pdf
再生铅烟气中铅尘环保治理湿电除尘超低排放设备.pdf
混合加料机构.pdf
漆包线脱漆机.pdf
灯具的透镜反光杯一体式结构.pdf
防坠网-栏一体装置.pdf
装修用管线固定结构.pdf
物料压平测高装置.pdf
氢溴酸山莨菪碱分析溶解仪器.pdf
粉料分散上料机.pdf
电磁屏蔽壳体以及新能源汽车.pdf
枫叶魔方.pdf
骨科助行装置.pdf
钣金件调直装置.pdf
节能建筑保温体系.pdf
超薄柔性关节镜通道.pdf
适用于有水状态下闸门井封堵后干地修复的施工方法.pdf
车辆充电桩用防尘防水枪头结构.pdf
棉浆粕生产用的蒸球.pdf
码垛机器人几何精度标定与辨识方法.pdf
钕铁硼自动化发蓝设备.pdf
钢管外涂塑粉末自清洁回收装置.pdf
货物出库管理方法、系统及介质.pdf
词汇标签的确定方法及系统和查询语句的生成方法及系统.pdf
用于电机定子的端部叠片和电机定子.pdf
船舶目标检测方法、系统、可读存储介质及计算机.pdf
具有圆柱螺旋形液流流道的磁流变阻尼器.pdf
基于AES算法的硬件安全方法、设备及存储介质.pdf
接口管理方法及接口管理系统.pdf
涂料的质量监测方法及系统.pdf
空气灰尘检测仪.pdf
高铬合金离心铸管机及其使用方法.pdf
机械成型鼓及其无缝锁定方法.pdf
一种新型连续结晶与提纯装置.pdf
尼伐地平新晶型及其制备方法.pdf
高强抗震钢筋及其制造方法.pdf
治疗与年龄相关的黄斑变性AMD的方法.pdf
用于操作具有液体喷射装置的干衣机的方法.pdf
三角形滤芯过滤器.pdf
含有聚缩醛的刷.pdf
香豆素35的制备方法.pdf
一种节能电动汽车充电器温度控制装置.pdf
高低压式缓冲系统.pdf
聚酰胺与聚烯烃的深染工艺.pdf
齿轮花键磨床.pdf
硅棒切割机双V形硅棒支承座.pdf
硅棒切割夹具的同步夹紧机构.pdf
一种压滤板.pdf
丁苯那嗪的拆分方法.pdf
挠性管.pdf
光涂料与其制备方法以及利用该涂料制得的感光膜.pdf