diff options
Diffstat (limited to 'tksao/frame/box.C')
| -rw-r--r-- | tksao/frame/box.C | 424 |
1 files changed, 424 insertions, 0 deletions
diff --git a/tksao/frame/box.C b/tksao/frame/box.C new file mode 100644 index 0000000..4a97e48 --- /dev/null +++ b/tksao/frame/box.C @@ -0,0 +1,424 @@ +// Copyright (C) 1999-2016 +// Smithsonian Astrophysical Observatory, Cambridge, MA, USA +// For conditions of distribution and use, see copyright notice in "copyright" + +#include <tk.h> + +#include "box.h" +#include "fitsimage.h" + +Box::Box(Base* p, const Vector& ctr, const Vector& seg, double ang) + : BaseBox(p, ctr, ang) +{ + numAnnuli_ = 1; + annuli_ = new Vector[1]; + annuli_[0] = seg; + + strcpy(type_,"box"); + numHandle = 4; + + updateBBox(); +} + +Box::Box(Base* p, const Vector& ctr, + const Vector& seg, + double ang, + const char* clr, int* dsh, + int wth, const char* fnt, const char* txt, + unsigned short prop, const char* cmt, + const List<Tag>& tg, const List<CallBack>& cb) + : BaseBox(p, ctr, ang, clr, dsh, wth, fnt, txt, prop, cmt, tg, cb) +{ + numAnnuli_ = 1; + annuli_ = new Vector[1]; + annuli_[0] = seg; + + strcpy(type_,"box"); + numHandle = 4; + + updateBBox(); +} + +Box::Box(const Box& a) : BaseBox(a) {} + +void Box::editBegin(int h) +{ + switch (h) { + case 1: + return; + case 2: + annuli_[0] = Vector(-annuli_[0][0],annuli_[0][1]); + return; + case 3: + annuli_[0] = -annuli_[0]; + return; + case 4: + annuli_[0] = Vector(annuli_[0][0],-annuli_[0][1]); + return; + } + + doCallBack(CallBack::EDITBEGINCB); +} + +void Box::edit(const Vector& v, int h) +{ + Matrix mm = bckMatrix(); + Matrix nn = mm.invert(); + + // This annuli_s about the opposite node + Vector ov = annuli_[0]/2 * nn; + annuli_[0] = (annuli_[0]/2) - (v*mm); + Vector nv = annuli_[0]/2 * nn; + center -= nv-ov; + + updateBBox(); + doCallBack(CallBack::EDITCB); + doCallBack(CallBack::MOVECB); +} + +void Box::editEnd() +{ + annuli_[0] = annuli_[0].abs(); + + updateBBox(); + doCallBack(CallBack::EDITENDCB); +} + +void Box::analysis(AnalysisTask mm, int which) +{ + switch (mm) { + case HISTOGRAM: + if (!analysisHistogram_ && which) { + addCallBack(CallBack::MOVECB, analysisHistogramCB_[0], + parent->options->cmdName); + addCallBack(CallBack::EDITCB, analysisHistogramCB_[0], + parent->options->cmdName); + addCallBack(CallBack::ROTATECB, analysisHistogramCB_[0], + parent->options->cmdName); + addCallBack(CallBack::DELETECB, analysisHistogramCB_[1], + parent->options->cmdName); + } + if (analysisHistogram_ && !which) { + deleteCallBack(CallBack::MOVECB, analysisHistogramCB_[0]); + deleteCallBack(CallBack::EDITCB, analysisHistogramCB_[0]); + deleteCallBack(CallBack::ROTATECB, analysisHistogramCB_[0]); + deleteCallBack(CallBack::DELETECB, analysisHistogramCB_[1]); + } + + analysisHistogram_ = which; + break; + case PLOT3D: + if (!analysisPlot3d_ && which) { + addCallBack(CallBack::MOVECB, analysisPlot3dCB_[0], + parent->options->cmdName); + addCallBack(CallBack::EDITCB, analysisPlot3dCB_[0], + parent->options->cmdName); + addCallBack(CallBack::ROTATECB, analysisPlot3dCB_[0], + parent->options->cmdName); + addCallBack(CallBack::DELETECB, analysisPlot3dCB_[1], + parent->options->cmdName); + } + if (analysisPlot3d_ && !which) { + deleteCallBack(CallBack::MOVECB, analysisPlot3dCB_[0]); + deleteCallBack(CallBack::EDITCB, analysisPlot3dCB_[0]); + deleteCallBack(CallBack::ROTATECB, analysisPlot3dCB_[0]); + deleteCallBack(CallBack::DELETECB, analysisPlot3dCB_[1]); + } + + analysisPlot3d_ = which; + break; + case STATS: + if (!analysisStats_ && which) { + addCallBack(CallBack::MOVECB, analysisStatsCB_[0], + parent->options->cmdName); + addCallBack(CallBack::EDITCB, analysisStatsCB_[0], + parent->options->cmdName); + addCallBack(CallBack::ROTATECB, analysisStatsCB_[0], + parent->options->cmdName); + addCallBack(CallBack::UPDATECB, analysisStatsCB_[0], + parent->options->cmdName); + addCallBack(CallBack::DELETECB, analysisStatsCB_[1], + parent->options->cmdName); + } + if (analysisStats_ && !which) { + deleteCallBack(CallBack::MOVECB, analysisStatsCB_[0]); + deleteCallBack(CallBack::EDITCB, analysisStatsCB_[0]); + deleteCallBack(CallBack::UPDATECB, analysisStatsCB_[0]); + deleteCallBack(CallBack::ROTATECB, analysisStatsCB_[0]); + deleteCallBack(CallBack::DELETECB, analysisStatsCB_[1]); + } + + analysisStats_ = which; + break; + default: + // na + break; + } +} + +void Box::analysisHistogram(char* xname, char* yname, int num) +{ + double* x; + double* y; + Matrix mm = Rotate(angle) * Translate(center); + + // during resize, annuli_ can be negative + Vector vv = annuli_[0].abs(); + BBox bb(-vv * mm); + bb.bound( vv * mm); + bb.bound(Vector( vv[0],-vv[1]) * mm); + bb.bound(Vector(-vv[0], vv[1]) * mm); + + parent->markerAnalysisHistogram(this, &x, &y, bb, num); + analysisXYResult(xname, yname, x, y, num+1); +} + +void Box::analysisPlot3d(char* xname, char* yname, + Coord::CoordSystem sys, + Marker::AnalysisMethod method) +{ + double* x; + double* y; + Matrix mm = Rotate(angle) * Translate(center); + + // during resize, annuli_ can be negative + Vector vv = annuli_[0].abs(); + BBox bb(-vv * mm); + bb.bound( vv * mm); + bb.bound(Vector( vv[0],-vv[1]) * mm); + bb.bound(Vector(-vv[0], vv[1]) * mm); + + int num = parent->markerAnalysisPlot3d(this, &x, &y, bb, sys, method); + analysisXYResult(xname, yname, x, y, num); +} + +void Box::analysisStats(Coord::CoordSystem sys, Coord::SkyFrame sky) +{ + ostringstream str; + Matrix mm = Rotate(angle) * Translate(center); + + // during resize, annuli_ can be negative + Vector vv = annuli_[0].abs(); + BBox bb(-vv * mm); + bb.bound( vv * mm); + bb.bound(Vector( vv[0],-vv[1]) * mm); + bb.bound(Vector(-vv[0], vv[1]) * mm); + + parent->markerAnalysisStats(this, str, bb, sys, sky); + str << ends; + Tcl_AppendResult(parent->interp, str.str().c_str(), NULL); +} + +// list + +void Box::list(ostream& str, Coord::CoordSystem sys, Coord::SkyFrame sky, + Coord::SkyFormat format, int conj, int strip) +{ + FitsImage* ptr = parent->findFits(sys,center); + listPre(str, sys, sky, ptr, strip, 0); + + switch (sys) { + case Coord::IMAGE: + case Coord::PHYSICAL: + case Coord::DETECTOR: + case Coord::AMPLIFIER: + listNonCel(ptr, str, sys); + break; + default: + if (ptr->hasWCSCel(sys)) { + Vector rr = ptr->mapLenFromRef(annuli_[0],sys,Coord::ARCSEC); + double aa = parent->mapAngleFromRef(angle,sys,sky); + switch (format) { + case Coord::DEGREES: + { + Vector vv = ptr->mapFromRef(center,sys,sky); + str << type_ << '(' + << setprecision(10) << vv << ',' + << setprecision(3) << fixed << setunit('"') << rr << ','; + str.unsetf(ios_base::floatfield); + str << setprecision(8) << radToDeg(aa) << ')'; + } + break; + case Coord::SEXAGESIMAL: + listRADEC(ptr,center,sys,sky,format); + str << type_ << '(' << setprecision(8) + << ra << ',' << dec << ',' + << setprecision(3) << fixed << setunit('"') << rr << ','; + str.unsetf(ios_base::floatfield); + str << setprecision(8) << radToDeg(aa) << ')'; + break; + } + } + else + listNonCel(ptr, str, sys); + } + + listPost(str, conj, strip); +} + +void Box::listNonCel(FitsImage* ptr, ostream& str, Coord::CoordSystem sys) +{ + Vector vv = ptr->mapFromRef(center,sys); + Vector rr = ptr->mapLenFromRef(annuli_[0],sys); + double aa = parent->mapAngleFromRef(angle,sys); + str << type_ << '(' << setprecision(8) << vv << ',' << rr << ',' + << radToDeg(aa) << ')'; +} + +void Box::listXML(ostream& str, Coord::CoordSystem sys, Coord::SkyFrame sky, + Coord::SkyFormat format) +{ + FitsImage* ptr = parent->findFits(sys,center); + + XMLRowInit(); + XMLRow(XMLSHAPE,type_); + + XMLRowCenter(ptr,sys,sky,format); + XMLRowRadius(ptr,sys,annuli_[0]); + XMLRowAng(sys,sky); + + XMLRowProps(ptr,sys); + XMLRowEnd(str); +} + +void Box::listCiao(ostream& str, Coord::CoordSystem sys, int strip) +{ + FitsImage* ptr = parent->findFits(); + listCiaoPre(str); + + // radius is always in image coords + switch (sys) { + case Coord::IMAGE: + case Coord::PHYSICAL: + case Coord::DETECTOR: + case Coord::AMPLIFIER: + { + Vector vv = ptr->mapFromRef(center,Coord::PHYSICAL); + Vector rr = ptr->mapLenFromRef(annuli_[0],Coord::PHYSICAL); + str << type_ << '(' << setprecision(8) << vv << ',' << rr << ',' + << radToDeg(angle) << ')'; + } + break; + default: + if (ptr->hasWCSCel(sys)) { + listRADEC(ptr,center,sys,Coord::FK5,Coord::SEXAGESIMAL); + Vector rr = ptr->mapLenFromRef(annuli_[0],sys,Coord::ARCMIN); + str << type_ << '(' + << ra << ',' << dec << ',' + << setprecision(5) << fixed << setunit('\'') << rr << ','; + str.unsetf(ios_base::floatfield); + str << setprecision(8) << radToDeg(angle) << ')'; + } + } + + listCiaoPost(str, strip); +} + +void Box::listSAOtng(ostream& str, Coord::CoordSystem sys, Coord::SkyFrame sky, + Coord::SkyFormat format, int strip) +{ + FitsImage* ptr = parent->findFits(); + listSAOtngPre(str, strip); + + // radius is always in image coords + switch (sys) { + case Coord::IMAGE: + case Coord::PHYSICAL: + case Coord::DETECTOR: + case Coord::AMPLIFIER: + { + Vector vv = ptr->mapFromRef(center,Coord::IMAGE); + Vector rr = ptr->mapLenFromRef(annuli_[0],Coord::IMAGE); + str << type_ << '(' << setprecision(8) << vv << ',' << rr << ',' + << radToDeg(angle) << ')'; + } + break; + default: + if (ptr->hasWCSCel(sys)) { + Vector rr = ptr->mapLenFromRef(annuli_[0],Coord::IMAGE); + switch (format) { + case Coord::DEGREES: + { + Vector vv = ptr->mapFromRef(center,sys,sky); + str << type_ << '(' + << setprecision(10) << vv << ',' + << setprecision(8) << rr << ',' + << setprecision(8) << radToDeg(angle) << ')'; + } + break; + case Coord::SEXAGESIMAL: + listRADEC(ptr,center,sys,sky,format); + str << type_ << '(' + << ra << ',' << dec << ',' + << setprecision(8) << rr << ',' + << setprecision(8) << radToDeg(angle) << ')'; + break; + } + } + } + + listSAOtngPost(str,strip); +} + +void Box::listPros(ostream& str, Coord::CoordSystem sys, Coord::SkyFrame sky, + Coord::SkyFormat format, int strip) +{ + FitsImage* ptr = parent->findFits(); + + switch (sys) { + case Coord::IMAGE: + case Coord::DETECTOR: + case Coord::AMPLIFIER: + sys = Coord::IMAGE; + case Coord::PHYSICAL: + { + Vector vv = ptr->mapFromRef(center,sys); + Vector rr = ptr->mapLenFromRef(annuli_[0],Coord::IMAGE); + coord.listProsCoordSystem(str,sys,sky); + str << "; "<< type_ << ' ' << setprecision(8) << vv << ' ' << rr << ' ' + << radToDeg(angle); + } + break; + default: + if (ptr->hasWCSCel(sys)) { + Vector rr = ptr->mapLenFromRef(annuli_[0],sys,Coord::ARCSEC); + switch (format) { + case Coord::DEGREES: + { + Vector vv = ptr->mapFromRef(center,sys,sky); + coord.listProsCoordSystem(str,sys,sky); + str << "; " << type_ << ' ' + << setprecision(10) << setunit('d') << vv << ' ' + << setprecision(3) << fixed << setunit('"') << rr << ' '; + str.unsetf(ios_base::floatfield); + str << setprecision(8) << radToDeg(angle); + } + break; + case Coord::SEXAGESIMAL: + listRADECPros(ptr,center,sys,sky,format); + coord.listProsCoordSystem(str,sys,sky); + str << "; " << type_ << ' ' + << ra << ' ' << dec << ' ' + << setprecision(3) << fixed << setunit('"') << rr << ' '; + str.unsetf(ios_base::floatfield); + str << setprecision(8) << radToDeg(angle); + break; + } + } + } + + listProsPost(str, strip); +} + +void Box::listSAOimage(ostream& str, int strip) +{ + FitsImage* ptr = parent->findFits(); + listSAOimagePre(str); + + Vector vv = ptr->mapFromRef(center,Coord::IMAGE); + str << type_ << '(' << setprecision(8) << vv << ',' << annuli_[0] << ',' + << radToDeg(angle) << ')'; + + listSAOimagePost(str, strip); +} + |