(57)【要約】【課題】ドキュメントの一部が切り取られて一部の基準点が失われても、電子透かしの同期回復を行うことができるとともに、基準点を規則的に配置して、秘匿性の劣化を招くことなく、電子透かしの同期回復を行うことができるようにする。【解決手段】複数の基準点の対応関係に基づいてオリジナルの紙面画像と印刷紙面画像間の幾何学的変化を示す幾何学変換パラメータを算出する幾何学変換パラメータ算出部１０を設け、印刷紙面画像補正部１１が幾何学変換パラメータ算出部１０により算出された幾何学変換パラメータを用いて、印刷紙面画像を幾何学的に変換することで、印刷紙面画像を幾何学的に補正する。【選択図】図１
2010103925.tif 000002
【特許請求の範囲】
【請求項１】
電子透かしの埋め込み位置の特定用として、オリジナルの紙面画像に埋め込まれている複数の基準点の位置を取得するオリジナル紙面画像基準点位置取得手段と、上記オリジナルの紙面画像が印刷された紙面上の印刷紙面画像に埋め込まれている複数の基準点の位置を算出する印刷紙面画像基準点位置算出手段と、上記オリジナル紙面画像基準点位置取得手段により取得された基準点間を結ぶ線分の長さと上記印刷紙面画像基準点位置算出手段により算出された基準点間を結ぶ線分の長さとの比である線分比を算出する線分比算出手段と、上記オリジナル紙面画像基準点位置取得手段により取得された基準点間を結ぶ線分がなす方向角と上記印刷紙面画像基準点位置算出手段により算出された基準点間を結ぶ線分がなす方向角との差分である方向角差分を算出する方向角差分算出手段と、上記線分比算出手段により算出された線分比と上記方向角差分算出手段により算出された方向角差分から、上記オリジナルの紙面画像に埋め込まれている複数の基準点と上記印刷紙面画像に埋め込まれている複数の基準点との対応関係を決定する対応関係決定手段と、上記対応関係決定手段により決定された複数の基準点の対応関係に基づいて上記オリジナルの紙面画像と上記印刷紙面画像間の幾何学的変化を示す幾何学変換パラメータを算出する幾何学変換パラメータ算出手段と、上記幾何学変換パラメータ算出手段により算出された幾何学変換パラメータを用いて、上記印刷紙面画像を幾何学的に変換することで、上記印刷紙面画像を幾何学的に補正する印刷紙面画像補正手段とを備えた電子透かしの同期回復装置。
【請求項２】
オリジナルの紙面画像に対して、複数の基準点がランダムに配置されていることを特徴とする請求項１記載の電子透かしの同期回復装置。
【請求項３】
電子透かしの埋め込み位置の特定用として、オリジナルの紙面画像に埋め込まれている複数の基準点の位置を取得するオリジナル紙面画像基準点位置取得処理手順と、上記オリジナルの紙面画像が印刷された紙面上の印刷紙面画像に埋め込まれている複数の基準点の位置を算出する印刷紙面画像基準点位置算出処理手順と、上記オリジナル紙面画像基準点位置取得処理手順により取得された基準点間を結ぶ線分の長さと上記印刷紙面画像基準点位置算出処理手順により算出された基準点間を結ぶ線分の長さとの比である線分比を算出する線分比算出処理手順と、上記オリジナル紙面画像基準点位置取得処理手順により取得された基準点間を結ぶ線分がなす方向角と上記印刷紙面画像基準点位置算出処理手順により算出された基準点間を結ぶ線分がなす方向角との差分である方向角差分を算出する方向角差分算出処理手順と、上記線分比算出処理手順により算出された線分比と上記方向角差分算出処理手順により算出された方向角差分から、上記オリジナルの紙面画像に埋め込まれている複数の基準点と上記印刷紙面画像に埋め込まれている複数の基準点との対応関係を決定する対応関係決定処理手順と、上記対応関係決定処理手順により決定された複数の基準点の対応関係に基づいて上記オリジナルの紙面画像と上記印刷紙面画像間の幾何学的変化を示す幾何学変換パラメータを算出する幾何学変換パラメータ算出処理手順と、上記幾何学変換パラメータ算出処理手順により算出された幾何学変換パラメータを用いて、上記印刷紙面画像を幾何学的に変換することで、上記印刷紙面画像を幾何学的に補正する印刷紙面画像補正処理手順とをコンピュータに実行させるための電子透かしの同期回復プログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
この発明は、ドキュメントに埋め込まれている複数の基準点の位置を検出して、電子透かしの埋め込み位置の特定を可能にする電子透かしの同期回復装置及び電子透かしの同期回復プログラムに関するものである。
【背景技術】
【０００２】
従来より、例えば、印刷物などのドキュメントにうすくパターンを重畳し、そのパターンの微視的な形状に情報を秘匿することにより、印刷物のセキュリティを向上させる手法が存在する。
例えば、以下の特許文献１には、ドットでパターンを形成し、特定位置のドットの有無によって、バイナリ情報を秘匿する手法が開示されている。
また、以下の特許文献２には、情報ユニットと呼ばれる複数種類の波形を有するブロックをドキュメント上に配置して印刷し、どの種類のブロックが特定位置に存在するかによって情報を埋め込む手法が開示されている。この情報を埋め込んだパターンは、しばしば電子透かしと呼ばれる。
【０００３】
これらの手法は、ドキュメントの位置と特定のドット（または、ブロック）を対応させて情報を埋め込んでいるので、埋め込んだ情報の検出を行うには、埋め込み時の位置情報が回復されなければならない。
即ち、埋め込み時と検出時では、同じ座標軸を用いて、埋め込み処理と検出処理を行う必要がある。これを電子透かしの同期回復と称される。
【０００４】
電子透かしの同期回復を行う手法として、多くの手法が提案されているが、代表的な手法は、同期回復用のパターンを電子透かしと同時に埋め込む方法である。
以下の特許文献１には、ドットの集合で矩形の情報ブロックを構成し、この情報ブロックの四隅（または、その一部）に同期回復用の基準パターンを配置する方法が開示されている。
この情報ブロックは、ドキュメント上に繰り返し配置されるので、基準パターンも情報ブロックと同様に周期的に配置される。
従来の電子透かしの同期回復装置は、パターンマッチングの手法を用いて、周期的に配置されている基準パターンを抽出し、その基準パターンを基準にして座標を定めることによって電子透かしの同期回復を行うようにしている。
【０００５】
【特許文献１】特開平１１−１１２７８７号公報（段落番号［００１２］から［００１３］、図１３）
【特許文献１】特開２００３−２０９６７６号公報（段落番号［００１０］から［００１２］、図１）
【発明の開示】
【発明が解決しようとする課題】
【０００６】
従来の電子透かしの同期回復装置は以上のように構成されているので、ドキュメントの一部が切り取られた場合には、一部の基準パターンを抽出しても、どの部分が切り取られたかを特定することができず、電子透かしの同期回復を行うことができないことがある課題があった。
また、基準パターンを情報ブロックと同様に周期的に配置する必要があるが、人間の目は規則性に敏感であるため、周期的な基準パターンは目につきやすく、電子透かしの秘匿性が劣化することがある課題があった。
【０００７】
この発明は上記のような課題を解決するためになされたもので、ドキュメントの一部が切り取られて一部の基準点が失われても、電子透かしの同期回復を行うことができるとともに、基準点を規則的に配置して、秘匿性の劣化を招くことなく、電子透かしの同期回復を行うことができる電子透かしの同期回復装置及び電子透かしの同期回復プログラムを得ることを目的とする。
【課題を解決するための手段】
【０００８】
この発明に係る電子透かしの同期回復装置は、オリジナル紙面画像基準点位置取得手段により取得された基準点間を結ぶ線分の長さと印刷紙面画像基準点位置算出手段により算出された基準点間を結ぶ線分の長さとの比である線分比を算出する線分比算出手段と、オリジナル紙面画像基準点位置取得手段により取得された基準点間を結ぶ線分がなす方向角と印刷紙面画像基準点位置算出手段により算出された基準点間を結ぶ線分がなす方向角との差分である方向角差分を算出する方向角差分算出手段と、線分比算出手段により算出された線分比と方向角差分算出手段により算出された方向角差分から、オリジナルの紙面画像に埋め込まれている複数の基準点と印刷紙面画像に埋め込まれている複数の基準点との対応関係を決定する対応関係決定手段と、対応関係決定手段により決定された複数の基準点の対応関係に基づいてオリジナルの紙面画像と印刷紙面画像間の幾何学的変化を示す幾何学変換パラメータを算出する幾何学変換パラメータ算出手段とを設け、印刷紙面画像補正手段が幾何学変換パラメータ算出手段により算出された幾何学変換パラメータを用いて、印刷紙面画像を幾何学的に変換することで、印刷紙面画像を幾何学的に補正するようにしたものである。
【発明の効果】
【０００９】
この発明によれば、オリジナル紙面画像基準点位置取得手段により取得された基準点間を結ぶ線分の長さと印刷紙面画像基準点位置算出手段により算出された基準点間を結ぶ線分の長さとの比である線分比を算出する線分比算出手段と、オリジナル紙面画像基準点位置取得手段により取得された基準点間を結ぶ線分がなす方向角と印刷紙面画像基準点位置算出手段により算出された基準点間を結ぶ線分がなす方向角との差分である方向角差分を算出する方向角差分算出手段と、線分比算出手段により算出された線分比と方向角差分算出手段により算出された方向角差分から、オリジナルの紙面画像に埋め込まれている複数の基準点と印刷紙面画像に埋め込まれている複数の基準点との対応関係を決定する対応関係決定手段と、対応関係決定手段により決定された複数の基準点の対応関係に基づいてオリジナルの紙面画像と印刷紙面画像間の幾何学的変化を示す幾何学変換パラメータを算出する幾何学変換パラメータ算出手段とを設け、印刷紙面画像補正手段が幾何学変換パラメータ算出手段により算出された幾何学変換パラメータを用いて、印刷紙面画像を幾何学的に変換することで、印刷紙面画像を幾何学的に補正するように構成したので、ドキュメントの一部が切り取られて一部の基準点が失われても、電子透かしの同期回復を行うことができるとともに、基準点を規則的に配置して、秘匿性の劣化を招くことなく、電子透かしの同期回復を行うことができる効果がある。
【発明を実施するための最良の形態】
【００１０】
実施の形態１．
図１はこの発明の実施の形態１による電子透かしの同期回復装置を示す構成図である。
図１において、電子透かし埋め込み装置１は例えばパソコンやプリンタなどから構成されており、電子透かしをオリジナルの紙面画像に埋め込むとともに、その電子透かしの埋め込み位置の特定用として、配置が既知である複数の基準点をオリジナルの紙面画像に埋め込み、そのオリジナルの紙面画像を紙面に印刷する処理を実施する。
ここで、基準点の配置については実施の形態２で説明するが、例えば、上記の特許文献２に開示されているような基準パターンを配置するようにしてもよい（図４を参照）。
【００１１】
オリジナル紙面画像基準点位置取得部２は例えば通信器やＵＳＢポートなどの情報入力インタフェースを備えており、電子透かし埋め込み装置１からオリジナルの紙面画像に埋め込まれている複数の基準点の位置を取得する処理を実施する。なお、オリジナル紙面画像基準点位置取得部２はオリジナル紙面画像基準点位置取得手段を構成している。
因みに、オリジナルの紙面画像が文書である場合、電子透かし埋め込み装置１により基準点が電子透かしとして埋め込まれる。この場合、電子透かし埋め込み装置１から電子透かしとして埋め込まれている基準点の位置を示す位置情報を取得する。あるいは、電子透かし埋め込み装置１からオリジナルの紙面画像の画像データ（電子データ）を受け取り、その画像データから基準点の位置を算出する。
オリジナルの紙面画像がデジタル化された自然物を撮影している自然画像などである場合には、基準点が電子透かしとして埋め込まれていなくてもよい。基準点が電子透かしとして埋め込まれていない場合、適当なフィルタ等を使用して、自然画像等の画像信号から特徴点を抽出し、この特徴点を基準点として用いるようにする。
【００１２】
オリジナル紙面画像ベクトル算出部３は例えばＭＰＵなどの演算装置から構成されており、オリジナル紙面画像基準点位置取得部２により取得された複数の基準点の中から任意の２つの基準点を選択して、２つの基準点を結ぶベクトル（線分）を決定するとともに、そのベクトルの長さを算出する処理を実施する。
【００１３】
印刷紙面画像取得部４は例えばデジタルスキャナなどから構成されており、電子透かし埋め込み装置１により印刷された紙面を走査して、紙面上の印刷紙面画像を取得する処理を実施する。
印刷紙面画像基準点位置算出部５は例えば画像検出器やＭＰＵなどの演算装置から構成されており、例えばパターンマッチング法を実施して、印刷紙面画像取得部４により取得された印刷紙面画像に埋め込まれている複数の基準点を検出し、複数の基準点の位置を算出する処理を実施する。
なお、印刷紙面画像取得部４及び印刷紙面画像基準点位置算出部５から印刷紙面画像基準点位置算出手段が構成されている。
【００１４】
印刷紙面画像ベクトル算出部６は例えばＭＰＵなどの演算装置から構成されており、印刷紙面画像基準点位置算出部５により検出された複数の基準点の中から任意の２つの基準点を選択して、２つの基準点を結ぶベクトル（線分）を決定するとともに、そのベクトルの長さを算出する処理を実施する。
線分比算出部７は例えばＭＰＵなどの演算装置から構成されており、オリジナル紙面画像ベクトル算出部３により算出されたベクトルの長さと、印刷紙面画像ベクトル算出部６により算出されたベクトルの長さの比である線分比を算出する処理を実施する。
なお、オリジナル紙面画像ベクトル算出部３、印刷紙面画像ベクトル算出部６及び線分比算出部７から線分比算出手段が構成されている。
【００１５】
方向角差分算出部８は例えばＭＰＵなどの演算装置から構成されており、オリジナル紙面画像ベクトル算出部３により決定されたベクトルがなす方向角と印刷紙面画像ベクトル算出部６によりにより決定されたベクトルがなす方向角との差分である方向角差分を算出する処理を実施する。
なお、オリジナル紙面画像ベクトル算出部３、印刷紙面画像ベクトル算出部６及び方向角差分算出部８から方向角差分算出手段が構成されている。
【００１６】
基準点対応関係決定部９は例えばＭＰＵなどの演算装置から構成されており、線分比算出部７により算出された線分比と、方向角差分算出部８により算出された方向角差分とから、オリジナルの紙面画像に埋め込まれている複数の基準点と印刷紙面画像に埋め込まれている複数の基準点との対応関係を決定する処理を実施する。なお、基準点対応関係決定部９は対応関係決定手段を構成している。
【００１７】
幾何学変換パラメータ算出部１０は例えばＭＰＵなどの演算装置から構成されており、基準点対応関係決定部９により決定された複数の基準点の対応関係に基づいて、オリジナルの紙面画像と印刷紙面画像間の幾何学的変化（オリジナルの紙面画像と印刷紙面画像との間で生じた幾何学変換）を示す幾何学変換パラメータを算出する処理を実施する。なお、幾何学変換パラメータ算出部１０は幾何学変換パラメータ算出手段を構成している。
【００１８】
印刷紙面画像補正部１１は例えばＭＰＵなどの演算装置から構成されており、幾何学変換パラメータ算出部１０により算出された幾何学変換パラメータを用いて、印刷紙面画像を幾何学的に変換することで、その印刷紙面画像を幾何学的に補正する処理を実施する。
なお、印刷紙面画像補正部１１は印刷紙面画像補正手段を構成している。
【００１９】
図１の例では、電子透かしの同期回復装置の構成要素であるオリジナル紙面画像基準点位置取得部２、オリジナル紙面画像ベクトル算出部３、印刷紙面画像取得部４、印刷紙面画像基準点位置算出部５、印刷紙面画像ベクトル算出部６、線分比算出部７、方向角差分算出部８、基準点対応関係決定部９、幾何学変換パラメータ算出部１０及び印刷紙面画像補正部１１のそれぞれが専用のハードウェアで構成されているものを示しているが、電子透かしの同期回復装置がコンピュータで構成される場合、オリジナル紙面画像基準点位置取得部２、オリジナル紙面画像ベクトル算出部３、印刷紙面画像取得部４、印刷紙面画像基準点位置算出部５、印刷紙面画像ベクトル算出部６、線分比算出部７、方向角差分算出部８、基準点対応関係決定部９、幾何学変換パラメータ算出部１０及び印刷紙面画像補正部１１の処理内容を記述している電子透かしの同期回復プログラムを当該コンピュータのメモリに格納し、当該コンピュータのＣＰＵが当該メモリに格納されている電子透かしの同期回復プログラムを実行するようにしてもよい。
図２はこの発明の実施の形態１による電子透かしの同期回復装置の処理内容を示すフローチャートである。
【００２０】
次に動作について説明する。
この実施の形態１における電子透かしの同期回復装置の概略を説明すると、オリジナルの紙面画像上の２つの特徴点を結んだベクトルの長さ及び方向と、印刷紙面画像上の２つの特徴点を結んだベクトルの長さ及び方向とを比較する一般化ハフ変換を実施して、オリジナル紙面画像上の特徴点と印刷紙面画像上の特徴点との対応関係を決定することで、幾何学変換のパラメータを決定し、さらに、その幾何学変換のパラメータを用いる幾何学変換によって印刷紙面画像を補正することで、電子透かしの同期回復を行うものである。
【００２１】
ここで、ハフ変換は、画像の中から特定の図形を検出する手法であり、基本的には、パターンマッチングの手法であるが、投票によって対応点を決定するため、総当り法と比べて、処理効率が高い特徴がある。
ハフ変換は、当初、直線を検出する手法として提案されたものであるが、例えば、円など、少数のパラメータで定式化される図形の検出にも利用することができる。
一般化ハフ変換は、輪郭の各点の座標と、各点における輪郭の方向を登録することによって、任意形状の図形を検出できるようにハフ変換を拡張したものである。
この実施の形態１では、この手法をオリジナルの紙面画像に埋め込まれた電子透かしの基準点と、そのオリジナルの紙面画像が印刷された紙の印刷物をスキャンして取得した電子データである印刷紙面画像から検出した基準点との対応関係の決定に利用している。
【００２２】
以下、電子透かしの同期回復装置の処理内容を具体的に説明する。
まず、オリジナル紙面画像基準点位置取得部２は、電子透かし埋め込み装置１からオリジナルの紙面画像に埋め込まれている複数の基準点の位置を示す位置情報を取得する（ステップＳＴ１）。このとき、電子透かし埋め込み装置１からオリジナルの紙面画像の画像データを受け取り、その画像データから基準点の位置を算出するようにしてもよい。
ここでは、説明の便宜上、Ｍ個の基準点Ｐk１（ｋ1＝０，１，２，…，Ｍ−１）の位置を示す位置情報を取得するものとする。
【００２３】
オリジナル紙面画像ベクトル算出部３は、オリジナル紙面画像基準点位置取得部２がＭ個の基準点Ｐk１の位置を取得すると、Ｍ個の基準点Ｐk１の中から任意の２つの基準点Ｐpi＝（ｘpi，ｙpi），Ｐpj＝（ｘpj，ｙpj）を選択して、２つの基準点Ｐpi，Ｐpjを結ぶベクトルＰpiＰpjを決定し、下記の式（１）に示すように、そのベクトルＰpiＰpjの長さupipjを算出する（ステップＳＴ２）。図３は２つの基準点間を結ぶ線分の方向と線分の中点の座標を示す説明図である。
ベクトル              ：ＰpiＰpj＝（ｘpi−ｘpj，ｙpi−ｙpj）＝(px,py)
ベクトルＰpiＰpjの長さ：upipj＝(（ｐx)2＋(ｐy)2）1/2
（１）
なお、オリジナル紙面画像ベクトル算出部３は、全ての基準点Ｐk1（ｋ1＝０，１，２，…，Ｍ−１）の組み合わせについてベクトルの長さを算出する。ただし、同じ２つの基準点の組み合わせに対しては、１つのベクトルの長さのみを算出し、逆方向のベクトルの長さは算出しない。
【００２４】
印刷紙面画像取得部４は、電子透かし埋め込み装置１により印刷された紙面を走査して、紙面上の印刷紙面画像を示す電子データを取得する（ステップＳＴ３）。
印刷紙面画像基準点位置算出部５は、印刷紙面画像取得部４が印刷紙面画像を示す電子データを取得すると、例えば、パターンマッチング法を実施して、その印刷紙面画像に埋め込まれている複数の基準点を検出し、複数の基準点の位置を算出する（ステップＳＴ４）。
ここでは、説明の便宜上、Ｎ個の基準点Ｑk2（ｋ2＝０，１，２，…，Ｎ−１）の位置を算出するものとする。
因みに、印刷紙面画像基準点位置算出部５が基準点Ｑk2の位置を算出する際の座標系は暫定的に決定すればよいが、一般に、ハフ変換では、ドキュメントの中心に原点を取ると性能が向上する。
【００２５】
ここで、電子透かし埋め込み装置１によりオリジナルの紙面画像が紙に印刷され、その紙面を走査して、紙面上の印刷紙面画像を示す電子データを取得する場合、オリジナルの紙面画像と印刷紙面画像との間で幾何学的変化が生じる。
この実施の形態１では、幾何学的変化を示す幾何学変換パラメータを（ａ，ｂ，ｃ，ｄ，ｅ，ｆ）で表すものとする。
このとき、オリジナル紙面画像基準点位置取得部２により取得されたＭ個の基準点Ｐk1の位置（ｘk1，ｙk1）と、印刷紙面画像基準点位置算出部５により算出されたＮ個の基準点Ｑk2の位置（ｘk2，ｙk2）は、下記のように表される。
ｘk1＝ａ×ｘk2＋ｂ×ｙk2＋ｃ
ｙk1＝ｄ×ｘk2＋ｅ×ｙk2＋ｄ
なお、印刷紙面画像基準点位置算出部５によるＮ個の基準点Ｑk2の検出には一定量の誤りが含まれていてもよく、そのような誤りは、以下のステップＳＴ９における基準点の対応関係を決定する処理の段階で排除される。
【００２６】
印刷紙面画像ベクトル算出部６は、印刷紙面画像基準点位置算出部５がＮ個の基準点Ｑk2を検出すると、Ｎ個の基準点Ｑk2の中から任意の２つの基準点Ｑqi＝（ｘqi，ｙqi），Ｑqj＝（ｘqj，ｙqj）を選択して、２つの基準点Ｑqi，Ｑqjを結ぶベクトルＱqiＱqjを決定し、下記の式（２）に示すように、そのベクトルＱqiＱqjの長さｖqiqjを算出する（ステップＳＴ５）。
ベクトル              ：ＱqiＱqj＝（ｘqi-ｘqj，ｙqi-ｙqj）＝(qx,qy)
ベクトルＱqiＱqjの長さ：ｖqiqj＝(（qx)2＋(qy)2）1/2
（２）
なお、印刷紙面画像ベクトル算出部６は、全ての基準点Ｑk2（ｋ2＝０，１，２，…，Ｎ−１）の組み合わせについてベクトルの長さを算出する。ただし、同じ２つの基準点の組み合わせに対しては、１つのベクトルの長さのみを算出し、逆方向のベクトルの長さは算出しない。
【００２７】
線分比算出部７は、オリジナル紙面画像ベクトル算出部３がベクトルＰpiＰpjの長さupipjを算出し、印刷紙面画像ベクトル算出部６がベクトルＱqiＱqjの長さｖqiqjを算出すると、下記の式（３）に示すように、そのベクトルＰpiＰpjの長さupipjとベクトルＱqiＱqjの長さｖqiqjの比である線分比ＳPQを算出する（ステップＳＴ６）。
線分比：ＳPQ＝ｖqiqj／upipj                                    （３）
【００２８】
方向角差分算出部８は、オリジナル紙面画像ベクトル算出部３がベクトルＰpiＰpjを決定し、印刷紙面画像ベクトル算出部６がベクトルＱqiＱqjを決定すると、下記の式（４）に示すように、そのベクトルＰpiＰpjがなす方向角と、そのベクトルＱqiＱqjがなす方向角との差分である方向角差分θPQを算出する（ステップＳＴ７）。
内積＝ベクトルＰpiＰpj×ベクトルＱqiＱqj
＝upipj×ｖqiqjcosθPQ
＝px×qx＋py×qy
∴  θPQ＝cos-1((px×qx＋py×qy)/(upipj×ｖqiqj)
（４）
【００２９】
基準点対応関係決定部９は、線分比算出部７が線分比ＳPQを算出し、方向角差分算出部８が方向角差分θPQを算出すると、その線分比ＳPQと方向角差分θPQから、オリジナルの紙面画像に埋め込まれている複数の基準点と印刷紙面画像に埋め込まれている複数の基準点との対応関係を決定する（ステップＳＴ８）。
以下、基準点対応関係決定部９におけるオリジナルの紙面画像上の基準点Ｐk1（ｋ１＝０，１，２，…，Ｍ−１）と、印刷紙面画像上の基準点Ｑk2（ｋ２＝０，１，２，…，Ｎ−１）との対応関係の決定処理を具体的に説明する。
【００３０】
まず、基準点対応関係決定部９は、ハフ変換において投票に用いる２次元のバッファｈ（ｓ，θ）をＮ−１面分用意する。
この２次元のバッファｈは、ｋ２＝０〜Ｎ−１の各ｋ２に対して、幾何学変換パラメータｓ，θの取り得る値の（ｓ，θ）の組み合わせ分だけ用意する。従って、バッファｈは（ｋ２，ｓ，θ）の３次元バッファとみなすこともできる。なお、ｓは上記の線分比ＳPQを意味し、θは上記の方向角差分θPQを意味する。
【００３１】
次に、ｋ２＝０〜Ｎ−１のすべてのｋ２に対する２次元のバッファｈ（ｓ，θ）の値を全て“０”に初期化する。
オリジナルの紙面画像上の基準点Ｐk1（ｋ１＝０，１，２，…，Ｍ−１）において、ｋ１＝ｐｉである基準点Ｐpiから、同紙面画像上のｋ１＝ｐｉ以外のすべての基準点Ｐpj（ｐｊ＝０，１，２，…，Ｍ−１）へのベクトルＰpiＰpj（ｐｉ＝０，１，２，…，Ｍ−１であり、ｐｊはｐｉを除く、ｋ１＝０，１，２，…，Ｍ−１のいずれかである。合計（Ｍ−１）×（Ｍ−１）個）は、以下のようになる。
【００３２】
ベクトル:           Ｐ0Ｐ1，Ｐ0Ｐ2， …，Ｐ0Ｐpj， …，Ｐ0ＰM-2，Ｐ0ＰM-1←M-1個
ベクトル:  Ｐ1Ｐ0，         Ｐ1Ｐ2， …，Ｐ1Ｐpj， …，Ｐ1ＰM-2，Ｐ1ＰM-1←M-1個
ベクトル:  Ｐ2Ｐ0，  Ｐ2Ｐ1，        …，Ｐ2Ｐpj， …，Ｐ2ＰM-2，Ｐ2ＰM-1←M-1個
ベクトル:  Ｐ3Ｐ0，  Ｐ3Ｐ1，Ｐ3Ｐ2， …，Ｐ3Ｐpj， …，Ｐ3ＰM-2，Ｐ3ＰM-1←M-1個
ベクトル:  ……………………………………ＰpiＰpj…………              ←M-1個
ベクトル: ＰM-1Ｐ0，ＰM-1Ｐ1，ＰM-1Ｐ2，…，ＰM-1Ｐpj，…，ＰM-1ＰM-2       ←M-1個
【００３３】
次に、印刷紙面画像上の基準点Ｑk2（ｋ２＝０，１，２，…，Ｎ−１）において、ｋ２＝ｑｉである点Ｑqiから、同紙面画像上のｋ２＝ｑｉ以外のすべての基準点Ｑqj（ｑｊ＝０，１，２，…，Ｎ−１）へのベクトルＱqiＱqj（ｑｉ＝０，１，２，…，Ｎ−１であり、ｑｊはｑｉを除く、ｋ２＝０，１，２，…，Ｎ−１のいずれかである。合計（Ｎ−１）×（Ｎ−１）個）は、以下のようになる。
ベクトル:           Ｑ0Ｑ1，Ｑ0Ｑ2，  …，Ｑ0Ｑqj， …，Ｑ0ＱN-2，Ｑ0ＱN-1←N-1個
ベクトル:   Ｑ1Ｑ0，        Ｑ1Ｑ2，  …，Ｑ1Ｑqj， …，Ｑ1ＱN-2，Ｑ1ＱN-1←N-1個
ベクトル:   Ｑ2Ｑ0， Ｑ2Ｑ1，         …，Ｑ2Ｑqj， …，Ｑ2ＱN-2，Ｑ2ＱN-1←N-1個
ベクトル:   Ｑ3Ｑ0， Ｑ3Ｑ1，Ｑ3Ｑ2，  …，Ｑ3Ｑqj， …，Ｑ3ＱN-2，Ｑ3ＱN-1←N-1個
ベクトル:   ……………………………………ＱqiＱqj…………              ←N-1個
ベクトル: ＱN-1Ｑ0，ＱN-1Ｑ1，ＱN-1Ｑ2，…，ＱN-1Ｑqj， …，ＱN-1ＱN-2     ←N-1個
【００３４】
ここで、ｋ２をｑｉに固定し、ｑｊを０〜Ｎ−１（ｑｉを除く）で変化させたＮ−１個のベクトルＱqiＱqjに対し、ｋ１（＝ｐｉ）を０〜Ｎ−１に変化させた場合のそれぞれにおいて、ｐｊを０〜Ｍ−１（ｐｉを除く）に変化させた（Ｍ−１）×（Ｍ−１）個のベクトルＰpiＰpjとの間で得られた線分比ＳPQと方向角差分θPQとに相当する２次元のバッファｈの値を“１”だけインクリメントする（図７を参照）。
【００３５】
ｑｉを固定し、ｑｊを０〜Ｎ−１（ｑｉを除く）で変化させたＮ−１個のベクトルＱqiＱqjに対し、上記線分比ＳPQと方向角差分θPQの算出をすべての（Ｍ−１）×（Ｍ−１）個のベクトルＰpiＰpjに対して行った後、ｋ１＝０〜Ｍ−１のそれぞれの２次元バッファｈ（ｓ，θ）の値を比較し、最も大きな値を含む２次元バッファｈを持つｋ１の値がｐｉであるとき、印刷紙面画像上の基準点Ｑk2（ｋ２＝ｑｉ）は、オリジナル紙面画像上の基準点Ｐk1（ｋ１＝ｐｉ）と対応するものとする。
これを０〜Ｎ−１のすべての基準点Ｑk2に行うことで、印刷紙面画像上で検出された全ての基準点Ｑk2に対応するオリジナル紙面画像上の基準点Ｐk1を決定することができる。
【００３６】
上記の例では、単純に“最も大きな値を含む２次元バッファｈ”をもつｋ１の値を探して対応関係を決定するものについて示したが、例えば、バッファｈの値に対して、ある閾値を設定し、その閾値以上のバッファｈが見つからない場合、印刷紙面画像上の基準点Ｑk2（ｋ２＝ｑｉ）に対応するオリジナル紙面画像上の基準点Ｐk1（ｋ１＝ｐｉ）が見つからないものとしてもよい。
【００３７】
幾何学変換パラメータ算出部１０は、基準点対応関係決定部９がオリジナルの紙面画像上の基準点Ｐk1（ｋ１＝０，１，２，…，Ｍ−１）と印刷紙面画像上の基準点Ｑk2（ｋ２＝０，１，２，…，Ｎ−１）との対応関係を決定すると、それらの基準点の対応関係に基づいて、ドキュメントが通信路で生じている幾何学的変化を示す幾何学変換パラメータ（ａ，ｂ，ｃ，ｄ，ｅ，ｆ）を算出する（ステップＳＴ９）。
【００３８】
即ち、幾何学変換パラメータ算出部１０は、オリジナルの紙面画像上の特徴点（ｘk1，ｙk1）と、印刷紙面画像上の特徴点（ｘk2，ｙk2）との対応関係が、以下のアフィン変換によって近似できると仮定して、幾何学的変化を示す幾何学変換パラメータ（ａ，ｂ，ｃ，ｄ，ｅ，ｆ）を決定する。
ｘk1＝ａ×ｘk2＋ｂ×ｙk2＋ｃ
ｙk1＝ｄ×ｘk2＋ｅ×ｙk2＋ｄ
変換前座標と変換後座標の二乗和の誤差は、以下のように表される。
Ｓx＝Σi=1n（ａｘk2＋ｂｙk2＋ｃ−ｘk1）2
Ｓy＝Σi=1n（ｃｘk2＋ｄｙk2＋ｅ−ｘk1）2
ただし、ｎはオリジナルの紙面画像と印刷紙面画像との間で対応がとれている特徴点の数である。
幾何学変換パラメータ算出部１０は、これらの二乗和の誤差Ｓx，Ｓyが最小になるように、最小二乗法の計算を実施して、幾何学変換パラメータ（ａ，ｂ，ｃ，ｄ，ｅ，ｆ）を決定する。
【００３９】
印刷紙面画像補正部１１は、幾何学変換パラメータ算出部１０が幾何学変換パラメータ（ａ，ｂ，ｃ，ｄ，ｅ，ｆ）を決定すると、その幾何学的パラメータ（ａ，ｂ，ｃ，ｄ，ｅ，ｆ）をアフィン変換のパラメータとして、印刷紙面画像のアフィン変換を実行することで、その印刷紙面画像を幾何学的に補正する（ステップＳＴ１０）。
また、印刷紙面画像補正部１１は、補正後の印刷紙面画像が、オリジナルの紙面画像上のどこに対応するかを決定する。
【００４０】
図示せぬ後段の電子透かし検出装置では、印刷紙面画像に埋め込まれている電子透かしを検出する際、印刷紙面画像補正部１１により補正された印刷紙面画像を参照して、電子透かしの同期回復を行う。
【００４１】
以上で明らかなように、この実施の形態１によれば、オリジナル紙面画像基準点位置取得部２により取得された基準点間を結ぶ線分の長さと印刷紙面画像基準点位置算出部５により算出された基準点間を結ぶ線分の長さとの比である線分比を算出する線分比算出部７と、オリジナル紙面画像基準点位置取得部２により取得された基準点間を結ぶ線分がなす方向角と印刷紙面画像基準点位置算出部５により算出された基準点間を結ぶ線分がなす方向角との差分である方向角差分を算出する方向角差分算出部８と、線分比算出部７により算出された線分比と方向角差分算出部８により算出された方向角差分から、オリジナルの紙面画像に埋め込まれている複数の基準点と印刷紙面画像に埋め込まれている複数の基準点との対応関係を決定する基準点対応関係決定部９と、基準点対応関係決定部９により決定された複数の基準点の対応関係に基づいてオリジナルの紙面画像と印刷紙面画像間の幾何学的変化を示す幾何学変換パラメータを算出する幾何学変換パラメータ算出部１０とを設け、印刷紙面画像補正部１１が幾何学変換パラメータ算出部１０により算出された幾何学変換パラメータを用いて、印刷紙面画像を幾何学的に変換することで、印刷紙面画像を幾何学的に補正するように構成したので、ドキュメントの一部が切り取られて一部の基準点が失われても、電子透かしの同期回復を行うことができるとともに、基準点を規則的に配置して、秘匿性の劣化を招くことなく、電子透かしの同期回復を行うことができる効果を奏する。
【００４２】
また、基準点を結ぶベクトルの方向と、そのベクトルの長さにより、一般化ハフ変換を用いて電子透かしで埋め込んだ基準点の対応関係を決定するようにしているので、基準点を規則的に配置することなく、電子透かしの同期回復を行うことができる効果を奏する。
また、幾何学変換パラメータ（ａ，ｂ，ｃ，ｄ，ｅ，ｆ）の算出が主な基準点の一致に基づいて行われるので、一部の基準点が失われている場合や、誤りがある基準点が存在する場合でも、安定的に正しい幾何学変換パラメータ（ａ，ｂ，ｃ，ｄ，ｅ，ｆ）が得られる効果を奏する。
【００４３】
なお、この実施の形態１では、基準点対応関係決定部９が、最大値のバッファｈ（ｓ，θ）に対応するオリジナルの紙面画像上の基準点Ｐk1と、印刷紙面画像上の基準点Ｑk2が対応するものとして対応関係を決定するものについて示したが、印刷紙面画像基準点位置算出部５により検出された基準点Ｑk2の位置に誤差が存在する場合、投票が真の値の周囲に分散してしまって、検出精度が低下することがある。
この検出精度が低下する問題は、２次元のバッファｈ（ｓ，θ）をスムージングすることによって軽減することができる。
さらに、２次元のバッファｈ（ｓ，θ）の±１近傍のパラメータの平均を取って、最終的な２次元のバッファｈの値とすることで、検出精度の低下を軽減することができる。
【００４４】
実施の形態２．
上記実施の形態１では、オリジナルの紙面画像に対する基準点の配置について特に言及していないが、オリジナルの紙面画像に対して、複数の基準点がランダムに配置されているようにしてもよい。
即ち、電子透かし埋め込み装置１によりオリジナルの紙面画像に埋め込まれる複数の基準点の位置を、例えば、乱数などを用いて決定することにより、オリジナルの紙面画像に対して、複数の基準点がランダムに配置されるようにしてもよい。
【００４５】
図５はランダムに配置された基準点を示す説明図である。
図５において、×は基準点の位置を示しており、複数の基準点が平面内にランダムに配置され、基準点の一部分の配置が他の部分と重なる可能性が低くなるように配置されている。
このため、印刷物の一部分が発見されたとき、その一部分がオリジナルの紙面画像のどの部分に相当するものであるかを、基準点の配置だけから高い確率で推定することができる。
【００４６】
例えば、図５の点線で囲まれた中央領域だけが切り出されたとき、その内部の基準点の配置が他の部分と重なる可能性が低い。
これに対して、図６に示すように、基準点が規則的に配置されている場合には、その一部の基準点の配置が他の部分と重なり、切り取られた断片が元のドキュメントのどの部分に相当するのかを、基準点から断定することができない。
【００４７】
以上で明らかなように、この実施の形態２によれば、オリジナルの紙面画像に対して、複数の基準点がランダムに配置されているので、ドキュメントの一部が切り取られた場合でも、電子透かしの基準点の対応関係を決定することができる効果を奏する。
【００４８】
実施の形態３．
上記実施の形態１では、電子透かし埋め込み装置１からオリジナルの紙面画像に埋め込まれている複数の基準点の位置を示す位置情報を取得するものについて示したが、電子透かし埋め込み装置１又は電子透かし埋め込み装置１の保管メモリなどから、電子データであるオリジナルの紙面画像を例えば通信路経由で受け取り、パターンマッチング法などを実施することにより、そのオリジナルの紙面画像から複数の基準点を検出するようにしてもよい。
また、電子透かし埋め込み装置１が複数の基準点を埋め込む際に実施する基準点の埋め込み位置決定アルゴリズムと同じアルゴリズムを実施することで、オリジナルの紙面画像に埋め込まれている複数の基準点の位置を決定するようにしてもよい。
【図面の簡単な説明】
【００４９】
【図１】この発明の実施の形態１による電子透かしの同期回復装置を示す構成図である。
【図２】この発明の実施の形態１による電子透かしの同期回復装置の処理内容を示すフローチャートである。
【図３】２つの基準点間を結ぶ線分の方向と線分の中点の座標を示す説明図である。
【図４】基準パターンの配置を示す説明図である。
【図５】ランダムに配置された基準点を示す説明図である。
【図６】規則的に配置された基準点を示す説明図である。
【図７】オリジナルの紙面画像上の基準点と印刷紙面画像上の基準点と対応関係の決定方法を示す説明図である。
【符号の説明】
【００５０】
１  電子透かし埋め込み装置、２  オリジナル紙面画像基準点位置取得部（オリジナル紙面画像基準点位置取得手段）、３  オリジナル紙面画像ベクトル算出部（線分比算出手段、方向角差分算出手段）、４  印刷紙面画像取得部（印刷紙面画像基準点位置算出手段）、５  印刷紙面画像基準点位置算出部（印刷紙面画像基準点位置算出手段）、６  印刷紙面画像ベクトル算出部（線分比算出手段、方向角差分算出手段）、７  線分比算出部（線分比算出手段）、８  方向角差分算出部（方向角差分算出手段）、９  基準点対応関係決定部（対応関係決定手段）、１０  幾何学変換パラメータ算出部（幾何学変換パラメータ算出手段）、１１  印刷紙面画像補正部（印刷紙面画像補正手段）。
【図１】
2010103925.tif 000003
【図２】
2010103925.tif 000004
【図３】
2010103925.tif 000005
【図４】
2010103925.tif 000006
【図５】
2010103925.tif 000007
【図６】
2010103925.tif 000008
【図７】
2010103925.tif 000009