CalDist.py 6.1 KB

123456789101112131415161718192021222324252627282930313233343536373839404142434445464748495051525354555657585960616263646566676869707172737475767778798081828384858687888990919293949596979899100101102103104105106107108109110111112113114115116117118119120121122123124125126127
  1. from math import radians, cos, sin, asin, sqrt
  2. import pandas as pd
  3. import numpy as np
  4. from datetime import datetime
  5. from datetime import timedelta
  6. from GpsRank import *
  7. from ProcessDfBms import *
  8. from ProcessDfGps import *
  9. from LIB.BACKEND import DBManager
  10. import DBManager
  11. #####################################配置环境分割线#################################################
  12. def GetDistInfo(input_sn,input_starttime,input_endtime):
  13. #####################################配置参数分割线#################################################
  14. dbManager = DBManager.DBManager()
  15. data_raw = dbManager.get_data(sn=input_sn, start_time=input_starttime,
  16. end_time=input_endtime)
  17. #拆包预处理
  18. df_bms_raw=data_raw['bms']
  19. df_gps_raw=data_raw['gps']
  20. df_bms=preprocess_Df_Bms(df_bms_raw)
  21. df_gps=preprocess_Df_Gps(df_gps_raw)
  22. #####################################数据预处理分割线#################################################
  23. # mode: 0:正常取数; 1:7255 取数
  24. if input_sn[0:2] == 'UD' or input_sn[0:2] == 'MG':
  25. mode = 1
  26. else:
  27. mode = 0
  28. #获取状态表,mode默认为0,mode=1放电时电流为负,mode=0充电时电流为正
  29. df_bms_drive_timetable=get_bms_drive_timetable(df_bms,mode)
  30. df_gps_drive_cycle_accum=pd.DataFrame()
  31. if len(df_bms_drive_timetable)>0:
  32. for index in range(len(df_bms_drive_timetable)):
  33. #筛选drivecycle数据
  34. drive_start_time=df_bms_drive_timetable.loc[index,'drive_start_time']#开始时间
  35. drive_end_time=df_bms_drive_timetable.loc[index,'drive_end_time']#结束时间
  36. time_condition=(df_gps['time']>drive_start_time)&(df_gps['time']<drive_end_time)#时间判断条件
  37. df_gps_drive_cycle=df_gps.loc[time_condition,:].copy()
  38. df_gps_drive_cycle=df_gps_drive_cycle.reset_index(drop=True)#重置index
  39. #计算drivecycle GPS累计里程,存入表格
  40. condition_a=df_gps_drive_cycle['deltatime']>60*3
  41. condition_b=(df_gps_drive_cycle['deltatime']>90*1)&(df_gps_drive_cycle['distance']>1000)
  42. drive_cycle_dist_array=df_gps_drive_cycle.loc[~(condition_a|condition_b),'distance'].values
  43. drive_cycle_dist_array=drive_cycle_dist_array[np.where((drive_cycle_dist_array>=1)&(drive_cycle_dist_array<1000))[0]]
  44. gps_dist=drive_cycle_dist_array.sum()
  45. df_bms_drive_timetable.loc[index,'gps_dist']=gps_dist#得到GPS路径
  46. #计算头-尾的空缺时间段对应的预估SOC
  47. if len(df_gps_drive_cycle)>2:
  48. gps_starttime=df_gps_drive_cycle.loc[1,'time']#gps开始时间
  49. gps_endtime=df_gps_drive_cycle.loc[df_gps_drive_cycle.index[-1],'time']#gps结束时间
  50. #从drive_start_time到gps开始时间,使用SOC计算的里程
  51. #gps结束时间到drive_end_time,使用SOC计算的里程
  52. unrecorded_odo_head=cal_deltasoc(df_bms,drive_start_time,gps_starttime)
  53. unrecorded_odo_tail=cal_deltasoc(df_bms,gps_endtime,drive_end_time)
  54. else:
  55. #计算数据丢失行unrecordeodo
  56. unrecorded_odo_head=cal_deltasoc(df_bms,drive_start_time,drive_end_time)
  57. unrecorded_odo_tail=0
  58. #计算中间的预估SOC
  59. predict_dist=cal_unrecorded_gps(df_gps_drive_cycle,df_bms)
  60. #计算总的预估SOC
  61. totaldist=predict_dist+unrecorded_odo_head+ unrecorded_odo_tail#得到GPS路径
  62. df_bms_drive_timetable.loc[index,'predict_dist']=totaldist
  63. else :
  64. pass
  65. #####################################统计行驶里程End#################################################
  66. #打印输出结果#
  67. index_list=list(range(len(df_bms_drive_timetable)))
  68. dist_gps=0
  69. dist_predict=0
  70. day_start_time=''#当日开始时间
  71. day_end_time=''#当日结束时间
  72. day_start_soc=0#当日开始soc
  73. day_end_soc=0#当日结束soc
  74. day_min_soc=101#当日最低soc
  75. drive_accum_soc=0#累计使用SOC
  76. if len(df_bms_drive_timetable)>0:
  77. #开始行
  78. day_start_soc=df_bms_drive_timetable.loc[1,'drive_start_soc']#开始soc
  79. day_start_time=df_bms_drive_timetable.loc[1,'drive_start_time']#开始时间
  80. #结束行
  81. day_end_time=df_bms_drive_timetable.loc[len(df_bms_drive_timetable)-1,'drive_end_time']#结束时间
  82. day_end_soc=df_bms_drive_timetable.loc[len(df_bms_drive_timetable)-1,'drive_end_soc']#结束soc
  83. for index in index_list:
  84. '''汇总里程'''
  85. dist_gps+=df_bms_drive_timetable.loc[index,'gps_dist']/1000#计算GPS里程
  86. dist_predict+=df_bms_drive_timetable.loc[index,'predict_dist']#计算预估里程
  87. drive_start_soc=df_bms_drive_timetable.loc[index,'drive_start_soc']#驾驶周期开始的soc
  88. drive_end_soc=df_bms_drive_timetable.loc[index,'drive_end_soc']#驾驶周期结束的soc
  89. day_min_soc=min(day_min_soc,drive_start_soc,drive_end_soc)
  90. delta_soc=drive_start_soc-drive_end_soc#驾驶周期SOC变化量
  91. drive_accum_soc+=abs(delta_soc)#所有drive cycle累计消耗的soc
  92. # gps_score=get_df_gps_score(input_starttime,input_endtime,df_gps)
  93. # gps_score=round(gps_score,1)
  94. #计算总里程
  95. dist_gps=round(dist_gps,3)
  96. dist_predict=round(dist_predict,3)
  97. dist_all=round(dist_gps+dist_predict,3)
  98. #输出统计结果
  99. # print ('为您查询到,从'+input_starttime+'到'+input_endtime+'时间段内:')
  100. # print('SOC变化量:'+str(df_bms['bmspacksoc'].max()-df_bms['bmspacksoc'].min())+' %')
  101. # print('行驶总里程:'+str(dist_all)+' km')
  102. return {'SN':input_sn,'range':dist_all,'accum_soc':drive_accum_soc,'day_start_soc':day_start_soc,
  103. 'day_end_soc':day_end_soc,'day_start_time':day_start_time,'day_end_time':day_end_time,
  104. 'day_min_soc':day_min_soc}
  105. # print('其中GPS信号在线时里程:'+str(dist_gps)+' km')
  106. # print('其中GPS信号掉线时预估里程:'+str(dist_predict)+' km')
  107. # print('GPS信号质量评分为:'+str(gps_score),'分\n')
  108. #####################################打印结果End#################################################