Pileus Git - ~andy/rsl/blob - examples/cappi_image.c

   1 /*
   2  * For RSL version 0.28 and higher.
   3  *
   4  * v1.0 Began 2/16/94 by John Merritt.
   5  *
   6  * Demonstrates reading NEXRAD files and loading the Radar structure.
   7  */
   8
   9 #include <stdio.h>
  10 #include <stdlib.h>
  11 #include <string.h>
  12 #include "rsl.h"
  13
  14 usage(int argc, char **argv)
  15 {
  16   fprintf(stderr,"Usage: %s infile [firstfile | callid]\n", argv[0]);
  17   exit(-1);
  18 }
  19
  20 process_args(int argc, char **argv,
  21                          char **in_file, char **out_file)
  22 {
  23
  24   if (argc == 2) {
  25         *in_file = strdup(argv[1]);
  26         *out_file = NULL;
  27   } else if (argc == 3) {
  28         *in_file = strdup(argv[1]);
  29         *out_file = strdup(argv[2]);
  30   } else {
  31         usage(argc, argv);
  32   }
  33 }
  34
  35
  36 main(int argc, char **argv)
  37 {
  38   char *infile, *firstfile;
  39   char outfile[100];
  40
  41   Radar *radar;
  42   Cappi *cappi;
  43
  44   int i, j;
  45
  46 /* 1. Process the arguments. */
  47   process_args(argc, argv, &infile, &firstfile); /* malloc for in/outfile */
  48
  49
  50   i = DZ_INDEX;
  51
  52   RSL_radar_verbose_on();
  53   if ((radar = RSL_anyformat_to_radar(infile, firstfile)) == NULL) exit(-1);
  54
  55   printf("Radar date: %2.2d/%2.2d/%2.2d\n", radar->h.month, radar->h.day, radar->h.year);
  56   printf("Radar time: %2.2d:%2.2d:%f\n", radar->h.hour, radar->h.minute, radar->h.sec);
  57   printf("Radar file: %s\n", infile);
  58   printf("Radar site: %c%c%c%c\n",
  59                  radar->h.name[0],
  60                  radar->h.name[1],
  61                  radar->h.name[2],
  62                  radar->h.name[3]);
  63   printf("Radar date: %2.2d/%2.2d/%2.2d\n", radar->h.month, radar->h.day, radar->h.year);
  64   printf("Radar time: %2.2d:%2.2d:%f\n", radar->h.hour, radar->h.minute, radar->h.sec);
  65
  66
  67   if (i == DZ_INDEX) RSL_load_refl_color_table();
  68   if (i == VR_INDEX) RSL_load_vel_color_table();
  69   if (i == SW_INDEX) RSL_load_sw_color_table();
  70   if (i == VR_INDEX) RSL_rebin_velocity_volume(radar->v[i]); /* Modifies v[i]. */
  71
  72   for (j=1; j<=5; j++) {
  73         cappi = RSL_cappi_at_h( radar -> v[i], (float)j/2.0, 200.0);
  74         sprintf(outfile, "cappi_%2.2d.gif", j);
  75         printf("Writing %s.\n", outfile);
  76         RSL_sweep_to_gif(cappi->sweep, outfile, 400, 400, 200.0);
  77         RSL_free_cappi(cappi);
  78   }
  79
  80
  81   exit(0);
  82
  83 }