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

sam_doggy

                               
登录/注册后可看大图

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我有些区域插值报错，遂用最新版本计算的方法
有空我写个总结

hzwjy

1. 
   
2. def fy4disk(rawfile,dim):
   
3.     sz = np.fromfile(rawfile, dtype=float, count=dim*dim*2)
   
4.     latlon = np.reshape(sz,(dim,dim,2))
   
5. 
   
6.     lat = latlon[:,:,0]
   
7.     lon = latlon[:,:,1]
   
8. 
   
9.     lat[lat > 100] = -9999.
   
10.     lon[lon < 0  ] = lon[lon < 0  ] + 360.
    
11.     lon[lon > 361] = -9999.
    
12. 
    
13.     return lon, lat
    
14. 
    
15. if __name__ == '__main__':
    
16. 
    
17.     rawfile = r'FullMask_Grid_4000.raw'
    
18. 
    
19.     # dim = 5496 # 2km
    
20.     dim = 2748 # 4km
    
21. 
    
22.     # Two ways of reading an binary file
    
23.     # lon, lat = fy4raw(rawfile,dim)
    
24.     lon, lat = fy4disk(rawfile,dim)
    

复制代码

sam_doggy

hzwjy 发表于 2019-4-9 10:30

谢谢！这个注释很详细

slddbs

hzwjy 发表于 2019-4-9 10:30

用您的方法将raw数据读取，但为2748*2748，下载的FY4A数据为nc格式，查看变量为2748*1108这个怎么解，感谢专家啦

slddbs

用您的方法将raw数据读取，但为2748*2748，下载的FY4A数据为nc格式，查看变量为2748*1108这个怎么解，感谢专家啦

sam_doggy

slddbs 发表于 2019-7-10 06:03
用您的方法将raw数据读取，但为2748*2748，下载的FY4A数据为nc格式，查看变量为2748*1108这个怎么解，感谢 ...

2748*2748 是4000M分辨率的某一时次的网格的行列数,应该不会变.不太懂1108是什么.

一阵妖风

slddbs 发表于 2019-7-10 06:03
用您的方法将raw数据读取，但为2748*2748，下载的FY4A数据为nc格式，查看变量为2748*1108这个怎么解，感谢 ...

应该是全圆盘和区域数据的区别

一阵妖风

hzwjy 发表于 2019-4-9 10:30

你好，请问有matlab读取的代码吗？谢谢了！

一阵妖风

hzwjy 发表于 2019-4-9 10:30

能不能问一下，他的说明文件里说的是经度在前，但为什么读出来的两个数组第一个是纬度第二个是经度呢？有点困惑，谢谢了！

slddbs

一阵妖风 发表于 2019-9-30 10:41
能不能问一下，他的说明文件里说的是经度在前，但为什么读出来的两个数组第一个是纬度第二个是经度呢？有 ...

你输出数值，看数值范围基本就明白了