问个小白问题，python如何读取.grd数据文件呢？

burf

比如data.grd这样的二进制格式数据文件，用python怎么读取呢？ 我尝试了一下open('file','rb')没有成功。。。读出来不知道是什么。。。
有知道的朋友不？求助一下，谢谢！

阿阿飞飞

不清楚你的文件格式，有没有头文件等信息，如果没有头文件信息，就是一个二进制数组的话就很好读。
给你看一个例子吧：

1. import numpy as np
   
2. 
   
3. #打开一个二进制文件
   
4. f=open('c:\\u.bin','r')
   
5. #使用numpy的fromfile来读取，注意dtype必须要对应，要不然读取出的数据不对
   
6. a=np.fromfile(f,dtype=np.float32)
   
7. f.close()
   
8. b=a.reshape(181,360,43)

复制代码

最后这个数组b就是我要的，一个三维数组。

jeffggg

阿阿飞飞 发表于 2013-11-18 21:45

                               
登录/注册后可看大图

不清楚你的文件格式，有没有头文件等信息，如果没有头文件信息，就是一个二进制数组的话就很好读。
给你看 ...

阿飞你好 很想学python 但是我看了一些教程感觉不怎么好学，对大妈里面的讲解不够 ，难以理解，你有好的建议吗？

贫道敬孔

多谢各位分享{:5_213:}

schliezer

你好，你这个问题解决了吗？我用楼下介绍的方法，读国家气象信息中心的格点化分析降水产品:
ydef 440 linear 15.05  0.1
xdef 700 linear 70.05  0.1
纬度起始：15.05,共440格点；0.1
经度起始：70.05,共700格点；0.1
但是读出来降水的值除了0就是XXe-38,XXe-40这种怪异值，不知是哪里出了问题

aifei

请问找到方法了，我试过用grads可是数据有问题

Alexander

二楼厉害啊{:5_275:}{:5_275:}{:5_275:}{:5_275:}{:5_275:}{:5_275:}

nsw704406795

schliezer 发表于 2018-9-21 17:19
你好，你这个问题解决了吗？我用楼下介绍的方法，读国家气象信息中心的格点化分析降水产品:
ydef 440 line ...

同问，大神有结果了吗！我读出来的更离谱呢[-999.-999.-999...0.0.0] 有点晕

Masterpiece

nsw704406795 发表于 2019-4-15 00:02
同问，大神有结果了吗！我读出来的更离谱呢[-999.-999.-999...0.0.0] 有点晕

-999是缺省值，国界以外的格点都是-999

陈三青

grib文件可以先用grads打开，转化成nc文件，然后再用python读取nc文件