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