R には、同じ経度と緯度の座標を持つ 216 のグリッドまたはセルを持つ 2 つのラスター オブジェクトがあります。最初のラスターには 183 のレイヤーがあり、以下の情報で表されます。
> CMAWRF[[1]]
class : SpatRaster
dimensions : 12, 18, 183 (nrow, ncol, nlyr)
resolution : 0.25, 0.25 (x, y)
extent : 3.125, 7.625, 50.625, 53.625 (xmin, xmax, ymin, ymax)
coord. ref. : lon/lat WGS 84
source : IPSL-IPSL-CM5A-MR_r1i1p1_IPSL-WRF381P1980.nc
varname : sfcWindAdjust (Bias-Adjust Near-Surface Wind Speed)
names : sfcWi~ust_1, sfcWi~ust_2, sfcWi~ust_3, sfcWi~ust_4, sfcWi~ust_5, sfcWi~ust_6, ...
unit : m s-1, m s-1, m s-1, m s-1, m s-1, m s-1, ...
time : 1980-01-01 10:30:00 to 1980-12-31 10:30:00
2 番目のラスターにはレイヤーが 1 つしかなく、範囲も同じです。お気づきかもしれませんが、緯度座標は 3.125 から 7.625 で、経度座標は 50.625 から 53.625 です。これらの座標は、オランダをカプセル化します。
したがって、私の目標は、人口ラスターの各 I,J (行の場合は I、列の場合は J) セル値を、各レイヤーの風速ラスターの各 I,J セル値で乗算することです。ただし、オランダの国境内にあるセル値に対してのみこれを実行したいので、空間ラスターでカバーされる地域全体ではありません。
Rでこれを行う方法を知っている人はいますか?