用grads读取txt中经纬度画站点图

最新推荐文章于 2021-10-15 21:10:45 发布

转载最新推荐文章于 2021-10-15 21:10:45 发布 · 5.6k 阅读

GrADS 专栏收录该内容

30 篇文章

订阅专栏

本文介绍如何使用GRADS软件从TXT文件中读取站点的经纬度数据，并将其应用于绘制站点图的过程。具体步骤包括打开格点资料、读取TXT文件中的站点信息、解析经纬度数据并根据不同条件绘制标记。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

用grads读取txt中经纬度，用于画站点图
　　
　　txt中的站点信息格式如下：
　　
　　stn_num lon lat data1 data2 .... （每个变量用空格格开一行为一个记录）
　　
　　****************draw grid**************************************
　　
　　 'open ******.ctl' *打开一个格点资料用于画底图的值，譬如全球或者区域海温分布，有值就好
　　……
　　
　　……
　　
　　 'd var' *画出底图
　　
　　***********************draw station********************************
　　
　　cnt=1000
　　
　　i=1
　　
　　while(i　　
　　aa=read('westpacific.txt')
　　 aa1=sublin(aa,2)
　　 lat=subwrd(aa1,5)
　　 lon=subwrd(aa1,4)
　　
　　 data＝subwrd(aa1,6)
　　
　　'q w2xy 'lon' 'lat
　　
　　x = subwrd(result,3)
　　y = subwrd(result,6)
　　
　　if(data<=……) *判读控制mark的值，从而体现data的差别
　　'draw mark 8 'x' 'y' 0.01'
　　
　　endif
　　
　　if(data<=……)
　　
　　'draw mark 8 'x' 'y' 0.03'
　　
　　endif
　　
　　if(data<=……)
　　
　　'draw mark 8 'x' 'y' 0.06'
　　
　　endif
　　
　　………
　　i=i+1
　　endwhile
　　ret=close(‘westpacific.txt')
　　'cbarn 1 0 4.9 0.3'
　　'print'
　　'disable print'
　　'reinit'
　　