Function main

tests/test_reconstruct_api.cpp:48–171

void print_version() { std::cout << "roofer test_reconstruct_api\n"; }

int main(int argc, const char* argv[]) {
  auto cmdl = argh::parser();

  cmdl.parse(argc, argv);
  std::string program_name = cmdl[0];

  bool verbose = cmdl[{"-v", "--verbose"}];
  bool version = cmdl[{"-V", "--version"}];

  if (cmdl[{"-h", "--help"}]) {
    print_help(program_name);
    return EXIT_SUCCESS;
  }
  if (version) {
    print_version();
    return EXIT_SUCCESS;
  }

  auto& logger = roofer::logger::Logger::get_logger();

  if (verbose) {
    logger.set_level(roofer::logger::LogLevel::debug);
  } else {
    logger.set_level(roofer::logger::LogLevel::warning);
  }
  std::string path_pointcloud =
      "data/wippolder/objects/503100000030812/crop/"
      "503100000030812_pointcloud.las";
  std::string path_footprint =
      "data/wippolder/objects/503100000030812/crop/503100000030812.gpkg";
  float floor_elevation = -0.16899998486042023;

  // read inputs
  auto pj = roofer::misc::createProjHelper();
  auto PointReader = roofer::io::createPointCloudReaderLASlib(*pj);
  auto VectorReader = roofer::io::createVectorReaderOGR(*pj);

  VectorReader->open(path_footprint);
  std::vector<roofer::LinearRing> footprints;
  VectorReader->readPolygons(footprints);

  PointReader->open(path_pointcloud);
  logger.info("Reading pointcloud from {}", path_pointcloud);
  roofer::vec1i classification;
  roofer::PointCollection points, points_ground, points_roof;
  roofer::AttributeVecMap attributes;
  PointReader->readPointCloud(points, &classification);
  logger.info("Read {} points", points.size());

  // split into ground and roof points
  for (size_t i = 0; i < points.size(); ++i) {
    if (2 == classification[i]) {
      points_ground.push_back(points[i]);
    } else if (6 == classification[i]) {
      points_roof.push_back(points[i]);
    }
  }
  logger.info("{} ground points and {} roof points", points_ground.size(),