Tweet
Tweet

Slide 1

Slide 1 text

関西大学総合情報学部 浅野　晃 画像情報処理 2024年度秋学期 第3部・ＣＴスキャナ — 投影からの画像の 再構成 ／ 第11回 逆投影法による再構成

Slide 2

Slide 2 text

20 2

Slide 3

Slide 3 text

20 2 投影切断面定理の復習

Slide 4

Slide 4 text

19 2024年度秋学期　画像情報処理 ／　関西大学総合情報学部　浅野　晃 投影切断面定理 3 投影群から２次元関数を再構成する fx F(fx, fy) fy θ ξ スライス F(ξcosθ, ξsinθ) ξ Gθ(ξ) 等しい x y θ s s g(s, θ) u 物体 投 影 物体の２次元 フーリエ変換 投影の１次元 フーリエ変換

Slide 5

Slide 5 text

19 2024年度秋学期　画像情報処理 ／　関西大学総合情報学部　浅野　晃 投影切断面定理 3 投影群から２次元関数を再構成する fx F(fx, fy) fy θ ξ スライス F(ξcosθ, ξsinθ) ξ Gθ(ξ) 等しい x y θ s s g(s, θ) u 物体 投 影 物体の２次元 フーリエ変換 投影の１次元 フーリエ変換

Slide 6

Slide 6 text

19 2024年度秋学期　画像情報処理 ／　関西大学総合情報学部　浅野　晃 投影切断面定理 3 投影群から２次元関数を再構成する fx F(fx, fy) fy θ ξ スライス F(ξcosθ, ξsinθ) ξ Gθ(ξ) 等しい x y θ s s g(s, θ) u 物体 投 影 物体の２次元 フーリエ変換 投影の１次元 フーリエ変換

Slide 7

Slide 7 text

19 2024年度秋学期　画像情報処理 ／　関西大学総合情報学部　浅野　晃 投影切断面定理 3 投影群から２次元関数を再構成する fx F(fx, fy) fy θ ξ スライス F(ξcosθ, ξsinθ) ξ Gθ(ξ) 等しい x y θ s s g(s, θ) u 物体 投 影 物体の２次元 フーリエ変換 投影の１次元 フーリエ変換

Slide 8

Slide 8 text

19 2024年度秋学期　画像情報処理 ／　関西大学総合情報学部　浅野　晃 投影切断面定理 3 投影群から２次元関数を再構成する fx F(fx, fy) fy θ ξ スライス F(ξcosθ, ξsinθ) ξ Gθ(ξ) 等しい x y θ s s g(s, θ) u 物体 投 影 物体の２次元 フーリエ変換 投影の１次元 フーリエ変換

Slide 9

Slide 9 text

19 2024年度秋学期　画像情報処理 ／　関西大学総合情報学部　浅野　晃 投影切断面定理 3 投影群から２次元関数を再構成する 「断面」がすべてそろえば，２次元逆フーリエ変換で２次元関数が再構成できる fx F(fx, fy) fy θ ξ スライス F(ξcosθ, ξsinθ) ξ Gθ(ξ) 等しい x y θ s s g(s, θ) u 物体 投 影 物体の２次元 フーリエ変換 投影の１次元 フーリエ変換

Slide 10

Slide 10 text

19 2024年度秋学期　画像情報処理 ／　関西大学総合情報学部　浅野　晃 フーリエ変換法による再構成の問題点 4 ２次元フーリエ変換の「すべての断面」を求めることはできない fx F(fx, fy) fy θ ξ スライス F(ξcosθ, ξsinθ) ξ Gθ(ξ) 等しい x y θ s s g(s, θ) u 物体 投 影 物体の２次元 フーリエ変換 投影の１次元 フーリエ変換

Slide 11

Slide 11 text

19 2024年度秋学期　画像情報処理 ／　関西大学総合情報学部　浅野　晃 フーリエ変換法による再構成の問題点 4 ２次元フーリエ変換の「すべての断面」を求めることはできない fx fy fx F(fx, fy) fy θ ξ スライス F(ξcosθ, ξsinθ) ξ Gθ(ξ) 等しい x y θ s s g(s, θ) u 物体 投 影 物体の２次元 フーリエ変換 投影の１次元 フーリエ変換

Slide 12

Slide 12 text

19 2024年度秋学期　画像情報処理 ／　関西大学総合情報学部　浅野　晃 フーリエ変換法による再構成の問題点 4 ２次元フーリエ変換の「すべての断面」を求めることはできない fx fy ひとつの投影＝ひとつの断面 fx F(fx, fy) fy θ ξ スライス F(ξcosθ, ξsinθ) ξ Gθ(ξ) 等しい x y θ s s g(s, θ) u 物体 投 影 物体の２次元 フーリエ変換 投影の１次元 フーリエ変換

Slide 13

Slide 13 text

19 2024年度秋学期　画像情報処理 ／　関西大学総合情報学部　浅野　晃 フーリエ変換法による再構成の問題点 4 ２次元フーリエ変換の「すべての断面」を求めることはできない fx fy ひとつの投影＝ひとつの断面 fx F(fx, fy) fy θ ξ スライス F(ξcosθ, ξsinθ) ξ Gθ(ξ) 等しい x y θ s s g(s, θ) u 物体 投 影 物体の２次元 フーリエ変換 投影の１次元 フーリエ変換

Slide 14

Slide 14 text

19 2024年度秋学期　画像情報処理 ／　関西大学総合情報学部　浅野　晃 フーリエ変換法による再構成の問題点 4 ２次元フーリエ変換の「すべての断面」を求めることはできない fx fy ひとつの投影＝ひとつの断面 fx F(fx, fy) fy θ ξ スライス F(ξcosθ, ξsinθ) ξ Gθ(ξ) 等しい x y θ s s g(s, θ) u 物体 投 影 物体の２次元 フーリエ変換 投影の１次元 フーリエ変換

Slide 15

Slide 15 text

19 2024年度秋学期　画像情報処理 ／　関西大学総合情報学部　浅野　晃 フーリエ変換法による再構成の問題点 4 ２次元フーリエ変換の「すべての断面」を求めることはできない fx fy ひとつの投影＝ひとつの断面 fx F(fx, fy) fy θ ξ スライス F(ξcosθ, ξsinθ) ξ Gθ(ξ) 等しい x y θ s s g(s, θ) u 物体 投 影 物体の２次元 フーリエ変換 投影の１次元 フーリエ変換

Slide 16

Slide 16 text

19 2024年度秋学期　画像情報処理 ／　関西大学総合情報学部　浅野　晃 フーリエ変換法による再構成の問題点 4 ２次元フーリエ変換の「すべての断面」を求めることはできない fx fy ひとつの投影＝ひとつの断面 fx F(fx, fy) fy θ ξ スライス F(ξcosθ, ξsinθ) ξ Gθ(ξ) 等しい x y θ s s g(s, θ) u 物体 投 影 物体の２次元 フーリエ変換 投影の１次元 フーリエ変換

Slide 17

Slide 17 text

19 2024年度秋学期　画像情報処理 ／　関西大学総合情報学部　浅野　晃 フーリエ変換法による再構成の問題点 4 ２次元フーリエ変換の「すべての断面」を求めることはできない fx fy ひとつの投影＝ひとつの断面 有限個の投影では，２次元フーリエ変換を埋め尽くすことはできない fx F(fx, fy) fy θ ξ スライス F(ξcosθ, ξsinθ) ξ Gθ(ξ) 等しい x y θ s s g(s, θ) u 物体 投 影 物体の２次元 フーリエ変換 投影の１次元 フーリエ変換

Slide 18

Slide 18 text

19 2024年度秋学期　画像情報処理 ／　関西大学総合情報学部　浅野　晃 フーリエ変換法による再構成の問題点 4 ２次元フーリエ変換の「すべての断面」を求めることはできない fx fy ひとつの投影＝ひとつの断面 有限個の投影では，２次元フーリエ変換を埋め尽くすことはできない →補間を行う fx F(fx, fy) fy θ ξ スライス F(ξcosθ, ξsinθ) ξ Gθ(ξ) 等しい x y θ s s g(s, θ) u 物体 投 影 物体の２次元 フーリエ変換 投影の１次元 フーリエ変換

Slide 19

Slide 19 text

19 2024年度秋学期　画像情報処理 ／　関西大学総合情報学部　浅野　晃 フーリエ変換法による再構成の問題点 5 fx fy 断面は極座標 周波数空間の誤差は，画像全体にひろがる アーティファクトを生む 補間を行う。が，コンピュータで計算する限りは「離散的」 ２次元フーリエ変換は正方座標 コンピュータの能力が低かった時代は 精密な計算が難しかった　→さてどうした？ fx F(fx, fy) fy θ ξ スライス F(ξcosθ, ξsinθ) ξ Gθ(ξ) 等しい x y θ s s g(s, θ) u 物体 投 影 物体の２次元 フーリエ変換 投影の１次元 フーリエ変換

Slide 20

Slide 20 text

20 6

Slide 21

Slide 21 text

20 6 逆投影法💡💡

Slide 22

Slide 22 text

19 2024年度秋学期　画像情報処理 ／　関西大学総合情報学部　浅野　晃 素朴な再構成 ー 逆投影 7 x y f(x, y)

Slide 23

Slide 23 text

19 2024年度秋学期　画像情報処理 ／　関西大学総合情報学部　浅野　晃 素朴な再構成 ー 逆投影 7 ２次元関数の任意の点 f(x, y) での値は x y f(x, y)

Slide 24

Slide 24 text

19 2024年度秋学期　画像情報処理 ／　関西大学総合情報学部　浅野　晃 素朴な再構成 ー 逆投影 7 ２次元関数の任意の点 f(x, y) での値は x y f(x, y)

Slide 25

Slide 25 text

19 2024年度秋学期　画像情報処理 ／　関西大学総合情報学部　浅野　晃 素朴な再構成 ー 逆投影 7 ２次元関数の任意の点 f(x, y) での値は f(x, y) を通るすべての投影（線積分）が わかれば求められる x y f(x, y)

Slide 26

Slide 26 text

19 2024年度秋学期　画像情報処理 ／　関西大学総合情報学部　浅野　晃 素朴な再構成 ー 逆投影 7 ２次元関数の任意の点 f(x, y) での値は f(x, y) を通るすべての投影（線積分）が わかれば求められる x y f(x, y)

Slide 27

Slide 27 text

19 2024年度秋学期　画像情報処理 ／　関西大学総合情報学部　浅野　晃 素朴な再構成 ー 逆投影 7 ２次元関数の任意の点 f(x, y) での値は f(x, y) を通るすべての投影（線積分）が わかれば求められる x y f(x, y) なら，それらの線積分をすべて合計してみれば？

Slide 28

Slide 28 text

19 2024年度秋学期　画像情報処理 ／　関西大学総合情報学部　浅野　晃 素朴な再構成 ー 逆投影 7 ２次元関数の任意の点 f(x, y) での値は f(x, y) を通るすべての投影（線積分）が わかれば求められる x y f(x, y) なら，それらの線積分をすべて合計してみれば？ どの線積分にも f(x, y) は含まれているのだから， 合計したら f(x, y) が強調される？

Slide 29

Slide 29 text

19 2024年度秋学期　画像情報処理 ／　関西大学総合情報学部　浅野　晃 逆投影法 8 x y θ s 軸s g(s, θ) u 物体 投 影 0 g(0, θ) s

Slide 30

Slide 30 text

19 2024年度秋学期　画像情報処理 ／　関西大学総合情報学部　浅野　晃 逆投影法 8 g(s,θ) x y θ s 軸s g(s, θ) u 物体 投 影 0 g(0, θ) s

Slide 31

Slide 31 text

19 2024年度秋学期　画像情報処理 ／　関西大学総合情報学部　浅野　晃 逆投影法 8 g(s,θ) x y θ s 軸s g(s, θ) u 物体 投 影 0 g(0, θ) s この線上では

Slide 32

Slide 32 text

19 2024年度秋学期　画像情報処理 ／　関西大学総合情報学部　浅野　晃 逆投影法 8 g(s,θ) x y θ s 軸s g(s, θ) u 物体 投 影 0 g(0, θ) s この線上では x cos θ + y sin θ − s = 0 　 　

Slide 33

Slide 33 text

19 2024年度秋学期　画像情報処理 ／　関西大学総合情報学部　浅野　晃 逆投影法 8 g(s,θ) x y θ s 軸s g(s, θ) u 物体 投 影 0 g(0, θ) s この線上では x cos θ + y sin θ − s = 0 　 　 s = x cos θ + y sin θ つまり

Slide 34

Slide 34 text

19 2024年度秋学期　画像情報処理 ／　関西大学総合情報学部　浅野　晃 逆投影法 8 g(s,θ) x y θ s 軸s g(s, θ) u 物体 投 影 0 g(0, θ) s この線上では x cos θ + y sin θ − s = 0 　 　 s = x cos θ + y sin θ つまり よって投影は

Slide 35

Slide 35 text

19 2024年度秋学期　画像情報処理 ／　関西大学総合情報学部　浅野　晃 逆投影法 8 g(s,θ) x y θ s 軸s g(s, θ) u 物体 投 影 0 g(0, θ) s この線上では x cos θ + y sin θ − s = 0 　 　 s = x cos θ + y sin θ つまり よって投影は g(x cos θ + y sin θ, θ)

Slide 36

Slide 36 text

19 2024年度秋学期　画像情報処理 ／　関西大学総合情報学部　浅野　晃 逆投影法 8 g(s,θ) x y θ s 軸s g(s, θ) u 物体 投 影 0 g(0, θ) s この線上では これを半周分足し合わせたのが逆投影 x cos θ + y sin θ − s = 0 　 　 s = x cos θ + y sin θ つまり よって投影は g(x cos θ + y sin θ, θ)

Slide 37

Slide 37 text

19 2024年度秋学期　画像情報処理 ／　関西大学総合情報学部　浅野　晃 逆投影法 8 g(s,θ) x y θ s 軸s g(s, θ) u 物体 投 影 0 g(0, θ) s この線上では これを半周分足し合わせたのが逆投影 x cos θ + y sin θ − s = 0 　 　 s = x cos θ + y sin θ つまり よって投影は g(x cos θ + y sin θ, θ) b(x, y) = π 0 g(x cos θ + y sin θ, θ)dθ

Slide 38

Slide 38 text

19 2024年度秋学期　画像情報処理 ／　関西大学総合情報学部　浅野　晃 逆投影法 8 g(s,θ) x y θ s 軸s g(s, θ) u 物体 投 影 0 g(0, θ) s この線上では これを半周分足し合わせたのが逆投影 x cos θ + y sin θ − s = 0 　 　 s = x cos θ + y sin θ つまり よって投影は g(x cos θ + y sin θ, θ) b(x, y) = π 0 g(x cos θ + y sin θ, θ)dθ

Slide 39

Slide 39 text

19 2024年度秋学期　画像情報処理 ／　関西大学総合情報学部　浅野　晃 逆投影法 8 g(s,θ) x y θ s 軸s g(s, θ) u 物体 投 影 0 g(0, θ) s この線上では これを半周分足し合わせたのが逆投影 x cos θ + y sin θ − s = 0 　 　 s = x cos θ + y sin θ つまり よって投影は g(x cos θ + y sin θ, θ) b(x, y) = π 0 g(x cos θ + y sin θ, θ)dθ 復元できているのか？ f(x, y) とどれほど違うのだろうか？

Slide 40

Slide 40 text

19 2024年度秋学期　画像情報処理 ／　関西大学総合情報学部　浅野　晃 復元できているのだろうか？ 9 b(x, y) = π 0 g(x cos θ + y sin θ, θ)dθ 逆投影

Slide 41

Slide 41 text

19 2024年度秋学期　画像情報処理 ／　関西大学総合情報学部　浅野　晃 復元できているのだろうか？ 9 b(x, y) = π 0 g(x cos θ + y sin θ, θ)dθ g(s, θ) = ∞ −∞ f(x, y)δ(x cos θ + y sin θ − s)dxdy Radon変換 逆投影

Slide 42

Slide 42 text

19 2024年度秋学期　画像情報処理 ／　関西大学総合情報学部　浅野　晃 復元できているのだろうか？ 9 b(x, y) = π 0 g(x cos θ + y sin θ, θ)dθ g(s, θ) = ∞ −∞ f(x, y)δ(x cos θ + y sin θ − s)dxdy Radon変換 逆投影

Slide 43

Slide 43 text

19 2024年度秋学期　画像情報処理 ／　関西大学総合情報学部　浅野　晃 復元できているのだろうか？ 9 s = x cos θ + y sin θ b(x, y) = π 0 g(x cos θ + y sin θ, θ)dθ g(s, θ) = ∞ −∞ f(x, y)δ(x cos θ + y sin θ − s)dxdy Radon変換 逆投影

Slide 44

Slide 44 text

19 2024年度秋学期　画像情報処理 ／　関西大学総合情報学部　浅野　晃 復元できているのだろうか？ 9 s = x cos θ + y sin θ b(x, y) = π 0 g(x cos θ + y sin θ, θ)dθ g(s, θ) = ∞ −∞ f(x, y)δ(x cos θ + y sin θ − s)dxdy Radon変換 逆投影

Slide 45

Slide 45 text

19 2024年度秋学期　画像情報処理 ／　関西大学総合情報学部　浅野　晃 復元できているのだろうか？ 9 s = x cos θ + y sin θ b(x, y) = π 0 g(x cos θ + y sin θ, θ)dθ g(s, θ) = ∞ −∞ f(x, y)δ(x cos θ + y sin θ − s)dxdy Radon変換 逆投影 b(x, y) = π 0 ∞ −∞ f(x , y )δ(x cos θ + y sin θ − (x cos θ + y sin θ))dx dy dθ = ∞ −∞ f(x , y ) π 0 δ((x − x) cos θ + (y − y) sin θ))dθ dx dy

Slide 46

Slide 46 text

19 2024年度秋学期　画像情報処理 ／　関西大学総合情報学部　浅野　晃 復元できているのだろうか？ 9 s = x cos θ + y sin θ b(x, y) = π 0 g(x cos θ + y sin θ, θ)dθ g(s, θ) = ∞ −∞ f(x, y)δ(x cos θ + y sin θ − s)dxdy Radon変換 逆投影 b(x, y) = π 0 ∞ −∞ f(x , y )δ(x cos θ + y sin θ − (x cos θ + y sin θ))dx dy dθ = ∞ −∞ f(x , y ) π 0 δ((x − x) cos θ + (y − y) sin θ))dθ dx dy

Slide 47

Slide 47 text

19 2024年度秋学期　画像情報処理 ／　関西大学総合情報学部　浅野　晃 復元できているのだろうか？ 9 s = x cos θ + y sin θ b(x, y) = π 0 g(x cos θ + y sin θ, θ)dθ g(s, θ) = ∞ −∞ f(x, y)δ(x cos θ + y sin θ − s)dxdy Radon変換 逆投影 b(x, y) = π 0 ∞ −∞ f(x , y )δ(x cos θ + y sin θ − (x cos θ + y sin θ))dx dy dθ = ∞ −∞ f(x , y ) π 0 δ((x − x) cos θ + (y − y) sin θ))dθ dx dy

Slide 48

Slide 48 text

19 2024年度秋学期　画像情報処理 ／　関西大学総合情報学部　浅野　晃 復元できているのだろうか？ 9 s = x cos θ + y sin θ b(x, y) = π 0 g(x cos θ + y sin θ, θ)dθ g(s, θ) = ∞ −∞ f(x, y)δ(x cos θ + y sin θ − s)dxdy Radon変換 逆投影 b(x, y) = π 0 ∞ −∞ f(x , y )δ(x cos θ + y sin θ − (x cos θ + y sin θ))dx dy dθ = ∞ −∞ f(x , y ) π 0 δ((x − x) cos θ + (y − y) sin θ))dθ dx dy

Slide 49

Slide 49 text

19 2024年度秋学期　画像情報処理 ／　関西大学総合情報学部　浅野　晃 復元できているのだろうか？ 9 s = x cos θ + y sin θ b(x, y) = π 0 g(x cos θ + y sin θ, θ)dθ g(s, θ) = ∞ −∞ f(x, y)δ(x cos θ + y sin θ − s)dxdy Radon変換 逆投影 b(x, y) = π 0 ∞ −∞ f(x , y )δ(x cos θ + y sin θ − (x cos θ + y sin θ))dx dy dθ = ∞ −∞ f(x , y ) π 0 δ((x − x) cos θ + (y − y) sin θ))dθ dx dy

Slide 50

Slide 50 text

19 2024年度秋学期　画像情報処理 ／　関西大学総合情報学部　浅野　晃 復元できているのだろうか？ 10 b(x, y) = π 0 ∞ −∞ f(x , y )δ(x cos θ + y sin θ − (x cos θ + y sin θ))dx dy dθ = ∞ −∞ f(x , y ) π 0 δ((x − x) cos θ + (y − y) sin θ))dθ dx dy

Slide 51

Slide 51 text

19 2024年度秋学期　画像情報処理 ／　関西大学総合情報学部　浅野　晃 復元できているのだろうか？ 10 b(x, y) = π 0 ∞ −∞ f(x , y )δ(x cos θ + y sin θ − (x cos θ + y sin θ))dx dy dθ = ∞ −∞ f(x , y ) π 0 δ((x − x) cos θ + (y − y) sin θ))dθ dx dy

Slide 52

Slide 52 text

19 2024年度秋学期　画像情報処理 ／　関西大学総合情報学部　浅野　晃 復元できているのだろうか？ 10 h(θ)が有限個のθk でしか0にならないとき b(x, y) = π 0 ∞ −∞ f(x , y )δ(x cos θ + y sin θ − (x cos θ + y sin θ))dx dy dθ = ∞ −∞ f(x , y ) π 0 δ((x − x) cos θ + (y − y) sin θ))dθ dx dy δ[h(θ)] = k 1 |h (θk)| δ[θ − θk]

Slide 53

Slide 53 text

19 2024年度秋学期　画像情報処理 ／　関西大学総合情報学部　浅野　晃 復元できているのだろうか？ 10 h(θ)が有限個のθk でしか0にならないとき b(x, y) = π 0 ∞ −∞ f(x , y )δ(x cos θ + y sin θ − (x cos θ + y sin θ))dx dy dθ = ∞ −∞ f(x , y ) π 0 δ((x − x) cos θ + (y − y) sin θ))dθ dx dy δ[h(θ)] = k 1 |h (θk)| δ[θ − θk] (x − x) cos θ + (y − y) sin θ = (x − x)2 + (y − y)2 sin(θ + α), α = cos−1 y −y √ (x −x)2+(y −y)2 = sin−1 x −x √ (x −x)2+(y −y)2 三角関数を合成

Slide 54

Slide 54 text

19 2024年度秋学期　画像情報処理 ／　関西大学総合情報学部　浅野　晃 復元できているのだろうか？ 10 h(θ)が有限個のθk でしか0にならないとき b(x, y) = π 0 ∞ −∞ f(x , y )δ(x cos θ + y sin θ − (x cos θ + y sin θ))dx dy dθ = ∞ −∞ f(x , y ) π 0 δ((x − x) cos θ + (y − y) sin θ))dθ dx dy δ[h(θ)] = k 1 |h (θk)| δ[θ − θk] (x − x) cos θ + (y − y) sin θ = (x − x)2 + (y − y)2 sin(θ + α), α = cos−1 y −y √ (x −x)2+(y −y)2 = sin−1 x −x √ (x −x)2+(y −y)2 三角関数を合成 θ = π − α のときだけ0

Slide 55

Slide 55 text

19 2024年度秋学期　画像情報処理 ／　関西大学総合情報学部　浅野　晃 復元できているのだろうか？ 10 h(θ)が有限個のθk でしか0にならないとき b(x, y) = π 0 ∞ −∞ f(x , y )δ(x cos θ + y sin θ − (x cos θ + y sin θ))dx dy dθ = ∞ −∞ f(x , y ) π 0 δ((x − x) cos θ + (y − y) sin θ))dθ dx dy δ[h(θ)] = k 1 |h (θk)| δ[θ − θk] (x − x) cos θ + (y − y) sin θ = (x − x)2 + (y − y)2 sin(θ + α), α = cos−1 y −y √ (x −x)2+(y −y)2 = sin−1 x −x √ (x −x)2+(y −y)2 δ((x − x) cos θ + (y − y) sin θ)) = 1 (x − x)2 + (y − y)2 cos(π) δ(θ − (π − α)) 三角関数を合成 θ = π − α のときだけ0

Slide 56

Slide 56 text

19 2024年度秋学期　画像情報処理 ／　関西大学総合情報学部　浅野　晃 復元できているといえなくもないが… 11 b(x, y) = π 0 ∞ −∞ f(x , y )δ(x cos θ + y sin θ − (x cos θ + y sin θ))dx dy dθ = ∞ −∞ f(x , y ) π 0 δ((x − x) cos θ + (y − y) sin θ))dθ dx dy

Slide 57

Slide 57 text

19 2024年度秋学期　画像情報処理 ／　関西大学総合情報学部　浅野　晃 復元できているといえなくもないが… 11 b(x, y) = π 0 ∞ −∞ f(x , y )δ(x cos θ + y sin θ − (x cos θ + y sin θ))dx dy dθ = ∞ −∞ f(x , y ) π 0 δ((x − x) cos θ + (y − y) sin θ))dθ dx dy

Slide 58

Slide 58 text

19 2024年度秋学期　画像情報処理 ／　関西大学総合情報学部　浅野　晃 復元できているといえなくもないが… 11 b(x, y) = π 0 ∞ −∞ f(x , y )δ(x cos θ + y sin θ − (x cos θ + y sin θ))dx dy dθ = ∞ −∞ f(x , y ) π 0 δ((x − x) cos θ + (y − y) sin θ))dθ dx dy δ((x − x) cos θ + (y − y) sin θ)) = 1 (x − x)2 + (y − y)2 cos(π) δ(θ − (π − α)) 　 　 　

Slide 59

Slide 59 text

19 2024年度秋学期　画像情報処理 ／　関西大学総合情報学部　浅野　晃 復元できているといえなくもないが… 11 b(x, y) = π 0 ∞ −∞ f(x , y )δ(x cos θ + y sin θ − (x cos θ + y sin θ))dx dy dθ = ∞ −∞ f(x , y ) π 0 δ((x − x) cos θ + (y − y) sin θ))dθ dx dy δ((x − x) cos θ + (y − y) sin θ)) = 1 (x − x)2 + (y − y)2 cos(π) δ(θ − (π − α)) 　 　 　

Slide 60

Slide 60 text

19 2024年度秋学期　画像情報処理 ／　関西大学総合情報学部　浅野　晃 復元できているといえなくもないが… 11 積分すると １ b(x, y) = π 0 ∞ −∞ f(x , y )δ(x cos θ + y sin θ − (x cos θ + y sin θ))dx dy dθ = ∞ −∞ f(x , y ) π 0 δ((x − x) cos θ + (y − y) sin θ))dθ dx dy δ((x − x) cos θ + (y − y) sin θ)) = 1 (x − x)2 + (y − y)2 cos(π) δ(θ − (π − α)) 　 　 　

Slide 61

Slide 61 text

19 2024年度秋学期　画像情報処理 ／　関西大学総合情報学部　浅野　晃 復元できているといえなくもないが… 11 積分すると １ b(x, y) = π 0 ∞ −∞ f(x , y )δ(x cos θ + y sin θ − (x cos θ + y sin θ))dx dy dθ = ∞ −∞ f(x , y ) π 0 δ((x − x) cos θ + (y − y) sin θ))dθ dx dy b(x, y) = ∞ −∞ f(x , y ) 1 (x − x)2 + (y − y)2 dx dy = f(x, y) ∗ 1 x2 + y2 δ((x − x) cos θ + (y − y) sin θ)) = 1 (x − x)2 + (y − y)2 cos(π) δ(θ − (π − α)) 　 　 　 よって

Slide 62

Slide 62 text

19 2024年度秋学期　画像情報処理 ／　関西大学総合情報学部　浅野　晃 復元できているといえなくもないが… 11 積分すると １ b(x, y) = π 0 ∞ −∞ f(x , y )δ(x cos θ + y sin θ − (x cos θ + y sin θ))dx dy dθ = ∞ −∞ f(x , y ) π 0 δ((x − x) cos θ + (y − y) sin θ))dθ dx dy b(x, y) = ∞ −∞ f(x , y ) 1 (x − x)2 + (y − y)2 dx dy = f(x, y) ∗ 1 x2 + y2 δ((x − x) cos θ + (y − y) sin θ)) = 1 (x − x)2 + (y − y)2 cos(π) δ(θ − (π − α)) 　 　 　 よって コンヴォリューションになっている

Slide 63

Slide 63 text

19 2024年度秋学期　画像情報処理 ／　関西大学総合情報学部　浅野　晃 復元するには 12 b(x, y)= f(x, y) ∗ 1 x2 + y2

Slide 64

Slide 64 text

19 2024年度秋学期　画像情報処理 ／　関西大学総合情報学部　浅野　晃 復元するには 12 フーリエ変換すると，コンヴォリューション→積 b(x, y)= f(x, y) ∗ 1 x2 + y2

Slide 65

Slide 65 text

19 2024年度秋学期　画像情報処理 ／　関西大学総合情報学部　浅野　晃 復元するには 12 フーリエ変換すると，コンヴォリューション→積 b(x, y)= f(x, y) ∗ 1 x2 + y2 FT[b(x, y)] = FT[f(x, y)] × FT 1 x2 + y2 ∴ FT[f(x, y)] = FT[b(x, y)]/FT 1 x2 + y2

Slide 66

Slide 66 text

19 2024年度秋学期　画像情報処理 ／　関西大学総合情報学部　浅野　晃 復元するには 12 フーリエ変換すると，コンヴォリューション→積 b(x, y)= f(x, y) ∗ 1 x2 + y2 FT[b(x, y)] = FT[f(x, y)] × FT 1 x2 + y2 ∴ FT[f(x, y)] = FT[b(x, y)]/FT 1 x2 + y2 FT 1 x2 + y2 = 1 f2 x + f2 y なので

Slide 67

Slide 67 text

19 2024年度秋学期　画像情報処理 ／　関西大学総合情報学部　浅野　晃 復元するには 12 フーリエ変換すると，コンヴォリューション→積 b(x, y)= f(x, y) ∗ 1 x2 + y2 FT[b(x, y)] = FT[f(x, y)] × FT 1 x2 + y2 ∴ FT[f(x, y)] = FT[b(x, y)]/FT 1 x2 + y2 FT[f(x, y)] = f2 x + f2 y × FT[b(x, y)] FT 1 x2 + y2 = 1 f2 x + f2 y なので

Slide 68

Slide 68 text

19 2024年度秋学期　画像情報処理 ／　関西大学総合情報学部　浅野　晃 これでいいのでしょうか？ 13 FT[f(x, y)] = f2 x + f2 y × FT[b(x, y)]

Slide 69

Slide 69 text

19 2024年度秋学期　画像情報処理 ／　関西大学総合情報学部　浅野　晃 これでいいのでしょうか？ 13 FT[f(x, y)] = f2 x + f2 y × FT[b(x, y)] fx = fy = 0 　 　

Slide 70

Slide 70 text

19 2024年度秋学期　画像情報処理 ／　関西大学総合情報学部　浅野　晃 これでいいのでしょうか？ 13 FT[f(x, y)] = f2 x + f2 y × FT[b(x, y)] fx = fy = 0 　 　 周波数0の成分は

Slide 71

Slide 71 text

19 2024年度秋学期　画像情報処理 ／　関西大学総合情報学部　浅野　晃 これでいいのでしょうか？ 13 FT[f(x, y)] = f2 x + f2 y × FT[b(x, y)] fx = fy = 0 　 　 FT[f(x, y)] = 0 　 　 周波数0の成分は

Slide 72

Slide 72 text

19 2024年度秋学期　画像情報処理 ／　関西大学総合情報学部　浅野　晃 これでいいのでしょうか？ 13 FT[f(x, y)] = f2 x + f2 y × FT[b(x, y)] fx = fy = 0 　 　 FT[f(x, y)] = 0 　 　 周波数0の成分は

Slide 73

Slide 73 text

19 2024年度秋学期　画像情報処理 ／　関西大学総合情報学部　浅野　晃 これでいいのでしょうか？ 13 FT[f(x, y)] = f2 x + f2 y × FT[b(x, y)] fx = fy = 0 　 　 FT[f(x, y)] = 0 　 　 周波数0の成分は 全体の平均が0

Slide 74

Slide 74 text

19 2024年度秋学期　画像情報処理 ／　関西大学総合情報学部　浅野　晃 これでいいのでしょうか？ 13 つまりf(x, y)は FT[f(x, y)] = f2 x + f2 y × FT[b(x, y)] fx = fy = 0 　 　 FT[f(x, y)] = 0 　 　 周波数0の成分は 全体の平均が0

Slide 75

Slide 75 text

19 2024年度秋学期　画像情報処理 ／　関西大学総合情報学部　浅野　晃 これでいいのでしょうか？ 13 つまりf(x, y)は FT[f(x, y)] = f2 x + f2 y × FT[b(x, y)] fx = fy = 0 　 　 FT[f(x, y)] = 0 　 　 周波数0の成分は 0 全体の平均が0

Slide 76

Slide 76 text

19 2024年度秋学期　画像情報処理 ／　関西大学総合情報学部　浅野　晃 これでいいのでしょうか？ 13 つまりf(x, y)は FT[f(x, y)] = f2 x + f2 y × FT[b(x, y)] fx = fy = 0 　 　 FT[f(x, y)] = 0 　 　 周波数0の成分は 0 全体の平均が0

Slide 77

Slide 77 text

19 2024年度秋学期　画像情報処理 ／　関西大学総合情報学部　浅野　晃 これでいいのでしょうか？ 13 つまりf(x, y)は FT[f(x, y)] = f2 x + f2 y × FT[b(x, y)] fx = fy = 0 　 　 FT[f(x, y)] = 0 　 　 周波数0の成分は 0 こうなっていることになる 全体の平均が0

Slide 78

Slide 78 text

19 2024年度秋学期　画像情報処理 ／　関西大学総合情報学部　浅野　晃 これでいいのでしょうか？ 13 つまりf(x, y)は FT[f(x, y)] = f2 x + f2 y × FT[b(x, y)] fx = fy = 0 　 　 FT[f(x, y)] = 0 　 　 周波数0の成分は 0 こうなっていることになる おかしい。 吸収率なんだから 全体の平均が0

Slide 79

Slide 79 text

19 2024年度秋学期　画像情報処理 ／　関西大学総合情報学部　浅野　晃 これでいいのでしょうか？ 13 つまりf(x, y)は FT[f(x, y)] = f2 x + f2 y × FT[b(x, y)] fx = fy = 0 　 　 FT[f(x, y)] = 0 　 　 周波数0の成分は 0 こうなっていることになる おかしい。 吸収率なんだから 0 全体の平均が0

Slide 80

Slide 80 text

19 2024年度秋学期　画像情報処理 ／　関西大学総合情報学部　浅野　晃 これでいいのでしょうか？ 13 つまりf(x, y)は FT[f(x, y)] = f2 x + f2 y × FT[b(x, y)] fx = fy = 0 　 　 FT[f(x, y)] = 0 　 　 周波数0の成分は 0 こうなっていることになる おかしい。 吸収率なんだから 0 全体の平均が0

Slide 81

Slide 81 text

19 2024年度秋学期　画像情報処理 ／　関西大学総合情報学部　浅野　晃 これでいいのでしょうか？ 13 つまりf(x, y)は FT[f(x, y)] = f2 x + f2 y × FT[b(x, y)] fx = fy = 0 　 　 FT[f(x, y)] = 0 　 　 周波数0の成分は 0 こうなっていることになる おかしい。 吸収率なんだから 0 値は正のはず 全体の平均が0

Slide 82

Slide 82 text

19 2024年度秋学期　画像情報処理 ／　関西大学総合情報学部　浅野　晃 これでいいのでしょうか？ 13 つまりf(x, y)は FT[f(x, y)] = f2 x + f2 y × FT[b(x, y)] fx = fy = 0 　 　 FT[f(x, y)] = 0 　 　 周波数0の成分は 0 こうなっていることになる おかしい。 吸収率なんだから 0 値は正のはず そもそも 全体の平均が0

Slide 83

Slide 83 text

19 2024年度秋学期　画像情報処理 ／　関西大学総合情報学部　浅野　晃 これでいいのでしょうか？ 13 つまりf(x, y)は FT[f(x, y)] = f2 x + f2 y × FT[b(x, y)] FT[b(x, y)] = FT[f(x, y)] × FT 1 x2 + y2 fx = fy = 0 　 　 FT[f(x, y)] = 0 　 　 周波数0の成分は 0 こうなっていることになる おかしい。 吸収率なんだから 0 値は正のはず そもそも 全体の平均が0

Slide 84

Slide 84 text

19 2024年度秋学期　画像情報処理 ／　関西大学総合情報学部　浅野　晃 これでいいのでしょうか？ 13 つまりf(x, y)は FT[f(x, y)] = f2 x + f2 y × FT[b(x, y)] FT[b(x, y)] = FT[f(x, y)] × FT 1 x2 + y2 FT 1 x2 + y2 = 1 f2 x + f2 y fx = fy = 0 　 　 FT[f(x, y)] = 0 　 　 周波数0の成分は 0 こうなっていることになる おかしい。 吸収率なんだから 0 値は正のはず そもそも 全体の平均が0

Slide 85

Slide 85 text

19 2024年度秋学期　画像情報処理 ／　関西大学総合情報学部　浅野　晃 これでいいのでしょうか？ 13 つまりf(x, y)は FT[f(x, y)] = f2 x + f2 y × FT[b(x, y)] FT[b(x, y)] = FT[f(x, y)] × FT 1 x2 + y2 FT 1 x2 + y2 = 1 f2 x + f2 y fx = fy = 0 　 　 FT[f(x, y)] = 0 　 　 周波数0の成分は 0 こうなっていることになる おかしい。 吸収率なんだから 0 値は正のはず そもそも fx = fy = 0 で発散している 全体の平均が0

Slide 86

Slide 86 text

19 2024年度秋学期　画像情報処理 ／　関西大学総合情報学部　浅野　晃 フィルタ補正逆投影法 14 再び投影切断面定理 fx F(fx, fy) fy θ ξ 断面 F(ξcosθ, ξsinθ) ξ Gθ(ξ) 等しい x y θ s s g(s, θ) u 物体 投 影 物体の２次元 フーリエ変換 投影の１次元 フーリエ変換

Slide 87

Slide 87 text

19 2024年度秋学期　画像情報処理 ／　関西大学総合情報学部　浅野　晃 フィルタ補正逆投影法 14 再び投影切断面定理 fx F(fx, fy) fy θ ξ 断面 F(ξcosθ, ξsinθ) ξ Gθ(ξ) 等しい x y θ s s g(s, θ) u 物体 投 影 物体の２次元 フーリエ変換 投影の１次元 フーリエ変換 復元

Slide 88

Slide 88 text

19 2024年度秋学期　画像情報処理 ／　関西大学総合情報学部　浅野　晃 フィルタ補正逆投影法 14 再び投影切断面定理 fx F(fx, fy) fy θ ξ 断面 F(ξcosθ, ξsinθ) ξ Gθ(ξ) 等しい x y θ s s g(s, θ) u 物体 投 影 物体の２次元 フーリエ変換 投影の１次元 フーリエ変換 復元 ２次元逆フーリエ変換

Slide 89

Slide 89 text

19 2024年度秋学期　画像情報処理 ／　関西大学総合情報学部　浅野　晃 フィルタ補正逆投影法 14 再び投影切断面定理 fx F(fx, fy) fy θ ξ 断面 F(ξcosθ, ξsinθ) ξ Gθ(ξ) 等しい x y θ s s g(s, θ) u 物体 投 影 物体の２次元 フーリエ変換 投影の１次元 フーリエ変換 f(x, y) = ∞ −∞ F(fx, fy) exp(i2π(fxx + fyy))dfxdfy 復元 ２次元逆フーリエ変換

Slide 90

Slide 90 text

19 2024年度秋学期　画像情報処理 ／　関西大学総合情報学部　浅野　晃 フィルタ補正逆投影法 15 f(x, y) = ∞ −∞ F(fx, fy) exp(i2π(fxx + fyy))dfxdfy

Slide 91

Slide 91 text

19 2024年度秋学期　画像情報処理 ／　関西大学総合情報学部　浅野　晃 フィルタ補正逆投影法 15 極座標に変換 f(x, y) = ∞ −∞ F(fx, fy) exp(i2π(fxx + fyy))dfxdfy

Slide 92

Slide 92 text

19 2024年度秋学期　画像情報処理 ／　関西大学総合情報学部　浅野　晃 フィルタ補正逆投影法 15 極座標に変換 f(x, y) = ∞ −∞ F(fx, fy) exp(i2π(fxx + fyy))dfxdfy fx = ξ cos θ, fy = ξ sin θ

Slide 93

Slide 93 text

19 2024年度秋学期　画像情報処理 ／　関西大学総合情報学部　浅野　晃 フィルタ補正逆投影法 15 極座標に変換 f(x, y) = ∞ −∞ F(fx, fy) exp(i2π(fxx + fyy))dfxdfy fx = ξ cos θ, fy = ξ sin θ

Slide 94

Slide 94 text

19 2024年度秋学期　画像情報処理 ／　関西大学総合情報学部　浅野　晃 フィルタ補正逆投影法 15 極座標に変換 f(x, y) = ∞ −∞ F(fx, fy) exp(i2π(fxx + fyy))dfxdfy fx = ξ cos θ, fy = ξ sin θ f(x, y) = 2π 0 ∞ 0 F(ξ cos θ, ξ sin θ) exp(i2πξ(x cos θ + y sin θ))ξdξdθ

Slide 95

Slide 95 text

19 2024年度秋学期　画像情報処理 ／　関西大学総合情報学部　浅野　晃 フィルタ補正逆投影法 15 極座標に変換 f(x, y) = ∞ −∞ F(fx, fy) exp(i2π(fxx + fyy))dfxdfy fx = ξ cos θ, fy = ξ sin θ f(x, y) = 2π 0 ∞ 0 F(ξ cos θ, ξ sin θ) exp(i2πξ(x cos θ + y sin θ))ξdξdθ

Slide 96

Slide 96 text

19 2024年度秋学期　画像情報処理 ／　関西大学総合情報学部　浅野　晃 フィルタ補正逆投影法 15 極座標に変換 f(x, y) = ∞ −∞ F(fx, fy) exp(i2π(fxx + fyy))dfxdfy fx = ξ cos θ, fy = ξ sin θ f(x, y) = 2π 0 ∞ 0 F(ξ cos θ, ξ sin θ) exp(i2πξ(x cos θ + y sin θ))ξdξdθ dfx dfy 正方座標の微小面積

Slide 97

Slide 97 text

19 2024年度秋学期　画像情報処理 ／　関西大学総合情報学部　浅野　晃 フィルタ補正逆投影法 15 極座標に変換 f(x, y) = ∞ −∞ F(fx, fy) exp(i2π(fxx + fyy))dfxdfy fx = ξ cos θ, fy = ξ sin θ f(x, y) = 2π 0 ∞ 0 F(ξ cos θ, ξ sin θ) exp(i2πξ(x cos θ + y sin θ))ξdξdθ dfx dfy ξ dξ ξdθ dθ 正方座標の微小面積 極座標の微小面積

Slide 98

Slide 98 text

19 2024年度秋学期　画像情報処理 ／　関西大学総合情報学部　浅野　晃 フィルタ補正逆投影法 15 極座標に変換 f(x, y) = ∞ −∞ F(fx, fy) exp(i2π(fxx + fyy))dfxdfy fx = ξ cos θ, fy = ξ sin θ f(x, y) = 2π 0 ∞ 0 F(ξ cos θ, ξ sin θ) exp(i2πξ(x cos θ + y sin θ))ξdξdθ dfx dfy ξ dξ ξdθ dθ 正方座標の微小面積 極座標の微小面積

Slide 99

Slide 99 text

19 2024年度秋学期　画像情報処理 ／　関西大学総合情報学部　浅野　晃 フィルタ補正逆投影法 15 極座標に変換 f(x, y) = ∞ −∞ F(fx, fy) exp(i2π(fxx + fyy))dfxdfy fx = ξ cos θ, fy = ξ sin θ f(x, y) = 2π 0 ∞ 0 F(ξ cos θ, ξ sin θ) exp(i2πξ(x cos θ + y sin θ))ξdξdθ Gθ(ξ) = F(ξ cos θ, ξ sin θ) 投影切断面定理より dfx dfy ξ dξ ξdθ dθ 正方座標の微小面積 極座標の微小面積

Slide 100

Slide 100 text

19 2024年度秋学期　画像情報処理 ／　関西大学総合情報学部　浅野　晃 フィルタ補正逆投影法 15 極座標に変換 f(x, y) = ∞ −∞ F(fx, fy) exp(i2π(fxx + fyy))dfxdfy fx = ξ cos θ, fy = ξ sin θ f(x, y) = 2π 0 ∞ 0 F(ξ cos θ, ξ sin θ) exp(i2πξ(x cos θ + y sin θ))ξdξdθ Gθ(ξ) = F(ξ cos θ, ξ sin θ) f(x, y) = 2π 0 ∞ 0 Gθ(ξ) exp(i2πξ(x cos θ + y sin θ))ξdξdθ 投影切断面定理より dfx dfy ξ dξ ξdθ dθ 正方座標の微小面積 極座標の微小面積

Slide 101

Slide 101 text

19 2024年度秋学期　画像情報処理 ／　関西大学総合情報学部　浅野　晃 フィルタ補正逆投影法 16 f(x, y) = 2π 0 ∞ 0 Gθ(ξ) exp(i2πξ(x cos θ + y sin θ))ξdξdθ 　 　 　

Slide 102

Slide 102 text

19 2024年度秋学期　画像情報処理 ／　関西大学総合情報学部　浅野　晃 フィルタ補正逆投影法 16 積分区間を変換 f(x, y) = 2π 0 ∞ 0 Gθ(ξ) exp(i2πξ(x cos θ + y sin θ))ξdξdθ 　 　 　

Slide 103

Slide 103 text

19 2024年度秋学期　画像情報処理 ／　関西大学総合情報学部　浅野　晃 フィルタ補正逆投影法 16 積分区間を変換 f(x, y) = 2π 0 ∞ 0 Gθ(ξ) exp(i2πξ(x cos θ + y sin θ))ξdξdθ 　 　 　 f(x, y) = π 0 ∞ −∞ Gθ(ξ) exp(i2πξ(x cos θ + y sin θ))|ξ|dξdθ 　 　 　

Slide 104

Slide 104 text

19 2024年度秋学期　画像情報処理 ／　関西大学総合情報学部　浅野　晃 フィルタ補正逆投影法 16 積分区間を変換 f(x, y) = 2π 0 ∞ 0 Gθ(ξ) exp(i2πξ(x cos θ + y sin θ))ξdξdθ 　 　 　 f(x, y) = π 0 ∞ −∞ Gθ(ξ) exp(i2πξ(x cos θ + y sin θ))|ξ|dξdθ 　 　 　

Slide 105

Slide 105 text

19 2024年度秋学期　画像情報処理 ／　関西大学総合情報学部　浅野　晃 フィルタ補正逆投影法 16 積分区間を変換 f(x, y) = 2π 0 ∞ 0 Gθ(ξ) exp(i2πξ(x cos θ + y sin θ))ξdξdθ 　 　 　 f(x, y) = π 0 ∞ −∞ Gθ(ξ) exp(i2πξ(x cos θ + y sin θ))|ξ|dξdθ 　 　 　

Slide 106

Slide 106 text

19 2024年度秋学期　画像情報処理 ／　関西大学総合情報学部　浅野　晃 フィルタ補正逆投影法 16 積分区間を変換 f(x, y) = 2π 0 ∞ 0 Gθ(ξ) exp(i2πξ(x cos θ + y sin θ))ξdξdθ 　 　 　 f(x, y) = π 0 ∞ −∞ Gθ(ξ) exp(i2πξ(x cos θ + y sin θ))|ξ|dξdθ 　 　 　 s = x cos θ + y sin θ

Slide 107

Slide 107 text

19 2024年度秋学期　画像情報処理 ／　関西大学総合情報学部　浅野　晃 フィルタ補正逆投影法 16 積分区間を変換 f(x, y) = 2π 0 ∞ 0 Gθ(ξ) exp(i2πξ(x cos θ + y sin θ))ξdξdθ 　 　 　 f(x, y) = π 0 ∞ −∞ Gθ(ξ) exp(i2πξ(x cos θ + y sin θ))|ξ|dξdθ 　 　 　 f(x, y) = π 0 ∞ −∞ |ξ|Gθ(ξ) exp(i2πsξ)dξ dθ s = x cos θ + y sin θ

Slide 108

Slide 108 text

19 2024年度秋学期　画像情報処理 ／　関西大学総合情報学部　浅野　晃 フィルタ補正逆投影法 16 積分区間を変換 f(x, y) = 2π 0 ∞ 0 Gθ(ξ) exp(i2πξ(x cos θ + y sin θ))ξdξdθ 　 　 　 f(x, y) = π 0 ∞ −∞ Gθ(ξ) exp(i2πξ(x cos θ + y sin θ))|ξ|dξdθ 　 　 　 f(x, y) = π 0 ∞ −∞ |ξ|Gθ(ξ) exp(i2πsξ)dξ dθ s = x cos θ + y sin θ

Slide 109

Slide 109 text

19 2024年度秋学期　画像情報処理 ／　関西大学総合情報学部　浅野　晃 フィルタ補正逆投影法 16 積分区間を変換 f(x, y) = 2π 0 ∞ 0 Gθ(ξ) exp(i2πξ(x cos θ + y sin θ))ξdξdθ 　 　 　 f(x, y) = π 0 ∞ −∞ Gθ(ξ) exp(i2πξ(x cos θ + y sin θ))|ξ|dξdθ 　 　 　 f(x, y) = π 0 ∞ −∞ |ξ|Gθ(ξ) exp(i2πsξ)dξ dθ s = x cos θ + y sin θ

Slide 110

Slide 110 text

19 2024年度秋学期　画像情報処理 ／　関西大学総合情報学部　浅野　晃 フィルタ補正逆投影法 16 積分区間を変換 f(x, y) = 2π 0 ∞ 0 Gθ(ξ) exp(i2πξ(x cos θ + y sin θ))ξdξdθ 　 　 　 f(x, y) = π 0 ∞ −∞ Gθ(ξ) exp(i2πξ(x cos θ + y sin θ))|ξ|dξdθ 　 　 　 f(x, y) = π 0 ∞ −∞ |ξ|Gθ(ξ) exp(i2πsξ)dξ dθ s = x cos θ + y sin θ の逆フーリエ変換 |ξ|Gθ(ξ)

Slide 111

Slide 111 text

19 2024年度秋学期　画像情報処理 ／　関西大学総合情報学部　浅野　晃 フィルタ補正逆投影法 17 f(x, y) = π 0 ∞ −∞ |ξ|Gθ(ξ) exp(i2πsξ)dξ dθ の逆フーリエ変換 |ξ|Gθ(ξ)

Slide 112

Slide 112 text

19 2024年度秋学期　画像情報処理 ／　関西大学総合情報学部　浅野　晃 フィルタ補正逆投影法 17 f(x, y) = π 0 ∞ −∞ |ξ|Gθ(ξ) exp(i2πsξ)dξ dθ の逆フーリエ変換 |ξ|Gθ(ξ) ˆ g(s, θ) ≡ ˆ g(x cos θ + y sin θ, θ) ≡ ∞ −∞ |ξ|Gθ(ξ) exp(i2πsξ)dξ

Slide 113

Slide 113 text

19 2024年度秋学期　画像情報処理 ／　関西大学総合情報学部　浅野　晃 フィルタ補正逆投影法 17 f(x, y) = π 0 ∞ −∞ |ξ|Gθ(ξ) exp(i2πsξ)dξ dθ の逆フーリエ変換 |ξ|Gθ(ξ) ˆ g(s, θ) ≡ ˆ g(x cos θ + y sin θ, θ) ≡ ∞ −∞ |ξ|Gθ(ξ) exp(i2πsξ)dξ とおくと，

Slide 114

Slide 114 text

19 2024年度秋学期　画像情報処理 ／　関西大学総合情報学部　浅野　晃 フィルタ補正逆投影法 17 f(x, y) = π 0 ∞ −∞ |ξ|Gθ(ξ) exp(i2πsξ)dξ dθ の逆フーリエ変換 |ξ|Gθ(ξ) ˆ g(s, θ) ≡ ˆ g(x cos θ + y sin θ, θ) ≡ ∞ −∞ |ξ|Gθ(ξ) exp(i2πsξ)dξ f(x, y) = π 0 ˆ g(x cos θ + y sin θ, θ)dθ とおくと，

Slide 115

Slide 115 text

19 2024年度秋学期　画像情報処理 ／　関西大学総合情報学部　浅野　晃 フィルタ補正逆投影法 17 f(x, y) = π 0 ∞ −∞ |ξ|Gθ(ξ) exp(i2πsξ)dξ dθ の逆フーリエ変換 |ξ|Gθ(ξ) ˆ g(s, θ) ≡ ˆ g(x cos θ + y sin θ, θ) ≡ ∞ −∞ |ξ|Gθ(ξ) exp(i2πsξ)dξ f(x, y) = π 0 ˆ g(x cos θ + y sin θ, θ)dθ とおくと， 投影を，「周波数空間で |ξ| 倍するフィルタ」を適用してから，逆投影

Slide 116

Slide 116 text

19 2024年度秋学期　画像情報処理 ／　関西大学総合情報学部　浅野　晃 コンヴォリューション逆投影法 18 ˆ g(s, θ) ≡ ˆ g(x cos θ + y sin θ, θ) ≡ ∞ −∞ |ξ|Gθ(ξ) exp(i2πsξ)dξ

Slide 117

Slide 117 text

19 2024年度秋学期　画像情報処理 ／　関西大学総合情報学部　浅野　晃 コンヴォリューション逆投影法 18 ˆ g(s, θ) ≡ ˆ g(x cos θ + y sin θ, θ) ≡ ∞ −∞ |ξ|Gθ(ξ) exp(i2πsξ)dξ

Slide 118

Slide 118 text

19 2024年度秋学期　画像情報処理 ／　関西大学総合情報学部　浅野　晃 コンヴォリューション逆投影法 18 ˆ g(s, θ) ≡ ˆ g(x cos θ + y sin θ, θ) ≡ ∞ −∞ |ξ|Gθ(ξ) exp(i2πsξ)dξ ある角度 θ での投影を 周波数空間で |ξ| 倍

Slide 119

Slide 119 text

19 2024年度秋学期　画像情報処理 ／　関西大学総合情報学部　浅野　晃 コンヴォリューション逆投影法 18 ˆ g(s, θ) ≡ ˆ g(x cos θ + y sin θ, θ) ≡ ∞ −∞ |ξ|Gθ(ξ) exp(i2πsξ)dξ ある角度 θ での投影を 周波数空間で |ξ| 倍 実空間では FT−1[|ξ|] とのコンヴォリューション

Slide 120

Slide 120 text

19 2024年度秋学期　画像情報処理 ／　関西大学総合情報学部　浅野　晃 コンヴォリューション逆投影法 18 ˆ g(s, θ) ≡ ˆ g(x cos θ + y sin θ, θ) ≡ ∞ −∞ |ξ|Gθ(ξ) exp(i2πsξ)dξ ˆ g(s, θ) = g(s, θ) ∗ FT−1[|ξ|] ある角度 θ での投影を 周波数空間で |ξ| 倍 実空間では FT−1[|ξ|] とのコンヴォリューション

Slide 121

Slide 121 text

19 2024年度秋学期　画像情報処理 ／　関西大学総合情報学部　浅野　晃 最初からこれを定義しておけば，フーリエ変換の必要はない コンヴォリューション逆投影法 18 ˆ g(s, θ) ≡ ˆ g(x cos θ + y sin θ, θ) ≡ ∞ −∞ |ξ|Gθ(ξ) exp(i2πsξ)dξ ˆ g(s, θ) = g(s, θ) ∗ FT−1[|ξ|] ある角度 θ での投影を 周波数空間で |ξ| 倍 実空間では FT−1[|ξ|] とのコンヴォリューション

Slide 122

Slide 122 text

19 2024年度秋学期　画像情報処理 ／　関西大学総合情報学部　浅野　晃 最初からこれを定義しておけば，フーリエ変換の必要はない コンヴォリューション逆投影法 18 ˆ g(s, θ) ≡ ˆ g(x cos θ + y sin θ, θ) ≡ ∞ −∞ |ξ|Gθ(ξ) exp(i2πsξ)dξ ˆ g(s, θ) = g(s, θ) ∗ FT−1[|ξ|] ある角度 θ での投影を 周波数空間で |ξ| 倍 実空間では FT−1[|ξ|] とのコンヴォリューション 極座標→直交座標の変換は実空間で行うので， 誤差が画面全体に拡散することはない

Slide 123

Slide 123 text

19 2024年度秋学期　画像情報処理 ／　関西大学総合情報学部　浅野　晃 フィルタ関数 19 ξ H(ξ) ξmax – ξmax ξ H(ξ) ξ H(ξ) ξmax – ξmax (a) (b) (c)

Slide 124

Slide 124 text

19 2024年度秋学期　画像情報処理 ／　関西大学総合情報学部　浅野　晃 フィルタ関数 19 ξ H(ξ) ξmax – ξmax ξ H(ξ) ξ H(ξ) ξmax – ξmax (a) (b) (c) |ξ|

Slide 125

Slide 125 text

19 2024年度秋学期　画像情報処理 ／　関西大学総合情報学部　浅野　晃 フィルタ関数 19 ξ H(ξ) ξmax – ξmax ξ H(ξ) ξ H(ξ) ξmax – ξmax (a) (b) (c) |ξ| H(ξ) = |ξ|rect ξ 2ξmax

Slide 126

Slide 126 text

19 2024年度秋学期　画像情報処理 ／　関西大学総合情報学部　浅野　晃 フィルタ関数 19 ξ H(ξ) ξmax – ξmax ξ H(ξ) ξ H(ξ) ξmax – ξmax (a) (b) (c) |ξ| H(ξ) = |ξ|rect ξ 2ξmax H(ξ) = |ξ|sinc ξ 2ξmax rect ξ 2ξmax

Slide 127

Slide 127 text

19 2024年度秋学期　画像情報処理 ／　関西大学総合情報学部　浅野　晃 フィルタ関数 19 ξ H(ξ) ξmax – ξmax ξ H(ξ) ξ H(ξ) ξmax – ξmax (a) (b) (c) |ξ| H(ξ) = |ξ|rect ξ 2ξmax H(ξ) = |ξ|sinc ξ 2ξmax rect ξ 2ξmax 高い周波数成分を増幅すると ノイズを強調してしまうので，抑える

Slide 128

Slide 128 text

19 2024年度秋学期　画像情報処理 ／　関西大学総合情報学部　浅野　晃 フィルタ関数 19 ξ H(ξ) ξmax – ξmax ξ H(ξ) ξ H(ξ) ξmax – ξmax (a) (b) (c) |ξ| H(ξ) = |ξ|rect ξ 2ξmax H(ξ) = |ξ|sinc ξ 2ξmax rect ξ 2ξmax 高い周波数成分を増幅すると ノイズを強調してしまうので，抑える ※現代のCTスキャナでは， 　初期状態の物体から計算で投影を求める→実際の投影と比較して，物体を修正する 　という操作を繰り返すことで，実際の物体に近づけていく，という方法（逐次近似法）も 　用いられています