解像度と範囲が異なる2つのラスターがあります。
> res(Elevation)
[1] 0.002083333 0.002083333
> res(Ann_precip)
[1] 0.008333333 0.008333333
> extent(Elevation)
class : Extent
xmin : -15.07722
xmax : -7.641806
ymin : 7.193611
ymax : 12.67694
> extent(Ann_precip)
class : Extent
xmin : -15.075
xmax : -7.641667
ymin : 7.191667
ymax : 12.675
私の質問は、これらの2つのラスターが一致する解像度と範囲を持つためには、
A)raster::aggregate
関数を使用する
> 0.008333333/0.002083333
[1] 4
Elevation_res<-aggregate(Elevation, fact=4, fun=mean)
そしてraster::extend
機能
Elevation_res<-extend(Elevation_res, Ann_precip, values=NA)
(ただし、ここでも異なるが非常に似た範囲と解像度が得られます):
> res(Elevation_res)
[1] 0.008333333 0.008333333
> res(Ann_precip)
[1] 0.008333333 0.008333333
> res(Elevation_res)==res(Ann_precip)
[1] FALSE FALSE
> extent(Elevation_res)
class : Extent
xmin : -15.07722
xmax : -7.635556
ymin : 7.193611
ymax : 12.67694
> extent(Ann_precip)
class : Extent
xmin : -15.075
xmax : -7.641667
ymin : 7.191667
ymax : 12.675
または
b)raster::resample
関数を使用する
Elevation_res<-resample(Elevation, Ann_precip, method="bilinear")
> res(Elevation_res)==res(Ann_precip)
[1] TRUE TRUE
> extent(Elevation_res)==extent(Ann_precip)
[1] TRUE
私が読んだので、私はこれを聞いてるのよWegmann(2016)ら大幅にリサンプリングするピクセル値に影響を与え、そのこと(P110)(私が正しく理解している場合)aggregate()
、extend()
およびcrop()
代わりに使用する必要があります。私の場合、解像度と範囲の違いは非常に小さいので、リサンプリングによって生じるバイアスはここでは最小限になると思いますか?