重要

翻訳はあなたが参加できるコミュニティの取り組みです。このページは現在 98.88% 翻訳されています。

22. GRASS GIS の統合

GRASS統合は、GRASS GISデータベースと機能へのアクセスを提供します（文献とWeb参照中のGRASS-PROJECTを参照）。統合は、プロバイダとプラグインの二つの部分から構成されています。プロバイダでは、GRASSラスタおよびベクタレイヤを閲覧、管理、視覚化できます。プラグインを使用すると、新しいGRASS locationとmapsetを作成し、GRASS領域を変更し、ベクタレイヤを作成・編集し、400の以上のGRASSモジュールでGRASS 2-Dと3-Dデータを分析できます。このセクションではプロバイダやプラグイン機能を導入し、GRASSデータを管理したりそれで作業するいくつかの例をあげます。

プロバイダはGRASSバージョン6および7をサポートし、プラグインはGRASS 6および7をサポートしています（QGIS 2.12以降）。QGISの配布には、GRASS 6またはGRASS 7、若しくは同時に両方のバージョン（バイナリは異なるファイル名を持ちます）のプロバイダ/プラグインを含んでいる場合があります。ただし、プロバイダ/プラグインのバージョンは実行時に1つだけ読み込めます。

22.1. デモデータセット 

As an example, we will use the QGIS Alaska dataset (see section サンプルデータのダウンロード). It includes a small sample GRASS LOCATION with three vector layers and one raster elevation map. Create a new folder called grassdata, download the QGIS 'Alaska' dataset qgis_sample_data.zip from https://download.qgis.org//downloads/data/ and unzip the file into grassdata.

More sample GRASS LOCATIONs are available at the GRASS website.

22.2. GRASSラスタおよびベクタレイヤを読み込む 

プロバイダが QGIS に読み込まれている場合、GRASS アイコンの付いたロケーションアイテムが、GRASS ロケーションを含む各フォルダアイテムの下のブラウザツリーに追加されます。フォルダ grassdata に移動し、ロケーション alaska とマップセット demo を展開します。

他のレイヤーと同様に、ブラウザからレイヤー項目をダブルクリックするか、マップキャンバスや凡例にドラッグ＆ドロップするかによって、任意のGRASSラスタレイヤとベクタレイヤを読み込めます。

Tip

GRASSデータの読み込み

GRASSロケーションアイテムが表示されない場合は、GRASSベクタプロバイダが読み込まれているかどうか、ヘルプ ► QGISについて ► プロバイダを検証してください。

22.3. ドラッグ＆ドロップでGRASSロケーションへデータをインポートする 

このセクションでは、GRASS mapsetにラスタとベクタデータをインポートする方法を例示します。

QGISのブラウザで、データをインポートしたいマップセットに移動します。
QGISのブラウザで、GRASSにインポートするレイヤを探します。ソースデータがツリー内のマップセットから離れすぎている場合は、ブラウザの別のインスタンス（ブラウザパネル（2））を開くことができることを覚えておいてください。
レイヤをドラッグし、インポート先のマップセットにドロップします。大きなレイヤのインポートには時間がかかることがあり、インポートが完了するまで、新たなレイヤ項目の前にアニメーションアイコンが表示されます。

ラスタデータが異なるCRSにある場合、それらは Approximate （高速）または Exact （正確な）変換を使用して再投影できます。元ラスタへのリンクが作成された場合（ r.external を使用して）、元データは同じCRSにあり、フォーマットがGDALにはよく知られている、元データのCRSが使用されます。 GRASSオプションのブラウザタブでこれらのオプションを設定できます。

元ラスタがより多くのバンドを持っている場合、 .<band number> 接尾辞を持つ新しいGRASS地図が各レイヤについて作成され、すべてのマップのグループアイコンが作成されます。外部ラスタは別のアイコンを持ちます。

22.4. QGISブラウザでGRASSデータを管理する 

マップをコピーする：GRASSマップは、同じロケーション内のマップセット間でドラッグアンドドロップを使ってコピーできます。
マップを削除する：GRASSマップ上で右クリックしてコンテキストメニューから削除を選択します。
マップの名前を変更する：GRASSマップ上で右クリックしてコンテキストメニューから名前を変更を選択します。

22.5. GRASSオプション 

GRASSオプションは GRASS オプションダイアログ中で設定できます。このダイアログはブラウザでロケーションまたはマップセットの項目を右クリックした後、 GRASS オプションを選択すると開きます。

22.6. GRASSプラグインを起動する 

QGISでGRASSの機能を使用するには、プラグインマネージャを使ってGRASSプラグインを選択し、読み込む必要があります。これをするには、プラグイン ► プラグインの管理とインストール... メニューに行き、 GRASS を選択し、 :guilabel:`OK`をクリックします。

GRASSプラグインを起動すると、以下の主な機能がGRASSメニュー（プラグイン ► GRASS ）で提供されています：

^{Mapsetを開く}
^新規Mapset
^{Mapsetを閉じる}
^{GRASS ツールを開く}
^{現在のGRASS領域を表示}
^{GRASSオプション}

22.7. GRASS mapsetを開く 

プラグインでGRASSツールへのアクセスを取得するためにGRASS mapsetを開く必要があります（mapsetが何も開かれていない場合はツールが無効になっています）。mapsetはブラウザから開くことができます：mapset項目を右クリックして、 Mapsetを開くをコンテキストメニューから選択します。

22.8. GRASS[場所]と[地図セット]

GRASSのデータはGISDBASEと呼ばれるディレクトリに保存されます。このディレクトリは grassdata と呼ばれることが多く、QGISでGRASSプラグインを使用する前に作成する必要があります。このディレクトリの中で、GRASS GISデータは LOCATIONs と呼ばれるサブディレクトリに格納されているプロジェクトごとに整理されています。各 LOCATION は座標系、地図投影法、地理的境界によって定義されます。各 LOCATION は複数の MAPSETs (LOCATION のサブディレクトリ) を持つことができ、それらはプロジェクトを異なるトピックやサブ領域に細分化したり、個々のチームメンバーのワークスペースとして使用されます (ref:literature_and_web の Neteler & Mitasova 2008 を参照)。GRASSモジュールでベクタやラスタのレイヤを分析するには、一般的にGRASSの LOCATION にインポートする必要があります。(これは厳密には正しくありません。GRASSモジュールの r.external と v.external を使えば、インポートせずに外部のGDALサポートデータセットへの読み取り専用リンクを作成することができます。これは初心者がGRASSを使うときの一般的な方法ではないので、ここでは説明しません)。

../../../_images/grass_location.png — 図 22.1 アラスカLOCATIONのGRASSデータ

22.9. GRASS LOCATIONへデータをインポートする 

データを簡単にブラウザにドラッグ＆ドロップしてインポートできる方法を見つけるために、ドラッグ＆ドロップでGRASSロケーションへデータをインポートするセクションを参照してください。

このセクションでは、標準GRASSモジュールを使用して、伝統的な方法でQGIS「アラスカ」データセットによって提供される「アラスカ」GRASS LOCATION に、ラスタおよびベクタデータをインポートする方法の例を示します。したがって、QGIS「アラスカ」データセットから土地被覆ラスター地図 landcover.img とベクターGMLファイル lakes.gml を使用します（サンプルデータのダウンロードを参照）。

QGISを起動し、GRASSプラグインがロードされていることを確認します。
GRASSツールバーで、 ^{Mapsetを開く} アイコンをクリックして、 Mapset ウィザードを起動します。
GRASSデータベースとして、QGISアラスカデータセットの grassdata フォルダを、 ロケーション として 'alaska'、 マップセット として 'demo' を選択し、 OK をクリックします。
ここで ^{GRASSツールを開く} アイコンをクリックします。GRASSツールボックスダイアログが表示されます（セクション GRASSツールボックスを参照）。
ラスター地図 landcover.img をインポートするには、モジュールタブの r.in.gdal モジュールをクリックして下さい。この GRASSモジュールは GDAL がサポートしているラスタファイルをGRASS LOCATION にインポートします。 r.in.gdal モジュールダイアログが表示されます。
QGIS「アラスカ」データセット内の raster フォルダを参照して、ファイル landcover.img を選択します。
出力ラスタマップ名として landcover_grass を定義し、実行をクリックします。guilabel:Output`タブに、現在実行中のGRASSコマンド ``r.in.gdal -o input=/path/to/landcover.img output=landcover_grass` が表示されます。
成功しました と表示されたら、出力を見るをクリックします。これで landcover_grass ラスタレイヤがGRASSにインポートされ、QGISキャンバスに表示されます。
ベクタGMLファイル lakes.gml をインポートするにはモジュールツリータブのモジュール v.in.ogr をクリックします。このGRASSモジュールを使用すると、GRASS LOCATION にOGRサポートのベクターファイルをインポートできます。 v.in.ogr のモジュールダイアログが表示されます。
QGIS「アラスカ」データセット内の gml フォルダを参照し、OGRファイルとして lakes.gml ファイルを選択します。
ベクタ出力名として lakes_grass を定義し、実行をクリックします。この例では他のオプションを気にする必要はありません。出力タブには、現在実行中のGRASSコマンド v.in.ogr -o dsn=/path/to/lakes.gml output=lakes_grass が表示されます。
成功しました と表示されたら、出力を見るをクリックします。lakes_grass ベクタレイヤがGRASSにインポートされ、QGISキャンバスに表示されます。

22.9.1. 新しいGRASS[場所]を作成する 

例として、単位としてフィートを使い、アルバース等積投影で投影された、サンプルGRASS [場所]アラスカ です。このサンプルGRASS [場所]アラスカ は、以下のGRASS関連のセクションのすべての例と演習のために使用されます。コンピュータ上のデータセットをダウンロードしてインストールすると便利です（サンプルデータのダウンロードを参照）。

QGISを起動し、GRASSプラグインがロードされていることを確認します。
QGISアラスカデータセット（サンプルデータのダウンロード参照）からの alaska.shp シェープファイル（セクションファイルからレイヤを読み込む参照）を可視化します。
GRASSツールバーで、 ^{地図セット} アイコンをクリックして新しい地図セットウィザードを起動してください。
既存のGRASSデータベース(GISDBASE)フォルダ grassdata を選択するか、コンピュータのファイルマネージャを使って新しい LOCATION 用のフォルダを作成します。そして次へをクリックします。
このウィザードを使用して、既存の LOCATION の中に新しい MAPSET を作成したり（新しい[地図セット]を追加するのセクションを参照）、完全に新しい LOCATION を作成することができます。新規ロケーションを作成を選択します（参照図 22.2）。
LOCATION の名前 - ここでは'alaska'とした - を入力し、:guilabel:`次へ`をクリックします。
投影リストを有効にするラジオボタン投影法をクリックして投影を定義します。
私たちはアルバース等積アラスカ（フィート）投影を使用しています。それはEPSG ID 2964によって表されることを知っているので、検索ボックスに入力します。（注意：このプロセスを別の 場所 と投影法で繰り返したいが、EPSG IDを記憶していない場合は、ステータスバーの右下隅にある ^{CRS Status} アイコンをクリックしてください（投影法の操作セクションを参照））。
フィルタに 2964 と入力して投影法を選択して下さい.
次へをクリックします。
デフォルト区域を定義するために、東西南北の LOCATION の境界を入力しなければなりません。ここで、ボタン現在のQGISの範囲を設定をクリックすると、読み込まれたレイヤ alaska.shp の範囲が GRASS のデフォルト領域として適用されます。
次へをクリックします。
新しい :file:‘LOCATION` 内に MAPSET を定義する必要もあります（新しい LOCATION を作成するときはこれが必要です）。それには好きに名前を付けられます - ここでは'demo'を使いました。GRASSは PERMANENT と呼ばれる特別な MAPSET を自動的に作成します。それはプロジェクトのコアデータ、そのデフォルトの空間的範囲および座標系の定義を格納するように設計されています（文献とWeb参照中のNeteler＆Mitasova 2008を参照）。
要約が正しいことを確認し、完了をクリックします。
新しい LOCATION 「アラスカ」と２つの MAPSET 'demo'と'PERMANENT'が作られました. 現在開かれているワーキングセットはあなたが定義した 'demo' です。
GRASSツールバーのそれまでは利用できなかったいくつかのツールが利用可能になっています.

../../../_images/create_grass_location.png — 図 22.2 QGIS内で新しいGRASS[場所]または新しい[地図セット]を作成する

手順が多いように感じたかもしれませんが、実際はそれほど悪くなく、LOCATION を素早く作成することができます。これで LOCATION 'alaska' はデータをインポートする準備が整いました（GRASS LOCATIONへデータをインポートするセクションを参照）。また、QGIS 'Alaska' データセットサンプルデータのダウンロードに含まれるサンプルのGRASS LOCATION 'alaska'に既に存在するベクタデータとラスタデータを使って、セクション GRASSベクタデータモデルに進むこともできます。

22.9.2. 新しい[地図セット]を追加する 

ユーザーは、自分で作成したGRASS MAPSET への書き込みアクセス権を持ちます。これは、自分の MAPSET へのアクセス以外にも、他のユーザーの MAPSET 中の地図を（そして他のユーザーはあなたのものを）読むことができることを意味しますが、変更したり削除できるの自分の MAPSET 中の地図だけです。

すべての MAPSET には WIND ファイルが含まれ,そこには現在の領域の座標値と選択されているラスタの解像度が格納されています (Neteler & Mitasova 2008 文献とWeb参照, セクション GRASS領域ツール参照).

QGISを起動し、GRASSプラグインがロードされていることを確認します。
GRASSツールバーで、 ^{地図セット} アイコンをクリックして新しい地図セットウィザードを起動してください。
さらに「テスト」という MAPSET を追加するため、 LOCATION 「アラスカ」の GRASS データベース (GISDBASE) フォルダ grassdata を選択して下さい.
次へをクリックします。
このウィザードを使用して、既存の LOCATION 内に新しい MAPSET を作成したり、全く新しい LOCATION を作成することができます。ラジオボタン場所を選択をクリックし（図 22.2 を参照）、次へをクリックします。
新しい MAPSET に test という名前を入力します。ウィザードの下の方に、既存の MAPSET とそれに対応する所有者のリストが表示されます。
次へをクリックし、要約が正しいことを確認し、完了をクリックします。

22.10. GRASSベクタデータモデル 

デジタイズに先立って GRASSベクタデータモデルを理解することが重要です。一般的には、GRASSは、トポロジ的ベクタモデルを使用しています。これは、領域が閉じたポリゴンとしてではなく、一つ以上の境界で表現されることを意味します。二つの隣接する領域間の境界は一度だけデジタイズされ、それは、両方の領域によって共有されています。境界は隙間なく接続され閉じている必要があります。領域は、領域の重心によって識別（およびラベル付け）されます。

境界と重心のほかに、ベクタマップにはポイントとラインも含めることができます。これらのすべてのジオメトリ要素は、一つのベクタで混合できる一つのGRASSベクタマップ内の異なるいわゆる「レイヤ」で表されます。だからGRASSにおいては、レイヤはベクタまたはラスタ地図ではなく、ベクタレイヤ内部の１つのレベルです。これは慎重に区別することが重要です。（ジオメトリ要素を混合することは可能ですが、それは異常なことであり、GRASSにおいてさえベクタネットワーク解析などの特別な場合に使用されるだけです。通常は、異なるジオメトリ要素は異なるレイヤに格納することをお勧めします。）

1つのベクタデータセット内には複数の「レイヤー」を格納できます。例えば、田畑、森林、湖沼は、一つのベクタに格納できます。隣接する森林と湖沼には同じ境界を共有できますが、それらは別々の属性テーブルを持っています。境界に属性を付けることも可能です。例として、湖沼と森林の境界が道路なので、異なる属性テーブルを持つ場合があります。

地物の「レイヤ」はGRASS内部で「レイヤ」によって示されます。「レイヤ」は（そのジオメトリが森か湖かなど）データセット内に複数のレイヤがあるかどうかを示す数値です。今のところ、それは数にだけできます。将来的には、GRASSは、ユーザーインターフェイス内のフィールドとして名前をサポートします。

属性はGRASS 場所 内部にdBase、SQLite3として、または外部データベーステーブル、例えば、PostgreSQL、MySQL、Oracleなどに格納できます。

データベーステーブル中の属性は「カテゴリ」値を使用してジオメトリ要素とリンクされます。

「カテゴリ」（キーID）はジオメトリプリミティブに付随した整数であり、それはデータベーステーブル内の1つのキー列へのリンクとして使用されます。

Tip

GRASSベクタモデルについて調べる

GRASSベクタモデルとその機能を学ぶ最良の方法は、ベクタモデルがより深く説明されている多くのGRASSチュートリアルの一つをダウンロードすることです。より詳しい情報、書籍、チュートリアルは https://grass.osgeo.org/learn/manuals/ をご覧ください。

22.11. 新しいGRASSベクタレイヤを作成する 

新しいGRASSベクタレイヤを作るには、ブラウザで地図セットコンテキストメニューから次の項目のいずれかを選択します：

点レイヤの新規作成
線レイヤの新規作成
ポリゴンレイヤの新規作成

そして、ダイアログに名前を入力します。新しいベクタマップが作成され、レイヤはキャンバスに追加され、編集が開始されます。レイヤの種類を選択することは、ベクタマップにデジタイズできるジオメトリタイプを制限するものではありません。GRASSでは、1つのベクタマップですべてのジオメトリタイプ（ポイント、ライン、ポリゴン）をまとめることができます。QGISはひとつのレイヤが特定のタイプを持っている必要があるので、種類はレイヤをキャンバスに追加するのに使われるだけです。

既存のベクタマップのコンテキストメニューから、上記項目のいずれかを選択し、既存のベクタマップにレイヤを追加することもできます。

GRASSはトポロジ的ベクタモデルを使用しているので、ひとつのレイヤでジオメトリタイプの全種類（ポイント、ラインおよびエリア）をまとめることができ、新しいGRASSベクタを作成するときにジオメトリタイプを選択する必要はありません。これがQGISでのシェープファイルの作成と異なっているのは、シェープファイルが単純地物ベクタモデルを使用しているためです（新しいベクタレイヤを作成するの項を参照）。

22.12. GRASSベクタレイヤをデジタイズして編集する 

GRASSベクタレイヤは、標準のQGISデジタイズツールを使用してデジタイズできます。しかし以下のためにいくつかの特殊性があり、これらは知っておくべきでしょう：

GRASS トポロジモデル　対　QGIS単純地物
GRASS モデルの複雑さ
- 単独のマップに複数のレイヤ
- 単独のマップに複数のジオメトリタイプ
- 複数のレイヤからの複数の地物によってジオメトリを共有します

その特殊性はこの後のセクションで説明します。

保存する、変更を破棄する、元に戻す、やり直す

警告

編集中に行われたすべての変更は、すぐにベクタマップと関連属性テーブルに書き込まれます。

変更は、各操作の後に書き込まれますが、編集を終了するときに、元に戻す/やり直す、すべての変更を破棄することができます。元に戻すまたは変更を破棄する場合、元の状態がベクタマップと属性テーブルに再び書き込まれます。

この動作には主な理由が2つあります：

これは、ユーザーは自分がしていることがやりたいことであり、作業が突然中断されたとき（例えば停電）には、データが保存されている方が良い、という信念に由来するGRASSベクタの性質です。
トポロジ的なデータを効果的に編集するために必要なのは、トポロジ的に正しいかどうかの情報を視覚化することであり、そのような情報は、地図に変更が書き込まれた場合にのみ、GRASSのベクタ地図から得ることができます。

ツールバー

GRASSレイヤが編集されるときに「デジタイジングツールバー」にはいくつかの特定のツールがあります：

アイコン	ツール	目的
	点を追加	新しい点をデジタイズする
	線を追加	新しいラインをデジタイズする
	境界線を追加	新しい境界線をデジタイズする
	重心を追加	新しい重心をデジタイズする (ラベルのあるエリア)
	閉じた境界線を追加	新しい閉境界をデジタイズする

表 GRASS デジタイジング : GRASS デジタイジングツール

Tip

GRASSでポリゴンをデジタイズする

GRASSでポリゴンを作成する場合は、最初にポリゴンの境界線をデジタイズします。その後、閉じた境界線に重心（ラベルポイント）を追加します。この理由は、トポロジ的ベクタモデルはポリゴンの属性情報を、境界に対してではなく、常に重心に対してリンクしているためです。

カテゴリ

しばしばcatと呼ばれるカテゴリは、IDの一種です。その名前はGRASSベクターのみが単独で「カテゴリ」属性持っていた時代から来ています。カテゴリは、ジオメトリと属性の間のリンクとして使用されています。単一のジオメトリは、複数のカテゴリを持ち、したがって異なるレイヤに複数の地物を表すことができます。現在はQGIS編集ツールを使ってレイヤごとに1つのカテゴリだけを割り当てることができます。新しい地物は、境界線を除き、自動的に新しいユニークなカテゴリを割り当てます。境界線は通常は領域を形成し、線形の地物を表すものではありません。しかし、たとえば異なるレイヤで、後から境界線の属性を定義することができます。

新しいカテゴリは、常に現在編集中のレイヤにのみ作成されます。

QGIS編集を使ってジオメトリに複数のカテゴリを割り当てることはできません。そのようなデータは複数の地物として適切に表され、個々の地物は、異なるレイヤからであっても、削除できます。

属性

現在編集中のレイヤの属性のみが変更できます。ベクタマップにより多くのレイヤが含まれている場合、他のレイヤの地物は、その属性が編集できないことを警告するため、すべての属性が '<not editable (layer #)>' に設定されます。その理由は、QGISがレイヤ毎に１つの固定されたフィールドセットしかサポートしないのに対して、他のレイヤが、異なるフィールドセットを持っことができ、通常は持っているからです。

ジオメトリプリミティブにカテゴリが割り当てられていない場合は、新しいユニークなカテゴリが自動的に割り当てられ、そのジオメトリの属性が変更されたとき、属性テーブルに新しいレコードが作成されます。

Tip

テーブル内の属性の一括更新を、たとえば「フィールド計算機」（フィールド計算機を使用する）を使って、行いたい場合で、更新したくないカテゴリを持たない地物（通常は境界線）がある場合、「高度なフィルタ」 cat is not null を設定することによりフィルタで除外できます。

編集スタイル

トポロジ的シンボルは、トポロジ的データを効果的に編集するために不可欠です。開始を編集する場合は、専門の「GRASS編集」レンダラーをレイヤに自動的に設定され、編集が閉じているときに、元のレンダラーが復元されます。スタイルは、レイヤのプロパティ「スタイル」タブでカスタマイズできます。スタイルもプロジェクトファイルまたは任意の他のスタイルとして個別のファイルに保存できます。スタイルをカスタマイズする場合は名前を変更しないでください。名前は編集が再び開始されるときスタイルをリセットするために使用されます。

Tip

レイヤが編集されたときに、プロジェクトファイルを保存しないでください、レイヤは「編集スタイル」で格納されますが、レイヤが編集されていない場合には意味がありません。

スタイルは一時的にフィールド「topo_symbol」として属性テーブルに追加されたトポロジ情報に基づいています。フィールドは編集が閉じられたときに自動的に削除されます。

Tip

属性テーブルから「topo_symbol」フィールドを削除しないでください。レンダラーはその列に基づいているため、地物が見えなくなるでしょう。

スナップ

エリアを形成するためには、接続された境界の頂点は正確に同じ座標を持っている必要があります。これは、キャンバスとベクタマップが同じCRSを持っている場合にのみ、スナップツールを使用して達成できます。そうでない場合は、地図座標からキャンバスへのそしてその逆の正しい変換で、座標が表現誤差とCRS変換のためにわずかに異なってくることがあります。

Tip

編集するときはキャンバスについてもレイヤのCRSを使用してください。

制限

同時に同じベクタ内の複数のレイヤを同時編集することはサポートされていません。これは主に、単一のデータソースのスタックを元に戻す複数の処理が不可能であることに起因します。

LinuxおよびMac OS上でだけGRASSレイヤを一度に一つだけ編集できます。これは、ランダムな順序でデータベースドライバを閉じられないというGRASSのバグによるものです。これは、GRASSの開発者によって解決されてようとしています。

Tip

GRASS編集権限

編集を行いたい場合は GRASS MAPSET の所有者になることが必要です。所有する以外の MAPSET に属するレイヤは、ファイルに書き込み権限を持っていても編集できません。

22.13. GRASS領域ツール 

GRASSの領域定義（空間的な作業ウィンドウを設定）は、ラスタレイヤを操作するために重要です。ベクタ分析はデフォルトで任意の定義された領域の定義に限定されていません。しかし、すべての新しく作成されたラスタは、元の範囲と解像度にかかわらず、現在定義されているGRASS領域の空間的な範囲と解像度を持つことになります。現在のGRASS領域は $LOCATION/$MAPSET/WIND ファイルに格納されており、それが北、南、東と西の境界、列と行の数、水平方向と垂直方向の空間解像度を定義します。

使用してスイッチを入れるとQGISキャンバスでGRASS領域の可視化をオフすることが可能です ^{現在のGRASSのregion} ボタンを表示します。

領域は「GRASSツール」ドックウィジェットの「領域」タブで変更できます。新しい領域の境界と解像度を入力し、適用をクリックします。キャンバス上のドラッグで領域を指定するをクリックすると、QGISキャンバス上でマウスを使って矩形をドラッグし、対話的に新しい領域を選択することができます。

GRASSモジュール g.region はラスタ分析のための適切な領域の広がりと解像度を定義するための多くのパラメータを提供します。 GRASSツールボックスセクションで説明した、GRASSツールボックスでこれらのパラメータを使用できます。

22.14. GRASSツールボックス 

^{GRASSツールを開く} ボックスは、選択されたGRASS LOCATION と MAPSET 内部のデータを操作するためのGRASSモジュールの機能を提供します。GRASSツールボックスを使用するには、書き込み権限を持っている LOCATION と MAPSET を開く必要があります（ MAPSET を作成した場合、通常は権限が付与されています）。分析中に作成された新しいラスターまたはベクターレイヤーは、現在選択されている LOCATION と MAPSET に書き込む必要があるため、これは必要です。

../../../_images/grass_toolbox_moduletree.png — 図 22.3 GRASSツールボックスおよびモジュールツリー

22.14.1. GRASSモジュールの操作 

GRASSツールボックスのGRASSシェルではほとんどすべての(300以上)GRASSモジュールをコマンドラインインターフェイスで提供しています。もっとユーザーフレンドリな作業環境としてグラフィカルダイアログでは約200のGRASSモジュールと機能が提供されています。

QGISバージョン 3.40 のグラフィカルツールボックスで利用可能なGRASSモジュールの完全な一覧は、GRASS wikiのhttps://grasswiki.osgeo.org/wiki/GRASS-QGIS_relevant_module_list で利用できます。

GRASSツールボックスの内容をカスタマイズすることもできます. この手順はセクション GRASSツールボックスのカスタマイズに記述されています.

図 22.3 に示すように、テーマごとにグループ化されたモジュールツリーや検索可能なモジュールリストタブを使用して、適切なGRASSモジュールを探すことができます。

グラフィカルモジュールアイコンをクリックすると、新しいタブがツールボックス]ダイアログボックスに追加され、オプション、出力とマニュアルの3つの新しいサブタブを提供します。

オプション

オプションタブでは、通常、QGISキャンバスに可視化されたラスターまたはベクターレイヤーを選択し、モジュールを実行するために、さらにモジュール固有のパラメータを入力でき、簡略化モジュールダイアログを提供します。

../../../_images/grass_module_option.png — 図 22.4 GRASSツールボックスモジュールのオプション

提供されているモジュールパラメータは、ダイアログをシンプルに保つために完全でないことがよくあります。さらにモジュールパラメータとフラグを使いたい場合は、GRASSシェルを起動し、コマンドラインでモジュールを実行する必要があります。

QGIS 1.8以降の新機能は、オプションタブの単純化されたモジュール]ダイアログの下にある詳細オプションを表示しますボタンです。現時点では、使用例としてモジュール v.in.ascii にのみ追加されていますが、QGISの将来のバージョンではそれはおそらくより多くの、またはすべての、GRASSツールボックス内のモジュールの一部になります。これはGRASSシェルに切り替える必要なく、完全なGRASSモジュールオプションを使用できます。

出力

../../../_images/grass_module_output.png — 図 22.5 GRASSツールボックスモジュール出力

出力タブは、モジュールの出力状況に関する情報を提供します。実行ボタンをクリックすると、モジュールは出力タブに切り替わり、解析プロセスに関する情報が表示されます。すべてがうまくいくと、最後に 成功しました というメッセージが表示されます。

マニュアル

../../../_images/grass_module_manual.png — 図 22.6 GRASSツールボックスモジュールマニュアル

マニュアルタブはGRASSモジュールのHTMLヘルプページを示しています。それを使って他のモジュールパラメータとフラグを調べたり、モジュールの目的についてのより深い知識を得ることができます。各モジュールのマニュアルページの最後には、メインヘルプ索引 、 テーマ別索引 と 全索引 へのさらなるリンクがあります。これらのリンクは、モジュール g.manual と同じ情報を提供します。

Tip

結果をすぐ表示

計算結果をすぐにマップキャンバスに表示したい場合、モジュールタブの一番下にある '出力を見る' ボタンを利用できます.

22.14.2. GRASSモジュールの例 

以下の例はいくつかの GRASS モジュールの力を示すものです。

22.14.2.1. 等高線を作成する 

最初の例では、標高ラスタ（DEM）からベクタ等高線地図を作成します。セクション GRASS LOCATIONへデータをインポートするで説明したように、ここでアラスカ LOCATION を設定していることを想定しています。

まず、 ^{Mapsetを開く} ボタンをクリックし、アラスカロケーションを選択してロケーションを開きます。
ここで ^{GRASSツールを開く} ボタンでツールボックスを開きます。
ツールカテゴリのリストで、 Raster（ラスタ） ► Surface management（表面管理） ► Generate vector contour lines（ベクタ等高線を生成）をダブルクリックします。
ツール r.contour をシングルクリックすると、上述したように、ツールダイアログが開きます（ GRASSモジュールの操作を参照）。
入力ラスタマップ名に gtopo30 を入力します。
等高線水準間の増加分に値100を入力します（これで等高線が100メートルの間隔で作成されます。）
出力するベクタマップ名に ctour_100 と入力します。
実行`をクリックして処理を開始します。出力ウィンドウに ``成功しました` というメッセージが表示されるまでしばらく待ちます。そして出力を見ると閉じるをクリックします。

これは大きな領域ですので、表示されるまでに数分かかるでしょう。レンダリングを完了したら、レイヤのプロパティウィンドウを開き、ベクタプロパティダイアログのように、等高線が標高ラスタの上に明確に表示されるように線の色を変更できます。

次に、アラスカ州の中心部にある小さな、山岳地帯にズームインします。近くにズームインすると等高線が鋭い角を持っていることがわかります。GRASSではそれらの全体の形状を維持しながら、少しベクタマップを変更する v.generalize ツールを提供しています。このツールは、異なる目的を持ついくつかの異なるアルゴリズムを使用しています。アルゴリズムの一部（すなわち、ダグラス・ペウカーおよび頂点削減）は、頂点の一部を除去することにより、ラインを簡素化します。得られたベクタは、より速く読み込まれます。非常に詳細なベクタを持っているが、非常に小縮尺の地図を作成しているので詳細が不要なとき、このプロセスが便利です。

Tip

単純化ツール

QGISにはちょうどGRASS v.generalize Douglas-Peukerアルゴリズムのように動作するベクタ ► ジオメトリツール ► ジオメトリを簡素化ツールがあることに注意してください。

しかし、この例の目的は異なっています。 r.contour によって作成された等高線には鋭い角があって滑らかにする必要があります。 v.generalize アルゴリズムの中には、ちょうどそれをするChaikenのもの（エルミートスプラインも）があります。これらのアルゴリズムはベクターに追加の頂点を追加し、読み込みがさらに遅くなることに注意してください。

GRASSツールボックスを開き、カテゴリ Vector（ベクタ） ► Develop map（地図を展開） ► Generalization（簡素化）をダブルクリックし、それから v.generalize モジュールをクリックしてそのオプションウィンドウを開いてください。
入力ベクタの名前として「ctour_100」ベクタが表示されていることを確認してください。
アルゴリズムのリストからChaiken'sを選択します。他のオプションはデフォルトのままにしておき、最後の行までスクロールして出力ベクタマップ名フィールドに 'ctour_100_smooth' と入力し、実行をクリックします。
処理にはしばらくかかります。出力ウィンドウに 成功しました が表示されたら、出力を見るをクリックし、それから閉じるをクリックします。
ラスタの背景に明確に表示するとともに元の等高線と対比するため、ベクタの色を変更することができます。新たな等高線は、元の全体的な形状に忠実でありながら、元々のものより角が滑らかになっていることが分かります。

../../../_images/grass_toolbox_vgeneralize.png — 図 22.7 ベクタマップを平滑化するGRASSモジュールv.generalize

Tip

その他にr.contourも使えます

上記の手順は、他の同等の状況で使用できます。降水量データのラスタマップを持っている場合は、例えば、同じ方法が等雨量（一定の降雨量）線のベクタマップを作成するために使用されます。

22.14.2.2. 陰影3D効果の作成 

いくつかの方法は、標高レイヤを表示し、地図に3-D効果を与えるために使用されます。上記のような等高線の使用は、多くの場合、地形図を生成するために選択される1つの一般的な方法です。3-D効果を表示するための別の方法は、陰影起伏によるものです。陰影起伏効果は最初、空での太陽の位置をシミュレートし、各セルに反射率値を与え、その後、各セルの傾きと傾斜方向を計算することにより、DEM（標高）ラスタから作成されます。したがって、太陽に面した斜面が点灯取得します。（影で）太陽から離れて面する斜面を暗くしています。

gtopo30 標高ラスタを読み込むことによって、この例を始めます。GRASSツールボックスを起動し、ラスタカテゴリの下で、ダブルクリックして Spatial analysis（空間分析） ► Terrain analysis（地形分析）を開きます。
そして r.shaded.relief をクリックしてモジュールを開きます。
方位角を 270から315に変更します。
新しい陰影ラスタのために ``gtopo30_shade``と入力し、実行をクリックします。
プロセスが完了すると陰影ラスタが地図に追加されます。これはグレースケールで表示されます。
gtopo30 の陰影起伏と色の両方を一緒に表示するには、目次中の gtopo30 地図の下に陰影起伏図を移動し、次に gtopo30 のプロパティウィンドウを開き、透明度タブに切り替えて約25％の透明度レベルを設定`します。

これでカラーマップと透明度を設定した gtopo30` 標高がグレースケール陰影起伏図の上に表示されるはずです。陰影の視覚効果を見るには、 gtopo30_shade マップをオフにして、それから再びオンにしてください。

GRASSシェルを使用する

QGISのGRASSプラグインは、GRASSを初めて使い、すべてのモジュールやオプションに精通していないユーザ向けに設計されています。そのため、ツールボックスに表示されるモジュールの中には、オプションのすべてが表示されていないものや、まったく表示されていないものがあります。GRASSシェル（またはコンソール）は、ツールボックスのツリーに表示されない追加のGRASSモジュールへのアクセスをユーザーに提供し、また、ツールボックスにあるモジュールで最も簡単なデフォルトパラメーターを持つものに追加のオプションを提供します。この例では、上に出てきた r.shaded.relief モジュールの追加オプションの使い方を示します。

../../../_images/grass_toolbox_shell.png — 図 22.8 GRASSシェル、r.shaded.reliefモジュール

モジュール r.shaded.relief は、陰影起伏効果がより顕著になるように、X-Y座標単位に対して標高値に乗算するパラメータ、zmult を取ることができます。

上述したように gtopo30 標高ラスタを読み込んでからGRASSツールボックスを起動し、GRASSシェルをクリックします。シェルウィンドウでコマンド r.shaded.relief map=gtopo30 shade=gtopo30_shade2 azimuth=315 zmult=3 を入力し、 Enter を押します。
プロセスが終了した後、ブラウズタブにシフトし、新しい gtopo30_shade2 ラスタをダブルクリックしてQGISでそれを表示します。
上で説明したように、目次で陰影起伏ラスタを gtopo30 ラスタの下に移動し、それから、着色された gtopo30 レイヤの透明性を確認します。3-D効果がより強く最初の陰影起伏図と比べて際立っていることがわかります。

../../../_images/grass_toolbox_shadedrelief.png — 図 22.9 GRASSモジュールr.shaded.reliefで作成した陰影を表示します

22.14.2.3. ベクタマップによるラスタ統計 

次の例はGRASS モジュールがラスタデータを集計してベクタマップのそれぞれのポリゴンの列に統計値を追加するものです。

再びアラスカのデータを使って、 GRASS LOCATIONへデータをインポートするを参照して、 shapefiles/trees.shp ファイルをGRASSにインポートします。
ここで、中間のステップが必要です：重心をインポートされた樹木図に追加して（境界と重心の両方を含む）完全なGRASSエリアベクタを作る必要があります。
ツールボックスで Vector（ベクタ） ► Manage features（地物を管理）を選択して v.centroids モジュールを開きます。
出力するベクトルマップ名として 'forest_areas' と入力してモジュールを実行して下さい.
次に、 forest_areas ベクタを読み込んで、森林の種類 - 落葉、常緑、混合 - を異なる色で表示します：レイヤのプロパティウィンドウのシンボロジタブの中で、凡例タイプから「ユニークな値」選び、分類フィールドに「VEGDESC」を設定します。（ベクタセクションのシンボロジプロパティ中のシンボロジタブの説明を参照してください。）
次に GRASS ツールボックスを再び開いてベクター ► ベクター更新を他の地図で開いて下さい.
v.rast.stats モジュールをクリックして gtopo30 と forest_areas と入力して下さい.
追加で必要なパラメータは1つだけです： column prefix elev と入力し、実行をクリックします。これは計算量が多く、長時間（おそらく2時間程度）実行される。
最後に、 forest_areas 属性テーブルを開き、各森林のポリゴンのために、いくつかの新しい列が、 elev_min 、 elev_max 、 elev_mean などを含め、追加されていることを確認してください。

22.14.3. GRASSツールボックスのカスタマイズ 

ほぼすべてのGRASSモジュールは、GRASSツールボックスに追加できます。XMLインターフェイスは、ツールボックス内のモジュールの外観とパラメーターを設定する非常に単純なXMLファイルを解析するために提供されます。

モジュール v.buffer （v.buffer.qgm）を生成するためのサンプルXMLファイルは次のようになります。

<?xml version="1.0" encoding="UTF-8"?>
<!DOCTYPE qgisgrassmodule SYSTEM "http://mrcc.com/qgisgrassmodule.dtd">

<qgisgrassmodule label="Vector buffer" module="v.buffer">
        <option key="input" typeoption="type" layeroption="layer" />
        <option key="buffer"/>
        <option key="output" />
</qgisgrassmodule>

The parser reads this definition and creates a new tab inside the Toolbox when you select the module.