QGIS 3.2 でJGD2011認識してくれない

[TJ]

ArcGIS で作成した JGD2011 の Shapefile を QGIS 3.2 に追加する、CRS が 「EPSG:4019 - Unknown datum based upon the GRS 1980 ellipsoid - Geographic」 になります。QGIS 3.2 の CRS 一覧には 「EPSG:6668 - JGD2011 - Geographic」ちゃんとあります。毎回手動で CRS を変更しています。
※ QGIS 2.18.22 で確認しましたが、QGIS 3.2 と同じく「EPSG:4019 - Unknown datum based upon the GRS 1980 ellipsoid - Geographic」 になります。

また手動で CRS を JGD 2011 に変更しても、何等かのプロセシングツール（Select by location など）を実行すると、出力したファイルの CRS が 「EPSG:4019 - Unknown datum based upon the GRS 1980 ellipsoid - Geographic」 に戻ってしまいます。これはかなり不便です。

ちなみに、JGD 2000 で定義した Shapefile は CRS は「EPSG:4612 - JGD2000 - Geographic」と正しく認識してくれます。

QGIS で JGD2011 を自動的に認識させる方法があれば教えていただきたいです。

[福岡]

TJ　様

ArcGISを持っていないので、見当違いでしたらご容赦ください。m(__)m

えーと、今回のトラブルですがShapeファイルに原因があるような気がします。

ShapeではCRSに関する定義をprjという拡張子を持ったファイルで行っています。ArcGISで作成したShapeのprjファイルが一般的なJGD2011になっていないのではないでしょうか？

推測ですがJGD2011の定義を手入力で行って、EPSGを4019にしているように思います。

そうするとprjファイルの修正が大事な作業になると思います。

prjファイルを修正するのには、2.18以前では「ベクタ － データマネジメントツール － 現在の投影法の定義」で行います。

レイヤのプロパティでもCRSを修正できますが、実はこの方法ではprjの修正まではしない便宜的なものですので、この方法では修正しないでください。

3.2ではプロセッシングツールの「レイヤの投影法の定義」が同等のツールかと思います。

こちらを使ってShapeファイルを修正するのが良いかと思いました。

[TJ]

福岡様

早速なご返答ありがとうございます。

prjファイルの中身を確認したら、確かに微妙に違っていました。

ArcGISで出力したprj
GEOGCS["GCS_JGD_2011",DATUM["D_JGD_2011",SPHEROID["GRS_1980",6378137.0,298.257222101]],PRIMEM["Greenwich",0.0],UNIT["Degree",0.0174532925199433]]

QGISで出力したprj
GEOGCS["JGD2011",DATUM["D_Japanese_Geodetic_Datum_2011",SPHEROID["GRS_1980",6378137,298.257222101]],PRIMEM["Greenwich",0],UNIT["Degree",0.017453292519943295]]

QGISで出力したprjの場合、QGISに追加Shapeを追加したら、JGD 2011と正しく認識されます。
今後QGISでJGD 2011のShapeを処理する場合、QGISで一回投影法を定義し直してからにします。

本当に助かりました。

[sakuya_izayoi]

こちらizayoiです。

試に、QGIS3.2でJGD2011(EPSG6666)に設定したレイヤーをスクラッチレイヤーで作り、図形を描いてEPSG6666で保存していった結果、このshpを読み込ませたらEPSG3822(TWD97)になってしまう時があります。

どういう時にprjを誤認識するのかちょっと判断できません。

3.2で自分で作ったshpが誤認識されるという事態が起きた方は、他にお見えになりませんでしょうか？

[sakuya_izayoi]

こちらizayoiです。座標系を教えてください。

JGD2000という時は

EPSGコードが4612だとずっと思っていたのですが

QGISにはJGD2000で地理参照系にEPSGコードが6696というのもあります。

またgeocent(地心座標系と訳せばよいかしら？）という参照系ではEPSG:4946とEPSG:4354がありました。

JGD2011にも同じように座標参照系が4つあります。

これらの使い分けって普段どうしていらっしゃいますか？

お恥ずかしい話ですが、JGD2000はEPSG:4612　JGD2011はEPSG:6668

であると思い込んでいたために、これらの知識が欠損しております。分かりやすい説明があればリンク先等をよろしくお願いいたします。

下記画面はJGD2000の場合のQGIS3.2のキャプチャです。

[TJ]

実は同じ疑問を持っていました。

調べましたが、その辺の情報があまりなかったです。

但し、１つヒントを見つけました。QGIS 3.2では、座標系を選択すると、座標系の範囲がダイアログの一番下に赤い塗りつぶしで示してくれます。また、その左側に座標系範囲の情報が表示されます。それぞれ比較してみると範囲の違い分かる。例えば、EPSG:4612は沖縄諸島も含まれているの対し、EPSG:6696には含まれていません。それの違いかなと思います。ただ、EPSG:4946とEPSG:4354の範囲情報は全く一緒だったので、違いが正直分かりません。普段使う分、JGD2000はEPSG:4612、JGD2011はEPSG:6668にすれば間違いないと思います。他の方もし詳細情報があれば共有していただければと思います。

EPSG:4612
Extent: 122.38, 17.09, 157.65, 46.05
Proj4: +proj=longlat +ellps=GRS80 +towgs84=0,0,0,0,0,0,0 +no_defs

EPSG:6696
Extent: 129.30, 30.94, 145.87, 45.54
Proj4: +proj=longlat +ellps=GRS80 +towgs84=0,0,0,0,0,0,0 +vunits=m +no_defs

EPSG:4946
Extent: 122.38, 17.09, 157.65, 46.05
Proj4: +proj=geocent +ellps=GRS80 +units=m +no_defs

EPSG:4354
Extent: 122.38, 17.09, 157.65, 46.05
Proj4: +proj=geocent +ellps=GRS80 +units=m +no_defs

[福岡]

皆様へ

いろいろと検索を繰り返しながら、混乱の真っ最中です。(^_^;)アセアセ

ただ直感的には、4612は「+no_defs」、6696は「+vunits=m+no_defs」となっているので、6696は3次元を視野に入れたものなのか？などと想像しています。

また、空間参照では地心座標系や地心直交座標系というのもあるそうで、3.2で出てくるgeocentという項目にあるEPSG:4946や EPSG:4354はそれに該当しそうな感触です。

[sakuya_izayoi]

TJさま

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

geocent（地心座標系)については全くのノーマークでした。(EPSG:4946とEPSG:4354の事です。)

http://www.geod.jpn.org/web-text/part2/2-1/

のリンク先の2-1-5. 3次元直交座標系

が分かりやすく書いてありました。

そしてEPSG:4354の方ですが2002年に制定されたようですが

https://epsg.io/4354

には推奨されないとあり、EPSG:4946の方が2014年に定義されたみたいですね。

https://epsg.io/4946

地心座標系を使いたい場合はEPSG4946の方が良さそうだという事なのですかね？

ちょっと自信がないため、情報をお待ちしております。

[sakuya_izayoi]

福岡さま

書いている間に更新されていたことに気が付かず、申し訳ないです。

そうですね、3次元空間座標系つまり、空間情報を持つ点群データなどに使用する座標系なのかなという印象を持ちます。

[TJ]

福岡様
izayoi様

いろいろ調査お疲れ様でした。
言われてみると3D表現するために、3次元空間座標系も必要ですね。確かにizayoi様添付のページの定義見るとZ軸の定義も記載されています。私もいろいろ調べてみます。