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
30 static int get_param_int(char *buf)
32 /* Returns an integer parameter from a header line. */
37 substr = index(buf, ':');
38 sscanf(substr, ": %d", &value);
42 static float get_param_float(char *buf)
44 /* Returns a floating point parameter from a header line. */
49 substr = index(buf, ':');
50 sscanf(substr, ": %f", &value);
55 static char *get_param_string(char *buf)
57 /* Returns a string parameter from a header line. */
59 static char string[20];
62 substr = index(buf, ':');
63 sscanf(substr, ": %s", string);
67 void A_label(Rainbow_hdr *rainbow_header, char* buf)
72 sscanf(buf, "%c%d", &label, &labelnum);
75 rainbow_header->az_step = get_param_float(buf);
76 else if (labelnum == 9)
77 rainbow_header->nsweeps = get_param_int(buf);
80 void F_label(Rainbow_hdr *rainbow_header, char* buf)
82 int labelnum, day, month, year, sec, minute, hour;
85 sscanf(buf, "%*c%d", &labelnum);
89 rainbow_header->compressed = get_param_int(buf);
92 sscanf(buf, "%*c%*d : %f %f", &lon, &lat);
93 rainbow_header->lon = lon;
94 rainbow_header->lat = lat;
97 sscanf(buf, "%*c%*d : %d %d %d", &day, &month, &year);
98 rainbow_header->month = month;
99 rainbow_header->day = day;
100 rainbow_header->year = year;
103 sscanf(buf, "%*c%*d : %d %d %d", &sec, &minute, &hour);
104 rainbow_header->hour = hour;
105 rainbow_header->minute = minute;
106 rainbow_header->sec = sec;
109 rainbow_header->datatype = get_param_int(buf);
114 void H_label(Rainbow_hdr *rainbow_header, char* buf)
119 sscanf(buf, "%c%d", &label, &labelnum);
122 rainbow_header->filetype = get_param_int(buf);
123 else if (labelnum == 8)
124 strcpy(rainbow_header->radarname, get_param_string(buf));
127 void P_label(Rainbow_hdr *rainbow_header, char* buf)
132 sscanf(buf, "%c%d", &label, &labelnum);
135 rainbow_header->range_start= get_param_float(buf);
136 else if (labelnum == 4)
137 rainbow_header->range_stop = get_param_float(buf);
138 else if (labelnum == 5)
139 rainbow_header->range_step = get_param_float(buf);
142 void R_label(Rainbow_hdr *rainbow_header, char* buf)
146 sscanf(buf, "%*c%d", &labelnum);
149 rainbow_header->nbins = get_param_int(buf);
150 else if (labelnum == 2)
151 rainbow_header->bin_resolution = get_param_float(buf);
152 else if (labelnum == 8)
153 rainbow_header->nvalues = get_param_int(buf);
156 void W_label(Rainbow_hdr *rainbow_hdr, char* buf)
158 int labelnum, az_start, az_stop, pw_code, prf_high, prf_low, zdata;
159 int vdata, wdata, unfolding, cdata, ddata, uzdata;
160 float elev, az_step, az_rate, range_stop;
163 sscanf(buf, "%*c%d : %f %d %d %f %f %d %d %d %d %d %d %d %f %d %d %d",
164 &labelnum, &elev, &az_start, &az_stop, &az_step,
165 &az_rate, &pw_code, &prf_high, &prf_low, &zdata, &vdata,
166 &wdata, &unfolding, &range_stop, &cdata, &ddata, &uzdata);
168 /* Note: Only need to collect parameters 1, 5, 7, 8, and 13 for each
169 * elevation. Parameters 2, 3, and 4 are fixed at 0, 359, and 1 for volume
170 * scan. The remaining parameters can be taken once from label number 1
171 * (first elevation). Also, don't need az_step; got it from A3.
175 rainbow_hdr->az_start = az_start;
176 rainbow_hdr->az_stop = az_stop;
177 rainbow_hdr->pulse_width_code = pw_code;
178 rainbow_hdr->zdata = zdata;
179 rainbow_hdr->vdata = vdata;
180 rainbow_hdr->wdata = wdata;
181 rainbow_hdr->unfolding = unfolding;
182 rainbow_hdr->cdata = cdata;
183 rainbow_hdr->ddata = ddata;
184 rainbow_hdr->uzdata = uzdata;
186 rainbow_hdr->elev_params[labelnum-1] =
187 (struct elev_params *) malloc(sizeof(struct elev_params));
188 rainbow_hdr->elev_params[labelnum-1]->elev_angle = elev;
189 rainbow_hdr->elev_params[labelnum-1]->az_rate = az_rate;
190 rainbow_hdr->elev_params[labelnum-1]->prf_high = prf_high;
191 rainbow_hdr->elev_params[labelnum-1]->prf_low = prf_low;
192 rainbow_hdr->elev_params[labelnum-1]->maxrange = range_stop;
195 /**********************************************************/
199 /**********************************************************/
201 static int read_hdr_line(char *buf, int maxchars, FILE *fp)
203 /* Read a line from the Rainbow file header into character buffer.
204 * Function returns the first character in buffer (the "label") if
205 * the line was successfully read, -1 otherwise.
207 * Note: the following control characters (defined in rainbow.h) are used
208 * in the Rainbow header:
209 * CR - Carriage Return: end of a line.
210 * ETB - End of Transmission Block: divides header into sections (we can
212 * ETX - End of Text: end of header.
219 while ((c = getc(fp)) != CR && c != ETX) {
221 c = getc(fp); /* Read past both <ETB> and the */
222 if (c == CR) c = getc(fp); /* combination <ETB><CR>. */
234 fprintf(stderr,"A header line exceeded buffer size %d.\n",maxchars);
235 fprintf(stderr,"Did not find end-of-line character 0x%02x.\n",CR);
236 buf[maxchars - 1] = '\0'; /* Make it a legal string anyway. */
242 if (c != ETX) c = (int) buf[0];
246 /**********************************************************/
248 /* read_rainbow_header */
250 /**********************************************************/
254 void read_rainbow_header(Rainbow_hdr *rainbow_header, FILE *fp)
256 /* Reads parameters from Rainbow file header into a rainbow header
262 /* Read each line of the header and extract parameters according to the
263 * label. The label is a single alphabetic character at the beginning
264 * of the line which indicates a category of parameters.
267 while ((label = read_hdr_line(buf, BUFSIZE, fp)) != ETX && label > 0) {
269 case 'H': H_label(rainbow_header, buf);
271 case 'P': P_label(rainbow_header, buf);
273 case 'W': W_label(rainbow_header, buf);
275 case 'F': F_label(rainbow_header, buf);
277 case 'R': R_label(rainbow_header, buf);
279 case 'A': A_label(rainbow_header, buf);
281 case 'N': rainbow_header->product = get_param_int(buf);
283 case 'X': /* X_label(rainbow_header, buf); */