From 57b2b895c60eae69eb3b47e4dffcf4a2d67dbb3d Mon Sep 17 00:00:00 2001 From: William Joye Date: Fri, 3 Mar 2017 12:51:34 -0500 Subject: fix position angle for linear WCS --- tksao/frame/frmap.C | 32 ++++++++------------------------ 1 file changed, 8 insertions(+), 24 deletions(-) diff --git a/tksao/frame/frmap.C b/tksao/frame/frmap.C index feb267a..42c467e 100644 --- a/tksao/frame/frmap.C +++ b/tksao/frame/frmap.C @@ -6,7 +6,8 @@ #include "framebase.h" #include "fitsimage.h" -double Base::mapAngleFromRef(double angle, Coord::CoordSystem sys, Coord::SkyFrame sky) +double Base::mapAngleFromRef(double angle, Coord::CoordSystem sys, + Coord::SkyFrame sky) { double rr = angle; FitsImage* ptr = currentContext->cfits; @@ -20,23 +21,15 @@ double Base::mapAngleFromRef(double angle, Coord::CoordSystem sys, Coord::SkyFra case Coord::AMPLIFIER: break; default: - switch (ptr->getWCSOrientation(sys,sky)) { - case Coord::NORMAL: - rr += ptr->getWCSRotation(sys,sky); - break; - case Coord::XX: - rr = -(angle + ptr->getWCSRotation(sys,sky) + M_PI); - break; - case Coord::YY: - case Coord::XY: - break; - } + rr += ptr->getWCSRotation(sys,sky); + break; } return zeroTWOPI(rr); } -double Base::mapAngleToRef(double angle, Coord::CoordSystem sys, Coord::SkyFrame sky) +double Base::mapAngleToRef(double angle, Coord::CoordSystem sys, + Coord::SkyFrame sky) { double rr = angle; FitsImage* ptr = currentContext->cfits; @@ -50,17 +43,8 @@ double Base::mapAngleToRef(double angle, Coord::CoordSystem sys, Coord::SkyFrame case Coord::AMPLIFIER: break; default: - switch (ptr->getWCSOrientation(sys,sky)) { - case Coord::NORMAL: - rr -= ptr->getWCSRotation(sys,sky); - break; - case Coord::XX: - rr = -(angle + ptr->getWCSRotation(sys,sky) + M_PI); - break; - case Coord::YY: - case Coord::XY: - break; - } + rr -= ptr->getWCSRotation(sys,sky); + break; } return zeroTWOPI(rr); -- cgit v0.12