基于Matlab实现MTF计算

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

---

一、核心代码实现

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

---

二、使用方法

1. 输入图像要求

   • 包含明显边缘（如黑白分界线）
   • 建议分辨率≥1024×768
   • 格式支持：JPG/PNG/BMP

2. 调用示例

   mtf_analysis('test_image.jpg');
   

---

三、关键步骤解析

1. 边缘检测优化

   • 使用Canny算法自动检测边缘
   • 选取最长连续边缘段（避免噪声干扰）

2. ESF插值处理

   esf = interp1(raw_edge_points, x(raw_edge_points), x, 'spline');
   

   • 三次样条插值提高采样精度

3. 频谱计算要点

   F = fft(lsf,N);
   mtf = abs(F(1:N/2)) / max(abs(F(1:N/2)));
   

   • 仅取前半频谱（对称性）
   • 归一化处理

---

四、扩展功能实现

1. 多方向MTF计算

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

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

2. 噪声抑制处理

% 中值滤波降噪
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/contentcnp/135774.html

六、注意事项

1. 采样率要求

   • 图像分辨率应≥2倍于系统截止频率（Nyquist准则）

2. 边缘选择标准

   • 对比度>20%
   • 边缘长度>100像素
   • 避免包含混叠现象

3. 坐标系转换

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

---

七、完整工程文件

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