測量系統97跟98的差別 - 土木

也跳出來補充一下好了(汗)
(以下調整刪文)

※ 引述《AHELF (工作消失中 QQ)》之銘言：
: ※ 引述《biceso (terry)》之銘言：
: : 因為工地用的是97系統
: : 跟要上報檢查的98系統不一樣

不知道業界98系統是什麼. @@"

: 土木工程測量學應該有教的:
: (以下引自WIKI網頁資料)

這個.
http://wiki.osgeo.org/wiki/Taiwan_datums

＊我也是從積丹尼台灣電力座標入門 TWD67 與 97 的轉換的。

: 台灣地區現有轉換參數與精度:
: (TWD67與TWD97相互轉換)專業人員通常以小區域轉換,
: 以取得極高(公分級)的精度, 不以全島為轉換區域.
: 最小曲率法: 官方版台灣本島轉換程式MTWDCON宣告精度為30cm.
: http://www.nlsc.gov.tw/websites/up_page/image/c-images/95-101.gif
: http://www.nlsc.gov.tw/websites/up_page/image/c-images/95-102.gif
: 使用畫面(我找不到下載網頁)
: PS.這就是原文n大推文的網形轉換法

bingo~ 我真的找不到這個程式. /_\ google MTWDCON 下載

: 參考網頁 http://wsserver.epa.gov.tw/gis/coordtrans.aspx
: 行政院環保署有弄成網頁版(不過拿掉地籍坐標部分)
: 二參數轉換: TWD67橫座標=TWD97橫座標-828公尺,
: TWD67縱座標=TWD97縱座標+207公尺,
: 誤差約在5公尺以內.
: (一般人只記這個.但如果業主要求更高精度...)
: PS.這是原文推文g大提的公式

關於這個 97_x+(-828) 97_y+(+207) ... 還得再看這個表
http://wiki.osgeo.org/index.php?title=Taiwan_datums/Test_points&uselang=zh-tw
在
Taiwan datums 最下面的 外部連結 External links
台灣三角點 TWD67 座標表 Taiwan Triangulation Points' Coordinates
最下面的 "另見 See also"
測試點 （Test points）

＊貼到 excel 用資料剖析

: 平面四參數轉換:(土木測量唯一有教的轉換方式)
: 利用兩個已知點(TWD67 和 TWD97都須知道)求解出4個參數
: 再寫出該區域轉換公式
: PS.所以才會各地區參數不同...因為各地區引用的已知點不同
: 利用各地區不同引用點可使誤差降低
: 以下公式僅適用台灣本島，最大誤差約2公尺
: X67=X97-807.8-A*X97-B*Y97; Y67=Y97+248.6-A*Y97-B*X97
: X97=X67+807.8+A*X67+B*Y67; Y97=Y67-248.6+A*Y67+B*X67
: (A=0.00001549, B=0.000006521)

成功大學水工所資地組大地坐標轉換測試程式
http://gis.thl.ncku.edu.tw/coordtrans/coordtrans.aspx

開始講誤差一事：
採用的轉換程式、網頁
中研院 wgs84_tm2 (可批次)
上河文化
成大水工所
PCTrans425 (可批次)

我主要是用上述的 test points 表
假設基準 84 經緯貼在 GoogleEarth 上，
然後用表內的 67tm2 橫x,縱y, (67高 , 應該是水準高H?)
用不同的轉換程式轉成 84經緯，貼回 GoogleEarth 比較。
（表內97xy 轉 wgs84 誤差可忽略不計，各家轉換都幾乎無差的)

以七星山為例。 :P <-- 等一下會說明

wgs84_tm2（單機，可批次）
http://www.ascc.sinica.edu.tw/gis/ISTIS/tools.html
第 10 個, v1.6
= 中研院web轉換（web 版不能輸入高程）
http://webgis.sinica.edu.tw/cgi-tran/webtrans.htm

所使用的三參數 dx dy dz = -752 -358 -179

很令我驚訝的是，上河文化網頁居然也 ≒ 中研院web
http://www.sunriver.com.tw/taiwanmap/grid_tm2_convert.php

以上 (-2.16 , -2.27) （大概。用量的。 :P （被飛踢））
h=0 (-2.24 , -2.21) 偏西北17cm
.....................................↙ 3.12 公尺

後來不知道在哪邊看到的，有人說試著改用新的三參數
wgs84_tm2 三參數改用 -685 -470 -237
*啊想起來了： mobile01 minstrelsy

-> (-0.61 , 0.75)
h=0 (-0.77 , 0.8) 再往西北西偏個 16cm
.....................................↖ 0.93 ~ 1.08 公尺，好很多了?


PCTrans425 採用 bursa-wolf 7para 轉換: (4.2.5 版已內建 TWD67 datum)
功能 2.1 /平面轉經緯 (無關高程)
-> (-0.53 , 0.49) ...................↖ 0.53 公尺
（傳說的 7para 誤差 1m ???)

有趣的是，從公式上看，會先平面轉經緯 67tm2xy -> 67ll(h) -> xyz 7para...
還是可以輸入 67 橢球高，不過表上沒有。 :P
且用 84橢球高代入，後來發現不論84橢球高和67高程代入都一樣 (畢竟才差 20m)

->(-0.67 , 0.60) ...................↖ 0.88 公尺
(多考慮反而差更多? = =a)

成大水工所web
http://gis.thl.ncku.edu.tw/coordtrans/coordtrans.aspx
->(-0.15 , 0.53)
...............北北西 0.55 公尺
（這個四參數不錯呀... )

MTWDCON
-> ( 0 , 0 )
是的。據說這個程式是用內建的三角點來檢核平差，七星山是檢核點。

於是我試著算某地67座標 (某處的電線杆電力座標 :P )
照 test points 表上找出距離該點最近的三角點，採用該點的 dx , dy
97x = 67x - dx
97y = 67y - dy
得到 97xy 後轉 WGS84 ，結果跟 MTWDCON 轉換結果完全一樣... so...


的確，一些論文也提及（以下講述個大概，很多不太記得了）
這些參數轉換法，多少都會產生某方向的偏差或以某處為中心的放射狀偏差
詳見
"TWD67 與TWD97 大地基準轉換方法之研究"
（因為我是不殺狼七參數教徒，所以看到有比七參數好的會腦羞 :P ）
（會因為台灣山區地形起伏之故，不太可能一個參數到底）

所以我在想業界要 67 轉 97 應該還是要用 MTWDCON 來轉吧... @@"
大概就這樣，很多都記不太清楚了，資料或有誤。

--