r - 更改 netcdf 的投影

Question

我有每日温度和降水数据的 netcdf 文件如何更改 netcdf 的投影？

我曾尝试使用 raster 函数来读取文件并使用 projectraster 重新投影它们？

代码_使用

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”的继承方法另外：警告消息：在min（dim（to）[1：2]）中： min 没有非缺失参数；返回 Inf

score 3 · Accepted Answer

的帮助projectRaster给出了这种用法：

Usage:

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

如果您使用两个未命名的参数调用它，它们将与fromand匹配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包中具有相同名称的函数使用 .

1 回答 1