(57)【要約】【課題】  視覚上電子透かしが目立つことなく画像に埋め込むことができるデータの情報量を増やす。【解決手段】  右画像に右画像を識別するための電子透かしと第１のデータＤ１を表す電子透かしを挿入し、左画像に左画像を識別するための電子透かしと第２のデータＤ２を表す電子透かしを挿入する。撮影画像を撮影右画像と撮影左画像に分割し、右画像を識別するための電子透かしが挿入されている撮影半画像（撮影右画像又は撮影左画像）から第１のデータを検出し、左画像を識別するための電子透かしが挿入されている撮影半画像（撮影左画像又は撮影右画像）から第２のデータを検出する。【選択図】        図２
2006050069.tif 000002
【特許請求の範囲】
【請求項１】
周波数領域において、円上に並ぶ第１の所定数の候補位置のうちの挿入すべき第ｊ（ｊ＝１、２、・・・、Ｎ）のデータにより相互に組み合わされる位置であって前記第１の所定数より少ない第２の所定数のものに非ゼロ周波数成分を生成するサブステップと、
前記周波数領域において、画像の回転角度を検出し、前記画像の第ｊ領域を識別するための非ゼロ周波数成分を生成するサブステップと、
生成された前記非ゼロ周波数成分を合わせて前記周波数領域から空間領域に変換するサブステップと、
原画像の第ｊ領域の全部又は一部のブロックに前記変換により得た画像を加算又は減算するサブステップと、
をＮ個のｊについて繰り返すステップを備えることを特徴とする電子透かし挿入方法。
【請求項２】
周波数領域において、画像の回転角度を検出し、前記画像の第ｊ１（ｊ１＝１、２、・・・、Ｎ１）領域を識別するための非ゼロ周波数成分を生成する第１生成サブステップと、
該第１生成サブステップで生成された前記非ゼロ周波数成分を前記周波数領域から空間領域に変換する第１変換サブステップと、
原画像の第ｊ１領域の全部又は一部のブロックに前記第１変換サブステップで得た画像を加算又は減算するサブステップと、
をＮ１個のｊ１について繰り返すステップと、
前記周波数領域において、円上に並ぶ第１の所定数の候補位置のうちの挿入すべき第ｊ２（ｊ２＝１、２、・・・、Ｎ２）のデータにより相互に組み合わされる位置であって前記第１の所定数より少ない第２の所定数のものに非ゼロ周波数成分を生成する第２生成サブステップと、
該第２生成サブステップで生成された前記非ゼロ周波数成分を前記周波数領域から空間領域に変換する第２変換サブステップと、
前記原画像の第ｊ２領域の全部又は一部のブロックに前記第２変換サブステップで得た画像を加算又は減算するサブステップと、
をＮ２個のｊ２について繰り返すステップと、
を備えることを特徴とする電子透かし挿入方法。
【請求項３】
周波数領域において、円上に並ぶ第１の所定数の候補位置のうちの挿入すべき第ｊ（ｊ＝１、２、・・・、Ｎ）のデータにより相互に組み合わされる位置であって前記第１の所定数より少ない第２の所定数のものに非ゼロ周波数成分を生成する第１生成サブステップと、
前記第１生成サブステップで生成された前記非ゼロ周波数成分を前記周波数領域から空間領域に変換する第１変換サブステップと、
原画像の第ｊ領域の全部又は一部のブロックに前記第１変換サブステップで得た画像を加算又は減算するサブステップと、
をＮ個のｊについて繰り返すステップと、
前記周波数領域において、画像の回転角度を検出するための非ゼロ周波数成分を生成する第２生成ステップと、
前記第２生成ステップで生成された前記非ゼロ周波数成分を前記周波数領域から空間領域に変換する第２変換ステップと、
原画像の全部又は一部のブロックに前記第２変換ステップで得た画像を加算又は減算するステップと、
を備えることを特徴とする電子透かし挿入方法。
【請求項４】
請求項１乃至３の何れか１項に記載の電子透かし挿入方法により生成された画像が印刷された印刷物。
【請求項５】
画像の第ｋ（ｋ＝１、２、・・・、Ｎ）領域より１以上の画像ブロックを得るサブステップと、
前記１以上の画像ブロックのデータを空間領域から周波数領域に変換するサブステップと、
前記１以上の画像ブロックの各周波数成分を基に、各周波数成分の振幅を得るサブステップと、
前記周波数領域において、前記各周波数成分の振幅を基に、前記画像の回転角度を検出し、前記画像の第ｊ（ｊ＝１、２、・・・Ｎ）領域を識別するための非ゼロ周波数成分の位置を検出することにより前記画像の回転角度を検出し、前記画像の第ｊ領域を識別するサブステップと、
前記周波数領域において、円周上の第１の所定数の位置であって検出された前記回転角度を回転基準としたもののうちの前記第１の所定数よりも少ない第２の所定数の何れかの位置に非ゼロ周波数成分があるか否かを異なる半径の各円毎に検出するサブステップと、
前記第２の所定数の何れかの位置にあることが検出された前記非ゼロ周波数成分の位置の組合せを基に、前記画像に挿入されている第ｊのデータを検出するサブステップと、
をＮ個のｋについて繰り返すステップと、
Ｎ個の第ｊのデータ及び順位ｊを基にデータを復元するステップと、
を備えることを特徴とする電子透かし検出方法。
【請求項６】
画像の第ｋ（ｋ＝１、２、・・・、Ｎ）領域より１以上の画像ブロックを得るサブステップと、
前記１以上の画像ブロックのデータを空間領域から周波数領域に変換するサブステップと、
前記１以上の画像ブロックの各周波数成分を基に、各周波数成分の振幅を得るサブステップと、
前記周波数領域において、前記各周波数成分の振幅を基に、前記画像の回転角度を検出し、前記画像の第ｊ（ｊ＝１、２、・・・Ｎ）領域を識別するための非ゼロ周波数成分の位置を検出し、該検出された非ゼロ周波数成分の位置及び領域番号ｋを基に、データ番号ｊを検出するサブステップと、
検出された前記非ゼロ周波数成分の位置を基に画像の回転角度を検出するサブステップと、
前記周波数領域において、円周上の第１の所定数の位置であって検出された前記回転角度を回転基準としたもののうちの前記第１の所定数よりも少ない第２の所定数の何れかの位置に非ゼロ周波数成分があるか否かを異なる半径の各円毎に検出するサブステップと、
前記第２の所定数の何れかの位置にあることが検出された前記非ゼロ周波数成分の位置の組合せを基に、前記画像に挿入されている第ｊのデータを検出するサブステップと、
をＮ個のｋについて繰り返すステップと、
Ｎ個の第ｊのデータを基にデータを復元するステップと、
を備えることを特徴とする電子透かし検出方法。
【請求項７】
画像の第ｋ（ｋ＝１、２、・・・、Ｍ×Ｎ）領域より１以上の画像ブロックを得るサブステップと、
前記１以上の画像ブロックのデータを空間領域から周波数領域に変換するサブステップと、
前記１以上の画像ブロックの各周波数成分を基に、各周波数成分の振幅を得るサブステップと、
前記周波数領域において、前記各周波数成分の振幅を基に、前記画像の回転角度を検出し、前記画像の第ｊ（ｊ＝１、２、・・・Ｎ）領域を識別するための非ゼロ周波数成分の位置を検出することにより前記画像の小回転角度を検出するサブステップと、
をＭ×Ｎ個のｋについて繰り返すステップと、
前記画像の回転角度を検出し、前記画像の第ｊ（ｊ＝１、２、・・・、Ｎ）領域を識別するための非ゼロ周波数成分のＭ×Ｎ個のｋに亘る分布に従って、ｊとｋとの対応関係及び大回転角度を求めるステップと、
前記小回転角度及び前記大回転角度を基に画像の回転角度を求めるステップと、
ｊの値が等しいＭ個の領域の間で各周波数成分の振幅を足し合わせるステップと、
ｊの値が等しいＭ個の領域毎に、足し合わされた各周波数成分の振幅を基に、前記周波数領域において、円周上の第１の所定数の位置であって検出された前記回転角度を回転基準としたもののうちの前記第１の所定数よりも少ない第２の所定数の何れかの位置に非ゼロ周波数成分があるか否かを異なる半径の各円毎に検出するステップと、
ｊの値が等しいＭ個の領域毎に、前記第２の所定数の何れかの位置にあることが検出された前記非ゼロ周波数成分の位置の組合せを基に、前記画像に挿入されている第ｊのデータを検出するステップと、
Ｎ個の第ｊのデータ及び順位ｊを基にデータを復元するステップと、
を備えることを特徴とする電子透かし検出方法。
【請求項８】
画像の第ｋ（ｋ＝１、２、・・・、Ｎ）領域より１以上の画像ブロックを得るサブステップと、
前記１以上の画像ブロックのデータを空間領域から周波数領域に変換するサブステップと、
前記１以上の画像ブロックの各周波数成分を基に、各周波数成分の振幅を得るサブステップと、
前記周波数領域において、前記各周波数成分の振幅を基に、前記画像の回転角度を検出し、前記画像の第ｓ（ｓ＝１、２、・・・Ｎ／Ｍ）領域を識別するための非ゼロ周波数成分の位置を検出することにより前記画像の小回転角度を検出するサブステップと、
をＮ個のｋについて繰り返すステップと、
前記画像の回転角度を検出し、前記画像の第ｓ（ｓ＝１、２、・・・、Ｎ）領域を識別するための非ゼロ周波数成分のＮ個のｋに亘る分布に従って、ｋとのデータ番号ｊとの対応関係及び大回転角度を求めるステップと、
前記小回転角度及び前記大回転角度を基に画像の回転角度を求めるステップと、
各周波数成分の振幅を基に、前記周波数領域において、円周上の第１の所定数の位置であって検出された前記回転角度を回転基準としたもののうちの前記第１の所定数よりも少ない第２の所定数の何れかの位置に非ゼロ周波数成分があるか否かを異なる半径の各円毎に検出するサブステップと、
前記第２の所定数の何れかの位置にあることが検出された前記非ゼロ周波数成分の位置の組合せを基に、前記画像に挿入されている第ｊのデータを検出するサブステップと、
をＮ個のｋについて繰り返すステップと、
Ｎ個の第ｊのデータ及び順位ｊを基にデータを復元するステップと、
を備えることを特徴とする電子透かし検出方法。
【請求項９】
画像より１以上の画像ブロックを得るステップと、
前記１以上の画像ブロックのデータを空間領域から周波数領域に変換するステップと、
前記１以上の画像ブロックの各周波数成分を基に、各周波数成分の振幅を得るステップと、
前記周波数領域において、前記各周波数成分の振幅を基に、前記画像の回転角度を検出するための非ゼロ周波数成分の位置を検出することにより前記画像の回転角度を検出するステップと、
前記回転角度により前記画像の回転角度を補正するステップと、
を備え、更に、
回転角度が補正された画像の第ｊ（ｊ＝１、２、・・・、Ｎ）領域より１以上の画像ブロックを得るサブステップと、
前記１以上の画像ブロックのデータを空間領域から周波数領域に変換するサブステップと、
前記１以上の画像ブロックの各周波数成分を基に、各周波数成分の振幅を得るサブステップと、
前記周波数領域において、円周上の第１の所定数の位置のうちの前記第１の所定数よりも少ない第２の所定数の何れかの位置に非ゼロ周波数成分があるか否かを異なる半径の各円毎に検出するサブステップと、
前記第２の所定数の何れかの位置にあることが検出された前記非ゼロ周波数成分の位置の組合せを基に、前記画像に挿入されている第ｊのデータを検出するサブステップと、
をＮ個のｊについて繰り返すステップと、
Ｎ個の第ｊのデータ及び順位ｊを基にデータを復元するステップと、
を備えることを特徴とする電子透かし検出方法。
【請求項１０】
周波数領域において、円上に並ぶ第１の所定数の候補位置のうちの挿入すべき第ｊ（ｊ＝１、２、・・・、Ｎ）のデータにより相互に組み合わされる位置であって前記第１の所定数より少ない第２の所定数のものに非ゼロ周波数成分を生成するステップと、
前記周波数領域において、画像の回転角度を検出し、前記画像の第ｊ領域を識別するための非ゼロ周波数成分を生成するステップと、
生成された前記非ゼロ周波数成分を合わせて前記周波数領域から空間領域に変換するステップと、
原画像の第ｊ領域の全部又は一部のブロックに前記変換により得た画像を加算又は減算するステップと、
をＮ個のｊについて繰り返す手段を備えることを特徴とする電子透かし挿入装置。
【請求項１１】
周波数領域において、画像の回転角度を検出し、前記画像の第ｊ１（ｊ１＝１、２、・・・、Ｎ１）領域を識別するための非ゼロ周波数成分を生成する第１生成ステップと、
該第１生成ステップで生成された前記非ゼロ周波数成分を前記周波数領域から空間領域に変換する第１変換ステップと、
原画像の第ｊ１領域の全部又は一部のブロックに前記第１変換ステップで得た画像を加算又は減算するステップと、
をＮ１個のｊ１について繰り返す手段と、
前記周波数領域において、円上に並ぶ第１の所定数の候補位置のうちの挿入すべき第ｊ２（ｊ２＝１、２、・・・、Ｎ２）のデータにより相互に組み合わされる位置であって前記第１の所定数より少ない第２の所定数のものに非ゼロ周波数成分を生成する第２生成ステップと、
該第２生成ステップで生成された前記非ゼロ周波数成分を前記周波数領域から空間領域に変換する第２変換ステップと、
前記原画像の第ｊ２領域の全部又は一部のブロックに前記第２変換ステップで得た画像を加算又は減算するステップと、
をＮ２個のｊ２について繰り返す手段と、
を備えることを特徴とする電子透かし挿入装置。
【請求項１２】
周波数領域において、円上に並ぶ第１の所定数の候補位置のうちの挿入すべき第ｊ（ｊ＝１、２、・・・、Ｎ）のデータにより相互に組み合わされる位置であって前記第１の所定数より少ない第２の所定数のものに非ゼロ周波数成分を生成する生成ステップと、
前記生成ステップで生成された前記非ゼロ周波数成分を前記周波数領域から空間領域に変換する変換ステップと、
原画像の第ｊ領域の全部又は一部のブロックに前記変換ステップで得た画像を加算又は減算するステップと、
をＮ個のｊについて繰り返す手段と、
前記周波数領域において、画像の回転角度を検出するための非ゼロ周波数成分を生成する生成手段と、
前記生成手段で生成された前記非ゼロ周波数成分を前記周波数領域から空間領域に変換する変換手段と、
原画像の全部又は一部のブロックに前記変換手段で得た画像を加算又は減算する手段と、
を備えることを特徴とする電子透かし挿入装置。
【請求項１３】
画像の第ｋ（ｋ＝１、２、・・・、Ｎ）領域より１以上の画像ブロックを得るステップと、
前記１以上の画像ブロックのデータを空間領域から周波数領域に変換するステップと、
前記１以上の画像ブロックの各周波数成分を基に、各周波数成分の振幅を得るステップと、
前記周波数領域において、前記各周波数成分の振幅を基に、前記画像の回転角度を検出し、前記画像の第ｊ（ｊ＝１、２、・・・Ｎ）領域を識別するための非ゼロ周波数成分の位置を検出することにより前記画像の回転角度を検出し、前記画像の第ｊ領域を識別するステップと、
前記周波数領域において、円周上の第１の所定数の位置であって検出された前記回転角度を回転基準としたもののうちの前記第１の所定数よりも少ない第２の所定数の何れかの位置に非ゼロ周波数成分があるか否かを異なる半径の各円毎に検出するステップと、
前記第２の所定数の何れかの位置にあることが検出された前記非ゼロ周波数成分の位置の組合せを基に、前記画像に挿入されている第ｊのデータを検出するステップと、
をＮ個のｋについて繰り返す手段と、
Ｎ個の第ｊのデータ及び順位ｊを基にデータを復元する手段と、
を備えることを特徴とする電子透かし検出装置。
【請求項１４】
画像の第ｋ（ｋ＝１、２、・・・、Ｎ）領域より１以上の画像ブロックを得るステップと、
前記１以上の画像ブロックのデータを空間領域から周波数領域に変換するステップと、
前記１以上の画像ブロックの各周波数成分を基に、各周波数成分の振幅を得るステップと、
前記周波数領域において、前記各周波数成分の振幅を基に、前記画像の回転角度を検出し、前記画像の第ｊ（ｊ＝１、２、・・・Ｎ）領域を識別するための非ゼロ周波数成分の位置を検出し、該検出された非ゼロ周波数成分の位置及び領域番号ｋを基に、データ番号ｊを検出するステップと、
検出された前記非ゼロ周波数成分の位置を基に画像の回転角度を検出するステップと、
前記周波数領域において、円周上の第１の所定数の位置であって検出された前記回転角度を回転基準としたもののうちの前記第１の所定数よりも少ない第２の所定数の何れかの位置に非ゼロ周波数成分があるか否かを異なる半径の各円毎に検出するステップと、
前記第２の所定数の何れかの位置にあることが検出された前記非ゼロ周波数成分の位置の組合せを基に、前記画像に挿入されている第ｊのデータを検出するステップと、
をＮ個のｋについて繰り返す手段と、
Ｎ個の第ｊのデータを基にデータを復元する手段と、
を備えることを特徴とする電子透かし検出装置。
【請求項１５】
画像の第ｋ（ｋ＝１、２、・・・、Ｍ×Ｎ）領域より１以上の画像ブロックを得るステップと、
前記１以上の画像ブロックのデータを空間領域から周波数領域に変換するステップと、
前記１以上の画像ブロックの各周波数成分を基に、各周波数成分の振幅を得るステップと、
前記周波数領域において、前記各周波数成分の振幅を基に、前記画像の回転角度を検出し、前記画像の第ｊ（ｊ＝１、２、・・・Ｎ）領域を識別するための非ゼロ周波数成分の位置を検出することにより前記画像の小回転角度を検出するステップと、
をＭ×Ｎ個のｋについて繰り返す手段と、
前記画像の回転角度を検出し、前記画像の第ｊ（ｊ＝１、２、・・・、Ｎ）領域を識別するための非ゼロ周波数成分のＭ×Ｎ個のｋに亘る分布に従って、ｊとｋとの対応関係及び大回転角度を求める手段と、
前記小回転角度及び前記大回転角度を基に画像の回転角度を求める手段と、
ｊの値が等しいＭ個の領域の間で各周波数成分の振幅を足し合わせる手段と、
ｊの値が等しいＭ個の領域毎に、足し合わされた各周波数成分の振幅を基に、前記周波数領域において、円周上の第１の所定数の位置であって検出された前記回転角度を回転基準としたもののうちの前記第１の所定数よりも少ない第２の所定数の何れかの位置に非ゼロ周波数成分があるか否かを異なる半径の各円毎に検出する手段と、
ｊの値が等しいＭ個の領域毎に、前記第２の所定数の何れかの位置にあることが検出された前記非ゼロ周波数成分の位置の組合せを基に、前記画像に挿入されている第ｊのデータを検出する手段と、
Ｎ個の第ｊのデータ及び順位ｊを基にデータを復元する手段と、
を備えることを特徴とする電子透かし検出装置。
【請求項１６】
画像の第ｋ（ｋ＝１、２、・・・、Ｎ）領域より１以上の画像ブロックを得るステップと、
前記１以上の画像ブロックのデータを空間領域から周波数領域に変換するステップと、
前記１以上の画像ブロックの各周波数成分を基に、各周波数成分の振幅を得るステップと、
前記周波数領域において、前記各周波数成分の振幅を基に、前記画像の回転角度を検出し、前記画像の第ｓ（ｓ＝１、２、・・・Ｎ／Ｍ）領域を識別するための非ゼロ周波数成分の位置を検出することにより前記画像の小回転角度を検出するステップと、
をＮ個のｋについて繰り返す手段と、
前記画像の回転角度を検出し、前記画像の第ｓ（ｓ＝１、２、・・・、Ｎ）領域を識別するための非ゼロ周波数成分のＮ個のｋに亘る分布に従って、ｋとのデータ番号ｊとの対応関係及び大回転角度を求める手段と、
前記小回転角度及び前記大回転角度を基に画像の回転角度を求める手段と、
各周波数成分の振幅を基に、前記周波数領域において、円周上の第１の所定数の位置であって検出された前記回転角度を回転基準としたもののうちの前記第１の所定数よりも少ない第２の所定数の何れかの位置に非ゼロ周波数成分があるか否かを異なる半径の各円毎に検出するステップと、
前記第２の所定数の何れかの位置にあることが検出された前記非ゼロ周波数成分の位置の組合せを基に、前記画像に挿入されている第ｊのデータを検出するステップと、
をＮ個のｋについて繰り返す手段と、
Ｎ個の第ｊのデータ及び順位ｊを基にデータを復元する手段と、
を備えることを特徴とする電子透かし検出装置。
【請求項１７】
画像より１以上の画像ブロックを得る手段と、
前記１以上の画像ブロックのデータを空間領域から周波数領域に変換する手段と、
前記１以上の画像ブロックの各周波数成分を基に、各周波数成分の振幅を得る手段と、
前記周波数領域において、前記各周波数成分の振幅を基に、前記画像の回転角度を検出するための非ゼロ周波数成分の位置を検出することにより前記画像の回転角度を検出する手段と、
前記回転角度により前記画像の回転角度を補正する手段と、
回転角度が補正された画像の第ｊ（ｊ＝１、２、・・・、Ｎ）領域より１以上の画像ブロックを得るステップと、
前記１以上の画像ブロックのデータを空間領域から周波数領域に変換するステップと、
前記１以上の画像ブロックの各周波数成分を基に、各周波数成分の振幅を得るステップと、
前記周波数領域において、円周上の第１の所定数の位置のうちの前記第１の所定数よりも少ない第２の所定数の何れかの位置に非ゼロ周波数成分があるか否かを異なる半径の各円毎に検出するステップと、
前記第２の所定数の何れかの位置にあることが検出された前記非ゼロ周波数成分の位置の組合せを基に、前記画像に挿入されている第ｊのデータを検出するステップと、
をＮ個のｊについて繰り返す手段と、
Ｎ個の第ｊのデータ及び順位ｊを基にデータを復元する手段と、
を備えることを特徴とする電子透かし検出装置。
【請求項１８】
コンピュータに請求項１乃至３の何れか１項に記載の電子透かし挿入方法を行わせるためのプログラム。
【請求項１９】
コンピュータに請求項５乃至９の何れか１項に記載の電子透かし検出方法を行わせるためのプログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、画像データに電子透かしを挿入するための電子透かし挿入方法、その装置並びにそのプログラム、電子透かしが挿入された画像データから電子透かしを検出するための電子透かし検出方法、その装置及びそのプログラムに関する。
【背景技術】
【０００２】
近年、画像データに電子透かしを挿入する技術が開発され、著作権等を有する者等を特定する情報、著作物等の利用を許諾する場合の利用方法及び条件に関する情報等を電子透かしの形態で画像データに挿入することが実用化されてきている。
【０００３】
通常は電子透かしは電子データの形態で挿入され、検出されるが、印刷にも対応した電子透かしも登場してきている。これは、電子透かしが挿入された画像データが印刷された場合であっても、印刷された画像データから電子透かしを検出できるようにしたものである。
【０００４】
ところが、従来からの印刷にも対応した電子透かしは、スキャナを用いて印刷物を読み取ることを前提としている。このような前提のもとでは、読み取られた画像の縮尺は予め知られており、画像が回転していることはない。また、スキャナは６４００ｄｐｉ程度の高解像度を有するので高周波成分に埋め込まれている電子透かしを検出することができる。更に、カメラで読み取る場合と異なり手振れが生ずる場合もない。更に、スキャナにより一定の光源により画像が照射されるので、信号対雑音比、コントラストもよい。
【０００５】
一方、近年カメラを内蔵し、電子メール機能及びＷＷＷ(World Wide Web)ブラウザ機能を有する携帯電話が普及してきている。電子透かしが挿入された画像をこのような携帯電話のカメラで読み取り、携帯電話からそのような画像を電子メールの添付ファイル又はホームページのフォーム情報として送れば、ＷＷＷサーバ側は電子透かしを検出し、それに基づいてデータベースから情報を検索し、携帯電話はＷＷＷサーバ、メールサーバ等からその様な情報の提供等の様々なサービスを受けることができる。
【０００６】
従って、携帯電話の内蔵カメラ等で読み取った画像から電子透かしを検出することを可能とする電子透かし挿入方法、その装置及びそのプログラム並びに電子透かし検出方法、その装置及びそのプログラムの登場が望まれていた。
【０００７】
しかしながら、携帯電話の内蔵カメラ等で読み取った画像は、画像とカメラの間の距離等が変動要因となり、縮尺が不定である。また、携帯電話の内蔵カメラ等で読み取った画像は、画像とカメラの間の相対回転角等が変動要因となり、少なからず回転している場合が多い。更に、大多数の携帯電話の内蔵カメラで読み取った画像は、解像度が低い。例えば、３ｃｍ×４ｃｍの画像を６４００ｄｐｉのスキャナで読み取った場合の画素数は、約７６００万であるのに対し、現在普及している大多数の携帯電話の内蔵カメラの画素数は３０万程度であることより明らかである。更に、携帯電話の内蔵カメラ等で読み取った画像は、手振れが生じている場合がある。更に、携帯電話の内蔵カメラ等で読み取った画像は、一般の外光のもとで読み取られるので、信号対雑音比及びコントラストが十分でない場合が多い。
【０００８】
そこで、縮尺が不定であり、回転していて、解像度が低く、手振れが生じていて、信号対雑音比が悪く、コントラストも悪い画像であっても電子透かしの形態で埋め込まれたデータを安定して検出することを可能とするロバストな電子透かし挿入方法、その装置及びそのプログラム並びに電子透かし検出方法、その装置およびそのプログラムを提供することを目的とした発明が発明され、特許文献１で公開されている。
【特許文献１】特開２００４−１５３９６号公報
【発明の開示】
【発明が解決しようとする課題】
【０００９】
しかしながら、特許文献１で公開されている発明によると、画像に埋め込むことができるデータの情報量がさほど多くない。また、多くの情報量のデータを画像に埋め込もうとすると、データを表す電子透かしの量を増やさなければならない。
【００１０】
そこで、本発明は、視覚上電子透かしが目立つことなく画像に埋め込むことができるデータの情報量を増やすことを可能とする電子透かし挿入方法、その装置及びそのプログラム並びに電子透かし検出方法、その装置及びそのプログラムを提供することを目的とする。
【課題を解決するための手段】
【００１１】
本発明によれば、周波数領域において、円上に並ぶ第１の所定数の候補位置のうちの挿入すべき第ｊ（ｊ＝１、２、・・・、Ｎ）のデータにより相互に組み合わされる位置であって前記第１の所定数より少ない第２の所定数のものに非ゼロ周波数成分を生成するサブステップと、前記周波数領域において、画像の回転角度を検出し、前記画像の第ｊ領域を識別するための非ゼロ周波数成分を生成するサブステップと、生成された前記非ゼロ周波数成分を合わせて前記周波数領域から空間領域に変換するサブステップと、原画像の第ｊ領域の全部又は一部のブロックに前記変換により得た画像を加算又は減算するサブステップと、をＮ個のｊについて繰り返すステップを備えることを特徴とする電子透かし挿入方法が提供される。
【００１２】
また、本発明によれば、周波数領域において、画像の回転角度を検出し、前記画像の第ｊ１（ｊ１＝１、２、・・・、Ｎ１）領域を識別するための非ゼロ周波数成分を生成する第１生成サブステップと、該第１生成サブステップで生成された前記非ゼロ周波数成分を前記周波数領域から空間領域に変換する第１変換サブステップと、原画像の第ｊ１領域の全部又は一部のブロックに前記第１変換サブステップで得た画像を加算又は減算するサブステップと、をＮ１個のｊ１について繰り返すステップと、前記周波数領域において、円上に並ぶ第１の所定数の候補位置のうちの挿入すべき第ｊ２（ｊ２＝１、２、・・・、Ｎ２）のデータにより相互に組み合わされる位置であって前記第１の所定数より少ない第２の所定数のものに非ゼロ周波数成分を生成する第２生成サブステップと、該第２生成サブステップで生成された前記非ゼロ周波数成分を前記周波数領域から空間領域に変換する第２変換サブステップと、前記原画像の第ｊ２領域の全部又は一部のブロックに前記第２変換サブステップで得た画像を加算又は減算するサブステップと、をＮ２個のｊ２について繰り返すステップと、を備えることを特徴とする電子透かし挿入方法が提供される。
【００１３】
更に、本発明によれば、周波数領域において、円上に並ぶ第１の所定数の候補位置のうちの挿入すべき第ｊ（ｊ＝１、２、・・・、Ｎ）のデータにより相互に組み合わされる位置であって前記第１の所定数より少ない第２の所定数のものに非ゼロ周波数成分を生成する第１生成サブステップと、前記第１生成サブステップで生成された前記非ゼロ周波数成分を前記周波数領域から空間領域に変換する第１変換サブステップと、原画像の第ｊ領域の全部又は一部のブロックに前記第１変換サブステップで得た画像を加算又は減算するサブステップと、をＮ個のｊについて繰り返すステップと、前記周波数領域において、画像の回転角度を検出するための非ゼロ周波数成分を生成する第２生成ステップと、前記第２生成ステップで生成された前記非ゼロ周波数成分を前記周波数領域から空間領域に変換する第２変換ステップと、原画像の全部又は一部のブロックに前記第２変換ステップで得た画像を加算又は減算するステップと、を備えることを特徴とする電子透かし挿入方法が提供される。
【００１４】
更に、本発明によれば、画像の第ｋ（ｋ＝１、２、・・・、Ｎ）領域より１以上の画像ブロックを得るサブステップと、前記１以上の画像ブロックのデータを空間領域から周波数領域に変換するサブステップと、前記１以上の画像ブロックの各周波数成分を基に、各周波数成分の振幅を得るサブステップと、前記周波数領域において、前記各周波数成分の振幅を基に、前記画像の回転角度を検出し、前記画像の第ｊ（ｊ＝１、２、・・・Ｎ）領域を識別するための非ゼロ周波数成分の位置を検出することにより前記画像の回転角度を検出し、前記画像の第ｊ領域を識別するサブステップと、前記周波数領域において、円周上の第１の所定数の位置であって検出された前記回転角度を回転基準としたもののうちの前記第１の所定数よりも少ない第２の所定数の何れかの位置に非ゼロ周波数成分があるか否かを異なる半径の各円毎に検出するサブステップと、前記第２の所定数の何れかの位置にあることが検出された前記非ゼロ周波数成分の位置の組合せを基に、前記画像に挿入されている第ｊのデータを検出するサブステップと、をＮ個のｋについて繰り返すステップと、Ｎ個の第ｊのデータ及び順位ｊを基にデータを復元するステップと、を備えることを特徴とする電子透かし検出方法が提供される。
【００１５】
更に、本発明によれば、画像の第ｋ（ｋ＝１、２、・・・、Ｎ）領域より１以上の画像ブロックを得るサブステップと、前記１以上の画像ブロックのデータを空間領域から周波数領域に変換するサブステップと、前記１以上の画像ブロックの各周波数成分を基に、各周波数成分の振幅を得るサブステップと、前記周波数領域において、前記各周波数成分の振幅を基に、前記画像の回転角度を検出し、前記画像の第ｊ（ｊ＝１、２、・・・Ｎ）領域を識別するための非ゼロ周波数成分の位置を検出し、該検出された非ゼロ周波数成分の位置及び領域番号ｋを基に、データ番号ｊを検出するサブステップと、検出された前記非ゼロ周波数成分の位置を基に画像の回転角度を検出するサブステップと、前記周波数領域において、円周上の第１の所定数の位置であって検出された前記回転角度を回転基準としたもののうちの前記第１の所定数よりも少ない第２の所定数の何れかの位置に非ゼロ周波数成分があるか否かを異なる半径の各円毎に検出するサブステップと、前記第２の所定数の何れかの位置にあることが検出された前記非ゼロ周波数成分の位置の組合せを基に、前記画像に挿入されている第ｊのデータを検出するサブステップと、をＮ個のｋについて繰り返すステップと、Ｎ個の第ｊのデータを基にデータを復元するステップと、を備えることを特徴とする電子透かし検出方法が提供される。
【００１６】
更に、本発明によれば、画像の第ｋ（ｋ＝１、２、・・・、Ｍ×Ｎ）領域より１以上の画像ブロックを得るサブステップと、前記１以上の画像ブロックのデータを空間領域から周波数領域に変換するサブステップと、前記１以上の画像ブロックの各周波数成分を基に、各周波数成分の振幅を得るサブステップと、前記周波数領域において、前記各周波数成分の振幅を基に、前記画像の回転角度を検出し、前記画像の第ｊ（ｊ＝１、２、・・・Ｎ）領域を識別するための非ゼロ周波数成分の位置を検出することにより前記画像の小回転角度を検出するサブステップと、をＭ×Ｎ個のｋについて繰り返すステップと、前記画像の回転角度を検出し、前記画像の第ｊ（ｊ＝１、２、・・・、Ｎ）領域を識別するための非ゼロ周波数成分のＭ×Ｎ個のｋに亘る分布に従って、ｊとｋとの対応関係及び大回転角度を求めるステップと、前記小回転角度及び前記大回転角度を基に画像の回転角度を求めるステップと、ｊの値が等しいＭ個の領域の間で各周波数成分の振幅を足し合わせるステップと、ｊの値が等しいＭ個の領域毎に、足し合わされた各周波数成分の振幅を基に、前記周波数領域において、円周上の第１の所定数の位置であって検出された前記回転角度を回転基準としたもののうちの前記第１の所定数よりも少ない第２の所定数の何れかの位置に非ゼロ周波数成分があるか否かを異なる半径の各円毎に検出するステップと、ｊの値が等しいＭ個の領域毎に、前記第２の所定数の何れかの位置にあることが検出された前記非ゼロ周波数成分の位置の組合せを基に、前記画像に挿入されている第ｊのデータを検出するステップと、Ｎ個の第ｊのデータ及び順位ｊを基にデータを復元するステップと、を備えることを特徴とする電子透かし検出方法が提供される。
【００１７】
更に、本発明によれば、画像の第ｋ（ｋ＝１、２、・・・、Ｎ）領域より１以上の画像ブロックを得るサブステップと、前記１以上の画像ブロックのデータを空間領域から周波数領域に変換するサブステップと、前記１以上の画像ブロックの各周波数成分を基に、各周波数成分の振幅を得るサブステップと、前記周波数領域において、前記各周波数成分の振幅を基に、前記画像の回転角度を検出し、前記画像の第ｓ（ｓ＝１、２、・・・Ｎ／Ｍ）領域を識別するための非ゼロ周波数成分の位置を検出することにより前記画像の小回転角度を検出するサブステップと、をＮ個のｋについて繰り返すステップと、前記画像の回転角度を検出し、前記画像の第ｓ（ｓ＝１、２、・・・、Ｎ）領域を識別するための非ゼロ周波数成分のＮ個のｋに亘る分布に従って、ｋとのデータ番号ｊとの対応関係及び大回転角度を求めるステップと、前記小回転角度及び前記大回転角度を基に画像の回転角度を求めるステップと、各周波数成分の振幅を基に、前記周波数領域において、円周上の第１の所定数の位置であって検出された前記回転角度を回転基準としたもののうちの前記第１の所定数よりも少ない第２の所定数の何れかの位置に非ゼロ周波数成分があるか否かを異なる半径の各円毎に検出するサブステップと、前記第２の所定数の何れかの位置にあることが検出された前記非ゼロ周波数成分の位置の組合せを基に、前記画像に挿入されている第ｊのデータを検出するサブステップと、をＮ個のｋについて繰り返すステップと、Ｎ個の第ｊのデータ及び順位ｊを基にデータを復元するステップと、を備えることを特徴とする電子透かし検出方法が提供される。
【００１８】
更に、本発明によれば、画像より１以上の画像ブロックを得るステップと、前記１以上の画像ブロックのデータを空間領域から周波数領域に変換するステップと、前記１以上の画像ブロックの各周波数成分を基に、各周波数成分の振幅を得るステップと、前記周波数領域において、前記各周波数成分の振幅を基に、前記画像の回転角度を検出するための非ゼロ周波数成分の位置を検出することにより前記画像の回転角度を検出するステップと、前記回転角度により前記画像の回転角度を補正するステップと、を備え、更に、回転角度が補正された画像の第ｊ（ｊ＝１、２、・・・、Ｎ）領域より１以上の画像ブロックを得るサブステップと、
前記１以上の画像ブロックのデータを空間領域から周波数領域に変換するサブステップと、前記１以上の画像ブロックの各周波数成分を基に、各周波数成分の振幅を得るサブステップと、前記周波数領域において、円周上の第１の所定数の位置のうちの前記第１の所定数よりも少ない第２の所定数の何れかの位置に非ゼロ周波数成分があるか否かを異なる半径の各円毎に検出するサブステップと、前記第２の所定数の何れかの位置にあることが検出された前記非ゼロ周波数成分の位置の組合せを基に、前記画像に挿入されている第ｊのデータを検出するサブステップと、をＮ個のｊについて繰り返すステップと、Ｎ個の第ｊのデータ及び順位ｊを基にデータを復元するステップと、を備えることを特徴とする電子透かし検出方法が提供される。
【発明の効果】
【００１９】
本発明によれば、データを複数に分割し、各分割データを電子透かしの形態で各分割画像に挿入するので、各分割画像に埋め込む電子透かしを少なく維持したまま、画像に埋め込むことができるデータの情報量を増やすことが可能となる。
【００２０】
また、本発明によれば、画像の回転角度を検出し、画像の各領域を識別するための電子透かしにより、どの分割画像がどのような角度でどこにあるのかを検出することができ、従って、分割データを基に元のデータを復元することができる。
【発明を実施するための最良の形態】
【００２１】
以下、図面を参照して本発明を実施するための最良の形態について詳細に説明する。
【００２２】
図１は、本発明の実施形態で用いる電子透かしの周波数領域における振幅を示す。
【００２３】
図１を参照すると、電子透かしは２種類あり、第１の電子透かしが画像の左半分（以下、「左画像」という。）に挿入され、第２の電子透かしが画像の右半分（以下、「右画像」という）に挿入される。左画像に挿入する電子透かしの周波数領域における振幅を見ると、第１の電子透かしは、外円６０１に分布する非ゼロ周波数成分及び内円６０２に分布する非ゼロ周波数成分を備える。同様に、右画像に挿入する電子透かしの周波数領域における振幅を見ると、第２の電子透かしは、外円６０３に分布する非ゼロ周波数成分及び内円６０４に分布する非ゼロ周波数成分を備える。
【００２４】
外円に分布する非ゼロ周波数成分は外円を２４等分して得られる２４個の等分点のうちのｘ軸と重なる２つの点及びｙ軸と重なる２つの点を除いた２０個の点のうちの１０個の点に配置される。非ゼロ周波数成分の振幅の分布は原点について対象であるので、２０個の点から１０個の点を選択する場合の数は、１０個の点から５個の点を選択する場合の数と等しい。従って、左画像には、１０Ｃ５＝２５２通りの第１データを挿入することができる。同様に、右画像には、１０Ｃ５＝２５２通りの第２データを挿入することができる。従って、合わせて２５２＾２＝６３５０４通りのデータを画像に挿入することができる。
【００２５】
内円６０２には、画像の回転角度を検出するための非ゼロ周波数成分が配置される。同様に、内円６０４にも、画像の回転角度を検出するための非ゼロ周波数成分が配置される。図１に示す例では、内円６０２に６個の非ゼロ周波数成分が配置され、内円６０４に６個の非ゼロ周波数成分が配置されていて、内円６０２に配置される非ゼロ周波数成分の分布と内円６０４に配置される非ゼロ周波数成分の分布は異なる。各内円に配置される非ゼロ周波数成分の振幅は、原点について点対称であるが、これは、必然的なものであり、如何に非ゼロ周波数成分を配置させても点対称性を崩すことはできず、従って、１つの円に分布する非ゼロ周波数成分だけでは画像が１８０度観点していることを検出することはできない。一般には、１８０度以外の角度だけ回転させた場合の非ゼロ周波数成分の分布が回転させない場合の非ゼロ周波数成分の分布と一致せず、且つ、内円６０２に配置される非ゼロ周波数成分の分布と内円６０４に配置される非ゼロ周波数成分の分布が一致しないという２つの条件を満たすように１又は２以上の非ゼロ周波数成分を内円６０２、６０４の各々に配置すれば、画像の１８０度の範囲の回転角度を検出できると共に、画像が逆さまであるか否か（１８０度回転しているか否か）を検出することができる。
【００２６】
符号６０５が付されている図は、回転角度を検出するための非ゼロ周波数成分を重心で表した図であり、符号６０６が付されている図は、データを記号で表した図である。以下の説明は、このような図を用いて行う。
【００２７】
図２は、図１に示すような電子透かしの印刷時における配置及び向きとカメラによる撮影時における配置及び向きを示す。
【００２８】
図２を参照すると、印刷時には、左画像の回転角度検出用電子透かしの非ゼロ周波数成分の重心は、１３５度の位置及び３１５度の位置にあり、右画像の回転角度検出用電子透かしの非ゼロ周波数成分の重心は、４５度の位置及び２２５度の位置にあり、左画像には、第２データの電子透かしが挿入され、右画像には第１データの電子透かしが挿入される。
【００２９】
カメラにより０度付近の回転角度で撮影された時には、印刷時と同様に、左画像の回転角度検出用電子透かしの非ゼロ周波数成分の重心は、１３５度の位置及び３１５度の位置にあり、右画像の回転角度検出用電子透かしの非ゼロ周波数成分の重心は、４５度の位置及び２２５度の位置にあり、左画像には、第２データの電子透かしが挿入されており、右画像には第１データの電子透かしが挿入されている。
【００３０】
一方、カメラにより１８０度付近の回転角度で撮影された時には、印刷時と異なり、左画像の回転角度検出用電子透かしの非ゼロ周波数成分の重心は、４５度の位置及び２２５度の位置にあり、右画像の回転角度検出用電子透かしの非ゼロ周波数成分の重心は、１３５度の位置及び３１５度の位置にあり、左画像には、第１データの電子透かしが挿入されており、右画像には第２データの電子透かしが挿入されている。
【００３１】
従って、回転角度に依存しないで、第１データは、重心が４５度及び２２５度の位置にある回転角度検出用電子透かしに対応し、第２データは、重心が１３５度及び２２５度の位置にある回転角度検出用電子透かしの非ゼロ周波数成分に対応する。従って、カメラで撮影した画像を左右に分割し、左画像及び右画像別々に処理を行い、各画像について、データ及び回転角度検出用電子透かしを検出して、第１データ又は第２データの値及び第１データであるのか第２データであるのかの別を識別することにより、最終的な全体のデータを得ることができる。具体的には、例えば、２０個の位置のうちの１０個の位置に非ゼロ周波数成分を配置することにより各分割画像のデータを挿入しているのであれば、第１データに２５２を乗じた値に第２データを加算した値が最終的な全体のデータである。
【００３２】
図２は、画像を左右に２分割する場合の電子透かし挿入及び電子透かし検出の様子を示すものであったが、図３は、画像を上下左右に４分割する場合の電子透かし挿入及び電子透かし検出の様子を示すものである。
【００３３】
図３を参照すると、印刷時には、右上画像の回転角度検出用電子透かしの非ゼロ周波数成分の重心は、４５度の位置及び２２５度の位置にあり、左上画像の回転角度検出用電子透かしの非ゼロ周波数成分の重心は、９０度の位置及び２７０度の位置にあり、左下画像の回転角度検出用電子透かしの非ゼロ周波数成分の重心は、１３５度の位置及び３１５度の位置にあり、右下画像の回転角度検出用電子透かしの非ゼロ周波数成分の重心は、０度の位置及び１８０度の位置にあり、右上画像には、第１データの電子透かしが挿入され、左上画像には第２データの電子透かしが挿入され、左下画像には、第３データの電子透かしが挿入され、右下画像には第４データの電子透かしが挿入される。
【００３４】
図３に示す例では、例えば、２０個の位置のうちの１０個の位置に非ゼロ周波数成分を配置することにより各分割画像にデータを挿入しているのであれば、２５２＾４＝４，０３２，７５８，０１６通りのデータを表すことができる。
【００３５】
カメラにより０度付近の回転角度で撮影された時には、印刷時と同様に、右上画像の回転角度検出用電子透かしの非ゼロ周波数成分の重心は、４５度の位置及び２２５度の位置にあり、左上画像の回転角度検出用電子透かしの非ゼロ周波数成分の重心は、９０度の位置及び２７０度の位置にあり、左下画像の回転角度検出用電子透かしの非ゼロ周波数成分の重心は、１３５度の位置及び３１５度の位置にあり、右下画像の回転角度検出用電子透かしの非ゼロ周波数成分の重心は、０度の位置及び１８０度の位置にあり、右上画像には、第１データの電子透かしが挿入されており、左上画像には第２データの電子透かしが挿入されており、左下画像には、第３データの電子透かしが挿入されており、右下画像には第４データの電子透かしが挿入されている。
【００３６】
一方、カメラにより９０度付近の回転角度で撮影された時には、印刷時と異なり、右上画像の回転角度検出用電子透かしの非ゼロ周波数成分の重心は、９０度の位置及び２７０度の位置にあり、左上画像の回転角度検出用電子透かしの非ゼロ周波数成分の重心は、１３５度の位置及び３１５度の位置にあり、左下画像の回転角度検出用電子透かしの非ゼロ周波数成分の重心は、０度の位置及び１８０度の位置にあり、右下画像の回転角度検出用電子透かしの非ゼロ周波数成分の重心は、４５度の位置及び２２５度の位置にあり、右上画像には、第４データの電子透かしが挿入されており、左上画像には第１データの電子透かしが挿入されており、左下画像には、第２データの電子透かしが挿入されており、右下画像には第３データの電子透かしが挿入されている。
【００３７】
また、カメラにより１８０度付近の回転角度で撮影された時には、印刷時と異なり、右上画像の回転角度検出用電子透かしの非ゼロ周波数成分の重心は、１３５度の位置及び３１５度の位置にあり、左上画像の回転角度検出用電子透かしの非ゼロ周波数成分の重心は、０度の位置及び１８０度の位置にあり、左下画像の回転角度検出用電子透かしの非ゼロ周波数成分の重心は、４５度の位置及び２２５度の位置にあり、右下画像の回転角度検出用電子透かしの非ゼロ周波数成分の重心は、９０度の位置及び２７０度の位置にあり、右上画像には、第３データの電子透かしが挿入されており、左上画像には第４データの電子透かしが挿入されており、左下画像には、第１データの電子透かしが挿入されており、右下画像には第２データの電子透かしが挿入されている。
【００３８】
更に、カメラにより２７０度付近の回転角度で撮影された時には、印刷時と異なり、右上画像の回転角度検出用電子透かしの非ゼロ周波数成分の重心は、０度の位置及び１８０度の位置にあり、左上画像の回転角度検出用電子透かしの非ゼロ周波数成分の重心は、４５度の位置及び２２５度の位置にあり、左下画像の回転角度検出用電子透かしの非ゼロ周波数成分の重心は、９０度の位置及び２７０度の位置にあり、右下画像の回転角度検出用電子透かしの非ゼロ周波数成分の重心は、１３５度の位置及び３１５度の位置にあり、右上画像には、第２データの電子透かしが挿入されており、左上画像には第３データの電子透かしが挿入されており、左下画像には、第４データの電子透かしが挿入されており、右下画像には第１データの電子透かしが挿入されている。
【００３９】
図２に示す左右２分割の場合には、他の分割画像の回転角度検出用非ゼロ周波数成分の重心を考慮せず、当該分割画像の回転角度検出用非ゼロ周波数成分の重心のみを基に、当該分割画像のデータは第１データであるのか第２データであるのかの別を判断することができた。これに対し、図３に示す上下左右４分割の場合にはそのようは判断はできない。すなわち、例えば、回転角度が０度である場合には、４５度及び２２５度の回転角度検出用非ゼロ周波数成分の重心に対応するデータは第１データであるのに対し、回転角度が９０度である場合には、４５度及び２２５度の回転角度検出用非ゼロ周波数成分の重心に対応するデータは第３データであり、回転角度が１８０度である場合には、４５度及び２２５度の回転角度検出用非ゼロ周波数成分の重心に対応するデータは第１データであり、回転角度が２７０度である場合には、４５度及び２２５度の回転角度検出用非ゼロ周波数成分の重心に対応するデータは第３データである。しかし、当該分割画像がどこの分割画像であるのかということ、及び、当該分割画像からどの重心を有する回転角度検出用電子透かしが検出されたのかを検出することにより、当該分割画像には何番目のデータが挿入されているのかを検出することができる。具体的には、右上分割画像を例に取ると、当該分割画像が右上分割画像であるということ、及び当該分割画像に４５度及び２２５度の重心を有する観点角度検出用電子透かしが検出されたことより、当該分割画像には第１データが挿入されていることを検出することができ、当該分割画像が右上分割画像であるということ、及び当該分割画像に９０度及び２７０度の重心を有する観点角度検出用電子透かしが検出されたことより、当該分割画像には第４データが挿入されていることを検出することができ、当該分割画像が右上分割画像であるということ、及び当該分割画像に１３５度及び３１５度の重心を有する観点角度検出用電子透かしが検出されたことより、当該分割画像には第３データが挿入されていることを検出することができ、当該分割画像が右上分割画像であるということ、及び当該分割画像に０度及び１８０度の重心を有する観点角度検出用電子透かしが検出されたことより、当該分割画像には第２データが挿入されていることを検出することができる。
【００４０】
図３に示す方法では、画像の回転角度が０度付近であるのか、９０度付近であるのか、１８０度付近であるのか、又は、２７０度付近であるのかを検出することができるのに対し、図２に示す方法では、画像の回転角度が０度付近であるのか、又は、１８０度付近であるのかを検出することができるが、画像の回転角度が９０度付近であること及び画像の回転角度が２７０度付近であることを検出することができない。しかし、カメラを横にした状態で画像を撮影する場合のみならず、カメラを縦にした状態で画像を撮影する場合もあるので、画像の回転角度が９０度付近であること及び画像の回転角度が２７０度付近であることも検出できることが好ましい。図４は、このようなことを可能とする方法を説明するための図である。
【００４１】
この方法は、図２に示す場合と同様に、画像を左右に２分割して、電子透かしを挿入することを前提としながらも、画像を上下左右に４分割して電子透かしを検出する。
【００４２】
カメラにより０度付近の回転角度で撮影された時には、右上画像の回転角度検出用電子透かしの非ゼロ周波数成分の重心は、４５度の位置及び２２５度の位置にあり、左上画像の回転角度検出用電子透かしの非ゼロ周波数成分の重心は、１３５度の位置及び３１５度の位置にあり、左下画像の回転角度検出用電子透かしの非ゼロ周波数成分の重心は、１３５度の位置及び３１５度の位置にあり、右下画像の回転角度検出用電子透かしの非ゼロ周波数成分の重心は、４５度の位置及び２２５度の位置にあり、右上画像には、第１データの電子透かしが挿入されており、左上画像には第２データの電子透かしが挿入されており、左下画像には、第２データの電子透かしが挿入されており、右下画像には第１データの電子透かしが挿入されている。
【００４３】
一方、カメラにより９０度付近の回転角度で撮影された時には、右上画像の回転角度検出用電子透かしの非ゼロ周波数成分の重心は、１３５度の位置及び３１５度の位置にあり、左上画像の回転角度検出用電子透かしの非ゼロ周波数成分の重心は、１３５度の位置及び３１５度の位置にあり、左下画像の回転角度検出用電子透かしの非ゼロ周波数成分の重心は、４５度の位置及び２２５度の位置にあり、右下画像の回転角度検出用電子透かしの非ゼロ周波数成分の重心は、４５度の位置及び２２５度の位置にあり、右上画像には、第１データの電子透かしが挿入されており、左上画像には第１データの電子透かしが挿入されており、左下画像には、第２データの電子透かしが挿入されており、右下画像には第２データの電子透かしが挿入されている。
【００４４】
また、カメラにより１８０度付近の回転角度で撮影された時には、右上画像の回転角度検出用電子透かしの非ゼロ周波数成分の重心は、１３５度の位置及び３１５度の位置にあり、左上画像の回転角度検出用電子透かしの非ゼロ周波数成分の重心は、４５度の位置及び２２５度の位置にあり、左下画像の回転角度検出用電子透かしの非ゼロ周波数成分の重心は、４５度の位置及び２２５度の位置にあり、右下画像の回転角度検出用電子透かしの非ゼロ周波数成分の重心は、１３５度の位置及び３１５度の位置にあり、右上画像には、第２データの電子透かしが挿入されており、左上画像には第１データの電子透かしが挿入されており、左下画像には、第１データの電子透かしが挿入されており、右下画像には第２データの電子透かしが挿入されている。
【００４５】
更に、カメラにより２７０度付近の回転角度で撮影された時には、右上画像の回転角度検出用電子透かしの非ゼロ周波数成分の重心は、４５度の位置及び２２５度の位置にあり、左上画像の回転角度検出用電子透かしの非ゼロ周波数成分の重心は、４５度の位置及び２２５度の位置にあり、左下画像の回転角度検出用電子透かしの非ゼロ周波数成分の重心は、１３５度の位置及び３１５度の位置にあり、右下画像の回転角度検出用電子透かしの非ゼロ周波数成分の重心は、１３５度の位置及び３１５度の位置にあり、右上画像には、第２データの電子透かしが挿入されており、左上画像には第２データの電子透かしが挿入されており、左下画像には、第１データの電子透かしが挿入されており、右下画像には第１データの電子透かしが挿入されている。
【００４６】
図４に示す方法の場合には、他の分割画像の回転角度検出用非ゼロ周波数成分の重心を考慮せず、当該分割画像の回転角度検出用非ゼロ周波数成分の重心のみを基に、当該分割画像のデータは第１データであるのか第２データであるのかの別を判断することはできない。すなわち、例えば、回転角度が０度である場合には、第１象限において４５度及び２２５度の回転角度検出用非ゼロ周波数成分の重心に第１データが対応し、第４象限において４５度及び２２５度の回転角度検出用非ゼロ周波数成分の重心に第１データが対応するのに対し、回転角度が９０度である場合には、第１象限において１３５度及び３１５度の回転角度検出用非ゼロ周波数成分の重心に第１データが対応し、第２象限において１３５度及び３１５度の回転角度検出用非ゼロ周波数成分の重心に第１データが対応し、回転角度が１８０度である場合には、第３象限において４５度及び２２５度の回転角度検出用非ゼロ周波数成分の重心に第１データが対応し、第４象限において４５度及び２２５度の回転角度検出用非ゼロ周波数成分の重心に第１データが対応し、回転角度が２７０度である場合には、第３象限において１３５度及び３１５度の回転角度検出用非ゼロ周波数成分の重心に第１データが対応し、第４象限において１３５度及び３１５度の回転角度検出用非ゼロ周波数成分の重心に第１データが対応する。
【００４７】
従って、図４に示す挿入時左右２分割検出時上下左右４分割の場合には、４つの分割画像における回転角度検出用非ゼロ周波数成分の重心を総合的にみることにより、画像の回転角度が０度付近であるのか、９０度付近であるのか、１８０度付近であるのか、又は、２７０度付近であるのかを判断する。
【００４８】
この様にして、各分割画像における回転角度検出用非ゼロ周波数成分の重心を総合的にみることにより画像の回転角度を検出し、その回転角度を基に、各分割画像から検出されたデータが第何番目のデータであるのかを判断し、各分割画像から検出されたデータとその番号を基に最終的な全体のデータを求めることができる。具体的には、例えば、２０個の位置のうちの１０個の位置に非ゼロ周波数成分を配置することにより各分割画像のデータを挿入しているのであれば、第１データに２５２＾２を乗じた値に第２データを加算することにより得た値が最終的な全体のデータである。
【００４９】
上記の何れの方法においても、カメラで撮影した画像を分割してから、各分割画像において回転角度検出用の電子透かしを検出した。これに対し、図５を参照して説明する方法は、カメラで撮影した画像を分割する前に、画像全体に挿入した回転角度検出用の電子透かしを検出し、これを基に画像の回転角度を検出し、この回転角度により画像の回転補正を行った後に、画像を分割し、各分割画像においてデータの電子透かしを検出するものである。
【００５０】
図５を参照すると、符号６０７で示すように、画像全体に画像の回転角度を検出するための電子透かしを挿入する。また、例えば、データを２分割して、２５２＾２＝６３５０４通りのデータを表す場合には、符号６０８で示すように、画像を２分割し、第ｉ割画像に第ｉデータを表す電子透かしを挿入する（ｉ＝１、２）。同様に、例えば、データを４分割して、２５２＾４＝４０３２７５８０１６通りのデータを表す場合には、符号６０９で示すように、画像を４分割し、第ｉ割画像に第ｉデータを表す電子透かしを挿入する（ｉ＝１〜４）。例えば、データを６分割して、２５２＾６＝２．５６１×１０＾１４通りのデータを表す場合には、符号６０９で示すように、画像を６分割し、第ｉ割画像に第ｉデータを表す電子透かしを挿入する（ｉ＝１〜６）。例えば、データを９分割して、２５２＾９＝４．０９８×１０＾２１通りのデータを表す場合には、符号６０９で示すように、画像を９分割し、第ｉ割画像に第ｉデータを表す電子透かしを挿入する（ｉ＝１〜９）。
【００５１】
次に、上述した電子透かし挿入方法及び電子透かし検出方法を行う電子透かし挿入装置及び電子透かし検出装置について説明する。
【００５２】
図６は、本発明の実施形態による電子透かし挿入装置の構成を示すブロック図である。
【００５３】
図６を参照すると、本発明の実施形態による電子透かし挿入装置は、回転角度検出用パターン生成部１０１、データ保持部１０２、データパターン生成部１０３、合成部１０４、逆離散フーリエ変換部１０５、二値化部１０６及び加算部（減算部）１０７を備える。なお、これらの部分は、ハードウェアによって実現してもよいが、コンピュータをこれらの部分として機能させるためのプログラムをコンピュータが読み込んで実行することにより実現してもよい。
【００５４】
次に、各分割画像毎に異なる回転角度検出用の電子透かし及び部分データを各分割画像毎に挿入する場合の電子透かし挿入装置により行われる電子透かし挿入方法について、図６及び７を参照して説明する。
【００５５】
この電子透かし挿入方法においては、第１の分割画像から第Ｎの分割画像についてステップＳ１５２乃至Ｓ１５６を繰り返す（ステップＳ１５１）。各繰返しにおいては以下のことを行う。
【００５６】
まず、回転角度検出用パターン生成部１０１は、第ｊの回転角度検出用パターンを生成する（ステップＳ１５１）。このパターンは、非ゼロ周波数成分より構成される。
【００５７】
次に、データパターン生成部１０３は、データ保持部１０２に保持されている全体データを読み込み、現在扱っているｊの値に応じて第ｊのデータを求め、第ｊのデータに応じて第ｊのデータパターンを生成する（ステップＳ１５２）。このパターンも非ゼロ周波数成分より構成される。第ｊのデータパターンが図１に示すようなものである場合には、第ｊのデータＤｊは、全体データをＤとすると、
Ｄｊ＝ｍｏｄ（Ｄ，２５２＾（Ｎ＋１−ｊ））／２５２＾（Ｎ−ｊ）
により表される。
【００５８】
次に、合成部１０４は、ステップＳ１５１で生成した第ｊの回転角度検出用パターンとステップＳ１５２で生成した第ｊのデータパターンを足し合わせることにより合成パターンを生成する（ステップＳ１５３）。
【００５９】
次に、逆離散フーリエ変換部１０５は、ステップＳ１５３で合成したパターンに対し、逆離散フーリエ変換をかけることにより、空間領域パターンを得る（ステップＳ１５４）。
【００６０】
次に、二値化部１０６は、ステップＳ１５４で得た空間領域パターンを二値化する（ステップＳ１５５）。すなわち、各画素について、ステップＳ１５４で得た空間領域パターンの画素値が０以上であれば、所定の値（例えば、１０）とし、ステップＳ１５４で得た空間パターンの画素値が０未満であれば、所定の値（例えば−１０）とする。
【００６１】
次に、加算部（減算部）１０７は、ステップＳ１０６で二値化された空間領域パターンを原画像の黄成分の第ｊ領域の全部又は一部のブロックに加算又は減算する（ステップＳ１５６）。なお、二値化部１０６が二値化を行うのは、電子透かし検出装置が電子透かしを検出することを容易にするためであり、省略してもよい。この場合には、逆離散フーリエ変換部１０５が出力する空間パターンを二値化せずに直接原画像の黄成分の第ｊ領域に加算又は減算する。
【００６２】
その後、上記の方法により電子透かしが挿入された画像は印刷物に印刷される。
【００６３】
次に、分割画像間で共通な（画像全体で共通な）回転角度検出用の電子透かしを画像全体に挿入する場合の電子透かし挿入装置により行われる電子透かし挿入方法について、図６及び８を参照して説明する。
【００６４】
まず、まず、回転角度検出用パターン生成部１０１は、分割画像間で共通な回転角度検出用パターンを生成する（ステップＳ１５１）。このパターンは、非ゼロ周波数成分より構成される。
【００６５】
次に、第１の分割画像から第Ｎの分割画像についてステップＳ５５３乃至Ｓ５５７を繰り返す（ステップＳ５５２）。各繰返しにおいては以下のことを行う。
【００６６】
まず、データパターン生成部１０３は、データ保持部１０２に保持されている全体データを読み込み、現在扱っているｊの値に応じて第ｊのデータを求め、第ｊのデータに応じて第ｊのデータパターンを生成する（ステップＳ５５３）。このパターンも非ゼロ周波数成分より構成される。第ｊのデータパターンが図１に示すようなものである場合には、第ｊのデータＤｊは、全体データをＤとすると、
Ｄｊ＝ｍｏｄ（Ｄ，２５２＾（Ｎ＋１−ｊ））／２５２＾（Ｎ−ｊ）
により表される。
【００６７】
次に、合成部１０４は、ステップＳ５５１で生成した共通回転角度検出用パターンとステップＳ５５３で生成した第ｊのデータパターンを足し合わせることにより合成パターンを生成する（ステップＳ５５４）。
【００６８】
次に、逆離散フーリエ変換部１０５は、ステップＳ５５４で合成したパターンに対し、逆離散フーリエ変換をかけることにより、空間領域パターンを得る（ステップＳ５５５）。
【００６９】
次に、二値化部１０６は、ステップＳ５５５で得た空間領域パターンを二値化する（ステップＳ５５６）。すなわち、各画素について、ステップＳ５５５で得た空間領域パターンの画素値が０以上であれば、所定の値（例えば、１０）とし、ステップＳ５５５で得た空間パターンの画素値が０未満であれば、所定の値（例えば−１０）とする。
【００７０】
次に、加算部（減算部）１０７は、ステップＳ５５６で二値化された空間領域パターンを原画像の黄成分の第ｊ領域の全部又は一部のブロックに加算又は減算する（ステップＳ５５７）。なお、二値化部１０６が二値化を行うのは、電子透かし検出装置が電子透かしを検出することを容易にするためであり、省略してもよい。この場合には、逆離散フーリエ変換部１０５が出力する空間パターンを二値化せずに直接原画像の黄成分の第ｊ領域に加算又は減算する。
【００７１】
その後、上記の方法により電子透かしが挿入された画像は印刷物に印刷される。
【００７２】
図９は、図２を参照して説明した検出方法（データ２分割・回転角度検出用電子透かし２分割・検出時２分割）を行うための電子透かし検出装置の構成を示すブロック図である。
【００７３】
図９を参照すると、この電子透かし検出装置は、空間フィルタ２０１、ブロック化部２０２、離散フーリエ変換部２０３、振幅計算部２０４、振幅積算部２０５、回転角度検出・領域番号識別部２０６、部分データ検出部２０７及びデータ検出部２０８を備える。なお、これらの部分はハードウェアによって実現してもよいが、コンピュータをこれらの部分として機能させるためのプログラムをコンピュータが読み込んで実行することによって実現してもよい。
【００７４】
次に、図９に示す電子透かし検出装置により行われる電子透かし検出方法について、図９及び１０を参照して説明する。
【００７５】
この電子透かし検出方法においては、第１の分割画像から第２の分割画像についてステップＳ２５２乃至Ｓ２５９を繰り返す（ステップＳ２５１）。各繰返しにおいては以下のことを行う。
【００７６】
まず、カメラにより撮影した画像の青成分の現在の分割画像（ｋ番目の分割画像）に対して空間フィルタ２０１を適用する（ステップＳ２５２）。この空間フィルタは、例えば、図１１に示すようなタップ係数を有する。タップ係数から明らかなように、この空間フィルタは帯域通過の特性を有する。従って、画像にこの空間フィルタ２０１を適用することにより、原画像に起因する低周波数成分や高周波数成分を抑圧することができる。
【００７７】
次に、ブロック化部２０２は、空間フィルタが適用された分割画像を離散フーリエ変換を行うためのサイズの１又は２以上のブロックに分割する（ステップＳ２５３）。
【００７８】
次に、離散フーリエ変換部２０３は、各ブロックの画像に対して離散フーリエ変換を行い、周波数成分を得る（ステップＳ２５４）。
【００７９】
次に、振幅計算部２０４は、各ブロックの各周波数成分の振幅を求める（ステップＳ２５５）。
【００８０】
次に、振幅積算部２０５は、分割画像に含まれるブロックの間で同一の周波数の周波数成分を積算することにより、分割画像についての各周波数成分の積算振幅を得る（ステップＳ２５６）。
【００８１】
次に、回転角度検出・領域番号識別部２０６は、各周波数成分の積算振幅を入力し、図１の内円６０２に分布する非ゼロ周波数成分の分布の様子をみることにより、現在の分割画像が、第１分割画像であるのか第２分割画像であるのかを判断する（ステップＳ２５７）。すなわち、図２を参照すると、重心が４５度の位置と２２５度の位置にある場合には、現在の分割画像が第１分割画像であると判断し、重心が１３５度の位置と３１５度の位置にある場合には、現在の分割画像が第２分割画像であると判断する。
【００８２】
なお、回転角度を検出するための非ゼロ周波数成分の２つの重心を求めることは、例えば、任意の１つの非ゼロ周波数成分を選択し、この非ゼロ周波数成分からの距離が相対的に短い（半数−１）の数の非ゼロ周波数成分を含めて、半数の非ゼロ周波数成分を抽出し、その半数の非ゼロ周波数成分が第１グループに属するとし、残りの半数が第２グループに属するとし、各グループの重心を求めることにより可能となる。
【００８３】
次に、回転角度検出・領域番号識別部２０６は、各周波数成分の積算振幅を入力し、図１の内円６０２に分布する非ゼロ周波数成分の分布の様子をみて、更に、ステップＳ２５７で判断した現在の分割画像の番号をみることにより、画像の回転角度を検出する（ステップＳ２５８）。すなわち、現在の分割画像が左画像であり且つ第２分割画像であると判断されたのであれば、画像が回転していない場合の重心が１３５度の位置にある非ゼロ周波数成分を基準として、重心が１３５度付近の位置にある実際の非ゼロ周波数成分の位置を基に、０度付近の回転角度を検出し、現在の分割画像が右画像であり且つ第１分割画像であると判断されたのであれば、画像が回転していない場合の重心が４５度の位置にある非ゼロ周波数成分を基準として、重心が４５度付近の位置にある実際の非ゼロ周波数成分の位置を基に、０度付近の回転角度を検出し、現在の分割画像が左画像であり且つ第１分割画像であると判断されたのであれば、画像が回転していない場合の重心が４５度の位置にある非ゼロ周波数成分を基準として、重心が４５度付近の位置にある実際の非ゼロ周波数成分の位置を基に、１８０度付近の回転角度を検出し、現在の分割画像が右画像であり且つ第２分割画像であると判断されたのであれば、画像が回転していない場合の重心が１３５度の位置にある非ゼロ周波数成分を基準として、重心が１３５度付近の位置にある実際の非ゼロ周波数成分の位置を基に１８０度付近の回転角度を検出する。
【００８４】
次に、部分データ検出部２０６は、各周波数成分の積算振幅を入力し、図１の外円６０１に分布する非ゼロ周波数成分の分布の様子をみることにより、データを検出し、これを第ｊデータとする（ステップＳ２５９）。ここで、番号ｊは、ステップＳ２５７で検出された分割画像の番号のことである。
【００８５】
ステップＳ２５１の繰返しが終了したならば、ステップＳ２５９をＮ回行うことによって得た第１データから第Ｎデータを基に全体のデータを検出する（ステップＳ２６０）。
【００８６】
ステップＳ２５７及びＳ２５８では、内円６０２に分布する回転角データを検出するための各非ゼロ周波数成分の周波数領域における極座標でみた角度を検出する必要があるが、そのためには、以下の方法を用いる。
【００８７】
すなわち、図１２を参照すると、まず、パラメータである半径をゼロに初期化し（ステップＳ９０１）、現在の半径の円周上にＮ１個（図１の例では６個）の非ゼロ周波数成分があるか否かについて判断する（ステップＳ９０２）。
【００８８】
ステップＳ９０２においては、以下のようなことをする。すなわち、図１の例では、円周を１６個に等分割する角度毎に現在の半径の円上にウィンドウを設け、全周波数成分をそのウィンドウを通した後にしきい値を超える振幅の非ゼロ周波数成分がＮ１個検出できるか否かをウィンドウの角度をずらしながらみていく。検出できた場合には、そのＮ１個の非ゼロ周波数成分が検出しようとしている非ゼロ周波数成分である。
【００８９】
無ければ（ステップＳ９０２でＮＯ）、半径を周波数領域における単位長さΔＲだけ増加させ（ステップＳ９０３）、ステップＳ９０２に戻る。
【００９０】
あれば（ステップＳ９０２でＹＥＳ）、現在の半径の円上に分布する回転角度を検出するための各非ゼロ周波数成分の周波数領域における極座標でみた角度を検出する（ステップＳ９０４）。
【００９１】
図１２には示していないが、半径が最大半径になってもＮ１個の非ゼロ周波数成分が検出できない場合には、ウィンドウの半径方向の幅を広げてステップＳ９０２及びステップＳ９０３を繰り返す。それでもどの半径でもＮ１個の非ゼロ周波数成分が検出できない場合には、ウィンドウの半径方向の幅を戻した上で上記しきい値を下げて、ステップＳ９０２及びＳ９０３を繰り返す。それでもどの半径でもＮ１個の非ゼロ周波数成分が検出できない場合には、再度ウィンドウの半径方向の幅を広げてステップＳ９０２及びＳ９０３を繰り返す。以下同様なことを繰り返す。
【００９２】
また、ステップＳ２５９では、外円６０１に分布するデータを表す各非ゼロ周波数成分の周波数領域における極座標でみた角度を検出する必要があるが、そのためには以下の方法を用いる。
【００９３】
すなわち、図１３を参照すると、まず、ステップ２５８で検出した回転角度により全ての周波数成分の振幅の位置を補正する（ステップＳ９２１）。
【００９４】
次に、パラメータである半径をステップ２５７で検出した内円の半径に初期化し（ステップＳ９２２）、現在の半径の円周上にＮ２個（図１の例では１０個）の非ゼロ周波数成分があるか否かについて判断する（ステップＳ９２３）。
【００９５】
ステップＳ９２３においては、以下のようなことをする。すなわち、図１に示す例であれば、現在の半径の円上のデータを表す非ゼロ周波数成分が存在する可能性がある２０の点に所定の広さを有するウィンドウを設け、全周波数成分をそのウィンドウを通した後にしきい値を超える振幅の非ゼロ周波数成分が１０個検出できるか否かをみていく。検出できた場合には、その１０個の非ゼロ周波数成分が検出しようとしている非ゼロ周波数成分である。
【００９６】
無ければ（ステップＳ９２３でＮＯ）、半径を周波数領域における単位長さΔＲだけ増加させ（ステップＳ９２４）、ステップＳ９２３に戻る。
【００９７】
あれば（ステップＳ９２３でＹＥＳ）、以下の方法により第ｊデータを検出する（ステップＳ９２５）。すなわち、現在の半径の円上に分布するデータを表す各非ゼロ周波数成分の周波数領域における極座標でみた角度を検出する。そして、非ゼロ周波数成分の分布により非ゼロ周波数成分が表すデータの値を求める。そして、ステップＳ２５７で検出した番号ｊに従って、そのデータを第ｊデータとする。
【００９８】
図１３には示していないが、半径が最大半径になっても１０個の非ゼロ周波数成分が検出できない場合には、ウィンドウの半径方向の幅を広げてステップＳ９２３及びステップＳ９２４を繰り返す。それでもどの半径でも１０個の非ゼロ周波数成分が検出できない場合には、ウィンドウの半径方向の幅を戻した上で上記しきい値を下げて、ステップＳ９２３及びＳ９２４を繰り返す。それでもどの半径でも１０個の非ゼロ周波数成分が検出できない場合には、再度ウィンドウの半径方向の幅を広げてステップＳ９２３及びＳ９２４を繰り返す。以下同様なことを繰り返す。
【００９９】
図１４は、図３を参照して説明した検出方法（データ４分割・回転角度検出用電子透かし４分割・検出時４分割）を行うための電子透かし検出装置の構成を示すブロック図である。
【０１００】
図１４を参照すると、この電子透かし検出装置は、空間フィルタ２０１、ブロック化部２０２、離散フーリエ変換部２０３、振幅計算部２０４、振幅積算部２０５、回転角度検出・パターン番号識別部３０１、部分データ番号検出部３０２、部分データ検出部３０３及びデータ検出部３０４を備える。なお、これらの部分はハードウェアによって実現してもよいが、コンピュータをこれらの部分として機能させるためのプログラムをコンピュータが読み込んで実行することによって実現してもよい。
【０１０１】
次に、図１４に示す電子透かし検出装置により行われる電子透かし検出方法について、図１４及び１５を参照して説明する。
【０１０２】
この電子透かし検出方法においては、まず、第１の分割画像から第４の分割画像についてステップＳ３５２乃至Ｓ３６０を繰り返す（ステップＳ３５１）。各繰返しにおいては以下のことを行う。
【０１０３】
まず、カメラにより撮影した画像の青成分の現在の分割画像（ｋ番目の分割画像）に対して空間フィルタ２０１を適用する（ステップＳ３５２）。この空間フィルタは、例えば、図１１に示すようなタップ係数を有する。タップ係数から明らかなように、この空間フィルタは帯域通過の特性を有する。従って、画像にこの空間フィルタ２０１を適用することにより、原画像に起因する低周波数成分や高周波数成分を抑圧することができる。
【０１０４】
次に、ブロック化部２０２は、空間フィルタが適用された分割画像を離散フーリエ変換を行うためのサイズの１又は２以上のブロックに分割する（ステップＳ３５３）。
【０１０５】
次に、離散フーリエ変換部２０３は、各ブロックの画像に対して離散フーリエ変換を行い、周波数成分を得る（ステップＳ３５４）。
【０１０６】
次に、振幅計算部２０４は、各ブロックの各周波数成分の振幅を求める（ステップＳ３５５）。
【０１０７】
次に、振幅積算部２０５は、分割画像に含まれるブロックの間で同一の周波数の周波数成分を積算することにより、分割画像についての各周波数成分の積算振幅を得る（ステップＳ３５６）。
【０１０８】
次に、回転角度検出・パターン番号識別部３０１は、各周波数成分の積算振幅を入力し、図１の内円６０２に分布する非ゼロ周波数成分の分布の様子をみることにより、パターン番号ｐを識別する（ステップＳ３５７）。すなわち、図３を参照すると、重心が０度の位置の近傍と１８０度の位置の近傍にある場合には、現在の分割画像に含まれているパターン番号は１であると判断し、重心が４５度の位置の近傍と２２５度の位置の近傍にある場合には、現在の分割画像に含まれているパターン番号は２であると判断し、重心が９０度の位置の近傍と２７０度の位置の近傍にある場合には、現在の分割画像に含まれているパターン番号は３であると判断し、重心が１３５度の位置の近傍と３１５度の位置の近傍にある場合には、現在の分割画像に含まれているパターン番号は４であると判断する。
【０１０９】
次に、部分データ番号検出部３０２は、現在の繰返し番号ｋと一致する分割画像の番号ｋと第ｋの分割画像を処理した時に検出されたパターン番号ｐとの対応関係全体を基に、現在処理している分割画像に何番目の部分データが対応するのか（すなわち、撮影分割画像番号ｋと部分データ番号ｊとの対応関係）を検出する（ステップＳ３５８）。すなわち、図３を参照すると、撮影分割画像の番号ｋを象限の番号であるとすると、ｋ＝１の繰返しにおいては、ｋとｐとの対応関係が、（１、２）である場合には、ｋとｊとの対応関係が（１、１）であることを検出し、ｋとｐとの対応関係が、（１、３）である場合には、ｋとｊとの対応関係が（１、４）であることを検出し、ｋとｐとの対応関係が、（１、４）である場合には、ｋとｊとの対応関係が（１、３）であることを検出し、ｋとｐとの対応関係が、（１、１）である場合には、ｋとｊとの対応関係が（１、２）であることを検出する。ｋ＝２の繰返しにおいては、ｋとｐとの対応関係が、（２、３）である場合には、ｋとｊとの対応関係が（２、２）であることを検出し、ｋとｐとの対応関係が、（２、４）である場合には、ｋとｊとの対応関係が（２、１）であることを検出し、ｋとｐとの対応関係が、（２、１）である場合には、ｋとｊとの対応関係が（２、４）であることを検出し、ｋとｐとの対応関係が、（２、２）である場合には、ｋとｊとの対応関係が（２、３）であることを検出する。ｋ＝３の繰返しにおいては、ｋとｐとの対応関係が、（３、４）である場合には、ｋとｊとの対応関係が（３、３）であることを検出し、ｋとｐとの対応関係が、（３、１）である場合には、ｋとｊとの対応関係が（３、２）であることを検出し、ｋとｐとの対応関係が、（３、２）である場合には、ｋとｊとの対応関係が（３、１）であることを検出し、ｋとｐとの対応関係が、（３、３）である場合には、ｋとｊとの対応関係が（３、４）であることを検出する。ｋ＝４の繰返しにおいては、ｋとｐとの対応関係が、（４、１）である場合には、ｋとｊとの対応関係が（４、４）であることを検出し、ｋとｐとの対応関係が、（４、２）である場合には、ｋとｊとの対応関係が（４、３）であることを検出し、ｋとｐとの対応関係が、（４、３）である場合には、ｋとｊとの対応関係が（４、２）であることを検出し、ｋとｐとの対応関係が、（４、４）である場合には、ｋとｊとの対応関係が（４、１）であることを検出する。
【０１１０】
なお、データ２分割に対応した図１０を参照して説明した方法では、回転角度検出用の非ゼロ周波数成分の重心のみから直ちに部分データの番号を検出することができた。これに対し、データ４分割に対応した当該方法では、回転角度検出用の非ゼロ周波数成分の重心のみからでは部分データの番号を検出することはできず、回転角度検出用の非ゼロ周波数成分の重心と現在処理している分割画像の象限の両者との対応関係から部分データの番号を検出することができる。
【０１１１】
ステップＳ３５８の次に、回転角度検出用の電子透かしを基に、画像の回転角度を検出する（ステップＳ３５９）。すなわち、図３を参照すると、現在の繰返しのステップＳ３５７において回転角度検出用電子透かしのパターンがパターン１乃至４のうちのどちらであるのかを検出したので、画像が回転していない場合のその番号のパターンの角度を基準として区分内回転角度を検出する。すなわち、画像が回転していない場合のその番号のパターンの角度と検出されたパターンの角度との差を小回転角度として検出する。更に、現在の分割画像の番号ｋとパターン番号ｐとの対応関係により、画像の回転角度が０度付近であるのか、９０度付近であるのか、１８０度付近であるのか、又は、２７０度付近であるのかを検出し、それに応じて、０度、９０度、１８０度又は２７０度を大回転角度として検出する。そして、小回転角度と大回転角度を合わせたものを回転角度とする。
【０１１２】
次に、部分データ検出部３０３は、当該分割画像から部分データを検出する（ステップＳ３６０）。検出方法は、図１３を参照して説明した方法と同様である。この際、ステップＳ３５９で検出した回転角度をステップＳ９２１で補正角度として利用する。また、ステップＳ３５８で求めた番号ｊを部分データの番号として利用する。
【０１１３】
ステップＳ３５１の繰返しが終了したならば、データ検出部３０４は、全ての繰返しにおけるステップＳ３６０で求めた第１部分データ乃至第４部分データを基にデータを検出する（ステップＳ３６２）。
【０１１４】
なお、空間フィルタ適用（ステップＳ３５２）、ブロック化（ステップＳ３５３）、離散フーリエ変換（ステップＳ３５４）、振幅計算（ステップＳ３５５）は、ステップＳ３５１の繰返しの内部において行う代わりにステップＳ３５１の繰返しに入る前に画像全体に対して行ってもよい。
【０１１５】
また、図１のパターンに対して、図１４に示した電子透かし検出装置及び図１５を参照して説明した電子透かし検出方法を適用することもできる。
【０１１６】
図１６は、図４を参照して説明した検出方法（データ２分割・回転角度検出用電子透かし２分割・検出時４分割）を行うための電子透かし検出装置の構成を示すブロック図である。
【０１１７】
図１６を参照すると、この電子透かし検出装置は、空間フィルタ２０１、ブロック化部２０２、離散フーリエ変換部２０３、振幅計算部２０４、振幅積算部２０５、回転角度検出・パターン番号識別部４０１、部分データ番号検出部４０２、振幅合算部４０３、部分データ検出部４０４及びデータ検出部４０５を備える。なお、これらの部分はハードウェアによって実現してもよいが、コンピュータをこれらの部分として機能させるためのプログラムをコンピュータが読み込んで実行することによって実現してもよい。
【０１１８】
次に、図１６に示す電子透かし検出装置により行われる電子透かし検出方法について、図１６及び１７を参照して説明する。
【０１１９】
この電子透かし検出方法においては、まず、第１の分割画像から第４の分割画像についてステップＳ４５２乃至Ｓ４５８を繰り返す（ステップＳ４５１）。各繰返しにおいては以下のことを行う。
【０１２０】
まず、カメラにより撮影した画像の青成分の現在の分割画像（ｋ番目の分割画像）に対して空間フィルタ２０１を適用する（ステップＳ４５２）。この空間フィルタは、例えば、図１１に示すようなタップ係数を有する。タップ係数から明らかなように、この空間フィルタは帯域通過の特性を有する。従って、画像にこの空間フィルタ２０１を適用することにより、原画像に起因する低周波数成分や高周波数成分を抑圧することができる。
【０１２１】
次に、ブロック化部２０２は、空間フィルタが適用された分割画像を離散フーリエ変換を行うためのサイズの１又は２以上のブロックに分割する（ステップＳ４５３）。
【０１２２】
次に、離散フーリエ変換部２０３は、各ブロックの画像に対して離散フーリエ変換を行い、周波数成分を得る（ステップＳ４５４）。
【０１２３】
次に、振幅計算部２０４は、各ブロックの各周波数成分の振幅を求める（ステップＳ４５５）。
【０１２４】
次に、振幅積算部２０５は、分割画像に含まれるブロックの間で同一の周波数の周波数成分を積算することにより、分割画像についての各周波数成分の積算振幅を得る（ステップＳ４５６）。
【０１２５】
次に、回転角度検出・パターン番号識別部４０１は、各周波数成分の積算振幅を入力し、図１の内円６０２に分布する非ゼロ周波数成分の分布の様子をみることにより、パターン番号ｐを識別する（ステップＳ４５７）。すなわち、図４を参照すると、重心が４５度の位置の近傍と２２５度の位置の近傍にある場合には、現在の分割画像に含まれているパターン番号は１であると判断し、重心が１３５度の位置の近傍と３１５度の位置の近傍にある場合には、現在の分割画像に含まれているパターン番号は２であると判断する。
【０１２６】
次に、回転角度検出用の電子透かしを基に、画像の小回転角度を検出する（ステップＳ４５８）。すなわち、図４を参照すると、現在の繰返しのステップＳ４５７において回転角度検出用電子透かしのパターンがパターン１又は２のうちのどちらであるのかを検出したので、画像が回転していない場合のその番号のパターンの角度を基準として小回転角度を検出する。すなわち、画像が回転していない場合のその番号のパターンの角度と検出されたパターンの角度との差を小回転角度として検出する。
【０１２７】
ステップＳ４５１の繰返しが終了したならば、部分データ番号検出部４０４は、分割画像の番号ｋと第ｋの分割画像を処理した時に検出されたパターン番号ｐとの対応関係全体を基に、何番目の分割画像に何番目の部分データが対応するのか（すなわち、撮影分割画像番号ｋと部分データ番号ｊとの対応関係）を検出する（ステップＳ４５９）。すなわち、図４を参照すると、第１象限で検出されたパターン番号乃至第４象限で検出されたパターン番号が（１、２、２、１）である場合には、第１の分割画像及び第４の分割画像に第１の部分データが対応し、第２の分割画像及び第３の分割画像に第２の部分データが対応し、第１象限で検出されたパターン番号乃至第４象限で検出されたパターン番号が（２、２、１、１）である場合には、第１の分割画像及び第２の分割画像に第１の部分データが対応し、第３の分割画像及び第４の分割画像に第２の部分データが対応し、第１象限で検出されたパターン番号乃至第４象限で検出されたパターン番号が（２、１、１、２）である場合には、第２の分割画像及び第３の分割画像に第１の部分データが対応し、第１の分割画像及び第４の分割画像に第２の部分データが対応し、第１象限で検出されたパターン番号乃至第４象限で検出されたパターン番号が（１、１、２、２）である場合には、第３の分割画像及び第４の分割画像に第１の部分データが対応し、第１の分割画像及び第２の分割画像に第２の部分データが対応する。
【０１２８】
なお、データ２分割に対応し、１８０度単位のカメラアングルのバリエーションに対応した図１０を参照して説明した方法では、回転角度検出用の非ゼロ周波数成分の重心のみから直ちに部分データの番号を検出することができ、データ４分割に対応した図１５を参照して説明した方法では、回転角度検出用の非ゼロ周波数成分の重心と現在処理している分割画像の象限の両者との対応関係から部分データの番号を検出することができたが、データ２分割に対応し、９０度単位のカメラアングルのバリエーションに対応した当該方法では、回転角度検出用の非ゼロ周波数成分の重心と現在処理している分割画像の象限の両者との対応関係から部分データの番号を検出することはできず、全ての分割画像についての回転角度検出用の非ゼロ周波数成分の重心の全体からどの分割画像にどの部分データが対応しているのかを検出することができる。
【０１２９】
次に、分割画像の番号ｋと第ｋの分割画像を処理した時に検出されたパターン番号ｐとの対応関係全体を基に、画像の回転角度が０度付近であるのか、９０度付近であるのか、１８０度付近であるのか、又は、２７０度付近であるのかを検出し、それに応じて、０度、９０度、１８０度又は２７０度を大回転角度として検出する（ステップＳ４６０）。
【０１３０】
次に、小回転角度と大回転角度を合わせたものを回転角度とする（ステップＳ４６１）。
【０１３１】
ステップＳ４６１の次に、第１の部分データから第２の部分データについてステップＳ４６３及びＳ４６４を繰り返す（ステップＳ４６２）。
【０１３２】
まず、現在の番号の部分データに対応する分割画像全てについての周波数毎の振幅積算値を合算する（ステップＳ４６３）。すなわち、現在の番号が１であり、且つ、回転角度が０度近傍であれば、第１の部分データに対応する第１象限及び第４象限についての周波数毎の振幅積算値を合算し、現在の番号が１であり、且つ、回転角度が９０度近傍であれば、第１の部分データに対応する第１象限及び第２象限についての周波数毎の振幅積算値を合算し、現在の番号が１であり、且つ、回転角度が１８０度近傍であれば、第１の部分データに対応する第２象限及び第３象限についての周波数毎の振幅積算値を合算し、現在の番号が１であり、且つ、回転角度が２７０度近傍であれば、第１の部分データに対応する第３象限及び第４象限についての周波数毎の振幅積算値を合算する。また、現在の番号が２であり、且つ、回転角度が０度近傍であれば、第２の部分データに対応する第２象限及び第３象限についての周波数毎の振幅積算値を合算し、現在の番号が２であり、且つ、回転角度が９０度近傍であれば、第２の部分データに対応する第３象限及び第４象限についての周波数毎の振幅積算値を合算し、現在の番号が２であり、且つ、回転角度が１８０度近傍であれば、第２の部分データに対応する第１象限及び第４象限についての周波数毎の振幅積算値を合算し、現在の番号が２であり、且つ、回転角度が２７０度近傍であれば、第２の部分データに対応する第１象限及び第２象限についての周波数毎の振幅積算値を合算する。
【０１３３】
次に、部分データ検出部４０３は、合算した周波数の振幅から部分データを検出する（ステップＳ４６４）。検出方法は、図１３を参照して説明した方法と同様である。この際、ステップＳ４６１で計算した回転角度をステップＳ９２１で補正角度として利用する。また、現在の繰返し番号ｊを部分データの番号として利用する。
【０１３４】
ステップＳ４６２の繰返しが終了したならば、データ検出部４０５は、全ての繰返しにおけるステップＳ４６４で求めた第１部分データ乃至第２部分データを基にデータを検出する（ステップＳ４６５）。
【０１３５】
なお、空間フィルタ適用（ステップＳ４５２）、ブロック化（ステップＳ４５３）、離散フーリエ変換（ステップＳ４５４）及び振幅計算（ステップＳ４５５）は、第１の繰返しの内部において行う代わりに第１の繰返しに入る前に画像全体に対して行ってもよい。
【０１３６】
図１８に示すように、回転角度を検出するための第１の電子透かしを左画像に挿入し、回転角度を検出するための第２の電子透かしを右画像に挿入し、第１部分データを右上画像に挿入し、第２部分データを左上画像に挿入し、第３部分データを左下画像に挿入し、第４部分データを右下画面に挿入した場合にも、図１６に示す電子透かし検出装置及び図１７を参照して説明した電子透かし検出方法を適用することができる。但し、この場合には、ステップＳ４６１において、１つの部分データに対応する部分画面の数は１つのみである。
【０１３７】
図１８に示すような電子透かしを挿入する方法について図１９を参照して説明する。
【０１３８】
この電子透かし挿入方法においては、第１の分割画像から第Ｎ１の分割画像についてステップＳ６５２乃至Ｓ６５５を繰り返す（ステップＳ６５１）。各繰返しにおいては以下のことを行う。
【０１３９】
まず、回転角度検出用パターン生成部１０１は、第ｊ１の回転角度検出用パターンを生成する（ステップＳ６５２）。このパターンは、非ゼロ周波数成分より構成される。
【０１４０】
次に、逆離散フーリエ変換部１０５は、ステップＳ６５２で生成したパターンに対し、逆離散フーリエ変換をかけることにより、空間領域パターンを得る（ステップＳ６５３）。
【０１４１】
次に、二値化部１０６は、ステップＳ６５３で得た空間領域パターンを二値化する（ステップＳ６５４）。すなわち、各画素について、ステップＳ６５３で得た空間領域パターンの画素値が０以上であれば、所定の値（例えば、１０）とし、ステップＳ６５３で得た空間パターンの画素値が０未満であれば、所定の値（例えば−１０）とする。
【０１４２】
次に、加算部（減算部）１０７は、ステップＳ６５４で二値化された空間領域パターンを原画像の黄成分の第ｊ領域の全部又は一部のブロックに加算又は減算する（ステップＳ６５５）。なお、二値化部１０６が二値化を行うのは、電子透かし検出装置が電子透かしを検出することを容易にするためであり、省略してもよい。この場合には、逆離散フーリエ変換部１０５が出力する空間パターンを二値化せずに直接原画像の黄成分の第ｊ領域に加算又は減算する。
【０１４３】
次に、第１の分割画像から第Ｎ２の分割画像についてステップＳ６５７乃至Ｓ６７０を繰り返す（ステップＳ１５１）。ここで、画像をＮ１個に分割する場合の各分割領域と画像をＮ２個に分割する場合の各分割領域は一般に異なる。各繰返しにおいては以下のことを行う。
【０１４４】
まず、データパターン生成部１０３は、データ保持部１０２に保持されている全体データを読み込み、現在扱っているｊ２の値に応じて第ｊ２のデータを求め、第ｊ２のデータに応じて第ｊ２のデータパターンを生成する（ステップＳ６５７）。このパターンは非ゼロ周波数成分より構成される。第ｊ２のデータパターンが図１８に示すようなものである場合には、第ｊのデータＤｊは、全体データをＤとすると、
Ｄｊ＝ｍｏｄ（Ｄ，２５２＾（Ｎ＋１−ｊ））／２５２＾（Ｎ−ｊ）
により表される。
【０１４５】
次に、逆離散フーリエ変換部１０５は、ステップＳ６５７で生成したパターンに対し、逆離散フーリエ変換をかけることにより、空間領域パターンを得る（ステップＳ６５８）。
【０１４６】
次に、二値化部１０６は、ステップＳ６５８で得た空間領域パターンを二値化する（ステップＳ６５９）。すなわち、各画素について、ステップＳ６５８で得た空間領域パターンの画素値が０以上であれば、所定の値（例えば、１０）とし、ステップＳ１５４で得た空間パターンの画素値が０未満であれば、所定の値（例えば−１０）とする。
【０１４７】
次に、加算部（減算部）１０７は、ステップＳ６５９で二値化された空間領域パターンを原画像の黄成分の第ｊ領域の全部又は一部のブロックに加算又は減算する（ステップＳ６７０）。なお、二値化部１０６が二値化を行うのは、電子透かし検出装置が電子透かしを検出することを容易にするためであり、省略してもよい。この場合には、逆離散フーリエ変換部１０５が出力する空間パターンを二値化せずに直接原画像の黄成分の第ｊ領域に加算又は減算する。
【０１４８】
その後、上記の方法により電子透かしが挿入された画像は印刷物に印刷される。
【０１４９】
なお、ステップＳ６５４、Ｓ６５５、Ｓ６５９及びＳ６７０を削除し、ステップＳ６５６の繰り返しの後に、画像全体に亘り回転角度検出用の電子透かしとデータ用電子透かしを足し合わせ、足し合わせの結果を画像全体に加算又は減算してもよい。加算又は減算の前に必要に応じて足し合わせの結果を二値化してもよい。
【０１５０】
図１８に示す電子透かしは次に説明するような方法によっても検出することができる。
【０１５１】
図２０を参照すると、この電子透かし検出方法においては、まず、第１の分割画像から第４の分割画像についてステップＳ７５２乃至Ｓ７５８を繰り返す（ステップＳ７５１）。各繰返しにおいては以下のことを行う。
【０１５２】
まず、カメラにより撮影した画像の青成分の現在の分割画像（ｋ番目の分割画像）に対して空間フィルタ２０１を適用する（ステップＳ７５２）。この空間フィルタは、例えば、図１１に示すようなタップ係数を有する。タップ係数から明らかなように、この空間フィルタは帯域通過の特性を有する。従って、画像にこの空間フィルタ２０１を適用することにより、原画像に起因する低周波数成分や高周波数成分を抑圧することができる。
【０１５３】
次に、ブロック化部２０２は、空間フィルタが適用された分割画像を離散フーリエ変換を行うためのサイズの１又は２以上のブロックに分割する（ステップＳ７５３）。
【０１５４】
次に、離散フーリエ変換部２０３は、各ブロックの画像に対して離散フーリエ変換を行い、周波数成分を得る（ステップＳ７５４）。
【０１５５】
次に、振幅計算部２０４は、各ブロックの各周波数成分の振幅を求める（ステップＳ７５５）。
【０１５６】
次に、振幅積算部２０５は、分割画像に含まれるブロックの間で同一の周波数の周波数成分を積算することにより、分割画像についての各周波数成分の積算振幅を得る（ステップＳ７５６）。
【０１５７】
次に、回転角度検出・パターン番号識別部４０１は、各周波数成分の積算振幅を入力し、図１の内円６０２に分布する非ゼロ周波数成分の分布の様子をみることにより、パターン番号ｐを識別する（ステップＳ４５７）。すなわち、図１８を参照すると、重心が４５度の位置の近傍と２２５度の位置の近傍にある場合には、現在の分割画像に含まれているパターン番号は１であると判断し、重心が１３５度の位置の近傍と３１５度の位置の近傍にある場合には、現在の分割画像に含まれているパターン番号は２であると判断する。
【０１５８】
次に、回転角度検出用の電子透かしを基に、画像の小回転角度を検出する（ステップＳ７５８）。すなわち、図１８を参照すると、現在の繰返しのステップＳ７５７において回転角度検出用電子透かしのパターンがパターン１又は２のうちのどちらであるのかを検出したので、画像が回転していない場合のその番号のパターンの角度を基準として小回転角度を検出する。すなわち、画像が回転していない場合のその番号のパターンの角度と検出されたパターンの角度との差を小回転角度として検出する。
【０１５９】
ステップＳ７５１の繰返しが終了したならば、部分データ番号検出部４０４は、分割画像の番号ｋと第ｋの分割画像を処理した時に検出されたパターン番号ｐとの対応関係全体を基に、何番目の分割画像に何番目の部分データが対応するのか（すなわち、撮影分割画像番号ｋと部分データ番号ｊとの対応関係）を検出する（ステップＳ７５９）。すなわち、図４を参照すると、第１象限で検出されたパターン番号乃至第４象限で検出されたパターン番号が（１、２、２、１）である場合には、第１の分割画像に第１の部分データが対応し、第２の分割画像に第２の部分データが対応し、第３の分割画像に第３の部分データが対応し、第４の分割画像に第４の部分データが対応し、第１象限で検出されたパターン番号乃至第４象限で検出されたパターン番号が（２、２、１、１）である場合には、第１の部分画像に第４の部分データが対応し、第２の部分画像に第１の部分データが対応し、第３の部分画像に第２の部分データが対応し、第４の部分画像に第３の部分データが対応し、第１象限で検出されたパターン番号乃至第４象限で検出されたパターン番号が（２、１、１、２）である場合には、第１の分割画像に第３の部分データが対応し、第２の分割画像に第４の部分データが対応し、第３の分割画像に第１の部分データが対応し、第４の分割画像に第２の部分データが対応し、第１象限で検出されたパターン番号乃至第４象限で検出されたパターン番号が（１、１、２、２）である場合には、第１の分割画像に第４の部分データが対応し、第２の分割画像に第３の部分データが対応し、第３の分割画像に第２の部分データが対応し、第４の分割画像に第１の部分データが対応する。
【０１６０】
次に、分割画像の番号ｋと第ｋの分割画像を処理した時に検出されたパターン番号ｐとの対応関係全体を基に、画像の回転角度が０度付近であるのか、９０度付近であるのか、１８０度付近であるのか、又は、２７０度付近であるのかを検出し、それに応じて、０度、９０度、１８０度又は２７０度を大回転角度として検出する（ステップＳ７６０）。
【０１６１】
次に、小回転角度と大回転角度を合わせたものを回転角度とする（ステップＳ７６１）。
【０１６２】
ステップＳ７６１の次に、第１の部分データから第４の部分データについてステップＳ７６３を繰り返す（ステップＳ４６２）。
【０１６３】
ステップＳ４６２では、部分データ検出部４０３は、ステップＳ７５６で求めた振幅から部分データを検出する。検出方法は、図１３を参照して説明した方法と同様である。この際、ステップＳ７６１で計算した回転角度をステップＳ９２１で補正角度として利用する。また、現在の繰返し番号ｊを部分データの番号として利用する。
【０１６４】
ステップＳ７６２の繰返しが終了したならば、データ検出部４０５は、全ての繰返しにおけるステップＳ７６３で求めた第１部分データ乃至第４部分データを基にデータを検出する（ステップＳ７６４）。
【０１６５】
なお、空間フィルタ適用（ステップＳ７５２）、ブロック化（ステップＳ７５３）、離散フーリエ変換（ステップＳ７５４）及び振幅計算（ステップＳ７５５）は、第１の繰返しの内部において行う代わりに第１の繰返しに入る前に画像全体に対して行ってもよい。
【０１６６】
図２１は、図５を参照して説明した検出方法（データ検出前画像全体回転補正）を行うための電子透かし検出装置の構成を示すブロック図である。
【０１６７】
図２１を参照すると、この電子透かし検出装置は、空間フィルタ２０１、ブロック化部２０２、離散フーリエ変換部２０３、振幅計算部２０４、振幅積算部２０５、回転角度検出部５０１、回転補正部５０３、空間フィルタ５０４、ブロック化部５０５、離散フーリエ変換部５０６、振幅計算部５０７、振幅積算部５０８、部分データ検出部５０９及びデータ検出部５１０を備える。なお、これらの部分はハードウェアによって実現してもよいが、コンピュータをこれらの部分として機能させるためのプログラムをコンピュータが読み込んで実行することによって実現してもよい。
【０１６８】
次に、図２１に示す電子透かし検出装置により行われる電子透かし検出方法について、図２１及び２２を参照して説明する。
【０１６９】
まず、カメラにより撮影した画像の青成分全体に対して空間フィルタ２０１を適用する（ステップＳ５５１）。この空間フィルタは、例えば、図１１に示すようなタップ係数を有する。タップ係数から明らかなように、この空間フィルタは帯域通過の特性を有する。従って、画像にこの空間フィルタ２０１を適用することにより、原画像に起因する低周波数成分や高周波数成分を抑圧することができる。
【０１７０】
次に、ブロック化部２０２は、空間フィルタが適用された画像全体を離散フーリエ変換を行うためのサイズの１又は２以上のブロックに分割する（ステップＳ５５２）。
【０１７１】
次に、離散フーリエ変換部２０３は、各ブロックの画像に対して離散フーリエ変換を行い、周波数成分を得る（ステップＳ５５３）。
【０１７２】
次に、振幅計算部２０４は、各ブロックの各周波数成分の振幅を求める（ステップＳ５５４）。
【０１７３】
次に、振幅積算部２０５は、画像全体に含まれるブロックの間で同一の周波数の周波数成分を積算することにより、画像全体についての各周波数成分の積算振幅を得る（ステップＳ５５５）。
【０１７４】
次に、回転角度検出部５０１は、各周波数成分の積算振幅を入力し、図１の内円６０２に分布する非ゼロ周波数成分の分布の様子をみることにより、画像全体の回転角度を検出する（ステップＳ５５６）。回転角度を検出するためには図１２を参照して説明した方法を利用することができる。
【０１７５】
次に、回転補正部５０３は、ステップＳ５５６で検出した回転角度を利用して画像全体の回転角度を補正する（ステップＳ５５７）。
【０１７６】
次に、第１の部分データから第Ｎの部分データについてステップＳ５５９乃至Ｓ５６４を繰り返す（ステップＳ５５８）。各繰返しにおいては以下のことを行う。
【０１７７】
まず、回転補正された画像の青成分の現在の分割画像（ｊ番目の分割画像）に対して空間フィルタ５０４を適用する（ステップＳ５５９）。この空間フィルタは、例えば、図１１に示すようなタップ係数を有する。タップ係数から明らかなように、この空間フィルタは帯域通過の特性を有する。従って、画像にこの空間フィルタ２０１を適用することにより、原画像に起因する低周波数成分や高周波数成分を抑圧することができる。
【０１７８】
次に、ブロック化部５０５は、空間フィルタが適用された分割画像を離散フーリエ変換を行うためのサイズのブロックに分割する（ステップＳ５６０）。
【０１７９】
次に、離散フーリエ変換部５０６は、各ブロックの画像に対して離散フーリエ変換を行い、周波数成分を得る（ステップＳ５６１）。
【０１８０】
次に、振幅計算部５０７は、各ブロックの各周波数成分の振幅を求める（ステップＳ５６２）。
【０１８１】
次に、振幅積算部５０８は、分割画像に含まれるブロックの間で同一の周波数の周波数成分を積算することにより、分割画像についての各周波数成分の積算振幅を得る（ステップＳ５６３）。
【０１８２】
次に、各周波数成分の積算振幅から第ｊの部分データを検出する（ステップＳ５６４）。第ｊの部分データを検出するためには、図１３を参照して説明した検出方法を利用することができる。但し、既にステップＳ５５７で回転補正をしているので、ステップＳ９２１を省略する。
【０１８３】
ステップＳ５５８の繰返しが終了したならば、各繰返しのステップＳ５６４で検出した第１部分データ乃至第Ｎ部分データを基に全体データを検出する（ステップＳ５６５）。
【産業上の利用可能性】
【０１８４】
本発明は、印刷物に印刷される画像に電子透かしを挿入し、その画像をカメラ等で撮影し、撮影された画像から電子透かしを検出することにより、電子透かしに対応するデータを取得するために利用することができる。
【図面の簡単な説明】
【０１８５】
【図１】本発明の実施形態で利用するデータ２分割の場合の電子透かしを示す図である。
【図２】本発明の実施形態で利用するデータ２分割・回転角度検出用電子透かし２分割・検出時２分割の場合の印刷画像と撮影画像における電子透かしを示す図である。
【図３】本発明の実施形態で利用するデータ４分割・回転角度検出用電子透かし４分割・検出時４分割の場合の印刷画像と撮影画像における電子透かしを示す図である。
【図４】本発明の実施形態で利用するデータ２分割・回転角度検出用電子透かし２分割・検出時４分割の場合の印刷画像と撮影画像における電子透かしを示す図である。
【図５】本発明の実施形態で利用するデータ検出前画像全体回転補正の場合の印刷画像における電子透かしを示す図である。
【図６】本発明の実施形態による電子透かし挿入装置の構成を示すブロック図である。
【図７】本発明の実施形態による電子透かし挿入方法を示すフローチャートである。
【図８】本発明の他の実施形態による電子透かし挿入方法を示すフローチャートである。
【図９】本発明のデータ２分割・回転角度検出用電子透かし２分割・検出時２分割の場合の実施形態による電子透かし検出装置の構成を示すブロック図である。
【図１０】本発明のデータ２分割・回転角度検出用電子透かし２分割・検出時２分割の場合の実施形態による電子透かし検出方法を示すフローチャートである。
【図１１】本発明の実施形態で利用する空間フィルタのタップ係数を示す図である。
【図１２】本発明の実施形態による回転角度検出用電子透かしを検出する方法を示すフローチャートである。
【図１３】本発明の実施形態によるデータを表す電子透かしを検出する方法を示すフローチャートである。
【図１４】本発明のデータ４分割・回転角度検出用電子透かし４分割・検出時４分割の場合の実施形態による電子透かし検出装置の構成を示すブロック図である。
【図１５】本発明のデータ４分割・回転角度検出用電子透かし４分割・検出時４分割の場合の実施形態による電子透かし検出方法を示すフローチャートである。
【図１６】本発明のデータ２分割・回転角度検出用電子透かし２分割・検出時４分割の場合の実施形態による電子透かし検出装置の構成を示すブロック図である。
【図１７】本発明のデータ２分割・回転角度検出用電子透かし２分割・検出時４分割の場合の実施形態による電子透かし検出方法を示すフローチャートである。
【図１８】本発明のデータ４分割・回転角度検出用電子透かし２分割・検出時４分割の場合の印刷画像と撮影画像における電子透かしを示す図である。
【図１９】本発明の更に他の実施形態による電子透かし挿入方法を示すフローチャートである。
【図２０】本発明のデータ４分割・回転角度検出用電子透かし２分割・検出時４分割の場合の実施形態による電子透かし検出方法を示すフローチャートである。
【図２１】本発明のデータ検出前画像全体回転補正の場合の実施形態による電子透かし検出装置の構成を示すブロック図である。
【図２２】本発明のデータ検出前画像全体回転補正の場合の実施形態による電子透かし検出方法を示すフローチャートである。
【符号の説明】
【０１８６】
１０１  回転角度検出用パターン生成部
１０２  データ保持部
１０３  データパターン生成部
１０４  合成部
１０５  逆離散フーリエ変換部
１０６  二値化部
１０７  加算部（減算部）
２０１  空間フィルタ
２０２  ブロック化部
２０３  離散フーリエ変換部
２０４  振幅計算部
２０５  振幅積算部
２０６  回転角度検出・領域番号識別部
２０７  部分データ検出部
２０８  データ検出部
３０１  回転角度検出・パターン番号識別部
３０２  部分データ番号検出部
３０３  部分データ検出部
３０４  データ検出部
４０１  回転角度検出・パターン番号識別部
４０２  部分データ番号検出部
４０３  振幅合算部
４０４  部分データ検出部
４０５  データ検出部
５０１  回転角度検出部
５０３  回転補正部
５０４  空間フィルタ
５０５  ブロック化部
５０６  離散フーリエ変換部
５０７  振幅計算部
５０８  振幅積算部
５０９  部分データ検出部
５１０  データ検出部
【図１】
2006050069.tif 000003
【図２】
2006050069.tif 000004
【図３】
2006050069.tif 000005
【図４】
2006050069.tif 000006
【図５】
2006050069.tif 000007
【図６】
2006050069.tif 000008
【図７】
2006050069.tif 000009
【図８】
2006050069.tif 000010
【図９】
2006050069.tif 000011
【図１０】
2006050069.tif 000012
【図１１】
2006050069.tif 000013
【図１２】
2006050069.tif 000014
【図１３】
2006050069.tif 000015
【図１４】
2006050069.tif 000016
【図１５】
2006050069.tif 000017
【図１６】
2006050069.tif 000018
【図１７】
2006050069.tif 000019
【図１８】
2006050069.tif 000020
【図１９】
2006050069.tif 000021
【図２０】
2006050069.tif 000022
【図２１】
2006050069.tif 000023
【図２２】
2006050069.tif 000024