python用salem.roi失败

bzz

今天画图发现这样问题：
正常应该是这样：

                               
登录/注册后可看大图

换了一个shp以后变成这样：

                               
登录/注册后可看大图

网上似乎没有好的解决办法，
特来问问大家，非常感谢！
下面是shp和代码：


http://bbs.06climate.com/forum.php?mod=attachment&aid=MTExNjYwfGViMDk2YThhZjE0NDJhZDM1MGU0ZmQwY2Y2NmIwOWE2fDE3NTczOTc5NDc%3D&request=yes&_f=.shp，http://bbs.06climate.com/forum.php?mod=attachment&aid=MTExNjg5fGU3Y2RiZWIxOGVmOWIwNjQ3OTAxZDQ2MDMwZDdlNDkxfDE3NTczOTc5NDc%3D&request=yes&_f=.dbf，http://bbs.06climate.com/forum.php?mod=attachment&aid=MTExNjkwfDc0Y2I0NWU4OGY3ZTNkZTQwNTVhYmNiZTNmZjc4YWE4fDE3NTczOTc5NDc%3D&request=yes&_f=.prj，http://bbs.06climate.com/forum.php?mod=attachment&aid=MTExNjkxfDMyYmY4YTlhM2FjZDlmY2NhYTlkZTc1YzMxZjBlZmZlfDE3NTczOTc5NDc%3D&request=yes&_f=.sbn，http://bbs.06climate.com/forum.php?mod=attachment&aid=MTExNjkyfDY0MjgyNTU2ZWRjZGI0NGQ2ODI1NjZiY2UyMWVkODQ3fDE3NTczOTc5NDc%3D&request=yes&_f=.sbx。

1. def cut_shp(data):
   
2.     import geopandas as gpd    ;   import salem
   
3.     shp = "DBATP\\DBATP_Line.shp"
   
4.     shpfile=gpd.read_file(shp)
   
5.     data=data.salem.roi(shape=shpfile)
   
6.     return data
   
7. 
   
8. def pr_data_select():
   
9.     data_path='Pre_1961_2022_monthly.nc'
   
10.     data_all = xr.open_dataset(data_path).pre
    
11.     data_all = cut_shp(data_all)
    
12.     return data_all

复制代码

灭火器

你只传一个shp，不传shx和dbf其他人是读不了的。
不过根据文件名带line，猜一个是因为文件是Polyline类型，需要尝试在geopandas里手动转换成Polygon

Jinnee

楼上说得对，你这个可能是线文件，得线转面

bzz

灭火器 发表于 2024-3-9 10:30
你只传一个shp，不传shx和dbf其他人是读不了的。
不过根据文件名带line，猜一个是因为文件是Polyline类型 ...

非常感谢指点，问题已经解决了