大地测量反演模型准外部检验的进一步探讨

反演模型质量的好坏对反演结果及后期的反演解释影响很大.本文给出了大地测量反演模型准外部检验的基本理论,提出了观测数据的划分方法和反演模型准外部质量控制的标准,当没有先验信息对模型进行检验和约束且数据量很大、很复杂时,该法是一种很有效的检验方法.     

顾及协方差函数的自适应四叉树InSAR数据压缩算法

利用InSAR数据进行形变机理反演时,由于InSAR数据点众多,且含有较多的误差乃至粗差点,严重制约了反演的效率和可靠性。为此,本文提出了顾及InSAR数据物理空间相关性来设立协方差函数,并依据协方差函数确定四叉树象限分解阈值和最大象限大小的自适应四叉树分解InSAR数据压缩算法。本算法能够在尽可能保留形变信号特征细节信息的同时,极大地降低InSAR数据量。本文以西安地区地面沉降InSAR结果为例进行了实验分析,验证了该算法的有效性,结果表明该方法能够在不损失形变信号特征的同时,有效的实现InSAR数据压缩和噪声消除的目的。     

基于S-RVoG模型的PolInSAR森林高度非线性复数最小二乘反演算法

针对经典的PolInSAR森林高度三阶段几何反演算法在单基线条件容易受到地体幅度比假设以及地形坡度影响的问题,从测量平差角度提出了基于S-RVoG模型的PolInSAR非线性复数最小二乘森林高度反演算法。该算法不再需要假设某一个极化通道地体幅度比为零,且采用考虑地形坡度影响的S-RVoG模型作为平差模型。为了验证算法,本文采用欧空局BioSAR2008项目提供的3景P波段极化干涉SAR数据进行两组单基线森林高度反演试验。结果表明,在单基线条件下,基于RVoG模型的非线性复数最小二乘算法反演结果优于三阶段几何反演算法,而基于S-RVoG模型的非线性复数最小二乘算法进一步提高反演精度,对于坡度较大区域（坡度>10°）,精度平均提高了18.48%。     

粗差探测的Bayes方法

讨论粗差探测的Bayes方法。首先根据Bayes统计推断的基本原理,建立判断粗差的Bayes方法———后验概率法,然后针对测量平差实际,分别给出非等权独立观测条件下基于均值漂移模型和方差膨胀模型的后验概率的计算公式,最后结合一边角网算例,验证本文方法的效果。     

附不等式约束平差模型的一种快速搜索算法

大地测量中常存在一些先验不等式约束信息,充分利用它们可以保证参数解的唯一性和稳定性。然而,现有的不等式约束平差算法主要是基于优化理论,算法通常比较复杂,需要选取有效约束或建立罚函数。在最小二乘平差准则基础上,把不等式约束看成是一个可行域,借助Fisher函数在可行域中快速搜索使误差平方和达到最小的最优解,推导出了可行解为最优解的充分必要条件。建立了基于Wolfe-Powell算法的非精确快速搜索算法,从而减小了搜索算法的计算量,得到了一种新的不等式约束平差计算方法。该算法的平差准则与最小二乘平差准则一致,不需要矩阵求逆运算,可适用于维数较大的平差问题解算。     

SBAS-InSAR方法支持下的矿区地表沉降监测及参数反演

煤矿开采极易诱发各类地质灾害,从而对矿区人民的生命财产安全带来严重影响。因此,对矿区地表进行形变监测及对其开采面参数进行反演获得形变规律,可有效防止矿区地质灾害的发生。本文选取山西省的店坪煤矿和汾源煤业作为研究区域,基于Sentinel-1A卫星影像数据,使用小基线集（SBAS）方法提取研究区的形变结果。使用Okada模型对5-101工作面参数进行反演和分析,并将反演结果与采工图实测数据进行对比验证。结果表明：①煤矿开采引起的地表形变程度与煤矿开采长度成正比;②急倾斜厚煤层开采引发的地表形变程度要比近水平薄煤层引起的地表形变程度大,更容易造成地表塌陷;③Okada模型反演较短近水平薄煤层的开采参数精度较高,说明该模型适用于此类型煤矿开采参数的反演。     

2016年意大利阿马特里切Mw 6.2地震震源机制InSAR反演

2016年8月24日，意大利中部阿马特里切（Amatrice）地区发生Mw 6.2地震。采用ALOS-2条带模式和SENTINEL-1A宽幅模式的合成孔径雷达（synthetic aperture radar，SAR）数据分别进行SAR差分干涉测量处理，获取了该地震的同震形变场。结果显示，本次地震造成意大利中部地区发生明显的地壳形变，在雷达视线向最大沉降量达19.6 cm。基于合成孔径雷达干涉测量（interferometry synthetic aperture radar，InSAR）和GPS同震形变场数据对此次地震的发震断层进行联合反演，通过改进倾角和平滑系数获取方法，得到了最优滑动分布模型。通过使用单断层模型和双断层模型进行反演可知，双断层模型反演结果优于单断层反演结果，两种模型下反演模型相关系数分别为0.85和0.89，发震断层走向分别为160°和158°，倾角分别为44°和46°，倾滑分布主要位于地下5~7 km，平均倾滑角为-80°，最大倾滑量0.9 m位于地壳深度5 km处，该发震断层是亚平宁冲断带的一部分，为NW-SE向延伸的正断层，断层长约20 km。综合使用地震同震形变场和GPS数据对震源机制进行反演、模拟和分析，获取了高精度的震源参数，可以为分析地震危险性和断层破裂参数等提供数据支持。     

大地测量和地震数据联合反演地震震源破裂过程研究进展

综述了近10多年来大地测量和地震数据联合反演地震震源破裂过程的研究成果及最新进展,特别是大地测量和地震数据联合反演地震震源破裂过程的模型、算法及典型震例。展望了利用GPSI、nSAR和地震波数据联合反演震源破裂过程中顾及先验信息、附有不等式约束时反演模型的建立方法,以及相对权比的确定、全局优化方法的应用等问题。     

综合地表与深部位移监测数据的滑坡多目标加权位移反分析方法

针对滑坡参数反演存在的多目标优化问题,同时为弥补滑坡深部位移监测点位稀疏的不足,提出了一种综合地表与深部位移监测数据的滑坡多目标加权位移反分析方法。该方法利用地表与深部位移监测数据构建多目标位移反分析模型,采用抗差Helmert方差分量估计法计算各类位移观测量的抗差验后随机模型,优化反演模型的权参数。试验结果表明:深部位移信息量不足会导致位移反分析结果出现严重偏差,综合地表与深部位移信息的反演结果更为准确;基于抗差Helmert方差分量估计的多目标加权位移反分析方法不仅可以合理确定不同类型观测数据的权重,还能够有效抵制异常粗差对反演结果的影响,提高了反演计算精度。     

Bootstrap方法在遥感线性BRDF模型反演中的应用

在遥感反演中,通常假设反演参数和模型误差的先验分布服从正态分布,这个假设通常不太符合实际。为此,本文提出由Bootstrap方法估计反演参数和模型误差的先验分布的方案。同时对先验数据按照地物分类,统计假设检验表明将先验知识分类的合理性。最后,以RossThick-LiTransit核组合的线性核驱动BRDF模型为例,用NOAA-AVHRR观测数据对使用Bootstrap方法的反演算法进行试验,并与正态假设下的Tikhonov正则化反演和Bayes反演结果比较,说明对先验知识分类和使用Bootstrap方法的遥感反演方法能明显减小参数反演结果的不确定性,提高其可信度。