SNAP：哨兵2 L1C数据大气校正、植被参数反演及缺少波段问题

NaRi~

已于 2023-10-20 22:32:33 修改

阅读量3.4k

点赞数 7

分类专栏：图像预处理文章标签： sentinel 图像处理

于 2023-05-12 16:06:56 首次发布

本文链接：https://blog.csdn.net/qq_41358707/article/details/130642915

版权

图像预处理专栏收录该内容

5 篇文章 0 订阅

订阅专栏

1、介绍

S2波段设置。

使用SNAP Toolbox内的内置生物物理处理器，生物物理处理器使用八个反射带（B3、B4、B5、B6、B7、B8A、B11和B12），以及使用辐射传输模型（RTM）（即PROSAIL和神经网络算法）观测天顶、太阳天顶和相对方位角。

（资料来源：《Retrieval of Evapotranspiration from Sentinel-2:Comparison of Vegetation Indices, Semi-Empirical Models and SNAP Biophysical Processor Approach》LCC=CCC/LAI

PS：SNAP反演LAI及冠层叶绿素含量CCC

SNAP帮助文档算法描述：

Biophysical Processor Algorithm Specifications

对于Sentinel 2产品，已经实现了两种不同的神经网络架构，NNET20m和NNET10m。两者使用下表中所示的不同输入频带组合作为神经网络的输入。

NNET10m：仅使用10m波段，因此它产生10m分辨率的输出产品，但它只能计算LAI、FAPAR和FVC指数。（对应SNAP工具：Biophysical Processor S2_10m）

NNET20m：The Biophysical Processor S2 utilizes 8 bands at 20m resolution and it is capable of computing the LAI, Cab, CWC, FAPAR and FVC indexes. （对应SNAP工具：Biophysical Processor S2）

Biophysical Processor Description

The Biophysical Processor S2 utilizes 8 bands at 20m resolution and it is capable of computing the LAI, Cab, CWC, FAPAR and FVC indexes.

SNAP中Biophysical Processor反演的参数为冠层叶绿素含量，而非上述描述的叶片叶绿素含量，见以下帮助文档

2、大气校正

Sen2Cor插件下载地址： Sen2Cor – STEP

1、下载Sen2cor压缩包（本人下载的为Sen2Cor_v2.10版本），下载后解压到自定义文件夹下。

2、添加环境变量。此电脑-属性-高级系统设置-环境变量-PATH

3、此目录下运行cmd，执行L2A_Process --help命令，查看是否配置成功。

4、将下载解压好的L1C数据，进行大气校正（注意：L1C文件夹名称或路径不要太长，会报错）

执行命令：for /D %s in (S2A_MSIL1C路径*) do L2A_process %s

5、大气校正后为L2A级数据的文件夹（上为校正前，下为校正后）

3、重采样

3.1 方式一：下载安装Sen2Res，可以将20m\60m波段数据升尺度到10m的插件。

下载地址：Sen2Res – STEP

按照以下步骤对SNAP进行更新，重启

这种方法耗时较长，采样精度高

安装好插件按照以下步骤进行重采样

设置面板参数

保存类型为SANP格式

其他参数设置默认即可。

3.2 方式二：

1、打开L2A级数据，xml文件

2、进行重采样

重采样参数以B2为基准（10m），存储格式选择为BEAN-DIMAP。（若选择ENVI，则输出结果为单波段img，不便于进行LAI及Cab反演），重采样结果默认在SNAP中打开

按照自身需求修改分辨率，若只反演LAI等参数，不反演Cab，10m分辨率即可；若需要反演Cab，重采样至20m。

重采样结果如下：

3.3 方式三（最推荐，Sentinel2 重采样工具）

后续需要反演参数，保存为SNAP默认格式

4、LAI反演

4.1 10m分辨率反演

设置要反演的参数及数据源

4.2 20m分辨率反演

若选择S2工具，可反演LAI、Cab等参数

影像原图及LAI反演图

5、缺少波段B01、B10问题

用户文档如下：

2A级处理包括应用于大气顶部（TOA）1C级正射图像产品的场景分类和大气校正。2A级的主要输出是经大气校正的正射图像、地表反射率产品。
额外的输出是气溶胶光学厚度（AOT）图、水汽（WV）图和场景分类（SCL）图，以及60米分辨率下云和雪概率的质量指标（QI）。对于所有波段（10米、20米或60米），2A级输出图像产品被重新采样并以相等的空间分辨率生成。

10米：包含光谱带2、3、4、8、真彩色图像（TCI）以及从20米重新采样的AOT和WVP图。
20米：包含光谱波段1-7、波段8A、11和12、真彩色图像（TCI）、场景分类（SCL）图以及AOT和WVP图。频带B8被省略，因为B8A提供了更精确的光谱信息。
60米：包含重新采样到60米的20米产品的所有波段，以及带1和9、真彩色图像（TCI）、场景分类（SCL）图以及AOT和WVP图。卷云波段B10被省略，因为它不包含地表信息，若需要此波段，可在L1C中获取。

因此，大气校正后缺少的B10波段为60m卷云波段，不包含地表信息，不参与大气校正。

但20m文件夹中缺少B01波段，缺少的波段在60m文件夹中