個人積分：545分
文章編號：5164067
個人積分：545分
文章編號：5167676
Wei_1144 wrote:
關於資料的檢驗與勘誤，我有個想法:
也許把不同來源的基點資料，分別轉成KML，在Google Earth裡面同時觀看，會比較簡單明暸一些
不同版本的基點 用不同顏色的標記來表示，但是都以統一的標點編號來標示，這樣就可以知道各版本基點的位置差異
用Google Earth來觀看有一個好處，可以切換 顯示或不顯示，這樣我們可以一次只看同一個基點的各種版本資料，逐個基點來比較，不需要全部的基點都顯示，那樣會 太"花"
如果一次只顯示特定基點的不同版本資料時，發現不同顏色(不同版本)的標點差異太大，我們就知道其中有個版本的資料是有問題的；反之，如果都很接近，資料大致上就是OK的；要整合或實地勘測時，應該以版本差異大的基點優先，效率與效益會比較好一點...
合併7777、小花、內政部 資料，簡單做一下差異比較:
不同版本資料中，基點名稱相同的基點會放在一起
. 紅色三角標: 7777
. 黃色三角標: 小花
. 綠色三角標: 內政部 資料
每個 基點名稱 資料夾 右方括號內的數字，代表同名基點中，(X_min, Y_min) 與 (X_max, Y_max) 的距離(公尺)
X_min = Min(X1, X2, X3...)，Y_min = Min(Y1, Y2, Y3...)，同法求最大值
KML 檔案: 基點_不同版本比較.kml
小花的資料是由.gdb轉出，經緯度精確到"度"小數點下七位
因為是要做版本差異比較，所以:
. 不含 7777資料中 無 點名(山名) 者
. 不含 內政部資料中 無 精確經緯度 者
. 所列出者 一定有不同版本的資料；無不同版本資料者 不列出；只在單一版本(例如 7777)中重複名稱者 不列出
個人積分：339分
文章編號：5228985
個人積分：545分
文章編號：5229331
個人積分：5分
文章編號：5232051
millerliu wrote:
http://www.geodetic.gov.hk/smo/gsi/programs/tc/geo.htm看到香港測繪處的網站, 將測點的詳細資料全部數字公布, http://www.geodetic.gov.hk/smo/gsi/data/satellite/HKFN.pdf看了真讓人羨慕, 大概香港人不認為這是國家機密吧!
謝謝miller兄提供的資訊，我也順便看了一下香港的測量座標系統
官網上有簡介說明，還有轉換公式：
http://www.geodetic.gov.hk/smo/gsi/data/pdf/explanatorynotes_c.pdf
還有使用手持Gpsr時，如Garmin，自訂Datum和grid的參數：
http://www.geodetic.gov.hk/smo/gsi/programs/tc/grid.htm
香港的官網還提供DGPS的功能
發現香港本來使用UTM，但小小的香港，竟然剛好要分兩個UTM的分帶 (剛好在東經114度，兩個分帶中間)，可能會有點不方便
後來1978年重測時，創了香港自己的Grid，不過參數複雜了點，False N, False E, 中央經線都不是整數，不過要是轉換程式設定正確，其實也是沒有問題的。測量基點的資料，是一個點一個pdf檔來供查詢，假如想自己Key-in到excel，工程也不小。
個人積分：545分
文章編號：5233086
個人積分：339分
文章編號：5233170
Wei_1144 wrote:
發現越遠離虎子山，偏差量會越多，有沒有辦法修正算法或參數? 高手給個意見
...(恕刪)
http://www.afasi.gov.tw/download/download02.htm這個轉換程式的說明已經提到具體作法了, 當然與我們用單一參數來轉換全島座標是不同的方式, 但是保證精準!!
個人積分：5分
文章編號：5237449
Wei_1144 wrote:
我用米的不落格中的程式碼來轉換全國三角點(TWD67TM2→TWD97DMS)、內政部(TWD97TM2→TWD97DMS)的資料，但是發現越遠離虎子山，偏差量會越多，有沒有辦法修正算法或參數?
我看了一下米的不落果裡程式碼，它用的dX,dY,dZ, rx,ry,rz,s 七參數和內政部86年控制點平差工作報告裡的參數有所不同，可能有些偏差是來自此。
前者為：
dX=-752.32
dY=-361.32
dZ=-180.51
rx=-0.00000117
ry=0.00000184
rz=0.00000098
s=0.00002329
而後者內政部86年：
dX=750.739
dY=359.515
dZ=180.510
rx=0.00003863
ry=0.00001721
rz=0.00000197
s=0.99998180
(以上是原始在網路上找到的資料，請注意前者的rx,ry,rz是徑度，內政部的rx,ry,rz是角秒，單位不同，前者的s是縮放率的差，後者的s是縮放率)
基本上我搞不清楚米的程式碼的七參數是怎麼來的，程式的參數甚至不在米自己整理的轉換參數列表中：
http://blog.minstrel.idv.tw/2004/06/taiwan-datum-parameter.html
程式碼中的註也提到參數的不確定。
另外您提到修正算法和參數：
以我所知，因為座標系統可能會因地形複雜而有變形，所以多項式多參數回歸的轉換法對此有變形的座標系是比七參數要更符合轉換的點，但多項式就是需要很多實際真實的點來求得參數，而且反轉換回來沒有小誤差的對稱性，要再求得另一組參數，這麼多參數，不適合在以前還沒電腦的時候來計算和使用，而且也不適合一般人用手持機時來自訂參數，所以在無高精度要求時還是用比較基本的轉換方法。
現在我們已有TWD97的官方精確座標，假設7777資料庫裡的也是官方精確座標(TWD67)，理論上我們自己就可以算出多項式的轉換公式的參數，以供自己轉換使用。
個人積分：5分
文章編號：5237929
內政部版本的參數, 記得之前有試過, 不過好像不盡理想, 所以garmin才又自己弄了一套出來.
而且TWD67的七參數部份, 就我所知, 每家地圖廠商都有自己的參數. 這部份很頭痛...
喔, 對了, 補充一點, 我自己所用的TWD67七參數, 是以921之前的數據比對出來的. 921之後, 虎子山自己也偏掉了不少, 當時為了要不要重做一次TWD67的偏移修正爭論了很久. 在檢驗參數的同時, 可能要考慮三角點的資料是什麼時候取得的.
個人積分：545分
文章編號：5239999
