利用python和Sen2cor对Sentinel2进行批量大气校正
利用python和Sen2cor对Sentinel2进行批量大气校正
p.s本文只适合大量的批量校正,一般两三幅影像的那种就直接搜批量校正就好了。
废话不多说,直接上代码
import os
import shutil
import subprocessfile_path = "E:\\ESRIRS\\Sentinel2Imag\\"
year_list = ["2016", "2017", "2018", "2019", "2020", "2021", "2022"]
satellite_list = ["Sentinel2A", "Sentinel2B"]def readNameWriteInTxt():for year in year_list:for satellite in satellite_list:file_path2 = file_path + year + "\\" + satellitetxt_file_path = file_path + satellite + year + ".txt"names = os.listdir(file_path2)if names is None:names += "Null"else:with open(txt_file_path, "w+") as log:log.write(file_path2 + "\n")for name in names:name = name + "\n"log.write(name)def AutomationAtmosphericCorrectionC(satellite_listA, year_listA, bat_sen2cor):for satellite in satellite_listA:for year in year_listA:txt_file_path = foder_path + satellite + year + ".txt"with open(txt_file_path, "r+") as txt_file:txt_first_line = next(txt_file)lines = txt_file.readlines()#index = 1for line in lines:if line != "":cmd_args = [bat_sen2cor, txt_first_line.strip("\n") + "\\" + line.strip("\n")]print("{} processing begin!".format(line.strip("\n")))subprocess.call(cmd_args)print("{} processing finished!\n".format(line.strip("\n")))line = ""files = open(txt_file_path, "w+")files.write(txt_first_line)line = files.writelines(line)files.close()
readNameWriteInTxt是读取file_path 目录下所有文件夹,并将其保存在txt文件中。txt第一行是影像父文件夹的绝对路径,除第一行外是影像的文件名。.
AutomationAtmosphericCorrectionC是逐行读取之前创建的txt文件,并对影像进行校正,每校正一次影像就会删除一条影像的文件名,这样第一天没校正完的位置也会被记录,第二天再运行脚本即可继续校正。
bat_sen2cor是Sen2cor的绝对路径
楼主电脑是九代i5标压版,一般跑十一二个小时能校正完二十幅影像。
p.s2016及以前的sentinel影像的校正好像有问题,不能用05版sne2cor,目前搜到的时要使用10版的,不过楼主没试验过。
利用python和Sen2cor对Sentinel2进行批量大气校正相关推荐
- 使用Acolite批量大气校正Sentinel-2 L1C数据
目录 1.前言 2.使用方法 3.批量校正程序 建了一个QQ群,大家可以在里边聊聊水色遥感数据下载和数据处理方面的事情:1087024529 1.前言 接下来的几篇博客,会介绍几种适用于浑浊水域的大气 ...
- 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 ...
- 利用Python在NASA Earthdata Search 上批量下载MODIS数据
1.下载网址:https://search.earthdata.nasa.gov 2.第一步:注册账号,右上角Earthdata Login,如实注册,身份可以选择教育.用途可以选择科研团队,记得添加 ...
- 利用Python调用ECMWF欧洲中心API批量下载数据
前段时间由于需要下载ECMWF(欧洲中期天气预报中心)的再分析数据,学习了如何利用Python调用ECMWF欧洲中心API进行批量下载.这种下载ECMWF数据的方法在官网上有非常详细的介绍.我只是对这 ...
- 使用RSD对高分1号卫星数据进行批量大气校正
高分数据处理可以有下述5中选项 MSS数据的批量正射校正 MSS数据批量正射校正和大气校正 PAN 数据的批量正射校正 MSS与PAN批量正射校正和融合 MSS与PAN批量正射校正.MSS的大气校正和 ...
- 使用SeaDas对Sentinel-3(OLCI)、Sentinel-2(MSI)进行大气校正
目录 1.前言 2.程序安装 3.处理事例 建了一个QQ群,大家可以在里边聊聊水色遥感数据下载和数据处理方面的事情:1087024529 1.前言 最近,有小伙伴询问我关于使用SeaDas-OCSSW ...
- 利用python读取excel中邮箱进行批量群发,使用中国移动内部邮箱或139邮箱-html编辑的正文并且带附件。
简介 此项目是通过读取excel文件内容(包含公司名称,公司类型,邮箱-每个公司对应多个邮箱)自动批量发送邮件到目标公司.亲测可用. 完整代码 import smtplib from email.mi ...
最新文章
- C#模糊查询绑定datagridview
- python调用dll传指针参数
- 运动检测(前景检测)之(二)混合高斯模型GMM
- C# Socket编程
- 【转】解决多文件编译中连接错误 multiple definition of*****
- servu无法显示远程文件夹_win10软件:解决文件资源管理器无法显示文件夹大小...
- 安卓调用系统拍照功能:1、启动拍照返回图片,2、启动拍照,图片存储在指定路径下
- C++ OpenCV 开发笔记整理分类合集【所有的相关记录,都整理在此】
- 禁掉win2003/2008服务 提高系统运行速度
- JAVA设计模式-11-代理模式(动态)(一)
- Papi酱: 左手奶瓶右手尿布× 左手papitube右手娱乐圈√
- 【其他】神舟精盾K470-i3 D3更换支持5GHz的双频无线网卡
- 【金融风控-贷款违约预测】数据挖掘学习:5.模型融合
- Second-Order Cone Programming(SOCP) 二阶锥规划
- Learning Memory-guided Normality for Anomaly Detection 代码解析
- Python爬虫实战: 爬取网易云歌单
- 《金刚经说什么》诗词辑录
- 用html写一个简单课表
- zend studio自动换行
- PHPExcel大文件导入数据库
热门文章
- 全球程序员收入出炉!北京收入排入全球第十
- DDS的verilog 实现个人总结
- 《蜘蛛侠3》5月2日中国提前上映
- access筛选末位数字为1或5_数字鸿沟理论不会用?北大必考常考,我们以快手短视频为例...
- ubuntu使用deepin-wine安装微信出现版本过低问题的解决
- 小程序利用canvas 绘制图案 (生成海报, 生成有特色的头像)
- Linux pv显示进度条命令
- Android入门小项目--微信登录界面源码(倒计时、fragement切换、activity信息传递)
- python爬京东延迟加载_python大规模爬取京东
- 推流是什么,直播为什么要推流