Initial revision.

diff --git a/Tools/process-dem.pl b/Tools/process-dem.pl
new file mode 100755 (executable)
index 0000000..55d870a
--- /dev/null
+++ b/Tools/process-dem.pl
@@ -0,0 +1,218 @@
+#!/usr/bin/perl
+
+#---------------------------------------------------------------------------
+# Toplevel script to automate DEM file processing and conversion
+#
+# Written by Curtis Olson, started January 1998.
+#
+# Copyright (C) 1997  Curtis L. Olson  - curt@infoplane.com
+#
+# This program is free software; you can redistribute it and/or modify
+# it under the terms of the GNU General Public License as published by
+# the Free Software Foundation; either version 2 of the License, or
+# (at your option) any later version.
+#
+# This program is distributed in the hope that it will be useful,
+# but WITHOUT ANY WARRANTY; without even the implied warranty of
+# MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
+# GNU General Public License for more details.
+#
+# You should have received a copy of the GNU General Public License
+# along with this program; if not, write to the Free Software
+# Foundation, Inc., 675 Mass Ave, Cambridge, MA 02139, USA.
+#
+# $Id$
+# (Log is kept at end of this file)
+#---------------------------------------------------------------------------
+
+
+$| = 1;                         # flush buffers after every write
+
+
+# return the file name root (ending at last ".")
+sub file_root {
+    my($file) = @_;
+    my($pos);
+
+    $pos = rindex($file, ".");
+    return substr($file, 0, $pos);
+}
+
+
+# set the FG_ROOT environment variable if it hasn't already been set.
+if ( $ENV{FG_ROOT} eq "" ) {
+    # look for a file called fgtop as a place marker
+    if ( -e "fgtop" ) {
+        $ENV{FG_ROOT} = ".";
+    } elsif ( -e "../fgtop" ) {
+        $ENV{FG_ROOT} = "..";
+    }
+}
+
+
+# 1.  Start with file.dem
+
+$dem_file = shift(@ARGV);
+$error = shift(@ARGV);
+$error += 0.0;
+
+print "Source file = $dem_file  Error tolerance = $error\n";
+
+if ( $error < 0.5 ) {
+    die "I doubt you'll be happy with an error tolerance as low as $error.\n";
+}
+
+# 2.  dem2node $FG_ROOT dem_file tolerance^2 (meters)
+# 
+#     - dem2node .. dem_file 160000
+#
+#     splits dem file into 64 file.node's which contain the
+#     irregularly fitted vertices
+
+$command = "./Dem2node/demfit $ENV{FG_ROOT} $dem_file $error";
+
+print "Running '$command'\n";
+
+open(OUT, "$command |");
+while ( <OUT> ) {
+    print $_;
+    if ( m/Scenery/ ) {
+       $subdir = $_;
+       $subdir =~ s/Dir = //;
+    }
+}
+close(OUT);
+
+
+# 3.  triangle -q file (Takes file.node and produces file.1.node and
+#                      file.1.ele)
+
+$subdir = "../Scenery/w120n030/w111n033";
+print "Subdirectory for this dem file is $subdir\n";
+
+@FILES = `ls $subdir`;
+foreach $file ( @FILES ) {
+    print $file;
+    chop($file);
+    if ( ($file =~ m/\.node$/) && ($file !~ m/\.\d\.node$/) ) {
+       $command = "./Triangle/triangle -q $subdir/$file";
+       print "Running '$command'\n";
+       open(OUT, "$command |");
+       while ( <OUT> ) {
+           print $_;
+       }
+       close(OUT);
+
+       # remove input file.node
+       unlink("$subdir/$file");
+    }
+}
+
+
+# 4.  fixnode file.dem subdir
+#
+#     Take the original .dem file (for interpolating Z values) and the
+#     subdirecotry containing all the file.1.node's and replace with
+#     fixed file.1.node
+
+$command = "./FixNode/fixnode $dem_file $subdir";
+print "Running '$command'\n";
+open(OUT, "$command |");
+while ( <OUT> ) {
+    print $_;
+}
+close(OUT);
+
+
+# 5.  tri2obj file (.1.node)
+#
+#     Take the file.1.node and file.1.ele and produce file.1.obj
+
+@FILES = `ls $subdir`;
+foreach $file ( @FILES ) {
+    chop($file);
+    if ( $file =~ m/\.1\.node$/ ) {
+       $file =~ s/\.node$//;  # strip off the ".node"
+
+       $command = "./Tri2obj/tri2obj $subdir/$file";
+       print "Running '$command'\n";
+       open(OUT, "$command |");
+       while ( <OUT> ) {
+           print $_;
+       }
+       close(OUT);
+
+       unlink("$subdir/$file.node");
+       unlink("$subdir/$file.node.orig");
+       unlink("$subdir/$file.ele");
+
+    }
+}
+
+
+# 6.  strip file.1.obj
+# 
+#     Strip the file.1.obj's
+#
+# 7.  cp bands.d file.2.obj
+#
+#     strips produces a file called "bands.d" ... copy this to file.2.obj
+
+@FILES = `ls $subdir`;
+foreach $file ( @FILES ) {
+    chop($file);
+    if ( $file =~ m/\.1\.obj$/ ) {
+       $command = "./Stripe_u/strips $subdir/$file";
+       print "Running '$command'\n";
+       open(OUT, "$command |");
+       while ( <OUT> ) {
+           print $_;
+       }
+       close(OUT);
+
+       # copy to destination file
+       $newfile = $file;
+       $newfile =~ s/\.1\.obj$//;
+       print "Copying to $subdir/$newfile.2.obj\n";
+       open(IN, "<bands.d");
+       open(OUT, ">$subdir/$newfile.2.obj");
+       while ( <IN> ) {
+           print OUT $_;
+       }
+       close(IN);
+       close(OUT);
+
+       unlink("$subdir/$file");
+    }
+}
+
+
+# 8.  fixobj file-new
+#
+#     Sort file.2.obj by strip winding
+
+@FILES = `ls $subdir`;
+foreach $file ( @FILES ) {
+    chop($file);
+    if ( $file =~ m/\.2\.obj$/ ) {
+       $newfile = $file;
+       $newfile =~ s/\.2\.obj$/.obj/;
+
+       $command = "./FixObj/fixobj $subdir/$file $subdir/$newfile";
+       print "Running '$command'\n";
+       open(OUT, "$command |");
+       while ( <OUT> ) {
+           print $_;
+       }
+       close(OUT);
+
+       unlink("$subdir/$file");
+    }
+}
+
+
+#---------------------------------------------------------------------------
+# $Log$
+# Revision 1.1  1998/01/09 23:06:46  curt
+# Initial revision.
+#