(57)【要約】【課題】電子透かしをオリジナルの処理対象情報から生成して利用することができると共に、処理対象情報への改竄を高精度に検知することができる情報処理装置等を提供する。【解決手段】周波数分解が可能な処理対象情報に対して周波数変換を行う変換処理部１０と、生成された複数の高周波成分ごとに当該各成分を複数のブロックに区画し、区画された複数のブロック領域ごとにグループ化し、各グループにおいて所定の条件を満たすブロックを特定するブロック特定部１１と、特定された前記ブロックに所定の値を設定して透かし成分を生成する透かし生成部１２と、生成された前記透かし成分を前記周波数変換された任意の１又は複数の成分に埋め込む透かし埋込部１３と、前記変換処理部１０の逆変換を行って前記透かし成分が埋め込まれた前記処理対象物を生成する逆変換処理部１４とを備える。【選択図】図５
2014175852.tif 000002
【特許請求の範囲】
【請求項１】
周波数分解が可能な処理対象情報に対して周波数変換を行う変換手段と、
生成された複数の高周波成分ごとに当該各成分を複数のブロックに区画し、区画された複数のブロック領域ごとにグループ化し、各グループにおいて所定の条件を満たすブロックを特定するブロック特定手段と、
特定された前記ブロックに所定の値を設定して透かし成分を生成する透かし生成手段と、
生成された前記透かし成分を前記周波数変換された任意の１又は複数の成分に埋め込む透かし埋込手段と、
前記変換手段の逆変換を行って前記透かし成分が埋め込まれた前記処理対象物を生成する逆変換手段とを備えることを特徴とする情報処理装置。
【請求項２】
請求項１に記載の情報処理装置において、
前記ブロック特定手段が、前記各グループにおける成分値が最大となるブロックを特定することを特徴とする情報処理装置。
【請求項３】
請求項１又は２に記載の情報処理装置において、
前記透かし成分が埋め込まれた処理対象物に対して周波数変換を行う第２変換手段と、
前記透かし成分が埋め込まれた成分から当該透かし成分を抽出する透かし抽出手段と、
前記第２変換手段により生成された複数の高周波成分ごとに当該各成分を複数のブロックに区画し、区画された複数のブロック領域ごとにグループ化し、各グループにおいて最大値となるブロックを特定する推定ブロック特定手段と、
特定された前記ブロックに所定の値を設定して推定透かし成分を生成する推定透かし生成手段と、
前記透かし抽出手段により抽出された前記透かし成分と、前記推定透かし生成手段により生成された前記推定透かし成分との差分を演算する差分演算手段と、
前記差分演算手段による演算結果に基づいて、前記処理対象情報の改竄を特定する改竄特定手段とを備えることを特徴とする情報処理装置。
【請求項４】
請求項１ないし３のいずれかに記載の情報処理装置において、
処理対象物に対して周波数変換を行う場合に、区間演算による周波数変換を行うことを特徴とする情報処理装置。
【請求項５】
請求項１ないし４のいずれかに記載の情報処理装置において、
前記周波数変換の処理が、二重ツリー複素ウェーブレット変換又はダイアディックウェーブレット変換であることを特徴とする情報処理装置。
【請求項６】
周波数分解が可能な処理対象情報に対して周波数変換を行う変換手段と、
生成された複数の高周波成分に、生成された低周波成分を所定の定数倍して加算する低周波成分加算手段と、
前記周波数変換された低周波成分、又は、前記低周波成分が加算された前記高周波成分のいずれか１又は複数の成分に透かしを埋め込む透かし埋込手段と、
前記変換手段の逆変換を行って前記透かしが埋め込まれた前記処理対象物を生成する逆変換手段とを備えることを特徴とする情報処理装置。
【請求項７】
周波数分解が可能な処理対象情報に対して周波数変換を行う変換手段、
生成された複数の高周波成分ごとに当該各成分を複数のブロックに区画し、区画された複数のブロック領域ごとにグループ化し、各グループにおいて最大値となるブロックを特定するブロック特定手段、
特定された前記ブロックに所定の値を設定して透かし成分を生成する透かし生成手段、
生成された前記透かし成分を前記周波数変換された任意の１又は複数の成分に埋め込む透かし埋込手段、
前記変換手段の逆変換を行って前記透かし成分が埋め込まれた前記処理対象物を生成する逆変換手段としてコンピュータを機能させる情報処理プログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、処理対象情報に対して透かし成分を埋め込む情報処理装置等に関する。
【背景技術】
【０００２】
近年、画像編集ツール等の普及により、誰でも簡単に画像編集が可能である。そのため、例えば、ドライブレコーダなどの記録メディアで保存されたデジタルコンテンツの証拠能力が疑問視されている。デジタルコンテンツの証拠能力を保証するために、デジタルコンテンツに対する改竄を検知する技術を開発する必要がある。
【０００３】
電子透かし法に関する技術として、発明者らにより、例えば非特許文献１に示す技術が開示されている。非特許文献１に示す技術は、区間演算と二重ツリー複素数離散ウェーブレット変換に基づいた電子透かし法である。
【先行技術文献】
【特許文献】
【０００４】
【非特許文献１】大浦龍二，皆本晃弥、”二重ツリー複素離散ウェーブレット変換と区間演算に基づく電子透かし法”、電気関係学会九州支部連合大会講演論文集、２０１１年９月１６日、６４巻、Ｎｏ．３−１Ａ−１０
【発明の概要】
【発明が解決しようとする課題】
【０００５】
しかしながら、非特許文献１に示す技術は、改竄を検知する具体的な手段が記載されておらず、改竄検知に関して不十分な技術であると共に、電子透かしを事前に別途用意しなければならず、透かしの埋め込み及び認証の際に煩わしい処理が必要となってしまうという課題を有する。また、仮に改竄の有無を検知することができたとしても、その位置を正確に特定することは難しく誤検知も多く発生してしまう。
【０００６】
本発明は、電子透かしをオリジナルの処理対象情報から生成して利用することができると共に、処理対象情報への改竄を高精度に検知することができる情報処理装置等を提供する。
【課題を解決するための手段】
【０００７】
本発明に係る情報処理装置は、周波数分解が可能な処理対象情報に対して周波数変換を行う変換手段と、生成された複数の高周波成分ごとに当該各成分を複数のブロックに区画し、区画された複数のブロック領域ごとにグループ化し、各グループにおいて最大値となるブロックを特定するブロック特定手段と、特定された前記ブロックに所定の値を設定して透かし成分を生成する透かし生成手段と、生成された前記透かし成分を前記周波数変換された任意の１又は複数の成分に埋め込む透かし埋込手段と、前記変換手段の逆変換を行って前記透かし成分が埋め込まれた前記処理対象物を生成する逆変換手段とを備えるものである。
【０００８】
このように、本発明に係る情報処理装置においては、周波数変換されて得られた高周波成分から透かし成分を生成し、当該透かし成分をオリジナルの処理対象情報に埋め込むため、事前に電子透かしを用意する必要がなく、処理の手間を省くことができるという効果を奏する。
【０００９】
また、オリジナルの処理対象情報に基づいて透かし成分が生成されるため、透かし成分がオリジナルに馴染み易く、オリジナルの情報を壊し難くなるという効果を奏する。
【００１０】
本発明に係る情報処理装置は、前記ブロック特定手段が、前記各グループにおける成分値が最大となるブロックを特定するものである。
【００１１】
このように、本発明に係る情報処理装置においては、各グループにおける成分値が最大となるブロックを特定して透かし成分が生成されるため、オリジナルの情報に馴染み易い透かしを生成することができるという効果を奏する。
【００１２】
また、成分値が最大となるブロックを特定することで、改竄がされやすい箇所（例えば、画像の場合は空や海のような特徴がない平坦な箇所ではなく、特徴的な成分が多く含まれる箇所）に透かし成分が埋め込まれるため、改竄の検出精度を上げることができるという効果を奏する。
【００１３】
本発明に係る情報処理装置は、前記透かし成分が埋め込まれた処理対象物に対して周波数変換を行う第２変換手段と、前記透かし成分が埋め込まれた成分から当該透かし成分を抽出する透かし抽出手段と、前記第２変換手段により生成された複数の高周波成分ごとに当該各成分を複数のブロックに区画し、区画された複数のブロック領域ごとにグループ化し、各グループにおいて最大値となるブロックを特定する推定ブロック特定手段と、特定された前記ブロックに所定の値を設定して推定透かし成分を生成する推定透かし生成手段と、前記透かし抽出手段により抽出された前記透かし成分と、前記推定透かし生成手段により生成された前記推定透かし成分との差分を演算する差分演算手段と、前記差分演算手段による演算結果に基づいて、前記処理対象情報の改竄を特定する改竄特定手段とを備えるものである。
【００１４】
このように、本発明に係る情報処理装置においては、抽出した透かし成分と、埋め込みの際に生成した方法と同様の方法で得られた推定透かし成分とを比較し、その差分に応じて改竄を特定するため、事前に透かし画像等を用意することなく、確実に改竄を検知することができるという効果を奏する。
【００１５】
本発明に係る情報処理装置は、処理対象物に対して周波数変換を行う場合に、区間演算による周波数変換を行うものである。
【００１６】
このように、本発明に係る情報処理装置においては、処理対象物に対して周波数変換を行う場合に、区間演算による周波数変換を行うため、低周波成分を含む高周波成分を演算工程において自然に得ることができ、高画質を維持しつつ、ロバストな透かし成分の埋め込みが可能になるという効果を奏する。
【００１７】
本発明に係る情報処理装置は、前記周波数変換の処理が、二重ツリー複素ウェーブレット変換又はダイアディックウェーブレット変換であるものである。
【００１８】
このように、本発明に係る情報処理装置においては、二重ツリー複素ウェーブレット変換又はダイアディックウェーブレット変換を利用することで、回転やサイズ変更といった幾何学的な攻撃に対してもロバストな透かし成分の埋め込みが可能となり、改竄の有無と共に改竄箇所も正確に特定することができるという効果を奏する。
【００１９】
本発明に係る情報処理装置は、周波数分解が可能な処理対象情報に対して周波数変換を行う変換手段と、生成された複数の高周波成分に、生成された低周波成分を所定の定数倍して加算する低周波成分加算手段と、前記周波数変換された低周波成分、又は、前記低周波成分が加算された前記高周波成分のいずれか１又は複数の成分に透かしを埋め込む透かし埋込手段と、前記変換手段の逆変換を行って前記透かしが埋め込まれた前記処理対象物を生成する逆変換手段とを備えるものである。
【００２０】
このように、本発明に係る情報処理装置においては、高周波成分に低周波成分を加算することで、高周波成分に透かし成分を埋め込んだ場合であっても、高画質を維持し攻撃に対してロバストな埋め込み情報を得ることができるという効果を奏する。
【図面の簡単な説明】
【００２１】
【図１】ＤＴ−ＣＤＷＴの処理を示す模式図である。
【図２】ＩＡに基づくＤＴ−ＣＤＷＴの処理結果の一例を示す図である。
【図３】ＩＡを用いた場合と用いない場合との画質を比較した結果を示す図である。
【図４】ＤＷＴベースの透かし法とＤＴ−ＣＤＷＴベースの透かし法との比較結果を示す図である。
【図５】第１の実施形態に係る情報処理装置における透かしの埋め込み処理に関する機能ブロック図である。
【図６】第１の実施形態に係る情報処理装置における透かしの埋め込み処理に関する動作を示すフローチャートである。
【図７】ＷＴＭＭの処理を示す模式図である。
【図８】第１の実施形態に係る情報処理装置の認証処理に関する機能ブロック図である。
【図９】第１の実施形態に係る情報処理装置の認証処理に関する動作を示すフローチャートである。
【図１０】実験に用いた画像サンプルを示す図である。
【図１１】圧縮を施した透かし入り画像から透かしを抽出した場合の赤、緑、青、輝度、青の色差、赤の色差の各成分に関する認証結果を示す図である。
【図１２】いかなる操作も施していない場合の認証結果を示す図である。
【図１３】コントラストを調整した場合の認証結果を示す図である。
【図１４】画像の一部が削除された場合の認証結果を示す第１の図である。
【図１５】画像の一部が削除された場合の認証結果を示す第２の図である。
【図１６】画像の一部が追加された場合の認証結果を示す図である。
【図１７】改竄領域が６４ピクセル以下の場合の認証結果を示す図である。
【図１８】いかなる操作も施していない場合の認証結果及びコントラストを調整した場合の認証結果を示す図である。
【図１９】画像の一部が削除又は追加された場合の認証結果を示す図である。
【図２０】いかなる操作も施していない場合の認証結果を示す図である。
【図２１】コントラストを調整した場合の認証結果を示す図である。
【図２２】画像の一部が削除された場合の認証結果を示す図である。
【図２３】画像の一部が追加された場合の認証結果を示す図である。
【発明を実施するための最良の形態】
【００２２】
以下、本発明の実施の形態を説明する。本実施形態の全体を通して同じ要素には同じ符号を付けている。
【００２３】
（本発明の第１の実施形態）
本実施形態に係る情報処理装置について、図１ないし図９を用いて説明する。本実施形態においては、処理対象情報としてドライブレコーダで多く採用されているＪＰＥＧ画像を用い、ドライブレコーダに対する改竄を検知する場合について説明する。また、周波数変換を区間演算（ＩＡ：Interval  Arithmetic）に基づく二重ツリー複素離散ウェーブレット変換（ＤＴ−ＣＤＷＴ：Dual-Tree  Complex  Discrete  Wavelet  Transform）で行うものとする。
【００２４】
近年、ドライブレコーダシステムを搭載した車両が増えている。ドライブレコーダシステムは、大きくシングルカメラシステムとダブルカメラシステムとに分類される。シングルカメラシステムは、ドライバの正面の光景のみを視覚的に記録するのに対して、ダブルカメラシステムは、車内のドライバや搭乗者も記録する。ここでは、より広く普及しているシングルカメラシステムの情報を用いるものとする。例えば、衝突事故の場合、ドライブレコーダは、法的手段や保険補償に関してドライバを保護するために利用される。これらの記録情報が現実的な証拠として認められるためには、情報が改竄されていないことを証明する確実な方法が必要である。
【００２５】
なお、以下に詳細に説明する本発明に係る情報処理装置は、上記ドライブレコーダの記録情報以外にも、例えば著作権の保護や、その他デジタル情報の照合、認証に利用することが可能である。また、画像情報に限らず、周波数特定を示す音声情報等にも、本発明に係る情報処理装置を適用することができる。
【００２６】
本実施形態において用いるＤＴ−ＣＤＷＴは、完全シフト不変性（ＰＴＩ：Perfect  Translation  Invariance）を特徴とし、ウェーブレット係数のエネルギー分布が、空間領域における変化の影響を受けない。このため、改竄された部位を特定するのに有用である。また、多くのドライブレコーダシステムでは、動画を記録するのにモーションＪＰＥＧ（Ｍ−ＪＰＥＧ）フォーマットを採用している。Ｍ−ＪＰＥＧでは、それぞれのビデオフレームやデジタルビデオシーケンスのインターレース領域が個別にＪＰＥＧ画像として圧縮される。すなわち、これらの動画は静止画の連続で成り立っていることとなる。したがって、本実施形態においては、静止画像への透かしの埋め込み及び認証について考慮すればよい。
【００２７】
ＤＴ−ＣＤＷＴについて説明する。対象のデジタル信号｛ｆｌ｝は、実数及び虚数のスケーリング関数φＲ（ｔ−ｋ）及びφＩ（ｔ−ｋ）（ただしｋ∈Ｚ）を用いて、次式で表される。
【００２８】
【数１】
2014175852.tif 000003
【００２９】
ここで、関数ｆ（ｔ）は、対象のデジタル信号｛ｆｌ｝を補間する。すなわち、ｆｎ＝ｆ（ｎ）（ただし、ｎ∈Ｚ）。（￣φ（ｔ））（上線の表現は、このように括弧書きで括る。以下同様。）はφ（ｔ）の共役複素数、Ｚは整数の集合を示す。ここで、ＤＴ−ＣＤＷＴは次の分解アルゴリズムを用いて算出される。
【００３０】
【数２】
2014175852.tif 000004
【００３１】
ここで、｛ａｎＲ｝と｛ｂｎＲ｝は実数の分解シーケンス、｛ａｎＩ｝と｛ｂｎＩ｝は虚数の分解シーケンスである。｛ａｎＲ｝、｛ａｎＩ｝は低域フィルタ、｛ｂｎＲ｝、｛ｂｎＩ｝は広域フィルタとして機能する。上記シーケンスを取りまとめると、
（１）シーケンス｛ｃｊ，ｋＲ｝と｛ｃｊ，ｋＩ｝は、それぞれレベルｊの実数及び虚数のスケーリング係数である。
（２）シーケンス｛ｄｊ，ｋＲ｝と｛ｄｊ，ｋＩ｝は、それぞれレベルｊの実数及び虚数のウェーブレット係数である。
（３）ｊの値は周波数に対応する。より具体的には、その成分の周波数は、対応するｊの値が小さくなるにつれて大きくなる。この点から、ｊ＝０を設定し、それぞれの係数に対してレベル０と−１のみを考える。
（４）式（２．１）より、対象のデジタル信号｛ｆｌ｝からレベル０の実数及び虚数のスケーリング係数｛ｃ０，ｋＲ｝、｛ｃ０，ｋＩ｝を求める。
（５）式（２．１）で求めた｛ｃ０，ｋＲ｝、｛ｃ０，ｋＩ｝を用いて、レベル−１の実数及び虚数のスケーリング係数｛ｃ−１，ｋＲ｝、｛ｃ−１，ｋＩ｝、式（２．２）の分解アルゴリズムを用いてレベル−１の実数及び虚数のウェーブレット係数｛ｄ−１，ｋＲ｝、｛ｄ−１，ｋＩ｝を求める。
（６）２次元のＤＴ−ＣＤＷＴであるため、演算がより複雑になり、式（２．２）に記載の入力シーケンスは必ずしも｛ｃ０，ｋＲ｝、｛ｃ０，ｋＩ｝にはならない。演算処理については、図１を用いて後述する。
【００３２】
ＤＴ−ＣＤＷＴの逆変換は、次の再構成アルゴリズムによって算出される。
【００３３】
【数３】
2014175852.tif 000005
【００３４】
ここで、｛ｇｎＲ｝と｛ｈｎＲ｝は、実数の再構成シーケンス、｛ｇｎＩ｝と｛ｈｎＩ｝は、虚数の再構成シーケンスである。式（２．３）、（２．１）及びｎ∈Ｚにおけるｆｎ＝ｆ（ｎ）の関係から、元の離散的信号｛ｆｎ｝が求められる。
【００３５】
次に、ＩＡについて説明する。区間Ａは実数の集合Ｒの連結部分集合である。［ａ１，ａ２］で示される閉区間は、｛ｔ｜ａ１≦ｔ≦ａ２，ａ１，ａ２∈Ｒ｝で得られる実数の組である。本実施形態においては、「区間」という言葉は閉区間を指すこととする。区間Ａの下限及び上限をそれぞれｉｎｆ（Ａ）＝ａ１及びｓｕｐ（Ａ）＝ａ２で表し、いずれかの非空区間Ａは、ω（Ａ）＝ａ２−ａ１で定義される。２つの区間Ａ＝［ａ１，ａ２］及びＢ＝［ｂ１，ｂ２］について、基本的な四則演算は次のように定義される。
【００３６】
【数４】
2014175852.tif 000006
【００３７】
特に、式（２．４）のγ＝ｂ１＝ｂ２を設定する際のスカラー乗法は次式で定義される。
【００３８】
【数５】
2014175852.tif 000007
【００３９】
その成分が区間を構成する区間ベクトル及びマトリクスには、各成分にこれらの演算を当てはめる。
【００４０】
ただし、一般に区間演算を数値的に厳密な方式で行うことは不可能である。特に、丸め誤差のせいで、式（２．４）に定義された区間を正確に生成することは不可能である。例えば、任意のＡとＢをＡ＋Ｂ（／∈）（斜線の表現は、このように括弧書きで括る。以下同様。）ｆｌ（Ａ＋Ｂ）とするのは可能である。ここで、左辺は正確な合計、右辺は一般的な浮動小数点で数値的に算出したものである。数値的に算出された区間が実際の区間を含まないという状態を避けるために、計算する際に、区間の下端点を常に切り捨て、上端点を切り上げるという丸めルールを用いる。例えば、２つの区間を合わせるとき、ｎｕｍ（Ａ＋Ｂ）＝［▽（ａ１＋ｂ１），△（ａ２＋ｂ２）］として計算を行う。ただし、ｎｕｍ（Ａ＋Ｂ）は数値的に算出された合計を表し、▽及び△は切り捨て及び切り上げの演算をそれぞれ表す。すなわち、Ａ＋Ｂ⊂ｎｕｍ（Ａ＋Ｂ）の関係は常に満たされる。
【００４１】
式（２．４）から式（２．５）に定義された演算によると、一般に、区間の幅は計算した数に比例して大きくなる。この現象は区間拡張と呼ばれ、コンピュータを利用して数学的に厳密な処理を行いたい場合には、極めて不都合な要素であるが、対照的に本実施形態においては、オリジナルの情報から余剰部分を作り出す有用なツールとして利用できる。
【００４２】
ＩＡに基づくＤＴ−ＣＤＷＴについて説明する。ＩＡに基づくＤＴ−ＣＤＷＴを次式で定義する。
【００４３】
【数６】
2014175852.tif 000008
【００４４】
ここで、Ｉ（Δｋ）＝［１−Δｋ，１＋Δｋ］、Δｋは正の実数である。式（２．６）に含まれる全ての演算はＩＡを用いて行われるため、ｃｊ−１，ｎＲ∈Ｉ（ｃｊ−１，ｎＲ）、ｄｊ−１，ｎＲ∈Ｉ（ｄｊ−１，ｎＲ）、ｃｊ−１，ｎＩ∈Ｉ（ｃｊ−１，ｎＩ）、ｄｊ−１，ｎＩ∈Ｉ（ｄｊ−１，ｎＩ）が得られ、実際にその関係は、式（２．４）及び式（２．２）によって、
【００４５】
【数７】
2014175852.tif 000009
【００４６】
となる。すなわち、ＩＡに基づくＤＴ−ＣＤＷＴは、元のＤＴ−ＣＤＷＴを含んでいることとなる。つまり、ＩＡに基づくＤＴ−ＣＤＷＴは、元の信号から余剰の信号を作り出している。画像の場合、式（２．２）及び式（２．６）を各方向、すなわち垂直方向及び水平方向に当てはめる。この手順を図１に示す。図１において、Ｃ、Ｄ、Ｅ及びＦは、それぞれ垂直、水平、対角線方向の低周波成分及び高周波成分を示す。上付きの文字Ｒ及びＩは、それぞれ「実数」及び「虚数」を示す。例えば、図１のＤ−１ＲＩは、式（２．２）に記載の垂直方向の低域フィルタ｛ａｎＲ｝（実数の分解シーケンス）で求めた係数の集合、及び、水平方向の広域フィルタ｛ｂｎＩ｝（虚数の分解シーケンス）で求めた係数の集合を示す。式（２．６）に基づくＤ−１ＲＩに相当する区間の成分は、Ｉ（Ｄ−１ＲＩ）で示される。Ｄ−１ＲＩ（ｍ，ｎ）（０≦ｍ≦Ｎｘ／２，  ０≦ｎ≦Ｎｙ／２，  ｍ，ｎ∈Ｚ）は、（ｍ，ｎ）位置の係数を表す。ここで、ｍ、ｎはそれぞれ垂直、水平方向の位置、Ｎｘ、Ｎｙはそれぞれ画像の垂直、水平方向のサイズを示す。
【００４７】
デジタル画像透かし法は、透かし入り画像の堅牢性と画質とのトレードオフの関係がある。低周波成分は画像を表す重要な情報の大部分を含んでいるため、低周波成分に透かしを埋め込むと堅牢性が増すが、透かし入り画像を劣化させてしまう。画質を維持するためには、高周波成分に透かしを埋め込む必要があるが、画像の堅牢性が低減してしまう。ＩＡに基づくＤＴ−ＣＤＷＴにより得られた高周波成分は、事実上低周波成分を含んでいることがわかる。この特性により、画像の重要な情報の大部分を含んでいる低周波成分を高周波成分のみから得ることができる。
【００４８】
図２は、この状態を示している。図２（Ａ）はＤ−１ＲＩ、図２（Ｂ）はｓｕｐ（Ｉ（Ｄ−１ＲＩ））、図２（Ｃ）はｉｎｆ（Ｉ（Ｄ−１ＲＩ））、図２（Ｄ）はｓｕｐ（Ｉ（Ｄ−１ＲＩ））−ｉｎｆ（Ｉ（Ｄ−１ＲＩ））である。また、図２（Ｅ）はｓｕｐ（Ｉ（Ｄ−１ＲＲ））−ｉｎｆ（Ｉ（Ｄ−１ＲＲ））、図２（Ｆ）はｓｕｐ（Ｉ（Ｄ−１ＲＩ））−ｉｎｆ（Ｉ（Ｄ−１ＲＩ））、図２（Ｇ）はｓｕｐ（Ｉ（Ｄ−１ＩＲ））−ｉｎｆ（Ｉ（Ｄ−１ＩＲ））、図２（Ｈ）はｓｕｐ（Ｉ（Ｄ−１ＩＩ））−ｉｎｆ（Ｉ（Ｄ−１ＩＩ））である。図２に示すように、例えば、ｓｕｐ（Ｉ（Ｄ−１ＲＩ））とｉｎｆ（Ｉ（Ｄ−１ＲＩ））は、どちらも少量の低周波成分を含んでおり、ｓｕｐ（Ｉ（Ｄ−１ＲＩ））−ｉｎｆ（Ｉ（Ｄ−１ＲＩ））の値を正規化することで、これらの低周波成分を抽出することができる。
【００４９】
実際に、ＩＡに基づくＤＴ−ＣＤＷＴを用いた透かし法は、画像圧縮に対してＩＡを用いないＤＴ−ＣＤＷＴに比べて堅牢である。図３はこの状態を示している。図３（Ａ）は１８．８％の圧縮率でＪＰＥＧ圧縮した透かし入り画像及び抽出した透かし、図３（Ｂ）は１８．８８％の圧縮率でＪＰＥＧ圧縮した透かし入り画像及び抽出した透かしを示しており、図３（Ａ）はＩＡを用いないＤＴ−ＣＤＷＴベースの透かし法であり、図３（Ｂ）はＩＡを用いたＤＴ−ＣＤＷＴベースの透かし法である。図３（Ｃ）は１２．６５％の圧縮率でＪＰＥＧ２０００圧縮した透かし入り画像及び抽出した透かし、図３（Ｄ）は１２．６１％の圧縮率でＪＰＥＧ２０００圧縮した透かし入り画像及び抽出した透かしを示しており、図３（Ｃ）はＩＡを用いないＤＴ−ＣＤＷＴベースの透かし法であり、図３（Ｄ）はＩＡを用いたＤＴ−ＣＤＷＴベースの透かし法である。
【００５０】
図３で実証されているように、二値の透かしには「Ｋ」の文字が４つ表示されており、ＩＡを用いた場合の透かしの画質は、ＩＡを用いない場合の透かしの画質に比べて優れていることが明らかである。
【００５１】
ここで、離散ウェーブレット変換（ＤＷＴ：Discrete  Wavelet  transform）とＤＴ−ＣＤＷＴの相違について述べる。ＤＷＴに基づく透かし法は数多く提案されている。このアプローチの難点は、シフト不変性がないことである。この特性のせいで、ＤＷＴベースの透かし法は、回転やサイズ変更などの幾何学的な攻撃に対してロバストではない。これに対して、ＤＴ−ＣＤＷＴはシフト不変であるため、幾何学的な攻撃に対してロバストである。
【００５２】
図４に、ＤＷＴベースの透かし法とＤＴ−ＣＤＷＴベースの透かし法との比較を示す。図４（Ａ）は３０度回転した場合の透かし入り画像から抽出した透かし、図４（Ｂ）は０．９のスケーリングで縮小した場合の透かし入り画像から抽出した透かしであり、いずれもＤＷＴベースの透かし法を用いている。図４（Ｃ）は３０度回転した場合の透かし入り画像から抽出した透かし、図４（Ｄ）は０．９のスケーリングで縮小した場合の透かし入り画像から抽出した透かしであり、いずれもＤＴ−ＣＤＷＴベースの透かし法を用いている。ＤＷＴベースの透かし法の場合は、抽出した透かしが大きくダメージを受けており、透かしの存在をほとんど確認することができない。これに対して、ＤＴ−ＣＤＷＴベースの透かし法の場合は、抽出した透かしが多少劣化しているものの、透かしの存在は一目で確認することができる。画像認証には改竄されたデータの区域を正確に検出できることが重要であるから、認証のための透かし法を開発する際には、ＤＴ−ＣＤＷＴがＤＷＴより優れている。
【００５３】
次に、本実施形態に係る情報処理装置のアルゴリズムを説明する。まず、透かしの埋め込み処理について説明する。図５は、本実施形態に係る情報処理装置の埋め込み処理に関する機能ブロック図である。情報処理装置１は、オリジナル画像を通常の方法でＤＴ−ＣＤＷＴ演算を行う第１変換部１０ａ、及び、オリジナル画像をＩＡを用いてＤＴ−ＣＤＷＴ演算を行う第２変換部１０ｂからなる変換処理部１０と、第１変換部１０ａで生成された複数の高周波成分ごとに各成分を、例えば８×８ピクセルのブロックに分割し、分割された複数のブロック領域ごと（例えば、隣接する４ブロックごと）にグループ化し、各グループにおいて最大となるブロックを特定するブロック特定部１１と、特定されたブロックに１又は−１を設定して二値の透かしを生成する透かし生成部１２と、第２変換部１０ｂで生成された１６の区間成分のうち、１又は複数の区間成分に透かしを埋め込む透かし埋込部１３と、逆ＤＴ−ＣＤＷＴを行って画像を再構成する逆変換処理部１４とを備える。
【００５４】
次に埋め込み処理の動作について説明する。図６は、本実施形態に係る情報処理装置の埋め込み処理に関する動作を示すフローチャートである。まず、第１変換部１０ａが、オリジナル画像に対して通常の方法でＤＴ−ＣＤＷＴを適用して、図１に示す１６の成分（Ｃ−１ＲＲ，Ｄ−１ＲＲ，・・・，Ｆ−１ＩＩ）を得る（Ｓ１１）。ブロック特定部１１が、Ｓ１１で得られた各成分を８×８ピクセルのオーバーラップしないブロックに分割し、以下の演算によりＷＴＭＭ（Wavelet  Transform  Modulus  maxima）を求める。
【００５５】
【数８】
2014175852.tif 000010
【００５６】
式（３．２）に示すように、平方根の項は１２の高周波成分からなる。各ブロックのＷＴＭＭ（ｍ，ｎ）の合計を算出し、図７に示すように４つの隣接するブロック（非オーバーラップ）のうち、成分値が最大となるブロックを特定する（Ｓ１２）。透かし生成部１２が、特定されたブロックにＷ（ｍ，ｎ）＝１又はＷ（ｍ，ｎ）＝−１を設定して透かしＷを生成する（Ｓ１３）。第２変換部１０ｂが、ＩＡを用いたＤＴ−ＣＤＷＴをオリジナル画像に適用して、図１のＣ−１ＲＲ，Ｄ−１ＲＲ，・・・，Ｆ−１ＩＩに相当する１６の区間成分を求める（Ｓ１４）。これらは、Ｉ（Ｃ−１ＲＲ），Ｉ（Ｄ−１ＲＲ），・・・，Ｉ（Ｆ−１ＩＩ）で表される。１６の区間成分から、透かしを埋め込む１又は複数の成分を選び、Ｉ（Ｓｉ）（ｉ＝１，２，・・・，１６，  １≦Ｎ≦１６）とする（Ｓ１５）。その量Ｓ’ｉは、区間Ｉ（Ｓｉ）内の点であればどこでもよいが、ここでは仮にＳ’ｉ＝ｓｕｐ（Ｉ（Ｓｉ））で定義する。他の区間成分を浮動小数点に置き換える。サイズ（２Ｋ＋１）×（２Ｌ＋１）のブロック（ブロック形式のスライディング・ウインドウ）を用いて、以下の演算により平滑化を行う（Ｓ１６）。
【００５７】
【数９】
2014175852.tif 000011
【００５８】
ここで、Ｋ及びＬは固定の自然数であり、ｓｇｎ（ａ）は実数ａの通常の符号関数である。以下の演算により、透かしを埋め込む（Ｓ１７）。
【００５９】
【数１０】
2014175852.tif 000012
【００６０】
なお、０＜α＜１は、堅牢性を調整する因子となる。~Ｓｉ及び浮動小数点に置き換えた成分を用いて、逆ＤＴ−ＣＤＷＴにより画像を再構成し（Ｓ１８）、透かし入り画像~Ｃ０を得る。
【００６１】
削除、追加、差し替えなどの悪意のある操作を検出できることは極めて重要である。また、ＪＰＥＧやＪＰＥＧ２０００圧縮、コントラスト調整などの悪意のない操作を認証できることも必要である。全ての高周波成分には、形状や輪郭など画像の重要な特徴に関する情報を含んでいるため、本実施形態においては全ての高周波成分を用いて透かしを生成している。悪意のある操作は、高周波成分に対して変化を加えると考えられるため、本実施形態の透かしを用いて、それらを簡単に検出することが可能となる。
【００６２】
本実施形態に係る情報処理装置では、認証において必要となる情報は透かしを埋め込んだ成分に関する情報だけである。この情報は、透かしを抽出するのに不可欠となる。以下、認証処理について説明する。図８は、本実施形態に係る情報処理装置の認証処理に関する機能ブロック図である。情報処理装置１は、透かし入り画像を通常の方法でＤＴ−ＣＤＷＴ演算を行う変換処理部１０と、変換された成分から透かし画像を抽出する透かし抽出部２１と、埋め込みの処理の場合と同様の方法で、ブロック分割及びグループ化を行って、最大となるブロックを特定する推定ブロック特定部２２と、埋め込み処理の場合と同様の方法で、特定されたブロックに１又は−１を設定して推定透かしを生成する推定透かし生成部２３と、抽出された透かしと推定透かしとを比較して差分を演算する差分演算部２４と、差分が演算された結果に基づいて、改竄の部位を特定する改竄特定部２５とを備える。
【００６３】
次に認証処理の動作について説明する。図９は、本実施形態に係る情報処理装置の認証処理に関する動作を示すフローチャートである。まず、変換処理部１０が、透かし入り画像~Ｃ０に対して通常の方法でＤＴ−ＣＤＷＴを適用して、図１に示す１６の成分（~Ｃ−１ＲＲ，~Ｄ−１ＲＲ，・・・，~Ｆ−１ＩＩ）を得る（Ｓ２１）。透かし抽出部２１が、透かしを入れた成分~Ｓｉから、符号関数を用いて、~Ｗｉ＝ｓｇｎ（｜~Ｓｉ｜−｜~Ｓｉ｜）及び~Ｗｅ＝Σｉ＝１Ｎ~Ｗｉを計算し、~Ｗｅを閾値化して二値の透かし~Ｗを抽出する（Ｓ２２）。推定ブロック特定部２２が、Ｓ２１で得られた各成分を８×８ピクセルのオーバーラップしないブロックに分割し、以下の演算によりＷＴＭＭを求める。
【００６４】
【数１１】
2014175852.tif 000013
【００６５】
各ブロックのＷＴＭＭ’（ｍ，ｎ）の合計を算出し、４つの隣接するブロック（非オーバーラップ）のうち、成分値が最大となる推定ブロックを特定する（Ｓ２３）。推定透かし生成部２３が、特定された推定ブロックにＷ’（ｍ，ｎ）＝１又はＷ’（ｍ，ｎ）＝−１を設定して推定透かしＷ’を生成する（Ｓ２４）。差分演算部２４が、抽出した透かし~Ｗと推定透かしＷ’との差分を演算し（Ｓ２５）、改竄特定部２５が、選択したブロックにおいて、閾値γ及び前述の差分ｄ＝Σｍ，ｎ｜Ｗ’（ｍ，ｎ）−~Ｗ（ｍ，ｎ）｜を用いて、改竄された部位を特定する（Ｓ２６）。選択したブロックから透かしが抽出できない場合は、ｄの値が大きくなる傾向がある。そこで、選択したブロックを含む隣接したブロックの集合は改竄されたものであると判断する。
【００６６】
このように、本実施形態に係る情報処理装置においては、オリジナル画像から透かし成分を生成し、当該透かし成分をオリジナル画像に埋め込むため、事前に透かしを用意する必要がなく、処理の手間を省くことができる。また、オリジナルの画像に基づいて透かし成分が生成されるため、透かし成分がオリジナルに馴染み易く、オリジナルの情報を壊し難くなる。さらに、各グループにおける成分値が最大となるブロックを特定して透かし成分が生成されるため、オリジナルの情報に馴染み易い透かしを生成することができる。さらにまた、成分値が最大となるブロックを特定することで、改竄がされやすい箇所（例えば、画像の場合は空や海のような特徴がない平坦な箇所ではなく、特徴的な成分が多く含まれる箇所）に透かし成分が埋め込まれるため、改竄の検出精度を上げることができる。
【００６７】
（本発明の第２の実施形態）
本実施形態に係る情報処理装置について説明する。本実施形態におては、第１の実施形態で用いたＤＴ−ＣＤＷＴの代わりに、ダイアディックウェーブレット変換（ＤＹＷＴ：Dyadic  Wavelet  Transform）を用いるものである。
なお、本実施形態において前記第１の実施形態と重複する説明は省略する。
【００６８】
ＩＡに基づくＤＹＷＴについて説明する。Ｃ０［ｍ，ｎ］をオリジナルの画像とする。次式で示すＤＹＷＴが一般的によく知られている。
【００６９】
【数１２】
2014175852.tif 000014
【００７０】
ここで、Ｃｋ，ｌｊ［ｍ，ｎ］＝Ｃｊ［ｍ＋２ｊｋ，ｎ＋２ｊｌ］，ｈは、ローパスフィルタ、ｇはハイパスフィルタである。より正確には、Ｃｊ［ｍ，ｎ］，Ｄｊ［ｍ，ｎ］，Ｅｊ［ｍ，ｎ］，Ｆｊ［ｍ，ｎ］は、それぞれ低周波成分、水平方向の高周波成分、垂直方向の高周波成分、及び、両方向の高周波成分である。ｍ及びｎは、それぞれ水平方向、垂直方向の位置を示す。また、逆ＤＹＷＴは、以下のアルゴリズムによって演算される。
【００７１】
【数１３】
2014175852.tif 000015
ここで、
【００７２】
【数１４】
2014175852.tif 000016
【００７３】
であり、~ｈは双対ローパスフィルタで、~ｇが双対ハイパスフィルタである。式（４．１）及び式（４．２）を用いてオリジナルの画像Ｃ０が得られる。次に、式（２．４）及び式（４．１）に基づいて、以下のように２次元のＩＡに基づくＤＹＷＴを定義付ける。
【００７４】
【数１５】
2014175852.tif 000017
ここで、
【００７５】
【数１６】
2014175852.tif 000018
【００７６】
であり、δｋ,δｌは正の実数である。演算は全てＩＡを用いて実行され、Ｃｊ＋１⊂ＩＣｊ＋１、Ｄｊ＋１⊂ＩＤｊ＋１、Ｅｊ＋１⊂ＩＥｊ＋１及びＦｊ＋１⊂ＩＦｊ＋１の関係が保たれる。すなわち、ＩＡに基づくＤＹＷＴで得られた成分は、通常のＤＹＷＴで得られた成分を含む。
【００７７】
次に、アルゴリズムについて説明する。アルゴリズムは基本的には、前記第１の実施形態と同じである。異なるのは、ＤＴ−ＣＤＷＴに代わってＤＹＷＴを用いているため、分解される周波数成分が図１に示すような１６成分ではなく、上記に示したＣ、Ｄ、Ｅ及びＦの４つの成分となる。すなわち、ＷＴＭＭを求める際には、高周波成分であるＤ、Ｅ及びＦの３つの項から求めることとなる。また、透かしを埋め込む成分もＣ、Ｄ、Ｅ及びＦのうちの１又は複数の成分から選ばれることとなる。その他具体的な処理については、第１の実施形態におけるアルゴリズムと同じである。
【００７８】
このように、本実施形態に係る情報処理装置においては、ＤＹＷＴを用いて変換処理を行うため、オリジナルの画像サイズと分解後の画像サイズが同値となり、より細かい透かしの生成及び埋め込みが可能となる。すなわち、改竄箇所の検知をより細かい範囲で特定することが可能となり、ピンポイントに小さな改竄を特定することができる。
【００７９】
なお、上記第１の実施形態及び第２の実施形態においては、ＩＡを用いることで高周波成分に低周波成分が含まれるようにしたが、ＩＡを用いずに、意図的に高周波成分に低周波成分の定数倍の成分を加算するようにしてもよい。そうすることで、ＩＡによる弊害（区間拡張やオリジナル画像の劣化）等を防止して、高画質且つ堅牢性を維持しながら、簡単な処理で透かしの埋め込み及び認証を行うことが可能となる。
【００８０】
また、上記第１の実施形態及び第２の実施形態に係る情報処理装置の処理は、ＲＯＭやＲＡＭ等のメモリに記憶されたプログラムをＣＰＵが読み込んで実行することで実現される。また、情報処理装置１は、必要に応じてハードディスクやフラッシュメモリ等の大容量の記憶装置にデータを保存してもよい。さらに、情報処理装置１は、利用者からの操作情報等を入力したり、情報を出力するための入出力インターフェースを備えると共に、外部機器との通信を行うための通信インターフェースを備える構成としてもよい。
【実施例】
【００８１】
本発明に係る情報処理装置の能力を評価するために、以下の実験を行った。本実験には、ドライブレコーダで記録した１８０×２８８ピクセルの２４ビットカラー画像を採用した。周波数変換はＩＡを用いて行った。実験においては、透かしを４つの成分Ｉ（Ｄ−１ＲＲ）、Ｉ（Ｄ−１ＲＩ）、Ｉ（Ｄ−１ＩＲ）、Ｉ（Ｄ−１ＩＩ）に埋め込んだ。図１０に、オリジナル画像（上段）、ＷＴＭＭ画像（中段）、生成した透かし（下段）を示す。ＲＧＢ（Ｒｅｄ：赤、Ｇｒｅｅｎ：緑、Ｂｌｕｅ：青）モデルを採用し、透かしをＧ成分に埋め込んだ。
【００８２】
本実験では、式（２．１）に示した手順は省略した。すなわち、式（２．１）を用いる代わりに、ｃ０，ｋＲ＝ｆｋ／２及びｃ０，ｋＩ＝ｆｋ／２を設定し、元のフィルタ長は８１だが、そのうち主要な３６フィルタだけを用いて演算時間を短縮した。ＤＴ−ＣＤＷＴは完全シフト不変性を持たないが、おおよそのシフト不変性を持っており、以下のような良好な結果を得ることができた。
【００８３】
パラメータ値は、ピーク信号対雑音比（ＳＮ比：ＰＳＮＲ）の値が著しく低下しないように選択した。ＳＮ比は以下の式で求められる。
【００８４】
【数１７】
2014175852.tif 000019
【００８５】
ここで、Ｎｘ及びＮｙは、それぞれ水平方向及び垂直方向の画像サイズであり、Ｃ０（ｉ，ｊ，ｋ）と~Ｃ０（ｉ，ｊ，ｋ）における変数ｋは、ＲＧＢの３成分に対応する数を表している。
【００８６】
本実験では、３５を上回るＰＳＮＲを必要とし、α及びΔｋは、ＪＰＥＧ及びＪＰＥＧ２０００圧縮に対する堅牢性を最大化するため、できる限り大きくした。以下の表１〜４は、本実験で使用するパラメータを決めるための予備の実験結果を示す。表１〜４は、赤（Ｒ）、緑（Ｇ）、青（Ｂ）、輝度（Ｙ）、青の色差（Ｃｒ）、赤の色差（Ｃｂ）の各成分に関するＰＳＮＲを示しており、表１はＫ＝Ｌ＝１０、Δｋ＝０．００６、表２はＫ＝Ｌ＝５、Δｋ＝０．００６、表３はＫ＝Ｌ＝１０、α＝０．９、表４はＫ＝Ｌ＝１０、α＝０．７の場合の結果を示している。
【００８７】
【表１】
2014175852.tif 000020
【００８８】
【表２】
2014175852.tif 000021
【００８９】
【表３】
2014175852.tif 000022
【００９０】
【表４】
2014175852.tif 000023
【００９１】
表１及び表２によると、Ｋ及びＬはＰＳＮＲに目立った影響を及ぼさない。一方、表３及び表４から、Δｋ及びαがＰＳＮＲに強い影響を与えていることを示している。これらの結果にしたがって、Δｋ＝０．００６、Ｋ＝Ｌ＝１０、α＝０．９として実験を行った。
【００９２】
次に、表１〜４によると、透かしをどの成分に埋め込んだとしてもＰＳＮＲはごく僅かしか変わらないことがわかる。また、図１１に、圧縮を施した透かし入り画像から透かしを抽出した場合の赤（Ｒ）、緑（Ｇ）、青（Ｂ）、輝度（Ｙ）、青の色差（Ｃｒ）、赤の色差（Ｃｂ）の各成分に関する認証結果を示す。図１１（Ａ）はＪＰＥＧ圧縮（ＱＦ＝６０）、図１１（Ｂ）はＪＰＥＧ２０００圧縮（ＱＦ＝７）の場合である。これらの結果から、圧縮攻撃に対する堅牢性を実現するために、Ｇ成分に透かしを埋め込むことで、ＪＰＥＧ圧縮及びＪＰＥＧ２０００圧縮（悪意のない通常の操作）が施された画像を通過させられることを示している。さらに、１つの８×８ピクセルのブロックが２５％を超えて改竄されている場合に改竄された部位を特定できるように、閾値γ＝３２とした。
【００９３】
図１２にいかなる操作も施していない場合の認証結果、図１３にコントラストを調整した場合の認証結果を示す。それぞれの図において、図１２（Ａ）及び図１３（Ａ）は透かし入り画像、図１２（Ｂ）及び図１３（Ｂ）は抽出した透かし、図１２（Ｃ）及び図１３（Ｃ）は推定透かし、図１２（Ｄ）及び図１３（Ｄ）は差分画像、図１２（Ｅ）はいかなる操作も施していない場合の認証結果の画像、図１３（Ｅ）はコントラストを調整した場合の認証結果の画像を示している。これらはどの部分も改竄されていないことを示している。図１２及び図１３に示すように、コントラスト調整のような悪意のない操作については、差分画像に差異が表れず、認証を通過することが可能であることが示された。
【００９４】
図１４及び図１５に画像の一部が削除された場合の認証結果、図１６に画像の一部が追加された場合の認証結果を示す。図１４において、上段は交通標識が削除された場合、中段は建物が削除された場合、下段は人と自動車が削除された場合の認証結果を示しており、右側の画像の斜線部分は攻撃された部位の位置を示している。図１５において、上段は近づいてくる車両が削除された場合、中段は街灯が削除された場合、下段は道路脇の溝が削除された場合の認証結果を示しており、右側の画像の斜線部分は攻撃された部位の位置を示している。図１６において、上段は歩行者が追加された場合、中段は横断歩道が追加された場合、下段は歩行者が追加された場合の認証結果を示しており、右側の画像の斜線部分は攻撃された部位の位置を示している。
【００９５】
図１４ないし図１６で示されるように、本発明の情報処理装置では、悪意のある攻撃を確実に検出し、且つ、その攻撃部位を正確に特定できることが明らかである。
【００９６】
図１７に改竄領域が６４ピクセル以下の場合の認証結果を示す。すなわち、本実験においては、改竄される部位が１ブロック８×８ピクセルより大きいことを想定してブロックを設定しているため、８×８ピクセルよりも小さいサイズの改竄について検討した。図１７において、上段は改竄領域が６４ピクセル、中段は改竄領域が４０ピクセル、下段は改竄領域が５６ピクセルの場合を示す。図１７に示されるように、改竄領域が１ブロック未満のサイズであっても、改竄の有無及び改竄部位を正確に特定することができることが明らかとなった。
【００９７】
なお、この１ブロックのサイズ設定は、使用態様に応じて任意に設定することができる。例えば、コンテンツの利用者に透かしの存在を認識させたい場合には、ブロックのサイズを大きく設定することで、透かしが入っていることを意識させることができる。また、逆に、透かしによる画質の劣化を防止したい場合には、ブロックサイズを小さく設定することで、透かしを意識することなくコンテンツを利用することができる。
【００９８】
次に、ＤＹＷＴを行った場合の実験結果を示す。表５ないし表７に本実験で使用するパラメータを決めるための予備の実験結果を示す。表５〜７は、赤（Ｒ）、緑（Ｇ）、青（Ｂ）の各成分に関するＰＳＮＲを示しており、表５はＫ＝Ｌ＝１０、γ＝３、表６はＫ＝Ｌ＝２０、γ＝３、表７はＫ＝Ｌ＝２０、γ＝２の場合の結果を示している。
【００９９】
【表５】
2014175852.tif 000024
【０１００】
【表６】
2014175852.tif 000025
【０１０１】
【表７】
2014175852.tif 000026
【０１０２】
表５及び表６によると、Ｋ及びＬはＰＳＮＲに目立った影響を及ぼさない。一方、表７は、δｋｌがＰＳＮＲに影響を及ぼすことが示されている。これらの結果にしたがって、γ＝３のときはδ＝δｋ＝δｌ＝０．０２００、Ｋ＝Ｌ＝２０とし、γ＝２のときはδ＝δｋ＝δｌ＝０．０３５０、Ｋ＝Ｌ＝２０とした。また、上記と同様の理由で、圧縮攻撃に対する堅牢性を実現するために、Ｇ成分に透かしを埋め込む。
【０１０３】
図１８にいかなる操作も施していない場合の認証結果及びコントラストを調整した場合の認証結果を示す。図１８（Ａ）、（Ｅ）は透かし入り画像、図１８（Ｂ）、（Ｆ）は抽出した透かし、図１８（Ｃ）、（Ｇ）は推定透かし、図１８（Ｄ）はいかなる操作も施していない場合の認証結果の画像、図１８（Ｈ）はコントラストを調整した場合の認証結果の画像を示している。図１８に示すように、コントラスト調整のような悪意のない操作については、認証を通過することが可能であることが示された。
【０１０４】
図１９に画像の一部が削除又は追加された場合の認証結果を示す。図１９（Ａ）、（Ｅ）は透かし入り画像、図１９（Ｂ）、（Ｆ）は抽出した透かし、図１９（Ｃ）、（Ｇ）は推定透かし、図１９（Ｄ）は街灯が削除された場合の認証結果の画像、図１９（Ｈ）は歩行者が追加された場合の認証結果の画像を示している。なお、右側の画像の斜線部分は攻撃された部位の位置を示している。これらの結果から、悪意のある改竄に対して、その改竄の有無及び改竄位置を正確に特定することが可能であることが示された。
【０１０５】
図２０ないし図２３に、改竄がない場合とある場合との認証結果を示す。図２０にいかなる操作も施していない場合の認証結果、図２１にコントラストを調整した場合の認証結果を示す。それぞれの図において、上段は透かし入り画像、下段は抽出した透かしと推定透かしとの差分画像である。いずれの場合も、悪意のない操作で認証を通過することができる。図２２に標識を削除した場合の認証結果、図２３に横断歩道を追加した場合の認証結果を示す。それぞれの図において、上段は透かし入り画像、下段は抽出した透かしと推定透かしとの差分画像である。図２２及び図２３から、上記のＤＴ−ＣＤＷＴを用いた場合と比べて、全体的に改竄領域をより正確に細かく特定していることがわかる。ＤＹＷＴによって得られた各周波数成分のサイズは、オリジナルの画像サイズと同じであるため、ＤＴ−ＣＤＷＴと比べてより小さいサイズのブロックを用いることができる。すなわち、改竄領域をより正確に細かく特定することが可能となる。
【符号の説明】
【０１０６】
１  情報処理装置
１０  変換処理部
１０ａ  第１変換部
１０ｂ  第２変換部
１１  ブロック特定部
１２  透かし生成部
１３  透かし埋込部
１４  逆変換処理部
２１  透かし抽出部
２２  推定ブロック特定部
２３  推定透かし生成部
２４  差分演算部
２５  改竄特定部
【図１】
2014175852.tif 000027
【図５】
2014175852.tif 000028
【図６】
2014175852.tif 000029
【図７】
2014175852.tif 000030
【図８】
2014175852.tif 000031
【図９】
2014175852.tif 000032
【図１１】
2014175852.tif 000033
【図２】
2014175852.tif 000034
【図３】
2014175852.tif 000035
【図４】
2014175852.tif 000036
【図１０】
2014175852.tif 000037
【図１２】
2014175852.tif 000038
【図１３】
2014175852.tif 000039
【図１４】
2014175852.tif 000040
【図１５】
2014175852.tif 000041
【図１６】
2014175852.tif 000042
【図１７】
2014175852.tif 000043
【図１８】
2014175852.tif 000044
【図１９】
2014175852.tif 000045
【図２０】
2014175852.tif 000046
【図２１】
2014175852.tif 000047
【図２２】
2014175852.tif 000048
【図２３】
2014175852.tif 000049