在Matlab中实现Hough变换检测直线

mytomorrow

, Q+ s# ^: y' @5 z4 ZHough变换是图像处理中从图像中识别几何形状的基本方法之一。Hough变换的基本原理在于利用点与线的对偶性，将原始图像空间的给定的曲线通过曲线表达形式变为参数空间的一个点。这样就把原始图像中给定曲线的检测问题转化为寻找参数空间中的峰值问题。也即把检测整体特性转化为检测局部特性。比如直线、椭圆、圆、弧线等。
8 @5 K" C. T4 B; a1 |6 N' o4 x下面给出在matlab中用Hough变化检测直线的实现，文章末尾提供相关资源的下载。
其中先用log(laplacian of gaussian)算子检测图像边界：

) W1 F" h0 l+ r5 E" E* i- t; B! K01        sourcePic=imread('D: matlab程序Hough变换检测直线curve.bmp');
( \( Q. B8 \. g2 K# l* I02        [height,width,l]=size(sourcePic);
2 D' r* y* S1 U" g, ?, v! w03        grayPic=zeros(height,width);
04        for i=1:height %转换成灰度图象
1 S. i7 C, c6 `' f: `' O( L1 O. k05            for j=1:width
1 B: f$ ]5 u0 |4 H06                grayPic(i,j)=(sourcePic(i,j,1) sourcePic(i,j,2) sourcePic(i,j,3))/3;
07            end
08        end
09        logNum=0; %log算子的计算值
10        edgePic=zeros(height,width);
11        for i=3:height-2 %log算子提取边界
12            for j=3:width-2
13                logNum=16*grayPic(i,j)-grayPic(i-2,j)-grayPic(i-1,j-1)-2*grayPic(i-1,j)-grayPic(i-1,j 1)-grayPic(i,j-2)-2*grayPic(i,j-1)-2*grayPic(i,j 1)-grayPic(i,j 2)-grayPic(i 1,j-1)-2*grayPic(i 1,j)-grayPic(i 1,j 1)-grayPic(i 2,j);
) a, N8 I2 E) ^* K1 g6 J14                if(logNum > 0)
15                    edgePic(i,j)=255;
16                else
; ?: g' |# p( [0 y6 h' A. o17                    edgePic(i,j)=0;
18                end
19            end
: q( _( T( L  ~& q. j4 U+ U20        end
0 r9 R! b; i5 g( e21        % Hough变换检测直线，g(x)=(a,p)为边界点对应的平面
* ]" Z! n" ?0 z& J$ ?3 I22        ma=180; %a的值为0到180度
23        mp=round(sqrt(height^2 width^2)); %对应P的最大值
24        npc=zeros(ma,2*mp); %用于记录(a,p)对应的点的个数
( F* S, _4 |" E25        npp=cell(ma,2*mp); %用于记录(a,p)对应的点的坐标
26        for i=1:height %计算(a,p)的值，并做相应记录
27            for j=1:width
28                if(edgePic(i,j)==255)
29                    for k=1:ma
30                        p = round(i*cos(pi*k/180) j*sin(pi*k/180));
' [. g) T4 Z8 d7 w4 ]9 k+ c+ V, r31                        if(p > 0)
32                            %npc(k,mp*2)=npc(k,mp*2) 1;
% u9 q1 b) W5 e33                            %npp{k,mp*2}=[npp{k,mp*2},[i,j]'];
34                             10.01.05修正
2 u' a$ m: ?& v4 N35                            npc(k,mp p)=npc(k,mp p) 1;
# d) ?, e+ }+ N# x) V$ K36                            npp{k,mp p}=[npp{k,mp p},[i,j]'];
37                        else
38                            ap=abs(p) 1;
39                            npc(k,ap)=npc(k,ap) 1;
40                            npp{k,ap}=[npp{k,ap},[i,j]'];
41                        end
& k; e* K; M3 l' S) \  U42                    end
( r' T* Y. U9 h2 M  X# C8 v43                end
* K6 z7 [2 c' G7 N, w  f& D. C44            end
45        end
1 a4 k- k: X6 v, x( w* T% q46        for i=1:ma %根据(a,p)对应的点的个数，用红色标出相应直线
& M3 l; x2 r% K: P, U, W- I! @47            for j=1:mp*2
! N. a; i; {& X( H; r  }48                if(npc(i,j) > 90) %将提取直线的阈值设为90
  `! \9 c9 f( C0 ^; J2 i49                    lp=npp{i,j};
/ c5 ?) x$ J- f) W6 G8 k( }50                    for k=1:npc(i,j)
51                        sourcePic(lp(1,k),lp(2,k),1)=255;
52                        sourcePic(lp(1,k),lp(2,k),2)=0;
53                        sourcePic(lp(1,k),lp(2,k),3)=0;
) Q9 W: S. I0 x1 P- i& J1 g% v5 K54                    end
# b- Z" K( v. q' g5 E55                end
" c+ d1 i! ^9 h) {  \56            end
57        end
58        figure,imshow(sourcePic);

NingW

在Matlab中实现Hough变换检测直线