overlaps.cpp

Go to the documentation of this file.

//=============================================================================
//    Copyright (C) 2021-2024 Wageningen University - All Rights Reserved
//                     Author: Gonzalo Mier
//                        BSD-3 License
//=============================================================================
 
#include "fields2cover/objectives/sg_obj/overlaps.h"
 
namespace f2c::obj {
 
double Overlaps::computeCost(const F2CCell& poly, const F2CSwaths& swaths) {
  return computeCost(F2CCells(poly), swaths);
}
 
double Overlaps::computeCost(const F2CCells& poly, const F2CSwaths& swaths) {
  double area_overlaped {0.0};
  std::vector<F2CCells> areas;
  for (int i = swaths.size() - 1; i >= 0; --i) {
    areas.emplace_back(swaths[i].areaCovered(poly));
    for (int j = i + 1; j < areas.size(); ++j) {
      if (areas[i].intersects(areas[j])) {
        area_overlaped += F2CCells(areas[i]).intersection(areas[j]).area();
      }
    }
  }
  return area_overlaped / poly.area();
}
 
 
}  // namespace f2c::obj