（更新）最最最准确的中国地图（附带南海，精确到市界）

尽头的尽头

subtropical 发表于 2018-7-24 18:29
哇哦 大神果然厉害 谢谢，还想请教一个问题，我把站点降水资料插值成格点画图，结果如图所示，我要怎么把 ...

mask一下，http://www.ncl.ucar.edu/Applications/shapefiles.shtml，参考21,22,23例子

subtropical

尽头的尽头 发表于 2018-7-25 08:36
mask一下，http://www.ncl.ucar.edu/Applications/shapefiles.shtml，参考21,22,23例子

谢谢大神，我看了这个例子，然后仿照这样来做，首先根据官网说的shapefile_mask_data[size=13.3333px] function in [size=13.3333px]shapefile_utils.ncl[size=13.3333px] 下载了shapefile_utils.ncl这个脚本，但是报错了说这个下载的脚本里面有句出问题了，我就看了下但是这个怎么改呢？麻烦您帮忙看下，谢谢您了



我的部分脚本：
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/shea_util.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/cnmap/cnmap.ncl"
load "./shapefile_utils.ncl"
begin

......
;---------------------------------使用cressman插值函数
rscan=(/1,0.5,0.25/)
r_final=obj_anal_ic_Wrap(lon, lat, rian, olon, olat, rscan, False)
data_mask = shapefile_mask_data(r_final,"D:/cygwin/app/ncl/lib/ncarg/nclscripts/cnmap/cnmap.shp",True)    ;---shapefile_mask_data

尽头的尽头

subtropical 发表于 2018-7-25 09:57
谢谢大神，我看了这个例子，然后仿照这样来做，首先根据官网说的shapefile_mask_data function in [/back ...

还需要一些属性设置的

subtropical

尽头的尽头 发表于 2018-7-25 15:56
还需要一些属性设置的

哦哦 那需要什么属性设置？ 可以教我一下吗？ 谢谢您了{:eb511:}

尽头的尽头

subtropical 发表于 2018-7-25 16:46
哦哦 那需要什么属性设置？ 可以教我一下吗？ 谢谢您了

例子里面有啊

subtropical

尽头的尽头 发表于 2018-7-25 22:28
例子里面有啊

您是说ncl官网第21那个例子吗？我好好看看学习，谢谢

subtropical

尽头的尽头 发表于 2018-7-25 22:28
例子里面有啊

大神您好，我参考例子21改了点，加载的家园的cnmap.shp.还有load "./shapefile_utils.ncl"
我的资料是站点资料sta_id=f_data1(:,0)
lon=f_data1(:,1)
lat=f_data1(:,2)
rian=f_data1(:,4)
这次没有那个357行的问题了，但是说shapefile_mask_data这个函数不存在，这个函数不是在加载的shapefile_utils.ncl里面吗。。。这该怎么办？麻烦您帮忙看下，谢谢您了

尽头的尽头

subtropical 发表于 2018-7-26 14:59
大神您好，我参考例子21改了点，加载的家园的cnmap.shp.还有load "./shapefile_utils.ncl"
我的资料是站 ...

是shapefile的变量名不对，不是这个函数不对，你检查一下

subtropical

尽头的尽头 发表于 2018-7-26 15:41
是shapefile的变量名不对，不是这个函数不对，你检查一下

哦哦 好的，我看了下shapefile_utils.ncl 里面有这个函数。。。。，难道是那个shapefile_name不对吗？是写cnmap.shp吗？

尽头的尽头

subtropical 发表于 2018-7-26 18:33
哦哦 好的，我看了下shapefile_utils.ncl 里面有这个函数。。。。，难道是那个shapefile_name不对吗？是 ...

查一下shapefile文件属性就知道了