使用python GDAL生成COG(Cloud Optimized GeoTIFF)

原创

已于 2023-08-15 17:23:13 修改 · 3.1k 阅读

9 ·

CC 4.0 BY-SA版权

文章标签：

#python

于 2022-03-24 17:23:40 首次发布

本文介绍了GDAL 3.1版本新增的COG（Cloud-Optimized GeoTIFF）栅格驱动，使得在warp和translate中直接输出COG格式成为可能，简化了之前的操作流程。通过gdalwarp和gdal_translate命令示例展示了如何生成COG文件，并提供了相应的Python函数实现。内容包括设置nodata值、压缩选项以及重采样方法等，帮助用户高效地转换和处理GeoTIFF数据。

更新

gdal 3.1版本更新了COG栅格driver，现在可以在warp或translate里直接用-of COG输出了，之前那一大堆都不需要了
官网的示例：

gdalwarp src1.tif src2.tif out.tif -of COG
gdal_translate world.tif world_webmerc_cog.tif -of COG -co TILING_SCHEME=GoogleMapsCompatible -co COMPRESS=JPEG

官网文档：https://gdal.org/drivers/raster/cog.html#raster-cog

参考资料：

https://trac.osgeo.org/gdal/wiki/CloudOptimizedGeoTIFF#HowtogenerateitwithGDAL
https://gdal.org/programs/gdal_translate.html
https://gdal.org/drivers/raster/cog.html
https://geoexamples.com/other/2019/02/08/cog-tutorial.html/

几个主要函数

def translateToCOG(in_ds, out_path):
    """
    将dataset转为cog文件
    :param in_ds: 输入dataset
    :param out_path: 输出路径
    :return:
    """
    im_bands = in_ds.RasterCount
    for i in range(im_bands):
        # 获取nodata和波段统计值
        nodataVal = in_ds.GetRasterBand(i + 1).GetNoDataValue()
        maxBandValue =