From 3cdfeccc7a80884676a3f364e406c4d7a47d1e98 Mon Sep 17 00:00:00 2001
From: =?UTF-8?q?Jostein=20B=C3=B8=20Fl=C3=B8ystad?=
 <jostein.floystad@gmail.com>
Date: Fri, 5 Jul 2013 13:12:34 +0200
Subject: [PATCH] iradon_sart: flip signs to reflect changes to radon.

The changes to radon were introduced in gh-596.
---
 skimage/transform/_radon_transform.pyx | 16 ++++++++--------
 1 file changed, 8 insertions(+), 8 deletions(-)

diff --git a/skimage/transform/_radon_transform.pyx b/skimage/transform/_radon_transform.pyx
index 3f865a93..91b943b9 100644
--- a/skimage/transform/_radon_transform.pyx
+++ b/skimage/transform/_radon_transform.pyx
@@ -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