读取风云4号卫星FullMask_Grid.raw文件

读取风云4号卫星FullMask_Grid.raw文件

1、下载FullMask_Grid.raw文件
下载网址：国家卫星气象中心（http://satellite.nsmc.org.cn/PortalSite/StaticContent/DocumentDownload.aspx?TypeID=3）

2、读取FullMask_Grid.raw文件

import numpy as np
import h5py

def fy4disk(rawfile,dim):
    """
    FY-4A数据行列号和经纬度查找表2km、4km，读取raw文件，存成hdf.
    :param rawfile: raw文件全路径
    :param dim: 行列数（2km:5496,4km:2748)
    :return:
    """
    sz = np.fromfile(rawfile, dtype=float, count=dim*dim*2)
    latlon = np.reshape(sz,(dim,dim,2))

    lat = latlon[:,:,0]
    lon = latlon[:,:,1]

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

    return lon, lat

if __name__ == '__main__':

    rawfile = r'D:\work\fire_point\FullMask_Grid_2000\FullMask_Grid_2000.raw'

    dim = 5496 # 2km
    #dim = 2748 # 4km

    # Two ways of reading an binary file
    # lon, lat = fy4raw(rawfile,dim)
    lon, lat = fy4disk(rawfile,dim)
    with h5py.File(r"D:\work\fire_point\FullMask_Grid_2000\FullMask_Grid_2000.hdf5",'w') as f:
        f["Lon"] = lon
        f["Lat"] = lat

参考网址：http://bbs.06climate.com/forum.php?mod=viewthread&tid=90074&highlight=raw