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

hwsu315

hwsu315
個人積分：5分
文章編號：5447097

5分

191樓

2008-03-25 20:11

Wei_1144 wrote:
Excel 表是現成...(恕刪)

Well done! 結果很不錯呢:)
您可以試試以最新算出得到的參數再作一次七個參數的最佳計算，看看是否能讓平均誤差更小。
另外我認為參考點應該多取才正確，少取的話可能只有這些少數的點十分吻合，但其他點反而誤差大了，您可以實驗看看。

millerliu wrote:

成大水工所網頁程式轉換誤差:
TWD67--TWD97(WGS84)誤差是 ( 0.352m, 0.3936m) = 0.3247m

我一直覺得成大的公式在本島上真的是"便宜又大碗"^^

中研院有個工具程式連結：GIS 工具程式
其中第31項是個CGI的轉換程式，對中研院轉換法有興趣也對CGI binary有研究的人不妨可參考看看

而第26項，用在商業GIS軟體上的Project定義檔，我發現TWD67它上面用的1/f是GRS67: 298.247167427
並不是298.25。
不曉得中研院網站上用的是哪一個

Wei_1144

Wei_1144
個人積分：540分
文章編號：5448315

540分

192樓

2008-03-25 21:51

hwsu315 wrote:
您可以試試以最新算出得到的參數再作一次七個參數的最佳計算，看看是否能讓平均誤差更小...(恕刪)

好主意
跑了10個循環，結果是:

dx 750.6095019
dy 358.2376998
dz 182.8246414
rx 6.74221E-07
ry 3.00371E-07
rz 2.0955E-09
M 0.999018561

差異_距離_平均 = 2.15

millerliu

millerliu
個人積分：339分
文章編號：5448806

339分

樓主

2008-03-25 22:22

Wei_1144 wrote:
沒問題啦...(恕刪)

等你的gpx檔囉, 順便試試七參數中的平移參數瞧瞧是否適用Garmin三參數? 可惜, Garmin只接收整數, 還是Magellan慷慨, 可以輸入小數點一位.

基石資料庫csv檔: 200803/mobile01-a9272e06937889f75ea918a267352dca.zip

millerliu

millerliu
個人積分：339分
文章編號：5448968

339分

樓主

2008-03-25 22:31

hwsu315 wrote:
我一直覺得成大的公式在本島上真的是"便宜又大碗"^^...(恕刪)

可惜不能進行批次作業, 外島座標181筆我還是一筆一筆輸入, then複製結果出來, 但是它可以計算以121為中央經線澎湖地區座標, 可惜無法用MTWDCON驗證外島誤差到底是多少? anyway, 對資料少的user來說, 成大水工所網頁還是一級棒的!

不過我倒是非常失望, 中研院提供軟體竟然這麼遜, 誤差到了2公尺!

Wei_1144

Wei_1144
個人積分：540分
文章編號：5449047

540分

195樓

2008-03-25 22:36

Wei_1144 wrote:
好主意跑了10個循環...(恕刪)

跑了100個循環的結果:

dx 749.0779102
dy 371.1820498
dz 189.5234563
rx 6.74221E-07
ry 3.00371E-07
rz 2.0955E-09
M 0.992637506

差異_距離_平均 = 1.94

好像越來越...不對勁

不管，有時間再跑它個1000次

Wei_1144

Wei_1144
個人積分：540分
文章編號：5452491

540分

196樓

2008-03-26 6:52

跑了1000個循環的結果:

dx 741.6152757
dy 434.0795152
dz 223.0937678
rx 6.74221E-07
ry 3.00371E-07
rz 2.09615E-09
M 0.962913744

差異_距離_平均 = 0.99

Wei_1144

Wei_1144
個人積分：540分
文章編號：5452520

540分

197樓

2008-03-26 7:09

不好意思，這些都是錯的，因為發現Excel少設定一個地方，重跑一次再修正結果

Wei_1144

Wei_1144
個人積分：540分
文章編號：5453789

540分

198樓

2008-03-26 9:57

millerliu wrote:
等你的gpx檔囉, ...(恕刪)

轉好了，沒有標高資料的點位，暫以0m帶入

附加壓縮檔: 200803/mobile01-6fe0c9b3fe865602eb32ae98c383cdcd.zip

Wei_1144

Wei_1144
個人積分：540分
文章編號：5454953

540分

199樓

2008-03-26 11:22

Wei_1144 wrote:
不好意思，這些都是錯的，因為發現Excel少設定一個地方，重跑一次再修正結果...(恕刪)

原來的設定沒有考慮高程的差異，所以把高程也納入考量

跑了100個循環的結果:

差異_距離_XYZ_平均 = 3.93

100個循環過程中 dx、dy、dz 與差異_距離_XYZ_平均的變化，發現其實只要算5個循環就夠了

millerliu

millerliu
個人積分：339分
文章編號：5455124

339分

樓主

2008-03-26 11:35

內政部另一個程式的用途剛才一個朋友提供一些指示, 供大家參考: (如下文)

Geoidint 可能是內插程式, 此程式可將大地起伏網格值內插任一位置（輸入十進
位之經緯度值），而參考之大地起伏網格檔（twgeoido1.grd3）亦由黃金維教授所
提供，為 30〞×30〞之解析度（內政部，2001），此網格與實測大地起伏比較，在
西部平原可達公分及之精度（Hwang, 1997）。將橢球高減去內插計算所得之大地
起伏值，即可得正高。

原文在
空載雷射掃瞄與航空測量數據高程精度比較之研究
http://nhmrc.nctu.edu.tw/People/tyshih/Publications/SE2003/51.pdf 有興趣 DEM資料的，可參考看看。