请选择 进入手机版 | 继续访问电脑版
爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

新浪微博登陆

只需一步, 快速开始

QQ登录

只需一步,快速开始

搜索
查看: 67|回复: 1

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

[复制链接] |关注本帖

新浪微博达人勋

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

您需要 登录 才可以下载或查看,没有帐号?立即注册 新浪微博登陆

x
本帖最后由 nsw704406795 于 2019-4-15 10:15 编辑

本人气象小白,麻烦大神们帮看看!数据是中国气象数据网下载的中国自动站与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.脚本运行后出现错误

图1.jpg
数组大小616000,而不是CTL给的700*440,这是怎么回事呢?
4.我把数组大小改了一下运行成功后读出来的数组错误
图2.jpg
616000=700*440*2,我把数组改成700*880,运行成功后读出来这个数组...
恳请大神们帮忙看看,给学生指点指点,感激不尽!


密码修改失败请联系qq:937062711

新浪微博达人勋

 楼主| 发表于 4 天前 | 显示全部楼层 |取消关注该作者的回复
3.png
数组全显示后...这是错误数据吗
密码修改失败请联系qq:937062711
回复 支持 反对

使用道具 举报

您需要登录后才可以回帖 登录 | 立即注册 新浪微博登陆

本版积分规则

Copyright ©2011-2014 bbs.06climate.com All Rights Reserved.  Powered by Discuz! (京ICP-10201084)

本站信息均由会员发表,不代表气象家园立场,禁止在本站发表与国家法律相抵触言论

快速回复 返回顶部 返回列表