資料-2

精度管理への不確かさの導入

不確かさの表現ガイドの概要
Guide to the Expression of Uncertainty in Measurement:
GUM

日本測量機器工業会編「ISO17123‐1 測量機器の現場試験理論改定概要」
より引用

従来の誤差評価の問題点• 誤差評価では、「誤差＝測定値 真の値」と定義。実際に「真の値」を求めることは困難。
• 誤差評価においては、偶然誤差や系統誤差を総合的に評価するための合理的な合成方法が定められていない。
• 測定結果を表現するための用語の定義が国・地域や専門分野によって異なる場合がある。

GUMのISO標準への導入• 測量分野では「ISO17123-1測量機器の現場試験手順：理論（第2版）」で導入• 精度の表現：不確かさ• 標準不確かさのタイプA評価• 標準不確かさのタイプB評価• 不確かさ評価バジェット• 統計的検定法

用語と定義① 不確かさ• 計測のばらつきを示す用語として、不確かさ（Uncertainty）を導入し、従来の規格で扱っていた偶然誤差だけでなく、系統誤差を含む誤差の包括的な表現として定義する。

② 標準不確かさ• 正規分布する偶然誤差のばらつきの尺度である標準偏差を、適当な分布を与えることで系統誤差にまで拡張したもの。不確かさバジェットの要素となる。評価法によってタイプAとタイプBに分けられる。

③ 合成標準不確かさ• 測定結果にまつわる全ての標準不確かさ（誤差）を誤差の伝搬則によって足し合わせたもので、測量成果の不確かさを表す。

④ 拡張不確かさ• 標準偏差（σ）に対して2σあるいは3σなどで表現される誤差範囲を、標準不確かさに適用したもの。

⑤ 不確かさバジェット• 測定結果にまつわる全ての不確かさ要素を表形式に整理したもので、最終的に合成標準不確かさあるいは拡張不確かさが見積もられる。

標準不確かさ• 標準不確かさのタイプA評価• タイプAの標準不確かさ（u ）は、統計的評価を用いて求められ、実験標準偏差s で表される。

• 標準不確かさのタイプB評価• タイプBの標準不確かさ（u）は、統計的評価以外の方法で求められるもので、以下のような様々な情報を利用した科学的判断に基づいて評価される。
• 以前の測定データ• 一般的知識及び経験• 製造業者の仕様• 校正データ• 想定される様々な誤差分布に応じた標準偏差• 一様（矩形）分布：u＝a／√12, aは矩形の幅

合成標準不確かさと拡張不確かさ• 合成標準不確かさ• 標準不確かさの2乗和の平方根で計算• 拡張不確かさ• 合成標準不確かさに包含係数kを掛けて計算• k=2は従来の概念の誤差範囲2σに相当：約95%の信頼水準

平面位置座標の不確かさの表現方法

平面位置座標の不確かさをどう表すか• 測量、測位分野において平面位置の「精度」の表示方法が不明
確という問題• 各次元（X, Y）ごとに1次元の標準偏差で表示 ←• 各次元の標準偏差の合成で表示

明確• 合成の式が不明確• 合成した標準偏差を誤差円と仮定しているか不明• 誤差円と仮定している場合でも、誤差円内の確率分布に無頓着

誤差円の考え方• 1次元Gauss分布の場合、確率変数が±1σ内に存在する確率は68％、±2σ内に存在する確率は95％• 2次元Gauss分布の場合、確率変数が1σの円内に存在する確率は39％、1.510σでは68％、2.4477σでは95％
誤差円の半径誤差円内に落ちる確率（%）

1σ391.1774σ

50

√2σ

63

1.510σ

682.1460σ

902.4477σ

9510

作業規程の準則における表し方• 各水平座標成分ごとの標準偏差: Mx, My• 1次元の標準偏差→68％信頼区間• 鉛直成分は1次元なので同じく68％信頼区間• 座標の標準偏差: Ms =√(Mx2+My2) ←誤差円とは書いていない• 2次元の標準偏差であり、誤差円の半径を表していると考えると→ Mx≒ MyのときMs =√2Mxなので信頼水準63%の円• 信頼水準95%の円の半径は、1.7308Ms

FGDC（米国連邦地理データ委員会）のGeospatialPositioning Accuracy Standards (1998)• 水平座標の位置正確度 positional accuracy for horizontal coordinates• 真値（理論値）の95%がその中に落ちる不確かさの円の半径• 鉛直座標の位置正確度 positional accuracy for vertical coordinates• 真値（理論値）の95%がその+/-の幅の中に落ちる不確かさの区間

測角の視準の不確かさ
機器検定データから導いた結果

13

セット間較差の標準偏差と視準の標準偏差の関係• １視準の測定値θiの標準偏差：𝜎 （観測方向、望遠鏡正反、対回、セットによらないと仮定） の 目• １夾角(θ2 –θ1)当たりの標準偏差： 2𝜎• １夾角1対回観測（望遠鏡正反）平均値の標準偏差：

方盛
向ゼロ=𝜎

• １夾角N対回観測平均値の標準偏差：• 観測平均値の2セット間の較差の標準偏差：

θ1θ2

≡𝜎14

セット間較差から推定した1視準当たりの標準偏差• 検定データのセット間較差は絶対値のみで表される
１級• 絶対値のセット間較差の平均 最小目盛値 10"読み5"読み1"読み TS
値𝑥̅ を用いて（符号付）セッ セット間較差
1.16 1.10 0.74 0.59
ト間較差の標準偏差𝜎 の推定 平均値
値𝑠 =𝑥̅ を求める。

• （符号付）セット間較差の標準偏差から1視準当たりの標準偏差σsを推定

セット間較差によるσsの推定値

1.8

1.7

1.10.9≡𝜎
15

1視準当たりの測角誤差の構成要素• 目盛分解能 → 最小目盛値• 視準誤差• 望遠鏡の分解能 → レーリーの分解能• 肉眼の分解能• 目標物に合わせる個人的能力

16

最小目盛値と測定の不確かさの関連付け• 最小目盛1区間の確率分布を一様（矩形分布）とする（タイプBの標準不確かさ）

xi-1

xi+1
最小目盛値：αr=xi-xi-1xi• 最小目盛値に起因する不確かさ：ur＝αr／√12• αr=1”のときur ＝0.289”• αr=5”のときur ＝1.44”17望遠鏡分解能と測定の不確かさの関係付け
αA
αA/2• 視準の不確かさはレーリーの分解能よりも小さいはずだが、その関係は？
• レーリーの分解能（エアリーディスクの半径）：αA＝1.22λ／Dから、D=45mmのときαA＝3.1” （λ=555nm）• 視準誤差がαA/2を超えることはまれと想定• ±αA/2の範囲内の視準誤差の確率分
布は矩形分布を仮定• その仮定の下でαA=3.1”のときuA＝ αA /√12 =0.89”18

肉眼の分解能の影響• 肉眼の分解能は約60秒と言われており、倍率30倍の望遠鏡を通すと約2秒に拡大される
望遠鏡の分解能3.1秒よりも小さいので不確かさバジェットには加算しない
望遠鏡分解能と肉眼分解能は2つの独立した誤差ではなく、2つのシステム関数の畳み込みで評価。
エアリーディスクの大きさはシステム関数の畳み込みでは変わらない。

• 目標物に合わせる個人的能力の差は検定
データでは不明

19

目盛分解能と望遠鏡分解能の合成から計算した不確かさとセット間較差から推定した視準の不確かさの比較
最小目盛値

目盛値起因の 分解能起因の合成不確かさ
不確かさ
不確かさ

αr (")ur (")

検定データから求めた不確かさ

ur+A=√(ur2+uA2) ss,S (")

uA (")

10.29

0.89

0.94

1.1

5

1.44

0.89

1.69

1.7

• 目盛分解能と望遠鏡分解能を合成した不確かさは、
1秒読みTS：セット間較差から推定した不確かさより小さい。他の不確かさ要因がある。
5秒読みTS：セット間較差から推定した不確かさとほぼ等しい（大小について判定できない）。目盛分解能と望遠鏡分解能が主要因。
20

測角の視準の不確かさ
野外測定実験データから導いた結果

21

夾角の測定の標準偏差• 基線場測定実験から得られた夾角の測定データを処理すること
で、セット間較差からの推定ではなく、夾角測定の標準偏差を直接計算できる。

• 2方向×正反観測（2回）×3対回×2セット×10回測定×3台＝720標本• 機器（TS）を交換すると夾角の値が変わるので、各機器ごとに各方向ごとの夾角の平均値を求め、その平均値からの偏差を分析• 推定した1夾角の観測の標準偏差は3.5秒である。

22

夾角の偏差の分布（全720標本）
集計
250
200
150
100
500集計

-10-8-6-4-2024681039808763197

116

135

2802

夾角の偏差 (秒)標準偏差：3.50秒
23夾角測定の標準偏差から推定される測定値
野外測定実験に基づく1夾角測定値の標準偏差である3.5秒をもとに推定すると、• 1視準当たりの測定の標準偏差は 3.5”/√2=2.5”• 検定データの絶対値セット間較差（ 1170標本）から推定した
1視準当たりの測定の標準偏差は1.7秒であり、野外測定から求めた値の2.5秒は明らかにそれより大きい。

24• 両者の差が測定者による目標の視準の不確かさに起因すると仮定すると、その不確かさの大きさは約1.8”である。
つまり、 √(2.52-1.72)=1.8• 機器検定においては目標物が細線であり、レチクルとの合致が容易だが、野外測定では目標物が反射鏡であり、レチクルとの合致が困難であることが原因と推測。
• 2重レチクル線の中間に目標物を視準する精度について、研究論文が１篇(*)しか見つからず、視準の不確かさが1.8”との推測の妥当性は確認できない。
(*)田中・矢野(1979):「精密測定における2重レチクル線合致の誤差と判断機構」, 人間工学, 15, 315-32325

目標の見え方の違い
十字線との合致前十字線との合致後機器検定：目標（コリメータの十字線）が明瞭

野外測定：目標（11m先の反射鏡中心）が不明瞭
26

倍角差の許容範囲を定める根拠

1夾角測定の標準偏差から倍角差の制限
値を推定• 倍角差Bは３つの倍角（bj, j=1,2,3）の範囲Rangeである• bj = rj +lj ，rj：j対回目の望遠鏡正（右）での夾角，lj ： j対回目の望遠鏡反（左）での夾角
• 𝜎 𝑟 = 𝜎 𝑙 = 2𝜎 , 𝜎 𝑏 = 2𝜎𝜎 は1視準の標準偏差• bj|max : j個の倍角のうちの最大値。 bj|min : j個の倍角のうちの最小値。
• B= bj|max - bj|min• bjの任意抽出ではないので、𝜎 𝐵 = 2𝜎(𝑏)とはならない• 範囲によるQCにはJIS Z 9020-2（シューハート管理図）を適用• 標準偏差 𝜎の測定値の範囲R の期待値はE 𝑅 = 𝑑 ∗ 𝜎• 範囲Rの上方管理限界は𝑈 𝑅 = 𝐷 ∗ 𝜎• Rの平均値にR の標準偏差の3倍を加えたものに相当： R の上限値とみなす
28

検定における倍角差の許容範囲•𝑈

𝑅 = 𝐷 ∗ 𝜎の関係を適用• 倍角の標準偏差𝜎(𝑏)から倍角差の許容範囲𝑈
定める

𝐵 = 𝐷 ∗ 𝜎(𝑏)を• 検定においては𝜎 = 1.7“、𝜎 𝑏 = 3.4“、検定における対回数n=3のとき、D2 = 4.358• 𝑈 𝐵 ≑ 15"• 機器検定から得られた倍角差の最大値20”• 倍角差15”以内が99.8%
準則の制限値30秒は過大
倍角差の許容範囲を15”とすることは妥当
29

基準点測量における倍角差の許容範囲• 野外測定においては𝜎 𝑠 = 2.5“、𝜎 𝑏 = 5.0“、基準点測量の対回数n=2ならばD2 =3.686 、またはn=3ならばD2 = 4.358• 𝑈 𝐵 ≑ 18“ または𝑈 𝐵 ≑ 22“
基線場測定から得られた倍角差の最大値15”（全体の約6%）基準点測量における観測値の倍角差の許容範囲30” は過大
𝑈 𝐵 に近い20”とするのが妥当。

• 較差の標準偏差と観測差の許容範囲についても同様に計算可能• 教科書や準則の示すところでは観測差は倍角差より小さいが、検定及び野外実験のデータでは両者に差があるという結果は得られなかった• 品質管理に観測差と倍角差の両方が必要か、未検討

30