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