风云4A/4B卫星行列号和经纬度查找表文件下载及读取方式

下载地址

https://satellite.nsmc.org.cn/PortalSite/StaticContent/DocumentDownload.aspx?TypeID=15

如图：

点击绿色小图标下载。

下载完了后是压缩包，解压缩后进入文件夹获取raw格式的文件

注意！！

注意133°E对应的卫星文件名里的1330E，
注意105°E对应的卫星文件名里的1050E.

读取文件

解压缩后拿到数据文件如图

时间部分信息不重要。解析代码如下：

#用PYTHON读取RAW文件
import numpy as np
rawfile=r"H:\temp_estimate\FY4A-readanddraw\FullMask_Grid_4000.raw"
dim=2748
data=np.fromfile(rawfile,dtype=float,count=dim*dim*2)
latlon=np.reshape(data,(dim,dim,2))
lat=latlon[:,:,0]
lon=latlon[:,:,1]

打印lat和lon,结果如下：

性能对比

官网上还给了一种方式是公式计算。对比一下两者的效率。
公式计算封装成函数，用时0.8s左右；查表用时0.08s左右。查表速度快十倍。