在柵格列表上執行循環

Question

需要解決方案，非常感謝您的幫助。

在下面的代碼中，我創建了三個柵格。然後我在這個光柵上創建一個點位置的隨機number，我正在接收一個三個矩陣的列表，這些矩陣的座標是那些稱爲samples的隨機位置。然後我將這些位置和樣本柵格值接收到samplevalues。

我想改變的是我想創建一組100,150,200和250個隨機點位置（numberv）。因此，在生成這些位置並接收到位置列表之後，每個柵格將被採樣length(numberv)次（在這種情況下爲4次）。因爲我有三個柵格，那麼我想要獲得一個列表，其中包含從我的三個柵格中抽取100個樣本值的第一個元素，第二個元素抽樣的柵格值爲150次，等等。該列表將包含length(numberv)個元素。然後，我將使用這些位置來獲取這些位置的柵格值。

我只用一個樣本粘貼了一個簡單案例的乾淨代碼（使用了1個元素矢量number），希望它有幫助。

y <- matrix(1:300,100,3) 
mv <- c(1,2,3) 
rep = 200 

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) 

number = 100      # current number for random sample points number 
numberv = c(100,150,200,250)  # sample number vector i want to use 

# samples below will hold only coordinate values: 
samples <- setNames(lapply(rasters, function(x) randomPoints(raster(x), 
                  n=number)), 
        paste0('sample',1:length(mv))) 

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

for (i in 1:ncol(y)) { 
     samplevalues[[i]] <- data.frame(samples[[i]],extract(rasters[[i]], 
                  samples[[i]])) 
} 

Answer 1

這是行嗎？

# 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) 

# Access sample 1 of number 200 
summary(sample.values[[3]][[1]])