ホーム » Taiki Bessho の投稿 (ページ 13)
作者アーカイブ: Taiki Bessho
ノイズの計算 (準備)
ノイズについてより深く説明します.
#説明しながら数学も少し使います.
キーワード
- 測定
- ノイズ量を表す指標
- 分散,標準偏差
まとめ
- 測定値に含まれる誤差がノイズ
- ノイズの量は分散と標準偏差によって表される.
- 天体写真がザラザラするのは標準偏差がシグナルに対して大きく,S/Nが小さいから
一度ノイズをより深く理解すると,非常に有益な多くの実用的な知識が芋づる式に分かります.例を挙げればきりがありませんが,例えば以下のことが分かるようになります.
- S/Nについて理解が深まるため,天体写真でコンポジットを行うことが何故有効なのか
- フラット画像を作る際は非常に多くの枚数をコンポジットする必要がある
- ダーク減算をすると実は天体写真の背景ノイズが大きくなる場合がある
ノイズを理解するために,まずノイズの量を表す指標である分散と標準偏差について説明します.
いきなり天体写真の場合について考えると難しいので,鉛筆の長さを定規で測る例を使って考えていきます.
長さが150.0mmだと既に分かっている鉛筆をものさしで計ったところ,誤差により150.1mmと計れてしまったとしましょう.
この誤差0.1mmは「ものさしの当て方」,「目盛りを読み取るときのズレ」などが主な原因でしょう.これらは偶然発生したもので,再現性はないので偶然誤差(ノイズ)です.
再現性がないため次にもう一度計ると誤差の量が変わって今度は150.3mmや149.8mmに計れることもあるでしょう.
この例では鉛筆の長さ(真の値)が150.0mmであることがあらかじめ分かっていましたが,普通何かを測定する場合は真の値が分からないことがほとんどです.したがって測定値から真の値がどのあたりにありそうかを推定することになります.
そのためには測定値にどのくらいの誤差が含まれているか知る必要があります.
測定値に含まれる誤差の大きさを表す指標として分散と標準偏差というものがあります.それぞれ順に説明していきます.
(ギリシャ文字は「シグマ」と読みます.)
先ほどの鉛筆をものさしで連続回計ることを考えましょう.第回目の測定値をとします.つまりはと何回目の測定かを表します.回の測定の平均値をと書きます.つまりです.
(ギリシャ文字は「ミュー」と読みます).
このとき分散は次で定義されます:
(1)
そして標準偏差はその平方根で定義されます:
(2)
文字のままだとイメージがわきにくいので,具体的に3回測定した場合()で計算してみましょう.
最初の測定値が, 2回目が, 3回目がだったとします.このとき平均値はです.
これから分散を求めると:
(3)
と計算出来ます.
(計算式中で平均値の値として150.1ではなく150.07を用いたのは計算の精度を上げるためです)
標準偏差を求めるにはこの平方根を取って:
(4)
となります.
測定値とその誤差を表す際に今求めた標準偏差を使って
と書けます.
真の値150.0mmがちゃんと標準偏差の範囲内で含まれているので,この3回の測定は悪くない測定だったと言うことが出来ます.
さて,数式や具体的な計算をしながら分散と標準偏差を説明しました.ここで特に覚えておいて欲しいのは
「測定値に含まれる誤差の大きさを表すのに分散と標準偏差が使われる」
ということです.
天体写真では写真の滑らかさを表すのにS/Nが使われますが,なのです.
写真がノイズでザラザラに見えるのはシグナルに対して標準偏差が大きいからだと言えます.
ここで説明しているのは厳密には「標本平均」「標本分散」「標本の標準偏差」です.
本来ならば有限回の試行回数でものさしの測定精度を見積もるためには「不偏分散」「標準偏差の不偏推定量」を考えるべきですが,ここでは話を簡単にするため「標本分散」「標本の標準偏差」で代用しました.)
最近のコメント