より複雑なオブジェクト
raster パッケージは RasterLayer オブジェクトを提供しますが、さらに複雑なオブジェクトとして RasterStack と RasterBrick もあります。これら2つのオブジェクトは、同じ範囲と次元を持つ多数のラスタ(いわゆる「マルチバンド」または「マルチレイヤー」ラスタ)を保存するために設計されています。
RasterLayer は行列のように考えられますが、RasterStack と RasterBrick は三次元配列に近いイメージです。これらを扱うためにもう1つ知っておくべきことは、特定のレイヤーをどのように指定するかです。
RasterStack や RasterBrick に対して $ または [[ でサブセット化すると、1つのレイヤーを取り出して新しい RasterLayer オブジェクトを返せます。たとえば、x が RasterStack の場合、x$layer_name や x[["layer_name"]] は、layer_name という名前のレイヤーだけを含む RasterLayer を返します。
ここでは、pop と同じ範囲をカバーし、人口をいくつかの年齢層に分けたレイヤーを含む pop_by_age という RasterStack オブジェクトを見ていきます。
この演習はコースの一部です
R で地理空間データを可視化する
演習の手順
pop_by_ageを表示します。すべてのレイヤー名が見えますか?[[サブセットを使ってunder_1レイヤーを抽出します。- 前の指示で書いたコードをそのまま
plot()に渡して、under_1レイヤーをプロットします。
実践的なインタラクティブ演習
このサンプルコードを完成させて、この演習に挑戦してみましょう。
# Print pop_by_age
# Subset out the under_1 layer using [[
# Plot the under_1 layer