14.2. Lesson: 地図をジオリファレンスする

林業においてよくある作業は、森林地域についての情報の更新でしょう。ある地域についての前回の情報が何年も前のもので、アナログで(つまり紙で)集められていたり、デジタイズされたが残っているのはその目録データの紙版しかなかったり、といったことがあり得ます。

その情報を、例えば後の目録と比較するために、GISで使用したいと思うことはよくあります。これは、GISソフトウェアを使って手元の情報をデジタイズする必要があることを意味します。しかしデジタイズを始める前に、紙の地図をスキャンしジオリファレンスするという重要な最初のステップがあります。

このレッスンの目標: QGISでジオリファレンスサツールを使用する方法を学ぶ。

14.2.1. basic 地図をスキャンする

最初にしなければならない仕事は、地図をスキャンすることです。地図が大きすぎる場合は、分割してスキャンできますが、各部分について前処理とジオリファレンスの作業を繰り返す必要があることに留意してください。ですから、可能な限り少ない分割で地図をスキャンします。

このマニュアルで提供されるものと異なる地図を使用する場合は、ご自分のスキャナを使い、300 DPIの解像度で地図を画像ファイルとしてスキャンしてください。地図が色付きの場合は、カラーで画像をスキャンし、後でそれらの色を使用して地図からの情報を別々のレイヤ(例えば、林分、等高線、道路...)に分離できるようにしておきます。

この演習では、以前にスキャンした地図を使用します。この地図はデータフォルダ exercise_data/forestryrautjarvi_map.tif として格納されています

14.2.2. basic Follow Along: スキャンした地図をジオリファレンスする

  1. QGISを開き、 プロジェクト ► プロパティ ► CRS で、プロジェクトのCRSを EPSG:3067 - ETRS89 / TM35FIN(E,N) に設定します。これは現在フィンランドで使用されているCRSです。

    ../../../_images/f_1.png
  2. QGISプロジェクトを map_digitizing.qgs として保存します。

QGISのジオリファレンスツールである ジオリファレンサ を使用します。地図をジオリファレンスするには:

  1. レイヤ ► georefRun ジオリファレンサ... でジオリファレンスツールを開きます。

    ../../../_images/georeferencer.png
  2. 地図画像ファイル、 rautjarvi_map.tif をジオリファレンスする画像として追加します: ファイル ► ラスタを開く

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

次に地図をジオリファレンスするための変換設定を定義する必要があります:

  1. 設定 ► 変換を設定 を開きます。

  2. 変換型線形 に、 リサンプリング方法最近傍 に設定します。

  3. 変換先CRS オプションの隣にある setProjection CRSを選択 ボタンを押し、 EPSG:2392 - KKJ / Finland zone 2 CRS を選択します。これはこの地図が作られた1994年当時のフィンランドで使われていたCRSです。

  4. 出力ファイル ボックスの横にあるアイコンをクリックし、フォルダに移動して exercise_data\forestry\digitizing フォルダを作成し、ファイル名を rautjarvi_georef.tif にします。

  5. checkbox 完了後にプロジェクトに読み込む をチェックします

  6. 残りのパラメータはデフォルトのままにします。

    ../../../_images/Clipboard10.png
  7. OK をクリックします。

この地図にいくつかある、座標を示す十字線を使って画像をジオリファレンスします。QGISで通常行うズームやパンのツールを使って、ジオリファレンサ・ウインドウに表示された画像を確認することができます。

  1. 地図の左下隅に拡大して、前述したように KKJ / Finland zone 2 CRSによるXとYの座標ペアを持つ十字線があることに注意してください。この点は地図のジオリファレンスに使用する最初の地上基準点(グランドコントロールポイント)として使用します。

  2. 点を追加 ツールを選択し、(必要に応じてパンとズームして)十字線の交点をクリックします。

  3. 地図座標の入力 ダイアログに、地図上に表示される座標(X: 2557000, Y: 6786000)と、そのCRS(EPSG:2392 - KKJ / Finland zone 2)を記入します

    ../../../_images/fillcoordinates.png
  4. OK をクリックします。

    ジオリファレンスのための最初の座標の準備ができました。

  5. 画像内の他の黒い十字線を探します。それらは北と東の方向にお互いに1000メートル離れています。これらの点の座標は、最初の点との関係で計算できるでしょう。

  6. 画像を縮小して、他の十字線を見つけるまで右か上に移動し、何キロメートル移動したか推定してください。地上基準点はお互いにできるだけ離すようにしてください。

  7. 少なくともあと3つの地上基準点を、最初のものと同じ方法でデジタイズします。次のようなものが出来上がるはずです:

    ../../../_images/Clipboard09.png

すでに3つのデジタイズした地上基準点では、ジオリファレンスエラーを点から出る赤い線として見ることができます。ピクセル単位の誤差は、 GCP テーブルdX (ピクセル)dY (ピクセル) 列にも見ることができます。ピクセル単位の誤差が10ピクセルより大きくなってはいけません。もしそうであれば、デジタイズした点と入力した座標を見直して、何が問題なのかを探す必要があります。上の画像を参考にするとよいでしょう。

基準点が決まったら、後で使用するために保存しておきます:

  1. ファイル ► 名前をつけて保存... に移動します。

  2. フォルダ exercise_data\forestry\digitizing で、ファイルの名前を rautjarvi_map.tif.points にします。

最期に地図をジオリファレンスします:

  1. ファイル ► ジオリファレンスを開始 に進みます。

  2. ジオリファレンスの設定を編集したときに、ファイル名を rautjarvi_georef.tif にしたことに注意してください。

    これでQGISプロジェクトにジオリファレンスされたラスタとして地図が表示されました。ラスタが少し回転しているように見えますが、これは単にデータが KKJ / Finland zone 2 で、プロジェクトが ETRS89 / TM35FIN(E,N) であるためです。

  3. データが正しくジオリファレンスされているかどうかを確認するには、以下の方法があります

    1. exercise_data\forestry フォルダにある、 rautjarvi_aerial.tif という名前の航空写真を開きます。

    2. あなたの地図とこの画像は良く合うはずです: 地図の透明度を50%に設定し、航空写真と比較してみてください。

    ../../../_images/Clipboard14.png
  4. QGISプロジェクトへの変更を保存してください。次のレッスンはこの時点から続きます。

14.2.3. In Conclusion

これで紙地図をジオリファレンスし、QGISで地図レイヤとして利用できるようになりました。

14.2.4. What's Next?

次のレッスンでは、地図上の林分をポリゴンとしてデジタイズし、目録データを追加します。