如何通过经纬度坐标限制Geopandas绘图的显示范围？

当然可以！你完全可以通过指定经纬度范围来实现Geopandas地图的局部放大，这里有两种实用的方法，根据你的需求选择即可：

方法1：直接设置绘图轴的范围 #

Geopandas的plot()方法会返回一个matplotlib的Axes对象，你可以利用这个对象的set_xlim()和set_ylim()方法直接限定经度（x轴）和纬度（y轴）的显示范围。

示例代码：

import geopandas as gpd
import matplotlib.pyplot as plt

# 读取shapefile
map_gdf = gpd.read_file('path_to_shapefile')

# 目标区域的经纬度范围（示例：纽约市区大致范围）
min_lon, max_lon = -74.05, -73.9
min_lat, max_lat = 40.6, 40.8

# 绘图并获取轴对象
fig, ax = plt.subplots(figsize=(10, 8))
map_gdf.plot(ax=ax)

# 设置经纬度范围
ax.set_xlim([min_lon, max_lon])
ax.set_ylim([min_lat, max_lat])

plt.show()

这种方法的好处是保留所有要素，只是只显示指定区域，适合你需要快速查看局部区域的场景。

方法2：先裁剪GeoDataFrame再绘图 #

如果你希望只绘制目标区域内的要素（提升绘图效率，尤其是处理大型shapefile时），可以用Geopandas的cx索引器直接裁剪GeoDataFrame，只保留落在指定经纬度范围内的要素，再进行绘图。

示例代码：

import geopandas as gpd

# 读取shapefile
map_gdf = gpd.read_file('path_to_shapefile')

# 目标区域的经纬度范围
min_lon, max_lon = -74.05, -73.9
min_lat, max_lat = 40.6, 40.8

# 裁剪GeoDataFrame：cx[经度范围, 纬度范围]
map_gdf_cropped = map_gdf.cx[min_lon:max_lon, min_lat:max_lat]

# 绘制裁剪后的地图
map_gdf_cropped.plot(figsize=(10, 8))

重要注意事项 #

• 坐标系匹配：上述方法仅当你的GeoDataFrame使用的是**WGS84坐标系（EPSG:4326）**时有效（也就是常用的经纬度格式）。如果你的shapefile是其他坐标系（比如UTM），需要先转换坐标系：

  map_gdf = map_gdf.to_crs(epsg=4326)
  

• 经纬度范围的顺序：注意cx索引器是先写经度范围（x轴），再写纬度范围（y轴），不要搞反哦！

内容的提问来源于stack exchange，提问作者Jason