clear;clc;
figure;
yt=imread（‘H:\大學\數字圖象處理\課程設計2\dogOriginal.bmp‘）;
yt_noise=imread（‘H:\大學\數字圖象處理\課程設計2\dogDistorted.bmp‘）;
subplot（431）;imshow（yt）;title（‘原圖‘）;
subplot（433）;imhist（yt）;title（‘原圖‘）;
a=fftshift（fft2（yt））;
subplot（432）;imshow（mat2gray（log（abs（a））））;title（‘原圖頻譜圖‘）;
subplot（434）;imshow（yt_noise）;title（‘噪聲圖‘）;
subplot（436）;imhist（yt_noise）;title（‘噪聲圖‘）;
a=fftshift（fft2（yt_noise））;
subplot（435）;imshow（mat2gray（log（abs（a））））;title（‘噪聲圖頻譜圖‘）;
A=fspecial（‘average‘）;
qsjnoise=filter2（Ayt_noise）/255;?
subplot（437）;imshow（qsjnoise）;title（‘去隨機噪聲‘）;
subplot（439）;imhist（qsjnoise）;title（‘去隨機噪聲‘）;
a=fftshift（fft2（qsjnoise））;
subplot（438）;imshow（mat2gray（log（abs（a））））;title（‘去隨機噪聲后頻譜圖‘）;
x=double（qsjnoise）;
x1=fftshift（fft2（qsjnoise））;
x2=x1;
[mn]=size（x1）;
nn=2;
d0=30;
m0=fix（m/2）;
n0=fix（n/2）;
for?i=1:m
????for?j=1:n
????????d=sqrt（（i-m0）^2+（j-n0）^2）;
????????h=1*exp（-1/2*（d^2/d0^2））;?
????????x2（ij）=h*x1（ij）;
????end
end
qnoise=abs（ifft2（x2））;
qnoise=mat2gray（qnoise）;
subplot（4310）;imshow（qnoise）;title（‘去噪聲‘）;
subplot（4312）;imhist（qnoise）;title（‘去噪聲‘）;
a=fftshift（fft2（qnoise））;
subplot（4311）;imshow（mat2gray（log（abs（a））））;title（‘去噪聲后頻譜圖‘）;

cha=yt-uint8（qnoise）;
he=sum（sum（cha））;
avg=he/（m*n）;
wucha=cha-avg;
result1=（sum（sum（wucha.^2）））/（m*n）;
if?result1==0
disp（‘dog圖均方誤差‘）;result2=0
else
disp（‘dog圖均方誤差‘）;result2=sqrt（result1）
end




figure;
yt=imread（‘H:\大學\數字圖象處理\課程設計2\swanOriginal.bmp‘）;
yt_noise=imread（‘H:\大學\數字圖象處理\課程設計2\swanNoise.bmp‘）;
subplot（431）;imshow（yt）;title（‘原圖‘）;
subplot（433）;imhist（yt）;title（‘原圖‘）;
a=fftshift（fft2（yt））;
subplot（432）;imshow（mat2gray（log（abs（a））））;title（‘原圖頻譜圖‘）;
subplot（434）;imshow（yt_noise）;title（‘噪聲圖‘）;
subplot（436）;imhist（yt_noise）;title（‘噪聲圖‘）;
a=fftshift（fft2（yt_noise））;
subplot（435）;imshow（mat2gray（log（abs（a））））;title（‘噪聲圖頻譜圖‘）;
qjynoise=medfilt2（yt_noise[33]）;
subplot（437）;imshow（qjynoise）;title（‘去隨機噪聲‘）;
subplot（439）;imhist（qjynoise）;title（‘去隨機噪聲‘）;
a=fftshift（fft2（qjynoise））;
subplot（438）;imshow（mat2gray（log（abs（a））））;title（‘去隨機噪聲后頻譜圖‘）;
x=double（qjynoise）;
x1=fftshift（fft2（qjynoise））;
x2=x1;
[mn]=size（x1）;
nn=2;
d0=30;
m0=fix（m/2）;
n0=fix（n/2）;
for?i=1:m
????for?j=1:n
????????d=sqrt（（i-m0）^2+（j-n0）^2）;
????????h=1/（1+0.414*（d/d0）^（2*nn））;
????????x2（ij）=h*x1（ij）;
????end
end
qnoise=abs（ifft2（x2））;

subplot（4310）;imshow（qnoise[]）;title（‘去噪聲‘）;
subplot（4312）;imhist（mat2gray（qnoise））;title（‘去噪聲‘）;
a=fftshift（fft2（qnoise））;
subplot（4311）;imshow（mat2gray（log（abs（a））））;title（‘去噪聲后頻譜圖‘）;

cha=yt-uint8（qnoise）;
he=sum（sum（cha））;
avg=he/（m*n）;
wucha=cha-avg;
result1=（sum（sum（wucha.^2）））/（m*n）;
if?result1==0
disp（‘swan圖均方誤差‘）;result2=0
else
disp（‘swan圖均方誤差‘）;result2=sqrt（result1）
end







%%%%帶阻濾波器，線性濾波器

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----

?????文件?????867328??2018-05-18?07:58??Matlab圖像增強\圖像處理課程設計2.doc

?????文件?????132442??2018-05-18?07:59??Matlab圖像增強\dogDistorted.bmp

?????文件?????132442??2018-05-18?07:59??Matlab圖像增強\dogOriginal.bmp

?????文件???????3041??2018-06-20?20:24??Matlab圖像增強\keshe2.m

?????文件?????188578??2018-05-18?07:59??Matlab圖像增強\swanNoise.bmp

?????文件?????188578??2018-05-18?07:59??Matlab圖像增強\swanOriginal.bmp

?????文件????1480409??2018-06-15?09:36??Matlab圖像增強\圖像處理課程設計.docx

?????目錄??????????0??2018-10-14?23:22??Matlab圖像增強

-----------?---------??----------?-----??----

??????????????2992818????????????????????8