(57)【要約】【課題】高い画像品質を保持しかつ予期される最悪のノイズ・レベルに適応できるウォーターマーク埋め込み手段を提供する。【解決手段】変倍ビン符号化方式と呼ばれる新符号化方式を利用する変換ドメイン・ウォーターマーク埋め込み技法を提供する。該方式に従って、高い画像品質を保持しながら予期されるノイズ・レベルに適応するように変換係数を修正することによって変換係数に埋め込まれるべきウォーターマーク・メッセージが符号化される。ウォーターマークの取り出しは、変換係数の統計値、画像処理操作や悪意ある攻撃によって係数に持ち込まれるノイズの最悪ケース統計モデルに基づいて最尤法を使用する復号法によって実行される。【選択図】図２
2004507177.tif 000002
【特許請求の範囲】
【請求項１】
関連するドメイン表現を有するデジタル画像データへ電子透かしを埋め込む方法であって、
誤り訂正コードを使用して電子透かしデータを符号化して符号化済み電子透かしデータを生成するステップと、
変倍ビン符号化方式を使用して、符号化済み電子透かしデータを画像のドメイン表現のエレメントの第１のサブセットに埋め込み、該ドメイン表現の電子透かし埋め込みサブセットを生成するステップと、
電子透かし埋め込みエレメント・サブセットを画像の該ドメイン表現の電子透かしの埋め込まれてない残りのエレメント部分と組み合わせて、画像の電子透かし埋め込みドメイン表現を生成するステップと、
を含む方法。
【請求項２】
符号化済み電子透かしデータを画像のドメイン表現のエレメントの第１のサブセットに埋め込む前に、デジタル画像データを空間ドメイン表現から変換ドメイン表現へ変換するステップと、
画像の電子透かし埋め込みドメイン表現を逆変換して、デジタル画像データおよび電子透かしデータを含む電子透かし埋め込みデジタル画像の空間ドメイン表現を生成するステップと、
を更に含む請求項１に記載の方法。
【請求項３】
変倍ビン符号化方式が、
画像データの統計値、予期ノイズ・レベルおよび許容歪みモデルを表す変倍パラメータで画像のドメイン表現のサブセットの各エレメントを変倍するステップと、
変倍された各エレメントに埋め込まれるべき対応する符号化済み電子透かしデータ・ビットの論理状態に応じて非対称離散マッピングのペアの１つに変倍された各エレメントをマップするステップと、
変倍された各エレメントとその対応するオリジナル・エレメントとの間の差を取得するステップと、
前記差をその対応するオリジナル・エレメントに加算して、電子透かしの埋め込まれた各エレメントを取得するステップと、
を含む、請求項２に記載の方法。
【請求項４】
前記画像の空間ドメイン表現を埋め込まれノイズのある電子透かしから、該電子透かしを取り出すプロセスを含み、
該プロセスが、
ノイズのある電子透かし埋め込み画像の空間ドメイン表現を変換ドメイン表現に変換するステップと、
ノイズのある電子透かし埋め込み画像の変換ドメイン表現から前記第１のエレメント・サブセットに対応する第２のエレメント・サブセットを選択するステップと、
該第２のエレメント・サブセットの統計パラメータを決定するステップと、
該統計パラメータ、予期ノイズと許容歪みとのモデルおよび前記変倍ビン符号化方式に関連する変倍パラメータを使用する確率型復号方式を用いて該第２のエレメント・サブセットから電子透かしの埋め込まれた符号化データを抽出するステップと、
誤り訂正コードに基づく誤り訂正復号器を使用して抽出された電子透かし埋め込み符号化データを復号して、電子透かしデータを生成するステップと、
を含む、請求項２に記載の方法。
【請求項５】
埋め込まれた電子透かしデータを抽出する前記ステップが、
第２のサブセットのエレメントの各々についてかつ離散マッピングの各々について、前記第１のエレメント・サブセットの範囲内の対応するオリジナル・エレメントを符号化するため所与のマッピング値が使用された第１の確率および前記所与のマッピング値が使用された前記第１の確率を所与として前記第２のエレメント・サブセットの各エレメントを観察する第２の確率を決定するステップと、
前記第１および第２の確率から決定される可能な電子透かしデータ符号ワードの各々についてスコア値を決定するステップと、
スコア値を最大にする電子透かしデータ符号ワードを選択して、選択した電子透かしデータ符号ワードを推定電子透かしデータとするステップと、
を含む、請求項４に記載の方法。
【請求項６】
変換ドメイン表現が関連係数を有するＤＣＴ変換ドメインに対応する、請求項５に記載の方法。
【請求項７】
対角線平行斜線の順序に従ってＤＣＴドメイン表現の第１の変換係数サブセットを配列して、対角線平行斜線に対応する第１の変換係数サブセットを選択するステップを更に含み、
該第１の変換係数サブセットの係数グループが共通の歪みおよび分散特性を有する、
請求項６に記載の方法。
【請求項８】
第１の変換係数サブセットを選択する前記ステップの前に、変換ドメイン表現係数を配列するステップを更に含み、
第２のサブセットの係数が第１のサブセットの係数と同じ変換ドメイン位置を有する、請求項７に記載の方法。
【請求項９】
電子透かし符号化データを抽出するステップがアダマール復号方式を使用する、請求項８に記載の方法。
【請求項１０】
デジタル画像データへ電子透かしを埋め込むシステムであって、
デジタル画像データを変換ドメイン表現に変換し、第１の変換係数を生成し、
誤り訂正コードを使用して電子透かしデータを符合化して符号化済み電子透かしデータを生成し、
変倍ビン符合化方式を使用することによって符号化済み電子透かしデータを第１の変換係数の第１のサブセットに埋め込み、電子透かし埋め込み変換係数サブセットを生成し、
電子透かしが埋め込まれてない残りの第１の変換係数と共に電子透かしの埋め込まれた第１の係数サブセットを逆変換して、デジタル画像データおよび電子透かしデータを含む電子透かし埋め込みデジタル画像データを生成する、
第１のデジタル画像プロセッサと、
電子透かし埋め込み画像データを第２の変換係数の変換ドメイン表現に変換し、
第１の変換係数サブセットに対応する第２の変換係数の第２のサブセットを選択し、
第２の係数サブセットの統計を使用して第１の変換係数サブセットの統計パラメータを推定し、
あらかじめ定められた変倍パラメータ、変倍およびマッピング・ステップの既知の諸側面、予期ノイズと許容歪みのモデルおよび推定された統計パラメータを使用する確率型復号方式を用いて第２の係数サブセットから埋め込まれた電子透かしデータを抽出し、
前記誤り訂正コードに基づく誤り訂正復号器を使用して抽出された電子透かしデータを復号して電子透かしデータを生成する、
第２のデジタル画像プロセッサと、
を備えるシステム。
【発明の詳細な説明】
【０００１】
【発明の属する技術分野】
本発明は、デジタル・データの電子透かし（ｗａｔｅｒｍａｒｋ、ウォーターマーク）埋め込みに関し、特に、デジタル画像データに電子透かしデータを埋め込み、オリジナルのデジタル画像データを使用せずに電子透かし埋め込みデジタル画像のノイズ・バージョンから埋め込まれた電子透かしデータを取り出す技法に関する。
【０００２】
【従来の技術】
デジタル情報技術の急速な開発が、（例えば、オープンにアクセス可能なウェブ・サイトを使用する）インターネットのような通信ネットワークを介するデジタル・データの複製および流通の容易化とともに、コンテンツ・プロバイダをしてその（画像、映像、音声、テキストなどの）電子データを保護する課題に直面させている。このような技術の急速な開発は、著作権保護のための高度な電子透かし埋め込みおよび情報秘匿手法の設計および研究に向けた多くの開発努力を刺激してきた。
【０００３】
例えば画像に秘密情報を埋め込み、この情報が秘密キーにアクセスせずに除去または判読されることができないようにするため、電子透かし埋め込み手法が使用される。一方、電子透かしがオリジナル画像に埋め込まれたとしても、品質を維持するためには、電子透かし埋め込み画像はオリジナル画像と知覚上同一でなければならない。電子透かし埋め込みの主なアプリケーションは、データの所有権の提供および偽造者に対する防護にある。従って、画像が役に立たなくなるほどに画像に損傷を与えることなく悪意を持つ攻撃者が電子透かしを除去することができないことを保証することが極めて重要である。更に、電子透かし埋め込み技法は、プリント、スキャン、有損失圧縮（例えばＪＰＥＧ）、フィルタリング等々のように画像に対して実行される標準的操作に対して強力であることが必要とされる。
【０００４】
従って、電子透かし埋め込み手法の設計において、攻撃と画像処理操作に対する堅固性、高品質および埋め込まれるべき秘密情報の相対的量（すなわち符号化率）という３つの主要なパラメータの観点からトレードオフが考慮される。
【０００５】
従来の画像電子透かし埋め込み手法は、電子透かしが埋め込まれる機能セットのタイプに従って２つのグループに分けられる。すなわち、電子透かしは、画像の空間ドメイン表現における輝度の強度値または画像の変換ドメイン表現（例えばＤＣＴ、ＤＷＴ）における変換係数に埋め込まれる。電子透かし信号を検出するために使用されるアルゴリズムもまた、電子透かし検出プロセスの際にオリジナルの画像を使用するか否かに基づいて２つのカテゴリに分類される。一般的には、電子透かしを変換ドメインの係数に埋め込みオリジナルの画像を必要とせずに電子透かしを検出するという電子透かし埋め込み手法は、標準画像処理操作および悪意のある攻撃に対する堅固性において多くの利点を有する。
【０００６】
変換ドメインの範囲内において電子透かしデータを符合化することによって画像Ｉに電子透かしメッセージまたはデータを埋め込む基本手法は、画像Ｉを初期的に変換して画像の変換ドメイン（例えばＤＣＴ、ＤＦＴその他）表現を取得することによって実行される。次に、変換係数のサブセットが選択され、変換係数値のサブセットを若干変更することによって電子透かしメッセージが符合化され、その結果、電子透かしを埋め込まれた変換係数が生成される。次に、電子透かし埋め込み係数が、電子透かしを埋め込まれてない係数と組み合わせられ、画像の電子透かし埋め込み変換表現が生成される。次に、電子透かし埋め込み変換表現は、逆変換され、電子透かし埋め込み画像
2004507177.tif 000003
が生成される。
【０００７】
2004507177.tif 000004
がオリジナルの画像Ｉと視覚的に見分けがつかないようにするため、巧妙な係数修正が行われる。
【０００８】
電子透かし埋め込み画像
2004507177.tif 000005
から電子透かしを取り出すため、復号器は、
2004507177.tif 000006
の歪められたまたはノイズのあるバージョンを受け取る。すなわち、この画像バージョンは、一定の画像処理操作または悪意のある攻撃を受けた可能性を持つ。復号器は、受け取ったバージョンを適切な変換ドメイン表現に変換し、電子透かし信号が符合化されたものと同じ係数サブセットを選択する。復号器は、次に、復号手順を使用して係数サブセットから電子透かしメッセージを抽出する。
【０００９】
オリジナルの画像に頼ることなく復号するため、変換ドメインにおける電子透かし埋め込みを行う２つの方法、すなわち、係数摂動およびディザリング変調が最近提案されている。
【００１０】
係数摂動を使用する１つの従来技術方法においては、有意なＤＣＴ（すなわち変換）係数の値を摂動することによって電子透かしデータがオリジナルの画像に付加される。例えば、ｓがＤＣＴ係数であるとすれば、係数をｓ＋εに変更することによって電子透かしデータの”ゼロ”ビットが符号化され、係数をｓ−εに変更することによって電子透かしデータの”１”ビットが符号化される。ここでεは小さい定数である。電子透かし埋め込み画像
2004507177.tif 000007
から電子透かしメッセージの取り出しは、適切なＤＣＴ係数を電子透かしメッセージに相関させることによって実行される。このような電子透かし埋め込み手法のバリエーションとして、例えば、ＤＦＴまたはＤＷＴのような異なる変換ドメインの使用や異なる摂動方式（例えば、ｓ（１±ε）またはｓ±ε｜ｓ｜）などが提案されている。
【００１１】
ディザリング変調は、変換ドメイン係数を量子化することに基づく方式である。この場合、選択された係数サブセットへの電子透かしメッセージの埋め込みは、電子透かしメッセージに依存する形態でこれら係数を量子化された値に置き換えることに基づく。例えば、”ゼロ”ビットは、量子化因子ｑ０で係数を量子化することによって符合化され、”１”ビットは、量子化因子ｑ１で係数を量子化することによって符合化される。
【００１２】
ディザー変調電子透かし埋め込み方式の１つのバリエーションにおいては、”ゼロ”ビットおよび”１”ビットの値は、２つの”自己ノイズ抑止”関数ｆ０およびｆ１を使用して符合化される。これら関数は、ディザー変調方式において使用される量子化因子ｑ０およびｑ１を若干修正したものに基づく。
【００１３】
【発明が解決しようとする課題】
現在まで、係数摂動およびディザー変調を含む従来技術電子透かし埋め込み方式は、定義済み歪みレベルまたはノイズ基準に適応するように明示的に設計されてない。その結果、これらの方式は、比較的低率の情報埋め込みを行うか、さもなければ、可変的ノイズ条件の下で高いエラー率に苦しむ。更に重要な点であるが、電子透かし埋め込み画像データに対する悪意のある攻撃者による攻撃の試みが未知のノイズ条件という形式を持ち込むので、これらの従来技術方式においては、電子透かし埋め込みの主な目的の１つであるデータのセキュリティに対する攻撃を防護することが容易でない。
【００１４】
【課題を解決するための手段】
本発明は、変倍ビン符号化と呼ばれる新しい符合化方式に基づく電子透かし埋め込み手法を提示する。本発明に従ったこの新符号化方式は、高い画像品質（すなわち低い歪みレベル）を保持しかつ予期される（すなわち一層悪いケース）のノイズ・レベルに適応するように画像値を修正することによって電子透かしデータを画像データに符号化する。
【００１５】
電子透かしが意図的または非意図的ノイズにさらされた後電子透かし埋め込み画像から電子透かしデータを取り出す処理は、確率型プロシージャを使用する復号方法（例えば最尤法復号法）を介して実行される。この方法は、オリジナル画像値の推定された統計および画像処理操作または攻撃ノイズによって画像に導入されるノイズの予期統計モデルに基づくものであり、それによって、強固で高品質の電子透かし埋め込みシステムおよび方法が提供される。
【００１６】
本発明は、デジタル画像データへデジタル・電子透かしデータを埋め込み、オリジナルのデジタル画像データを使用することなくかつ電子透かし埋め込みデジタル画像データに導入されたノイズに影響されることなく、電子透かしの埋め込まれたデジタル・データから埋め込まれた電子透かしを取り出す方法およびシステムを提供する。本発明のシステムおよび方法は、オリジナルのデジタル・データに電子透かしを埋め込み、電子透かし埋め込みデジタル・データを取得するプロセスおよび電子透かし埋め込みデジタル・データを復号してオリジナルの電子透かし埋め込みデータを取得するプロセスという２つのプロセスによって実行される。電子透かしデータは、変倍ビン符号化と呼ばれる符合化方法を使用してデジタル画像データに埋め込まれる。埋め込まれた電子透かしデータは、確率型復号方式を使用して電子透かし埋め込みデータから取り出される。１つの側面において、確率型復号方式は最尤法復号方式である。
【００１７】
本発明の１つの側面において、電子透かしデータの埋め込みは、オリジナルのデジタル画像データを第１の変換係数に変換するステップ、誤り訂正コードを使用して電子透かしデータを符合化するステップ、変倍ビン符号化方式を使用して、符号化済み電子透かしデータを前記第１の変換係数の第１のサブセットに埋め込み、埋め込まれた変換係数のサブセットを生成するステップ、および、電子透かしが埋め込まれてない残りの第１の変換係数と共に電子透かしの埋め込まれた第１の係数サブセットを逆変換して、オリジナルのデジタル画像データおよび電子透かしデータを含む電子透かし埋め込みデジタル画像データを生成するステップによって達成される。
【００１８】
上記側面によれば、変倍ビン符号化方式は、予期ノイズ・レベルおよび許容歪みモデルを表すあらかじめ定められた変倍パラメータで変換係数サブセットの各係数を変倍するステップ、変倍された各係数に埋め込まれるべき対応する符号化済み電子透かしデータ・ビットの論理状態に応じて非対称離散マッピングのペアの１つに変倍された各係数をマップするステップ、変倍された各係数とその対応するマップされ変倍された係数との間の差を取得するステップ、および、前記差をその対応するオリジナルの（すなわち変倍されてなくマップされてない）係数に加算して、電子透かしの埋め込まれた各変換係数を取得するステップによって実行される。
【００１９】
本発明のシステムおよび方法の別の１つの側面によれば、電子透かしの符号化の間に使用されたものと同じ変換を使用して電子透かしの埋め込まれた画像を第２の変換係数に変換するステップ、変換係数の第１のサブセットに対応する電子透かし埋め込み変換係数の第２のサブセットを選択するステップ、第２の電子透かし埋め込み係数サブセットの画像統計を使用して第１の変換係数サブセットの統計パラメータを推定するステップ、あらかじめ定められた変倍パラメータ、変倍およびマッピング・ステップの既知の諸側面、予期ノイズと許容歪みのモデルおよび推定された統計パラメータを使用する確率型復号方式を介して第２の係数サブセットから埋め込まれた電子透かしデータを抽出するステップ、および、電子透かし符号化の際に使用された誤り訂正コードに基づく誤り訂正復号器を使用して抽出された電子透かしデータを復号して推定オリジナル・電子透かしを取得するステップによって、電子透かしデータが電子透かしの埋め込まれたデジタル画像データから取り出される。１つの側面において、最尤法復号方式が確率型復号方式として使用される。
【００２０】
電子透かし埋め込みシステムの１つの側面は、デジタル画像データへ電子透かしを埋め込むため、デジタル画像データを変換ドメイン表現に変換する画像変換器、誤り訂正コードを使用して電子透かしデータを符合化して電子透かしデータを符合化する符号化器、変倍ビン符合化方式を使用することによって符号化済み電子透かしデータを第１の変換係数の第１のサブセットに埋め込み、電子透かし埋め込み変換係数サブセットを生成する電子透かし符号化器、および、電子透かしが埋め込まれてない残りの第１の変換係数と共に電子透かしの埋め込まれた第１の係数サブセットを逆変換して、デジタル画像データおよび電子透かしデータを含む電子透かし埋め込みデジタル画像データを生成する逆画像変換器を備える。
【００２１】
電子透かし埋め込みシステムの１つの側面は、電子透かし埋め込みデジタル画像から電子透かしデータを取り出すため、電子透かし埋め込み画像データを変換ドメイン表現に変換する画像変換器、係数の第２のサブセットの統計パラメータを決定する手段、あらかじめ定められた変倍パラメータ、変倍およびマッピング・ステップの既知の諸側面、予期ノイズと許容歪みのモデルおよび決定された統計パラメータを使用する確率型復号方式を介して第２の係数サブセットから埋め込まれた電子透かしデータを抽出する確率型電子透かし復号器、および、 オリジナルの誤り訂正コードに基づく誤り訂正復号器を使用して抽出された電子透かしデータを復号して電子透かしデータを取得する復号器を備える。
【００２２】
【発明の実施の形態】
電子透かし埋め込みプロセスは、電子透かし符号化段階１０、攻撃チャネル段階１１および電子透かし復号段階１２という３つの段階を有するものとして表されることができる。
【００２３】
電子透かし符号化段階１０は、埋め込まれるべき電子透かしデータまたはメッセージｍおよび電子透かしメッセージを埋め込むための埋め込み対象画像ｉという２つの入力を受け取る。埋め込み対象画像ｉは、テキスト・データ、ビデオ・データ、オーディオ・データまたはこれらタイプのデータの組合せからなるデジタル・データを含むことができる。符号化器の出力は、電子透かしを埋め込まれた画像ｉ’であり、これは視覚上オリジナルの画像ｉに類似している。攻撃チャネル１１は、予期される／悪意のない信号処理操作（例えば圧縮、フィルタリング、その他）および／あるいは電子透かし埋め込み画像に対する意図的攻撃の複合を表す。チャネルの電子透かし埋め込み画像出力ｉ＊は、入力された電子透かし埋め込み画像ｉ’の歪められたバージョンであり、ノイズ信号ｒから生じるものと仮定される。最後に、復号器１２は入力として画像ｉ＊を受け取り、オリジナルの埋め込み対象画像ｉに頼ることなく埋め込まれたメッセージｍを推定する。
【００２４】
変倍ビン符号化方式を使用する実施形態
図２の（Ａ）は、埋め込み対象画像に電子透かしメッセージを埋め込む電子透かし埋め込みシステムの１つの実施形態を示す。画像データｉへの電子透かしメッセージｍの埋め込みは変換ドメインにおいて実行される。画像ｉは、画像変換ブロック２０によって変換係数セットｔ＝Ｆ（ｉ）に初期的に変換される。画像変換ブロック２０は、デジタル信号処理の分野において既知の多数の方法で実施されることができる。例えば、画像変換ブロックは、離散コサイン変換（ＤＣＴ）表現または離散フーリェ変換（ＤＦＴ）表現に画像を変換するように実施されることができる。また、埋め込みの前に、電子透かしメッセージｍは誤り訂正コード・ブロック２１によって誤り訂正符号形式ｂに符合化される。電子透かしデータを埋め込むため係数サブセットｓが画像係数ｔから選択される。変倍ビン符号化ブロック２２は、特定の符号化方式を使用して、選択された画像データ変換係数ｓにメッセージｂを埋め込み、電子透かしの埋め込まれた係数ｘを生成する。この場合特定符号化方式においては、予期ノイズ・モデルと許容歪みモデルとの関数であるあらかじめ定められた変倍パラメータで変換係数サブセットの各係数が変倍され、変倍された各係数に埋め込まれるべき対応する符号化済み電子透かしデータ・ビットの論理状態に応じて非対称離散マッピングのペアの１つに変倍された各係数がマップされ、変倍された各係数とその対応するマップされた変倍係数との間の差が取得され、前記差がその対応するオリジナルの（すなわち変倍されてなくマップされてない）係数に加算されて電子透かしの埋め込まれた各変換係数が取得される。
【００２５】
電子透かし埋め込み係数ｘは対応する非電子透かし埋め込み係数ｓを置き換え、残りの非電子透かし埋め込み係数と組み合わされる。組み合わされた係数が画像逆変換ブロック２３によって逆変換され、電子透かし埋め込み画像ｉ’が生成される。
【００２６】
ここで理解されるべき点であるが、ビン符号化は、画像それ自体と符号化済み画像とに応じて所与のメッセージを符号化するため符号ワードが選択される手法であると一般に定義される。本発明のシステムおよび方法に従った変倍ビン符号化は、変換係数を崩壊させると予想されるほどのノイズ・レベルおよび許容歪みに符号化率（すなわちビンの数）を適応させるように適切な変倍パラメータをビン符号化方式に導入する技法として一般的に定義される。従って、本発明に従った変倍ビン符号化方式に基づく電子透かし埋め込みシステムおよび方法は、画像処理操作および悪意ある攻撃によって導入されるノイズに対して強固である。
【００２７】
画像変換器２０、符号化器２１、変倍ビン符号化器２２および画像逆変換器２３を含む（図２（Ａ）に示されている）システムは、コンピュータまたは処理システムによって実行されるソフトウェア実施形態で実現されることができるし、あるいは、ＡＳＩＣのような専用ハードウェア実施形態またはソフトウェアとハードウェアの任意の組み合わせでも実現することができる。
【００２８】
図２（Ｂ）は、図２（Ａ）に示されるシステムに対応する変倍ビン符合化方式を使用する電子透かし埋め込み方法の１つの実施形態を示す。図２（Ｂ）に示される方法は、画像データｉを変換係数に変換するステップ（ステップ２４）、電子透かしデータｍを誤り訂正コード形式に符合化して符号化済み電子透かしデータｂを生成するステップ（ステップ２５）、変倍ビン符号化方式を使用して符合化された電子透かしデータｂを変換係数サブセットｓに埋め込むステップ（ステップ２６）、および、電子透かし埋め込み係数ｘおよび残りの非電子透かし埋め込み係数を逆変換して電子透かし埋め込み画像ｉ’を生成するステップ（ステップ２７）を含む。
【００２９】
注意されるべき点であるが、画像への電子透かしデータの変倍ビン符号化は、画像の所与の変換ドメインの代わりに、画像の空間（すなわちピクセル）ドメインの範囲内で実行されることもできる。従って、電子透かし埋め込みシステムの１つの代替実施形態においては、画像変換ブロック２０および画像逆変換ブロック２３（図２（Ａ））は必要でなく、電子透かし埋め込み方法の１つの代替実施形態においては、変換および逆変換のステップ２４、２７は必要とされない。
【００３０】
１つの実施形態において、変倍ビン符号化方式を使用する埋め込みは、以下の許容歪みモデルおよび予期ノイズ・モデル基準を厳守することによって電子透かし埋め込み画像がオリジナル画像に視覚的に類似していることを保証するように実行される。
【００３１】
１．歪み制約ｄ（ｘ，ｓ）＜Ｄと整合するスカラー歪みレベル
2004507177.tif 000008
を両者が有するスカラー離散マッピング・ペアｑ０およびｑ１。
２． その値が攻撃チャネルにおける予期ノイズの相対的パワーに依存するスカラー・パラメータ０＜α＜１。
【００３２】
このように、Ｎ個ビットのシーケンスｂ＝（ｂ１，．．．，ｂＮ）がＮ個変換係数のシーケンスｓ＝（ｓ１，．．．，ｓＮ）に埋め込まれる電子透かし埋め込みシステムおよび方法の別の１つの実施形態において、変倍ビン符号化方式を使用する埋め込みは、シーケンスｓをαによって変倍し、次に、ｉ番目のコンポーネントがｂｉの信号に従ってｑ０またはｑ１によってマップされるようにシーケンスｂによって離散マッピングｑによって定義されたコンポーネント方式を使用してシーケンスｓを埋め込むことによって実行される。これを公式で記述すれば、
【数１】
2004507177.tif 000009
と表される。
【００３３】
次に、αｓとそのマップされたバージョンとの間の差がオリジナルのシーケンスｓに加算され、電子透かしの埋め込まれた信号
ｘ＝ｓ＋（ｑ（αｓ）−αｓ）
が取得される。
【００３４】
図３は、予期ノイズ・レベルの関数である変倍パラメータを使用し、許容歪みモデルの関数である離散的マッピングを使用する変倍ビン符号化方式を実行するプロセスを示す。この方法によれば、サブセットｓの各係数ｓｉがあらかじめ定められた変倍パラメータαで変倍され（ステップ３０）、次に、変倍された係数ｓｉに埋め込まれるべき対応する符号化された電子透かしデータ・ビットｂｉの論理状態に応じて非対称離散マッピング（ｑ０、ｑ１）の１つに各係数がマップされる（ステップ３１）。変倍された各係数とその対応するマップされ変倍された変倍係数との間の差が決定され（ステップ３２）、その対応するオリジナルの（すなわち変倍されてなくマップされてない）係数に上記の差が加算され、電子透かしの埋め込まれた各変換係数ｘｉが取得される（ステップ３３）。
【００３５】
確率型復号方式を使用する復号
電子透かし埋め込みシステムおよび方法の別の１つの実施形態に従って、電子透かし復号段階１３（図１）は、オリジナル画像の変換係数に関する統計情報を推定するパラメトリック統計的推定プロシージャおよび攻撃チャネル・ノイズに関する予期される（すなわち一層悪いケースの）モデルを含む確率型復号方式を使用して実施される。確率型復号方式は、オリジナル画像と電子透かし埋め込み画像とは統計的類似性を持ち、従って、電子透かし埋め込み画像を使用してオリジナル画像の推定統計データを取得することが可能とされるという事実を利用する。
【００３６】
確率型復号方式を使用する電子透かし埋め込みシステムの１つの実施形態（図４（Ａ））において、画像変換ブロック４０がノイズのある画像ｉ＊を受け取ってそれを変換表現に変換する。変換ブロック４０によって実行される変換のタイプは、電子透かし符号化プロセスの間に実行されたものと同じである。画像変換ブロック４０は、電子透かしが埋め込まれかつ変換された係数のサブセットｙを電子透かし復号器４１に提供する。電子透かしが埋め込まれ変換された係数のサブセットｙが、電子透かし符号化プロセスに使用されたオリジナル係数サブセット（すなわちｓ）に対応する。電子透かし復号器４１は、確率型復号方式を使用して、埋め込まれた電子透かしデータｂ＊をサブセットから抽出する。抽出された電子透かしデータｂ＊は、誤り訂正コード（ＥＣＣ）復号器４２に送られ、そこで電子透かしデータｂ＊は復号され、推定電子透かしデータｍ＊として出力される。
【００３７】
図４（Ｂ）は、確率型復号方式を使用する電子透かし埋め込み方法の１つの実施形態を示している。この方法に従えば、歪められたまたはノイズのある画像バージョンｉ＊から符合化された電子透かしメッセージｍが推定される。初期的に、ノイズのある電子透かし埋め込み画像ｉ＊が変換ドメイン表現Ｆ（ｉ＊）に変換され（ステップ４３）、電子透かし符号化プロセスにおいて使用された第１のサブセットｓ（ｓ１，．．．，ｓＮ）と同じ変換ドメイン位置を有する第２の変換係数サブセットが選択される（ステップ４４）。第２の変換係数はｙ＝（ｙ１，．．．，ｙＮ）と標記される。次に、歪められた係数ｙ＝（ｙ１，．．．，ｙＮ）の統計を決定することによってオリジナルの係数サブセットｓ（ｓ１，．．．，ｓＮ）の統計パラメータが推定される（ステップ４５）。変倍パラメータを含む変倍ビン符号化プロセスの既知の諸側面、予期ノイズと許容歪みのモデルおよび推定統計パラメータを利用する確率型復号方式を使用して、符号化された電子透かしデータ・ワードｂ＊がｙから抽出される（ステップ４６）。電子透かしデータｂ＊が決定されると、それを復号して、推定復号メッセージｍ＊が取得される（ステップ４７）。
【００３８】
画像変換器４０、電子透かし復号器４１およびＥＣＣ復号器４２を含むシステム（図４（Ａ））は、コンピュータまたは処理システムによって実行されるソフトウェア実施形態で実現されることができるし、あるいは、ＡＳＩＣのような専用ハードウェア実施形態またはソフトウェアとハードウェアの任意の組み合わせでも実現することができる。
【００３９】
変換ドメインの範囲内の係数の”位置”が変換から生成される係数アレイの範囲内の位置を指す点は理解されるべきことである。更に留意すべき点であるが、位置は、変換から生成される係数のアレイの次元に従って１、２、３、あるいはいかなる多次元の位置でもよい。
【００４０】
確率型復号方式が最尤法（ＭＬ法）である１つの実施形態においては、変倍ビン符号化プロシージャによって使用されたものと同じパラメータ（α、ｑ０およびｑ１）、予期ノイズと許容歪みのモデル、および、オリジナルのサブセットｓ（ｓ１，．．．，ｓＮ）の推定統計パラメータが使用される。
【００４１】
図５は、ＭＬ復号プロシージャの方法の１つの実施形態を示す。変倍ビン符号化パラメータ５１、オリジナルのマッピングの既知の諸側面５２、予期ノイズと許容歪みのモデル５３および推定統計パラメータ５４を所与として、以下の確率が決定される（ステップ５５、５６）。
【００４２】
＊対応するオリジナルの係数を符合化する際にポイントｑ０（ｊ）が使用された確率Ｐ０（ｊ）；
＊対応するオリジナルの係数を符合化する際にポイントｑ０（ｊ）が使用されたと仮定した場合に歪められた係数値ｙｉを観察する条件つき確率Ｐ０（ｙｉ｜ｊ）；
＊対応するオリジナルの係数を符合化する際にポイントｑ１（ｊ）が使用された確率Ｐ１（ｊ）；および
＊対応するオリジナルの係数を符合化する際にポイントｑ１（ｊ）が使用されたと仮定した場合に歪められた係数値ｙｉを観察する条件つき確率Ｐ１（ｙｉ｜ｊ）。
【００４３】
（ただし、ｊは有限インターバル｛−Ｊ，Ｊ｝にわたる離散的マッピング・ポイントの整数インデックスである）。
【００４４】
上記確率の結果を使用し、可能なコード・ワードｂの各々に関して（ただしｂは２ｋ個の可能な誤り訂正コード・ワードの範囲にある）、次式のスコアが決定される（ステップ５７）。
【００４５】
【数２】
2004507177.tif 000010
【００４６】
次に、スコアが評価され、スコアＳｃ（ｂ）を最大にする符号ワードｂが選択され、これが推定符合化電子透かしデータ・ワードｂ＊に相当する（ステップ５８）。
【００４７】
図６の（Ａ）および（Ｂ）は、変倍ビン符号化方式（図６（Ａ））およびＭＬ復号方式（図６（Ｂ））を使用する電子透かし埋め込み方法の１つの実施形態を示す。この方法によれば、電子透かしデータの符号化はノイズのある画像ｉ＊から電子透かしデータを抽出する時復号を単純化するように実行される。図６の（Ａ）に示されているように、電子透かしが埋め込まれるべき画像ｉは、初期的に、そのＤＣＴ変換表現ｔ＝ＤＣＴ（ｉ）に変換される（ステップ６０）。
【００４８】
次に、変換係数は、ＤＣＴドメインにおいて、図６（Ｃ）に示されているように、左上隅（＃１）から始めて右下隅（＃Ｍ）まで対角線平行斜線の順序に従って順に並べられる（ステップ６１）。ＤＣＴドメインが画像の周波数領域表現に対応している点に留意する必要がある。図６の（Ｃ）を参照すれば、係数は、左上隅１に位置している最低周波数（すなわちＤＣ）から右下隅Ｍの最高周波数へと斜めに並べられる。対角線平行斜線の各々（例えばＬ１、Ｌｎ）は、その領域における係数のグループ（すなわち対角線平行斜線）を含む周波数の１つの領域を表す。
【００４９】
次に、Ｎ個の変換係数サブセットｓが、Ｌ１からＬｈ（ただし１＜Ｌ１＜Ｌｈ＜Ｍ）までインデックスづけされた対角線平行斜線に対応する配列された係数から選択される（ステップ６２）。これらの係数は、電子透かしメッセージまたはデータがそこに符号化されるべき信号ｓ＝（ｓ１，．．．，ｓＮ）として使用される。領域［Ｌｉ，Ｌｈ］にインデックスを持つ対角線平行斜線の各々ごとに（ｉは係数のインデックスである）、この対角線平行斜線に対応する係数について統計パラメータが計算される（ステップ６３）。この場合、統計パラメータは、ｉ番目の係数の対角線平行斜線上にある係数集合の推定分散Ｑである。特定の対角線平行斜線に対応する係数（複数）は同じ推定分散を共有する。
【００５０】
数値集合｛αｉ｝ｋｉ＝１の推定分散の公式は、統計学の分野において周知のように、次の公式で表される。
【００５１】
【数３】
2004507177.tif 000011
【００５２】
対角線平行斜線ごとに、均一な離散マッピング・ペアが計算される（ステップ６４）。離散マッピングは、許容スカラー歪みレベルに依存する。係数ｓｉに関するこの歪みレベルは次式によって定義される。
【００５３】
【数４】
2004507177.tif 000012
ただし、εは視覚歪みレベルを制御する大域パラメータである。この歪み項は、人間視覚システムのモデルによって誘導されるもので、他の歪み関数と置き換えることもできる。留意されるべき点であるが、同じ対角線平行斜線に対応するすべての係数は、共通の歪みレベルＤｉを持ち、従って、同じ離散マッピングを共有しているので、マッピング計算を大幅に縮減させることができる。このようにして、各歪みレベル・パラメータＤｉに関して、均一なスカラー離散マッピング・ペアｑ０とｑ１が次式によって決定されることができる。
【００５４】
【数５】
2004507177.tif 000013
ただし、ｒｏｕｎｄ（）は最隣接整数値を実際の入力へ返す関数である。
【００５５】
マッピングｑｉ０、ｑｉ１は次式によって表される離散集合における値をそれぞれ取る。
【００５６】
【数６】
2004507177.tif 000014
【００５７】
上記離散集合の対応するエレメントはｑｉ０（ｊ）およびｑｉ１（ｊ）（ただしｊは整数インデックス）として標記される。
【００５８】
変換順序付けのタイプに応じて選択される符号を使用して、後段の復号ステップを最適化するため、電子透かしメッセージｍはＮビットのシーケンスｂ＝（ｂ１，．．．，ｂＮ）として符合化される（ステップ６５）。この場合、ＭＬ復号プロセスの間の確率計算を加速するため、誤り訂正コード・ワードとしてアダマール行列行を使用して、メッセージが符合化される。ｋ個のメッセージ・ビットの各ブロックは、次数ｋのアダマール行列の長さｎ＝２ｋの該当する行にマップされる。ここで、ｋは、対角線平行斜線の関数としてまた予期ノイズ・パワーの関数として変わり得るものである。本実施形態において、Ｎ＝ｎＲ（ただしＲは整数）と仮定されている。
【００５９】
電子透かし埋め込み係数ｘを取得するため、変倍ビン符号化方式を使用して、バイナリ・シーケンスｂ＝（ｂ１，．．．，ｂｎ）が係数サブセットｓ＝（ｓ１，．．．，ｓＮ）に埋め込まれる（ステップ６６）。サブセットｓ＝（ｓ１，．．．，ｓＮ）を使い切るまで、ｎ個の係数の各ブロックに関して変倍ビン符号化が実行される。係数は、初期的に、式α＝ε／（ε＋η）によって定義される変倍パラメータを使用して変倍される。ただし、εは前述のように定義され、ηは、予期される（すなわち一層悪いケースの）ノイズ・パワーと埋め込み対象信号パワーとの間の比率である。比率パラメータηは、多くの画像および種々のノイズ・タイプに関してその値を測定することによって実験的に決定される。ノイズが強くなるほどηの値は高い。電子透かし埋め込み方式がどの程度強固であるべきかを決定することによって、ノイズの最大レベルが選択され、対応するηの値が決定される。代替実施形態において、パラメータαはインデックスｉに依存し、従って固定的ではないようにすることもできる点留意されるべきである。次に、（ステップ６４において決定された）ベクトル離散マッピングｑＤｎを使用して、次式によって定義されるように変倍係数がマップされる。
【００６０】
【数７】
2004507177.tif 000015
【００６１】
次に、αｓとそのマップされたバージョンとの間の差がオリジナルのシーケンスｓから減算され、以下に標記されるような電子透かし埋め込み信号ｘが取得される。
【００６２】
ｘ＝ｓ＋（ｑ（αｓ）−αｓ）
ＤＣＴ表現ｔにおける係数シーケンスｓが電子透かし埋め込み係数シーケンスｘと置き換えられ、その結果の電子透かし埋め込み変換表現がｔ’として表される（ステップ６７）。次に、電子透かし埋め込みＤＣＴドメイン表現が逆変換され、電子透かし埋め込み画像ｉ’＝ＤＣＴ−１（ｔ’）が取得される（ステップ６８）。
【００６３】
図６の（Ｂ）を参照すれば、電子透かし画像ｉ’が何らかの形式の意図的または非意図的ノイズにさらされ、歪められたまたはノイズのある画像バージョンｉ＊となった後、復号器は電子透かし画像を受け取る。復号器は、最尤法復号方式を使用して、ノイズのある画像ｉ＊から埋め込まれた電子透かしデータを推定抽出する。初期的には、ノイズのある画像ｉ＊のＤＣＴ変換は、ｔ＊ ＝ ＤＣＴ（ｉ＊）と決定される（ステップ６９）。次に、変倍ビン符号化プロセスの場合と同様に、ノイズのある画像の変換係数ｔ＊が順に並べられる（ステップ７０）。この場合、係数は、左上隅（＃１）から始めて右下隅（＃Ｍ）までの対角線平行斜線に従って配列される。一旦係数が順に並べられると、オリジナルのサブセット（すなわち、ＬｉからＬｈ（ただし１＜Ｌｉ＜Ｌｈ＜Ｍ）までインデックスづけされた対角線平行斜線）に対応するＮ個の変換係数サブセットｙ＝（ｙ１，．．．，ｙＮ）が選択される（ステップ７１）。
【００６４】
本発明の方法に従えば、歪められた画像ｉ＊はオリジナルの画像ｉに視覚的に類似していると仮定され、従って、ノイズ・バージョン係数ｙに関して統計パラメータ（例えば分散推定Ｑｉ＊）を決定することによって、オリジナルの電子透かし埋め込みサブセットの統計パラメータを推定することができる。統計パラメータＱｉ＊は、サブセットｓについて決定した場合と同様に、サブセットｙの範囲内のグループ（すなわち対角線平行斜線）に関して決定される（ステップ７２）。
【００６５】
次に、変倍ビン復号プロセスに関して計算された場合と同じように対角線平行斜線の各々について均一離散マッピングのペアｑ０＊およびｑ１＊が計算される（ステップ７３）。
【００６６】
電子透かし埋め込みサブセットの統計パラメータが決定されると、符合化された電子透かし・デーダｂ＊がＭＬ法を使用して抽出される（ステップ７４）。具体的には、対応するオリジナル係数ｓｉを符合化する際にポイントｑｉ０＊（ｊ）およびｑｉ１＊（ｊ）が使用されたと仮定すれば、観察値ｙｉの各々について、値ｙｉを観察する確率Ｐ０（ｙｉ｜ｊ）およびＰ１（ｙｉ｜ｊ）が推定される。ＤＣＴ係数の連結分布に関する独立ガウス・モデル、許容電子透かし歪みおよび攻撃チャネルにおける予期ノイズを使用して、上記確率は次式によって与えられる。
【００６７】
【数８】
2004507177.tif 000016
【００６８】
対応するオリジナル係数ｓｉを符合化する際に対応するポイントｑｉ０＊（ｊ）およびｑｉ１＊（ｊ）を使用する事前確率Ｐ０（ｊ）およびＰ１（ｊ）は、ｙｉと無関係であり、次式によって推定されることができる。
【００６９】
【数９】
2004507177.tif 000017
【００７０】
各アダマール符号ワードｂ（ｂは２ｋ個の可能な誤り訂正符号ワードにわたる）に関して、次式のスコアが計算される。
【００７１】
【数１０】
2004507177.tif 000018
ただし、次式［数１１］で表される差があらゆるｉおよびｂ∈｛０，１｝について無視できるほど小さいように上式のＪは選択される。
【００７２】
【数１１】
2004507177.tif 000019
【００７３】
留意されるべき点であるが、上記定義のスコア関数は迅速アダマール変換を使用して効率的に計算される。このように、迅速アダマール変換を使用して確率計算を実施することができるように電子透かし符号化の間に対角線平行斜線係数を配列および選択することによって電子透かし復号動作の効率が向上される。
【００７４】
スコアのセットはすべてのアダマール符号ワードについて計算される。スコア関数を最大にする符号ワードｂが、推定された電子透かし埋め込み符号ワードｂ＊である。アダマール行列行のインデックスを使用して電子透かし埋め込み符号ワードｂ＊が復号され、推定復号メッセージｍ＊が取得される（ステップ７５）。留意されるべき点であるが、図６の（Ａ）および（Ｂ）に示された方法に従って、アダマール復号使用のＳｃ（ｂ）決定を容易にするため、係数は対角線平行斜線順序づけに従って配列および選択される。しかしなら、更に留意されるべき点であるが、ＭＬ復号の間に実行される確率計算を最小にするため、他の復号方式と組み合わされた他の変換順序づけを使用して電子透かし技法を同様に最適化することもできる。
【００７５】
以上、本発明の完全な理解を得るため、特定の変換、符号化、複合化および順序づけのタイプなど多数の特定の詳細が記述された。しかしながら、これら特定の詳細は本発明を実現する上でかならずしも必要とされない点は当業者に明らかであろう。また、周知の画像処理ステップは、本発明を不必要に不明瞭にすることを避けるため、詳細の記述はされてない。
【００７６】
更に、本発明の構成要素が特定の実施形態に関連して記述されているが、本発明を他の種々の形態で実施することが可能である点は認められることであろう。従って、上述された特定の実施形態は例示の目的に記述されたものであり、本発明をそれに制限するように意図されたものでない点は理解されるべきことである。
【図面の簡単な説明】
【図１】電子透かし埋め込みプロセスの基本段階を示すブロック図である。
【図２】本発明に従った変倍ビン符号化方式を使用した電子透かし埋め込みシステムおよび方法の第１の実施形態を示すブロック図である。
【図３】本発明に従った変倍ビン符号化方式の第１の実施形態を示す流れ図である。
【図４】本発明に従った確率型復号方式を使用した電子透かし埋め込みシステムおよび方法の第１の実施形態を示すブロック図である。
【図５】確率型復号方式を使用して電子透かしデータを復号する電子透かし埋め込み方法の１つの実施形態を示す流れ図である。
【図６】変倍ビン符号化方式および最尤法復号方式を使用する電子透かし埋め込み方法の１つの実施形態を示すブロック図である。
【図１】
2004507177.tif 000020
【図２Ａ】
2004507177.tif 000021
【図２Ｂ】
2004507177.tif 000022
【図３】
2004507177.tif 000023
【図４Ａ】
2004507177.tif 000024
【図４Ｂ】
2004507177.tif 000025
【図５】
2004507177.tif 000026
【図６Ａ】
2004507177.tif 000027
【図６Ｂ】
2004507177.tif 000028
【図６Ｃ】
2004507177.tif 000029
【国際公開パンフレット】
2004507177.tif 000030
2004507177.tif 000031
2004507177.tif 000032
2004507177.tif 000033
2004507177.tif 000034
2004507177.tif 000035
2004507177.tif 000036
2004507177.tif 000037
2004507177.tif 000038
2004507177.tif 000039
2004507177.tif 000040
2004507177.tif 000041
2004507177.tif 000042
2004507177.tif 000043
2004507177.tif 000044
2004507177.tif 000045
2004507177.tif 000046
2004507177.tif 000047
2004507177.tif 000048
2004507177.tif 000049
2004507177.tif 000050
2004507177.tif 000051
2004507177.tif 000052
2004507177.tif 000053
2004507177.tif 000054
2004507177.tif 000055
2004507177.tif 000056
2004507177.tif 000057
2004507177.tif 000058
2004507177.tif 000059
2004507177.tif 000060
2004507177.tif 000061
2004507177.tif 000062
【国際公開パンフレット（コレクトバージョン）】
2004507177.tif 000063
2004507177.tif 000064
【国際公開パンフレット（コレクトバージョン）】
2004507177.tif 000065
2004507177.tif 000066
2004507177.tif 000067
2004507177.tif 000068
2004507177.tif 000069
2004507177.tif 000070
2004507177.tif 000071
2004507177.tif 000072
2004507177.tif 000073
2004507177.tif 000074
2004507177.tif 000075
2004507177.tif 000076
2004507177.tif 000077
2004507177.tif 000078
2004507177.tif 000079
2004507177.tif 000080
2004507177.tif 000081
2004507177.tif 000082
2004507177.tif 000083
2004507177.tif 000084
2004507177.tif 000085
2004507177.tif 000086
2004507177.tif 000087
2004507177.tif 000088
2004507177.tif 000089
2004507177.tif 000090
2004507177.tif 000091
2004507177.tif 000092
2004507177.tif 000093
2004507177.tif 000094
2004507177.tif 000095
2004507177.tif 000096
2004507177.tif 000097
2004507177.tif 000098
【国際調査報告】
2004507177.tif 000099
2004507177.tif 000100