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

subtropical

hzwjy 发表于 2012-12-20 20:46
这个不错。
如果用NCL 6.1，用shapefile来添加中国国境线、省界什么的都是很简单的。

您好，这个shapefile可以分享下吗 谢谢 我把站点插值到格点上出现数值超出边界怎么改

subtropical

感谢大神的分享 我用这个cnmap画出了中国地图如下图，我想问怎么只显示中国的地图，旁边的那些怎么去掉？
我的脚本如下：
res=True
res@gsnDraw                        = False
res@gsnFrame                       = False
res@mpMinLatF         = min_lat
res@mpMaxLatF         = max_lat
res@mpMinLonF         = min_lon
res@mpMaxLonF         = max_lon
res@pmTickMarkDisplayMode="Always"

res@mpDataSetName               = "Earth..4"
res@mpDataBaseVersion           = "Ncarg4_1"
res@mpFillDrawOrder             = "PreDraw"
res@mpFillOn                    = True
res@mpLandFillColor             = "white"
res@gsnLeftString="station"
res@mpAreaMaskingOn=True
res@mpOutlineOn=True
map=gsn_csm_map(wks,res)
;;;;;;;;;绘制站点
mkres=True
mkres@gsMarkerColor="green"
mkres@gsMarkerSizeF=0.003
mkres@gsMarkerIndex=16
  dum=gsn_add_polymarker(wks,map,lon,lat,mkres)
;;;;添加中国地图和河流
cnres           = True
cnres@china     = True       ;draw china map or not
cnres@river     = True       ;draw changjiang&huanghe or not
cnres@province  = False     ;draw province boundary or not
cnres@nanhai    = False     ;draw nanhai or not
cnres@diqu      = False       ; draw diqujie or not
chinamap = add_china_map(wks,map,cnres)
draw(map)
frame(wks)
end

谢谢您了

尽头的尽头

subtropical 发表于 2018-7-24 16:53
感谢大神的分享 我用这个cnmap画出了中国地图如下图，我想问怎么只显示中国的地图，旁边的那些怎么去掉？
...

res@mpOutlineOn=True这个设置成False试试

subtropical

尽头的尽头 发表于 2018-7-24 17:48
res@mpOutlineOn=True这个设置成False试试

哇哦 大神果然厉害 谢谢，还想请教一个问题，我把站点降水资料插值成格点画图，结果如图所示，我要怎么把超出边界的去掉？谢谢您

尽头的尽头

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那个例子吗？我好好看看学习，谢谢