close?all
clear
clc

%原始圖像的讀取與顯示
im=imread（‘lenna.bmp‘）;
imshow（im）;
title（‘原始圖像‘）;

%?%高斯低通濾波得到模糊圖像
%?h=fspecial（‘gaussian‘[33]1）;%高斯低通濾波器（采用3*3的模板，標(biāo)準(zhǔn)差為1（默認(rèn)的為3*3模板，標(biāo)準(zhǔn)差為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次迭代后的效果圖‘）;

?屬性????????????大小?????日期????時(shí)間???名稱(chēng)
-----------?---------??----------?-----??----

?????文件??????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