基于相似值的SAR图像相干斑非局部均值抑制方法技术

本发明专利技术提出了一种基于相似值的SAR图像相干斑非局部均值抑制方法，用于解决现有技术中存在的相干斑噪声抑制能力较弱的技术问题，实现步骤为：对SAR图像I

【技术实现步骤摘要】
基于相似值的SAR图像相干斑非局部均值抑制方法
本专利技术属于雷达图像处理
，涉及一种SAR图像相干斑抑制方法，具体涉及一种基于相似值的SAR图像相干斑非局部均值抑制方法，可用于对SAR图像进行分类、目标检测和识别。
技术介绍
合成孔径雷达(SyntheticApertureRadar，SAR)具有全天时，全天候工作，穿透能力强的优点，被广泛应用于目标检测识别等领域。但是由于SAR相干成像技术往往使获得的SAR图像具有大量的相干斑噪声，严重影响了后续的SAR图像处理与解译。因此，相干斑抑制是SAR图像预处理过程中的重要步骤，具有重要的研究意义与实用价值。评价SAR图像相干斑抑制能力的指标主要有三个，第一个是等效视数(EquivalentNumbersofLooks,ENL)，用来评价SAR图像匀质区域的光滑程度的一个指标，ENL值越高，表示图像上的匀质区域光滑度越好；第二个是边缘保持指数(EdgePreservationIndex,EPI)，用来评价滤波处理后对图像的水平/垂直方向的边缘保持能力，EPI值越接近1，表示边缘保持能力本文档来自技高网...

【技术保护点】
1.一种基于相似值的SAR图像相干斑非局部均值抑制方法，其特征在于，包括以下步骤：/n(1)对SAR图像I

【技术特征摘要】
1.一种基于相似值的SAR图像相干斑非局部均值抑制方法，其特征在于，包括以下步骤：
(1)对SAR图像I1进行预处理：
计算大小为H×W的SAR图像I1中所有像素点像素值的最大值并通过对图像I1进行归一化，得到取值在[0,1]的归一化SAR图像I2；
(2)计算SAR图像I2的等效视数L和噪声标准差
计算SAR图像I2中所有像素点像素值的均值μ以及所有像素点像素值的方差σ2，并根据μ和σ2计算SAR图像I2的等效视数L，同时对SAR图像I2进行小波变换，并根据小波变换得到的对角高频子带小波系数HH计算SAR图像I2的噪声标准差
(3)计算SAR图像I2的邻域窗边长w1和搜索窗边长w2：
(3a)定义对等效视数L变换后范围在[3,10]的函数为f(L)，对L变换后范围在[1,6]的函数为g(L)：






其中，表示向上取整；
(3b)根据f(L)和g(L)确定SAR图像I2的邻域窗边长w1和搜索窗边长w2：
w1＝f(L),w1∈[3,10]
w2＝w1+g(L),w2∈[w1+1,w1+6]；
(4)对SAR图像I2进行区域划分：
通过以w1为边长的邻域窗，并按照先行后列的顺序对SAR图像I2进行逐像素滑动，得到图像块集合U1＝{up|1≤p≤(H-w1+1)×(W-w1+1)}，其中，up表示第p个图像块；
(5)获取每个图像块up的邻域：
选取以每个图像块up的中心为中心，以边长为w2的搜索窗内除up以外的所有图像块作为up的邻域集合U2＝{uq|1≤q≤(w2-w1+1)×(w2-w1+1)}，其中，uq表示up的邻域中的第q个图像块；
(6)基于相似值计算每个图像块up与其每个邻域图像块uq之间的滤波权值ω(up,uq)：
(6a)计算每个图像块up与其每个邻域图像块uq之间的像素强度差异值d1(up,uq)；
(6b)计算每个图像块up与其每个邻域图像块uq之间的结构差异值d2(up,uq)：
(6b1)按照如下公式对SAR图像I2进行滤波，得到差分SAR图像I3：



其中，I2(i,j)表示SAR图像I2第i行第j列像素的像素值，I3(ρ1,ρ2)表示SAR图像I3第ρ1行第ρ2列像素的像素值；
(6b2)通过I3中所有像素点像素值的最大值对I3进行归一化，得到取值在[0,1]的归一化差分图像I4，并通过I4构造有关像素差值的积分图Ot：



其中,ot(z)表示图像I4有关像...

【专利技术属性】
技术研发人员：刘峥，王梦白，谢荣，冉磊，
申请(专利权)人：西安电子科技大学，
类型：发明
国别省市：陕西;61