国土数値情報の世界測地系データの変換

[TDS OKNYM]

質問がございます。

JGD2011EPSG:6698　という座標系で作成しているマップに国土数値情報の神奈川の土砂災害警戒区域のデータをオーバレイさせたいのですが、なかなか、まったくうまくいきません。

どのようにすれば変換されオーバレイできますでしょうか。

ちなみに国土数値情報の土砂災害データは下記のURLになります

https://nlftp.mlit.go.jp/ksj/gml/datalist/KsjTmplt-A33-v1_4.html

[福岡]

TDS OKNYM 様

こんにちは

QGISは、データが正しければ自動的に重ねてくれるようになっています。
そしてデータがうまく重ならない場合は、CRSが誤って設定されていることが原因となっているケースが多いです。

このような投稿があった場合には、Webマップにデータを重ねて確認することをお勧めします。
CRSに関する簡単な説明とWebマップの利用について、Qiitaに記事を書いていますので、そちらをご覧くださればと思います。

参照座標系のざっくりとしたお話し
https://qiita.com/fukuoka1960/items/2318898e47ec1e749e1c

事前に作成していたマップについても、ご確認ください。

ダウンロードしたデータが重ならないということなので、直感的には事前に作成していたマップの方に不具合があるような気がします。

2022年2月22日火曜日 11:20:48 UTC+9 TDS OKNYM:

[TDS OKNYM]

元の顧客から受領したデータ自体がJGDで座標が、(-63000,-94000)というふうな座標（神奈川西部の某市町村）で、上記の土砂災害データは世界測地系で座標が、(139.000,36.000)というふうな緯度経度のような座標になります。

正規の方法では、福岡先生のおっしゃる通りにするべきなのでしょうけど、状況的にかなり作業をしてしまい、工期も近いため最小限の手順で解決したいと思っています。

そのような虫のいい方法はないでしょうか。

2022年2月22日火曜日 11:54:38 UTC+9 福岡:

[福岡]

TDS OKNYM 様

こんにちは

受領したデータのお示しくださった値から推測すると、CRSは「JGD2011/Japan Plane Rectangular CS IX（EPSG:6677）」ではないかと思います。
レイヤプロパティの「ソース」タブの「設定されたCRS」を変更してみてください。

それでうまく行った場合ですが、新規プロジェクトに受領ファイルを読み込んでみてください。
レイヤのプロパティからCRSを確認して「JGD2011/Japan Plane Rectangular CS IX（EPSG:6677）」なら、強引に誤ったCRSで読み込んでいます。
その際は、上記のとおりの「設定されたCRS」の変更でOKです。

もしJGD2011（EPSG:6668）であった場合には、受領ファイル修正までを行うのが根本的な対処となります。
念のために、私なら行う対処方法を記載します。データセットの提出まで求められるのなら、ここまでやった方が良いかと思います。

プリント出力のみを求められている場合などは、成果提出後に試してみてください。

【データ形式がShapeファイル】
プロセッシングツールの「ベクタ一般 － シェープファイルの投影法の定義」でCRSを再設定できます。

【データ形式がShapeファイル以外】
レイヤパネルの対象レイヤ上で右クリックして「エクスポート － 新規ファイルに地物を保存」で正しいファイルを書き出します。

2022年2月22日火曜日 14:20:21 UTC+9 TDS OKNYM:

[福岡]

すみません。追記です。

受領ファイルがラスタである可能性について、失念していました。

ラスタの場合は、メニューの「ラスタ － 投影法 － 投影法の割り当て」を使って、修正してください。

2022年2月22日火曜日 16:27:31 UTC+9 福岡:

[TDS OKNYM]

両データの隔たりが非常に大きく何か根本的におかしいように思います。

片方は日本の緯度経度3桁2桁の座標になっており、もう片方は桁が5桁の座標です。

自分の理解が悪いと思いますが、JGD2011（EPSG:6668）をJGD2011（EPSG:6677）に変換するのではほとんど変化はありません。

分かりやすい話、地球と月位離れている感じなのです。

再度ご教示いただけますでしょうか。

ちなみに全てシェイプファイルで、ベクターデータです。

2022年2月22日火曜日 18:47:49 UTC+9 福岡:

[TDS OKNYM]

補足いたします。拙い説明でわかりにくかったと思います。

つまり、平面直角座標系で設定されているプロジェクトファイルに対して、土砂災害データの緯度経度座標のシェイプファイルを重ね合わせたいということを

伝えたかったです。

2022年2月22日火曜日 19:47:14 UTC+9 TDS OKNYM:

[ありた]

こんにちは

両レイヤのレイヤプロパティを開き、情報タブを確認してください。

そのうち「座標参照系」と「領域」の値がそれぞれどのように表示されているかを

記載してください。

福岡さんがおっしゃってるように、本来の座標系に対し、それを読み取るための

座標系が適切に設定されていれば GIS 上で正しい位置に表示され、重なり合います。

おそらくその設定が不適切な状態であると予想されますが、具体的な情報がないため

予測しながらの回答となっているため、何度もやりとりが生じてしまっています。

データの状況についてなど、提示できる範囲でもっと提示するようにしてください。

[福岡]

ありた 様

フォローしてくださって、ありがとうございます。m(__)m
１人のみの投稿よりも説得力が違うので、大変助かります。

TDS OKNYM 様

ありたさんが補足してくださったように、現在のQGISはプロジェクトのCRSに合わせて、自動的に各レイヤを表示してくれる仕様になっています。
例えば各レイヤのCRSが地理座標と投影座標が混在していても、プロジェクトのCRSに合わせた見た目で重ね合わせてくれます。

つまり重ならないということは、レイヤのCRSが適切に設定されていないことが原因です。

例えばですが、土砂災害警戒区域でダウンロードしたファイルは３種類あると思うのですが、どれをお使いでしょうか？
GeoJSONをお使いの場合は、デフォルトではプロジェクトのCRS（この場合は平面直角座標9系）が適用されると思います。
正しくはJGD2011（EPSG:6668、緯度経度）ですので、手作業で設定してあげる必要が生じます。

上記のケースである可能性は結構あると思いますが、その他にも色々なケースが考えられます。
そこで私は自己判断が容易にできるWebマップとの重ね合わせをお勧めし、ありたさんは適切な情報提供をお願いしているわけです。

CRSの設定でお困りのようなので、おそらくはGISを使い始めたばかりかと推察いたします（間違っていたらごめんなさい）。
お手数をお掛けしますが、よろしくお願いいたします。

2022年2月22日火曜日 21:38:56 UTC+9 ありた:

[TDS OKNYM]

ありた様、福岡様、質問の仕方が適切でない点、理解度が低い点、申し訳ないです。

また、ご教示いただき誠にありがとうございます。

まずは、お詫びとお礼を述べさせていただきました。

Webマップというものがどこにあるのか探し見見ます。

15年くらい前に2年位ArcGISを使って研究していたことがあり、QGISは始めて3ヶ月位です。不勉強でご迷惑おかけしております。

ありた様

下記にデータのソースを貼り付けました。

【土砂災害警戒区域データの情報】

名前
A33-19_14Polygon
パス
W:\***********\A33-19_14_GML\A33-19_14Polygon.shp　　　（****は顧客名が入るので伏せ字にしました）
ストレージ
ESRI Shapefile
コメント

文字コード
Shift_JIS
ジオメトリ
Polygon (MultiPolygon)
座標参照系（CRS）
EPSG:6668 - JGD2011 - 地理学
領域
138.9841368923265463,35.1314768256272387 : 139.7439965655160847,35.6609433741549608
単位
度
地物数
49,518

【業務データ】

名前
***********
パス
W:\21年度\ ***********  \ ***********  .管渠.矢印線.shp
ストレージ
ESRI Shapefile
コメント

文字コード
Shift_JIS
ジオメトリ
Line (MultiLineString)
座標参照系（CRS）

領域
-69848.2260000030510128,-94548.9240007079206407 : -64066.9690000026894268,-92635.3350004213862121
単位

地物数
34

2022年2月23日水曜日 8:48:24 UTC+9 福岡:

[福岡]

TDS OKNYM 様

ご返信ありがとうございます。

Webマップですが、私の最初の投稿でリンクを張っておいた記事に、プラグインを使った表示方法を説明しています。
もう１度リンクを張っておきますので「５.Webマップに重ねてみよう」をご覧ください。

参照座標系のざっくりとしたお話し
https://qiita.com/fukuoka1960/items/2318898e47ec1e749e1c

さてご提示頂いた情報を拝見しました。
貼り付けミスでなければ「業務データ」のCRSが空欄になっているのが気になります。
領域の数値からすると、CRSは「JGD2011/Japan Plane Rectangular CS IX（EPSG:6677）」のようです。
プロセッシングールの「ベクタ一般 － シェープファイルの投影法の定義」でこのCRSを指定して実行すると解決するかもしれません。

念のため、ファイルのバックアップを取ってからお試しください。

2022年2月23日水曜日 18:17:41 UTC+9 TDS OKNYM:

[ありた]

こんにちは

情報のご提示ありがとうございます。

福岡さんも回答されていますが、少し補足いたします。

提示していただいた情報を見ると、業務データ（管渠矢印線データ）は、読み取るための

座標系情報が正しく設定されていません。

レイヤプロパティのソースタブをみると「不正な投影法」となっていると思います。

ESRI Shapefile 形式は、位置情報や属性情報が別々のファイルで保存された形式です。

そして座標系情報は *.prj ファイルとして用意されていれば、その設定が適用されます。

もし *.prj ファイルが同梱していないデータを QGIS で開こうとしたときには、座標情報を

どの座標系として認識したらよいのかわからないため、本来の位置に描画することができません。

（正確には、データの座標系情報が不明な場合、設定しない（未知の座標系／不正な座標系）、

　確認ダイアログを開く、プロジェクト CRS にあわせる、デフォルト CRS を適用する。の

　いずれかの動作となります。この動作は QGIS の設定により変えることができます。）

レイヤプロパティのソースタブで適切な座標系を指定することで、正しく読み取ることが

できるようになります。

（ただし、本プロジェクト内の設定であるため、再度読み込む際には不明なまま）

プロセシングツールの「シェープファイルの投影法の定義」であれば、レイヤの設定のほか

*.prj も作成してくれるようなので、次回以降は（今回定義した）座標系として読み取ることが

できるようになります。

[TDS OKNYM]

福岡様

ご回答ありがとうございます。

ご教授いただいたやりかたで、重ね合わせることができました。

大変ありがとうございました。

参照座標系のざっくりとしたお話しですが、

ちなみに、手順通り進めてみたのですが添付画面のところで、使用できるWebサービスの選択肢にGoogleやその他のサービスが入ってきていません。

どのようにすれば、あるべき姿になりますでしょうか。

2022年2月23日水曜日 20:06:32 UTC+9 福岡:

[TDS OKNYM]

ありた様

ご回答ありがとうございます。

＞レイヤプロパティのソースタブをみると「不正な投影法」となっていると思います。

たしかにそのとおりでした。

ご教授いただいたとおりプロセシングツールの「シェープファイルの投影法の定義」で重ね合わせることができました。

今はテスト用のフォルダで業務データを1つに絞ってテストしましたが

このあと、実際の業務プロジェクトで、使っているシェイプに同じ作業をしてまいります。

手順と考え方を教えていただきありがとうございました。

2022年2月23日水曜日 23:27:14 UTC+9 ありた:

[福岡]

TDS OKNYM 様

こんにちは。
お示しくださった画像で説明いたします。

上の方にタブが並んでいることが確認できると思います。
画像では「Visiblity」タブの内容を表示しています。

右隣の「More services」タブの上でクリックして、このタブへ移動してください。
移動したら、下の方に「Get contributed pack」というボタンがありますので、これをクリックします。

contributed packの取得に成功したら「Visiblity」タブをクリックし、このタブへ戻ると選択できるWebマップが増えているはずです。

2022年2月24日木曜日 10:06:37 UTC+9 TDS OKNYM:

[TDS OKNYM]

福岡様

こんにちは。ご回答ありがとうございます。

Webマップを増やすことができました。

ありがとうございました。

2022年2月24日木曜日 11:43:14 UTC+9 福岡: