Расчет расстояния между городами и поиск близлежащих городов на основе GeoPT, в Python на Google App Engine

У меня есть определенная модель городов, которая сохраняет geoname_idи местоположение(как GeoPt) город. Есть две вещи, которых я хочу достичь.

1. Я хочу получить все города в радиусе 500 кмот заданного города.
2. Я хочу рассчитать расстояние в кммежду двумя заданными городами.

Как лучше всего это сделать, учитывая, что у меня очень большая база данных городов, и я не хочу сильно жертвовать факторами производительности. Любая помощь или совет приветствуются.