Python'da en yakın enlem ve boylam noktaları

0

Soru

Liste 1 ve Liste 2'deki istasyonların bir listesi var Liste 2'den Liste 1'e en yakın istasyonları nasıl bulabilirim ?

liste 1 ve 2

Liste 1

SS Enlem Boylam Yok 977 23.141747 53.796469 946 23.398398 55.422916 742 23.615732 53.717952 980 23.633077 55.567046 660 23.6504 54.4007

Liste 2

SS Enlem Boylam Yok 962 23.657571 53.703683 745 23.671971 52.955976 743 23.766849 53.770344 978 23.847163 52.809653 748 23.942166 52.16236 744 23.955817 52.790424 760 23.984592 55.55764 945 24.030256 55.844842 894 24.03511 53.891547 856 24.741601 55.80063 893 24.04123 53.899958 387 24.059988 51.748138 675 24.061578 53.417912 664 24.063978 51.76195

Bunu Powerbi'de eşleştirerek manuel olarak yapabilirim, ancak ölçeklenebilir bir çözüm arıyorum ve Python'u tercih et.

geolocation python
2021-11-16 11:26:01
2

En iyi cevabı

0

geopy senin arkadaşın. Koordinat çiftleri arasındaki mesafeyi hesaplayan fonksiyonlara sahiptir. İşte bir yaklaşım:

from geopy import distance

s = {
    977: (23.141747, 53.796469),
    946: (23.398398, 55.422916),
    # etc etc
}

d = {
    962: (23.657571, 53.703683),
    745: (23.671971, 52.955976),
    743: (23.766849, 53.770344),
    # etc etc
}

for (ss, a) in s.items():
    best = None
    dist = None
    for (dd, b) in d.items():
        km = distance.distance(a, b).km
        if dist is None or km < dist:
            best = dd
            dist = km

    print(f'{ss} is nearest {best}: {dist} km')

Örnek verilerinizle çalıştırırsam elde ederim:

977 is nearest 962: 57.909274476145846 km
946 is nearest 760: 66.3613771288571 km
742 is nearest 962: 4.857141728990575 km
980 is nearest 760: 38.94400553289674 km
660 is nearest 743: 65.56437761273963 km

Liste 1'in her birinden liste 2'nin her birine olan mesafelere ihtiyacınız varsa, deneyebilirsiniz

sorted([[ss, dd, distance.distance(a, b).km] for (ss, a) in s.items()
        for (dd, b) in d.items()])

Bu, önce liste 1 öğesi tarafından sipariş edilen listelerin bir listesini verir:

[[660, 387, 273.98088337893], 
 [660, 664, 272.6633222300461], 
 [660, 675, 109.98235440892797], 
 [660, 743, 65.56437761273963],

 # etc

Daha fazla sıralama veya gruplama, okuyucu için bir alıştırma olarak bırakılmıştır.

2021-11-17 12:11:25

Teşekkürler! ama birden fazla istasyonu nasıl yaklaştırabilirim ? Liste 2'deki en iyi 2 veya 3 en yakın istasyonun liste 1'e olmasını ister miyim?
Raul V

Örneğin: 977 962'ye en yakın,: 57.909274476145846 km 760,: 57.909274476145846 km xyz: 57.909274476145846 km a
Raul V

Sana bir başlangıç yapmak için cevabımı düzenledim.
xpqz

sen bir canavarsın. Çok teşekkür ederim.
Raul V
0

Bu önceki soruya oldukça benzer
Enlem / boylam temelinde iki nokta arasındaki mesafeyi alma
yani, kopya olarak talep edilebilir.
Her neyse, Kurt Peek'in cevabını takiben, yapabilirsin:

import geopy.distance

def get_distnace(coords_1, coords_2):
    return geopy.distance.vincenty(coords_1, coords_2).km

mesafeyi 'km'olarak döndüren bir işlev elde etmek için. Daha sonra, sizinki gibi birkaç liste verildiğinde, formda

list1 = [[stat_name_1, lat_1, lon_1], [stat_name_1, lat_1, lon_1], ... ]

list2 = [[...], ... ]

muhtemelen yapabilirsin.:

min_stat = get_distance(list1[0][1:], list2[0][1:])
for stat_1 in list1:
    coord_1 = stat_1[1:]
    for stat_2 in list2:
        coord_2 = stat_2[1:]
        min_stat = min(min_stat, get_distance(coord_1, coord_2)

Herhangi bir uygulamadan önce hata ayıklanması ve test edilmesi için taslak, fikir olarak alın.

2021-11-16 12:03:54

Diğer dillerde

Bu sayfa diğer dillerde

Русский
..................................................................................................................
Italiano
..................................................................................................................
Polski
..................................................................................................................
Română
..................................................................................................................
한국어
..................................................................................................................
हिन्दी
..................................................................................................................
Français
..................................................................................................................
Česk
..................................................................................................................
Português
..................................................................................................................
ไทย
..................................................................................................................
中文
..................................................................................................................
Español
..................................................................................................................
Slovenský
..................................................................................................................