flightgear/Tools/process-dem.pl

#!/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);


# 4.1 findedges file (.1.node) (.1.ele)

#     Extract the edge vertices (in original geodetic coordinates) and
#     normals (in cartesian coordinates) and save them in something
#     very close to the .obj format as file.north, file.south,
#     file.east file.west.

@FILES = `ls $subdir`;
foreach $file ( @FILES ) {
    chop($file);
    if ( $file =~ m/\.1\.node$/ ) {
	$file =~ s/\.node$//;  # strip off the ".node"

	$command = "./FindEdges/findedges $subdir/$file";
	print "Running '$command'\n";
	open(OUT, "$command |");
	while ( <OUT> ) {
	    print $_;
	}
	close(OUT);
    }
}


exit(1);


# 4.2 read in tri files (node/ele) skipping edges, read edges out of
#     edge files, save including proper shared edges (as node/ele)
#     files.  If my edge and adjacent edge both exist, use other,
#     delete mine.  If only mine exists, use it.


# 4.3 Retriangulate fixed up files (without -q option)


# 5.  tri2obj file (.1.node) (.1.ele)
#
#     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.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.
#