星点法测MTF图像matlab代码

以下是使用星点法测量MTF图像的MATLAB代码示例：

% 读取MTF图像
img = imread('MTF_Image.png');

% 显示MTF图像
figure;
imshow(img);
title('MTF Image');

% 提取灰度值通道
img = rgb2gray(img);

% 计算星点位置
corners = detectHarrisFeatures(img);

% 显示星点
figure;
imshow(img);
hold on;
plot(corners.selectStrongest(10));
title('Harris Corners');

% 计算星点对应的模糊尺寸
blur_sizes = zeros(size(corners,1),1);
for i=1:size(corners,1)
    [blur_sizes(i), ~] = measure_blur(img, corners.Location(i,:));
end

% 显示模糊尺寸分布
figure;
histogram(blur_sizes);
xlabel('Blur Size');
ylabel('Count');
title('Blur Size Distribution');

% 计算MTF
mtf = 1 - blur_sizes / max(blur_sizes);

% 显示MTF曲线
figure;
plot(mtf);
xlabel('Spatial Frequency');
ylabel('MTF');
title('MTF Curve');

% 计算MTF50
mtf50 = find(mtf < 0.5, 1, 'last');

% 显示MTF50位置
hold on;
plot([mtf50 mtf50], [0 0.5], 'r--');
text(mtf50, 0.55, ['MTF50 = ' num2str(mtf50)]);

其中`measure_blur`函数用于计算给定位置的模糊尺寸，可以使用以下代码实现：

function [blur_size, blur_image] = measure_blur(img, center, patch_size)
% 默认取15x15的局部区域
if nargin < 3
    patch_size = 15;
end

% 提取局部区域
patch = imcrop(img, [center-patch_size/2 patch_size-1]);

% 计算Laplacian值
laplacian = abs(del2(patch));

% 取局部区域的平均Laplacian值作为模糊度量
blur_size = mean2(laplacian);

% 可选输出模糊区域图像
if nargout > 1
    blur_image = laplacian > blur_size;
end
end

需要注意的是，该方法的精度受到图像分辨率和星点密度的影响，因此需要根据实际情况进行调整。