FY-4B是中国风云四号静止轨道气象卫星系列的第二颗卫星 [1]，是在轨运行的风云四号A星（FY-4A）的姊妹星。在FY-4A的基础上，FY-4B进行了系统升级优化，各方面技术更为成熟，将成为中国第一颗第二代静止气象卫星，FY-4B_L1数据经纬度查找表数据说明-133E-V2 ### FY-4B-L1数据经纬度查找表数据说明 #### 一、FY-4B卫星简介 FY-4B是中国风云四号静止轨道气象卫星系列中的第二颗卫星，作为FY-4A的姊妹星，在FY-4A的基础上进行了系统升级优化，各方面技术更加成熟。FY-4B不仅继承了FY-4A的功能特性，还在多个方面进行了改进和增强，使其成为我国第一颗第二代静止气象卫星。 #### 二、FY-4B-L1经纬度查找表概述 ##### 1. 查找表功能 FY-4B L1数据经纬度查找表用于辅助FY-4B卫星的成像仪和快速成像仪L1数据处理，通过该表可以将图像的行列号转换为对应的地理坐标（经纬度），从而实现对特定地理区域的精确识别和分析。 ##### 2. 分辨率 查找表根据不同的分辨率生成，具体包括： - **250米** - **500米** - **1000米** - **2000米** - **4000米** 每种分辨率对应的查找表文件都是独立的，用户可以根据实际需求选择合适的分辨率进行数据处理。 ##### 3. 文件命名规则 查找表文件的命名规则为：“FY4B-_DISK_1330E_GEO_NOM_LUT_YYYYMMDDHHMMSS_RRRRR_V0001.raw”，其中： - **YYYYMMDDHHMMSS**：表示生成该查找表文件的时间戳。 - **RRRRR**：代表分辨率，例如“0250M”表示250米分辨率。 ##### 4. 文件格式 - **布局**：文件从北向南按行（从西向东）排列，每个数据单元存储对应的纬度和经度值。 - **单位**：经度和纬度均以度为单位。 - **类型**：数据类型为double，采用小端格式存储。 - **填充值**：未定义区域的经纬度值用999999.9999表示。 ##### 5. 标称网格大小 不同分辨率的查找表文件具有不同的行数和列数，以匹配其空间分辨率。具体对应关系如下： | 分辨率 | 行数 | 列数 | | --- | --- | --- | | 250米 | 43968 | 43968 | | 500米 | 21984 | 21984 | | 1000米 | 10992 | 10992 | | 2000米 | 5496 | 5496 | | 4000米 | 2748 | 2748 | ##### 6. 观测经纬度范围 不同分辨率下的观测范围有所不同，主要体现在可观察到的地理坐标范围上。以星下点133.0°E为中心，具体范围如下： | 分辨率 | 经度范围 | 纬度范围 | | --- | --- | --- | | 250米 | 51.711°~214.289°E (-145.711°W) | -81.295°~81.295° | | 500米 | 51.722°~214.278°E (-145.722°W) | -81.222°~81.222° | | 1000米 | 51.755°~214.245°E (-145.755°W) | -81.187°~81.187° | | 2000米 | 51.776°~214.224°E (-145.776°W) | -81.052°~81.052° | | 4000米 | 51.820°~214.180°E (-145.820°W) | -80.883°~80.883° | #### 三、读取示例 以4000米分辨率的查找表为例，以下是使用Matlab读取文件的示例代码： ```matlab clear all; clc; fid = fopen('FY4B-_DISK_1330E_GEO_NOM_LUT_20220317000000_4000M_V0001.raw','r'); data = fread(fid, [5496, 2748],'double','l'); fclose(fid); lat = data(1:2:end,:); %纬度 lon = data(2:2:end,:); %经度 ``` 通过上述代码，可以方便地从查找表文件中提取纬度和经度信息，并进一步应用于气象数据分析或地图可视化等场景中。这些信息对于理解FY-4B卫星的数据结构以及如何利用这些数据进行地理定位至关重要。

风云四（FY-4）数据读取代码-读取试例 包括HDF、HRIT以及NC三种格式数据

IDL对FY-4A辐射定标，写的不好，但可以用

风云4卫星数据格式说明

限制，只能处理FY-4的可见光与短波红外波段且是2km分辨率数据，这是一个处理FY-4数据的程序，主要进行两个处理，辐射定标与几何校正，定标为表观反射率。底层使用C++语言实现，主要用到HDF5与GDAL两个库，保存格式是GEOTIFF，坐标系为WGS-84地理坐标系。可以处理两个分辨率，2KM与4KM，由于没有进行分块处理，所以2KM数据占内存1.5G左右。