Update WMOreg: case island + coastal stations

Stations/WMOref.py

@@ -36,8 +36,8 @@ class WMOREG:
         """
         # Point is (x,y)-(lon,lat) 
         p1 = Point(lon, lat)
-        out = np.array([p1.within(poly) for poly in self.polygonWMO])
-        try : 
+        out = np.array([poly.contains(p1) for poly in self.polygonWMO])
+        try :
             index = np.where(out)[0][0]
             stReg, stSubreg = self.attributesPolygonWMO[index]
             if output=="num":
@@ -46,5 +46,19 @@ class WMOREG:
                 return [stReg, stSubreg]
         except:
             print("lon: {0}; lat: {1}".format(lon,lat))
-            print("No valid WMO region has been found, verify the lon lat")
+            print("Inexact location : finding the closest subregion.")
+            out = np.array([self.get_distance(poly,p1) for poly in self.polygonWMO])
+            index = np.argmin(out)
+            stReg, stSubreg = self.attributesPolygonWMO[index]
+            if output=="num":
+                return stReg*100+stSubreg
+            elif output == "list":
+                return [stReg, stSubreg]
+            
+    def get_distance(self, poly,p1):
+        try:
+           d = poly.exterior.distance(p1)
+        except:
+           d = np.min([p.exterior.distance(p1) for p in poly])
+        return d