手把手教你用Python调用天地图WMS/WFS服务，5分钟获取地理数据并可视化

Python实战5分钟调用天地图WMS/WFS服务实现地理数据可视化第一次接触地理信息服务时我被各种缩写搞晕了——WMS、WFS、WMTS这些OGC标准接口文档读起来像天书。直到用Python实际调用了一次天地图服务才发现原来获取专业地图数据可以如此简单。本文将带你用最少的代码快速实现从数据获取到可视化的完整流程。1. 环境准备与基础概念在开始编码前我们需要明确几个关键概念。WMSWeb Map Service提供的是地图图片适合作为底图WFSWeb Feature Service则返回矢量地理要素数据适合做进一步分析。天地图作为国内权威的地理信息服务其接口符合OGC标准但需要注意坐标系统一使用CGCS2000。安装必要的Python库pip install requests geopandas matplotlib folium核心工具说明requests用于发送HTTP请求获取服务数据geopandas处理地理空间数据的利器matplotlib基础可视化工具folium生成交互式Leaflet地图注意所有示例基于福州天地图公开接口实际使用时请遵守其服务条款。2. 获取WMS底图并显示WMS服务最常用的操作是GetMap我们需要构造包含以下关键参数的URLlayers指定地图图层bbox地图范围左下角经度,左下角纬度,右上角经度,右上角纬度size图片尺寸width x heightsrs坐标参考系统EPSG代码import requests from PIL import Image from io import BytesIO # 福州天地图WMS服务地址 wms_url http://tdt.fuzhou.gov.cn/serviceaccess/WMS/DMDZ params { service: WMS, request: GetMap, version: 1.1.1, layers: 1, styles: default, bbox: 119.28,26.08,119.30,26.10, # 福州某区域 width: 800, height: 600, srs: EPSG:4490, # CGCS2000坐标系 format: image/png } response requests.get(wms_url, paramsparams) img Image.open(BytesIO(response.content)) img.show()常见问题处理图片显示为空白检查bbox参数是否超出服务范围坐标系统错误确保srs参数与服务支持的坐标系一致请求被拒绝可能需要添加token等认证参数3. 调用WFS服务获取矢量数据WFS的GetFeature操作可以获取地理要素的详细信息。以下示例获取福州市某区域的地名点数据import geopandas as gpd wfs_url http://tdt.fuzhou.gov.cn/serviceaccess/WFS/DMDZ_WFS-G params { service: WFS, request: GetFeature, version: 1.0.0, typeName: DMDZ:地名, bbox: 119.28,26.08,119.30,26.10, srsname: EPSG:4490, outputFormat: application/json } response requests.get(wfs_url, paramsparams) features gpd.read_file(response.text) print(features.head())关键参数说明typeName要素类型名称格式为命名空间:要素类型outputFormat指定返回格式GML/JSONmaxFeatures限制返回要素数量避免数据量过大4. 数据叠加可视化实战将WMS底图与WFS矢量数据结合展示是地理分析的常见需求。我们使用matplotlib实现import matplotlib.pyplot as plt from matplotlib.offsetbox import OffsetImage, AnnotationBbox # 创建画布 fig, ax plt.subplots(figsize(10, 8)) # 显示WMS底图 ax.imshow(img, extent[119.28, 119.30, 26.08, 26.10]) # 叠加WFS数据 features.plot(axax, colorred, markersize50, alpha0.7) # 添加标注 for idx, row in features.iterrows(): ax.annotate(row[名称], (row.geometry.x, row.geometry.y), xytext(3,3), textcoordsoffset points, bboxdict(boxstyleround, fcw)) plt.title(福州市地名点分布) plt.xlabel(经度) plt.ylabel(纬度) plt.grid() plt.show()进阶技巧使用folium创建交互式地图import folium # 创建底图 m folium.Map(location[26.09, 119.29], zoom_start15) # 添加WMS图层 folium.WmsTileLayer( urlwms_url, name福州底图, fmtimage/png, layers1, transparentTrue, attribution天地图 ).add_to(m) # 添加WFS要素 for idx, row in features.iterrows(): folium.Marker( location[row.geometry.y, row.geometry.x], popuprow[名称], iconfolium.Icon(colorred) ).add_to(m) # 添加图层控制 folium.LayerControl().add_to(m) m.save(fuzhou_map.html)5. 性能优化与实用技巧实际项目中我们还需要考虑以下问题缓存策略优化from requests_cache import install_cache install_cache(ogc_cache, expire_after3600) # 缓存1小时坐标系转换示例from pyproj import Transformer # CGCS2000转WGS84 transformer Transformer.from_crs(EPSG:4490, EPSG:4326) lon, lat transformer.transform(119.29, 26.09)参数自动编码处理from urllib.parse import quote bbox 119.28,26.08,119.30,26.10 encoded_bbox quote(bbox) # 处理特殊字符错误处理最佳实践try: response requests.get(wms_url, paramsparams, timeout10) response.raise_for_status() except requests.exceptions.RequestException as e: print(f请求失败: {e}) if hasattr(e, response) and e.response is not None: print(f错误详情: {e.response.text})最后分享一个实用技巧天地图服务通常有QPS限制在循环请求时建议添加延迟import time time.sleep(0.1) # 每次请求间隔0.1秒