如何旋转Matplotlib地图？

Question

从https://s3.amazonaws.com/nyc-tlc/misc/taxi_zones.zip获取的shapefile，我希望绘制曼哈顿区，并为每个出租车区域绘制轮廓线。

这段代码将分别旋转每个单独的出租车区域，而不是一次性旋转所有的出租车区域。

import geopandas as gpd
from matplotlib import pyplot as plt


fname = "path_to_shapefile.shp"
df = gpd.read_file(fname)
df = df[df['borough'] == "Manhattan"]
glist = gpd.GeoSeries([g for g in df['geometry']])
glist = glist.rotate(90)
glist.plot()

[编辑] 我已经进一步完善了这个程序，使其能够以编程方式旋转图像。然而，如果我添加一个图例，那么图例也会旋转，这是不可取的。仍在寻找更好的解决方案。 请注意，还有这篇stackoverflow帖子（如何通过90度旋转matplotlib图来旋转图形？），但是，只能对图形进行90度旋转的解决方案不能对图像进行旋转。

import geopandas as gpd
from matplotlib import pyplot as plt

import numpy as np
from scipy import ndimage
from matplotlib import transforms


fname = "path_to_shapefile.shp"
df = gpd.read_file(fname)
df = df[df['borough'] == "Manhattan"]
df.plot()
plt.axis("off")
plt.savefig("test.png")

img = plt.imread('test.png')

rotated_img = ndimage.rotate(img, -65)
plt.imshow(rotated_img, cmap=plt.cm.gray)
plt.axis('off')
plt.show()

[编辑2]

@PMende下面给出的答案需要进行简单修改即可解决问题。

df = gpd.read_file(fname)
df = df[df['borough'] == "Manhattan"]
glist = gpd.GeoSeries([g for g in df['geometry']])
glist = glist.rotate(-65, origin=(0,0))
glist.plot()

关键在于将所有对象围绕一个单一点旋转，而不是围绕它们各自的起点旋转。

[EDIT 3] 如果有人正在尝试这样做，并且需要将旋转后的地理系列保存到数据框中（例如，基于附加列对几何图形进行着色），则需要创建一个新的数据框，只需简单写入即可。

df['geometry'] = glist

不起作用。我目前不确定原因。然而，下面的代码适用于我。

new_dataframe = gpd.GeoDataFrame(glist)
new_dataframe = new_dataframe.rename(columns={0:'geometry'}).set_geometry('geometry')
new_dataframe.plot()

- Connor

---

请记住，您可能不希望旋转几何图形，因为这会使其失去地理意义。您只需要旋转 matplotlib 绘图即可。这可能有助于搜索。 - martinfleis

@martinfleis 谢谢您的建议 - 那确实有帮助。我找到了一种旋转整个图像的方法，但如果我添加一个图例，那也会被旋转，这是一个问题。 - Connor

我知道这是一个老问题，但是否有一种方法可以添加底图（例如使用Contextily）？我在这里有一个类似的问题：https://stackoverflow.com/questions/65498106/how-to-rotate-a-contextily-basemap-in-matplotlib-and-jupyter-notebook - a11

Answer 1

如果我正确理解了GeoPandas的文档，您可以指定每个几何体旋转的起始点（默认为每个几何体的中心）。要获得所需的行为，您可以将每个形状绕相同的起始点旋转。

例如：

import geopandas as gpd
from matplotlib import pyplot as plt


fname = "path_to_shapefile.shp"
df = gpd.read_file(fname)
df = df[df['borough'] == "Manhattan"]

center = df["geometry"].iloc[0].centroid()
glist = gpd.GeoSeries([g for g in df['geometry']])
glist = glist.rotate(90, origin=center)

glist.plot()

我自己无法进行测试，但它应该可以帮助你朝着正确的方向开始。

（尽管我也同意 @martinfeleis 的观点，不一定要旋转几何图形，而是旋转绘图。）

- PMende

---

你甚至不需要找到中心 -- 实际上，你只需要把原点(0,0)作为旋转的中心点就可以了，这样也可以正常工作。谢谢你的帮助！ - Connor

我知道这是一个老问题，但是否有一种方法可以添加基础地图（例如使用Contextily）？我在这里有一个类似的问题：https://stackoverflow.com/questions/65498106/how-to-rotate-a-contextily-basemap-in-matplotlib-and-jupyter-notebook - a11