3 This is the TRMM Office Radar Software Library.
6 George Mason University
9 This library is free software; you can redistribute it and/or
10 modify it under the terms of the GNU Library General Public
11 License as published by the Free Software Foundation; either
12 version 2 of the License, or (at your option) any later version.
14 This library is distributed in the hope that it will be useful,
15 but WITHOUT ANY WARRANTY; without even the implied warranty of
16 MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
17 Library General Public License for more details.
19 You should have received a copy of the GNU Library General Public
20 License along with this library; if not, write to the Free Software
21 Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA 02111-1307 USA
29 static int get_param_int(char *buf)
31 /* Returns an integer parameter from a header line. */
36 substr = index(buf, ':');
37 sscanf(substr, ": %d", &value);
41 static float get_param_float(char *buf)
43 /* Returns a floating point parameter from a header line. */
48 substr = index(buf, ':');
49 sscanf(substr, ": %f", &value);
54 static char *get_param_string(char *buf)
56 /* Returns a string parameter from a header line. */
58 static char string[20];
61 substr = index(buf, ':');
62 sscanf(substr, ": %s", string);
66 void A_label(Rainbow_hdr *rainbow_header, char* buf)
71 sscanf(buf, "%c%d", &label, &labelnum);
74 rainbow_header->az_step = get_param_float(buf);
75 else if (labelnum == 9)
76 rainbow_header->nsweeps = get_param_int(buf);
79 void F_label(Rainbow_hdr *rainbow_header, char* buf)
81 int labelnum, day, month, year, sec, minute, hour;
84 sscanf(buf, "%*c%d", &labelnum);
88 rainbow_header->compressed = get_param_int(buf);
91 sscanf(buf, "%*c%*d : %f %f", &lon, &lat);
92 rainbow_header->lon = lon;
93 rainbow_header->lat = lat;
96 sscanf(buf, "%*c%*d : %d %d %d", &day, &month, &year);
97 rainbow_header->month = month;
98 rainbow_header->day = day;
99 rainbow_header->year = year;
102 sscanf(buf, "%*c%*d : %d %d %d", &sec, &minute, &hour);
103 rainbow_header->hour = hour;
104 rainbow_header->minute = minute;
105 rainbow_header->sec = sec;
108 rainbow_header->datatype = get_param_int(buf);
113 void H_label(Rainbow_hdr *rainbow_header, char* buf)
118 sscanf(buf, "%c%d", &label, &labelnum);
121 rainbow_header->filetype = get_param_int(buf);
122 else if (labelnum == 8)
123 strcpy(rainbow_header->radarname, get_param_string(buf));
126 void P_label(Rainbow_hdr *rainbow_header, char* buf)
131 sscanf(buf, "%c%d", &label, &labelnum);
134 rainbow_header->range_start= get_param_float(buf);
135 else if (labelnum == 4)
136 rainbow_header->range_stop = get_param_float(buf);
137 else if (labelnum == 5)
138 rainbow_header->range_step = get_param_float(buf);
141 void R_label(Rainbow_hdr *rainbow_header, char* buf)
145 sscanf(buf, "%*c%d", &labelnum);
148 rainbow_header->nbins = get_param_int(buf);
149 else if (labelnum == 2)
150 rainbow_header->bin_resolution = get_param_float(buf);
151 else if (labelnum == 8)
152 rainbow_header->nvalues = get_param_int(buf);
155 void W_label(Rainbow_hdr *rainbow_hdr, char* buf)
157 int labelnum, az_start, az_stop, pw_code, prf_high, prf_low, zdata;
158 int vdata, wdata, unfolding, cdata, ddata, uzdata;
159 float elev, az_step, az_rate, range_stop;
162 sscanf(buf, "%*c%d : %f %d %d %f %f %d %d %d %d %d %d %d %f %d %d %d",
163 &labelnum, &elev, &az_start, &az_stop, &az_step,
164 &az_rate, &pw_code, &prf_high, &prf_low, &zdata, &vdata,
165 &wdata, &unfolding, &range_stop, &cdata, &ddata, &uzdata);
167 /* Note: Only need to collect parameters 1, 5, 7, 8, and 13 for each
168 * elevation. Parameters 2, 3, and 4 are fixed at 0, 359, and 1 for volume
169 * scan. The remaining parameters can be taken once from label number 1
170 * (first elevation). Also, don't need az_step; got it from A3.
174 rainbow_hdr->az_start = az_start;
175 rainbow_hdr->az_stop = az_stop;
176 rainbow_hdr->pulse_width_code = pw_code;
177 rainbow_hdr->zdata = zdata;
178 rainbow_hdr->vdata = vdata;
179 rainbow_hdr->wdata = wdata;
180 rainbow_hdr->unfolding = unfolding;
181 rainbow_hdr->cdata = cdata;
182 rainbow_hdr->ddata = ddata;
183 rainbow_hdr->uzdata = uzdata;
185 rainbow_hdr->elev_params[labelnum-1] =
186 (struct elev_params *) malloc(sizeof(struct elev_params));
187 rainbow_hdr->elev_params[labelnum-1]->elev_angle = elev;
188 rainbow_hdr->elev_params[labelnum-1]->az_rate = az_rate;
189 rainbow_hdr->elev_params[labelnum-1]->prf_high = prf_high;
190 rainbow_hdr->elev_params[labelnum-1]->prf_low = prf_low;
191 rainbow_hdr->elev_params[labelnum-1]->maxrange = range_stop;
194 /**********************************************************/
198 /**********************************************************/
200 static int read_hdr_line(char *buf, int maxchars, FILE *fp)
202 /* Read a line from the Rainbow file header into character buffer.
203 * Function returns the first character in buffer (the "label") if
204 * the line was successfully read, -1 otherwise.
206 * Note: the following control characters (defined in rainbow.h) are used
207 * in the Rainbow header:
208 * CR - Carriage Return: end of a line.
209 * ETB - End of Transmission Block: divides header into sections (we can
211 * ETX - End of Text: end of header.
218 while ((c = getc(fp)) != CR && c != ETX) {
220 c = getc(fp); /* Read past both <ETB> and the */
221 if (c == CR) c = getc(fp); /* combination <ETB><CR>. */
233 fprintf(stderr,"A header line exceeded buffer size %d.\n",maxchars);
234 fprintf(stderr,"Did not find end-of-line character 0x%02x.\n",CR);
235 buf[maxchars - 1] = '\0'; /* Make it a legal string anyway. */
241 if (c != ETX) c = (int) buf[0];
245 /**********************************************************/
247 /* read_rainbow_header */
249 /**********************************************************/
253 void read_rainbow_header(Rainbow_hdr *rainbow_header, FILE *fp)
255 /* Reads parameters from Rainbow file header into a rainbow header
261 /* Read each line of the header and extract parameters according to the
262 * label. The label is a single alphabetic character at the beginning
263 * of the line which indicates a category of parameters.
266 while ((label = read_hdr_line(buf, BUFSIZE, fp)) != ETX && label > 0) {
268 case 'H': H_label(rainbow_header, buf);
270 case 'P': P_label(rainbow_header, buf);
272 case 'W': W_label(rainbow_header, buf);
274 case 'F': F_label(rainbow_header, buf);
276 case 'R': R_label(rainbow_header, buf);
278 case 'A': A_label(rainbow_header, buf);
280 case 'N': rainbow_header->product = get_param_int(buf);
282 case 'X': /* X_label(rainbow_header, buf); */