ホーム » 2016 » 9月 (ページ 3)

月別アーカイブ: 9月 2016

ノイズの計算（準備）

2016年9月28日 / コメントする

ノイズについてより深く説明します．
＃説明しながら数学も少し使います．

キーワード

測定
ノイズ量を表す指標
分散 $\sigma^2$ ，標準偏差 $\sigma$

まとめ

測定値に含まれる誤差がノイズ
ノイズの量は分散 $\sigma^2$ と標準偏差 $\sigma$ によって表される．
天体写真がザラザラするのは標準偏差 $\sigma$ がシグナル $S$ に対して大きく，S/Nが小さいから

一度ノイズをより深く理解すると，非常に有益な多くの実用的な知識が芋づる式に分かります．例を挙げればきりがありませんが，例えば以下のことが分かるようになります．

S/Nについて理解が深まるため，天体写真でコンポジットを行うことが何故有効なのか
フラット画像を作る際は非常に多くの枚数をコンポジットする必要がある
ダーク減算をすると実は天体写真の背景ノイズが大きくなる場合がある

ノイズを理解するために，まずノイズの量を表す指標である分散 $\sigma ^2$ と標準偏差 $\sigma$ について説明します．

いきなり天体写真の場合について考えると難しいので，鉛筆の長さを定規で測る例を使って考えていきます．

アセット 1 長さが150.0mmだと既に分かっている鉛筆をものさしで計ったところ，誤差により150.1mmと計れてしまったとしましょう．
この誤差0.1mmは「ものさしの当て方」，「目盛りを読み取るときのズレ」などが主な原因でしょう．これらは偶然発生したもので，再現性はないので偶然誤差(ノイズ)です．
再現性がないため次にもう一度計ると誤差の量が変わって今度は150.3mmや149.8mmに計れることもあるでしょう．

この例では鉛筆の長さ(真の値)が150.0mmであることがあらかじめ分かっていましたが，普通何かを測定する場合は真の値が分からないことがほとんどです．したがって測定値から真の値がどのあたりにありそうかを推定することになります．

そのためには測定値にどのくらいの誤差が含まれているか知る必要があります．

測定値に含まれる誤差の大きさを表す指標として分散 $\sigma ^2$ と標準偏差 $\sigma$ というものがあります．それぞれ順に説明していきます．
(ギリシャ文字 $\sigma$ は「シグマ」と読みます．)

先ほどの鉛筆をものさしで連続 $N$ 回計ることを考えましょう．第 $i$ 回目の測定値を $X_i$ とします．つまり $i$ は $i=1,2,\dots,N$ と何回目の測定かを表します． $N$ 回の測定の平均値を $\mu$ と書きます．つまり $\mu = \frac{X_1+X_2+\cdots+X_N}{N}$ です．
(ギリシャ文字 $\mu$ は「ミュー」と読みます)．

このとき分散 $\sigma^2$ は次で定義されます：

(1) $\begin{align*} \sigma ^2 = \frac{(X_1-\mu)^2 + (X_2-\mu)^2 + \cdots + (X_N -\mu)^2}{N} \end{align*}$

そして標準偏差 $\sigma$ はその平方根で定義されます：

(2) $\begin{align*} \sigma = \frac{\sqrt{(X_1-\mu)^2 + (X_2-\mu)^2 + \cdots + (X_N -\mu)^2}}{\sqrt{N}} \end{align*}$

文字のままだとイメージがわきにくいので，具体的に３回測定した場合( $N=3$ )で計算してみましょう．
最初の測定値が $X_1 = 150.1\text{mm}$ , ２回目が $X_2 = 150.3\text{mm}$ , ３回目が $X_3 = 149.8\text{mm}$ だったとします．このとき平均値は $\mu=\frac{150.1+150.3+149.8}{3}=150.067 \fallingdotseq 150.1\text{mm}$ です．
これから分散 $\sigma^2$ を求めると：