台灣電力座標定位: coortranscheck

轉換檢核

應該可基於測試點來做檢查

而所用的轉換參數，是參考這篇

其中 MOI、OSG1、OSG4 (其他不認得) 得注意，分為

A: MOI 與 OSG4 是用於 proj4js ，上河根據結果應該是用OSG4

MOI的東北角是OSG4 約 0.68m,朝向28.4度(北部)

B: PCTrans 與我的則是用 OSG1的參數(Miller Liu, 2008)

而 A的2.8m，朝向30度則是B

分別是我的、成大、台電在 0.6m內，朝向98度(東南東)

不過奇怪的是，proj4js 查看原始碼發現 towgs84 應該是

m,m,m,asec,asec,asec,ppm

其中 asec 讀入後是一律乘以 SEC_TO_RAD ，換言之讀入的是 asec

故 MOI 與 OSG4 不知道為何記載是這樣

過小的 rxyz 幾乎與 0 無異，sf 過小給下去則等同於1

最終 OSG1 若用在 proj4js ，其中 rxyz 負負得正，則可轉換到跟B一樣的位置

應該是 PCTrans/我的是 PVR，proj4js 是 CFR 的緣故(來源待查)

PVR: Position Vector Rotation
CFR: Coordinate Frame Rotation

參考:

原來以前就已經寫過找最近基石點的方式來轉換

將資料納入頁面只需100kb

平面轉經緯應該公式沒有問題。

經緯轉平面可能得參考 (U)TM Projection 。

https://github.com/jjimenezshaw/Leaflet.UTM

proj4js 平面轉經緯轉平面，結果是 0

我的和 PCTrans 結果一樣??? 97平面 dx = 0.003m, dy = -0.005m

其他:

PCTrans 4.2.8 (自4.2.5版後新增了 Miller Liu 的 TWD67 7參數定義)

台灣電力座標定位