FY-4A遥感影像（Disk/Regc区域）的几何校正
思路
-
选择区域经纬度范围，根据公式计算算出对应的行列号，并赋值
-
构建经纬度查找表，进行GLT几何校正
相较法1，法2速度快，且构建的GLT可以重复利用。故本节主要介绍构建经纬度查找表进行GLT几何校正的方法。
1. N_DISK全圆盘（4km）
有两种生成经纬度查找表的方法
一、直接读取官网提供的角度数据集raw，生成lon/lat数据集。
二、读取FY4A-AGRI_N_DISK_1047E_L1-GEO-MULT_NOM数据集，其中有圆盘的行列号数据集，无值区域为-1/有值区域为对应的行列号，然后根据公式由行列号计算出对应的经纬度，生成经纬度数据集。
步骤
-
生成经纬度查找表数据集（法2）
根据官网下载的经纬度查找表（4km），利用python读取raw文件（包含lon/lat数据集），生成longitude和latitude的hdf数据集。参考自（https://blog.csdn.net/tk20190411/article/details/115907130?utm_medium=distribute.pc_relevant.none-task-blog-2~default~BlogCommendFromMachineLearnPai2~default-3.control&depth_1-utm_source=distribute.pc_relevant.none-task-blog-2~default~BlogCommendFromMachineLearnPai2~default-3.control）
latitude
longitude
-
选定ROI区域，裁剪经纬度查找表和风云4A影像
subset data fromROIs
注意：裁剪区域必须在有值区域内，即不能包含圆盘以外的9999的无效值，否则ENVI在构建GLT的过程中会卡死。
-
Build GLT
注意：先选择lon，再选择lat，最后的Rotation选择0
-
Georeference from GLT
依次选择GLT文件、待校正数据，进行GLT几何校正。
结果展示
几何校正后的FY-4A裁剪区域
2. N_REGC中国区域（4km）
由于官网没有提供REGC的角度数据集，而只提供了REGC的GEO文件（FY4A-_AGRI--N_REGC_1047E_L1-GEO-_MULT_NOM_20210608071500_20210608071917_4000M_V0001），所以这里需要读取行列号计算经纬度，生成经纬度数据集。
步骤
以下步骤均相同。
-
选取ROI区域，裁剪经纬度查找表和风云4A影像
-
Build GLT
-
Georeference from GLT
结果展示
展望
-
可通过gdalwarp.exe构建VRT的方式完成批量几何校正，待研究。
参考李民录老师的csdn博客
-
裁剪ROI的目的是为了确保构建GLT的部分是有值区域（数组没有无效值（0/nan）和错值（lon>360，lat not in [-90, 90]），否则ENVI在创建GLT的过程中会报错。
其实也可不裁剪，对经度/纬度数据集进行插值，插补圆盘无效值区域，将圆盘边角全部填满。目前暂未成功。