利用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进行批量大气校正相关推荐

  1. 使用Acolite批量大气校正Sentinel-2 L1C数据

    目录 1.前言 2.使用方法 3.批量校正程序 建了一个QQ群,大家可以在里边聊聊水色遥感数据下载和数据处理方面的事情:1087024529 1.前言 接下来的几篇博客,会介绍几种适用于浑浊水域的大气 ...

  2. GEE|在GEE对Sentinel-2、Landsat7、Landsat8进行批量大气校正、去云,并进行均值融合、裁剪、显示和数据下载

    前言: 由于我要使用S2在2016-2020年时序数据进行作物分类,但GEE仅提供2018年后的S2 L2A数据('COPERNICUS/S2_SR'),因此需要进行大气校正,同时由于需要的数据量大, ...

  3. GF1-WFV批量正射校正和批量大气校正

    RSD提供的GF1-WFV批量处理只包括下述两项 MSS数据的批量正射校正 MSS数据批量正射校正和大气校正 检查正射校正使用的DEM是否准备好 见博文<使用RSD批量融合高分1号(GF1)卫星 ...

  4. python大气校正_sen2cor批量大气校正Sentinel2数据——python代码

    上个博客介绍了几种Sentinel数据快速下载的方法,没看到的biubiubiu→→→→→→→→→→→→→→→→→→→→ 接下来分享一下如何在python下调用cmd使用sen2cor批量对Senti ...

  5. 利用Python在NASA Earthdata Search 上批量下载MODIS数据

    1.下载网址:https://search.earthdata.nasa.gov 2.第一步:注册账号,右上角Earthdata Login,如实注册,身份可以选择教育.用途可以选择科研团队,记得添加 ...

  6. 利用Python调用ECMWF欧洲中心API批量下载数据

    前段时间由于需要下载ECMWF(欧洲中期天气预报中心)的再分析数据,学习了如何利用Python调用ECMWF欧洲中心API进行批量下载.这种下载ECMWF数据的方法在官网上有非常详细的介绍.我只是对这 ...

  7. 使用RSD对高分1号卫星数据进行批量大气校正

    高分数据处理可以有下述5中选项 MSS数据的批量正射校正 MSS数据批量正射校正和大气校正 PAN 数据的批量正射校正 MSS与PAN批量正射校正和融合 MSS与PAN批量正射校正.MSS的大气校正和 ...

  8. 使用SeaDas对Sentinel-3(OLCI)、Sentinel-2(MSI)进行大气校正

    目录 1.前言 2.程序安装 3.处理事例 建了一个QQ群,大家可以在里边聊聊水色遥感数据下载和数据处理方面的事情:1087024529 1.前言 最近,有小伙伴询问我关于使用SeaDas-OCSSW ...

  9. 利用python读取excel中邮箱进行批量群发,使用中国移动内部邮箱或139邮箱-html编辑的正文并且带附件。

    简介 此项目是通过读取excel文件内容(包含公司名称,公司类型,邮箱-每个公司对应多个邮箱)自动批量发送邮件到目标公司.亲测可用. 完整代码 import smtplib from email.mi ...

最新文章

  1. C#模糊查询绑定datagridview
  2. python调用dll传指针参数
  3. 运动检测(前景检测)之(二)混合高斯模型GMM
  4. C# Socket编程
  5. 【转】解决多文件编译中连接错误 multiple definition of*****
  6. servu无法显示远程文件夹_win10软件:解决文件资源管理器无法显示文件夹大小...
  7. 安卓调用系统拍照功能:1、启动拍照返回图片,2、启动拍照,图片存储在指定路径下
  8. C++ OpenCV 开发笔记整理分类合集【所有的相关记录,都整理在此】
  9. 禁掉win2003/2008服务 提高系统运行速度
  10. JAVA设计模式-11-代理模式(动态)(一)
  11. Papi酱: 左手奶瓶右手尿布× 左手papitube右手娱乐圈√
  12. 【其他】神舟精盾K470-i3 D3更换支持5GHz的双频无线网卡
  13. 【金融风控-贷款违约预测】数据挖掘学习:5.模型融合
  14. Second-Order Cone Programming(SOCP) 二阶锥规划
  15. Learning Memory-guided Normality for Anomaly Detection 代码解析
  16. Python爬虫实战: 爬取网易云歌单
  17. 《金刚经说什么》诗词辑录
  18. 用html写一个简单课表
  19. zend studio自动换行
  20. PHPExcel大文件导入数据库

热门文章

  1. 全球程序员收入出炉!北京收入排入全球第十
  2. DDS的verilog 实现个人总结
  3. 《蜘蛛侠3》5月2日中国提前上映
  4. access筛选末位数字为1或5_数字鸿沟理论不会用?北大必考常考,我们以快手短视频为例...
  5. ubuntu使用deepin-wine安装微信出现版本过低问题的解决
  6. 小程序利用canvas 绘制图案 (生成海报, 生成有特色的头像)
  7. Linux pv显示进度条命令
  8. Android入门小项目--微信登录界面源码(倒计时、fragement切换、activity信息传递)
  9. python爬京东延迟加载_python大规模爬取京东
  10. 推流是什么,直播为什么要推流