corrected frwd projector output and added sinogram segment definition…

… to invaux

niftypet/nipet/invaux.py

@@ -100,7 +100,26 @@ def axial_lut(Cnt):
         z += Cnt['AXR']
         rng[i, 1] = z
 
-
+    # > sinogram segments
+    if Cnt['SPN']>1:
+        SPN = Cnt['SPN']
+        MNRD = [-(SPN//2),]
+        MXRD = [SPN//2,]
+        for i in range(SPN//2+1, NRNG, SPN):
+            s = 2*(NRNG-i)-1
+
+            if i>Cnt['MRD']: break
+
+            MNRD.append(-i-(SPN-1))
+            MXRD.append(-i)
+
+            MNRD.append(i)
+            MXRD.append(i+(SPN-1))
+
+            #print(2*s, i, (i+2))
+
+        Cnt['MNRD'] = np.array(MNRD)
+        Cnt['MXRD'] = np.array(MXRD)
 
     #---------------------------------------------------------------------
     # > Michelogram for single slice rebinning

niftypet/nipet/prj_inv/prj_inv.py

@@ -123,7 +123,8 @@ def frwd_prj(im, scanner_params, isub=ISUB_DEFAULT, dev_out=False, attenuation=F
 
     # put the gaps back to form displayable sinogram
     if not dev_out and fullsino_out:
-        sino = mmraux.putgaps(sino, txLUT, Cnt)
+        sino = np.reshape(sino, (Cnt['NSBINS'], Cnt['NSANGLES'], nsinos))
+        sino = np.transpose(sino, (2,1,0))
 
     return sino