用星点法测图像MTF曲线matlab代码

以下是一个简单的Matlab代码，用于计算和绘制图像的MTF曲线。其中，我们使用星点法来测量MTF曲线。

% 读取图像
img = imread('test_image.jpg');

% 转换成灰度图
img_gray = rgb2gray(img);

% 计算傅里叶变换
fft_img = fft2(double(img_gray));

% 计算幅度谱
amp_spec = abs(fft_img);

% 计算相位谱
phase_spec = angle(fft_img);

% 计算MTF曲线
[ny, nx] = size(img_gray);
x = 0:nx/2-1;
y = 0:ny/2-1;
[X, Y] = meshgrid(x, y);
r = sqrt(X.^2 + Y.^2);
r = round(r(:));
amp_spec = round(amp_spec(:));
mtf = accumarray(r+1, amp_spec) ./ accumarray(r+1, ones(size(amp_spec)));
mtf = mtf(1:nx/2);

% 绘制MTF曲线
figure;
plot(x, mtf, 'o-');
xlabel('Spatial Frequency (lp/mm)');
ylabel('MTF');
title('MTF Curve using Star Chart Method');

代码中，我们首先读取图像并将其转换为灰度图像。然后，我们计算图像的傅里叶变换，并计算幅度谱和相位谱。接着，我们使用星点法测量MTF曲线，并绘制结果。