PostGISを使って平面直角座標系から世界測地系に変換する方法

こんな感じ。

SELECT AsText(Transform(GeomFromText('POINT(1 -1)', 2443), 4326))

"2443"のところは、平面直角座標系として日本のどこの座標原点を使っているかによって変わってくる。
それぞれのSRIDについてはYaskey GIS基礎知識を参照。

計算結果がの答え合わせは、緯度、経度への換算と比べる？（「度・分・秒」形式なので換算が必要）

PostGISで使う SRID

自分に必要なのはこれだけかな。

select * from spatial_ref_sys where srid in (4301, 4326)

結果を要約すると、

4301 → Tokyo （いわゆる日本測地系。日本の地図でけっこう使われている）
4326 → WGS 84 （いわゆる世界測地系。Google Maps等で使われている）

参考：Bubble://ちずろぐ/別巻/ SRIDList

PostGIS調査メモ

参考になったリンクや、後で見るかもしれないリンクの集。それに加えて気づいたことのメモ。

• PostGISを利用してみる。 - 黒い猫の開発日記
• 初心者には一番分かりやすそう
• 他の記事も参考になるかも → 「postgis」の検索結果 - 黒い猫の開発日記
• PostGISの関数
• PostGIS関数サンプル集
• PostGISで日本測地系(Tokyo)から世界測地系(WGS84)への変換ができない場合
• 自分でも現時点で最新版のPostGIS1.4.0-2で試したが、やはり変換されなかった。Tokyo（4301）のproj4textを見ると下記のようになっていた。変換値がデフォルトでは未設定ということか。
  
  

  +proj=longlat +ellps=bessel +no_defs
  

  
  
• 他のデータを参考にすると、下記のようにするのが正しい？（要検証。）最後に一応他のデータに合わせて半角スペースを入れてある。
  
  

  +proj=longlat +ellps=bessel +towgs84=-147.54,507.26,680.47,0,0,0,0 +no_defs
  

  
  
• projパラメータリスト
• PostGISについて
• PostGISとPostgreSQL幾何データ型の比較
• PostGIS : Home （本家）
  

インストール

• Windowsの PostgreSQL 8.4 / 8.3 にPostGISをインストールする方法
  
• PostgreSQL/PostGISのダウンロードとインストール

Shapeファイル関連

• Shapeファイルのインポート／エクスポート
• Bubble://ちずろぐ/別巻/ shp2pgsql
• shp2pgsqlで"シフトJIS"からの文字コード変換がうまくいかなかったのは仕方ない
• Shapeファイルからのロード時に文字コードとして「WindowsのShift_JIS」を指定する方法
• クローズでないリング発見
• shp2pgsqlのオプションに-wを指定すると人間が緯度経度を比較的見やすいINSERT文が作成される
  
• PostGISへShapefileの投入と表示確認 shp2pgsql pgsql2shp pgAdminⅢ
• シェープファイルについて
• 全国市区町村界データ | 対応データ | ESRIジャパン株式会社
• Shapeファイルの仕様を作ったESRI社の日本法人がShapeファイルを無償配布している
  

マニュアル和訳

• PostGIS 1.3.6 マニュアル 日本語訳
• PostGIS Manual （バージョン不明？）