JP2007333647A - 画像処理方法、画像処理装置および画像処理プログラム - Google Patents

Abstract

【課題】画像相関法による変位、変形、歪みを測定する画像処理方法、画像処理装置および画像処理プログラムの提供。
【解決手段】対象物Ｔの変化前後の画像データ取得過程と微小領域４の大きさと解析点３の移動量解析範囲５回転角解析範囲６および歪み解析範囲７の数値データのパラメータ取得過程と微小領域４をアフィン変換によりマスク８を生成する回転処理過程と各マスク８を走査しながら変化後画像データ２との相関値を計算する第一相関値算出過程と、最大相関値となる位置および角度を検出する位置角度検出過程と、最大相関値を示すマスク８をアフィン変換により伸縮させた二次マスク９を生成する伸縮処理過程と各二次マスク９と変化後画像データ２との第二相関値算出過程と最大相関値となる伸縮率を検出する伸縮率検出過程と最大相関値を与える位置、角度および伸縮率から解析点の移動量、回転角および歪みを決定する状態決定過程からなる画像処理方法。
【選択図】図４

Description

本発明は、対象物の変化前と変化後の二枚の画像の相関値を求め最大となった点を移動先とみなす画像相関法により、対象物の変位、変形および歪みを測定する画像処理方法、画像処理装置および画像処理プログラムに関するものである。

従来、物体の変位・変形を画像処理によって求める方法として、画像相関法が用いられている。画像相関法は、変化前と変化後の二枚の画像を用い、一方の画像データ上で他方の画像データの部分領域（基準画像）を走査させながら相関値を計算し、最大相関値が得られる点に基準画像が存在するとみなすものであり、二次元画像の変位量測定では、基準画像を二次元方向に走査させることで二次元の変位量を導き出すことができる。このような画像による測定は、接触式による測定が不可能な測定対象や、測定範囲や測定点の数が多いときに有効である。
ところで、一般に画像相関法による測定では、測定対象が並進運動する場合を前提に解析が行われていた（例えば、特許文献１）。そのため、測定対象が大きく変形し、回転やせん断変形、伸縮等の変形を伴っている場合には解析を行うことが困難であった。
特開平８−１６６２０８号公報

そこで発明者らは、画像相関法において基準画像を幾何学変換により回転させる手法を提案している（非特許文献１）。この手法は、基準画像をある角度ずつ回転させその都度走査させることを繰り返すもので、従来の並進運動のみを考慮する手法では対象が10度程度回転すると解析不能になっていたところ、この手法では30度回転した場合でもおおむね解析が可能である。
また、画像相関法による処理を行う際に、対象画像にアフィン変換を施す発明が提案されている（特許文献２）。ここでアフィン変換とは、平行移動、回転、左右反転、拡大／縮小、せん断の組み合わせによる座標変換である。ただし特許文献２においては、平行移動、回転、拡大／縮小（ｘ方向およびｙ方向に同じ倍率）のみを考慮している。この発明は、対象物の画像データ上を走査してあらかじめ登録されたモデル画像との相関値を求めるもので、アフィン変換のパラメータを変化させながら走査を所定回数繰り返し、最大相関値となる点を導き出す。さらに、何度もアフィン変換することを避けるため、対象を二値化処理して主軸角を求め、それに基づきアフィン変換のパラメータを定める。
廣瀬、田代、園田、野村、神谷、「2次元幾何学変換を用いた画像処理による変位測定に関する研究」、精密工学会秋季大会学術講演会講演論文集、2004年、p.1083-1084 特開平７−２２００２６号公報

しかしながら、非特許文献１の手法においては、基準画像の変換に際して回転しか考慮されていなかったため、連続体を対象とした場合、変形を伴わない移動や回転および連続体の微小部分は歪まないとみなせる小変形については解析可能であったが、微小部分の歪みまで考慮して解析することは不可能であった。
また、特許文献２の発明は、あらかじめ登録されたモデル画像により対象物の認識および位置の計測を行うものであり、連続体の変形を対象とすることはできなかった。さらに、主軸角からアフィン変換のパラメータを定める方法も、連続体の変形を計測するような場合には一つの主軸角を定めることができず、用いることができなかった。

本発明は、上記事情を鑑みたものであり、主に連続体を対象とし、画像相関法によりその変位、変形および歪みを測定することができる画像処理方法、画像処理装置および画像処理プログラムを提供することを目的とする。

本発明のうち請求項１の発明は、画像データ取得手段が対象物の変化前画像データおよび変化後画像データを受け付ける画像データ取得過程と、パラメータ取得手段が変化前画像データ上の解析点の座標および解析点を中心とする微小領域の大きさならびに解析点の移動量解析範囲、回転角解析範囲および歪み解析範囲の数値データを受け付けるパラメータ取得過程と、回転処理手段が微小領域をアフィン変換により回転角解析範囲において所定角度ずつ回転させたマスクを生成する回転処理過程と、第一相関値算出手段が変化後画像データ上の移動量解析範囲において各マスクを走査しながら変化後画像データとの相関値を計算する第一相関値算出過程と、位置角度検出手段が最大相関値となる位置および角度を検出する位置角度検出過程と、伸縮処理手段が最大相関値を示すマスクをアフィン変換により歪み解析範囲において所定伸縮率ずつ伸縮させた二次マスクを生成する伸縮処理過程と、第二相関値算出手段が各二次マスクと変化後画像データとの相関値を計算する第二相関値算出過程と、伸縮率検出手段が最大相関値となる伸縮率を検出する伸縮率検出過程と、状態決定手段が最大相関値を与える位置、角度および伸縮率から、それぞれ解析点の移動量、回転角および歪みを決定する状態決定過程と、を備えることを特徴とする。

また、請求項２の発明は、前記微小領域の形状が、円形であることを特徴とする。

さらに、請求項３の発明は、連続体上の複数点についての変位、変形および歪みの解析を行う場合において、解析範囲決定手段が解析しようとする現解析点の移動量、回転角および歪みのうち少なくとも一つについて、現解析点の近傍の点ですでに解析した前解析点の移動量、回転角および歪みの値を基準とした相対的な解析範囲を定める解析範囲決定過程を備えることを特徴とする。

また、請求項４の発明は、画像処理装置であって、変化前および変化後の対象物を撮像し、変化前画像データおよび変化後画像データとして出力する撮像手段と、請求項１、２又は３記載の画像処理方法を実行する画像処理手段と、を備えることを特徴とする。

さらに、請求項５の発明は、画像処理プログラムであって、コンピュータに、請求項１、２又は３記載の画像処理方法を実行させることを特徴とする。

本発明のうち請求項１の発明によれば、対象物に接触することなく、変化前画像データおよび変化後画像データのみから対象物の変位、変形および歪みを測定することができる。特に、画像相関の際に微小領域に対してアフィン変換による伸縮を適用することで、従来の方法では不可能であった微小部分の歪みを伴う大きな変形についても解析可能であり、また従来可能であった解析についてもより高い精度で行うことができる。さらに、歪みから各点の応力を直接的に把握することも可能である。なお、本発明の具体的な用途としては、金属の大変形の測定を始めとし、橋やビル等の構造体の変形、地滑り等の土砂災害の予兆システム等の解析が考えられる。

本発明のうち請求項２の発明によれば、微小領域の形状を円形とすることで、微小領域をアフィン変換により回転させた場合でも領域内に含まれる計算点数が変化せず、回転させない場合と同じ条件で解析することができる。

本発明のうち請求項３の発明によれば、連続体上のある点の解析範囲をその近傍の点の解析結果から定めることで、解析範囲を限定し、計算量を低減することができる。

本発明のうち請求項４の発明によれば、使用者は、撮像手段により対象物を撮像し、画像処理手段により画像データを処理することで、本発明の画像処理方法による測定を行うことができる。

本発明のうち請求項５の発明によれば、使用者は、汎用のパーソナルコンピュータ等によって、本発明の画像処理方法を実施することができる。

本発明は、対象物の変化前と変化後の二枚の画像の相関値を求め最大となった点を移動先とみなす画像相関法により、対象物の変位、変形および歪みを測定するものである。そこでまず、本発明において用いられる画像のアフィン変換について式により説明する。

アフィン変換による座標変換式は、変換前の座標を（ｘ，ｙ）、変換後の座標を（ｘ’，ｙ’）、回転中心の座標を（ｑ_ｘ，ｑ_ｙ）、回転角をθ、ｘ方向およびｙ方向の伸縮率を（ｓ_ｘ，ｓ_ｙ）としてとして式（１）で表される。
なお、このような座標変換を行う方法には、順変換とよばれる方法と逆変換とよばれる方法がある。式（１）で表される順変換は、変換前画像の座標（ｘ，ｙ）の光強度の数値を変換後の座標（ｘ’，ｙ’）に内挿する方法であるのに対し、以下の式（２）で表される逆変換は、変換後画像の座標（ｘ’，ｙ’）に内挿する光強度を変換前画像の座標（ｘ，ｙ）から参照する方法である。逆変換による方法では内挿する光強度を補間法により計算することができるため、ここでは逆変換を用いる。また、補間法としては最近隣内挿法、線形補間法、三次畳み込み内挿法等があるが、特に三次畳み込み内挿法が高精度であり適当である。

式（２）の逆行列を計算すると式（３）となる。
さらに式（３）を展開すると式（４）が求められる。

次に、本発明で用いる画像相関の手法について図および式により説明する。図１に示すのは、本発明の画像処理方法の概念図である。ここでは対象物Ｔとして、片持ち梁を想定する。

まず、変化前画像データ１上において、解析点３の座標と、解析点３を中心とする微小領域４の形状および大きさを定める。この微小領域４を変化後画像データ２上において１画素ずつシフトさせながら相関値を計算し、最大相関値が得られる点に解析点３が存在するとみなす。

微小領域４の大きさは、小さいほど空間分解能が高くなるが、測定の精度は低下しエラーが生じやすくなる。よってそのバランスを考慮して実施者が決定すべきものであるが、ここでは１５×１５画素とする。解析点３の変位がある程度予想できる場合は、変位する方向に短く、変位しない方向に長い形状としてもよい。

また、微小領域４をシフトさせる範囲（移動量解析範囲５）については、対象物Ｔの変位がまったく予想できない場合は変化後画像データ２上全域となるが、ある程度予想できる場合は、その予想範囲内でシフトさせればよい。例えば図１に示すような場合には、梁は必ず下方に変位すると考えられるので、移動量解析範囲５を解析点３よりも下方の範囲のみに限定してもよい。

ここで、ある点の相関値は、変化前後の画像において、当該点１画素とその近傍８画素の計９画素の差分をとり、その差分値によって求めることにする。式（５）において、変化前画像の点（ｘ，ｙ）における光強度の差分値をｇ（ｘ，ｙ）、変化後画像の点（ｘ，ｙ）における光強度の差分値をｇ’（ｘ，ｙ）として、相関値Ｓ（Δｘ，Δｙ）を求める。最大相関値を示すときの変位量がΔｘ，Δｙとなる。

なお、ここに挙げた手法は一例であり、画像データ相互の類似度を判断できる手法であればどのようなものであってもよい。

次に、微小領域４の回転について説明する。対象物Ｔの変位が平行移動又は微小な回転のみである場合には上記の方法のみで対応できるが、大きな回転変位を伴う場合、上記の方法では適切な点で最大相関値が得られない。そこで、微小領域４をアフィン変換により回転させて相関をとる必要がある。ここでは、式（４）において伸縮の項を省いた以下の式（６）を用いる。

式（６）により所定角度ずつ回転させたそれぞれの微小領域４を、ここではマスク８とよぶ。ここで、回転中心は微小領域４の中心であり、すなわち解析点３の座標となる。このように生成されるマスク８を、一枚ずつ変化後画像データ２上においてシフトさせながら相関値を計算し、最大相関値が得られるΔｘ，Δｙおよびθが、それぞれ解析点３の移動量および回転角とみなされる。

ここで、微小領域４を回転させる範囲（回転角解析範囲６）については、対象物Ｔの回転がまったく予想できない場合は0°〜360°の全域となるが、ある程度予想できる場合は、その予想範囲内で回転させればよい。また、初めから細かい角度毎（例えば、1°毎）のマスク８を生成してもよいが、粗い角度毎（例えば、5°毎）のマスク８を生成してそれぞれの相関値を計算し、最大相関値が得られた角度の近傍で細かい角度毎のマスク８を生成して再び相関値を計算すると、全体のマスク８生成数が減り、解析時間を短縮することができる。

なお、このように微小領域４を回転させる場合、微小領域４の形状が問題となる。すなわち、微小領域４の形状が正方形の場合、図２（ａ）に示すように、回転させることで変形後画像と相関をとるための計算点数が減少してしまい（図では25点→17点）、同じ条件を保つことができない。そこで、図２（ｂ）に示すように、微小領域４を円形とすれば、回転させても計算点数が減少しない。

次に、微小領域４の伸縮について説明する。上記の手法に加えて、さらに微小領域４の伸縮を考慮することにより、解析点３の歪みを求めることができる。

上記の手法により移動量および回転角が定められた微小領域４について、式（４）により所定伸縮率ずつ伸縮させたそれぞれの微小領域４を、ここでは二次マスク９とよぶ。このように生成される二次マスク９について、一枚ずつ変化後画像データ２との相関値を計算し、最大相関値が得られるｓ_ｘ，ｓ_ｙが、ｘ方向およびｙ方向の歪みとみなされる。

ここで、微小領域４を伸縮させる範囲（歪み解析範囲７）については、変化前画像データ１および変化後画像データ２を見比べて定めてもよいし、対象物Ｔの材質の歪みの限界値（破壊に至る直前の歪みの値）がわかる場合にはその値から定めてもよい。

以上より、解析点３の移動量、回転角および歪みが求められた。この一連の解析を対象物Ｔの全域に渡って行えば、対象物全体の変位、変形および歪みの状態を把握することができる。

なお、上記の手法においては、変化前画像データ１上において解析点３および微小領域４を定め、それを変化後画像データ２上においてシフトさせたが、これとは逆に、変化後画像データ２上において解析点３および微小領域４を定め、それを変化前画像データ１上においてシフトさせてもよい。前者の方法によれば、変化前の対象物の各点の移動先が把握できるので、ある物体に力を加えたときにどのように変形するか、といったことを検討する場合に好適である。一方、後者の方法によれば、変化後の対象物Ｔの各点の移動元が把握できるので、地滑りの際にどの部分がどのように変位したのか、といったことを検討する場合に好適である。

実際に本発明の画像処理方法を実施するには、以下に示す画像処理装置および画像処理プログラムを利用する。

本発明の画像処理プログラムは、キーボードやマウス等からなる入力装置、ディスプレイ等からなる出力装置、本プログラムの命令を順番に実行するＣＰＵ、本プログラムおよび本プログラムの実行に必要なデータおよび計算結果等を保存しておく記憶装置を構成要素とする標準的なコンピュータにより実行され、コンピュータが各種の手段として機能する。また、本発明の画像処理装置における撮像手段２００は、一般的なデジタルカメラからなり、画像処理手段１００は、前記プログラムが実行されるコンピュータにより構成される。図３に示すのは、本発明の画像処理装置の第一実施形態の構成を示すブロック図であり、図４に示すのは、本発明の画像処理装置の第一実施形態が実行する画像処理方法の流れを示すフローチャートである。

最初の過程は画像データ取得過程であり、撮像手段２００が、対象物の変化前画像および変化後画像を撮像し、それを画像データ取得手段１０１が、変化前画像データおよび変化後画像データとして取得する。画像データの解像度は、両者が同じであれば特に限定されないが、もちろん高解像度であるほうが解析の精度は高くなる。また、この段階で画像データ上の明らかに対象物でない部分（背景に相当する部分）を省いてしまってもよい。

次はパラメータ取得過程であり、パラメータ取得手段１０２が、使用者による入力手段３００からのパラメータの入力を受け付ける。入力手段３００は一般的なキーボードやマウス等の入力装置からなる。ここでは、解析点の座標および解析点を中心とする微小領域の大きさならびに解析点の移動量解析範囲（例えば、ｘ方向およびｙ方向に±80画素）、回転角解析範囲（例えば、±30°）および歪み解析範囲（例えば、ｘ方向およびｙ方向に0.9〜1.1）を設定し、さらに１画素に対応する実寸を設定する。解析点の座標は複数指定することもできる。また、微小領域の形状も指定することができ、形状に合わせて、矩形であれば縦横の長さ、円形であれば半径を指定する。

次はフィルタリング過程であり、フィルタリング手段１５１が、変化前画像データおよび変化後画像データに対してフィルタ処理を行い、画像を先鋭化する。フィルタには種々のものがあるが、ここではラプラシアンフィルタが好適である。このような処理を行うことで、画像のエッジが強調され、測定精度が向上する。

次は回転処理過程であり、回転処理手段１０３が、解析点を中心とした微小領域をアフィン変換により回転させたマスクを生成する。例えば、回転角解析範囲を±30°、回転単位を1°とすれば、60枚のマスクが生成される。

次は第一相関値算出過程であり、第一相関値算出手段１０４が、各マスクをｘ方向およびｙ方向へ１画素ずつシフトさせながら、変形後画像との相関値を計算する。シフトさせる範囲は、先に設定した移動量解析範囲であり、設定した範囲内におけるすべてのΔｘ，Δｙ，θの組み合わせについての相関値が求められる。

次は位置角度検出過程であり、位置角度検出手段１０５が、第一相関値算出過程において得られた相関値の中から最大相関値を検出し、最大相関値を与えるΔｘ，Δｙおよびθを解析点の移動量および回転角とみなす。

次は伸縮処理過程であり、伸縮処理手段１０６が、位置角度検出過程において最大相関値が得られたマスクについて、アフィン変換により伸縮させた二次マスクを生成する。例えば、伸縮率解析範囲をｘ方向およびｙ方向にそれぞれ0.9〜1.1、伸縮単位を0.01とすれば、440枚のマスクが生成される。

次は第二相関値算出過程であり、第二相関値算出手段１０７が、各マスクと変形後画像との相関値を計算する。設定した範囲内におけるすべてのｓ_ｘ，ｓ_ｙの組み合わせについての相関値が求められる。

次は伸縮率検出過程であり、伸縮率検出手段１０８が、第二相関値算出過程において得られた相関値の中から最大相関値を検出し、最大相関値を与えるｓ_ｘ，ｓ_ｙを解析点の歪みとみなす。

次は状態決定過程であり、状態決定手段１０９が、位置角度検出過程において最大相関値を与えるΔｘ，Δｙおよびθならびに伸縮率検出過程において最大相関値を与えるｓ_ｘ，ｓ_ｙを、解析点の移動量、回転角および歪みとみなし、変化後の解析点の状態を決定する。

以上の過程により、ある解析点についての移動量、回転角および歪みが求められる。解析点を複数指定した場合、その数の分だけ回転処理過程〜状態決定過程を繰り返す。

すべての解析点についての移動量、回転角および歪みが求められたあと、最後は出力過程であり、出力手段４００が、解析結果を出力・表示する。出力手段４００は一般的なディスプレイ等の出力装置からなる。出力の形式は、各解析点の移動量、回転角および歪みの数値データをリストにして表示するものや、変位をベクトルで表示するもの、等変位線を色の濃淡やカラーで表示するものなど、どのようなものであってもよい。

図５に示すのは、本発明の画像処理装置の第二実施形態の構成を示すブロック図であり、図６に示すのは、本発明の画像処理装置の第二実施形態が実行する画像処理方法の流れを示すフローチャートである。

ここで、第二実施形態の考え方を示す。解析範囲の設定に際して、第一実施形態においては、全解析点について同じ解析範囲で計算を行うか、解析点毎に解析範囲を設定するが、連続体の変位・変形について考える場合には、さらに効率的に解析範囲を定めることができる。すなわち、連続体上のある点の移動量、回転角および歪みの値は、その点の近傍の点の各値に極めて近いものと考えられる。そこで、ある点の解析範囲を定めるにあたり、その近傍の点ですでに解析した点の移動量、回転角および歪みの値を基準とした相対的な解析範囲を定めることができる。この方法によれば、連続体上の各解析点において相対的な解析範囲を狭くすることができるので、解析時間を短縮することができる。以下において、第一実施形態と異なる部分のみ説明する。

第一実施形態では、パラメータ取得過程において移動量解析範囲、回転角解析範囲および歪み解析範囲が設定されていたが、これらはこれから解析しようとする現解析点を基準に定められていた（例えば、「現解析点から±80画素」）。一方、第二実施形態では、解析範囲決定過程において、解析範囲決定手段１１０が、パラメータ取得手段１０２が受け付けた使用者の相対的解析範囲の入力値と、状態決定手段１０９が決定した現解析点の近傍の点ですでに解析した前解析点の移動量、回転角および歪みのデータとを読み込み、それらの値から解析範囲を設定する（例えば、使用者の相対的解析範囲の入力値が±5画素、前解析点の移動量が＋10画素であれば、「現解析点から＋5画素〜＋15画素」）。この方法による場合、現解析点の移動量、回転角および歪みは、前解析点のそれらと非常に近いと考えられるため、各解析範囲を狭く設定することができ、結果として計算量が少なくなる。

なお、画像相関法を用いる際には、対象物の表面に特徴ある模様や印が付いている必要がある。対象物表面が元々ランダムな模様であればそれを使えばよいが、模様がない場合や不適切な模様がある場合は、ランダムなドットを塗布するとよい。この際、ドットの数や分布および大きさ等が測定精度に影響するが、白地に黒の円形ドットを塗布する場合、同じ直径のドットよりランダムな直径のドットを塗布するほうが、測定精度がよい。また、黒いドット部分より白い地の部分の比率が高い方が、測定精度がよい。さらに、種々の色（例えば、赤、緑、青）のドットを塗布してもよい。

また、画像相関法の解析による変位量測定の最小単位は１画素であり、基本的にはそれ以下の感度で測定することができない。しかし、二次元サブピクセル処理によって１画素の1/10程度の感度で測定することができる。この方法は、最大相関点とその近傍８点の相関値による近似曲面を最小二乗法により作成し、真の最大相関位置を１画素以下の単位で求めることにより、１画素以下の変位を得る方法である。これにより、変位が微小な場合でも適切な結果が得られる。

上記実施形態ではプログラムによりコンピュータを各手段として機能させているが、各手段の一部又は全部をそれぞれ専用のハードウェアにより構成するものであってもよい。また、各過程の順序は、同等の結果を得られる範囲で上記実施形態と異なるものであってもよく、二以上の過程が並列的に実行されるものであってもよい。

以下において、本発明の画像処理方法により実際に物体の変形を測定した解析例を示す。測定対象は金属製の防振ダンパで、これを油圧シリンダで加力して変形させたものであり、解析にはその一部分を拡大した画像を用いた。図７に示すのが解析に用いた画像であり、（ａ）は変化前画像、（ｂ）は変化後画像である。画像中の防振ダンパの下部は固定されており、上部は右方へ変位している。

まず、比較対照のために、測定対象の並進運動のみを考慮し、回転を考慮しない手法による解析を行う。図８に示すのがその解析結果で、（ａ）は各点の移動量を色の濃淡で表したものであり（色が濃いほど移動量が大きい）、（ｂ）は各点の移動の様子をベクトルにより表したものである。（ａ）においては、所々に色の濃淡が不連続に変化している部分があり、正確な解析結果が得られていないことがわかる。また（ｂ）においても、所々に周囲と異なる向きのベクトルが存在する。

次に、本発明の手法による解析を行う。本発明の手法では、並進および回転を考慮した解析により測定対象の移動量および回転角が求められ、さらに伸縮を考慮した解析により測定対象の歪みが求められる。まず図９に示すのは移動量の解析結果で、（ａ）および（ｂ）が表すものは図８の場合と同様である。回転も考慮して解析すると、（ａ）においては色の濃淡の不連続な部分がほとんど見られず、正確な解析結果が得られていることがわかる。また（ｂ）においても、周囲と異なる向きのベクトルは見られない。

さらに、伸縮を考慮した解析について、ここでは解析点を図１０のように三点選んだ。解析点Ａは歪みが少ないと考えられる点、解析点Ｂは歪みがわずかに生じていると考えられる点、解析点Ｃは大きな歪みが生じていると考えられる点である。解析にあたっては、伸縮率解析範囲をｘ方向およびｙ方向にそれぞれ0.960〜1.040、伸縮単位を0.005とした。図１１の（ａ）、（ｂ）、（ｃ）はそれぞれ、Ａ、Ｂ、Ｃの各点における相関値の計算結果を示すものである。これらより最大相関値が求められ、それを与える伸縮率が各点の歪みとなる。

表１に、各点の座標と、移動量、回転角および歪みの解析結果を示す。座標は画像の左上を原点とし、右方向をｘ方向、下方向をｙ方向とした。表より、解析点Ａは歪みがほとんどなく、解析点Ｂはｙ方向に２％の歪みがあり、解析点Ｃはｘ方向およびｙ方向にそれぞれ２％の歪みがあることがわかる。

以上のように、本発明の手法によれば、測定対象の移動量および回転角を高い精度で求めることができる。さらに、従来の手法では不可能であった歪みの解析も可能であり、対象物の状態をより詳細に把握することができる。

本発明の画像処理方法において用いられる画像相関の概念図。 微小領域の形状と回転の関係を示す説明図。 本発明の画像処理装置の第一実施形態の構成を示すブロック図。 本発明の画像処理装置の第一実施形態が実行する画像処理方法の流れを示すフローチャート。 本発明の画像処理装置の第二実施形態の構成を示すブロック図。 本発明の画像処理装置の第二実施形態が実行する画像処理方法の流れを示すフローチャート。 解析に用いた実際の画像。 従来の手法による移動量の解析結果を示す図。 本発明の手法による移動量の解析結果を示す図。 解析点の位置を示す説明図。 歪みの解析における相関値の計算結果を示すグラフ。

符号の説明

１ 変化前画像データ
２ 変化後画像データ
３ 解析点
４ 微小領域
５ 移動量解析範囲
６ 回転角解析範囲
７ 歪み解析範囲
８ マスク
９ 二次マスク
Ｔ 対象物
１００ 画像処理手段
１０１ 画像データ取得手段
１０２ パラメータ取得手段
１０３ 回転処理手段
１０４ 第一相関値算出手段
１０５ 位置角度検出手段
１０６ 伸縮処理手段
１０７ 第二相関値算出手段
１０８ 伸縮率検出手段
１０９ 状態決定手段
１１０ 解析範囲決定手段
２００ 撮像手段

Claims (5)

1. 画像データ取得手段（１０１）が対象物（Ｔ）の変化前画像データ（１）および変化後画像データ（２）を受け付ける画像データ取得過程と、
   パラメータ取得手段（１０２）が変化前画像データ（１）上の解析点（３）の座標および解析点（３）を中心とする微小領域（４）の大きさならびに解析点（３）の移動量解析範囲（５）、回転角解析範囲（６）および歪み解析範囲（７）の数値データを受け付けるパラメータ取得過程と、
   回転処理手段（１０３）が微小領域（４）をアフィン変換により回転角解析範囲（６）において所定角度ずつ回転させたマスク（８）を生成する回転処理過程と、
   第一相関値算出手段（１０４）が変化後画像データ（２）上の移動量解析範囲（５）において各マスク（８）を走査しながら変化後画像データ（２）との相関値を計算する第一相関値算出過程と、
   位置角度検出手段（１０５）が最大相関値となる位置および角度を検出する位置角度検出過程と、
   伸縮処理手段（１０６）が最大相関値を示すマスク（８）をアフィン変換により歪み解析範囲（７）において所定伸縮率ずつ伸縮させた二次マスク（９）を生成する伸縮処理過程と、
   第二相関値算出手段（１０７）が各二次マスク（９）と変化後画像データ（２）との相関値を計算する第二相関値算出過程と、
   伸縮率検出手段（１０８）が最大相関値となる伸縮率を検出する伸縮率検出過程と、
   状態決定手段（１０９）が最大相関値を与える位置、角度および伸縮率から、それぞれ解析点の移動量、回転角および歪みを決定する状態決定過程と、
   を備えることを特徴とする画像処理方法。
2. 前記微小領域（４）の形状が、円形であることを特徴とする請求項１記載の画像処理方法。
3. 連続体上の複数点についての変位、変形および歪みの解析を行う場合において、
   解析範囲決定手段（１１０）が解析しようとする現解析点の移動量、回転角および歪みのうち少なくとも一つについて、現解析点の近傍の点ですでに解析した前解析点の移動量、回転角および歪みの値を基準とした相対的な解析範囲を定める解析範囲決定過程を備えることを特徴とする請求項１又は２記載の画像処理方法。
4. 変化前および変化後の対象物（Ｔ）を撮像し、変化前画像データ（１）および変化後画像データ（２）として出力する撮像手段（２００）と、
   請求項１、２又は３記載の画像処理方法を実行する画像処理手段（１００）と、
   を備えることを特徴とする画像処理装置。
5. コンピュータに、
   請求項１、２又は３記載の画像処理方法を実行させることを特徴とする画像処理プログラム。