(57)【特許請求の範囲】
【請求項１】
信号を量子化する方法であって、
ある入力データ範囲を有する入力データ信号を受信するステップと、
前記入力データ範囲を所定数の量子化段階に分割するステップと、
前記所定数の量子化段階を複数の量子化状態およびレベルに割り当てるステップと、
前記量子化された状態およびレベルに誤差最小化アルゴリズムを実行するステップと、
前記入力データ信号から各量子化状態およびレベルの重心値を適応的に計算するステップと、
前記の計算された重心に基づき前記量子化状態およびレベルを更新するステップと、
重心の所定の最適化が達成されるまで、前記分割するステップから前記更新するステップまでのステップを反復するステップと、
ひとたび重心の所定の最適化が達成されると、前記の更新された重心を量子化出力として送信するステップと
を備える方法。
【請求項２】
前記量子化された状態およびレベルに誤差最小化アルゴリズムを実行する前記ステップが、
サンプリング・データの各入力に関して前記量子化状態およびレベルを示す指標を格納するステップと、
各量子化状態およびレベルに関して前記の格納された指標位置の行列を形成するステップと
を含む請求項１記載の方法。
【請求項３】
前記最小誤差アルゴリズムが最小平均二乗（ＬＭＳ）アルゴリズムから成る請求項１記載の方法。
【請求項４】
前記重心の所定の最適化が所定の誤差閾値に基づく請求項１記載の方法。
【請求項５】
信号対雑音比（ＳＮＲ）値または値の変化が前記の計算された重心に基づいて決定される請求項１記載の方法。
【請求項６】
前記重心の所定の最適化が前記ＳＮＲ値を所定の閾値と比較することに基づいており、
前記方法が更にＳＮＲ値又は当該値の変化が所定のスレッショルドに到達するまで、最小誤差アルゴリズムを再実行し、重心値を適応的に計算し、前記量子化状態及びレベルを更新し、及びＳＮＲ値を獲得する各ステップを繰り返すステップを有する
請求項５記載の方法。
【請求項７】
前記ＬＭＳアルゴリズムがビタービ・アルゴリズムである請求項３記載の方法。
【請求項８】
デジタル処理システムにおいて、ある入力データ範囲を有する入力データ信号を量子化する方法であって、
前記入力データ信号を受信するステップと、
前記入力データ信号にトレリス経路の最適化を実行して、最小誤差経路の履歴を生成するステップと、
前記最小誤差経路に基づき経路行列を形成するステップと、
重心を前記経路行列に基づいて計算するステップと、
前記トレリス経路最適化を再度実行して、前記入力データ範囲のそれぞれの部分が前記量子化状態およびレベルに不均一に割り当てられる量子化状態およびレベルの不均一分布を生成するステップと、
前記の計算された重心を量子化出力として送信するステップと
を備える方法。
【請求項９】
前記データ信号がデシメートされたデータ信号である請求項８記載の方法。
【請求項１０】
前記トレリス経路の最適化がその鏡像である別の構成で実行される請求項８記載の方法。
【発明の詳細な説明】
（関連出願）
本出願は、モトローラ社に譲渡された１９９７年３月２０日にYolanda Prieto出願の米国出願番号第０８／８２２，４０３号「Adaptive Filtering For Use With Data Compression and Signal Reconstruction」および１９９７年３月２０日にYolanda Prieto出願の米国出願番号第０８／８２２，４０４号「Data Compression System, Method, and Apparatus」に関連する。
【０００１】
（産業上の利用分野）
本発明は、一般にデータ圧縮に関し、さらに詳しくは、デジタル・データの圧縮に関する。とりわけ、本発明は、トレリス符号化量子化器を利用するシステムにおいて改善された量子化を行う方法および装置に関する。
【０００２】
（従来の技術）
テレビ会議およびデジタル画像の格納などに関する技術やサービスの出現に伴い、デジタル信号処理の分野では大きな進歩がなされた。当業者には明白であろうが、デジタル信号処理は、通常は、サンプリングされたデータ信号を生成し、格納および／または送信を行うために信号を圧縮し、その後で被圧縮信号から元のデータを再構築するシステム，装置および方法に関する。効率が良く、費用効果が優れたデジタル信号処理システムに重要なことは、圧縮を行うために用いる方法である。
【０００３】
当技術においては周知の如く、データ圧縮とは、元のデータ信号をチャネル上での通信や適切な媒体での格納に適したビット・ストリーム内にマッピングする（割り付ける）ために実行されるステップを指す。元のデータを表現し回復するために必要な情報量を最小限に抑えることのできる方法が、演算の複雑性を軽減し、コストを低く抑えるために望ましい。コストに加えて、最小限の遅延をもって品質の良いデータ再生を行うことのできるハードウェアおよびソフトウェア設備は簡単であることが同様に望ましい。
【０００４】
量子化とは、データのビット・ストリームを取り出して、後で再生するためにそれを圧縮する技術を指す。データ信号の圧縮を行うために用いることのできるアルゴリズムはいくつか存在する。最も初歩的な方法では、信号サンプルを取り出し、各信号サンプルをより小さい数のレベルの１つに量子化する。周知の量子化器には、最近隣最小二乗誤差の原則（nearest neighbor, minimum square erro (MSE) rules）を用いて入力を指定された数のレベルに分類するロイド−マックス最適量子化器（Lloyd-Max optimal quantizer）がある。
【０００５】
量子化は、ベクトル，スカラー，丸め，切り捨ておよび特にトレリス符号化量子化（TCQ: trellis coded quantization）を含む種々の他の方法を用いて実行することができる。ベクトル量子化は、サンプル群を取り出して、ただ１つのシンボルによりその群を表現して、圧縮を改善し、高い信号対雑音比（SNR: signal-to-noise）を生成するので、効率的な量子化手段と考えられる。しかし、ベクトル量子化は、計算が集約的でルックアップ・テーブルを必要とする（すなわちメモリを要する）。スカラー量子化は、１シンボルにつき１サンプルを表現して、ベクトル量子化よりも計算負担が少ない。しかし、スカラー量子化には圧縮速度が遅いという欠点がある。量子化の代替法として出力信号の丸めまたは切り捨てがあるが、これらは分解能の問題を招くことがある。
【０００６】
現在のところ、静止画像のためのJPEG2000（Joint Photographic Experts Group）システムの次の標準は、ウェーブレット（wavelet）を用いて入力信号の分解を行うアルゴリズムと、信号圧縮のためのトレリス符号化量子化器とを提案する。
【０００７】
トレリス符号化量子化は、最小誤差経路を生成するための検索を行うビタービ・アルゴリズムに基づく。一般に、ビタービは、格子（lattice）量子化器（均一格子）を利用する。格子量子化は、演算効率が良いがその結果は歪みにより最適ではない。圧縮率が高くなると（ビット／サンプルが低い）、SNRが悪くなる。そのため、データ信号圧縮中に発されるノイズ（量子化ノイズ）を低減し、データの再構築の際に達成可能な信号対雑音比（品質）を改善するように動作する改善されたTCQ方法および装置に対する必要性が依然として存在する。
【０００８】
従って、特にデジタル信号処理システムに適応される改善されたトレリス符号化量子化（TCQ）法が必要である。
【０００９】
（好適な実施例の説明）
本明細書は、新規と見なされる本発明の特徴を定義する請求項を結論とするが、以下の説明と図面とを関連して考察することにより本発明がより良く理解頂けると信ずるものである。図面内では、同様の参照番号が継続して用いられる。
【００１０】
ここに開示する本発明は、データ・コデックの一部として実現することのできる量子化段階を改善する。図１を参照して、本発明によるデータ・コデック・ブロック図１００が図示される。このデータ・コデック１００は、ビデオなどの三次元（3-D）システム，静止画像などの二次元（2-D）システムならびに音声などの一次元（1-D）システムに適応することができる。
【００１１】
データ・コデック１００は、エンコーダ１１０とデコーダ１２０とを具備する。要するに、所定のビット／サンプル分解能を有する入力データ信号S(n)１０２が分解段１１２で変換を受ける。被変換信号x(n)１１３は、次に量子化段１１４を通過して、本発明によりさらなる圧縮を受ける。本発明により、量子化段１１４は、重心準拠量子化（centroid based quantization）を行って不均一な間隔の量子化器を生成する。被量子化信号１１５は、次に好ましくは周知の種々の符号化法の１つを用いる符号化段１１６で符号化される。被符号化信号１１８は、次にチャネル１３０上で通信されるか格納される。受信側において、デコーダ１２０は逆符号化段１３２，逆量子化段１３３および逆分解段１３４を用いて、入力データ信号S(n)を出力データ信号S'(n)１３６として再現する。
【００１２】
分解段１１２は、基本的には、ウェーブレット，FFT（Fast Fourier Transform: 高速フーリエ変換），DCT（Discrete Cosine Transform: 離散コサイン変換），DFT（Discrete Fourier Transform: 離散フーリエ変換），DST（Discrete Sine Transform: 離散サイン変換），KLT（Karhunen-Loewe Transform: カールネン−ロエベ変換），WHT（Walsh-Hadamard Transform: ウォルシュ−アダマール変換）など種々の周知の方法の１つによりデータ変換を実行する周波数領域アナライザである。ここで用いる変換とは、いわゆる部分帯域符号化（subband coding）である。部分帯域符号化においては、1-Dであろうと2-Dであろうと3-Dであろうと、入力データ信号１０２は、反復的に、低域通過および高域通過フィルタのフィルタ・バンクを通じて分割され、且つデシメート（ｄｅｃｉｍａｔｅ）され、元の信号の部分帯域被選択係数の出力を生成する。被変換信号x(n)１１３は、このため被変換係数の信号とも呼ばれる。また、被変換信号x(n)１１３は、デシメートされた信号でもある。本発明の好適な実施例においては、分解段１１２は、JPEG2000システムの次標準において提案されるようなウェーブレットにより構成されることもある。ウェーブレットは、一定の多重分解能機能を達成することを助け、画像データの場合のブロック効果を最小限に抑える。
【００１３】
符号化段１１６と逆符号化段１３２は、算術，ハフマン（Huffman）やその他の周知の符号化および解読法などを含む種々の符号化法を通じて動作することができる。
【００１４】
チャネル１３０は、たとえば無線周波数（RF）または光ファイバ伝送経路などのワイヤレス，ワイヤライン，音響または光学的リンク経路を含む。チャネル１３０は、同様に、周知の、あるいは将来開発される格納媒体とすることもできる。
【００１５】
逆量子化段１３３は、トレリス符号化量子化，スカラー，ベクトルなどの種々の逆量子化法または適応重心準拠量子化器１１４を逆にしたものを通じて動作することができる。逆分解段１３４は、好ましくは、逆ウェーブレット変換などの周知の手段またはフィルタ・バンクの集合により制御される。
【００１６】
本発明により、量子化段１１４は、係数の全範囲（部分帯域の全部または一部を表す）を取り、状態内で量子化器レベルの最適不均一割当を生成する重心準拠量子化器である。図２を参照して、本発明の好適な実施例による４状態トレリス符号化量子化（TCQ）段１１４のブロック図を示す。要するに、本発明により、トレリス経路最適化２０２が、被変換係数信号x(n)１１３に関して実行されて、経路行列２０６として構築される出力が生成され、これに関して重心演算が段２０８で実行される。重心は、この後で、トレリス経路最適段２０２に帰還されて、経路行列の別の反復および生成と更新された重心演算のための量子化レベルの新しい適応中心として用いられる。経路行列２０６は、ベクトル要素「nk,j」で構成され、ｋはＫ個の状態の総集合からの量子化状態を表し、ｊは１状態あたりＪ個のレベルの総集合からのレベルを表す。本発明により、トレリス経路最適化２０２，経路行列構築２０６および重心の更新２０８のステップは、適切なSNRまたは他の所望の誤差測定値，エントロピまたはビット速度が達成される２１０まで適応的に反復される。
【００１７】
本発明により、トレリス経路最適化２０２は、好ましくは、誤差が最小になる経路を得るために量子化器間の移行を選択する「最小平均二乗（LMS：least mean square）」アルゴリズムを通じて行われる。しかし、最小−最大などのその他の最適化法を用いても最小誤差経路を得ることができる。
【００１８】
経路最適化２０２の基本的実行例を好適な実施例２０９に示す。ただしこのとき、状態（各状態がＪ個のレベルを有する。）と呼ばれるＫ個の量子化器（量子）の集合に関して、量子（状態およびレベル）は、入力データの全範囲をカバーする。状態およびレベルは、より大きな量子化器を、初期に均一に分布する部分集合に分割した後で得られる。各状態の量子化レベルは、中点により、あるいは重心により、レベル間の量子化間隔を定義する。
【００１９】
入力データ・サンプルに対応する最も良いレベルと状態とを探す過程で、最適化を行うことにより、データ・サンプルが入力されて、歪み，誤差値，エントロピまたはビット速度を最小限にする状態とレベルの最良のシーケンスが発見される（２１０）。トレリスは特定のデータ集合または集合群に関して最適化され、いくつかの状態を伴う適応量子化器となる。これらは各々がより小型で最適化される量子化器となる。状態シーケンスと特定のデータ・シーケンスとの間には一致が見られる。最適化の後で異なるデータ集合が入力されると、確率密度関数（pdf：probability density function）またはヒストグラムは同一であっても、歪み削減の多くが消えて、新たな最適化が実行され、適応手順となる。Ｋ状態のトレリス量子化器に関して、各着信データ・サンプルにつき使用可能なＪレベルのうちから１つだけを選択する。これにより、Ｍがトレリス内の各状態に関して量子化段階の全体数を表すとすると、Ｍ／Ｋのレベルが割り当てられる。第１入力データ・サンプルに関する第１状態およびレベルは、誤差が最も低くなる状態を発見することにより、独立して決定される。続いて、各着信入力値が図２に示されるテーブル２０９などの移行テーブルに従い全部でＫ個の状態の２f状態のうちの１つに割り当てられる。従って、デコーダに対してたどるべき移行経路を伝えるには「ｆ」ビットかかる。たとえば、２０９に示されるように、ｆ＝１であると、Q０から、Q０（またはQ1）により次のサンプル入力を量子化することができ、さらにどのレベルを用いるかをデコーダに伝えるには一定のビットを要する。また、４状態トレリス経路最適化２０９または他のＫ状態トレリスを、基本的に鏡像である別の構成で実行することができる。
【００２０】
例として、ビタービ・アルゴリズムとして知られるLMSアルゴリズムが説明される最小誤差アルゴリズムとなる。ビタービ・アルゴリズムは、全範囲のデータ、すなわちx(n)１１３を取り出し（ただしｎは０からN-1とする）、このデータを第１回反復の一部として量子化して、データを各状態に関して異なるレベルに割り当てることができる。量子化器間にあるデータであれば、最も近いレベルに割り当てられる。
【００２１】
ビタービ・アルゴリズムは、量子化器間のすべての異なる経路に関して最小誤差を発見する「最小平均二乗（LMS）」法である。ビタービ・アルゴリズムは、トレリスを通る経路を最適化する。すなわち量子化器集合｛Qi｝をもつｉ番目の反復においては経路Πi（Qi）となり、許容できる移行に従う。たとえば、図２に示すビタービを参照して、t=0のサンプルx(n)がQ0(0)に量子化され次にt=1で移行はQ0(1)またはQ1(1)のいずれか一方にしか進めない。LMSアルゴリズムの最後に、最小誤差を生むためにとられる経路が簡潔に導かれ、三次元経路行列２０６の形で表現される。
【００２２】
経路行列２０６から、段２０８で重心が計算され、それを用いて、特定の量子化状態およびレベルに入るデータ・サンプルを表す。重心演算は、以下の等式により決まる平均値に基づく：
【００２３】
【数１】
0003657517.tif 000002
ただし
【数２】
0003657517.tif 000003
は、ベクトルであり、
【数３】
0003657517.tif 000004
が空のときは次のようになる：
【００２４】
【数４】
0003657517.tif 000005
新たに計算され更新される重心Ck,j(i+1)がトレリス経路最適化段２０２に戻され、そこでｋ番目の状態のｊ番目のレベルに関する新たな中心として用いられ、次の(i+1)番目の反復が開始される。これを図２の帰還矢印により示す。
【００２５】
一例として、x(n)３０４の入力データ列が１０個のサンプルからなるとする：
入力データ＝  x(n)  ＝  ［１０．０，１０．２，０．５，２．１８，１．９，１３．７５，１０．４０，０．１，１１．２，２．５］
0003657517.tif 000006
また、［０，１，２，３，．．．９］はそれぞれ特定のデータ・サンプル値を指す指標を表すものとする。（データ・サンプルの実際の値は、被変換係数x(n)であることに注意。）
次に、被変換係数入力x(n)１１３がトレリス２０２を用いて量子化される。この例では、トレリス２０２は４つの状態［Q0，Q1，Q2，Q3］を有し、各状態が４つのレベルを有するよう構築される（従って、１６の量子化段階があり、計算すべき重心は１６ある）。４状態４レベルのトレリスについて図示および説明されるが、トレリス２０２は、２つ以上の状態で構築され各状態がいくつかのレベルを持つ様々な構造に形成することができることは当業者には理解頂けよう。しかし、一般にはトレリス量子化器は、最も単純なケースである４状態を利用する。
【００２６】
この場合も、トレリス２０２の４状態は、量子化器Q0，Q1，Q2，Q3であり、各状態に４つのレベルが割り当てられる。第１回反復に関して、ビタービ・アルゴリズムを用いて最小誤差経路を量子化割当において得るものとするが、他のLMSアルゴリズムも同様に利用することができる。第１回反復に関して、レベルの均一な分布が用いられるが、第１回反復は、不均一なレベル分布を用いても同様に実行することができることを当業者には理解頂けよう。たとえば、量子化状態およびレベルを最良に割り当てる入力データの分布を用いる初期量子化器割当の場合、不均一な分布のレベルのほうが、より適している場合がある。たとえば、ラプラシアン型分布（すなわち鋭いピークを持つゼロ平均分布）を有する入力サンプル範囲については、縮小因子（shrink factor）を量子化段階に適用して、不均一な分布のレベルを得ることができる。
【００２７】
次に図３を参照して、i=0の第１回反復を行う４量子化状態の各々のレベル分布をグラフに示す。グラフ３００は、１６個の量子化段階に分かれたデータ入力範囲の４量子化状態Q0，Q1，Q2，Q3対４レベルの第１回反復の可能な例を表す。
【００２８】
以下の表は、入力データの様々な範囲と、特定の範囲に関する量子化状態およびレベルの割り当てを示す（ただし”［”は両端の値を含めることを表し、”（”は含めない）。
【００２９】
【表１】
0003657517.tif 000007
図２のビタービ・トレリスに戻って、トレリスはサンプルx(ni)を取り、サンプルが時刻t=tiにおいて状態Q0（またはQ2）の４レベルのいずれかに量子化された場合は、時刻t=ti+1における次のサンプルx(ni+1)だけを、状態Q0またはQ1の４レベルのいずれかに量子化することができる。同様に、サンプルx(ni)が時刻t=tiにおいて状態Q1（またはQ3）の４レベルのいずれかに量子化されると、時刻t=ti+1における次のサンプルx(ni+1)だけが状態Q2またはQ3の４レベルのいずれかに量子化することができる。
【００３０】
そのため、上記の入力データx(n)に関しては、データはサンプリングされたデータ値のみならず、以前の量子化器の状態にも基づいてLMS段２０２を通じて量子化される。経路行列２０６は、入力データを指示する指標をn(state, level)の種々の行列位置内に格納することによって構築される。次の表は、サンプリングされたデータがどのように量子化され、指標が行列位置に割り当てられるかを示す。
【００３１】
【表２】
0003657517.tif 000008
各行列位置に格納される指標は、元のデータのLMS（最小平均誤差）（すなわち最小誤差の経路）を表す。そのため、（この例では均一分布から始まる）第１回反復を進めると、行列２０６は以下の被格納指標で構成される：
【００３２】
【表３】
0003657517.tif 000009
格納される指標の図を図４に示し、経路行列２０６が実際に量子化器状態，レベルおよび指標の3-D構造を表すのを示す。各入力データ点は、サンプル入力点が起こる位置を示す指標により表される。図２の経路行列２０６において、各要素 は、その要素が状態ｋ，レベルｊに量子化された特定の入力データ・サンプルの指標であるベクトルによって構成さる。状態ｋ，レベルｊに量子化される入力サンプルがなかった場合、この行列要素
【数５】
0003657517.tif 000010
は空になる。
【００３３】
i=0の第１回反復を進めて、次のステップは重心を計算することである。上記の等式（１），（２）に示されるように、重心の計算は平均値に基づく。このため、この例の第１回反復（t=i）の重心演算により、次が得られる：
n00に関して  →  c00= x[7] = 0.1
および
c01＝ 以前の設定通り
c02＝ 以前の設定通り
n03に関して  →  c03= x[8] = 11.2
およびn10に関して  →  c10=1/2[x(3)+x(9)]=1/2(2.18+2.5)=2.34
（この反復(t=i-1）前は、重心c10は１．５であったことに注意）
c11＝ 以前の設定通り
n12に関して  →  c12= x[1] = 10.2
n13に関して  →  c13= x[5] = 13.75
c13＝ 以前の設定通り
および
n20に関して  →  c20=1/2[x(2)+x(4)]=1/2(1.5+1.5)=1.2
（この反復(t=i-1）前は、重心c20は２．５であったことに注意）
c21＝ 以前の設定通り
および
n22に関して  →  c22=1/2[x(0)+x(6)]=1/2(10.0+10.4)=10.2
（この反復(t=i-1）前は、重心c22は１０．５であったことに注意）
c23＝ 以前の設定通り
c30＝ 以前の設定通り
c31＝ 以前の設定通り
c32＝ 以前の設定通り
c32＝ 以前の設定通り
c33＝ 以前の設定通り
かくして、図３の例に示されるように、以前の重心がすべて０．５を中心にある場合、量子化器の各状態およびレベル（現在は不均一な）について新しい重心が得られる。同じ集合のデータx(n)が第２回反復において新たな重心をもって再分類され、この反復プロセスは、特定の入力信号に関する誤差測定値が所定の所望値に達するまで繰り返される。図５は、１回または数回の適応反復の後で、種々の状態およびレベルが均一分布から、入力データに最適に適合する分布に変化する様子の任意の例を示す。本発明により説明される方法は、各反復においてより低い誤差への収束を保証する。
【００３４】
このように、最小平均二乗アルゴリズムによるトレリス量子化のステップを通じてサンプリングされたデータ係数x(n)を取り出して行列を形成し、各量子化状態に関して更新された重心を計算することによって、プルーニング（pruning）アルゴリズムがもたらされる。このアルゴリズムは、所望のSNRを生成するのに適した収束が得られるまで再度繰り返される。新しい部分帯域または入力データの新しい集合の場合、前回のデータ実行から計算された重心値を用いて第１回反復を開始するように量子化器を調整することもできる。この方法では、次の入力データ範囲内で同様の入力データが予想される場合は、最適化されるのにはるかに近い重心値で量子化を開始することができ、そのためにSNRを全面的に最適化するのに必要な反復数が少なくなる。
【００３５】
従来技術の場合は、量子化法は基本的にはビタービ・アルゴリズム段で停止した。これは、量子化レベル間のデルタが均一であることを前提としたためである。本発明により説明される出願人の量子化法１１４を適用することにより、重心演算は量子化レベル間の距離が変わるとダイナミックに調整される（図３および図５を比較することでわかる）。このため、従来技術は均一に間隔をとるシステム上で動作したが、本発明の量子化法では、不均一性と入力データに基づく適応が可能になる。
【００３６】
本発明の量子化段１１４は、更新された重心を用いて量子化応対のレベルを適応的に割り当てる。このため、量子化器はデータ分布に最適化されて、SNRを改善する。量子化器の出力１１５は、段１１６で符号化され、チャネル１３０に送信される。これにより量子化器出力１１５は、着信信号に基づいて適応的に割り当てられる。本発明により説明される適応重心に基づく量子化器は、画像信号などの入力信号を、入力信号に基づいて量子化を適応的に割り当てることにより最適化される。量子化ブロック１１４は、テレビ会議中の顔の画像などの特定の種類の信号に関して、量子化のステップをある所定の状態およびレベル割当において開始して処理時間を短くするように調整することもできる。
【００３７】
以下のステップは、本発明による適応重心準拠量子化法をまとめる：
i)  入力データ範囲が所定数の量子化段階に分割される。上記に説明する例では、４ビットのトレリスについて１６の量子化段階がある。
【００３８】
ii)  量子化段階が状態の行列に形成される。上記の例では、量子化段階１，５，９，１３が状態１に進み、量子化段階２，６，１０，１４が状態２に進むなど。これが初期状態割当である。
【００３９】
iii)  入力データが、ラプラシアン分布など所定の分布をもつと予想される場合は、前述の要領で縮小因子を導入することができる。
【００４０】
【数６】
0003657517.tif 000011
iv)  LMSアルゴリズム、好ましくはビタービ・アルゴリズムを伴うトレリスを実行する。
【００４１】
v)  平均二乗誤差を計算することにより最初の生SNR値または他の誤差測定値を求める。
【００４２】
【数７】
0003657517.tif 000012
vi)  データ・サンプルの重心値を、状態毎に各レベルに１つ計算する。
【００４３】
vii)  新たな平均二乗誤差を計算することにより、精密なより高いSNR値を求める。
【００４４】
viii)  ステップ(iv)からの更新重心値を用いて、LMSアルゴリズムを伴うトレリスの２回目を実行する。
【００４５】
viv)  重心を再計算し、SNRの新たな値を求める。
【００４６】
SNRの最後の値は最初のものよりも良くなっている。SNRは、重心の更新により各反復サイクル毎に収穫逓減的に増加する。
【００４７】
x)  所定のSNR規準が満たされたら反復を停止する。
【００４８】
本発明により開示される適応重心準拠量子化法の利点には、収束だけでなく反復処理時間が（傾斜法に比べて）早いことがあげられる。
【００４９】
上記のステップで１つのオプション・パラメータ、すなわち縮小因子が導入された。入力データがゼロ平均ラプラシアンなどの鋭い分布を有する場合、最外側レベルにはデータ・サンプルはほとんどない。従って、縮小因子を適用して、すべてのレベルを入力範囲のより活動的な中心に向かって集中させることができる。その結果得られるSNRは、反復のたびに単調に増加して、最適なSNRレベルに収束する。縮小因子の効果は、最適なSNR値を得るための総反復回数を減らすことである。
【００５０】
開示されるシステム明細は、さらに、変換段，量子化段，符号化，チャネル，逆符号化および逆分解段を含む。本発明の別の実施例においては、量子化１１４に先立ってベクトル化（2Dを1Dデータ・ストリームに変換）と閾値段とが実行されて、反復データの長いストリングを１つの符号で表し、非ゼロ1Dデータ・ストリームを得られるようにする。図６には、データx(n)が段６０２において1-Dにベクトル化されてベクトル係数６０４を生成する実行例が提案される。ベクトル係数６０４は、次に量子化に先立って、閾値６０６と比較される。基本的には、ゼロのデータ・サンプルおよび／またはゼロに近い値は閾値６０６によって切り捨てられて、非ゼロ要素の修正された一次元入力データ・アレイを適応重心準拠量子化器１１４の入力として得る。
【００５１】
図７は、完全なコデック・システム６００の一部としての1-Dデータ・ストリーム変換ステップを示す。量子化されると、好ましくは符号化段１１６において符号化される反復データがあり、さらに圧縮される。これにより、送信時間と計算時間とが短縮される。この場合も、ダイナミック重心準拠量子化１１４を通じて得られる利点を、システム７００の実施例で実現することができる。これには被符号化信号を通じて反復データの効率的な送信という更なる利点も伴う。
【００５２】
本発明の好適な実施例が図示および説明されたが、本発明がそれに制限されないことは明らかである。添付の請求項に定義される本発明の精神および範囲から逸脱せずに数多くの修正，変更，変換，置き換えおよび同等物が当業者には可能であろう。
【図面の簡単な説明】
【図１】  本発明の好適な実施例によるデータ・コデックのブロック図である。
【図２】  本発明による適応重心準拠量子化器である。
【図３】  本発明による図２の量子化器に関して、i=0の第１回反復を行う４つの量子化器状態のうち個々の状態に関するレベル分布例をグラフに示す。
【図４】  本発明による図２の量子化器に関して第１回反復から格納される指標を図に表す例である。
【図５】  本発明による図２の量子化器に関して、任意の反復を行う４つの量子化器状態のうち個々の状態に関するレベル分布例をグラフに示す。
【図６】  本発明の代替実施例による量子化前の一次元データ・ストリーム変換の提案される実行例である。
【図７】  本発明の代替実施例による完全なデータ・コデック・システムの一部としての1-Dデータ・ストリーム変換の段を示す。
【図１】
0003657517.tif 000013
【図２】
0003657517.tif 000014
【図３】
0003657517.tif 000015
【図４】
0003657517.tif 000016
【図５】
0003657517.tif 000017
【図６】
0003657517.tif 000018
【図７】
0003657517.tif 000019