基于Gamma滤波的SAR图像超像素分割方法技术

本发明专利技术提出了一种基于Gamma滤波的SAR图像超像素分割方法，用于解决现有基于滤波的超像素分割方法中存在的超像素分割结果准确度低的技术问题，实现步骤为：采用均值滤波对输入的原始SAR图像进行滤波，得到像素灰度的平均值和标准差；分别计算像素灰度平均值的方差系数、相干斑噪声的方差系数和SAR图像的变化系数；判断各系数的关系，确定是否进行Gamma滤波；采用Gamma滤波方法对输入的原始SAR图像进行相干斑噪声降噪；对降噪后的SAR图像进行超像素分割，得到多个超像素块并输出。本发明专利技术降低了SAR图像中相干斑噪声的影响，提高了SAR图像超像素分割的准确度，可用于对SAR图像的目标检测、识别及分类。

【技术实现步骤摘要】
基于Gamma滤波的SAR图像超像素分割方法
本专利技术属于图像处理
，涉及一种SAR图像超像素分割方法，具体涉及一种基于Gamma滤波的SAR图像超像素分割方法，可用于对SAR图像的目标检测、识别及分类。
技术介绍
合成孔径雷达(Syntheticapertureradar,SAR)是一种工作在微波波段的相干成像雷达，它以其高分辨率和全天候、全天时、大面积的数据获取能力而成为当前遥感观测的重要手段，在资源、环境、考古以及军事等方面得到广泛的应用。地面目标对雷达发射电磁波的后向散射形成SAR图像。随着SAR图像的大量获取，智能的图像理解与解译技术成为当今的研究热点。作为图像理解与解译的关键性工作，图像的超像素分割是将一副输入的SAR图像按照不同的属性、内容、特征等将相邻的像素划分为图像块。对一幅图像进行有效而精确的超像素分割，能够提高进一步对图像理解的准确度，获取更多的有效信息，并很大程度上降低后续SAR图像处理任务的复杂度。然而SAR图像固有的相干成像机理使SAR图像存在无法消除的相干斑噪声，导致每个像素的灰度值发生随机的变化，从而影响SAR图像的超像素分割效果及后续的SAR图像处理结果。目前，大多数学者主要采用以下三种方法来解决上述问题：1、基于统计模型的超像素分割方法；2、基于尺度的超像素分割方法。3、基于滤波的超像素分割方法。第一种方法是对SAR图像中的相干斑噪声进行概率统计，假设概率分布符合Rayleigh分布，Gamma分布，K分布等，将概率统计模型融合到常规的超像素分割算法的模型中，从而得到适于SAR图像的超像素分割结果，该方法的不足之处是，虽然利用了SAR图像的相干斑噪声信息来进行SAR图像的超像素分割，但是其基于统计信息来进行计算，从而增加了计算过程的复杂度，降低了算法的运行效率。第二种方法的核心思想是将SAR图像看作不同尺度的纹理，根据不同类别的目标呈现不同纹理这一性质，将对SAR图像的超像素分割，转化为对SAR图像中不同纹理的识别。但是该类方法实现的过程复杂，得到的结果相比于其他的超像素分割方法，并没有明显的优势。第三种方法是首先采用滤波器对SAR图像中的相干斑噪声进行滤除，降低相干斑噪声对SAR图像处理的影响，从而得到降噪后的SAR图像，然后对降噪后的SAR图像进行超像素分割。该方法实现简单，在复杂度比前两种方法的复杂度小的情况下得到的超像素分割的准确度相同。在此类方法中，关键是对SAR图像进行相干斑噪声的滤除，而此技术已经非常成熟，其中包括有采用中值滤波器，局域滤波器，Lee滤波器，Sigma滤波器，Frost滤波器等对SAR图像进行预处理去除相干斑噪声。但这些滤波器对SAR图像的相干斑噪声进行滤除时，虽然都能够滤除相干斑噪声，但是图像中的边缘信息、纹理信息及线性特征信息等不能得到很好的保持。
技术实现思路
本专利技术的目的在于克服上述现有技术存在的不足，提出了一种基于Gamma滤波的SAR图像超像素分割方法，在降低SAR图像中相干斑噪声对其干扰的情况下，有效保持了图像的边缘信息、纹理信息以及线性特征信息，用于解决现有基于滤波的超像素分割方法中存在的超像素分割准确度低的技术问题。为实现上述目的，本专利技术采取的技术方案包括如下步骤：(1)输入原始SAR图像；(2)采用均值滤波对输入的原始SAR图像进行滤波，得到原始SAR图像像素灰度的平均值和标准差σ(i,j)；(3)利用原始SAR图像像素灰度的平均值和标准差σ(i,j)，计算原始SAR图像像素灰度平均值的方差系数CI：(4)计算原始SAR图像中相干斑噪声的方差系数Cu：其中，L表示原始SAR图像的视数；(5)计算原始SAR图像的变化系数Cmax:(6)当Cu≤CI≤Cmax时，采用Gamma滤波方法，对输入的原始SAR图像进行相干斑噪声降噪，得到降噪后的SAR图像IGam(i,j)：其中，α为异质参数，且(7)对降噪后的SAR图像IGam(i,j)进行超像素分割，得到多个超像素块并输出。本专利技术与现有技术相比，具有以下优点：本专利技术由于在对SAR图像进行超像素分割之前首先采用Gamma滤波滤除原始SAR图像中的相干斑噪声，得到降噪后的SAR图像再对其进行超像素分割，与现有技术相比，通过本专利技术能够在降低了相干斑噪声对SAR图像处理的影响时，有效的保持了图像的边缘信息、纹理信息以及线性特征信息，提高了超像素分割的准确度。附图说明图1为本专利技术的实现流程框图；图2为本专利技术模拟SAR图像产生的超像素图像；图3为本专利技术含有纹理的模拟SAR图像产生的超像素图像；图4为本专利技术真实SAR图像产生的超像素图像。具体实施方式下面结合附图和实施例，对本专利技术作进一步的详细描述。实施例1参照图1：本专利技术包括如下步骤：步骤1：输入原始SAR图像。原始SAR图像分为模拟SAR图像和真实SAR图像两大类，在本实施例中，原始SAR图像采用不含纹理的模拟SAR图像。步骤2：采用均值滤波对输入的原始SAR图像进行滤波，得到原始SAR图像像素灰度的平均值和标准差σ(i,j)；均值滤波是典型的线性滤波算法，在进行均值滤波时，需要一个滤波模板，用模板中的全体像素的平均值来代替原来的像素值。对于滤波模板的选择，一般包括有3*3的模板，5*5的模板，或者更大窗口的模板，在本专利技术中对于模板的选择不受限定，但是由于窗口越大，对于滤波后的结果会有所影响，在本专利技术的实施例中选择3*3的权系数矩阵模板，其表达式如下：其中，s为模板矩阵表达式的行，t为模板矩阵表达式的列；(2a)采用权系数矩阵模板w(s,t)，对输入的原始SAR图像进行滤波，得到原始SAR图像像素灰度的平均值(2b)计算原始SAR图像像素灰度的标准差σ(i,j)：其中，N表示输入的原始SAR图像的横坐标上的像素个数，M表示输入的原始SAR图像的纵坐标上的像素个数。步骤3：利用原始SAR图像像素灰度的平均值和标准差σ(i,j)，计算原始SAR图像像素灰度平均值的方差系数CI：步骤4：计算原始SAR图像中相干斑噪声的方差系数Cu：其中，L表示原始SAR图像的视数；步骤5：计算原始SAR图像的变化系数Cmax:步骤6：判断原始SAR图像像素灰度平均值的方差系数CI与原始SAR图像中相干斑噪声的方差系数Cu及原始SAR图像变化系数Cmax的关系，当Cu≤CI≤Cmax时，采用Gamma滤波方法，对输入的原始SAR图像进行相干斑噪声降噪，得到降噪后的SAR图像IGam(i,j)：其中，α为异质参数，且当CI＜Cu时，则SAR图像的像素灰度值为原始SAR图像像素灰度的平均值当CI＞Cmax时，则SAR图像的像素灰度值为原始SAR图像的像素灰度值I(i,j)。而对于这两种情况，不对其进行Gamma滤波。步骤7：对降噪后的SAR图像IGam(i,j)进行超像素分割，得到多个超像素块并输出。本步骤的具体实现如下：(7a)计算降噪后的SAR图像IGam(i,j)水平方向上相邻像素的像素值的差W(i,j)：W(i,j)＝|IGam(i,j)-IGam(i-1,j)|其中，IGam(i-1,j)表示像素值IGam(i,j)水平方向上相邻像素的像素值；(7b)计算降噪后的SAR图像IGam(i,j)垂直方向上相邻像素的像素值的差V(i,j)：V本文档来自技高网...

【技术保护点】
一种基于Gamma滤波的SAR图像超像素分割方法，包括如下步骤：(1)输入原始SAR图像；(2)采用均值滤波对输入的原始SAR图像进行滤波，得到原始SAR图像像素灰度的平均值

【技术特征摘要】
1.一种基于Gamma滤波的SAR图像超像素分割方法，包括如下步骤：(1)输入原始SAR图像；(2)采用均值滤波对输入的原始SAR图像进行滤波，得到原始SAR图像像素灰度的平均值和标准差σ(i,j)；(3)利用原始SAR图像像素灰度的平均值和标准差σ(i,j)，计算原始SAR图像像素灰度平均值的方差系数CI：(4)计算原始SAR图像中相干斑噪声的方差系数Cu：其中，L表示原始SAR图像的视数；(5)计算原始SAR图像的变化系数Cmax:(6)当Cu≤CI≤Cmax时，采用Gamma滤波方法，对输入的原始SAR图像进行相干斑噪声降噪，得到降噪后的SAR图像IGam(i,j)：其中，α为异质参数，且(7)对降噪后的SAR图像IGam(i,j)进行超像素分割，得到多个超像素块并输出。2.根据权利要求1所述的基于Gamma滤波的SAR图像超像素分割方法，其特征在于：步骤(2)中所述的采用均值滤波对输入的原始SAR图像进行滤波，实现步骤如下：(2a)采用权系数矩阵模板w(s,t)，对输入的原始SAR图像进行滤波，得到原始SAR图像像素灰度的平均值其中，s为模板矩阵表达式的行，t为模板矩阵表达式的列；(2b)计算原始SAR图像像素灰度...

【专利技术属性】
技术研发人员：冯冬竹，余航，袁晓光，戴浩，范琳琳，高飞飞，孙景荣，
申请(专利权)人：西安电子科技大学，
类型：发明
国别省市：陕西,61