一种基于SAR图像的山区地表微形变信息提取方法技术

本发明专利技术公开了一种基于SAR图像的山区地表微形变信息提取方法，通过对SAR影像组进行精配准并提取振幅与相位信息，然后根据像元点的振幅离差指数阈值法和归一化同极化相位差阈值法选取永久散射点；同时利用SAR影像组和高程数据进行差分干涉处理得到差分相位图；再利用二次滤波算法选取分布式散射点，并对分布式散射点的差分相位值进行优化；再将永久散射点和分布式散射点合并为稳定散射点，构造稳定散射点的Delaunay三角网；最后利用网络平差法对稳定散射点的差分相位进行求解，最终得到地表微形变信息。

【技术实现步骤摘要】

本专利技术属于图像处理
，更为具体地讲，涉及一种基于SAR图像的山区地表微形变信息提取方法。
技术介绍
永久散射体差分干涉测量技术(PS-InSAR)是一种利用时间序列SAR影像组提取地表形变信息的技术。PS-InSAR技术能够全天时全天候对地表进行监测而不受天气影响，因而为大范围的地表形变监测提供了快速、有效且经济的手段。中国地域辽阔，自然环境复杂，PS-InSAR技术能够在观测条件有限的情况下实现地表信息提取。相对于传统的位移计、水准仪等实地测量，该技术能够极大的提高覆盖范围、降低成本、提高空间分辨率。PS-InSAR提取地表形变信息的步骤包括：配准、差分干涉处理、永久散射点的提取、相位分解等。由于处理过程环环相扣，任意一个环节的误差都将直接被带入下一步处理中，因此，需要每一步都达到较高的精度，才能保证最后结果的准确性。PS-InSAR技术是一种对永久散射体进行差分干涉处理，从而提取出整幅影像地表形变信息的技术。永久散射体指的是研究区域中具有高强度稳定散射特性的地物目标，可能是人工建筑、裸露岩石等。在SAR图像上，永久散射体所对应的像元具有较高的亮度和信噪比，在时间轴上呈现出很好的稳定性。可以通过对时间序列SAR影像组进行分析从而选取出永久散射点，永久散射点选取算法的准确度和精度将直接影响PS-InSAR的技术的精度。目前关于永久散射点选取算法虽然较多，但存在着许多不足之处。传统的永久散射点选取算法受到地表类型限制，人工建筑、裸露岩石砂砾等地面目标具有较强的相关性因而能够被选取为永久散射点，然而植被覆盖较多的山区往往提取不到足够的永久散射点。同时算法对时间序列影像的数量要求高，传统PS-InSAR算法往往需要十几甚至二十几幅图像才能保证永久散射点的可靠性，因此，必须对现有的技术进行改进，以实现提高山区稳定散射
点的数量，同时降低对SAR图像数量的要求。
技术实现思路
本专利技术的目的在于克服现有技术的不足，提供一种基于SAR图像的山区地表微形变信息提取方法，通过提取稳定散射点，以实现山区的地表微形变信息提取。为实现上述专利技术目的，本专利技术为一种基于SAR图像的山区地表微形变信息提取方法，其特征在于，包括以下步骤：(1)、对时间序列SAR影像组进行预处理在预先采集的N幅时间序列SAR影像组中，根据综合相关函数选取一幅SAR图像作为公共主影像，其余的为从影像；再利用sarscape软件对每幅从影像按主影像进行配准，并提取每幅SAR图像的振幅与相位信息，再为每幅SAR图像各生成一幅振幅图与相位图；(2)、差分干涉将时间序列SAR影像组进行干涉处理并去平，得到干涉相位图；同时将外部高程数据转化为高程相位分量；然后将干涉相位图的相位与高程相位分量进行差分处理，得到差分干涉相位图；(3)、利用振幅离差阈值法与归一化同极化相位差阈值法选取永久散射点(3.1)、根据振幅图，计算每个像元点的振幅离差指数DA(i,j)： D A ( i , j ) = σ A ( i , j ) m A ( i , j ) ]]>其中，(i,j)表示像元点的坐标，mA(i,j)、σA(i,j)分别为振幅的均值与方差： m A ( i , j ) = 1 N Σ k = 1 N A k ( i , j ) , σ A = 1 N Σ k = 1 N ( A k ( i , j ) - m A ( i , j ) ) 2 ]]>其中，Ak(i,j)表示第k幅振幅图的振幅值，N表示振幅图或相位图的总个数；(3.2)、根据相位图，计算每个像元点的归一化同极化相位差χ(i,j)：其中分别为第k相位图的HH极化与VV极化相位值；(3.3)、根据振幅离差指数DA(i,j)和归一化同极化相位差χ(i,j)选取永久散射点；(4)、利用双重自适应性滤波算法选取分布式散射点(4.1)、选取中心像元点根据振幅图，计算每个像元点的同极化振幅幅比ρ(i,j)： ρ ( i , j ) = 1 N Σ k = 1 N A k H H ( i , j ) A k V V ( i 本文档来自技高网...

【技术保护点】
一种基于SAR图像的山区地表微形变信息提取方法，其特征在于，包括以下步骤：(1)、对时间序列SAR影像组进行预处理在预先采集的N幅时间序列SAR影像组中，根据综合相关函数选取一幅SAR图像作为公共主影像，其余的为从影像；再利用sarscape软件对每幅从影像按主影像进行配准，并提取每幅SAR图像的振幅与相位信息，再为幅SAR图像各生成一幅振幅图与相位图；(2)、差分干涉将时间序列SAR影像组进行干涉处理并去平，得到干涉相位图；同时将外部高程数据转化为高程相位分量；然后将干涉相位图的相位与高程相位分量进行差分处理，得到差分干涉相位图；(3)、利用振幅离差阈值法与归一化同极化相位差阈值法选取永久散射点(3.1)、根据振幅图，计算每个像元点的振幅离差指数DA(i,j)：DA(i,j)=σA(i,j)mA(i,j)]]>其中，(i,j)表示像元点的坐标，mA(i,j)、σA(i,j)分别为振幅的均值与方差：mA(i,j)=1NΣk=1NAk(i,j),σA=1NΣk=1N(Ak(i,j)-mA(i,j))2]]>其中，Ak(i,j)表示第k幅振幅图的振幅值，N表示振幅图或相位图的总个数；(3.2)、根据振幅图，计算每个像元点的归一化同极化相位差χ(i,j)：其中分别为第k相位图的HH极化与VV极化相位值；(3.3)、根据振幅离差指数DA(i,j)和归一化同极化相位差χ(i,j)选取永久散射点；(4)、利用双重自适应性滤波算法选取分布式散射点(4.1)、选取中心像元点根据振幅图，计算每个像元点的同极化振幅比ρ(i,j)：ρ(i,j)=1NΣk=1NAkHH(i,j)AkVV(i,j)]]>其中，表示第k幅振幅图在HH极化振幅值，表示第k幅振幅图在VV极化振幅值；设定振幅离差指数阈值阈值D'T、同极化振幅比阈值将振幅离差指数满足DA(i,j)＜D'T且同极化振幅比值满足的像元点设置为中心像元点；(4.2)、利用统计学中的KS检验以中心像元点为中心、滤波窗口内，判定中心像元点的同质点；(4.3)、确定同质点族：在滤波窗口内，舍去没有与中心像元点直接或间接相连的同质点，将中心像元点与其直接或间接相连的同质点作为一个同质点族；(4.4)、判断同质点族大小是否超过阈值T，如果同质点族大小超过阈值T，则将该同质点归入一次候选分布式散射点集CDS1中；(4.5)、计算同质点族的振幅均值，并作为CDS1的振幅值；(4.6)、对CDS1重复步骤(4.2)～(4.4)，再次判定同质点族，当同质点族大小超过阈值T时，归入二次候选分布式散射点集CDS2中；(4.7)、根据差分干涉相位图计算出CDS2同质点族的差分相位均值，并作为CDS2的差分相位值；同时计算CDS2中各像元点的振幅离差指数，将振幅离差指数DA(i,j)小于预设阈值D'T'的像元点设置为分布式散射点；(5)、形变信息提取将永久散射点与分布式散射点合并为稳定散射点集，建立Delaunay三角网，利用网络平差法对稳定散射点的差分相位进行相位分解，从而提取出山区地表微形变信息。...

【技术特征摘要】
1.一种基于SAR图像的山区地表微形变信息提取方法，其特征在于，包括以下步骤：(1)、对时间序列SAR影像组进行预处理在预先采集的N幅时间序列SAR影像组中，根据综合相关函数选取一幅SAR图像作为公共主影像，其余的为从影像；再利用sarscape软件对每幅从影像按主影像进行配准，并提取每幅SAR图像的振幅与相位信息，再为幅SAR图像各生成一幅振幅图与相位图；(2)、差分干涉将时间序列SAR影像组进行干涉处理并去平，得到干涉相位图；同时将外部高程数据转化为高程相位分量；然后将干涉相位图的相位与高程相位分量进行差分处理，得到差分干涉相位图；(3)、利用振幅离差阈值法与归一化同极化相位差阈值法选取永久散射点(3.1)、根据振幅图，计算每个像元点的振幅离差指数DA(i,j)： D A ( i , j ) = σ A ( i , j ) m A ( i , j ) ]]>其中，(i,j)表示像元点的坐标，mA(i,j)、σA(i,j)分别为振幅的均值与方差： m A ( i , j ) = 1 N Σ k = 1 N A k ( i , j ) , σ A = 1 N Σ k = 1 N ( A k ( i , j ) - m A ( i , j ) ) 2 ...

【专利技术属性】
技术研发人员：陈彦，陈浩，童玲，
申请(专利权)人：电子科技大学，
类型：发明
国别省市：四川;51