grads绘制分省地图问题

lilispring

建议你好好看grads手册

breeze_li

四叶草 发表于 2016-8-25 12:54
回归书本，回归基础。楼主，建议你好好看grads手册。ctl文件的编写，手册上有讲的。看后不明白的大家再 ...

手册我看了，站点和格点的ctl文件我是根据数据编写的，然后用gs进行格点到站点的插值，绘制插值后的图片。这里因为我是要用python代开ctl文件，所以想问这里能不能把插值后数据编写一个ctl，然后在打开ctl文件，直接d变量出图？因为不是学这个的，而是有一个这样的需求，所以理解不对的地方，还请纠正。上面这个gmp文件是干什么的，跟map文件有什么区别？

四叶草

breeze_li 发表于 2016-8-25 15:20
手册我看了，站点和格点的ctl文件我是根据数据编写的，然后用gs进行格点到站点的插值，绘制插值后的图片 ...

要实现你所说的，按照格点数据的ctl编写一个就行啦

breeze_li

四叶草 发表于 2016-8-25 17:20
要实现你所说的，按照格点数据的ctl编写一个就行啦

格点数据的ctl不是只是一个网格么，没有数据，display出来就是一个空白网格。我这有一个例子model.ctl和model.gmp,在grads中代开model.ctl，然后d（ps/ta/ua/va...）就能出图
这个是model.ctl文件：
dset ^model.grb
title "Sample Model Data for lats4d Tutorial"
undef 1e+20
dtype grib
index ^model.gmp
xdef 72 linear 0.000000 5.000000
ydef 46 linear -90.000000 4.000000
zdef 7 levels
1000 850 700 500 300 200 100
tdef 5 linear 0Z1jan1987 1dy
vars 8
ps        0    1,  1,  0,  0 Surface pressure [hPa]
ua        7   33,100 Eastward wind [m/s]
va        7   34,100 Northward wind [m/s]
zg        7    7,100 Geopotential height [m]
ta        7   11,100 Air Temperature [K]
hus       7   51,100 Specific humidity [kg/kg]
ts        0   11,105,  2 Surface (2m) air temperature [K]
pr        0   59,  1,  0,  0 Total precipitation rate [kg/(m^2*s)]
endvars
model.gmp文件我这打不开，这个如果你有兴趣看的话就只能下载了，不好意思。

后来我尝试用python和grads的接口，在python中写了gs文件，可以成功绘图，不过想知道上面这种是怎么实现的。

breeze_li

四叶草 发表于 2016-8-25 17:20
要实现你所说的，按照格点数据的ctl编写一个就行啦

我更新了回帖奖励，附件的5个金币应该可以赚回去

四叶草

breeze_li 发表于 2016-8-25 17:44
我更新了回帖奖励，附件的5个金币应该可以赚回去

.gmp文件我就不清楚啦，平时绘图，我只用到.ctl，.map，.grd三个文件就行。你所说的python这个我也不懂

breeze_li

四叶草 发表于 2016-8-25 18:05
.gmp文件我就不清楚啦，平时绘图，我只用到.ctl，.map，.grd三个文件就行。你所说的python这个我也不懂

好吧，那这个例子中这种用map能实现么

四叶草

breeze_li 发表于 2016-8-25 18:18
好吧，那这个例子中这种用map能实现么

DSET  E:/project/text2/grid.grd  
TITLE Grid Data Sample
UNDEF -9.99E33  
XDEF 3 LINEAR 105.29  2.5  
YDEF 4 LINEAR 31.42  2.5  
ZDEF  1 LEVELS 1000   
TDEF   1 LINEAR 01may2011 1mo  
VARS      1
g     0  99  Grid Data
ENDVARS
你先把你例子的数据通过fwrite提取出来，作为单独的数据文件，与上面的ctl联合起来就能实现你说的结果。ctl里的数据名称应与你存放数据的名称一致。

wavelet

顶一下，求高手帮忙解答

breeze_li

wavelet 发表于 2016-8-25 18:51
顶一下，求高手帮忙解答

。。好人