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
JP2500654B2 - Three-dimensional shape relative position estimation device - Google Patents
[go: Go Back, main page]

JP2500654B2 - Three-dimensional shape relative position estimation device - Google Patents

Three-dimensional shape relative position estimation device

Info

Publication number
JP2500654B2
JP2500654B2 JP5288034A JP28803493A JP2500654B2 JP 2500654 B2 JP2500654 B2 JP 2500654B2 JP 5288034 A JP5288034 A JP 5288034A JP 28803493 A JP28803493 A JP 28803493A JP 2500654 B2 JP2500654 B2 JP 2500654B2
Authority
JP
Japan
Prior art keywords
relative position
range data
connection error
partial
estimating
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Expired - Fee Related
Application number
JP5288034A
Other languages
Japanese (ja)
Other versions
JPH06241731A (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.)
NEC Corp
Original Assignee
Nippon Electric Co Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Nippon Electric Co Ltd filed Critical Nippon Electric Co Ltd
Priority to JP5288034A priority Critical patent/JP2500654B2/en
Publication of JPH06241731A publication Critical patent/JPH06241731A/en
Application granted granted Critical
Publication of JP2500654B2 publication Critical patent/JP2500654B2/en
Anticipated expiration legal-status Critical
Expired - Fee Related legal-status Critical Current

Links

Landscapes

  • Length Measuring Devices By Optical Means (AREA)
  • Image Processing (AREA)
  • Image Analysis (AREA)

Description

【発明の詳細な説明】Detailed Description of the Invention

【0001】[0001]

【産業上の利用分野】この発明は、レンジデータ測定装
置により多視点から物体の形状を測定し、それぞれの視
点から得られた物体の部分的な形状(以下部分形状と呼
ぶ)を接続することにより物体の形状を再構成する装置
において、テクスチャを利用して部分形状間の相対位置
を推定する装置に関する。
BACKGROUND OF THE INVENTION 1. Field of the Invention The present invention measures the shape of an object from multiple points of view by means of a range data measuring device and connects the partial shapes of the object (hereinafter referred to as partial shapes) obtained from the respective points of view. The present invention relates to an apparatus for reconstructing the shape of an object by estimating the relative position between partial shapes using texture.

【0002】[0002]

【従来の技術】主なレンジデータ測定法として、ステレ
オ法、レンジファインダを用いる方法がある(尾上他編
「画像処理ハンドブック」昭晃堂(1987)。いずれ
の方法でも一回に取得できるのは物体の一部分(部分形
状)だけであり、物体全体の形状を再構成するには、さ
まざまな方向から物体の形状を測定し、その結果を統合
しなければならない。形状再構成には、部分形状間の相
対位置が必要であるが、ステレオカメラやレンジファイ
ンダの相対位置を求めることが困難であることにより、
相対位置を形状データから推定する方法が開発されてい
る。これらはいずれも隣接する部分形状が重なり領域
(共通領域)を持つことを前提として、重なり領域から
相対位置を推定する。
2. Description of the Related Art As a main range data measuring method, there are a method using a stereo method and a range finder (Ogami et al., "Image Processing Handbook", Shokoido (1987). Either method can be acquired at one time. It is only a part (partial shape) of the object, and in order to reconstruct the shape of the whole object, it is necessary to measure the shape of the object from various directions and integrate the results. Although it is necessary to find the relative position between the two, it is difficult to find the relative position of the stereo camera or range finder.
Methods have been developed for estimating relative position from shape data. In each of these, the relative position is estimated from the overlapping area on the assumption that adjacent partial shapes have an overlapping area (common area).

【0003】例えばレンジファインダにより得られた形
状を三角パッチで近似的に表現し、三角パッチの対応の
妥当な組み合せを求めることにより相対位置を推定する
方法(Robert Bergevin他:Estim
ating the 3DRigid Transfo
rmation Between Two Range
Views of a Complex Objec
t,ICPR’92)、また、レンジデータを分割し、
分割されてできた領域間のマッチングにより相対位置を
推定する方法(河井他:”多視点レンジデータからの3
次元形状復元”、電子情報通信学会パターン認識・理解
研究会技術報告91−19(1991)、物体表面上の
テクスチャを手がかりにして部分形状間の重なり領域を
求め、重なり領域の位置・向きを一致させることで相対
位置を求める方法(丸家、根本、高島:3次元形状再構
成のための部分形状接続法、1992年電子情報通信学
会春季全国大会)などがある。
For example, a method of approximating a shape obtained by a range finder by triangular patches and estimating a relative combination of the triangular patches to estimate the relative position (Robert Bergevin et al .: Estim).
aating the 3D RigTrans
rmation Between Two Range
Views of a Complex Objec
t, ICPR'92), and also divide the range data,
A method of estimating relative position by matching between divided regions (Kawai et al .: 3 from multi-view range data)
Dimensional shape restoration ", IEICE Technical Report 91-19 (1991), Pattern Recognition and Understanding Research Society of the Institute of Electronics, Information and Communication Engineers, using the texture on the surface of an object as a clue to find the overlapping area between partial shapes and match the position and orientation of the overlapping area. There is a method of obtaining a relative position by doing so (Maruya, Nemoto, Takashima: partial shape connection method for three-dimensional shape reconstruction, 1992 IEICE Spring National Convention).

【0004】[0004]

【発明が解決しようとする課題】しかしながら、形状に
関する情報だけで相対位置を推定するのは困難な場合が
ある。例えば、重なり領域が、曲率が急激に変化する、
あるいは不連続が存在するなどの特徴を含まない限り、
重なり領域が一意に決まらないか、あるいは決まったと
してもレンジデータに含まれているノイズの影響を受け
ている可能性が高い。
However, it may be difficult to estimate the relative position based only on the information about the shape. For example, the overlap region has a sharp change in curvature,
Or unless it includes features such as the presence of discontinuities,
There is a high possibility that the overlapping area is not uniquely determined or even if it is, it is affected by the noise included in the range data.

【0005】一方テクスチャを用いる方式では、相対位
置が一意に決まらないという問題は起こりにくい。しか
しながら前述の方法では、ステレオカメラの光軸が、常
にある水平面内になければならないという制約があっ
た。
On the other hand, in the method using the texture, the problem that the relative position is not uniquely determined is unlikely to occur. However, the method described above has a constraint that the optical axis of the stereo camera must always be in a certain horizontal plane.

【0006】本発明の目的は、テクスチャを用いること
で、安定な相対位置推定が行え、かつ位置、向きに関す
る制約の無い相対位置推定装置を提供することにある。
An object of the present invention is to provide a relative position estimating device which can perform stable relative position estimation by using a texture and has no restriction on the position and orientation.

【0007】[0007]

【課題を解決するための手段】上述した問題点を解決す
るため、第1の発明の3次元形状相対位置推定装置は、
レンジデータとテクスチャデータ、及びレンジデータ測
定装置の傾斜角を入力するデータ入力手段と、部分形状
を表現する座標系を変換してレンジデータ測定装置の傾
きの影響を解消する座標変換手段と、水平面を設定し、
部分形状と水平面との交線を求める交線計算手段と、部
分形状間の縦ずれを設定する縦ずれ設定手段と、前記交
線計算手段で得られた交線間の重なりを推定する重なり
幅推定手段と、部分形状間の水平面上の相対位置を推定
する相対位置推定手段と、接続誤差を計算する接続誤差
計算手段と、接続誤差が最小となる縦ずれを選択する最
小接続誤差選択手段と、を備えている。
In order to solve the above-mentioned problems, the three-dimensional shape relative position estimating device of the first invention is
Data input means for inputting range data and texture data, and an inclination angle of the range data measuring device, coordinate conversion means for converting a coordinate system expressing a partial shape to eliminate the influence of the inclination of the range data measuring device, and a horizontal plane. Set
An intersection line calculating means for obtaining an intersection line between a partial shape and a horizontal plane, a vertical deviation setting means for setting a vertical deviation between the partial shapes, and an overlapping width for estimating an overlap between the intersection lines obtained by the intersection line calculating means. An estimating means; a relative position estimating means for estimating a relative position between the partial shapes on a horizontal plane; a connection error calculating means for calculating a connection error; and a minimum connection error selecting means for selecting a vertical deviation that minimizes the connection error. , Are provided.

【0008】第2の発明の3次元形状相対位置推定装置
は、レンジデータとテクスチャデータ、及びレンジデー
タ測定装置の傾斜角を入力するデータ入力手段と、部分
形状を表現する座標系を変換してレンジデータ測定装置
の傾きの影響を解消する座標変換手段と、水平面を設定
し、部分形状と水平面との交線を求める交線計算手段
と、部分形状間の縦ずれを設定する縦ずれ設定手段と、
交線間の重なり幅の候補を抽出する交線重なり幅候補抽
出手段と、上下の重なり幅候補同士を照合し、整合性の
高い重なり幅を選択する重なり幅選択手段と、部分形状
間の水平面上の相対位置を推定する相対位置推定手段
と、接続誤差を計算する接続誤差計算手段と、接続誤差
が最小となる縦ずれを選択する最小接続誤差選択手段
と、を備えている。
The three-dimensional shape relative position estimating apparatus of the second invention converts the range data and the texture data, the data input means for inputting the inclination angle of the range data measuring apparatus, and the coordinate system expressing the partial shape. Coordinate conversion means for eliminating the influence of the tilt of the range data measuring device, intersection line calculation means for setting a horizontal plane to find the intersection line between the partial shape and the horizontal plane, and vertical deviation setting means for setting the vertical deviation between the partial shapes. When,
Intersection line overlap width candidate extraction means for extracting overlap width candidates between intersection lines, overlap width selection means for collating upper and lower overlap width candidates and selecting an overlap width with high consistency, and a horizontal plane between partial shapes It comprises a relative position estimating means for estimating the above relative position, a connection error calculating means for calculating a connection error, and a minimum connection error selecting means for selecting a vertical deviation that minimizes the connection error.

【0009】[0009]

【実施例】次に、本発明について図面を参照して説明す
る。
Next, the present invention will be described with reference to the drawings.

【0010】第1の発明について図面を参照して説明す
る。図1は、第1の発明の一実施例を示すブロック図で
ある。図1を参照すると、第1の発明は、データ入力手
段11と、部分形状生成手段12と、座標変換手段13
と、交線計算手段14と、縦ずれ設定手段15と、重な
り幅推定手段16と、相対位置推定手段17と、接続誤
差計算手段18と、最小接続誤差選択手段19から構成
される。
The first invention will be described with reference to the drawings. FIG. 1 is a block diagram showing an embodiment of the first invention. Referring to FIG. 1, the first aspect of the invention is a data input unit 11, a partial shape generation unit 12, and a coordinate conversion unit 13.
The intersection line calculation means 14, the vertical deviation setting means 15, the overlap width estimation means 16, the relative position estimation means 17, the connection error calculation means 18, and the minimum connection error selection means 19.

【0011】次に、図1から図8までを参照して、本実
施例の動作を説明する。本実施例において、入力された
データは、図1に示す矢印の順番で処理される。
Next, the operation of this embodiment will be described with reference to FIGS. 1 to 8. In the present embodiment, the input data is processed in the order of arrows shown in FIG.

【0012】データ入力手段11では物体のレンジデー
タ及びテクスチャデータを様々な視点測定し入力する。
その際、レンジデータ測定時のレンジデータ測定装置の
傾斜角も併せて入力する。図3に、レンジデータ測定装
置としてステレオカメラ31を用いた場合の例を示す。
図3において、レンジデータ測定装置としてのステレ
オカメラ31に、カメラの光軸が水平面となす角βを測
定する傾斜計32と、カメラの光軸を中心とした回転角
αを測定する傾斜計33が取り付けられている。なお、
このカメラを用いて得られるレンジデータは、カメラを
原点とした視点座標系(x1,y1,z1)で表現され
る。
The data input means 11 measures and inputs various range data and texture data of the object.
At that time, the tilt angle of the range data measuring device at the time of measuring the range data is also input. FIG. 3 shows an example in which a stereo camera 31 is used as the range data measuring device.
In FIG. 3, a stereo camera 31 as a range data measuring device includes an inclinometer 32 for measuring an angle β formed by the optical axis of the camera and a horizontal plane, and an inclinometer 33 for measuring a rotation angle α about the optical axis of the camera. Is attached. In addition,
Range data obtained using this camera is represented by a viewpoint coordinate system (x1, y1, z1) with the camera as the origin.

【0013】部分形状生成手段12では、レンジデータ
から部分形状を生成する。そのために、まずレンジデー
タの不連続な部分などを手がかりにして、背景部分のレ
ンジデータを消去する。そして、残ったレンジデータか
ら面を生成し、部分形状とする。
The partial shape generating means 12 generates a partial shape from the range data. Therefore, the range data in the background part is erased by using the discontinuous part of the range data as a clue. Then, a surface is generated from the remaining range data to form a partial shape.

【0014】座標変換手段13では、レンジデータ測定
時の、角αと角βの測定点毎の違いを解消するために、
カメラを原点とした視点座標系(x1,y1,z1)で
表現されている部分形状を、レンジデータ測定時の角α
と角βがゼロであった場合の座標系(z2,y2,z
2)に変換する。図3に、座標系(x1,y1,z1)
と(x2,y2,z2)の関係を示す。この変換を次に
述べる。
In the coordinate conversion means 13, in order to eliminate the difference between the measurement points of the angle α and the angle β when measuring the range data,
The partial shape represented by the viewpoint coordinate system (x1, y1, z1) with the camera as the origin is the angle α at the time of range data measurement.
And the angle β is zero, the coordinate system (z2, y2, z
Convert to 2). In Fig. 3, the coordinate system (x1, y1, z1)
And (x2, y2, z2) are shown. This conversion is described next.

【0015】図4において、oが原点である。座標系
(x1,y1,z1)では、水平面は、x1軸に対して
α、z1軸に対してβだけ傾いた面(abco)として
観測される。座標変換により、この面abcoを水平に
するには、まずz1軸を中心にして角−αだけ回転さ
せ、次に、x1軸を角γだけ回転させればよい(反時計
まわりの回転を正とする)。従って座標変換は、
In FIG. 4, o is the origin. In the coordinate system (x1, y1, z1), the horizontal plane is observed as a plane (abco) inclined by α with respect to the x1 axis and β with respect to the z1 axis. To make the plane abco horizontal by coordinate transformation, first rotate the z1 axis by an angle −α and then rotate the x1 axis by an angle γ (counterclockwise rotation is positive). And). Therefore, the coordinate transformation is

【0016】[0016]

【数1】 [Equation 1]

【0017】となる。## EQU1 ##

【0018】なお、γは次のように求められる。図4に
おいて、a’は、aをz1軸を中心に−α回転して得ら
れる点(sinα・tanβ,cosα・tanβ,
1)であり、b’は、bをz1軸を中心に−αだけ回転
して得られる点である。同様に、c’は、cをz1軸を
中心に−αだけ回転して得られる点である。さらに、平
面a’b’c’oがx1軸を中心にしたγの回転によ
り、(x1,z1)平面上に移ればよいのであるが、
a’b’と(y1,z1)平面との交点をdとすると、
dの位置は(0,cosα・tanβ,1)、線分od
の長さは、(1+cos2α・tan2β)1/2 と
なる。従って、 sinγ=cosα・tanβ/(1+cos2α・t
an2β)1/2 cosγ=1/(1+cos2α・tan2β)1/2 である。
Note that γ is obtained as follows. In FIG. 4, a ′ is a point (sinα · tanβ, cosα · tanβ, obtained by rotating α by −α around the z1 axis.
1), and b ′ is a point obtained by rotating b by −α around the z1 axis. Similarly, c ′ is a point obtained by rotating c by −α about the z1 axis. Furthermore, the plane a′b′c′o may be moved to the (x1, z1) plane by the rotation of γ around the x1 axis.
If the intersection of a'b 'and the (y1, z1) plane is d,
The position of d is (0, cos α · tan β, 1), the line segment od
Has a length of (1 + cos2α · tan2β) 1/2. Therefore, sinγ = cosα · tanβ / (1 + cos2α · t
an2β) 1/2 cosγ = 1 / (1 + cos2α · tan2β) 1/2.

【0019】レンジデータ測定時の角αと角βの違いを
以上の座標変換により解消すると、部分形状間の相対位
置の内、不明な成分は、垂直な方向に関する相対位置
h、水平な方向に関する相対位置p,q、垂直軸を中心
とする回転θとなる(図5)。
When the difference between the angle α and the angle β at the time of range data measurement is eliminated by the above coordinate conversion, an unknown component among the relative positions between the partial shapes is related to the vertical relative position h and the horizontal direction. The relative position is p, q, and the rotation is θ about the vertical axis (FIG. 5).

【0020】交線計算手段14では、水平面を等間隔に
複数設定し、部分形状と水平面との交線を求める(図
6)。水平面の数が多い方が正確な結果が得られる。
The intersecting line calculating means 14 sets a plurality of horizontal planes at equal intervals and obtains an intersecting line between the partial shape and the horizontal plane (FIG. 6). Accurate results can be obtained with a large number of horizontal planes.

【0021】縦ずれ設定手段15では、部分形状間の垂
直方向に関する位置ずれhを、あらかじめ決められた範
囲内で一つ設定する。縦ずれの単位は、水平面の間隔と
する。
The vertical deviation setting means 15 sets one positional deviation h between the partial shapes in the vertical direction within a predetermined range. The unit of vertical deviation is the interval between horizontal planes.

【0022】重なり幅推定手段16では、縦ずれ設定手
段15で設定した縦ずれにおいて、同じ高さとなる交線
間の重なり幅を推定する。そして、相対位置推定手段1
7で、推定された重なり領域が場所的に一致するような
相対位置を推定する。重なり幅推定手段16及び相対位
置推定手段17としては、前述の方法(丸家、根本、高
島:3次元形状再構成のための部分形状接続法、199
2年電子情報通信学会春季全国大会)が使える。
The overlap width estimating means 16 estimates the overlap width between intersecting lines having the same height in the vertical deviation set by the vertical deviation setting means 15. Then, the relative position estimating means 1
At 7, the relative position is estimated such that the estimated overlap areas are spatially coincident. The overlapping width estimating means 16 and the relative position estimating means 17 are the above-mentioned methods (Maruya, Nemoto, Takashima: partial shape connection method for three-dimensional shape reconstruction, 199.
2 years IEICE Spring National Convention) can be used.

【0023】具体的に述べると、まず、交線に沿って等
間隔(δt)にサンプル点を設定し、そこでの形状デー
タ及びテクスチャデータ(輝度値)をサンプリングす
る。次に同一水平面上にあるセグメント間の重なり幅を
δtずつ変化させながら、重なり部のテクスチャデータ
の相互相関を計算する(図7)。そして、相互相関が最
も高くなる重なり幅を求め、サンプル点同士を対応付け
るすべての水平面上の交線に対して行なう。図8に、重
なり幅が4、水平面が一つしかない場合の例を示す。図
では、a1とb1、a2とb2、a3とb3、a4とb
4が対応付られている。
More specifically, first, sample points are set at equal intervals (δt) along the intersection line, and the shape data and texture data (luminance values) there are sampled. Next, the cross-correlation of the texture data of the overlapping portion is calculated while changing the overlapping width between the segments on the same horizontal plane by δt (FIG. 7). Then, the overlapping width that maximizes the cross-correlation is obtained, and this is performed for all intersection lines on the horizontal plane that associate the sample points with each other. FIG. 8 shows an example in which the overlapping width is 4 and there is only one horizontal plane. In the figure, a1 and b1, a2 and b2, a3 and b3, a4 and b
4 is associated.

【0024】次に、対応付けられたサンプル点間の距離
の平均自乗が最も小さくなるような、部分形状の相対位
置(p,q,θ)を求める。
Next, the relative position (p, q, θ) of the partial shape is calculated so that the mean square of the distance between the associated sample points becomes the smallest.

【0025】接続誤差計算手段18では、部分形状の重
なった領域でのすきま(接続誤差と呼ぶ)を水平面毎に
計算する。まず、相対位置推定手段17で推定された
(p,q,θ)に従って、部分形状に対して平行移動、
回転を行なう。その時のサンプル点間の距離の平均自乗
を、その縦ずれにおける接続誤差とする。ここで、設定
された範囲内のすべての縦ずれについて接続誤差を計算
したかどうかを判定する。未計算の部分があれば、縦ず
れhを一定の値(水平面の間隔)だけ変化させて、縦ず
れ設定手段15にもどる。すべて接続誤差計算が終了し
ていれば、最小接続誤差選択手段19に進む。
The connection error calculating means 18 calculates the clearance (referred to as connection error) in the region where the partial shapes overlap, for each horizontal plane. First, in accordance with (p, q, θ) estimated by the relative position estimation means 17, a parallel movement with respect to the partial shape,
Rotate. The mean square of the distance between the sample points at that time is the connection error in the vertical shift. Here, it is determined whether or not the connection error has been calculated for all the vertical deviations within the set range. If there is an uncalculated portion, the vertical shift h is changed by a constant value (the interval between horizontal planes), and the process returns to the vertical shift setting means 15. If all the connection error calculations have been completed, the process proceeds to the minimum connection error selection means 19.

【0026】最小接続誤差選択手段19では、接続誤差
が最小となる縦ずれを検出し、その時の縦ずれh、およ
びその縦ずれ位置での(p,q,θ)を部分形状間の相
対位置であると推定する。
The minimum connection error selection means 19 detects the vertical deviation that minimizes the connection error, and determines the vertical deviation h at that time and (p, q, θ) at that vertical deviation position as relative positions between the partial shapes. Presumed to be

【0027】第2の発明について図面を参照して説明す
る。図2は、第2の発明の一実施例を示すブロック図で
ある。図2を参照すると、第2の発明は、データ入力手
段11と、部分形状生成手段12と、座標変換手段13
と、交線計算手段14と、縦ずれ設定手段15と、重な
り候補抽出手段21と、重なり幅選択手段22と、相対
位置推定手段17と、接続誤差計算手段18と、最小接
続誤差選択手段19から構成される。入力されたデータ
は、図1に示す矢印の順番で処理される。第1の発明と
の相違は、重なり候補抽出手段21と、重なり幅選択手
段22であるので、この二つについて説明する。
The second invention will be described with reference to the drawings. FIG. 2 is a block diagram showing an embodiment of the second invention. Referring to FIG. 2, the second aspect of the invention is a data input unit 11, a partial shape generation unit 12, and a coordinate conversion unit 13.
, Intersection line calculation means 14, vertical deviation setting means 15, overlap candidate extraction means 21, overlap width selection means 22, relative position estimation means 17, connection error calculation means 18, and minimum connection error selection means 19. Composed of. The input data is processed in the order of arrows shown in FIG. The difference from the first invention lies in the overlap candidate extracting means 21 and the overlap width selecting means 22, so these two will be explained.

【0028】第2の発明の重なり候補抽出手段21は、
交線毎に、テクスチャデータの相互相関が極大になる重
なり幅すべてを、重なり幅候補として抽出する。そし
て、重なり幅選択手段22において、上下の交線間で重
なり幅を照合し、最も整合性の高い重なり幅を選択す
る。 重なり幅選択手段22での処理を詳細に述べる。
まず、上下に隣接する交線上のサンプル点間で、対応付
けを行なう。お互いに最も近いサンプル点を1対1で対
応づけるが、通常交線の長さが異なるため、交線の端で
は対応先がないサンプル点が生じる。図9に、交線91
とそのすぐ下の交線92との対応付け、及び交線93と
そのすぐ下の交線94のとの対応付けの例を示す。ここ
で、w1は交線91と93での重なり幅候補、w2は交
線92と94での重なり幅候補である。また、S1は、
交線93の左端での相手先のないサンプル点数であり、
S2は、交線92の右端での相手先のないサンプル点数
である。この場合において、w1+S1 = w2−S
2の関係が成り立っていれば、重なり幅候補w1とw2
は整合性が高いと判定される。この整合性チェックをす
べての水平面上の交線で行い、より多くの水平面で整合
性かあると判定された重なり幅を選択する。
The overlapping candidate extracting means 21 of the second invention is
For each line of intersection, all overlapping widths at which the cross-correlation of texture data is maximized are extracted as overlapping width candidates. Then, the overlapping width selecting means 22 compares the overlapping widths between the upper and lower intersecting lines, and selects the overlapping width having the highest consistency. The processing in the overlapping width selection means 22 will be described in detail.
First, the sample points on the line of intersection that are vertically adjacent to each other are associated with each other. The sample points closest to each other are associated with each other in a one-to-one manner, but the lengths of the intersecting lines are usually different, so that there are sample points with no corresponding points at the ends of the intersecting lines. In FIG. 9, the intersection line 91
An example of the correspondence between the intersection line 92 and the intersection line 92 immediately below and the correspondence line 93 and the intersection line 94 immediately below the intersection line 93 is shown. Here, w1 is an overlapping width candidate at the intersecting lines 91 and 93, and w2 is an overlapping width candidate at the intersecting lines 92 and 94. Also, S1 is
It is the number of sample points at the left end of the intersection line 93 with no partner,
S2 is the number of sample points at the right end of the intersecting line 92 with no other party. In this case, w1 + S1 = w2-S
If the relationship of 2 holds, the overlapping width candidates w1 and w2
Are determined to be highly consistent. This consistency check is performed on intersection lines on all horizontal planes, and the overlapping width determined to be consistent on more horizontal planes is selected.

【0029】[0029]

【発明の効果】第1の発明の形状相対位置推定装置で
は、部分形状上のテクスチャを用いて部分形状の重なり
を検出し、2つの部分形状の共通な部分を検出する。共
通部分の形状を合わせることにより2つの部分形状の相
対位置を推定する。
According to the shape relative position estimating apparatus of the first invention, the texture of the partial shapes is used to detect the overlapping of the partial shapes, and the common portion of the two partial shapes is detected. The relative positions of the two partial shapes are estimated by matching the shapes of the common parts.

【0030】そのため、曲率が急激に変化する、あるい
は不連続が存在するなどの特徴がなくても、相対位置を
決めることができる。また、レンジデータ測定装置の位
置、向きに関する制約が無いため、レンジデータ測定の
死角及び手間が大幅に削減される。
Therefore, the relative position can be determined even if there is no characteristic such as a rapid change in curvature or the presence of discontinuity. Further, since there is no restriction on the position and orientation of the range data measuring device, the blind spot and time and effort for measuring the range data are significantly reduced.

【0031】第2の発明の形状相対位置推定装置では、
重なり幅推定において整合性のチェックを行なうので、
より高精度な推定が行える。
In the shape relative position estimating device of the second invention,
Since the consistency is checked in the overlap width estimation,
More accurate estimation can be performed.

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

【図1】第1の発明の一実施例を示すブロック図であ
る。
FIG. 1 is a block diagram showing an embodiment of a first invention.

【図2】第2の発明の一実施例を示すブロック図であ
る。
FIG. 2 is a block diagram showing an embodiment of the second invention.

【図3】ステレオカメラと、ステレオカメラに対する傾
斜計の取り付け方を示す図である。
FIG. 3 is a diagram showing a stereo camera and how to attach an inclinometer to the stereo camera.

【図4】座標系間の関係を示す図である。FIG. 4 is a diagram showing a relationship between coordinate systems.

【図5】部分形状間の相対位置を表わすパラメータを示
す図である。
FIG. 5 is a diagram showing parameters representing relative positions between partial shapes.

【図6】部分形状と水平面との関係を示す図である。FIG. 6 is a diagram showing a relationship between a partial shape and a horizontal plane.

【図7】重なり幅推定を示す図である。FIG. 7 is a diagram showing overlap width estimation.

【図8】セグメント間の相対位置を表わすパラメータを
示す図である。
FIG. 8 is a diagram showing parameters representing relative positions between segments.

【図9】サンプル点間の上下の対応付けを示す図であ
る。
FIG. 9 is a diagram showing upper and lower correspondences between sample points.

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

11 データ入力手段 12 部分形状生成手段 13 座標変換手段 14 交線計算手段 15 縦ずれ設定手段 16 重なり幅推定手段 17 相対位置推定手段 18 接続誤差計算手段 19 最小接続誤差選択手段 21 重なり幅候補抽出手段 22 重なり幅選択手段 91 交線 92 交線 93 交線 94 交線 11 data input means 12 partial shape generation means 13 coordinate conversion means 14 intersection line calculation means 15 vertical deviation setting means 16 overlap width estimation means 17 relative position estimation means 18 connection error calculation means 19 minimum connection error selection means 21 overlap width candidate extraction means 22 Overlap width selection means 91 Intersection line 92 Intersection line 93 Intersection line 94 Intersection line

Claims (2)

(57)【特許請求の範囲】(57) [Claims] 【請求項1】 レンジデータ測定装置により、複数の視
点から物体のレンジデータ、テクスチャデータを入力
し、そこから得られた物体の部分形状間の相対位置を推
定する3次元形状相対位置推定装置において、レンジデ
ータとテクスチャデータ、及びレンジデータ測定装置の
傾斜角を入力するデータ入力手段と、部分形状を表現す
る座標系を変換してレンジデータ測定装置の傾きの影響
を解消する座標変換手段と、水平面を設定し、部分形状
と水平面との交線を求める交線計算手段と、部分形状間
の縦ずれを設定する縦ずれ設定手段と、前記交線計算手
段で得られた交線間の重なりを推定する重なり幅推定手
段と、部分形状間の水平面上の相対位置を推定する相対
位置推定手段と、接続誤差を計算する接続誤差計算手段
と、接続誤差が最小となる縦ずれを選択する最小接続誤
差選択手段と、を備えたことを特徴とする3次元形状相
対位置推定装置。
1. A three-dimensional shape relative position estimation device for inputting range data and texture data of an object from a plurality of viewpoints by a range data measuring device and estimating relative positions between partial shapes of the object obtained therefrom. Data input means for inputting the range data and texture data, and the inclination angle of the range data measuring device, and coordinate conversion means for converting the coordinate system expressing the partial shape to eliminate the influence of the inclination of the range data measuring device, An intersection line calculation means for setting a horizontal plane and obtaining an intersection line between the partial shape and the horizontal plane, a vertical deviation setting means for setting a vertical deviation between the partial shapes, and an overlap between the intersection lines obtained by the intersection line calculation means. Of the overlap width, a relative position estimating means for estimating the relative position between the partial shapes on the horizontal plane, a connection error calculating means for calculating the connection error, and a minimum connection error. And a minimum connection error selecting means for selecting the vertical deviation.
【請求項2】 前記重なり幅推定手段に代わって、交線
間の重なり幅の候補を抽出する交線重なり幅候補抽出手
段と、 上下交線間で重なり幅候補同士を照合し、整合性の高い
重なり幅を選択する重なり幅選択手段と、を備えたこと
を特徴とする請求項1記載の3次元形状相対位置推定装
置。
2. An intersection line overlap width candidate extraction unit that extracts candidates for an overlap width between intersection lines instead of the overlap width estimation unit, and an overlap width candidate between upper and lower intersection lines are collated to check consistency. 3. The three-dimensional shape relative position estimating device according to claim 1, further comprising an overlapping width selecting unit that selects a high overlapping width.
JP5288034A 1992-11-20 1993-11-17 Three-dimensional shape relative position estimation device Expired - Fee Related JP2500654B2 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP5288034A JP2500654B2 (en) 1992-11-20 1993-11-17 Three-dimensional shape relative position estimation device

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
JP4-311841 1992-11-20
JP31184192 1992-11-20
JP5288034A JP2500654B2 (en) 1992-11-20 1993-11-17 Three-dimensional shape relative position estimation device

Publications (2)

Publication Number Publication Date
JPH06241731A JPH06241731A (en) 1994-09-02
JP2500654B2 true JP2500654B2 (en) 1996-05-29

Family

ID=26556990

Family Applications (1)

Application Number Title Priority Date Filing Date
JP5288034A Expired - Fee Related JP2500654B2 (en) 1992-11-20 1993-11-17 Three-dimensional shape relative position estimation device

Country Status (1)

Country Link
JP (1) JP2500654B2 (en)

Families Citing this family (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP3869876B2 (en) 1995-12-19 2007-01-17 キヤノン株式会社 Image measuring method and image measuring apparatus
JP3712847B2 (en) * 1997-05-20 2005-11-02 株式会社リコー Three-dimensional shape measurement method, three-dimensional shape measurement device, and posture detection device for imaging means
JP2003269928A (en) 2002-03-12 2003-09-25 Nec Corp Three-dimensional shape measuring method and apparatus and program
JP4985213B2 (en) * 2007-08-20 2012-07-25 日本電気株式会社 Three-dimensional shape measuring method, apparatus and program
WO2012126522A1 (en) * 2011-03-23 2012-09-27 Metaio Gmbh Method for registering at least one part of a first and second image using a collineation warping function

Also Published As

Publication number Publication date
JPH06241731A (en) 1994-09-02

Similar Documents

Publication Publication Date Title
Scharstein View synthesis using stereo vision
Garcia et al. A combined temporal tracking and stereo-correlation technique for accurate measurement of 3D displacements: application to sheet metal forming
KR102424135B1 (en) Structured light matching of a set of curves from two cameras
CN110044374B (en) Image feature-based monocular vision mileage measurement method and odometer
JP2009041972A (en) Image processing apparatus and method
CN112771573A (en) Depth estimation method and device based on speckle images and face recognition system
JP2953154B2 (en) Shape synthesis method
JPH08136220A (en) Article position detection method and apparatus
Yang et al. Error analysis of 3D shape construction from structured lighting
Jaw et al. Registration of ground‐based LiDAR point clouds by means of 3D line features
JP5462662B2 (en) Position / orientation measurement apparatus, object identification apparatus, position / orientation measurement method, and program
JP2500654B2 (en) Three-dimensional shape relative position estimation device
JP2006214893A (en) Computer software program for measuring the three-dimensional shape of an object using the object measuring method and computer system
JPH07103715A (en) Method and apparatus for recognizing three-dimensional position and attitude based on visual sense
Maas Automatic DEM generation by multi-image feature based matching
Belhaoua et al. Error evaluation in a stereovision-based 3D reconstruction system
JP4935769B2 (en) Plane region estimation apparatus and program
JP6773280B2 (en) Tunnel excavation management method and management equipment
Markiewicz Evaluation of 2d affine—hand-crafted detectors for feature-based tls point cloud registration
CN101939666B (en) Method for the computer-aided calculation of the movement of an object using sensor data
CN113140031A (en) Three-dimensional image modeling system and method and oral cavity scanning equipment applying same
JP6835665B2 (en) Information processing equipment and programs
JP7420231B2 (en) Elevator 3D data processing equipment
JP2023009981A (en) Facility recognition system and facility recognition method
JP4153322B2 (en) Method and apparatus for associating measurement points in photogrammetry

Legal Events

Date Code Title Description
A01 Written decision to grant a patent or to grant a registration (utility model)

Free format text: JAPANESE INTERMEDIATE CODE: A01

Effective date: 19960123

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

Free format text: PAYMENT UNTIL: 20080313

Year of fee payment: 12

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

Free format text: PAYMENT UNTIL: 20090313

Year of fee payment: 13

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

Free format text: PAYMENT UNTIL: 20090313

Year of fee payment: 13

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

Free format text: PAYMENT UNTIL: 20100313

Year of fee payment: 14

LAPS Cancellation because of no payment of annual fees