4

ラスターの解像度と範囲を変更したいと考えています。私はさまざまな方法で結合しようとしましたが、うまくいきませんでした... 以下は、適切な解像度と範囲を取得するために書いたコードですが、結果として、合計ピクセル値が変更されます (そして、差はかなり大きくなる可能性があります。ある場合...)。ラスターの範囲と解像度に関するオンラインの投稿にはかなりの数がありますが、私は自分の問題を解決できませんでした。総ピクセル値 (cellStats(r, sum)) の変化を最小限に抑えながら、解像度と範囲を変更する方法はありますか? どうもありがとうございました。

library(raster)

#sample rasters
r <- raster(extent(-180,179.9999, -55.91668 , 83.58326))
res(r)<-c(0.5/6, 0.5/6)
r <- setValues(r, runif(ncell(r)))
cellStats(r, sum) #3615229

temp_extent <- raster(extent(-180,180, -90 , 90))
res(temp_extent)<-c(0.5, 0.5)

# to get resolution (0.5, 0.5)
r1 <- aggregate(r, fact=6, fun=sum)
cellStats(r1, sum) #3615229
r1
# to get extent (-180,180,-90,90)
r2 <- resample(x= r1, y=temp_extent, method="bilinear") # what I am hoping is for NA to fill the added pixels
cellStats(r2, sum) #problem: gives different total value: 3615916
r2
4

1 に答える 1