save backprojections

Author: Georg Schramm

python/02_reconstruct_petsird.py

@@ -325,6 +325,8 @@
 
 non_tof_backproj = proj.adjoint(xp.ones(coords0.shape[0], dtype="float32"))
 
+xp.save(Path(fname).with_suffix(".non_tof_backproj.npy"), non_tof_backproj)
+
 #### HACK assumes same TOF parameters for all module type pairs
 sigma_tof = scanner_info.tof_resolution[0][0] / 2.35
 tof_bin_edges = scanner_info.tof_bin_edges[0][0].edges
@@ -341,6 +343,8 @@
 
 tof_backproj = proj.adjoint(xp.ones(coords0.shape[0], dtype="float32"))
 
+xp.save(Path(fname).with_suffix(".tof_backproj.npy"), tof_backproj)
+
 del proj
 
 # %%
@@ -402,6 +406,7 @@
             end="\r",
         )
         lm_exp = effs[sl] * lm_subset_projs[i_subset](res_model(recon))
+
         tmp = num_subsets * res_model.adjoint(
             lm_subset_projs[i_subset].adjoint(effs[sl] / lm_exp)
         )