Question

我对这段代码大汗淋漓。我之前收到过帮助来构建它here。简而言之，我在这里做的是我有三个栅格的列表，我随机抽样numberv次。因此，输出是四个列表的列表，每个列表有三个栅格。获取随机点位置后，我会在此位置获取栅格值。我想解决的问题是我想采用第二个样本位置，即sample.set[[1]][2]并从rasters[1]获取栅格值。然后我想取sample.set[[1]][3]并从rasters[2]获取栅格值。然后sample.set[[2]][2]并从rasters[1]和sample.set[[2]][3]获取栅格值并从rasters[2]获取栅格值等。结果将是4个列表的列表，每个列表包含2个元素样本xy值（位置）和先前的栅格值。非常感谢帮助。

y <- matrix(1:150,50,3)
mv <- c(1,2,3)
rep = 20

valuematrix <- vector("list",ncol(y))

for (i in 1:ncol(y)) {
        newmatrix <- replicate(rep,y[,i])
        valuematrix[[i]] <- newmatrix
}

library(sp)
library(raster)

rasters <- setNames(lapply(valuematrix, function(x) raster(x)), 
                    paste0('raster',1:length(mv)))

# Create a loop that will sample the rasters
library(dismo)

numberv = c(10,12,14,16)        # sample number vector 

# Function to sample using a given number (returns list of three)
sample.number <- function(x) {
        rps <- lapply(rasters, function(y) randomPoints(raster(y),n=x))
        setNames(rps,paste0('sample',1:length(mv)))
}

# Apply sample.number() to your numberv list
sample.set <- lapply(numberv,sample.number)

# Function to extract values from a given sample
sample.extract <- function(x) {
        lapply(1:length(x),function(y) data.frame(x[[y]],
                                                  extract(rasters[[y]],x[[y]])))
}

# Apply sample.extract() to the set of samples (returns list of four lists)
sample.values <- lapply(sample.set,sample.extract)

现在我想使用列表sample.set的第二个元素中的样本值来对列表rasters中的第一个栅格进行采样我尝试了这个但没有成功：

sample.extract.prev <- function(x) {
        lapply(1:length(x),function(y) data.frame(x[[y]],
                                                  extract(rasters[[y]],x+1[[y]])))
}


sample.values.prev <- lapply(sample.set,sample.extract.prev)

Answer 1

管理解决这个问题（对我来说是高五;）不幸的是，我设法用一个循环来做，很高兴看到一个函数的例子。

samplevaluesnext <- vector("list",length(sample.set))

## Look up values
for (j in 1:length(sample.set)) {
        for (i in 1:(length(rasters)-1)) {
                samplevaluesnext[[j]][[i]] <- data.frame(sample.set[[j]][[i+1]],
                                                         extract(rasters[[i]], 
                                                                 as.data.frame(sample.set[[j]][i+1])))
        }
}

使用R中的函数对栅格和位置列表进行采样

1 个答案: