Question

    #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）：＆＃39;数据＆＃39;必须是矢量类型，是＆＃39; NULL＆＃39;电话：重新分类...... .rasterObjectFromFile - ＆gt; .rasterFromRasterFile - ＆gt; readIniFile - ＆gt; 矩阵另外：有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 个答案: