图像椭圆定位算法故障求助：代码修复与中心偏移计算需求

原代码问题分析 #

1. 滤波函数冗余且过度模糊：原getfilter中imgaussfilt(im,10)的sigma值过大，会严重模糊椭圆边缘；两次imresize操作（下采样再上采样）会导致图像像素化，丢失椭圆细节，不利于后续检测。
2. 检测函数逻辑错误：
   • 用imm5=int16(im>1);覆盖了之前的二值化结果im2bw(im,Thr)，逻辑完全错误，应该基于二值化后的图像处理。
   • 过时的hist函数应替换为histcounts，避免版本兼容性问题。
   • 未对检测到的区域进行筛选，无法区分椭圆和噪点；循环调整结构元素大小的逻辑容易陷入死循环或误判。
   • 参数itrun未被使用，属于冗余定义。

修复后的完整代码 #

1. 优化后的滤波函数 #

function [im_new] = getfilter(im)
    % 适度高斯模糊抑制噪点，保留椭圆边缘
    imm1 = imgaussfilt(im, 2);
    % 中值滤波去除椒盐噪点
    imm2 = medfilt2(imm1, [3, 3]);
    im_new = imm2;
end

2. 修复后的椭圆检测函数 #

function [s, Thr] = findeclipse(im)
    % 用histcounts替代过时的hist函数
    [counts, bins] = histcounts(im(:), 0:255);
    Thr = otsuthresh(counts);
    % 正确二值化图像
    imm5 = imbinarize(im, Thr);
    % 形态学开运算去除小噪点（先腐蚀后膨胀）
    SE = strel('disk', 2);
    BW3 = imopen(imm5, SE);
    
    % 提取区域属性，增加面积筛选排除噪点
    s = regionprops(BW3, ...
        'Centroid', ...
        'MajorAxisLength', ...
        'MinorAxisLength', ...
        'Orientation', ...
        'Area');
    
    % 筛选符合椭圆特征的区域（面积足够大、长短轴比例合理）
    valid_idx = [];
    for i = 1:length(s)
        % 面积阈值可根据图像调整，这里针对目标图像设置
        if s(i).Area > 1000 && (s(i).MajorAxisLength / s(i).MinorAxisLength) < 3
            valid_idx = [valid_idx, i];
        end
    end
    
    if isempty(valid_idx)
        error('未检测到椭圆');
    elseif length(valid_idx) > 1
        error('检测到多个疑似椭圆区域');
    else
        s = s(valid_idx);
    end
end

3. 主调用脚本（计算偏移量） #

% 读取图像
im = imread('cAR5q.png');
% 预处理滤波
im_filtered = getfilter(im);
% 检测椭圆
[s, ~] = findeclipse(im_filtered);
% 计算图像中心
img_center = [size(im, 2)/2, size(im, 1)/2];
% 计算椭圆中心与图像中心的偏移量
offset = s.Centroid - img_center;

% 输出结果
fprintf('椭圆中心坐标：(%.2f, %.2f)\n', s.Centroid(1), s.Centroid(2));
fprintf('图像中心坐标：(%.2f, %.2f)\n', img_center(1), img_center(2));
fprintf('偏移量（水平, 垂直）：(%.2f, %.2f)\n', offset(1), offset(2));

% 可视化结果（可选）
figure; imshow(im); hold on;
plot(s.Centroid(1), s.Centroid(2), 'r+', 'MarkerSize', 15);
plot(img_center(1), img_center(2), 'bo', 'MarkerSize', 10);
legend('椭圆中心', '图像中心');
title('椭圆定位与中心偏移');

修复说明 #

• 滤波函数调整了高斯模糊的sigma值和中值滤波窗口大小，在抑制噪点的同时保留椭圆边缘细节。
• 检测函数修复了二值化逻辑错误，替换了过时函数，增加了区域筛选条件（面积、长短轴比例），确保只检测目标椭圆。
• 主脚本实现了偏移量的计算，并添加了可视化功能便于验证结果。

内容的提问来源于stack exchange，提问作者Nir