Deprecated: The each() function is deprecated. This message will be suppressed on further calls in /home/zhenxiangba/zhenxiangba.com/public_html/phproxy-improved-master/index.php on line 456
JP4046282B2 - Range finder - Google Patents
[go: Go Back, main page]

JP4046282B2 - Range finder - Google Patents

Range finder Download PDF

Info

Publication number
JP4046282B2
JP4046282B2 JP2003110466A JP2003110466A JP4046282B2 JP 4046282 B2 JP4046282 B2 JP 4046282B2 JP 2003110466 A JP2003110466 A JP 2003110466A JP 2003110466 A JP2003110466 A JP 2003110466A JP 4046282 B2 JP4046282 B2 JP 4046282B2
Authority
JP
Japan
Prior art keywords
camera
pixel
incident
light beam
image
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Expired - Fee Related
Application number
JP2003110466A
Other languages
Japanese (ja)
Other versions
JP2004317246A (en
Inventor
千秋 青山
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Honda Motor Co Ltd
Original Assignee
Honda Motor Co Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Honda Motor Co Ltd filed Critical Honda Motor Co Ltd
Priority to JP2003110466A priority Critical patent/JP4046282B2/en
Priority to US10/806,547 priority patent/US7627197B2/en
Publication of JP2004317246A publication Critical patent/JP2004317246A/en
Application granted granted Critical
Publication of JP4046282B2 publication Critical patent/JP4046282B2/en
Anticipated expiration legal-status Critical
Expired - Fee Related legal-status Critical Current

Links

Images

Landscapes

  • Length Measuring Devices By Optical Means (AREA)
  • Measurement Of Optical Distance (AREA)

Description

【0001】
【発明の属する技術分野】
本発明は、非接触で対象物の形状を測定するレンジファインダに関する。
【0002】
【従来の技術】
近年、CCDカメラやコンピュータ画像の処理の発展に伴って画像を使用した3次元計測が一般的になってきている。CCDカメラとコンピュータ画像処理を用いた3次元計測装置として、レンジファインダがある。レンジファインダは、対象物に対し、投光器からスリット光などの平行光線を照射し、対象物をカメラで撮像し、その撮像した画像から対象物の形状を測定する装置である。その原理は、投光器とカメラの位置関係と、撮像した画像中の平行光線の位置とから、三角測量に基づいて対象物の表面位置を計算する(例えば、特許文献1参照)。
【0003】
【特許文献1】
特開平9−287927号公報(段落0025、図2等)
【0004】
【発明が解決しようとする課題】
ところで、従来のレンジファインダでは、投光器が、固定された一点から平行光線を照射し、カメラに入射される光線は、一点に収束することを前提として形状を測定している。つまり、カメラのレンズ系がピンホールカメラモデルに基づいていることが前提となっている。このピンホールカメラモデルとは、図7に示すように、基点位置(針穴:ピンホールH)を通して入射する光(入射光)のみが、画像面上に到達して、三次元空間が、画像面上の二次元空間に対応付けられるモデルを言う。このようにピンホールカメラモデルは、入射光線が一点のピンホールを通って画像が形成されることを想定している。
【0005】
しかし、ピンホールカメラではなく、ガラス等からなるレンズ系を有するカメラでは入射光線を延長しても1点に収束しない。そのため、レンズ系を有するカメラで撮像した画像には、非線形の歪みが存在し、周辺視野ほどその歪みが大きい。
したがって、従来のレンジファインダでは、この歪みの分だけ対象物の位置測定結果に誤差を生じていた。この誤差は、対象物との距離にも依存し、かつ対象物の距離による誤差の変化が非線形であった。そのため、対象物の距離がわかっていれば、位置測定結果と、実際の位置との関係をとって較正することは可能ではあるが、対象物の距離は測定するまで未知であるので、完全に較正することはできなかった。
【0006】
本発明は、以上のような問題点に鑑みてなされたものであり、従来のレンジファインダにおいて、ピンホールカメラモデルを前提とした誤差の根本原因を除去し、カメラと対象物の距離によらず正確な形状測定をすることが可能なレンジファインダを提供することを目的とする。
【0007】
【課題を解決するための手段】
本発明は、前記目的を達成するために創案されたものであり、まず、請求項1に記載のレンジファインダは、投光器により対象物に平行光線を照射し、カメラにより前記平行光線が照射された対象物を撮像し、前記投光器の発光位置と前記カメラの光学中心との位置関係、及び撮像された画像における前記平行光線の照射位置から、前記対象物の形状を測定するレンジファインダであって、前記カメラで撮像した画像を入力する画像入力手段と、この画像入力手段で入力された画像において、照射された平行光線の画素位置を検出する画素位置検出手段と、前記カメラで撮像される画像の画素位置と、その画素への入射光線の方向及び基準位置から入射光線への変位量とを関連づけた較正テーブルを蓄積した較正テーブル蓄積手段と、前記画素位置検出手段で検出された画素位置に基づき、前記較正テーブルを参照して、その画素位置に対応する入射光線の方向及び変位量を取得し、この入射光線の方向及び変位量を利用して前記照射位置を計算する位置算出手段とを備えることを特徴とする。
【0008】
このようなレンジファインダは、まず、対象物に対し、投光器により平行光線が照射される。この平行光線は、投光器から広がることなく照射される光であれば良く、1点で照射されるスポット光や、帯状の光であるスリット光や、特定のパターンを有するパターン光などを単独で、またはこれらを組み合わせて使用することができる。
そして、画像入力手段が、カメラで撮像した画像をレンジファインダ内に入力し、画素位置検出手段がこの画像から照射された平行光線の画素位置を認識する。さらに位置算出手段が、較正テーブルを参照して前記画素位置に対応する入射光線の方向及び変位量を取得し、この入射光線の方向及び基準位置からの変位量を利用して照射位置、すなわち対象物の表面位置を計算する。
そのため、レンズに入射してくる入射光線の正確な方向に基づき、対象物の形状を測定できるので、対象物の形状を正確に測定することができる。なお、学中心(opticalcenter)とは、レンズ系の中央位置を示し、ピンホールカメラモデルのレンズ系においては、ピンホールの位置に該当するものである。
【0009】
また、請求項2に記載のレンジファインダは、請求項1に記載のレンジファインダにおいて、前記基準位置が、前記カメラと一定関係にある固定点を表し、前記変位量が、前記基準位置から入射光線上の任意の一点へ向かうベクトルで表され、前記任意の一点が、前記カメラから、視線方向のベクトルを延ばす基点となり、前記視線方向が、前記画素への入射光線の方向を示し、前記較正テーブルが、前記基準位置に対して、前記画素位置ごとに前記変位量が関連づけられていることを特徴とする。
【0011】
【発明の実施の形態】
次に、図面を参照しながら本発明の実施形態について説明する。
まず、一般にピンホールカメラモデルと呼ばれるカメラにおいて、画像の歪みの原因となる、入射光線が一点で交わらないカメラの非ピンホール性について説明し、その非ピンホール性を有するカメラの特性を数値化したキャリブレーションデータについて説明する。そして、カメラの撮像画素毎のキャリブレーションデータを測定して較正テーブルを生成する方法について説明し、この較正テーブルを有するレンジファインダについて順次説明する。
【0012】
[カメラの非ピンホール性について]
まず、図8を参照して、レンズ系を有するカメラで撮像した画像において、歪みが発生する原因について説明する。図8は、レンズ系を有するカメラの模式図である。ここでは、説明を簡略化するため、レンズ系を板ガラスGとし、ンホールHが生成されているものとする。このカメラCの板ガラスGに垂直に入射する入射光線r1は、ピンホールHを通って撮像面I上の画素R1に撮像される。また、板ガラスGに斜めに入射した入射光線r2及びr3は、板ガラスG内で屈折した後にピンホールHを通って撮像面I上の画素R2及びR3に撮像される。
【0013】
しかし、このカメラCは、板ガラスGを通過する前の入射光線r2及びr3の延長線であるr2´及びr3´と、入射光線r1とは、一点では交わらず、実際にはピンホールカメラモデルとはなっていないことがわかる。このため、撮像面I上の画素R3には、ピンホールカメラモデルで想定している入射光線rrとは、距離D分だけずれた入射光線r3が撮像されることになる。
【0014】
このように、レンズ系(ここでは板ガラスG)に入射される入射光線によって像を撮像するカメラは、ピンホール性が崩れていることになる(非ピンホール性)。以下、レンズ系を有するカメラを「非ピンホールカメラ」と呼ぶこととする。
【0015】
[キャリブレーションデータについて]
次に、図9を参照して、非ピンホールカメラの特性を数値化したキャリブレーションデータについて説明する。図9は、キャリブレーションデータの内容を説明するための説明図である。図9に示すように、レンズlに入射する入射光線Rは2点で特定することができる。ここでは、第1の光源位置P1と、第2の光源位置P2とから発光される光が同一の画素(図示せず)に撮像されたときに、入射光線Rがその画素に対応する入射光線であると特定する。
【0016】
ここで、すべての入射光線との距離の自乗和が最小となる点を光学中心Oと定義し、各撮像画素に対応する入射光線Rと光学中心Oとの距離が最小となる点を、その入射光線Rの入射光基点Kと定義する。
【0017】
すなわち、光学中心O(x0,y0,z0)は、すべての入射光線において、光源位置P1(x1,y1,z1)と光源位置P2(x2,y2,z2)とで特定される、入射光線Rからの距離dの自乗((1)式)和が最小になる位置を、最小自乗法によって求めた位置となる。
【0018】
d2=-(A2/B)+C ・・・(1)
【0019】
ただし、
A=(x2-x1)(x1-x0)+(y2-y1)(y1-y0)+(z2-z1)(z1-z0)
B=(x2-x1)2+(y2-y1)2+(z2-z1)2
C=(x1-x0)2+(y1-y0)2+(z1-z0)2
とする。
【0020】
これによって、入射光線の方向(光源位置P1及びP2で特定される方向)と、光学中心Oから入射光基点Kへの変位量(3次元ベクトルVD(dx,dy,dz)で表現)とを対応付けたデータをキャリブレーションデータとすることで、非ピンホールカメラの特性を数値化することができる。
【0021】
なお、キャリブレーションデータは、これに限定されるものではない。例えば、前記した例では、光学中心Oを基準位置とし、光学中心Oから入射光線に降ろした垂線の足までのベクトルを変位量VDとしているが、基準位置は、光学中心に限らず、カメラと一定関係にある固定点であれば、どのような点でも構わない。そして、変位量VDは、基準位置から入射光線上の任意の一点へ向かうベクトルであればよく、基準位置から入射光線へ降ろした垂線の足へ向かうベクトルには限られない。
【0022】
[較正テーブルの生成方法]
次に、図10を参照して、非ピンホールカメラの特性を数値化したキャリブレーションデータをカメラの撮影画素(投影画素に相当)毎に関連付けた較正テーブルを生成する方法について説明する。図10は、較正テーブルを生成する方法の原理を示す概念図である。なお、図10(a)は、特定の入射光線に対して、カメラのパン及びチルトを変化させることで、キャリブレーションデータを測定する原理を示す概念図であり、図10(b)は、固定したカメラに対して、入射光線を変化させることで、キャリブレーションデータを測定する原理を示す概念図である。
【0023】
図10(a)に示すように、キャリブレーションデータを撮影画素毎に関連付けた較正テーブルを生成するには、非ピンホール性を有するカメラCに対して、光源位置をP1及びP2の1方向に移動させ(1軸移動)、光源位置P1及びP2で特定される入射光線Rを決定し、光源位置P1及びP2から発光される光が共に、測定を行う撮像画素に入射されるように、カメラCのパン及びチルトを調整(2軸回転)することで、カメラCの撮像画素毎に入射される入射光線Rの方向を特定する。
【0024】
また、図10(b)に示すように、カメラCは固定しておき、光源位置P1及びP2の2点で発光した入射光線Rが、測定画素に入射されるように、光源位置P1及びP2をXYZ方向に移動(3軸移動)させることで、その測定画素に入射される光源位置P1及びP2の2点で定まる入射光線Rの方向を特定することとしてもよい。
【0025】
図10(a)又は(b)で、測定を行った各撮像画素毎に特定された入射光線Rに基づいて、入射光線Rの方向と、光学中心Oから入射光基点K(図9参照)への変位量とをキャリブレーションデータとして撮像画素(投影画素)毎に関連付けることで較正テーブルを生成することができる。
【0026】
[レンジファインダの構成]
次に、図1を参照して、レンジファインダについて説明を行う。図1は、本発明の実施の形態であるレンジファインダの構成を示したブロック図である。図1に示したレンジファインダ1は、レーザ照射装置50で対象物OBに対し照射した光(平行光線)を、カメラCで撮像し、得られた画像から、対象物OBの3次元位置を検出するものである。そして、レーザ照射装置50の光束の方向を上下、左右に振って、多数の点で3次元位置を測定することにより、対象物OBの形状を測定する。
ここでは、レンジファインダ1を、レーザ照射装置(投光器)50と、カメラCと、画像入力手段10と、画素位置検出手段20と、較正テーブル蓄積手段30と、位置算出手段40とを備えて構成した。
【0027】
レーザ照射装置50は、図2に示すブロック図のように、レーザ光を発光するレーザ発光部51と、レーザ発光部51から照射されたレーザ光を集光して細いビームにする集光レンズ52と、集光レンズ52で集光されたレーザ光を複数のビームに分ける回折格子53とを備えている。回折格子53は、紙面に垂直な方向へビームを分けるものである。
また、回折格子の後段には、分かれたビームを1方向に拡散してスリット光を生成するビーム拡散レンズ54が設けられている。ビーム拡散レンズ54は、シリンドリカルレンズ等で構成される。このビーム拡散レンズ54によって、複数のビームのそれぞれは、60°の角度で拡散される。
図3に、レーザ光の照射状態を図示する。図3に示すように、レーザ照射装置50から発せられたレーザ光LBは、円錐面状に拡散され、平面PL上で二次曲線の軌跡RLを形成し、反射・散乱している。この反射・散乱が、特許請求の範囲にいう照射位置を示すことになる。レーザ光LBは、回折格子53によって、5つのビーム(レーザ光LB)に分けられており、さらにビーム拡散レンズ54によって、5つのビームが60°に拡散されている。これらレーザ光LBは、5つに分かれたビームごとに異なる原点OLから放射状に広がるように発せられている。
なお、図3では、見やすくするためにビームの数を「5」としているが、実用上は、より細かい間隔でレーザ光を分ける。たとえば、図3に示す角度Bが32°であり、隣り合うレーザ光の角度Cを1.6°とする。すなわち、ビームの数は「21」となる。
【0028】
また、レーザ照射装置50は、レンズ系で回折・拡散するものに限られず、一点のスポット光を発するレーザ光を、上方向(x方向)、横方向(y方向)に機械的に回動させてもよい。その場合には、レーザ光が、一点から放射状に発せられるように回動機構を構成すると、対象物の位置算出が容易になる。
【0029】
画像入力手段10は、カメラCで対象物OBを撮像した画像を入力するものである。入力する画像は、カラー画像や、多階調の白黒画像として入力する。画像入力手段10には、カメラCで撮像された各画像を一時的に記憶しておく、図示していないメモリを備えており、このメモリに記憶された画像を、後記する画素位置検出手段20が参照するものとする。
【0030】
画素位置検出手段20は、画像入力手段10で入力された画像において、レーザ光の照射位置を検出し、その画素位置を特定するものである。具体的には、画像上の明度分布を見て、明度がピークとなる画素を検出する。ここで検出された各画像の画素位置は、位置算出手段40へ入力される。
【0031】
較正テーブル蓄積手段30は、ハードディスク等の一般的な記憶媒体であって、カメラCの撮像画素毎にキャリブレーションデータを対応付けた較正テーブル31を蓄積したものである。較正テーブル31は、例えば、図4のように画素のx座標、y座標の組合せに対し、その画素に入射してくる光線を特定する情報として、光学中心Oからの変位量VD(dx,dy,dz)及び角度α、βを対応づけて記憶されている。図5に示すように、角度α、βは、レーザ照射装置50のレーザ光LBが照射される原点OLを座標軸の原点とし、レーザ照射装置50の向きを基準にX軸、Y軸、Z軸を設定し、図5のXZ平面に対する入射光線の角度をα、図5のXY平面に対する入射光線の角度をβとしている。
【0032】
位置算出手段40は、画素位置検出手段20で検出された各画像の画素位置に対応した、較正テーブル30のキャリブレーションデータ(入射光線の方向及び変位量VD)に基づいて、対照物OBの位置(3次元位置)を算出するものである。
【0033】
[レンジファインダの動作]
次に、図5及び図6を参照して、レンジファインダ1の動作について説明する。図6は、レンジファインダ1の動作を示すフローチャートである。
まず、レンジファインダ1は、対象物OBにレーザ光LBを照射し(ステップS1)、カメラCで対象物OBを撮像する(ステップS2)。そして、画像入力手段10によって、撮像した画像を入力する(ステップS3)。
次に、レンジファインダ1は、画像入力手段10で入力された各画像において、画素位置検出手段20により、明度のピークを検出して、レーザ光LBの対象物OB上の照射位置に対応する画素位置を検出する(ステップS4)。
【0034】
そして、レンジファインダ1は、位置算出手段40によって、較正テーブル31を参照して、前記画素位置に対応するキャリブレーションデータを取得し(ステップS5)、そのキャリブレーションデータに基づいて、カメラ基点位置PC(Cx,Cy,Cz)を算出する(ステップS6)。カメラ基点位置は、図5に示すように、カメラCから、視線方向のベクトルを延ばす基点となる位置である。図5の座標系で、カメラCの光学中心Oの座標を、O(Ox,Oy,Oz)とすると、カメラ基点位置PC(Cx,Cy,Cz)は、
Cx=Ox+dx
Cy=Oy+dy
Cz=Oz+dz
により表される。なお、本実施の形態では、較正テーブル31に、画素位置と、dx,dy,dzを対応させ、Cx,Cy,Czを別途計算により算出しているが、較正テーブル31に、画素位置とCx,Cy,Czを対応させて記憶させても構わない。その場合には、特許請求の範囲にいう基準位置は、図5の座標系の原点、すなわち、レーザ光LBが照射される原点OLとなる。
【0035】
そして、カメラ基点位置PC(Cx,Cy,Cz)と、較正テーブル31から取得した角度α、β、及びレーザ光の方向を特定する角度ω、角度φに基づき、照射位置を計算する。なお、角度ωは、予め測定する対象までの距離が限定されている場合には、画像上の一定の範囲内には所定のレーザ光しか存在しないため、画素位置からレーザ光を決定でき対応する角度も決定できることになる。また、スポット光を機械的に振るようにレーザ照射装置を構成した場合には、そのスポット光が向けられた方向を機械的又は電気的に測定し、角度ω、角度φを画像と同期して位置算出手段40に取り込み、計算に使用すればよい。
【0036】
位置の計算は、例えば次のようにして行う。
図5のように、レーザ光BLの原点OLから照射位置P0(X,Y,Z)までの距離をr、カメラCのカメラ基点位置PCから照射位置P0(X,Y,Z)までの距離をRとする。
そこで、照射位置P0(X,Y,Z)は、レーザ光BLの方向ベクトルおよび距離rとの関係では、幾何学的に次の(2)〜(4)式
X=rcosωcos φ ・・・(2)
Y=rsinω ・・・(3)
Z=‐rcosωsinφ・・・(4)
の関係があり、視線方向ベクトルおよび距離Rとの関係では、次の(5)〜(7)式
X=Rcosαcosβ + Cx・・・(5)
Y=Rsinαcosβ + Cy・・・(6)
Z=-Rsinβ + Cz ・・・(7)
の関係がある。
これらの式(2)〜(7)から、r、R、を求めると、
(0≦Cy and 0≦ω) or (Cy<0 and ω<0)のとき
R=(-B + √(B2-4AC))/2A
(0≦Cy and ω<0) or (Cy<0 and 0≦ω)のとき
R=(-B - √(B2-4AC))/2A
ω0のとき
r=(Rsinαcosβ + Cy)/sinω
ただし、
A=cos2αcos2β+sin2β ‐ cos2ωsin2αcos2β/sin2ω
B=cosαcosβ Cx - sinβ Cz ‐ cos2ωsinαcosβ Cy/sin2ω
C=Cx2 + Cz2 ‐ cos2ω Cy2/sin2ω
ω= 0 のとき
R=-Cy/sin α cos β
となる。
【0037】
このようにして、レンジファインダ1は、カメラCの非ピンホール性を補正して、正確に照射位置を検出することが可能になる。そして、すべてのレーザ光LBの照射位置について3次元位置を計算すれば、対象物OBの3次元形状を測定することができる。
なお、このレンジファインダ1を、移動ロボット、自動車等に組み込んで用いることも可能である。例えば、移動ロボットに本発明を適用し、移動ロボットが、床の形状を検出することで、床の凹凸を正確に認識することができ、移動ロボットが安定した歩行を行うことが可能になる。
【0038】
【発明の効果】
以上説明したとおり、本発明に係るレンジファインダによれば、カメラの非ピンホール性を考慮して、光学中心と実際にカメラのレンズ系に入射する光のズレを補正し、対象物の正確な位置・形状を測定することができる。
【図面の簡単な説明】
【図1】実施形態に係るレンジファインダの構成を示すブロック図である。
【図2】レーザ照射装置のブロック図である。
【図3】レーザの照射状態を示す斜視図である。
【図4】較正テーブルの一例を示す図である。
【図5】照射位置の計算方法を説明する図である。
【図6】実施形態に係るレンジファインダの動作を示すフローチャートである。
【図7】ピンホールカメラモデルの概念を説明する図である。
【図8】レンズ系を有するカメラの模式図である。
【図9】キャリブレーションデータの内容を説明する図である。
【図10】キャリブレーションデータの生成方法を説明する図である。
【符号の説明】
1 レンジファインダ
10 画像入力手段
20 画素位置検出手段
30 較正テーブル蓄積手段
31 較正テーブル
40 位置算出手段
50 レーザ照射装置
C カメラ
[0001]
BACKGROUND OF THE INVENTION
The present invention relates to a range finder that measures the shape of an object in a non-contact manner.
[0002]
[Prior art]
In recent years, with the development of CCD camera and computer image processing, three-dimensional measurement using images has become common. There is a range finder as a three-dimensional measuring apparatus using a CCD camera and computer image processing. The range finder is a device that irradiates a target with parallel light beams such as slit light from a projector, images the target with a camera, and measures the shape of the target from the captured image. The principle is that the surface position of the object is calculated based on triangulation from the positional relationship between the projector and the camera and the position of the parallel rays in the captured image (see, for example, Patent Document 1).
[0003]
[Patent Document 1]
Japanese Patent Laid-Open No. 9-287927 (paragraph 0025, FIG. 2, etc.)
[0004]
[Problems to be solved by the invention]
By the way, in the conventional range finder, the projector measures the shape on the premise that the light beam irradiates a parallel light beam from a fixed point and the light beam incident on the camera converges to the single point. That is, it is assumed that the camera lens system is based on a pinhole camera model. As shown in FIG. 7, this pinhole camera model is such that only light (incident light) incident through the base point position (needle hole: pinhole H) reaches the image plane, and the three-dimensional space is an image. A model associated with a two-dimensional space on a surface. Thus, the pinhole camera model assumes that an incident light beam passes through a single pinhole to form an image.
[0005]
However, not a pinhole camera but a camera having a lens system made of glass or the like does not converge to one point even if the incident light beam is extended. Therefore, non-linear distortion exists in an image captured by a camera having a lens system, and the distortion is larger in the peripheral visual field.
Therefore, in the conventional range finder, an error is generated in the position measurement result of the object by this distortion. This error also depends on the distance to the object, and the change in the error due to the distance of the object is non-linear. Therefore, if the distance of the object is known, it is possible to calibrate based on the relationship between the position measurement result and the actual position, but since the distance of the object is unknown until measurement, It could not be calibrated.
[0006]
The present invention has been made in view of the above problems, and in the conventional range finder, eliminates the root cause of errors based on the pinhole camera model, regardless of the distance between the camera and the object. An object of the present invention is to provide a range finder capable of performing accurate shape measurement.
[0007]
[Means for Solving the Problems]
The present invention was devised to achieve the above object, and first, the range finder according to claim 1 irradiates an object with parallel rays by a projector and the parallel rays with a camera . A range finder that images a target object and measures the shape of the target object from the positional relationship between the light emission position of the projector and the optical center of the camera, and the irradiation position of the parallel rays in the captured image, An image input unit that inputs an image captured by the camera, a pixel position detection unit that detects a pixel position of an irradiated parallel ray in the image input by the image input unit, and an image captured by the camera Calibration table storage means for storing a calibration table associating a pixel position with the direction of the incident light beam to the pixel and the amount of displacement from the reference position to the incident light beam; Based on the pixel position detected by the position detection means, referring to the calibration table, the direction and displacement amount of the incident light beam corresponding to the pixel position are obtained, and the direction and displacement amount of the incident light beam are used to And a position calculating means for calculating the irradiation position .
[0008]
In such a range finder , first, a parallel light beam is irradiated onto an object by a projector. This parallel light may be light that is irradiated from the projector without spreading, and spot light that is irradiated at one point, slit light that is band-shaped light, pattern light having a specific pattern, and the like alone. Or they can be used in combination.
Then, the image input means inputs the image captured by the camera into the range finder, and the pixel position detection means recognizes the pixel position of the parallel light beam emitted from this image. Further, the position calculating means refers to the calibration table to acquire the direction and displacement amount of the incident light beam corresponding to the pixel position, and uses the direction of the incident light beam and the displacement amount from the reference position to obtain the irradiation position, that is, the target. Calculate the surface position of the object.
Therefore, since the shape of the object can be measured based on the exact direction of the incident light incident on the lens, the shape of the object can be accurately measured. Incidentally, it with the optical center (opticalcenter) indicates the center position of the lens system, the lens system of the pinhole camera model, it corresponds to the position of the pinhole.
[0009]
The range finder according to claim 2 is the range finder according to claim 1, wherein the reference position represents a fixed point having a fixed relationship with the camera, and the amount of displacement is incident light from the reference position. The calibration table is represented by a vector directed to an arbitrary point on the line, the arbitrary point becomes a base point extending a vector of the line-of-sight direction from the camera, and the line-of-sight direction indicates the direction of incident light on the pixel. However, the displacement amount is associated with the reference position for each pixel position .
[0011]
DETAILED DESCRIPTION OF THE INVENTION
Next, embodiments of the present invention will be described with reference to the drawings.
First, we will explain the non-pinhole characteristics of a camera that does not intersect incident light at a single point, which causes image distortion in a camera that is generally called a pinhole camera model, and quantifies the characteristics of the non-pinhole camera. The calibration data thus made will be described. Then, a method for generating calibration tables by measuring calibration data for each imaging pixel of the camera will be described, and a range finder having this calibration table will be sequentially described.
[0012]
[Non-pinhole characteristics of the camera]
First, with reference to FIG. 8, the cause of distortion in an image captured by a camera having a lens system will be described. FIG. 8 is a schematic diagram of a camera having a lens system. Here, in order to simplify the description, the lens system and the glass sheet G, it is assumed that the pinhole H is generated. Incident light beam r1 incident perpendicularly on the glass plate G of the camera C passes through the pinhole H and is imaged on the pixel R1 on the imaging surface I. In addition, incident rays r2 and r3 incident obliquely on the plate glass G are refracted in the plate glass G, pass through the pinhole H, and are imaged on the pixels R2 and R3 on the imaging surface I.
[0013]
However, in this camera C, the incident rays r2 ′ and r3 ′, which are the extension lines of the incident rays r2 and r3 before passing through the glass sheet G, and the incident ray r1 do not intersect at one point, but actually are a pinhole camera model. You can see that it is not. For this reason, the incident light beam r3 that is shifted by the distance D from the incident light beam rr assumed in the pinhole camera model is imaged on the pixel R3 on the imaging surface I.
[0014]
Thus, a camera that captures an image with incident light rays incident on the lens system (here, the glass sheet G) has a pinhole property (non-pinhole property). Hereinafter, a camera having a lens system is referred to as a “non-pinhole camera”.
[0015]
[About calibration data]
Next, calibration data obtained by quantifying the characteristics of the non-pinhole camera will be described with reference to FIG. FIG. 9 is an explanatory diagram for explaining the contents of the calibration data. As shown in FIG. 9, the incident light ray R incident on the lens 1 can be specified by two points. Here, when light emitted from the first light source position P1 and the second light source position P2 is captured by the same pixel (not shown), the incident light R corresponds to the pixel. To be identified.
[0016]
Here, the point at which the sum of squares of the distances to all incident rays is minimized is defined as the optical center O, and the point at which the distance between the incident ray R and the optical center O corresponding to each imaging pixel is minimized. This is defined as an incident light base point K of the incident light R.
[0017]
That is, the optical center O (x 0 , y 0 , z 0 ) has a light source position P1 (x 1 , y 1 , z 1 ) and a light source position P2 (x 2 , y 2 , z 2 ) for all incident rays. The position where the sum of the squares of the distance d from the incident light ray R (equation (1)) is minimized is determined by the least square method.
[0018]
d 2 =-(A 2 / B) + C (1)
[0019]
However,
A = (x 2 -x 1 ) (x 1 -x 0 ) + (y 2 -y 1 ) (y 1 -y 0 ) + (z 2 -z 1 ) (z 1 -z 0 )
B = (x 2 -x 1 ) 2 + (y 2 -y 1 ) 2 + (z 2 -z 1 ) 2
C = (x 1 -x 0 ) 2 + (y 1 -y 0 ) 2 + (z 1 -z 0 ) 2
And
[0020]
As a result, the direction of the incident light (the direction specified by the light source positions P1 and P2), the amount of displacement from the optical center O to the incident light base point K (represented by a three-dimensional vector V D (dx, dy, dz)), and By using the data associated with the as the calibration data, the characteristics of the non-pinhole camera can be digitized.
[0021]
Note that the calibration data is not limited to this. For example, in the above-described example, the optical center O is set as the reference position, and the vector from the optical center O to the foot of the perpendicular dropped to the incident ray is set as the displacement amount V D. However, the reference position is not limited to the optical center, Any point can be used as long as it is a fixed point that is in a fixed relationship. The displacement amount V D may be a vector that goes from the reference position to an arbitrary point on the incident light beam, and is not limited to a vector that goes from the reference position to the perpendicular foot dropped to the incident light beam.
[0022]
[Calibration table generation method]
Next, a method for generating a calibration table in which calibration data obtained by quantifying the characteristics of a non-pinhole camera is associated with each shooting pixel (corresponding to a projection pixel) of the camera will be described with reference to FIG. FIG. 10 is a conceptual diagram showing the principle of a method for generating a calibration table. FIG. 10A is a conceptual diagram showing the principle of measuring calibration data by changing the pan and tilt of the camera with respect to a specific incident ray, and FIG. 10B is a fixed diagram. It is a conceptual diagram which shows the principle which measures calibration data by changing incident light rays with respect to the camera which did.
[0023]
As shown in FIG. 10A, in order to generate a calibration table in which calibration data is associated with each imaging pixel, the light source position is set in one direction of P1 and P2 with respect to the camera C having non-pinhole characteristics. The camera is moved (one-axis movement) to determine the incident light ray R specified by the light source positions P1 and P2, and the light emitted from the light source positions P1 and P2 is incident on the imaging pixel to be measured. By adjusting the pan and tilt of C (biaxial rotation), the direction of the incident light beam R incident on each imaging pixel of the camera C is specified.
[0024]
Further, as shown in FIG. 10B, the camera C is fixed, and the light source positions P1 and P2 are set so that incident light rays R emitted at two points of the light source positions P1 and P2 are incident on the measurement pixel. Is moved in the XYZ directions (triaxial movement), and the direction of the incident light ray R determined by two points of the light source positions P1 and P2 incident on the measurement pixel may be specified.
[0025]
In FIG. 10 (a) or (b), based on the incident light ray R specified for each imaging pixel that has been measured, the direction of the incident light ray R and the incident light base point K from the optical center O (see FIG. 9). The calibration table can be generated by associating the amount of displacement to the image data (projection pixel) with each other as calibration data.
[0026]
[Range finder configuration]
Next, the range finder will be described with reference to FIG. FIG. 1 is a block diagram showing a configuration of a range finder according to an embodiment of the present invention. The range finder 1 shown in FIG. 1 picks up the light (parallel rays) irradiated to the object OB by the laser irradiation device 50 with the camera C, and detects the three-dimensional position of the object OB from the obtained image. To do. Then, the shape of the object OB is measured by measuring the three-dimensional position at a large number of points by swinging the direction of the light beam of the laser irradiation device 50 up and down and left and right.
Here, the range finder 1 includes a laser irradiation device (projector) 50, a camera C, an image input unit 10, a pixel position detection unit 20, a calibration table storage unit 30, and a position calculation unit 40. did.
[0027]
As shown in the block diagram of FIG. 2, the laser irradiation device 50 includes a laser light emitting unit 51 that emits laser light, and a condensing lens 52 that condenses the laser light emitted from the laser light emitting unit 51 into a thin beam. And a diffraction grating 53 for dividing the laser beam condensed by the condenser lens 52 into a plurality of beams. The diffraction grating 53 divides the beam in a direction perpendicular to the paper surface.
In addition, a beam diffusing lens 54 for diffusing the divided beams in one direction to generate slit light is provided at the subsequent stage of the diffraction grating. The beam diffusion lens 54 is configured by a cylindrical lens or the like. By the beam diffusion lens 54, each of the plurality of beams is diffused at an angle of 60 °.
FIG. 3 illustrates a laser light irradiation state. As shown in FIG. 3, the laser beam LB emitted from the laser irradiation device 50 is diffused in a conical surface, forms a quadratic locus RL on the plane PL, and is reflected and scattered. This reflection / scattering indicates the irradiation position in the claims. The laser beam LB is divided into five beams (laser beam LB) by the diffraction grating 53, and the five beams are diffused by 60 ° by the beam diffusion lens 54. These laser beam LB is emitted as radiating from different origin O L in each beam which is divided into five.
In FIG. 3, the number of beams is “5” for easy viewing, but in practice, the laser light is divided at finer intervals. For example, the angle B shown in FIG. 3 is 32 °, and the angle C between adjacent laser beams is 1.6 °. That is, the number of beams is “21”.
[0028]
Further, the laser irradiation device 50 is not limited to the one that diffracts and diffuses in the lens system, but mechanically rotates the laser light that emits a single spot light in the upward direction (x direction) and the lateral direction (y direction). May be. In that case, if the rotation mechanism is configured such that the laser light is emitted radially from one point, the position of the object can be easily calculated.
[0029]
The image input means 10 inputs an image obtained by capturing the object OB with the camera C. The input image is input as a color image or a multi-tone monochrome image. The image input unit 10 includes a memory (not shown) that temporarily stores each image captured by the camera C. The image position stored in the memory is a pixel position detection unit 20 to be described later. Shall refer to.
[0030]
The pixel position detection unit 20 detects the irradiation position of the laser beam in the image input by the image input unit 10 and specifies the pixel position. Specifically, by looking at the brightness distribution on the image, the pixel having the peak brightness is detected. The pixel position of each image detected here is input to the position calculation means 40.
[0031]
The calibration table storage means 30 is a general storage medium such as a hard disk, and stores a calibration table 31 in which calibration data is associated with each imaging pixel of the camera C. For example, as shown in FIG. 4, the calibration table 31 has a displacement amount V D (dx, dx, d) from the optical center O as information for specifying a light ray incident on the pixel with respect to a combination of the x coordinate and the y coordinate of the pixel. dy, dz) and angles α and β are stored in association with each other. As shown in FIG. 5, the angle alpha, beta is the origin O L in which the laser beam LB of the laser irradiation apparatus 50 is irradiated as the origin of the coordinate axes, X-axis relative to the orientation of the laser irradiation device 50, Y-axis, Z An axis is set, and the angle of incident light with respect to the XZ plane in FIG. 5 is α, and the angle of incident light with respect to the XY plane in FIG. 5 is β.
[0032]
The position calculation means 40 is based on the calibration data (the direction of the incident light and the displacement amount V D ) of the calibration table 30 corresponding to the pixel position of each image detected by the pixel position detection means 20. The position (three-dimensional position) is calculated.
[0033]
[Rangefinder operation]
Next, the operation of the range finder 1 will be described with reference to FIGS. FIG. 6 is a flowchart showing the operation of the range finder 1.
First, the range finder 1 irradiates the object OB with the laser beam LB (step S1), and images the object OB with the camera C (step S2). Then, the captured image is input by the image input means 10 (step S3).
Next, in the range finder 1, the pixel corresponding to the irradiation position on the object OB of the laser beam LB is detected by the pixel position detecting unit 20 in each image input by the image input unit 10. The position is detected (step S4).
[0034]
Then, the range finder 1 acquires calibration data corresponding to the pixel position by the position calculation means 40 with reference to the calibration table 31 (step S5), and based on the calibration data, the camera base position P C (Cx, Cy, Cz) is calculated (step S6). As shown in FIG. 5, the camera base point position is a position that becomes a base point for extending a vector in the viewing direction from the camera C. In the coordinate system of FIG. 5, when the coordinates of the optical center O of the camera C are O (Ox, Oy, Oz), the camera base position P C (Cx, Cy, Cz) is
Cx = Ox + dx
Cy = Oy + dy
Cz = Oz + dz
Is represented by In the present embodiment, the pixel position and dx, dy, dz are associated with the calibration table 31 and Cx, Cy, Cz are calculated separately, but the pixel position and Cx are calculated in the calibration table 31. , Cy, Cz may be stored in association with each other. In that case, the reference position is referred to in the claims, the coordinate system of the origin of the 5, i.e., the origin O L in which the laser beam LB is irradiated.
[0035]
Then, the irradiation position is calculated based on the camera base position P C (Cx, Cy, Cz), the angles α and β acquired from the calibration table 31, and the angle ω and the angle φ that specify the direction of the laser beam. Note that when the distance to the object to be measured is limited in advance, the angle ω corresponds to the laser beam that can be determined from the pixel position because there is only a predetermined laser beam within a certain range on the image. The angle can also be determined. In addition, when the laser irradiation apparatus is configured to mechanically shake the spot light, the direction in which the spot light is directed is measured mechanically or electrically, and the angles ω and φ are synchronized with the image. What is necessary is just to take in the position calculation means 40 and to use for calculation.
[0036]
The position is calculated as follows, for example.
As shown in FIG. 5, the distance from the origin OL of the laser beam BL to the irradiation position P 0 (X, Y, Z) is r, and the camera base point PC of the camera C to the irradiation position P 0 (X, Y, Z). The distance to is R.
Therefore, the irradiation position P 0 (X, Y, Z) is geometrically expressed by the following equations (2) to (4) in relation to the direction vector of the laser beam BL and the distance r.
X = rcosω cos φ・ ・ ・ (2)
Y = rsinω (3)
Z = ‐rcosωsin φ・ ・ ・ (4)
In relation to the line-of-sight direction vector and the distance R, the following equations (5) to (7)
X = Rcosαcosβ + Cx (5)
Y = Rsinαcosβ + Cy (6)
Z = -Rsinβ + Cz (7)
There is a relationship.
From these equations (2) to (7), r and R are obtained.
(0 ≦ Cy and 0 ≦ ω) or (Cy <0 and ω <0)
R = (-B + √ (B 2 -4AC)) / 2A
(0 ≦ Cy and ω <0) or (Cy <0 and 0 ≦ ω)
R = (-B-√ (B 2 -4AC)) / 2A
When ω 0
r = (Rsinαcosβ + Cy) / sinω
However,
A = cos 2 αcos 2 β + sin 2 β‐cos 2 ωsin 2 αcos 2 β / sin 2 ω
B = cosαcosβ Cx-sinβ Cz ‐ cos 2 ωsinαcosβ Cy / sin 2 ω
C = Cx 2 + Cz 2 ‐ cos 2 ω Cy 2 / sin 2 ω
When ω = 0
R = -Cy / sin α cos β
It becomes.
[0037]
In this way, the range finder 1 can correct the non-pinhole property of the camera C and accurately detect the irradiation position. And if the three-dimensional position is calculated about the irradiation position of all the laser beams LB, the three-dimensional shape of the target object OB can be measured.
The range finder 1 can be used by being incorporated in a mobile robot, an automobile, or the like. For example, when the present invention is applied to a mobile robot and the mobile robot detects the shape of the floor, the unevenness of the floor can be recognized accurately, and the mobile robot can stably walk.
[0038]
【The invention's effect】
As described above, according to the range finder according to the present invention, in consideration of the non-pinhole property of the camera, the deviation of the optical center and the light actually incident on the lens system of the camera is corrected, and the object is accurately detected. The position and shape can be measured.
[Brief description of the drawings]
FIG. 1 is a block diagram illustrating a configuration of a range finder according to an embodiment.
FIG. 2 is a block diagram of a laser irradiation apparatus.
FIG. 3 is a perspective view showing a laser irradiation state.
FIG. 4 is a diagram illustrating an example of a calibration table.
FIG. 5 is a diagram illustrating a method for calculating an irradiation position.
FIG. 6 is a flowchart showing an operation of the range finder according to the embodiment.
FIG. 7 is a diagram illustrating the concept of a pinhole camera model.
FIG. 8 is a schematic diagram of a camera having a lens system.
FIG. 9 is a diagram illustrating the contents of calibration data.
FIG. 10 is a diagram illustrating a calibration data generation method.
[Explanation of symbols]
DESCRIPTION OF SYMBOLS 1 Range finder 10 Image input means 20 Pixel position detection means 30 Calibration table storage means 31 Calibration table 40 Position calculation means 50 Laser irradiation apparatus C Camera

Claims (2)

投光器により対象物に平行光線を照射し、カメラにより前記平行光線が照射された対象物を撮像し、前記投光器の発光位置と前記カメラの光学中心との位置関係、及び撮像された画像における前記平行光線の照射位置から、前記対象物の形状を測定するレンジファインダであって、
前記カメラで撮像した画像を入力する画像入力手段と、
この画像入力手段で入力された画像において、照射された平行光線の画素位置を検出する画素位置検出手段と、
前記カメラで撮像される画像の画素位置と、その画素への入射光線の方向及び基準位置から入射光線への変位量とを関連づけた較正テーブルを蓄積した較正テーブル蓄積手段と、
前記画素位置検出手段で検出された画素位置に基づき、前記較正テーブルを参照して、その画素位置に対応する入射光線の方向及び変位量を取得し、この入射光線の方向及び変位量を利用して前記照射位置を算出する位置算出手段とを備えることを特徴とするレンジファインダ。
The object is irradiated with a parallel light beam by the projector, and the target object irradiated with the parallel light is imaged by the camera. The positional relationship between the light emission position of the projector and the optical center of the camera, and the parallel in the captured image. A range finder for measuring the shape of the object from a light irradiation position,
Image input means for inputting an image captured by the camera;
In the image input by the image input means, pixel position detection means for detecting the pixel position of the irradiated parallel light beam,
Calibration table accumulating means for accumulating a calibration table in which the pixel position of the image captured by the camera, the direction of the incident light beam to the pixel and the amount of displacement from the reference position to the incident light beam are associated;
Based on the pixel position detected by the pixel position detecting means, the calibration table is referred to obtain the direction and displacement amount of the incident light beam corresponding to the pixel position, and use the direction and displacement amount of the incident light beam. A position finder for calculating the irradiation position.
前記基準位置は、前記カメラと一定関係にある固定点を表し、The reference position represents a fixed point that is in a fixed relationship with the camera,
前記変位量は、前記基準位置から入射光線上の任意の一点へ向かうベクトルで表され、The amount of displacement is represented by a vector from the reference position toward an arbitrary point on the incident ray,
前記任意の一点は、前記カメラから、視線方向のベクトルを延ばす基点となり、The one arbitrary point is a base point for extending a vector in the sight line direction from the camera
前記視線方向は、前記画素への入射光線の方向を示し、The line-of-sight direction indicates the direction of light rays incident on the pixel,
前記較正テーブルは、前記基準位置に対して、前記画素位置ごとに前記変位量が関連づけられていることを特徴とする請求項1に記載のレンジファインダ。The range finder according to claim 1, wherein the displacement amount is associated with the reference position for each pixel position in the calibration table.
JP2003110466A 2003-04-07 2003-04-15 Range finder Expired - Fee Related JP4046282B2 (en)

Priority Applications (2)

Application Number Priority Date Filing Date Title
JP2003110466A JP4046282B2 (en) 2003-04-15 2003-04-15 Range finder
US10/806,547 US7627197B2 (en) 2003-04-07 2004-03-23 Position measurement method, an apparatus, a computer program and a method for generating calibration information

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2003110466A JP4046282B2 (en) 2003-04-15 2003-04-15 Range finder

Publications (2)

Publication Number Publication Date
JP2004317246A JP2004317246A (en) 2004-11-11
JP4046282B2 true JP4046282B2 (en) 2008-02-13

Family

ID=33471316

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2003110466A Expired - Fee Related JP4046282B2 (en) 2003-04-07 2003-04-15 Range finder

Country Status (1)

Country Link
JP (1) JP4046282B2 (en)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP5136108B2 (en) * 2008-02-18 2013-02-06 トヨタ自動車株式会社 3D shape measuring method and 3D shape measuring apparatus
CN108362215B (en) * 2017-10-09 2021-07-16 同济大学 Automatic deformation monitoring system and method for various measurement robots to freely set up stations

Also Published As

Publication number Publication date
JP2004317246A (en) 2004-11-11

Similar Documents

Publication Publication Date Title
JP5140761B2 (en) Method for calibrating a measurement system, computer program, electronic control unit, and measurement system
US7627197B2 (en) Position measurement method, an apparatus, a computer program and a method for generating calibration information
US7075661B2 (en) Apparatus and method for obtaining three-dimensional positional data from a two-dimensional captured image
US6741279B1 (en) System and method for capturing document orientation information with a digital camera
US9188430B2 (en) Compensation of a structured light scanner that is tracked in six degrees-of-freedom
JP5943547B2 (en) Apparatus and method for non-contact measurement
DK2993450T3 (en) Method and arrangement for recording acoustic and optical information as well as a corresponding computer program and a corresponding computer-readable storage medium
US20130194569A1 (en) Substrate inspection method
WO2004044522A1 (en) Three-dimensional shape measuring method and its device
JP2008241643A (en) Three-dimensional shape measuring device
CN101387596A (en) Method and device for horizontal adjustment of optical mechanism
JP4969279B2 (en) Position detection method and position detection apparatus
CN103676487A (en) Workpiece height measuring device and correcting method thereof
JP2012026816A (en) Dimension measuring method and device
JP4046282B2 (en) Range finder
JP4837538B2 (en) End position measuring method and dimension measuring method
JP2017198470A (en) Measurement device, measurement method, system, and goods manufacturing method
JP2000205821A (en) Instrument and method for three-dimensional shape measurement
JP2007093412A (en) Three-dimensional shape measuring device
US12439024B2 (en) Provision of real world and image sensor correspondence points for use in calibration of an imaging system for three dimensional imaging based on light triangulation
US20040263862A1 (en) Detecting peripheral points of reflected radiation beam spots for topographically mapping a surface
JP2018044863A (en) Measuring device, measuring method, system, and article manufacturing method
CN110702378B (en) Optical axis pointing detection device and detection method
CN109073371B (en) Apparatus and method for tilt detection
JP4741943B2 (en) Inspection apparatus and inspection method

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20051130

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20070815

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20070822

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20071003

TRDD Decision of grant or rejection written
A01 Written decision to grant a patent or to grant a registration (utility model)

Free format text: JAPANESE INTERMEDIATE CODE: A01

Effective date: 20071114

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20071116

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20101130

Year of fee payment: 3

R150 Certificate of patent or registration of utility model

Ref document number: 4046282

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150

Free format text: JAPANESE INTERMEDIATE CODE: R150

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20101130

Year of fee payment: 3

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20111130

Year of fee payment: 4

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20111130

Year of fee payment: 4

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20121130

Year of fee payment: 5

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20131130

Year of fee payment: 6

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

LAPS Cancellation because of no payment of annual fees