相关疑难解决方法(0)

我有一个栅格时间序列存储在多个GeoTIFF文件（*.tif）中，我想将其转换为单个NetCDF文件。数据是uint16。

我可能会使用以下gdal_translate命令将每个图像转换为netcdf：

 gdal_translate -of netcdf -co FORMAT=NC4 20150520_0164.tif foo.nc

然后使用一些脚本NCO从文件名中提取日期然后进行连接，但是我想知道是否可以在Python中使用xarray它更有效地进行此操作，并且它是新的rasterio后端。

我可以轻松读取文件：

import glob
import xarray as xr
f = glob.glob('*.tif')
da = xr.open_rasterio(f[0]) 
da

哪个返回

<xarray.DataArray (band: 1, y: 5490, x: 5490)>
[30140100 values with dtype=uint16]
Coordinates:
  * band     (band) int64 1
  * y        (y) float64 5e+05 5e+05 5e+05 5e+05 5e+05 4.999e+05 4.999e+05 ...
  * x        (x) float64 8e+05 8e+05 8e+05 8e+05 8.001e+05 8.001e+05 …