基于SSIM校正聚类稀疏表示的SAR图像降斑方法技术

本发明专利技术公开了一种基于SSIM校正聚类的稀疏表示SAR图像降斑方法，主要解决现有技术中不能保留更多的点目标和纹理细节的问题。其实现步骤为：(1)输入图像；(2)估计输入SAR图像中噪声的方差；(3)对输入图像进行分类；(4)对分类后图像进行初始化字典；(5)更新分类后图像的稀疏表示矩阵和字典；(6)利用更新后的稀疏表示矩阵和字典重构图像，得到降斑后的SAR图像。本发明专利技术不仅能够提高SAR图像中斑点噪声的抑制效果，而且能很好保留图像中点目标和纹理细节信息，可用于SAR图像分割和目标识别的预处理。

【技术实现步骤摘要】
基于SSIM校正聚类稀疏表示的SAR图像降斑方法
本专利技术属于图像处理
，更进一步涉及合成孔径雷达SAR图像的降斑方法，可用于SAR图像分割和目标识别的预处理。
技术介绍
合成孔径雷达SAR具有全天候、全天时、高分辨率和穿透力强等特点，因此SAR在军事方面和民用方面都得到了广泛的应用。由于SAR系统是微波相干成像，因此SAR图像在获取的过程中容易受到斑点噪声的影响。斑点噪声的存在大大降低了SAR图像的分辨率，影响了SAR图像后续的处理与解译，如何有效地抑制SAR图像中的斑点噪声变得尤为重要。近年来，稀疏表示被广泛的应用到图像去噪中。然而，标准的稀疏表示算法假定稀疏矩阵中非零系数的出现是随机的，并没有考虑到图像块之间的相似性。因此，一些学者提出了将图像的几何结构相似性与稀疏表示相结合的算法。其中，董伟生等人提出的基于聚类的稀疏表示CSR算法认为图像块之间的稀疏表示系数的分布不是随机的，而是在某种程度上与图像信号的自相似性有关。西安电子科技大学在其专利申请“基于稀疏域噪声分布约束的SAR图像抑斑方法”(专利申请号：201410229349.2，公开号：CN103971346A)中提出了一种基于稀疏域噪声分布约束的SAR图像抑斑方法。该方法首先对SAR图像进行对数变换，利用K均值聚类的方法对图像块进行聚类。针对每类图像块，采用主成分分析的方法求出对应的子字典。然后，利用非局部均值算法和软阈值收缩算法得到最终的稀疏表示系数，通过反稀疏变换和反对数变换得到最终的SAR图像。该方法虽然解决了现有降斑方法不能同时保持图像纹理细节与辐射特性的问题，能对图像块进行聚类，但是不能保证所有图像块都得到了最好的划分。西安电子科技大学在其专利申请“基于SAR图像局部统计特性的KSVD相干斑抑制方法”(专利申请号：201110318457.3，公开号：CN102509263A)中提出了一种基于SAR图像局部统计特性的KSVD相干斑抑制方法。该方法从图像中随机抽取训练样本，利用SAR_KSVD算法对字典进行训练，得到最终的训练字典。然后，利用SAR_OMP算法进行稀疏编码的过程，得到稀疏编码系数。最后，根据冗余稀疏表示图像噪声抑制理论，得到最终的降斑图像。该方法能够很好的解决现有相干斑抑制方法中的边缘和纹理等细节信息被模糊的问题。但是由于该方法在提取训练图像块时并没有将图像块进行详细的划分，而是通过随机抽取的方法进行选取，导致训练得到的字典不能充分的包含图像中的纹理细节信息。
技术实现思路
本专利技术针对上述现有技术的不足，提出一种基于SSIM校正聚类稀疏表示的SAR图像降斑方法，以在保留点目标和纹理细节信息的前提下，提高SAR图像中的斑点噪声抑制效果。为实现上述目的，本专利技术的技术方案包括如下步骤：(1)输入任选的一幅合成孔径雷达SAR图像；(2)估计输入SAR图像中噪声的方差：(2a)对SAR图像进行方向波变换，得到低频子带系数矩阵和高频子带系数矩阵；(2b)按照下式，采用基于方向波域的斑点噪声估计方法，估计SAR图像中噪声的方差：其中，σ表示SAR图像中噪声的标准差，C表示噪声的归一化标准差，Ψ表示控制参数，μy表示SAR图像的均值，y表示SAR图像矩阵，表示SAR图像的高频子带系数矩阵的方差，Dy表示SAR图像的高频子带系数矩阵；(3)图像块分类：(3a)采用KSVD算法，对SAR图像进行预滤波，得到滤波后的SAR图像y′；(3b)以滤波后的SAR图像y′中任意像素为中心，抽取大小为8×8的图像块，得到SAR图像像素图像块集合B1，以SAR图像的低频子带系数矩阵中任意系数为中心，抽取大小为8×8的图像块，得到SAR图像低频子带图像块集合B2，合并SAR图像像素图像块集合B1和SAR图像低频子带图像块集合B2，得到总图像块集合B；(3c)计算总图像块集合B中每个图像块的方差，根据得到的图像块的方差大小，将总图像块集合B分为光滑图像块集合E和非光滑图像块集合F；(4)初始化字典：(4a)采用基于SSIM校对的K均值聚类算法，对非光滑图像块集合F中的所有图像块进行聚类，得到N类相似图像块集合Sk，1≤k≤N，N为K均值聚类所得到的类别总数；(4b)采用主成分分析方法，求解相似图像块集合Sk对应的64个特征向量，将这64个特征向量按列排列，得到Sk对应的字典Dk，1≤k≤N；(4c)将字典Dk作为相似图像块集合Sk的初始字典，将大小为64×256的离散余弦变换DCT字典作为光滑图像块集合E的初始字典；(5)根据基于结构聚类的稀疏表示方法和字典更新方法，对相似图像块集合Sk和光滑图像块集合E进行稀疏表示矩阵更新和字典更新，得到最终的稀疏表示矩阵和最终的字典(6)根据上述得到的总图像块集合B、最终的稀疏表示矩阵和最终的字典利用下式得到降斑后的SAR图像x：其中，λ表示拉格朗日因子，I表示与SAR图像大小一样的单位矩阵，R表示总图像块集合B中任意图像块矩阵，RT表示图像块矩阵R的转置。本专利技术与现有的技术相比具有以下优点：第一，由于本专利技术采用了方向波变换来获取不同尺度下的图像，克服了现有技术中不能充分结合不同尺度下图像块的自相似性而导致训练得到的字典包含信息量少的缺点，使得本专利技术充分的利用了图像的细节信息。第二，由于本专利技术采用了基于方向波域的斑点噪声估计方法，克服了现有技术中不能更精确的估计SAR图像中斑点噪声方差而导致降斑结果不理想的问题，使得本专利技术能够更好的去除图像中的斑点噪声。第三，由于本专利技术采用了基于SSIM校正的K均值聚类方法，克服了现有技术中利用K均值聚类方法对图像块进行聚类时而导致的图像块划分不准确的问题，使得本专利技术能够更好的保留图像中的细节信息。附图说明图1是本专利技术的实现流程图；图2为本专利技术与现有技术对haigang合成孔径雷达SAR图像的降斑效果对比图；图3为本专利技术与现有技术对field合成孔径雷达SAR图像的降斑效果对比图。具体实施方式下面结合附图对本专利技术做进一步详细描述。参照图1，本专利技术的具体步骤如下：步骤1：输入图像。输入任选的一幅合成孔径雷达SAR图像。本专利技术实施例中使用的合成孔径雷达SAR图像分别如附图2(a)和附图3(a)所示。其中，附图2(a)是合成孔径雷达SAR图像海港图haigang，大小为256×256，附图3(a)是合成孔径雷达SAR图像田野图field，大小为512×512。步骤2：估计输入SAR图像中的噪声方差。(2a)对SAR图像进行方向波变换，得到低频子带系数矩阵和高频子带系数矩阵；(2b)按照下式，采用基于方向波域的斑点噪声估计方法，估计SAR图像中噪声的方差：其中，σ表示SAR图像中噪声的标准差，C表示噪声的归一化标准差，Ψ表示控制参数，μy表示SAR图像的均值，y表示SAR图像矩阵，表示SAR图像的高频子带系数矩阵的方差，Dy表示SAR图像的高频子带系数矩阵。步骤3：图像块分类。(3a)采用KSVD算法，对SAR图像进行预滤波，得到滤波后的SAR图像y′，其中，KSVD算法的具体步骤如下：(3a1)以SAR图像y中任意像素为中心，抽取大小为8×8的图像块，得到的SAR图像像素图像块集合M，初始化字典D为离散余弦变换DCT字典，其大小为64×256；(3a2)按照下式，更新稀疏表本文档来自技高网...

【技术保护点】
一种基于SSIM校正聚类稀疏表示的SAR图像降斑方法，包括如下步骤：(1)输入任选的一幅合成孔径雷达SAR图像；(2)估计输入SAR图像中噪声的方差：(2a)对SAR图像进行方向波变换，得到低频子带系数矩阵和高频子带系数矩阵；(2b)按照下式，采用基于方向波域的斑点噪声估计方法，估计SAR图像中噪声的方差：σ2=C2(Ψμy2+σDy2)(1+C2),]]>其中，σ表示SAR图像中噪声的标准差，C表示噪声的归一化标准差，Ψ表示控制参数，μy表示SAR图像的均值，y表示SAR图像矩阵，表示SAR图像的高频子带系数矩阵的方差，Dy表示SAR图像的高频子带系数矩阵；(3)图像块分类：(3a)采用KSVD算法，对SAR图像进行预滤波，得到滤波后的SAR图像y′；(3b)以滤波后的SAR图像y′中任意像素为中心，抽取大小为8×8的图像块，得到滤波后的SAR图像y′像素图像块集合B1，以SAR图像的低频子带系数矩阵中任意系数为中心，抽取大小为8×8的图像块，得到SAR图像低频子带图像块集合B2，合并SAR图像像素图像块集合B1和SAR图像低频子带图像块集合B2，得到总图像块集合B；(3c)计算总图像块集合B中每个图像块的方差，根据得到的图像块的方差大小，将总图像块集合B分为光滑图像块集合E和非光滑图像块集合F；(4)初始化字典：(4a)采用基于SSIM校对的K均值聚类算法，对非光滑图像块集合F中的所有图像块进行聚类，得到N类相似图像块集合Sk，1≤k≤N，N为K均值聚类所得到的类别总数；(4b)采用主成分分析方法，求解相似图像块集合Sk对应的64个特征向量，将这64个特征向量按列排列，得到Sk对应的字典Dk，1≤k≤N；(4c)将字典Dk作为相似图像块集合Sk的初始字典，将大小为64×256的离散余弦变换DCT字典Dl作为光滑图像块集合E的初始字典；(5)根据基于结构聚类的稀疏表示方法和字典更新方法，对相似图像块集合Sk和光滑图像块集合E进行稀疏表示矩阵更新和字典更新，得到最终的稀疏表示矩阵和最终的字典(6)根据上述得到的总图像块集合B、最终的稀疏表示矩阵和最终的字典利用下式得到降斑后的SAR图像x：x=(λI+ΣRTR)-1(λy′+ΣRTD^a^),]]>其中，λ表示拉格朗日因子，I表示与SAR图像大小一样的单位矩阵，R表示总图像块集合B中任意图像块矩阵，RT表示图像块矩阵R的转置。...

【技术特征摘要】
1.一种基于SSIM校正聚类稀疏表示的SAR图像降斑方法，包括如下步骤：(1)输入任选的一幅合成孔径雷达SAR图像；(2)估计输入SAR图像中噪声的方差：(2a)对SAR图像进行方向波变换，得到低频子带系数矩阵和高频子带系数矩阵；(2b)按照下式，采用基于方向波域的斑点噪声估计方法，估计SAR图像中噪声的方差：其中，σ表示SAR图像中噪声的标准差，C表示噪声的归一化标准差，Ψ表示控制参数，μy表示SAR图像的均值，y表示SAR图像矩阵，表示SAR图像的高频子带系数矩阵的方差，Dy表示SAR图像的高频子带系数矩阵；(3)图像块分类：(3a)采用KSVD算法，对SAR图像进行预滤波，得到滤波后的SAR图像y′；(3b)以滤波后的SAR图像y′中任意像素为中心，抽取大小为8×8的图像块，得到滤波后的SAR图像y′像素图像块集合B1，以SAR图像的低频子带系数矩阵中任意系数为中心，抽取大小为8×8的图像块，得到SAR图像低频子带图像块集合B2，合并SAR图像像素图像块集合B1和SAR图像低频子带图像块集合B2，得到总图像块集合B；(3c)计算总图像块集合B中每个图像块的方差，根据得到的图像块的方差大小，将总图像块集合B分为光滑图像块集合E和非光滑图像块集合F；(4)初始化字典：(4a)采用基于SSIM校对的K均值聚类算法，对非光滑图像块集合F中的所有图像块进行聚类，得到N类相似图像块集合Sk，1≤k≤N，N为K均值聚类所得到的类别总数；所述的基于SSIM校对的K均值聚类算法，对非光滑图像块集合F中的所有图像块进行聚类，按如下步骤进行：(4a1)利用K均值聚类算法，对非光滑图像块集合F中的所有图像块进行N类聚类，得到每个图像块的类别和每个类别的聚类中心；(4a2)计算每个图像块与其所属类别的聚类中心的结构相似度SSIM值，判断该SSIM值是否小于0.85，若是，则执行步骤(4a3)，否则，执行步骤(4a4)；(4a3)计算该图像块与其他类别的聚类中心的SSIM值，将该图像块划分到其最大的SSIM值所对应的类别中；(4a4)将所属类别为k的所有图像块聚集起来，得到相似图像块集合Sk，1≤k≤N,N为K均值聚类所得到的类别总数；(4b)采用主成分分析方法，求解相似图像块集合Sk对应的64个特征向量，将这64个特征向量按列排列，得到Sk对应的字典Dk，1≤k≤N；(4c)将字典Dk作为相似图像块集合Sk的初始字典，将大小为64×256的离散余弦变换DCT字典Dl作为光滑图像块集合E的初始字典；(5)根据基于结构聚类的稀疏表示方法和字典更新方法，对相似图像块集合Sk和光滑图像块集合E进行稀疏表示矩阵更新和字典更新，得到最终的稀疏表示矩阵和最终的字典所述的稀疏表示矩阵更新和字典更新按如下步骤进行：(5a)按照下式，更新稀疏表示矩阵和聚类中心所对应的稀疏表示矩阵，得到更新后的稀疏表示矩阵和聚类中心所对应的稀疏表示矩阵

【专利技术属性】
技术研发人员：白静，焦李成，张钊，韩雪云，王爽，马晶晶，马文萍，侯彪，
申请(专利权)人：西安电子科技大学，
类型：发明
国别省市：陕西;61