python读取grads格式的ctl和dat格式的数据

樱桃小丸子

工作需要将grads可读的ctl加dat文件提取部分数据输出成txt格式，后续再进行处理 ，因为grads用的不是特别6，所以想用python实现，我知道python有grads的包，但网上相关数据少，没查到怎么用，也没查到是不是要linux才有那个包，想问问有没有大神知道怎么用  或者有没有其他python读取dat数据的方法，之前直接用读取二进制的方法打开数据不对

weinihou

不建议用 PyGrADS，这个是用Python调用GrADS，Python只是作为接口，没学过或用过GrADS的话相当于白瞎

我之前直接读取二进制是没错的，比如
data = numpy.fromfile(filename,dtype='f')
楼主可以再去试试，不要轻易 放 弃  治     疗嘛

不过这样对单变量的文件好一些，多变量的需要根据说明文件再后续加些调整

樱桃小丸子

weinihou 发表于 2017-12-20 20:21
https://stackoverflow.com/questions/14245094/how-to-read-part-of-binary-file-with-numpy

所以说 ...

非常感谢1111111

樱桃小丸子

weinihou 发表于 2017-12-14 14:05
fang qi zhi liao   竟然是会被屏蔽的

你好  我试着用你说的方法，这个文件2个g左右，python会报错   应该是文件太大 我试着把count设的小一点  是可以读出来的  但是全部就没办法读出来  我就没法下一步reshape了，还想问问你遇到过类似的情况吗

MeteoInfo

可以考虑用MeteoInfoLab，脚本语言是Jython，和Python语法一样。可以方便读取各类气象数据格式（包括GrADS），参考这里：http://www.meteothinker.com/docs ... ata/asciiwrite.html

樱桃小丸子

MeteoInfo 发表于 2017-12-14 11:49
可以考虑用MeteoInfoLab，脚本语言是Jython，和Python语法一样。可以方便读取各类气象数据格式（包括GrADS ...

因为后面转换后的txt数据还要再用python进行处理  所以需要整个流程全部用python

weinihou

fang qi zhi liao   竟然是会被屏蔽的

平流层的萝卜

4楼地板是正解，读取后对data.reshape(t维数，z维数，y维数，x维数)，就得到了相应的四维数组

weinihou

樱桃小丸子 发表于 2017-12-18 14:48
你好  我试着用你说的方法，这个文件2个g左右，python会报错   应该是文件太大 我试着把count设的小一点  ...

没搞过这么大的！
2个G的你都想一次性读进来啊，考虑到之后还要做一些操作和运算，就算都进来了也容易崩吧

我觉着lz你就一段一段搞好了…………

樱桃小丸子

weinihou 发表于 2017-12-19 10:54
没搞过这么大的！
2个G的你都想一次性读进来啊，考虑到之后还要做一些操作和运算，就算都进来了也容易崩 ...

那麻烦再请教一下  一段一段读的话data=numpy.fromfile(filename,dtype='f',count=？？？)  是不是可以一层一层的读  那比如说count第一段取前2000个，那第2000-4000个我用count怎么写呢？我在numpy的官方文档也没找到  或者有什么其他方法可以切成一段一段的吗

weinihou

樱桃小丸子 发表于 2017-12-20 10:46
那麻烦再请教一下  一段一段读的话data=numpy.fromfile(filename,dtype='f',count=？？？)  是不是可以一 ...

https://stackoverflow.com/questi ... ary-file-with-numpy

所以说要多用搜索引擎诶……