From 28de2f978abf5fa7f9a1c4d71548d36325f1da44 Mon Sep 17 00:00:00 2001 From: =?UTF-8?q?Jostein=20B=C3=B8=20Fl=C3=B8ystad?= Date: Tue, 18 Jun 2013 18:23:09 +0200 Subject: [PATCH] transform.radon: Remove unneccesary matrix inverse. --- 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 59e6c338..e3f46f30 100644 --- a/skimage/transform/radon_transform.py +++ b/skimage/transform/radon_transform.py @@ -98,17 +98,16 @@ def radon(image, theta=None, circle=False): [0, 0, 1]]) def build_rotation(theta): - T = -np.deg2rad(theta) + T = np.deg2rad(theta) - R = np.array([[np.cos(T), -np.sin(T), 0], - [np.sin(T), np.cos(T), 0], + R = np.array([[np.cos(T), np.sin(T), 0], + [-np.sin(T), np.cos(T), 0], [0, 0, 1]]) return shift1.dot(R).dot(shift0) for i in range(len(theta)): - rotated = _warp_fast(padded_image, - np.linalg.inv(build_rotation(-theta[i]))) + rotated = _warp_fast(padded_image, build_rotation(theta[i])) out[:, i] = rotated.sum(0)[::-1]