Geopy / как найти наиболее близкую точку из списка?

Question

[Сергей Ильин]

есть два словаря:

income = {'': 'Error', 'Kaliningrad (Konigsberg), ': (54.710454, 20.512733), 'Krasnojarsk, Тельмана, 36': (56.05031, 92.967305), 'Moscow > Moskva, ': (55.755814, 37.617635), 'Moscow > Moskva, Багратионовский пр. д.7 кор.20Б': (55.743137, 37.504232)}
# в income могут быть строки!

etalon = {'115035, г. Москва, М-35, ул. Балчуг, 2': (55.747459, 37.625271), '308000, г. Белгород, пр-т Славы, 74': (50.598079, 36.585839), '241050, г. Брянск, ул. Горького, 34': (53.244984, 34.365707), '600000, г. Владимир, ул. Большая Московская, 29': (56.129745, 40.406446)}

x_values = income.values()
y_values = etalon.values()

хочу циклом пройтись по income, взять геопару (широта, долгота) и сравнить с каждой геопарой в etalon. результатом должен быть список вида:
[адрес в income, адрес в etalon]
или словарь:
{адрес в income: адрес в etalon}

from geopy.distance import geodesic

for x in x_values:
    if type(x) == str:
        None
    else:
        distance = 10000000 # выбираем какую-то от балды дистанцию побольше
        for y in y_values():
            try:
                distance_new = geodesic(x, y, ellipsoid='WGS-84')
                distance_new = geodesic(x_values, y_values, ellipsoid='WGS-84')
                    if distance_new < distance:
                        distance = distance_new
                        total = {income.keys(): etalon.keys()}
            except:
                print('Error')

не работает ((

Comments

[Сергей Ильин]

Итоговое решение стало таким:

total_all = []
errors = []

for x, val_x in income.items():
    distance = 10000000
    if type(val_x) == str:
        print(x,val_x)
        errors.append((x,val_x))
    else:
        compare_addresses = ()
        for i in range(len(income)):
            for y, val_y in etalon.items():
                distance_new = geodesic(val_x, val_y, ellipsoid='WGS-84').km
                if distance_new < distance:
                    distance = distance_new
                    compare_addresses = (x, y, distance)

        print(compare_addresses)
        total_all.append(compare_addresses)

итого в errors мы получим список с отбракованными адресами, в total_all - список вида

[адрес на вход, ближайший найденный адрес из списка, расстояние]

Всем спасибо.
Правда, работает долго ((

Answer 1

[Дмитрий]

dict.keys() возвращает Dictionary view objects (dict_keys), а этот объект мутабельный и значит не может быть ключом в словаре.

Answer 2

[zexer]

1. Переменная y_values, а в цикле for вы пишите y_values()
2. Блок try имеет сильно большой отступ
3. Покажите ошибку.

Comments

[Сергей Ильин]

виноват, поправил.
(чистил цикл перед тем, как запостить сюда)

чтобы показать ошибку, убрал try / except:

for x in x_values:
    if type(x) == str:
        None
    else:
        distance = 10000000
        for y in y_values:
            distance_new = geodesic(x, y, ellipsoid='WGS-84')
            if distance_new < distance:
                distance = distance_new
                total = {income.keys(): etalon.keys()}
print(total)

TypeError: unhashable type: 'dict_keys'
total = {dict_x.keys(): dict_y.keys()}

очевидно, я не правильно возвращаю .keys()

[zexer]

Насколько я вижу, в словаре income множество ключей, однако вы пытаетесь создать новый словарь, и ключом в нем делаете множество ключей из словаря income. Ключом в новом словаре должно быть единственное значение, а не множество ключей из другого словаря.

[Сергей Ильин]

zexer,
то есть, обход по словарю мне нужен в духе
for k in dic:
print(k)
print(dic[k])

или

for k, val in dic.items():
print(k)
print(val)

т.е. не ссылаться сразу на dic.values(), а брать либо dic и в нем dic[key], пару "ключ-значение" в dic.items()?

[zexer]

Сергей Ильин, Получается, что вы правильно понимаете)

Answer 3

[freeExec]

Не стоит изобретать велосипед. Давно существуют пространственные индексы и поиск ближайшего.
Гуглить:

K nearest neighbors (kNN)

Comments

[Сергей Ильин]

kNN - ok