逐像素乾燥/溼潤頻率RasterBrickř

Question

鑑於standardized rainfall一個rasterbricka與monthly時間序列數據爲大約60年，我怎樣才能計算上rasterbrick以下，而不必訴諸於數據塌陷成dataframe：

Q = negative/positive number of severely dry/wet months per total months of severely dry/wet events × 100 

其中severely wet包括在所述範圍內的值1.50 to 1.99和severely dry包括在所述範圍內的值−1.99 to −1.50

我更喜歡Q爲raster格式。下載數據集並執行：

dd=spei03_df 
dd[1:2]<-dd[2:1]#swap lat and lon 
a=rasterFromXYZ(dd) 

有什麼想法？謝謝。

Example data set can be found here

Answer 1

您可以使用光柵代數從raster包。這裏是一個計算嚴重潮溼的例子。最終輸出a_wet_perc是一個柵格圖層，單元值範圍爲0到100（以百分比表示）。您可以按照相同的邏輯創建總結嚴重乾燥條件的柵格圖層。

library(raster) 

# Calculate the Severe Wet 
a_wet <- a 

# Replace vlaues >= 1.5 and <= 1.99 to 100 
# (Does not need to be 100. A very large number should be good) 
a_wet[a_wet >= 1.5 & a_wet <= 1.99] <- 100 
# Replace other values to 0 
a_wet[a_wet < 100] <- 0 
# Replace 100 to 1 
a_wet[a_wet == 100] <- 1 

# Calculate the percentage of severe wet 
a_wet_perc <- mean(a_wet) * 100