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
JP3296099B2 - Corresponding point position measurement method and apparatus for correlation matching - Google Patents
[go: Go Back, main page]

JP3296099B2 - Corresponding point position measurement method and apparatus for correlation matching - Google Patents

Corresponding point position measurement method and apparatus for correlation matching

Info

Publication number
JP3296099B2
JP3296099B2 JP17802294A JP17802294A JP3296099B2 JP 3296099 B2 JP3296099 B2 JP 3296099B2 JP 17802294 A JP17802294 A JP 17802294A JP 17802294 A JP17802294 A JP 17802294A JP 3296099 B2 JP3296099 B2 JP 3296099B2
Authority
JP
Japan
Prior art keywords
point
correlation value
correlation
equation
model
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
JP17802294A
Other languages
Japanese (ja)
Other versions
JPH0843024A (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.)
Meidensha Corp
Original Assignee
Meidensha 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 Meidensha Corp filed Critical Meidensha Corp
Priority to JP17802294A priority Critical patent/JP3296099B2/en
Publication of JPH0843024A publication Critical patent/JPH0843024A/en
Application granted granted Critical
Publication of JP3296099B2 publication Critical patent/JP3296099B2/en
Anticipated expiration legal-status Critical
Expired - Fee Related legal-status Critical Current

Links

Landscapes

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

Description

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

【0001】[0001]

【産業上の利用分野】本発明は、いわゆる画像処理装置
における相関マッチング処理により、入力画像中のモデ
ルとの対応点の位置(座標)を決定する技術に関し、特
には、入力画像の画素の分解能よりも精密に位置を計測
できる方法及び装置に関する。
BACKGROUND OF THE INVENTION 1. Field of the Invention The present invention relates to a technique for determining a position (coordinate) of a corresponding point with a model in an input image by a correlation matching process in a so-called image processing apparatus, and more particularly, to a resolution of a pixel of the input image. The present invention relates to a method and an apparatus capable of measuring a position more precisely.

【0002】[0002]

【従来の技術】画像処理装置を用いて或る対象物の検査
を行う場合、対象物の画像中での位置(座標)を計測す
る必要がある。このような目的のために、画像処理装置
においては、相関マッチング処理が広く利用されてい
る。相関マッチング処理では、 (1)まず、検査前の予備段階として、対象物のサンプ
ルの画像の一部を切り出してその部分の濃度に関するパ
ターンを、モデルとして記憶しておく。 (2)次に検査時に、入力画像から1画素ずれた各点毎
に順次モデルと同じ形状で且つ同じ大きさの濃度に関す
るパターンを切り出し、各点毎にモデルと切り出したパ
ターンとの相関値を計算し、相関の最も強いパターンが
切り出された点の位置をモデルとの対応点の位置とす
る。なお、相関値の大小と相関の強弱は相関値を求める
計算式により異なり、相関値が最小の場合に相関が最強
となる計算式と、相関値が最大の場合に相関が最強とな
る計算式とがある。
2. Description of the Related Art When an object is inspected using an image processing apparatus, it is necessary to measure a position (coordinate) of the object in an image. For such a purpose, correlation matching processing is widely used in image processing apparatuses. In the correlation matching process, (1) First, as a preliminary stage before inspection, a part of an image of a sample of an object is cut out, and a pattern relating to the density of the part is stored as a model. (2) Next, at the time of inspection, a pattern relating to density having the same shape and the same size as the model is sequentially cut out for each point shifted by one pixel from the input image, and the correlation value between the model and the cut out pattern is obtained for each point Calculation is performed, and the position of the point where the pattern having the strongest correlation is cut out is set as the position of the corresponding point with the model. The magnitude of the correlation value and the strength of the correlation are different depending on the calculation formula for calculating the correlation value. The calculation formula that gives the strongest correlation when the correlation value is the minimum and the calculation formula that gives the strongest correlation when the correlation value is the maximum. There is.

【0003】下式(9)は2つのパターン間の相関値を
求める計算式の一例であり、式中でPi はモデルのi番
目の画素の濃度であり、Qi は切り出したi番目の画素
の濃度であり、この式(9)の場合は、切り出した各パ
ターンのうち、相関値Cが最小となるパターンを切り出
した点の位置を、対応点の位置として決定することにな
る。
The following equation (9) is an example of an equation for calculating a correlation value between two patterns. In the equation, P i is the density of the i-th pixel of the model, and Q i is the i-th pixel cut out. In the case of Expression (9), the position of a point at which the pattern with the smallest correlation value C is extracted from the extracted patterns is determined as the position of the corresponding point.

【数9】 (Equation 9)

【0004】また、下式(10)は相関値の計算式の他
の例であり、この場合は、切り出したパターンのうち、
相関値Cが最大になるパターンを切り出した点の位置
を、対応点の位置として決定することになる。
The following equation (10) is another example of the equation for calculating the correlation value. In this case, among the extracted patterns,
The position of the point at which the pattern with the maximum correlation value C is cut out is determined as the position of the corresponding point.

【数10】 (Equation 10)

【0005】ここでパターンを切り出す点と、モデルと
の対応点について、図8を例にして説明する。図8に示
した例では、家を入力画像から検出する対象物とし、家
の位置を屋根の頂点の位置で代表するものと仮定して説
明する。即ち、図8において、この例のモデル100は
家101の画像を含む四角形の濃度パターンであり、モ
デル100内の各点の位置はそのうち任意の1点を基準
にとれば、その基準点からの位置で表現できる。そこ
で、屋根の頂点102が家101を代表するので、この
頂点102を基準点とする。これに対し、入力画像20
0からモデル100と同じ形状で且つ同じ大きさのパタ
ーン201を或る点202で切り出す際は、パターン2
01と点202との位置関係が、モデル100とその基
準点(屋根の頂点102)との位置関係に一致する。そ
して、順次1画素ずつ点202をずらしてパターンの切
り出しを行った後、或るパターン201Aでモデル10
0との相関が最も強いことが判れば、そのパターン20
1Aを切り出した点202Aがモデル100の屋根の頂
点102に対応するから、モデル100との対応点とな
る。従って、相関が最も強いパターン201Aを切り出
した点202Aの位置が入力画像200におけるモデル
との対応点の位置となり、入力画像200の中の家20
3の位置を代表する。
[0005] Here, a point where a pattern is cut out and a corresponding point with a model will be described with reference to FIG. 8 as an example. In the example illustrated in FIG. 8, a description is given on the assumption that a house is an object to be detected from an input image, and that the position of the house is represented by the position of the top of the roof. That is, in FIG. 8, the model 100 of this example is a rectangular density pattern including an image of the house 101, and the position of each point in the model 100 is determined from any one of the points. Can be expressed by position. Accordingly, since the vertex 102 of the roof represents the house 101, the vertex 102 is set as a reference point. In contrast, the input image 20
When cutting out a pattern 201 having the same shape and the same size as the model 100 from 0 at a certain point 202, a pattern 2
The positional relationship between 01 and the point 202 matches the positional relationship between the model 100 and its reference point (the roof vertex 102). Then, after the point 202 is sequentially shifted by one pixel to cut out the pattern, the model 10 is extracted with a certain pattern 201A.
If the correlation with 0 is found to be the strongest, the pattern 20
Since a point 202A obtained by cutting 1A corresponds to the vertex 102 of the roof of the model 100, the point 202A is a corresponding point with the model 100. Accordingly, the position of the point 202A obtained by cutting out the pattern 201A having the highest correlation is the position of the corresponding point with the model in the input image 200, and the house 20 in the input image 200
3 position.

【0006】[0006]

【発明が解決しようとする課題】上述した相関マッチン
グ処理は、常識的で広く使用されているが、厳密に言う
と以下の点で問題がある。即ち、1画素ずつずれた各点
毎にパターンを切り出していることから、式(9)また
は式(10)の計算から得られる対応点の位置は離散的
なものであり、場合によっては真の対応点から最大で±
0.5画素分ずれた位置になる危険性がある。この0.
5画素のずれが実際の対象物上でどのような寸法になる
かは、撮影に用いるビデオカメラの分解能等に依存する
ので、従来の相関マッチング処理では対処することがで
きないものである。しかしながら、実用的には0.5画
素のずれは影響が大きく、精密な位置決めには使用する
ことができない。例えば、ビデオカメラの分解能が水平
方向に512画素で視角が60°、ビデオカメラと対象
物間の距離が2mである場合は、±0.5画素のずれ
は、約5mmとなる。
The above-described correlation matching processing is common sense and widely used, but strictly speaking, has the following problems. That is, since the pattern is cut out for each point shifted by one pixel, the position of the corresponding point obtained from the calculation of Expression (9) or Expression (10) is discrete, and in some cases, the true position is obtained. ± at maximum from corresponding point
There is a risk that the position will be shifted by 0.5 pixel. This 0.
The size of the displacement of the five pixels on the actual object depends on the resolution of the video camera used for photographing, and cannot be dealt with by the conventional correlation matching processing. However, in practice, a shift of 0.5 pixel has a large effect and cannot be used for precise positioning. For example, when the resolution of the video camera is 512 pixels in the horizontal direction and the viewing angle is 60 °, and the distance between the video camera and the object is 2 m, the deviation of ± 0.5 pixels is about 5 mm.

【0007】本発明は上述した従来技術に鑑み、入力画
像の分解能よりも精密に対応点の位置を計測することが
できる方法及び装置を提供することを目的とする。
SUMMARY OF THE INVENTION In view of the above-mentioned prior art, an object of the present invention is to provide a method and an apparatus capable of measuring the position of a corresponding point more precisely than the resolution of an input image.

【0008】[0008]

【課題を解決するための手段】上述の目的を達成する第
1の発明に係る相関マッチングの対応点位置計測方法
は、入力画像の各点毎に、モデルと同じ形状で且つ同じ
大きさのパターンを切り出し、モデルと切り出したパタ
ーンとの相関値を相関が強いほど相関値が小さくなる計
算式で計算すること;各点での相関値の計算結果から、
相関が最も強い第1の点の位置A0 と、同第1の点での
相関値C0 と、第1の点を含む一次元上で第1の点の直
前に位置する第2の点での相関値C-1と、同一次元上で
第1の点の直後に位置する第3の点での相関値C+1とを
求めること;第2の点での相関値C-1と第3の点での相
関値C+1との大小関係を判定すること;判定結果から、
-1>C+1の場合は下式(11)により、C-1<C+1
場合は下式(12)により、入力画像におけるモデルと
の対応点の位置Ar を補間演算で求めること;により、
画像分解能よりも精密に対応点の位置を計測すること
特徴とする相関マッチングの対応点位置計測方法であ
る。
According to a first aspect of the present invention, there is provided a method for measuring a corresponding point position of a correlation matching, wherein each pattern of an input image has a pattern having the same shape and the same size as a model. Is calculated, and the correlation value between the model and the cut pattern is calculated by a calculation formula in which the stronger the correlation, the smaller the correlation value is. From the calculation result of the correlation value at each point,
The position A 0 of the first point having the strongest correlation, the correlation value C 0 at the first point, and the second point located immediately before the first point on one dimension including the first point a correlation value C -1 at a second point; the correlation values C -1 in, obtaining a correlation value C +1 at the third point located immediately after the first point on the same dimension Determining the magnitude relationship with the correlation value C + 1 at the third point; from the determination result,
When C -1 > C +1 , the position Ar of the corresponding point with the model in the input image is interpolated by the following equation (11), and when C -1 <C +1 , the following equation (12). by; obtaining
This is a corresponding point position measuring method for correlation matching, characterized in that the position of the corresponding point is measured more precisely than the image resolution .

【0009】[0009]

【数11】 Ar =A0 −(1/2)(C-1−C+1)/(C0 −C-1) …式(11)A r = A 0 − (1/2) (C −1 −C +1 ) / (C 0 −C −1 ) Equation (11)

【0010】[0010]

【数12】 Ar =A0 +(1/2)(C+1−C-1)/(C0 −C+1) …式(12)A r = A 0 + (1 /) (C +1 −C −1 ) / (C 0 −C +1 ) Equation (12)

【0011】また、第2の発明に係る相関マッチングの
対応点位置計測方法は、入力画像の各点毎に、モデルと
同じ形状で且つ同じ大きさのパターンを切り出し、モデ
ルと切り出したパターンとの相関値を相関が強いほど相
関値が大きくなる計算式で計算すること;各点での相関
値の計算結果から、相関が最も強い第1の点の位置A0
と、同第1の点での相関値C0 と、第1の点を含む一次
元上で第1の点の直前に位置する第2の点での相関値C
-1と、同一次元上で第1の点の直後に位置する第3の点
での相関値C+1とを求めること;第2の点での相関値C
-1と第3の点での相関値C+1との大小関係を判定するこ
と;判定結果から、C-1<C+1の場合は前記式(11)
により、C-1>C+1の場合は前記式(12)により、入
力画像におけるモデルとの対応点の位置Ar を補間演算
で求めること;により、画像分解能よりも精密に対応点
の位置を計測することを特徴とする相関マッチングの対
応点位置計測方法である。
Further, in the method for measuring the corresponding point position of the correlation matching according to the second invention, a pattern having the same shape and the same size as the model is cut out for each point of the input image, and the pattern is compared with the model. Calculating the correlation value using a calculation formula in which the correlation value increases as the correlation increases; from the calculation result of the correlation value at each point, the position A 0 of the first point having the strongest correlation is calculated.
And a correlation value C 0 at the first point and a correlation value C 2 at a second point located immediately before the first point on one dimension including the first point.
-1 and a correlation value C + 1 at a third point located immediately after the first point on the same dimension; a correlation value C at the second point
Determining the magnitude relationship between -1 and the correlation value C +1 at the third point; from the determination result, if C -1 <C +1 , the above equation (11) is used.
Accordingly, C -1> C according to the equation in the case of +1 (12), determining the position A r of the corresponding point of the model in the input image by interpolation calculation; by, precisely corresponding points than the image resolution
Is a corresponding point position measuring method of correlation matching characterized by measuring the position of a corresponding point.

【0012】更に、第3の発明に係る相関マッチングの
対応点位置計測装置は、入力画像の各点毎に、モデルと
同じ形状で且つ同じ大きさのパターンを切り出し、モデ
ルと切り出したパターンとの相関値を相関が強いほど相
関値が小さくなる計算式で計算する手段と;各点での相
関値の計算結果から、相関が最も強い第1の点の位置A
0 と、同第1の点での相関値C0 と、第1の点を含む一
次元上で第1の点の直前に位置する第2の点での相関値
-1と、同一次元上で第1の点の直後に位置する第3の
点での相関値C+1とを求める手段と;前記2つの式(1
1)及び式(12)を格納している手段と;第2の点で
の相関値C-1と第3の点での相関値C+1との大小関係を
判定する手段と;判定結果から、C-1>C+1の場合は前
記式(11)を選択し、C-1<C+1の場合は前記式(1
2)を選択する手段と;選択された式により、入力画像
におけるモデルとの対応点の位置Ar を補間演算で求め
る手段と;を備え、画像分解能よりも精密に対応点の位
置を計測することを特徴とする相関マッチングの対応点
位置計測装置である。
Further, the corresponding point position measuring apparatus for correlation matching according to the third invention cuts out a pattern having the same shape and the same size as the model for each point of the input image, and compares the pattern with the model. Means for calculating the correlation value by a calculation formula in which the stronger the correlation is, the smaller the correlation value is; and calculating the position A of the first point having the strongest correlation from the calculation result of the correlation value at each point.
0 , a correlation value C 0 at the first point, and a correlation value C −1 at a second point immediately before the first point on one dimension including the first point, in the same dimension Means for obtaining a correlation value C +1 at a third point located immediately after the first point above;
Means for storing 1) and equation (12); means for determining the magnitude relationship between the correlation value C -1 at the second point and the correlation value C +1 at the third point; determination results Therefore, when C -1 > C +1, the formula (11) is selected. When C -1 <C +1 , the formula (1) is selected.
Means for selecting 2); means for obtaining the position Ar of the corresponding point with the model in the input image by interpolation by the selected formula; and determining the position of the corresponding point more precisely than the image resolution.
A corresponding point position measuring device for correlation matching characterized by measuring the position.

【0013】また、第4の発明に係る相関マッチングの
対応点位置計測装置は、入力画像の各点毎に、モデルと
同じ形状で且つ同じ大きさのパターンを切り出し、モデ
ルと切り出したパターンとの相関値を相関が強いほど相
関値が大きくなる計算式で計算する手段と;各点での相
関値の計算結果から、相関が最も強い第1の点の位置A
0 と、同第1の点での相関値C0 と、第1の点を含む一
次元上で第1の点の直前に位置する第2の点での相関値
-1と、同一次元上で第1の点の直後に位置する第3の
点での相関値C+1とを求める手段と;前記2つの式(1
1)及び式(12)を格納している手段と;第2の点で
の相関値C-1と第3の点での相関値C+1との大小関係を
判定する手段と;判定結果から、C-1<C+1の場合は前
記式(11)を選択し、C-1>C+1の場合は前記式(1
2)を選択する手段と;選択された式により、入力画像
におけるモデルとの対応点の位置Ar を補間演算で求め
る手段と;を備え、画像分解能よりも精密に対応点の位
置を計測することを特徴とする相関マッチングの対応点
位置計測装置である。
Further, in the apparatus for measuring the corresponding point position of the correlation matching according to the fourth invention, a pattern having the same shape and the same size as the model is cut out for each point of the input image. Means for calculating the correlation value using a calculation formula in which the stronger the correlation is, the larger the correlation value is; and calculating the position A of the first point having the strongest correlation from the calculation result of the correlation value at each point.
0 , a correlation value C 0 at the first point, and a correlation value C −1 at a second point immediately before the first point on one dimension including the first point, in the same dimension Means for obtaining a correlation value C +1 at a third point located immediately after the first point above;
Means for storing 1) and equation (12); means for determining the magnitude relationship between the correlation value C -1 at the second point and the correlation value C +1 at the third point; determination results Therefore, when C -1 <C +1 , the above equation (11) is selected. When C -1 > C +1 , the above equation (1) is selected.
Means for selecting 2); means for obtaining the position Ar of the corresponding point with the model in the input image by interpolation by the selected formula; and determining the position of the corresponding point more precisely than the image resolution.
A corresponding point position measuring device for correlation matching characterized by measuring the position.

【0014】[0014]

【作用】入力画像中でモデルとの真の対応点の近傍で
は、図3に示すように、真の対応点Ar を含む一次元即
ち直線11上の前後で、相関の強弱は真の対応点Ar
中心に良い近似で対称に変化する。即ち、直線11と直
交する軸12上に相関値Cをとれば、真の対応点Ar
り前側での相関値の変化直線13と、後側での相関値の
変化直線14とは、互いに傾きの符号が逆で、その絶対
値が同一の傾きである。なお、図3では相関値が小さい
ほど相関が強いものとしているが、逆の場合でも、同様
である。
In the vicinity of the true corresponding point with the model in the input image, as shown in FIG. 3, one-dimensional including the true corresponding point Ar, that is, before and after on the straight line 11, the strength of the correlation is the true correspondence. It changes symmetrically with a good approximation around the point Ar . That is, taking the correlation value C on the axis 12 perpendicular to the straight line 11, the change in linear 14 correlation values in the true and change linearly 13 of the correlation value in front of the corresponding points A r, rear, together The signs of the slopes are opposite, and their absolute values are the same. In FIG. 3, it is assumed that the smaller the correlation value is, the stronger the correlation is.

【0015】上述した対称性を利用して補間すれば、真
の対応点Ar に実質的に等しい精密な対応点を求めるこ
とができる。この補間には、前述した式(11)及び式
(12)を用いると、これらは比較的単純な式なので、
補間演算に要する時間が短いという利点がある。
By performing interpolation using the above-described symmetry, a precise corresponding point substantially equal to the true corresponding point Ar can be obtained. If the above-described equations (11) and (12) are used for this interpolation, these are relatively simple equations.
There is an advantage that the time required for the interpolation calculation is short.

【0016】<式(11)及び(12)について>次
に、図4〜図7を参照して式(11)及び式(12)の
導出過程と、その使い分けを説明する。
<Equations (11) and (12)> Next, the process of deriving the equations (11) and (12) and their proper use will be described with reference to FIGS.

【0017】まず、前記式(9)など、相関が強いほど
相関値Cが小さくなる計算式を用いる場合、AとCとの
直交座標系を導入すると、図4または図5に示すよう
に、パターンを切り出した点の位置とその点での相関値
との対(A0 ,C0 ),(A-1,C-1),(A+1
+1)をプロットすることかできる。但し、A0 は相関
が最も強い第1の点の位置、C0 は第1の点での相関
値、A-1は第1の点を含む一次元A上で第1の点の直前
に位置する第2の点の位置、C-1は第2の点での相関
値、A+1は同一次元A上で第1の点の直後に位置する第
3の点の位置、C+1は第3の点での相関値をそれぞれ示
している。
First, when using a calculation formula such as the above formula (9) in which the correlation value C becomes smaller as the correlation becomes stronger, if an orthogonal coordinate system of A and C is introduced, as shown in FIG. 4 or FIG. Pairs (A 0 , C 0 ), (A −1 , C −1 ), (A +1 ,
C +1 ) can be plotted. Here, A 0 is the position of the first point having the strongest correlation, C 0 is the correlation value at the first point, and A −1 is the one-dimensional A including the first point immediately before the first point. The position of the second point located, C -1 is the correlation value at the second point, A +1 is the position of the third point located immediately after the first point on the same dimension A, C +1 Indicates the correlation value at the third point.

【0018】図4はC-1>C+1となる場合の例であり、
この場合は、(A0 ,C0 )と(A -1,C-1)を通る実
線の直線15を引き、この直線15と符号が反対で絶対
値が等しい傾きを持ち且つ(A+1,C+1)を通る二点鎖
線の直線16との交点を求めれば、そのA座標が対応点
の位置Ar となる。この場合、直線15は次式(13)
で表わされる。従って、直線16は次式(14)で表さ
れる。両直線の交点を求めるために式(13),(1
4)からCを消去すると、次式(15)が得られる。式
(15)中のAは対応点の位置Ar にほかならないか
ら、式(15)を展開すると、次式(16)を得る。そ
して、パターンを切り出す点は1画素ずつずれるので、
0 −A-1=1画素、A+1−A0 =1画素、(A-1+A
+1)/2=A0 であることから、式(16)を整理する
と、前記の式(11)が得られる。即ち、相関値を計算
式(9)等で求める場合は、C-1>C+1であれば、前式
(11)を用いた補間演算によりA0 と、C0 ,C-1
びC+1とを用いて精密な対応点の位置Ar を求めること
ができる。
FIG. 4 shows C-1> C+1This is an example where
In this case, (A0, C0) And (A -1, C-1) Passing through
Draw a straight line 15 of the line
Values have equal slopes and (A+1, C+1) Two-point chain passing
If the intersection of the line with the straight line 16 is found, the A coordinate is the corresponding point.
Position ArBecomes In this case, the straight line 15 is expressed by the following equation (13).
Is represented by Therefore, the straight line 16 is expressed by the following equation (14).
It is. Equations (13) and (1) are used to determine the intersection of the two straight lines.
By eliminating C from 4), the following equation (15) is obtained. formula
A in (15) is the position A of the corresponding pointrIs it nothing but
By expanding the equation (15), the following equation (16) is obtained. So
Then, since the point where the pattern is cut out is shifted by one pixel,
A0-A-1= 1 pixel, A+1-A0= 1 pixel, (A-1+ A
+1) / 2 = A0Therefore, formula (16) is rearranged.
And the above equation (11) is obtained. That is, calculate the correlation value
In the case of obtaining by equation (9) or the like, C-1> C+1If
A is calculated by interpolation using (11).0And C0, C-1Passing
And C+1And the exact position A of the corresponding point usingrSeeking
Can be.

【0019】[0019]

【数13】 (C−C-1)= (C0 −C-1)/(A0 −A-1) (A−A-1) …式(13)(C−C −1 ) = (C 0 −C −1 ) / (A 0 −A −1 ) (A−A −1 ) Equation (13)

【0020】[0020]

【数14】 (C−C+1)=− (C0 −C-1)/(A0 −A-1) (A−A+1)…式(14)(C−C +1 ) = − (C 0 −C −1 ) / (A 0 −A −1 ) (A−A +1 ) Equation (14)

【0021】[0021]

【数15】 (C-1−C+1)=− (C0 −C-1)/(A0 −A-1) (2A−A-1−A+1) …式(15)(C −1 −C +1 ) = − (C 0 −C −1 ) / (A 0 −A −1 ) (2A−A −1 −A +1 ) Equation (15)

【0022】[0022]

【数16】 Ar =(1/2) (A-1+A+1)−((C-1−C+1)/(C0−C-1))(A0−A-1) =A0 − (1/2)(C-1−C+1)/(C0 −C-1) …式(16)A r = (1/2) (A −1 + A +1 ) − ((C −1 −C +1 ) / (C 0 −C −1 )) (A 0 −A −1 ) = A 0 − (1/2) (C −1 −C +1 ) / (C 0 −C −1 ) Equation (16)

【0023】図5はC-1<C+1となる場合の例であり、
この場合は、(A0 ,C0 )と(A +1,C+1)を通る実
線の直線17を引き、この直線17と符号が反対で絶対
値が等しい傾きを持ち且つ(A-1,C-1)を通る二点鎖
線の直線18との交点を求めれば、そのA座標が対応点
の位置Ar となる。この場合、直線17は次式(17)
で表わされる。従って、直線18は次式(18)で表さ
れる。両直線の交点を求めるために式(17),(1
8)からCを消去すると、次式(19)が得られる。式
(19)中のAは対応点の位置Ar にほかならないか
ら、式(19)を展開すると、次式(20)を得る。こ
の場合も、パターンを切り出す点は1画素ずつずれるの
で、A0 −A-1=1画素、A+1−A0 =1画素、(A-1
+A+1)/2=A0 であることから、式(20)を整理
すると、前記の式(12)が得られる。即ち、相関値を
計算式(9)等で求める場合は、C-1<C+1であれば、
前式(12)を用いた補間演算によりA0 と、C0 ,C
-1及びC+1とを用いて精密な対応点の位置Ar を求める
ことができる。
FIG.-1<C+1This is an example where
In this case, (A0, C0) And (A +1, C+1) Passing through
Draw a straight line 17 of the line, and sign
Values have equal slopes and (A-1, C-1) Two-point chain passing
If the intersection of the line with the straight line 18 is found, the A coordinate is the corresponding point.
Position ArBecomes In this case, the straight line 17 is expressed by the following equation (17).
Is represented by Therefore, the straight line 18 is expressed by the following equation (18).
It is. Equations (17) and (1) are used to determine the intersection of the two straight lines.
Eliminating C from 8) gives the following equation (19). formula
A in (19) is the position A of the corresponding pointrIs it nothing but
Then, when Expression (19) is expanded, the following Expression (20) is obtained. This
In the case of, the point where the pattern is cut out is shifted by one pixel.
And A0-A-1= 1 pixel, A+1-A0= 1 pixel, (A-1
+ A+1) / 2 = A0Equation (20) is rearranged.
Then, the above equation (12) is obtained. That is, the correlation value
In the case of using the calculation formula (9) or the like, C-1<C+1If,
A is calculated by interpolation using the above equation (12).0And C0, C
-1And C+1And the exact position A of the corresponding point usingrAsk for
be able to.

【0024】[0024]

【数17】 (C−C+1)= (C+1−C0)/(A+1−A0) (A−A+1) …式(17)(C−C +1 ) = (C + 1− C 0 ) / (A + 1− A 0 ) (A−A +1 ) Equation (17)

【0025】[0025]

【数18】 (C−C-1)=− (C+1−C0)/(A+1−A0) (A−A-1) …式(18)(C−C −1 ) = − (C + 1− C 0 ) / (A + 1− A 0 ) (A−A −1 ) Equation (18)

【0026】[0026]

【数19】 (C+1−C-1)=− (C+1−C0)/(A+1−A0) (2A−A+1−A-1) …式(19)Equation 19] (C +1 -C -1) = - (C +1 -C 0) / (A +1 -A 0) (2A-A +1 -A -1) ... (19)

【0027】[0027]

【数20】 Ar =(1/2) (A+1+A-1)−((C+1−C-1)/(C+1−C0))(A+1−A0)] =A0 + (1/2)(C+1−C-1)/(C0 −C+1) …式(20)Equation 20] A r = (1/2) (A +1 + A -1) - ((C +1 -C -1) / (C +1 -C 0)) (A +1 -A 0)] = A 0 + (1/2) (C +1 −C −1 ) / (C 0 −C +1 ) Equation (20)

【0028】次に、前記式(10)など、相関が強いほ
ど相関値Cが大きくなる計算式を用いる場合、AとCと
の直交座標系を導入すると、図6または図7に示すよう
に、パターンを切り出した点とその点での相関値との対
(A0 ,C0 ),(A-1,C -1),(A+1,C+1)をプ
ロットすることかできる。但し、A0 は相関が最も強い
第1の点の位置、C0 は第1の点での相関値、A-1は第
1の点を含む一次元A上で第1の点の直前に位置する第
2の点の位置、C-1は第2の点での相関値、A +1は同一
次元A上で第1の点の直後に位置する第3の点の位置、
+1は第3の点での相関値をそれぞれ示している。
Next, when the correlation is strong, as in the above equation (10),
When a calculation formula that increases the correlation value C is used, A and C
When the rectangular coordinate system is introduced, as shown in FIG. 6 or FIG.
The pair of the point at which the pattern was cut out and the correlation value at that point
(A0, C0), (A-1, C -1), (A+1, C+1)
Can be lotted. Where A0Has the strongest correlation
Location of the first point, C0Is the correlation value at the first point, A-1Is the
The first position on the one-dimensional A including the point 1 located immediately before the first point
Position of point 2, C-1Is the correlation value at the second point, A +1Are the same
The position of a third point on the dimension A immediately after the first point;
C+1Indicates the correlation value at the third point.

【0029】図6はC-1<C+1となる場合の例であり、
この場合は、図4と同じく(A0 ,C0 )と(A-1,C
-1)を通る実線の直線15を引き、この直線15と符号
が反対で絶対値が等しい傾きを持ち且つ(A+1,C+1
を通る二点鎖線の直線16との交点を求めれば、そのA
座標が対応点の位置Ar となる。この場合、直線15は
前記式(13)で表わされる。従って、直線16も前記
式(14)で表される。両直線の交点を求めるために式
(13),(14)からCを消去すると、前式(15)
が得られる。式(15)中のAは対応点の位置Ar にほ
かならないから、式(15)を展開すると、前式(1
6)を得る。そして、パターンを切り出す点は1画素ず
つずれるので、A0 −A-1=1画素、A+1−A0 =1画
素、(A-1+A+1)/2=A0 であることから、式(1
6)を整理すると、前記の式(11)が得られる。即
ち、相関値を計算式(10)等で求める場合は、C-1
+1であれば、前式(11)を用いた補間演算によりA
0 と、C0 ,C-1及びC+1とを用いて精密な対応点の位
置Ar を求めることができる。
FIG. 6 shows an example in which C -1 <C +1 .
In this case, (A 0 , C 0 ) and (A −1 , C
-1 ), a solid line 15 is drawn, the sign of which is opposite to that of this straight line 15, the slopes are equal, and (A +1 , C +1 )
Is obtained by finding the intersection of the two-dot chain line with the straight line 16 passing through
Coordinate is a position A r of the corresponding points. In this case, the straight line 15 is represented by the equation (13). Therefore, the straight line 16 is also represented by the equation (14). Eliminating C from equations (13) and (14) to find the intersection of both straight lines yields equation (15)
Is obtained. Since the A in the formula (15) nothing but positions A r of the corresponding points, Expanding equation (15), Equation (1
6) is obtained. Since the point at which the pattern is cut is shifted by one pixel, A 0 −A −1 = 1 pixel, A +1 −A 0 = 1 pixel, and (A −1 + A +1 ) / 2 = A 0 , Equation (1)
By rearranging 6), the above equation (11) is obtained. That is, when the correlation value is obtained by the calculation formula (10) or the like, C -1 <
If C +1 , A is calculated by interpolation using equation (11).
0, it is possible to determine the position A r of precise corresponding points by using the C 0, C -1 and C +1.

【0030】図7はC-1>C+1となる場合の例であり、
この場合は、図5と同じく(A0 ,C0 )と(A+1,C
+1)を通る実線の直線17を引き、この直線17と符号
が反対で絶対値が等しい傾きを持つ且つ(A-1,C-1
を通る二点鎖線の直線16との交点を求めれば、そのA
座標が対応点の位置Ar となる。この場合、直線17は
前記式(17)で表わされる。従って、直線18も前記
式(18)で表される。従って、両直線の交点を求める
ために式(17),(18)からCを消去すると、前式
(19)が得られる。式(19)中のAは対応点の位置
r にほかならないから、式(19)を展開すると、前
式(30)を得る。この場合も、パターンを切り出す点
は1画素ずつずれるので、A0 −A-1=1画素、A+1
0 =1画素、(A-1+A+1)/2=A0 であることか
ら、式(20)を整理すると、前記の式(12)が得ら
れる。即ち、相関値を計算式(10)等で求める場合
は、C-1>C+1であれば、前式(12)を用いた補間演
算によりA0 と、C0 ,C-1及びC+1とを用いて精密な
対応点の位置Ar を求めることができる。
FIG. 7 shows an example in which C -1 > C +1 .
In this case, (A 0 , C 0 ) and (A +1 , C
+1 ), a straight line 17 is drawn, the sign of which is opposite to that of this straight line 17 and whose absolute value is equal to that of (A −1 , C −1 ).
Is obtained by finding the intersection of the two-dot chain line with the straight line 16 passing through
Coordinate is a position A r of the corresponding points. In this case, the straight line 17 is represented by the equation (17). Therefore, the straight line 18 is also represented by the equation (18). Therefore, if C is eliminated from equations (17) and (18) in order to find the intersection of the two straight lines, equation (19) is obtained. Since the A in the formula (19) nothing but positions A r of the corresponding points, obtained by expanding the formula (19), Equation (30). Also in this case, the point at which the pattern is cut is shifted by one pixel, so that A 0 −A −1 = 1 pixel and A +1
Since A 0 = 1 pixel and (A -1 + A +1 ) / 2 = A 0 , the above equation (12) is obtained by rearranging equation (20). That is, when the correlation value is obtained by the calculation formula (10) or the like, if C −1 > C +1 , A 0 , C 0 , C −1, and C 0 are obtained by interpolation using the previous formula (12). it is possible to obtain the position a r of precise corresponding points by using the +1.

【0031】[0031]

【実施例】以下、図面を参照して本発明の実施例を説明
する。図1は相関マッチングの対応点位置計測装置の構
成例を示し、この例の装置は相関値計算部1と、モデル
格納部2と、データ抽出部3と、大小関係判定部4と、
補間式選択部5と、補間式格納部6と、補間演算部7と
から構成されている。
Embodiments of the present invention will be described below with reference to the drawings. FIG. 1 shows an example of a configuration of a corresponding point position measuring device for correlation matching. The device of this example includes a correlation value calculating unit 1, a model storage unit 2, a data extracting unit 3, a magnitude relation determining unit 4,
It is composed of an interpolation formula selection unit 5, an interpolation formula storage unit 6, and an interpolation calculation unit 7.

【0032】相関値計算部1はモデル格納部2に予め記
憶させてある濃度パターンのモデルの情報2aを用い
て、入力画像8の1画素ずつすらした各点毎に、モデル
と同じ形状で且つ同じ大きさのパターンを切り出し、前
出の式(9)や式(10)等の計算式によりモデルと切
り出したパターンとの相関値を計算し、各点での相関値
の計算結果1aをデータ抽出部3に与える。
The correlation value calculating unit 1 uses the density pattern model information 2a stored in the model storage unit 2 in advance, and for each point of the input image 8 which is even one pixel at a time, has the same shape as the model. A pattern having the same size is cut out, and a correlation value between the model and the cut out pattern is calculated according to the calculation formulas (9) and (10) described above, and a calculation result 1a of a correlation value at each point is calculated as data. It is given to the extraction unit 3.

【0033】データ抽出部3は相関値計算部1から与え
られた各点での相関値の計算結果1aから、相関が最も
強い第1の点の位置A0 と、同第1の点での相関値C0
と、第1の点を含む一次元上で第1の点の直前に位置す
る第2の点での相関値C-1と、同一次元上で第1の点の
直後に位置する第3の点での相関値C+1とを求め、これ
らのデータを大小関係判定部4と補間演算部7に与え
る。
The data extracting unit 3 calculates the position A 0 of the first point having the strongest correlation and the position A 0 of the first point from the calculation result 1a of the correlation value at each point given from the correlation value calculating unit 1. Correlation value C 0
And a correlation value C -1 at a second point located immediately before the first point on one dimension including the first point, and a third value located immediately after the first point on the same dimension. A correlation value C +1 at the point is obtained, and these data are given to the magnitude relationship determination unit 4 and the interpolation calculation unit 7.

【0034】大小関係判定部4は本実施例では、データ
抽出部3から与えられる第2の点での相関値C-1と第3
の点での相関値C+1との大小関係を、第1の点での相関
値C 0 との差の絶対値|C-1−C0 |,|C+1−C0
を用いて判定することとし、その判定結果4aを補間式
選択部5に与える。
In the present embodiment, the magnitude relation determining unit 4 determines the data
Correlation value C at the second point given from extraction unit 3-1And the third
Correlation value C at point+1The magnitude relationship with the first point correlation
Value C 0The absolute value of the difference with | C-1-C0|, | C+1-C0|
And the result of determination 4a is calculated using an interpolation formula.
This is given to the selection unit 5.

【0035】補間式選択部5は大小関係判定部4から与
えられた判定結果に応じて、補間式格納部6に予め記憶
させてある前式(11)及び前式(12)の一方を選択
する。具体的には、(1)|C-1−C0 |>|C+1−C
0 |の場合は、式(11)を選択し、(2)|C-1−C
0 |<|C+1−C0 |の場合は、式(12)を選択す
る。
The interpolation formula selection unit 5 selects one of the previous formulas (11) and (12) stored in the interpolation formula storage unit 6 in advance according to the judgment result given from the magnitude relation judgment unit 4. I do. Specifically, (1) | C -1 -C 0 |> | C +1 -C
0 |, equation (11) is selected, and (2) | C −1 −C
0 | <| C +1 -C 0 | For selects equation (12).

【0036】補間演算部7は補間式選択部6により選択
した式(11)または式(12)を用いて、データ抽出
部3から与えられたA0 ,C0 ,C-1及びC+1から、精
密な対応点の位置Ar を求める。
The interpolation operation unit 7 uses the equations (11) or (12) selected by the interpolation equation selection unit 6 to provide A 0 , C 0 , C -1 and C +1 given from the data extraction unit 3. from obtains the position a r of the precision corresponding points.

【0037】なお、|C-1−C0 |=|C+1−C0 |の
場合は、第1の点が真の対応点そのものであるから、補
間計算をする必要なくA0 が対応点の位置となる。しか
し、処理の繁雑さを避けることから、判定条件の一方に
等号を入れて、例えば下記のようにしても良い。(1)
|C-1−C0 |≧|C+1−C0 |の場合は、式(11)
を選択し、(2)|C-1−C0 |<|C+1−C0 |の場
合は、式(12)を選択する。
In the case of | C −1 −C 0 | = | C +1 −C 0 |, since the first point is the true corresponding point, A 0 corresponds without the need for interpolation calculation. The position of the point. However, in order to avoid the complexity of the processing, an equal sign may be included in one of the determination conditions, for example, as described below. (1)
When | C −1 −C 0 | ≧ | C + 1 −C 0 |, the equation (11) is used.
And (2) If | C −1 −C 0 | <| C +1 −C 0 |, the equation (12) is selected.

【0038】ところで、大小関係判定部4で|C-1−C
0 |と|C+1−C0 |との間の大小関係を判定するよう
にしたのは、相関値計算部1で式(9)のように相関が
強いほど相関値が小さくなる計算式と式(10)のよう
に相関が強いほど相関値が大きくなる計算式いずれを用
いても、それに影響されずに適切な補間式(11)また
は(12)を選択できるようにしたためである。即ち、 (1)相関値を計算式(9)で求める場合、|C-1−C
0 |>|C+1−C0 |はC-1>C+1と等価であり式(1
1)を選択し、|C-1−C0 |<|C+1−C0|はC-1
<C+1と等価であり式(12)を選択する。 (2)逆に、相関値を計算式(10)で求める場合は、
|C-1−C0 |>|C +1−C0 |がC-1<C+1と等価で
あり式(11)を選択し、|C-1−C0 |<|C+1−C
0 |がC-1>C+1と等価であり式(12)を選択する。
By the way, | C-1-C
0| And | C+1-C0To determine the magnitude relationship between |
The reason is that the correlation is calculated by the correlation value calculation unit 1 as shown in Expression (9).
Equation (10), where the correlation value decreases as the value increases.
Whichever formula is used, the stronger the correlation, the larger the correlation value.
Even if it is, the appropriate interpolation formula (11) or
Is because (12) can be selected. That is, (1) When calculating the correlation value by the calculation formula (9), | C-1-C
0| > | C+1-C0| Is C-1> C+1Equation (1)
1) Select | C-1-C0| <| C+1-C0| Is C-1
<C+1Equation (12) is selected. (2) On the contrary, when the correlation value is calculated by the formula (10),
| C-1-C0| > | C +1-C0| Is C-1<C+1Is equivalent to
Expression (11) is selected and | C-1-C0| <| C+1-C
0| Is C-1> C+1Equation (12) is selected.

【0039】もちろん、相関値計算部1が計算式(9)
など相関が強いほど相関値が小さくなる計算式しか使わ
ない場合は、単に大小関係判定部4でC-1とC+1の大小
関係を判定し、補間式選択部5ではC-1≧C+1の場合に
式(11)を選択し、C-1<C+1の場合に式(12)を
選択するようにしても良い。
Of course, the correlation value calculator 1 calculates the equation (9).
For example, when using only a calculation formula in which the correlation value becomes smaller as the correlation becomes stronger, the magnitude relationship between C −1 and C +1 is simply determined by the magnitude relationship determination unit 4, and the interpolation formula selection unit 5 determines C −1 ≧ C select the expression (11) in the case of +1 may be selected equation (12) in the case of C -1 <C +1.

【0040】同様に、相関値計算部1が計算式(10)
など相関が強いほど相関値が大きくなる計算式しか使わ
ない場合は、単に大小関係判定部4でC-1とC+1の大小
関係を判定し、補間式選択部5ではC-1≦C+1の場合に
式(11)を選択し、C-1>C+1の場合に式(12)を
選択するようにしても良い。
Similarly, the correlation value calculation unit 1 calculates the equation (10)
If, for example, only a calculation formula that increases the correlation value as the correlation becomes stronger is used, the magnitude relationship between C −1 and C +1 is simply determined by the magnitude relationship determination unit 4, and the interpolation formula selection unit 5 determines C −1 ≦ C select the expression (11) in the case of +1 may be selected equation (12) in the case of C -1> C +1.

【0041】図2に、図1に示した実施例装置の動作の
流れを示す。
FIG. 2 shows an operation flow of the embodiment apparatus shown in FIG.

【0042】<二次元、三次元等への拡張>上述した説
明は一次元における対応点の位置計測についてである
が、二次元または三次元等における対応点の位置を計測
するには次のようにすれば良い。 (1)二次元の場合:この場合は図1に示した装置のう
ち、二点鎖線で囲む部分9を2つ用意し、XY直交座標
系を考えて、一方ではAをX、A0 をX0 、A r をXr
と読み替え、各点での相関値の計算結果から、相関が最
も強い第1の点位置X0 と、同第1の点での相関値C0
と、第1の点を含むX軸上で第1の点の直前に位置する
第2の点での相関値C-1と、同X軸上で第1の点の直後
に位置する第3の点での相関値C+1とを求め、第2の点
での相関値C-1と第3の点での相関値C+1との大小関係
を判定し、判定結果に応じて式(11)または式(1
2)を選択して、入力画像におけるモデルとの対応点の
X軸上の位置Xr を求める。また、他方ではAをY、A
0 をY0 、Ar をYr と読み替え、各点での相関値の計
算結果から、相関が最も強い第1の点の位置Y0 と、同
第1の点での相関値C0 と、第1の点を含むY軸上で第
1の点の直前に位置する第2の点での相関値C-1と、同
Y軸上で第1の点の直後に位置する第3の点での相関値
+1を求め、第2の点での相関値C-1と第3の点での相
関値C+1との大小関係を判定し、判定結果に応じて式
(11)または式(12)を選択して、入力画像におけ
るモデルとの対応点のY軸上の位置Yr を求める。これ
により、(Xr ,Yr )が二次元での対応点の位置とし
て求まる。 (2)三次元の場合:この場合は二点鎖線で囲んだ部分
9を更に1つ追加し、AをZ、A0 をZ0 、Ar をZr
と読み替え、各点での相関値の計算結果から、相関が最
も強い第1の点の位置Z0 と、同第1の点での相関値C
0 と、第1の点を含むZ軸上で第1の点の直前に位置す
る第2の点での相関値C-1と、同Z軸上で第1の点の直
後に位置する第3の点での相関値C+1を求め、第2の点
での相関値C-1と第3の点での相関値C+1との大小関係
を判定し、判定結果に応じて式(11)または式(1
2)を選択し、入力画像におけるモデルとの対応点のZ
軸上の位置Zr を求める。これにより、(Xr ,Yr
r )が三次元での対応点の位置として求まる。 (3)上記2つの例では並列処理であるが、二点鎖線で
囲んだ部分9は1つだけとし、順次Xr ,Yr ,Zr
求めるようにしても良い。
<Extension to 2D, 3D, etc.>
Akira is about measuring the position of corresponding points in one dimension
Measures the position of corresponding points in 2D or 3D
To do so, do the following. (1) Two-dimensional case: In this case, the apparatus shown in FIG.
Then, two portions 9 surrounded by a two-dot chain line are prepared, and XY rectangular coordinates are prepared.
Considering the system, on the other hand, A is X, A0To X0, A rTo Xr
From the calculation result of the correlation value at each point.
Strong first point position X0And the correlation value C at the first point0
And on the X-axis including the first point, immediately before the first point
Correlation value C at the second point-1And immediately after the first point on the same X axis
The correlation value C at the third point located at+1And the second point
Correlation value C at-1And the correlation value C at the third point+1Relationship with
Is determined, and equation (11) or equation (1) is determined according to the determination result.
2) Select the corresponding point of the input image
Position X on X axisrAsk for. On the other hand, A is Y, A
0Is Y0, ArIs YrAnd calculate the correlation value at each point.
From the calculation result, the position Y of the first point having the strongest correlation0And the same
Correlation value C at the first point0And on the Y axis including the first point
Correlation value C at a second point located immediately before point 1-1And the same
Correlation value at a third point located immediately after the first point on the Y axis
C+1, And the correlation value C at the second point-1And the phase in the third point
Function value C+1Judge the magnitude relationship with
Select (11) or (12) to select
Y on the Y-axis of the corresponding point with the modelrAsk for. this
Gives (Xr, Yr) Is the position of the corresponding point in two dimensions
Is determined. (2) Three-dimensional case: In this case, the part surrounded by the two-dot chain line
9 is added, A is Z, A0To Z0, ArTo Zr
From the calculation result of the correlation value at each point.
Strong first point position Z0And the correlation value C at the first point
0And on the Z axis including the first point, immediately before the first point.
The correlation value C at the second point-1And the first point on the Z axis
Correlation value C at the third point located later+1Find the second point
Correlation value C at-1And the correlation value C at the third point+1Relationship with
Is determined, and equation (11) or equation (1) is determined according to the determination result.
2) is selected and Z of the corresponding point with the model in the input image is selected.
Position Z on axisrAsk for. Thus, (Xr, Yr,
Zr) Is obtained as the position of the corresponding point in three dimensions. (3) In the above two examples, parallel processing is used.
Enclosed part 9 is only one, and Xr, Yr, ZrTo
You may ask for it.

【0043】[0043]

【発明の効果】以上の説明から判るように、従来はモデ
ルとの対応点の位置を画像の分解能に相当する精度でし
か計測できなかったが、本発明によれば補間演算を行う
ことにより、画像の分解能よりも精密に対応点の位置を
計測することができる。また、補間演算に用いる式は比
較的単純であるから、演算に要する時間が短いという利
点がある。
As can be understood from the above description, the position of the corresponding point with the model can be measured only with an accuracy equivalent to the resolution of the image in the past. The position of the corresponding point can be measured more precisely than the resolution of the image. Further, since the equation used for the interpolation calculation is relatively simple, there is an advantage that the time required for the calculation is short.

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

【図1】本発明の一実施例の装置構成を示す図。FIG. 1 is a diagram showing an apparatus configuration according to an embodiment of the present invention.

【図2】同装置の動作の流れを示す図。FIG. 2 is a diagram showing a flow of operation of the apparatus.

【図3】本発明の原理となる相関値変化の対称性を示す
図。
FIG. 3 is a diagram showing a symmetry of a correlation value change which is a principle of the present invention.

【図4】補間式導出の説明図。FIG. 4 is an explanatory diagram of derivation of an interpolation expression.

【図5】補間式導出の説明図。FIG. 5 is an explanatory diagram of derivation of an interpolation formula.

【図6】補間式導出の説明図。FIG. 6 is an explanatory diagram of derivation of an interpolation formula.

【図7】補間式導出の説明図。FIG. 7 is an explanatory diagram of derivation of an interpolation equation.

【図8】相関マッチング処理及び対応点の説明図。FIG. 8 is an explanatory diagram of a correlation matching process and corresponding points.

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

1 相関値計算部 2 モデル格納部 3 データ抽出部 4 大小関係判定部 5 補間式選択部 6 補間式格納部 7 補間演算部 8 入力画像 Reference Signs List 1 correlation value calculation unit 2 model storage unit 3 data extraction unit 4 magnitude relation determination unit 5 interpolation formula selection unit 6 interpolation formula storage unit 7 interpolation calculation unit 8 input image

───────────────────────────────────────────────────── フロントページの続き (58)調査した分野(Int.Cl.7,DB名) G01B 11/00 - 11/30 102 G06T 1/00 G06T 7/00 G06T 7/60 ──────────────────────────────────────────────────続 き Continuation of front page (58) Fields investigated (Int. Cl. 7 , DB name) G01B 11/00-11/30 102 G06T 1/00 G06T 7/00 G06T 7/60

Claims (4)

(57)【特許請求の範囲】(57) [Claims] 【請求項1】 入力画像の各点毎に、モデルと同じ形状
で且つ同じ大きさのパターンを切り出し、モデルと切り
出したパターンとの相関値を相関が強いほど相関値が小
さくなる計算式で計算すること; 各点での相関値の計算結果から、相関が最も強い第1の
点の位置A0 と、同第1の点での相関値C0 と、第1の
点を含む一次元上で第1の点の直前に位置する第2の点
での相関値C-1と、同一次元上で第1の点の直後に位置
する第3の点での相関値C+1とを求めること; 第2の点での相関値C-1と第3の点での相関値C+1との
大小関係を判定すること; 判定結果から、C-1>C+1の場合は下式(1)により、
-1<C+1の場合は下式(2)により、入力画像におけ
るモデルとの対応点の位置Ar を補間演算で求めるこ
と;により、画像分解能よりも精密に対応点の位置を計測す
ること を特徴とする相関マッチングの対応点位置計測方
法。 【数1】 Ar =A0 −(1/2)(C-1−C+1)/(C0 −C-1) …式(1) 【数2】 Ar =A0 +(1/2)(C+1−C-1)/(C0 −C+1) …式(2)
1. A pattern having the same shape and the same size as a model is cut out for each point of an input image, and a correlation value between the model and the cut out pattern is calculated using a calculation formula in which the stronger the correlation, the smaller the correlation value. From the calculation result of the correlation value at each point, the position A 0 of the first point having the strongest correlation, the correlation value C 0 at the first point, and the one-dimensional position including the first point To calculate a correlation value C −1 at a second point located immediately before the first point and a correlation value C +1 at a third point located immediately after the first point on the same dimension. Determining the magnitude relationship between the correlation value C −1 at the second point and the correlation value C +1 at the third point; from the determination result, if C −1 > C +1 , By (1),
In the case of C -1 <C + 1 , the position Ar of the corresponding point with the model in the input image is obtained by interpolation by the following equation (2); thereby, the position of the corresponding point is measured more precisely than the image resolution. You
The corresponding point position measuring method of the correlation matching characterized by Rukoto. A r = A 0- (1/2) (C -1 -C +1 ) / (C 0 -C -1 ) Equation (1) Ar = A 0 + (1) / 2) (C +1 −C −1 ) / (C 0 −C +1 ) Equation (2)
【請求項2】 入力画像の各点毎に、モデルと同じ形状
で且つ同じ大きさのパターンを切り出し、モデルと切り
出したパターンとの相関値を相関が強いほど相関値が大
きくなる計算式で計算すること; 各点での相関値の計算結果から、相関が最も強い第1の
点の位置A0 と、同第1の点での相関値C0 と、第1の
点を含む一次元上で第1の点の直前に位置する第2の点
での相関値C-1と、同一次元上で第1の点の直後に位置
する第3の点での相関値C+1とを求めること; 第2の点での相関値C-1と第3の点での相関値C+1との
大小関係を判定すること; 判定結果から、C-1<C+1の場合は下式(3)により、
-1>C+1の場合は下式(4)により、入力画像におけ
るモデルとの対応点の位置Ar を補間演算で求めるこ
と;により、画像分解能よりも精密に対応点の位置を計測す
ること を特徴とする相関マッチングの対応点位置計測方
法。 【数3】 Ar =A0 −(1/2)(C-1−C+1)/(C0 −C-1) …式(3) 【数4】 Ar =A0 +(1/2)(C+1−C-1)/(C0 −C+1) …式(4)
2. A pattern having the same shape and the same size as the model is cut out for each point of the input image, and the correlation value between the model and the cut out pattern is calculated by a calculation formula in which the stronger the correlation, the larger the correlation value. From the calculation result of the correlation value at each point, the position A 0 of the first point having the strongest correlation, the correlation value C 0 at the first point, and the one-dimensional position including the first point To calculate a correlation value C −1 at a second point located immediately before the first point and a correlation value C +1 at a third point located immediately after the first point on the same dimension. Determining the magnitude relationship between the correlation value C -1 at the second point and the correlation value C +1 at the third point; from the determination result, if C -1 <C +1 , (3)
In the case of C -1 > C + 1 , the position Ar of the corresponding point with the model in the input image is obtained by interpolation by the following equation (4); thereby, the position of the corresponding point is measured more precisely than the image resolution. You
The corresponding point position measuring method of the correlation matching characterized by Rukoto. Equation 3] A r = A 0 - (1/2 ) (C -1 -C +1) / (C 0 -C -1) ... Equation (3) [Expression 4] A r = A 0 + (1 / 2) (C +1 −C −1 ) / (C 0 −C +1 ) Equation (4)
【請求項3】 入力画像の各点毎に、モデルと同じ形状
で且つ同じ大きさのパターンを切り出し、モデルと切り
出したパターンとの相関値を相関が強いほど相関値が小
さくなる計算式で計算する手段と; 各点での相関値の計算結果から、相関が最も強い第1の
点の位置A0 と、同第1の点での相関値C0 と、第1の
点を含む一次元上で第1の点の直前に位置する第2の点
での相関値C-1と、同一次元上で第1の点の直後に位置
する第3の点での相関値C+1とを求める手段と; 下記2つの式(5)及び式(6)を格納している手段
と; 第2の点での相関値C-1と第3の点での相関値C+1との
大小関係を判定する手段と; 判定結果から、C-1>C+1の場合は下式(5)を選択
し、C-1<C+1の場合は下式(6)を選択する手段と; 選択された式により、入力画像におけるモデルとの対応
点の位置Ar を補間演算で求める手段と; を備え、画像分解能よりも精密に対応点の位置を計測す
ることを特徴とする相関マッチングの対応点位置計測装
置。 【数5】 Ar =A0 −(1/2)(C-1−C+1)/(C0 −C-1) …式(5) 【数6】 Ar =A0 +(1/2)(C+1−C-1)/(C0 −C+1) …式(6)
3. A pattern having the same shape and the same size as the model is cut out for each point of the input image, and the correlation value between the model and the cut out pattern is calculated by a calculation formula in which the stronger the correlation, the smaller the correlation value becomes. A position A 0 of the first point having the strongest correlation, a correlation value C 0 at the first point, and a one-dimensional pattern including the first point based on the calculation result of the correlation value at each point. A correlation value C -1 at a second point located immediately before the first point above and a correlation value C +1 at a third point located immediately after the first point on the same dimension are calculated as follows. Means for obtaining; means for storing the following two equations (5) and (6); magnitude of correlation value C -1 at the second point and correlation value C +1 at the third point Means for judging the relationship; means for selecting the following equation (5) if C -1 > C +1 and selecting the following equation (6) if C -1 <C +1 from the judgment results: The selected expression More, means for determining the interpolation calculation the position A r of the corresponding point of the model in the input image; equipped with, than the image resolution to measure the position of the precisely corresponding points
The corresponding point position measuring device of the correlation matching characterized by Rukoto. A r = A 0- (1/2) (C -1 -C +1 ) / (C 0 -C -1 ) Equation (5) Ar = A 0 + (1) / 2) (C +1 -C -1 ) / (C 0 -C +1 ) Equation (6)
【請求項4】 入力画像の各点毎に、モデルと同じ形状
で且つ同じ大きさのパターンを切り出し、モデルと切り
出したパターンとの相関値を相関が強いほど相関値が大
きくなる計算式で計算する手段と; 各点での相関値の計算結果から、相関が最も強い第1の
点の位置A0 と、同第1の点での相関値C0 と、第1の
点を含む一次元上で第1の点の直前に位置する第2の点
での相関値C-1と、同一次元上で第1の点の直後に位置
する第3の点での相関値C+1とを求める手段と; 下記2つの式(7)及び式(8)を格納している手段
と; 第2の点での相関値C-1と第3の点での相関値C+1との
大小関係を判定する手段と; 判定結果から、C-1<C+1の場合は下式(7)を選択
し、C-1>C+1の場合は下式(8)を選択する手段と; 選択された式により、入力画像におけるモデルとの対応
点の位置Ar を補間演算で求める手段と; を備え、画像分解能よりも精密に対応点の位置を計測す
ることを特徴とする相関マッチングの対応点位置計測装
置。 【数7】 Ar =A0 −(1/2)(C-1−C+1)/(C0 −C-1) …式(7) 【数8】 Ar =A0 +(1/2)(C+1−C-1)/(C0 −C+1) …式(8)
4. A pattern having the same shape and the same size as the model is cut out for each point of the input image, and the correlation value between the model and the cut out pattern is calculated by a calculation formula in which the stronger the correlation, the larger the correlation value. A position A 0 of the first point having the strongest correlation, a correlation value C 0 at the first point, and a one-dimensional pattern including the first point based on the calculation result of the correlation value at each point. A correlation value C -1 at a second point located immediately before the first point above and a correlation value C +1 at a third point located immediately after the first point on the same dimension are calculated as follows. Means for obtaining; means for storing the following two equations (7) and (8); magnitude of the correlation value C -1 at the second point and the correlation value C +1 at the third point Means for judging the relationship; means for selecting the following equation (7) if C -1 <C +1 and selecting the following equation (8) if C -1 > C +1 from the judgment results: The selected expression More, means for determining the interpolation calculation the position A r of the corresponding point of the model in the input image; equipped with, than the image resolution to measure the position of the precisely corresponding points
The corresponding point position measuring device of the correlation matching characterized by Rukoto. A r = A 0 − (1 /) (C −1 −C +1 ) / (C 0 −C −1 ) Equation (7) Ar = A 0 + (1) / 2) (C +1 −C −1 ) / (C 0 −C +1 ) Equation (8)
JP17802294A 1994-07-29 1994-07-29 Corresponding point position measurement method and apparatus for correlation matching Expired - Fee Related JP3296099B2 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP17802294A JP3296099B2 (en) 1994-07-29 1994-07-29 Corresponding point position measurement method and apparatus for correlation matching

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP17802294A JP3296099B2 (en) 1994-07-29 1994-07-29 Corresponding point position measurement method and apparatus for correlation matching

Publications (2)

Publication Number Publication Date
JPH0843024A JPH0843024A (en) 1996-02-16
JP3296099B2 true JP3296099B2 (en) 2002-06-24

Family

ID=16041217

Family Applications (1)

Application Number Title Priority Date Filing Date
JP17802294A Expired - Fee Related JP3296099B2 (en) 1994-07-29 1994-07-29 Corresponding point position measurement method and apparatus for correlation matching

Country Status (1)

Country Link
JP (1) JP3296099B2 (en)

Also Published As

Publication number Publication date
JPH0843024A (en) 1996-02-16

Similar Documents

Publication Publication Date Title
US7746377B2 (en) Three-dimensional image display apparatus and method
Zhu et al. Generalized parallel-perspective stereo mosaics from airborne video
US8208029B2 (en) Method and system for calibrating camera with rectification homography of imaged parallelogram
US20190311497A1 (en) Methods and computer program products for calibrating stereo imaging systems by using a planar mirror
EP4411627A1 (en) Photogrammetry method, apparatus and device, and storage medium
US6993179B1 (en) Strapdown system for three-dimensional reconstruction
CN112288826B (en) Calibration method and device of binocular camera and terminal
JP2002517045A (en) How to determine the exact position of the template match point
JPH09231373A (en) Device for measuring three-dimensional position
JP3075023B2 (en) 3D shape data processing method
JPH0773344A (en) Method and apparatus for three- dimensional point in two-dimensional graphic display
JP3282516B2 (en) 3D image measurement method
JP2000171214A (en) Corresponding point retrieving method and three- dimensional position measuring method utilizing same
US11385053B2 (en) Distance measuring camera
JP3296099B2 (en) Corresponding point position measurement method and apparatus for correlation matching
JP2802034B2 (en) 3D object measurement method
JPH06195472A (en) Image-processing system
JP3912638B2 (en) 3D image processing device
JP4444440B2 (en) Stereo image measuring device
JP3277105B2 (en) Method and apparatus for creating partial solid model
KR100914857B1 (en) Method and apparatus for detecting of workpiece position and matching of CAD data using 3 dimensional vision camera system
JP4153322B2 (en) Method and apparatus for associating measurement points in photogrammetry
CN115409693A (en) Two-dimensional positioning method based on pipeline foreign matters in three-dimensional image
JP2001124519A (en) Recurring correspondent point survey method, three- dimensional position measuring method thereof, these devices and recording medium
JP3385783B2 (en) 3D position measurement method in stereo measurement

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: 20020312

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

Free format text: PAYMENT UNTIL: 20080412

Year of fee payment: 6

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

Free format text: PAYMENT UNTIL: 20090412

Year of fee payment: 7

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

Free format text: PAYMENT UNTIL: 20090412

Year of fee payment: 7

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

Free format text: PAYMENT UNTIL: 20100412

Year of fee payment: 8

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

Free format text: PAYMENT UNTIL: 20100412

Year of fee payment: 8

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

Free format text: PAYMENT UNTIL: 20110412

Year of fee payment: 9

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

Free format text: PAYMENT UNTIL: 20130412

Year of fee payment: 11

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

Free format text: PAYMENT UNTIL: 20140412

Year of fee payment: 12

LAPS Cancellation because of no payment of annual fees