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()