単純な i、j、k 位置で識別される 3D グリッドにいくつかのデータがあります (実際の空間情報はありません)。これらのデータは現在 RasterStack にあります。
b <- stack(system.file("external/rlogo.grd", package="raster"))
# add more layers
b <- stack(b,b)
# dimensions
dim(b)
[1] 77 101 6
77 行、101 列、6 層になります。
# upscale by 2
up <- aggregate(b,fact=2)
dim(up)
[1] 39 51 6
39 行、51 列、6 層になります。
期待される動作: 3 つのレイヤー。
各レイヤー内で集約するという現在の動作に加えて、レイヤー間で集約する方法を探しています。私は他のデータ構造にもオープンですが、自分で書いたものよりも既存のアップスケーリング/リサンプリング/集計アルゴリズムを好みます。
関連する可能性があるのは、http: //quantitative-advice.gg.mq.edu.au/t/fast-way-to-grid-and-sum-coordinates/110/5または時空間パッケージです。より複雑になります。