1 2 3 4 5 6 7 8 9 10 11 12
/* ./funcalc -f wcs.in $S foo.fits */ local double dval1, dval2; void *wcs; int init=0; end if( !init++ ) FunInfoGet(fun, FUN_WCS, &wcs, 0); pix2wcs(wcs, (double)cur->x, (double)cur->y, &dval1, &dval2); cur->RA:D = dval1; cur->Dec:D = dval2;