// (Log is kept at end of this file)
+#ifdef HAVE_CONFIG_H
+# include <config.h>
+#endif
+
#include <ctype.h> // isspace()
+#include <stdlib.h> // atoi()
#include <math.h> // rint()
#include <stdio.h>
-#include <stdlib.h> // atoi()
#include <string.h>
#include <sys/stat.h> // stat()
#include <unistd.h> // stat()
#include <Include/fg_constants.h>
+#define MAX_EX_NODES 10000
+
#ifdef WIN32
# define MKDIR(a) mkdir(a,S_IRWXU) // I am just guessing at this flag (NHV)
#endif // WIN32
// read and parse DEM "A" record
-void fgDEM::read_a_record( void ) {
+int fgDEM::read_a_record( void ) {
int i, inum;
double dnum;
char name[144];
inum = next_int(fd);
printf(" DEM level code = %d\n", inum);
+ if ( inum > 3 ) {
+ return(0);
+ }
+
// Pattern code, 1 indicates a regular elevation pattern
inum = next_int(fd);
printf(" Pattern code = %d\n", inum);
// number of profiles
dem_num_profiles = cols = next_int(fd);
printf(" Expecting %d profiles\n", dem_num_profiles);
+
+ return(1);
}
int fgDEM::parse( void ) {
int i;
- cur_row = 0;
+ cur_col = 0;
- read_a_record();
+ if ( !read_a_record() ) {
+ return(0);
+ }
for ( i = 0; i < dem_num_profiles; i++ ) {
// printf("Ready to read next b record\n");
printf(" Done parsing\n");
- return(0);
+ return(1);
}
// Use least squares to fit a simpler data set to dem data
-void fgDEM::fit( char *fg_root, double error, struct fgBUCKET *p ) {
+void fgDEM::fit( double error, fgBUCKET *p ) {
double x[DEM_SIZE_1], y[DEM_SIZE_1];
double m, b, ave_error, max_error;
double cury, lasty;
// printf("Please hit return: "); gets(junk);
}
- outputmesh_output_nodes(fg_root, p);
+ // outputmesh_output_nodes(fg_root, p);
}
}
-// Write out a node file that can be used by the "triangle" program
-void fgDEM::outputmesh_output_nodes( char *fg_root, struct fgBUCKET *p ) {
+// Write out a node file that can be used by the "triangle" program.
+// Check for an optional "index.node.ex" file in case there is a .poly
+// file to go along with this node file. Include these nodes first
+// since they are referenced by position from the .poly file.
+void fgDEM::outputmesh_output_nodes( char *fg_root, fgBUCKET *p ) {
+ double exnodes[MAX_EX_NODES][3];
+ double junk1, junk2, junk3;
struct stat stat_buf;
- char base_path[256], dir[256], file[256];
+ char base_path[256], dir[256], file[256], exfile[256];
#ifdef WIN32
char tmp_path[256];
#endif
FILE *fd;
long int index;
int colmin, colmax, rowmin, rowmax;
- int i, j, count, result;
+ int i, j, count, excount, result;
// determine dimensions
colmin = p->x * ( (cols - 1) / 8);
index = fgBucketGenIndex(p);
sprintf(file, "%s/%ld.node", dir, index);
+ // get (optional) extra node file name (in case there is matching
+ // .poly file.
+ strcpy(exfile, file);
+ strcat(exfile, ".ex");
+
+ // load extra nodes if they exist
+ excount = 0;
+ if ( (fd = fopen(exfile, "r")) != NULL ) {
+ fscanf(fd, "%d %d %d %d", &excount, &junk1, &junk2, &junk3);
+
+ if ( excount > MAX_EX_NODES - 1 ) {
+ printf("Error, too many 'extra' nodes, increase array size\n");
+ exit(-1);
+ } else {
+ printf(" Expecting %d 'extra' nodes\n", excount);
+ }
+
+ for ( i = 1; i <= excount; i++ ) {
+ fscanf(fd, "%d %lf %lf %lf\n", &junk1,
+ &exnodes[i][0], &exnodes[i][1], &exnodes[i][2]);
+ printf("(extra) %d %.2f %.2f %.2f\n",
+ i, exnodes[i][0], exnodes[i][1], exnodes[i][2]);
+ }
+ }
+ fclose(fd);
+
printf("Creating node file: %s\n", file);
fd = fopen(file, "w");
- // first count nodes to generate header
+ // first count regular nodes to generate header
count = 0;
for ( j = rowmin; j <= rowmax; j++ ) {
for ( i = colmin; i <= colmax; i++ ) {
}
// printf(" count = %d\n", count);
}
- fprintf(fd, "%d 2 1 0\n", count);
+ fprintf(fd, "%d 2 1 0\n", count + excount);
+
+ // now write out extra node data
+ for ( i = 1; i <= excount; i++ ) {
+ fprintf(fd, "%d %.2f %.2f %.2f\n",
+ i, exnodes[i][0], exnodes[i][1], exnodes[i][2]);
+ }
- // now write out actual node data
- count = 1;
+ // write out actual node data
+ count = excount + 1;
for ( j = rowmin; j <= rowmax; j++ ) {
for ( i = colmin; i <= colmax; i++ ) {
if ( output_data[i][j] > -9000.0 ) {
fgDEM::~fgDEM( void ) {
// printf("class fgDEM DEstructor called.\n");
+ free(dem_data);
+ free(output_data);
}
// $Log$
+// Revision 1.11 1998/07/20 12:46:11 curt
+// When outputing to a .node file, first check for an optional
+// "index.node.ex" file in case there is a .poly file to go along with this
+// node file. Include these nodes first since they are referenced by position
+// from the .poly file. This is my first pass at adding an area "cutout"
+// feature to the terrain generation pipeline.
+//
+// Revision 1.10 1998/07/13 20:58:02 curt
+// .
+//
+// Revision 1.9 1998/07/13 15:29:49 curt
+// Added #ifdef HAVE_CONFIG_H
+//
+// Revision 1.8 1998/07/04 00:47:18 curt
+// typedef'd struct fgBUCKET.
+//
+// Revision 1.7 1998/06/05 18:14:39 curt
+// Abort out early when reading the "A" record if it doesn't look like
+// a proper DEM file.
+//
+// Revision 1.6 1998/05/02 01:49:21 curt
+// Fixed a bug where the wrong variable was being initialized.
+//
// Revision 1.5 1998/04/25 15:00:32 curt
// Changed "r" to "rb" in gzopen() options. This fixes bad behavior in win32.
//