【已解决】python画青藏高原阴影报错

又又又加零

参考家园里这位老师的帖子青藏高原阴影教程，在中国地图基础上进行叠加青藏高原shp文件增加阴影时发生了报错：

1. fea_country = cfeat.ShapelyFeature(shapereader.Reader(country_shp).geometries(),proj,edgecolor='k', facecolor='ivory')
   
2. ax.add_feature(fea_country, linewidth=0.5, alpha=0.9)
   
3. 
   
4. fea_pro = cfeat.ShapelyFeature(pro_shp.geometries(), proj,edgecolor='darkolivegreen', facecolor='white')
   
5. ax.add_feature(fea_pro, linewidth=0.3, alpha=1)

复制代码

报错显示如下：

----> 4 fea_pro = cfeat.ShapelyFeature(pro_shp.geometries(), proj,edgecolor='darkolivegreen', facecolor='white')AttributeError: 'str' object has no attribute 'geometries'

看报错应该是shp文件的问题，但是用该shp文件是可以在ncl里进行添加青藏高原阴影的，实在不知道该怎么改，有好心人可以解答吗？不胜感激！使用的shp文件是该贴的文件青藏高原shp

【解决方法】按照上贴的链接教程的代码写即可，也不知为何我按照回复里的老师写的赋属性就不行，如果后续有老师能够解惑就好了~

千与千叶

这个错误是因为 pro_shp 是一个字符串对象(string)，并不是一个具有 geometries 属性的对象。请检查代码，确保 pro_shp 是一个具有 geometries 属性的对象，比如 shapely.geometry.Polygon 或 MultiPolygon 对象。如果 pro_shp 是一个字符串对象，则需要将其转换成具有 geometries 属性的对象。常见的方法是使用 geopandas 库中的 read_file() 函数读取 shapefile 文件，并指定 geometry 列名，例如：

import geopandas as gpd

pro_shp = gpd.read_file('path/to/pro.shp')
pro_shp = pro_shp.set_geometry('geometry')
这将返回一个具有 geometries 属性的 GeoDataFrame 对象，可以直接传递给 cfeat.ShapelyFeature() 函数。

myj0219

谢谢分享！！！！！

又又又加零

千与千叶 发表于 2023-4-6 00:20
这个错误是因为 pro_shp 是一个字符串对象(string)，并不是一个具有 geometries 属性的对象。请检查代码， ...

非常感谢您的解答，原来是这样，但尝试用您的方法修改后，依旧报错了。报错依旧是这条，是否是因为我使用的shp文件属性不完整呢？

----> 4 fea_pro = cfeat.ShapelyFeature(pro_shp.geometries(), proj,edgecolor='darkolivegreen', facecolor='white')
-> 5902 return object.__getattribute__(self, name)AttributeError: 'GeoDataFrame' object has no attribute 'geometries'

查看shp文件，是下图的形式。

清宵逝水

又又又加零 发表于 2023-4-6 09:55
非常感谢您的解答，原来是这样，但尝试用您的方法修改后，依旧报错了。报错依旧是这条，是否是因为我使用 ...

如果已经进行了pro_shp = pro_shp.set_geometry('geometry')的操作，那直接cfeat.ShapelyFeature(pro_shp) 就好了

又又又加零

清宵逝水 发表于 2023-4-7 15:21
如果已经进行了pro_shp = pro_shp.set_geometry('geometry')的操作，那直接cfeat.ShapelyFeature(pro_shp ...

已经解决了，感谢您的帮助