首页 > 其他分享 >【RS】Sentinel-2哨兵二号数据预处理(大气校正、重采样、波段组合)

【RS】Sentinel-2哨兵二号数据预处理(大气校正、重采样、波段组合)

时间:2024-03-06 18:03:20浏览次数:43  
标签:采样 L2A RS SNAP ENVI 波段 Sentinel 数据 预处理

        刚分享过新版欧空局的数据下载教程,正好把哨兵2号预处理的教程也分享一下,主要就是使用官方插件Sen2or对L1C级数据进行大气校正,使用SNAP软件对L2A级数据进行重采样并导出ENVI可以打开的格式,最后使用ENVI对重采样后的数据进行波段组合,以便于后期的定量分析。

1.软件下载

1.1 Sen2Cor下载

用途:L1C级数据大气校正生成L2A级数据

网站:Sen2Cor – STEP (esa.int)

百度云盘:Sen2Cor– RS迷途小书童

1.2 SNAP下载

用途:欧空局官方公布的哨兵数据处理软件,这里主要使用其波段重采样以及波段导出ENVI格式。

网站:SNAP Download – STEP (esa.int)

百度云盘:SNAP– RS迷途小书童

2.Sen2or大气校正

2.1 单文件处理

        使用命令提示符,Win+R,输入cmd(命令提示符)。使用命令cd /d 路径,将目录切换至Sen2Cor的目录。如下:

        输入命令L2A_Process.bat 影像路径,注意这里的影像路径必须子目录就是文件,而不是还有一层,解压的时候可能会生成两级目录,切记!!!还有就是不能有中文路径!!!如下:

        等待处理成功,返回Successfully即可。

2.2 批量处理

        使用命令提示符,Win+R,输入cmd(命令提示符)。使用命令cd /d 路径,将目录切换至Sen2Cor的目录。(与单文件一样,目的是调用L2A_Process插件)

        输入命令for /D %s in (G:\1\S2A_MSIL1C*) do L2A_process %s,其中G:\1\是影像的路径,S2A_MSIL1C*是指路径中所有以S2A_MSIL1C开头的文件,如果是S2B数据,就改下S2B_MSIL1C*即可。如下:

3.SNAP重采样

3.1 打开影像

        在SNAP中,点击左上角的File,使用Open product打开L2A数据文件夹中的MTD_MSIL2A.xml文件即可。

3.2 重采样

        打开Optical>Geometric>S2 Resampling Processor工具,设置导出格式为ENVI,修改导出路径,修改重采样的分辨率为10m,run即可。

4.ENVI波段组合

4.1 加载波段

        我们将重采样完成后文件夹中的.img波段文件拖入ENVI打开。

4.2 波段组合

        ENVI5.6以上版本使用Build Layer Stack工具,ENVI5.3使用Layer Stacking工具,将波段按顺序依次倒入,坐标系通过影像读取,输出即可!

5.总结

        以上即为Sentinel2数据的预处理流程,如果是从官网下载的L2A级数据,就直接重采样+波段组合即可,无需进行大气校正!此外ENVI5.6过后支持打开Sentinel2数据,所以SNAP导出ENVI格式主要是用来ENVI5.3的,我们这里使用SNAP是用来重采样。

        有什么问题直接留言沟通即可,后续我可能会更新使用Python一键批量预处理L1C级数据。感兴趣的话可以关注一下。

标签:采样,L2A,RS,SNAP,ENVI,波段,Sentinel,数据,预处理
From: https://www.cnblogs.com/RSran/p/18057204

相关文章

  • Caused by: org.apache.ibatis.reflection.ReflectionException: Error parsing prope
    问题描述1、编译过程是不会被发现的;2、运行时报错:Causedby:org.apache.ibatis.reflection.ReflectionException:Errorparsingpropertyname'lambda$pages$4d14352e$1'.Didn'tstartwith'is','get'or'set'。问题分析MyBatis-Plus的条件构造器常用......
  • CF1935E Distance Learning Courses in MAC
    CF1935EDistanceLearningCoursesinMAC题目大意给定\(n\)个变量\(z_i\in[x_i,y_i]\),你可以在范围内任意指定\(z_i\)的值。\(q\)次查询,每次查询给定区间\([l_i,r_i]\),求用这些变量得到的二进制或最大值。思路选择\(z\in[x,y]\),贡献分为两部分(1)\([x,y]\)的......
  • PostgerSql 相较于MSSQL 要注意的地方
    十年河东,十年河西,莫欺少年穷学无止境,精益求精1、大小写敏感传统的SQLSERVER数据库默认大小写不敏感,但PgSql是大小写敏感的,如下:2、varChar和Text其实没你想的那么简单MSSQL中用得最多的的文本类型是NVARCHAR,这是一个带长度限制的文本类型,对应地,PG中有VARCHAR,这样用没问题,......
  • 【RS】最新欧空局Sentinel-2卫星数据下载(哨兵1、2、3、5P、6系列)
    ​        之前分享过Sentinel2数据下载的方法,但是有粉丝反应欧空局的网站更新了,之前的网站已经不能用了。所以自己抽空研究了一下新版的欧空局网站,今天就和大家分享一下如何使用新版的欧空局网站下载哨兵系列的卫星数据,本文以Sentinel2数据下载为例。1.前言   ......
  • [转帖]linux-windows文件实时同步:Rsync使用教程
    http://luomuren.top/articles/2021/04/06/1617641017252.html#:~:text=linux-windows%E6%96%87%E4%BB%B6%E5%AE%9E%E6%97%B6%E5%90%8C%E6%AD%A5%EF%BC%9ARsync%E4%BD%BF%E7%94%A8%E6%95%99%E7%A8%8B%201%20%E4%B8%80%20%E3%80%81%E4%BB%80%E4%B9%88%E6%98%AFrync%20%EF%BC%9F......
  • DW_apb_timers介绍 --20240305
    DW_apb_timers是新思设计的基于APB总线的一个IP DW_apb_timers框图如下:DW_apb_timers特性:1、多达8个可编程定时器;2、定时器宽度可设:8至32位;3、支持两种运行模式:freerunning和user-mode;4、支持定时器的独立计时;5、每个中断可配置极性:低电平(下降沿),高电平(上升沿);6、单个或......
  • m基于FPGA的RS+卷积级联编译码实现,RS用IP核实现,卷积用verilog实现,包含testbench测
    1.算法仿真效果Vivado2019.2仿真结果如下:  2.算法涉及理论知识概要        级联码是一种通过将两种或多种纠错码结合使用来提高纠错能力的编码方案。在RS+卷积级联编码中,通常首先使用卷积码对原始数据进行编码,以增加冗余并提供一定的纠错能力。然后,将卷积码的输......
  • CF622F The Sum of the k-th Powers 题解
    原式为\(k+1\)次多项式,所以需要\(k+2\)个点确定。然后转化,前缀和。\[\begin{equation}n=k+2\\\end{equation}\]\[\begin{equation}f(x)=\sum\limits_{i=0}^{n}y_i\prod\limits_{j=0,j\nei}^{n}\frac{x-x_j}{x_i-x_j}\end{equation}\]\[\begin{equation}x_0=......
  • Git - error:you need to resolve your current index first 解决方案
    场景:从dev-test分支上拉取dev分支上的代码(意外操作,本应该拉取dev-test分支)相当于从一个分支A,切换到分支B,对B分支进行了pull的操作错误提示:error:youneedtoresolveyourcurrentindexfirst原因:在执行pull操作时,实际是执行了:fetch+merge两个操作。由于分支B很久未......
  • 无人直播,云SRS,Docker一键搭建安装教程
    如果有服务器并已经安装Docker,可以直接从第四步开始。第一步:选购轻量应用服务器并安装Ubuntu22或者20系统。境外服务器以香港和日本、新加坡为佳。24小时无人直播流量需求比较大,请按照视频码率自己算好,超出后流量非常昂贵。【以视频码率5000为例:5000/8/1024*3600*24*31/1024/102......