JP2008015416A - 音響信号に対する情報の埋め込み装置および音響信号からの情報の抽出装置 - Google Patents
音響信号に対する情報の埋め込み装置および音響信号からの情報の抽出装置 Download PDFInfo
- Publication number
- JP2008015416A JP2008015416A JP2006189196A JP2006189196A JP2008015416A JP 2008015416 A JP2008015416 A JP 2008015416A JP 2006189196 A JP2006189196 A JP 2006189196A JP 2006189196 A JP2006189196 A JP 2006189196A JP 2008015416 A JP2008015416 A JP 2008015416A
- Authority
- JP
- Japan
- Prior art keywords
- frame
- spectrum
- acoustic
- information
- component
- 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
Links
Images
Landscapes
- Signal Processing For Digital Recording And Reproducing (AREA)
Abstract
【課題】 従来と同等以上の抽出精度を維持しつつ、より多くの情報を埋め込むことが可能な音響信号に対する情報の埋め込み装置、音響信号からの情報の抽出装置を提供する。
【解決手段】 音響信号を所定区間に区分し、埋め込むべき２ビットおよびこの２ビットの付加情報の１ワード内の位置によりとり得る８値に応じて、前記所定区間の低周波成分の状態を変更する。具体的には、低周波成分をさらにより高い低周波成分とより低い低周波成分に分けるとともに、時間的に前と後に分け、得られる４つの領域を埋め込むべき８値に応じて変化させる。
【選択図】 図２
【解決手段】 音響信号を所定区間に区分し、埋め込むべき２ビットおよびこの２ビットの付加情報の１ワード内の位置によりとり得る８値に応じて、前記所定区間の低周波成分の状態を変更する。具体的には、低周波成分をさらにより高い低周波成分とより低い低周波成分に分けるとともに、時間的に前と後に分け、得られる４つの領域を埋め込むべき８値に応じて変化させる。
【選択図】 図２
Description
本発明は、ＣＤ・ＤＶＤ等を用いた民生・業務用途における鑑賞用のパッケージ音楽分野、放送事業者等が商業目的で配信する放送・ネットワーク音楽配信分野における音楽著作権の保護（不正コピーの監視）および音楽属性情報の提供（楽曲タイトル検索サービス）分野、ミュージアム、イベント会場における展示説明ナレーションに連動した文字情報の提供サービス分野、放送番組やＣＤ／ＤＶＤパッケージの音声信号からＵＲＬなどの情報を抽出し、携帯電話を用いて所定のコンテンツに関連するｗｅｂサイトにアクセスして詳細情報を抽出したり、アンケートに回答したりする非接触なインターネットのゲートウェイサービス分野に関する。
最近、流れている音楽のタイトル等を知ることができる楽曲属性情報の提供サービスとして、放送された音楽に対して日時と地域を放送局に照会したり、携帯電話で流れている音楽断片を録音してデータベースに登録されているメロディーと照合したりするサービスが実用化されている（例えば、特許文献１、２参照）。
特許文献１、２に記載の発明では、録音した音楽の断片と、データベースに登録されているメロディーと照合するため、データベースに登録される楽曲が増えると、処理負荷が増え、類似したメロディーを誤判定する可能性が増える。そこで、曲名やアーチスト情報などの楽曲属性情報を不可聴な電子透かしとして音響信号に埋め込む手法も提案されている（例えば、特許文献３〜６参照）。
特許文献１〜６に記載の手法では、埋め込み可能な情報量が少なく、音質が少なからず劣化し、各種信号処理により透かし情報が消失し、またアナログコピーに対しては、透かし検出が困難であるという問題がある。そこで、本出願人は、複数チャンネル有する音響信号の低周波数成分の比率を属性情報のビット値に応じて変更することにより、属性情報（付加情報）を埋め込む手法を提案した（特許文献７参照）。
特許文献７に記載の発明は、ステレオ音響信号にのみ対応したものであったが、本出願人は、モノラル音響信号にも対応可能なように改良を加え、低周波成分のうち、より高い低周波成分とより低い低周波成分を利用して付加情報を埋め込む手法を提案した（特許文献８参照）。
特許文献８に記載の発明では、低い方の低周波数帯が、一般的なマイクロフォンの感度範囲に含まれないことになる場合が多く、埋め込まれた情報を正しく識別することができない場合が生じるため、本出願人は、音響信号の所定区間における低周波成分の先頭部分と後部部分のスペクトル強度の割合を変更することにより、付加情報を埋め込む手法を提案した（特許文献９参照）。
上記特許文献７〜９に記載の発明では、いずれも周波数４４．１ｋＨｚでサンプリングした音響信号について、４０９６サンプルを１音響フレームとし、この１音響フレームについて１ビットの情報を埋め込む処理を行っている。実際には、区切りを示す情報や、埋め込みエラーが生じた場合に対処する音響フレームも必要となり、全ての音響フレームに埋め込み対象とする有効な情報を埋め込むことができるわけではないため、現実には、毎秒１０ビット程度の情報量が埋め込み可能となっている。
特開２００２−２５９４２１号公報
特開２００３−１５７０８７号公報
特開平１１−１４５８４０号公報
特開平１１−２１９１７２号公報
特許第３３２１７６７号公報
特開２００３−９９０７７号公報
特願２００５−５１５７号
特願２００５−５８８２４号
特願２００５−１４７７４３号
しかしながら、付加情報として例えば“http:www.○○○.co.jp”といったＵＲＬを埋め込もうとすると、最低でも１８秒程度を要することになり、抽出時にも同様の時間を要するため、実用的に問題がある。埋め込み可能な情報量を増やすためには、単純に上記１音響フレームを構成するサンプル数を１／２にするという手法がある。
しかし、特許文献７〜９では、音響フレームに周波数変換を行うことにより処理するものであるため、１音響フレームのサンプル数を１／２にすると、周波数変換により得られる周波数成分の要素数も１／２となり、１音響フレームにおける解析のための情報量は、現状に比べ１／４となってしまう。このため、特に抽出時において、正確な情報を抽出することが困難となる。また、通常、付加情報は所定ビット数を１ワードとした場合に複数ワードで構成されるが、処理単位である１ワードの区切りを示す情報を埋め込むために１音響フレームを費やさなければならず、埋め込み効率を悪化させるという問題もある。
そこで、本発明は、従来と同等以上の抽出精度を維持しつつ、より多くの情報を埋め込むことが可能な音響信号に対する情報の埋め込み装置、音響信号からの情報の抽出装置を提供することを課題とする。
上記課題を解決するため、本発明では、時系列のサンプル列で構成される音響信号に対して、付加情報を聴取不能な状態で埋め込む装置であって、前記音響信号より、所定数のサンプルを音響フレームとして読み込む音響フレーム読込手段と、前記音響フレームに対して第１窓関数、第２窓関数、第３窓関数を用いてそれぞれ周波数変換を行い、前記第１窓関数に対応するスペクトルである第１窓スペクトル、前記第２窓関数に対応するスペクトルである第２窓スペクトル、前記第３窓関数に対応するスペクトルである第３窓スペクトルを生成する周波数変換手段と、前記生成された第１窓スペクトルから互いに重複しない２つの低周波数範囲のスペクトル集合ＳＰ１ＵとＳＰ１Ｄを抽出するとともに、前記第３窓スペクトルから互いに重複しない２つの低周波数範囲のスペクトル集合ＳＰ３ＵとＳＰ３Ｄを抽出し、前記埋め込むべき付加情報の情報配列および当該情報配列の付加情報内の１単位における位置に基づいて、ＳＰ１Ｕ、ＳＰ１Ｄ、ＳＰ３Ｕ、ＳＰ３Ｄの割合を変更すると共に、中央の第２窓関数に対応する低周波成分を除去する低周波成分変更手段と、前記変更された低周波スペクトルを含む各窓スペクトルに対して周波数逆変換を行って、改変音響フレームを生成する周波数逆変換手段と、前記生成された改変音響フレームを順次出力する改変音響フレーム出力手段を有する音響信号に対する情報の埋め込み装置を提供する。
本発明によれば、音響信号の低周波領域から互いに重複しない２つの低周波成分を抽出するとともに、音響信号における所定区間の低周波成分を前部部分と後部部分に分け、所定区間の低周波領域において分けた４成分のスペクトル強度の割合を、埋め込むべき付加情報の情報配列および当該情報配列の付加情報内の１単位における位置に基づいて変更するようにしたので、従来と同等以上の抽出精度を維持しつつ、より多くの情報を埋め込むことが可能となるという効果を奏する。
以下、本発明の実施形態について図面を参照して詳細に説明する。
（１．本発明の基本概念）
まず、本発明の基本概念について、特許文献８、９に記載の発明と比較して説明する。本発明、特許文献８、９に記載の発明は、いずれも音響信号から所定数（例えば４０９６個）のサンプルを１音響フレームとして抽出し、この音響フレームに対して、周波数変換の後、その低周波成分の状態を、埋め込むべき付加情報の値に応じて変更するものである。図１に、従来の手法による低周波成分の状態を１音響フレーム単位で示す。図１に示す各音響フレームにおいて、横軸は時間方向、縦軸は周波数方向を示している。また、網掛けされた部分は、低周波成分が存在している部分を示し、網掛けが濃いほど成分強度が強いことを示している。
（１．本発明の基本概念）
まず、本発明の基本概念について、特許文献８、９に記載の発明と比較して説明する。本発明、特許文献８、９に記載の発明は、いずれも音響信号から所定数（例えば４０９６個）のサンプルを１音響フレームとして抽出し、この音響フレームに対して、周波数変換の後、その低周波成分の状態を、埋め込むべき付加情報の値に応じて変更するものである。図１に、従来の手法による低周波成分の状態を１音響フレーム単位で示す。図１に示す各音響フレームにおいて、横軸は時間方向、縦軸は周波数方向を示している。また、網掛けされた部分は、低周波成分が存在している部分を示し、網掛けが濃いほど成分強度が強いことを示している。
特許文献８、９に記載の発明については、図１（ａ）（ｅ）に示すような原状態の低周波成分に対して、符号０を埋め込む場合、Ｌ-ch（左チャンネル）の信号の低周波成分を全て除去し、除去した成分と同等の成分をＲ-ch（右チャンネル）の信号の低周波成分に加算して図１（ｂ）（ｆ）に示すような状態に変更するのは同じである。
特許文献８に記載の発明では、低周波成分のうちより高い成分とより低い成分の分布により埋め込む符号の区別を行う。したがって、符号１を埋め込む場合、図１（ｃ）に示すように、Ｌ-chの信号のより高い低周波成分を除去し、除去した成分と同等の成分をＲ-chの信号に加算し、逆に、Ｒ-chの信号のより低い低周波成分を除去し、除去した成分と同等の成分をＬ-chの信号に加算する。符号２を埋め込む場合は、図１（ｄ）に示すように、Ｌ-chの信号のより低い低周波成分を除去し、除去した成分と同等の成分をＲ-chの信号に加算し、逆に、Ｒ-chの信号のより高い低周波成分を除去し、除去した成分と同等の成分をＬ-chの信号に加算する。
特許文献９に記載の発明では、低周波成分のうち時間的に前部（過去）の成分と後部（未来）の成分の分布により埋め込む符号の区別を行う。したがって、符号１を埋め込む場合、図１（ｇ）に示すように、Ｌ-chの信号の後部の低周波成分を除去し、除去した成分と同等の成分をＲ-chの信号に加算し、逆に、Ｒ-chの信号の前部の低周波成分を除去し、除去した成分と同等の成分をＬ-chの信号に加算する。符号２を埋め込む場合は、図１（ｈ）に示すように、Ｌ-chの信号の前部の低周波成分を除去し、除去した成分と同等の成分をＲ-chの信号に加算し、逆に、Ｒ-chの信号の後部の低周波成分を除去し、除去した成分と同等の成分をＬ-chの信号に加算する。
特許文献８に記載の発明では、図１（ｂ）〜（ｄ）に示すような３つの状態に低周波成分を変更することにより、情報の埋め込みを行い、特許文献９に記載の発明では、図１（ｆ）〜（ｈ）に示すような３つの状態に低周波成分を変更することにより、情報の埋め込みを行っている。３つの状態のうち、１つは、埋め込むべき付加情報を構成するワード間の区切りやエラー処理等に用いるものであるので、実質埋め込むことができるのは２つの状態、すなわち、１ビットに相当する。したがって、特許文献８、９に記載の発明では、１音響フレームに対して１ビットの情報が埋め込み可能となる。
図２に、本発明による低周波成分の状態を１音響フレーム単位で示す。図２に示す各音響フレームにおいても、横軸は時間方向、縦軸は周波数方向を示している。本発明では、特許文献８、９に記載の発明のそれぞれの特徴を組み合わせ低周波成分のうちより高い成分とより低い成分の分布と、低周波成分のうち時間的に前部の成分と後部の成分の分布を組み合わせて埋め込む符号の区別を行う。
本発明においても、図２（ａ）に示すような原状態の低周波成分に対して、符号０を埋め込む場合、Ｌ-chの信号の低周波成分を全て除去し、除去した成分と同等の成分をＲ-chの信号の低周波成分に加算して図２（ｂ）に示すような状態に変更するのは、特許文献８、９に記載の発明と同じである。なお、図２においては、図１と異なり、縦軸の周波数方向が５つに区分されているが、これは、最下部にギャップ領域を設けたためである。このギャップ領域は、抽出対象であるＬ-chの信号においては、常に成分が除去され、除去された成分と同等の成分がＲ-chの信号の低周波成分に加算される。最下部にギャップ領域を設けることにより、埋め込み可能な低周波成分の上限を高くすることができる。さらに、５つに区分された低周波成分のうち最上部と上から３番目にもギャップ領域を設けている。上から２番目は上位低周波領域であり、上から４番目は下位低周波領域である。図２（ｂ）に示すように、上位低周波領域の前部におけるスペクトルをＳＰ１Ｕ、上位低周波領域の後部におけるスペクトルをＳＰ３Ｕ、下位低周波領域の前部におけるスペクトルをＳＰ１Ｄ、下位低周波領域の後部におけるスペクトルをＳＰ３Ｄで表現することとする。
本発明で、符号１を埋め込む場合、図２（ｃ）に示すように、Ｌ-chの信号の後部の低周波成分を除去し、除去した成分と同等の成分をＲ-chの信号に加算し、逆に、Ｒ-chの信号の前部の低周波成分を除去し、除去した成分と同等の成分をＬ-chの信号に加算する。符号２を埋め込む場合は、図２（ｄ）に示すように、Ｌ-chの信号の前部の低周波成分を除去し、除去した成分と同等の成分をＲ-chの信号に加算し、逆に、Ｒ-chの信号の後部の低周波成分を除去し、除去した成分と同等の成分をＬ-chの信号に加算する。
符号３を埋め込む場合、図２（ｅ）に示すように、Ｌ-chの信号の前部の上位低周波領域の成分である上位低周波成分と、後部の下位低周波領域の成分である下位低周波成分を除去し、除去した成分と同等の成分をＲ-chの信号に加算し、逆に、Ｒ-chの信号の前部の下位低周波成分と後部の上位低周波成分を除去し、除去した成分と同等の成分をＬ-chの信号に加算する。符号４を埋め込む場合は、図２（ｆ）に示すように、Ｌ-chの信号の前部の下位低周波成分と後部の上位低周波成分を除去し、除去した成分と同等の成分をＲ-chの信号に加算し、逆に、Ｒ-chの信号の前部の上位低周波成分と後部の下位低周波成分を除去し、除去した成分と同等の成分をＬ-chの信号に加算する。
符号５を埋め込む場合、図２（ｇ）に示すように、Ｌ-chの信号の前部の上位低周波成分と、後部の低周波成分を除去し、除去した成分と同等の成分をＲ-chの信号に加算し、逆に、Ｒ-chの信号の前部の下位低周波成分を除去し、除去した成分と同等の成分をＬ-chの信号に加算する。符号６を埋め込む場合は、図２（ｈ）に示すように、Ｌ-chの信号の前部の低周波成分と後部の上位低周波成分を除去し、除去した成分と同等の成分をＲ-chの信号に加算し、逆に、Ｒ-chの信号の後部の下位低周波成分を除去し、除去した成分と同等の成分をＬ-chの信号に加算する。
符号７を埋め込む場合、図２（ｉ）に示すように、Ｌ-chの信号の前部の低周波成分と、後部の下位低周波成分を除去し、除去した成分と同等の成分をＲ-chの信号に加算し、逆に、Ｒ-chの信号の後部の上位低周波成分を除去し、除去した成分と同等の成分をＬ-chの信号に加算する。符号８を埋め込む場合は、図２（ｊ）に示すように、Ｌ-chの信号の前部の下位低周波成分と後部の低周波成分を除去し、除去した成分と同等の成分をＲ-chの信号に加算し、逆に、Ｒ-chの信号の前部の上位低周波成分を除去し、除去した成分と同等の成分をＬ-chの信号に加算する。
本発明では、図２（ｂ）〜（ｊ）に示すような９つの状態に低周波成分を変更することにより、情報の埋め込みを行っている。９つの状態のうち、１つは、埋め込むことができない状態におけるエラー処理に用いるものであるので、実質埋め込むことができるのは８つの状態、すなわち、３ビットに相当する。ただし、３ビットのうち１ビットは埋め込むべき付加情報を構成するワード間の区切りを示すために用いるので、埋め込むべき付加情報については、本発明では、１音響フレームに対して２ビットの情報が埋め込み可能となる。ただし、区切りを示す情報を埋め込むためだけに１音響フレームを専有することはないので、全体としては付加情報の埋め込み効率が高まることになる。
（２．音響信号に対する情報の埋め込み装置）
図３は、本発明に係る音響信号に対する情報の埋め込み装置の構成を示す機能ブロック図である。図３において、１０は音響フレーム読込手段、２０は周波数変換手段、３０は低周波成分変更手段、４０は周波数逆変換手段、５０は改変音響フレーム出力手段、６０は記憶手段、６１は音響信号記憶部、６２は付加情報記憶部、６３は改変音響信号記憶部、７０は付加情報読込手段である。なお、図３に示す装置は、ステレオ音響信号、モノラル音響信号の両方に対応可能であるが、ここでは、ステレオ音響信号に対して処理を行う場合について説明していく。
図３は、本発明に係る音響信号に対する情報の埋め込み装置の構成を示す機能ブロック図である。図３において、１０は音響フレーム読込手段、２０は周波数変換手段、３０は低周波成分変更手段、４０は周波数逆変換手段、５０は改変音響フレーム出力手段、６０は記憶手段、６１は音響信号記憶部、６２は付加情報記憶部、６３は改変音響信号記憶部、７０は付加情報読込手段である。なお、図３に示す装置は、ステレオ音響信号、モノラル音響信号の両方に対応可能であるが、ここでは、ステレオ音響信号に対して処理を行う場合について説明していく。
音響フレーム読込手段１０は、付加情報の埋め込み対象とする元のステレオ音響信号の各チャンネルから所定数のサンプルを１フレームとして読み込む機能を有している。周波数変換手段２０は、音響フレーム読込手段１０が読み込んだ音響信号のフレームをフーリエ変換等により周波数変換してフレームスペクトルを生成する機能を有している。低周波成分変更手段３０は、生成されたフレームスペクトルから所定の低周波数範囲に相当するスペクトル集合を複数抽出し、付加情報記憶部６２から抽出した付加情報に基づいて、低周波強度データのスペクトル集合の状態を変更する機能を有している。周波数逆変換手段４０は、変更された低周波強度データを含む複数のフレームスペクトルに対して周波数逆変換を行うことにより、改変音響フレームを生成する機能を有している。改変音響フレーム出力手段５０は、生成された改変音響フレームを順次出力する機能を有している。記憶手段６０は、付加情報を埋め込む対象とするステレオ音響信号を記憶した音響信号記憶部６１と、ビット配列として構成され、ステレオ音響信号に埋め込まれる付加情報を記憶した付加情報記憶部６２と、付加情報埋め込み後の改変音響信号を記憶する改変音響信号記憶部６３を有しており、その他処理に必要な各種情報を記憶するものである。付加情報読込手段７０は、付加情報記憶部６２から付加情報を抽出する機能を有している。なお、付加情報とは、音響情報に付加して埋め込むべき情報であり、タイトルやアーティスト名等の属性情報、および属性情報以外の他の情報を含むものである。本実施形態では、付加情報は８ビットを１ワードとし、この１ワードを１単位として構成される。図３に示した各構成手段は、現実にはコンピュータおよびその周辺機器等のハードウェアに専用のプログラムを搭載することにより実現される。すなわち、コンピュータが、専用のプログラムに従って各手段の内容を実行することになる。
（３．埋め込み装置の処理動作）
次に、図３に示した音響信号に対する情報の埋め込み装置の処理動作について説明する。ここでは、音響信号として、Ｌ（左）、Ｒ（右）の２チャンネルを有するステレオ音響信号に対して処理を行う場合について説明していく。音響フレーム読込手段１０は、音響信号記憶部６１に記憶されたステレオ音響信号の左右の各チャンネルから、それぞれ所定数のサンプルを１音響フレームとして読み込む。音響フレーム読込手段１０が読み込む１音響フレームのサンプル数は、適宜設定することができるが、サンプリング周波数が４４．１ｋＨｚの場合、４０９６サンプル程度とすることが望ましい。したがって、音響フレーム読込手段１０は、左チャンネル、右チャンネルについてそれぞれ４０９６サンプルずつ、順次音響フレームとして読み込んでいくことになる。
次に、図３に示した音響信号に対する情報の埋め込み装置の処理動作について説明する。ここでは、音響信号として、Ｌ（左）、Ｒ（右）の２チャンネルを有するステレオ音響信号に対して処理を行う場合について説明していく。音響フレーム読込手段１０は、音響信号記憶部６１に記憶されたステレオ音響信号の左右の各チャンネルから、それぞれ所定数のサンプルを１音響フレームとして読み込む。音響フレーム読込手段１０が読み込む１音響フレームのサンプル数は、適宜設定することができるが、サンプリング周波数が４４．１ｋＨｚの場合、４０９６サンプル程度とすることが望ましい。したがって、音響フレーム読込手段１０は、左チャンネル、右チャンネルについてそれぞれ４０９６サンプルずつ、順次音響フレームとして読み込んでいくことになる。
音響フレームとしては、ＡタイプとＢタイプが存在する。Ａタイプの音響フレーム、Ｂタイプの音響フレームは、それぞれ同タイプの先行する音響フレームの最後のサンプルの次のサンプルを先頭サンプルとして設定される。そして、ＡタイプとＢタイプの音響フレームは互いに所定数（本実施形態では２０４８）のサンプルを重複して設定される。例えば、Ａタイプの音響フレームを先頭からＡ１、Ａ２、Ａ３…とし、Ｂタイプの音響フレームを先頭からＢ１、Ｂ２、Ｂ３…とすると、Ａ１はサンプル１〜４０９６、Ａ２はサンプル４０９７〜８１９２、Ａ３はサンプル８１９３〜１２２８８、Ｂ１はサンプル２０４９〜６１４４、Ｂ２はサンプル６１４５〜１０２４０、Ｂ３はサンプル１０２４１〜１４３３６となる。なお、ＡタイプとＢタイプは相対的なものであるので、どちらが先であっても良い。すなわち、上記とは逆にＡ１がサンプル２０４９〜６１４４、Ａ２がサンプル６１４５〜１０２４０、Ａ３がサンプル１０２４１〜１４３３６、Ｂ１がサンプル１〜４０９６、Ｂ２がサンプル４０９７〜８１９２、Ｂ３がサンプル８１９３〜１２２８８であっても良い。
周波数変換手段２０は、音響フレーム読込手段１０が読み込んだ音響フレームに対して、周波数変換を行って、その音響フレームのスペクトルであるフレームスペクトルを得る。具体的には、窓関数を利用して周波数変換を行う。周波数変換としては、フーリエ変換、ウェーブレット変換その他公知の種々の手法を用いることができる。本実施形態では、フーリエ変換を用いた場合を例にとって説明する。
一般に、所定の信号に対してフーリエ変換を行う場合、信号を所定の長さに区切って行う必要があるが、この場合、所定長さの信号に対してそのままフーリエ変換を行うと、擬似高調波成分が発生する。そこで、一般にフーリエ変換を行う場合には、ハニング窓と呼ばれる窓関数を用いて、信号の値を変化させた後、変化後の値に対してフーリエ変換を実行する。
本発明では、擬似高調波成分の発生を防止するためだけではなく、特許文献９に記載の発明と同様、一つの音響フレームから、情報を埋め込むための複数の状態を作り出すため、複数の窓関数を用意し、一つの音響フレームに対して、各窓関数を利用してフーリエ変換を行い、複数のスペクトルを得る。複数の窓関数として、特許文献９に示すものと同一のものを用いても本発明による効果は得られるが、より良い効果を得るため、本発明では、図４（ｂ）〜（ｄ）に示したような第１窓関数Ｗ（１，ｉ）、第２窓関数Ｗ（２，ｉ）、第３窓関数Ｗ（３，ｉ）を用意し、抽出側で認識し易いようにした。第１窓関数Ｗ（１，ｉ）は、音響フレームの前部を抽出するためのものであり、図４（ｂ）に示すように前部の所定のサンプル番号ｉの位置において、最大値１をとり、後部においては、最小値０をとるように設定されている。どのサンプル番号の場合に最大値をとるかについては、窓関数Ｗ（１，ｉ）の設計によって異なってくるが、本実施形態では、後述する〔数式１〕で定義される。窓関数Ｗ（１，ｉ）を乗じることにより、図４（ａ）に示すような音響フレームの信号波形は、図４（ｆ）に示すように、前部に信号成分が残り、後部の信号成分が削除されたものとなり、これがフーリエ変換対象となる。
また、第２窓関数Ｗ（２，ｉ）は、音響フレームの中央部を抽出するためのものであり、図４（ｃ）に示すように、中央部の所定のサンプル番号ｉの位置において、最大値１をとり、前部、後部においては、最小値０をとるように設定されている。どのサンプル番号の場合に最大値をとるかについては、窓関数Ｗ（２，ｉ）の設計によって異なってくるが、本実施形態では、後述する〔数式２〕で定義される。窓関数Ｗ（２，ｉ）を乗じることにより、図４（ａ）に示すような音響フレームの信号波形は、図４（ｇ）に示すように、中央部に信号成分が残り、前部と後部の信号成分が除去されたものとなり、これがフーリエ変換対象となる。
また、第３窓関数Ｗ（３，ｉ）は、音響フレームの後部を抽出するためのものであり、図４（ｄ）に示すように、前部においては最小値０をとり、後部の所定のサンプル番号ｉの位置において、最大値１をとるように設定されている。どのサンプル番号の場合に最大値をとるかについては、窓関数Ｗ（３，ｉ）の設計によって異なってくるが、本実施形態では、後述する〔数式３〕で定義される。窓関数Ｗ（３，ｉ）を乗じることにより、図４（ａ）に示すような音響フレームの信号波形は、図４（ｈ）に示すように、前部の信号成分が除去され、後部に信号成分が残ったものとなり、これがフーリエ変換対象となる。このように前部、中央部、後部を抽出した後、フーリエ変換を実行するため、前部、中央部、後部に対応したスペクトルが得られることになる。１つの音響フレームにビット値を埋め込むためには、本来、前部と後部の２つに分けられれば良いのであるが、抽出側においては、必ずしも、信号を同期して読み込むことができるとは限らず、したがって、前部と後部をはっきりと区別するため、本発明では、埋め込み時に中央部の信号成分を常に削除し、前部と後部を時間的に分離することとしている（ただし、抽出時は前部と後部だけを解析すればよく、中央部は無視してよい）。本発明において用いる窓関数は、窓関数Ｗ（１，ｉ）と窓関数Ｗ（３，ｉ）が左右非対称であるため、抽出時において、埋め込まれていた情報の誤認識が起こりにくくなる。
また、本発明では、音響フレームを重複させて読み込み、奇数フレーム（または偶数フレーム）については、窓関数Ｗ（１，ｉ）、Ｗ（２，ｉ）、Ｗ（３，ｉ）を用い、偶数フレーム（または奇数フレーム）については、図４（ｅ）に示したような第４窓関数Ｗ（４，ｉ）を用いるようにした。
なお、本発明においては、音響フレームは重複して読み込まれる。すなわち、奇数番目の音響フレームと偶数番目の音響フレームは、所定数のサンプルを重複して読み込む。上記のように、奇数フレームと偶数フレームでは、用いられる窓関数が異なるが、奇数フレームと偶数フレームは単に奇数か偶数かの違いだけであるため、どちらに対して処理をどちらの処理を行っても良い。したがって、本明細書では、奇数フレーム、偶数フレームの一方をＡタイプフレーム、他方をＢタイプフレームと呼ぶことにする。本実施形態では、奇数フレームをＡタイプフレーム、偶数フレームをＢタイプフレームとして説明するが、逆に偶数フレームをＡタイプフレーム、奇数フレームをＢタイプフレームとしても良い。
本実施形態では、窓関数Ｗ（１，ｉ）〜Ｗ（４，ｉ）は、以下の〔数式１〕〜〔数式４〕で定義される。なお、図４において、横軸は時間軸（ｉ）である。ｉは、後述するように、各音響フレーム内のＮ個のサンプルに付した通し番号であるため時刻ｔに比例している。また、図４（ａ）（ｆ）（ｇ）（ｈ）（ｉ）において縦軸は信号の振幅値（レベル）を示す。図４（ｂ）〜（ｅ）において縦軸は窓関数Ｗ（１，ｉ）、Ｗ（２，ｉ）、Ｗ（３，ｉ）、Ｗ（４，ｉ）の値を示しており、Ｗ（１，ｉ）、Ｗ（２，ｉ）、Ｗ（３，ｉ）、Ｗ（４，ｉ）の最大値はいずれも１である。
〔数式１〕
ｉ≦３Ｎ／８のとき、Ｗ（１，ｉ）＝０．５−０．５ｃｏｓ（８πｉ／（３Ｎ））
３Ｎ／８＜ｉ≦Ｎ／２のとき、Ｗ（１，ｉ）＝０．５−０．５ｃｏｓ（８π（ｉ−Ｎ／４）／Ｎ）
ｉ＞Ｎ／２のとき、Ｗ（１，ｉ）＝０．０
ｉ≦３Ｎ／８のとき、Ｗ（１，ｉ）＝０．５−０．５ｃｏｓ（８πｉ／（３Ｎ））
３Ｎ／８＜ｉ≦Ｎ／２のとき、Ｗ（１，ｉ）＝０．５−０．５ｃｏｓ（８π（ｉ−Ｎ／４）／Ｎ）
ｉ＞Ｎ／２のとき、Ｗ（１，ｉ）＝０．０
〔数式２〕
ｉ≦３Ｎ／８のとき、Ｗ（２，ｉ）＝０．０
３Ｎ／８＜ｉ≦Ｎ／２のとき、Ｗ（２，ｉ）＝０．５−０．５ｃｏｓ（４π（ｉ−Ｎ／４）／Ｎ）
ｉ＞３Ｎ／４のとき、Ｗ（２，ｉ）＝０．０
ｉ≦３Ｎ／８のとき、Ｗ（２，ｉ）＝０．０
３Ｎ／８＜ｉ≦Ｎ／２のとき、Ｗ（２，ｉ）＝０．５−０．５ｃｏｓ（４π（ｉ−Ｎ／４）／Ｎ）
ｉ＞３Ｎ／４のとき、Ｗ（２，ｉ）＝０．０
〔数式３〕
ｉ≦Ｎ／２のとき、Ｗ（３，ｉ）＝０．０
ｉ＞Ｎ／２のとき、Ｗ（３，ｉ）＝０．５−０．５ｃｏｓ（４π（ｉ−Ｎ／２）／Ｎ）
ｉ≦Ｎ／２のとき、Ｗ（３，ｉ）＝０．０
ｉ＞Ｎ／２のとき、Ｗ（３，ｉ）＝０．５−０．５ｃｏｓ（４π（ｉ−Ｎ／２）／Ｎ）
〔数式４〕
ｉ≦Ｎ／４のとき、Ｗ（４，ｉ）＝０．０
Ｎ／４＜ｉ≦Ｎ／２のとき、Ｗ（４，ｉ）＝０．５−０．５ｃｏｓ（４π（ｉ−Ｎ／４）／Ｎ）
Ｎ／２＜ｉ≦７Ｎ／８のとき、Ｗ（４，ｉ）＝０．５−０．５ｃｏｓ（８π（ｉ−Ｎ／８）／（３Ｎ））
ｉ＞７Ｎ／８のとき、Ｗ（４，ｉ）＝０．０
ｉ≦Ｎ／４のとき、Ｗ（４，ｉ）＝０．０
Ｎ／４＜ｉ≦Ｎ／２のとき、Ｗ（４，ｉ）＝０．５−０．５ｃｏｓ（４π（ｉ−Ｎ／４）／Ｎ）
Ｎ／２＜ｉ≦７Ｎ／８のとき、Ｗ（４，ｉ）＝０．５−０．５ｃｏｓ（８π（ｉ−Ｎ／８）／（３Ｎ））
ｉ＞７Ｎ／８のとき、Ｗ（４，ｉ）＝０．０
なお、図４および上記〔数式１〕〜〔数式４〕から明らかなように、窓関数Ｗ（１，ｉ）とＷ（３，ｉ）は、互いに非対称な形状である。これは、後述する抽出側において、両者の識別を容易にするためである。また、窓関数Ｗ（１，ｉ）、Ｗ（２，ｉ）、Ｗ（３，ｉ）は、ｉが所定の値のときに最大値１をとり、ｉがその他の値をとる場合には、ｉの値に応じて単調増加、または単調減少する窓関数を分割したものであるため、窓関数Ｗ（１，ｉ）とＷ（３，ｉ）が定まると、窓関数Ｗ（２，ｉ）も必然的に定まる。このため、窓関数Ｗ（２，ｉ）は左右非対称の形状となっている。
本発明においては、奇数フレームと偶数フレームを、所定サンプルずつ重複して読み込むため、情報の埋め込みを行った後、音響信号に復元する際に、窓関数を乗じた奇数フレームと、窓関数を乗じた偶数フレームの重複サンプルを加算した場合に、ほぼ元の値に戻るようにしなければならない。このため、窓関数Ｗ（４，ｉ）の形状は、窓関数Ｗ（１，ｉ）、Ｗ（２，ｉ）、Ｗ（３，ｉ）の値に応じて必然的に定まる。すなわち、奇数フレームと偶数フレームの重複部分において、窓関数Ｗ（１，ｉ）、Ｗ（２，ｉ）、Ｗ（３，ｉ）、Ｗ（４，ｉ）を加算すると、全区間固定値１になるように定義されている。
周波数変換手段２０が、Ａタイプの音響フレームに対してフーリエ変換を行う場合は、左チャンネル信号Ｘｌ（ｉ）、右チャンネル信号Ｘｒ（ｉ）（ｉ＝０，…，Ｎ−１）に対して、３つの窓関数Ｗ（１，ｉ）、Ｗ（２，ｉ）、Ｗ（３，ｉ）を用いて、以下の〔数式５〕に従った処理を行い、左チャンネルに対応する変換データの実部Ａｌ（１，ｊ）、Ａｌ（２，ｊ）、Ａｌ（３，ｊ）、虚部Ｂｌ（１，ｊ）、Ｂｌ（２，ｊ）、Ｂｌ（３，ｊ）、右チャンネルに対応する変換データの実部Ａｒ（１，ｊ）、Ａｒ（２，ｊ）、Ａｒ（３，ｊ）、虚部Ｂｒ（１，ｊ）、Ｂｒ（２，ｊ）、Ｂｒ（３，ｊ）を得る。なお、窓関数Ｗ（１，ｉ）、Ｗ（２，ｉ）、Ｗ（３，ｉ）は、それぞれ音響フレームの前部（先頭）付近、中央付近、後部付近において値が大きくなる関数となっている。
〔数式５〕
Ａｌ（１，ｊ）＝Σi=0,…,N-1Ｗ（１，ｉ）・Ｘｌ（ｉ）・ｃｏｓ（２πｉｊ／Ｎ）
Ｂｌ（１，ｊ）＝Σi=0,…,N-1Ｗ（１，ｉ）・Ｘｌ（ｉ）・ｓｉｎ（２πｉｊ／Ｎ）
Ａｌ（２，ｊ）＝Σi=0,…,N-1Ｗ（２，ｉ）・Ｘｌ（ｉ）・ｃｏｓ（２πｉｊ／Ｎ）
Ｂｌ（２，ｊ）＝Σi=0,…,N-1Ｗ（２，ｉ）・Ｘｌ（ｉ）・ｓｉｎ（２πｉｊ／Ｎ）
Ａｌ（３，ｊ）＝Σi=0,…,N-1Ｗ（３，ｉ）・Ｘｌ（ｉ）・ｃｏｓ（２πｉｊ／Ｎ）
Ｂｌ（３，ｊ）＝Σi=0,…,N-1Ｗ（３，ｉ）・Ｘｌ（ｉ）・ｓｉｎ（２πｉｊ／Ｎ）
Ａｒ（１，ｊ）＝Σi=0,…,N-1Ｗ（１，ｉ）・Ｘｒ（ｉ）・ｃｏｓ（２πｉｊ／Ｎ）
Ｂｒ（１，ｊ）＝Σi=0,…,N-1Ｗ（１，ｉ）・Ｘｒ（ｉ）・ｓｉｎ（２πｉｊ／Ｎ）
Ａｒ（２，ｊ）＝Σi=0,…,N-1Ｗ（２，ｉ）・Ｘｒ（ｉ）・ｃｏｓ（２πｉｊ／Ｎ）
Ｂｒ（２，ｊ）＝Σi=0,…,N-1Ｗ（２，ｉ）・Ｘｒ（ｉ）・ｓｉｎ（２πｉｊ／Ｎ）
Ａｒ（３，ｊ）＝Σi=0,…,N-1Ｗ（３，ｉ）・Ｘｒ（ｉ）・ｃｏｓ（２πｉｊ／Ｎ）
Ｂｒ（３，ｊ）＝Σi=0,…,N-1Ｗ（３，ｉ）・Ｘｒ（ｉ）・ｓｉｎ（２πｉｊ／Ｎ）
Ａｌ（１，ｊ）＝Σi=0,…,N-1Ｗ（１，ｉ）・Ｘｌ（ｉ）・ｃｏｓ（２πｉｊ／Ｎ）
Ｂｌ（１，ｊ）＝Σi=0,…,N-1Ｗ（１，ｉ）・Ｘｌ（ｉ）・ｓｉｎ（２πｉｊ／Ｎ）
Ａｌ（２，ｊ）＝Σi=0,…,N-1Ｗ（２，ｉ）・Ｘｌ（ｉ）・ｃｏｓ（２πｉｊ／Ｎ）
Ｂｌ（２，ｊ）＝Σi=0,…,N-1Ｗ（２，ｉ）・Ｘｌ（ｉ）・ｓｉｎ（２πｉｊ／Ｎ）
Ａｌ（３，ｊ）＝Σi=0,…,N-1Ｗ（３，ｉ）・Ｘｌ（ｉ）・ｃｏｓ（２πｉｊ／Ｎ）
Ｂｌ（３，ｊ）＝Σi=0,…,N-1Ｗ（３，ｉ）・Ｘｌ（ｉ）・ｓｉｎ（２πｉｊ／Ｎ）
Ａｒ（１，ｊ）＝Σi=0,…,N-1Ｗ（１，ｉ）・Ｘｒ（ｉ）・ｃｏｓ（２πｉｊ／Ｎ）
Ｂｒ（１，ｊ）＝Σi=0,…,N-1Ｗ（１，ｉ）・Ｘｒ（ｉ）・ｓｉｎ（２πｉｊ／Ｎ）
Ａｒ（２，ｊ）＝Σi=0,…,N-1Ｗ（２，ｉ）・Ｘｒ（ｉ）・ｃｏｓ（２πｉｊ／Ｎ）
Ｂｒ（２，ｊ）＝Σi=0,…,N-1Ｗ（２，ｉ）・Ｘｒ（ｉ）・ｓｉｎ（２πｉｊ／Ｎ）
Ａｒ（３，ｊ）＝Σi=0,…,N-1Ｗ（３，ｉ）・Ｘｒ（ｉ）・ｃｏｓ（２πｉｊ／Ｎ）
Ｂｒ（３，ｊ）＝Σi=0,…,N-1Ｗ（３，ｉ）・Ｘｒ（ｉ）・ｓｉｎ（２πｉｊ／Ｎ）
周波数変換手段２０が、Ｂタイプの音響フレームに対してフーリエ変換を行う場合は、左チャンネル信号Ｘｌ（ｉ）、右チャンネル信号Ｘｒ（ｉ）（ｉ＝０，…，Ｎ−１）に対して、窓関数Ｗ（４，ｉ）を用いて、以下の〔数式６〕に従った処理を行い、左チャンネルに対応する変換データの実部Ａｌ（４，ｊ）、虚部Ｂｌ（４，ｊ）、右チャンネルに対応する変換データの実部Ａｒ（４，ｊ）、虚部Ｂｒ（４，ｊ）を得る。
〔数式６〕
Ａｌ（４，ｊ）＝Σi=0,…,N-1Ｗ（４，ｉ）・Ｘｌ（ｉ）・ｃｏｓ（２πｉｊ／Ｎ）
Ｂｌ（４，ｊ）＝Σi=0,…,N-1Ｗ（４，ｉ）・Ｘｌ（ｉ）・ｓｉｎ（２πｉｊ／Ｎ）
Ａｒ（４，ｊ）＝Σi=0,…,N-1Ｗ（４，ｉ）・Ｘｒ（ｉ）・ｃｏｓ（２πｉｊ／Ｎ）
Ｂｒ（４，ｊ）＝Σi=0,…,N-1Ｗ（４，ｉ）・Ｘｒ（ｉ）・ｓｉｎ（２πｉｊ／Ｎ）
Ａｌ（４，ｊ）＝Σi=0,…,N-1Ｗ（４，ｉ）・Ｘｌ（ｉ）・ｃｏｓ（２πｉｊ／Ｎ）
Ｂｌ（４，ｊ）＝Σi=0,…,N-1Ｗ（４，ｉ）・Ｘｌ（ｉ）・ｓｉｎ（２πｉｊ／Ｎ）
Ａｒ（４，ｊ）＝Σi=0,…,N-1Ｗ（４，ｉ）・Ｘｒ（ｉ）・ｃｏｓ（２πｉｊ／Ｎ）
Ｂｒ（４，ｊ）＝Σi=0,…,N-1Ｗ（４，ｉ）・Ｘｒ（ｉ）・ｓｉｎ（２πｉｊ／Ｎ）
上記〔数式５〕〔数式６〕において、ｉは、各音響フレーム内のＮ個のサンプルに付した通し番号であり、ｉ＝０，１，２，…Ｎ−１の整数値をとる。また、ｊは周波数の値について、値の小さなものから順に付した通し番号であり、ｉと同様にｊ＝０，１，２，…Ｎ／２−１の整数値をとる。サンプリング周波数が４４．１ｋＨｚ、Ｎ＝４０９６の場合、ｊの値が１つ異なると、周波数が１０．８Ｈｚ異なることになる。
上記〔数式５〕〔数式６〕に従った処理を実行することにより、各音響フレームの信号成分を周波数に対応した成分であるスペクトルで表現されたフレームスペクトルが得られる。続いて、低周波成分変更手段３０が、生成されたフレームスペクトルから３つの所定周波数範囲のスペクトル集合を抽出する。人間の聴覚は、２００〜３００Ｈｚ程度までの低周波成分については、方向性を感知しにくくなっていることが知られている（コロナ社１９９０年１０月３０日発行「音響工学講座１．基礎音響工学、日本音響学会編」ｐ２４７図９・２６参照）。したがって、特許文献７〜９に記載の発明においては、低周波成分を２００〜３００Ｈｚ程度以下としていた。しかし、本出願人による研究の結果、最低周波数帯における成分を削除した場合は、４００Ｈｚ付近においても方向性を感知しにくくなることがわかった。そこで、本実施形態では、０〜１５０Ｈｚ付近の最低周波数帯の成分を除去し、１５０〜４００Ｈｚ程度の低周波数帯に情報を埋め込むことにしている。さらに、本実施形態では、埋め込み対象とする低周波数帯を比較的高い低周波数帯と比較的低い低周波数帯に分けて、それぞれに対して加工を行うようにしている。
ここで、本発明の処理による周波数成分全体の状況を図５に示す。図５において、上下方向は周波数方向であり、ｊ＝２０４７は２２．０５ｋＨｚに対応している。図５（ａ）は特許文献９に示した従来方式、図５（ｂ）は本発明による方式を示す。図５（ｂ）は、最低周波数帯をｊ＝０〜１４（０〜約１５０Ｈｚ）、埋め込み対象とする低周波数帯のうち比較的高い低周波数帯をｊ＝２７〜３５（約２９０〜３８０Ｈｚ）、埋め込み対象とする低周波数帯のうち比較的低い低周波数帯をｊ＝１５〜２３（約１６０〜２５０Ｈｚ）とした場合を示している。また、網掛けの濃さは、信号成分の強度を概念的に示したものである。従来方式も本発明もｊ＝４０（従来はｊ＝３０）以上の高周波成分に対しては、変更を加えず、原音を維持する点については同じである。また、原音部との境を明確にするために、数成分（本実施形態では、ｊ＝３６〜３９の４成分）を左（Ｌ）信号から右（Ｒ）信号に移動させているのも、従来方式と同じである。
本発明が従来方式と異なるのは、従来方式がｊ＝０を原音維持とし、ｊ＝１〜２９の成分を原則左（Ｌ）信号から右（Ｒ）信号に移動させ、ｊ＝１〜２６の成分に対しては付加情報に応じて選択的に右（Ｒ）信号から左（Ｌ）信号に移動させることにより、ｊ＝１〜２６を埋め込み領域としているのに対して、本発明ではｊ＝０〜３９の成分を原則左（Ｌ）信号から右（Ｒ）信号に移動させ、ｊ＝１５〜２３、２７〜３５の成分に対しては付加情報に応じて選択的に右（Ｒ）信号から左（Ｌ）信号に移動させることにより、ｊ＝１５〜３５を埋め込み領域としている点である。なお、図５中のｊ＝１５〜２３、２７〜３５の領域については、埋め込まれる情報により、ＬＲのどちらかに成分が移動されることを示している。
低周波成分変更手段３０は、Ａタイプの音響フレームについて、付加情報読込手段７０が読み込んだビット値に応じて、低周波成分の割合を変更する処理を行う。本発明では、付加情報を２ビットずつ読み込み、１音響フレームに対して２ビットの情報を埋め込む。埋め込まれる２ビットの値は、“００”“０１”“１０”“１１”の４通りがある。本実施形態では、これらを値１〜値４と定義する。この際、“００”“０１”“１０”“１１”の４通りのうち、いずれを値１〜値４と定義しても良い。抽出側において、埋め込み側で埋め込まれた２ビットの配列が特定できれば良いためである。したがって、この定義は、埋め込み側と抽出側で一致している必要がある。
また、本発明では、区切りを示す情報を付加情報と別に埋め込むようにしていないため、１ワードを構成する最後の２ビットを埋め込む際に、１ワードの終わりであることを示す情報を追加する。そのため、１ワード中の最後の２ビットを埋め込む際には、値１〜値４にさらに値４を追加し、値５〜値８とする処理を行う。そして、低周波成分変更手段３０は、この値１〜値８に応じて、低周波成分を状態１〜状態８のいずれかに変更する。具体的には、まず、低周波成分変更手段３０は、上記〔数式５〕により求めた左チャンネルの実部Ａｌ（１，ｊ）、Ａｌ（３，ｊ）、虚部Ｂｌ（１，ｊ）、Ｂｌ（３，ｊ）、右チャンネルの実部Ａｒ（１，ｊ）、Ａｒ（３，ｊ）、虚部Ｂｒ（１，ｊ）、Ｂｒ（３，ｊ）を利用して、以下の〔数式７〕により、合算値Ｅ1〜Ｅ4を算出する。
〔数式７〕
Ｅ1＝Σj=m,m+P-4｛Ａｌ（１，ｊ）2＋Ｂｌ（１，ｊ）2＋Ａｒ（１，ｊ）2＋Ｂｒ（１，ｊ）2｝
Ｅ2＝Σj=m,m+P-4｛Ａｌ（３，ｊ）2＋Ｂｌ（３，ｊ）2＋Ａｒ（３，ｊ）2＋Ｂｒ（３，ｊ）2｝
Ｅ3＝Σj=m+p,m+2P-4｛Ａｌ（１，ｊ）2＋Ｂｌ（１，ｊ）2＋Ａｒ（１，ｊ）2＋Ｂｒ（１，ｊ）2｝
Ｅ4＝Σj=m+p,m+2P-4｛Ａｌ（３，ｊ）2＋Ｂｌ（３，ｊ）2＋Ａｒ（３，ｊ）2＋Ｂｒ（３，ｊ）2｝
Ｅ1＝Σj=m,m+P-4｛Ａｌ（１，ｊ）2＋Ｂｌ（１，ｊ）2＋Ａｒ（１，ｊ）2＋Ｂｒ（１，ｊ）2｝
Ｅ2＝Σj=m,m+P-4｛Ａｌ（３，ｊ）2＋Ｂｌ（３，ｊ）2＋Ａｒ（３，ｊ）2＋Ｂｒ（３，ｊ）2｝
Ｅ3＝Σj=m+p,m+2P-4｛Ａｌ（１，ｊ）2＋Ｂｌ（１，ｊ）2＋Ａｒ（１，ｊ）2＋Ｂｒ（１，ｊ）2｝
Ｅ4＝Σj=m+p,m+2P-4｛Ａｌ（３，ｊ）2＋Ｂｌ（３，ｊ）2＋Ａｒ（３，ｊ）2＋Ｂｒ（３，ｊ）2｝
上記〔数式７〕において、ｍは低周波数帯の下限の成分の番号、Ｍは低周波数帯の上限の成分の番号であり、Ｐ＝（Ｍ−ｍ）／２である。なお、Ｐが整数にならない場合、切捨て処理を行う。ｍ＝１５、Ｍ＝４０と設定すると、Ｐ＝１２となり、図５に例示したスケールと同じになる。上記〔数式７〕により算出されたＥ1〜Ｅ4は音響フレームの所定の範囲におけるスペクトル集合の成分強度の合算値を示すことになる。続いて、この合算値Ｅ1〜Ｅ4がレベル下限値Ｌｅｖより大きいかどうかの判定を行う。レベル下限値Ｌｅｖは、音響信号Ｘｌ（ｉ）、Ｘｒ（ｉ）の振幅最大値が１に正規化されており、Ｍ＝４０に設定されている場合、０．２５に設定する。このＬｅｖ＝０．２５という値は、経験的にアナログ変換への耐性が維持できるレベルであり、低周波成分が少ない場合は適宜下げることになるが、その場合は、アナログ変換により検出精度も低下することになる。
合算値Ｅ1〜Ｅ4がレベル下限値Ｌｅｖより大きいかどうかを判断するのは、信号の強度が小さいと、信号を変化させても、その変化を抽出側で検出することができないためである。すなわち、各低周波スペクトル成分を十分な大きさとすることができるかどうかを判断することになる。
Ｅ1およびＥ3が共にレベル下限値Ｌｅｖより大きく、埋め込むべき情報が“値１”である場合、以下の〔数式８〕に従った処理を実行することにより、低周波成分の状態を“状態１”、すなわち、図２（ｃ）に示したような状態に変更する。
〔数式８〕
ｊ＝ｍ〜ｍ＋Ｐ−４、ｍ＋Ｐ〜ｍ＋２Ｐ−４の各成分に対して
Ａｌ´（３，ｊ）＝０
Ｂｌ´（３，ｊ）＝０
Ｅ（３，ｊ）＝｛Ａｌ（３，ｊ）2＋Ｂｌ（３，ｊ）2＋Ａｒ（３，ｊ）2＋Ｂｒ（３，ｊ）2｝1/2
Ａｒ´（３，ｊ）＝Ａｒ（３，ｊ）・Ｅ（３，ｊ）／｛Ａｒ（３，ｊ）2＋Ｂｒ（３，ｊ）2｝1/2
Ｂｒ´（３，ｊ）＝Ｂｒ（３，ｊ）・Ｅ（３，ｊ）／｛Ａｒ（３，ｊ）2＋Ｂｒ（３，ｊ）2｝1/2
Ａｒ´（１，ｊ）＝０
Ｂｒ´（１，ｊ）＝０
Ｅ（１，ｊ）＝｛Ａｌ（１，ｊ）2＋Ｂｌ（１，ｊ）2＋Ａｒ（１，ｊ）2＋Ｂｒ（１，ｊ）2｝1/2
Ａｌ´（１，ｊ）＝Ａｌ（１，ｊ）・Ｅ（１，ｊ）／｛Ａｌ（１，ｊ）2＋Ｂｌ（１，ｊ）2｝1/2
Ｂｌ´（１，ｊ）＝Ｂｌ（１，ｊ）・Ｅ（１，ｊ）／｛Ａｌ（１，ｊ）2＋Ｂｌ（１，ｊ）2｝1/2
ｊ＝ｍ〜ｍ＋Ｐ−４、ｍ＋Ｐ〜ｍ＋２Ｐ−４の各成分に対して
Ａｌ´（３，ｊ）＝０
Ｂｌ´（３，ｊ）＝０
Ｅ（３，ｊ）＝｛Ａｌ（３，ｊ）2＋Ｂｌ（３，ｊ）2＋Ａｒ（３，ｊ）2＋Ｂｒ（３，ｊ）2｝1/2
Ａｒ´（３，ｊ）＝Ａｒ（３，ｊ）・Ｅ（３，ｊ）／｛Ａｒ（３，ｊ）2＋Ｂｒ（３，ｊ）2｝1/2
Ｂｒ´（３，ｊ）＝Ｂｒ（３，ｊ）・Ｅ（３，ｊ）／｛Ａｒ（３，ｊ）2＋Ｂｒ（３，ｊ）2｝1/2
Ａｒ´（１，ｊ）＝０
Ｂｒ´（１，ｊ）＝０
Ｅ（１，ｊ）＝｛Ａｌ（１，ｊ）2＋Ｂｌ（１，ｊ）2＋Ａｒ（１，ｊ）2＋Ｂｒ（１，ｊ）2｝1/2
Ａｌ´（１，ｊ）＝Ａｌ（１，ｊ）・Ｅ（１，ｊ）／｛Ａｌ（１，ｊ）2＋Ｂｌ（１，ｊ）2｝1/2
Ｂｌ´（１，ｊ）＝Ｂｌ（１，ｊ）・Ｅ（１，ｊ）／｛Ａｌ（１，ｊ）2＋Ｂｌ（１，ｊ）2｝1/2
上記〔数式８〕においては、ｊ＝ｍ〜ｍ＋Ｐ−４、ｍ＋Ｐ〜ｍ＋２Ｐ−４において、Ａｌ´（３，ｊ）、Ｂｌ´（３，ｊ）を共に０としている。これは、図２（ｃ）の上段に示すように、L-chにおいて、SP3U、SP3D内の各成分を０にすることを示しているが、“状態１”はSP1U、SP1Ｄとの差を明確にすることができれば十分であるため、必ずしも０にする必要はなく、小さな値であれば良い。
Ｅ2およびＥ4が共にレベル下限値Ｌｅｖより大きく、埋め込むべき情報が“値２”である場合は、以下の〔数式９〕に従った処理を実行することにより、低周波成分の状態を“状態２” 、すなわち、図２（ｄ）に示したような状態に変更する。
〔数式９〕
ｊ＝ｍ〜ｍ＋Ｐ−４、ｍ＋Ｐ〜ｍ＋２Ｐ−４の各成分に対して
Ａｌ´（１，ｊ）＝０
Ｂｌ´（１，ｊ）＝０
Ｅ（１，ｊ）＝｛Ａｌ（１，ｊ）2＋Ｂｌ（１，ｊ）2＋Ａｒ（１，ｊ）2＋Ｂｒ（１，ｊ）2｝1/2
Ａｒ´（１，ｊ）＝Ａｒ（１，ｊ）・Ｅ（１，ｊ）／｛Ａｒ（１，ｊ）2＋Ｂｒ（１，ｊ）2｝1/2
Ｂｒ´（１，ｊ）＝Ｂｒ（１，ｊ）・Ｅ（１，ｊ）／｛Ａｒ（１，ｊ）2＋Ｂｒ（１，ｊ）2｝1/2
Ａｒ´（３，ｊ）＝０
Ｂｒ´（３，ｊ）＝０
Ｅ（３，ｊ）＝｛Ａｌ（３，ｊ）2＋Ｂｌ（３，ｊ）2＋Ａｒ（３，ｊ）2＋Ｂｒ（３，ｊ）2｝1/2
Ａｌ´（３，ｊ）＝Ａｌ（３，ｊ）・Ｅ（３，ｊ）／｛Ａｌ（３，ｊ）2＋Ｂｌ（３，ｊ）2｝1/2
Ｂｌ´（３，ｊ）＝Ｂｌ（３，ｊ）・Ｅ（３，ｊ）／｛Ａｌ（３，ｊ）2＋Ｂｌ（３，ｊ）2｝1/2
ｊ＝ｍ〜ｍ＋Ｐ−４、ｍ＋Ｐ〜ｍ＋２Ｐ−４の各成分に対して
Ａｌ´（１，ｊ）＝０
Ｂｌ´（１，ｊ）＝０
Ｅ（１，ｊ）＝｛Ａｌ（１，ｊ）2＋Ｂｌ（１，ｊ）2＋Ａｒ（１，ｊ）2＋Ｂｒ（１，ｊ）2｝1/2
Ａｒ´（１，ｊ）＝Ａｒ（１，ｊ）・Ｅ（１，ｊ）／｛Ａｒ（１，ｊ）2＋Ｂｒ（１，ｊ）2｝1/2
Ｂｒ´（１，ｊ）＝Ｂｒ（１，ｊ）・Ｅ（１，ｊ）／｛Ａｒ（１，ｊ）2＋Ｂｒ（１，ｊ）2｝1/2
Ａｒ´（３，ｊ）＝０
Ｂｒ´（３，ｊ）＝０
Ｅ（３，ｊ）＝｛Ａｌ（３，ｊ）2＋Ｂｌ（３，ｊ）2＋Ａｒ（３，ｊ）2＋Ｂｒ（３，ｊ）2｝1/2
Ａｌ´（３，ｊ）＝Ａｌ（３，ｊ）・Ｅ（３，ｊ）／｛Ａｌ（３，ｊ）2＋Ｂｌ（３，ｊ）2｝1/2
Ｂｌ´（３，ｊ）＝Ｂｌ（３，ｊ）・Ｅ（３，ｊ）／｛Ａｌ（３，ｊ）2＋Ｂｌ（３，ｊ）2｝1/2
上記〔数式９〕においては、ｊ＝ｍ〜ｍ＋Ｐ−４、ｍ＋Ｐ〜ｍ＋２Ｐ−４において、Ａｌ´（１，ｊ）、Ｂｌ´（１，ｊ）を共に０としている。これは、図２（ｄ）の上段に示すように、L-chにおいて、SP1U、SP1Ｄ内の各成分を０にすることを示しているが、“状態２”はSP3U、SP3Dとの差を明確にすることができれば十分であるため、必ずしも０にする必要はなく、小さな値であれば良い。
Ｅ1およびＥ4が共にレベル下限値Ｌｅｖより大きく、埋め込むべき情報が“値３”である場合、以下の〔数式１０〕に従った処理を実行することにより、低周波成分の状態を“状態３” 、すなわち、図２（ｅ）に示したような状態に変更する。
〔数式１０〕
ｊ＝ｍ〜ｍ＋Ｐ−４の各成分に対して
Ａｌ´（３，ｊ）＝０
Ｂｌ´（３，ｊ）＝０
Ｅ（３，ｊ）＝｛Ａｌ（３，ｊ）2＋Ｂｌ（３，ｊ）2＋Ａｒ（３，ｊ）2＋Ｂｒ（３，ｊ）2｝1/2
Ａｒ´（３，ｊ）＝Ａｒ（３，ｊ）・Ｅ（３，ｊ）／｛Ａｒ（３，ｊ）2＋Ｂｒ（３，ｊ）2｝1/2
Ｂｒ´（３，ｊ）＝Ｂｒ（３，ｊ）・Ｅ（３，ｊ）／｛Ａｒ（３，ｊ）2＋Ｂｒ（３，ｊ）2｝1/2
ｊ＝ｍ＋Ｐ〜ｍ＋２Ｐ−４の各成分に対して
Ａｌ´（１，ｊ）＝０
Ｂｌ´（１，ｊ）＝０
Ｅ（１，ｊ）＝｛Ａｌ（１，ｊ）2＋Ｂｌ（１，ｊ）2＋Ａｒ（１，ｊ）2＋Ｂｒ（１，ｊ）2｝1/2
Ａｒ´（１，ｊ）＝Ａｒ（１，ｊ）・Ｅ（１，ｊ）／｛Ａｒ（１，ｊ）2＋Ｂｒ（１，ｊ）2｝1/2
Ｂｒ´（１，ｊ）＝Ｂｒ（１，ｊ）・Ｅ（１，ｊ）／｛Ａｒ（１，ｊ）2＋Ｂｒ（１，ｊ）2｝1/2
ｊ＝ｍ〜ｍ＋Ｐ−４の各成分に対して
Ａｒ´（１，ｊ）＝０
Ｂｒ´（１，ｊ）＝０
Ｅ（１，ｊ）＝｛Ａｌ（１，ｊ）2＋Ｂｌ（１，ｊ）2＋Ａｒ（１，ｊ）2＋Ｂｒ（１，ｊ）2｝1/2
Ａｌ´（１，ｊ）＝Ａｌ（１，ｊ）・Ｅ（１，ｊ）／｛Ａｌ（１，ｊ）2＋Ｂｌ（１，ｊ）2｝1/2
Ｂｌ´（１，ｊ）＝Ｂｌ（１，ｊ）・Ｅ（１，ｊ）／｛Ａｌ（１，ｊ）2＋Ｂｌ（１，ｊ）2｝1/2
ｊ＝ｍ＋Ｐ〜ｍ＋２Ｐ−４の各成分に対して
Ａｒ´（３，ｊ）＝０
Ｂｒ´（３，ｊ）＝０
Ｅ（３，ｊ）＝｛Ａｌ（３，ｊ）2＋Ｂｌ（３，ｊ）2＋Ａｒ（３，ｊ）2＋Ｂｒ（３，ｊ）2｝1/2
Ａｌ´（３，ｊ）＝Ａｌ（３，ｊ）・Ｅ（３，ｊ）／｛Ａｌ（３，ｊ）2＋Ｂｌ（３，ｊ）2｝1/2
Ｂｌ´（３，ｊ）＝Ｂｌ（３，ｊ）・Ｅ（３，ｊ）／｛Ａｌ（３，ｊ）2＋Ｂｌ（３，ｊ）2｝1/2
ｊ＝ｍ〜ｍ＋Ｐ−４の各成分に対して
Ａｌ´（３，ｊ）＝０
Ｂｌ´（３，ｊ）＝０
Ｅ（３，ｊ）＝｛Ａｌ（３，ｊ）2＋Ｂｌ（３，ｊ）2＋Ａｒ（３，ｊ）2＋Ｂｒ（３，ｊ）2｝1/2
Ａｒ´（３，ｊ）＝Ａｒ（３，ｊ）・Ｅ（３，ｊ）／｛Ａｒ（３，ｊ）2＋Ｂｒ（３，ｊ）2｝1/2
Ｂｒ´（３，ｊ）＝Ｂｒ（３，ｊ）・Ｅ（３，ｊ）／｛Ａｒ（３，ｊ）2＋Ｂｒ（３，ｊ）2｝1/2
ｊ＝ｍ＋Ｐ〜ｍ＋２Ｐ−４の各成分に対して
Ａｌ´（１，ｊ）＝０
Ｂｌ´（１，ｊ）＝０
Ｅ（１，ｊ）＝｛Ａｌ（１，ｊ）2＋Ｂｌ（１，ｊ）2＋Ａｒ（１，ｊ）2＋Ｂｒ（１，ｊ）2｝1/2
Ａｒ´（１，ｊ）＝Ａｒ（１，ｊ）・Ｅ（１，ｊ）／｛Ａｒ（１，ｊ）2＋Ｂｒ（１，ｊ）2｝1/2
Ｂｒ´（１，ｊ）＝Ｂｒ（１，ｊ）・Ｅ（１，ｊ）／｛Ａｒ（１，ｊ）2＋Ｂｒ（１，ｊ）2｝1/2
ｊ＝ｍ〜ｍ＋Ｐ−４の各成分に対して
Ａｒ´（１，ｊ）＝０
Ｂｒ´（１，ｊ）＝０
Ｅ（１，ｊ）＝｛Ａｌ（１，ｊ）2＋Ｂｌ（１，ｊ）2＋Ａｒ（１，ｊ）2＋Ｂｒ（１，ｊ）2｝1/2
Ａｌ´（１，ｊ）＝Ａｌ（１，ｊ）・Ｅ（１，ｊ）／｛Ａｌ（１，ｊ）2＋Ｂｌ（１，ｊ）2｝1/2
Ｂｌ´（１，ｊ）＝Ｂｌ（１，ｊ）・Ｅ（１，ｊ）／｛Ａｌ（１，ｊ）2＋Ｂｌ（１，ｊ）2｝1/2
ｊ＝ｍ＋Ｐ〜ｍ＋２Ｐ−４の各成分に対して
Ａｒ´（３，ｊ）＝０
Ｂｒ´（３，ｊ）＝０
Ｅ（３，ｊ）＝｛Ａｌ（３，ｊ）2＋Ｂｌ（３，ｊ）2＋Ａｒ（３，ｊ）2＋Ｂｒ（３，ｊ）2｝1/2
Ａｌ´（３，ｊ）＝Ａｌ（３，ｊ）・Ｅ（３，ｊ）／｛Ａｌ（３，ｊ）2＋Ｂｌ（３，ｊ）2｝1/2
Ｂｌ´（３，ｊ）＝Ｂｌ（３，ｊ）・Ｅ（３，ｊ）／｛Ａｌ（３，ｊ）2＋Ｂｌ（３，ｊ）2｝1/2
上記〔数式１０〕においては、ｊ＝ｍ〜ｍ＋Ｐ−４においてＡｌ´（３，ｊ）、Ｂｌ´（３，ｊ）を共に０とし、ｊ＝ｍ＋Ｐ〜ｍ＋２Ｐ−４においてＡｌ´（１，ｊ）、Ｂｌ´（１，ｊ）を共に０としている。これは、図２（ｅ）の上段に示すように、L-chにおいて、SP3Ｄ、SP1U内の各成分を０にすることを示しているが、“状態３”はSP3U、SP1Dとの差を明確にすることができれば十分であるため、必ずしも０にする必要はなく、小さな値であれば良い。
Ｅ2およびＥ3が共にレベル下限値Ｌｅｖより大きく、埋め込むべき情報が“値４”である場合、以下の〔数式１１〕に従った処理を実行することにより、低周波成分の状態を“状態４” 、すなわち、図２（ｆ）に示したような状態に変更する。
〔数式１１〕
ｊ＝ｍ〜ｍ＋Ｐ−４の各成分に対して
Ａｌ´（１，ｊ）＝０
Ｂｌ´（１，ｊ）＝０
Ｅ（１，ｊ）＝｛Ａｌ（１，ｊ）2＋Ｂｌ（１，ｊ）2＋Ａｒ（１，ｊ）2＋Ｂｒ（１，ｊ）2｝1/2
Ａｒ´（１，ｊ）＝Ａｒ（１，ｊ）・Ｅ（１，ｊ）／｛Ａｒ（１，ｊ）2＋Ｂｒ（１，ｊ）2｝1/2
Ｂｒ´（１，ｊ）＝Ｂｒ（１，ｊ）・Ｅ（１，ｊ）／｛Ａｒ（１，ｊ）2＋Ｂｒ（１，ｊ）2｝1/2
ｊ＝ｍ＋Ｐ〜ｍ＋２Ｐ−４の各成分に対して
Ａｌ´（３，ｊ）＝０
Ｂｌ´（３，ｊ）＝０
Ｅ（３，ｊ）＝｛Ａｌ（３，ｊ）2＋Ｂｌ（３，ｊ）2＋Ａｒ（３，ｊ）2＋Ｂｒ（３，ｊ）2｝1/2
Ａｒ´（３，ｊ）＝Ａｒ（３，ｊ）・Ｅ（３，ｊ）／｛Ａｒ（３，ｊ）2＋Ｂｒ（３，ｊ）2｝1/2
Ｂｒ´（３，ｊ）＝Ｂｒ（３，ｊ）・Ｅ（３，ｊ）／｛Ａｒ（３，ｊ）2＋Ｂｒ（３，ｊ）2｝1/2
ｊ＝ｍ〜ｍ＋Ｐ−４の各成分に対して
Ａｒ´（３，ｊ）＝０
Ｂｒ´（３，ｊ）＝０
Ｅ（３，ｊ）＝｛Ａｌ（３，ｊ）2＋Ｂｌ（３，ｊ）2＋Ａｒ（３，ｊ）2＋Ｂｒ（３，ｊ）2｝1/2
Ａｌ´（３，ｊ）＝Ａｌ（３，ｊ）・Ｅ（３，ｊ）／｛Ａｌ（３，ｊ）2＋Ｂｌ（３，ｊ）2｝1/2
Ｂｌ´（３，ｊ）＝Ｂｌ（３，ｊ）・Ｅ（３，ｊ）／｛Ａｌ（３，ｊ）2＋Ｂｌ（３，ｊ）2｝1/2
ｊ＝ｍ＋Ｐ〜ｍ＋２Ｐ−４の各成分に対して
Ａｒ´（１，ｊ）＝０
Ｂｒ´（１，ｊ）＝０
Ｅ（１，ｊ）＝｛Ａｌ（１，ｊ）2＋Ｂｌ（１，ｊ）2＋Ａｒ（１，ｊ）2＋Ｂｒ（１，ｊ）2｝1/2
Ａｌ´（１，ｊ）＝Ａｌ（１，ｊ）・Ｅ（１，ｊ）／｛Ａｌ（１，ｊ）2＋Ｂｌ（１，ｊ）2｝1/2
Ｂｌ´（１，ｊ）＝Ｂｌ（１，ｊ）・Ｅ（１，ｊ）／｛Ａｌ（１，ｊ）2＋Ｂｌ（１，ｊ）2｝1/2
ｊ＝ｍ〜ｍ＋Ｐ−４の各成分に対して
Ａｌ´（１，ｊ）＝０
Ｂｌ´（１，ｊ）＝０
Ｅ（１，ｊ）＝｛Ａｌ（１，ｊ）2＋Ｂｌ（１，ｊ）2＋Ａｒ（１，ｊ）2＋Ｂｒ（１，ｊ）2｝1/2
Ａｒ´（１，ｊ）＝Ａｒ（１，ｊ）・Ｅ（１，ｊ）／｛Ａｒ（１，ｊ）2＋Ｂｒ（１，ｊ）2｝1/2
Ｂｒ´（１，ｊ）＝Ｂｒ（１，ｊ）・Ｅ（１，ｊ）／｛Ａｒ（１，ｊ）2＋Ｂｒ（１，ｊ）2｝1/2
ｊ＝ｍ＋Ｐ〜ｍ＋２Ｐ−４の各成分に対して
Ａｌ´（３，ｊ）＝０
Ｂｌ´（３，ｊ）＝０
Ｅ（３，ｊ）＝｛Ａｌ（３，ｊ）2＋Ｂｌ（３，ｊ）2＋Ａｒ（３，ｊ）2＋Ｂｒ（３，ｊ）2｝1/2
Ａｒ´（３，ｊ）＝Ａｒ（３，ｊ）・Ｅ（３，ｊ）／｛Ａｒ（３，ｊ）2＋Ｂｒ（３，ｊ）2｝1/2
Ｂｒ´（３，ｊ）＝Ｂｒ（３，ｊ）・Ｅ（３，ｊ）／｛Ａｒ（３，ｊ）2＋Ｂｒ（３，ｊ）2｝1/2
ｊ＝ｍ〜ｍ＋Ｐ−４の各成分に対して
Ａｒ´（３，ｊ）＝０
Ｂｒ´（３，ｊ）＝０
Ｅ（３，ｊ）＝｛Ａｌ（３，ｊ）2＋Ｂｌ（３，ｊ）2＋Ａｒ（３，ｊ）2＋Ｂｒ（３，ｊ）2｝1/2
Ａｌ´（３，ｊ）＝Ａｌ（３，ｊ）・Ｅ（３，ｊ）／｛Ａｌ（３，ｊ）2＋Ｂｌ（３，ｊ）2｝1/2
Ｂｌ´（３，ｊ）＝Ｂｌ（３，ｊ）・Ｅ（３，ｊ）／｛Ａｌ（３，ｊ）2＋Ｂｌ（３，ｊ）2｝1/2
ｊ＝ｍ＋Ｐ〜ｍ＋２Ｐ−４の各成分に対して
Ａｒ´（１，ｊ）＝０
Ｂｒ´（１，ｊ）＝０
Ｅ（１，ｊ）＝｛Ａｌ（１，ｊ）2＋Ｂｌ（１，ｊ）2＋Ａｒ（１，ｊ）2＋Ｂｒ（１，ｊ）2｝1/2
Ａｌ´（１，ｊ）＝Ａｌ（１，ｊ）・Ｅ（１，ｊ）／｛Ａｌ（１，ｊ）2＋Ｂｌ（１，ｊ）2｝1/2
Ｂｌ´（１，ｊ）＝Ｂｌ（１，ｊ）・Ｅ（１，ｊ）／｛Ａｌ（１，ｊ）2＋Ｂｌ（１，ｊ）2｝1/2
上記〔数式１１〕においては、ｊ＝ｍ〜ｍ＋Ｐ−４においてＡｌ´（１，ｊ）、Ｂｌ´（１，ｊ）を共に０とし、ｊ＝ｍ＋Ｐ〜ｍ＋２Ｐ−４においてＡｌ´（３，ｊ）、Ｂｌ´（３，ｊ）を共に０としている。これは、図２（ｆ）の上段に示すように、L-chにおいて、SP1D、SP3U内の各成分を０にすることを示しているが、“状態４”はSP3D、SP1Uとの差を明確にすることができれば十分であるため、必ずしも０にする必要はなく、小さな値であれば良い。
Ｅ1がレベル下限値Ｌｅｖより大きく、埋め込むべき情報が“値５”である場合、以下の〔数式１２〕に従った処理を実行することにより、低周波成分の状態を“状態５” 、すなわち、図２（ｇ）に示したような状態に変更する。
〔数式１２〕
ｊ＝ｍ＋Ｐ〜ｍ＋２Ｐ−４の各成分に対して
Ａｌ´（１，ｊ）＝０
Ｂｌ´（１，ｊ）＝０
Ｅ（１，ｊ）＝｛Ａｌ（１，ｊ）2＋Ｂｌ（１，ｊ）2＋Ａｒ（１，ｊ）2＋Ｂｒ（１，ｊ）2｝1/2
Ａｒ´（１，ｊ）＝Ａｒ（１，ｊ）・Ｅ（１，ｊ）／｛Ａｒ（１，ｊ）2＋Ｂｒ（１，ｊ）2｝1/2
Ｂｒ´（１，ｊ）＝Ｂｒ（１，ｊ）・Ｅ（１，ｊ）／｛Ａｒ（１，ｊ）2＋Ｂｒ（１，ｊ）2｝1/2
ｊ＝ｍ〜ｍ＋Ｐ−４、ｍ＋Ｐ〜ｍ＋２Ｐ−４の各成分に対して
Ａｌ´（３，ｊ）＝０
Ｂｌ´（３，ｊ）＝０
Ｅ（３，ｊ）＝｛Ａｌ（３，ｊ）2＋Ｂｌ（３，ｊ）2＋Ａｒ（３，ｊ）2＋Ｂｒ（３，ｊ）2｝1/2
Ａｒ´（３，ｊ）＝Ａｒ（３，ｊ）・Ｅ（３，ｊ）／｛Ａｒ（３，ｊ）2＋Ｂｒ（３，ｊ）2｝1/2
Ｂｒ´（３，ｊ）＝Ｂｒ（３，ｊ）・Ｅ（３，ｊ）／｛Ａｒ（３，ｊ）2＋Ｂｒ（３，ｊ）2｝1/2
ｊ＝ｍ〜ｍ＋Ｐ−４の各成分に対して
Ａｒ´（１，ｊ）＝０
Ｂｒ´（１，ｊ）＝０
Ｅ（１，ｊ）＝｛Ａｌ（１，ｊ）2＋Ｂｌ（１，ｊ）2＋Ａｒ（１，ｊ）2＋Ｂｒ（１，ｊ）2｝1/2
Ａｌ´（１，ｊ）＝Ａｌ（１，ｊ）・Ｅ（１，ｊ）／｛Ａｌ（１，ｊ）2＋Ｂｌ（１，ｊ）2｝1/2
Ｂｌ´（１，ｊ）＝Ｂｌ（１，ｊ）・Ｅ（１，ｊ）／｛Ａｌ（１，ｊ）2＋Ｂｌ（１，ｊ）2｝1/2
ｊ＝ｍ＋Ｐ〜ｍ＋２Ｐ−４の各成分に対して
Ａｌ´（１，ｊ）＝０
Ｂｌ´（１，ｊ）＝０
Ｅ（１，ｊ）＝｛Ａｌ（１，ｊ）2＋Ｂｌ（１，ｊ）2＋Ａｒ（１，ｊ）2＋Ｂｒ（１，ｊ）2｝1/2
Ａｒ´（１，ｊ）＝Ａｒ（１，ｊ）・Ｅ（１，ｊ）／｛Ａｒ（１，ｊ）2＋Ｂｒ（１，ｊ）2｝1/2
Ｂｒ´（１，ｊ）＝Ｂｒ（１，ｊ）・Ｅ（１，ｊ）／｛Ａｒ（１，ｊ）2＋Ｂｒ（１，ｊ）2｝1/2
ｊ＝ｍ〜ｍ＋Ｐ−４、ｍ＋Ｐ〜ｍ＋２Ｐ−４の各成分に対して
Ａｌ´（３，ｊ）＝０
Ｂｌ´（３，ｊ）＝０
Ｅ（３，ｊ）＝｛Ａｌ（３，ｊ）2＋Ｂｌ（３，ｊ）2＋Ａｒ（３，ｊ）2＋Ｂｒ（３，ｊ）2｝1/2
Ａｒ´（３，ｊ）＝Ａｒ（３，ｊ）・Ｅ（３，ｊ）／｛Ａｒ（３，ｊ）2＋Ｂｒ（３，ｊ）2｝1/2
Ｂｒ´（３，ｊ）＝Ｂｒ（３，ｊ）・Ｅ（３，ｊ）／｛Ａｒ（３，ｊ）2＋Ｂｒ（３，ｊ）2｝1/2
ｊ＝ｍ〜ｍ＋Ｐ−４の各成分に対して
Ａｒ´（１，ｊ）＝０
Ｂｒ´（１，ｊ）＝０
Ｅ（１，ｊ）＝｛Ａｌ（１，ｊ）2＋Ｂｌ（１，ｊ）2＋Ａｒ（１，ｊ）2＋Ｂｒ（１，ｊ）2｝1/2
Ａｌ´（１，ｊ）＝Ａｌ（１，ｊ）・Ｅ（１，ｊ）／｛Ａｌ（１，ｊ）2＋Ｂｌ（１，ｊ）2｝1/2
Ｂｌ´（１，ｊ）＝Ｂｌ（１，ｊ）・Ｅ（１，ｊ）／｛Ａｌ（１，ｊ）2＋Ｂｌ（１，ｊ）2｝1/2
上記〔数式１２〕においては、ｊ＝ｍ＋Ｐ〜ｍ＋２Ｐ−４においてＡｌ´（１，ｊ）、Ｂｌ´（１，ｊ）を共に０とし、ｊ＝ｍ〜ｍ＋Ｐ−４、ｍ＋Ｐ〜ｍ＋２Ｐ−４においてＡｌ´（３，ｊ）、Ｂｌ´（３，ｊ）を共に０としている。これは、図２（ｇ）の上段に示すように、L-chにおいて、SP1U、SP3U、SP3D内の各成分を０にすることを示しているが、“状態５”はSP1Dとの差を明確にすることができれば十分であるため、必ずしも０にする必要はなく、小さな値であれば良い。
Ｅ2がレベル下限値Ｌｅｖより大きく、埋め込むべき情報が“値６”である場合、以下の〔数式１３〕に従った処理を実行することにより、低周波成分の状態を“状態６” 、すなわち、図２（ｈ）に示したような状態に変更する。
〔数式１３〕
ｊ＝ｍ〜ｍ＋Ｐ−４、ｍ＋Ｐ〜ｍ＋２Ｐ−４の各成分に対して
Ａｌ´（１，ｊ）＝０
Ｂｌ´（１，ｊ）＝０
Ｅ（１，ｊ）＝｛Ａｌ（１，ｊ）2＋Ｂｌ（１，ｊ）2＋Ａｒ（１，ｊ）2＋Ｂｒ（１，ｊ）2｝1/2
Ａｒ´（１，ｊ）＝Ａｒ（１，ｊ）・Ｅ（１，ｊ）／｛Ａｒ（１，ｊ）2＋Ｂｒ（１，ｊ）2｝1/2
Ｂｒ´（１，ｊ）＝Ｂｒ（１，ｊ）・Ｅ（１，ｊ）／｛Ａｒ（１，ｊ）2＋Ｂｒ（１，ｊ）2｝1/2
ｊ＝ｍ＋Ｐ〜ｍ＋２Ｐ−４の各成分に対して
Ａｌ´（３，ｊ）＝０
Ｂｌ´（３，ｊ）＝０
Ｅ（３，ｊ）＝｛Ａｌ（３，ｊ）2＋Ｂｌ（３，ｊ）2＋Ａｒ（３，ｊ）2＋Ｂｒ（３，ｊ）2｝1/2
Ａｒ´（３，ｊ）＝Ａｒ（３，ｊ）・Ｅ（３，ｊ）／｛Ａｒ（３，ｊ）2＋Ｂｒ（３，ｊ）2｝1/2
Ｂｒ´（３，ｊ）＝Ｂｒ（３，ｊ）・Ｅ（３，ｊ）／｛Ａｒ（３，ｊ）2＋Ｂｒ（３，ｊ）2｝1/2
ｊ＝ｍ〜ｍ＋Ｐ−４の各成分に対して
Ａｒ´（３，ｊ）＝０
Ｂｒ´（３，ｊ）＝０
Ｅ（３，ｊ）＝｛Ａｌ（３，ｊ）2＋Ｂｌ（３，ｊ）2＋Ａｒ（３，ｊ）2＋Ｂｒ（３，ｊ）2｝1/2
Ａｌ´（３，ｊ）＝Ａｌ（３，ｊ）・Ｅ（３，ｊ）／｛Ａｌ（３，ｊ）2＋Ｂｌ（３，ｊ）2｝1/2
Ｂｌ´（３，ｊ）＝Ｂｌ（３，ｊ）・Ｅ（３，ｊ）／｛Ａｌ（３，ｊ）2＋Ｂｌ（３，ｊ）2｝1/2
ｊ＝ｍ〜ｍ＋Ｐ−４、ｍ＋Ｐ〜ｍ＋２Ｐ−４の各成分に対して
Ａｌ´（１，ｊ）＝０
Ｂｌ´（１，ｊ）＝０
Ｅ（１，ｊ）＝｛Ａｌ（１，ｊ）2＋Ｂｌ（１，ｊ）2＋Ａｒ（１，ｊ）2＋Ｂｒ（１，ｊ）2｝1/2
Ａｒ´（１，ｊ）＝Ａｒ（１，ｊ）・Ｅ（１，ｊ）／｛Ａｒ（１，ｊ）2＋Ｂｒ（１，ｊ）2｝1/2
Ｂｒ´（１，ｊ）＝Ｂｒ（１，ｊ）・Ｅ（１，ｊ）／｛Ａｒ（１，ｊ）2＋Ｂｒ（１，ｊ）2｝1/2
ｊ＝ｍ＋Ｐ〜ｍ＋２Ｐ−４の各成分に対して
Ａｌ´（３，ｊ）＝０
Ｂｌ´（３，ｊ）＝０
Ｅ（３，ｊ）＝｛Ａｌ（３，ｊ）2＋Ｂｌ（３，ｊ）2＋Ａｒ（３，ｊ）2＋Ｂｒ（３，ｊ）2｝1/2
Ａｒ´（３，ｊ）＝Ａｒ（３，ｊ）・Ｅ（３，ｊ）／｛Ａｒ（３，ｊ）2＋Ｂｒ（３，ｊ）2｝1/2
Ｂｒ´（３，ｊ）＝Ｂｒ（３，ｊ）・Ｅ（３，ｊ）／｛Ａｒ（３，ｊ）2＋Ｂｒ（３，ｊ）2｝1/2
ｊ＝ｍ〜ｍ＋Ｐ−４の各成分に対して
Ａｒ´（３，ｊ）＝０
Ｂｒ´（３，ｊ）＝０
Ｅ（３，ｊ）＝｛Ａｌ（３，ｊ）2＋Ｂｌ（３，ｊ）2＋Ａｒ（３，ｊ）2＋Ｂｒ（３，ｊ）2｝1/2
Ａｌ´（３，ｊ）＝Ａｌ（３，ｊ）・Ｅ（３，ｊ）／｛Ａｌ（３，ｊ）2＋Ｂｌ（３，ｊ）2｝1/2
Ｂｌ´（３，ｊ）＝Ｂｌ（３，ｊ）・Ｅ（３，ｊ）／｛Ａｌ（３，ｊ）2＋Ｂｌ（３，ｊ）2｝1/2
上記〔数式１３〕においては、ｊ＝ｍ〜ｍ＋Ｐ−４、ｍ＋Ｐ〜ｍ＋２Ｐ−４においてＡｌ´（１，ｊ）、Ｂｌ´（１，ｊ）を共に０とし、ｊ＝ｍ＋Ｐ〜ｍ＋２Ｐ−４においてＡｌ´（３，ｊ）、Ｂｌ´（３，ｊ）を共に０としている。これは、図２（ｈ）の上段に示すように、L-chにおいて、SP1U、SP3U、SP1D内の各成分を０にすることを示しているが、“状態６”はSP3Dとの差を明確にすることができれば十分であるため、必ずしも０にする必要はなく、小さな値であれば良い。
Ｅ4がレベル下限値Ｌｅｖより大きく、埋め込むべき情報が“値７”である場合、以下の〔数式１４〕に従った処理を実行することにより、低周波成分の状態を“状態７” 、すなわち、図２（ｉ）に示したような状態に変更する。
〔数式１４〕
ｊ＝ｍ〜ｍ＋Ｐ−４、ｍ＋Ｐ〜ｍ＋２Ｐ−４の各成分に対して
Ａｌ´（１，ｊ）＝０
Ｂｌ´（１，ｊ）＝０
Ｅ（１，ｊ）＝｛Ａｌ（１，ｊ）2＋Ｂｌ（１，ｊ）2＋Ａｒ（１，ｊ）2＋Ｂｒ（１，ｊ）2｝1/2
Ａｒ´（１，ｊ）＝Ａｒ（１，ｊ）・Ｅ（１，ｊ）／｛Ａｒ（１，ｊ）2＋Ｂｒ（１，ｊ）2｝1/2
Ｂｒ´（１，ｊ）＝Ｂｒ（１，ｊ）・Ｅ（１，ｊ）／｛Ａｒ（１，ｊ）2＋Ｂｒ（１，ｊ）2｝1/2
ｊ＝ｍ〜ｍ＋Ｐ−４の各成分に対して
Ａｌ´（３，ｊ）＝０
Ｂｌ´（３，ｊ）＝０
Ｅ（３，ｊ）＝｛Ａｌ（３，ｊ）2＋Ｂｌ（３，ｊ）2＋Ａｒ（３，ｊ）2＋Ｂｒ（３，ｊ）2｝1/2
Ａｒ´（３，ｊ）＝Ａｒ（３，ｊ）・Ｅ（３，ｊ）／｛Ａｒ（３，ｊ）2＋Ｂｒ（３，ｊ）2｝1/2
Ｂｒ´（３，ｊ）＝Ｂｒ（３，ｊ）・Ｅ（３，ｊ）／｛Ａｒ（３，ｊ）2＋Ｂｒ（３，ｊ）2｝1/2
ｊ＝ｍ＋Ｐ〜ｍ＋２Ｐ−４の各成分に対して
Ａｒ´（３，ｊ）＝０
Ｂｒ´（３，ｊ）＝０
Ｅ（３，ｊ）＝｛Ａｌ（３，ｊ）2＋Ｂｌ（３，ｊ）2＋Ａｒ（３，ｊ）2＋Ｂｒ（３，ｊ）2｝1/2
Ａｌ´（３，ｊ）＝Ａｌ（３，ｊ）・Ｅ（３，ｊ）／｛Ａｌ（３，ｊ）2＋Ｂｌ（３，ｊ）2｝1/2
Ｂｌ´（３，ｊ）＝Ｂｌ（３，ｊ）・Ｅ（３，ｊ）／｛Ａｌ（３，ｊ）2＋Ｂｌ（３，ｊ）2｝1/2
ｊ＝ｍ〜ｍ＋Ｐ−４、ｍ＋Ｐ〜ｍ＋２Ｐ−４の各成分に対して
Ａｌ´（１，ｊ）＝０
Ｂｌ´（１，ｊ）＝０
Ｅ（１，ｊ）＝｛Ａｌ（１，ｊ）2＋Ｂｌ（１，ｊ）2＋Ａｒ（１，ｊ）2＋Ｂｒ（１，ｊ）2｝1/2
Ａｒ´（１，ｊ）＝Ａｒ（１，ｊ）・Ｅ（１，ｊ）／｛Ａｒ（１，ｊ）2＋Ｂｒ（１，ｊ）2｝1/2
Ｂｒ´（１，ｊ）＝Ｂｒ（１，ｊ）・Ｅ（１，ｊ）／｛Ａｒ（１，ｊ）2＋Ｂｒ（１，ｊ）2｝1/2
ｊ＝ｍ〜ｍ＋Ｐ−４の各成分に対して
Ａｌ´（３，ｊ）＝０
Ｂｌ´（３，ｊ）＝０
Ｅ（３，ｊ）＝｛Ａｌ（３，ｊ）2＋Ｂｌ（３，ｊ）2＋Ａｒ（３，ｊ）2＋Ｂｒ（３，ｊ）2｝1/2
Ａｒ´（３，ｊ）＝Ａｒ（３，ｊ）・Ｅ（３，ｊ）／｛Ａｒ（３，ｊ）2＋Ｂｒ（３，ｊ）2｝1/2
Ｂｒ´（３，ｊ）＝Ｂｒ（３，ｊ）・Ｅ（３，ｊ）／｛Ａｒ（３，ｊ）2＋Ｂｒ（３，ｊ）2｝1/2
ｊ＝ｍ＋Ｐ〜ｍ＋２Ｐ−４の各成分に対して
Ａｒ´（３，ｊ）＝０
Ｂｒ´（３，ｊ）＝０
Ｅ（３，ｊ）＝｛Ａｌ（３，ｊ）2＋Ｂｌ（３，ｊ）2＋Ａｒ（３，ｊ）2＋Ｂｒ（３，ｊ）2｝1/2
Ａｌ´（３，ｊ）＝Ａｌ（３，ｊ）・Ｅ（３，ｊ）／｛Ａｌ（３，ｊ）2＋Ｂｌ（３，ｊ）2｝1/2
Ｂｌ´（３，ｊ）＝Ｂｌ（３，ｊ）・Ｅ（３，ｊ）／｛Ａｌ（３，ｊ）2＋Ｂｌ（３，ｊ）2｝1/2
上記〔数式１４〕においては、ｊ＝ｍ〜ｍ＋Ｐ−４、ｍ＋Ｐ〜ｍ＋２Ｐ−４においてＡｌ´（１，ｊ）、Ｂｌ´（１，ｊ）を共に０とし、ｊ＝ｍ〜ｍ＋Ｐ−４においてＡｌ´（３，ｊ）、Ｂｌ´（３，ｊ）を共に０としている。これは、図２（ｉ）の上段に示すように、L-chにおいて、SP1U、SP1D、SP3D内の各成分を０にすることを示しているが、“状態７”はSP3Uとの差を明確にすることができれば十分であるため、必ずしも０にする必要はなく、小さな値であれば良い。
Ｅ3がレベル下限値Ｌｅｖより大きく、埋め込むべき情報が“値８”である場合、以下の〔数式１５〕に従った処理を実行することにより、低周波成分の状態を“状態８” 、すなわち、図２（ｊ）に示したような状態に変更する。
〔数式１５〕
ｊ＝ｍ〜ｍ＋Ｐ−４の各成分に対して
Ａｌ´（１，ｊ）＝０
Ｂｌ´（１，ｊ）＝０
Ｅ（１，ｊ）＝｛Ａｌ（１，ｊ）2＋Ｂｌ（１，ｊ）2＋Ａｒ（１，ｊ）2＋Ｂｒ（１，ｊ）2｝1/2
Ａｒ´（１，ｊ）＝Ａｒ（１，ｊ）・Ｅ（１，ｊ）／｛Ａｒ（１，ｊ）2＋Ｂｒ（１，ｊ）2｝1/2
Ｂｒ´（１，ｊ）＝Ｂｒ（１，ｊ）・Ｅ（１，ｊ）／｛Ａｒ（１，ｊ）2＋Ｂｒ（１，ｊ）2｝1/2
ｊ＝ｍ〜ｍ＋Ｐ−４、ｍ＋Ｐ〜ｍ＋２Ｐ−４の各成分に対して
Ａｌ´（３，ｊ）＝０
Ｂｌ´（３，ｊ）＝０
Ｅ（３，ｊ）＝｛Ａｌ（３，ｊ）2＋Ｂｌ（３，ｊ）2＋Ａｒ（３，ｊ）2＋Ｂｒ（３，ｊ）2｝1/2
Ａｒ´（３，ｊ）＝Ａｒ（３，ｊ）・Ｅ（３，ｊ）／｛Ａｒ（３，ｊ）2＋Ｂｒ（３，ｊ）2｝1/2
Ｂｒ´（３，ｊ）＝Ｂｒ（３，ｊ）・Ｅ（３，ｊ）／｛Ａｒ（３，ｊ）2＋Ｂｒ（３，ｊ）2｝1/2
ｊ＝ｍ＋Ｐ〜ｍ＋２Ｐ−４の各成分に対して
Ａｒ´（１，ｊ）＝０
Ｂｒ´（１，ｊ）＝０
Ｅ（１，ｊ）＝｛Ａｌ（１，ｊ）2＋Ｂｌ（１，ｊ）2＋Ａｒ（１，ｊ）2＋Ｂｒ（１，ｊ）2｝1/2
Ａｌ´（１，ｊ）＝Ａｌ（１，ｊ）・Ｅ（１，ｊ）／｛Ａｌ（１，ｊ）2＋Ｂｌ（１，ｊ）2｝1/2
Ｂｌ´（１，ｊ）＝Ｂｌ（１，ｊ）・Ｅ（１，ｊ）／｛Ａｌ（１，ｊ）2＋Ｂｌ（１，ｊ）2｝1/2
ｊ＝ｍ〜ｍ＋Ｐ−４の各成分に対して
Ａｌ´（１，ｊ）＝０
Ｂｌ´（１，ｊ）＝０
Ｅ（１，ｊ）＝｛Ａｌ（１，ｊ）2＋Ｂｌ（１，ｊ）2＋Ａｒ（１，ｊ）2＋Ｂｒ（１，ｊ）2｝1/2
Ａｒ´（１，ｊ）＝Ａｒ（１，ｊ）・Ｅ（１，ｊ）／｛Ａｒ（１，ｊ）2＋Ｂｒ（１，ｊ）2｝1/2
Ｂｒ´（１，ｊ）＝Ｂｒ（１，ｊ）・Ｅ（１，ｊ）／｛Ａｒ（１，ｊ）2＋Ｂｒ（１，ｊ）2｝1/2
ｊ＝ｍ〜ｍ＋Ｐ−４、ｍ＋Ｐ〜ｍ＋２Ｐ−４の各成分に対して
Ａｌ´（３，ｊ）＝０
Ｂｌ´（３，ｊ）＝０
Ｅ（３，ｊ）＝｛Ａｌ（３，ｊ）2＋Ｂｌ（３，ｊ）2＋Ａｒ（３，ｊ）2＋Ｂｒ（３，ｊ）2｝1/2
Ａｒ´（３，ｊ）＝Ａｒ（３，ｊ）・Ｅ（３，ｊ）／｛Ａｒ（３，ｊ）2＋Ｂｒ（３，ｊ）2｝1/2
Ｂｒ´（３，ｊ）＝Ｂｒ（３，ｊ）・Ｅ（３，ｊ）／｛Ａｒ（３，ｊ）2＋Ｂｒ（３，ｊ）2｝1/2
ｊ＝ｍ＋Ｐ〜ｍ＋２Ｐ−４の各成分に対して
Ａｒ´（１，ｊ）＝０
Ｂｒ´（１，ｊ）＝０
Ｅ（１，ｊ）＝｛Ａｌ（１，ｊ）2＋Ｂｌ（１，ｊ）2＋Ａｒ（１，ｊ）2＋Ｂｒ（１，ｊ）2｝1/2
Ａｌ´（１，ｊ）＝Ａｌ（１，ｊ）・Ｅ（１，ｊ）／｛Ａｌ（１，ｊ）2＋Ｂｌ（１，ｊ）2｝1/2
Ｂｌ´（１，ｊ）＝Ｂｌ（１，ｊ）・Ｅ（１，ｊ）／｛Ａｌ（１，ｊ）2＋Ｂｌ（１，ｊ）2｝1/2
上記〔数式１５〕においては、ｊ＝ｍ〜ｍ＋Ｐ−４においてＡｌ´（１，ｊ）、Ｂｌ´（１，ｊ）を共に０とし、ｊ＝ｍ〜ｍ＋Ｐ−４、ｍ＋Ｐ〜ｍ＋２Ｐ−４においてＡｌ´（３，ｊ）、Ｂｌ´（３，ｊ）を共に０としている。これは、図２（ｊ）の上段に示すように、L-chにおいて、SP1D、SP3U、SP3D内の各成分を０にすることを示しているが、“状態８”はSP1Uとの差を明確にすることができれば十分であるため、必ずしも０にする必要はなく、小さな値であれば良い。
上記〔数式８〕〜〔数式１５〕に従った処理を実行し、低周波成分の状態を“状態１”〜“状態８”に変更した場合、いずれの場合であっても、さらに以下の〔数式１６〕に従った処理を実行し、最低周波数帯におけるギャップ領域Ｇ１Ｄ、Ｇ３Ｄ、上位低周波数帯と下位低周波数帯の間にギャップ領域Ｇ１Ｍ、Ｇ３Ｍ、上位低周波数帯の高周波側にギャップ領域Ｇ１Ｕ、Ｇ３Ｕを設ける。
〔数式１６〕
ｊ＝０〜ｍ−１、ｍ＋Ｐ−３〜ｍ＋Ｐ−１、ｍ＋２Ｐ−３〜Ｍ−１の各成分に対して
Ａｌ´（１，ｊ）＝０
Ｂｌ´（１，ｊ）＝０
Ａｌ´（３，ｊ）＝０
Ｂｌ´（３，ｊ）＝０
Ｅ（１，ｊ）＝｛Ａｌ（１，ｊ）2＋Ｂｌ（１，ｊ）2＋Ａｒ（１，ｊ）2＋Ｂｒ（１，ｊ）2｝1/2
Ａｒ´（１，ｊ）＝Ａｒ（１，ｊ）・Ｅ（１，ｊ）／｛Ａｒ（１，ｊ）2＋Ｂｒ（１，ｊ）2｝1/2
Ｂｒ´（１，ｊ）＝Ｂｒ（１，ｊ）・Ｅ（１，ｊ）／｛Ａｒ（１，ｊ）2＋Ｂｒ（１，ｊ）2｝1/2
Ｅ（３，ｊ）＝｛Ａｌ（３，ｊ）2＋Ｂｌ（３，ｊ）2＋Ａｒ（３，ｊ）2＋Ｂｒ（３，ｊ）2｝1/2
Ａｒ´（３，ｊ）＝Ａｒ（３，ｊ）・Ｅ（３，ｊ）／｛Ａｒ（３，ｊ）2＋Ｂｒ（３，ｊ）2｝1/2
Ｂｒ´（３，ｊ）＝Ｂｒ（３，ｊ）・Ｅ（３，ｊ）／｛Ａｒ（３，ｊ）2＋Ｂｒ（３，ｊ）2｝1/2
ｊ＝０〜Ｍ−１の各成分に対して
Ａｌ´（２，ｊ）＝０
Ｂｌ´（２，ｊ）＝０
Ｅ（２，ｊ）＝｛Ａｌ（２，ｊ）2＋Ｂｌ（２，ｊ）2＋Ａｒ（２，ｊ）2＋Ｂｒ（２，ｊ）2｝1/2
Ａｒ´（２，ｊ）＝Ａｒ（２，ｊ）・Ｅ（２，ｊ）／｛Ａｒ（２，ｊ）2＋Ｂｒ（２，ｊ）2｝1/2
Ｂｒ´（２，ｊ）＝Ｂｒ（２，ｊ）・Ｅ（２，ｊ）／｛Ａｒ（２，ｊ）2＋Ｂｒ（２，ｊ）2｝1/2
ｊ＝０〜ｍ−１、ｍ＋Ｐ−３〜ｍ＋Ｐ−１、ｍ＋２Ｐ−３〜Ｍ−１の各成分に対して
Ａｌ´（１，ｊ）＝０
Ｂｌ´（１，ｊ）＝０
Ａｌ´（３，ｊ）＝０
Ｂｌ´（３，ｊ）＝０
Ｅ（１，ｊ）＝｛Ａｌ（１，ｊ）2＋Ｂｌ（１，ｊ）2＋Ａｒ（１，ｊ）2＋Ｂｒ（１，ｊ）2｝1/2
Ａｒ´（１，ｊ）＝Ａｒ（１，ｊ）・Ｅ（１，ｊ）／｛Ａｒ（１，ｊ）2＋Ｂｒ（１，ｊ）2｝1/2
Ｂｒ´（１，ｊ）＝Ｂｒ（１，ｊ）・Ｅ（１，ｊ）／｛Ａｒ（１，ｊ）2＋Ｂｒ（１，ｊ）2｝1/2
Ｅ（３，ｊ）＝｛Ａｌ（３，ｊ）2＋Ｂｌ（３，ｊ）2＋Ａｒ（３，ｊ）2＋Ｂｒ（３，ｊ）2｝1/2
Ａｒ´（３，ｊ）＝Ａｒ（３，ｊ）・Ｅ（３，ｊ）／｛Ａｒ（３，ｊ）2＋Ｂｒ（３，ｊ）2｝1/2
Ｂｒ´（３，ｊ）＝Ｂｒ（３，ｊ）・Ｅ（３，ｊ）／｛Ａｒ（３，ｊ）2＋Ｂｒ（３，ｊ）2｝1/2
ｊ＝０〜Ｍ−１の各成分に対して
Ａｌ´（２，ｊ）＝０
Ｂｌ´（２，ｊ）＝０
Ｅ（２，ｊ）＝｛Ａｌ（２，ｊ）2＋Ｂｌ（２，ｊ）2＋Ａｒ（２，ｊ）2＋Ｂｒ（２，ｊ）2｝1/2
Ａｒ´（２，ｊ）＝Ａｒ（２，ｊ）・Ｅ（２，ｊ）／｛Ａｒ（２，ｊ）2＋Ｂｒ（２，ｊ）2｝1/2
Ｂｒ´（２，ｊ）＝Ｂｒ（２，ｊ）・Ｅ（２，ｊ）／｛Ａｒ（２，ｊ）2＋Ｂｒ（２，ｊ）2｝1/2
埋め込まれる側の低周波成分の強度が小さすぎると、信号を変化させても、その変化を抽出側で検出することができない。そこで、このような場合には、付加情報が埋め込まれている音響フレームとの違いを明確にするため、以下の〔数式１７〕に従った処理を実行し、その低周波成分を除去する。この状態を“状態０”と呼ぶことにする。これは、図２（ｂ）に示すような状態である。
〔数式１７〕
ｊ＝０〜Ｍ−１の各成分に対して
Ａｌ´（１，ｊ）＝０
Ｂｌ´（１，ｊ）＝０
Ａｌ´（２，ｊ）＝０
Ｂｌ´（２，ｊ）＝０
Ａｌ´（３，ｊ）＝０
Ｂｌ´（３，ｊ）＝０
Ｅ（１，ｊ）＝｛Ａｌ（１，ｊ）2＋Ｂｌ（１，ｊ）2＋Ａｒ（１，ｊ）2＋Ｂｒ（１，ｊ）2｝1/2
Ａｒ´（１，ｊ）＝Ａｒ（１，ｊ）・Ｅ（１，ｊ）／｛Ａｒ（１，ｊ）2＋Ｂｒ（１，ｊ）2｝1/2
Ｂｒ´（１，ｊ）＝Ｂｒ（１，ｊ）・Ｅ（１，ｊ）／｛Ａｒ（１，ｊ）2＋Ｂｒ（１，ｊ）2｝1/2
Ｅ（２，ｊ）＝｛Ａｌ（２，ｊ）2＋Ｂｌ（２，ｊ）2＋Ａｒ（２，ｊ）2＋Ｂｒ（２，ｊ）2｝1/2
Ａｒ´（２，ｊ）＝Ａｒ（２，ｊ）・Ｅ（２，ｊ）／｛Ａｒ（２，ｊ）2＋Ｂｒ（２，ｊ）2｝1/2
Ｂｒ´（２，ｊ）＝Ｂｒ（２，ｊ）・Ｅ（２，ｊ）／｛Ａｒ（２，ｊ）2＋Ｂｒ（２，ｊ）2｝1/2
Ｅ（３，ｊ）＝｛Ａｌ（３，ｊ）2＋Ｂｌ（３，ｊ）2＋Ａｒ（３，ｊ）2＋Ｂｒ（３，ｊ）2｝1/2
Ａｒ´（３，ｊ）＝Ａｒ（３，ｊ）・Ｅ（３，ｊ）／｛Ａｒ（３，ｊ）2＋Ｂｒ（３，ｊ）2｝1/2
Ｂｒ´（３，ｊ）＝Ｂｒ（３，ｊ）・Ｅ（３，ｊ）／｛Ａｒ（３，ｊ）2＋Ｂｒ（３，ｊ）2｝1/2
ｊ＝０〜Ｍ−１の各成分に対して
Ａｌ´（１，ｊ）＝０
Ｂｌ´（１，ｊ）＝０
Ａｌ´（２，ｊ）＝０
Ｂｌ´（２，ｊ）＝０
Ａｌ´（３，ｊ）＝０
Ｂｌ´（３，ｊ）＝０
Ｅ（１，ｊ）＝｛Ａｌ（１，ｊ）2＋Ｂｌ（１，ｊ）2＋Ａｒ（１，ｊ）2＋Ｂｒ（１，ｊ）2｝1/2
Ａｒ´（１，ｊ）＝Ａｒ（１，ｊ）・Ｅ（１，ｊ）／｛Ａｒ（１，ｊ）2＋Ｂｒ（１，ｊ）2｝1/2
Ｂｒ´（１，ｊ）＝Ｂｒ（１，ｊ）・Ｅ（１，ｊ）／｛Ａｒ（１，ｊ）2＋Ｂｒ（１，ｊ）2｝1/2
Ｅ（２，ｊ）＝｛Ａｌ（２，ｊ）2＋Ｂｌ（２，ｊ）2＋Ａｒ（２，ｊ）2＋Ｂｒ（２，ｊ）2｝1/2
Ａｒ´（２，ｊ）＝Ａｒ（２，ｊ）・Ｅ（２，ｊ）／｛Ａｒ（２，ｊ）2＋Ｂｒ（２，ｊ）2｝1/2
Ｂｒ´（２，ｊ）＝Ｂｒ（２，ｊ）・Ｅ（２，ｊ）／｛Ａｒ（２，ｊ）2＋Ｂｒ（２，ｊ）2｝1/2
Ｅ（３，ｊ）＝｛Ａｌ（３，ｊ）2＋Ｂｌ（３，ｊ）2＋Ａｒ（３，ｊ）2＋Ｂｒ（３，ｊ）2｝1/2
Ａｒ´（３，ｊ）＝Ａｒ（３，ｊ）・Ｅ（３，ｊ）／｛Ａｒ（３，ｊ）2＋Ｂｒ（３，ｊ）2｝1/2
Ｂｒ´（３，ｊ）＝Ｂｒ（３，ｊ）・Ｅ（３，ｊ）／｛Ａｒ（３，ｊ）2＋Ｂｒ（３，ｊ）2｝1/2
上記のように、Ａタイプの音響フレームに対しては、埋め込むべきビット値に応じて、低周波成分の変更を行う必要があるため、低周波成分変更手段３０は、上記〔数式８〕〜〔数式１７〕に従った処理を実行する。ところが、Ｂタイプの音響フレームは、Ａタイプの音響フレームのみの場合に生じる両端部分の不連続性を防止するために用いられるものであるので、ビット値に応じて低周波成分の変更を行う必要はない。そこで、低周波成分変更手段３０は、Ｂタイプの音響フレームについては、以下の〔数式１８〕に従った処理を実行し、常にその低周波成分を除去する。
〔数式１８〕
ｊ＝０〜Ｍ−１の各成分に対して
Ａｌ´（４，ｊ）＝０
Ｂｌ´（４，ｊ）＝０
Ｅ（４，ｊ）＝｛Ａｌ（４，ｊ）2＋Ｂｌ（４，ｊ）2＋Ａｒ（４，ｊ）2＋Ｂｒ（４，ｊ）2｝1/2
Ａｒ´（４，ｊ）＝Ａｒ（４，ｊ）・Ｅ（４，ｊ）／｛Ａｒ（４，ｊ）2＋Ｂｒ（４，ｊ）2｝1/2
Ｂｒ´（４，ｊ）＝Ｂｒ（４，ｊ）・Ｅ（４，ｊ）／｛Ａｒ（４，ｊ）2＋Ｂｒ（４，ｊ）2｝1/2
ｊ＝０〜Ｍ−１の各成分に対して
Ａｌ´（４，ｊ）＝０
Ｂｌ´（４，ｊ）＝０
Ｅ（４，ｊ）＝｛Ａｌ（４，ｊ）2＋Ｂｌ（４，ｊ）2＋Ａｒ（４，ｊ）2＋Ｂｒ（４，ｊ）2｝1/2
Ａｒ´（４，ｊ）＝Ａｒ（４，ｊ）・Ｅ（４，ｊ）／｛Ａｒ（４，ｊ）2＋Ｂｒ（４，ｊ）2｝1/2
Ｂｒ´（４，ｊ）＝Ｂｒ（４，ｊ）・Ｅ（４，ｊ）／｛Ａｒ（４，ｊ）2＋Ｂｒ（４，ｊ）2｝1/2
周波数逆変換手段４０は、上記のようにして、低周波成分の状態が変更されたフレームスペクトルを周波数逆変換して改変音響フレームを得る処理を行う。この周波数逆変換は、当然のことながら、周波数変換手段２０が実行した手法に対応していることが必要となる。本実施形態では、周波数変換手段２０において、フーリエ変換を施しているため、周波数逆変換手段４０は、フーリエ逆変換を実行することになる。
具体的には、Ａタイプの音響フレームに対しては、周波数逆変換手段４０は、上記〔数式８〕〜〔数式１７〕のいずれかにより得られたスペクトルの左チャンネルの実部Ａｌ´（１，ｊ）等、虚部Ｂｌ´（１，ｊ）等、右チャンネルの実部Ａｒ´（１，ｊ）等、虚部Ｂｒ´（１，ｊ）等を用いて、以下の〔数式１９〕に従った処理を行い、Ｘｌ´（ｉ）、Ｘｒ´（ｉ）を算出する。なお、上記〔数式８〕〜〔数式１７〕において改変されていない周波数成分については、Ａｌ´（１，ｊ）等として、元の周波数成分であるＡｌ（１，ｊ）等を用いる。
〔数式１９〕
Ｘｌ´（ｉ）＝１／Ｎ・｛ΣjＡｌ´（１，ｊ）・ｃｏｓ（２πｉｊ／Ｎ）−ΣjＢｌ´（１，ｊ）・ｓｉｎ（２πｉｊ／Ｎ）｝＋１／Ｎ・｛ΣjＡｌ´（２，ｊ）・ｃｏｓ（２πｉｊ／Ｎ）−ΣjＢｌ´（２，ｊ）・ｓｉｎ（２πｉｊ／Ｎ）｝＋１／Ｎ・｛ΣjＡｌ´（３，ｊ）・ｃｏｓ（２πｉｊ／Ｎ）−ΣjＢｌ´（３，ｊ）・ｓｉｎ（２πｉｊ／Ｎ）｝＋Ｘｌｐ（ｉ＋Ｎ／２）
Ｘｒ´（ｉ）＝１／Ｎ・｛ΣjＡｒ´（１，ｊ）・ｃｏｓ（２πｉｊ／Ｎ）−ΣjＢｒ´（１，ｊ）・ｓｉｎ（２πｉｊ／Ｎ）｝＋１／Ｎ・｛ΣjＡｒ´（２，ｊ）・ｃｏｓ（２πｉｊ／Ｎ）−ΣjＢｒ´（２，ｊ）・ｓｉｎ（２πｉｊ／Ｎ）｝＋１／Ｎ・｛ΣjＡｒ´（３，ｊ）・ｃｏｓ（２πｉｊ／Ｎ）−ΣjＢｒ´（３，ｊ）・ｓｉｎ（２πｉｊ／Ｎ）｝＋Ｘｒｐ（ｉ＋Ｎ／２）
Ｘｌ´（ｉ）＝１／Ｎ・｛ΣjＡｌ´（１，ｊ）・ｃｏｓ（２πｉｊ／Ｎ）−ΣjＢｌ´（１，ｊ）・ｓｉｎ（２πｉｊ／Ｎ）｝＋１／Ｎ・｛ΣjＡｌ´（２，ｊ）・ｃｏｓ（２πｉｊ／Ｎ）−ΣjＢｌ´（２，ｊ）・ｓｉｎ（２πｉｊ／Ｎ）｝＋１／Ｎ・｛ΣjＡｌ´（３，ｊ）・ｃｏｓ（２πｉｊ／Ｎ）−ΣjＢｌ´（３，ｊ）・ｓｉｎ（２πｉｊ／Ｎ）｝＋Ｘｌｐ（ｉ＋Ｎ／２）
Ｘｒ´（ｉ）＝１／Ｎ・｛ΣjＡｒ´（１，ｊ）・ｃｏｓ（２πｉｊ／Ｎ）−ΣjＢｒ´（１，ｊ）・ｓｉｎ（２πｉｊ／Ｎ）｝＋１／Ｎ・｛ΣjＡｒ´（２，ｊ）・ｃｏｓ（２πｉｊ／Ｎ）−ΣjＢｒ´（２，ｊ）・ｓｉｎ（２πｉｊ／Ｎ）｝＋１／Ｎ・｛ΣjＡｒ´（３，ｊ）・ｃｏｓ（２πｉｊ／Ｎ）−ΣjＢｒ´（３，ｊ）・ｓｉｎ（２πｉｊ／Ｎ）｝＋Ｘｒｐ（ｉ＋Ｎ／２）
上記〔数式１９〕においては、式が繁雑になるのを防ぐため、Σj=0,…,N-1をΣjとして示している。上記〔数式１９〕における第１式の“＋Ｘｌｐ（ｉ＋Ｎ／２）”、第２式の“＋Ｘｒｐ（ｉ＋Ｎ／２）”の項は、直前に改変された改変音響フレームのデータＸｌｐ（ｉ）、Ｘｒｐ（ｉ）が存在する場合に、時間軸上Ｎ／２サンプル分重複することを考慮して加算するためのものである。上記〔数式１９〕によりＡタイプの改変音響フレームの左チャンネルの各サンプルＸｌ´（ｉ）、右チャンネルの各サンプルＸｒ´（ｉ）、が得られることになる。
Ｂタイプの音響フレームに対しては、周波数逆変換手段４０は、上記〔数式１８〕により得られたスペクトルの左チャンネルの実部Ａｌ´（４，ｊ）、虚部Ｂｌ´（４，ｊ）、右チャンネルの実部Ａｒ´（４，ｊ）、虚部Ｂｒ´（４，ｊ）を用いて、以下の〔数式２０〕に従った処理を行い、Ｘｌ´（ｉ）、Ｘｒ´（ｉ）を算出する。なお、上記〔数式１８〕において改変されていない周波数成分については、以下の〔数式２０〕においてはＡｌ´（４，ｊ）、Ｂｌ´（４，ｊ）、Ａｒ´（４，ｊ）、Ｂｒ´（４，ｊ）として、元の値であるＡｌ（４，ｊ）、Ｂｌ（４，ｊ）、Ａｒ（４，ｊ）、Ｂｒ（４，ｊ）を用いる。
〔数式２０〕
Ｘｌ´（ｉ）＝１／Ｎ・｛ΣjＡｌ´（４，ｊ）・ｃｏｓ（２πｉｊ／Ｎ）−ΣjＢｌ´（４，ｊ）・ｓｉｎ（２πｉｊ／Ｎ）｝＋Ｘｌｐ（ｉ＋Ｎ／２）
Ｘｒ´（ｉ）＝１／Ｎ・｛ΣjＡｒ´（４，ｊ）・ｃｏｓ（２πｉｊ／Ｎ）−ΣjＢｒ´（４，ｊ）・ｓｉｎ（２πｉｊ／Ｎ）｝＋Ｘｒｐ（ｉ＋Ｎ／２）
Ｘｌ´（ｉ）＝１／Ｎ・｛ΣjＡｌ´（４，ｊ）・ｃｏｓ（２πｉｊ／Ｎ）−ΣjＢｌ´（４，ｊ）・ｓｉｎ（２πｉｊ／Ｎ）｝＋Ｘｌｐ（ｉ＋Ｎ／２）
Ｘｒ´（ｉ）＝１／Ｎ・｛ΣjＡｒ´（４，ｊ）・ｃｏｓ（２πｉｊ／Ｎ）−ΣjＢｒ´（４，ｊ）・ｓｉｎ（２πｉｊ／Ｎ）｝＋Ｘｒｐ（ｉ＋Ｎ／２）
上記〔数式２０〕によりＢタイプの改変音響フレームの左チャンネルの各サンプルＸｌ´（ｉ）、右チャンネルの各サンプルＸｒ´（ｉ）、が得られることになる。
改変音響フレーム出力手段５０は、周波数逆変換手段４０の処理により得られたＡタイプの改変音響フレーム、Ｂタイプの改変音響フレームを順次出力ファイルに出力する。
次に、図３に示した音響信号に対する情報の埋め込み装置の処理の全体的な流れを、図６のフローチャートに従って説明する。図３に示した装置を構成する各構成要素は、連携して図６に従った処理を実行する。図６は、付加情報１ワード分の処理に対応したものとなっている。１ワードとしては、任意のビット数に設定することができるが、通常１バイト（８ビット）に設定する。また、情報の埋め込みは、Ａタイプの音響フレームに対して行われるので、図６は、Ａタイプの音響フレームについての説明となっている。Ｂタイプの音響フレームについては、Ａタイプの音響フレームと並行して、音響フレーム読込手段１０により読み込まれ、窓関数Ｗ（４，ｉ）を利用して周波数変換手段２０により周波数変換された後、低周波成分変更手段３０により低周波成分が除去され、周波数逆変換手段４０により周波数逆変換された後、改変音響フレーム出力手段５０により出力される。
図６においては、まず、付加情報読込手段７０は、付加情報記憶部６２から付加情報を１ワード単位で読み込む（Ｓ１０１）。具体的には、音響信号に対する情報の埋め込み装置として用いられるコンピュータ内のレジスタに１ワード読み込むことになる。
次に、低周波成分変更手段３０が、レジスタに保持された１ワードから２ビットを読み込む処理を行う（Ｓ１０２）。そして、読み込む２ビットが末端２ビットでない場合は、読み込んだ２ビットの上位に１ビットを加える処理を行う。これは、１０進数に換算した場合に４を加えることに相当する（Ｓ１０３）。
続いて、音響フレーム読込手段１０が、音響信号記憶部６１に記憶されたステレオ音響信号の左右の各チャンネルから、それぞれ所定数のサンプルを１音響フレームとして読み込む（Ｓ１０４）。
続いて、周波数変換手段２０は、読み込んだ音響フレームに対して、周波数変換を行って、その音響フレームのスペクトルであるフレームスペクトルを得る（Ｓ１０９）。具体的には、各音響フレームについて、３つの窓関数Ｗ（１，ｉ）、Ｗ（２，ｉ）、Ｗ（３，ｉ）を用いて、上記〔数式５〕に従った処理を行う。
続いて、低周波成分変更手段３０は、上記〔数式７〕によりＥ1〜Ｅ4を算出する。そして、Ｅ1〜Ｅ4がレベル下限値Ｌｅｖ以下である場合、低周波成分を“状態０”とする処理を行う（Ｓ１１２）。具体的には、上記〔数式１７〕に従った処理を実行することになる。
上記〔数式１７〕に従った処理を実行することにより、左チャンネルのフレームスペクトルの低周波数成分は、上位低周波領域前部のスペクトルＳＰ１Ｕ、上位低周波領域後部のスペクトルＳＰ３Ｕ、下位低周波領域前部のスペクトルＳＰ１Ｄ、下位低周波領域後部のスペクトルＳＰ３Ｄ全てにおいて、“０” で同一となる。この“状態０”のパターンは、合算値Ｅ1〜Ｅ4が下限値Ｌｅｖ以下であることにより、当該フレームに値１〜値８の記録が行なわれなかったことを示す情報となる。なお、上記〔数式１７〕においては、ｊ＝０〜Ｍ−１の各成分に対してＡｌ´（ｊ）＝Ｂｌ´（ｊ）＝０としているが、抽出側において、“状態１”〜“状態８”のいずれでもないことを認識可能とすることを目的としているため、十分小さな値であれば、必ずしも０とする必要はない。また、必ずしも４つのスペクトル成分の値が同一である必要はなく、差が小さければ良い。
一方、低周波成分変更手段３０は、合算値Ｅ1〜Ｅ4がレベル下限値Ｌｅｖより大である場合、埋め込むべき情報が値１〜値８の場合に、低周波成分の状態を“状態１”〜“状態８”のいずれかの状態に変更する処理を行う（Ｓ１１１）。すなわち、Ｓ１１１においては、上記〔数式８〕〜〔数式１１〕に従った処理を実行する。
次に、周波数逆変換手段４０が、上記Ｓ１１１、Ｓ１１２の処理により各窓成分のスペクトル集合間の割合が変更されたフレームスペクトルを周波数逆変換して改変音響フレームを得る処理を行う（Ｓ１１４）。この周波数逆変換は、当然のことながら、周波数変換手段２０がＳ１０９において実行した手法に対応していることが必要となる。本実施形態では、周波数変換手段２０において、フーリエ逆変換を施しているため、周波数逆変換手段４０は、フーリエ逆変換を実行することになる。具体的には、上記〔数式８〕〜〔数式１７〕のいずれかにより得られたスペクトルの左チャンネルの実部Ａｌ´（１，ｊ）等、虚部Ｂｌ´（１，ｊ）等、右チャンネルの実部Ａｒ´（１，ｊ）等、虚部Ｂｒ´（１，ｊ）等を用いて、上記〔数式１９〕に従った処理を行い、Ｘｌ´（ｉ）、Ｘｒ´（ｉ）を算出する。
改変音響フレーム出力手段５０は、得られた改変音響フレームを順次出力ファイルに出力する。こうして１つの音響フレームに対する処理を終えたら、状態０に設定した場合には、音響フレーム読込手段１０が、次の音響フレームを読み込む（Ｓ１０４）。一方、状態１〜状態８のいずれかに設定した場合には、低周波成分変更手段３０が付加情報のビット配列中の次の２ビットを読み込む（Ｓ１０２）。以上のような処理を音響信号の両チャンネルの全サンプルに渡って実行していく。すなわち、所定数のサンプルを音響フレームとして読み込み、音響信号から読み込むべき音響フレームがなくなったら（Ｓ１０４）、処理を終了する。なお、Ｓ１０１において読み込んだ１ワードのデータの各ビットに対応する処理を終えた場合、Ｓ１０２からＳ１０１に戻り、付加情報の次のワードを読み込む処理をすることになる。付加情報の全ワードに対して処理が終了した場合は、付加情報の先頭ワードに戻って処理を行う。この結果、全ての音響フレームに対して処理を行った全ての改変音響フレームが出力ファイルに記録されて、改変音響信号として得られる。得られた改変音響信号は、記憶手段６０内の改変音響信号記憶部６３に出力され、記憶される。
以上の処理による左チャンネル信号の変化の様子について、１音響フレームに１ビットを埋め込む従来の手法と比較して説明する。図７に、従来の手法によるビット埋め込みの手順を示す。図７において、図面左右方向は、時間軸方向である。また、図中多数存在する矩形は、改変音響フレームの第１成分、第２成分の存在（削除されていないこと）を示している。図７の例では、固定長８ビットの連続する４ワードの付加情報を埋め込む場合を示している。図７に示すように、最初に区切りを示す情報が埋め込まれた後は、ビット値を示す情報が埋め込まれる。低周波成分が小さいために埋め込みができないフレームが存在しない場合には、連続して１ワード分の８ビットが埋め込まれ、その後、再び、区切りを示す情報が埋め込まれ、続いて、ビット値を示す情報が埋め込まれる。低周波成分が小さいために埋め込みができないフレームが存在した場合には、図７に示す３ワード目のように、２ビット目の後、区切りを示す情報が埋め込まれ、その後、３ビット目から埋め込まれる。
これに対して、本発明では、１音響フレームに１ビットを埋め込むため、より多くの情報を埋め込むことができる。図８に、本発明によるビット埋め込みの手順を示す。図８において、第１成分と第２成分を結ぶ破線は、音響フレーム長を示しており、矩形は、改変音響フレームの第１成分、第２成分の存在を示している。図８の例でも、固定長８ビットの連続する６ワード強の付加情報を埋め込む場合を示している。図８に示すように、２ビットの値を示す情報が順次埋め込まれる。すなわち“状態５”〜“状態８”のいずれかの状態に変更される。４つ目の音響フレームは、１ワードの最後の２ビットを埋め込むものであるため、“状態１”〜“状態４”のいずれかの状態に変更される。これにより１ワードの終了を示すことができる。低周波成分が小さいために埋め込みができないフレームが存在しない場合には、連続して１ワード分の８ビットが埋め込まれ、その後、再び、区切りを示す情報が埋め込まれ、続いて、２ビットの値を示す情報が埋め込まれる。低周波成分が小さいために埋め込みができないフレームが存在した場合には、図８に示す５ワード目のように、２ビット目の後、“状態０”に変更される。なお、本実施形態で、低周波成分が小さいために埋め込みができないフレームが存在した場合とは、上記〔数式７〕により算出された合算値Ｅ1〜Ｅ4のうち埋め込むべき値に対応するものが、レベル下限値Ｌｅｖ以下となる音響フレームが存在する場合を示している。
図７、図８は、従来手法と、本発明における処理の一例であるが、図７に示すように、従来の手法では、３６個の音響フレームを利用して３０ビットしか埋め込むことができないのに対し、本発明の手法では、３１個の音響フレームを利用して５０ビット埋め込むことができ、埋め込み効率が高まる。ワード間の区切りを１つの音響フレームに埋め込むようにした場合は、図８の例では、さらに４つの音響フレームを必要とすることになり、埋め込み効率は低くなる。したがって、ワード間の区切りを４つ目の音響フレームにおいて表現した本発明の手法では、埋め込み効率がその分だけさらに高くなることを示している。
なお、本実施形態では、付加情報の１ワードを１バイトとした場合について説明したが、抽出側と取り決めがある限り、付加情報の１ワードを任意のビット数単位で記録することが可能である。
上記のようにして得られた改変音響信号の左チャンネルのうち、付加情報が埋め込まれている部分については、低周波成分は、上記状態０〜状態８の９通りの分布しかないことになる。しかし、高周波成分については、元の音響信号のままであるので、制作者の設定に基づいた種々な分布になる。また、上記の例で示したように、ステレオ音響信号を利用した場合には、左チャンネルにおいて変化させられた低周波成分は、上記〔数式１１〕〜〔数式１７〕の処理からも明らかなように、必ず右チャンネルの低周波成分に付加されている。したがって、右チャンネルが左チャンネルにおいて削除された成分を補っているため、両チャンネル全体として見ると、信号の劣化がない。人間の聴覚は、高周波成分については、方向性を感知し易いが、低周波成分については、方向性を感知しにくくなっている。したがって、低周波成分が一方に偏っていても、聴いている人にとっては、通常の音響信号と変わりなく聴こえることになる。
（４．音響信号からの情報の抽出装置）
次に、本発明に係る音響信号からの情報の抽出装置について説明する。図９は、本発明に係る音響信号からの情報の抽出装置の一実施形態を示す構成図である。図９において、１００は音響信号入力手段、１１０は基準フレーム獲得手段、１２０は位相変更フレーム設定手段、１３０は周波数変換手段、１４０は符号判定パラメータ算出手段、１５０は符号出力手段、１６０は付加情報抽出手段、１７０は音響フレーム保持手段である。
次に、本発明に係る音響信号からの情報の抽出装置について説明する。図９は、本発明に係る音響信号からの情報の抽出装置の一実施形態を示す構成図である。図９において、１００は音響信号入力手段、１１０は基準フレーム獲得手段、１２０は位相変更フレーム設定手段、１３０は周波数変換手段、１４０は符号判定パラメータ算出手段、１５０は符号出力手段、１６０は付加情報抽出手段、１７０は音響フレーム保持手段である。
音響信号入力手段１００は、流れている音声をデジタル音響信号として取得し、入力する機能を有している。現実には、マイクロフォンおよびＡ／Ｄ変換器により実現される。マイクロフォンとしては、低周波成分が検出可能なものであれば、モノラル無指向性のものであっても、ステレオ指向性のものであっても使用可能である。ステレオ指向性のものであっても一方のチャンネルだけ利用すれば良い。また、図３に示した装置で情報の埋め込みを行った場合には、特別精度の高いものでなく、一般的な精度のマイクロフォンを用いても情報の抽出が可能となる。基準フレーム獲得手段１１０は、入力されたデジタルのモノラル音響信号（あるいはステレオ音響信号の１チャンネル）から所定数のサンプルで構成される音響フレームを基準フレームとして読み込む機能を有している。位相変更フレーム設定手段１２０は、基準フレームと所定サンプルずつ移動させることにより位相を変更した音響フレームを位相変更フレームとして設定する機能を有している。周波数変換手段１３０は、図３に示した周波数変換手段２０と同様の機能を有している。
符号判定パラメータ算出手段１４０は、生成されたフレームスペクトルから所定の周波数以下に相当する各低周波強度データを抽出し、各成分に対応した各低周波強度データの合算値ＥC1〜ＥC4を以下の〔数式２１〕に基づいて算出し、この合算値ＥC1〜ＥC4を符号判定パラメータとし、この符号判定パラメータＥC1〜ＥC4の比率等に基づいて、所定の状態であると判断する機能を有している。以下の〔数式２１〕は上記〔数式７〕において右チャンネル成分を削除したもので、抽出時には右チャンネル成分を参照しないためである。
〔数式２１〕
ＥC1＝Σj=m,…,m+P-4｛Ａｌ（１，ｊ）2＋Ｂｌ（１，ｊ）2｝
ＥC2＝Σj=m,…,m+P-4｛Ａｌ（３，ｊ）2＋Ｂｌ（３，ｊ）2｝
ＥC3＝Σj=m+P,…,m+2P-4｛Ａｌ（１，ｊ）2＋Ｂｌ（１，ｊ）2｝
ＥC4＝Σj=m+P,…,m+2P-4｛Ａｌ（３，ｊ）2＋Ｂｌ（３，ｊ）2｝
ＥC1＝Σj=m,…,m+P-4｛Ａｌ（１，ｊ）2＋Ｂｌ（１，ｊ）2｝
ＥC2＝Σj=m,…,m+P-4｛Ａｌ（３，ｊ）2＋Ｂｌ（３，ｊ）2｝
ＥC3＝Σj=m+P,…,m+2P-4｛Ａｌ（１，ｊ）2＋Ｂｌ（１，ｊ）2｝
ＥC4＝Σj=m+P,…,m+2P-4｛Ａｌ（３，ｊ）2＋Ｂｌ（３，ｊ）2｝
符号出力手段１５０は、１つの基準フレームに対応する音響フレーム（基準フレームおよび位相変更フレーム）の中から最適な位相であると判断されるものを判断し、その音響フレームの状態に対応する符号を出力する機能を有している。
符号判定パラメータ算出手段１４０と符号出力手段１５０により符号化手段が構成される。付加情報抽出手段１６０は、符号出力手段１５０により出力された９値の配列を、所定の規則により変換して意味のある付加情報として抽出する機能を有している。音響フレーム保持手段１７０は、連続する２個の基準フレームを保持可能なバッファメモリである。図９に示した各構成手段は、現実には情報処理機能を有する小型のコンピュータおよびその周辺機器等のハードウェアに専用のプログラムを搭載することにより実現される。特に、本発明の目的をより簡易に達成するためには、携帯型端末装置をハードウェアとして用いることが望ましい。
符号判定パラメータ算出手段１４０と符号出力手段１５０により符号化手段が構成される。付加情報抽出手段１６０は、符号出力手段１５０により出力された９値の配列を、所定の規則により変換して意味のある付加情報として抽出する機能を有している。音響フレーム保持手段１７０は、連続する２個の基準フレームを保持可能なバッファメモリである。図９に示した各構成手段は、現実には情報処理機能を有する小型のコンピュータおよびその周辺機器等のハードウェアに専用のプログラムを搭載することにより実現される。特に、本発明の目的をより簡易に達成するためには、携帯型端末装置をハードウェアとして用いることが望ましい。
（５．抽出装置の処理動作）
次に、図９に示した音響信号からの情報の抽出装置の処理動作について図１０のフローチャートに従って説明する。まず、本装置では、平均符号レベルＨＬ１〜ＨＬ８、位相判定テーブルＳ（ｐ）、非符号カウンタＮｎを初期化する（Ｓ２００）。これらについて説明する。平均符号レベルＨＬ１〜ＨＬ８は、値１〜値８が埋め込まれていたと判断される音響フレーム（以下、有効フレームと呼ぶことにする）についての、上記〔数式２１〕で算出される低周波成分の合算値ＥC1〜ＥC4を利用して得られるＥC1＋ＥC3、ＥC2＋ＥC4、ＥC1＋ＥC4、ＥC2＋ＥC3、２・ＥC1、２・ＥC2、２・ＥC4、２・ＥC3それぞれの平均値、すなわち、過去の有効フレームにおけるこれらの値の平均値で与えられるものであり、初期値は、０．１に設定されている。位相判定テーブルＳ（ｐ）は、位相を判定するためのテーブルであり、ｐは０〜５の整数値をとる。初期値はＳ（ｐ）＝０に設定されている。非符号カウンタＮｎは、信号レベルが低く、非符号（埋め込みエラーを示す情報）であると判断されるフレーム数のカウンタであり、初期状態では、Ｎｎ＝０に設定される。
次に、図９に示した音響信号からの情報の抽出装置の処理動作について図１０のフローチャートに従って説明する。まず、本装置では、平均符号レベルＨＬ１〜ＨＬ８、位相判定テーブルＳ（ｐ）、非符号カウンタＮｎを初期化する（Ｓ２００）。これらについて説明する。平均符号レベルＨＬ１〜ＨＬ８は、値１〜値８が埋め込まれていたと判断される音響フレーム（以下、有効フレームと呼ぶことにする）についての、上記〔数式２１〕で算出される低周波成分の合算値ＥC1〜ＥC4を利用して得られるＥC1＋ＥC3、ＥC2＋ＥC4、ＥC1＋ＥC4、ＥC2＋ＥC3、２・ＥC1、２・ＥC2、２・ＥC4、２・ＥC3それぞれの平均値、すなわち、過去の有効フレームにおけるこれらの値の平均値で与えられるものであり、初期値は、０．１に設定されている。位相判定テーブルＳ（ｐ）は、位相を判定するためのテーブルであり、ｐは０〜５の整数値をとる。初期値はＳ（ｐ）＝０に設定されている。非符号カウンタＮｎは、信号レベルが低く、非符号（埋め込みエラーを示す情報）であると判断されるフレーム数のカウンタであり、初期状態では、Ｎｎ＝０に設定される。
このように、初期値が設定されている状態で、利用者が流れている音楽について、その楽曲名等の属性情報を知りたいと思った場合、まず、抽出装置に対して、抽出装置としての起動の指示を行う。これは、例えば、抽出装置を携帯電話機等の携帯端末で実現している場合は、所定のボタンを操作することにより実行できる。抽出装置は、指示が入力されると、音響信号入力手段１００が、流れている音楽を録音し、デジタル化してデジタル音響信号として入力する。具体的には、無指向性マイクロフォン（または指向性マイクロフォンの一方のチャンネル）から入力される音声を、Ａ／Ｄ変換器によりデジタル化する処理を行うことになる。
続いて、基準フレーム獲得手段１１０が、音響信号入力手段１００から入力された音響信号から、所定数のサンプルで構成される音響フレームを基準フレームとして抽出する（Ｓ２０１）。具体的には、基準フレームを抽出して音響フレーム保持手段１７０に読み込むことになる。基準フレーム獲得手段１１０が基準フレームとして読み込む１音響フレームのサンプル数は、図３に示した音響フレーム読込手段１０で設定されたものと同一にする必要がある。したがって、本実施形態の場合、基準フレーム獲得手段１１０は、４０９６サンプルずつ、順次基準フレームとして読み込んでいくことになる。音響フレーム保持手段１７０には、上述のように２個の基準フレームが格納可能となっており、新しい基準フレームが読み込まれると、古い基準フレームを破棄するようになっている。したがって、音響フレーム保持手段１７０には、常に基準フレーム２個分（連続する８１９２サンプル）が格納されていることになる。
埋め込み装置で処理する音響フレームは、先頭から途切れることなく隣接して設定される基準フレームと、この基準フレームと位相を変更した位相変更フレームとに分けることができる。基準フレームについては、最初の基準フレームとしてサンプル番号１からサンプル番号４０９６までを設定したら、次の基準フレームは、サンプル番号４０９７からサンプル番号８１９２、さらに次の基準フレームは、サンプル番号８１９３からサンプル番号１２２８８、というように途切れることなく設定される。そして、各基準フレームについて、１／６フレーム（約６８３サンプル）ずつ移動した５個の位相変更フレームを設定する。例えば、最初の基準フレームについては、サンプル番号６８３、１３６６、２０４９、２７３２、３４１３から始まる４０９６のサンプルで構成される５個の位相変更フレームが設定されることになる。続いて、周波数変換手段１３０、符号判定パラメータ算出手段１４０が、読み込んだ各音響フレームから、埋め込まれている情報を判定し、対応する符号を出力する（Ｓ２０２）。出力される情報の形式は、付加情報が埋め込まれる場合に対応する値１〜値８の８値、および値０の計９値の形式となる。
ここで、ステップＳ２０２の符号判定処理の詳細を図１１のフローチャートに従って説明する。まず、周波数変換手段１３０が、読み込んだ各音響フレームに対して、周波数変換を行ってフレームスペクトルを得る（Ｓ４０１）。この処理は、図３に示した周波数変換手段２０における処理と同様である。ただし、抽出に用いるのは、左チャンネルだけであるので、上記〔数式５〕に従った処理を行い、左チャンネルに対応する変換データの実部Ａｌ（１，ｊ）等、虚部Ｂｌ（１，ｊ）等を得る。
上記周波数変換手段１３０における処理により、周波数に対応した成分であるスペクトルで表現されたフレームスペクトルが得られる。続いて、符号判定パラメータ算出手段１４０は、平均符号レベルＨＬ１〜ＨＬ８の算出を行う（Ｓ４０２）。具体的には、過去“状態１”であると判断された音響フレームについての合算値ＥC1と合算値ＥC3の和の積算値であるｖ１を、過去“状態１”であると判断された音響フレームの数であるｎ１で除算することによりＨＬ１を算出し、過去“状態２”であると判断された音響フレームについての合算値ＥC2と合算値ＥC4の和の積算値であるｖ２を、過去“状態２”であると判断された音響フレームの数であるｎ２で除算することによりＨＬ２を算出し、過去“状態３”であると判断された音響フレームについての合算値ＥC1と合算値ＥC4の和の積算値であるｖ３を、過去“状態３”であると判断された音響フレームの数であるｎ３で除算することによりＨＬ３を算出し、過去“状態４”であると判断された音響フレームについての合算値ＥC2と合算値ＥC3の和の積算値であるｖ４を、過去“状態４”であると判断された音響フレームの数であるｎ４で除算することによりＨＬ４を算出する。
さらに、過去“状態５”であると判断された音響フレームについての合算値ＥC1の２倍の積算値であるｖ５を、過去“状態５”であると判断された音響フレームの数であるｎ５で除算することによりＨＬ５を算出し、過去“状態６”であると判断された音響フレームについての合算値ＥC2の２倍の積算値であるｖ６を、過去“状態６”であると判断された音響フレームの数であるｎ６で除算することによりＨＬ６を算出し、過去“状態７”であると判断された音響フレームについての合算値ＥC3の２倍の積算値であるｖ７を、過去“状態７”であると判断された音響フレームの数であるｎ７で除算することによりＨＬ７を算出し、過去“状態８”であると判断された音響フレームについての合算値ＥC4の２倍の積算値であるｖ８を、過去“状態８”であると判断された音響フレームの数であるｎ８で除算することによりＨＬ８を算出する。したがって、平均符号レベルＨＬ１〜ＨＬ４は、過去対応する状態であると判断された音響フレームの低周波強度データの合算値の平均値となる。平均符号レベルＨＬ５〜ＨＬ８は、過去対応する状態であると判断された音響フレームの低周波強度データの合算値の２倍の平均値となる。
続いて、符号判定パラメータ算出手段１４０は、候補符号テーブルの初期化を行う（Ｓ４０３）。候補符号テーブルは、１つの基準フレームおよび５個の位相変更フレームを特定する０〜５の位相番号および、この６個の音響フレームの状態から得られる９値の符号を記録するものである。
続いて、符号判定パラメータ算出手段１４０は、合算値ＥC1〜合算値ＥC4を用いて、低周波成分の状態がどのような状態であるか、すなわち、２ビットの値としてどのような値が埋め込まれていたかを判断する処理を行う（Ｓ４０４）。具体的には、まず、以下の４群の判定処理を実行し、各群における全ての条件を満たす場合に、対応する状態であったと判断し、対応する値を出力する。
ＥC1＞０．０００５・ＨＬ１、かつＥC3＞０．０００５・ＨＬ３、かつＥC1＞２・ＥC2、かつＥC3＞２・ＥC4の場合、“状態１”であると判断し、値１を出力する。
ＥC2＞０．０００５・ＨＬ２、かつＥC4＞０．０００５・ＨＬ４、かつＥC2＞２・ＥC1、かつＥC4＞２・ＥC3の場合、“状態２”であると判断し、値２を出力する。
ＥC1＞０．０００５・ＨＬ１、かつＥC4＞０．０００５・ＨＬ４、かつＥC1＞２・ＥC2、かつＥC4＞２・ＥC3の場合、“状態３”であると判断し、値３を出力する。
ＥC2＞０．０００５・ＨＬ２、かつＥC3＞０．０００５・ＨＬ３、かつＥC2＞２・ＥC1、かつＥC3＞２・ＥC4の場合、“状態４”であると判断し、値４を出力する。
上記各群のいずれも満たさない場合は、符号判定パラメータ算出手段１４０は、さらに以下の各群における判定処理を実行して状態を判断し、対応する値を出力する。
ＥC1＞０．０００５・ＨＬ１、かつＥC1＞２・ＥC2、かつＥC1＞ＥC3、かつＥC1＞ＥC4の場合、“状態５”であると判断し、値５を出力する。
ＥC2＞０．０００５・ＨＬ２、かつＥC2＞２・ＥC1、かつＥC2＞ＥC3、かつＥC2＞ＥC4の場合、“状態６”であると判断し、値６を出力する。
ＥC4＞０．０００５・ＨＬ４、かつＥC4＞２・ＥC3、かつＥC4＞ＥC1、かつＥC4＞ＥC2の場合、“状態７”であると判断し、値７を出力する。
ＥC3＞０．０００５・ＨＬ３、かつＥC3＞２・ＥC4、かつＥC3＞ＥC1、かつＥC3＞ＥC2の場合、“状態８”であると判断し、値８を出力する。
符号判定パラメータ算出手段１４０は、各音響フレーム単位で、上記判定結果に応じて値１〜値８のいずれかを出力する（Ｓ４０７）。上記いずれの群の条件も満たさない場合は、“状態０”であると判断し、値０を出力する（Ｓ４０８）。この“状態０”は、合算値Ｅ1〜Ｅ4が下限値Ｌｅｖ以下であることにより、当該フレームに値１〜値８の記録が行なわれなかったことを示す埋め込みエラーを示すものとなる。
また、上記判定の結果、値１〜値８のいずれかを出力した場合には、さらに、以下の〔数式２２〕に従って位相判定テーブルＳ（ｐ）の更新を行う（Ｓ４０９）。
〔数式２２〕
“状態１”であると判断し、値１を出力した場合、Ｓ（ｐ）←Ｓ（ｐ）＋ＥC1＋ＥC3
“状態２”であると判断し、値２を出力した場合、Ｓ（ｐ）←Ｓ（ｐ）＋ＥC2＋ＥC4
“状態３”であると判断し、値３を出力した場合、Ｓ（ｐ）←Ｓ（ｐ）＋ＥC1＋ＥC4
“状態４”であると判断し、値４を出力した場合、Ｓ（ｐ）←Ｓ（ｐ）＋ＥC2＋ＥC3
“状態５”であると判断し、値５を出力した場合、Ｓ（ｐ）←Ｓ（ｐ）＋２・ＥC1
“状態６”であると判断し、値６を出力した場合、Ｓ（ｐ）←Ｓ（ｐ）＋２・ＥC2
“状態７”であると判断し、値７を出力した場合、Ｓ（ｐ）←Ｓ（ｐ）＋２・ＥC4
“状態８”であると判断し、値８を出力した場合、Ｓ（ｐ）←Ｓ（ｐ）＋２・ＥC3
“状態１”であると判断し、値１を出力した場合、Ｓ（ｐ）←Ｓ（ｐ）＋ＥC1＋ＥC3
“状態２”であると判断し、値２を出力した場合、Ｓ（ｐ）←Ｓ（ｐ）＋ＥC2＋ＥC4
“状態３”であると判断し、値３を出力した場合、Ｓ（ｐ）←Ｓ（ｐ）＋ＥC1＋ＥC4
“状態４”であると判断し、値４を出力した場合、Ｓ（ｐ）←Ｓ（ｐ）＋ＥC2＋ＥC3
“状態５”であると判断し、値５を出力した場合、Ｓ（ｐ）←Ｓ（ｐ）＋２・ＥC1
“状態６”であると判断し、値６を出力した場合、Ｓ（ｐ）←Ｓ（ｐ）＋２・ＥC2
“状態７”であると判断し、値７を出力した場合、Ｓ（ｐ）←Ｓ（ｐ）＋２・ＥC4
“状態８”であると判断し、値８を出力した場合、Ｓ（ｐ）←Ｓ（ｐ）＋２・ＥC3
続いて、符号判定パラメータ算出手段１４０は、候補符号テーブルに、最適位相となる候補を保存する（Ｓ４１０）。具体的には、位相判定テーブルに記録されているＳ（ｐ）の値が最大となる位相番号ｐの値、前記Ｓ４０７、Ｓ４０８により判定された９値のいずれかの符号、その音響フレームについての上記〔数式２１〕に従った処理を実行することにより算出した、低周波数成分に対応するＥC1＋ＥC3、ＥC2＋ＥC4、ＥC1＋ＥC4、ＥC2＋ＥC3、２・ＥC1、２・ＥC2、２・ＥC4、２・ＥC3のいずれかを最適位相の候補として候補符号テーブルに保存する。
続いて、全ての位相番号ｐに対応する処理を終えたかどうかを判定する（Ｓ４１１）。これは、ある基準フレームに対して全ての位相変更フレームの処理を行ったかどうかを判定している。本実施形態では、ｐが０〜５までの値をとるので、６回分処理していない場合は、処理していた音響フレームから所定サンプル数ずらして、位相の異なる音響フレームを設定し、Ｓ４０４に戻って処理を繰り返す。なお、ｐ＝０の場合が基準フレームであり、ｐ＝１〜５の場合が位相変更フレームである。全ての位相番号ｐに対応する処理を終えた場合は、候補保存テーブルに記録されている位相番号ｐに対応する位相が最適位相であると判定し、候補保存テーブルに記録されている符号を出力する（Ｓ４１２）。
再び図１０のフローチャートに戻って説明する。Ｓ２０２による処理の結果、値１〜値８に相当する符号が出力された場合には、平均符号レベルのパラメータの更新処理を行う（Ｓ２０３）。具体的には、平均符号レベルＨＬ１〜ＨＬ８算出の際の分子となる積算値ｖ１〜ｖ８に、それぞれＥC1＋ＥC3、ＥC2＋ＥC4、ＥC1＋ＥC4、ＥC2＋ＥC3、２・ＥC1、２・ＥC2、２・ＥC4、２・ＥC3のいずれかを加算して積算値ｖ１〜ｖ８を更新し、分母となるフレーム数ｎ１〜ｎ８にそれぞれ１を加算してフレーム数ｎ１〜ｎ８を更新することにより行う。ｎ１〜ｎ８はあらかじめ設定したｎmax（例えば１００）を超えないように設定する。
ここからは、出力された値が値１〜４であるか、値５〜８であるかにより異なる。値５〜８である場合、非符号カウンタの初期化処理を行う（Ｓ２０４）。具体的には、上記Ｓ２００の初期化処理における処理と同様、非符号カウンタＮｎ＝０と設定する。続いて、出力された値から４減じた値、すなわち値１〜４に対応する２ビットの値をバッファに保存する（Ｓ２０５）。続いて、ビットカウンタを“２”だけカウントアップする（Ｓ２０６）。そして、Ｓ２０１に戻って、次の基準フレームを抽出する処理を行う。
Ｓ２０２において出力された値が値１〜４である場合、出力された２ビットの値をバッファに保存する（Ｓ２０７）。値１〜４は、１ワードの終了を示しているので、バッファに記録された１ワード分のデータを、付加情報抽出手段１６０が出力する（Ｓ２０８）。そして、ビットカウンタを０に初期化する（Ｓ２０９）。そして、Ｓ２０１に戻って、次の基準フレームを抽出する処理を行う。
Ｓ２０２による処理の結果、値０が出力された場合には、非符号カウンタのカウントアップ処理を行う（Ｓ２１０）。具体的には、非符号カウンタＮｎの値に１を加算する。そして、非符号カウンタＮｎの値がｎmax以上である場合は、Ｓ２００に戻って初期化処理を行う。非符号カウンタＮｎの値がｎmax未満である場合は、Ｓ２０１に戻って符号抽出処理を継続する。図１０に示す処理を各基準フレームに対して実行することにより、付加情報が抽出されることになる。Ｓ２０１において全ての基準フレームが抽出されたと判断された場合には、処理を終了する。
上記Ｓ２０８の処理において、付加情報抽出手段１６０は、まず、符号出力手段１５０により出力された値を、所定の規則により変換して意味のある付加情報として抽出する。所定の規則としては、情報を埋め込む者が意図した情報が受け取った者に認識可能な状態とできるものであれば、さまざまな規則が適用できるが、本実施形態では、文字情報として認識するための規則としている。すなわち、付加情報抽出手段１６０は、符号判定パラメータ算出手段１４０が判定し、符号出力手段１５０から出力される符号から得られる２ビット値の配列を１バイト（８ビット）単位で認識し、これを設定されたコード体系に従って文字情報を認識する。このようにして得られた文字情報は、表示装置（図示省略）の画面に表示出力される。
従って、埋め込み装置により音響信号に、その楽曲の曲名やアーチスト等の属性情報を文字情報として埋め込んでおけば、利用者は、その音楽が流れているのを聞いて、その曲名やアーチストを知りたいと思ったときに、抽出装置として機能する自身の携帯端末に所定の操作を行えば、自身の携帯端末の画面に曲名やアーチスト等の属性情報が文字情報として表示されることになる。
（６．位相補正処理について）
上記のように、抽出時には、埋め込み時に埋め込んだ音響フレームに対応して、音響信号を読み込むことができるとは限らない。そこで、音響フレームの位相をずらして複数通り（本実施形態では６通り）で読み込み、その中で最適な位相を決定し、その位相で特定される音響フレームに対応する符号を出力することにしている。例えば、６通りで読み込む場合、先頭の音響フレームは、本来サンプル番号１〜４０９６のサンプルであるが、サンプル番号１、６８３、１３６６、２０４９、２７３２、３４１３から始まる４０９６のサンプルで構成される６個の各音響フレームに対して処理を行い、最適な音響フレームに対応する符号を出力することになる。この位相補正処理は、Ｓ４０３、Ｓ４０９、Ｓ４１０、Ｓ４１１、Ｓ４１２における処理を中心として行われることになる。
上記のように、抽出時には、埋め込み時に埋め込んだ音響フレームに対応して、音響信号を読み込むことができるとは限らない。そこで、音響フレームの位相をずらして複数通り（本実施形態では６通り）で読み込み、その中で最適な位相を決定し、その位相で特定される音響フレームに対応する符号を出力することにしている。例えば、６通りで読み込む場合、先頭の音響フレームは、本来サンプル番号１〜４０９６のサンプルであるが、サンプル番号１、６８３、１３６６、２０４９、２７３２、３４１３から始まる４０９６のサンプルで構成される６個の各音響フレームに対して処理を行い、最適な音響フレームに対応する符号を出力することになる。この位相補正処理は、Ｓ４０３、Ｓ４０９、Ｓ４１０、Ｓ４１１、Ｓ４１２における処理を中心として行われることになる。
（７．下限閾値補正処理について）
信号レベルが小さい場合には、窓成分の強度の大小が判定できず、抽出側で誤判断することが多くなる。そこで、合算値ＥC1〜ＥC4が所定の閾値以下のフレームについては、無効なフレームであると判断するようにしているが、この際の閾値を過去の有効フレームについての低周波強度の積算値を利用して補正する処理を行っている。このように閾値を変動させることにより、信号レベルが変動しても無効なフレームであるか、有効なフレームであるかを正確に判断することが可能となる。この下限閾値補正処理は、Ｓ４０２、Ｓ２０３における処理を中心として行われることになる。
信号レベルが小さい場合には、窓成分の強度の大小が判定できず、抽出側で誤判断することが多くなる。そこで、合算値ＥC1〜ＥC4が所定の閾値以下のフレームについては、無効なフレームであると判断するようにしているが、この際の閾値を過去の有効フレームについての低周波強度の積算値を利用して補正する処理を行っている。このように閾値を変動させることにより、信号レベルが変動しても無効なフレームであるか、有効なフレームであるかを正確に判断することが可能となる。この下限閾値補正処理は、Ｓ４０２、Ｓ２０３における処理を中心として行われることになる。
（８．モノラル音響信号の場合）
上記実施形態においては、埋め込み装置、抽出装置のいずれにおいても、左右のチャンネルを有するステレオ音響信号の左チャンネル信号に付加情報を埋め込む場合を例にとって説明したが、逆に右チャンネル信号に付加情報を埋め込むようにしても良い。本発明は、左右の特性には無関係だからである。また、１つのチャンネルしかないモノラル音響信号に対して処理を行う場合は、上記実施形態において、左チャンネル信号に対して行った処理を行うことになる。本発明は、１つのチャンネル信号に対して付加情報を埋め込み、また抽出を行うので、モノラル音響信号であってもステレオ音響信号であっても同様に行うことができる。
上記実施形態においては、埋め込み装置、抽出装置のいずれにおいても、左右のチャンネルを有するステレオ音響信号の左チャンネル信号に付加情報を埋め込む場合を例にとって説明したが、逆に右チャンネル信号に付加情報を埋め込むようにしても良い。本発明は、左右の特性には無関係だからである。また、１つのチャンネルしかないモノラル音響信号に対して処理を行う場合は、上記実施形態において、左チャンネル信号に対して行った処理を行うことになる。本発明は、１つのチャンネル信号に対して付加情報を埋め込み、また抽出を行うので、モノラル音響信号であってもステレオ音響信号であっても同様に行うことができる。
（９．信号成分が小さくても情報の埋め込みを可能とする手法）
ここまで、説明してきた処理においては、上記Ｅ1〜Ｅ4が所定値より大きいことが必要となり、上記Ｅ1〜Ｅ4が所定値以下の場合には、情報の埋め込みを行うことができない。そこで、以下、上記Ｅ1〜Ｅ4が所定値以下であっても、信号の埋め込みを可能とする手法について説明する。
ここまで、説明してきた処理においては、上記Ｅ1〜Ｅ4が所定値より大きいことが必要となり、上記Ｅ1〜Ｅ4が所定値以下の場合には、情報の埋め込みを行うことができない。そこで、以下、上記Ｅ1〜Ｅ4が所定値以下であっても、信号の埋め込みを可能とする手法について説明する。
この場合、図３に示した埋め込み装置における情報の埋め込み処理は、図６のフローチャートに従って行われるが、合算値Ｅ1〜Ｅ4がレベル下限値Ｌｅｖ以下となることがなくなるため、Ｓ１０９において低周波成分変更手段３０がレベルの判定を行わない。これは、この処理では、信号レベルが小さくても強制的に情報を埋め込むこととしているため、情報の埋め込みが不可能な信号レベルが小さい部分が存在するかどうかを判断し、区切りモードに設定する必要がないからである。
したがって、Ｓ１１１における状態１〜状態８に設定する処理としては、まず、以下の〔数式２３〕に従って算出される固定値Ｖを、低周波成分の強度として設定する。
〔数式２３〕
Ｖ＝｛０．５・Ｌｅｖ／（２Ｐ−６）｝1/2
Ｖ＝｛０．５・Ｌｅｖ／（２Ｐ−６）｝1/2
そして、状態１とする場合は、上記〔数式８〕に従った処理を実行した後、以下の〔数式２４〕に従った処理を実行する。
〔数式２４〕
ｊ＝ｍ〜ｍ＋Ｐ−４、ｍ＋Ｐ〜ｍ＋２Ｐ−４の各成分に対して
Ａｌ´（１，ｊ）＝Ａｌ（１，ｊ）・Ｖ／｛Ａｌ（１，ｊ）2＋Ｂｌ（１，ｊ）2｝1/2
Ｂｌ´（１，ｊ）＝Ｂｌ（１，ｊ）・Ｖ／｛Ａｌ（１，ｊ）2＋Ｂｌ（１，ｊ）2｝1/2
ｊ＝ｍ〜ｍ＋Ｐ−４、ｍ＋Ｐ〜ｍ＋２Ｐ−４の各成分に対して
Ａｌ´（１，ｊ）＝Ａｌ（１，ｊ）・Ｖ／｛Ａｌ（１，ｊ）2＋Ｂｌ（１，ｊ）2｝1/2
Ｂｌ´（１，ｊ）＝Ｂｌ（１，ｊ）・Ｖ／｛Ａｌ（１，ｊ）2＋Ｂｌ（１，ｊ）2｝1/2
状態２とする場合は、上記〔数式９〕に従った処理を実行した後、以下の〔数式２５〕に従った処理を実行する。
〔数式２５〕
ｊ＝ｍ〜ｍ＋Ｐ−４、ｍ＋Ｐ〜ｍ＋２Ｐ−４の各成分に対して
Ａｌ´（３，ｊ）＝Ａｌ（３，ｊ）・Ｖ／｛Ａｌ（３，ｊ）2＋Ｂｌ（３，ｊ）2｝1/2
Ｂｌ´（３，ｊ）＝Ｂｌ（３，ｊ）・Ｖ／｛Ａｌ（３，ｊ）2＋Ｂｌ（３，ｊ）2｝1/2
ｊ＝ｍ〜ｍ＋Ｐ−４、ｍ＋Ｐ〜ｍ＋２Ｐ−４の各成分に対して
Ａｌ´（３，ｊ）＝Ａｌ（３，ｊ）・Ｖ／｛Ａｌ（３，ｊ）2＋Ｂｌ（３，ｊ）2｝1/2
Ｂｌ´（３，ｊ）＝Ｂｌ（３，ｊ）・Ｖ／｛Ａｌ（３，ｊ）2＋Ｂｌ（３，ｊ）2｝1/2
状態３とする場合は、上記〔数式１０〕に従った処理を実行した後、以下の〔数式２６〕に従った処理を実行する。
〔数式２６〕
ｊ＝ｍ〜ｍ＋Ｐ−４の各成分に対して
Ａｌ´（１，ｊ）＝Ａｌ（１，ｊ）・Ｖ／｛Ａｌ（１，ｊ）2＋Ｂｌ（１，ｊ）2｝1/2
Ｂｌ´（１，ｊ）＝Ｂｌ（１，ｊ）・Ｖ／｛Ａｌ（１，ｊ）2＋Ｂｌ（１，ｊ）2｝1/2
ｊ＝ｍ＋Ｐ〜ｍ＋２Ｐ−４の各成分に対して
Ａｌ´（３，ｊ）＝Ａｌ（３，ｊ）・Ｖ／｛Ａｌ（３，ｊ）2＋Ｂｌ（３，ｊ）2｝1/2
Ｂｌ´（３，ｊ）＝Ｂｌ（３，ｊ）・Ｖ／｛Ａｌ（３，ｊ）2＋Ｂｌ（３，ｊ）2｝1/2
ｊ＝ｍ〜ｍ＋Ｐ−４の各成分に対して
Ａｌ´（１，ｊ）＝Ａｌ（１，ｊ）・Ｖ／｛Ａｌ（１，ｊ）2＋Ｂｌ（１，ｊ）2｝1/2
Ｂｌ´（１，ｊ）＝Ｂｌ（１，ｊ）・Ｖ／｛Ａｌ（１，ｊ）2＋Ｂｌ（１，ｊ）2｝1/2
ｊ＝ｍ＋Ｐ〜ｍ＋２Ｐ−４の各成分に対して
Ａｌ´（３，ｊ）＝Ａｌ（３，ｊ）・Ｖ／｛Ａｌ（３，ｊ）2＋Ｂｌ（３，ｊ）2｝1/2
Ｂｌ´（３，ｊ）＝Ｂｌ（３，ｊ）・Ｖ／｛Ａｌ（３，ｊ）2＋Ｂｌ（３，ｊ）2｝1/2
状態４とする場合は、上記〔数式１１〕に従った処理を実行した後、以下の〔数式２７〕に従った処理を実行する。
〔数式２７〕
ｊ＝ｍ〜ｍ＋Ｐ−４の各成分に対して
Ａｌ´（３，ｊ）＝Ａｌ（３，ｊ）・Ｖ／｛Ａｌ（３，ｊ）2＋Ｂｌ（３，ｊ）2｝1/2
Ｂｌ´（３，ｊ）＝Ｂｌ（３，ｊ）・Ｖ／｛Ａｌ（３，ｊ）2＋Ｂｌ（３，ｊ）2｝1/2
ｊ＝ｍ＋Ｐ〜ｍ＋２Ｐ−４の各成分に対して
Ａｌ´（１，ｊ）＝Ａｌ（１，ｊ）・Ｖ／｛Ａｌ（１，ｊ）2＋Ｂｌ（１，ｊ）2｝1/2
Ｂｌ´（１，ｊ）＝Ｂｌ（１，ｊ）・Ｖ／｛Ａｌ（１，ｊ）2＋Ｂｌ（１，ｊ）2｝1/2
ｊ＝ｍ〜ｍ＋Ｐ−４の各成分に対して
Ａｌ´（３，ｊ）＝Ａｌ（３，ｊ）・Ｖ／｛Ａｌ（３，ｊ）2＋Ｂｌ（３，ｊ）2｝1/2
Ｂｌ´（３，ｊ）＝Ｂｌ（３，ｊ）・Ｖ／｛Ａｌ（３，ｊ）2＋Ｂｌ（３，ｊ）2｝1/2
ｊ＝ｍ＋Ｐ〜ｍ＋２Ｐ−４の各成分に対して
Ａｌ´（１，ｊ）＝Ａｌ（１，ｊ）・Ｖ／｛Ａｌ（１，ｊ）2＋Ｂｌ（１，ｊ）2｝1/2
Ｂｌ´（１，ｊ）＝Ｂｌ（１，ｊ）・Ｖ／｛Ａｌ（１，ｊ）2＋Ｂｌ（１，ｊ）2｝1/2
状態５とする場合は、上記〔数式１２〕に従った処理を実行した後、以下の〔数式２８〕に従った処理を実行する。
〔数式２８〕
ｊ＝ｍ〜ｍ＋Ｐ−４の各成分に対して
Ａｌ´（１，ｊ）＝Ａｌ（１，ｊ）・Ｖ／｛Ａｌ（１，ｊ）2＋Ｂｌ（１，ｊ）2｝1/2
Ｂｌ´（１，ｊ）＝Ｂｌ（１，ｊ）・Ｖ／｛Ａｌ（１，ｊ）2＋Ｂｌ（１，ｊ）2｝1/2
ｊ＝ｍ〜ｍ＋Ｐ−４の各成分に対して
Ａｌ´（１，ｊ）＝Ａｌ（１，ｊ）・Ｖ／｛Ａｌ（１，ｊ）2＋Ｂｌ（１，ｊ）2｝1/2
Ｂｌ´（１，ｊ）＝Ｂｌ（１，ｊ）・Ｖ／｛Ａｌ（１，ｊ）2＋Ｂｌ（１，ｊ）2｝1/2
状態６とする場合は、上記〔数式１３〕に従った処理を実行した後、以下の〔数式２９〕に従った処理を実行する。
〔数式２９〕
ｊ＝ｍ〜ｍ＋Ｐ−４の各成分に対して
Ａｌ´（３，ｊ）＝Ａｌ（３，ｊ）・Ｖ／｛Ａｌ（３，ｊ）2＋Ｂｌ（３，ｊ）2｝1/2
Ｂｌ´（３，ｊ）＝Ｂｌ（３，ｊ）・Ｖ／｛Ａｌ（３，ｊ）2＋Ｂｌ（３，ｊ）2｝1/2
ｊ＝ｍ〜ｍ＋Ｐ−４の各成分に対して
Ａｌ´（３，ｊ）＝Ａｌ（３，ｊ）・Ｖ／｛Ａｌ（３，ｊ）2＋Ｂｌ（３，ｊ）2｝1/2
Ｂｌ´（３，ｊ）＝Ｂｌ（３，ｊ）・Ｖ／｛Ａｌ（３，ｊ）2＋Ｂｌ（３，ｊ）2｝1/2
状態７とする場合は、上記〔数式１４〕に従った処理を実行した後、以下の〔数式３０〕に従った処理を実行する。
〔数式３０〕
ｊ＝ｍ＋Ｐ〜ｍ＋２Ｐ−４の各成分に対して
Ａｌ´（３，ｊ）＝Ａｌ（３，ｊ）・Ｖ／｛Ａｌ（３，ｊ）2＋Ｂｌ（３，ｊ）2｝1/2
Ｂｌ´（３，ｊ）＝Ｂｌ（３，ｊ）・Ｖ／｛Ａｌ（３，ｊ）2＋Ｂｌ（３，ｊ）2｝1/2
ｊ＝ｍ＋Ｐ〜ｍ＋２Ｐ−４の各成分に対して
Ａｌ´（３，ｊ）＝Ａｌ（３，ｊ）・Ｖ／｛Ａｌ（３，ｊ）2＋Ｂｌ（３，ｊ）2｝1/2
Ｂｌ´（３，ｊ）＝Ｂｌ（３，ｊ）・Ｖ／｛Ａｌ（３，ｊ）2＋Ｂｌ（３，ｊ）2｝1/2
状態８とする場合は、上記〔数式１５〕に従った処理を実行した後、以下の〔数式３１〕に従った処理を実行する。
〔数式３１〕
ｊ＝ｍ＋Ｐ〜ｍ＋２Ｐ−４の各成分に対して
Ａｌ´（１，ｊ）＝Ａｌ（１，ｊ）・Ｖ／｛Ａｌ（１，ｊ）2＋Ｂｌ（１，ｊ）2｝1/2
Ｂｌ´（１，ｊ）＝Ｂｌ（１，ｊ）・Ｖ／｛Ａｌ（１，ｊ）2＋Ｂｌ（１，ｊ）2｝1/2
ｊ＝ｍ＋Ｐ〜ｍ＋２Ｐ−４の各成分に対して
Ａｌ´（１，ｊ）＝Ａｌ（１，ｊ）・Ｖ／｛Ａｌ（１，ｊ）2＋Ｂｌ（１，ｊ）2｝1/2
Ｂｌ´（１，ｊ）＝Ｂｌ（１，ｊ）・Ｖ／｛Ａｌ（１，ｊ）2＋Ｂｌ（１，ｊ）2｝1/2
〔数式２４〕〜〔数式３１〕により算出されたＡｌ´、Ｂｌ´の値を、上記〔数式７〕におけるＡｌ、Ｂｌに代えて用いると、Ｅ1〜Ｅ4＝Ｌｅｖとなるため、Ｓ１０９において、Ｌｅｖとの大小関係の判断を行う必要が不要となる。
上記のように、周波数成分が小さい場合に情報を埋め込んだ場合であっても、抽出側の、音響信号からの情報の抽出装置の構成は図９と同一であり、処理動作は図１０のフローチャートに従ったものと同一である。
１０・・・音響フレーム読込手段
２０・・・周波数変換手段
３０・・・低周波成分変更手段
４０・・・周波数逆変換手段
５０・・・改変音響フレーム出力手段
６０・・・記憶手段
６１・・・音響信号記憶部
６２・・・付加情報記憶部
６３・・・改変音響信号記憶部
７０・・・付加情報読込手段
１００・・・音響信号入力手段
１１０・・・基準フレーム獲得手段
１２０・・・位相変更フレーム設定手段
１３０・・・周波数変換手段
１４０・・・符号判定パラメータ算出手段
１５０・・・符号出力手段
１６０・・・付加情報抽出手段
１７０・・・音響フレーム保持手段
２０・・・周波数変換手段
３０・・・低周波成分変更手段
４０・・・周波数逆変換手段
５０・・・改変音響フレーム出力手段
６０・・・記憶手段
６１・・・音響信号記憶部
６２・・・付加情報記憶部
６３・・・改変音響信号記憶部
７０・・・付加情報読込手段
１００・・・音響信号入力手段
１１０・・・基準フレーム獲得手段
１２０・・・位相変更フレーム設定手段
１３０・・・周波数変換手段
１４０・・・符号判定パラメータ算出手段
１５０・・・符号出力手段
１６０・・・付加情報抽出手段
１７０・・・音響フレーム保持手段
Claims (10)
- 時系列のサンプル列で構成される音響信号に対して、付加情報を聴取不能な状態で埋め込む装置であって、
前記音響信号より、所定数のサンプルを音響フレームとして読み込む音響フレーム読込手段と、
前記音響フレームに対して第１窓関数、第２窓関数、第３窓関数を用いてそれぞれ周波数変換を行い、前記第１窓関数に対応するスペクトルである第１窓スペクトル、前記第２窓関数に対応するスペクトルである第２窓スペクトル、前記第３窓関数に対応するスペクトルである第３窓スペクトルを生成する周波数変換手段と、
前記生成された第１窓スペクトルから互いに重複しない２つの低周波数範囲のスペクトル集合ＳＰ１ＵとＳＰ１Ｄを抽出するとともに、前記第３窓スペクトルから互いに重複しない２つの低周波数範囲のスペクトル集合ＳＰ３ＵとＳＰ３Ｄを抽出し、前記埋め込むべき付加情報の情報配列および当該情報配列の付加情報内の１単位における位置に基づいて、ＳＰ１Ｕ、ＳＰ１Ｄ、ＳＰ３Ｕ、ＳＰ３Ｄの割合を変更すると共に、中央の第２窓関数に対応する低周波成分を除去する低周波成分変更手段と、
前記変更された低周波スペクトルを含む各窓スペクトルに対して周波数逆変換を行って、改変音響フレームを生成する周波数逆変換手段と、
前記生成された改変音響フレームを順次出力する改変音響フレーム出力手段と、
を有することを特徴とする音響信号に対する情報の埋め込み装置。 - 請求項１において、
前記低周波成分変更手段は、所定の固定値Ｖを利用して、ＳＰ１Ｕ、ＳＰ１Ｄ、ＳＰ３Ｕ、ＳＰ３Ｄのいずれかに所定の強度を与えることにより、前記ＳＰ１Ｕ、ＳＰ１Ｄ、ＳＰ３Ｕ、ＳＰ３Ｄのスペクトル強度の割合を変更するものであることを特徴とする音響信号に対する情報の埋め込み装置。 - 請求項１または請求項２において、
前記音響フレーム読込手段は、先行する音響フレームと所定数のサンプルを重複させて読み込み、読み込んだ音響フレーム全体に所定の窓関数を乗じて前記周波数変換手段に渡すものであり、
前記改変音響フレーム出力手段は、前記生成された改変音響フレームを先行する改変音響フレームと連結させて出力するものであることを特徴とする音響信号に対する情報の埋め込み装置。 - 請求項１から請求項３のいずれかにおいて、
前記低周波成分変更手段は、前記低周波数範囲をＦ１以上およびＦ２以下に設定し、前記スペクトル集合ＳＰ１ＤおよびＳＰ３Ｄの下限周波数はＦ１に一致させて設定し、前記スペクトル集合ＳＰ１ＤおよびＳＰ３Ｄの上限周波数はＳＰ１ＵおよびＳＰ３Ｕの下限周波数より所定の幅だけ小さく設定し、前記ＳＰ１ＵおよびＳＰ３Ｕの上限周波数はＦ２より所定の幅だけ小さく設定し、直流成分から前記ＳＰ１ＤおよびＳＰ３Ｄの下限周波数までのギャップ領域をＧ１ＤおよびＧ３Ｄとし、ＳＰ１ＤおよびＳＰ３Ｄの上限周波数とＳＰ１ＵおよびＳＰ３Ｕの下限周波数までのギャップ領域をＧ１ＭおよびＧ３Ｍとし、ＳＰ１ＵおよびＳＰ３Ｕの上限周波数からＦ２までのギャップ領域をＧ１ＵおよびＧ３Ｕとするとき、前記全てのギャップ領域の成分を常に除去するものであることを特徴とする音響信号に対する情報の埋め込み装置。 - 請求項１から請求項４のいずれかにおいて、
前記低周波成分変更手段は、前記Ｆ１を１５０Ｈｚ、前記Ｆ２を４００Ｈｚに設定するものであることを特徴とする音響信号に対する情報の埋め込み装置。 - 請求項１または請求項４において、
前記音響信号が左右２チャンネルの時系列のサンプル列で構成されるステレオ音響信号であって、
前記音響フレーム読込手段は、各チャンネルに対応する音響フレームをそれぞれ読み込み、
前記周波数変換手段は、各チャンネルの音響フレームに対して周波数変換を行い、各チャンネル別の第１窓スペクトル、第２窓スペクトル、第３窓スペクトルを生成するものであり、
前記低周波成分変更手段は、一方のチャンネルの窓スペクトルから抽出した所定の低周波数範囲に対応する各スペクトル集合に対して、前記埋め込むべき付加情報の情報配列および当該情報配列の付加情報内の１単位における位置に基づいて、ＳＰ１Ｕ、ＳＰ１Ｄ、ＳＰ３Ｕ、ＳＰ３Ｄの割合を変更するにあたり、ＳＰ１Ｕ、ＳＰ１Ｄ、ＳＰ３Ｕ、ＳＰ３Ｄのいずれかの成分を除去する際、前記一方のチャンネルにおいて行われた変更により除去された成分を補足するよう他方のチャンネルの前記スペクトル集合間の割合を変更するものであり、
前記周波数逆変換手段は、各チャンネルについて、変更されたスペクトル集合を含むフレームスペクトルに対して周波数逆変換を行って改変音響フレームを生成し、前記改変音響フレーム出力手段は、各チャンネルについて、生成された改変音響フレームを順次出力するものであることを特徴とする音響信号に対する情報の埋め込み装置。 - コンピュータを、請求項１から請求項６のいずれかに記載の音響信号に対する情報の埋め込み装置として、機能させるためのプログラム。
- 音響信号から、あらかじめ聴取不能な状態で埋め込まれた付加情報を抽出する装置であって、
前記音響信号の所定区間をデジタル化して、所定数のサンプルで構成される音響フレームを獲得する音響フレーム獲得手段と、
前記音響フレームに対して第１窓関数、第３窓関数を用いてそれぞれ周波数変換を行い、前記第１窓関数に対応するスペクトルである第１窓スペクトル、前記第３窓関数に対応するスペクトルである第３窓スペクトルを生成する周波数変換手段と、
前記生成された第１窓スペクトルから互いに重複しない２つの低周波数範囲のスペクトル集合ＳＰ１ＵとＳＰ１Ｄを抽出するとともに、前記第３窓スペクトルから互いに重複しない２つの低周波数範囲のスペクトル集合ＳＰ３ＵとＳＰ３Ｄを抽出し、４種類のスペクトル集合ＳＰ１Ｕ、ＳＰ１Ｄ、ＳＰ３Ｕ、ＳＰ３Ｄのスペクトル強度を算出し、その各スペクトル強度の大きさおよび２種類のスペクトル強度の大小関係に基づいて、埋め込まれていた付加情報の情報配列および当該情報配列の付加情報内の１単位における位置を示す所定の符号を出力する符号化手段と、
前記出力された符号に対応する情報配列を、所定の規則により変換して付加情報を抽出する付加情報抽出手段と、
を有することを特徴とする音響信号からの情報の抽出装置。 - 請求項８において、
前記音響フレーム獲得手段は、前記音響信号から、所定数のサンプルで構成される音響フレームを基準フレームとして獲得する基準フレーム獲得手段と、前記基準フレームと所定サンプルずつ移動させることにより位相を変更して設定される複数の音響フレームを位相変更フレームとして設定する位相変更フレーム設定手段により構成され、
前記符号化手段は、前記抽出したスペクトル集合に基づいて、符号判定パラメータを算出する符号判定パラメータ算出手段と、基準フレームが異なる過去の同位相の音響フレームにおいて算出された符号判定パラメータに基づいて、前記基準フレームおよび複数の位相変更フレームのうち１つの音響フレームを位相が最適なものであると判断し、当該最適な位相の音響フレームについて判断された前記符号判定パラメータに基づいて、所定の符号を出力する符号出力手段を有するものであることを特徴とする音響信号からの情報の抽出装置。 - コンピュータを、請求項８または請求項９に記載の音響信号からの情報の抽出装置として、機能させるためのプログラム。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2006189196A JP5011855B2 (ja) | 2006-07-10 | 2006-07-10 | 音響信号に対する情報の埋め込み装置および音響信号からの情報の抽出装置 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2006189196A JP5011855B2 (ja) | 2006-07-10 | 2006-07-10 | 音響信号に対する情報の埋め込み装置および音響信号からの情報の抽出装置 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JP2008015416A true JP2008015416A (ja) | 2008-01-24 |
| JP5011855B2 JP5011855B2 (ja) | 2012-08-29 |
Family
ID=39072464
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP2006189196A Expired - Fee Related JP5011855B2 (ja) | 2006-07-10 | 2006-07-10 | 音響信号に対する情報の埋め込み装置および音響信号からの情報の抽出装置 |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JP5011855B2 (ja) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2010026242A (ja) * | 2008-07-18 | 2010-02-04 | Yamaha Corp | 電子透かし情報の埋め込みおよび抽出を行う装置、方法およびプログラム |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2002044408A (ja) * | 2000-06-29 | 2002-02-08 | Internatl Business Mach Corp <Ibm> | 電子透かし方法およびそのシステム |
| JP2008009104A (ja) * | 2006-06-29 | 2008-01-17 | Dainippon Printing Co Ltd | 音響信号に対する情報の埋め込み装置および音響信号からの情報の抽出装置 |
-
2006
- 2006-07-10 JP JP2006189196A patent/JP5011855B2/ja not_active Expired - Fee Related
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2002044408A (ja) * | 2000-06-29 | 2002-02-08 | Internatl Business Mach Corp <Ibm> | 電子透かし方法およびそのシステム |
| JP2008009104A (ja) * | 2006-06-29 | 2008-01-17 | Dainippon Printing Co Ltd | 音響信号に対する情報の埋め込み装置および音響信号からの情報の抽出装置 |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2010026242A (ja) * | 2008-07-18 | 2010-02-04 | Yamaha Corp | 電子透かし情報の埋め込みおよび抽出を行う装置、方法およびプログラム |
Also Published As
| Publication number | Publication date |
|---|---|
| JP5011855B2 (ja) | 2012-08-29 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| JP4899416B2 (ja) | ネットワーク接続装置 | |
| JP4839775B2 (ja) | 音響信号に対する情報の埋め込み装置、方法、プログラム | |
| JP4660275B2 (ja) | 音響信号に対する情報の埋め込み装置および方法 | |
| JP5011849B2 (ja) | 音響信号に対する情報の埋め込み装置および音響信号からの情報の抽出装置 | |
| JP4867765B2 (ja) | 音響信号に対する情報の埋め込み装置および音響信号からの情報の抽出装置 | |
| JP4629495B2 (ja) | 音響信号に対する情報の埋め込み装置および方法 | |
| JP2007292827A (ja) | 音響信号検索装置 | |
| JP4531653B2 (ja) | 音響信号からの情報の抽出装置 | |
| JP5082257B2 (ja) | 音響信号検索装置 | |
| JP5011855B2 (ja) | 音響信号に対する情報の埋め込み装置および音響信号からの情報の抽出装置 | |
| JP4831335B2 (ja) | 音響信号に対する情報の埋め込み装置および音響信号からの情報の抽出装置 | |
| JP5011872B2 (ja) | 音響信号に対する情報の埋め込み装置および音響信号からの情報の抽出装置 | |
| JP4831333B2 (ja) | 音響信号に対する情報の埋め込み装置および音響信号からの情報の抽出装置 | |
| JP4713181B2 (ja) | 音響信号に対する情報の埋め込み装置、音響信号からの情報の抽出装置、および音響信号再生装置 | |
| JP4713180B2 (ja) | 音響信号からの情報の抽出装置 | |
| JP5011865B2 (ja) | 音響信号に対する情報の埋め込み装置および音響信号からの情報の抽出装置 | |
| JP4831334B2 (ja) | 音響信号に対する情報の埋め込み装置および音響信号からの情報の抽出装置 | |
| JP4876978B2 (ja) | 音響信号に対する情報の埋め込み装置および音響信号からの情報の抽出装置 | |
| JP5003164B2 (ja) | 音響信号からの情報の抽出装置 | |
| JP2008129198A (ja) | 音響信号に対する情報の埋め込み装置、音響信号からの情報の抽出装置 | |
| JP4760540B2 (ja) | 音響信号に対する情報の埋め込み装置 | |
| JP4760539B2 (ja) | 音響信号に対する情報の埋め込み装置 | |
| JP2008139491A (ja) | 音響信号に対する情報の埋め込み装置および音響信号からの情報の抽出装置 | |
| JP2009075332A (ja) | 音響信号に対する情報の埋め込み装置および音響信号からの情報の抽出装置 | |
| JP4910921B2 (ja) | 音響信号に対する情報の埋め込み装置および音響信号からの情報の抽出装置 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20090625 |
|
| A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20110412 |
|
| A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20110609 |
|
| A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20110804 |
|
| 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: 20120508 |
|
| 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: 20120521 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20150615 Year of fee payment: 3 |
|
| R150 | Certificate of patent or registration of utility model |
Ref document number: 5011855 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R150 Free format text: JAPANESE INTERMEDIATE CODE: R150 |
|
| LAPS | Cancellation because of no payment of annual fees |