clear?all;
x=imread（‘020.bmp‘）;%輸入文件
g=double（rgb2gray（x））;
[growgcol]=size（g）;
Dij（11:8）=0;
?
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%旋濾波部分的程序，分成8個方向，取方差最小為濾波方向。
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
for?i=4:grow-4
????for?j=4:gcol-4

???????B（1:81:7）=[g（i-3j）?g（i-2j）?g（i-1j）g（ij）?g（i+1j）?g（i+2j）?g（i+3j）;
???????g（i-3j+2）?g（i-2j+1）?g（i-1j+1）g（ij）??g（i+1j-1）?g（i+2j-1）??g（i+3j-2）;
???????g（i-3j+3）??g（i-2j+2）?g（i-1j+1）g（ij）?g（i+1j-1）?g（i+2j-2）?g（i+3j-3）;
???????g（i-1j+3）??g（i-1j+2）?g（ij+1）?g（ij）g（ij-1）?g（i+1j-2）?g（i+1j-3）;
???????g（ij-3）?g（ij-2）?g（ij-1）?g（ij）g（ij+1）?g（ij+2）?g（ij+3）;
???????g（i-1j-3）?g（i-1j-2）??g（ij-1）?g（ij）g（ij+1）?g（i+1j+2）?g（i+1j+3）;
???????g（i-3j-3）?g（i-2j-2）?g（i-1j-1）?g（ij）g（i+1j+1）??g（i+2j+2）?g（i+3j+3）;
???????g（i-3j-2）?g（i-2j-1）?g（i-1j+1）?g（ij）g（i+1j+1）?g（i+2j+1）?g（i+3j+2）];
???????B=double（B‘）;
???????Dij=Dij+std（B）;
%end
%end
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%根據方向，確定濾波方向，進行中值濾波，輸入新的矩陣g1中。
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
[kkk]=find（Dij==min（Dij））;

Direction（ij）=kk;

??switch?（kk（1））%根據方向確定入口
??case?1
????
?????????????????B1=[g（ij-3）?g（ij-2）?g（ij-1）?g（ij）g（ij+1）?g（ij+2）?g（ij+3）];
?????????????????B1=double（B1）;
?????????????????g1（ij）=median3（B1）;%進行中值濾波
???????????

?case?2
???????
?????????????????B1=[g（ij-3）?g（ij-2）?g（ij-1）?g（ij）g（ij+1）?g（ij+2）?g（ij+3）];
?????????????????B1=double（B1）;
?????????????????g1（ij）=median3（B1）;
??????????????

?case?3
??????
?????????????????B1=[g（i-3j+3）??g（i-2j+2）?g（i-1j+1）g（ij）?g（i+1j-1）?g（i+2j-2）?g（i+3j-3）];
?????????????????B1=double（B1）;
?????????????????g1（ij）=median3（B1）;
????????????????

?case?4
???????
?????????????????B1=[g（i-1j+3）??g（i-1j+2）?g（ij+1）?g（ij）g（ij-1）?g（i+1j-2）?g（i+1j-3）];
?????????????????B1=double（B1）;
?????????????????g1（ij）=median3（B1）;
?????????????????

?case?5
???????
??????????????????B1=[g（ij-3）?g（ij-2）?g（ij-1）?g（ij）g（ij+1）?g（ij+2）?g（ij+3）];
??????????????????B1=double（B1）;
??????????????????g1（ij）=median3（B1）;
?????????????????

?case?6
??????
????
??????????????????B1=[g（ij-3）?g（ij-2）?g（ij-1）?g（ij）g（ij+1）?g（ij+2）?g（ij+3）];
??????????????????B1=double（B1）;
??????????????????g1（ij）=median3（B1）;
????????????????

?case?7
??????
??????????????????B1=[g（ij-3）?g（ij-2）?g（ij-1）?g（ij）g（ij+1）?g（ij+2）?g（ij+3）];
??????????????????B1=double（B1）;
??????????????????g1（ij）=median3（B1）;
????????????

?case?8
??????
??????????????????B1=[g（ij-3）?g（ij-2）?g（ij-1）?g（ij）g（ij+1）?g（ij+2）?g（ij+3）];
??????????????????B1=double（B1）;
??????????????????g1（ij）=median3（B1）;
???????????????????
??????
?end

end
end

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%將文件重新存回原目錄
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

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

?????文件???????3197??2008-02-10?22:20??spinfilter.m

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

?????????????????3197????????????????????1