如何使用GDAL获取NetCDF与GeoTIFF文件中的波段（变量）名称及存储顺序

我完全理解你想精简依赖、只用GDAL解决这个问题的需求——确实，GDAL的Dataset.GetMetaData()确实不直接返回波段/变量名称列表，不过我们可以通过GDAL的其他原生API绕开这个限制，针对NetCDF和GeoTIFF分别处理：

首先要纠正你当前代码里的一个小问题：你现在打开的是NetCDF里的单个变量（路径里加了:var），这样只能操作这个单独的变量，没法获取整个文件里的所有变量。我们需要先打开整个数据集，再提取变量/波段信息。

1. 处理NetCDF文件 #

NetCDF在GDAL里是以**子数据集（SubDatasets）**的形式存储变量的，我们可以用Dataset.GetSubDatasets()来获取所有变量的信息：

from osgeo import gdal

def get_netcdf_variables(ds: gdal.Dataset) -> dict:
    """从NetCDF数据集提取变量名称及存储顺序"""
    variables = {}
    sub_datasets = ds.GetSubDatasets()
    # 遍历子数据集，按存储顺序索引（从0开始）
    for idx, (sub_path, sub_desc) in enumerate(sub_datasets):
        # 从子数据集路径中提取变量名，格式为NETCDF:"path.nc":varname
        var_name = sub_path.split(':')[-1].strip('"')
        variables[idx] = var_name
    return variables

# 示例用法：打开整个NetCDF文件
input_dir = "/your/input/dir"
file_name = "your_file"
ds = gdal.Open(f"{input_dir}/{file_name}.nc")
if ds:
    nc_vars = get_netcdf_variables(ds)
    print("NetCDF变量及顺序：")
    for order, name in nc_vars.items():
        print(f"顺序 {order}: {name}")
    ds = None  # 释放资源

2. 处理GeoTIFF文件 #

GeoTIFF的波段名称通常存储在每个波段的元数据或描述中，我们可以遍历每个波段，用GetDescription()或读取波段元数据的NAME字段来获取：

def get_geotiff_band_names(ds: gdal.Dataset) -> dict:
    """从GeoTIFF数据集提取波段名称及存储顺序"""
    bands = {}
    band_count = ds.RasterCount
    # 遍历波段，GDAL波段索引从1开始
    for band_idx in range(1, band_count + 1):
        band = ds.GetRasterBand(band_idx)
        # 优先尝试从元数据获取NAME字段
        band_meta = band.GetMetadata()
        band_name = band_meta.get("NAME")
        # 如果元数据里没有，用波段描述；都没有则用默认命名
        if not band_name:
            band_name = band.GetDescription() or f"Band_{band_idx}"
        # 存储顺序保留GDAL的1-based索引，也可以改成0-based
        bands[band_idx] = band_name
        band = None  # 释放资源
    return bands

# 示例用法：打开GeoTIFF文件
tiff_ds = gdal.Open("/your/path/your_file.tif")
if tiff_ds:
    tiff_bands = get_geotiff_band_names(tiff_ds)
    print("GeoTIFF波段及顺序：")
    for order, name in tiff_bands.items():
        print(f"顺序 {order}: {name}")
    tiff_ds = None

通用封装函数 #

可以把两种格式的处理封装成一个函数，自动判断文件类型：

def get_dataset_variables(ds: gdal.Dataset) -> dict:
    """自动判断数据集类型，返回变量/波段名称及存储顺序"""
    driver_name = ds.GetDriver().ShortName
    if driver_name == "netCDF":
        return get_netcdf_variables(ds)
    elif driver_name == "GTiff":
        return get_geotiff_band_names(ds)
    else:
        # 其他格式默认返回波段索引和默认名称，可按需扩展
        bands = {}
        for idx in range(1, ds.RasterCount +1):
            bands[idx] = f"Band_{idx}"
        return bands

# 通用用法
ds = gdal.Open("/your/file/path")
if ds:
    vars_info = get_dataset_variables(ds)
    print("变量/波段信息：")
    for order, name in vars_info.items():
        print(f"存储顺序 {order}: {name}")
    ds = None

这个方案完全基于GDAL原生API，不需要引入xarray、rasterio等额外依赖，完美符合你精简脚本的需求。

内容的提问来源于stack exchange，提问作者Philipe Riskalla Leal