Python 将 NC 格式降水数据转换为栅格数据的代码

以下是将 NC 格式降水数据转栅格的 Python 代码，需要使用 netCDF4 和 numpy 库：

import netCDF4 as nc
import numpy as np
import gdal

# 打开 NC 文件
nc_file = nc.Dataset('path/to/nc/file.nc')

# 获取数据
data = nc_file.variables['precipitation'][:]

# 获取经纬度坐标
lon = nc_file.variables['lon'][:]
lat = nc_file.variables['lat'][:]

# 定义输出文件名和路径
output_file = 'path/to/output/raster.tif'

# 定义栅格的元数据
geotransform = (lon[0], lon[1]-lon[0], 0, lat[-1], 0, lat[-1]-lat[-2])
projection = 'WGS84'

# 创建栅格文件
driver = gdal.GetDriverByName('GTiff')
raster = driver.Create(output_file, len(lon), len(lat), 1, gdal.GDT_Float32)
raster.SetGeoTransform(geotransform)
raster.SetProjection(projection)

# 写入数据到栅格文件
band = raster.GetRasterBand(1)
band.WriteArray(np.flipud(data))
band.FlushCache()

注意：

1. 在运行代码之前，需要先安装 netCDF4、numpy 和 gdal 库。
2. 代码中的路径和文件名需要根据实际情况进行修改。