计算图像Spatial Frequence的Matlab程序SF.m

可以计算图像的行与列的空间频率以及总频率

%{
@ 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;

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

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

----------- ---------  ---------- -----  ----

                  548                    1