Aqui está meu primeiro geodatframe:
!pip install geopandas
import pandas as pd
import geopandas
city1 = [{'City':"Buenos Aires","Country":"Argentina","Latitude":-34.58,"Longitude":-58.66},
{'City':"Brasilia","Country":"Brazil","Latitude":-15.78 ,"Longitude":-70.66},
{'City':"Santiago","Country":"Chile ","Latitude":-33.45 ,"Longitude":-70.66 }]
city2 = [{'City':"Bogota","Country":"Colombia ","Latitude":4.60 ,"Longitude":-74.08},
{'City':"Caracas","Country":"Venezuela","Latitude":10.48 ,"Longitude":-66.86}]
city1df = pd.DataFrame(city1)
city2df = pd.DataFrame(city2)
gcity1df = geopandas.GeoDataFrame(
city1df, geometry=geopandas.points_from_xy(city1df.Longitude, city1df.Latitude))
gcity2df = geopandas.GeoDataFrame(
city2df, geometry=geopandas.points_from_xy(city2df.Longitude, city2df.Latitude))Cidade1
City Country Latitude Longitude geometry
0 Buenos Aires Argentina -34.58 -58.66 POINT (-58.66000 -34.58000)
1 Brasilia Brazil -15.78 -47.91 POINT (-47.91000 -15.78000)
2 Santiago Chile -33.45 -70.66 POINT (-70.66000 -33.45000)e meu segundo geodataframe: City2:
City Country Latitude Longitude geometry
1 Bogota Colombia 4.60 -74.08 POINT (-74.08000 4.60000)
2 Caracas Venezuela 10.48 -66.86 POINT (-66.86000 10.48000)eu gostaria do terceiro dataframe com a cidade mais próxima de city1 a city2 com a distância como:
City Country Latitude Longitude geometry Nearest Distance
0 Buenos Aires Argentina -34.58 -58.66 POINT (-58.66000 -34.58000) Bogota 111 KmAqui está minha solução atual usando geodjango e dict (mas é muito longo):
from django.contrib.gis.geos import GEOSGeometry
result = []
dict_result = {}
for city01 in city1 :
dist = 99999999
pnt = GEOSGeometry('SRID=4326;POINT( '+str(city01["Latitude"])+' '+str(city01['Longitude'])+')')
for city02 in city2:
pnt2 = GEOSGeometry('SRID=4326;POINT('+str(city02['Latitude'])+' '+str(city02['Longitude'])+')')
distance_test = pnt.distance(pnt2) * 100
if distance_test < dist :
dist = distance_test
result.append(dist)
dict_result[city01['City']] = city02['City']Aqui estão minhas tentativas:
from shapely.ops import nearest_points
# unary union of the gpd2 geomtries
pts3 = gcity2df.geometry.unary_union
def Euclidean_Dist(df1, df2, cols=['x_coord','y_coord']):
return np.linalg.norm(df1[cols].values - df2[cols].values,
axis=1)
def near(point, pts=pts3):
# find the nearest point and return the corresponding Place value
nearest = gcity2df.geometry == nearest_points(point, pts)[1]
return gcity2df[nearest].City
gcity1df['Nearest'] = gcity1df.apply(lambda row: near(row.geometry), axis=1)
gcity1dfaqui :
City Country Latitude Longitude geometry Nearest
0 Buenos Aires Argentina -34.58 -58.66 POINT (-58.66000 -34.58000) Bogota
1 Brasilia Brazil -15.78 -70.66 POINT (-70.66000 -15.78000) Bogota
2 Santiago Chile -33.45 -70.66 POINT (-70.66000 -33.45000) BogotaSaudações
Olá e bem-vindo ao StackOverflow! Você parece ter a impressão de que o StackOverflow é um site em que você publica um problema e recebe algum código em troca. De fato, não é esse o caso. Sua pergunta provavelmente será encerrada ou mesmo excluída em breve. Para impedir que isso aconteça no futuro, faça o tour e dê uma olhada no centro de ajuda . Em particular, torne-se familiar quanto ao que é considerado tópico por aqui
—
azro
Além disso, quando você postar sobre DF, poste código pyhton com o conteúdo do DF, para todas as pessoas que gostariam de ajudá-lo a não escrever por conta própria
—
azro
@ azro eu editei e adicionei minha solução ao problema e aos meus dados iniciais.
—
user462794
são suas cidades apenas na América do Sul? Caso contrário, a que distância eles podem estar um do outro? Quantas cidades podem estar na cidade1 e quantas na cidade2? É importante encontrar a solução mais rápida ou uma solução mais simples, executada em um tempo razoável, OK? Se este for o caso, qual seria um tempo razoável?
—
Walter Tross
@ WalterTross minha cidade está em todo o mundo e estou procurando a solução fastet. Graças
—
user462794