В GDAL пришел к этому:
Код: Выделить всё
mask=gdal.Open(path1,gdal.GA_ReadOnly)
src=gdal.Open(path2,gdal.GA_ReadOnly)
dst=gdal.GetDriverByName('NetCDF').Create(pathdst,xsize=mask.RasterXSize,ysize=mask.RasterYSize)
dst.SetGeoTransform(mask.GetGeoTransform())
dst.SetProjection(mask.GetProjection())
optt=gdal.WarpOptions(format='NetCDF',xRes=mask.RasterXSize,yRes=mask.RasterYSize)
dst=gdal.Warp(dst,src,options=optt)
dst=None
src=None
mask=None
Добился результата с помощью Reprojectimage, но там обрабатывается только 1 слой, а я хочу как-то для всех сразу сделать warp. С Translate та же история (это логично).
Итак, вопрос: как с помощью Python и GDAL сделать перепроецирование растра в сетку другого растра для всех каналов (слоев времени в данном случае) сразу?
Если кто-то найдет в интернете ссылку по этому вопросу, которую я ещё не видел - тому мой огромный респект. Но я пересмотрел уже все что только можно, в том числе документацию. Оно не работает.