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开始计算,本文中就是0123);如下图所示。

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

至此,大功告成。

相关推荐
tan180°6 小时前
MySQL表的操作(3)
linux·数据库·c++·vscode·后端·mysql
优创学社27 小时前
基于springboot的社区生鲜团购系统
java·spring boot·后端
why技术7 小时前
Stack Overflow,轰然倒下!
前端·人工智能·后端
幽络源小助理7 小时前
SpringBoot基于Mysql的商业辅助决策系统设计与实现
java·vue.js·spring boot·后端·mysql·spring
烛阴7 小时前
简单入门Python装饰器
前端·python
好开心啊没烦恼8 小时前
Python 数据分析:numpy,说人话,说说数组维度。听故事学知识点怎么这么容易?
开发语言·人工智能·python·数据挖掘·数据分析·numpy
面朝大海,春不暖,花不开8 小时前
使用 Python 实现 ETL 流程:从文本文件提取到数据处理的全面指南
python·etl·原型模式
ai小鬼头8 小时前
AIStarter如何助力用户与创作者?Stable Diffusion一键管理教程!
后端·架构·github
简佐义的博客8 小时前
破解非模式物种GO/KEGG注释难题
开发语言·数据库·后端·oracle·golang