网形稳健度不同计算法之统计等值性检定.pptVIP

2005/10/7 網形穩健度不同計算法之統計等值性檢定 學生：林雅婷 指導教授：許榮欣 老師 流程 前言 變形指標 變形向量及穩健度計算 等值性檢定 實驗成果及分析 結論 前言(1/2) 測量所得之觀測量中難免有粗差存在而未被偵測出來，傳統上以Baarda的可靠度理論來分析，未被偵測出來的粗差對測量網所造成之影響(Baarda，1968)。 可靠度理論未探討沒被偵測出的粗差如何影響個別網點位置的估計，因此必須藉由穩健度分析理論來進一步分析。 Vaní?ek等(1991，2001)、陶本藻(1992)、許榮欣(2004，2005)皆有針對穩健度理論做過相關之研究。 前言(2/2) 本文將分別介紹Vaní?ek (1991，2001)與陶本藻(1992)所提之穩健度理論。 因為採用的方法不同，分析出之各個網點的三種變形指標值亦不相同，文中以平面大地網為研究資料，利用Vaní?ek與陶本藻的方法進行穩健度分析，則每個網點都會有兩組不同的變形指標值，最後透過統計方法中的等值性檢定，檢驗這兩種方法所計算出的變形指標，對同一個測量網而言，在統計上是否具有等值性。 變形指標(1/4) 觀測量中所含之粗差將導致測量網內點位產生移位，透過點位之移位量梯度可用以衡量點位之變形。 之移位量為：(許榮欣等，2005) 之變形矩陣為(Vaní?ek et al. 1991，2001)： 變形指標(2/4) 變形指標(Vaní?ek et al. 1991，2001)： 1.平均應變(Mean strain) 代表尺度變形(Deformation in scale)。 變形指標(3/4) 2.總剪應變(Total shear) 總剪應變代表形狀變形(Deformation in shape)。 變形指標(4/4) 3.局部微旋轉( Local differential rotation) 研究點對Z軸之微旋轉(Differential rotation)為 上式可再分解成區塊旋轉(block rotation) 與局部微旋轉(local differential rotation) 。 吾人常藉由局部微旋轉( Local differential rotation)來描述各個網點之局部扭轉(Local twisting)狀況。 變形向量及穩健度計算 變形向量之基本公式 Vaní?ek計算法 陶氏計算法 Vaní?ek計算法與陶氏計算法之比較 變形向量之基本公式(1/3) 點位之變形向量計算，除了考慮研究點 外，其周圍與之有連結或距離某半徑範圍內之t個 點，都必須列入一起考慮。 令 變形向量之基本公式(2/3) 變形向量可表示為(Hsu and Li，2004) 將移位向量 擴及全網，透過引入粗差向量 的概念，將變形矩陣改寫為：(Hsu and Li，2004) 變形向量之基本公式(3/3) 變形向量之基本公式 ，代表因第k個觀測量中未被偵測出之粗差所造成之移位向量。 Vaní?ek計算法(1/1) 每一個網點皆存有一個變形向量 ，變形向量的組成受到所有觀測量之影響。 每一點位之變形向量皆不相同，計算時皆須重複執行n次的迴圈(Do-Loop)來求解(許榮欣等，2005)。 透過各點之變形向量可決定各點之三種變形指標，並選取其絕對值最大者 為該研究點之穩健度指標。 對於每個點位而言，每個觀測量都可產生一個 ，則共有n個變形矩陣，而 可以產生3個變形指標，所以每一點位上可得3n個變形指標，即 各有n個。 當絕對值愈小則表示點位愈穩健愈不易受粗差影響。 陶氏計算法(1/1) 各網點之變形向量 由最大移位向量 組成 Vaní?ek計算法與陶氏計算法之比較 (1/1) 都是計算變形指標。 Vaní?ek計算法使用各點之變形向量 來計算每個點之三種變形指標，變形指標受到所有觀測量之影響。 Tao’s計算法僅使用最大的變形向量 來計算變形指標，所有的變形指標均來自 ，又稱為最大移位向量法(Hsu，2005)。 計算效率上以陶氏法較為簡便，然Vaní?ek計算法考慮所有觀測量中未被偵測出之粗差對於點位估計之影響，而陶氏法僅認為點位估計只受到最大移位向量中之2m個元素影響，因此就理論而言應以Vaní?ek計算法較為嚴謹。 等值性檢定(1/2) 等值性檢