### Python 批量转换 NC 文件到 TIFF 格式的实现 为了完成从 NetCDF (NC) 到 GeoTIFF 的批量转换，可以采用 `netCDF4` 和 `rasterio` 库来读取和写入地理空间数据。下面是一个完整的解决方案[^1]： ```python import netCDF4 as nc import numpy as np import rasterio from rasterio.transform import from_origin import os def read_nc_file(file_path): dataset = nc.Dataset(file_path, 'r') lats = dataset.variables['lat'][:] lons = dataset.variables['lon'][:] var_data = dataset.variables[list(dataset.variables.keys())[-1]][:] return lats, lons, var_data def write_tif(output_filename, data_array, latitudes, longitudes): transform = from_origin(longitudes[0], latitudes[0], abs(longitudes[1]-longitudes[0]), abs(latitudes[1]-latitudes[0])) with rasterio.open( output_filename, 'w', driver='GTiff', height=data_array.shape[1], width=data_array.shape[2], count=1, dtype=str(data_array.dtype), crs='+proj=latlong', transform=transform, ) as dst: dst.write(data_array.astype(rasterio.float32), 1) def batch_convert_nc_to_tif(input_dir, output_dir): 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) lats, lons, var_data = read_nc_file(file_path) time_steps = var_data.shape[0] base_name = os.path.splitext(filename)[0] for t in range(time_steps): out_tif = f"{output_dir}/{base_name}_time_{t}.tif" single_time_slice = var_data[t,:,:].reshape(1,var_data.shape[1],var_data.shape[2]) write_tif(out_tif, single_time_slice, lats, lons) input_directory = "path/to/your/nc/files" # 替换为实际路径 output_directory = "desired/output/path" # 替换为期望保存位置 batch_convert_nc_to_tif(input_directory, output_directory) ``` 此脚本定义了三个主要函数：用于读取 `.nc` 文件的数据；创建并保存对应的 `.tif` 文件；以及遍历指定目录中的所有 `.nc` 文件执行上述两个过程。