Sentinel-2L1C级数据是经正射校正和亚像元级几何精校正后的大气表观反射率(TOA)产品,在我们需要使用辐射率(Radiance)数据时需要一些转换,所以该数据的辐射定标不在是简单的乘以一个Gain,加上一个Offset那么简单。而是利用公式关系反算辐射率数据。其详细计算公式参见(ENVI官博)http://blog.sina.com.cn/s/blog_764b1e9d0102y7g4.html。下面就ENVI官博中各个参数如何获取做一点说明。
Sentinel-2L1C级数据根目录中包含一个名为MTD_MSIL1C.xml的文件,该文件为该景影像的元数据。
:大气表观反射率,这里其实是指各个波段的影像文件,打开MTD_MSIL1C.xml,在Granule_List节点下可以找到影像文件存放的相对路径。如GRANULE/L1C_T51RTQ_A015101_20180514T024123/IMG_DATA/T51RTQ_20180514T023551_B01。
d:日地距离,打开MTD_MSIL1C.xml文件,在Product_Info节点中有个名称为U的节点,该节点数值即为日地距离。如0.980958599408787。
:太阳辐照度,打开MTD_MSIL1C.xml文件,在Product_Info节点中有个名称为Solar_Irradiance_List的节点,即为每个波段对应的太阳辐照度值。如1884.69。
:太阳高度角,这个值不在MTD_MSIL1C.xml文件中,在相对路径如GRANULE/L1C_T51RTQ_A015101_20180514T024123下包含另一个名为MTD_TL的xml文件,打开这个xml文件,找到节点Mean_Sun_Angle,在此节点下包含两个角度,即ZENITH_ANGLE(天顶角),AZIMUTH_ANGLE(方位角)。我们知道太阳高度角与天顶角互余,因此天顶角的余角即为太阳高度角。如xml中读取ZENITH_ANGLE值为19.3760450927414,那么太阳高度角为70.6239549072586,单位为度。