measure_filter.h 5.3 KB

123456789101112131415161718192021222324252627282930313233343536373839404142434445464748495051525354555657585960616263646566676869707172737475767778798081828384858687888990919293949596979899100101102103104105106107108109110111112113114115116117118119120121122123124125126127128129130
  1. /*
  2. * @Description: 测量结果滤波器
  3. * @Author: yct
  4. * @Date: 2021-03-15 14:41:46
  5. * @LastEditTime: 2021-07-30 10:33:47
  6. * @LastEditors: yct
  7. */
  8. #ifndef MEASURE_FILTER_HH
  9. #define MEASURE_FILTER_HH
  10. #include "../tool/singleton.h"
  11. #include "../tool/common_data.h"
  12. #include "../error_code/error_code.h"
  13. #include <map>
  14. #include <deque>
  15. #include <chrono>
  16. #include <vector>
  17. #include <algorithm>
  18. #include "glog/logging.h"
  19. #define FILTER_SIZE 6
  20. #define MAX_QUEUE_SIZE 12
  21. #define MAX_TIME_INTERVAL_MILLI 9000
  22. class Measure_filter : public Singleton<Measure_filter>
  23. {
  24. // 子类必须把父类设定为友元函数,这样父类才能使用子类的私有构造函数。
  25. friend class Singleton<Measure_filter>;
  26. public:
  27. // 必须关闭拷贝构造和赋值构造,只能通过 get_instance 函数来进行操作唯一的实例。
  28. Measure_filter(const Measure_filter &) = delete;
  29. Measure_filter &operator=(const Measure_filter &) = delete;
  30. ~Measure_filter() = default;
  31. // 更新测量数据回调
  32. void update_data(int terminal_id, Common_data::Car_wheel_information_stamped data)
  33. {
  34. if(!data.wheel_data.correctness)
  35. return;
  36. // LOG(INFO) << data.wheel_data.to_string();
  37. // 未创建队列
  38. if (m_measure_results_map.find(terminal_id) == m_measure_results_map.end())
  39. {
  40. std::deque<Common_data::Car_wheel_information_stamped> t_deque;
  41. t_deque.push_back(data);
  42. m_measure_results_map.insert(std::pair<int, std::deque<Common_data::Car_wheel_information_stamped>>(terminal_id, t_deque));
  43. }
  44. else//已创建
  45. {
  46. m_measure_results_map[terminal_id].push_back(data);
  47. // 剔除超时数据
  48. while(std::chrono::duration_cast<std::chrono::milliseconds>(std::chrono::system_clock::now()-m_measure_results_map[terminal_id].front().measure_time).count() > MAX_TIME_INTERVAL_MILLI)
  49. {
  50. m_measure_results_map[terminal_id].pop_front();
  51. }
  52. // 维持队列长度
  53. while(m_measure_results_map[terminal_id].size() > MAX_QUEUE_SIZE)
  54. {
  55. m_measure_results_map[terminal_id].pop_front();
  56. }
  57. }
  58. }
  59. // 获取滤波后最新测量数据
  60. Error_manager get_filtered_wheel_information(int terminal_id, Common_data::Car_wheel_information& result)
  61. {
  62. // 检查数据量
  63. if(m_measure_results_map.find(terminal_id) == m_measure_results_map.end() || m_measure_results_map.find(terminal_id)->second.size() < FILTER_SIZE)
  64. {
  65. return Error_manager(WJ_FILTER_LACK_OF_RESULT, MINOR_ERROR,
  66. (std::string("缺少足够用于滤波的结果")+
  67. std::to_string(terminal_id)+
  68. std::string(m_measure_results_map.find(terminal_id) == m_measure_results_map.end()?
  69. "end":
  70. std::to_string(m_measure_results_map.find(terminal_id)->second.size()))).c_str());
  71. }
  72. // 填充待滤波结果到数组
  73. std::deque<Common_data::Car_wheel_information_stamped>* tp_deque = &m_measure_results_map.find(terminal_id)->second;
  74. std::vector<Common_data::Car_wheel_information_stamped> t_result_vec;
  75. Common_data::Car_wheel_information_stamped t_avg_result;
  76. for (std::deque<Common_data::Car_wheel_information_stamped>::reverse_iterator t_iter = tp_deque->rbegin(); t_iter != tp_deque->rend(); t_iter++)
  77. {
  78. t_result_vec.push_back(*t_iter);
  79. t_avg_result += *t_iter;
  80. if (t_result_vec.size() >= FILTER_SIZE)
  81. break;
  82. }
  83. if(t_result_vec.size()<FILTER_SIZE)
  84. return Error_manager(WJ_FILTER_LACK_OF_RESULT, MINOR_ERROR, "结果缺失");
  85. // 各结果计算平均,并打分
  86. t_avg_result /= t_result_vec.size();
  87. std::vector<float> t_score_vec, t_sorted_score_vec;
  88. for (size_t i = 0; i < t_result_vec.size(); i++)
  89. {
  90. t_score_vec.push_back((t_result_vec[i] - t_avg_result).calc_score());
  91. t_sorted_score_vec.push_back((t_result_vec[i] - t_avg_result).calc_score());
  92. }
  93. // 排序,丢掉最高两个(即与均值偏差最大两个,剩下平均获得最终值)
  94. std::sort(t_sorted_score_vec.begin(), t_sorted_score_vec.end());
  95. Common_data::Car_wheel_information_stamped t_final_result;
  96. int t_result_count = 0;
  97. for (size_t i = 0; i < t_score_vec.size(); i++)
  98. {
  99. if(t_score_vec[i] > t_sorted_score_vec[t_sorted_score_vec.size()-3])
  100. {
  101. continue;
  102. }else{
  103. t_final_result += t_result_vec[i];
  104. t_result_count++;
  105. }
  106. }
  107. if(t_result_count<=0)
  108. return Error_manager(WJ_FILTER_FLUCTUATING, MINOR_ERROR, "结果波动");
  109. t_final_result /= t_result_count;
  110. result = t_final_result.wheel_data;
  111. // LOG(INFO) << "\navg: \n\t"<<t_avg_result.wheel_data.to_string()<<"final: \t"<<t_final_result.wheel_data.to_string();
  112. return SUCCESS;
  113. }
  114. private:
  115. // 父类的构造函数必须保护,子类的构造函数必须私有。
  116. Measure_filter() = default;
  117. // 各终端测量数据队列
  118. std::map<int, std::deque<Common_data::Car_wheel_information_stamped> > m_measure_results_map;
  119. };
  120. #endif // !MEASURE_FILTER_HH