JP2001119559A - 画像処理装置及び方法及び記憶媒体 - Google Patents

画像処理装置及び方法及び記憶媒体

Info

Publication number
JP2001119559A
JP2001119559A JP29816199A JP29816199A JP2001119559A JP 2001119559 A JP2001119559 A JP 2001119559A JP 29816199 A JP29816199 A JP 29816199A JP 29816199 A JP29816199 A JP 29816199A JP 2001119559 A JP2001119559 A JP 2001119559A
Authority
JP
Japan
Prior art keywords
image data
information
embedding
mask
embedded
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
JP29816199A
Other languages
English (en)
Other versions
JP3869983B2 (ja
Inventor
Junichi Hayashi
淳一 林
Tomochika Murakami
友近 村上
Keiichi Iwamura
恵市 岩村
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Canon Inc
Original Assignee
Canon Inc
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Canon Inc filed Critical Canon Inc
Priority to JP29816199A priority Critical patent/JP3869983B2/ja
Priority to US09/688,108 priority patent/US6826290B1/en
Publication of JP2001119559A publication Critical patent/JP2001119559A/ja
Application granted granted Critical
Publication of JP3869983B2 publication Critical patent/JP3869983B2/ja
Anticipated expiration legal-status Critical
Expired - Fee Related legal-status Critical Current

Links

Landscapes

  • Image Analysis (AREA)
  • Image Processing (AREA)
  • Editing Of Facsimile Originals (AREA)
  • Facsimile Image Signal Circuits (AREA)

Abstract

(57)【要約】 【課題】 電子透かし情報が埋め込まれた画像に幾何的
な歪みが与えられた場合においても、この歪みを的確に
判断できる様にし、結果的に後段の電子透かし情報の抽
出を正確に行える様にする。 【解決手段】 画像データの幾何学的な状態を示すレジ
ストレーション信号、及び前記画像データに電子透かし
情報が予め人間に目視しにくく埋め込まれているはずの
画像データから、該電子透かし情報を抽出する為の画像
処理装置であって、前記画像データを周波数変換する変
換手段と、前記変換手段により得られる信号に微分処理
を施す微分手段と、 前記微分手段で得られた信号と、
前記レジストレーション信号に関する情報とに基づい
て、前記画像データの幾何学的な補正を行う補正手段と
を有する。

Description

【発明の詳細な説明】
【0001】
【発明の属する技術分野】本発明は、電子透かし情報が
埋め込まれた画像データからこの情報を正確に抽出でき
る為の処理を行う画像処理装置及び方法及びこの方法を
記憶した記憶媒体に関するものである。
【0002】
【従来の技術】近年のコンピュータ及びネットワークの
発達は著しく、文字データ、画像データ、音声データ等
のマルチメディアデータがコンピュータ、ネットワーク
上で頻繁に扱われる様になっている。
【0003】これらマルチメディアデータはディジタル
データである為、詳細で良質なデータの複製を容易に作
成することが可能である。この様なデータの著作権を保
護する為に、元のデータ中に著作権情報や利用者情報を
電子透かし情報として埋め込む技術が利用されることが
ある。
【0004】例えば、著作権情報や利用者情報が埋め込
まれたマルチメディアデータを配布した場合には、これ
ら各マルチメディアデータの使用状況と、このデータか
ら抽出される著作権情報や利用者情報とを照らし合わせ
ることにより、正しく使用されているか検証することも
可能である。
【0005】この様な電子透かし技術においては、画像
の品質と攻撃に対する耐性と埋め込み可能な情報量の3
つのパラメータが密接に関わっている。
【0006】画像の品質とは、オリジナルの画像に対す
る、電子透かし情報が埋め込まれた画像の視覚的な劣化
の度合いに相当する。電子透かし技術では、オリジナル
の画像データを微量だけ操作する為、電子透かし情報の
埋め込み操作をすると、信号レベルでは必ずオリジナル
の画像に対して劣化が生じる。しかし、この劣化を人間
の視覚が知覚できない様にすれば、オリジナルの画像の
持っている価値は失われずに使用可能となる。電子透か
し情報の埋め込みによって、視覚的な画質劣化は生じな
いことが好ましい。
【0007】攻撃に対する耐性とは、電子透かし情報が
埋め込まれた画像が種々の攻撃を受けた場合でも、その
後電子透かし情報が誤りなく抽出可能である度合いに相
当する。なお上記攻撃には、拡大、縮小、位置ずらし等
の幾何学的歪みを与える処理、或いはJPEG圧縮の様な非
可逆圧縮を行う処理、或いは部分切り取り、ノイズ付
加、濃度変換等が存在する。また、電子透かし情報が埋
め込まれた画像をプリンタで紙等の記録媒体に出力し、
その後スキャナで入力した後に電子透かし情報を抽出す
る一連の処理についても、画質が微量でも変化する可能
性が有れば、これら各処理或いは全部を攻撃であると考
えることができる。これらの攻撃を受けても電子透かし
情報は正しく抽出されることが好ましい。
【0008】埋め込み可能な情報量とは、電子透かし情
報として埋め込むことができる情報量である。これはで
きるだけ大きなほうが好ましい。
【0009】即ち、電子透かし情報を埋め込むことによ
る視覚的な画質劣化が無く、種々の攻撃に対しては電子
透かし情報の耐性が有り、埋め込み可能な情報量が大き
くなる様な、電子透かし情報の埋め込みアルゴリズムが
要求される。
【0010】上記3つのパラメータ間にはトレードオフ
の関係があると言われている。例えば、攻撃に対する耐
性を上げようとすると、画質劣化が大きくなったり、埋
め込み可能な情報量が少なくなったりする。一方、画質
劣化を小さくしようとすると、攻撃に対する耐性が小さ
くなったり、埋め込み可能な情報量も少なくなったりす
る。一方、埋め込み可能な情報量を多くしようとする
と、画質劣化が大きくなったり攻撃に対する耐性が小さ
くなったりする。
【0011】よって、電子透かし技術の開発は、電子透
かし情報の埋め込みが行われる状況に応じて、この状況
に最適な方式を検討する必要がある。
【0012】
【発明が解決しようとする課題】従来、電子透かし情報
が埋め込まれた画像データに対して、幾何的な歪みを与
える攻撃が施された場合、その後の画像データを電子透
かし情報を正確に抽出可能とする為の処理方法は未だ確
立されていなかった。
【0013】本発明は上記従来例に鑑みて成されたもの
であり、電子透かし情報が埋め込まれた画像に幾何的な
歪みが与えられた場合においても、この歪みを的確に判
断できる様にし、結果的に後段の電子透かし情報の抽出
を正確に行える様にすることを目的とする。
【0014】
【課題を解決するための手段】上記課題を解決するため
に、本発明の画像処理装置によれば、画像データの幾何
学的な状態を示すレジストレーション信号、及び前記画
像データに電子透かし情報が予め人間に目視しにくく埋
め込まれているはずの画像データから、該電子透かし情
報を抽出する為の画像処理装置であって、前記画像デー
タを周波数変換する変換手段と、前記変換手段により得
られる信号に微分処理を施す微分手段と、前記微分手段
で得られた信号と、前記レジストレーション信号に関す
る情報とに基づいて、前記画像データの幾何学的な補正
を行う補正手段とを有することを特徴とする。
【0015】
【発明の実施の形態】[1 電子透かし埋め込み装置]
以下、本実施の形態における電子透かし埋め込み装置の
概要について図を用いて説明する。
【0016】図1に、本実施の形態における電子透かし
埋め込み装置を示す。図1に示す様に、電子透かし埋め
込み装置は、色成分抽出手段0101、レジストレーシ
ョン信号埋め込み手段0102、埋め込み位置決定手段
0103、付加情報埋め込み手段0104、色成分合成
手段0105から構成される。
【0017】電子透かし埋め込み装置には、画像データ
Iが入力される。これは1画素当たり所定の複数ビット
が割り当てられた多値画像データである。また本実施の
形態では、入力される画像データIがグレイスケール画
像データであってもカラー画像データであっても対応可
能である。グレイスケール画像データは、1画素当たり
1種類の要素から構成されているものであり、カラー画
像データは、1画素当たり3種類の要素から構成されて
いるものとする。この3種類の要素とは本実施の形態で
は赤色成分、青色成分、緑色成分である。しかしなが
ら、別の色成分の組み合わせにも本発明は適用可能であ
る。
【0018】電子透かし埋め込み装置に入力された画像
データIは、まず色成分抽出手段0101に入力され
る。
【0019】入力された画像データIがカラー画像デー
タである場合には、色成分抽出手段0101においてこ
のカラー画像データから青色成分だけを分離し、後段の
レジストレーション信号埋め込み手段0102に出力す
る。
【0020】一方その他の色成分は後段の色成分合成手
段0105に出力する。即ちここでは、電子透かし情報
が埋め込まれるべき色成分のみが分離されて電子透かし
の処理系に送出される。
【0021】本実施の形態においては青色成分に対して
電子透かし情報の埋め込みが行われる。これは赤色成
分、青色成分、緑色成分の中で、人間の視覚にとっては
青色成分が最も鈍感であることによる。よって青色成分
に電子透かし情報を埋め込むことは、他の色成分に対し
て電子透かし情報を埋め込むことに比べて、電子透かし
情報による画質劣化が人間の目に知覚しにくくなる効果
がある。
【0022】また、入力された画像データIがグレイス
ケール画像データである場合には、色成分抽出手段01
01はグレイスケール画像データを一旦疑似的なカラー
画像データに変換する。ここで疑似的なカラー画像デー
タとは、1画素当たり3種類の要素から構成されるカラ
ー画像データであるが、この場合、3種類の要素の値が
全て等しい画像データである。グレイスケール画像デー
タは上記疑似的なカラー画像データに変換され、このカ
ラー画像データの中で青色成分を分離し、レジストレー
ション信号埋め込み手段0102へ出力される。
【0023】一方その他の色成分は後段の色成分合成手
段0105に出力する。この様にして、上述したカラー
画像データの場合と同様に、青色成分に対して電子透か
し情報が埋め込まれる。
【0024】なお、以後の説明では、画像データIがカ
ラー画像データである場合とグレースケール画像データ
である場合とは、できるだけ区別しない様にして説明す
る。即ち、カラー画像データと擬似的なカラー画像デー
タとは区別せずに説明する。
【0025】次に、レジストレーション信号埋め込み手
段0102について説明する。ここで、レジストレーシ
ョン信号とは、電子透かし情報の抽出の前処理として幾
何的な補正を実行する為に必要とされる信号である。
【0026】レジストレーション埋め込み手段0102
には、色成分抽出手段0101で得られた青色成分の画
像データが入力される。レジストレーション信号埋め込
み手段0102では、画像データに対して一種の電子透
かしの技術を用いてレジストレーション信号が埋め込ま
れる。即ち、レジストレーション信号を埋め込まれた画
像データにおいて、人間の視覚はレジストレーション信
号を知覚することができない。このレジストレーション
信号の埋め込み方法等、詳細は後述する。
【0027】レジストレーション信号埋め込み手段01
02はレジストレーション信号が埋め込まれた画像デー
タを出力する。
【0028】次の埋め込み位置決定手段0103は、レ
ジストレーション信号埋め込み手段0102から入力さ
れる画像データにおける、付加情報Infの埋め込み位置
を決定する。
【0029】埋め込み位置決定手段0103は、入力さ
れた画像データと共に、画像中に付加情報Infを埋め込
む位置を表す制御データを付加情報埋め込み手段010
4に出力する。
【0030】付加情報埋め込み手段0104は、上記画
像データと制御データに加えて、付加情報Inf(複数の
ビット情報)を入力する。この付加情報Infは、青色成
分の画像データにおける上記決定された埋め込み位置
に、電子透かし技術を用いて埋め込まれる。この電子透
かしの技術を用いた付加情報Infの埋め込みについても
後述する。
【0031】付加情報埋め込み手段0104からは付加
情報Infが埋め込まれた画像データが出力され、色成分
合成手段0105に入力される。
【0032】色成分合成手段0105は、前段(付加情
報埋め込み手段0104)までに処理された青色成分
と、色成分抽出手段0101から直接入力される赤色成
分及び緑色成分とを用いて、通常のカラー画像データの
形態に合成する。
【0033】以上の処理により、レジストレーション信
号及び付加情報Infが電子透かし技術によって埋め込ま
れた画像データwIが出力される。
【0034】なお、本実施の形態ではこの画像データwI
に種々の幾何的な歪みを生じる攻撃が加えられるものと
して説明する。例えば、ユーザーによる意図的な画像編
集であり、或いは画像データwIを印刷した後、印刷物を
スキャナ等でスキャンニングする動作である。攻撃を受
けた画像データは後述する図2のwI'となる。
【0035】以上説明した各手段による全体の流れを、
図31に示すフローチャートを用いて説明する。
【0036】まずステップ3402において画像データ
Iが色成分抽出手段0101に入力される。これは写真
や印刷物をスキャナなどで読み取り、画像データを生成
する工程も含む。更に青色成分を分離し、後段のレジス
トレーション信号入力に使用する。
【0037】次に、ステップ3403においてレジスト
レーション信号を発生し、このレジストレーション信号
がステップ3404で埋め込まれる。このステップ34
04におけるレジストレーション信号埋め込み処理は、
図1におけるレジストレーション信号埋め込み手段01
02の内部で実行される処理に相当し、詳細な説明は後
述する。
【0038】更に、ステップ3405によりマスクが作
成され、作成されたマスクは、ステップ3406に入力
され、埋め込みビット情報と埋め込み位置の関係を規定
する。ステップ3407ではマスクを拡大マスクに拡張
する。このマスク・パターン配列対応手段についても詳
細な説明は後述する。
【0039】ステップ3403,3404においてレジ
ストレーション信号が埋め込まれた画像データに対し
て、ステップ3408において付加情報Infが埋め込ま
れる。この付加情報埋め込み処理はマクロブロック単位
で付加情報Infを画像全体に繰り返し埋め込む。この処
理は後述の図10で詳細に述べる。ここで、マクロブロ
ックとは、最小埋め込み単位のことを示し、このマクロ
ブロックに相当する画像領域に1つの完結した付加情報
Infの全情報が埋め込まれる。
【0040】上記画像データに対して付加情報Infが埋
め込まれた後に、ステップ3409により、電子透かし
情報が埋め込まれた画像データwIが出力される。
【0041】上述した様に、この画像データwIは、後述
する図32の電子透かし抽出開始処理が実行される前
に、種々の幾何的な歪みを生じる攻撃が加えられるかも
しれない。
【0042】[2 電子透かし抽出装置]次に、本実施
の形態における電子透かし抽出装置の概要について説明
する。
【0043】図2に、本実施の形態における電子透かし
抽出装置を示す。図2に示す様に、電子透かし抽出装置
は、色成分抽出手段0201、レジストレーション手段
0202、及び付加情報抽出手段0203から構成され
る。
【0044】電子透かし抽出装置には、画像データwI'
が入力される。画像データwI'は画像データwIに対して
種々の幾何的な歪みを生じる攻撃を受けている可能性を
有する画像データである。攻撃には、JPEG圧縮などの非
可逆圧縮、スケーリング、回転、印刷&スキャニングな
ど種々のものを含み、更にこれらの組合せも考えられ
る。
【0045】従って、画像データwI'とwIとは内容が同
一であることが理想的ではあるが、実際にはこの2つの
画像データは内容が著しく異なっていることが多い。
【0046】色成分抽出手段0201は、画像データw
I'を入力し、青色成分を抽出した後、この青色成分の画
像データを後段のレジストレーション手段0202へ出
力する。画像データwI'のうち、青色成分以外の赤色成
分、及び緑色成分は必要無いのでここで廃棄される。
【0047】レジストレーション手段0202には、色
成分抽出手段0201で得られた青色成分の画像データ
wI1'が入力される。そしてこの青色成分の画像データwI
1'を用いて、幾何学的な歪みが補正された画像データwI
2'を生成し、付加情報埋め込み開始点pを判別する。
【0048】上述した様に、画像データwI'は画像デー
タwIとは異なるスケールである可能性を有するのに対
し、画像データwI2'は画像データwIと必ず同一のスケー
ルになる。この理由、及び画像データwI2'を画像データ
wIと同じスケールにする処理の詳細は後述する。
【0049】ここで付加情報埋め込み開始点pとは、レ
ジストレーション手段0202の後段の付加情報抽出手
段0203が実行する抽出処理の開始点を示す情報であ
る。これに関しても詳細は後述する。
【0050】レジストレーション手段0202は、画像
データwI2'及び付加情報埋め込み開始点pの情報を、付
加情報抽出手段0203に出力する。
【0051】付加情報抽出手段0203は、付加情報埋
め込み手段0103での埋め込み方式に合わせた所定の
処理を施こすことによって、画像データwI2'に埋め込ま
れている電子透かし情報を抽出することが可能であり、
抽出された付加情報Infを出力する。
【0052】以上説明した各手段による全体の流れを、
図32のフローチャートを用いて説明する。まずステッ
プ3502では画像データwI'が入力される。この画像
データwI'は、画像データwIであると予想される画像デ
ータをネットワークやメモリから読み出したり、画像デ
ータwIに基づいて印刷した物をスキャナ等でスキャンニ
ングしたりすることで得られる。後者の場合には画像デ
ータwI'が画像データwIとは著しく異なっている可能性
が非常に高い。
【0053】また、この画像データwI'の青色成分のみ
を抽出され、次のステップに使用される。
【0054】次に、ステップ3503では、入力された
青色成分の画像データwI1'のスケールが補正される。
【0055】更にステップ3504では、入力された青
色成分の画像データのオフセットwI1'が補正される。こ
のスケール合わせ処理、及びオフセット合わせ処理は、
図2のレジストレーション手段0202の内部で実行さ
れる処理であり、詳細な説明は後述する。
【0056】次に第1のパターン配列を用いた抽出処理
がステップ3506において、第2のパターン配列を用
いた抽出処理がステップ3505において、スケールと
オフセットが既に補正された画像データwI2'から、夫々
埋め込まれている付加情報Infが抽出される。
【0057】統計検定ステップ3507では、上記抽出
された電子透かし情報の確からしさを算出、判定され、
正しい電子透かし情報ではないと判定されればステップ
3502に戻り、電子透かし情報が埋め込まれていると
考えられる画像を入力しなおす。一方、充分正確な電子
透かし情報であると判定された場合には、ステップ35
08の比較処理により電子透かし情報(付加情報Inf)
の抽出を行う。また、ステップ3510では上記確から
しさを示す情報を後述する信頼性指標Dとして表示す
る。
【0058】以上の第1のパターン配列を用いた抽出処
理、及び第2のパターン配列を用いた抽出処理、及び統
計検定処理、及び比較処理は、図2の付加情報抽出手段
0203の内部で実行される処理であり、詳細な説明は
後述する。
【0059】[3 各部の詳細説明]次に、各部の詳細
な説明をする。
【0060】まず、電子透かし抽出側におけるレジスト
レーション手段0202、ステップ3503、ステップ
3504が実行するレジストレーションという処理につ
いて説明する。
【0061】レジストレーション処理とは電子透かし情
報の抽出の際に、電子透かし抽出装置に入力された画像
データwI'から電子透かし情報が抽出できる様にする前
処理である。以下では、まず印刷系の処理が施された画
像データがどの様な変化を受けるかを考える。そして、
その様な変化に対するレジストレーション処理を考察
し、印刷系に対するレジストレーション処理を考える。
【0062】電子透かし抽出装置には、電子透かし埋め
込み装置から出力された画像データwIがそのまま入力さ
れるとは限らない。
【0063】例えば、画像データwIがYMCKインクジ
ェットプリンタによって印刷され、この印刷物をスキャ
ナでスキャンニングする場合を挙げて考察する。
【0064】この時、プリンタによる出力解像度と、ス
キャナによる入力解像度が異なる場合には、元のカラー
画像データwIとスキャンニングで得られた画像データの
スケールは異なってしまう。よって、得られた画像デー
タwI'からはそのまま正確に電子透かし情報を抽出でき
る可能性は少ない。従って、これらのスケールの違いを
補正できる手段を備えておくことが必要である。
【0065】本実施の形態では、入力解像度と出力解像
度の両方が既知である為、これらの比からスケールの比
を算出できる。例えば、出力解像度が600dpiであり
入力解像度が300dpiである場合には、印刷前の画像
とスキャンニング後の画像のスケールの比は2倍であ
る。よって、算出されたスケールの比にしたがって、適
当なスケーリングアルゴリズムを用いて画像データwI'
に対してスケーリングを施す。これにより画像データwI
と画像データwI'が表す画像サイズを同じスケールにで
きる。
【0066】しかし、全ての場合に出力と入力の解像度
が既知であるとは限らない。双方の解像度が既知でない
場合には、上述の方法は使用できない。この場合には、
スケールの違いを補正する手段に加えて、更にスケール
の比を知る手段が必要である。
【0067】また、画像データwIに印刷系の処理が施さ
れた画像は、スキャナのスキャンニングにより入力され
た後には図3に示す様な画像となる。図3においては0
301全体が画像データwI'が表す画像である。この画
像データ0301は、画像データwIにより表される本来
の画像0302と白色の余白部分0303から構成され
る。この様な余白は、ユーザーがマウスなどによって切
り取った場合には不正確になってしまう。
【0068】以上挙げた点が、印刷系を介して得られた
画像データwI'を表す画像には、必ず生じていると考え
られ、画像データwIに印刷系の処理が施される場合であ
れば、これらを解決する必要がある。
【0069】以上、電子透かしを抽出する前に印刷系の
処理を少なくとも1度介してから画像データを得る場合
について説明したが、この様な状況は人為的な編集によ
っても起こり得る。
【0070】レジストレーション処理とは上述した様な
問題を解決する為の処理であり、付加情報抽出処理の前
処理として定義される。
【0071】以下、入出力解像度の比が不明である場合
を想定し、上記問題点を解決する為に備えられた、レジ
ストレーション信号埋め込み手段及びレジストレーショ
ン手段について説明する。
【0072】[3−1 レジストレーション信号埋め込
み処理]まず、レジストレーション信号埋め込み手段0
102(ステップ3404)について詳細を説明する。
【0073】レジストレーション信号埋め込み手段01
02は、付加情報埋め込み手段0104の前段に位置す
るものである。この手段0102は、図2のレジストレ
ーション手段における画像データwI'のレジストレーシ
ョンに参照されるレジストレーション信号を、元の画像
データに予め埋め込むものである。このレジストレーシ
ョン信号は、電子透かし情報として人間の目に見えにく
く画像データ(本実施の形態ではカラー画像データの青
色成分)に埋め込まれる。
【0074】図4にレジストレーション信号埋め込み手
段0102の内部構成を示す。レジストレーション信号
埋め込み手段0102は、図4のブロック分割手段04
01、フーリエ変換手段0402、加算手段0403、
逆フーリエ変換手段0404、ブロック合成手段040
5から構成される。以下、各手段の詳細について解説す
る。
【0075】ブロック分割手段0401は、入力された
画像データを互いに重ならない複数のブロックに分割す
る。このブロックのサイズは本実施の形態においては2
のべき乗に定める。実際はこれ以外のサイズでも適応可
能であるが、ブロックのサイズが2のべき乗である場合
は、ブロック分割手段0401の後に結合されるフーリ
エ変換手段0402において高速の処理を行うことが可
能である。
【0076】ブロック分割手段0401により分割され
たブロックは、二つの集合I1とI2に分けられ、このうち
I1は後段のフーリエ変換手段0402に入力され、I2
後段のブロック合成手段0405に入力される。本実施
の形態では、I1としては、ブロック分割手段0401に
より得られた各ブロックのうち、画像データIの中で最
も中心付近に位置する一つのブロックが選択され、残り
のブロックは全てI2として選択される。
【0077】これは、本実施の形態が少なくとも一つの
ブロックを用いることで実現可能であり、ブロックの数
が少ない方が処理時間を短縮できることによる。しかし
ながら本発明はこれに限らず、I1として二つ以上のブロ
ックを選択する場合も範疇に含む。
【0078】また、どの様なサイズでブロックに分割さ
れ、どのブロックがレジストレーション信号の埋め込み
対象に選択されるかという情報は、電子透かし埋め込み
装置と、電子透かし抽出装置で共有する必要がある。
【0079】ブロック分割手段0401の分割で得られ
た画像データの一部I1は、フーリエ変換手段0402に
入力される。
【0080】次にフーリエ変換手段0402は、入力さ
れた画像データI1に対してフーリエ変換を施す。入力さ
れた画像データI1の元のデータ形態を空間領域と呼ぶの
に対して、フーリエ変換された後のデータ形態を周波数
領域と呼ぶ。フーリエ変換は、入力された全てのブロッ
クに対して施される。なお本実施の形態では、入力され
るブロックのサイズが2のべき乗であるので、処理の高
速化の為に高速フーリエ変換を用いる。
【0081】高速フーリエ変換とは、フーリエ変換がn
×n回の演算量を必要とするのに対して、(n/2)log
2(n)の演算量で実行可能な変換アルゴリズムである。た
だしここでnは正の整数である。高速フーリエ変換とフ
ーリエ変換は、演算結果を得る為の速度が違うだけであ
り、両者からは等しい結果が得られる。よって本実施の
形態では、高速フーリエ変換とフーリエ変換を区別して
説明しない。
【0082】尚、フーリエ変換により得られた周波数領
域の画像データは振幅スペクトルと位相スペクトルによ
って表現される。このうち振幅スペクトルだけが加算手
段0403に入力される。一方で、位相スペクトルは逆
フーリエ変換装置0404に入力される。
【0083】次に、加算手段0403について説明す
る。加算手段0403には、上記振幅スペクトルと共
に、レジストレーション信号と呼ばれる信号rが別に入
力される。レジストレーション信号の例としては、図5
に示す様なインパルス信号が挙げられる。
【0084】図5では、フーリエ変換されて得られた2
次元空間周波数成分のうちの振幅スペクトルが示されて
いる。中心が低周波成分、周囲が高周波数成分である。
0501は本来の画像成分が持っている信号成分の振幅
スペクトルであり、写真等の自然画像に相当する信号で
は低域に多くの大きな信号が集中する。一方で高域には
ほとんど信号が存在しない。
【0085】なお、本実施の形態では自然画像に一連の
処理を施すことを想定して説明するが、本発明はこれに
限らず、文書画像、CG画像等も同様に処理しても良
い。ただし、本実施の形態は中間濃度が比較的多い自然
画像を処理する場合に特に有効である。
【0086】図5は、自然画像の本来持つ信号0501
に対して周波数領域の信号の水平垂直ナイキスト周波数
成分にインパルス信号0502、0503、0504、
0505を加えた本実施の形態の一例である。この例の
様に、レジストレーション信号はインパルス性の信号で
あることが望ましい。これは後述する電子透かし抽出装
置において、レジストレーション信号だけを抽出しやす
いからである。
【0087】図5では入力信号のナイキスト周波数成分
に対してインパルス信号を加えているが、本発明はこれ
に限定されることはない。即ち、電子透かし情報が埋め
込まれた画像が攻撃を受けた場合にもレジストレーショ
ン信号が除去されないものであれば良い。上述した様
に、JPEG圧縮などの非可逆圧縮方式はローパスフィルタ
的な効果がある。よって、ここでの情報圧縮の対象とな
る高周波成分へインパルス信号を埋めこんでも、圧縮・
伸張処理により除去されてしまう可能性がある。
【0088】一方で、低周波成分へのインパルスの埋め
込みは、高周波成分への埋め込みと比べて、人間の視覚
特性からノイズとして知覚されやすいという欠点があ
る。よって、本実施の形態では人間の視覚に認識しにく
い第1の周波数以上であり、かつ非可逆圧縮・伸張処理
により容易に除去されない第2の周波数以下である中間
レベルの周波数へインパルス信号を埋め込むものとす
る。また、このレジストレーション信号は、加算手段0
403に入力された各ブロック(本実施の形態では1ブ
ロック)に対して加えられる。
【0089】加算手段0403は、周波数領域の画像デ
ータの振幅スペクトルにレジストレーション信号が加え
られた信号を逆フーリエ変換手段0404に出力する。
【0090】逆フーリエ変換手段0404は、入力され
た周波数領域の画像データに対して逆フーリエ変換を施
す。この逆フーリエ変換は、入力された全てのブロック
に対して施される。上述のフーリエ変換手段0402の
場合と同様に、入力されたブロックのサイズが2のべき
乗であるので処理の高速化の為に高速フーリエ変換を用
いる。逆フーリエ変換手段0404へ入力された周波数
領域の信号は、逆フーリエ変換されることにより空間領
域の信号に変換され出力される。
【0091】逆フーリエ変換手段0404から出力され
た空間領域の画像データは、ブロック結合手段0405
に入力される。
【0092】ブロック結合手段0405は、ブロック分
割手段0405で行われた分割と逆の処理を行う。ブロ
ック結合手段0405の処理の結果、画像データ(青色
成分)が再構成されて出力される。
【0093】以上、図1に示したレジストレーション信
号埋め込み手段0102の詳細を述べた。
【0094】図4では、レジストレーション信号をフー
リ変換領域において埋め込む方式を説明した。一方で、
レジストレーション信号を空間領域において埋め込む方
式も考えられる。この方式を図29を用いて説明する。
【0095】図29は、ブロック分割手段3201、加
算手段3202、ブロック合成手段3203、逆フーリ
エ変換手段3204から構成される。
【0096】ブロック分割手段3201及びブロック合
成手段3203は、図4におけるブロック分割手段04
01及びブロック合成手段0405と同様の動作をす
る。レジストレーション信号埋め込み手段0102に入
力される画像データはまずブロック分割手段3201に
入力され、分割される。ここで得られたブロックは加算
手段3202に入力される。一方で、レジストレーショ
ン信号rは逆フーリエ変換手段3204に入力され、逆
フーリエ変換処理により信号r'に変換される。ここで、
レジストレーション信号rは、図5に示したものと同様
に周波数領域上の信号である。加算手段3202には、
ブロック分割手段3201からのブロックと逆フーリエ
変換手段3204からの信号r'が入力され、夫々加算さ
れる。加算手段3202から出力された信号はブロック
合成手段3203に入力され、画像データ(青色成分)
が再構成されて出力される。
【0097】以上図29の手段構成は、図4の手段構成
と同様の処理を空間領域で行うものである。図4の手段
構成に比べて、フーリエ変換手段を必要としないので高
速な処理を行うことが可能である。
【0098】更に、図29の中で、信号r'は入力画像デ
ータIに対して独立な信号である。よって信号r'の算
出、即ち逆フーリエ変換手段3204の処理は入力画像
データIが入力される毎に実行される必要はなく、予め
r'を生成しておくことが可能である。この場合には、図
29の手段構成から逆フーリエ変換手段3204を除去
することもでき、更に高速にレジストレーション信号を
埋め込むことができる。このレジストレーション信号を
参照するレジストレーション処理については後述する。
【0099】≪ パッチワーク法 ≫本実施の形態では
付加情報Infの埋め込みの為にパッチワーク法と呼ばれ
る原理を用いている。そこで、まずパッチワーク法の原
理を説明する。
【0100】パッチワーク法では画像に対して統計的偏
りを生じさせることによって付加情報Infの埋め込みを
実現している。
【0101】これを図30を用いて説明する。図30に
おいては、3301、3302は各々画素の部分集合、
3303は画像全体である。画像全体3303から二つ
の部分集合A 3301とB 3302を選択する。
【0102】この二つの部分集合の選択方法は、互いに
重ならならなければ本実施の形態におけるパッチワーク
法による付加情報Infの埋め込みが実行可能である。た
だし、この二つの部分集合の大きさや選択方法は、この
パッチワーク法によって埋め込まれた付加情報Infの耐
性、即ち画像データwIが攻撃を受けた際に付加情報Inf
を失わない為の強度に大きく影響を及ぼす。これについ
ては後述する。
【0103】今、選択した部分集合Aの要素の持つ値を
{an}、部分集合Bの要素の持つ値を{bn}とする。{an}、
{bn}は具体的には、各部分集合に含まれる各画素の値
(本実施の形態ではカラー画像データ中の青色成分の値
に相当)である。
【0104】ここで、次の指標dを定義する。 d = 1/N Σ (ai − bi) これは、二つの集合の画素値の差の期待値を示してい
る。
【0105】一般的な自然画像に対して、適当な部分集
合Aと部分集合Bを選択し、指標dを定義すると、 d ≒ 0 となる性質がある。以降ではdを信頼度距離と呼ぶ。
【0106】一方で、付加情報Infを構成する各ビット
の埋め込み操作として、 a'i = ai + c b'i = bi − c という操作を行う。これは部分集合Aの要素全てに対し
て値cを加え、部分集合Bの要素全てに対してcを減ずる
という操作である。
【0107】ここで、先程の場合と同様に、付加情報In
fが埋め込まれた画像から部分集合Aと部分集合Bを選択
し、指標dを計算する。
【0108】すると、 d = 1/N Σ (a'i - b'i) = 1/N Σ {(ai + c) - (bi - c)} = 1/N Σ (ai - bi) + 2c = 2c となり0にはならない。
【0109】即ち、ある画像が与えられた時に、画像に
対して信頼度距離dを算出することによって、d≒0なら
ば付加情報Infは埋め込まれておらず、一方でdが0から
一定量以上離れた値であるなら付加情報Infが埋め込ま
れていると判断できる。
【0110】以上がパッチワーク法の基本的な考え方で
ある。
【0111】このパッチワーク法の原理を応用し、本実
施の形態では複数のビットの情報を埋め込む。本方法に
おいては、部分集合Aと部分集合Bの選択の方法について
もパターン配列によって定義している。
【0112】上述の方法においては、原画像の所定の要
素に対してパターン配列の要素を加えたり減じたりする
ことによって、付加情報Infの埋め込みを実現してい
る。
【0113】簡単なパターン配列の例を図9に示す。図
9は、1ビットを埋め込む為に8×8画素を参照する場
合の、元の画像からの画素値の変更量を示すパターン配
列である。図9の様に、パターン配列は正の値を持つ配
列要素、負の値を持つ配列要素、及び0の値を持つ配列
要素から構成される。
【0114】図9のパターンにおいて、+cの配列要素
で示される位置は対応位置の画素値をcだけ上昇させる
位置を示し、上述した部分集合Aに相当する位置であ
る。一方−cの配列要素で示される位置は対応位置の画
素値をc減少させる位置を示し、上述した部分集合Bに
相当する位置である。また0で示される位置は上述した
部分集合A,B以外の位置であることを示す。
【0115】本実施の形態では、画像の全体的な濃度を
変化させない為にも正の値を持つ配列要素の個数と負の
値を持つ配列要素の個数を等しくしている。即ち、1つ
のパターン配列において全ての配列要素の和が0になっ
ている。なお、後述する付加情報Infの抽出操作の時に
はこの条件が必須である。
【0116】以上の様なパターン配列を用いて付加情報
Infを構成する各ビット情報の埋め込み操作を行う。
【0117】本実施の形態では、図9のパターンを元の
画像データにおける互いに異なる領域に複数回配置して
画素値を上昇/減少させることにより、複数のビット情
報、即ち付加情報Infを埋め込む。言い換えれば、1つ
の画像の互いに異なる領域に、部分集合AとBの組み合わ
せだけでなく、部分集合A’とB’、部分集合A”とB”、
…という複数の組み合わせを想定することで、複数のビ
ットからなる付加情報Infを埋め込む。
【0118】なお、本実施の形態では元の画像データが
大きい場合には、繰り返し付加情報Infを埋め込むこと
になる。これはパッチワーク法が統計的性質を利用して
いるものであることから、統計的性質が現れるのに充分
な数を必要とすることによる。
【0119】また、本実施の形態では複数ビットを埋め
込む際にパターン配列を用いて画素値を変更する領域が
重ならない様にする為、予め互いのビット同志でパター
ン配列を使用する相対位置を決定する。即ち、付加情報
Infを構成する1ビット目の情報を埋め込む為のパター
ン配列の位置と、2ビット目の情報を埋め込む為のパタ
ーン配列の位置の関係は適切に定められる。
【0120】例えば付加情報Infが16ビットで構成さ
れていれば、1〜16ビット目の夫々の8×8画素のパ
ターン配列の位置関係は、32×32画素よりも大きい
サイズの領域上で画質劣化が少なくなる様に相対的に与
えられる。
【0121】更に、上記付加情報Inf(これを構成する
各ビット情報)は、画像データが大きい場合には、でき
るだけ多くの回数繰り返して埋め込む。これは、付加情
報Infの各ビットを正しく抽出可能とすることが目的で
ある。特に本実施の形態では、同一の付加情報Infが繰
り返し埋め込まれていることを利用した統計的な計測を
行うので、上記繰り返しは重要である。
【0122】以上説明した様な埋め込み位置の選択は、
図1における埋め込み位置決定手段0103において実
行される。次に、この埋め込み位置決定手段の動作につ
いて説明する。
【0123】[3−2埋め込み位置決定処理]図11は
埋め込み位置決定手段0113の内部構成を記してい
る。
【0124】図11のマスク作成手段1101は、付加
情報Infを構成する各ビット情報の埋め込み位置を規
定する為のマスクの作成を行う。マスクとは、各ビット
情報に対応するパターン配列(図9参照)の相対的な配
置方法を規定する位置情報を備えたマトリクスである。
【0125】図17の1701にマスクの一例を示す。
マスクの内部には夫々係数値が割り当てられており、各
々の係数値はマスク内で等しい出現頻度を有している。
このマスクを用いる場合であれば、最大16ビットから
なる付加情報Infを埋め込むことが可能である。
【0126】次にマスク参照手段1102は、マスク作
成手段1101で作成したマスクを読みこみ、マスク内
の各係数値と、各ビット情報が何ビット目かの情報を対
応付けて、各ビット情報を埋め込む為のパターン配列の
配置方法を決定する。
【0127】更にマスク・パターン配列対応手段110
3は、マスク内の各係数値の位置に各パターン配列の配
列要素(8×8サイズ)を展開する。即ち、図17の1
701に示されるマスクの各係数値(1マス)を、同図
1702の様に8×8倍にして各パターン配列の埋め込
み位置として参照可能とする。
【0128】後述の付加情報埋め込み手段0104は、
図17の埋め込み先頭座標1702を参照して、パター
ンサイズを用いて各ビット情報を埋め込むことになる。
【0129】なお、本実施の形態ではマスク作成手段1
101に、画像データ(青色成分)を入力する毎に上記
マスクを作成する。よって、大きいサイズの画像データ
を入力する場合には、複数回繰り返して同一の付加情報
Infを埋め込むことになる。
【0130】上記の方法では画像から付加情報Infを抽
出する場合に、上記マスクの構成(係数値の配列)が鍵
の役割を果たす。即ち、鍵の所有者だけが情報の抽出を
行えるという効果がある。
【0131】なお本発明は、リアルタイムにマスクを作
成せずに、予め作成しておいたマスクをマスク作成手段
1101の内部記憶手段などに記憶させておき、必要時
に呼び出す場合も範疇に含む。この場合には高速に後段
の処理に移行することができる。
【0132】次に、埋め込み位置決定手段0103の中
で実行される各処理の詳細を説明する。
【0133】[3−2−1 マスク作成手段]はじめに、
マスク作成手段1101について説明する。
【0134】パッチワーク法を用いた付加情報Infの埋
め込みにおいて、攻撃耐性を強める為に画素値に大きな
操作を加えて情報を埋め込んだ場合(例えばパターン配
列のCの値を大きく設定した場合)には、元の画像デー
タが表す画像において画素値の急激な変化の有るいわゆ
るエッジ部分では画質の劣化は比較的目立ちにくいが、
画素値の変化の少ない平坦部では画素値を操作した部分
がノイズとして目立ってしまう。
【0135】図13に人間の目で知覚する空間周波数特
性を示す。横軸は空間周波数(Radial Frequency)を示
し、縦軸は視覚の応答値である。画素値を操作し、情報
を埋めこんだ場合、人間の目が敏感に知覚できる低周波
数領域では画質劣化が目立つことが図13から分かる。
【0136】その為本実施の形態では、多値画像の2値
化処理に通常使用されているブルーノイズマスクやコー
ンマスクの特性を考慮し、各ビットに対応するパターン
の配置を行う。
【0137】次にブルーノイズマスク、コーンマスクの
特性について簡単な説明を行う。
【0138】まず初めにブルーノイズマスクの特性につ
いて説明する。
【0139】ブルーノイズマスクはどの様な閾値で2値
化してもブルーノイズパターンとなる特性を有する。こ
のブルーノイズパターンとは空間周波数が高周波領域に
偏った周波数特性を示すパターンである。
【0140】図37に、あるブルーノイズマスクの一部
を示す。
【0141】また図14の1401には、閾値10で2
値化したブルーノイズマスクの空間周波数特性の概略図
を示す。
【0142】1401の横軸はRadial Frequencyであ
り、ブルーノイズマスクをフーリエ変換したときの原点
(直流成分)からの距離を示している。縦軸はPower sp
ectrumであり、横軸Radial Frequencyの示す距離にある
振幅成分の2乗和をとり平均化した値である。なお、同
図は画像の2次元周波数特性を1次元グラフ化し視覚的
に分かりやすくしたものである。
【0143】図13と比較してみると、ブルーノイズマ
スクは高周波成分に偏りがある為、人間の目に知覚され
にくいことが分かる。従って、インクジェットプリンタ
等では、ドットを用いた面積階調で多値画像の階調を表
現する際に、ブルーノイズマスクを用いることで、空間
周波数成分を高周波に偏らせ、人の目に目立つことなく
面積階調を表現できることが知られている。
【0144】次にブルーノイズマスクの生成過程の一例
を以下に示す。 1. ホワイトノイズを生成する 2. 階調gの2値画像Pgl(初期値はホワイトノイズマス
ク)にローパスフィルタリングを行い、多値画像P'gl
生成 3. 階調 g (初期値:127)の画像とローパスフィルタ
リング画像P'gl(多値)を比較し、誤差の大きい順に、
2値画像Pgの白黒の画素を反転させ、2値画像Pgl+1
得る。 4. 誤差が最小になるまで、2,3の操作を繰り返し、
2値画像Pgl(初期値はホワイトノイズマスク)を少しづ
つ階調 g(初期値:127) の2値画像Pg(ブルーノイズ
マスク)に変えていく。 5. Pg画像に階調g+1(g-1)の2値の黒(白)の点を
ランダムな位置に与え、2,3の操作を繰り返し、P
g+1(Pg-1)を得る。
【0145】以上の操作を繰り返すことにより、全ての
階調でのブルーノイズマスクを作成し、ディザマトリク
スを生成する。
【0146】例えば32×32のブルーノイズマスクで
は1階調ごとに4点増加(減少)する。
【0147】ただし、このとき256階調持たせる為に
前の階調gで決まった黒(白)のビットは反転できない
為、低または高階調では制限条件が厳しくなり、一様性
に乏しいランダムパターンしか得られないという欠点が
ある。
【0148】図12にブルーノイズマスクを構成する各
係数の出現頻度分布(ヒストグラム)1201を示す。
図12では0〜255の全ての値(係数)がマスク内に
同数存在している。
【0149】上記ブルーノイズマスクが多値画像の2値
化に用いられる技術は良く知られており、例えば、「J.
Opt.Soc.Am A/Vol.9, No.11/November 1992 Digital ha
lftoning technique using a blue-noise mask Tehopha
no Mitsa, Kevin J.Parker」等に詳しく開示されてい
る。
【0150】次にコーンマスクの特性について説明す
る。
【0151】コーンマスクは、このマスクに含まれる各
係数を2値化した場合に、ここで得られた2値情報を表
す空間周波数領域上で、図14の1402に示される様
に、周期的または擬似周期的なピークが発生することを
1つの特徴とする。ただし、低周波領域ではピークが立
たない様に設計されている。
【0152】図38に、あるコーンマスクの係数配列の
一部を示す。
【0153】コーンマスクをどの閾値で2値化した場合
にも、ドット間で適度な距離が保たれる為、低周波領域
でのピークは立たない。
【0154】図14の1402に、コーンマスクの閾値
10で2値化した場合の空間周波数特性の概略図を示
す。1401のブルーノイズマスクの空間周波数特性と
同じく、1402の特性でも低周波成分が少ないことが
分かる。
【0155】コーンマスクの場合は、閾値が低い場合で
あっても高い場合であっても、ブルーノイズマスクが持
つ低域周波数より高い周波数からピークが発生する為、
ブルーノイズマスクに比べて埋め込み位置に密集した部
分が少なくなる。その為、付加情報Infを埋め込んだ際
に生じる埋め込みノイズは、ブルーノイズに比べて更に
目立たなくなる利点がある。
【0156】またコーンマスクを構成する係数の使用頻
度も、ブルーノイズマスクの時と同じく、図12の12
01で示す出現頻度分布(ヒストグラム)になる。
【0157】従って、このマスクの係数に対応付けて、
付加情報Infを構成する各ビット情報に対応するパター
ンを、画像データに埋め込む様にするならば、この画像
データ中には各ビット情報に対応するパターンを同数配
置することができ、結果的に付加情報Infをバランス良
く埋め込むことができる。
【0158】本実施の形態では以上の利点から埋め込み
参照マスクにコーンマスクを用いることとする。
【0159】[3−2−2 マスク参照手段]マスク作成
手段1101で作成されたマスク(コーンマスク)は、
マスク参照手段1102に入力される。
【0160】マスク参照手段1102では、画像に埋め
込むNビットの情報の埋め込み位置とマスクの番号(画
素値)を対応付けて、埋め込み位置を決定する。
【0161】マスク参照手段1102で行う埋め込み位
置決定方法を説明する。
【0162】本実施の形態では上述したコーンマスクを
用いるが、ここでは説明を分かりやすくする為、図15
の1501に示す4×4マスクを用いて説明する。
【0163】図15のマスクは、4×4個の係数を有
し、0〜15までの係数値を1つずつ配置しているマス
クである。この4×4マスクを用い付加情報Infの埋め
込み位置の参照を行う。この説明で用いるマスクの場合
には最大16ビットで構成される付加情報Infを埋め込
むことができるが、以下8ビットの付加情報Infを埋め
込む場合を説明する。
【0164】まず付加情報Infの構成を、図36を用い
て説明する。同図の様に付加情報Infは、スタートビッ
トInf1と利用情報Inf2から構成される。
【0165】このスタートビットInf1は、理想位置から
実際の付加情報Infが埋め込まれている位置がずれてい
ることを認識し、これに合わせて電子透かし(付加情報
Inf)の抽出開始位置を補正する為に、電子透かし抽出
装置側に含まれるオフセット合わせ手段で用いられる。
詳細については後述する。
【0166】また、利用情報Inf2は、本来の付加情報、
即ち実際に画像データIの付加的情報として利用される
情報である。この情報には、例えば画像データwIの不正
利用の際に原因を追跡することを目的とするならば、図
1に示す装置のID或いはユーザーのID等が含まれ
る。また、画像データwIの印刷物をコピー禁止にするの
であれば、コピーが禁止であることを示す制御情報が含
まれる。
【0167】本実施の形態では、スタートビットは5ビ
ットとし、「11111」というビット列を用いる。し
かし、本発明はこれに限らず、付加情報Infのうち5ビ
ット以外のビット数をスタートビットとして用いること
も可能であり、同様に「11111」というビット列以
外を用いることも可能である。ただし、スタートビット
のビット数とビット系列は電子透かし埋め込み装置と電
子透かし抽出装置で共有しておく必要がある。
【0168】上述した様な4×4個の係数からなるコー
ンマスクを用いて、スタートビット5ビット、利用情報
3ビットの合計8ビットの付加情報Infを埋め込む簡単
な例について説明をする。
【0169】しかしながら、本発明はこれに限定される
ものではない。例えば32×32のコーンマスクを用い
て、スタートビット5bit、利用情報64bitからなる合
計69bitの付加情報Infを埋め込み場合等にも適用可能
である。
【0170】付加情報Infはスタートビットが5ビット
「11111」で利用情報が3ビット「010」とす
る。夫々1番目が1、2番目が1、3番目が1、4番目
が1、5番目が1、6番目が0、7番目が1、8番目が
0のビット情報を持つ。
【0171】これら各ビットに相当するパターン(図9
参照)の夫々をコーンマスクの係数の1つずつに対応す
る位置に割り当て、この位置関係に従って元の画像デー
タの各画素値を±c変更する。これにより、1つのコー
ンマスクに相当するサイズの元の画像データに対して、
1つの付加情報Infが埋め込まれることとなる。
【0172】本実施の形態では、付加情報Infを埋めこ
む為に最低限必要なビット数に基づいて、ある閾値を決
定し、上記コーンマスクにおける、この閾値以下の係数
が配置されている位置に、対応する各ビット情報を埋め
こむ。これにより付加情報Infのビット数に関わらず1
つのコーンマスクには1つの付加情報Infが埋め込まれ
る。
【0173】なお、本発明は上記方法に限らず、ある閾
値以上の係数が配置されている位置に、対応する各ビッ
ト情報が埋め込まれる様にし、これを前提に閾値を決定
する様にしても良い。
【0174】次に本実施の形態では、埋め込みに用いる
閾値以下の係数の個数が、マスク全体の係数の個数にお
いて占める割合を、埋め込み充填率と呼ぶことにする。
【0175】8ビットの付加情報Infを整数回正しく埋
め込むには、図15のマスク1501内においてどの係
数を埋め込み参照位置に用いるかを決定する為の閾値
は、8或いは16である必要が有る。この閾値は、耐性
及び画質への影響を考慮して最適なものが決定される。
【0176】ここで、マスクの閾値を8とした場合に
は、埋め込み充填率は50%となる。即ち、上記マスク
と照らし合わされる元の画像データのうち50%が、図
9のパターン配列を用いた処理の対象になる。
【0177】各ビット情報とマスク内の係数の対応関係
の一例を対応表1に示す。
【0178】
【表1】 ここで、S1〜S5まではオフセット合わせ処理装置で用
いる位置合わせの為のビット情報(スタートビット)で
ある。1〜5は5ビットの利用情報である。
【0179】対応表1の対応に従うと、図16の160
1で表現された係数(0〜7)の位置に相当する入力画
像データの画素の位置に各ビット情報がパターン(図9
参照)を用いて埋め込まれる。埋め込むビット情報の順
序とマスク内の係数値の対応関係は鍵情報の一部であ
り、この対応関係を知ることなしに各ビット情報の抽出
を行うことはできない。本実施の形態では説明を簡単に
する為、対応表1の様に0から閾値までの係数値に順に
S1〜S5と利用情報の3ビットを対応させることとす
る。
【0180】次に、実際に32×32のサイズのコーン
マスクを用いて埋めこむ場合の充填率についても少し述
べておく。なお、処理の手順はマスク1501を用いた
場合と同じである。
【0181】まず最初に、埋め込み時の画質の劣化等を
考慮し、付加情報Infの埋め込みを整数回正しく行う為
に必要な閾値を決定する。
【0182】更に、付加情報Infを構成する各ビットが
等しい繰り返し数で埋め込まれる様にする為、閾値以下
の係数の数を付加情報Infを構成するビット数Nで割り、
1つのマスクサイズで各ビットを何回埋め込めるかを決
定する。
【0183】例えば、0〜255までの係数値に相当す
る元の画像データに、上述した一例のスタートビット5
ビットと利用情報64ビットからなる、69ビットの付
加情報Infを埋めこむ場合には、例えば閾値を137に
設定する。
【0184】この場合、マスク内の有効な係数値の数は
138個となる。1つの付加情報Infを表すのに必要な
ビット数は69個であるから、1つのマスクサイズにお
いては各ビット情報は138/69=2回ずつ埋め込む
ことができる。
【0185】なお、コーンマスクを用いて埋め込み位置
を決定する時に、ある閾値以下の係数値を持つ全ての点
に対して埋め込みを行う理由は、空間周波数の低周波成
分にピークが立たないというコーンマスクの特性を生か
す為である。
【0186】上述した様に埋め込み位置を決定した結
果、埋め込み充填率50%、埋め込み情報量69ビット
となる場合には、付加情報Infを構成する各ビット情報
とコーンマスクを構成する各係数値との関係は対応表2
の様になる。
【0187】
【表2】 ここでS1〜S5はスタートビットであり、オフセット合
わせ処理装置で用いる位置合わせの為のビット情報であ
る。1〜64は利用情報である。
【0188】ただし本発明はこの対応関係に限らず、0
から閾値(または閾値から255)までの係数の位置の
全てに、図9のパターンを使用して各ビット情報を順に
埋め込むならば、各ビット情報と各係数値との対応関係
は別のものであっても良い。
【0189】32×32のコーンマスクの場合には、1
つのマスク内に同じ係数を持つ位置が夫々4個存在す
る。
【0190】全ての係数に上記対応表2に基づいて、元
の画像データに各ビット情報を埋め込む場合、32×3
2,64×64等の大きいサイズのコーンマスクなら
ば、付加情報Infを構成する各ビット情報はほぼ等しい
回数埋め込まれる。また、元の画像データにおいて同一
のビット情報が拡散して埋め込まれる。
【0191】パッチワーク法では、従来各ビット情報に
対応するパターン(図9参照)が互いに重なり合わない
様にランダム的に埋め込み位置を選んでいたが、本実施
の形態では、上記コーンマスクを参照することにより同
様の効果を得る事ができ、更に画質劣化も少ない。
【0192】以上の結果、マスク参照手段1102で
は、各ビット情報に対応する埋め込み位置の座標(x,y)
を得る。
【0193】配列S[bit][num]=(x,y) でその情報を表
現すると、bitは対応表1の場合、スタートビットS1〜
S5と利用情報1〜3ビットを表す。またnumはコーンマ
スク内で繰り返し表れる各係数に付けられる順序であ
る。(x,y)にはマスク内での相対座標が収められる。
【0194】以上の操作をマスク参照手段1102で行
う。
【0195】[3−2−3 マスク・パターン配列対応
手段]マスク参照手段1102で得られた各ビット情報
のコーンマスク内での埋め込み位置は、マスク・パター
ン配列対応手段1103に入力される。
【0196】マスク参照手段1102で決定した埋め込
み位置は夫々のビット情報のパターンの位置(8×8画
素分)であるので、パッチワーク法では更に図9に示さ
れる加算領域(+c)と減算領域(−c)とそれ以外
(0)とを割り当てる必要がある。この為、マスク参照
手段1102で参照したコーンマスクの全ての位置に図
9に相当する8×8サイズのパターン配列を展開する操
作をマスク・パターン配列対応手段1103で行う。
【0197】具体的にはマスク参照手段1102で得ら
れた配列S[bit][num]=(x,y)の座標に対して、x座標に
はパターン配列の横方向のサイズを掛け、y座標にはパ
ターン配列の縦方向のサイズを掛ける操作を行う。その
結果、図17のマスク内の座標1701は、マスク内の
1画素が1パターン配列に拡大された先頭座標1702
となる。
【0198】この先頭座標から図19に示されるパター
ン配列を用いると、パターン配列の大きさを持つ領域1
703に重なり合わず埋め込みを行うことが可能にな
る。
【0199】座標(x,y)は座標(x',y')に変化するが、配
列S[bit][num]のbitとnumは変化しない。
【0200】従って、配列S[bit][num]のbitに相当する
付加情報Infを(x',y')をパターン配列を埋め込む先頭位
置とし、複数のビット情報の埋め込みが可能になる。
【0201】なお、マスク・パターン配列対応手段11
03によりコーンマスクの各係数が、8×8のパターン
配列に展開(拡大)された大きなマスクを、拡大マスク
と呼ぶ。
【0202】上記拡大マスクのサイズは、(32×8)
×(32×8)サイズとなり、このサイズが付加情報In
fを少なくとも1つ埋め込む為に必要な最低限の画像単
位(マクロブロックとする)となる。
【0203】以上がマスク・パターン配列対応手段11
03で行われる操作である。
【0204】一般的に小さなマスクは大きなマスクより
作成時にドットの配置位置に自由度が少なく、コーンマ
スクの様な所望の特性を持つマスクの作成が難しい。例
えば小さなマスクを画像データの全体に繰り返し割り当
てることで付加情報Infを埋め込む場合には、小さなマ
スクの持つ空間周波数が画像データ全体に出現する。
【0205】その一方で、一つのマスクからは完結した
付加情報Infが抽出されるので、マスクのサイズを大き
く設定することによって、切り抜き耐性(部分的な画像
データwI’から付加情報Infを抽出できる可能性)は小
さくなってしまう。ゆえに、上記切り抜き耐性と画質劣
化のバランスを考慮し、マスクのサイズを決定する必要
がある。
【0206】以上が図1の埋め込み位置決定手段010
3で行われる処理である。
【0207】[3−3付加情報埋め込み処理]以上の様に
して決定された、画像データへの各ビット情報の埋め込
み位置を参照して、図1の付加情報埋め込み手段010
4は実際に付加情報Infを埋め込む。
【0208】図10に付加情報Infを繰り返して埋め込
む処理の動作の流れを示す。
【0209】図10に示す方式では、画像全体に割り当
て可能なマクロブロックを複数個割り当て、更に、これ
ら全マクロブロックに対してまず1ビット目のビット情
報を繰り返して埋め込み、続いて2ビット目、3ビット
目…という様に繰り返して埋め込む。これは未だ埋め込
み処理が行われていないビット情報が有れば、未処理の
マクロブロックの全てに1001〜1003の埋め込み
処理を施す手順により構成されている。
【0210】しかしながら、本発明はこの順序に限らず
2つのループ処理の内外関係を逆にしても良い。即ち、
未処理のマクロブロックが有れば、これに対して未だ埋
め込んでいないビット情報を全て埋め込む手順に変更し
ても良い。
【0211】具体的には付加情報Infの埋め込みは、埋
め込まれる各ビット情報が“1”の場合には図9のパタ
ーン配列を加える。また埋め込むビットが“0”の場合
には図9のパターン配列を減ずる、即ち図9の正負符号
を反転したものを加算する。
【0212】上記加減算の処理は、図10における切り
替え手段1001を、埋め込むビット情報に応じて切り
替え制御することによって実現される。即ち、埋め込ま
れるビット情報が“1”の時は加算手段1002に接続
し、上記ビット情報が“0”の時は減算手段1003に
接続する。これら1001〜1003の処理はビット情
報とパターン配列の情報を参照しながら行われる。
【0213】ここで上記ビット情報の1つが埋め込まれ
る様子を図19に示す。同図では埋め込まれるビット情
報が“1”、即ちパターン配列を加える場合の例を示
す。
【0214】図19に示す例では、I(x,y)が元の画像、
P(x,y)が8×8のパターン配列である。8×8のパター
ン配列を構成する各係数は、このパターン配列と同じサ
イズの元の画像データ(青色成分)に重ね合わせられ、
同位置の値同志が加減算される。その結果、I'(x,y)が
算出され、ビット情報が埋め込まれた青色成分の画像デ
ータとして図1の色成分合成手段0104へ出力され
る。
【0215】上述した8×8のパターン配列を用いた加
減算処理を、上記対応表2で決定された埋め込み位置
(各ビット情報を埋め込む為のパターン配列を割り当て
た位置)の全てに対して繰り返し行う。
【0216】次に、図10の内部のループ処理が行われ
る様子を図18に示す。
【0217】図18では、各ビット情報を繰り返し埋め
込む為に、画像データ全体1801(1803)全体に、
マクロブロック1802を左上から右下までラスター順
に繰り返し割り当てて埋め込み(図10の1001〜1
003)を行う。
【0218】付加情報埋め込み手段0104により以上
の操作が行われ、付加情報Infが画像全体に埋め込まれ
る。
【0219】以上の処理により、画像データ中に付加情
報Infが埋め込まれる。この付加情報Infが埋め込まれた
画像データの各画素が、十分少ないドット数で表現され
るのであれば、パターン配列の大きさも十分小さくなる
ので、これらパターン配列の1つ1つ非常に小さい点で
しか人間の目には知覚されない。よって、コーンマスク
の空間周波数特性も維持され、人間の目に見えにくい。
【0220】[3−4 レジストレーション処理]次
に、電子透かし抽出装置側に備えられる、図2における
レジストレーション手段0202について詳細を説明す
る。
【0221】レジストレーション手段0202とは、付
加情報抽出手段0203の前段に位置する手段であり、
付加情報Inf抽出処理の前処理である。レジストレーシ
ョン手段0202には、前段の色成分抽出手段0201
により抽出された青色成分の画像が入力される。
【0222】レジストレーション手段0202では、電
子透かし埋め込み装置から出力された画像データwIと、
電子透かし抽出装置に入力された画像データwI’のスケ
ールの違いが補正される。
【0223】レジストレーション手段0202の詳細を
図7に示す。図7に示す様に、レジストレーション手段
0202は、ブロック分割手段0701、フーリエ変換
手段0702、インパルス抽出手段0703、スケーリ
ング率算出手段0704、スケーリング手段0705か
ら構成される。
【0224】ブロック分割手段0701では、上述した
レジストレーション信号埋め込み手段0102(ブロッ
ク分割手段0401)と同様のブロック分割処理が行わ
れる。この処理によって、一般的にはレジストレーショ
ン信号埋め込み手段0102(ブロック分割手段040
1)と同様のブロックを抽出することは困難である。こ
れは電子透かし情報が埋め込まれた画像データwIが印刷
系の処理を施されることによって、大きさが変化し、更
に位置がずれることによるものである。
【0225】しかし、このブロックの抽出はある程度違
っていても問題はない。これは、電子透かし埋め込み装
置において、レジストレーション信号が画像データのう
ち振幅スペクトルに埋め込まれているからである。振幅
スペクトルは、画像データの空間領域における位置ずれ
には影響されないという性質がある。よって、電子透か
し埋め込み装置と電子透かし抽出装置の夫々において、
各々のブロック分割手段により分割されたブロックが、
空間領域で多少の位置ずれが生じていても問題はない。
【0226】ブロック分割手段0701はブロック分割
した画像データをフーリエ変換手段0702に出力す
る。フーリエ変換手段0702は上述したレジストレー
ション信号埋め込み手段0102の場合と同様に、空間
領域の画像データを周波数領域の画像データに変換す
る。フーリエ変換された周波数領域の画像データは振幅
スペクトルと位相スペクトルによって表現される。この
うち振幅スペクトルだけがインパルス抽出手段0703
に入力される。一方で、位相スペクトルは破棄される。
【0227】周波数領域に変換された画像データはイン
パルス抽出手段0703に入力される。インパルス抽出
手段0703では、周波数領域に変換された画像データ
からインパルス性の信号だけを抽出する。即ち、既に画
像データに埋め込まれている図5の0502,0503,
0504,0505を抽出する。
【0228】これは公知の画像処理技術を用いて行うこ
とが可能である。例えば、周波数領域に変換された画像
データを閾値処理することで実現できる。この例を図8
(a)に示す。図8(a)にはインパルス抽出手段0703に
入力された振幅スペクトル0801を閾値0802によ
って閾値処理する様子を示す。なお説明の為、図8にお
いて変換された画像データを1次元で表現した。適当な
閾値0802を選択することによってインパルス信号を
抽出することが可能である。しかしながら、低域に存在
するインパルス信号と同じ程度の大きさを持つ様な本来
の画像データも同時に抽出してしまう。
【0229】この問題を解決した本実施の形態の方式を
図8(b)に示す。周波数領域に変換された画像データ0
801に対して2次微分処理を施す。これはラプラシア
ンフィルタ等を施すことに等しい。周波数領域に変換さ
れた画像データ0801に対して2次微分を施したもの
を0803に示す。このデータ0803に対して適当な
閾値0804を選択し閾値処理を施すことによって、イ
ンパルス信号を抽出可能である。
【0230】このインパルス信号の抽出に関してもう少
し詳細な原理を図26を用いて説明する。なお、この図
では上述したレジストレーション信号埋め込み側の処理
も記載している。
【0231】レジストレーション信号埋め込み手段01
02においては、空間領域の画像データ2601が周波
数領域に変換され画像データ2602となり、周波数領
域においてインパルス信号2603が加えられる。
【0232】インパルス信号(レジストレーション信
号)2603が加えられた周波数領域の画像データは、
逆周波数変換されることによって空間領域の信号260
1’へ再び戻る。再び空間領域へ戻された画像データ2
601’にはインパルス信号が付加された影響があるは
ずであるが、人間の目には知覚しにくく、実質的には画
像データ2601と画像データ2601’は同一物に見
える。これは、周波数領域で加えられたインパルス信号
2603が逆フーリエ変換によって、画像データ全体に
小さな振幅で分布するからである。
【0233】図26の2603の様なインパルス信号が
加えられた場合は、ある一定の周波数成分を持つ画像デ
ータが空間領域に加えられたことに等しくなる。この加
えられたインパルス信号が人間の知覚できる周波数より
も大きく、更に振幅が人間の知覚できる限度以下である
ならば、加えられたインパルス信号は人間の目には見え
ない。よって、上記レジストレーション信号の埋め込み
自体も一種の電子透かし処理であると言える。
【0234】なお、本実施の形態では画像データ260
1にレジストレーション信号2603が埋め込まれ、更
に実際に埋め込むべき付加情報Infが埋め込まれた後、
空間領域の信号2601’を復元する。
【0235】図26の様に埋め込まれたレジストレーシ
ョン信号は、抽出の際に再びフーリエ変換が施される。
これによって、空間領域では一度画像データ全体に拡散
されたレジストレーション信号2603が、周波数領域
に変換され再びインパルス信号として現れる。
【0236】電子透かし情報が埋め込まれた画像がJPEG
圧縮等の非可逆圧縮などの攻撃を受けた場合、このイン
パルスは振幅が小さくなる可能性が高い。一方で、スケ
ーリングなどの幾何的な攻撃を受けた場合、このインパ
ルスはその位置が移動する。いずれの場合も上述した様
な適当なインパルス抽出処理を施すことによって、イン
パルス信号は抽出可能であり、元の画像データからの変
化を推測できる。この変化を補正すれば本実施の形態で
埋め込まれる付加情報Infを確実に抽出できる状態を作
り出せる。
【0237】以上の処理により、図7のインパルス抽出
手段0703からは上述したインパルス信号が出力さ
れ、スケーリング率算出手段0704に入力される。ス
ケーリング率算出手段0704は、入力されたインパル
ス信号の座標を用いてどの様なスケーリングが施された
かを算出する。
【0238】本実施の形態の電子透かし抽出装置側で
は、予めどの周波数成分にインパルス信号を埋め込んだ
かを知っているものとする。この場合には、この予め埋
め込まれた周波数と、インパルスが検出された周波数の
比によりスケーリング率を算出することが可能である。
例えば、予めインパルス信号が埋め込まれた周波数を
a、検出されたインパルス信号の周波数をbとすると、a/
b倍のスケーリングが施されていることがわかる。これ
は良く知られたフーリエ変換の性質である。以上の処理
により、スケーリング率算出手段0704からはスケー
リング率が出力される。
【0239】しかしながら本発明はこれに限らず、電子
透かし埋め込み装置側から必要に応じてレジストレーシ
ョン信号を埋め込んだ位置(周波数)の情報を受信する
様にしても良い。例えばこの位置情報は暗号化信号とし
て受信し、上記スケーリング率の算出処理を行う形態も
本発明の範疇に含まれる。こうすることによって、付加
情報Infを正しく抽出できるのはレジストレーション信
号を知っている人だけとなる。この場合、レジストレー
ション信号を付加情報Infを抽出する為の鍵として使用
できる。
【0240】スケーリング率算出手段0704から出力
されたスケーリング率は、スケーリング手段0705に
入力される。スケーリング手段0705には画像データ
wI1'も入力され、画像データwI1'は入力されたスケーリ
ング率によってスケーリング処理が施される。スケーリ
ング処理はバイリニア補間やバイキュービック補間など
種々のものが適応可能である。そして、スケーリング手
段0705からはスケーリング処理が施された画像デー
タwI2'が出力される。
【0241】[3−5 付加情報抽出処理]次に、図1
の付加情報埋め込み手段0103で付加情報Infが埋め
込まれた画像データwI’の青色成分からこの付加情報In
fを抽出する図2の付加情報抽出手段0203の動作に
ついて述べる。
【0242】この付加情報Infの抽出処理のブロック図
を図20に示す。
【0243】[3−5−1 埋め込み位置決定処理]図
20に示す様に、まず埋め込み位置決定手段2001に
おいて、画像データwI2'(青色成分)中のどの領域から
付加情報Infを抽出するかを決定する。この埋め込み位
置決定手段2001によってなされる動作は、前述した
埋め込み位置決定手段0103と同じであり、その為、
0103と2001によって決定される領域は同一のも
のとなる。
【0244】決定された領域から、前述した対応表2を
用い、更に図9に示されるパターン配列を用いて付加情
報Infが抽出される。
【0245】ここで、付加情報Infの抽出は、決定され
た領域に対してパターン配列を畳み込むことによって実
現される。
【0246】[3−5−2 信頼度距離演算手段]信頼度
距離dは埋め込んだ情報を抽出する際に必要となる計算
値である。
【0247】各ビット情報に対応する信頼度距離dを求
める方法を図6に示す。
【0248】まず始めに、図中の畳み込み演算手段06
01で行う処理を図21及び図22を用いて説明する。
【0249】図21及び図22に、付加情報Infを構成
する1ビットの情報を抽出する例を示す。
【0250】図21は付加情報Infを構成するある1ビ
ット情報が埋め込まれた画像データ(青色成分)I''(x,
y)に対してこの1ビット情報の抽出処理を行った例、そ
して図22は上記1ビット情報が埋め込まれていない画
像データI''(x,y)に対して1ビット情報の抽出処理を行
おうとした例である。
【0251】図21において、I''(x,y)が1ビット情報
が埋め込まれた画像データ、P(x,y)が畳み込み処理に用
いられる8×8のパターン配列(付加情報Inf抽出用の
パターン配列)である。この8×8のパターン配列を構
成する各要素(0,±c)は、入力画像データI''(x,y)
の同位置に配置されている画素値に積算され、更に各積
算値の和が算出される。即ち、I''(x,y)に対してP(x,y)
が畳み込まれる。ここで、I''(x,y)は、画像データI'
(x,y)が攻撃を受けた場合の画像を含んだ表現である。
攻撃を受けていない場合には、I''(x,y)=I'(x,y)であ
る。I''(x,y)に1ビット情報が埋め込まれている画像で
ある場合には、畳み込みの結果、図21に示す様に非零
の値が得られる可能性が非常に高い。特にI''(x,y)=I'
(x,y)の時には畳み込みの結果は32c2となる。
【0252】なお、本実施の形態では、埋め込みに用い
るパターン配列と抽出に用いるパターン配列は同様のも
のを用いている。しかしながら、これは本発明において
限定されるものではない。一般的には、埋め込みに用い
るパターン配列をP(x,y)、抽出に用いるパターン配列を
P'(x,y)とした場合には、 P'(x,y)=aP(x,y) という関係に変形できる。ここでaは任意の実数であ
り、本実施の形態では、簡単の為、a=1の場合について
説明する。
【0253】一方、図22に示す例では、上述の演算と
同様の演算が1ビット情報が埋め込まれていない画像デ
ータI''(x,y)に対して施されている。原画像(画像デー
タIに相当)からは畳み込み演算の結果、図22に示す
様に零の値が得られる。
【0254】以上、図21及び図22を用いて1ビット
情報の抽出方法を説明した。しかし、以上の説明は、付
加情報Infが埋め込まれる対象の画像データIにおいて畳
み込み演算の結果が0である場合であり、非常に理想的
な場合である。一方で、実際の画像データIの8×8の
パターン配列に相当する領域においては畳み込み演算の
結果が0であることはなかなか少ない。
【0255】即ち、原画像(画像データI)における8
×8のパターン配列に相当する領域について、図9のパ
ターン配列(コーンマスクも配置情報として参照)を用
いて畳み込み演算を行った場合、理想と異なり、非零の
値が算出されることもある。逆に、付加情報Infが埋め
込まれた画像(画像データwI)における8×8のパター
ン配列に相当する領域について、同じく畳み込み演算を
行った結果が“32c2“でなく“0”になってしまう
こともある。
【0256】しかしながら、付加情報Infを構成するビ
ット情報の夫々は、通常、元の画像データに複数回埋め
込まれている。即ち付加情報Infが画像に複数回埋め込
まれている。従って、各ビット情報が埋め込まれている
n個のマクロブロックにおいて、8×8パターン配列毎
の畳み込み演算を行い、ここで得られる各ビット情報毎
のn個の畳み込み演算結果に基づいて“各ビット情報が
埋め込まれているのか否か”或いは“各ビット情報が1
か0か”等を統計的に判断すれば良い。これら統計的な
判断の方法については後述する。
【0257】よって畳み込み演算手段0601は、付加
情報Infを構成する各ビット情報に1ついて、夫々複数
の畳み込み演算結果の和を求める。例えば、付加情報In
fが8ビットであれば、8個の和が得られる。この各ビ
ット情報に対応する和は平均計算手段0602に入力さ
れ、夫々が全マクロブロックの数nで割られて平均化さ
れる。この平均値が信頼度距離dである。即ち、この信
頼度距離dは、図21の“32c2”と“0”のどちらに
類似しているかを多数決的に生成した値である。
【0258】ただし、信頼度距離dは、先のパッチワー
ク法の説明ではd = 1/N Σ(ai-bi)と定義していたの
で、厳密には信頼度距離dは、P'(x,y) = 1/c P(x,y)を
用いて畳み込み演算を行った結果の平均値である。しか
しながら、P'(x,y) = aP(x,y)を用いて畳み込み演算を
行っても、畳み込み演算結果の平均値は、上記信頼度距
離dの実数倍になっているだけであり、本質的には同様
の効果が得られる。よって本発明には、P'(x,y) = aP
(x,y)を用いた畳み込み演算結果の平均値を信頼度距離d
に用いることも十分可能である。
【0259】求められた信頼度距離dは0603の記憶
媒体に蓄えられる。
【0260】畳み込み演算手段0601は、付加情報In
fを構成する各ビットについて上記信頼度距離dを繰り返
し生成し、順次記憶媒体0603に格納する。
【0261】この演算値に関してもう少し詳細な説明を
する。元の画像データIに対して図9のパターン配列
(コーンマスクも配置情報として参照)を用いて算出さ
れる信頼度距離dは理想的には0である。しかしながら
実際の画像データIにおいては、この値は非常に0に近
くはあるが非零の値が多い。各ビット情報について発生
する信頼度距離dの頻度分布を調べると、図23の様に
なる。
【0262】図23において、横軸は、各ビット情報毎
に発生する信頼度距離dの値であり、縦軸はその信頼度
距離dを生じる畳み込みが行われたビット情報の数(信
頼度距離dの出現頻度)を示している。図を見ると正規
分布に類似していることがわかる。また、元の画像デー
タIにおいては信頼度距離dは必ずしも0ではないが、そ
の平均値は0(或はそれに非常に近い値)である。
【0263】一方、元の画像データIではなく、図19
の様にビット情報“1”を埋め込んだ後の画像データ
(青色成分)をI’(x,y)に上記畳み込みを行った場合に
は、信頼度距離dは図24に示す様な頻度分布となる。
即ち、図の様に図23の分布形状を保ったまま、右方向
にシフトしている。この様に、付加情報Infを構成する
ある1ビットを埋め込んだ後の画像データは、信頼度距
離dが必ずしもcという訳ではないが、その平均値はc
(或はそれに非常に近い値)となる。
【0264】なお、図24ではビット情報“1”を埋め
込んだ例を示したが、ビット情報“0”を埋め込んだ場
合は図23に示した頻度分布が、左にシフトすることに
なる。
【0265】以上説明した様に、パッチワーク法を用い
て付加情報Inf(各ビット情報)を埋め込む場合には、
埋め込むビット数(パターン配列の使用回数)を出来る
だけ多くした方が、図23及び図24に示す様な統計的
分布が正確に現れやすい。即ち、付加情報Infを構成す
る各ビット情報が埋め込まれているか否か、或いは埋め
込まれているビット情報が“1”か“0”かを検出でき
る精度が高くなる。
【0266】[3−5−3 オフセット合わせ処理]次
にオフセット合わせ手段2002の構成について解説す
る。
【0267】オフセット合わせ手段2002には、適切
なスケーリングを施された後の画像データwI2'が入力さ
れる。この後、図6の信頼度距離演算を用い、スタート
ビットを検出する。なお、オフセット合わせ手段200
2はスタートビットInf1の5ビット分に対応する5つの
信頼度距離だけを生成する。スタートビットInf1とは、
図36に示す様に、付加情報埋め込み手段0104にお
いて予め埋め込んである付加情報Infの一部であり、本
実施の形態では5ビット分である。
【0268】なお、これらスタートビットInf1は、概念
的には最初の5ビット分であるが、付加情報Infが埋め
込まれている画像においては隣接、密集して存在する訳
ではなく、むしろ点在している。これは、対応表2のコ
ーンマスクを構成する各係数値に対応付けて順に埋め込
まれる為である。
【0269】オフセット合わせ手段2002の処理のフ
ローチャートを図28に示す。以下の解説は図28のフ
ローチャートの流れに沿って行う。
【0270】オフセット合わせ手段2002では、入力
された画像データwI2'に対して、ステップ2801によ
り、まず最も左上の座標を埋め込み開始座標と仮定す
る。同時に、最大値MAXを0に設定する。そして、ステ
ップ2802により、図6の信頼度距離演算手段を用い
て、スタートビットの検出を試みる。
【0271】ここで得られた1〜5番目のビット情報
が、ステップ2803により正しいスタートビット「1
1111」であるかどうかを判定する。この点が正しい
埋め込み開始座標であれば検出結果として5つの連続し
た正の信頼度距離dが検出されるが、そうでない場合に
は正の信頼度距離dが5つ連続しないことが多い。上記
判断を順次行い、正しいスタートビットInf1が検出でき
る位置を、埋め込み開始座標であると決定すれば良い。
【0272】しかしながら、実際のところ埋め込み開始
座標以外の点でも正しいスタートビットInf1が検出され
てしまう場合も有り得る。この原因を図27を用いて説
明する。
【0273】図27は、本実施の形態で用いるパッチワ
ーク法で埋め込まれた付加情報Infを抽出する為、付加
情報Infの埋め込み時に用いたものと同一のパターン配
列(2702、2704)(コーンマスクも配置情報と
して参照)を用いて、畳み込みを行いながら本来のマク
ロブロックの位置(2701、2703、2704)を
探索する様子を示している。左の図から右の図へ向かっ
て探索が連続的に進んでいるものとする。
【0274】図27では、簡単の為、画像データwI2'の
一部である1つのマクロブロック(付加情報Infが抽出
できる最小単位)に注目している。この図の1マスは1
ビット情報を埋め込む為のパターン配列の大きさの概念
を示している。
【0275】図27左において、2701と2702の
関係を有する場合、即ち2702が実際のマクロブロッ
ク2701よりも左上に位置して場合、元の画像と付加
情報Inf抽出用のパターン配列の位置は、斜線領域のみ
で重なっている。
【0276】また同図中央には、更に探索を進め、探索
中の位置と実際のマクロブロックの位置が完全に一致し
ている場合が示されている。この状態では、畳み込み対
象のパターン配列とマクロブロックが最大面積重なって
いる。
【0277】また同図右においては、探索中の位置が実
際に付加情報Infが埋め込まれているマクロブロックの
位置よりも右下に位置している。この状態では、この状
態では、畳み込み対象のパターン配列とマクロブロック
は斜線領域のみで重なる。
【0278】図27の全ての場合において、畳み込み対
象のパターン配列とマクロブロックが十分に重なってい
れば正しいスタートビットInf1を抽出することが可能で
ある。ただし、これら3つの場合は重なる面積が異なっ
ているので信頼度距離dが異なる。
【0279】上記重なる面積は前述した信頼度距離dに
置き換えて考えることができる。即ち、畳み込み対象の
パターン配列とマクロブロックの位置関係が完全に一致
していれば、各ビット情報共に信頼度距離dは上述した
±32cに非常に近くなる。
【0280】よって、本実施の形態においては、図28
の様に、ステップ2803において正しいスタートビッ
トInf1でないと判定された場合には、ステップ2807
によりラスタ順で次の探索点に移動する。一方で、正し
いスタートビットInf1であると判定された場合には、ス
テップ2804により、信頼度距離dが最大値MAXより大
きいかどうかを判定する。最大値MAXより小さな場合に
は、ステップ2807によりラスタ順で次の探索点に移
動する。一方で、信頼度距離dが最大値MAXよりも大きな
場合には、最大値MAXを現在の信頼度距離dに更新し、同
時に現在の探索点を埋め込み開始点として記憶する。そ
して、ステップ2806において全ての探索点を探索し
たかどうかを判定し、全て終了していない場合には、ス
テップ2807によりラスタ順で次の探索点に移動す
る。一方で、全て終了している場合には、その時記憶さ
れている埋め込み開始点を出力し処理を終了する。
【0281】以上の一連の処理により、本実施の形態に
おけるオフセット合わせ手段2002は、スタートビッ
トInf1を検出し、正しいスタートビットInf1が得られた
座標の中で、最も信頼度距離dの大きな座標の情報を、
付加情報Infの埋め込み開始点であると判断し、埋め込
み開始座標として後段へ出力する。
【0282】[3−5−4 利用情報抽出手段]利用情報
抽出手段2003は、前段のオフセット合わせ手段20
02から埋め込み開始座標、及び付加情報Infが埋め込
まれた画像データを入力し、図6で説明した動作を同じ
く用いて、ここでは利用情報Inf2を構成する各ビット情
報のみについて信頼度距離dを算出し、これらビット情
報に対する信頼度距離d1を後段の統計検定手段2006
に出力する。
【0283】なお、利用情報Inf2を構成する各ビット情
報に相当する信頼度距離d1を得ることは、実質的には、
埋め込まれた利用情報Inf2の各ビットを抽出することに
相当する。これについては後述する。
【0284】ここでは、上記探索により判別された埋め
込み開始座標に基づいて、各信頼度距離dを算出するの
みであり、スタートビットInf1の5ビット分については
抽出しない。
【0285】[3−6統計検定処理]統計検定手段200
6では、図20の利用情報抽出手段2003で得られる
信頼度距離d1の信頼性を判定する。この判定は、付加情
報Inf(利用情報Inf2)の抽出に用いた第1のパターン
配列とは異なる第2のパターン配列を用いて信頼度距離
d2を生成し、この信頼度距離d2の出現頻度分布を参照し
て信頼性指標Dを生成することで行われる。
【0286】ここで信頼度距離d1は利用情報抽出手段2
003において利用情報Inf2を抽出するために、第1の
パターン配列(コーンマスクも配置情報として参照)を
用いて得られる信頼度距離であり、信頼度距離d2は第1
のパターン配列をは異なる後述する第2のパターン配列
を用いて得られる信頼度距離である。第1のパターン配
列は、通常付加情報Inf(スタートビットInf1,利用情
報Inf2)を埋め込む際に用いた図9のパターン配列であ
る。
【0287】なお第2のパターン配列、信頼性指標D等
についての詳細は後述する。
【0288】[3―6―1 第2のパターン配列による
抽出処理] ≪ 中心極限定理 ≫{an},{bn}はn個の要素からなる画素
値の集合で、夫々図30に示される様な部分集合Aと部
分集合Bの要素の持つ画素値とする。
【0289】信頼度距離d(Σ(ai-bi) / n)は,十分な
数 nで{an}, {bn}を取る場合、画素の値anとbnには相関
がなく、信頼度距離dの期待値は0になる。また中心極
限定理より信頼度距離dの分布は独立な正規分布をと
る。
【0290】ここで中心極限定理について簡単に説明す
る。
【0291】平均値mc、標準偏差σcの母集団(正規分
布でなくても良い)から大きさncの任意標本を抽出した
時、標本平均値Scの分布はncが大きくなるにつれて正規
分布N(mc,(σc/√nc)^2)に近づくことを示す定理であ
る。
【0292】一般には母集団の標準偏差σcは不明なこ
とが多いが、サンプル数ncが十分大きく、母集団の数Nc
がサンプル数ncに比べてさらに十分大きいときは標本の
標準偏差scをσcの代わりに用いても実用上ほとんど差
し支えない。
【0293】本実施の形態に戻って説明する。まず利用
情報抽出手段2003で求められた信頼度距離d1の出現
頻度分布は、利用情報Inf2を正しく抽出できたか否かで
大きく異なる。
【0294】例えば、スタートビットInf1の検出に誤り
があった場合(オフセット合わせに失敗した場合)等に
は、利用情報Inf2が埋め込まれているはずの位置には実
際にはビット情報が埋め込まれていないので、信頼度距
離d1の出現頻度分布は図25の正規分布2501の様に
なる。
【0295】一方、正しく抽出できている場合には、利
用情報Inf2を構成するビット情報“1”に対応する各信
頼度距離d1が正規分布2502の位置に累積され、利用
情報Inf2を構成するビット情報“0”に対応する各信頼
度距離d1が正規分布2503の位置に累積される。よっ
て、この場合には2つの“山”が現れる。この2つの
“山”の大きさの比は、利用情報Inf2を構成するビット
情報“1”と“0”の比とほぼ等しい。
【0296】ただし、これは付加情報が埋め込まれてい
ない元の画像に対して第1のパターン配列で畳み込み処
理を行って得られる信頼度距離d1が、正規分布2501
の様になることを前提としたものである。
【0297】従って、現実的には、元の画像の状態を知
らない限り、正しく抽出できているか否かの判断を行う
ことは出来ない。
【0298】よって本実施の形態では付加情報が埋め込
まれていても元の画像の状態を十分判別できる、いわゆ
る第2のパターン配列を用いて、信頼度距離d2の正規分
布を生成し、この正規分布を2501として考えること
によって、利用情報Inf2が正しく抽出できているか否か
の判断を行う。
【0299】例えば、信頼度距離d2で作成した正規分布
2501を構成する斜線部分(中心から95%までの構
成要素)より外側に信頼度距離d1の出現頻度分布が存在
すれば、対象となっている画像に統計的偏りが存在し、
利用情報Inf2が埋め込まれていると確定できる。考える
ことができ、利用情報Inf2の確からしさを統計的に判断
することができる。この詳しい方法については後述す
る。
【0300】次に、付加情報Inf(利用情報Inf2)が埋
め込まれている画像データを用いて、付加情報Infが埋
め込まれる前の信頼度距離d1の出現頻度分布に類似する
もの(図25の様な正規分布2501)を生成する方法
を説明する。
【0301】本実施の形態では、第2のパターン配列に
よる抽出手段2005を用いて、正規分布2501に類
似する分布を構成する信頼度距離d2を求める。
【0302】第2のパターン配列による抽出手段200
5は、利用情報抽出手段2003に用いた第1のパター
ン配列と“直交する”第2のパターン配列を用いて、信
頼度距離d2を求める手段であり、畳み込み処理を行う点
等、利用情報抽出手段0203と動作自体はほぼ同じで
ある。
【0303】なお、対比説明の為、利用情報抽出装置2
003で用いた図9のパターン配列、及びこのパターン
配列を配置する位置を参照する為のマスク(コーンマス
ク)を、夫々「第1のパターン配列」、及び「第1の位
置参照マスク」と呼び、第1のパターン配列に“直交す
る”パターン配列、及びこのパターン配列を配置する位
置を参照する為のマスクを、夫々「第2のパターン配
列」、及び「第2の位置参照マスク」と呼ぶ。
【0304】第2のパターン配列による抽出手段200
5に、まず、オフセット合わせ手段2002から埋め込
み開始座標を入力し、上述した図6の信頼度距離演算を
用いて信頼度距離d2の計算も行う。
【0305】この時、図6の信頼度距離演算で用いるパ
ターン配列は埋め込みに用いた図9のパターン配列09
01ではなく、このパターン配列0901に“直交す
る”パターン配列3601或いは3602を用いる。
【0306】この理由は、図33のパターン配列360
1及び3602を用いて計算される信頼度距離d2には、
付加情報Infの埋め込みに用いた図9のパターン配列0
901で操作した影響がまったく反映されない為であ
る。
【0307】図34に示す様に、図9のパターン配列0
901とこれにこれに“直交する”上記パターン配列3
601とを畳み込み処理した結果は0である。これはパ
ターン配列3602についても同様である。即ち、第
1、第2のパターン配列の畳み込み結果は0である。従
って、元の画像の濃度が第1のパターン配列を用いて変
更されていたとしても、第2のパターン配列を用いて畳
み込み処理を行って得られる信頼度距離dには全く影響
が無い。
【0308】よって、付加情報Infが埋め込まれている
画像に対して上記第2のパターン配列を用いた畳み込み
処理を施して得られる信頼度距離d2の出現頻度分布は、
図25の正規分布2501とほぼ同様のものになる。従
って上記出現頻度分布を正規分布2501とみなす。
【0309】ここで得られる正規分布2501は、図3
2の3507の統計検定処理に必要な判断基準となる。
【0310】第2のパターン配列による抽出処理200
5は、上述の様に図33の3601、3602の様な
「第1のパターンとは“直交する”パターン配列」と、
図35の3802に示す第2の位置参照マスクを用い
て、信頼度距離d2の正規分布を生成するである。
【0311】なお、上記「第1のパターンとは直交する
パターン配列」の条件を以下に示すと、 (1)図33に示す様に、図9の0901と同じサイズ
であること (2)パターン配列3601、3602の様に、付加情
報Infの埋め込み時に用いた図9のパターン配列090
1との畳み込み処理の結果が0になることである。
【0312】また、図34に示す畳み込み処理は、図2
1及び図22に示される畳み込み処理と同じである。
【0313】本実施の形態では、畳み込みの結果が0に
なることを、ベクトルの内積が直交する場合に0になっ
ていることになぞらえ、「互いのパターン配列が“直交
している”」と呼ぶ。従って図33の3601、360
2は「図9のパターン配列0901に“直交する”パタ
ーン配列」である。
【0314】付加情報Infの埋め込み時に用いたパター
ン配列に“直交する”パターン配列を信頼度距離d2の計
算に用いる理由は、信頼度距離d2の分布に統計的な偏り
を存在させない、即ち0を中心の出現頻度分布を生成す
る為である。
【0315】また、「第1のパターンとは“直交する”
パターン配列」は、 (3)利用情報抽出処理2003に用いたパターン配列
の非零の要素と等しい数の非零の要素を持ち、正と負の
要素の数が夫々等しいことも必要な条件である。これは
同一の演算条件で、信頼度距離d1と信頼度距離d2が抽出
される様にする為である。
【0316】次に、本実施の形態では「第2の位置参照
マスク」は、付加情報Infの埋め込み時に用いた380
1とは別のパターンを有し、かつ3801とは異なるサ
イズの、図35の3802に示される参照マスクを用い
る。
【0317】以上、上記第1と第2のパターン配列が異
なっていれば信頼度距離d2の出現頻度分布はほぼ正規分
布2501となる。
【0318】しかしながら、スタートビットの検出位置
が完全でない場合等には、第2のパターン配列を用いて
畳み込みを行ったにも拘わらず統計的な偏りが検出され
てしまう可能性も有る。本実施の形態ではこの可能性も
考慮して、第1と第2の位置参照マスクの大きさを異な
らせることで、周期的な要素を打ち消す様にする。或い
はマスク内の各パターン配列の配置方法を異ならせるこ
とで、同一領域での畳み込みを行わない様にする。
【0319】また、この場合には「第2の位置参照マス
ク」は、これを構成する各係数がランダムに分布してい
れば良く、コーンマスクでなくとも構わない。
【0320】もし「第2の埋めこみ位置参照マスク」
が、「第1の埋めこみ位置参照マスク」と異なる様に設
定する場合には、「第2の埋めこみ位置参照マスク」は
図20の埋めこみ位置決定手段2004で作成すること
とする。
【0321】一般的には、上述した切り抜き耐性を考慮
して、第1の位置参照マスク(コーンマスク)は、付加
情報Infの埋め込み対象となる画像データ全体に対して
それほど大きなサイズを取ることは考えられない。よっ
て、「第2の位置参照マスク」は比較的大きなものを用
いると良い。本実施の形態では、付加情報Infを埋め込
む時に参照する第1のマスクよりも、付加情報Inf側で
信頼度距離d1を計算する際に用いる第2のマスクのサイ
ズが大きくなる様に設定することとする。
【0322】しかしながら本発明はこれに限らず、互い
のマスクサイズが等しくてもある程度の効果を奏する。
従って、「第2の位置参照マスク」は図20の埋めこみ
位置決定手段2001で作成されるものでも良い。
【0323】互いのマスクの最低限の条件としては、互
いのマスクに適用される付加情報Infを構成する各ビッ
トの繰り返し数が、同一サイズの画像領域内で等しいこ
とが必要である。
【0324】なお、第2のパターン配列による抽出処理
で十分な結果が得られない場合は、上述した条件を備え
る別の第2のパターン配列や第2の位置参照マスクを用
いて、再度信頼度距離d2を計算することにより、理想的
な出現頻度分布である図25の2501を生成できる可
能性も有る。
【0325】次に第2のパターン配列による抽出手段2
005の具体的な操作を示す。
【0326】本実施の形態では、第1の位置参照マスク
が32×32のコーンマスクで、第2の位置参照マスク
が64×64のコーンマスクとし、2つのマスクにおい
て、各係数の相対的な配列は全く異なっているとする。
【0327】まず、第2のパターン配列による抽出手段
2005では、抽出位置の決定は以下の対応表3に従っ
て行うこととする。
【0328】
【表3】
【0329】第2の位置参照マスクでは、同じ値の係数
がマスク内に各16個存在する。一方、32×32の第
1の位置参照マスクは、先の対応例2でマスクの参照を
行っている場合、32×32の中で同一係数の繰り返し
数は4個である。即ち、同一サイズの画像データにおい
ては、第1の位置参照マスクも第2の位置参照マスクも
同じ値の係数は同数存在する。
【0330】本実施の形態では、上記対応例3の規則に
従った位置関係に第2のパターン配列を割り当て、順次
畳み込み処理を行い、各ビット情報に対応する69個の
信頼度距離d2を算出する。
【0331】[3―6―2 信頼性指標D]第2のパター
ン配列による抽出手段2005にて生成される信頼度距
離d2は、ほぼ正規分布2501と同一の分布で出現する
が、正規分布においては、一般的に以下の式(25.1)
の範囲で95%のサンプル(信頼度距離d2)が出現する
ことが知られている。
【0332】 m-1.96σ < d2 < m+1.96σ … 式(25.1) ここで、σは上記信頼度距離d2についての標準偏差であ
り、mは平均である。
【0333】なお上記場合の範囲のことを“95%の信
頼区間”と呼ぶ。
【0334】m-1.96σ , m+1.96σ は、第2のパター
ン配列による抽出手段2005で信頼度距離d2が得られ
た後、これを用いて計算される。
【0335】利用情報抽出手段2003から統計検定手
段2006に入力される信頼度距離d1の出現頻度分布
は、ビット情報が“1”の場合は図25の正規分布25
02になり、ビット情報が“0”の場合は正規分布25
03になるので、利用情報Inf2に対応する信頼度距離d1
は、第2のパターン配列による抽出手段2005で求め
られる95%の信頼区間(図25の斜線部分)の外に存
在する確率が非常に高い。
【0336】ところで、オフセット合わせ手段2002
の処理時点で、この処理の対象となる画像に利用情報In
f2が存在しない場合には、信頼度距離d1の出現頻度分布
も正規分布2501の様になる。
【0337】この場合、利用情報Inf2に対応する64個
の信頼度距離d1の全てが式(25.1)の信頼区間に含ま
れない確率は、(1−0.95)の64乗と非常に小さ
い。
【0338】従って、信頼度距離d2に基づいて正規分布
2501を求めておけば、この正規分布の大半を占める
範囲に、信頼度距離d1に基づいて求められた出現頻度分
布が含まれるか否かを考えることにより、付加情報Inf
(利用情報Inf2)が埋め込まれているか否かをほぼ確実
に判断できる。
【0339】統計検定手段2006では、上述した様な
性質を用いて付加情報Inf(利用情報Inf2)が埋め込ま
れていることの信頼度を判断する。
【0340】本実施の形態では、利用情報Infが埋め込
まれていることの信頼度を、信頼性指標Dとして扱う。
【0341】この信頼性指標Dは、利用情報抽出手段2
003で生成する全ての信頼度距離d1における、式(2
5.1)の範囲の外に存在する信頼度距離d1の個数の割合
で定義される。
【0342】統計検定手段2006は、この信頼性指標
Dが閾値αより大きければ、信頼度距離d1の総合的な出
現頻度分布は図25の2502や2503の様な位置に
人為的に偏らされているもの、即ち利用情報Inf2が確実
に埋め込まれている画像であると判断する。
【0343】従って、ここでの判定に使用された信頼度
距離d1自体が、信頼性の有る情報であると考え、この信
頼度距離d1を更に後段の比較手段2007へ転送するこ
とを許可する。
【0344】なお、信頼性指標Dは、図32の信頼性表
示ステップ3510に示される様に、利用情報Inf2の信
頼性指標D、或いは指標Dに基づくメッセージをモニタ
等に表示しても良い。
【0345】例えば、信頼性指標Dが閾値αより大きく
ない場合は、「利用情報Inf2は正確に抽出できていませ
ん」との趣旨のメッセージを表示させ、図32の統計検
定ステップ3507から画像を再度入力するステップ3
502に戻る。
【0346】[3−7 比較処理]図20の比較手段2
007は、利用情報抽出手段2003と統計検定手段2
006とを経て出力された信頼度距離d1の値を入力す
る。ここに入力される信頼度距離d1は信頼性の高い情報
であるので、ここでは信頼度距離d1に対応する各ビット
情報が“1”と“0”の何れであるかを単純に判定する
だけで良い。
【0347】具体的には、利用情報Inf2を構成するある
ビット情報の信頼度距離d1が正であれば、このビット情
報が“1”であると判定し、信頼度距離d1が負の場合は
このビット情報が“0”であると判定する。
【0348】上記判定により得られた利用情報Inf2は、
ユーザーの参照情報、或いは制御信号にする為の最終的
なデータとして出力される。
【0349】以上で、付加情報の埋め込みから抽出まで
の一連の処理の説明を終わる。
【0350】(変形例)以上の実施の形態において、付
加情報Inf(利用情報Inf2)には誤り訂正符号化された
ものを用いることも可能であり、そうする事によって、
更に抽出された利用情報Inf2の信頼性が向上する。
【0351】なお、本発明は、複数の機器(例えばホス
トコンピュータ、インタフェース機器、リーダ、プリン
タ等)から構成されるシステムの1部として適用して
も、1つの機器(たとえば複写機、ファクシミリ装置)
からなるものの1部に適用してもよい。
【0352】また、本発明は上記実施の形態を実現する
為の装置及び方法のみに限定されるものではなく、上記
システム又は装置内のコンピュータ(CPUあるいはMPU)
に、上記実施の形態を実現する為のソフトウエアのプロ
グラムコードを供給し、このプログラムコードに従って
上記システムあるいは装置のコンピュータが上記各種デ
バイスを動作させることにより上記実施の形態を実現す
る場合も本発明の範疇に含まれる。
【0353】またこの場合、前記ソフトウエアのプログ
ラムコード自体が上記実施の形態の機能を実現すること
になり、そのプログラムコード自体、及びそのプログラ
ムコードをコンピュータに供給する為の手段、具体的に
は上記プログラムコードを格納した記憶媒体は本発明の
範疇に含まれる。
【0354】この様なプログラムコードを格納する記憶
媒体としては、例えばフロッピー(登録商標)ディス
ク、ハードディスク、光ディスク、光磁気ディスク、CD
-ROM、磁気テープ、不揮発性のメモリカード、ROM等を
用いることができる。
【0355】また、上記コンピュータが、供給されたプ
ログラムコードのみに従って各種デバイスを制御するこ
とにより、上記実施の形態の機能が実現される場合だけ
ではなく、上記プログラムコードがコンピュータ上で稼
働しているOS(オペレーティングシステム)、あるいは他
のアプリケーションソフト等と共同して上記実施の形態
が実現される場合にもかかるプログラムコードは本発明
の範疇に含まれる。
【0356】更に、この供給されたプログラムコード
が、コンピュータの機能拡張ボードやコンピュータに接
続された機能拡張ユニットに備わるメモリに格納された
後、そのプログラムコードの指示に基づいてその機能拡
張ボードや機能格納ユニットに備わるCPU等が実際の処
理の一部または全部を行い、その処理によって上記実施
の形態が実現される場合も本発明の範疇に含まれる。
【0357】なお、上記実施の形態ではコーンマスクを
用いて電子透かし情報を埋め込む場合について説明した
が、本発明はこれに限らない。特にブルーノイズマスク
を用いて電子透かし情報を埋め込む場合も本発明の範疇
に含まれる。
【0358】また、上述した種々の特徴点の少なくとも
1つを含む構成であれば本発明の範疇に含まれる。
【0359】
【発明の効果】以上説明したように本発明によれば、電
子透かし情報が埋め込まれた画像に幾何的な歪みが与え
られた場合においても、この歪みを的確に判断できる様
にし、結果的に後段の電子透かし情報の抽出を正確に行
える。
【図面の簡単な説明】
【図1】電子透かし埋め込み装置の全体構成ブロック図
【図2】電子透かし抽出装置の全体構成ブロック図
【図3】印刷系処理において抽出側で生成された画像デ
ータの一例を示す図
【図4】レジストレーション信号埋め込み手段のブロッ
ク図
【図5】レジストレーション信号を説明する図
【図6】信頼度距離演算手段を示す図
【図7】スケール合わせ手段のブロック図
【図8】レジストレーション信号の抽出を説明する図
【図9】付加情報の埋め込み時及び抽出時に用いるパタ
ーン配列を示す図
【図10】付加情報埋め込み手段のブロック図
【図11】埋め込み位置決定手段のブロック図
【図12】コーンマスク及びブルーノイズマスクの概念
【図13】人間の視覚の空間周波数特性を示す図
【図14】ブルーノイズマスク、コーンマスクの空間周
波数特性を示す図
【図15】位置参照マスクを説明する図
【図16】位置参照マスク内の埋め込み位置を示す概念
【図17】図16のマスクに各パターン配列を展開する
様子を示す図
【図18】全画像において付加情報Infの埋め込みに必
要な領域のを示す図
【図19】付加情報Infを埋め込む演算を説明する図
【図20】付加情報抽出手段を説明する図
【図21】付加情報Infを抽出する様子を説明する図
【図22】付加情報Infが存在しないにも拘わらず、抽
出しようとした様子を示す図
【図23】原画像から信頼度距離dを抽出した場合の理
想的な出現頻度分布を示す図
【図24】電子透かしが埋め込まれた画像から信頼度距
離dを抽出した場合を示す図
【図25】本実施の形態における信頼度距離d1,d2の出
現頻度分布の例を説明する図
【図26】レジストレーション信号の埋め込みと抽出の
原理を説明する図
【図27】レジストレーションを行う様子を示した図で
ある
【図28】レジストレーション処理を説明するフローチ
ャート
【図29】空間領域におけるレジストレーション信号埋
め込み手段のブロック図
【図30】パッチーワーク法における二つの集合を説明
する図
【図31】電子透かし埋め込み処理の全体を説明するフ
ローチャート
【図32】電子透かし抽出処理の全体を説明するフロー
チャート
【図33】図9のパターンに直交するパターン配列の例
を示す図
【図34】“直交する”パターン配列を説明する図
【図35】第1、第2の位置参照マスクを示した図
【図36】付加情報Infの構成を示す図
【図37】ブルーノイズマスク内の各係数の一例を示す
【図38】コーンマスクの画素値の各係数の一例を示す
【符号の説明】
0101 色成分抽出手段 0102 レジストレーション信号埋め込み手段 0103 埋め込み位置決定手段 0104 付加情報埋め込み手段 0105 色成分合成手段
【手続補正書】
【提出日】平成11年11月26日(1999.11.
26)
【手続補正1】
【補正対象書類名】明細書
【補正対象項目名】0026
【補正方法】変更
【補正内容】
【0026】レジストレーション信号埋め込み手段01
02には、色成分抽出手段0101で得られた青色成分
の画像データが入力される。レジストレーション信号埋
め込み手段0102では、画像データに対して一種の電
子透かしの技術を用いてレジストレーション信号が埋め
込まれる。即ち、レジストレーション信号を埋め込まれ
た画像データにおいて、人間の視覚はレジストレーショ
ン信号を知覚することができない。このレジストレーシ
ョン信号の埋め込み方法等、詳細は後述する。
【手続補正2】
【補正対象書類名】明細書
【補正対象項目名】0047
【補正方法】変更
【補正内容】
【0047】レジストレーション手段0202には、色
成分抽出手段0201で得られた青色成分の画像データ
wI1'が入力される。そしてこの青色成分の画像データwI
1'を用いて、幾何学的な歪みが補正された画像データwI
2'を生成する。
【手続補正3】
【補正対象書類名】明細書
【補正対象項目名】0049
【補正方法】削除
【手続補正4】
【補正対象書類名】明細書
【補正対象項目名】0050
【補正方法】変更
【補正内容】
【0050】レジストレーション手段0202は、画像
データwI2'を、付加情報抽出手段0203に出力する。
【手続補正5】
【補正対象書類名】明細書
【補正対象項目名】0051
【補正方法】変更
【補正内容】
【0051】付加情報抽出手段0203は、付加情報埋
め込み手段0104での埋め込み方式に合わせた所定の
処理を施こすことによって、画像データwI2'に埋め込ま
れている電子透かし情報を抽出することが可能であり、
抽出された付加情報Infを出力する。
【手続補正6】
【補正対象書類名】明細書
【補正対象項目名】0053
【補正方法】変更
【補正内容】
【0053】また、この画像データwI'の青色成分のみ
が抽出され、次のステップで使用される。
【手続補正7】
【補正対象書類名】明細書
【補正対象項目名】0055
【補正方法】変更
【補正内容】
【0055】更にステップ3504では、入力された青
色成分の画像データwI1'のオフセットが補正される。
【手続補正8】
【補正対象書類名】明細書
【補正対象項目名】0057
【補正方法】変更
【補正内容】
【0057】統計検定ステップ3507では、上記抽出
された電子透かし情報の確からしさが算出、判定され、
正しい電子透かし情報ではないと判定されればステップ
3502に戻り、電子透かし情報が埋め込まれていると
考えられる画像を入力しなおす。一方、充分正確な電子
透かし情報であると判定された場合には、ステップ35
08の比較処理により電子透かし情報(付加情報Inf)
の抽出を行う。また、ステップ3510では上記確から
しさを示す情報を後述する信頼性指標Dとして表示す
る。
【手続補正9】
【補正対象書類名】明細書
【補正対象項目名】0058
【補正方法】削除
【手続補正10】
【補正対象書類名】明細書
【補正対象項目名】0060
【補正方法】変更
【補正内容】
【0060】まず、電子透かし抽出側におけるレジスト
レーション手段0202、ステップ3503が実行する
レジストレーションという処理について説明する。
【手続補正11】
【補正対象書類名】明細書
【補正対象項目名】0070
【補正方法】削除
【手続補正12】
【補正対象書類名】明細書
【補正対象項目名】0094
【補正方法】変更
【補正内容】
【0094】図4では、レジストレーション信号をフー
リエ変換領域において埋め込む方式を説明した。一方
で、レジストレーション信号を空間領域において埋め込
む方式も考えられる。この方式を図29を用いて説明す
る。
【手続補正13】
【補正対象書類名】明細書
【補正対象項目名】0123
【補正方法】変更
【補正内容】
【0123】[3−2埋め込み位置決定処理]図11は埋
め込み位置決定手段0103の内部構成を記している。
【手続補正14】
【補正対象書類名】明細書
【補正対象項目名】0127
【補正方法】変更
【補正内容】
【0127】更にマスク・パターン配列対応手段110
3は、マスク内の各係数値の位置に各パターン配列の配
列要素(8×8サイズ)を展開する。即ち、図17の1
701に示されるマスクの各係数値(1マス)を、同図
1703の様に8×8倍にして各パターン配列の埋め込
み位置として参照可能とする。
【手続補正15】
【補正対象書類名】明細書
【補正対象項目名】0128
【補正方法】変更
【補正内容】
【0128】後述の付加情報埋め込み手段0104は、
図17の埋め込み先頭座標1702を参照して、パター
ン配列を用いて各ビット情報を埋め込むことになる。
【手続補正16】
【補正対象書類名】明細書
【補正対象項目名】0178
【補正方法】変更
【補正内容】
【0178】
【表1】 ここで、S1〜S5まではオフセット合わせ処理装置で用
いる位置合わせの為のビット情報(スタートビット)で
ある。1〜3は3ビットの利用情報である。
【手続補正17】
【補正対象書類名】明細書
【補正対象項目名】0198
【補正方法】変更
【補正内容】
【0198】この先頭座標から図9に示されるパターン
配列を用いると、パターン配列の大きさを持つ領域17
03に重なり合わず埋め込みを行うことが可能になる。
【手続補正18】
【補正対象書類名】明細書
【補正対象項目名】0214
【補正方法】変更
【補正内容】
【0214】図19に示す例では、I(x,y)が元の画像、
P(x,y)が8×8のパターン配列である。8×8のパター
ン配列を構成する各係数は、このパターン配列と同じサ
イズの元の画像データ(青色成分)に重ね合わせられ、
同位置の値同士が加減算される。その結果、I'(x,y)が
算出され、ビット情報が埋め込まれた青色成分の画像デ
ータとして図1の色成分合成手段0105へ出力され
る。
【手続補正19】
【補正対象書類名】明細書
【補正対象項目名】0219
【補正方法】変更
【補正内容】
【0219】以上の処理により、画像データ中に付加情
報Infが埋め込まれる。この付加情報Infが埋め込まれた
画像データの各画素が、十分少ないドット数で表現され
るのであれば、パターン配列の大きさも十分小さくなる
ので、これらパターン配列の1つ1つは非常に小さい点
でしか人間の目には知覚されない。よって、コーンマス
クの空間周波数特性も維持され、人間の目に見えにく
い。
【手続補正20】
【補正対象書類名】明細書
【補正対象項目名】0241
【補正方法】変更
【補正内容】
【0241】[3−5 付加情報抽出処理]次に、図1
の付加情報埋め込み手段0104で付加情報Infが埋め
込まれた画像データwI’の青色成分からこの付加情報In
fを抽出する図2の付加情報抽出手段0203の動作に
ついて述べる。
【手続補正21】
【補正対象書類名】明細書
【補正対象項目名】0250
【補正方法】変更
【補正内容】
【0250】図21は付加情報Infを構成するある1ビ
ット情報が埋め込まれた画像データ(青色成分)I''(x,
y)に対してこの1ビット情報の抽出処理を行った例、そ
して図22は上記1ビット情報が埋め込まれていない画
像データI(x,y)に対して1ビット情報の抽出処理を行お
うとした例である。
【手続補正22】
【補正対象書類名】明細書
【補正対象項目名】0253
【補正方法】変更
【補正内容】
【0253】一方、図22に示す例では、上述の演算と
同様の演算が1ビット情報が埋め込まれていない画像デ
ータI(x,y)に対して施されている。原画像(画像データ
Iに相当)からは畳み込み演算の結果、図22に示す様
に零の値が得られる。
【手続補正23】
【補正対象書類名】明細書
【補正対象項目名】0257
【補正方法】変更
【補正内容】
【0257】よって畳み込み演算手段0601は、付加
情報Infを構成する各ビット情報に1ついて、夫々複数
の畳み込み演算結果の和を求める。例えば、付加情報In
fが8ビットであれば、8個の和が得られる。この各ビ
ット情報に対応する和は平均計算手段0602に入力さ
れ、夫々が全マクロブロックの数nで割られて平均化さ
れる。この平均値が信頼度距離dである。即ち、この信
頼度距離dは、図21の“32c2”と“0”のどちらに
類似しているかを多数決的に生成した値である。
【手続補正24】
【補正対象書類名】明細書
【補正対象項目名】0273
【補正方法】変更
【補正内容】
【0273】図27は、本実施の形態で用いるパッチワ
ーク法で埋め込まれた付加情報Infを抽出する為、付加
情報Infの埋め込み時に用いたものと同一のパターン配
列(2702、2705)(コーンマスクも配置情報と
して参照)を用いて、畳み込みを行いながら本来のマク
ロブロックの位置(2701、2703、2704)を
探索する様子を示している。左の図から右の図へ向かっ
て探索が連続的に進んでいるものとする。
【手続補正25】
【補正対象書類名】明細書
【補正対象項目名】0279
【補正方法】変更
【補正内容】
【0279】上記重なる面積は前述した信頼度距離dに
置き換えて考えることができる。即ち、畳み込み対象の
パターン配列とマクロブロックの位置関係が完全に一致
していれば、各ビット情報共に信頼度距離dは上述した
±32c2に非常に近くなる。
【手続補正26】
【補正対象書類名】明細書
【補正対象項目名】0286
【補正方法】変更
【補正内容】
【0286】ここで信頼度距離d1は利用情報抽出手段2
003において利用情報Inf2を抽出するために、第1の
パターン配列(コーンマスクも配置情報として参照)を
用いて得られる信頼度距離であり、信頼度距離d2は第1
のパターン配列とは異なる後述する第2のパターン配列
を用いて得られる信頼度距離である。第1のパターン配
列は、通常付加情報Inf(スタートビットInf1,利用情
報Inf2)を埋め込む際に用いた図9のパターン配列であ
る。
【手続補正27】
【補正対象書類名】明細書
【補正対象項目名】0299
【補正方法】変更
【補正内容】
【0299】例えば、信頼度距離d2で作成した正規分布
2501を構成する斜線部分(中心から95%までの構
成要素)より外側に信頼度距離d1の出現頻度分布が存在
すれば、対象となっている画像に統計的偏りが存在し、
利用情報Inf2が埋め込まれていると考えることができ、
利用情報Inf2の確からしさを統計的に判断することがで
きる。この詳しい方法については後述する。
【手続補正28】
【補正対象書類名】明細書
【補正対象項目名】0302
【補正方法】変更
【補正内容】
【0302】第2のパターン配列による抽出手段200
5は、利用情報抽出手段2003に用いた第1のパター
ン配列と“直交する”第2のパターン配列を用いて、信
頼度距離d2を求める手段であり、畳み込み処理を行う点
等、利用情報抽出手段2003と動作自体はほぼ同じで
ある。
【手続補正29】
【補正対象書類名】明細書
【補正対象項目名】0314
【補正方法】変更
【補正内容】
【0314】付加情報Infの埋め込み時に用いたパター
ン配列に“直交する”パターン配列を信頼度距離d2の計
算に用いる理由は、信頼度距離d2の分布に統計的な偏り
を存在させない、即ち0が中心の出現頻度分布を生成す
る為である。
【手続補正30】
【補正対象書類名】明細書
【補正対象項目名】0330
【補正方法】変更
【補正内容】
【0330】本実施の形態では、上記対応表3の規則に
従った位置関係に第2のパターン配列を割り当て、順次
畳み込み処理を行い、各ビット情報に対応する69個の
信頼度距離d2を算出する。
【手続補正31】
【補正対象書類名】明細書
【補正対象項目名】図面の簡単な説明
【補正方法】変更
【補正内容】
【図面の簡単な説明】
【図1】電子透かし埋め込み装置の全体構成ブロック図
【図2】電子透かし抽出装置の全体構成ブロック図
【図3】印刷系処理において抽出側で生成された画像デ
ータの一例を示す図
【図4】レジストレーション信号埋め込み手段のブロッ
ク図
【図5】レジストレーション信号を説明する図
【図6】信頼度距離演算手段を示す図
【図7】スケール合わせ手段のブロック図
【図8】レジストレーション信号の抽出を説明する図
【図9】付加情報の埋め込み時及び抽出時に用いるパタ
ーン配列を示す図
【図10】付加情報埋め込み手段のブロック図
【図11】埋め込み位置決定手段のブロック図
【図12】コーンマスク及びブルーノイズマスクの係数
値の出現頻度分布を示す図
【図13】人間の視覚の空間周波数特性を示す図
【図14】ブルーノイズマスク、コーンマスクの空間周
波数特性を示す図
【図15】位置参照マスクを説明する図
【図16】位置参照マスク内の埋め込み位置を示す概念
【図17】図16のマスクに各パターン配列を展開する
様子を示す図
【図18】全画像において付加情報Infの埋め込みに必
要な領域を示す図
【図19】付加情報Infを埋め込む演算を説明する図
【図20】付加情報抽出手段を説明する図
【図21】付加情報Infを抽出する様子を説明する図
【図22】付加情報Infが存在しないにも拘わらず、抽
出しようとした様子を示す図
【図23】原画像から信頼度距離dを抽出した場合の理
想的な出現頻度分布を示す図
【図24】電子透かしが埋め込まれた画像から信頼度距
離dを抽出した場合を示す図
【図25】本実施の形態における信頼度距離d1,d2の出
現頻度分布の例を説明する図
【図26】レジストレーション信号の埋め込みと抽出の
原理を説明する図
【図27】オフセット合わせ処理を行う様子を示した図
【図28】オフセット合わせ処理を説明するフローチャ
ート
【図29】空間領域におけるレジストレーション信号埋
め込み手段のブロック図
【図30】パッチーワーク法における二つの集合を説明
する図
【図31】電子透かし埋め込み処理の全体を説明するフ
ローチャート
【図32】電子透かし抽出処理の全体を説明するフロー
チャート
【図33】図9のパターンに直交するパターン配列の例
を示す図
【図34】“直交する”パターン配列を説明する図
【図35】第1、第2の位置参照マスクを示した図
【図36】付加情報Infの構成を示す図
【図37】ブルーノイズマスク内の各係数の一例を示す
【図38】コーンマスクの画素値の各係数の一例を示す
【符号の説明】 0101 色成分抽出手段 0102 レジストレーション信号埋め込み手段 0103 埋め込み位置決定手段 0104 付加情報埋め込み手段 0105 色成分合成手段
【手続補正32】
【補正対象書類名】図面
【補正対象項目名】図25
【補正方法】変更
【補正内容】
【図25】
───────────────────────────────────────────────────── フロントページの続き (51)Int.Cl.7 識別記号 FI テーマコート゛(参考) H04N 1/40 Z (72)発明者 岩村 恵市 東京都大田区下丸子3丁目30番2号キヤノ ン株式会社内 Fターム(参考) 5B057 AA20 CA01 CA08 CA12 CA16 CD12 DA07 DB02 DB06 DB09 DC02 5C076 AA14 AA40 BA06 5C077 LL14 5L096 AA02 AA06 BA08 EA13 EA14 FA23 FA34 FA69 GA03 GA10 GA19 GA38 MA03

Claims (12)

    【特許請求の範囲】
  1. 【請求項1】 画像データの幾何学的な状態を示すレジ
    ストレーション信号、及び前記画像データに電子透かし
    情報が予め人間に目視しにくく埋め込まれているはずの
    画像データから、該電子透かし情報を抽出する為の画像
    処理装置であって、 前記画像データを周波数変換する変換手段と、 前記変換手段により得られる信号に微分処理を施す微分
    手段と、 前記微分手段で得られた信号と、前記レジストレーショ
    ン信号に関する情報とに基づいて、前記画像データの幾
    何学的な補正を行う補正手段とを有することを特徴とす
    る画像処理装置。
  2. 【請求項2】 前記幾何学的な補正はスケーリングであ
    ることを特徴とする請求項1に記載の画像処理装置。
  3. 【請求項3】 前記電子透かし情報には、前記電子透か
    しの基準位置を示す位置情報が含まれることを特徴とす
    る請求項1に記載の画像処理装置。
  4. 【請求項4】 前記電子透かし情報は、前記位置情報と
    利用情報を含み、該利用情報には前記装置のID或いは
    ユーザーIDが含まれる請求項3に記載の画像処理装
    置。
  5. 【請求項5】 前記電子透かし情報は、前記位置情報と
    利用情報を含み、該利用情報には前記装置を制御する為
    の情報が含まれる請求項3に記載の画像処理装置。
  6. 【請求項6】 更に、前記補正が施された画像データか
    ら前記電子透かし情報を抽出する抽出手段を有すること
    を特徴とする請求項1に記載の画像処理装置。
  7. 【請求項7】 前記レジストレーション信号は、前記画
    像データの周波数領域における所定位置に配置されるイ
    ンパルス信号であることを特徴とする請求項1に記載の
    画像処理装置。
  8. 【請求項8】 前記微分処理は二次微分処理であること
    を特徴とする請求項1に記載の画像処理装置。
  9. 【請求項9】 更に、前記変換手段は、前記画像データ
    の一部の色成分に対してのみ周波数変換を行うことを特
    徴とする請求項1に記載の画像処理装置。
  10. 【請求項10】 前記一部の色成分は青色成分であるこ
    とを特徴とする請求項1に記載の画像処理装置。
  11. 【請求項11】 画像データの幾何学的な状態を示すレ
    ジストレーション信号、及び前記画像データに電子透か
    し情報が予め人間に目視しにくく埋め込まれているはず
    の画像データから、該電子透かし情報を抽出する為の画
    像処理方法であって、 前記画像データを周波数変換する変換ステップと、 前記変換ステップで得られる信号に微分処理を施す微分
    ステップと、 前記微分ステップで得られた信号と、前記レジストレー
    ション信号に関する情報とに基づいて、前記画像データ
    の幾何学的な補正を行う補正ステップとを有することを
    特徴とする画像処理方法。
  12. 【請求項12】 画像データの幾何学的な状態を示すレ
    ジストレーション信号、及び前記画像データに電子透か
    し情報が予め人間に目視しにくく埋め込まれているはず
    の画像データから、該電子透かし情報を抽出する為の画
    像処理プログラムであって、 前記画像データを周波数変換する変換ステップと、 前記変換ステップで得られる信号に微分処理を施す微分
    ステップと、 前記微分ステップで得られた信号と、前記レジストレー
    ション信号に関する情報とに基づいて、前記画像データ
    の幾何学的な補正を行う補正ステップとを有することを
    特徴とする画像処理プログラムを、コンピュータから読
    み取り可能に記憶した記憶媒体。
JP29816199A 1999-10-20 1999-10-20 画像処理装置及び方法及び記憶媒体 Expired - Fee Related JP3869983B2 (ja)

Priority Applications (2)

Application Number Priority Date Filing Date Title
JP29816199A JP3869983B2 (ja) 1999-10-20 1999-10-20 画像処理装置及び方法及び記憶媒体
US09/688,108 US6826290B1 (en) 1999-10-20 2000-10-16 Image processing apparatus and method and storage medium

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP29816199A JP3869983B2 (ja) 1999-10-20 1999-10-20 画像処理装置及び方法及び記憶媒体

Publications (2)

Publication Number Publication Date
JP2001119559A true JP2001119559A (ja) 2001-04-27
JP3869983B2 JP3869983B2 (ja) 2007-01-17

Family

ID=17855998

Family Applications (1)

Application Number Title Priority Date Filing Date
JP29816199A Expired - Fee Related JP3869983B2 (ja) 1999-10-20 1999-10-20 画像処理装置及び方法及び記憶媒体

Country Status (1)

Country Link
JP (1) JP3869983B2 (ja)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR100454668B1 (ko) * 2001-12-27 2004-11-05 한국전자통신연구원 템플레이트를 포함하는 워터마크 삽입 및 검출 시스템
KR100837839B1 (ko) * 2001-10-25 2008-06-13 주식회사 케이티 속도변경 공격으로 훼손된 오디오 워터마크 복구 장치 및그 방법
US11134167B2 (en) * 2019-04-10 2021-09-28 Beyontech, Inc. Method of embedding and extracting watermark data using a correction image comparison with a noise-based image in the frequency domain

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR100837839B1 (ko) * 2001-10-25 2008-06-13 주식회사 케이티 속도변경 공격으로 훼손된 오디오 워터마크 복구 장치 및그 방법
KR100454668B1 (ko) * 2001-12-27 2004-11-05 한국전자통신연구원 템플레이트를 포함하는 워터마크 삽입 및 검출 시스템
US11134167B2 (en) * 2019-04-10 2021-09-28 Beyontech, Inc. Method of embedding and extracting watermark data using a correction image comparison with a noise-based image in the frequency domain

Also Published As

Publication number Publication date
JP3869983B2 (ja) 2007-01-17

Similar Documents

Publication Publication Date Title
JP4218920B2 (ja) 画像処理装置及び画像処理方法並びに記憶媒体
JP2001148776A (ja) 画像処理装置及び方法及び記憶媒体
US7142689B2 (en) Image processing apparatus for determining specific images
US6741758B2 (en) Image processor and image processing method
US7068809B2 (en) Segmentation in digital watermarking
JP2001218006A (ja) 画像処理装置、画像処理方法および記憶媒体
US6826290B1 (en) Image processing apparatus and method and storage medium
Kacker et al. Joint halftoning and watermarking
JP2004140764A (ja) 画像処理装置及び画像処理方法
US6993148B1 (en) Image processing apparatus and method, and storage medium
Amiri et al. Robust watermarking against print and scan attack through efficient modeling algorithm
JP4595014B2 (ja) 電子透かし埋め込み装置および検出装置
He et al. A practical print-scan resilient watermarking scheme
Chan et al. A framework of reversible color-to-grayscale conversion with watermarking feature
JP2003143388A (ja) 画像処理方法およびその装置ならびに印刷物
Chen et al. High-quality blind watermarking in halftones using random toggle approach
JP2008536380A (ja) 量子化・透かし入れの方法
JP3809310B2 (ja) 画像処理装置及び方法及び記憶媒体
JP4311698B2 (ja) 画像処理装置及び画像処理方法並びに記録媒体
JP3884891B2 (ja) 画像処理装置及び方法及び記憶媒体
JP2001119558A (ja) 画像処理装置及び方法及び記憶媒体
JP3869983B2 (ja) 画像処理装置及び方法及び記憶媒体
CN112767227B (zh) 一种可抗屏幕拍摄的图像水印方法
CN110570343B (zh) 一种基于自适应特征点提取的图像水印嵌入方法及装置
JP2004147253A (ja) 画像処理装置及び画像処理方法

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20041111

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20060622

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20060627

A521 Written amendment

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20060828

TRDD Decision of grant or rejection written
A01 Written decision to grant a patent or to grant a registration (utility model)

Free format text: JAPANESE INTERMEDIATE CODE: A01

Effective date: 20060926

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20061016

R150 Certificate of patent or registration of utility model

Free format text: JAPANESE INTERMEDIATE CODE: R150

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20091020

Year of fee payment: 3

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20101020

Year of fee payment: 4

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20101020

Year of fee payment: 4

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20111020

Year of fee payment: 5

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20111020

Year of fee payment: 5

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20121020

Year of fee payment: 6

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20131020

Year of fee payment: 7

LAPS Cancellation because of no payment of annual fees