我正在尝试为某些站点数据分配值，这些数据恰好位于我拥有天气数据的区域之外。我正在尝试根据最近的单元格值提取，如果可能的话，提取 40 公里内的单元格值。
我的光栅 (r) 如下所示：
class(r)
class : RasterBrick
dimensions : 201, 464, 93264, 23376 (nrow, ncol, ncell, nlayers)
resolution : 0.25, 0.25 (x, y)
extent : -40.5, 75.5, 25.25, 75.5 (xmin, xmax, ymin, ymax)
coord. ref. : +proj=longlat +datum=WGS84
data source : \\ueahome\eressci5\zuw13bqu\data\NTProfile\Desktop\EOBS European data\rr_0.25deg_reg_v10.0.nc
names : X1950.01.01, X1950.01.02, X1950.01.03, X1950.01.04, X1950.01.05, X1950.01.06, X1950.01.07, X1950.01.08, X1950.01.09, X1950.01.10, X1950.01.11, X1950.01.12, X1950.01.13, X1950.01.14, X1950.01.15, ...
Date : 1950-01-01, 2013-12-31 (min, max)
varname : rr
我使用以下代码根据纬度和经度数据提取
vals <- extract(r, matrix(c(issues[22,3], issues[22,2]), …我有一个包含年份的向量 A：例如。
A <- c(1978, 1979, 1980, 1981, 1984, 1987,1988,1989,1990,1991, 1992)
我希望能够在错过一年之前计算年份的顺序。所以我会寻找我的答案：4, 1, 6
我知道rle用于重复数字的序列，但不确定在这里做什么。