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

Added DBP,MBP,SBP,PP at each outlet

Browse Source
This commit is contained in:
Aloma Blanch 2020-07-16 20:30:27 -05:00
parent 0a949ffd3d
commit 2c64d160a2
2 changed files with 19 additions and 3 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

20
functions.py
View File

@ -7,6 +7,7 @@ from tkinter import Tk
from tkinter.filedialog import askopenfilename, asksaveasfile, askdirectory from tkinter.filedialog import askopenfilename, asksaveasfile, askdirectory
import pandas as pd import pandas as pd
import tkinter as tk import tkinter as tk
from statistics import mean
from scipy.signal import find_peaks from scipy.signal import find_peaks
@ -81,15 +82,30 @@ def pressure(folder,N_ts,T_cyc,dt,n_cyc):
Nc = round(T_cyc/dt) Nc = round(T_cyc/dt)
time = np.linspace(0,T_cyc,Nc) time = np.linspace(0,T_cyc,Nc)
fig, ax = plt.subplots() fig, ax = plt.subplots()
# SBP = []
# DBP = []
# MBP = []
SBP = np.empty(pressure.shape[1])
DBP = np.empty(pressure.shape[1])
MBP = np.empty(pressure.shape[1])
for i in range(0,pressure.shape[1]): for i in range(0,pressure.shape[1]):
ax.plot(time,pressure[N_ts-Nc:N_ts,-i]/1333.22,label='ROI-'+str(i+2)) ax.plot(time,pressure[N_ts-Nc:N_ts,i]/1333.22,label='ROI-'+str(i+2))
# SBP.append(np.amax(pressure[N_ts-Nc:N_ts,i]/1333.22))
# DBP.append(np.amin(pressure[N_ts-Nc:N_ts,i]/1333.22))
# MBP.append(mean(pressure[N_ts-Nc:N_ts,i]/1333.22))
SBP[i] = (np.amax(pressure[N_ts-Nc:N_ts,i]/1333.22))
DBP[i] = (np.amin(pressure[N_ts-Nc:N_ts,i]/1333.22))
MBP[i] = (mean(pressure[N_ts-Nc:N_ts,i]/1333.22))
PP = SBP-DBP
ax.set(xlabel='time [s]', ylabel='Pressure [mmHg]', ax.set(xlabel='time [s]', ylabel='Pressure [mmHg]',
title='Pressure @ each outlet') title='Pressure @ each 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)
plt.show() plt.show()
return (DBP,MBP,SBP,PP)

2
main.py
View File

@ -49,5 +49,5 @@ error_plot(folder,dt,rc,False)
periodicity(project,folder,dt,T_cyc,n_cyc) periodicity(project,folder,dt,T_cyc,n_cyc)
# Pressure # Pressure
pressure(folder,N_ts,T_cyc,dt,n_cyc) (DBP,MBP,SBP,PP) = pressure(folder,N_ts,T_cyc,dt,n_cyc)