如何在Python中使用gdal_translate将netCDF转为COG格式

一、subprocess调用的错误修正 #

你的脚本存在几个关键问题，导致COG文件无法生成：

1. 子数据集名称不匹配
   OSGeo4W命令中使用的子数据集是monthly_prcp_norm，但脚本里写的是mlyprcp_norm，名称不一致会导致GDAL找不到对应变量，直接执行失败。

2. subprocess命令格式问题
   你将整个命令放在列表的单个元素中，虽shell=True支持这种写法，但更规范且不易出错的方式是拆分参数为列表元素；同时你未等待进程执行完成，也未捕获错误输出，无法得知命令执行失败的原因。

3. 冗余的-of参数
   原命令同时写了-of GTiff和-of COG，GDAL中COG是基于GTiff的特殊格式，仅需-of COG即可，重复参数可能引发未知问题。

修正后的subprocess代码：

import os
import subprocess

input_folder_path = r"C:\MyTest\data"
output_folder_path = r"C:\MyTest\output"
gdal_translate = r"C:\Users\my_user\.conda\envs\geospatial\Library\bin\gdal_translate.exe"

test_input_nc = os.path.join(input_folder_path, 'my-precip-data.nc')
test_output_tif = os.path.join(output_folder_path, 'my_precip-monthly-normal.tif')

# 拆分参数为列表（推荐，shell=False更安全）
command = [
    gdal_translate,
    f"NETCDF:{test_input_nc}:monthly_prcp_norm",
    test_output_tif,
    "-of", "COG",
    "-co", "COMPRESS=LZW"
]

# 等待进程完成并捕获输出
process = subprocess.Popen(command, stdout=subprocess.PIPE, stderr=subprocess.PIPE, text=True)
stdout, stderr = process.communicate()

# 检查执行状态
if process.returncode != 0:
    print(f"执行错误: {stderr}")
else:
    print(f"执行成功: {stdout}")

二、使用gdal.translate()直接转COG #

无需调用subprocess，直接用GDAL Python绑定的gdal.Translate()方法即可，核心是指定format='COG'，并设置创建选项，同时正确指定netCDF的子数据集路径：

from osgeo import gdal
import os

input_folder_path = r"C:\MyTest\data"
output_folder_path = r"C:\MyTest\output"

test_input_nc = os.path.join(input_folder_path, 'my-precip-data.nc')
# 构造netCDF子数据集路径
input_ds_path = f"NETCDF:{test_input_nc}:monthly_prcp_norm"
test_output_tif = os.path.join(output_folder_path, 'my_precip-monthly-normal.tif')

# 设置COG的创建选项
creation_options = [
    'COMPRESS=LZW',
    # 可选：添加其他COG相关参数，比如'PREDICTOR=2'、'BIGTIFF=IF_SAFE'等
]

# 调用Translate方法
ds = gdal.Translate(
    destName=test_output_tif,
    srcDS=input_ds_path,
    format='COG',
    creationOptions=creation_options
)

# 释放数据集资源
if ds:
    ds = None
    print("COG文件生成成功")
else:
    print("生成失败，请检查输入路径和子数据集名称")

关键说明： #

• 子数据集路径格式必须为NETCDF:{nc文件路径}:{变量名}，和命令行格式一致。
• format='COG'告诉GDAL直接生成云优化GeoTIFF，无需额外设置format='GTiff'。
• creationOptions对应命令行的-co参数，可根据需求添加压缩、分块等设置。

内容的提问来源于stack exchange，提问作者gwydion93