一种基于非本地分类稀疏表示的SAR图像降噪方法技术

本发明专利技术公开了一种基于非本地分类稀疏表示的SAR图像降噪方法。属于数字图像处理技术领域。它是一种基于分类系数估计与稀疏表示相结合的SAR图像降噪方法。首先按照图像的同质区和异质区进行分类，而后分别采用非本地均值方法和3D变换域硬阈值收缩方法实现对同质区和异质区的降噪，并由降噪结果获得非本地系数。为进一步提高降噪性能，利用所估计出的非本地系数对稀疏表示模型中的系数进行约束，以更好的逼近真实图像的系数。最后,利用最终估计的系数重构出降噪后的图像。本发明专利技术对图像的同质区和异质区采用不同方法进行处理，不仅可有效的实现同质区抑斑，同时能保留异质区中存在的细节，可用于SAR图像降噪。

【技术实现步骤摘要】
一种基于非本地分类稀疏表示的SAR图像降噪方法
本专利技术属于数字图像处理
，它特别涉及基于图像非本地相似性与稀疏表示相结合的降噪方法，用于SAR图像降噪处理。
技术介绍
合成孔径雷达SAR被广泛的应用于多个领域，并在国民经济发展与国防安全保卫等领域中发挥着重要的作用。但SAR图像的相干成像机理使其成像后的图像中存在严重的相干斑噪声，对SAR图像的后处理，如图像分割、边缘检测和特征提取等带来了一定影响，不利于SAR图像的解译。因此，对相干斑噪声的有效抑制便成为SAR图像处理的一个重要环节。信号的稀疏表示能有效提取信号的本质特征，基于稀疏表示的图像降噪因其良好的降噪效果而得到国内外学者广泛关注。其中一些方法通过将稀疏表示与其他降噪方法相结合来提高降噪的性能，如将稀疏降噪与均值滤波，全微分相结合等。尤其是将稀疏降噪与非本地冗余特性相结合的降噪方法取得了较好的结果，如BM3D及其改进方法在降噪的整体效果方面达到了领先水平。
技术实现思路
本专利技术的目在于针对现有SAR图像降噪中存在的不足，提出一种基于非本地相似性与稀疏表示相结合的SAR降噪方法。该方法充分考虑了SAR图像中不同区域内的特征与噪声分布的情况，首先对图像分类，然后不同的方法对不同的图像块进行降噪以估计出每一图像块对应的非本地稀疏系数，最后利用非本地稀疏系数调节本地稀疏系数使最终的稀疏解更接近于真实图像的稀疏系数。因此，该方法可以提高稀疏降噪的性能。包括以下步骤：步骤一、图像内容的分类对于图像中的每个像素点，计算以该点为中心的图像块的方差，然后将其与预先设置好的阈值相比较，若大于阈值，则该中心像素点被标记为异质区，若小于阈值，则标记为同质区。步骤二、分类估计非本地稀疏系数首先，根据目标图像块分类后的属性设定相应的图像块大小与搜索区域。然后，由式(2)计算出目标图像块与搜索范围内候选图像块的相似度：其中yi和yj分别表示目标图像块和候选图像块，ω(yi,yj)为yi和yj之间的相似度，k为图像块内像素，L为SAR图像视数，ζ为权重调节参数，表示上一次降噪结果，初始降噪时为0。在求得候选图像块与目标图像块之间的相似度后，从中选择相似度最大的K个图像块作为目标图像块的相似图像块。对于异质区目标图像块，将其本身及相似图像块联合构造为3D矩阵，并对其进行小波变换与硬阈值操作，再经过小波逆变换得到目标图像块的降噪结果；对于同质区的目标图像块，将其对应的相似图像块进行加权平均以实现目标图像块的降噪处理。为获得稀疏表示中更适合该图像的字典，分别对同质区和异质区的目标图像块进行Kmeans聚类，并对每一类图像块使用主成分分析法训练相应的自适应学习字典D，最后将降噪后的目标图像块在该字典上投影，获得同质区目标图像块的非本地稀疏系数ηi和异质区非本地稀疏系数γi。步骤三、结合非本地稀疏系数与稀疏表示的SAR图像降噪为了提高图像降噪性能，利用非本地稀疏系数调节稀疏解使其更接近真实解，即将非本地稀疏系数代入以下稀疏模型进行求解：其中y表示含噪图像，D为字典，N为图像所包含的像素点个数，ki为分类后像素所属类型的标记，αi为本地稀疏系数，γi和ηi为非本地稀疏系数，λ1，λ2和λ3为权重调节参数；为了求出最终的降噪图像，首先将目标图像块在字典D上投影得到本地稀疏系数αi，即：αi＝DTyi式(4)然后利用式(5)和式(6)求解式(3)中的权重参数λ1和λ2或λ1和λ3：其中σi表示αi的标准差，δi和χi分别表示αi-ηi和αi-γi的标准差，||y||2表示图像块的2范数。在D、αi、λ1、λ2、λ3、ηi、γi确定后，利用迭代收缩算法求解式(3),得到最终的稀疏解最后，使用式(7)将重构得到降噪后的图像：将降噪后图像与初始含噪图像进行凸组合后的结果作为新的含噪图像，然后重复步骤二和步骤三的降噪处理，当相邻两次降噪所获图像间的平均绝对偏差MAD小于设置的阈值时结束循环。本专利技术的创新点是利用图像的非本地冗余信息提高图像稀疏降噪的性能；利用不同的方法实现对不同类型图像块的降噪以获得非本地稀疏系数；并利用估计出的非本地稀疏系数对稀疏表示模型中的稀疏系数进行约束，以进一步的提高降噪性能，并将该方法用于SAR图像降噪。本专利技术的有益效果：稀疏降噪的过程中结合图像的非本地冗余特性，一方面可以训练出更适合图像表示的字典，另一方面可调节本地稀疏系数使其更接近于干净图像的稀疏系数，从本质上提高稀疏降噪的性能；在求解非本地稀疏系数时，根据图像块内容的不同而选择不同的参数和方法以提高非本地稀疏系数的估计精度，从而更好的调节本地稀疏系数；从稀疏表示降噪原理出发，获得乘性噪声下的稀疏降噪模型及参数，以更好的保持SAR图像的固有特征，有利于SAR图像的后续处理。本专利技术主要采用仿真实验的方法进行验证，所有步骤、结论都在MATLAB8.0上验证正确。附图说明图1是本专利技术的工作流程框图；图2是本专利技术仿真中使用的真实含噪的田园SAR图像；其中白色矩形区域为选择的同质区；图3是PPB方法对图2的降噪结果图；图4是BM3D方法对图2的降噪结果图；图5是本专利技术方法对图2的降噪结果图。具体实施方式参照图1，本专利技术基于非本地分类稀疏表示的SAR图像降噪方法，具体步骤包括如下：步骤一、图像内容的分类首先，对SAR图像进行简单的高斯空域滤波，方差为1.2；然后计算8×8大小图像块的方差c＝σ/μ。设置阈值其中σn为相干斑的标准差，n为图像块内像素点个数。若c大于T，则该中心像素点属于异质区，并设ki＝1；反之，属于同质区，并设ki＝0。步骤二、分类估计非本地稀疏系数对属于同质区的目标图像块，设置图像块大小为9×9，搜索范围的大小为15×15，当目标图像块属于异质区时设置图像块的大小为3×3，搜索范围的大小为39×39；首先通过式(2)计算出搜索范围内K个候选图像块的之间相似度，其中ζ＝1，K＝16，并将这些候选图像块设置为初始相似图像块的集合，对相似度大小进行排序，然后逐个计算搜索范围内剩余候选图像块的相似度，若小于当前相似图像块集合中候选图像块的相似度则丢弃，否则将该图像块插入相似图像块集合并更新相似图像块集合。对异质区目标图像块及其相似图像块构成的3D矩阵采用3D变换域阈值收缩方法来进行非本地估计，如式(8)所示：其中Yi表示相似图像块组成的3D矩阵，Γ3D和分别表示小波变换和小波逆变换，τ表示阈值处理，且阈值的大小由局部噪声方差决定。当目标图像块属于同质区时，通过加权平均对其进行估计，如式(9)所示：其中表示图像块yi的非本地估计。利用Kmeans聚类算法，分别将同质区和异质区的目标图像块分为15类，针对每类图像块，利用主成分分析法训练自适应字典D，将降噪后的图像块在字典D上投影得到非本地稀疏系数ηi或γi。步骤三、结合估计稀疏系数与稀疏表示的SAR图像降噪首先利用式(4)求解目标图像块在字典D上投影的本地稀疏系数αi，采用式(5)和式(6)求解式(3)中的权重参数λ1，λ2和λ3，并结合步骤二中估计出的非本地稀疏系数ηi或γi，并由式(10)计算出式(3)的解其中λi表示ki确定后，式(3)右端非本地约束项的权重参数，βi为对应的非本地稀疏系数，的表达式如下：在得到稀疏解后，利用式(7)重构稀疏解得到降噪本文档来自技高网...

【技术保护点】
一种基于非本地分类稀疏表示的SAR图像降噪方法，其特征在于具体步骤如下：步骤一、图像内容的分类SAR图像中不同区域内噪声分布不同，并且不同区域间图像本身的特性也有所不同，为实现对不同图像内容采用不同方法降噪，将图像块的方差与预设的阈值进行比较，当方差大于阈值时视为异质区，小于阈值时视为同质区；步骤二、分类估计非本地系数结合分类后SAR图像中同质区和异质区的特点，选择不同的方法对其进行降噪处理；首先以图像的每个像素点为中心，依次抽取图像块，若图像块属于异质区，则根据异质区图像块包含奇异结构和细节边缘而难以在像素点附近范围内找到真正相似图像块的特点，选择较小的块尺寸并在较大的邻域范围内搜索与该目标块相似的图像块，并在得到相似图像块后使用3D变换域硬阈值收缩方法实现对目标图像块的降噪；对于同质区目标块，由于同质区内的像素点所包含的真实信号服从相同的分布，因此设置较大的块尺寸和较小的搜索范围以寻找相似块，并利用加权平均的方法实现对同质区目标块的降噪；然后利用主成分分析法训练出自适应学习字典，最后将降噪后的图像块在对应的自适应学习字典上投影得到图像块的非本地稀疏系数；步骤三、结合非本地系数与稀疏表示的SAR图像降噪为进一步提高图像降噪的性能，利用步骤二得到的非本地系数调节稀疏解使其更接近于真实解，即将非本地系数代入式(1)求解：α^=arg min||y-Dα||22+λ1ΣN||αi||1+λ2ΣNki||αi-γi||1+λ3ΣN(1-ki)||αi-ηi||1]]>  式(1)其中y表示含噪图像，D为字典，N为图像所包含的像素点个数，ki为分类后像素所属类型的标记，αi为本地稀疏系数，γi和ηi为非本地稀疏系数，λ1，λ2和λ3为权重调节参数；式(1)右端前两项为图像的本地稀疏表示，后两项为使本地稀疏系数αi更接近于真实解的约束项；由于图像内容仅分为同质区和异质区，ki仅为1或0，因此针对每个图像块，求解式(1)时均会从最后两项中选择出一项对进行αi约束，此时的优化问题转化为二范数求解问题；在确定了参数λ1，λ2和λ3后，先将目标图像块在字典上投影得到本地稀疏系数αi，然后利用非本地系数γi或ηi调节αi，得到最终稀疏解并将重构得到降噪后图像；将重构后的降噪图像与原含噪图像非凸组合后继续执行步骤二和步骤三，当前后两次降噪结果的平均绝对偏差MAD小于设定的阈值时停止循环，得到最终降噪图像。...

【技术特征摘要】
1.一种基于非本地分类稀疏表示的SAR图像降噪方法，其特征在于具体步骤如下：步骤一、图像内容的分类SAR图像中不同区域内噪声分布不同，并且不同区域间图像本身的特性也有所不同，为实现对不同图像内容采用不同方法降噪，将每个目标图像块的方差与预设的阈值进行比较，当方差大于阈值时目标图像块中心像素点归属为异质区，当方差小于阈值时目标图像块中心像素点归属为同质区；步骤二、分类估计非本地稀疏系数结合分类后SAR图像中同质区和异质区的特点，选择不同的方法对其进行降噪处理；首先以图像的每个像素点为中心，依次抽取目标图像块，若目标图像块属于异质区，则根据异质区图像块包含奇异结构和细节边缘而难以在像素点附近范围内找到真正相似图像块的特点，选择较小的图像块尺寸并在较大的邻域范围内搜索与该目标图像块相似的图像块，并在得到相似图像块后使用3D变换域硬阈值收缩方法实现对目标图像块的降噪；对于同质区目标图像块，由于同质区内的像素点所包含的真实信号服从相同的分布，因此设置较大的图像块尺寸和较小的搜索范围以寻找相似图像块，并利用加权平均的方法实现对同质区目标图像块的降噪；然后利用主成分分析法训练出自适应学习字典，最后将降噪后的目标...

【专利技术属性】
技术研发人员：刘书君，吴国庆，张新征，杨婷，徐礼培，
申请(专利权)人：重庆大学，
类型：发明
国别省市：重庆;85