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)
32 # error This library requires C++
38 #include <Bucket/bucketutils.h>
40 //#include <zlib/zlib.h>
41 #include <Misc/fgstream.hxx>
45 #define DEM_SIZE_1 1201
49 // file pointer for input
53 // coordinates (in arc seconds) of south west corner
54 double originx, originy;
56 // number of columns and rows
59 // Distance between column and row data points (in arc seconds)
60 double col_step, row_step;
62 // pointers to the actual mesh data allocated here
63 float (*dem_data)[DEM_SIZE_1];
64 float (*output_data)[DEM_SIZE_1];
66 // Current "A" Record Information
67 char dem_description[80], dem_quadrangle[80];
68 double dem_x1, dem_y1, dem_x2, dem_y2, dem_x3, dem_y3, dem_x4, dem_y4;
69 double dem_z1, dem_z2;
70 int dem_resolution, dem_num_profiles;
72 // Current "B" Record Information
73 int prof_col, prof_row;
74 int prof_num_cols, prof_num_rows;
75 double prof_x1, prof_y1;
78 // temporary values for the class to use
83 // return next token from input stream
86 // return next integer from input stream
89 // return next double from input stream
92 // return next exponential num from input stream
100 // open a DEM file (use "-" if input is coming from stdin)
101 int open ( const string& file );
109 // read and parse DEM "A" record
112 // read and parse DEM "B" record
113 void read_b_record();
115 // Informational methods
116 double info_originx( void ) { return(originx); }
117 double info_originy( void ) { return(originy); }
119 // return the current altitude based on mesh data. We should
120 // rewrite this to interpolate exact values, but for now this is
122 double interpolate_altitude( double lon, double lat );
124 // Use least squares to fit a simpler data set to dem data
125 void fit( double error, fgBUCKET *p );
127 // Initialize output mesh structure
128 void outputmesh_init( void );
130 // Get the value of a mesh node
131 double outputmesh_get_pt( int i, int j );
133 // Set the value of a mesh node
134 void outputmesh_set_pt( int i, int j, double value );
136 // Write out a node file that can be used by the "triangle" program
137 void outputmesh_output_nodes( const string& fg_root, fgBUCKET *p );
148 // Revision 1.8 1998/09/19 17:59:46 curt
149 // Use c++ streams (fg_gzifstream). Also converted many character arrays to
152 // Revision 1.7 1998/07/04 00:47:19 curt
153 // typedef'd struct fgBUCKET.
155 // Revision 1.6 1998/06/05 18:14:40 curt
156 // Abort out early when reading the "A" record if it doesn't look like
157 // a proper DEM file.
159 // Revision 1.5 1998/04/22 13:14:46 curt
160 // Fixed a bug in zlib usage.
162 // Revision 1.4 1998/04/21 17:03:41 curt
163 // Prepairing for C++ integration.
165 // Revision 1.3 1998/04/18 03:53:06 curt
166 // Added zlib support.
168 // Revision 1.2 1998/04/14 02:43:28 curt
169 // Used "new" to auto-allocate large DEM parsing arrays in class constructor.
171 // Revision 1.1 1998/04/08 22:57:23 curt
172 // Adopted Gnu automake/autoconf system.
174 // Revision 1.2 1998/03/23 20:35:42 curt
175 // Updated to use FG_EPSILON
177 // Revision 1.1 1998/03/19 02:54:47 curt
178 // Reorganized into a class lib called fgDEM.
180 // Revision 1.1 1998/03/19 01:46:29 curt