From 060a4df5c2993b9a0a4d5b1dd2614b71c14014ef Mon Sep 17 00:00:00 2001
From: 3011357 <3011357@stud.hs-mannheim.de>
Date: Sat, 2 Nov 2024 22:35:31 +0100
Subject: [PATCH] Histogramm des CT-Bilds erstellt + beschriftet

---
 main.py | 14 ++++++++++++++
 1 file changed, 14 insertions(+)

diff --git a/main.py b/main.py
index a7a0f15..6f75a36 100644
--- a/main.py
+++ b/main.py
@@ -1,5 +1,6 @@
 import SimpleITK as sitk
 import numpy as np
+import matplotlib.pyplot as plt
 
 img = sitk.ReadImage("Thorax1.0B40f_Downsampled_by_2.nrrd")
 
@@ -20,12 +21,25 @@ print(f"3) Die Methode GetSize() gibt den Typ tuple zurück. \n")
 #Sollte man kontrollieren
 print(f"4) Der Bildausschnitt in die drei Raumrichtungen ist X: {xSp*xSi} mm, Y: {ySp*xSi} mm, Z: {zSp*xSi} mm\n")
 
+# CT-Bild der Lunge erstellen
 seedList = [(120,140,100)]
 lungs = sitk.ConnectedThreshold(img, seedList , -1100, -200)
 sitk.WriteImage(lungs, "lungs.nrrd")
 
+# Volumen der Lunge berechnen
 lungsArray = np.array(sitk.GetArrayFromImage(lungs))
 voxelInLungs = lungsArray[lungsArray > 0].size * volumePerVoxel
 lungsInLitre = voxelInLungs / 1000000
 #print(lungsInLitre)
 
+# Histogramm des CT-Bilds
+thoraxArray = np.array(sitk.GetArrayFromImage(img))
+thoraxArray1D = thoraxArray.flatten()
+plt.hist(thoraxArray1D, 100)
+plt.title("Histogramm eines Thorax")
+plt.xlabel("Helligkeitswert in Hounsfield-Einheiten (HU)")
+plt.ylabel("Häufigkeit in Millionen (mio)")
+plt.show()
+plt.clf()
+
+