- 积分
- 279
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2016-11-29
- 最后登录
- 1970-1-1
|
登录后查看更多精彩内容~
您需要 登录 才可以下载或查看,没有帐号?立即注册
x
本渣有时想加上一些自己矢量化的shp,比如青藏高原等等,但是总是报错,有高人知道自己矢量化的底图有什么要求吗?比如坐标之类的。
谢谢啦!
附上代码/自己矢量化的shp:
begin
f= addfile ("ua_FGOALS.nc","r")
u = f->ua
u700 = u(:,3,:,:)
uTim = dim_avg_n_Wrap(u700,0)
q = addfile ("va_FGOALS.nc","r")
v = q->va
v700 = v(:,3,:,:)
vTim = dim_avg_n_Wrap(v700,0)
wks = gsn_open_wks("pdf","stream_all")
res = True
res@gsnFrame = False
res@gsnMaximize = True
res@mpOutlineOn = True
res@mpDataBaseVersion = "MediumRes"
res@gsnAddCyclic = False
res@mpMinLonF = 80
res@mpMaxLonF = 120
res@mpMinLatF = 15
res@mpMaxLatF = 50
res@gsnMajorLatSpacing = 10
res@gsnMajorLonSpacing = 10
res@mpLandFillColor = "tan"
res@mpOceanFillColor = "LightBue"
res@mpInlandWaterFillColor = "LightBlue"
res@mpOceanFillColor = "LightBlue"
res@vcRefMagnitudeF = 1.
res@vcRefLengthF = 0.045
res@vcGlyphStyle = "CurlyVector"
res@vcMinDistanceF = 0.006
res@vcLineArrowColor = "blue"
plot = gsn_csm_vector_map(wks,uTim,vTim,res)
plres = True ; resource list for polylines
plres@gsLineColor = "black"
pgres = True
pgres@gsLineColor = "yellow"
psres = True
psres@gsLineColor = "black"
pfres = True
pfres@gsFillColor = "white"
id = gsn_add_shapefile_polylines(wks,plot,"bou1_4l.shp",plres)
id1 = gsn_add_shapefile_polylines(wks,plot,"desert.shp",pgres)
id2 = gsn_add_shapefile_polylines(wks,plot,"sandy_land.shp",psres)
id3 = gsn_add_shapefile_polygons(wks,plot,"Qinghai-Tibet_Plateau.shp",pfres)
draw(plot)
frame(wks)
end
|
|