一种基于GVF-Snake模型边界探测的相位解缠算法
|
摘要
针对矿区沉降监测中低相干和强噪声下出现相位不连续、条纹间断等导致的InSAR(Interferometric SAR)监测可靠性不高等问题,提出了一种基于GVF-Snake(gradient vector flow snake)模型干涉条纹边界探测的相位解缠方法.该算法利用GVF-Snake模型对沉降漏斗区的干涉条纹进行边缘检测获取干涉条纹的边界,对检测边界内的像元按照相关准则进行解缠绕,通过改进的中值滤波方法消除孤立点.以济宁矿区为实验区,采用日本ALOS PALSAR真实雷达数据进行矿区沉降监测以验证新解缠算法的有效性,利用两组实验数据分别采用基于GVF-Snake模型的边界探测解缠算法与最小费用流等其它5种解缠算法进行对比实验与分析.结果表明:所提算法在干涉条纹质量较好的情况下与其它解缠算法高度一致,以MCF(minimum cost flow)最小费用流解缠方法作为评价标准,新方法解缠精度可以达到±0.005 8rad;在干涉条纹间断、噪声影响严重的情况下,可以实现InSAR矿区沉降漏斗的正确相位解缠,从而保证InSAR矿区大量级沉降监测的可靠性及精度.
In view of the low reliability caused by phase discontinuity and fringe disconnection under the low coherence and strong noises when monitoring the mining subsidence using InSAR technique,a new phase unwrapping algorithm based on the GVF-Snake model of edge detection was proposed.In the new algorithm,the boundary of the interferometric fringes in the subsidence funnel area was firstly obtained through edge detection using GVF-Snake model.Then the pixels surrounded by the detected boundary were unwrapped according to some rules.Finally an improved median filtering method was used to eliminate isolated points.Taking Ji'ning mining area(China)as the experimental plot,two Japanese ALOS PALSAR data for monitoringthe mining subsidence were used to verify the effectiveness and the performance of the new algorithm.Two groups of experimental data derived from the interferogram were used to carry out the comparative experiment and analysis for the new method based on Gradient Vector Flow Snake(GVF-Snake)model of edge detection algorithm and other five kinds of phase unwrapping algorithms,such as minimum cost flow(MCF)method,branch cut method proposed by Goldstein,and so on.The results show that the new algorithm under the condition of a better quality of interferogram ties up with other unwrapping algorithms with high consistency.Taking MCF method as the evaluation standard,the accuracy of new method can reach±0.0058 rad.On the other hand,under the condition of interferogram which is severely affected by the interferometric fringe disconnection and noise,the new method can get correct unwrapping result for the funnel area when monitoring the mining subsidence using InSAR technique.So when monitoring the mining subsidence using InSAR technique,the reliability and accuracy can be ensured.
In view of the low reliability caused by phase discontinuity and fringe disconnection under the low coherence and strong noises when monitoring the mining subsidence using InSAR technique,a new phase unwrapping algorithm based on the GVF-Snake model of edge detection was proposed.In the new algorithm,the boundary of the interferometric fringes in the subsidence funnel area was firstly obtained through edge detection using GVF-Snake model.Then the pixels surrounded by the detected boundary were unwrapped according to some rules.Finally an improved median filtering method was used to eliminate isolated points.Taking Ji'ning mining area(China)as the experimental plot,two Japanese ALOS PALSAR data for monitoringthe mining subsidence were used to verify the effectiveness and the performance of the new algorithm.Two groups of experimental data derived from the interferogram were used to carry out the comparative experiment and analysis for the new method based on Gradient Vector Flow Snake(GVF-Snake)model of edge detection algorithm and other five kinds of phase unwrapping algorithms,such as minimum cost flow(MCF)method,branch cut method proposed by Goldstein,and so on.The results show that the new algorithm under the condition of a better quality of interferogram ties up with other unwrapping algorithms with high consistency.Taking MCF method as the evaluation standard,the accuracy of new method can reach±0.0058 rad.On the other hand,under the condition of interferogram which is severely affected by the interferometric fringe disconnection and noise,the new method can get correct unwrapping result for the funnel area when monitoring the mining subsidence using InSAR technique.So when monitoring the mining subsidence using InSAR technique,the reliability and accuracy can be ensured.
引文
[1]程璞,许才军,王华.InSAR相位解缠算法研究[J].大地测量与地球动力学,2007,27(3):50-55.CHENG Pu,XU Caijun,WANG Hua.Analysis of phase unwrapping algorithm[J].Journal of Geodesy and Geodynamics.2007,27(3):50-55.
[2]FAN Hongdong,GU Wei,QIN Yong,et al.A model for extracting large deformation mining subsidence using D-InSAR technique and probality integral method[J].Transactions of Nonferrous Metals Society of China,2014,24:1242-1247.
[3]王志勇,张继贤,黄国满.基于InSAR的济宁矿区沉降精细化监测与分析[J].中国矿业大学学报,2014,43(1):169-174.WANG Zhiyong,ZHANG Jixian,HUANG Guoman.Precise monitoring and analysis of the land subsidence in Jining coal mining area based on InSAR technique[J].Journal of China University of Mining&Technology,2014,43(1):169-174.
[4]ALEX H M N,GE Linlin,LI Xiaojing.Assessments of land subsidence in the Gippsland Basin of Australia using ALOS PALSAR data[J].Remote Sensing of Environment,2015,159:86-101.
[5]GOLDSTEIN R M,ZEBKER H A,WENER C L.Satellite radar interferometry:Two-dimensional phase unwrapping[J].Radio Science,1988,23(4):713-720.
[6]GHIGLIA D C,ROMERO L A.Robust two-dimensional weighted and unweighted phase unwrapping that uses fast transforms and iterative methods[J].Journal of the Optical Society of America A,1994,11(1):107-117.
[7]COSTANTINI M.A novel phase unwrapping method based on network programming[J].IEEE Transactions on Geoscience and Remote Sensing,1998,36(3):813-821.
[8]于向明,孙学宏,刘丽萍,等.优化枝切线设置的孔径雷达相位解缠方法[J].测绘科学,2015,40(9):43-47.YU Xiangming,SUN Xuehong,LIU Liping,et al.A phase unwrapping method of interferometric SAR based on optimal allocation of branch-cuts[J].Science of Surveying and Mapping,2015,40(9):43-47.
[9]ZHAO M,HUANG L,ZHANG Q,et al.Qualityguided phase unwrapping technique:comparison of quality maps and guiding strategies[J].Applied Optics,2011,50(33):6214-6224.
[10]陈强,杨莹辉,刘国祥,等.基于边界探测的InSAR最小二乘整周相位解缠方法[J].测绘学报,2012,41(3):441-448.CHEN Qiang,YANG Yinghui,LIU Guoxiang,et al.InSAR phase unwrapping using least squares method with integer ambiguity resolution and edge detection[J].Acta Geodaetica et Cartographica Sinica,2012,41(3):441-448.
[11]CHEN C W,ZEBKER H A.Network approaches to two-Dimensional phase unwrapping;intractability and two new algorithms[J].Journal of the Optical Society of America A,2000,17(3):401-414.
[12]赵争,张继贤,张过.遗传算法在InSAR相位解缠中的应用[J].测绘科学,2002,27(3):37-39.ZHAO Zheng,ZHANG Jixian,ZHANG Guo.The application of genetic algorithms in InSAR phase unwrapping[J].Science of Surveying and Mapping,2002,27(3):37-39.
[13]魏志强,金亚秋.基于蚁群算法的InSAR相位解缠算法[J].电子与信息学报,2008,30(3):518-523.WEI Zhiqiang,JIN Yaqiu.InSAR phase unwrapping algorithm based on ant colony algorithm[J].Journal of Electronics&Information Technology,2008,30(3):518-523.
[14]刘国林,郝华东,陶秋香.卡尔曼滤波相位解缠及其与其他方法的对比分析[J].武汉大学学报(信息科学版),2010,35(10):1174-1178.LIU Guolin,HAO Huadong,TAO Qiuxiang.Kalman filter phase unwrapping algorithm and comparison and analysis with other methods[J].Geomatics and Information Science of Wuhan University,2010,35(10):1174-1178.
[15]詹总谦,钱俊,舒宁.基于区域分割和编码的相位解缠方法[J].武汉大学学报(信息科学版),2002,27(3):316-320.ZHAN Zongqian,QIAN Jun,SHU Ning.A method of phase unwrapping based on region segmentation and encoding[J].Geomatics and Information Science of Wuhan University,2002,27(3):316-320.
[16]GE L,CHANG H C,RIZOS C,Mine subsidence monitoring using multi-source satellite SAR images[J].Photogrammetric Engineering and Remote Sensing,2007,73(3):259-266.
[17]ROSEN P A,HENSLEY S,JOUGHIN I R,et al.Synthetic aperture radar interferometry[J].Proceedings of the IEEE,2000,88(3):333-382.
[18]王志勇,孟庆颖.InSAR相位解缠算法在矿区沉降监测中的对比与分析[J].遥感信息,2014,29(6):47-51.WANG Zhiyong,MENG Qingying.Comparison and analysis of phase unwrapping algorithms in mining subsidence monitoring[J].Remote Sensing Information,2014,29(6):47-51.
[19]赵凌君,贾承丽,匡纲要.SAR图像边缘检测方法综述[J].中国图象图形学报,2007,12(12):2042-2049.ZHAO Lingjun,JIA Chengli,KUANG Gangyao.Overview of edge detection in SAR images[J].Journal of Image and Graphics,2007,12(12):2042-2049.
[20]QIN Lunming,ZHU Ce,ZHAO Yao,et al.Generalized gradient vector flow for snakes:new observations,analysis and improvement[J].IEEE Transactions on Circuits and Systems for Video Technology,2013,23(5):883-897.
[21]KASS M,WITKIN A,TERZOPOULOS D.Snakes:active contour models[J].International Journal of Computer Vision,1988,1(4):321-331.
[22]XU Chenyang,PRINCE J L.Gradient vector flow:A new external force model for snakes[C]//IEEE Computer Society Conference on Computer Vision and Pattern Recognition.San Juan:Puerto Rico,1997:66-71.
[2]FAN Hongdong,GU Wei,QIN Yong,et al.A model for extracting large deformation mining subsidence using D-InSAR technique and probality integral method[J].Transactions of Nonferrous Metals Society of China,2014,24:1242-1247.
[3]王志勇,张继贤,黄国满.基于InSAR的济宁矿区沉降精细化监测与分析[J].中国矿业大学学报,2014,43(1):169-174.WANG Zhiyong,ZHANG Jixian,HUANG Guoman.Precise monitoring and analysis of the land subsidence in Jining coal mining area based on InSAR technique[J].Journal of China University of Mining&Technology,2014,43(1):169-174.
[4]ALEX H M N,GE Linlin,LI Xiaojing.Assessments of land subsidence in the Gippsland Basin of Australia using ALOS PALSAR data[J].Remote Sensing of Environment,2015,159:86-101.
[5]GOLDSTEIN R M,ZEBKER H A,WENER C L.Satellite radar interferometry:Two-dimensional phase unwrapping[J].Radio Science,1988,23(4):713-720.
[6]GHIGLIA D C,ROMERO L A.Robust two-dimensional weighted and unweighted phase unwrapping that uses fast transforms and iterative methods[J].Journal of the Optical Society of America A,1994,11(1):107-117.
[7]COSTANTINI M.A novel phase unwrapping method based on network programming[J].IEEE Transactions on Geoscience and Remote Sensing,1998,36(3):813-821.
[8]于向明,孙学宏,刘丽萍,等.优化枝切线设置的孔径雷达相位解缠方法[J].测绘科学,2015,40(9):43-47.YU Xiangming,SUN Xuehong,LIU Liping,et al.A phase unwrapping method of interferometric SAR based on optimal allocation of branch-cuts[J].Science of Surveying and Mapping,2015,40(9):43-47.
[9]ZHAO M,HUANG L,ZHANG Q,et al.Qualityguided phase unwrapping technique:comparison of quality maps and guiding strategies[J].Applied Optics,2011,50(33):6214-6224.
[10]陈强,杨莹辉,刘国祥,等.基于边界探测的InSAR最小二乘整周相位解缠方法[J].测绘学报,2012,41(3):441-448.CHEN Qiang,YANG Yinghui,LIU Guoxiang,et al.InSAR phase unwrapping using least squares method with integer ambiguity resolution and edge detection[J].Acta Geodaetica et Cartographica Sinica,2012,41(3):441-448.
[11]CHEN C W,ZEBKER H A.Network approaches to two-Dimensional phase unwrapping;intractability and two new algorithms[J].Journal of the Optical Society of America A,2000,17(3):401-414.
[12]赵争,张继贤,张过.遗传算法在InSAR相位解缠中的应用[J].测绘科学,2002,27(3):37-39.ZHAO Zheng,ZHANG Jixian,ZHANG Guo.The application of genetic algorithms in InSAR phase unwrapping[J].Science of Surveying and Mapping,2002,27(3):37-39.
[13]魏志强,金亚秋.基于蚁群算法的InSAR相位解缠算法[J].电子与信息学报,2008,30(3):518-523.WEI Zhiqiang,JIN Yaqiu.InSAR phase unwrapping algorithm based on ant colony algorithm[J].Journal of Electronics&Information Technology,2008,30(3):518-523.
[14]刘国林,郝华东,陶秋香.卡尔曼滤波相位解缠及其与其他方法的对比分析[J].武汉大学学报(信息科学版),2010,35(10):1174-1178.LIU Guolin,HAO Huadong,TAO Qiuxiang.Kalman filter phase unwrapping algorithm and comparison and analysis with other methods[J].Geomatics and Information Science of Wuhan University,2010,35(10):1174-1178.
[15]詹总谦,钱俊,舒宁.基于区域分割和编码的相位解缠方法[J].武汉大学学报(信息科学版),2002,27(3):316-320.ZHAN Zongqian,QIAN Jun,SHU Ning.A method of phase unwrapping based on region segmentation and encoding[J].Geomatics and Information Science of Wuhan University,2002,27(3):316-320.
[16]GE L,CHANG H C,RIZOS C,Mine subsidence monitoring using multi-source satellite SAR images[J].Photogrammetric Engineering and Remote Sensing,2007,73(3):259-266.
[17]ROSEN P A,HENSLEY S,JOUGHIN I R,et al.Synthetic aperture radar interferometry[J].Proceedings of the IEEE,2000,88(3):333-382.
[18]王志勇,孟庆颖.InSAR相位解缠算法在矿区沉降监测中的对比与分析[J].遥感信息,2014,29(6):47-51.WANG Zhiyong,MENG Qingying.Comparison and analysis of phase unwrapping algorithms in mining subsidence monitoring[J].Remote Sensing Information,2014,29(6):47-51.
[19]赵凌君,贾承丽,匡纲要.SAR图像边缘检测方法综述[J].中国图象图形学报,2007,12(12):2042-2049.ZHAO Lingjun,JIA Chengli,KUANG Gangyao.Overview of edge detection in SAR images[J].Journal of Image and Graphics,2007,12(12):2042-2049.
[20]QIN Lunming,ZHU Ce,ZHAO Yao,et al.Generalized gradient vector flow for snakes:new observations,analysis and improvement[J].IEEE Transactions on Circuits and Systems for Video Technology,2013,23(5):883-897.
[21]KASS M,WITKIN A,TERZOPOULOS D.Snakes:active contour models[J].International Journal of Computer Vision,1988,1(4):321-331.
[22]XU Chenyang,PRINCE J L.Gradient vector flow:A new external force model for snakes[C]//IEEE Computer Society Conference on Computer Vision and Pattern Recognition.San Juan:Puerto Rico,1997:66-71.