Perona_Malik 图像的各向异性扩散

图像的各向异性扩散，利用matlab实现Perona_Malik算法，梯度阈值K通过对图像求解2-范数的方法迭代获得

close all
clear
clc

%原始图像的读取与显示
im=imread（‘lenna.bmp‘）;
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;
imshow（imB）;
title（‘含噪图像‘）;


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

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

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

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件      66616  2010-07-04 10:47  Perona_Malik\lenna.bmp

     文件        816  2010-07-07 18:39  Perona_Malik\main.asv

     文件        816  2010-07-07 18:39  Perona_Malik\main.m

     文件       1612  2010-07-07 15:25  Perona_Malik\PM.m

     文件        201  2010-07-07 18:34  Perona_Malik\SNR.m

     目录          0  2010-07-09 16:15  Perona_Malik

----------- ---------  ---------- -----  ----

                70061                    6