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
JP5034852B2 - Image measuring method and image measuring apparatus - Google Patents
[go: Go Back, main page]

JP5034852B2 - Image measuring method and image measuring apparatus - Google Patents

Image measuring method and image measuring apparatus Download PDF

Info

Publication number
JP5034852B2
JP5034852B2 JP2007263306A JP2007263306A JP5034852B2 JP 5034852 B2 JP5034852 B2 JP 5034852B2 JP 2007263306 A JP2007263306 A JP 2007263306A JP 2007263306 A JP2007263306 A JP 2007263306A JP 5034852 B2 JP5034852 B2 JP 5034852B2
Authority
JP
Japan
Prior art keywords
image
template
calculated
template image
pixel
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
JP2007263306A
Other languages
Japanese (ja)
Other versions
JP2009093412A (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.)
Fujitsu Ltd
Original Assignee
Fujitsu 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 Fujitsu Ltd filed Critical Fujitsu Ltd
Priority to JP2007263306A priority Critical patent/JP5034852B2/en
Publication of JP2009093412A publication Critical patent/JP2009093412A/en
Application granted granted Critical
Publication of JP5034852B2 publication Critical patent/JP5034852B2/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

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

Description

本発明はFA(Factory Automation)分野の製造ラインなどでの部品・マークなどの位置、形状の計測に係り、特に極小の対象物を撮影し、対象物の基準テンプレートと撮影画像とのマッチングにより対象物の位置・形状の計測を精度良く行なう画像計測方法及び画像計測装置に関するものである。   The present invention relates to the measurement of the position and shape of parts, marks, and the like on a production line in the FA (Factor Automation) field, and in particular, it captures a very small object and matches it with a reference template of the object and a captured image. The present invention relates to an image measurement method and an image measurement apparatus that accurately measure the position and shape of an object.

FA分野の製造ラインなどでは、ロボットなどによりライン上の部品を掴み取り上げて、装置に装着あるいはプリント板の指定位置に挿入する。あるいはマークの形状を把握し、マークの方向を合わせて対象物の指定位置に貼り付けることが行なわれている。しかしながら、作今の装置などの小型化の進展に伴ない、対象の部品、マークなどは非常に小さく、また、複雑な形状の部品も多くなっている。このため対象物の位置・形状を探索し計測する必要がある。方策として、対象物を撮影した画像(探索画像)を対象物の基準テンプレート(誤解し易い場合はテンプレート画像と表現)とマッチング(テンプレートマッチング)させ、形状・位置を計測して目的の装置への装着、プリント板への挿入が行なわれる。   In a production line in the FA field, a part on the line is picked up by a robot or the like and is mounted on the apparatus or inserted into a designated position on the printed board. Alternatively, the shape of the mark is grasped, and the direction of the mark is aligned and pasted to the designated position of the object. However, along with the progress of miniaturization of devices such as devices, the target parts and marks are very small, and the number of parts with complicated shapes is also increasing. For this reason, it is necessary to search and measure the position and shape of the object. As a measure, the image (search image) obtained by photographing the object is matched (template matching) with the reference template of the object (if it is easily misunderstood, it is expressed as a template image), and the shape / position is measured to obtain the target device. Mounting and insertion into the printed board are performed.

図10は従来のテンプレートマッチングによるサブピクセル位置推定方法(その1)を示す図である。従来のテンプレートマッチングにより探索物を画像の基本単位である画素(ピクセル)以下の精度で計測を行なう方法を示している。なお、このピクセル以下の精度での計測をサブピクセル位置推定と表現する。1)テンプレート、2)探索画像とマッチング領域、3)フィッティング関数とピーク値を示している。
ア.1)で示すテンプレートを2)の探索画像上の左上を起点として、1ピクセル単位のピッチで位置をずらしながら(図の矢印の方向)探索画像とテンプレート画像とのマッチングを行ない、ピクセル単位の位置での相関値を算出する。2つの図形間の相関値の算出については一般に知られているので説明は省略する。
イ.ピクセル単位の座標位置を引数として算出した相関値を近似するフィッティング関数を算出する。フィッティング関数がピーク値となる相関ピーク位置を推定する(ここでは説明上1次元としXの関数として2次曲線で表す)。この例ではX座標の(319)と(320)の間にピークがあり、このX座標位置「x」より探索物の位置を推定する。位置を推定することを「探索物を計測する」と表現する。
FIG. 10 is a diagram illustrating a conventional subpixel position estimation method (No. 1) based on template matching. This shows a method of measuring a search object with a precision equal to or lower than a pixel (pixel) which is a basic unit of an image by conventional template matching. In addition, the measurement with the accuracy below this pixel is expressed as sub-pixel position estimation. 1) Template, 2) Search image and matching area, 3) Fitting function and peak value.
A. The template shown in 1) starts from the upper left of the search image in 2), and the search image and the template image are matched while shifting the position at a pitch of 1 pixel unit (in the direction of the arrow in the figure). The correlation value at is calculated. Since the calculation of the correlation value between two figures is generally known, a description thereof will be omitted.
A. A fitting function that approximates the correlation value calculated using the coordinate position of the pixel unit as an argument is calculated. The correlation peak position at which the fitting function becomes the peak value is estimated (here, it is assumed to be one-dimensional for the sake of explanation and represented by a quadratic curve as a function of X). In this example, there is a peak between X coordinates (319) and (320), and the position of the search object is estimated from the X coordinate position “x 0 ”. Estimating the position is expressed as “measuring a search object”.

しかしながら、テンプレート画像と探索画像が同一(サブピクセル移動量がゼロ)の場合においても、「±1」ピクセルのサンプリング位置における探索画像のマッチング領域が完全には一致しないため、「±1」ピクセルの位置での相関値は僅かに異なったものとなる。その結果、フィッティング関数近似によって得られたピーク位置がゼロにはならず、ずれてしまう。任意の探索画像での計測位置にもこのようなずれが加わることになり、サブピクセル以下の精度での探索物計測の誤差の要因となっている。   However, even when the template image and the search image are the same (subpixel movement amount is zero), the search image matching region at the sampling position of “± 1” pixels does not completely match, so The correlation values at the positions are slightly different. As a result, the peak position obtained by fitting function approximation does not become zero, but shifts. Such a shift is also added to a measurement position in an arbitrary search image, which causes a search object measurement error with sub-pixel accuracy or less.

図11は従来のテンプレートマッチングによるサブピクセル位置推定方法(その2)を示す図である。テンプレートとして対象物の一部を撮影したテンプレート画像を採用し、図10で述べたテンプレートマッチングを実行する。1)部分テンプレート、2)部分テンプレートによるマッチングを示している。   FIG. 11 is a diagram showing a conventional subpixel position estimation method (No. 2) based on template matching. A template image obtained by photographing a part of an object is adopted as a template, and the template matching described in FIG. 10 is executed. 1) Partial template, 2) Matching by partial template is shown.

図11の1)で示す探索物の一部のテンプレートと2)で示す探索画像とのマッチングを行ない、フィッティング関数を算出し、相関ピーク位置よりX座標の位置を推定する。   A part of the template of the search object shown in 1) of FIG. 11 and the search image shown in 2) are matched, a fitting function is calculated, and the position of the X coordinate is estimated from the correlation peak position.

この部分マッチングにより、相関値を算出した場合、各々のマッチング領域の画像の差異が大きいため、この誤差は更に強く現れる。この誤差のため製造ラインでの部品実装ができなくなる、あるいはマーク位置がずれるなどの問題が発生する。   When the correlation value is calculated by this partial matching, this error appears even more strongly because there is a large difference between the images in each matching region. Due to this error, problems such as the inability to mount components on the production line or the shift of the mark position occur.

上述したように、従来のテンプレートマッチングによるサブピクセル位置の推定では、テンプレート画像の内部構造(輝度分布)によらず、常に同一のフィッティング関数によりピーク位置推定を行っていた。しかし、実際にはテンプレート画像の輝度分布に応じて相関値の分布形状は微妙に異なるため、フィッティング関数近似によるピーク推定位置にずれが生じ、サブピクセルの計測誤差となっていた。   As described above, in the estimation of the subpixel position by the conventional template matching, the peak position is always estimated by the same fitting function regardless of the internal structure (luminance distribution) of the template image. However, in actuality, the distribution shape of the correlation value is slightly different depending on the luminance distribution of the template image, so that the peak estimation position is shifted due to the fitting function approximation, resulting in a subpixel measurement error.

特許文献1では複数のマッチング基本画像の各々に対しマッチング対象画像の位置をずらしてマッチングを行ない、各マッチングピークが得られた複数のマッチング度とマッチング位置よりマッチング位置偏差を補正して偏差を除いたマッチング位置及びマッチング度を算出する技術が開示されている。   In Patent Document 1, matching is performed by shifting the position of the matching target image with respect to each of a plurality of matching basic images, and the matching position deviation is corrected based on a plurality of matching degrees and matching positions at which each matching peak is obtained, and the deviation is removed. A technique for calculating a matching position and a matching degree is disclosed.

特許文献2では計測対象物のモデルの撮影画像上で、モデルの各々異なる特徴を示す複数のモデル画像を生成し、各々最大となる相関値を算出し対象画像の中心位置を算出する技術が開示されている。
特開平09−81729号公報 特開平10−153407号公報
Patent Document 2 discloses a technique for generating a plurality of model images showing different features of a model on a captured image of a model of a measurement target, calculating a maximum correlation value, and calculating a center position of the target image. Has been.
JP 09-81729 A JP-A-10-153407

解決しようとする課題は、FA分野の製造ラインで極小部品、マークなどの探索対象物を撮影した探索画像とテンプレート画像とのマッチングにより探索対象物の位置・形状を計測するテンプレートマッチングにおいて、撮影画像の基本単位の1ピクセル以下(サブピクセル)の計測ができない問題である。   The problem to be solved is that, in template matching in which the position and shape of a search object are measured by matching a search image obtained by shooting a search object such as a minimal part or a mark on a manufacturing line in the FA field with a template image, This is a problem that it is impossible to measure 1 pixel or less (subpixel) of the basic unit.

本発明は、テンプレートマッチングによる探索画像の計測において、ピクセル以下のサブピクセルの精度で探索対象物の位置・形状計測を実現する画像計測方法及び装置を提供することを目的とする。   An object of the present invention is to provide an image measurement method and apparatus that realizes the position / shape measurement of a search object with subpixel accuracy below a pixel in measurement of a search image by template matching.

第1の発明は、探索対象物の標準テンプレート画像と前記探索対象物を撮影した探索画像のテンプレートマッチングによる相関値の分布をフィッティングすることにより前記探索対象物の位置・形状を計測するテンプレートマッチングによる画像計測方法である。 A first aspect of the present invention is by template matching to measure the position and shape of the search object by fitting the distribution of the correlation value by template matching of the search image taken the search object as the standard template image in the search object This is an image measurement method.

前記標準テンプレート画像より前記探索対象物の領域を切り出したテンプレート画像の原点を前記標準テンプレート画像の原点と一致させた基準位置より、前記テンプレート画像を前記基準位置よりサブピクセル移動した位置と、前記サブピクセル移動位置より+1ピクセル移動した位置と前記サブピクセル移動位置より−1ピクセル移動した位置と、で前記テンプレート画像と前記標準テンプレート画像との自己相関値を各々算出し、前記算出した3つの自己相関値より曲線近似した関数を算出し、前記算出した近似関数の自己相関値が最大となる位置を算出し、前記算出した自己相関値最大位置とサブピクセル移動位置との差を補正量として算出し、前記サブピクセル移動位置を複数の異なる位置に移動させた位置で前記補正量を算出し、前記算出した前記複数補正量の内、前記補正量が一番小さいサブピクセル移動位置の補正量で前記探索対象物の位置を補正して前記探索対象物の位置・形状を計測する。 A position obtained by moving the template image by a sub pixel from the reference position from a reference position in which the origin of the template image obtained by cutting out the region of the search object from the standard template image is coincident with the origin of the standard template image; Autocorrelation values of the template image and the standard template image are respectively calculated from a position moved by +1 pixel from the pixel movement position and a position moved by -1 pixel from the subpixel movement position, and the calculated three autocorrelations A function that approximates a curve from the value is calculated, a position where the autocorrelation value of the calculated approximate function is maximum is calculated, and a difference between the calculated maximum position of the autocorrelation value and the subpixel movement position is calculated as a correction amount. The correction amount is calculated at a position where the sub-pixel movement position is moved to a plurality of different positions. Among the plurality correction amount the calculated, the correction amount is to measure the position and shape of the search target position compensation by the said search object in the correction amount of the smallest sub-pixel movement position.

第1の発明により、テンプレート画像が有する誤差を複数のサブピクセル移動位置毎の補正量の内、探索対象物の位置に一番近い移動位置の補正量で補正することができ、探索対象物の位置・形状をサブピクセルの精度で計測できる。 According to the first invention, the error of the template image can be corrected with the correction amount of the movement position closest to the position of the search target object among the correction amounts for the plurality of sub-pixel movement positions . The position and shape can be measured with subpixel accuracy.

第2の発明は、第1の発明において前記探索物の位置を補正する補正量は、切り出す前記探索対象物の領域を複数の領域で切り出し、前記複数の切り出した領域のテンプレート画像より前記補正量を算出し、前記算出した複数の補正量の内、最小となる領域での補正量で前記探索対象物の位置を補正する。 In the second invention, the correction amount for correcting the position of the search object in the first invention is that the region of the search object to be cut out is cut out in a plurality of regions, and the correction amount is obtained from a template image of the plurality of cut out regions. And the position of the search object is corrected with the correction amount in the minimum region among the plurality of calculated correction amounts.

第2の発明により、テンプレート画像が有する誤差を複数の領域で切り出した領域のテンプレート画像から算出した複数の補正量の内、最小となる領域での補正量を算出することができ、探索対象物の位置・形状の計測を精度良くできる。 According to the second invention, it is possible to calculate the correction amount in the minimum region among the plurality of correction amounts calculated from the template image of the region obtained by cutting out the error of the template image in the plurality of regions, and the search object The position and shape of the can be accurately measured.

第3の発明は、第1の発明において、相関値の分布によるフィッティングは前記テンプレート画像を複数位置に移動して移動位置毎に自己相関値を算出し、前記各移動位置と該移動位置での自己相関値をフィッティング関数で近似して算出する。 In a third aspect based on the first aspect, the fitting based on the correlation value distribution is performed by moving the template image to a plurality of positions and calculating an autocorrelation value for each movement position. The autocorrelation value is approximated with a fitting function.

第3の発明により、相関値の分布をフィッティング関数近似することにより、探索対象物の位置・形状の計測を精度良くできる。 According to the third aspect of the present invention, the position / shape of the search target can be measured with high accuracy by approximating the distribution of correlation values to the fitting function.

第4の発明は、探索対象物の標準テンプレート画像と前記探索対象物を撮影した探索画像のテンプレートマッチングによる相関値の分布をフィッティングすることにより前記探索対象物の位置・形状を計測するテンプレートマッチングによる画像計測装置である。
前記標準テンプレート画像より前記探索対象物の領域を切り出したテンプレート画像の原点を前記標準テンプレート画像の原点と一致させた基準位置より、前記テンプレート画像を前記基準位置よりサブピクセル移動した位置と、前記サブピクセル移動位置より+1ピクセル移動した位置と前記サブピクセル移動位置より−1ピクセル移動した位置と、で前記テンプレート画像と前記標準テンプレート画像との自己相関値を各々算出する手段と、
前記算出した3つの自己相関値より曲線近似した関数を算出する手段と、前記算出した近似関数の自己相関値が最大となる位置を算出する手段と、前記算出した自己相関値最大位置とサブピクセル移動位置との差を補正量として算出する手段と、前記サブピクセル移動位置を複数の異なる位置に移動させた位置で前記補正量を算出する手段と、
前記算出した前記複数補正量の内、前記補正量が一番小さいサブピクセル移動位置の補正量で前記探索対象物の位置を補正して前記探索対象物の位置・形状を計測する手段と、
を有する。
The fourth invention is based on template matching that measures the position and shape of the search object by fitting a correlation value distribution by template matching between the standard template image of the search object and the search image obtained by photographing the search object. This is an image measurement device.
A position obtained by moving the template image by a sub pixel from the reference position from a reference position in which the origin of the template image obtained by cutting out the region of the search object from the standard template image is coincident with the origin of the standard template image; Means for calculating autocorrelation values of the template image and the standard template image respectively at a position moved by +1 pixel from the pixel movement position and a position moved by -1 pixel from the sub-pixel movement position;
Means for calculating a function that approximates a curve from the calculated three autocorrelation values; means for calculating a position at which the autocorrelation value of the calculated approximation function is maximum; and the calculated autocorrelation value maximum position and subpixel Means for calculating a difference from a movement position as a correction amount; means for calculating the correction amount at a position obtained by moving the sub-pixel movement position to a plurality of different positions;
Means for measuring the position / shape of the search object by correcting the position of the search object with the correction amount of the sub-pixel movement position with the smallest correction amount among the calculated plurality of correction amounts ;
Have

第4の発明により、テンプレート画像が有する誤差を複数のサブピクセル移動位置毎の補正量の内、探索対象物の位置に一番近い移動位置の補正量で補正することができ、探索対象物の位置・形状をサブピクセルの精度で計測できる画像計測装置を実現できる。 According to the fourth invention, the error of the template image can be corrected with the correction amount of the movement position closest to the position of the search target object among the correction amounts for the plurality of sub-pixel movement positions. An image measuring device that can measure the position and shape with sub-pixel accuracy can be realized.

本発明により、個々のテンプレート画像が有している計測誤差の要因に併せてその誤差量を補正量として算出し、探索画像の計測結果を補正することにより、サブピクセルの精度で探索物の位置・形状を計測することが可能であり、計測誤差を個別に補正して計測精度を向上させることができる。   According to the present invention, the position of the search object is calculated with subpixel accuracy by calculating the amount of error as a correction amount in accordance with the measurement error factor of each template image and correcting the measurement result of the search image. -The shape can be measured, and measurement errors can be individually corrected to improve measurement accuracy.

(実施例1)
図1は本発明の一実施形態の画像位置計測装置ブロック構成を示す図である。
Example 1
FIG. 1 is a diagram showing a block configuration of an image position measuring apparatus according to an embodiment of the present invention.

画像位置計測装置1は画像入力部10、テンプレート画像処理部20、テンプレートマッチング処理部30、計測結果出力部40で構成する。以下に各部の機能処理内容を述べる。
1)画像入力部10
探索対象物の画像(以下、標準テンプレート画像と表現)、標準テンプレート画像から切り出し探索画像とのマッチングを行なうために登録する画像(以下、テンプレート画像と表現)、探索画像を入力し、登録部20、テンプレートマッチング処理部30に送る。
2)テンプレート画像処理部20
画像登録部21、補正量算出部22で構成する。
3)画像登録部21
画像入力部10から標準テンプレート画像を受信し、テンプレート画像として登録する。
4)補正量算出部22
テンプレート画像よりテンプレート自身の誤差を算出し、探索画像の計測結果を補正する補正量として保持する。テンプレート自身の持つフィッティング誤差の算出(補正量の算出)方法については後述する。(図2、図5)
5)テンプレートマッチング処理部30
相関演算部31、サブピクセル推定部32、推定位置補正部33で構成し、画像入力部10から探索画像、テンプレート画像処理部20よりテンプレート画像、補正量を受信し、探索画像の計測と補正を行い、計測結果を計測結果出力部40に出力する。
6)相関演算部31
探索画像に対し、テンプレート画像を移動させ、移動した各位置での相関演算を行なう。移動位置とその位置で相関値より、フィッティング関数で近似する処理を行なう。フィッティング関数の最大値とその位置を算出する。
7)サブピクセル位置推定部32
相関演算部31で算出したフィッティング関数の相関最大値とその位置より探索画像の位置をサブピクセル精度で推定する。
8)推定位置補正部33
サブピクセル推定部32での探索画像の位置推定結果とテンプレート画像処理部20からのテンプレート画像の誤差の補正量より推定した探索物の位置を補正し探索物の位置計測を行なう。
9)計測結果出力部40
テンプレートマッチング処理部30からの探索物の計測結果を出力する。
The image position measuring apparatus 1 includes an image input unit 10, a template image processing unit 20, a template matching processing unit 30, and a measurement result output unit 40. The function processing contents of each part are described below.
1) Image input unit 10
An image of a search object (hereinafter referred to as a standard template image), an image to be registered for matching with a search image cut out from the standard template image (hereinafter referred to as a template image), and a search image are input, and a registration unit 20 And sent to the template matching processing unit 30.
2) Template image processing unit 20
The image registration unit 21 and the correction amount calculation unit 22 are included.
3) Image registration unit 21
A standard template image is received from the image input unit 10 and registered as a template image.
4) Correction amount calculation unit 22
The error of the template itself is calculated from the template image, and is stored as a correction amount for correcting the measurement result of the search image. A method for calculating the fitting error of the template itself (calculation of correction amount) will be described later. (Fig. 2, Fig. 5)
5) Template matching processing unit 30
A correlation calculation unit 31, a sub-pixel estimation unit 32, and an estimated position correction unit 33 are configured to receive a search image from the image input unit 10, a template image and a correction amount from the template image processing unit 20, and measure and correct the search image. The measurement result is output to the measurement result output unit 40.
6) Correlation calculation unit 31
The template image is moved with respect to the search image, and correlation calculation is performed at each moved position. A process of approximating with a fitting function is performed from the moving position and the correlation value at that position. The maximum value of the fitting function and its position are calculated.
7) Subpixel position estimation unit 32
The position of the search image is estimated with sub-pixel accuracy from the maximum correlation value of the fitting function calculated by the correlation calculation unit 31 and its position.
8) Estimated position correction unit 33
The position of the search object is estimated by correcting the position of the search object estimated from the position image estimation result of the subpixel estimation unit 32 and the template image error correction amount from the template image processing unit 20.
9) Measurement result output unit 40
The measurement result of the search object from the template matching processing unit 30 is output.

図2はテンプレート画像の補正量算出手順(その1)を示す図である。
S1:テンプレート用画像(標準テンプレート画像)より対象物の領域を切り出し、テンプレート画像として登録する。
S2:標準テンプレート画像上に設定した基準位置でのテンプレート画像と標準テンプレート画像との自己相関値を算出する。
S3:テンプレート画像を標準テンプレート画像上で基準位置より「−1」及び「+1」ピクセル移動し、その位置での標準テンプレート画像との自己相関値を算出する。
FIG. 2 is a diagram showing a template image correction amount calculation procedure (part 1).
S1: A region of an object is cut out from a template image (standard template image) and registered as a template image.
S2: An autocorrelation value between the template image at the reference position set on the standard template image and the standard template image is calculated.
S3: The template image is moved by “−1” and “+1” pixels from the reference position on the standard template image, and an autocorrelation value with the standard template image at that position is calculated.

なお、基準位置とはテンプレート画像と標準テンプレート画像の各原点を一致させ(同一画像に相当)、ピクセル単位では画像が移動していない(サブピクセル移動量ゼロ)位置である。
S4:算出した3つの位置とその相関値より2次曲線で近似したフィッティング関数を算出する。
S5:フィッティング関数の自己相関値が最大となる位置を算出し、その自己相関最大位置と基準位置との差をテンプレートの補正量「H」とする。
Note that the reference position is a position where the origins of the template image and the standard template image coincide with each other (corresponding to the same image), and the image does not move in pixel units (zero sub-pixel movement amount).
S4: A fitting function approximated by a quadratic curve is calculated from the calculated three positions and their correlation values.
S5: A position where the autocorrelation value of the fitting function is maximized is calculated, and a difference between the maximum autocorrelation position and the reference position is set as a template correction amount “H”.

S4で3点の相関値から2次曲線の算出と最大値の算出について述べる。ここでは1次元とし位置xに対し、自己相関値をyとし、y=f(x)=ax+bx+c=a(x+b2/4a)+(c−b/4a)とする。[x]の位置と[x]より[−1]、[+1]ピクセル移動させた位置[x0−1]、[x0+1]とその自己相関値、(x0−1,y)、(x,y)、(x0+1,y)、より係数「a」、「b」、「c」は各々、a=(y+y−2*y)/2、b=(y−y)−a*(2*x0−1)と算出でき、自己相関値が最大値となる「xmax」はxmax=(−b/2a)となる。 In S4, calculation of the quadratic curve and calculation of the maximum value from the three correlation values will be described. Here, the one-dimensional position is assumed to be y with respect to the position x, and y = f (x) = ax 2 + bx + c = a (x + b 2 / 4a) 2 + (c−b 2 / 4a). To do. The position of the [x 0] [x 0] from [-1], [+ 1] positions were pixels moved [x 0-1], [x 0 + 1] and its autocorrelation values, (x 0-1 , Y 1 ), (x 0 , y 2 ), (x 0 + 1 , y 3 ), and the coefficients “a”, “b”, “c” are respectively a = (y 1 + y 3 −2 * y 2 ) / 2, b = (y 2 −y 1 ) −a * (2 * x 0−1 ), and “x max ” at which the autocorrelation value becomes the maximum value is x max = (− b / 2a).

なお、ここではフィッティング関数として二次曲線による近似について説明したが、テンプレートマッチングによる自己相関値の分布がより複雑となり、2次曲線では近似できないと判断した場合は、自己相関値分布に適した近似関数で行なえば良い。ここでは説明を省略する。また、後述する(図2、図3、図5、図8)でのフィッティング関数近似も同様に2次曲線での近似として説明する。   Although the approximation by the quadratic curve has been described here as the fitting function, if the autocorrelation value distribution by template matching becomes more complicated and cannot be approximated by the quadratic curve, the approximation suitable for the autocorrelation value distribution This can be done with a function. The description is omitted here. The fitting function approximation described later (FIGS. 2, 3, 5, and 8) is also described as an approximation with a quadratic curve.

また、自己相関によるテンプレート自身の誤差算出として、対象物の標準画像からテンプレート画像の撮影を繰り返し撮影した平均画像をテンプレート画像として用いることによって、ランダムノイズの影響を抑えることもできる。詳細な説明は省略する。   In addition, as an error calculation of the template itself due to autocorrelation, the influence of random noise can be suppressed by using an average image obtained by repeatedly capturing a template image from a standard image of an object as the template image. Detailed description is omitted.

図3は探索画像のサブピクセル位置推定手順(その1)を示す図である。
S11:計測する対象物の画像(探索画像)を撮影する。
S12:撮影した探索画像とテンプレート画像との相関値をテンプレート初期位置「E0」より順次所定位置移動して各位置で算出する。
S13:算出した相関値が最大となった位置とその位置より「−1」及び「+1」ピクセル移動した3点での相関値を2次曲線で近似したフィッティング関数を算出する。
S14:フィッティング関数の相関値が最大となるサブピクセル位置を算出し、そのサブピクセル位置を補正量「H」で補正し、探索画像の推定位置とする。
FIG. 3 is a diagram showing a subpixel position estimation procedure (part 1) of the search image.
S11: An image (search image) of an object to be measured is taken.
S12: The correlation value between the photographed search image and the template image is sequentially moved from the template initial position “E0” by a predetermined position and calculated at each position.
S13: A fitting function is calculated by approximating a correlation value at a point where the calculated correlation value is maximized and three points shifted by “−1” and “+1” pixels from the position with a quadratic curve.
S14: A sub-pixel position where the correlation value of the fitting function is maximized is calculated, and the sub-pixel position is corrected with the correction amount “H” to obtain an estimated position of the search image.

なお、S13、S14でのフィッティング関数の算出、フィッティング関数の相関値が最大となる位置の算出は図2の最後に述べた手順と同一の方法で行なう。   The calculation of the fitting function in S13 and S14 and the calculation of the position where the correlation value of the fitting function is maximized are performed by the same method as the procedure described at the end of FIG.

図4はテンプレート誤差の補正による探索画像位置推定を説明する図である。テンプレート自身の持つフィッティング誤差を算出し、算出した誤差より、探索画像の位置計測結果の補正について説明している。1)テンプレートの誤差、2)テンプレート誤差より探索画像の計測誤差補正を示している。
図4のアはテンプレートの自己相関値から算出したテンプレート自身の誤差であり、イはアの誤差を補正するテンプレート誤差補正量である(逆ベクトルとなる)。ウで示す探索画像の補正前の推定位置をイの誤差補正量で補正し、エで示す補正後の推定位置を得る。
(実施例2)
実施例1では、テンプレート自身の誤差の補正のために自己相関マッチングによる補正量算出として、サブピクセル移動の無い同一テンプレート画像から得た補正量を一律に利用する方法について説明した。それに対し、ここでは、サブピクセルの移動を仮定したテンプレート画像を生成し、それぞれとのマッチングを行って補正量を細かく求め、サブピクセル推定位置に応じて適切な補正を加える方法について説明する。
FIG. 4 is a diagram for explaining search image position estimation by correcting the template error. The fitting error of the template itself is calculated, and correction of the position measurement result of the search image is described based on the calculated error. 1) Template error, 2) Measurement error correction of search image based on template error.
4A is an error of the template itself calculated from the autocorrelation value of the template, and A is a template error correction amount for correcting the error of (a). The estimated position before correction of the search image indicated by C is corrected with the error correction amount of A to obtain the corrected estimated position indicated by D.
(Example 2)
In the first embodiment, the method of uniformly using the correction amount obtained from the same template image without sub-pixel movement as the correction amount calculation by autocorrelation matching for correcting the error of the template itself has been described. On the other hand, here, a description will be given of a method of generating a template image on the assumption that the subpixel is moved, performing matching with each of the template images, finely obtaining a correction amount, and applying an appropriate correction according to the subpixel estimated position.

図5はテンプレート画像の補正量算出手順(その2)を示す図である。
S21:テンプレート用画像(標準テンプレート画像)より対象物の領域を切り出し、テンプレート画像として登録する。
S22:テンプレート画像を標準テンプレート画像上の基準位置より1ピクセル以下(サブピクセル)移動させた位置「Ri」(i=0)に設定する。
S23:標準テンプレート画像とテンプレート画像の移動位置での相関値を算出する。
S24:テンプレート画像を標準テンプレート画像上でさらに「−1」及び「+1」ピクセル移動し、その2つの位置で標準テンプレート画像との相関値を算出する。
S25:設定した3つの位置とその位置での相関値を2次曲線で近似したフィッティング関数を算出する。
S26:フィッティング関数の相関値が最大となる位置「Ris」を算出し、「Ri」位置での補正値「Hri」(Hri=Ris−Ri)を算出する。
S27:i=i+1とする。
S28:設定されている次の移動位置があるか確認し、ある場合はS29に進み、そうで無い場合は終了となる。
S29:次のテンプレートのサブピクセル移動位置「Ri」を設定し、S22に戻る。
FIG. 5 is a diagram showing a template image correction amount calculation procedure (part 2).
S21: The object region is cut out from the template image (standard template image) and registered as a template image.
S22: The template image is set to a position “Ri” (i = 0) obtained by moving the template image by one pixel or less (subpixel) from the reference position on the standard template image.
S23: A correlation value at the movement position of the standard template image and the template image is calculated.
S24: The template image is further moved by “−1” and “+1” pixels on the standard template image, and a correlation value with the standard template image is calculated at the two positions.
S25: A fitting function is calculated by approximating the set three positions and the correlation values at the positions with a quadratic curve.
S26: A position “Ris” at which the correlation value of the fitting function is maximized is calculated, and a correction value “Hri” (Hri = Ris−Ri) at the “Ri” position is calculated.
S27: i = i + 1.
S28: It is confirmed whether there is a next set movement position. If there is, the process proceeds to S29, and if not, the process ends.
S29: The sub-pixel movement position “Ri” of the next template is set, and the process returns to S22.

なお、S25、26のフィッティング関数の算出、フィッティング関数の相関値が最大となる位置の算出は前述した。   The calculation of the fitting function in S25 and S26 and the calculation of the position where the correlation value of the fitting function is maximized have been described above.

図6はテンプレート領域の選択例(その1)を示す図である。図5で算出した標準テンプレート画像に対するテンプレート画像の位置より1つのテンプレート領域を選択する例を示している。1)テンプレート領域、2)自己相関による推定位置を示している。   FIG. 6 is a diagram showing a template area selection example (No. 1). FIG. 5 shows an example in which one template region is selected from the position of the template image with respect to the standard template image calculated in FIG. 1) A template area, 2) an estimated position by autocorrelation.

図5のテンプレート画像の位置「Ri」として1)の「A」、「B」、「C」の3点でのテンプレートマッチングによる誤差を算出し、誤差が一番小さい位置(ここでは「B」の位置)例えば、「B」の位置でのテンプレート画像を選択する。   As a position “Ri” of the template image in FIG. 5, an error due to template matching at three points “A”, “B”, and “C” of 1) is calculated, and the position where the error is the smallest (here “B”). For example, the template image at the position “B” is selected.

図7はテンプレート領域の選択例(その2)を示す図である。テンプレート画像がテンプレート群として機能する場合のテンプレート領域の選択方法を示している。部分テンプレートとして、部分テンプレート1、部分テンプレート2、部分テンプレート3を採り上げる(図の部分テンプレートとして括弧で囲まれた画像が部分テンプレート画像である)。なお、部分テンプレートは部分テンプレート画像を移動させる場合、複数の部分テンプレートが連動して移動する。   FIG. 7 is a diagram showing a template area selection example (No. 2). A method for selecting a template area when a template image functions as a template group is shown. As the partial template, the partial template 1, the partial template 2, and the partial template 3 are taken up (the image enclosed in parentheses as the partial template in the figure is the partial template image). Note that when a partial template image is moved, a plurality of partial templates are moved in conjunction with each other.

図6の例と同様に、複数の部分テンプレートを連動して移動させ、部分テンプレートの個数や領域を変えながら図5で述べた自己相関による計測位置をそれぞれ推定し、その推定位置のベクトル和を求め、そのベクトル和の大きさが小さい部分テンプレートの組み合わせ(個数と領域)をテンプレート群として選択する。   As in the example of FIG. 6, a plurality of partial templates are moved in conjunction with each other, the measurement positions based on the autocorrelation described in FIG. 5 are estimated while changing the number and area of the partial templates, and the vector sum of the estimated positions is calculated. Then, a combination (number and region) of partial templates having a small vector sum is selected as a template group.

図8は探索画像のサブピクセル位置推定手順(その2)を示す図である。
S31:計測する対象物の画像(探索画像)を撮影する。
S32:撮影した探索画像とテンプレート画像との相関値をテンプレート初期位置「E0」より順次所定位置移動して各位置で算出する。
S33:算出した相関値が最大となった位置とその位置より「−1」及び「+1」ピクセル移動した3点での相関値を2次曲線で近似したフィッティング関数を算出する。
S34:フィッティング関数の相関値が最大となるサブピクセル位置を算出し、そのサブピクセル位置に一番近い図5で算出したサブピクセル移動位置(Ri)の補正量「Hri」で補正し、探索画像の推定位置とする。
FIG. 8 is a diagram showing a sub-pixel position estimation procedure (part 2) of the search image.
S31: An image of an object to be measured (search image) is taken.
S32: The correlation value between the photographed search image and the template image is sequentially moved from the template initial position “E0” by a predetermined position and calculated at each position.
S33: A fitting function is calculated by approximating the correlation value at the three points shifted by “−1” and “+1” pixels from the position where the calculated correlation value is maximized, and a quadratic curve.
S34: The subpixel position where the correlation value of the fitting function is maximized is calculated, corrected with the correction amount “Hri” of the subpixel movement position (Ri) calculated in FIG. 5 closest to the subpixel position, and the search image Estimated position.

なお、S33のフィッティング関数の算出、フィッティング関数の相関値が最大となる位置の算出は前述した。   The calculation of the fitting function in S33 and the calculation of the position where the correlation value of the fitting function is maximized have been described above.

図9はサブピクセル位置に対応した補正を示す図である。図5で述べたテンプレート画像の位置として、「1/8」ピクセル移動した位置で算出したサブピクセル位置とその補正量を示している。   FIG. 9 is a diagram showing the correction corresponding to the sub-pixel position. As the position of the template image described in FIG. 5, the subpixel position calculated at the position moved by “1/8” pixel and the correction amount are shown.

図9中の「0.0」の位置が同一画像(サブピクセル移動量ゼロ)における位置補正量に相当する。実施例1の方法は、この一定の補正量を一律に加えるやり方に相当する。しかし、実際にはピーク推定位置のずれは一定ではなく、補正すべき量が点線のような変化を示すことが多い。従って、この例のように細かなサブピクセル(ここでは「1/8」ピクセル)移動したテンプレート画像を生成し、その画像とのテンプレートマッチングを各サブピクセル点で位置推定を行う。これにより、探索画像の位置補正をきめ細かく適切に補正することが可能となる。
(付記1)
テンプレート画像と探索対象物を撮影した探索画像のテンプレートマッチングによる相関値の分布をフィッティングすることにより前記探索対象物の位置・形状を計測するテンプレートマッチングによる画像計測方法であって、
前記計測対象物のテンプレート画像同士の自己相関値を算出し、
前記算出した自己相関値より前記テンプレート画像の位置を推定して前記テンプレート画像の誤差を算出し、
前記誤差より前記テンプレート画像の誤差を補正する補正量を算出し、
前記探索画像と前記テンプレート画像とのマッチングにより前記探索対象物の位置を算出し、
前記算出した探索対象物の位置を前記補正量で補正して前記探索対象物の位置・形状を計測することを特徴とするテンプレートマッチングによる画像計測方法。
(付記2)
付記1記載における前記探索物の位置を補正する補正量は、前記テンプレートを複数の位置に移動して複数の補正量を算出し、
前記算出した複数の補正量の内、最小となる前記移動位置での前記テンプレート画像とその補正量で前記探索対象物の位置を補正することを特徴とする付記1記載の画像計測方法。
(付記3)
付記1記載の画像計測方法において、前記テンプレート画像をサブピクセル量移動した位置に移動して自己相関値を算出し、
前記自己相関値より前記テンプレートの前記移動位置での誤差として補正量を算出し、
前記サブピクセル移動した移動位置で算出する補正量を複数のサブピクセル移動位置で算出し、
前記探索画像と前記テンプレート画像とのマッチングにより算出した位置を前記探索画像の位置に一番近い前記サブピクセル移動位置の補正量で補正して前記探索対象物の位置・形状を計測することを特徴とする付記1記載の画像計測方法。
(付記4)
付記1記載における相関値の分布によるフィッティングは前記テンプレート画像を複数位置に移動して移動位置毎に自己相関値を算出し、
前記各移動位置と該移動位置での自己相関値をフィッティング関数で近似して算出することを特徴とする付記1記載の画像計測方法。
(付記5)
テンプレート画像と探索対象物を撮影した探索画像のテンプレートマッチングによる相関値の分布をフィッティングすることにより前記探索対象物の位置・形状を計測するテンプレートマッチングによる画像計測装置であって、
前記画像計測装置は、
前記計測対象物のテンプレート画像同士の自己相関値を算出する手段と、
前記算出した自己相関値より前記テンプレート画像の位置を推定する手段と、
前記推定した位置より前記テンプレート画像の誤差を算出して前記テンプレート画像の誤差を補正する補正量を算出する手段と、
前記探索画像と前記テンプレート画像とのテンプレートマッチングにより前記探索対象物の位置を算出する手段と、で構成し、
前記算出した探索対象物の位置を前記補正量で補正して前記探索対象物の位置・形状を計測することを特徴とするテンプレートマッチングによる画像計測装置。
(付記6)
付記1記載における前記探索物の位置を補正する補正量は、前記テンプレート画像を複数撮影し、前記複数撮影の画像を各々テンプレート画像とし、
前記複数のテンプレート画像の自己相関値から各々補正量を算出し、
前記算出した複数補正量の平均値を前記テンプレート画像の補正量とすることを特徴とする付記1記載の画像計測方法。
(付記7)
付記3記載における前記サブピクセル移動した移動位置で算出する補正量は、
前記移動位置と前記移動位置より±1ピクセル移動した位置での前記テンプレート画像の自己相関値が最大となる位置を算出し、
前記最大となる位置より前記テンプレートの誤差補正量を算出し、
前記誤差補正量を前記移動位置での誤差補正量とすることを特徴とする付記3記載の画像計測方法。
The position of “0.0” in FIG. 9 corresponds to the position correction amount in the same image (subpixel movement amount zero). The method of Embodiment 1 corresponds to a method of uniformly adding this constant correction amount. However, in practice, the deviation of the peak estimation position is not constant, and the amount to be corrected often shows a change like a dotted line. Therefore, as shown in this example, a template image that is moved with fine sub-pixels (here, “1/8” pixel) is generated, and template matching with the image is performed at each sub-pixel point. This makes it possible to finely and appropriately correct the position correction of the search image.
(Appendix 1)
An image measuring method by template matching that measures the position and shape of the search object by fitting a distribution of correlation values by template matching of the search image obtained by photographing the template image and the search object,
Calculate an autocorrelation value between the template images of the measurement object,
An error of the template image is calculated by estimating the position of the template image from the calculated autocorrelation value,
A correction amount for correcting the error of the template image is calculated from the error,
Calculating the position of the search object by matching the search image and the template image;
An image measurement method based on template matching, wherein the position and shape of the search object are measured by correcting the calculated position of the search object with the correction amount.
(Appendix 2)
The correction amount for correcting the position of the search object according to appendix 1 is to calculate a plurality of correction amounts by moving the template to a plurality of positions,
The image measurement method according to claim 1, wherein the position of the search target is corrected with the template image at the minimum moving position and the correction amount among the calculated plurality of correction amounts.
(Appendix 3)
In the image measurement method according to appendix 1, the autocorrelation value is calculated by moving the template image to a position moved by a sub-pixel amount,
From the autocorrelation value, calculate a correction amount as an error at the movement position of the template,
A correction amount calculated at the movement position moved by the sub-pixel is calculated at a plurality of sub-pixel movement positions,
The position / shape of the search object is measured by correcting the position calculated by matching the search image and the template image with the correction amount of the sub-pixel movement position closest to the position of the search image. The image measurement method according to appendix 1.
(Appendix 4)
The fitting based on the distribution of the correlation value according to appendix 1 calculates the autocorrelation value for each movement position by moving the template image to a plurality of positions.
The image measuring method according to claim 1, wherein the moving position and the autocorrelation value at the moving position are calculated by approximation with a fitting function.
(Appendix 5)
An image measuring device by template matching that measures the position and shape of the search object by fitting a distribution of correlation values by template matching of the search image obtained by photographing the template image and the search object,
The image measuring device is
Means for calculating an autocorrelation value between the template images of the measurement object;
Means for estimating the position of the template image from the calculated autocorrelation value;
Means for calculating a correction amount for correcting an error of the template image by calculating an error of the template image from the estimated position;
Means for calculating the position of the search object by template matching between the search image and the template image,
An image measuring apparatus based on template matching, wherein the position and shape of the search object are measured by correcting the calculated position of the search object with the correction amount.
(Appendix 6)
The correction amount for correcting the position of the search object according to appendix 1 is that a plurality of the template images are photographed, and the plurality of photographed images are each used as a template image.
Calculating a correction amount from the autocorrelation values of the plurality of template images,
The image measuring method according to claim 1, wherein the average value of the calculated plurality of correction amounts is used as the correction amount of the template image.
(Appendix 7)
The correction amount calculated at the movement position moved by the sub-pixel according to Supplementary Note 3 is:
Calculating a position at which the autocorrelation value of the template image is maximized at the movement position and a position moved by ± 1 pixel from the movement position;
Calculate an error correction amount of the template from the maximum position,
The image measurement method according to claim 3, wherein the error correction amount is an error correction amount at the moving position.

図1は本発明の一実施形態の画像位置計測装置ブロック構成を示す図である。FIG. 1 is a diagram showing a block configuration of an image position measuring apparatus according to an embodiment of the present invention. 図2はテンプレート画像の補正量算出手順(その1)を示す図である。FIG. 2 is a diagram showing a template image correction amount calculation procedure (part 1). 図3は探索画像のサブピクセル位置推定手順(その1)を示す図である。FIG. 3 is a diagram showing a subpixel position estimation procedure (part 1) of the search image. 図4はテンプレート誤差の補正による探索画像位置推定を説明する図である。FIG. 4 is a diagram for explaining search image position estimation by correcting the template error. 図5はテンプレート画像の補正量算出手順(その2)を示す図である。FIG. 5 is a diagram showing a template image correction amount calculation procedure (part 2). 図6テンプレート領域の選択例(その1)を示す図である。6 is a diagram showing a template region selection example (part 1). 図7テンプレート領域の選択例(その2)を示す図である。7 is a diagram showing a template area selection example (part 2). 図8は探索画像のサブピクセル位置推定手順(その2)を示す図である。FIG. 8 is a diagram showing a sub-pixel position estimation procedure (part 2) of the search image. 図9はサブピクセル位置に対応した補正を示す図である。FIG. 9 is a diagram showing the correction corresponding to the sub-pixel position. 図10は従来のテンプレートマッチングによるサブピクセル位置推定方法(その1)を示す図である。FIG. 10 is a diagram illustrating a conventional subpixel position estimation method (No. 1) based on template matching. 図11は従来のテンプレートマッチングによるサブピクセル位置推定方法(その2)を示す図である。FIG. 11 is a diagram showing a conventional subpixel position estimation method (No. 2) based on template matching.

符号の説明Explanation of symbols

1 画像位置計測装置
10 画像入力部
20 テンプレート画像処理部
21 画像登録部
22 補正量算出部
30 テンプレートマッチング処理部
31 相関演算部
32 サブピクセル推定部
33 推定位置補正部
40 計測結果出力部
DESCRIPTION OF SYMBOLS 1 Image position measuring device 10 Image input part 20 Template image process part 21 Image registration part 22 Correction amount calculation part 30 Template matching process part 31 Correlation calculation part 32 Subpixel estimation part 33 Estimated position correction part 40 Measurement result output part

Claims (4)

探索対象物の標準テンプレート画像と前記探索対象物を撮影した探索画像のテンプレートマッチングによる相関値の分布をフィッティングすることにより前記探索対象物の位置・形状を計測するテンプレートマッチングによる画像計測方法であって、
前記標準テンプレート画像より前記探索対象物の領域を切り出したテンプレート画像の原点を前記標準テンプレート画像の原点と一致させた基準位置より、前記テンプレート画像を前記基準位置よりサブピクセル移動した位置と、前記サブピクセル移動位置より+1ピクセル移動した位置と前記サブピクセル移動位置より−1ピクセル移動した位置と、で前記テンプレート画像と前記標準テンプレート画像との自己相関値を各々算出し、
前記算出した3つの自己相関値より曲線近似した関数を算出し、
前記算出した近似関数の自己相関値が最大となる位置を算出し、
前記算出した自己相関値最大位置とサブピクセル移動位置との差を補正量として算出し、
前記サブピクセル移動位置を複数の異なる位置に移動させた位置で前記補正量を算出し、
前記算出した前記複数補正量の内、前記補正量が一番小さいサブピクセル移動位置の補正量で前記探索対象物の位置を補正して前記探索対象物の位置・形状を計測することを特徴とするテンプレートマッチングによる画像計測方法。
An image measuring method according to template matching for measuring the position and shape of the search object by fitting the distribution of the correlation values by the template matching of the standard template image and the search image taken the search object of the search object ,
A position obtained by moving the template image by a sub pixel from the reference position from a reference position in which the origin of the template image obtained by cutting out the region of the search object from the standard template image is coincident with the origin of the standard template image; Autocorrelation values of the template image and the standard template image are respectively calculated from a position moved by +1 pixel from the pixel movement position and a position moved by -1 pixel from the sub-pixel movement position;
Calculate a function that approximates a curve from the calculated three autocorrelation values,
Calculate the position where the autocorrelation value of the calculated approximate function is maximum,
Calculating the difference between the calculated autocorrelation value maximum position and the subpixel movement position as a correction amount;
Calculating the correction amount at a position where the sub-pixel movement position is moved to a plurality of different positions;
Among the plurality correction amount the calculated, characterized in that to measure the position and shape of the correction amount is the search object correct for position of the search object with the correction amount of the smallest sub-pixel movement position An image measurement method using template matching.
請求項1記載の切り出す前記探索対象物の領域を複数の領域で切り出し、
前記複数の切り出した領域のテンプレート画像より前記補正量を算出し、
前記算出した複数の補正量の内、最小となる領域での補正量で前記探索対象物の位置を補正することを特徴とする請求項1記載の画像計測方法。
The region of the search object to be cut out according to claim 1 is cut out in a plurality of regions,
Calculating the correction amount from a template image of the plurality of clipped regions;
The image measurement method according to claim 1, wherein the position of the search target is corrected with a correction amount in a minimum region among the plurality of calculated correction amounts.
請求項1記載における相関値の分布によるフィッティングは前記テンプレート画像を複数位置に移動して移動位置毎に自己相関値を算出し、前記各移動位置と該移動位置での自己相関値をフィッティング関数で近似して算出することを特徴とする請求項1記載の画像計測方法。 The fitting by the correlation value distribution according to claim 1 is performed by moving the template image to a plurality of positions, calculating an autocorrelation value for each movement position, and calculating the autocorrelation value at each movement position and the movement position using a fitting function. 2. The image measurement method according to claim 1, wherein the calculation is performed by approximation . 探索対象物の標準テンプレート画像と前記探索対象物を撮影した探索画像のテンプレートマッチングによる相関値の分布をフィッティングすることにより前記探索対象物の位置・形状を計測するテンプレートマッチングによる画像計測装置であって、
前記標準テンプレート画像より前記探索対象物の領域を切り出したテンプレート画像の原点を前記標準テンプレート画像の原点と一致させた基準位置より、前記テンプレート画像を前記基準位置よりサブピクセル移動した位置と、前記サブピクセル移動位置より+1ピクセル移動した位置と前記サブピクセル移動位置より−1ピクセル移動した位置と、で前記テンプレート画像と前記標準テンプレート画像との自己相関値を各々算出する手段と、
前記算出した3つの自己相関値より曲線近似した関数を算出する手段と、
前記算出した近似関数の自己相関値が最大となる位置を算出する手段と、
前記算出した自己相関値最大位置とサブピクセル移動位置との差を補正量として算出する手段と、
前記サブピクセル移動位置を複数の異なる位置に移動させた位置で前記補正量を算出する手段と、
前記算出した前記複数補正量の内、前記補正量が一番小さいサブピクセル移動位置の補正量で前記探索対象物の位置を補正して前記探索対象物の位置・形状を計測する手段と、
を有する事を特徴とするテンプレートマッチングによる画像計測装置。
An image measuring device by template matching that measures the position and shape of the search object by fitting a distribution of correlation values by template matching between a standard template image of the search object and a search image obtained by photographing the search object. ,
A position obtained by moving the template image by a sub pixel from the reference position from a reference position in which the origin of the template image obtained by cutting out the region of the search object from the standard template image is coincident with the origin of the standard template image; Means for calculating autocorrelation values of the template image and the standard template image respectively at a position moved by +1 pixel from the pixel movement position and a position moved by -1 pixel from the sub-pixel movement position;
Means for calculating a function that approximates a curve from the calculated three autocorrelation values;
Means for calculating a position where the calculated autocorrelation value of the approximate function is maximum;
Means for calculating a difference between the calculated autocorrelation value maximum position and the subpixel movement position as a correction amount;
Means for calculating the correction amount at a position where the sub-pixel movement position is moved to a plurality of different positions;
Means for measuring the position / shape of the search object by correcting the position of the search object with the correction amount of the sub-pixel movement position with the smallest correction amount among the calculated plurality of correction amounts ;
An image measuring device by template matching characterized by having
JP2007263306A 2007-10-09 2007-10-09 Image measuring method and image measuring apparatus Expired - Fee Related JP5034852B2 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2007263306A JP5034852B2 (en) 2007-10-09 2007-10-09 Image measuring method and image measuring apparatus

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2007263306A JP5034852B2 (en) 2007-10-09 2007-10-09 Image measuring method and image measuring apparatus

Publications (2)

Publication Number Publication Date
JP2009093412A JP2009093412A (en) 2009-04-30
JP5034852B2 true JP5034852B2 (en) 2012-09-26

Family

ID=40665343

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2007263306A Expired - Fee Related JP5034852B2 (en) 2007-10-09 2007-10-09 Image measuring method and image measuring apparatus

Country Status (1)

Country Link
JP (1) JP5034852B2 (en)

Families Citing this family (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP5353597B2 (en) * 2009-09-24 2013-11-27 日本電気株式会社 Pattern recognition device
JP5278347B2 (en) * 2010-02-18 2013-09-04 株式会社島津製作所 Material testing machine
EP2734912A1 (en) 2011-07-18 2014-05-28 MultiTouch Oy Correction of touch screen camera geometry
CN104732518B (en) * 2015-01-19 2017-09-01 北京工业大学 A kind of PTAM improved methods based on intelligent robot terrain surface specifications
JP6878908B2 (en) 2017-01-23 2021-06-02 セイコーエプソン株式会社 Encoders, robots and printers
JP7129667B2 (en) * 2017-10-31 2022-09-02 サン電子株式会社 NAIL DETERMINATION METHOD, NAIL DETERMINATION DEVICE, AND PROGRAM
WO2022009318A1 (en) 2020-07-07 2022-01-13 三菱電機株式会社 Visual feedback control device
WO2022118482A1 (en) * 2020-12-04 2022-06-09 日本電信電話株式会社 Image matching device, image matching method, and image matching program
CN116559170B (en) * 2022-01-27 2026-04-24 华为技术有限公司 A product quality testing method and related system

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH10153407A (en) * 1996-11-25 1998-06-09 Omron Corp Measurement processing device

Also Published As

Publication number Publication date
JP2009093412A (en) 2009-04-30

Similar Documents

Publication Publication Date Title
JP5034852B2 (en) Image measuring method and image measuring apparatus
CN107270810B (en) The projector calibrating method and device of multi-faceted projection
US8265343B2 (en) Apparatus, method and program for distance measurement
CN107633536A (en) A kind of camera calibration method and system based on two-dimensional planar template
US20090244299A1 (en) Image processing device, computer-readable storage medium, and electronic apparatus
US20120308126A1 (en) Corresponding image processing method for compensating colour
CN107507246A (en) A kind of camera marking method based on improvement distortion model
WO2015198930A1 (en) Distance measurement device, and distance measurement correction device using correction parameter
CN111080716A (en) Camera calibration target based on color coding phase shift stripes and calibration point extraction method
CN109373901B (en) Method for calculating center position of hole on plane
WO2017130304A1 (en) Measurement device, method and display device
CN107360377B (en) Vehicle-mounted video image stabilization method
KR20160105322A (en) Focus position detection device, focus position detection method and a computer program for focus position detection
Zhang et al. A flexible calibration method for large-range binocular vision system based on state transformation
JP2009301181A (en) Image processing apparatus, image processing program, image processing method and electronic device
CN108510463A (en) A kind of antidote and device of fault image
JP4887376B2 (en) A method for obtaining a dense parallax field in stereo vision
CN114219866B (en) Binocular structured light three-dimensional reconstruction method, reconstruction system and reconstruction equipment
JP2009239515A (en) Image processing apparatus, image processing method, and image processing program
JP6532328B2 (en) Image processing apparatus, control method therefor, and control program
CN116818129B (en) Temperature estimation and thermal distortion correction method applied to structured light reconstruction
JPH1153549A (en) Image processing apparatus, image processing method, and transmission medium
JP5875248B2 (en) Image processing apparatus, image processing method, and program
JP4102386B2 (en) 3D information restoration device
KR101142279B1 (en) An apparatus for aligning images in stereo vision system and the method thereof

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20100715

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20111027

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20111115

A521 Written amendment

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20120113

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

Free format text: JAPANESE INTERMEDIATE CODE: A01

Effective date: 20120605

A01 Written decision to grant a patent or to grant a registration (utility model)

Free format text: JAPANESE INTERMEDIATE CODE: A01

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20120618

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

Free format text: PAYMENT UNTIL: 20150713

Year of fee payment: 3

R150 Certificate of patent or registration of utility model

Free format text: JAPANESE INTERMEDIATE CODE: R150

LAPS Cancellation because of no payment of annual fees