CMORPH数据处理

枳忧9822

【新人求助】
请问各位大神，如何利用读取空间分辨率8km,时间分辨率为30min的CMORPH数据呢？找到了它的ctl文件如下：
DSET ../CRT/8km-30min/%y4/%y4%m2/CMORPH_V1.0_CRT_8km-30min_%y4%m2%d2%h2
OPTIONS template little_endian
UNDEF  -999.0
TITLE  Precipitation estimates
XDEF 4948 LINEAR   0.036378335 0.072756669
YDEF 1649 LINEAR -59.963614    0.072771377
ZDEF   01 LEVELS 1
TDEF 999999 LINEAR  00z01jan1998 30mn
VARS 1
cmorph   1  99  cmorph [ mm/hr ]
ENDVARS

微笑的橙子

with open(CMORPH_path, 'rb') as f:
    cmorph = np.fromfile(f, dtype=np.float64)
cmorph = np.reshape(cmorph, (4948, 1649))
cmorph = cmorph[::-1]
cmorph[cmorph == -999] = np.nan

枳忧9822

微笑的橙子 发表于 2019-11-2 16:46
with open(CMORPH_path, 'rb') as f:
    cmorph = np.fromfile(f, dtype=np.float64)
cmorph = np.resh ...

非常感谢您的回复，已经解决啦

Bonnenuitzz

请问能共享一下怎么解决的吗

gaoyanbo22

微笑的橙子 发表于 2019-11-2 16:46
with open(CMORPH_path, 'rb') as f:
    cmorph = np.fromfile(f, dtype=np.float64)
cmorph = np.resh ...

这个reshape(4949,1649)还是（1649，4948）啊？第一个数不应该是行数吗？