わが国では，次の空間参照系がよく使われます．

日本測地系は2002年3月31日まで使用されていたものです．世界測地系は現在使用されているものです．WGS84は米軍がGPS用に管理している測地系です．

世界測地系とWGS84とはあまり差がなく，メートル単位の誤差を許容する場合には，そのまま使用しても差し支えない程度です．

世界測地系と日本測地系は，投影変換を要します．

UTMは，1万分の1 から 20万分の1 の地形図に用いられます．ゾーンは地球全体を経度6度ごとにゾーンを定めています．よってゾーン番号は経度で決まります．ゾーン番号と経度範囲の関係は次の通りです．

ゾーンをまたぐ場合は，通常は，より広くインタセクトしている方のゾーンを用います．

平面直角座標系は，1万分の1 よりも大縮尺(詳細)の地図に用いられます．国内陸域を19個の系に分けています．概ね1つの系に対して複数の都府県が対応し，都府県が複数の系に分かれることは少ないですが，離島部を持つ都県は複数の系をまたぎますし，北海道は3つの系に分かれています．

詳細については平面直角座標系(平成14年国土交通省告示第9号)(国土地理院サイト)をご覧ください．

また，地図関連情報ガジェット(Finds.jp)にて，系番号を簡易的に確認できます．

ご注意: PostGISで投影変換を行う場合には，空間データベースにあるspatial_ref_sysを更新しなければなりません．

# Tokyo
<4301> +proj=longlat +ellps=bessel +no_defs  <>

なお，+no_defsは，デフォルトファイルを読み込まないことを指示しています．このエントリの記述通りに動作することを保障するために必要です．

<>(終端文字列)より前に，たとえば，+towgs84=-146.336,506.832,680.254,0,0,0,0を挿入します． 次のようになります．

	パラメータには，いくつかのバリエーションがあります．

# Tokyo
<4301> +proj=longlat +ellps=bessel +towgs84=-146.336,506.832,680.254,0,0,0,0 +no_defs  <>

4301のほか，30161-30179, 3092-3096についても同じように改訂することで，数メートル程度の誤差に収まるようになります．ただし，なお誤差を含みますので，十分にご注意ください．

spatial_ref_sysの中身を見てみましょう．

db=# select * from spatial_ref_sys where srid=4301;
 srid | auth_name | auth_srid |                                                                                                                                srtext                                                                                                                                |               proj4text               
------+-----------+-----------+----------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------+---------------------------------------
 4301 | EPSG      |      4301 | GEOGCS["Tokyo",DATUM["Tokyo",SPHEROID["Bessel 1841",6377397.155,299.1528128,AUTHORITY["EPSG","7004"]],AUTHORITY["EPSG","6301"]],PRIMEM["Greenwich",0,AUTHORITY["EPSG","8901"]],UNIT["degree",0.01745329251994328,AUTHORITY["EPSG","9122"]],AUTHORITY["EPSG","4301"]] | +proj=longlat +ellps=bessel +no_defs 
(1 row)

このうち，問題になるのは，proj4textカラムです．+proj=longlat +ellps=bessel +no_defsとなっていて，Proj.4と同じくtowgs84パラメータを持っていません．このままで実行すると，経度緯度が変化しません．

select ST_AsEWKT(ST_Transform('SRID=4612;POINT(135 35)'::geometry,4301));
        st_asewkt        
-------------------------
 SRID=4301;POINT(135 35)
(1 row)

proj4textを更新すると，変換結果が変わります．

db=# UPDATE spatial_ref_sys SET proj4text='+proj=longlat +ellps=bessel +towgs84=-146.336,506.832,680.254,0,0,0,0 +no_defs' WHERE srid=4301;
UPDATE 1

db=# select ST_AsEWKT(ST_Transform('SRID=4612;POINT(135 35)'::geometry,4301));
                     st_asewkt                      
----------------------------------------------------
 SRID=4301;POINT(135.002792603664 34.9968021008188)
(1 row)

	パラメータには，いくつかのバリエーションがあります．