PostGIS

前言 最近因為工作需求碰到地理空間搜尋與儲存的技術，其中包含 PostGIS 這個 PostgreSQL 實作的地理物件函式庫，裏頭涵蓋各地區座標系統的 SRID(Spatial Reference Identifier)，令人欣慰的是台灣地區常用的 TWD97 與 TWD67 兩種座標系統也涵蓋如此，但在實作過程發現 TWD67 的座標點位有高達數百公尺的誤差，追溯源頭後意外發現，原來是 PostGIS 內的參數設定有誤，解決之後就順手紀錄一下踩坑的過程。 座標參考系格式 - WKT vs proj.4 vs SRID 開始以前要先提一下地理知識，因為地球是圓的不是平的，所以要將一個橢圓形的物件投影到紙面上勢必會有投影扭曲的問題，而為了降低投影扭曲的影響，學術上有很多轉換方式，包含等面積投影、等角投影(麥卡托投影)…，但不管選用哪一種投影都只能降低扭曲的幅度，並不存在絕對好的投影方式。 為了描述座標參考系(CRS, Coordinate Reference System)，訂了幾種常用的描述格式，包含: WKT(Well-Known Text): 是一種描述點、線、多邊形、幾何形狀的文本格式，例如: POINT(121.5 22.5)，相較於其他格式更 human-readable。 舉例來說，Google map 使用的坐標系 WGS84 的 WKT 格式: GEOGCS["GCS_WGS_1984",DATUM["D_WGS_1984",SPHEROID["WGS_1984",6378137,298.257223563]],PRIMEM["Greenwich",0],UNIT["Degree",0.017453292519943295]] 在 PostGIS 中被應用在 spatial_ref_sys 表內的 srtext 欄位中 proj.4: 是另一種描述投影的文本格式，例如: +proj=tmerc +lat_0=0 +lon_0=121 +k=0.9999 +x_0=250000 +y_0=0 +ellps=aust_SA +units=m +no_defs，用幾個參數如橢圓角度、基準、長度單位…來定義投影的方式。 Proj.4 是開源社群常用的地圖投影資料庫，據我所知 PostGIS 和 GDAL/OGR 這些軟體都直接或間接的使用到 Proj.4 資料庫 在 PostGIS 中被應用在 spatial_ref_sys 表內的 proj4text 欄位中 SRID:...