大地测量学与导航

短基线集形变模型反演的正则化解算方法

  • 姜兆英 ,
  • 刘国林 ,
  • 陶秋香
展开
  • 1. 山东科技大学测绘科学与工程学院, 山东 青岛 266590;
    2. 青岛农业大学理学与信息科学学院, 山东 青岛 266109
姜兆英(1981-),女,博士生,研究方向为InSAR数据处理和模型算法。E-mail: 15964287012@163.com

收稿日期: 2015-03-17

  修回日期: 2016-01-04

  网络出版日期: 2016-05-30

基金资助

国家自然科学基金(41274007;41404003);山东省自然科学基金(ZR2012DM001);高等学校博士学科点专项科研基金(20123718110001);山东省泰山学者建设工程专项经费(TSXZ201509)

Regularization solution of Small Baseline Subset Deformation Model Inversion

  • JIANG Zhaoying ,
  • LIU Guolin ,
  • TAO Qiuxiang
Expand
  • 1. Geomatics College, Shandong University of Science and Technology, Qingdao 266590, China;
    2. College of Science and Information, Qingdao Agricultural University, Qingdao 266109, ChinaAbstract

Received date: 2015-03-17

  Revised date: 2016-01-04

  Online published: 2016-05-30

Supported by

The National Natural Science Foundation of China (Nos. 41274007;41404003);Shandong Province Natural Science Foundation of China (No. ZR2012DM001);Specialized Research Fund for the Doctoral Program of Higher Education (No. 20123718110001);Shandong Taishan Scholar Construction Project under Special Funding (No. TSXZ201509)

摘要

针对短基线集形变模型反演中法方程系数矩阵呈病态的问题,提出一种正则化稳健解算方法。该方法基于Tikhonov正则化理论,将形变速率求解问题转化为极小化问题,根据L-曲线法选取正则化参数,考虑最小二乘残差各个分量间的关系选取正则化矩阵,实现短基线集形变模型反演的稳健解算。分别采用LS法、岭估计法和Tikhonov正则化法对覆盖北京地区的29景ENVISAT ASAR数据进行处理,反演出研究区沉降速率图。通过对代表不同沉降情况的21个点的均方误差值和时间相干值、整个研究区的均方误差图等的对比分析,表明本文提出的短基线集形变模型反演的正则化稳健解算方法可获取更可靠的形变监测结果。

本文引用格式

姜兆英 , 刘国林 , 陶秋香 . 短基线集形变模型反演的正则化解算方法[J]. 测绘学报, 2016 , 45(5) : 566 -573 . DOI: 10.11947/j.AGCS.2016.20150143

Abstract

For the coefficient matrix of the normal equation is ill-conditioned during inverting deformation model of small baseline subset (SBAS) InSAR technique, a regularization robust method is proposed. Based on Tikhonov regularization theory, this method converts the problem of how to solve the deformation rate into minimization problem. According to L-curve method to choose regularization parameter, considering the relationship between the individual components of least-squares residuals to choose regularization matrix, thus it achieves robust solution of SBAS deformation model inversion. We adopt respectively least-squares estimation, ridge estimation and Tikhonov regularization method to deal with 29 ENVISAT ASAR dataset relevant to the Beijing area, achieving the subsidence rate map of the study area. Through comparative analysis among the mean square error (MSE) of 21 points on behalf of the different subsidence, temporal coherence values and MSE maps of the entire study area, we confirm that Tikhonov regularization robust method in inverting SBAS deformation model can obtain more reliable results of deformation monitoring.

参考文献

[1] BERARDINO P, FORNARO G, LANARI R, et al. A New Algorithm for Surface Deformation Monitoring Based on Small Baseline Differential SAR Interferograms[J]. IEEE Transactions on Geoscience and Remote Sensing, 2002, 40(11): 2375-2383.
[2] TIZZANI P, BERARDINO P, CASU F, et al. Surface Deformation of Long Valley Caldera and Mono Basin, California, Investigated with the SBAS-InSAR Approach[J]. Remote Sensing of Environment, 2007, 108(3): 277-289.
[3] 尹宏杰, 朱建军, 李志伟, 等. 基于SBAS的矿区形变监测研究[J]. 测绘学报, 2011, 40(1): 52-58. YIN Hongjie, ZHU Jianjun, LI Zhiwei, et al. Ground Subsidence Monitoring in Mining Area Using DInSAR SBAS Algorithm[J]. Acta Geodaetica et Cartographica Sinica, 2011, 40(1): 52-58.
[4] 李珊珊, 李志伟, 胡俊, 等. SBAS-InSAR技术监测青藏高原季节性冻土形变[J]. 地球物理学报, 2013, 56(5): 1476-1486. LI Shanshan, LI Zhiwei, HU Jun, et al. Investigation of the Seasonal Oscillation of the Permafrost over Qinghai-Tibet Plateau with SBAS-InSAR Algorithm[J]. Chinese Journal of Geophysics, 2013, 56(5): 1476-1486.
[5] 吴宏安, 张永红, 陈晓勇, 等. 基于小基线DInSAR技术监测太原市2003-2009年地表形变场[J]. 地球物理学报, 2011, 54(3): 673-680. WU Hongan, ZHANG Yonghong, CHEN Xiaoyong, et al. Ground Deformation Monitoring Using Small Baseline DInSAR Technique: A Case Study in Taiyuan City from 2003 to 2009[J]. Chinese Journal of Geophysics, 2011, 54(3): 673-680.
[6] 许文斌, 李志伟, 丁晓利, 等. 利用InSAR短基线技术估计洛杉矶地区的地表时序形变和含水层参数[J]. 地球物理学报, 2012, 55(2): 452-461. XU Wenbin, LI Zhiwei, DING Xiaoli, et al. Application of Small Baseline Subsets D-InSAR Technology to Estimate the Time Series Land Deformation and Aquifer Storage Coefficients of Los Angeles Area[J]. Chinese Journal of Geophysics, 2012, 55(2): 452-461.
[7] 张永红,张继贤,龚文瑜,等.基于SAR干涉点目标分析技术的城市地表形变监测[J]. 测绘学报, 2009, 38(6): 482-487, 493. DOI: 10.3321/j.issn:1001-1595.2009.06.003. ZHANG Yonghong, ZHANG Jixian, GONG Wenyu, et al. Monitoring Urban Subsidence Based on SAR Interferometric Point Target Analysis[J]. Acta Geodaetica et Cartographica Sinica, 2009, 38(6): 482-487, 493. DOI: 10.3321/j.issn:1001-1595.2009.06.003.
[8] 张永红, 吴宏安, 孙广通. 时间序列InSAR技术中的形变模型研究[J]. 测绘学报, 2012, 41(6): 864-869, 876. ZHANG Yonghong, WU Hongan, SUN Guangtong. Deformation Model of Time Series Interferometric SAR Techniques[J]. Acta Geodaetica et Cartographica Sinica, 2012, 41(6): 864-869, 876.
[9] FERRETTI A, PRATI C, ROCCA F. Nonlinear Subsidence Rate Estimation Using Permanent Scatterers in Differential SAR Interferometry[J]. IEEE Transactions on Geoscience and Remote Sensing, 2000, 38(5): 2202-2212.
[10] 李永生, 张景发, 李振洪, 等. 利用短基线集干涉测量时序分析方法监测北京市地面沉降[J]. 武汉大学学报(信息科学版), 2013, 38(11): 1374-1377. LI Yongsheng, ZHANG Jingfa, LI Zhenhong, et al. Land Subsidence in Beijing City from InSAR Time Series Analysis with Small Baseline Subset[J]. Geomatics and Information Science of Wuhan University, 2013, 38(11): 1374-1377.
[11] LANARI R, MORA O, MANUNTA M, et al. A Small-baseline Approach for Investigating Deformations on Full-resolution Differential SAR Interferograms[J]. IEEE Transactions on Geoscience and Remote Sensing, 2004, 42(7): 1377-1386.
[12] BONANO M, MANUNTA M, MARSELLA M, et al. Long-term ERS/ENVISAT Deformation Time-series Generation at Full Spatial Resolution via the Extended SBAS Technique[J]. International Journal of Remote Sensing, 2012, 33(15): 4756-4783.
[13] AKBARI V, MOTAGH M. Improved Ground Subsidence Monitoring Using Small Baseline SAR Interferograms and a Weighted Least Squares Inversion Algorithm[J]. IEEE Geoscience and Remote Sensing Letters, 2012, 9(3): 437-441.
[14] LAUKNES T R, ZEBKER H A, LARSEN Y. InSAR Deformation Time Series Using an L1-norm Small-baseline Approach[J]. IEEE Transactions on Geoscience and Remote Sensing, 2011, 49(1): 536-546.
[15] 王家映. 地球物理反演理论[M]. 北京: 高等教育出版社, 2002. WANG Jiaying. Inverse Theory in Geophysics[M]. Beijing: Higher Education Press, 2002.
[16] TIKHONOV A N, ARSENIN V Y. Solutions of Ill-posed Problems[M]. New York: John Wiley and Sons, 1977.
[17] HE Guoqiang, LIU Linxian. A Kind of Implicit Iterative Methods for Ill-posed Operator Equation[J]. Journal of Computational Mathematics, 1999, 17(3): 275-284.
[18] STEIHAUG T.The Conjugate Gradient Method and Trust Regions in Large Scale Optimization[J]. SIAM Journal on Numerical Analysis, 1983, 20(3): 626-637.
[19] 郭建峰. 测量平差系统病态性的诊断与处理[D]. 郑州: 信息工程大学, 2002. GUO Jianfeng. Diagnostics and Processing of Ill-posed in Surveying Adjustment System[D]. Zhengzhou: Information Engineering University, 2002.
[20] 王松桂, 贾忠贞. 矩阵论中不等式[M]. 合肥: 安徽教育出版社, 1994. WANG Songgui, JIA Zhongzhen. Inequality in Matrix Theory[M]. Hefei: Anhui Education Press, 1994.
[21] 王彦飞. 反演问题的计算方法及其应用[M]. 北京: 高等教育出版社, 2007. WANG Yanfei. Computational Methods for Inverse Problems and Their Applications[M]. Beijing: Higher Education Press, 2007.
[22] 王振杰. 测量中不适定问题的正则化解法[M].北京:科学出版社, 2006. WANG Zhenjie. Regularization Solutions of Ill-posed Problems in Geodesy[M]. Beijing: Science Press, 2006.
[23] 杨文采. 地球物理反演和地震层析成像[M]. 北京: 地质出版社, 1989. YANG Wencai. Geophysical Inversion and Seismic Tomography[M]. Beijing: Geological Press, 1989.
[24] 张方仁, 汪晓庆. 平差参数的岭估计和压缩估计[J]. 武汉测绘科技大学学报, 1989, 14(3): 46-58. ZHANG Fangren,WANG Xiaoqing.The Ridge Estimation and the Stein Estimation of Adjustment Parameters[J]. Journal of Wuhan Technical University of Surveying and Mapping, 1989, 14(3): 46-58.
[25] 陶秋香, 刘国林. PS InSAR公共主影像优化选取的一种新方法[J]. 武汉大学学报(信息科学版), 2011, 36(12): 1456-1460. TAO Qiuxiang, LIU Guolin. A New Method to Optimize and Select Common Master Images in PS InSAR[J]. Geomatics and Information Science of Wuhan University, 2011, 36(12): 1456-1460.
[26] PEPE A, LANARI R. On the Extension of the Minimum Cost Flow Algorithm for Phase Unwrapping of Multitemporal Differential SAR Interferograms[J]. IEEE Transactions on Geoscience and Remote Sensing, 2006, 44(9): 2374-2383.
文章导航

/