(57)【要約】受信された信号における時間オフセットを補正して、受信された信号のフレーム系列を送信されたシンボル系列に正しく配列するための方法及び装置が説明される。それぞれのシンボルは、ＴSの信号サンプルにわたり延びている。受信された信号は、はじめに、長さＴSのフレーム系列に分割され、次いでそれぞれのフレームは、Ｎbのサブフレームの倍数に分割される。続いて、Ｎbの値の系列が形成され、それぞれの系列における各連続する値は、それぞれ連続するフレーム内の対応するサブフレームから導出される。Ｎbの系列のそれぞれは、正しく配列された送信シンボル系列の予測値である。【特許請求の範囲】
【請求項１】
受信された信号におけるオフセットを補償する方法であって、該受信された信号は、それぞれのシンボルがＴSの信号サンプルにわたり延びるシンボルの系列により変更され、
該受信された信号をフレームに分割するステップと、
それぞれのフレームを複数のＮbのサブフレームに分割するステップと、
それぞれのフレーム内の対応するサブフレームから導出される値のＮbの系列を形成するステップと、
該Ｎbの系列を該シンボル系列に正しく配列された連続するフレーム系列の予測値として取得するステップと、
を有する方法。
【請求項２】
それぞれのフレームは、予め決定された長さＴSからなる、
請求項１記載の方法。
【請求項３】
フレーム間のオーバラップが存在する、
請求項１記載の方法。
【請求項４】
それぞれのサブフレームは、隣接するサブフレームとオーバラップする、
請求項１記載の方法。
【請求項５】
該Ｎbは２から８の範囲にある、
請求項１記載の方法。
【請求項６】
該シンボルの系列は、ＬWシンボルを有しており、該受信された信号は、ＬFフレームに分割され、該ＬFは該ＬWの整数倍である、
請求項１記載の方法。
【請求項７】
該シンボルの系列は、帯域制限された周波数応答及び平滑化された時間応答を有する窓成形関数で畳み込みされる値の系列を有する、
請求項１記載の方法。
【請求項８】
該窓成形関数は、対称な時間応答及び反対称な時間応答を有する、
請求項７記載の方法。
【請求項９】
該シンボルの系列は、レイズドコサイン関数又はバイフェーズ関数のうちの少なくとも一方の系列を有する、
請求項１記載の方法。
【請求項１０】
該オフセットは時間のずれである、
請求項１記載の方法。
【請求項１１】
どの予測値が最良の予測値であるかを判定するため、それぞれの予測値が正しく配列されたフレーム系列であるかのように該それぞれの予測値を処理するステップをさらに有する、
請求項１記載の方法。
【請求項１２】
該最良の予測値は、処理されたときに、１以上の予め決定された条件を超える第一の予測値であるとされ、該処理するステップは、該最良の予測値がひとたび決定されたときに停止する、
請求項１１記載の方法。
【請求項１３】
該予測値のそれぞれを該シンボルの系列に対応する基準値と相関をとるステップと、最大の相関ピーク値をもつ予測値を該最良の予測値として取得するステップとをさらに有する、
請求項１記載の方法。
【請求項１４】
第一の信号又は信号の一部について第一の最良の予測値がひとたび決定されると、該方法は、更に受信された信号又は信号の一部について繰り返され、該更なる信号からの予測値は、該第一の最良の予測値に依存する順序で処理される、
請求項１１乃至１３のいずれか記載の方法。
【請求項１５】
請求項１記載の方法を実行するために調整されるコンピュータプログラム。
【請求項１６】
請求項１５記載のコンピュータプログラムを有する記録キャリア。
【請求項１７】
請求項１５記載のコンピュータプログラムをダウンロードするために利用可能にされる方法。
【請求項１８】
受信された信号におけるオフセットを補償するために構成される装置であって、該受信された信号は、それぞれのシンボルがＴSの信号サンプルにわたり延びるシンボルの系列により変更され、
受信された信号をフレームに分割するために配置される分割手段と、
それぞれのフレームを複数のＮbのサブフレームに分割するために配置される分割手段と、
それぞれのフレーム内の対応するサブフレームから導出される値のＮbの系列を形成し、該Ｎbの系列を該シンボル系列に正しく配列された連続するフレーム系列の予測値として取得するために配置されるプロセッサと、
を有する装置。
【請求項１９】
該Ｎbの系列を記憶するために配置されるバッファをさらに有する、
請求項１８記載の装置。
【請求項２０】
請求項１８記載の装置を有するデコーダ。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、オーディオ、ビデオ又はデータ信号のような情報信号に埋め込まれた情報をデコードするための装置及び方法に関する。
【背景技術】
【０００２】
情報信号のウォーターマーキングは、更なるデータを情報信号と共に伝送するための技術である。たとえば、ウォーターマーキング技術は、著作権及びコピー制御情報をオーディオ信号に埋め込むために使用することができる。
【０００３】
ウォーターマーキングスキームの主要な要件は、信号からウォーターマークを取り除くための攻撃（たとえば、ウォーターマークを取り除くことは信号にダメージを与える）に対してロバストでありつつ、ウォーターマークが観察できない（すなわち、オーディオ信号のケースでは可聴ではない）ことである。ウォーターマークのロバスト性は、通常、ウォーターマークが埋め込まれる信号の品質に対してトレードオフにあることが理解される。たとえば、ウォーターマークがオーディオ信号に強固に埋め込まれている（及び取り除くことが難しい）場合、オーディオ信号の品質が低下される可能性が高い。
【０００４】
様々なタイプのオーディオのウォーターマーキングスキームが提案されており、それぞれのスキームは、それ自身の優れた点及び劣る点を有している。あるタイプのオーディオウォーターマーキングは、所望のデータ（たとえば、著作権情報）をオーディオ信号に埋め込むために時間的な補正技法を使用することである。この技術は、二次方程式を解くことでエコーの強度が決定される効果的なエコー隠しアルゴリズムである。二次方程式は、一方の位置がτに等しい遅延であり、他方の位置が０に等しい遅延である、２つの位置での自己相関値により生成される。デコーダでは、２つの遅延の位置での自己相関関数の比を決定することで、ウォーターマークが抽出される。
【０００５】
国際特許出願WO  00/00969は、Aris  Technologies社によるものであり、補助信号を情報ホスト又はカバー信号に埋め込み又はエンコードするための技術を開示している。カバー信号のレプリカ、又は特定のドメイン（時間、周波数又は空間）におけるカバー信号の一部は、ステゴ鍵に従って生成され、このステゴ鍵は、カバー信号のパラメータへの修正値を規定している。次いで、レプリカ信号は、埋め込まれるべき情報に対応する補助信号により修正され、ステゴ信号を形成するようにカバー信号に挿入される。
【０００６】
デコーダでは、オリジナルの補助データを抽出するため、オリジナルのカバー信号のレプリカと同様のやり方でステゴ信号のレプリカが生成される。次いで、結果的に得られるレプリカが受信されたステゴ信号と相関が取られ、補助信号が抽出される。補助信号の抽出は、比較的複雑であり、エンコーダ（又はエンベッダ）及びデコーダ（又はディテクタ）の両者でステゴ鍵を必要とする。さらに、ディテクタで補助信号に同期するために総当りのサーチが必要とされる。
【０００７】
さらに、ペイロード抽出の性能は、補助信号をどの位良好に予測することができるかに依存する。補助信号におけるペイロードのビットに関して高く期待されるエラーレートをもつシステムでは、これは達成することが難しい。このソリューションにより、非常に複雑な誤り訂正方法、又は情報のキャパシティの大幅な制限を招くことになる。
【発明の開示】
【発明が解決しようとする課題】
【０００８】
本発明の目的は、従来技術の少なくとも１つの問題に実質的に対処する、ウォーターマークのデコーディングスキームについて時間オフセットの補償を提供することにある。
【課題を解決するための手段】
【０００９】
第一の態様では、本発明は、受信された信号におけるオフセットを補償する方法を提供するものであり、該受信された信号はシンボル系列により変更される。それぞれのシンボルは、ＴSの信号サンプルにわたり延びている。本方法は、（ａ）受信された信号を予め決定された長さＴSのフレームに分割するステップ、（ｂ）それぞれのフレームを複数のＮbのサブフレームに分割するステップ、（ｃ）それぞれのフレーム内の対応するサブフレームから導出される値のＮbの系列を形成するステップ、及び（ｄ）該Ｎbの系列を、該シンボル系列に正しく配列された（オフセットがない）連続するフレーム系列の予測値として取得するステップを含んでいる。
【００１０】
好ましくは、それぞれのフレームは、隣接するフレームとオーバラップする。  
好ましくは、それぞれのサブフレームは、隣接するサブフレームとオーバラップする。  
好ましくは、Ｎbは、２から８の範囲にある。  
好ましくは、シンボル系列は、ＬWシンボルを備え、受信された信号は、ＬFフレームに分割され、この場合、ＬFは、ＴS・ＬWの整数倍である。
【００１１】
好ましくは、該系列のそれぞれのシンボルは、ＴSサンプルにわたり延びる。  
好ましくは、該シンボルは、帯域制限された周波数変化を有し、通過帯域で滑らかな（グレースフルな）時間変化を有する窓成形関数で成形される。さらに、窓成形関数は、対称的又は反対称的な時間変化を有していることが好ましい。  
好ましくは、該窓成形関数は、レイズドコサイン関数（raised  cosine  function）又はバイフェーズ関数（bi-phase  function）を有することが好ましい。
【００１２】
好ましくは、該オフセットは、受信された信号と送信された信号との間の時間のずれである。  
好ましくは、本方法は、どの予測値が最良の予測値であるかを決定するように、あたかもそれぞれの予測値が正しく配列されたフレーム系列であるかのように、ステップ（ｄ）で生成されたそれぞれの予測値を処理するステップをさらに有している。
【００１３】
本方法は、該予測値のそれぞれを該シンボル系列に対応する基準シーケンスと相関をとるステップ、及び最大の相関ピーク値をもつ予測値を最良の予測値として取得するステップをさらに有している。  
好ましくは、最良の予測値は、処理されたときに、１以上の予め決定された条件を超える最初の予測値であると考えられる。好ましくは、予測値の処理は、ワーキング予測値がひとたび決定されると停止する。  
好ましくは、第一の信号又は信号の一部についてワーキング予測値がひとたび決定されると、本方法は、更に受信された信号又は信号の一部について繰り返され、該更なる信号からの予測値は、該第一の最良の予測値に依存する順序で処理される。このように、本方法は、オフセットについて適応的に変更することができる。
【００１４】
更なる態様では、本発明は、先の方法を実行するために調整されるコンピュータプログラムを提供する。  
別の態様では、本発明は、先のコンピュータプログラムを有する記録キャリアを提供する。  
更なる態様では、本発明は、先のコンピュータプログラムをダウンロードするために利用可能にされる方法を提供する。
【００１５】
別の態様では、本発明は、受信された信号におけるオフセットを補償するために構成される装置を提供するものであって、該受信された信号は、シンボル系列により変更され、それぞれのシンボルは、ＴSの信号サンプルにわたり延びている。本装置は、受信された信号を好ましくは長さＴSのフレームに分割するために配置され、それぞれのフレームを複数のＮbのサブフレームに分割するために配置される分割手段、及びそれぞれのフレーム内の対応するサブフレームから導出されるＮbの値の系列を形成し、該Ｎbの系列をシンボル系列に正しく配列された（オフセットのない）連続するフレーム系列の予測値として取得するプロセッサを有している。
【００１６】
本発明に関する最良の理解のため、及び類似の実施の形態がどのように実施されるかを示すため、例を通して添付図面を参照しながら説明が行われる。
【発明を実施するための最良の形態】
【００１７】
図１は、マルチビットのペイロードのウォーターマークｗをホスト信号ｘに埋め込むためのデジタル信号処理を実行するために必要とされる装置のブロック図である。ホスト信号ｘは、装置の入力１２で提供される。ホスト信号ｘは、加算器２２を介して出力１４の方向に通過される。しかし、ホスト信号ｘのレプリカ（入力８）は、ウォーターマーク情報を伝送するために乗算器１８の方向に分割される。
【００１８】
ウォーターマーク信号ｗCは、ペイロードエンベッダ（embedder）及びウォーターマークコンディショニング装置６から得られ、ペイロードエンベッダ及びウォーターマークコンディショニング装置に入力される基準となる有限長のランダム系列ｗSから導出される。乗算器１８は、ウォーターマーク信号ｗCとレプリカオーディオ信号ｘとの積を計算するために利用される。次いで、結果的に得られる積ｗC・ｘは、ゲインコントローラ２４を介して加算器２２に通過される。ゲインコントローラ２４は、ゲインファクタαにより信号を増幅又は減衰するために使用される。
【００１９】
ゲインファクタαは、可聴性とウォーターマークのロバスト性との間のトレードオフを制御する。ゲインファクタは、定数であるか、又は時間、周波数及び空間の少なくとも１つに関する変数である場合がある。図１における装置は、αが変数であるとき、ホスト信号ｘの特性に基づいて、信号分析ユニット２６を介して自動的に調節することができる。好ましくは、オーディオ信号の場合にＨＡＳ（Human  Auditory  System）の精神音響モデル（psycho-acoustic  model）のような、適切に選択される知覚のコスト関数（perceptibility  cost  function）に従って信号品質への影響を最小にするように、利得αが自動的に調節される。かかるモデルは、たとえば、E.  Zwickerによる文献“Audio  Engineering  and  Psychoacoustics  :  Matching  signals  to  the  final  receiver,  Human  Auditory  System”Journal  of  the  Audio  Engineering  Society,  Vol.  39,  pp.  Vol.115-126,  March  1991に記載されている。
【００２０】
以下では、本発明のこの実施の形態を説明するために、例示のみを介して、オーディオウォーターマークが利用される。
【００２１】
その後、結果的に得られるウォーターマークオーディオ信号ｙは、ｗCとｘとの積の適切にスケーリングされたバージョンをホスト信号に加えることで、埋め込み装置１０の出力１４で得られる。
【数１】
2005522079.tif 000002
【００２２】
好ましくは、ウォーターマークｗCは、ｘで乗算されたとき、ｘの短時間のエンベロープを支配的に変更するように選択される。
【００２３】
図２は、１つの好適な実施の形態を示しており、この実施の形態では、図１における乗算器１８への入力８は、フィルタリングユニット１５におけるフィルタＨを使用してホスト信号ｘのレプリカをフィルタリングすることにより得られる。フィルタリング出力がｘbにより示される場合、この好適な実施の形態によれば、ウォーターマーク信号は、ｘbとウォーターマークｗCとの積をホスト信号ｘに加えることで生成される。
【数２】
2005522079.tif 000003
ここで、
（外１）
2005522079.tif 000004
は、
（外２）
2005522079.tif 000005
とｙbが、
（外３）
2005522079.tif 000006
であるように定義されるように定義され、ウォーターマーク信号ｙのエンベロープが変調された部分ｙbは、以下のように与えられる。
【数３】
2005522079.tif 000007
【００２４】
好ましくは、図３に示されるように、フィルタＨは、その下側の遮断周波数ｆL及び上側の遮断周波数ｆHにより特徴付けされる線形位相帯域通過フィルタである。図３（ｂ）でわかるように、フィルタＨは、通過帯域（ＢＷ）内で周波数ｆに関して、線形位相応答を有している。したがって、フィルタＨが帯域通過フィルタであるとき、ｘbはホスト信号の通過帯域成分の範囲内にあり、
（外４）
2005522079.tif 000008
は通過帯域成分の範囲外にある。最適なパフォーマンスのために、信号ｘb及び
（外５）
2005522079.tif 000009
は同相であることが好ましい。このことは、フィルタＨにより生じる位相歪みを適切に補償することで達成される。線形位相フィルタのケースでは、歪みは簡単な時間遅延である。
【００２５】
図４では、ペイロードエンベッダ及びウォーターマークコンディショニングユニット６の詳細が示されている。このユニットでは、最初の基準となるランダム系列ｗSは、マルチビットのウォーターマーク信号ｗCに変換される。
【００２６】
はじめに、以下ウォーターマークシード信号とも呼ばれる、有限長の、好ましくは、ゼロ平均かつ一様に分布されるランダム系列ｗSは、初期シードＳをもつランダム数発生器を使用して生成される。後に理解されるように、ウォーターマーク信号のコピーが補償処理のために検出器で生成することができるように、この初期シードＳがエンベッダ及び検出器の両者に知られていることが好ましい。これは、長さＬWの系列となる。
【数４】
2005522079.tif 000010
【００２７】
なお、ある応用では、シードは、代替的なチャネルを介してディテクタに送信することができるか、又はある予め決定されたプロトコルを使用して受信された信号から導出することができる。
【００２８】
その後、系列ＷSは、巡回シフトユニット３０を使用して量ｄ1及びｄ2だけ巡回的にシフトされ、ランダム系列ｗd1及びｗd2がそれぞれ取得される。これらの２つの系列（ｗd1及びｗd2）は、実際上は第一の系列及び第二の系列であり、第二の系列は、第一の系列に関して巡回的にシフトされている。それぞれの系列ｗdi、ｉ＝１，２，は、続いて乗算ユニット４０においてそれぞれの符号ビットｒiで乗算される。ここで、ｒi＝＋１又は−１である。それぞれのｒ1及びｒ2の値は、定数のままであり、ウォーターマークのペイロードが変化したときにのみ変化する。その後、それぞれの系列は、図４に示されるウォーターマークコンディショニング回路２０により長さＬW・ＴSの周期的、緩やかに変化する狭帯域信号ｗiに変換される。最後に、緩やかに変化する狭帯域信号ｗ1及びｗ2は、相対的な遅延Ｔr（この場合、Ｔr＜ＴS）を介して加算され、マルチビットのペイロードのウォーターマーク信号ｗCが与えられる。このことは、遅延ユニット４５を使用して量Ｔrだけ信号ｗ2をはじめに遅延し、続いて、この遅延された信号ｗ2をｗ1に加算ユニット５０で加えることで達成される。
【００２９】
図５は、ペイロードエンベッダ及びウォーターマークコンディショニング装置６で使用されるウォーターマークコンディショニング装置２０をより詳細に示している。ウォーターマークシード信号ｗSは、コンディショニング装置２０に入力される。
【００３０】
便宜上、唯一の系列Ｗdiの変更は、図５に示されているが、それぞれの系列は類似のやり方で変更され、結果は加算されてウォーターマーク信号ｗCが取得される。
【００３１】
図５に示されるように、それぞれのウォーターマーク信号の系列ｗdi[k]、ｉ＝１，２は、サンプルリピータ１８０の入力に印加される。チャート１８１は、＋１と−１の間のランダム数の値の系列として系列ｗdiのうちの１つを例示しており、系列は長さＬWからなる。サンプルリピータは、ウォーターマークシード信号系列内でそれぞれの値をＴS回繰返し、矩形のパルス列の信号を生成する。ＴSは、ウォーターマークシンボル周期と呼ばれ、オーディオ信号におけるウォーターマークシンボルのスパンを表している。チャート１８３は、チャート１８１で例示される信号がサンプルリピータ１８０をひとたび通過したときの該信号の結果を表している。
【００３２】
次いで、レイズドコサイン窓のような窓成形関数ｓ[n]は、ｗd1及びｗd2から導出された矩形のパルス関数を緩やかに変化するウォーターマーク系列の関数ｗ1[n]及びｗ2[n]にそれぞれ変換するために印加される。
【００３３】
チャート１８４は、典型的なレイズドコサイン窓成形関数を示しており、この関数もまたスパンＴSである。
【００３４】
次いで、生成されたウォーターマーク系列ｗ1[n]及びｗ2[n]は、相対的な遅延Ｔr（この場合、Ｔr＜ＴS）により加算され、マルチビットのペイロードのウォーターマーク信号ｗC[n]が与えられる。すなわち、
【数５】
2005522079.tif 000011
Ｔrの値は、ｗ1のゼロ交差がｗ2の最大の振幅ポイントに一致し、逆にｗ2のゼロ交差がｗ1の最大の振幅ポイントに一致するように選択される。したがって、レイズドコサイン窓成形関数についてＴr＝ＴS／２であり、バイフェーズ窓成形関数についてＴr＝ＴS／４である。他の窓成形関数について、Ｔrの他の値も可能である。
【００３５】
以下の説明により理解されるように、検出の間、ｗC[n]の相関は、（図１２に見ることができるように）ｐＬ’により分離される２つの相関ピークを生成する。ｐＬ’は、ペイロードの一部である、ｗd1とｗd2との間の巡回的なシフトｐＬの予測値であり、以下のように定義される。
【数６】
2005522079.tif 000012
【００３６】
ｐＬに加えて、埋め込まれたウォーターマークの相対的な符号を変えることで、余分の情報をエンコードすることができる。
【００３７】
検出器では、これは、相関ピークの間で相対的な符号ｒsignとして見ることができる。この相対的な符号は、以下のように定義される。
【数７】
2005522079.tif 000013
この場合、ρ1＝ｓｉｇｎ（ｃＬ1）及びρ2＝ｓｉｇｎ（ｃＬ2）は、図４の符号ビットｒ1（入力８０）及びｒ2（入力９０）のそれぞれの予測値であり、ｃＬ1及びｃＬ2は、ｗd1及びｗd2のそれぞれに対応する相関ピークの値である。全体のウォーターマークのペイロードｐＬWは、エラーのない検出のためにｒsignとｐＬとの結合として与えられる。
【数８】
2005522079.tif 000014
【００３８】
長さＬWのウォーターマーク系列により搬送することができる最大の情報（Ｉmax）は、ビット数として、以下により与えられる。
【数９】
2005522079.tif 000015
【００３９】
かかるスキームでは、ペイロードは、エンベッダと検出器との間の相対的なオフセットに対して免疫性があり、可能性のある時間スケールの変化に対しても免疫性がある。
【００４０】
窓成形関数は、本発明のウォーターマーキングスキームのロバスト性及び可聴性の応答を制御する主要なパラメータのうちの１つとして識別される。図６ａ及び図６ｂに例示されるように、本実施の形態では、窓成形関数の２つの可能な例である、レイズドコサイン関数及びバイフェーズ関数が説明される。
【００４１】
擬似ＤＣフリーなウォーターマーク信号を得るように、レイズドコサイン窓関数の代わりに、バイフェーズ窓関数を使用することが好ましい。これは、図７ａ及び図７ｂに例示されており、レイズドコサイン窓成形関数及びバイフェーズ窓成形関数のそれぞれで条件付けされるウォーターマーク系列（このケースでは、ｗdi[k]＝{1,1,-1,1,-1,-1}の系列）に対応する周波数スペクトルを示している。理解できるように、レイズドコサインで条件付けされたウォーターマーク系列の周波数スペクトルは、周波数ｆ＝０で最大値を有し、バイフェーズ成形されたウォーターマーク系列の周波数スペクトルは、周波数ｆ＝０で最小値を有しており、すなわち、非常に小さなＤＣ成分を有する。
【００４２】
有効な情報は、ウォーターマークのＤＣ成分ではない成分にのみ含まれている。結果的に、同じく加えられたウォーターマークのエネルギーについて、バイフェーズ窓で条件付けされたウォーターマークは、レイズドコサイン窓により条件付けされたウォーターマークよりも有効な情報を搬送する。結果として、バイフェーズ窓は、同じロバスト性について優れた可聴性のパフォーマンスを提供し、逆に、同じ可聴性の品質について良好なロバスト性を可能にする。
【００４３】
かかるバイフェーズ関数は、他のウォーターマーキングスキームについて窓成形関数として利用することができる。言い換えれば、バイフェーズ関数は、別の信号に組み込まれることになる（ウォーターマークのような）信号のＤＣ成分を低減するために印加される。
【００４４】
図８は、ウォーターマーク検出器２００，３００，４００のブロック図を示している。検出器は、３つの主要なステージから構成される。（ａ）ウォーターマークシンボル抽出ステージ２００、（ｂ）バッファリング及び補間ステージ３００、及び（ｃ）相関及び判定ステージ４００。
【００４５】
シンボル抽出ステージ２００では、受信されたウォーターマーク信号ｙ’[n]は、ウォーターマークが埋め込まれた系列の多数（Ｎb）の予測値を生成するために処理される。これらのウォーターマーク系列の予測値は、エンベッダと検出器との間に存在する場合がある時間オフセットを解決するために必要とされ、ウォーターマーク検出器は、ホスト信号に挿入されたウォーターマーク系列に同期することができる。
【００４６】
バッファリング及び補間ステージ３００では、これらの予測値は、Ｎb個の個別のバッファに分解され、補間は、生じる場合がある時間スケールの変更を解決するためにそれぞれのバッファに印加される。たとえば、サンプリング（クロック）周波数におけるドリフトは、時間領域の信号における伸張又は収縮となる場合がある（すなわち、ウォーターマークが伸張又は収縮される場合がある）。
【００４７】
相関及び判定ステージ４００では、それぞれのバッファの内容は、基準のウォーターマークと相関が取られ、最大の相関ピークは、ウォーターマークが受信された信号ｙ’[n]内に確かに埋め込まれたかに関する可能性を判定するために閾値と比較される。
【００４８】
ウォーターマーク検出の精度を最大にするため、ウォーターマーク検出処理は、ウォーターマーク系列の長さの３倍から４倍の長さである受信された信号ｙ’[n]にわたり典型的に実行される。したがって、検出されるべきそれぞれのウォーターマークシンボルは、該シンボルの幾つかの予測値の平均を得ることで構築することができる。この平均化処理は、平滑化と呼ばれ、平均化が行われる回数は、スムージングファクタｓfと呼ばれる。ＬDを検出窓の長さとし、ウォーターマーク検出の真理値が報告される（サンプル数での）オーディオセグメントの長さとして定義される。このとき、ＬD＝ｓf・ＬW・ＴSであり、この場合、ＴSはシンボル周期であり、ＬWは、ウォーターマーク系列内でのシンボル数である。シンボル抽出の間、ファクタＴSの決定は、エネルギー計算ステージで行われる。このように、バッファリング及び補間ステージ内でのそれぞれのバッファ（Ｂｕｆｆｅｒ）３２０の長さ（Ｌb）はＬb＝ｓf・ＬWである。
【００４９】
図８に示されるウォーターマークのシンボル抽出ステージ２００では、到来するウォーターマーク信号ｙ’[n]は、任意の信号コンディショニングフィルタＨb２１０に入力される。このフィルタ２１０は、典型的には帯域通過フィルタであり、図２に示される対応するフィルタＨ１５と同じ応答を有している。フィルタＨbの出力は、ｙ’b[n]であり、伝送媒体での線形性を仮定して、式（１）及び（３）から以下に従う。
【数１０】
2005522079.tif 000016
【００５０】
なお、上式において、エンベッダと検出器との間で可能性のある時間オフセットは暗黙的に無視される。一般的なウォーターマーキングスキームの原理の説明の容易さのため、以下では、エンベッダと検出器との間に完全な同期が存在する（すなわち、オフセットがない）ことが仮定される。本発明に従って時間オフセットがどのように補正されるかに関して図１１を参照して説明が与えられる。
【００５１】
なお、エンベッダにおいてフィルタが使用されないとき（すなわち、Ｈ＝１であるとき）、検出器におけるＨbもまた省略することができるか、又は検出のパフォーマンスを改善するためになお含めることができる。Ｈbが省略される場合、式（１０）におけるｙbはｙで置き換えられる。残りの処理は同じである。
【００５２】
オーディオ信号は長さＴSのフレームに分割され、かつｙ’b,m[n]はｍ番目のフィルタリングされたフレーム信号のうちのｎ番目のサンプルであることが仮定される。ｍ番目のフレームに対応するエネルギーＥ[m]は、以下に示される。
【数１１】
2005522079.tif 000017
【００５３】
この式（１１）を式（１０）と組み合わせて、以下のようになる。
【数１２】
2005522079.tif 000018
ここで、ｗe[m]は、ｍ番目に抽出されたウォーターマークシンボルであり、埋め込まれたウォーターマーク系列に関するＮb回乗算された予測値を含んでいる。式（１２）におけるｗe[m]について解き、αの高次の項を無視すると、以下の近似が得られる。
【数１３】
2005522079.tif 000019
【００５４】
図８に示されるウォーターマーク抽出ステージ２００では、フィルタＨbの出力ｙ’b[n]は、フレームディバイダ２２０の入力に供給され、このフレームディバイダは、オーディオ信号を長さＴSのフレーム、すなわちｙ’b,m[n]に分割し、次いで、エネルギー計算ユニット２３０は、式（１２）により、フレーム分割された信号のそれぞれに対応するエネルギーを計算する。次いで、このエネルギー計算ユニット２３０の出力は、ホワイトニングステージＨW２４０の入力として供給され、このホワイトニングステージは、出力ｗe[m]を供給するように、式（１３）で示される機能を実行する。図９及び図１０には、このホワイトニングステージの代替的な実現２４０Ａ，２４０Ｂが例示されている。
【００５５】
式（１３）の分母は、ホスト（オリジナル）信号ｘの情報を必要とする項を含んでいることが理解される。信号ｘが検出器で利用可能ではないので、ｗe[m]を計算するために、式（１３）の分母が予測される必要があることを意味している。
【００５６】
以下には、かかる予測が２つの説明された窓成形関数（レイズドコサイン窓成形関数及びバイフェーズ窓成形関数）についてどのように達成することができるかが説明されるが、この教示は他の窓成形関数に対しても拡張することできることを理解されたい。
【００５７】
図６ａに示されるレイズドコサイン窓成形関数に関して、ウォーターマークにより誘発されるオーディオエンベロープは、エネルギー関数Ｅ[m]の雑音のある部分にのみ寄与することが実現される。緩やかに変化する部分（すなわち、低周波成分）は、オリジナルのオーディオ信号ｘのエンベローブの寄与のために支配的である。このように、式（１３）は、以下により近似される場合がある。
【数１４】
2005522079.tif 000020
ここで、“ｌｏｗｐａｓｓ（．）”とは、低域通過フィルタ関数である。したがって、関数におけるレイズドコサイン窓成形のためのホワイトニングフィルタＨWは、図９に示されるように実現することができることが理解される。
【００５８】
理解することができるように、かかるホワイトニングフィルタＨW２４０は、信号Ｅ[m]を受信するための入力２４２Ａを備えている。その後、この信号の一部は、低域通過フィルタ２４７Ａを通して通過され、低域通過フィルタリングされたエネルギー信号ＥLP[m]が生成され、次に、このエネルギー信号は、関数Ｅ[m]と共に計算ステージ２４８Ａへの入力として供給される。次いで、計算ステージ２４８Ａは、抽出されたウォーターマークシンボルｗe[m]を計算するために、Ｅ[m]をＥLP[m]で割る。
【００５９】
バイフェーズ窓関数は、エンベッダのウォーターマークコンディショニングステージにおいて利用されるとき、オリジナルのオーディオのエンベロープを予測して、しかるにｗe[m]を計算するために、異なるアプローチが利用されるべきである。
【００６０】
図６ｂに示されるバイフェーズ窓関数を調べることで、オーディオエンベロープがかかる窓関数で変調されるとき、フレームの前半と後半が反対の方向でスケーリングされることが理解される。検出器では、この特性は、ホスト信号ｘのエンベロープエネルギーを予測するために利用される。
【００６１】
結果的に、検出器内で、それぞれのオーディオフレームは、はじめに半分に細分割される。前半のフレーム及び後半のフレームに対応するエネルギー関数は、したがって以下によりそれぞれ与えられる。
【数１５】
2005522079.tif 000021
及び
【数１６】
2005522079.tif 000022
オリジナル信号のエンベロープは、２つのサブフレーム内で反対の方向で変調されるので、オリジナルのオーディオエンベロープは、Ｅ1[m]及びＥ2[m]の平均として近似することができる。
【００６２】
さらに、これら２つの関数の間の差として、瞬間的な変調の値を取ることができる。したがって、バイフェーズ窓関数について、ウォーターマークｗe[m]は以下のように近似することができる。
【数１７】
2005522079.tif 000023
【００６３】
結果として、バイフェーズ窓成形関数について図８におけるホワイトニングフィルタＨW２４０Ｂは、図１０に示されるように実現することができる。入力２４２Ｂ及び２４３Ｂのそれぞれは、前半のフレームＥ1[m]及び後半のフレーム及びＥ2[m]のエネルギー関数を受ける。それぞれのエネルギー関数は、次いで２つに分割され、加算器２４５Ｂ及び２４６Ｂに供給され、この加算器２４５Ｂ及び２４６Ｂは、Ｅ1[m]−Ｅ2[m]、及びＥ1[m]＋Ｅ2[m]をそれぞれ計算する。これら計算された関数の両者は、次いで計算ユニット２４Ｂに通過され、この計算ユニット２４Ｂは、式（１７）に従って埋め込まれたウォーターマーク系列に関してＮb回乗算された予測値を含んでいるｗe[m]を計算するように、加算器２４５Ｂからの値を加算器２４６Ｂからの値で割る。
【００６４】
次いで、この出力ｗe[m]は、バッファリング及び補間ステージ３００に通過され、このステージでは、信号は、デマルチプレクサ（Ｄｅｍｕｘ）３１０により多重化が分解され、長さＬbのバッファ（Ｂｕｆｆｅｒ）３２０で緩衝記憶され、エンベッダと検出器との間の同期の欠陥が解決され、補間ユニット（Ｉｎｔｅｒｐ）３３０内で補間されてエンベッダと検出器との間の時間スケールの変更が補償される。
【００６５】
ウォーターマークの可能性のあるロバスト性を最大にするため、ウォーターマーキングシステムがエンベッダと検出器との間の時間オフセットに対して免疫性があることを確かめることが重要である。言い換えれば、ウォーターマーク検出器は、ホスト信号に挿入されるウォーターマーク系列に同期することが可能である必要がある。
【００６６】
図１１は、オフセット問題を解決するためにバッファリング及び補間ステージ３００により実行される処理を例示している。説明される例は、レイズドコサイン窓成形関数がウォーターマーク埋め込み処理で使用されたとき、オフセットを解決するための処理を例示している。しかし、原理的に、バイフェーズ窓成形関数が使用されたとき、同じ技術を適用することができる。
【００６７】
図１１を参照して、フィルタＨb  ２１０によるフィルタリングの後、到来するオーディオ信号ストリームｙ’b[n]は、フレームディバイダ２２０により、効果的な長さＴSのオーバラッピングフレーム３０２に分割されることが好ましい。
【００６８】
好ましくは、エンベッダと検出器との間の可能性のあるオフセットを解決するため、それぞれのフレームは、Ｎbのサブフレーム３０４ａ，３０４ｂ，３０４ｘに分割され、先の計算（式（１２）から式（１７））がサブフレーム毎に適用される。
【００６９】
好ましくは、それぞれのサブフレームは、隣接するサブフレームとオーバラップする。図示される例では、それぞれのサブフレーム３０４ａ，３０４ｂ，．．．，３０４ｘに関して５０％のオーバラップ（ＴS／Ｎb）が存在し、それぞれのサブフレームは、長さ２ＴS／Ｎbからなる。オーバラップするサブフレームが考慮されるとき、メインフレームは、図１１に示されるようなフレーム間のオーバラップを可能とするために、シンボル周期ＴSよりも長いことが好ましい。
【００７０】
その後、オーディオのエネルギーは、ホワイトニングステージ２４０によりそれぞれのサブフレームについて計算され、結果として得られる値は、デマルチプレクサ３１０によりＮbのバッファ３２０に分解される。バッファ３２０のうちのそれぞれ連続する１つＢ1，Ｂ2，．．．，ＢNbは、値の系列を含んでおり、第一のバッファＢ1は、それぞれのフレーム内に第一のサブフレームに対応する値の系列を含んでおり、第二のバッファＢ2は、それぞれのフレーム内の第二のサブフレームに対応する値の系列を含んでいる等である。
【００７１】
ｗDiがｉ番目のバッファのコンテンツである場合、以下を示すことができる。
【数１８】
2005522079.tif 000024
ここで、Ｌbはバッファ長である。
【００７２】
レイズドコサインウィンドウ成形関数について、埋め込まれたウォーターマークのエネルギーは、フレームの中央の近くに集中され、フレームの中央と最良に配列されるサブフレームが他の全てのサブフレームよりも、明らかに良好な埋め込まれたウォーターマークのシンボルの予測値となる。実際に、それぞれのバッファは、シンボル系列の予測値を含んでおり、系列に対応する予測値は、異なる時間オフセットを有している。
【００７３】
フレームの中央と最良に配列されるサブフレーム（すなわち、正しく配列されたフレームの最良の予測値）は、それぞれのバッファのコンテンツを基準のウォーターマークの系列と相関を取ることで決定される。最大の相関ピーク値をもつ系列は、正しく配列されたフレームの最良の予測値として選択される。対応する信頼レベルは、以下に説明されるように、検出の真理値を決定するために使用される。好ましくは、定義された閾値を超える相関ピークをもつ予測されたウォーターマークの系列がひとたび発見されると、相関プロセスは停止される。
【００７４】
典型的には、それぞれのバッファの長さは、ウォーターマーク系列の長さＬWの３倍と４倍との間であり、典型的には２０４８シンボルと８１９２シンボルとの間の長さからなり、Ｎbは、典型的に２から８の範囲内にある。
【００７５】
バッファは、それぞれのウォーターマークシンボルが該シンボルの幾つかの予測値の平均を取ることで構築することができるように、通常、ウォーターマーク系列の３倍から４倍の長さである。この平均処理は平滑化と呼ばれ、平均が行われる回数は、スムージングファクタｓfと呼ばれる。このように、バッファの長さＬb及びウォーターマーク系列の長さＬWが与えられると、スムージングファクタｓfは、以下のようになる。
【００７６】
【数１９】
2005522079.tif 000025
【００７７】
別の好適な実施の形態では、検出器は、前のサーチステップの結果に基づいてオフセットのサーチで使用されるパラメータをリファインする。たとえば、バッファＢ3に記憶されている結果が情報信号の最良の予測値を提供することを第一の予測値の系列が示す場合、（同じ受信された信号、又は次の検出窓の間に受信された信号のいずれかに関する）次のオフセットのサーチは、最良の予測値のサブフレームの位置に向けてサブフレームの位置をシフトすることでリファインされる。このようにして、ゼロのオフセットを有する系列の予測値は、繰返し改善することができる。
【００７８】
図８に示されるように、バッファリングステージからの出力ｗD1，ｗD2，．．．，ｗDNbは、補間ステージに通過され、補間の後、このステージの出力ｗI1，ｗI2，．．．，ｗINbは、ウォーターマーク信号における可能性のある時間スケールの変更を解決するために必要とされるものであり、相関及び判定ステージに通過される。異なる可能性のあるオフセット値に対応するウォーターマークの全ての予測値ｗI1，ｗI2，．．．，ｗINbは、相関及び判定ステージ４００に通過される。
【００７９】
相関器４１０は、基準となるウォーターマークの系列ｗC[k]に関してそれぞれの予測値ｗIj，ｊ＝１，．．．，Ｎbの相関を計算する。次いで、それぞれの予測値に対応するそれぞれの相関出力は、最大値の検出ユニット４２０に印加され、この検出ユニットは、その２つの予測値が最大の相関ピーク値を与えるかを判定する。これらの予測値は、基準となるウォーターマークの巡回的にシフトされたバージョンｗd1及びｗd2に最も良くフィットする予測値として選択される。これら予測された系列の相関値は、閾値検出器及びペイロード抽出ユニット４３０に通過される。
【００８０】
検出器で使用される基準のウォーターマーク系列ｗSは、ホスト信号に印加される（おそらく巡回的にシフトされたバージョンの）オリジナルのウォーターマーク系列に対応する。たとえば、エンベデッダ内でシードＳをもつランダム数発生器を使用してウォーターマーク信号が計算された場合、同様に、検出器は、ウォーターマーク信号を決定するように、同じランダム数発生アルゴリズム及び同じ初期シードＳを使用して同じランダム数の系列を計算することができる。代替的に、基準は単にいずれか予め決定された系列であるので、ウォーターマーク信号はエンベッダに最初に供給され、検出器により利用される。
【００８１】
図１２は、典型的な相関関数の形状を相関器４１０からの出力として示している。水平のスケールは、（系列のサンプルの観点での）相関遅延を示している。左側の垂直スケール（信頼レベルｃＬとも呼ばれる）は、正規分布された相関関数の標準偏差に関して正規化された相関ピークの値を表している。
【００８２】
理解できるように、典型的な相関は、ｃＬに関して比較的平坦であり、ｃＬ＝０に関して中央に置かれる。しかし、関数は、２つのピークを含んでおり、これら２つのピークは、ｐＬだけ離れており（式（６）参照）、ウォーターマークが存在するときに閾値を超えるｃＬ値の方向に延びている。相関ピークが負であるとき、先の説明はそれらの絶対値にも当てはまる。
【００８３】
（図ではｃＬ＝８．７として設定されて示される）水平ラインは、検出閾値を表している。検出閾値は、誤報率を表している。検出閾値は、誤報率を制御する。
【００８４】
２種類の誤報が存在する。偽陽性率は、ウォーターマークされていないアイテムでウォーターマークを検出する確率として定義され、偽陰性率は、ウォーターマークされたアイテムでウォーターマークを検出しない確率として定義される。一般に、偽陽性率の要件は、偽陰性率の要件よりも厳格である。図１１の右側のスケールは、偽陽性率ｂの確率を説明している。図示される例に見ることができるように、偽陽性率ｂ＝１０-12の確率は、閾値ｃＬ＝８．７に等価であり、ｂ＝１０-87はｃＬ＝２０に等価である。
【００８５】
それぞれの検出のインターバルの後、検出器は、オリジナルのウォーターマークが存在するか、存在しないかを判定し、これに基づいて“ｙｅｓ”又は“ｎｏ”の判定を出力する。望まれる場合、この判定実行処理を改善するため、多数の判定窓が考慮される場合がある。かかる例では、偽陽性率は、所望の基準に依存して、考慮されるそれぞれの検出窓に個々の確率の組み合わせである。たとえば、相関関数が３つの検出インターバルのうちのいずれかの２つに関して閾値ｃＬ＝７を超える２つのピークを有する場合、ウォーターマークが存在するものと思われることが判定される。かかる検出基準は、ウォーターマーク信号の所望の使用に依存して、並びにオリジナルのホスト信号の品質のようなファクタ及び通常の伝送の間に信号がどの位粗悪に破壊される可能性があるかを考慮して変更することができる。
【００８６】
ペイロード抽出ユニット４３０は、検出されたウォーターマーク信号からペイロード（たとえば、情報コンテンツ）を抽出するために続いて利用される場合がある。ユニットが検出の閾値を超える２つの相関ピークｃＬ1及びｃＬ2をひとたび予測すると、（式（６）で定義される）巡回シフトｃＬの予測値ｃＬ’は、ピーク間の距離として導出される。つぎに、相関ピークの符号ρ1及びρ2が決定され、したがってｒSIGNが式（７）から計算される。次いで、全体のウォーターマークのペイロードは、式（８）を使用して計算される場合がある。
【００８７】
たとえば、図１２において、ｐＬが２つのピーク間の相対的な距離であることを理解することができる。両方のピークは正、すなわちρ1＝＋１及びρ2＝＋１である。式（７）から、ｒSIGN＝３である。結果的に、ペイロードｐＬW＝＜３，ｐＬ＞である。
【００８８】
図８に記載されるシンボル抽出及びバッファリングステージは、図１３に示される装置５００により効率的に実現することができる。ここで、オフセットの補償は、余分の計算なしに達成される。また、多重化の分解も遅延及び間引きブロックのシンプルなセットで実現されることが理解される。
【００８９】
第一に、到来するフレーム信号ｙb,mは、長さＴS／ＮbのＮbのオーバラップしないサブフレームに細分割され、それぞれのサブフレームのエネルギーは、エネルギー計算ユニット２３０を使用して計算される。第二に、ホワイトニングフィルタＨWは、ホワイトニングユニット２４０に適用される。遅延ユニット５１０と加算ユニット５２０の組み合わせは、隣接するサブフレームの間で５０％のオーバラップを効果的に実現する。ウォーターマークシンボル系列ｗe[m]が加算ユニット５２０の出力で生成された後、この系列は、遅延セット５１２とダウンサンプリングセット５３０の組み合わせを使用してＮbのバッファユニット３２０にわたり続いて分散される。
【００９０】
これは、それぞれのバッファがｗe[k]の到来するサンプルにおけるそれぞれのＮbについて１つの値を得るように行われる。たとえば、全てのバッファが満たされるまで、第一のサンプルがｗD1に進み、第二のサンプルがｗD2に進み、第三のサンプルがｗD3に進み、．．．第ＮbのサンプルがｗDNbに進む場合、第（Ｎb＋１）のサンプルは、ｗD1に戻る等となる。このように、ｉ番目のバッファ入力ｗDi[k]は、以下のように表される。
【数２０】
2005522079.tif 000026
【００９１】
ｗDi[k]のサンプリング周波数は、ｗe[m]のサンプリング周波数の１／Ｎbであることが理解される。この間引きは、図１３における間引きセット５３２を介して達成される。
【００９２】
エネルギー計算ユニット２３０ではオーバラップしていないフレームが考慮されるので、Ｎbの系列を生成するために必要とされる全体の計算量は、全体のフレームにわたり延びるシンボルをもつ唯一の系列が計算される場合に必要とされた計算量と同じである。
【００９３】
当業者であれば、特に説明されていない様々な実現は、本発明の範囲に含まれるものとして理解されたい。たとえば、検出装置の機能のみが説明されたが、本装置は、デジタル回路、アナログ回路、コンピュータプログラム、又はそれらの組み合わせとして実現されることを理解されたい。
【００９４】
同様に、先の実施の形態は、オーディオ信号を参照しながら説明されたが、本発明は、たとえば、ビデオ及びデータ信号のような情報又はマルチメディア信号といった情報を他のタイプの信号に加算することに適用することができる。
【００９５】
さらに、本発明は、１つのみのウォーターマーキング系列を含むウォーターマーキングスキーム（すなわち、１ビットスキーム）、又は多数のウォーターマーキング系列を含むウォーターマーキングスキームに適用することができることが理解される。かかる多数の系列は、ホスト信号内に同時又は連続的に埋め込むことができる。
【００９６】
同様にして、先のウォーターマークの検出は、それぞれの予測値は相関がとられるものとして説明されたが、相関の手順は、ウォーターマークの正の検出がひとたび行われると停止するように調整することができることが理解される。さらに、デコーダは、バッファの再順序付け（バッファが相関付けされる順序）により、時間オフセットを適応的に補償するように調整することができ、これにより、現在の検出窓において最良に配列されるバッファが、次の検出窓で相関がとられるべき最初のバッファとなる。
【００９７】
本明細書では、用語“comprising”は、他のエレメント又はステップを排除するものではなく、“a”又は“and”は、複数であることを排除するものではなく、単一のプロセッサ又は他のユニットが請求項で列挙される幾つかの手段の機能を実現する場合があることが理解される。
【図面の簡単な説明】
【００９８】
【図１】ウォーターマークの埋め込み装置を説明する図である。
【図２】信号部分抽出フィルタＨを示す図である。
【図３】図３ａは、図２に示されるフィルタＨの周波数の関数としての典型的な振幅応答を示すグラフ、図３ｂは、図２に示されるフィルタＨの周波数の関数としての典型的な位相応答を示すグラグである。
【図４】図１に示される装置のペイロードエンベデッド及びウォーターマークコンディショニングステージを示す図である。
【図５】それぞれのステージで関連する信号のチャートを含む、図４のウォーターマークコンディショニング装置ＨCの詳細を説明する図である。
【図６】図６ａは、レイズドコサイン関数の形式での好適な代替の窓成形関数を示すグラグ、図６ｂは、バイフェーズ関数の形式での好適な代替の窓成形関数を示すグラフである。
【図７】図７ａは、レイズドコサイン窓成形関数によりコンディショニングされたウォーターマーク系列の周波数スペクトルを示すグラフ、図７ｂは、バイフェーズ窓成形関数によりコンディショニングされたウォーターマーク系列の周波数スペクトルを示すグラフである。
【図８】本発明の実施の形態に係るウォーターマーク検出器を説明する図である。
【図９】レイズドコサイン成形窓関数との使用のための、図８のホワイトニングフィルタＨWを示す図である。
【図１０】バイフェーズ窓成形関数との使用のための、図８のホワイトニングフィルタＨWを示す図である。
【図１１】本発明の実施の形態に係るウォーターマークシンボル抽出及びバッファリング処理の詳細を示す図である。
【図１２】図８に示されるウォーターマーク検出器の相関器からの相関関数出力の典型的な形状を示すグラフである。
【図１３】シンボル抽出及びバッファリングステージの１つの好適な実現の例を示す図である。
【図１】
2005522079.tif 000027
【図２】
2005522079.tif 000028
【図３ａ】
2005522079.tif 000029
【図３ｂ】
2005522079.tif 000030
【図４】
2005522079.tif 000031
【図５】
2005522079.tif 000032
【図６】
2005522079.tif 000033
【図７】
2005522079.tif 000034
【図８】
2005522079.tif 000035
【図９】
2005522079.tif 000036
【図１０】
2005522079.tif 000037
【図１１】
2005522079.tif 000038
【図１２】
2005522079.tif 000039
【図１３】
2005522079.tif 000040
【国際調査報告】
2005522079.tif 000041
2005522079.tif 000042
2005522079.tif 000043