FixNode/main.c

   1 // main.c -- read in a .node file and fix the z values of the interpolated
   2 //           points
   3 //
   4 // Written by Curtis Olson, started November 1997.
   5 //
   6 // Copyright (C) 1997  Curtis L. Olson  - curt@infoplane.com
   7 //
   8 // This program is free software; you can redistribute it and/or modify
   9 // it under the terms of the GNU General Public License as published by
  10 // the Free Software Foundation; either version 2 of the License, or
  11 // (at your option) any later version.
  12 //
  13 // This program is distributed in the hope that it will be useful,
  14 // but WITHOUT ANY WARRANTY; without even the implied warranty of
  15 // MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
  16 // GNU General Public License for more details.
  17 //
  18 // You should have received a copy of the GNU General Public License
  19 // along with this program; if not, write to the Free Software
  20 // Foundation, Inc., 675 Mass Ave, Cambridge, MA 02139, USA.
  21 //
  22 // $Id$
  23 // (Log is kept at end of this file)
  24 //
  25
  26
  27 #include <dirent.h>
  28 #include <stdio.h>
  29 #include <string.h>
  30 #include <sys/types.h>
  31
  32 #include "../DEM/dem.h"
  33 #include "fixnode.h"
  34 #include "triload.h"
  35
  36
  37 // Storage for the original DEM data which is used to interpolate z values
  38 static fgDEM dem;
  39 static float dem_data[DEM_SIZE_1][DEM_SIZE_1];
  40
  41 // Node list
  42 static double nodes[MAX_NODES][3];
  43
  44
  45 // find all the matching files in the specified directory and fix them
  46 void process_files(char *root_path) {
  47     DIR *d;
  48     struct dirent *de;
  49     char file_path[256];
  50     char *ptr;
  51     int len;
  52
  53     if ( (d = opendir(root_path)) == NULL ) {
  54         printf("cannot open directory '%s'.", root_path);
  55         exit(-1);
  56     }
  57
  58     while ( (de = readdir(d)) != NULL ) {
  59         len = strlen(de->d_name);
  60         if ( len > 7 ) {
  61             ptr = de->d_name;
  62             ptr += (len - 7);
  63             // printf("--> %s \n", ptr);
  64
  65             if ( strcmp(ptr, ".1.node") == 0 ) {
  66                 strcpy(file_path, root_path);
  67                 strcat(file_path, "/");
  68                 strcat(file_path, de->d_name);
  69                 printf("File = %s\n", file_path);
  70
  71                 // load the input data files
  72                 triload(file_path, nodes);
  73
  74                 fixnodes(file_path, dem, dem_data, nodes);
  75             }
  76         }
  77     }
  78 }
  79
  80
  81 // main
  82 int main(int argc, char **argv) {
  83     char demfile[256], root_path[256];
  84
  85     if ( argc != 3 ) {
  86         printf("Usage %s demfile root_path\n", argv[0]);
  87         exit(-1);
  88     }
  89
  90     strcpy(demfile, argv[1]);
  91     strcpy(root_path, argv[2]);
  92
  93     // load the corresponding dem file so we can interpolate elev values
  94     dem.open(demfile);
  95     dem.parse(dem_data);
  96     dem.close();
  97
  98     // process all the *.1.node files in the specified directory
  99     process_files(root_path);
 100
 101     return(0);
 102 }
 103
 104
 105 // $Log$
 106 // Revision 1.6  1998/04/06 21:09:44  curt
 107 // Additional win32 support.
 108 // Fixed a bad bug in dem file parsing that was causing the output to be
 109 // flipped about x = y.
 110 //
 111 // Revision 1.5  1998/03/19 02:50:20  curt
 112 // Updated to support -lDEM class.
 113 //
 114 // Revision 1.4  1998/03/03 16:00:58  curt
 115 // More c++ compile tweaks.
 116 //
 117 // Revision 1.3  1998/01/09 23:03:08  curt
 118 // Restructured to split 1deg x 1deg dem's into 64 subsections.
 119 //
 120 // Revision 1.2  1997/12/02 13:12:07  curt
 121 // Updated to fix every node.
 122 //
 123 // Revision 1.1  1997/11/27 00:17:34  curt
 124 // Initial revision.
 125 //
 126 //