ncl中shapefile_mask_data提取中国区域内数据问题

Akihabara

想利用ncl提供的shapefile_mask_data提取中国区域内格点数据，但是print出的数据几乎都是缺测，mask_in也都是缺测，原数据无缺测，不知道哪里出了问题  TT 谢谢各位

tidiff是插值后的格点数据（无缺测）
Coordinates:
lat: [17..57]
lon: [72..136]
脚本如下：
shp_filename="$NCARG_ROOT/lib/ncarg/nclscripts/cnmap/cnmap.shp"
mask_in = shapefile_mask_data(tdiff(:,:),shp_filename,True)
mask_io = where(ismissing(mask_in), 0, 1)
erai_mask  = tdiff                                                        
erai_mask  = mask (erai_mask, mask_io, 1)  
copy_VarCoords(tdiff ,erai_mask)                 
print( erai_mask )

Akihabara

shp文件的问题