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

[Bravo Astro]

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")#シェープファイル出力

誤りがあれば、ご指摘ください。

よろしくお願いします。


-------------- next part --------------
HTML¤ÎźÉÕ¥Õ¥¡¥¤¥ë¤¬½üµî¤µ¤ì¤Þ¤·¤¿.
URL: http://lists.osgeo.org/pipermail/osgeojapan-discuss/attachments/20100814/ec10509f/attachment.html