摘要
针对不连续地形难以精确重建地形高程的问题,将干涉图残差点和相邻像素点差异比较相结合的坏点判断方法与基于全变分能量函数的最大后验估计算法结合起来,该文提出一种改进的多基线InSAR高程重建方法。该方法利用干涉图残差点和相邻像素点差异比较对不连续点计算原理的不同,对基于全变分能量函数的最大后验估计算法重建的预估地形高程图进行坏点判断,有效判断出高程图不连续点中由于噪声或计算误差引起的坏点,并对预估地形高程图中坏点进行适当加权更新获得最终估计的地形高程图。该文算法在保留基于全变分能量函数的最大后验估计算法简便的优点的同时提高了目标地形高程的估计精度。最后通过2种不同类型地形的仿真实验,表明该方法简单有效。
This paper proposes an improved multi-baseline In SAR height reconstruction approach to reconstruct height models of discontinuous terrains through combining with the Maximum A Posteriori(MAP) estimation based on Total Variation(TV) model energy function and a defective pixel judgment method composed of interferogram residues and comparison of adjacent pixels differences. The proposed method utilizes the difference discontinuous pixels between interferogram residues and comparison of adjacent pixels differences, to make defective pixels judgment from reconstructed terrain graphs obtained by the MAP estimation based on TV model energy function, which estimates efficiently defective pixels caused by errors or noise from discontinuous pixels, and then updates relatively defective pixels to obtain the final terrain height estimation. The presented algorithm not only keeps the convenience of MAP estimation algorithm but also improves estimation precision for height models of objective terrains. Experiences of two different types of terrain show that the proposed approach is efficient as well as valid.
引文
[1]XU G,XING M D,XIA X G,et al.Sparse regularization of interferometric phase and amplitude for In SAR image formation based on Bayesian representation[J].IEEEGeoscience and Remote Sensing Society,2015,53(4):2123-2136.doi:10.1109/TGRS.2014.2355592.
[2]VEHMAS R,JYLHA J,VAILA M,et al.Analysis and comparison of multichannel SAR imaging algorithms[C].IEEE Radar Conference,Seattle,USA,2017:340-345.doi:10.1109/RADAR.2017.7944224.
[3]YU H W,LAN Y,XU J Y,et al.Large-scale L0-norm and L1-norm 2D phase unwrapping[J].IEEE Geoscience and Remote Sensing Society,2017,55(8):4712-4728.doi:10.1109/TGRS.2017.2698452.
[4]袁志辉,邓云凯,李飞,等.改进的基于最大似然估计的多通道In SAR高程重建方法[J].电子与信息学报,2013,35(9):2161-2167.doi:10.3724/SP.J.1146.2012.01530.YUAN Zhihui,DENG Yunkai,LI Fei,et al.Improved multichannel In SAR height reconstruction method based on maximum likelihood estimation[J].Journal of Electronics&Information Technology,2013,35(9):2161-2167.doi:10.3724/SP.J.1146.2012.01530.
[5]XU W,CHANG E,KWOH L,et al.Phase-unwrapping of SAR interferogram with multi-frequency or multi-baseline[C].Proceedings of IEEE International Geoscience and Remote Sensing Symposium,Pasadena Califormia,USA,1994:730-732.
[6]THOMPSON D G,ROBERTSON A E,ARNOLD D V,et al.Multi-baseline interferometric SAR for iterative height estimation[C].Proceedings of IEEE International Geoscience and Remote Sensing Symposium,Hamburg,Germany,1999:251-253.
[7]PASCAZIO V and SCHIRINZI G.Multifrequency In SARheight reconstruction through maximum likelihood estimation of local planes parameters[J].IEEE Transactions on Image Processing,2002,11(12):1478-1489.doi:10.1109/TIP.2002.804274.
[8]BASELICE F,FERRAIOLI G,PASCAZIO V,et al.Contextual information-based multichannel synthetic aperture radar interferometry:Addressing DEMreconstruction using contextual information[J].IEEE Signal Processing Magazine,2014,31(4):59-68.doi:10.1109/MSP.2014.2312282.
[9]FERRAIUOLO G,PASCAZIO V,and SCHIRINZI G.Maximum a posteriori estimation of height profiles in In SARimaging[J].IEEE Geoscience and Remote Sensing Letters,2004,1(2):66-70.doi:10.1109/LGRS.2003.822882.
[10]FERRAIOLI G,SHABOU A,TUPIN F,et al.Multichannel phase unwrapping with graph cuts[J].IEEE Geoscience and Remote Sensing Letters,2009,6(3):562-566.doi:10.1109/LGRS.2009.2021165.
[11]SHABOU A,BASELICE F,and FERRAIOLI G.Urban digital elevation model reconstruction using very high resolution multichannel In SAR data[J].IEEE Transactions on Geoscience and Remote Sensing,2012,50(11):4748-4758.doi:10.1109/TGRS.2012.2191155.
[12]HONG F,TANG J W,and LU P P.Multichannel DEMreconstruction method based on Markov random fields for bistatic SAR[J].Science China Information Sciences,2015,58(6):1-14.doi:10.1007/s11432-015-5320-z.
[13]LI Zhengfang,BAO Zheng,and SUO Zhiyong.A joint image coregistration,phase noise suppression,and phase unwrapping method based on subspace projection for multibaseline In SAR systems[J].IEEE Transactions on Geoscience and Remote Sensing,2007,45(3):584-591.doi:10.1109/TGRS.2006.888143.
[14]张红敏,靳国旺,徐青,等.中国余数定理在双基线In SAR相位解缠中的应用[J].测绘学报,2011(6):770-777.ZHANG Hongmin,JIN Guowang,XU Qing,et al.Application of Chinese remainder theorem in phase unwrapping for dual-baseline In SAR[J].Acta Geodaetica Et Cartographica Sinica,2011,40(6):770-777.
[15]于瀚雯,保铮.利用L1范数的多基线In SAR相位解缠绕技术[J].西安电子科技大学学报(自然科学版),2013,40(4):37-41.doi:10.3969/j.issn.1001-2400.2013.04.006.YU Hanwen and BAO Zheng.L1-norm method for multibaseline In SAR phase unwrapping[J].Journal of Xidian University,2013,40(4):37-41.doi:10.3969/j.issn.1001-2400.2013.04.006.
[16]索志勇,李真芳,吴建新,等.联合像素多基线In SAR降维处理方法[J].西安电子科技大学学报(自然科学版),2009,36(3):385-389.doi:10.3969/j.issn.1001-2400.2009.03.001.SUO Zhiyong,LI Zhenfang,WU Jianxin,et al.Reduceddimension method for joint-pixel multi-baseline In SARprocessing[J].Journal of Xidian University,2009,36(3):385-389.doi:10.3969/j.issn.1001-2400.2009.03.001.
[17]YUAN Z,DENG Y,LI F,et al.Multichannel In SAR DEMreconstruction through improved closed-form robust Chinese remainder theorem[J].IEEE Geoscience and Remote Sensing,2013,10(6):1314-1318.doi:10.1109/LGRS.2013.2238886.
[18]刘会涛,邢孟道,保铮.利用L∞+L1范数的多基线相位解缠绕方法[J].电子与信息学报,2015,37(5):1111-1115.doi:10.11999/JEIT140955.LIU Huitao,XING Mengdao,and BAO Zheng.L∞+L1-norm method for multi-baseline phase unwrapping[J].Journal of Electronics&Information Technology,2015,37(5):1111-1115.doi:10.11999/JEIT140955.
[19]XIE X M.Enhanced multi-baseline unscented kalman filtering phase unwrapping algorithm[J].Journal of Systems Engineering and Electronics,2016,27(2):343-351.doi:10.1109/JSEE.2016.00035.
[20]YU H W and LAN Y.Robust two-dimensional phase unwrapping for multibaseline SAR interferograms:A twostage programming approach[J].IEEE Transactions on Geoscience and Remote Sensing,2016,54(9):1-9.doi:10.1109/TGRS.2016.2558541.
[21]LEE J S,HOPPEL K W,MANGO S,et al.Intensity and phase statistics of multilook polarimetric and interferometric SAR imagery[J].IEEE Transactions on Geoscience and Remote Sensing,1994,32(5):1017-1028.doi:10.1109/36.312890.
[22]DELEDALLE C A,DENIS L,FERRAIOLI G,et al.Combining patch-based estimation and total variation regularization for 3D In SAR reconstruction[C].Proceedings of IEEE International Geoscience and Remote Sensing Symposium,Milan,Italy,2015:2485-2488.doi:10.1109/IGARSS.2015.7326314.
[23]ISHIKAWA H.Exact optimization for markov random fields with convex priors[J].IEEE Transactions on Pattern Analysis and Machine Intelligence,2003,25(10):1333-1336.doi:10.1.1.106.8689.