int main(int argc, char **argv){
  std::string analysis = PrintRichdemHeader(argc, argv);

  if(argc!=4){
    std::cout<<"Syntax: "<<argv[0]<<" gentest <SIZE> <OUTPUT NAME>\n";
    std::cout<<"\tGenerates a test DEM\n\n";
    std::cout<<"Syntax: "<<argv[0]<<" <ALGORITHM> <INPUT DEM> <OUTPUT NAME>\n";
    std::cout<<"Algorithms:\n";
    std::cout<<"\t1: Barnes et al. (2014) Flat Resolution Algorithm\n";
    std::cout<<"\t2: Garbtecht and Martz (1997) Flat Resolution Algorithm\n";
    return -1;
  }

  if(argv[1]==std::string("gentest")){
    GenerateDEM(std::stoi(argv[2]),argv[3],analysis);
  } else {
    Router(argv[2],argv[1][0],argv[2],argv[3],analysis);    
  }

  return 0;
}
Exemple #2
0
int main(int argc, char **argv){
  std::string analysis = PrintRichdemHeader(argc,argv);
  
  if(argc!=2){
    std::cerr<<argv[0]<<" <Input>"<<std::endl;
    return -1;
  }

  int width;
  int height;
  GDALDataType dtype;
  std::vector<double> geotransform(6);

  getGDALDimensions(argv[1],height,width,dtype,geotransform.data());

  std::cout<<"Geotransform for '"<<argv[1]<<"': ";
  for(const auto x: geotransform)
    std::cout<<x<<" ";
  std::cout<<width<<" "<<height<<std::endl;

  return 0;
}