Use a faster csv parser to parse the point sources

bootstrap.py

@@ -36,14 +36,10 @@
                     "c:/DEV/bld"  # avoid long path issues by using a short build path
                 )
 
-        overlay_ports = os.path.abspath(
-            os.path.join(os.path.dirname(__file__), "deps", "overlay-ports")
-        )
-
         if args.clean:
             vcpkg.clean(triplet=triplet)
         else:
-            vcpkg.bootstrap(ports_dir=os.path.join(".", "deps"), triplet=triplet, build_root=build_root, install_root=install_root, overlay_ports=overlay_ports, clean_after_build=args.clean_after_build)
+            vcpkg.bootstrap(ports_dir=os.path.join(".", "deps"), triplet=triplet, build_root=build_root, install_root=install_root, clean_after_build=args.clean_after_build)
     except KeyboardInterrupt:
         print("\nInterrupted")
         sys.exit(-1)

logic/configurationparser.cpp

@@ -14,6 +14,7 @@ namespace emap {
 
 using namespace inf;
 using namespace std::string_view_literals;
+namespace gdal = inf::gdal;
 
 static EmissionDestination emission_destination_from_string(std::string_view str)
 {

logic/debugtools.cpp

@@ -18,6 +18,7 @@ namespace emap {
 
 using namespace inf;
 using namespace std::string_literals;
+namespace gdal = inf::gdal;
 
 class VectorBuilder
 {

logic/emissioninventory.cpp

@@ -15,6 +15,7 @@ namespace emap {
 
 using namespace inf;
 using namespace date::literals;
+namespace gdal = inf::gdal;
 
 static fs::path throw_if_not_exists(const fs::path& path)
 {

logic/emissionscollector.cpp

@@ -77,9 +77,7 @@ void EmissionsCollector::add_emissions(const CountryCellCoverage& countryInfo, c
     }
 
     for (auto& entry : pointEmissions) {
-        if (entry.value().amount() > 0.0) {
-            _outputBuilder->add_point_output_entry(entry);
-        }
+        _outputBuilder->add_point_output_entry(entry);
     }
 
     if (diffuseEmissions.empty() && !pointEmissions.empty()) {

logic/gridprocessing.cpp

@@ -32,6 +32,7 @@ namespace emap {
 
 using namespace inf;
 using namespace std::string_literals;
+namespace gdal = inf::gdal;
 
 gdal::VectorDataSet transform_vector(const fs::path& vectorPath, const GeoMetadata& destMeta)
 {

logic/include/emap/emissioninventory.h

@@ -1,6 +1,7 @@
 #pragma once
 
 #include "emap/emissions.h"
+#include "infra/math.h"
 
 namespace emap {
 
@@ -223,6 +224,26 @@ class EmissionCollection
         throw inf::RuntimeError("No emission found with id: {}", id);
     }
 
+    const TEmission& emission_with_id_at_coordinate(const EmissionIdentifier& id, Coordinate coord) const
+    {
+        auto iter = std::find_if(_emissions.begin(), _emissions.end(), [&id, &coord](const TEmission& em) {
+            if (em.id() == id) {
+                if (auto coordOpt = em.coordinate(); coordOpt.has_value()) {
+                    return inf::math::approx_equal(coord.x, coordOpt->x, 1e-4) &&
+                           inf::math::approx_equal(coord.y, coordOpt->y, 1e-4);
+                }
+            }
+
+            return false;
+        });
+
+        if (iter == _emissions.end()) {
+            throw inf::RuntimeError("No emission found with id: {} at coordinate {}", id, coord);
+        }
+
+        return *iter;
+    }
+
     std::optional<TEmission> try_emission_with_id(const EmissionIdentifier& id) const noexcept
     {
         auto emissionIter = find_sorted(id);
@@ -243,6 +264,23 @@ class EmissionCollection
         return result;
     }
 
+    std::vector<TEmission> emissions_with_id_at_coordinate(const EmissionIdentifier& id, Coordinate coord) const
+    {
+        std::vector<TEmission> result;
+        std::copy_if(_emissions.begin(), _emissions.end(), std::back_inserter(result), [&id, &coord](const TEmission& em) {
+            if (em.id() == id) {
+                if (auto coordOpt = em.coordinate(); coordOpt.has_value()) {
+                    return inf::math::approx_equal(coord.x, coordOpt->x, 1e-4) &&
+                           inf::math::approx_equal(coord.y, coordOpt->y, 1e-4);
+                }
+            }
+
+            return false;
+        });
+
+        return result;
+    }
+
     size_t empty() const noexcept
     {
         return _emissions.empty();