这里提供一个基于GDAL库的Python代码示例,可以将netCDF格式的降水数据文件转换为栅格数据格式:

import gdal

# 输入文件路径
in_file = r"C:\Users\y\Desktop\pr_Amon_BCC-CSM2-MR_ssp585_r1i1p1f1_gn_201501-210012.nc"
# 输出栅格文件路径
out_file = r"C:\Users\y\Desktop\pr.tif"

# 打开netCDF文件
ds = gdal.Open("NETCDF:{}:pr".format(in_file))
# 获取数据集的投影信息和地理变换信息
proj = ds.GetProjection()
trans = ds.GetGeoTransform()

# 获取数据集的第一个波段
band = ds.GetRasterBand(1)
# 将波段数据写入栅格文件
gdal.Translate(out_file, band, outputType=gdal.GDT_Float32, format="GTiff", outputSRS=proj, outputBounds=trans)

执行此代码后,将生成一个名为“pr.tif”的栅格文件,包含了输入netCDF文件中的降水数据。如果需要处理多个变量或时间步长的数据,可以在循环中逐个读取和转换多个文件

用python将这个路径:CUsersyDesktoppr_Amon_BCC-CSM2-MR_ssp585_r1i1p1f1_gn_201501-210012nc的降水文件进行NC转栅格处理

原文地址: https://www.cveoy.top/t/topic/ccRo 著作权归作者所有。请勿转载和采集!

免费AI点我,无需注册和登录