grads绘制分省地图问题

breeze_li · 发表于 2016-8-19 12:19:22

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

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

x

本帖最后由 breeze_li 于 2016-8-25 17:43 编辑

看了清风老师的grads作图详解，我试着绘制分省地图，前面的站点数据正常显示。在插值的时候用province-basemap绘制陕西地图的时候数据显示不出来，而用cnbasemap数据可以显示。应该是province-basemap中设置问题，请各位帮忙看看问题出在哪。下面是我的截图与报错： ~3W@6716NBMPH}EF~1UF47U.png

这是运行gs文件的报错，
sta.gs文件内容：
'reinit'
'open e:/project/text2/grid.ctl'
'open e:/project/text2/sta.ctl'
'set lon 105.29 111.15'
'set lat 31.42 39.35'
'set mpdset shannx'
'define  rgrid=oacres(g,temper.2)'
'set gxout shaded'
'd rgrid'
'set gxout contour'
'province-basemap rgrid'
'printim e:/project/text2/sta_grid.png white'
;
province-basemap路径：
say args
  if (args='')
  say usage1 ;say usage2
  return
  endif
  map = subwrd(args,1)
* map is provincename
  varname = subwrd(args,2)
  color = subwrd(args,3)
if (color = '') ; color = 0  ; endif
file = 'D:/software/opengrads/Contents/Cygwin/Versions/2.0.a9.oga.1/out/shannx_out.txt';say "read file name="file
'set mpdset'map
  'd 'varname
grid.ctl文件
DSET  E:/project/text2/grid.grd
TITLE Grid Data Sample
UNDEF -9.99E33
XDEF 3 LINEAR 105.29  2.5
YDEF 4 LINEAR 111.15  2.5
ZDEF  1 LEVELS 1000
TDEF 1 LINEAR 01may2011 1mo
VARS    1
g    0  99  Grid Data
ENDVARS
sta.ctl文件：
DSET E:/project/text2/sta.grd
DTYPE  station
STNMAP E:/project/text2/sta.map
UNDEF  -32766
TITLE  Temperature Data Sample
TDEF 1 linear 01may2011 1mo
VARS 1
temper 0  99  Temper Data
ENDVARS
站点数据截图：

插值后province-basemap（左）与cnbasemap（右）截图：

请各位帮忙看看

四叶草 · 发表于 2016-8-19 16:22:49

breeze_li 发表于 2016-8-19 15:32
两个参数？还缺个什么，我换成cnbasemap可以，province-basemap应该怎么改？

换成'province-basemap shannx rgrid' 看看。三个参数的说明在脚本里面有讲

breeze_li · 发表于 2016-8-19 12:20:32

下面是我的数据：
站点经度纬度小时温度
53947  109.25 35.5667 21.6
53845  110.1167 32.8167 26.1
57127  109.3167 36.8833 20.1
57016  110.0833 34.4833 18
53646  110.1833 36.0667 21.3
57036  108.2333 34.55 26.5

breeze_li · 发表于 2016-8-19 13:37:08

顶一下，求高手帮忙解答

四叶草 · 发表于 2016-8-19 15:00:06

楼主，估计是你的'province-basemap rgrid' 这个命令参数不对，至少需要2个参数

breeze_li · 发表于 2016-8-19 15:32:02

四叶草发表于 2016-8-19 15:00
楼主，估计是你的'province-basemap rgrid' 这个命令参数不对，至少需要2个参数

两个参数？还缺个什么，我换成cnbasemap可以，province-basemap应该怎么改？

breeze_li · 发表于 2016-8-19 16:30:04

四叶草发表于 2016-8-19 16:22
换成'province-basemap shannx rgrid' 看看。三个参数的说明在脚本里面有讲

可以了，非常感谢，绘制出来的图与cbarn可以自定义颜色吗，哪里有相关资料

四叶草 · 发表于 2016-8-19 17:09:43

breeze_li 发表于 2016-8-19 16:30
可以了，非常感谢，绘制出来的图与cbarn可以自定义颜色吗，哪里有相关资料

颜色可以自己定义的，你可以自己看看手册或是在家园看看

breeze_li · 发表于 2016-8-19 17:24:21

四叶草发表于 2016-8-19 17:09
颜色可以自己定义的，你可以自己看看手册或是在家园看看

中科院的那个？嗯，谢啦，你有没用过python跟grads结合绘图？

四叶草 · 发表于 2016-8-19 18:12:33

breeze_li 发表于 2016-8-19 17:24
中科院的那个？嗯，谢啦，你有没用过python跟grads结合绘图？

嗯嗯。Python我不会

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

[图形美化] grads绘制分省地图问题

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

回帖奖励 +2 金钱