summaryrefslogtreecommitdiffstats
diff options
context:
space:
mode:
authorWilliam Joye <wjoye@cfa.harvard.edu>2017-11-10 21:52:56 (GMT)
committerWilliam Joye <wjoye@cfa.harvard.edu>2017-11-10 21:52:56 (GMT)
commit15f5271a23904beb595c8269b6d723b2743ef3bb (patch)
tree658860a2697874d834a334d95626e60257e96186
parent9e71ccb6ed107120b1e0a5f324a92c6529143d7e (diff)
downloadblt-15f5271a23904beb595c8269b6d723b2743ef3bb.zip
blt-15f5271a23904beb595c8269b6d723b2743ef3bb.tar.gz
blt-15f5271a23904beb595c8269b6d723b2743ef3bb.tar.bz2
update AST WCS
-rw-r--r--tksao/frame/fitsimage.C46
-rw-r--r--tksao/frame/fitsimage.h2
-rw-r--r--tksao/frame/fitsmap.C4
-rw-r--r--tksao/frame/grid25d.C2
-rw-r--r--tksao/frame/grid2d.C2
-rw-r--r--tksao/frame/grid3d.C2
6 files changed, 29 insertions, 29 deletions
diff --git a/tksao/frame/fitsimage.C b/tksao/frame/fitsimage.C
index 61d6aee..8a13d8f 100644
--- a/tksao/frame/fitsimage.C
+++ b/tksao/frame/fitsimage.C
@@ -1440,13 +1440,13 @@ void FitsImage::match(const char* xxname1, const char* yyname1,
int ss1 = sys1-Coord::WCS;
if (!(ss1>=0 && ast_ && ast_[ss1]))
return;
- if (!astWCSIsASkyFrame(ast_[ss1]))
+ if (!wcsIsASkyFrame(ast_[ss1]))
return;
int ss2 = sys2-Coord::WCS;
if (!(ss2>=0 && ast_ && ast_[ss2]))
return;
- if (!astWCSIsASkyFrame(ast_[ss2]))
+ if (!wcsIsASkyFrame(ast_[ss2]))
return;
// get doubles
@@ -1591,7 +1591,7 @@ void FitsImage::match(const char* xxname1, const char* yyname1,
Tcl_GetDoubleFromObj(interp_, objyy2[ii], yy2+ii);
setAstWCSSystem(newast_, sys1);
- if (!astWCSIsASkyFrame(newast_))
+ if (!wcsIsASkyFrame(newast_))
return;
setAstWCSSkyFrame(newast_, sky1);
for (int ii=0; ii<nxx1; ii++) {
@@ -1600,7 +1600,7 @@ void FitsImage::match(const char* xxname1, const char* yyname1,
}
setAstWCSSystem(newast_, sys2);
- if (!astWCSIsASkyFrame(newast_))
+ if (!wcsIsASkyFrame(newast_))
return;
setAstWCSSkyFrame(newast_, sky2);
for (int ii=0; ii<nxx2; ii++) {
@@ -2930,7 +2930,7 @@ double FitsImage::getWCSPixelSize(Coord::CoordSystem sys)
pt2[1] = wyy[2];
double out = (astDistance(newast_,pt0,pt1)+astDistance(newast_,pt0,pt2))/2.;
- if (astWCSIsASkyFrame(newast_))
+ if (wcsIsASkyFrame(newast_))
return radToDeg(out);
else
return out;
@@ -2967,7 +2967,7 @@ double FitsImage::getWCSPixelArea(Coord::CoordSystem sys)
double ll = astDistance(newast_,pt0,pt1);
double mm = astDistance(newast_,pt0,pt2);
- if (astWCSIsASkyFrame(newast_))
+ if (wcsIsASkyFrame(newast_))
return radToDeg(ll)*radToDeg(mm);
else
return ll*mm;
@@ -3037,7 +3037,7 @@ Coord::Orientation FitsImage::getWCSOrientation(Coord::CoordSystem sys,
Coord::Orientation rr = Coord::NORMAL;
if (!(isnan(ang)||isinf(ang)||(ang == -DBL_MAX)||(ang == DBL_MAX))) {
- if (astWCSIsASkyFrame(newast_))
+ if (wcsIsASkyFrame(newast_))
rr = ang>=0 ? Coord::NORMAL : Coord::XX;
else
rr = ang<=0 ? Coord::NORMAL : Coord::XX;
@@ -3126,7 +3126,7 @@ Vector FitsImage::pix2wcs(Vector in, Coord::CoordSystem sys,
int ss = sys-Coord::WCS;
if (ss>=0 && ast_ && ast_[ss]) {
- if (astWCSIsASkyFrame(ast_[ss])) {
+ if (wcsIsASkyFrame(ast_[ss])) {
setAstWCSSkyFrame(ast_[ss],sky);
Vector out = wcsTran(ast_[ss], in, 1);
if (astOK && checkWCS(out))
@@ -3160,7 +3160,7 @@ Vector FitsImage::pix2wcs(Vector in, Coord::CoordSystem sys,
double yy =0;
wcsTran(newast_, 1, in.v, in.v+1, 1, &xx, &yy);
if (astOK && checkWCS(xx,yy)) {
- if (astWCSIsASkyFrame(newast_))
+ if (wcsIsASkyFrame(newast_))
return Vector(radToDeg(xx),yy*180./M_PI);
else
return Vector(xx,yy);
@@ -3183,7 +3183,7 @@ char* FitsImage::pix2wcs(Vector in, Coord::CoordSystem sys,
int ss = sys-Coord::WCS;
if (ss>=0 && ast_ && ast_[ss]) {
ostringstream str;
- if (astWCSIsASkyFrame(ast_[ss])) {
+ if (wcsIsASkyFrame(ast_[ss])) {
setAstWCSSkyFrame(ast_[ss],sky);
Vector out = wcsTran(ast_[ss], in, 1);
if (!(astOK && checkWCS(out))) {
@@ -3260,7 +3260,7 @@ char* FitsImage::pix2wcs(Vector in, Coord::CoordSystem sys,
wcsTran(newast_, 1, in.v, in.v+1, 1, &xx, &yy);
if (astOK && checkWCS(xx,yy)) {
- if (astWCSIsASkyFrame(newast_)) {
+ if (wcsIsASkyFrame(newast_)) {
switch (format) {
case Coord::DEGREES:
xx =radToDeg(xx); // 0 to 360
@@ -3317,7 +3317,7 @@ Vector FitsImage::wcs2pix(Vector in, Coord::CoordSystem sys,
int ss = sys-Coord::WCS;
if (ss>=0 && ast_ && ast_[ss]) {
- if (astWCSIsASkyFrame(ast_[ss])) {
+ if (wcsIsASkyFrame(ast_[ss])) {
setAstWCSSkyFrame(ast_[ss],sky);
Vector out = wcsTran(ast_[ss], in.degToRad(), 0);
if (astOK && checkWCS(out))
@@ -3350,7 +3350,7 @@ Vector FitsImage::wcs2pix(Vector in, Coord::CoordSystem sys,
double xx =0;
double yy =0;
- if (astWCSIsASkyFrame(newast_))
+ if (wcsIsASkyFrame(newast_))
in *= M_PI/180.;
wcsTran(newast_, 1, in.v, in.v+1, 0, &xx, &yy);
@@ -3372,7 +3372,7 @@ double FitsImage::getWCSDist(Vector a, Vector b, Coord::CoordSystem sys)
astClearStatus; // just to make sure
double rr=0;
- if (astWCSIsASkyFrame(ast_[ss])) {
+ if (wcsIsASkyFrame(ast_[ss])) {
Vector aa = a*M_PI/180.;
Vector bb = b*M_PI/180.;
rr = astDistance(ast_[ss], aa.v, bb.v) *180./M_PI;
@@ -3392,7 +3392,7 @@ double FitsImage::getWCSDist(Vector aa, Vector bb, Coord::CoordSystem sys)
setAstWCSSystem(newast_,sys);
double rr=0;
- if (astWCSIsASkyFrame(newast_)) {
+ if (wcsIsASkyFrame(newast_)) {
aa *= M_PI/180.;
bb *= M_PI/180.;
rr = astDistance(newast_, aa.v, bb.v) *180./M_PI;
@@ -3417,7 +3417,7 @@ int FitsImage::hasWCSEqu(Coord::CoordSystem sys)
int ss = sys-Coord::WCS;
if (ss>=0 && ast_ && ast_[ss])
- if (astWCSIsASkyFrame(ast_[ss])) {
+ if (wcsIsASkyFrame(ast_[ss])) {
// special case of xLON/xLAT
char* bb = &(wcs_[ss]->c1type[1]);
if (!strncmp(bb,"LON",3) || !strncmp(bb,"LAT",3)) {
@@ -3449,7 +3449,7 @@ int FitsImage::hasWCSCel(Coord::CoordSystem sys)
int ss = sys-Coord::WCS;
if (ss>=0 && ast_ && ast_[ss])
- if (astWCSIsASkyFrame(ast_[ss]))
+ if (wcsIsASkyFrame(ast_[ss]))
return 1;
return 0;
@@ -3766,7 +3766,7 @@ void FitsImage::astinit(int ss, FitsHead* hd, FitsHead* prim)
#endif
// set default skyframe
- if (astWCSIsASkyFrame(ast_[ss]))
+ if (wcsIsASkyFrame(ast_[ss]))
setAstWCSSkyFrame(ast_[ss],Coord::FK5);
}
@@ -3823,7 +3823,7 @@ void FitsImage::astinit(FitsHead* hd, FitsHead* prim)
break;
}
- if (astWCSIsASkyFrame(newast_))
+ if (wcsIsASkyFrame(newast_))
setAstWCSSkyFrame(newast_,Coord::FK5);
}
#endif
@@ -3841,7 +3841,7 @@ void FitsImage::astinit0(int ss, FitsHead* hd, FitsHead* prim)
return;
// set default skyframe
- if (astWCSIsASkyFrame(ast_[ss]))
+ if (wcsIsASkyFrame(ast_[ss]))
setAstWCSSkyFrame(ast_[ss],Coord::FK5);
}
@@ -3877,7 +3877,7 @@ void FitsImage::setAstWCSFormat(AstFrameSet* aa, int id, const char* format)
void FitsImage::setAstWCSSkyFrame(AstFrameSet* ast, Coord::SkyFrame sky)
{
// is sky frame
- if (!astWCSIsASkyFrame(ast))
+ if (!wcsIsASkyFrame(ast))
return;
// is it already set?
@@ -3965,12 +3965,12 @@ void FitsImage::setAstWCSSystem(AstFrameSet* ast, Coord::CoordSystem sys)
#endif
#ifndef NEWWCS
-int FitsImage::astWCSIsASkyFrame(void* ast)
+int FitsImage::wcsIsASkyFrame(void* ast)
{
return astIsASkyFrame(astGetFrame(ast,AST__CURRENT));
}
#else
-int FitsImage::astWCSIsASkyFrame(void* ast)
+int FitsImage::wcsIsASkyFrame(void* ast)
{
int naxes = astGetI(ast,"Naxes");
switch (naxes) {
diff --git a/tksao/frame/fitsimage.h b/tksao/frame/fitsimage.h
index 72d25ee..f2d5222 100644
--- a/tksao/frame/fitsimage.h
+++ b/tksao/frame/fitsimage.h
@@ -398,7 +398,7 @@ class FitsImage {
double getWCSPixelArea(Coord::CoordSystem);
#endif
- int astWCSIsASkyFrame(void*);
+ int wcsIsASkyFrame(void*);
#ifndef NEWWCS
Vector wcsTran(AstFrameSet*, Vector&, int);
void wcsTran(AstFrameSet*, int, Vector*, int, Vector*);
diff --git a/tksao/frame/fitsmap.C b/tksao/frame/fitsmap.C
index 6f0e975..f9b220b 100644
--- a/tksao/frame/fitsmap.C
+++ b/tksao/frame/fitsmap.C
@@ -195,7 +195,7 @@ double FitsImage::mapLenFromRef(double dd, Coord::CoordSystem sys,
pt1[1] = wyy[1];
double out = astDistance(newast_,pt0,pt1);
- if (astWCSIsASkyFrame(newast_)) {
+ if (wcsIsASkyFrame(newast_)) {
out = radToDeg(out);
switch (dist) {
case Coord::DEGREE:
@@ -295,7 +295,7 @@ double FitsImage::mapLenToRef(double dd, Coord::CoordSystem sys,
maperr =0;
double rdd = dd;
- if (astWCSIsASkyFrame(newast_)) {
+ if (wcsIsASkyFrame(newast_)) {
rdd = degToRad(dd);
switch (dist) {
case Coord::DEGREE:
diff --git a/tksao/frame/grid25d.C b/tksao/frame/grid25d.C
index ec34d32..24fc7a6 100644
--- a/tksao/frame/grid25d.C
+++ b/tksao/frame/grid25d.C
@@ -67,7 +67,7 @@ int Grid25d::doit(RenderMode rm)
// Get 2D SkyFrame
AstFrameSet* ast = (AstFrameSet*)astCopy(fits->getAST(system_));
- if (fits->astWCSIsASkyFrame(ast))
+ if (fits->wcsIsASkyFrame(ast))
fits->setAstWCSSkyFrame(ast, sky_);
// Record the index of the current Frame
diff --git a/tksao/frame/grid2d.C b/tksao/frame/grid2d.C
index 7494868..336cd20 100644
--- a/tksao/frame/grid2d.C
+++ b/tksao/frame/grid2d.C
@@ -72,7 +72,7 @@ int Grid2d::doit(RenderMode rm)
AstFrameSet* ast = (AstFrameSet*)astCopy(fits->getAST(system_));
// set desired skyformat
- if (fits->astWCSIsASkyFrame(ast))
+ if (fits->wcsIsASkyFrame(ast))
fits->setAstWCSSkyFrame(ast, sky_);
// add wcs to frameset
diff --git a/tksao/frame/grid3d.C b/tksao/frame/grid3d.C
index 13f1991..774a393 100644
--- a/tksao/frame/grid3d.C
+++ b/tksao/frame/grid3d.C
@@ -100,7 +100,7 @@ int Grid3d::doit(RenderMode rm)
// Get 2D SkyFrame
AstFrameSet* ast = (AstFrameSet*)astCopy(fits->getAST(system_));
- if (fits->astWCSIsASkyFrame(ast))
+ if (fits->wcsIsASkyFrame(ast))
fits->setAstWCSSkyFrame(ast, sky_);
// Record the index of the current Frame