fly/flightgear
1
0
Fork 0
Code Activity

Initial revision.

Browse source
This commit is contained in:
curt 1998-01-09 23:06:46 +00:00
parent 0e813d878b
commit 35b9b8a360
1 changed files with 218 additions and 0 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

218
Tools/process-dem.pl Executable file
View file

@ -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.
#