Изменение проекции netcdf

Question

У меня есть файлы netcdf с данными о дневной температуре и осадках. Как можно изменить проекцию netcdf?

Я пытался сделать это с помощью растровой функции, чтобы прочитать файлы и перепроектировать их с помощью Projectraster?

Code_used

a <-raster(file.nc)
cr1<-"+proj=longlat +datum=NAD83 +no_defs+ellps=GRS80 +towgs84=0,0,0"
 projectRaster(a, cr1,res =  0.04166667)

ОШИБКА_

Ошибка в (функция (классы, fdef, mtable): невозможно найти унаследованный метод для функции 'res' для подписи '"символ"'. Дополнительно: предупреждающее сообщение: в мин. (Dim(to)[1:2]): нет не пропущенных аргументов для min; возвращает Inf

Answer 1

Помощь для projectRaster дает это использование:

Usage:

     projectRaster(from, to, res, crs, method="bilinear", 
                  alignOnly=FALSE, over=FALSE, filename="", ...) 

если вы называете это с двумя безымянными аргументами, они соответствуют from а также to, а также to является:

  to: Raster* object with the parameters to which 'from' should be
      projected

если вместо этого вы называете второй аргумент, это должно сработать:

  a_project = projectRaster(a, crs = cr1, res =  0.04166667)

потому что теперь второй аргумент сопоставляется как crs аргумент.

НО, если вы просто хотите преобразовать gdal-совместимый набор данных, вы можете использовать gdaltransformлибо из командной строки, либо через функцию с тем же именем из gdalUtils пакет.