python批量解析CLDAS-V2.0的.nc数据文件

原创

已于 2023-04-18 18:23:00 修改 · 5.5k 阅读

50 ·

CC 4.0 BY-SA版权

文章标签：

#NetCDF #气象数据 #python

于 2020-12-16 15:23:00 首次发布

这篇博客介绍了如何使用Python的NetCDF4和gdal库将CLDAS-V2.0的.nc数据文件（以四川省为例）转换为栅格图像。首先通过ArcMap的多维工具预览数据，然后通过编写Python脚本读取.nc文件，利用NetCDF4获取变量，并用gdal创建栅格数据，设置地理信息和投影，最后将数据写入.tif文件。脚本适用于批量处理多个.nc文件，并可在ArcMap中正常显示结果。

该文章已生成可运行项目，

处理CLDAS-V2.0的.nc数据文件——四川省为例

CLDAS-V2.0是中国气象局陆面数据同化系统实时产品数据集，可以从中国气象数据网下载在这里插入图片描述
以四川省为例，下载了一段时间的数据，数据格式为.nc格式（NETCDF DATA），如下

可以使用ArcMap的多维工具的“创建NetCDF栅格图层”的工具来打开，叠加一个全国的省界，可以看到这个数据是四川的某个数据
在这里插入图片描述
现在使用NetCDF4和gdal来处理这个文件夹中的459个.nc文件，转为栅格图像

使用NetCDF4读取一个nc文件，可以看到是这样的

其中variables是它的值，有经纬度（坐标）和一个TAIR的值（灰度）
下面pip好NetCDF4和gdal两个包，开始写一个处理的py文件

import netCDF4
from osgeo import gdal
import numpy as np
import os

# 读取一个nc文件转换为栅格
def read_nc_to_tif(nc_path, out_path):
    # 读取nc数据保存为字典
    dataset = netCDF4.Dataset(nc_path)  # 打开一个nc文件
    keys = dataset.variables.keys()  # 获取它的值的列表
    dicts = {
   
   }
    for key in keys:
        dicts[key] = np.array(</