Matlab输出二进制再用Grads画图的错误问题

酷睿314 · 发表于 2014-4-25 19:44:40

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

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

x

由于要把多年SST方差分布图画出来，我通过MATLAB把ＮＣ文件读取，并计算，得到一个矩阵（行列由于原本由NC读取，故没有进行变换），而后用它生成dat并编写ctl，可是用grads能够读取，读进去后却只画一点的分布图，而不是整个格点。。。试了很多别的方法都是这样。另外，我在MATLAB得到的方差分布图中将缺测值0设定问1E20。。求教。。

river · 发表于 2014-4-25 19:48:59

直接用grads计算方差，然后画图就好了

酷睿314 · 发表于 2014-4-25 19:50:35

river 发表于 2014-4-25 19:48
直接用grads计算方差，然后画图就好了

怎么用grads计算呢？能大概告诉一下gs怎么写么？

river · 发表于 2014-4-25 20:32:32

酷睿314 发表于 2014-4-25 19:50
怎么用grads计算呢？能大概告诉一下gs怎么写么？

没什么特别的，计算方差的公式书上有，把变量代进去就可以了，类似于fortran编程

Aires · 发表于 2014-4-25 21:40:20

楼主这样说太模糊了吧，把具体的脚本粘上来。

沙颖凯 · 发表于 2014-4-25 22:13:14

用MATLAB输出GrADS使用的数据一个是要注意维数的存放，还有一个是fwrite要选float选项

kongfeng0824 · 发表于 2014-4-26 13:20:21

就用matlab画吧

酷睿314 · 发表于 2014-5-3 22:53:00

最后直接用grads写了个脚本算方差搞定了。。

泰达巴神 · 发表于 2015-4-23 21:25:31

沙颖凯发表于 2014-4-25 22:13
用MATLAB输出GrADS使用的数据一个是要注意维数的存放，还有一个是fwrite要选float选项

你好，我在matlab中计算了一个矩阵a(29,17,8)，分别代表x,y,z方向，存储的语句为：
fid = fopen('D:\f\a.dat','w');
count = fwrite(fid,a,'float32');
fclose(fid);
接着用grads打开，当z=1时是可以画图的，但set z 2....也就是高层的时候，出现了错误
Entire grid contents are set to missing data
Grid limits of file: X = 1 29 Y = 1 17 Z = 1 8 T = 1 1 E = 1 1
Grid limits of request: X = 1 29 Y = 1 17 Z = 2 2 T = 1 1 E = 1 1
Warning issued for variable = v
Z的层次那里出问题了，请问该怎么办？

血嗜子木 · 发表于 2015-10-14 10:33:03

我也很知道这门技术

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

[脚本编辑] Matlab输出二进制再用Grads画图的错误问题

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