群集

时间:2016-08-19 13:05:09

标签: r raster r-raster

    #load a .nc file
    temp=brick(temp.nc)
temp
temp
class       : RasterBrick
dimensions  : 180, 360, 64800, 2928  (nrow, ncol, ncell, nlayers)
resolution  : 1, 1  (x, y)
extent      : 0, 360, -90, 90  (xmin, xmax, ymin, ymax)
coord. ref. : +proj=longlat +datum=WGS84
data source : /fhgfs/data/work/crct/bfb16057/pgf/1.0deg/3hourly/tas_3hourly_1948-1948.nc
names       : X1948.01.01.00.00.00, X1948.01.01.03.00.00, X1948.01.01.06.00.00, X1948.01.01.09.00.00, X1948.01.01.12.00.00, X1948.01.01.15.00.00, X1948.01.01.18.00.00, X1948.01.01.21.00.00, X1948.01.02.00.00.00, X1948.01.02.03.00.00, X1948.01.02.06.00.00, X1948.01.02.09.00.00, X1948.01.02.12.00.00, X1948.01.02.15.00.00, X1948.01.02.18.00.00, ...
Date/time   : 1948-01-01 00:00:00, 1948-12-31 21:00:00 (min, max)
varname     : tas
level       : 1

我正在使用包栅格在群集上处理大型netcdf文件。首先我遇到了一些自定义函数的问题,所以我不得不重写它们并使用getValues和setValues来对这些rasterBricks执行计算,从而解决了这个问题。

但现在我在主脚本中调用reclassify

reclasstemp=reclassify(temp,rcl)

R崩溃了。如果我在本地机器上运行它,一切都有效。是否有解决方法或解决方案?

  

矩阵错误(unlist(ini),ncol = 2,byrow = TRUE):'数据'必须   是矢量类型,是' NULL'电话:重新分类......   .rasterObjectFromFile - > .rasterFromRasterFile - > readIniFile - >   矩阵另外:有50个或更多警告(使用警告()   见第50页

我无法使用get / setValues重新分类,因为重新分类需要一个栅格,因此不会起作用。否则,我将不得不在v=getValues(temp)

中重新分类矩阵
> str(v)
 num [1:64800, 1:2928] 222 222 222 222 221 ...
 - attr(*, "dimnames")=List of 2
  ..$ : NULL
  ..$ : chr [1:2928] "X1948.01.01.00.00.00" "X1948.01.01.03.00.00" "X1948.01.01.06.00.00" "X1948.01.01.09.00.00" ...

0 个答案:

没有答案