%{
@?calculate?the?spatial?freuency?of?an?image
%}
function?t?=?SF（A）

%[M?N]?:?image?size
%RF:?row?frequency
%CR:?col?frequency

A?=?double（A）;
[MN]?=?size（A）;

%?Calculate?RF
RF?=?0.0;
for?i=1:M
????for?j=2:N
????????RF?=?RF?+?（A（ij）?-?A（ij-1））^2;
????end
end
RF?=?sqrt（RF/（M*N））;
%RF?=?（RF/（M*N））;

%?Calculate?CF
CF=0.0;
for?i=2:M
????for?j=1:N
????????CF?=?CF??+?（A（ij）?-?A（i-1j））^2;
????end
end
CF?=?sqrt（CF/（M*N））;
%CF?=?（CF/（M*N））;

%?Calculate?SF
SF?=?sqrt（RF^2?+?CF^2）;
%SF?=?（RF^2?+?CF^2）;
t?=?SF;

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

?????文件????????548??2008-10-29?13:04??SF.m

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

??????????????????548????????????????????1