FixNode/fixnode.cxx

   1 /* fixnode.c -- traverse the node file and fix the elevation of all the new
   2  *              interpolated 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 <stdio.h>
  28 #include <string.h>
  29 #include <unistd.h>
  30
  31 #include "fixnode.hxx"
  32 #include "triload.hxx"
  33
  34
  35 /* load the node information */
  36 void fixnodes( char *filename, fgDEM dem,
  37                float dem_data[DEM_SIZE_1][DEM_SIZE_1],
  38                double nodes[MAX_NODES][3] )
  39 {
  40     char toname[256];
  41     FILE *fd;
  42     int i;
  43
  44     printf("Fixing up node elevations\n");
  45
  46     /* we could just fix the new nodes as the first "for" statement
  47      * would do, but that leads to funny results (I need to figure out
  48      * why.)  So, let's try fixing all of them */
  49
  50     /* for ( i = origcount + 1; i <= nodecount; i++ ) { */
  51     for ( i = 1; i <= nodecount; i++ ) {
  52         /* printf("Current: %d %.2f %.2f %.2f\n", i, nodes[i][0],
  53                nodes[i][1], nodes[i][2]); */
  54
  55         nodes[i][2] =
  56             dem.interpolate_altitude(nodes[i][0], nodes[i][1]);
  57
  58         /* printf("Fixed: %d %.2f %.2f %.2f\n", i, nodes[i][0],
  59                nodes[i][1], nodes[i][2]); */
  60     }
  61
  62
  63     sprintf(toname, "%s.orig", filename);
  64     printf("Moving %s to %s\n", filename, toname);
  65     rename(filename, toname);
  66
  67     printf("Saving new node file: %s\n", filename);
  68
  69     fd = fopen(filename, "w");
  70
  71     fprintf(fd, "%d 2 1 0\n", nodecount);
  72
  73     for ( i = 1; i <= nodecount; i++ ) {
  74         fprintf(fd, "%d %.2f %.2f %.2f 0\n", i, nodes[i][0],
  75                nodes[i][1], nodes[i][2]);
  76     }
  77
  78     fclose(fd);
  79 }
  80
  81
  82 /* $Log$
  83 /* Revision 1.2  1998/04/14 02:26:03  curt
  84 /* Code reorganizations.  Added a Lib/ directory for more general libraries.
  85 /*
  86  * Revision 1.1  1998/04/08 23:05:56  curt
  87  * Adopted Gnu automake/autoconf system.
  88  *
  89  * Revision 1.5  1998/03/19 02:50:19  curt
  90  * Updated to support -lDEM class.
  91  *
  92  * Revision 1.4  1998/03/03 16:00:57  curt
  93  * More c++ compile tweaks.
  94  *
  95  * Revision 1.3  1998/01/09 23:03:08  curt
  96  * Restructured to split 1deg x 1deg dem's into 64 subsections.
  97  *
  98  * Revision 1.2  1997/12/02 13:12:07  curt
  99  * Updated to fix every node.
 100  *
 101  * Revision 1.1  1997/11/27 00:17:33  curt
 102  * Initial revision.
 103  *
 104  */