InSAR 数据处理误差的减弱措施

目录

  • 1.失相干误差
  • 2.基线误差
  • 3.DEM 误差
  • 4.大气误差
  • 5.解缠误差
  • 6.地理编码误差

在这里插入图片描述

本文由CSDN点云侠原创,爬虫网站请自重。

  InSAR 获取的干涉相位通常可表示为:
φ i n t ( η , ξ ) = φ d e f ( η , ξ ) + φ a t m ( η , ξ ) + φ t o p o ( η , ξ ) + φ o r b i t ( η , ξ ) + φ n o i s e ( η , ξ ) \varphi^{int}(\eta,\xi)=\varphi^{def}(\eta,\xi)+\varphi^{atm}(\eta,\xi)+\varphi^{topo}(\eta,\xi)+\varphi^{orbit}(\eta,\xi)+\varphi^{noise}(\eta,\xi) φint(η,ξ)=φdef(η,ξ)+φatm(η,ξ)+φtopo(η,ξ)+φorbit(η,ξ)+φnoise(η,ξ)
  其中, φ d e f \varphi^{def} φdef为形变相位, φ a t m \varphi^{atm} φatm为大气延迟误差相位, φ t o p o \varphi^{topo} φtopo为地形误差相位, φ o r b i t \varphi^{orbit} φorbit为轨道误差相位, φ n o i s e \varphi^{noise} φnoise为噪声相位。另外,在获取的 InSAR 相位中,还常常含有解缠误差。为获取高精度的形变监测结果,必须对其中的这些误差项进行有效的估计和剔除。

1.失相干误差

  对重复轨道干涉测量来说,其常常受到很多的失相干因素的影响,这些影响因素会在干涉处理中以偶然误差的形式传播到干涉形变结果中。如两幅影像获取期间随着地表的植被、土壤湿度等的变化,会在地表分辨单元内产生随机性的附加相位延迟变化,即时间失相干变化;另外,重复轨道监测中,偏移的轨道间隔即空间基线又会导致雷达的侧视角度不同,从而使得两次成像时刻会在同一分辨单元内产生不同的附加相位分量,这种影响在空间上也表现为随机性的偶然误差,这就是所谓的几何(基线)失相干。同样的,在沿轨方向即方位向上,因雷达两次成像的照射方向即多普勒中心不一致,也会引起相位的失相干,这就是所谓的多普勒中心失相关。这些随机性的失相干误差常常会引起相位信号的严重失相关,从而会给后续的相位解缠操作带来极大的困难,甚至还可能导致干涉分析无法进行下去。通常可通过滤波的方法来消除这种随机性的失相干误差对形变结果的影响。

2.基线误差

  基线是 InSAR 干涉测量的关键参数,它既是 SAR 干涉成像的基础,也是导致 InSAR干涉相位失相干的主要因素之一。在 InSAR 数据处理中,常用到的基线有两种:时间基线和空间基线。时间基线指的是用于进行干涉处理的两个 SAR 影像获取的时间间隔,它是导致 InSAR 干涉时间失相干的根源;而空间基线则指成像时刻两个卫星间的空间距离,它是产生 InSAR 干涉相位几何失相干的根源。这里,我们讨论的基线误差主要是针对空间基线。
  在InSAR干涉测量中,干涉基线一般是利用轨道的星历数据和状态矢量确定的。当卫星的轨道数据不准确时,就会导致估计的基线存在误差,其表现在干涉图中常为一些明显的系统性条纹,因此,基线误差对形变结果的影响属于系统误差,这常常使得我们很难将干涉相位中的形变分量、大气延迟分量和残差干涉条纹信息分离开来,从而直接影响了最终获取的形变的精度和可靠性。若要完全有效的去除这种基线误差导致的残差干涉条纹,往往要求精密轨道数据的绝对精度要优于1mm,就当前而言,该精度还远远不能满足。
  另外,基线估计算法的精度还和SAR数据的应用相关。在地震、冰川等差分形变监测应用中,ERS和Envisat精密轨道的精度已可满足要求;但是在地面沉降等高精度的地表形变监测应用中,轨道基线的精度要求达到亚厘米级,而现有的轨道数据是不满足的,因此必须提供精确客观的基线估计值。目前常用的基线估计和精化方法可分为基于卫星轨道状态矢量的估计方法、基于外部信息如GCP的估计方法和基于干涉图本身的估计方法。

3.DEM 误差

  无论是在两轨差分 DInSAR 还是时序 InSAR 数据处理中,辅助的外部 DEM 数据精度对形变结果的误差影响也是不容忽视的。DEM 误差对 InSAR 形变监测结果的影响主要表现在两个方面:一是外部 DEM 数据自身的精度,可通过高程模糊度将误差转换到形变信息中;二是在利用外部 DEM 数据与主影像配准时产生的配准误差影响。垂直基线越小,高程模糊度就越大,则由 DEM 误差导致的形变误差就越小。 另外,这种由 DEM 误差所引起的地形误差对形变结果的影响有时表现为系统误差,有时又表现为偶然误差。其中,随机性的地形误差会导致干涉图的相关系数显著降低,但系统性的地形误差几乎不会影响干涉图的相干性,但它会严重的扭曲地表的真实形变。针对 DEM 误差,其最直接的减弱措施即是选择高精度的 DEM;另外,还可以在干涉图较多的时候,将 DEM 作为一个系统参数估计出来,这是大地测量中常用到的一种含有系统误差的平差方法。

4.大气误差

  对星载 InSAR 干涉测量来说,大气误差比较复杂,它主要是由于雷达信号在大气层中传播时对流层折射率变化和电离层电子含量变化所引起的。其中,电离层延迟误差在空间范围内变化较小,可利用差分消除;而对流层延迟误差则在时间和空间上均会产生较大的变化,甚至其在短时间、小范围内也可能产生较大差异。因此,我们这里主要讨论对流层延迟误差对形变结果的影响。通常,对流层延迟误差在空间域内表现为较低频率的相关性,而在时间域上则呈现为高频特性。因此,大气误差对形变结果的影响有时候表现为系统误差,有时候又表现为偶然误差。其中,系统性的大气误差主要是和地形相关的:地形越高,距离卫星的距离越短,大气影响越小;反之,地形越低,距离卫星的距离越远,大气的影响就越大。因此,我们一般可以通过建立一个与地形相关的模型即采用参数估计的方法对这种系统性的大气误差进行消除。而当大气误差表现为偶然性时,也可以通过前面的干涉图滤波消除,但更多的则是在有多个干涉图的情况下,利用时间域的高通滤波来对其进行分离。
  总的来说,大气误差是 InSAR 处理中很关键的一项影响因素,特别的由于其具有空间上的低频特性,往往很难将其和地表形变相位区分开来,从而造成伪形变相位,这种伪形变相位会严重的影响形变结果的解译。目前,有关大气对流层延迟误差的改正方法可大致分为两类:一类是通过获取区域内的多幅影像,分析信号的统计特性进行大气延迟改正,如干涉图的线性叠加法、永久散射体方法等;另一类是利用外部独立数据辅助进行大气相位的消除。首先借助外部独立气象数据如 GPS、MERIS、MODIS 等计算得到大气延迟量,然后再将其从InSAR 干涉相位图中去除以达到大气相位的延迟改正。

5.解缠误差

  相位解缠是 InSAR 数据处理中的难点和重点,它的准确性将直接影响 InSAR 形变监测结果的精度和可靠性。目前,有关 InSAR 相位解缠算法可大致分为 3 类:第一类是基于路径跟踪的算法。路径跟踪法是基于逐个像素搜索,利用相位数据的局部信息,通过一定的方法得到全局最优解,故属于局部算法。路径跟踪算法能有效的抑制局部相位误差在整个积分区域的传播,可在相干性好的区域获得较精确的解缠结果,保证解缠值的可靠性;但是当干涉图含有较大噪声时,该类算法易造成误差传递或无法解缠的孤立区块。第二类是基于最小范数的算法。相对于路径积分算法,最小范数法可归为全局算法,它是通过在解缠结果中加入某种约束即一定的数学模型,使其与
原始相位达到某种程度的最优逼近,如最小二乘相位解缠,它不需识别残差点,具有较好的稳定性,但是正由于其是“穿过”而不“绕过”残差点,易导致误差传递。第三类是基于最优估计的算法,最优估计法是通过使解缠相位的导数和缠绕相位的导数之差最小化实现的,不再需要进行残差点识别,通常都能较好的抑制误差传递,但是其往往理论相对复杂,计算难度大,如常用的网络流模型等。
   然而,我们知道,InSAR 数据处理中受到许多的失相干因素的影响,如干涉噪声、地形起伏引起的叠掩、阴影以及形变梯度过大等,这些失相干因素往往导致我们获取的干涉相位不连续,相位的周期性和趋势性遭到破坏,从而大大的增加了相位解缠的难度,甚至可能造成错误的相位解缠,直接影响了 InSAR 形变监测结果的精度。特别的,这种错误的相位解缠,更多的则体现为解缠粗差,其表现在空间上常为一片一片的出现。通过干涉图滤波,能在一定的程度上提高干涉图的信噪比,提高相位解缠的精度,但更多的还需要通过多干涉图组合的方法来进行相位解缠误差的剔除。

6.地理编码误差

  地理编码是将 InSAR 获取的雷达 SAR 坐标系下的结果转换到地理坐标系下的过程。地理编码误差主要是由于所用的外部 DEM 数据的平面坐标不准确导致的。其表现为在地理坐标系下,当我们将 InSAR 形变结果与其它基础图层叠加分析时,形变点的位置常常会产生一个系统性的偏移,故其对形变结果的影响属于系统误差,一般可通过添加地面控制点来对其进行校正。

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:/a/115476.html

如若内容造成侵权/违法违规/事实不符,请联系我们进行投诉反馈qq邮箱809451989@qq.com,一经查实,立即删除!

相关文章

TCP 协议

文章目录 协议格式1面向连接:1.1三次握手(建立连接)1.2包序管理1.2四次挥手(断开连接) 2可靠传输:一。保证数据可靠有序的到达对端:确认应答机制超时重传机制 二。提高传输效率:1.提升自身发送数据量滑动窗口机制 rwnd滑动窗口丢包…

对5款驱动软件个人感受

1、基于安全考虑,个人建议优先选用官网可以下载到的版本。 2、拓展: 在线版:普通版,需要安装,且只能联网使用。离线版:把所有驱动都包含在内,无需联网也能安装驱动,一般从大小就可区…

教你怎么用Python每天自动给女朋友免费发短信

今天的教程就是教大家怎么发送免费短信给女朋友。 发送短信接口,我知道的常见的有两个平台,一个是 twilio,可以免费发短信 500 条,可发任意信息,一个是腾讯云,可以免费发短信 100 条,需要申请短…

直播间讨论区需要WebSocket,简单了解下

由于 http 存在一个明显的弊端(消息只能有客户端推送到服务器端,而服务器端不能主动推送到客户端),导致如果服务器如果有连续的变化,这时只能使用轮询,而轮询效率过低,并不适合。于是 WebSocket…

联想服务器-HTTP boot安装Linux系统

HTTP boot与传统PXE的主要差异 HTTP不再需要使用UDP协议的tftp服务(连接不可靠、不支持大文件)了,只需要dhcp 和http 两个服务即可,支持较稳定的大文件传输。 实验环境 ThinkSystem服务器SR650V2 SR660V2 通过HTTP boot安装Cen…

【C++基础知识学习笔记】精华版(复习专用)

常用语法 函数重载(Overload) 规则: 函数名相同 参数个数不同、参数类型不同、参数顺序不同 注意: 返回值类型与函数重载无关 调用函数时,实参的隐式类型转换可能会产生二义性 默认参数 C++ 允许函数设置默认参数,在调用时可以根据情况省略实参。规则如下: 默认参数只能…

Aop自定义注解生成日志

Aop自定义注解生成日志 1.编写自定义注解 //表示此注解可以标注在方法上 Target(ElementType.METHOD) //运行时生效 Retention(RetentionPolicy.RUNTIME) public interface OpetionLog {//定义一个变量,可以接收参数String value() default "";}2.Cont…

【移远QuecPython】EC800M物联网开发板的内置GNSS定位的恶性BUG(目前没有完全的解决方案)

【移远QuecPython】EC800M物联网开发板的内置GNSS定位的恶性BUG(目前没有完全的解决方案) GNSS配置如下: 【移远QuecPython】EC800M物联网开发板的内置GNSS定位获取(北斗、GPS和GNSS) 测试视频(包括BUG复…

智慧建筑工地管理平台源码

智慧工地是聚焦工程施工现场,紧紧围绕人、机、料、法、环等关键要素,综合运用物联网、云计算、大数据、移动计算和智能设备等软硬件信息技术,与施工生产过程相融合。 智慧工地管理平台充分运用数字化技术,聚焦施工现场岗位一线&am…

乌班图 Linux 系统 Ubuntu 23.10.1 发布更新镜像

Ubuntu 团队在其官网上发布了Ubuntu 23.10.1 版本,这是目前较新的 Ubuntu 23.10(Focal Fossa)操作系统系列的第一个发行版,旨在为社区提供最新的安装媒体。Ubuntu 22.04 LTS(Focal Fossa)操作系统系列于 2022 年 4 月 21 日发布。 Ubuntu 23.10 LTS(长期支持版本)可用…

数字人IP为何成家电品牌年轻化营销黑马?

伴随着数字人概念的出现,家电品牌逐渐通过3D虚拟数字人定制,让数字人成为内容、变现一体的IP,形成一定影响力的品牌效应,利用长线内容沉淀粉丝,使品牌实现年轻化营销。 *图片源于网络 如近日在海尔智家旗下品牌发布会上…

springboot--外部环境配置

外部环境配置 前言1、配置优先级配置文件优先级如下(后面的覆盖前面的)测试 2、外部配置3、导入配置4、属性占位符 前言 场景:线上应用如何快速修改配置,并引用最新配置? springBoot 使用配置优先级外部配置 简化配置…

重定向-缓冲区

1.重定向 文件描述符对应的分配规则是什么? 尝试用这个代码 关闭0,1,2文件描述符,看看有什么现象?关闭哪个,你打开的文件fd应该就是哪个 结论: 从0下标开始,寻找最小的没有没使用的数组位置,它…

基于php+thinkphp+vue的学生公寓管理系统-宿舍管理-寝室管理系统

运行环境 开发语言:PHP 数据库:MYSQL数据库 应用服务:apache服务器 使用框架:ThinkPHPvue 开发工具:VScode/Dreamweaver/PhpStorm等均可 项目简介 本系统结合计算机系统的结构、概念、模型、原理、方法,在计算机各种优势的情况下,采用PHP语…

自动曝光算法(第二讲)

序言 第一章说了,自动曝光算法的目的:已知当前raw图亮度、当前曝光时间、当前增益和目标亮度,当环境光发生变化的时候,是通过控制增益、曝光时间和光圈使raw图的亮度,保持在目标亮度附近。本章想讲一下目标亮度的相关…

美观且可以很方便自定义的MATLAB绘图颜色

函数介绍 主函数是draw_test,用于测试函数。 draw_h是函数,用于给Matlab提供美观且可以很方便自定义的绘图颜色。 draw_h函数介绍 这是一个带输入输出的函数,输入1/2/3,输出下面三种颜色库的配色,每种库均有五种颜色…

Find My手机保护壳|苹果Find My与手机保护壳结合,智能防丢,全球定位

随着科技水平的快速发展,科技美容这一行业做为新型产业新生而出。时尚IT品牌随着市场的多元化发展。针对手机品牌和功能的增加而呈多样化,将手机保护壳按质地分有PC壳,皮革 ,硅胶,布料,硬塑,皮套…

手机上有哪些支持设置农历日期提醒的工具

很多人的生日都是按照农历日期来安排的,而农历日期和公历日期相错的日子很多,在手机上如果想要记录农历生日提醒,需要借助一些支持设定农历日期的工具来实现。 手机上有哪些支持设置农历日期提醒的工具呢?敬业签是一款可以在手机…

C++:类和对象(中)

1.类的6个默认成员函数: 如果一个类中什么成员都没有,简称为空类。 空类中真的什么都没有吗?并不是,任何类在什么都不写时,编译器会自动生成以下6个默认成员函数。 默认成员函数:用户没有显式实现&#xff…

全志H6-LicheePi调试记录

LPDDR3:K4E6E304ED-EGCG和K4E6E304ED-EGCF 最高速率不一样,CG是2133MHz,CF是1866MHz