通過python讀取txt文件和繪制柱形圖的實(shí)現(xiàn)代碼
目的
臨床數(shù)據(jù)的記錄時(shí)間和對(duì)應(yīng)標(biāo)簽(逗號(hào)后面的數(shù)字)記錄在txt文件里,要把標(biāo)簽轉(zhuǎn)換為3類標(biāo)簽,并且計(jì)算出每個(gè)標(biāo)簽的分別持續(xù)時(shí)間,然后繪制成柱形圖方便查閱。
小難點(diǎn)分析:
(1)txt的切割讀取對(duì)應(yīng)內(nèi)容
(2)時(shí)間差計(jì)算
txt文件如圖:

使用效果
首先將原始txt轉(zhuǎn)換為 左列新標(biāo)簽 右列持續(xù)時(shí)間

繪制為柱形圖

為了直觀,每次只最多顯示 2個(gè)小時(shí),同時(shí)橫坐標(biāo)還是按照臨床的記錄時(shí)間順序。
代碼實(shí)現(xiàn)
# -*- coding: utf-8 -*-
from datetime import datetime
import matplotlib.pyplot as plt
from matplotlib.font_manager import FontProperties
font_set = FontProperties(fname=r"c:\windows\fonts\simsun.ttc", size=12)
address = "D:/XXXXXX/"
patient = "XXXXX"
year = '2020' # or 2020
txt_address = address + patient + "/timetable.txt"
filename = address + patient + "/newTimeTable.txt" # 保存的txt地址
# 讀取txt的相應(yīng)內(nèi)容
def showLabel():
content = []
label = []
rawtime = []
with open(txt_address, 'r', encoding="utf-8") as f:
lines = f.readlines() # 返回的是列表
print("文件的行數(shù):%d" % len(lines))
for line in lines:
content_value = line.split(sep=',') # 分為[0]和[1]
time_value = content_value[0].split(sep=year)
content_value = content_value[1].split(sep='\n')
temp_time = year + time_value[-1]
rawtime.append(temp_time)
content.append(content_value[0])
# 去掉第一行的line和轉(zhuǎn)換label
for i in range(1, len(lines)): # range是左閉右開
number = int(content[i])
if number == 4:
label_value = 1 # label 1 表示癲癇發(fā)作
elif number == 1:
label_value = 0 # label 0 表示數(shù)據(jù)缺失
else:
label_value = -1 # label -1 表示非癲癇發(fā)作
label.append(label_value)
del rawtime[0] # 去掉第一行的line
length = len(lines)
return label, rawtime, length
# 計(jì)算時(shí)間差
def getTimeSpan(new_label,time,length):
time_span = []
time_label = []
temp = datetime.strptime(time[0], "%Y-%m-%d %H:%M:%S")
for i in range(0, length - 2):
if new_label[i] != new_label[i + 1]:
end_time = datetime.strptime(time[i + 1], "%Y-%m-%d %H:%M:%S")
diff = end_time - temp
temp = end_time
# time_span.append(diff) 標(biāo)準(zhǔn)化時(shí)間
time_span.append(round(diff.total_seconds() / 60, 2)) # 保留兩位小數(shù)
time_label.append(new_label[i])
return time_span, time_label
# 畫圖,按照標(biāo)簽內(nèi)容給圖上色
def showTimeandLabel(time_span,time_label):
Yanse = []
for i in time_label:
if i == 0:
yanse = 'blue'
elif i == -1:
yanse = 'gray'
else:
yanse = 'red'
Yanse.append(yanse)
plt.bar(range(len(time_span)), time_span, width=0.8, tick_label=time_label, color=Yanse)
plt.xlabel('標(biāo)簽', fontproperties=font_set)
plt.ylabel("持續(xù)時(shí)間(分鐘)", fontproperties=font_set)
plt.title("patient: " + patient, size=20)
plt.ylim(0, 120)
plt.show()
# 保存到txt
def txtSave(time_span, time_label):
# w:向文件中寫入內(nèi)容時(shí),會(huì)先清空原文件中的內(nèi)容,
with open(filename, 'w', encoding="utf-8") as f:
line_number = len(time_label)
for i in range(line_number):
f.write(str(time_label[i])+" "+str(time_span[i]))
f.write("\n")
def main():
new_label, time, length = showLabel()
time_span, time_label = getTimeSpan(new_label, time, length)
txtSave(time_span, time_label)
showTimeandLabel(time_span, time_label)
if __name__ == '__main__':
main()
到此這篇關(guān)于通過python讀取txt文件和繪制柱形圖的文章就介紹到這了,更多相關(guān)python繪制柱形圖內(nèi)容請(qǐng)搜索腳本之家以前的文章或繼續(xù)瀏覽下面的相關(guān)文章希望大家以后多多支持腳本之家!
相關(guān)文章
pytorch教程resnet.py的實(shí)現(xiàn)文件源碼分析
torchvision.models這個(gè)包中包含alexnet、densenet、inception、resnet、squeezenet、vgg等常用的網(wǎng)絡(luò)結(jié)構(gòu),并且提供了預(yù)訓(xùn)練模型,可以通過簡(jiǎn)單調(diào)用來讀取網(wǎng)絡(luò)結(jié)構(gòu)和預(yù)訓(xùn)練模型2021-09-09
Python pytest裝飾器總結(jié)(實(shí)例詳解)
這篇文章主要介紹了Python pytest裝飾器總結(jié),本文給大家介紹的非常詳細(xì),對(duì)大家的學(xué)習(xí)或工作具有一定的參考借鑒價(jià)值,需要的朋友可以參考下2021-04-04
django如何根據(jù)現(xiàn)有數(shù)據(jù)庫(kù)表生成model詳解
這篇文章主要給大家介紹了關(guān)于django如何根據(jù)現(xiàn)有數(shù)據(jù)庫(kù)表生成model的相關(guān)資料,文中通過示例代碼介紹的非常詳細(xì),對(duì)大家學(xué)習(xí)或者使用Django具有一定的參考學(xué)習(xí)價(jià)值,需要的朋友可以參考下2022-08-08
PyCharm無(wú)法識(shí)別PyQt5的2種解決方法,ModuleNotFoundError: No module named
這篇文章主要介紹了PyCharm無(wú)法識(shí)別PyQt5的兩種解決辦法,ModuleNotFoundError: No module named 'pyqt5',需要的朋友可以參考下2020-02-02
使用python進(jìn)行廣告點(diǎn)擊率的預(yù)測(cè)的實(shí)現(xiàn)
這篇文章主要介紹了使用python進(jìn)行廣告點(diǎn)擊率的預(yù)測(cè)的實(shí)現(xiàn),文中通過示例代碼介紹的非常詳細(xì),對(duì)大家的學(xué)習(xí)或者工作具有一定的參考學(xué)習(xí)價(jià)值,需要的朋友們下面隨著小編來一起學(xué)習(xí)學(xué)習(xí)吧2019-07-07

