FixNode/main.cxx

   1 // main.cxx -- 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 <sys/types.h>
  28 #include <dirent.h>
  29 #include <stdio.h>
  30 #include <string.h>
  31 #include <string>
  32
  33 #ifdef HAVE_STDLIB_H
  34 #  include <stdlib.h>
  35 #endif // HAVE_STDLIB_H
  36
  37 #include <DEM/dem.hxx>
  38
  39 #include "fixnode.hxx"
  40 #include "triload.hxx"
  41
  42
  43 // Storage for the original DEM data which is used to interpolate z values
  44 fgDEM dem;
  45
  46 // Node list
  47 static double nodes[MAX_NODES][3];
  48
  49
  50 // find all the matching files in the specified directory and fix them
  51 void process_files(const string& root_path) {
  52     DIR *d;
  53     struct dirent *de;
  54     string file_path;
  55     char *ptr;
  56     int len;
  57
  58     if ( (d = opendir( root_path.c_str() )) == NULL ) {
  59         cout << "cannot open directory " + root_path + "\n";
  60         exit(-1);
  61     }
  62
  63     while ( (de = readdir(d)) != NULL ) {
  64         len = strlen(de->d_name);
  65         if ( len > 7 ) {
  66             ptr = de->d_name;
  67             ptr += (len - 7);
  68             // printf("--> %s \n", ptr);
  69
  70             if ( strcmp(ptr, ".1.node") == 0 ) {
  71                 file_path =  root_path + "/" + de->d_name;
  72                 cout << "File = " + file_path + "\n";
  73
  74                 // load the input data files
  75                 triload(file_path.c_str(), nodes);
  76
  77                 fixnodes(file_path.c_str(), &dem, nodes);
  78             }
  79         }
  80     }
  81 }
  82
  83
  84 // main
  85 int main(int argc, char **argv) {
  86     string demfile, root_path;
  87
  88     if ( argc != 3 ) {
  89         printf("Usage %s demfile root_path\n", argv[0]);
  90         exit(-1);
  91     }
  92
  93     cout << "Starting fixnode\n";
  94
  95     demfile = argv[1];
  96     root_path = argv[2];
  97
  98     // load the corresponding dem file so we can interpolate elev values
  99     dem.open(demfile);
 100     dem.parse();
 101     dem.close();
 102
 103     // process all the *.1.node files in the specified directory
 104     process_files(root_path);
 105
 106     return(0);
 107 }
 108
 109
 110 // $Log$
 111 // Revision 1.6  1998/09/19 18:01:27  curt
 112 // Support for changes to libDEM.a
 113 //
 114 // Revision 1.5  1998/07/22 21:46:41  curt
 115 // Fixed a bug that was triggering a seg fault.
 116 //
 117 // Revision 1.4  1998/06/27 16:55:24  curt
 118 // Changed include order for <sys/types.h>
 119 //
 120 // Revision 1.3  1998/04/26 05:02:06  curt
 121 // Added #ifdef HAVE_STDLIB_H
 122 //
 123 // Revision 1.2  1998/04/14 02:26:04  curt
 124 // Code reorganizations.  Added a Lib/ directory for more general libraries.
 125 //
 126 // Revision 1.1  1998/04/08 23:05:57  curt
 127 // Adopted Gnu automake/autoconf system.
 128 //
 129 // Revision 1.6  1998/04/06 21:09:44  curt
 130 // Additional win32 support.
 131 // Fixed a bad bug in dem file parsing that was causing the output to be
 132 // flipped about x = y.
 133 //
 134 // Revision 1.5  1998/03/19 02:50:20  curt
 135 // Updated to support -lDEM class.
 136 //
 137 // Revision 1.4  1998/03/03 16:00:58  curt
 138 // More c++ compile tweaks.
 139 //
 140 // Revision 1.3  1998/01/09 23:03:08  curt
 141 // Restructured to split 1deg x 1deg dem's into 64 subsections.
 142 //
 143 // Revision 1.2  1997/12/02 13:12:07  curt
 144 // Updated to fix every node.
 145 //
 146 // Revision 1.1  1997/11/27 00:17:34  curt
 147 // Initial revision.
 148 //
 149 //