• 41

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

我也來貼圖

成大轉換


86年內政部七參數轉換


米轉換


米轉換內政部參數


手持機GPS中研院三參數


手持機GPS Garmin三參數
我再貼

以成大轉換的差異大小排列,將差異大於1公尺的基點資料先刪掉,剩下373筆資料
再計算 ((距離差異5 米轉換內政部參數) - (距離差異4 米轉換)),並檢視其分布圖,期望能看出不同參數的影響
這個很奇怪,同樣的參數,但差異很大,應該是算法不一樣

(距離差異3 86年內政部七參數轉換) - (距離差異5 米轉換內政部參數)


與樓上的分佈圖比較,台灣東、西部點位的"差異的差異" 有同樣的傾向(圖形看起來往右傾斜),會不會是算法上的系統性差異?

如果是的話,可能要先修正算法,再求改善參數
Wei_1144 wrote:
(距離差異3 86年內政部七參數轉換)


我一直沒有搞清楚, 這個86年內政部七參數轉換是有現成的轉換程式嗎? 還是是以內政部的參數套用哪裡的程式?

Wei_1144 wrote:
與樓上的分佈圖比較,台灣東、西部點位的"差異的差異" 有同樣的傾向(圖形看起來往右傾斜),會不會是算法上的系統性誤差?


我覺得可能是算法的細部不同所造成的, 應該是差別在七參數的旋轉部份. 159樓講到的, 我發現我用的不是Bursa-Wolf. 或許改寫那一段, 算法就一致了.

Wei_1144 wrote:
如果是的話,可能要先修正算法,再求改善參數


修正算法. 唔. 七參數平移的部份沒什麼好修的, 旋轉的部份似乎有好幾種做法. 或者乾脆不用七參數?

... 我一直沒有搞清楚, 這個86年內政部七參數轉換是有現成的轉換程式嗎? 還是是以內政部的參數套用哪裡的程式?


這個要麻煩 hwsu315 解說一下
以上是以差異距離來做分析,再把距離的差異拆解成X、Y軸來看

我再貼

米轉換 的差異中
X軸的成分 與 分佈:


X軸的 成分 與 分佈,右旋90度來看:


Y軸的 成分 與 分佈:


Y軸的 成分 與 分佈,右旋90度來看:



米轉換內政部參數 的差異中
X軸的成分 與 分佈:


X軸的 成分 與 分佈,右旋90度來看:


Y軸的 成分 與 分佈:


Y軸的 成分 與 分佈,右旋90度來看:



86年內政部七參數轉換 的差異中
X軸的成分 與 分佈:


X軸的 成分 與 分佈,右旋90度來看:


Y軸的 成分 與 分佈:


Y軸的 成分 與 分佈,右旋90度來看:

再來看 差異的差異

((距離差異5 米轉換內政部參數) - (距離差異4 米轉換)) 之中
X軸的成分 與 分佈:


X軸的 成分 與 分佈,右旋90度來看:


Y軸的 成分 與 分佈:


Y軸的 成分 與 分佈,右旋90度來看:



((距離差異3 86年內政部七參數轉換) - (距離差異5 米轉換內政部參數)) 之中
X軸的成分 與 分佈:


X軸的 成分 與 分佈,右旋90度來看:


Y軸的 成分 與 分佈:


Y軸的 成分 與 分佈,右旋90度來看:





當然,這些都是假設7777與內政部的點位是一樣的

這個球到底該怎麼轉呢?
minstrelsy wrote:
我一直沒有搞清楚, 這個86年內政部七參數轉換是有現成的轉換程式嗎? 還是是以內政部的參數套用哪裡的程式?

真抱歉沒寫清楚
這裡的 "86年內政部七參數轉換" 是用我自己的excel來轉換的
目前好像沒有內政部七參數現成有轉換程式
我的excel裡有 TM2 <-> DMS, DMS <-> XYZ ,還有成大的轉換,手持機的Abridged Molodensky轉換 和 七參數Helmert轉換
所參考的公式資料連結:
TM2 <-> DMS 這個其實稍改一下就可換成TM3或TM6(UTM) 而改一下中央經線就可用在澎湖
DMS <-> XYZ
Abridged Molodensky和Helmert來源同上
在這期間,我也看到了不同出處的公式資料,大同小異,
本討論串最開始的某樓我附了我的 TM2 <-> DMS excel,和成大的結果是幾乎相同,但後來Excel愈加愈多,所以變得很亂,我也就沒再整理好分享上來,而且我還心想著應該寫個VB才不會excel中那麼多sheet這麼亂。

hwsu315 wrote:
本討論串最開始的某樓我附了我的 TM2 <-> DMS excel,和成大的結果是幾乎相同,但後來Excel愈加愈多,所以變得很亂,我也就沒再整理好分享上來,而且我還心想著應該寫個VB才不會excel中那麼多sheet這麼亂。


稍整理了一下我的excel
第一個附加檔案是Helmert 七參數轉換 由TWD97到TWD67
第二個附加檔案是Abridged Molodensky 由TWD67到TWD97
若要反向轉換, 交換Datum參數及 轉換參數變號 應該即可 誤差不大
裡面已有前幾樓的582筆資料 若要轉換更多筆 請copy & paste每個sheet的row來延伸
有興趣的人看了我的excel請別見笑 若有強者能改善 希望能再分享出來 這樣就最感激了
附加壓縮檔: 200803/mobile01-f3278a1011c6a53f13f23e1d9735e97f.zip
附加壓縮檔: 200803/mobile01-e8decb874364127e76ed4fe59a54f185.zip
三位高手是不是都在以內政部衛星控制站及7777三角點來推算轉換參數呢? 昨天自己終於找到內政部MTWDCON程式的指令格式了, 將溪子口山, 及虎子山的T97 TM2輸入來轉換, 發現與7777的資料所差無幾, 例如:

7777是 2765266.653, 304554.669 高度 142.649
MTWDCON轉出來是 2765266.6726, 304554.7205 高度是142.6281

如果直接將以MTWDCON程式為轉換依據是否足夠精確了呢? 上述誤差幾乎都在5cm上下. 如同林務局的Window轉換程式一樣, 便是直接將內政部程式產出的對照點作為轉換依據, 不知高手們看法如何?

同理, 將上述這兩筆資料全部以內政部程式轉換成"T97度"做成gpx檔, 這樣已經是接近毫無偏差的輸入Garmin GPS中了, 剩下的就只有驗證我用的三參數到底是否是最小誤差的參數組合了? 可惜, Garmin只能以整數輸入三參數! 當然, 921地震造成的三角點偏離原址也會造成資料有誤.

  • 41
內文搜尋
X
評分
評分
複製連結
請輸入您要前往的頁數(1 ~ 41)
Mobile01提醒您
您目前瀏覽的是行動版網頁
是否切換到電腦版網頁呢?