iradon_sart: redefine projection center.

This reflects the bugfixes related to issue gh-592 implemented in
gh-596.
This commit is contained in:
Jostein Bø Fløystad
2013-07-05 13:14:39 +02:00
parent 06857d3b92
commit f2790d658c
+2 -2
View File
@@ -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