]> git.mxchange.org Git - flightgear.git/blobdiff - Tools/process-dem.pl
Makefile tweaks
[flightgear.git] / Tools / process-dem.pl
index 55d870a0b9125119cacb6f37aa5fa1c4c36f349f..fe7cc3284489bd28185511255b48ec4bfcfba0c3 100755 (executable)
 
 $| = 1;                         # flush buffers after every write
 
+$do_dem2node =     0;
+$do_triangle_1 = 0;
+$do_fixnode =    0;
+$do_splittris =  0;
+$do_assemtris =  0;
+$do_triangle_2 = 0;
 
-# return the file name root (ending at last ".")
-sub file_root {
-    my($file) = @_;
-    my($pos);
-
-    $pos = rindex($file, ".");
-    return substr($file, 0, $pos);
-}
+$do_tri2obj =    0;
+$do_strips =     0;
+$do_fixobj =     1;
 
 
 # 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} = "..";
-    }
+    die "You must remember to set the FG_ROOT environment variable!\n";
 }
 
 
-# 1.  Start with file.dem
+if ( $#ARGV < 1 ) {
+    die "Usage: $0 <error^2> dem-file1 [ dem-file2 dem-file3 ... ]\n";
+}
+
+# Start with file.dem
 
-$dem_file = shift(@ARGV);
 $error = shift(@ARGV);
 $error += 0.0;
 
-print "Source file = $dem_file  Error tolerance = $error\n";
+while ( $dem_file = shift(@ARGV) ) {
+    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";
+    }
+
+
+    if ( $do_dem2node ) {
+       dem2node() ;
+    } else {
+       $subdir = "../Scenery/w120n030/w111n033";
+       print "WARNING:  Hardcoding subdir = $subdir\n";
+    }
+
+    triangle_1() if ( $do_triangle_1 );
+    fixnode() if ( $do_fixnode );
+    splittris() if ( $do_splittris );
+    assemtris() if ( $do_assemtris );
+    triangle_2() if ( $do_triangle_2);
+    tri2obj() if ( $do_tri2obj );
+    strips() if ( $do_strips );
+    fixobj() if ( $do_fixobj );
+}
+
+
+# exit normally
+exit(0);
 
-if ( $error < 0.5 ) {
-    die "I doubt you'll be happy with an error tolerance as low as $error.\n";
+
+# fix command to work with windoze, replaces first "/" with "\\"
+sub fix_command {
+    my($in) = @_;
+
+    $system = `uname -s`;
+    chop($system);
+
+    if ( $system =~ m/CYGWIN32/ ) { 
+       $in =~ s/\//\\\\/;
+    }
+
+    return($in);
+}
+
+
+# return the file name root (ending at last ".")
+sub file_root {
+    my($file) = @_;
+    my($pos);
+
+    $pos = rindex($file, ".");
+    return substr($file, 0, $pos);
 }
 
-# 2.  dem2node $FG_ROOT dem_file tolerance^2 (meters)
+
+# 1.  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 = //;
+sub dem2node {
+    $command = "Dem2node/dem2node $ENV{FG_ROOT} $dem_file $error";
+    $command = fix_command($command);
+    print "Running '$command'\n";
+
+    open(OUT, "$command |");
+    while ( <OUT> ) {
+       print $_;
+       if ( m/^Dir = / ) {
+           $subdir = $_;
+           $subdir =~ s/^Dir = //;
+           chop($subdir);
+       }
     }
-}
-close(OUT);
+    close(OUT);
+} 
 
 
-# 3.  triangle -q file (Takes file.node and produces file.1.node and
+# 2.  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 $_;
+sub triangle_1 {
+    @FILES = `ls $subdir`;
+    foreach $file ( @FILES ) {
+       print $file;
+       chop($file);
+       if ( ($file =~ m/\.node$/) && ($file !~ m/\.\d\.node$/) ) {
+           $command = "Triangle/triangle -q $subdir/$file";
+           $command = fix_command($command);
+           print "Running '$command'\n";
+           open(OUT, "$command |");
+           while ( <OUT> ) {
+               print $_;
+           }
+           close(OUT);
+
+           # remove input file.node
+           unlink("$subdir/$file");
        }
-       close(OUT);
-
-       # remove input file.node
-       unlink("$subdir/$file");
     }
 }
 
 
-# 4.  fixnode file.dem subdir
+# 3.  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 $_;
+sub fixnode {
+    $command = "FixNode/fixnode $dem_file $subdir";
+    $command = fix_command($command);
+    print "Running '$command'\n";
+    open(OUT, "$command |");
+    while ( <OUT> ) {
+       print $_;
+    }
+    close(OUT);
 }
-close(OUT);
 
 
-# 5.  tri2obj file (.1.node)
-#
-#     Take the file.1.node and file.1.ele and produce file.1.obj
+# 4.1 splittris file (.1.node) (.1.ele)
+
+#     Extract the corner, edge, and body vertices (in original
+#     geodetic coordinates) and normals (in cartesian coordinates) and
+#     save them in something very close to the .obj format as file.se,
+#     file.sw, file.nw, file.ne, file.north, file.south, file.east,
+#     file.west, and file.body.  This way we can reconstruct the
+#     region using consistant edges and corners.  
+
+#     Arbitration rules: If an opposite edge file already exists,
+#     don't create our matching edge.  If a corner already exists,
+#     don't create ours.  Basically, the early bird gets the worm and
+#     gets to define the edge verticies and normals.  All the other
+#     adjacent tiles must use these.
+
+sub splittris {
+    @FILES = `ls $subdir`;
+    foreach $file ( @FILES ) {
+       chop($file);
+       if ( $file =~ m/\.1\.node$/ ) {
+           $file =~ s/\.node$//;  # strip off the ".node"
+       
+           $command = "SplitTris/splittris $subdir/$file";
+           $command = fix_command($command);
+           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");
+       }
+    }
+}
+
+
+# 4.2 read in the split of version of the tiles, reconstruct the tile
+#     using the proper shared corners and edges.  Save as a node file
+#     so we can retriangulate.
+
+sub assemtris {
+    @FILES = `ls $subdir`;
+    foreach $file ( @FILES ) {
+       chop($file);
+       if ( $file =~ m/\.1\.body$/ ) {
+           $file =~ s/\.body$//;  # strip off the ".body"
+       
+           $command = "AssemTris/assemtris $subdir/$file";
+           $command = fix_command($command);
+           print "Running '$command'\n";
+           open(OUT, "$command |");
+           while ( <OUT> ) {
+               print $_;
+           }
+           close(OUT);
+       }
+       unlink("$subdir/$file.body");
+    }
+}
+
 
-@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 $_;
+# 4.3 Retriangulate reassembled files (without -q option) so no new
+#     nodes are generated.
+
+sub triangle_2 {
+    @FILES = `ls $subdir`;
+    foreach $file ( @FILES ) {
+       print $file;
+       chop($file);
+       if ( ($file =~ m/\.node$/) && ($file !~ m/\.\d\.node$/) ) {
+           $base = $file;
+           $base =~ s/\.node$//;
+           print("Test for $subdir/$base.q\n");
+           if ( -r "$subdir/$base.q" ) {
+
+               # if triangle hangs, we can create a filebase.q for
+               # the file it hung on.  Then, we test for that file
+               # here which causes the incremental algorithm to run
+               # (which shouldn't ever hang.)
+
+               $command = "Triangle/triangle -i $subdir/$file";
+           } else {
+               $command = "Triangle/triangle $subdir/$file";
+           }
+           $command = fix_command($command);
+           print "Running '$command'\n";
+           open(OUT, "$command |");
+           while ( <OUT> ) {
+               print $_;
+           }
+           close(OUT);
+
+           # remove input file.node
+           unlink("$subdir/$file");
        }
-       close(OUT);
+    }
+}
 
-       unlink("$subdir/$file.node");
-       unlink("$subdir/$file.node.orig");
-       unlink("$subdir/$file.ele");
 
+# 5.  tri2obj file (.1.node) (.1.ele)
+#
+#     Take the file.1.node and file.1.ele and produce file.1.obj
+#
+#     Extracts normals out of the shared edge/vertex files, and uses
+#     the precalcuated normals for these nodes instead of calculating
+#     new ones.  By sharing normals as well as vertices, not only are
+#     the gaps between tiles eliminated, but the colors and lighting
+#     transition smoothly across tile boundaries.
+
+sub tri2obj {
+    @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";
+           $command = fix_command($command);
+           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
+#     Strip the file.1.obj's.  Note, strips doesn't handle the minimal
+#     case of striping a square correctly.
 #
 # 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 $_;
+sub strips {
+    @FILES = `ls $subdir`;
+    foreach $file ( @FILES ) {
+       chop($file);
+       if ( $file =~ m/\.1\.obj$/ ) {
+           $command = "Stripe_u/strips $subdir/$file";
+           $command = fix_command($command);
+           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");
        }
-       close(IN);
-       close(OUT);
-
-       unlink("$subdir/$file");
     }
 }
 
@@ -191,28 +363,89 @@ foreach $file ( @FILES ) {
 #
 #     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 $_;
+sub fixobj {
+    @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";
+           $command = fix_command($command);
+           print "Running '$command'\n";
+           open(OUT, "$command |");
+           while ( <OUT> ) {
+               print $_;
+           }
+           close(OUT);
+
+           # unlink("$subdir/$file");
        }
-       close(OUT);
-
-       unlink("$subdir/$file");
     }
 }
 
 
 #---------------------------------------------------------------------------
 # $Log$
+# Revision 1.18  1998/05/20 20:55:40  curt
+# Makefile tweaks
+#
+# Revision 1.17  1998/04/28 01:23:25  curt
+# Added a work around so that we can get past the "triangle" program
+# hanging, by aborting and rerunning with that tile marked to use the "-i"
+# option.
+#
+# Revision 1.16  1998/04/18 03:57:53  curt
+# Added zlib library support.
+#
+# Revision 1.15  1998/04/08 23:24:07  curt
+# Adopted Gnu automake/autoconf system.
+#
+# Revision 1.14  1998/04/06 21:09:38  curt
+# Additional win32 support.
+# Fixed a bad bug in dem file parsing that was causing the output to be
+# flipped about x = y.
+#
+# Revision 1.13  1998/03/19 02:52:52  curt
+# Updated to reflect some minor tool reorganization and the creation of class
+# to handle DEM processing needs.
+#
+# Revision 1.12  1998/03/19 01:48:35  curt
+# Added gpc-2.01 (generic polygon clipping library)
+#
+# Revision 1.11  1998/03/03 03:36:57  curt
+# Cumulative tweaks.
+#
+# Revision 1.10  1998/02/01 03:42:26  curt
+# Modifications to handle compressed dem files.
+#
+# Revision 1.9  1998/01/27 18:36:54  curt
+# Lots of updates to get back in sync with changes made over in .../Src/
+#
+# Revision 1.8  1998/01/21 17:59:05  curt
+# Uncomment lines to remove several intermediate files.
+#
+# Revision 1.7  1998/01/19 19:51:06  curt
+# A couple final pre-release tweaks.
+#
+# Revision 1.6  1998/01/15 21:33:33  curt
+# Assembling triangles and building a new .node file with the proper shared
+# vertices now works.  Now we just have to use the shared normals and we'll
+# be all set.
+#
+# Revision 1.5  1998/01/15 02:50:08  curt
+# Tweaked to add next stage.
+#
+# Revision 1.4  1998/01/14 15:55:34  curt
+# Finished splittris, started assemtris.
+#
+# Revision 1.3  1998/01/14 02:15:52  curt
+# Updated front end script to keep plugging away on tile fitting.
+#
+# Revision 1.2  1998/01/12 20:42:08  curt
+# Working on fitting tiles together in a seamless manner.
+#
 # Revision 1.1  1998/01/09 23:06:46  curt
 # Initial revision.
 #