JP4429316B2 - 時間ドメインから周波数ドメインへ及びそれとは逆にデジタル信号のドメイン変換を実行する装置及び媒体 - Google Patents

時間ドメインから周波数ドメインへ及びそれとは逆にデジタル信号のドメイン変換を実行する装置及び媒体 Download PDF

Info

Publication number
JP4429316B2
JP4429316B2 JP2006532242A JP2006532242A JP4429316B2 JP 4429316 B2 JP4429316 B2 JP 4429316B2 JP 2006532242 A JP2006532242 A JP 2006532242A JP 2006532242 A JP2006532242 A JP 2006532242A JP 4429316 B2 JP4429316 B2 JP 4429316B2
Authority
JP
Japan
Prior art keywords
domain
matrix
lift
receive
transformation
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Expired - Fee Related
Application number
JP2006532242A
Other languages
English (en)
Other versions
JP2007507789A5 (ja
JP2007507789A (ja
Inventor
ハイビン ファン，
シャオ リン，
スザント ラハルジャ，
ロンシャン ユー，
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Agency for Science Technology and Research Singapore
Original Assignee
Agency for Science Technology and Research Singapore
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Agency for Science Technology and Research Singapore filed Critical Agency for Science Technology and Research Singapore
Publication of JP2007507789A publication Critical patent/JP2007507789A/ja
Publication of JP2007507789A5 publication Critical patent/JP2007507789A5/ja
Application granted granted Critical
Publication of JP4429316B2 publication Critical patent/JP4429316B2/ja
Anticipated expiration legal-status Critical
Expired - Fee Related legal-status Critical Current

Links

Images

Classifications

    • G—PHYSICS
    • G06—COMPUTING OR CALCULATING; COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10—Complex mathematical operations
    • G06F17/14—Fourier, Walsh or analogous domain transformations, e.g. Laplace, Hilbert, Karhunen-Loeve, transforms
    • G—PHYSICS
    • G06—COMPUTING OR CALCULATING; COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10—Complex mathematical operations
    • G06F17/14—Fourier, Walsh or analogous domain transformations, e.g. Laplace, Hilbert, Karhunen-Loeve, transforms
    • G06F17/147—Discrete orthonormal transforms, e.g. discrete cosine transform, discrete sine transform, and variations therefrom, e.g. modified discrete cosine transform, integer transforms approximating the discrete cosine transform
    • G—PHYSICS
    • G06—COMPUTING OR CALCULATING; COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10—Complex mathematical operations
    • G06F17/14—Fourier, Walsh or analogous domain transformations, e.g. Laplace, Hilbert, Karhunen-Loeve, transforms
    • G06F17/141—Discrete Fourier transforms
    • G—PHYSICS
    • G06—COMPUTING OR CALCULATING; COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10—Complex mathematical operations
    • G06F17/14—Fourier, Walsh or analogous domain transformations, e.g. Laplace, Hilbert, Karhunen-Loeve, transforms
    • G06F17/148—Wavelet transforms
    • G—PHYSICS
    • G10—MUSICAL INSTRUMENTS; ACOUSTICS
    • G10L—SPEECH ANALYSIS TECHNIQUES OR SPEECH SYNTHESIS; SPEECH RECOGNITION; SPEECH OR VOICE PROCESSING TECHNIQUES; SPEECH OR AUDIO CODING OR DECODING
    • G10L19/00—Speech or audio signals analysis-synthesis techniques for redundancy reduction, e.g. in vocoders; Coding or decoding of speech or audio signals, using source filter models or psychoacoustic analysis
    • G10L19/02—Speech or audio signals analysis-synthesis techniques for redundancy reduction, e.g. in vocoders; Coding or decoding of speech or audio signals, using source filter models or psychoacoustic analysis using spectral analysis, e.g. transform vocoders or subband vocoders
    • G10L19/0212—Speech or audio signals analysis-synthesis techniques for redundancy reduction, e.g. in vocoders; Coding or decoding of speech or audio signals, using source filter models or psychoacoustic analysis using spectral analysis, e.g. transform vocoders or subband vocoders using orthogonal transformation
    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04N—PICTORIAL COMMUNICATION, e.g. TELEVISION
    • H04N19/00—Methods or arrangements for coding, decoding, compressing or decompressing digital video signals
    • H04N19/60—Methods or arrangements for coding, decoding, compressing or decompressing digital video signals using transform coding
    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04N—PICTORIAL COMMUNICATION, e.g. TELEVISION
    • H04N19/00—Methods or arrangements for coding, decoding, compressing or decompressing digital video signals
    • H04N19/60—Methods or arrangements for coding, decoding, compressing or decompressing digital video signals using transform coding
    • H04N19/63—Methods or arrangements for coding, decoding, compressing or decompressing digital video signals using transform coding using sub-band based transform, e.g. wavelets
    • G—PHYSICS
    • G10—MUSICAL INSTRUMENTS; ACOUSTICS
    • G10L—SPEECH ANALYSIS TECHNIQUES OR SPEECH SYNTHESIS; SPEECH RECOGNITION; SPEECH OR VOICE PROCESSING TECHNIQUES; SPEECH OR AUDIO CODING OR DECODING
    • G10L19/00—Speech or audio signals analysis-synthesis techniques for redundancy reduction, e.g. in vocoders; Coding or decoding of speech or audio signals, using source filter models or psychoacoustic analysis
    • G10L19/0017—Lossless audio signal coding; Perfect reconstruction of coded audio signal by transmission of coding error
    • G—PHYSICS
    • G10—MUSICAL INSTRUMENTS; ACOUSTICS
    • G10L—SPEECH ANALYSIS TECHNIQUES OR SPEECH SYNTHESIS; SPEECH RECOGNITION; SPEECH OR VOICE PROCESSING TECHNIQUES; SPEECH OR AUDIO CODING OR DECODING
    • G10L19/00—Speech or audio signals analysis-synthesis techniques for redundancy reduction, e.g. in vocoders; Coding or decoding of speech or audio signals, using source filter models or psychoacoustic analysis
    • G10L19/04—Speech or audio signals analysis-synthesis techniques for redundancy reduction, e.g. in vocoders; Coding or decoding of speech or audio signals, using source filter models or psychoacoustic analysis using predictive techniques
    • G10L19/16—Vocoder architecture
    • G10L19/18—Vocoders using multiple modes
    • G10L19/24—Variable rate codecs, e.g. for generating different qualities using a scalable representation such as hierarchical encoding or layered encoding

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Mathematical Physics (AREA)
  • Data Mining & Analysis (AREA)
  • Theoretical Computer Science (AREA)
  • Computational Mathematics (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Pure & Applied Mathematics (AREA)
  • Software Systems (AREA)
  • Databases & Information Systems (AREA)
  • Algebra (AREA)
  • General Engineering & Computer Science (AREA)
  • Multimedia (AREA)
  • Signal Processing (AREA)
  • Discrete Mathematics (AREA)
  • Health & Medical Sciences (AREA)
  • Computational Linguistics (AREA)
  • Spectroscopy & Molecular Physics (AREA)
  • Audiology, Speech & Language Pathology (AREA)
  • Human Computer Interaction (AREA)
  • Acoustics & Sound (AREA)
  • Compression, Expansion, Code Conversion, And Decoders (AREA)
  • Complex Calculations (AREA)
  • Compression Or Coding Systems Of Tv Signals (AREA)
  • Magnetic Resonance Imaging Apparatus (AREA)

Description

関連出願へのクロスレファレンス
本出願は、２００３年９月２９日に出願された米国プロビジョナル特許出願第６０／５０７，２１０号、及び２００３年９月２９に出願された米国プロビジョナル特許出願第６０／５０７，４４０号の優先権の利益を請求するもので、これら各々の内容は、その全体を全ての目的で参考としてここに援用する。
更に、以下の共通に所有する特許出願は、本出願と同日に出願されたもので、その全体をここに援用する。
「Methodfor Transforming a Digital Signal from the Time Domain into the FrequencyDomain and Vice Versa」、代理人ドケット第Ｐ１００４４４号；及び
「Processand Device for Determining a Transforming Element for a Given TransformationFunction, Method and Device for Transforming a Digital Signal from the TimeDomain into the Frequency Domain and Vice Versa and Computer Readable Medium」、代理人ドケット第Ｐ１００４５２号。
背景
本発明は、一般に、デジタル信号処理に係り、より詳細には、デジタルフォーマットされた信号のドメイン変換を実行するシステム及び方法に係る。
ドメイン変換、例えば、離散的コサイン変換（ＤＣＴ）は、近代的な信号処理産業において広く使用されている。近年、ＤＣＴの変形で、整数ＤＣＴと称されるものが、ロスのないコード化アプリケーションにおいて重要な役割を果たすことから、多大な研究上の関心を集めている。「ロスのない」という語は、デコーダが、エンコードされたビットストリームからソース信号の厳密なコピーを発生できることを意味する。
ＤＣＴは、実数値型ブロック変換である。入力ブロックが整数のみで構成されても、ＤＣＴの出力ブロックは、非整数成分を含み得る。便宜上、入力ブロックは、入力ベクトルと称され、出力ブロックは、出力ベクトルと称される。ベクトルが整数成分しか含まない場合には、それが整数ベクトルと称される。ＤＣＴとは対照的に、整数ＤＣＴは、整数入力ベクトルから整数出力ベクトルを発生する。同じ整数入力ベクトルに対して、整数ＤＣＴの整数出力ベクトルは、ＤＣＴの実数出力ベクトルを厳密に近似する。従って、整数ＤＣＴは、スペクトル分析においてＤＣＴの全ての良好なプロパティを保持する。
整数ＤＣＴの重要なプロパティは、可逆性である。可逆性とは、整数逆ＤＣＴ（ＩＤＣＴ）が存在し、整数ＤＣＴが入力ベクトルｘから出力ベクトルｙを発生する場合に、整数ＩＤＣＴがベクトルｙからベクトルｘを回復できることを意味する。時々、整数ＤＣＴは、順方向変換とも称され、一方、整数ＩＤＣＴは、逆方向又は逆変換とも称される。
整数変更型離散的コサイン変換（ＩｎｔＭＤＣＴ）と称される変換が近年提案されて、ＩＳＯ／ＩＥＣ ＭＰＥＧ−４オーディオ圧縮に使用されている。このＩｎｔＭＤＣＴは、そのプロパティ、即ち変更型離散的コサイン変換（ＭＤＣＴ）から導出することができる。「Signal Processing with Lapped Transforms」、Artech House、１９９２年におけるH. S. Malvar著の開示は、ギブンズ(Givens)回転の列をＤＣＴ−ＩＶブロックと共にカスケード構成とすることでＭＤＣＴの効率的な実現を与える。ギブンズ回転は、整数対整数のマッピングを行うための３つのリフトステップへと因数分解できることが良く知られている。例えば、R. Geiger, T. Sporer, J. Koller, K. Brandenburg著の「Audio Coding based onInteger Transforms」、ＡＥＳ第１１１回コンベンション、米国、ニューヨーク、２００１年９月、を参照されたい。整数変換は、各ギブンズ回転を３つのリフトステップに置き換えることによりそれらのプロートタイプから直接コンバートすることができる。各リフトステップには１つの丸め演算があるので、整数変換の全丸め回数は、プロートタイプ変換のギブンズ回転回数の３倍である。離散的三角法変換（例えば、離散的フーリエ変換（ＤＦＴ）又は離散的コサイン変換（ＤＣＴ））の場合に、それに伴うギブンズ回転の回数は、通常、Ｎｌｏｇ２Ｎレベルにあり、ここで、Ｎは、ブロックのサイズであり、即ちデジタル信号が分割される各ブロックに含まれるデータ記号の量である。従って、全丸め回数も、直接コンバートされる整数変換のファミリーに対して、Ｎｌｏｇ２Ｎレベルにある。丸めのために、整数変換は、そのフローティングポイントプロートタイプしか近似しない。丸めの回数と共に近似エラーが増加する。
それ故、より少数の丸め演算を使用してドメイン変換を実行し、これにより、より正確で且つより計算の少ない集中的変換を生じさせるための方法及びシステムの実施が要望される。
発明の概要
本発明は、時間ドメインから周波数ドメインへ及びそれとは逆にデジタル信号のドメイン変換を実行するための方法であって、丸め演算の回数を著しく減少した方法を提供するという課題を解決する。
一実施形態において、時間ドメインから周波数ドメインへ及びそれとは逆にデジタル信号のドメイン変換を実行するための方法は、複数のリフトステージを含む変換エレメントにより変換を実行するステップを備え、変換は、変換マトリクスに対応する。少なくとも１つのリフトステージは、変換マトリクスそれ自体であるか又は低い次元の各変換マトリクスである少なくとも１つの補助的な変換マトリクスを備えている。更に、各リフトステージは、丸めユニットも備えている。この方法は、更に、各補助的な変換マトリクスによる変換の後に各丸めユニットにより信号を処理するステップも備えている。
特定の実施形態の詳細説明
図１２は、本発明の一実施形態によりデジタル信号のドメイン変換を実行するための方法を示す。最初に、１２１０において、複数のリフトステージを備えた変換エレメントによりデジタル信号の変換が実行される。この変換は、変換マトリクスに対応し、リフトステージの少なくとも１つは、補助的な変換マトリクス及び丸めユニットを備え、この補助的なマトリクスは、変換マトリクスそれ自体、又は低い次元の対応変換マトリクスで構成される。ドメイン変換に使用できるプロセス１２１０及びドメイン変換の実施例は、以下に詳細に述べる。
次いで、１２２０において、補助的なマトリクスによる変換の後に丸めユニットにより信号の丸め演算が実行される。
本発明の方法の好ましい実施形態によれば、デジタル信号のデータ記号がデータベクトルとして変換エレメントに与えられる。各リフトステージにおいて、データベクトル又はデータベクトルの一部分がドメイン変換装置により変換され、その変換結果が、次いで、整数ベクトルに丸められる。これは、丸めプロセスがデジタル信号の各ブロックにおいて個々のエレメント又はデータ記号に対して変換プロセス内で実行される従来技術の方法とは対照的である。従って、本発明による方法では丸め演算の回数が著しく減少される。丸め演算の回数が減少されたために、本発明による方法は、大きな計算時間及びコンピュータリソースを必要としない。更に、整数ドメイン変換の近似エラーも著しく減少できる。
好ましい実施形態では、本発明は、整数の形式ＩＶ ＤＣＴの変換を実現するための方法を提供する。本発明による方法は、従来の方法に比して、著しく減少された回数の丸め演算しか必要としない。その結果、近似エラーが著しく減少され、ＤＣＴ−ＩＶの場合に、これを、通常のＮｌｏｇ２Ｎレベルから一実施形態では１．５Ｎほどの低いものに、又、別の実施形態では２．５Ｎほどの低いものに減少することができ、ここで、Ｎは、デジタル信号のブロックサイズである。本発明による方法は、計算の複雑さが低く、且つ構造がモジュール式である。
本発明による方法及び装置は、オーディオ、イメージ又はビデオ信号のような任意の形式のデジタル信号に対して使用できる。物理的な測定信号に対応するデジタル化された信号であるデジタル信号は、対応するアナログ信号の少なくとも１つの特性特徴（例えば、ビデオ信号のルミナンス及びクロミナンス値、アナログサウンド信号の振幅、又はセンサからのアナログ感知信号）を走査することにより発生されてもよい。デジタル信号は、複数のデータ記号を含む。デジタル信号のデータ記号は、ブロックへとグループ編成され、各ブロックは、対応するアナログ信号のサンプリングレートに基づき同じ所定数のデータ記号を有する。
本発明による方法は、整数値を表わす入力デジタル信号を、これも整数値を表わす出力信号へと変換するのに使用できる。本発明による変換方法は、可逆である。出力信号は、本発明による変換方法を実行することにより元の入力信号へと変換して戻すことができる。本発明の方法による変換のこのような可逆プロパティは、出力信号が元の入力信号と同じでなければならないロスのないコード化に使用することができる。
本発明による信号のこのような整数変換は、多数のアプリケーション及びシステム、例えば、ＭＰＥＧオーディオ、イメージ及びビデオ圧縮、ＪＰＥＧ２０００、又はスペクトル分析（赤外線、紫外線又は核磁気放射線信号を分析するための）に使用することができる。又、これは、固定ポイントデジタル信号プロセッサ（ＤＳＰ）のようなハードウェアシステムにおいて、実数値信号変換の場合にオーバーフローのようなファクタを考慮する必要なく、容易に実施することができる。
以下に述べる本発明の好ましい実施形態は、モノ及びステレオの両アプリケーションに適している。モノのアプリケーションにおいては、２つの連続するサンプルブロックがグループ分けされて、一緒に処理される。これは、単一ブロック処理に比して１ブロック長さの信号遅延を導入する。しかしながら、ステレオのアプリケーションでは、左右のチャンネルからの同時サンプルブロックがグループ分けされて一緒に処理された場合に、この余計なブロック遅延を防止することができる。
本発明の方法によれば、デジタル信号は、変換エレメントにより周波数ドメインへ変換される。
好ましくは、変換エレメントは、複数のリフト(lifting)ステージを備えている。
変換エレメントは、リフトラダーのモデルに基づいて説明することができる。リフトラダーモデルは、データ記号の２つのグループの１つを各々受け取る２つのサイドピースを有する。これら２つのサイドピース間には、２つ以上のカスケード状のリフトステージが設けられる。各リフトステージは、一端（入力端）に信号を受信し、加算ユニットを経て他端（出力端）に信号を出力する。丸めユニットは、出力端に配列される。リフトステージは、隣接するリフトステージの出力（又は入力）端が異なるサイドピースに接続されるように、サイドピース間に交互に配列される。
変換エレメントは、リフトラダーモデルの形態で説明するが、これは、変換エレメントの変換経路の例示に過ぎないことに注意されたい。しかしながら、本発明は、このようなラダーモデルに限定されるものではない。
離散的コサイン変換、離散的サイン変換、離散的フーリエ変換、又は離散的Ｗ変換は、本発明の方法により実行できる変換の一例である。
本発明の１つの態様において、各リフトステージは、リフトマトリクスに対応し、これは、２つの可逆の整数マトリクスを１つの対角方向にもつ４つのサブマトリクスで構成されたブロック三角形マトリクスである。リフトマトリクスの一例は、次の通りである。
Figure 0004429316
置換マトリクスとは、別のマトリクスにおいてエレメントの位置を変化させるマトリクスである。又、置換マトリクスは、同一でもよい。又、可逆の整数マトリクスは、例えば、次のような別の仕方でリフトマトリクスに配列されてもよい。
Figure 0004429316
好ましくは、可逆の整数マトリクスは、１又は−１の成分しか含まない対角マトリクスである。
図１は、本発明の実施形態によるオーディオエンコーダ１００の構成を示す図である。このオーディオエンコーダ１００は、変更型離散的コサイン変換（ＭＤＣＴ）に基づく従来の知覚ベースレイヤーコーダと、整数変更型離散的コサイン変換（ＩｎｔＭＤＣＴ）に基づくロスなし改善コーダとを備えている。
例えば、マイクロホン１１０により発生されてアナログ／デジタルコンバータ１１１によりデジタル化されたオーディオ信号１０９が、オーディオエンコーダ１００へ供給される。このオーディオ信号１０９は、複数のデータ記号を含む。このオーディオ信号１０９は、複数のブロックに分割され、各ブロックは、デジタル信号の複数のデータ記号を含み、又、各ブロックは、変更型離散的コサイン変換（ＭＤＣＴ）装置１０１により変換される。ＭＤＣＴ係数は、知覚モジュール１０２の助けで量子化装置１０３により量子化される。知覚モデルは、量子化エラーにより生じる可聴歪が低くなるように量子化装置１０３を制御する。量子化されたＭＤＣＴ係数は、その後、ビットストリームエンコーダ１０４によりエンコードされ、このエンコーダは、ロスのある知覚コード化出力ビットストリーム１１２を発生する。
ビットストリームエンコーダ１０４は、その入力を、ハフマンコード化又はランレングスコード化のような標準的な方法により、ロスのないように圧縮して、その入力より平均ビットレートが低い出力を発生する。又、入力オーディオ信号１０９は、ＩｎｔＭＤＣＴ装置１０５にも供給され、この装置は、ＩｎｔＭＤＣＴ係数を発生する。量子化装置１０３の出力である量子化されたＭＤＣＴ係数は、ＩｎｔＭＤＣＴ係数を予想するのに使用される。量子化されたＭＤＣＴ係数は、逆量子化装置１０６へ供給され、その出力（回復され又は非量子化のＭＤＣＴ係数）は、丸めユニット１０７へ供給される。
丸めユニットは、ＭＤＣＴ係数を整数値へ丸め、整数値ＭＤＣＴとＩｎｔＭＤＣＴ係数との間の差である残余ＩｎｔＭＤＣＴ係数は、エントロピーコーダ１０８によりエントロピーコード化される。このエントロピーエンコーダは、ビットストリームエンコーダ１０４と同様のもので、その入力の平均ビットレートをロスなしに減少して、ロスなしの改善ビットストリーム１１３を発生する。このロスなしの改善ビットストリーム１１３は、知覚コード化ビットストリーム１１２と一緒に、入力オーディオ信号１０９を再構成するに必要な情報を最小限のエラーで搬送する。
図２は、図１に示すオーディオコーダ１００に対応する本発明の実施形態によるオーディオデコーダ２００の構成を示す。知覚コード化ビットストリーム２０７は、ビットストリームデコーダ２０１へ供給され、このデコーダは、図１のビットストリームエンコーダ１０４の動作に対して逆の動作を実行して、デコードされたビットストリームを発生する。デコードされたビットストリームは、逆量子化装置２０２へ供給され、その出力（回復されたＭＤＣＴ係数）は、逆ＭＤＣＴ装置２０３へ供給される。従って、再構成された知覚コード化オーディオ信号２０９が得られる。
ロスなしの改善ビットストリーム２０８は、エントロピーデコーダ２０４へ供給され、このデコーダは、図１のエントロピーエンコーダ１０８の動作とは逆の動作を実行し、対応する残余ＩｎｔＭＤＣＴ係数を発生する。逆量子化装置２０２の出力は、丸め装置２０５によって丸められて、整数値ＭＤＣＴ係数を発生する。この整数値ＭＤＣＴ係数は、残余ＩｎｔＭＤＣＴ係数に加算されて、ＩｎｔＭＤＣＴ係数を発生する。最終的に、逆ＩｎｔＭＤＣＴは、逆ＩｎｔＭＤＣＴ装置２０６によりＩｎｔＭＤＣＴ係数に適用されて、再構成されたロスなしのコード化オーディオ信号２１０を発生する。
図３は、ＤＣＴ−ＩＶを変換として使用すると共に、３つのリフトステージ、即ち第１リフトステージ３０１、第２リフトステージ３０２及び第３リフトステージ３０３を使用する本発明による方法の実施形態を示すフローチャート３００である。この方法は、好ましくは、図１のＩｎｔＭＤＣＴ装置１０５及び図２の逆ＩｎｔＭＤＣＴ装置２０６に使用されて、各々、ＩｎｔＭＤＣＴ及び逆ＩｎｔＭＤＣＴを実施する。図３において、ｘ １ 及びｘ ２ は、各々、デジタル信号の第１及び第２ブロックである。ｚは、中間信号であり、そしてｙ １ 及びｙ ２ は、各々、デジタル信号の第１及び第２ブロックに対応する出力信号である。
上述したように、ＤＣＴ−ＩＶアルゴリズムは、ロスなしのオーディオコード化において重要な役割を演じる。
Figure 0004429316
従って、この実施形態においては、変換エレメントに対応するマトリクスは、次の通りである。
Figure 0004429316
リフトマトリクスの数、ひいては、変換エレメントにおけるリフトステージの数は、ＤＣＴ−ＩＶが変換関数である本発明のこの実施形態では、３である。
Ｎポイント実数入力シーケンスｘ（ｎ）のＤＣＴ−ＩＶは、次のように定義される。
Figure 0004429316
ここで、ｘ １ 、ｘ ２ を２つの整数Ｎｘ１列ベクトルとする。これら列ベクトルｘ １ 、ｘ ２ は、本発明では１つの変換エレメントにより変換されるデジタル信号の２つのブロックに対応する。ｘ １ 、ｘ ２ のＤＣＴ−ＩＶ変換は、各々、ｙ １ 、ｙ ２ である。
Figure 0004429316
（８）及び（９）を合成すると、次のようになる。
Figure 0004429316
前記対角マトリクスは、変換エレメントが対応しているブロック対角マトリクスである。
これは、前記式が次のような簡単な代数変更により変化された場合にも、本発明の範囲内である。
Figure 0004429316
式（１３）は、（７）におけるＤＣＴ−ＩＶプロパティを使用して容易に照合することができる。式（１３）を用いると、式（１１）は、次のように表わすことができる。
Figure 0004429316
式（１４）から、１つの変換エレメントで２つの整数ＤＣＴ−ＩＶを計算する次の整数ＤＣＴ−ＩＶアルゴリズムが導出される。
図４は、ＤＣＴ−ＩＶを変換関数として使用する本発明による方法の実施形態を示す。この実施形態は、ＩｎｔＭＤＣＴを実施するために図１に示すオーディオコーダ１００に使用される。図３と同様に、ｘ １ 及びｘ ２ は、入力デジタル信号の２つのブロックであり、ｚは、中間信号であり、そしてｙ １ 及びｙ ２ は、出力信号の対応ブロックである。
図４に示す３つのリフトステージは、式（１４）における３つのリフトマトリクスに対応する。
図４に示すように、時間対周波数ドメインの整数変換は、次のように決定される。
第１ステージ４０１において、ｘ ２ がＤＣＴ−ＩＶ変換４０２により変換され、次いで、ＤＣＴ−ＩＶ係数が丸められる（４０３）。丸められたＤＣＴ−ＩＶ係数は、次いで、ｘ １ ４０４に加算される。従って、中間信号ｚが発生される。従って、中間信号ｚは、次の式を満足する。
Figure 0004429316
第２ステージ４０５において、ｚがＤＣＴ−ＩＶ変換４０６により変換され、ＤＣＴ−ＩＶ係数が丸められる（４０７）。次いで、丸められたＤＣＴ−ＩＶ係数から、ｘ ２ が減算される。従って、出力信号ｙ１が発生される。従って、出力信号ｙ１は、次の式を満足する。
Figure 0004429316
第３ステージ４０９において、ｙ １ がＤＣＴ−ＩＶ変換４１０により変換され、次いで、ＤＣＴ−ＩＶ係数が丸められる（４１１）。丸められたＤＣＴ−ＩＶ係数が、次いで、ｚから減算される。従って、出力信号ｙ ２ が発生される。従って、出力信号ｙ ２ は、次の式を満足する。
Figure 0004429316
図５は、ＤＣＴ−ＩＶを変換関数として使用する本発明の方法の実施形態に基づいて逆変換を行うためのアルゴリズムを示す。この実施形態は、逆ＩｎｔＭＤＣＴを実施するために図２に示すオーディオデコーダ２００に使用される。図５に示すアルゴリズムは、図４に示すアルゴリズムの逆である。異なる信号ｙ １ 、ｙ ２ 、ｘ １ 、ｘ ２ 、ｚに対する表示は、図４の表示に対応して選択されたものである。
図５に示すように、周波数対時間ドメインの整数変換は、次のように決定される。
第１ステージ５０１において、ｙ １ がＤＣＴ−ＩＶ変換５０２により変換され、ＤＣＴ−ＩＶ係数が丸められる（５０３）。丸められたＤＣＴ−ＩＶ係数は、次いで、ｙ ２ に加算される（５０４）。従って、中間信号ｚが発生される。従って、中間信号ｚは、次の式を満足する。
Figure 0004429316
第２ステージ５０５において、ｚがＤＣＴ−ＩＶ変換５０６により変換され、ＤＣＴ−ＩＶ係数が丸められる。次いで、その丸められたＤＣＴ−ＩＶ係数から、ｙ １ が減算される。従って、信号ｘ ２ が発生される。従って、信号ｘ ２ は、次の式を満足する。
Figure 0004429316
第３ステージ５０９において、ｘ ２ がＤＣＴ−ＩＶ変換５１０により変換され、ＤＣＴ−ＩＶ係数が丸められる（５１１）。この丸められたＤＣＴ−ＩＶ係数が、次いで、ｚから減算される。従って、信号ｘ １ が発生される。従って、ｘ １ は、次の式を満足する。
Figure 0004429316
式（１６ａ）から（１６ｃ）に基づくアルゴリズムは、式（１５ａ）から（１５ｃ）に基づくアルゴリズムに対して逆であることが明らかである。従って、図１及び２に示すエンコーダ及びデコーダに使用した場合には、これらアルゴリズムは、ロスのないオーディオコード化のための方法及び装置をもたらす。
式（１５ａ）から（１５ｃ）及び（１６ａ）から（１６ｃ）は、更に、２つのＮｘＮ整数ＤＣＴ−ＩＶを比較するために、３つのＮｘＮＤＣＴ−ＩＶ、３つのＮｘ１丸め、及び３つのＮｘ１加算が必要であることを示している。それ故、１つのＮｘＮ整数ＤＣＴ−ＩＶについては、平均が次のようになる。
Figure 0004429316
但し、ＲＣ（．）は、全丸め回数であり、ＡＣ（．）は、全演算回数である。直接コンバートされる整数ＤＣＴ−ＩＶアルゴリズムに比して、ここに提案する整数ＤＣＴ−ＩＶアルゴリズムは、ＲＣをレベルＮｌｏｇ２ＮからＮへ減少する。
式（１８）で示されたように、ここに提案する整数ＤＣＴ−ＩＶアルゴリズムの演算の複雑さは、ＤＣＴ−ＩＶアルゴリズムより約５０％も高い。しかしながら、ＲＣも考慮した場合には、ここに提案するアルゴリズムの合成複雑さ（ＡＣ＋ＲＣ）は、直接コンバート式の整数アルゴリズムを著しく越えることはない。アルゴリズムの複雑さの厳密な分析は、使用するＤＣＴ−ＩＶアルゴリズムに依存する。
図４及び５に示すように、ここに提案する整数ＤＣＴ−ＩＶアルゴリズムは、構造が簡単でモジュール状である。これは、何らかの既存のＤＣＴ−ＩＶアルゴリズムをそのＤＣＴ−ＩＶ計算ブロックに使用することができる。ここに提案するアルゴリズムは、ＩｎｔＭＤＣＴを要求するアプリケーション、例えば、ＭＰＥＧ−４オーディオ拡張３参照モデル０に適している。
図６は、５つのリフトステージ、即ち第１リフトステージ６０１、第２リフトステージ６０２、第３リフトステージ６０３、第４リフトステージ６０４、及び第５リフトステージ６０５を使用する本発明による方法の実施形態を示すフローチャート６００である。この方法は、図１のＩｎｔＭＤＣＴ装置１０５及び図２の逆ＩｎｔＭＤＣＴ装置２０６に使用して、各々、ＩｎｔＭＤＣＴ及び逆ＩｎｔＭＤＣＴを実施することができる。図６において、ｘ １ 及びｘ ２ は、各々、デジタル信号の第１及び第２ブロックである。ｚ １ 、ｚ ２ 、ｚ ３ は、中間信号であり、そしてｙ １ 及びｙ ２ は、各々、デジタル信号の第１及び第２ブロックに対応する出力信号である。
図７は、変換関数がＤＣＴ−ＩＶ変換関数である本発明による方法の実施形態を示すフローチャートである。
Figure 0004429316
Figure 0004429316
式（２０）は、更に、次のように書き表すことができる。
Figure 0004429316
変換エレメントは、式（２９）の５つのリフトマトリクスに対応する５つのリフトステージを備えている。
更に、変換エレメントは、置換マトリクスＰ co に対応するデータシャフルステージも備えている。
図７において、第１リフトステージの入力は、デジタル信号ｘ １ 及びｘ ２ の２つのブロックであり、ｚ １ 、ｚ ２ 及びｚ ３ は、中間信号であり、そしてｙ １ 及びｙ ２ は、各々、デジタル信号の第１及び第２ブロックに対応する出力信号である。
変換エレメントｘの入力、及び変換エレメントｘ １ 及びｘ ２ の第１リフトステージの２つの入力ブロックは、次の式を満足する。
Figure 0004429316
以下、リフトマトリクスＴ ３ に対応するリフトステージである第１リフトステージ７０１について説明する。
Figure 0004429316
式（２２）により与えられるＴ ３ の定義を使用すると、式（２９）は、次のように書き直される。
Figure 0004429316
この実施形態では、整数ＤＣＴ−ＩＶに対する可逆アルゴリズムが与えられるので、整数値への丸めが含まれる。従って、式（３０）に基づき、第１リフトステージ７０１の第１ステップ７０６において、ｘ １ にＫ ３ が乗算される。この乗算の結果が、ステップ７０７において、整数値に丸められる。丸められた値は、次いで、ステップ７０８において、ｘ ２ に加算される。従って、中間信号ｚ １ は、次の式を満足する。
Figure 0004429316
図７に示された変換エレメントの他の４つのリフトステージ７０２、７０３、７０４、７０５であって、マトリクスＴ ２ 、Ｓ、Ｒ ２ 及びＲ １ に各々対応するリフトステージは、第１リフトステージ７０１と同一の構造を有するので、それらの説明は省略する。第２のリフトステップ７０２の加算ステップ７０９では、ｘ １ が、Ｔ ２ の定義に基づき、−Ｄ N/2 で乗算されることだけ注意されたい。
以下、逆変換の変換エレメントのリフトステージについて、図８を参照して説明する。
図８は、図７に示す変換の逆変換を行う変換エレメントのリフトステージを示す。
図８において、第１リフトステージの入力は、デジタル信号ｙ １ 及びｙ ２ の２つのブロックであり、ｚ １ 、ｚ ２ 、ｚ ３ は、中間信号であり、そしてｘ １ 及びｘ ２ は、各々、デジタル信号の第１及び第２ブロックに対応する出力信号である。
図８に示す最後のリフトステージ８０５は、図７に示す第１のリフトステージ７０１に対して逆である。従って、最後のリフトステージ８０５の第１ステップ８０６において、ｘ １ がＫ ３ で乗算される。この乗算の結果がステップ８０７において整数値へ丸められる。丸められた値は、ステップ８０８において、ｚ １ から減算される。従って、信号ｘ ２ は、次の式を満足する。
Figure 0004429316
図８に示された変換エレメントの他の４つのリフトステージ８０１、８０２、８０３、８０４であって、図７に示された変換エレメントのリフトステージ７０５、７０４、７０３及び７０２に対して逆のリフトステージは、最後のリフトステージ８０５と同一の構造を有するので、それらの説明は省略する。第４のリフトステップ８０４の加算ステップ８０９の後に、この加算ステップ８０９の結果が−Ｄ N/2 で乗算されて、ｘ １ を生じることだけ注意されたい。
図８のリフトステージ８０５、８０４、８０３、８０２及び８０１は、各々、図７のリフトステージ７０１から７０５に対して逆であることが明らかであろう。又、マトリクスＰ co に対応する入力信号の順列を逆にできると共に、それに応じたデータシャフルステージが逆の変換エレメントにより構成されるので、個々に提供される方法は可逆であり、従って、図１及び２に示すオーディオエンコーダ１００及びオーディオデコーダ２００に使用された場合には、ロスのないオーディオコード化を行う方法及び装置が提供される。
この実施形態に使用される丸めの回数の分析は、本発明の説明の終りで行う。
図９は、本発明の実施形態によるイメージアーカイブシステムの構成を示す。
図９において、イメージソース９０１、例えば、カメラは、アナログイメージ信号を発生する。このイメージ信号は、アナログ／デジタルコンバータ９０２により処理されて、それに対応するデジタルイメージ信号を発生する。このデジタルイメージ信号は、時間ドメインから周波数ドメインへの変換を含むロスなしのイメージエンコーダ９０３によりロスなしにエンコードされる。この実施形態では、時間ドメインは、イメージの座標スペースに対応する。ロスなしのコード化イメージ信号が、記憶装置９０４、例えば、ハードディスク又はＤＶＤに記憶される。イメージが必要とされるときには、ロスなしのコード化イメージ信号が、記憶装置９０４からフェッチされて、ロスなしのイメージデコーダ９０５へ送られ、このデコーダは、ロスなしのコード化イメージ信号をデコードして、データロスのない元のイメージ信号を再構成する。
イメージ信号のこのようなロスなしのアーカイブは、例えば、イメージが半導体ウェハのエラーマップであって、後で分析するためにこれを記憶しなければならない場合に、重要である。
以下、本発明により時間ドメインから周波数ドメインへ及びそれとは逆にデジタル信号を変換するための方法の更に別の実施形態を説明する。この実施形態は、図９に示すイメージアーカイブシステムのロスなしのイメージエンコーダ９０３及びロスなしのイメージデコーダ９０５に使用されるのが好ましい。
図１０は、ＤＷＴ−ＩＶを変換関数として使用する本発明による方法の実施形態を示す。
Ｎポイント実数入力シーケンスｘ（ｎ）のＤＷＴ−ＩＶは、次のように定義される。
Figure 0004429316
ＤＷＴ−ＩＶマトリクスは、次の式へと因数分解される。
Figure 0004429316
Ｉ Ｎ／２ は、次数Ｎ／２の恒等マトリクスである。Ｊ Ｎ／２ は、次数Ｎ／２の対抗恒等マトリクスであり、即ち、次のようになる。
Figure 0004429316
Ｔは、次式で与えられるＮｘＮマトリクスである。
Figure 0004429316
Ｄ Ｎ／２ は、次式で与えられる次数Ｎ／２の対角マトリクスである。
Figure 0004429316
Ｒ Ｎ 及びＴは、リフトマトリクスの積へと更に因数分解することができる。
Figure 0004429316
従って、式（３８）は、次の式で書き表すことができる。
Figure 0004429316
リフトマトリクスＲ ３ 及びＴ １ は、リフトマトリクスＳへと集合させることができる。
Figure 0004429316
式（４９）及び（５０）から、ＤＷＴ−ＩＶマトリクスに対する次の因数分解式を得ることができる。
Figure 0004429316
更に、変換エレメントは、置換マトリクスＰ Ｎ に対応するデータシャフルステージを備えている。このデータシャフルステージは、各入力データブロックにおける成分の順序を再配列する。Ｐ Ｎ に基づいて、入力データベクトルは、次のように再配列される。即ち、ベクトルの第１の半分は、不変のままとし、ベクトルの第２の半分は、ひっくり返す。即ち、次のようになる。
Figure 0004429316
図１０において、第１リフトステージの入力は、デジタル信号ｘ １ 及びｘ ２ の２つのブロックであり、ｚ １ 、ｚ ２ 及びｚ ３ は、中間信号であり、そしてｙ １ 及びｙ ２ は、各々、デジタル信号の第１及び第２ブロックに対応する出力信号である。
変換エレメントｘの入力、及び変換エレメントｘ １ 及びｘ ２ の第１リフトステージの２つの入力ブロックは、次の式を満足する。
Figure 0004429316
以下、リフトマトリクスＴ ３ に対応するリフトステージである第１リフトステージ１００１について説明する。
Figure 0004429316
式（４５）により与えられるＴ ３ の定義を使用すると、式（５１）は、次のように書き直される。
Figure 0004429316
この実施形態では、整数ＤＷＴ−ＩＶに対する可逆アルゴリズムが与えられるので、整数値への丸めが含まれる。従って、式（５５）に基づき、第１リフトステージ１００１の第１ステップ１００６において、ｘ １ にＫ ３ が乗算される。この乗算の結果が、ステップ１００７において、整数値に丸められる。丸められた値は、次いで、ステップ１００８において、ｘ ２ に加算される。従って、中間信号ｚ １ は、次の式を満足する。
Figure 0004429316
図１０に示された変換エレメントの他の４つのリフトステージ１００２、１００３、１００４、１００５であって、マトリクスＴ ２ 、Ｓ、Ｒ ２ 及びＲ １ に各々対応するリフトステージは、第１リフトステージ１００１と同一の構造を有するので、それらの説明は省略する。第２のリフトステップ１００２の加算ステップ１００９では、ｘ １ が、Ｔ ２ の定義に基づき、Ｄ N/2 で乗算されることだけ注意されたい。
以下、逆変換の変換エレメントのリフトステージについて図１１を参照して説明する。
図１１は、図１０に示す変換の逆変換を行う変換エレメントのリフトステージを示す。
図１１において、第１リフトステージの入力は、デジタル信号ｙ １ 及びｙ ２ の２つのブロックであり、ｚ １ 、ｚ ２ 、ｚ ３ は、中間信号であり、そしてｘ １ 及びｘ ２ は、各々、デジタル信号の第１及び第２ブロックに対応する出力信号である。
図１１に示す最後のリフトステージ１１０５は、図１０に示す第１のリフトステージ１００１に対して逆である。従って、最後のリフトステージ１１０５の第１ステップ１１０６において、ｘ １ がＫ ３ で乗算される。この乗算の結果がステップ１１０７において整数値へ丸められる。丸められた値は、ステップ１１０８において、ｚ １ から減算される。従って、信号ｘ ２ は、次の式を満足する。
Figure 0004429316
図１１に示された変換エレメントの他の４つのリフトステージ１１０１、１１０２、１１０３、１１０４であって、リフトステージ１００５、１００４、１００３及び１００２に対して逆のリフトステージは、最後のリフトステージ１１０５と同一の構造を有するので、それらの説明は省略する。第４のリフトステップ１１０４の加算ステップ１１０９の後に、この加算ステップ１１０９の結果がＤ N/2 で乗算されて、ｘ １ を生じることだけ注意されたい。
図１１のリフトステージ１１０５、１１０４、１１０３、１１０２及び１１０１は、各々、図１０のリフトステージ１００１から１００５に対して逆であることが明らかであろう。又、マトリクスＰ Ｎ に対応する入力信号の順列を逆にできると共に、それに応じたデータシャフルステージが逆の変換エレメントにより構成されるので、全体的に提供される方法は可逆であり、従って、図９に示すロスなしのイメージエンコーダ９０３及びロスなしのイメージデコーダ９０５に使用された場合には、ロスのないイメージコード化を行う方法及び装置が提供される。
ここに説明する実施形態では、ＤＣＴ−ＩＶのための本発明による方法の実施形態がオーディオコード化に使用されると共に、ＤＷＴ−ＩＶのための本発明による方法の実施形態がイメージコード化に使用されたが、ＤＣＴ−ＩＶのための本発明による方法の実施形態をイメージコード化に使用することもでき、又、ＤＷＴ−ＩＶのための本発明による方法の実施形態をオーディオコード化に使用することもでき、更に、その全てをビデオ信号のような他のデジタル信号のコード化に使用することもできる。
式（３４）及び（３５）について考えると、各リフトステージにはＮ／２回の丸めがあることが明らかである。それ故、式（２６）について考えると、本発明による図７及び８に示すＤＣＴ−ＩＶアルゴリズムの実施形態の変換エレメントに対する合計丸め回数は、Ｎ／２の５倍、即ち２．５Ｎであり、これは、現状技術によるＮｌｏｇ２Ｎより著しく低い。
Figure 0004429316
整数ＤＷＴ−ＩＶ変換関数についても同様の結論を引き出すことができる。
以下、本発明により時間ドメインから周波数ドメインへ及びそれとは逆にデジタル信号を変換するための方法の更に別の実施形態を説明する。
本発明のこの実施形態では、ドメイン変換がＤＣＴ変換であり、従って、ブロックサイズＮはある整数であると共に、入力ベクトルが２つのサブベクトルを含む。
Figure 0004429316
Figure 0004429316
各リフトマトリクスＬ２Ｎに対して、リフトステージの可逆の整数対整数マッピングは、ここに援用する参照文献「Factoring Wavelet Transforms into Lifting Steps」、Tech. Report、I. Daubechies及びW. Sweldens著、Bell Laboratories, LucentTechnologies、１９９６年、に説明された２ｘ２リフトステップと同様に実現される。唯一の相違は、単一の変数ではなく、ベクトルに丸め演算が適用されることである。
他の実施形態の前記説明において、リフトマトリクスに対してリフトステージをいかに実現するかは既に述べた。従って、リフトマトリクスに対応するリフトステージの説明は、以下では省略する。
Figure 0004429316
前記式の右側を構成しているマトリクスは、次のように説明される。
Ｐ１は、次式で与えられる第１の置換マトリクスである。
Figure 0004429316
Ｐ２は、第２の置換マトリクスであり、その一例が、次のＭＡＴＬＡＢスクリプトにより発生される。
Figure 0004429316
一例として、Ｎが４であるときには、Ｐ２が次のような８ｘ８マトリクスである。
Figure 0004429316
Ｐ３は、第３の置換マトリクスであり、その一例が、次のＭＡＴＬＡＢスクリプトにより発生される。
Figure 0004429316
一例として、Ｎが４であるときには、Ｐ３が次のような８ｘ８マトリクスである。
Figure 0004429316
Figure 0004429316
Figure 0004429316
Figure 0004429316
Figure 0004429316
Figure 0004429316
Figure 0004429316
Figure 0004429316
Figure 0004429316
Figure 0004429316
式（８４）から、次元Ｎｘ１の２つの入力信号に対して整数ＤＣＴを計算することができる。
式（８４）は、ＤＣＴ−ＩＶ変換ドメインを記述するリフトマトリクス因数分解を与えるので、そのリフトマトリクスを、ここに示すように使用して、与えられた入力信号のドメイン変換を計算することができる。
式（８４）は、次のようにして導出することができる。
次の分解は、Wang、Zhongde著の「On Computing the Discrete Fourier and Cosine Transforms」、IEEE Transactions onAcoustics, Speech and Signal Processing、第ＡＳＳＰ−３３巻、第４号、１９８５年１０月、からの開示を使用して導出することができる。
Figure 0004429316
Figure 0004429316
Figure 0004429316
Figure 0004429316
Figure 0004429316
この実施形態では、ドメイン変換を計算するときに、４Ｎ回の丸め演算しか必要とされず、これについて、以下に述べる。
ここで、各々、α（＊）を実数加算の回数とし、μ（＊）を実数乗算の回数とし、γ（＊）を実数丸めの回数とする。ここに提案するＩｎｔＤＣＴアルゴリズムの場合、次のようになる。
α（ＩｎｔＤＣＴ）＝１１Ｎ＋３α（ＤＣＴ−ＩＶ）
μ（ＩｎｔＤＣＴ）＝９Ｎ＋３μ（ＤＣＴ−ＩＶ）
γ（ＩｎｔＤＣＴ）＝８Ｎ
前記結果は、データサンプルの２つのブロックに対するものである。というのは、ここに提案するＩｎｔＤＣＴアルゴリズムは、それらを一緒に処理するからである。従って、データサンプルの１つのブロックに対して、計算の回数が半減し、次のようになる。
α１（ＩｎｔＤＣＴ）＝５．５Ｎ＋１．５α（ＤＣＴ−ＩＶ）
μ１（ＩｎｔＤＣＴ）＝４．５Ｎ＋１．５μ（ＤＣＴ−ＩＶ）
γ１（ＩｎｔＤＣＴ）＝４Ｎ
但し、α１、μ１及びγ１は、各々、１つのサンプルブロックに対する実数加算の回数、実数乗算の回数、及び実数丸めの回数である。
ＤＣＴ−ＩＶ計算の場合には、ここに援用する参照文献「Signal Processing with lapped Transforms」、H. S. Malvar著、マサチューセッツ州ノルウッド、アーテック・ハウス、１９９２年、第１９９−２０１ページ、に掲載されたＦＦＴベースのアルゴリズムを使用することができ、これについては、
α（ＤＣＴ−ＩＶ）＝１．５Ｎｌｏｇ２Ｎ
μ（ＤＣＴ−ＩＶ）＝０．５Ｎｌｏｇ２Ｎ＋Ｎ
その結果、次のようになる。
α１（ＩｎｔＤＣＴ）＝２．２５Ｎｌｏｇ２Ｎ＋５．５Ｎ
μ１（ＩｎｔＤＣＴ）＝０．７５Ｎｌｏｇ２Ｎ＋６Ｎ
以下、時間ドメインから周波数ドメインへ及びそれとは逆にデジタル信号を変換する方法の更に別の実施形態を説明する。
この実施形態では、離散的高速フーリエ変換（ＦＦＴ）がドメイン変換として使用される。
Ｆを、正規化されたＦＦＴのＮｘＮ変換マトリクスとする。
Figure 0004429316
但し、Ｎは、変換サイズで、ある正の整数であり、一実施形態では、２の累乗、即ちＮ＝２ｉであり、そしてｉは、ゼロより大きな整数である。ｍ及びｎは、マトリクス指数である。
この実施形態のもとでは、次元ＮｘＮの置換マトリクスＰは、指数０又は１を含むマトリクスである。これにＮｘ１ベクトル（入力信号のマトリクス表示）を乗算した後、ベクトルにおけるエレメントの順序が変更される。
この実施形態では、リフトマトリクスは、次の式の２Ｎｘ２Ｎマトリクスとして定義される。
Figure 0004429316
但し、Ｐ１及びＰ２は、２つの置換マトリクスであり、Ｏは、ＮｘＮのゼロマトリクスであり、Ａは、任意のＮｘＮマトリクスである。リフトマトリクスＬについては、Ｉ．ダウベッチーズ氏のここに援用する前記参照文献における２ｘ２リフトステップと同様に実現される。しかしながら、上述したように、単一の変数ではなく、丸めがベクトルに適用される。Ｌの転置、ＬＴも、リフトマトリクスであることが明らかであろう。
更に、Ｔを２Ｎｘ２Ｎ変換マトリクスとする。
Figure 0004429316
従って、変更型変換マトリクスＴ（従って、ドメイン変換それ自体）は、リフトマトリクス因数分解として表現することができる。
Figure 0004429316
式（９６）において、かぎカッコ内のブランクスペースは、全てゼロのマトリクスエレメントを表わす。
式（９４）の右辺のリフトマトリクスは、この場合に変換マトリクスＦそれ自体である補助的な変換マトリクスを備えている。
式（９４）から明らかなように、リフトマトリクスの因数分解は、ここに述べる方法を使用して２つのＮｘ１複素数ベクトルに対する整数ＦＦＴを計算するのに使用できる。
この実施形態のもとでは、ドメイン変換の計算は、３Ｎ回の丸め演算しか必要とせず、これについて、以下に述べる。
α（＊）を実数加算の回数とする。
μ（＊）を実数乗算の回数とする。
γ（＊）を実数丸め演算の回数とする。
ここに提案するＩｎｔＦＦＴアルゴリズムの場合に、次のものが得られる。
α（ＩｎｔＦＦＴ）＝６Ｎ＋３α（ＦＦＴ）
μ（ＩｎｔＦＦＴ）＝３μ（ＦＦＴ）
γ（ＩｎｔＦＦＴ）＝６Ｎ
前記結果は、データサンプルの２つのブロックに対するものである。というのは、ここに提案するＩｎｔＦＦＴアルゴリズムは、それらを一緒に処理するからである。従って、データサンプルの１つのブロックに対して、計算の回数が半減し、次のようになる。
α１（ＩｎｔＦＦＴ）＝３Ｎ＋１．５α（ＦＦＴ）
μ１（ＩｎｔＦＦＴ）＝１．５μ（ＦＦＴ）
γ１（ＩｎｔＦＦＴ）＝３Ｎ
但し、α１、μ１及びγ１は、各々、１つのサンプルブロックに対する実数加算の回数、実数乗算の回数、及び実数丸め演算の回数である。
ＦＦＴ計算に対して、スプリット−ラディックスＦＦＴ（ＳＲＦＦＴ）アルゴリズムを使用することができ、これについては、
α（ＳＲＦＦＴ）＝３Ｎｌｏｇ２Ｎ−３Ｎ＋４
μ（ＳＲＦＦＴ）＝Ｎｌｏｇ２Ｎ−３Ｎ＋４
その結果、次のようになる。
α１（ＩｎｔＦＦＴ）＝４．５Ｎｌｏｇ２Ｎ−１．５Ｎ＋６
μ１（ＩｎｔＦＦＴ）＝１．５Ｎｌｏｇ２Ｎ−４．５Ｎ＋６
図１３は、上述したＤＣＴ変換技術及び上述したＦＦＴドメイン変換の変換精度を評価するのに使用される順方向及び逆方向変換コーダを示す。テストは、ここに援用する「Coding of Moving Pictures and Audio: Work plan for Evaluation ofInteger MDCT for FGS to Lossless Experimentation Framework」、ＩＳＯ／ＩＥＣ ＪＴＣ １／ＳＣ ２９／ＷＧ １１ Ｎ５５７８、Pattaya、Thailand、２００３年３月、に記載されたＭＰＥＧ−４ロスなしオーディオコードグループにより提案された評価規格に基づいて変換の平均２乗エラー（ＭＳＥ）を測定することを含む。
より詳細には、ＩｎｔＤＣＴ及び整数逆ＤＣＴ（ＩｎｔＩＤＣＴ）に対するＭＳＥは、次のように表わされる。
Figure 0004429316
但し、エラー信号ｅは、図１の場合と同様に、ＩｎｔＤＣＴについてはｅｆ、及びＩｎｔＩＤＣＴについてはｅｔである。Ｋは、評価に使用されるサンプルブロックの合計数である。
ＩｎｔＦＦＴ及び整数逆ＦＦＴ（ＩｎｔＩＦＦＴ）に対するＭＳＥは、次のように表わされる。
Figure 0004429316
両ドメイン変換の場合に、１５個の異なるタイプの音楽ファイルで合計４５０秒が４８ｋＨｚ／１６ビットテストセットに使用された。表１は、テスト結果を示す。
表１から明らかなように、本発明のシステム及び方法を使用して生じたＭＳＥは、非常に僅かであり、そして従来のシステムとは異なり、処理部とっくサイズとは実質的に独立したものである。ＤＣＴ−ＩＶドメイン変換を参照すれば、ＭＳＥは、４０９６ビットまでのブロックザイズＮの増加と共に僅かに増加するだけである。ＦＦＴのＭＳＥは、いっそう良好で、４０９６ビットまでのブロックサイズに対して０．４の一定ＭＳＥを示す。本発明の実証された性能を、より長いブロックサイズに対する現在の能力及び増加するニーズの観点から見ると、本発明の効果が明らかとなろう。
Figure 0004429316
援用する参照文献
以下の文書を、参考としてここに援用する。
H.S. Malvar著、「SignalProcessing with Lapped Transforms」、Artech House、１９９２年；
R.Geiger, T. Sporer, J. Koller, K. Brandenburg著、「Audio Coding based on Integer Transforms」、ＡＥＳ第１１１回コンベンション、米国、ニューヨーク、２００１年９月；
Wang,Zhongde著、「OnComputing the Discrete Fourier and Cosine Transforms」、IEEE Transactions onAcoustics, Speech and Signal Processing、第ＡＳＳＰ−３３巻、第４号、１９８５年１０月；
I.Daubechies及びW.Sweldens著、「Factoringwavelet transforms into lifting steps」、Tech. Report, Bell Laboratories, Lucent Technologies、１９９６年；
S.Oraintara, Y. J. Chen及びT. Q. Nguyen著、「Integer fast Fourier transform」、IEEE Trans. Signal Processing、第５０巻、第３号、２００２年３月、第６０７−６１８ページ；
P.Hao及びQ. Shi著、「Matrix factorization forreversible integer mapping」、IEEE Trans. Signal Processing、第４９巻、第１０号、２００１年１０月、第２３１４−２３２４ページ；
G.Plonka及びM.Tasch著、「Invertible integer DCTalgorithms」、Appl.Comput. Harmon. Anal. 12: 70-88、２００３年；
Y.H. Zeng、L. Z. Cheng、G. A. Bi及びAlex C. Kot著、「Integer DCTs and fastalgorithms」、IEEETrans. Signal Processing、第４９巻、第１１号、２００１年１１月、第２７７４−２７８２ページ；
J.Wang, J. Sun及びS.Yu著、「1-D and2-D transforms from integers to integers」、in Proc. Int. Conf. Acoustics, Speech and Signal Processing、香港、２００３年、第ＩＩ巻、第５４９−５５２ページ；及び
「Codingof Moving Pictures and Audio: Work plan for Evaluation of Integer MDCT for FGSto Lossless Experimentation Framework」、ISO/IEC JTC 1/SC 29/WG 11 N5578, Pattaya, Thailand、２００３年３月。
本発明の実施形態によるオーディオエンコーダの構成を示す図である。 図１に示すオーディオコーダに対応する、本発明の実施形態によるオーディオデコーダの構成を示す図である。 本発明による方法の実施形態を示すフローチャートである。 ＤＣＴ−ＩＶを変換関数として使用する本発明による方法の実施形態を示す図である。 図４に示す本発明の方法の実施形態に基づく逆変換のアルゴリズムを示す図である。 本発明による方法の実施形態を示すフローチャートである。 ＤＣＴ−ＩＶを変換関数として使用する本発明による方法の実施形態を示す図である。 図７に示す本発明の方法の実施形態に基づく逆変換のアルゴリズムを示す図である。 本発明の実施形態によるイメージアーカイブシステムの構成を示す図である。 ＤＷＴ−ＩＶを変換関数として使用する本発明による方法の実施形態を示す図である。 図１０に示す本発明の方法の実施形態に基づく逆変換のアルゴリズムを示す図である。 本発明の実施形態によるデジタル信号のドメイン変換を実行するための方法を示す図である。 本発明により与えられるＤＣＴ及びＥＦＴ変換方法の変換精度を評価するのに使用される順方向及び逆方向変換コーダを示す図である。
符号の説明
１００・・・オーディオエンコーダ、１０１・・・変更型離散的コサイン変換（ＭＤＣＴ）、１０２・・・知覚モジュール、１０３・・・量子化装置、１０４・・・ビットストリームエンコーダ、１０５・・・ＩｎｔＭＤＣＴ装置、１０６・・・逆量子化装置、１０７・・・丸めユニット、１０８・・・エントロピーコーダ、１０９・・・オーディオ信号、１１０・・・マイクロホン、１１１・・・アナログ／デジタルコンバータ、１１２・・・知覚コード化出力ビットストリーム、１１３・・・ロスなしの改善ビットストリーム、２００・・・オーディオでコーダ、２０１・・・ビットストリームデコーダ、２０２・・・逆量子化装置、２０３・・・逆ＭＤＣＴ装置、２０４・・・エントロピーデコーダ、２０５・・・丸め装置、２０６・・・逆ＩｎｔＭＤＣＴ装置、２０７・・・知覚コード化ビットストリーム、２０８・・・ロスなしの改善ビットストリーム、２０９・・・再構成された知覚コード化オーディオ信号、３０１・・・第１リフトステージ、３０２・・・第２リフトステージ、３０３・・・第３リフトステージ、４０１・・・第１ステージ、４０２・・・ＤＣＴ−ＩＶ変換、４０３・・・丸め、４０４・・・加算、４０５・・・第２ステージ、４０６・・・ＤＣＴ−ＩＶ変換、４０７・・・丸め、４０９・・・第３ステージ、４１０・・・ＤＣＴ−ＩＶ変換、４１１・・・丸め、５０１・・・第１ステージ、５０２・・・ＤＣＴ−ＩＶ変換、５０３・・・丸め、５０４・・・加算、５０５・・・第２ステージ、５０６・・・ＤＣＴ−ＩＶ変換、５０７・・・丸め、５０９・・・第３ステージ、５１０・・・ＤＣＴ−ＩＶ変換、５１１・・・丸め

Claims (11)

  1. 時間ドメインから周波数ドメインへ及びそれとは逆にデジタル信号のドメイン変換を実行する装置において、
    前記デジタル信号に対する前記ドメイン変換を実行するための変換エレメントを有する変換ユニットを備え、
    前記デジタル信号のデータ記号は前記変換エレメントにデータベクトルとして提供され、前記変換エレメントは前記ドメイン変換を複数のリフトステージで実行し、各リフトステージは前記データベクトルの少なくとも一部に対する前記ドメイン変換の一ステップを実行し、前記ドメイン変換は離散コサイン変換（ＤＣＴ）−ＩＶの変換マトリクスに対応し、該変換マトリクスに従って前記データベクトルが変換され、
    各リフトステージは、少なくとも１つの補助的な変換マトリクスを備え、該補助的な変換マトリクスは、前記変換マトリクスそれ自体であるか、又は低い次元の各変換マトリクスであり、各リフトステージにおいて前記データベクトルの少なくとも一部が前記少なくとも１つの補助的な変換マトリクスによって変換され、前記データベクトルの前記少なくとも一部はリフトステージにおいて前記少なくとも１つの補助的な変換マトリクスとの乗算に基づいて処理され、更に、各リフトステージは、丸め演算も含み、
    これにより、前記データベクトル又は前記データベクトルの前記一部が各リフトステージにおいて前記各補助的な変換マトリクスによる前記ドメイン変換の後に前記各丸め演算により処理され、
    前記変換ユニットは、
    複数のデータブロックを受け取るように結合されると共に、各データブロックをＭＤＣＴ係数へとドメイン変換するように構成された変更型離散コサイン変換装置と、
    前記ＭＤＣＴ係数の各々を受け取るように結合された量子化装置であって、それに応答して、量子化されたＭＤＣＴ係数を発生するように動作できる量子化装置と、
    前記量子化されたＭＤＣＴ係数を受け取るように結合されたビットストリームエンコーダであって、それに応答して、知覚的にコード化されたビットストリームを発生するビットストリームエンコーダと、
    前記量子化されたＭＤＣＴ係数を受け取るように結合された逆量子化装置であって、前記ＭＤＣＴ係数を非量子化状態に回復させるように動作できる逆量子化装置と、
    前記回復されたＭＤＣＴ係数を受け取るように結合されると共に、整数値のＭＤＣＴ係数を発生するように動作できる丸めユニットと、
    を備える、
    装置。
  2. 前記変換ユニットは、更に、
    前記データブロックを受け取るように結合されると共に、それに応答して、ＩｎｔＭＤＣＴ係数を発生するように動作できる逆変更型離散的コサイン変換装置と、
    各ＩｎｔＭＤＣＴ係数と整数値のＭＤＣＴ係数との間の差を計算して、各残余ＭＤＣＴ係数を発生するための手段と、
    前記残余ＭＤＣＴ係数を受け取るように結合されると共に、それに応答して、ロスなしの改善ビットストリームを発生するように動作できるエントロピーコーダと、
    を備える請求項１に記載の装置。
  3. 前記変換ユニットは、更に、
    前記知覚コード化ビットストリームを受け取るように結合されると共に、それに応答して、デコードされたビットストリームを出力するように動作できるビットストリームデコーダと、
    前記デコードされたビットストリームを受け取るように結合され、それに応答して、回復されたＭＤＣＴ係数を発生するための逆量子化装置と、
    前記回復されたＭＤＣＴ係数を受け取るように結合されると共に、各ＭＤＣＴ係数を整数値へと丸めるように動作できる丸めユニットと、
    前記回復されたＭＤＣＴストリームを受け取るように結合され、それに応答して、前記知覚コード化信号の再構成されたコピーを発生するための逆ＭＤＣＴ装置と、
    を備える請求項２に記載の装置。
  4. 前記変換ユニットは、更に、
    前記ロスなしのビットストリームを受け取るように結合されると共に、それに応答して、残余ＩｎｔＭＤＣＴ係数を発生するように動作できるエントロピーデコーダと、
    前記残余ＩｎｔＭＤＣＴ係数を整数値のＭＤＣＴ係数に加算して、ＩｎｔＭＤＣＴ係数を発生するための手段と、
    前記整数値のＭＤＣＴ係数及び前記ＩｎｔＭＤＣＴ係数の加算を受け取るように結合されて、ロスなしにコード化されたオーディオ信号の再構成されたコピーを発生するための逆ＩｎｔＭＤＣＴ装置と、
    を備える請求項３に記載の装置。
  5. 各リフトステージはリフトマトリクスに対応し、
    該リフトマトリクスは、前記変換の一ステップを実行するために使用され、４つのサブマトリクスと、これらサブマトリクスの２つとして２つの可逆整数マトリクスを１つの対角方向に伴うようなブロック三角マトリクスであり、
    前記４つのサブマトリクスのうち少なくとも一つは前記補助的な変換マトリクスである、請求項１に記載の装置。
  6. 各リフトマトリクスにおける前記可逆整数マトリクスは、１又は−１のいずれかである成分を伴う対角マトリクスである、請求項５に記載の装置。
  7. 前記変換エレメントは、３つのリフトステージを備える、請求項１及び５〜６のいずれか一項に記載の装置。
  8. 前記変換エレメントは、５つのリフトステージを備える、請求項１及び５〜６のいずれか一項に記載の装置。
  9. 前記変換エレメントは、８つのリフトステージを備える、請求項１及び５〜６のいずれか一項に記載の装置。
  10. オーディオ信号又はビデオ信号が前記デジタル信号として使用される、請求項１及び５〜９のいずれか一項に記載の装置。
  11. プログラムが記録されたコンピュータ読み取り可能な媒体であって、前記プログラムは、時間ドメインから周波数ドメインへ及びそれとは逆にデジタル信号のドメイン変換を実行する方法をコンピュータのプロセッサで実行させるように適応されたものであり、
    前記プログラムは、前記プロセッサに、前記デジタル信号に対する変換エレメントにより複数のリフトステージに従って前記ドメイン変換を実行させるコードを備え、前記デジタル信号のデータ記号は前記変換エレメントにデータベクトルとして提供され、各リフトステージは前記データベクトルの少なくとも一部に対する前記ドメイン変換の一ステップを実行し、前記ドメイン変換は離散コサイン変換（ＤＣＴ）−ＩＶの変換マトリクスに対応し、該変換マトリクスに従って前記データベクトルが変換され、前記複数のリフトステージのうちの少なくとも１つのリフトステージは、少なくとも１つの補助的な変換マトリクスと、丸めとを含み、前記補助的な変換マトリクスは、前記変換マトリクスそれ自体で構成されるか、又は低い次元の対応変換マトリクスで構成され、各リフトステージにおいて前記データベクトルの少なくとも一部が前記少なくとも１つの補助的な変換マトリクスによって変換され、前記データベクトルの前記少なくとも一部はリフトステージにおいて前記少なくとも１つの補助的な変換マトリクスとの乗算に基づいて処理され、
    前記プログラムは、前記プロセッサに、前記補助的な変換マトリクスによる前記ドメイン変換の後に前記丸めにより前記データベクトル又は前記データベクトルの少なくとも一部の丸め演算を実行させるコードを備え、
    前記プロセッサに前記ドメイン変換を実行させる前記プログラムの前記コードは、
    前記プロセッサに、複数のデータブロックを受け取らせ、各データブロックをＭＤＣＴ係数へとドメイン変換させるコードと、
    前記プロセッサに、前記ＭＤＣＴ係数の各々を受け取らせ、それに応答して、量子化されたＭＤＣＴ係数を発生させるコードと、
    前記プロセッサに、前記量子化されたＭＤＣＴ係数を受け取らせ、それに応答して、知覚的にコード化されたビットストリームを発生させるコードと、
    前記プロセッサに、前記量子化されたＭＤＣＴ係数を受け取らせ、前記ＭＤＣＴ係数を非量子化状態に回復させるコードと、
    を備え、
    前記プロセッサに前記丸め演算を実行させる前記プログラムの前記コードは、前記プロセッサに、前記回復されたＭＤＣＴ係数を受け取らせ、整数値のＭＤＣＴ係数を発生させるコードを備える、
    コンピュータ読み取り可能な媒体。
JP2006532242A 2003-09-29 2004-05-06 時間ドメインから周波数ドメインへ及びそれとは逆にデジタル信号のドメイン変換を実行する装置及び媒体 Expired - Fee Related JP4429316B2 (ja)

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
US50744003P 2003-09-29 2003-09-29
US50721003P 2003-09-29 2003-09-29
PCT/SG2004/000120 WO2005031595A1 (en) 2003-09-29 2004-05-06 Method for performing a domain transformation of a digital signal from the time domain into the frequency domain and vice versa

Publications (3)

Publication Number Publication Date
JP2007507789A JP2007507789A (ja) 2007-03-29
JP2007507789A5 JP2007507789A5 (ja) 2007-06-28
JP4429316B2 true JP4429316B2 (ja) 2010-03-10

Family

ID=34396342

Family Applications (3)

Application Number Title Priority Date Filing Date
JP2006532243A Pending JP2007507790A (ja) 2003-09-29 2004-05-06 時間ドメインから周波数ドメインへ及びそれとは逆にデジタル信号を変換する方法
JP2006532242A Expired - Fee Related JP4429316B2 (ja) 2003-09-29 2004-05-06 時間ドメインから周波数ドメインへ及びそれとは逆にデジタル信号のドメイン変換を実行する装置及び媒体
JP2009156963A Expired - Fee Related JP4942793B2 (ja) 2003-09-29 2009-07-01 時間ドメインから周波数ドメインへ及びそれとは逆にデジタル信号を変換する方法

Family Applications Before (1)

Application Number Title Priority Date Filing Date
JP2006532243A Pending JP2007507790A (ja) 2003-09-29 2004-05-06 時間ドメインから周波数ドメインへ及びそれとは逆にデジタル信号を変換する方法

Family Applications After (1)

Application Number Title Priority Date Filing Date
JP2009156963A Expired - Fee Related JP4942793B2 (ja) 2003-09-29 2009-07-01 時間ドメインから周波数ドメインへ及びそれとは逆にデジタル信号を変換する方法

Country Status (8)

Country Link
US (3) US20070276894A1 (ja)
EP (3) EP1668533A4 (ja)
JP (3) JP2007507790A (ja)
KR (2) KR100885437B1 (ja)
MX (2) MXPA06003509A (ja)
MY (2) MY139788A (ja)
TW (3) TW200527224A (ja)
WO (3) WO2005031596A1 (ja)

Families Citing this family (19)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
DE10331803A1 (de) * 2003-07-14 2005-02-17 Fraunhofer-Gesellschaft zur Förderung der angewandten Forschung e.V. Vorrichtung und Verfahren zum Umsetzen in eine transformierte Darstellung oder zum inversen Umsetzen der transformierten Darstellung
EP1668533A4 (en) * 2003-09-29 2013-08-21 Agency Science Tech & Res METHOD FOR CARRYING OUT A DOMAIN TRANSFORMATION OF A DIGITAL SIGNAL FROM THE TIME DOMAIN INTO THE FREQUENCY DOMAIN AND IN VICE
DE10345996A1 (de) * 2003-10-02 2005-04-28 Fraunhofer Ges Forschung Vorrichtung und Verfahren zum Verarbeiten von wenigstens zwei Eingangswerten
US7574064B2 (en) * 2005-04-12 2009-08-11 The Aerospace Corporation Fast lifting lossless wavelet transform
US7577307B2 (en) * 2005-04-12 2009-08-18 The Aerospace Corporation Fast adaptive lifting lossless wavelet transform
US8392176B2 (en) 2006-04-10 2013-03-05 Qualcomm Incorporated Processing of excitation in audio coding and decoding
US20100292986A1 (en) * 2007-03-16 2010-11-18 Nokia Corporation encoder
US20080288568A1 (en) * 2007-05-14 2008-11-20 Hou Hsieh S Low power Fast Hadamard transform
US8428957B2 (en) 2007-08-24 2013-04-23 Qualcomm Incorporated Spectral noise shaping in audio coding based on spectral dynamics in frequency sub-bands
US8548815B2 (en) 2007-09-19 2013-10-01 Qualcomm Incorporated Efficient design of MDCT / IMDCT filterbanks for speech and audio coding applications
EP2211335A1 (en) * 2009-01-21 2010-07-28 Fraunhofer-Gesellschaft zur Förderung der angewandten Forschung e.V. Apparatus, method and computer program for obtaining a parameter describing a variation of a signal characteristic of a signal
US8885701B2 (en) * 2010-09-08 2014-11-11 Samsung Electronics Co., Ltd. Low complexity transform coding using adaptive DCT/DST for intra-prediction
US8995532B2 (en) 2010-09-30 2015-03-31 Texas Instruments Incorporated Low complexity large transform
WO2015164575A1 (en) 2014-04-25 2015-10-29 Dolby Laboratories Licensing Corporation Matrix decomposition for rendering adaptive audio using high definition audio codecs
US10068577B2 (en) 2014-04-25 2018-09-04 Dolby Laboratories Licensing Corporation Audio segmentation based on spatial metadata
CN105895109B (zh) * 2016-05-10 2019-02-26 信阳师范学院 一种基于dwt和dct的数字语音取证和篡改恢复方法
KR20180089858A (ko) * 2017-02-01 2018-08-09 엘지전자 주식회사 레이어드 기븐스 변환을 이용하여 변환을 수행하는 방법 및 장치
EP3616196A4 (en) 2017-04-28 2021-01-20 DTS, Inc. Audio coder window and transform implementations
US12278630B2 (en) * 2021-09-24 2025-04-15 Altera Corporation Fast fourier transform (FFT) based digital signal processing (DSP) engine

Family Cites Families (21)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4031038A (en) * 1975-06-16 1977-06-21 The Dow Chemical Company Water insoluble chelate exchange resins having a crosslinked polymer matrix and pendant thereto a plurality of methyleneaminopyridine groups
US5339265A (en) * 1992-08-31 1994-08-16 University Of Maryland At College Park Optimal unified architectures for the real-time computation of time-recursive discrete sinusoidal transforms
US5523847A (en) * 1992-10-09 1996-06-04 International Business Machines Corporation Digital image processor for color image compression
JP3773260B2 (ja) * 1994-11-09 2006-05-10 コーニンクレッカ フィリップス エレクトロニクス エヌ ヴィ 画像又は音声信号を伝送するシステム及び方法
US5999656A (en) * 1997-01-17 1999-12-07 Ricoh Co., Ltd. Overlapped reversible transforms for unified lossless/lossy compression
US6058215A (en) * 1997-04-30 2000-05-02 Ricoh Company, Ltd. Reversible DCT for lossless-lossy compression
US6073153A (en) * 1998-06-03 2000-06-06 Microsoft Corporation Fast system and method for computing modulated lapped transforms
US6421464B1 (en) * 1998-12-16 2002-07-16 Fastvdo Llc Fast lapped image transforms using lifting steps
WO2000055757A1 (en) * 1999-03-17 2000-09-21 The Johns Hopkins University A fast multiplierless transform
US6681052B2 (en) * 2000-01-15 2004-01-20 Sony Corporation Methods and systems for performing inverse quantization and inverse weighting of DV video
WO2002045005A1 (en) * 2000-12-01 2002-06-06 Lizardtech, Inc. Method for lossless encoding of image data by approximating linear transforms and preserving selected properties
US6934676B2 (en) * 2001-05-11 2005-08-23 Nokia Mobile Phones Ltd. Method and system for inter-channel signal redundancy removal in perceptual audio coding
DE10129240A1 (de) * 2001-06-18 2003-01-02 Fraunhofer Ges Forschung Verfahren und Vorrichtung zum Verarbeiten von zeitdiskreten Audio-Abtastwerten
JP3971135B2 (ja) * 2001-07-11 2007-09-05 株式会社テクノマセマティカル Ｄｃｔ行列分解方法及びｄｃｔ装置
US7082450B2 (en) * 2001-08-30 2006-07-25 Nokia Corporation Implementation of a transform and of a subsequent quantization
JP3796432B2 (ja) 2001-10-31 2006-07-12 キヤノン株式会社 フィルタ処理装置およびフィルタ処理方法
DE10217297A1 (de) * 2002-04-18 2003-11-06 Fraunhofer Ges Forschung Vorrichtung und Verfahren zum Codieren eines zeitdiskreten Audiosignals und Vorrichtung und Verfahren zum Decodieren von codierten Audiodaten
JP3902990B2 (ja) * 2002-07-02 2007-04-11 キヤノン株式会社 アダマール変換処理方法及びその装置
DE10236694A1 (de) * 2002-08-09 2004-02-26 Fraunhofer-Gesellschaft zur Förderung der angewandten Forschung e.V. Vorrichtung und Verfahren zum skalierbaren Codieren und Vorrichtung und Verfahren zum skalierbaren Decodieren
DE10331803A1 (de) * 2003-07-14 2005-02-17 Fraunhofer-Gesellschaft zur Förderung der angewandten Forschung e.V. Vorrichtung und Verfahren zum Umsetzen in eine transformierte Darstellung oder zum inversen Umsetzen der transformierten Darstellung
EP1668533A4 (en) * 2003-09-29 2013-08-21 Agency Science Tech & Res METHOD FOR CARRYING OUT A DOMAIN TRANSFORMATION OF A DIGITAL SIGNAL FROM THE TIME DOMAIN INTO THE FREQUENCY DOMAIN AND IN VICE

Also Published As

Publication number Publication date
EP1668534A1 (en) 2006-06-14
JP2007507790A (ja) 2007-03-29
JP2009266250A (ja) 2009-11-12
EP1668534A4 (en) 2013-08-21
EP1668533A4 (en) 2013-08-21
EP1668533A1 (en) 2006-06-14
MXPA06003509A (es) 2007-01-25
MY137920A (en) 2009-03-31
KR20070026325A (ko) 2007-03-08
KR100885437B1 (ko) 2009-02-24
MXPA06003508A (es) 2007-01-25
WO2005031596A1 (en) 2005-04-07
WO2005031595A1 (en) 2005-04-07
TW200524290A (en) 2005-07-16
KR20070015363A (ko) 2007-02-02
US8126950B2 (en) 2012-02-28
US20080030385A1 (en) 2008-02-07
JP4942793B2 (ja) 2012-05-30
TWI343534B (en) 2011-06-11
MY139788A (en) 2009-10-30
TW200527224A (en) 2005-08-16
US20070276894A1 (en) 2007-11-29
JP2007507789A (ja) 2007-03-29
TW200529040A (en) 2005-09-01
US8126951B2 (en) 2012-02-28
TWI348282B (en) 2011-09-01
EP1687736A1 (en) 2006-08-09
WO2005031597A1 (en) 2005-04-07
US20070276893A1 (en) 2007-11-29
KR100885438B1 (ko) 2009-02-24

Similar Documents

Publication Publication Date Title
JP4942793B2 (ja) 時間ドメインから周波数ドメインへ及びそれとは逆にデジタル信号を変換する方法
CN1860527B (zh) 一种处理具有离散值序列的信号的设备和方法
US8195730B2 (en) Apparatus and method for conversion into a transformed representation or for inverse conversion of the transformed representation
CN100435137C (zh) 用于处理至少两个输入值的装置和方法
KR20060069254A (ko) 손실 및 무손실 2차원 데이터 압축을 위한 가역 변환
US20040220805A1 (en) Method and device for processing time-discrete audio sampled values
CN100520765C (zh) 确定给定变换函数的变换元素的过程和设备,数字信号变换方法和设备及计算机可读介质
Oraintara et al. The integer MDCT and its application in the MPEG layer III audio

Legal Events

Date Code Title Description
A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20070507

A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20070507

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20080828

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20080909

A601 Written request for extension of time

Free format text: JAPANESE INTERMEDIATE CODE: A601

Effective date: 20081209

A602 Written permission of extension of time

Free format text: JAPANESE INTERMEDIATE CODE: A602

Effective date: 20081216

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20090107

RD03 Notification of appointment of power of attorney

Free format text: JAPANESE INTERMEDIATE CODE: A7423

Effective date: 20090107

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20090303

A601 Written request for extension of time

Free format text: JAPANESE INTERMEDIATE CODE: A601

Effective date: 20090527

A602 Written permission of extension of time

Free format text: JAPANESE INTERMEDIATE CODE: A602

Effective date: 20090603

A601 Written request for extension of time

Free format text: JAPANESE INTERMEDIATE CODE: A601

Effective date: 20090703

A602 Written permission of extension of time

Free format text: JAPANESE INTERMEDIATE CODE: A602

Effective date: 20090710

A601 Written request for extension of time

Free format text: JAPANESE INTERMEDIATE CODE: A601

Effective date: 20090803

A602 Written permission of extension of time

Free format text: JAPANESE INTERMEDIATE CODE: A602

Effective date: 20090810

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20090902

TRDD Decision of grant or rejection written
A01 Written decision to grant a patent or to grant a registration (utility model)

Free format text: JAPANESE INTERMEDIATE CODE: A01

Effective date: 20091124

A01 Written decision to grant a patent or to grant a registration (utility model)

Free format text: JAPANESE INTERMEDIATE CODE: A01

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20091215

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20121225

Year of fee payment: 3

R150 Certificate of patent or registration of utility model

Ref document number: 4429316

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150

Free format text: JAPANESE INTERMEDIATE CODE: R150

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20131225

Year of fee payment: 4

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

LAPS Cancellation because of no payment of annual fees