IPP Software Navigation Tools IPP Links Communication Pan-STARRS Links

Changeset 14574


Ignore:
Timestamp:
Aug 20, 2007, 3:44:48 PM (19 years ago)
Author:
eugene
Message:

do not require photcode or nstars

File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/Ohana/src/getstar/src/ReadImageHeader.c

    r14496 r14574  
    44int ReadImageHeader (Header *header, Image *image) {
    55
    6   int Nastro, ccdnum, hour, min, Nx, Ny, photcode;
     6  int Nastro, ccdnum, hour, min, Nx, Ny;
    77  double tmp, sec, Cerror;
    88  char *c, photname[64], line[80];
     
    6666  image[0].cerror = tmp * 50.0;
    6767 
    68   /* get photcode from header */
    69   /*** XXX is photcode needed?? ***/
    70   if (!gfits_scan (header, "PHOTCODE", "%s", 1, photname)) {
    71     fprintf (stderr, "photcode not supplied in header\n");
    72     return (FALSE);
    73   }
    74   photcode = GetPhotcodeCodebyName (photname);
    75   if (photcode == 0) {
    76     fprintf (stderr, "no valid photcode is supplied\n");
    77     return (FALSE);
    78   }
    79   image[0].photcode = photcode;
    80 
    81   // XXX we don't really need these
    82   // gfits_scan (header, ExptimeKeyword,  "%lf", 1, &tmp);
    83   // image[0].exptime = tmp;
    84  
    8568  /*** why are we no longer using APMIFIT?? ***/
    8669  tmp = 0;
     
    10891  image[0].fwhm_y = tmp * 25.0 * image[0].coords.cdelt1 * 3600.0;
    10992
    110   // XXX we don't really need the following values
    111 # if (0) 
    112   if (STKeyword[0]) {
    113     /* get ST (used for airmass calculation) */
    114     gfits_scan (header, STKeyword, "%s", 1, line);
    115     /* remove ':' characters */
    116     for (c = strchr (line, ':'); c != (char *) NULL; c = strchr (line, ':')) { *c = ' '; }
    117     sscanf (line, "%d %d %lf", &hour, &min, &sec);
    118     image[0].sidtime = hour + min/60.0 + sec/3600.0;
    119   } else {
    120     double jd;
    121     jd = sec_to_jd (image[0].tzero);
    122     image[0].sidtime  = ohana_lst (jd, Longitude);
    123   }
    124   image[0].latitude = Latitude;
    125 
    126   tmp = 0;
    127   gfits_scan (header, "TRATE",   "%lf", 1, &tmp);
    128   image[0].trate = 10000 * tmp;
    129 
    130   tmp = 0;
    131   gfits_scan (header, AirmassKeyword, "%lf", 1, &tmp);
    132   image[0].secz = MIN (NO_MAG, tmp);
    133 
    134   if (!gfits_scan (header, CCDNumKeyword, "%d", 1, &ccdnum)) {
    135     image[0].ccdnum = 0xff;
    136   } else {
    137     image[0].ccdnum = ccdnum;
    138   }
    139 
    140   gfits_scan (header, "ZERO_PT", "%lf", 1, &ZeroPt);
    141   SetZeroPoint (ZeroPt);
    142 # endif
     93  image[0].trate = 0.0;
     94  image[0].secz = NO_MAG;
     95  image[0].ccdnum = 0;
    14396
    14497  /* secz is in units milli-airmass */
     
    148101  bzero (image[0].dummy, sizeof(image[0].dummy));
    149102
    150   /* find expected number of stars */
    151   if (!gfits_scan (header, "NSTARS", "%d", 1, &image[0].nstar)) {
    152     fprintf (stderr, "can't get NSTARS from header\n");
    153     return (FALSE);
    154   }
    155 
    156103  return (TRUE);
    157104}
Note: See TracChangeset for help on using the changeset viewer.