JP2008014882A - Three-dimensional measuring device - Google Patents
Three-dimensional measuring device Download PDFInfo
- Publication number
- JP2008014882A JP2008014882A JP2006188705A JP2006188705A JP2008014882A JP 2008014882 A JP2008014882 A JP 2008014882A JP 2006188705 A JP2006188705 A JP 2006188705A JP 2006188705 A JP2006188705 A JP 2006188705A JP 2008014882 A JP2008014882 A JP 2008014882A
- Authority
- JP
- Japan
- Prior art keywords
- equation
- sensor
- light
- coordinates
- dimensional
- 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.)
- Pending
Links
Images
Landscapes
- Length Measuring Devices By Optical Means (AREA)
- Length Measuring Devices With Unspecified Measuring Means (AREA)
Abstract
Description
本発明は、三次元計測装置に関し、さらに詳しくは、大型構造物の形状と、その構造物に存在する亀裂の三次元計測をリアルタイムに行なう三次元計測装置に関するものである。 The present invention relates to a three-dimensional measurement apparatus, and more particularly to a three-dimensional measurement apparatus that performs three-dimensional measurement of a shape of a large structure and a crack existing in the structure in real time.
従来からトンネルや大型下水管等の保守点検作業では、定期的にそれらの構造物を目視により観察して亀裂の有無や状態の現状把握を行なっていた。特に、トンネルなどの亀裂の状態を記録する場合、従来は構造物壁面に直接チョークでマス目を描き、目視により読み取った形状を方眼紙に記入する方法が採られていた。しかし、この方法では、作業員の個人差や測定誤差が大きく、形状も三次元であるために電子化が困難であるといった問題がある。 Conventionally, in maintenance inspection work such as tunnels and large sewage pipes, these structures are regularly observed visually to determine the presence and state of cracks. In particular, when recording the state of a crack in a tunnel or the like, conventionally, a method has been employed in which squares are drawn directly on the wall surface of the structure with chalk and the shape read visually is written on a graph paper. However, this method has a problem that it is difficult to digitize because there are large individual differences among workers and measurement errors, and the shape is also three-dimensional.
また、現在、商品化されている三次元画像計測装置では、比較的小さな構造物の計測を対象としている装置が多く、大型構造物の形状を計測するためには、多くの時間と労力を要した。即ち、大型構造物を計測する場合は、CCDカメラを移動させながら計測し、計測されたデータを繋ぐ作業が必要であった。 In addition, currently commercially available 3D image measuring devices are often used for measuring relatively small structures, and it takes a lot of time and labor to measure the shape of large structures. did. That is, when measuring a large structure, it is necessary to perform measurement while moving the CCD camera and to connect the measured data.
また、従来の三次元形状計測装置として特許文献1には、被測定物に対し、上方からこの被測定物を含む計測範囲内を縦横走査しながらレーザ光を照射し、このレーザ光及び上記被測定物上で反射されるレーザ光を外部に送出し、上記計測範囲内に複数の測定点を設定し、かつ、送出された上記被測定物に照射されるレーザ光と上記被測定物上で反射されるレーザ光との位相差を検出して上記複数の測定点各々に対する奥行き情報を求め、上記被測定物の平面形状を撮影するカメラと、このカメラによって撮影された画像情報から上記被測定物の二次元情報を求め、上記二次元情報と形状計算器から出力される上記奥行き情報とを組み合わせて上記被測定物の三次元形状情報を求める三次元形状計測装置について開示されている。
しかしながら、特許文献1に開示されている従来技術は、二次元情報と形状計算器から出力される奥行き情報とを組み合わせて被測定物の三次元形状情報を求めるため、二次元情報と奥行き情報を合成する処理が必要となり、多くの処理時間を要するばかりでなく、画像情報の画質が二次元情報に影響を与えるため、結果的に三次元形状の精度にも影響が及ぼされるといった問題がある。 However, the conventional technique disclosed in Patent Document 1 uses the two-dimensional information and the depth information to obtain the three-dimensional shape information of the object to be measured by combining the two-dimensional information and the depth information output from the shape calculator. There is a problem that not only the processing to be combined is required and a lot of processing time is required, but also the image quality of the image information affects the two-dimensional information, and consequently the accuracy of the three-dimensional shape is also affected.
本発明は、かかる課題に鑑み、撮像カメラとレーザ投光器に三次元位置と方向を検出できる三次元磁気センサを取り付け、撮像カメラとレーザ投光器をそれぞれ独立して動作させながら被測定対象物全体を計測すると共に、レーザ投光器をスポット光に切り替えて、被測定対象物に存在する亀裂を計測ことにより、形状と亀裂の両方をリアルタイムに計測し、且つ同一画面上に合成して表示することが可能な三次元計測装置を提供することを目的とする。 In view of such problems, the present invention attaches a three-dimensional magnetic sensor capable of detecting a three-dimensional position and direction to an imaging camera and a laser projector, and measures the entire object to be measured while operating the imaging camera and the laser projector independently. At the same time, by switching the laser projector to spot light and measuring the cracks present in the object to be measured, both the shape and cracks can be measured in real time and combined and displayed on the same screen. An object is to provide a three-dimensional measuring apparatus.
本発明はかかる課題を解決するために、請求項1は、物体の三次元情報に基づいて前記物体の三次元形状を計測する三次元計測装置において、前記物体に光を照射する投光手段と、該投光手段により前記物体表面に照射された光を撮像する撮像手段と、前記投光手段に備えられ該投光手段の三次元位置情報及び姿勢情報を生成する第1のセンサと、前記撮像手段に備えられ該撮像手段の三次元位置情報及び姿勢情報を生成する第2のセンサと、所定のエリアに磁界ベクトルを形成するトランスミッタと、前記物体に照射された光の平面方程式を算出する方程式算出手段と、前記投光手段のビーム形状を切り替えるビーム形状切替手段と、を備え、前記投光手段と撮像手段がそれぞれ独立に動作することにより前記第1のセンサ及び第2のセンサより得られる三次元位置情報及び姿勢情報に基づいて、前記方程式算出手段が前記物体に照射された光の平面方程式を算出する場合と、前記物体全体の平面方程式を算出する場合と、該物体に存在する亀裂の平面方程式を算出する場合とに応じて、前記ビーム形状切替手段により前記ビーム形状を切り替えるように構成したことを特徴とする。 In order to solve this problem, the present invention provides a three-dimensional measuring apparatus that measures a three-dimensional shape of the object based on three-dimensional information of the object, and a light projecting unit that irradiates the object with light. Imaging means for imaging light irradiated on the object surface by the light projecting means; a first sensor provided in the light projecting means for generating three-dimensional position information and posture information of the light projecting means; A second sensor that is provided in the imaging unit and generates three-dimensional position information and posture information of the imaging unit, a transmitter that forms a magnetic field vector in a predetermined area, and a plane equation of light irradiated on the object are calculated. An equation calculating means and a beam shape switching means for switching the beam shape of the light projecting means, and the light projecting means and the imaging means operate independently of each other, whereby the first sensor and the second sensor. Based on the obtained three-dimensional position information and posture information, the equation calculation means calculates a plane equation of light applied to the object, calculates a plane equation of the entire object, It is characterized in that the beam shape is switched by the beam shape switching means in accordance with the case of calculating the plane equation of an existing crack.
本発明では、投光手段と撮像手段を独立に動作するようにして、それぞれに三次元位置情報及び姿勢情報を生成するセンサを備え、それぞれが独立に動作することにより得られた三次元位置情報及び姿勢情報を基に平面方程式を算出するものである。また、測定対象物に存在する亀裂を計測する場合は、ビーム形状を切り替えるようにして測定するものである。これにより、大型の測定対象物の形状を正確に測定できるばかりでなく、細かい亀裂を正確に測定することができる。 In the present invention, the light projecting means and the imaging means are operated independently, and each is provided with a sensor for generating three-dimensional position information and posture information, and the three-dimensional position information obtained by operating each independently. The plane equation is calculated based on the posture information. Moreover, when measuring the crack which exists in a measuring object, it measures by switching a beam shape. Thereby, not only the shape of a large measurement object can be measured accurately, but also a fine crack can be measured accurately.
請求項2は、前記第1のセンサ及び第2のセンサは、前記トランスミッタにより形成された磁界ベクトルを受信することにより、前記第1のセンサ及び第2のセンサの三次元位置情報及びロール角、ピッチ角、及びヨー角の姿勢情報を取得することを特徴とする。
投光手段と撮像手段の位置情報と姿勢情報を得るには、各種の方法が考えられるが、光の光路を邪魔せず、且つ視界を妨げない方法が必要条件である。その点では磁界による方法が最適である。本発明では半球状の磁界を発生するトランスミッタを備え、その磁界の磁力線に対するベクトルを計算することにより三次元位置情報及び姿勢情報を取得するものである。従って、物体と投光手段及び撮像手段の距離はこの磁界のエリア内に限定される。
According to a second aspect of the present invention, the first sensor and the second sensor receive a magnetic field vector formed by the transmitter, thereby obtaining three-dimensional position information and a roll angle of the first sensor and the second sensor. Pitch angle and yaw angle attitude information is acquired.
Various methods are conceivable for obtaining the position information and the posture information of the light projecting means and the imaging means, but a method that does not disturb the optical path of light and does not disturb the field of view is a necessary condition. In that respect, the method using a magnetic field is optimal. The present invention includes a transmitter that generates a hemispherical magnetic field, and obtains three-dimensional position information and posture information by calculating a vector with respect to the magnetic field lines of the magnetic field. Accordingly, the distance between the object, the light projecting means, and the imaging means is limited to the area of the magnetic field.
請求項3は、前記方程式算出手段は、前記第1及び第2のセンサから得られた位置を(xow,yow,zow)、姿勢を(φ,θ,ψ)とした場合、前記各センサのセンサ座標に原点を置いたセンサ座標上の点(xr,yr,zr)は、出力された前記各センサのデータを基に、
ここで
C:Cos、S:Sinとすると、前記センサ座標上で光平面上の任意の3点の座標(xr,yr,zr)を前記式(1)及び(2)により前記トランスミッタを原点としたワールド座標に変換し、前記ワールド座標(xrw,yrw,zrw)に変換することにより、
前記式(3)及び式(4)として光の平面方程式を取得し、前記撮像手段の焦点から受像面までの距離をfと置くと、前記受像面上の座標(u,v)は前記撮像手段の焦点を原点とした座標系において、
となり、Z=λとおき線形化すると、
と表わされ、該式(6)を前記センサ座標で表すために、回転・平行移動を考慮して、
式(7)として表し、(ここでk11〜k33のパラメータには前記撮像手段の位置や姿勢などをはじめとする計測対象と、前記撮像手段の位置を表すデータが全て含まれているものとする)該式(7)が前記センサ座標(x,y,z)と前記受像面上の座標(u,v)の関係式となり、前記式(7)を展開し、整理すると式(8)を得て、
前記式(8)の2平面の交線で表される直線は前記撮像手段の焦点から計測点に向かう直線として
式(9)のレシーバ座標を求め、前記式(1)と同様にワールド座標に変換して前記撮像手段の焦点から計測点に向かう直線である式(10)を求め、
(ここで、(wx,wy,wz)はトランスミッタを原点としたレシーバのワールド座標である)
前記式(10)とレーザ平面上の方程式である前記式(4)を連立させてワールド座標(X,Y,Z)を求めることを特徴とする。
According to a third aspect of the present invention, in the equation calculation means, when the position obtained from the first and second sensors is (x ow , y ow , z ow ) and the posture is (φ, θ, ψ), The point (x r , y r , z r ) on the sensor coordinates with the origin at the sensor coordinates of each sensor is based on the output data of each sensor.
here
Assuming that C: Cos and S: Sin, the coordinates (x r , y r , z r ) of any three points on the optical plane on the sensor coordinates are the origins of the transmitter according to the equations (1) and (2). To the world coordinates (x rw , y rw , z rw ),
When the plane equation of light is acquired as the equations (3) and (4), and the distance from the focal point of the imaging means to the image receiving surface is set to f, the coordinates (u, v) on the image receiving surface are the imaging In the coordinate system with the focal point of the means as the origin,
And when linearizing with Z = λ,
In order to express the equation (6) in the sensor coordinates, in consideration of rotation / translation,
Expressed as Equation (7), (wherein k 11 and the measurement object, including position and orientation of the imaging means to the parameters to k 33, those data representing the position of the imaging means contains all The expression (7) becomes a relational expression between the sensor coordinates (x, y, z) and the coordinates (u, v) on the image receiving surface. When the expression (7) is developed and rearranged, the expression (8) )
The straight line represented by the intersection of the two planes in the equation (8) is a straight line from the focal point of the imaging means to the measurement point.
Obtain the receiver coordinates of Equation (9), convert to world coordinates as in Equation (1), and obtain Equation (10), which is a straight line from the focal point of the imaging means to the measurement point,
(Where (wx, wy, wz) are the world coordinates of the receiver with the transmitter as the origin)
World coordinates (X, Y, Z) are obtained by combining the above expression (10) and the above expression (4) which is an equation on the laser plane.
請求項4は、前記方程式算出手段は、前記平面方程式をリアルタイムに算出することにより前記投光手段の独立走査を可能としたことを特徴とする。
物体に照射された光の平面方程式は方程式算出手段により計算される。このとき計算の処理速度が遅い場合、一旦投光手段により照射された平面情報を記憶してその後バッチ処理により方程式の計算処理を行う。このように計算処理がリアルタイムに行われないと、投光手段を物体の任意の場所に操作することが不可能となり、撮像手段と独立にした意味がなくなってしまう。そこで本発明では、平面方程式をリアルタイムに計算できるようにするために、第1及び第2のセンサから得られた位置情報と姿勢情報を独立に取得することにより、演算速度高めることを可能とした。
According to a fourth aspect of the present invention, the equation calculating means enables independent scanning of the light projecting means by calculating the plane equation in real time.
The plane equation of the light irradiated to the object is calculated by the equation calculating means. If the calculation processing speed is low at this time, the plane information once irradiated by the light projecting means is temporarily stored, and then the equation calculation processing is performed by batch processing. Thus, if the calculation process is not performed in real time, it becomes impossible to operate the light projecting means at an arbitrary place on the object, and the meaning independent of the imaging means is lost. Therefore, in the present invention, in order to be able to calculate the plane equation in real time, it is possible to increase the calculation speed by independently acquiring the position information and the posture information obtained from the first and second sensors. .
請求項5は、前記ビーム形状切替手段は、前記物体全体の平面方程式を算出する場合は、該ビーム形状をスリット光に切り替え、前記物体に存在する亀裂の平面方程式を算出する場合は、前記ビーム形状をスポット光に切り替えることを特徴とする。
物体全面の平面方程式を算出する場合は、投光手段のビーム形状をスリット状にして走査することにより、大型構造物の表面の情報を正確に且つ迅速に得ることができる。一方、亀裂のような細い複雑な形状を測定するには、投光手段のビーム形状をスポット状にしてその亀裂に沿って走査することが好ましい。
The beam shape switching means may switch the beam shape to slit light when calculating a plane equation of the entire object, and calculate the plane equation of a crack existing in the object. The shape is switched to spot light.
When calculating the plane equation of the entire surface of the object, information on the surface of the large structure can be obtained accurately and quickly by scanning with the beam shape of the light projecting means being slit. On the other hand, in order to measure a thin and complicated shape such as a crack, it is preferable to make the beam shape of the light projecting means a spot shape and scan along the crack.
請求項6は、前記投光手段の光源はレーザ光若しくはLEDであることを特徴とする。
光は光源から照射されてから物体に到達するまでできるだけ拡散しないことが好ましい。その点でレーザ光は優れている。しかし、光源と物体との距離が至近の場合は、光の拡散は少ないのでLEDの光でも可能である。
According to a sixth aspect of the present invention, the light source of the light projecting means is a laser beam or an LED.
It is preferable that light is not diffused as much as possible until it reaches the object after being irradiated from the light source. In this respect, the laser beam is excellent. However, when the distance between the light source and the object is very close, the light diffusion is small, so LED light is also possible.
本発明によれば、物体に光を照射する投光手段と、この投光手段により物体表面に照射された光を撮像する撮像手段と、投光手段に備えられこの投光手段の三次元位置情報及び姿勢情報を生成する第1のセンサと、撮像手段に備えられこの撮像手段の三次元位置情報及び姿勢情報を生成する第2のセンサと、所定のエリアに磁界ベクトルを形成するトランスミッタと、物体に照射された光の平面方程式を算出する方程式算出手段と、投光手段のビーム形状を切り替えるビーム形状切替手段と、を備え、投光手段と撮像手段がそれぞれ独立に動作することにより第1のセンサ及び第2のセンサより得られる三次元位置情報及び姿勢情報に基づいて、方程式算出手段が物体に照射された光の平面方程式を算出する場合、物体全体の平面方程式を算出する場合と、物体に存在する亀裂の平面方程式を算出する場合とで、ビーム形状切替手段によりビーム形状を切り替えるようにしたので、測定対象物に存在する亀裂を精度良く計測することができ、且つ測定対象物内に亀裂の画像を自動的に合成することができる。 According to the present invention, the light projecting means for irradiating the object with light, the image capturing means for capturing the light irradiated on the object surface by the light projecting means, and the three-dimensional position of the light projecting means provided in the light projecting means A first sensor that generates information and orientation information; a second sensor that is provided in the imaging means and generates three-dimensional position information and orientation information of the imaging means; a transmitter that forms a magnetic field vector in a predetermined area; Equation calculating means for calculating a plane equation of light irradiated to the object, and beam shape switching means for switching the beam shape of the light projecting means, and the light projecting means and the imaging means operate independently, respectively. When the equation calculation means calculates the plane equation of the light applied to the object based on the three-dimensional position information and posture information obtained from the second sensor and the second sensor, the plane equation of the entire object is calculated. Since the beam shape is switched by the beam shape switching means in the case of calculating the plane equation of the crack existing in the object, the crack existing in the measurement object can be accurately measured, and An image of a crack can be automatically synthesized in the measurement object.
以下、本発明を図に示した実施形態を用いて詳細に説明する。但し、この実施形態に記載される構成要素、種類、組み合わせ、形状、その相対配置などは特定的な記載がない限り、この発明の範囲をそれのみに限定する主旨ではなく単なる説明例に過ぎない。 Hereinafter, the present invention will be described in detail with reference to embodiments shown in the drawings. However, the components, types, combinations, shapes, relative arrangements, and the like described in this embodiment are merely illustrative examples and not intended to limit the scope of the present invention only unless otherwise specified. .
図1は本発明の実施形態に係る三次元計測装置を使用して対象物を計測している斜視図である。この三次元計測装置100は、対象物1にレーザスリット光またはレーザスポット光を照射するレーザ投光器(投光手段)20と、このレーザ投光器20により対象物1の表面に照射された光を撮像する撮像装置(撮像手段)30と、所定のエリアに磁界ベクトルを形成するトランスミッタ40と、図示しないパーソナルコンピュータ(PC)を備えて構成される。そして、レーザ投光器20はT字型の架台2と、スリット状のレーザ光またはレーザスポット光を出射するレーザ光源3と、レーザ投光器20の位置情報と姿勢情報を生成する磁気センサ4と、レーザ投光器20のビーム形状を切り替えるビーム形状切替器(ビーム形状切替手段)44と、を備えて構成され、レーザ光源3への電源と磁気センサ4の情報はケーブル5によりPCに伝えられる。また、撮像装置30はT字型の架台6と、対象物1に照射されたレーザ光を撮像するCCDカメラ7と、撮像装置30の位置情報と姿勢情報を生成する磁気センサ8とを備えて構成され、CCDカメラ7の信号と磁気センサ8の情報はケーブル9によりPCに伝えられる。尚、架台2、6は磁界の影響を磁気センサに与えないために、木製、プラスチック、ゴム等の材質により形成される。また、光源としてレーザ光を使用したが、LED或いは他の可視光を使用しても構わない。また、本実施形態ではレーザ投光器20と撮像装置30の位置情報と姿勢情報を生成する手段として、トランスミッタ40からの磁界ベクトルを受信するために磁気センサを使用したが、他の手段により生成しても構わない。ここで、投光手段は主としてレーザ投光器20により構成され、撮像手段は主として撮像装置30により構成され、第1のセンサは主として磁気センサ4により構成され、第2のセンサは主として磁気センサ8により構成され、方程式算出手段は主としてPCにより構成される。
FIG. 1 is a perspective view of measuring an object using a three-dimensional measuring apparatus according to an embodiment of the present invention. The three-
次に本実施形態の三次元計測装置100の概略動作について説明する。尚、対象物としてトンネルを例にして説明する。まず、トランスミッタ40をトンネル1の近傍に配置し、電源を投入してトランスミッタ40の磁界をレーザ投光器20及び撮像装置30が検出できるようにトンネル1との位置関係を設定する。これは図示しないPCの画面を見ながら調整を行う。次に、ビーム形状切替器44をスリット光側に切り替える。そしてレーザ投光器20と撮像装置30を手に持ち、PCの画面を見ながらCCDカメラ7がトンネル1の全体画像を捉えられるように配置する。そしてレーザ光源3の電源を投入して、スリット光がトンネル1の走査開始位置に来るようにしてPCに格納されたプログラムを開始する。そして、レーザ投光器20をトンネル1の走査開始位置から順次移動し、トンネル1の全体の走査が完了するまで行う。これらの走査の過程でCCDカメラ7から読み込まれたスリット光の平面方程式がリアルタイムに計算される。また、走査の過程でレーザ投光器20及び撮像装置30を移動しても、磁気センサ4、磁気センサ8から位置情報とロール角、ピッチ角、ヨー角の姿勢情報がPCに取り込まれリアルタイムに処理されてトンネル1の平面方程式に反映される。次にビーム形状切替器44をスポット光側に切り替え、そのスポット光を亀裂45に沿って走査する。
そしてトンネル1と亀裂45の走査が完了するとPCに取り込まれたデータに基づいて、トンネル1と亀裂45の三次元画像を再生する。この三次元画像はあらゆる角度からのデータにより構成されているので、対象物を回転して任意の角度から観察することも可能となる。
Next, a schematic operation of the three-
When the scanning of the tunnel 1 and the
図2は図1の斜視図に基づいて本実施形態の三次元計測装置100の全体構成を表す模式図である。同じ構成要素には同じ参照番号が付されているので、重複する説明は省略する。ここで、CCDカメラ7の信号はケーブル9によりイメージプロセッサ17に接続され、磁気センサ8の信号はケーブル9によりレーザ位置検出回路11に接続され、磁気センサ4はケーブル5によりレーザ位置検出回路11に接続され、トランスミッタ40はケーブル41によりレーザ位置検出回路11に接続されている。またレーザ位置検出回路11の出力及びイメージプロセッサ17の出力信号はPC12に入力され、PC12を制御するプログラムはROM(Read Only Memory)13に格納されている。またPC12からは処理された三次元画像を表示するモニタ18が接続されている。尚、イメージプロセッサ17はCCDカメラ7により撮像された画像データをPC12が処理し易いように変換する機能があり、CCDカメラ7からの画像信号からレーザ輝線の座標をリアルタイムに演算するために回路をFPGA(Field Programmable Gate Array)により構成している。またレーザ位置検出回路11は磁気センサ4、磁気センサ8の信号から三次元位置情報及び姿勢情報を検出回路である。
また説明の都合上CCDカメラ7の光軸15とレーザ光源3の光軸16のなす角度をαとすると、トンネル1の表面形状は滑らかであるので、角度αを大きくして一度に広い範囲の情報を取得して処理するものである。
FIG. 2 is a schematic diagram showing the overall configuration of the three-
For convenience of explanation, if the angle formed by the optical axis 15 of the
図3は本発明の磁気センサ4、磁気センサ8が位置情報と姿勢情報を生成する原理を説明する図である。図3(a)は磁界ベクトルと各磁気センサの位置関係を表す図であり、トランスミッタ40から放射された磁界ベクトル42が半球状に形成され、その磁界ベクトル42上に例えば磁気センサ4、磁気センサ8が存在し、対象物46が磁界ベクトル42の内側に含まれるものとする。図3(b)は磁界ベクトルを説明する模式図である。例えば磁気センサ4ではトランスミッタ40の位置における磁界の強さおよび方向をAとし、磁界の強さおよび方向AのX軸方向磁界の強さ、Y軸方向磁界の強さ、Z軸方向磁界の強さをそれぞれAx、Ay、Azとし、cosα、cosβ、cosγを磁界の強さおよび方向Aの方向余波とすれば、Ax=Acosα、Ay=Acosβ、Az=Acosγであり、磁気センサ4のX軸方向の磁界検出コイルからはAxの出力が、Y軸方向の磁界検出コイルからはAyの出力が、Z軸方向の磁界検出コイルからはAzの出力が送出される。そして磁界の強さおよび方向Aは、A=(Ax2+Ay2+Az2)1/2で与えられる。同じく磁気センサ8の磁界の強さおよび方向Bは、B=(Bx2+By2+Bz2)1/2で与えられる。
FIG. 3 is a diagram for explaining the principle by which the magnetic sensor 4 and the magnetic sensor 8 of the present invention generate position information and posture information. FIG. 3A is a diagram showing the positional relationship between the magnetic field vector and each magnetic sensor. The magnetic field vector 42 radiated from the
次にレーザ平面方程式の算出方法について説明する。
図4は各座標系の位置関係を表す図である。カメラ座標(画像モニタ上の座標u,v)50、CCDカメラ7に取り付けられた磁気センサ8のレシーバ座標51、レーザ光源3に取り付けられた磁気センサ4のレシーバ座標54、トランスミッタ40の中心を基準とした実座標であるワールド座標53をそれぞれ設定する。レーザ平面52の位置や傾きを自由に設定することを可能にするために、レーザ光源3に取り付けられた磁気センサ4が検知した情報(投光器20の位置と姿勢)を用いて、随時、レーザ平面52の方程式を算出することを行う。また、図7のように受像面55上にある計測点の位置Pとカメラの焦点Fを結ぶ直線lが磁気センサ8が検知した情報、(撮像装置30の位置と姿勢)を基に決定される。また、直線lとレーザ平面52の方程式の交点を求めることで、計測点のワールド座標53(X,Y,Z)を算出する。
Next, a method for calculating the laser plane equation will be described.
FIG. 4 is a diagram showing the positional relationship of each coordinate system. The camera coordinates (coordinates u, v on the image monitor) 50, the receiver coordinates 51 of the magnetic sensor 8 attached to the
図5はレシーバ座標とレーザ平面の関係を表す図である。
即ち、レーザ投光器20と撮像装置30取り付けられた磁気センサ4、8を用いてレーザ平面方程式をリアルタイムに算出することでレーザスリット平面の独立走査が可能となる。三次元磁気センサからはセンサ自体の位置(xow,yow,zow)と姿勢(φ,θ,ψ)のデータが得られる。三次元磁気センサのレシーバ座標に原点を置いたレシーバ座標上の点(xr,yr,zr)は、出力された三次元磁気センサのデータを基に(1)式によって、トランスミッタ40を原点としたワールド座標に変換される。
FIG. 5 is a diagram showing the relationship between the receiver coordinates and the laser plane.
That is, the laser slit plane can be independently scanned by calculating the laser plane equation in real time using the
ここで
C:Cos、S:Sinとする。
here
C: Cos, S: Sin.
レシーバ座標上でレーザ平面52上の任意の3点の座標(xrw,yrw,zrw)を式(1)に代入し、ワールド座標に変換することでレーザ平面52の方程式を得ることができる。つまり、
が得られるレーザ平面52の方程式である。
By substituting the coordinates (x rw , y rw , z rw ) of any three points on the
Is the equation of the
図6はカメラ座標とレシーバ座標の関係を示す図である。また、図7に二次元化したカメラ焦点と受像面及び計測点の関係を示す。カメラの焦点20から受像面までの距離をfと置くと、受像面上の座標(カメラ座標50)はカメラの焦点20を原点とした座標系において、
となる。
z=λとおき、線形化すると、
と表される。
FIG. 6 is a diagram showing the relationship between camera coordinates and receiver coordinates. FIG. 7 shows the relationship between the two-dimensional camera focus, the image receiving surface, and the measurement points. When the distance from the camera
It becomes.
When z = λ and linearized,
It is expressed.
ここで、式(6)を三次元磁気センサのレシーバを原点とした座標系(レシーバ座標54)で表すために、回転・平行移動を考慮して、以下のような行列で表すことができる。
ここでk11〜k33のパラメータにはCCDカメラ7の位置や姿勢などをはじめとする計測対象と、CCDカメラ7の位置を表すデータが全て含まれている。従って、式(7)がレシーバ座標(x,y,z)54とカメラ位置(u,v)の関係式になる。
Here, since Expression (6) is expressed by a coordinate system (receiver coordinates 54) with the receiver of the three-dimensional magnetic sensor as the origin, it can be expressed by the following matrix in consideration of rotation and translation.
Here, the parameters k 11 to k 33 include all the measurement objects including the position and orientation of the
式(7)を展開し、整理すると次式のように表せる。
尚、11個の未知数k11〜k33は、既知のワールド座標の基凖点(x,y,z)と、それに対応するカメラ座標21の点(u,v)の組み合わせを式(8)に代入し、連立方程式を解くことで求めることができる。
When formula (7) is expanded and arranged, it can be expressed as the following formula.
The eleven unknowns k 11 to k 33 are obtained by combining the known world coordinate base point (x, y, z) and the corresponding point (u, v) of the camera coordinate 21 with the formula (8). It can be obtained by substituting into and solving simultaneous equations.
また式(8)の2平面の交線で表される直線はカメラの焦点から計測点に向かう直線で以下のような媒介変数で表せる。
Further, the straight line represented by the intersection of the two planes in the equation (8) is a straight line from the focal point of the camera to the measurement point and can be represented by the following parameter.
式(9)はレシーバ座標54であるので、式(1)と同様にワールド座標53に変換する。
ここで、(wx,wy,wz)はトランスミッタ40を原点とした磁気センサ8のワールド座標である。
カメラの焦点から計測点に向かう直線である式(10)とレーザ平面上の方程式である式(4)を連立させ計測点21のワールド座標(X,Y,Z)を求めることができる。
Since the equation (9) is the receiver coordinate 54, it is converted into the world coordinate 53 as in the equation (1).
Here, (wx, wy, wz) is the world coordinates of the magnetic sensor 8 with the
The world coordinates (X, Y, Z) of the
図8はキャリブレーションの構成を示す図である。50mmおきに目盛りが書かれたキャリブレーションボード70、71をレシーバ72のXY平面と平行になるように設置し撮影する。そして、画面上での座標(u,v)をクリックし読み取る。また、対応するキャリブレーションボード上の座標(x,y,z)も読み取る。同じz軸方向距離で4点の座標を読み取り、それをz軸方向に50mm移動しながら複数回(6回以上)繰り返す。
式(8)に画面上での座標(u,v)と対応するキャリブレーションボード上の座標(x,y,z)の組み合わせを代入し、以下の連立方程式を組み立てる。
FIG. 8 is a diagram showing the configuration of calibration.
The following simultaneous equations are assembled by substituting a combination of coordinates (u, v) on the screen and coordinates (x, y, z) on the calibration board corresponding to equation (8).
ここで
here
そこで11個の係数
は擬似逆行列
を用いて算出する。
この行列式を解き係数
が求まると、カメラ座標とレシーバ座標の関係式(7)が決定される。
So 11 coefficients
Is the pseudoinverse
Calculate using.
Solve this determinant with coefficients
Is obtained, the relational expression (7) between the camera coordinates and the receiver coordinates is determined.
図9は本発明の本実施形態の三次元計測装置の動作を説明するフローチャートである。尚、対象物としてトンネルを例にして説明する。まず、トランスミッタ40をトンネル1の近傍に配置し(S1)、電源を投入してトランスミッタ40の磁界をレーザ投光器20及び撮像装置30が検出できるようにトンネル1との位置関係を設定する(S2)。これはPC12の画面を見ながら調整を行う。次に、ビーム形状切替器44をスリット光側に切り替える(S3)。そしてレーザ投光器20と撮像装置30を手に持ち、PC12の画面を見ながらCCDカメラ7がトンネル1の全体画像を捉えられるように配置する(S4)。そしてレーザ光源3の電源を投入して、スリット光がトンネル1の走査開始位置に来るようにしてPC12に格納されたプログラムを開始する(S5)。そして、レーザ投光器20をトンネル1の走査開始位置から順次移動する(S6)。これらの走査の過程でCCDカメラ7から読み込まれたスリット光の平面方程式がリアルタイムに計算される(S7)。そしてこの動作がトンネル1の全体の走査が完了するまで行う(S8)。また、走査の過程でレーザ投光器20及び撮像装置30を移動しても、磁気センサ4、磁気センサ8から位置情報とロール角、ピッチ角、ヨー角の姿勢情報がPCに取り込まれリアルタイムに処理されてトンネル1の平面方程式に反映される。
FIG. 9 is a flowchart for explaining the operation of the three-dimensional measuring apparatus according to this embodiment of the present invention. A description will be given by taking a tunnel as an example of the object. First, the
次にビーム形状切替器44をスポット光側に切り替え(S9)、そしてレーザ投光器20と撮像装置30を手に持ち、PC12の画面を見ながらCCDカメラ7が亀裂45の画像を捉えられるように配置する(S10)。そしてスポット光を亀裂45の走査開始位置に来るようにセットする(S11)。そして、レーザ投光器20を亀裂45に沿って順次移動する(S12)。これらの走査の過程でCCDカメラ7から読み込まれたスポット光の平面方程式がリアルタイムに計算される(S13)。そしてこの動作が亀裂45の全体の走査が完了するまで行う(S14)。走査が完了するとトンネル1と亀裂45の三次元形状が計算されてモニタ18上に三次元画像として再生される(S15)。
Next, the
以上の通り本発明によれば、物体に光を照射するレーザ投光器20と、このレーザ投光器20により物体表面に照射された光を撮像する撮像装置30と、レーザ投光器20に備えられこのレーザ投光器20の三次元位置情報及び姿勢情報を生成する磁気センサ4と、撮像装置30に備えられこの撮像装置30の三次元位置情報及び姿勢情報を生成する磁気センサ8と、所定のエリアに磁界ベクトルを形成するトランスミッタ40と、物体に照射された光の平面方程式を算出する方程式算出手段と、レーザ投光器20のビーム形状を切り替えるビーム形状切替手段44と、を備え、レーザ投光器20と撮像装置30がそれぞれ独立に動作することにより磁気センサ4及び磁気センサ8より得られる三次元位置情報及び姿勢情報に基づいて、方程式算出手段が物体に照射された光の平面方程式を算出する場合、物体全体の平面方程式を算出する場合と、物体に存在する亀裂の平面方程式を算出する場合とで、ビーム形状切替手段44によりビーム形状を切り替えるようにしたので、トンネル1に存在する亀裂45を精度良く計測することができ、且つトンネル1内に亀裂45の画像を自動的に合成することができる。
As described above, according to the present invention, the
また、撮像装置30の焦点から計測点に向かう直線である式(10)とレーザ平面上の方程式である式(4)を連立させ、計測点のワールド座標(X,Y,Z)を求めるので、座標変換の考え方を取り入れて平面方程式をリアルタイムに計算できるようにし、処理速度の高速化を図りレーザ投光器20の独立走査を行うことができる。
また、磁気センサ4と8から得られた位置情報と姿勢情報を基に、光の平面方程式をリアルタイムに算出することが可能となり、レーザ投光器20の独立走査を行うことができる。
Further, since the equation (10) that is a straight line from the focal point of the
Moreover, it becomes possible to calculate the plane equation of light in real time based on the position information and posture information obtained from the magnetic sensors 4 and 8, and the
また、レーザ投光器20の光源はレーザ光若しくはLEDを使用するので、光源と物体との距離及び装置のコストとのバランスにより何れかの光源を任意に選択することができる。
また、所定のエリアに磁界ベクトルを形成するトランスミッタ40を備えたので、各センサが磁界ベクトルを計算することにより、位置情報及び姿勢情報を正確に取得することができる。
Further, since the light source of the
Further, since the
1 トンネル、2、6 架台、3 レーザ光源、4、8 磁気センサ、7 CCDカメラ、11 レーザ位置検出回路、12 PC、13 ROM、17 イメージプロッセサ、18 モニタ、20 レーザ投光器、30 撮像装置、40 トランスミッタ DESCRIPTION OF SYMBOLS 1 Tunnel, 2 and 6 mounts, 3 Laser light source, 4 and 8 Magnetic sensor, 7 CCD camera, 11 Laser position detection circuit, 12 PC, 13 ROM, 17 Image processor, 18 Monitor, 20 Laser projector, 30 Imaging device, 40 Transmitter
Claims (6)
前記物体に光を照射する投光手段と、該投光手段により前記物体表面に照射された光を撮像する撮像手段と、前記投光手段に備えられ該投光手段の三次元位置情報及び姿勢情報を生成する第1のセンサと、前記撮像手段に備えられ該撮像手段の三次元位置情報及び姿勢情報を生成する第2のセンサと、所定のエリアに磁界ベクトルを形成するトランスミッタと、前記物体に照射された光の平面方程式を算出する方程式算出手段と、前記投光手段のビーム形状を切り替えるビーム形状切替手段と、を備え、
前記投光手段と撮像手段がそれぞれ独立に動作することにより前記第1のセンサ及び第2のセンサより得られる三次元位置情報及び姿勢情報に基づいて、前記方程式算出手段が前記物体に照射された光の平面方程式を算出する場合と、前記物体全体の平面方程式を算出する場合と、該物体に存在する亀裂の平面方程式を算出する場合とに応じて、前記ビーム形状切替手段により前記ビーム形状を切り替えるように構成したことを特徴とする三次元計測装置。 In a three-dimensional measurement apparatus that measures the three-dimensional shape of the object based on the three-dimensional information of the object,
Projecting means for irradiating the object with light, imaging means for imaging the light irradiated on the object surface by the projecting means, and three-dimensional position information and attitude of the projecting means provided in the projecting means A first sensor that generates information; a second sensor that is provided in the imaging unit and generates three-dimensional position information and orientation information of the imaging unit; a transmitter that forms a magnetic field vector in a predetermined area; and the object An equation calculating means for calculating a plane equation of the light irradiated to the light beam, and a beam shape switching means for switching the beam shape of the light projecting means,
The equation calculation unit is irradiated on the object based on the three-dimensional position information and the posture information obtained from the first sensor and the second sensor by the light projecting unit and the imaging unit operating independently. The beam shape is changed by the beam shape switching means according to a case where a plane equation of light is calculated, a case where a plane equation of the entire object is calculated, and a case where a plane equation of a crack existing in the object is calculated. A three-dimensional measuring apparatus characterized by being configured to be switched.
ここで
C:Cos、S:Sinとすると、
前記センサ座標上で光平面上の任意の3点の座標(xr,yr,zr)を前記式(1)及び(2)により前記トランスミッタを原点としたワールド座標に変換し、前記ワールド座標(xrw,yrw,zrw)に変換することにより、
前記式(3)及び式(4)として光の平面方程式を取得し、前記撮像手段の焦点から受像面までの距離をfと置くと、前記受像面上の座標(u,v)は前記撮像手段の焦点を原点とした座標系において、
となり、Z=λとおき線形化すると、
と表わされ、該式(6)を前記センサ座標で表すために、回転・平行移動を考慮して、
式(7)として表し、(ここでk11〜k33のパラメータには前記撮像手段の位置や姿勢などをはじめとする計測対象と、前記撮像手段の位置を表すデータが全て含まれているものとする)該式(7)が前記センサ座標(x,y,z)と前記受像面上の座標(u,v)の関係式となり、前記式(7)を展開し、整理すると式(8)を得て、
前記式(8)の2平面の交線で表される直線は前記撮像手段の焦点から計測点に向かう直線として
式(9)のレシーバ座標を求め、前記式(1)と同様にワールド座標に変換して前記撮像手段の焦点から計測点に向かう直線である式(10)を求め、
(ここで、(wx,wy,wz)はトランスミッタを原点としたレシーバのワールド座標である)
前記式(10)とレーザ平面上の方程式である前記式(4)を連立させてワールド座標(X,Y,Z)を求めることを特徴とする請求項1又は2に記載の三次元計測装置。 The equation calculation means, when the position data obtained from the first and second sensors is (x ow , y ow , z ow ) and the posture data is (φ, θ, ψ), The point (x r , y r , z r ) on the sensor coordinates with the origin at the sensor coordinates of the sensor is based on the output data of each sensor.
here
C: Cos, S: Sin,
The coordinates (x r , y r , z r ) of any three points on the optical plane on the sensor coordinates are converted into world coordinates with the transmitter as the origin by the equations (1) and (2), and the world By converting to coordinates (x rw , y rw , z rw ),
When the plane equation of light is acquired as the equations (3) and (4), and the distance from the focal point of the imaging means to the image receiving surface is set to f, the coordinates (u, v) on the image receiving surface are the imaging In the coordinate system with the focal point of the means as the origin,
And when linearizing with Z = λ,
In order to express the equation (6) in the sensor coordinates, in consideration of rotation / translation,
Expressed as Equation (7), (wherein k 11 and the measurement object, including position and orientation of the imaging means to the parameters to k 33, those data representing the position of the imaging means contains all The expression (7) becomes a relational expression between the sensor coordinates (x, y, z) and the coordinates (u, v) on the image receiving surface. When the expression (7) is developed and rearranged, the expression (8) )
The straight line represented by the intersection of the two planes in the equation (8) is a straight line from the focal point of the imaging means to the measurement point.
Obtain the receiver coordinates of Equation (9), convert to world coordinates as in Equation (1), and obtain Equation (10), which is a straight line from the focal point of the imaging means to the measurement point,
(Where (wx, wy, wz) are the world coordinates of the receiver with the transmitter as the origin)
The three-dimensional measuring apparatus according to claim 1 or 2, wherein the world coordinates (X, Y, Z) are obtained by simultaneously combining the expression (10) and the expression (4) which is an equation on a laser plane. .
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2006188705A JP2008014882A (en) | 2006-07-07 | 2006-07-07 | Three-dimensional measuring device |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2006188705A JP2008014882A (en) | 2006-07-07 | 2006-07-07 | Three-dimensional measuring device |
Publications (1)
Publication Number | Publication Date |
---|---|
JP2008014882A true JP2008014882A (en) | 2008-01-24 |
Family
ID=39072027
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2006188705A Pending JP2008014882A (en) | 2006-07-07 | 2006-07-07 | Three-dimensional measuring device |
Country Status (1)
Country | Link |
---|---|
JP (1) | JP2008014882A (en) |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2009192483A (en) * | 2008-02-18 | 2009-08-27 | Toyota Motor Corp | Three dimensional shape measuring method and three dimensional shape measuring device |
JP2010281811A (en) * | 2009-05-08 | 2010-12-16 | Univ Of Miyazaki | Device for making temperature image three dimensional |
CN102589530A (en) * | 2012-02-24 | 2012-07-18 | 合肥工业大学 | Method for measuring position and gesture of non-cooperative target based on fusion of two dimension camera and three dimension camera |
KR101217186B1 (en) | 2011-11-23 | 2013-01-02 | 한국과학기술원 | Displacement measurement system of structure and method of the same |
CN103471519A (en) * | 2013-09-27 | 2013-12-25 | 国家电网公司 | Method for measuring deformation of power transmission and transformation tower by adoption of prism-free photoelectric total station |
CN106019398A (en) * | 2016-05-26 | 2016-10-12 | 武汉固德超前高新科技研发有限公司 | Detecting device, detecting system and detecting method for exploring internal structure of cave |
CN106600604A (en) * | 2016-11-14 | 2017-04-26 | 电子科技大学 | Threshold optimization method of magneto-optic imaging defect detection method |
CN108802043A (en) * | 2017-05-11 | 2018-11-13 | 成都中信华瑞科技有限公司 | Tunnel detector, detecting system and tunnel defect information extracting method |
-
2006
- 2006-07-07 JP JP2006188705A patent/JP2008014882A/en active Pending
Cited By (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2009192483A (en) * | 2008-02-18 | 2009-08-27 | Toyota Motor Corp | Three dimensional shape measuring method and three dimensional shape measuring device |
JP2010281811A (en) * | 2009-05-08 | 2010-12-16 | Univ Of Miyazaki | Device for making temperature image three dimensional |
KR101217186B1 (en) | 2011-11-23 | 2013-01-02 | 한국과학기술원 | Displacement measurement system of structure and method of the same |
WO2013077623A1 (en) * | 2011-11-23 | 2013-05-30 | 한국과학기술원 | Structure displacement measurement system and method |
CN102589530A (en) * | 2012-02-24 | 2012-07-18 | 合肥工业大学 | Method for measuring position and gesture of non-cooperative target based on fusion of two dimension camera and three dimension camera |
CN103471519A (en) * | 2013-09-27 | 2013-12-25 | 国家电网公司 | Method for measuring deformation of power transmission and transformation tower by adoption of prism-free photoelectric total station |
CN103471519B (en) * | 2013-09-27 | 2016-09-07 | 国家电网公司 | A kind of method applied without the deformation of prism photoelectricity total station survey power transmission and transformation pylon |
CN106019398A (en) * | 2016-05-26 | 2016-10-12 | 武汉固德超前高新科技研发有限公司 | Detecting device, detecting system and detecting method for exploring internal structure of cave |
CN106600604A (en) * | 2016-11-14 | 2017-04-26 | 电子科技大学 | Threshold optimization method of magneto-optic imaging defect detection method |
CN106600604B (en) * | 2016-11-14 | 2019-06-14 | 电子科技大学 | A kind of threshold values optimization method of magneto-optic image defects detection method |
CN108802043A (en) * | 2017-05-11 | 2018-11-13 | 成都中信华瑞科技有限公司 | Tunnel detector, detecting system and tunnel defect information extracting method |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
EP2132523B1 (en) | Method and device for exact measurement of objects | |
JP4885584B2 (en) | Rangefinder calibration method and apparatus | |
US9500469B2 (en) | Laser line probe having improved high dynamic range | |
JP2008014882A (en) | Three-dimensional measuring device | |
KR20150128300A (en) | method of making three dimension model and defect analysis using camera and laser scanning | |
JP5260175B2 (en) | Non-contact measuring method and measuring system for measuring coordinates of target surface | |
JP2007120993A (en) | Object shape measuring device | |
JP2006064690A (en) | In-tubing profile measuring device | |
US20150085108A1 (en) | Lasergrammetry system and methods | |
JP5522630B2 (en) | 3D system | |
JP2007322162A (en) | Three-dimensional shape measuring apparatus and three-dimensional shape measuring method | |
JP2013064644A (en) | Shape-measuring device, shape-measuring method, system for manufacturing structures, and method for manufacturing structures | |
JP2001148025A5 (en) | ||
JP2014153149A (en) | Shape measurement device, structure manufacturing system, shape measurement method and program | |
CN113658270B (en) | Method, device, medium and system for multi-vision calibration based on workpiece hole center | |
JP4591103B2 (en) | X-ray CT inspection apparatus and X-ray CT inspection method | |
JP3991040B2 (en) | Three-dimensional measuring apparatus and three-dimensional measuring method | |
WO2019230171A1 (en) | Laser calibration device, calibration method therefor, and image input device including laser calibration device | |
JP4778855B2 (en) | Optical measuring device | |
JP2000039310A (en) | Method and device for measuring shape | |
JP2007033040A (en) | Method and device for calibrating optical head part in three-dimensional shape measuring instrument by optical cutting method | |
JP2021527217A (en) | Alignment method for the beam alignment unit of the interference measuring device and the measuring device for performing the interference measurement using laser irradiation. | |
JP2011138096A (en) | Measuring microscope | |
JP2009264898A (en) | Workpiece position/attitude measurement method and measuring device | |
JP2012013593A (en) | Calibration method for three-dimensional shape measuring machine, and three-dimensional shape measuring machine |