### 如何使用 Python 将多个 NC 文件批量转换为 TIFF 文件 以下是实现这一功能的一种方法，主要依赖于 `netCDF4` 和 `rasterio` 库来处理 `.nc` 文件并将其保存为 `.tiff` 格式的栅格数据。 #### 安装必要的库 在运行代码之前，请确保安装了所需的第三方库： ```bash pip install netCDF4 rasterio numpy ``` #### 示例代码 以下是一个完整的脚本示例： ```python import os from netCDF4 import Dataset import rasterio from rasterio.transform import from_origin import numpy as np def nc_to_tif(nc_file, tif_folder, variable_name="data", time_index=0): """ 将单个 .nc 文件中的指定变量提取出来并保存为 GeoTIFF。 参数: nc_file (str): 输入的 NetCDF (.nc) 文件路径。 tif_folder (str): 输出的 GeoTIFF 文件夹路径。 variable_name (str): 需要提取的变量名称，默认为 'data'。 time_index (int): 时间维度索引，默认为第一个时间片。 """ # 加载NetCDF文件 dataset = Dataset(nc_file) # 获取变量数组 var_array = dataset.variables[variable_name][time_index] # 获取地理坐标信息 lon = dataset.variables['lon'][:] lat = dataset.variables['lat'][:] # 计算像素大小和仿射变换矩阵 pixel_size_x = abs(lon[1] - lon[0]) pixel_size_y = abs(lat[1] - lat[0]) transform = from_origin( min(lon), max(lat), pixel_size_x, pixel_size_y ) # 创建GeoTIFF文件名 base_name = os.path.basename(nc_file).split('.')[0] output_path = os.path.join(tif_folder, f"{base_name}_{time_index}.tif") # 使用RasterIO写入GeoTIFF with rasterio.open( output_path, 'w', driver='GTiff', height=var_array.shape[0], width=var_array.shape[1], count=1, dtype=str(var_array.dtype), crs='+proj=latlong', # 假设WGS84投影 transform=transform, ) as dst: dst.write(var_array.astype(rasterio.float32), indexes=1) def batch_nc_to_tif(input_folder, output_folder, variable_name="data"): """ 批量将输入文件夹中的所有.nc文件转换为.tiff文件 参数: input_folder (str): 存放 .nc 文件的文件夹路径。 output_folder (str): 输出 .tiff 文件的目标文件夹路径。 variable_name (str): 提取的变量名称，默认为 'data'。 """ if not os.path.exists(output_folder): os.makedirs(output_folder) for file in os.listdir(input_folder): if file.endswith(".nc"): nc_filepath = os.path.join(input_folder, file) nc_to_tif(nc_filepath, output_folder, variable_name=variable_name) # 设置参数 input_dir = "/path/to/your/input/nc_files" output_dir = "/path/to/your/output/tif_files" batch_nc_to_tif(input_dir, output_dir, variable_name="temperature") # 替换为实际变量名 ``` 上述代码实现了两个核心函数： - **`nc_to_tif`**: 处理单一 `.nc` 文件并将其中的一个时间切片的数据存储到 `.tiff` 中。 - **`batch_nc_to_tif`**: 对整个目录下的所有 `.nc` 文件执行批量化操作[^1]。 #### 关键点说明 1. 函数支持自定义变量名 (`variable_name`) 和时间索引 (`time_index`) 的设置，以便灵活选择所需的时间片段或特定变量。 2. 地理空间信息通过 `rasterio.transform.from_origin()` 构建，基于纬度(`lat`)和经度(`lon`)计算得到像元尺寸以及原点位置。 3. 默认假设投影为 WGS84（即经纬度），如果原始 `.nc` 数据具有其他 CRS 投影，则需调整 `crs` 参数以匹配实际情况。 #### 注意事项 - 如果目标 `.nc` 文件包含多维复杂结构或者缺失值标记，请提前确认其内部布局，并适当修改读取逻辑。 - ArcGIS 工具可能用于辅助验证最终生成的 `.tiff` 是否符合预期的空间范围与分辨率要求。