请选择 进入手机版 | 继续访问电脑版
爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

新浪微博登陆

只需一步, 快速开始

搜索
查看: 4859|回复: 3

有关ncl叠图的问题

[复制链接]

新浪微博达人勋

发表于 2019-4-19 10:49:39 | 显示全部楼层 |阅读模式

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

您需要 登录 才可以下载或查看,没有帐号?立即注册 新浪微博登陆

x
请大家帮忙看看 程序运行没有错误 这个程序是为了在一张图上画出30个台风路径 然而我设置了先画后画的顺序 可是出来的图地图部分把路径盖住了 不知哪里有错 请大神指点
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"
wks=gsn_open_wks("png","track")
number=(/"01","02","03","04","05","06","07","08","09","10","11","12","13","14","15","16","17","18","19","20","21","22","23","24","25","26","27","28","29","30"/)
row=(/38,38,14,34,16,27,9,51,47,18,41,17,39,24,6,24,21,30,46,28,28,42,27,44,44,18,9,49,20,14/)
res=True
res@mpDataBaseVersion     = "MediumRes"     
;res@pmTickMarkDisplayMode = "Always"        ; better map tickmarks
res@gsnDraw=False
res@gsnFrame=False                       ;don't advance frame
res@mpMinLatF=0                     
res@mpMaxLatF=40                    
res@mpMinLonF=100
res@mpMaxLonF=160
res@mpFillDrawOrder="PreDraw"
res@mpLandFillColor="grey90"
res@tfPolyDrawOrder="PostDraw"
plot=gsn_csm_map(wks,res)               
resn=True
resn@gsLineColor=60
resp=True
do i=0,29
typhoon:=asciiread("/cygdrive/c/Users/Administrator/Desktop/1997 best track/"+number(i)+".txt",-1,"float")
a=row(i)
use:=onedtond(typhoon(8:),(/a,6/))
y:=use(:,2)*0.1
x:=use(:,3)*0.1
c=a-1
do b=0,c
x1:=use(b,3)*0.1
y1:=use(b,2)*0.1
gsn_polyline(wks,plot,x,y,resn)
  resp@gsMarkerColor=use(b,1)*10
resp@gsMarkerIndex=1
resp@gsMarkerSizeF=0.01
resp@tfPolyDrawOrder="PostDraw"
gsn_polymarker(wks,plot,x1,y1,resp)
end do
end do
ress=True
ress@txFontHeightF=0.015
ress@txPerimOn=True
ress@tfPolyDrawOrder="PostDraw"
xue=gsn_add_text(wks,plot,"",x(c),y(c),ress)
draw(plot)
frame(wks)



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

新浪微博达人勋

发表于 2019-4-19 13:04:59 | 显示全部楼层
res@mpFillDrawOrder="PreDraw"   这个去掉试试呢
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

新浪微博达人勋

发表于 2019-4-19 15:15:50 | 显示全部楼层
先draw(plot), 然后在gsn_polymarker(wks,plot,x1,y1,resp)
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

新浪微博达人勋

 楼主| 发表于 2019-4-19 16:31:36 | 显示全部楼层
还是一样呢
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

您需要登录后才可以回帖 登录 | 立即注册 新浪微博登陆

本版积分规则

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

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

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