用NCL画中国地图

天行健 · 发表于 2011-8-25 20:37:05

登录后查看更多精彩内容~

您需要登录才可以下载或查看，没有帐号？立即注册

x

欢迎使用NCL的同仁加入NCL讨论QQ群：123189536

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/wrf/WRFUserARW.ncl"
;************************************************
begin
;f=addfile("hgt.2002.nc","r")
a=asciiread ("barmes.txt", (/25,29,3/) , "float")
b=asciiread ("hgt.txt", (/25,29/) , "float")
lat=fspan(10,70,25)
lon=fspan(60,130,29)
lat@units="degrees_north"
lon@units="degrees_east"
a!0="lat"
a!1="lon"
a&lat=lat
a&lon=lon

wks = gsn_open_wks("X11","barnes")
res=True
;;;;;;;;;;;限定地图区域绘图;;;;;;地图和数据要分别设置，互相独立;;;;;;;;;;;;;
res@gsnAddCyclic       = False
    res@mpLimitMode    = "LatLon"
    res@mpMinLatF       = 10
    res@mpMaxLatF       = 70
    res@mpMinLonF       = 60
res@mpMaxLonF       = 130
;;;;;;;;;;;;;只显示中国区域的图;;;;;;;;;;;;;;;;;;;;;;
res@mpDataBaseVersion="Ncarg4_1" ;中等分辨率
res@mpDataSetName="Earth..4" ;第4版地图，有中国边界数据
res@mpOutlineOn          = True
res@mpOutlineSpecifiers=(/"China:states","Taiwan"/) ;中国边界有问题，缺藏南、台湾
res@mpOutlineBoundarySets ="NoBoundaries"
;;;;;;;;;;;填充颜色;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
res@cnFillDrawOrder  = "PreDraw"    ;设置画图顺序，先填充颜色，后把海洋覆盖掉
res@mpAreaMaskingOn = True    ;使能填充覆盖
   res@mpMaskAreaSpecifiers = (/"China:states","Taiwan"/)
res@mpOceanFillColor = 1      ;用白色填充海洋0是colormap的索引值
res@mpInlandWaterFillColor = 0

res@cnFillOn    = True
    res@cnLinesOn    = False    ;等值线不显示
    res@cnLineLabelsOn = False

plot = gsn_csm_contour_map(wks,a(:,:,0),res)
end

中国各省市边界名：

993 3 2 7 Land : Eurasia : China
994 4 2 6 Land : Eurasia : China : Anhui
995 4 2 7 Land : Eurasia : China : Beijing Shi
996 4 2 6 Land : Eurasia : China : Chongqing Shi
997 4 2 7 Land : Eurasia : China : Fujian
998 4 2 7 Land : Eurasia : China : Fujian . Haitan Dao
999 4 2 7 Land : Eurasia : China : Fujian . Nanri Dao
1000 4 2 7 Land : Eurasia : China : Fujian . Xiamen Dao (Amoy)
1001 4 2 5 Land : Eurasia : China : Gansu
1002 4 2 4 Land : Eurasia : China : Guandong
1003 4 2 4 Land : Eurasia : China : Guandong . Donghai Dao
1004 4 2 4 Land : Eurasia : China : Guandong . Gaolan Dao
1005 4 2 4 Land : Eurasia : China : Guandong . Hailing Dao
1006 4 2 4 Land : Eurasia : China : Guandong . Nansan Dao
1007 4 2 4 Land : Eurasia : China : Guandong . Naozhou Dao
1008 4 2 4 Land : Eurasia : China : Guandong . Shangchuan Dao
1009 4 2 4 Land : Eurasia : China : Guandong . Xiachuan Dao
1010 4 2 4 Land : Eurasia : China : Guandong . Xinliao Dao
1011 4 2 7 Land : Eurasia : China : Guangxi
1012 4 2 7 Land : Eurasia : China : Guangxi . Weizhou Dao
1013 4 2 7 Land : Eurasia : China : Guangxi . Xieyang Dao
1014 4 2 5 Land : Eurasia : China : Guizhou
1015 4 2 5 Land : Eurasia : China : Hainan
1016 4 2 6 Land : Eurasia : China : Hebei
1017 4 2 5 Land : Eurasia : China : Heilongjiang
1018 4 2 8 Land : Eurasia : China : Henan
1019 4 2 7 Land : Eurasia : China : Hong Kong
1020 4 2 7 Land : Eurasia : China : Hong Kong . Dayushan
1021 4 2 7 Land : Eurasia : China : Hubei
1022 4 2 8 Land : Eurasia : China : Hunan
1023 4 2 5 Land : Eurasia : China : Jiangsu
1024 4 2 5 Land : Eurasia : China : Jiangxi
1025 4 2 6 Land : Eurasia : China : Jilin
1026 4 2 5 Land : Eurasia : China : Liaoning
1027 4 2 5 Land : Eurasia : China : Liaoning . Changxing Dao
1028 4 2 5 Land : Eurasia : China : Liaoning . Dachangshan Dao
1029 4 2 5 Land : Eurasia : China : Liaoning . Guanglu Dao
1030 4 2 5 Land : Eurasia : China : Liaoning . Haiyang Dao
1031 4 2 5 Land : Eurasia : China : Liaoning . Juhua Dao
1032 4 2 5 Land : Eurasia : China : Liaoning . Shicheng Dao
1033 4 2 5 Land : Eurasia : China : Liaoning . Xizhong Dao
1034 4 2 5 Land : Eurasia : China : Liaoning . Zhangzi Dao
1035 4 2 7 Land : Eurasia : China : Nei Mongol
1036 4 2 6 Land : Eurasia : China : Ningxia Huizu
1037 4 2 4 Land : Eurasia : China : Qinghai
1038 4 2 4 Land : Eurasia : China : Shaanxi
1039 4 2 7 Land : Eurasia : China : Shandong
1040 4 2 7 Land : Eurasia : China : Shandong . Daqin Dao
1041 4 2 7 Land : Eurasia : China : Shandong . Lingshan Dao
1042 4 2 7 Land : Eurasia : China : Shandong . Nanchangshan Dao
1043 4 2 7 Land : Eurasia : China : Shandong . Tuoji Dao
1044 4 2 7 Land : Eurasia : China : Shanghai Shi
1045 4 2 7 Land : Eurasia : China : Shanghai Shi . Chongming Dao
1046 4 2 5 Land : Eurasia : China : Shanxi
1047 4 2 7 Land : Eurasia : China : Sichuan
1048 4 2 4 Land : Eurasia : China : Tianjin Shi
1049 4 2 7 Land : Eurasia : China : Xinjiang Uygur
1050 4 2 7 Land : Eurasia : China : Xinjiang Uygur . (Sinkiang)
1051 4 2 8 Land : Eurasia : China : Xizang
1052 4 2 8 Land : Eurasia : China : Xizang . (Tibet)
1053 4 2 4 Land : Eurasia : China : Yunnan
1054 4 2 4 Land : Eurasia : China : Zhejiang
1055 4 2 4 Land : Eurasia : China : Zhejiang . Dai Shan
1056 4 2 4 Land : Eurasia : China : Zhejiang . Daqu Shan
1057 4 2 4 Land : Eurasia : China : Zhejiang . Liuheng Dao
1058 4 2 4 Land : Eurasia : China : Zhejiang . Sijiao Shan
1059 4 2 4 Land : Eurasia : China : Zhejiang . Taohua Dao
1060 4 2 4 Land : Eurasia : China : Zhejiang . Zhoushan Dao
1061 4 2 4 Land : Eurasia : China : Zhejiang . Zhujia Jian

2221 3 2 8 Land : Eurasia : Taiwan
2222 3 2 8 Land : Eurasia : Taiwan . Jinmen Dao (Quemoy)
2223 3 2 8 Land : Eurasia : Taiwan . Lan Yu
2224 3 2 8 Land : Eurasia : Taiwan . Liuch'iu Yu
2225 3 2 8 Land : Eurasia : Taiwan . Lu Tao
2226 3 2 8 Land : Eurasia : Taiwan . Mienhua Yu
2227 3 2 8 Land : Eurasia : Taiwan . P'engchia Yu
2228 3 2 8 Land : Eurasia : Taiwan . P'enghu Ch'untao (Pescadores)

善人/jw · 发表于 2011-8-25 21:17:47

抢个沙发~~

athlon8048 · 发表于 2011-8-25 21:34:47

沙发没了，坐板凳！

支持一下楼主！

artmunich · 发表于 2011-8-25 21:45:34

虽然暂时不懂，但貌似很厉害~~

易小凯 · 发表于 2011-8-26 09:22:58

load "$NCARG_ROOT/lib/ncarg/nclscripts/wrf/WRFUserARW.ncl"
这个是NCL里面自带的？

天行健 · 发表于 2011-8-27 09:42:17

易小凯发表于 2011-8-26 09:22

登录/注册后可看大图

load "$NCARG_ROOT/lib/ncarg/nclscripts/wrf/WRFUserARW.ncl"
这个是NCL里面自带的？

en
en
en
en

ningwentao · 发表于 2011-9-2 09:59:34

天哥，我是小宁，问一下，既然ncl的图有问题，用国家地理信息中心的图怎么画？

archtim · 发表于 2012-2-27 16:40:24

请问这几个文本文件内容如何呢？如barmes.txt

吴龙刚 · 发表于 2012-5-26 10:54:12

支持天哥！

刘昊炎 · 发表于 2012-5-30 08:23:58

同问barmes.txt和hgt.txt

		自动登录	找回密码
密码			立即注册

[经验总结] 用NCL画中国地图

登录后查看更多精彩内容~

评分