fly/terragear
1
0
Fork 0
Code Activity

Added support for GSHHS lakes, islands, and ponds.

Browse source
This commit is contained in:
curt 2000-05-06 20:01:58 +00:00
parent 3c558ab7e6
commit 34db12893e
1 changed files with 17 additions and 9 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

26
src/Prep/GSHHS/main.cxx
View file

@ -55,9 +55,9 @@ int main( int argc, char **argv ) {
fglog().setLogLevels( FG_ALL, FG_DEBUG ); fglog().setLogLevels( FG_ALL, FG_DEBUG );
if ( argc < 3 ) { if ( argc < 4 ) {
FG_LOG( FG_GENERAL, FG_ALERT, "Usage: " << argv[0] FG_LOG( FG_GENERAL, FG_ALERT, "Usage: " << argv[0]
<< " <gshhs_file> <work_dir>" ); << " <gshhs_file> <work_dir> <level> [ area_type ]" );
exit(-1); exit(-1);
} }
@ -66,12 +66,20 @@ int main( int argc, char **argv ) {
string command = "mkdir -p " + work_dir; string command = "mkdir -p " + work_dir;
system( command.c_str() ); system( command.c_str() );
// get the specified data level
int target_level = atoi(argv[3]);
if ( target_level < 1 || target_level > 4) {
FG_LOG( FG_GENERAL, FG_ALERT, argv[0] <<
": level must be 1 (land), 2 (lake), 3 (island), or 4 (pond)");
exit(-1);
}
// allow us to override the area type from the command line. All // allow us to override the area type from the command line. All
// polygons in the processed shape file will be assigned this area // polygons in the processed shape file will be assigned this area
// type // type
string force_area_type = ""; string force_area_type = "";
if ( argc == 4 ) { if ( argc == 5 ) {
force_area_type = argv[3]; force_area_type = argv[4];
} }
// initialize persistant polygon counter // initialize persistant polygon counter
@ -157,13 +165,13 @@ int main( int argc, char **argv ) {
} }
max_east = 180000000; /* Only Eurasiafrica needs 270 */ max_east = 180000000; /* Only Eurasiafrica needs 270 */
if ( h.level != 1 ) { if ( h.level != target_level ) {
continue; continue;
} }
if ( h.id != 2 ) { // if ( h.id != 2 ) {
continue; // continue;
} // }
FG_LOG( FG_GENERAL, FG_INFO, " record = " << h.id << " size = " << FG_LOG( FG_GENERAL, FG_INFO, " record = " << h.id << " size = " <<
shape.contour_size( 0 ) ); shape.contour_size( 0 ) );
@ -181,7 +189,7 @@ int main( int argc, char **argv ) {
// cout << "waiting 2 seconds ..." << endl; // cout << "waiting 2 seconds ..." << endl;
// sleep(2); // sleep(2);
if ( shape.contour_size( 0 ) > 10000 ) { if ( shape.contour_size( 0 ) > 16000 ) {
// use my crude line clipping routine to whittle down the // use my crude line clipping routine to whittle down the
// huge polygons into something more managable // huge polygons into something more managable