1 /* fixnode.c -- traverse the node file and fix the elevation of all the new
4 * Written by Curtis Olson, started November 1997.
6 * Copyright (C) 1997 Curtis L. Olson - curt@infoplane.com
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.
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.
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.
23 * (Log is kept at end of this file)
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] )
44 printf("Fixing up node elevations\n");
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 */
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]); */
56 dem.interpolate_altitude(dem_data, nodes[i][0], nodes[i][1]);
58 /* printf("Fixed: %d %.2f %.2f %.2f\n", i, nodes[i][0],
59 nodes[i][1], nodes[i][2]); */
63 sprintf(toname, "%s.orig", filename);
64 printf("Moving %s to %s\n", filename, toname);
65 rename(filename, toname);
67 printf("Saving new node file: %s\n", filename);
69 fd = fopen(filename, "w");
71 fprintf(fd, "%d 2 1 0\n", nodecount);
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]);
83 /* Revision 1.5 1998/03/19 02:50:19 curt
84 /* Updated to support -lDEM class.
86 * Revision 1.4 1998/03/03 16:00:57 curt
87 * More c++ compile tweaks.
89 * Revision 1.3 1998/01/09 23:03:08 curt
90 * Restructured to split 1deg x 1deg dem's into 64 subsections.
92 * Revision 1.2 1997/12/02 13:12:07 curt
93 * Updated to fix every node.
95 * Revision 1.1 1997/11/27 00:17:33 curt