读awx定量格点产品

pixi

参照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. 
   
7. file=r'E:\work\metplot\data\sat\cta\CTA_MLT_OTG_20170917_1030_FY2E.AWX'
   
8. f=open(file,'rb')
   
9. #一级头信息共40 byte
   
10. name,=unpack('12s',f.read(12))
    
11. name=name.decode()
    
12. 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))
    
13. format_instruct=format_instruct.decode()
    
14. 
    
15. #二级头信息共80 byte
    
16. sat_name,=unpack('8s',f.read(8))
    
17. para,para_byte,base_value,scalor,valid_time=unpack('5h',f.read(10))
    
18. 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))
    
19. 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))
    
20. f.read(24)
    
21. 
    
22. grid_num=x_grids*y_grids
    
23. format=str(grid_num)+'b'
    
24. data_field=unpack(format,f.read(grid_num))
    
25. data=np.array(data_field).reshape(x_grids,y_grids)
    
26. img=Image.fromarray(data)
    
27. img.show()

复制代码

但是出来的图错位了，左边的一块被移到了右边，没想明白，求助。
另外读到的文件名后面带\x00，怎样把\x00去掉？（b'FY2E\x00\x00\x00\x00'）

pixi

对比了下fy2g和fy2e的，数据区域范围不同，但都是偏离了12个经度，只好用笨办法把右边的区域拼接到左边了

18393810639

可以采用 np.split()对字符串进行拆解

dazhui121

请问楼主有没有处理静止卫星的呀