【图像去噪】非线性扩散PM算法图像去噪【含Matlab源码 2130期】

⛄一、PM模型图像降噪简介

1 引言
数字图像在获取、存储和传输中总会受到噪声的影响，因此图像去噪一直是备受关注的研究问题。作为图像处理技术的一大分支，基于偏微分方程的方法在图像去噪领域发挥着功不可没的作用。其中，PM (Perona-Malik, PM）模型是一种典型的非线性扩散模型，

其中，k表示噪声阈值，用来控制沿着梯度方向的扩散量。

PM模型的显式离散解为：

PM模型利用图像的梯度模对扩散系数进行调节，在平坦区域进行大的扩散，在边沿区域主要沿着边沿的切线方向进行扩散，具有很好保持边沿的作用。然而，PM去噪后图像却存在明显的阶梯效应，且容易引入斑点噪声。

2 本文算法
在PM模型中，扩散系数直接影响去噪效果。经典PM模型中，扩散系数直接依赖于待求像素点的梯度模，而受到噪声的影响，尤其随着噪声水平的提高，梯度的计算将产生很大的误差。文献[6]提出利用噪声和图像边沿在邻域范围内的不同结构特征，使用待求像素点的邻域像素值来构造水平和垂直两个方向的扩散函数，达到更好的保护边沿的作用，然而该方法仍然会在去噪图像中引入斑点噪声。受到文献[6]的启发，本文从水平、垂直，45°对角和135°对角4个方向来构造扩散函数：