diff options
author | William Joye <wjoye@cfa.harvard.edu> | 2017-12-07 22:08:42 (GMT) |
---|---|---|
committer | William Joye <wjoye@cfa.harvard.edu> | 2017-12-07 22:08:42 (GMT) |
commit | 5606c9ac2deadb681b5aaf91f9045b646d9f2729 (patch) | |
tree | 99bebd97e86803ef020735ae1af3d0da575dd037 /tksao/frame/fitsmap.C | |
parent | 3fbb8a09c8befe6cc8678961155dfcd6447907b5 (diff) | |
download | blt-5606c9ac2deadb681b5aaf91f9045b646d9f2729.zip blt-5606c9ac2deadb681b5aaf91f9045b646d9f2729.tar.gz blt-5606c9ac2deadb681b5aaf91f9045b646d9f2729.tar.bz2 |
update AST WCS
Diffstat (limited to 'tksao/frame/fitsmap.C')
-rw-r--r-- | tksao/frame/fitsmap.C | 24 |
1 files changed, 12 insertions, 12 deletions
diff --git a/tksao/frame/fitsmap.C b/tksao/frame/fitsmap.C index 59924f4..9e93eeb 100644 --- a/tksao/frame/fitsmap.C +++ b/tksao/frame/fitsmap.C @@ -208,17 +208,17 @@ double FitsImage::mapLenFromRef(double dd, Coord::CoordSystem sys, if (!hasWCS(sys)) return 0; - setWCSSystem(newast_, sys); - setWCSSkyFrame(newast_, Coord::FK5); + setWCSSystem(ast_, sys); + setWCSSkyFrame(ast_, Coord::FK5); Vector in[2]; Vector out[2]; in[0] = center(); in[1] = center()+Vector(0,dd); - wcsTran(newast_, 2, in, 1, out); - double rr = wcsDistance(newast_,out[0],out[1]); + wcsTran(ast_, 2, in, 1, out); + double rr = wcsDistance(ast_,out[0],out[1]); - if (wcsIsASkyFrame(newast_)) { + if (wcsIsASkyFrame(ast_)) { rr = radToDeg(rr); switch (dist) { case Coord::DEGREE: @@ -310,10 +310,10 @@ double FitsImage::mapLenToRef(double dd, Coord::CoordSystem sys, return 0; astClearStatus; // just to make sure - setWCSSystem(newast_, sys); + setWCSSystem(ast_, sys); double rdd = dd; - if (wcsIsASkyFrame(newast_)) { + if (wcsIsASkyFrame(ast_)) { rdd = degToRad(dd); switch (dist) { case Coord::DEGREE: @@ -328,12 +328,12 @@ double FitsImage::mapLenToRef(double dd, Coord::CoordSystem sys, } Vector cc = center(); - Vector wcc = wcsTran(newast_,cc,1); + Vector wcc = wcsTran(ast_,cc,1); Vector wpp = wcc+Vector(0,rdd); - Vector pp = wcsTran(newast_,wpp,0); - astInvert(newast_); - double rr = wcsDistance(newast_,cc,pp); - astInvert(newast_); + Vector pp = wcsTran(ast_,wpp,0); + astInvert(ast_); + double rr = wcsDistance(ast_,cc,pp); + astInvert(ast_); return rr; } |