From f2790d658c0e6f94bc1a479c4eed162907bf5a05 Mon Sep 17 00:00:00 2001 From: =?UTF-8?q?Jostein=20B=C3=B8=20Fl=C3=B8ystad?= Date: Fri, 5 Jul 2013 13:02:03 +0200 Subject: [PATCH] iradon_sart: redefine projection center. This reflects the bugfixes related to issue gh-592 implemented in gh-596. --- skimage/transform/_radon_transform.pyx | 4 ++-- 1 file changed, 2 insertions(+), 2 deletions(-) diff --git a/skimage/transform/_radon_transform.pyx b/skimage/transform/_radon_transform.pyx index d6112a57..3f865a93 100644 --- a/skimage/transform/_radon_transform.pyx +++ b/skimage/transform/_radon_transform.pyx @@ -33,7 +33,7 @@ cpdef bilinear_ray_sum(cnp.double_t[:, :] image, cnp.double_t theta, """ theta = theta / 180. * M_PI cdef cnp.double_t radius = image.shape[0] // 2 - 1 - cdef cnp.double_t projection_center = image.shape[0] // 2 - 1 + cdef cnp.double_t projection_center = image.shape[0] // 2 cdef cnp.double_t rotation_center = image.shape[0] // 2 # (s, t) is the (x, y) system rotated by theta cdef cnp.double_t t = ray_position - projection_center @@ -118,7 +118,7 @@ cpdef bilinear_ray_update(cnp.double_t[:, :] image, deviation = 0. theta = theta / 180. * M_PI cdef cnp.double_t radius = image.shape[0] // 2 - 1 - cdef cnp.double_t projection_center = image.shape[0] // 2 - 1 + cdef cnp.double_t projection_center = image.shape[0] // 2 cdef cnp.double_t rotation_center = image.shape[0] // 2 # (s, t) is the (x, y) system rotated by theta cdef cnp.double_t t = ray_position - projection_center