查看原文
其他

全球10米遥感数据免费下载及预处理(附软件下载)

小助手 GIS前沿 2022-12-04


关于各种数据下载的方法我们分享过很多:各种实用航测遥感数据数据免费获取,速来领取这次看标题估计大家都猜出来了是Sentinel-2(哨兵2),没错,就是Sentinel-2。(这次重新整理,更加完善)


今天我们就来看看10米分辨率的Sentinel-2下载及简单的处理。





数据下载


Sentinel-2(哨兵2)系列卫星的主要载荷多光谱成像仪拥有13个谱段,从可见光到近红外再到短波红外,空间分辨率从10m到60m。由Sentinel 2A 和Sentinel 2B 两颗卫星组成,图像采集时间分辨率是每颗卫星10天,两颗互补,重访周期为5天。


下载Sentinel-2数据网站有很多,欧空局官网、USGS、地理空间数据云等一些国内外的网站都可以下载。


我们就以欧空局官网讲解。


下载地址为欧空局哥白尼数据中心:https://scihub.copernicus.eu/dhus/#/home(建议选择谷歌浏览器)


下载前我们先了解下Sentinel-2的数据


目前我们可下载的数据有,Sentinel-2的L1B,L1C和L2A数据。


具体可看下产品级别介绍:



L1C级的数据未经大气校正,需要用户自行进行大气校正处理。L2A级数据则是在L1C基础上经过大气校正的大气底层反射率数据,我们不再需要进行大气校正。


不过,后面有L1C级数据预处理流程。

首先你需要注册账户。


欧空局注册非常快,只要填写正确的邮箱基本都是立马给你发送邮件确认。



注册登录后可在右侧切换图层,选择矢量框选。



框选你所需要下载的区域点击按钮可以选择拍摄时间,含云量等等。



选择Sentinel-2

含云量可以自己设置,我们一般设置[0 TO 5]


如果,我们只想下载L2A数据,则可在Product Type中选着S2MSI2A,这样就是筛选后的L2A数据。



如果不做数据筛选,也可根据数据名来进行下载。



可点击预览数据情况。


最后开始下载,最好科学上网。



下载速度还是可以的,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程序)


一、安装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工具,指定输出路径,其他参数默认,默默等待重采样完成。


 

第二种方法(利用Sen2Cor插件)



打开SNAP,Tools—>Plugins,切换到Available Plugins选项,搜索Sen2;选择下图中的280版本安装。安装完后,选择重启软件选项。


 

重启软件后。先打开需要处理的数据


 

再按照下图打开Sen2Cor280工具

 


指定分辨率All,直接生成L2A级产品。


 

重采样步骤同第一种方法一样。


软件下载

进入本公众号后台

发送数字

1019

获取下载链接

文件名:Sentinel 2预处理


- END -


影像三维地形图制作(附练习数据)

龚健雅院士:遥感大数据与智能解译

遥感卫星影像处理流程

ArcGIS土地利用转移矩阵(附土地利用数据下载)

干掉卫星、无人机?只用气球就能获取10厘米分辨率遥感影像

您可能也对以下帖子感兴趣

文章有问题?点此查看未经处理的缓存