结合方向聚集性的SAR图像相干斑抑制方法技术

本发明专利技术公开了一种结合图像方向聚集性的SAR图像相干斑抑制方法，主要解决现有技术的细节信息丢失的问题。其相干斑抑制过程为：(1)利用Primal?Sketch稀疏表示模型中的稀疏编码方法，提取SAR图像由线段组成的边脊草图来描述图像的奇异信息；(2)根据边脊草图，将SAR图像划分为均匀区域和非均匀区域；(3)对非均匀区域，采用结合图像结构方向信息和聚集性的SAR图像相干斑抑制方法进行处理；(4)对均匀区域，采用窗口为5*5的Lee滤波方法进行处理；(5)将处理后非均匀区域和均匀区域进行合并，得到SAR图像的相干斑抑制结果。本发明专利技术实现了SAR图像良好的相干斑抑制效果，且有效地保持了图像的细节信息，可用于SAR图像相干斑抑制。

【技术实现步骤摘要】

本专利技术属于图像处理
，涉及SAR图像相干斑抑制方法，可用于SAR图像相干斑抑制。
技术介绍
合成孔径雷达SAR是一种高分辨率成像雷达，属于主动式遥感系统。SAR不受光照，天气等条件的影响，具有全天候的成像能力，在测绘、资源环境监测、地质、地震、火山监测、海洋、军事等领域得到了广泛的应用。SAR图像是通过孔径合成的方式来获得一个较大的雷达天线，利用该大天线发出的电磁波经目标散射元后向散射所成的图像。在SAR的成像过程中，因其利用相干方式来实现孔径合成，故往往会伴随着后向散射电磁波之间的相干现象。这使得在SAR图像中总会出现由于相干而形成的斑点信息。这些相干斑信息是由不同散射源或具有不同传播路径的散射电磁波相干而形成的。从一定意义上来说，这些相干斑信息反映了散射源的一些特性，但同时也给SAR图像的理解和解译带来了很大的困难。在SAR图像中，奇异性信息的保持对于图像信息的理解和解译是非常重要的。在空域的SAR图像的相干斑抑制方法中，有的方法是直接把数字图像处理的方法应用到SAR图像的斑点抑制中，并没有考虑到SAR图像噪声的统计特性，如均值滤波和中值滤波方法；有的方法是以SAR图像斑点噪声满足完全发展相干斑的假设条件，利用局域统计特性对中心像素的真实值进行估计，如Lee滤波和Frost滤波方法，这类方法由于没有考虑到图像不同区域不同特性，因而导致图像细节信息的模糊或丢失。针对这一问题，有学者提出了增强滤波的SAR图像相干斑抑制方法，如增强Lee滤波方法，这类方法利用局域信息将图像进行划分，并针对不同区域采用不同的相干斑抑制方法，能够较好地实现图像奇异性信息保持与斑点噪声抑制的平衡。但是，上述方法都采用局域统计的方法对中心像素进行估计，并不能充分反映图像的结构信息，这在很大程度上会导致一些细节信息丢失。
技术实现思路
本专利技术的目的在于克服上述已有方法的不足，提出了一种结合方向聚集性的SAR 图像相干斑抑制方法，以提升SAR图像相干斑抑制的效果。实现本专利技术的技术方案是利用Primal Sketch稀疏表示模型中稀疏编码的方法提取SAR图像边脊草图，并根据该边脊草图把SAR图像划分为均勻和非均勻区域，对均勻区域采用Lee滤波方法处理，对非均勻区域采用结合了图像结构方向信息和聚集性的SAR图像相干斑抑制方法处理，能够有效地保持和体现SAR图像中一些奇异性区域存在的细节信息，提高了 SAR图像的相干斑抑制效果。具体步骤如下(1)利用I^rimal Sketch稀疏表示模型中的稀疏编码方法，提取SAR图像由线段组成的边脊草图来描述图像的奇异信息；(2)根据边脊草图，将SAR图像划分为非均勻区域和均勻区域；以边脊草图中所标记的线段为中心，将其沿垂直方向两边扩展成具有2*d+l个像素宽的区域，所有这些区4域共同组成SAR图像中含奇异信息的非均勻区域，剩余的区域即为不含奇异信息的均勻区域，其中d取值为1 3;(3)对非均勻区域，采用结合图像结构方向信息和聚集性的SAR图像相干斑抑制方法进行处理(3a)将边脊草图中的线段作为单位，以线段上的每一点为中心，k为半径在线段上选择邻域，计算SAR图像中与该邻域相应的像素的均值μ，再依次计算SAR图像中该邻域内每个像素的值与该均值μ的平均值m，并将平均值m作为该像素的估计值，其中k取值为 1 3 ；(3b)将步骤(3a)中的线段沿其垂直方向的两边逐渐平移d个像素，得到2*d条线段，对这2*d条线段按步骤(3a)进行处理，得到与这2*d条线段相应像素的估计值；(3c)重复步骤(3a)-C3b)，计算出边脊草图中的所有线段，即得到SAR图像非均勻区域的估计值；(4)对于均勻区域，采用窗口为5*5的Lee滤波方法进行处理；(5)将处理后非均勻区域和均勻区域进行合并，得到SAR图像最终的相干斑抑制结果。本专利技术与现有技术相比具有如下优点本专利技术由于利用I^rimal Sketch稀疏表示模型中的稀疏编码方法提取SAR图像的边脊草图，并根据该边脊草图把SAR图像划分为均勻区域和非均勻区域，并分别对均勻区域和非均勻区域采用不同方法进行相干斑抑制处理对均勻区域采用Lee滤波方法进行处理，可较好地实现相干斑抑制和奇异信息保持之间的平衡；对非均勻区域，采用结合图像结构方向信息和聚集性的SAR相干斑抑制方法，该方法简单易行，操作方便，充分利用了局部非均勻区域的结构信息和方向聚集性，尤其能够有效地保持和体现SAR图像中非均勻区域存在的一些细节信息；因此，本专利技术提高了 SAR图像的相干斑抑制效果。附图说明图1是本专利技术对SAR图像相干斑的抑制流程图；图2是本专利技术中基于I^rimal Sketch稀疏表示模型提取SAR图像边脊草图的结果图；图3是用本专利技术与现有方法对DRA SAR X波段分辨率为3米的Bedfordshire图像相干斑抑制结果图。具体实施例方式参照图1，本专利技术的实施过程如下步骤1，利用Primal Sketch稀疏表示模型中的稀疏编码方法，提取输入图像如图 2(a)由线段组成的边脊草图来描述图像线状奇异信息。(1.1)将如图2(a)所示的输入图像与N个尺度和M个方向上的扩展高斯一阶导滤波器和高斯二阶导滤波器进行卷积，以检测图像的边和脊，其中N取值为3 5，M取值为 18 ；(1. 2)计算输入图像每个像素的联合响应，得到联合响应最大值图，并对联合响应最大值图用非极大抑制方法处理，将处理的结果作为建议草图S ；(1. 3)根据建议草图S中最大联合响应的位置，把建议草图S中与该位置连通的点连接成线段，生成一个原始边/脊模型；(1.4)在边/脊模型中添加线段，评价图像的编码长度增益八L，若AL< ε，ε是阈值取值为25，则拒绝接受该线段；否则在边/脊模型中添加该线段，并搜索满足该线段末端与其余像素在平均拟合误差内的分割线作为建议线段，若搜索到该线段转至步骤(1. 4)， 否则结束，即得到了如图2(b)所示的边脊草图。步骤2，根据边脊草图，将SAR图像划分为均勻区域和非均勻区域。以边脊草图中所标记的线段为中心，将其沿垂直方向两边扩展成具有2*d+l个像素宽的区域，所有这些区域共同组成SAR图像中含奇异信息的非均勻区域，剩余的区域即为不含奇异信息的均勻区域，其中d取值为1 3 ；步骤3，对非均勻区域，采用提出的结合图像结构方向信息和聚集性的SAR图像相干斑抑制方法。(3a)将边脊草图中的线段作为单位，以线段上的每一点为中心，k为半径在线段上选择邻域，计算SAR图像中与该邻域相应的像素的均值μ，再依次计算SAR图像中该邻域内每个像素的值与该均值μ的平均值m，并将平均值m作为该像素的估计值，其中k取值为 1 3 ；(3b)将步骤(3a)中的线段沿其垂直方向的两边逐渐平移d个像素，得到2*d条线段，对这2*d条线段按步骤(3a)进行处理，得到与这2*d条线段相应像素的估计值；(3c)重复步骤(3a)-C3b)，计算出边脊草图中的所有线段，即得到SAR图像非均勻区域的估计值。步骤4，对于均勻区域，采用窗口为5*5的Lee滤波方法进行处理。对于非均勻区域，用窗口为5*5的Lee滤波方法对邻域中心为进行估计，公式如下R = R + k*(I-R)其中，本文档来自技高网...

【技术保护点】
１．一种结合方向聚集性的ＳＡＲ图像相干斑抑制方法，包括如下步骤：（１）利用Ｐｒｉｍａｌ　Ｓｋｅｔｃｈ稀疏表示模型中的稀疏编码方法，提取ＳＡＲ图像由线段组成的边脊草图来描述图像的奇异信息；（２）根据边脊草图，将ＳＡＲ图像划分为非均匀区域和均匀区域；以边脊草图中所标记的线段为中心，将其沿垂直方向两边扩展成具有２＊ｄ＋１个像素宽的区域，所有这些区域共同组成ＳＡＲ图像中含奇异信息的非均匀区域，剩余的区域即为不含奇异信息的均匀区域，其中ｄ取值为１～３；（３）对非均匀区域，采用结合图像结构方向信息和聚集性的ＳＡＲ图像相干斑抑制方法进行处理：（３ａ）将边脊草图中的线段作为单位，以线段上的每一点为中心，ｋ为半径在线段上选择邻域，计算ＳＡＲ图像中与该邻域相应的像素的均值μ，再依次计算ＳＡＲ图像中该邻域内每个像素的值与该均值μ的平均值ｍ，并将平均值ｍ作为该像素的估计值，其中ｋ取值为１～３；（３ｂ）将步骤（３ａ）中的线段沿其垂直方向的两边逐渐平移ｄ个像素，得到２＊ｄ条线段，对这２＊ｄ条线段按步骤（３ａ）进行处理，得到与这２＊ｄ条线段相应像素的估计值；（３ｃ）重复步骤（３ａ）－（３ｂ），计算出边脊草图中的所有线段，即得到ＳＡＲ图像非均匀区域的估计值；（４）对于均匀区域，采用窗口为５＊５的Ｌｅｅ滤波方法进行处理；（５）将处理后非均匀区域和均匀区域进行合并，得到ＳＡＲ图像最终的相干斑抑制结果。...

【技术特征摘要】
1.一种结合方向聚集性的SAR图像相干斑抑制方法，包括如下步骤(1)利用Sketch稀疏表示模型中的稀疏编码方法，提取SAR图像由线段组成的边脊草图来描述图像的奇异信息；(2)根据边脊草图，将SAR图像划分为非均勻区域和均勻区域；以边脊草图中所标记的线段为中心，将其沿垂直方向两边扩展成具有2*d+l个像素宽的区域，所有这些区域共同组成SAR图像中含奇异信息的非均勻区域，剩余的区域即为不含奇异信息的均勻区域，其中d取值为1 3 ；(3)对非均勻区域，采用结合图像结构方向信息和聚集性的SAR图像相干斑抑制方法进行处理(3a)将边脊草图中的线段作为单位，以线段上的每一点为中心，k为半径在线段上选择邻域，计算SAR图像中与该邻域相应的像素的均值μ，再依次计算SAR图像中该邻域内每个像素的值与该均值μ的平均值m，并将平均值m作为该像素的估计值，其中k取值为1 3 ；(3b)将步骤(3a)中的线段沿其垂直方向的两边逐渐平移d个像素，得到2*d条线段， 对这2*d条线段按步骤(3a)进行处理，得到与这2*d条线段相应像素的估计值；(3c)重复步骤(3a)-C3b)，计算出边脊草图中的所有线段，即得到SAR图像非均勻区域的估计值；(4)对于均勻区域，采用窗口为5*5的Lee滤波方法进行处理；(5)将处理后非均勻区域和均勻区域进行合并，得到SAR图像最终的相干斑抑制结果。2.根据权利要求1所述的...

【专利技术属性】
技术研发人员：刘芳，焦李成，王爽，戚玉涛，武杰，侯彪，修冬阳，郝红侠，马文萍，蒋文梅，尚荣华，
申请(专利权)人：西安电子科技大学，
类型：发明
国别省市：87