Impressed how fast GeoPandas spatial join based on distance (GeoDataFrame.sjoin nearest) is. 1000 locations tested against 90k locations finished in 9 seconds on my Macbook. Now to get ChatGPT to write this code more elegantly.