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

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

新浪微博登陆

只需一步, 快速开始

搜索
查看: 12327|回复: 3

[求助] 读awx定量格点产品

[复制链接]

新浪微博达人勋

发表于 2017-10-20 10:50:13 | 显示全部楼层 |阅读模式

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

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

x
本帖最后由 pixi 于 2017-10-20 12:13 编辑

参照Awx的格式说明,自己写了段代码处理云量产品
  1. from struct import *
  2. import matplotlib.pyplot as plt
  3. from PIL import Image
  4. import numpy as np
  5. from math import *   

  6. file=r'E:\work\metplot\data\sat\cta\CTA_MLT_OTG_20170917_1030_FY2E.AWX'
  7. f=open(file,'rb')
  8. #一级头信息共40 byte
  9. name,=unpack('12s',f.read(12))
  10. name=name.decode()
  11. byte_seq,first_head_len,second_head_len,fulfill_len,rec_len,head_record_num,product_data_record_num,product_type,compression,format_instruct,quality=unpack('9h8sh',f.read(28))
  12. format_instruct=format_instruct.decode()

  13. #二级头信息共80 byte
  14. sat_name,=unpack('8s',f.read(8))
  15. para,para_byte,base_value,scalor,valid_time=unpack('5h',f.read(10))
  16. begin_year,begin_mon,begin_day,begin_hour,begin_min,end_year,end_mon,end_day,end_hour,end_min=unpack('10h',f.read(20))
  17. upper_left_lat,upper_left_lon,lower_right_lat,lower_right_lon,grid_unit,x_int,y_int,x_grids,y_grids=unpack('9h',f.read(18))
  18. f.read(24)

  19. grid_num=x_grids*y_grids
  20. format=str(grid_num)+'b'
  21. data_field=unpack(format,f.read(grid_num))
  22. data=np.array(data_field).reshape(x_grids,y_grids)
  23. img=Image.fromarray(data)
  24. img.show()
复制代码
但是出来的图错位了,左边的一块被移到了右边,没想明白,求助。
另外读到的文件名后面带\x00,怎样把\x00去掉?(b'FY2E\x00\x00\x00\x00')

自己画的图

自己画的图

meteoinfo出的图

meteoinfo出的图
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2017-10-21 09:15:27 | 显示全部楼层
对比了下fy2g和fy2e的,数据区域范围不同,但都是偏离了12个经度,只好用笨办法把右边的区域拼接到左边了
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2018-11-26 15:20:26 | 显示全部楼层
可以采用 np.split()对字符串进行拆解
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2020-9-29 10:28:28 | 显示全部楼层
请问楼主有没有处理静止卫星的呀
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

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

本版积分规则

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

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

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