Allow querying an existing station via interpolate/summarize

earthobservations · Nov 29, 2022 · d598d97 · d598d97
1 parent d54e0f0
commit d598d97
Show file tree

Hide file tree

Showing 2 changed files with 4 additions and 1 deletion.
diff --git a/CHANGELOG.rst b/CHANGELOG.rst
@@ -4,6 +4,9 @@ Changelog
 Development
 ***********
 
+- Interpolation/Summary: Now the queried point can be an existing station laying on the border of the polygon that it's
+  being checked against
+
 0.49.0 (28.11.2022)
 *******************
 

diff --git a/wetterdienst/core/scalar/interpolate.py b/wetterdienst/core/scalar/interpolate.py
@@ -158,7 +158,7 @@ def get_valid_station_groups(stations_dict: dict, utm_x: float, utm_y: float):
     for station_group in combinations(stations_dict.keys(), 4):
         coords = [(stations_dict[s][0], stations_dict[s][1]) for s in station_group]
         pol = Polygon(coords)
-        if pol.contains(point):
+        if pol.covers(point):
             valid_groups.put(station_group)
 
     return valid_groups