3 // dem.h -- DEM management class
5 // Written by Curtis Olson, started March 1998.
7 // Copyright (C) 1998 Curtis L. Olson - curt@me.umn.edu
9 // This program is free software; you can redistribute it and/or
10 // modify it under the terms of the GNU General Public License as
11 // published by the Free Software Foundation; either version 2 of the
12 // License, or (at your option) any later version.
14 // This program is distributed in the hope that it will be useful, but
15 // WITHOUT ANY WARRANTY; without even the implied warranty of
16 // MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
17 // General Public License for more details.
19 // You should have received a copy of the GNU General Public License
20 // along with this program; if not, write to the Free Software
21 // Foundation, Inc., 675 Mass Ave, Cambridge, MA 02139, USA.
24 // (Log is kept at end of this file)
33 #include <Scenery/Bucket/bucketutils.h>
37 #define DEM_SIZE_1 1201
41 // file pointer for input
44 // coordinates (in arc seconds) of south west corner
45 double originx, originy;
47 // number of columns and rows
50 // Distance between column and row data points (in arc seconds)
51 double col_step, row_step;
53 // the actual mesh data allocated here
54 // float dem_data[DEM_SIZE_1][DEM_SIZE_1];
55 // float output_data[DEM_SIZE_1][DEM_SIZE_1];
57 // Current "A" Record Information
58 char dem_description[80], dem_quadrangle[80];
59 double dem_x1, dem_y1, dem_x2, dem_y2, dem_x3, dem_y3, dem_x4, dem_y4;
60 double dem_z1, dem_z2;
61 int dem_resolution, dem_num_profiles;
63 // Current "B" Record Information
64 int prof_col, prof_row;
65 int prof_num_cols, prof_num_rows;
66 double prof_x1, prof_y1;
69 // temporary values for the class to use
76 // Constructor (opens a DEM file)
79 // open a DEM file (use "-" if input is coming from stdin)
80 int open ( char *file );
86 int parse( float dem_data[DEM_SIZE_1][DEM_SIZE_1] );
88 // read and parse DEM "A" record
89 void read_a_record( );
91 // read and parse DEM "B" record
92 void read_b_record( float dem_data[DEM_SIZE_1][DEM_SIZE_1] );
94 // Informational methods
95 double info_originx( void ) { return(originx); }
96 double info_originy( void ) { return(originy); }
98 // return the current altitude based on mesh data. We should
99 // rewrite this to interpolate exact values, but for now this is
101 double interpolate_altitude( float dem_data[DEM_SIZE_1][DEM_SIZE_1],
102 double lon, double lat);
104 // Use least squares to fit a simpler data set to dem data
105 void fit( float dem_data[DEM_SIZE_1][DEM_SIZE_1],
106 float output_data[DEM_SIZE_1][DEM_SIZE_1],
107 char *fg_root, double error, struct fgBUCKET *p );
109 // Initialize output mesh structure
110 void outputmesh_init( float output_data[DEM_SIZE_1][DEM_SIZE_1] );
112 // Get the value of a mesh node
113 double outputmesh_get_pt( float output_data[DEM_SIZE_1][DEM_SIZE_1],
116 // Set the value of a mesh node
117 void outputmesh_set_pt( float output_data[DEM_SIZE_1][DEM_SIZE_1],
118 int i, int j, double value );
120 // Write out a node file that can be used by the "triangle" program
121 void outputmesh_output_nodes( float output_data[DEM_SIZE_1][DEM_SIZE_1],
122 char *fg_root, struct fgBUCKET *p );
133 // Revision 1.1 1998/04/08 22:57:23 curt
134 // Adopted Gnu automake/autoconf system.
136 // Revision 1.2 1998/03/23 20:35:42 curt
137 // Updated to use FG_EPSILON
139 // Revision 1.1 1998/03/19 02:54:47 curt
140 // Reorganized into a class lib called fgDEM.
142 // Revision 1.1 1998/03/19 01:46:29 curt