我目前有一个函数:
def cross_country(contiguous_usa, full_geo_data):
full_geo_data['Coordinates'] = full_geo_data[['longitude', 'latitude']].values.tolist()
full_geo_data['Coordinates'] = full_geo_data['Coordinates'].apply(Point)
full_geo_data = gpd.GeoDataFrame(full_geo_data, geometry='Coordinates')
fig = plt.figure(figsize=(10,15))
ax1 = plt.subplot(212, projection=gcrs.AlbersEqualArea(central_latitude=-98, central_longitude=39.5))
gplt.kdeplot(full_geo_data[full_geo_data['speedkmh'] == 0],
projection=gcrs.AlbersEqualArea(), cmap="cool", clip=contiguous_usa.geometry, ax=ax1)
gplt.polyplot(contiguous_usa, projection=gcrs.AlbersEqualArea(), ax=ax1)
plt.title("Test")
cross_country(contiguous_usa, full_geo_data)然而,它工作得很好,当我运行它时,地图是这样的:

我知道这是一件微不足道的事情,但我已经查阅了文档,我找不到任何与改变基本方向有关的东西,因为它只是简单地旋转绘图。
发布于 2019-03-13 04:13:23
最简单的答案可能涉及到您正在使用的投影或crs。然而,如果你不能让它工作,你可以使用shapely来修改geodataframe的单独行。
def rotator(row):
row['geometry'] = shapely.affinity.rotate(row['geometry'], -90)
return row
full_geo_data = full_geo_data.apply(rotator, axis = 1)https://stackoverflow.com/questions/54962758
复制相似问题