GISGIS実習オープン教材の学習(その3)今回は、課題の第3回として主に座標系の変換を学習します。参考URLは、以下になります。今回の課題は、3つあります。1.空間座標の変換(緯度経度から平面直角座標への変換)2.ジオリファレンス(画像データとベクトル位置データとの関連付け)3.ジオコーディング(背景地図上に位置情報を持つテキストデータの重ね合わせ)早速見ていきましょう。課題1の内容は、以下となります。東京23区コンビニのデータを世界測地系の平面直角座標系へ投影変換し、地図をレイアウトしてください。次に、地理座標系と投影座標系の違いについて、WEBで調べて300字程度でまとめてください。課題3作業手順は以下となります。地図をレイアウト(整飾)地理座標系と投影座標系の違いについて、300字程度でまとめる1.東京23区コンビニデータのダウンロード下記URLからデータをダウンロードしましょう。https://github.com/gis-oer/datasets/raw/master/tasks/tokyo23ku-cvs_task.zipデータは、属性付きのシェープファイルで、座標系は世界測地の緯度経度(JGD2000、EPSG:4612)の点データが3,853件収容されています。データに添付されているREADME.mdには、JGD2011と記載されていますが、ogrinfoでShapeファイルの内容を直接確認するとJGD2000と記述されています。課題の指示は、世界測地系の平面直角座標(EPSG:2451)として先ほどのシェープファイルを読み込みなさいということですので、QGISのデフォルトの測地系を2451にしておいて、ダウンロードしたデータをインポートします。QGISのデフォルト測地系の設定手順は、メニュー[設定]から[オプション]を選択し表示されたダイアログから[CRS]を選択します。そこで、”新しいプロジェクトのCRS”の設定を、”デフォルトのCRSを使用”に指定し、リストからEPSG:2451を選択してOKボタンを押下します。(下図参照)なお、”新しいレイヤのCRS”は、”CRSダイアログを表示”にしておいてください。略で、参照座標系を意味しています。次にダウンロードしたデータをインポートします。手順は、メニュー[レイヤ]から[データソースマネージャ]を起動してダウンロードしたシェープファイルを指定します。エンコーディングは、”UTF-8”としてください。インポートが成功すると地図が表示されますが、赤枠で示した座標が緯度経度ではないことを確認してください。インポートしたシェープファイルには、座標系情報(.prjファイル)が添付されていますので、インポート時に自動的にEPSG:2451に座標変換が行われます。EPSG:2451(平面直角9系)の原点は千葉県野田市ですから、Y座標がすべてマイナスになっていることを確認してください。メニュー[プロジェクト]から[新規プリントレイアウト]を選択して、レイアウトを作成していきます。完成イメージは以下になります。4.地理座標系と投影座標系の違いについて、300字程度でまとめる例えば、以下のような感じで書きます。次元である地球上の位置を緯度と経度で表現する座標系です。地理座標系では、地球の形状(準拠楕円体)や重心によってモデルを定義しています。日本測地系では準拠楕円体にベッセル楕円体を採用していましたが、現在の測量法ではGRS80(WGS84とほぼ同じ)を採用しています。座標で表現する座標系です。日本国内でよく使われている平面直角座標系を例にすると、国内の19ヵ所に基準点を設け、特定の楕円体モデルを前提に平面投影します。その他全世界を表現できる例としては、UTM座標系やWebメルカトル(EPSG:3857)などがあります。まとめ今回は、少し長くなってしまったので、課題1のみとさせてください。課題2をアップします。作業を進めるうえでのTips(お役立ち情報)を以下に示します。1.GDALの導入:以下のURLを参考にするとよいと思います。https://www.kkaneko.jp/tools/win/gisinternals.html2.マークダウンファイル(.md)の表示:Typoraというソフトを使いました。ベータですがWindows版もあります。ogrinfoにてUTF-8ファイルの文字化けを防ぐ方法:コマンドプロンプトにて、以下のコマンドを打ちます。650014.手順2で設定した”新しいプロジェクトのCRS”設定この設定は元に(”最初のレイヤのCRSを使う”)戻した方がよいと思います。5.EPSGコード一覧:以下がよいと思いました。technology.MENU。