aloma/SimVascular_report_PostProcess
1
0
Fork 0
Code Issues Pull Requests Releases Wiki Activity

Changed aimed.csv imput file for aimed.txt

Browse Source
This commit is contained in:
Aloma Blanch 2020-07-24 08:22:33 -05:00
parent 070428663c
commit c3c570764b
1 changed files with 7 additions and 7 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

14
functions.py
View File

@ -42,7 +42,7 @@ def cycle(folder,dt,save_path):
ax.plot(time,pressure[:,-1]/1333.22,'b') ax.plot(time,pressure[:,-1]/1333.22,'b')
ax.plot(t_pk_loc, P_peaks/1333.22, "or",label='peak location') ax.plot(t_pk_loc, P_peaks/1333.22, "or",label='peak location')
ax.set(xlabel='time [s]', ylabel='Pressure [mmHg]', ax.set(xlabel='time [s]', ylabel='Pressure [mmHg]',
title='Pressure at last outlet') title='Pressure at last outlet')
ax.spines['right'].set_visible(False) ax.spines['right'].set_visible(False)
ax.spines['top'].set_visible(False) ax.spines['top'].set_visible(False)
ax.legend(loc=0) ax.legend(loc=0)
@ -221,9 +221,11 @@ def rcr(project_folder):
return Rc_C_Rd return Rc_C_Rd
def barPlot(project_folder,DBP,MBP,SBP,PP,Q_avg,save_path): def barPlot(project_folder,DBP,MBP,SBP,PP,Q_avg,save_path):
import csv aimed_all = [1] #No li agrada fer append en una llista buida, no recordo perquè i segur que es pot fer millor
with open(project_folder+'/aimed.csv', 'r') as f: f = open(project_folder+'/aimed.txt', 'r') #obrir el fitxer
aimed_all = list(csv.reader(f, delimiter=';')) for line in f:
aimed_all = line.split() #fas una llista amb cada valor de la línia carregan-te els espais buits
del aimed_all[0] #esborrar el primer valor que has posat perquè no estigues buida la llista
def plot_bar(CFD,aimed,name,save_path,unit): def plot_bar(CFD,aimed,name,save_path,unit):
labels = [] labels = []
@ -292,6 +294,4 @@ def barPlot(project_folder,DBP,MBP,SBP,PP,Q_avg,save_path):
# DBP # DBP
CFD = [round(float(i),1) for i in Q_avg] CFD = [round(float(i),1) for i in Q_avg]
aimed = [round(float(i),1) for i in aimed_all[4]] aimed = [round(float(i),1) for i in aimed_all[4]]
plot_bar(CFD,aimed,'Q_avg',save_path,' [mL/s]') plot_bar(CFD,aimed,'Q_avg',save_path,' [mL/s]')