• 沒有找到結果。

5-1 結論

(1)個別投影坐標系統之變形量具有強烈之地區性,因此在小區域內個別 系統自身仍維持著相當程度之相對精度,本實驗透過一個「地段」為 範圍進行「小規模」之轉換工作,由實驗結果發現可達到相當高之精 度。

(2)共同點之點數及品質是整個坐標轉換精度好壞的重要因素之ㄧ,本實 驗區外圍同時具有 TWD97 與 TWD67 之共同點清查完後僅有 4 點,經計 算後其共同點之改正數過大,且自由度過少,因此本實驗另採用均勻 分佈於實驗區四個象限之四等控制點,該 10 個點位雖僅有 TWD67 坐 標,吾人藉當年辦理重測作業時該 10 點與上述 4 點基本控制點 GPS 聯測之原始觀測資料,強制附合至該 4 點基本控制點上,重新解算出 10 點控制點之 TWD97 坐標,以這 10 個四等控制點為共同點,並以實 驗區內 125 個圖根點為檢核點,實驗結果發現以最小二乘配置法結合 四、六參數轉換中均僅有 3 個點位誤差大於 2 公分,122 個點位誤差 小於 2 公分,佔所有點數 97.6%;可見以最小二乘配置法結合四參數 或六參數兩種模式進行坐標轉換,於 125 個圖根點精度上沒有明顯的 差異,且較單純四、六參數轉換模式好。

(3)以最小二乘配置法結合六參數模式進行坐標轉換獲得以下成果:宗地 之界址點最大之誤差為 1.3 公分,平均點位誤差為 0.5 公分;每段邊 長平均增加 0.001 公尺,增減最大值為 0.027 公尺,最小值為-0.001 公尺,99.66%邊長增減在 1 公分以內;面積增加最多之比率為 0.228

%,減少最多的為 0.685%,平均每筆面積增加 0.007%,宗地面積增

加最大者為 1.23 平方公尺,減少最多的為 0.15 平方公尺,99.6%宗 地面積增減均在 0.1%內;由以上實驗成果可知坐標轉換之成果均可符 合地籍測量實施規則之要求,本坐標轉換可視為無誤差,並可直接以 轉換後之 TWD97 坐標進行複丈作業。

5-2 建議

(1)地籍圖是民眾土地產權界址重要依據,如需進行坐標系統轉換建議以 小區域為轉換單元,並配合最小二乘配置法來配賦控制點的殘餘誤 差,因其控制點坐標於轉換過程中不變,可兼顧法律面與實務面之執 行。另坐標系統經轉換後其邊長、形狀及面積檢查甚為重要,建議應 訂定嚴格之標準作業程序來檢核轉換後之成果,以確保民眾之產權。

(2)一般小區域坐標轉換,最麻煩的地方就是共同點數量難以掌握,本實 驗主要採用當年重測作業所設置之四等控制點為共同點,這些點位可 能因年久失修損毀或遺失,但由於當年觀測之原始觀測資料保持完 整,藉這些原始觀測資料重新化算,才使得坐標轉換作業可順利完成,

且轉換後之精度證實可視為無誤差,因此其他類似本實驗區之地段建 議可以本實驗方法進行坐標轉換,減少外業時間。

(3)坐標系統整合後,各種不同圖籍資料可藉地理資訊系統(GIS)予以整 合,除可降低因系統不同所引起之測量問題外,並可提供決策者下達 政策前輔助工具,因此如何將不同圖籍資料結合以方便相關圖籍管理 及廣泛運用於地理資訊系統 GIS 中,為目前重要議題之一。

參考文獻

1.內政部,2003,地籍測量實施規則。

2.于宗壽、魯林成,1982,測量平差基礎。

3.王蜀嘉,1984,最小二乘排位法在含有隨機特性之坐標轉換上的應用,

台南。

4.尤瑞哲,1998,三度空間 GPS 坐標和台灣地區二度 TM 坐標轉換可行性分 析之研究,NSC87-2211-E006-36。

5.黃華尉,2001,TWD97 與 TWD67 二度 TM 坐標轉換之研究,台南。

6.許皓寧、蔡榮得,2004,台北市地籍資料 TWD67 與 TWD97 坐標轉換之比 較研究,地籍測量,第 23 卷、第 1 期,中華民國地籍測量學會會刊,

P39~63。

7.戴翰國、余致義、曾清涼,2002,利用六參數平面轉換與最小二乘配置 進行小區域 TWD67 與 TWD97 之地籍資料坐標轉換-以台北市大安通化段為 例,第五屆 GPS 衛星科技研討會,P66~71。

8.邱仲銘、洪慧齡,1999,土地測量成果坐標整合之研究",第十八屆測量 學術及應用研討會論文集。

9.李祺滄、高書屏,2005,整合地籍圖、航測地形圖與都市計畫圖套疊後 精度分析之研究,國立中興大學碩士學位論文。

10.曾清涼、儲慶美,1999,GPS 衛星測量原理與應用。

11.楊名,2003,TWD67 與 TWD97 坐標系統轉換。

12.尤瑞哲,2006,最小二乘配置法坐標轉換電腦程式 COLLOCATION 使用手 冊。

13.Erdas Inc.,1999,Erdas Field Guide Version 8.5.,USA。

14. Vanieek, P. and Steeves, R. R.,1996. “ Tramsformation of Coordinates between Two Horizontal Geodetic Datums", Journal of Geodesy,Vol 70,pp 740-745.

15.Vincenty,T.,1980. “ Height-Controlled Three-Dimensional Adjuestment of Horizontal Network " , Bull. Geod. , Vol.54 , pp.37-43

附 表

四參數 四參數+

四參數 四參數+

四參數 四參數+

附表 B 界址點協助指界內容

TWD67 TWD97 較差

附表 C 邊界樁距離與方位角度差明細表

起始邊界樁 結束邊界樁 距離差(公分) 方位角度差(秒)

S101 S102 0.2 -1.9

S102 S103 0.1 0.4

S103 S104 0.2 0.7

S104 S105 -0.3 -0.1

S105 S106 0.3 -6.7

S106 S107 0.7 -6.3

S122 S123 0.3 -0.4

S123 S124 0.3 0.8

S124 S125 0.6 1.2

S125 S126 1.3 1.2

S126 S127 0.2 0.4

S127 S128 0.3 1.5

S128 S129 0.1 0.1

S129 S130 0.3 1.0

S130 S131 0.1 -0.4

S131 S132 0.3 1.0

S132 S133 0.3 -0.5

S133 S134 0.2 0.2

S134 S135 0.4 -0.1

S135 S136 0.1 -0.2

S136 S137 0.1 -0.8

S137 S138 0.2 -0.7

S138 S139 0.1 0.0

S139 S140 0.2 0.1

S140 S141 0.1 -0.5

S141 S142 0.2 0.4

S142 S143 0.3 0.0

S143 S144 0.1 2.4

S144 S145 -0.2 0.9

S145 S146 -0.1 0.3

附表 D 控制測量成果基線成果比較表

網型尺度比 = 1.000002179175

=============================================================

後視 前視 觀測距離 平差後距離 較差 精度 中誤差 (m) (m) (m) (1/ppm) (m) 3 A13 A12 1029.618 1029.619 0.001 1272705 0.007 6 A13 F031 2672.499 2672.492 -0.007 406772 0.018 9 A13 A17 1073.398 1073.400 0.001 943233 0.009 12 A13 1096 4817.381 4817.380 -0.001 4807764 0.009 15 A13 5122 1803.802 1803.807 0.005 383462 0.009 18 A13 S117 868.689 868.684 -0.005 166001 0.008 21 S117 A12 1571.597 1571.590 -0.007 240415 0.007 24 S117 F031 3100.859 3100.854 -0.006 540218 0.019 27 S117 A17 1876.858 1876.855 -0.003 620038 0.010 30 S117 1096 5671.113 5671.108 -0.005 1083720 0.010 33 S117 5122 1900.861 1900.867 0.005 355103 0.010 36 A12 F031 1642.886 1642.888 0.002 671938 0.019 39 A12 A17 1593.057 1593.061 0.003 466353 0.010 42 A12 1096 4749.516 4749.516 -0.000 9000000 0.010 45 A12 5122 2813.197 2813.203 0.006 469572 0.008 48 A17 F031 3019.238 3019.247 0.009 344270 0.021 51 A17 1096 3813.874 3813.873 -0.001 3080673 0.012 54 A17 5122 1706.083 1706.087 0.004 421672 0.009 57 F031 1096 5080.094 5080.101 0.007 725107 0.021 60 F031 5122 4441.948 4441.955 0.007 661596 0.020 63 5122 1096 5015.115 5015.118 0.004 1432073 0.011 66 A13 TS11 1017.040 1017.037 -0.003 361036 0.013 69 A13 S012 1291.610 1291.613 0.003 426557 0.015 72 A13 A12 1029.619 1029.619 -0.000 2105559 0.014 75 A13 F031 2672.503 2672.492 -0.011 248743 0.015

78 A13 A17 1073.402 1073.400 -0.002 439558 0.023

180 S117 5122 1900.867 1900.867 -0.000 9000000 0.011

282 5122 1096 5015.135 5015.118 -0.017 298412 0.013

384 GA02 GA07 484.126 484.134 0.008 59454 0.017

486 S119 F718 4511.043 4511.045 0.001 3478061 0.011 489 S119 1071 3291.667 3291.677 0.009 354287 0.014 492 S012 F718 4780.564 4780.567 0.003 1785793 0.019 495 S012 1071 3196.766 3196.770 0.004 764229 0.021 498 1071 F718 3533.267 3533.262 -0.004 834497 0.011

=============================================================

*本控制網共有 166 條基線,其中有 0 條基線低於 1/20,000

附表 E 參數轉換成果

一、以四個共同點坐標轉換

(一)四參數坐標轉換

(1)輸入測區四個已知點 TWD97、TWD67 坐標

1096,2737408.677,322610.178,2737204.077,323441.302 D029,2740508.157,331755.179,2740303.652,332586.571 F031,2741114.792,326084.547,2740910.183,326915.793 F037,2738863.382,328848.101,2738658.828,329679.520 -99,

(2)輸出坐標轉換成果

區域中心座標:N0= 2739473.752 E0= 327324.50125 E(X) = A * (E - E0) + B * (N - N0) + C

N(Y) =- B * (E - E0) + A * (N - N0) + D A= 1.00002847713642

B=-3.78439410957165E-06 C= 328155.7965

D= 2739269.185

OFFSET E(X)=831.295 N(Y)=-204.567

旋轉 : - 0 度 00 分 00.7806 秒 (順正逆負) Scale : 1.00002847714358

轉換後區域中心座標 2739269.1850 328155.7965

點名 原座標(N) 原座標(E) => 轉換座標(N) 轉換座標(E) DeltaN DeltaE 點位誤差

1096 2737408.677 322610.178 2737204.033 323441.347 -0.044 0.045 0.063 D029 2740508.157 331755.179 2740303.636 332586.597 -0.016 0.026 0.030 F031 2741114.792 326084.547 2740910.267 326915.801 0.084 0.008 0.084 F037 2738863.382 328848.101 2738658.803 329679.442 -0.025 -0.078 0.082

(二)、六參數坐標轉換

(1)輸入測區四個已知點 TWD97、TWD67 坐標

1096,2737408.677,322610.178,2737204.077,323441.302 D029,2740508.157,331755.179,2740303.652,332586.571 F031,2741114.792,326084.547,2740910.183,326915.793 F037,2738863.382,328848.101,2738658.828,329679.520 -99,

(2)輸出坐標轉換成果

區域中心座標:N0= 2739473.752 E0= 327324.50125 E = A * (E - E0) + B * (N - N0) + E

N = C * (E - E0) + D * (N - N0) + F A= 1.0000356842439

B=-1.28311098208534E-05 C= 1.40688195112507E-05 D= 0.999987582032065 E= 328155.7965

F= 2739269.185

點名 原座標(N) 原座標(E) => 轉換座標(N) 轉換座標(E) DeltaN DeltaE 點位誤差

================================================================================================

1096 2737408.677 322610.178 2737204.069 323441.332 -0.008 0.030 0.031 D029 2740508.157 331755.179 2740303.639 332586.619 -0.013 0.048 0.050 F031 2741114.792 326084.547 2740910.187 326915.777 0.004 -0.016 0.017 F037 2738863.382 328848.101 2738658.844 329679.458 0.016 -0.062 0.064

一、以十個共同點坐標轉換

OFFSET E(X)=831.402 N(Y)=-204.618

旋轉 : - 0 度 00 分 03.8068 秒 (順正逆負) Scale : 1.00002890556789

轉換後區域中心座標 2737676.2999 329777.4735

點名 原座標(N) 原座標(E) => 轉換座標(N) 轉換座標(E) DeltaN DeltaE 點位誤差

(二)六參數坐標轉換

以四參數及最小二乘配置坐標轉換

(一)四參數坐標轉換

(1)輸入測區四個已知點 TWD97、TWD67 坐標 宜蘭網

GA01 2737282.168 328859.683 2737077.527 329691.086 GA02 2737421.264 329399.206 2737216.643 330230.623 GA03 2737516.898 329998.788 2737312.288 330830.219 GA04 2738375.472 329709.288 2738170.894 330540.705 GA05 2738603.385 328861.188 2738398.810 329692.578 GA06 2738110.469 328674.851 2737905.873 329506.242 GA07 2737861.922 329198.712 2737657.308 330030.123 TS012 2737877.277 328236.476 2737672.635 329067.866 TS11 2738318.784 328187.55 2738114.140 329018.954 TS119 2737441.54 328334.969 2737236.881 329166.339 -99

N039 2744578.720 324412.922 2744374.507 325244.357

(2)輸出坐標轉換成果

VARIANCE-COVARIANCE MATRIX OF TRANSFORMATION PARAMETERS =======================================================

.4546E-10 .0000E+00 .1270E-19 .2381E-20 .4546E-10 -.2381E-20 .1270E-19 .2458E-04 .0000E+00 .2458E-04

VARIANCE AND COVARIANCE =======================

2 .000 -1. .0001786

-1 : DETERMINED AUTOMATICALLY BY PROGRAM

COVARIANCE FUNCTION ===================

COVARIANCE FUNCTION X--X : .0002146455 - .000000021465 * D +- .000000 +- .000000000000 COVARIANCE FUNCTION Y--Y : .0001786284 - .000000017863 * D +- .000000 +- .000000000000

FINAL SETTINGS OF CONTROL PARAMETERS ======================================

REFERENCE POINTS ================

PERCENTAGE OF FILTERED NOISES AT YE-COMPONENT = .00 %

INTERPOLATED POINTS ===================

NAME O-XON O-YOE C-XNN C-YNE SX SY --- B01 2738540.913 329031.656 2738336.337 329863.051 .022 .003 B02 2738517.897 329179.116 2738313.321 330010.516 .019 .002 B03 2738532.476 329406.984 2738327.902 330238.390 .017 .003 B04 2738455.659 329551.955 2738251.083 330383.366 .014 .002 B05 2738116.424 329845.748 2737911.837 330677.171 .007 -.001 B06 2737893.881 329884.275 2737689.285 330715.701 .004 -.004 B07 2737630.544 329932.522 2737425.938 330763.951 .001 -.007 B08 2737380.605 329924.260 2737175.990 330755.691 -.001 -.009 B09 2737375.621 329790.639 2737171.004 330622.067 .000 -.009 B10 2737358.853 329623.869 2737154.233 330455.292 .000 -.008 B11 2737351.208 329540.447 2737146.586 330371.868 .001 -.008 B12 2737335.630 329273.270 2737131.002 330104.685 .000 -.007

  (其餘資料 略)

相關文件