15.2. ラスタ解析

ラスタ計算機
- 例
ラスタを揃える

15.2.1. ラスタ計算機 

ラスタメニューのラスタ計算機は、既存のラスタのピクセル値に基づいた計算を実行します（図 15.17 参照）。結果はGDALがサポートする形式の新規ラスタレイヤとして書き出されます。

../../../_images/raster_calculator1.png — 図 15.17 ラスタ計算機

バンド のリストには、使用可能なロードされたラスタレイヤがすべて含まれています。ラスタ計算機の式フィールドにラスタを追加するには、フィールドリスト内で名前をダブルクリックします。そこに演算子ボタンを使用したりボックス内にキーボードで演算子を入力して、計算式を作成します。

出力レイヤ セクションでは、出力レイヤを定義する必要があります。入力ラスタレイヤに基づいて、またはX,Y座標とカラム数と行数に基づいて、計算領域の範囲を定義し出力レイヤの解像度を設定します。出力レイヤの解像度が入力レイヤと異なる場合には、最近傍アルゴリズムを使用して値がリサンプリングされます。

演算子 セクションには利用可能なすべての演算子があります。演算子をラスタ計算機の式ボックスに追加するには、それに応じたボタンをクリックします。算術計算（ + 、 - 、 * など）や三角関数（ sin 、 cos 、 tan など）が利用できます。条件式（ = 、 != 、 < 、 >= など）は、偽の場合は0、真の場合は1を返すので、他の演算子や関数と共に使うことができます。

結果をプロジェクトに追加するチェックボックスにチェックを入れると、結果レイヤは自動的に凡例エリアに追加されて表示されます。

ヒント

ラスタ計算機のアルゴリズムについても参照してください。

15.2.1.1. 例 

標高値をメートルからフィートに変換する

メートル単位の標高ラスタからフィート単位のラスタを作成するには、メートルからフィートへの換算係数 3.28 を掛けます。式は次のとおりです：

"elevation@1" * 3.28

マスクを使用する

例えば、標高0メートル以上の部分だけに興味があるなど、ラスタの一部をマスクしたい場合には、マスクの作成とその結果のラスタへの適用が以下の式を使用して一度に行えます。

("elevation@1" >= 0) * "elevation@1"

これはつまり、標高値が0以上の各セルは条件式の評価が1になり、元の値に1を掛けるので値が維持されます。そうでない場合には条件式の評価が0になり、ラスタ値が0になります。これにより、標高値に応じたマスクが作成されます。

ラスタを例えば2つの標高クラスに分類したい場合には、次の式を使用すると、1つのステップで2種類の値 1と2を持つラスタを作成することができます。

("elevation@1" < 50) * 1 + ("elevation@1" >= 50) * 2

これはつまり、標高値が50未満の各セルには値1が設定され、50以上の各セルには値2が設定されます。

15.2.2. ラスタを揃える 

このツールは複数のラスタを入力として、それらの位置を完全に合わせることが可能です。これはつまり、

同じCRSに再投影し、
同じセルサイズや同じグリッドオフセットにリサンプリングし、
関心のある領域で切り抜きを行い、
必要に応じて値を再スケールします。

ラスタはすべて、別のファイルに保存されます。

まず、ラスタ ► ラスタを揃える... からツールを開き、 ^{新しいラスタを追加} ボタンをクリックして、QGISの既存のラスタを1つ選択します。揃えたあとのラスタを保存する出力ファイルやリサンプリング方法を選択し、ツールで必要ならばセルサイズに応じてスケールにチェックを入れます。リサンプリング方法は以下が選択できます（図 15.18 参照）：

最近傍（Nearest Neighbor）
バイリニア（2x2）
キュービック（4x4 ） ：三次の畳込み内挿近似
キュービックBスプライン（4x4 ） ：キュービックBスプライン近似
ランチョス（6x6） ：ランチョス窓によるsinc補間
平均（Average） ：NODATAでない寄与ピクセルすべての平均を計算
モード ：すべてのサンプリングポイントの中で最も頻繁に出現する値を選択
最大値 、 最小値 、 中央値 、 第１四分位（Q1） 、 第３四分位（Q3） ：それぞれNODATAでない寄与ピクセルすべての最大値、最小値、中央値、第１四分位値、第３四分位値を選択