ウェーブレット

ウェーブレット（英: wavelet）やマザーウェーブレット（英: mother wavelet）とは、数学において、局在する波、つまり、有限の長さの波もしくは速やかに減衰する波の事。ファーザーウェーブレット（英: father wavelet）とは、多重解像度解析にて使われる、マザーウェーブレット関数とセットで使われるスケーリング関数の事。waveletはwave（波）とlet（小さい）の合成語である。 ウェーブレット変換・ウェーブレット解析とは、ウェーブレットを用いて変換・解析する事。信号表現は入力信号に合致するようなウェーブレット波形の拡大縮小（スケーリング）・平行移動（シフト）により行われる。より正確には、この信号表現はウェーブレット系列と呼ばれ、これは2乗可積分関数のヒルベルト空間における完全正規直交系の基底関数集合（正規直交基底）を用いた線形基底展開である。 ウェーブレットの発展は、20世紀初頭のハンガリー人数学者ハール・アルフレッドによるいくつかの断片的な考察に基づく。ガーボル・デーネシュによるその後の研究でガボール・アトム（英語版）が得られた。ガボール・アトムはウェーブレットと似た形で構成され、似た目的に応用された。ウェーブレット理論への大きな貢献のひとつは、1975年のジョージ・ツワイクによる連続ウェーブレット変換（英語版）（初期には cochlear 変換と呼ばれていて、耳の音に対する反応を研究していたときに発見された。） の発見である。 ウェーブレットの概念は、1975年にエルフで石油探査をしていたフランス人地球物理学者ジャン・モーレー（英語版）が発見した。1981年、モーレーは クロアチア系フランス人物理学者アレックス・グロスマン（英語版）との共同研究から連続ウェーブレット変換（英語版）の定式化（Goupillaud）を行なった。彼らはフランス語で"小さい波"を意味するondeletteという言葉を用いていたが、少し後に英語に翻訳された際に"onde"は"wave"と訳されてウェーブレット("wavelet")という用語が誕生した。 その後のウェーブレット理論における大きな貢献には、Strombergによる離散ウェーブレット変換における初期研究(1983年)、イングリッド・ドブシーによるコンパクト台を持つ直交ウェーブレット(1988年)、Mallatによる多重解像度解析に関する提案(1989年)、Delpratによる連続ウェーブレット変換における時間-周波数変換(1991年)、Newlandによるハーモニックウェーブレット変換（英語版）など、枚挙にいとまがない。 1909年: 最初のウェーブレット (Haarによるハールウェーブレット) 1950年代以降: Jean Morlet , Alex Grossman 1980年代以降: Yves Meyer, Stephane Mallat, Ingrid Daubechies, Ronald Coifman, Victor Wickerhauser ウェーブレット理論は、いくつかの異なる目的で応用される。 全てのウェーブレット変換は、時間周波数表現（英語版）であると考えられるが、調和解析とも関係がある。 ウェーブレット変換は、大きく連続ウェーブレット変換（英語版）(CWT)と離散ウェーブレット変換(DWT)に分類される。これらの違いは、連続ウェーブレット変換では可能な全てのスケールとシフトが用いられるのに対して、離散ウェーブレット変換では一部分のみが使われることにある。 連続ウェーブレット変換（英語版）は、ハイゼンベルクの不確定性原理に支配されている。同様に、離散ウェーブレットにおいても不確定性原理は考慮されなければならない。 多くの場合に有用である離散ウェーブレット変換は、有限インパルス応答(FIR)フィルタで構成されるフィルタバンクである。 ウェーブレット変換は３つに分類されることが多い。連続ウェーブレット変換（英語版）、離散ウェーブレット変換、多重解像度解析(MRA)による離散ウェーブレット変換である。以下、この３つについて解説する。 連続ウェーブレット変換（英語版）においては、有限なエネルギーを持った信号は、連続な周波数バンドの群（もしくは L 2 ( R ) {\displaystyle L^{2}(\mathbb {R} )} 関数空間の一部) として投影される。得られた周波数成分は、適切な積分によって元の信号を再構成することができる。 部分空間の群は、スケール1の部分空間を拡大縮小（スケール）して生成されたものである。この部分空間は、１つの関数すなわちマザーウェーブレット ψ ∈ L 2 ( R ) {\displaystyle \psi \in L^{2}(\mathbb {R} )} をシフトすることによって生成される。 一般的なマザーウェーブレットの例は以下のとおりである。 スケールaの部分空間は、以下の式で生成される。(これはベビーウェーブレットと呼ばれることがあるがあまり一般的ではない) ψ a , b ( t ) = 1 a ψ ( t − b a ) {\displaystyle \psi _{a,b}(t)={\frac {1}{\sqrt {a}}}\psi \left({\frac {t-b}{a}}\right)} ただし、aは正の実数でありスケールを決定する。bは任意の実数でありシフトを決定する。(a,b)のペアは、 R + × R {\displaystyle \mathbb {R} _{+}\times \mathbb {R} } の上半面において定義される。 関数xをスケールaの部分空間へ投影すると、以下の式で示される。 x a ( t ) = ∫ R W T ϕ { x } ( a , b ) ⋅ ψ a , b ( t ) d b {\displaystyle x_{a}(t)=\int _{\mathbb {R} }WT_{\phi }\{x\}(a,b)\cdot \psi _{a,b}(t)\,db} 但し、WTはウェーブレット係数である。 W T ϕ { x } ( a , b ) = ⟨ x , ψ a , b ⟩ = ∫ R x ( t ) ψ a , b ( t ) ¯ d t {\displaystyle WT_{\phi }\{x\}(a,b)=\langle x,\psi _{a,b}\rangle =\int _{\mathbb {R} }x(t){\overline {\psi _{a,b}(t)}}\,dt} . 信号xの解析のためには、ウェーブレット係数をスケーログラム（英語版）にする。 全てのウェーブレット係数を使って信号を解析することは実用上不可能である。信号を対応するウェーブレット係数から再構成することは，上半面の離散部分集合さえ取り出せば十分可能だと思うだろう。その一つとして実数パラメータa>1，b>0によるアフィン系がある。対応する半面の離散部分集合は，全ての点 ( a m , n a m b ) {\displaystyle (a^{m},n\,a^{m}b)} を含む（ m , n ∈ Z {\displaystyle m,n\in \mathbb {Z} } ）。対応するベビーウェーブレットは以下で与えられる。 ψ m , n ( t ) = a − m / 2 ψ ( a − m t − n b ) {\displaystyle \psi _{m,n}(t)=a^{-m/2}\psi (a^{-m}t-nb)} . 式 x ( t ) = ∑ m ∈ Z ∑ n ∈ Z ⟨ x , ψ m , n ⟩ ⋅ ψ m , n ( t ) {\displaystyle x(t)=\sum _{m\in \mathbb {Z} }\sum _{n\in \mathbb {Z} }\langle x,\,\psi _{m,n}\rangle \cdot \psi _{m,n}(t)} による有限エネルギーを持つ任意の信号xの再構成のための十分条件は，関数群 { ψ m , n : m , n ∈ Z } {\displaystyle \{\psi _{m,n}:m,n\in \mathbb {Z} \}} が L 2 ( R ) {\displaystyle L^{2}(\mathbb {R} )} のタイトフレーム（英語版）を形作ることである。 各種あるウェーブレット変換の離散化の全ての方法において、上半面上の各有界矩形領域は有限個の係数のみを持つ。しかし、各係数を求めるためには積分の評価が必要となる。このような数値的な複雑さを避けるために、スケーリング関数（ファーザーウェーブレット）と呼ばれる補助関数 ϕ ∈ L 2 ( R ) {\displaystyle \phi \in L^{2}(\mathbb {R} )} が利用される。このとき a は整数でなければならない。例えば典型的な係数として a=2、b=1 が用いられる。最も有名なファーザー・マザーウェーブレットの組としてドブシーの4タップウェーブレットがある。 マザー・ファーザーそれぞれのウェーブレットから部分空間 V m = span ⁡ ( ϕ m , n : n ∈ Z ) {\displaystyle V_{m}=\operatorname {span} (\phi _{m,n}:n\in \mathbb {Z} )} , where ϕ m , n ( t ) = 2 − m / 2 ϕ ( 2 − m t − n ) {\displaystyle \phi _{m,n}(t)=2^{-m/2}\phi (2^{-m}t-n)} と W m = span ⁡ ( ψ m , n : n ∈ Z ) {\displaystyle W_{m}=\operatorname {span} (\psi _{m,n}:n\in \mathbb {Z} )} , where ψ m , n ( t ) = 2 − m / 2 ψ ( 2 − m t − n ) {\displaystyle \psi _{m,n}(t)=2^{-m/2}\psi (2^{-m}t-n)} . が構成される。これらより、系列 { 0 } ⊂ ⋯ ⊂ V 1 ⊂ V 0 ⊂ V − 1 ⊂ ⋯ ⊂ L 2 ( R ) {\displaystyle \{0\}\subset \dots \subset V_{1}\subset V_{0}\subset V_{-1}\subset \dots \subset L^{2}(\mathbb {R} )} は L 2 ( R ) {\displaystyle L^{2}(\mathbb {R} )} の多重解像度解析を形成することになり、また部分空間 … , W 1 , W 0 , W − 1 , … … {\displaystyle \dots ,W_{1},W_{0},W_{-1},\dots \dots } は上の系列の直交する差分、つまり W m {\displaystyle W_{m}} は V m − 1 {\displaystyle V_{m-1}} 中にある V m {\displaystyle V_{m}} の直交補空間となる。サンプリング定理と同様に、sampling distance 2 m {\displaystyle 2^{m}} の空間 V m {\displaystyle V_{m}} は 0 から 2 − m − 1 {\displaystyle 2^{-m-1}} の周波数帯域をほぼカバーすることになる。また W m {\displaystyle W_{m}} は直交補空間として帯域 [ 2 − m − 1 , 2 − m ] {\displaystyle [2^{-m-1},2^{-m}]} を大まかにカバーする。 このような包含と直交の関係より，2組の恒等式 h n = ⟨ ϕ 0 , 0 , ϕ 1 , n ⟩ {\displaystyle h_{n}=\langle \phi _{0,0},\,\phi _{1,n}\rangle } ϕ ( t ) = 2 ∑ n ∈ Z h n ϕ ( 2 t − n ) {\displaystyle \phi (t)={\sqrt {2}}\sum _{n\in \mathbb {Z} }h_{n}\phi (2t-n)} と g n = ⟨ ψ 0 , 0 , ϕ 1 , n ⟩ {\displaystyle g_{n}=\langle \psi _{0,0},\,\phi _{1,n}\rangle } ψ ( t ) = 2 ∑ n ∈ Z g n ϕ ( 2 t − n ) {\displaystyle \psi (t)={\sqrt {2}}\sum _{n\in \mathbb {Z} }g_{n}\phi (2t-n)} を満たす系列 h = { h n } n ∈ Z {\displaystyle h=\{h_{n}\}_{n\in \mathbb {Z} }} と g = { g n } n ∈ Z {\displaystyle g=\{g_{n}\}_{n\in \mathbb {Z} }} が存在することになる。 2番目の恒等式はファーザーウェーブレット ϕ {\displaystyle \phi } の洗練条件（英語版）と呼ばれる。これらの恒等関係は高速ウェーブレット変換（英語版）アルゴリズムの土台となっている。 実応用での効率性を考えると、マザー（プロトタイプ）ウェーブレット（関数）はコンパクトサポートの連続微分可能関数であることが望ましい。しかし、（連続ウェーブレット変換における）解析的であることの要求と、理論的な理由から、一般的にウェーブレット関数は 空間 L 1 ( R ) ∩ L 2 ( R ) {\displaystyle L^{1}(\mathbb {R} )\cap L^{2}(\mathbb {R} )} の部分空間から選ばれる。これは絶対値積分可能かつ2乗積分可能な可測関数の空間である。 ∫ − ∞ ∞ | ψ ( t ) | d t < ∞ {\displaystyle \int _{-\infty }^{\infty }|\psi (t)|\,dt<\infty } and ∫ − ∞ ∞ | ψ ( t ) | 2 d t < ∞ {\displaystyle \int _{-\infty }^{\infty }|\psi (t)|^{2}\,dt<\infty } . この関数空間では必ずゼロ平均と二乗ノルムの条件が定式化できる。 ∫ − ∞ ∞ ψ ( t ) d t = 0 {\displaystyle \int _{-\infty }^{\infty }\psi (t)\,dt=0} （ゼロ平均の条件） ∫ − ∞ ∞ | ψ ( t ) | 2 d t = 1 {\displaystyle \int _{-\infty }^{\infty }|\psi (t)|^{2}\,dt=1} （二乗ノルム正規の条件） ψ {\displaystyle \psi } が 連続ウェーブレット変換（英語版）（正確な議論はリンク先参照）のウェーブレットであるためには、マザーウェーブレットは安定な逆変換を持つための許容性の規範（簡単に言うとこれは半微分可能性のようなもの）を満たさなければならない。 離散ウェーブレット変換における最低限満たさなければならない条件として、ウェーブレット系列はLp空間 L 2 ( R ) {\displaystyle L^{2}(\mathbb {R} )} 中の単位元でなければならない。離散ウェーブレット変換のほとんどの構成は多重解像度解析を用いており、この場合ウェーブレットはスケール関数により決定される。このスケール関数自体が汎関数方程式である。 多くの場合において ψ {\displaystyle \psi } をvanishing moments を表すより大きい数字 M の連続関数、つまり全ての整数 m<M について以下の式を満たす関数に限定することは有用である。 ∫ − ∞ ∞ t m ψ ( t ) d t = 0 {\displaystyle \int _{-\infty }^{\infty }t^{m}\,\psi (t)\,dt=0} マザーウェーブレットは、 a {\displaystyle a} の因数による拡大縮小（スケール）と、 b {\displaystyle b} の因数による平行移動（シフト）により、（Morlet によるオリジナルの定式化のように）以下のように与えられる。 ψ a , b ( t ) = 1 a ψ ( t − b a ) {\displaystyle \psi _{a,b}(t)={1 \over {\sqrt {a}}}\psi \left({{t-b} \over a}\right)} . 連続ウェーブレット変換では、(a,b) のペアは全半平面 R + × = R {\displaystyle \mathbb {R} _{+}\times =R} 上で変化する。また離散ウェーブレット変換では、このペアは、affine群とも呼ばれる離散部分集合上で変化する。 これらの関数はたびたび（連続）変換の基底関数という誤った捉え方をされる。事実、連続フーリエ変換にあるような基底は、連続ウェーブレット変換にはみあたらない。時間周波数解釈では少し違う定式化が使用される（Delpratによる）。 ウェーブレット変換は、三角関数の級数表現のフーリエ変換としばしば比較される。主な違いは、ウェーブレット変換は時間と周波数の両方の成分を局在化するが、標準的なフーリエ変換は周波数成分だけを局在化することである。短時間フーリエ変換も時間と周波数の両方を局在化できるが、時間周波数分解能に問題がある。一方、ウェーブレット変換ではしばしば多重解像度解析という、より良い表現が用いられる。 また、離散ウェーブレット変換の計算量はO(N)であり、高速フーリエ変換のO(N log N)に比べて小さい(ここで、Nはデータの大きさである)。 ウェーブレット（およびウェーブレット族）の定義の仕方には様々な方法がある。 ウェーブレットはもっぱらスケーリングフィルタ g により定義される。これは長さ 2N で和が 1 となる有限インパルス応答(FIR) の低域通過フィルタである。双直交ウェーブレットでは、分解と再合成のフィルタが別々に決定される。 分析では、高域通過フィルタは低域通過フィルタの QMF として計算され、再合成フィルタは分解フィルタの時間反転である。 例えば、ドブシー・ウェーブレット（英語版）とSymletウェーブレット（チェコ語版）は、スケーリングフィルタで定義することができる。 ウェーブレットは、ウェーブレット関数 ψ ( t ) {\displaystyle \psi (t)} (マザーウェーブレット)とスケーリング関数 ϕ ( t ) {\displaystyle \phi (t)} (ファーザーウェーブレット)から定義される。 実際のところウェーブレット関数は帯域通過フィルターであり、それぞれの水準の半分の帯域幅でスケールされている。これによって、全てのスペクトルを扱うために無限の水準が必要となる問題が生じる。スケーリング関数を用いれば、最低限の水準で全てのスペクトルを扱うことができる。詳細な説明は[1]にある。 コンパクトサポートをもつウェーブレットでは、 ϕ ( t ) {\displaystyle \phi (t)} は有限長であり、スケーリングフィルターgと同等である。 例えば、Meyerウェーブレットは、スケーリング関数で定義することができる。 ウェーブレットはウェーブレット関数 ψ ( t ) {\displaystyle \psi (t)} のような時間領域表現をただ一つ持つ。 例えば、Mexican hat waveletは、ウェーブレット関数で定義することができる。 いくつかの連続ウェーブレット（英語版）のリストを参照。 大まかに、離散ウェーブレット変換はデータ圧縮に使われる一方で、連続ウェーブレット変換は信号解析に使われる。その結果として、離散ウェーブレット変換は工学と計算機科学において一般的に使われ、連続ウェーブレット変換は科学研究においてもっともよく使われている。ウェーブレット変換は、現在非常に多くの様々な用途に、しばしば従来のフーリエ変換を置き換えて使用されている。分子動力学、第一原理計算、宇宙物理学、密度行列局在、地震地球物理学、光学、乱流そして量子力学を含む、物理学の多くの分野でこのパラダイムシフトが起こった。この変化が起こった他の分野は画像処理、血圧、心拍やECGの解析、DNA解析、タンパク質解析、気候学、一般的な信号処理、音声認識、コンピュータグラフィックスそしてマルチフラクタル解析（英語版）である。コンピュータビジョンや画像処理において、尺度空間（英語版）表現やガウス微分オペレータの概念は正規化された多重解像度表現の一つであると考えられている。 ウェーブレットはデータ圧縮の分野でも用いられる。デジタル信号処理における他の時間-周波数変換と同様、ウェーブレット変換は(たとえば画像などの)圧縮されていないデータに対し適用でき、その後圧縮処理がなされることで、結果として効果的なデータ圧縮を実現できる。JPEG 2000はウェーブレットを利用した画像形式の一つである。ウェーブレットを利用したデータ圧縮についてはウェーブレット圧縮を参照されたい。 JPEG 2000で使われているウェーブレットは双直交ウェーブレットであり、ウェーブレット系列の座標は異なる2つの基底関数集合を用いて計算されるため、注意を要する。 異なる用途に応じて、多くのウェーブレット変換が存在する。以下にその一例を列挙するが、すべてのウェーブレット変換についてのリストはウェーブレット変換の一覧（英語版）を参照されたい。 連続ウェーブレット変換（英語版） (CWT) 離散ウェーブレット変換 (DWT) 高速ウェーブレット変換（英語版） (FWT) ウェーブレットパケット分解（英語版） (WPD) 定常ウェーブレット変換（英語版） (SWT) Beylkin (18) Coiflet (6, 12, 18, 24, 30) ドブシー・ウェーブレット（英語版） (2, 4, 6, 8, 10, 12, 14, 16, 18, 20) Cohen-Daubechies-Feauveauウェーブレット（英語版） (CDF N/PまたはDaubechies双直交ウェーブレットとして参照されることがある) ハールウェーブレット Symletウェーブレット（チェコ語版） Mexican hat wavelet Hermitian wavelet Hermitian hat wavelet Complex mexican hat wavelet Morlet wavelet Modified Morlet wavelet Beta wavelet Hilbert-Hermitian wavelet Chirplet変換（英語版） Curvelet（英語版） ウェーブレット木（英語版） - データ構造 フィルタバンク 非整数次フーリエ変換 多重解像度解析 短時間フーリエ変換 尺度空間（英語版） 超広帯域ラジオはウェーブレットを送信する。