爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

新浪微博登陆

只需一步, 快速开始

搜索
楼主: 尽头的尽头

[作图] (更新)最最最准确的中国地图(附带南海,精确到市界)

    [复制链接]

新浪微博达人勋

 楼主| 发表于 2018-7-25 08:36:11 | 显示全部楼层
subtropical 发表于 2018-7-24 18:29
哇哦 大神果然厉害 谢谢,还想请教一个问题,我把站点降水资料插值成格点画图,结果如图所示,我要怎么把 ...

mask一下,http://www.ncl.ucar.edu/Applications/shapefiles.shtml,参考21,22,23例子
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2018-7-25 09:57:05 | 显示全部楼层
本帖最后由 subtropical 于 2018-7-25 09:59 编辑

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

2.png
我的部分脚本:
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


密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2018-7-25 15:56:03 | 显示全部楼层
subtropical 发表于 2018-7-25 09:57
谢谢大神,我看了这个例子,然后仿照这样来做,首先根据官网说的shapefile_mask_data function in [/back ...

还需要一些属性设置的
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2018-7-25 16:46:54 | 显示全部楼层
尽头的尽头 发表于 2018-7-25 15:56
还需要一些属性设置的

哦哦 那需要什么属性设置? 可以教我一下吗? 谢谢您了{:eb511:}
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2018-7-25 22:28:15 | 显示全部楼层
subtropical 发表于 2018-7-25 16:46
哦哦 那需要什么属性设置? 可以教我一下吗? 谢谢您了

例子里面有啊
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2018-7-26 11:18:23 | 显示全部楼层

您是说ncl官网第21那个例子吗?我好好看看学习,谢谢
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2018-7-26 14:59:35 | 显示全部楼层
本帖最后由 subtropical 于 2018-7-26 15:16 编辑

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

新浪微博达人勋

 楼主| 发表于 2018-7-26 15:41:26 | 显示全部楼层
subtropical 发表于 2018-7-26 14:59
大神您好,我参考例子21改了点,加载的家园的cnmap.shp.还有load "./shapefile_utils.ncl"
我的资料是站 ...

是shapefile的变量名不对,不是这个函数不对,你检查一下
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2018-7-26 18:33:46 | 显示全部楼层
尽头的尽头 发表于 2018-7-26 15:41
是shapefile的变量名不对,不是这个函数不对,你检查一下

哦哦 好的,我看了下shapefile_utils.ncl 里面有这个函数。。。。,难道是那个shapefile_name不对吗?是写cnmap.shp吗? 6.png 5.png
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2018-7-26 20:51:09 | 显示全部楼层
subtropical 发表于 2018-7-26 18:33
哦哦 好的,我看了下shapefile_utils.ncl 里面有这个函数。。。。,难道是那个shapefile_name不对吗?是 ...

查一下shapefile文件属性就知道了
密码修改失败请联系微信:mofangbao
您需要登录后才可以回帖 登录 | 立即注册 新浪微博登陆

本版积分规则

Copyright ©2011-2014 bbs.06climate.com All Rights Reserved.  Powered by Discuz! (京ICP-10201084)

本站信息均由会员发表,不代表气象家园立场,禁止在本站发表与国家法律相抵触言论

快速回复 返回顶部 返回列表