clear?all;
tic;
a?=imread（‘cameraman.tif‘）;
[mn]=size（a）;
a=imnoise（a‘salt?&?pepper‘0.02）;???????????????%加椒鹽噪聲
b=im2double（a）?;?????????????????????????????????%由于matlab中讀入的圖像是uint8的，在求平均值時(shí)很容易溢出，因此將其轉(zhuǎn)換成double型進(jìn)行運(yùn)算

for?i=2:m-1
????for?j=2:n-1
??????I1=b（i-1j-1）-b（ij）;??????????????????????%求鄰域和中心像素的差值
??????I2=b（i-1j）-b（ij）;
??????I3=b（i-1j+1）-b（ij）;
??????I4=b（ij-1）-b（ij）;
??????I6=b（ij+1）-b（ij）;
??????I7=b（i+1j-1）-b（ij）;
??????I8=b（i+1j）-b（ij）;
??????I9=b（i+1j+1）-b（ij）;
?????
??????q=（I1+I2+I3+I4+I7+I6+I8+I9+9*b（ij））/9;????%將中心像素周圍的灰度值相加求平均值
??????b（ij）=q;
????end
end
c=im2uint8（b）;??????????????????????????????????%將double型數(shù)據(jù)轉(zhuǎn)換成uint8
figuresubplot（121）imshow（a）;
title（‘原圖‘）;
imwrite（a‘原圖.jpg‘）;
subplot（122）imshow（c）;
title（‘平均濾波‘）;
imwrite（c‘平均濾波.jpg‘）;
toc;

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

?????文件????????918??2011-12-04?10:09??第一題\average_filter.m

?????文件???????1718??2011-12-04?19:01??第一題\direction.m

?????文件???????1870??2011-12-04?10:08??第一題\k_neighborhood.m

?????文件???????9727??2011-12-04?19:08??第一題\K鄰域圖像.jpg

?????文件??????98816??2011-12-04?19:08??第一題\三種方法對(duì)比+編程解釋.doc

?????文件??????19798??2011-12-04?19:08??第一題\原圖.jpg

?????文件??????10087??2011-12-04?19:08??第一題\平均濾波.jpg

?????文件??????14347??2011-12-04?19:08??第一題\方向平滑.jpg

?????目錄??????????0??2011-12-05?17:01??第一題

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

???????????????157281????????????????????9