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
|
// Wrapping of WCSTools wcs.c
//
// Not wrapping the functions that return WCS coords as string because can convert using decimal2hms() etc. in astCoords
%feature("autodoc", "0");
%module (package="PyWCSTools") wcs
%include "carrays.i"
%include "typemaps.i"
%array_functions(double, doubleArray)
/* %feature("autodoc", "1") */
%{
#include <wcstools/wcs.h>
%}
/* WORKING */
struct WorldCoor *wcsinit(char *hstring); /* set up a WCS structure from a FITS image header */
struct WorldCoor *wcsxinit(double cra, double cdec, double secpix, int xrpix, int yrpix, int nxpix, int nypix, double rotate, int equinox, double epoch, char *proj); /* set up a WCS structure from arguments */
struct WorldCoor *wcskinit(int nxpix, int nypix, char *ctype1, char *ctype2, double crpix1, double crpix2, double crval1, double crval2, double *cd, double cdelt1, double cdelt2, double crota, double equinox, double epoch); /* set up a WCS structure from keyword values */
int iswcs(struct WorldCoor *wcs); /* Return 1 if WCS structure is filled, else 0 */
int nowcs(struct WorldCoor *wcs); /* Return 0 if WCS structure is filled, else 1 */
%apply double *OUTPUT { double *xpix, double *ypix, int *offscl};
void wcs2pix(struct WorldCoor *wcs, double xpos, double ypos, double *xpix, double *ypix, int *offscl); /* Convert World Coordinates to pixel coordinates */
%apply double *OUTPUT { double *xpos, double *ypos};
void pix2wcs(struct WorldCoor *wcs, double xpix, double ypix, double *xpos, double *ypos); /* Convert pixel coordinates to World Coordinates */
void wcscent(struct WorldCoor *wcs); /* Print the image center and size in WCS units */
char *getradecsys(struct WorldCoor *wcs); /* Return current value of coordinate system */
void wcsoutinit(struct WorldCoor *wcs, char *coorsys); /* Initialize WCS output coordinate system for use by pix2wcs */
void wcsininit(struct WorldCoor *wcs, char *coorsys); /* Initialize WCS input coordinate system for use by wcs2pix */
char *getwcsout(struct WorldCoor *wcs); /* Return current value of WCS output coordinate system */
char *getwcsin(struct WorldCoor *wcs); /* Return current value of WCS input coordinate system */
%apply double *OUTPUT {double *cra, double *cdec, double *dra, double *ddec};
void wcssize(struct WorldCoor *wcs, double *cra, double *cdec, double *dra, double *ddec); /* Return RA and Dec of image center, size in RA and Dec */
%apply double *OUTPUT {double *width, double *height};
void wcsfull(struct WorldCoor *wcs, double *cra, double *cdec, double *width, double *height); /* Return RA and Dec of image center, size in degrees */
struct WorldCoor {
double xref; /* X reference coordinate value (deg) */
double yref; /* Y reference coordinate value (deg) */
double xrefpix; /* X reference pixel */
double yrefpix; /* Y reference pixel */
double xinc; /* X coordinate increment (deg) */
double yinc; /* Y coordinate increment (deg) */
double rot; /* rotation around axis (deg) (N through E) */
double cd[4]; /* rotation matrix */
double dc[4]; /* inverse rotation matrix */
double equinox; /* Equinox of coordinates default to 1950.0 */
double epoch; /* Epoch of coordinates default to equinox */
double nxpix; /* Number of pixels in X-dimension of image */
double nypix; /* Number of pixels in Y-dimension of image */
double plate_ra; /* Right ascension of plate center */
double plate_dec; /* Declination of plate center */
double plate_scale; /* Plate scale in arcsec/mm */
double x_pixel_offset; /* X pixel offset of image lower right */
double y_pixel_offset; /* Y pixel offset of image lower right */
double x_pixel_size; /* X pixel_size */
double y_pixel_size; /* Y pixel_size */
double ppo_coeff[6]; /* pixel to plate coefficients for DSS */
double x_coeff[20]; /* X coefficients for plate model */
double y_coeff[20]; /* Y coefficients for plate model */
double xpix; /* X (RA) coordinate (pixels) */
double ypix; /* Y (dec) coordinate (pixels) */
double zpix; /* Z (face) coordinate (pixels) */
double xpos; /* X (RA) coordinate (deg) */
double ypos; /* Y (dec) coordinate (deg) */
double crpix[9]; /* Values of CRPIXn keywords */
double crval[9]; /* Values of CRVALn keywords */
double cdelt[9]; /* Values of CDELTn keywords */
double pc[81]; /* Values of PCiiijjj keywords */
double projp[10]; /* Constants for various projections */
int pvfail; /* If non-zero, significant inaccuracy likely to occur in projection */
double projppv[2*MAXPV]; /* SCAMP constants for the PV coordinates */
struct poly *inv_x; /* SCAMP projection correction polynom in x */
struct poly *inv_y; /* SCAMP projection correction polynom in y */
double longpole; /* Longitude of North Pole in degrees */
double latpole; /* Latitude of North Pole in degrees */
double rodeg; /* Radius of the projection generating sphere */
double imrot; /* Rotation angle of north pole */
double pa_north; /* Position angle of north (0=horizontal) */
double pa_east; /* Position angle of east (0=horizontal) */
double radvel; /* Radial velocity (km/sec away from observer)*/
double zvel; /* Radial velocity (v/c away from observer)*/
double zpzd; /* Colat of FIP (degs) */
double zpr; /* Radius of FIP (degs) */
int imflip; /* If not 0, image is reflected around axis */
int prjcode; /* projection code (-1-32) */
int latbase; /* Latitude base 90 (NPA), 0 (LAT), -90 (SPA) */
int ncoeff1; /* Number of x-axis plate fit coefficients */
int ncoeff2; /* Number of y-axis plate fit coefficients */
int zpnp; /* ZP polynomial order (0-9) */
int changesys; /* 1 for FK4->FK5, 2 for FK5->FK4 */
/* 3 for FK4->galactic, 4 for FK5->galactic */
int printsys; /* 1 to print coordinate system, else 0 */
int ndec; /* Number of decimal places in PIX2WCST */
int degout; /* 1 to always print degrees in PIX2WCST */
int tabsys; /* 1 to put tab between RA & Dec, else 0 */
int rotmat; /* 0 if CDELT, CROTA; 1 if CD */
int coorflip; /* 0 if x=RA, y=Dec; 1 if x=Dec, y=RA */
int offscl; /* 0 if OK, 1 if offscale */
int wcson; /* 1 if WCS is set, else 0 */
int naxis; /* Number of axes in image (for WCSLIB 3.0) */
int naxes; /* Number of axes in image */
int wcsproj; /* WCS_OLD: AIPS worldpos() and worldpix()
WCS_NEW: Mark Calabretta's WCSLIB subroutines
WCS_BEST: WCSLIB for all but CAR,COE,NCP
WCS_ALT: AIPS for all but CAR,COE,NCP */
int linmode; /* 0=system only, 1=units, 2=system+units */
int detector; /* Instrument detector number */
char instrument[32]; /* Instrument name */
char ctype[9][9]; /* Values of CTYPEn keywords */
char c1type[9]; /* 1st coordinate type code:
RA--, GLON, ELON */
char c2type[9]; /* 2nd coordinate type code:
DEC-, GLAT, ELAT */
char ptype[9]; /* projection type code:
SIN, TAN, ARC, NCP, GLS, MER, AIT, etc */
char units[9][32]; /* Units if LINEAR */
char radecsys[32]; /* Reference frame: FK4, FK4-NO-E, FK5, GAPPT*/
char radecout[32]; /* Output reference frame: FK4,FK5,GAL,ECL */
char radecin[32]; /* Input reference frame: FK4,FK5,GAL,ECL */
double eqin; /* Input equinox (match sysin if 0.0) */
double eqout; /* Output equinox (match sysout if 0.0) */
int sysin; /* Input coordinate system code */
int syswcs; /* WCS coordinate system code */
int sysout; /* Output coordinate system code */
/* WCS_B1950, WCS_J2000, WCS_ICRS, WCS_GALACTIC,
* WCS_ECLIPTIC, WCS_LINEAR, WCS_ALTAZ */
char center[32]; /* Center coordinates (with frame) */
struct wcsprm wcsl; /* WCSLIB main projection parameters */
struct linprm lin; /* WCSLIB image/pixel conversion parameters */
struct celprm cel; /* WCSLIB projection type */
struct prjprm prj; /* WCSLIB projection parameters */
struct IRAFsurface *lngcor; /* RA/longitude correction structure */
struct IRAFsurface *latcor; /* Dec/latitude correction structure */
int distcode; /* Distortion code 0=none 1=SIRTF */
struct Distort distort; /* SIRTF distortion coefficients */
char *command_format[10]; /* WCS command formats */
/* where %s is replaced by WCS coordinates */
/* where %f is replaced by the image filename */
/* where %x is replaced by image coordinates */
double ltm[4]; /* Image rotation matrix */
double ltv[2]; /* Image offset */
int idpix[2]; /* First pixel to use in image (x, y) */
int ndpix[2]; /* Number of pixels to use in image (x, y) */
struct WorldCoor *wcs; /* WCS upon which this WCS depends */
struct WorldCoor *wcsdep; /* WCS depending on this WCS */
char *wcsname; /* WCS name (defaults to NULL pointer) */
char wcschar; /* WCS character (A-Z, null, space) */
int logwcs; /* 1 if DC-FLAG is set for log wavelength */
};
|
