非线性扩散--PM算法完整matlab代码

很容易理解的PM模型代码，输入原始图像，多次迭代PM算法得到平滑后的图像。非线性扩散方法在图像处理方面的应用越来越广泛，因其对于图像边缘的增强和保护作用十分明显，所以对该种方法的深入研究是十分有必要的。

close all
clear
clc

%原始图像的读取与显示
%im=imread（‘lenna.bmp‘）;
im=imread（‘VV.jpg‘）;

imshow（im）;
title（‘原始图像‘）;

% %高斯低通滤波得到模糊图像
% h=fspecial（‘gaussian‘[33]1）;%高斯低通滤波器（采用3*3的模板，标准差为1（默认的为3*3模板，标准差为0.5））
% imA=imfilter（imh）;
% figure;
% imshow（uint8（imA））;
% title（‘模糊图像‘）;

%添加高斯白噪声
imB=imnoise（im‘gaussian‘00.003）;
SNR（imimB）
%imB=imA+randn（size（imA））*5;
figure（1）;
imshow（imB）;
title（‘含噪图像‘）;

PM_image=PM（imB2000.022）;
SNR（imPM_image）
figure（2）;
imshow（uint8（PM_image））;
title（‘50次迭代后的效果图‘）;

PM_image=PM（imB3000.022）;
SNR（imPM_image）
figure（3）;
imshow（uint8（PM_image））;
title（‘100次迭代后的效果图‘）;

PM_image=PM（imB5000.022）;
SNR（imPM_image）
figure（4）;
imshow（uint8（PM_image））;
title（‘200次迭代后的效果图‘）;

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     目录           0  2014-05-15 21:21  Perona_Malik\
     文件        1699  2014-05-14 22:09  Perona_Malik\PM.m
     文件         201  2010-07-07 18:34  Perona_Malik\SNR.m
     文件      304774  2014-05-15 21:21  Perona_Malik\VV.jpg
     文件       66614  2014-05-12 20:22  Perona_Malik\aniso-synimgn21.bmp
     文件       66614  2014-05-12 20:22  Perona_Malik\aniso-synimgn22.bmp
     文件       66614  2014-05-12 20:22  Perona_Malik\aniso-synimgn23.bmp
     文件       66614  2014-05-12 20:22  Perona_Malik\aniso-synimgn24.bmp
     文件       66616  2010-07-04 10:47  Perona_Malik\lenna.bmp
     文件         816  2010-07-07 18:39  Perona_Malik\main.asv
     文件         856  2014-05-15 21:24  Perona_Malik\main.m