wassname/scikit-image
1
0
Fork 0
mirror of https://github.com/wassname/scikit-image.git synced 2026-07-25 13:30:51 +08:00
Code Issues Packages Projects Releases Wiki Activity

iradon_sart: flip signs to reflect changes to radon.

Browse Source 
The changes to radon were introduced in gh-596.
This commit is contained in:
Jostein Bø Fløystad
2013-07-05 13:14:39 +02:00
parent af47941435
commit 3cdfeccc7a
1 changed files with 8 additions and 8 deletions
Download Patch File Download Diff File Expand all files Collapse all files
skimage/transform/_radon_transform.pyx
+8 -8
View File
@@ -50,11 +50,11 @@ cpdef bilinear_ray_sum(cnp.double_t[:, :] image, cnp.double_t theta,
if Ns > 0:
# step length between samples
ds = 2 * s0 / Ns
dx = ds * cos(theta)
dy = ds * sin(theta)
dx = -ds * cos(theta)
dy = -ds * sin(theta)
# point of entry of the ray into the reconstruction circle
x0 = -s0 * cos(theta) + t * sin(theta)
y0 = -s0 * sin(theta) - t * cos(theta)
x0 = s0 * cos(theta) - t * sin(theta)
y0 = s0 * sin(theta) + t * cos(theta)
for k in range(Ns+1):
x = x0 + k * dx
y = y0 + k * dy
@@ -134,11 +134,11 @@ cpdef bilinear_ray_update(cnp.double_t[:, :] image,
if Ns > 0:
# Step length between samples
ds = 2 * s0 / Ns
dx = ds * cos(theta)
dy = ds * sin(theta)
dx = -ds * cos(theta)
dy = -ds * sin(theta)
# Point of entry of the ray into the reconstruction circle
x0 = -s0 * cos(theta) + t * sin(theta)
y0 = -s0 * sin(theta) - t * cos(theta)
x0 = s0 * cos(theta) - t * sin(theta)
y0 = s0 * sin(theta) + t * cos(theta)
for k in range(Ns+1):
x = x0 + k * dx
y = y0 + k * dy