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
JP3208900B2 - Method and apparatus for recognizing three-dimensional position and orientation based on vision - Google Patents
[go: Go Back, main page]

JP3208900B2 - Method and apparatus for recognizing three-dimensional position and orientation based on vision - Google Patents

Method and apparatus for recognizing three-dimensional position and orientation based on vision

Info

Publication number
JP3208900B2
JP3208900B2 JP04924093A JP4924093A JP3208900B2 JP 3208900 B2 JP3208900 B2 JP 3208900B2 JP 04924093 A JP04924093 A JP 04924093A JP 4924093 A JP4924093 A JP 4924093A JP 3208900 B2 JP3208900 B2 JP 3208900B2
Authority
JP
Japan
Prior art keywords
error
transformation matrix
perspective transformation
target object
unit
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 - Lifetime
Application number
JP04924093A
Other languages
Japanese (ja)
Other versions
JPH06258028A (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.)
Denso Corp
Original Assignee
Denso Corp
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 Denso Corp filed Critical Denso Corp
Priority to JP04924093A priority Critical patent/JP3208900B2/en
Priority to US08/206,828 priority patent/US5499306A/en
Publication of JPH06258028A publication Critical patent/JPH06258028A/en
Application granted granted Critical
Publication of JP3208900B2 publication Critical patent/JP3208900B2/en
Anticipated expiration legal-status Critical
Expired - Lifetime legal-status Critical Current

Links

Landscapes

  • Length Measuring Devices By Optical Means (AREA)
  • Manipulator (AREA)
  • Image Processing (AREA)
  • Control Of Position Or Direction (AREA)
  • Image Analysis (AREA)
  • Measurement Of Optical Distance (AREA)
  • Control Of Position, Course, Altitude, Or Attitude Of Moving Bodies (AREA)

Description

【発明の詳細な説明】DETAILED DESCRIPTION OF THE INVENTION

【0001】[0001]

【産業上の利用分野】本発明は、自律移動ロボットなど
に搭載され単眼視カメラにより撮影した既知対象の画像
情報に基いて単眼視カメラの三次元位置および姿勢を認
識するようにした視覚に基く三次元位置および姿勢の認
識方法とその装置に関する。
BACKGROUND OF THE INVENTION 1. Field of the Invention The present invention is based on vision that recognizes the three-dimensional position and orientation of a monocular camera based on image information of a known object captured by a monocular camera mounted on an autonomous mobile robot or the like. The present invention relates to a three-dimensional position and orientation recognition method and apparatus.

【0002】[0002]

【従来の技術】近年、例えば、自律移動ロボット等のよ
うに、形状および絶対座標位置が既知である対象物体に
対して接近するように移動し、アームにより対象物体を
把持するといった能動的操作を行うようにしたものがあ
る。このような自律移動ロボット等においては、CCD
カメラなどの撮像手段を搭載してその二次元的な画像情
報に基いて対象物体に対する自己の三次元的な相対位置
を認識することが必要になる。この場合、位置認識を行
うための装置においては、高精度且つ迅速に演算を行っ
て三次元的な位置のデータを求めることが要求され、し
かも安価に実現できる構成が望まれている。
2. Description of the Related Art In recent years, active operations such as moving an object such as an autonomous mobile robot so as to approach a target object having a known shape and absolute coordinate position and grasping the target object by an arm are performed. There's something we did. In such autonomous mobile robots, CCDs
It is necessary to mount imaging means such as a camera and recognize its own three-dimensional relative position with respect to the target object based on the two-dimensional image information. In this case, in an apparatus for performing position recognition, it is required to perform calculations with high accuracy and speed to obtain three-dimensional position data, and a configuration that can be realized at low cost is desired.

【0003】そこで、従来では、比較的簡単なものとし
て、例えば、特開平3−166072号公報あるいは特
開平3−166073号公報に示されるようなものが考
えられている。すなわち、これらにおいては、対象装置
に固定された形状および寸法が既知の特殊な幾何学形状
(マーカと呼ぶ)をカメラにより撮影し、そのマーカの
二次元的な画像情報に基いて、その重心を計算する等の
方法により対象装置に対するカメラの相対的な三次元の
位置関係を求めるようにしたものである。
[0003] Conventionally, as a relatively simple device, for example, a device as disclosed in Japanese Patent Application Laid-Open No. Hei 3-166072 or Japanese Patent Application Laid-Open No. Hei 3-166073 has been considered. That is, in these, a special geometric shape (referred to as a marker) having a known shape and dimensions fixed to the target device is photographed by a camera, and the center of gravity is determined based on the two-dimensional image information of the marker. The relative three-dimensional positional relationship of the camera with respect to the target device is obtained by a method such as calculation.

【0004】[0004]

【発明が解決しようとする課題】しかしながら、上述の
ような従来方法によると、対象装置にマーカを設ける必
要があるため、その設置のための準備をその都度行う必
要があると共に、マーカに対してカメラの光軸を常に垂
直に設定する必要があるため、任意の姿勢をとるカメラ
に対しては適用が困難となる不具合がある。
However, according to the above-described conventional method, since it is necessary to provide a marker on the target device, it is necessary to prepare for the installation each time, and the marker must be provided for the marker. Since it is necessary to always set the optical axis of the camera vertically, there is a problem that it is difficult to apply the camera to an arbitrary posture.

【0005】このような不具合を解消し得る理論的手法
として、一般には、透視n点問題において最小二乗推定
した透視変換行列の成分からカメラの位置および姿勢を
表すパラメータを逆算する方法がコンピュータビジョン
の手法として存在する。この方法によれば、カメラの光
軸と対象物体を記述する絶対座標の原点との間の関係が
未知の場合でも、原理的にはカメラの位置および姿勢を
表すパラメータを求めることが可能となる。
As a theoretical method capable of solving such a problem, a method of inversely calculating parameters representing the position and orientation of a camera from components of a perspective transformation matrix estimated by least squares in a perspective n-point problem is generally used in computer vision. It exists as a method. According to this method, even when the relationship between the optical axis of the camera and the origin of the absolute coordinates describing the target object is unknown, it is possible in principle to obtain parameters representing the position and orientation of the camera. .

【0006】ところが、このようなコンピュータビジョ
ンの手法においては、実際にコンピュータを用いて明示
的に計算する場合に、適用する式が非常に複雑になるこ
とに加えてその演算量が膨大なものとなるため、迅速な
演算処理が困難になり実用的には適用が難しくなるもの
である。
However, in such a computer vision method, when an explicit calculation is actually performed using a computer, the applied formula becomes extremely complicated and the amount of calculation is enormous. Therefore, rapid arithmetic processing becomes difficult, and practical application becomes difficult.

【0007】本発明は、上記事情に鑑みてなされたもの
で、その目的は、カメラの光軸と既知形状の対象物体を
記述する絶対座標の原点との間の関係が未知の場合で
も、カメラにより撮影した二次元の画像情報に基いて、
比較的簡単かつ迅速ににカメラの位置と姿勢を表すパラ
メータとを求めることができるようにした視覚情報に基
く三次元位置と姿勢の認識方法とその装置を提供するこ
とにある。
SUMMARY OF THE INVENTION The present invention has been made in view of the above circumstances, and an object of the present invention is to provide a camera which can be used even when the relationship between the optical axis of the camera and the origin of absolute coordinates describing a target object having a known shape is unknown. Based on the two-dimensional image information taken by
It is an object of the present invention to provide a method and apparatus for recognizing a three-dimensional position and orientation based on visual information, which can relatively easily and quickly obtain a parameter representing the position and orientation of a camera.

【0008】[0008]

【課題を解決するための手段】本発明の視覚に基く三次
元位置および姿勢の認識方法は、形状を示す三次元絶対
座標が既知である対象物体を撮像手段により撮影したと
きの画像情報に基いてその対象物体の形状の特徴を示す
二次元の画像特徴点を抽出する画像処理工程と、前記撮
像手段の概略位置と姿勢とに基いて、前記画像処理工程
で抽出した画像特徴点と記憶手段にあらかじめ記憶され
た前記対象物体の見え方のパターンとの間の特徴点の対
応リストを選択すると共に、その対応リストに基いて前
記撮像手段の現在位置と姿勢とにおける前記対象物体の
見え方を表すカメラ行列の最小二乗推定を行って推定透
視変換行列を求める最小二乗推定工程と、前記撮像手段
と前記対象物体との概略的な推定距離に応じて必要な演
算精度に対応する密度の量子化空間点を有する探索領域
を設定する探索領域設定工程と、前記探索領域の量子化
空間点において前記撮像手段により前記対象物体を撮影
したときに得られるべき二次元の画像特徴点を算出する
ための透視変換行列を、あらかじめ記憶手段に記憶され
た前記対象物体の三次元絶対座標に対する前記量子化空
間点の推定位置と前記撮像手段の撮像姿勢を示す姿勢パ
ラメータとに基いて演算する透視変換行列演算工程と、
この透視変換行列演算工程にて求められた前記透視変換
行列と前記最小二乗推定工程において求められた推定透
視変換行列との間の適合度を示す評価誤差値を所定の誤
差評価関数に基いて演算する誤差演算工程と、この誤差
演算工程にて求められた前記評価誤差値を前記探索領域
に対応して設定された要求誤差精度と比較し、前記評価
誤差値が前記要求誤差精度よりも小さくなることを判定
条件として、前記探索領域設定工程,前記透視変換行列
演算工程および前記誤差演算工程を繰り返し実行させる
誤差判定工程と、この誤差判定工程において前記判定条
件を満たしたときに、そのときの量子化空間点の三次元
絶対座標と前記撮像手段の姿勢パラメータとを出力する
データ出力工程とから構成したところに特徴を有する。
According to the present invention, there is provided a method for recognizing a three-dimensional position and orientation based on vision, based on image information obtained by photographing a target object whose shape has known three-dimensional absolute coordinates by an image pickup means. An image processing step of extracting two-dimensional image feature points indicating characteristics of the shape of the target object; and an image feature point extracted in the image processing step based on the approximate position and orientation of the imaging means, and a storage means. A correspondence list of feature points between the pattern of the appearance of the target object stored in advance is selected, and based on the correspondence list, the appearance of the target object at the current position and orientation of the imaging unit is determined. A least-squares estimation step of obtaining an estimated perspective transformation matrix by performing a least-squares estimation of the camera matrix to be represented, and corresponding to a required calculation accuracy according to a rough estimated distance between the imaging unit and the target object. A search area setting step of setting a search area having a quantization space point of a degree, and a two-dimensional image feature point to be obtained when the target object is photographed by the imaging means at the quantization space point of the search area. A perspective transformation matrix for calculation is calculated based on an estimated position of the quantization space point with respect to the three-dimensional absolute coordinates of the target object stored in advance in the storage unit and a posture parameter indicating an imaging posture of the imaging unit. A perspective transformation matrix operation step,
An evaluation error value indicating the degree of conformity between the perspective transformation matrix calculated in the perspective transformation matrix calculation step and the estimated perspective transformation matrix calculated in the least square estimation step is calculated based on a predetermined error evaluation function. Error calculating step, and comparing the evaluation error value obtained in the error calculating step with a required error precision set corresponding to the search area, and the evaluation error value becomes smaller than the required error precision. As a determination condition, an error determination step in which the search area setting step, the perspective transformation matrix calculation step, and the error calculation step are repeatedly executed, and when the determination condition is satisfied in the error determination step, the quantum at that time is determined. It is characterized in that it comprises a data output step of outputting the three-dimensional absolute coordinates of the image space point and the attitude parameter of the imaging means.

【0009】また、本発明の視覚に基く三次元位置およ
び姿勢の認識装置は、画像情報を出力する撮像手段と、
形状を示す三次元絶対座標が既知である対象物体を前記
撮像手段により撮影したときの画像情報に基いてその対
象物体の形状の特徴を示す二次元の画像特徴点を抽出す
る画像処理手段と、前記撮像手段の概略位置と姿勢とに
基いて、前記画像処理手段により抽出した画像特徴点と
記憶手段にあらかじめ記憶された前記対象物体の見え方
のパターンとの間の特徴点の対応リストを選択すると共
に、その対応リストに基いて前記撮像手段の現在位置と
姿勢とにおける前記対象物体の見え方を表すカメラ行列
の最小二乗推定を行って推定透視変換行列を求める最小
二乗推定手段と、前記撮像手段と前記対象物体との概略
的な推定距離に応じて必要な演算精度に対応する密度の
量子化空間点を有する探索領域を設定する探索領域設定
手段と、前記探索領域の量子化空間点において前記撮像
手段により前記対象物体を撮影したときに得られるべき
二次元の画像特徴点を算出するための透視変換行列を、
あらかじめ記憶手段に記憶された前記対象物体の三次元
絶対座標に対する前記量子化空間点の推定位置と前記撮
像手段の撮像姿勢を示す姿勢パラメータとに基いて演算
する透視変換行列演算手段と、この透視変換行列演算手
段により求められた前記透視変換行列と前記最小二乗推
定手段により求められた推定透視変換行列との間の適合
度を示す評価誤差値を所定の誤差評価関数に基いて演算
する誤差演算手段と、この誤差演算手段にて求められた
前記評価誤差値を前記探索領域に対応して設定された要
求誤差精度と比較し、前記評価誤差値が前記要求誤差精
度よりも小さくなることを判定条件として、前記探索領
域設定手段,前記透視変換行列演算手段および前記誤差
演算手段に前記演算を繰り返し実行させる誤差判定手段
と、この誤差判定手段により前記判定条件が満たされた
ときに、そのときの量子化空間点の三次元絶対座標と前
記撮像手段の姿勢パラメータとを出力するデータ出力手
段とから構成したところに特徴を有する。
In addition, the apparatus for recognizing a three-dimensional position and orientation based on vision according to the present invention comprises: an image pickup means for outputting image information;
Image processing means for extracting a two-dimensional image feature point indicating a characteristic of the shape of the target object based on image information when the target object whose three-dimensional absolute coordinates indicating the shape is known is photographed by the imaging means, A correspondence list of feature points between the image feature points extracted by the image processing means and the pattern of the appearance of the target object stored in advance in the storage means is selected based on the approximate position and orientation of the imaging means. And a least-squares estimating means for performing a least-squares estimation of a camera matrix representing the appearance of the target object at the current position and orientation of the imaging means based on the correspondence list to obtain an estimated perspective transformation matrix; Search area setting means for setting a search area having a quantized space point having a density corresponding to a required calculation accuracy according to a rough estimated distance between the means and the target object; The perspective transformation matrix for calculating the two-dimensional image feature points to be obtained when photographing the object by the imaging unit in the quantization space point range,
Perspective transformation matrix computing means for performing computation based on the estimated position of the quantization space point with respect to the three-dimensional absolute coordinates of the target object stored in advance in the storage means and a posture parameter indicating the imaging posture of the imaging means; Error calculation for calculating an evaluation error value indicating a degree of conformity between the perspective transformation matrix obtained by the transformation matrix calculation means and the estimated perspective transformation matrix obtained by the least-squares estimation means based on a predetermined error evaluation function Means, and comparing the evaluation error value obtained by the error calculation means with a required error accuracy set corresponding to the search area, and determining that the evaluation error value is smaller than the required error accuracy. Error determination means for causing the search area setting means, the perspective transformation matrix calculation means, and the error calculation means to repeatedly execute the calculation; When the judgment condition is satisfied by the step, characterized in was constructed from the data output means for outputting a three-dimensional absolute coordinates of the quantized spatial points at that time and attitude parameters of the imaging means.

【0010】[0010]

【作用】本発明の視覚に基く三次元位置および姿勢の認
識方法によれば、対象物体に対して移動接近して作業を
行う装置等に搭載した撮像手段により、三次元絶対座標
が既知である対象物体が撮影されると、画像処理工程に
おいては、その画像情報に基いて対象物体の形状の特徴
を示す二次元の画像特徴点を抽出する。
According to the method for recognizing a three-dimensional position and orientation based on vision of the present invention, the three-dimensional absolute coordinates are known by the imaging means mounted on a device or the like which moves and approaches the target object. When the target object is photographed, in the image processing step, a two-dimensional image feature point indicating the shape characteristic of the target object is extracted based on the image information.

【0011】続いて、最小二乗推定工程においては、前
記撮像手段の概略位置と姿勢とに基いて、前記画像処理
工程で抽出した画像特徴点と記憶手段にあらかじめ記憶
された前記対象物体の見え方のパターンとの間の特徴点
の対応リストを選択し、この後、その対応リストに基い
て前記撮像手段の現在位置と姿勢とにおいて前記対象物
体の見え方を表す変換行列としてのカメラ行列の最小二
乗推定を行って推定透視変換行列を求める。
Subsequently, in the least squares estimation step, the image feature points extracted in the image processing step and the appearance of the target object stored in the storage means in advance based on the approximate position and orientation of the imaging means. A correspondence list of feature points between the pattern and the corresponding pattern is selected. Thereafter, based on the correspondence list, the minimum of a camera matrix as a transformation matrix representing the appearance of the target object at the current position and orientation of the imaging unit is selected. Square estimation is performed to obtain an estimated perspective transformation matrix.

【0012】次に、探索領域設定工程において、撮像手
段とその対象物体との概略的な推定距離に応じて必要な
演算精度に対応する密度の量子化空間点を有する探索領
域を設定し、続く透視変換行列演算工程において、その
探索領域の量子化空間点に前記撮像手段を配置して前記
対象物体を撮影したときに得られるべき二次元の画像特
徴点を算出するための透視変換行列を、あらかじめ記憶
手段に記憶された前記対象物体の三次元絶対座標に対す
る前記量子化空間点の推定位置と前記撮像手段の撮像姿
勢を示す姿勢パラメータとに基いて演算する。
Next, in a search area setting step, a search area having quantized spatial points of a density corresponding to a required calculation accuracy is set according to the approximate estimated distance between the imaging means and the target object, and then the search area is set. In the perspective transformation matrix calculation step, a perspective transformation matrix for calculating a two-dimensional image feature point to be obtained when the imaging unit is arranged at the quantization space point of the search area and the target object is photographed, The calculation is performed based on the estimated position of the quantization space point with respect to the three-dimensional absolute coordinates of the target object stored in advance in the storage unit and a posture parameter indicating the imaging posture of the imaging unit.

【0013】そして、誤差演算工程において、透視変換
行列演算工程にて求められた前記透視変換行列と前記最
小二乗推定工程において求められた推定透視変換行列と
の間の適合度を示す評価誤差値を所定の誤差評価関数に
基いて演算する。これは、各行列における対応する要素
のずれを演算するもので、評価誤差値が小さい程両行列
間の適合度合が高いことを示している。
In the error calculation step, an evaluation error value indicating the degree of conformity between the perspective transformation matrix obtained in the perspective transformation matrix calculation step and the estimated perspective transformation matrix obtained in the least square estimation step is calculated. The calculation is performed based on a predetermined error evaluation function. This calculates the shift of the corresponding element in each matrix, and indicates that the smaller the evaluation error value, the higher the degree of matching between the two matrices.

【0014】次に、誤差判定工程において、上述の誤差
演算工程にて求められた前記評価誤差値を前記探索領域
に対応して設定された要求誤差精度と比較し、前記評価
誤差値が前記要求誤差精度よりも小さくなることを判定
条件として、前記探索領域設定工程,前記透視変換行列
演算工程および前記誤差演算工程を繰り返し実行させ、
最終的に、誤差判定工程にて上記条件に適合するように
なったときには、データ出力工程において、そのときの
量子化空間点の三次元絶対座標と前記撮像手段の姿勢パ
ラメータとを出力するようになり、実際の撮像手段の三
次元絶対座標に最も近い量子化空間点を近似点として求
めることができると共に、撮像手段の光軸のずれ等を表
す姿勢パラメータを検出することができるようになり、
これらのデータに基いて対象物体に対する装置の移動制
御などを適確に行うことができるようになるものであ
る。
Next, in the error judging step, the evaluation error value obtained in the above-described error calculating step is compared with a required error accuracy set corresponding to the search area, and the evaluation error value The search area setting step, the perspective transformation matrix calculation step, and the error calculation step are repeatedly executed, with the determination condition that the difference is smaller than the error precision,
Finally, when the above condition is met in the error determination step, in the data output step, the three-dimensional absolute coordinates of the quantization space point at that time and the posture parameter of the imaging unit are output. Thus, the quantization space point closest to the three-dimensional absolute coordinates of the actual imaging unit can be obtained as an approximate point, and the posture parameter representing the shift of the optical axis of the imaging unit can be detected.
Based on these data, movement control of the apparatus with respect to the target object can be performed accurately.

【0015】また、本発明の視覚に基く三次元位置およ
び姿勢の認識装置によれば、上述の各工程を実施する手
段により、同様の動作を行なわせることができ、実際の
撮像手段の三次元絶対座標に最も近い量子化空間点を近
似点として求めることができると共に、撮像手段の光軸
のずれ等を表す姿勢パラメータを検出することができる
ようになり、これらのデータに基いて対象物体に対する
装置の移動制御などを適確に行うことができるようにな
るものである。
Further, according to the apparatus for recognizing a three-dimensional position and orientation based on vision of the present invention, the same operation can be performed by means for performing the above-mentioned steps, and the three-dimensional position of the actual imaging means can be controlled. The quantization space point closest to the absolute coordinates can be obtained as an approximate point, and the posture parameter representing the shift of the optical axis of the imaging means can be detected. This makes it possible to accurately control the movement of the device.

【0016】[0016]

【実施例】以下、本発明を工場内を走行する自律移動ロ
ボットに搭載する三次元位置認識装置に適用した場合の
一実施例について図面を参照しながら説明する。全体構
成の概略を示す図2において、自律移動ロボット1は、
例えば、工場内などの床面に設けられた軌道に沿って作
業エリア2間を移動するようになっており、本発明の三
次元位置認識装置(図1参照)が搭載されている。自律
移動ロボット1の上部にはアーム3が設けられ、そのア
ーム3の先端部には撮像手段としてのカメラ4および対
象物体としてのワーク5を把持する等の作業を行うため
のハンド6が配設されている。この場合、作業エリア2
は壁面から窪んだ状態に形成された凹部に設けられ、所
定形状をなすワーク5はその作業エリア2の底面2aの
所定位置に載置されている。
DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS An embodiment in which the present invention is applied to a three-dimensional position recognition apparatus mounted on an autonomous mobile robot running in a factory will be described below with reference to the drawings. In FIG. 2 showing an outline of the entire configuration, an autonomous mobile robot 1
For example, it moves between work areas 2 along a track provided on a floor surface in a factory or the like, and is equipped with a three-dimensional position recognition device (see FIG. 1) of the present invention. An arm 3 is provided on the upper part of the autonomous mobile robot 1, and a camera 4 as an imaging means and a hand 6 for performing work such as grasping a work 5 as a target object are provided at a tip of the arm 3. Have been. In this case, work area 2
The work 5 having a predetermined shape is placed at a predetermined position on the bottom surface 2 a of the work area 2.

【0017】さて、三次元位置認識装置は図1に示すよ
うに構成される。画像処理手段としての画像処理部7
は、前処理計算部8および特徴点抽出部9からなる。前
処理計算部8は、カメラ4が撮影したワーク5を含むシ
ーンの画像信号を入力して前処理を施し、特徴点抽出部
9は、前処理計算部8により計算された画像情報に基い
て、後述のように、ワーク5の特徴点に対応する二次元
の画像特徴点を抽出するようになっている。
Now, the three-dimensional position recognition apparatus is configured as shown in FIG. Image processing unit 7 as image processing means
Comprises a pre-processing calculation unit 8 and a feature point extraction unit 9. The pre-processing calculation unit 8 inputs an image signal of a scene including the work 5 captured by the camera 4 and performs pre-processing, and the feature point extraction unit 9 performs processing based on the image information calculated by the pre-processing calculation unit 8. As will be described later, a two-dimensional image feature point corresponding to the feature point of the work 5 is extracted.

【0018】最小二乗推定手段としての最小二乗推定部
10は、画像処理部7にて抽出された画像特徴点が与え
られると、これに基いて後述するようにして、カメラ4
によりワーク5を撮影したときにそのカメラ4の位置お
よび姿勢に対応して得られるべき見え方のパターンを選
択すると共に、カメラ4の現在位置と姿勢とに対するカ
メラ行列の最小二乗推定を行い、この結果得られる推定
透視変換行列Msを演算して求める。
When the image feature points extracted by the image processing unit 7 are given, the least squares estimation unit 10 as a least squares estimation unit, based on the image feature points,
, The appearance pattern to be obtained corresponding to the position and orientation of the camera 4 when the work 5 is photographed is selected, and the least square estimation of the camera matrix for the current position and orientation of the camera 4 is performed. The resulting estimated perspective transformation matrix Ms is calculated and obtained.

【0019】量子化空間探索部11は、探索領域設定手
段としての探索領域設定部12,透視変換行列演算手段
としての透視変換行列演算部13,誤差演算手段として
の誤差演算部14,誤差判定手段としての誤差判定部1
5およびデータ出力手段としてのデータ出力部16から
構成されている。
The quantization space search section 11 includes a search area setting section 12 as search area setting means, a perspective transformation matrix operation section 13 as perspective transformation matrix operation means, an error operation section 14 as error operation means, and an error determination means. Error determination unit 1 as
5 and a data output unit 16 as data output means.

【0020】この場合、探索領域設定部12は、探索領
域設定工程として、作業エリア2のワーク5に対する自
律移動ロボット1の概略位置情報に基いた密度の複数の
量子化空間点Eを有する空間探索領域Sを設定する。透
視変換行列演算部13は、透視変換行列演算工程とし
て、設定された空間探索領域Sの量子化空間点Eにおい
てカメラ4によりワーク5を撮影したときに画像信号か
ら得られる画像内における画像特徴点の座標を求めるた
めの透視変換行列Mを、自律移動ロボット1の概略位置
情報およびカメラ4の姿勢パラメータに応じて演算す
る。
In this case, as the search area setting step, the search area setting section 12 performs a space search having a plurality of quantization space points E having a density based on the approximate position information of the autonomous mobile robot 1 with respect to the work 5 in the work area 2. An area S is set. The perspective transformation matrix calculation unit 13 performs, as a perspective transformation matrix calculation step, image feature points in an image obtained from an image signal when the camera 4 photographs the work 5 at the quantization space point E in the set space search area S. Is calculated in accordance with the approximate position information of the autonomous mobile robot 1 and the posture parameter of the camera 4.

【0021】誤差演算部14は、誤差演算工程として、
最小二乗推定部10にて求められた推定透視変換行列M
sと透視変換行列演算部13にて演算された透視変換行
列Mとの間の誤差を各行列要素の二乗誤差和として求め
る誤差評価関数に基いて演算する。誤差判定部15は、
誤差判定工程として、演算された二乗誤差和が要求誤差
精度εよりも小さいか否かを判定し、小さいときにはそ
のときのカメラ4の位置および姿勢のデータを出力端子
17に出力し、大きいときには再び探索領域設定部12
に戻って上述の各部において演算を実行させる。そし
て、データ出力部16は、データ出力工程として、上述
の判定をした結果に基いてそのときの空間探索領域Sに
おける量子化空間点Eを出力する。
The error calculating section 14 performs an error calculating step as follows.
Estimated perspective transformation matrix M obtained by least-squares estimator 10
The calculation is performed based on an error evaluation function that calculates an error between s and the perspective transformation matrix M computed by the perspective transformation matrix computation unit 13 as a sum of square errors of each matrix element. The error determination unit 15
As an error determination step, it is determined whether or not the calculated sum of squared errors is smaller than the required error precision ε. If smaller, the data of the position and orientation of the camera 4 at that time is output to the output terminal 17. Search area setting unit 12
Then, the calculation is executed in each of the above-described units. Then, as a data output step, the data output unit 16 outputs the quantization space point E in the space search area S at that time based on the result of the above determination.

【0022】また、量子化空間探索部11には、後述す
る各種データが記憶される記憶手段としてのメモリ18
が接続され、その各種データの書き込みおよび読み出し
動作が行えるようになっていると共に、外部情報入力端
子19が設けられ、自律移動ロボット1の走行制御を行
う図示しない制御部あるいは図示しない位置検出センサ
などから自律移動ロボット1の概略位置を示す信号が入
力されるようになっている。この場合、メモリ18にあ
らかじめ記憶されるデータとしては、特徴点の対応リス
ト,ワーク5および作業エリア2の三次元絶対座標デー
タ,カメラ4のパラメータ,量子化空間点Eのリストと
その座標データなどである。
The quantization space search unit 11 has a memory 18 as storage means for storing various data described later.
Are connected, and write and read operations of the various data can be performed. An external information input terminal 19 is provided, and a control unit (not shown) for controlling the traveling of the autonomous mobile robot 1 or a position detection sensor (not shown) is provided. , A signal indicating the approximate position of the autonomous mobile robot 1 is input. In this case, the data stored in advance in the memory 18 include a correspondence list of feature points, three-dimensional absolute coordinate data of the work 5 and the work area 2, parameters of the camera 4, a list of quantization space points E and their coordinate data, and the like. It is.

【0023】次に、本実施例の作用について図3ないし
図10をも参照して説明する。三次元位置認識装置は、
その動作の流れとして図3に示すようなステップにした
がって自律移動ロボット1の存在位置を認識する。ま
ず、図2の状態において、自律移動ロボット1が軌道上
を移動してきて作業エリア2に差し掛かると、カメラ4
の撮影シーン内に作業エリア2が入って撮影されると、
画像処理部7は、カメラ4の画像信号を入力するように
なる(ステップT1)。この場合、自律移動ロボット1
は、所定軌道上を所定方向に移動しており、また移動を
制御する制御部の内部情報などにより概略の停止位置が
わかっているので、その停止予定位置において対象物体
としてのワーク5がカメラ4の視野内に入るように設定
することができる。
Next, the operation of this embodiment will be described with reference to FIGS. The three-dimensional position recognition device
The existence position of the autonomous mobile robot 1 is recognized in accordance with the steps shown in FIG. First, in the state of FIG. 2, when the autonomous mobile robot 1 moves on the track and approaches the work area 2, the camera 4
When the work area 2 enters the shooting scene and is shot,
The image processing section 7 receives the image signal of the camera 4 (step T1). In this case, the autonomous mobile robot 1
Is moving in a predetermined direction on a predetermined trajectory, and a rough stop position is known from internal information of a control unit for controlling the movement. Can be set so as to fall within the field of view.

【0024】次に、カメラ4により撮影されたワーク5
の画像信号は、画像処理部7の前処理計算部8にて前処
理としての雑音除去処理(ステップT2)およびワーク
5の画像の切り出し処理(ステップT3)を実行し、こ
の後、特徴点抽出部9にてその画像信号に基いて次のよ
うにしてm個の画像特徴点Qjの集合Qを求めて出力す
るようになる(ステップT4)。ここで、画像特徴点Q
jの集合Qは次のように定義される。 Q={Qj;j=1,2,…,m} …(1) Qj=(uj,vj) …(2) u,v;表示画面上における画素単位の二次元座標の各
要素値を示す
Next, a work 5 photographed by the camera 4
Of the image signal of the image processing unit 7, the preprocessing calculation unit 8 of the image processing unit 7 executes a noise removal process (step T2) as a preprocess and a cutout process of the image of the work 5 (step T3), and thereafter, the feature point extraction The unit 9 calculates and outputs a set Q of m image feature points Qj based on the image signal as follows (step T4). Here, the image feature point Q
The set Q of j is defined as follows. Q = {Qj; j = 1, 2,..., M} (1) Qj = (uj, vj) (2) u, v; each element value of two-dimensional coordinates in pixel units on the display screen

【0025】この画像特徴点Qjの抽出に際しては、外
部情報として入力端子19から与えられる概略位置情報
に基いて次のように演算を行う。すなわち、概略位置情
報として、例えば、自律移動ロボット1の移動経路を記
憶した作業計画,走行制御により移動した空間移動の履
歴情報などが入力されるので、これらの外部情報に基い
て、自律移動ロボット1の中心位置の概略的な三次元絶
対座標として、絶対座標で表した概略位置A0、 A0=(X0 ,Y0 ,Z0 ) …(3) を得ることができる。なお、概略位置A0の三次元絶対
座標は、作業エリア2の所定の位置を原点Opとして表
現したものである。
In extracting the image feature point Qj, the following calculation is performed based on the approximate position information provided from the input terminal 19 as external information. That is, as the approximate position information, for example, a work plan in which the moving path of the autonomous mobile robot 1 is stored, history information of the spatial movement moved by the travel control, and the like are input. As approximate three-dimensional absolute coordinates of the center position of No. 1, approximate positions A0, A0 = (X0, Y0, Z0) (3) expressed in absolute coordinates can be obtained. Note that the three-dimensional absolute coordinates of the approximate position A0 represent a predetermined position in the work area 2 as the origin Op.

【0026】次に、メモリ18にあらかじめ記憶されて
いる作業エリア2に載置されたワーク5の種々の見え方
パターン(図4(a)ないし(f)参照)の中から、概
略位置A0の値に対応する見え方パターンを選択する
(ステップT5)。いま、例えば、概略位置A0に対応
する見え方パターンが図4(b)に対応しているときに
は、その見え方パターンに現れているワーク5の三次元
の特徴点Piの集合、 P={Pi;i=1,2,…,n} …(4) に対応する画像特徴点Qjを求めて順次抽出する。抽出
した画像特徴点Qjはメモリ18に記憶される。これ
は、例えば、図9(a)に示すように、作業エリア2と
ワーク5との特徴点Piの座標が記述されたメモリ18
の座標データと、同図(b)に示す画像信号により得ら
れた画像特徴点Qjとの間の対応関係を上述の見え方パ
ターン(図4(b)参照)に対応させて選択し、これに
基いて特徴点の対応リストを選択するのである。
Next, from various appearance patterns (see FIGS. 4A to 4F) of the work 5 placed on the work area 2 stored in the memory 18 in advance, the approximate position A0 is determined. The appearance pattern corresponding to the value is selected (step T5). Now, for example, when the appearance pattern corresponding to the approximate position A0 corresponds to FIG. 4B, a set of three-dimensional feature points Pi of the work 5 appearing in the appearance pattern, P = {Pi I = 1, 2,..., N} (4) Image feature points Qj corresponding to The extracted image feature points Qj are stored in the memory 18. This is, for example, as shown in FIG. 9A, the memory 18 in which the coordinates of the feature points Pi of the work area 2 and the work 5 are described.
Is selected in correspondence with the above-mentioned appearance pattern (see FIG. 4 (b)), and the correspondence between the coordinate data of (1) and the image feature point Qj obtained by the image signal shown in FIG. The correspondence list of the feature points is selected based on the.

【0027】さて、最小二乗推定部10においては、上
述のように選択された対応リストに基いてカメラ4の現
在位置と姿勢に対するカメラ行列Cの最小二乗推定を行
う(ステップT6)。この場合、カメラ行列Cは、次式
(5),(6),(7)に基いて設定されるもので、
(u,v,w)は斉次座標系における画像表示面の1点
を表し、(U,V)は画像表示面の座標軸で表した座標
値であり、(x,y,z)は三次元絶対座標空間におい
て対応する1点の座標である。つまり、カメラ行列Cは
三次元絶対空間の座標と画像表示面上との間の変換を行
うためのパラメータである。なお、wはカメラ4による
スケールファクタを表している。
The least-squares estimating unit 10 performs least-squares estimation of the camera matrix C for the current position and orientation of the camera 4 based on the correspondence list selected as described above (step T6). In this case, the camera matrix C is set based on the following equations (5), (6), and (7).
(U, v, w) represents one point on the image display surface in the homogeneous coordinate system, (U, V) is a coordinate value represented by a coordinate axis of the image display surface, and (x, y, z) is a tertiary coordinate. These are the coordinates of one corresponding point in the original absolute coordinate space. That is, the camera matrix C is a parameter for performing conversion between the coordinates in the three-dimensional absolute space and the image display surface. Note that w represents a scale factor by the camera 4.

【0028】[0028]

【数1】 (Equation 1)

【0029】そして、式(5),(6)から、 u−Uw=0 , v−Vw=0 …(8) が導かれ、また、斉次座標の性質を考慮すると、式
(7)から次式(9),(10),(11)が得られ
る。
From the equations (5) and (6), u-Uw = 0, v-Vw = 0 (8) is derived. Further, considering the properties of the homogeneous coordinates, the following equation is obtained. The following equations (9), (10), and (11) are obtained.

【0030】[0030]

【数2】 (Equation 2)

【0031】次に、式(9),(10),(11)を式
(8)に代入すると、次式(12),(13)が得られ
る。
Next, when the equations (9), (10) and (11) are substituted into the equation (8), the following equations (12) and (13) are obtained.

【0032】[0032]

【数3】 これらの式から次式(14),(15)のようにスケー
ルファクタとしてのwを消去して変形することができ
る。
(Equation 3) From these equations, w as the scale factor can be deleted and transformed as in the following equations (14) and (15).

【0033】[0033]

【数4】 この結果から、次式(16),(17)で示される評価
関数を導出することができ、Eu,Evを最小化するよ
うな次式(18),(19)で示される条件のもとでカ
メラ行列Cの各要素を係数とする線形方程式を得て最終
的にカメラ行列Cを決定することができる。
(Equation 4) From this result, an evaluation function represented by the following equations (16) and (17) can be derived, and under the conditions represented by the following equations (18) and (19) that minimize Eu and Ev. Then, a linear equation using each element of the camera matrix C as a coefficient is obtained to finally determine the camera matrix C.

【0034】[0034]

【数5】 そして、このようにして得られたカメラ行列C(広義の
透視変換行列)の推定結果を推定透視変換行列Msとし
て記憶する。
(Equation 5) Then, the estimation result of the camera matrix C (perspective transformation matrix in a broad sense) obtained as described above is stored as the estimated perspective transformation matrix Ms.

【0035】次に、量子化空間探索部11においては、
探索領域設定部12により、概略位置A0(X0 ,Y0
,Z0 )を基準とした空間探索領域Sc、 Sc(c=1,2,…,C) …(20) を設定する(ステップT7)。この空間探索領域Sc
は、図5に示すように、所定の密度で格子状に量子化し
て配列された複数の量子化空間点E(k,c)、 E(k,c)(k=1,2,…,K、c=1,2,…,C) …(21) を含んで構成される。
Next, in the quantization space search unit 11,
The search area setting unit 12 sets the approximate position A0 (X0, Y0
, Z0) are set as the space search areas Sc, Sc (c = 1, 2,..., C) (20) (step T7). This space search area Sc
Are, as shown in FIG. 5, a plurality of quantized spatial points E (k, c), E (k, c) (k = 1, 2,... K, c = 1, 2,..., C) (21)

【0036】この場合、作業エリア2の原点Opを基準
とした絶対座標系におけるX軸方向(作業エリア2に対
向する方向)の座標値X=XL1の平面内で、カメラ4は
ワーク5を向いていると仮定できる。ここで、XL1はカ
メラ4の視点の位置を示すX座標値の推定値であり、こ
の値は自律移動ロボット1の概略位置A0のX方向の座
標値X0 とカメラ4のロボット上の相対座標から決定で
きる。また、本実施例における自律移動ロボット1の場
合には、概略位置A0の高さ方向の座標値Z0は、一定
の値から大きく外れることがないと想定されるので、カ
メラ4のZ方向の座標を求めるための探索範囲は実質的
に狭い範囲に絞り込むことができる。
In this case, the camera 4 faces the work 5 in the plane of the coordinate value X = XL1 in the X-axis direction (the direction facing the work area 2) in the absolute coordinate system based on the origin Op of the work area 2. Can be assumed. Here, XL1 is an estimated value of the X coordinate value indicating the position of the viewpoint of the camera 4, and this value is obtained from the coordinate value X0 of the approximate position A0 of the autonomous mobile robot 1 in the X direction and the relative coordinates of the camera 4 on the robot. Can decide. Further, in the case of the autonomous mobile robot 1 in the present embodiment, since the coordinate value Z0 in the height direction of the approximate position A0 is assumed not to largely deviate from a certain value, the coordinate value of the camera 4 in the Z direction is assumed. Can be narrowed down to a substantially narrow range.

【0037】次に、カメラ4の各種パラメータの探索範
囲を設定する(ステップT8)。ここで、カメラ4のパ
ラメータとしては、図6に対応関係を示すように、 (a)カメラ4の姿勢パラメータCP(α,β,γ)、 (b)カメラ4の焦点距離特性f (c)カメラ4の位置(r,θ,φ) がある。そして、この場合において、(b)の焦点距離
特性fはカメラ4の特性によりあらかじめ決まっている
ので、例えば、メモリ18内にデータとして記憶させて
おくことができる。また、(c)のカメラ4の位置
(r,θ,φ)の値は以降の演算工程にて決めることが
できる値である。
Next, a search range of various parameters of the camera 4 is set (step T8). Here, as the parameters of the camera 4, as shown in a correspondence relationship in FIG. 6, (a) a posture parameter CP (α, β, γ) of the camera 4, (b) a focal length characteristic f (c) of the camera 4 There is a position (r, θ, φ) of the camera 4. In this case, since the focal length characteristic f in (b) is determined in advance by the characteristic of the camera 4, it can be stored as data in the memory 18, for example. Further, the value of the position (r, θ, φ) of the camera 4 in (c) is a value that can be determined in the subsequent calculation steps.

【0038】したがって、ここでは、カメラ4のワーク
5の中心位置からの光軸のずれを表している姿勢パラメ
ータCP(α,β,γ)の探索範囲を決定すれば良い。
すなわち、カメラ4の姿勢パラメータCPの探索範囲
は、例えば図7に示すように、γ=0の場合に対して、
ワーク5の所定位置としての原点Opの画像位置Ops
(画像情報に含まれていない場合には含まれている特徴
点を原点とする座標系に変換する)が画像の表示画面上
における画素単位の二次元座標(u,v)において、例
えば(Na,Nb)であるとすると、次式よりαおよび
βの推定値が計算できる。
Therefore, here, the search range of the posture parameter CP (α, β, γ) representing the deviation of the optical axis from the center position of the work 5 of the camera 4 may be determined.
That is, the search range of the posture parameter CP of the camera 4 is, for example, as shown in FIG.
Image position Ops of origin Op as predetermined position of work 5
(If it is not included in the image information, it is converted into a coordinate system with the feature point included as the origin.) In pixel-based two-dimensional coordinates (u, v) on the image display screen, for example, (Na , Nb), the estimated values of α and β can be calculated from the following equations.

【0039】 α=2×arctan{(Na/NH )× tan(wH /2)} …(22) β=2×arctan{(Nb/NV )× tan(wV /2)} …(23) ただし、NH ;画像の表示画面上における水平方向の画
像サイズ NV ;画像の表示画面上における垂直方向の画像サイズ wH ;カメラの水平画角 wV ;カメラの垂直画角 である。
Α = 2 × arctan {(Na / NH) × tan (wH / 2)} (22) β = 2 × arctan {(Nb / NV) × tan (wV / 2)} (23) , NH; the horizontal image size of the image on the display screen NV; the vertical image size of the image on the display screen wH; the horizontal field angle of the camera wV; the vertical field angle of the camera.

【0040】この場合、要求誤差精度があまり高くない
場合には、線形近似として、 α=(Na/NH )×wH …(24) β=(Nb/NV )×wV …(25) という式で代用することもできる。
In this case, if the required error accuracy is not very high, then, as a linear approximation, α = (Na / NH) × wH (24) β = (Nb / NV) × wV (25) It can be substituted.

【0041】そして、γ=0でない場合には、ワーク5
の中心の画像位置Opsを含む円環状の領域Sa(Op
s)内にγの補正をした場合のOpscが存在すると考
えられるので、Sa(Ops)内の量子化空間点Ops
cについて式(22),(23)あるいは式(24),
(25)に基いてα,βを定め、 γ=γpsc(OpsとOpscとのなす角度) として、以降の探索工程を行うようにする。なお、本実
施例においては、自律移動ロボット1は、略水平な床の
上を移動することと、概略位置A0の座標値を用いて高
精度にカメラ4の回転角を調整できるので、カメラ4の
光軸ずれに対する姿勢パラメータ(α,β,γ)のう
ち、最も広い探索能力を要求されるのは水平方向の角度
ずれを表すαである。したがって、β,γについては微
小範囲の探索を行えば良い。
If γ = 0 is not satisfied, the work 5
Annular region Sa (Op including the image position Ops at the center of
Since it is considered that Opsc when γ is corrected exists in s), the quantization space point Ops in Sa (Ops) is considered.
Equation c (22), (23) or equation (24)
Based on (25), α and β are determined, and γ = γpsc (the angle between Ops and Opsc), and the subsequent search step is performed. In the present embodiment, the autonomous mobile robot 1 can move on a substantially horizontal floor and adjust the rotation angle of the camera 4 with high accuracy using the coordinate value of the approximate position A0. Among the attitude parameters (α, β, γ) for the optical axis shift, the one that requires the widest search capability is α representing the horizontal angle shift. Therefore, it is sufficient to search for a minute range for β and γ.

【0042】次に、透視変換行列演算部13にて、空間
探索領域S1内の初期探索点として量子化空間点E
(k,1)(k=1,2,…,K)のうちの1点である
量子化空間点E(1,1)を選択して(ステップT9)
その三次元絶対座標と、カメラ4の姿勢パラメータCP
(α,β,γ)の初期値とを用いて、 透視変換行列M(k,c)(k=1,2,…,K、c=1,2,…,C) …(26) を演算する(ステップT10)。
Next, the perspective transformation matrix operation unit 13 sets the quantized spatial points E as initial search points in the spatial search area S1.
A quantization space point E (1,1) which is one point of (k, 1) (k = 1, 2,..., K) is selected (step T9).
The three-dimensional absolute coordinates and the posture parameter CP of the camera 4
Using the initial values of (α, β, γ), the perspective transformation matrix M (k, c) (k = 1, 2,..., K, c = 1, 2,..., C) (26) An operation is performed (step T10).

【0043】この場合、透視変換行列Mは、次のように
して演算される。すなわち、まず、ワーク5の原点Op
に対して設定されている三次元絶対座標(X,Y,Z)
を基準として、初期探索領域S1を概略位置A0のX座
標値であるX=X0 の平面内にとり、この初期探索領域
S1内の三次元絶対座標(X,Y,Z)を三次元極座標
(r,θ,φ)に変換する透視変換行列Mを求める。カ
メラ4の姿勢パラメータCP(α,β,γ)と、焦点距
離特性fから求めたカメラ4の視点位置と射影平面に対
する距離dとに基いて、図6に示す関係に基いて透視変
換行列M(k,c)を求める。
In this case, the perspective transformation matrix M is calculated as follows. That is, first, the origin Op of the work 5
Three-dimensional absolute coordinates (X, Y, Z) set for
, The initial search area S1 is set within the plane of X = X0, which is the X coordinate value of the approximate position A0, and the three-dimensional absolute coordinates (X, Y, Z) in the initial search area S1 are converted to the three-dimensional polar coordinates (r , Θ, φ). Based on the posture parameter CP (α, β, γ) of the camera 4 and the viewpoint position of the camera 4 obtained from the focal length characteristic f and the distance d to the projection plane, the perspective transformation matrix M is obtained based on the relationship shown in FIG. (K, c) is obtained.

【0044】ここで、透視変換行列M(k,c)は次の
ように表すことができる。 M(k,c)=T・M1・M2・M3 …(27) ただし、 T ;回転,並行移動,拡大,縮小等の変換行列(ここ
では、カメラ4の姿勢パラメータCPで表される回転の
演算) M1;撮影しているワークの射影空間内への透視変換 M2;射影空間内の平面から射影平面への透視変換 M3;射影平面から画像の表示面への座標変換 このうち、画像の表示面とワーク5との間の透視変換を
演算するための透視変換行列M1とM2との積および座
標変換行列M3は、それぞれ、次式(28),(29)
のように表すことができる。
Here, the perspective transformation matrix M (k, c) can be expressed as follows. M (k, c) = T ・ M1 ・ M2 ・ M3 (27) where T: a transformation matrix for rotation, parallel movement, enlargement, reduction, etc. (here, the rotation matrix represented by the posture parameter CP of the camera 4) Calculation) M1: Perspective transformation of the work being photographed into the projection space M2: Perspective transformation from the plane in the projection space to the projection plane M3: Coordinate transformation from the projection plane to the display surface of the image The product of the perspective transformation matrices M1 and M2 for calculating the perspective transformation between the surface and the work 5 and the coordinate transformation matrix M3 are given by the following equations (28) and (29), respectively.
Can be expressed as

【0045】[0045]

【数6】 また、カメラ4の光軸ずれに対応する姿勢パラメータC
P(α,β,γ)に応じた回転の演算を行うための変換
行列Tは、前述の図6と図8に示す関係とに基いて次式
(30),(31),(32)のように表すことができ
る。なお、このとき、光軸ずれの変換においては、、ワ
ーク5との距離が変化しないと仮定すると、T0 =T1
とおくことができるので、次式(33)の関係が成立す
る。
(Equation 6) Also, the posture parameter C corresponding to the optical axis shift of the camera 4
The transformation matrix T for performing the rotation operation in accordance with P (α, β, γ) is given by the following equations (30), (31), and (32) based on the relationship shown in FIGS. Can be expressed as At this time, in the conversion of the optical axis shift, assuming that the distance from the workpiece 5 does not change, T0 = T1.
Therefore, the following equation (33) holds.

【0046】[0046]

【数7】 次に、誤差演算部14において、このようにして演算に
より求めた透視変換行列M(k,r)と、最小二乗推定
部10で求めた推定透視変換行列Msとの間の誤差評価
関数D(1,1)を計算する(ステップT11)。この
場合、誤差評価関数D(k,c)は透視変換行列Mと推
定透視変換行列Msとの各行列要素間の二乗誤差の和を
演算するものである。
(Equation 7) Next, the error calculation unit 14 calculates an error evaluation function D () between the perspective transformation matrix M (k, r) obtained by the calculation in this way and the estimated perspective transformation matrix Ms obtained by the least square estimation unit 10. (1, 1) is calculated (step T11). In this case, the error evaluation function D (k, c) calculates the sum of the square errors between the matrix elements of the perspective transformation matrix M and the estimated perspective transformation matrix Ms.

【0047】続いて、現在の量子化空間点E(1,1)
において姿勢パラメータCP(α,β,γ)のうちγを
そのままにした状態でα,βを変更し、上述と同様にし
て誤差評価関数D(nE ,1)(nE =1,2,…,N
E )を演算する(ステップT12,13,10,1
1)。この後、これらの誤差評価関数D(nE ,1)か
ら最適値としてDE1(αE ,βE ,γE )を算出する
(ステップT14)。続いて、誤差判定部15におい
て、演算された誤差評価関数D(k,1)の値から得ら
れた最適値DE1の演算値が、あらかじめ設定されている
検出条件としての要求誤差精度εの値に対して、 DE k<ε …(34) とした判定条件を満たすか否かを判断する(ステップT
15)。
Subsequently, the current quantization space point E (1, 1)
In the posture parameters CP (α, β, γ), α and β are changed while γ remains unchanged, and the error evaluation function D (nE, 1) (nE = 1, 2,. N
E) (Steps T12, 13, 10, 1)
1). Thereafter, DE1 (αE, βE, γE) is calculated as an optimum value from these error evaluation functions D (nE, 1) (step T14). Subsequently, in the error determination unit 15, the calculated value of the optimum value DE1 obtained from the value of the calculated error evaluation function D (k, 1) is used as the value of the required error accuracy ε as a preset detection condition. It is determined whether or not the following determination condition is satisfied: DE k <ε (34)
15).

【0048】そして、上述の判定条件式(37)を満た
していない場合には、以上の演算を初期探索領域S1内
における他の量子化空間点E(k,1)についても同様
に実行する(ステップT16,T17を経てステップT
9〜T15を実行)。さらに、他の量子化空間点E
(k,1)について判定条件式(34)を満たさない場
合には、そのときの最適値DE kのなかから要求誤差精
度εに一番近い値が得られた量子化空間点E(k,1)
の位置に基づいて、新たな空間探索領域S2を設定し
(ステップT18,T19)、前述の演算工程を再び実
行する(ステップT5〜T15)。
If the above-mentioned determination condition expression (37) is not satisfied, the above operation is similarly performed for the other quantization space points E (k, 1) in the initial search area S1 ( After steps T16 and T17, step T
9 to T15). Further, other quantization space points E
If the determination condition expression (34) is not satisfied for (k, 1), the quantization space point E (k, 1) at which a value closest to the required error accuracy ε is obtained from the optimum value DE k at that time. 1)
A new space search area S2 is set based on the position (steps T18 and T19), and the above-described calculation process is executed again (steps T5 to T15).

【0049】一方、判定条件式(34)の条件を満たし
た場合には(ステップT15)、データ出力部16によ
り、そのときの最適値DE kに相当する量子化空間点E
(k,c)の位置(rE ,θE ,φE )の値およびカメ
ラ4の姿勢パラメータ(αE,βE ,γE )の値を検出
データとして出力端子17に出力し(ステップT2
0)、探索プログラムを終了する。
On the other hand, when the condition of the criterion condition (34) is satisfied (step T15), the data output unit 16 outputs the quantization space point E corresponding to the optimum value DE k at that time.
The value of the position (rE, θE, φE) of (k, c) and the value of the posture parameter (αE, βE, γE) of the camera 4 are output to the output terminal 17 as detection data (step T2).
0), end the search program.

【0050】また、最終の空間探索領域SCの設定にも
拘らず判定条件式(34)を満たさない場合には(ステ
ップT18)、要求誤差精度εを満たす量子化空間点E
(k,c)が存在しないことをメッセージ出力した(ス
テップT21)後、それまでの間に演算された最適値D
E kのなかから要求誤差精度εに一番近い値が得られた
量子化空間点E(k,c)の位置(rE ,θE ,φE )
の値およびカメラ4の姿勢パラメータ(αE ,βE ,γ
E )の値を検出データとして出力端子17に出力して
(ステップT20)探索プログラムを終了する。
If the determination condition expression (34) is not satisfied despite the setting of the final space search area SC (step T18), the quantization space point E satisfying the required error accuracy ε is obtained.
After outputting a message that (k, c) does not exist (step T21), the optimum value D calculated up to that time is output.
The position (rE, θE, φE) of the quantization space point E (k, c) at which the value closest to the required error precision ε is obtained from Ek
And the posture parameters of the camera 4 (αE, βE, γ
The value of E) is output to the output terminal 17 as detection data (step T20), and the search program ends.

【0051】なお、上述の場合には、作業エリア2のワ
ーク5の原点Opから一定の距離(X=XL1)におい
て、要求誤差精度εを満たす量子化空間点E(k,1)
を求めるように、空間探索領域Sの設定をS1,S2,
…と設定する探索過程を示しているが、例えば、図10
に示すように、自律移動ロボット1の進行によるワーク
5への近接(X<XL1)するのに応じてその位置に応じ
た空間探索領域S(k,c)を設定する探索過程を実行
することもできる。
In the above case, at a certain distance (X = XL1) from the origin Op of the work 5 in the work area 2, the quantization space point E (k, 1) satisfying the required error accuracy ε is obtained.
Is set to S1, S2,
.. Are shown. For example, FIG.
Executing a search process of setting a spatial search area S (k, c) according to the position of the autonomous mobile robot 1 as it approaches (X <XL1) to the work 5 as the autonomous mobile robot 1 advances. Can also.

【0052】すなわち、例えば、自律移動ロボット1の
進行過程において、原点Opからの距離Xに対応して一
つの空間探索領域S1を設定し、そのなかの量子化空間
点E(k,1)のうち要求誤差精度εに一番近いものを
カメラ4の位置および姿勢パラメータを示す仮の検出デ
ータとして求め、その検出データに基いて自律移動ロボ
ット1の移動量を制御し、原点Opからの距離Xがさら
に小さくなったときに新たに空間探索領域S2を設定し
てさらに要求誤差精度εに一番近い量子化空間点E
(k,2)を求めるようにし、以下、同様にして、ワー
ク5に自律移動ロボット1が最も近接したときに、最終
的に要求誤差精度εを満たすように探索過程を設定する
ようにしても良い。
That is, for example, in the course of the movement of the autonomous mobile robot 1, one space search area S1 is set corresponding to the distance X from the origin Op, and the quantization space point E (k, 1) in the space search area S1 is set. Among them, the one closest to the required error accuracy ε is obtained as temporary detection data indicating the position and orientation parameters of the camera 4, the amount of movement of the autonomous mobile robot 1 is controlled based on the detection data, and the distance X from the origin Op is calculated. Becomes smaller, a new space search area S2 is set, and the quantization spatial point E closest to the required error accuracy ε is further set.
(K, 2) is calculated, and similarly, when the autonomous mobile robot 1 is closest to the workpiece 5, the search process may be set so as to finally satisfy the required error accuracy ε. good.

【0053】このように、本実施例においては、既知形
状の対象物体としての作業エリア2およびワーク5をカ
メラ4により撮像してその画像情報から画像特徴点Qj
を抽出し、最小二乗推定部10により、抽出した画像特
徴点とワーク5の見え方のパターンとの間の特徴点の対
応リストを選択してカメラ4の現在位置と姿勢とにおけ
るワーク5の見え方を表すカメラ行列の最小二乗推定を
行って推定透視変換行列Msを求め、一方、量子化した
空間点を複数個有する空間探索領域Scを設定すると共
に、カメラ4の姿勢パラメータCP(α,β,γ)を設
定してワーク5をその量子化空間点E(k,c)におい
て撮像したときに得られるべきワーク5の画像特徴点を
求めるための透視変換行列M(k,c)を演算し、それ
らの対応関係を求めて誤差評価関数D(k,c)により
要求誤差精度εを満たすか否かを判定してマッチングを
とることによりカメラ4の位置および姿勢パラメータC
Pを推定するようにしたので、以下のような効果を得る
ことができる。
As described above, in the present embodiment, the work area 2 and the work 5 as a target object having a known shape are imaged by the camera 4 and the image feature points Qj are obtained from the image information.
Is extracted by the least-squares estimating unit 10, and a correspondence list of feature points between the extracted image feature points and the appearance pattern of the work 5 is selected, and the appearance of the work 5 at the current position and orientation of the camera 4 is selected. The estimated perspective transformation matrix Ms is obtained by performing least-squares estimation of the camera matrix representing the two directions, while a spatial search area Sc having a plurality of quantized spatial points is set, and the posture parameters CP (α, β , Γ), and calculates a perspective transformation matrix M (k, c) for obtaining an image feature point of the work 5 to be obtained when the work 5 is imaged at the quantization space point E (k, c). The position and orientation parameters C of the camera 4 are obtained by determining whether or not the required error accuracy ε is satisfied by the error evaluation function D (k, c) by obtaining the correspondence between them and performing matching.
Since P is estimated, the following effects can be obtained.

【0054】すなわち、第1に、従来と異なり、特殊マ
ーカを設ける必要がなくなり、マーカ認識用の移動時間
が不要となって演算時間を短縮して高速化を図れる。第
2に、透視変換行列の成分からパラメータを算出する方
法ではなく、所定密度で配置される量子化空間点Eを有
する空間探索領域Scを設定して透視変換行列を求め、
推定透視変換行列との間の適合度を評価するので、演算
量が膨大となることなく迅速に演算処理を実行できる。
第3に、カメラ4の姿勢パラメータCP(α,β,γ)
が未知の場合でも適用できる。第4に、カメラ4の位置
を求めるための演算の精度を、量子化空間点Eの取り方
で調整することができる。つまり、例えば、作業内容や
環境の実情を考慮して遠距離においては粗い精度で高速
演算を行わせ、そしてワーク5に近接したときには細か
い精度で正確な演算を行なわせるように設定するという
ように、情況に応じた適切な精度に設定することができ
る。第5に、このように演算を実行するので、対象物体
としてのワーク5に対して、能動的にカメラ4が移動す
る際に、遠距離での光軸合わせ(粗い推定)から近距離
での高精度推定まで、移動系と連動して多段的に用いる
ことができる。
First, unlike the related art, there is no need to provide a special marker, and no moving time for marker recognition is required, so that the calculation time can be reduced and the speed can be increased. Second, instead of calculating parameters from the components of the perspective transformation matrix, a spatial search area Sc having quantized spatial points E arranged at a predetermined density is set to obtain a perspective transformation matrix,
Since the degree of conformity with the estimated perspective transformation matrix is evaluated, arithmetic processing can be executed quickly without an enormous amount of arithmetic.
Third, the posture parameter CP (α, β, γ) of the camera 4
Is applicable even if is unknown. Fourth, the accuracy of the calculation for obtaining the position of the camera 4 can be adjusted by how to obtain the quantization space point E. That is, for example, in consideration of the work content and the actual circumstances of the environment, it is set so that high-speed calculation is performed with coarse accuracy at long distances, and accurate calculation is performed with fine accuracy when approaching the work 5. It can be set to an appropriate accuracy according to the situation. Fifth, since the calculation is performed in this manner, when the camera 4 is actively moving with respect to the workpiece 5 as the target object, the optical axis alignment at a long distance (coarse estimation) can be used for a short distance. Until highly accurate estimation, it can be used in multiple stages in conjunction with the mobile system.

【0055】[0055]

【発明の効果】以上説明したように、本発明の視覚に基
く三次元位置および姿勢の認識方法によれば、画像処理
工程により、形状を示す三次元絶対座標が既知である対
象物体を撮像手段により撮影したときの画像情報に基い
てその対象物体の形状の特徴を示す二次元の画像特徴点
を抽出し、最小二乗推定工程において、撮像手段の概略
位置と姿勢とに基いて、抽出した画像特徴点と記憶手段
に記憶された対象物体の見え方パターンとの間の対応リ
ストを選択し、カメラ行列の最小二乗推定により推定透
視変換行列を求め、探索領域設定工程,透視変換行列演
算工程,誤差演算工程,誤差判定工程およびデータ出力
工程により、撮像手段の概略的な推定位置に応じた探索
領域を設定してその探索領域における複数の量子化され
た空間点における対象物体に対する透視変換行列を演算
して前記推定透視変換行列との誤差を誤差評価関数によ
り評価して要求誤差精度を満たすときにそのときの量子
化空間点を撮像手段の三次元絶対座標位置として推定す
ると共に、撮像手段の姿勢パラメータを姿勢として出力
するようにしたので、撮像手段の光軸ずれがある場合で
も、特殊マーカなどを用いることなく、迅速且つ高精度
に撮像手段の位置と姿勢を認識することができ、しか
も、要求誤差精度の設定の仕方により演算速度と演算精
度との兼ね合いを距離に応じて適切に設定することがで
き、これらのデータに基いて対象物体に対する装置の移
動制御などを適確に行うことができるという優れた効果
を奏するものである。
As described above, according to the method for recognizing a three-dimensional position and orientation based on vision according to the present invention, the image processing means captures an object whose three-dimensional absolute coordinates indicating its shape are known by the image processing step. The two-dimensional image feature points indicating the characteristics of the shape of the target object are extracted based on the image information obtained when the image is captured, and in the least square estimation step, the extracted image is extracted based on the approximate position and orientation of the imaging means. A correspondence list between the feature points and the appearance pattern of the target object stored in the storage means is selected, an estimated perspective transformation matrix is obtained by least-squares estimation of the camera matrix, a search area setting step, a perspective transformation matrix calculation step, An error calculation step, an error determination step, and a data output step set a search area corresponding to the approximate estimated position of the imaging means, and set a plurality of quantized spatial points in the search area. When the perspective transformation matrix for the elephant object is calculated and the error with the estimated perspective transformation matrix is evaluated by an error evaluation function and the required error accuracy is satisfied, the quantization space point at that time is used as the three-dimensional absolute coordinate position of the imaging means. In addition to the estimation, the posture parameter of the imaging unit is output as the posture. Therefore, even if there is an optical axis deviation of the imaging unit, the position and posture of the imaging unit can be quickly and accurately determined without using a special marker. Recognition can be performed, and the balance between the calculation speed and the calculation accuracy can be appropriately set according to the distance depending on the method of setting the required error accuracy. Based on these data, the movement control of the apparatus with respect to the target object can be performed. And the like can be accurately performed.

【0056】また、本発明の視覚に基く三次元位置およ
び姿勢の認識装置によれば、上述の各工程を実行する手
段により、同様の動作を行なわせることができ、実際の
撮像手段の三次元絶対座標に最も近い量子化空間点を近
似点として求めることができると共に、撮像手段の光軸
のずれ等を表す姿勢パラメータを検出することができる
ようになり、これらのデータに基いて対象物体に対する
装置の移動制御などを的確に行うことができるようにな
るという優れた効果を奏する。
Further, according to the apparatus for recognizing three-dimensional position and orientation based on vision of the present invention, the same operation can be performed by the means for executing the above-described steps, and the three-dimensional operation of the actual imaging means can be performed. The quantization space point closest to the absolute coordinates can be obtained as an approximate point, and the posture parameter representing the shift of the optical axis of the imaging means can be detected. There is an excellent effect that the movement control of the device can be accurately performed.

【図面の簡単な説明】[Brief description of the drawings]

【図1】本発明の一実施例を示す機能的構成のブロック
FIG. 1 is a block diagram of a functional configuration showing an embodiment of the present invention.

【図2】自律移動ロボットと作業エリアとの位置関係を
示す図
FIG. 2 is a diagram showing a positional relationship between an autonomous mobile robot and a work area.

【図3】プロセスの流れを説明するフローチャートFIG. 3 is a flowchart illustrating a process flow.

【図4】種々の見え方パターンFIG. 4 Various appearance patterns

【図5】空間探索領域と作業エリアとの位置関係を示す
作用説明図
FIG. 5 is an operation explanatory view showing a positional relationship between a space search area and a work area.

【図6】絶対座標軸とカメラおよび画像の座標軸との関
係を示す図
FIG. 6 is a diagram showing a relationship between absolute coordinate axes and coordinate axes of a camera and an image.

【図7】画像表示部上のワークとカメラの姿勢パラメー
タとの関係を示す説明図
FIG. 7 is an explanatory diagram showing a relationship between a workpiece on an image display unit and a posture parameter of a camera.

【図8】透視変換行列を演算するための座標軸の関係を
示す作用説明図
FIG. 8 is an operation explanatory diagram showing a relationship between coordinate axes for calculating a perspective transformation matrix.

【図9】(a)作業エリアとワークの特徴点の定義と
(b)画像信号による画像特徴点との対応関係を示す作
用説明図
FIGS. 9A and 9B are operation explanatory diagrams showing a correspondence relationship between (a) a definition of a feature point of a work area and a work and (b) an image feature point based on an image signal.

【図10】探索領域をワークとの距離に応じて設定する
場合の作用説明図
FIG. 10 is an operation explanatory diagram when a search area is set according to a distance from a workpiece;

【符号の説明】[Explanation of symbols]

1は自律移動ロボット、2は作業エリア、4はカメラ
(撮像手段)、5はワーク(対象物体)、7は画像処理
部、8は前処理計算部、9は特徴点抽出部、10は最小
二乗推定部(最小二乗推定手段)、11は量子化空間探
索部、12は探索領域設定部(探索領域設定手段)、1
3は透視変換行列演算部(透視変換行列演算手段)、1
4は誤差演算部(誤差演算手段)、15は誤差判定部
(誤差判定手段)、16はデータ出力部(データ出力手
段)、18はメモリ(記憶手段)、19は外部情報入力
端子である。
1 is an autonomous mobile robot, 2 is a work area, 4 is a camera (imaging means), 5 is a work (target object), 7 is an image processing unit, 8 is a pre-processing calculation unit, 9 is a feature point extraction unit, and 10 is a minimum. A square estimating unit (least square estimating means), 11 is a quantization space searching unit, 12 is a search area setting unit (search area setting means), 1
3 is a perspective transformation matrix computing unit (perspective transformation matrix computing means), 1
Reference numeral 4 denotes an error calculation unit (error calculation unit), 15 denotes an error determination unit (error determination unit), 16 denotes a data output unit (data output unit), 18 denotes a memory (storage unit), and 19 denotes an external information input terminal.

フロントページの続き (58)調査した分野(Int.Cl.7,DB名) G01B 11/00 - 11/30 102 B25J 9/10 - 9/22 B25J 13/00 - 13/08 G06T 1/00 Continuation of the front page (58) Field surveyed (Int.Cl. 7 , DB name) G01B 11/00-11/30 102 B25J 9/10-9/22 B25J 13/00-13/08 G06T 1/00

Claims (2)

(57)【特許請求の範囲】(57) [Claims] 【請求項1】 形状を示す三次元絶対座標が既知である
対象物体を撮像手段により撮影したときの画像情報に基
いてその対象物体の形状の特徴を示す二次元の画像特徴
点を抽出する画像処理工程と、 前記撮像手段の概略位置と姿勢とに基いて、前記画像処
理工程で抽出した画像特徴点と記憶手段にあらかじめ記
憶された前記対象物体の見え方のパターンとの間の特徴
点の対応リストを選択すると共に、その対応リストに基
いて前記撮像手段の現在位置と姿勢とにおける前記対象
物体の見え方を表すカメラ行列の最小二乗推定を行って
推定透視変換行列を求める最小二乗推定工程と、 前記撮像手段と前記対象物体との概略的な推定距離に応
じて必要な演算精度に対応する密度の量子化空間点を有
する探索領域を設定する探索領域設定工程と、 前記探索領域の量子化空間点において前記撮像手段によ
り前記対象物体を撮影したときに得られるべき二次元の
画像特徴点を算出するための透視変換行列を、あらかじ
め記憶手段に記憶された前記対象物体の三次元絶対座標
に対する前記量子化空間点の推定位置と前記撮像手段の
撮像姿勢を示す姿勢パラメータとに基いて演算する透視
変換行列演算工程と、 この透視変換行列演算工程にて求められた前記透視変換
行列と前記最小二乗推定工程において求められた推定透
視変換行列との間の適合度を示す評価誤差値を所定の誤
差評価関数に基いて演算する誤差演算工程と、 この誤差演算工程にて求められた前記評価誤差値を前記
探索領域に対応して設定された要求誤差精度と比較し、
前記評価誤差値が前記要求誤差精度よりも小さくなるこ
とを判定条件として、前記探索領域設定工程,前記透視
変換行列演算工程および前記誤差演算工程を繰り返し実
行させる誤差判定工程と、 この誤差判定工程において前記判定条件を満たしたとき
に、そのときの量子化空間点の三次元絶対座標と前記撮
像手段の姿勢パラメータとを出力するデータ出力工程と
からなる視覚に基く三次元位置および姿勢の認識方法。
1. An image for extracting a two-dimensional image feature point indicating a feature of a shape of a target object based on image information obtained by photographing the target object having a known three-dimensional absolute coordinate indicating a shape by an imaging unit. Processing step, based on the approximate position and orientation of the imaging means, the feature points between the image feature points extracted in the image processing step and the pattern of the appearance of the target object stored in advance in the storage means A least-squares estimation step of selecting a correspondence list and performing a least-squares estimation of a camera matrix representing the appearance of the target object at the current position and orientation of the imaging means based on the correspondence list to obtain an estimated perspective transformation matrix A search area setting step of setting a search area having a quantized spatial point having a density corresponding to a required calculation accuracy according to a rough estimated distance between the imaging unit and the target object; A perspective transformation matrix for calculating a two-dimensional image feature point to be obtained when the target object is photographed by the imaging unit at the quantization space point of the search area is stored in the storage unit in advance. A perspective transformation matrix calculation step of calculating based on the estimated position of the quantization space point with respect to the three-dimensional absolute coordinates and an attitude parameter indicating the imaging attitude of the imaging means; An error calculation step of calculating an evaluation error value indicating the degree of conformity between the perspective transformation matrix and the estimated perspective transformation matrix obtained in the least square estimation step based on a predetermined error evaluation function; Compare the obtained evaluation error value with the required error accuracy set corresponding to the search area,
An error determining step of repeatedly executing the search area setting step, the perspective transformation matrix calculating step, and the error calculating step, on the condition that the evaluation error value is smaller than the required error accuracy, A method for recognizing a three-dimensional position and orientation based on vision, comprising a data output step of outputting the three-dimensional absolute coordinates of the quantized space point at that time and the orientation parameter of the imaging means when the determination condition is satisfied.
【請求項2】 画像情報を出力する撮像手段と、 形状を示す三次元絶対座標が既知である対象物体を前記
撮像手段により撮影したときの画像情報に基いてその対
象物体の形状の特徴を示す二次元の画像特徴点を抽出す
る画像処理手段と、 前記撮像手段の概略位置と姿勢とに基いて、前記画像処
理手段により抽出した画像特徴点と記憶手段にあらかじ
め記憶された前記対象物体の見え方のパターンとの間の
特徴点の対応リストを選択すると共に、その対応リスト
に基いて前記撮像手段の現在位置と姿勢とにおける前記
対象物体の見え方を表すカメラ行列の最小二乗推定を行
って推定透視変換行列を求める最小二乗推定手段と、 前記撮像手段と前記対象物体との概略的な推定距離に応
じて必要な演算精度に対応する密度の量子化空間点を有
する探索領域を設定する探索領域設定手段と、 前記探索領域の量子化空間点において前記撮像手段によ
り前記対象物体を撮影したときに得られるべき二次元の
画像特徴点を算出するための透視変換行列を、あらかじ
め記憶手段に記憶された前記対象物体の三次元絶対座標
に対する前記量子化空間点の推定位置と前記撮像手段の
撮像姿勢を示す姿勢パラメータとに基いて演算する透視
変換行列演算手段と、 この透視変換行列演算手段により求められた前記透視変
換行列と前記最小二乗推定手段により求められた推定透
視変換行列との間の適合度を示す評価誤差値を所定の誤
差評価関数に基いて演算する誤差演算手段と、 この誤差演算手段にて求められた前記評価誤差値を前記
探索領域に対応して設定された要求誤差精度と比較し、
前記評価誤差値が前記要求誤差精度よりも小さくなるこ
とを判定条件として、前記探索領域設定手段,前記透視
変換行列演算手段および前記誤差演算手段に前記演算を
繰り返し実行させる誤差判定手段と、 この誤差判定手段により前記判定条件が満たされたとき
に、そのときの量子化空間点の三次元絶対座標と前記撮
像手段の姿勢パラメータとを出力するデータ出力手段と
からなる視覚に基く三次元位置および姿勢の認識装置。
2. An image pickup means for outputting image information, and a characteristic of the shape of the target object based on image information obtained by photographing the target object whose three-dimensional absolute coordinates indicating the shape is known by the image pickup means. Image processing means for extracting two-dimensional image feature points; image feature points extracted by the image processing means based on the approximate position and orientation of the imaging means; and the appearance of the target object stored in advance in storage means. And selecting a correspondence list of feature points between the two patterns and performing a least-squares estimation of a camera matrix representing the appearance of the target object at the current position and orientation of the imaging unit based on the correspondence list. A least-squares estimating unit for obtaining an estimated perspective transformation matrix; and a search having a quantized spatial point having a density corresponding to a required calculation accuracy according to a rough estimated distance between the imaging unit and the target object. Search area setting means for setting a search area; and a perspective transformation matrix for calculating a two-dimensional image feature point to be obtained when the target object is photographed by the imaging means at a quantization space point of the search area. A perspective transformation matrix calculating unit that calculates based on an estimated position of the quantization space point with respect to the three-dimensional absolute coordinates of the target object stored in advance in the storage unit and a posture parameter indicating an imaging posture of the imaging unit; An error for calculating an evaluation error value indicating a degree of conformity between the perspective transformation matrix calculated by the perspective transformation matrix calculation unit and the estimated perspective transformation matrix calculated by the least square estimation unit based on a predetermined error evaluation function. Calculating means, comparing the evaluation error value obtained by the error calculating means with a required error accuracy set corresponding to the search area,
An error determining unit that causes the search area setting unit, the perspective transformation matrix calculating unit, and the error calculating unit to repeatedly execute the calculation, based on a determination condition that the evaluation error value is smaller than the required error accuracy; When the determination condition is satisfied by the determination means, the three-dimensional absolute coordinate of the quantization space point at that time and the data output means for outputting the attitude parameter of the imaging means. Recognition device.
JP04924093A 1993-03-08 1993-03-10 Method and apparatus for recognizing three-dimensional position and orientation based on vision Expired - Lifetime JP3208900B2 (en)

Priority Applications (2)

Application Number Priority Date Filing Date Title
JP04924093A JP3208900B2 (en) 1993-03-10 1993-03-10 Method and apparatus for recognizing three-dimensional position and orientation based on vision
US08/206,828 US5499306A (en) 1993-03-08 1994-03-07 Position-and-attitude recognition method and apparatus by use of image pickup means

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP04924093A JP3208900B2 (en) 1993-03-10 1993-03-10 Method and apparatus for recognizing three-dimensional position and orientation based on vision

Publications (2)

Publication Number Publication Date
JPH06258028A JPH06258028A (en) 1994-09-16
JP3208900B2 true JP3208900B2 (en) 2001-09-17

Family

ID=12825353

Family Applications (1)

Application Number Title Priority Date Filing Date
JP04924093A Expired - Lifetime JP3208900B2 (en) 1993-03-08 1993-03-10 Method and apparatus for recognizing three-dimensional position and orientation based on vision

Country Status (1)

Country Link
JP (1) JP3208900B2 (en)

Families Citing this family (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP4595212B2 (en) * 2001-02-20 2010-12-08 コニカミノルタビジネステクノロジーズ株式会社 POSITION INFORMATION SETTING DEVICE AND ENVIRONMENTAL INFORMATION ACQUIRING DEVICE
KR101040957B1 (en) * 2009-07-27 2011-06-17 전자부품연구원 Location recognition method and robot
JP5494427B2 (en) * 2010-11-12 2014-05-14 富士通株式会社 Image processing program and image processing apparatus
JP6227950B2 (en) * 2013-09-19 2017-11-08 株式会社Ihi Sensor orientation control method and apparatus
CN106553195B (en) * 2016-11-25 2018-11-27 中国科学技术大学 Object 6DOF localization method and system during industrial robot crawl
JP7076097B2 (en) * 2018-04-24 2022-05-27 大成建設株式会社 Image leveling device and its program, and drawing generation system
CN114440776B (en) * 2022-01-28 2024-07-19 上海交途科技有限公司 Automatic displacement measurement method and system based on machine vision
CN118332708B (en) * 2024-06-17 2024-10-18 成都飞机工业(集团)有限责任公司 Large-size component assemblability evaluation method based on characteristic point distance error matrix

Also Published As

Publication number Publication date
JPH06258028A (en) 1994-09-16

Similar Documents

Publication Publication Date Title
US7283661B2 (en) Image processing apparatus
JP3525896B2 (en) Three-dimensional object recognition method and bin picking system using the same
CN110116407A (en) Flexible robot&#39;s pose measuring method and device
JP6922348B2 (en) Information processing equipment, methods, and programs
CN114022560A (en) Calibration method and related device and equipment
EP0157299A1 (en) Image processing apparatus
JP3208900B2 (en) Method and apparatus for recognizing three-dimensional position and orientation based on vision
JP2730457B2 (en) Three-dimensional position and posture recognition method based on vision and three-dimensional position and posture recognition device based on vision
JP3326851B2 (en) Method and apparatus for recognizing three-dimensional position and orientation based on vision
CN114742883B (en) An automated assembly method and system based on planar workpiece positioning algorithm
US20090297036A1 (en) Object detection on a pixel plane in a digital image sequence
Das et al. Sensor fusion in autonomous vehicle using LiDAR and camera Sensor
JP2003136465A (en) Method for determining 3D position / posture of detection target and visual sensor for robot
JP2001076128A (en) Obstacle detection apparatus and method
JP2778430B2 (en) Three-dimensional position and posture recognition method based on vision and three-dimensional position and posture recognition device based on vision
Viéville et al. Experimenting with 3D vision on a robotic head
Lange et al. Predictive vision based control of high speed industrial robot paths
CN113701766B (en) Robot map construction method, robot positioning method and device
JP3170345B2 (en) 3D information extraction method
WO2022141262A1 (en) Object detection
Brandt et al. Controlled active exploration of uncalibrated environments
CN119806157A (en) A robot recharging method, device, equipment and storage medium
Motai et al. SmartView: hand-eye robotic calibration for active viewpoint generation and object grasping
CN119169073A (en) Anchor hole positioning method and system
CN110877339B (en) Automatic control method and device in robot operation process

Legal Events

Date Code Title Description
FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20110713

Year of fee payment: 10

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

Free format text: PAYMENT UNTIL: 20120713

Year of fee payment: 11

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

Free format text: PAYMENT UNTIL: 20120713

Year of fee payment: 11

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

Free format text: PAYMENT UNTIL: 20130713

Year of fee payment: 12

EXPY Cancellation because of completion of term