算法（十一）：imglib之sobel

Sobel是一个边缘检测的算子，也可以对8位图像进行边缘检测。

由于imglib提供了好几个sobel的算法，这里我们使用的是IMG_sobel_3x3_8()

IMG_sobel_3x3_8函数并不是用C语言写的，而是用DSP的汇编语言写的，执行效率还是很高的。

1. 
   
2.                 .if     __TI_EABI__
   
3.                     .asg    IMG_sad_8x8,  _IMG_sad_8x8
   
4.                 .endif
   
5. 
   
6.                 .sect ".text:optimized"
   
7.                 .global _IMG_sad_8x8
   
8. 
   
9. _IMG_sad_8x8:  
   
10. * ===================== SYMBOLIC REGISTER ASSIGNMENTS ===================== *
    
11.         .asg            A0,         A_p
    
12.         .asg            A1,         A_i
    
13.         .asg            A2,         A_s0
    
14.         .asg            A3,         A_k1
    
15.         .asg            A4,         A_srcImg
    
16.         .asg            A4,         A_retval
    
17.         .asg            A5,         A_sad
    
18.         .asg            A6,         A_pitch
    
19.         .asg            A6,         A_s3210
    
20.         .asg            A7,         A_d7654
    
21.         .asg            A7,         A_s7654
    
22.         .asg            A8,         A_r3210
    
23.         .asg            A8,         A_s1
    
24.         .asg            A9,         A_d3210
    
25.         .asg            A9,         A_r7654
    
26.         .asg            B3,         B_retaddr
    
27.         .asg            B4,         B_refImg
    
28.         .asg            B5,         B_pitch
    
29.         .asg            B6,         B_r3210
    
30.         .asg            B7,         B_d7654
    
31.         .asg            B7,         B_r7654
    
32.         .asg            B8,         B_d3210
    
33.         .asg            B8,         B_s3210
    
34.         .asg            B9,         B_s7654
    
35.         .asg            B16,        B_sad
    
36.         .asg            B17,        B_srcImg
    
37.         .asg            B18,        B_s2
    
38.         .asg            B19,        B_s3
    
39. * ========================================================================= *
    
40. * =========================== PIPE LOOP PROLOG ============================ *
    
41.         LDNDW   .D2T2   *B_refImg,              B_r7654:B_r3210 ;[ 1,1]
    
42. ||      B       .S2     loop
    
43. ||      MVK     .S1     0x0101,     A_k1
    
44. ||      ADD     .L2X    A_srcImg,   8,          B_srcImg
    
45. 
    
46.         LDDW    .D2T1   *B_srcImg++(16),        A_s7654:A_s3210 ;[ 2,1]
    
47. ||      LDDW    .D1T2   *A_srcImg++(16),        B_s7654:B_s3210 ;[ 2,1]
    
48. ||      MV      .L2X    A_pitch,    B_pitch
    
49. ||      PACK2   .L1     A_k1,       A_k1,       A_k1
    
50. 
    
51.         LDNDW   .D2T1   *++B_refImg(B_pitch),   A_r7654:A_r3210 ;[ 3,1]
    
52. ||      ZERO    .L2     B_sad
    
53. ||      ZERO    .L1     A_sad
    
54. ||      MVK     .S1     3,          A_i
    
55. 
    
56. * =========================== PIPE LOOP KERNEL ============================ *
    
57. loop:
    
58.   [！A_i]ADD     .L1     A_sad,      A_s1,       A_sad           ;[13,1]
    
59. ||[！A_i]ADD     .S2     B_sad,      B_s3,       B_sad           ;[13,1]
    
60. ||[ A_i]B       .S1     loop                                    ;[10,2]
    
61. ||      DOTPU4  .M1     A_d3210,    A_k1,       A_s0            ;[10,2]
    
62. ||      SUBABS4 .L2     B_s7654,    B_r7654,    B_d7654         ;[ 7,3]
    
63. ||      LDNDW   .D2T2   *++B_refImg(B_pitch),   B_r7654:B_r3210 ;[ 1,5]
    
64. 
    
65.   [！A_i]ADD     .S2     B_sad,      B_s2,       B_sad           ;[14,1]
    
66. ||      DOTPU4  .M2X    B_d7654,    A_k1,       B_s3            ;[ 8,3]
    
67. ||      SUBABS4 .L1     A_s7654,    A_r7654,    A_d7654         ;[ 8,3]
    
68. ||      SUBABS4 .L2     B_s3210,    B_r3210,    B_d3210         ;[ 8,3]
    
69. ||      LDDW    .D2T1   *B_srcImg++(16),        A_s7654:A_s3210 ;[ 2,5]
    
70. ||      LDDW    .D1T2   *A_srcImg++(16),        B_s7654:B_s3210 ;[ 2,5]
    
71. 
    
72.   [ A_i]SUB     .D1     A_i,        1,          A_i             ;[15,1]
    
73. ||[！A_i]ADD     .S1     A_sad,      A_s0,       A_sad           ;[15,1]
    
74. ||      DOTPU4  .M1     A_d7654,    A_k1,       A_s1            ;[ 9,3]
    
75. ||      DOTPU4  .M2X    B_d3210,    A_k1,       B_s2            ;[ 9,3]
    
76. ||      SUBABS4 .L1     A_s3210,    A_r3210,    A_d3210         ;[ 9,3]
    
77. ||      LDNDW   .D2T1   *++B_refImg(B_pitch),   A_r7654:A_r3210 ;[ 3,5]
    
78. 
    
79. * =========================== PIPE LOOP EPILOG ============================ *
    
80.         ADD     .L1     A_sad,      A_s1,       A_sad           ;[13,4]
    
81. ||      ADD     .S2     B_sad,      B_s3,       B_sad           ;[13,4]
    
82. ||      DOTPU4  .M1     A_d3210,    A_k1,       A_s0            ;[10,5]
    
83. 
    
84.         ADD     .L2     B_sad,      B_s2,       B_sad           ;[14,4]
    
85. ||      RET     .S2     B_retaddr
    
86. 
    
87.         ADD     .S1     A_sad,      A_s0,       A_sad           ;[15,4]
    
88. 
    
89.         ADD     .L1     A_sad,      A_s1,       A_sad           ;[13,5]
    
90. ||      ADD     .S2     B_sad,      B_s3,       B_sad           ;[13,5]
    
91. 
    
92.         ADD     .S2     B_sad,      B_s2,       B_sad           ;[14,5]
    
93. 
    
94.         ADD     .S1     A_sad,      A_s0,       A_sad           ;[15,5]
    
95. 
    
96.         ADD     .S1X    A_sad,      B_sad,      A_retval
    

函数并不需要用户指定mask，而是默认使用以下mask;

IMG_sobel_3x3_8()函数有四个参考，说明如下：

In：bmp图像的数据区指针

Out：运算后的结果

Cols：图像的宽度

Rows：图像的高度

所以对于一个576*448的图像，这个函数的用法就是：

IMG_sobel_3x3_8(In,out,bmpWidth,bmpHeight);

将out指向bmp文件的数据区，重新生成bmp文件后就可以看到sobel边缘检测的效果了。

以下是对8位灰度图像进行sobel边缘检测的结果

以下是对8位二值化图像进行sobel边缘检测的结果

最后和源码与二进制文件会在最后的结项中。