使用 geopandas 和 geoplot 在纽约地图上绘制位置

编程问答 2022-05-22

问题描述

我有一个商店位置文件，我试图将其绘制到纽约市地图上。到目前为止，我一直在使用以下两个链接作为我的指南，但我一直无法使代码正常工作。

import geopandas
import geoplot
import pandas as pd
from shapely.geometry import Point,polygon

store_locDF=pd.read_csv('stores.csv') #Import file
cols_to_keep=['store_longitude','store_latitude']
store_locDF = store_locDF[store_locDF.columns.intersection(cols_to_keep)] #Keep only long and lat

#Convert long and lat to "points"
geometry = [Point(xy) for xy in zip( store_locDF["store_longitude"],store_locDF["store_latitude"])]
#Store geometry in dataframe
store_locDF['geometry']=geometry

#print(store_locDF.head())

boroughs = geopandas.read_file(geoplot.datasets.get_path('nyc_boroughs'))

ax=geoplot.polyplot(boroughs,edgecolor='white',facecolor='lightgray',figsize=(12,8))

geoplot.polyplot(geometry,ax=ax) #<--Didnt work
#geoplt.pointplot(geometry,ax=ax) <--Also didnt work
''''


The error I am getting is: 
ValueError                                Traceback (most recent call last)
<ipython-input-12-3c7d64229a6c> in <module>
     19 ax=geoplot.polyplot(boroughs,8))
     20 
---> 21 geoplot.polyplot(geometry,ax=ax)

~\anaconda3\envs\IST 652\lib\site-packages\geoplot\geoplot.py in polyplot(df,projection,extent,figsize,ax,**kwargs)
    917             return ax
    918 
--> 919     plot = polyPlot(df,figsize=figsize,ax=ax,extent=extent,projection=projection,**kwargs)
    920     return plot.draw()
    921 

~\anaconda3\envs\IST 652\lib\site-packages\geoplot\geoplot.py in __init__(self,df,**kwargs)
    881     class polyPlot(Plot):
    882         def __init__(self,**kwargs):
--> 883             super().__init__(df,**kwargs)
    884 
    885         def draw(self):

~\anaconda3\envs\IST 652\lib\site-packages\geoplot\geoplot.py in __init__(self,**kwargs)
    614             # self when it is asked for its geometry property,and so it will never be the source
    615             # of this error.
--> 616             raise ValueError(
    617                 'The input GeoDataFrame does not have a "geometry" column set.'
    618             )

ValueError: The input GeoDataFrame does not have a "geometry" column set.


Any help is appreciated

解决方法

正如评论者“Paul H”所说，问题在于我提供的是列表而不是地理数据框。以下是工作代码：

import geopandas
import geoplot
import pandas as pd
from shapely.geometry import Point,Polygon
import matplotlib.pyplot as plt

store_locDF=pd.read_csv('stores.csv')
cols_to_keep=['store_longitude','store_latitude']

crs={'init':'epsg:4326'}

boroughs = geopandas.read_file(geoplot.datasets.get_path('nyc_boroughs'))

geo_df=geopandas.GeoDataFrame(store_locDF,crs=crs,geometry=geopandas.points_from_xy(store_locDF["store_longitude"],store_locDF["store_latitude"]))



fig,ax=plt.subplots(figsize=(15,15))
boroughs.plot(ax=ax,alpha=0.4,color="grey")
geo_df[geo_df['active_yn']=="Y"].plot(ax=ax,markersize=200,color="green",label="Active Store")
geo_df[geo_df['active_yn']=="N"].plot(ax=ax,color="red",label="Inactive Stores")
plt.legend()
geopandasgeoplotpointpointpythonshapely