JP7488451B2

JP7488451B2 - 三次元計測装置

Info

Publication number: JP7488451B2
Application number: JP2020066152A
Authority: JP
Inventors: 太一丸谷
Original assignee: Denso Wave Inc
Current assignee: Denso Wave Inc
Priority date: 2020-04-01
Filing date: 2020-04-01
Publication date: 2024-05-22
Anticipated expiration: 2040-04-01
Also published as: JP2021162510A

Description

位相シフト法により物体の三次元形状を計測する三次元計測装置に関する。

特許文献１に開示されているように、計測対象物の三次元形状を計測する三次元計測装置として、位相シフト法を用いる装置が知られている。位相シフト法では位相をずらした複数枚の縞パターン画像を投影し、縞パターン画像が投影された計測対象物を撮影する。そして、計測対象物の各点の輝度をもとに、三角測量の原理で計測対象物の各点の座標を決定する。

特開２０１４－２０２４９２号公報

計測対象物が鏡面物体であると、計測対象物を撮影した画像が白飛びする可能性がある。撮影画像が白飛びしてしまうと、位相シフト法による三次元形状の測定精度が低下する。

また、計測対象物が半透明物体であると、表面で光が反射せず、内部反射が検出されることがある。内部反射光に基づいて三次元形状を計測してしまうと、三次元形状の測定精度が低下してしまう。

本開示は、この事情に基づいて成されたものであり、その目的は、精度よく計測対象物の三次元形状を計測できる三次元計測装置を提供することにある。

上記目的は独立請求項に記載の特徴の組み合わせにより達成され、また、下位請求項は更なる有利な具体例を規定する。特許請求の範囲に記載した括弧内の符号は、一つの態様として後述する実施形態に記載の具体的手段との対応関係を示すものであって、開示した技術的範囲を限定するものではない。

上記目的を達成するための三次元計測装置は、
縞パターン画像を投影するプロジェクタ（２０）と、
計測対象物（５）に投影された縞パターン画像を撮影するカメラ（４０、１４０）とを備え、
カメラが撮影した画像の輝度値に基づいて、位相シフト法により計測対象物の三次元形状を計測する三次元計測装置であって、
プロジェクタが投影した光を直線偏光にする偏光子（３０）と、
縞パターン画像が表された光が計測対象物によって反射された反射光から、少なくとも４つの振動方向の光を選択的に通過させる検光子（４１、１４１）と、
検光子を通過してカメラにより検出された光の輝度値を取得する輝度値取得部（Ｓ１２）と、
輝度値取得部が取得した輝度値に基づき、検光子が通過させる光の振動方向の違いによる輝度値の変動幅と、検光子が通過させる光の振動方向の違いによらない輝度値のベース値とを算出する輝度成分算出部（Ｓ１３）と、
変動幅とベース値とを比較して、小さい方を形状計算用の輝度値として用いて計測対象物の三次元形状を算出する形状算出部（Ｓ１４）と、を備える。

この三次元計測装置は、直線偏光を計測対象物に投影する。したがって、計測対象物からの鏡面反射成分の振動方向は特定の方向になる。また、この三次元計測装置は、反射光から、少なくとも４つの振動方向の光を選択的に通過させる検光子を備えている。検光子が、振動方向が異なる光を選択的に通過させるので、検光子を通過した光のうち、鏡面反射成分は、検光子が通過させる光の振動方向に応じて輝度が相違する。

一方、計測対象物からの反射光には拡散反射成分が含まれる。拡散反射成分には、電場の振動方向が種々の方向を持つ光が含まれている。換言すれば、拡散反射成分は、光の振動方向によらず、ほぼ一定と考えることができる。よって、反射光の拡散反射成分は、検光子が通過させる光の振動方向によらず、ほぼ同一になる。

これらのことから、４つ以上の振動方向について検光子を通過した輝度値を比較した場合の変動幅が、反射光の鏡面反射成分の大きさであると考えることができる。また、４つ以上の振動方向について検光子を通過した輝度値を比較した場合の、検光子が通過させる振動方向によらない部分であるベース値が、反射光の拡散反射成分であると考えることができる。

計測対象物が鏡面物体であれば、鏡面反射成分を使ってしまうと、その値が飽和していることがあるので、輝度値が三次元形状を反映していない。そこで、拡散反射成分の輝度値を表すベース値を使う。鏡面物体であれば、鏡面反射成分を表す輝度値の変動幅よりも、拡散反射成分を表す輝度値のベース値の方が小さい。

また、計測対象物が半透明物体である場合、計測対象物の内部で拡散反射が生じることがある。したがって、計測対象物が半透明物体である場合、輝度値の拡散反射成分は三次元形状を正確に表していない。そこで、計測対象物が半透明物体であれば、鏡面反射成分を使うことが好ましい。鏡面反射成分は変動幅であり、計測対象物が半透明物体である場合、変動幅はベース値よりも小さい。

次に、計測対象物が鏡面物体ではないが半透明物体でもない物体（以下、通常物体）である場合を考える。通常物体は、物体表面での拡散反射が反射光の主成分となる物体である。通常物体については、三次元形状を計測する輝度値に、反射光の鏡面反射成分を使っても、拡散反射成分を使ってもよい。

ただし、拡散反射成分は、表面拡散成分であるか内部拡散反射成分であるかは区別できない。鏡面反射成分を使えば、その区別を必要としないので、通常物体については、鏡面反射成分を使えばよいことになる。通常物体では、鏡面反射成分のほうが拡散反射成分よりも小さいと考えられる。

以上をまとめると、計測対象物が、鏡面物体であっても、半透明物体であっても、通常物体であっても、変動幅とベース値を比較して、小さい方を形状計算用の輝度値として用いればよいことになる。

形状算出部は、変動幅とベース値とを比較して、小さい方を形状計算用の輝度値として用いて計測対象物の三次元形状を算出する。よって、計測対象物の表面性状によらず、精度よく三次元形状を算出することができる。

また、三次元計測装置において、カメラ（４０）は、互いに通過させる光の振動方向が相違する４種類の検光子（４１）を備えた偏光カメラとすることができる。偏光カメラを備える場合、カメラよりも光路の手前に検光子を別途配置する必要がない。

また、三次元計測装置において、検光子（１４１）は、反射光の光路においてカメラ（１４０）の前に配置されていてもよい。このようにすれば、偏光カメラではないカメラを用いることができる。

実施形態の三次元計測装置１の構成を示す図。縞パターン画像を示す図。カメラ４０が備える検光子４１の繰り返し単位を示す図。検光子４１の配列の一部を示す図。カメラ４０が備えるカラーフィルタを示す図。第１実施形態における三次元形状を計測する処理を示す図。偏光角度と光検出素子が検出する光強度との関係を示す図。水平座標（ｘｍ、ｙｍ）の算出方法を説明する図である。第２実施形態における三次元形状を計測する処理を示す図。第３実施形態における三次元形状を計測する処理を示す図。第４実施形態の三次元計測装置１００の構成を示す図。

＜第１実施形態＞
以下、実施形態を図面に基づいて説明する。図１は、第１実施形態の三次元計測装置１の構成を示す図である。三次元計測装置１は、制御装置１０と、プロジェクタ２０と、偏光子３０と、カメラ４０とを備えている。三次元計測装置１は、作業台２の上に置かれた計測対象物５の三次元形状を位相シフト法により計測する。作業台２の上面は平面であり、作業台２の任意の位置に計測対象物５が位置する。三次元計測装置１は、たとえば、ロボットにピッキング、組付け作業、製品検査等を行わせる際のロボットの目として利用する。

制御装置１０は、コンピュータを備えたものとすることができる。制御装置１０は、プロジェクタ２０が投影する画像のデータとなる画像データを生成してプロジェクタ２０へ出力する。プロジェクタ２０が投影する画像には、縞パターン画像がある。

また、制御装置１０は、プロジェクタ２０から縞パターン画像が計測対象物５に投影された状態で、カメラ４０が撮影した画像を表す画像データを取得する。そして、その画像データをもとに位相シフト法により、計測対象物５の三次元形状を計測する。

図２に縞パターン画像を示す。本実施形態の縞パターン画像は、カラーの縞パターン画像である。なお、これとは異なり、単色の縞パターン画像を投影してもよい。図２に示す縞パターン画像は、詳しくは、赤、緑、青ともに、０から２５５までの輝度範囲で輝度を変化させた３色の単色縞パターン画像を合成した合成縞パターン画像である。

単色縞パターン画像は、赤、緑、青のいずれか１色の輝度が画像の一方向には正弦波状に変化し、その一方向と直交する方向は輝度が一定である画像である。合成縞パターン画像は、３色の単色縞パターン画像の位相が所定の角度だけずれている。

一例としては、赤色の単色縞パターン画像の位相が最も進んでおり、緑色の単色縞パターン画像の位相がそれよりも２π／３遅れている。青色の単色縞パターン画像は、緑色の単色縞パターン画像よりもさらに２π／３だけ位相が遅れている。縞パターン画像は、赤、緑、青の単色縞パターン画像が均等に含まれる画像であるため、ｘ画素座標の変化に伴い虹状に色が変化する。

プロジェクタ２０は、カラー画像を投影可能なプロジェクタである。プロジェクタ２０が投影する縞パターン画像は、偏光子３０を通過して計測対象物５に投影される。偏光子３０は、プロジェクタ２０が投影した光を直線偏光にする。偏光子３０による偏光方向は特に制限はない。

カメラ４０は、カラー画像を撮影可能なデジタルカメラであり、フォトダイオードなどの光検出素子（すなわち画素）を受光面に縦横に多数備えている。カメラ４０は偏光カメラになっている。偏光カメラであるカメラ４０は、撮像素子よりも光路手前に、図３に示す検光子４１を備える。図３には、検光子４１の繰り返し単位を示している。検光子４１は、詳しくは、互いに通過する光の振動方向が異なる４種類の検光子４１ａ、４１ｂ、４１ｃ、４１ｄに分かれる。これら４種類の検光子４１ａ、４１ｂ、４１ｃ、４１ｄを区別しないときは検光子４１と記載する。

本実施形態の４種類の検光子４１ａ、４１ｂ、４１ｃ、４１ｄは、通過する光の振動方向の角度差が、いずれも４５度で等しくなっている。検光子４１ａを通過する光の振動方向を０度とすると、検光子４１ｂを通過する光の振動方向は４５度、検光子４１ｃを通過する光の振動方向は９０度、検光子４１ｄを通過する光の振動方向は１３５度である。

図４には検光子４１の配列の一部を示す。図４に示すように、４種類の検光子４１を各１つずつ含む繰り返し単位が、縦横に連続して配列されている。この検光子４１と光検出素子の間に図５に示すＲＧＢカラーフィルタが配置されている。

ＲＧＢカラーフィルタは、赤と緑と青のいずれかのカラーフィルタが各光検出素子よりも光到来方向に配置されたものである。赤と緑と青のカラーフィルタの配列は、一般にベイヤ配列に従っている。図５において、Ｒは赤色フィルタ、Ｇは緑色フィルタ、Ｂは青色フィルタを意味する。

各色のフィルタの最小単位は１つの検光子４１に対応する。各色のフィルタは、その最小単位が縦横に２つずつ配列されている。したがって、各色のフィルタとも、４種類すべての検光子４１を通過した光が入射する。フィルタは、最小単位４つ分の大きさを１つの繰り返し単位としてベイヤ配列に従い配列されている。

［三次元形状を計測する処理］
次に、三次元形状を計測する処理を説明する。図６に三次元形状を計測する処理を示している。図６に示す処理は、ユーザの操作に基づき、制御装置１０が実行する。ステップ（以下、ステップを省略）Ｓ１１では、縞パターン画像を計測対象物５に投影し、カメラ４０により、そのときの計測対象物５の画像を撮影する。

Ｓ１２では、赤、緑、青の３色の色別の撮影画像を、検光子４１の種類別に取得する。これは、カラーフィルタの色別かつ検光子４１の種類別に、画素が検出した輝度値を取得することを意味する。このＳ１２は輝度値取得部に相当する。

Ｓ１３は輝度成分算出部に相当する。Ｓ１３では、色別かつ検光子４１の種類別の撮影画像を用いて偏光解析を行う。偏光解析は、図７に示すＩｍａｘ－Ｉｍｉｎと、Ｉｍｉｎを算出するものである。

図７の横軸は偏光角度である。偏光角度は４種類の検光子４１が通過させる光の振動方向を意味する。縦軸は各画素が検出する光強度、すなわち、輝度値である。したがって、Ｉｍａｘ－Ｉｍｉｎは、検光子４１の違いによる輝度値の変動幅であり、Ｉｍｉｎは検光子４１によらない輝度値のベース値である。

ここで、偏光角度を連続的に変化させて輝度値を検出したと仮定すると、図７に示すように、偏光角度に対する輝度値の変化は正弦波状になるはずである。したがって、画素が検出する輝度値をＩ、偏光角度をθ、偏光角度の初期値をφとすると、式１が成立する。

式１において、（１＋ｓｉｎ（２×（θ＋φ）））／２は、最大値が１、最小値が０であり、「θ＋φ」は図７に示す正弦波の各偏光角度における位相を意味する。変動幅であるＩｍａｘ－Ｉｍｉｎに、（１＋ｓｉｎ（２×（θ＋φ）））／２を乗じた右辺第１項は、各偏光角度における変動成分の大きさを意味する。

式１においてＩｍａｘ、Ｉｍｉｎ、θ、φが未知数であり、Ｉは各画素が検出した輝度値である。４つの偏光角度、０度、４５度、９０度、１３５度ではＩを取得できる。したがって、４つの未知数に対して４つの等式が得られるので連立方程式を解くことによりＩｍａｘ、Ｉｍｉｎ、θ、φを求めることができる。そして、Ｉｍａｘ、Ｉｍｉｎを求めたら、それらから輝度値の変動幅すなわちＩｍａｘ－Ｉｍｉｎを算出する。

ここで、変動幅とベース値が表す輝度成分の意味を説明する。図７にも記載しているように、変動幅は反射光の鏡面反射成分を表しており、ベース値は反射光の拡散反射成分を表している。また、図７から分かるように、鏡面反射成分は反射光の交流成分であり、拡散反射成分は反射光の直流成分である。

次に、変動幅が反射光の鏡面反射成分を表している理由を説明する。本実施形態の三次元計測装置１は偏光子３０を備えており、直線偏光を計測対象物５に投影する。鏡面反射する場合、反射光も直線偏光が維持される。検光子４１の偏光角度が、反射光の直線偏光の角度と一致しているほど、その検光子４１に対応する画素が検出する輝度値は大きくなる。つまり、各画素の手前に配置されている検光子４１の偏光角度が異なることにより、各画素に検出される反射光の鏡面反射成分の大きさは変化する。それ故、変動幅は、反射光の鏡面反射成分を表しているのである。

次に、ベース値が反射光の拡散反射成分を表している理由を説明する。計測対象物５に投影される光が直線偏光であっても、拡散反射成分は全方位光になる。全方位光であるので、偏光角度がどの方向であっても拡散反射成分の大きさはほぼ同じになる。したがって、ベース値が反射光の拡散反射成分を表していると考えることができるのである。

Ｓ１４は形状算出部に相当する。Ｓ１４は、具体的にはＳ１５以下である。Ｓ１５では形状計算用の輝度値を色別に決定する。形状計算用の輝度値は、Ｓ１３で算出した変動幅とベース値とを比較して小さい方の値である。

この理由を説明する。計測対象物５が鏡面物体である場合には、鏡面反射成分を使ってしまうと、その値が飽和していることがあるので、輝度値が三次元形状を表していない。一方、鏡面物体でも、通常、少しは反射光に拡散反射成分が含まれる。鏡面物体からの反射光に含まれる拡散反射成分であれば、飽和はしていない。したがって、計測対象物５が鏡面物体である場合、拡散反射成分を使うことが好ましい。鏡面物体からの反射光において、鏡面反射成分と拡散反射成分とを比較すると、拡散反射成分の方が小さい。

次に、計測対象物５が半透明物体である場合を考える。半透明物体では、計測対象物５の内部で拡散反射が生じることがある。したがって、計測対象物５が半透明物体である場合、輝度値の拡散反射成分は三次元形状を正確に表していない。そこで、計測対象物５が半透明物体であれば、鏡面反射成分を使うことが好ましい。鏡面反射は物体の表面で生じる反射だからである。半透明物体での鏡面反射成分は大きくない。したがって、半透明物体での鏡面反射成分は拡散反射成分よりも小さい。

次に、計測対象物５が通常物体である場合を考える。通常物体は、鏡面物体でも半透明物体でもない物体である。通常物体であれば、鏡面反射成分を採用してもよいし、拡散反射成分を採用してもよい。ただし、任意の物体の表面形状を計測する場合には、計測対象物５が半透明物体である可能性もある。拡散反射は、物体内部での拡散反射であるか、物体表面での拡散反射であるかを区別することが難しい。そこで、通常物体であれば、鏡面反射成分を使うことにすればよい。通常物体では、拡散反射成分よりも鏡面反射成分の方が小さい。

以上より、物体が鏡面物体でも、半透明物体でも、通常物体でも、拡散反射成分と鏡面反射成分を比較して小さい方を採用すればよいことになる。

Ｓ１６では、Ｓ１５で決定した形状計算用の輝度値を色別に正規化する。Ｓ１７では、各座標（ｘ、ｙ）の色別に正規化した輝度値をもとに、式２から、各座標（ｘ、ｙ）における位相θ（ｘ、ｙ）を算出する。なお、図５を用いて説明したように、本実施形態では、色毎に縦横に２つずつ、合計４つの画素が１つの繰り返し単位となっている。ここでの座標は、その繰り返し単位ごとの座標である。たとえば、４つの画素の中心座標を、ここでの座標とする。

式２において、Ｎは位相シフト総回数、ｎは色別に取得した撮影画像の位相シフト回数である。縞パターン画像において最も早い位相とした色のｎが０、次に位相が早い色のｎが１、最も位相が遅い色のｎが２である。位相シフト総回数Ｎは３である。また、ａ（ｘ、ｙ）は輝度振幅、ｂ（ｘ、ｙ）は背景輝度、θ（ｘ、ｙ）はｎ＝０での位相θである。

式２において、未知数は、ａ（ｘ、ｙ）、ｂ（ｘ、ｙ）、θ（ｘ、ｙ）の３つである。したがって、Ｓ１６で色別に決定した３つの撮像画像についての各座標（ｘ、ｙ）の輝度値を用いれば、位相θを含む、３つの未知数、ａ（ｘ、ｙ）、ｂ（ｘ、ｙ）、θ（ｘ、ｙ）を算出することができる。

Ｓ１８では、Ｓ１７で算出した各座標（ｘ、ｙ）の位相θ（ｘ、ｙ）から、座標計測点Ｐの高さ座標ｚｍを決定する。座標計測点Ｐは、計測対象物５あるいは作業台２の表面上の点である。

高さ座標ｚｍは、プロジェクタ２０とカメラ４０とを含む平面から物体までの距離である。高さ座標ｚｍは、位相θと高さ座標ｚｍとの関係を示すグラフと、Ｓ１７で算出した位相θとを用いて決定する。位相θと高さ座標ｚｍとの関係を示すグラフは、プロジェクタ２０の座標、カメラ４０の座標、高さ座標ｚｍ、基準面における縞パターン画像の１周期分の長さが分かれば作成することができる。なお、基準面は、作業台２の表面である投影面に平行であって、プロジェクタ２０およびカメラ４０までの距離がｚｍとなっている面である。

プロジェクタ２０とカメラ４０を固定すれば、プロジェクタ２０の座標、カメラ４０の座標は既知になる。また、基準面までの高さ座標ｚｍは与える値である。さらに、基準面までの高さ座標ｚｍが決まれば、その基準面における縞パターン画像の１周期分の長さも決まる。よって、位相θと高さ座標ｚｍとの関係を示すグラフは事前に求めることができる。

位相θと高さ座標ｚｍとの関係を示すグラフを事前に求めておき、Ｓ１８では、事前に求めた上記グラフに、Ｓ１７で算出した位相θ（ｘ、ｙ）を当てはめて、各座標計測点Ｐの高さ座標ｚｍを決定する。なお、位相θが何周期目であるかが不明だと、高さ座標ｚｍも決定することができない。しかし、ある座標計測点Ｐにおける高さ座標ｚｍは、その座標計測点Ｐに隣接する座標計測点Ｐの高さ座標ｚｍに対して連続的な変化をする。したがって、位相θが何周期目であるかが不明でも計測対象物５の三次元形状を計測することはできる。また、作業台２の高さ座標ｚｍは既知であるので、作業台２の高さ座標ｚｍと比較をすることで、座標計測点Ｐの高さ座標ｚｍを決定してもよい。

Ｓ１９では、Ｓ１８で高さ座標ｚｍを決定した座標計測点Ｐについて、水平座標（ｘｍ、ｙｍ）を決定する。Ｓ１８において決定した高さ座標ｚｍは、画素上の座標には対応付けられている。画素上の座標が決まると、カメラ４０に対する方向（αｘ、αｙ）は定まる。なお、αｘは、図８に示すように、カメラ４０から座標計測点Ｐに向かう方向のうち、ｘｍｚｍ平面におけるｚｍ軸との間の角度である。αｙは図８には図示していないが、αｙはカメラ４０から座標計測点Ｐに向かう方向のうち、ｙｍｚｍ平面におけるｚｍ軸との間の角度である。図８から分かるように、水平座標（ｘｍ、ｙｍ）は、高さ座標ｚｍとαｍ、αｙから幾何学計算により算出することができる。

Ｓ１４の処理を各座標（ｘ、ｙ）に対して実行することで、計測対象物５の三次元形状を計測することができる。

［実施形態のまとめ］
以上、説明した本実施形態では、偏光子３０により直線偏光とした縞パターン画像を計測対象物５に投影する。また、カメラ４０は、偏光角度が互いに相違する４つの検光子４１を備える。この構成により、輝度値の変動幅（すなわち反射光の鏡面反射成分の大きさ）と、輝度値のベース値（すなわち拡散反射成分の大きさ）を決定することができる（Ｓ１３）。そして、変動幅とベース値とを比較して小さい方を形状計算用の輝度値とすることで（Ｓ１５）、計測対象物５の性状によらず、計測対象物５の三次元形状を精度よく算出することができる。

また、カメラ４０として偏光カメラを用いているので、カメラ４０よりも手前に検光子を配置する必要がない。

＜第２実施形態＞
次に、第２実施形態を説明する。この第２実施形態以下の説明において、それまでに使用した符号と同一番号の符号を有する要素は、特に言及する場合を除き、それ以前の実施形態における同一符号の要素と同一である。また、構成の一部のみを説明している場合、構成の他の部分については先に説明した実施形態を適用できる。

第２実施形態は、計測対象物５が鏡面物体であることが事前に分かっている場合に適用できる実施形態である。図９に第２実施形態において図６に代えて実行する処理を示す。図９では、図６と同じＳ１１、Ｓ１２を実行して、赤、緑、青の３色の色別の撮影画像を、検光子４１の種類別に取得する。

続いて、図６のＳ１３、Ｓ１４に代えてＳ２３、Ｓ２４を実行する。Ｓ２３は、Ｓ１３と同じ計算をして、Ｉｍａｘ、Ｉｍｉｎ、θ、φを求める。Ｓ１３との違いは、Ｉｍａｘ－Ｉｍｉｎは算出しない点である。

Ｓ２４は具体的には、Ｓ２６、Ｓ２７、Ｓ２８、Ｓ２９の処理である。これらＳ２６、Ｓ２７、Ｓ２８、Ｓ２９は、それぞれ、Ｓ１６、Ｓ１７、Ｓ１８、Ｓ１９と同じあるいは類似の処理である。Ｓ１４にはＳ１５が含まれていたが、Ｓ２４にはＳ１５に相当する処理はない。計測対象物５が鏡面物体である場合には、形状計算用の輝度値にはベース値であるＩｍｉｎを使えばよいことが分かっている。そのため、輝度値の変動幅の大きさとベース値の大きさを比較する必要がないからである。

Ｓ２６では、第２実施形態における形状計算用の輝度値である各座標（ｘ、ｙ）のＩｍｉｎを色別に正規化する。Ｓ２７、Ｓ２８、Ｓ２９は、それぞれ、Ｓ１７、Ｓ１８、Ｓ１９と同じである。したがって、Ｓ２７では、各座標（ｘ、ｙ）の色別に正規化した輝度値をもとに、式２から、各座標（ｘ、ｙ）における位相θ（ｘ、ｙ）を算出する。Ｓ２８では、Ｓ２７で算出した各座標（ｘ、ｙ）の位相θ（ｘ、ｙ）から、座標計測点Ｐの高さ座標ｚｍを決定する。Ｓ２９では、Ｓ２８で高さ座標ｚｍを決定した座標計測点Ｐについて、水平座標（ｘｍ、ｙｍ）を、高さ座標ｚｍとαｍ、αｙから幾何学計算により算出する。

計測対象物５が鏡面物体であることが分かっている場合には、この第２実施形態のように、輝度値のベース値と変動幅の比較を省略し、輝度値のベース値であるＩｍｉｎを形状計算用の輝度値として用いて計測対象物５の三次元形状を計測することができる。

＜第３実施形態＞
第３実施形態は、計測対象物５が半透明物体であることが事前に分かっている場合に適用できる実施形態である。図１０に第３実施形態において図６に代えて実行する処理を示す。図１０では、図６と同じＳ１１、Ｓ１２を実行して、赤、緑、青の３色の色別の撮影画像を、検光子４１の種類別に取得する。

続くＳ３３は図６のＳ１３と同じであり、Ｉｍａｘ、Ｉｍｉｎ、θ、φを求め、さらに、Ｉｍａｘ、ＩｍｉｎからＩｍａｘ－Ｉｍｉｎを算出する。

Ｓ３４は具体的には、Ｓ３６、Ｓ３７、Ｓ３８、Ｓ３９の処理である。これらＳ３６、Ｓ３７、Ｓ３８、Ｓ３９は、それぞれ、Ｓ１６、Ｓ１７、Ｓ１８、Ｓ１９と同じあるいは類似の処理である。Ｓ１４にはＳ１５が含まれていたが、Ｓ３４にはＳ１５に相当する処理はない。計測対象物５が半透明物体である場合には、形状計算用の輝度値には変動幅であるＩｍａｘ－Ｉｍｉｎを使えばよいことが分かっている。そのため、輝度値の変動幅の大きさとベース値の大きさを比較する必要がないからである。

Ｓ３６では、第３実施形態における形状計算用の輝度値である各座標（ｘ、ｙ）のＩｍａｘ－Ｉｍｉｎを色別に正規化する。Ｓ３７、Ｓ３８、Ｓ３９は、それぞれ、Ｓ１７、Ｓ１８、Ｓ１９と同じである。したがって、Ｓ３７では、各座標（ｘ、ｙ）における位相θ（ｘ、ｙ）を算出する。Ｓ３８では、座標計測点Ｐの高さ座標ｚｍを決定する。Ｓ３９では、座標計測点Ｐの水平座標（ｘｍ、ｙｍ）を算出する。

計測対象物５が半透明物体であることが分かっている場合、この第３実施形態のように、輝度値のベース値と変動幅の比較を省略し、輝度値の変動幅であるＩｍａｘ－Ｉｍｉｎを形状計算用の輝度値として用いて計測対象物５の三次元形状を計測することができる。

＜第４実施形態＞
図１１に第４実施形態の三次元計測装置１００の構成を示す。三次元計測装置１００は、偏光カメラではないカメラ１４０を備える。そして、このカメラ１４０よりも、反射光の光路においてカメラ１４０の前となる位置に、検光子１４１が１つ配置されている。検光子１４１は、たとえば円形薄板状であり、厚さ方向に貫通する中心軸周りに回転可能である。

検光子１４１は、その中心軸周りの任意の回転角度で固定させることができる。したがって、ある回転角度を０度とすると、０度、４５度、９０度、１３５度の回転角度で固定することができる。これら４つの回転角度で固定した状態は、それぞれ、互いに異なる４つの振動方向の光を選択的に通過させる状態である。

三次元計測装置１００を用いて計測対象物５の三次元形状を計測する場合、Ｓ１１において、検光子１４１の回転角度を４種類以上の回転角度として、それぞれの回転角度で、縞パターン画像の投影と、計測対象物５の撮影を行えばよい。

この三次元計測装置１００は、偏光カメラではないカメラ１４０を用いることができる。また、全部の画素が同じ偏光角度の輝度を検出する。したがって、三次元形状の分解能を向上することができる。

以上、実施形態を説明したが、開示した技術は上述の実施形態に限定されるものではなく、次の変形例も開示した範囲に含まれ、さらに、下記以外にも要旨を逸脱しない範囲内で種々変更して実施できる。

＜変形例１＞
実施形態では、カラーの縞パターン画像を投影していた。しかし、互いに位相が異なる３つの単色の縞パターン画像を投影してもよい。

＜変形例２＞
第４実施形態の検光子１４１は、回転可能に構成されていた。しかし、回転可能に構成することに代えて、所定の取付部材に対して検光子１４１が着脱可能になっており、偏光方向を異ならせた検光子１４１を、順次、取付部材に装着してもよい。

１：三次元計測装置２：作業台５：計測対象物１０：制御装置２０：プロジェクタ３０：偏光子４０：カメラ４１：検光子１００：三次元計測装置１４０：カメラ１４１：検光子Ｓ１２：輝度値取得部Ｓ１３：輝度成分算出部Ｓ１４：形状算出部Ｓ２３：輝度成分算出部Ｓ２４：形状算出部Ｓ３３：輝度成分算出部Ｓ３４：形状算出部

Claims

縞パターン画像を投影するプロジェクタ（２０）と、
計測対象物（５）に投影された前記縞パターン画像を撮影するカメラ（４０、１４０）とを備え、
前記カメラが撮影した画像の輝度値に基づいて、位相シフト法により前記計測対象物の三次元形状を計測する三次元計測装置であって、
前記プロジェクタが投影した光を直線偏光にする偏光子（３０）と、
前記縞パターン画像が表された光が前記計測対象物によって反射された反射光から、少なくとも４つの振動方向の光を選択的に通過させる検光子（４１、１４１）と、
前記検光子を通過して前記カメラにより検出された光の輝度値を取得する輝度値取得部（Ｓ１２）と、
前記輝度値取得部が取得した輝度値に基づき、前記検光子が通過させる光の振動方向の違いによる前記輝度値の変動幅と、前記検光子が通過させる光の振動方向の違いによらない前記輝度値のベース値とを算出する輝度成分算出部（Ｓ１３）と、
前記変動幅と前記ベース値とを比較して、小さい方を形状計算用の輝度値として用いて前記計測対象物の三次元形状を算出する形状算出部（Ｓ１４）と、を備える三次元計測装置。
請求項１に記載の三次元計測装置であって、
前記カメラ（４０）は、互いに通過させる光の振動方向が相違する４種類の前記検光子（４１）を備えた偏光カメラである、三次元計測装置。
請求項１に記載の三次元計測装置であって、
前記検光子（１４１）は、前記反射光の光路において前記カメラ（１４０）の前に配置されている三次元計測装置。