災害・自然破壊による環境劣化のモニタリング・評価システムの開発
環境劣化評価 QGIS プラグイン 利用ガイド

1. 地球地図データの用意(1) ダウンロード
地球地図データの入手先は国土地理院の以下のサイトにて紹介されています。
http://www.gsi.go.jp/kankyochiri/gm_data.html#gm_data
各国の地球地図データが、地球地図国際運営委員会(ISCGM)のサイトにて提供されております。一方、日本国内のデータであれば、国土地理院のサイトからもダウンロードが可能です。
http://www.gsi.go.jp/kankyochiri/gm_jpn.html(2) データフォーマットの変換
地球地図データは、BIL・GML・JPG・Shape・TIFF 等の形式で提供されています。この
うち、QGIS で扱ううえで注意の必要な BIL・GML について説明します。
BIL の場合

QGIS のデータ変換ライブラリに使われている GDAL は、BIL フォーマットをサポートしています。したがって、QGIS で「ラスタレイヤの追加」メニューを実行し、拡張子.bil のファイルを選択すれば、レイヤとして追加することが可能です。
ただし、土地利用図における海域のように特定の値が大半を占めるような領域（日本国内であれば、沖縄や小笠原近辺）を取り扱う際には注意が必要です。QGIS のデフォルトのス
タイル設定では、累積カウント数で 2.0 – 98.0 % の範囲に収まる値しか採用されないため、次のように一面が黒く塗りつぶされた画像が表示されることがあります。

このような場合は、スタイル設定の「最小値/最大値のロード」を「最小/最大」に変更して読み込みボタンをクリックしてください。

これにより、割合の低いセルの値も表示されるようになります。

GML の場合

QGIS の「ベクタレイヤの追加」メニューからそのまま GML を読み込むと、緯度・経度が反転して表示されてしまいます。

OGR の config オプションを利用して、一度正常な緯度・経度の shape を作成してから QGISに読み込ませるようにします。実行する OGR コマンドは次の通りです。
> ogr2ogr -f "ESRI Shapefile" --config GML_INVERT_AXIS_ORDER_IF_LAT_LONG YES--config GML_CONSIDER_EPSG_AS_URN YES 出力ファイル(shp) 入力ファイル(gml)

出力された shp ファイルをベクタレイヤとして追加すると、緯度・経度の反転が解消されていることが分かります。

2. インストール(1) プラグインのコピー
環境劣化評価プラグイン(EnviromentalEvaluation)を QGIS のプラグインフォルダにフォルダごとコピーします。QGIS のプラグインフォルダは以下の通りです。
・Windows
C:¥Users¥ユーザ名¥.qgis2¥python¥plugins※「ユーザ名」はログインしているアカウントのフォルダ名になります・Linux$HOME/.qgis2/python/plugins(2) プラグインの登録
QGIS に環境劣化評価プラグインを登録します。
① QGIS を起動し、メニューから「プラグイン」→「プラグインの管理とインストール」を選択して「プラグインマネージャ」のダイアログを表示させます。

② 「プラグイン＞インストール済」のリストを表示し、リストから「Evaluation ofenvironmental deterioration」プラグインを選択し、チェックボックスにチェックを入れます。

③ 正常にインストールされると、「プラグイン」メニューに「環境劣化評価」および「オーバーレイ＆マトリクス集計」、「環境劣化影響分析」が追加されます。また、プラグインアイコンとして、２つのアイコン（緑地に OMS,EDA の文字）も追加されます。

3. 利用方法
環境劣化評価プラグインは以下の手順で利用します。
(1) 入力データの読み込みと表示
計算に使用するデータを読み込みます。計算に用いるデータは以下の通りです。
① ラスタデータ
土地利用区分等のラスタデータで、単一バンド、整数型のデータが利用可能です。
② ポリゴンデータ
行政区分等のポリゴンデータで、このデータの各ポリゴンに対してラスタデータの値毎にカウントし、集計します。
以下に読み込み例を示します。lu レイヤがラスタデータ、AdmArea レイヤがポリゴンデータになります。

(2) オーバーレイ＆マトリクス集計
各ポリゴンと重なるラスタデータの値をカウントし、集計します。集計結果はシェープファイルとして出力します。
① 「オーバーレイ＆マトリクス集計」プラグインの起動
メニューから「プラグイン」→「環境劣化評価」→「オーバーレイ＆マトリクス集計」を選択し、ダイアログを表示させます。

② 入力データの選択
入力データを選択します。ラスタ、ポリゴンとも地図上に表示されているレイヤから選択することができます。
③ 出力シェープファイルの指定
計算結果を出力するシェープファイルを指定します。「ブラウズ」ボタンを押すとファイル選択のダイアログが表示されます。

既存のファイルを上書きする場合は、確認のダイアログが表示されます。ただし、入力ポリゴンと同じファイルを指定することは出来ません（実行時にエラーとなります）。

④ 計算実行
入出力データの選択、指定が完了したら実行ボタンを押して計算を行います。計算はラスタデータのポリゴンへの変換、データの集計の２ステップで行います。それぞれ進行状況がプログレスバーに表示されます。

計算が完了すると「結果を地図に表示しますか？」というダイアログが表示されます。
「はい」を選択すると地図上に出力シェープファイルで指定したポリゴンが表示されます。

⑤ 計算結果
計算結果は、出力したシェープファイルの属性情報として登録されています。ラスタデータの値毎にフィールドを作成し、そのフィールド毎にカウントしたラスタの数が登録されます。結果を登録するフィールド名は「C + ラスタの値」となります。例えば、ラスタ値が ”10” であった場合、対応するフィールド名は “C10” となります。このフィールドのデータ型は整数型となります。

(3) 環境劣化影響分析「オーバーレイ＆マトリクス集計」の結果を使って、トータルコストを計算します。
① 「オーバーレイ＆マトリクス集計」プラグインの起動
メニューから「プラグイン」→「環境劣化評価」→「環境劣化影響分析」を選択し、ダイアログを表示させます。

② 入力データ（ポリゴン）の選択「オーバーレイ＆マトリクス集計」の計算結果のポリゴンを選択します。地図上に表示されているポリゴンレイヤから選択することができます。
③ コストパラメータプファイルの指定計算を行うためのコストパラメータファイルを指定します。コストパラメータファイル
は、計算を行うフィールド名（文字型）,乗じる係数（コスト単価、実数型）のレコードからなるテキスト（CSV）ファイルです。以下に例を示します。
C10,1.24
C20,2.2
C30,3.3
C40,5.1
C50,8.678
C60,10.45
C70,20.2
C80,30.1
C90,50④ 出力シェープファイルの指定
計算結果を出力するシェープファイルを指定します。「ブラウズ」ボタンを押すとファイル選択のダイアログが表示されます。

⑤ 計算実行
入出力データの選択、指定が完了したら実行ボタンを押して計算を行います。進行状況がプログレスバーに表示されます。
計算が完了すると「結果を地図に表示しますか？」というダイアログが表示されます。
「はい」を選択すると地図上に出力シェープファイルで指定したポリゴンが表示されます。

⑥ 計算結果
計算結果は、出力したシェープファイルの “TOTAL” というフィールドに登録されます。