"入手画像を周波数領域の信号に変換し、その周波数領域の信号とオリジナル画像に埋め込まれている位置合わせ信号を構成している複数の位置合わせデータが示す周波数平面上の座標とを照合して、上記オリジナル画像に対して、上記入手画像がどのように拡大縮小、回転、反転されているかを特定する拡大縮小等特定手段と、上記拡大縮小等特定手段の特定結果を用いて、上記複数の位置合わせデータの中から、上記周波数領域の信号に対応する位置合わせデータを特定し、上記位置合わせデータから、上記オリジナル画像に対して、上記入手画像がどのように平行移動しているかを特定し、その特定結果から上記オリジナル画像と上記入手画像の同期位置を特定する同期位置特定手段と、上記同期位置特定手段により特定された同期位置と上記位置合わせ信号を用いて、上記拡大縮小又は上記回転の少なくとも一方の微調整量を算出する微調整量算出手段とを備えた画像位置合わせ装置。入手画像の一部を構成している部分画像を切り出す画像切出手段を設け、拡大縮小等特定手段は、上記画像切出手段により切り出された部分画像を周波数領域の信号に変換し、その周波数領域の信号と位置合わせ信号を構成している複数の位置合わせデータが示す周波数平面上の座標とを照合して、オリジナル画像に対して、上記入手画像がどのように拡大縮小、回転、反転されているかを特定することを特徴とする請求項1記載の画像位置合わせ装置。拡大縮小等特定手段は、画像切出手段により切り出された部分画像をフーリエ変換して、その変換結果であるフーリエ係数を出力するフーリエ変換部と、上記フーリエ変換部から出力されたフーリエ係数と位置合わせ信号を構成している複数の位置合わせデータである周波数スペクトルの周波数平面上の座標とを照合して、オリジナル画像に対して、入手画像がどのように拡大縮小、回転、反転されているかを特定する拡大縮小回転反転特定部とから構成されており、同期位置特定手段は、上記拡大縮小回転反転特定部の特定結果を用いて、上記位置合わせ信号を構成している複数の周波数スペクトルの中から、上記フーリエ変換部から出力されたフーリエ係数に対応する周波数スペクトルを特定し、上記周波数スペクトルの位相から、上記オリジナル画像に対して、上記入手画像がどのように平行移動しているかを特定することを特徴とする請求項2記載の画像位置合わせ装置。微調整量算出手段により拡大縮小の微調整量が算出された場合、上記微調整量を用いて、上記拡大縮小を補正する拡大縮小補正手段を設け、同期位置特定手段は、拡大縮小等特定手段により特定された拡大縮小の代わりに、上記拡大縮小補正手段により補正された拡大縮小を用いて、周波数領域の信号に対応する位置合わせデータを特定することを特徴とする請求項1から請求項3のうちのいずれか1項記載の画像位置合わせ装置。微調整量算出手段により回転の微調整量が算出された場合、上記微調整量を用いて、上記回転を補正する回転補正手段を設け、同期位置特定手段は、拡大縮小等特定手段により特定された回転の代わりに、上記回転補正手段により補正された回転を用いて、周波数領域の信号に対応する位置合わせデータを特定することを特徴とする請求項1から請求項3のうちのいずれか1項記載の画像位置合わせ装置。入手画像を周波数領域の信号に変換し、その周波数領域の信号とオリジナル画像に埋め込まれている位置合わせ信号を構成している複数の位置合わせデータが示す周波数平面上の座標とを照合して、上記オリジナル画像に対して、上記入手画像がどのように拡大縮小、回転、反転されているかを特定する拡大縮小等特定処理手順と、上記拡大縮小等特定処理手順での特定結果を用いて、上記複数の位置合わせデータの中から、上記周波数領域の信号に対応する位置合わせデータを特定し、上記位置合わせデータから、上記オリジナル画像に対して、上記入手画像がどのように平行移動しているかを特定し、その特定結果から上記オリジナル画像と上記入手画像の同期位置を特定する同期位置特定処理手順と、上記同期位置特定処理手順により特定された同期位置と上記位置合わせ信号を用いて、上記拡大縮小又は上記回転の少なくとも一方の微調整量を算出する微調整量算出処理手順とをコンピュータに実行させるための画像位置合わせプログラム。入手画像を周波数領域の信号に変換し、その周波数領域の信号とオリジナル画像に埋め込まれている位置合わせ信号を構成している複数の位置合わせデータが示す周波数平面上の座標とを照合して、上記オリジナル画像に対して、上記入手画像がどのように拡大縮小、回転、反転されているかを特定する拡大縮小等特定手段と、上記拡大縮小等特定手段の特定結果を用いて、上記複数の位置合わせデータの中から、上記周波数領域の信号に対応する位置合わせデータを特定し、上記位置合わせデータから、上記オリジナル画像に対して、上記入手画像がどのように平行移動しているかを特定し、その特定結果から上記オリジナル画像と上記入手画像の同期位置を特定する同期位置特定手段と、上記同期位置特定手段により特定された同期位置と上記位置合わせ信号を用いて、上記拡大縮小又は上記回転の少なくとも一方の微調整量を算出する微調整量算出手段と、上記微調整量算出手段により算出された拡大縮小の微調整量を用いて、上記拡大縮小を補正するとともに、上記微調整量算出手段により算出された回転の微調整量を用いて、上記回転を補正する拡大縮小回転補正手段と、上記拡大縮小回転補正手段により補正された拡大縮小及び回転、上記反転及び上記同期位置を用いて、上記入手画像を補正する画像補正手段と、上記画像補正手段により補正された入手画像から電子透かしを検出する電子透かし検出手段とを備えた電子透かし検出装置。従来の画像位置合わせ装置は以上のように構成されているので、拡大縮小等の位置合わせの結果に誤差が生じていても、その誤差を補償する手段がなく、画像の位置合わせの精度が低下してしまう課題があった。また、拡大縮小等の位置合わせの結果に誤差があると、その誤差の影響が平行移動に関する位置合わせに及んで、平行移動に関する位置合わせの結果に大きな誤差が発生してしまうことがある課題もあった。この発明は上記のような課題を解決するためになされたもので、拡大縮小等の位置合わせの結果に生じている誤差を補償して、画像の位置合わせ精度を高めることができる画像位置合わせ装置及び画像位置合わせプログラムを得ることを目的とする。また、この発明は、画像の位置合わせ精度を高めて、入手画像から正確に電子透かしを検出することができる電子透かし検出装置を得ることを目的とする。この発明に係る画像位置合わせ装置は、入手画像を周波数領域の信号に変換し、その周波数領域の信号とオリジナル画像に埋め込まれている位置合わせ信号を構成している複数の位置合わせデータが示す周波数平面上の座標とを照合して、オリジナル画像に対して、入手画像がどのように拡大縮小、回転、反転されているかを特定する拡大縮小等特定手段と、拡大縮小等特定手段の特定結果を用いて、複数の位置合わせデータの中から、上記周波数領域の信号に対応する位置合わせデータを特定し、その位置合わせデータから、オリジナル画像に対して、入手画像がどのように平行移動しているかを特定し、その特定結果からオリジナル画像と入手画像の同期位置を特定する同期位置特定手段とを設け、微調整量算出手段が、同期位置特定手段により特定された同期位置と位置合わせ信号を用いて、上記拡大縮小又は上記回転の少なくとも一方の微調整量を算出するようにしたものである。この発明によれば、拡大縮小等特定手段の特定結果を用いて、複数の位置合わせデータの中から、周波数領域の信号に対応する位置合わせデータを特定し、その位置合わせデータから、オリジナル画像に対して、入手画像がどのように平行移動しているかを特定し、その特定結果からオリジナル画像と入手画像の同期位置を特定する同期位置特定手段を設け、微調整量算出手段が、同期位置特定手段により特定された同期位置と位置合わせ信号を用いて、オリジナル画像に対する入手画像の拡大縮小又は回転の少なくとも一方の微調整量を算出するように構成したので、オリジナル画像に対する入手画像の拡大縮小や回転の誤差を補償して、画像の位置合わせ精度を高めることができる効果がある。実施の形態1.図1はこの発明の実施の形態1による画像位置合わせ装置を示す構成図である。図1において、部分画像切出部1は拡大縮小、回転、反転、切取などの編集が施されている可能性があるデジタル画像非圧縮輝度画像であって、位置合わせ信号が電子透かしとして埋め込まれているデジタル画像以下、入手画像と称するを入手すると、その入手画像の一部を構成している部分画像を切り出し、その部分画像をフーリエ変換部2に出力する処理を実施する。なお、部分画像切出部1は画像切出手段を構成している。図1の例では、部分画像切出部1が画像位置合わせ装置に実装されているが、例えば、入手画像のサイズが、後段の処理部で取扱いが可能なサイズであるような場合には、部分画像切出部1が画像位置合わせ装置に実装されていなくてもよい。フーリエ変換部2は部分画像切出部1により切り出された部分画像をフーリエ変換して、その変換結果であるフーリエ係数を拡大縮小回転反転特定部3及び同期位置特定部4に出力する処理を実施する。拡大縮小回転反転特定部3はオリジナル画像拡大縮小等が施されていない入手画像に埋め込まれている位置合わせ信号を構成している複数の周波数スペクトル位置合わせデータを入力すると、フーリエ変換部2から出力されたフーリエ係数と、複数の周波数スペクトルの周波数平面上の座標とを照合して、オリジナル画像に対して、入手画像がどのように拡大縮小、回転、反転されているかを特定する。即ち、オリジナル画像に対する入手画像の拡大率、回転角及び反転状況を特定する処理を実施する。なお、フーリエ変換部2及び拡大縮小回転反転特定部3から拡大縮小等特定手段が構成されている。同期位置特定部4は拡大縮小回転反転特定部3により特定された拡大率、回転角及び反転状況を用いて、その位置合わせ信号を構成している複数の周波数スペクトルの中から、フーリエ変換部2から出力されたフーリエ係数に対応する周波数スペクトルを特定し、その周波数スペクトルの位相から、オリジナル画像に対して、入手画像がどのように平行移動しているかを特定する。即ち、オリジナル画像に対する入手画像の平行移動量を算出する処理を実施する。また、同期位置特定部4はオリジナル画像に対する入手画像の平行移動量からオリジナル画像と入手画像の同期位置を特定する処理を実施する。なお、同期位置特定部4は同期位置特定手段を構成している。微調整量算出部5は同期位置特定部4により特定された同期位置と位置合わせ信号を用いて、拡大縮小回転反転特定部3により特定された拡大率拡大縮小及び回転角回転の微調整量を算出する処理を実施する。なお、微調整量算出部5は微調整量算出手段を構成している。この実施の形態1では、微調整量算出部5が拡大率の微調整量と、回転角の微調整量の双方を算出する例を説明するが、拡大率又は回転角の少なくとも一方の微調整量を算出するようにしてもよい。虹斑模様図1の例では、画像位置合わせ装置の構成要素である部分画像切出部1、フーリエ変換部2、拡大縮小回転反転特定部3、同期位置特定部4及び微調整量算出部5のそれぞれが専用のハードウェア例えば、CPUを実装している半導体集積回路、あるいは、ワンチップマイコンで構成されているものを想定しているが、画像位置合わせ装置がコンピュータで構成されていてもよい。画像位置合わせ装置がコンピュータで構成されている場合、部分画像切出部1、フーリエ変換部2、拡大縮小回転反転特定部3、同期位置特定部4及び微調整量算出部5の処理内容を記述している画像位置合わせプログラムをコンピュータのメモリに格納し、当該コンピュータのCPUが当該メモリに格納されているプログラムを実行するようにすればよい。図2はこの発明の実施の形態1による画像位置合わせ装置の処理内容を示すフローチャートである。次に動作について説明する。最初に、画像位置合わせ装置による画像の位置合わせの原理について説明する。即ち、前提となる2次元離散フーリエ変換以下、単にフーリエ変換と称するによる位置合わせの理論について説明する。まず、オリジナル画像に埋め込まれる位置合わせ信号は、横Mピクセル、縦Mピクセルの正方形のパターンであり、s個の周波数スペクトルSii∈[1,s]を含むパターンであるとする。位置合わせ信号のパターンは、一般的には長方形でもよいが、実装の多くで正方形が用いられるので、ここでは、正方形であるとする。また、フーリエ変換サイズはMであるとし、周波数スペクトルSiSiは複素数は周波数平面上でXi,YiXi∈[1,M/2-1]の位置に存在しているものとする。オリジナル画像に足し合わせる位置合わせ信号のパターンは実数であるから、必然的に周波数平面上の-Xi,-Yiの位置に周波数スペクトルSiの共役複素数の周波数スペクトルが配置される。即ち、位置合わせ信号のパターンをfx,yx∈[-M/2,M/2-1]y∈[-M/2,M/2-1]とすると、fx,yのフーリエ変換FX,YX∈[-M/2,M/2-1],Y∈[-M/2,M/2-1]は、下記のようになる。FX,Y=SiX,Y∈{Xi,Yi|i∈[1,s]}Siの共役複素数X,Y∈{-Xi,-Yi|i∈[1,s]}0その他のX,Yここで、フーリエ変換FX,Y及び逆フーリエ変換fx,yは、下記のように表される。この実施の形態1では、横Mピクセル、縦Mピクセルの位置合わせ信号のパターンfx,yが、循環的にオリジナル画像に加算されているものとする。画像の特性に応じて、検出精度や主観画質を向上させるために、画素毎に重みを調節しながら加算されているものであってもよい。この実施の形態1では、位置合わせ信号のパターンfx,yの左上画素を、同期位置と呼ぶことにする。この場合、同期位置は、左上から横方向に第M×p1画素、縦方向に第M×q1画素に現れるp,qは非負整数。図3はオリジナル画像に埋め込まれている位置合わせ信号を示す説明図であり、上記の内容を図解したものである。画像の位置合わせは、入手画像の任意の横Mピクセル、縦Mピクセルの正方形部分を対象に行われる。その正方形部分をフーリエ変換して、周波数スペクトルの絶対値振幅を調べることで、位置合わせ用の周波数スペクトルの位置を検出することができる。この検出処理は、周波数スペクトルの絶対値マップから、値が突出している座標をピックアップする作業であり、公知の画像処理技術によって十分可能であるため、その方法の詳細については割愛する。検出された周波数スペクトルの数がt、周波数スペクトルがTkk∈[1,t]であるとする。ここで、入手画像がある角度だけ回転されていれば、s個の周波数スペクトルSiの周波数平面上での位置が、原点を中心に同じ角度だけ回転される。また、入手画像がある対称線を軸に反転されていれば、s個の周波数スペクトルSiの周波数平面上での位置が、同じ角度で原点を通る対称線を軸に反転される。また、入手画像がある拡大率で拡大されていれば、s個の周波数スペクトルSiの周波数平面上での位置が、原点を中心にその拡大率の逆数の拡大率で拡大される周波数スペクトルSiの周波数平面上での原点からの距離は、画像の拡大率の逆数倍になる。この性質を利用することで、入手画像がオリジナル画像に対して、どのように拡大縮小、回転、反転されているかを判断することができる。これらの方法の詳細については、上記の特許文献1,2に開示されている。s個の周波数スペクトルSiの周波数平面上での位置から判断される拡大率、回転角及び反転軸の角度に基づいて入手画像を補正することで、平行移動の不確定性だけが残された状態まで補正することができる。ここでは、平行移動の不確定性だけが残された状態まで補正された入手画像を半補正画像と称することにする。半補正画像は、オリジナル画像からx軸方向にmピクセルm∈[-M/2,M/2-1]、y軸方向にnピクセルn∈[-M/2,M/2-1]だけ平行移動されているものとする。このm,nを求めるのが、平行移動の位置合わせである。拡大縮小、回転、反転が施されている場合の他、拡大縮小、回転、反転が施されていない場合であっても、切取トリミングが施されている場合は、平行移動の位置合わせが必要になる。ここで、周波数スペクトルTkに対応する周波数スペクトルがSiTkとSiの対応関係は、拡大縮小、回転、反転に対する位置合わせの手続で明らかになっているであるとし、その周波数スペクトルSiの位相をPiSi=|Si|expjPi、その周波数スペクトルTkの位相をQkTk=|Tk|expjQkとすると、下記の式が成立するmとnを選択する。フーリエ変換サイズMのフーリエ変換に基づく平行移動の位置合わせでは、真の平行移動量を得ることはできず、得ることができるのは、真の平行移動量のMを法とする値のみである。これは、オリジナル画像に横Mピクセル、縦Mピクセルの位置合わせ信号のパターンfx,yを循環的に適用しているため、横方向のMの倍数ピクセル、縦方向にMの倍数ピクセルだけ移動した位置の位置合わせの信号パターンが全て同じになるからである。通常、フーリエ変換を利用する位置合わせでは、真の平行移動量のMを法とする値を得ることが目的であり、真の平行移動量の算出は目的外である。即ち、画像の位置合わせは、オリジナル画像から入手画像に至るまでになされた編集の、拡大、回転、反転についてのパラメータ拡大率、回転角等を明らかにすることと、Mを法とする平行移動量を明らかにすることである。また、周期的に現れる同期位置の座標を明らかにすることと考えることもできる。これらは等価であって、いずれの解釈であっても、オリジナル画像において、縦横共にMの倍数ピクセルの平行移動の不確定性があることを除けば、オリジナル画像と入手画像との対応関係を明らかにしたことになる。次に、画像位置合わせ装置の処理内容を具体的に説明する。まず、部分画像切出部1は、拡大縮小、回転、反転、切取などの編集が施されている可能性がある入手画像を入力すると、その入手画像の一部を構成している部分画像を切り出し、その部分画像をフーリエ変換部2に出力する図2のステップST1。この実施の形態1では、横Mピクセル、縦Mピクセルの正方形の部分画像を切り出すものとする。部分画像の切り出し位置は、ランダムに選択してもよいし、左上隅を選択してもよいし、様々な実装が考えられる。また、部分画像は、入手画像に完全に含まれていてもよいし、一部が入手画像に含まれて、残りの部分が一定画素値でパディングされていてもよい。また、入手画像が部分画像のサイズに満たない場合には、必然的にパディングによって横Mピクセル、縦Mピクセルとすることが必要になる。ここで、図4はオリジナル画像の一部を含む入手画像から4つの部分画像が切り出されている例を示す説明図である。フーリエ変換部2は、部分画像切出部1から横Mピクセル、縦Mピクセルの部分画像を受けると、その部分画像に対して、M×Mの2次元離散フーリエ変換を実施して、その変換結果であるフーリエ係数を拡大縮小回転反転特定部3及び同期位置特定部4に出力するステップST2。拡大縮小回転反転特定部3は、オリジナル画像に埋め込まれている位置合わせ信号を構成している複数の周波数スペクトル位置合わせデータを入力すると、フーリエ変換部2から出力されたフーリエ係数と、複数の周波数スペクトルの周波数平面上の座標とを照合して、オリジナル画像に対する入手画像の拡大率、回転角及び反転状況を特定し、拡大率を示す拡大縮小情報、回転角を示す回転情報及び反転状況を示す反転情報を出力するステップST3。フーリエ係数と複数の周波数スペクトルの周波数平面上の座標とを照合して、拡大率、回転角及び反転状況を特定する処理自体は公知の技術であり、この技術は特許文献1,2にも開示されているので、詳細な説明は省略する。拡大縮小回転反転特定部3から出力された拡大縮小情報、回転情報及び反転情報は、後段の同期位置特定部4に入力されるが、その拡大縮小情報には、オリジナル画像に対する入手画像の拡大率の他に、拡大縮小の有無や、拡大縮小の中心となる位置の座標を示す情報などを含めることができる。また、回転情報には、オリジナル画像に対する入手画像の回転角の他に、回転の有無や、回転の中心となる位置の座標を示す情報などを含めることができる。また、反転情報には、オリジナル画像に対する入手画像の反転の有無の他に、反転軸の勾配や、反転軸が通る座標を示す情報などを含めることができる。拡大縮小情報、回転情報及び反転情報は、入手画像の拡大率、回転角及び反転状況を特定することができれば、どのようなパラメータ群によって構成されていてもよい。例えば、拡大縮小が施されていない場合、拡大率を1で表すことで、拡大縮小の有無に関するパラメータを省略することができる。また、回転が施されていない場合、回転角を0で表すことで、回転の有無に関するパラメータを省略することができる。また、反転軸は、反転軸が通る2点の座標をパラメータとすることもできるし、反転軸が通る1点の座標と反転軸の勾配をパラメータとすることもできる。この実施の形態1では、説明の簡単化のため、拡大縮小情報はオリジナル画像に対する入手画像の拡大率だけを含み、回転情報はオリジナル画像に対する入手画像の回転角だけを含み、反転情報はオリジナル画像に対する入手画像の反転の有無だけを含んでいるものとして説明する。これは、オリジナル画像から入手画像への編集モデルを限定することによって可能となる最小のパラメータ構成である。この実施の形態1では、オリジナル画像から入手画像への編集モデルは、平行移動、特定点を中心とする拡大縮小、特定点を中心とする回転、特定軸による左右反転の順になされる編集であるものとする。このように編集モデルを限定しても一般性を失うことはない平行移動、拡大縮小、回転、反転の組み合わせによる任意の編集を表現することができる。また、この実施の形態1では、フーリエ変換に基づき生成された位置合わせ信号がオリジナル画像に適用されて流通していることを前提とする。同期位置特定部4は、拡大縮小回転反転特定部3から拡大縮小情報、回転情報及び反転情報を受けると、その拡大縮小情報が示す拡大率、回転情報が示す回転角及び反転情報が示す回転の有無を用いて、その位置合わせ信号を構成している複数の周波数スペクトルの中から、フーリエ変換部2から出力されたフーリエ係数に対応する周波数スペクトルを特定する。即ち、同期位置特定部4は、例えば、位置合わせ信号を構成している周波数スペクトルの周波数平面上での座標を、拡大縮小情報が示す拡大率の逆数の拡大率だけ拡大し、回転情報が示す回転角だけ回転し、反転情報が反転有であれば、反転する処理を実施する。そして、上記処理後のオリジナル画像に埋め込まれている位置合わせ信号を構成している複数の周波数スペクトルの中から、フーリエ変換部2から出力されたフーリエ係数に対応する周波数スペクトルを特定するフーリエ係数の周波数平面上でのピークと、位置合わせ信号における周波数スペクトルとの対応関係を明らかにする。同期位置特定部4は、位置合わせ信号を構成している複数の周波数スペクトルの中から、フーリエ変換部2から出力されたフーリエ係数の周波数平面上でのピークに対応する周波数スペクトルを特定すると、その周波数スペクトルの位相から、オリジナル画像に対する入手画像の平行移動量を算出するステップST4。フーリエ係数と周波数スペクトルの位相から平行移動量を算出する技術自体は公知の技術であり、この技術は、上記の特許文献1にも開示されているため、詳細な説明は省略する。なお、同期位置特定部4により算出される平行移動量は、平行移動の大きさだけでなく、平行移動の向きなどを表すパラメータも含んでいるものとする。ただし、この平行移動量を示すパラメータは、オリジナル画像に対する入手画像の平行移動を特定することができれば、どのようなパラメータ群によって構成されていてもよい。例えば、平行移動の向きを表すベクトルの仰角と平行移動の大きさをパラメータとすることもできるし、横方向移動量と縦方向移動量をパラメータとすることもできる。この実施の形態1では、オリジナル画像に対する入手画像の平行移動量は、右向きを正とする横方向移動量を示すパラメータと、下向きを正とする縦方向移動量を示すパラメータから構成されているものとする。同期位置特定部4は、オリジナル画像に対する入手画像の平行移動量を算出すると、その平行移動量からオリジナル画像と入手画像の同期位置を特定するステップST5。オリジナル画像と入手画像の同期位置を示すパラメータについては、どのようなパラメータ群によって構成されていてもよい。この実施の形態1では、オリジナル画像と入手画像の同期位置は、半補正画像平行移動の不確定性だけが残された状態まで補正された入手画像に対して、右向きにx軸、下向きにy軸を設定した際の座標で表すものとする。同期位置を示す座標は、相対的な位置関係しか意味を持たないため、座標軸を設定する際の原点の位置は任意でよい。また、ここでは、同期位置を座標軸の設定による同期位置座標としているが、同期位置を特定することができれば、任意のパラメータ群を使用してよい。なお、平行移動量から同期位置を示す座標を求める際は、同期位置を示す座標が前提とする座標軸に、半補正画像を当て嵌め、部分画像の起点から平行移動量だけ移動した点を同期位置の座標とすればよい。ここで、図5は同期位置座標が前提とする座標軸において、部分画像、平行移動量、同期位置の関係の一例を示す説明図である。微調整量算出部5は、同期位置特定部4が幾つかの部分画像について同期位置を特定するとステップST6、それらの同期位置と位置合わせ信号を用いて、拡大縮小回転反転特定部3により特定された拡大率及び回転角の微調整量を算出するステップST7。なお、同期位置特定部4は、部分画像毎に同期位置を特定する必要があるが、拡大縮小回転反転特定部3は、部分画像毎に拡大率等を特定する処理を実施するようにしてもよいし、最初の部分画像についてだけ拡大率等を特定する処理を実施するようにしてもよい。拡大縮小回転反転特定部3が、最初の部分画像についてだけ拡大率等を特定する場合、同期位置特定部4では、他の部分画像について同期位置を特定する際、最初の部分画像について特定された拡大率等を流用する。図5では、部分画像1から部分画像4について同期位置の座標を算出している例を示している。仮に、位置合わせが完璧であるとすれば、同期位置を示す座標P1,P2,P3,P4は、Mピクセル間隔で整然と並ぶはずである。ところが実際は、算出誤差があるため、少しずれている。微調整量算出部5では、M図5の例では、M=256ピクセル間隔で整然と並ぶべき同期位置のずれは、拡大縮小回転反転特定部3により特定された拡大率及び回転角の誤差に起因すると考えて、その拡大率及び回転角の微調整量を算出するようにしている。以下、微調整量算出部5による微調整量の算出処理を具体的に説明する。最初に、拡大率の微調整量を算出する場合について説明する。拡大率の微調整量を算出する場合、P1とP2間の長さと、M=256との比を微調整量とすることが考えられる。図5の例では、P1とP2間の長さは約258.008であり、M=256の1.008倍である。この場合、半補正画像は、原画像の1.008倍であると考え、1.008倍を拡大率の微調整量とする。また、精度を高めるために、P1とP3間の長さや、P3とP4間の長さを考慮することもできる。P1とP3間の長さは約258.002であり、M=256の1.000008倍である。また、P3とP4間の長さは約257.008であり、M=256の1.004倍である。このことから、例えば、1.008倍、1.000008倍、1.004倍の平均をとって、約1.004倍を拡大率の微調整量とすることができる。このように、より多くのデータの平均をとるなどして、精度を高めることができる。一方、拡大率の微調整量を1倍とするべきケースも考えられる。例えば、P1とP2間の長さ、P1とP3間の長さ、P3とP4間の長さが、M=256より大きかったり、小さかったりで、まちまちの場合、拡大率の誤差に起因すると考えるより、平行移動の誤差に起因して同期位置がばらついていると考えるべきである。平行移動の誤差に起因して同期位置がばらついていると考えられる場合、拡大率を微調整しない判断を行うことで、位置合わせの精度を高めることができる。ばらつきの判断は、例えば、分散のように、ばらつきを示す統計量を計算し、閾値との比較によって行うことができる。ここで、図6は拡大率の微調整量の算出処理を示すフローチャートである。図6の例では、微調整量算出部5は、複数の隣接する同期位置間の距離として、例えば、P1とP2間の長さや、P1とP3間の長さなどを算出しステップST11、複数の隣接する同期位置間の距離のばらつき例えば、分散を算出するステップST12。微調整量算出部5は、複数の隣接する同期位置間の距離のばらつきと所定の閾値を比較し、そのばらつきが閾値未満であればステップST13、複数の隣接する同期位置間の距離の平均値を算出しステップST14、その平均値とMとの比を拡大率の微調整量とするステップST15。一方、そのばらつきが閾値以上であればステップST13、拡大率の微調整量を1とするステップST16。次に、回転角の微調整量を算出する場合について説明する。回転角の微調整を算出する場合、P1からP2を見る向きと、x軸正の向きとの成す角を回転角の微調整とすることが考えられる。図5の例では、ベクトルP1P2は2582であり、x軸正の向きとの成す角が約0.008ラジアンである。この場合、半補正画像は、原画像に対して0.008ラジアン回転していると考え、0.008ラジアンを回転角の微調整量とする。また、精度を高めるために、ベクトルP1P3や、ベクトルP3P4を考慮することもできる。ベクトルP1P3はy軸正の向きから0.004ラジアン回転しており、ベクトルP3P4はx軸正の向きから約0.008ラジアン回転している。このことから、例えば、0.008ラジアン、0.004ラジアン、0.008ラジアンの平均をとって、約0.0067ラジアンを回転角の微調量とすることができる。このように、より多くのデータの平均をとるなどして、精度を高めることができる。一方、回転角の微調整量を1倍とするべきケースも考えられる。例えば、上記の成す角が0ラジアンより大きかったり、小さかったりで、まちまちの場合、回転角の誤差に起因すると考えるより、平行移動の誤差に起因して同期位置がばらついていると考えるべきである。平行移動の誤差に起因して同期位置がばらついていると考えられる場合、回転角を微調整しない判断を行うことで、位置合わせの精度を高めることができる。ばらつきの判断は、例えば、分散のように、ばらつきを示す統計量を計算し、閾値との比較によって行うことができる。ここで、図7は回転角の微調整量の算出処理を示すフローチャートである。図7の例では、微調整量算出部5は、複数の隣接する同期位置間を結ぶベクトル例えば、ベクトルP1P2や、ベクトルP1P3などと、x軸又はy軸との成す角を算出すしステップST21、複数の成す角のばらつき例えば、分散を算出するステップST22。微調整量算出部5は、複数の成す角のばらつきと所定の閾値を比較し、そのばらつきが閾値未満であればステップST23、複数の成す角の平均値を算出しステップST24、その平均値を回転角の微調整量とするステップST25。一方、そのばらつきが閾値以上であればステップST23、回転角の微調整量を0とするステップST26。ここまでは、拡大率の微調整量の算出処理と回転角の微調整量の算出処理を説明したが、拡大率の微調整と回転角の微調整を組み合わせて、拡大率及び回転角の微調整量を出力するようにしてもよい。上述したように、同期位置特定部4により特定された複数の同期位置の関係から判断して、拡大率と回転角のそれぞれについて個別に判断し、微調整が必要である考えられる方についてのみ微調整量を算出して出力することもできる。また、微調整量算出部5は、更に洗練された方法として、最小二乗法を適用して、拡大率及び回転角の微調整量を算出するようにしてもよい。最小二乗法に基づいて計算すると、拡大率の微調整量r倍と回転角の微調整量θラジアンは、下記のようになる。即ち、同期位置の誤差が専ら平行移動の誤差に起因する場合、上式では0に近いθや、1に近いrを返す一方、同期位置の誤差の多くが拡大率や回転角の誤差に起因する場合、最も確からしい微調整量としてのθやrを得ることができる。勿論、統計量や閾値を用いた微調整の要否判断をした後に、上記の最小二乗法に基づく微調整量の算出処理を行ってもよい。以下、上式を用いて、拡大率の微調整量r倍と回転角の微調整量θラジアンを算出する例を説明する。Pは578.25279である。Pi’は真の同期位置であるため、ここではP1’00、P2’2560、P3’0256、P4’256256とする。実際には、検出された同期位置であるP1は、オリジナル画像の左上端から右に256ピクセルの位置であるため、P1’2560が正解であるが、256の倍数ピクセルの不確定性については、この実施の形態1では問題としないため、P1’からP4’を上記のように設定することに何ら問題はない。こうすると、P’は128128である。この場合、r,θは、下記のように算出される。P1’-P’=-128-128P2’-P’=128-128P3’-P’=-128128P4’-P’=128128P1-P=-128.25-129P2-P=129.75-127P3-P=-129.25127P4-P=127.75129P1’-P’×P1-P=96P2’-P’×P2-P=352P3’-P’×P3-P=288P4’-P’×P4-P=160P1’-P’P1-P=32928P2’-P’P2-P=32864P3’-P’P3-P=32800P4’-P’P4-P=32864半補正画像は、原画像に対して約1.003倍に拡大されていて、約0.006816ラジアン回転されていると推定することができる。したがって、微調整量算出部5は、1.003倍を拡大率の微調整量として出力し、0.006816ラジアンを回転角の微調整量として出力する。なお、r、θの算出後に閾値処理を行い、rが1より閾値以上離れている場合や、θが0より閾値以上離れている場合などに限定して、微調整が必要であると判断して微調整量を出力し、さもなければ、微調整を不要とする判断をすることもできる。ここでは、拡大率の微調整量を拡大率で表現し、回転角の微調整量を回転角で表現しているが、それぞれ拡大縮小と回転の微調整を特定できる任意のパラメータ群としてもよい。この実施の形態1では、本発明の主な目的である拡大縮小や回転の位置合わせを微調整する部分に限定し、拡大率の微調整量や回転角の微調整量を出力するものを示したが、実際にオリジナル画像と入手画像の対応関係を算出するには、拡大率の微調整量や回転角の微調整量を得た後に、拡大率の微調整量を用いて、拡大縮小回転反転特定部3により特定された拡大率を補正し、回転角の微調整量を用いて、拡大縮小回転反転特定部3により特定された回転角を補正する必要がある。この実施の形態1では、拡大縮小情報及び拡大率の微調整量は拡大率を示し、回転情報及び回転角の微調整量は回転角を示しているので、単純に乗算あるいは加算することで、微調整済の拡大縮小情報及び微調整済の回転情報を得ることができる。即ち、拡大縮小回転反転特定部3により特定された拡大率と微調整量算出部5により算出された拡大率の微調整量とを乗算することで、微調整済の拡大縮小情報を求め、拡大縮小回転反転特定部3により特定された回転角と微調整量算出部5により算出された回転角の微調整量を加算することで、微調整済の回転情報を求めることができる。このようにして得られた微調整済の拡大縮小情報拡大縮小の微調整を行わない場合は拡大縮小回転反転特定部3から出力された拡大縮小情報及び微調整済の回転情報回転角の微調整を行わない場合は拡大縮小回転反転特定部3から出力された回転情報と、拡大縮小回転反転特定部3から出力された反転情報とを組み合わせることで、オリジナル画像から入手画像に至る過程で拡大縮小、回転、反転がどのようになされたかを特定することができる。更に、同期位置特定部4により特定された同期位置を用いることで、どのように平行移動されたかも特定することができるため、オリジナル画像から入手画像に至る過程の平行移動、拡大縮小、回転、反転の全てが特定できたことになり、オリジナル画像と入手画像の位置的対応関係が明確になる。以上で明らかなように、この実施の形態1によれば、拡大縮小回転反転特定部3により特定された拡大率、回転角及び反転状況を用いて、位置合わせ信号を構成している複数の周波数スペクトルの中から、フーリエ変換部2から出力されたフーリエ係数に対応する周波数スペクトルを特定し、その周波数スペクトルの位相から、オリジナル画像に対する入手画像の平行移動量を算出して、その平行移動量からオリジナル画像と入手画像の同期位置を特定する同期位置特定部4を設け、微調整量算出部5が、同期位置特定部4により特定された同期位置と位置合わせ信号を用いて、拡大縮小回転反転特定部3により特定された拡大率及び回転角の微調整量を算出するように構成したので、オリジナル画像に対する入手画像の拡大率や回転角の誤差を補償して、画像の位置合わせ精度を高めることができる効果を奏する。実施の形態2.上記実施の形態1では、微調整量算出部5が、同期位置特定部4により特定された同期位置と位置合わせ信号を用いて、拡大縮小回転反転特定部3により特定された拡大率及び回転角の微調整量を算出するものを示したが、拡大率及び回転角の微調整量を用いて、拡大縮小回転反転特定部3により特定された拡大率及び回転角を補正するようにしてもよい。図8はこの発明の実施の形態2による画像位置合わせ装置を示す構成図であり、図において、図1と同一符号は同一または相当部分を示すので説明を省略する。拡大率補正部11は拡大縮小回転反転特定部3により特定された拡大率と微調整量算出部5により算出された拡大率の微調整量とを乗算することで、微調整済の拡大縮小情報を求め、微調整済の拡大縮小情報を同期位置特定部4に出力する処理を実施する。なお、拡大率補正部11は拡大縮小補正手段を構成している。回転角補正部12は拡大縮小回転反転特定部3により特定された回転角と微調整量算出部5により算出された回転角の微調整量を加算することで、微調整済の回転情報を求め、微調整済の回転情報を同期位置特定部4に出力する処理を実施する。なお、回転角補正部12は回転補正手段を構成している。図8の例では、画像位置合わせ装置の構成要素である部分画像切出部1、フーリエ変換部2、拡大縮小回転反転特定部3、同期位置特定部4、微調整量算出部5、拡大率補正部11及び回転角補正部12のそれぞれが専用のハードウェア例えば、CPUを実装している半導体集積回路、あるいは、ワンチップマイコンで構成されているものを想定しているが、画像位置合わせ装置がコンピュータで構成されていてもよい。画像位置合わせ装置がコンピュータで構成されている場合、部分画像切出部1、フーリエ変換部2、拡大縮小回転反転特定部3、同期位置特定部4、微調整量算出部5、拡大率補正部11及び回転角補正部12の処理内容を記述している画像位置合わせプログラムをコンピュータのメモリに格納し、当該コンピュータのCPUが当該メモリに格納されているプログラムを実行するようにすればよい。図9はこの発明の実施の形態2による画像位置合わせ装置の処理内容を示すフローチャートである。次に動作について説明する。拡大率補正部11及び回転角補正部12が実装されている点以外は、上記実施の形態1と同様であるため、主に拡大率補正部11及び回転角補正部12の処理内容を説明する。拡大率補正部11は、微調整量算出部5が拡大率の微調整量を算出すると、その微調整量を拡大縮小回転反転特定部3により特定された拡大率に乗算することで、微調整済の拡大縮小情報を求め拡大縮小回転反転特定部3により特定された拡大率を補正する、補正後の拡大率を示す微調整済の拡大縮小情報を同期位置特定部4に出力するステップST8,ST9。回転角補正部12は、微調整量算出部5が回転角の微調整量を算出すると、その微調整量を拡大縮小回転反転特定部3により特定された回転角に加算することで、微調整済の回転情報を求め拡大縮小回転反転特定部3により特定された回転角を補正する、補正後の回転角を示す微調整済の回転情報を同期位置特定部4に出力するステップST8,ST9。同期位置特定部4は、拡大率補正部11から微調整済の拡大縮小情報を受け、回転角補正部12から微調整済の回転情報を受けると、微調整済の拡大縮小情報が示す補正後の拡大率と、微調整済の回転情報が示す補正後の回転角と、拡大縮小回転反転特定部3から出力された反転情報とを用いて、再度、オリジナル画像に対する入手画像の平行移動量を算出して、その平行移動量からオリジナル画像と入手画像の同期位置を特定するステップST4,ST5。この実施の形態2では、同期位置特定部4がオリジナル画像と入手画像の同期位置を2回特定することになるが、1回目の特定処理で使用する拡大率及び回転角より、2回目の特定処理で使用する拡大率及び回転角の方が精度が高いため、1回目の特定処理より2回目の特定処理の方が、同期位置の特定精度が向上する。また、微調整量算出部5が、同期位置特定部4の2回目の特定処理で得られた同期位置を用いて、拡大率及び回転角の微調整量を算出することで、その微調整量の算出精度も高めることができる。ここでは、同期位置の特定処理と微調整量の算出処理を2回行うものを示したが、所定の回数だけ、同期位置の特定処理と微調整量の算出処理を繰り返し行うようにしてもよい。この繰り返しの回数は、予め指定するようにしてもよいし、拡大率及び回転角の微調整量が所定の閾値を下回るまで繰り返すようにしてもよい。上記実施の形態1,2の画像位置合わせ装置では、入手画像がデジタル非圧縮輝度画像である物を示したが、例えば、YCbCr形式やRGB形式に代表されるカラー画像であってもよい。このような場合でも、例えば、YやRなどの一つの成分を選択し、その1成分からなる画像データを非圧縮輝度画像と見立てて、実施の形態1,2の画像位置合わせ装置を適用することができる。ただし、事前に、どの画像成分に位置合わせ信号が埋め込まれているかは知っておかなければならない。例えば、RGB形式のG成分に位置合わせ信号が埋め込まれている場合、YCbCr形式の画像の位置合わせを行うには、まず、その画像をRGB形式に変換し、その上でG成分に対して、実施の形態1,2の画像位置合わせ装置を適用するG成分だけからなる画像を入手画像として、実施の形態1,2の画像位置合わせ装置を適用することで、画像の位置合わせが可能である。また、入手画像がJPEGやGIFなど、何らかの符号化方式圧縮方式によって符号化圧縮されている場合、復号伸張してから対象となる非圧縮成分に着目して、実施の形態1,2の画像位置合わせ装置を適用するようにすればよい。このように、実施の形態1,2の画像位置合わせ装置は、色の有無や色空間形式、画像の圧縮/非圧縮に依らずに適用することができる。上記実施の形態1,2では、フーリエ変換に基づく位置合わせ信号が埋め込まれた画像について、フーリエ変換に基づく位置合わせを行うものを示したが、フーリエ変換に基づく位置合わせに限らず、他の位置合わせ技術に基づく構成としてもよい。例えば、フーリエ変換以外にも離散コサイン変換のように、画像処理分野で多用される周波数領域への変換方式を用いてもよい。離散コサイン変換に基づく拡大縮小等の位置合わせを行う場合、拡大縮小回転反転特定部3では、位置合わせデータとして離散コサイン変換係数の周波数平面上での座標を用いる。また、フーリエ変換部2の代わりに、部分画像を離散コサイン変換する変換部を設け、その変換結果である係数に存在するピークの周波数平面上での座標と、位置合わせデータを構成している周波数スペクトルを照合することで、拡大縮小情報、回転情報及び反転情報を算出するようにする。また、平行移動の位置合わせについては、部分画像と位置合わせ信号のパターンとの相関演算による位置合わせも可能である。この場合、同期位置特定部4では、位置合わせデータとして位置合わせ信号のパターンそのものそのパターンを生成するための周波数スペクトル等の情報であってもよいを使用するもとし、部分画像と位置合わせ信号のパターンとを少しずつずらしながら相関値を算出し、相関値の絶対値が最大となるような部分画像と位置合わせ信号のパターンとの位置関係を明らかにすることで平行移動量を算出し、その平行移動量から同期位置を特定する。このように、拡大縮小回転反転特定部3と同期位置特定部4は、どのような技術による位置合わせ方式であっても、その技術に合わせて構成することで、適切に動作する画像位置合わせ装置とすることができる。即ち、拡大縮小回転反転特定部3と同期位置特定部4は、それぞれ拡大縮小等の位置合わせと平行移動の位置合わせが可能な任意の構成とすることができる。実施の形態3.図10はこの発明の実施の形態3による電子透かし検出装置を示す構成図である。図10において、画像位置合わせ装置21は図8の画像位置合わせ装置である。画像補正部22は画像位置合わせ装置21の拡大率補正部11による微調整済の拡大縮小情報、回転角補正部12による微調整済の回転情報、拡大縮小回転反転特定部3から出力された反転情報及び同期位置特定部4により特定された同期位置を示す座標を用いて、入手画像を補正する処理を実施する。なお、画像補正部22は画像補正手段を構成している。電子透かしデータ検出部23は画像補正部22により補正された入手画像から電子透かしデータを検出する処理を実施する。なお、電子透かしデータ検出部23は電子透かし検出手段を構成している。図10の例では、電子透かし検出装置の構成要素である画像位置合わせ装置21、画像補正部22及び電子透かしデータ検出部23のそれぞれが専用のハードウェア例えば、CPUを実装している半導体集積回路、あるいは、ワンチップマイコンで構成されているものを想定しているが、電子透かし検出装置がコンピュータで構成されていてもよい。電子透かし検出装置がコンピュータで構成されている場合、画像位置合わせ装置21、画像補正部22及び電子透かしデータ検出部23の処理内容を記述しているプログラムをコンピュータのメモリに格納し、当該コンピュータのCPUが当該メモリに格納されているプログラムを実行するようにすればよい。図11はこの発明の実施の形態3による電子透かし検出装置の処理内容を示すフローチャートである。次に動作について説明する。画像位置合わせ装置21は、上記実施の形態2と同様に、画像の位置合わせを実施し、微調整済の拡大縮小情報、微調整済の回転情報、反転情報及び同期位置を示す座標を出力するステップST11。画像補正部22は、画像位置合わせ装置21から微調整済の拡大縮小情報、微調整済の回転情報、反転情報及び同期位置を示す座標を受けると、これらの情報を用いて、入手画像を幾何変形することで入手画像を補正し、補正後の入手画像を出力するステップST22。画像補正部22が入手画像を補正する際に必要となる計算は、拡大率を指定した拡大縮小、回転角を指定した回転、左右反転、同期位置を示す座標から明らかとなる平行移動量分の平行移動であり、これらは、画像処理の公知技術で十分可能であるため詳細な説明を省略する。なお、平行移動について、横方向にMの倍数ピクセル、縦方向にMの倍数ピクセルの不確定性が残されるが、この実施の形態3の電子透かしデータ検出装置では、問題とならない。電子透かしデータ検出部23は、画像補正部22が入手画像を補正すると、補正後の入手画像から電子透かしデータを検出するステップST13。電子透かしデータは、公知の電子透かし技術によって埋め込まれており、位置補正を行った画像から検出することができる。電子透かしデータを埋め込む方法は、どのような方法であってもよいが、位置合わせ信号と干渉しないようにしなければならない。例えば、位置合わせデータが、フーリエ変換に基づく周波数スペクトルであり、画像の輝度成分に埋め込まれているとすれば、電子透かしデータを埋め込む方法は、フーリエ変換に基づく周波数スペクトルを用いた方法であるならば、色差成分に対して埋め込むなど埋め込む対象成分を別にするか、あるいは、輝度成分に埋め込むならば、位置合わせデータである周波数スペクトルに悪影響を与えない、別の方法でなければならない。周波数スペクトルによらない方法としては、特定少数画素について、その画素値を変更する方法や、特定の近隣画素集合に対して画素値を変更する方法など、様々な方法があり、電子透かしデータ検出部23で使用することができる。具体的な一例としては、例えば、特願2010-121354号に記載されている方法で電子透かしデータを埋め込むことができる。この方法は、周波数スペクトルを使用しないため、周波数スペクトルを使用する位置合わせ信号と干渉しない。また、当該文献におけるデータ埋込対象画像のWdataとHdataを、フーリエ変換サイズMと同じに値に設定して、電子透かしデータを埋め込んでおけば、縦横Mピクセル毎に、全く同じ埋込方法が適用されることになるため、画像位置合わせ装置21で残される横Mの倍数ピクセル、縦Mの倍数ピクセルの不確定性の影響を受けずに、電子透かしデータを検出することができる。勿論、Mピクセル毎に繰り返される埋込方法によって埋め込めば、他の電子透かし埋込方法によって埋め込まれている場合であっても、画像位置合わせ装置21が残す平行移動の不確定性の影響を受けずに、電子透かしデータを検出することができる。電子透かしデータ検出方法の詳細については、適用する電子透かし技術によるので詳細な説明を省略する。電子透かしデータ検出部23は、電子透かしデータの検出に成功すればステップST14、検出した電子透かしデータを出力して、一連の処理を終了する。一方、電子透かしデータの検出に失敗した場合ステップST14、入手画像の補正をやり直してステップST15,ST11、電子透かしデータの検出を試みることによって、電子透かしデータの検出成功率を高めることができる。電子透かしデータの検出成否の判断は、電子透かしデータに誤り検出データを付加することで実現することができる。誤り検出データを付加するものに限るものではなく、電子透かしデータ埋込技術特有の方法によって検出成否を判断するようにしてもよい。また、入手画像の補正のやり直しの際には、画像位置合わせ装置21から、前回出力された補正パラメータ以外の補正パラメータを出力させ、その出力された補正パラメータに基づいて、画像補正部22が入手画像を補正し、電子透かしデータ検出部23が当該入手画像から電子透かしデータを検出し直すことで行われる。この際、画像位置合わせ装置21は、最も確からしい補正パラメータ以外にも、スコアリングで順位付けするなどして、複数の補正パラメータの候補を保持しておき、スコアの高い補正パラメータから順に出力するなどの処理を行う。スコアリングせずに、可能性のある補正パラメータを順次出力する構成も可能であるが、スコアリングを行う構成とすることで、電子透かしデータ検出成功までの繰り返し回数を短縮できるというメリットがある。可能性のある補正パラメータを全て試した上で、電子透かしデータの検出に失敗した場合や、電子透かしデータの検出に成功する前に、繰り返し回数が予め設定された繰り返し回数の上限に達した場合は、最終的に電子透かしデータの検出に失敗したものとして、一連の処理を終了する。この実施の形態3の電子透かしデータ検出装置は、画像位置合わせ装置21を使用することで、高精度な画像位置合わせが可能になるため、電子透かしデータの検出成功率が高まる効果が得られる。なお、本願発明はその発明の範囲内において、各実施の形態の自由な組み合わせ、あるいは各実施の形態の任意の構成要素の変形、もしくは各実施の形態において任意の構成要素の省略が可能である。"