狠狠撸

多次元信号処理の基礎
と画像処理のための
二次元変換技術
平成２３年９月３０日（金）
新潟大学工学部電気電子工学科
准教授村松正吾

自己紹介
（むらまつしょうご）
? 所属
? 新潟大学工学部電気電子工学科
? 画像?映像信号処理の教育研究に従事
? 学部：「プログラミング」，「画像情報工学」
? 修士：「画像処理特論」，博士：「多次元信号処理論」
? 学会：IEEE，電子情報通信学会，映像情報メディア学会
? 著書など
? 「マルチメディア技術の基礎DCT入門」
（CQ出版社, 1997年）
? 「MATLABによる画像&映像信号処理」
（CQ出版社, 2007年）

講演内容
? 一次元信号の分析と合成
? 画像変換の基礎
? 重複変換とウェーブレット変換
? 非分離変換処理
? 応用と今後の発展

信号処理とは？
? 音や光等の物理的変量を伝送?記録?解析?加工する技術
電気信号（電波?光信号）
センシング
音?光?超音波?地震?X線など
再生
（必要に応じて）

一次元信号処理
? 例：音（モノラル）
? 一変量の関数とみなせる
?(?)
? 音の標本化
? 数列（サンプル列）が得られる
?[?]
? 音の解析
? 「どのような成分から構成されているか？」
フーリエ解析など
１秒間に440回振動

多次元信号処理
? 複数の変量を持つ信号
? 画像?映像
? 多チャンネル信号
? 画像の標本化（モノクロ）
? 連続分布 → 配列
? ?0, ?1 → ?[?0, ?1]
? 画像の解析
? 「どのような成分から構成されているか？」
多次元フーリエ解析など
高
高

準備：遅延器
? ???
：?サンプル分時間シフト
3
1
3
1
5
3
3
1
3
1
5
3
3
1
3
1
5
3
3
1
3
1
5
3
1
1
??1 ??1 ??2
0 0 0
1 1 2

信号の解析(一次元）
? 近似成分と詳細成分に分解
3
1
3
1
5
3
2 2
4
2
3
1 1 1
-1 -2
3
1
×1/2
×1/2
+
4 4 4 6 8
+
隣同士を足す
右隣を引く
近似成分
詳細成分
2 -2 2 -4 2
-1
4
-

信号の合成
? 分解した信号を足し合わせると…
+
元に戻る！
2 2
4
2
3
1 1 1
-1
-2
近似成分
詳細成分
3
1
3
1
5
3
-1
4
足し算、引き算は
信号の情報を逃さない
立派な成分分析法！
低周波成分
高周波成分

フィルタバンク表現
??1
+
1/2
1/2
遅延器
??1
+
-1/2
1/2
+
入力出力
近似成分
詳細成分
分析器合成器
群（集まり）の意
?0 ? =
1
2
1 + ??1
?1 ? =
1
2
?1 + ??1
? ? = ?0 ? + ?1 ? = ??1
全体の伝達関数?(?)が単なる遅延

分析と合成の交換
入力
分析器
+
出力
合成器
??1
+
1/2
1/2
+
-1/2
1/2
各処理を合成器に移動しても元に戻る
??1
?0 ? =
1
2
1 + ??1
?1 ? =
1
2
?1 + ??1
? ? = ?0 ? + ?1 ? = ??1
全体の伝達関数?(?)が単なる遅延

その他の構成
??1
+
1/2
1/2
??1
+
-1/2
1/2
入力
近似成分
詳細成分
分析器
+
出力
合成器
同じ分析器に対し異なる合成器も有り得る
（合成器）（分析器）
??1
+
1/2
1/2
+
1/2 ??1
-1/2
?0 ? =
1
2
1 + ??1
?1 ? = ±
1
2
?1 + ??1
?0 ? =
1
2
1 + ??1
?1 ? = ?
1
2
?1 + ??1
? ? = ?0 ? ?0 ? + ?1 ? ?1 ?
=
1
4
1 + 2??1 + ??2
±1 ?1
遅延！

信号の分析と合成
入力出力
遅延???
になればよい！
分析器合成器
絵の具に例えると…
サンプリング密度が増加する（冗長となる）
量が増える！？

無駄のない信号の解析
? 無駄を省くための間引きも構成も可能
3
1
3
1
5
3
2 2
4
1 1 1
3
1
×1/2
×1/2
+
4 4 4 6 8
－
縮小近似成分
縮小詳細成分
2 -2 2 -4 2
絵の具に例えると…
2点に１点棄却
右隣を引く
隣同士を足す

無駄のない信号の合成
? 分解した信号を補間して足し合わせると…
+
元に戻る！
2 2
4
1 1 1
縮小近似成分
縮小詳細成分
3
1
3
1
5
3
低周波成分
高周波成分
2 2
4
2 2
4
1 1 1
そのまま補間
符号反転して補間 -1 -1 -1

最大間引きフィルタバンク表現
??1
+
1/2
1/2
??1
+
-1/2
1/2
+
入力出力
縮小近似成分
縮小詳細成分
分析器合成器
↓2
↓2
↑2
↑2
無駄がない（非冗長）
ダウンサンプラーアップサンプラー
??1
+
1
1
+
1
-1 ??1

? 定義
? 例
ダウンサンプラ－の定義
↓?
? ?nx ? ?my ? ? ? ?Mmxmy ?
?サンプル毎に１サンプルを出力
? ? 1個のサンプルを棄却
↓2
2 2
4
2
3
4
2 2
4
×
×
×

アップサンプラ－の定義
? 定義
? 例
↑?
? ?nx ? ?my
? ? ? ?
??
?
?
? ???
?
その他0
2,,0 ?MMm
M
mx
my
各サンプル間に
? ? 1個の零値を挿入
↑2
2 2
4
2 2
4
0 0 0

補間処理の例
? 繰り返し補間（ゼロ次ホールド）
? 符号反転補間
2 2
4
1 1 1
2 2
4
2 2
4
1 1 1
-1 -1 -1
↑2
??1
+
1
+
1
-1 ??1
2 2
4
1
↑2
1 1 1 左隣を引く
隣同士を足す

分析と合成の行列表現
（等価変換）
??1
+
1/2
1/2
??1
+
-1/2
1/2
+
入力出力
分析器合成器
↓2
↓2
↑2
↑2
??1
+
1
1
+
1
-1 ??1
本当に単なる遅延になるのか？
（完全再構成できるのか？）

（等価変換）
??1
+
1/2
1/2
??1
+
-1/2
1/2
+
入力
出力
分析器合成器
↓2
↓2
↑2
↑2
??1
+
1
1
+
1
-1 ??1

（等価変換）
??1
+
1/2
1/2
??1
+
-1/2
1/2
+
入力
出力
分析器合成器
↓2
↓2
↑2
↑2
??1
+
1
1
+
-1
1 ??1
遅延器を一つにまとめられる遅延器を一つにまとめられる

（等価変換）
+
1/2
1/2
??1
+
-1/2
1/2
+
入力
出力
分析器合成器
↓2
↓2
↑2
↑2
??1
+
1
1
+
-1
1
二つの内積演算→行列演算二つの内積演算→行列演算

（等価変換）
??1
+
入力
出力
分析器合成器
↓2
↓2
↑2
↑2
??11
2
１１
1 ?1
１１
1 ?1
捨てる値の計算は不要零値の計算も不要
1
3
2
-1
0
0
0
0

（等価変換）
??1
+
入力
出力
分析器合成器
↓2
↓2
↑2
↑2
??1
? =
1
2
１１
1 ?1
? =
１１
1 ?1
?0
?1
?0
?1
?0
?1
?0
?1
= ?
?0
?1
= ??
?0
?1
=
１ 0
0 1
?0
?1
=
?0
?1
予め
捨てる
後で零値
挿入
等しい

ブロック化とデブロック化
??1
+
入力
出力
↓2
↓2
↑2
↑2
??1
3
1
3
1
5
3
本当に単なる遅延になるのか？
（完全再構成できるのか？）
↑
0
3
1
3
1
5
3
↑
0
1 1 3
3 3
5
1 1 3
3 3
5
1 1 3
3 3
1 1 3
5
↑
0
単なる遅延！

分析処理のブロック処理表現
（最大間引きフィルタバンク）
? 2点ブロック毎の行列演算に帰着
3
1
3
1
5
3
2 2
4
1 1 1
縮小近似成分
縮小詳細成分
2
1
=
1
2
1 1
1 ?1
3
1
2
1
=
1
2
1 1
1 ?1
3
1
4
1
=
1
2
1 1
1 ?1
5
1
?
?
?
?
?
?

間引き構成のブロック処理表現
（合成処理）
? 係数ブロック毎の逆行列演算に帰着
元に戻る！
2 2
4
1 1 1
縮小近似成分
縮小詳細成分
3
1
3
1
5
33
1
=
1 1
1 ?1
2
1
3
1
=
1 1
1 ?1
2
1
5
3
=
1 1
1 ?1
4
1
?
?
?
?
?
?

信号変換の応用例
（冗長系も非冗長系も）
? 圧縮符号化
?JPEG, JPEG2000
?MPEG, H.264/AVC
? ノイズ除去
?ウェーブレット縮退法
? 特徴抽出
?Viola-Jonesの顔検出
分析合成
：：
縮退
分析合成
：：
量子化出力入力
入力出力

音声圧縮符号化
? ATRAC3(MD)
? MP3/AAC(iPod)
1. 人の耳に聞こえにくい音を分析
2. 音の成分を分別して捨てる！
約1/10など
CD MD
※http://www.sony.co.jp/Products/ATRAC3/wa.htmより引用
周波数
大
き
さ
約1,400kbps 132kbpsなど

画像圧縮符号化
? JPEG/JPEG2000（静止画）
? MPEG2, H.264/AVC（動画）
1. 人の目に見えにくい成分を分析
2. 画像の成分を分別して捨てる！
信号解析
1
2
3
4
5
6
7
8
0
2
4
6
8
-200
0
200
400
600
800
1000
1200
放送カメラ DVD/Blu-ray

1秒辺りのビット数R(bps) ＝
画素数×画面数/秒×ビット数/画素
? ハイビジョン品質（HDTV）
R = (1920×1080)×30×24 ≒ 1.5 Gbps
? アナログ放送品質（SDTV）
R = (720×480)×30×24 ≒ 248 Mbps
地デジ
1.5Gbps
248Mbps
約14Mbps
圧縮前
地デジ
圧縮前
約4Mbps

自然画像の性質
??
?
?
??
?
?
?
1
0
x
x
x
1x
隣接画素の散布図

信号変換の効果
? 隣接値との相関を減らすことができる
? スパース表現が得られる
??
?
?
??
?
?
?
1
0
x
x
x
0x 1x
0x
1x 1y
0y
??
?
?
??
?
?
??
?
?
??
?
?
?
???
?
?
??
?
?
1
0
1
0
11
11
2
1
x
x
y
y
1y
0y
広がり広がり広がり
相関が減る
スパース化される

変換行列の例
? ハール変換（π/4回転）行列
??
?
?
??
?
?
?
?
11
11
2
1
A ??
?
?
??
?
? ?
??
11
11
2
11
A
IAA ??1

信号変換の基礎
? 基底ベクトル
? ?110
1
,,, ?
?
? MbbbA ?
111100 ?????? MMyyy bbbx ?
信号xは基底ベクトルの線形結合として表現できる
逆行列の各列ベクトルを基底ベクトルと呼ぶ

? 基底ベクトルの例
? 係数は順変換により求められる
??
?
?
??
?
?
??
?
?
?
?
?
?
?
??
?
?
??
?
?
?
???
?
?
??
?
?
?
0
2
2
2
11
11
2
1
1
0
y
y
y
2 2 2 2
22
? ? ?
?
?
?
?
?
?
?
???
?
?
??
?
??
????
?
?
??
?
?
????
?
?
??
?
?
??
?
?
??
?
? ?
???
?
?
??
?
?
?
2
2
1
1
2
1
0
1
1
2
1
2
11
11
2
1
101
0
1
0
yyy
y
x
x
x
0b 1b
基底ベクトルの線形結合である
22

? 座標上での関係
1x
0x
0y
1y
2
2
22
2
3
0
-1
??
?
?
??
?
??
????
?
?
??
?
?
???
?
?
?
?
?
?
?
1
1
2
1
1
1
1
2
1
3
22
2
? ? ??
?
?
??
?
??
?????
?
?
??
?
?
???
?
?
?
?
?
?
?
1
1
2
1
1
1
1
2
1
3
2
22
0b1b
1
??
?
?
??
?
??
????
?
?
??
?
?
???
?
?
?
?
?
?
?
1
1
2
1
0
1
1
2
1
2
2
2
変換係数は新たな軸上での座標

基底ベクトルの調べ方
? 変換係数を一つのみ１、残りを０として逆変換
? ハール変換の例
,
0
0
1
1
0
?
?
?
?
?
?
?
?
?
?
?
?
?
?
? ?
?
Ab ,
0
1
0
1
1
?
?
?
?
?
?
?
?
?
?
?
?
?
?
? ?
?
Ab
?
?
?
?
?
?
?
?
?
?
?
?
?
?
? ?
?
1
0
0
, 1
1
?
? AbM
,
1
1
2
1
0
1
11
11
2
1
0 ??
?
?
??
?
?
???
?
?
??
?
?
??
?
?
??
?
? ?
?b ??
?
?
??
?
??
???
?
?
??
?
?
??
?
?
??
?
? ?
?
1
1
2
1
1
0
11
11
2
1
1b

基底ベクトルの調べ方
+↑2
↑2
??1
A?1 =
1
2
１ ?１
1 1
+↑2
↑2
等価
??1
+
1 2
??1
+
1 2
?1 2
1 2
1
2
1
2
1
2
1
21
0
1
0
1
2
1
2
?0 ?
?1 ?
合成フィルタ??(?)のインパルス応答は，?番目基底ベクトルの要素である
インパルス応答

画像変換の基礎
? 基底画像（ベクトルから画像へ）
画像Xは基底画像の線形結合として表現できる
1,11,11,01,00,00,0 1010 ???????? MMMMyyy BBBX ?
基底画像

? 基底画像の例
1/2 1/2
1/2 1/2
０２
４６
-1/2 1/2
-1/2 1/2
-1/2 -1/2
1/2 1/2
1/2 -1/2
-1/2 1/2
=
6 × 2×
4× 0×
1,01,0 B?y0,00,0 B?y
X
0,10,1 B?y 1,11,1 B?y
変換係数基底画像

? 分離処理
? 順変換?逆変換の分離処理の行列表現
T
AXAY ?
T??
? YAAX 1
順
逆
yn
xn
変換係数を求める
元の配列に戻す
yk
xn
? = ??
? = ?? ? ?
= ?? ?
※画像処理では
直軸を下向きにとる

? の2次元ハール変換
? 逆変換は以下のとおり
??
?
?
??
?
?
???
?
?
??
?
? ?
??
?
?
??
?
?
??
?
?
??
?
?
?
???
?
?
??
?
?
??
04
26
11
11
64
20
11
11
2
1
64
20 TT
AAAXAY
????????????????????
1,10,11,00,0
10
00
0
01
00
4
00
10
2
00
01
6
04
26
BBBB
AAAAAAAA
AAYAA
??
?
?
??
?
?
????
?
?
??
?
?
????
?
?
??
?
?
????
?
?
??
?
?
??
??
?
?
??
?
?
?
TTTT
TT
０２
４６

基底画像の調べ方
? 変換係数を一つのみ１、残りを０として逆変換
??
?
?
??
?
?
?
?
???
?
?
??
?
?
???
?
?
??
?
?
??
?
?
??
?
? ?
???
?
?
??
?
?
?
??
?
?
??
?
? ??
???
?
?
??
?
?
???
?
?
??
?
?
??
?
?
??
?
? ?
???
?
?
??
?
?
?
??
?
?
??
?
?
?
?
???
?
?
??
?
?
???
?
?
??
?
?
??
?
?
??
?
? ?
???
?
?
??
?
?
?
??
?
?
??
?
?
???
?
?
??
?
?
???
?
?
??
?
?
??
?
?
??
?
? ?
???
?
?
??
?
?
?
11
11
2
1
11
11
10
00
11
11
2
1
10
00
11
11
2
1
11
11
01
00
11
11
2
1
01
00
11
11
2
1
11
11
00
10
11
11
2
1
00
10
11
11
2
1
11
11
00
01
11
11
2
1
00
01
1,1
0,1
1,0
0,0
AAB
AAB
AAB
AAB
T
T
T
T

水平シフタ?垂直シフタ
? ?y
??：?画素分下シフト
? ?y
+?：?画素分上シフト
? ?x
??： ?画素分右シフト
? ?x
+?： ?画素分左シフト
?y
?1
?y
+1
?x
?1
?x
+1

画像変換のフィルタバンク表現
?y
?1
+
?y
?1
+
入力
↓2×1
1 2
1 2
?1 2
1 2
↓2×1
?x
?1
+
?x
?1
+
↓1×2
1 2
1 2
?1 2
1 2
↓1×2
?x
?1
+
?x
?1
+
↓1×2
1 2
1 2
?1 2
1 2
↓1×2
垂直ダウン
サンプラー水平ダウン
サンプラー

画像逆変換のフィルタバンク表現
+
出力
?y
?1
+
1 2
1 2
+
?1 2
1 2 ?y
?1
↑2×1
↑2×1+
?x
?1
+
1 2
1 2
+
?1 2
1 2 ? ?
?1
↑1×2
↑1×2
+
?x
?1
+
1 2
1 2
+
?1 2
1 2 ?x
?1
↑1×2
↑1×2
水平アップ
サンプラー
垂直アップ
サンプラー

基底画像の調べ方
↑2×1+
↑1×2 ?0 ?x
↑1×2 ?1 ?x
↑1×2 ?0 ?x
↑1×2 ?1 ?x
+ ↑2×1
?0 ?y
?1 ? ?
+
合成フィルタ??(? ?) ??(?x)のインパルス応答は，?, ?番目基底画像の要素である
1
2
1
1
2
1
2
1
2 1
2
1
2

? × ?行列ブロック変換処理
（一次元）
??1
+
入力
出力
分析器合成器
↓?
?
?0
?1
?0
?1
完全再構成条件：??=?
↓?
↓?
?
?
?
??1
?
?
? ?
?
?
? ??1
? ??1
??1
??1
?
?
?
↑?
↑?
↑?
?
?
?
?
?0
?1
? ??1
+

ブロック変換のイメージ
? ブロック処理手順（４×４行列の例）
? ?nx ? ?nx?
E
E
E
R
R
R

離散コサイン変換(DCT)
? DCT(Type-II)の特徴
?変換核が余弦（コサイン）関数で与えられる
?直交変換である
?基底ベクトルが対称性をもつ
?フーリエ解析と密接な関係がある
?高速アルゴリズムが存在する
?準最適変換である
JPEG,MPEG-1/2, H.264/AVCに採用されている

? DCT-IIの定義（一次元）
? ? ? ? ? ? 1,,1,0,
2
12
cos
2 1
0
???
?
?
?
?
? ?
? ?
?
?
Mm
M
mn
nxk
M
my
N
n
m ?
?
? ? ? ? ? ? 1,,1,0,
2
12
cos
2 1
0
???
?
?
?
?
? ?
? ?
?
?
Mn
M
mn
myk
M
nx
N
m
m ?
?
順変換
逆変換 ??
?
?
?
??
?
?
1,,2,11
0
2
1
Mm
m
km
?

? 基底ベクトルと基底画像（８点DCT）
0 2 4 6 8
-0.5
0
0.5
p0
0 2 4 6 8
-0.5
0
0.5
p2
0 2 4 6 8
-0.5
0
0.5
p4
0 2 4 6 8
-0.5
0
0.5
p6
0 2 4 6 8
-0.5
0
0.5
p1
0 2 4 6 8
-0.5
0
0.5
p3
0 2 4 6 8
-0.5
0
0.5
p5
0 2 4 6 8
-0.5
0
0.5
p7

? 画像の離散コサイン変換
サイズ８×８のブロック処理例

JPEGの概要
? ベースライン方式の符号化手順
RGB→YUV
ブロック化
2D-DCT
量子化
可逆
DPCM
ジグザグ
スキャン
ハフマン符号化
二次元
ハフマン符号化
多重化
DC
AC
前処理変換符号化エントロピー符号化ビット列生成
量子化テーブルハフマンテーブル
入
力出
力
符
号

? 正方行列（DCTなど）の問題点
ブロックノイズモスキートノイズ
変換係数を棄却?量子化したときの問題

より一般的な?等分割最大間引き
フィルタバンク（一次元）
+
入力
出力
分析器合成器
?0
?1
↓?
↓?
↓?
?
?
?
?
?
?
?
?
?
? ??1
?
?
?
+
↑?
↑?
↑?
?
?
?
?0 ?
?1 ?
? ??1 ?
?0 ?
?1 ?
? ??1 ?
インパルス応答の長さは?点を超えてもよい

4タップ２等分割最大間引き
フィルタバンク
? ?30h
? ?10h
? ?00h
? ?20h
2
? ?31h
? ?11h
? ?01h
? ?21h
? ?00f
? ?01f
? ?20f
? ?30f
? ?10f
? ?21f
? ?31f
? ?11f
2
2
2 2
2
2
2
? ? ? ? ? ?
? ? ? ? ? ?
? ? ? ? ? ?
? ? ? ? ? ?
??
135
024
113
202
xxx
xxx
xxx
xxx
?
?
? ? ? ? ? ??? 012 xxx
? ? ? ? ? ?
? ? ? ? ? ?
??
101
101
111
000
?
?
yyy
yyy
2x4行列演算 4x2行列演算 ? ?
? ?1
1
3
2
?
?
u
u
?
? ?
? ?
? ?
? ?0
0
0
0
3
2
1
0
u
u
u
u
? ?
? ?
?
1
1
1
0
u
u
? ? ? ? ? ??? 012 xxx
?
?
重複保持重複加算
? ?
? ?
? ?
? ?3
2
1
0
x
x
x
x
?
?
1?
z
1?
z
1?
z
1?
z
1?
z
1?
z
隣の隣の隣まで計算

重複ブロック変換のイメージ
? 重複ブロック処理手順
? ?nx ? ?nx?
P
P
P
Q
Q
Q

重複ブロック処理の例
（一次元，4タップ２分割）
??1
+
入力
出力
分析器合成器
↓2
?0
?1
?0
?1
完全再構成は可能か？
↓2
↓2
??1
?3
??1
↑2
↑2
↑2
?0
?1
?3
??1
↓2
?2
? ? ??1
+
↑2 +
??1
?2
2x4行列演算 4x2行列演算

??1
入力
分析器合成器
↓2
?0
?1
↓2
↓2
??1
?0
?1
?3
??1
↓2
?2
?0
?1
+
出力
?0
?1
?3
??1
↑2
↑2
↑2
??1
+
↑2 +
??1
?2
?0
?1
+
+
? = ?0 ?1 ? =
?1
?0

??1
入力
↓2
??1
?0
?1
??1
↓2
?0
?1
+
出力
?0
?1
??1
↑2
↑2
??1
??1
?0
?1
?0
?1
+
+
+
+
分析器合成器

??1
入力
↓2
?0
?1
↓2
E ?
=?0+?1 ??1
+
出力
?0
?1
??1
↑2
↑2
R ?
=?0+?1 ??1
?0
?1
分析器合成器
記憶をもった行列演算となっている
? ? ? ? = ?0+?1 ??1
?0+?1 ??1
= ???
?完全再構成条件：

4タップ2分割フィルタバンクの例
? 2次Daubechies直交フィルタバンク(DB2)
1 1.5 2 2.5 3 3.5 4
-0.2
0
0.2
0.4
0.6
0.8
1
1.2
Reconstruction low-pass filter
1 1.5 2 2.5 3 3.5 4
-0.5
0
0.5
1
Reconstruction high-pass filter
The synthesis filters for db2
基底ベクトル
P =
0.4830 0.8365 0.2241 -0.1294
-0.1294 -0.2241 0.8365 -0.4830
Q =
0.4830 -0.1294
0.8365 -0.2241
0.2241 0.8365
-0.1294 -0.4830

? × ?行列重複ブロック変換処理
（一次元）
??1
+
入力
出力
分析器合成器
↓?
?(?)
?0
?1
?0
?1
完全再構成条件：? ? ? ? =??? ?
↓?
↓?
?
?
?
??1
?
?
? ?
?
?
? ??1
? ??1
??1
??1
?
?
?
↑?
↑?
↑?
?
?
?
?(?)
?0
?1
? ??1
+

重複直交変換（LOT/GenLOT)
? 16タップ8等分割LOT
0 5 10 15
-0.5
0
0.5
n (samples)
振幅
Impulse 応答
0 5 10 15
-0.5
0
0.5
n (samples)
振幅
Impulse 応答
0 5 10 15
-0.5
0
0.5
n (samples)
振幅
Impulse 応答
0 5 10 15
-0.5
0
0.5
n (samples)
振幅
Impulse 応答
0 5 10 15
-0.5
0
0.5
n (samples)
振幅
Impulse 応答
0 5 10 15
-0.5
0
0.5
n (samples)
振幅
Impulse 応答
0 5 10 15
-0.5
0
0.5
n (samples)
振幅
Impulse 応答
0 5 10 15
-0.5
0
0.5
n (samples)
振幅
Impulse 応答
基底ベクトル基底画像

離散ウェーブレット変換（DWT）
? ツリー構成（不等分割最大間引きフィルタバンク）
? ?
? ? ? ?
? ? ? ?
? ?42
1
21
0
0
0 zHzHzH
? ?
? ?zH 0
1
? ?
? ?zF 0
1
8?
8?
2?
8?
8?
2?
? ?
? ? ? ?
? ?21
1
0
0 zHzH 4? 4?
? ?
? ? ? ?
? ? ? ?
? ?42
0
21
0
0
0 zHzHzH
? ?
? ? ? ?
? ? ? ?
? ?zFzFzF 0
0
21
0
42
1
? ?
? ? ? ?
? ?zFzF 0
0
21
1
? ?
? ? ? ?
? ? ? ?
? ?zFzFzF 0
0
21
0
42
0
等価
? ?
? ?zH 0
0
? ?
? ?zH 0
1
2?
2?
? ?
? ?zH 1
0
? ?
? ?zH 1
1
2?
2?
? ?
? ?zH 2
0
? ?
? ?zH 2
1
2?
2?
? ?
? ?zF 0
0
? ?
? ?zF 0
1
2?
2?
? ?
? ?zF 1
0
? ?
? ?zF 1
1
2?
2?
? ?
? ?zF 2
0
? ?
? ?zF 2
1
2?
2?

DCTとDWTの比較
? 多重解像度解析
DCT(JPEGなど) DWT(JPEG2000など)
垂直周波数
水平周波数
低
低
高
高

DCTとDWTの比較
? DWTがノイズに強い理由
周
波
数
低
高
DCT ウェーブレット(DWT)
周波数に応じて
ブロックサイズが異なる
ブロックサイズが
互いに重複している
耐ブロックノイズ
耐モスキートノイズ

多次元信号処理の基礎と応用
非分離
自由度
大
簡便性大小
小
可分離
可分離では直交性、対称性、重複性を同時に満たす２ｘ２等分割は存在しない
可分離システムでは斜め方向の指向性は可分離では実現できない
LOT DirLOT

フーリエ解析の復習
? 信号は（複素）正弦波の足し合わせ
? ?tx
t
t
t
t
t
どの様な正弦波で構成されているか？
? ??jX
?
MAX?
（複素）正弦波の周波数分布
最大角周波数
周波数スペクトル

? 標本化（連続時間）
一次元標本化とスペクトル
t
? ?tx? ?txT
t t
＝ ×
? ??
?
???
?
n
nTt?
TT
t?? t?
? ??jXT
?
? ??jX 標本化 ? ?
T
jX ?
?0 t?t??
：大T
エイリアジング

? 標本化（連続時空間）
二次元標本化とスペクトル
1?
0?
? ?10 , jΩjΩXV
0
2
?
?
1
2
?
?
スペクトルは周期的となる
1?
0?
? ?10 , jΩjΩX
0
2
?
?
1
2
?
?

寛容な多次元標本化
? ?ΩjX
0
2
P
?
1
2
P
?
1?
0
2
P
?
1
2
P
?
1?
0?
? ?ΩU jX
0?
標本化前標本化後
??
?
?
??
?
?
?
1
0
0
0
P
P
V
?
?
?
?
?
?
?
?
?
?
?
?
?
1
0
2
0
0
2
P
P
?
?
U
各次元独立に標本化定理を
適用してはならない！
単純な
直交標本化

エッジ画像のスペクトル例
画像
振幅特性
フィルタリングには周波数サポートの考慮が欠かせない
局所的方向に適した基底画像の必要性がいえる
0P
?
1P
?

可分離システムと非分離システム
↑2×1+
↑1×2 ?0 ?x
↑1×2 ?1 ?x
↑1×2 ?0 ?x
↑1×2 ?1 ?x
+ ↑2×1
?0 ?y
?1 ? ?
+
+
↑2×2 ?0 ?y, ?x
↑2×2 ?1 ?y, ?x
↑2×2 ?0 ?y, ?x
↑2×2 ?1 ?y, ?x
+
+非分離
可分離
直接２ｘ２分割を行う

指向性非分離変換
? 非分離指向性フィルタの導入
水平
DCT
垂直
DCT
垂直
DCT
垂直
DCT
垂直
DCT
垂直
IDCT
垂直
IDCT
垂直
IDCT
垂直
IDCT
水平
IDCT
帯域形状
水平
垂直
変換
水平
垂直
逆変換
帯域形状
可分離基底
非分離基底

非分離/指向性GenLOT
入力
出力
分析器合成器
↓?
?(?y, ?x)
?0
?1
?0
?1
完全再構成条件：? ?y, ?x ? ?y, ?x =???
?
↓?
↓?
?
?
?
?
?
?
? ??1
? ??1
↑?
↑?
↑?
?
?
?
?(?y, ?x)
?0
?1
? ??1

ノイズ除去応用
? ウェーブレット縮退法
分析合成
：：
縮退入力出力
原画像ノイズ画像可分離DWT 指向性ODWT

冗長系とスパース表現
? アトムと辞書
合成
：
出力
原画像ノイズ画像可分離DWT 指向性DWT
合成
：
合成
：

画像処理のための冗長系の例
? 無指向性変換の例
?非間引き（シフト不変）ハール変換／DCT
? 指向性変換の例
?Curvelet
?非間引きContorlet(NSCT)
?混成DirLOT

まとめ
? 一次元信号の分析と合成
? 画像変換の基礎
? 重複変換とウェーブレット変換
? 非分離変換処理
? 応用と今後の発展
?冗長系
?局所的適応制御

狠狠撸

多次元信号処理の基础と画像処理のための二次元変换技术

Recommended

More Related Content

What's hot (20)

Similar to 多次元信号処理の基础と画像処理のための二次元変换技术 (20)

More from Shogo Muramatsu (6)

多次元信号処理の基础と画像処理のための二次元変换技术