(57)【要約】【課題】検査対象物の表面をより高精度に検査すること。【解決手段】実施形態の画像検査装置は、円周方向に所定周期で空間的に光の強度が変化するように照明部を制御する照明制御部と、前記照明部から出射された光の検査対象物での反射光を撮影して画像を取得するように撮像部を制御する撮影制御部と、空間的な光の強度の変化に応じた前記画像の輝度値と空間的な光の強度の変化に応じて前記所定周期で変化する複素参照信号との乗算値の、光の強度変化の円周方向一周分の積分値に基づく複素画像を算出する複素画像算出部と、前記複素画像から各画素の前記輝度値の位相を算出する位相算出部と、前記位相に基づいて前記検査対象物の表面を検査する検査部と、を備える。【選択図】図１
2018112478.tif 000002
【特許請求の範囲】
【請求項１】
円周方向に所定周期で空間的に光の強度が変化するように照明部を制御する照明制御部と、
前記照明部から出射された光の検査対象物での反射光を撮影して画像を取得するように撮像部を制御する撮影制御部と、
空間的な光の強度の変化に応じた前記画像の輝度値と空間的な光の強度の変化に応じて前記所定周期で変化する複素参照信号との乗算値の、光の強度変化の円周方向一周分の積分値に基づく複素画像を算出する複素画像算出部と、
前記複素画像から各画素の前記輝度値の位相を算出する位相算出部と、
前記位相に基づいて前記検査対象物の表面を検査する検査部と、
を備える画像検査装置。
【請求項２】
前記複素画像は複素位相画像を含む、請求項１に記載の画像検査装置。
【請求項３】
前記検査部は、前記位相に基づいて前記検査対象物の表面の異常の有無を判定する、請求項１または請求項２に記載の画像検査装置。
【請求項４】
前記検査部は、前記位相に基づいて、位相限定ラプラシアンを用いた処理を行うことによって、前記検査対象物の表面の異常の有無を判定する、請求項３に記載の画像検査装置。
【請求項５】
前記検査部は、位相限定ラプラシアンを用いた処理によって得たデータに対して、ソーベルフィルタを用いた処理を行うことによって、前記検査対象物の表面の異常の有無を判定する、請求項４に記載の画像検査装置。
【請求項６】
前記検査部による検査の結果を示す画像を表示するように表示部を制御する表示制御部を、さらに備える、請求項１〜５のいずれか一項に記載の画像検査装置。
【請求項７】
前記照明部は、全体が円環状に構成され、円周方向に複数の発光部が配置されており、
前記照明制御部は、前記複数の発光部のうち向かい合う発光部同士を円周方向に順番に点灯させるように前記照明部を制御し、
前記撮影制御部は、点灯する前記発光部が変わるたびに画像を取得するように前記撮像部を制御する、請求項１〜６のいずれか一項に記載の画像検査装置。
【請求項８】
円周方向に所定周期で空間的に光の強度が変化するように照明部を制御する照明制御ステップと、
前記照明部から出射された光の検査対象物での反射光を撮影して画像を取得するように撮像部を制御する撮影制御ステップと、
空間的な光の強度の変化に応じた前記画像の輝度値と空間的な光の強度の変化に応じて前記所定周期で変化する複素参照信号との乗算値の、光の強度変化の円周方向一周分の積分値に基づく複素画像を算出する複素画像算出ステップと、
前記複素画像から各画素の前記輝度値の位相を算出する位相算出ステップと、
前記位相に基づいて前記検査対象物の表面を検査する検査ステップと、
を含む画像検査方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、画像検査装置および画像検査方法に関する。
【背景技術】
【０００２】
従来、検査対象物に光を照射し、当該検査対象物の表面からの反射光を撮影して画像データとして取得し、当該画像データの輝度値の変化等に基づいて検査対象物の表面を検査する技術が、提案されている。
【先行技術文献】
【特許文献】
【０００３】
【特許文献１】特開２００８−１８５４７９号公報
【発明の概要】
【発明が解決しようとする課題】
【０００４】
この種の技術では、検査対象物の表面をより高精度に検査することができれば有意義である。
【課題を解決するための手段】
【０００５】
実施形態の画像検査装置は、円周方向に所定周期で空間的に光の強度が変化するように照明部を制御する照明制御部と、前記照明部から出射された光の検査対象物での反射光を撮影して画像を取得するように撮像部を制御する撮影制御部と、空間的な光の強度の変化に応じた前記画像の輝度値と空間的な光の強度の変化に応じて前記所定周期で変化する複素参照信号との乗算値の、光の強度変化の円周方向一周分の積分値に基づく複素画像を算出する複素画像算出部と、前記複素画像から各画素の前記輝度値の位相を算出する位相算出部と、前記位相に基づいて前記検査対象物の表面を検査する検査部と、を備える。
【図面の簡単な説明】
【０００６】
【図１】図１は、実施形態の画像検査システムの模式的な構成図である。
【図２】図２は、実施形態の画像検査システムにおける照明装置の構成を示した模式的な例示図である。
【図３】図３（ａ）は、実施形態の検査対象物の撮像画像である。図３（ｂ）は、図３（ａ）の撮像画像の一部を拡大した図である。
【図４】図４は、実施形態で用いる座標系である円座標の説明図である。
【図５】図５は、実施形態の検査対象物の一部に光を照射した様子を示す図である。
【図６】図６は、図５の点ａ、ｂ、ｃにおける輝度値の変化の様子を示す模式図である。
【図７】図７は、実施形態の画像検査装置によって得られた強度画像である。
【図８】図８は、座標にしたがって変位する位相を示す図である。
【図９】図９は、実施形態の画像検査装置によって得られた処理結果画像である。
【図１０】図１０は、実施形態の画像検査装置による処理を示すフローチャートである。
【図１１】図１１は、実施形態の画像検査装置によって得られたツールマーク方向画像である。
【発明を実施するための形態】
【０００７】
以下、本発明の例示的な実施形態が開示される。以下に示される実施形態の構成や制御、ならびに当該構成や制御によってもたらされる作用および結果（効果）は、一例である。本発明は、以下の実施形態に開示される構成や制御以外によっても実現可能であるとともに、基本的な構成や制御によって得られる種々の効果（派生的な効果も含む）を得ることが可能である。
【０００８】
まず、実施形態の画像検査システムの構成について説明する。図１は、実施形態の画像検査システムＳの模式的な構成図である。画像検査システムＳは、台５に載せられた検査対象物４（例えば金属製品）の被検査面４ａを検査する。画像検査システムＳは、被検査面４ａを撮影した画像（データ）に対する演算処理により、被検査面４ａにおける凹凸等の異常の有無を判定する。なお、以下において、画像とは、画像データを指し、具体的には、離散化された二次元の各点（各座標、各画素）での輝度値等のデータ群である。
【０００９】
ここで、本実施形態で例にとる検査対象物４について説明する。図３（ａ）は、実施形態の検査対象物４の撮像画像である。図３（ｂ）は、図３（ａ）の撮像画像の一部を拡大した図である。図３（ａ）に示すように、検査対象物４は、上面視で円形の被検査面４ａを有している。被検査面４ａの中央の黒い部分は、検査対象外の領域である。以下、被検査面４ａとは、その検査対象外の領域以外を指すものとする。
【００１０】
被検査面４ａには、同心円状の多数のツールマーク（例えば旋盤等の加工機による切削痕）、および、ツールマークと異なる方向のキズ（図３（ｂ）の領域１００１内の右上がりの白い部分）がある。本実施形態では、当該キズを検出する場合を例にとる。なお、キズの方向がツールマークの方向と同じでないことを前提とする。また、キズは、引っ掻いたことにより生じたものでもよいし、あるいは、打痕でもよい。
【００１１】
図１に戻って、画像検査システムＳは、画像検査装置１と、撮像部２と、照明装置３と、を備える。照明装置３について、図２も合わせて参照して説明する。図２（図１のＡ−Ａ矢視方向で見た図）は、実施形態の画像検査システムＳにおける照明装置３の構成を示した模式的な例示図である。照明装置３は、全体が円環状に構成されており、１２個の発光部Ｌ１〜Ｌ１２（以下、区別しないときは「発光部Ｌ」という。）を有する。照明装置３は、中心周りに開口部が設けられた円環状に構成されている。
【００１２】
図２に示すように、照明装置３は、１２個の発光部（照明領域）が円環状に配置されることで構成されている。各発光部Ｌ１〜Ｌ１２は、照明装置３における検査対象物４に面した部分が、円環の中心に対して等角度間隔（３０度間隔）で１２等分された領域に対応する。発光部Ｌ１〜Ｌ１２の各々につき、同じ数（１つまたは複数）の光源が設けられている。これにより、実施形態では、各発光部Ｌ１〜Ｌ１２が、それぞれ同一の光量の光を照射することが可能なように構成されている。これにより、発光部Ｌは、検査対象物４の方向に光を出射できる。
【００１３】
また、各発光部Ｌは、画像検査装置１の制御部１４の制御にしたがって、個別に点灯消灯できるようになっている。
【００１４】
各発光部Ｌが光を出射すると、検査対象物４における対応する部分に光が照射される。例えば、発光部Ｌ２、Ｌ３、Ｌ４、Ｌ８、Ｌ９、Ｌ１０が光を出射すると、図５に示すように、検査対象物４におけるそれらに対応する部分に光が照射される。各発光部Ｌは、例えば、複数のＬＥＤ（Light  Emitting  Diode）により構成される。なお、各発光部Ｌは、ＯＥＬＤ（Organic  Electro  Luminescent  Display）のシートによって構成されてもよい。
【００１５】
図１に戻って、撮像部２は、照明装置３の中央の開口部を通して検査対象物４に臨み、照明装置３から出射され被検査面４ａで反射した光を撮影する。撮像部２は、例えば、ディジタルカメラ（ディジタルスチルカメラ）であり、ディジタル画像、すなわち２次元の画素毎の輝度値（データ）を取得する。撮像部２は、例えば、ＣＭＯＳ（Complementary  Metal-Oxide-Semiconductor）イメージセンサやＣＣＤ（Charge-Coupled  Device）イメージセンサ等の撮像素子を有する。
【００１６】
次に、画像検査装置１について説明する。画像検査装置１は、コンピュータ装置であり、記憶部１１、入力部１２、表示部１３、制御部１４を有する。画像検査装置１は、外部装置（撮像部２、照明装置３等）との通信インタフェースも備えるが、図１での図示や説明を省略する。
【００１７】
記憶部１１は、情報を記憶する手段であり、例えば、ＲＡＭ（Random  Access  Memory）、ＲＯＭ（Read  Only  Memory）、ＨＤＤ（Hard  Disk  Drive）、ＳＳＤ（Solid  State  Drive）等である。入力部１２は、ユーザが情報を入力する手段であり、例えば、キーボードやマウスである。表示部１３は、情報を表示する手段であり、例えば、ＬＣＤ（Liquid  Crystal  Display）である。
【００１８】
制御部１４は、例えばＣＰＵ（Central  Processing  Unit）や制御コントローラ等として構成される。制御部１４は、照明制御部１４１、撮影制御部１４２、撮像画像取得部１４３、強度算出部１４４、複素復調演算部１４５、位相算出部１４６、振幅算出部１４７、異常判定部１４８、表示制御部１４９を備える。
【００１９】
照明制御部１４１は、円周方向に所定周期で空間的に光の強度が変化するように照明装置３（照明部）を制御する。照明制御部１４１は、例えば、照明装置３における１２個の発光部Ｌ（図２）のうち６個ずつを次の点灯順１〜６で点灯させる。
（点灯順１）発光部Ｌ１２、Ｌ１、Ｌ２、Ｌ６、Ｌ７、Ｌ８
（点灯順２）発光部Ｌ１、Ｌ２、Ｌ３、Ｌ７、Ｌ８、Ｌ９
（点灯順３）発光部Ｌ２、Ｌ３、Ｌ４、Ｌ８、Ｌ９、Ｌ１０（図５）
（点灯順４）発光部Ｌ３、Ｌ４、Ｌ５、Ｌ９、Ｌ１０、Ｌ１１
（点灯順５）発光部Ｌ４、Ｌ５、Ｌ６、Ｌ１０、Ｌ１１、Ｌ１２
（点灯順６）発光部Ｌ５、Ｌ６、Ｌ７、Ｌ１１、Ｌ１２、Ｌ１
【００２０】
撮影制御部１４２は、照明装置３から出射された光の検査対象物４での反射光を撮影して画像を取得するように撮像部２を制御する。具体的には、撮影制御部１４２は、照明制御部１４１が、上述の点灯順で発光部Ｌを点灯させた場合ごとに撮影を行うように撮像部２を制御する。なお、本実施形態の照明制御部１４１は、時間の遷移にしたがって、円周方向で、明暗が変化するように、点灯制御を行うこととした。
【００２１】
ここで、実施形態で用いる座標系について説明する。図４は、実施形態で用いる円座標の説明図である。本実施形態では、座標系として、動径方向をＹとし、角度方向をＸとする円座標を用いる。
【００２２】
次に、検査対象物４の撮像画像における各点（位置。座標）における輝度値の変化について説明する。図５は、実施形態の検査対象物４の一部に光を照射した様子を示す図である。図６は、図５の点ａ、ｂ、ｃにおける輝度値の変化の様子を示す模式図である。
【００２３】
図５に示すように点ａ、ｂ、ｃをとったものとする。その場合、照明装置３における発光部Ｌを点灯順１〜６で点灯させると、点ａ、ｂ、ｃそれぞれの輝度値は、図６に示すように変化する。例えば、点ａの輝度値の変化は、点灯順２で最高となり点灯順５で最低となるような正弦曲線（サインカーブ。余弦曲線と同様）として近似できる。また、点ｂの輝度値の変化は、点灯順３で最高となり点灯順６で最低となるような正弦曲線として近似できる。また、点ｃの輝度値の変化は、点灯順４で最高となり点灯順１で最低となるような正弦曲線として近似できる。
【００２４】
なお、ツールマークとキズはいずれも被検査面４ａにおける溝（凹凸）という点で共通するが、ツールマークの方向は緩やかにしか変化しないのに対して、キズの方向はツールマークの方向と異なっている。したがって、画素単位で見た場合、キズがある部分の輝度値の変化の正弦曲線は、キズがない部分の輝度値の変化の正弦曲線に比べて、位相と、振幅の点で差が生じる。
【００２５】
また、ツールマークの深さ（凹凸の大きさ）は、加工機ごとに異なり、さらに、同じ加工機であっても刃先状態（刃先の鋭利さ）によって異なる。そして、ツールマークの深さは、輝度値の変化の正弦曲線の振幅に影響するが、位相には影響しない。
【００２６】
したがって、画素単位で輝度値の変化の正弦曲線の位相が局所的に変化している部分を特定することで、ツールマークの深さに関係なく（ツールマークがない場合も含む。）、キズを検出することができる。以下、そのための処理について詳述する。なお、キズを検出する場合、強度、位相、振幅のいずれか、あるいはそれらの２つ以上の組み合わせによって検出処理を構築することで、検出性能向上や過検出低減を図ることができる。本実施形態では、一例として、位相を使用する場合について説明する。撮像画像取得部１４３〜表示制御部１４９の説明は、後述する。
【００２７】
まず、照明装置３を点灯順１〜６で点灯させながら撮影された画像データに基づく画像処理について説明する。上述のように、照明装置３を点灯順１〜６で点灯させると、撮像画像における各点における輝度値の変化は、正弦曲線として近似できる。
【００２８】
上述したように、本実施形態では、照明装置３の光の強度（輝度値）は、円周方向に所定周期で空間的に変化している。このように空間的に光の周期的な強度変化が与えられた環境において、本実施形態では、点灯順ごとに撮像部２による撮影が実行される。
【００２９】
また、座標（ｍ，ｎ）の各画素における強度Ｉ（ｍ，ｎ）は、以下の式（１）で定義される。
【数１】
2018112478.tif 000003
【００３０】
ここで、「ｍ」および「ｎ」は、それぞれ、２次元の画像中に配置された画素の位置座標を示すパラメータである。Ａ＋Ｂｃｏｓ（ωｘ＋φ）は、各画素の輝度値（関数）である。φは、位相の遅れである。強度Ｉ（ｍ，ｎ）による画像を強度画像と呼ぶ。図７は、実施形態の画像検査装置１によって得られた強度画像である。図７に示す強度画像では、図３と同様、同心円状の多数のツールマーク、および、ツールマークと異なる方向のキズ（領域１００２内の右上がりの白い部分）がある。
【００３１】
また、各画素における複素画像Ｆ（ｍ，ｎ）は、以下の式（２）で定義される。
【数２】
2018112478.tif 000004
【００３２】
ここで、「ｊ」は、虚数単位（√（−１））である。「ｅ−ｊωｘ」は、複素参照信号（複素指数関数）である。「ωｘ」は、複素参照信号の位相である。αは、任意に設定される分割数で、本実施形態では発光部Ｌによる点灯順１〜６の数に対応して「６」である。Ｔは、分割数αに応じて定まる増分である。なお、Ｔは、露光時間に比例する。
【００３３】
オイラーの公式により、複素参照信号「ｅ−ｊωｘ」の実部はｃｏｓωｘで、虚部は−ｓｉｎωｘである。複素画像Ｆ（ｍ，ｎ）は、各画素の輝度値（Ａ＋Ｂｃｏｓ（ωｘ＋φ））と複素参照信号ｅ−ｊωｘとの乗算値の、複素参照信号の一周期分の積分値に相当し、輝度値（Ａ＋Ｂｃｏｓ（ωｘ＋φ））の複素直交復調信号や、複素直交検波信号等とも称されうる。
【００３４】
式（２）において、各画素における複素画像Ｆ（ｍ，ｎ）の実部Ｆｒ（ｍ，ｎ）は、式（３）となり、虚部Ｆｉ（ｍ，ｎ）は、式（４）となる。
【数３】
2018112478.tif 000005
【数４】
2018112478.tif 000006
【００３５】
以下、各画素における複素画像Ｆ（ｍ，ｎ）に関して、虚部と実部の２乗和の平方根を「振幅」と称する。各画素における複素画像Ｆ（ｍ，ｎ）の振幅Ｓ（ｍ，ｎ）は、式（５）と定義され、位相Ｐ（ｍ，ｎ）は、式（６）と定義される。なお、振幅Ｓ（ｍ，ｎ）は、複素振幅画像とも称され、位相Ｐ（ｍ，ｎ）は、複素位相画像とも称されうる。
【数５】
2018112478.tif 000007
【数６】
2018112478.tif 000008
【００３６】
ここで、式（３），（４）から、以下の式（７），（８）が得られる。
【数７】
2018112478.tif 000009
【数８】
2018112478.tif 000010
【００３７】
よって、式（５），（６）から、以下の式（９），（１０）が得られる。
【数９】
2018112478.tif 000011
【数１０】
2018112478.tif 000012
【００３８】
式（１０）から、複素画像Ｆ（ｍ，ｎ）の位相Ｐ（ｍ，ｎ）が輝度値Ｌ（ｍ，ｎ）の位相ωｘの遅れφに対応していることがわかる。つまり、複素画像Ｆ（ｍ，ｎ）を算出し、さらにその位相Ｐ（ｍ，ｎ）を得ることにより、輝度値Ｌ（ｍ，ｎ）の位相ωｘの遅れφを把握することができる。以下、位相Ｐ（ｍ，ｎ）を位相画像データともいう。
【００３９】
また、本実施形態では直接使わないが、式（９）から、複素画像Ｆ（ｍ，ｎ）の振幅Ｓ（ｍ，ｎ）が輝度値Ｌ（ｍ，ｎ）の振幅Ｂに対応していることがわかる。つまり、複素画像Ｆ（ｍ，ｎ）を算出し、さらにその振幅Ｓ（ｍ，ｎ）を得ることにより、輝度値Ｌ（ｍ，ｎ）の振幅Ｂを把握することができる。
【００４０】
次に、位相限定ラプラシアンを用いて、位相の分布の勾配に関する特徴を算出することで異常を検出する方法について説明する。
【００４１】
上記の式（１０）からわかるように、位相画像データの各画素値（Ｐ（ｍ，ｎ））は、−π〜πの範囲に折りたたまれる。したがって、位相画像データの各画素値は、−πからπに、またはπから−πに不連続に変化しうる（位相ジャンプ）。このような位相ジャンプは、被検査面４ａが、凹凸などの局所的な異常部分（キズ）を含まない平坦になっている場合でも、位相画像データ上にエッジとして現れる。
【００４２】
キズ等の異常部分における画素値の変化は、エッジにおける画素値の変化よりも小さい。したがって、位相画像データに対して閾値などを用いた通常の検出処理を行うだけでは、エッジがノイズとなって、異常部分を検出することが容易でない。
【００４３】
上述したように、ツールマークの方向は緩やかにしか変化しないので、ツールマークしかない部分では、隣接する画素ごとの位相の変化は一次関数的に緩やかにしか変化しない。ここで、図８は、座標にしたがって変位する位相を示す図である。ツールマークと異なる方向のキズがある部分では、図８の領域１１００で示すように、位相が急峻に変化する。したがって、定常的な位相の変化を無視することができれば、位相が急峻に変化する異常部分（キズ）のみを容易に検出することができる。
【００４４】
そこで、本実施形態においては、位相限定ラプラシアンを用いた処理を行うことにより、上記の位相ジャンプと緩やかな変化の影響を無視した位相の分布の勾配に関する特徴を算出する。ここで、位相限定ラプラシアンとは、振幅および位相を含む複素数で表現される画像データ（以下、単に「画像データ」という。）のうち、振幅部分を無視した位相部分にのみ２階微分を施すための演算式であり、位相の定常的な変化を無視し、位相の急峻な変化を検出するための演算式である。
【００４５】
一例として、以下に示すラプラシアンフィルタ（Ａ）に対応する位相限定ラプラシアンについて説明する。
【表１】
2018112478.tif 000013
【００４６】
このラプラシアンフィルタ（Ａ）は、通常の実数の画素値を有する画像データにおいて、処理対象の画素の画素値と、処理対象の画素の周囲に隣接する８個の画素の画素値との差分をとるための、いわゆる８近傍ラプラシアンフィルタである。より具体的には、このラプラシアンフィルタ（Ａ）は、通常の実数の画素値を有する画像データにおいて、処理対象の画素の画素値を８倍した値から、処理対象の画素の周囲に隣接する８個の画素の画素値の総和を引いた値を、処理対象の画素の新たな画素値とするためのフィルタである。
【００４７】
撮像部２によって得られる画像データの各画素値は、振幅および位相を含む複素指数関数で表現される。したがって、画像データの位相部分にこのラプラシアンフィルタ（Ａ）に対応する処理を施したい場合、まず、処理対象の画素の画素値（ｇω（ｉ，ｊ）とする）の位相を８倍するために、処理対象の画素の画素値を８乗する必要がある。
【００４８】
そして、処理対象の画素の周囲に隣接する８個の画素の画素値ｇω（ｉ−１，ｊ−１），ｇω（ｉ−１，ｊ），ｇω（ｉ−１，ｊ＋１），ｇω（ｉ，ｊ−１），ｇω（ｉ，ｊ＋１），ｇω（ｉ＋１，ｊ−１），ｇω（ｉ＋１，ｊ），ｇω（ｉ＋１，ｊ＋１）の位相の総和を求めるために、これら８個の画素の画素値を全て掛け合わせる必要がある。
【００４９】
そして、処理対象の画素の画素値の位相を８倍にしたものと、処理対象の画素の周囲に隣接する８個の画素の位相の総和と、の差を求めるために、処理対象の画素の画素値を８乗したものを、処理対象の画素の周囲に隣接する８個の画素の画素値を掛け合わせたもので割る必要がある。
【００５０】
上記の３つの演算を式で表すと、下記の式（１１）〜（１３）のようになる。
【００５１】
【数１１】
2018112478.tif 000014
【数１２】
2018112478.tif 000015
【数１３】
2018112478.tif 000016
【００５２】
ここで、式（１１）においてｇω（ｉ，ｊ）を８乗ではなく９乗している理由は、式（１３）における乗算にｇω（ｉ，ｊ）が含まれているからである。これにより、式（１１）の値を式（１２）の値で割った値の指数部分が、処理対象の画素の画素値ｇω（ｉ，ｊ）の位相を８倍にしたものと、処理対象の画素の周囲に隣接する８個の画素の位相の総和と、の差に対応する。なお、式（１３）では、式（１１）の値を式（１２）の値で割った後、位相のみ（ａｒｇ）を取り出している。
【００５３】
以上の式（１１）〜（１３）で示した位相限定ラプラシアンを用いた処理を、画像データの全画素について実行することにより、画像データから、２階微分が施された位相画像データを算出することができる。
【００５４】
なお、ラプラシアンフィルタは、ほかに、処理対象の画素と、処理対象の画素の上下および左右に隣接する４個の画素との差分をとるための、いわゆる４近傍ラプラシアンフィルタであってもよい。また、３×３のラプラシアンフィルタとは異なる、５×５のラプラシアンフィルタであってもよい。一般的なラプラシアンフィルタに対応する位相限定ラプラシアンは、下記の式（１４）〜（１６）で表される。
【００５５】
【数１４】
2018112478.tif 000017
【数１５】
2018112478.tif 000018
【数１６】
2018112478.tif 000019
【００５６】
式（１４）は、フィルタ係数が正となる画素の集合Ｐのそれぞれの画素値をフィルタ係数で重み付けしたものを掛け合わせる、すなわちフィルタ係数で重み付けした位相の総和をとるための式である。この式（１４）において、Ｍ（ｋ，ｌ）は、処理対象の画素の位置を（ｉ，ｊ）とした場合における、位置（ｉ＋ｋ，ｊ＋ｌ）の画素に対応するフィルタ係数（正の値）である。
【００５７】
また、式（１５）は、フィルタ係数が負となる画素の集合Ｎのそれぞれの画素値をフィルタ係数で重み付けしたものを掛け合わせる、すなわちフィルタ係数で重み付けした位相の総和をとるための式である。この式（１５）において、Ｎ（ｋ，ｌ）は、処理対象の画素の位置を（ｉ，ｊ）とした場合における、位置（ｉ＋ｋ，ｊ＋ｌ）の画素に対応するフィルタ係数（負の値）の符号を反転させたものである。
【００５８】
そして、式（１６）は、式（１４）の値を式（１５）の値で割ることにより、フィルタ係数が正となる画素の集合Ｐのフィルタ係数で重み付けした位相の総和と、フィルタ係数が負となる画素の集合Ｎのフィルタ係数で重み付けした位相の総和との差分をとり、その差分を取り出すための式である。
【００５９】
このようにして、位相画像データに対して位相限定ラプラシアンを用いることで、位相が急峻に変化する異常部分（キズ）のみを容易に検出することができる。
【００６０】
なお、位相限定ラプラシアンを用いた処理によって得たデータに対して、ソーベルフィルタを用いることで、異常部分をより強調することができる。ソーベルフィルタとは、１次微分を計算し、輪郭を検出するフィルタである。例えば、上方向のソーベルフィルタは、次の通りである。
【表２】
2018112478.tif 000020
【００６１】
同様にして、下方向、右方向、左方向、右下方向、右上方向、左下方向、左上方向のソーベルフィルタを用いることができる。また、上下左右の４方向同時のソーベルフィルタや、８方向同時のソーベルフィルタを用いてもよい。
【００６２】
位相画像データに対して位相限定ラプラシアンとソーベルフィルタによる処理を行うことで、図９のような画像を得ることができる。図９は、実施形態の画像検査装置１によって得られた処理結果画像である。図９の処理結果画像では、図３や図７と比較してわかるように、同心円状の多数のツールマークは現れず、ツールマークと異なる方向のキズ（領域１００３内の右上がりの白い部分）だけが現れている。
【００６３】
図１に戻って、画像検査装置１の制御部１４において、撮像画像取得部１４３は、撮像部２で撮影した画像から、各位置における輝度値Ｌ（ｍ，ｎ）を取得し、記憶部１１に記憶する。
【００６４】
強度算出部１４４は、輝度値Ｌ（ｍ，ｎ）から強度Ｉ（ｍ，ｎ）を算出（設定）する。また、複素復調演算部１４５（複素画像算出部）は、輝度値Ｌ（ｍ，ｎ）および複素参照信号に基づいて、複素画像Ｆ（ｍ，ｎ）を算出する。より詳細には、複素復調演算部１４５は、空間的な光の強度の変化に応じた画像の輝度値と空間的な光の強度の変化に応じて所定周期で変化する複素参照信号との乗算値の、光の強度変化の円周方向一周分の積分値に基づく複素画像を算出する。複素参照信号は、記憶部１１に記憶されているか、あるいは記憶部１１に記憶されたデータに基づいて算出される。位相算出部１４６は、複素画像Ｆ（ｍ，ｎ）から位相Ｐ（ｍ，ｎ）を算出する。振幅算出部１４７は、複素画像Ｆ（ｍ，ｎ）から振幅Ｓ（ｍ，ｎ）を算出する。
【００６５】
異常判定部１４８（検査部）は、位相Ｐ（ｍ，ｎ）に基づいて検査対象物４の表面を検査する。例えば、異常判定部１４８（検査部）は、位相Ｐ（ｍ，ｎ）に基づいて検査対象物４の表面の異常の有無を判定する。その際、例えば、異常判定部１４８（検査部）は、は、位相Ｐ（ｍ，ｎ）に基づいて、位相限定ラプラシアンを用いた処理を行うことによって、検査対象物４の表面の異常の有無を判定する。また、異常判定部１４８（検査部）は、位相限定ラプラシアンを用いて得たデータに対して、ソーベルフィルタを用いた処理を行うことによって、検査対象物４の表面の異常の有無を判定する。表示制御部１４９は、異常判定部１４８による異常の判定の結果を示す画像を表示するよう、表示部１３を制御する。
【００６６】
なお、制御部１４による上述した演算処理や制御は、ソフトウエアによって実行されてもよいし、ハードウエアによって実行されてもよいし、それらの併用により実行されてもよい。ソフトウエアによる処理の場合にあっては、制御部１４は、ＲＯＭ、ＨＤＤ、ＳＳＤ等に記憶されたプログラムを読み出して実行する。この場合、プログラムには、制御部１４に含まれる各部に対応するモジュールが含まれる。制御部１４の全部あるいは一部がハードウエアによって構成される場合、制御部１４には、ＦＰＧＡ（Field  Programmable  Gate  Array）やＡＳＩＣ（Application  Specific  Integrated  Circuit）等が含まれる。
【００６７】
図１０は、実施形態の画像検査装置１による処理を示すフローチャートである。画像検査装置１による検査は、図１０に示されるような手順で実行される。まず、照明制御部１４１は、照明装置３における所定の発光部Ｌ（例えば、点灯順１に対応する各発光部Ｌ）が点灯するように、照明装置３を制御する（ステップＳ１）。これを受けて、照明装置３における所定の発光部Ｌが光を出射する。
【００６８】
撮影制御部１４２は、撮影を行うように撮像部２を制御する（ステップＳ２）。これを受けて、撮像部２は、照明装置３から出射され被検査面４ａで反射した光を撮影する。次に、撮像画像取得部１４３は、撮像部２で撮影した画像から、各点における輝度値Ｌ（ｍ，ｎ）（画像データ）を取得し、記憶部１１に記憶する（ステップＳ３）。
【００６９】
１周期分の処理が終了していない場合（ステップＳ４でＮｏ）、照明制御部１４１は、照明装置３における次の発光部Ｌ（例えば、点灯順２に対応する各発光部Ｌ）が点灯するように、照明装置３を制御する（ステップＳ５）。これを受けて、照明装置３における次の発光部Ｌが光を出射する。ステップＳ５の後、ステップＳ２に戻る。
【００７０】
１周期分の処理が終了した場合（ステップＳ４でＹｅｓ）、ステップＳ６に進む。１周期分の処理が終了した場合とは、例えば、点灯順１〜６における撮影が終了した場合である。
【００７１】
次に、強度算出部１４４は、輝度値Ｌ（ｍ，ｎ）から強度Ｉ（ｍ，ｎ）を算出（設定）する（ステップＳ６）。次に、複素復調演算部１４５は、輝度値Ｌ（ｍ，ｎ）および複素参照信号に基づいて、複素画像Ｆ（ｍ，ｎ）を算出する（ステップＳ７）。次に、位相算出部１４６は、複素画像Ｆ（ｍ，ｎ）から位相Ｐ（ｍ，ｎ）を算出する（ステップＳ８）。
【００７２】
次に、振幅算出部１４７は、複素画像Ｆ（ｍ，ｎ）から振幅Ｓ（ｍ，ｎ）を算出する（ステップＳ９）。次に、異常判定部１４８は、位相Ｐ（ｍ，ｎ）に基づいて、被検査面４ａにおける異常の有無を判定する（ステップＳ１０）。異常判定部１４８は、例えば、上述した位相限定ラプラシアンとソーベルフィルタによる処理を用いて異常を判定する。次に、表示制御部１４９は、異常判定部１４８による異常の判定結果を示す画像を表示するよう、表示部１３を制御する（ステップＳ１１）。
【００７３】
このように、本実施形態の画像検査装置１によれば、照明装置３から出射する光の強度を円周方向に所定周期で空間的に変化させ、空間的に光の強度が変化するたびに撮像部２で取得した画像の輝度値と複素参照信号との乗算値の積分値に基づく複素画像を算出し、複素画像から各画素の輝度値の位相を算出し、その位相に基づいて検査対象物４の表面をより高精度に検査することができる。
【００７４】
例えば、上述したように、複素画像に含まれる複素位相画像（Ｐ（ｍ，ｎ））を得ることにより、輝度値Ｌ（ｍ，ｎ）の位相ωｘの遅れφを把握することができ、それによって検査対象物４の表面をより高精度に検査することができる。
【００７５】
また、異常判定部１４８（検査部）は、輝度値の位相に基づいて、位相限定ラプラシアンを用いた処理を行うことによって、位相の定常的な変化を無視し、位相の急峻な変化を異常として検出することができる。
【００７６】
また、異常判定部１４８（検査部）は、位相限定ラプラシアンを用いた処理によって得たデータに対して、ソーベルフィルタを用いた処理を行うことによって、異常部分（キズ）をより強調することができる。
【００７７】
また、異常判定部１４８（検査部）による検査の結果を示す画像を表示部１３に表示することで、ユーザは、キズ等を視認することができる。
【００７８】
また、全体が円環状で１２個の発光部Ｌを有し、その１２個の発光部Ｌの点灯消灯を個別に制御できる照明装置３を用いて、向かい合う発光部Ｌ同士を円周方向に順番（例えば点灯順１〜６の順番）に点灯させることで、検査対象物４の被検査面４ａにおける各点に対し、光の入射角や入射方向、光源からの距離等の点でバランスよく光を照射することができ、良好な撮像画像を得ることができる。
【００７９】
従来技術では、検査対象物にツールマークとキズの両方がある場合、例えば、両者の輝度差によってキズのみを検出しようとするが、その方法では、両者にある程度以上の輝度差がないとキズを高精度で検出できない。一方、本実施形態の画像検査装置１では、両者の輝度差がなくても、ツールマークの方向は緩やかにしか変化しないことを利用して、ツールマークと異なる方向のキズを高精度で検出できる。
【００８０】
なお、画像検査装置１を用いて、ツールマークの方向を表すツールマーク方向画像を作成し、表示することもできる。図１１は、実施形態の画像検査装置１によって得られたツールマーク方向画像である。光の照射方向の水平方向成分（例えば、図２の発光部Ｌ１の場合、真下よりやや左の方向）と平行な方向のツールマークは明暗の差が小さく、光の照射方向の水平方向成分と垂直な方向のツールマークは明暗の差が大きいので、ツールマークの方向は、位相の遅れφに９０度を足したものとして定義できる。図１１のツールマーク方向画像では、ツールマーク（およびキズ）の方向が同じものには同じ色を与えるようにして表示している。そして、図１１の領域１００４には、ツールマークと異なる方向のキズ（右上がりの白い部分）が現れている。
【００８１】
なお、上述した実施形態の画像検査装置１で実行される検査プログラムは、インストール可能な形式又は実行可能な形式のファイルでＣＤ−ＲＯＭ、フレキシブルディスク（ＦＤ）、ＣＤ−Ｒ、ＤＶＤ（Digital  Versatile  Disk）等のコンピュータで読み取り可能な記録媒体に記録されて提供可能である。
【００８２】
また、上述した実施形態の画像検査装置１で実行される検査プログラムを、インターネット等のネットワークに接続されたコンピュータ上に格納し、ネットワーク経由でダウンロードさせることにより提供するように構成してもよい。また、上述した実施形態の画像検査装置１で実行される検査プログラムをインターネット等のネットワーク経由で提供または配布するように構成してもよい。
【００８３】
本発明のいくつかの実施形態及び変形例を説明したが、これらの実施形態及び変形例は、例として提示したものであり、発明の範囲を限定することは意図していない。これら新規な実施形態及び変形例は、その他の様々な形態で実施されることが可能であり、発明の要旨を逸脱しない範囲で、種々の省略、置き換え、変更を行うことができる。これら実施形態やその変形例は、発明の範囲や要旨に含まれるとともに、特許請求の範囲に記載された発明とその均等の範囲に含まれる。
【００８４】
例えば、上述の実施形態では、ツールマークがある場合のキズの検出について説明したが、ツールマークがなくても、同様の処理でキズを検出することができる。
【００８５】
また、ツールマークの形状は、同心円状に限定されず、線状、うろこ状などであってもよい。ただし、ツールマーク同士が交差していないことが好ましい。
【００８６】
また、上述したように、ツールマークの深さ（凹凸の大きさ）は輝度値の変化の正弦曲線の振幅に影響する。したがって、画像検査装置１は、振幅を用いて、キズなどの異常判定ではなく、検査対象物４を加工した加工機の刃先状態（刃先の鋭利さ）を評価することもできる。
【００８７】
また、上述の実施形態では、画像検査装置１において、制御部１４と撮像部２を別々に設けたが、制御部１４と撮像部２を別々に設ける必要は無く、制御部１４の少なくとも一部が撮像部２と一体に設けられたカメラを有する構成としてもよい。
【００８８】
また、照明装置３の１２個の発光部Ｌのうち６個（向かい合う３個×２）ずつを所定の点灯順で点灯させるようにしたが、これに限定されず、４個（向かい合う２個×２）ずつや、８個（向かい合う４個×２）ずつを、所定の点灯順で点灯させるようにしてもよい。
【００８９】
また、照明装置３が有する発光部Ｌの数は、１２個に限定されず、１８個、２４個等、別の個数であってもよい。
【００９０】
また、検査対象物４は、金属製品に限定されず、樹脂製品や、塗装を施した製品等、光をある程度反射するものであれば、何でもよい。
【００９１】
また、検査対象物４の表面に段差がある場合でも、段差の影響を受けずにキズを検出できることが、実験で判明している。
【符号の説明】
【００９２】
１…画像検査装置、２…撮像部、３…照明装置、４…検査対象物、４ａ…被検査面、５…台、１１…記憶部、１２…入力部、１３…表示部、１４…制御部、１４１…照明制御部、１４２…撮影制御部、１４３…撮像画像取得部、１４４…強度算出部、１４５…複素復調演算部（複素画像算出部）、１４６…位相算出部、１４７…振幅算出部、１４８…異常判定部（検査部）、１４９…表示制御部、Ｓ…画像検査システム。
【図１】
2018112478.tif 000021
【図２】
2018112478.tif 000022
【図３】
2018112478.tif 000023
【図４】
2018112478.tif 000024
【図５】
2018112478.tif 000025
【図６】
2018112478.tif 000026
【図７】
2018112478.tif 000027
【図８】
2018112478.tif 000028
【図９】
2018112478.tif 000029
【図１０】
2018112478.tif 000030
【図１１】
2018112478.tif 000031