厳密高低網平均の計算方法
このページは、網平均計算について興味のある方のために具体的な計算方法を掲載しました。
「厳密高低網平均の計算方法」
(上記ファイル名をクリックするとダウンロード可ですが、厳密高低網平均の計算方法エクセルの改良はご遠慮ください)
(Space Netを操作する上では、本ページ内容と「厳密高低網平均の計算方法」エクセルを理解する必要はありません)
測量計算ソフトSpace Net につきましては以下のページをご覧ください。
「Space Net XY・高低網」、「Space Net 三次元網」、「Space Net 水準網」、
「Space Net 自由網(XY・高低)」、「Space Net 自由網(三次元)」
シート名 | 種類 | 説明 | エクセルに 使用した主な計算式 | |||
計算の内容 | 計算式 | 記号の説明 | エクセル使用箇所 | |||
入力1 | 入力 | ・厳密高低網平均に必要なデータを入力します。 ・与点標高と新点標高近似値を入力します。 ・新点標高近似値(H2)は右の計算式より求めます。 |
新点標高の近似値 | H2=(H2正方向+H2反方向)/2 H2正方向=H1+Dsinα1+i1-f2+K H2反方向=H1-Dsinα2-i2+f1-K |
H1:既知点標高 H2:新点標高近似値 i:器械高 f:目標高 D:測定距離(斜距離) S:基準面上の距離 Z:観測した鉛直角 α:高低角 α=90°-Z K:両差(気差及び球差) K=(1-k)S2/(2R) k:屈折係数(0.133) R:平均曲率半径6370000 |
本仕組みファイルは、左記の計算式の結果を入力しますので左記の計算式は使用していません。 |
入力2 | 〃 | ・厳密高低網平均に必要なデータを入力します。 ・各測線について、厳密高低網平均結果の基準面上の球面距離を入力します。 |
||||
入力3 | 〃 | ・厳密高低網平均に必要なデータを入力します。 ・各測線(正・反の両方)について、観測値(器械高、目標高、高低角)を入力します。 |
||||
1高低角補正 | 計算 | ・高低角の補正をしています。 ・入力3シートで入力した高低角を標石上面の高低角に換算します。最初に正反観測値別々に、器械高・目標高・基準面上の距離から、観測高低角を標石上面の高低角に換算します。次に正と反の値を平均した標石上面の高低角(網平均計算に使用)を求めます。 |
高低角の補正 | @正の高低角に対する補正量 dα1=tan-1((f2-i1) cosα1/(S/cosα1-(f2-i1)sinα1)) A反の高低角に対する補正量 dα2=tan-1((f1-i2) cosα2/(S/cosα2-(f1-i2)sinα2)) B補正した観測高低角(標石上面) α’1(正)=α1-dα1 α’2(反)=α2-dα2 C観測高低角(正反平均) α’m=(α1-α2)/2 |
α:観測した高低角 dα:Aに対する補正量 α’:補正後の観測高低角 α’m:補正後の観測高低角(正反平均) i:器械高 f:目標高 S:基準面上の距離 |
@H〜X列 AAD〜AT列 BY〜AC列(正) AU〜AY列(反) CAZ〜BE列 |
2定数項・係数 | 〃 | ・観測方程式の定数項と係数を求めています。 ・各測線毎に、定数項と係数を求めます。 ・定数項は、与点標高と近似標高と基準面上の距離等から求めた近似高低角と、正反平均の標石上面の観測高低角との差です。 ・係数C1とC2は、近似高低角と、基準面上の距離と標高から求めます。 |
定数項 係数C1、C2 |
@近似標高により求めた高低角 A’=tan-1((H’2-H’1)/S・(1-(H’1+H’2)/(2R))) A観測方程式の定数項 L=A’-α’m B観測方程式の係数 C1=cos2α’/S・(1-H’1/R)ρ″ C2=cos2α’/S・(1-H’2/R)ρ″ |
A’:近似高低角 H’:近似標高 R:平均曲率半径6370000 L:定数項 C:係数 ρ″:206265 |
@G〜Q列 AS〜X列 BZ〜AA列 |
3方程式 | 〃 | ・観測方程式の係数行列を求めています。 ・各測線毎に、該当する係数C1・C2を「2定数項・係数」シートから呼び込んでいます。出発点は係数に-1をかけます。行列の縦が各測線、横が新点です。与点も新点と同様に扱います。 ・係数C1、C2を呼び込むため、各測線の出発点、到着点の区別が重要になります。 |
観測方程式の係数行列 | V(α)=-C1凾1+C2凾2-L 上記が観測方程式ですが、本シートでは、係数行列のみ表示しています(定数項Lは「2定数項・係数」シート参照) |
V(α):高低角の観測方程式 凾:未知量(近似標高の補正量) |
D〜J列 |
4行列 | 〃 | ・行列計算により最小二乗法で近似標高の補正量等を求めています。 ・観測方程式の各要素(定数項ベクトル、係数行列、重み行列)から最小二乗法による行列計算で近似標高の補正量、観測値の残差を求めます。(ベクトルは1列) ・行列計算A’PVが0であれば正しく計算されたことになります。(点検) |
行列計算 | A’PA=N A’PL=U -N-1・U=W AW+L=V |
A:係数行列 N:正規行列 U:定数項ベクトル W:補正量ベクトル V:残差ベクトル P:重み対角行列 -1:逆行列 |
全体 |
5標高最確値・標準偏差 | 〃 | ・標高の最確値、観測値と標高の標準偏差を求めています。 ・新点標高近似値に行列計算で求めた補正量を足して標高の最確値を求めます。 ・観測値の残差、重み、方程式数、未知数から単位重みの標準偏差を求めます。 ・正規行列の逆行列N-1の対角要素と単位重みの標準偏差から、新点標高の標準偏差を求めます。(測量精度の一尺度) |
@標高の最確値 A観測値の標準偏差 B標高の標準偏差 |
@H=H’+凾 Aσ0=√(1/(m-n)ΣPVV) Bσ0・√q |
H:標高最確値 H’:標高近似値 凾:近似標高の補正量 σ0:単位重みの標準偏差 V:観測値の残差 P:重み=1 m:観測方程式数 n:未知数(新点数) q:正規方程式の逆行列の対角要素 |
@G〜I列 AB〜D列 BJ〜K列 |
6高低角最確値・残差 | 〃 | ・高低角の最確値と、高低角の残差を求めています。 ・最確値の標高と距離から高低角の最確値を求めます。 ・最確値高低角から、観測高低角を引いて高低角の残差を求めます。(行列計算の残差ベクトルと一致します) |
@高低角の最確値 A高低角の残差 |
@高低角の最確値=tan-1(最確値高低差/平均標高上の最確値距離) A高低角の残差=高低角の最確値-観測高低角(α’m) |
@E〜P列 AQ〜Z列 |