[OSGeoJapan-discuss]日本測地系からWGS 84への変換方法

[Nobusuke Iwasaki]

Bravo Astroさま，みなさま

岩崎です。
QGISの場合，「カスタムCRS」を使って設定することが出来ます。
メニューの「設定」から「カスタムCRS」を選択して下さい。

名称とパラメータが表示されると思いますので，緯度経度座標系であれば，例えば
名称：Tokyo_LatLon
パラメータ：+proj=longlat +ellps=bessel
+towgs84=-146.336,506.832,680.254,0,0,0,0 +no_defs

と入力して，フロッピーディスクマークをクリックして保存してもらえればOKで
す。あとは，旧測地系のファイルを開く場合に，指定してもらえればOKです。
（UTM座標系の場合，またちょっと設定方法が違いますが）

ただし，注意が必要なのですが，この方法では日本測地系と世界測地系の間の正
確な変換は出来ないようです。これは，変換ライブラリであるProj.4が「異なる
楕円体間の変換」の変換をサポートしていないためのようです。たぶんRの変換
もproj.4を使っていると思われますので，同様の問題があるかもしれません。

どなたか，内容について補足や，「正確な」変換が出来るライブラリ等について
情報をご存じでしたら，よろしくお願いします。

#オートデスクさんがなんかそういうライブラリを持っているというお話を聞い
た気もするのですが・・・。


On Sat, 14 Aug 2010 14:58:32 +0900
Bravo Astro <mnara-260116 ＠ kamakuranet.ne.jp> wrote:

> OSGeoの皆様
> Bravo Astroです。QGIS1.5を使用しています。
> 初歩的質問で申し訳ありません。
> 
> 日本測地系シェープファイルからWGS 84シュープファイルへの変換方法を教えてください。
> ここギコさんのHP
> http://kokogiko.net/m/archives/002049.html
> によりますと、単に楕円体をベッセルからWGS84に変換すだけでは不十分で
> 楕円体重心のオフセット
> +towgs84=-146.336,506.832,680.254
> も考慮する必要があるということです。
> この+towgs84をQGISのどこで設定すればよいのかわかりません。
> ちなみに、Rでの変換は次のようにしています。（shift　JISのシェープファイルを使用）
> library(maptools)
> 
> library(rgdal)
> 
> pj<-CRS("+proj=longlat +ellps=bessel +towgs84=-146.336,506.832,680.254 
> +no_defs")
> 
> map<-readShapePoly("kanagawa_murakami.shp",proj4string=pj)#シェープファイルを日本測地系で読み込む
> 
> pj3<-CRS("+proj=longlat +ellps=WGS84 +no_defs")#WGS84を定義する。
> 
> map3<-spTransform(map,pj3)# WGS84へ変換
> 
> writePolyShape(map3,"kanagawa_murakami_wgs84.shp")#シェープファイル出力
> 
> 誤りがあれば、ご指摘ください。
> 
> よろしくお願いします。
> 
>