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

mofangbao

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

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

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

绘制脚本：

1. 'reinit'
   
2. 'open d:\test.ctl'
   
3. 'set lon 70 140'
   
4. 'set lat 15 55'
   
5. 'set mpdset cnbase'
   
6. 'set map 1 1 5.0'
   
7. 'set grads off'
   
8. 'set grid off'
   
9. 'set xlopts 1 4 0.15'
   
10. 'set ylopts 1 4 0.15'
    
11. 'set xlint 10'
    
12. 'set ylint 10'
    
13. 'draw map'
    
14. 'river'
    
15. 'd maskout(r,r-30)'
    
16. 'set gxout print'
    
17. 'd b'
    
18. res = result
    
19. siz.1=0.12
    
20. siz.2=0.09
    
21. siz.3=0.06
    
22. siz.4=0.06
    
23. siz.5=0.09
    
24. siz.6=0.12
    
25. siz.7=0.15
    
26. i=2
    
27. while(1)
    
28.     ln=sublin(res,i)
    
29.     if(ln='');break;endif
    
30.     i=i+1
    
31.     lat1=subwrd(ln,3)
    
32.     lon1=subwrd(ln,2)
    
33.     data=sublin(res,i)
    
34.     i=i+1
    
35. 'q w2xy 'lon1' 'lat1''
    
36. x1=subwrd(result,3)
    
37. y1=subwrd(result,6)
    
38. data=subwrd(data,1)*1
    
39. 'set line 6'
    
40. if(data<-5)
    
41.     'draw mark 5 'x1' 'y1' 'siz.4''
    
42. endif
    
43. 
    
44. if(data>=-5 & data<-3)
    
45.     'draw mark 2 'x1' 'y1' 'siz.5''
    
46. endif
    
47. 'set line 1'
    
48. if(data>=-3 & data<0)
    
49.     'draw mark 6 'x1' 'y1' 'siz.7''
    
50. endif
    
51. 
    
52. if(data>=0 & data<3)
    
53.     'draw mark 1 'x1' 'y1' 'siz.6''
    
54. endif
    
55. 'set line 9'
    
56. if(data>=3 & data<5)
    
57.     'draw mark 6 'x1' 'y1' 'siz.3''
    
58. endif
    
59. 
    
60. if(data>=5 & data<8)
    
61.     'draw mark 2 'x1' 'y1' 'siz.1''
    
62. endif
    
63. 
    
64. if(data>=8)
    
65.     'draw mark 5 'x1' 'y1' 'siz.2''
    
66. endif
    
67. endwhile
    
68. 'set line 1'
    
69. 'set string 1 l 5 0'
    
70. 'set strsiz 0.2 0.2'
    
71. 'set font 1'
    
72. 'draw string 1.0 7.0 (1)'
    
73. 'set font 0'
    
74. 'set string 1 l 5 0'
    
75. 'set strsiz 0.12 0.15'
    
76. 
    
77. 'draw string 1.3 3.0  <-5'
    
78. 'draw string 1.3 2.7  -5~-3'
    
79. 'draw string 1.3 2.4  -3~0'
    
80. 'draw string 1.3 2.1  0~3'
    
81. 'draw string 1.3 1.8  3~5'
    
82. 'draw string 1.3 1.5  5~8'
    
83. 'draw string 1.3 1.2  >8'
    
84. 
    
85. 'set line 6'
    
86. 'draw mark 5 1 3.0 'siz.1''
    
87. 'draw mark 2 1 2.7 'siz.2''
    
88. 'set line 1'
    
89. 'draw mark 6 1 2.4 'siz.3''
    
90. 
    
91. 'draw mark 1 1 2.1 'siz.4''
    
92. 'set line 9'
    
93. 'draw mark 6 1 1.8 'siz.5''
    
94. 'draw mark 2 1 1.5 'siz.6''
    
95. 'draw mark 5 1 1.2 'siz.7''
    
96. 
    
97. 'southsea'
    
98. 'enable print D:\test.gmf'
    
99. 'print'
    
100. 'disable print'
     
101. ;
     

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

绘制代码如下：

才仔

抢楼{:5_235:}

river

抢到地板也不错，支持！

lqouc

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

clare

3floor~

曾记否

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

愉悦

前几天还在求呢~

愉悦

前几天还在求呢~

wk1984

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

mofangbao

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

                               
登录/注册后可看大图

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

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