ivdriver_lidar_rs16.cpp 4.1 KB

123456789101112131415161718192021222324252627282930313233343536373839404142434445464748495051525354555657585960616263646566676869707172737475767778798081828384858687888990919293949596979899100101102103104105106107108109110111112113114115116117118119120121122123124125126127128129130131132133134135136137138139
  1. #include "ivdriver_lidar_rs16.h"
  2. namespace iv {
  3. ivdriver_lidar_rs16::ivdriver_lidar_rs16()
  4. {
  5. }
  6. int ivdriver_lidar_rs16::processudpData(QByteArray ba, pcl::PointCloud<pcl::PointXYZI>::Ptr point_cloud)
  7. {
  8. static double fAngle_Total = 0;
  9. static double fAngle_Last = 0;
  10. int nrtn = 0;
  11. float Ang = 0;
  12. float Range = 0;
  13. int Group = 0;
  14. int pointi = 0;
  15. float wt = 0;
  16. int wt1 = 0;
  17. float V_theta[16] = {-15,-13,-11,-9,-7,-5,-3,-1,15,13,11,9,7,5,3,1};
  18. char * pstr = ba.data();
  19. wt1 = ((pstr[2 + Group * 100] *256) + pstr[ 3 + Group * 100]) ;
  20. wt = wt1/ 100.0;
  21. double fAngX = wt;
  22. if(fabs(fAngX-fAngle_Last)>300)
  23. {
  24. fAngle_Total = fAngle_Total + fabs(fabs(fAngX-fAngle_Last)-360);
  25. }
  26. else
  27. {
  28. fAngle_Total = fAngle_Total + fabs(fabs(fAngX-fAngle_Last));
  29. }
  30. fAngle_Last = fAngX;
  31. if(fAngle_Total > 360)
  32. {
  33. nrtn = 1;
  34. fAngle_Total = 0;
  35. }
  36. for (Group = 0; Group <= 11; Group++)
  37. {
  38. wt1 = ((pstr[2 + Group * 100] *256) + pstr[ 3 + Group * 100]) ;
  39. wt = wt1/ 100.0;
  40. for (pointi = 0; pointi <16; pointi++)
  41. {
  42. // Ang = (0 - wt - w * T[pointi] - H_BETA[pointi]+213) / 180.0 * Lidar_Pi;
  43. Ang = (0 - wt) / 180.0 * M_PI;
  44. Range = ((pstr[ Group * 100 + 4 + 3 * pointi] << 8) + pstr[Group * 100 + 5 + 3 * pointi]);
  45. unsigned char intensity = pstr[ Group * 100 + 6 + 3 * pointi];
  46. Range=Range* 5.0/1000.0;
  47. if(Range<150)
  48. {
  49. pcl::PointXYZI point;
  50. point.x = Range*cos(V_theta[pointi] / 180 * M_PI)*cos(Ang + mfrollang);
  51. point.y = Range*cos(V_theta[pointi] / 180 * M_PI)*sin(Ang + mfrollang);
  52. point.z = Range*sin(V_theta[pointi] / 180 * M_PI);
  53. if(mbinclix)
  54. {
  55. double y,z;
  56. y = point.y;
  57. z = point.z;
  58. point.y = y*mcos_finclinationang_xaxis +z*msin_finclinationang_xaxis;
  59. point.z = z*mcos_finclinationang_xaxis - y*msin_finclinationang_xaxis;
  60. }
  61. if(mbincliy)
  62. {
  63. double z,x;
  64. z = point.z;
  65. x = point.x;
  66. point.z = z*mcos_finclinationang_yaxis + x*msin_finclinationang_yaxis;
  67. point.x = x*mcos_finclinationang_yaxis - z*msin_finclinationang_yaxis;
  68. }
  69. point.intensity = intensity;
  70. point_cloud->points.push_back(point);
  71. ++point_cloud->width;
  72. }
  73. }
  74. wt = wt + 0.18;
  75. for (pointi = 0; pointi < 16; pointi++)
  76. {
  77. Ang = (0 - wt) / 180.0 * M_PI;
  78. // Ang = Ang+angdiff;
  79. Range = ((pstr[ Group * 100 + 52 + 3 * pointi] << 8) + pstr[Group * 100 + 53 + 3 * pointi]);
  80. unsigned char intensity = pstr[ Group * 100 + 54 + 3 * pointi];
  81. Range=Range* 5.0/1000.0;
  82. if(Range<150)
  83. {
  84. pcl::PointXYZI point;
  85. point.x = Range*cos(V_theta[pointi] / 180 * M_PI)*cos(Ang + mfrollang);
  86. point.y = Range*cos(V_theta[pointi] / 180 * M_PI)*sin(Ang + mfrollang);
  87. point.z = Range*sin(V_theta[pointi] / 180 * M_PI);
  88. if(mbinclix)
  89. {
  90. double y,z;
  91. y = point.y;z = point.z;
  92. point.y = y*mcos_finclinationang_xaxis +z*msin_finclinationang_xaxis;
  93. point.z = z*mcos_finclinationang_xaxis - y*msin_finclinationang_xaxis;
  94. }
  95. if(mbincliy)
  96. {
  97. double z,x;
  98. z = point.z;x = point.x;
  99. point.z = z*mcos_finclinationang_yaxis + x*msin_finclinationang_yaxis;
  100. point.x = x*mcos_finclinationang_yaxis - z*msin_finclinationang_yaxis;
  101. }
  102. point.intensity = intensity;
  103. point_cloud->points.push_back(point);
  104. ++point_cloud->width;
  105. }
  106. }
  107. }
  108. return nrtn;
  109. }
  110. }