NCL中如何添加自定义掩膜？

雨里1998

这是百度了某个网友给出方法的代码，但报错了，而且我也在官网上确实没找到shapefile_mask_data函数

1. begin
   
2.     cata="d:/paper/climate dynamics/"
   
3.     shpfile="tibetan/tibetan.shp"
   
4.     file1="ERA5_dataset/surface_pressure1.nc"
   
5.     f=addfile(cata+file1,"r")
   
6.     pres=short2flt(f->sp)
   
7.     opt=True
   
8.     opt@return_mask=True
   
9.     p_mask=shapefile_mask_data(pres,cata+shpfile,opt)
   
10.     p_land_mask=where(p_mask .eq. 1,pres,pres@_FillValue)
    
11.     copy_VarMeta(pres,p_land_mask)
    
12. end

复制代码

雨里1998

已经解决了，要导入官网提供的shapefile_utils.ncl脚本