《SAR影像斑点噪声的抑制方法》由会员分享，可在线阅读，更多相关《SAR影像斑点噪声的抑制方法（4页珍藏版）》请在装配图网上搜索。
1、二、SAR影像斑点噪声的抑制方法1、均值滤波均值滤波器是采用滤波窗口内所有像素灰度值的平均值来代替中心像素的值，均值滤波器具 有很好的噪声平滑能力，噪声标准差按窗口内像元数的均方根降低。然而，均值滤波器进行 平滑时对噪声和边缘信息不加区分，从而导致边缘信息临近区域分辨率下降，为了减少这一 问题，通常采用3x3或5x5的小滤波窗口。也正是由于它滤波时对噪声和信号不加以区 分的特点，对非噪声像素即信号也进行了平滑，从而不可避免地导致了影像的整体模糊和分 辨率下降。2、中值滤波中值滤波器是采用滤波窗口内所有像素的中值来代替中心像素的值，它能有效地孤立斑 点噪声。然而，这种滤波器存在边缘模糊，消除细的
2、线性特征以及目标形状扭曲等常见问题。 经中值滤波滤波后的影像失真度较大，纹理等细节信息损失较严重。3、局部滤波局域滤波基于方位把活动窗口分为8块(北、南、东、西、西北、东北、西南及东南) 对于每个区域，方差可由下式来计算：该算法比较窗口中心点附近8个局域方差值，窗口中心点的值就被区域中具有最小方差区 域像元的平均值所代替。方差平均值小的区域像元被认为受干涉的干扰小，这一点很像窗口 中心点的像元，方差小的区域相对于周围的区域来说受干涉影响小。4、Lee滤波在缺乏信号x的精确模型的情况下，使用影像本身从3x3或其它的滤波窗口内的局域均值 z和局域方差var(z)来估计信号的先验均值和方差。根据乘性
3、噪声模型，信号x的先验均 值和方差可以这样来估算：，假设线性滤波器的形式为，这里，。要注意的是必须确保var(x)为非负，如果为负则置 var(x)为0,否则可能在影像上引入认为的噪声成分。这一滤波方法的直观解释是，在均匀区域var(x)=0,滤波后的像素值(窗口内像素的平 均值)；对于高反差区域(或边缘)，var(x)较大，(像素本身的值)。该滤波器存在一个 问题是边缘区域的噪声并没有被平滑。5、Sigma 滤波Sigma滤波是基于高斯分布的Sigma概率，它通过对滤波窗口内落在中央像素的两个 Sigma范围内的像素进行平均来滤除影像噪声。高斯分布的两个Sigma概率是0.955， 即高斯分
4、布随机样本的95.5%都落在其均值的两个标准偏差范围内。对于乘性噪声模型而 言，两个Sigma范围是：.事先计算出所有灰度级(例如256个灰度级)的Sigma范围，并存储在数组中。对滤波 窗口内的中央像素，从数组中提取出Sigma范围值，将窗口内像素与这些上下限进行比较， 对落在上下限内的像素进行平均，并用平均值来替代中央像素的值。落在这两个Sigma范 围之外的像素将被忽略。如果没有其它窗口像素落在两个Sigma范围内时，引入一个阈值KS，如果落在Sigma范 围内的像素总数小于或等于KS时，就用中间像素的四个最近的相邻像素的平均值来替代。6、Frost 滤波Frost滤波器是特定大小窗口的
5、像素值和指数脉冲响应m卷积的Wiener自适应滤波器： 其中K是滤波器参数，代表中心像素的位置，|t|是距 的距离。这种响应是由目标反射率 的自回归指数模型得到的。Frost滤波器采用的斑点噪声模型采用的形式如下：。这里 是系统响应函数，*为卷积 算子。尽管该算法适用于任何系统响应函数，但在通常的应用中，一般假定为delta函数(例 如假定的功率谱密度在感兴趣的波段宽度上是不变的)。最小均方滤波器形式如下：。这里t对应于空间域中像素之间的距离。选择脉冲函数m(t)， 使下式最小：。按照频率域中Wiener滤波器的推导，可以容易地找到上式的解：，0衰减常数o的大小取决于x，var(x)和a。在应
6、用中，a取作一个常数，尽管它应当是与具 体图像有关的。其他两个量则通过5x5窗口内像素的局域均值和方差来估计。7、Gamma MAP滤波对于多视SAR图像，通过假设影像概率密度函数pdf为Gamma分布，应用最大概率(MAP) 滤除相干斑噪声，可得到Gamma MAP滤波器：其中，异质参数为。图像的局域方差系数。和 分别表示滤波窗口内图像像素的标准差和均值，相干斑的局域 方差系数，时视数。值得注意的是，下式只适用于多视SAR图像，处理单视图像时会存在估计偏差，必须对下 式作无偏修正。单视图像的无偏MAP估计为：8、基于小波分解的阈值滤波小波阈值主要有两种方法：一种是软阈值(Soft Thres
7、holding)方法，就是将小于某一 阈值的小波系数以0来代替，并将大于阈值的小波系数减去阈值作为新的小波系数的值； 另一种是硬阈值(Haard Thresholding)方法，就是直接将小于阈值的小波系数用0代替， 而大于阈值的值不做处理。软阈值去噪公式为：硬阈值去噪公式为：其中，软阈值的方法适用于小尺度，这是因为噪声集中在小尺度上，采用软阈值可以使其收 缩(变小)，减小其在重建的图像中所占的比例，这样就最大限度的减小了噪声的影响。反之， 因为在大尺度由于噪声成分相对较少，所以大于阈值的小波系数可以保留不变，这样图像的 特征就不会因为小波收缩而被削弱。1) 信号和噪声的交点作为阈值将近似部分
8、经多层分解所得的近似部分视为无噪的信号部分，将细节部分所得的细节部 分视为噪声成分。具体步骤为：将影像进行小波分解，做出近似部分经多层分解所得的近似部分和细节部分所 得的细节部分的统计曲线图，其交点所对应的横坐标视为阈值；然后进行软、硬阈值。下图为一信号噪声统计曲线图，其中绿色为噪声信号，红色为无噪信号。2) 基于小波系数的SOT结构进行滤波如上图，影像通过小波变换后，按其频带从低到高形成一个树状结构，树根是最低频带的结 点，它有3个子女，分别位于次低频子带的相应位置；其余子带（最高频子带除外）的结点 都有4个子女，位于高一级子带的相应位置，这样图1中的3级小波分解就形成深度为4 的树。基于小
9、波域中SOT结构的相干斑噪声抑制的算法步骤： 选取参数，将原始SAR影像进行3次分解，设置相同的阈值。 对分解的所有小波系数，基于SOT在处设置阈值。即：以的系数为根节点，对、和 设置阈值（一对一关系）；以、和的系数为根节点，分别对、和设置阈值（一对四关系）；以、和的系数为根节点，分别对、和设置阈值（一对四关系）； 进行阈值处理，即对9块系数分别进行阈值处理。 进行逆小波变换。二、SAR影像滤波效果质量评价1、平均值图像像素值的均值定义为：滤波后影像的均值越接近原始影像，说明影像的辐射特性保持的越好。2、标准差标准差定义为：滤波后影像的标准差越接近原始影像，则影像的灰度动态范围保持的越好。3、
10、辐射分辨率辐射分辨率的定义为：辐射分辨率是SAR反映地物目标微波散射特性精度的衡量。4、有效视数（ENL）有效视数又称等效视数，其计算公式为：式中，、分别为图像灰度4次方的均值、平方的均值。有效视数反映去斑点噪声的能力，ENL越大，去斑点噪声的能力越强。5、信噪比（SNR）信噪比的定义为：信噪比值越高，则影像的去噪效果越好。6、峰值信噪比（PSNR）峰值信噪比的定义为：峰值信噪比值越高，则影像的去噪效果越好。7、边缘保持（ESI）边缘保持的定义为滤波后影像的边缘与滤波前之比。边缘保持系数越接近1，则边缘保 持的效果越好。经典算法结果比较：七种常用经典滤波算法（均值滤波、中值滤波、局部滤波、Le
11、e滤波、Lee-Sigma滤波、 Frost滤波和Gamma-MAP滤波）中，Lee滤波和Gamma-MAP滤波的效果最好。均值滤波和中值滤波保持边缘的能力较差；局部滤波虽然有较好的保持边缘能力，但是信噪 比和峰值信噪比很低；Lee-Sigma滤波和Frost滤波的有效视数较低，而且对原始影像均 值的改变较大。Lee滤波和Gamma-MAP滤波在各方面都有较好的效果。非经典算法结果比较：Wiener滤波、一级db32小波分解的Wiener滤波、一级haar小波分解的Wiener滤波、 一级db32小波分解的对数变换Wiener滤波和一级db32小波分解的软阈值滤波算法中， 使用haar小波基的
12、效果最差，目视效果就出现很大的模糊，信噪比的数值也相当低；一级db32小波分解的Wiener滤波算法对于原始影像是否进行对数变换的效果相当，在统 计量上的表现也是如此；直接Wiener滤波信噪比和峰值信噪比数值较高，边缘保持能力也比较强，但是有效视数偏 低；这五个算法中，一级db32小波分解的软阈值滤波效果最好，整幅影像的均值保持很好，信 噪比、峰值信噪比和有效视数都比较高，边缘的保持能力也比较强。基于交点的软硬阈值滤波总体说来，效果不如经典滤波算法，它降低了影像的信噪比和有效 视数。经典算法和非经典算法的比较：一级db32小波分解的软阈值滤波算法与经典Lee滤波、Gamma-MAP滤波相比，
13、影像均 值保持、边缘保持效果比较好；信噪比、峰值信噪比优于Gamma-MAP滤波，但是不及 Lee滤波；辐射分辨率和有效视数有待于进一步的提高。对于小波分解方法的总结：基于小波分解的方法是在不同分解级别上对低频图像进行低通滤波去噪，并加上相应级别的 高频图像的边缘信息，实现在充分抑制图像斑点噪声的同时，保留不同尺度的结构和纹理特 征。但是基于小波分解的方法所得的结果并不理想。最主要的原因：没有将噪声和信号进行有效 的分离，使得滤得的影像既包含一些噪声信号，又缺失了一部分无噪信号。应该根据傅立叶 频谱设计有效的方法进行频段分离，继而在相应的某层小波系数上设计滤波器滤波，这也是 一个难点。以上所讨论的各种方法，都在不同程度上抑制了斑点噪声，保留了边缘信息。但是效果都不 是很理想，故需要引入影像融合，以不仅滤除噪声等非有效信息的干扰，而且较大程度上保 持了图像的细节信息
SAR影像斑点噪声的抑制方法
