video_to_photo.py 14 KB

123456789101112131415161718192021222324252627282930313233343536373839404142434445464748495051525354555657585960616263646566676869707172737475767778798081828384858687888990919293949596979899100101102103104105106107108109110111112113114115116117118119120121122123124125126127128129130131132133134135136137138139140141142143144145146147148149150151152153154155156157158159160161162163164165166167168169170171172173174175176177178179180181182183184185186187188189190191192193194195196197198199200201202203204205206207208209210211212213214215216217218219220221222223224225226227228229230231232233234235236237238239240241242243244245246247248249250251252253254255256257258259260261262263264265266267268269270271272273274275276277278279280281282283284285286287288289290291292293294295296297298299300301302303304305306307308309310311312313314315316317318319320321322323324325326327328329330331332333334335336337338339340341342343344345346347348349350351352353354355356357358359360361362363364365366367368369370371372373374375376377378379380381382383384385386387388389390391392393394395396397398399400401402403404405406407408409410411412413414415416417418419420421422423424425426427428429430431432433434435436437438439440441442443444445446447448449450451452453454455456457458459460461462463464465466467468469470471472473474475476477478479480481482483484485486487488489490491492493494495496
  1. import os
  2. import math
  3. import cv2
  4. from demo_func import run_openpose_for_normal
  5. from PIL import Image
  6. from analyse_jump import analyse_npy_side_jump,analyse_npy_side_pingban,analyse_npy_side_juanfu,analyse_npy_side_gaotaitui,analyse_npy_side_shendun
  7. import glob
  8. import time
  9. import os
  10. def split_video(path):
  11. video_path = path
  12. os.path.split(path)
  13. tur=os.path.split(path)
  14. print(tur[len(tur)-1].split('.')[0])
  15. name=tur[len(tur)-1].split('.')[0]
  16. if os.path.exists("./capture_image/"+name)==False:
  17. os.mkdir("./capture_image/"+name)
  18. cap = cv2.VideoCapture(video_path)
  19. FPS = cap.get(5)
  20. print(FPS)
  21. c = 1
  22. frameRate = 1
  23. if os.path.exists("./capture_image/"+name)==False:
  24. os.mkdir("./capture_image/"+name)
  25. while (True):
  26. ret, frame = cap.read()
  27. if ret:
  28. if (c % frameRate == 0):
  29. cv2.imwrite("./capture_image/"+name+"/capture_image" + str(c) + '.png', frame)
  30. c += 1
  31. cv2.waitKey(0)
  32. else:
  33. break
  34. cap.release()
  35. return c
  36. def process_photo(name,c):
  37. print(name)
  38. for i in range(1,c):
  39. run_openpose_for_normal(i,name)
  40. def gather_video(path):
  41. video_path = path
  42. os.path.split(path)
  43. tur=os.path.split(path)
  44. print(tur[len(tur)-1].split('.')[0])
  45. name=tur[len(tur)-1].split('.')[0]
  46. cap = cv2.VideoCapture(video_path)
  47. fps = cap.get(5)
  48. num_frames = len(os.listdir(r'./capture_image/'+name))
  49. print(num_frames)
  50. img_array = []
  51. im=Image.open("./capture_image/"+name+"/capture_image_result1.png")
  52. #img_width = 720
  53. #img_height = 1280
  54. k=0
  55. for i in range(num_frames + 1):
  56. filename = "./capture_image/"+name+"/capture_image_result" + str(i) + ".png"
  57. k+=1
  58. img = cv2.imread(filename)
  59. if img is None:
  60. continue
  61. img_array.append(img)
  62. print(k)
  63. save_name="./video/video_a/"+name+".mp4"
  64. out = cv2.VideoWriter(save_name, cv2.VideoWriter_fourcc(*"mp4v"), fps, im.size)
  65. for i in range(len(img_array)):
  66. out.write(img_array[i])
  67. out.release()
  68. im.close()
  69. def do_analysis(a_path,video_path,mode):
  70. path = video_path
  71. os.path.split(path)
  72. tur=os.path.split(path)
  73. name=tur[len(tur)-1].split('.')[0]
  74. if os.path.exists("./cap_file/"+name)==False:
  75. os.mkdir("./cap_file/"+name)
  76. #ticks=time.time()
  77. source_npy_side=a_path
  78. output=open("./cap_file/"+name+"/file2"+".txt","w")
  79. output.close()
  80. result={}
  81. p='./capture_image/'+name+"/"
  82. lenF=int(len(glob.glob(p+'*.png'))/2)
  83. for i in range(1,lenF):#run_openpose_for_normal(i)
  84. if(mode=="jump"):
  85. result_side = analyse_npy_side_jump(name,i,source_npy_side.format(i))
  86. elif (mode=="pingban"):
  87. result_side = analyse_npy_side_pingban(name,i,source_npy_side.format(i))
  88. elif(mode=="juanfu"):
  89. result_side = analyse_npy_side_juanfu(name,i,source_npy_side.format(i))
  90. elif(mode=="gaotaitui"):
  91. result_side =analyse_npy_side_gaotaitui(name,i,source_npy_side.format(i))
  92. elif(mode=="shendun"):
  93. result_side =analyse_npy_side_shendun(name,i,source_npy_side.format(i))
  94. if(result_side!=False):
  95. result[i]=result_side
  96. output=open("./cap_file/"+name+"/file1.txt","w")
  97. for re in result:
  98. print(re,file=output)
  99. for(name,value) in result[re].items():
  100. print(name,value,sep=',',file=output)
  101. print("kv-over",file=output)
  102. output.close()
  103. def read2(fileName="./cap_file/file1.txt"):
  104. file=open(fileName)
  105. db={}
  106. key=file.readline().strip()
  107. while(key != "kv-over\n"):
  108. if not key:
  109. break
  110. person={}
  111. field=file.readline().strip()
  112. while field != "kv-over":
  113. name, value= field.split(",")
  114. person[name]=value.strip()
  115. field=file.readline().strip()
  116. if(field=="kv-over\n"):
  117. break
  118. db[key]=person
  119. key = file.readline().strip()
  120. print("db==== " , db)
  121. return db
  122. def ignore_data_pingban(result):
  123. result_ig=result
  124. k="0"
  125. last={}
  126. max_index=1
  127. for re in result:
  128. if k=="0":
  129. k=re
  130. else:
  131. if(max_index<int(re)):
  132. max_index=int(re)
  133. max_angle=abs(float(result_ig[k]["s_e_angle"])-float(result_ig[k]["e_w_angle"]))
  134. max_num=int(k)
  135. less_angle=100
  136. less_num=1
  137. max_less_a=1
  138. max_less_num=k
  139. for i in result_ig:
  140. if (int(i)>=max_index/3 and int(i)<=max_index/3*2):
  141. val=math.fabs( math.fabs(float(result_ig[i]["hip_knee_angle"]))-math.fabs(float(result_ig[i]["knee_ankle_angle"])))
  142. an_=abs(float(result_ig[i]["s_e_angle"])-float(result_ig[i]["e_w_angle"]))
  143. if(math.fabs( math.fabs(float(result_ig[i]["hip_knee_angle"]))-math.fabs(float(result_ig[i]["knee_ankle_angle"])))<less_angle):
  144. less_angle=val
  145. less_num=int(i)
  146. if(an_<max_angle):
  147. max_angle=an_
  148. max_num=i
  149. if(math.fabs( math.fabs(float(result_ig[i]["hip_knee_angle"]))-math.fabs(float(result_ig[i]["knee_ankle_angle"])))>max_less_a):
  150. max_less_an=val
  151. max_less_num=int(i)
  152. re={}
  153. re["1"]={}
  154. re["1"]["num"]=max_num
  155. re["1"]["val"]=max_angle
  156. re["2"]={}
  157. re["2"]["num"]=less_num
  158. re["2"]["val"]=less_angle
  159. re["3"]={}
  160. re["3"]["num"]=max_less_num
  161. re["3"]["val"]=max_less_a
  162. #1表示手肘和肩膀最小角度
  163. #2表示最小的身体角度差
  164. #3表示最大的身体角度差
  165. #1用来说明手部力量 3-2说明稳定性 2说明最好的状态 3说明最差的状态
  166. print(re)
  167. return result_ig,re
  168. def ignore_data_juanfu(result):
  169. result_ig=result
  170. k="0"
  171. last={}
  172. max_index=1
  173. for re in result:
  174. if k=="0":
  175. k=re
  176. else:
  177. if(max_index<int(re)):
  178. max_index=int(re)
  179. max_angle=result[k]["nose_neck_angle"]
  180. max_num=int(k)
  181. less_h_k_a=result[k]["hip_knee_angle"]
  182. less_h_k_num=k
  183. max_h_k_a=result[k]["hip_knee_angle"]
  184. max_h_k_num=k
  185. less_k_a_a=result[k]["knee_ankle_angle"]
  186. less_k_a_num=k
  187. max_k_a_a=result[k]["knee_ankle_angle"]
  188. max_k_a_num=k
  189. max_neck_hip_a=result[k]["neck_hip_angle"]
  190. max_neck_hip_num=k
  191. for i in result_ig:
  192. if(float(result_ig[i]["nose_neck_angle"])>float(max_angle)):
  193. max_num=int(i)
  194. max_angle=result_ig[i]["nose_neck_angle"]
  195. if(float(result_ig[i]["hip_knee_angle"])>float(max_h_k_a)):
  196. max_h_k_num=int(i)
  197. max_h_k_a=result_ig[i]["hip_knee_angle"]
  198. if(float(result_ig[i]["knee_ankle_angle"])>float(max_k_a_a)):
  199. max_k_a_num=int(i)
  200. max_k_a_a=result_ig[i]["knee_ankle_angle"]
  201. if(float(result_ig[i]["neck_hip_angle"])>float(max_neck_hip_a)):
  202. max_neck_hip_num=int(i)
  203. max_neck_hip_a=result_ig[i]["neck_hip_angle"]
  204. if(float(result_ig[i]["hip_knee_angle"])<float(less_h_k_a)):
  205. less_h_k_num=int(i)
  206. less_h_k_a=result_ig[i]["hip_knee_angle"]
  207. if(float(result_ig[i]["knee_ankle_angle"])<float(less_k_a_a)):
  208. less_k_a_num=int(i)
  209. less_k_a_a=result_ig[i]["knee_ankle_angle"]
  210. re={}
  211. re["1"]={}
  212. re["1"]["num"]=max_num
  213. re["1"]["val"]=max_angle
  214. re["2"]={}
  215. re["2"]["num"]=max_neck_hip_num
  216. re["2"]["val"]=max_neck_hip_a
  217. re["3"]={}
  218. re["3"]["num"]=float(max_h_k_a)-float(less_h_k_a)
  219. re["3"]["num1"]=max_h_k_num
  220. re["3"]["num2"]=less_h_k_num
  221. re["3"]["val1"]=max_h_k_a
  222. re["3"]["val2"]=less_h_k_a
  223. re["4"]={}
  224. re["4"]["num"]=float(max_k_a_a)-float(less_k_a_a)
  225. re["4"]["num1"]=max_k_a_num
  226. re["4"]["num2"]=less_k_a_num
  227. re["4"]["val1"]=max_k_a_a
  228. re["4"]["val2"]=less_k_a_a
  229. #1表示头部和脖子角度
  230. #2表示脖子和臀部角度
  231. #3表示hip-knee的极差 稳定性
  232. #4表示knee-ankle的极差 稳定性
  233. #3的数据误差太大 待用性待考察
  234. print(re)
  235. return result_ig,re
  236. def ignore_data_gaotaitui(result):
  237. result_ig=result
  238. k="0"
  239. last={}
  240. max_index=1
  241. for re in result:
  242. if k=="0":
  243. k=re
  244. else:
  245. if(max_index<int(re)):
  246. max_index=int(re)
  247. max_leg_knee_r_angle=result_ig[k]["knee_ankle_r_angle"]
  248. max_leg_knee_l_angle=result_ig[k]["knee_ankle_l_angle"]
  249. max_leg_knee_angle_num=int(k)
  250. panduan={}
  251. for i in result_ig:
  252. if(float(result_ig[i]["l_d"])<30):
  253. if(float(result_ig[i]["knee_ankle_r_angle"])>120 and float(result_ig[i]["hip_knee_r_angle"])<90):
  254. panduan[i]=result_ig[i]["r_d"]
  255. f=0
  256. #ma=0
  257. ma=panduan[k]
  258. key_k=k
  259. for key in panduan:
  260. if f==0:
  261. ma=panduan[key]
  262. f=1
  263. key_k=key
  264. elif float(ma)>float(panduan[key]):
  265. key_k=key
  266. ma=panduan[key]
  267. re={}
  268. re["1"]={}
  269. re["1"]["num"]=key_k
  270. re["1"]["val"]=ma
  271. re["1"]["hip_knee_r_angle"]=result_ig[key_k]["hip_knee_r_angle"]
  272. re["1"]["knee_ankle_r_angle"]=result_ig[key_k]["knee_ankle_r_angle"]
  273. return result_ig,re
  274. def ignore_data_jump(result):
  275. result_ig=result
  276. k="0"
  277. last={}
  278. max_index=1
  279. for re in result:
  280. if k=="0":
  281. k=re
  282. else:
  283. if(max_index<int(re)):
  284. max_index=int(re)
  285. print(result_ig)
  286. max_leg_knee_angle=result_ig[k]["knee_ankle_angle"]
  287. max_leg_knee_angle_num=int(k)
  288. less_angle=100
  289. less_num=1
  290. ankle=result_ig[k]["leg_heigh"]
  291. ankle_num=int(k)
  292. max_leg_knee_hip_knee=math.fabs(float(result_ig[k]["hip_knee_angle"])-float(result_ig[k]["knee_ankle_angle"]))
  293. max_l_h_num=1
  294. hip_hei_init=float(result_ig[k]["hip_heigh"])
  295. #print("max")
  296. #print(max_index)
  297. #print(hip_hei_init)
  298. for i in result_ig:
  299. if int(i)<=max_index/5:
  300. if float(max_leg_knee_angle)<float(result_ig[i]["knee_ankle_angle"]):
  301. max_leg_knee_angle_num=int(i)
  302. max_leg_knee_angle=result_ig[i]["knee_ankle_angle"]
  303. if int(i)<=max_index/2 and int(i)>max_index/3:
  304. val=math.fabs( math.fabs(float(result_ig[i]["hip_knee_angle"]))-math.fabs(float(result_ig[i]["knee_ankle_angle"])))
  305. if(math.fabs( math.fabs(float(result_ig[i]["hip_knee_angle"]))-math.fabs(float(result_ig[i]["knee_ankle_angle"])))<less_angle):
  306. less_angle=val
  307. less_num=int(i)
  308. if int(i)<=max_index and int(i)>(max_index/3):
  309. if(float(result_ig[i]["hip_heigh"]))<(float(hip_hei_init)) and max_l_h_num!=1:
  310. break
  311. if (math.fabs(math.fabs(float(result_ig[i]["hip_knee_angle"]))-math.fabs(float(result_ig[i]["knee_ankle_angle"]))))>max_leg_knee_hip_knee and math.fabs(float(result_ig[i]["leg_heigh"])-float(ankle))<=20:
  312. max_leg_knee_hip_knee=math.fabs(float(result_ig[i]["hip_knee_angle"])-math.fabs(float(result_ig[i]["knee_ankle_angle"])))
  313. max_l_h_num=int(i)
  314. # ankle=result_ig[i]["leg_heigh"]
  315. # ankle_num=int(i)
  316. #elif ankle!=1 and int(i)>=(le/3*2):
  317. # break
  318. print(max_leg_knee_angle)
  319. print(max_leg_knee_angle_num)
  320. print(less_num)
  321. print(less_angle)
  322. #print(ankle)
  323. #print(ankle_num)
  324. print(max_leg_knee_hip_knee)
  325. print(max_l_h_num)
  326. re={}
  327. re["1"]={}
  328. re["1"]["num"]=max_leg_knee_angle_num
  329. re["1"]["val"]=max_leg_knee_angle
  330. re["2"]={}
  331. re["2"]["num"]=less_num
  332. re["2"]["val"]=less_angle
  333. re["3"]={}
  334. re["3"]["num"]=max_l_h_num
  335. re["3"]["val"]=max_leg_knee_hip_knee
  336. print(re)
  337. return result_ig,re
  338. def ignore_data_shendun(result):
  339. result_ig=result
  340. k="0"
  341. last={}
  342. max_index=1
  343. for re in result:
  344. if k=="0":
  345. k=re
  346. else:
  347. if(max_index<int(re)):
  348. max_index=int(re)
  349. max_d_num=k
  350. max_d=result[k]["r_d"]
  351. for i in result_ig:
  352. if(float(result_ig[i]["r_d"])>float(max_d)):
  353. max_d_num=int(i)
  354. max_d= result_ig[i]["r_d"]
  355. re={}
  356. re["1"]={}
  357. re["1"]["num"]=max_d_num
  358. re["1"]["val"]=max_d
  359. #1最大的深蹲角度
  360. return result_ig,re
  361. def write1(result,name):
  362. output=open(name,"w")
  363. for re in result:
  364. print(re,file=output)
  365. for(name,value) in result[re].items():
  366. print(name,value,sep=',',file=output)
  367. print("kv-over",file=output)
  368. output.close()
  369. return name
  370. def run(video_path,mode,need_split):
  371. path=video_path
  372. os.path.split(path)
  373. tur=os.path.split(path)
  374. if(mode=="pingban"):
  375. print(need_split)
  376. name=tur[len(tur)-1].split('.')[0]
  377. if(need_split==1):
  378. c=split_video(video_path)
  379. process_photo(name,c)
  380. gather_video(video_path)
  381. do_analysis("./capture_image/"+name+"/capture_image{}-1.png.npy",video_path,mode)
  382. switch = {'jump': ignore_data_jump, # 注意此处不要加括号
  383. 'pingban': ignore_data_pingban,
  384. 'juanfu': ignore_data_juanfu,
  385. 'shendun':ignore_data_shendun,
  386. 'gaotaitui':ignore_data_gaotaitui,
  387. }
  388. choice = mode # 获取选择
  389. result,re=switch.get(choice, ignore_data_jump)(read2("./cap_file/"+name+"/file1.txt"))
  390. # 执行对应的函数,如果没有就执行默认的函数
  391. #result,re=ignore_data_juanfu(read2("./cap_file/"+name+"/file1.txt"))
  392. write1(result,"./cap_file/"+name+"/file3.txt")
  393. write1(re,"./cap_file/"+name+"-file.txt")
  394. #write1()
  395. return re
  396. if __name__ == '__main__':
  397. video_path = "./video/shendun.mp4"
  398. #path=video_path
  399. #os.path.split(path)
  400. #tur=os.path.split(path)
  401. #name=tur[len(tur)-1].split('.')[0]
  402. #c=split_video(video_path)
  403. #process_photo(name,c)
  404. #gather_video(video_path)
  405. #do_analysis("./capture_image/"+name+"/capture_image{}-1.png.npy",video_path,"shendun")
  406. #result,re=ignore_data_shendun(read2("./cap_file/"+name+"/file1.txt"))
  407. #write1(result,"./cap_file/"+name+"/file3.txt")
  408. #write1(re,"./cap_file/"+name+"-file.txt")
  409. #write1()
  410. #print(re)
  411. re=run(video_path,"shendun",0)
  412. #参数一 video_path
  413. #参数二 shendun jump gaotaitui juanfu pingban
  414. #第三个参数为1 表示需要视频切割