• 検索結果がありません。

4.講義内容

N/A
N/A
Protected

Academic year: 2021

シェア "4.講義内容"

Copied!
20
0
0

読み込み中.... (全文を見る)

全文

(1)

2次元フーリエ変換

講義内容

講義内容

空間周波数の概念

2次元フーリエ変換

代表的な2次元フーリエ変換対

2次元離散フーリエ変換

(2)

フーリエ変換と逆変換

u

v

F.T.

∫ ∫

−∞∞ ∞ ∞ − − + = f x y j ux vy dxdy v u F( , ) ( , )exp{ 2π ( )}

連続系

連続系

離散系

離散系

∑∑

− = − = + − = 1 0 1 0 } / ) ( 2 exp{ ) , ( 1 ) , ( N x N y N vy ux j y x f N v u F π

x

y

)

,

(

x

y

f

)

,

( v

u

F

I. F.T. ) , ( vu F ただし,ここで は絶対値を とって画像化

∑∑

− = − = + = 1 0 1 0 } / ) ( 2 exp{ ) , ( 1 ) , ( N x N y N vy ux j v u F N y x f π

順変換

順変換

逆変換

逆変換

(3)

2次元フーリエ変換の具体的なイメージ

∑∑

− = − =

+

=

1 0 1 0

}

/

)

(

2

exp{

)

,

(

1

)

,

(

N i N j

N

vy

ux

j

y

x

f

N

v

u

F

π

}

/

)

(

2

exp{

j

π

ux

+

vy

N

)

,

(

x

y

f

対応する画素ごとに積をとって 最後に総和をとる. 対応する画素ごとに積をとって 最後に総和をとる.

はどんなパターンか?

それでは

exp{

j

2

π

(

ux

+

vy

)

/

N

}

離散系での説明

(4)

2次元フーリエ変換の具体的なイメージ

)

(

2

sin

)

(

2

cos

)}

(

2

exp{

j

π

ux

+

vy

=

π

ux

+

vy

j

π

ux

+

vy

に注目して考える. のうち,実部cos2π (ux + vy) を与える. この直線は なる. の直線は以下のように 1 2 cos ,... ,..., 2 , 1 , 0 = = + n n vy ux π

x

y

u / 1 v / 1 v / 2 v / 3 『空間周波数』と呼ば れる. を与える. は空間的な波の周波数 ⇒ ) , (u v 方向の周波数成分 方向の周波数成分 y v x u : : 「間隔が大きい」 が小さい」 「 となる. で 軸上に注目すると), (すなわち とおくと において, ⇔ =⇔ = = = = + u ux u u x ux x y n vy ux 1 ) cos( ,... / 2 , / 1 , 0 ... 2 , 1 , 0 0 ,... ,..., 2 , 1 , 0 u / 2

(5)

空間周波数の例

) ( 2 cos π ux + vy ,... 2 , , 0 ,... 2 , 1 , 0 0 / D y x D D x vy ux + = + = ⇔ = x y 例1) D 2 ) 0 , / 1 ( ) , (u v = D D x y D ) 0 , / 2 ( ) , (u v = D 例2) ,... 2 / 3 , , 2 / , 0 ,... 2 , 1 , 0 0 / 2x D y x D D D vy ux + = + = ⇔ = D 2

(6)

演習

) ( 2 cos π ux + vy x y u v D / 1 D A B 例題2 上図A,B,Cの位置に対応する空間周波 数のパターン(余弦波)をスケッチ しなさい. 例題1 下の図に対応する余弦関数を式で書き なさい.ただし黒い線は1の値をもち, 余弦関数の最大値を描いているものと する. 5 / D D / 1 D / 2 C

(7)

フーリエ変換演算のまとめ

One-comonent Image

x

y

u

v

x

y

0 1 2 3

0

1

2

3

u

v

x

y

x

y

x

y

∑∑

− = − =

+

=

1 0 1 0

}

/

)

(

2

exp{

)

,

(

1

)

,

(

N i N j

N

vy

ux

j

y

x

f

N

v

u

F

π

(8)

フーリエの合成のデモ

順次,高周波数成分を 追加していく. Manhattan distanceで Dm=3のスペクトル

u

v

u

v

F.T. F(u,v)

(9)

フーリエの合成のデモ(つづき)

Dm=3まで Dm=10まで Dm=6まで u v u v u v u v

(10)

2次元フーリエ変換

講義内容

講義内容

空間周波数の概念

2次元フーリエ変換

代表的な2次元フーリエ変換対

2次元離散フーリエ変換

(11)

代表的な2次元フーリエ変換対(1)

1 ) , ( ) , ( ) , (x y = x yF u v = f δ x u ) , ( ) , (x y x y f = δ 1 ) , (u v = F 0の関数. で無限大になり,他で 0 , 0 : ) , (x y x = y = δ 2変数のデルタ関数: 0の関数. で無限大になり,他で b y a x b y a x − , − ): = , = ( δ y v

(12)

代表的な2次元フーリエ変換対(2)

) ( sinc ) ( sinc ) , ( ) ( rect ) ( rect ) , (x y x y F u v u v f = ⇔ = x u y v x u 0 u 0 u − 0 0 v 0 vv y u / 1 2/u v / 1 v / 2 v / 3 v / 4 )} , ( ) , ( { 2 1 ) , ( )] ( 2 cos[ ) , (x y u0x v0y F u v u u0 v v0 u u0 v v0 f = π + ⇔ = δ − − +δ + −

(13)

代表的な2次元フーリエ変換対(3)

2 2 ) ( ) , ( r x y d r circ y x f = = + u v J1: ベッセル関数 x y d x u )] ( exp[ ] exp[ ) , ( 2 2 2 y x r y x f + − = − = π π y v 2 2 1 2 , ) ( ) , ( u v d d J d v u F = ρ = + ρ π ρ π π )] ( exp[ ] exp[ ) , ( 2 2 2 v u v u F + − = − = π πρ Gauss関数

(14)

2次元フーリエ変換の計算例

-矩形1-

) ( sinc ) ( sinc ) , ( ) ( rect ) ( rect ) , ( F u v au bv b y a x y x f = ⇔ = 6 , 12 = = b a

(15)

2次元フーリエ変換の計算例

-矩形1-

) ( sinc ) ( sinc ) , ( ) ( rect ) ( rect ) , ( F u v au bv b y a x y x f = ⇔ = 24 , 6 = = b a 24 , 6 = = b a 64 , 6 = = b a 64 , 6 = = b a

(16)

2次元フーリエ変換の計算例

-円形1-

2 2 ) ( ) , ( r x y d r circ y x f = = + ( , ) 2 1( ), u2 v2 d d J d v u F = ρ = + ρ π ρ π π

(17)

2次元フーリエ変換

講義内容

講義内容

空間周波数の概念

2次元フーリエ変換

代表的な2次元フーリエ変換対

2次元離散フーリエ変換

(18)

一般に,赤枠のように,原点が中央になるよう に配列し直して表示する方がわかりやすい. x y u u 2次元フーリエ変換 および振幅(絶対値) の対数変換表示 v v 2DFFTの結果は図のように原 点を端として切り出されたスペク トルと解釈できる.

2次元離散フーリエ変換

(19)

2次元離散フーリエ変換のデータの並び

N-1 0 N-1 0 N/2 Nyquist freq. N/2 u v

(20)

参照

関連したドキュメント

内 容 受講対象者 受講者数 研修月日

令和元年度予備費交付額 267億円 令和2年度第1次補正予算額 359億円 令和2年度第2次補正予算額 2,048億円 令和2年度第3次補正予算額 4,199億円 令和2年度予備費(

[r]

・条例第 37 条・第 62 条において、軽微なものなど規則で定める変更については、届出が不要とされ、その具 体的な要件が規則に定められている(規則第

区分 授業科目の名称 講義等の内容 備考.. 文 化

授業科目の名称 講義等の内容 備考

必修 幼二種 単位 ディプロマポリシーとの関連性

年平均値から日平均値(日平均値の2%除外値又は日平均値の年間98%値)への変換