111 年 國立成功大學測量及空間資訊學系碩士班《測量學》
第 1 題10 分
為檢驗某電子秤,以此電子秤量測 1kg 砝碼 6 次得觀測數據如下表:
觀測值(g) 995 991 992 994 990 992
(1). 請計算這組觀測值之觀測誤差、平均誤差及均方根誤差。
(2). 請計算這組觀測值之平均值、殘差及樣本標準偏差。
(3). 請進一步以精密度(precision)及精確度(accuracy)的概念,說明這組觀測資料的品質。
登入後即可作答並保存紀錄。
本題主要在考驗觀測數據的基本統計分析,包括觀測誤差、平均誤差、均方根誤差、平均值、殘差、樣本標準偏差,以及對數據品質的判斷。
核心概念:
- 觀測誤差 (Observed Error):指單一觀測值與真實值之間的差異。在實際測量中,真實值通常未知,我們常用平均值來近似真實值。
- 平均誤差 (Mean Error, ME):觀測誤差的平均值。理論上,平均誤差應趨近於零。
- 均方根誤差 (Root Mean Square Error, RMSE):也稱為標準差 (Standard Deviation),是衡量觀測值離散程度的指標。
- 平均值 (Mean):所有觀測值的總和除以觀測次數。
- 殘差 (Residual):觀測值與其平均值之間的差異。
- 樣本標準偏差 (Sample Standard Deviation):用來估計母體標準差的統計量,計算公式為 。
- 精密度 (Precision):指觀測值之間的一致性或分散程度。精密度越高,觀測值越接近彼此。
- 精確度 (Accuracy):指觀測值與真實值之間的接近程度。精確度越高,觀測值越接近真實值。
解題步驟:
首先,我們將觀測值整理如下:
g
g
g
g
g
g
觀測次數 。
(1). 計算觀測誤差、平均誤差及均方根誤差
由於真實值未知,我們通常以平均值來近似真實值。
計算平均值 ():
g
計算殘差 ():
觀測誤差 (Error):在實際情況下,我們無法直接得知真實值,因此無法計算精確的觀測誤差。我們通常用殘差來近似觀測誤差。
例如,第一個觀測值的近似觀測誤差為 g。
平均誤差 (Mean Error, ME):理論上,平均誤差是所有觀測誤差的平均值。在實際應用中,我們計算殘差的平均值,它應該趨近於零。
ME = g。
這表示在這一組數據中,觀測值的偏差是相互抵銷的。
均方根誤差 (RMSE):這是衡量數據分散程度的重要指標,通常就是樣本標準偏差。
首先計算殘差平方和:
第 2 題10 分
以捲尺觀測兩點間之水平距離共 6 次,得觀測數據如下表,若以觀測量品管之信心
水準為99%(即E99 = 2.5σ)之容許誤差檢視觀測量,則此組觀測量中是否有錯誤觀
測量?(10%)
觀測值(m) 45.120 45.118 45.121 45.126 45.222 45.129
登入後即可作答並保存紀錄。
本題考驗如何使用統計方法來檢核觀測數據中是否存在異常值 (outlier),並應用品管中的容許誤差概念。
核心概念:
- 信心水準 (Confidence Level):表示我們對觀測結果的信賴程度。
- 容許誤差 (Allowable Error):在給定的信心水準下,允許觀測值與真實值之間的 maximum 差異。
- 異常值 (Outlier):顯著偏離其他觀測值的數據點,可能由儀器故障、人為錯誤或特殊現象引起。
- 標準差 (Standard Deviation, σ):衡量數據分散程度的指標。
- E99:在 99% 信心水準下的容許誤差。題目中給定 E99 = 2.5σ。
解題步驟:
首先,我們需要計算這組觀測數據的平均值和標準偏差。
觀測值 ():45.120, 45.118, 45.121, 45.126, 45.222, 45.129 (m)
觀測次數 () = 6
1. 計算平均值 ():
m
2. 計算殘差 ():
第 3 題10 分
某直角三角形土地如下圖,觀測了C邊距離得 c=35.623±0.012m,並觀測 θ
角得 θ = 24°20′35″±30”,請計算a及b邊距離及其中誤差 σα及σb。(10%)
🖼️【此處有附圖,請對照原卷】
(圖示為一個直角三角形,其中斜邊為 c,夾角 θ 的對邊為 b,鄰邊為 a。)
登入後即可作答並保存紀錄。
本題考驗如何利用三角學關係以及誤差傳播定律來計算邊長的中誤差。
核心概念:
- 三角學關係:在直角三角形中,邊長與角度之間的關係。
- 誤差傳播定律 (Law of Propagation of Errors):當一個量是其他幾個量的函數時,這些量的誤差如何影響該量的誤差。
- 中誤差 (Mean Error, ME):一般情況下,我們用標準差 (σ) 來表示中誤差。
- 角度單位轉換:將度分秒轉換為十進位度或弧度,以便於計算。
解題步驟:
首先,我們需要將給定的角度轉換為十進位度或弧度。題目中的角度誤差是 ±30”,通常在誤差傳播中,角度誤差是以弧度為單位。
1. 角度單位轉換:
的誤差
將角度誤差轉換為弧度:
radians
將角度 轉換為十進位度:
轉換為弧度:
radians
2. 三角學關係:
從圖中可知,直角三角形的邊長與角度關係為:
3. 誤差傳播定律應用:
計算邊長 a 的中誤差 :
計算偏導數:
代入誤差傳播公式:
計算邊長 b 的中誤差 :
計算偏導數:
代入誤差傳播公式:
4. 代入數值計算:
已知:
m, m
(或 0.42489 radians)
radians
第 4 題10 分
設採用A方法觀測某筆土地面積,共觀測4次,得觀測數據如下表,請依據最小二
乘平差法:(10%)
A方法觀測值(公頃) 2.5531 2.5523 2.5549 2.5566
(1). 列出所有的觀測方程式
(2). 列出最小二乘之解算條件式(即殘差之平方合為最小)
(3). 解算此筆土地面積之最或是值。
登入後即可作答並保存紀錄。
本題考驗最小二乘平差法的基本應用,特別是針對單一未知數(此處為土地面積)的觀測情況。
核心概念:
- 最小二乘法 (Least Squares Method):一種統計方法,用於找到最佳擬合模型,使得觀測值與模型預測值之間的殘差平方和最小。
- 觀測方程式 (Observation Equation):描述觀測值與未知數之間的數學關係。
- 解算條件式 (Condition Equation):最小二乘法的核心,即最小化殘差平方和。
- 最或然值 (Most Probable Value, MPV):在最小二乘法中,指對未知數的最佳估計值。
解題步驟:
已知觀測了某筆土地面積 4 次,得到四個觀測值。我們假設土地面積為一個單一未知數 。
觀測值 ():2.5531, 2.5523, 2.5549, 2.5566 (公頃)
觀測次數 () = 4
(1). 列出所有的觀測方程式
觀測方程式描述觀測值與未知數之間的關係。在此情況下,我們假設每一次的觀測都試圖測量同一個未知數 。
若設 為土地面積的最或然值,則觀測方程式可以寫成:
其中 是第 次的觀測值, 是我們要求的最或然值,而 是第 次觀測的殘差。
所以,四個觀測方程式為:
(2). 列出最小二乘之解算條件式(即殘差之平方合為最小)
根據最小二乘法的原理,我們需要最小化所有殘差的平方和,即 。
從觀測方程式,我們可以得到殘差的表示式:
第 5 題10 分
續上題,若再以B方法觀測這筆土地面積,共觀測3次,得觀測數據如下表,設A
方法之先驗中誤差為±0.0015公頃,而B方法之先驗中誤差為±0.0006公頃:
B方法觀測值(公頃) 2.5561 2.5569 2.5563
(1). 請以加權最小二乘平差法列出所有的觀測方程式,並依據先驗中誤差釐定權值
(2). 請整合A及B方法之觀測值,以加權平差方式解算此筆土地面積之最或是值
(3). 請計算殘差、後驗單位權中誤差及最或是值中誤差
登入後即可作答並保存紀錄。
核心觀念
本題考查加權間接觀測平差(Weighted Observation Equation Adjustment)或多組不等精度觀測值之加權平均值解算。核心概念與公式包含:
- 權(Weight)與中誤差(Standard Error)之關係:
其中 為任選之先驗單位權中誤差, 為各觀測群之先驗中誤差。權與誤差平方(變異數)成反比。 - 觀測方程式(Observation Equations):
設單一未知數(面積最或是值)為 ,各項觀測值為 ,殘差為 ,則觀測方程式為:
- 加權最或是值(Weighted Mean):
- 殘差與後驗單位權中誤差(Posterior Reference Standard Error):
其中 為總觀測數, 為獨立未知數個數(本題 ), 為多餘觀測數(自由度)。 - 最或是值中誤差(Standard Error of the Weighted Mean):
前題(第 4 題)觀測資料回顧
本題題幹標明「續上題」,成大測量所 111 年第 4 題中,A 方法對該筆土地進行了 4 次獨立觀測,觀測值分別為:
結合第 5 題 B 方法的 3 次觀測值:
總觀測個數 。
解題方法與詳細計算
(1) 釐定權值與列出觀測方程式
-
依據先驗中誤差釐定權值:
已知 A 方法先驗中誤差 ,B 方法先驗中誤差 。
選取 (或任意便利常數,亦可取 之公倍數):
即 A 方法每測次之權 ,B 方法每測次之權 (兩者權之比值為 )。 -
列出觀測方程式:
設該筆土地面積之最或是值為 (公頃)。
為計算簡潔,取概略值 ,設改正數為 。各觀測值減去 後化為以 為單位的常數項 :- A 方法觀測方程式:
- B 方法觀測方程式:
- A 方法觀測方程式:
(2) 整合 A 及 B 方法觀測值解算最或是值
總權值 :
採用相對概略值 之化簡數值 :
計算加權和 :
解算改正數 :
故面積最或是值 為:
(3) 計算殘差、後驗單位權中誤差及最或是值中誤差
第 6 題20 分
已知 AB 兩樁位點之平面(E, N) 坐標為 A (168500.123 m, 2545003.361 m) 及 B
(168534.981 m, 2545123.334 m),P 點為未知點。若於 P 點架設全站儀觀測得水
平距離 LAP = 120.534 m,並觀測得水平角∠APB = 63°45′22": (20%)
(1). 請依題目的敘述,繪製點位及觀測量關係圖
(2). 請依據關係圖計算∠OBAP
(3). 請計算方位角QAB及ΦΑΡ
(4). 請計算 P 點坐標
(5). 設 AB 點坐標無誤差,而此儀器之先驗測距精度為±(5mm+10ppm),先驗測角
精度為±30”,請計算 P 點坐標中誤差
登入後即可作答並保存紀錄。
核心觀念
本題利用三角形 的幾何關係:
- 坐標已知,可求 長度及方位角。
- 已知 與 ,利用正弦定理求 。
- 由 方位角及 ,求 方位角。
- 以極坐標轉換計算 點坐標。
- 坐標中誤差則利用距離與方向角誤差傳播公式計算。
題目中的「」依幾何敘述解讀為 。
(1) 點位及觀測量關係圖
已知:
關係圖如下:
B
●
/ \
/ \
AB / \ BP
/ \
●---------●
A AP P
∠APB
其中:
- 為已知基線;
- 為全站儀觀測之水平距離;
- 為在 點觀測之水平角;
- 為待求點。
(2) 計算
先求
坐標差為:
因此:
在三角形 中, 對應角 , 對應角 。依正弦定理:
因此:
代入:
得到:
所以:
(3) 計算方位角 及
方位角
方位角以北方向為 ,順時針增加,因此:
故:
方位角
由幾何關係, 與 的夾角為 。
若依關係圖取 位於 的右側,則:
由於題目文字沒有明確指定 位於 哪一側,幾何上另有鏡射解:
換算為 :
(4) 計算 點坐標
方位角以北為基準時,坐標增量公式為:
採用關係圖所示的 :
因此:
第 7 題10 分
設某一點位之平面 EN坐標中誤差為σE = ±18mm 及σN = ±12 mm:(10%)
(1). 假設兩坐標誤差不相關(即σEN = 0),則其坐標值之協變方矩陣為何?請計算
誤差橢圓之長短軸及長軸傾角,並繪圖表示之(須標註橢圓外框、長短軸及長
軸傾角)。(誤差橢圓參考公式如下)
2σεν συ = ±√osin20 + 25encososine + oßcos20
參考公式:tan(20) = σ-σε σν = ±√2cos20 – 20ɛncososine + osin20
(2). 又若兩坐標誤差之協變方為σEN = -120mm²,則其坐標值之協變方矩陣為
何?其點位誤差相關性為何?請計算誤差橢圓之長短軸及長軸傾角,並繪圖表
示之(須標註橢圓外框、長短軸及長軸傾角)
登入後即可作答並保存紀錄。
本題考驗誤差橢圓的計算與應用,涉及到協方差矩陣、特徵值與特徵向量的計算,以及誤差傳播定律在多變量情況下的應用。
核心概念:
- 誤差橢圓 (Error Ellipse):在二維坐標系統中,表示一個點位置誤差的統計分佈。橢圓的形狀和方向反映了誤差的大小和相關性。
- 協方差矩陣 (Covariance Matrix):描述多個隨機變量之間協方差的矩陣。對於坐標誤差 (E, N),其協方差矩陣為 。
- 長短軸 (Semi-major and Semi-minor Axes):誤差橢圓的半徑,分別代表最大和最小的標準差。
- 長軸傾角 (Orientation Angle):誤差橢圓長軸與坐標軸的夾角。
- 特徵值與特徵向量 (Eigenvalues and Eigenvectors):協方差矩陣的特徵值對應於誤差橢圓長短軸的平方,特徵向量對應於長短軸的方向。
解題步驟:
已知 mm, mm。
注意單位:mm。
(1). 假設兩坐標誤差不相關(即σEN = 0)
協方差矩陣 (Covariance Matrix):
當誤差不相關時,協方差 。
誤差橢圓的長短軸及長軸傾角:
對於對角協方差矩陣,誤差橢圓的長短軸與坐標軸平行。
長軸的平方是較大的特徵值,短軸的平方是較小的特徵值。
在此情況下,特徵值直接就是對角線上的元素。
長軸的平方 (因為 )
短軸的平方
長軸 mm
短軸 mm
長軸傾角:
當 時,長軸與坐標軸平行。
長軸沿著 E 軸(如果 )或 N 軸(如果 )。
在此情況下,長軸沿著 E 軸,短軸沿著 N 軸。
所以,長軸傾角為 0°。
繪圖表示:
繪製一個坐標系,標示 E 和 N 軸。
以原點為中心,繪製一個橢圓。
長軸沿 E 軸,長度為 2a = 36 mm。
短軸沿 N 軸,長度為 2b = 24 mm。
橢圓外框應標示長軸 (36 mm) 和短軸 (24 mm),以及長軸傾角 (0°)。
【答案】
協方差矩陣:。
誤差橢圓長軸 mm,短軸 mm。
長軸傾角為 0°。
(2). 兩坐標誤差之協變方為σEN = -120mm²
協方差矩陣 (Covariance Matrix):
點位誤差相關性:
由於 是負值,表示 E 和 N 坐標誤差之間存在負相關性。
當 E 坐標的誤差偏向正值時,N 坐標的誤差傾向於偏向負值,反之亦然。
誤差橢圓的長短軸及長軸傾角:
長短軸的平方是協方差矩陣的特徵值。
特徵值 滿足 。
第 8 題10 分
何謂衛星靜態基線測量?若應用兩部接收儀進行兩個地面點位A及B點之靜態基線
測量,觀測計算得這兩個點位間之基線分量,請說明此基線分量之定義(請繪圖並
以文字描述之)。若其中A點之地心(地固)坐標為已知點,B點則為未知點,請說明
如何利用所得之基線分量計算B點坐標。(10%)
登入後即可作答並保存紀錄。
本題考驗對衛星靜態基線測量的理解,包括基線分量的定義、計算方法以及應用於未知點坐標確定。
核心概念:
- 衛星靜態基線測量 (Satellite Static Baseline Measurement):使用兩台或多台 GNSS 接收機在固定點位上長時間連續接收衛星信號,用於測量兩個點位之間的高精度相對位置(基線向量)。
- 基線向量 (Baseline Vector):連接兩個測量點的向量,其分量表示兩個點在三維空間中的相對位置差。
- 基線分量 (Baseline Components):基線向量在不同坐標系下的分量(例如,當地水平坐標系或大地坐標系)。
- 地心坐標系 (Geocentric Coordinate System):以地球質心為原點,X 軸指向本初子午線與赤道的交點,Z 軸指向北極,Y 軸垂直於 XZ 平面構成右手系。
- 地固坐標系 (Earth-Fixed Coordinate System):與地球剛體一起旋轉的坐標系,通常與地心坐標系相同。
- 大地坐標系 (Geodetic Coordinate System):以大地水準面為參考橢球體的坐標系,通常用經度 (λ)、緯度 (φ) 和大地高 (h) 表示。
- 相對定位 (Relative Positioning):利用兩個接收機之間的觀測數據來確定它們之間的相對位置。
解題步驟:
(1). 何謂衛星靜態基線測量?
衛星靜態基線測量是指,在測量過程中,兩台或多台 GNSS 接收機的天線被固定在地面上的已知或未知點位上,並進行長時間(通常為數小時或數天)連續的衛星信號觀測。這種測量方法主要用於獲取兩個點位之間的相對高精度位置信息,即基線向量。由於接收機是靜止的,可以積累大量的觀測數據,通過後處理軟件進行數據解算,消除或減弱隨機誤差和部分系統誤差,從而獲得非常精確的基線結果。這種方法是大地測量、工程測量和科學研究中測量遠距離高精度基線的標準方法。
(2). 請說明此基線分量之定義(請繪圖並以文字描述之)。
基線向量 表示從點 A 到點 B 的相對位置向量。
其中 和 分別是點 A 和點 B 的位置向量。
基線分量是指這個基線向量在某個特定坐標系下的投影。最常見的是在大地坐標系(經度 λ, 緯度 φ, 大地高 h)或地心直角坐標系 (X, Y, Z) 下表示。
-
在地心直角坐標系 (X, Y, Z) 下的基線分量:
設點 A 的地心坐標為 ,點 B 的地心坐標為 。
則基線向量在 X, Y, Z 方向的分量為:
基線向量 。 -
在大地坐標系 (λ, φ, h) 下的基線分量:
這種表示方式更直觀。基線向量可以分解為三個方向:
a. 徑向方向 (Radial direction):垂直於參考橢球面的法線方向。
b. 經向方向 (Meridional direction):沿著參考橢球面的子午線方向。
c. 緯向方向 (Latitudinal direction):沿著參考橢球面的卯酉圈方向。然而,更常用的是將基線向量分解為:
- 東西方向 (East, dE):與當地經線垂直,指向東方的方向。
- 南北方向 (North, dN):沿著當地經線方向,指向北方的方向。
- 天頂方向 (Up, dU):垂直於當地水平面,指向天頂的方向。
如果點 A 的大地坐標為 ,點 B 的大地坐標為 。
則基線分量 (近似值,對於短基線) 可以表示為:
(需要嚴謹的轉換公式)
第 9 題10 分
我國基本地形圖的圖幅分幅方式為何?此種圖幅經台灣TM2投影後之圖面圖幅框線
之形狀為何?若圖幅編號為9418、9418-11、9418-II-NE及9418-11-010,請說明這
些圖幅編號是屬於哪一種比例尺的基本地形圖圖幅編號?(10%)
登入後即可作答並保存紀錄。
本題考驗對我國基本地形圖圖幅分幅系統、投影方式以及圖幅編號規則的掌握。
核心概念:
- 基本地形圖 (Basic Topographic Map):國家統一編制,反映地表自然與人文要素的比例尺地圖。
- 圖幅分幅方式 (Map Sheet Division System):將大範圍區域劃分為若干圖幅的規則。
- 投影方式 (Map Projection):將地球表面曲面信息轉換到平面地圖上的數學方法。
- TM2 投影 (Transverse Mercator Projection, Zone 2):一種橫軸墨卡托投影,台灣地區常用的投影方式。
- 圖幅編號 (Map Sheet Numbering):用來識別和定位特定圖幅的代碼。
解題步驟:
(1). 我國基本地形圖的圖幅分幅方式為何?
我國基本地形圖主要採用 經緯度分幅 和 方里網分幅 相結合的方式。
- 經緯度分幅:以經緯度線為界,將全國劃分為不同大小的圖幅。
- 方里網分幅:基於經緯度分幅的基礎上,再細分為更小的圖幅,並採用方里網(例如 10 公里 x 10 公里)進行標示。
最常見的分幅系統是基於 1:25,000 的基本地形圖,其圖幅範圍約為 7.5 分鐘經度 x 5 分鐘緯度。
更大比例尺的地圖(如 1:5,000, 1:10,000)則基於 1:25,000 的圖幅進一步細分。
更小比例尺的地圖(如 1:50,000, 1:100,000, 1:250,000, 1:500,000, 1:1,000,000)則是由 1:25,000 的圖幅合併而成。
目前我國基本地形圖的主體是 1:25,000 的圖幅。
其分幅規則通常是:
- 1:1,000,000 地圖:劃分為 100 張圖幅。
- 1:500,000 地圖:由 1:1,000,000 圖幅細分。
- 1:250,000 地圖:由 1:500,000 圖幅細分。
- 1:100,000 地圖:由 1:250,000 圖幅細分。
- 1:50,000 地圖:由 1:25,000 圖幅細分。
- 1:25,000 地圖:作為基本圖幅,通常是 7.5' 經度 x 5' 緯度。
【答案】
我國基本地形圖採用經緯度分幅與方里網細分結合的方式。以 1:25,000 圖幅為基本單元,其範圍約為 7.5 分鐘經度 x 5 分鐘緯度。
(2). 此種圖幅經台灣TM2投影後之圖面圖幅框線之形狀為何?
台灣 TM2 投影 (Transverse Mercator Projection, Zone 2):
台灣地區的 TM2 投影採用 橫軸墨卡托投影。
橫軸墨卡托投影的基本特徵是:
- 將地球投影到一個與其赤道面相切的圓柱面上,然後將圓柱面展開成平面。
- 對於橫軸墨卡托投影,圓柱體的軸與地球的自轉軸垂直,圓柱體的母線與地球的自轉軸平行。
- 它能保持投影帶內的形變最小,特別是對於長度、面積和角度的變形。
圖幅框線形狀:
經過橫軸墨卡托投影後,由於投影方式的特性,圖幅的邊界(框線)在投影後的平面上通常是 矩形。
然而,由於地球是球體(或橢球體),而地圖是平面,即使是橫軸墨卡托投影,其投影帶外的圖幅邊緣也會有一定程度的變形。但對於標準的投影帶內的圖幅,其邊緣(經緯線)在展開後會近似於矩形。
更精確地說,由於投影帶的邊界是經線,而圖幅的頂部和底部邊界是緯線,在展開成平面時,經線變為平行線,緯線也變為平行線,但經線之間的距離在緯線上會隨緯度變化,緯線之間的距離在經線上也會隨經度變化。
在 TM2 投影帶內,經線投影後是平行直線,緯線投影後也是平行直線。因此,圖幅的框線(由經緯線構成)在投影後近似為 矩形。
在實際的地圖上,圖幅的邊緣通常會被繪製成矩形。