Question

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

我尝试使用光栅功能来执行此操作以读取文件并用projectraster重新投影？

已使用代码

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”的继承方法另外：警告消息：在min（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。