Берёшь точку.
Перебираешь твой словарь.
Переводишь значения из строки в два float.
Вычисляешь расстояние между твоей точкой и текущей точкой в словаре.
Пишешь в какой-нибудь список кортеж из результата функции (расстояние между точками) и ключа текущей точки.
Сортируешь список по первому значению кортежа.
Берёшь слайс первых четырёх значений.
PROFIT.
Функция вот.
from math import radians, cos, sin, asin, sqrt
def haversine(lat1, lon1, lat2, lon2):
"""
Вычисляет расстояние в километрах между двумя точками, учитывая окружность Земли.
https://en.wikipedia.org/wiki/Haversine_formula
"""
# convert decimal degrees to radians
lon1, lat1, lon2, lat2 = map(radians, (lon1, lat1, lon2, lat2))
# haversine formula
dlon = lon2 - lon1
dlat = lat2 - lat1
a = sin(dlat / 2) ** 2 + cos(lat1) * cos(lat2) * sin(dlon / 2) ** 2
c = 2 * asin(sqrt(a))
km = 6367 * c
return km