Python中ArcPy切割栅格文件为四等分

  本文介绍基于Python 中的ArcPy模块,基于一个大文件夹 ,遍历其中每一个子文件夹 中所有的遥感影像栅格文件 ,并将原本的每一景 遥感影像文件四等分 切割,或裁剪为其他指定个数的小块的方法。

  首先明确一下我们的需求。现有一个大文件夹 ,其中具有多个子文件夹,如下图所示。

  其中,每一个子文件夹 中,都含有大量的栅格遥感影像文件 (我们这里就以.tif格式的栅格图像文件为例来介绍);例如,随意打开上图中任意一个子文件夹,都是如下图所示的情况。

  我们希望实现的是,将每一个子文件夹 中的每一景遥感影像 ,都加以裁剪、切割,将原本的一景遥感影像切分为新的4份,相当于是四等分 。这里大家想将原本的图像分为几部分都可以 ,即可以将原来的一景图像分为N * M个部分。

  了解了需求,我们就可以开始代码的撰写。本文所用代码,其实和我们之前的文章# ArcPy批量切割遥感影像为大量小矩形,以及# Python ArcPy依据矢量面将栅格图像分割为多个小块中提到的代码思路比较类似,但是在文件读取、裁剪参数设置等方面有所不同;大家如果有需要,可以先查阅以上这两篇文章。

  本文所用代码如下。

python 复制代码
# -*- coding: utf-8 -*-
"""
Created on Mon Aug 28 22:39:52 2023

@author: fkxxgis
"""

import os
import arcpy

tif_file_path = r"E:/02_Project/202307_NDVIProduce/Beijing_Preprocessing/Original"
result_file_path = r"E:/02_Project/202307_NDVIProduce/Beijing_Preprocessing/Four_Result"
arcpy.env.parallelProcessingFactor = 0

for root, dirs, files in os.walk(tif_file_path):
    for dir_name in dirs:
        dir_path = os.path.join(root, dir_name)
        arcpy.env.workspace = dir_path
        tif_file_list = arcpy.ListRasters("*", "tif")

        for tif_file in tif_file_list:
            arcpy.SplitRaster_management(tif_file,
                                         result_file_path,
                                         tif_file.split(".tif")[0] + "_",
                                         "NUMBER_OF_TILES",
                                         "TIFF",
                                         "BILINEAR",
                                         "2 2",
                                         "#",
                                         "#",
                                         "PIXELS",
                                         "#",
                                         "#",
                                         "#",
                                         "#",
                                         "#",
                                         "0")
        print(dir_path)

  首先,我们需要设置输入和输出文件夹路径;其中,tif_file_path表示存放原始遥感影像的大文件夹 路径,而result_file_path则表示存放拆分后结果的文件夹路径。随后,这里需要设置一下ArcPy的环境参数,通过arcpy.env.parallelProcessingFactor = 0这句代码禁用并行处理,以确保在处理期间不会产生问题------关于这一设置的原理,大家参考文章# ArcMap分割栅格Split Raster后得不到结果的解决方法即可。

  随后,我们使用os.walk,遍历源文件夹中的所有子文件夹 ;对于每个子文件夹 ,设置ArcPy工作环境为该子文件夹 的路径,从而使得我们可以使用 arcpy.ListRasters获取这一个子文件夹 中的所有.tif格式的遥感影像文件。

  接下来,对于每个遥感影像,使用arcpy.SplitRaster_management函数对其进行切分。其中,第一个参数tif_file是要拆分的遥感影像的路径,第二个参数result_file_path是拆分后结果保存的文件夹路径,第三个参数tif_file.split(".tif")[0] + "_"是输出文件的前缀,这里使用原始文件名去掉.tif后缀,并在末尾添加下划线;接下来的参数用于设置拆分方式、输出格式、插值方法等。这里我们选择将遥感影像拆分为2 * 2块(大家如果需要将遥感影像拆分为其他的数量,修改这里就可以),插值方法为BILINEAR,输出格式为TIFF。最后一个参数0表示将原本遥感影像中的0值作为切割后小图像的NoData值。

  运行上述代码,我们即可在结果文件夹中,看到每一景遥感影像都已经被切分为了4份,且这里每一份的文件后都带有一个数字后缀(数字后缀从0开始计算,本文中就是0、1、2与3);如下图所示。

  接下来,如果我们希望将这些得到的小遥感影像 按照文件名的特征,复制、剪切到其他文件夹中,可以参考文章# Python将多个文件夹内的文件按照文件名特征复制到不同目标路径,以及# Python提取指定名称的文件并批量复制到其他文件夹中提到的代码加以自动实现。

  至此,大功告成。

相关推荐
Jo乔戈里10 小时前
免费本地搜图软件
图像处理·python·搜索引擎·ai
for_ever_love__10 小时前
字符串处理——f-string、切片与正则,清洗文本的第一把刀
python·大模型·虚拟环境
夜晚回家10 小时前
Python 零基础入门 | 输入的数字为什么不能直接算
开发语言·python
打工仔折腾 AI10 小时前
把AI Agent托管在家用电脑:UU远程终端与端口映射实测记录
人工智能·后端·python·langchain·ai agent 实战
yivifu10 小时前
中文古籍电子书注释集成
前端·javascript·python·beautifulsoup·epub
零基础12310 小时前
Agent 的 Memory 怎么做科研:以中医诊断场景为例
人工智能·经验分享·python·语言模型
vx_Biye_Design11 小时前
springboot高校选课系统82776-计算机课程设计、毕业设计
java·vue.js·spring boot·后端·spring·课程设计·express
databook11 小时前
面向数据工程师的正则表达式:从日志清洗到字段提取
python·正则表达式·数据分析
ss27312 小时前
AI全栈实战 | 3.2-01 Python 基础:四大数据容器怎么选,推导式为什么是 Pythonic 的灵魂
开发语言·人工智能·python