全球10米遥感数据免费下载及预处理(附软件下载)
看标题估计大家都猜出来了是Sentinel-2(哨兵2),没错,就是Sentinel-2。
之前我们推过关于高清谷歌数据的下载方法:谷歌卫星影像免费下载?!来看这些软件,今天我们就来看看10米分辨率的Sentinel-2下载及简单的处理。
Sentinel-2(哨兵2)系列卫星的主要载荷多光谱成像仪拥有13个谱段,从可见光到近红外再到短波红外,空间分辨率从10m到60m。由Sentinel 2A 和Sentinel 2B 两颗卫星组成,图像采集时间分辨率是每颗卫星10天,两颗互补,重访周期为5天。
下载地址为欧空局哥白尼数据中心:https://scihub.copernicus.eu/dhus/#/home(建议选择谷歌浏览器)
首先你需要注册账户。
欧空局注册非常快,只要填写正确的邮箱基本都是立马给你发送邮件确认。
注册登录后可在右侧切换图层,选择矢量框选。
框选你所需要下载的区域点击
选择Sentinel-2
含云量可以自己设置,我们一般设置[0 TO 5]
可以点击
最后开始下载,最好科学上网。
下载速度还是可以的,500kb/s左右(如果下载时卡住,可取消重新下载)
▐ 数据预处理
Sentinel 2数据预处理就是对下载的L1C级产品进行辐射定标和大气校正生成L2A级产品。目前有二种方法:一种利用独立的Sen2Cor程序进行处理;另一种直接利用snap软件里Sen2Cor插件。推荐第一种做法。Sen2Rec插件主要用于将所有波段重采样至统一分辨率,当然也可以用SNAP自带的重采样工具,文中有具体介绍。
SNAP下载地址:http://step.esa.int/main/download/snap-download/
Sen2Cor下载地址:http://step.esa.int/main/snap-supported-plugins/sen2cor/sen2cor-v2-9/
Sen2Res下载地址:https://step.esa.int/main/snap-supported-plugins/sen2res/
一、安装Sen2Cor
下载后解压,双击运行L2A_Process.bat批处理文件,闪现窗口,不用管;
接着WIN+R输入“cmd”进入“命令提示符”,定位到Sen2Cor解压后的目录:比如我的输入:cd C:\Users\DELL\Sen2Cor-02.09.00-win64;
输入L2A_Process --help,查看是否安装成功。
如果未出现。删除解压后的文件。重新解压到电脑用户文件夹下(比如:C:\Users\DELL\)。重复上面操作直至成功。
二、数据预处理
对L1C级产品进行辐射定标和大气校正生成L2A级产品。
开始数据处理:继续刚才的命令提示符,输入“L2A_Process +(下载的sentinel 2号数据目录)”注意目录到.SAFE后缀,比如我的处理:L2A_Process E:\S2A_MSIL1C_20211130T030101_N0301_R032_T50TMK_20211130T050208.SAFE(目录路径不要有中文和空格)
处理完成后会在同样目录下生成一个处理后的L2A级数据,如果你要指定生成目录就在刚才的命令后加”--output_dir=路径“,正如刚才帮助中列出的提示。
如果是批处理,先将所有数据放在一个目录。
直接输入:
“for /D %s in (你的数据目录,比如我的:E: \S2A_MSIL1C*) do L2A_process %s”
到这里预处理已经完成,我们可以打开安装的SNAP(第一次打开提示更新插件选择是),分别加载L2A和L1C产品进行对比。
打开2个级别影像后,横向对比(B8波段中植被的光谱曲线)。
已经完成辐射校正和大气校正。
三、统一分辨率(重采样)
由于波段分辨率并不一致,所以你需要重采样至10m分辨率,让其一致以便后续分析。SNAP有2个自带的重采样功能:
推荐第一个,性价比最高,记得选择分辨率和输出路径,重采样后,右击重采样的产品打开RGB影像,你会发现可选择项变多了,统一分辨率都可以融合了。
当然如果你对重采样有更高要求,可以用更专业的插件Sen2Res,它的效果更好,但是需要更多的处理时间和空间,进入第三步操作。
四、加载Sen2Res插件
添加Sen2Res插件:安装SNAP好后,打开SNAP,Tools——>Plugins.加载后,勾选安装,安装好后重新启动SNAP软件。
打开Optical—>Sentinel-2 super resolution工具,指定输出路径,其他参数默认,默默等待重采样完成。
打开SNAP,Tools—>Plugins,切换到Available Plugins选项,搜索Sen2;选择下图中的280版本安装。安装完后,选择重启软件选项。
重启软件后。先打开需要处理的数据
再按照下图打开Sen2Cor280工具
指定分辨率All,直接生成L2A级产品。
重采样步骤同第一种方法一样。
授权转载:GIS前沿
END