基于现代测量平差的InSAR三维形变估计理论与方法

【摘要】合成孔径雷达干涉测量(InSAR)是二十世纪六十年代发展起来一项空间大地测量与现代遥感技术,可以大范围、高精度和空间连续地监测地球表面发生的形变,为人类了解地球运动机制、防御自然及人为地质灾害等提供了重要的手段。然而,传统的单轨InSAR形变测量结果只能反映真实的三维地表形变在雷达视线方向(LOS)上的一维投影,容易引起对地表形变监测信息的误判或漏判,使得该技术的发展和应用受到了极大的阻碍。通过融合多颗SAR卫星平台获取的多个轨道数据,可以将InSAR一维形变测量拓展至三维。但是现有方法普遍缺乏严密的数据处理理论,在处理多平台、多轨道和多时域InSAR测量值时,无法充分考虑这些多源异质数据在时间尺度和测量精度上的差异,而且不能很好地顾及InSAR多源误差的影响,从而极大的限制了三维形变测量的精度和时间分辨率。现代测量平差是一种以误差理论和经典测量平差为基础的数据处理技术,并在不同层面上扩充、发展形成新理论和新方法,目前已被广泛应用于“3S”及其集成的数据处理中。InSAR三维形变测量本质上可以视为一个平差问题,因此现代测量平差理论和方法在该领域中具有良好的应用潜力。本文以InSAR三维形变估计的函数模型、随机模型和估计准则为核心研究内容。针对InSAR三维形变估计理论中的关键技术问题,借鉴和研究了滤波、插值、最小二乘、方差分量估计、卡尔曼滤波和总体最小二乘等现代测量平差方法的原理和思想,系统性地研究了基于现代测量平差的InSAR三维形变估计理论与方法,为高精度三维地表形变的获取提供了新的途径。本文的主要贡献和创新之处在于：(1)建立了基于最小二乘的InSAR三维地表形变估计的函数模型,并针对模型的观测值设计了InSAR趋势误差改正算法,显著改善了InSAR三维形变估计结果的精度。本文深入分析了从InSAR的LOS方向一维测量值求解三维地表形变的可能性和精度,系统地研究了多源异质InSAR(包括D-InSAR、Offset-Tracking和MAI)测量值与三维地表形变之间的函数关系,并通过理论推导将其扩展为适用于融合多源异质InSAR测量值和GPS测量值估计三维地表形变的函数模型。同时,为了进一步提高MAI技术获取的方位向形变测量值的精度,本文通过研究电离层误差和形变信号在空间上的不同分布规律,发展了基于方向性滤波和插值的MAI电离层趋势误差改正方法；而针对D-InSAR技术在毗邻轨道获取的形变结果的不一致问题,本文利用在毗邻轨道的公共区域提取的连接点作为附加约束,发展了基于整体同步最小二乘的多轨道D-InSAR轨道趋势误差改正模型。在此基础上,本文分别精确获取了2010年新西兰Darfield地震和2011年日本Tohoku-Oki地震的第一个同震三维地表形变场,精度可达到cm级。(2)提出了基于方差分量估计的InS(?)R三维地表形变速率估计算法,实现了InSAR测量值随机模型的后验估计,显著改善了InSAR三维形变速率估计结果的精度。在融合多源异质InSAR及GPS测量值估计三维地表形变速率时,除了需要建立精确的观测函数模型之外,还需要获取InSAR及GPS测量值的随机模型(即方差),然而目前没有一个公认的可以精确确定多源异质InSAR测量值先验方差的方法。针对这一难题,本文首先根据多源异质InSAR及GPS测量噪声的统计特性实现了多源异质InSAR和GPS测量值的最优分组；然后基于线性模型思想建立了最小二乘残差和多源异质InSAR及GPS测量值的单位权中误差之间的函数模型,通过迭代运算得到了多源异质InSAR和GPS测量值的后验方差(或权重)。此外,考虑到在实际中InSAR多余观测量较少,本文在顾及三维形变速率估计精度和计算量的前提下,通过大量实验研究得到了方差分量估计中InSAR测量值的最优配置。模拟和美国南加州地区的真实数据实验结果表明,本文方法的三维形变速率估计精度比现有方法有了明显改善。随后,将该方法应用于青藏高原冬克玛底冰川的监测中,揭示了青藏高原山岳冰川的第一个三维运动速率场,为研究青藏高原的冰水质量迁移和转换提供了一个新的视角。(3)提出了基于卡尔曼滤波的InSAR三维地表形变序列动态估计函数模型,显著提高了InSAR三维地表形变测量结果的时间分辨率和精度。传统的InSAR三维形变测量方法受限于不同平台、不同轨道SAR数据的时态差异,只能对多源异质InSAR测量值进行静态平差求解,获取地表的瞬时三维形变或三维形变速率,而无法估计三维形变序列。本文首先深入分析了融合多平台、多轨道和多时域InSAR测量值求解三维地表形变序列的局限性；然后通过联合时序InSAR测量值的成像几何建立空间联系和引入三维形变速率矢量作为时态约束,建立了三维地表形变的观测方程和状态方程,从而实现了三维形变序列的动态估计。为了提高三维形变序列的估计精度,本文引入方差扩大模型来抑制InSAR测量值中解缠误差等粗差的影响。而在实际数据处理中,利用基于最小二乘的InSAR时序测量值处理方法得到的三维低频形变和地形残差来确定滤波初始状态矢量和协助InSAR干涉相位解缠,从而进一步提高了卡尔曼滤波方程组的解算精度。模拟和美国南加州地区的真实数据实验结果表明,本文方法适用于线性或非线性形变监测,实现了三维地表形变序列的准实时估计,提高了计算效率的同时降低了数据的存储量,可以得到时间分辨率和精度明显优于现有方法的三维形变序列结果。(4)提出了基于总体最小二乘的InSAR三维地表形变估计方法,解决了传统的基于高斯-马尔可夫(GM)模型的估计准则无法考虑DEM误差影响的问题,提高了InSAR三维形变估计结果的精度。在利用InSAR技术监测地表平行位移的三维形变时,通常都借助于DEM数据提供的坡度角信息来降低对InSAR观测量的要求,这使得InSAR三维形变监测模型的系数矩阵被DEM误差所污染。传统的基于GM模型的估计方法(如最小二乘)只能考虑观测值含误差的情况,本文方法则以变量含误差(EIV)模型为核心,在满足观测值误差平方和最小的同时,使得系数矩阵误差平方和也达到最小,因此可以得到更加优化的三维地表形变估值。针对InSAR测量最为常见的配置情况,本文首先设计了两种地表平行位移的InSAR三维形变监测模型,一种是基于升降轨InSAR的LOS向测量值,另一种是基于单轨InS(?)R的LOS向和方位向测量值；然后分别研究了DEM误差(或坡度角方差)在这两种模型中的传播方式,并通过理论推导得到了相应的系数矩阵方差,进而实现了三维地表形变的总体最小二乘估计。模拟和河南大寺坡真实数据实验结果表明,总体最小二乘在估计三维地表形变时比现有方法表现更加稳健,提高了三维形变结果的估计精度。
【作者】胡俊；
【导师】朱建军；李志伟；丁晓利；
【作者基本信息】中南大学，测绘科学与技术，2013，博士
【关键词】合成孔径雷达干涉测量(InSAR)；现代测量平差；三维形变估计；函数模型；随机模型；

【参考文献】
[1]高鸣.人羊膜修复腱鞘缺损的实验研究[D].河北医科大学,外科学,2013,硕士.
[2]霍崇卓.抵押权实现程序研究[D].安徽大学,法律,2014,硕士.
[3]汤永清.基于共性分析的MIS设计与实现[D].华侨大学,电工理论与新技术,2003,硕士.
[4]杜庆波.数据专线的测试[J].铁道通信信号.2005(01)
[5]石琴.基于现代设计理论的车身结构设计方法研究[D].合肥工业大学,2006.
[6]张吉人.SBS改性沥青防水卷材屋面工程施工及质量控制[J].建筑技术.1997(06)
[7]何永祥.论法律裁决的不谬性[D].西南政法大学,法学理论,2012,硕士.
[8]张颖.大连邮电设计院实施项目经理制改造研究[D].大连理工大学,工商管理,2004,硕士.
[9]陈德顺.Mg-Li-Al-Si合金的制备及其组织性能研究[D].重庆大学,材料工程（专业学位）,2014,硕士.
[10]曹庆连.国际中转型枢纽陆港发展思路与策略研究[D].长安大学,物流工程与管理,2013,硕士.
[11]张蓉.民用飞机系统故障诊断测试性分析与评估技术研究[D].南京航空航天大学,载运工具运用工程,2012,硕士.
[12]邓安燕.农村初中历史有效教学策略研究[D].云南师范大学,历史课程与教学论,2013,硕士.
[13]张利军.责任编辑:图书质量的守望者[J].出版广角,2014,11:33-35.
[14]付静.水环境重金属检测的电化学传感器的研究[D].浙江大学,2007.
[15]王恒宝.冠状动脉旁路移植术治疗冠状动脉三支血管病变的临床分析[D].吉林大学,外科学,2014,硕士.
[16]田忠.产业集聚视角下的就业效应研究[D].贵州财经大学,劳动经济学,2013,硕士.
[17]黄河.我国网球方向硕士研究生培养状况的研究[D].北京体育大学,体育教育训练学,2013,硕士.
[18]毛博.基于RFID的深井油田注水控制系统设计与实现[D].华中科技大学,机械电子工程,2013,硕士.
[19]周炳华.炼厂多网络集成优化方法研究[D].东北大学,系统工程,2010,硕士.
[20]钟杰.静电火花试验系统及其静电点火能量计算方法研究[D].重庆大学,电气工程,2014,硕士.
[21]马杰.印度梨形孢诱导烟草促生、抗病、抗逆作用及其机理的初步研究[D].浙江大学,植物保护学（专业学位）,2012,硕士.
[22]任怀强.基于驱动层的计算机数据安全防护系统的设计与实现[D].西安电子科技大学,计算机软件与理论,2011,硕士.
[23]游玉霞.苏联的中东政策研究（1964-1982年）[D].黑龙江大学,世界史,2013,硕士.
[24]赖俊辉.基于iPad平台的在线办公系统设计与实现[D].北京交通大学,软件工程（专业学位）,2013,硕士.
[25]段鹏.基于粗糙集和神经网络的不完备信息系统数据挖掘研究[D].云南师范大学,2004.
[26]范雪.论省级电台的创新与发展[D].渤海大学,新闻学,2014,硕士.
[27]何雅.论银行信托理财业务的法律规制[D].辽宁大学,经济法学,2012,硕士.
[28]石淼.股份有限公司中小股东权益保护制度研究[D].苏州大学,法律,2003,硕士.
[29]王磊.钢丝绳综合实验系统的结构设计与研究[D].太原科技大学,机械设计及理论,2013,硕士.
[30]丁然.未来通信中的语音图像数据格式设计与仿真[D].成都理工大学,电子与通信工程,2013,硕士.
[31]张程.色散管理孤子掺铒光纤锁模激光器脉冲压缩特性的研究[D].天津大学,2013.
[32]沈丽芳,冯志仙,赵雪红,徐红,沈丽娜,徐巧君,黄丽华.急性心肌梗死单病种护理质量评价标准的研究[J].护理与康复,2014,04:310-313.
[33]王玉丽.基于石墨烯的电化学传感器的制备及应用[D].安徽师范大学,2011.
[34]周威.多菌灵、百菌清、毒死蜱在大棚和露地蔬菜及其土壤中的残留特征[D].浙江大学,植物保护,2013,硕士.
[35]耿精媛.带有放弃的On-off源流体队列模型[D].长安大学,应用数学,2014,硕士.
[36]苗世洪,马帅,尚亚男,姜臻,侯俊贤,于之虹.基于割点和路径搜索的输电断面快速识别方法[J].电力系统自动化,2014,02:39-45.
[37]闻千.网络化运营条件下城市轨道交通应急指挥管理评价方法研究[D].西南交通大学,交通运输规划与管理,2013,硕士.
[38]王邦军.中国软件企业管理创新研究[D].武汉理工大学,管理科学与工程,2004,硕士.
[39]栾靖.小学生田径类基础动作技能发展水平与评价研究[D].南京师范大学,体育教育训练学,2012,硕士.
[40]罗依儒.英语教师对话智慧研究[D].宁波大学,课程与教学论,2013,硕士.
[41]鲁烈琴.三维叠前深度偏移技术在西部复杂地区的应用与效果[J].石油地球物理勘探,1997,05:703-708+750.
[42]钱旭培.基于DFL的多Agent在对抗环境和非对抗环境中的协调模型及应用研究[D].苏州大学,计算机应用技术,2003,硕士.
[43]王志斌.虚拟环境下基于物理属性的线缆布局与装配仿真技术[D].北京理工大学,机械工程,2014,硕士.
[44]袁光.CucurbitacinⅠ诱导胶质母细胞瘤凋亡和保护性自噬的机制研究[D].山东大学,外科学,2014,博士.
[45]瞿绍生.农村信用社内部责任机制研究[D].湘潭大学,法律,2011,硕士.
[46]苏学灵,纪昌明,黄小锋,贾东旭,郭希海.混合式抽水蓄能电站在梯级水电站群中的优化调度[J].电力系统自动化,2010,04:29-33.
[47]陈璇.《红楼梦》续书研究[D].苏州大学,中国古代文学,2003,硕士.
[48]金玲玲.带外场的Vlasov-Poisson-BGK方程的全局解[D].华中科技大学,应用数学,2013,硕士.
[49]李淑华.区域房地产市场分析框架及技术研究[D].华东师范大学,企业管理,2004,硕士.
[50]张帆.无线公钥基础设施—WPKI设计与实现研究[D].西安电子科技大学,计算机系统结构,2004,硕士.