买专利卖专利找龙图腾，真高效！ 查专利查商标用IPTOP,全免费！专利年费监控用IP管家,真方便！

申请/专利权人：中国科学院生态环境研究中心

摘要：本发明提供了一种监测山体现今抬升速度的雷达时序差分干涉测量方法，包括：第一步骤：主被动遥感数据协同处理；第二步骤：执行雷达数据差分干涉处理；第三步骤：提取永久散射体的累积形变；第四步骤：针对提取的永久散射体的累积形变执行多尺度滤波；第五步骤：在多尺度滤波后执行山脊范围提取；第六步骤：根据滤波后永久散射体在时间序列上的绝对高程变化反演山体抬升形变。本发明提供了一种能够在宏观尺度上精准监测山体现今抬升速度的雷达时序差分干涉测量方法。

主权项：1.一种监测山体现今抬升速度的雷达时序差分干涉测量方法，其特征在于包括：第一步骤：主被动遥感数据协同处理；第二步骤：执行雷达数据差分干涉处理；第三步骤：提取永久散射体的累积形变；第四步骤：针对提取的永久散射体的累积形变执行多尺度滤波；其中，第四步骤对永久散射体提取结果执行多尺度滤波，在空间域滤除稳定地表、植被变化及人类扰动引起的地表形变，在时间域滤除气候变化引起的地表形变，使监测对象锚定于地体抬升形变；其中，多尺度滤波分为空间域滤波和时间域滤波；空间域滤波采用相位标准差低通滤波滤除植被影响，采用监督式高通滤波滤除稳定地表，采用监督式低通滤波滤除新建设施，采用土地利用分类掩膜人类活动区域；时间域滤波为永久散射体在时间序列上的畸变滤波以滤除突发性地表形变，以及结合物候特征的时间域非线性滤波，滤除季节性地表形变；其中，时间域非线性滤波公式为：Pyi∣a,xi,σ2∽Naxi+b,σ2；式中P为地体累积抬升量，N为地体累积抬升条件概率，a为观测频数，b为地体起始形变值，yi表示第i个地体累积抬升量，xi表示第i个地体绝对高程观测值，σ2为随机扰动项的方差，假设随机扰动项εi服从期望为0，方差为σ2的正态分布；令Y表示地体累积抬升量y1～yi的列向量，ε表示残差的列向量，X为观测值矩阵，β为a的列向量，则区域地体抬升形变的矩阵表达形式如公式：Y＝Xβ+ε；第五步骤：在多尺度滤波后执行山脊范围提取；其中，在多尺度滤波后执行山脊范围提取包括：基于DEM数据使用分水岭方法圈定山脊范围，依次进行山谷填洼、流向计算、流量计算、累积流量分割、分水岭处理，实现流域范围圈定，进而提取分水岭，随后通过设置累积流量分割阈值，对流域范围进行掩膜，将监测窗口锚定于山脊范围；第六步骤：根据滤波后永久散射体在时间序列上的绝对高程变化反演山体抬升形变。

全文数据：

权利要求：

百度查询： 中国科学院生态环境研究中心 一种监测山体现今抬升速度的雷达时序差分干涉测量方法