DemChop/demchop.cxx

   1 // demchop.cxx -- chop up a dem file into it's corresponding pieces and stuff
   2 //                them into the workspace directory
   3 //
   4 // Written by Curtis Olson, started March 1999.
   5 //
   6 // Copyright (C) 1997  Curtis L. Olson  - curt@flightgear.org
   7 //
   8 // This program is free software; you can redistribute it and/or modify
   9 // it under the terms of the GNU General Public License as published by
  10 // the Free Software Foundation; either version 2 of the License, or
  11 // (at your option) any later version.
  12 //
  13 // This program is distributed in the hope that it will be useful,
  14 // but WITHOUT ANY WARRANTY; without even the implied warranty of
  15 // MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
  16 // GNU General Public License for more details.
  17 //
  18 // You should have received a copy of the GNU General Public License
  19 // along with this program; if not, write to the Free Software
  20 // Foundation, Inc., 675 Mass Ave, Cambridge, MA 02139, USA.
  21 //
  22 // $Id$
  23 // (Log is kept at end of this file)
  24
  25
  26 #include <Include/compiler.h>
  27
  28 #include STL_STRING
  29
  30 #include <Debug/logstream.hxx>
  31 #include <Bucket/newbucket.hxx>
  32 #include <DEM/dem.hxx>
  33
  34 #include "point2d.hxx"
  35
  36 FG_USING_STD(string);
  37
  38
  39 int main(int argc, char **argv) {
  40     /*
  41     fgDEM dem;
  42     FGBucket p;
  43     string fg_root;
  44     string filename;
  45     double error;
  46     int i, j;
  47     */
  48
  49     fglog().setLogLevels( FG_ALL, FG_DEBUG );
  50
  51     if ( argc != 4 ) {
  52         FG_LOG( FG_GENERAL, FG_ALERT,
  53                 "Usage " << argv[0] << " <dem_file> <work_dir>" );
  54         exit(-1);
  55     }
  56
  57     string dem_name = argv[1];
  58     string work_dir = argv[2];
  59     string command = "mkdir -p " + work_dir;
  60     system( command.c_str() );
  61
  62     fgDEM dem(dem_name);
  63     dem.parse();
  64     dem.close();
  65
  66     point2d min, max;
  67     min.x = dem.get_originx() / 3600.0;
  68     min.y = dem.get_originy() / 3600.0;
  69     FGBucket b_min( min.x, min.y );
  70
  71     max.x = min.x + ( dem.get_cols() * dem.get_col_step() ) / 3600.0;
  72     max.y = min.y + ( dem.get_rows() * dem.get_row_step() ) / 3600.0;
  73     FGBucket b_max( max.x, max.y );
  74
  75     if ( b_min == b_max ) {
  76         dem.write_area( b_min );
  77     } else {
  78         FGBucket b_cur;
  79         int dx, dy, i, j;
  80
  81         fgBucketDiff(b_min, b_max, &dx, &dy);
  82         cout << "DEM file spans tile boundaries" << endl;
  83         cout << "  dx = " << dx << "  dy = " << dy << endl;
  84
  85         if ( (dx > 20) || (dy > 20) ) {
  86             cout << "somethings really wrong!!!!" << endl;
  87             exit(-1);
  88         }
  89
  90         for ( j = 0; j <= dy; j++ ) {
  91             for ( i = 0; i <= dx; i++ ) {
  92                 b_cur = fgBucketOffset(min.x, min.y, i, j);
  93                 dem.write_area( b_cur );
  94             }
  95         }
  96     }
  97
  98     return 0;
  99 }
 100
 101
 102 // $Log$
 103 // Revision 1.2  1999/03/10 16:09:44  curt
 104 // Hacking towards the first working version.
 105 //
 106 // Revision 1.1  1999/03/10 01:02:54  curt
 107 // Initial revision.
 108 //