哨兵2号波段_Sentinel2 哨兵二号数据下载及处理教程

本文介绍了哨兵2号卫星的基本信息,详细讲解了如何从官方网站下载数据,以及如何使用Sen2Cor进行L1C到L2A的处理。还提到了使用SNAP软件进行数据格式转换,包括在ENVI5.1和ENVI5.3中处理哨兵2号数据的方法。
摘要由CSDN通过智能技术生成

一 、Sentinel-2 哨兵二号简介

Sentinel-2 是高分辨率多光谱成像卫星,一颗卫星的重访周期为10天,两颗互补,重访周期为5天。分为2A和2B两颗卫星。2A于2015年6月23日01:52 UTC以“织女星”运载火箭发射升空。2B于2017年3月07日北京时间9时49分 UTC以“织女星”运载火箭发射升空。两者同时进入运行状态后,每5天可完成一次对地球赤道地区的完整成像,而对于纬度较高的地区,这一周期仅需3天。

Sentinel-2 卫星携带一枚多光谱仪器(MSI),可覆盖13个光谱波段,地面分辨率分别有10m、20m和60m:

Sensor

Band number

Band name

Sentinel-2A

Sentinel-2B

Resolution (meters)

Central wavelength (nm)

Bandwidth (nm)

Central wavelength (nm)

Bandwidth (nm)

MSI

1

Coastal aerosol

443.9

20

442.3

20

60

MSI

2

Blue

496.6

65

492.1

65

10

MSI

3

Green

560.0

35

559

35

10

MSI

4

Red

664.5

30

665

30

10

MSI

5

Vegetation Red Edge

703.9

15

703.8

15

20

MSI

6

Vegetation Red Edge

740.2

15

739.1

15

20

MSI

7

Vegetation Red Edge

782.5

20

779.7

20

20

MSI

8

NIR

835.1

115

833

115

10

MSI

8b

Narrow NIR

864.8

20

864

20

20

MSI

9

Water vapour

945.0

20

943.2

20

60

MSI

10

SWIR –  Cirrus

1373.5

30

1376.9

30

60

MSI

11

SWIR

1613.7

90

1610.4

90

20

MSI

12

SWIR

2202.4

180

2185.7

180

20

Level-1C 是经过正射校正和几何精校正的大气表观反射率产品,并没有进行大气校正。

Level-2A 主要包含经过大气校正的大气底层反射率数据(Bottom-of-Atmosphere corrected reflectance),但这个L2A数据需要用户根据需求自行生产,为此,ESA发布了专门生产L2A级数据的插件Sen2cor。

二、 哨兵数据下载

https://scihub.copernicus.eu/  官网下载,需要先注册

(浏览器建议选择用google浏览器或火狐浏览器)

fe03968128f9cf014f13aea0f5568740.png

544fcdd41f94c442d283b39327a6f72d.png

通过鼠标左键和滚轮将区域缩放到感兴趣区域(我们以保定区域为例)

e0a615a4f869cf34fe321d7dbb75e9c2.png

绘制感兴趣区域, 点击绘图工具,在浏览器端绘制出感兴趣区域

4f2a0e2cfe269e0e15c5eafa93dbb152.png

高级搜索:输入要检索的时间和传感器类型,点击搜索

cd4f3bdb2f7170984ffaa06259e8355d.png

搜索结果

b08c877e2a95a78d16a03a5286c3b5e8.png

我们一般选择云量较少的影像

点击,可以对影像进行预览

bbe20da66ff7cc62fc08ac47e43e325d.png

a6ec6e6cc62e43f786f7eb1029d9aa90.png

点击 下载 进行影像数据下载(下载时间段最好是凌晨4-8点钟,白天其他时间段下载影像数据非常慢,速度只有50KB/s左右,而且只能同时下载2幅影像,同时下载多幅影像会出错!),

4f0333740dc95f821e1d355a33fc7e19.png

e9d989fb3a366d1914792c48e3c142d1.png

三、影像预处理

下载后的哨兵2号数据产品级别为Level-1C数据。要转换成L2A级产品需要用户自己进行处理生产。(L2A级数据主要包含经过辐射定标和大气校正的大气底层反射率数据(Bottom-of-Atmosphere corrected reflectance))(当然也可以下载L2A级别影像,只不过数量较少)。

1e1e7c3955532d9400cab92b98395007.png

1. L1C->L2A

哨兵数据预处理(L1C级产品处理为L2A级)

http://step.esa.int/main/third-party-plugins-2/sen2cor/sen2cor_v2-8/

65aecdd7f2e82349d275faa948ecf748.png

因为我的电脑是Windows系统,所以选择 第一个,下载完毕后解压文件

生成一个Sen2Cor-02.08.00-win64文件夹

ef85d31276f16dad8c2ec524b655bcb8.png

将这个文件夹复制拷贝到 c:\User文件夹下面

5f9c09d11ee3f5e05cc8c49e92192a80.png

将下载的哨兵-2影像进行解压

2e17433536913a892ad20c5f36153b4e.png

打开命令行(Win+R)

切换到C:\Users\Sen2Cor-02.08.00-win64,输入 L2A_Process.bat   --help

如果输出一下信息,则表明工具可用

087d6301780e34d2306906081c0502bf.png

在命令行中输入 L2A_Process.bat ,并将要处理的文件夹拖入命令行中,如下图所示,回车,进行处理。。。

37021c68bb2db211cad146c8d24f684f.png

等进度显示100%的时候表示运行结束(这个过程预计持续20-30分钟),文件同一目录会产生

d6b7cc90a391b34ba86e9f086b116a52.png

0c0367fbde6ee0bda6caeb9f52e09a62.png

S2B_MSIL2A的文件夹,即为预处理后的影像(L1C->L2A级别)。

2 哨兵数据格式转换

哨兵的格式转换主要通过SNAP软件进行处理

http://step.esa.int/main/download/

下载后进行安装

e9499249347babaaf368e277773c37c6.png

44405f19864949da763f7f50e2a5d267.png

打开以MTD开头的以.XML结尾的文件

0f930d21a53c30f50a34ddb1e525de44.png

在Product Explorer中可以看到相关数据信息

dc3f63561aafe7f5c47f957737c3eb68.png

右键 导入的文件—>Open RGB Image Window,

afafcf43cc3f8863546947ac17585dad.png

c8bb06342f74adaa2ea02aa9abf842c8.png

2cfd40017bf12cf3458af6725d3db8ce.png

257b4db3d6f29bd98d32da971e2c0621.png

影像分辨率还是非常高的。

接下来进行文件转换

2.1 ENVI5.1

ENVI5.1无法直接打开Sentinel数据,但大家对ENVI操作会更加熟练,这里补充说明下如何用ENVI5.1打开Sentinel数据,主要就是重采样并转换格式。

打开重采样工具

b7489b05b36f1f9dab239c3a3cfef073.png

8a5b91eadcfaaaf7d7311b4353bc6404.png

854db32c985359960b1898b1da3e1bc8.png

这里建议以10m重采样。若以B1波段重采样的话,结果所有波段都会为60m分辨率,数据质量会下降。而设置为10m或者以B2波段重采样,结果波段的空间分辨率变为10m,这对于非10m分辨率的波段,虽然分辨率提高了,但实质并无变化。

重采样后的结果就可以导入ENVI5.1中显示了,但其存储方式是一个波段一个文件,在ENVI中融合一下(ENVI中搜索layerstacking工具)就可以得到多波段文件了。

2.2 ENVI5.3

ENVI5.3 可以直接打开哨兵的影像(jp2格式)

打开

F:\2020_11_29\S2A\S2A_MSIL2A_20200919T030551_N0214_R075_T50SLH_20200919T060633\S2A_MSIL2A_20200919T030551_N0214_R075_T50SLH_20200919T060633.SAFE\GRANULE\L2A_T50SLH_A027385_20200919T031126\IMG_DATA\R10m

edb1b8813836c2e77e609437e2016079.png

里面存储的为R10M分辨率的波段信息

c602a15049474d481be69c5a89b86e52.png

将JP2格式的波段文件,另存为ENVI支持的.dat文件

通过Layer stacking工具进行合成成新的文件,进行后期数据处理和分析工作.

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值