[OSGeoJapan-discuss] QGISのプロジェクト座標系の設定について

2025年 1月 10日 (金) 20:39:05 PST

吉中さん

確かに、 EPSG:4301 <-> JGD:6668 の表示に関しては変換処理がスキップされますね。

一応、 EPSG(PROJ) 上では tky2jgd.gsb があれば EPSG:6740 で変換するはずですが、
あってもうまくいきませんね。
（ EPSG:6740 は TKY2JGD.gsb しか使っていないので結果は実質 JGD2000 であり、
touhokutaiheiyouoki2011.gsb も用いる EPSG:6714 が選択肢に出てこないのはなぜ。。。
プロセシングツールの「ベクタレイヤを再投影」では選択できますが）

% projinfo -k operation -o PROJ -s EPSG:4301 -t EPSG:6668
Candidate operations found: 2
Note: using '--spatial-test intersects' would bring more results (3)
-------------------------------------
Operation No. 1:

EPSG:6740, Tokyo to JGD2011 (2), 0.2 m, Japan - onshore, excluding northern
prefectures of 'main province' (see remarks).

PROJ string:
+proj=pipeline
  +step +proj=axisswap +order=2,1
  +step +proj=unitconvert +xy_in=deg +xy_out=rad
  +step +proj=hgridshift +grids=tky2jgd.gsb
  +step +proj=unitconvert +xy_in=rad +xy_out=deg
  +step +proj=axisswap +order=2,1

-------------------------------------
Operation No. 2:

unknown id, Ballpark geographic offset from Tokyo to JGD2011, unknown
accuracy, World, has ballpark transformation

PROJ string:
+proj=noop

上記の No.2 の手法により noop (No Operation) が実行されている感じです。
以前は変換できないが、今は EPSG:6668 の towgs84 パラメータ付与によって
楕円体変換レベルの低精度変換は行われると勘違いしていました。すみません


無難なのは、一旦 JGD2000 への変換を挟んだのちに JGD2011 への再投影でしょうか。


2025年1月11日(土) 12:33 yosinaka teruhiko <teruhiko ＠ keroro.keroro.com>:

> Naoki Arita様
>
> 吉中です。今年もよろしくお願いします。
> ご返信および対応策についてありがとうございます。
>
> 経緯をもう少し細かく申しますと以下の通りです。
>
> 旧日本座標系で計測された経度・緯度情報があったとします。
> ・計測された場所は、「とある高校」地点で1999年に計測された旧日本測地系データです（もちろんダミー）。
>
> ・QGIS上で上記地点を読み込む際に、レイヤ座標系を「EPSG4301」で読み込ませてやると、添付PDFのP1赤星のように、想定通りの地点にプロットされます。
> ・仮に、誤ってレイヤ座標系をJGD2011で読み込んだ場合、添付PDFのP1緑星のように、想定とはかなりずれた地点にプロットされる。
>
> ここで、先日投稿させてもらった時の状況なのですが、知り合いは
>
> 「QGIS3.0以降は、オンザフライ投影状態なので、プロジェクト座標系を現在主流のJGD2011で設定したところ、添付PDFのP2のように、旧日本測地系レイヤとJGD2011レイヤが重なってしまった。P3のように、プロジェクト座標系を旧日本座標系にしたところ、これもこれで想定と違う。」
> ということになっておりました。
> JGD2000で投影しておかないと、旧日本座標系レイヤとJDG2011レイヤが混在していると、変換がうまくいかないのを知っていればすぐ気が付けますが（これもザックリとはという注釈付きますが）、慣れてないと厳しいのではないかなと。
>
>
> そこで、普段皆様上級ユーザーの方々は、旧日本座標系とJGD2011が混ざっている場合に、どのように対応されているのかなということをご教示いただきたいという次第でありました。
>
> 確かに、変換うまくいってませんなどの警告が出されても、対応までできるかというとここはここで難しいですが・・・。
>
>
> ＿＿＿＿＿＿＿＿＿＿＿＿＿＿＿＿＿＿＿＿＿＿＿＿＿＿＿＿＿＿＿＿＿＿＿＿＿＿
> 吉中輝彦
> teruhiko ＠ keroro.keroro.com
> (Gmailの場合は送受信エラーが発生するため、teruhikoyoshi ＠ outlook.jp までお願いします。）
>
> ------------------------------
> *差出人:* OSGeoJapan-discuss <osgeojapan-discuss-bounces ＠ lists.osgeo.org> が
> Naoki Arita (OSGeoJapan-discuss 経由) <osgeojapan-discuss ＠ lists.osgeo.org>
> の代理で送信
> *送信日時:* 2025年1月10日 22:18
> *宛先:* osgeojapan-discuss ＠ lists.osgeo.org <
> osgeojapan-discuss ＠ lists.osgeo.org>
> *件名:* Re: [OSGeoJapan-discuss] QGISのプロジェクト座標系の設定について
>
> 吉中さんこんにちは
> 今年もよろしくお願いします。
>
> データの具体的な内容と「意図した場所」というのがどのような場所かわかりませんが、
> まず大前提として、データがどのような座標系であっても「基本的には」適切に表示されます。
> すなわち読み込むデータが緯度経度であっても直角平面座標系であっても UTM であっても、
> 旧測地系であっても、「座標値＋その座標値がどのような座標系か」の情報があれば
> 同じ位置を表すデータは同じ位置で表示されるはずで、表示だけならばデータ座標系は
> どんな座標系であっても問題ありません。
>
>
> ただし、ここで日本の旧測地系 (Tokyo) のようにズレが不規則な場合、変換には
> (1) 日本全域の平均的なズレ情報を使う方法（低精度）
> (2) 各地点のズレを保存されたファイルを利用する方法（高精度）
> のいずれかの方法があります。 Even さんが示しているように (2) は別途ファイルが
> インストール（規定のフォルダにコピー）されている必要があります。
> https://github.com/tohka/JapanGridShift では PROJ (QGIS) で使用可能な形式である
> グリッドファイルを変換してまとめています。
> https://github.com/tohka/JapanGridShift/blob/master/README.ja.md
>
>
> (2) のグリッドファイルが存在しない場合には、 (1) の方法で変換せざるを得ず、
> 低精度でしか変換できないため、旧測地系 (Tokyo) を表示するときは若干の誤差が
> ある状態で表示されます。
> 逆に (2) のグリッドファイルが存在している場合には、高精度で変換されるため、
> 旧測地系 (Tokyo) のデータであっても、正しい位置に表示することができます。
>
>
>
> また近年は解決しておりますが、以前の EPSG データセットに登録されている JGD2011 には
> towgs84 パラメータが付与されていなかった関係で、 JGD2011 への変換が行われない
> 問題がありました。（パラメータがない場合、変換スキップは PROJ の仕様）
>
> https://proj.org/en/5.0/faq.html#why-do-i-get-different-results-with-4-5-0-and-4-6-0
>
> このため、当時は Tokyo -> JGD2011 の変換（表示のための変換含む）においては
> 座標変換処理が行われず、約450m のズレが存在する状態で表示されていました。
> https://groups.google.com/g/qgisshitumon01/c/9afZj2IDzZY/m/mleSVPysAAAJ
>
> https://sakaik.hateblo.jp/entry/20181218/MySQL8013_ST_Transform_first_impression
> https://qiita.com/tohka383/items/43434990e30ad76f6651
>
>
>
> 上記のように、グリッドファイルを用意すれば高精度に変換し適切に表示できますが、
> ある意味各自の環境に依ってしまう部分がありますので、無難な処理としては TKY2JGD.gsb と
> touhokutaiheiyouoki2011.gsb をインストールしたうえで、データを JGD2011 に再投影しておくと
> 確実ではあるかと思います。
> （正確には touhokutaiheiyouoki2011.gsb はあくまで東日本大震災での変動情報であり、
> JGD2000 <-> JGD2011 間の変換パラメータではありませんし、大震災で非常に複雑な
> 地殻変動が起きた地域はパラメータによる変換は適切でないため、変換パラメータを
> 欠損させた箇所も含まれていますが）
>
>
>
> なお「変換ができていない場合」ではありませんが、 QGIS のメニューバーの [設定] - [オプション] -
> [座標変換] - [複数の座標変換を利用できる場合は尋ねる] にチェックがあるか確認してみてください。
>
> ありた
>
> 2025年1月10日(金) 14:16 yosinaka teruhiko (OSGeoJapan-discuss 経由)
> <osgeojapan-discuss ＠ lists.osgeo.org>:
> >
> > English version（I apologize for my basic level of English.）
> >
> > Dear Even-san,
> > Thank you for your reply.
> > I understand that the reason there is no warning when the coordinate
> transformation fails in QGIS is because it is intentionally designed this
> way. Indeed, if there is no transformation method available, it would be
> confusing for the user if they were simply told it's an error.
> > I believe that if something doesn't seem to be working, users can
> resolve the issue by consulting with others in their community. I also
> think that the QGIS forum and this mailing list, where we can discuss and
> receive advice from experienced users and experts, are wonderful resources.
> > Thank you for your continued support.
> >
> >
> > 日本語版
> >
> > Even san
> >  ご返信ありがとうございます。
> > QGISで座標の変換がうまくいってない場合に警告を出さないのは、あえてそのようにしているということなんですね。
> > 確かに、変換方法がない場合に、エラーですと言われてもユーザー側はどうすればいいか困ってしまいます。
> > うまくいってないなと思った場合には、身近なユーザー間での相談で解決となるのかなと思います。
> >
> またQGISのフォーラムでディスカッションして経験者や有識者から助言をいただけるこのメーリングリストも素晴らしいと思ってます。今後ともよろしくお願いします。
> >
> >
> > 吉中輝彦
> > teruhiko ＠ keroro.keroro.com
> > (Gmailの場合は送受信エラーが発生するため、teruhikoyoshi ＠ outlook.jp までお願いします。）
> >
> >
> > ________________________________
> > 差出人: Even Rouault <even.rouault ＠ spatialys.com>
> > 送信日時: 2025年1月10日 11:45
> > 宛先: osgeojapan-discuss ＠ lists.osgeo.org <
> osgeojapan-discuss ＠ lists.osgeo.org>
> > CC: yosinaka teruhiko <teruhiko ＠ keroro.keroro.com>
> > 件名: Re: [OSGeoJapan-discuss] QGISのプロジェクト座標系の設定について
> >
> >
> > 吉中さん、
> >
> > (英語から Google Translate で翻訳。英語のオリジナルメッセージは後ほど)
> >
> > QGIS は、座標変換を実行するために PROJ ソフトウェア ライブラリ (https://proj.org) に依存しています (私は
> PROJ 開発者の 1 人です)。PROJ ソフトウェア ライブラリは、座標参照システムと座標操作に EPSG データベース (
> https://epsg.org) を使用し、測地基準点間の変換を参照します。JGD2000 から JGD2011 への特定のケースでは、EPSG
> データベースには、touhokutaiheiyouoki2011.gsb というグリッド ファイルを参照するレコード (
> https://epsg.org/transformation_6713/JGD2000-to-JGD2011-1.html)
> があります。残念ながら、このグリッド ファイルは、PROJ が使用できると認識しているグリッド ファイルのコレクションである PROJ データ
> リポジトリ (https://github.com/OSGeo/PROJ-data/) にまだ含まれていません。 EPSG
> データセットは、グリッドファイルの名前のみを示しており、ダウンロード元は示していません。また、ダウンロード元がわかっていても、ファイルのライセンスが適切でない可能性があります。そのため、そのグリッドがないと、PROJ、ひいては
> QGIS は、その場合座標操作を実行できません。PROJ データ データベースでグリッドが参照されていない場合、警告は発行されません。これは、グリッド
> ファイルを見つける実用的な方法がない場合 (アクセス制限、有料コンテンツなど)
> があるため、ダウンロードできないファイルに関する警告を発すると、ユーザーを困らせるだけだからです。
> >
> > チケット https://github.com/OSGeo/PROJ-data/issues/36 および
> https://github.com/OSGeo/PROJ-data/pull/37
> では、過去に日本の測地グリッドに関する議論がありましたが、それらを含める上での唯一の障害は、現在のライセンス条件であり、これは明らかにオープン データ
> ライセンスではありません。
> > とはいえ、https://github.com/tohka/JapanGridShift リポジトリの
> https://github.com/tohka/JapanGridShift/raw/refs/heads/master/gsb_files/touhokutaiheiyouoki2011.gsb
> にグリッド ファイルがあります。これをダウンロードして QGIS
> インストールの「proj」サブディレクトリにインストールすると、jp_gsi_gsigeo2011.tif (2011 ジオイド
> ファイルのライセンスを取得するための 1 回限りの成功した取り組みの結果) などの他のグリッド ファイルも見つかります。そうすれば、JGD2000
> から JGD2011 への変換は問題なく機能するはずです。
> > 日本のコミュニティが国土地理院 (GSI)
> と協力してこれらのライセンス問題を解決できれば素晴らしいと思います。理想的には、これらのグリッドは、何らかの申請を必要とする「互換ライセンス」ではなく、クリエイティブ
> コモンズ アトリビューション 4.0 (CC-BY-4.0) の条件に基づいてライセンスされることを望みます。
> >
> > フランスより、よろしくお願いいたします。
> >
> > Even
> >
> >
> >
> > English version:
> >
> >
> > Yosinaka-san,
> >
> > QGIS relies on the PROJ software library (https://proj.org) to perform
> coordinate transformation (I'm one of the PROJ developers). The PROJ
> software library uses the EPSG database (https://epsg.org) for coordinate
> reference system and coordinate operation, which reference transformations
> between geodetic datums. For the particular case of JGD2000 to JGD2011, the
> EPSG database has a record (
> https://epsg.org/transformation_6713/JGD2000-to-JGD2011-1.html) which
> references a grid file called touhokutaiheiyouoki2011.gsb . Unfortunately
> this grid file is not yet included into the PROJ-data repository (
> https://github.com/OSGeo/PROJ-data/), which is the collection of grid
> files PROJ knows it can use. The EPSG dataset just gives the names of the
> grid files, not where to download them from. And even when the location
> where to download them from is known, the file might not be appropriately
> licensed. So without that grid, PROJ, and thus QGIS, cannot perform the
> coordinate operation in that case. They don't issue a warning when a grid
> is not referenced in the PROJ-data database, because, in some cases there
> is no practical way to find the grid file (restricted access, paywalls,
> etc.), so emitting warnings about files that cannot be downloaded would
> just annoy users.
> > In ticket https://github.com/OSGeo/PROJ-data/issues/36 and
> https://github.com/OSGeo/PROJ-data/pull/37, we had a discussion in past
> years about that topic of Japanese geodetic grids, and the only obstacle
> for their inclusion is their current licensing terms, which are not clearly
> an open data license.
> > That said, in the https://github.com/tohka/JapanGridShift repository,
> you can find the grid file at
> https://github.com/tohka/JapanGridShift/raw/refs/heads/master/gsb_files/touhokutaiheiyouoki2011.gsb
> . If you download it and install it in your QGIS installation in the "proj"
> sub-directory, where you can find other grid files such as
> jp_gsi_gsigeo2011.tif (which was the result of a one-time successful effort
> to get a license for the 2011 geoid file), then the JGD2000 to JGD2011
> transformation should work fine.
> > It would be great if the Japanese community managed to resolve with
> Geospatial Information Authority of Japan (GSI) those licensing issues.
> Ideally we would want those grids to be licensed under the Creative Commons
> Attribution 4.0 (CC-BY-4.0) terms, and not a "compatible license" with it
> which requires to fill some application.
> >
> > Best regards from France,
> >
> > Even
> >
> >
> >
> > Le 10/01/2025 à 03:01, yosinaka teruhiko (OSGeoJapan-discuss 経由) a écrit
> :
> >
> > 皆様
> > FOSS4G TOKAIで活動させていただいております吉中です。
> > 本年もよろしくお願いします。
> >
> > QGISのプロジェクト座標系の設定についてです。
> > 旧日本座標系ファイルとJGD2011の座標系のファイルが混在しているプロジェクトがある場合です。
> >
> 先ほど知り合いが、QGISのプロジェクト座標をJGD2011に設定しているのに、旧日本座標系ファイルとJGD2011座標系ファイルが同じ場所に表示されるのはなんでかという質問が私のところに来ました。
> >
> 「旧日本座標系→JGD2011へは直接変換できない」というのが私の認識であり、プロジェクトをJGD2000で表示させればおそらく解決するということで実施してもらったところ、ちゃんと意図した場所に表示されるようになりました。
> > 皆様、旧日本座標系とJGD2011座標系が混在するデータを操るときにどのようにされてますか？
> > ①旧日本座標系をJGD2000形式に変換する
> > ②QGISのプロジェクト座標をJGD2000で作業する。
> > ③その他
> > このあたり、正解はないとは思いますが、自分はこうやると便利だと思うというのがあれば、ぜひ教えてもらいたいと思います。
> >
> > 追伸
> >
> QGISで上記のようにプロジェクト座標系で、座標変換がうまくいかない場合、対応していない場合に「変換ができてない可能性があります」などの警告が出るように設定することはできるのでしょうか？
> >
> > ＿＿＿＿＿＿＿＿＿＿＿＿＿＿＿＿＿＿＿＿＿＿＿＿＿＿＿＿＿＿＿＿＿
> > 吉中輝彦
> > teruhiko ＠ keroro.keroro.com
> > (Gmailの場合は送受信エラーが発生するため、teruhikoyoshi ＠ outlook.jp までお願いします。）
> >
> >
> > _______________________________________________
> > OSGeoJapan-discuss mailing list
> > OSGeoJapan-discuss ＠ lists.osgeo.org
> > https://lists.osgeo.org/mailman/listinfo/osgeojapan-discuss
> >
> > --
> > http://www.spatialys.com
> > My software is free, but my time generally not.
> > Butcher of all kinds of standards, open or closed formats. At the end,
> this is just about bytes.
> >
> > _______________________________________________
> > OSGeoJapan-discuss mailing list
> > OSGeoJapan-discuss ＠ lists.osgeo.org
> > https://lists.osgeo.org/mailman/listinfo/osgeojapan-discuss
>
>
>
> --
> 有田 直樹 <ToNaokiArita ＠ gmail.com>
> _______________________________________________
> OSGeoJapan-discuss mailing list
> OSGeoJapan-discuss ＠ lists.osgeo.org
> https://lists.osgeo.org/mailman/listinfo/osgeojapan-discuss
>


-- 
有田 直樹 <ToNaokiArita ＠ gmail.com>
-------------- next part --------------
HTML$B$NE:IU%U%!%$%k$rJ]4I$7$^$7$?(B...
URL: <http://lists.osgeo.org/pipermail/osgeojapan-discuss/attachments/20250111/755c9c85/attachment-0001.htm>