From b3746b90902dbb47d5da2d8500973c6f00392ea9 Mon Sep 17 00:00:00 2001 From: =?UTF-8?q?Jostein=20B=C3=B8=20Fl=C3=B8ystad?= Date: Sat, 6 Jul 2013 18:29:31 +0200 Subject: [PATCH] iradon: Cleanup by locating related code in one place. --- skimage/transform/radon_transform.py | 5 ++--- 1 file changed, 2 insertions(+), 3 deletions(-) diff --git a/skimage/transform/radon_transform.py b/skimage/transform/radon_transform.py index 2f2a18e7..dfac0748 100644 --- a/skimage/transform/radon_transform.py +++ b/skimage/transform/radon_transform.py @@ -226,9 +226,6 @@ def iradon(radon_image, theta=None, output_size=None, [X, Y] = np.mgrid[0.0:x, 0.0:y] xpr = X - int(output_size) // 2 ypr = Y - int(output_size) // 2 - if circle: - radius = (output_size - 1) // 2 - reconstruction_circle = (xpr**2 + ypr**2) < radius**2 # Reconstruct image by interpolation for i in range(len(theta)): @@ -243,6 +240,8 @@ def iradon(radon_image, theta=None, output_size=None, backprojected = interpolant(t) reconstructed += backprojected if circle: + radius = (output_size - 1) // 2 + reconstruction_circle = (xpr**2 + ypr**2) < radius**2 reconstructed[~reconstruction_circle] = 0. return reconstructed * np.pi / (2 * len(th))