diff options
author | William Joye <wjoye@cfa.harvard.edu> | 2017-11-09 21:15:34 (GMT) |
---|---|---|
committer | William Joye <wjoye@cfa.harvard.edu> | 2017-11-09 21:15:34 (GMT) |
commit | b0f7b82f117089df97bc53d6e4220c2893f09116 (patch) | |
tree | 50a5e652ad8cd45a6febed2f3122e1a0a499c9f3 /tksao/frame | |
parent | 092bfbf6d9154407ad8f36b505b07230f022d2d3 (diff) | |
download | blt-b0f7b82f117089df97bc53d6e4220c2893f09116.zip blt-b0f7b82f117089df97bc53d6e4220c2893f09116.tar.gz blt-b0f7b82f117089df97bc53d6e4220c2893f09116.tar.bz2 |
update AST WCS
Diffstat (limited to 'tksao/frame')
-rw-r--r-- | tksao/frame/fitsimage.C | 12 |
1 files changed, 5 insertions, 7 deletions
diff --git a/tksao/frame/fitsimage.C b/tksao/frame/fitsimage.C index 1cbba07..8e45946 100644 --- a/tksao/frame/fitsimage.C +++ b/tksao/frame/fitsimage.C @@ -4211,8 +4211,6 @@ void FitsImage::astWCSTran(AstFrameSet* ast, int npoint, int forward, double* xout, double* yout) { - // astTran2(ast, npoint, xin, yin, forward, xout, yout); - int naxes = astGetI(ast,"Naxes"); switch (naxes) { case 1: @@ -4228,7 +4226,7 @@ void FitsImage::astWCSTran(AstFrameSet* ast, int npoint, ptr_in[1] = (double*)yin; ptr_in[2] = new double[npoint]; for (int kk=0; kk<npoint; kk++) - ptr_in[2][kk] = 1; + ptr_in[2][kk] = forward ? context_->slice(2) : 0; double* ptr_out[3]; @@ -4252,8 +4250,8 @@ void FitsImage::astWCSTran(AstFrameSet* ast, int npoint, ptr_in[2] = new double[npoint]; ptr_in[3] = new double[npoint]; for (int kk=0; kk<npoint; kk++) { - ptr_in[2][kk] = 1; - ptr_in[3][kk] = 1; + ptr_in[2][kk] = forward ? context_->slice(2) : 0; + ptr_in[3][kk] = forward ? context_->slice(3) : 0; } double* ptr_out[4]; @@ -4266,10 +4264,10 @@ void FitsImage::astWCSTran(AstFrameSet* ast, int npoint, if (ptr_in[2]) delete [] ptr_in[2]; - if (ptr_in[3]) - delete [] ptr_in[3]; if (ptr_out[2]) delete [] ptr_out[2]; + if (ptr_in[3]) + delete [] ptr_in[3]; if (ptr_out[3]) delete [] ptr_out[3]; } |