给我一段空间频率评价灰度图像的matlab代码

以下是一个可以评价灰度图像的空间频率的MATLAB代码示例：

function freq_score = compute_spatial_frequency(img)
    % 计算图像的傅里叶变换
    img_fft = fft2(double(img));
    img_fft_shift = fftshift(img_fft);

    % 计算图像的频率谱
    freq_spectrum = abs(img_fft_shift);
    
    % 计算频率谱的中心位置
    [M, N] = size(freq_spectrum);
    cx = round(N/2);
    cy = round(M/2);
    
    % 计算频率谱的平均值和标准差
    freq_mean = mean2(freq_spectrum);
    freq_std = std2(freq_spectrum);
    
    % 计算频率谱的方差
    freq_var = sum(sum((freq_spectrum - freq_mean).^2)) / (M * N);
    
    % 计算空间频率得分
    freq_score = freq_var / freq_std;
end

在这个函数中，我们首先对输入图像进行傅里叶变换，然后计算其频率谱。接着，我们计算频率谱的平均值、标准差和方差，并将方差除以标准差得到空间频率得分。空间频率得分是一种用于评价图像的空间频率特征的指标，其值越小表示图像的频率特征越窄，即图像的空间频率越集中。