asc转tif python,NC转TIF数据处理（python）

在进行气候分析时，很多人都会用到ERA5数据。关于ERA5数据的下载，很多人都给出了相应的下载方法，官方也有相应的处理代码，这里简单说一下怎么把NC数据转成GIS分析常用的栅格TIF数据。代码如下

import numpy as np

import netCDF4 as nc

from osgeo import gdal,osr,ogr

import os

import glob

def NC_to_tiffs(data,Output_folder):

nc_data_obj = nc.Dataset(data)

Lon = nc_data_obj.variables['longitude'][:]

Lat = nc_data_obj.variables['latitude'][:]

tp_arr = np.asarray(nc_data_obj.variables['tp'])  #将tp数据读取为数组

tp_arr_float = tp_arr.astype(float)*1000

LonMin,LatMax,LonMax,LatMin = [Lon.min(),Lat.max(),Lon.max(),Lat.min()]

N_Lat = len(Lat)

N_Lon = len(Lon)

Lon_Res = (LonMax - LonMin) /(float(N_Lon)-1)

Lat_Res = (LatMax - LatMin) / (float(N_Lat)-1)

for i in range(len(tp_arr[:])):

#创建.tif文件

driver =