台灣一等一級水準網之重新平差計算.pptVIP

台灣一等一級水準網之重新平差計算 指導教授：許榮欣 學生：鄞守毅 報告流程 前言 資料介紹 IAUE 定相關法 Theil model 實驗 結論 資料介紹 本研究所計算之資料，是由內政部委託國立成功大學衛星中心所進行的『台灣一等一級水準網測量督導查核工作』所用之資料。 根據「一等一級水準網測量督導查核工作總報告書」，此水準網共有1033筆觀測高差、1020個水準點，以基隆市民族英雄公墓內之K999水準點(5.66883m)為基準點。 資料介紹 本研究將水準網分為26條測線，17個待定高程，並以基隆市民族公墓內之K999為固定點。 簡圖與整理表如下頁所示： 路線簡圖 資料介紹 上頁中各測線先驗方差m是以下式計算得的 上式中 是各測段閉合差， 是各測段長，i=1~n是測段編號。 IAUE 迭代近乎無偏估計法(Iterated Almost Unbiased Estimation)，即IAUE，是一種估計最佳方差的方法，此法將根據網形幾何結構與觀測值品質估計各測線的方差及測線間的協方差，進而得到最佳的權矩陣。 。 IAUE計算流程(1) 首先應該要列出觀測方程式，並計算出各測線的先驗方差，然後進行分組，將所有觀測量依某種標準分為k組，並將同一組的先驗方差統一為該組方差的平均數。 IAUE計算流程(2) 。 計算流程(3) ，其中各個對角線元素就是各觀測量的多餘觀測數 計算流程(4） 計算各組的scale factor (f數) ，其中 是各組成員r的總和。 重複2～6步驟直到各組scale factor和平差後之後驗單位權方差都收斂到1為止。 IAUE 依照上面步驟進行到最後將會得到各觀測量的最佳方差估計，待方差估計完畢後，尤其所組成的權矩陣會是最佳的權矩陣，直接拿來做平差計算的結果也將是各待測點位的近乎無偏估計。 定相關法 定相關法也是期望為平差方法估計出最佳的權模式，讓平差結果更理想。 不同於IAUE的是它是直接求最佳權，而不是先求各測線之最佳方差，再來組成權。 相關係數R可以根據lag-one autocorrelation，也可以根據任何lag的autocorrelation來估計，而平差結果也會隨相關係數的選擇而異。 計算流程(1) 利用各測線間測段之閉合差，計算各測線的R（相關係數） 其中 代表各測段閉合差，i代表測線編號， 代表該測線的測段數，k代表lag數。 . 計算流程(2) 代入平差計算，得各點高程 定相關流程圖 Theil model 測量計畫中，有時會遇到加密的問題，此類問題，需先將某些點加以固定，以進行加密動作，而那些固定點資料就稱為約制條件。 加密問題可分兩種方法處理，一為不將已知點資料給權之平差，即不含虛擬觀測之平差；一為將已知點資料加權之平差，即加入虛擬觀測之平差。 Theil model 所謂虛擬觀測就是將已知點資料視為一已觀測之資料（實際上我們並沒有去觀測)，一起加入平差計算，藉此可以減少對網形的約制，避免網形為迎合約制條件而變形，或讓秩虧的設計矩陣得以滿秩。 Theil model 已知點資料並不是絕對精確的，凡是觀測或多或少必有誤差，只是我們有沒有將它反映出來爾，若為虛擬觀測加權，就是將它視為有誤差之量，再進行平差計算，則稱之為加權的約制/加密平差；若無給予加權，即視為無誤差之已知量，則稱為不加權的約制/加密平差。本研究即以加入虛擬觀測之Theil Model做加密計算。 計算流程(1) 將含有誤差之已知約制參數與未知的待求加密未知數分開，分別為 的 和 的 ，其中 ， 是已知的約制參數數目， 是待求的加密點數目 計算流程(2) Theil Model觀測方程式可寫成 : 計算流程(3) 將(7)代入(6)可得新虛擬觀測方程式： 其中 ， Theil model流程圖 IAUE實驗 IAUE實驗分成依： 先驗方差 地形 距離固定點距離 這三種標準來分組 依先驗方差分組 依地形分組 依距離分組