FY-4A图像数据的辐射定标和几何校正

1.1 FY-4A AGRI图像数据（辐射定标、几何校正）

（1）辐射定标

辐射定标就是把星上扫描辐射仪输出的计数值转换为要探测的物理量的绝对值。对于可见光一般转换为灰度和反照率，对于红外和水汽通道转换为灰度等级、辐射能量和绝对温度。 科学数据集中主要包括两部分，如图1。第一部分是图像数据层“NOMChannelX”，X表示14个通道编号，分别存储了14个通道的图像数据，即数字量化值（DN），没有物理意义。第7通道的有效数据范围是[0,65534]，其余通道的有效数据范围为[0，4095]，无效填充值为65535。第二部分是定标表“CALChannelX”，用于将数字量化值处理为反照率或亮温。数据类型为浮点型，第7通道的数据维度为65536×1，其余通道为4096×1。可见光/近红外波段有效值范围为[0,1.5]，DN值处理后反照率，红外波段的有效范围为[100,500]，DN值处理后为亮温。定标表中按DN值的取值范围，存放了各个有效DN值对应的反照率或亮温数据，形成一一对应的查找表，表的索引就是DN值。

辐射定标的步骤：

1. 读取多个波段的图像数据层的像元像素值；
2. 利用定标表对数据进行定标；
3. 将原始的数字信号转换为反射率等实际物理量；

 图1 风云四号多通道扫描成像辐射计L1 4KM数据列表

（2）几何校正

通过经纬度查找表校正法，将静止轨道卫星的标称投影转化为矢量数据的投影。也就是将图像数据的行列号转换成经纬度。主要的步骤如下：

a.   FY-4A数据行列号和经纬度查找表4km，文件格式为raw。需编辑头文件方可打开文件。

图2 4km行列号和经纬度查找表的头文件设置

b.  在图像与raw文件都需要做研究区域的roi。建立GTL，把画好的raw文件，band2为X，band1为Y。但会报错，错误原因是envi_build_glt_doit:An error has occured during processing Error:"Attempt to subscript DXS withj is out range." The result may be invalid.没有尝试成功。

图3  GLT文件投影设置

 c.  根据国家卫星气象中心给出的转换关系，进行数据行列号的经纬度相互转换，得到经纬度信息。【FY-4A卫星采用CGMS LRIT/HRIT全球规范定义的静止轨道标称投影，地理坐标基于WGS84参考椭球计算得到。】

图4 风云四号中国区图像数据