(57)【要約】  本発明は、変位前後の位相特異点から変位を検出する変位検出方法、及び、変位検出装置、変位検出プログラム、並びに、特徴点マッチング処理方法、特徴点マッチングプログラムに関し、変位前後の位相特異点から変位を検出する変位検出方法であって、位相特異点の位相構造から所定の要素を取得し、要素に基づいて位相特異点を特定し、変位後の位相特異点の位置を検出することにより、変位前後の位相特異点を確実に特定し、その変位を容易、かつ、確実に、さらには、回転変位を含む種々の変位を検出可能とすることを特徴とする。
2007088759.tif 000002
【特許請求の範囲】
【請求項１】
変位前後の位相特異点から変位を検出する変位検出方法であって、
前記位相特異点の位相構造から所定の要素を取得し、前記要素に基づいて前記位相特異点を特定する位相特異点特定手順と、
前記位相特異点特定手順で取得された要素に基づいて変位後の位相特異点の位置を検出する変位検出手順とを有することを特徴とする変位検出方法。
【請求項２】
前記位相特異点特定手順で取得される要素は、前記位相特異点近傍の実部と虚部の値を平面で補間し、その結果、実部がゼロとなる線分と虚部がゼロとなる線分とが交差する角度であることを特徴とする請求項１記載の変位検出方法。
【請求項３】
前記位相特異点特定手順で取得される要素は、前記位相特異点近傍の強度分布の形状の離心率であることを特徴とする請求項１記載の変位検出方法。
【請求項４】
前記位相特異点特定手順で取得される要素は、前記位相特異点近傍の強度の勾配であることを特徴とする請求項１記載の変位検出方法。
【請求項５】
前記位相特異点特定手順は、前記位相特異点近傍の実部と虚部の値を平面で補間し、その結果、実部がゼロとなる線分と虚部がゼロとなる線分とが交差する角度と、前記位相特異点近傍の強度分布の形状の離心率と、前記位相特異点近傍の強度の勾配である渦度とを含むことを特徴とする請求項１記載の変位検出方法。
【請求項６】
前記変位検出手順は、前記位相特異点特定手順で特定された位相特異点の変位前後の位置を前記要素に基づいて特定する変位前後位置検出手順と、
前記変位前後位置検出手順で検出された変位前後の位相特異点の位置に基づいて前記位相特異点の変位を測定する変位測定手順とを有することを特徴とする請求項１記載の変位検出方法。
【請求項７】
前記変位測定手順は、前記変位前後位置検出手順で検出された変位前後の位相特異点の位置の座標の差分の二乗平方根を位相特異点の変位量として算出することを特徴とする請求項６記載の変位検出方法。
【請求項８】
前記変位測定手順は、前記変位前後位置検出手順で検出された複数組の変位前後の位相特異点の位置を結ぶ線分の複数の垂直二等分線の距離の和が最小となる点を回転中心として求める回転中心取得手順と、
前記変位前後位置検出手順で検出された変位前後の位相特異点の位置と前記回転中心取得手順で取得した前記回転中心とからなる三角形の前記回転中心を挟む辺の角度を回転角として求める回転角度取得手順とを有することを特徴とする請求項６記載の変位検出方法。
【請求項９】
取得した画像の強度信号をフーリエ変換するフーリエ変換手順と、
前記フーリエ変換手順でフーリエ変換された信号の複素共役の対称性を崩すようなフィルタリング処理を行い、擬似位相情報を取得するフィルタリング処理手順と、
前記フィルタリング処理手順で取得した擬似位相情報に基づいて位相特異点を取得する位相特異点取得手順を有することを特徴とする請求項１記載の変位検出方法。
【請求項１０】
前記フィルタリング処理手順は、前記フーリエ変換手順でフーリエ変換された信号に対して中心に位相特異点を持つらせん状の位相特性を持ち、かつ、ドーナツ状の振幅特性を持つラゲールガウスビームフィルタ特性により前記擬似位相情報を取得することを特徴とする請求項９記載の変位検出方法。
【請求項１１】
前記画像の強度信号の平均値を前記画像の各強度信号から減算した画像を前記画像とする前処理手順を有することを特徴とする請求項９記載の変位検出方法。
【請求項１２】
前記フィルタリング処理手順は、前記フーリエ変換手順でフーリエ変換された信号に対して中心に位相特異点を持つらせん状の位相特性を持つフィルタ特性、又は、ヒルベルト変換により前記擬似位相情報を取得することを特徴とする請求項９記載の変位検出方法。
【請求項１３】
変位前の特徴点の特徴パラメータと変位後の特徴点の特徴パラメータとの差異に基づいて変位前の特徴点に対応する変位後の特徴点を探しだし、その特徴点の地理的な位置関係の類似性から変位前の特徴点に対応する変位後の特徴点のマッチングを行うマッチング手順を、更に、有することを特徴とする変位検出方法。
【請求項１４】
前記マッチング手順は、変位前の複数の位相特異点間の特徴パラメータの総和距離と変位後の複数の位相特異点間の特徴パラメータの総和距離との差分とが最小となる特徴点の組を変位前と変位後とで互いにマッチングする特徴点として特定することを特徴とする請求項１３記載の変位検出方法。
【請求項１５】
前記マッチング手順は、変位前の所定の位相特異点の特徴パラメータをＴi1、Ｔi2、変位後の所定の位相特異点の特徴パラメータをＰw(i1)、Ｐw(i2)、αを縮尺としたとき、
【数１】
2007088759.tif 000027
が最小となる位相特異点の組を変位前後でマッチングする位相特異点の組として抽出することを特徴とする請求項１３記載の変位検出方法。
【請求項１６】
変位前後の位相特異点から変位を検出する変位検出装置であって、
前記位相特異点の位相構造から所定の要素を取得し、前記要素に基づいて前記位相特異点を特定する位相特異点特定手段と、
前記位相特異点特定手段で取得された要素に基づいて変位後の位相特異点の位置を検出する変位検出手段とを有することを特徴とする変位検出装置。
【請求項１７】
前記位相特異点特定手段で取得される要素は、前記位相特異点近傍の実部と虚部の値を平面で補間し、その結果、実部がゼロとなる線分と虚部がゼロとなる線分とが交差する角度であることを特徴とする請求項１６記載の変位検出装置。
【請求項１８】
前記位相特異点特定手段で取得される要素は、前記位相特異点近傍の強度分布の形状の離心率であることを特徴とする請求項１６記載の変位検出装置。
【請求項１９】
前記位相特異点特定手段で取得される要素は、前記位相特異点近傍の強度の勾配であることを特徴とする請求項１６記載の変位検出装置。
【請求項２０】
前記位相特異点特定手段は、前記位相特異点近傍の実部と虚部の値を平面で補間し、その結果、実部がゼロとなる線分と虚部がゼロとなる線分とが交差する角度と、前記位相特異点近傍の強度分布の形状の離心率と、前記位相特異点近傍の強度の勾配である渦度とを含むことを特徴とする請求項１６記載の変位検出装置。
【請求項２１】
前記変位検出手段は、前記位相特異点特定手段で特定された位相特異点の変位前後の位置を前記要素に基づいて特定する変位前後位置検出手段と、
前記変位前後位置検出手段で検出された変位前後の位相特異点の位置に基づいて前記位相特異点の変位を測定する変位測定手段とを有することを特徴とする請求項１６記載の変位検出装置。
【請求項２２】
前記変位測定手段は、前記変位前後位置検出手段で検出された変位前後の位相特異点の位置の座標の差の二乗平方根を位相特異点の変位量として算出することを特徴とする請求項２１記載の変位検出装置。
【請求項２３】
前記変位測定手段は、前記変位前後位置検出手段で検出された複数組の変位前後の位相特異点の位置を結ぶ線分の複数の垂直二等分線の距離の和が最小となる点を回転中心として求める回転中心取得手段と、
前記変位前後位置検出手段で検出された変位前後の位相特異点の位置と前記回転中心取得手段で取得した前記回転中心とからなる三角形の前記回転中心を挟む辺の角度を回転角として求める回転角度取得手段とを有することを特徴とする請求項２１記載の変位検出装置。
【請求項２４】
取得した画像の強度信号をフーリエ変換するフーリエ変換手段と、
前記フーリエ変換手段でフーリエ変換された信号の複素共役の対称性を崩すようなフィルタリング処理を行い、擬似位相情報を取得するフィルタリング処理手段と、
前記フィルタリング処理手段で取得した擬似位相情報に基づいて位相特異点を取得する位相特異点取得手段を有することを特徴とする請求項１６記載の変位検出装置。
【請求項２５】
前記フィルタリング処理手段は、前記フーリエ変換手段でフーリエ変換された信号に対して中心に位相特異点を持つらせん状の位相特性を持ち、かつ、ドーナツ状の振幅特性を持つラゲールガウスビームフィルタ特性により前記擬似位相情報を取得することを特徴とする請求項２４記載の変位検出装置。
【請求項２６】
前記画像の強度信号の平均値を前記画像の各強度信号から減算した画像を前記画像とする前処理手段を有することを特徴とする請求項２４記載の変位検出装置。
【請求項２７】
前記フィルタリング処理手段は、前記フーリエ変換手段でフーリエ変換された信号に対して中心に位相特異点を持つらせん状の位相特性を持つフィルタ特性、又は、ヒルベルト変換により前記擬似位相情報を取得することを特徴とする請求項２４記載の変位検出装置。
【請求項２８】
変位前の特徴点の特徴パラメータと変位後の特徴点の特徴パラメータとの差異に基づいて変位前の特徴点に対応する変位後の特徴点を探しだし、その特徴点の地理的な位置関係の類似性から変位前の特徴点に対応する変位後の特徴点のマッチングを行うマッチング手段を、更に、有することを特徴とする変位検出装置。
【請求項２９】
前記マッチング手段は、変位前の複数の位相特異点間の特徴パラメータの総和距離と変位後の複数の位相特異点間の特徴パラメータの総和距離との差分とが最小となる特徴点の組を変位前と変位後とで互いにマッチングする特徴点として特定することを特徴とする請求項２８記載の変位検出装置。
【請求項３０】
前記マッチング手段は、変位前の所定の位相特異点の特徴パラメータをＴi1、Ｔi2、変位後の所定の位相特異点の特徴パラメータをＰw(i1)、Ｐw(i2)、αを縮尺としたとき、
【数２】
2007088759.tif 000028
が最小となる位相特異点の組を変位前後でマッチングする位相特異点の組として抽出することを特徴とする請求項２８記載の変位検出装置。
【請求項３１】
コンピュータに、
位相特異点の位相構造から所定の要素を取得し、前記要素に基づいて前記位相特異点を特定する位相特異点特定手順と、
前記位相特異点特定手順で取得された要素に基づいて変位後の位相特異点の位置を検出する変位検出手順とを実行させる変位検出プログラム。
【請求項３２】
前記位相特異点特定手順で取得される要素は、前記位相特異点近傍の実部と虚部の値を平面で補間し、その結果、実部がゼロとなる線分と虚部がゼロとなる線分とが交差する角度であることを特徴とする請求項３１記載の変位検出プログラム。
【請求項３３】
前記位相特異点特定手順で取得される要素は、前記位相特異点近傍の強度分布の形状の離心率であることを特徴とする請求項３１記載の変位検出プログラム。
【請求項３４】
前記位相特異点特定手順で取得される要素は、前記位相特異点近傍の強度の勾配であることを特徴とする請求項３１記載の変位検出プログラム。
【請求項３５】
前記位相特異点特定手順は、前記位相特異点近傍の実部と虚部の値を平面で補間し、その結果、実部がゼロとなる線分と虚部がゼロとなる線分とが交差する角度と、前記位相特異点近傍の強度分布の形状の離心率と、前記位相特異点近傍の強度の勾配である渦度とを含むことを特徴とする請求項３１記載の変位検出プログラム。
【請求項３６】
前記変位検出手順は、前記位相特異点特定手順で特定された位相特異点の変位前後の位置を前記要素に基づいて特定する変位前後位置検出手順と、
前記変位前後位置検出手順で検出された変位前後の位相特異点の位置に基づいて前記位相特異点の変位を測定する変位測定手順とを有することを特徴とする請求項３１記載の変位検出プログラム。
【請求項３７】
前記変位測定手順は、前記変位前後位置検出手順で検出された変位前後の位相特異点の位置座標の差分の二乗平方根を位相特異点の変位量として算出することを特徴とする請求項３５記載の変位検出プログラム。
【請求項３８】
前記変位測定手順は、前記変位前後位置検出手順で検出された複数組の変位前後の位相特異点の位置を結ぶ線分の複数の垂直二等分線の距離の和が最小となる点を回転中心として求める回転中心取得手順と、
前記変位前後位置検出手順で検出された変位前後の位相特異点の位置と前記回転中心取得手順で取得した前記回転中心とからなる三角形の前記回転中心を挟む辺の角度を回転角として求める回転角度取得手順とを有することを特徴とする請求項３６記載の変位検出プログラム。
【請求項３９】
取得した画像の強度信号をフーリエ変換するフーリエ変換手順と、
前記フーリエ変換手順でフーリエ変換された信号の複素共役の対称性を崩すようなフィルタリング処理を行い、擬似位相情報を取得するフィルタリング処理手順と、
前記フィルタリング処理手順で取得した擬似位相情報に基づいて位相特異点を取得する位相特異点取得手順を有することを特徴とする請求項３１記載の変位検出プログラム。
【請求項４０】
前記フィルタリング処理手順は、前記フーリエ変換手順でフーリエ変換された信号に対して中心に位相特異点を持つらせん状の位相特性を持ち、かつ、ドーナツ状の振幅特性を持つラゲールガウスビームフィルタ特性により前記擬似位相情報を取得することを特徴とする請求項３９記載の変位検出プログラム。
【請求項４１】
前記画像の強度信号の平均値を前記画像の各強度信号から減算した画像を前記画像とする前処理手順を有することを特徴とする請求項３９記載の変位検出プログラム。
【請求項４２】
前記フィルタリング処理手順は、前記フーリエ変換手順でフーリエ変換された信号に対して中心に位相特異点を持つらせん状の位相特性を持つフィルタ特性、又は、ヒルベルト変換により前記擬似位相情報を取得することを特徴とする請求項３９記載の変位検出プログラム。
【請求項４３】
変位前の位相特異点の特徴パラメータと変位後の位相特異点の特徴パラメータとの差に基づいて前記変位前に特定された位相特異点と変位後に特定された位相特異点とをマッチングを行うマッチング手順を、更に、有することを特徴とする請求項３１記載の変位検出プログラム。
【請求項４４】
前記マッチング手順は、変位前の複数の位相特異点間の特徴パラメータの総和距離と変位後の複数の位相特異点間の特徴パラメータの総和距離との差分とが最小となる特徴点の組を変位前と変位後とで互いにマッチングする特徴点として特定することを特徴とする請求項４３記載の変位検出プログラム。
【請求項４５】
前記マッチング手順は、変位前の所定の位相特異点の特徴パラメータをＴi1、Ｔi2、変位後の所定の位相特異点の特徴パラメータをＰw(i1)、Ｐw(i2)、αを縮尺としたとき、
【数３】
2007088759.tif 000029
が最小となる位相特異点の組を変位前後でマッチングする位相特異点の組として抽出することを特徴とする請求項４３記載の変位検出プログラム。
【請求項４６】
変位前に特定された特徴点と変位後に特定された特徴点とをマッチングを行うポイントマッチング処理方法であって、
変位前の特徴点の特徴パラメータと変位後の特徴点の特徴パラメータとの差異に基づいて変位前の特徴点に対応する変位後の特徴点を探しだし、その特徴点の地理的な位置関係の類似性から変位前の特徴点に対応する変位後の特徴点のマッチングを行うことを特徴とすることを特徴とする特徴点マッチング処理方法。
【請求項４７】
変位前の複数の特徴点間の特徴パラメータの総和距離と変位後の複数の特徴点間の特徴パラメータの総和距離との差分とが最小となる特徴点の組を変位前と変位後とで互いにマッチングする特徴点として特定することを特徴とする請求項３７記載の特徴点マッチング処理方法。
【請求項４８】
変位前の所定の特徴点の特徴パラメータをＴi1、Ｔi2、変位後の所定の特徴点の特徴パラメータをＰw(i1)、Ｐw(i2)、αを縮尺としたとき、
【数４】
2007088759.tif 000030
が最小となる特徴点の組を変位前後でマッチングする特徴点の組として抽出することを特徴とする請求項３８記載の特徴点マッチング処理方法。
【請求項４９】
前記特徴点は、画像の強度情報に基づいて取得される位相特異点であることを特徴とする請求項３７記載の特徴点マッチング処理方法。
【請求項５０】
前記位相特異点は、前記特徴パラメータとして位置情報、離心率、及び、渦度、並びに、実軸と虚軸とが交わる角度を有することを特徴とする請求項４０記載の特徴点マッチング処理方法。
【請求項５１】
コンピュータに、
変位前の特徴点の特徴パラメータと変位後の特徴点の特徴パラメータとの差異に基づいて変位前の特徴点に対応する変位後の特徴点を探しだし、その特徴点の地理的な位置関係の類似性から変位前の特徴点に対応する変位後の特徴点のマッチングを行うことを特徴とするコンピュータにより読み取り可能な特徴点マッチング処理プログラム。
【請求項５２】
変位前の複数の特徴点間の特徴パラメータの総和距離と変位後の複数の特徴点間の特徴パラメータの総和距離との差分とが最小となる特徴点の組を変位前と変位後とで互いにマッチングする特徴点として特定することを特徴とする請求項５１記載の特徴点マッチング処理プログラム。
【請求項５３】
変位前の所定の特徴点の特徴パラメータをＴi1、Ｔi2、変位後の所定の特徴点の特徴パラメータをＰw(i1)、Ｐw(i2)、縮尺をαとしたとき、
【数５】
2007088759.tif 000031
が最小となる特徴点の組を変位前後でマッチングする特徴点の組として抽出することを特徴とする請求項５２記載の特徴点マッチング処理プログラム。
【請求項５４】
前記特徴点は、画像の強度情報に基づいて取得される位相特異点であることを特徴とする請求項５１記載の特徴点マッチング処理プログラム。
【請求項５５】
前記位相特異点は、前記特徴パラメータとして位置情報、離心率、及び、渦度、並びに、実軸と虚軸とが交わる角度を有することを特徴とする請求項５４記載の特徴点マッチング処理プログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は変位検出方法、及び、変位検出装置、変位検出プログラム、並びに、特徴点マッチング処理方法、特徴点マッチングプログラムに係り、特に、変位前後の位相特異点から変位を検出する変位検出方法、及び、変位検出装置、変位検出プログラム、並びに、特徴点マッチング処理方法、特徴点マッチングプログラムに関する。
【背景技術】
【０００２】
近年、レーザースペックルパターンやランダムドットパターンなどのように空間的にランダムな構造をもつパターンやテクスチャを指標として物体の微小変位を非接触で計測する技術は非破壊検査や材料強度試験など産業応用上の重要な位置を占めている。パターンやテクスチャを用いて微小変位を検出する際の信号処理技術として相関法が知られている。【０００３】
図３８、図３９は相関法の動作を説明するため図を示す。図３８（Ａ）は場の変位、図３８（Ｂ）相関法の概念図、図３８（Ｃ）は本測定法の概念図を示している。また、図３９（Ａ）は物体Ａの一方向への変位、図３９（Ｂ）は物体Ａの回転変位を説明するための図を示している。【０００４】
相関法を用いて変位を検出する場合、図３８（Ａ）に矢印で示すように場の変位が位置により異なるときには図３８（Ｂ）に示すように各領域に分けて計算する必要があり、各点での相関値を決定するごとに２次元相関積分計算または２次元フーリエ変換・逆変換を行う必要があったため、１ピクセル以下での変位検出の分解能を上げようとしてピクセル内部での相関値の計算点数を増やすと計算時間が膨大になっていた。【０００５】
また、従来の相関法では、図３９（Ａ）に示すように一方向の変位に対しては相関をとることにより変位量を求めることができるが、図３９（Ｂ）に示すように回転に対しては適用できなかった。【０００６】
なお、位相情報を利用する相関関数の計算方法としては、位相限定相関法が知られている（特許文献１参照）。この相関法は、複素関数である合成スペクトルの振幅を一定化又は対数関数等により抑制し、位相情報のみからなる振幅限定複素合成スペクトルを作成してこれをフーリエ逆変換して相互相関関数を求めるものである。
【特許文献１】特許第３０３５６５４号
【発明の開示】
【発明が解決しようとする課題】
【０００７】
しかるに、従来の変位検出方法では、計算時間が膨大になるとともに、回転変位の検出が容易に行えないなどの問題点があった。【０００８】
本発明は上記の点に鑑みてなされたもので、位相特異点の構造を特徴付ける要素から所望の位相特異点を特定することにより変位前後の位相特異点を確実に特定し、その変位を容易、かつ、確実に、さらには、回転変位を含む種々の変位を検出できる変位検出方法、及び、変位検出装置、変位検出プログラム、並びに、特徴点マッチング処理方法、特徴点マッチングプログラムを提供することを目的とする。
【課題を解決するための手段】
【０００９】
本発明は、変位前後の位相特異点から変位を検出する変位検出方法であって、位相特異点の位相構造から所定の要素を取得し、要素に基づいて位相特異点を特定し、変位後の位相特異点の位置を検出することを特徴とする。【００１０】
このとき、取得される要素は、位相特異点近傍の実部と虚部の値を平面で補間し、その結果、実部がゼロとなる線分と虚部がゼロとなる線分とが交差する角度、位相特異点近傍の強度分布の形状の離心率、位相特異点近傍の強度の勾配であることを特徴とする。【００１１】
また、特定された位相特異点の変位前後の位置を要素に基づいて特定し、特定された変位前後の位相特異点の位置に基づいて位相特異点の変位を測定することを特徴とする。このとき、変位前後の位相特異点の位置の座標の差分の二乗平方根を位相特異点の変位量として算出することを特徴とする。【００１２】
複数組の変位前後の位相特異点の位置を結ぶ線分の複数の垂直二等分線の距離の和が最小となる点を回転中心として求め、変位前後の位相特異点の位置と回転中心とからなる三角形の回転中心を挟む辺の角度を回転角として求めることを特徴とする。【００１３】
さらに、本発明は、取得した画像の強度信号をフーリエ変換し、フーリエ変換された信号の複素共役の対称性を崩すようなフィルタリング処理を行い、擬似位相情報を取得し、取得した擬似位相情報に基づいて位相特異点を取得することを特徴とする。【００１４】
このとき、フーリエ変換された信号に対して中心に位相特異点を持つらせん状の位相特性を持ち、かつ、ドーナツ状の振幅特性を持つラゲールガウスビームフィルタ特性によりフィルタリングすることを特徴とする。また、画像の強度信号の平均値を画像の各強度信号から減算した画像を元画像とすることを特徴とする。【００１５】
さらに、フーリエ変換された信号に対して中心に位相特異点を持つらせん状の位相特性を持つフィルタ特性、又は、ヒルベルト変換によりフィルタリングを行い、擬似位相情報を取得することを特徴とする。【００１６】
また、本発明は、変位前に特定された特徴点と変位後に特定された特徴点とをマッチングを行うポイントマッチング処理方法であって、変位前の特徴点の特徴パラメータと変位後の特徴点の特徴パラメータとの差異に基づいて変位前の特徴点に対応する変位後の特徴点を探しだし、その特徴点の地理的な位置関係の類似性から変位前の特徴点に対応する変位後の特徴点のマッチングを行うことを特徴とする。【００１７】
また、このとき、本発明は、変位前の複数の特徴点間の特徴パラメータの総和距離と変位後の複数の特徴点間の特徴パラメータの総和距離との差分とが最小となる特徴点の組を変位前と変位後とで互いにマッチングする特徴点として特定することを特徴とする。【００１８】
さらに、本発明は、変位前の所定の特徴点の特徴パラメータをＴi1、Ｔi2、変位後の所定の特徴点の特徴パラメータをＰw(i1)、Ｐw(i2)、αを縮尺としたとき、【００１９】
【数１】
2007088759.tif 000003
が最小となる特徴点の組を変位前後でマッチングする特徴点の組として抽出することを特徴とする。【００２０】
なお、特徴点は、画像の強度情報に基づいて取得される位相特異点であることを特徴とし、位相特異点は特徴パラメータとして位置情報、離心率、及び、渦度、並びに、実軸と虚軸とが交わる角度を有することを特徴とする。
【発明の効果】
【００２１】
本発明によれば、位相特異点の位相構造から所定の要素を取得し、要素に基づいて位相特異点を特定することにより、変位前後の各位相特異点を一意に特定することができ、これによって、変位後の位相特異点を追跡できるなどの効果を得られる。
【図面の簡単な説明】
【００２２】
【図１】本発明の一実施例のシステム構成図である。
【図２】変位検出プログラムの処理フローチャートである。
【図３】位相特異点取得処理の処理フローチャートである。
【図４】位相特異点取得処理の動作説明図である。
【図５】位相特異点取得処理の動作説明図である。
【図６】位相特異点の位置をサブピクセルで決定する動作を説明するための図である。
【図７】位相特異点の位置をサブピクセルで決定する動作を説明するための図である。
【図８】位相特異点の位置をサブピクセルで決定する動作を説明するための図である。
【図９】位相特異点特定処理の処理フローチャートである。
【図１０】離心率ｅ、渦度ωを取得する動作を説明するための図である。
【図１１】離心率ｅ、渦度ωを取得する動作を説明するための図である。
【図１２】離心率ｅ、渦度ωを取得する動作を説明するための図である。
【図１３】全位相特異点の位置の差を求めｘ方向，ｙ方向の変位をヒストグラムに表した図である。
【図１４】回転変位の検出動作を説明するための図である。
【図１５】シミュレーション実験を説明するための図である。
【図１６】回転前後の位相特異点の数をヒストグラムで表した図である。
【図１７】回転前後の位相特異点の回転角θiをヒストグラムで表した図である。
【図１８】回転前後の角度Δθの差をヒストグラムで表した図である。
【図１９】回転前後の対応する位相特異点の離心率ｅをヒストグラムで表した図である。
【図２０】回転前後の対応する位相特異点の渦度ωの差をヒストグラムで表した図である。
【図２１】ポイントマッチング処理の動作説明図である。
【図２２】ポイントマッチング処理のフローチャートである。
【図２３】マッチングアルゴリズムの処理動作を説明するための図である。
【図２４】マッチングアルゴリズムの処理動作を説明するための図である。
【図２５】マッチングアルゴリズムの処理動作を説明するための図である。
【図２６】実験・検証に用いたＬＧフィルタの振幅特性図である。
【図２７】マッチング処理の動作説明図である。
【図２８】マッチング処理の動作説明図である。
【図２９】実験・検証結果を説明するための図である。
【図３０】実験・検証結果を説明するための図である。
【図３１】実験・検証結果を説明するための図である。
【図３２】実験・検証結果を説明するための図である。
【図３３】実験・検証結果を説明するための図である。
【図３４】実験・検証結果を説明するための図である。
【図３５】実験・検証結果を説明するための図である。
【図３６】点間総和距離の差によりマッチングを行う動作を説明するための図である。
【図３７】点間総和距離の差によりマッチングを行う動作を説明するための図である。
【図３８】相関法の動作を説明するため図である。
【図３９】相関法の動作を説明するため図である。
【符号の説明】
【００２３】
１０１  変位検出装置、１０２  測定対象
１１１  画像検出部、１１２  信号処理部
１２１  光源、１２２  撮像装置
１３１  フレームメモリ、１３２  処理部、１３３  フーリエ変換部
１３４  ＲＯＭ
１３５  ハードディスクドライブ、１３６  ディスプレイ、１３７  入力装置
【発明を実施するための最良の形態】
【００２４】
〔システム構成〕
図１は本発明の一実施例のシステム構成図を示す。【００２５】
本実施例の変位検出装置１０１は、画像検出部１１１及び信号処理部１１２から構成されている。【００２６】
画像検出部１１１は、光源１２１及び撮像装置１２２から構成されている。画像検出部１１１は、光源１２１から測定対象１０２に光を照射し、撮像装置１２２により測定対象１０２の表面画像を撮像する。【００２７】
画像検出部１１１で撮像された画像は、信号処理部１１２に供給される。信号処理部１１２は、例えば、コンピュータシステムであり、フレームメモリ１３１、処理部１３２、フーリエ変換部１３３、ＲＯＭ１３４、ハードディスクドライブ１３５、ディスプレイ１３６，入力装置１３７から構成されている。【００２８】
フレームメモリ１３１は、画像検出部１１１からの画像をフレーム毎に記憶する。処理部１３２は、ハードディスクドライブ１３５に予めインストールされた変位検出プログラムに基づいて測定対象１０２の微小変位を検出するための処理を実行する。【００２９】
フーリエ変換部１３３は、画像検出部１１１で撮像された画像に対してフーリエ変換を行う。メモリ１３４は、処理部１３２の作業用記憶領域として用いられる。ハードディスクドライブ１３５は、変位検出プログラムや画像検出部１１１で検出された画像データ、あるいは、変位検出プログラムでの途中結果、変位検出結果などのデータが記憶される。【００３０】
ディスプレイ１３６は、ＬＣＤ、ＣＲＴなどから構成されており、画像、変換画像、途中結果、変位検出結果などの処理部１３２での処理結果を表示する。入力装置１３７は、キーボード、マウスなどから構成されており、データ入力や各種コマンドの入力するために用いられる。【００３１】
本実施例の変位検出装置１０１は、インストールされた変位検出プログラムに基づいて測定対象１０２から取得した画像を解析して、その変位を検出する。【００３２】
〔変位検出プログラム〕
図２は変位検出プログラムの処理フローチャートを示す。【００３３】
処理部１３２は、まず、変位検出プログラムに従って、ステップＳ１−１で変位前後の対象画像をハードディスクドライブ１３５からメモリ１３４などに読み出す。【００３４】
次に処理部１３２は、ステップＳ１−２で、読み出された変位前後の対象画像に対して各々位相特異点を取得するための位相特異点取得処理を実行する。【００３５】
次に処理部１３２は、ステップＳ１−３で変位前後の対象画像から取得した位相特異点が持つ、実部と虚部とのゼロクロス角度、離心率、渦度などの要素により位相特異点を特定する位相特異点特定処理を実行する。【００３６】
次に処理部１３２は、ステップＳ１−４で、実部と虚部とのゼロクロス角度、離心率、渦度などの位相特異点に特有の要素により変位前後で対応する位相特異点を特定し、変位前後の位相特異点の位置から対象画像の変位を検出する変位検出処理を実行する。【００３７】
位相特異点によって、対象画像の変位を検出することで対象画像の変位を正確に特定することが可能となる。【００３８】
〔複素解析信号取得処理〕
次にステップＳ１−２の位相特異点取得処理について説明する。【００３９】
図３は位相特異点取得処理の処理フローチャート、図４、図５は位相特異点取得処理の動作説明図を示す。【００４０】
ステップＳ１−２の位相特異点取得処理は、取得した画像の強度信号をフーリエ変換し、フーリエ変換された信号の複素共役の対称性を崩すようなフィルタリング処理を行い、擬似位相情報を取得し、取得した擬似位相情報に基づいて位相特異点を取得する処理である。【００４１】
処理部１３２は、まず、ステップＳ２−１で前処理を行う。ステップＳ２−１の前処理は、対象画像の全体の強度情報の平均値を求め、求めた全体画像の強度情報の平均値を対象画像の強度情報から減算する処理を行う処理であり、この処理によって、画像データ強度情報中の直流成分が除去される。直流成分が除去されることにより、位相特異点が強調される。【００４２】
次に、処理部１３２は、ステップＳ２−２でフーリエ変換部１３３を制御して、強度情報に対して２次元高速フーリエ変換を行い、空間周波数スペクトルを得る。【００４３】
例えば、図４において、画像検出部１１１で得られた物体のテクスチャの強度情報２１１の関数２０１を【００４４】
【数２】
2007088759.tif 000004
とする。この強度情報を２次元フーリエ変換することにより、空間周波数スペクトル【００４５】
【数３】
2007088759.tif 000005
を求めることができる。なお、図４において２１２は空間周波数スペクトルの振幅分布、２１３は位相分布を示している。空間周波数スペクトルは、実部のみをもつ強度情報をフーリエ変換したものなので、周波数原点に対して複素共役の対称性、すなわちエルミート対称性をもつ。ＦＦＴを関数Ｆとすると、例えば、式２０２のように表せる。【００４６】
そこで、処理部１３２は、ステップＳ２−３で実部を基に虚部を作り出して複素解析信号を得るために、図４に２２１〜２２５で示すような分布を持つ、周波数領域上の複素共役の対称性を崩すようなフィルタリング処理を行う。フィルタリング処理では、例えば、ヒルベルト変換、スパイラル位相フィルタ、ＬＧ（ラゲールガウス）フィルタなどが用いられる。なお、図４において２２１はヒルベルト変換の振幅分布、２２２はスパイラル位相フィルタの位相分布、２２４はＬＧフィルタの振幅分布、２２５はＬＧフィルタの位相分布を示している。【００４７】
ヒルベルト変換の演算式は、例えば、図４、２０３で表すことができる。また、スパイラル位相フィルタの演算式は、例えば、図４、２０４で表すことができる。ＬＧフィルタの振幅透過分布２２４の演算式は、図４、２０５（Ａ）で表すことができる。ＬＧフィルタは、中心を対称に複素共役を位相上でπだけずらすことにより，ドーナツ状の振幅透過率分布が得られる。なお、ＬＧフィルタの位相分布２２５の演算式は、例えば、図４、２０５（Ａ）に角度φを導入することにより図４、２０５（Ｂ）に示すように表される。【００４８】
以上のように、演算式２０２に演算式２０３、２０４、２０５などを乗算することにより、本実施例で必要となる解析信号を得ることができる。また、エラーの原因となる高周波数成分を除去でき、直流成分を自動的、かつ、連続的にゼロにできる。さらに、ドーナツ状の強度分布により低周波数成分を取り除くことにより、位相特異点の数を増やすことができる。【００４９】
さらに、本実施例ではランダムなパターンの強度情報を用いて解析を行っており、測定対象１０２の構造などによってスペックという丸い斑点のようなパターンが多数存在することになる。このスペックルは、本実施例で指標としている位相特異点に対応関係がある。【００５０】
このスペックルのサイズによりフーリエ変換したときの周波数領域に位相特異点の情報を多く含む周波数が変わる。例えば、スペックルのサイズが小さいと位相特異点の情報を含む周波数成分は高くなり、逆にサイズが大きいと位相特異点の情報を含む周波数成分は低くなる。つまり、位相特異点の情報を高精度に求めるにはドーナツ状の強度分布をもつＬＧフィルタが最適であり、ドーナツ状の幅のサイズをスペックルのサイズに合わせて調節することにより、より位相特異点の情報を高精度に求めることができる。一方、ローパスフィルタでは、低周波数成分を強調し、位相特異点の情報をうまく取り出すことができない。【００５１】
以上のように、本実施例における解析信号の生成過程ではＬＧフィルタを用いることが最適となる。【００５２】
なお、図４において２１４はフィルタリング処理後の振幅分布、２１５はフィルタリング処理後の位相分布を示しており、その演算式は、例えば、図４、２０６で表すことができる。【００５３】
なお、２２１、２２３に示すようなヒルベルト変換またはスパイラル位相フィルタをかけても同種の解析信号を得ることはできるが、ヒルベルト変換またはスパイラル位相フィルタを使うためには予め直流成分を除去しておかなければならない。【００５４】
処理部１３２は、ステップＳ２−４でフィルタをかけた空間周波数スペクトルをフーリエ逆変換して複素解析信号【００５５】
【数４】
2007088759.tif 000006
を得る。なお、図４において２１６はフーリエ逆変換後の振幅分布、２１７はフーリエ逆変換後の位相特性を示しており、その演算式は、例えば、図４、２０７で表すことができる。【００５６】
図５（Ａ）はフーリエ逆変換後の振幅分布を拡大したものであり、図５（Ｂ）は図５（Ａ）の所定の部分を更に拡大したものである。【００５７】
この複素解析信号の擬似位相情報から図５（Ｂ）に示すような位相特異点Ｐを取得することができる。【００５８】
次に処理部１３２は、ステップＳ２−５で位相特異点の位置を１ピクセルより小さい分解能であるサブピクセルの分解能で決定する処理を行う。【００５９】
図６、図７、図８は位相特異点の位置をサブピクセルで決定する動作を説明するための図を示す。【００６０】
図６（Ａ）は擬似位相情報、図６（Ｂ）は位相特異点周辺画素の拡大図、図６（Ｃ）は実部、図６（Ｄ）は虚部を示している。また、図７（Ａ）は位相特異点近傍の擬似位相情報、図７（Ｂ）は位相特異点近傍の擬似位相情報の虚部、図７（Ｃ）は位相特異点近傍の擬似位相情報の実部、図７（Ｄ）は位相特異点近傍の擬似位相情報を直線補間したもの、図７（Ｅ）は位相特異点近傍の擬似位相情報の虚部を直線補間したもの、図７（Ｆ）は位相特異点近傍の擬似位相情報の実部を直線補間したものを示している。【００６１】
擬似位相情報は(-π，π)までの値をとり、図６、図７において黒から白になるにしたがって値が大きくなる。【００６２】
位相特異点Ｐは、その点を囲む閉経路上の位相勾配の積分が±２πの整数倍になり、また、複素数の実部と虚部が共にゼロとなるという性質を持っている。【００６３】
図６（Ｂ）、図７（Ｃ）に示すような互いに隣接しあう４画素を結ぶ正方形閉経路に対して位相勾配の積分が±２πの整数倍になるという性質を用いて位相特異点の存在している位置を決定していた。このため、１ピクセル以上の分解能をもって位相特異点の位置を決定することができない。【００６４】
本実施例では、複素数の実部と虚部が共にゼロとなるという性質を用いて１ピクセルの分解能より小さいサブピクセルの分解能で位相特異点の位置を決定するようにしている。本実施例では、図７（Ｂ）、図７（Ｃ）に示す特異点近傍の離散化された実部と虚部を平面で線形補間することにより図７（Ｅ）、図７（Ｆ）に示すように実部及び虚部を共に滑らかな構造として、実部がゼロとなる線分Ｌre0、虚部がゼロとなる線分Ｌim0を求める。図７（Ｄ）に示すように位相特異点の位置は実部がゼロとなる線分Ｌre0と虚部がゼロとなる線分Ｌim0の交わる点を位相特異点Ｐ0として設定することによりサブピクセルの分解能をもって位相特異点を決定することができる．
図８（Ａ）に示すように位相勾配の積分が±２πの整数倍になる性質を用いる方法では、１ピクセルの面積を持つ正方形の閉経路内のどこかに位相特異点が存在することはわかってもその位置を特定することができなかった。これに対して図８（Ｂ）に示すように複素数の実部と虚部が共にゼロとなるという性質を用いて実部と虚部の内挿と実部と虚部のゼロの交差点から位相特異点を求めることにより１ピクセルより小さいサブピクセルの分解能をもって位相特異点の位置を決定できている。【００６５】
〔位相特異点特定処理〕
図９は位相特異点特定処理の処理フローチャートを示す。【００６６】
ステップＳ１−３の位相特異点特定処理は、位相特異点の位相構造から所定の要素を取得し、要素に基づいて位相特異点を特定する処理である。【００６７】
処理部１３２は、ステップＳ３−１で位相特異点を特定する要素としてゼロクロス角を取得する。ゼロクロス角は、図７（Ｄ）、図８（Ｂ）に示す実部がゼロとなる線分Ｌre0と虚部がゼロとなる線分Ｌim0との交差する角度Δθである。【００６８】
また、処理部１３２は、ステップＳ３−２で位相特異点を特定する要素として離心率ｅを取得し、ステップＳ３−３で渦度ωを取得する。【００６９】
図１０、図１１、図１２は離心率ｅ、渦度ωを取得する動作を説明するための図を示す。図１０（Ａ）は擬似位相情報、図１０（Ｂ）は擬似位相情報の位相特異点近傍の振幅分布、図１０（Ｃ）は位相分布を示している。【００７０】
図１０（Ｂ）に示すように振幅分布は位相特異点近傍に近づくにつれて暗くなっていて、これは値がゼロに近づいていることを示している。【００７１】
また、図１１は特異点近傍の振幅分布の構造図を示している。図１１（Ａ）は振幅を等高線で表示したものを示しており、図１１（Ｂ）は三次元表示したものを示す。【００７２】
図１１（Ａ）に示すように位相特異点近傍の等振幅の分布は略楕円形状をしている。一般に位相特異点近傍の振幅分布は略楕円形をしており、その形状は位相特異点に応じて異なっており、位相特異点を特定する要素として用いることができる。この楕円形状の特徴を二次曲線の特徴を表す値である離心率によって示すことによって、位相特異点を特定することが可能となる。【００７３】
例えば、図１２に示すような楕円の場合、一般に数式で、【００７４】
【数５】
2007088759.tif 000007
で表せる。図１２に示す楕円の離心率ｅはａ：長軸半径，ｂ：短軸半径，ｆ：焦点距離とすると、離心率ｅは【００７５】
【数６】
2007088759.tif 000008
で表させる。【００７６】
図１１（Ｂ）に示される振幅を三次元で表示した図を見ると位相特異点近傍の形状は楕円形の円錐のような構造をしていることが分かる。したがって、本実施例では、各位相特異点を特定する要素の一つとして位相特異点近傍の渦度を求める。渦度とは微小部分の局所的な回転の強さと回転の方向を表す量である。図１１（Ａ）では振幅の等高線の間隔が狭いほど渦度は大きくなり、図１１（Ｂ）では三次元表示した位相特異点近傍の勾配（変化の度合い）が急になるほど渦度は大きくなる。すなわち、渦度は位相特異点近傍の振幅の大きさに依存する。【００７７】
渦度ωは次式で表される。【００７８】
【数７】
2007088759.tif 000009
以上のように位相特異点は位置（ｘ、ｙ），ゼロクロス角Δθ，離心率e，渦度ωの要素を持っていて、これらの要素を用いて各位相特異点を一意に特定することができる。このとき、複数の要素を用いることにより確実に位相特異点を特定することができる。【００７９】
〔直線変位検出〕
次にステップＳ１−４の位相特異点変位検出処理について説明する。【００８０】
位相変位検出処理では、変位前後の位相特異点を上記の特定方法により特定することにより変位量を検出する。【００８１】
このとき、まず、変位前後で位相特異点をそれぞれ一意に特定するためのパラメータとして前述のように求めたゼロクロス角θ，離心率ｅ，渦度ωを用いて位相特異点を特定する。その方法は変位前の各パラメータをθ1、ｅ1、ω1、変位後の各パラメータをθ2、ｅ2、ω2とする。【００８２】
上記パラメータ（Δθ1、ｅ1、ω1）、（Δθ2、ｅ2、ω2）のそれぞれの差（Δθ2−Δθ1、ｅ2−ｅ1、ω2−ω1）の絶対値がある値より小さい時に同一の位相特異点であるとみなす。各位相特異点をそれぞれ一意に特定できるため、変位前後で各位相特異点がどこに移動したのか分かる。このとき、変位前の位相特異点の位置（ｘ1、ｙ1）と変位後の位相特異点の位置（ｘ2、ｙ2）の差から変位量を求めることができる。【００８３】
実際の変位量ΔＬは、Δｘ＝（ｘ2−ｘ1）、Δｙ＝（ｙ2−ｙ1）とすると、【００８４】
【数８】
2007088759.tif 000010
で表すことができる。【００８５】
なお、図１３は、全位相特異点の位置の差を求めｘ方向，ｙ方向の変位をヒストグラムに表した図を示す。図１３（Ａ）はｘ方向の変位、図１３（Ｂ）はｙ方向の変位を示している。【００８６】
図１３に示すように全位相特異点の位置がｘ方向に略平行に変位したことがわかる。【００８７】
〔回転変位検出〕
また、本実施例では、各位相特異点を回転変位前後で特定できるため、回転の検出も可能である。【００８８】
図１４は回転変位の検出動作を説明するための図を示す。【００８９】
図１４（Ａ）に示すように回転前の２組の位相特異点Ｐ11、Ｐ12のｘｙ座標をそれぞれ
（ｘ1、ｙ1）、（ｘ2、ｙ2）
とし，回転後の対応する位相特異点Ｐ21、Ｐ22のｘｙ座標をそれぞれ
（ｘ'1、ｙ'1）、（ｘ'2、ｙ'2）
とすると、回転の中心
（ｘc、ｙc）
は図１４（Ａ）に示すように回転前後の位相特異点Ｐ11、Ｐ12と位相特異点Ｐ21、Ｐ22の垂直二等分線【００９０】
【数９】
2007088759.tif 000011
の交わる点となる。【００９１】
一般的に回転の中心を求めるには２組の回転前後の位相特異点が分かれば十分であるが、実際は多数の位相特異点が存在するので回転前後の各位相特異点から得られる垂直二等分線が図１４（Ｂ）に示すに必ずしも１点では交わらない。このような場合には、回転の中心（ｘc、ｙc）と各垂直二等分線との距離の和が最小となるような回転の中心（ｘc、ｙc）を最小二乗法により求めるようにする。【００９２】
以上により回転中心を決定できる。【００９３】
次に図１４（Ｃ）に示すように回転前後の各位相特異点Ｐ11、Ｐ12と回転中心（ｘc、ｙc）からなる三角形に対して余弦定理を用いることによって、回転角θを求めることができる。【００９４】
以上により、本実施例によれば、回転変位の検出を行うことが可能となる。【００９５】
〔シミュレーション実験結果〕
実験で得た強度情報をコンピュータ内部で９０度だけ回転させて，回転させる前後の強度情報を用いてシミュレーション実験を行った。【００９６】
図１５はシミュレーション実験を説明するための図を示す。なお、図１５は、画素数が１０２４×１０２４で、回転させる前後の強度情報からＬＧフィルタを用いて２つの解析信号を得たときの振幅分布と位相分布を表示したものである。図１５（Ａ）は回転変位前の擬似位相情報、図１５（Ｂ）は回転変位前の擬似位相情報の振幅分布、図１５（Ｃ）は回転変位前の擬似位相情報の位相分布、図１５（Ｄ）は回転変位後の擬似位相情報、図１５（Ｅ）は回転変位後の擬似位相情報の振幅分布、図１５（Ｆ）は回転変位後の擬似位相情報の位相分布を示している。【００９７】
この回転前後の位相特異点の数は両方とも正の位相特異点が322個，負の位相区異点321個の合計643個で同数の位相特異点が存在していた。図１５（Ｂ）及び図１５（Ｅ）に示される振幅分布からも明らかなように回転前の三角で囲った２つのリング状の振幅が回転後には回転の中心（全画素の中心）に対して９０°回転している。【００９８】
次に位相特異点を特定するための要素であるRe=0とIm=0の直線の交わる角度：Δθ，離心率（eccentricity）：ｅ，渦度（vorticity）：ωを回転前後の全位相特異点に対して求める。回転前後で対応する位相特異点を探す条件は以下の条件を用いた。【００９９】
条件１：Re=0とIm=0の直線の交わる角度Δθの差が±0.5度以下
条件２：回転前後の離心率ｅの差が±0.005以下
条件３：回転前後の渦度ωの差が±0.05以下
回転前の位相特異点に対応する回転後の位相特異点が１対１の対応になっているのかを調べた。【０１００】
図１６は回転前の位相特異点に対応する回転後の位相特異点の数をヒストグラムで表した図を示す。【０１０１】
回転前後でほぼ全ての位相特異点が一意に特定できた。
シミュレーションで得られた測定結果の回転の中心は
（ｘc、ｙc）＝（511.5、511.5）
となった。これは設定した回転の中心
（ｘc、ｙc）＝（511.5、511.5）
と一致する。【０１０２】
図１７は回転前の位相特異点に対応する回転後の位相特異点と回転の中心から求めたそれぞれの回転角θiをヒストグラムで表した図を示す。【０１０３】
シミュレーション結果である回転角のヒストグラムは±0.001度の精度で全位相特異点が９０度回転していることを示している。シミュレーション実験では、測定対象１０２を９０度回転させており、図１７に示す測定結果も９０度回転していることを示しているので、この測定結果から本発明が回転に対しても有効に計測できることを確認できる。【０１０４】
図１８は回転前後の対応する位相特異点の角度Δθの差をヒストグラムで表した図、図１９は回転前後の対応する位相特異点の離心率ｅをヒストグラムで表した図、図２０は回転前後の対応する位相特異点の渦度ωの差をヒストグラムで表した図を示す。【０１０５】
図１８、図１９、図２０に示すようにゼロ付近に固まっており、誤差が少なく、位相特異点を確実に特定できていることがわかる。【０１０６】
なお、上記実施例では、位相特異点を用いて回転変位を検出する処理について説明したが縮尺、移動、回転を簡単な処理により検出することも可能である。【０１０７】
〔ポイントマッチング処理〕
ポイントマッチング処理は、任意の複数点の位置関係と、その点の空間構造の類似性を用いて変位前のパターンから変位後のパターンを検出する処理である。【０１０８】
すなわち、位相構造が類似する点を候補点として、その候補点を含んだパターン内の点間の距離が最小となる点を変位後のパターンとして決定するというものである。【０１０９】
以下にポイントマッチング処理の手順を説明する。【０１１０】
まず、変位前後の２つの強度分布画像に対して各々前述の位相特異点を抽出するためのフィルタ演算を行うことにより位相特異点を抽出する。抽出した位相特異点近傍にある実部及び虚部を線形補間することによって、サブピクセルの精度で全位相特異点の位置を特定する。更に、位置が特定された各位相特異点に対して離心率ｅ、実軸と虚軸との角Δθ、渦度ωを求め、ラベリングを行う。【０１１１】
図２１は、ポイントマッチング処理の動作説明図を示す。図２１（Ａ）は変位前の位相特異点の状態、図２１（Ｂ）は変位後の位相特異点の状態を示している。【０１１２】
まず、図２１（Ａ）に示す位相特異点から黒丸で示す位相特異点を選択する。選択する位相特異点はその特徴パラメータである離心率ｅ、実軸と虚軸との角Δθ、渦度ωが他の多くの位相特異点とは明らかにことなるものを抽出することが望ましい。これによって、変位後の位相特異点を特定しやすくなる。【０１１３】
図２１（Ａ）に示す変位前の画像から選択した位相特異点の特徴パラメータと位置と図２１（Ｂ）に示す変位後の画像から選択した位相特異点の特徴パラメータと位置とから図２１（Ｂ）に黒丸で示すように図２１（Ａ）に黒丸で示す位相特異点に対応する位相特異点を検出する。【０１１４】
位相特異点の特徴パラメータ、離心率e、渦度ω、実軸と虚軸との角Δθは、変位前後で同じであるか、または、ほとんど変化しない。このため、変位前後の位相特異点を正確に特定することによって移動、回転、縮尺などを正確に測定することができる。【０１１５】
図２２はポイントマッチング処理のフローチャートを示す。【０１１６】
まず、処理部１３２は、ステップＳ４−１で、図３とともに説明した位相特異点取得処理によって変位前後の位相特異点を取得し、ラベリングする。【０１１７】
次に、処理部１３２は、ステップＳ４−２で変位前の各位相特異点にほぼ同じ、あるいは、所定の範囲内に近似した特徴パラメータを持つ位相特異点を変位後の位相特異点から探索する。【０１１８】
次に処理部１３２は、ステップＳ４−３で変位前の位相特異点から所定の２点を選択し、それを組み合わせて組み合わせＭとする。次に、処理部１３２は、ステップＳ４−４で変数ｋに「３」が代入される。【０１１９】
次に、処理部１３２は、ステップＳ４−５で変数ｋが定数ｎより小さいか否かを判定する。なお、定数ｎは、探索する位相特異点の数である。【０１２０】
次に処理部１３２は、ステップＳ４−６で選択された複数の位相特異点から次の探索対象となる位相特異点を抽出する。【０１２１】
次に処理部１３２は、ステップＳ４−７で新たに抽出された位相特異点に対して特徴パラメータが近似する変位後の位相特異点を探索し、組み合わせＭとして追加する。【０１２２】
次に処理部１３２は、ステップＳ４−８で組み合わせＭの位相特異点の特徴パラメータに対して下記の式（１）によって計算を行い、閾値以下のものを変位前の位相特異点と変位後の位相特異点との有力な組み合わせの候補として残す。【０１２３】
【数１０】
2007088759.tif 000012
式（１）においてｔｉとｐw(i)は変位前と変位後の点を表している。なお、ｔｉはｐw(i)と対応している。
Ｔ＝｛ｔ１，ｔ２，ｔ３，・・・,ｔｎ｝
は変位前の位相特異点の組み合わせを表している。そして、
Ｐ＝｛ｐ１，ｐ２，ｐ３，・・・,ｐｎ｝
は変位後の位相特異点の組み合わせを表している。【０１２４】
また、｜ｔi1ｔi2｜は、点ｔi1と点ｔi2との間の距離，αは縮尺のスケール(度合い)を表している。さらに、変数ｗはＴ→Ｐの対応関係を表している。【０１２５】
なお、上記式（１）の関数をコスト関数（cost function）と呼ぶ。【０１２６】
処理部１３２は、ステップＳ４−９で変数ｋに（ｋ＋１）を代入し、ステップＳ４−５に戻る。処理部１３２は、上記ステップＳ４−６〜Ｓ４−９をステップＳ４−５で変数ｋが定数ｎより大きくなるまで繰り返す。【０１２７】
次に処理部１３２は、ステップＳ４−１０で残った組み合わせＭの位相特異点の特徴パラメータ間の距離が最小となるものをマッチングパターンとして抽出する。【０１２８】
次に、処理部１３２は、ステップＳ４−１１で抽出した特徴パラメータ間の距離が最小となる位相特異点の変換パラメータ、例えば、平行移動(shift),回転(rotation)、縮尺(scale)を計算する。【０１２９】
上記ポイントマッチング処理により、目標はパターン内の複数の点の位置関係において変位前後で最も良く一致するものを見つけ出すことができる。【０１３０】
次に、上記マッチングアルゴリズムの処理動作を、図面を用いてさらに詳細に説明する。【０１３１】
図２３、図２４、図２５はマッチングアルゴリズムの処理動作を説明するための図を示す。【０１３２】
図２３（Ａ）は変位前の位相特異点の中から複数の位相特異点の組を選んだものであり、これをＴ空間とし、ｔi（ｉ＝１、...,ｎ）、Ｔi（Ｘi,Ｙi,ｅi,ωi,Δθi）で表す。図２３（Ｂ）は変位後の位相特異点でＰ空間とし、ｐj（ｊ＝１、...,ｋ）、Ｐj（Ｘj,Ｙj,ｅj,ωj,Δθj）で表される。する。【０１３３】
まず、図２３（Ａ）に黒丸で示すように、変位前の位相特異点から任意の個数の位相特異点を選ぶ。【０１３４】
変位前のパターン内のそれぞれの位相特異点ｔiの特徴パラメータである離心率ｅ、ゼロクロス角θ、渦度ωと変位後の位相特異点ｐjの特徴パラメータである離心率ｅ、ゼロクロス角θ、渦度ωとの差を閾値と比較する。それは、例えば、次式で表せる。【０１３５】
【数１１】
2007088759.tif 000013
なお、式（２）で２番目の分母の｜ωTi＋ωPj｜は正規化するために付け加えたものである。実際のプログラム内ではこの形で使われている。また、ｅthは離心率の閾値、θthはゼロクロス角の閾値、ωthは渦度の閾値を示している。【０１３６】
次にＴ空間の位相特異点の組とＰ空間の位相特異点の組とのマッチングを行う。【０１３７】
図２４において、実線で示す丸印は正の位相特異点を示しており、破線で示した丸印は負の位相特異点を点線で示している。【０１３８】
まず、図２４（Ａ）に示すＴ空間の２個の位相特異点の極性、正か負か符号と、その位置関係とに基づいて図２４（Ｂ）に破線の直線で示すようにＰ空間に存在しうる可能性のある２個の位相特異点の組を探索する。【０１３９】
このマッチングでは移動(shift)，回転(rotation)，だけでなく縮尺(scale)も考慮しているので、Ｔ空間での初めの２つの組をＰ空間の中から２つを選ぶ組み合わせは、この例では正の位相特異点と負の位相特異点の組み合わせ全てであり、正が７個、負が６個なので４２通り存在する。【０１４０】
なお、ここで、正の位相特異点と負の位相特異点というのは特徴パラメータの渦度(vorticity)の値が正のものを正の位相特異点、また、負のものを"負の位相特異点としている。【０１４１】
この得られた位相特異点の組み合わせに対して，位相特異点の特徴パラメータの条件に合わない組み合わせを切り捨てる。このとき、各特徴パラメータは実験環境により変化しやすいことを考え、条件を緩くとることが好ましい。ここで言う条件とは、式（２）の閾値ｅth、θth、ωthのことである。【０１４２】
次に、ステップＳ４−３で空間Ｔ内の２つの位相特異点ｔi1、ｔi2を選ぶ。この位相特異点ｔi1、ｔi2は少なくとも１つのマッチング候補点を持っている。位相特異点ｔi1はｎ1個のマッチング候補の点を、位相特異点ｔi2はｎ2個のマッチング候補点を持つ。したがって、変位後の空間Ｐ内の候補点の組み合わせはｎ1ｎ２通りあることがわかる。それをテンポラルマッチング（temporal matching）と呼ぶことにする。【０１４３】
次に、図２５のようにＴ空間の位相特異点を１つ増やし、Ｐ空間の中から３点目として可能性のある点を探す。これにより３点目を含む組み合わせが得られる。なお、特徴パラメータを用いて明らかに間違った点は除外される。このような操作をＴ空間で選択される位相特異点がｎ個になるまで繰り返す。【０１４４】
つまり、新しいマッチングポイントを加えつつ、上記テンポラルマッチングを最終的にすべての位相特異点のマッチング点が決定されるまで繰り返す。【０１４５】
この操作を行うためには、新しく追加されていく点の地理的な位置関係が変位前と変位後で類似性があるかどうか調べる必要がある。【０１４６】
式（１）では、テンポラルマッチングの類似性を調べると同時に距離の閾値は間違って検出されたマッチングポイントを除外する操作も行っていることになる。【０１４７】
最終的にで得られた位相特異点の組み合わせＭの中から式（１）の評価式によりＴ空間の全位相特異点間の距離と対応するＰ空間の位相特異点間の空間的な距離及び特徴点の距離の差の総和が最小になるものを選び出してマッチングを行っている。【０１４８】
さらに、図面を用いて説明を続ける。【０１４９】
図２６、図２７はマッチング処理の動作説明図を示す。図２６（Ａ）〜（Ｃ）、図２７（Ａ）〜（Ｃ）は変位前の位相特異点の分布、すなわち、Ｔ空間における位相特異点の分布、図２６（Ｄ）〜（Ｆ）、図２７（Ｄ）〜（Ｆ）は変位後の位相特異点の分布、すなわち、Ｐ空間における位相特異点の分布を示す。【０１５０】
まず、図２６（Ａ）、（Ｄ）に示すようにＴ空間及びＰ空間のそれぞれの位相特異点の特徴パラメータを計算し、注目ポイントを選択する。図２６（Ａ）では、ｔ1、ｔ2、ｔ4の組み合わせを選択した。【０１５１】
次に、図２６（Ｂ）、（Ｅ）に示すようにＴ空間内のそれぞれの点の特徴パラメータと似たＰ空間の候補点Ｐ1〜Ｐ4を探索する。【０１５２】
次に、図２６（Ｃ）、（Ｅ）に示すようにＴ空間内の２つの位相特異点ｔ1、ｔ2を選択して、Ｐ空間内の対応する２つの位相特異点ｐ1、ｐ2を選択し、組み合わせＭとする。【０１５３】
次に、図２７（Ａ）、（Ｄ）に示すように既に選択した位相特異点以外で候補となる位相特異点ｔ3、ｔ4及びそれに対応するＰ空間内の位相特異点ｐ3、ｐ4を探索する。【０１５４】
次に、図２７（Ｂ）、（Ｆ）に示すようにＴ空間の位相特異点ｔ1、ｔ2、ｔ4とＰ空間の位相特異点ｐ1、ｐ2、ｐ4の特徴パラメータを式（１）に代入して、式（１）の値を算出する。同様に、Ｔ空間の注目位相特異点ｔ1、ｔ2、ｔ4とＰ空間の位相特異点ｐ1、ｐ2、ｐ3の特徴パラメータを式（１）に代入して、式（１）の値を算出する。【０１５５】
算出結果のうち閾値より小さい位相特異点の組を抽出する。例えば、ここでは、組み合わせＭ１（ｐ1、ｐ2、ｐ4）及びＭ2（ｐ1、ｐ2、ｐ3）が抽出されたとする。抽出された位相特異点の組Ｍ1、Ｍ2をマッチングした組み合わせとして抽出する。【０１５６】
次に図２７（Ｃ）、（Ｇ）に示すようにマッチングした組Ｍ1、Ｍ2のうち最もコスト、すなわち、式（１）の値が低い組Ｍ1をマッチングした組として抽出する。【０１５７】
以上は、説明を簡単にするために注目する位相特異点の数を３点とした場合について説明したが、３点以上の位相特異点に注目して同様な処理を行うことも可能である。【０１５８】
次に、ステップＳ４−１１で算出される変換パラメータについて説明する。すなわち、位相特異点の移動が平行移動、回転、縮尺の３つのパラメータに分解できることについて説明する。【０１５９】
一般に、流体力学の観点からどんな複雑な動きもそれを平行移動(transition)，回転(rotation)，縮尺(scale)で表すことができることが知られており、変位前後座標の回転関係式は次式で表すことができる。【０１６０】
【数１２】
2007088759.tif 000014
ここで、（ｘi、ｙi）は変位前の位相特異点の座標、（ｘi'、ｙi'）は変位後の位相特異点の座標、αは縮尺、θは回転角である。【０１６１】
【数１３】
2007088759.tif 000015
は変位前の位相特異点の重心、【０１６２】
【数１４】
2007088759.tif 000016
は変位後の位相特異点の重心を表している。
まず、はじめに、マッチングを表すのに適した関数について説明する。【０１６３】
この関数はそれぞれ変位前後で対応している位相特異点の差を足し合わせた形で書ける。【０１６４】
【数１５】
2007088759.tif 000017
ここで、Ｃ＝αcosθ、Ｓ＝αsinθを表している。Ｅの値の最小値を見つけるためにはそれぞれのパラメータをまず計算しなければならない。例えば、【０１６５】
【数１６】
2007088759.tif 000018
などである。【０１６６】
ここで、【０１６７】
【数１７】
2007088759.tif 000019
とおくと、【０１６８】
【数１８】
2007088759.tif 000020
から次式が得られる。【０１６９】
【数１９】
2007088759.tif 000021
これを解くと次のようになる。【０１７０】
【数２０】
2007088759.tif 000022
同様にして【０１７１】
【数２１】
2007088759.tif 000023
から次式が得られる。【０１７２】
【数２２】
2007088759.tif 000024
これを解くと次式が得られる。【０１７３】
【数２３】
2007088759.tif 000025
ＳとＣの定義より以下のパラメータが得られる。【０１７４】
【数２４】
2007088759.tif 000026
このように、位相特異点の移動は、平行移動Δｘ、Δｙ、回転θ、縮尺αの３つのパラメータで表せることがわかる。この３つのパラメータを変換パラメータと呼び、ステップＳ４−１１で算出される。【０１７５】
〔実験・検証〕
次に、上記マッチング処理方法の実験・検証結果について説明する。【０１７６】
図２８は実験・検証に用いられたＬＧフィルタの振幅特性図、図２９〜図３５は実験・検証結果を説明するための図を示す。図２９は実験画像を示す図であり、図２９（Ａ）は変位前の画像、図２９（Ｂ）は変位後の画像を示している。図３０は図２９に示す画像から抽出された位相特異点を示す画像であり、図２９（Ａ）は変位前の位相特異点の分布、図２９（Ｂ）は変位後の位相特異点の分布を示している。図３１は選択した位相特異点を拡大した図であり、図３１（Ａ）は変位前、図３１（Ｂ）は変位後の拡大図を示している。【０１７７】
まず、今回の実験では、対象としてみどりふぐを使用し、高速ＣＣＤカメラを用いて秒間３０枚、連続撮影した画像を用いてマッチング処理を検証している。【０１７８】
撮影した画像に対してＬＧフィルタ演算を行い、位相特異点を得る。このとき、ＬＧフィルタのフィルタ幅は、ω2＝３０としている。なお、ＬＧフィルタは、spiral phase filterにバンドパスフィルタの特性を加えたものであり、図２８に示すような振幅特性とした。なお、ＬＧフィルタの振幅幅ωを広くとれば、出てくる位相特異点の数は多くなり、逆に狭くすれば少なくなる。【０１７９】
図２９に示す変位前後の画像に対して今回提案したマッチングアルゴリズムを適用した例について説明する。図２９に示す画像の変位前後の位相特異点の分布は図３０に示すようになる。図３０（Ａ）に示す変位前の位相特異点から黒丸で示される６つの位相特異点を選択してマッチングを行った。【０１８０】
図３０（Ｂ）において、黒丸で示す位相特異点は変位前の黒丸で示す位相特異点に対応しており、実際にマッチングアルゴリズムにより検出されたパターンである。【０１８１】
図３１は変位前に選んだ位相特異点の組み合わせと、ポイントマッチングアルゴリズムによって検出された変位後の位相特異点の組み合わせを拡大図示したものである。このように変位前，変位後のパターンを見比べると移動，回転，縮尺を組み合わせた変形が起こっているがマッチングアルゴリズムによって検証された。【０１８２】
異なる２枚の画像に対して異なったパターンに着目して、本実施例のポイントマッチング処理を採用した場合の結果についても検証した。【０１８３】
図３２は異なる２枚の画像に対して異なったパターンに着目して、本実施例のポイントマッチング処理を採用した場合の変位前後の位相特異点の分布を示す図であり、図３２（Ａ）は変位前の分布、図３２（Ｂ）は変位後の分布を示している。図３３は図３２に示した位相特異点の変位を示す図であり、図３３（Ｂ）は選択された位相特異点を拡大した図である。【０１８４】
この検証例においても図３２（Ａ）に示すように変位前の４つの位相特異点が変位後、図３２（Ｂ）に示す位置に移動し、図３３に示すように注目した４つの位相特異点が変位していることが検証できた。【０１８５】
さらに同じ作業をふぐの連続撮影画像１００枚に対して行った。【０１８６】
図３４は連続した画像について実験・検証例を行った場合の画像を示す図である。図３４（Ａ）は最初の画像、図３４（Ｂ）は最後の画像を示している。図３５は図３４に示す画像から６つの位相特異点に注目し、その軌跡を追跡したものである。【０１８７】
このように、本実施例のポイントマッチング処理を用いることにより、位相特異点を連続して追跡できる。【０１８８】
さらに、任意の位相特異点の組み合わせに対して位相特異点が持つ位相構造と他の位相特異点との距離とによって対応する位相特異点の組を検出し、それらの位相特異点の組を追跡することによってランダムな動きを追跡できることがわかる。特徴パラメータが変化してしまう理由で計測不可能であった大きな変位を伴った物体の運動変位計測が可能となった。【０１８９】
また、本実施例では、位相特異点の組のうち式（１）が最小となる組み合わせを検出することにより、マッチングを行ったが、位相特異点の点間総和距離の差を算出し、変位前のパターンと変位後のパターンとで点間総和距離の差が最小となる位相特異点を検出することによりマッチングを行うようにしてもよい。【０１９０】
図３６、図３７は点間総和距離の差によりマッチングを行う動作を説明するための図を示す。【０１９１】
図３６（Ａ）、図３７（Ａ）は変位前の位相特異点の分布、図３６（Ｂ）、図３７（Ａ）は変位後の位相特異点の分布を示す。【０１９２】
まず、変位前の画像から位相特異点ｐ1〜ｐ4を抽出するとともに、変位後の画像から位相特異点ｐ1'、ｐ2'、ｐ4'、ｐ5'を抽出し、変位前の位相特異点ｐ1〜ｐ4ｔと変位後の位相特異点ｐ1'、ｐ2'、ｐ4'、ｐ5'との点間総和距離の差Ｓ1を求める。点間総和距離の差Ｓ1は、
Ｓ1＝（｜ｐ1ｐ3｜−｜ｐ1'ｐ5'｜）2＋（｜ｐ1ｐ2｜−｜ｐ1'ｐ2'｜）2
＋（｜ｐ2ｐ3｜−｜ｐ2'ｐ5'｜）2＋（｜ｐ4ｐ3｜−｜ｐ4'ｐ5'｜）2
＋（｜ｐ1ｐ4｜−｜ｐ1'ｐ4'｜）2＋（｜ｐ2ｐ4｜−｜ｐ2'ｐ4'｜）2
で求められる。【０１９３】
また、変位後の画像から位相特異点ｐ1'、ｐ2'、ｐ4'、ｐ5'とは異なる位相特異点ｐ1"、ｐ2"、ｐ4"、ｐ5"を抽出し、変位前の位相特異点ｐ1〜ｐ4ｔと変位後の位相特異点ｐ1"、ｐ2"、ｐ4"、ｐ5"との点間総和距離の差Ｓ2を求める。点間総和距離の差Ｓ2は、
Ｓ2＝（｜ｐ1ｐ4｜−｜ｐ1"ｐ4"｜）2＋（｜ｐ1ｐ2｜−｜ｐ1"ｐ2"｜）2
＋（｜ｐ2ｐ3｜−｜ｐ2"ｐ3"｜）2＋（｜ｐ3ｐ4｜−｜ｐ3"ｐ4"｜）2
＋（｜ｐ1ｐ3｜−｜ｐ1"ｐ3"｜）2＋（｜ｐ2ｐ4｜−｜ｐ2"ｐ4"｜）2
で求められる。【０１９４】
このとき、例えば、Ｓ2<Ｓ1であれば、変位前の画像から位相特異点ｐ1〜ｐ4は、位相特異点ｐ1"、ｐ2"、ｐ4"、ｐ5"に対応すると判定する。【０１９５】
なお、ここでは、説明を簡単にするために４点について説明したが、２点以上であれば、点間総和距離の差によりマッチングを行うことが可能であることは言うまでもない。【０１９６】
本実施例のポイントマッチング処理は、複数の位相特異点の位置関係を用いるので縮尺、移動そして、回転により変位前後で完全にその位置関係が一致しなくてもポイントマッチングを行うことで一番近い組を選び出すことができるので、本発明はあらゆる変位を計測することができる。【０１９７】
なお、本実施例では、特徴点として位置情報、離心率、渦度、実軸と虚軸との交叉角度であるゼロクロス角を特徴パラメータとして持つ位相特異点を例に説明したが、特徴点はこれに限定されるものではなく、例えば、位置情報だけであってもよい。【０１９８】
以上、本実施例によれば、干渉計なしで空間的なランダムパターンの強度情報を位相情報に対応づけ、位相情報の中に存在する全位相特異点の特徴を基に変位前後で各位相特異点がどこに移動したか特定することにより、計測物体の全体の変位計測に対しても、物体の局所領域における変位計測のいずれに対しても適用することができ、さらに場の回転の計測をも実現した変位検出装置を提供することできる。【０１９９】
本国際出願は、２００６年２月１日に出願した日本国特許出願２００６−０２４８４６号、及び、２００６年１１月８日に出願した日本国特許出願２００６−３０３２３８号に基づく優先権を主張するものであり、日本国特許出願２００６−０２４８４６号、及び、日本国特許出願２００６−３０３２３８号の全内容を本国際出願に援用する。
【図１】
2007088759.tif 000032
【図２】
2007088759.tif 000033
【図３】
2007088759.tif 000034
【図４】
2007088759.tif 000035
【図５】
2007088759.tif 000036
【図６】
2007088759.tif 000037
【図７】
2007088759.tif 000038
【図８】
2007088759.tif 000039
【図９】
2007088759.tif 000040
【図１０】
2007088759.tif 000041
【図１１】
2007088759.tif 000042
【図１２】
2007088759.tif 000043
【図１３】
2007088759.tif 000044
【図１４】
2007088759.tif 000045
【図１５】
2007088759.tif 000046
【図１６】
2007088759.tif 000047
【図１７】
2007088759.tif 000048
【図１８】
2007088759.tif 000049
【図１９】
2007088759.tif 000050
【図２０】
2007088759.tif 000051
【図２１】
2007088759.tif 000052
【図２２】
2007088759.tif 000053
【図２３】
2007088759.tif 000054
【図２４】
2007088759.tif 000055
【図２５】
2007088759.tif 000056
【図２６】
2007088759.tif 000057
【図２７】
2007088759.tif 000058
【図２８】
2007088759.tif 000059
【図２９】
2007088759.tif 000060
【図３０】
2007088759.tif 000061
【図３１】
2007088759.tif 000062
【図３２】
2007088759.tif 000063
【図３３】
2007088759.tif 000064
【図３４】
2007088759.tif 000065
【図３５】
2007088759.tif 000066
【図３６】
2007088759.tif 000067
【図３７】
2007088759.tif 000068
【図３８】
2007088759.tif 000069
【図３９】
2007088759.tif 000070
【手続補正書】
【提出日】平成19年11月30日(2007.11.30)
【手続補正１】
【補正対象書類名】特許請求の範囲
【補正対象項目名】全文
【補正方法】変更
【補正の内容】
【特許請求の範囲】
【請求項１】
取得した画像の強度信号をフーリエ変換するフーリエ変換手順と、
前記フーリエ変換手順でフーリエ変換された信号のエルミート対称性を崩すようなフィルタリング処理を行う手順と、
前記フィルタリング処理の後にフーリエ逆変換することにより複素解析信号を得る手順と、
前記複素解析信号の実部と虚部の値が共にゼロとなる位相特異点を取得する位相特異点取得手順と、
前記位相特異点の位相構造から所定の特徴パラメータを取得し、前記特徴パラメータに基づいて前記位相特異点を特定する位相特異点特定手順と、
前記位相特異点特定手順で取得された特徴パラメータに基づいて変位後の位相特異点の位置を検出することにより変位を検出する変位検出手順とを有することを特徴とする変位検出方法。
【請求項２】
前記位相特異点特定手順で取得される特徴パラメータは、前記位相特異点近傍の前記複素解析信号の実部と虚部の値を平面で補間し、その結果、実部がゼロとなる線分と虚部がゼロとなる線分とが交差する角度であることを特徴とする請求項１記載の変位検出方法。
【請求項３】
前記位相特異点特定手順で取得される特徴パラメータは、前記位相特異点近傍の前記複素解析信号の強度分布の形状の離心率であることを特徴とする請求項１記載の変位検出方法。
【請求項４】
前記位相特異点特定手順で取得される特徴パラメータは、前記位相特異点近傍の前記複素解析信号の実部の勾配と前記複素解析信号の虚部の勾配とのベクトル積であることを特徴とする請求項１記載の変位検出方法。
【請求項５】
前記位相特異点特定手順は、前記位相特異点近傍の前記複素解析信号の実部と虚部の値を平面で補間し、その結果、実部がゼロとなる線分と虚部がゼロとなる線分とが交差する角度と、前記位相特異点近傍の前記複素解析信号の強度分布の形状の離心率と、前記位相特異点近傍の前記複素解析信号の実部の勾配と前記複素解析信号の虚部の勾配とのベクトル積である渦度とを含むことを特徴とする請求項１記載の変位検出方法。
【請求項６】
前記変位検出手順は、前記位相特異点特定手順で特定された位相特異点の変位前後の位置を前記特徴パラメータに基づいて特定する変位前後位置検出手順と、
前記変位前後位置検出手順で検出された変位前後の位相特異点の位置に基づいて前記位相特異点の変位を測定する変位測定手順とを有することを特徴とする請求項１記載の変位検出方法。
【請求項７】
前記変位測定手順は、前記変位前後位置検出手順で検出された変位前後の位相特異点の位置の座標の差分の二乗平方根を位相特異点の変位量として算出することを特徴とする請求項６記載の変位検出方法。
【請求項８】
前記変位測定手順は、前記変位前後位置検出手順で検出された複数組の変位前後の位相特異点の位置を結ぶ線分の複数の垂直二等分線の距離の和が最小となる点を回転中心として求める回転中心取得手順と、
前記変位前後位置検出手順で検出された変位前後の位相特異点の位置と前記回転中心取得手順で取得した前記回転中心とからなる三角形の前記回転中心を挟む辺の角度を回転角として求める回転角度取得手順とを有することを特徴とする請求項６記載の変位検出方法。
【請求項９】
（削除）
【請求項１０】
前記フィルタリング処理は、前記フーリエ変換手順でフーリエ変換された信号に対して中心に位相特異点を持つらせん状の位相特性を持ち、かつ、ドーナツ状の振幅特性を持つラゲールガウスビームフィルタ特性による処理であることを特徴とする請求項１記載の変位検出方法。
【請求項１１】
前記画像の強度信号の平均値を前記画像の各強度信号から減算した画像を前記画像とする前処理手順を有することを特徴とする請求項１記載の変位検出方法。
【請求項１２】
前記フィルタリング処理は、前記フーリエ変換手順でフーリエ変換された信号に対して中心に位相特異点を持つらせん状の位相特性を持つフィルタ特性、又は、ヒルベルト変換による処理であることを特徴とする請求項１記載の変位検出方法。
【請求項１３】
変位前の位相特異点の特徴パラメータと変位後の位相特異点の特徴パラメータとの差異に基づいて変位前の位相特異点に対応する変位後の位相特異点を探しだし、その位相特異点の地理的な位置関係の類似性から変位前の位相特異点に対応する変位後の位相特異点のマッチングを行うマッチング手順を、更に、有することを特徴とする請求項１記載の変位検出方法。
【請求項１４】
前記マッチング手順は、変位前の複数の位相特異点間の特徴パラメータの総和距離と変位後の複数の位相特異点間の特徴パラメータの総和距離との差分とが最小となる特徴点の組を変位前と変位後とで互いにマッチングする特徴点として特定することを特徴とする請求項１３記載の変位検出方法。
【請求項１５】
前記マッチング手順は、変位前の所定の位相特異点の特徴パラメータをＴｉ１、Ｔｉ２、変位後の所定の位相特異点の特徴パラメータをＰｗ（ｉ１）、Ｐｗ（ｉ２）、αを縮尺としたとき、
【数１】
2007088759.tif 000071
が最小となる位相特異点の組を変位前後でマッチングする位相特異点の組として抽出することを特徴とする請求項１３記載の変位検出方法。
【請求項１６】
取得した画像の強度信号をフーリエ変換するフーリエ変換手段と、
前記フーリエ変換手段でフーリエ変換された信号の複素共役の対称性を崩すようなフィルタリング処理を行う手段と、
前記フィルタリング処理の後にフーリエ逆変換することにより複素解析信号を得る手段と、
前記複素解析信号の実部と虚部の値が共にゼロとなる位相特異点を取得する位相特異点取得手段と、
前記位相特異点の位相構造から所定の特徴パラメータを取得し、前記特徴パラメータに基づいて前記位相特異点を特定する位相特異点特定手段と、
前記位相特異点特定手段で取得された特徴パラメータに基づいて変位後の位相特異点の位置を検出することにより変位を検出する変位検出手段とを有することを特徴とする変位検出装置。
【請求項１７】
前記位相特異点特定手段で取得される特徴パラメータは、前記位相特異点近傍の前記複素解析信号の実部と虚部の値を平面で補間し、その結果、実部がゼロとなる線分と虚部がゼロとなる線分とが交差する角度であることを特徴とする請求項１６記載の変位検出装置。
【請求項１８】
前記位相特異点特定手段で取得される特徴パラメータは、前記位相特異点近傍の前記複素解析信号の強度分布の形状の離心率であることを特徴とする請求項１６記載の変位検出装置。
【請求項１９】
前記位相特異点特定手段で取得される特徴パラメータは、前記位相特異点近傍の前記複素解析信号の実部の勾配と前記複素解析信号の虚部の勾配とのベクトル積であることを特徴とする請求項１６記載の変位検出装置。
【請求項２０】
前記位相特異点特定手段は、前記位相特異点近傍の前記複素解析信号の実部と虚部の値を平面で補間し、その結果、実部がゼロとなる線分と虚部がゼロとなる線分とが交差する角度と、前記位相特異点近傍の前記複素解析信号の強度分布の形状の離心率と、前記位相特異点近傍の前記複素解析信号の実部の勾配と前記複素解析信号の虚部の勾配とのベクトル積である渦度とを含むことを特徴とする請求項１６記載の変位検出装置。
【請求項２１】
前記変位検出手段は、前記位相特異点特定手段で特定された位相特異点の変位前後の位置を前記特徴パラメータに基づいて特定する変位前後位置検出手段と、
前記変位前後位置検出手段で検出された変位前後の位相特異点の位置に基づいて前記位相特異点の変位を測定する変位測定手段とを有することを特徴とする請求項１６記載の変位検出装置。
【請求項２２】
前記変位測定手段は、
前記変位前後位置検出手段で検出された変位前後の位相特異点の位置の座標の差の二乗平方根を位相特異点の変位量として算出することを特徴とする請求項２１記載の変位検出装置。
【請求項２３】
前記変位測定手段は、前記変位前後位置検出手段で検出された複数組の変位前後の位相特異点の位置を結ぶ線分の複数の垂直二等分線の距離の和が最小となる点を回転中心として求める回転中心取得手段と、
前記変位前後位置検出手段で検出された変位前後の位相特異点の位置と前記回転中心取得手段で取得した前記回転中心とからなる三角形の前記回転中心を挟む辺の角度を回転角として求める回転角度取得手段とを有することを特徴とする請求項２１記載の変位検出装置。
【請求項２４】
（削除）
【請求項２５】
前記フィルタリング処理は、前記フーリエ変換手段でフーリエ変換された信号に対して中心に位相特異点を持つらせん状の位相特性を持ち、かつ、ドーナツ状の振幅特性を持つラゲールガウスビームフィルタ特性による処理であることを特徴とする請求項１６記載の変位検出装置。
【請求項２６】
前記画像の強度信号の平均値を前記画像の各強度信号から減算した画像を前記画像とする前処理手段を有することを特徴とする請求項１６記載の変位検出装置。
【請求項２７】
前記フィルタリング処理は、前記フーリエ変換手段でフーリエ変換された信号に対して中心に位相特異点を持つらせん状の位相特性を持つフィルタ特性、又は、ヒルベルト変換による処理であることを特徴とする請求項１６記載の変位検出装置。
【請求項２８】
変位前の位相特異点の特徴パラメータと変位後の位相特異点の特徴パラメータとの差異に基づいて変位前の位相特異点に対応する変位後の位相特異点を探しだし、その位相特異点の地理的な位置関係の類似性から変位前の位相特異点に対応する変位後の位相特異点のマッチングを行うマッチング手段を、更に、有することを特徴とする請求項１６記載の変位検出装置。
【請求項２９】
前記マッチング手段は、変位前の複数の位相特異点間の特徴パラメータの総和距離と変位後の複数の位相特異点間の特徴パラメータの総和距離との差分とが最小となる特徴点の組を変位前と変位後とで互いにマッチングする特徴点として特定することを特徴とする請求項２８記載の変位検出装置。
【請求項３０】
前記マッチング手段は、変位前の所定の位相特異点の特徴パラメータをＴｉ１、Ｔｉ２、変位後の所定の位相特異点の特徴パラメータをＰｗ（ｉ１）、Ｐｗ（ｉ２）、αを縮尺としたとき、
【数２】
2007088759.tif 000072
が最小となる位相特異点の組を変位前後でマッチングする位相特異点の組として抽出することを特徴とする請求項２８記載の変位検出装置。
【請求項３１】
コンピュータに、
取得した画像の強度信号をフーリエ変換するフーリエ変換手順と、
前記フーリエ変換手順でフーリエ変換された信号のエルミート対称性を崩すようなフィルタリング処理を行う手順と、
前記フィルタリング処理の後にフーリエ逆変換することにより複素解析信号を得る手順と、
前記複素解析信号の実部と虚部の値が共にゼロとなる位相特異点を取得する位相特異点取得手順と、
前記位相特異点の位相構造から所定の特徴パラメータを取得し、前記特徴パラメータに基づいて前記位相特異点を特定する位相特異点特定手順と、
前記位相特異点特定手順で取得された特徴パラメータに基づいて変位後の位相特異点の位置を検出することにより変位を検出する変位検出手順とを実行させる変位検出プログラム。
【請求項３２】
前記位相特異点特定手順で取得される特徴パラメータは、前記位相特異点近傍の前記複素解析信号の実部と虚部の値を平面で補間し、その結果、実部がゼロとなる線分と虚部がゼロとなる線分とが交差する角度であることを特徴とする請求項３１記載の変位検出プログラム。
【請求項３３】
前記位相特異点特定手順で取得される特徴パラメータは、前記位相特異点近傍の前記複素解析信号の強度分布の形状の離心率であることを特徴とする請求項３１記載の変位検出プログラム。
【請求項３４】
前記位相特異点特定手順で取得される特徴パラメータは、前記位相特異点近傍の前記複素解析信号の実部の勾配と前記複素解析信号の虚部の勾配とのベクトル積であることを特徴とする請求項３１記載の変位検出プログラム。
【請求項３５】
前記位相特異点特定手順は、前記位相特異点近傍の前記複素解析信号の実部と虚部の値を平面で補間し、その結果、実部がゼロとなる線分と虚部がゼロとなる線分とが交差する角度と、前記位相特異点近傍の前記複素解析信号の強度分布の形状の離心率と、前記位相特異点近傍の前記複素解析信号の実部の勾配と前記複素解析信号の虚部の勾配とのベクトル積である渦度とを含むことを特徴とする請求項３１記載の変位検出プログラム。
【請求項３６】
前記変位検出手順は、前記位相特異点特定手順で特定された位相特異点の変位前後の位置を前記特徴パラメータに基づいて特定する変位前後位置検出手順と、
前記変位前後位置検出手順で検出された変位前後の位相特異点の位置に基づいて前記位相特異点の変位を測定する変位測定手順とを有することを特徴とする請求項３１記載の変位検出プログラム。
【請求項３７】
前記変位測定手順は、前記変位前後位置検出手順で検出された変位前後の位相特異点の位置座標の差分の二乗平方根を位相特異点の変位量として算出することを特徴とする請求項３６記載の変位検出プログラム。
【請求項３８】
前記変位測定手順は、前記変位前後位置検出手順で検出された複数組の変位前後の位相特異点の位置を結ぶ線分の複数の垂直二等分線の距離の和が最小となる点を回転中心として求める回転中心取得手順と、
前記変位前後位置検出手順で検出された変位前後の位相特異点の位置と前記回転中心取得手順で取得した前記回転中心とからなる三角形の前記回転中心を挟む辺の角度を回転角として求める回転角度取得手順とを有することを特徴とする請求項３６記載の変位検出プログラム。
【請求項３９】
（削除）
【請求項４０】
前記フィルタリング処理は、前記フーリエ変換手順でフーリエ変換された信号に対して中心に位相特異点を持つらせん状の位相特性を持ち、かつ、ドーナツ状の振幅特性を持つラゲールガウスビームフィルタ特性による処理であることを特徴とする請求項３１記載の変位検出プログラム。
【請求項４１】
前記画像の強度信号の平均値を前記画像の各強度信号から減算した画像を前記画像とする前処理手順を有することを特徴とする請求項３１記載の変位検出プログラム。
【請求項４２】
前記フィルタリング処理は、前記フーリエ変換手順でフーリエ変換された信号に対して中心に位相特異点を持つらせん状の位相特性を持つフィルタ特性、又は、ヒルベルト変換による処理であることを特徴とする請求項３１記載の変位検出プログラム。
【請求項４３】
変位前の位相特異点の特徴パラメータと変位後の位相特異点の特徴パラメータとの差に基づいて前記変位前に特定された位相特異点と変位後に特定された位相特異点とをマッチングを行うマッチング手順を、更に、有することを特徴とする請求項３１記載の変位検出プログラム。
【請求項４４】
前記マッチング手順は、変位前の複数の位相特異点間の特徴パラメータの総和距離と変位後の複数の位相特異点間の特徴パラメータの総和距離との差分とが最小となる特徴点の組を変位前と変位後とで互いにマッチングする特徴点として特定することを特徴とする請求項４３記載の変位検出プログラム。
【請求項４５】
前記マッチング手順は、変位前の所定の位相特異点の特徴パラメータをＴｉ１、Ｔｉ２、変位後の所定の位相特異点の特徴パラメータをＰｗ（ｉ１）、Ｐｗ（ｉ２）、αを縮尺としたとき、
【数３】
2007088759.tif 000073
が最小となる位相特異点の組を変位前後でマッチングする位相特異点の組として抽出することを特徴とする請求項４３記載の変位検出プログラム。
【請求項４６】
取得した画像の強度信号をフーリエ変換するフーリエ変換手順と、
前記フーリエ変換手順でフーリエ変換された信号のエルミート対称性を崩すようなフィルタリング処理を行う手順と、
前記フィルタリング処理の後にフーリエ逆変換することにより複素解析信号を得る手順と、
前記複素解析信号の実部と虚部の値が共にゼロとなる位相特異点を取得する位相特異点取得手順と、
前記位相特異点の位相構造から所定の特徴パラメータを取得し、前記特徴パラメータに基づいて前記位相特異点を特定する位相特異点特定手順と、
前記位相特異点特定手順で取得された特徴パラメータに基づいて変位後の位相特異点の位置を検出する手順とを有し、
変位前の位相特異点の特徴パラメータと変位後の位相特異点の特徴パラメータとの差異に基づいて変位前の位相特異点に対応する変位後の位相特異点を探しだし、その位相特異点の地理的な位置関係の類似性から変位前の位相特異点に対応する変位後の位相特異点のマッチングを行うことを特徴とする位相特異点マッチング処理方法。
【請求項４７】
変位前の複数の位相特異点間の特徴パラメータの総和距離と変位後の複数の位相特異点間の特徴パラメータの総和距離との差分とが最小となる位相特異点の組を変位前と変位後とで互いにマッチングする位相特異点として特定することを特徴とする請求項４６記載の位相特異点マッチング処理方法。
【請求項４８】
変位前の所定の位相特異点の特徴パラメータをＴｉ１、Ｔｉ２、変位後の所定の位相特異点の特徴パラメータをＰｗ（ｉ１）、Ｐｗ（ｉ２）、αを縮尺としたとき、
【数４】
2007088759.tif 000074
が最小となる位相特異点の組を変位前後でマッチングする位相特異点の組として抽出することを特徴とする請求項４６記載の位相特異点マッチング処理方法。
【請求項４９】
（削除）
【請求項５０】
前記位相特異点は、前記特徴パラメータとして、位置情報、前記位相特異点近傍の前記複素解析信号の強度分布の形状の離心率、前記位相特異点近傍の前記複素解析信号の実部の勾配と前記複素解析信号の虚部の勾配とのベクトル積である渦度、及び前記位相特異点近傍の前記複素解析信号の実部がゼロとなる線分と虚部がゼロとなる線分とが交わる角度を有することを特徴とする請求項４６記載の位相特異点マッチング処理方法。
【請求項５１】
コンピュータに、
取得した画像の強度信号をフーリエ変換するフーリエ変換手順と、
前記フーリエ変換手順でフーリエ変換された信号のエルミート対称性を崩すようなフィルタリング処理を行う手順と、
前記フィルタリング処理の後にフーリエ逆変換することにより複素解析信号を得る手順と、
前記複素解析信号の実部と虚部の値が共にゼロとなる位相特異点を取得する位相特異点取得手順と、
前記位相特異点の位相構造から所定の特徴パラメータを取得し、前記特徴パラメータに基づいて前記位相特異点を特定する位相特異点特定手順と、
前記位相特異点特定手順で取得された特徴パラメータに基づいて変位後の位相特異点の位置を検出する手順とを有し、
変位前の位相特異点の特徴パラメータと変位後の位相特異点の特徴パラメータとの差異に基づいて変位前の位相特異点に対応する変位後の位相特異点を探しだし、その位相特異点の地理的な位置関係の類似性から変位前の位相特異点に対応する変位後の位相特異点のマッチングを行うことを特徴とするコンピュータにより読み取り可能な位相特異点マッチング処理プログラム。
【請求項５２】
変位前の複数の位相特異点間の特徴パラメータの総和距離と変位後の複数の位相特異点間の特徴パラメータの総和距離との差分とが最小となる位相特異点の組を変位前と変位後とで互いにマッチングする位相特異点として特定することを特徴とする請求項５１記載の位相特異点マッチング処理プログラム。
【請求項５３】
変位前の所定の位相特異点の特徴パラメータをＴｉ１、Ｔｉ２、変位後の所定の位相特異点の特徴パラメータをＰｗ（ｉ１）、Ｐｗ（ｉ２）、縮尺をαとしたとき、
【数５】
2007088759.tif 000075
が最小となる位相特異点の組を変位前後でマッチングする位相特異点の組として抽出することを特徴とする請求項５１記載の位相特異点マッチング処理プログラム。
【請求項５４】
（削除）
【請求項５５】
前記位相特異点は、前記特徴パラメータとして、位置情報、前記位相特異点近傍の前記複素解析信号の強度分布の形状の離心率、前記位相特異点近傍の前記複素解析信号の実部の勾配と前記複素解析信号の虚部の勾配とのベクトル積である渦度、及び前記位相特異点近傍の前記複素解析信号の実部がゼロとなる線分と虚部がゼロとなる線分とが交わる角度を有することを特徴とする請求項５１記載の位相特異点マッチング処理プログラム。
【手続補正３】
【補正対象書類名】明細書
【補正対象項目名】０００１
【補正方法】変更
【補正の内容】
【０００１】
技術分野
［０００１］
本発明は変位検出方法、及び、変位検出装置、変位検出プログラム、並びに、位相特異点マッチング処理方法、位相特異点マッチング処理プログラムに係り、特に、変位前後の位相特異点から変位を検出する変位検出方法、及び、変位検出装置、変位検出プログラム、並びに、位相特異点マッチング処理方法、位相特異点マッチング処理プログラムに関する。
背景技術
［０００２］
近年、レーザースペックルパターンやランダムドットパターンなどのように空間的にランダムな構造をもつパターンやテクスチャを指標として物体の微小変位を非接触で計測する技術は非破壊検査や材料強度試験など産業応用上の重要な位置を占めている。パターンやテクスチャを用いて微小変位を検出する際の信号処理技術として相関法が知られている。
［０００３］
図３８、図３９は相関法の動作を説明するため図を示す。図３８（Ａ）は場の変位、図３８（Ｂ）相関法の概念図、図３８（Ｃ）は本測定法の概念図を示している。また、図３９（Ａ）は物体Ａの一方向への変位、図３９（Ｅ）は物体Ａの回転変位を説明するための図を示している。
［０００４］
相関法を用いて変位を検出する場合、図３８（Ａ）に矢印で示すように場の変位が位置により異なるときには図３８（Ｂ）に示すように各領域に分けて計算する必要があり、各点での相関値を決定するごとに２次元相関積分計算または２次元フーリエ変換・逆変換を行う必要があったため、１ピクセル以下での変位検出の分解能を上げようとしてピクセル内部での相関値の計算点数を増やすと計算時間が膨大になっていた。
［０００５］
また、従来の相関法では、図３９（Ａ）に示すように一方向の変位に対しては相関をとることにより変位量を求めることができるが、図３９（Ｂ）に示すように回転に対しては適用できなかった。
［０００６］
なお、位相情報を利用する相関関数の計算方法としては、位相限定相関法が知ら
【手続補正４】
【補正対象書類名】明細書
【補正対象項目名】０００２
【補正方法】変更
【補正の内容】
【０００２】
れている（特許文献１参照）。この相関法は、複素関数である合成スペクトルの振幅を一定化又は対数関数等により抑制し、位相情報のみからなる振幅限定複素合成スペクトルを作成してこれをフーリエ逆変換して相互相関関数を求めるものである。
特許文献１：特許第３０３５６５４号
発明の開示
発明が解決しようとする課題
［０００７］
しかるに、従来の変位検出方法では、計算時間が膨大になるとともに、回転変位の検出が容易に行えないなどの問題点があった。
［０００８］
本発明は上記の点に鑑みてなされたもので、位相特異点の構造を特徴付ける特徴パラメータから所望の位相特異点を特定することにより変位前後の位相特異点を確実に特定し、その変位を容易、かつ、確実に、さらには、回転変位を含む種々の変位を検出できる変位検出方法、及び、変位検出装置、変位検出プログラム、並びに、位相特異点マッチング処理方法、位相特異点マッチング処理プログラムを提供することを目的とする。
課題を解決するための手段
［０００９］
本発明は、取得した画像の強度信号をフーリエ変換するフーリエ変換手順と、フーリエ変換手順でフーリエ変換された信号のエルミート対称性を崩すようなフィルタリング処理を行う手順と、フィルタリング処理の後にフーリエ逆変換することにより複素解析信号を得る手順と、複素解析信号の実部と虚部の値が共にゼロとなる位相特異点を取得する位相特異点取得手順と、位相特異点の位相構造から所定の特徴パラメータを取得し、特徴パラメータに基づいて前記位相特異点を特定する位相特異点特定手順と、位相特異点特定手順で取得された特徴パラメータに基づいて変位後の位相特異点の位置を検出することにより変位を検出する変位検出手順とを有することを特徴とする。
［００１０］
このとき、取得される特徴パラメータは、位相特異点近傍の実部と虚部の値を平面で補間し、その結果、実部がゼロとなる線分と虚部がゼロとなる線分とが交差する角度、位相特異点近傍の強度分布の形状の離心率、位相特異点近傍の強度の勾配であることを特徴とする。
［００１１］
また、特定された位相特異点の変位前後の位置を特徴パラメータに基づいて特定し、特定された変位前後の位相特異点の位置に基づいて位相特異点の変位を測定することを特徴とする。このとき、変位前後の位相特異点の位置の座標の差分の二乗平方根を位相特異点の変位量として算出することを特徴とする。
［００１２］
複数組の変位前後の位相特異点の位置を結ぶ線分の複数の垂直二等分線の距離の和が最小となる点を回転中心として求め、変位前後の位相特異点の位置と回転中心とからなる三角形の回転中心を挟む辺の角度を回転角として求めることを特徴と
【手続補正５】
【補正対象書類名】明細書
【補正対象項目名】０００３
【補正方法】変更
【補正の内容】
【０００３】
する。
［００１３］
さらに、本発明は、取得した画像の強度信号をフーリエ変換し、フーリエ変換された信号の複素共役の対称性を崩すようなフィルタリング処理を行い、擬似位相情報を取得し、取得した擬似位相情報に基づいて位相特異点を取得することを特徴とする。
［００１４］
このとき、フーリエ変換された信号に対して中心に位相特異点を持つらせん状の位相特性を持ち、かつ、ドーナツ状の振幅特性を持つラゲールガウスビームフィルタ特性によりフィルタリングすることを特徴とする。また、画像の強度信号の平均値を画像の各強度信号から減算した画像を元画像とすることを特徴とする。
［００１５］
さらに、フーリエ変換された信号に対して中心に位相特異点を持つらせん状の位相特性を持つフィルタ特性、又は、ヒルベルト変換によりフィルタリングを行い、擬似位相情報を取得することを特徴とする。
［００１６］
また、本発明は、変位前に特定された位相特異点と変位後に特定された位相特異点とをマッチングを行うポイントマッチング処理方法であって、変位前の位相特異点の特徴パラメータと変位後の位相特異点の特徴パラメータとの差異に基づいて変位前の位相特異点に対応する変位後の位相特異点を探しだし、その位相特異点の地理的な位置関係の類似性から変位前の位相特異点に対応する変位後の位相特異点のマッチングを行うことを特徴とする。
［００１７］
また、このとき、本発明は、変位前の複数の位相特異点間の特徴パラメータの総和距離と変位後の複数の位相特異点間の特徴パラメータの総和距離との差分とが最小となる位相特異点の組を変位前と変位後とで互いにマッチングする特徴点として特定することを特徴とする。
［００１８］
さらに、本発明は、変位前の所定の位相特異点の特徴パラメータをＴｉ１、Ｔｉ２、変位後の所定の位相特異点の特徴パラメータをＰｗ（ｉ１）、Ｐｗ（ｉ２）、αを縮尺としたとき、
［００１９］
［数１］
2007088759.tif 000076
【手続補正６】
【補正対象書類名】明細書
【補正対象項目名】０００４
【補正方法】変更
【補正の内容】
【０００４】
が最小となる位相特異点の組を変位前後でマッチングする位相特異点の組として抽出することを特徴とする。
［００２０］
なお、位相特異点は、画像の強度情報に基づいて取得される位相特異点であることを特徴とし、位相特異点は特徴パラメータとして位置情報、離心率、及び、渦度、並びに、実軸と虚軸とが交わる角度を有することを特徴とする。
発明の効果
［００２１］
本発明によれば、位相特異点の位相構造から所定の特徴パラメータを取得し、特徴パラメータに基づいて位相特異点を特定することにより、変位前後の各位相特異点を一意に特定することができ、これによって、変位後の位相特異点を追跡できるなどの効果を得られる。
図面の簡単な説明
［００２２］
［図１］本発明の一実施例のシステム構成図である。
［図２］変位検出プログラムの処理フローチャートである。
［図３］位相特異点取得処理の処理フローチャートである。
［図４］位相特異点取得処理の動作説明図である。
［図５］位相特異点取得処理の動作説明図である。
［図６］位相特異点の位置をサブピクセルで決定する動作を説明するための図である。
［図７］位相特異点の位置をサブピクセルで決定する動作を説明するための図である。
［図８］位相特異点の位置をサブピクセルで決定する動作を説明するための図である。
［図９］位相特異点特定処理の処理フローチャートである。
［図１０］離心率ｅ、渦度ωを取得する動作を説明するための図である。
［図１１］離心率ｅ、渦度ωを取得する動作を説明するための図である。
［図１２］離心率ｅ、渦度ωを取得する動作を説明するための図である。
［図１３］全位相特異点の位置の差を求めｘ方向，ｙ方向の変位をヒストグラムに表した図である。
［図１４］回転変位の検出動作を説明するための図である。
［図１５］シミュレーション実験を説明するための図である。
［図１６］回転前後の位相特異点の数をヒストグラムで表した図である。
［図１７］回転前後の位相特異点の回転角θｉをヒストグラムで表した図である。
【手続補正７】
【補正対象書類名】明細書
【補正対象項目名】０００７
【補正方法】変更
【補正の内容】
【０００７】
［００３２］
〔変位検出プログラム〕
図２は変位検出プログラムの処理フローチャートを示す。
［００３３］
処理部１３２は、まず、変位検出プログラムに従って、ステップＳ１−１で変位前後の対象画像をハードディスクドライブ１３５からメモリ１３４などに読み出す。
［００３４］
次に処理部１３２は、ステップＳ１−２で、読み出された変位前後の対象画像に対して各々位相特異点を取得するための位相特異点取得処理を実行する。
［００３５］
次に処理部１３２は、ステップＳ１−３で変位前後の対象画像から取得した位相特異点が持つ、実部と虚部とのゼロクロス角度、離心率、渦度などの特徴パラメータにより位相特異点を特定する位相特異点特定処理を実行する。
［００３６］
次に処理部１３２は、ステップＳ１−４で、実部と虚部とのゼロクロス角度、離心率、渦度などの位相特異点に特有の特徴パラメータにより変位前後で対応する位相特異点を特定し、変位前後の位相特異点の位置から対象画像の変位を検出する変位検出処理を実行する。
［００３７］
位相特異点によって、対象画像の変位を検出することで対象画像の変位を正確に特定することが可能となる。
［００３８］
〔複素解析信号取得処理〕
次にステップＳ１−２の位相特異点取得処理について説明する。
［００３９］
図３は位相特異点取得処理の処理フローチャート、図４、図５は位相特異点取得処理の動作説明図を示す。
［００４０］
ステップＳ１−２の位相特異点取得処理は、取得した画像の強度信号をフーリエ変換し、フーリエ変換された信号の複素共役の対称性を崩すようなフィルタリング処理を行い、擬似位相情報を取得し、取得した擬似位相情報に基づいて位相特異点を取得する処理である。
［００４１］
処理部１３２は、まず、ステップＳ２−１で前処理を行う。ステップＳ２−１の前処理は、対象画像の全体の強度情報の平均値を求め、求めた全体画像の強度情報の平均値を対象画像の強度情報から減算する処理を行う処理であり、この処理によって、画像データ強度情報中の直流成分が除去される。直流成分が除去されることにより、位相特異点が強調される。
【手続補正８】
【補正対象書類名】明細書
【補正対象項目名】００１１
【補正方法】変更
【補正の内容】
【００１１】
［００６１］
擬似位相は（−π，π）までの値をとり、図６、図７において黒から白になるにしたがって値が大きくなる。
［００６２］
位相特異点Ｐは、その点を囲む閉経路上の位相勾配の積分が±２πの整数倍になり、また、複素解析信号の実部と虚部が共にゼロとなるという性質を持っている。
［００６３］
図６（Ｂ）、図７（Ｃ）に示すような互いに隣接しあう４画素を結ぶ正方形閉経路に対して位相勾配の積分が±２πの整数倍になるという性質を用いて位相特異点の存在している位置を決定していた。このため、１ピクセル以上の分解能をもって位相特異点の位置を決定することができない。
［００６４］
本実施例では、複素解析信号の実部と虚部が共にゼロとなるという性質を用いて１ピクセルの分解能より小さいサブピクセルの分解能で位相特異点の位置を決定するようにしている。本実施例では、図７（Ｂ）、図７（Ｃ）に示す特異点近傍の離散化された実部と虚部を平面で線形補間することにより図７（Ｅ）、図７（Ｆ）に示すように実部及び虚部を共に滑らかな構造として、実部がゼロとなる線分Ｌｒｅ０、虚部がゼロとなる線分Ｌｉｍ０を求める。図７（Ｄ）に示すように位相特異点の位置は実部がゼロとなる線分Ｌｒｅ０と虚部がゼロとなる線分Ｌｉｍ０の交わる点を位相特異点Ｐ０として設定することによりサブピクセルの分解能をもって位相特異点を決定することができる．
図８（Ａ）に示すように位相勾配の積分が±２πの整数倍になる性質を用いる方法では、１ピクセルの面積を持つ正方形の閉経路内のどこかに位相特異点が存在することはわかってもその位置を特定することができなかった。これに対して図８（Ｂ）に示すように複素解析信号の実部と虚部が共にゼロとなるという性質を用いて実部と虚部の内挿と実部と虚部のゼロの交差点から位相特異点を求めることにより１ピクセルより小さいサブピクセルの分解能をもって位相特異点の位置を決定できている。
［００６５］
〔位相特異点特定処理〕
図９は位相特異点特定処理の処理フローチャートを示す。
［００６６］
ステップＳ１−３の位相特異点特定処理は、位相特異点の位相構造から所定の特徴パラメータを取得し、特徴パラメータに基づいて位相特異点を特定する処理である。
［００６７］
処理部１３２は、ステップＳ３−１で位相特異点を特定する特徴パラメータとしてゼロクロス角を取得する。ゼロクロス角は、図７（Ｄ）、図８（Ｂ）に示す実部がゼロとなる線分Ｌｒｅ０と虚
【手続補正９】
【補正対象書類名】明細書
【補正対象項目名】００１２
【補正方法】変更
【補正の内容】
【００１２】
部がゼロとなる線分Ｌｉｍ０との交差する角度Δθである。
［００６８］
また、処理部１３２は、ステップＳ３−２で位相特異点を特定する特徴パラメータとして離心率ｅを取得し、ステップＳ３−３で渦度ωを取得する。
［００６９］
図１０、図１１、図１２は離心率ｅ、渦度ωを取得する動作を説明するための図を示す。図１０（Ａ）は擬似位相情報、図１０（Ｂ）は擬似位相情報の位相特異点近傍の振幅分布、図１０（Ｃ）は位相分布を示している。
［００７０］
図１０（Ｂ）に示すように振幅分布は位相特異点近傍に近づくにつれて暗くなっていて、これは値がゼロに近づいていることを示している。
［００７１］
また、図１１は特異点近傍の振幅分布の構造図を示している。図１１（Ａ）は振幅を等高線で表示したものを示しており、図１１（Ｂ）は三次元表示したものを示す。
［００７２］
図１１（Ａ）に示すように位相特異点近傍の等振幅の分布は略楕円形状をしている。一般に位相特異点近傍の振幅分布は略楕円形をしており、その形状は位相特異点に応じて異なっており、位相特異点を特定する特徴パラメータとして用いることができる。この楕円形状の特徴を二次曲線の特徴を表す値である離心率によって示すことによって、位相特異点を特定することが可能となる。
［００７３］
例えば、図１２に示すような楕円の場合、一般に数式で、
［００７４］
［数５］
2007088759.tif 000077
で表せる。図１２に示す楕円の離心率ｅはａ：長軸半径，ｂ：短軸半径，ｆ：焦点距離とすると、離心率ｅは
［００７５］
［数６］
【手続補正１０】
【補正対象書類名】明細書
【補正対象項目名】００１３
【補正方法】変更
【補正の内容】
【００１３】
2007088759.tif 000078
で表させる。
［００７６］
図１１（Ｂ）に示される振幅を三次元で表示した図を見ると位相特異点近傍の形状は楕円形の円錐のような構造をしていることが分かる。したがって、本実施例では、各位相特異点を特定する特徴パラメータの一つとして位相特異点近傍の渦度を求める。渦度とは微小部分の局所的な回転の強さと回転の方向を表す量である。図１１（Ａ）では振幅の等高線の間隔が狭いほど渦度は大きくなり、図１１（Ｂ）では三次元表示した位相特異点近傍の勾配（変化の度合い）が急になるほど渦度は大きくなる。すなわち、渦度は位相特異点近傍の振幅の大きさに依存する。
［００７７］
渦度ωは次式で表される。
［００７８］
［数７］
2007088759.tif 000079
以上のように位相特異点は位置（ｘ、ｙ），ゼロクロス角Δθ，離心率ｅ，渦度ωの特徴パラメータを持っていて、これらの特徴パラメータを用いて各位相特異点を一意に特定することができる。このとき、複数の特徴パラメータを用いることにより確実に位相特異点を特定することができる。
［００７９］
〔直線変位検出〕
次にステップＳ１−４の位相特異点変位検出処理について説明する。
［００８０］
位相変位検出処理では、変位前後の位相特異点を上記の特定方法により特定す
【手続補正１１】
【補正対象書類名】明細書
【補正対象項目名】００１４
【補正方法】変更
【補正の内容】
【００１４】
ることにより変位量を検出する。
［００８１］
このとき、まず、変位前後で位相特異点をそれぞれ一意に特定するためのパラメータとして前述のように求めたゼロクロス角Δθ，離心率ｅ，渦度ωを用いて位相特異点を特定する。その方法は変位前の各パラメータをΔθ１、ｅ１、ω１、変位後の各パラメータをΔθ２、ｅ２、ω２とする。
［００８２］
上記パラメータ（Δθ１、ｅ１、ω１）、（Δθ２、ｅ２、ω２）のそれぞれの差（Δθ２−Δθ１、ｅ２−ｅ１、ω２−ω１）の絶対値がある値より小さい時に同一の位相特異点であるとみなす。各位相特異点をそれぞれ一意に特定できるため、変位前後で各位相特異点がどこに移動したのか分かる。このとき、変位前の位相特異点の位置（ｘ１、ｙ１）と変位後の位相特異点の位置（ｘ２、ｙ２）の差から変位量を求めることができる。
［００８３］
実際の変位量ΔＬは、Δｘ＝（ｘ２−ｘ１）、Δｙ＝（ｙ２−ｙ１）とすると、
［００８４］
［数８］
2007088759.tif 000080
で表すことができる。
［００８５］
なお、図１３は、全位相特異点の位置の差を求めｘ方向，ｙ方向の変位をヒストグラムに表した図を示す。図１３（Ａ）はｘ方向の変位、図１３（Ｂ）はｙ方向の変位を示している。
［００８６］
図１３に示すように全位相特異点の位置がｘ方向に略平行に変位したことがわかる。
［００８７］
〔回転変位検出〕
また、本実施例では、各位相特異点を回転変位前後で特定できるため、回転の検出も可能である。
［００８８］
図１４は回転変位の検出動作を説明するための図を示す。
［００８９］
図１４（Ａ）に示すように回転前の２組の位相特異点Ｐ１１、Ｐ１２のｘｙ座標をそれぞれ
【手続補正１２】
【補正対象書類名】明細書
【補正対象項目名】００１６
【補正方法】変更
【補正の内容】
【００１６】
１０２４×１０２４で、回転させる前後の強度情報からＬＧフィルタを用いて２つの解析信号を得たときの振幅分布と位相分布を表示したものである。図１５（Ａ）は回転変位前の擬似位相情報、図１５（Ｂ）は回転変位前の擬似位相情報の振幅分布、図１５（Ｃ）は回転変位前の擬似位相情報の位相分布、図１５（Ｄ）は回転変位後の擬似位相情報、図１５（Ｅ）は回転変位後の擬似位相情報の振幅分布、図１５（Ｆ）は回転変位後の擬似位相情報の位相分布を示している。
［００９７］
この回転前後の位相特異点の数は両方とも正の位相特異点が３２２個，負の位相区異点３２１個の合計６４３個で同数の位相特異点が存在していた。図１５（Ｂ）及び図１５（Ｅ）に示される振幅分布からも明らかなように回転前の三角で囲った２つのリング状の振幅が回転後には回転の中心（全画素の中心）に対して９０°回転している。
［００９８］
次に位相特異点を特定するための特徴パラメータであるＲｅ＝０とＩｍ＝０の直線の交わる角度：Δθ，離心率（ｅｃｃｅｎｔｒｉｃｉｔｙ）：ｅ，渦度（ｖｏｒｔｉｃｉｔｙ）：ωを回転前後の全位相特異点に対して求める。回転前後で対応する位相特異点を探す条件は以下の条件を用いた。
［００９９］
条件１：Ｒｅ＝０とＩｍ＝０の直線の交わる角度Δθの差が±０．５度以下
条件２：回転前後の離心率ｅの差が±０．００５以下
条件３：回転前後の渦度ωの差が±０．０５以下
回転前の位相特異点に対応する回転後の位相特異点が１対１の対応になっているのかを調べた。
［０１００］
図１６は回転前の位相特異点に対応する回転後の位相特異点の数をヒストグラムで表した図を示す。
［０１０１］
回転前後でほぼ全ての位相特異点が一意に特定できた。
シミュレーションで得られた測定結果の回転の中心は
（ｘｃ、ｙｃ）＝（５１１．５、５１１．５）
となった。これは設定した回転の中心
（ｘｃ、ｙｃ）＝（５１１．５、５１１．５）
と一致する。
［０１０２］
図１７は回転前の位相特異点に対応する回転後の位相特異点と回転の中心から求めたそれぞれの回転角θｉをヒストグラムで表した図を示す。
【手続補正１３】
【補正対象書類名】明細書
【補正対象項目名】００１７
【補正方法】変更
【補正の内容】
【００１７】
［０１０３］
シミュレーション結果である回転角のヒストグラムは±０．００１度の精度で全位相特異点が９０度回転していることを示している。シミュレーション実験では、測定対象１０２を９０度回転させており、図１７に示す測定結果も９０度回転していることを示しているので、この測定結果から本発明が回転に対しても有効に計測できることを確認できる。
［０１０４］
図１８は回転前後の対応する位相特異点の角度Δθの差をヒストグラムで表した図、図１９は回転前後の対応する位相特異点の離心率ｅの差をヒストグラムで表した図、図２０は回転前後の対応する位相特異点の渦度ωの差をヒストグラムで表した図を示す。
［０１０５］
図１８、図１９、図２０に示すようにゼロ付近に固まっており、誤差が少なく、位相特異点を確実に特定できていることがわかる。
［０１０６］
なお、上記実施例では、位相特異点を用いて回転変位を検出する処理について説明したが縮尺、移動、回転を簡単な処理により検出することも可能である。
［０１０７］
〔ポイントマッチング処理〕
ポイントマッチング処理は、任意の複数点の位置関係と、その点の空間構造の類似性を用いて変位前のパターンから変位後のパターンを検出する処理である。
［０１０８］
すなわち、位相構造が類似する点を候補点として、その候補点を含んだパターン内の点間の距離が最小となる点を変位後のパターンとして決定するというものである。
［０１０９］
以下にポイントマッチング処理の手順を説明する。
［０１１０］
まず、変位前後の２つの強度分布画像に対して各々前述の位相特異点を抽出するためのフィルタ演算を行うことにより位相特異点を抽出する。抽出した位相特異点近傍にある実部及び虚部を線形補間することによって、サブピクセルの精度で全位相特異点の位置を特定する。更に、位置が特定された各位相特異点に対して離心率ｅ、実軸と虚軸との角Δθ、渦度ωを求め、ラベリングを行う。
［０１１１］
図２１は、ポイントマッチング処理の動作説明図を示す。図２１（Ａ）は変位前の位相特異点の状態、図２１（Ｂ）は変位後の位相特異点の状態を示している。
［０１１２］
まず、図２１（Ａ）に示す位相特異点から黒丸で示す位相特異点を選択する。選択する位相特異点はその特徴パラメータである離心率ｅ、実軸と虚軸との角Δθ、渦度ωが他の多くの位相特異点とは明らかにことなるものを抽出することが望ましい。これ
【手続補正１４】
【補正対象書類名】明細書
【補正対象項目名】００１９
【補正方法】変更
【補正の内容】
【００１９】
［０１２３］
［数１０］
2007088759.tif 000081
式（１）においてＴｉとＰｗ（ｉ）は変位前と変位後の点を表している。なお、ＴｉはＰｗ（ｉ）と対応している。
Ｔ＝｛Ｔ１，Ｔ２，Ｔ３，・・・，Ｔｎ｝
は変位前の位相特異点の組み合わせを表している。そして、
Ｐ＝｛Ｐ１，Ｐ２，Ｐ３，・・・，Ｐｎ｝
は変位後の位相特異点の組み合わせを表している。
［０１２４］
また、｜Ｔｉ１−Ｔｉ２｜は、点Ｔｉ１と点Ｔｉ２との間の距離，αは縮尺のスケール（度合い）を表している。さらに、変数ｗはＴ→Ｐの対応関係を表している。
［０１２５］
なお、上記式（１）の関数をコスト関数（ｃｏｓｔ  ｆｕｎｃｔｉｏｎ）と呼ぶ。
［０１２６］
処理部１３２は、ステップＳ４−９で変数ｋに（ｋ＋１）を代入し、ステップＳ４−５に戻る。処理部１３２は、上記ステップＳ４−６〜Ｓ４−９をステップＳ４−５で変数ｋが定数ｎより大きくなるまで繰り返す。
［０１２７］
次に処理部１３２は、ステップＳ４−１０で残った組み合わせＭの位相特異点の特徴パラメータ間の距離が最小となるものをマッチングパターンとして抽出する。
［０１２８］
次に、処理部１３２は、ステップＳ４−１１で抽出した特徴パラメータ間の距離が最小となる位相特異点の変換パラメータ、例えば、平行移動（ｓｈｉｆｔ），回転（ｒｏｔａｔｉｏｎ）、縮尺（ｓｃａｌｅ）を計算する。
［０１２９］
上記ポイントマッチング処理により、目標はパターン内の複数の点の位置関係において変位前後で最も良く一致するものを見つけ出すことができる。
［０１３０］
次に、上記マッチングアルゴリズムの処理動作を、図面を用いてさらに詳細に説明する。
［０１３１］
図２３、図２４、図２５はマッチングアルゴリズムの処理動作を説明するための図を示
【手続補正１５】
【補正対象書類名】明細書
【補正対象項目名】００２０
【補正方法】変更
【補正の内容】
【００２０】
す。
［０１３２］
図２３（Ａ）は変位前の位相特異点の中から複数の位相特異点の組を選んだものであり、これをＴ空間とし、ｔｉ（ｉ＝１、．．．，ｎ）、Ｔｉ（Ｘｉ，Ｙｉ，ｅｉ，ωｉ，Δθｉ）で表す。図２３（Ｂ）は変位後の位相特異点でＰ空間とし、ｐｊ（ｊ＝１、．．．，ｋ）、Ｐｊ（Ｘｊ，Ｙｊ，ｅｊ，ωｊ，Δθｊ）で表される。
［０１３３］
まず、図２３（Ａ）に黒丸で示すように、変位前の位相特異点から任意の個数の位相特異点を選ぶ。
［０１３４］
変位前のパターン内のそれぞれの位相特異点ｔｉの特徴パラメータである離心率ｅｔｉ、ゼロクロス角Δθｔｉ、渦度ωｔｉと変位後の位相特異点ｐｊの特徴パラメータである離心率ｅｐｊ、ゼロクロス角Δθｐｊ、渦度ωｐｊとの差を閾値と比較する。それは、例えば、次式で表せる。
［０１３５］
［数１１］
2007088759.tif 000082
なお、式（２）で２番目の分母の｜ωｔｉ＋ωｐｊ｜は正規化するために付け加えたものである。実際のプログラム内ではこの形で使われている。また、ｅｔｈは離心率の閾値、Δθｔｈはゼロクロス角の閾値、ωｔｈは渦度の閾値を示している。
［０１３６］
次にＴ空間の位相特異点の組とＰ空間の位相特異点の組とのマッチングを行う。
［０１３７］
図２４において、実線で示す丸印は正の位相特異点を示しており、破線で示した丸印は負の位相特異点を点線で示している。
［０１３８］
まず、図２４（Ａ）に示すＴ空間の２個の位相特異点の極性、正か負か符号と、その位置関係とに基づいて図２４（Ｂ）に破線の直線で示すようにＰ空間に存在しうる可能性のある２個の位相特異点の組を探索する。
［０１３９］
このマッチングでは移動（ｓｈｉｆｔ），回転（ｒｏｔａｔｉｏｎ），だけでなく縮尺（ｓｃａｌｅ）も考慮しているので、Ｔ空間での初めの２つの組をＰ空間の中から２つを選ぶ組み合わせは、この例では正の位相特異点と負の位相特異点の組み合わせ全てであり、正が７個、負が
【手続補正１６】
【補正対象書類名】明細書
【補正対象項目名】００２１
【補正方法】変更
【補正の内容】
【００２１】
６個なので４２通り存在する。
［０１４０］
なお、ここで、正の位相特異点と負の位相特異点というのは特徴パラメータの渦度（ｖｏｒｔｉｃｉｔｙ）の値が正のものを正の位相特異点、また、負のものを”負の位相特異点としている。
［０１４１］
この得られた位相特異点の組み合わせに対して，位相特異点の特徴パラメータの条件に合わない組み合わせを切り捨てる。このとき、各特徴パラメータは実験環境により変化しやすいことを考え、条件を緩くとることが好ましい。ここで言う条件とは、式（２）の閾値ｅｔｈ、Δθｔｈ、ωｔｈのことである。
［０１４２］
次に、ステップＳ４−３で空間Ｔ内の２つの位相特異点ｔｉ１、ｔｉ２を選ぶ。この位相特異点ｔｉ１、ｔｉ２は少なくとも１つのマッチング候補点を持っている。位相特異点ｔｉ１はｎ１個のマッチング候補の点を、位相特異点ｔｉ２はｎ２個のマッチング候補点を持つ。したがって、変位後の空間Ｐ内の候補点の組み合わせはｎ１ｎ２通りあることがわかる。それをテンポラルマッチング（ｔｅｍｐｏｒａｌ  ｍａｔｃｈｉｎｇ）と呼ぶことにする。
［０１４３］
次に、図２５のようにＴ空間の位相特異点を１つ増やし、Ｐ空間の中から３点目として可能性のある点を探す。これにより３点目を含む組み合わせが得られる。なお、特徴パラメータを用いて明らかに間違った点は除外される。このような操作をＴ空間で選択される位相特異点がｎ個になるまで繰り返す。
［０１４４］
つまり、新しいマッチングポイントを加えつつ、上記テンポラルマッチングを最終的にすべての位相特異点のマッチング点が決定されるまで繰り返す。
［０１４５］
この操作を行うためには、新しく追加されていく点の地理的な位置関係が変位前と変位後で類似性があるかどうか調べる必要がある。
［０１４６］
式（１）では、テンポラルマッチングの類似性を調べると同時に距離の閾値は間違って検出されたマッチングポイントを除外する操作も行っていることになる。
［０１４７］
最終的にで得られた位相特異点の組み合わせＭの中から式（１）の評価式によりＴ空間の全位相特異点間の距離と対応するＰ空間の位相特異点間の空間的な距離及び特徴点の距離の差の総和が最小になるものを選び出してマッチングを行っている。
［０１４８］
さらに、図面を用いて説明を続ける。
【手続補正１７】
【補正対象書類名】明細書
【補正対象項目名】００２３
【補正方法】変更
【補正の内容】
【００２３】
［０１５９］
一般に、流体力学の観点からどんな複雑な動きもそれを平行移動（ｔｒａｎｓｉｔｉｏｎ），回転（ｒｏｔａｔｉｏｎ），縮尺（ｓｃａｌｅ）で表すことができることが知られており、変位前後座標の回転関係式は次式で表すことができる。
［０１６０］
［数１２］
2007088759.tif 000083
ここで、（ｘｉ、ｙｉ）は変位前の位相特異点の座標、（ｘｉ’、ｙｉ’）は変位後の位相特異点の座標、αは縮尺、θは回転角である。
［０１６１］
［数１３］
2007088759.tif 000084
は変位前の位相特異点の重心、
［０１６２］
［数１４］
2007088759.tif 000085
は変位後の位相特異点の重心を表している。
まず、はじめに、マッチングを表すのに適した関数について説明する。
［０１６３］
この関数はそれぞれ変位前後で対応している位相特異点の差を足し合わせた形で書ける。
【手続補正１８】
【補正対象書類名】明細書
【補正対象項目名】００２８
【補正方法】変更
【補正の内容】
【００２８】
示した位相特異点の変位を示す図であり、図３３（Ｂ）は選択された位相特異点を拡大した図である。
［０１８４］
この検証例においても図３２（Ａ）に示すように変位前の４つの位相特異点が変位後、図３２（Ｂ）に示す位置に移動し、図３３に示すように注目した４つの位相特異点が変位していることが検証できた。
［０１８５］
さらに同じ作業をふぐの連続撮影画像１００枚に対して行った。
［０１８６］
図３４は連続した画像について実験・検証例を行った場合の画像を示す図である。図３４（Ａ）は最初の画像、図３４（Ｂ）は最後の画像を示している。図３５は図３４に示す画像から６つの位相特異点に注目し、その軌跡を追跡したものである。
［０１８７］
このように、本実施例のポイントマッチング処理を用いることにより、位相特異点を連続して追跡できる。
［０１８８］
さらに、任意の位相特異点の組み合わせに対して位相特異点が持つ位相構造と他の位相特異点との距離とによって対応する位相特異点の組を検出し、それらの位相特異点の組を追跡することによってランダムな動きを追跡できることがわかる。特徴パラメータが変化してしまう理由で計測不可能であった大きな変位を伴った物体の運動変位計測が可能となった。
［０１８９］
また、本実施例では、位相特異点の組のうち式（１）が最小となる組み合わせを検出することにより、マッチングを行ったが、位相特異点の点間総和距離の差を算出し、変位前のパターンと変位後のパターンとで点間総和距離の差が最小となる位相特異点を検出することによりマッチングを行うようにしてもよい。
［０１９０］
図３６、図３７は点間総和距離の差によりマッチングを行う動作を説明するための図を示す。
［０１９１］
図３６（Ａ）、図３７（Ａ）は変位前の位相特異点の分布、図３６（Ｂ）、図３７（Ｂ）は変位後の位相特異点の分布を示す。
［０１９２］
まず、変位前の画像から位相特異点ｐ１〜ｐ４を抽出するとともに、変位後の画像から位相特異点ｐ１’、ｐ２’、ｐ３’、ｐ４’を抽出し、変位前の位相特異点ｐ１〜ｐ４と変位後の位相特異点ｐ１’、ｐ２’、ｐ３’、ｐ４’との点間総和距離の差Ｓ１を求める。点間総和距離の差Ｓ１は、
【手続補正１９】
【補正対象書類名】明細書
【補正対象項目名】００２９
【補正方法】変更
【補正の内容】
【００２９】
Ｓ１＝（｜ｐ１ｐ３｜−｜ｐ１’ｐ３’｜）２＋（｜ｐ１ｐ２｜−｜ｐ１’ｐ２’｜）２
＋（｜ｐ２ｐ３｜−｜ｐ２’ｐ３’｜）２＋（｜ｐ３ｐ４｜−｜ｐ３’ｐ４’｜）２
＋（｜ｐ１ｐ４｜−｜ｐ１’ｐ４’｜）２＋（｜ｐ２ｐ４｜−｜ｐ２’ｐ４’｜）２
で求められる。
ここで、｜ｐｉｐｊ｜は、点ｐｉと点ｐｊとの間の距離を表している。
［０１９３］
また、変位後の画像から位相特異点ｐ１’、ｐ２’、ｐ３’、ｐ４’とは異なる位相特異点ｐ１″、ｐ２″、ｐ３″、ｐ４″を抽出し、変位前の位相特異点ｐ１〜ｐ４と変位後の位相特異点ｐ１″、ｐ２″、ｐ３″、ｐ４″との点間総和距離の差Ｓ２を求める。点間総和距離の差Ｓ２は、
Ｓ２＝（｜ｐ１ｐ３｜−｜ｐ１″ｐ３″｜）２＋（｜ｐ１ｐ２｜−｜ｐ１″ｐ２″｜）２
＋（｜ｐ２ｐ３｜−｜ｐ２″ｐ３″｜）２＋（｜ｐ３ｐ４｜−｜ｐ３″ｐ４″｜）２
＋（｜ｐ１ｐ４｜−｜ｐ１″ｐ４″｜）２＋（｜ｐ２ｐ４｜−｜ｐ２″ｐ４″｜）２
で求められる。
［０１９４］
このとき、例えば、Ｓ２＜Ｓ１であれば、変位前の画像から位相特異点ｐ１〜ｐ４は、位相特異点１″、ｐ２″、ｐ３″、ｐ４″に対応すると判定する。
［０１９５］
なお、ここでは、説明を簡単にするために４点について説明したが、２点以上であれば、点間総和距離の差によりマッチングを行うことが可能であることは言うまでもない。
［０１９６］
本実施例のポイントマッチング処理は、複数の位相特異点の位置関係を用いるので縮尺、移動そして、回転により変位前後で完全にその位置関係が一致しなくてもポイントマッチングを行うことで一番近い組を選び出すことができるので、本発明はあらゆる変位を計測することができる。
［０１９７］
なお、本実施例では、特徴点として位置情報、離心率、渦度、実部がゼロとなる線分Ｌｒｅ０と虚部がゼロとなる線分Ｌｉｍ０との交叉角度であるゼロクロス角を特徴パラメータとして持つ位相特異点を例に説明したが、特徴パラメータはこれに限定されるものではなく、例えば、位置情報だけであってもよい。
［０１９８］
以上、本実施例によれば、干渉計なしで空間的なランダムパターンの強度情報を位相情報に対応づけ、位相情報の中に存在する全位相特異点の特徴を基に変位前後で各位相特異点がどこに移動したか特定することにより、計測物体の全体の変位計測に対しても、物体の局所領域における変位計測のいずれに対しても適用することができ、さらに場の回転の計測をも実現した変位検出装置を提供することできる。
［０１９９］
本国際出願は、２００６年２月１日に出願した日本国特許出願２００６−０２４８４６号、
【手続補正２０】
【補正対象書類名】図面
【補正対象項目名】図１４
【補正方法】変更
【補正の内容】
【図１４】
2007088759.tif 000086
【手続補正２１】
【補正対象書類名】図面
【補正対象項目名】図２２
【補正方法】変更
【補正の内容】
【図２２】
2007088759.tif 000087
【国際調査報告】
2007088759.tif 000088
2007088759.tif 000089
2007088759.tif 000090