int fits_update_key_from_str(fitsfile* fptr, int datatype, char* name, char* valstr, char* comment, int* status) {
if (*status > 0) return *status;
if (datatype == UNDEFINED_TYPE)
datatype = fits_detect_key_type(valstr);
float cf[2];
double cd[2];
char keyval[FLEN_VALUE];
switch(datatype) {
case TSTRING:
removequotes(valstr, keyval);
fits_update_key_str(fptr, name, keyval, comment, status);
break;
case TBYTE: case TSBYTE: case TUSHORT: case TSHORT: case TINT: case TLONG:
fits_update_key_lng(fptr, name, atol(valstr), comment, status);
break;
case TUINT: case TULONG:
fits_update_key_lng(fptr, name, atof(valstr), comment, status);
break;
case TLOGICAL:
fits_update_key_log(fptr, name, strcasecmp(valstr,"TRUE") == 0 ? TRUE: FALSE , comment, status);
break;
case TFLOAT:
fits_update_key_flt(fptr, name, atof(valstr), -7, comment, status);
break;
case TDOUBLE:
fits_update_key_dbl(fptr, name, atof(valstr), -15, comment, status);
break;
/* For the complex numbers, we assume the string format: 1+3i or 1+3j */
case TCOMPLEX:
if (sscanf(valstr,"(%f,%f)", cf, cf+1) == 2)
fits_update_key_fixcmp(fptr, name, cf, -7, comment, status);
else
fprintf(stderr, " Error in %s : wrong complex format in %s", __func__, valstr);
break;
case TDBLCOMPLEX:
if (sscanf(valstr,"(%lf,%lf)", cd, cd+1) == 2)
fits_update_key_dblcmp(fptr, name, cd, -15, comment, status);
else
fprintf(stderr, " Error in %s : wrong complex format in %s", __func__, valstr);
break;
default: *status = BAD_DATATYPE;
}
return *status;
}
int HPXhdr(fitsfile *fptr, struct healpix *hpxdat)
{
char comment[64], cval[16], *ctype1, *ctype2, *descr1, *descr2, *pcode;
int status;
float cos45, crpix1, crpix2, crval1, crval2, lonpole;
double cdelt1, cdelt2;
status = 0;
fits_update_key_log(fptr, "EXTEND", 0,
"No FITS extensions are present", &status);
fits_write_date(fptr, &status);
/* Set pixel transformation parameters. */
if (hpxdat->layout == 0) {
crpix1 = (5 * hpxdat->nside + 1) / 2.0f;
} else {
crpix1 = (4 * hpxdat->nside + 1) / 2.0f;
}
crpix2 = crpix1;
fits_write_key(fptr, TFLOAT, "CRPIX1", &crpix1,
"Coordinate reference pixel", &status);
fits_write_key(fptr, TFLOAT, "CRPIX2", &crpix2,
"Coordinate reference pixel", &status);
cos45 = (float)sqrt(2.0) / 2.0f;
if (hpxdat->layout == 0) {
fits_write_key_flt(fptr, "PC1_1", cos45, -8,
"Transformation matrix element", &status);
fits_write_key_flt(fptr, "PC1_2", cos45, -8,
"Transformation matrix element", &status);
fits_write_key_flt(fptr, "PC2_1", -cos45, -8,
"Transformation matrix element", &status);
fits_write_key_flt(fptr, "PC2_2", cos45, -8,
"Transformation matrix element", &status);
}
cdelt1 = -90.0 / hpxdat->nside / sqrt(2.0);
cdelt2 = -cdelt1;
fits_write_key_dbl(fptr, "CDELT1", cdelt1, -8,
"[deg] Coordinate increment", &status);
fits_write_key_dbl(fptr, "CDELT2", cdelt2, -8,
"[deg] Coordinate increment", &status);
/* Celestial transformation parameters. */
if (hpxdat->layout == 0) {
pcode = "HPX";
} else {
pcode = "XPH";
}
if (hpxdat->crdsys == 'G') {
/* Galactic. */
ctype1 = "GLON";
ctype2 = "GLAT";
descr1 = "Galactic longitude";
descr2 = "Galactic latitude";
} else if (hpxdat->crdsys == 'E') {
/* Ecliptic, who-knows-what. */
ctype1 = "ELON";
ctype2 = "ELAT";
descr1 = "Ecliptic longitude";
descr2 = "Ecliptic latitude";
} else if (hpxdat->crdsys == 'Q') {
/* Equatorial, who-knows-what. */
ctype1 = "RA--";
ctype2 = "DEC-";
descr1 = "Right ascension";
descr2 = "Declination";
} else {
/* Unknown. */
ctype1 = "XLON";
ctype2 = "XLAT";
descr1 = "Longitude";
descr2 = " Latitude";
}
sprintf(cval, "%s-%s", ctype1, pcode);
sprintf(comment, "%s in an %s projection", descr1, pcode);
fits_write_key_str(fptr, "CTYPE1", cval, comment, &status);
sprintf(cval, "%s-%s", ctype2, pcode);
sprintf(comment, "%s in an %s projection", descr2, pcode);
fits_write_key_str(fptr, "CTYPE2", cval, comment, &status);
crval1 = 0.0f + 90.0f * hpxdat->quad;
if (hpxdat->layout == 0) {
crval2 = 0.0f;
} else if (hpxdat->layout == 1) {
crval1 += 180.0f;
crval2 = 90.0f;
} else {
crval1 += 180.0f;
crval2 = -90.0f;
}
if (360.0f < crval1) crval1 -= 360.0f;
sprintf(comment, "[deg] %s at the reference point", descr1);
fits_write_key(fptr, TFLOAT, "CRVAL1", &crval1, comment, &status);
sprintf(comment, "[deg] %s at the reference point", descr2);
fits_write_key(fptr, TFLOAT, "CRVAL2", &crval2, comment, &status);
if (hpxdat->layout) {
lonpole = 180.0f;
sprintf(comment, "[deg] Native longitude of the celestial pole");
fits_write_key(fptr, TFLOAT, "LONPOLE", &lonpole, comment, &status);
}
if (hpxdat->layout == 0) {
fits_write_key_lng(fptr, "PV2_1", (LONGLONG)4,
"HPX H parameter (longitude)", &status);
fits_write_key_lng(fptr, "PV2_2", (LONGLONG)3,
"HPX K parameter (latitude)", &status);
}
/* Commentary. */
fits_write_record(fptr, " ", &status);
if (hpxdat->layout == 0) {
fits_write_comment(fptr,
"Celestial map with FITS-standard HPX coordinate system generated by",
&status);
} else {
fits_write_comment(fptr,
"Celestial map with XPH coordinate system (polar HPX) generated by",
&status);
}
fits_write_comment(fptr,
"'HPXcvt' which reorganises HEALPix data without interpolation as",
&status);
fits_write_comment(fptr,
"described in \"Mapping on the HEALPix grid\" by Mark Calabretta and",
&status);
fits_write_comment(fptr,
"Boud Roukema. See http://www.atnf.csiro.au/people/Mark.Calabretta",
&status);
return status;
}
