智能车摄像头开源—1.1核心代码:自顺应八向迷宫(上) ...

打印 上一主题 下一主题

主题 994|帖子 994|积分 2982

目录
一、开源项目简介
二、前期探索与算法确定
1.常用二值化算法对比(个人经验)
        大津法
        卷积核
        差比和
        自顺应迷宫(八领域)
三、自顺应八向迷宫


一、开源项目简介

        2024.8.19,决战前夕,正式决定将十个多月来的心血进行开源。本在区赛打完之后,就筹划如果国一,就部分开源,若是国二,就全部开源。但显然,在封车的那一刻,结局就已经注定,已经无望国一。以是此行怀揣失意而来,注定满载遗憾而去,但不甘心于自己全部的成果随着战败的我而隐去,因此在此将代码全部开源。
        筹划将整个摄像头代码进行详细讲述,并在后续连续更新,做成一套完备的系列,供新入门智能车图像处理的学弟学妹提供参考。同时也谨慎声明,此系列成果与经验均由本人实践得出,并不保证完全精确,仅供入门学习与参考。
        本人将智能车摄像头图像处理大致分为四部分:图像基础处理与边线提取(底层),信息提取(中层),元素判断与元素处理(上层),以及速度决议与控制相干代码(链接层)。
        本在考虑到底从那边讲起,但显然,图像的预处理与边线的提取为整个摄像头图像处理的基石,良好的边线提取算法能够快速、准确且高效的完成使命,同时对强光以及明暗不均的情况有较好的顺应性,不占用大量的运算资源。因此,第一篇就直接开门见山,将最关键的先行给出,后续再对其他处理进行补充阐明。
二、前期探索与算法确定

1.常用二值化算法对比(个人经验)

        大津法

        想必大多数摄像头位入门后接触的第一个图像算法便是大津法,大津法简单易用,同时网上流传的开源代码大多是大津法,因此可以快速的移植并直接使用。(若还没接触过大津法的铁柱或金发妹可以先行去了解原理与算法的实现,这里不过多解说原理,究竟有非常多的开源项目可供参考)。但使用大津二值化图像,淹灭时间较长,同时使用单阈值处理整幅图像,使其对局部较亮或较暗的情况表现并不抱负。
        针对于对光线顺应性不佳的题目,网上又给出了分块大津含糊大津等优化算法,在一定水平上进步了对光线的顺应性,也优化了计算时间,但仍无法到达我满足的结果,且计算时长对主控的占用较长,使一些运算性能较低的主控压力拉满。若不想在此处淹灭时间,大津法还是比较推荐的,可以直接移植现成的代码后便可直接使用,得到二值化后的图像,但显然,这并不是我心目中的抱负算法。
        这里给出常规大津和三分割大津法:(算法泉源于学长传承)
  1. /************************************************************************
  2. *  函数名称:short GetOSTU (unsigned char tmImage[LCDH][LCDW])
  3. *  功能说明:大津法求阈值大小
  4. *  参数说明:tmImage : 图像数据
  5. *  函数返回:无
  6. *  修改时间:2011年10月28日
  7. *  备    注:  GetOSTU(Image_Use);//大津法阈值
  8. Ostu方法又名最大类间差方法,通过统计整个图像的直方图特性来实现全局阈值T的自动选取,其算法步骤为:
  9. 1) 先计算图像的直方图,即将图像所有的像素点按照0~255共256个bin,统计落在每个bin的像素点数量
  10. 2) 归一化直方图,也即将每个bin中像素点数量除以总的像素点
  11. 3) i表示分类的阈值,也即一个灰度级,从0开始迭代 1
  12. 4) 通过归一化的直方图,统计0~i 灰度级的像素(假设像素值在此范围的像素叫做前景像素) 所占整幅图像
  13. 的比例w0,        并统计前景像素的平均灰度u0;统计i~255灰度级的像素(假设像素值在此范围的像素叫做背
  14. 景像素)  * 所占整幅图像的比例w1,并统计背景像素的平均灰度u1;
  15. 5) 计算前景像素和背景像素的方差 g = w0*w1*(u0-u1) (u0-u1)
  16. 6) i++;转到4),直到i为256时结束迭代
  17. 7) 将最大g相应的i值作为图像的全局阈值
  18. 缺陷:OSTU算法在处理光照不均匀的图像的时候,效果会明显不好,因为利用的是全局像素信息。
  19. ************************************************************************/
  20. short GetOSTU (unsigned char tmImage[LCDH][LCDW])
  21. {
  22.   signed short i, j;
  23.   unsigned long Amount = 0;
  24.   unsigned long PixelBack = 0;
  25.   unsigned long PixelshortegralBack = 0;
  26.   unsigned long Pixelshortegral = 0;
  27.   signed long PixelshortegralFore = 0;
  28.   signed long PixelFore = 0;
  29.   float OmegaBack, OmegaFore, MicroBack, MicroFore, SigmaB, Sigma; // 类间方差;
  30.   signed short MinValue, MaxValue;
  31.   signed short Threshold = 0;
  32.   unsigned char HistoGram[256] = {0}; //初始化灰度直方图
  33.   for (j = 0; j < LCDH; j++)
  34.   {
  35.     for (i = 0; i < LCDW; i++)
  36.     {
  37.       HistoGram[tmImage[j][i]]++; //统计灰度级中每个像素在整幅图像中的个数
  38.     }
  39.   }
  40.   for (MinValue = 0; MinValue < 256 && HistoGram[MinValue] == 0; MinValue++);        //获取最小灰度的值
  41.   for (MaxValue = 255; MaxValue > MinValue && HistoGram[MinValue] == 0; MaxValue--); //获取最大灰度的值
  42.   if (MaxValue == MinValue)
  43.     return MaxValue;         // 图像中只有一个颜色
  44.   if (MinValue + 1 == MaxValue)
  45.     return MinValue;        // 图像中只有二个颜色
  46.   for (j = MinValue; j <= MaxValue; j++)
  47.     Amount += HistoGram[j];        //  像素总数
  48.   Pixelshortegral = 0;
  49.   for (j = MinValue; j <= MaxValue; j++)
  50.   {
  51.     Pixelshortegral += HistoGram[j] * j;        //灰度值总数
  52.   }
  53.   SigmaB = -1;
  54.   for (j = MinValue; j < MaxValue; j++)
  55.   {
  56.     PixelBack = PixelBack + HistoGram[j];     //前景像素点数
  57.     PixelFore = Amount - PixelBack;           //背景像素点数
  58.     OmegaBack = (float) PixelBack / Amount;   //前景像素百分比
  59.     OmegaFore = (float) PixelFore / Amount;   //背景像素百分比
  60.     PixelshortegralBack += HistoGram[j] * j;  //前景灰度值
  61.     PixelshortegralFore = Pixelshortegral - PixelshortegralBack;  //背景灰度值
  62.     MicroBack = (float) PixelshortegralBack / PixelBack;   //前景灰度百分比
  63.     MicroFore = (float) PixelshortegralFore / PixelFore;   //背景灰度百分比
  64.     Sigma = OmegaBack * OmegaFore * (MicroBack - MicroFore) * (MicroBack - MicroFore);   //计算类间方差
  65.     if (Sigma > SigmaB)                    //遍历最大的类间方差g //找出最大类间方差以及对应的阈值
  66.     {
  67.       SigmaB = Sigma;
  68.       Threshold = j;
  69.     }
  70.   }
  71. //  Threshold += 40;
  72.   return Threshold;                        //返回最佳阈值;
  73. }
  74. /************************************************************************
  75. *  函数名称:void Get_Bin_Image_seg (unsigned char mode)
  76. *  功能说明:分割大津法,将图片水平分成三部分,分别通过大津法取阈值,判断每部分的阈值与全局阈值的差别,如果相差太大 ,使用全局阈值。
  77. *  参数说明:
  78. *  函数返回:无
  79. *  修改时间:2022年7月28日
  80. *  备    注:  Get_Bin_Image(0); //使用大津法二值化
  81. ************************************************************************/
  82. void Get_Bin_Image_seg()
  83. {
  84.     uint8 Thresholds[3]={0};    //存储三部分阈值
  85.     uint8 threshold1=0;
  86.     uint8 (*p_image)[94] = &OLED_Bin_data[0];
  87.     //  计算全局阈值,在下面分成三个部分,分别计算各部分的阈值,如果和全局阈值偏差过大,则使用全局阈值
  88.     threshold1= SimplifyGetOSTU(OLED_Bin_data[0],LCDW,LCDH);
  89.     for(k=0; k<3; k++)
  90.     {
  91.         maxvar = 0;
  92.         w0 = 0;
  93.         u = 0;
  94.         gray_hh = 0;
  95.         u0tmp = 0;
  96.         var = 0;
  97.         Thresholds[k] = 0;
  98.         // 初始化灰度比例直方图
  99.          for (i = 0; i < 256; i++)
  100.          {
  101.             bin_float[i] = 0;
  102.          }
  103.          uint32 gray_sum=0;
  104.          for (i = nchu_point[k].y0; i <= nchu_point[k].y1; i++)
  105.          {
  106.              for (j = nchu_point[k].x0; j <= nchu_point[k].x1; j++)
  107.              {
  108.                      ++bin_float[*(*(p_image + i) + j)];
  109.                      gray_sum+=*(*(p_image + i) + j);
  110.              }
  111.          }
  112.          size = (nchu_point[k].y1 - nchu_point[k].y0 + 1) * (nchu_point[k].x1 - nchu_point[k].x0 + 1);
  113.          for (i = 0; i < 256; i++)
  114.          {
  115.              bin_float[i] = bin_float[i] / size; // 每个灰度值的像素点所占比例
  116.              u += i * bin_float[i];
  117.          }
  118.          //创建比例灰度直方图
  119.          for (i = 0; i < 256; i++)
  120.          {
  121.              /**
  122.              w0 += bin_float[i];
  123.              gray_hh += i * bin_float[i];             //灰度和
  124.              u0 = gray_hh / w0;
  125.              var = (u0 - u) * (u0 - u) * w0 / (1 - w0);
  126.              if (var > maxvar)
  127.              {
  128.                      maxgray = gray_hh;
  129.                      maxbin = w0;
  130.                      maxvar = var;
  131.                      Thresholds[k] = (uint8)i;
  132.              }**/
  133.              w0 += bin_float[i];//背景部分每个灰度值的像素点所占比例之和
  134.              u0tmp += i * bin_float[i];             //灰度和
  135.              w1 = 1-w0;
  136.              u1tmp = gray_sum/size - u0tmp;
  137.              u0 = u0tmp / w0;
  138.              u1 = u1tmp / w1;
  139.              u = u0tmp+ u1tmp;
  140.              var = (u0 - u) * (u0 - u) * w0 / (1 - w0);
  141.              if (var > maxvar)
  142.              {
  143.                      maxgray = u0tmp;
  144.                      maxbin = w0;
  145.                      maxvar = var;
  146.                      Thresholds[k] = (uint8)i;
  147.              }
  148.          }
  149.          if (k == 0)
  150.          {
  151.              if (gray_hh > 15 && gray_hh <= 33)
  152.              {
  153.                  if (maxbin < 0.9f)
  154.                  {
  155.                          Thresholds[k] = (uint8)(Thresholds[1] - 3);
  156.                  }
  157.              }
  158.              else if (gray_hh > 41 && gray_hh <= 47)
  159.              {
  160.                  if (maxbin < 0.64f || maxbin > 0.76f)
  161.                  {
  162.                          Thresholds[k] = (uint8)(Thresholds[1] - 3);
  163.                  }
  164.              }
  165.              else if (gray_hh > 50 && gray_hh <= 60)
  166.              {
  167.                  if (maxbin < 0.42f || maxbin > 0.58f)
  168.                  {
  169.                          Thresholds[k] = (uint8)(Thresholds[1] - 3);
  170.                  }
  171.              }
  172.              if (fabs(threshold1 - Thresholds[k]) >= 30)//相差太大直接用标准值
  173.              {
  174.                  Thresholds[k] = threshold1;
  175.              }
  176.              else
  177.              {
  178.                  Thresholds[k] = (uint8)(Thresholds[k] + 0.5f * (threshold1 - Thresholds[k]));//这里是取了平均,可以通过调整占比适应不同的光照条件
  179.              }
  180.          }
  181.          else if (k == 1)
  182.          {
  183.              if(gray_hh > 69 && gray_hh < 80)
  184.              {
  185.                  if (maxbin > 0.15f)
  186.                  {
  187.                      Thresholds[k] = (uint8)(Thresholds[0] + 3);
  188.                  }
  189.              }
  190.              if (fabs(threshold1 - Thresholds[k]) >= 30)//相差太大直接用标准值
  191.              {
  192.                  Thresholds[k] = threshold1;
  193.              }
  194.              else
  195.              {
  196.                  Thresholds[k] = (uint8)(Thresholds[k] + 0.5f * (threshold1 - Thresholds[k]));//这里是取了平均,可以通过调整占比适应不同的光照条件
  197.              }
  198.          }
  199.          else if (k == 2)
  200.          {
  201.              if (maxbin < 0.85f && gray_hh < 28)
  202.              {
  203.                  Thresholds[k] = (uint8)(Thresholds[1] - 3);
  204.              }
  205.              else if(gray_hh > 69 && gray_hh < 79)
  206.              {
  207.                  if (maxbin < 0.5f || maxbin > 0.15f)
  208.                  {
  209.                      Thresholds[k] = (uint8)(Thresholds[1] - 3);
  210.                  }
  211.              }
  212.               if (fabs(threshold1 - Thresholds[k]) >= 30)//相差太大直接用标准值
  213.              {
  214.                  Thresholds[k] = threshold1;
  215.              }
  216.              else
  217.              {
  218.                  Thresholds[k] = (uint8)(Thresholds[k] + 0.5f * (threshold1 - Thresholds[k]));//这里是取了平均,可以通过调整占比适应不同的光照条件
  219.              }
  220.          }
  221.          for (i = 0; i < Thresholds[k]; i++)
  222.          {
  223.              Bin_Array[i] = 0;
  224.          }
  225.          for (i = Thresholds[k]; i < 256; i++)
  226.          {
  227.              Bin_Array[i] = 255;
  228.          }
  229.          for (i = nchu_point[k].y0; i <= nchu_point[k].y1; i++)
  230.          {
  231.              for (j = nchu_point[k].x0; j <= nchu_point[k].x1; j++)
  232.              {
  233.                  Pixle[i][j] = Bin_Array[*(*(p_image + i) + j)];
  234.              }
  235.          }
  236.     }
  237. }
  238. /************************************************************************
  239. *  函数名称:void Get_Bin_Image (unsigned char mode)
  240. *  功能说明:图像二值化到Bin_Image[][]
  241. *  参数说明:mode  :
  242. *    0:使用大津法阈值
  243. *    1:使用平均阈值
  244. *    2: sobel 算子改进型  手动阈值,同时输出改为提取边沿的图像
  245. *    3:sobel 算子改进型   动态阈值,同时输出改为提取边沿的图像
  246. *  函数返回:无
  247. *  修改时间:2020年10月28日
  248. *  备    注:  Get_Bin_Image(0); //使用大津法二值化
  249. ************************************************************************/
  250. void Get_Bin_Image (unsigned char mode)
  251. {
  252.   unsigned short i = 0, j = 0;
  253.   unsigned short Threshold = 0;
  254.   unsigned long tv = 0;
  255.   if (mode == 0)
  256.   {
  257.     Threshold = GetOSTU(OLED_Bin_data);  //大津法阈值
  258.   }
  259.   if (mode == 1)
  260.   {
  261.     //累加
  262.     for (i = 0; i < LCDH; i++)
  263.     {
  264.       for (j = 0; j < LCDW; j++)
  265.       {
  266.         tv += OLED_Bin_data[i][j];   //累加
  267.       }
  268.     }
  269.     Threshold =(unsigned short)(tv / LCDH / LCDW);   //求平均值,光线越暗越小,全黑约35,对着屏幕约160,一般情况下大约100
  270.     Threshold = Threshold + 20;      //此处阈值设置,根据环境的光线来设定
  271.   }
  272.   // 二值化 //
  273.   for (i = 0; i < LCDH-1; i++)
  274.   {
  275.     for (j = 0; j < LCDW-1; j++)
  276.     {
  277.       if (OLED_Bin_data[i][j] > Threshold) //数值越大,显示的内容越多,较浅的图像也能显示出来
  278.         Pixle[i][j] = 1;//白
  279.       else
  280.         Pixle[i][j] = 0;//黑
  281.     }
  282.   }
  283. }
复制代码
        卷积核

        使用卷积核提取边沿进而二值化图像不失为一个抱负的选择。常用的卷积核包括Sobel卷积核,拉普拉斯卷积核等,龙邱在其官方例程中给出了基于Sobel的二值化算法,代码如下:
  1. void sobelAutoThreshold (u8 source[MT9V03X_DVP_H/2][MT9V03X_DVP_W],u8 target[MT9V03X_DVP_H/2][MT9V03X_DVP_W])
  2. {
  3.     /** 卷积核大小 */
  4.     short KERNEL_SIZE = 3;
  5.     short xStart = KERNEL_SIZE / 2;
  6.     short xEnd = MT9V03X_DVP_W - KERNEL_SIZE / 2;
  7.     short yStart = KERNEL_SIZE / 2;
  8.     short yEnd = MT9V03X_DVP_H - KERNEL_SIZE / 2;
  9.     short i, j, k;
  10.     short temp[3];
  11.     for (i = yStart; i < yEnd; i++)
  12.     {
  13.         for (j = xStart; j < xEnd; j++)
  14.         {
  15.             /* 计算不同方向梯度幅值  */
  16.             temp[0] = -(short) source[i - 1][j - 1] + (short) source[i - 1][j + 1]            // {-1, 0, 1},
  17.             - (short) source[i][j - 1] + (short) source[i][j + 1]                                          // {-1, 0, 1},
  18.             - (short) source[i + 1][j - 1] + (short) source[i + 1][j + 1];                       // {-1, 0, 1};
  19.             temp[1] = -(short) source[i - 1][j - 1] + (short) source[i + 1][j - 1]             // {-1, -1, -1},
  20.             - (short) source[i - 1][j] + (short) source[i + 1][j]                                  // { 0,  0,  0},
  21.             - (short) source[i - 1][j + 1] + (short) source[i + 1][j + 1];                           // { 1,  1,  1};
  22.             temp[2] = -(short) source[i - 1][j] + (short) source[i][j - 1]                                //  {0, -1, -1},
  23.             - (short) source[i][j + 1] + (short) source[i + 1][j]                                       //  {1,  0, -1},
  24.             - (short) source[i - 1][j + 1] + (short) source[i + 1][j - 1];                            //  {1,  1,  0};
  25.             
  26.             temp[3] = -(short) source[i - 1][j] + (short) source[i][j + 1]                        // {-1, -1,  0},
  27.             - (short) source[i][j - 1] + (short) source[i + 1][j]                                        // {-1,  0,  1},
  28.             - (short) source[i - 1][j - 1] + (short) source[i + 1][j + 1];                             //  {0,  1,  1};
  29.             temp[0] = abs(temp[0]);
  30.             temp[1] = abs(temp[1]);
  31.             temp[2] = abs(temp[2]);
  32.             temp[3] = abs(temp[3]);
  33.             /* 找出梯度幅值最大值  */
  34.             for (k = 1; k < 3; k++)
  35.             {
  36.                 if (temp[0] < temp[k])
  37.                 {
  38.                     temp[0] = temp[k];
  39.                 }
  40.             }
  41.             /* 使用像素点邻域内像素点之和的一定比例    作为阈值  */
  42.             temp[3] =
  43.                     (short) source[i - 1][j - 1] + (short) source[i - 1][j] + (short) source[i - 1][j + 1]
  44.                     + (short) source[i][j - 1] + (short) source[i][j] + (short) source[i][j + 1]
  45.                     + (short) source[i + 1][j - 1] + (short) source[i + 1][j] + (short) source[i + 1][j + 1];
  46.             if (temp[0] > (temp[3] / 12.0f))
  47.             {
  48.                 target[i][j] = 0xFF;
  49.             }
  50.             else
  51.             {
  52.                 target[i][j] = 0x00;
  53.             }
  54.         }
  55.     }
  56. }
  57. ————————————————
  58.                             版权声明:本文为博主原创文章,遵循 CC 4.0 BY-SA 版权协议,转载请附上原文出处链接和本声明。
  59.                         
  60. 原文链接:https://blog.csdn.net/fgfgfdg8/article/details/126403306
复制代码

        卷积核在图像处理中有着非常普遍的使用。其中Sobel重要用于边沿检测,其原理是使用卷积核(即一个3*3的数字矩阵)计算每一个中心点灰度值变革最大的梯度方向。例如一个3*3区域内,第一行三个点的灰度值分别为:200, 200, 200,而最下面一行的三个点灰度值分别为:50, 50, 50,那么可如许计算:ABS((a1*200 + a2*200 + a3*200)-  (b1*50 + b2*50 + b3*50)),得出的结果与一个阈值进行比较(a1,a2,a3,b1,b2,b3为系数)。在这种情况下,灰度值由上到下有极其剧烈的变革,那么这个地方一个较大的梯度方向便是由上往下,即可以判断存在边沿(举例子,实际上并不会出现这么好的数值)。同样的,卷积核一般计算四个方向,分别为左右、上下和两个对角的梯度变革。这里简单解说原理,可以搜刮卷积核的资料进行详细学习。
        回到龙邱代码中,其函数输入第一参数为要处理的图像,第二参数为用于存储处理后的图像。改变最后一步 if 语句中temp[3]所除的系数,即可改变对边线提取的精致水平(即改变梯度对比的阈值,梯度值大于这个阈值即判断存在边界,越精致会导致存在的噪点更多,需要后续更多的降噪处理,不宜灵敏度过大,否则噪点会严重影响边线提取),而改变最后一步中 if 语句中的“ > ”为 “ < ”可将黑白反向,当读懂代码后,黑白反向原理就很简单了,不需过多赘述。
        卷积核自己具备对光线较好的顺应性,但显然运算时间也较长,与大津相差不大,相比于大津,若不想在边线提取深究,建议选用卷积核将图像二值化,个人以为在一定水平上结果要优于大津法。
        其中也曾偶尔得到优化后的Sobel算法,但并未进行验证,这里直接给出,结果自行查验:
  1. //imag[cam_w][cam_h]原图像数组
  2. //imag1[cam_w][cam_h]处理后图像数组
  3. //
  4. //180mhz 主频lpcdan片机处理60 * 120图像实测7.2ms
  5. //如果需要更快的处理速度可选择局部处理提高速度
  6. * /
  7. //#define cam_h  60   //高度(行)
  8. //#define cam_w  120  //宽度(列)
  9. void sobel(uint8_t imag[cam_h][cam_w], uint8_t imag1[cam_h][cam_w])
  10. {
  11.     int tempx = 0, tempy = 0, temp = 0, i = 0, j = 0;
  12.     for (i = 1; i < cam_h - 1; i++)
  13.     {
  14.         for (j = 1; j < cam_w - 1; j++)
  15.         {
  16.             tempx = (-imag[i - 1][j - 1]) + (-2 * imag[i][j - 1]) + (-imag[i + 1][j - 1])
  17.                   + (imag[i - 1][j + 1])  + (2 * imag[i][j + 1])  + (imag[i + 1][j + 1]);
  18.             if (tempx < 0)
  19.             {
  20.                 tempx = -tempx;
  21.             }
  22.             tempy = (imag[i + 1][j - 1])  + (2 * imag[i + 1][j])  + (imag[i + 1][j + 1])
  23.                   + (-imag[i - 1][j - 1]) + (-2 * imag[i - 1][j]) + (-imag[i - 1][j + 1]);
  24.             if (tempy < 0)
  25.             {
  26.                 tempy = -tempy;
  27.             }
  28.             temp = tempx + tempy;
  29.             if (temp > 255)
  30.             {
  31.                 temp = 255;
  32.             }
  33.             imag1[i][j] = temp;
  34.         }
  35.     }
  36. }
  37. temp = sqrt(tempx * tempx + tempy * tempy);
  38. //这里加点我的理解:
  39. //1.可将上述函数内:temp = tempx + tempy 替换为 temp = sqrt(tempx * tempx + tempy * tempy) 以优化
  40. //2.算法整体好像并没有对图像二值化,只是将梯度结果限幅后丢给图像。至于二值化可以参考上述龙邱的方法
  41. //3.龙邱方法是比对四个方向内的梯度最大值后与阈值比较,但此处算法不同,需读者仔细斟酌,至于怎么处理就得依靠读者发挥了
复制代码
        差比和

        不久前逐飞官方的摄像头培训中给出了差比和求取边线的方法,原理为:从最长白直列开始,向左(向右)依次寻找边线点,对于左侧找边线,用当前点的灰度值减去右侧相距五个点的灰度值,再将差值除二者的和值并取绝对值,再将结果向左移7位(乘128倍,使用移位比直接乘更节省算力和时间,故不使用乘100),将最后结果与一个阈值进行比较,当大于阈值时,阐明当前点灰度值变革剧烈,即可判断为赛道黑边界。右侧边线同理,只是对比的点位于当前点的左侧。
        此算法对光线变革表现出较好的顺应性能,同时在我们双核200Mhz的主控上,单核内处理十字一张188 * 120图像到达了0.6ms的运算速度,极其之快,显然是一个非常不错的抱负算法,可以快速地提取出赛道的一维边线数组。代码这里不列出,可自行在Bilibili中搜刮逐飞科技找到代码实现方法并进行后续优化。(注:逐飞并没有将算法进行顺应性优化,且给出的对比阈值固定,任何好的算法终极都需要非常多的优化才能表现出较好的顺应性和结果,想使用逐飞方案,显然需要自行对算法进行进一步处理)。至于我为什么没给出来,因为懒的写,而且原理很简单,如果真想学真想用,大可不必纠结没有完善的开源代码,最好的代码就在你手里。
        有些开源方案中,对a和b差值取绝对值,但我们只需要由白赛道对黑线的对比,不需由黑到白的对比,显然不取绝对值可有用防止误判。
        在本人算法中也有使用差比和的代码部分,这里给出源码:
  1. uint8 Compare_Value = 20;
  2. //差比和   ****
  3. //第一形参a为灰度值较大的点,第二形参为灰度值较低的点
  4. int16 Compare_Num(int16 a, int16 b, int8 compare_value)               //****
  5. {
  6.     if((((a - b) << 7) / (a + b)) > compare_value)
  7.     {
  8.         return 1;
  9.     }
  10.     else
  11.     {
  12.         return 0;
  13.     }
  14. }
复制代码
        自顺应迷宫(八领域)

          在算法的探索中,以上算法或是速度,或是对光线的顺应能力,都不能兼具,(除差比和方案外,但显然,差比和提取的信息量较少,每行左右边线各只有一个点,即一维边线,不利于后续的处理)以至于我始终不能确定终极的算法。但在网上寻找思路的过程中,接触到了上交AuTo战队开源出的自顺应迷宫算法,使我豁然开朗。自顺应是选用局部5×5区域内的像素灰度值均值作为当前区域分割黑白的阈值,而迷宫是一种爬线算法,可用八领域替代(固然,迷宫不输八邻域),上征战队将边线点的求取与爬线算法相结合。一方面,小局部阈值对明暗不均的图像表现出极好的顺应性,同时边爬线边求取边线点,避免了对赛道外以及赛道内无关图像的处理(只对边界进行提取),详细资料可以去网上找上交开源代码,细致研读后,定会受益匪浅。
        同时,b站上名为“__苏格拉没有底___”的开源博主也给出了自顺应八领域直接处理灰度图求出边线的源代码,原理与上交方案雷同,都是直接处理灰度图像,只是将爬线算法更换为八领域,但二者相差不大。也可自行搜刮获取。
       在明白透彻二者的方案后,我意识到这种直接处理灰度图的方式具备极高的上限,也赞叹于开源作者的能力与灵感。意识到这就是我一直渴求的边线提取算法,因此我毅然决然的选定这种算法作为我摄像头边线提取的核心算法。


三、自顺应八向迷宫

       在确定终极算法思路后,我又对代码进行更深的分析与明白,意识到仍然存在优化空间,因此将迷宫的速度上风与八邻域的方向上风相结合,给出自己的八向迷宫,对每一个点的阈值计算进行进一步优化,同时扩充爬线时获取的信息量,用于后续拐点判断和元素处理。经实测,在我们双核200MHz主控单核处理一张188 * 120的图像(包括压缩),到达了0.3ms左右的极速(其中还包罗一系列处理,会在后续给出),同时,算法对于差别光线场景均表现出很好的顺应性。对于本人来说,已是极为抱负的摄像头算法,固然,还需后来者评判。
       本人以为对于此算法,在运算速度方面的优化已经无限靠近于极致,固然,可能仍然存在优化空间等待后来者的探索,本人能力有限,已无法给出更多的优化思路。
       最快速的明白方式,便是对着代码的运行思路,绘图一步一步明白,自顺应八向迷宫对于新接触到此算法的车友有一些难度,但实在真正明白后会发现很简单,此算法的特点就是简单粗暴,以至于大道至简。
       因篇幅缘故原由,这里先将源代码给出,在下篇文章中再进行详细讲述:
  1. //**************方向计算****************
  2. //记录方向时,要将每个坐标变为单一数字,便于简化后续算法。例如八邻域的方向(逆时针    顺时针)
  3. //3  2  1               1  2  3
  4. //4     0               0     4
  5. //5  6  7               7  6  5
  6. /
  7. //{-1,-1},{0,-1},{+1,-1},
  8. //{-1, 0},       {+1, 0},
  9. //{-1,+1},{0,+1},{+1,+1},
  10. //迷宫缺点是无法像八邻域一样逐步记录方向,但八邻域无非就是以上八个点,且与中心点差值固定,不会随迷宫算法的朝向和移动方向而改变,因此先将每个横坐标乘 3
  11. //{-3,-1},{0,-1},{+3,-1},
  12. //{-3, 0},       {+3, 0},
  13. //{-3,+1},{0,+1},{+3,+1},
  14. //再将横坐标减去纵坐标
  15. // -2, 1, 4
  16. // -3,    3
  17. // -4,-1, 2
  18. //由此可以得到一个八邻方向坐标。只需在每次移动后,在方向数组里记录对应数字,就可确定生长方向
  19. //此算法无任何原理,只是为了得到八个不一样的值可以用来判定方向,横坐标可以乘大于 2 的任意值,2以内会出现重复
  20. //************************************
  21. //******************自适应方向迷宫参数************************
  22. const int8 L_Face_Dir[4][2] = {{0,-1},{1,0},{0,1},{-1,0}};  //左侧迷宫面向
  23. //  0
  24. //3   1
  25. //  2
  26. const int8 L_Face_Dir_L[4][2] = {{-1,-1},{1,-1},{1,1},{-1,1}};  //左侧面向的左前方
  27. //0   1
  28. //
  29. //3   2
  30. const int8 R_Face_Dir[4][2] = {{0,-1},{1,0},{0,1},{-1,0}};  //右侧迷宫面向
  31. //  0
  32. //3   1
  33. //  2
  34. const int8 R_Face_Dir_R[4][2] = {{1,-1},{1,1},{-1,1},{-1,-1}};  //右侧面向的右前方
  35. //3   0
  36. //
  37. //2   1
  38. const int8 Square_0[25][2] = {              //一个5 * 5的矩阵,用来求中心点周围的局部阈值(局部阈值)
  39. {-2,-2},{-1,-2},{0,-2},{+1,-2},{+2,-2},
  40. {-2,-1},{-1,-1},{0,-1},{+1,-1},{+2,-1},
  41. {-2,-0},{-1, 0},{0, 0},{+1, 0},{+2,-0},
  42. {-2,+1},{-1,+1},{0,+1},{+1,+1},{+2,+1},
  43. {-2,+2},{-1,+2},{0,+2},{+1,+2},{+2,+2}
  44. };
  45. //迷宫单侧停止爬线标志位
  46. uint8 L_Stop_Flag = 0;
  47. uint8 R_Stop_Flag = 0;
  48. //**********************************************************
  49. /******
  50. * 函数功能:      求取赛道二维数组边线
  51. * 特殊说明:      基于上交代码的自适应迷宫优化后的自适应八向迷宫
  52. * 形  参:        uint16 Break_Flag         最大循环次数,防止卡死程序,一般为3~4倍图像宽度
  53. *                uint8(*image)[Image_X]     提取边线的图像
  54. *                uint8(*l_line)[2]          存放左侧边线的二维数组
  55. *                uint8(*r_line)[2]          存放右侧边线的二维数组
  56. *                int8 *l_dir                存放左侧边线每个点的生长方向
  57. *                int8 *r_dir                存放右侧边线每个点的生长方向
  58. *                uint16 *l_stastic          记录左侧边线点的个数
  59. *                uint16 *r_stastic          记录右侧边线点的个数
  60. *                uint8 *x_meet              记录左右两侧爬线相遇点的X坐标
  61. *                uint8 *y_meet              记录左右两侧爬线相遇点的Y坐标
  62. *                uint8 l_start_x            左侧爬线起始点的X坐标
  63. *                uint8 l_start_y            左侧爬线起始点的Y坐标
  64. *                uint8 r_start_x            右侧爬线起始点的X坐标
  65. *                uint8 r_start_y            右侧爬线起始点的Y坐标
  66. *                uint8 clip_value           计算每个阈值时相加的经验值,一般为-5 ~ 5,避免强行分割,可直接设为0
  67. *
  68. * 示例:         Dir_Labyrinth_5((uint16)Use_Num, Find_Line_Image, Adaptive_L_Line, Adaptive_R_Line, Adaptive_L_Grow_Dir, Adaptive_R_Grow_Dir, &Adaptive_L_Statics, &Adaptive_R_Statics, &Adaptive_X_Meet, &Adaptive_Y_Meet,
  69.                    Adaptive_L_Start_Point[0], Adaptive_L_Start_Point[1], Adaptive_R_Start_Point[0], Adaptive_R_Start_Point[1], 0);
  70. * 返回值:        无
  71. */
  72. void Dir_Labyrinth_5(uint16 Break_Flag, uint8(*image)[Image_X], uint8(*l_line)[2], uint8(*r_line)[2], int8 *l_dir, int8 *r_dir, uint16 *l_stastic, uint16 *r_stastic, uint8 *x_meet, uint8 *y_meet,
  73.                      uint8 l_start_x, uint8 l_start_y, uint8 r_start_x, uint8 r_start_y, uint8 clip_value)
  74. {
  75.     uint8 j = 0;
  76.     L_Stop_Flag = 0;
  77.     R_Stop_Flag = 0;
  78. //左边变量
  79.     uint8  L_Center_Point[2] = {0};     //存放每次找到的XY坐标
  80.     uint16 L_Data_Statics = 0;          //统计左边找到的边线点的个数
  81.     uint8  L_Front_Value = 0;           //左侧 面向的前方点的灰度值
  82.     uint8  L_Front_L_Value = 0;         //左侧 面向的左前方点的灰度值
  83.     uint8  L_Dir = 0;                   //此参数用于转向
  84.     uint8  L_Turn_Num = 0;              //记录转向次数,若中心点前后左右都是黑色像素,就会在一个点转向四次,记录到四次时退出循环防止卡死程序
  85.     uint16 L_Pixel_Value_Sum = 0;       //中心点与周围24个点的像素值和
  86.     float L_Thres = 0;                 //局部阈值,即L_Pixel_Value_Sum / 25
  87. //右边变量
  88.     uint8  R_Center_Point[2] = {0};     //存放每次找到的XY坐标
  89.     uint16 R_Data_Statics = 0;          //统计右边找到的边线点的个数
  90.     uint8  R_Front_Value = 0;           //右侧 面向的前方点的灰度值
  91.     uint8  R_Front_R_Value = 0;         //右侧 面向的左前方点的灰度值
  92.     uint8  R_Dir = 0;                   //此参数用于转向
  93.     uint8  R_Turn_Num = 0;              //记录转向次数,若中心点前后左右都是黑色像素,就会在一个点转向四次,记录到四次时退出循环防止卡死程序
  94.     uint16 R_Pixel_Value_Sum = 0;       //中心点与周围24个点的像素值和
  95.     float R_Thres = 0;                 //局部阈值
  96. //第一次更新坐标点  将找到的起点值传进来
  97.     L_Center_Point[0] = l_start_x + 1;//x
  98.     L_Center_Point[1] = l_start_y;//y
  99.     R_Center_Point[0] = r_start_x - 1;//x
  100.     R_Center_Point[1] = r_start_y;//y
  101.     //开启方向迷宫循环
  102.     while (Break_Flag--)
  103.     {
  104.          //左边
  105.         //判定出死区后,挂出停止标志位,单侧爬线停止。
  106.         if(L_Stop_Flag == 0)
  107.         {
  108.             l_line[L_Data_Statics][0] = L_Center_Point[0];  //找到的中心点X坐标计入左边线数组
  109.             l_line[L_Data_Statics][1] = L_Center_Point[1];  //找到的中心点Y坐标计入左边线数组
  110.             if(L_Data_Statics != 0)
  111.             {
  112.                 switch(l_dir[L_Data_Statics - 1])  //下面这一坨可以根据上一个点的生长方向大幅优化爬线时间
  113.                 {
  114.                     //从第二个点开始,第一个点的阈值要25个点全部加一遍
  115.                     //当横向或纵向生长时,将原先的25次加法运算简化为十次计算
  116.                     //当斜向生长时,将原先的25次加法计算简化为十八次运算
  117.                     //可以画出图来,跟着代码走几遍,就可以很快速的理解
  118.                     case 1:
  119.                     {
  120.                         L_Pixel_Value_Sum = L_Pixel_Value_Sum - image[L_Center_Point[1] + 3][L_Center_Point[0] + 2] - image[L_Center_Point[1] + 3][L_Center_Point[0] + 1]
  121.                                                               - image[L_Center_Point[1] + 3][L_Center_Point[0] + 0] - image[L_Center_Point[1] + 3][L_Center_Point[0] - 1]
  122.                                                               - image[L_Center_Point[1] + 3][L_Center_Point[0] - 2]
  123.                                                               + image[L_Center_Point[1] - 2][L_Center_Point[0] + 2] + image[L_Center_Point[1] - 2][L_Center_Point[0] + 1]
  124.                                                               + image[L_Center_Point[1] - 2][L_Center_Point[0] + 0] + image[L_Center_Point[1] - 2][L_Center_Point[0] - 1]
  125.                                                               + image[L_Center_Point[1] - 2][L_Center_Point[0] - 2];
  126.                         break;
  127.                     }
  128.                     case -2:
  129.                     {
  130.                         L_Pixel_Value_Sum = L_Pixel_Value_Sum - image[L_Center_Point[1] - 1][L_Center_Point[0] + 3] - image[L_Center_Point[1] - 0][L_Center_Point[0] + 3]
  131.                                                               - image[L_Center_Point[1] + 1][L_Center_Point[0] + 3] - image[L_Center_Point[1] + 2][L_Center_Point[0] + 3]
  132.                                                               - image[L_Center_Point[1] + 3][L_Center_Point[0] + 3] - image[L_Center_Point[1] + 3][L_Center_Point[0] + 2]
  133.                                                               - image[L_Center_Point[1] + 3][L_Center_Point[0] + 1] - image[L_Center_Point[1] + 3][L_Center_Point[0] - 0]
  134.                                                               - image[L_Center_Point[1] + 3][L_Center_Point[0] - 1]
  135.                                                               + image[L_Center_Point[1] + 2][L_Center_Point[0] - 2] + image[L_Center_Point[1] + 1][L_Center_Point[0] - 2]
  136.                                                               + image[L_Center_Point[1] + 0][L_Center_Point[0] - 2] + image[L_Center_Point[1] - 1][L_Center_Point[0] - 2]
  137.                                                               + image[L_Center_Point[1] - 2][L_Center_Point[0] - 2] + image[L_Center_Point[1] - 2][L_Center_Point[0] - 1]
  138.                                                               + image[L_Center_Point[1] - 2][L_Center_Point[0] - 0] + image[L_Center_Point[1] - 2][L_Center_Point[0] + 1]
  139.                                                               + image[L_Center_Point[1] - 2][L_Center_Point[0] + 2];
  140.                         break;
  141.                     }
  142.                     case -3:
  143.                     {
  144.                         L_Pixel_Value_Sum = L_Pixel_Value_Sum - image[L_Center_Point[1] - 2][L_Center_Point[0] + 3] - image[L_Center_Point[1] - 1][L_Center_Point[0] + 3]
  145.                                                               - image[L_Center_Point[1] + 0][L_Center_Point[0] + 3] - image[L_Center_Point[1] + 1][L_Center_Point[0] + 3]
  146.                                                               - image[L_Center_Point[1] + 2][L_Center_Point[0] + 3]
  147.                                                               + image[L_Center_Point[1] - 2][L_Center_Point[0] - 2] + image[L_Center_Point[1] - 1][L_Center_Point[0] - 2]
  148.                                                               + image[L_Center_Point[1] - 0][L_Center_Point[0] - 2] + image[L_Center_Point[1] + 1][L_Center_Point[0] - 2]
  149.                                                               + image[L_Center_Point[1] + 2][L_Center_Point[0] - 2];
  150.                         break;
  151.                     }
  152.                     case -4:
  153.                     {
  154.                         L_Pixel_Value_Sum = L_Pixel_Value_Sum - image[L_Center_Point[1] - 3][L_Center_Point[0] - 1] - image[L_Center_Point[1] - 3][L_Center_Point[0] + 0]
  155.                                                               - image[L_Center_Point[1] - 3][L_Center_Point[0] + 1] - image[L_Center_Point[1] - 3][L_Center_Point[0] + 2]
  156.                                                               - image[L_Center_Point[1] - 3][L_Center_Point[0] + 3] - image[L_Center_Point[1] - 2][L_Center_Point[0] + 3]
  157.                                                               - image[L_Center_Point[1] - 1][L_Center_Point[0] + 3] - image[L_Center_Point[1] + 0][L_Center_Point[0] + 3]
  158.                                                               - image[L_Center_Point[1] + 1][L_Center_Point[0] + 3]
  159.                                                               + image[L_Center_Point[1] - 2][L_Center_Point[0] - 2] + image[L_Center_Point[1] - 1][L_Center_Point[0] - 2]
  160.                                                               + image[L_Center_Point[1] + 0][L_Center_Point[0] - 2] + image[L_Center_Point[1] + 1][L_Center_Point[0] - 2]
  161.                                                               + image[L_Center_Point[1] + 2][L_Center_Point[0] - 2] + image[L_Center_Point[1] + 2][L_Center_Point[0] - 1]
  162.                                                               + image[L_Center_Point[1] + 2][L_Center_Point[0] - 0] + image[L_Center_Point[1] + 2][L_Center_Point[0] + 1]
  163.                                                               + image[L_Center_Point[1] + 2][L_Center_Point[0] + 2];
  164.                         break;
  165.                     }
  166.                     case -1:
  167.                     {
  168.                         L_Pixel_Value_Sum = L_Pixel_Value_Sum - image[L_Center_Point[1] - 3][L_Center_Point[0] - 2] - image[L_Center_Point[1] - 3][L_Center_Point[0] - 1]
  169.                                                               - image[L_Center_Point[1] - 3][L_Center_Point[0] + 0] - image[L_Center_Point[1] - 3][L_Center_Point[0] + 1]
  170.                                                               - image[L_Center_Point[1] - 3][L_Center_Point[0] + 2]
  171.                                                               + image[L_Center_Point[1] + 2][L_Center_Point[0] - 2] + image[L_Center_Point[1] + 2][L_Center_Point[0] - 1]
  172.                                                               + image[L_Center_Point[1] + 2][L_Center_Point[0] + 0] + image[L_Center_Point[1] + 2][L_Center_Point[0] + 1]
  173.                                                               + image[L_Center_Point[1] + 2][L_Center_Point[0] + 2];
  174.                         break;
  175.                     }
  176.                     case 2:
  177.                     {
  178.                         L_Pixel_Value_Sum = L_Pixel_Value_Sum - image[L_Center_Point[1] + 1][L_Center_Point[0] - 3] - image[L_Center_Point[1] + 0][L_Center_Point[0] - 3]
  179.                                                               - image[L_Center_Point[1] - 1][L_Center_Point[0] - 3] - image[L_Center_Point[1] - 2][L_Center_Point[0] - 3]
  180.                                                               - image[L_Center_Point[1] - 3][L_Center_Point[0] - 3] - image[L_Center_Point[1] - 3][L_Center_Point[0] - 2]
  181.                                                               - image[L_Center_Point[1] - 3][L_Center_Point[0] - 1] - image[L_Center_Point[1] - 3][L_Center_Point[0] + 0]
  182.                                                               - image[L_Center_Point[1] - 3][L_Center_Point[0] + 1]
  183.                                                               + image[L_Center_Point[1] - 2][L_Center_Point[0] + 2] + image[L_Center_Point[1] - 1][L_Center_Point[0] + 2]
  184.                                                               + image[L_Center_Point[1] - 0][L_Center_Point[0] + 2] + image[L_Center_Point[1] + 1][L_Center_Point[0] + 2]
  185.                                                               + image[L_Center_Point[1] + 2][L_Center_Point[0] + 2] + image[L_Center_Point[1] + 2][L_Center_Point[0] + 1]
  186.                                                               + image[L_Center_Point[1] + 2][L_Center_Point[0] + 0] + image[L_Center_Point[1] + 2][L_Center_Point[0] - 1]
  187.                                                               + image[L_Center_Point[1] + 2][L_Center_Point[0] - 2];
  188.                         break;
  189.                     }
  190.                     case 3:
  191.                     {
  192.                         L_Pixel_Value_Sum = L_Pixel_Value_Sum - image[L_Center_Point[1] + 2][L_Center_Point[0] - 3] - image[L_Center_Point[1] + 1][L_Center_Point[0] - 3]
  193.                                                               - image[L_Center_Point[1] - 0][L_Center_Point[0] - 3] - image[L_Center_Point[1] - 1][L_Center_Point[0] - 3]
  194.                                                               - image[L_Center_Point[1] - 2][L_Center_Point[0] - 3]
  195.                                                               + image[L_Center_Point[1] + 2][L_Center_Point[0] + 2] + image[L_Center_Point[1] + 1][L_Center_Point[0] + 2]
  196.                                                               + image[L_Center_Point[1] + 0][L_Center_Point[0] + 2] + image[L_Center_Point[1] - 1][L_Center_Point[0] + 2]
  197.                                                               + image[L_Center_Point[1] - 2][L_Center_Point[0] + 2];
  198.                         break;
  199.                     }
  200.                     case 4:
  201.                     {
  202.                         L_Pixel_Value_Sum = L_Pixel_Value_Sum - image[L_Center_Point[1] + 3][L_Center_Point[0] + 1] - image[L_Center_Point[1] + 3][L_Center_Point[0] - 0]
  203.                                                               - image[L_Center_Point[1] + 3][L_Center_Point[0] - 1] - image[L_Center_Point[1] + 3][L_Center_Point[0] - 2]
  204.                                                               - image[L_Center_Point[1] + 3][L_Center_Point[0] - 3] - image[L_Center_Point[1] + 2][L_Center_Point[0] - 3]
  205.                                                               - image[L_Center_Point[1] + 1][L_Center_Point[0] - 3] - image[L_Center_Point[1] - 0][L_Center_Point[0] - 3]
  206.                                                               - image[L_Center_Point[1] - 1][L_Center_Point[0] - 3]
  207.                                                               + image[L_Center_Point[1] + 2][L_Center_Point[0] + 2] + image[L_Center_Point[1] + 1][L_Center_Point[0] + 2]
  208.                                                               + image[L_Center_Point[1] - 0][L_Center_Point[0] + 2] + image[L_Center_Point[1] - 1][L_Center_Point[0] + 2]
  209.                                                               + image[L_Center_Point[1] - 2][L_Center_Point[0] + 2] + image[L_Center_Point[1] - 2][L_Center_Point[0] + 1]
  210.                                                               + image[L_Center_Point[1] - 2][L_Center_Point[0] + 0] + image[L_Center_Point[1] - 2][L_Center_Point[0] - 1]
  211.                                                               + image[L_Center_Point[1] - 2][L_Center_Point[0] - 2];
  212.                         break;
  213.                     }
  214.                 }
  215.             }
  216.             else
  217.             {
  218.                 for (j = 0; j < 25; j++)    //第一个阈值将25个点全部加一遍,后续点的阈值根据生长方向计算
  219.                 {
  220.                     L_Pixel_Value_Sum += image[L_Center_Point[1] + Square_0[j][1]][L_Center_Point[0] + Square_0[j][0]];
  221.                 }
  222.             }
  223.             L_Thres = (L_Pixel_Value_Sum + Thres_Interfere) / Thres_Num_Interfere;   //阈值为25个点灰度值的平均值
  224.             L_Thres -= clip_value;              //将得到的灰度阈值减去一个经验值,用来优化判定
  225.             //这里为反向平滑滤波(不知网上是否有这种算法,此处为本人原创),用于适应由较暗到(局部)高亮区域,或反之的光线情况。后续详细讲解原理
  226.             if(Thres_Filiter_Flag_1 == 1 || Thres_Filiter_Flag_2 == 1)
  227.             {
  228.                 if(L_Data_Statics > 3)
  229.                 {
  230.                     L_Thres = L_Thres * 1.3f - L_Thres_Record[L_Data_Statics - 1] * 0.2f - L_Thres_Record[L_Data_Statics - 2] * 0.1f;
  231.                 }
  232.             }
  233.             L_Thres_Record[L_Data_Statics] = L_Thres;
  234.             L_Data_Statics++;                   //每找到一个点统计个数+1
  235.             L_Judge_Again:    //L_Judge_Again 与 goto 配合使用
  236.             if(L_Stop_Flag == 0)
  237.             {
  238.                 L_Front_Value = image[L_Center_Point[1] + L_Face_Dir[L_Dir][1]][L_Center_Point[0] + L_Face_Dir[L_Dir][0]];          //记录面向的前方点的灰度值
  239.                 L_Front_L_Value = image[L_Center_Point[1] + L_Face_Dir_L[L_Dir][1]][L_Center_Point[0] + L_Face_Dir_L[L_Dir][0]];    //记录面向的左前方点的灰度值
  240.                 if((float)L_Front_Value < L_Thres)     //面向的前方点是黑色
  241.                 {
  242.                     L_Dir = (L_Dir + 1) % 4;    //需右转一次
  243.                     L_Turn_Num ++;
  244.                     if(L_Turn_Num == 4)        //死区处理
  245.                     {
  246.                         L_Stop_Flag = 1;       //当前后左右都是黑色时,进入死区,停止左侧爬线
  247.                     }
  248.                     goto L_Judge_Again;
  249.                 }
  250.                 else if((float)L_Front_L_Value < L_Thres)   //左前方点是黑色,前方点是白色
  251.                 {
  252.                     L_Center_Point[0] += L_Face_Dir[L_Dir][0];
  253.                     L_Center_Point[1] += L_Face_Dir[L_Dir][1];      //向前走一步
  254.                     l_dir[L_Data_Statics - 1] = (L_Face_Dir[L_Dir][0] * 3) - L_Face_Dir[L_Dir][1];
  255.                     L_Turn_Num = 0;
  256.                 }
  257.                 else        //左前方和前方都是白色点
  258.                 {
  259.                     L_Center_Point[0] += L_Face_Dir_L[L_Dir][0];
  260.                     L_Center_Point[1] += L_Face_Dir_L[L_Dir][1];        //向左前方走一步
  261.                     l_dir[L_Data_Statics - 1] = (L_Face_Dir_L[L_Dir][0] * 3) - L_Face_Dir_L[L_Dir][1];
  262.                     L_Dir = (L_Dir + 3) % 4;        //左转一次
  263.                     L_Turn_Num = 0;
  264.                 }
  265.                 if(L_Data_Statics >= 4)     //O环处理,即转了一圈后回到原处,也是一种死区,当立即停止爬线
  266.                 {
  267.                     if(l_line[L_Data_Statics][0] == l_line[L_Data_Statics - 4][0]&&
  268.                        l_line[L_Data_Statics][1] == l_line[L_Data_Statics - 4][1])
  269.                     {
  270.                         L_Stop_Flag = 1;
  271.                     }
  272.                 }
  273.             }
  274.         }
  275.         //右侧与左侧同理,代码也类似,理解左侧后右侧就很简单
  276.         if(R_Stop_Flag == 0)
  277.         {
  278.             r_line[R_Data_Statics][0] = R_Center_Point[0];
  279.             r_line[R_Data_Statics][1] = R_Center_Point[1];
  280.             if(R_Data_Statics != 0)
  281.             {
  282.                 switch(r_dir[R_Data_Statics - 1])
  283.                 {
  284.                     case 1:
  285.                     {
  286.                         R_Pixel_Value_Sum = R_Pixel_Value_Sum - image[R_Center_Point[1] + 3][R_Center_Point[0] + 2] - image[R_Center_Point[1] + 3][R_Center_Point[0] + 1]
  287.                                                               - image[R_Center_Point[1] + 3][R_Center_Point[0] + 0] - image[R_Center_Point[1] + 3][R_Center_Point[0] - 1]
  288.                                                               - image[R_Center_Point[1] + 3][R_Center_Point[0] - 2]
  289.                                                               + image[R_Center_Point[1] - 2][R_Center_Point[0] + 2] + image[R_Center_Point[1] - 2][R_Center_Point[0] + 1]
  290.                                                               + image[R_Center_Point[1] - 2][R_Center_Point[0] + 0] + image[R_Center_Point[1] - 2][R_Center_Point[0] - 1]
  291.                                                               + image[R_Center_Point[1] - 2][R_Center_Point[0] - 2];
  292.                         break;
  293.                     }
  294.                     case -2:
  295.                     {
  296.                         R_Pixel_Value_Sum = R_Pixel_Value_Sum - image[R_Center_Point[1] - 1][R_Center_Point[0] + 3] - image[R_Center_Point[1] - 0][R_Center_Point[0] + 3]
  297.                                                               - image[R_Center_Point[1] + 1][R_Center_Point[0] + 3] - image[R_Center_Point[1] + 2][R_Center_Point[0] + 3]
  298.                                                               - image[R_Center_Point[1] + 3][R_Center_Point[0] + 3] - image[R_Center_Point[1] + 3][R_Center_Point[0] + 2]
  299.                                                               - image[R_Center_Point[1] + 3][R_Center_Point[0] + 1] - image[R_Center_Point[1] + 3][R_Center_Point[0] - 0]
  300.                                                               - image[R_Center_Point[1] + 3][R_Center_Point[0] - 1]
  301.                                                               + image[R_Center_Point[1] + 2][R_Center_Point[0] - 2] + image[R_Center_Point[1] + 1][R_Center_Point[0] - 2]
  302.                                                               + image[R_Center_Point[1] + 0][R_Center_Point[0] - 2] + image[R_Center_Point[1] - 1][R_Center_Point[0] - 2]
  303.                                                               + image[R_Center_Point[1] - 2][R_Center_Point[0] - 2] + image[R_Center_Point[1] - 2][R_Center_Point[0] - 1]
  304.                                                               + image[R_Center_Point[1] - 2][R_Center_Point[0] - 0] + image[R_Center_Point[1] - 2][R_Center_Point[0] + 1]
  305.                                                               + image[R_Center_Point[1] - 2][R_Center_Point[0] + 2];
  306.                         break;
  307.                     }
  308.                     case -3:
  309.                     {
  310.                         R_Pixel_Value_Sum = R_Pixel_Value_Sum - image[R_Center_Point[1] - 2][R_Center_Point[0] + 3] - image[R_Center_Point[1] - 1][R_Center_Point[0] + 3]
  311.                                                               - image[R_Center_Point[1] + 0][R_Center_Point[0] + 3] - image[R_Center_Point[1] + 1][R_Center_Point[0] + 3]
  312.                                                               - image[R_Center_Point[1] + 2][R_Center_Point[0] + 3]
  313.                                                               + image[R_Center_Point[1] - 2][R_Center_Point[0] - 2] + image[R_Center_Point[1] - 1][R_Center_Point[0] - 2]
  314.                                                               + image[R_Center_Point[1] - 0][R_Center_Point[0] - 2] + image[R_Center_Point[1] + 1][R_Center_Point[0] - 2]
  315.                                                               + image[R_Center_Point[1] + 2][R_Center_Point[0] - 2];
  316.                         break;
  317.                     }
  318.                     case -4:
  319.                     {
  320.                         R_Pixel_Value_Sum = R_Pixel_Value_Sum - image[R_Center_Point[1] - 3][R_Center_Point[0] - 1] - image[R_Center_Point[1] - 3][R_Center_Point[0] + 0]
  321.                                                               - image[R_Center_Point[1] - 3][R_Center_Point[0] + 1] - image[R_Center_Point[1] - 3][R_Center_Point[0] + 2]
  322.                                                               - image[R_Center_Point[1] - 3][R_Center_Point[0] + 3] - image[R_Center_Point[1] - 2][R_Center_Point[0] + 3]
  323.                                                               - image[R_Center_Point[1] - 1][R_Center_Point[0] + 3] - image[R_Center_Point[1] + 0][R_Center_Point[0] + 3]
  324.                                                               - image[R_Center_Point[1] + 1][R_Center_Point[0] + 3]
  325.                                                               + image[R_Center_Point[1] - 2][R_Center_Point[0] - 2] + image[R_Center_Point[1] - 1][R_Center_Point[0] - 2]
  326.                                                               + image[R_Center_Point[1] + 0][R_Center_Point[0] - 2] + image[R_Center_Point[1] + 1][R_Center_Point[0] - 2]
  327.                                                               + image[R_Center_Point[1] + 2][R_Center_Point[0] - 2] + image[R_Center_Point[1] + 2][R_Center_Point[0] - 1]
  328.                                                               + image[R_Center_Point[1] + 2][R_Center_Point[0] - 0] + image[R_Center_Point[1] + 2][R_Center_Point[0] + 1]
  329.                                                               + image[R_Center_Point[1] + 2][R_Center_Point[0] + 2];
  330.                         break;
  331.                     }
  332.                     case -1:
  333.                     {
  334.                         R_Pixel_Value_Sum = R_Pixel_Value_Sum - image[R_Center_Point[1] - 3][R_Center_Point[0] - 2] - image[R_Center_Point[1] - 3][R_Center_Point[0] - 1]
  335.                                                               - image[R_Center_Point[1] - 3][R_Center_Point[0] + 0] - image[R_Center_Point[1] - 3][R_Center_Point[0] + 1]
  336.                                                               - image[R_Center_Point[1] - 3][R_Center_Point[0] + 2]
  337.                                                               + image[R_Center_Point[1] + 2][R_Center_Point[0] - 2] + image[R_Center_Point[1] + 2][R_Center_Point[0] - 1]
  338.                                                               + image[R_Center_Point[1] + 2][R_Center_Point[0] + 0] + image[R_Center_Point[1] + 2][R_Center_Point[0] + 1]
  339.                                                               + image[R_Center_Point[1] + 2][R_Center_Point[0] + 2];
  340.                         break;
  341.                     }
  342.                     case 2:
  343.                     {
  344.                         R_Pixel_Value_Sum = R_Pixel_Value_Sum - image[R_Center_Point[1] + 1][R_Center_Point[0] - 3] - image[R_Center_Point[1] + 0][R_Center_Point[0] - 3]
  345.                                                               - image[R_Center_Point[1] - 1][R_Center_Point[0] - 3] - image[R_Center_Point[1] - 2][R_Center_Point[0] - 3]
  346.                                                               - image[R_Center_Point[1] - 3][R_Center_Point[0] - 3] - image[R_Center_Point[1] - 3][R_Center_Point[0] - 2]
  347.                                                               - image[R_Center_Point[1] - 3][R_Center_Point[0] - 1] - image[R_Center_Point[1] - 3][R_Center_Point[0] + 0]
  348.                                                               - image[R_Center_Point[1] - 3][R_Center_Point[0] + 1]
  349.                                                               + image[R_Center_Point[1] - 2][R_Center_Point[0] + 2] + image[R_Center_Point[1] - 1][R_Center_Point[0] + 2]
  350.                                                               + image[R_Center_Point[1] - 0][R_Center_Point[0] + 2] + image[R_Center_Point[1] + 1][R_Center_Point[0] + 2]
  351.                                                               + image[R_Center_Point[1] + 2][R_Center_Point[0] + 2] + image[R_Center_Point[1] + 2][R_Center_Point[0] + 1]
  352.                                                               + image[R_Center_Point[1] + 2][R_Center_Point[0] + 0] + image[R_Center_Point[1] + 2][R_Center_Point[0] - 1]
  353.                                                               + image[R_Center_Point[1] + 2][R_Center_Point[0] - 2];
  354.                         break;
  355.                     }
  356.                     case 3:
  357.                     {
  358.                         R_Pixel_Value_Sum = R_Pixel_Value_Sum - image[R_Center_Point[1] + 2][R_Center_Point[0] - 3] - image[R_Center_Point[1] + 1][R_Center_Point[0] - 3]
  359.                                                               - image[R_Center_Point[1] - 0][R_Center_Point[0] - 3] - image[R_Center_Point[1] - 1][R_Center_Point[0] - 3]
  360.                                                               - image[R_Center_Point[1] - 2][R_Center_Point[0] - 3]
  361.                                                               + image[R_Center_Point[1] + 2][R_Center_Point[0] + 2] + image[R_Center_Point[1] + 1][R_Center_Point[0] + 2]
  362.                                                               + image[R_Center_Point[1] + 0][R_Center_Point[0] + 2] + image[R_Center_Point[1] - 1][R_Center_Point[0] + 2]
  363.                                                               + image[R_Center_Point[1] - 2][R_Center_Point[0] + 2];
  364.                         break;
  365.                     }
  366.                     case 4:
  367.                     {
  368.                         R_Pixel_Value_Sum = R_Pixel_Value_Sum - image[R_Center_Point[1] + 3][R_Center_Point[0] + 1] - image[R_Center_Point[1] + 3][R_Center_Point[0] - 0]
  369.                                                               - image[R_Center_Point[1] + 3][R_Center_Point[0] - 1] - image[R_Center_Point[1] + 3][R_Center_Point[0] - 2]
  370.                                                               - image[R_Center_Point[1] + 3][R_Center_Point[0] - 3] - image[R_Center_Point[1] + 2][R_Center_Point[0] - 3]
  371.                                                               - image[R_Center_Point[1] + 1][R_Center_Point[0] - 3] - image[R_Center_Point[1] - 0][R_Center_Point[0] - 3]
  372.                                                               - image[R_Center_Point[1] - 1][R_Center_Point[0] - 3]
  373.                                                               + image[R_Center_Point[1] + 2][R_Center_Point[0] + 2] + image[R_Center_Point[1] + 1][R_Center_Point[0] + 2]
  374.                                                               + image[R_Center_Point[1] - 0][R_Center_Point[0] + 2] + image[R_Center_Point[1] - 1][R_Center_Point[0] + 2]
  375.                                                               + image[R_Center_Point[1] - 2][R_Center_Point[0] + 2] + image[R_Center_Point[1] - 2][R_Center_Point[0] + 1]
  376.                                                               + image[R_Center_Point[1] - 2][R_Center_Point[0] + 0] + image[R_Center_Point[1] - 2][R_Center_Point[0] - 1]
  377.                                                               + image[R_Center_Point[1] - 2][R_Center_Point[0] - 2];
  378.                         break;
  379.                     }
  380.                 }
  381.             }
  382.             else
  383.             {
  384.                 for (j = 0; j < 25; j++)
  385.                 {
  386.                     R_Pixel_Value_Sum += image[R_Center_Point[1] + Square_0[j][1]][R_Center_Point[0] + Square_0[j][0]];
  387.                 }
  388.             }
  389.             R_Thres = (R_Pixel_Value_Sum + Thres_Interfere) / Thres_Num_Interfere;
  390.             R_Thres -= clip_value;
  391.             if(Thres_Filiter_Flag_1 == 1 || Thres_Filiter_Flag_2 == 1)
  392.             {
  393.                 if(R_Data_Statics > 3)
  394.                 {
  395.                     R_Thres = R_Thres * 1.3f - R_Thres_Record[R_Data_Statics - 1] * 0.2f - R_Thres_Record[R_Data_Statics - 2] * 0.1f;
  396.                 }
  397.             }
  398.             R_Thres_Record[R_Data_Statics] = R_Thres;
  399.             R_Data_Statics++;
  400.             R_Judgme_Again:
  401.             if(R_Stop_Flag == 0)
  402.             {
  403.                 R_Front_Value = image[R_Center_Point[1] + R_Face_Dir[R_Dir][1]][R_Center_Point[0] + R_Face_Dir[R_Dir][0]];
  404.                 R_Front_R_Value = image[R_Center_Point[1] + R_Face_Dir_R[R_Dir][1]][R_Center_Point[0] + R_Face_Dir_R[R_Dir][0]];
  405.                 if((float)R_Front_Value < R_Thres)
  406.                 {
  407.                     R_Dir = (R_Dir + 3) % 4;
  408.                     R_Turn_Num ++;
  409.                     if(R_Turn_Num == 4)
  410.                     {
  411.                         R_Stop_Flag = 1;
  412.                     }
  413.                     goto R_Judgme_Again;
  414.                 }
  415.                 else if((float)R_Front_R_Value < R_Thres)
  416.                 {
  417.                     R_Center_Point[0] += R_Face_Dir[R_Dir][0];
  418.                     R_Center_Point[1] += R_Face_Dir[R_Dir][1];
  419.                     r_dir[R_Data_Statics - 1] = R_Face_Dir[R_Dir][0] * 3 - R_Face_Dir[R_Dir][1];
  420.                     R_Turn_Num = 0;
  421.                 }
  422.                 else
  423.                 {
  424.                     R_Center_Point[0] += R_Face_Dir_R[R_Dir][0];
  425.                     R_Center_Point[1] += R_Face_Dir_R[R_Dir][1];
  426.                     r_dir[R_Data_Statics - 1] = R_Face_Dir_R[R_Dir][0] * 3 - R_Face_Dir_R[R_Dir][1];
  427.                     R_Dir = (R_Dir + 1) % 4;
  428.                     R_Turn_Num = 0;
  429.                 }
  430.                 if(R_Data_Statics >= 4)
  431.                 {
  432.                     if(r_line[R_Data_Statics][0] == r_line[R_Data_Statics - 4][0]&&
  433.                        r_line[R_Data_Statics][1] == r_line[R_Data_Statics - 4][1])
  434.                     {
  435.                         R_Stop_Flag = 1;
  436.                     }
  437.                 }
  438.             }
  439.         }
  440.         if(L_Stop_Flag == 0 && R_Stop_Flag == 0)
  441.         {
  442.             if ((My_ABS(r_line[R_Data_Statics - 1][0] - l_line[L_Data_Statics - 1][0]) <= 1)
  443.                 && (My_ABS(r_line[R_Data_Statics - 1][1] - l_line[L_Data_Statics - 1][1]) <= 1))        //两侧爬线相遇,退出循环,一张图像爬线结束
  444.             {
  445.                 *y_meet = (r_line[R_Data_Statics - 1][1] + l_line[L_Data_Statics - 1][1]) >> 1;  //记录相遇点Y
  446.                 *x_meet = (r_line[R_Data_Statics - 1][0] + l_line[L_Data_Statics - 1][0]) >> 1;  //记录相遇点X
  447.                 break;
  448.             }
  449.         }
  450.         //有一侧存在死区时,对相遇点的判定放宽松一些,防止实际相遇但没有判定出,导致爬线紊乱的情况
  451.         else
  452.         {
  453.             if ((My_ABS(r_line[R_Data_Statics - 1][0] - l_line[L_Data_Statics - 1][0]) <= 3)
  454.                 && (My_ABS(r_line[R_Data_Statics - 1][1] - l_line[L_Data_Statics - 1][1]) <= 3))        //两侧爬线相遇,退出循环,一张图像爬线结束
  455.             {
  456.                 *y_meet = (r_line[R_Data_Statics - 1][1] + l_line[L_Data_Statics - 1][1]) >> 1;  //记录相遇点Y
  457.                 *x_meet = (r_line[R_Data_Statics - 1][0] + l_line[L_Data_Statics - 1][0]) >> 1;  //记录相遇点X
  458.                 break;
  459.             }
  460.         }
  461.     }
  462.     L_Stop_Flag = 0;
  463.     R_Stop_Flag = 0;
  464.     *l_stastic = L_Data_Statics;    //记录左侧边线点个数
  465.     *r_stastic = R_Data_Statics;    //记录右侧边线点个数
  466. }
复制代码
        最后还是建议不要一味地拷代码,一定要自己一点一点敲出来,真正把算法明白透彻才可谈优化算法。车赛不是拷代码大赛,更多的是学新东西、学新思路。
        同时补充一点,核心算法的稳定性,决定了后续元素处理的方法;核心算法的速度,决定了团体代码的上限和部分控制上限;核心算法提取的信息量,决定了后续处理元素的难易水平。以是务必在此部分下功夫,确保明白透彻。

智能车摄像头开源—1.1 核心代码:自顺应八向迷宫(上)
智能车摄像头开源—1.2 核心算法:自顺应八向迷宫(下)
智能车摄像头开源—2 摄像头基础参数设置经验分享
智能车摄像头开源—3 图像基础处理、迭代优化与结果展示

免责声明:如果侵犯了您的权益,请联系站长,我们会及时删除侵权内容,谢谢合作!更多信息从访问主页:qidao123.com:ToB企服之家,中国第一个企服评测及商务社交产业平台。
回复

使用道具 举报

0 个回复

倒序浏览

快速回复

您需要登录后才可以回帖 登录 or 立即注册

本版积分规则

万有斥力

金牌会员
这个人很懒什么都没写!
快速回复 返回顶部 返回列表