Changeset 14574
- Timestamp:
- Aug 20, 2007, 3:44:48 PM (19 years ago)
- File:
-
- 1 edited
-
trunk/Ohana/src/getstar/src/ReadImageHeader.c (modified) (4 diffs)
Legend:
- Unmodified
- Added
- Removed
-
trunk/Ohana/src/getstar/src/ReadImageHeader.c
r14496 r14574 4 4 int ReadImageHeader (Header *header, Image *image) { 5 5 6 int Nastro, ccdnum, hour, min, Nx, Ny , photcode;6 int Nastro, ccdnum, hour, min, Nx, Ny; 7 7 double tmp, sec, Cerror; 8 8 char *c, photname[64], line[80]; … … 66 66 image[0].cerror = tmp * 50.0; 67 67 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 these82 // gfits_scan (header, ExptimeKeyword, "%lf", 1, &tmp);83 // image[0].exptime = tmp;84 85 68 /*** why are we no longer using APMIFIT?? ***/ 86 69 tmp = 0; … … 108 91 image[0].fwhm_y = tmp * 25.0 * image[0].coords.cdelt1 * 3600.0; 109 92 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; 143 96 144 97 /* secz is in units milli-airmass */ … … 148 101 bzero (image[0].dummy, sizeof(image[0].dummy)); 149 102 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 156 103 return (TRUE); 157 104 }
Note:
See TracChangeset
for help on using the changeset viewer.
