Изменение проекции netcdf - PullRequest
Новая
       29

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

0 голосов
/ 19 сентября 2018

У меня есть файлы 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)

ERROR_

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

1 Ответ

0 голосов
/ 20 сентября 2018

Справка для 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.

...