编写基于matlab的sfr函数

以下是基于MATLAB的SFR函数的示例代码：

function [resolution, sfr, freq] = calculate_sfr(image, threshold, fft_size, angle)
% 计算图像的分辨率
%
% 输入参数：
% image - 输入图像，为灰度图像
% threshold - Sobel算子的阈值，默认为0.1
% fft_size - 傅里叶变换的大小，默认为512
% angle - SFR算法的方向，0表示水平方向，90表示垂直方向，默认为0
%
% 输出参数：
% resolution - 图像分辨率
% sfr - SFR曲线
% freq - 空间频率

% 设置默认参数
if nargin < 4
    angle = 0;
end
if nargin < 3
    fft_size = 512;
end
if nargin < 2
    threshold = 0.1;
end

% 计算边缘
if angle == 0
    edge_image = edge(image, 'sobel', threshold, 'horizontal');
else
    edge_image = edge(image, 'sobel', threshold, 'vertical');
end

% 计算SFR
sfr = sum(edge_image, angle+1)';

% 对SFR进行傅里叶变换
sfr_fft = abs(fft(sfr, fft_size));
sfr_fft = sfr_fft(1:fft_size/2);

% 计算空间频率
freq = (0:(fft_size/2-1)) / fft_size;

% 找到SFR曲线的第一个峰值
[pks, locs] = findpeaks(sfr_fft);
[~, max_idx] = max(pks);
resolution = freq(locs(max_idx));

% 绘制SFR曲线
plot(freq, sfr_fft);
xlabel('空间频率');
ylabel('幅度');
title(sprintf('分辨率 = %.2f', resolution));

使用示例：

image = imread('lena.png');
gray_image = rgb2gray(image);
[resolution, sfr, freq] = calculate_sfr(gray_image);

需要注意的是，该函数仅适用于灰度图像。如果要对彩色图像进行分辨率计算，需要先将其转换为灰度图像。此外，SFR算法的参数需要根据图像的特性进行适当调整。