JP2009266250A - 時間ドメインから周波数ドメインへ及びそれとは逆にデジタル信号を変換する方法 - Google Patents

時間ドメインから周波数ドメインへ及びそれとは逆にデジタル信号を変換する方法 Download PDF

Info

Publication number
JP2009266250A
JP2009266250A JP2009156963A JP2009156963A JP2009266250A JP 2009266250 A JP2009266250 A JP 2009266250A JP 2009156963 A JP2009156963 A JP 2009156963A JP 2009156963 A JP2009156963 A JP 2009156963A JP 2009266250 A JP2009266250 A JP 2009266250A
Authority
JP
Japan
Prior art keywords
matrix
conversion
receive
dct
blocks
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.)
Granted
Application number
JP2009156963A
Other languages
English (en)
Other versions
JP4942793B2 (ja
Inventor
Haibin Huang
ハイビン ファン，
Xiao Lin
シャオ リン，
Susanto Rahardja
スザント ラハルジャ，
Rongshan Yu
ロンシャン ユー，
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 JP2009266250A publication Critical patent/JP2009266250A/ja
Application granted granted Critical
Publication of JP4942793B2 publication Critical patent/JP4942793B2/ja
Anticipated expiration legal-status Critical
Expired - Fee Related legal-status Critical Current

Links

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR 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; CALCULATING OR 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; CALCULATING OR 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; CALCULATING OR 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)
  • Computational Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Data Mining & Analysis (AREA)
  • Pure & Applied Mathematics (AREA)
  • Mathematical Optimization (AREA)
  • Mathematical Analysis (AREA)
  • Databases & Information Systems (AREA)
  • General Engineering & Computer Science (AREA)
  • Software Systems (AREA)
  • Algebra (AREA)
  • Signal Processing (AREA)
  • Multimedia (AREA)
  • Discrete Mathematics (AREA)
  • Health & Medical Sciences (AREA)
  • Computational Linguistics (AREA)
  • Audiology, Speech & Language Pathology (AREA)
  • Human Computer Interaction (AREA)
  • Spectroscopy & Molecular Physics (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)

Abstract

【課題】整数型離散コサイン変換の丸め回数を少なくし、近似エラーを減少する。
【解決手段】順方向変換及び逆方向変換において、デジタル信号は、複数のブロックへとグループ分けされるデータ記号で構成され、各ブロックが所定数のデータ記号を含み、デジタル信号の２つのブロックを１つの変換エレメントにより変換するプロセスを備え、変換エレメントは、２つのサブマトリクスを含むブロック対角マトリクスに対応し、各サブマトリクスは、変換マトリクスを備え、且つ変換エレメントは、複数のリフトステージを備え、更に、各リフトステージは、補助的な変換及び丸めユニットによるデジタル信号のブロックの処理を含む。
【選択図】図３

Description

関連出願へのクロスレファレンス
本出願は、２００３年９月２９日に出願された米国プロビジョナル特許出願第６０／５０７，２１０号、及び２００３年９月２９に出願された米国プロビジョナル特許出願第６０／５０７，４４０号の優先権の利益を請求するもので、これら各々の内容は、その全体を全ての目的で参考としてここに援用する。
更に、以下の共通に所有する特許出願は、本出願と同日に出願されたもので、その全体をここに援用する。
「Methodfor Performing a Domain Transformation of a Digital Signal from the Time Domaininto the Frequency Domain 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」、ＡＥＳ第１１１回コンベンション、米国、ニューヨーク、２００１年９月、を参照されたい。
それ故、ＩｎｔＭＤＣＴの実現は、整数ＤＣＴ−ＩＶの効率的な実施に依存する。整数変換は、各ギブンズ回転を３つのリフトステップに置き換えることによりそれらのプロートタイプから直接コンバートすることができる。各リフトステップには１つの丸め演算があるので、整数変換の全丸め回数は、プロートタイプ変換のギブンズ回転回数の３倍である。離散的三角法変換（例えば、離散的フーリエ変換（ＤＦＴ）又は離散的コサイン変換（ＤＣＴ））の場合に、それに伴うギブンズ回転の回数は、通常、Ｎｌｏｇ２Ｎレベルにあり、ここで、Ｎは、ブロックのサイズであり、即ちデジタル信号が分割される各ブロックに含まれるデータ記号の量である。従って、全丸め回数も、直接コンバートされる整数変換のファミリーに対して、Ｎｌｏｇ２Ｎレベルにある。丸めのために、整数変換は、そのフローティングポイントプロートタイプしか近似しない。丸めの回数と共に近似エラーが増加する。
それ故、デジタル信号をより効率的な仕方でドメイン変換するためのシステム及び方法が要望される。
本発明は、デジタル信号をドメイン変換するためのシステム及び方法であって、入力データの２つのブロックが同じオペレーションにおいて同時にドメイン変換されるようなシステム及び方法を提供する。この構成は、有効な丸め演算の回数、ひいては、近似エラーを減少させる。
本発明の一実施形態において、変換関数を使用してデジタル信号を時間ドメインから周波数ドメインへ及びそれとは逆に変換するための方法が提供される。変換関数は、変換マトリクスを含み、デジタル信号は、複数のブロックへとグループ分けされたデータ記号を含み、各ブロックは所定数のデータ記号を含む。この方法は、デジタル信号の２つのブロックを１つの変換エレメントで変換することを備え、この変換エレメントは、２つのサブマトリクスを含むブロック対角マトリクスに対応し、各サブマトリクスは、変換マトリクスを含むと共に、変換エレメントは、複数のリフトステージを含み、更に、各リフトステージは、補助的な変換及び丸めユニットによるデジタル信号のブロックの処理を含む。
本発明のこれら及び他の特徴は、添付図面及び特定実施形態の詳細な説明から良好に理解されよう。
本発明の実施形態によるオーディオエンコーダの構成を示す図である。 図１に示すオーディオコーダに対応する、本発明の実施形態によるオーディオデコーダの構成を示す図である。 本発明による方法の実施形態を示すフローチャートである。 ＤＣＴ−ＩＶを変換関数として使用する本発明による方法の実施形態を示す図である。 図４に示す本発明の方法の実施形態に基づく逆変換のアルゴリズムを示す図である。 本発明の実施形態によるイメージアーカイブシステムの構成を示す図である。 ここに提案するシステム及び方法の性能を評価するのに使用される順方向及び逆方向変換コーダを示す図である。
図１は、本発明の実施形態によるオーディオエンコーダ１００の構成を示す図である。このオーディオエンコーダ１００は、変更型離散的コサイン変換（ＭＤＣＴ）に基づく従来の知覚ベースレイヤーコーダと、整数変更型離散的コサイン変換（ＩｎｔＭＤＣＴ）に基づくロスなし改善コーダとを備えている。
例えば、マイクロホン１１０により発生されてアナログ／デジタルコンバータ１１１によりデジタル化されたオーディオ信号１０９が、オーディオエンコーダ１００へ供給される。このオーディオ信号１０９は、複数のデータ記号を含む。このオーディオ信号１０９は、複数のブロックに分割され、各ブロックは、デジタル信号の複数のデータ記号を含み、又、各ブロックは、変更型離散的コサイン変換（ＭＤＣＴ）装置１０１により変換される。ＭＤＣＴ係数は、知覚モジュール１０２の助けで量子化装置１０３により量子化される。知覚モデルは、量子化エラーにより生じる可聴歪が低くなるように量子化装置１０３を制御する。量子化されたＭＤＣＴ係数は、その後、ビットストリームエンコーダ１０４によりエンコードされ、このエンコーダは、ロスのある知覚コード化出力ビットストリーム１１２を発生する。
ビットストリームエンコーダ１０４は、その入力を、ハフマンコード化又はランレングスコード化のような標準的な方法により、ロスのないように圧縮して、その入力より平均ビットレートが低い出力を発生する。又、入力オーディオ信号１０９は、ＩｎｔＭＤＣＴ装置１０５にも供給され、この装置は、ＩｎｔＭＤＣＴ係数を発生する。量子化装置１０３の出力である量子化されたＭＤＣＴ係数は、ＩｎｔＭＤＣＴ係数を予想するのに使用される。量子化されたＭＤＣＴ係数は、逆量子化装置１０６へ供給され、その出力（回復され又は非量子化のＭＤＣＴ係数）は、丸めユニット１０７へ供給される。
丸めユニットは、供給されたＭＤＣＴ係数を整数値へ丸め、整数値ＭＤＣＴとＩｎｔＭＤＣＴ係数との間の差である残余ＩｎｔＭＤＣＴ係数は、エントロピーコーダ１０８によりエントロピーコード化される。このエントロピーエンコーダは、ビットストリームエンコーダ１０４と同様のもので、その入力の平均ビットレートをロスなしに減少して、ロスなしの改善ビットストリーム１１３を発生する。このロスなしの改善ビットストリーム１１３は、知覚コード化ビットストリーム１１２と一緒に、入力オーディオ信号１０９を再構成するに必要な情報を最小限のエラーで搬送する。
図２は、図１に示すオーディオコーダ１００に対応する本発明の実施形態によるオーディオデコーダ２００の構成を示す。知覚コード化ビットストリーム２０７は、ビットストリームデコーダ２０１へ供給され、このデコーダは、図１のビットストリームエンコーダ１０４の動作に対して逆の動作を実行して、デコードされたビットストリームを発生する。デコードされたビットストリームは、逆量子化装置２０２へ供給され、その出力（回復されたＭＤＣＴ係数）は、逆ＭＤＣＴ装置２０３へ供給される。従って、再構成された知覚コード化オーディオ信号２０９が得られる。
ロスなしの改善ビットストリーム２０８は、エントロピーデコーダ２０４へ供給され、このデコーダは、図１のエントロピーエンコーダ１０８の動作とは逆の動作を実行し、対応する残余ＩｎｔＭＤＣＴ係数を発生する。逆量子化装置２０２の出力は、丸め装置２０５によって丸められて、整数値ＭＤＣＴ係数を発生する。この整数値ＭＤＣＴ係数は、残余ＩｎｔＭＤＣＴ係数に加算されて、ＩｎｔＭＤＣＴ係数を発生する。最終的に、逆ＩｎｔＭＤＣＴは、逆ＩｎｔＭＤＣＴ装置２０６によりＩｎｔＭＤＣＴ係数に適用されて、再構成されたロスなしのコード化オーディオ信号２１０を発生する。
図３は、ＤＣＴ−ＩＶを変換として使用すると共に、３つのリフトステージ、即ち第１リフトステージ３０１、第２リフトステージ３０２及び第３リフトステージ３０３を使用する本発明による方法の実施形態を示すフローチャート３００である。この方法は、好ましくは、図１のＩｎｔＭＤＣＴ装置１０５及び図２の逆ＩｎｔＭＤＣＴ装置２０６に使用されて、各々、ＩｎｔＭＤＣＴ及び逆ＩｎｔＭＤＣＴを実施する。図３において、ｘ １ 及びｘ ２ は、各々、デジタル信号の第１及び第２ブロックである。ｚは、中間信号であり、そしてｙ １ 及びｙ ２ は、各々、デジタル信号の第１及び第２ブロックに対応する出力信号である。
上述したように、ＤＣＴ−ＩＶアルゴリズムは、ロスなしのオーディオコード化において重要な役割を演じる。
従って、この実施形態では、本発明による変換エレメントに対応する変換マトリクスは、次の通りである。
リフトマトリクスの数、ひいては、変換エレメントにおけるリフトステージの数は、ＤＣＴ−ＩＶが変換関数である本発明のこの実施形態では、３である。
Ｎポイント実数入力シーケンスｘ（ｎ）のＤＣＴ−ＩＶは、次のように定義される。
ここで、ｘ １ 、ｘ ２ を２つの整数Ｎｘ１列ベクトルとする。これら列ベクトルｘ １ 、ｘ ２ は、本発明では１つの変換エレメントにより変換されるデジタル信号の２つのブロックに対応する。ｘ １ 、ｘ ２ のＤＣＴ−ＩＶ変換は、各々、ｙ １ 、ｙ ２ である。
（５）及び（６）を合成すると、次のようになる。
前記対角マトリクスは、変換エレメントが、本発明により、対応しているブロック対角マトリクスである。
これは、前記式が次のような簡単な代数変更により変化された場合にも、本発明の範囲内である。
式（１０）は、（３）におけるＤＣＴ−ＩＶプロパティを使用して容易に照合することができる。式（１０）を用いると、式（８）は、次のように表わすことができる。
式（１１）における３つのリフトマトリクスは、図３に示す３つのリフトステージに対応する。
式（１１）から、１つの変換エレメントで２つの整数ＤＣＴ−ＩＶを計算する次の整数ＤＣＴ−ＩＶアルゴリズムが導出される。
図４は、ＤＣＴ−ＩＶを変換関数として使用する本発明による方法の実施形態を示す。この実施形態は、ＩｎｔＭＤＣＴを実施するために図１に示すオーディオコーダ１００に使用される。図３と同様に、ｘ １ 及びｘ ２ は、入力デジタル信号の２つのブロックであり、ｚは、中間信号であり、そしてｙ １ 及びｙ ２ は、出力信号の対応ブロックである。
図４に示す３つのリフトステージは、式（１１）における３つのリフトマトリクスに対応する。
図４に示すように、時間対周波数ドメインの整数変換は、次のように決定される。
第１ステージ４０１において、ｘ ２ がＤＣＴ−ＩＶ変換４０２により変換され、次いで、ＤＣＴ−ＩＶ係数が丸められる（４０３）。丸められたＤＣＴ−ＩＶ係数は、次いで、ｘ １ ４０４に加算される。従って、中間信号ｚが発生される。従って、中間信号ｚは、次の式を満足する。
第２ステージ４０５において、ｚがＤＣＴ−ＩＶ変換４０６により変換され、ＤＣＴ−ＩＶ係数が丸められる（４０７）。次いで、丸められたＤＣＴ−ＩＶ係数から、ｘ １ が減算される。従って、出力信号ｙ １ が発生される。従って、出力信号ｙ １ は、次の式を満足する。
第３ステージ４０９において、ｙ １ がＤＣＴ−ＩＶ変換４１０により変換され、次いで、ＤＣＴ−ＩＶ係数が丸められる（４１１）。丸められたＤＣＴ−ＩＶ係数が、次いで、ｚから減算される。従って、出力信号ｙ ２ が発生される。従って、出力信号ｙ ２ は、次の式を満足する。
図５は、ＤＣＴ−ＩＶを変換関数として使用する本発明の方法の実施形態に基づいて逆変換を行うためのアルゴリズムを示す。この実施形態は、逆ＩｎｔＭＤＣＴを実施するために図２に示すオーディオデコーダ２００に使用される。図５に示すアルゴリズムは、図４に示すアルゴリズムの逆である。異なる信号ｙ １ 、ｙ ２ 、ｘ １ 、ｘ ２ 、ｚに対する表示は、図４の表示に対応して選択されたものである。
図５に示すように、周波数対時間ドメインの整数変換は、次のように決定される。
第１ステージ５０１において、ｙ １ がＤＣＴ−ＩＶ変換５０２により変換され、ＤＣＴ−ＩＶ係数が丸められる（５０３）。丸められたＤＣＴ−ＩＶ係数は、次いで、ｙ ２ に加算される（５０４）。従って、中間信号ｚが発生される。従って、中間信号ｚは、次の式を満足する。
第２ステージ５０５において、ｚがＤＣＴ−ＩＶ変換５０６により変換され、ＤＣＴ−ＩＶ係数が丸められる。次いで、その丸められたＤＣＴ−ＩＶ係数から、ｙ １ が減算される。従って、信号ｘ ２ が発生される。従って、信号ｘ ２ は、次の式を満足する。
第３ステージ５０９において、ｘ ２ がＤＣＴ−ＩＶ変換５１０により変換され、ＤＣＴ−ＩＶ係数が丸められる（５１１）。この丸められたＤＣＴ−ＩＶ係数が、次いで、ｚから減算される。従って、信号ｘ １ が発生される。従って、ｘ １ は、次の式を満足する。
式（１３ａ）から（１３ｃ）に基づくアルゴリズムは、式（１２ａ）から（１２ｃ）に基づくアルゴリズムに対して逆であることが明らかである。従って、図１及び２に示すエンコーダ及びデコーダに使用した場合には、これらアルゴリズムは、ロスのないオーディオコード化のための方法及び装置をもたらす。
以下に述べる本発明の実施形態では、上述した方法がイメージアーカイブシステムに使用される。
式（１２ａ）から（１２ｃ）及び（１３ａ）から（１３ｃ）は、更に、２つのＮｘＮ整数ＤＣＴ−ＩＶを計算するために、３つのＮｘＮＤＣＴ−ＩＶ、３つのＮｘ１丸め、及び３つのＮｘ１加算が必要であることを示している。それ故、１つのＮｘＮ整数ＤＣＴ−ＩＶについては、平均が次のようになる。

但し、ＲＣ（．）は、全丸め回数であり、ＡＣ（．）は、全演算回数である。直接コンバートされる整数ＤＣＴ−ＩＶアルゴリズムに比して、ここに提案する整数ＤＣＴ−ＩＶアルゴリズムは、ＲＣをレベルＮｌｏｇ２ＮからＮへ減少する。
式（１５）で示されたように、ここに提案する整数ＤＣＴ−ＩＶアルゴリズムの演算の複雑さは、ＤＣＴ−ＩＶアルゴリズムより約５０％も高い。しかしながら、ＲＣも考慮した場合には、ここに提案するアルゴリズムの合成複雑さ（ＡＣ＋ＲＣ）は、直接コンバート式の整数アルゴリズムを著しく越えることはない。アルゴリズムの複雑さの厳密な分析は、使用するＤＣＴ−ＩＶアルゴリズムに依存する。
図４及び５に示すように、ここに提案する整数ＤＣＴ−ＩＶアルゴリズムは、構造が簡単でモジュール状である。これは、何らかの既存のＤＣＴ−ＩＶアルゴリズムをそのＤＣＴ−ＩＶ計算ブロックに使用することができる。ここに提案するアルゴリズムは、ＩｎｔＭＤＣＴを要求するアプリケーション、例えば、ＭＰＥＧ−４オーディオ拡張３参照モデル０に適している。
図６は、本発明の実施形態によるイメージアーカイブシステムの構成を示す。
図６において、イメージソース６０１、例えば、カメラは、アナログイメージ信号を発生する。このイメージ信号は、アナログ／デジタルコンバータ６０２により処理されて、それに対応するデジタルイメージ信号を発生する。このデジタルイメージ信号は、時間ドメインから周波数ドメインへの変換を含むロスなしのイメージエンコーダ６０３によりロスなしにエンコードされる。この実施形態では、時間ドメインは、イメージの座標スペースに対応する。ロスなしのコード化イメージ信号が、記憶装置６０４、例えば、ハードディスク又はＤＶＤに記憶される。イメージが必要とされるときには、ロスなしのコード化イメージ信号が、記憶装置６０４からフェッチされて、ロスなしのイメージデコーダ６０５へ送られ、このデコーダは、ロスなしのイメージエンコーダ６０３に対応するもので、ロスなしのコード化イメージ信号をデコードして、データロスを伴わずに元のイメージ信号を再構成する。
イメージ信号のこのようなロスなしのアーカイブは、例えば、イメージが半導体ウェハのエラーマップであって、後で分析するためにこれを記憶しなければならない場合に、重要である。
本発明のこの実施形態において、図３から５に示す方法の実施形態は、ロスなしのイメージエンコーダ６０３及びロスなしのイメージデコーダ６０５に使用される。上述したように、図３から５に示す方法の実施形態は、可逆の変換を提供し、従って、特に、ロスなしのイメージコード化方法を提供する。
本発明による方法は、オーディオ及びイメージ信号に限定されない。他のデジタル信号、例えば、ビデオ信号も、本発明の方法により変換することができる。
以下、本発明により時間ドメインから周波数ドメインへ及びそれとは逆にデジタル信号を変換する方法の更に別の実施形態を説明する。
本発明のこの実施形態では、ドメイン変換がＤＣＴ変換であり、従って、ブロックサイズＮは、何らかの整数である。一実施形態では、Ｎが２の累乗である。
各リフトマトリクスＬ２Ｎに対して、リフトステージの可逆の整数対整数マッピングは、ここに援用する参照文献「Factoring Wavelet Transforms into Lifting Steps」、Tech. Report、I. Daubechies及びW. Sweldens著、Bell Laboratories, LucentTechnologies、１９９６年、に説明された２ｘ２リフトステップと同様に実現される。唯一の相違は、単一の変数ではなく、ベクトルに丸め演算が適用されることである。
他の実施形態の前記説明において、リフトマトリクスに対してリフトステージをいかに実現するかは既に述べた。従って、リフトマトリクスに対応するリフトステージの説明は、以下では省略する。
前記式の右側を構成しているマトリクスは、次のように説明される。
Ｐ１は、次式で与えられる第１の置換マトリクスである。
Ｐ２は、第２の置換マトリクスであり、その一例が、次のＭＡＴＬＡＢスクリプトにより発生される。
一例として、Ｎが４であるときには、Ｐ２が次のような８ｘ８マトリクスである。
Ｐ３は、第３の置換マトリクスであり、その一例が、次のＭＡＴＬＡＢスクリプトにより発生される。
一例として、Ｎが４であるときには、Ｐ３が次のような８ｘ８マトリクスである。
式（４２）から、次元Ｎｘ１の２つの入力信号に対して整数ＤＣＴを計算することができる。
式（４２）は、ＤＣＴ−ＩＶ変換ドメインを記述するリフトマトリクス因数分解を与えるので、そのリフトマトリクスを、ここに示すように使用して、与えられた入力信号のドメイン変換を計算することができる。
式（４２）は、次のようにして導出することができる。
次の分解は、Wang、Zhongde著の「On Computing the Discrete Fourier and Cosine Transforms」、IEEE Transactions onAcoustics, Speech and Signal Processing、第ＡＳＳＰ−３３巻、第４号、１９８５年１０月、からの開示を使用して導出することができる。
この実施形態では、ドメイン変換を計算するときに、４Ｎ回の丸め演算しか必要とされず、これについて、以下に述べる。
ここで、各々、α（＊）を実数加算の回数とし、μ（＊）を実数乗算の回数とし、γ（＊）を実数丸めの回数とする。ここに提案するＩｎｔＤＣＴアルゴリズムの場合、次のようになる。
α（ＩｎｔＤＣＴ）＝１１Ｎ＋３α（ＤＣＴ−ＩＶ）
μ（ＩｎｔＤＣＴ）＝９Ｎ＋３μ（ＤＣＴ−ＩＶ）
γ（ＩｎｔＤＣＴ）＝８Ｎ
前記結果は、データサンプルの２つのブロックに対するものである。というのは、ここに提案するＩｎｔＤＣＴアルゴリズムは、それらを一緒に処理するからである。従って、データサンプルの１つのブロックに対して、計算の回数が半減し、次のようになる。
α１（ＩｎｔＤＣＴ）＝５．５Ｎ＋１．５α（ＤＣＴ−ＩＶ）
μ１（ＩｎｔＤＣＴ）＝４．５Ｎ＋１．５μ（ＤＣＴ−ＩＶ）
γ１（ＩｎｔＤＣＴ）＝４Ｎ
但し、α１、μ１及びγ１は、各々、１つのサンプルブロックに対する実数加算の回数、実数乗算の回数、及び実数丸めの回数である。
ＤＣＴ−ＩＶ計算の場合には、ここに援用する参照文献「Signal Processing with lapped Transforms」、H. S. Malvar著、マサチューセッツ州ノルウッド、アーテック・ハウス、１９９２年、第１９９−２０１ページ、に掲載されたＦＦＴベースのアルゴリズムを使用することができ、これについては、
α（ＤＣＴ−ＩＶ）＝１．５Ｎｌｏｇ２Ｎ
μ（ＤＣＴ−ＩＶ）＝０．５Ｎｌｏｇ２Ｎ＋Ｎ
その結果、次のようになる。
α１（ＩｎｔＤＣＴ）＝２．２５Ｎｌｏｇ２Ｎ＋５．５Ｎ
μ１（ＩｎｔＤＣＴ）＝０．７５Ｎｌｏｇ２Ｎ＋６Ｎ
以下、時間ドメインから周波数ドメインへ及びそれとは逆にデジタル信号を変換する方法の更に別の実施形態を説明する。
この実施形態では、離散的高速フーリエ変換（ＦＦＴ）がドメイン変換として使用される。
Ｆを、正規化されたＦＦＴのＮｘＮ変換マトリクスとする。

但し、Ｎは、変換サイズであり、ｍ及びｎは、マトリクス指数である。
この実施形態のもとでは、次元ＮｘＮの置換マトリクスＰは、指数０又は１を含むマトリクスである。これにＮｘ１ベクトル（入力信号のマトリクス表示）を乗算した後、ベクトルにおけるエレメントの順序が変更される。
この実施形態では、リフトマトリクスは、次の式の２Ｎｘ２Ｎマトリクスとして定義される。

但し、Ｐ１及びＰ２は、２つの置換マトリクスであり、Ｏは、ＮｘＮのゼロマトリクスであり、Ａは、任意のＮｘＮマトリクスである。リフトマトリクスＬについては、Ｉ．ダウベッチーズ氏のここに援用する前記参照文献における２ｘ２リフトステップと同様に実現される。しかしながら、上述したように、単一の変数ではなく、丸めがベクトルに適用される。Ｌの転置、ＬＴも、リフトマトリクスであることが明らかであろう。
更に、Ｔを２Ｎｘ２Ｎ変換マトリクスとする。
従って、変更型変換マトリクスＴ（従って、ドメイン変換それ自体）は、リフトマトリクス因数分解として表現することができる。
式（５３）において、かぎカッコ内のブランクスペースは、全てゼロのマトリクスエレメントを表わす。
式（５１）から明らかなように、リフトマトリクスの因数分解は、ここに述べる方法を使用して２つのＮｘ１複素数ベクトルに対する整数ＦＦＴを計算するのに使用できる。
この実施形態のもとでは、ドメイン変換の計算は、３Ｎ回の丸め演算しか必要とせず、これについて、以下に述べる。
α（＊）を実数加算の回数とする。
μ（＊）を実数乗算の回数とする。
γ（＊）を実数丸め演算の回数とする。
ここに提案するＩｎｔＦＦＴアルゴリズムの場合に、次のものが得られる。
α（ＩｎｔＦＦＴ）＝６Ｎ＋３α（ＦＦＴ）
μ（ＩｎｔＦＦＴ）＝３μ（ＦＦＴ）
γ（ＩｎｔＦＦＴ）＝６Ｎ
前記結果は、データサンプルの２つのブロックに対するものである。というのは、ここに提案するＩｎｔＦＦＴアルゴリズムは、それらを一緒に処理するからである。従って、データサンプルの１つのブロックに対して、計算の回数が半減し、次のようになる。
α１（ＩｎｔＦＦＴ）＝３Ｎ＋１．５α（ＦＦＴ）
μ１（ＩｎｔＦＦＴ）＝１．５μ（ＦＦＴ）
γ１（ＩｎｔＦＦＴ）＝３Ｎ
但し、α１、μ１及びγ１は、各々、１つのサンプルブロックに対する実数加算の回数、実数乗算の回数、及び実数丸め演算の回数である。
ＦＦＴ計算に対して、スプリット−ラディックスＦＦＴ（ＳＲＦＦＴ）アルゴリズムを使用することができ、これについては、
α（ＳＲＦＦＴ）＝３Ｎｌｏｇ２Ｎ−３Ｎ＋４
μ（ＳＲＦＦＴ）＝Ｎｌｏｇ２Ｎ−３Ｎ＋４
その結果、次のようになる。
α１（ＩｎｔＦＦＴ）＝４．５Ｎｌｏｇ２Ｎ−１．５Ｎ＋６
μ１（ＩｎｔＦＦＴ）＝１．５Ｎｌｏｇ２Ｎ−４．５Ｎ＋６
図７は、上述したＤＣＴ変換技術及び上述したＦＦＴドメイン変換の変換精度を評価するのに使用される順方向及び逆方向変換コーダを示す。テストは、ここに援用する「Coding of Moving Pictures and Audio: Work plan for Evaluation ofInteger MDCT for FGS to Lossless Experimentation Framework」、ＩＳＯ／ＩＥＣ ＪＴＣ １／ＳＣ ２９／ＷＧ １１ Ｎ５５７８、Pattaya、Thailand、２００３年３月、に記載されたＭＰＥＧ−４ロスなしオーディオコードグループにより提案された評価規格に基づいて変換の平均２乗エラー（ＭＳＥ）を測定することを含む。
より詳細には、ＩｎｔＤＣＴ及び整数逆ＤＣＴ（ＩｎｔＩＤＣＴ）に対するＭＳＥは、次のように表わされる。

但し、エラー信号ｅは、図１の場合と同様に、ＩｎｔＤＣＴについてはｅｆ、及びＩｎｔＩＤＣＴについてはｅｔである。Ｋは、評価に使用されるサンプルブロックの合計数である。
ＩｎｔＦＦＴ及び整数逆ＦＦＴ（ＩｎｔＩＦＦＴ）に対するＭＳＥは、次のように表わされる。
両ドメイン変換の場合に、１５個の異なるタイプの音楽ファイルで合計４５０秒が４８ｋＨｚ／１６ビットテストセットに使用された。表１は、テスト結果を示す。
表１から明らかなように、本発明のシステム及び方法を使用して生じたＭＳＥは、非常に僅かであり、そして従来のシステムとは異なり、処理部とっくサイズとは実質的に独立したものである。ＤＣＴ−ＩＶドメイン変換を参照すれば、ＭＳＥは、４０９６ビットまでのブロックザイズＮの増加と共に僅かに増加するだけである。ＦＦＴのＭＳＥは、いっそう良好で、４０９６ビットまでのブロックサイズに対して０．４の一定ＭＳＥを示す。本発明の実証された性能を、より長いブロックサイズに対する現在の能力及び増加するニーズの観点から見ると、本発明の効果が明らかとなろう。
援用する参照文献
以下の文書を、参考としてここに援用する。
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. プログラムが記録されたコンピュータ読み取り可能な媒体であって、前記プログラムは、離散コサイン変換（ＤＣＴ）−ＩＶの変換マトリクスを含む変換関数を使用して時間ドメインから周波数ドメインへ及びそれとは逆にデジタル信号を変換する方法をコンピュータのプロセッサに実行させるように適応され、前記デジタル信号は、データ記号で構成されて、複数のブロックへと分割され、各ブロックが同数のデータ記号を含んでおり、
    前記プログラムが、前記プロセッサに、前記デジタル信号の２つのブロックの変換を１つの変換エレメントにより実行させるコードを備え、前記変換エレメントは、２つのサブマトリクスを含むブロック対角マトリクスに対応し、各サブマトリクスは、前記変換マトリクスに対応し、前記変換エレメントは、前記変換を複数のリフトステージで実行し、各リフトステージは、前記二つのブロックに対する前記変換の一ステップを実行し、前記二つのブロックは、該二つのブロックのうちの一つとＤＣＴ−ＩＶマトリクスとの乗算に基づいてリフトステージにおいて処理され、
    前記プロセッサに前記変換を実行させる前記プログラムの前記コードは、
    前記プロセッサに、前記複数のデータブロックを受け取らせ、各データブロックをＭＤＣＴ係数へとドメイン変換させるためのコードと、
    前記プロセッサに、前記ＭＤＣＴ係数の各々を受け取らせ、それに応答して、量子化されたＭＤＣＴ係数を発生させるためのコードと、
    前記プロセッサに、前記量子化されたＭＤＣＴ係数を受け取らせ、それに応答して、知覚的にコード化されたビットストリームを発生させるためのコードと、
    前記プロセッサに、前記量子化されたＭＤＣＴ係数を受け取らせ、前記ＭＤＣＴ係数を非量子化状態に回復させるためのコードと、
    前記プロセッサに、前記回復されたＭＤＣＴ係数を受け取らせ、整数値のＭＤＣＴ係数を発生させるためのコードと、
    を含む、コンピュータ読み取り可能な媒体。
JP2009156963A 2003-09-29 2009-07-01 時間ドメインから周波数ドメインへ及びそれとは逆にデジタル信号を変換する方法 Expired - Fee Related JP4942793B2 (ja)

Applications Claiming Priority (4)

Application Number Priority Date Filing Date Title
US50744003P 2003-09-29 2003-09-29
US50721003P 2003-09-29 2003-09-29
US60/507,210 2003-09-29
US60/507,440 2003-09-29

Related Parent Applications (1)

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

Publications (2)

Publication Number Publication Date
JP2009266250A true JP2009266250A (ja) 2009-11-12
JP4942793B2 JP4942793B2 (ja) 2012-05-30

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 (2)

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 時間ドメインから周波数ドメインへ及びそれとは逆にデジタル信号のドメイン変換を実行する装置及び媒体

Country Status (8)

Country Link
US (3) US8126951B2 (ja)
EP (3) EP1687736A1 (ja)
JP (3) JP2007507790A (ja)
KR (2) KR100885437B1 (ja)
MX (2) MXPA06003508A (ja)
MY (2) MY137920A (ja)
TW (3) TW200527224A (ja)
WO (3) WO2005031597A1 (ja)

Families Citing this family (18)

* 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
JP2007507790A (ja) * 2003-09-29 2007-03-29 エージェンシー フォー サイエンス，テクノロジー アンド リサーチ 時間ドメインから周波数ドメインへ及びそれとは逆にデジタル信号を変換する方法
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
US10068577B2 (en) 2014-04-25 2018-09-04 Dolby Laboratories Licensing Corporation Audio segmentation based on spatial metadata
EP3134897B1 (en) 2014-04-25 2020-05-20 Dolby Laboratories Licensing Corporation Matrix decomposition for rendering adaptive audio using high definition audio codecs
CN105895109B (zh) * 2016-05-10 2019-02-26 信阳师范学院 一种基于dwt和dct的数字语音取证和篡改恢复方法
KR20180089858A (ko) * 2017-02-01 2018-08-09 엘지전자 주식회사 레이어드 기븐스 변환을 이용하여 변환을 수행하는 방법 및 장치
WO2018201113A1 (en) 2017-04-28 2018-11-01 Dts, Inc. Audio coder window and transform implementations

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2002103684A1 (de) * 2001-06-18 2002-12-27 Fraunhofer-Gesellschaft zur Förderung der angewandten Forschung e. V. Verfahren und vorrichtung zum verarbeiten von zeitdiskreten audio-abtastwerten
WO2003088212A1 (de) * 2002-04-18 2003-10-23 Fraunhofer-Gesellschaft zur Förderung der angewandten Forschung e. V Vorrichtung und verfahren zum codieren eines zeitdiskreten audiosignals und vorrichtung und verfahren zum decodieren von codierten audiodaten

Family Cites Families (19)

* 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
EP0739569B1 (en) * 1994-11-09 1999-07-28 Koninklijke Philips Electronics N.V. System and method for transmitting picture or sound signals
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
EP1311975A4 (en) * 1999-03-17 2003-05-21 Univ Johns Hopkins Fast transformation without multiplier
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
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 キヤノン株式会社 フィルタ処理装置およびフィルタ処理方法
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
JP2007507790A (ja) * 2003-09-29 2007-03-29 エージェンシー フォー サイエンス，テクノロジー アンド リサーチ 時間ドメインから周波数ドメインへ及びそれとは逆にデジタル信号を変換する方法

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2002103684A1 (de) * 2001-06-18 2002-12-27 Fraunhofer-Gesellschaft zur Förderung der angewandten Forschung e. V. Verfahren und vorrichtung zum verarbeiten von zeitdiskreten audio-abtastwerten
WO2003088212A1 (de) * 2002-04-18 2003-10-23 Fraunhofer-Gesellschaft zur Förderung der angewandten Forschung e. V Vorrichtung und verfahren zum codieren eines zeitdiskreten audiosignals und vorrichtung und verfahren zum decodieren von codierten audiodaten

Also Published As

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

Similar Documents

Publication Publication Date Title
JP4942793B2 (ja) 時間ドメインから周波数ドメインへ及びそれとは逆にデジタル信号を変換する方法
US8195730B2 (en) Apparatus and method for conversion into a transformed representation or for inverse conversion of the transformed representation
Liang et al. Fast multiplierless approximations of the DCT with the lifting scheme
KR101226094B1 (ko) 5-포인트 ｄｃｔ-ⅱ, ｄｃｔ-ⅳ, 및 ｄｓｔ-ⅳ 및 아키텍쳐들의 계산을 위한 고속 알고리즘들
RU2325708C2 (ru) Устройство и способ обработки сигнала, имеющего последовательность дискретных значений
JP5086274B2 (ja) 共通因数を用いる変換
JP4328805B2 (ja) 少なくとも２つの入力値を処理するための装置および方法
JP3814611B2 (ja) 時間離散オーディオサンプル値を処理する方法と装置
CN100520765C (zh) 确定给定变换函数的变换元素的过程和设备,数字信号变换方法和设备及计算机可读介质
Liang et al. Fast multiplierless approximation of the DCT with the lifting scheme
Huang et al. Integer fast modified cosine transform

Legal Events

Date Code Title Description
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: 20120207

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: 20120228

R150 Certificate of patent or registration of utility model

Ref document number: 4942793

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: 20150309

Year of fee payment: 3

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