matlab环境下狭缝法测图像MTF曲线代码

以下是一个简单的 Matlab 程序来计算图像的 MTF 曲线，使用了狭缝法：

% 狭缝法测图像MTF曲线
clc; clear; close all;

% 读入图像
img = imread('test.jpg');

% 灰度化
if size(img, 3) == 3
    img = rgb2gray(img);
end

% 显示原始图像
figure(1);
imshow(img);
title('Original Image');

% 水平方向狭缝宽度
slit_width = 20;

% 每个狭缝位置的像素值
slit_values = zeros(size(img, 1), size(img, 2) - slit_width);
for i = 1:size(img, 1)
    for j = 1:size(img, 2) - slit_width
        slit_values(i, j) = mean(img(i, j:j+slit_width-1));
    end
end

% 每个像素值的傅里叶变换
fft_values = abs(fftshift(fft2(slit_values)));

% 计算MTF
mtf = zeros(size(fft_values, 2), 1);
for i = 1:size(fft_values, 2)
    mtf(i) = mean(fft_values(:,i));
end
mtf = mtf / mtf(1);

% 画MTF曲线
figure(2);
plot(mtf);
title('MTF Curve');
xlabel('Spatial Frequency (cycles/pixel)');
ylabel('Contrast');

在这个程序中，我们首先读入图像并将其灰度化。然后，我们定义了狭缝的宽度，计算每个狭缝位置的平均像素值，并使用傅里叶变换计算每个像素值的频谱。最后，我们计算 MTF 并绘制 MTF 曲线。

请注意，这只是一个简单的程序，您可能需要根据您的具体应用程序进行修改。