Python实现栅格数据拼接、裁剪等操作

从地理空间数据云（http://www.gscloud.cn/）上下载了GDEMV2—30米分辨率的数据，检索按照行政区来的。

下载完的数据解压缩以后，如下：

1.首先把所有的文件夹下的dem数据路径读入到TXT文件中，代码：

此部分参考：https://blog.csdn.net/weixin_42236288/article/details/81810761

import os# 查询符合条件的文件
f = open(r'D:\ProfessionalProfile\DEMdata\StatisticsDEMfile.txt','w')
j,k=0,0
def file(root, ext):for i in os.listdir(root):  # os.listdir(root) 读取root目录下的文件和目录# os.path.splitext()将文件名和扩展名分开ext = os.path.splitext(root + '\\' + i)  # os.path.splitext()读取文件扩展名if ext[1]:  # 有扩展名时 为文件global jif i.count(ex) > 0:  # 包含特定字符pathDEM = root +'\\'+ if.write(pathDEM + '\n')f.flush()# print(pathDEM)  # 输出j += 1else:global kif i.count(ex) > 0:  # 包含特定字符pathDEM2 = root + r'\\' + i# print( pathDEM2)  # 输出k += 1root1 = root + '\\' + i  # 拼接目录file(root1, ext)  # 递归调用print('文件数量为：', j)# print('文件夹数量为：', k)
# f.close()# root = os.getcwd()  # 生成目录 root
# 包含的特定字符
ex = 'dem'
# 调用函数
file('D:\ProfessionalProfile\DEMdata\GDEMV2_30m_shandong','.tif')

运行后得到文件：

2.对该TXT文件里的影像进行拼接

参考：https://blog.csdn.net/qq_43177210/article/details/108402353?utm_medium=distribute.pc_relevant.none-task-blog-BlogCommendFromMachineLearnPai2-3.control&dist_request_id=2a4eb29a-3fae-44ad-8384-ce28a6a41f0a&depth_1-utm_source=distribute.pc_relevant.none-task-blog-BlogCommendFromMachineLearnPai2-3.control

# 如果图片大于4G，需要在输入图片路径前加入 -co BIGTIFF=YES
# subprocess 模块首先推荐使用的是它的 run 方法，更高级的用法可以直接使用 Popen 接口。
import subprocess# 路径
create_slope = r'''E:\r_python\python.exe E:\r_python\Scripts/gdal_merge.py -of GTiff -o '''
list = []
# 打开TXT文件路径
f = open(r'D:\ProfessionalProfile\DEMdata\StatisticsDEMfile.txt', 'r')
# 获取TXT文件中所有DEMTIF文件路径
# readlines()函数：1、一次性读取整个文件。2、自动将文件内容分析成一个行的列表。
lines = f.readlines()
# strip()函数：用于移除字符串头尾指定的字符（默认为空格或换行符）或字符序列。# 注意：该方法只能删除开头或是结尾的字符，不能删除中间部分的字符。
fileDir = [x.strip() for x in lines]
# 遍历所有TIF文件
list = fileDir
# print(list)# 生成文件路径
tarDir = r'D:\ProfessionalProfile\DEMdata'
#这句没懂
filename = ' '.join(list)
# l3 = len(filename)
# print(filename)
# subprocess.call执行指定的命令，返回命令执行状态，其功能类似于os.system(cmd)。详见：https://www.cnblogs.com/zhou2019/p/10582716.html
# subprocess.call(args, *, stdin=None, stdout=None, stderr=None, shell=False, timeout=None)
subprocess.call(create_slope + r'D:\ProfessionalProfile\DEMdata\demMosaic1224.tif' + ' ' + '-co COMPRESS=LZW ' + filename)

运行结果如下（ENVI中打开）：

此外，还学习到了很多，如：

subprocess模块的信息可参考：https://www.cnblogs.com/zhou2019/p/10582716.html

https://blog.csdn.net/m0_38051293/article/details/103218755?utm_medium=distribute.pc_relevant.none-task-blog-baidujs_title-7&spm=1001.2101.3001.4242

https://blog.csdn.net/weixin_40450867/article/details/104191939

https://blog.csdn.net/qq_36178899/article/details/83989031?utm_medium=distribute.pc_relevant_download.none-task-blog-baidujs-1.nonecase&depth_1-utm_source=distribute.pc_relevant_download.none-task-blog-baidujs-1.nonecase

3.裁剪

https://blog.csdn.net/weixin_40501429/article/details/114140680