From ac460a9777303df3d4df17a06c1ec70b53122f68 Mon Sep 17 00:00:00 2001 From: =?UTF-8?q?Jostein=20B=C3=B8=20Fl=C3=B8ystad?= Date: Fri, 5 Jul 2013 14:13:41 +0200 Subject: [PATCH] iradon: Use util.pad for array padding. This depends on util.pad accepting 0s in pad_width, which has a proposed solution in gh-634. --- skimage/transform/radon_transform.py | 9 ++++----- 1 file changed, 4 insertions(+), 5 deletions(-) diff --git a/skimage/transform/radon_transform.py b/skimage/transform/radon_transform.py index dcb13cae..93979977 100644 --- a/skimage/transform/radon_transform.py +++ b/skimage/transform/radon_transform.py @@ -109,11 +109,10 @@ def radon(image, theta=None, circle=False): def _sinogram_circle_to_square(sinogram): - size = int(np.ceil(np.sqrt(2) * sinogram.shape[0])) - sinogram_padded = np.zeros((size, sinogram.shape[1])) - pad = (size - sinogram.shape[0]) // 2 - sinogram_padded[pad:pad + sinogram.shape[0], :] = sinogram - return sinogram_padded + diagonal = int(np.ceil(np.sqrt(2) * sinogram.shape[0])) + pad = diagonal - sinogram.shape[0] + pad_width = ((pad // 2, pad - pad // 2), (0, 0)) + return util.pad(sinogram, pad_width, mode='constant', constant_values=0) def iradon(radon_image, theta=None, output_size=None,