added zoom to mMR recon

Author: pjmark

niftypet/nipet/__init__.py

@@ -79,8 +79,8 @@
 from . import sigaux
 from .prj_sig import prjsig
 
-from . import lm_sig
-from .lm_sig.hst_sig import lminfo_sig
+#from . import lm_sig
+#from .lm_sig.hst_sig import lminfo_sig
 
 # > SynchroPET
 from . import synaux

niftypet/nipet/include/scanner_mmr.h

@@ -117,6 +117,8 @@ struct Cnst {
   // float ICOSSTP;
 
   float ETHRLD;
+
+  float ZOOM; // zoom for changing the voxel size and dimensions
 };
 //<><><><><><><><><><><><><><><><><><><><><><><><><><><><><><><><><><><><><><><>
 

niftypet/nipet/prj/src/prj_module.cu

@@ -276,6 +276,8 @@ static PyObject *frwd_prj(PyObject *self, PyObject *args, PyObject *kwargs) {
   Cnt.LOG = (char)PyLong_AsLong(pd_log);
   PyObject *pd_devid = PyDict_GetItemString(o_mmrcnst, "DEVID");
   Cnt.DEVID = (char)PyLong_AsLong(pd_devid);
+  PyObject *pd_zoom = PyDict_GetItemString(o_mmrcnst, "zoom");
+  Cnt.ZOOM = (float)PyFloat_AsDouble(pd_zoom);
 
   /* Interpret the input objects as numpy arrays. */
   // axial LUTs:
@@ -471,6 +473,8 @@ static PyObject *back_prj(PyObject *self, PyObject *args, PyObject *kwargs) {
   Cnt.LOG = (char)PyLong_AsLong(pd_log);
   PyObject *pd_devid = PyDict_GetItemString(o_mmrcnst, "DEVID");
   Cnt.DEVID = (char)PyLong_AsLong(pd_devid);
+  PyObject *pd_zoom = PyDict_GetItemString(o_mmrcnst, "zoom");
+  Cnt.ZOOM = (float)PyFloat_AsDouble(pd_zoom);
 
   /* Interpret the input objects as numpy arrays. */
   // axial LUTs:

niftypet/nipet/prj/src/prjb.cu

@@ -221,10 +221,12 @@ void gpu_bprj(float *d_im, float *d_sino, float *li2rng, short *li2sn, char *li2
     nrng_c = NRINGS;
     nil2r_c = NLI2R;
   }
-  // voxels in axial direction
-  vz0 = 2 * Cnt.RNG_STRT;
-  vz1 = 2 * (Cnt.RNG_END - 1);
-  nvz = 2 * nrng_c - 1;
+
+  // voxels in axial direction (to accounting for the variable AFOV)
+  vz0 = (int)Cnt.ZOOM * (2 * Cnt.RNG_STRT);
+  vz1 = (int)Cnt.ZOOM * (2 * (Cnt.RNG_END - 1));
+  nvz = (int)Cnt.ZOOM * (2 * nrng_c - 1);
+
   if (Cnt.LOG <= LOGDEBUG) {
     printf("i> detector rings range: [%d, %d) => number of  sinos: %d\n", Cnt.RNG_STRT,
            Cnt.RNG_END, snno);

niftypet/nipet/prj/src/prjf.cu

@@ -238,10 +238,12 @@ void gpu_fprj(float *d_sn, float *d_im, float *li2rng, short *li2sn, char *li2no
     nrng_c = NRINGS;
     nil2r_c = NLI2R;
   }
+  
   // voxels in axial direction
-  vz0 = 2 * Cnt.RNG_STRT;
-  vz1 = 2 * (Cnt.RNG_END - 1);
-  nvz = 2 * nrng_c - 1;
+  vz0 = (int)Cnt.ZOOM * (2 * Cnt.RNG_STRT);
+  vz1 = (int)Cnt.ZOOM * (2 * (Cnt.RNG_END - 1));
+  nvz = (int)Cnt.ZOOM * (2 * nrng_c - 1);
+
   if (Cnt.LOG <= LOGDEBUG) {
     printf("i> detector rings range: [%d, %d) => number of  sinos: %d\n", Cnt.RNG_STRT,
            Cnt.RNG_END, snno);