(57)【要約】エコー、残響および／またはノイズを引き起こす音響経路を通じた透かし入れされたオーディオ信号の伝送後は、透かし検出相関結果ピーク値は、主たる相関結果ピークのまわりの限られた時間的範囲内に集中し、該限られた時間的範囲は全相関長Nよりずっと小さい。透かし検出は下記に基づく。マークされていないコンテンツについての相関結果値についての平均的なまたは期待される確率分布においてnp個の相関結果ピーク値v＝(v1,v2,…,vnp)が与えられたとき、前記限られた時間的範囲をカバーする窓内で、これらのピークvより大きいかそれに等しい、Nの相関結果値集合からのnp個以上の値がある確率が計算される。
2016515721.tif 000002
【特許請求の範囲】
【請求項１】
エコー、残響および／またはノイズを含みうる受領オーディオ信号中の透かしシンボルを判別する方法であって、
・前記受領オーディオ信号のセクションの、対応する異なる候補透かしシンボルisに関係した少なくとも二つの異なる候補参照パターン信号との相関を調べ、各候補透かしシンボルisについてそれぞれの場合にN個の相関結果値
【数６６】
2016515721.tif 000068
を生成するステップと、
・各候補透かしシンボルisについて、前記相関結果値cis内のM個のピーク値を決定するステップとを含んでおり、
・グローバル・ピーク値と記される前記M個のピーク値のそれぞれについておよび各候補透かしシンボルisについて、前記相関結果値cisから、長さLのM個の窓のそれぞれの中のnp個のピーク値の対応するベクトル
【数６７】
2016515721.tif 000069
を決定するステップであって、各長さLの窓は前記M個のグローバル・ピーク値のうちの一つを含み、長さLの窓全部の合計長さはNより小さく、k＝1,…,Mである、ステップと、
・長さLの前記M個の窓のそれぞれについて、ピーク値の前記ベクトルviskおよび透かし入れされていないオーディオ信号についての相関結果値についての平均のまたは期待される確率分布に基づいて、偽陽性確率値
【数６８】
2016515721.tif 000070
を計算するステップであって、前記偽陽性確率FP(visk)は、一回または複数回、そのような長さLの窓が、透かし入れされていないオーディオ信号コンテンツについての相関値についての平均または期待される確率分布において前記np個のピーク値より大きいまたはそれに等しいnp個以上の相関結果絶対値を含むを含む確率として定義される、ステップと、
・前記受領オーディオ信号の前記セクションについて、最小の偽陽性確率値minis(mink(FP(visk)))をもつ候補透かしシンボルisを、該最小の偽陽性確率値minis(mink(FP(visk)))が所定の閾値（Tmax）より小さい場合には、検出された透かしシンボルとして選択し、そうでない場合には透かしシンボルが検出されなかったと想定するステップとを特徴とする、
方法。
【請求項２】
前記偽陽性確率FP(visk)値の前記計算において、その値がさらなる所定の閾値（Tmin）より小さい場合には、正しい候補透かしシンボルisが見つかったと想定され、前記のさらなる処理が省略され、前記所定の閾値（Tmax）は前記さらなる所定の閾値（Tmin）より大きい、請求項１記載の方法。
【請求項３】
M＝3である、請求項１または２記載の方法。
【請求項４】
二つの異なる候補透かしシンボルがある、請求項１ないし３のうちいずれか一項記載の方法。
【請求項５】
エコー、残響および／またはノイズを含みうる受領オーディオ信号中の透かしシンボルを判別する装置であって、
・前記受領オーディオ信号のセクションの、対応する異なる候補透かしシンボルisに関係した少なくとも二つの異なる候補参照パターン信号との相関を調べ、各候補透かしシンボルisについてそれぞれの場合にN個の相関結果値
【数６９】
2016515721.tif 000071
を生成するよう適応された手段と、
・各候補透かしシンボルisについて、前記相関結果値cis内のM個のピーク値を決定し、各候補透かしシンボルisについての前記M個の相関結果ピーク値から導出されて、前記受領オーディオ信号の前記セクションについて、前記候補透かしシンボルisのうちの透かしシンボルを決定するよう適応された手段とを含んでおり、
グローバル・ピーク値と記される前記M個のピーク値のそれぞれについておよび各候補透かしシンボルisについて、前記相関結果値cisから、長さLのM個の窓のそれぞれの中のnp個のピーク値の対応するベクトル
【数７０】
2016515721.tif 000072
が決定され、各長さLの窓は前記M個のグローバル・ピーク値のうちの一つを含み、長さLの窓全部の合計長さはNより小さく、k＝1,…,Mであり、
長さLの前記M個の窓のそれぞれについて、ピーク値の前記ベクトルviskおよび透かし入れされていないオーディオ信号についての相関結果値についての平均のまたは期待される確率分布に基づいて、偽陽性確率値
【数７１】
2016515721.tif 000073
が計算され、ここで、前記偽陽性確率FP(visk)は、一回または複数回、そのような長さLの窓が、透かし入れされていないオーディオ信号コンテンツについての相関値についての平均または期待される確率分布において前記np個のピーク値より大きいまたはそれに等しいnp個以上の相関結果絶対値を含むを含む確率として定義され、
前記受領オーディオ信号の前記セクションについて、最小の偽陽性確率値minis(mink(FP(visk)))をもつ候補透かしシンボルisが、該最小の偽陽性確率値minis(mink(FP(visk)))が所定の閾値（Tmax）より小さい場合には、検出された透かしシンボルとして選択され、
そうでない場合には透かしシンボルが検出されなかったと想定される、
装置。
【請求項６】
前記偽陽性確率FP(visk)値の前記計算において、その値がさらなる所定の閾値（Tmin）より小さい場合には、正しい候補透かしシンボルisが見つかったと想定され、前記のさらなる処理が省略され、前記所定の閾値（Tmax）は前記さらなる所定の閾値（Tmin）より大きい、請求項５記載の装置。
【請求項７】
M＝3である、請求項５または６記載の装置。
【請求項８】
二つの異なる候補透かしシンボルがある、請求項５ないし７のうちいずれか一項記載の装置。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、エコー、残響および／またはノイズを含みうる受領オーディオ信号中の透かしシンボルを判別する方法および装置に関する。
【背景技術】
【０００２】
オーディオ透かし入れは、耳に聞こえない仕方で追加的情報をオーディオ信号中に埋め込むプロセスである。埋め込みは、たとえば擬似ランダムなノイズまたはエコーを加えることによってオーディオ信号を変化させることによって実行される。埋め込みを耳に聞こえないようにするために、埋め込みの強度は、オーディオ信号の音響心理学的解析によって制御される。特許文献１は、オーディオ信号、たとえばマイクロフォンによって受領されたラウドスピーカー音におけるエコー、残響および／またはノイズの存在時の透かし検出を記載している。これらのエコーは、透かしシンボル（すなわち参照信号）との長さNの相関結果値シーケンス内の複数のピークを与え、透かし検出の信頼性を改善するために使われる。この統計的検出器の基本的段階は次のようなものである。
・各候補透かしシンボルについて前記相関結果シーケンスにおけるピーク値v＝(v1,v2,…,vnp)を見出す。ここで、v1≧…≧vnpであり、npは考慮に入れられる相関結果ピーク値の数である。
・前記候補透かしシンボルが埋め込まれている偽陽性確率P(k)、k＝1,…,npを再帰的に計算する。
・最低のP(k)値を与える候補透かしシンボルを選択する。
【０００３】
P(k)は、候補透かしシンボルを誤って受け入れる確率である。これは、透かしが入っていない信号セクションからのk個以上の相関結果値が考慮対象の実際のk個のピーク値以上である確率を記述する。
【０００４】
この統計的検出器は次の問題を解決する。
ａ）いかにして確率P(k)を再帰的に評価するか。ここで、考慮されるピーク値の数kが徐々に増大させられることができる。
ｂ）すでに実行された計算を再利用することによりいかにして計算量を最小化するか。
【０００５】
この統計的検出器は、検出パフォーマンスを改善するためにいくつかの相関結果ピークを使う。特に、この改善は、透かし入れされたトラックが、エコーに起因するマルチパス検出につながる音響経路を通じて伝送される場合に有利である。P(k)を計算するために使われるnp個のピークv1≧…≧vnpは、長さNの相関結果値シーケンス全体から取られる。
【先行技術文献】
【特許文献】
【０００６】
【特許文献１】国際公開第2011/141292号
【発明の概要】
【発明が解決しようとする課題】
【０００７】
しかしながら、この種の処理は、物理的な現実を最適な仕方で反映するものではない。相関結果値シーケンスにおける追加的なピーク値が残響から生じ、よって主ピークのまわりに近接してグループ化されるからである。すなわち、これらの追加的なピーク値は、ラウドスピーカーと検出装置のマイクロフォンとの間の直接経路に対して時間遅延される。ただし、限られた時間期間内のみである。透かし入れ検出装置が源からの距離ds（サンプルs単位で測る）をもつ音響経路を通じて音を受領するとき、伝搬距離はdm＝dscTメートルである。ここで、cTは一サンプリング区間内に音が伝搬する距離であり、T＝1/fはサンプリング・レートであり、cは音速である。
【０００８】
たとえば、直接（direct）経路についてdmd＝4mであり、間接（indirect）経路についてdmi＝2dmd＝8mである場合、直接経路の主ピークと反射に起因する近隣ピークとの間のサンプル単位での距離は次のようになる。
【０００９】
【数１】
2016515721.tif 000003
典型的にはセッティングでは、相関長はN＝16kサンプルである（1k≡1024）。したがって、ピークは、主ピーク値または複数のピーク値のまわりの、たとえば
【００１０】
【数２】
2016515721.tif 000004
サンプルのサイズ（たとえばL≪N）の窓において探索され、対応する偽陽性確率が計算されることになる。主相関結果ピーク値はN個の相関結果値の現在の集合におけるどこかに位置するので、偽陽性確率は、サイズNのバッファ内のサイズLの窓のあらゆる可能なN−L＋1通りのシフトについて計算される。
【００１１】
本発明によって解決されるべき課題は、既知の統計的透かし検出器よりよく、マルチパス受信の物理を考慮に入れ、それにより偽陽性確率計算および透かし検出のパフォーマンスを改善することである。
【課題を解決するための手段】
【００１２】
この課題は、請求項１に開示される方法によって解決される。この方法を利用する装置は請求項２において開示される。
【００１３】
上述したように、エコー、残響および／またはノイズを引き起こす音響経路を通じた透かし入れされたオーディオ信号の伝送後は、透かし検出相関結果ピーク値は一般に、限られた時間的範囲内で主たる相関結果ピークのまわりに集中し、Lと表わす該時間的範囲の最大サイズは全相関長Nよりずっと小さい。透かし検出のための関係したタスクは、次のように定式化できる。
【００１４】
マークされていないオーディオ信号コンテンツについて相関結果値についての平均的なまたは期待される確率分布において、np個のピーク値v＝(v1,v2,…,vnp)が与えられたとき、これらのピークvより大きいかそれに等しい、現在の長さNの相関結果値集合からの長さLのスライディング窓内のnp個以上のピーク値がある確率PFP(np,L,v)をいかにして計算するか。
【００１５】
測定により、マークされていないコンテンツについての相関結果値についてのそのような平均または期待される確率分布がガウス分布に対応するまたはガウス分布と同様であることが示されている。
【００１６】
コメント：一サンプルずつ参照パターンを現在の入力オーディオ信号セクションを通じてシフトさせていくことにより相関を実行するとき、この相関のN値結果は正のピーク値のほか負のピーク値をもつことがあり、これらは一緒にして「ピーク量値」と表わされる。
【００１７】
以下の記述を簡単にするために、いくつかの記法を導入しておく。
【００１８】
【数３】
2016515721.tif 000005
は、N内の位置jで始まる長さLのj番目のスライディング窓内のL個の相関結果値を表わす。
【００１９】
【数４】
2016515721.tif 000006
は、ベクトルcjにおいて、マークされていないコンテンツについての相関結果値についての平均または期待される確率分布において前記np個のピーク値vより大きいまたはそれに等しいnp個以上の現在の相関結果値があることを示す。
【００２０】
【数５】
2016515721.tif 000007
は、ベクトルcjにおいて、マークされていないコンテンツについての相関結果値についての平均または期待される確率分布において前記np個のピーク値vより大きいまたはそれに等しい現在の相関結果値がnp個未満であることを意味する。
【００２１】
コメント：Prが確率を表わすとして、
【００２２】
【数６】
2016515721.tif 000008
という点で、
【００２３】
【数７】
2016515721.tif 000009
は
【００２４】
【数８】
2016515721.tif 000010
と相補的である。
【００２５】
長さLのスライディング窓をN個の相関結果値を通じてシフトさせる。上述したように、求められるFP（false  positive［偽陽性］）確率は、一回または複数回、スライディング窓が、透かし入れされていないオーディオ信号コンテンツについての相関結果値についての平均または期待される確率分布において前記np個のピーク値より大きいまたはそれに等しいnp個以上の相関結果量値を含む確率である。このFP確率についての相補的な場合は、透かし入れされていないオーディオ信号コンテンツについての相関結果値についての平均または期待される確率分布においてnp個のピーク値より大きいまたはそれに等しいnp個以上の相関結果量値を含むスライディング窓がない、つまり
【００２６】
【数９】
2016515721.tif 000011
というものである。よって、
【００２７】
【数１０】
2016515721.tif 000012
は偽陽性確率
【００２８】
【数１１】
2016515721.tif 000013
（一つまたは複数の窓がjによってインデックス付けされる）についての相補的確率である。
【００２９】
最終的なFP確率計算は次のように表現できる（詳細は後述）。
【００３０】
【数１２】
2016515721.tif 000014
L＝Nについては（この場合、
【００３１】
【数１３】
2016515721.tif 000015
である）、
【００３２】
【数１４】
2016515721.tif 000016
の定義における上記のコメントに鑑み、この一般的な公式は、特許文献１の統計的検出器における再帰的計算の場合
【００３３】
【数１５】
2016515721.tif 000017
に帰着する。
【００３４】
【数１６】
2016515721.tif 000018
の計算を、特許文献１の統計的検出器について記述しておく。現在の信号セクションについての相関結果値におけるピーク値に基づいて、候補シンボルのどれが現在の信号セクションに存在するかが、透かしシンボルのその種類の検出の偽陽性確率の関係した値を使って、検出される。ここで、偽陽性確率は再帰的な仕方で計算され、所与の数の相関結果ピーク値についての全偽陽性確率が、初期にはより少数の相関結果ピーク値についての偽陽性確率を使い、要求される検出信頼度に従って徐々に考慮される相関結果ピーク値の数を増していくことによって、評価される。
【００３５】
したがって、以下の記述における偽陽性確率の計算は
【００３６】
【数１７】
2016515721.tif 000019
を計算する問題にマッピングされる。所与の数npのピーク値について、この確率は、ピークi＝1,2で始めてnp−1個のピーク値についての確率から再帰的に計算されることができる（「発明を実施するための形態」の節参照）。
【００３７】
本発明は、偽陽性検出判断を有意に低減することによって特許文献１の統計的透かし検出処理の検出パフォーマンスを改善するが、本発明は特許文献１の透かし検出のすべての利点を保持する。
【００３８】
有利なことに、本発明は、ユーザーがテレビを視聴し、透かし入れされたテレビ音を介して透かし情報を得るというセカンドスクリーンのシナリオにおいて使用されることができる。ここで、前記透かし情報はたとえば、該透かし情報を使って現在のテレビ番組に関係した追加的情報をダウンロードして基本的に同期して呈示するタブレット・コンピュータのためである。
【００３９】
原理的には、本発明の方法は、エコー、残響および／またはノイズを含みうる受領オーディオ信号中の透かしシンボルを判別するために好適であり、本方法は、
・受領されたオーディオ信号の現在セクションの、対応する異なる候補透かしシンボルisに関係した少なくとも二つの異なる候補参照パターン信号との相関を調べ、各候補透かしシンボルisについてそれぞれの場合にN個の現在の相関結果値cisを生成するステップと、
・各候補透かしシンボルisについて、前記相関結果値cis内のM個のピーク値を決定するステップと、
・各候補透かしシンボルisについての前記M個の相関結果ピーク値から導出されて、前記受領オーディオ信号の前記現在セクションについて、前記候補透かしシンボルisのうちの透かしシンボルを決定するまたは透かしシンボルが検出されなかったと判定するステップとを含み、前記透かしシンボルの決定はさらに、
・グローバル・ピーク値と記される前記M個のピーク値のそれぞれについておよび各候補透かしシンボルisについて、現在の相関結果値cisから、長さLのM個の窓のそれぞれの中のnp個のピーク値の対応するベクトル
【００４０】
【数１８】
2016515721.tif 000020
を決定するステップであって、各長さLの窓は前記M個のグローバル・ピーク値のうちの一つを含み、長さLの窓全部の合計長さはNより小さく、k＝1,…,Mである、ステップと、
・長さLの前記M個の窓のそれぞれについて、現在の相関結果値cisから、一回または複数回、そのような長さLの窓が、透かし入れされていないオーディオ信号についての相関結果値についての平均または期待される確率分布においてnp個のピーク値より大きいまたはそれに等しいnp個以上の相関結果値を含むを含む確率に対応する偽陽性確率
【００４１】
【数１９】
2016515721.tif 000021
の値を計算するステップと、
・前記受領オーディオ信号の現在セクションについて、最小の偽陽性確率値
【００４２】
【数２０】
2016515721.tif 000022
をもつ候補透かしシンボルisを、検出された透かしシンボルとして選択するステップとを含む。
【００４３】
原理的には、本発明の装置は、エコー、残響および／またはノイズを含みうる受領オーディオ信号中の透かしシンボルを判別するために好適であり、本装置は、
・受領されたオーディオ信号の現在セクションの、対応する異なる候補透かしシンボルisに関係した少なくとも二つの異なる候補参照パターン信号との相関を調べ、各候補透かしシンボルisについてそれぞれの場合にN個の現在の相関結果値cisを生成するよう適応された手段と、
・各候補透かしシンボルisについて、前記相関結果値cis内のM個のピーク値を決定し、各候補透かしシンボルisについての前記M個の相関結果ピーク値から導出されて、前記受領オーディオ信号の前記現在セクションについて、前記候補透かしシンボルisのうちの透かしシンボルを決定するまたは透かしシンボルが検出されなかったと判定するよう適応された手段とを含み、
グローバル・ピーク値と記される前記M個のピーク値のそれぞれについておよび各候補透かしシンボルisについて、現在の相関結果値cisから、長さLのM個の窓のそれぞれの中のnp個のピーク値の対応するベクトル
【００４４】
【数２１】
2016515721.tif 000023
が決定され、各長さLの窓は前記M個のグローバル・ピーク値のうちの一つを含み、長さLの窓全部の合計長さはNより小さく、k＝1,…,Mであり、
長さLの前記M個の窓のそれぞれについて、現在の相関結果値cisから、一回または複数回、そのような長さLの窓が、透かし入れされていないオーディオ信号についての相関結果値についての平均または期待される確率分布においてnp個のピーク値より大きいまたはそれに等しいnp個以上の相関結果値を含むを含む確率に対応する偽陽性確率
【００４５】
【数２２】
2016515721.tif 000024
の値が計算され、
前記受領オーディオ信号の現在セクションについて、最小の偽陽性確率値
【００４６】
【数２３】
2016515721.tif 000025
をもつ候補透かしシンボルisが検出された透かしシンボルとして選択される。
【００４７】
本発明の有利な追加的実施形態はそれぞれの従属請求項において開示される。
【図面の簡単な説明】
【００４８】
本発明の例示的な実施形態が付属の図面を参照して記述される。
【図１】三つの相関結果ピーク値v3,v2,v1より大きいまたはそれに等しい三つ以上の相関値の分離した〔重なりのない〕可能性を示す図である。
【図２】c2Lが≧vである二つの値をもつ分離した可能性を、対応する代表ベクトル[0,1,1],[0,0,2],[m≧2,0],[m≧1,1,0],[m≧1,0,1]および零補償された代表ベクトル[[0,1,1]],[[0,1,1]],[[1,m≧1,0]],[m≧1,0],[m≧1,0,1]とともに示す図である。
【図３】「m≧」の場合についての一般的な値分布を示す図である。
【図４】本発明の処理についてのブロック図である。
【図５】本発明の処理についての流れ図である。
【図６】三つの長さLの窓が描かれている例示的な相関結果を示す図である。
【発明を実施するための形態】
【００４９】
〈Ａ．偽陽性（FP）確率の定義〉
相互相関に基づく透かし検出において、埋め込まれた透かし情報を判別するために単一の相関結果ピーク値を使うことはよく知られている。しかしながら、本発明では、複数の相関結果ピーク量値を用いる最適な透かし検出が記述される。FP確率を評価するために、次のことが想定される。
・任意のマークされていない信号セクション・コンテンツと参照パターンまたは透かし〔ウォーターマーク〕シンボルとの間の相関値が独立同分布のランダム変数であること、および
・特定の相関結果ピーク値vについて、対応する確率分布尾部確率
【００５０】
【数２４】
2016515721.tif 000026
をいかにして計算するかは知られていること。ここで、cはガウシアンなど特定の分布に従う相関値を表わす。
【００５１】
長さNの現在の相関結果における単一のピーク値vについては、FP確率は、マークされていないコンテンツについての相関結果値についての平均または期待される確率分布において、N個の相関値のうちそのピーク値vより小さくない一つまたは複数の値がある確率である。同様に、np個のピーク値v1≧v2≧…≧vnpについては、偽陽性確率は、マークされていないコンテンツについての相関結果値についての平均または期待される確率分布において、これらのピークv1,v2,…,vnpより大きいまたはそれに等しいnp個以上の相関結果値がある確率として定義される。ここで、対応する尾部確率は
【００５２】
【数２５】
2016515721.tif 000027
として定義される。
【００５３】
ソート後の特定の相関結果ベクトル・サンプル
【００５４】
【数２６】
2016515721.tif 000028
（下付のn1および上付きのn2を付けたxはn1からn2のインデックスを付けられたシーケンス値xを集めたベクトルを表わす）、すなわち、c1≧c2≧…≧cLが与えられると、それが
【００５５】
【数２７】
2016515721.tif 000029
より大きいまたはそれに等しいnp個以上の値をもつかどうかを決定できる。
【００５６】
nGT＝0;
for  i＝1:np
if  ci≧vi
nGT＝nGT＋1;
end
end
。
【００５７】
nGT＝npであれば、c1Lはv1npより大きいまたはそれに等しいnp個以上の値をもつ。そのことは、以下では
【００５８】
【数２８】
2016515721.tif 000030
〔またはc1L≧v1np〕と記される。c1L≧v1npとなる場合は分離した〔重なりのない〕イベントに分割される。最終的なFP確率は、個々に生起するイベントの確率の和である。
【００５９】
〈Ａ１．相関値分布、複数ピークとの比較〉
相関結果値を複数のピーク値と比較するために、相関結果量値の完全な範囲がnp＋1個の区間
[−∞,vnp)，[vnp,vnp-1)，[vnp-1,vnp-2)，…,[v1,＋∞)
に分割される。
【００６０】
次いで、相関値分布は、個々の区間内に何個の相関結果量値が位置しているかを数えることによって実行される。それは代表ベクトルによって記述できる。時によっていくつかの区間内のいくつかの値がFP確率評価にとっては重要でないので、代表ベクトルは異なる値をもつことがある。したがって、本稿では、代表ベクトルの最も右の要素が常に区間[v1,＋∞)に対応し、一方、最も左の要素は第一の要素と称される。
【００６１】
図１は、np＝3個以上の相関結果値が三つのピーク値v3≦v2≦v1より大きいまたはそれに等しくなる種々の場合についての相関値分布を示している。場合P1は、代表ベクトル[m≧3]によって記述できる。これは、区間[v1,＋∞)内に3個以上の相関結果値があり、残りのL−m個の相関結果値はピーク値v1未満であることを示す。一方では、区間[v3,v2)や[v2,v1)に何個の値があるかに関わりなく、これらの場合はみな、v1未満のL−m個の相関結果値があるという条件に含まれる。他方では、三つのピークとの比較のために、個々の区間[v3,v2)，[v2,v1)，[v1,＋∞)内に何個の相関結果値があるかを検査することが便利である。代表ベクトル[m≧3]については、最悪ケースは、区間[v3,v2)にも区間[v2,v1)にも相関結果値がないというものであり、代表ベクトル[0,0,m≧3]に対応する。
【００６２】
そのような場合、v1≧v2、v1≧v3のため、それでも、少なくとも、≧v1である一つの値、≧v2である一つの値および≧v3である一つの値がある。これは、区間[v3,v2)および[v2,v1)に関連付けられた二つの0が区間[v1,＋∞)におけるm≧3によって補償されていると解釈されることができる。簡単のため、[[…]]と記される「零補償された代表ベクトル」がこの場合、もとの代表ベクトル[m≧3]から導出されて、[[1,1,m≧1]]として定義できる。特に断わりのない限り、長さnpをもつ零補償された代表ベクトルが区間[vnp,＋∞)について使われる。
【００６３】
たとえば、場合P2では、区間[v1,＋∞)内に二つの相関結果値、区間[v2,v1)内に0個の相関結果値、区間[v3,v2)内にm≧1個の相関結果値があり、L−2−m個の相関結果値はv3未満である。
【００６４】
代表ベクトルa＝[an,…,a2,a1]を与えられたとき、その対応する零補償されたものは次のように導出される。
・np−n個の0を加えることによってベクトルaを長さnpに拡張する
【００６５】
【数２９】
2016515721.tif 000031
・{aj(t),1≦t≦T}として、a'内の「1」より大きな要素を見出す。
・t＝1,2,…,Tについて、零補償を実行する。aj(t)の左側に0があれば、右から左に、最大でaj(t)−1個の0が補償できる。零補償の間、0の要素および対応するaj(t)は逐次的に更新される。「0」が補償されるたびに、a'におけるその値が「1」に変更され、a'におけるaj(t)は1減らされる。
【００６６】
結果として得られるベクトルa'は、aについての零補償された代表ベクトルである。k≧nとして任意の区間[vk,＋∞)についての零補償された代表ベクトルが、npをkで置き換えることによって同様に得ることができることを注意しておく。
【００６７】
たとえば、図１に挙げた場合についての代表ベクトルおよびそれらの零補償された代表ベクトルはそれぞれ、[m≧3]，[m≧1,0,2]，[m≧1,2]，[m≧2,1]，[m≧1,1,1]および[[1,1,m≧1]]，[[m≧1,1,1]]，[[1,m≧1,1]]，[[1,m≧1,1]][[m≧1,1,1]]である。明らかに、代表ベクトルから対応する零補償された代表ベクトルへのマッピングは一意的である。一方、逆マッピングについてはこれは成り立たない可能性がある。
【００６８】
有利なことに、そのような零補償された代表ベクトルの導入により、複数のピークを含む相関結果値を比較することがずっと簡単になる。具体的には、零補償された代表ベクトルに0の値がなければ、
【００６９】
【数３０】
2016515721.tif 000032
と記されるN個の相関値を集めた相関ベクトルがnp個のピークより大きいまたはそれに等しいと想定できる。このことは、簡潔に
【００７０】
【数３１】
2016515721.tif 000033
〔またはc1N≧v1np〕と記される。さらに、零補償された代表ベクトルにk個の0があれば、c1Nはちょうどnp−k個の、≧v1npである値をもつと述べることができる。図１に示した値分布によれば、個々のイベントについての確率は次のように計算される。
【００７１】
【数３２】
2016515721.tif 000034
最終的なFP確率はPFP＝P1＋P2＋P3＋P4＋P5である。因子(pi−pi-1)は、範囲[vi,vi-1]内の値をもつ相関結果ピークを得る確率であり、一方、因子(1−pi)は、範囲(−∞,vi]内にピークを得る確率である。
【００７２】
〈Ｂ．限られた範囲内での相関結果ピーク値についてのFP確率〉
上述したように、音響経路を通じた信号伝送については、相関結果ピーク量値は一般に、相関長Nよりずっと小さい最大サイズLの限られた時間的範囲内に集中する。よって、N個の相関結果値のうちの長さLの現在の窓内のnp個のピーク量値について、マークされないコンテンツついての相関結果ピーク値についての平均または期待される確率分布において、この現在の長さLの窓内でnp個のピーク量値より大きいまたはそれに等しいnp個以上の相関結果量値がある確率が計算される。
【００７３】
〈Ｂ．１  発明のFP確率計算処理〉
【００７４】
【数３３】
2016515721.tif 000035
についての定義は上記で与えてある。
【００７５】
スライディング窓が一回または複数回、np個の期待されるピーク量値より大きいまたはそれに等しいnp個以上の値を含むことの相補的な場合は、np個の期待されるピーク量値より大きいまたはそれに等しいnp個以上の値を含むスライディング窓がない、すなわち
【００７６】
【数３４】
2016515721.tif 000036
ということである。よって、FP確率についての相補的な確率は次のようになる。
【００７７】
【数３５】
2016515721.tif 000037
連鎖律を使って、同時確率
【００７８】
【数３６】
2016515721.tif 000038
は、条件付き確率によって計算できる。
【００７９】
【数３７】
2016515721.tif 000039
相関ベクトルcjについては、直前のものcj-1のみが条件付き確率
【００８０】
【数３８】
2016515721.tif 000040
において有意である。該直前のものが、cjの一つの新たな要素を除いてすべての要素を含んでいるからである。
【００８１】
【数３９】
2016515721.tif 000041
したがって、同時確率
【００８２】
【数４０】
2016515721.tif 000042
は次のように書ける。
【００８３】
【数４１】
2016515721.tif 000043
さらに、Lが十分大きく、恒等式
【００８４】
【数４２】
2016515721.tif 000044
が成り立つので、L個の相関サンプルの部分集合がすでにN個のサンプル全部の代表的な集合を表わす。上記の恒等式が式(2)に用いられる。
【００８５】
【数４３】
2016515721.tif 000045
いかにして
【００８６】
【数４４】
2016515721.tif 000046
を評価するかは知られているので、FP確率計算は条件付き確率
【００８７】
【数４５】
2016515721.tif 000047
の評価に帰着する。
【００８８】
〈Ｂ．２  条件付き確率の評価〉
条件付き確率
【００８９】
【数４６】
2016515721.tif 000048
は、条件付き確率の定義を使って再定式化できる。二つのイベントAおよびBが与えられ、P(B)＞0とすると、Bが与えられたときのAの条件付き確率は次のように定義される。
【００９０】
【数４７】
2016515721.tif 000049
したがって、
【００９１】
【数４８】
2016515721.tif 000050
は次のように書くことができる。
【００９２】
【数４９】
2016515721.tif 000051
よって、FP確率は次のように評価できる。
【００９３】
【数５０】
2016515721.tif 000052
この一般的な公式は、L＝Nについては、特許文献１の統計的透かし検出処理における再帰的計算の場合（その場合、
【００９４】
【数５１】
2016515721.tif 000053
である）に帰着する。
【００９５】
【数５２】
2016515721.tif 000054
L≠Nの場合について式(6)を計算するために、次節で説明するように分割が実行される。
【００９６】
〈Ｂ．３  相関値分布に基づく同時確率評価〉
同時確率
【００９７】
【数５３】
2016515721.tif 000055
は次のように表現できる。
【００９８】
【数５４】
2016515721.tif 000056
c2Lがちょうど(np−1)個の≧vである値をもつ場合は、ここでもまた二つの分離した群に分割される。
【００９９】
・c2Lがちょうど(np−1)個の≧v1np-1である値を区間[vnp-1,＋∞)にもち、区間[vnp,vnp-1)には値をもたない。この場合、区間[vnp-1,＋∞)についての対応する零補償された代表ベクトルは0の要素をもたず、同時に「1」の要素のみをもつ。さもなければ、区間[vnp,vnp-1)における0が補償され、c2Lがちょうど(np−1)個の≧vである値をもつという制約条件に矛盾する。
【０１００】
・c2Lがちょうど(np−2)個の≧v1np-1である値を区間[vnp-1,＋∞)にもち、ちょうど(np−1)個の≧vである値を区間[vnp,＋∞)にもつ。これは、値分布についての何らかの区間における「m≧b」の場合に対応する。このことは次の例から明白である。
【０１０１】
np＝3の場合、c2Lはちょうどnp−2＝1個の≧v1np-1である値を区間[v2,＋∞)にもつ。これは、区間[v2,＋∞)についてのすべての零補償された代表ベクトルが一つの0要素をもつことを意味する。より具体的には、[[0,1]]および[[m≧1,0]]がnp＝3についての対応する零補償された代表ベクトルである。
【０１０２】
・[[0,1]]については、c2Lがちょうど(np−1)＝2個の≧vである値を区間[v3,＋∞)にもつためには、区間[v3,v2)内に一つまたは複数の値がある必要がある。よって、区間[v3,＋∞)についての零補償された代表ベクトルは[[m≧1,0,1]]となる。つまり、区間[v3,v2)においてm≧1個の値がある。
【０１０３】
・[[m≧1,0]]については、[[1,0]]または[[m≧2,0]]にさらに区別できる。[[1,0]]については、c2Lがちょうど(np−1)＝2個の≧vである値を区間[v3,＋∞)にもつよう、零補償された代表ベクトルは[[m≧1,1,0]]である必要がある。[[m≧2,0]]については、対応する代表ベクトルは[m≧2,0]であり、一方、c2Lがちょうど(np−1)＝2個の≧vである値をもつという条件を満たすために区間[v3,v2)内の値の数に対するさらなる要件は必要ない。
【０１０４】
図２は、c2Lはちょうど二つの、≧vである値をもつ場合を列挙している。P1およびP2は区間[v3,v2]においてc2Lが値をもたない場合に対応する。場合P3、P4およびP5は一つの区間において「m≧b」をもつ。図２の五つの場合について、すべての零補償された代表ベクトルは単一の0要素のみを含む。よって、cL+1をc2Lに加えるとc2≧vとなり、c1をc2Lに加えると
【０１０５】
【数５５】
2016515721.tif 000057
となるようcL+1およびc1についての値範囲を選ぶことはストレートにできる。
【０１０６】
たとえば、P2は零補償された代表ベクトル[[0,1,1]]に対応する。c1＜v3であれば、c1をc2Lに加えると
【０１０７】
【数５６】
2016515721.tif 000058
となる。他方、cL+1≧v3であれば、cL+1をc2Lに加えるとc2L+1≧v13となる。
【０１０８】
よって、個々の確率は次のように計算される。
【０１０９】
【数５７】
2016515721.tif 000059
一般に、P3、P4、P5確率についての和の項を計算することは、二項定理を用いることによって再定式化できる。
【０１１０】
【数５８】
2016515721.tif 000060
これはn1≪n2−n1であれば計算量の有意に軽減する。
【０１１１】
さらに、図２に示したP3、P4、P5のような「m≧b」の場合について、計算量を軽減するために、和についてのルックアップテーブルが構築されることができる。
【０１１２】
図３は、「m≧b」の場合についての一般的な値分布を示している。区間[vq,∞)における値の数はsで表わされ、区間[vq+1,vq)にはm個の値がある。よって、1≦q≦np−1かつ0≦s≦np−2かつb＝np−1−sである。ピークの数npが与えられると、qおよびsに依存するルックアップテーブルは次のように構築される。
【０１１３】
【数５９】
2016515721.tif 000061
c2Lがちょうど(np−1)個の≧vである値をもつすべての分離した場合について代表ベクトルを有すると、FP確率を評価することはストレートにできる。
【０１１４】
〈Ｂ．４  再帰的な代表ベクトル構築〉
代表ベクトルおよび零補償されたベクトルは、偽陽性確率
【０１１５】
【数６０】
2016515721.tif 000062
を計算するために使われる。これらのベクトルが既知であれば、確率P1ないしP5についての式が定式化されることができる。図２についての記述における例を参照。このように、FP(visk)の計算は、所与のピーク数npについてすべてのベクトルを構築する問題にマッピングされる。
【０１１６】
以下では、これらの代表ベクトルをいかにして再帰的に得るかを説明する。先に論じたように、図２におけるすべての場合は、対応する代表ベクトル[0,1,1]，[0,0,2]，[m≧2,0]，[m≧1,1,0]，[m≧1,0,1]によって表わせる。最初の二つのベクトルは再帰的に決定できる（上記参照）。最後の三つのベクトルの再帰的構築に関し、最後の三つのベクトルについての二つの属性は次のようなものである。
・「m≧b」の要素が常に第一の要素である。
・すべての要素の和がm≧2になる。
【０１１７】
換言すれば、代表ベクトルにおける「m≧b」要素を除いた残りの要素を与えられると、「m≧b」要素が推定できる。図３を参照するに、sが与えられると、b＝np−1−sである。よって、図２におけるP3、P4、P5の場合についての代表ベクトルは[0]，[1,0]，[0,1]に帰着する。そして問題は、np個のピークについていかにしてこれらの簡略化された代表ベクトルを見出すことに変えられる。
【０１１８】
図３に動機付けられて、これらの代表ベクトルはs∈{0,1,…,np−2}の値に従って異なる集合に分割される。sについて諸代表ベクトルが与えられると、s＋1についての対応する代表ベクトルは、拡張し、単位ベクトルを加えることによって構築できる。
【０１１９】
〈Ｂ．４．１  初期化〉
s＝0については代表ベクトルの集合は{[0]}であり、s＝1については代表ベクトルの集合は{[1,0],[1,0]}である。
【０１２０】
〈Ｂ．４．２  1≦s≦np−3についての更新〉
sについての集合における個々の代表ベクトルの長さljに依存して、s＋1についての現在の再帰の代表ベクトルは異なる仕方で構築される。
【０１２１】
・長さlj≦s＋1であれば、代表ベクトル中の最初の零でない要素を見出す。iが最初の零でない要素の位置であるとする。代表ベクトルに単位ベクトル
【０１２２】
【数６１】
2016515721.tif 000063
を加えて新しい代表ベクトルを得る。例示的な単位ベクトルはei＝(0,…,0,1,0,…,0)である。ここで、最初の0は位置「1」をもち、二番目の0は位置「i−1」をもち、「1」は位置「i」をもち、三番目の0は位置「i＋1」をもち、最後の0は位置「lj」をもつ。
【０１２３】
・長さlj＝s＋1であれば、先頭に一つの0を加えることによって代表ベクトルを拡張する。拡張された代表ベクトルにおいて最初の零でない値を見出し、単位ベクトルを加える。
【０１２４】
しかしながら、ある代表ベクトルについての零補償された代表ベクトルにおける「0」要素がその代表ベクトルにおける最初の零でない要素の右側にある場合には、拡張することなく、この代表ベクトルについて、上記の見出し、加える手順を実行する。
【０１２５】
〈Ｂ．４．３  例１〉
s＝1についていずれも長さlj＝s＋1＝2をもつ二つの代表ベクトル[1,0]および[0,1]を更新。ベクトル[0,1]については、最初の零でない要素の後には0要素はない。よって、新しい代表ベクトルは
【０１２６】
【数６２】
2016515721.tif 000064
として得られる。他方、ベクトル[1,0]については、その最初の零でない要素の後に一つの0要素があり、零補償後に一つの0が残る。よって、新しい代表ベクトルは、拡張なしでの単位ベクトルを加える
【０１２７】
【数６３】
2016515721.tif 000065
のほかに拡張して単位ベクトルを加える
【０１２８】
【数６４】
2016515721.tif 000066
がある。
【０１２９】
〈Ｂ．４．４  例２〉
s＝2について代表ベクトル集合{[0,0,2]，[0,1,1]，[1,0,1]，[0,2,0]，[1,1,0]，[2,0]}に関して更新。上記に従って、s＝3についての代表ベクトルは次のように得られる。
【０１３０】
【数６５】
2016515721.tif 000067
[0,1,0,2]については、これは[[0,1,1,1]]をその零補償された代表ベクトルとしており、零補償後はその最初の零でない要素の右側に0はない。一方、[1,0,0,2]についてはその零補償された代表ベクトルは[[1,0,1,1]]であり、零補償後でもその最初の零でない要素の右側に一つ0が残っている。よって、s＝4についての代表ベクトルを得るためには、[0,1,0,2]と[1,0,0,2]について異なる更新手順が実行される。
【０１３１】
図４における本発明の透かしデコーダのブロック図において、受領されたオーディオ信号は、取得または受領セクション・ステップまたはステージ４１においてサンプリングされ、その後、前処理ステップまたはステージ４２（これはたとえばスペクトル白色化を実行できる）を通過して、相関ステップまたはステージ４３に至り、ここで受領されたオーディオ信号の現在セクションが、参照パターン信号生成ステップまたはステージ４６からの、対応する異なる透かしシンボルに関係した、少なくとも二つの異なる候補参照パターン信号isと相関付けられ、それによりそれぞれの場合において、各候補透かしシンボルisについてN個の現在の相関結果値cisを生成する。合計nSymbols個の異なる候補参照パターン信号isがあり、たとえばnSymbols＝2である。秘密鍵を使って、ランダム位相生成ステップまたはステージ４５において対応する擬似ランダムな位相信号を生成し、その擬似ランダムな位相信号がステップ／ステージ４６において参照パターン信号生成を制御する。ステップ／ステージ４３における相関付けに続く透かしシンボル検出ステップまたはステージ４４では、対応する検出された透かしシンボルを決定し、出力するために、現在のN個の相関結果値cisのM個の主相関結果ピークのまわりで長さLの窓が検査され、上記および図５との関連で記述されるような処理が実行される。
【０１３２】
図５の流れ図では、ステップ５０で受領される入力相関ブロックcisについて、パラメータisが「1」に設定される。たとえば、is＝1は候補透かしシンボル「0」に対応し、is＝2は候補透かしシンボル「1」に対応する。各候補透かしシンボルisについて、対応する相関結果値の集合cisがステップ５１においてその大きさに従ってソートされ、最大のM個の量値cisk（すなわち、グローバル・ピーク量値）が決定される。たとえば、M＝3、k＝1,2,3である。ステップ５２では、各グローバル・ピーク量値k＝1,…,Mについておよび現在の候補透かしシンボルisについて、各グローバル・ピークのまわりの（すなわち、各グローバル・ピークを含む）もとの相関結果値の±L/2の範囲内（すなわち長さLの窓内）のnp個のピーク量値の対応するベクトルviskが決定される。比較ステップ５３によって制御されて、パラメータisは対応するループにおいて「1」から「nSymbols」まで走る。
【０１３３】
LはNより著しく小さいので、すべての長さLの窓の合計長さは長さNより小さい。実際上、L≪Nであり、LはNより少なくとも一桁小さい。すなわち、N/L＞10である。たとえば、N＝16kであり、L＝1kである。
【０１３４】
その後、is＝1からnSymbolsまで走る外側ループおよびk＝1からMまで走る内側ループにはいり、それぞれ比較ステップ５７および５６によって制御される。内側ループでは、偽陽性確率FP(visk)がステップ５４において現在の相関結果値から計算され、それに比較ステップ５５、kのインクリメント、比較ステップ５６および外側ループでのisのインクリメントおよび比較ステップ５７が続く。偽陽性確率FP(visk)値は、一回または複数回、そのような長さLの窓が、透かし入れされていないオーディオ信号についての相関結果量値についての平均または期待される確率分布においてnp個のピーク値より大きいまたはそれに等しいnp個以上の相関結果量値を含む確率に対応する。比較ステップ５７に続いて、候補透かしシンボルisがステップ５８において決定される。その候補透かしシンボルについて、最終的な偽陽性確率fpはminis(mink(FP(isk)))である。ステップ５９は、そのfp値をさらなる閾値Tmaxと比較することによって、透かしシンボルが検出されたかどうかを検査する。それが真であれば、検出された透かしシンボルが出力される。真でなければ、透かしシンボルは検出されなかったことになる。
【０１３５】
比較ステップ５５において、FP(visk)が所定の閾値Tminより小さければ、正しい候補透かしシンボルisがみつかったと想定され、計算時間節約のために両方のループを抜け、その透かしシンボルが出力される。
【０１３６】
図６は、長さNの例示的な相関結果および描かれているk＝1,2,3の三つの長さLの窓を示している。それぞれの長さLの窓はN個の相関結果値におけるピークのまわりに配置されている。N内の最初の可能な長さLの窓の位置はj＝1であり、最後の可能な長さLの窓の位置はj＝N−L＋1である。
【０１３７】
本発明の処理は、単一のプロセッサまたは電子回路によって、あるいは並列して動作するおよび／または本発明の異なる部分に対して動作するいくつかのプロセッサまたは電子回路によって実行されることができる。
【図１】
2016515721.tif 000074
【図２】
2016515721.tif 000075
【図３】
2016515721.tif 000076
【図４】
2016515721.tif 000077
【図５】
2016515721.tif 000078
【図６】
2016515721.tif 000079
【手続補正書】
【提出日】平成27年10月9日(2015.10.9)
【手続補正１】
【補正対象書類名】特許請求の範囲
【補正対象項目名】全文
【補正方法】変更
【補正の内容】
【特許請求の範囲】
【請求項１】
エコー、残響および／またはノイズを含みうる受領オーディオ信号中の透かしシンボルを判別する方法であって、
・前記受領オーディオ信号のセクションの、対応する異なる候補透かしシンボルisに関係した少なくとも二つの異なる候補参照パターン信号との相関を調べ、各候補透かしシンボルisについてそれぞれの場合にN個の相関結果値
【数６６】
2016515721.tif 000080
を生成するステップと、
・各候補透かしシンボルisについて、前記相関結果値cis内のM個のピーク値を決定するステップとを含んでおり、
・グローバル・ピーク値と記される前記M個のピーク値のそれぞれについておよび各候補透かしシンボルisについて、前記相関結果値cisから、長さLのM個の窓のそれぞれの中のnp個のピーク値の対応するベクトル
【数６７】
2016515721.tif 000081
を決定するステップであって、各長さLの窓は前記M個のグローバル・ピーク値のうちの一つを含み、長さLの窓全部の合計長さはNより小さく、k＝1,…,Mである、ステップと、
・長さLの前記M個の窓のそれぞれについて、ピーク値の前記ベクトルviskおよび透かし入れされていないオーディオ信号についての相関結果値についての平均のまたは期待される確率分布に基づいて、偽陽性確率値
【数６８】
2016515721.tif 000082
を計算するステップであって、前記偽陽性確率FP(visk)は、一回または複数回、そのような長さLの窓が、透かし入れされていないオーディオ信号についての相関値についての平均または期待される確率分布において前記np個のピーク値より大きいまたはそれに等しいnp個以上の相関結果絶対値を含むを含む確率として定義される、ステップと、
・前記受領オーディオ信号の前記セクションについて、最小の偽陽性確率値minis(mink(FP(visk)))をもつ候補透かしシンボルisを、該最小の偽陽性確率値minis(mink(FP(visk)))が所定の閾値（Tmax）より小さい場合には、検出された透かしシンボルとして選択し、そうでない場合には透かしシンボルが検出されなかったと判断するステップとを特徴とする、
方法。
【請求項２】
前記偽陽性確率FP(visk)値の前記計算において、その値がさらなる所定の閾値（Tmin）より小さい場合には、正しい候補透かしシンボルisが見つかったと想定され、前記のさらなる処理が省略され、前記所定の閾値（Tmax）は前記さらなる所定の閾値（Tmin）より大きい、請求項１記載の方法。
【請求項３】
M＝3である、請求項１または２記載の方法。
【請求項４】
二つの異なる候補透かしシンボルがある、請求項１ないし３のうちいずれか一項記載の方法。
【請求項５】
エコー、残響および／またはノイズを含みうる受領オーディオ信号中の透かしシンボルを判別する装置であって、
・前記受領オーディオ信号のセクションの、対応する異なる候補透かしシンボルisに関係した少なくとも二つの異なる候補参照パターン信号との相関を調べ、各候補透かしシンボルisについてそれぞれの場合にN個の相関結果値
【数６９】
2016515721.tif 000083
を生成するよう適応された手段と、
・各候補透かしシンボルisについて、前記相関結果値cis内のM個のピーク値を決定し、各候補透かしシンボルisについての前記M個の相関結果ピーク値から導出されて、前記受領オーディオ信号の前記セクションについて、前記候補透かしシンボルisのうちの透かしシンボルを決定するよう適応された手段とを含んでおり、
グローバル・ピーク値と記される前記M個のピーク値のそれぞれについておよび各候補透かしシンボルisについて、前記相関結果値cisから、長さLのM個の窓のそれぞれの中のnp個のピーク値の対応するベクトル
【数７０】
2016515721.tif 000084
が決定され、各長さLの窓は前記M個のグローバル・ピーク値のうちの一つを含み、長さLの窓全部の合計長さはNより小さく、k＝1,…,Mであり、
長さLの前記M個の窓のそれぞれについて、ピーク値の前記ベクトルviskおよび透かし入れされていないオーディオ信号についての相関結果値についての平均のまたは期待される確率分布に基づいて、偽陽性確率値
【数７１】
2016515721.tif 000085
が計算され、ここで、前記偽陽性確率FP(visk)は、一回または複数回、そのような長さLの窓が、透かし入れされていないオーディオ信号についての相関値についての平均または期待される確率分布において前記np個のピーク値より大きいまたはそれに等しいnp個以上の相関結果絶対値を含むを含む確率として定義され、
前記受領オーディオ信号の前記セクションについて、最小の偽陽性確率値minis(mink(FP(visk)))をもつ候補透かしシンボルisが、該最小の偽陽性確率値minis(mink(FP(visk)))が所定の閾値（Tmax）より小さい場合には、検出された透かしシンボルとして選択され、
そうでない場合には透かしシンボルが検出されなかったと判断される、
装置。
【請求項６】
前記偽陽性確率FP(visk)値の前記計算において、その値がさらなる所定の閾値（Tmin）より小さい場合には、正しい候補透かしシンボルisが見つかったと想定され、前記のさらなる処理が省略され、前記所定の閾値（Tmax）は前記さらなる所定の閾値（Tmin）より大きい、請求項５記載の装置。
【請求項７】
M＝3である、請求項５または６記載の装置。
【請求項８】
二つの異なる候補透かしシンボルがある、請求項５ないし７のうちいずれか一項記載の装置。
【国際調査報告】
2016515721.tif 000086
2016515721.tif 000087
2016515721.tif 000088