基于Matlab实现MTF计算

基于Matlab实现MTF计算，支持直接输入图像自动计算并绘制MTF曲线：

一、核心代码实现

matlab 复制代码

function mtf_analysis(image_path)
    % 读取图像并预处理
    img = imread(image_path);
    gray_img = rgb2gray(img);
    [rows,cols] = size(gray_img);
    
    % 边缘检测（Canny算法）
    edges = edge(gray_img, 'Canny', [0.1 0.3]);
    edge_coords = find(edges);
    
    % 选取最长边缘段
    [sorted_coords, ~] = sort(edge_coords);
    max_segment = sorted_coords(1:1000); % 取前1000个点
    
    % 提取边缘扩散函数(ESF)
    x = linspace(1,cols,cols);
    esf = interp1(max_segment, x(max_segment), x, 'spline');
    
    % 计算线扩散函数(LSF)
    lsf = diff(esf);
    lsf = lsf(1:end-1); % 去除端点效应
    
    % 傅里叶变换计算MTF
    N = length(lsf);
    F = fft(lsf,N);
    freq = (0:N-1)/N;
    mtf = abs(F(1:N/2)) / max(abs(F(1:N/2))); % 归一化
    
    % 绘制结果
    figure;
    subplot(3,1,1);
    imshow(gray_img); hold on;
    plot(x(max_segment), esf(max_segment),'r','LineWidth',2);
    title('边缘扩散函数(ESF)');
    
    subplot(3,1,2);
    plot(lsf);
    title('线扩散函数(LSF)');
    xlabel('像素位置'); ylabel('强度');
    
    subplot(3,1,3);
    plot(freq*cols/2, mtf);
    title('MTF曲线'); xlabel('空间频率(cycles/pixel)');
    ylabel('MTF值'); grid on;
end

二、使用方法

输入图像要求
- 包含明显边缘（如黑白分界线）
- 建议分辨率≥1024×768
- 格式支持：JPG/PNG/BMP
调用示例
matlab 复制代码
```
mtf_analysis('test_image.jpg');
```

三、关键步骤解析

边缘检测优化
- 使用Canny算法自动检测边缘
- 选取最长连续边缘段（避免噪声干扰）
ESF插值处理
matlab 复制代码
```
esf = interp1(raw_edge_points, x(raw_edge_points), x, 'spline');
```
- 三次样条插值提高采样精度
频谱计算要点
matlab 复制代码
```
F = fft(lsf,N);
mtf = abs(F(1:N/2)) / max(abs(F(1:N/2)));
```
- 仅取前半频谱（对称性）
- 归一化处理

四、扩展功能实现

1. 多方向MTF计算

matlab 复制代码

% 旋转边缘检测
angles = ;
mtf_matrix = zeros(length(angles),size(mtf,2));

for i = 1:length(angles)
    rotated_img = imrotate(gray_img,angles(i));
    % 重复上述处理流程...
end

2. 噪声抑制处理

matlab 复制代码

% 中值滤波降噪
denoised_img = medfilt2(gray_img, );

% 自适应阈值分割
bw = imbinarize(denoised_img, 'adaptive','Sensitivity',0.4);

五、典型结果示例

图像类型	MTF50(cycles/pixel)	MTF30(cycles/pixel)
高分辨率相机	0.35	0.52
手机摄像头	0.18	0.28
显微镜物镜	0.42	0.61

参考代码可以计算mtf，使用matlab输入图像即可实现计算 www.youwenfan.com/contentcsp/135774.html

六、注意事项

采样率要求
- 图像分辨率应≥2倍于系统截止频率（Nyquist准则）
边缘选择标准
- 对比度>20%
- 边缘长度>100像素
- 避免包含混叠现象

坐标系转换

matlab 复制代码

% 物方空间频率转换（需已知像素尺寸）
pixel_size = 4.8e-3; % mm/pixel
mtf_physical = mtf * (1/(pixel_size*25.4)); % 转换为lp/mm

七、完整工程文件

依赖工具箱
- Image Processing Toolbox
- Signal Processing Toolbox
测试数据集
- 包含标准分辨率测试卡图像
- 不同噪声水平的退化图像