首页 > 编程语言 >Python依据遥感影像的分幅筛选出对应的栅格文件

Python依据遥感影像的分幅筛选出对应的栅格文件

时间:2024-06-01 13:01:40浏览次数:31  
标签:条带 Python 栅格 遥感 文件夹 分幅 path 影像

  本文介绍基于Python语言,结合已知研究区域中所覆盖的全部遥感影像分幅条带号,从大量的遥感影像文件中筛选落在这一研究区域中遥感影像文件的方法。

  首先,先来明确一下本文所需实现的需求。现已知一个研究区域(四川省),且已知覆盖这一研究区域所需的全部遥感影像的分幅条带号(这里就以Sentinel-2数据所采用的分幅规则为例来介绍),如下图所示。

image

  其中,矢量图层上的标签Label)就是所需的遥感影像的分幅条带号;且这一矢量要素的属性表中,有具体每一个分幅条带号的具体字段(如上图红色方框所示的那一列)。

  同时,我们还下载好了大量的遥感影像——这些遥感影像数量众多,既有位于这一研究区域内的遥感影像,同时还有位于研究区域外的。其中,这里所有的遥感影像都保存在一个大文件夹下,而这一大文件夹下还有多个子文件夹,具体遥感影像就存在于子文件夹中;如下图所示。

  其中,红色框内的文件夹即为前述“大文件夹”,紫色框内的每一个文件夹即为前述“子文件夹”。此外,每一个子文件夹的命名也是和遥感影像的分幅条带号对应的,比如43LCC这一个子文件夹就是保存了UTM编号为43的网格内的全部遥感影像。这里关于Sentinel-2遥感影像分幅条带号的介绍,如果大家有需要,可以参考哨兵2号Sentinel-2分幅条带介绍与MGRS网格矢量文件获取这篇文章。

  同时,在每一个遥感影像文件文件名中,都有分幅条带号这一字段;例如,打开43LCC这一子文件夹,如下图所示。

  我们希望实现的,就是基于上图所示的研究区域内遥感影像的分幅条带号信息,在大量的遥感影像中筛选出这些分幅对应的遥感影像文件,并将其通过复制的方式放入到一个新的结果文件夹中。

  了解了具体需求,我们就可以开始代码的撰写。首先,本文所需的完整代码如下所示。

# -*- coding: utf-8 -*-
"""
Created on Fri Feb  3 18:33:11 2023

@author: fkxxgis
"""

import os
import xlrd
import shutil

xls_file = r"E:\02_Project\Sentinel.xls"
tif_path = r"G:\S2_CSI_LCC_product\LCC_10days_2020"
result_path = r"G:\Select"

book = xlrd.open_workbook(xls_file)
sh = book.sheet_by_index(0)
row_num = sh.nrows

for row in range(1, row_num):
    mgrs = sh.cell(row,6).value
    utm = mgrs[0:2]
    for root, dirs, files in os.walk(tif_path):
        for name in dirs:
            if name[0:2] == utm:
                for root_1, dirs_1, files_1 in os.walk(os.path.join(root, name)):
                    for name_1 in files_1:
                        if name_1.find(mgrs) != -1:
                            if not os.path.exists(os.path.join(result_path, mgrs)):
                                os.makedirs(os.path.join(result_path, mgrs))
                                print("Make " + mgrs + " path.")
                            shutil.copy(os.path.join(root, name, name_1), os.path.join(result_path, mgrs, name_1))

  其中,我们需要导入osxlrdshutil3Python库,分别实现系统文件遍历、Excel数据读取与研究区域内遥感影像数据的复制。其中,关于xlrd库的配置,如果大家是用的Anaconda环境,可以参考Python表格处理模块xlrd在Anaconda中的安装这篇文章。此外,我们首先需要将本文开头提及的表示遥感影像数据分幅条带号的矢量文件的属性表导出(通过“System Toolboxes”→“Conversion Tools.tbx”→“Excel”→“Table To Excel”选项,如下图所示),保存为.xls格式的表格文件,从而方便我们后期代码对数据的读取。

  在代码中,xls_file表示我们将属性表导出后得到的.xls格式的表格文件,tif_path表示存放所有遥感影像的大文件夹,result_path则表示找到我们所需的结果遥感影像文件后,希望将其复制到的路径。

  随后,我们通过xlrd.open_workbook()函数打开.xls格式的表格文件,读取其中的数据并获取其行数row_num;接下来,我们即可对这一表格文件的单元格数据加以遍历,也就是代码中的第1for循环——我们从第2行开始(因为第1行是表头),读取每一行第7列的数据(也就是存放有遥感影像分幅条带号的那一列),从而依次获取所需的遥感影像条带号。

  接下来,由于遥感影像的分幅条带号的前两位,以及前述子文件夹名称的前两位,都是表示UTM编号的两位数字,因此我们通过utm = mgrs[0:2]这句代码,截取当前遥感影像分幅条带号的前两位,并通过os.walk()函数进行遍历,通过一个if判断语句找到大文件夹下对应的子文件夹;随后,在这一子文件夹中同样通过os.walk()函数,进行遥感影像文件的遍历——也就是代码中第2个、第3个、第4for循环所做的事情。

  接下来,针对遍历得到的子文件夹中每一个遥感影像文件,我们通过.find()函数找到所有满足这一分幅条带号的遥感影像文件,也就是代码中第5for循环与第2if判断语句的工作。找到符合当前分幅条带号的遥感影像文件后,我们即可开始将其复制到目标文件夹中。

  这里还有一个需求,因为我们这里保存的是多时相遥感影像数据(即每一个分幅条带号对应着多个不同时相的遥感影像文件),因此我们希望在目标文件夹中,同样用各个分幅条带号作为名称,创建多个子文件夹;然后将当前分幅条带号对应的全部遥感影像数据放入这一文件夹中。基于此,我们通过os.path.exists()函数判断是否存在指定的子文件夹,如果不存在的话就新建这一文件夹;随后,通过shutil.copy()函数复制这些遥感影像文件。

  运行上述代码,最终我们得到的结果如下图所示;其中,每一个子文件夹都表示一个分幅条带号,子文件夹内即为这一分幅条带号所对应的全部时相的遥感影像文件。

  至此,大功告成。

标签:条带,Python,栅格,遥感,文件夹,分幅,path,影像
From: https://www.cnblogs.com/fkxxgis/p/18225883

相关文章

  • 基于Python+OpenCV高速公路行驶车辆的速度检测系统
    欢迎大家点赞、收藏、关注、评论啦,由于篇幅有限,只展示了部分核心代码。文章目录一项目简介二、功能三、系统四.总结一项目简介  一、项目背景与意义随着交通流量的增加和高速公路的快速发展,高速公路上的车辆速度管理成为了保障道路安全和提升通行效率的重......
  • SockJS Python 客户端
    SockJS是一个用于浏览器和服务器之间建立全双工连接的库,它允许在不支持原生WebSocket的浏览器中提供类似WebSocket的API。Python中也有很多SockJS客户端库,例如`python-socketio`和`sockjs-client-py`。以下是如何使用Python客户端(在这个例子中,我们将使用`pytho......
  • 【多进程并发笔记】Python-Multiprocess
    目录调用函数后,函数内的变量如何释放?python2.7怎么使用多线程加速forloop多进程进程池,函数序列化错误的处理Time模块计算程序运行时间使用多进程,Start()后,如何获得返回值?使用多进程并行,每个进程都将结果写入sqlite3数据库,可以么python创建进程池进程池的最大进程数怎么确......
  • 【计算机毕业设计】谷物识别系统Python+人工智能深度学习+TensorFlow+卷积算法网络模
    谷物识别系统,本系统使用Python作为主要编程语言,通过TensorFlow搭建ResNet50卷积神经算法网络模型,通过对11种谷物图片数据集('大米','小米','燕麦','玉米渣','红豆','绿豆','花生仁','荞麦','黄豆','黑米','黑豆')进行训练......
  • 为什么 python 会出现这种行为?
    我试图在Python中将数字动态追加到2D数组中。temp=[]arr=[tempforiinrange(2)]Arr[0].append("erg;erg)arr[0].append("ergse")print(arr)我得到的输出结果是......
  • 基于python美食网站的设计与实现论文
    目录摘要IAbstractII第1章绪论11.1项目研究的背景11.2开发目的和意义11.3国内外研究现状1第2章系统开发工具32.1Python编程语言32.2B/S模式32.3MySQL数据库42.4Django框架介绍42.5Vue开发技术52.6JavaScript简介5第3章系统分析73.......
  • 基于Python高校成绩分析的设计与实现论文
    目录1绪论11.1开发背景11.2国内外研究现状和发展趋势综述11.3开发设计的意义及研究方向12系统开发技术32.1Python可视化技术32.2Django框架32.3IDEA介绍32.4B/S架构42.5MySQL数据库介绍43系统分析53.1整体分析53.2功能需求分析53.3系......
  • python参数传递与Java参数传递的对比
    目录一、本文简介二、参数传递的类型三、python中参数传递的类型3.1可变对象与不可变对象3.2对象引用传递和引用传递的区别3.3如何理解python中列表的"+"操作和append操作的区别四、Java中参数传递的类型附录一、本文简介python参数传递的类型是“对象引用传递”......
  • Python While 循环语句
    Python编程中while语句用于循环执行程序,即在某条件下,循环执行某段程序,以处理需要重复处理的相同任务。其基本形式为:```while判断条件(condition):  执行语句(statements)……```执行语句可以是单个语句或语句块。判断条件可以是任何表达式,任何非零、或非空(null)的值......
  • Python 运算符
    相关学习资料什么是运算符?本章节主要说明Python的运算符。举个简单的例子 4+5=9 。例子中,4 和 5 被称为操作数,+ 称为运算符。Python语言支持以下类型的运算符:算术运算符比较(关系)运算符赋值运算符逻辑运算符位运算符成员运算符身份运算符运算符优先级接下来让......