Merge branch 'brendt/fbp' into brendt/de-astra-fy-extended

lanl · Oct 9, 2024 · 36e68c7 · 36e68c7
2 parents 69a0525 + 5d37894
commit 36e68c7
Show file tree

Hide file tree

Showing 2 changed files with 5 additions and 4 deletions.
diff --git a/scico/linop/xray/_xray.py b/scico/linop/xray/_xray.py
@@ -125,13 +125,14 @@ def back_project(self, y: ArrayLike) -> snp.Array:
         return XRayTransform2D._back_project(y, self.x0, self.dx, self.nx, self.y0, self.angles)
 
     def fbp(self, y: ArrayLike) -> snp.Array:
-        """Compute filtered back projection (FBP) inverse of projection.
+        r"""Compute filtered back projection (FBP) inverse of projection.
 
         Compute the filtered back projection inverse by filtering each
         row of the sinogram with the filter defined in (61) in
         :cite:`kak-1988-principles` and then back projecting. The
-        projection angles are assumed to be evenly spaced: poor results
-        may be obtained if this assumption is violated.
+        projection angles are assumed to be evenly spaced in
+        :math:`[0, \pi)`; reconstruction quality may be poor if
+        this assumption is violated.
 
         Args:
             y: Input projection, (num_angles, N).

diff --git a/scico/test/linop/xray/test_xray.py b/scico/test/linop/xray/test_xray.py
@@ -81,7 +81,7 @@ def test_fbp(dx, det_count_factor):
 
     det_count = int(det_count_factor * N)
     n_proj = 360
-    angles = np.linspace(0, np.pi, n_proj)
+    angles = np.linspace(0, np.pi, n_proj, endpoint=False)
     A = XRayTransform2D(x_gt.shape, angles, det_count=det_count, dx=dx)
     y = A(x_gt)
     x_fbp = A.fbp(y)