Commit 9b7565aa authored by Radek Galabov's avatar Radek Galabov
Browse files

Added script used for flow integration along time. Not important for the user.

parent 7ceb499d
Loading
Loading
Loading
Loading
+27 −0
Original line number Diff line number Diff line
import glob
import os
import numpy as np
import matplotlib.pyplot as plt


folder_path = r"C:\Users\galr\Desktop\pacienti_rozpracovani\G318_patient62_bodyPart2_visit0_autom_manual\my_flow" 

files = []
for filename in glob.glob(folder_path + os.sep + "**" + ".csv", recursive=True):
    files.append(np.genfromtxt(filename, delimiter=","))

flux_int = []
dist_int = []
for file in files:
    print("Time:",file[1][0])
    print("Flux_int",np.transpose(file)[:][4])
    flux_values = np.transpose(file)[:][4]
    dist_values = np.transpose(file)[:][3]
    print("Suma:", np.sum(flux_values[1:]))
    flux_int.append(flux_values)
    dist_int.append(dist_values)
    print("--------------")

print(np.sum(flux_int,0))
plt.plot(dist_int[0],np.sum(flux_int,0))
plt.show()