Added a small tool to rectify heights using cliff contours placed in
the height array directory.
This commit is contained in:
James.Hester
parent
990417c000
commit
09d61732cd
1 changed files with 131 additions and 0 deletions
131
src/Lib/Array/rectify_height.cxx
Normal file
131
src/Lib/Array/rectify_height.cxx
Normal file
|
|
@ -0,0 +1,131 @@
|
||||||
|
//rectify_height.cxx
|
||||||
|
#ifdef HAVE_CONFIG_H
|
||||||
|
# include <config.h>
|
||||||
|
#endif
|
||||||
|
|
||||||
|
#include <simgear/compiler.h>
|
||||||
|
#include <simgear/bucket/newbucket.hxx>
|
||||||
|
#include <simgear/misc/sg_path.hxx>
|
||||||
|
#include <simgear/debug/logstream.hxx>
|
||||||
|
|
||||||
|
#include <sys/types.h>
|
||||||
|
#include <iostream>
|
||||||
|
#include <string.h>
|
||||||
|
#include <stdlib.h>
|
||||||
|
#include "array.hxx"
|
||||||
|
|
||||||
|
/* This program will find all height files in the directory provided,
|
||||||
|
and if they have an associated cliff file, will adjust and then output
|
||||||
|
the heights. Single threaded to avoid reading/writing the same file. */
|
||||||
|
|
||||||
|
// display usage and exit
|
||||||
|
static void usage( const std::string name ) {
|
||||||
|
SG_LOG(SG_GENERAL, SG_ALERT, "Usage: " << name);
|
||||||
|
SG_LOG(SG_GENERAL, SG_ALERT, " --work-dir=<directory>");
|
||||||
|
SG_LOG(SG_GENERAL, SG_ALERT, " --height-dir=<directory>");
|
||||||
|
SG_LOG(SG_GENERAL, SG_ALERT, "[ --tile-id=<id>]");
|
||||||
|
SG_LOG(SG_GENERAL, SG_ALERT, " --min-lon=<degrees>");
|
||||||
|
SG_LOG(SG_GENERAL, SG_ALERT, " --max-lon=<degrees>");
|
||||||
|
SG_LOG(SG_GENERAL, SG_ALERT, " --min-lat=<degrees>");
|
||||||
|
SG_LOG(SG_GENERAL, SG_ALERT, " --max-lat=<degrees>");
|
||||||
|
SG_LOG(SG_GENERAL, SG_ALERT, "[ --min-dist=<float>]");
|
||||||
|
exit(-1);
|
||||||
|
}
|
||||||
|
|
||||||
|
int main( int argc, char **argv) {
|
||||||
|
std::string output_dir = ".";
|
||||||
|
std::string work_dir = ".";
|
||||||
|
std::string height_dir = "";
|
||||||
|
SGGeod min,max;
|
||||||
|
long tile_id = -1;
|
||||||
|
double bad_zone = 30; //distance in m from cliff to rectify
|
||||||
|
|
||||||
|
sglog().setLogLevels(SG_ALL,SG_INFO);
|
||||||
|
sglog().set_log_priority(SG_DEBUG);
|
||||||
|
|
||||||
|
//
|
||||||
|
// Parse the command-line arguments.
|
||||||
|
//
|
||||||
|
int arg_pos;
|
||||||
|
for (arg_pos = 1; arg_pos < argc; arg_pos++) {
|
||||||
|
std::string arg = argv[arg_pos];
|
||||||
|
|
||||||
|
if (arg.find("--work-dir=") == 0) {
|
||||||
|
work_dir = arg.substr(11);
|
||||||
|
} else if (arg.find("--height-dir=") == 0) {
|
||||||
|
height_dir = arg.substr(13);
|
||||||
|
} else if (arg.find("--tile-id=") == 0) {
|
||||||
|
tile_id = atol(arg.substr(10).c_str());
|
||||||
|
} else if ( arg.find("--min-lon=") == 0 ) {
|
||||||
|
min.setLongitudeDeg(atof( arg.substr(10).c_str() ));
|
||||||
|
} else if ( arg.find("--max-lon=") == 0 ) {
|
||||||
|
max.setLongitudeDeg(atof( arg.substr(10).c_str() ));
|
||||||
|
} else if ( arg.find("--min-lat=") == 0 ) {
|
||||||
|
min.setLatitudeDeg(atof( arg.substr(10).c_str() ));
|
||||||
|
} else if ( arg.find("--max-lat=") == 0 ) {
|
||||||
|
max.setLatitudeDeg(atof( arg.substr(10).c_str() ));
|
||||||
|
} else if ( arg.find("--min-dist=") == 0) {
|
||||||
|
bad_zone = atof(arg.substr(11).c_str());
|
||||||
|
} else if (arg.find("--") == 0) {
|
||||||
|
usage(argv[0]);
|
||||||
|
} else {
|
||||||
|
break;
|
||||||
|
}
|
||||||
|
}
|
||||||
|
|
||||||
|
SG_LOG(SG_GENERAL, SG_ALERT, "Working directory is " << work_dir);
|
||||||
|
SG_LOG(SG_GENERAL, SG_ALERT, "Heights are in " << height_dir);
|
||||||
|
SG_LOG(SG_GENERAL, SG_ALERT, "Rectification zone within " << bad_zone << "m of cliffs");
|
||||||
|
|
||||||
|
if ( tile_id > 0 ) {
|
||||||
|
SG_LOG(SG_GENERAL, SG_ALERT, "Tile id is " << tile_id);
|
||||||
|
} else {
|
||||||
|
if (min.isValid() && max.isValid() && (min != max))
|
||||||
|
{
|
||||||
|
SG_LOG(SG_GENERAL, SG_ALERT, "Longitude = " << min.getLongitudeDeg() << ':' << max.getLongitudeDeg());
|
||||||
|
SG_LOG(SG_GENERAL, SG_ALERT, "Latitude = " << min.getLatitudeDeg() << ':' << max.getLatitudeDeg());
|
||||||
|
} else
|
||||||
|
{
|
||||||
|
SG_LOG(SG_GENERAL, SG_ALERT, "Lon/Lat unset or wrong");
|
||||||
|
exit(1);
|
||||||
|
}
|
||||||
|
}
|
||||||
|
|
||||||
|
// Now generate the buckets
|
||||||
|
std::vector<SGBucket> bucketList;
|
||||||
|
if (tile_id == -1) {
|
||||||
|
// build all the tiles in an area
|
||||||
|
SG_LOG(SG_GENERAL, SG_ALERT, "Fixing heights within given bounding box");
|
||||||
|
|
||||||
|
SGBucket b_min( min );
|
||||||
|
SGBucket b_max( max );
|
||||||
|
|
||||||
|
if ( b_min == b_max ) {
|
||||||
|
bucketList.push_back( b_min );
|
||||||
|
} else {
|
||||||
|
SG_LOG(SG_GENERAL, SG_ALERT, " adjustment area spans tile boundaries");
|
||||||
|
sgGetBuckets( min, max, bucketList );
|
||||||
|
}
|
||||||
|
} else {
|
||||||
|
// adjust the specified tile
|
||||||
|
SG_LOG(SG_GENERAL, SG_ALERT, "Adjusting tile " << tile_id);
|
||||||
|
bucketList.push_back( SGBucket( tile_id ) );
|
||||||
|
}
|
||||||
|
|
||||||
|
// Finally, loop over the buckets adjusting heights
|
||||||
|
for (unsigned int i=0; i < bucketList.size(); i++) {
|
||||||
|
SGBucket bucket = bucketList[i];
|
||||||
|
TGArray array;
|
||||||
|
std::string base_path = bucket.gen_base_path();
|
||||||
|
std::string array_path = work_dir + "/" + height_dir + "/" + base_path + "/" + bucket.gen_index_str();
|
||||||
|
if (!array.open(array_path)) {
|
||||||
|
SG_LOG(SG_GENERAL,SG_DEBUG, "Failed to open array file " << array_path);
|
||||||
|
continue;
|
||||||
|
}
|
||||||
|
array.parse(bucket);
|
||||||
|
array.close();
|
||||||
|
array.remove_voids();
|
||||||
|
array.rectify_heights(bad_zone);
|
||||||
|
array.write_bin(work_dir + "/" + height_dir,true,bucket);
|
||||||
|
}
|
||||||
|
}
|
||||||
Loading…
Reference in a new issue