Merge pull request #211 from DiamondLightSource/fbp2dcorr

corrections to the fbp2d astra to flipud the result

Author: dkazanc

httomolibgpu/recon/algorithm.py

@@ -119,12 +119,14 @@ def FBP2d_astra(
 
     # loop over detY slices
     for slice_index in range(0, detY_size):
-        reconstruction[:, slice_index, :] = RecTools.FBP(
-            data[:, slice_index, :],
-            filter_type=filter_type,
-            filter_parameter=filter_parameter,
-            filter_d=filter_d,
-            recon_mask_radius=recon_mask_radius,
+        reconstruction[:, slice_index, :] = np.flipud(
+            RecTools.FBP(
+                data[:, slice_index, :],
+                filter_type=filter_type,
+                filter_parameter=filter_parameter,
+                filter_d=filter_d,
+                recon_mask_radius=recon_mask_radius,
+            )
         )
     return reconstruction
 

zenodo-tests/test_recon/test_algorithm.py

@@ -41,7 +41,7 @@ def test_reconstruct_FBP2d_astra_i12_dataset1(i12_dataset1):
         recon_mask_radius=0.9,
     )
     assert recon_data.flags.c_contiguous
-    assert_allclose(np.sum(recon_data), 84673.68, rtol=1e-07, atol=1e-6)
+    assert_allclose(np.sum(recon_data), 84672.84, atol=1e-2)
     assert recon_data.dtype == np.float32
     assert recon_data.shape == (2560, 50, 2560)