From 006adb3112cbfad7de7da0611861f6c4d9d70fc5 Mon Sep 17 00:00:00 2001
From: 2wenty1ne <sirdarkvic@gmail.com>
Date: Sun, 3 Nov 2024 18:53:40 +0100
Subject: [PATCH] works, but weak performance seperating all brightness values
 from the lung

---
 main.py | 18 ++++++++++++++++--
 1 file changed, 16 insertions(+), 2 deletions(-)

diff --git a/main.py b/main.py
index 6f75a36..f6a3c09 100644
--- a/main.py
+++ b/main.py
@@ -35,11 +35,25 @@ lungsInLitre = voxelInLungs / 1000000
 # Histogramm des CT-Bilds
 thoraxArray = np.array(sitk.GetArrayFromImage(img))
 thoraxArray1D = thoraxArray.flatten()
-plt.hist(thoraxArray1D, 100)
-plt.title("Histogramm eines Thorax")
+plt.hist(thoraxArray1D, bins=100)
+plt.title("Helligkeitswerte eines Thorax")
 plt.xlabel("Helligkeitswert in Hounsfield-Einheiten (HU)")
 plt.ylabel("Häufigkeit in Millionen (mio)")
 plt.show()
 plt.clf()
 
+# Alle Helligkeitswerte der Lunge seperieren
+lungBrightness = []
+lungsArray1D = lungsArray.flatten()
 
+for (ind, elem) in enumerate(lungsArray1D):
+    if elem > 0:
+        lungBrightness.append(thoraxArray1D[ind])
+
+# Histogram aller Helligkeitswerte in der Lunge 
+plt.hist(lungBrightness, bins=100)
+plt.title("Helligkeitswerte einer Lunge")
+plt.xlabel("Helligkeitswert in Hounsfield-Einheiten (HU)")
+plt.ylabel("Häufigkeit")
+plt.show()
+plt.clf()