diff options
author | William Joye <wjoye@cfa.harvard.edu> | 2019-10-11 16:44:40 (GMT) |
---|---|---|
committer | William Joye <wjoye@cfa.harvard.edu> | 2019-10-11 16:44:40 (GMT) |
commit | 54bcabe6518c97ba76f40c18c04229a96b2b6216 (patch) | |
tree | fcc97027842f60fdd8909e7fdf81689e2ca161a3 /tksao/frame/basecommand.C | |
parent | ccf5c66b117b2b3eb6d7356d3558af4f52002a04 (diff) | |
download | blt-54bcabe6518c97ba76f40c18c04229a96b2b6216.zip blt-54bcabe6518c97ba76f40c18c04229a96b2b6216.tar.gz blt-54bcabe6518c97ba76f40c18c04229a96b2b6216.tar.bz2 |
support higher fits dims
Diffstat (limited to 'tksao/frame/basecommand.C')
-rw-r--r-- | tksao/frame/basecommand.C | 16 |
1 files changed, 7 insertions, 9 deletions
diff --git a/tksao/frame/basecommand.C b/tksao/frame/basecommand.C index 071f5d9..a5905b9 100644 --- a/tksao/frame/basecommand.C +++ b/tksao/frame/basecommand.C @@ -1748,18 +1748,17 @@ void Base::getFitsSliceCmd(int ii) printInteger(currentContext->slice(ii)); } -void Base::getFitsSliceFromImageCmd(Coord::CoordSystem sys, Coord::SkyFrame sky) +void Base::getFitsSliceFromImageCmd(Coord::CoordSystem sys) { - getFitsSliceFromImageCmd(currentContext->slice(2), sys, sky); + getFitsSliceFromImageCmd(currentContext->slice(2), sys); } -void Base::getFitsSliceFromImageCmd(int ss, Coord::CoordSystem sys, - Coord::SkyFrame sky) +void Base::getFitsSliceFromImageCmd(int ss, Coord::CoordSystem sys) { if (currentContext->cfits) { FitsImage* ptr = currentContext->fits; Vector3d dd = Vector3d(ptr->center(),ss) * Translate3d(-.5,-.5,-.5); - Vector3d out = ptr->mapFromRef(dd,sys,sky); + Vector3d out = ptr->mapFromRef(dd,sys,Coord::FK5); // set precision high enough for plot in GHz freq ostringstream str; @@ -1770,14 +1769,13 @@ void Base::getFitsSliceFromImageCmd(int ss, Coord::CoordSystem sys, Tcl_AppendResult(interp, "1", NULL); } -void Base::getFitsSliceToImageCmd(double dd, Coord::CoordSystem sys, - Coord::SkyFrame sky) +void Base::getFitsSliceToImageCmd(double dd, Coord::CoordSystem sys) { if (currentContext->cfits) { FitsImage* ptr = currentContext->fits; Vector3d cc = Vector3d(ptr->center(),1) * Translate3d(-.5, -.5, -.5); - Vector3d wcc = ptr->mapFromRef(cc,sys,sky); - Vector3d oo = ptr->mapToRef(Vector3d(wcc[0],wcc[1],dd),sys,sky); + Vector3d wcc = ptr->mapFromRef(cc,sys,Coord::FK5); + Vector3d oo = ptr->mapToRef(Vector3d(wcc[0],wcc[1],dd),sys,Coord::FK5); Vector3d out = oo * Translate3d(.5, .5, .5); printInteger(out[2]); } |