From a4870242b63458ed9ee79db6a951007c94deebf4 Mon Sep 17 00:00:00 2001 From: =?UTF-8?q?Jostein=20B=C3=B8=20Fl=C3=B8ystad?= Date: Mon, 24 Jun 2013 10:18:54 +0200 Subject: [PATCH] iradon_sart: Test clip functionality and add a test for it. --- skimage/transform/radon_transform.py | 2 +- skimage/transform/tests/test_radon_transform.py | 3 +++ 2 files changed, 4 insertions(+), 1 deletion(-) diff --git a/skimage/transform/radon_transform.py b/skimage/transform/radon_transform.py index 0f26e343..ff3af596 100644 --- a/skimage/transform/radon_transform.py +++ b/skimage/transform/radon_transform.py @@ -413,5 +413,5 @@ def iradon_sart(radon_image, theta=None, image=None, projection_shifts=None, projection_shifts[angle_index]) image += relaxation * image_update if not clip is None: - image = clip(image, clip[0], clip[1]) + image = np.clip(image, clip[0], clip[1]) return image diff --git a/skimage/transform/tests/test_radon_transform.py b/skimage/transform/tests/test_radon_transform.py index 74d20a24..11fced56 100644 --- a/skimage/transform/tests/test_radon_transform.py +++ b/skimage/transform/tests/test_radon_transform.py @@ -347,6 +347,9 @@ def test_iradon_sart(): delta = np.mean(np.abs(reconstructed - image)) print('delta (2 iterations) =', delta) assert delta < 0.013 * error_factor + reconstructed = iradon_sart(sinogram, theta, clip=(0, 1)) + print('delta (1 iteration, clip) =', delta) + assert delta < 0.013 * error_factor np.random.seed(1239867) shifts = np.random.uniform(-3, 3, sinogram.shape[1])