Merge pull request #237 from NicoleJurjew/PSMR_24_PET_introductory

Bugfix in MLEM function of DIY_OSEM notebook in PET
SyneRBI · May 19, 2024 · 057e6b1 · 057e6b1
2 parents 7a97fde + ac04f57
commit 057e6b1
Showing 1 changed file with 2 additions and 0 deletions.
diff --git a/notebooks/PET/Solution_Snippets/DIY_OSEM_01_MLEM.py b/notebooks/PET/Solution_Snippets/DIY_OSEM_01_MLEM.py
@@ -6,4 +6,6 @@ def MLEM(acquired_data, acq_model, initial_image, num_iterations):
         quotient.fill(numpy.nan_to_num(quotient.as_array()))
         mult_update = acq_model.backward(quotient)/sensitivity         # A^t * quotient / A^t1
         mult_update.fill(numpy.nan_to_num(mult_update.as_array()))
+        estimated_image *= mult_update
+        estimated_image.maximum(0)
     return estimated_image