求路过高手解答画台风路径的问题~~

丹心侠骨

求路过高手解答画台风路径的问题~~，正在做毕业论文遇到的问题，出图不太对，范围一大低值中心很多就乱了，该怎么做呢？求大神指导啊？{:soso_e154:}附上gs文件

                               
登录/注册后可看大图

1.msl.gs
'reinit'
'!rm F:\\ec\\1323\\route\\track.txt'
'sdfopen F:\ec\1323\route\msl.nc'
'set grads off'
'set lon 110 124'
'set lat 20 35'
'set z 1'
i=20
while(i<=23)
'set t 'i''
'q time'
tt=subwrd(result,3)
year=substr(tt,9,4)
month=substr(tt,6,3)
day=substr(tt,4,2)
time=substr(tt,1,2)
name='F:\ec\1323\route\track.txt'
'define mpslv=minloc(min(msl,lat=20,lat=35),lon=110,lon=124)'
'd mpslv'
rcx=subwrd(result,4)
*x格点
'define mpslv=minloc(min(msl,lon=110,lon=124),lat=20,lat=35)'
'd mpslv'
rcy=subwrd(result,4)
*y格点
'set x 'rcx''
lon=subwrd(result,4)
'set y 'rcy''
lat=subwrd(result,4)
'd msl'
*minline=sublin(result,2)
*minvalue=subwrd(minline,4)
minvalue=subwrd(result,4)/100
'c'
record=year month day time ' ' lon ' ' lat ' ' minvalue ' '
rw=write(name,record,append)
i=i+1
endwhile

lat1=21
lat2=34
lon1=111
lat2=123

while(i<=32)
'set t 'i''
'q time'
tt=subwrd(result,3)
year=substr(tt,9,4)
month=substr(tt,6,3)
day=substr(tt,4,2)
time=substr(tt,1,2)
name='F:\ec\1323\route\track.txt'
'define mpslv=minloc(min(msl,lat-0.6='lat1',lat+0.6='lat1'),lon-0.6='lon1',lon+0.6='lon1')'
'd mpslv'
rcx=subwrd(result,4)
*x格点
'define mpslv=minloc(min(msl,lat-0.6='lat1',lat+0.6='lat1'),lon-0.6='lon1',lon+0.6='lon1')'
'd mpslv'
rcy=subwrd(result,4)
*y格点
'set x 'rcx''
lon=subwrd(result,4)
'set y 'rcy''
lat=subwrd(result,4)
lat1=lat
lon1=lon
'd msl'
*minline=sublin(result,2)
*minvalue=subwrd(minline,4)
minvalue=subwrd(result,4)/100
'c'
record=year month day time ' ' lon ' ' lat ' ' minvalue ' '
rw=write(name,record,append)
i=i+1
endwhile
;
2.ht.gs
'reinit' ;*第2部分--读出位置绘制台风路径
'sdfopen F:\ec\1323\route\msl.nc'
'set grads off'
'set mproj latlon'
'set mpdset hires'
'set lon 110 124'
'set lat 20 35'
'set z 1'
'set cmax -10000'
'set xlint 4' ;*定义横坐标的标记间隔
'set ylint 2' ;*定义纵坐标的标记间隔
'set xlopts 1 4 0.22' ;*横坐标 颜色 线宽 大小
'set ylopts 1 4 0.20' ;*纵坐标
'set clopts -1 -1 0.20' ;*等值线
'd msl'
*读取文件
read_file=read('F:\ec\1323\route\track.txt') ;*打开文件进行读取
while(1) ;*文件内容200行之内
xbefore=x0
ybefore=y0
read_file=read('F:\ec\1323\route\track.txt')
read_code=sublin(read_file,1) ;*读取文件打开代码
if(read_code>0);
break;
endif;
*读取正常为0，其他要么出错，要么结束
read_line=sublin(read_file,2) ;*读取文件行内容
say read_line
date0=subwrd(read_line,1) ;*取得第1部分
lon0=subwrd(read_line,2) ;*取得第2部分
lat0=subwrd(read_line,3) ;*取得第3部分
'q w2xy 'lon0' 'lat0'' ;*将经纬度坐标转换为GRADS坐标
x0=subwrd(result,3) ;*取得X坐标
y0=subwrd(result,6) ;*取得Y坐标
*say x0
*say y0
'draw title  bilis tracked by lsx'
'draw wxsym 40 'x0' 'y0' 0.3 2 6'
'draw line 'xbefore' 'ybefore' 'x0' 'y0''
*'draw string 'x0+0.1' 'y0' 'date0'';
endwhile
'gxprint F:\ec\1323\route\track.pdf white'
；

wk1021

丹心侠骨 发表于 2015-4-30 11:12
是菲特，好像登陆后，925hPa的被冷空气侵入，底层不明显了，谢谢啦

我说应该就是被其他系统影响的缘故，不然这种方法理论上是没错的，那现在你解决了问题没有？

丹心侠骨

自己顶一个

丹心侠骨

没有人帮忙吗

丹心侠骨

再顶一个，求解决

wk1021

是台风菲特的路径图么？

wk1021

不知道，是不是因为范围内还有其他的低压系统的干扰？或者你只画一小段，看看会不会乱？

丹心侠骨

wk1021 发表于 2015-4-29 17:12
不知道，是不是因为范围内还有其他的低压系统的干扰？或者你只画一小段，看看会不会乱？

是菲特，好像登陆后，925hPa的被冷空气侵入，底层不明显了，谢谢啦

丹心侠骨

wk1021 发表于 2015-4-30 13:34
我说应该就是被其他系统影响的缘故，不然这种方法理论上是没错的，那现在你解决了问题没有？

应该解决了，陆地上的路径不能画了

wk1021

丹心侠骨 发表于 2015-5-3 13:23
应该解决了，陆地上的路径不能画了

恩，加油~