求可根据经纬度范围输出县/州边界图像的Python工具

方法1：基于GeoPandas + Matplotlib + Contextily（适配你已用的工具链） #

你可以在GeoPandas基础上，通过数据裁剪+局部绘图+卫星底图叠加实现需求，无需换工具：

1. 加载边界数据：可以用US Census Bureau的TIGER/Line shapefile（县/州级），或GeoPandas内置的naturalearth_lowres（精度足够时使用）
2. 用四角坐标裁剪数据：利用GeoPandas的cx索引快速筛选指定区域内的边界
3. 转换坐标适配卫星底图：卫星底图通常用Web Mercator（EPSG:3857），需把WGS84（EPSG:4326）的边界数据转换过去
4. 绘图并保存为图像

示例代码：

import geopandas as gpd
import matplotlib.pyplot as plt
import contextily as ctx

# 1. 加载县边界数据（以美国为例，替换为你需要的区域数据）
counties = gpd.read_file(gpd.datasets.get_path('naturalearth_lowres'))
# 若需更高精度，替换为US Census的shapefile路径：counties = gpd.read_file('tl_2023_us_county.shp')

# 2. 定义四角经纬度（示例：西经100，北纬40，西经95，北纬35）
min_lon, max_lat, max_lon, min_lat = -100, 40, -95, 35

# 3. 裁剪区域内的边界
local_counties = counties.cx[min_lon:max_lon, min_lat:max_lat]

# 4. 转换为Web Mercator坐标（适配卫星底图）
local_counties_3857 = local_counties.to_crs(epsg=3857)

# 5. 绘图
fig, ax = plt.subplots(figsize=(10, 10))
local_counties_3857.boundary.plot(ax=ax, linewidth=2, color='red')

# 6. 添加卫星底图
ctx.add_basemap(ax, source=ctx.providers.Esri.WorldImagery, crs=local_counties_3857.crs)

# 7. 设置视图范围（确保完全显示目标区域）
ax.set_xlim(local_counties_3857.total_bounds[0], local_counties_3857.total_bounds[2])
ax.set_ylim(local_counties_3857.total_bounds[1], local_counties_3857.total_bounds[3])

# 8. 关闭坐标轴并保存图像
ax.set_axis_off()
plt.savefig('local_county_boundaries.png', dpi=300, bbox_inches='tight')

方法2：用Folium生成交互式地图后导出静态图像 #

如果需要快速预览并导出，Folium支持直接添加边界和卫星底图，再通过工具导出为静态图：

1. 初始化Folium地图，设置中心点和缩放级别（基于四角坐标计算）
2. 添加县/州边界图层（GeoJson格式）
3. 设置卫星底图（如Esri World Imagery）
4. 用第三方工具（如selenium）导出为PNG

示例代码（导出需额外安装selenium）：

import folium
from selenium import webdriver
import time
import geopandas as gpd

# 定义四角坐标
min_lon, max_lat, max_lon, min_lat = -100, 40, -95, 35
center_lat = (max_lat + min_lat)/2
center_lon = (max_lon + min_lon)/2

# 初始化地图，设置卫星底图
m = folium.Map(
    location=[center_lat, center_lon], 
    zoom_start=8, 
    tiles='https://server.arcgisonline.com/ArcGIS/rest/services/World_Imagery/MapServer/tile/{z}/{y}/{x}', 
    attr='Esri'
)

# 添加边界图层
counties = gpd.read_file(gpd.datasets.get_path('naturalearth_lowres'))
local_counties = counties.cx[min_lon:max_lon, min_lat:max_lat]
folium.GeoJson(
    local_counties, 
    style_function=lambda x: {'color': 'red', 'weight': 2}
).add_to(m)

# 保存为HTML
m.save('local_map.html')

# 用Selenium导出为PNG
driver = webdriver.Chrome()
driver.get('file:///path/to/local_map.html')  # 替换为实际文件路径
time.sleep(2)  # 等待地图加载完成
driver.save_screenshot('local_boundary_satellite.png')
driver.quit()

方法3：基于Rasterio + GeoPandas（专业栅格叠加） #

如果需要更精细的栅格控制（比如调整透明度、分辨率），可以用Rasterio处理卫星底图，再叠加矢量边界：

1. 下载指定区域的卫星影像（或用Rasterio从在线源获取）
2. 用GeoPandas读取边界数据并裁剪到影像范围
3. 用Rasterio或Matplotlib将边界绘制到影像上并保存

关键注意事项 #

• 边界数据精度：如果naturalearth_lowres精度不足，直接使用对应国家的官方shapefile（如美国用US Census数据）
• 坐标系统：确保矢量和栅格数据的坐标一致，避免出现偏移
• 图像分辨率：保存时调整dpi参数，满足后续叠加卫星图像的需求

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