### 将多个NetCDF文件批量转换为GeoTIFF 为了实现将多个NetCDF文件批量转换为GeoTIFF的功能，可以利用Python中的`osgeo.gdal`库以及`netCDF4`库来处理这些操作。下面提供一段完整的代码示例用于说明如何完成此任务。 #### 安装必要的依赖包 在执行具体的操作之前，确保已经安装了所需的Python库： ```bash pip install gdal netCDF4 numpy pandas ``` #### 批量转换脚本 这段代码会遍历指定目录下的所有`.nc`文件，并将其逐一转换成对应的`.tif`文件存储到目标路径下。 ```python import os from netCDF4 import Dataset import numpy as np from osgeo import gdal, osr def get_netcdf_info(nc_file_path): dataset = Dataset(nc_file_path) # 获取地理坐标系信息和其他元数据 spatial_ref = osr.SpatialReference() spatial_ref.ImportFromEPSG(4326) # 假设使用WGS84投影 return { 'dataset': dataset, 'spatial_reference': spatial_ref.ExportToWkt(), 'bands_count': len(dataset.variables), 'band_names': list(dataset.variables.keys()) } def write_geotiff(output_filename, array_data, geo_transform, projection): driver = gdal.GetDriverByName('GTiff') out_raster = driver.Create( output_filename, array_data.shape[1], array_data.shape[0], 1, gdal.GDT_Float32) out_raster.SetGeoTransform(geo_transform) out_band = out_raster.GetRasterBand(1) out_band.WriteArray(array_data) out_raster_srs = osr.SpatialReference(wkt=projection) out_raster.SetProjection(out_raster_srs.ExportToWkt()) out_band.FlushCache() input_dir = './nc_files/' # 输入NetCDF文件夹位置 output_dir = './tifs/' # 输出GeoTIFF文件夹位置 if not os.path.exists(output_dir): os.makedirs(output_dir) for filename in os.listdir(input_dir): if filename.endswith(".nc"): file_path = os.path.join(input_dir, filename) info_dict = get_netcdf_info(file_path) ds = info_dict['dataset'] for band_name in info_dict['band_names']: try: var = ds.variables[band_name][:] # 这里假设所有的变量都是二维数组；如果不是，则需要额外判断逻辑 # 创建输出文件名 base_name = os.path.splitext(os.path.basename(filename))[0] tif_output_path = f"{output_dir}{base_name}_{band_name}.tif" # 设置仿射变换参数（这里仅作为例子给出固定值） geotransform = (-180.0, 0.5, 0.0, 90.0, 0.0, -0.5) write_geotiff(tif_output_path, var.data, geotransform, info_dict['spatial_reference']) except Exception as e: print(f"Error processing {filename} with variable {band_name}: ", str(e)) ``` 上述代码实现了从NetCDF文件读取数据并保存为带有空间参考系统的GeoTIFF文件的过程[^2]。