最近在用STM32H7这块板子做机器视觉的学习，在这里记录下学习心得与过程。有不对的地方欢迎大家指正。（stm32H7开发板，lcd9332显示屏）
(一）边缘检测概念
: L# @+ E- [  l) i: j9 ?) y: Z边缘是图像的基本特征，图像的边缘通常意味着一个区域的结束与另一个区域的开始，对于图像后续处理与检测有重要作用。
在图像处理中，边缘通常被认为是周围的像素灰度发生阶跃变化或者屋顶变化的像素集合。
. N) r2 i* ]0 I/ X" S6 l边缘的特征是沿边缘走向的像素变化平缓，垂直与边缘走向的像素变化剧烈。
2 }- r& n0 S* R边缘检测就是采用某种算法来提取出图像中灰度剧烈变化的区域边界，并且将其来接起来的方法。
（二）Canny边缘检测
# a/ |$ [2 G! e: o, VCanny算子检测是目前使用较为广泛的边缘检测方法之一，本文也是通过canny算子来实现边缘检测的。
其主要步骤为：高斯滤波–>一阶差分算子计算梯度幅值及其方向–>对梯度幅值进行极大值抑制–>双阈值提取边缘。

1. 基于全局阈值的二值化处理
假设某一副灰度图有如下的直方图，该图像由暗色背景下的较亮物体组成，从背景中提取这一物体时，将阈值T作为分割点，分割后的图像g(x, y)由下述公式给出，称为全局阈值处理

' n4 G1 u. J/ ]- T; F: M; L5 b. i
8 O( j/ k2 W9 a0 c6 V基本全局阈值处理方法

为全局阈值T选择一个初始的估计值

用T分割图像，产生两组像素：G1由大于T的像素组成，G2由小于T的像素组成

0 n6 i$ Z% K* s) R1 j对G1和G2的像素分别计算平均灰度值m1和m2
; ?, m1 F6 I) }2 X
计算新的阈值T = 1/2 * (m1 + m2)
. u& V" {: J- f$ B- U  K; f# L
1 o. e0 i1 E# q. ]6 f3 {重复步骤2-4，直到连续迭代中的T值差小于一个预定义的参数ΔT
代码如下：

1. while(1){
   2 i3 t* f5 o9 C) W' {
2.                         T = 0.5 * (m1 + m2);
   
3.                         m1 = 0.0;
   
4.                         m2 = 0.0;
   ( y# l: k* M! T) `, |4 J
5.                         m1_num = 0;
   - g, r" t$ m0 X: w( U- |
6.                         m2_num = 0;
   
7. 
   
8.                         for (i = 0; i < len; i++){
   
9.                                 if (GetImageData<i> </i><= T){
   
10.                                                 m1 += GetImageData;
    ; ~; s; n: [5 e9 n& [% a% t2 @
11.                                                 m1_num++;
    ) x# ]7 g# g3 t" m3 G/ R2 M
12.                                         }
    
13.                                 else{
    
14.                                                 m2 += GetImageData;
    
15.                                                 m2_num++;
    " X# ^' y/ d2 S2 O
16.                                                 }
    9 \9 C; O# W' Z) F$ _+ [) A# V- F
17. //                                GetImageData+=1;
    ) C; w% V8 D3 {3 U
18. //                                DataNum++;
    . k) y: \1 K6 |) l+ ^; }
19.                         }
    5 s" l5 S1 |& m( }/ K: \# P; U
20.                         if (m1_num != 0)
    # o& r, f2 l" }! l* w
21.                                         m1 /= m1_num;
    
22.                         if (m2_num != 0)
    
23.                                         m2 /= m2_num;
    " u% c- i' l! h4 g) N9 z
24.                         if(fabs(T-(0.5*(m1 + m2)))<delt_t)
    9 J1 ?+ i: {$ J4 u7 I
25.                                 break;
    
26.         }
    , M% [) s$ {; ?: C
27.         GetImageData-=DataNum;
    
28.         for(j=0;j<len;j++)
    * G6 i+ M; l; O  w
29.         {
    # p2 O! g4 r) F# e
30.                 if(GetImageData[j]<=T)
    
31.                         GetImageData[j]=0x00;
    % r9 w/ p! q2 w& t( m* I
32.                 else
    9 C" Y5 ]$ x( Q1 R  F; T* Q' ~
33.                         GetImageData[j]=0xFF;        

复制代码

原图：
/ ~1 g$ S: H& F" L& d
1 S3 p  v3 \* D0 {7 B" R! |# W" ?

9 |/ \+ |( }% d4 u% W4 V
/ R+ y  {* g4 h7 p二值化之后的图片：

我在刚开始做二值化的时候直接给定阈值进行处理，每张图片处理效果都不一样，很多时候会出现去不掉阴影的情况，此处通过参考阈值处理后达到所需要的效果。
4 s& A$ U" E! U. h0 r; r; `
2. 高斯滤波
进行简单的二值化处理之后开始Canny边缘检测的第一部分：高斯滤波。
滤波的主要目的是降噪，一般的图像处理算法都需要先进行降噪。而高斯滤波主要使图像变得平滑。
高斯函数是一个类似与正态分布的中间大两边小的函数。
对于一个位置（x,y）的像素点，其灰度值（这里只考虑二值图）为f(x,y)。
那么经过高斯滤波后的灰度值将变为：

$ ^0 l) @  x% F2 J* D

$ ~# Z* ]; P; _; N& a简单说就是用一个高斯矩阵乘以每一个像素点及其邻域，取其带权重的平均值作为最后的灰度值。
4 v6 @9 x: S" T0 Q
如下图3x3的像素区域为例，此矩阵中，我们选择的原点（0，0）周围的八个点为其领域：
8 i6 X  R  T4 c) F

8 }. k/ w5 G- B  h4 C7 {
2 p' Q6 z: r3 t0 a" i" k' v! u我们设σ=1.5，将这个此区域的每个像素点做高斯变换之后就得到了以下的数据：



2 M; M8 R3 N* N9 E* ~% v此矩阵为其权重矩阵，这个区域的权重总和为0.4787147，我们需要让他们的权重和为1，于是每个点的值在除以0.4787147，此时我们重新得到了一个3x3的矩阵：
4 G: @7 y5 v+ `/ |

7 _6 `2 o5 Y% z0 H# K3 o9 M
有了相应的3x3矩阵的权重值，我们再对原图像素点进行卷积运算，简单来说就是原图中的点附近的3x3大小区域乘以高斯模板区域。
4 P; t( D" y4 I& \代码如下：

6 v2 b- C" D" A6 I1 E5 c( ]4 l

1. for(i=widthG+1;i<((heightG-1)*widthG)-2;i++)
   
2.         {
   
3.                         temp=(0.094*ImageNowg[i-widthG-1]+0.118*ImageNowg[i-widthG]+0.094*ImageNowg[i-widthG+1]
   2 u! s9 M& g% h: q/ p# v
4.                 +0.118*ImageNowg[i-1]+0.147*ImageNowg+0.118*ImageNowg[i+1]+0.094*ImageNowg[i+widthG-1]
   $ Q9 r* \2 |" a" m" p( s/ f2 V6 s
5.                 +0.118*ImageNowg[i+widthG]+0.094*ImageNowg[i+widthG+1]);
   
6.                         
   
7.                         gerytemp=(u8)temp;
   
8.         }
   - X8 f) J8 L, g7 E5 D, a& u& V( w
9.         for(j=widthG+1;j<((heightG-1)*widthG)-2;j++)
   
10.         {
    
11.                
    ; D, M; V% l. @" ]: C! k1 J
12.                         if(gerytemp[j]>150)
    # D6 x* Z# J+ O! H. K
13.                         {
    ' u: K- v+ I9 N$ P- n
14.                                         gerytemp[j]=255;        
    6 X' p9 _. D1 x: L- f6 q9 A! f
15.                         }
    4 l& g! I& t/ x1 i" D
16.                         ImageNowg[j]=gerytemp[j];
    
17.                         
    ( O8 I" N- d$ n* }
18.         }

复制代码

/ v! K' ~- y8 N; q7 k$ M0 m+ b) a. T- ^高斯模糊后的图片（由于显示像素我没有进行色彩转换，高斯模糊的部分不是灰色），其中蓝色的部分为高斯模糊的像素。

- \/ l, K7 b; F. V/ K- n( e' Z6 I3.Sobel算子与极大值抑制
6 o! V1 g9 B7 P! G, HSobel算子是图像处理中的基本算子之一，他是一离散差分算子，用于计算图像亮度函数的梯度近似值。
4 U- L7 H' _- y1 hSobel算子有两个，一个是水平算子，一个是垂直算子。由于Sobel算子是滤波算子公式，用于提取边缘，可以利用快速卷积函数，简单有效，在Canny边缘检测被作为常用算子之一。
Sobel算子也是一种梯度幅值，该算子包含两组3x3的矩阵，分别为横向与纵向，将之与 图像作平面卷积既可以得出横向及纵向的灰度差分近似值。假设原图像为I，横向边缘检测图像为G(x),纵向边缘检测图像为G(y)，模板卷积因子如下：

. f' O6 ]# {- d( [+ ~/ i

7 w& \; o( W4 Z
6 ^3 A% F0 V5 b3 P% @有了模板卷积因子，我们计算出每个像素点的横向纵向梯度近似值之后，需要知道梯度的幅值大小以及梯度放向，梯度幅值由以下公式得出：
4 e4 }0 w, v/ P" P+ L) qG=|G(x)+G(y)|
% L* N  p4 o0 M, p3 l" |梯度方向公式：
' {* {. r' r+ m& |" l' W0 h9 mΘ=artan(G(y)/G(x))
: P$ u# e+ q6 ?& `计算出梯度幅值之后，对图像进行门限处理，设定门限值为T，梯度幅值大于T时，规定该像素点的灰度值为0XFF，否则为0，即：
/ g9 j3 o; V" Z/ S) m7 l# y: \
& @; O& w+ X, O9 V

  N! \; A0 @& m非极大值抑制
在Sobel算子对图像进行处理之后，图像的边缘不能作为最佳边缘，边缘有且应当只有一个准确的响应。非极大值抑制作用是对于梯度不是最大值的点进行抑制，对梯度图像中每个像素进行非极大值抑制的算法是：
7 [. \$ g4 n: z$ v
1.将当前像素的梯度强度与沿正负梯度方向上的两个像素进行比较。
2.如果当前像素的梯度强度与另外两个像素相比最大，则该像素点保留为边缘点，否则该像素点将被抑制。
如下图所示，其中P点为当前像素点，P1,P2为领域像素梯度幅值比较点，在这三点当中，如果P点的幅值是最大值，则保留，如果不是则抑制：
2 J  r8 i+ S# P- z; c

为了方便计算，通常会在俩个相邻像素（G(i+1,j-1),G(i+1,j))，之间采用线性插值的方法计算P1，P2的梯度幅值。
具体计算方式如下：
$ h* D. T3 j" R% M* {- eP2=（1-tanΘ)xG(i+1,j)+tanΘxG(i+1,j-1);

4 q8 i- }; r6 a) zP1=（1-tanΘ)xG(i-1,j)+tanΘxG(i-1,j+1);
代码如下：

1. for(i=widthS+1;i<((widthS-1)*(heightS-1))-2;i++)
   
2.         {
   5 B" @: J  r; c% N- e' C7 t
3.                
   : O6 l7 n  _% V! c8 ?9 N4 U0 P1 |
4.                 i00=ImageNow[i-widthS-1];
   
5.                 i01=ImageNow[i-widthS];
   
6.                 i02=ImageNow[i-widthS+1];
   % G$ V9 ]( {7 c  ?4 W( [' P
7.                 i10=ImageNow[i-1];
   
8.                 i12=ImageNow[i+1];
   
9.                 i20=ImageNow[i+widthS-1];
   
10.                 i21=ImageNow[i+widthS];
    & O0 E/ ?, s, z  ?; |
11.                 i22=ImageNow[i+widthS+1];
    $ M. x  i5 n5 C( [' r
12.                
    
13.                
    ! W) I7 Q% S+ _+ M& C: W' ]. \
14.                 V=-i00-2*i01-i02+i20+2*i21+i22;
    
15.                 H=-i00+i02-2*i10+2*i12-i20+i22;
    
16.                
    5 Q2 X/ o( Y) j1 V1 t" W4 T. h9 Y
17.                 grey=(u8)abs(H)+(u8)abs(V);
    - g% x6 b* L7 i5 w- U; F
18.                
    & `! U& c9 R+ ?, W- U+ d
19.                         if(grey==0)
    6 p: r% h( d. f
20.                 {
    
21.                         ImageOperation=0x00;
    
22.                         continue;
    " [; |; k% Q+ I, ?& |3 X6 z6 I* }9 m
23.                 }
    
24.                 else
    7 r0 T) Q6 j# F
25.                 {
    / x' L! X, R; H
26.                         if(abs(H)>abs(V))///X///
    
27.                         {
    ! p, G6 e2 X2 I4 l2 ]; J* ]* L
28.                                 weight=abs(V)/abs(H);
    ! R8 S" ?$ m5 g' Q% Q4 U$ x& k
29.                                 if(H*V>0)
    
30.                                 {
    
31.                                         comp1=grey[i+1]*(1-weight)+grey[i-widthS+1]*weight;
    
32.                                         comp2=grey[i-1]*(1-weight)+grey[i+widthS-1]*weight;
    
33.                                 }
    % b) f# \/ Y$ @8 f
34.                                 else
    
35.                                 {
    ) @& Y  j& g# }1 J4 R% l
36.                                         comp1=grey[i+1]*(1-weight)+grey[i+widthS+1]*weight;
    
37.                                         comp2=grey[i-1]*(1-weight)+grey[i-widthS-1]*weight;
    4 Z# F7 y6 z# Y) A" [
38.                                 }
    
39.                         }
    
40.                         else/Y/
    4 q' J  o: n/ `; u; s
41.                         {
    
42.                                 weight=abs(H)/abs(V);
    
43.                                 if(H*V>0)
    5 {. Y  v( S7 Z% N% I
44.                                 {
    ) c, w+ w6 w4 F0 f7 j0 V$ r0 n
45.                                         comp1=grey[i-widthS]*(1-weight)+grey[i-widthS+1]*weight;
    
46.                                         comp2=grey[i+widthS]*(1-weight)+grey[i+widthS-1]*weight;
    " B+ e4 B: {+ j+ ]# W& F7 d
47.                                 }
    6 m& C! \: `: S, l* i5 c& {2 U
48.                                 else
    
49.                                 {
    ! v5 O, n- l' r& a# E
50.                                         comp1=grey[i+widthS]*(1-weight)+grey[i+widthS+1]*weight;
    5 m, l5 ^% {& g2 x) P
51.                                         comp2=grey[i-widthS]*(1-weight)+grey[i-widthS-1]*weight;
    % D  O+ G3 i9 Y
52.                                 }
    0 K' Z/ ~; _/ z+ `
53.                         }
    ( Y- X9 i/ }- [- m  w/ J- m
54. 
    
55.                                 if((grey>comp1)&&(grey>comp2))
    
56.                                         ImageOperation=0xff;
    
57.                                 
    5 J" ?2 R1 s; i
58.                                 else
    6 P2 F( M& |9 S( s, Z. x
59.                                 ImageOperation=0x00;
    - z0 B( K4 g- C  A3 i+ G; m+ v0 L0 L
60.                 }
    7 }  K- O" V& f4 n+ m4 f* D
61.                
    
62.         }
    3 D. D" i0 V& o2 O. o8 G# z" m( j
63.         for(j=widthS+1;j<((widthS-1)*(heightS-1))-2;j++)
    
64.         {
    4 h" L% c3 O- ~9 ^
65.                         ImageNow[j]=ImageOperation[j];
    
66.         }

复制代码

代码中有判断当前像素梯度方向落在哪个象限的代码，具体看代码。
效果图：


) g/ k: i/ ]4 S0 l4 I2 K8 d% o

9 ~) Q) g& U+ y/ x4.双阈值处理
6 M5 q# X! I9 n2 b  f' i在施加非极大值抑制之后，剩余的像素可以更准确地表示图像中的实际边缘。然而，仍然存在由于噪声和颜色变化引起的一些边缘像素。为了解决这些杂散响应，必须用弱梯度值过滤边缘像素，并保留具有高梯度值的边缘像素，可以通过选择高低阈值来实现。如果边缘像素的梯度值高于高阈值，则将其标记为强边缘像素；如果边缘像素的梯度值小于高阈值并且大于低阈值，则将其标记为弱边缘像素；如果边缘像素的梯度值小于低阈值，则会被抑制。阈值的选择取决于给定输入图像的内容。
" ]" w( W! b# ~具体实现代码如下：
5 f' |( O- T* F
% o* C& ^) P* L7 A, Z

1. 
   
2. nhist=(int*)malloc((400)*sizeof(int));
   
3.         
   
4.         for(i=0;i<400;i++)
   + }  h( T4 x' D; K2 q, H$ w  Q
5.                 nhist=0;
   
6.         for(i=0;i<widthS*heightS;i++)
   ; s& L( |& [9 w0 e' g" t" G
7.                 if(ImageNow==0xff)
   9 }- h' k6 x- b0 N
8.                         nhist[grey]++;
   ' N9 ^1 e! O" y+ t+ S% T+ l4 ^+ K
9.         sum=0;
   / o, o# N6 m  |" z$ w. e  s; _) J& Y
10.         for(i=0;i<400;i++)
    / w5 y7 W" k3 `( D3 i
11.                 sum+=nhist;
    
12.         sum2=0;
    
13.         for(i=0;i<400;i++)
    & Z: i- _! P/ i) s6 @( Y
14.                 {
    
15.                         sum2+=nhist;
    
16.                         if(sum2>(0.7*sum+0.5))
    
17.                         {
    
18.                                 high=i;                                                        //取出高阈值
    
19.                                 break;
    
20.                         }
    
21.                 }
    
22.         low=(int)(0.4*((float)high)+0.5);//取出低阈值
    % G3 C4 D0 D' V3 R: {8 \% E5 |
23.         free(nhist);
    
24.         
    4 j' [, d1 I, J4 x5 i
25.                 for(i=widthS+1;i<widthS*heightS;i++)双阈值检测处理
    
26.                 {
    5 X" V  f( H8 s: h/ U- ^
27.                         if((grey>=high)&&(ImageNow==0xff))
    
28.                                                         ImageNow=0xff;
    5 C+ p) o( n0 r6 k; d
29.                         else if(grey>low&&grey<0xff)
    
30.                         {
    : [( F0 s/ \  d2 y1 ~* e2 q
31.                                 judge=1;
    
32.                                 while(judge)
    
33.                                 {
    
34.                                         for(j=-widthS;j<widthS+1;j+=widthS)
      t5 \) U/ v3 F/ P
35.                                         {
    + J1 s6 t2 x9 p
36.                                                 for(k=-1;k<2;k++)
    2 I: @+ A7 g, @: m; g! H
37.                                                 {
    
38.                                                         if((grey[i+j+k-1]>low)&&(ImageNow==0xff))
    
39.                                                         {
    % [4 W( n+ R0 ?+ F" l! @" T2 V  o
40.                                                                 ImageNow[i+j+k-1]=0xff;
    3 m7 c9 ~, V$ N0 K1 B
41.                                                                 i=i+j+k;
    
42.                                                                 cnt=0;
    
43.                                                                 break;
    ! Z: M. P. X6 D+ ~) m. L
44.                                                         }
    
45.                                                         else
    
46.                                                                 cnt++;
    
47.                                                         if(cnt>=8)
    
48.                                                         {
    ( q6 J! j% f, X" Y9 h
49.                                                                 cnt=0;
    
50.                                                                 judge=0;
    
51.                                                                 break;
    
52.                                                         }
    
53.                                                 }
    
54.                                                 if(judge==0||cnt==0)
    & m" |1 T5 L8 d/ p2 h0 J( H
55.                                                         break;
    
56.                                         }
    9 p( y6 e0 m' g- [4 d0 Q* z
57.                                         if(j>widthS+1||cnt==0)
    
58.                                                 judge=0;
    ) h1 _0 I$ M- m
59.                                 }
    2 J% l: f. p# g( p2 g9 Y
60.                         }
    
61.                         else
    
62. <i>                                              </i>          ImageNow<span style="font-style: italic;"><span style="font-style: normal;">=0;</span>
    % J9 Q" F  t3 i$ B
63.                </span> }
    4 z8 g6 X7 k$ P, e1 k$ [/ I1 |6 [/ s: z

复制代码

处理效果图：

/ z' L: Q3 [6 W% n* Y3 s2 L: ~% m

1 H$ O* J5 I" ?% z
由于我的图片噪声较小，与上一步的处理结果看不出太大的差别。
' J  B3 n9 q3 W/ d; H) Q) R. J& v
stm32h7的Canny边缘检测算法基本完毕，在32上运行起来和Matlab上的差距很大，我觉得应该是我图片存储方式的问题，我这里是用ImageLcd软件将图片直接转为一维数组进行操作显示的。建议使用SD卡存储图片在读取运算（由于我手头没有SD卡只能用这种方式）。
9 V+ H8 R; j" X9 N0 d  j4 T7 T
0 T0 R- S0 A& d9 X6 q

& x8 A  p0 o! D