• 41

山友分享的全國基石及衛星控制點資料(98.11.18更新檔案)

hwsu315 wrote:
在說明我的想法前,我...(恕刪)


該基石資料datum都是TWD67, 因此沒有你想得這麼複雜! 至於你用GPSr進行檢驗這個作法其實還是偏離我說的觀點, 你的結果是在看接收機的誤差, 因此無法有效去檢定軟體甚至基點的誤差! 除非你拿出測量級GPS進行實測, 不然你只是將這問題更擴大了!

http://www.sunriver.com.tw/grid_tm2.htm#04是虎子山座標, 這才是最精確的數據, 全部都經過測量!

millerliu wrote:
該基石資料datum都是TWD67, 因此沒有你想得這麼複雜! 至於你用GPSr進行檢驗這個作法其實還是偏離我說的觀點, 你的結果是在看接收機的誤差, 因此無法有效去檢定軟體甚至基點的誤差! 除非你拿出測量級GPS進行實測, 不然你只是將這問題更擴大了!

http://www.sunriver.com.tw/grid_tm2.htm#04是虎子山座標, 這才是最精確的數據, 全部都經過測量!


的確,我的想法和您的有些不同,我著重在七千多個TM2座標資料是否是同個Datum,而對於不同的轉換程式,我皆可接受,因為全都在誤差10米內(garmin勉強),您提到資料全都是同個datum,那我就清楚沒問題了。

再來關於您說的哪一種程式的誤差最小,我的想法是找出同一點的精確的 TWD67 座標 和 TWD97 座標,再拿不同程式的轉換來比較。
謝謝您的這個連結,我們可以這個虎子山來試試:
由以上網址所得虎子山 TWD67 TM2: (2652130, 2652336), TWD97 TM2: (248171, 2652130)
用中研院和成大的網頁來轉換虎子山TWD67 TM2的座標到 TWD97 TM2
中研院和成大都是相同的 (248171,2652130)
結果符合預期,因為中研院和成大在台灣中部的轉換結果很相近。

再來就來看台灣其他地區的差別了,因為沒有高精度的GPS,所以目前看來最可靠的資料應該是內政部地政司衛星測量中心的資料。
很可惜的是公布的資料有意或無意不提供最精確的座標,末位數字統一為0,不知是四捨五入還是無條件消去。

內政部的資料,提供了三種不同表示座標的方式:TM2, 經緯 和 XYZ
其中經緯的座標誤差最大,完全無法用
TM2的南北及東西的誤差為+-5米,或+-10米 (看末位0是四捨五入或是無條件消失),勉強可一用
而第三個XYZ,在我計算轉換後發現誤差為TM2座標的1/2到1/3,也就是+-1~3米,我覺得是值得以此為準來探討的

假設山友TWD67 TM2基點資料是正確的

先以七星山為例,
內政部 : TWD97 (121.5534697, 25.1706954) or TM2 (305790, 2784799)

山友TWD67 TM2 (304957, 2785003)
轉換後:
中研院:TWD97 (121.5534236, 25.1706755) or TM2 (305786, 2784797)
成大 :TWD97 (121.5534425, 25.1707037) or TM2 (305788, 2784800)

再選南部一個點:枋寮 原三等三角點 點號537
內政部:TWD97 (120.5850662, 22.9893685) or TM2 (207461, 2543166)

山友TWD67 TM2 (206634, 2543372)

中研院:TWD97 (120.5850903, 22.9893737) or TM2 (207463, 2543167)
成大 :TWD97 (120.5850736, 22.9893523) or TM2 (207461, 2543164)

結果是滿有趣的,山友的這個資料其實滿準確的
不管是成大或中研院,作出來的轉換,都在內政部公布資料誤差內(+- 1~3m)
以這兩點來看,成大似乎稍稍準一點點,但,仍要多取幾個點才比較看得出來。
我用mac沒有mapsource可用,所以懶得列Garmin的結果。

另外,用內政部XYZ算出來的TM2和橢球高來看,我推測內部政末位數的0是四捨五入。
分享我的發現給大家。
hwsu315 wrote:
很可惜的是公布的資料有意或無意不提供最精確的座標,末位數字統一為0,不知是四捨五入還是無條件消去。...(恕刪)


我有取得內政部部份基點的精確座標(2566筆,到小數下三位),在 這邊

以七星山為例,內政部資料:

點名 n e lat lon
七星山 2784799.355 305787.7829 25:10:14 121:33:12

TM2 座標比較準
GPX檔是以TWD97 TM2 轉 經緯度 所製成的



millerliu wrote:
我用溪子口山來對照:

內政部TWD97資料是: 305380, 2765060, 基石資料庫則出現305396, 2765062, 前面的座標值十位數是錯的! 這是比較能肯定的推論!

經緯度部分: 24.59 30, 121.32 50, 基石資料庫則出現 24.59 33.1 , 121.32 55.6 看不出來誤差!


以溪子口這個來看
我從內政部 TWD97 的XYZ資料算得出的 TM2為 (305384, 2765064);

基石資料庫 TWD67 TM2 (304554.669, 2765266.653 )
用成大方法轉換
TWD97 TM2 (305385, 2765063)

東西南北差一公尺,這是另一個兩邊資料符合的例子
溪子口山,內政部資料:

點名 n e lat lon
溪子口山 2765062.499 305385.3217 24:59:33 121:32:55

Wei_1144 wrote:
我有取得內政部部份基點的精確座標(2566筆,到小數下三位),在 這邊


原來您有更精確的座標 :D
資料是從特別的管道取得嗎?
hwsu315 wrote:
原來您有更精確的座標 :D
資料是從特別的管道取得嗎?

不是
Wei_1144 wrote:
我有取得內政部部份基點的精確座標(2566筆,到小數下三位),在 這邊

以七星山為例,內政部資料:

點名 n e lat lon
七星山 2784799.355 305787.7829 25:10:14 121:33:12

TM2 座標比較準
GPX檔是以TWD97 TM2 轉 經緯度 所製成的


我已經轉換成gpi放在Vista HCx裡面了, 實在很驚訝你的超人能力! 可以提供excel檔讓我可以用各定位點的TM2及經緯度來看看GPS機子的表現嗎?

再告訴你一個好消息, 7777大大剛寄給我最新基石資料庫, 居然高達10585筆! 我試著自己轉換看看, 有困難可以向你求援嗎? 哇, 四年來又整理出近3000筆資料了!

我有些不解了, 內政部衛星測量資料, 與聯勤的三角點測量位置一樣, 怎麼高度會不同呢? 水準基點不是都是基隆平均海平面嗎?
關於計算的誤差,可以做個簡單的驗證:

假設所用的參數與算法是完美的,那麼任意給定一個座標,不管轉換多少次,應該還是能還原出原來的座標值

例如:

TWD67 TM2 → TWD67 經緯度 → TWD97 經緯度 → TWD97 TM2 → TWD97 經緯度 → TWD67 經緯度 → TWD67 TM2

總共做六次轉換,第一個 TWD67 TM2 和 最後一個 TWD67 TM2 應該是一樣的

誤差大小就看 參數與算法了



如果在意的是某一個部分的計算誤差,也可以就那部份多做幾次,看看誤差如何
例如: TWD67 經緯度 → TWD97 經緯度 → TWD67 經緯度
多做幾個循環,看看誤差增長的速度

Wei_1144 wrote:
如果在意的是某一個部分的計算誤差,也可以就那部份多做幾次,看看誤差如何
例如: TWD67 經緯度 → TWD97 經緯度 → TWD67 經緯度
多做幾個循環,看看誤差增長的速度
...(恕刪)

我好奇的是Vista HCx實際接收時的誤差, 而不是軟體計算上的轉換誤差! 尤其Garmin新機種的平均功能可以降低接收誤差到什麼程度, 況且接近中央經線與遠離中央經線的定位點數據應該也不一樣才是! 這也是我想要內政部原始數據的原因!
  • 41
內文搜尋
X
評分
評分
複製連結
請輸入您要前往的頁數(1 ~ 41)
Mobile01提醒您
您目前瀏覽的是行動版網頁
是否切換到電腦版網頁呢?