DIP数字图像处理（一）图像分割：霍夫变换、阈值处理及其Matlab手动实现（不调库）

饼干妹妹324

已于 2022-12-12 15:28:45 修改

阅读量475

点赞数

分类专栏： DIP 文章标签： matlab 算法图像处理计算机视觉

于 2022-12-07 14:50:13 首次发布

本文链接：https://blog.csdn.net/SolvingMaster/article/details/128214885

版权

DIP 专栏收录该内容

2 篇文章 0 订阅

订阅专栏

图像分割算法目录

点、线和边缘检测
- 霍夫变换
阈值处理

点、线和边缘检测

霍夫变换

考虑直线的法线表示 $xcos\theta+ysin\theta=\rho$ ，及其参数空间 $(\rho, \theta)$ . 其中 $-90^o\leq \theta\leq 90^o$ ， $-D\leq \rho\leq D$ ， $D$ 是图像对角最大距离。

算法步骤

将参数空间 $(\rho, \theta)$ 划分为有限数量的累加器，累加器的值初始化为0.

对于 $x y$ 空间的每个非背景点，令 $\theta$ 等于 $\theta$ 轴上每个允许的值，利用公式 $\rho=xcos\theta+ysin\theta$ 解出round( $\rho$ )

将得到的 $(\rho, \theta)$ 对应的adder+=1

最终，累加器 $(\rho, \theta)$ 的值等价于 $x y$ 空间中位于该直线上点的数量。算法复杂度 $O (n)$ .

Matlab实现的注意点

需要先对图像做阈值处理化为0-1图像

初始化参数空间 $(\rho, \theta)$ 时，注意矩阵大小应为img_hough = zeros(2*(D/dD)+1, 2*(90/dtheta)+1) 如下图所示

提取直线参数做下游任务时，注意取出二维矩阵最大元素索引的方法为

[val, idx] = max(img_hough);
[val2, idx_col] = max(val);
idx_row = idx(idx_col);
idx_col = idx_col;

霍夫变换代码

close all;clear;clc;
img = imread('./images/tilted.PNG');
imshow(img,[]);
[m, n, c] = size(img);

% thresholding
img_binary = zeros(m,n);

for i = 1:m
    for j = 1:n
        if img(i,j)<250
            img_binary(i,j)=1;
        else
            img_binary(i,j)=0;
        end
    end
end
imshow(uint8(img_binary),[]);

% hough
D = m+n;
dD = 1;
dtheta = 1;
img_hough = zeros(2*(D/dD)+1, 2*(90/dtheta)+1);
for i=1:m
for j=1:n
    if img_binary(i,j)>0.5
        for theta = -90:90
            theta_rad = theta/180*pi;
            rho = i*cos(theta_rad)+j*sin(theta_rad);
            rho = round(rho);
            index1 = rho+D+1;
            index2 = theta+90+1;
            img_hough(index1, index2)=img_hough(index1, index2)+1;
        end
    end
end
end
imshow(img_hough,[]);

% argmax and its index
[val, idx] = max(img_hough);
[val2, idx_col] = max(val);
idx_row = idx(idx_col);

% rotated angle
angle = idx_col-90-1;
img_final = imrotate(img, -angle, 'bilinear');
imshow(img_final, []);