使用Acolite批量大气校正Sentinel-2 L1C数据
目录
1.前言
2.使用方法
3.批量校正程序
建了一个QQ群,大家可以在里边聊聊水色遥感数据下载和数据处理方面的事情:1087024529
1.前言
接下来的几篇博客,会介绍几种适用于浑浊水域的大气校正方法如何使用,以及使用这些方法批量大气校正水色数据。
Acolite大气校正方法,在浑浊水域有着较高的大气校正精度。可用于Landsat (5/7/8) and Sentinel-2 (A/B) 数据。程序下载地址:
官网:http://odnature.naturalsciences.be/remsem/software-and-data/acolite
Github:https://github.com/acolite/acolite
在官网会有Acolite的使用说明书,建议大家从前到后看一遍操作说明。
2.使用方法
在官网下载好Acolite压缩包后,解压,可以看到有这四个文件。点击acolite.exe,就可以直接运行。
Acolite提供了两种方式运行,一种是使用GUI界面,一种是在命令行中运行。点击acolite.exe后,会跳出GUI界面
2.1.GUI界面中运行
使用GUI界面进行大气校正时,可以看到可以设置的参数较少。需要设置更多参数的,要修改setting文件,然后使用命令行调用setting文件进行处理。
关于setting文件,里边包含的是关于大气校正时需要的参数。setting文件的内容格式是怎么样的?可以首先使用GUI界面处理一景影像,就会生成setting文件,然后根据需要修改setting文件中的参数。都可以设置哪些参数,建议大家阅读官方说明文档。
以下只介绍GUI界面中的参数设置:
(1)输入,输出文件的路径;
(2)兴趣区范围,即最南、最北、最西、最东边的度数(十进制),可以不设置,不设置就输出原始影像的范围;
(3)L2W parameters: 即需要输出的L2W的参数。可以填写很多,详见官方说明文档;
(4)PNG outputs:一些生成的数据可以输出为PNG格式的图片;
(6)Save or restore settings: Save是可以保存上方的参数到一个setting文件。Restore是导入之前已经生成好的Setting文件(没有的话就不进行);
(7)Run processing就是开始执行,Exit就是退出程序。
2.2.命令行中运行
主要是调用Setting文件。然后进行处理。在命令行中输入:
> acolite安装路径 --cli --settings=setting文件路径
3.批量校正程序
主要的思路是:首先通过GUI界面生成一个Setting文件,然后根据要处理的数据,将Setting文件中的inputfile和outfile对应的路径名给更改掉。然后acolite根据setting文件进行批量处理。程序其实比较简单的。
#-*-coding: utf-8-*-
#@Author: zhaohaiyang
#@E-mail: hyzhao_rs@163.com
#使用Acolite进行大气校正批处理
#
#
#
################################################import os
from subprocess import check_call#--------------------------------------------------------
#需要修改的变量
example='G:/settingfileexample.txt' #一个Setting文件例子 #注意setting文件的第91行为输入文件路径(从0开始)#第92行为输出文件的路径。运行之前应该先检查,不是的话根据实际情况修改
s2dir='G:/Sentinel-2/' #原始数据的存放目录,注意加上最后的/
resultDir='G:/Result/' #Acolite大气校正后结果的存放路径。注意加上最后的/
acolitepath='D:/Program Files/acolite_py_win/dist/acolite/acolite.exe' #acolite的路径
#--------------------------------------------------------#读出例子Setting文件的内容
with open(example,'r') as ef:examplecon=ef.read().split('\n')[:-1]s2fileList=os.listdir(s2dir)
settingList=[] #存放各setting文件的路径for s2f in s2fileList:fp1=os.listdir(s2dir+s2f)[0]inputfile=s2dir+s2f+'/'+fp1outputfile=resultDir+s2fos.mkdir(outputfile)examplecon[91]='inputfile='+inputfileexamplecon[92]='output='+outputfilesettingFile=resultDir+s2f+'_setting.txt'settingList.append(settingFile)with open(settingFile,'w') as outsetting:for ec in examplecon:outsetting.write(ec+'\n')for stl in settingList:cmdline=acolitepath+' --cli --settings='+stlcallout=check_call(cmdline)
使用Acolite批量大气校正Sentinel-2 L1C数据相关推荐
- 利用python和Sen2cor对Sentinel2进行批量大气校正
利用python和Sen2cor对Sentinel2进行批量大气校正 p.s本文只适合大量的批量校正,一般两三幅影像的那种就直接搜批量校正就好了. 废话不多说,直接上代码 import os impo ...
- GEE|在GEE对Sentinel-2、Landsat7、Landsat8进行批量大气校正、去云,并进行均值融合、裁剪、显示和数据下载
前言: 由于我要使用S2在2016-2020年时序数据进行作物分类,但GEE仅提供2018年后的S2 L2A数据('COPERNICUS/S2_SR'),因此需要进行大气校正,同时由于需要的数据量大, ...
- GF1-WFV批量正射校正和批量大气校正
RSD提供的GF1-WFV批量处理只包括下述两项 MSS数据的批量正射校正 MSS数据批量正射校正和大气校正 检查正射校正使用的DEM是否准备好 见博文<使用RSD批量融合高分1号(GF1)卫星 ...
- python大气校正_sen2cor批量大气校正Sentinel2数据——python代码
上个博客介绍了几种Sentinel数据快速下载的方法,没看到的biubiubiu→→→→→→→→→→→→→→→→→→→→ 接下来分享一下如何在python下调用cmd使用sen2cor批量对Senti ...
- 使用RSD对高分1号卫星数据进行批量大气校正
高分数据处理可以有下述5中选项 MSS数据的批量正射校正 MSS数据批量正射校正和大气校正 PAN 数据的批量正射校正 MSS与PAN批量正射校正和融合 MSS与PAN批量正射校正.MSS的大气校正和 ...
- Sentinel-2批量大气校正_基于Sen2Cor
将sentinel-2数据放在一个文件夹下(S2A和S2B数据要分开处理) 1. 2.for /D %s in (E:\1\S2A_MSIL1C*) do L2A_process %s for /D ...
- 大气校正后的ndvi_Sentinel2 L1C下载、大气校正、重采样
点击蓝字关注我哦 1.基本信息(成像仪/重访周期/波段数/分辨率) 哨兵2号是高分辨率多光谱成像卫星,携带一枚多光谱成像仪(MSI),用于陆地监测,可提供植被.土壤和水覆盖.内陆水路及海岸区域等图像, ...
- 使用ACOLITE对OLCI影像进行大气校正+批处理
目录 下载软件 处理单幅OLCI影像 批处理思路及代码示例 建了一个QQ群,大家可以在里边聊聊水色遥感数据下载和数据处理方面的事情:1087024529 喜大普奔,ACOLITE可以大气校正OLCI传 ...
- sentinel-2 大气校正方法
参考文章: 哨兵2号(Sentinel-2)介绍.下载.预处理及批处理_GeoSuper的博客-CSDN博客_哨兵2号 Sentinel-2批量大气校正_基于Sen2Cor_GeoSuper的博客-C ...
最新文章
- 血泪史:阿里云+ubuntu+vnc+xfce4
- Android应用开发经常使用知识
- SharePoint Add-in Model (App Model) 介绍 – 概念、托管方式、开发语言
- 扩展欧几里得求解ax+by=c的特殊解(模板)
- 循环语句 for循环、while循环、do while循环
- cookies与session
- 算法概论第八章课后习题8.3
- hp t410微型计算机使用,HP 发表新款 t410 AIO Smart Zero 精简型电脑,仅需网络线即可作为电源驱动使用...
- 闪迪u盘量产工具万能版_加密保护,更高性能!闪迪酷奂CZ74 U盘新品果然够惊艳...
- java的小于等于符号怎么打_「小于符号」mybatis的一些特殊符号标识(大于,小于,等于,不等于) - seo实验室...
- AUTOCAD——成组命令
- 中国志愿者服务器注册,如何注册成为志愿者?中国志愿服务网注册流程
- 入行手游必知的8项基础知识
- 信息系统项目管理师-3项目立项管理
- 最大连续登录天数统计对应的用户数分布的三种实现方式(HQL+Spark Rdd)
- python设置excel单元格格式_Python帮你做Excel——格式设置与画图
- 简单易懂的计算机网络相关名词的解释--集线器篇(持续更新)
- IC 拔取器 rework station
- pytho lambda表达式
- 基于vue3+ts+scss的后台管理系统(一)
热门文章
- win10+安装Python 3.9.0b4+anaconda3
- 手机投屏时电视显示服务器有问题,小屏变大屏,手机投屏这几招你学废了么?...
- 刺激战场压枪——外挂
- 在线生成广西geojson数据+精确到乡镇街道级别手绘地图json
- HTML5 Day02 知识点总结
- php关闭warn,PHP warn、notice 有关问题的解决方案
- 编程程序,提示用户输入两个点(x1,y1)、(x2,y2),然后显示两点之间的距离。使用Math.pow()
- flex 实现菜单九宫格布局
- 英语语法快速入门2--并列句(附思维导图)
- APP上架实战干货腾讯《应用宝》