PostGIS の測地系変換

今ちょっと PostGIS をいじっているんだけど、日本測地系(Tokyo)と WGS 84 の間の変換がどうもうまくいかない。何度やっても、同じ数値が返ってきてしまう。

ということで、いろいろぐぐって調べてみた。どうやら spatial_ref_sys の定義がおかしいとのこと(このへんとかこのへんとかこのへんとか)。見つけたサイトの日付がちょいと古かったんだが、調べてみたら確かにそれが原因みたい。

ということで、

UPDATE spatial_ref_sys 
SET proj4text='+proj=longlat +ellps=bessel towgs84=-147.54,507.26,680.47'
WHERE srid=4301;

としてテーブルを修正して解決。