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像素。