6.1. Lesson: データを再投影および変換する

ここでもう一度、座標参照系(CRS)の話をしましょう。以前にも少し触れたことがありますが、実用上どのような意味があるのかについては触れていませんでした。

このレッスンの目標: ベクターデータセットの再投影および変換をします。

6.1.1. basic Follow Along: 投影法

現在、すべてのデータおよび地図自体のCRSは、 WGS84 と呼ばれています。これは、データを表現するための非常に一般的な地理座標系(GCS)です。しかし、これから見ていくように問題があります。

  1. 現在の地図を保存してください

  2. 次に、 exercise_data/world/world.qgs の下にある世界地図を開いてください

  3. 拡大 ツールを使って南アフリカを拡大表示します

  4. 画面下の ステータスバー にある 縮尺 フィールドにスケールを設定してみてください。南アフリカの上で、この値を 1:5 000 000 (1対500万)に設定してください。

  5. 縮尺 フィールドを見ながら、マップをパンします

スケールが変化していることに気づいたでしょうか?それは、画面の中心にある 1:5 000 000 でズームした1点から遠ざかっているからです。その点から離れたところでは、縮尺が異なります。

理由を理解するために、地球の球体を考えます。そこには北から南に複数の線が引かれています。これらの経線は、赤道では遠く離れていますが、極で出会っています。

GCSでは、この球体で作業をしているのですが、画面は平らです。球体を平面で表現しようとすると、テニスボールを切り開いて平らにしようとするのと同じような歪みが発生します。これは地図上では、経線は極点(合流するはずの場所)でも同じように離れていることを意味します。つまり、地図上で赤道から遠ざかるにつれて、見えるものの縮尺はどんどん大きくなっていくのです。このことは、現実的には、地図上に一定の縮尺が存在しないことを意味します!

この問題を解決するため、かわりに投影座標系(PCS)を使用してみましょう。PCSでは縮尺変更のための余裕を作り、それを修正する方法でデータを「投影」または変換します。そのため、一定の縮尺を維持するために、PCSを使用するために私たちのデータを投影変換する必要があります。

6.1.2. basic Follow Along: オンザフライ 再投影

デフォルトでは、QGISはデータを「オンザフライ」で再投影します。この意味は、データそのものが別のCRSにある場合でも、QGISはそれを任意のCRSにあるかのように投影することができるということです。

QGISの右下にある projectionEnabled 現在のCRS ボタンをクリックすることで、プロジェクトのCRSを変更することができます。

  1. 表示されたダイアログで、global という単語を フィルタ フィールドに入力します。その下の あらかじめ定義されたCRS フィールドに、いくつかの CRS が表示されるはずです。

  2. WGS 84 / NSIDC EASE-Grid 2.0 Global | EPSG:6933 のエントリをクリックして選択し、 OK をクリックします。

    南アフリカの形状が変化することに注意してください。すべて投影法の変更によって地球の見た目としての形状が変わります。

  3. 前と同様に、1:5 000 000 の縮尺に拡大します。

  4. 地図をパンニングします。

    縮尺は同じままであることに注意します!

「オンザフライ」再投影は異なるCRSのデータセットを組み合わせて使う際にも用いられます。

  1. 南アフリカのデータのみを含む別のベクタレイヤを地図に追加します。これは exercise_data/world/RSA.shp として見つかります。

  2. それを読み込みます。そのCRSを見る簡単な方法は、凡例のレイヤにマウスカーソルを置くことです。それは EPSG:3410 です。

    何に気づきますか?

    レイヤは continents と異なる CRS を持つ場合でも、表示されます。

6.1.3. moderate Follow Along: 他のCRSに設定したデータセットの保存

時には、既存のデータセットを別の CRS でエクスポートする必要があります。次のレッスンで説明するように、レイヤー上で距離計算をする必要がある場合、投影座標系でレイヤを持つことが常に良いことです。

「オンザフライ」再投影は、プロジェクト に関連するものであり、単一のレイヤに関連するものではないことに注意し てください。つまり、レイヤを 正しい 位置で見たとしても、プロジェクトとは異なるCRSを持つことがあります。

レイヤを別のCRSでエクスポートすることが簡単にできます。

  1. training_data.gpkg から buildings データセットを追加します

  2. レイヤ パネルで buildings レイヤを右クリックします

  3. 表示されたメニューから エクスポート ► 新規ファイルに地物を保存... を選択します。すると、 名前をつけてベクタレイヤを保存... ダイアログが表示されます。

  4. ファイル名 フィールドの隣にある ブラウズ ボタンをクリックします

  5. exercise_data/ に移動して、新しいレイヤの名前を buildings_reprojected.shp に指定します。

  6. CRS の値を変更します。最近使用されたCRSのみがドロップダウンメニューに表示されます。ドロップダウンメニューの横にある、 setProjection CRSの選択 ボタンをクリックします。

  7. CRSの選択 ダイアログが表示されます。その Filterフィルタ フィールドで 34S を検索します。

  8. リストから WGS 84 / UTM zone 34S | EPSG:32734 を選択します。

    ../../../_images/CRSselector.png
  9. 他のオプションは変えずにおきます。名前を付けてベクタレイヤを保存... ダイアログは次のようになります:

    ../../../_images/save_vector_dialog.png
  10. OK をクリックします

これでレイヤの新旧の投影を比較すると、2つの異なるCRSでありながら、重なっていることがわかります。

6.1.4. hard Follow Along: 独自の投影法の作成

投影はデフォルトでQGISに含まれるものだけよりも多くあります。自身の投影も作成できます。

  1. 新しい地図を始めます

  2. world/oceans.shp データセットを読み込みます

  3. 設定 ► カスタム投影法... を実行すると、このようなダイアログが表示されます。

    ../../../_images/custom_crs.png
  4. symbologyAdd CRSを追加 ボタンをクリックし、新しい投影法を作成します

  5. 面白い投影法として Van der Grinten I というのがあります。その名前を 名前 フィールドに入力してください。

    他のほとんどの投影がそうであるように、この投影は、長方形のものの代わりに円形フィールドに地球を表します。

  6. 形式 で、 WKT(推奨) を選択します

  7. パラメータ フィールドに次の文字列を追加します:

    PROJCRS["unknown",
        BASEGEOGCRS["unknown",
            DATUM["unknown",
                ELLIPSOID["unknown",6371000,0,
                    LENGTHUNIT["metre",1,
                        ID["EPSG",9001]]]],
            PRIMEM["Greenwich",0,
                ANGLEUNIT["degree",0.0174532925199433],
                ID["EPSG",8901]]],
        CONVERSION["unknown",
            METHOD["Van Der Grinten"],
            PARAMETER["Longitude of natural origin",0,
                ANGLEUNIT["degree",0.0174532925199433],
                ID["EPSG",8802]],
            PARAMETER["False easting",0,
                LENGTHUNIT["metre",1],
                ID["EPSG",8806]],
            PARAMETER["False northing",0,
                LENGTHUNIT["metre",1],
                ID["EPSG",8807]]],
        CS[Cartesian,2],
            AXIS["(E)",east,
                ORDER[1],
                LENGTHUNIT["metre",1,
                    ID["EPSG",9001]]],
            AXIS["(N)",north,
                ORDER[2],
                LENGTHUNIT["metre",1,
                    ID["EPSG",9001]]]]
    
    ../../../_images/new_crs_parameters.png
  8. OK をクリックします

  9. ステータスバーの右側にある projectionEnabled 現在のCRS ボタンをクリックし、プロジェクトの CRS を変更します

  10. 新しく定義した投影法を選択します( フィルタ フィールドでその名前を検索します)

  11. この投影法を適用するため地図は再投影され、こうなります:

    ../../../_images/van_grinten_projection.png

6.1.5. In Conclusion

異なる投影は、異なる目的のために有用です。正しい投影を選択することにより、地図上の地物が正確に表現されていることを確認できます。

6.1.6. Further Reading

このレッスンの Advanced セクションの教材は、 `この記事<https://anitagraser.com/2012/03/18/beautiful-global-projections-adding-custom-projections-to-qgis/>`_から引用したものです。

Coordinate Reference Systems に関する詳細な情報はこちらをご覧ください。

6.1.7. What's Next?

次のレッスンでは、QGISの様々なベクトル解析ツールを使って、ベクトルデータを解析する方法を学びます。