【GDAL-Python】6-在Python中使用GDAL合并栅格数据

最新推荐文章于 2025-03-29 20:29:56 发布

代码搬运工的逆袭

最新推荐文章于 2025-03-29 20:29:56 发布

阅读量1.6k

点赞数 23

分类专栏： GDAL教程文章标签： python 卫星影像三维重建 gdal gdal-python

本文链接：https://blog.csdn.net/weixin_44702962/article/details/137863624

版权

本文介绍了如何使用Python的GDAL库，通过创建虚拟栅格vrt和gdal_translate命令，将同一目录下的多幅栅格切片合并成单个GeoTIFF文件，同时讨论了合并要点和代码实现细节。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

1-介绍

（1）教程内容：介绍如何使用 Python 中的 GDAL 将同一目录中的一堆栅格切片合并到单个 geotiff 文件中；

（2）技术路线：使用GDAL合并栅格数据有两个选项：
1）在Python脚本中使用子进程调用gdal_merge.py （终端也可以）
2）创建一个虚拟栅格vrt，然后使用translate转换为geotiff
PS:在合并影像文件较多且数据量大时推荐2），1）方法会加载所有文件，对硬件的内存占用较高。

1.1 虚拟栅格vrt介绍

虚拟栅格是一种将一堆现有数据集组织到一个文件中的方法，也就是告诉影像文件所在位置，以merge.vrt文件为例，内存2.48k，里面内容如下所示

<VRTDataset rasterXSize="5024" rasterYSize="4851">
  <SRS dataAxisToSRSAxisMapping="2,1">GEOGCS["WGS 84",DATUM["WGS_1984",SPHEROID["WGS 84",6378137,298.257223563,AUTHORITY["EPSG","7030"]],AUTHORITY["EPSG","6326"]],PRIMEM["Greenwich",0],UNIT["degree",0.0174532925199433,AUTHORITY["EPSG","9122"]],AXIS["Latitude",NORTH],AXIS["Longitude",EAST],AUTHORITY["EPSG","4326"]]</SRS>
  <GeoTransform> -8.1622514831999993e+01,  3.3507924922972912e-06,  0.0000000000000000e+00,  3.0274203637999999e+01,  0.0000000000000000e+00, -3.3507923440937610e-06</GeoTransform>
  <VRTRasterBand dataType="UInt16" band="1">
    <NoDataValue>65535</NoDataValue>
    <ColorInterp

最低0.47元/天解锁文章