%Frost?filter?for?speckle?noise?reduction?
%Author?:?Jeny?Rajan?
function?[ft]?=?frost（I）?
%?I?is?the?noisy?input?image?
[x?y?z]?=?size（I）;?
I?=?double（I）;?
K?=?1;?
N?=?I;?
for?i?=?1:x?
????for?j?=?1:y???????????????????????????????
????????if?（i>1?&?i1?&?j????????????mat（1）?=?I（i-1j）;?
????????????mat（2）?=?I（i+1j）;?
????????????mat（3）?=?I（ij-1）;?
????????????mat（4）?=?I（ij+1）;?
????????????d（1）?=?sqrt（（i-（i-1））^2）;?
????????????d（2）?=?sqrt（（i-（i+1））^2）;?
????????????d（3）?=?sqrt（（j-（j-1））^2）;?
????????????d（4）?=?sqrt（（j-（j+1））^2）;?
????????????mn?=?mean（mean（mat））;?
????????????c?=?mat-mn;?
????????????c2?=?c.^2;?
????????????c3?=?c/（c2+.0000001）;?
????????????Cs?=?0.25*sum（sum（c3））;?
????????????m（1）?=?exp（-K*Cs*d（1））;?
????????????m（2）?=?exp（-K*Cs*d（2））;?
????????????m（3）?=?exp（-K*Cs*d（3））;?
????????????m（4）?=?exp（-K*Cs*d（4））;?
????????????ms?=?sum（sum（m））;?
????????????mp?=?m/ms;?
????????????N（ij）?=?sum（sum（mp.*mat））;?????????????????????
????????end?
?????end?
end?
ft?=?uint8（N）;
end

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

?????文件???????1054??2018-04-21?12:12??frost.m

?????文件???????2289??2018-04-21?12:10??kuan.m

?????文件???????2704??2018-04-21?12:03??lee.m

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

?????????????????6047????????????????????3