【新人求助】python读取中国自动站与CMORPH降水产品融合的逐时降水量数据时出现问题

nsw704406795

本人气象小白，麻烦大神们帮看看！数据是中国气象数据网下载的中国自动站与CMORPH降水产品融合的逐时降水量网格数据集，grd格式。
1.网上找的python脚本
import numpy as np
file = 'E:/2019030117.grd'  # 文件路径
with open(file, 'rb') as f:              # 用with open 打开文件更智能
    a = np.fromfile(f, dtype=np.float32)   # 把二进制文件转为可读的数字
#a = np.reshape(a, (1400,440))           # rashape
#a = a[::-1]                            # 翻转矩阵
print(a)
2.中国气象数据网给的CTL文件
DSET ^SEVP_CLI_CHN_MERGE_FY2_PRE_HOUR_GRID_0.10-%y4%m2%d2%h2.grd
UNDEF -999.0
OPTIONS   little_endian  template
TITLE  China Hourly Merged Precipitation Analysis
xdef  700 linear  70.05  0.10
ydef  440 linear  15.05  0.10
ZDEF     1 LEVELS 1  
TDEF 9999
LINEAR 00Z01Aug2010 1hr
VARS 2                           
crain      1 00  CH01   combined analysis (mm/Hour)
gsamp      1 00  CH02   gauge numbers
ENDVARS
3.脚本运行后出现错误

数组大小616000，而不是CTL给的700*440，这是怎么回事呢？

4.我把数组大小改了一下运行成功后读出来的数组错误

616000=700*440*2，我把数组改成700*880，运行成功后读出来这个数组...
恳请大神们帮忙看看，给学生指点指点，感激不尽！

nsw704406795

数组全显示后...这是错误数据吗

Were_Mercy

请问楼主解决了么，同问！

tulalang

http://bbs.06climate.com/forum.php?mod=viewthread&tid=90192

Were_Mercy

是不是数据类型写错了呢？应该是<f4

微笑的橙子

因为你这里是两个变量哇~所以要乘2

微笑的橙子

Were_Mercy 发表于 2019-6-30 13:28
是不是数据类型写错了呢？应该是

float32应该也可以？

Were_Mercy

微笑的橙子 发表于 2019-11-2 16:41
float32应该也可以？

可以的，但是<f4表示使用的Little Endian的数据排列方式

Ay_1

请问楼主解决问题了吗