如何结合`GeoDataFrame`和`add_feature()`来限定绘制范围？

结合`GeoDataFrame`和`add_feature()`来限定绘制范围的一般步骤如下：

1. 首先，你需要有一个包含地理位置信息（如经纬度）的`GeoDataFrame`。这通常来自于CSV文件或者一些地理空间数据库，例如`geopandas`库读取的数据。

2. 创建`GeoDataFrame`时，它的坐标系统默认可能会是WGS84或者其他投影系统。你可以使用`to_crs()`函数将其转换为与你的地图投影一致的坐标系统，以便更好地控制绘制范围。

3. 确定你要在地图上显示的具体经纬度范围。这通常是基于你的研究兴趣或者数据特点，例如`(start_longitude, start_latitude, end_longitude, end_latitude)`。

4. 使用`GeoDataFrame`的`bounds`属性获取其最小和最大坐标值，如果需要，可以调整这个范围以适应实际需求。

5. 利用`GeoAxes`的`set_extent()`方法，将经纬度范围作为参数传递进去，例如`ax1.set_extent(geo_df.bounds)`. 这会自动根据`GeoDataFrame`的经纬度范围更新地图视窗。

6. 最后，在`add_feature()`前调用`set_extent()`，即可保证只有指定范围内的地理特征会被绘制出来。

示例代码如下：

import geopandas as gpd

# 加载GeoDataFrame
gdf = gpd.read_file('your_shapefile.shp')

# 转换坐标系统并设置范围
gdf = gdf.to_crs(epsg=3857)  # 或者你的目标投影
min_lon, min_lat, max_lon, max_lat = gdf.total_bounds

# 更新绘制范围
ax1.set_extent([min_lon, max_lon, min_lat, max_lat])

# 添加地理特征
ax1.add_feature(gdf, ...)