(57)【要約】【課題】電子透かしの検出の信頼性を落とさず、また電子透かしを視覚的に目立たせることもなく、より多くの情報を埋め込むことのできる電子透かし埋め込み装置を提供する。【解決手段】電子透かし埋め込み装置１において、埋め込み系列生成部１０は、埋め込み情報を、変位によって変化するＫ種（Ｋは、３以上の整数）の異なる値に含まれるいずれかの予め対応付けられる値に置き換えて情報系列を生成し、拡散系列に基づいて、生成した情報系列を変位させて埋め込み系列を生成する。パターン重畳部４０は、埋め込み系列に基づいて入力信号に対して電子透かしを重畳して透かし埋め込み済み信号Ｉ’を出力する。【選択図】図１
2014022972.tif 000002
【特許請求の範囲】
【請求項１】
埋め込み情報を電子透かしとして人間の知覚に感知されないように入力信号に対して埋め込む電子透かし埋め込み装置であって、
前記埋め込み情報の全体、または前記埋め込み情報を任意の長さに分割した各部分を、変位によって変化するＫ種（Ｋは、３以上の整数）の異なる値に含まれるいずれかの予め対応付けられる値に置き換えて情報系列を生成し、拡散系列に基づいて、生成した前記情報系列を変位させて埋め込み系列を生成する埋め込み系列生成部と、
前記埋め込み系列に基づいて前記入力信号に対して電子透かしを重畳して電子透かし埋め込み済み信号を生成するパターン重畳部と、
を備えることを特徴とする電子透かし埋め込み装置。
【請求項２】
前記変位によって変化するＫ種（Ｋは、３以上の整数）の異なる値は、複素数であり、
前記埋め込み系列生成部は、
前記埋め込み情報の全体、または前記埋め込み情報を任意の長さに分割した各部分であるシンボルを、Ｋ種（Ｋは、３以上の整数）の異なる複素数の値に含まれるいずれかの予め対応付けられる複素数に置き換えて複素数の情報系列を生成し、生成した前記複素数の情報系列の各要素と複素数の拡散系列の各要素に基づいて前記複素数の情報系列の位相または偏角を変位させて埋め込み系列を生成する
ことを特徴とする請求項１に記載の電子透かし埋め込み装置。
【請求項３】
前記入力信号は、Ｎ（Ｎは、２以上の整数）以上の次元を有する信号であり、
前記埋め込み系列生成部が生成した前記埋め込み系列に基づいてＮ−１次元の埋め込みパターンを生成する埋め込みパターン生成部と、
前記Ｎ−１次元の埋め込みパターンに含まれる値に基づいて変調を行いＮ次元の透かし情報パターンを生成する変調部と、を備え、
前記パターン重畳部は、
前記入力信号に対して前記Ｎ次元の透かし情報パターンを重畳して電子透かし埋め込み済み信号を生成する
ことを特徴とする請求項１または２に記載の電子透かし埋め込み装置。
【請求項４】
前記入力信号は、Ｎ（Ｎは、２以上の整数）以上の次元を有する信号であり、
前記埋め込み系列生成部が生成した前記埋め込み系列に基づいてＮ−１次元の埋め込みパターンを生成する埋め込みパターン生成部を備え、
前記パターン重畳部は、
直交変換された前記入力信号に前記埋め込みパターンを重畳した信号を生成し、
前記重畳した信号を直交逆変換して電子透かし埋め込み済み信号を生成する
ことを特徴とする請求項１または２に記載の電子透かし埋め込み装置。
【請求項５】
埋め込み情報を変位によって変化するＫ種（Ｋは、３以上の整数）の異なる値の組に含まれるいずれかの予め対応付けられる値に置き換えることにより生成された情報系列と拡散系列とに基づいて生成された電子透かしを入力信号から検出する電子透かし検出装置であって、
前記入力信号から変位によって変化する値の検出系列を生成する検出系列抽出部と、
前記検出系列と前記拡散系列に基づいて相関値を算出する相関値計算部と、
前記相関値と、前記Ｋ種の異なる値のそれぞれを区分する境界とに基づいて埋め込まれている情報系列を判定する透かし判定部と、
を備えることを特徴とする電子透かし検出装置。
【請求項６】
前記変位によって変化するＫ種（Ｋは、３以上の整数）の異なる値は、複素数であり、
前記透かし判定部は、
前記異なる複素数の値ごとに算出される前記相関値を変位させて同じ方向に向かせて前記相関値の総和を算出し、算出した前記相関値の総和の位相に基づいて新たに定められる前記Ｋ種の異なる複素数の値を区分する境界に基づいて埋め込まれている情報系列を判定する
ことを特徴とする請求項５に記載の電子透かし検出装置。
【請求項７】
埋め込み情報を電子透かしとして人間の知覚に感知されないように入力信号に対して埋め込む電子透かし埋め込み方法であって、
埋め込み系列生成部が、前記埋め込み情報の全体、または前記埋め込み情報を任意の長さに分割した各部分を、変位によって変化するＫ種（Ｋは、３以上の整数）の異なる値に含まれるいずれかの予め対応付けられる値に置き換えて情報系列を生成し、拡散系列に基づいて、生成した前記情報系列を変位させて埋め込み系列を生成するステップと、
パターン重畳部が、前記埋め込み系列に基づいて前記入力信号に対して電子透かしを重畳して電子透かし埋め込み済み信号を生成するステップと、
を含むことを特徴とする電子透かし埋め込み方法。
【請求項８】
埋め込み情報を変位によって変化するＫ種（Ｋは、３以上の整数）の異なる値の組に含まれるいずれかの予め対応付けられる値に置き換えることにより生成された情報系列と拡散系列とに基づいて生成された電子透かしを入力信号から検出する電子透かし検出方法であって、
検出系列抽出部が、前記入力信号から変位によって変化する値の検出系列を生成するステップと、
相関値計算部が、前記検出系列と前記拡散系列に基づいて相関値を算出するステップと、
透かし判定部が、前記相関値と、前記Ｋ種の異なる値のそれぞれを区分する境界とに基づいて埋め込まれている情報系列を判定するステップと、
を含むことを特徴とする電子透かし検出方法。
【請求項９】
埋め込み情報を電子透かしとして人間の知覚に感知されないように入力信号に対して埋め込むコンピュータに、
前記埋め込み情報の全体、または前記埋め込み情報を任意の長さに分割した各部分を、変位によって変化するＫ種（Ｋは、３以上の整数）の異なる値に含まれるいずれかの予め対応付けられる値に置き換えて情報系列を生成し、拡散系列に基づいて、生成した前記情報系列を変位させて埋め込み系列を生成するステップと、
前記埋め込み系列に基づいて前記入力信号に対して電子透かしを重畳して電子透かし埋め込み済み信号を生成するステップと、
を実行させるためのプログラム。
【請求項１０】
埋め込み情報を変位によって変化するＫ種（Ｋは、３以上の整数）の異なる値の組に含まれるいずれかの予め対応付けられる値に置き換えることにより生成された情報系列と拡散系列とに基づいて生成された電子透かしを入力信号から検出するコンピュータに、
前記入力信号から変位によって変化する値の検出系列を生成するステップと、
前記検出系列と前記拡散系列に基づいて相関値を算出するステップと、
前記相関値と、前記Ｋ種の異なる値のそれぞれを区分する境界とに基づいて埋め込まれている情報系列を判定するステップと、
を実行させるためのプログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、画像、映像信号などの信号に対し、透かしを知覚されないように埋め込み、また、透かしの埋め込まれた信号からこの透かしを読み取るための電子透かし埋め込み装置、電子透かし検出装置、方法及びプログラムに関する。
【背景技術】
【０００２】
近年、デジタルコンテンツに対し電子透かしを埋め込むことで、デジタルコンテンツの著作権保護を行ったり、デジタルコンテンツに関する著作権情報などのメタデータを参照したりする技術が実現されている。また、電子透かしが埋め込まれたデジタルコンテンツが、テレビのコマーシャルメッセージあるいは印刷物などのアナログ媒体を通じて広告として提供され、このデジタルコンテンツをデジタルカメラなどで撮影し、電子透かしを読み取ることで広告に関連した情報をユーザに取得させたりするといったサービスも実現されている。
【０００３】
（静止画電子透かし）
静止画像に電子透かしを埋め込む方法としては、例えば、特許文献1に記載の技術が知られている。特許文献１においては、擬似乱数によって生成した埋め込み系列を画像の直交変換領域（例えば、フーリエ変換領域）の実部および虚部に埋め込んでおき、埋め込み系列と検出系列の相関を用いて電子透かしの検出を行うスペクトル拡散型の電子透かし方式が開示されている。
【０００４】
（動画電子透かし）
映像信号に電子透かしを埋め込む方法としては、例えば、特許文献２に記載の技術が知られている。特許文献２においては、擬似乱数によって生成した埋め込み系列を時間方向周期信号によって変調して映像に重畳しておき、埋め込み系列と検出系列の複素相関を用いて電子透かしの検出を行うスペクトル拡散型の電子透かし方式が開示されている。
【先行技術文献】
【特許文献】
【０００５】
【特許文献１】特開２００３−２１９１４８号公報
【特許文献２】国際公開第２００７／１０２４０３号
【発明の概要】
【発明が解決しようとする課題】
【０００６】
ところで、特許文献２には、ｎビットの埋め込み情報をｍ倍の長さに拡散し、拡散された各ビットをＮ−１次元の空間上の複素数の実部と虚部に対応させることで、Ｎ−１次元空間上のｍ／２個の複素数を用いて元の１ビットの情報を埋め込む技術が開示されている。この技術において、信号に埋め込む情報を増やすためには、Ｎ−１次元空間上で埋め込みに用いる位置を増やす必要がある。すなわち、Ｎ−１次元の複素配列の配列数を増加させて、埋め込みに用いる複素数を増やす必要がある。
【０００７】
しかしながら、埋め込みに用いる複素数を増やすということは、埋め込みパターンのＮ−１次元空間上での空間周波数のより高い周波数領域を用いることになる。高い周波数の領域を用いる場合、再撮影などをすることにより解像度が落ちた場合に電子透かしの検出が困難になり、また電子透かしが視覚的に目立ちやすくなってしまうなどの問題がある。
【０００８】
本発明は、上記問題を解決すべくなされたもので、その目的は、電子透かしを埋め込む際に、電子透かしの検出の信頼性を落とさず、また電子透かしを視覚的に目立たせることもなく、より多くの情報を埋め込むことのできる電子透かし埋め込み装置、電子透かし検出装置、方法及びプログラムを提供することにある。
【課題を解決するための手段】
【０００９】
上記問題を解決するために、本発明の一態様は、埋め込み情報を電子透かしとして人間の知覚に感知されないように入力信号に対して埋め込む電子透かし埋め込み装置であって、前記埋め込み情報の全体、または前記埋め込み情報を任意の長さに分割した各部分を、変位によって変化するＫ種（Ｋは、３以上の整数）の異なる値に含まれるいずれかの予め対応付けられる値に置き換えて情報系列を生成し、拡散系列に基づいて、生成した前記情報系列を変位させて埋め込み系列を生成する埋め込み系列生成部と、前記埋め込み系列に基づいて前記入力信号に対して電子透かしを重畳して電子透かし埋め込み済み信号を生成するパターン重畳部と、を備えることを特徴とする電子透かし埋め込み装置である。
【００１０】
また、本発明の一態様は、上記に記載の発明において、前記変位によって変化するＫ種（Ｋは、３以上の整数）の異なる値は、複素数であり、前記埋め込み系列生成部は、前記埋め込み情報の全体、または前記埋め込み情報を任意の長さに分割した各部分であるシンボルを、Ｋ種（Ｋは、３以上の整数）の異なる複素数の値に含まれるいずれかの予め対応付けられる複素数に置き換えて複素数の情報系列を生成し、生成した前記複素数の情報系列の各要素と複素数の拡散系列の各要素に基づいて前記複素数の情報系列の位相または偏角を変位させて埋め込み系列を生成することを特徴とする。
【００１１】
また、本発明の一態様は、上記に記載の発明において、前記入力信号は、Ｎ（Ｎは、２以上の整数）以上の次元を有する信号であり、前記埋め込み系列生成部が生成した前記埋め込み系列に基づいてＮ−１次元の埋め込みパターンを生成する埋め込みパターン生成部と、前記Ｎ−１次元の埋め込みパターンに含まれる値に基づいて変調を行いＮ次元の透かし情報パターンを生成する変調部と、を備え、前記パターン重畳部は、前記入力信号に対して前記Ｎ次元の透かし情報パターンを重畳して電子透かし埋め込み済み信号を生成することを特徴とする。
【００１２】
また、本発明の一態様は、上記に記載の発明において、前記入力信号は、Ｎ（Ｎは、２以上の整数）以上の次元を有する信号であり、前記埋め込み系列生成部が生成した前記埋め込み系列に基づいてＮ−１次元の埋め込みパターンを生成する埋め込みパターン生成部を備え、前記パターン重畳部は、直交変換された前記入力信号に前記埋め込みパターンを重畳した信号を生成し、前記重畳した信号を直交逆変換して電子透かし埋め込み済み信号を生成することを特徴とする。
【００１３】
また、本発明の一態様は、埋め込み情報を変位によって変化するＫ種（Ｋは、３以上の整数）の異なる値の組に含まれるいずれかの予め対応付けられる値に置き換えることにより生成された情報系列と拡散系列とに基づいて生成された電子透かしを入力信号から検出する電子透かし検出装置であって、前記入力信号から変位によって変化する値の検出系列を生成する検出系列抽出部と、前記検出系列と前記拡散系列に基づいて相関値を算出する相関値計算部と、前記相関値と、前記Ｋ種の異なる値のそれぞれを区分する境界とに基づいて埋め込まれている情報系列を判定する透かし判定部と、を備えることを特徴とする電子透かし検出装置である。
【００１４】
また、本発明の一態様は、上記に記載の発明において、前記変位によって変化するＫ種（Ｋは、３以上の整数）の異なる値は、複素数であり、前記透かし判定部は、前記異なる複素数の値ごとに算出される前記相関値を変位させて同じ方向に向かせて前記相関値の総和を算出し、算出した前記相関値の総和の位相に基づいて新たに定められる前記Ｋ種の異なる複素数の値を区分する境界に基づいて埋め込まれている情報系列を判定することを特徴とする。
【００１５】
また、本発明の一態様は、埋め込み情報を電子透かしとして人間の知覚に感知されないように入力信号に対して埋め込む電子透かし埋め込み方法であって、埋め込み系列生成部が、前記埋め込み情報の全体、または前記埋め込み情報を任意の長さに分割した各部分を、変位によって変化するＫ種（Ｋは、３以上の整数）の異なる値に含まれるいずれかの予め対応付けられる値に置き換えて情報系列を生成し、拡散系列に基づいて、生成した前記情報系列を変位させて埋め込み系列を生成するステップと、パターン重畳部が、前記埋め込み系列に基づいて前記入力信号に対して電子透かしを重畳して電子透かし埋め込み済み信号を生成するステップと、を含むことを特徴とする電子透かし埋め込み方法である。
【００１６】
また、本発明の一態様は、埋め込み情報を変位によって変化するＫ種（Ｋは、３以上の整数）の異なる値の組に含まれるいずれかの予め対応付けられる値に置き換えることにより生成された情報系列と拡散系列とに基づいて生成された電子透かしを入力信号から検出する電子透かし検出方法であって、検出系列抽出部が、前記入力信号から変位によって変化する値の検出系列を生成するステップと、相関値計算部が、前記検出系列と前記拡散系列に基づいて相関値を算出するステップと、透かし判定部が、前記相関値と、前記Ｋ種の異なる値のそれぞれを区分する境界とに基づいて埋め込まれている情報系列を判定するステップと、を含むことを特徴とする電子透かし検出方法である。
【００１７】
また、本発明の一態様は、埋め込み情報を電子透かしとして人間の知覚に感知されないように入力信号に対して埋め込むコンピュータに、前記埋め込み情報の全体、または前記埋め込み情報を任意の長さに分割した各部分を、変位によって変化するＫ種（Ｋは、３以上の整数）の異なる値に含まれるいずれかの予め対応付けられる値に置き換えて情報系列を生成し、拡散系列に基づいて、生成した前記情報系列を変位させて埋め込み系列を生成するステップと、前記埋め込み系列に基づいて前記入力信号に対して電子透かしを重畳して電子透かし埋め込み済み信号を生成するステップと、を実行させるためのプログラムである。
【００１８】
また、本発明の一態様は、埋め込み情報を変位によって変化するＫ種（Ｋは、３以上の整数）の異なる値の組に含まれるいずれかの予め対応付けられる値に置き換えることにより生成された情報系列と拡散系列とに基づいて生成された電子透かしを入力信号から検出するコンピュータに、前記入力信号から変位によって変化する値の検出系列を生成するステップと、前記検出系列と前記拡散系列に基づいて相関値を算出するステップと、前記相関値と、前記Ｋ種の異なる値のそれぞれを区分する境界とに基づいて埋め込まれている情報系列を判定するステップと、を実行させるためのプログラムである。
【発明の効果】
【００１９】
この発明によれば、電子透かしを埋め込む際に、電子透かしの検出の信頼性を落とさず、また電子透かしを視覚的に目立たせることもなく、より多くの情報を埋め込むことが可能となる。
【図面の簡単な説明】
【００２０】
【図１】本発明の第１実施形態による電子透かし埋め込み装置の構成を示すブロック図である。
【図２】同実施形態による電子透かし埋め込み装置の埋め込み系列生成部の内部構成を示すブロック図である。
【図３】同実施形態による電子透かし埋め込みの処理を示したフローチャートである。
【図４】同実施形態によるシンボルと複素数の対応付け（その１）を示した図である。
【図５】同実施形態によるシンボルと複素数の対応付け（その２）を示した図である。
【図６】同実施形態による埋め込み系列の配列への設定例を示した図である。
【図７】同実施形態による他のパターン重畳を行う電子透かし埋め込み装置の構成を示すブロック図である。
【図８】本発明の第１実施形態による電子透かし検出装置の構成を示すブロック図である。
【図９】同実施形態による電子透かし検出の処理を示したフローチャートである。
【図１０】同実施形態による相関値の分布の例（その１）を示した図である。
【図１１】同実施形態による相関値の分布の例（その２）を示した図である。
【図１２】同実施形態による相関値の位相同期の例を示した図である。
【図１３】第２実施形態による電子透かし埋め込み装置の構成を示すブロック図である。
【図１４】第２実施形態による電子透かし検出装置の構成を示すブロック図である。
【図１５】第３実施形態による電子透かし埋め込み装置の構成を示すブロック図である。
【図１６】第３実施形態による電子透かし検出装置の構成を示すブロック図である。
【図１７】第４実施形態による電子透かし埋め込み装置の構成を示すブロック図である。
【図１８】同実施形態による電子透かし埋め込み装置の埋め込み系列生成部の内部構成を示すブロック図である。
【図１９】第４実施形態による電子透かし検出装置の構成を示すブロック図である。
【発明を実施するための形態】
【００２１】
以下、本発明の実施形態について図面を参照して説明する。
（第１実施形態）
第１実施形態では、電子透かしとして入力信号に埋め込む埋め込み情報を分割した部分列をシンボルとする。そして、各々のシンボルに異なる複素数を割り当てることにより複素情報系列を生成し、生成した複素情報系列に基づいて電子透かしを埋め込みまた電子透かしの検出を行う構成を示す。
【００２２】
（第１実施形態の電子透かし埋め込み装置）
図１は、第１実施形態による電子透かし埋め込み装置１の概略ブロック図である。電子透かし埋め込み装置１は、埋め込み系列生成部１０と、埋め込みパターン生成部２０と、時間変調部３０と、パターン重畳部４０を備えている。
【００２３】
電子透かし埋め込み装置１において、埋め込み系列生成部１０は、埋め込み情報Ｕを入力し、埋め込み系列ｗを出力する。埋め込みパターン生成部２０は、埋め込み系列生成部１０が生成した埋め込み系列ｗを入力し、埋め込み系列ｗに基づいて生成した埋め込みパターンＶを出力する。時間変調部３０は、埋め込みパターンＶを入力し、埋め込みパターンＶに基づいて生成した透かし映像パターンＷを出力する。パターン重畳部４０は、透かし埋め込み対象映像Ｉと、時間変調部３０が生成した透かし映像パターンＷとを入力し、透かし埋め込み対象映像Ｉに透かし映像パターンＷを重畳した透かし埋め込み済み映像Ｉ’を生成して出力する。
【００２４】
図２は、埋め込み系列生成部１０の内部ブロックを示した図である。埋め込み系列生成部１０は、複素情報系列生成部１１と、埋め込み系列計算部１２を備えている。
埋め込み系列生成部１０において、複素情報系列生成部１１は、埋め込み情報Ｕを入力し、入力した埋め込み情報Ｕに基づいて生成した複素情報系列ｓを出力する。埋め込み系列計算部１２は、複素情報系列生成部１１が生成した複素情報系列ｓを入力し、複素情報系列ｓに基づいて算出した埋め込み系列ｗを出力する。
【００２５】
なお、埋め込み情報Ｕおよび透かし埋め込み対象映像Ｉの入力は、外部の記憶装置に記憶されているものを読み出して入力するか、あるいは外部の装置が出力した情報を入力するなどのように構成することができる。また、埋め込み系列ｗ、埋め込みパターンＶ，透かし映像パターンＷ、複素情報系列ｓは、各機能部の間で直接入出力する構成の他、メモリ等の記憶手段を介するような構成であってもよい。すなわち、それぞれの機能部によって、これらの情報が生成されたときに、生成の処理を行った機能部が電子透かし埋め込み装置１に備えられるメモリなどの記憶手段に生成した情報を書き込み、次段の機能部がそれを読み出すことによって処理されるようにしてもよい。
【００２６】
次に、図３を参照しつつ、電子透かし埋め込み装置１による電子透かし埋め込みの処理について説明する。
埋め込み系列生成部１０は、埋め込みを行う埋め込み情報Ｕを入力し、長さＬの複素数の系列である埋め込み系列ｗを生成する（ステップＳ１０１：埋め込み系列生成ステップ）。埋め込み系列生成ステップは、ステップＳ１０１１とステップＳ１０１２の２段階を有する。
【００２７】
（ステップＳ１０１１の処理：複素情報系列の生成）
複素情報系列生成部１１が埋め込み情報Ｕに基づいて複素情報系列ｓを生成して出力する。以下、ステップＳ１０１１における複素情報系列生成部１１による複素情報系列の生成について詳細に説明する。
まず、埋め込み情報Ｕが、ｎビットから構成される情報、例えば、ｂ１，ｂ２，…，ｂｉ，…，ｂｎであり、ｂｉ∈｛−１，１｝である場合について説明する。
【００２８】
（１）複素情報系列生成部１１が、入力したｎビットの埋め込み情報Ｕを、ｋ個のシンボルＳ１，Ｓ２，…，Ｓｋに分割する。ここで、シンボルとは、それぞれ埋め込み情報Ｕのうちの一部分の情報を表し、実際の電子透かし埋め込みの処理単位となる情報である。全てのシンボルが同じ数、例えば、ｇビットずつであってもよいし、各シンボルが異なるビット数で構成されていてもよい。例えば、埋め込み情報Ｕが６４ビット長の２進数値で与えられたとき、２ビットずつの情報に区切って２ビットの情報の各々をシンボルとしてもよい。
【００２９】
（２）次に、複素情報系列生成部１１が、ｋ個のシンボルＳ１，Ｓ２，…，Ｓｋに対し、それぞれ複素数を割り当て、ｋ個の複素数の列ｓ１，ｓ２，…，ｓｋを生成する。シンボルＳｉに割り当てる複素数ｓｉは、各シンボルがｇビットで構成されている場合、Ｓｉに対応するビット値ｂｈ，ｂｈ＋１，…，ｂｈ＋ｇ−１の組み合わせに応じたそれぞれ異なる位相を有する複素数となる。ｇビットを１シンボルで表している場合、２ｇ種の複素数をシンボル値として用いることになる。
【００３０】
例えば、埋め込み情報Ｕがｎ＝２ｋビットの場合、ビット値ｂ２ｉ−１，ｂ２ｉがシンボルＳｉに対応しているとき、Ｓｉに対応する複素数ｓｉを下記の式（１）のように定めることができる。
【００３１】
【数１】
2014022972.tif 000003
【００３２】
ここで、ｊは虚数単位である。式（１）に基づくシンボルと複素数の対応付けを複素平面でみた場合、シンボルＳｉに対応する複素数ｓｉは、図４（ａ）のように示されることになる。
【００３３】
他の複素数の定め方として、埋め込み情報Ｕが、ｎ＝３ｋビットの場合、ビット値ｂ３ｉ−２，ｂ３ｉ−１，ｂ３ｉがシンボルＳｉに対応しているとき、Ｓｉに対応する複素数ｓｉを図４（ｂ）のようにビット値の組に応じて複素数の位相を変えるようにして定めることができる。
【００３４】
なお、図４（ｂ）のように定めた場合、互いに位相が隣接するシンボルの複素数に対応するビット値が、高々１ビットしか異ならないように、グレイコードを用いて表現してもよい。このように表現した場合、検出時にエラーが発生した場合のビットエラーレートを低く抑えることが可能となる。
また、上記の例において、ｋ＝１としてｎビットの情報を１つのシンボルで表現するようにすることもできる。このとき、１つのシンボルで２ｎ種類の複素数が表現されることになる。
【００３５】
さらに、他の複素数の定め方として、埋め込み情報Ｕをｔ進数に変換し、ｔ進数の１桁を１シンボルに割り当てるようにしてもよい。例えば、埋め込み情報Ｕを３進数で表現する。３進数は、０、１、２のいずれかの数値で示されるので、これらの数値に図５に示すように３種類の複素数を割り当てておく。そして、３進数で表された埋め込み情報Ｕの一桁ｔｉ∈｛０，１，２｝を１シンボルを割り当てるようにすれば、３種類の複素数をシンボルに対応する複素数として定めることができる。
【００３６】
以上の手順により、ステップＳ１０１１において複素情報系列生成部１１が、｛ｓ１，ｓ２，…，ｓｋ｝を複素情報系列ｓとして生成して出力する。
なお、埋め込み情報Ｕは、上述したビット構成、すなわちｂ１，ｂ２，…，ｂｉ，…，ｂｎは、ｂｉ∈｛−１，１｝に限られるものではなく、他のビット構成、例えば、一般的な０と１からなるビット構成などであっても複素情報系列ｓを生成することは可能である。
【００３７】
（ステップＳ１０１２の処理：埋め込み系列の生成）
次に、埋め込み系列計算部１２が、複素情報系列生成部１１が出力した複素情報系列ｓを入力し、入力した複素情報系列ｓを擬似乱数列で拡散して埋め込み系列ｗを算出して出力する。以下、埋め込み系列計算部１２による埋め込み系列ｗの算出について詳細に説明する。
【００３８】
（１）埋め込み系列計算部１２は、入力した複素情報系列であるｓを複素数の擬似乱数列ＰＮ＝｛ＰＮ１，ＰＮ２，…，ＰＮｍｋ｝で拡散して埋め込み系列ｗを算出する。
具体的には、まず、長さｋの複素情報系列ｓ１，ｓ２，…，ｓｋから、各複素数をｍ回ずつ繰り返した式（２）によって表される系列ｓ（ｍ）を生成する。このとき、上述した埋め込み系列ｗの長さＬは、Ｌ＝ｍｋとなる。
【００３９】
【数２】
2014022972.tif 000004
【００４０】
次に、埋め込み系列計算部１２は、埋め込み系列ｗ＝｛ｗ１，ｗ２，…，ｗｍｋ｝をｓ（ｍ）とＰＮの要素ごとの積として式（３）のように算出する。
【００４１】
【数３】
2014022972.tif 000005
【００４２】
なお、複素数の擬似乱数列ＰＮは、Ｍ系列やＧＯＬＤ系列から得られる値を用いて生成することができる。また、擬似乱数列ＰＮは、処理の際に演算によって生成する構成の他、予め演算によって生成しておき、内部のメモリ等の記憶手段に予め記憶させておいたものを読み出すことによって利用する構成であってもよい。
【００４３】
（埋め込み系列の他の表現形式）
なお、埋め込み系列ｗは、上記の複素情報系列ｓおよび埋め込み系列ｗを複素数値の系列で表す形式に限られず、埋め込みを行う透かし情報に応じて定めることができる系列であれば他の形式で表すこともできる。例えば、各々の複素数の偏角を表す値の列で表すという形式がある。複素情報系列ｓ、埋め込み系列ｗ、擬似乱数列ＰＮのそれぞれの列が下式（４）のような複素数値で表されるとする。
【００４４】
【数４】
2014022972.tif 000006
【００４５】
複素数の偏角の値の列で表す複素情報系列＾ｓ、埋め込み系列＾ｗ、擬似乱数列＾ＰＮは、複素数の偏角、θ、φ、ψを用いて式（５）のように表されることになる。
【００４６】
【数５】
2014022972.tif 000007
【００４７】
なお、＾ｓ、＾ｗ、＾ＰＮは、それぞれ式（５）の左辺の文字に対応する文字であり、「＾」の記号は、本明細書中において「＾Ｃ」のように記載した場合、式（５）の左辺のように文字「Ｃ」の上に「＾」が記載されているものであると解釈するものとする。また、前述したφは、次の式（６）の文字と同じギリシャ文字のファイであり、本明細書中ではφとして記載するが、式中に示される（６）の文字と同じ意味として解釈するものとする。
【００４８】
【数６】
2014022972.tif 000008
【００４９】
複素数の偏角で表した場合、埋め込み系列＾ｗの算出は、次式（７）のように行われる。
【００５０】
【数７】
2014022972.tif 000009
【００５１】
上記のように、複素情報系列ｓとしては、シンボルごとの位相もしくは方向を表す値の列であればどのような値を用いてもよい。
また、擬似乱数列ＰＮは、複素情報系列ｓを拡散するための拡散系列であり、拡散系列の要素毎の位相もしくは方向をあらわす値の列であればどのような値を用いてもよい。
また、埋め込み系列ｗは、擬似乱数列ＰＮの位相もしくは方向を、複素情報系列ｓの位相もしくは方向によって要素毎に変化させることで得られる系列であればどのような値で表現されていてもよい。
【００５２】
（ステップＳ１０２の処理：埋め込みパターンの生成）
次に、図３に戻り、埋め込みパターン生成部２０が、埋め込み系列生成部１０が生成した埋め込み系列ｗに基づいて埋め込みパターンＶを生成する（ステップＳ１０２）。以下、ステップＳ１０２を３つの段階（ステップＳ１０２−１、Ｓ１０２−２、Ｓ１０２−３）に分けて説明する。
【００５３】
（ステップＳ１０２−１）
最初に、埋め込みパターン生成部２０は、大きさＸ１×Ｘ２×…×ＸＮ−１のＮ−１次元複素配列の領域を確保する。ここで、Ｘ１，Ｘ２，…，ＸＮ−１はあらかじめ定められた要素数である。
【００５４】
（ステップＳ１０２−２）
次に、埋め込みパターン生成部２０は、埋め込み系列ｗから、順に１つずつ、合計Ｍ個の値を読み出し、Ｎ−１次元複素配列のＭ個の配列の位置［０，０，…，０］から順に、読み出した値を配列の要素値として配列に値を書き込んで設定していく。すなわち、複素配列が、Ａ［ｐ１，ｐ２，…ｐｎ，…，ｐＮ−１］（ｐｎ≧０）であり、埋め込み系列ｗが、ｗ１，ｗ２，…，ｗＬである場合には、下式（８）のようになる。また、図６は、２次元の複素配列の場合の配列への値の設定例を示した図である。
【００５５】
【数８】
2014022972.tif 000010
【００５６】
（ステップＳ１０２−３）
埋め込みパターン生成部２０は、ステップＳ１０２−２において生成した複素配列を埋め込みパターンＶとして出力する。
なお、複素配列への埋め込み系列の値の設定は、上記の式（８）や図６で示した構成に限られるものではなく、複素配列の中で埋め込み系列の値を設定する要素を選択的に行うようにしてもよい。例えば、複素配列の中で値を設定する要素部位をあらかじめ選択しておき、その要素部位に対して埋め込み系列の値を順に並べて設定するようにし、それ以外の要素の値を０にするようにしてもよい。
【００５７】
また、埋め込み系列の値を複素配列へ設定する前、すなわちステップＳ１０２−２の前に、各埋め込み系列の順序を擬似乱数に基づいてランダムな順序に入れ替えるようにしてもよい。順序を入れ替えることにより、埋め込まれている情報の不正な解析や、書き換えなどの攻撃を困難にすることができるとともに、インタリーブ符号化としての効果や、透かし埋め込み済み映像Ｉに対する局所的な耐性の不均衡を防ぐという効果も奏する。この構成は、順序の入れ替えに使用する擬似乱数の種の値を電子透かし埋め込みの鍵として与え、同一の鍵を電子透かしの検出時に使用するなどの構成により実現することもできる。
【００５８】
また、埋め込み系列ｗの順序の入れ替えを行う代わりに、ステップＳ１０２−２で得られた複素配列、すなわち埋め込み系列ｗを設定した後の複素配列の要素を入れ替えるようにしてもよい。
また、上記のステップＳ１０２−３において、複素配列を空間方向、すなわちＮ−１次元の方向で周波数変換してから埋め込みパターンＶとして出力するようにしてもよい。このとき、ステップＳ１０２−２において埋め込み系列ｗの値を設定する要素を選択的に行った場合は、埋め込みに用いる空間方向での周波数を特定の周波数係数に限って、要素の選択を行うことになる。
【００５９】
（ステップＳ１０３の処理：時間変調）
次に、時間変調部３０は、埋め込みパターン生成部２０が生成した埋め込みパターンＶを時間軸方向で変調して透かし映像パターンＷを生成する（ステップＳ１０３）。時間変調の処理は、例えば、以下に説明する特許文献２に記載されている時間変調部の処理を適用することができる。
【００６０】
（１）時間変調部３０が、同一の基本周波数を有する直交する２つの周期信号を生成する。例えば、１つの周期信号を基にそれぞれ位相が９０°、すなわち１／４周期分異なるように２つの周期信号を生成する。
（２）時間変調部３０が、入力した埋め込みパターンＶの実部、虚部を（１）で生成した２つの周期信号によってそれぞれ時間方向に変調する。具体的には、埋め込みパターンＶの位置ごとの値の実部、虚部の値を、（１）で生成した周期信号を搬送波として、ＡＭ（Amplitude Modulation）変調してＮ次元のパターンに変換する。
（３）時間変調部３０が、変調した２つの信号を加算して透かし映像パターンＷを生成する。
なお、上記（１）で２つの周期信号を生成するのに用いる周期信号は、１周期分積分した結果が０となり、自己相関関数が鋭敏なピークを示さないものであればよく、具体的には、正弦波、三角波、矩形波などを適用することができる。
【００６１】
以下、Ｎ＝３の映像信号、すなわち３次元の入力信号に対し周期信号として正弦波ｆｒ（ｔ）、ｆｉ（ｔ）を用いて時間変調を行い、透かし映像パターンＷを生成する手順を説明する。正弦波ｆｒ（ｔ）、ｆｉ（ｔ）を次式（９）とする。
【００６２】
【数９】
2014022972.tif 000011
【００６３】
このとき、透かし映像パターンＷは、式（１０）として表される。
【００６４】
【数１０】
2014022972.tif 000012
【００６５】
ただし、式（１０）においてωは周期Ｔに対応する角速度で、ω＝２π／Ｔである。なお、式（１０）における記号（１１）と記号（１２）は、それぞれＣの実部、虚部を取り出す演算を示す。
【００６６】
【数１１】
2014022972.tif 000013
【００６７】
【数１２】
2014022972.tif 000014
【００６８】
なお、これらの処理は、２つの周期信号を表す関数として複素関数ｆ（ｔ）＝ｅｊωｔ＝ｃｏｓωｔ＋ｊｓｉｎωｔを用いて式（１３）のように行ってもよい。式（１３）において、ｆ（＊）（ｔ）は、ｆ（ｔ）の共役複素数である。
【００６９】
【数１３】
2014022972.tif 000015
【００７０】
また、他の時間変調の例として、例えば、特許文献２に示されている以下の２つ手法を適用することができる。
（手法１）手法１は、埋め込みパターンＶの値の偏角を用いて、１つの周期関数の位相を変化させるように、例えば、式（１４）のような演算を行う手法である。
【００７１】
【数１４】
2014022972.tif 000016
【００７２】
式（１４）において、ａｒｇ（Ｃ）は、Ｃの偏角の値を算出する演算を表す。この場合、埋め込みパターンＶを複素数としてではなく、複素数の偏角を表す実数値として記憶しておくようにしてもよい。
【００７３】
（手法２）手法２は、埋め込みパターンＶの値を映像信号におけるＮ次元目の軸（例えば、時間軸）の特定の周波数に対するフーリエ係数であると見なして、離散逆フーリエ変換して透かし映像パターンＷを生成する手法である。
【００７４】
この手法により、透かし映像パターンＷの生成を行うことで、埋め込みパターンＶ上の位置（例えば、空間方向の位置）によってＮ次元目の方向（例えば、時間方向）の位相がそれぞれ異なるように透かし映像パターンＷを生成することができる。
【００７５】
なお、埋め込みパターン生成部２０による埋め込みパターン生成ステップＳ１０２の処理と、時間変調部３０による時間変調ステップＳ１０３の処理は、一体として、埋め込み系列ｗに基づいて透かし映像パターンＷを生成するという透かし映像パターン生成の処理とみなすことができる。
【００７６】
（ステップＳ１０４の処理：パターンの重畳）
パターン重畳部４０は、時間変調部３０が生成したＮ次元の透かし映像パターンＷとＮ次元信号の透かし埋め込み対象映像Ｉとを入力する。そして、パターン重畳部４０は、Ｎ次元信号の透かし埋め込み対象映像ＩにＮ次元の透かし映像パターンＷを加算することで重畳し、重畳した結果のＮ次元信号を透かし埋め込み済み映像Ｉ’として出力する（ステップＳ１０４）。
【００７７】
パターン重畳部４０による処理は、例えば、特許文献２に記載の埋め込みパターン重畳部の構成を適用することができる。例えば、Ｎ次元の透かし埋め込み対象映像Ｉを、Ｉ（ｘ１，ｘ２，…，ｘＮ−１，ｔ）とし、透かし映像パターンＷを、Ｗ（ｘ１，ｘ２，…，ｘＮ−１，ｔ）とした場合、透かし埋め込み済み映像Ｉ’（ｘ１，ｘ２，…，ｘＮ−１，ｔ）は、次式（１５）として求められる。
【００７８】
【数１５】
2014022972.tif 000017
【００７９】
ここで、αは、あらかじめ定められた透かし埋め込みの強度パラメータである。パターン重畳部４０は、算出した透かし埋め込み済み映像Ｉ’を出力する（ステップＳ１０４）。
なお、埋め込みパラメータαは、人間による電子透かしの知覚されやすさに応じて変化するパラメータとして算出されていてもよい。
【００８０】
また、埋め込みを行う対象とする信号成分は、透かし埋め込み対象映像Ｉの１つの成分であってもよいし、複数の成分であってもよい。例えば、透かし埋め込み対象映像Ｉの複数の色成分の内の１つであってもよいし複数の成分であってもよい。具体的には、Ｒ，Ｇ，Ｂの中の１つまたは複数の成分値であってもよいし、Ｙ，Ｃｂ，Ｃｒの成分値であってもよく、この他の信号成分の何れかであってもよい。
【００８１】
また、上記の第１実施形態の構成において、透かし埋め込み対象映像Ｉの大きさが、透かし映像パターンＷの大きさよりも大きい場合、透かし映像パターンＷを、時間方向や空間方向に繰り返して加算するようにしてもよい。
【００８２】
また、上記の第１実施形態の構成において、ステップＳ１０４のパターン重畳部４０による処理の前に、透かし映像パターンＷを任意の大きさ、もしくは透かし埋め込み対象映像Ｉの大きさに一致させるように拡大しておいてもよい。当該拡大は、空間方向に拡大してもよいし、時間方向に拡大してもよく、どのようなアルゴリズムを適用してもよい。例えば、拡大した１ブロックに１つの値を対応させて、拡大後のブロックを全て対応する１つの値と同じ値にする手法や、線形補間やバイキュービックなどの公知の補間手法を適用することができる。また、重畳される信号値が量子化によって失われてしまうのを防ぐため、量子化されないような値の大きさを保つよう、空間方向で誤差拡散やディザなどの処理を行うように拡大を行ってもよい。
【００８３】
また、上記の第１実施形態の構成において、パターン重畳部４０による処理において、透かし埋め込み対象映像Ｉの信号値に直接透かし映像パターンＷの値を加算するのではなく、透かし埋め込み対象映像Ｉの所定の特徴量を、透かし映像パターンＷの値だけ変更する、あるいはそのスカラー倍だけ変更するように透かし埋め込み対象映像Ｉを変更することによって重畳を行うようにしてもよい。所定の特徴量としては、例えば、透かし埋め込み対象映像Ｉのブロックごとの信号値の平均値や周波数変換した周波数係数などがあり、これら以外に、映像、画像の画素の輝度値や色差、ＲＧＢの色信号値などを適用することができる。
【００８４】
（その他の透かし重畳手法）
また、第１実施形態においては、埋め込みパターンＶを時間変調した透かし映像パターンＷを生成した上で、透かし埋め込み対象画像Ｉに透かし映像パターンＷを重畳する例を示したが、当該構成に限られず、例えば、特許文献２に示されている以下の手順による手法を用いてもよい。
【００８５】
図７は、当該手順を行う電子透かし埋め込み装置１’の構成を示した概略ブロック図である。電子透かし埋め込み装置１’は、電子透かし埋め込み装置から時間変調部３０を除いた構成となっており、以下、異なる構成であるパターン重畳部４０’について説明する。
【００８６】
（１）パターン重畳部４０’は、入力する透かし埋め込み対象映像Ｉをいったん時間方向で周波数変換して変換済み信号を求める。
（２）パターン重畳部４０’は、埋め込みパターンＶを、変換済み信号の特定周波数係数に重畳して逆変換前信号を求める。すなわち、変換済み信号の一部のＮ−１次元平面に埋め込みパターンＶを重畳する。
（３）パターン重畳部４０’は、逆変換前信号を逆周波数変換して、透かし埋め込み済み映像Ｉ’を求める。ここで、周波数変換としては、直交変換を行う演算であればよく、例えば離散フーリエ変換が適用できる。
【００８７】
例えば、Ｎ＝３である透かし埋め込み対象映像Ｉを時間方向で周波数変換した変換済み信号に対して、時間方向周波数の１．０Ｈｚに対応する周波数係数を取り出すと、２次元の空間位置毎の周波数係数値が並んだ２次元配列が得られる。これに対し、埋め込みパターンＶを空間位置毎に加算することで逆変換前信号を得る。これを逆周波数変換して透かし埋め込み済み映像Ｉ’を求めることで、時間変調部３０において周期が１．０秒である時間方向周期信号を用いて時間変調を行った後、パターン重畳部４０により透かし映像パターンＷの重畳を行うのと同様の透かし埋め込み済み映像Ｉ’を得ることができる。
【００８８】
なお、時間変調部３０で用いる周期信号や、パターン重畳部４０’で用いる特定周波数は、透かし埋め込み済み映像Ｉ’における透かし信号の知覚されやすさを考慮して決定してもよい。すなわち、人の視覚特性により、知覚されやすい（知覚感度の高い）時間方向周波数と知覚されにくい（知覚感度の低い）時間方向周波数があることから、知覚感度の低い時間方向周波数に対して重畳を行うことで、視覚的に目立たない電子透かしを埋め込むことができる。
【００８９】
上記のように、図７に示す電子透かし埋め込み装置１’は、透かし埋め込み対象映像Ｉを周波数変換した変換済み信号に対して空間位置を指定して埋め込みパターンＶを埋め込むことができるため、例えば、知覚感度の低い時間方向周波数の選択を意識した埋め込み等を容易に行うことが可能となる。一方、図１に示す電子透かし埋め込み装置１は、透かし埋め込み対象映像Ｉを周波数変換する必要がないため、電子透かし埋め込み装置１’に比べて少ない演算量で電子透かしの埋め込みを行うことが可能となる。
【００９０】
また、パターン重畳部４０’が、透かし埋め込み対象映像Ｉに対する周波数の変換を行うのではなく、パターン重畳部４０’とは別に周波数変換部を電子透かし埋め込み装置１’が備え、当該周波数変換部が透かし埋め込み対象映像Ｉを周波数変換してパターン重畳部４０’に出力する構成であってもよい。
【００９１】
なお、１シンボルで表現できる複素数の数をＫとした場合、上記の実施形態の記載における例の関係は、それぞれ以下の通りの構成となる。すなわち、
Ｋ＝３：図５に示す３進数の１桁を１シンボルに割り当てる例、
Ｋ＝４：図４（ａ）に示す２ビットの１シンボルに割り当てる例、
Ｋ＝８：図４（ｂ）に示す３ビットを１シンボルに割り当てる例、
Ｋ＝２ｇ：ｇビットずつを１ビットに割り当てる例、
Ｋ＝ｔｌ：ｔ進数のｌ桁を１シンボルに割り当てる例、
となる。このように、３以上の種類の複素数を１シンボルとして割り当てて複素情報系列を生成することで、特許文献２に記載の技術に比べてより多くのビット値を１シンボルで表現することが可能となる。
【００９２】
上述した第１実施形態の電子透かし埋め込み装置１の構成により、埋め込み情報Ｕから複素情報系列ｓを生成し、生成した複素情報系列ｓと複素数の拡散系列との要素ごとの積を用いて埋め込み系列ｗを生成する構成とした。これにより、ｍ個の複素数の拡散系列に対して従来よりも多くのビット値を与えることができ、拡散系列の長さが同じ、すなわち複素数として埋め込みのできる位置の個数が同じ場合には、全体として多くの情報を埋め込むことが可能となる。また、拡散系列の全体の長さを同じとした場合、１シンボルに割り当てる拡散系列の長さｍを２倍にし、一方で１シンボルで表すビット値を２倍にすれば、ｍが大きくなると信頼性が増すことから、全体として従来と同じ情報量を信頼性高く埋め込むことが可能となる。すなわち、従来と同じ大きさの複素配列に対して、従来よりも多くの情報を埋め込むことができ、具体的には、映像信号に対して埋め込みの情報量を増加させることが可能となる。また、埋め込みの情報量を増加させずに維持した場合には、信頼性高く埋め込むことが可能となる。
【００９３】
（第１実施形態の電子透かし検出装置）
次に、第１実施形態による電子透かし検出装置２について図面を参照しつつ説明する。
図８は、第１実施形態による電子透かし検出装置２の概略ブロック図である。
電子透かし検出装置２は、時間復調部５０と、検出系列抽出部６０と、相関値計算部７０と、透かし判定部８０とを備え、透かし埋め込み済み映像Ｈを入力し、透かし検出結果Ｒを出力する。
【００９４】
時間復調部５０は、透かし埋め込み済み映像Ｈを入力し、透かし埋め込み済み映像Ｈに基づいて生成した複素配列Ｑを出力する。なお、透かし埋め込み済み映像Ｈは、電子透かし埋め込み装置１が生成した透かし埋め込み済み映像Ｉ’と同じ映像である場合の他、電子透かしが埋め込まれていない映像の場合も含む。検出系列抽出部６０は、時間復調部５０が生成した複素配列Ｑを入力し、複素配列Ｑに基づいて生成した検出系列ｄを出力する。相関値計算部７０は、検出系列抽出部６０が生成した検出系列ｄを入力し、検出系列ｄに基づいて生成した相関値ρを出力する。透かし判定部８０は、相関値ρを入力し、相関値ρに基づいて透かしの有無の判定を行い、シンボルを検出して透かし検出結果Ｒを出力する。
【００９５】
なお、透かし埋め込み済み映像Ｈの入力は、外部の記憶装置に記憶されているものを読み出して入力するか、あるいは外部の装置が出力した情報を入力するなどのように構成することができる。また、複素配列Ｑ、検出系列ｄ、相関値ρは、各機能部の間で直接入出力する構成の他、メモリ等の記憶手段を介するような構成であってもよい。すなわち、それぞれの機能部によって、これらの情報が生成されたときに、生成の処理を行った機能部が電子透かし検出装置２に備えられるメモリなどの記憶手段に生成した情報を書き込み、次段の機能部がそれを読み出すことによって処理されるようにしてもよい。
【００９６】
次に、図９を参照しつつ、電子透かしの検出の処理について説明する。
まず、最初に、時間復調部５０は、透かし埋め込み済み映像Ｈを入力し、この透かし埋め込み済み映像Ｈに対して時間復調の処理を行い複素配列Ｑを生成する（ステップＳ２０１）
時間復調部５０は、例えば、特許文献２に記載されている時間復調部の構成を適用することができ、以下のような手順で処理を行う。
【００９７】
（１）時間復調部５０は、同一の基本周波数を有する直交する２つの周期信号を生成する。例えば、１つの周期信号を基にそれぞれ位相が９０°、すなわち１／４周期分異なるように２つの周期信号を生成する。なお、周期信号については上述した時間変調部３０で用いたものを適用することができる。
（２）時間復調部５０は、入力した透かし埋め込み済み映像Ｈの時間方向成分を元に、（１）で生成した２つの周期信号のそれぞれを用いて復調し、２つのＮ−１次元信号を算出する。
（３）時間復調部５０は、（２）で算出した２つのＮ−１次元信号が、それぞれ実部、虚部となる複素数のＮ−１次元配列である複素配列Ｑを生成する。
【００９８】
以下に、Ｎ＝３の映像信号に対し周期信号として正弦波ｆｒ（ｔ），ｆｉ（ｔ）を用いて時間復調を行い、複素配列Ｑを生成する処理を説明する。まず、透かし埋め込み済み映像ＨがＨ（ｘ１，ｘ２，…，ｘＮ−１，ｔ）で表されているとする。正弦波ｆｒ（ｔ），ｆｉ（ｔ）を式（１６）とする。
【００９９】
【数１６】
2014022972.tif 000018
【０１００】
このとき、複素配列Ｑ（ｘ１，ｘ２，…，ｘＮ−１）は、次式（１７）によって算出される。なお、式（１７）においてｊは虚数単位である。
【０１０１】
【数１７】
2014022972.tif 000019
【０１０２】
なお、上記の構成では、積分演算に基づいて複素配列Ｑを算出する例を示したが、透かし埋め込み済み映像Ｈが離散信号として得られている場合は、積和演算により複素配列Ｑを算出することもできる。また、透かし埋め込み済み映像Ｈが連続信号として得られている場合であっても任意の標本化手段を用いて離散信号に変換してから処理することもできる。
【０１０３】
また、この他の時間復調の例として、例えば、特許文献２に記載されている手法を適用することができる。すなわち、時間復調部５０が、入力した透かし埋め込み済み映像Ｈから所定の区間を取りだし、Ｎ次元目の軸（例えば、時間方向）に一次元離散フーリエ変換し、所定の周波数係数を算出して複素配列Ｑとする手法である。
【０１０４】
次に、検出系列抽出部６０は、時間復調部５０が生成した複素配列Ｑから抽出して、それぞれ複素数の系列である検出系列ｄを生成する（ステップＳ２０２）。
（１）検出系列抽出部６０は、入力した複素配列Ｑから大きさＸ１×Ｘ２×…×ＸＮ−１のＮ−１次元複素配列を生成する。ただし、Ｘ１，Ｘ２，…，ＸＮ−１は、上述した電子透かし埋め込み装置１の埋め込みパターン生成部２０が埋め込みパターンＶの生成のために確保した配列と同じ要素数である。
【０１０５】
透かし埋め込み済み映像Ｈの大きさがＸ１×Ｘ２×…×ＸＮ−１と等しい場合は、複素配列Ｑを同じ大きさのＸ１×Ｘ２×…×ＸＮ−１のＮ−１次元複素配列とすることができる。また、透かし埋め込み済み映像Ｈの大きさがＸ１×Ｘ２×…×ＸＮ−１と異なる場合は、複素配列Ｑの大きさもＸ１×Ｘ２×…×ＸＮ−１とは異なってくる場合がある。その場合は、拡大、縮小などの任意の標本化手段を用いて標本化して大きさをＸ１×Ｘ２×…×ＸＮ−１に一致させてＮ−１次元の複素配列とすることができる。
【０１０６】
具体的には、複素配列ＱをＸ１×Ｘ２×…×ＸＮ−１のブロック数の領域に分割し、各ブロック領域に含まれる値の平均値を求めるようにして縮小する手法などを適用することができる。また、線形補間やバイキュービックなどの公知の補完手法を用いて拡大、縮小を行ってもよい。
【０１０７】
（２）検出系列抽出部６０は、手順（１）で生成したＮ−１次元複素配列Ｑから、順に１つずつ複素数値を読み出し、読み出した複素数値を並べて検出系列ｄを生成する。例えば、Ｎ−１次元複素配列Ｑが、Ｂ［ｐ１，ｐ２，…ｐｎ，…，ｐＮ−１］（ｐｎ≧０）、検出系列ｄが、ｄ１，ｄ２，…，ｄＬと表されているとき、式（１８）のように抽出されることになる。
【０１０８】
【数１８】
2014022972.tif 000020
【０１０９】
また、上述した埋め込みパターン生成部２０において、各埋め込み系列ｗの順序を擬似乱数を用いて、ランダムな順序に入れ替えて埋め込みが行われている場合、上記の手順（２）の後に、擬似乱数を用いて検出系列ｄの順序を埋め込みパターン生成部２０が入れ替えた場合とは逆の順序で入れ替えて、埋め込み系列ｗの順序に対応付けられる順序に戻しておく必要がある。このとき、順序の入れ替えに使用する擬似乱数の種の値は、電子透かし検出の鍵として、電子透かしの埋め込み時に用いたのと同一の鍵を与えることになる。
【０１１０】
次に、相関値計算部７０は、検出系列抽出部６０が生成した検出系列ｄを入力し、検出系列ｄと埋め込みの際に用いた擬似乱数列との相関値を、ｋ個のシンボルの各々に対して算出する（ステップＳ２０３）。具体的には以下の手順で処理を行う。
【０１１１】
（１）相関値計算部７０は、電子透かし埋め込み装置１において埋め込み系列ｗの生成の処理に用いたのと同じ複素数の擬似乱数列ＰＮ＝｛ＰＮ１，ＰＮ２，…，ＰＮｍｋ｝を用いて処理を行う。当該擬似乱数系列ＰＮの長さは、Ｌ＝ｍｋであり、電子透かし埋め込み装置１の場合と同じく、処理の際に演算によって生成する構成の他、予め演算によって生成しておき、内部のメモリ等の記憶手段に予め記憶させておいたものを読み出すことによって利用する構成などを適用することができる。
（２）相関値計算部７０は、検出系列ｄをｍ個ずつに分割する。このとき、ｊ番目の部分系列をｄ（ｊ）で表す。
（３）相関値計算部７０は、ｄ（ｊ）と擬似乱数列ＰＮの対応する部分系列ＰＮ（ｊ）との複素相関を算出する。すなわち、j 番目の部分に対する相関値をρ（ｊ）とすると、式（１９）のように算出することができる。
【０１１２】
【数１９】
2014022972.tif 000021
【０１１３】
ただし、式（１９）において「・」は数列をベクトルと見たときの内積演算を表し、Ｃ＊は複素数Ｃの共役複素数を表す。
また、後述する参考文献１や３に記載の検出信頼性の評価基準にそろえるため、例えばｄ（ｊ）およびＰＮ（ｊ）の各要素をあらかじめ平均０、実部、虚部の分散が１となるように正規化した、式（２０）で示される検出系列の部分系列、式（２１）で示される擬似乱数列の部分系列の値を求めるようにしてもよい。
【０１１４】
【数２０】
2014022972.tif 000022
【０１１５】
【数２１】
2014022972.tif 000023
【０１１６】
これらの正規化した検出系列及び擬似乱数列を用いて、式（２２）のように相関値の算出の際、定数項を乗じて演算を行うことにより評価基準がそろった相関値を算出することができる。
【０１１７】
【数２２】
2014022972.tif 000024
【０１１８】
次に、透かし判定部８０は、相関値計算部７０が算出した相関値ρを入力し、相関値ρに基づいて判定を行う（ステップＳ２０４）。電子透かし検出装置２が入力する透かし埋め込み済み映像Ｈに電子透かしが埋め込まれている場合と、埋め込まれていない場合では、相関値計算部７０が算出するｋ個の相関値ρ（ｋ）の絶対値に違いが現れる。すなわち、電子透かしが埋め込まれている場合、ｋ個の相関値ρ（ｋ）の絶対値のそれぞれの値が大きな値となり、逆に電子透かしが埋め込まれていない場合には、k個の相関値ρ（ｋ）の絶対値のそれぞれの値が小さな値となる。
【０１１９】
したがって、予め定められる所定の閾値を用いて相関値ρ（ｋ）の絶対値が閾値を超えているかどうかを判定することで、透かしが埋め込まれているか否かを判定することができる。なお、ｋ個の相関値ρ（ｋ）のそれぞれの絶対値で判定するのではなく、ｋ個のシンボル全体で判定を行ってもよく、その場合、式（２３）で示されるk個の相関値ρ（ｋ）の総和の絶対値を用いて判定を行うようにしてもよい。
【０１２０】
【数２３】
2014022972.tif 000025
【０１２１】
また、式（２４）で示されるk個の相関値ρ（ｋ）の絶対値の総和を用いて判定を行うようにしてもよい。
【０１２２】
【数２４】
2014022972.tif 000026
【０１２３】
なお、上述したように検出信頼性の評価基準をそろえるように相関値を正規化して算出していた場合、求められる誤検出率の上限値に基づいて検出信頼性を評価して閾値を定めておくこともできる。また、特許文献２に示される電子透かし検出装置と同様に、相関値の偏角を用いて時間方向の同期を行うようにしてもよい。
【０１２４】
上記の判定によって、透かし判定部８０が、電子透かしが埋め込まれていると判定した場合、透かし判定部８０は、次に、相関値の位相に基づいて、埋め込まれているシンボルの値を判定し、判定結果のシンボル値を並べて透かし検出結果Ｒを算出する。以下、シンボルの総数ｋ＝１の場合を例として透かし判定部８０によるシンボルの値の判定の具体的方法について説明する。
埋め込み系列ｗの１要素ｗｉについて取り出して考えると、前述の埋め込み系列生成部１０による埋め込み系列の生成と同じくｗｉは、式（２５）として表すことができる。
【０１２５】
【数２５】
2014022972.tif 000027
【０１２６】
これを元に生成された透かし映像パターンＷを、パターン重畳を行うステップＳ１０４で入力信号に重畳することは、式（２６）で示されるように埋め込み系列ｗｉにノイズ成分Ｎｉを加えた透かし入り信号ｆｉを得ることと等価である。
【０１２７】
【数２６】
2014022972.tif 000028
【０１２８】
また、透かし埋め込み済み映像Ｈに対し、Ｎ次元目の方向（すなわち時間方向）で異なるタイミングから検出を開始した場合は、これはｆｉの位相がずれることを意味するため、その結果、検出系列の要素ｄｉは、位相の変位量をφとすると式（２７）となる。
【０１２９】
【数２７】
2014022972.tif 000029
【０１３０】
上述したように、相関値計算部７０によるステップＳ２０３の処理では、式（２８）のように検出系列ｄと擬似乱数列ＰＮの複素相関を算出している。
【０１３１】
【数２８】
2014022972.tif 000030
【０１３２】
式（２８）を、一要素のみについてみると式（２９）として表すことができる。
【０１３３】
【数２９】
2014022972.tif 000031
【０１３４】
式（２９）の右辺の最後の式の第一項はシンボル値ｓｉを位相変位量φだけ複素平面状で回転させた値の実数倍となる。一方、擬似乱数列ＰＮが平均０であるように生成されていれば、第二項は平均０のランダムな値となり、式（２９）の値は、図１０（ａ）のような値の分布を取ることになる。
【０１３５】
相関値ρの算出において、式（２９）の各要素の総和をとると、第一項は共通の偏角を有するのに対して、第二項はランダムな偏角となり、結果として図１０（ｂ）のような分布となる。
さらに、特許文献２や後述する参考文献３に記載の電子透かし検出装置と同様に、相関値ρの偏角を用いて時間方向の同期を行えば、位相位置をキャンセルした位置に相関値ρを変換することができる。
【０１３６】
このように変換した相関値ρに基づいて、透かし判定部８０は、シンボルの判定を行う。例えば、複素情報系列ｓの取りうる値を図４（ａ）のような４値を取るようにしていた場合、図１１（ａ）に示す４つの領域のいずれかに分布することとなる。図１１（ａ）のように実軸、虚軸をシンボル値境界として設定して、どの領域に相関値ρが属しているかを判定することで、どのシンボル値が埋め込まれていたのかを判定することができる。
【０１３７】
また、複素情報系列ｓの取りうる値を図４（ｂ）のような８値を取るようにしていた場合、図１１（ｂ）のような８つの領域のいずれかに分布することとなる。図１１（ｂ）のような破線で示すシンボル値境界を設定し、相関値ρの位相を判定することで、どのシンボル値が埋め込まれていたのかを判定することができる。
【０１３８】
次に、シンボルの総数がｋ＞１の場合、上述の相関値ρはシンボル毎に求められ、シンボルｓｊに対して相関値ρ（ｊ）が求められる。ρ（ｊ）の分布は図１２（ａ）のように、あらかじめ埋め込まれていたシンボルの値に応じて複数の分布のいずれかを取ることになる。電子透かしの埋め込み時と検出時で、Ｎ次元目の方向（すなわち時間方向）の同期がずれていた場合、図１２（ｂ）のように位相変位φだけ回転していることになる。
【０１３９】
本実施形態の電子透かし検出装置２の透かし判定部８０は、１つのシンボルに対する相関値の位相に基づいて多種のシンボル値のいずれであるかを判定する。そのため、シンボル値の種類の数が増加する場合、誤った判定を防ぐため、より高い精度での位相同期を行う必要が生じる。位相の同期は、例えば、次のように行われる。相関値計算ステップで得られた各ρ（ｊ）に対し、以下の（ｉ）、（ｉｉ）、（ｉｉｉ）のルールを適用する。
【０１４０】
【数３０】
2014022972.tif 000032
【０１４１】
これにより、全てのρ（ｊ）を同じ方向を向けることができ、これらの総和である＾ρを算出する。ρ（ｊ）を各々−ａｒｇ（＾ρ）＋π／４の値だけ位相を変位させることで、位相の同期を行うことができる。ここで、ａｒｇ（Ｃ）、またはａｒｇＣは、複素数Ｃの偏角を求める演算を表す。
さらに、以下のように行うことで、より精度の高い位相同期を行うことができる。
【０１４２】
（１）上記の手順により総和である＾ρを求める。
（２）ρ（ｊ）に対して以下の（ｉ）、（ｉｉ）、（ｉｉｉ）のルールを適用し、適用した後に総和を求め＾ρ’とする。
【０１４３】
【数３１】
2014022972.tif 000033
【０１４４】
（３）ρ（ｊ）を各々−ａｒｇ（＾ρ）＋π／４の値だけ位相を変位させる。
これにより、より精度の高い位相の同期をとることができる。なお、手順（２）で得られた＾ρ’を＾ρに置き換え、手順（２）を繰り返すようにしてもよい。
【０１４５】
すなわち、（１）固定の位相の境界値を用いて相関値を判定して位相変位させた上で総和を求め、（２）総和の位相に基づき新たに得られる境界値を用いて相関値を判定して位相変位させた上で再度総和を求める、という手順を行うことでより精度の高い位相同期が可能となり、検出性能が向上する効果が得られる。
【０１４６】
上記の第１実施形態の構成では、４種のシンボル値を一例として説明したが、当該構成に限られず、この他の種類のシンボル値を取る場合でも同様の方法を用いることができる。具体的には、例えば、ｑ種のシンボル値を取る場合には、上記の手順（２）で、ρ（ｊ）について次式（３２）を満たすものについてはｅ−２πρ／ｑを乗じるようにすればよい。
【０１４７】
【数３２】
2014022972.tif 000034
【０１４８】
なお、シンボル値の誤判定に対する耐性を持たせるため、埋め込み情報Ｕを、あらかじめ誤り訂正符号を用いて符号化して用い、透かし検出結果Ｒを得た後で誤り訂正復号を行うようにしてもよい。例えば、シンボル値がｑ種の値である場合に、ｑ元のＢＣＨ（Bose-Chaudhuri-Hocquenghem）符号やリード・ソロモン符号を用いるようにしてもよい。バースト誤り耐性の高いｑ元の誤り訂正符号を用いることにより、電子透かしに対して付加されるノイズの特性によって発生する可能性があるシンボル単位での誤判定を軽減することが可能となる。また、相関値ρの絶対値に対する閾値を設け、閾値よりも小さい相関値ρの場合には、そのシンボルについては値が得られなかったとして処理し、例えば消失のある通信路を経た場合の誤り訂正符号等の方式を用いて誤り訂正を行うようにしてもよい。
【０１４９】
（第１実施形態の検出の信頼性について）
第１実施形態では、埋め込み系列ｗがランダムに生成された複素数列となっている。したがって、特許文献２や参考文献２に示されているように、埋め込む複素数列の分散が等しければ検出の信頼性は共通になることが導ける。
特許文献２や参考文献２に記載の技術を用いた参考文献３に記載の技術では、複素数の複素情報系列ｓを用いず、実数の１と−１の系列を用いて擬似乱数列を埋め込み系列に変換しているため、いわば、１つのシンボルとして１ビットの情報しか表していない。一方、本実施形態の構成では、複素数の複素情報系列ｓを用いて埋め込み系列ｗに変換している。
【０１５０】
したがって、透かしの埋め込まれていない信号から誤って透かしを検出してしまう確率として表される検出の信頼性を特許文献２や参考文献２に記載の技術と同じ小さな値に保ちつつ、Ｎ−１次元空間上で埋め込みに用いる位置を増やすことなく、１つのシンボルに対して複数ビットの情報を表現して電子透かしの埋め込みを行うことができる。
【０１５１】
すなわち、従来と同じ大きさの複素配列に対して、検出の信頼性を維持しつつ、従来よりも埋め込みの情報量を増加させることができる。例えば、図４（ａ）のように１つのシンボルに２ビットの情報を割り当てれば、参考文献３と比較して２倍の情報を埋め込むことが可能となる。また、逆に、従来と同じ大きさの複素配列に対して、従来と同じ情報量を埋め込む場合には、１つのシンボルあたりに対応する埋め込み位置の数、すなわち対応する拡散系列の長さを長くすることができるため、検出の信頼性を高めることができる。
【０１５２】
参考文献１：中村高雄、片山淳、山室雅司、曽根原登、「カメラ付き携帯電話機を用いたアナログ画像からの高速電子透かし検出方式」、電子通信学会和文論文誌D-II、2004、Vol. J87-D-II、 No. 12、pp. 2145―2155
【０１５３】
参考文献２：山本奏、中村高雄、片山淳、安野貴之、「単一周波数平面スペクトル拡散を利用した時間同期外し耐性を持つ動画電子透かし」，電子通信学会和文論文誌D、2007、Vol. J90-D, No.7, pp. 1755-1764
【０１５４】
参考文献３：中村高雄、山本奏、北原亮、片山淳、安野貴之、小池秀樹、曽根原登、「SFPSS 法に基づくリアルタイム検出可能な映像向けモバイル電子透かし」、情報処理学会論文誌、2008、Vol.49、No.6、pp.1885-1895
【０１５５】
上記の第１実施形態の構成により、複素数の複素情報系列を用いて、複素情報系列と擬似乱数系列の要素ごとの積を算出することで、擬似乱数系列の要素の値の位相を、複素情報系列の要素の値の位相で変位させた、埋め込み系列を生成することができる。これにより、検出の信頼性を保ったまま、映像信号に対して埋め込みの情報量を増加させた電子透かしの埋め込みが可能となる。また、埋め込みの情報量を増加させずに維持した場合には、検出の信頼性を高めることができるという効果が得られる。
【０１５６】
（第２実施形態）
次に、本発明の第２実施形態について図面を参照しつつ説明する。第１実施形態では、電子透かしを映像に埋め込む構成について示したが、第２実施形態では、電子透かしを静止画像に対して埋め込む構成を示す。
【０１５７】
図１３は、第２実施形態の電子透かし埋め込み装置１ａの概略ブロック図である。電子透かし埋め込み装置１ａは、埋め込み系列生成部１０ａ、埋め込みパターン生成部２０ａ、空間変調部３０ａ、パターン重畳部４０ａを備える。
次に、第２実施形態における電子透かしの埋め込みの処理について説明する。第２実施形態の電子透かし埋め込み装置１ａの構成は、以下の点で第１実施形態と異なり、それ以外の符号から「ａ」を除いた符号が第１実施形態の機能部の符号と同じ機能部の構成については第１実施形態と同様である。以下、異なる構成について説明する。
【０１５８】
埋め込みパターン生成部２０ａは、第１実施形態の埋め込みパターン生成部２０の構成においてＮ＝２とした構成となる。また、パターン重畳部４０ａについても、第１実施形態のパターン重畳部４０においてＮ＝２次元の画像を対象にパターン重畳を行う構成となる。
【０１５９】
空間変調部３０ａは、第１実施形態の時間変調部３０と異なり、埋め込みパターンＶ１に対して時間方向で変調する代わりに空間方向で変調する。空間方向での変調は、例えば、X 軸方向に並べたＮ−１＝１次元の埋め込みパターンＶ１を、Ｙ軸方向に周期信号を用いて式（３３）のように変調する。
【０１６０】
【数３３】
2014022972.tif 000035
【０１６１】
このように変調することにより、Ｎ−１＝１次元の埋め込みパターンＶ１上の位置（例えば、Ｘ軸方向の位置）によって、Ｎ次元目の方向（例えば、Ｙ軸方向）の位相がそれぞれ異なるように式（３３）の左辺のＷ（ｘ，ｙ）、すなわち透かし画像パターンＷ１を生成することができる。
【０１６２】
なお、空間方向での変調は、上記の構成に限られず、この他の変調方法であってもかまわない。例えば、Ｙ軸方向に並べた埋め込みパターンＶ１を、Ｘ軸方向に周期信号を用いて変調してもよいし、斜めの方向に並べた埋め込みパターンＶ１を別な斜めの方向に変調するようにしてもよい。斜め方向に変調した場合は、検出における空間復調においてその方向で復調を行うことになる。
【０１６３】
なお、第１実施形態と同じく、透かし埋め込み対象画像Ｉ１を一度、空間方向で周波数変換し、その特定周波数係数に埋め込みパターンＶ１を重畳した上で、埋め込みパターンＶ１が重畳された信号を逆周波数変換して透かし埋め込み済み画像Ｉ１’を求めるようにしてもよい。
【０１６４】
図１４は、第２実施形態における電子透かし検出装置２ａの概略ブロック図である。
電子透かし検出装置２ａは、空間復調部５０ａと、検出系列抽出部６０ａと、相関値計算部７０ａと、透かし判定部８０ａとを備え、透かし埋め込み済み画像Ｈ１を入力し、透かし検出結果Ｒ１を出力する。なお、透かし埋め込み済み画像Ｈ１は、電子透かし埋め込み装置１ａが生成した透かし埋め込み済み画像Ｉ１’と同じ映像である場合の他、電子透かしが埋め込まれていない画像の場合も含む。
【０１６５】
次に、第２実施形態における電子透かしを検出する処理について説明する。第２実施形態の電子透かし検出装置２ａの構成は、以下の点で第１実施形態と異なり、それ以外の符号から「ａ」を除いた符号が第１実施形態の機能部の符号と同じ機能部の構成については第１実施形態と同様である。以下、異なる構成について説明する。
空間復調部５０ａは、第１実施形態の時間復調部５０と異なり、透かし埋め込み済み画像Ｈ１を時間方向で復調する代わりに空間方向で復調する。
【０１６６】
空間方向での復調は例えば、Ｎ＝２次元の透かし埋め込み済み画像Ｈ１を、Ｙ軸方向の周期信号を用いて復調し、Ｘ軸方向に並んだＮ−１＝1 次元の複素配列Ｑ１を次式（３４）の演算により生成する。なお、式（３４）においてＱ（ｘ）が、複素配列Ｑ１であり、Ｈ（ｘ，ｙ）が、透かし埋め込み済み画像Ｈ１である。
【０１６７】
【数３４】
2014022972.tif 000036
【０１６８】
空間復調部５０ａによる空間方向での復調は、空間変調部３０ａでの変調方法に対応する復調となる。例えば、空間変調部３０ａで、Ｙ軸方向に並べた埋め込みパターンＶ１を、Ｘ軸方向に周期信号を用いて変調している場合は、空間復調部５０ａは、Ｘ軸方向の周期信号を用いて復調してＹ軸方向に並んだ複素配列Ｑ１を生成する。また、空間変調部３０ａで、斜めの方向に並べた埋め込みパターンＶ１を別な斜めの方向に変調するようにしている場合は、空間復調部５０ａは、２つ目の斜め方向で復調を行う。
【０１６９】
空間変調部３０ａが、Ｙ軸方向に並べた埋め込みパターンＶ１を、Ｘ軸方向に周期信号を用いて変調している場合、Ｘ軸方向の変位が行われた場合、例えば、画像の切り取りが行われた場合、位相変位が発生する。また、この他の方向に周期信号を用いて変調している場合は、当該方向の変位が行われた場合に位相変位が発生する。これらの位相変位が生じている場合の同期は、第１実施形態の透かし判定部８０による位相同期と同じ方式で行うことができる。
【０１７０】
上記の第２実施形態の構成により、静止画像に対する電子透かし埋め込みを行う場合にも、複素数の複素情報系列を用いて、複素情報系列と擬似乱数系列の要素ごとの積を算出することで、擬似乱数系列の要素の値の位相を、複素情報系列の要素の値の位相で変位させた、埋め込み系列を生成することができる。これにより、検出の信頼性を保ったまま、静止画像に対して埋め込みの情報量を増加させた電子透かしの埋め込みが可能となる。また、埋め込みの情報量を増加させずに維持した場合には、検出の信頼性を高めることができるという効果が得られる。
【０１７１】
（第３実施形態）
次に、本発明の第３実施形態について図面を参照しつつ説明する。第３実施形態では、第２実施形態と同じく電子透かしを静止画像に対して埋め込む構成を示す。
図１５は、第３実施形態の電子透かし埋め込み装置１ｂの概略ブロック図である。電子透かし埋め込み装置１ｂは、埋め込み系列生成部１０ｂ、埋め込みパターン生成部２０ｂ、逆周波数変換部３０ｂ、パターン重畳部４０ｂを備える。
【０１７２】
次に、第３実施形態における電子透かしの埋め込みの処理について説明する。第３実施形態の電子透かし埋め込み装置１ｂの構成は、以下の点で第１実施形態と異なり、それ以外の符号から「ｂ」を除いた符号が第１実施形態の機能部の符号と同じ機能部の構成については第１実施形態と同様である。以下、異なる構成について説明する。
【０１７３】
埋め込みパターン生成部２０ｂは、逆周波数変換部３０ｂでの処理の結果として得られる透かし画像パターンＷ２が実数になるように、原点対象位置の値の組が複素共役になるように埋め込みパターンＶ２を生成する。なお、埋め込みパターンＶ２は、Ｎ−１次元ではなくＮ次元パターンである。
【０１７４】
逆周波数変換部３０ｂは、埋め込みパターン生成部２０ｂが生成した埋め込みパターンＶ２を逆周波数変換し、透かし画像パターンＷ２を生成する。逆周波数変換としては、埋め込みパターンＶ２を周波数係数とみて実空間に変換する直交変換であればよく、例えば、逆離散フーリエ変換を適用してもよい。パターン重畳部４０ｂは、Ｎ＝２次元の画像を対象にパターン重畳を行う構成のみが第１実施形態のパターン重畳部４０と異なる。
【０１７５】
なお、第１実施形態と同じく、透かし埋め込み対象画像Ｉ２を一度、空間方向で周波数変換し、その特定周波数係数に埋め込みパターンＶ２を重畳した上で、埋め込みパターンＶ２が重畳された信号を逆周波数変換して透かし埋め込み済み画像Ｉ２’を求めるようにしてもよい。
【０１７６】
図１６は、第３実施形態における電子透かし検出装置２ｂの概略ブロック図である。
電子透かし検出装置２ｂは、周波数変換部５０ｂと、検出系列抽出部６０ｂと、相関値計算部７０ｂと、透かし判定部８０ｂとを備え、透かし埋め込み済み画像Ｈ２を入力し、透かし検出結果Ｒ２を出力する。なお、透かし埋め込み済み画像Ｈ２は、電子透かし埋め込み装置１ｂが生成した透かし埋め込み済み画像Ｉ２’と同じ映像である場合の他、電子透かしが埋め込まれていない画像の場合も含む。
【０１７７】
次に、第３実施形態における電子透かしを検出する処理について説明する。第３実施形態の電子透かし検出装置２ｂの構成は、以下の点で第１実施形態と異なり、それ以外の符号から「ｂ」を除いた符号が第１実施形態の機能部の符号と同じ機能部の構成については第１実施形態と同様である。以下、異なる構成について説明する。
【０１７８】
周波数変換部５０ｂは、透かし埋め込み済み画像Ｈ２を周波数変換し、複素配列Ｑ２を生成する。周波数変換としては、透かし埋め込み済み画像Ｈ２を実空間値とみて周波数空間に変換する直交変換であればよく、例えば、離散フーリエ変換を適用してもよい。
なお、第３実施形態では、第１実施形態や第２実施形態のように共通の位相変位が発生しないため、透かし判定部８０ｂにおいて位相同期は行わない。
【０１７９】
上記の第３実施形態の構成により、静止画像に対する電子透かし埋め込みを行う場合にも、複素数の複素情報系列を用いて、複素情報系列と擬似乱数系列の要素ごとの積を算出する。そして、擬似乱数系列の要素の値の位相を、複素情報系列の要素の値の位相で変位させた、埋め込み系列を生成する。これにより、検出の信頼性を保ったまま、静止画像に対して埋め込みの情報量を増加させた電子透かしの埋め込みが可能となる。また、埋め込みの情報量を増加させずに維持した場合には、検出の信頼性を高めることができるという効果が得られる。
【０１８０】
（第４実施形態）
次に、本発明の第４実施形態について図面を参照しつつ説明する。第４実施形態は、第１、２、３実施形態の埋め込み系列生成部１０、１０ａ、１０ｂが行う埋め込み系列の演算手法とそれに対応する相関値計算部７０、７０ａ、７０ｂにおいて別の手法を用いる構成である。
【０１８１】
図１７は、第４実施形態による電子透かし埋め込み装置１ｃの概略ブロック図である。電子透かし埋め込み装置１ｃは、埋め込み系列生成部１０ｃと、埋め込みパターン生成部２０ｃと、時間変調部３０ｃと、パターン重畳部４０ｃを備えている。
【０１８２】
図１８は、埋め込み系列生成部１０ｃの内部ブロックを示した図である。埋め込み系列生成部１０ｃは、複素情報系列生成部１１ｃと、埋め込み系列計算部１２ｃを備えている。
次に、第４実施形態における電子透かしの埋め込みの処理について説明する。第４実施形態の電子透かし埋め込み装置１ｃの構成は、以下の点で第１実施形態と異なり、それ以外の符号から「ｃ」を除いた符号が第１実施形態の機能部の符号と同じ機能部の構成については第１実施形態と同様である。以下、異なる構成について説明する。
【０１８３】
電子透かし埋め込み装置１ｃの埋め込み系列計算部１２ｃは、以下の手順により埋め込み系列ｗ３の算出を行う。
複素情報系列生成部１１ｃが生成した複素情報系列ｓ３を複素数の擬似乱数列ＰＮ＝｛ＰＮ１，ＰＮ２，…，ＰＮＬｋ｝で拡散し、長さＬの埋め込み系列ｗ３を算出する。ここで、ｋは、複素情報系列ｓ３ｉの個数（ｉ＝１〜ｋ）を示す。具体的には、複素情報系列ｓ３ｉに対して、式（３５）で示される埋め込み系列ｗ３（ｉ）を式（３６）のように算出する。なお、式（３５）において、ｗが埋め込み系列ｗ３に対応し、式（３６）のｓｉは、複素情報系列ｓ３ｉに対応する。
【０１８４】
【数３５】
2014022972.tif 000037
【０１８５】
【数３６】
2014022972.tif 000038
【０１８６】
さらに、埋め込み系列ｗ３、すなわちｗ＝｛ｗ１，ｗ２，…，ｗＬ｝とするとき、それぞれのｗをｗ（ｉ）の各ｉについての総和として式（３７）のように求める。
【０１８７】
【数３７】
2014022972.tif 000039
【０１８８】
なお、複素数の擬似乱数列ＰＮは、Ｍ系列やＧＯＬＤ系列から得られる値を用いて生成することができる。また、擬似乱数列ＰＮは、処理の際に演算によって生成する構成の他、予め演算によって生成しておき、内部のメモリ等の記憶手段に予め記憶させておいたものを読み出すことによって利用する構成であってもよい。
【０１８９】
図１９は、第４実施形態における電子透かし検出装置２ｃの概略ブロック図である。
電子透かし検出装置２ｃは、時間復調部５０ｃと、検出系列抽出部６０ｃと、相関値計算部７０ｃと、透かし判定部８０ｃとを備え、透かし埋め込み済み画像Ｈ３を入力し、透かし検出結果Ｒ３を出力する。なお、透かし埋め込み済み画像Ｈ３は、電子透かし埋め込み装置１ｃが生成した透かし埋め込み済み画像Ｉ３’と同じ映像である場合の他、電子透かしが埋め込まれていない画像の場合も含む。
【０１９０】
次に、第４実施形態における電子透かしを検出する処理について説明する。第３実施形態の電子透かし検出装置２ｃの構成は、以下の点で第１実施形態と異なり、それ以外の符号から「ｃ」を除いた符号が第１実施形態の機能部の符号と同じ機能部の構成については第１実施形態と同様である。以下、異なる構成について説明する。
電子透かし検出装置２ｃの相関値計算部７０ｃにおける相関値ρ３の算出の手順を以下に説明する。
【０１９１】
（１）電子透かし埋め込み装置１ｃにおいて埋め込み系列の生成の際に用いたものと同じ複素数の擬似乱数列ＰＮ＝｛ＰＮ１，ＰＮ２，…，ＰＮＬｋ｝を用いて処理を行う。当該擬似乱数系列ＰＮの長さは、Ｌｋであり、電子透かし埋め込み装置１ｃの場合と同じく、処理の際に演算によって生成する構成の他、予め演算によって生成しておき、内部のメモリ等の記憶手段に予め記憶させておいたものを読み出すことによって利用する構成などを適用することができる。
（２）擬似乱数列をＬ個ずつに分割し、ｊ番目の部分系列をＰＮ（ｊ）とする。
（３）検出系列ｄ３をｄ＝ｄ｛ｄ１，ｄ２，…，ｄＬ｝とし、擬似乱数列の部分系列ＰＮ（ｊ）との複素相関を計算する。すなわち、ｊ番目の部分に対する相関値をρ（ｊ）とすると式（３８）のように相関値を求めることができる。
【０１９２】
【数３８】
2014022972.tif 000040
【０１９３】
ただし、「・」は、数列をベクトルとみたときの内積演算を表し、Ｃ＊は複素数Ｃの共役複素数を表す。なお、検出信頼性の評価基準を一致させるため、上述した正規化の手法により正規化した上で演算を行うようにしてもよい。
なお、上記の第４実施形態については、第１実施形態との構成の相違の点について説明したが、第４実施形態の構成は、第２、３実施形態にも同様に適用することができる。
【０１９４】
第４実施形態の埋め込み系列の計算方法では、１シンボルが長さＬの埋め込み系列全体に拡散することになり、透かし埋め込み対象映像Ｉ３の埋め込み系列の一部に対応する成分が、バースト的に大きな値となって電子透かしに対するノイズとなる状況においても、検出の信頼性を高めることができる。
一方、第１実施形態の埋め込み系列の計算方法では、第４実施形態の方法と異なり、１シンボルは埋め込み系列の一部分に対してのみ拡散されることとなるが、埋め込み系列の１つ値は１つのシンボルの値を元に生成され、複数のシンボルの値が混在して表現されることがない。これにより、第４実施形態の方法と比較して、透かし埋め込み済み信号が量子化されるなど、微小な透かし成分が除去されるノイズが加えられた場合に対する耐性が高まるという特性がある。
【０１９５】
上記の第１から第４の実施形態の構成により、電子透かし埋め込み装置において、複素数のように位相もしくは方向を持つ値からなる情報系列を用いて、擬似乱数に基づく拡散系列の要素値の位相もしくは方向を、情報系列の値の位相もしくは方向の分だけ変位させ、埋め込み系列を求めてこれを電子透かしとして埋め込むことが可能となる。そして、電子透かし検出装置において、検出系列と拡散系列の複素相関値の位相もしくは方向を、所定のシンボル値境界と比較し、埋め込まれていた情報系列の値を判定することで電子透かしの検出を行うことにより、検出の信頼性を保ったまま、埋め込みの情報量を増加させることが可能となる。
【０１９６】
なお、本発明は、上記の実施形態の構成に限られず、映像、画像以外の信号を対象に同様の電子透かし埋め込みおよび電子透かし検出を行うようにしてもよい。例えば、多チャネルを持つ音声信号を対象に同様の電子透かし埋め込みおよび電子透かし検出を行うこともできる。
【０１９７】
また、上記の実施形態において複素数を用いている構成は、複素数に限られるものではなく、４元数、８元数、ｋ次元（ｋは２以上）のベクトルを用いるようにしてもよい。また、位相や偏角を用いている箇所はベクトルの方向や所定の軸（例えば、Ｘ、Ｙ、Ｚ軸など）からの角度を用いるようにしてもよい。また、相関値計算においては、４元数や８元数の場合、複素数と同様に共役な数との積和をとってもよい。任意のｋ次元ベクトルの場合は、複素数、４元数、８元数と同様の積演算を行ってもよいし、ベクトル間の内積を行うようにしてもよい。
また、上記の実施形態において、埋め込み系列生成部での埋め込み情報の処理に先立って誤り訂正符号を用いて埋め込み情報を符号化してもよく、逆に、検出結果の出力に先立って誤り訂正符号を復号するようにしても良い。
【０１９８】
また、上記の実施形態において、電子透かし埋め込み装置が出力した電子透かし埋め込み済信号を直接電子透かし検出装置が入力する場合についての記載も行っているが、本発明の構成は当該構成に限られるものではなく、埋め込み済信号を圧縮、符号化、配信、編集、改変などした上で電子透かし検出装置が入力するようにしてもよい。また、埋め込み済信号を磁気媒体、例えば、ＶＴＲ(video tape recorder)、ＤＶＤ(Digital Versatile disc)、フレキシブルディスク、ＣＤ（compact disc）、ＨＤＤ（Hard Disc Drive）等やその他の媒体、例えば、フィルム等に一旦記録し、それを読み出すこと、ネットワークを通じて伝送することも装置間の入出力に含まれるものとする。また、光学的なデバイスを用いて再生（例えば、映画としてスクリーンに映写する、ＣＲＴ（Cathode ray tube）や液晶、プラズマなどのディスプイで表示する等）したものをビデオカメラ、携帯電話のカメラ、フィルムを用いたカメラ等の撮影手段を用いて再撮影することなども装置間の入出力に含まれるものとする。
【０１９９】
なお、本発明における各処理部の機能を実現するためのプログラムをコンピュータ読み取り可能な記録媒体に記録して、この記録媒体に記録されたプログラムをコンピュータシステムに読み込ませ、実行することによりアセンブリプログラムの変換を行ってもよい。なお、ここでいう「コンピュータシステム」とは、ＯＳや周辺機器等のハードウェアを含むものとする。また、「コンピュータシステム」は、ホームページ提供環境（あるいは表示環境）を備えたＷＷＷシステムも含むものとする。また、「コンピュータ読み取り可能な記録媒体」とは、フレキシブルディスク、光磁気ディスク、ＲＯＭ、ＣＤ−ＲＯＭ等の可搬媒体、コンピュータシステムに内蔵されるハードディスク等の記憶装置のことをいう。さらに「コンピュータ読み取り可能な記録媒体」とは、インターネット等のネットワークや電話回線等の通信回線を介してプログラムが送信された場合のサーバやクライアントとなるコンピュータシステム内部の揮発性メモリ（ＲＡＭ）のように、一定時間プログラムを保持しているものも含むものとする。
【０２００】
また、上記プログラムは、このプログラムを記憶装置等に格納したコンピュータシステムから、伝送媒体を介して、あるいは、伝送媒体中の伝送波により他のコンピュータシステムに伝送されてもよい。ここで、プログラムを伝送する「伝送媒体」は、インターネット等のネットワーク（通信網）や電話回線等の通信回線（通信線）のように情報を伝送する機能を有する媒体のことをいう。また、上記プログラムは、前述した機能の一部を実現するためのものであっても良い。さらに、前述した機能をコンピュータシステムにすでに記録されているプログラムとの組み合わせで実現できるもの、いわゆる差分ファイル（差分プログラム）であっても良い。
【符号の説明】
【０２０１】
１  電子透かし埋め込み装置
１０  埋め込み系列生成部
１１  複素情報系列生成部
１２  埋め込み系列計算部
２０  埋め込みパターン生成部
３０  時間変調部
４０  パターン重畳部
２  電子透かし検出装置
５０  時間復調部
６０  検出系列抽出部
７０  相関値計算部
８０  透かし判定部
【図１】
2014022972.tif 000041
【図２】
2014022972.tif 000042
【図３】
2014022972.tif 000043
【図４】
2014022972.tif 000044
【図５】
2014022972.tif 000045
【図６】
2014022972.tif 000046
【図７】
2014022972.tif 000047
【図８】
2014022972.tif 000048
【図９】
2014022972.tif 000049
【図１０】
2014022972.tif 000050
【図１１】
2014022972.tif 000051
【図１２】
2014022972.tif 000052
【図１３】
2014022972.tif 000053
【図１４】
2014022972.tif 000054
【図１５】
2014022972.tif 000055
【図１６】
2014022972.tif 000056
【図１７】
2014022972.tif 000057
【図１８】
2014022972.tif 000058
【図１９】
2014022972.tif 000059