1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209
210
211
212
213
214
215
216
217
218
219
220
221
222
223
224
225
226
227
228
229
230
231
232
233
234
235
236
237
238
239
240
241
242
243
244
245
246
247
248
249
250
251
252
253
254
255
256
257
258
259
260
261
262
263
264
265
266
267
268
269
270
271
272
273
274
275
276
277
278
279
280
281
282
283
284
285
286
287
288
289
290
291
292
293
294
295
296
297
298
299
300
301
302
303
304
305
306
307
308
309
310
311
312
313
314
315
316
317
318
319
320
321
322
323
324
325
326
327
328
329
330
331
332
333
334
335
336
337
338
339
340
341
342
343
344
345
346
347
348
349
350
351
352
353
354
355
356
357
358
359
360
361
362
363
364
365
366
367
368
369
370
371
372
373
374
375
376
377
378
379
380
381
382
383
384
385
386
387
388
389
390
391
392
393
394
395
396
397
398
399
400
401
402
403
404
405
406
407
408
409
410
411
412
413
414
415
416
417
418
419
420
421
422
423
424
425
426
427
428
429
430
431
432
433
434
435
436
437
438
439
440
441
442
443
444
445
446
447
448
449
450
451
452
453
454
455
456
457
458
459
460
461
462
463
464
465
466
467
468
469
470
471
472
473
474
475
476
|
#ifndef wcslib_h_
#define wcslib_h_
/*=============================================================================
*
* WCSLIB - an implementation of the FITS WCS proposal.
* Copyright (C) 1995-2002, Mark Calabretta
*
* This library is free software; you can redistribute it and/or
* modify it under the terms of the GNU Lesser General Public
* License as published by the Free Software Foundation; either
* version 2 of the License, or (at your option) any later version.
*
* This library is distributed in the hope that it will be useful,
* but WITHOUT ANY WARRANTY; without even the implied warranty of
* MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
* Lesser General Public License for more details.
*
* You should have received a copy of the GNU Lesser General Public
* License along with this library; if not, write to the Free Software
* Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA 02111-1307 USA
*
* Correspondence concerning WCSLIB may be directed to:
* Internet email: mcalabre@atnf.csiro.au
* Postal address: Dr. Mark Calabretta,
* Australia Telescope National Facility,
* P.O. Box 76,
* Epping, NSW, 2121,
* AUSTRALIA
*
* Author: Mark Calabretta, Australia Telescope National Facility
* $Id: wcslib.h,v 1.2 2016/03/30 20:09:45 joye Exp $
*===========================================================================*/
#ifdef __cplusplus
extern "C" {
#endif
#if !defined(__STDC__) && !defined(__cplusplus)
#ifndef const
#define const
#endif
#endif
#define MAXPV 100
#define WCS_NGRIDPOINTS 12 /* Number of WCS grid points / axis */
#define WCS_NGRIDPOINTS2 (WCS_NGRIDPOINTS*WCS_NGRIDPOINTS)
#define WCS_INVMAXDEG 9 /* Maximum inversion polynom degree */
#define WCS_INVACCURACY 0.04 /* Maximum inversion error (pixels) */
#define WCS_NRANGEPOINTS 32 /* Number of WCS range points / axis */
#ifndef PI
#define PI 3.1415926535898 /* never met before? */
#endif
/* DEG/ARCSEC is now D2S and ARCSEC/DEG is S2D */
/* #define DEG (PI/180.0) 1 deg in radians */
/* #define ARCSEC (DEG/3600.0) 1 arcsec in radians */
#define NAXISPV 2
/* poly.h
*%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
* Part of: A program using polynomial fits
* Author: E.BERTIN (IAP)
* Contents: Include for poly.c
* Last modified: 03/03/2004
*%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
*/
#ifndef _POLY_H_
#define _POLY_H_
/*--------------------------------- constants -------------------------------*/
#define POLY_MAXDIM 4 /* Max dimensionality of polynom */
#define POLY_MAXDEGREE 10 /* Max degree of the polynom */
/*---------------------------------- macros ---------------------------------*/
/*--------------------------- structure definitions -------------------------*/
typedef struct poly
{
double *basis; /* Current values of the basis functions */
double *coeff; /* Polynom coefficients */
int ncoeff; /* Number of coefficients */
int *group; /* Groups */
int ndim; /* dimensionality of the polynom */
int *degree; /* Degree in each group */
int ngroup; /* Number of different groups */
} polystruct;
/*---------------------------------- protos --------------------------------*/
extern polystruct *poly_init(int *group,int ndim,int *degree,int ngroup);
extern double poly_func(polystruct *poly, double *pos);
extern int cholsolve(double *a, double *b, int n),
*poly_powers(polystruct *poly);
extern void poly_addcste(polystruct *poly, double *cste),
poly_end(polystruct *poly),
poly_fit(polystruct *poly, double *x, double *y,
double *w, int ndata, double *extbasis),
poly_solve(double *a, double *b, int n),
svdsolve(double *a, double *b, int m, int n,
double *vmat, double *wmat);
#endif
extern int npcode;
extern char pcodes[26][4];
struct prjprm {
char code[4];
int flag;
double phi0, theta0;
double r0;
double p[10];
double w[20];
int n;
int npv;
double ppv[2*MAXPV];
struct poly *inv_x;
struct poly *inv_y;
#if __STDC__ || defined(__cplusplus)
int (*prjfwd)(const double, const double,
struct prjprm *,
double *, double *);
int (*prjrev)(const double, const double,
struct prjprm *,
double *, double *);
#else
int (*prjfwd)();
int (*prjrev)();
#endif
};
#if __STDC__ || defined(__cplusplus)
int prjset(const char [], struct prjprm *);
int prjfwd(const double, const double, struct prjprm *, double *, double *);
int prjrev(const double, const double, struct prjprm *, double *, double *);
int azpset(struct prjprm *);
int azpfwd(const double, const double, struct prjprm *, double *, double *);
int azprev(const double, const double, struct prjprm *, double *, double *);
int szpset(struct prjprm *);
int szpfwd(const double, const double, struct prjprm *, double *, double *);
int szprev(const double, const double, struct prjprm *, double *, double *);
int tanset(struct prjprm *);
int tanfwd(const double, const double, struct prjprm *, double *, double *);
int tanrev(const double, const double, struct prjprm *, double *, double *);
int stgset(struct prjprm *);
int stgfwd(const double, const double, struct prjprm *, double *, double *);
int stgrev(const double, const double, struct prjprm *, double *, double *);
int sinset(struct prjprm *);
int sinfwd(const double, const double, struct prjprm *, double *, double *);
int sinrev(const double, const double, struct prjprm *, double *, double *);
int arcset(struct prjprm *);
int arcfwd(const double, const double, struct prjprm *, double *, double *);
int arcrev(const double, const double, struct prjprm *, double *, double *);
int zpnset(struct prjprm *);
int zpnfwd(const double, const double, struct prjprm *, double *, double *);
int zpnrev(const double, const double, struct prjprm *, double *, double *);
int zeaset(struct prjprm *);
int zeafwd(const double, const double, struct prjprm *, double *, double *);
int zearev(const double, const double, struct prjprm *, double *, double *);
int airset(struct prjprm *);
int airfwd(const double, const double, struct prjprm *, double *, double *);
int airrev(const double, const double, struct prjprm *, double *, double *);
int cypset(struct prjprm *);
int cypfwd(const double, const double, struct prjprm *, double *, double *);
int cyprev(const double, const double, struct prjprm *, double *, double *);
int ceaset(struct prjprm *);
int ceafwd(const double, const double, struct prjprm *, double *, double *);
int cearev(const double, const double, struct prjprm *, double *, double *);
int carset(struct prjprm *);
int carfwd(const double, const double, struct prjprm *, double *, double *);
int carrev(const double, const double, struct prjprm *, double *, double *);
int merset(struct prjprm *);
int merfwd(const double, const double, struct prjprm *, double *, double *);
int merrev(const double, const double, struct prjprm *, double *, double *);
int sflset(struct prjprm *);
int sflfwd(const double, const double, struct prjprm *, double *, double *);
int sflrev(const double, const double, struct prjprm *, double *, double *);
int parset(struct prjprm *);
int parfwd(const double, const double, struct prjprm *, double *, double *);
int parrev(const double, const double, struct prjprm *, double *, double *);
int molset(struct prjprm *);
int molfwd(const double, const double, struct prjprm *, double *, double *);
int molrev(const double, const double, struct prjprm *, double *, double *);
int aitset(struct prjprm *);
int aitfwd(const double, const double, struct prjprm *, double *, double *);
int aitrev(const double, const double, struct prjprm *, double *, double *);
int copset(struct prjprm *);
int copfwd(const double, const double, struct prjprm *, double *, double *);
int coprev(const double, const double, struct prjprm *, double *, double *);
int coeset(struct prjprm *);
int coefwd(const double, const double, struct prjprm *, double *, double *);
int coerev(const double, const double, struct prjprm *, double *, double *);
int codset(struct prjprm *);
int codfwd(const double, const double, struct prjprm *, double *, double *);
int codrev(const double, const double, struct prjprm *, double *, double *);
int cooset(struct prjprm *);
int coofwd(const double, const double, struct prjprm *, double *, double *);
int coorev(const double, const double, struct prjprm *, double *, double *);
int bonset(struct prjprm *);
int bonfwd(const double, const double, struct prjprm *, double *, double *);
int bonrev(const double, const double, struct prjprm *, double *, double *);
int pcoset(struct prjprm *);
int pcofwd(const double, const double, struct prjprm *, double *, double *);
int pcorev(const double, const double, struct prjprm *, double *, double *);
int tscset(struct prjprm *);
int tscfwd(const double, const double, struct prjprm *, double *, double *);
int tscrev(const double, const double, struct prjprm *, double *, double *);
int cscset(struct prjprm *);
int cscfwd(const double, const double, struct prjprm *, double *, double *);
int cscrev(const double, const double, struct prjprm *, double *, double *);
int qscset(struct prjprm *);
int qscfwd(const double, const double, struct prjprm *, double *, double *);
int qscrev(const double, const double, struct prjprm *, double *, double *);
int raw_to_pv(struct prjprm *prj, double x, double y, double *xo, double *yo);
#else
int prjset(), prjfwd(), prjrev();
int azpset(), azpfwd(), azprev();
int szpset(), szpfwd(), szprev();
int tanset(), tanfwd(), tanrev();
int stgset(), stgfwd(), stgrev();
int sinset(), sinfwd(), sinrev();
int arcset(), arcfwd(), arcrev();
int zpnset(), zpnfwd(), zpnrev();
int zeaset(), zeafwd(), zearev();
int airset(), airfwd(), airrev();
int cypset(), cypfwd(), cyprev();
int ceaset(), ceafwd(), cearev();
int carset(), carfwd(), carrev();
int merset(), merfwd(), merrev();
int sflset(), sflfwd(), sflrev();
int parset(), parfwd(), parrev();
int molset(), molfwd(), molrev();
int aitset(), aitfwd(), aitrev();
int copset(), copfwd(), coprev();
int coeset(), coefwd(), coerev();
int codset(), codfwd(), codrev();
int cooset(), coofwd(), coorev();
int bonset(), bonfwd(), bonrev();
int pcoset(), pcofwd(), pcorev();
int tscset(), tscfwd(), tscrev();
int cscset(), cscfwd(), cscrev();
int qscset(), qscfwd(), qscrev();
int raw_to_pv();
#endif
extern const char *prjset_errmsg[];
extern const char *prjfwd_errmsg[];
extern const char *prjrev_errmsg[];
#define PRJSET 137
struct celprm {
int flag;
double ref[4];
double euler[5];
};
#if __STDC__ || defined(__cplusplus)
int celset(const char *, struct celprm *, struct prjprm *);
int celfwd(const char *,
const double, const double,
struct celprm *,
double *, double *,
struct prjprm *,
double *, double *);
int celrev(const char *,
const double, const double,
struct prjprm *,
double *, double *,
struct celprm *,
double *, double *);
#else
int celset(), celfwd(), celrev();
#endif
extern const char *celset_errmsg[];
extern const char *celfwd_errmsg[];
extern const char *celrev_errmsg[];
#define CELSET 137
struct linprm {
int flag;
int naxis;
double *crpix;
double *pc;
double *cdelt;
/* Intermediates. */
double *piximg;
double *imgpix;
};
#if __STDC__ || defined(__cplusplus)
int linset(struct linprm *);
int linfwd(const double[], struct linprm *, double[]);
int linrev(const double[], struct linprm *, double[]);
int matinv(const int, const double [], double []);
#else
int linset(), linfwd(), linrev(), matinv();
#endif
extern const char *linset_errmsg[];
extern const char *linfwd_errmsg[];
extern const char *linrev_errmsg[];
#define LINSET 137
struct wcsprm {
int flag;
char pcode[4];
char lngtyp[5], lattyp[5];
int lng, lat;
int cubeface;
};
#if __STDC__ || defined(__cplusplus)
int wcssett(const int,
const char[][9],
struct wcsprm *);
int wcsfwd(const char[][9],
struct wcsprm *,
const double[],
const double[],
struct celprm *,
double *,
double *,
struct prjprm *,
double[],
struct linprm *,
double[]);
int wcsrevv(const char[][9],
struct wcsprm *,
const double[],
struct linprm *,
double[],
struct prjprm *,
double *,
double *,
const double[],
struct celprm *,
double[]);
int wcsmix(const char[][9],
struct wcsprm *,
const int,
const int,
const double[],
const double,
int,
double[],
const double[],
struct celprm *,
double *,
double *,
struct prjprm *,
double[],
struct linprm *,
double[]);
#else
int wcsset(), wcsfwd(), wcsrev(), wcsmix();
#endif
extern const char *wcsset_errmsg[];
extern const char *wcsfwd_errmsg[];
extern const char *wcsrev_errmsg[];
extern const char *wcsmix_errmsg[];
#define WCSSET 137
#if __STDC__ || defined(__cplusplus)
int sphfwd(const double, const double,
const double [],
double *, double *);
int sphrev(const double, const double,
const double [],
double *, double *);
#else
int sphfwd(), sphrev();
#endif
#ifdef PI
#undef PI
#endif
#ifdef D2R
#undef D2R
#endif
#ifdef R2D
#undef R2D
#endif
#ifdef SQRT2
#undef SQRT2
#endif
#ifdef SQRT2INV
#undef SQRT2INV
#endif
#ifdef D2S
#undef D2S
#endif
#ifdef S2D
#undef S2D
#endif
#define PI 3.141592653589793238462643
#define D2R PI/180.0
#define R2D 180.0/PI
#define S2D 1.0/3600.0
#define D2S 3600.0
#define SQRT2 1.4142135623730950488
#define SQRT2INV 1.0/SQRT2
#if !defined(__STDC__) && !defined(__cplusplus)
#ifndef const
#define const
#endif
#endif
#if __STDC__ || defined(__cplusplus)
double cosdeg(const double);
double sindeg(const double);
double tandeg(const double);
double acosdeg(const double);
double asindeg(const double);
double atandeg(const double);
double atan2deg(const double, const double);
#else
double cosdeg();
double sindeg();
double tandeg();
double acosdeg();
double asindeg();
double atandeg();
double atan2deg();
#endif
/* Domain tolerance for asin and acos functions. */
#define WCSTRIG_TOL 1e-10
#ifdef __cplusplus
}
#endif
#endif /* wcslib_h_ */
/* Feb 3 2000 Doug Mink - Make cplusplus ifdefs for braces all-inclusive
*
* Feb 15 2001 Doug Mink - Undefine math constants if already defined
* Sep 19 2001 Doug Mink - Update for WCSLIB 2.7, especially proj.h and cel.h
*
* Mar 12 2002 Doug Mink - Update for WCSLIB 2.8.2, especially proj.h
* Nov 29 2006 Doug Mink - Drop semicolon at end of C++ ifdef
* Jan 4 2007 Doug Mink - Drop extra declarations of SZP subroutines
*
* Mar 30 2011 Doug Mink - Add raw_to_pv() subroutine for SCAMP from Ed Los
*/
|
