xarray nc转tif

import xarray as xr

import rioxarray  # conda install -c conda-forge rioxarray

da = xr.open_Dataset(paht_file)  # 最好将原始dataset改为dataarray, 如下

ds = xr.Dataarray(da.prec.values, coords=[ds.lat, ds.lon], dims=['y', 'x'])  #  这里的维度名称最好用x ,y表示,如果输出tif文件,但在rio中识别不了其他的维度名称

ds.rio.write_crs("epsg:4236"), inplace=True) # 第一个参数是投影信息,第二个是将投影信息的名称等替换或者写入dataarray中

ds.rio.to_raster(out_path + '****.tif') # 文件可以被arcgis打开查看属性信息,不知道为啥看图软件打不开

如果要批量 就写个循环就行,根据时间 文件数量都行。也可以定义为一个函数。

更多信息 可以去rioxarray官网查看,该模块可以用来处理raster,https://corteva.github.io/rioxarray/stable/getting_started/manage_information_loss.html

 

 

你可能感兴趣的:(python,python)