SNAP：哨兵3 OLCI level-1B产品处理

1、    Sentinel-3A  OLCI level-1数据介绍

        Sentinel-3A卫星于2016年2月16日成功发射， 其搭载的OLCI传感器具有从可见光到近红外的 21 个通道，空间分辨率为 300 m，波段范围从 400~1020 nm。OLCI level-1B产品包含的辅助数据集有：经纬度文件、观测几何文件（太阳天顶角、太阳方位角、卫星天顶角、卫星方位角）、大气参量信息文件（水平风速、海平面压力、臭氧总量、相对湿度、参考压力水平、大气温度廓线、水汽总量）。

         OL_1_EFR产品包含21个波段的TOA数据、质量标志和结合点网格数据。每个波段提供的数据为反射辐射度(Oa*_radiance)、太阳通量(solar_flux_band_*)、中心波长(lambda0_band_*)和全宽半最大值(FWHM_band_*)。这些数据的坐标集是每像素地理坐标(纬度、经度和海拔)。

        联络点网格数据，包括太阳传感器目标几何(太阳和观测方位角和天顶角)和来自欧洲中期天气预报中心(ECMWF)数据集的气象元数据:风速、850 hPa的湿度、海平面压力、大气水汽和臭氧含量以及大气温度廓线。联络点数据具有不同的坐标(tp_纬度，tp_经度)和跨径(经度)方向的不同分辨率(联络点为77列，地理坐标为4865列)，导致间隔为16 km。扎点(TP)是图像网格的下采样版本;因此，TP和像素之间的关联应通过它们的图像坐标(具有适当的子采样因子(4865−1)/(77−1)= 64(对于FR))来完成，而不是通过不同的地理坐标来完成。TP为卫星视距与参考椭球的交点，不存在观测视差校正。

        因此直接使用ENVI中建立GLT无法提取逐像元观测几何信息。Sentinel-3 影像数据的空间参考为地理查找表，不能直接得到影像的投影信息和经纬度坐标，因此使用 SNAP 软件对 Sentinel-3 遥感影像进行几何校正。

        本文主要对OLCI level-1B产品（FR）进行几何校正及观测几何提取。数据信息如下。

 OL_1_EFR产品（OL_1_ERR表示分辨率降低）及介绍如下表所示。资料来源哨兵数据官方指南：Level-1B Products Description - Sentinel-3 OLCI Technical Guide - Sentinel Online - Sentinel Online (esa.int)

 

2、SNAP几何校正

 snap中打开数据，xml文件

投影

 设置输出格式为ENVI格式并设置相关路径

设置输出坐标系，默认WGS84即可

 

 run，进行处理，时间可能稍长...

 处理完成后，输出文件中文件如下，且Img数据均为wgs84坐标系

 将需处理的波段加载进ENVI进行处理，波段辐射数据进行波段组合，角度数据进行转存即可。

波段组合如下。

 

 角度数据也可使用Layer Stacking工具进行转存。

3、几何校正结果

 波段数据（21 bands）几何校正及波段组合结果。（假彩色显示为17 6 3 波段组合 ）

 角度数据，以太阳天顶角为例。

参考资料：

[1]王峰,汪小钦,丁宇.利用Sentinel-3A OLCI可见光通道反演台湾岛气溶胶光学厚度

[2]《Google Earth Engine Sentinel-3 OLCI Level-1 Dataset Deviates from the Original Data: Causes and Consequences》

[3]Course: CMDS Essentials Pack for Ocean Applications | self-paced learning, Section: Downloading SNAP for Visualising data. (eumetsat.int)

[4]Sentinel-3 OLCI - Technical Guide - Sentinel Online - Sentinel Online (esa.int)

（小白，欢迎批评指正！）