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

气象家园

 找回密码
 立即注册

新浪微博登陆

只需一步, 快速开始

QQ登录

只需一步,快速开始

搜索
查看: 208|回复: 2

[求助] 风云卫星FullMask_Grid.raw文件如何读取呢[已解决]

[复制链接] |关注本帖

新浪微博达人勋

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

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

x
本帖最后由 sam_doggy 于 2019-4-9 08:27 编辑


                               
登录/注册后可看大图

FY-4A数据行列号和经纬度查找表4km
这个文件,我用
  1. rawfile = np.fromfile(r'E:\climate bbs\FullMask_Grid_4000.raw')
复制代码
读取后发现无法使用reshape使得其变为2748*2748的矩阵,因为他的维数为(15103008L,)

https://blog.csdn.net/modabao/article/details/83834580 参考上面博客得知答案,注意是他的github的早期版本,后期改为计算了
然而他的早期版本插值我还有些问题,读进去raw我有些区域插值报错,遂用最新版本计算的方法
有空我写个总结
密码修改失败请联系qq:937062711

新浪微博达人勋

发表于 2019-4-9 10:30:05 | 显示全部楼层 |取消关注该作者的回复

  1. def fy4disk(rawfile,dim):
  2.     sz = np.fromfile(rawfile, dtype=float, count=dim*dim*2)
  3.     latlon = np.reshape(sz,(dim,dim,2))

  4.     lat = latlon[:,:,0]
  5.     lon = latlon[:,:,1]

  6.     lat[lat > 100] = -9999.
  7.     lon[lon < 0  ] = lon[lon < 0  ] + 360.
  8.     lon[lon > 361] = -9999.

  9.     return lon, lat

  10. if __name__ == '__main__':

  11.     rawfile = r'FullMask_Grid_4000.raw'

  12.     # dim = 5496 # 2km
  13.     dim = 2748 # 4km

  14.     # Two ways of reading an binary file
  15.     # lon, lat = fy4raw(rawfile,dim)
  16.     lon, lat = fy4disk(rawfile,dim)
复制代码
密码修改失败请联系qq:937062711
回复 支持 反对

使用道具 举报

新浪微博达人勋

 楼主| 发表于 2019-4-9 16:41:12 | 显示全部楼层 |取消关注该作者的回复

谢谢!这个注释很详细
密码修改失败请联系qq:937062711
回复 支持 反对

使用道具 举报

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

本版积分规则

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

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

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