Nejbližší zeměpisné šířky a délky body v pythonu

0

Otázka

Mám seznam stanic, Seznam 1 a Seznam 2 Jak mohu najít nejbližší stanice do Seznamu 1 nebo v Seznamu 2 ?

seznam 1 a 2

Seznam 1

SS Č. Latitude Longtitude 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

Seznam 2

SS Č. Latitude Longtitude 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

Mohu to udělat ručně tím, že jejich mapování na PowerBI, ale mám zájem o škálovatelné řešení, a Přednost Pythonu.

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

Nejlepší odpověď

0

geopy je váš přítel. To má funkce, které výpočet vzdálenosti mezi dvojicemi souřadnic. Zde je jeden přístup:

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')

Když jsem jej spustit s vaší příklad dat:

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

Pokud potřebujete vzdálenosti od každého z seznam 1, seznam 2, můžete zkusit

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

To dává seznam seznamů, objednat podle seznamu 1 položka první:

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

 # etc

Další objednávce nebo seskupení odešel jako exercice pro čtenáře.

2021-11-17 12:11:25

Díky!!! ale jak mohu získat více stanic poblíž ? Pokud chci, horní 2 nebo 3 nejbližší stanice v seznam 2 seznam 1?
Raul V

Například : 977 je nejbližší 962,: 57.909274476145846 km 760,: 57.909274476145846 km xyz: 57.909274476145846 km
Raul V

Upravil jsem svou odpověď, aby vám začít.
xpqz

ty jsi ZVÍŘE. Děkuji moc.
Raul V
0

To je docela podobná předchozí otázku
Získání vzdálenosti mezi dvěma body na základě zeměpisné šířky/délky
tak, to by mohla být uplatněna jako duplicitní.
Každopádně, následující Kurt Nahlédnout odpověď, jsi to mohl udělat:

import geopy.distance

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

získat funkci, která vrátí vzdálenost v "km". Pak, vzhledem k tomu, pár seznamů, jako je vaše, ve formě

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

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

to by se asi dalo:

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)

Berte to jako návrh, nápad, ladit a zkoušet, než jakékoliv aplikace.

2021-11-16 12:03:54

V jiných jazycích

Tato stránka je v jiných jazycích

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