Недавно я занимался некоторой статистикой на основе OpenStreetMaps. Когда я использовал эстакад для экспорта немецкого Autobahnen и Geopandas для расчета их длины, я обнаружил, что длины не соответствуют официальным цифрам. Перепускной запрос уже исключает соединительные рампы, выходы, (...).
сравнивая самые длинные дюжины из них, большинство из них явно переоценили, поскольку цифры были на 50/60% выше, чем официальные. (См. Image)
. />[out:json][timeout:100];
{{geocodeArea:Deutschland}}->.searchArea;
way["highway"="motorway"]["ref"="A 7"](area.searchArea);
out body;
>;
out skel qt;
code (python)
import geopandas as gpd
import pandas as pd
autobahn = "A7"
columns = ["ref", "length_m"]
df_total = pd.DataFrame(columns=columns)
gdf = gpd.read_file(f"https://fliessbaden.de/wp-content/uploads/A7.geojson")
gdf = gdf.to_crs(epsg=3857)
gdf["length_m"] = gdf.geometry.length
total_m = gdf['length_m'].sum()
total_km = total_m / 1000
print(autobahn + ": " + str(round(total_km, 2)) + " km")
Подробнее здесь: https://stackoverflow.com/questions/796 ... -distances
Геопандас переоценивает расстояния? ⇐ Python
-
- Похожие темы
- Ответы
- Просмотры
- Последнее сообщение