cartopy如何读取shapefile的投影信息

sam_doggy

cartopy在读取shapefile的时候不能默认加载其中的.prj文件，也就是说无法自动获取其投影，我现在想到的做法是用geopandas读取后转换投影，不知道有没有方法不依赖geopandas库读取呢（应该暴力破解正则匹配prj文件内容也行，但是我不太了解.prj文件的格式，也不想重复造轮子）

1. import geopandas as gpd
   
2. shdf = gpd.read_file(filename)
   
3. shdf = shdf.to_crs("EPSG:4326")
   
4. shape_feature = cfeature.ShapelyFeature(shdf.geometry, ccrs.LambertConformal(), edgecolor='r', facecolor='none')
   
5. ax.add_feature(shape_feature)

复制代码

edwardli

1、解析WKT文本的工作可以从数据格式入手，自己写解析代码
2、或者使用GDAL、osgeo，基本是一句话实现，转为Proj.4子串，进而解析
3、但是，第一步的成本是高昂的，第二步的成本和使用Fiona或者geopandas进行坐标转换是一样的。
综上所述：cartopy不是全能的，也没有任何一个包是全能的；
专业的事情得找专业的包干，这也是Python所以强大的原因：开源、丰富；
考虑开发/研究成本，能解决问题的代码就是好代码。

sam_doggy

edwardli 发表于 2021-8-14 21:49
1、解析WKT文本的工作可以从数据格式入手，自己写解析代码
2、或者使用GDAL、osgeo，基本是一句话实现，转 ...

我看了cartopy的项目库，发现是在待开发的功能当中，我本以为这个是自带的，因为prj明显是shp文件的一部分，检查了半天才发现是shp投影的问题，，感觉浪费了不少时间，看了一下开发进度发现比想象中的成本高很多，这个停滞了几年了，确实应该把专业的功能交给专业库。