The Swiss confederation provides a shapefile with the boundaries of cantons and municipals at:
The Swiss government has used the Swiss coordinate system (LV03 and LV95).. 🤦♂️
There are multiple rows defining different POLYGON for each canton, my guess is that they didn't know that MULTIPOLYGON exists and so instead of using this class they repeat the row for each POLYGON...
After having cast the coordinates to the classic coordinate system, We can simply re-use the code we have written for the taxi dataset to perform the spatial join and label our coordinates to the corresponding canton.
Number of WLAN observations for each canton
Using OpenStreetMaps we can define a rough Bounding Box of Lausanne to then grid the area
# Lausanne BB
min_long = 6.5472
max_long = 6.7023
min_lat = 46.4972
max_lat = 46.5570
The Area of the bounding box would be over 113 km^2
With grids of 500m x 500m we would have roughly 455 zones...
We would not use them all anyway, should we maybe use smaller grids and then remove the ones without any point inside?
The area delimited by the bounding box has been subdivided into grids 500m x 500m.