ORG+Cartopy完成shp文件筛选与显示

最新推荐文章于 2023-08-18 11:20:20 发布

黑眼圈的小白猫

最新推荐文章于 2023-08-18 11:20:20 发布

阅读量2.1k

点赞数 4

文章标签： python

本文链接：https://blog.csdn.net/qq_41850864/article/details/105168501

版权

在进行肝毕设时，使用Python编程绘制研究区概况图是常见任务。通常，利用basemap、cartopy.io.shapereader或shapefile库能轻松处理SHP文件。但遇到特殊SHP文件，会出现编码错误。通过尝试，发现GDAL库能有效解决此类问题，即使对复杂SHP文件也有强大读写能力。这是一个从错误中找到解决方案并记录下来的经验分享。

摘要由CSDN通过智能技术生成

肝毕设的过程总是充满艰辛，一个编程小白就此起步。
基本上来讲要做稍微宏观一点的调查研究，画一张好看的研究区概况图都是最基本的。作图中清晰不冗余的边界线和地理背景都是必不可少。这里用python编程来画张图，用得好确实节省时间，而且不装那些又大又沉重的软件感觉非常轻盈。
首先说一下手头的shp文件，有些shp文件确实非常友好，基本按照现有教程一遍就能跑出来，这样一般针对python的shp文件包无论是使用basemap，cartopy.io.shapereader还是shapefile都能解决。
比如cartpy的shp读取模块：

from cartopy.io.shapereader import Reader
import matplotlib.pyplot as plt

shp = Reader("filepath/file.shp")
region = [70,140,15,60]
proj = ccrs.PlateCarree()
plt.clf()
fig = plt.figure(figsize=(18,6),dpi=200)
ax = fig.add_subplot(1, 1, 1,projection = proj)
ax.set_extent(region,crs = ccrs.PlateCarree())
ax.add_geometries(shp.geometries(),crs = ccrs.PlateCarree(), edgecolor='m', linewidths=0.5,facecolor&#