Python:如何使用GDAL实现多波段遥感影像的批量拼接?

本打算自己写一个拼接函数,但是感觉太折腾了,查了一下gdal的实现方法,发现有封装的非常好的一个方法(gdal.Warp)可以解决,直接就拿来用了。

一行代码就可以解决,需要注意的

gdal.Warp(out_mosaic_path, [img_path_1, img_path_2], srcNodata=0.0, dstNodata=0.0)

out_mosaic_path:拼接好的影像的输出路径

[img_path_1, img_path_2, ······]:列表形式,传入所有需要进行拼接的tiff文件的路径;例如此处需要拼接的tiff文件有两个。

srcNodata:原始数据集(待拼接影像)的缺失值数值。(例如你的遥感影像中0.0为缺失值,那么就传入0.0)

dstNodata:设置拼接好的遥感影像的缺失值(这里就还是照旧0.0吧,如果说你觉得0.0可能是有意义的,那么你可以设为例如-9999等)

至于批量处理,例如我需要相同日期的影像进行拼接,那么可以如下操作:

in_dir = r'J:\Objects\TobaccoLeafRecognition\Data\VV_VH\TIFF\in_TIFF'
out_dir = r'J:\Objects\TobaccoLeafRecognition\Data\VV_VH\TIFF\in_TIFF'
for i in range(1, 13):
    date_str = 'S1A_IW_GRDH_1SDV_2019_{:02}*.tiff'.format(i)  # 每个月
    mosaic_paths = glob(os.path.join(in_dir, date_str))
    gdal.Warp(os.path.join(out_dir, 'S1_{:02}.tiff'.format(i)), mosaic_paths, srcNodata=0.0, dstNodata=0.0)

相关推荐

最近更新

  1. TCP协议是安全的吗?

    2024-03-17 13:38:04       18 阅读
  2. 阿里云服务器执行yum,一直下载docker-ce-stable失败

    2024-03-17 13:38:04       19 阅读
  3. 【Python教程】压缩PDF文件大小

    2024-03-17 13:38:04       19 阅读
  4. 通过文章id递归查询所有评论(xml)

    2024-03-17 13:38:04       20 阅读

热门阅读

  1. 基础练习题之函数

    2024-03-17 13:38:04       20 阅读
  2. ClickHouse副本节点数据损坏恢复

    2024-03-17 13:38:04       23 阅读
  3. Clickhouse MergeTree原理(二)—— 表和分区的维护

    2024-03-17 13:38:04       18 阅读
  4. Centos设置docker自启动,以及容器程序自启动

    2024-03-17 13:38:04       21 阅读
  5. Python:递归函数

    2024-03-17 13:38:04       22 阅读
  6. html导航栏+下拉菜单+表单验证

    2024-03-17 13:38:04       23 阅读
  7. HTML

    HTML

    2024-03-17 13:38:04      20 阅读
  8. 在CentOS 7系统下通过二进制方式安装MySQL 8.0.34

    2024-03-17 13:38:04       21 阅读
  9. Jtti:如何在CentOS中安装和配置Tomcat应用服务器

    2024-03-17 13:38:04       21 阅读