画素単位乾き/湿り周波数

Question

約60年間、standardized rainfallのmonthly時系列データとrasterbrickaを考えると、どのように私はrasterbrickで次のように計算することができます。

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 Iはを好むであろう範囲内の値−1.99 to −1.50

を含みますは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