日本放射線技術学会中四国部会 2014 夏季学術大会 画像情報研究会

ノイズの発生メカニズムと低減
および利用に関する画像技術

香川高等専門学校 電気情報工学科
本田道隆

今日お話する内容
１．画像に重畳する構造的ノイズとその発生原因、および一般的な対策
・

あまり本質的なものではないが、故障や調整不良のシステムで画像に障害陰影を与えるノイズと対策

２．ランダム性を有するノイズの発生原因と理論的なアプローチ（分析法）
・・

放射線像に本質的に重畳するノイズ、特に量子モトルと熱雑音
ノイズ理論の復習と推奨する計測法や分析法

３．ノイズ低減手法と利用法について
・・

システム側の低減手法の紹介、および画像処理による低減法
ノイズは低減すべきものだが、利用することもある１．画像に重畳する構造的ノイズとその発生原因、および一般的な対策

ラスターノイズ（古い映像装置での）走査線が画像に写り込む
走査型映像機器
電子銃

ターゲット

ラスターノイズ（古い映像装置での）走査線が画像に写り込む
走査型映像機器
電子銃

近年はビームスキャン方式が姿を消しており、このタイプ
ターゲットのノイズはほとんどなくなった。
しかし、万一古い装置などで低減させる必要があるときは、ノッチ型フィルタ（後述）がある程度効果的。

いわゆる不均一性ノイズ（近年の映像系に多い）

多重の読み出し機構などで各ラインのゲイン補正にズレが発生したとき
CCD素子など

偶数列用
ゲイン 100
ゲイン 99
このノイズは、現にFPDで発生する。
ゲイン 100
ゲイン 99
そこで、均一な被写体を撮影し、このノイズを収集して実奇数列用ブロックに分けた
際の画像にはその逆のパターンを補正として使用する。
読み出しや処理ノッチ型フィルタが有効な場合もある。

いわゆる不均一性ノイズ（近年の映像系に多い）

多重の読み出し機構などで各ラインのゲイン補正にズレが発生したとき
CCD素子など

偶数列用
ゲイン 100
ゲイン 99奇数列用

ゲイン 100

ゲイン 99ブロックに分けた
読み出しや処理いわゆる不均一性ノイズ（近年の映像系に多い）

多重の読み出し機構などで各ラインのゲイン補正にズレが発生したとき
CCD素子など

偶数列用
ゲイン 100
ゲイン 99
このノイズは、現にFPDで発生する。
ゲイン 100
ゲイン 99
そこで、均一な被写体を撮影し、このノイズを収集して実奇数列用ブロックに分けた
際の画像にはその逆のパターンを補正として使用する。
読み出しや処理ノッチ型フィルタが有効な場合もある。

ノッチ型フィルタとは？

実例

ある特定の周波数成分だけを画像上から除くように設定したフィルタ
周波数変調度

重み（実空間フィルタ係数）
1.0
0.8
0.6
0.4
0.2
0.0-0.2

1.0

最大周波数成分を除去

0.0-2.0-1.002.0[lp/mm]

周波数2 lp/mm の周期成分を除去

1.0-5 -4 -3 -2 -1 012345

画像上の相対画素位置

周波数変調度

重み（実空間フィルタ係数）

1.0

0.8

0.0-0.2

0.0-2.0-1.0

周波数01.0[lp/mm]

2.0-5 -4 -3 -2 -1 0123

画像上の相対画素位置

45

外来（飛び込み）ノイズ

走査中に外来ノイズが画像に混入する
走査型映像機器
電子銃

このノイズは、微小な信号検出の際に、モータなどの電磁ターゲット波が検出回路に影響を与える場合などに生じる。（現在はEMI対策などがルールになっており十分対策が進んでいる）
代表的対策としてメディアンフィルタがある。

故障によるノイズと緊急措置

平面検出器の場合、信号線やゲート線が縦横に走行している（故障時には線状影が発生しやすい）
信号線

何らかの異常ゲート線よほどの事故（検出器が壊れるほどの強い衝撃を受けるなど）がなければ現在の製造技術では考えにくいが、万一発生した場合には目立たせないような緊急措置は可能。
信号が正常に出てこない
これも、メディアンフィルタが最も単純で簡単な対応。

メディアンフィルタとは？

実例

メディアンフィルタ ＝ カーネル内の中央値を出力するフィルタ

カーネルサイズ：

3×3のとき上から順に９個の画素値を並べると{ 38, 31, 25, 24, 16, 15, 10, 9, 7 }となる。中央値は 5 番目に大きい値なので16がフィルタ出力（これまでの24を16に置き換える。）

29182513729381631930724153222192591033484302917

中央値なのでノイズが平滑化される。
エッジをぼかさない。
細線を消せる。（消えてしまう）
3 × 3 画素領域16が出力される
カーネル内に細い線があってもメディアンフィルタで除かれる(パーフェクトではないが）

２．ランダム性を有するノイズの発生原因と理論的なアプローチ（分析法）
・ノイズの発生原因と少しばかりの理論的復習
・回路系ノイズの基礎的な説明・計測法に関する補足ランダムノイズの２大発生原因

ショット雑音と熱雑音

ショット雑音（ショットノイズ）

ターゲット

ショットノイズ（電子がばらばらに発生）

透過確率e ‐L

熱雑音（ジョンソンノイズ）（サーマルノイズ）

電荷移動を阻止する作用がある抵抗器内における電子の不規則運動

熱雑音（いわゆる電子回路のノイズ）

電圧のばらつきショットノイズ = ポアソンノイズ (1)
入力粒子

L[cm]

No.１

5個のうち1個透過する確率は？
1個が透過する確率 ： P14個が透過しない確率 : (1-P)4

No.2
No.3

だから
145C1 × P × (1-P)

No.4
No.5

透過する確率 P

こういうのを「二項過程」や「ベルヌーイ試行」などと言います

この要領で、0個から5個まで計算すると
0.7
0.6
0.5
0.4
0.3
0.2
0.1
0

P=0.1
P=0.2

012

P=0.3

345

ショットノイズ = ポアソンノイズ (2)X線の場合は、多量の電子がターゲットに当たり、しかも電子 → X線となる確率も小さい （約１％）ターゲット

１％程度多くの電子（例えば 500個）

多くの電子X線

このような場合、出てくるX線はポアソン分布に従う
確率(ヒストグラム）

L[cm]

No.１

出てくるX線光子の平均 λ(この数値例の場合はλ=5)No.2No.499No.500X線になる確率0.01

出てくるX線光子数

多くの入力 (n 個）が 確率 p で 出力されるときの ポアソン分布発生した事象数が k 個である確率
期待値 λ = np

P (k ) ke k!

分散 σ2 = np

期待値と分散が等しいことは大事な事実
λ が大きくなれば正規分布に近づいてくることも重要な事実
たとえば透視の場合1画素にX線光子（フォトン数）が平均 9 個吸収されると、分散も 9 。
だから、標準偏差は 3 なのでノイズの大きさは 3 となる。
ただし、間接変換では少しボケるのでノイズはその分少なくなる。

画像の平均値からノイズが分かるので条件制御などで、インテリジェンスな制御もできる。

ショットノイズ の伝搬
被写体上の1×1mm面積に入るX線光子
平均 40,000個とすると標準偏差は 200個

SNR = 200

平均 400個
標準偏差は 20個

SNR = 20

平均 100個
標準偏差は 10個

SNR = 10

被写体（透過率 0.01）

検出膜（吸収率 25%）

間接変換ではボケ
平均 100個
標準偏差は 5 個

見かけ上のSNR = 20

平均 100個
標準偏差は 6 個

見かけ上のSNR = 16.7

読み出し回路のノイズ

ショットノイズ の伝搬
ポアソン分布に従う統計的性質
被写体上の1×1mm面積に入るX線光子
平均 40,000個とすると標準偏差は 200個
被写体（透過率 0.01）
平均 400個
標準偏差は 20個
検出膜（吸収率 25%）
平均 100個
標準偏差は 10個

SNR = 200
ここは二項過程による統計的作用が発生するが、ポアソン分布に従う性質は保たれる

SNR = 20
細かい被写体ではボケの影響を受けるので、周波数毎に DQE(f) でSNR = 10
評価されることが多い

間接変換ではボケ
平均 この3者をまとめて検出器の性能
100個
見かけ上のSNR= 20を量子検出効率(DQE) で表し、標準偏差は5個X線の利用率を数値化する
読み出し回路のノイズ
平均 100個
標準偏差は 6 個

見かけ上のSNR = 16.7

熱雑音 = 電子ノイズ = 回路ノイズ (1)電圧のばらつき抵抗器から発生する！

電圧のばらつき(ノイズ）V 

4 kTR  fk : ボルツマン定数
T : 絶対温度
R : 抵抗値
f : 周波数帯域

絶対温度が高く、抵抗値が大きいと熱雑音は多く発生する。

昔、冷却したカメラを使った映像装置なども試作されたことがありました。

熱雑音 = 電子ノイズ = 回路ノイズ (2)

知っておいていただきたい
回路ノイズの発生原因

熱雑音そのものは神の産物ともいえるUncontrollableなもの。
ここが完全に 0 [V] (グラウンド電位)
理想
これをいかにして画像に伝搬しないようにするかが技術。
であれば微小電荷でも完璧に検出できる
フォトダイオードa-Se
映像系技術に関しては、
チャージ・浮遊容量を減らすこと
アンプ・信号の量を極力増やすこと＜1画素＞などが著名な知見であり、鉄則ともいえる

実際

安定

抵抗から発生した電位変化のため電荷の一部がひっかかる

フォトダイオード
a-Se

＜1画素＞

信号線とGND線の浮遊容量
2V

ノイズあり
大きな抵抗(ノイズ源）

21

熱雑音 = 電子ノイズ = 回路ノイズ (3)浮遊容量を減らすと言っても・・
実寸大の検出器は苦労が多い

実寸サイズの長い配線
小面積の場合大面積の場合抵抗値大抵抗値小容量小

容量大どうしても配線がクロスになる
金属と金属が向き合うのでコンデンサを形成
データ（画素電荷）が流れる信号線
読み出しの選択線（ゲート線）

回路で発生するノイズの影響の変化
ランダムノイズ
刷毛ではいたようなノイズ

現在、臨床に使えるようになったが極低線量部は少し残ることがある

経験的に回路ノイズに「補正」はなく、本質的に低減しなければうまくいかない回路ノイズが影響する例（刷毛ではいたようなノイズ）
信号線ゲート線（読み出し行の選択線）

電荷の出入りOFF状態ではゲート線上のノイズは同じ波形

水色の信号電荷に同量のノイズ電荷（赤色）が各行の画素に一斉に重畳する現在は対策が進んでいるため、この例のような劣化画像を出力することはないが、極低線量の部分などで存在を指摘されることはあり得る。
画素濃度（濃度で正規化してないパワースペクトル）
600500垂直方向プロフィールのパワースペクトル

400300水平方向プロフィールのパワースペクトル200
100013579 11 13 15 17 19 21 23 25

空間周波数 [lp/cm]

ランダムノイズの計測について

・標準偏差・ウィナースペクトル（パワースペクトル）
この両者は品質管理などでは水平方向と垂直方向のそれぞれで計測するとより正確な管理ができる。
（画質や処理の研究でも場合によっては必要）
背景トレンド除去の代わりに、画像を2枚収集して減算する手法でもかなり良い結果が得られる。
（ノイズは 約 1.4 倍になっていることに注意）
標準偏差は簡便で大変有益な指標だが、これ以外にノイズの周波数特性も計測することを勧める。

ランダムノイズの周波数特性は画質評価で極めて重要

低周波成分多い
Contrast 0.5
CNR2.0高周波成分多い
Contrast 0.5
CNR2.0３．ノイズ低減手法と利用法について
・回路ノイズの低減は一般には Uncontrollable？
・画像処理によるノイズ低減を中心に・ノイズを利用する技術について

ランダムノイズ低減の基本手法
１．フレーム間平均効果による低減（リカーシブフィルタ）

出力画像

ｘ a-n

ｘ a-3

ｘ a-2

画像に重みをかけて加算平均するとノイズの分散は(a02 + a-12+ +a-n2) 倍になる。

ｘ a-1

ｘ a0+

例： 4枚の加算平均で 各重みを 0.25 とすると分散は 0.252 × 4 倍 = 0.25 倍になる。
ノイズ量 (標準偏差) では 0.5 倍になる。

1 フレーム遅延×透視像
Xk

1-a係数a＋×係数 a = 0.6 の場合

Yk = 0.6 Xk + 0.4Yk-1
Yk = 0.6 Xk + 0.24 Xk-1 + 0.096 Xk-2 + ・・・

ランダムノイズ低減の基本手法
２．画素同士の平均効果による低減（線形空間フィルタ）

平滑化されたノイズ

30

ランダムノイズ低減の基本手法
２．画素同士の平均効果による低減（線形空間フィルタ）

平滑化されたノイズ

31

線形フィルタのサンプル
原画

空間（平滑化）フィルタ

線形フィルタのサンプル
原画

リカーシブフィルタ（係数 0.6）

ランダムノイズ低減の基本手法の不満な点
リカーシブフィルタ（単純フレーム平均）線形空間フィルタ
（平滑化フィルタ）ランダムノイズは低減
するが、動く信号のコントラストも低下させてしまう

ランダムノイズは低減
するが、細かい信号のコントラストも低下させてしまう

信号とノイズを分けて個別に処理する必要がある
・先見情報に基づいてパターン認識などで信号と背景と弁別するタイプ・ノイズと信号の一般的性質の差を検知しながらフィルタを変えるタイプ
これらは、「エッジ保存型ノイズ低減処理」と呼ばれる

実際に応用されている処理例
観察対象を限定せず、信号と背景を局所の統計的性質で弁別する処理
コヒーレントフィルタ
原画

統計的性質が近い画素同士で加算平均する

単純なリカーシブフィルタ

35実際に応用されている処理例

１（フィルタの効果）

透視像のノイズ比較
ノイズ低減処理なし

新しいノイズ低減フィルタ

ノイズ（相対値） [%]120

画像の明るさによらずノイズを約40％にすることができた。

10080
60
40200Lung

Heart

Media

Abd

1

ノイズ低減処理なし
新しいノイズ低減フィルタ

MTF0 .8

このフィルタを適用してもMTFの劣化はなかった。

0 .60 .4
0 .2
000 .5

1

1 .5

22 .5

3

Spatial frequency [m m - 1 ]36

現在研究途上の処理例
線状信号が観察しにくい透視画像（シミュレーションで作成）

CNR :
線状信号

矩形領域のテンプレートによるマッチング信号

統計検定量 (T値)

1.3
背景領域：

CNR :

ノイズの平滑化

2.6

背景領域：

CNR :

1.3（原画のハイパス処理)×0.7

CNRは同じでも、背景処理の方法で画質改善する期待があることが確認された

現在研究途上の処理例２

原画

背景ノイズを高周波化してさらにノイズ量を約 70%に抑えた結果

なぜノイズが細かいと線状信号が見やすくなるのか？
~ Statistical Decision Theory Model（SDTM）からの説明 ~◆大きい信号を見る場合◆小さい信号を見る場合信号の周りのノイズは、広い範囲で平滑化して見ている。

信号の周りのノイズは、狭い範囲で平滑化し、細かいノイズとして見ている。

「信号を見るとき知覚の過程で、信号と同じ周波数特性をもつフィルタをノイズに自然にかけている」 と考えると大変うまく説明できる。
実際、この考え方は放射線画像分析でよく使われる。

Statistical Decision Theory Model （統計的決定理論モデル）
信号

視覚系の周波数特性

 S (f)  V (f)df 22

SNR SD 22 N(f)  V (f) S (f)  V (f) df
2

ノイズのパワースペクトル22

ここに信号特性がフィルタのように入る

SNRSD が同じ画像例

右と左のどちらの画像が線が見やすいでしょうか？

なぜノイズが速く動くと線状信号が見やすくなるのか？
ヒトの視覚系は速い動きを正しいレベルで感知できない。
つまり、時間方向に平滑化フィルタをかけて知覚する。

明滅１辺サイズ画像上のランダムノイズを低減させるには？

・ノイズ量（ノイズの標準偏差）を低減させる。
物理的な低減。
普通、ノイズ低減とはこのことを指す。

・ノイズの周波数特性を相対的に増強する。
知覚的な低減。
画像技術に通用するノイズ低減の一方式。

・ノイズをダイナミックに高速で動かす。
知覚的な低減。
画像技術に通用するノイズ低減の一方式。

ノイズの利用 （一般的なもの）

工学的には、量子化ノイズの低減法や、いわゆる心地よい 1/f ゆらぎを形成したインターフェース技術などが知られている。

同一信号が継続的に発生するような場合、ノイズがなければデジタル化するときに誤差が生じる。

ランダムノイズを重畳させて計測を何度も行って平均すると、量子化レベル以下の精度で信号がデジタル化できる。

ノイズの利用 （X線技術特有）

X線のフォトン数がポアソン分布に従うことを利用
期待値 λ = np

分散 σ2 = np

画像レベルを計測すると、ノイズ量が推定できる！
検出器

透視

ある場所の SNR = 10 だった。
MTFでノイズが少し減っていることを補正すると、X線光子数は 50個/画素
DSA では管電圧を下げてコントラストをあげ、さらに SNR を10 倍にしたい

X線管

DSA

70kVで光子数を 5000個 にすればよいので、テーブルをひいて mAsを決めれば良い！！

ノイズの利用 （画像処理の例）

パターン認識には誤認もあるアーチファクトを目立たせないようにして、オペレータに違和感を与えないようにすることは重要。

ノイズの少ない画像

適度なノイズがある画像

ご清聴ありがとうございました。