1 / 24

６．４　離散的コサイン変換 （ DCT : discrete cosine transform ） （１）ＤＣＴとは

６．４　離散的コサイン変換 （ DCT : discrete cosine transform ） （１）ＤＣＴとは. ＤＦＴの式. が偶関数であれば，ＤＦＴの式の虚数項を０として. ＤＣＴの解釈. 半周期分 が与えられたものとして， 偶関数 になるように拡張し （Ｓｉｎ成分，すなわち虚数項＝０） 周期的に接続 したＤＦＴとみなすことができる。. （２）偶関数への拡張の方法. ① 軸のところを中心とみなす方法 ②軸から０．５ずれたところで折り返す方法. ① ｎ = ０ を中心として偶関数に拡張. 元の信号.

chun
Télécharger la présentation

６．４　離散的コサイン変換 （ DCT : discrete cosine transform ） （１）ＤＣＴとは

An Image/Link below is provided (as is) to download presentation Download Policy: Content on the Website is provided to you AS IS for your information and personal use and may not be sold / licensed / shared on other websites without getting consent from its author. Content is provided to you AS IS for your information and personal use only. Download presentation by click this link. While downloading, if for some reason you are not able to download a presentation, the publisher may have deleted the file from their server. During download, if you can't get a presentation, the file might be deleted by the publisher.

E N D

Presentation Transcript

  1. ６．４　離散的コサイン変換（DCT : discrete cosine transform ）（１）ＤＣＴとは ＤＦＴの式 が偶関数であれば，ＤＦＴの式の虚数項を０として

  2. ＤＣＴの解釈 半周期分が与えられたものとして， 偶関数になるように拡張し （Ｓｉｎ成分，すなわち虚数項＝０） 周期的に接続したＤＦＴとみなすことができる。

  3. （２）偶関数への拡張の方法 ①軸のところを中心とみなす方法 ②軸から０．５ずれたところで折り返す方法 ① ｎ = ０ を中心として偶関数に拡張 元の信号 ① ｎ = －０．５ を中心として偶関数に拡張

  4. 軸から０．５ずらす場合の接続方法も２種類 ①周期２Ｎ－１として接続 ②周期２Ｎとして接続（第Ⅱ種ＤＣＴで採用） ①周期２Ｎ－１として周期的に接続 ②周期２Ｎとして周期的に接続

  5. （３）ＤＣＴの種類 偶関数への拡張方法，周期的な接続方法の違いにより，以下のような種類がある。

  6. （４）第Ⅱ種ＤＣＴ 第Ⅱ種は，ＪＰＥＧ，ＭＰＥＧなど 画像処理で最もよく使われる方法 （逆変換）第Ⅲ種と同じ

  7. VBAでの記述　①データ宣言とパルス信号の設定（色々な信号を設定してみよう）VBAでの記述　①データ宣言とパルス信号の設定（色々な信号を設定してみよう） Private Const Num = 128 Private Const PI = 6.28318530717959 / 2 Private X(Num) As Double Private Y(Num) As Double Private Z(Num) As Double Private Sub パルス() For K = 0 To Num - 1 X(K) = -1 If K >= 30 And K < 80 Then X(K) = 1 Next X(Num) = X(0) End Sub

  8. VBAでの記述　②結果設定とボタンＣｌｉｃｋイベントハンドラVBAでの記述　②結果設定とボタンＣｌｉｃｋイベントハンドラ Private Sub 結果設定() With Worksheets("Sheet1") .Cells(1, 1) = "No.": .Cells(1, 2) = "X" .Cells(1, 3) = "Y": .Cells(1, 4) = "Z" For K = 0 To Num .Cells(K + 2, 1) = K: .Cells(K + 2, 2) = X(K) .Cells(K + 2, 3) = Y(K): .Cells(K + 2, 4) = Z(K) Next End With End Sub Sub ボタン1_Click() パルス DCT IDCT 結果設定 End Sub

  9. VBAでの記述　③ＤＣＴ Private Sub DCT() 'DCT(離散的コサイン変換) Dim C0 As Double: Dim Cn0 As Double C0 = 1 / Sqr(Num): Cn0 = Sqr(2#) * C0 For K = 0 To Num - 1 Y(K) = 0 For j = 0 To Num - 1 Y(K) = Y(K) + X(j) * Cos((2 * j + 1) * K * PI / (2 * Num)) Next If K = 0 Then Y(K) = C0 * Y(K) Else Y(K) = Cn0 * Y(K) End If Next Y(Num) = 0 End Sub

  10. VBAでの記述　④ＩＤＣＴ Private Sub IDCT() 'IDCT(離散的逆変換) Dim C0 As Double: Dim C As Double C0 = 1# / Sqr(2#): C = Sqr(2# / Num) For j = 0 To Num - 1 Z(j) = C0 * Y(0) For K = 1 To Num - 1 Z(j) = Z(j) + Y(K) * Cos((2 * j + 1) * K * PI / (2 * Num)) Next Z(j) = C * Z(j) Next Z(Num) = Z(0) End Sub

  11. ＤＣＴの結果

  12. 色々なＤＣＴの結果（１）

  13. 色々なＤＣＴの結果（２）

  14. （５）高速ＤＣＴ（Ｃｈｅｎのアルゴリズム）（５）高速ＤＣＴ（Ｃｈｅｎのアルゴリズム） ＦＦＴを使ってＤＣＴを効率よく実行することも可能。 ここでは，Ｃｈｅｎのアルゴリズムを示す。 ＤＣＴはＮ＝４のとき次のように行列式で表現できる。 これを小行列に分解し，因数分解することで計算量を減らす。

  15. 計算の手順 ＪＰＥＧ，ＭＰＥＧで使われているＮ＝８の場合

  16. VBAでの記述　①データ宣言とパルス信号の設定VBAでの記述　①データ宣言とパルス信号の設定 Private Const Num = 8 Private Const PI = 3.12159265358979 Private X(Num) As Double '入力データ Private Y(Num) As Double '変換結果 Private Z(Num) As Double '逆変換結果 Private Sub パルス() For K = 0 To Num - 1 X(K) = -1 If K >= 3 And K <= 5 Then X(K) = 1 Next X(Num) = X(0) End Sub

  17. VBAでの記述　②結果設定とボタンのＣｌｉｃｋイベントハンドラVBAでの記述　②結果設定とボタンのＣｌｉｃｋイベントハンドラ Private Sub 結果設定() With Worksheets("Sheet2") .Cells(1, 1) = "No.": .Cells(1, 2) = "X" .Cells(1, 3) = "Y": .Cells(1, 4) = "Z" For K = 0 To Num .Cells(K + 2, 1) = K: .Cells(K + 2, 2) = X(K) .Cells(K + 2, 3) = Y(K): .Cells(K + 2, 4) = Z(K) Next End With End Sub Sub Sheet2_ボタン1_Click() パルス FDCT IDCT 結果設定 End Sub

  18. VBAでの記述　③Ｃｈｅｎのアルゴリズム（その１）VBAでの記述　③Ｃｈｅｎのアルゴリズム（その１） Private Sub FDCT() 'FDCT(Chenのアルゴリズム) '(既に入力XがYに複写されているものとする） Cos4 = Cos(PI / 4): Cos8 = Cos(PI / 8) Cos16 = Cos(PI / 16): Cos163 = Cos(3 * PI / 16) Sin8 = Sin(PI / 8): Sin16 = Sin(PI / 16) Sin163 = Sin(3 * PI / 16) For J = 0 To Num Y(J) = X(J) Next For J = 0 To 3 tmp = Y(J) + Y(7 - J): Y(7 - J) = Y(J) - Y(7 - J): Y(J) = tmp Next tmp = Y(0) + Y(3): Y(3) = Y(0) - Y(3): Y(0) = tmp tmp = Y(1) + Y(2): Y(2) = Y(1) - Y(2): Y(1) = tmp tmp = Cos4 * (-Y(5) + Y(6)) Y(6) = Cos4 * (Y(5) + Y(6)) Y(5) = tmp

  19. VBAでの記述　④Ｃｈｅｎのアルゴリズム（その２）VBAでの記述　④Ｃｈｅｎのアルゴリズム（その２） tmp = Cos4 * (Y(0) + Y(1)) Y(1) = Cos4 * (Y(0) - Y(1)) Y(0) = tmp tmp = Sin8 * Y(2) + Cos8 * Y(3) Y(3) = -Cos8 * Y(2) + Sin8 * Y(3) Y(2) = tmp tmp = Y(4) + Y(5): Y(5) = Y(4) - Y(5): Y(4) = tmp tmp = -Y(6) + Y(7) Y(7) = Y(6) + Y(7) Y(6) = tmp tmp = Sin16 * Y(4) + Cos16 * Y(7) Y(7) = -Cos16 * Y(4) + Sin16 * Y(7) Y(4) = tmp

  20. VBAでの記述　⑤Ｃｈｅｎのアルゴリズム（その３）VBAでの記述　⑤Ｃｈｅｎのアルゴリズム（その３） tmp = Cos163 * Y(5) + Sin163 * Y(6) Y(6) = -Sin163 * Y(5) + Cos163 * Y(6) Y(5) = tmp tmp = Y(1): Y(1) = Y(4): Y(4) = tmp tmp = Y(3): Y(3) = Y(6): Y(6) = tmp For J = 0 To Num - 1 Y(J) = 0.5 * Y(J) Next Y(Num) = Y(0) End Sub

  21. VBAでの記述　⑥確認用逆変換 Private Sub IDCT() 'IDT(離散的逆変換) Dim C0 As Double: Dim C As Double C0 = 1# / Sqr(2#): C = Sqr(2# / Num) For J = 0 To Num - 1 Z(J) = C0 * Y(0) For K = 1 To Num - 1 Z(J) = Z(J) + Y(K) * Cos((2 * J + 1) * K * PI / (2 * Num)) Next Z(J) = C * Z(J) Next Z(Num) = Z(0) End Sub

  22. ＦＤＣＴの結果

  23. 色々なＦＤＣＴの結果（１）

  24. 色々なＦＤＣＴの結果（２）

More Related