9.4. Lesson: 補足の実習

このレッスンでは、QGISでの完全なGIS解析を通して案内します。

注釈

LinfinitiとS Motala（ケープ半島工科大学）が開発したレッスン

9.4.1. 問題の状態

あなたはケープ半島とその周辺で、珍しいフィンボスの植物種に適した生息地である領域を見つける任務を課されています。ケープ半島内での調査地域の範囲は次のとおりです。Melkbosstrandの南、ストランドの西。問題の植物種は以下の嗜好を示すことを植物学者が示してくれています：

• 東向き斜面に生息します。

• 15％と60％の間の勾配の斜面に生息します。

• 年間総降水量1200ミリメートル超の領域に生息します。

• どの集落からも少なくとも250メートル離れたところで発見されます。

• それが発生した植生の面積は、面積が少なくとも6000平方メートルである必要があります。

ケープ自然保護ボランティアとして、あなたは自宅に最も近い適当な土地の一画でその植物を探すことに合意しました。GISのスキルを使用してどこを見に行くべきか決定しましょう。

9.4.2. ソリューションの概要

The data for this exercise can be found in the exercise_data/more_analysis folder. You are going to find the candidate area that is closest to a house of your choice. The solution will involve:

1. DEMを分析し、東向き斜面および該当する勾配を有する斜面を見つける。

2. 降雨ラスターを分析し、該当する降雨量の領域を見つける。

3. ゾーニングベクターレイヤーを分析し、集落から離れていて該当するサイズである領域を見つける。

9.4.3. 地図を設定する

1. 画面の一番右下にある「CRSステータス」ボタンをクリックします。表示される画面の CRS タブの下に ワールド座標参照系 ボックスが表示されるでしょう。

2. このボックスで、 投影座標系 ‣ ユニバーサル横メルカトル(UTM) に移動します。

3. Select the entry WGS 84 / UTM zone 33S (with the EPSG code 32733).

4. Click OK. The map is now in the UTM33S coordinate reference system.

5. プロジェクトに名前を付けて保存 ツールバーボタンをクリックして地図を保存するか、 ファイル ‣ プロジェクトに名前を付けて保存... メニュー項目を使用します。

6. Save the map in a directory called Rasterprac that you should create somewhere on your computer. You will save whatever layers you create in this directory as well.

9.4.4. 地図へデータを読み込む

データを処理するために、地図キャンバスに必要なレイヤー（通りの名前、ゾーン、降雨、DEM）をロードする必要があります。

9.4.4.1. For vectors...

1. Click on the ^{Open Data Source Manager} button and enable the Vector tab in the dialog that appears, or use the Layer ‣ Add Layer ‣ Add Vector Layer... menu item.

2. Ensure that the File radio button is selected.

3. Click on the ... button to browse for Vector dataset(s).

4. In the dialog that appears, open the exercise_data/more_analysis/streets directory.

5. Select the file Street_Names_UTM33S.shp.

6. 開く をクリックします。

   The dialog closes and shows the original dialog, with the file path specified in the text field next to the Vector dataset(s) button. This allows you to ensure that the correct file is selected. It is also possible to enter the file path in this field manually, should you wish to do so.

7. Click Add. The vector layer will load in your map. Its color is automatically assigned. It will be changed later.

8. Rename the layer to Streets.

   1. Right-click on it in the Layers panel (by default, the pane along the left-hand side of the screen).

   2. Click Rename in the dialog that appears and rename it, pressing the Enter key when done.

9. Repeat the vector adding process, but this time select the Generalised_Zoning_Dissolve_UTM33S.shp file in the Zoning directory.

10. Rename it to Zoning.

9.4.4.2. For rasters...

1. Click on the ^{Open Data Source Manager} button and enable the Raster tab in the dialog that appears, or use the Layer ‣ Add Layer ‣ Add Raster Layer... menu item.

2. Ensure that the File radio button is selected.

3. 該当するファイルに移動し、それを選択し、 開く をクリックしてください。

4. Do this for each of the two raster files. The files you want are DEM/reproject/DEM.tif and Rainfall/reprojected/rainfall.tif.

5. Rename the rainfall raster to Rainfall (with an initial capital). Initially when you load them, the images will be gray rectangles. Don't worry, this will be changed later.

6. 地図を保存します。

何が起こっているのかを適切に確認するために、レイヤーのためのシンボルを変更する必要があります。

9.4.5. ベクターレイヤーのシンボル体系を変更する

1. In the Layers panel, right-click on the Streets layer.

2. 表示されているメニューで プロパティ を選択します。

3. Switch to the Symbology tab in the dialog that appears.

4. 街路 レイヤーの現在の色を示す四角が付いた 変更 とラベルのついたボタンをクリックしてください。

5. 表示されるダイアログで新しい色を選択します。

6. OK をクリックします。

7. レイヤープロパティ ダイアログの OK をもう一度クリックします。これは 街路 レイヤーの色を変更します。

8. Zoning レイヤーについて同様のプロセスに従って、そのための適切な色を選択してください。

9.4.6. ラスターレイヤーのシンボル体系を変更する

ラスターレイヤーシンボル体系はいくらか違います。

1. Rainfall ラスターの プロパティ ダイアログを開きます。

2. Switch to the Symbology tab. You'll notice that this dialog is very different from the version used for vector layers.

3. 標準偏差を使用する が選択されていることを確認します。

4. Change the value in the associated box to 2.00 (it should be set to 0.00 by default).

5. コントラスト強調 見出しの下に、 Current ドロップダウンリストの値を 最大最小に引き延ばす に変更します。

6. OK をクリックします。「降水量」ラスターは、見える場合は、画素毎に異なる輝度値が見えるように、色を変更しなければなりません。

7. Repeat this process for the DEM, but set the standard deviations used for stretching to 4.00.

9.4.7. レイヤー順序を変更する

1. In the Layers panel, click and drag layers up and down to change the order they appear in on the map.

2. Newer versions of QGIS may have a Control rendering order checkbox beneath the Layers panel. Ensure that it is checked.

これですべてのデータがロードされ正常に表示されますので、分析を開始できます。クリップ操作が最初に行われるならば、最高です。使用するつもりのない領域の値を計算することに処理能力を無駄使いしないようにするためです。

9.4.8. 正しい地区の検索

1. Load the vector layer admin_boundaries/Western_Cape_UTM33S.shp into your map.

2. Rename it to Districts.

3. Right-click on the Districts layer in the Layers panel.

4. 表示されるメニューで、 クエリ... メニュー項目を選択します。 クエリビルダー ダイアログが表示されます。

今、地区の次のリストだけを選択するためのクエリを作成します：

• Bellville

• Cape

• Goodwood

• Kuils River

• Mitchells Plain

• Simons Town and

• Wynberg.

1. フィールド リスト中で NAME_2 フィールドをダブルクリックします。以下の SQLのwhere句 のテキストフィールドに表示されます。

2. Click the = button; an = sign is added to the SQL query.

3. （現在空） 値 リストの下の すべて ボタンをクリックします。短い遅延の後、これは選択したフィールドの値（NAME_2 ）を持つ 値 リストを投入します。

4. 値 リストの中で Bellville をダブルクリックします。前述のように、これはSQLクエリに追加されます。

In order to select more than one district, you'll need to use the OR boolean operator.

1. OR ボタンをクリックすると、それがSQLクエリに追加されます。

2. 上記と同様のプロセスを使用して、既存のSQLクエリに次の行を追加します。

   "NAME_2" = 'Cape'
   

3. 別の OR 演算子を追加し、その後、同様の方法で、上記の地区のリストを通して作業を続けます。

   最終的なクエリはこうであるべきです

   "NAME_2" = 'Bellville' OR "NAME_2" = 'Cape' OR "NAME_2" = 'Goodwood' OR
   "NAME_2" = 'Kuils River' OR "NAME_2" = 'Mitchells Plain' OR "NAME_2" =
   'Simons Town' OR "NAME_2" = 'Wynberg'
   

4. OK をクリックします。地図に表示される地区は現在、上記のリストのものに限定されています。

9.4.9. ラスターのクリップ

今、興味の領域を持っている、このエリアへのラスターをクリップできます。

1. 見えるレイヤーは DEM 、Rainfall 、Districts レイヤーだけであることを確認してください。

2. Districts は見えるように最上になければなりません。

3. ラスター ‣ 抽出 ‣ Clipper メニュー項目を選択することにより、クリッピングダイアログを開きます。

4.  入力ファイル（ラスター）レイヤー ドロップダウンリストで DEM レイヤーを選択します。

5. 選択... ボタンをクリックして、 出力ファイル テキストフィールドに出力先を指定します。

6. Rasterprac ディレクトリに移動します。

7. ファイル名を入力します。

8. ファイルを保存します。データ値なし チェックボックスはチェックなしのままにしておきます。

9. 正しいラジオボタンが選択されていることを確保することによって 範囲 クリップモードを使用してください。

10. 地区を含む領域が選択されるように、キャンバス内のエリアをクリックしてドラッグ。

11. Check the Open output file after running algorithm box.

12. Click Run.

13. クリップ操作が完了した後、クリップ ダイアログを閉じないでください。（そうすると、すでに定義されているクリップ領域を失う原因となります。）

14. 入力ファイル（ラスター） ドロップダウンリストで Rainfall ラスターを選択し、別の出力ファイル名を選択します。

15. 他のオプションを変更しないでください。以前に描いた既存のクリップ領域を変更しないでください。すべてのものを同じままに残し、実行 をクリックしてください。

16. 第二のクリップ操作が完了したら、クリップ ダイアログを閉じることができます。

17. 地図を保存します。

9.4.10. 地図をクリーンアップします。

1. レイヤー パネルから元の 降雨 と DEM レイヤーを削除します：

2. これらのレイヤー上で右クリックし、 削除 を選択します。

   注釈

   これは、記憶装置からデータを削除しません。地図から外すだけです。

3. 街路 レイヤーのラベルを無効にします。

   1. ラベル付け ボタンをクリックします。

   2. このレイヤーにラベル付け ボックスのチェックをオフにします。

   3. OK をクリックします。

4. 再度 街路 をすべて表示します。

   1. レイヤー パネルでレイヤーを右クリックします。

   2. クエリ を選択します。

   3. 表示された クエリ ダイアログで、 クリア ボタンをクリックし、 OK をクリックします。

   4. データがロードされるまで待ちます。すべての街は今、表示されます。

5. 以前のようにラスターシンボルを変更します（ *ラスターレイヤーのシンボル体系を変更する 参照）。

6. 地図を保存します。

7. 今 レイヤー パネル中のベクターレイヤーを、横にあるボックスをオフにして非表示できます 。これにより、地図のレンダリングが速くなりますし、時間の節約になります。

陰影起伏を作成するには、この目的のために書かれたアルゴリズムを使用する必要があります。

9.4.11. 陰影図の作成

1. In the Layers panel, ensure that the DEM is the active layer (i.e., it is highlighted by having been clicked on).

2. Click on the Raster ‣ Analysis ‣ Hillshade menu item to open the Hillshade dialog.

3. 出力レイヤーのための適切な場所を指定し、それを hillshade と呼びます。

4. Check the Open output file after running algorithm box.

5. Click Run.

6. 処理が完了するのを待ちます。

新しい hillshade レイヤーが レイヤーリスト に現れています。

1. レイヤーリスト 中で hillshade レイヤーを右クリックし、プロパティ ダイアログを呼び出します。

2. 透明度 タブをクリックして透明度スライダを 80% に設定します。

3. ダイアログで 実行 をクリックします。

4. 透明な陰影起伏がクリップされたDEM上に重ね合わされたときに効果に注意されたいです。

9.4.12. 傾斜

1. ラスター ‣ 地形分析 メニュー項目をクリックします。

2. クリップされたDEMを入力レイヤーとして 傾斜 分析タイプを選択します。

3. 出力の目的のために適切なファイル名と場所を指定します。

4. Check the Open output file after running algorithm box.

5. Click Run.

勾配画像が計算され、地図に追加されました。しかし、いつものように単なる灰色の四角形です。何が起こっているかを適切に確認するには、次のようにシンボルを変更します。

1. （いつものようにレイヤーの右クリックメニューを経由して） レイヤー Properties ダイアログを開きます。

2. Click on the Symbology tab.

3. ( カラーマップ ドロップダウンメニューで) グレースケール となっているところを 疑似カラー に変更します。

4. 標準偏差を使用 ラジオボタンが選択されていることを確認します。

9.4.13. 傾斜方位

傾斜を計算するのと同じアプローチを使用しますが、最初のダイアログボックスで 傾斜方向 を選択します。

定期的に地図を保存することを忘れないで下さい。

9.4.14. ラスターを再分類する

1. メニュー項目 ラスター ‣ ラスター計算機 をクリックしてください。

2. 出力レイヤーのための場所として Rasterprac ディレクトリを指定します。

3. Ensure that the Open output file after running algorithm box is selected.

左の ラスタバンド リストに、レイヤ パネル中のすべてのラスタレイヤが表示されます。傾斜レイヤが 傾斜 という名であれば、それは 傾斜@1 、と表示されます。

傾きが 15 と 60 度の間にあることが必要です。したがって 15 未満 60 超はすべて排除しなければなりません。

1. インターフェイスのリスト項目とボタンを使用して以下の式を作成します：

   ((slope@1 < 15) OR (slope@1 > 60)) = 0
   

2. 適切な場所とファイル名で 出力レイヤー フィールドを設定してください。

3. Click Run.

Now find the correct aspect (east-facing: between 45 and 135 degrees) using the same approach.

1. 次の式を組み立てます:

   ((aspect@1 < 45) OR (aspect@1 > 135)) = 0
   

2. Find the correct rainfall (greater than 1200mm) the same way. Build the following expression:

   (rainfall@1 < 1200) = 0
   

すべてのラスターを再分類した今、それらは（それらは正しく地図に追加されていると仮定して）自分の地図に灰色の長方形として表示が表示されます。（ 1 と 0 は真または偽を意味する）2つのクラスだけでラスターデータを正しく表示するには、それらのシンボルを変更する必要があります。

9.4.15. 再分類されたレイヤーにスタイルを設定する

1. Open the Symbology tab in the layer's Properties dialog as usual.

2. バンドから最小/最大値ロード 見出しの下の 実際（遅い） ラジオボタンを選択します。

3. 読み込み ボタンをクリックします:

今 カスタム最小/最大値 フィールドにはれぞれ 0 と 1 が取り込まれているはずです。（そうでない場合、データの再分類で間違いがあったので、その部分を再びやり直す必要があります。）

1. コントラスト強調 見出しの下で Current ドロップダウンリストを MinMaxへ引き延ばしする に設定します。

2. OK をクリックします。

3. 3つのすべてのためにこれを行うには、ラスターを再分類し、作業を保存することを忘れないでください！

残る基準は、その領域が都市地域から 250m 離れてなければならないことだけです。この要件をは、計算する領域が農村地域の端から 250m 以上であることを確実にすることで満たします。したがってまず、すべての農村地域を見つける必要があります。

9.4.16. 農村地域を検索する

1. レイヤー パネルですべてのレイヤーを非表示にします。

2. Zoning を再表示します。

3. その上で右クリックして クエリ ダイアログを起動します。

4. 次のクエリをビルドします:

   "Gen_Zoning" = 'Rural'
   

   行き詰まったら Streets クエリを構築するための以前の説明書を参照してください。

5. それを実行したら クエリ ダイアログを閉じます。

Zoning レイヤーからポリゴンのコレクションが表示されるはずです。新しいレイヤーファイルにこれらを保存する必要があります。

1. Zoning の右クリックメニューで 名前を付けて保存... を選択します。

2. Zoning ディレクトリの下にレイヤーを保存します。

3. 出力ファイルの名前を rural.shp とします。

4. OK をクリックします。

5. レイヤーを地図に追加します。

6. ベクター ‣ ジオプロセッシングツール ‣ ディゾルブ をクリックします。

7. 選択地物のみ使用 ボックスのチェックはないままで、入力ベクターレイヤーとして rural レイヤーを選択します。

8. Leave empty the Dissolve field(s) option to combine all selected features in a single one.

9. Zoning ディレクトリの下にレイヤーを保存します。

10. Check the Open output file after running algorithm box.

11. Click Run.

12. ディゾルブ ダイアログを閉じます。

13. rural と Zoning レイヤーを外します。

14. 地図を保存します。

ここで農村地域の端から 250m 内にある領域を除外する必要があります。これは以下に説明するように、負のバッファを作成することによって行います。

9.4.17. 負のバッファを作成する

1. メニュー項目 ベクター ‣ バッファ ‣ ジオプロセッシングツール をクリックしてください。

2. 表示されるダイアログで、入力ベクターレイヤーとして rural_dissolve レイヤーを選択します（ 選択された地物のみ使用 はチェックすべきではありません）。

3. バッファ距離 ボタンを選択し、関連するフィールドに値 -250 を入力してください。負の値は、バッファが内部バッファでなければならないことを意味しています。

4. バッファ結果をディゾルブする ボックスにチェックします。

5. 他の農村のベクターファイルと同じディレクトリに出力ファイルを設定します。

6. Name the output file rural_buffer.shp.

7. 保存 をクリックします。

8. OK をクリックして処理が完了するのを待ちます。

9. 現れたダイアログで はい を選択します。

10. バッファ ダイアログを閉じます。

11. rural_dissolve レイヤーを削除します。

12. 地図を保存します。

3つの既存のラスターと同じ分析に農村地帯を組み込むためには、それは同様にラスター化する必要があります。しかし、ラスターは、分析のために互換性があるようにするために、同じサイズにする必要があります。そのため、ラスター化する前に次の3つのラスターと同じエリアにベクターをクリップする必要があります。ベクターは別のベクターによってクリップできるのみなので、最初にラスターと同じサイズのバウンディングボックスポリゴンを作成する必要があります。

9.4.18. バウンディングボックスベクターを作成する

1. メニュー項目 レイヤー ‣ 新規 ‣ 新規シェープファイルレイヤー... をクリックします。

2. タイプ 見出しの下で ポリゴン ボタンを選択します。

3. Click Specify CRS and set the coordinate reference system WGS 84 / UTM zone 33S : EPSG:32733.

4. OKをクリックします。

5. 新しいベクターレイヤー ダイアログでも OK をクリックします。

6. Zoning ディレクトリにベクターを保存します。

7. Name the output file bbox.shp.

8. 新しい bbox レイヤーと再分類ラスターのレイヤーを除くすべてのレイヤーを非表示にします。

9. Ensure that the bbox layer is highlighted in the Layers panel.

10. 表示>ツールバー メニュー項目に移動し、 デジタル化 が選択されていることを確認します。その後、その上に鉛筆やコーキーと、ツールバーのアイコンが表示されるはずです。 編集を切り替え ボタンです。

11. 編集を切り替え ボタンをクリックして 編集モード に入ってください。これにより、ベクターレイヤーが編集できるようになります。

12. 地物を追加 ボタンをクリックします、編集を切り替え ボタンの近くにあるはずです。これは、二重の矢印ボタンの後ろに隠れてすることができます。 デジタイズ ツールバーの非表示のボタンもしそうなら、表示するために二重矢印をクリックしてください。

13. 地物を追加 ツールをアクティブにして、ラスターの角を左クリックしてください。それが正確であることを保証するために、マウスのホイールでズームインする必要があるかもしれません。このモードでは、地図をパンするには、マウスの中央ボタンかマウスホイールをクリックしたまま地図内でドラッグします。

14. 第四及び最終点に関しては、形状を確定するために右クリックします。

15. シェープIDの任意の番号を入力します。

16. OK をクリックします。

17. 編集を保存 ボタンをクリックします。

18. 編集を切り替え ボタンをクリックして、編集セッションを停止してください。

19. 地図を保存します。

これでバウンディングボックスができましたので、農村地域のバッファレイヤーをクリップするためにそれを使用できます。

9.4.19. ベクターレイヤーをクリップする

1. bbox と rural_buffer レイヤーだけが、後者を上にして、表示されていることを確認してください。

2. ベクター>ジオプロセッシングツール>クリップ メニュー項目をクリックしてください。

3. 表示されたダイアログで、入力ベクターレイヤーを rural_buffer に、クリップレイヤーを bbox に、両方とも 選択地物のみ使用 ボックスをチェックしないで、設定します。

4. Zoning ディレクトリの下に、出力ファイルを置きます。

5. Name the output file rural_clipped.

6. OK をクリックします。

7. TOCにレイヤーを追加するように求められたら、 はい をクリックしてください。

8. ダイアログを閉じます。

9. 3つのベクターを比較し、自分のために結果を参照してください。

10. bbox と rural_buffer レイヤーを削除し、地図を保存します。

これでラスター化する準備ができました。

9.4.20. ベクターレイヤーをラスター化する

作成した新しいラスターのピクセルサイズを指定する必要がありますので、まず既存のラスターの1つの大きさを知っておく必要があります。

1. 3つの既存のラスターのいずれかの プロパティ ダイアログを開きます。

2. メタデータ タブに切り替えます。

3. Make a note of the X and Y values under the heading Dimensions in the Metadata table.

4. プロパティ ダイアログを閉じます。

5. ラスター ‣ 変換 ‣ ラスター化 メニュー項目 をクリックしてください。データセットがサポートされていないことについて警告を受けることがあります。それをクリックして消し、無視します。

6. 入力レイヤーとして rural_clipped を選択します。

7. Zoning ディレクトリ内の出力ファイルの場所を設定します。

8. Name the output file rural_raster.tif.

9. Check the New size box and enter the X and Y values you made a note of earlier.

10. キャンバスにロードします ボタンをチェックします。

11. Click the pencil icon next to the text field which shows the command that will be run. At the end of the existing text, add a space and then the text -burn 1. This tells the Rasterize function to "burn" the existing vector into the new raster and give the areas covered by the vector the new value of 1 (as opposed to the rest of the image, which will automatically be 0).

12. OK をクリックします。

13. それが計算された後、新たなラスターが地図に表示されるはずです。

14. 新しいラスターはグレーの長方形のようになります - 再分類ラスターのために行ったように、表示スタイルを変更できます。

15. 地図を保存します。

Now that you have all four criteria each in a separate raster, you need to combine them to see which areas satisfy all the criteria. To do so, the rasters will be multiplied with each other. When this happens, all overlapping pixels with a value of 1 will retain the value of 1, but if a pixel has the value of 0 in any of the four rasters, then it will be 0 in the result. In this way, the result will contain only the overlapping areas.

9.4.21. ラスターを組み合わせる

1. ラスター ‣ ラスター計算機... メニュー項目をクリックします。

2. （それらを呼ばれるものに応じて、そのレイヤーのための適切な名前で）次の式を構築します。

   [Rural raster] * [Reclassified aspect] * [Reclassified slope] *
   [Reclassified rainfall]
   

3. 出力場所を Rasterprac に設定します。

4. Name the output raster cross_product.tif.

5. Ensure that the Open output file after running algorithm box is checked.

6. Click Run.

7. 他の再分類ラスターのスタイルを設定するのと同じ方法で、新しいラスターのシンボルを変更します。新しいラスターが正しく、すべての基準が満たされている領域を表示します。

To get the final result, you need to select the areas that are greater than 6000m^2. However, computing these areas accurately is only possible for a vector layer, so you will need to vectorize the raster.

9.4.22. ラスターをベクター化する

1. メニュー項目 ラスター ‣ 変換 ‣ ポリゴン化 をクリックしてください。

2. Select the cross_product.tif raster.

3. 出力場所を Rasterprac に設定します。

4. Name the file candidate_areas.shp.

5. Ensure that Open output file after running algorithm is checked.

6. Click Run.

7. 処理が完了したらダイアログを閉じます。

All areas of the raster have been vectorized, so you need to select only the areas that have a value of 1.

1. 新しいベクターのため クエリ ダイアログを開きます。

2. クエリのビルド:

   "DN" = 1
   

3. OK をクリックします。

4. Create a new vector file from the results by saving the candidate_areas vector after the query is complete (and only the areas with a value of 1 are visible). Use the Save as... function in the layer's right-click menu for this.

5. Save the file in the Rasterprac directory.

6. 候補領域のみ.shp と名前を付けます。

7. 地図を保存します。

9.4.23. 各ポリゴンの面積を計算する

1. 新規ベクターレイヤーの右クリックメニューを表示します。

2. 属性テーブルを開く を選択します。

3. テーブルの下に沿った 編集モード切り替え ボタンをクリックするか、Ctrlキー+E を押します。

4. テーブル底部に沿った フィールド計算機を開く ボタンをクリックするか、Ctrlキー+I を押します。

5. Under the New field heading in the dialog that appears, enter the field name area. The output field type should be an integer, and the field width should be 10.

6. フィールド演算式 に、タイプします:

   $area
   

   計算された値で、このフィールドの計算は、ベクターレイヤーの各ポリゴンの面積を計算し、次いで（ 面積 と呼ばれる）新しい整数列を投入することを意味します。

7. OK をクリックします。

8. id と呼ばれる別の新しいフィールドに同じことを行います。フィールド演算式 中にキーボード入力します：

   $id
   

   これは、各ポリゴンは識別目的のためのユニークなIDを持っていることを保証します。

9. 再び 編集モード切り替え をクリックし、編集内容を保存するようプロンプトが表示されたら、そうします。

9.4.24. 与えられたサイズの面積を選択する

今この面積が知られています：

1. Build a query (as usual) to select only the polygons larger than 6000m^2. The query is:

   "area" > 6000
   

2. Save the selection as a new vector layer called solution.shp.

今はもう解領域が出ているので、家に最も近いものを選択します。

9.4.25. あなたの家のデジタイズ

1. 以前のように新しいベクターレイヤーを作成しますが、今回は type 値が Point となるよう選択します。

2. それが正しいCRSであることを確認します！

3. Name the new layer house.shp.

4. 新規レイヤー作成を終了します。

5. （新規レイヤーが選択された状態で）編集モードにします。

6. 街路をガイドとして使用しながら、家または居住の他の現在の場所がある点をクリックします。家を見つけるのを助けるために他のレイヤーを開く必要があります。近くのどこにも住んでいない場合は、単に家が考えられるかもしれない街の中のどこかをクリックします。

7. シェープIDの任意の番号を入力します。

8. OK をクリックします。

9. 編集内容を保存し、編集モードを終了します。

10. 地図を保存します。

どれが家に最も近いかを決定するため、ソリューションエリアのポリゴンの重心（「質量中心」）を見つける必要があります。

9.4.26. ポリゴンの中心点の算出

1. Click on the Vector ‣ Geometry Tools ‣ Centroids menu item.

2. 入力レイヤーを solution.shp と指定します。

3. 出力場所を Rasterprac と指定します。

4. Call the destination file solution_centroids.shp.

5. Check Open output file after running algorithm to add the result to the TOC (Layers panel).

6. Click Run and close the dialog.

7. 新しいレイヤーを、見えるようレイヤーの順序の一番上にドラッグします。

9.4.27. あなたの家に最も近い重心の計算

1. メニュー項目 ベクター ‣ 分析ツール ‣ 距離行列 をクリックします。

2. The input layer should be your house, and the target layer solution_centroids. Both of these should use the id field as their unique ID field.

3. 出力マトリクス型は linear になるはずです。

4. 適切な出力場所と名前を設定します。

5. OK をクリックします。

6. テキストエディタでファイルを開きます(またはスプレッドシートにインポート)。最短 距離 に関連付けられたターゲットIDに注意してください。同距離に複数あることもあります。

7. Build a query in QGIS to select only the solution areas closest to your house (selecting it using the id field).

これは、研究の質問への最終的な答えです。

提出のために、お好みの魅力的なラスター (例えば DEM とか 傾斜 ラスター、など)の上に半透明の陰影起伏のレイヤーを、また、最も近いソリューション領域(複数可)のポリゴンだけでなく、家を含めます。 出力図を作成する際は地図作成のためのすべてのベストプラクティスに従ってください。