Garmin 手持機種精準度的種種觀察

millerliu

相信許多人一定也會很在意到底Garmin手持機種是否準確定位, 因此下面分享近幾個月的研究心得:

首先我們將精準度分為兩部分, 一個是機內轉換精確度(以下簡稱"機內誤差"), 也就是在不開機接收GPS衛星訊號的情況下, 觀察一個定點轉換到另一個座標時的誤差大不大, 另一個則是實際接收的精確度(以下簡稱"實測誤差"), 在戶外開機接收GPS後, 與內建定點座標比對後, 觀察顯示座標與實際座標誤差如何. 我們下面將依序談談機內誤差, 與實測誤差.

以上提到的定位點通常是有兩種座標格式, 例如: TWD67, TWD97(WGS84), 這樣的座標經過統計台灣澎湖地區有499個, 其中47個TWD67-TWD97座標誤差超出1m,主要是921大地震或資料有誤,因此還有452個三角點位置與內政部衛星控制點是同一位置(此資料見附件"共同點.zip", 所有的基準轉換都是以這些資料來運算, 因此相互轉換如果有任何誤差, 也就表示是轉換軟體造成的!

說明機內誤差時, 先要理解GPS手持機的資料處理方式, 一般來說機台儲存任何輸入的座標(不管你用哪一種座標格式, 或大地基準), 一律都會存成WGS84經緯度, 格式為121.00000, 25.000(比如:E121度, N25度), 也就是到小數點下五位數. 如果你要得到TWD67二度分帶座標(這是許多官方地形圖及山友常用的格式), Garmin機台便會利用內建的TWD67轉換機制顯示出二度分帶. 當然其他國家使用者選用非WGS84大地基準時也是一樣的運作邏輯.

反過來, 如果你直接輸入TWD67二度分帶座標, Garmin一樣是自動利用機內轉換機制變成WGS84經緯度存起來, 很可惜的是依照許多精確座標資料來比對, Garmin內建TWD67轉換是偏離了10公尺的, 如果你改用user datum(座標系統-->user), 將DX, DY, DY, DA, DF改設為-685, -470, -237, -23, -0.00081204的話, 這個機內誤差平均可以縮小到1.64公尺左右.

因此, 我將機內誤差簡短來說, 如果你是利用TWD97(台灣版的WGS84)的經緯度或二度分帶座標, 那你保證可以做到機內誤差是1公尺的理想狀態, 對於業餘者來說這樣其實是很容易在這一米內找到目標的! 但是若你要使用TWD67的經緯度或二度分帶, 便會受限於機內轉換機制的影響, 前面提到Garmin內建TWD67選項有高達10m誤差, 若使用user datum自行設定參數, 機內誤差平均即可降到1.64公尺.

接下來談談實測誤差, 我們知道機內誤差以WGS84較為精確, 我們便用內政部衛星測量中心的測量點來進行對照, 我在台北市文山區附近的溪子口山, 瑠公公園測量時幾乎都僅有1-2米誤差, 當地對空良好, 因此衛星訊號也極佳, 機台顯示估計誤差都是3m, 加上平均功能重複計算50-300次之後, 可說是已經達到民用機種的最精確定位表現了, 但是換到景美橋這個水準點時, 由於東側有10層樓遮住, 僅有一側才能接收訊號, 前面方法都用過了, 誤差依然高達10米以上, 因此要降低實測誤差只有找到天空開闊, 衛星分佈在45度以內, 45度以外時, 定位精準性才會提高!


如果這篇文章對你有益, 給點積分吧, 好久沒有吃大補丸了!!


========更新=============
台灣地區:
網友以程式驗證後發現user datum最佳轉換參數
DX, DY, DY, DA, DF為-685, -470, -237, -23, -0.00081204
這樣機內誤差平均可以縮小到1.64公尺左右.

澎湖地區:
網友以程式驗證後發現user datum最佳轉換參數
DX, DY, DY, DA, DF為-752, -349, -179, -23, -0.00081204
這樣機內誤差平均是 2.44公尺, 最大誤差是4.09公尺

同時也請一併設定user Grid(使用者自訂方格)讓澎湖金馬地區可以顯示正確的二度分帶座標格式: 中央經線--E119度, 尺度比:0.9999, 橫座標平移量:250000, 縱座標平移量: 0

順帶提醒: 金門馬祖地區由於缺乏檢核點, 無法推算轉換參數, 若以上述參數進行轉換, 誤差將高達(500, 160)m, 建議以WGS84進行定位, 不管是TM2或經緯度均可.

========6/10更新===========
近幾個月以共同點檢核座標精確時發現, 使用TWD97二度分帶來觀察誤差還是有發生錯誤的四捨五入的情形, 例如: xxxx00.56, yyyyy00.63座標居然會變成 xxxx00.0, yyyyy00.0(自動捨去!), 因此還是以經緯度(以度為單位)才是最精確的檢核方式, 如: 24.99253258, 121.54865689是溪子山(仙跡岩)精確WGS84座標, GPS應該顯示24.99253, 121.54866才是誤差最小的結果.

http://www.itc.nl/ilwis/downloads/tools/geodeticTools.asp  及 http://earth-info.nga.mil/GandG/geotrans/index.html 是與此篇轉換模式Abridged Molodensky Transformation有關的程式.
尤其第一個連結中的Inverse Molodensky軟體, 算出來的DX, DY, DZ居然很適合, 我又有懶惰的理由!

[ 本文最後由 millerliu 於 2008-6-10 13:54 編輯 ]

millerliu

前述附件檔可以讓大家檢驗我提供的DX, DY, DZ, DA, DZ是否較Garmin精確, 目前我正在利用"共同點.zip"的金門馬祖地區資料繼續計算轉換參數,  有興趣者不妨將此當作研究專題吧! 上述五個參數就是Abridged Molodensky transformation中的, 很容易找到轉換公式的! www.gmat.unsw.edu.au/snap/gps/clynch_pdfs/Datum_ii.pdf  中7-35兩頁就是.

台灣, 澎湖地區參數是Mobile01網友協助計算的, 目前建議給Garmin 何經理, 不清楚他們何時願意接受這套參數, 文件中有建議文字, 但是尚未得到回應!

[ 本文最後由 millerliu 於 2008-6-3 16:37 編輯 ]

millerliu

原文由 jason204 於 2008-6-6 17:28 發表

就我所知，金門及馬祖地區官方以往所用的坐標系統都不是TWD67坐標系統，金門地區用的是單純的平面坐標系統（X,Y），而馬祖地區則是南京坐標系統，採用的投影方法為UTM－6度分帶投影，投影分帶為50及51分帶， ...

http://www.geo.ntnu.edu.tw/facul ... rojection/index.htm 有日治時期至今的四種座標系統, 因此我相信沒有理由南京座標系統(即Garmin內建的Hu-Tzu-Shan系統)是一直使用到TWD97出現才停用. 此外金馬地區地形圖是管制品, 也讓我們對此保持模糊狀態, 但是官方程式卻是可以在外島轉換為TWD67座標, 因此我的說法我存疑! 下面是我找到的文字:

"聯勤401廠於民國88年起，為配合GPS之運用，於軍圖GRS-67系統下之UTM方格網，加印WGS-84之方格網（以紫色短線加印於圖幅邊上），此種軍圖目前已撥發部隊使用。使用者可同時運用兩種座標系，檢查PLGR內建之WGS-84座標系與USER1轉換GRS-67座標系之精度。"

至於您提到共同點剔除部分, 其實雖然TWD67共同點是民國69年的老資料, 但是你一定沒想到這些TWD67座標與近幾年公布的TWD97座標相差無幾, 這些是用官方轉換程式驗證過的, 超出1米以上誤差的共同點被我刪除, 剩下的450幾點就是我進行轉換參數的基本資料.

七參數(Molodensky-Badeka, 或Bursa Wolfe)雖然較為精準但是那是測量級等高階GPS內建的轉換法, Garmin民用機僅採用Abridged Molodensky Transformation, 如果可以將全部450個計算出的DX, DY, DZ參數算出來, 你可以隨意計算任何區域的轉換參數, 例如: 我上面提供的是適用全台灣的, 但是你能找出北台灣, 中台灣, 南台灣的DX, DY, DZ, 這樣轉換參數會比台灣全島用的更精確, 只是目前我還找不到如何用複回歸(Multiple Regression)方式進行最後階段的計算, 需要高手指導, 前面的參數是Mobile01網友算出來的, 可是我卻無法繼續問出具體計算步驟.

其間還將共同點座標換算成X, Y, Z三維座標來計算DX, DY, DZ, 結果發現僅能進行極微小區域的轉換, 如果用在全台灣地區誤差竟然有幾十公尺, 只有該共同點附近才能精確轉換. Jason, 可以教教我嗎?

[ 本文最後由 millerliu 於 2008-6-10 11:12 編輯 ]

millerliu

原文由 chihhsiungchen 於 2008-6-10 11:34 發表
手持機是否準確定位，和手持機接收和運算的硬體有關係。
和大地座標系統沒有關係。

只是現在國際標準，都是用WGS84溝通。（這是沒辦法的事）
要用虎子山（TWD67）地理座標系統，所造成的誤差那是因為轉換資 ...

mobile01網友算出什麼？我不知道。（沒有連結）
但是mobile01充斥太多人云亦云，或半路出師抄襲別人論文的假知識，不可不慎。那些轉換運算，真正精準到一個程度的校正，早有做出研究報告了。

你手中已經有Garmin GPS或相關軟體進行驗證了嗎? 是否精確無須用這樣以偏該全的結論來表達, 這與mobile01無關.  

轉換精準的研究報告有以Abriged Molodensky Transformaton模式進行嗎? 許多論文我都去翻閱過, 多數都是七參數, 最小曲面法, 最小二乘法等專業模式, 但是沒有適用Garmin gps可以採用的轉換參數--DX, DY, DZ.

另外, 手持機是否準確定位這個問題, 請你閱讀一下全文理解後再下結論不遲, 只要用邏輯來討論都歡迎, 你的用語與我也有差異!

millerliu

原文由 chihhsiungchen 於 2008-6-10 14:44 發表
抱歉，我要對樓主您說一下：
您把我最重要的一句話遺漏引用：
....料而來，初學者先不要被誤導了。
其他的話。。。。。
抱歉，既然您說我以偏概(該？)全，那也不需歡迎我討論了。

願意具體來討論這樣的技術問題都是我歡迎的對象!

millerliu

原文由 jason204 於 2008-6-10 21:50 發表
3.對於Abridged Molodensky Transformation不熟，說真的還是第一次聽到，個人比較常用的是Molodensky-7參數轉換(基準轉換)及經過2°TM地圖投影後的平面4參數Helmert轉換與6參數Affine轉換，不過我都是檢現成的程式在RUN，因為沒有寫電腦程式的經驗(其實是不會啦)，複雜的轉換公式也早就還給老師了。

4.對於TWD67座標與近幾年公布的TWD97座標轉換與共同點剔除，大概是應用及精度上的差異，個人的經驗因採用全區轉換不符合工作上的精度需求，所以大多只能用分區域的坐標轉換方法。

5.對於共同點座標換算成X, Y, Z三維座標來計算DX, DY, DZ, 結果發現僅能進行極微小區域的轉換，意思不太了解。但在該共同點附近才能精確轉換的部分，我想是否是因為共同點的二組坐標不符值是會被待轉換點吸收掉，.............，所以對於共同點的篩選與適度的剔除，在要取得較高精度的坐標轉換成果之實務上是相當重要的。

台灣光復初期直到58年間確實是用Hu-Tzu-Shan系統, 也就是你說的南京系統, 使用橢球1924, 現在還是有山友為了尋找當時官方記載的三角點還是用Hu-Tzu-Shan經緯度資料找出不少遺棄山林中的基石.

看來你應該是專業測量者, 才會整天使用七參數轉換, Abridged http://www.itc.nl/ilwis/downloads/tools/geodeticTools.asp  是Abridged Molodensky Transformation計算程式--Inverse Molodensky軟體可以直接算出來的DX, DY, DZ(經緯度以秒為單位輸入, DX, DY, DZ計算到以cm為單位), 對於Garmin使用者來說實在是很方便. 我用溪子山座標計算出來的向你說明好了.

我用上述座標以Inverse Molodensky軟體算出 -735, -388, -193, 經過台灣南北幾個座標來確認轉換TWD67誤差時, 發現這組參數非常精確, 甚至不亞於網友先前幫我算的 -685, -470, -237.   

興沖沖地接著計算澎湖地區得出-557, -274, -187, 似乎誤差就挺離譜的, 不理解這是怎麼回事? 應用在台灣地區OK, 用在離島居然荒腔走板! 看來可能還是要用公式手工計算一次來必對一下, 沒得偷懶!

至於我先前將經緯度轉換成卡式XYX座標, 將TWD67, TWD97的XYZ座標相減直接算出的DX, DY, DZ是(溪子山) -740, -381, -189(以上高度為0), -750, -365, -180(以上高度依照原始資料輸入), 但是僅能讓台北市南區座標精確, 但是離開這範圍誤差便非常大, 比起前兩段算出的參數都不理想, 根本還無法用在台北縣市區域呢? 正是因為如此才發現Abridge Molodensky Transformation才是Garmin使用的轉換模式, 目前驗證似乎就是如此, 就看如何找出正確的參數了!

[ 本文最後由 millerliu 於 2008-6-11 09:16 編輯 ]

millerliu

http://www.gartrip.de/long.htm  網路上類似我這樣的無聊業餘者還真不少, 這篇文章還比較了2000年前後有無SA的差異, 及有無EGNOS差異, 及不同接收時間(5min, 20min, 1hr, 3hr, 6hr, 12hr, 24hr, 36hr)的精度比較. 德國人的研究精神實在讓人佩服!