更改 netcdf 的投影

【问题标题】：Changing projection of a netcdf更改 netcdf 的投影
【发布时间】：2019-02-24 02:43:05
【问题描述】：

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

我曾尝试使用光栅函数来读取文件并用投影仪重新投影？

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

【问题讨论】：

【解决方案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 包中的同名函数。

【讨论】：