Запитання з тегом «geopandas»

Це здається досить простим питанням, але я не можу зрозуміти, як перетворити пандус DataFrame в GeoDataFrame для просторового з'єднання. Ось приклад того, як виглядають мої дані df.head(): Date/Time Lat Lon ID 0 4/1/2014 0:11:00 40.7690 -73.9549 140 1 4/1/2014 0:17:00 40.7267 -74.0345 NaN Насправді цей фрейм даних був створений з …

42 csv  geopandas  pandas 

У мене є 2 геодетарам: import geopandas as gpd from shapely.geometry import Point gpd1 = gpd.GeoDataFrame([['John',1,Point(1,1)],['Smith',1,Point(2,2)],['Soap',1,Point(0,2)]],columns=['Name','ID','geometry']) gpd2 = gpd.GeoDataFrame([['Work',Point(0,1.1)],['Shops',Point(2.5,2)],['Home',Point(1,1.1)]],columns=['Place','geometry']) і я хочу знайти ім'я найближчої точки в gpd2 для кожного рядка в gpd1: desired_output = Name ID geometry Nearest 0 John 1 POINT (1 1) Home 1 Smith 1 POINT …

20 python  distance  geopandas  pandas 

У мене простий DataPreme GeoPandas: Я хотів би завантажити цей GeoDataframe в таблицю PostGIS. У мене вже встановлена ​​база даних з розширенням PostGIS, але не можу додати цей Dataframe як таблицю. Я спробував таке: engine = <> meta = MetaData(engine) eld_test = Table('eld_test', meta, Column('id', Integer, primary_key=True), Column('key_comb_drvr', Text), Column('geometry', …

17 python  postgis  geopandas 

Враховуючи geopandas GeoDataFrameсерію полігонів, я хотів би отримати площу в км кв кожної функції в моєму списку. Це досить поширена проблема, і звичайним запропонованим рішенням у минулому було використання shapelyта pyprojбезпосередньо (наприклад, тут і тут ). Чи є спосіб зробити це по-чистому geopandas?

16 python  area  shapely  geopandas  pyproj 

У мене є список фреймів даних, всі однакові за структурою і всі є багатокутниками. Як я можу додати їх усі разом за допомогою геопанд?

16 python  geopandas 

Я намагався створити функцію, яка в основному робить те саме, що QGIS "розчиняє" функцію. Я подумав, що це буде дуже просто, але, мабуть, ні. Тож із того, що я зібрав навколо, найкращим варіантом тут має бути використання фіона з витонченою формою. Я тільки почав возитися з векторними файлами, так що …

15 python  shapely  fiona  geopandas 

Я намагаюся прив’язати дві лінії один до одного за допомогою Shapely / Geopandas, але результат оснащення дуже дивний. Я намагався : import geopandas as gpd from shapely.geometry import * from shapely.ops import snap lines1 = gpd.GeoDataFrame.from_file('lines1.shp') lines1 = lines1.to_crs({'init': 'epsg:2227'}) lines2 = gpd.GeoDataFrame.from_file('lines2.shp') lines2 = lines2.to_crs({'init': 'epsg:2227'}) res = lines1 …

14 shapely  snapping  geopandas 

Я використовую наведений нижче код, щоб знайти країну (а іноді й штат) для мільйонів GPS-точок. На даний момент код займає близько однієї секунди за точку, що надзвичайно повільно. Форма файлу - 6 Мб. Я читав, що геопанди використовують rtrees для просторових приєднань, що робить їх неймовірно ефективними, але це, здається, …

14 python  spatial-join  geopandas  rtree 

Чи можна читати необроблені дані в a geopandas GeoDataFrame, a la a pandas DataFrame? Наприклад, такі роботи: import pandas as pd import requests data = requests.get("https://data.cityofnewyork.us/api/geospatial/arq3-7z49?method=export&format=GeoJSON") pd.read_json(io.BytesIO(r.content)) Наступне не: import geopandas as gpd import requests data = requests.get("https://data.cityofnewyork.us/api/geospatial/arq3-7z49?method=export&format=GeoJSON") gpd.read_file(io.BytesIO(r.content)) Іншими словами, чи можна зчитувати геопросторові дані, що знаходяться в пам'яті, …

14 python  fiona  geopandas  pandas 

У мене є геопанди, GeoDataFrame що містять сотні фігурнихPolygon та MultiPolygonгеометрій. Багатокутники перекриваються в багатьох місцях. Я хотів би скласти нову геометрію, яка містить підрахунок кількості перекриттів. Щось на зразок цього: Хтось має ідеї, як підійти до цього? Я навіть не можу побачити вихід. Врешті-решт, я особливо хотів би мати …

13 python  shapely  geopandas  pandas 

це виглядає як дурне запитання, але я не можу знайти чітку відповідь на це: які одиниці геопанд / формально використовують розрахунок відстані / площі між об'єктами?

13 python  shapely  geopandas 

Я маю GeoPandas 0.2.1 py27_0 fiona 1.7.0 np110py27_1 Пітон 2.7.10 встановлений у середовищі Anaconda 2-4.1.1-Windows-x86. Я вмію конструювати GeoDataFrame, читаючи набори вхідних даних та керувати ними, але збереження вихідного набору даних не зберігає систему координат. import geopandas as gpd world = gpd.read_file(gpd.datasets.get_path('naturalearth_lowres')) world.to_file(driver='ESRI Shapefile',filename=r'C:\GIS\Temp\world_out.shp') Файл world_out.shpмістить особливості та схему вихідного …

13 python  fiona  geopandas  pandas  anaconda 

У мене є географічний каркас даних із ряду рядків рядків, у яких є деякі дані, пов'язані з кожною вершиною / точкою: Point_x = (Lat, Lon, Time, ID, Data1, Data2, Data3) Точки перетворюються у рядки на основі ідентифікатора та упорядковуються часом. Я хочу розбити рядки рядків там, де виконується певна умова. …

11 python  shapely  geopandas 

Я маю географічний кадр даних в EPSG: 4326, і я би створив новий кадр даних, що складається з усіх рядків, які потрапляють у певне обмежувальне поле. Спочатку я дістаю обмежувальне поле, про яке мені важливо (що насправді є обмежувальним вікном іншого фрейму даних): print df_sussex.total_bounds [ -1.57239292 50.57467674 0.14528384 51.27465152] …

11 python  geopandas 

Я намагаюся знайти, де кілька ліній перетинають полігон для двох різних геоданих кадрів: from shapely.geometry import Polygon, LineString import geopandas as gpd polygon = Polygon([(0, 0), (1, 0), (1, 1), (0, 1), (0, 0)]) line1 = LineString([(0.5, 0.5), (0.7, 0.7)]) line2 = LineString([(0.9, 0.9), (0.2, 0.6)]) poly_gdf = gpd.GeoDataFrame(geometry=[polygon]) line_gdf …

11 python  intersection  shapely  geopandas