Minimalistic geospatial searching solution for Python

Question

I am looking for a minimalistic solution for doing basic geospatial search in Python. We have a dataset of roughly 10 k locations and we need to solve the find the all locations within a radius of N kilometers from a given location. I am not looking for explicit database with geospatial support. I hope to get around another external solution. Is there something that would use Python only?

Answer 1

Shapely seems to be a good solution. Its description seems to correspond to what you're looking for :

[Shapely] It lets you do PostGIS-ish stuff outside the context of a database using Python.

It is based on GEOS, which a widely used C++ library.

Here is a link to the documentation

Answer 2

scipy.spatial has a kd-tree implementation that might be the most popular in Python.