哨兵2号预处理

1 下载（https://scihub.copernicus.eu/dhus/#/home）

 

   登录——选区域——选日期——选类型

PS：如果显示图像offline（时间太久远就会被暂时下架），需要申请上线（添加购物车，点下载后会出现提示栏，绿色说明已提交申请，红色是提交失败，需要过会儿再试试），大概两三个小时就可以下载了。

2 Sen2Cor里进行辐射定标和大气校正（L1C转L2A）

  打开cmd，输入Sen2Cor文件所在的绝对路径，输入下面代码进行批处理，加粗部分是L1C文件所在地址，其他不变。

for /D %s in (F:\L1C\S2A_MSIL1C*) do L2A_process %s

（这样说明处理完成）

（处理好的数据会和原始数据在一个文件夹内）

3 利用SNAP输出ENVI模式的波段

  在SNAP中输入哨兵的L2A文件，依此选择File-Import-Optical Sensors-Sentinel-2-S2-MSI L2A，然后选择MTD_MSIL2A.xml。

   然后选择Raster-Data Conversion-Band Select，波段一般选择B2（Blue）、B3（Green）、B4（Red）和B8（NIR），这四个波段分辨率均为10m。同时，将输出格式选择为ENVI，确定输出位置，点击Run。

 

4 利用ENVI进行波段融合（Layer Stacking）

  波段选择输出的波段是分开的，加载到ENVI里是单波段灰色，需要利用Layer Stacking模块进行波段融合，得到一个彩色图像。打开ENVI，在Toolbox选择Raster Management，找到Layer Stacking（波段融合），将想要融合的波段加载进来，点击Reorder Files将波段进行排序，选择投影坐标和输出路径，选择重采样方式，最后确定。

 （上图是四景波段融合后的影像)

 5 利用ENVI进行拼接

  在Toolbox选择Mosaicking-Seamless Mosaic，打开后点击绿色加号，把影像都加载进来，将Date Ignore Value值都改为零，选择一景影像为基准。

   点击Color Correction-Histogram Matching-Entirs Scene

   点击Seamlines/Feathering-Edge Feathering

   点击Seamlines-Auto Generate Seamlines选择接边线

   打开后选择ENVI格式，输出地址，背景值设为0，重采样选择双线性内插法（按需而定），然后点结束，开始镶嵌。

  效果不太理想，改进：接边线自己点线（先保存，再调用），右下的图像为基准（reference），羽化范围设成10像素。