GrADS绘制站点分布图-升级版

mofangbao · 发表于 2013-10-30 22:59:19

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

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

x

经常看到有同学问绘制站点分布图的各种问题，我想很多人都有自己的绘制方法，比如用MeteoInfo就能非常方便的进行绘制。除了其他工具，GrADS也是可以绘制出质量不错的图形的，GrADS绘制站点图的思路基本都是一致的，最基础的就是直接d出来，但是这样不能设置点的样式，用的更多的是我很久之前转帖（点击查看）中那样类似的思路：获取站点经纬度，然后绘制，不过如果需要根据数据值的大小来分别设置绘制点的样式，甚至颜色大小呢，比如下面这样的图：

（由于数据是别人要发论文的，因此我在绘图时对结果做了处理，图像也很乱，只是说明绘制方法）

绘制思路：
利用站点格式的数据，设置输出到屏幕，将N个站点数据保存在变量中，然后对数据进行循环，逐行取数据的经纬度以及数据值，根据数据值判断应该绘制什么样的点，最后绘制图例。
注：本脚本参考了一位“未知师姐”的绘制程序，我对数据处理做了优化，原来的程序需要额外提供一个txt的文本，现在直接从站点文件取数据。

绘制脚本：

'reinit'
'open d:\test.ctl'
'set lon 70 140'
'set lat 15 55'
'set mpdset cnbase'
'set map 1 1 5.0'
'set grads off'
'set grid off'
'set xlopts 1 4 0.15'
'set ylopts 1 4 0.15'
'set xlint 10'
'set ylint 10'
'draw map'
'river'
'd maskout(r,r-30)'
'set gxout print'
'd b'
res = result
siz.1=0.12
siz.2=0.09
siz.3=0.06
siz.4=0.06
siz.5=0.09
siz.6=0.12
siz.7=0.15
i=2
while(1)
ln=sublin(res,i)
if(ln='');break;endif
i=i+1
lat1=subwrd(ln,3)
lon1=subwrd(ln,2)
data=sublin(res,i)
i=i+1
'q w2xy 'lon1' 'lat1''
x1=subwrd(result,3)
y1=subwrd(result,6)
data=subwrd(data,1)*1
'set line 6'
if(data<-5)
'draw mark 5 'x1' 'y1' 'siz.4''
endif
if(data>=-5 & data<-3)
'draw mark 2 'x1' 'y1' 'siz.5''
endif
'set line 1'
if(data>=-3 & data<0)
'draw mark 6 'x1' 'y1' 'siz.7''
endif
if(data>=0 & data<3)
'draw mark 1 'x1' 'y1' 'siz.6''
endif
'set line 9'
if(data>=3 & data<5)
'draw mark 6 'x1' 'y1' 'siz.3''
endif
if(data>=5 & data<8)
'draw mark 2 'x1' 'y1' 'siz.1''
endif
if(data>=8)
'draw mark 5 'x1' 'y1' 'siz.2''
endif
endwhile
'set line 1'
'set string 1 l 5 0'
'set strsiz 0.2 0.2'
'set font 1'
'draw string 1.0 7.0 (1)'
'set font 0'
'set string 1 l 5 0'
'set strsiz 0.12 0.15'
'draw string 1.3 3.0 <-5'
'draw string 1.3 2.7 -5~-3'
'draw string 1.3 2.4 -3~0'
'draw string 1.3 2.1 0~3'
'draw string 1.3 1.8 3~5'
'draw string 1.3 1.5 5~8'
'draw string 1.3 1.2 >8'
'set line 6'
'draw mark 5 1 3.0 'siz.1''
'draw mark 2 1 2.7 'siz.2''
'set line 1'
'draw mark 6 1 2.4 'siz.3''
'draw mark 1 1 2.1 'siz.4''
'set line 9'
'draw mark 6 1 1.8 'siz.5''
'draw mark 2 1 1.5 'siz.6''
'draw mark 5 1 1.2 'siz.7''
'southsea'
'enable print D:\test.gmf'
'print'
'disable print'
;

'reinit' 
'open d:\test.ctl' 
'set lon 70 140' 
'set lat 15 55' 
'set mpdset cnbase' 
'set map 1 1 5.0' 
'set grads off' 
'set grid off' 
'set xlopts 1 4 0.15' 
'set ylopts 1 4 0.15' 
'set xlint 10' 
'set ylint 10' 
'draw map' 
'river' 
'd maskout(r,r-30)' 
'set gxout print' 
'd b' 
res = result 
siz.1=0.12 
siz.2=0.09 
siz.3=0.06 
siz.4=0.06 
siz.5=0.09 
siz.6=0.12 
siz.7=0.15 
i=2 
while(1) 
    ln=sublin(res,i) 
    if(ln='');break;endif 
    i=i+1 
    lat1=subwrd(ln,3) 
    lon1=subwrd(ln,2) 
    data=sublin(res,i) 
    i=i+1 
'q w2xy 'lon1' 'lat1'' 
x1=subwrd(result,3) 
y1=subwrd(result,6) 
data=subwrd(data,1)*1 
'set line 6' 
if(data<-5) 
    'draw mark 5 'x1' 'y1' 'siz.4'' 
endif 
 
if(data>=-5 & data<-3) 
    'draw mark 2 'x1' 'y1' 'siz.5'' 
endif 
'set line 1' 
if(data>=-3 & data<0) 
    'draw mark 6 'x1' 'y1' 'siz.7'' 
endif 
 
if(data>=0 & data<3) 
    'draw mark 1 'x1' 'y1' 'siz.6'' 
endif 
'set line 9' 
if(data>=3 & data<5) 
    'draw mark 6 'x1' 'y1' 'siz.3'' 
endif 
 
if(data>=5 & data<8) 
    'draw mark 2 'x1' 'y1' 'siz.1'' 
endif 
 
if(data>=8) 
    'draw mark 5 'x1' 'y1' 'siz.2'' 
endif 
endwhile 
'set line 1' 
'set string 1 l 5 0' 
'set strsiz 0.2 0.2' 
'set font 1' 
'draw string 1.0 7.0 (1)' 
'set font 0' 
'set string 1 l 5 0' 
'set strsiz 0.12 0.15' 
 
'draw string 1.3 3.0  <-5' 
'draw string 1.3 2.7  -5~-3' 
'draw string 1.3 2.4  -3~0' 
'draw string 1.3 2.1  0~3' 
'draw string 1.3 1.8  3~5' 
'draw string 1.3 1.5  5~8' 
'draw string 1.3 1.2  >8' 
 
'set line 6' 
'draw mark 5 1 3.0 'siz.1'' 
'draw mark 2 1 2.7 'siz.2'' 
'set line 1' 
'draw mark 6 1 2.4 'siz.3'' 
 
'draw mark 1 1 2.1 'siz.4'' 
'set line 9' 
'draw mark 6 1 1.8 'siz.5'' 
'draw mark 2 1 1.5 'siz.6'' 
'draw mark 5 1 1.2 'siz.7'' 
 
'southsea' 
'enable print D:\test.gmf' 
'print' 
'disable print' 
;

绘制代码如下：

才仔 · 发表于 2013-10-30 23:08:14

抢楼

{:5_235:}

river · 发表于 2013-10-31 06:06:22

抢到地板也不错，支持！

lqouc · 发表于 2013-10-31 08:05:27

虽然这种图我一般都用meteoinfo了，但是还是喜欢可以自己编写脚本的，用着放心啊。

clare · 发表于 2013-10-31 08:17:15

3floor~

曾记否 · 发表于 2013-10-31 08:43:32

又可以学习新方法了，赞一个！

愉悦 · 发表于 2013-10-31 09:16:48

前几天还在求呢~

愉悦 · 发表于 2013-10-31 09:17:51

前几天还在求呢~

wk1984 · 发表于 2013-10-31 09:58:31

我觉得这种手动设置分段比较麻烦
换一个数据又得重新改很多

mofangbao · 发表于 2013-10-31 11:18:57

wk1984 发表于 2013-10-31 09:58

登录/注册后可看大图

我觉得这种手动设置分段比较麻烦
换一个数据又得重新改很多

分段当然是根据自己的需求来设置的，GrADS做这个本来就不是最方便的，只是介绍下思路而已，具体用什么实现每个人都有自己的选择~

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

[分享资料] GrADS绘制站点分布图-升级版

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

评分

相关帖子