JP4817901B2 - Correlation evaluation system and correlation evaluation method - Google Patents

Correlation evaluation system and correlation evaluation method Download PDF

Info

Publication number
JP4817901B2
JP4817901B2 JP2006073249A JP2006073249A JP4817901B2 JP 4817901 B2 JP4817901 B2 JP 4817901B2 JP 2006073249 A JP2006073249 A JP 2006073249A JP 2006073249 A JP2006073249 A JP 2006073249A JP 4817901 B2 JP4817901 B2 JP 4817901B2
Authority
JP
Japan
Prior art keywords
correlation
value
destination
region
reliability
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.)
Active
Application number
JP2006073249A
Other languages
Japanese (ja)
Other versions
JP2007249682A (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.)
Subaru Corp
Original Assignee
Fuji Jukogyo KK
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 Fuji Jukogyo KK filed Critical Fuji Jukogyo KK
Priority to JP2006073249A priority Critical patent/JP4817901B2/en
Publication of JP2007249682A publication Critical patent/JP2007249682A/en
Application granted granted Critical
Publication of JP4817901B2 publication Critical patent/JP4817901B2/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Measurement Of Optical Distance (AREA)
  • Image Processing (AREA)
  • Image Analysis (AREA)

Description

本発明は、相関元領域の相関先として特定された相関先領域の位置的な信頼性を評価する相関評価システムおよび相関評価方法に関する。   The present invention relates to a correlation evaluation system and a correlation evaluation method for evaluating positional reliability of a correlation destination area specified as a correlation destination of a correlation source area.

従来より、相関元領域の相関先としての相関先領域の位置を詳細に特定する技術が知られている。例えば、特許文献1には、ステレオ画像における相関先領域の位置の特定をサブピクセルレベルで行う手法が開示されている。本手法では、相関先領域のシティブロック距離(以下CB距離という)と、相関先領域の直前または直後に位置的に隣接する相関先候補のCB距離とに基づいて、サブピクセル成分を算出する。このサブピクセル成分を用いて、相関元領域に対する相関先領域のピクセルレベル(データ領域の構成要素レベル)のズレ量を補間することによって、サブピクセルレベルでのズレ量(視差)が算出される。   2. Description of the Related Art Conventionally, a technique for specifying in detail the position of a correlation destination area as a correlation destination of a correlation source area is known. For example, Patent Document 1 discloses a technique for specifying the position of a correlation destination region in a stereo image at a subpixel level. In this method, the subpixel component is calculated based on the city block distance (hereinafter referred to as CB distance) of the correlation destination region and the CB distance of the correlation destination candidate that is positioned immediately before or after the correlation destination region. By using this sub-pixel component, the amount of deviation (parallax) at the sub-pixel level is calculated by interpolating the amount of deviation of the pixel level (component level of the data area) of the correlation destination area with respect to the correlation source area.

特開2003−97944号公報JP 2003-97944 A

しかしながら、従来の技術では、画像データ(データ領域)に乗る、外乱としてのノイズ等によって、算出されたサブピクセルレベルでの相関先領域の位置には、誤差が含まれている。この誤差は、相関先領域の位置を入力とした後段の各種処理の精度を低下させるおそれがある。   However, in the conventional technique, an error is included in the position of the correlation destination region at the calculated subpixel level due to noise as disturbance that is carried on the image data (data region). This error may reduce the accuracy of various subsequent processes using the position of the correlation destination region as an input.

本発明の目的は、相関元領域の相関先として特定された相関先領域の位置的な信頼性を評価することによって、相関先領域の位置を入力とした各種処理の精度向上を図ることである。   An object of the present invention is to improve the accuracy of various processes using the position of a correlation destination area as an input by evaluating the positional reliability of the correlation destination area specified as the correlation destination of the correlation source area. .

また、本発明の別の目的は、ノイズの影響を考慮した上で、相関先領域の位置的な信頼性を評価することによって、耐ノイズ性に優れた各種処理を可能にすることである。   Another object of the present invention is to enable various types of processing with excellent noise resistance by evaluating the positional reliability of the correlation destination region in consideration of the influence of noise.

かかる課題を解決するために、第1の発明は、相関元領域の相関先として特定された相関先領域の位置的な信頼性を評価する相関評価システムを提供する。この評価システムは、マッチング処理部と、補間処理部と、信頼性評価部と、データ補正部とを有する。マッチング処理部は、相関元領域の相関先の候補となる複数の相関先候補を所定の探索範囲内に設定し、複数の相関先候補のそれぞれに関して、相関元領域との相関の程度を示す相関値を算出する。そして、マッチング処理部は、相関先候補毎に算出された相関値同士を比較することによって、複数の相関先候補のいずれかを相関先領域として特定する。また、マッチング処理部は、相関先領域の位置を離散的な第1の座標値として算出する。補間処理部は、相関先領域の相関値と、相関先領域と位置的に隣接する2以上の相関先候補の相関値とを用いた補間によって、相関先領域の位置を前記第1の座標値よりも細かな分解能を有する第2の座標値として算出する。信頼性評価部は、相関先領域の相関値と、相関先領域と位置的に隣接した相関先候補の相関値との間の変化量に応じて、相関値の誤差の範囲を示す第1の誤差幅を相関先領域の位置を示す座標軸上に投影するとともに、当該投影によって座標軸上に生じた第2の誤差幅に基づいて、上記第1の座標値と上記第2の座標値とによって特定される相関先領域の位置に関する信頼性を評価する。データ補正部は、マッチング処理部の前段に設けられており、補正前のデータに対して所定の補正を施した上で、補正後のデータを出力する。この場合、信頼性評価部は、補正前のデータに対する補正後のデータの補正の度合いに応じて、第1の誤差幅を可変に設定する。 In order to solve such a problem, the first invention provides a correlation evaluation system that evaluates the positional reliability of a correlation destination area specified as a correlation destination of a correlation source area. The evaluation system includes a matching processing unit, an interpolation processing unit, a reliability evaluation unit, and a data correction unit . The matching processing unit sets a plurality of correlation destination candidates as correlation destination candidates in the correlation source area within a predetermined search range, and indicates the degree of correlation with the correlation source area for each of the plurality of correlation destination candidates. Calculate the value. Then, the matching processing unit identifies one of the plurality of correlation destination candidates as a correlation destination region by comparing the correlation values calculated for each correlation destination candidate. Further, the matching processing unit calculates the position of the correlation destination region as a discrete first coordinate value. The interpolation processing unit converts the position of the correlation destination region to the first coordinate value by interpolation using a correlation value of the correlation destination region and correlation values of two or more correlation destination candidates that are positionally adjacent to the correlation destination region. The second coordinate value having a finer resolution is calculated. The reliability evaluation unit includes a first value indicating an error range of the correlation value according to an amount of change between the correlation value of the correlation destination region and the correlation value of the correlation destination candidate that is positioned adjacent to the correlation destination region. The error width is projected on the coordinate axis indicating the position of the correlation destination area, and specified by the first coordinate value and the second coordinate value based on the second error width generated on the coordinate axis by the projection. The reliability regarding the position of the correlation destination area is evaluated. The data correction unit is provided in front of the matching processing unit, and performs a predetermined correction on the data before correction and outputs the corrected data. In this case, the reliability evaluation unit variably sets the first error width according to the degree of correction of the data after correction with respect to the data before correction.

第1の発明において、信頼性評価部は、相関先領域の相関値と、相関先領域と位置的に隣接した相関先候補の相関値との間の変化量が小さくなるほど、信頼性を低く評価することが望ましい。   In the first invention, the reliability evaluation unit evaluates the reliability lower as the amount of change between the correlation value of the correlation destination region and the correlation value of the correlation destination candidate positioned adjacent to the correlation destination region is smaller. It is desirable to do.

また、第1の発明において、相関元領域または探索範囲内に含まれるデータに基づいて、当該データを代表する代表値を算出する代表値算出部を設けてもよい。この場合、信頼性評価部は、代表値に応じて、第1の誤差幅を可変に設定する。   In the first invention, a representative value calculation unit that calculates a representative value representing the data based on data included in the correlation source region or the search range may be provided. In this case, the reliability evaluation unit variably sets the first error width according to the representative value.

第1の発明において、センサと、アンプとを設けてもよい。センサは、所定サイズのデータ群によって構成され、かつ、相関元領域または探索範囲として用いられる領域を含むフレーム領域を出力する。アンプは、マッチング処理部の前段に設けられており、センサから出力されたフレーム領域をゲインによって調整する。この場合、信頼性評価部は、ゲインに応じて、第1の誤差幅を可変に設定する。   In the first invention, a sensor and an amplifier may be provided. The sensor outputs a frame region that includes a data group having a predetermined size and includes a region used as a correlation source region or a search range. The amplifier is provided before the matching processing unit, and adjusts the frame region output from the sensor by the gain. In this case, the reliability evaluation unit variably sets the first error width according to the gain.

さらに、第1の発明において、周囲の温度を検出する温度センサをさらに設けてもよい。この場合、信頼性評価部は、温度に基づいて、第1の誤差幅を可変に設定する。   Furthermore, in the first invention, a temperature sensor for detecting the ambient temperature may be further provided. In this case, the reliability evaluation unit variably sets the first error width based on the temperature.

第2の発明は、相関元領域の相関先として特定された相関先領域の位置的な信頼性を評価する相関評価方法を提供する。この評価方法は、補正前のデータに対して所定の補正を施した上で、補正後のデータを出力するステップと、相関元領域の相関先の候補となる複数の相関先候補を所定の探索範囲内に設定し、複数の相関先候補のそれぞれに関して、相関元領域との相関の程度を示す相関値を算出するステップと、相関先候補毎に算出された相関値同士を比較することによって、複数の相関先候補のいずれかを相関先領域として特定するステップと、相関先領域の位置を離散的な第1の座標値として算出するステップと、相関先領域の相関値と、相関先領域と位置的に隣接する2以上の相関先候補の相関値とを用いた補間によって、相関先領域の位置を前記第1の座標値よりも細かな分解能を有する第2の座標値として算出するステップと、相関先領域の相関値と、相関先領域と位置的に隣接した相関先候補の相関値との間の変化量に応じて、相関値の誤差の範囲を示す第1の誤差幅を相関先領域の位置を示す座標軸上に投影するとともに、当該投影によって座標軸上に生じた第2の誤差幅に基づいて、上記第1の座標値と上記第2の座標値とによって特定される相関先領域の位置に関する信頼性を評価するステップとを有する。この信頼性を評価するステップでは、補正前のデータに対する補正後のデータの補正の度合いに応じて、第1の誤差幅を可変に設定する。

The second invention provides a correlation evaluation method for evaluating the positional reliability of the correlation destination area specified as the correlation destination of the correlation source area. In this evaluation method, a predetermined correction is performed on the uncorrected data, and then the corrected data is output, and a plurality of correlation destination candidates that are candidates for the correlation destination in the correlation source region are determined by a predetermined search. By setting the value within the range and calculating the correlation value indicating the degree of correlation with the correlation source region for each of the plurality of correlation destination candidates, by comparing the correlation values calculated for each correlation destination candidate, A step of specifying any one of a plurality of correlation destination candidates as a correlation destination region, a step of calculating a position of the correlation destination region as a discrete first coordinate value, a correlation value of the correlation destination region, a correlation destination region, Calculating a position of the correlation destination region as a second coordinate value having a finer resolution than the first coordinate value by interpolation using correlation values of two or more correlation destination candidates that are adjacent in position; Correlation of correlation destination area And the first error width indicating the error range of the correlation value on the coordinate axis indicating the position of the correlation destination area according to the amount of change between the correlation value and the correlation value of the correlation destination candidate that is positioned adjacent to the correlation destination area. And the reliability of the position of the correlation destination area specified by the first coordinate value and the second coordinate value is evaluated based on the second error width generated on the coordinate axis by the projection. A step of performing. In the step of evaluating the reliability, the first error width is variably set according to the degree of correction of the data after correction with respect to the data before correction.

本発明によれば、相関先領域の相関値を基準とした相関先領域と位置的に隣接した相関先候補の相関値の変化量に応じて、相関先領域の位置に関する信頼性が評価されることになる。この信頼性に関する評価の結果を、相関先領域の位置を入力とした各種処理に用いることによって、これら各種処理の精度向上を図ることができる。例えば、ノイズの影響を考慮して、上記相関値の変化量と信頼性との間の関係を設定することによって、耐ノイズ性に優れた各種処理が可能となる。   According to the present invention, the reliability regarding the position of the correlation destination region is evaluated according to the amount of change in the correlation value of the correlation destination candidate that is positioned adjacent to the correlation destination region based on the correlation value of the correlation destination region. It will be. By using the result of the evaluation related to reliability for various processes using the position of the correlation destination area as an input, the accuracy of these various processes can be improved. For example, by considering the influence of noise and setting the relationship between the change amount of the correlation value and the reliability, various processes with excellent noise resistance can be performed.

(第1の実施形態)
図1は、本実施形態に係る相関評価システムのブロック構成図である。本相関評価システムは、相関元領域の相関先として特定された相関先領域の位置的な信頼性を評価する。記憶部1aには、第1のデータ領域を規定するデータが格納されており、記憶部1bには、第2のデータ領域を規定するデータが格納されている。相関元領域は、第1のデータ領域の一部を構成し、相関先領域(および相関先候補)は、第2のデータ領域の一部を構成する。所定の入力部は、これら一対のデータ領域に相当する一対のデータ群を記憶部1a,1bに出力する。例えば、この一対のデータ群として、同時に撮像された一対のフレーム画像データ(ステレオ画像データ)が挙げられる。この場合、入力部には、ステレオカメラが用いられる。
(First embodiment)
FIG. 1 is a block configuration diagram of a correlation evaluation system according to the present embodiment. This correlation evaluation system evaluates the positional reliability of the correlation destination area specified as the correlation destination of the correlation source area. Data that defines the first data area is stored in the storage unit 1a, and data that defines the second data area is stored in the storage unit 1b. The correlation source area constitutes a part of the first data area, and the correlation destination area (and the correlation destination candidate) constitutes a part of the second data area. The predetermined input unit outputs a pair of data groups corresponding to the pair of data areas to the storage units 1a and 1b. For example, the pair of data groups includes a pair of frame image data (stereo image data) captured simultaneously. In this case, a stereo camera is used for the input unit.

マッチング処理部2は、相関元領域の相関先の候補となる複数の相関先候補を、第2のデータ領域上の所定の探索範囲内に設定する。これらの複数の相関先候補のそれぞれに関して、相関値Sが算出される。この相関値Sは、相関元領域と相関先領域との相関の程度を示す。マッチング処理部2は、相関先候補毎の相関値S同士を比較し、これによって、複数の相関先候補のいずれかを相関先領域として特定する。例えば、二次元平面である第1のデータ領域および第2のデータ領域のi−j座標系は、データ領域の左下隅を原点として、水平方向をi座標軸、垂直方向をj座標軸とする。この座標上で整数値をとる、格子点状のデータ領域の要素に、データが対応付けられる。また、一例として、相関元領域を、4×4要素(水平方向に4要素、垂直方向に4要素)の小領域とすることができる。この相関元領域とデータ特性の相関を有する相関先領域を、例えば横長状の128×4の探索範囲において特定する。探索範囲内では、相関先候補を、探索範囲の左端から右端へと水平方向に1要素ずつずらして更新しながら、相関値Sが算出される。相関先候補は、相関元領域との水平方向のズレ量kによって特定され、この相関値Sは、kの関数S(k)として表される。相関値S(k)を最小とする相関先候補(k=kp)は、相関先領域として特定される。相関先領域の水平方向への整数レベルの座標値kpが、算出され、相関値S(kp)とともに補間処理部3に出力される。また、相関先領域の左側(直前)に隣接する相関先候補の相関値S(kp−1)、および、相関先領域の右側(直後)に隣接する相関先候補の相関値S(kp+1)が、補間処理部3に出力される。ここでは、二次元データ群を用いることを前提に、探索範囲として、第1のデータ領域上の相関元領域を水平方向に移動させた領域に対応する、第2のデータ領域上の範囲を設定している。これに代えて、3次元以上のデータ群を用い、また、探索範囲として、相関元領域を任意の方向に移動させた領域に対応する範囲を設定してもよい。   The matching processing unit 2 sets a plurality of correlation destination candidates as correlation destination candidates in the correlation source area within a predetermined search range on the second data area. A correlation value S is calculated for each of the plurality of correlation destination candidates. The correlation value S indicates the degree of correlation between the correlation source area and the correlation destination area. The matching processing unit 2 compares the correlation values S for each correlation destination candidate, and thereby specifies any one of a plurality of correlation destination candidates as a correlation destination region. For example, in the ij coordinate system of the first data area and the second data area which are two-dimensional planes, the lower left corner of the data area is the origin, the horizontal direction is the i coordinate axis, and the vertical direction is the j coordinate axis. Data is associated with an element of a grid-point data area that takes an integer value on this coordinate. As an example, the correlation source region can be a small region of 4 × 4 elements (4 elements in the horizontal direction and 4 elements in the vertical direction). A correlation destination area having a correlation between the correlation source area and data characteristics is specified in, for example, a horizontally long 128 × 4 search range. In the search range, the correlation value S is calculated while updating the correlation destination candidates by shifting one element at a time in the horizontal direction from the left end to the right end of the search range. The correlation destination candidate is specified by a horizontal shift amount k with respect to the correlation source region, and the correlation value S is expressed as a function S (k) of k. The correlation destination candidate (k = kp) that minimizes the correlation value S (k) is specified as the correlation destination region. The coordinate value kp of the integer level in the horizontal direction of the correlation destination region is calculated and output to the interpolation processing unit 3 together with the correlation value S (kp). Further, the correlation value S (kp−1) of the correlation destination candidate adjacent to the left side (immediately before) of the correlation destination region, and the correlation value S (kp + 1) of the correlation destination candidate adjacent to the right side (immediately after) of the correlation destination region. Is output to the interpolation processing unit 3. Here, on the assumption that a two-dimensional data group is used, a range on the second data area corresponding to an area obtained by moving the correlation source area on the first data area in the horizontal direction is set as a search range. is doing. Instead of this, a three-dimensional or more data group may be used, and a range corresponding to an area obtained by moving the correlation source area in an arbitrary direction may be set as a search range.

相関値S(k)としては、一例として、絶対差分和を用いることができる。この絶対差分和を用いた相関値S(k)は、数式1のように、評価対象となる一対のデータ領域において、位置的に対応する要素同士のデータの差(絶対値)を算出し、この値を全要素分加算した総和値として表される。同数式において、i,jは、相関元領域の位置を特定する相関元領域の左下座標であり、kは、相関元領域から相関先候補への水平方向のズレ量である。Rは相関元領域内の各要素に対応するデータであり、Lは相関先候補内の各要素に対応するデータである。m,nは、相関元領域および相関先候補のサイズである。また、相関値S(k)として、差分自乗和等の他の周知の評価関数を使用することも可能である。差分自乗和を用いた相関値S(k)は、評価対象となる一対のデータ領域において、位置的に対応する要素同士の輝度差の自乗を算出し、この値を全要素分加算した総和値として表される。

Figure 0004817901
As the correlation value S (k), for example, an absolute difference sum can be used. The correlation value S (k) using the sum of absolute differences is calculated as a difference (absolute value) between data of elements corresponding to each other in a pair of data areas to be evaluated, as in Equation 1. This value is expressed as a total value obtained by adding all elements. In the equation, i and j are the lower left coordinates of the correlation source region that specifies the position of the correlation source region, and k is a horizontal shift amount from the correlation source region to the correlation destination candidate. R is data corresponding to each element in the correlation source region, and L is data corresponding to each element in the correlation destination candidate. m and n are the sizes of the correlation source region and the correlation destination candidate. Further, as the correlation value S (k), other known evaluation functions such as a sum of squared differences can be used. The correlation value S (k) using the sum of squared differences is a sum value obtained by calculating the square of the luminance difference between the elements corresponding to each other in the pair of data areas to be evaluated and adding this value for all the elements. Represented as:
Figure 0004817901

これらの相関評価手法を用いた場合、一対のデータ領域の相関が高くなるほど、換言すれば、両者のデータ特性が類似するほど、相関値S(k)が小さくなり、完全に一致するデータ領域に関しては、理論的には相関値S(k)が0になる。ただし、実際には、各種のノイズが外乱として作用するため、本来完全に一致するはずのデータ領域であっても、その相関値S(k)が0になることは殆どない。   When these correlation evaluation methods are used, the correlation value S (k) decreases as the correlation between a pair of data regions increases, in other words, the similarity between the data characteristics of the two data regions, and the data regions that are completely coincident with each other. Theoretically, the correlation value S (k) becomes zero. However, in actuality, various noises act as disturbances, and therefore, even in a data region that should originally be perfectly matched, the correlation value S (k) is hardly zero.

補間処理部3は、相関先領域の相関値S(kp)と、隣接相関先候補の相関値S(kp−1)およびS(kp+1)とを用いて補間する。この補間によって、相関先領域が属するi−j座標系における相関先領域の位置(ズレ量)が小数レベルの座標値として算出される。ここでは、この相関先領域の水平方向のズレ量を、つぎの数式2に示す値kr、つまり、kp(整数部分)とksp(小数部分)との和として特定する。ズレ量krと、相関値S(kp−1),S(kp),S(kp+1)とは信頼性評価部4に出力される。

Figure 0004817901
The interpolation processing unit 3 interpolates using the correlation value S (kp) of the correlation destination region and the correlation values S (kp−1) and S (kp + 1) of the adjacent correlation destination candidates. By this interpolation, the position (displacement amount) of the correlation destination area in the ij coordinate system to which the correlation destination area belongs is calculated as a decimal level coordinate value. Here, the horizontal shift amount of this correlation destination area is specified as the value kr shown in the following Equation 2, that is, the sum of kp (integer part) and ksp (decimal part). The deviation amount kr and the correlation values S (kp−1), S (kp), and S (kp + 1) are output to the reliability evaluation unit 4.
Figure 0004817901

具体的には、補間処理部3は、つぎの数式3,4の右辺を演算して、小数部分kspを得る。数式3は、相関値S(kp−1)が相関値S(kp+1)より大きい場合(換言すれば、相関値S(k)が座標値kpの後よりも前で大きく変化する場合)を示す。また、数式4は、相関値S(kp−1)が相関値S(kp+1)以下の場合(相関値S(k)が座標値kpの前よりも後で大きく変化する場合)を示す。正確には、数式4においては、kspは負となり、この場合、(kp−1)が相関値krの整数部分を示し、(1−ksp)がズレ量krの小数部分を示す。

Figure 0004817901
Figure 0004817901
Specifically, the interpolation processing unit 3 calculates the right side of the following equations 3 and 4 to obtain the decimal part ksp. Equation 3 shows the case where the correlation value S (kp−1) is larger than the correlation value S (kp + 1) (in other words, the case where the correlation value S (k) changes largely before the coordinate value kp). . Equation 4 shows a case where the correlation value S (kp−1) is equal to or smaller than the correlation value S (kp + 1) (a case where the correlation value S (k) changes greatly after the coordinate value kp). To be exact, in Equation 4, ksp is negative. In this case, (kp-1) indicates the integer part of the correlation value kr, and (1-ksp) indicates the decimal part of the shift amount kr.
Figure 0004817901
Figure 0004817901

例えば、ステレオ画像データにおける視差dについてのサブピクセル成分(小数成分)の補間については、特開2000−283753号公報に開示されているので、必要ならば参照されたい。また、これに代えて、特開2003−97944号公報に開示されているサブピクセル成分の補間手法を用いることも可能である。この補間手法においては、上述の3つの相関値のすべては用いられず、相関値S(kp)と、隣接相関値S(kp−1),S(kp+1)のうちのより小さな相関値とが用いられて、サブピクセルレベルの補間が行われる。   For example, the interpolation of the subpixel component (decimal component) for the parallax d in the stereo image data is disclosed in Japanese Patent Laid-Open No. 2000-283653, and should be referred to if necessary. Alternatively, a subpixel component interpolation method disclosed in Japanese Patent Laid-Open No. 2003-97944 can be used instead. In this interpolation method, not all of the above three correlation values are used, and the correlation value S (kp) and the smaller correlation value of the adjacent correlation values S (kp-1) and S (kp + 1) are obtained. Used to perform sub-pixel level interpolation.

信頼性評価部4は、まず、相関値S(kp)を基準とした隣接相関値S(kp−1)またはS(kp+1)の変化量と、誤差幅E1とに基づいて、誤差幅E2を推定する。誤差幅E1は、相関値S(k)の誤差の範囲を示し、誤差幅E2は、相関先領域のズレ量krの誤差の範囲を示す。誤差幅E1として、相関値S(k)に重畳するノイズのばらつきを推定した一定値を設定してよく、また、誤差幅E1として、その誤差を変動させるノイズ要因に応じた可変値を使用してもよい。誤差幅E1に可変値を用いる場合については、第2の実施形態で説明する。つぎに、信頼性評価部4は、誤差幅E2に基づいて、相関先領域のズレ量krに関する信頼性を評価する。   First, the reliability evaluation unit 4 calculates the error width E2 based on the variation amount of the adjacent correlation value S (kp-1) or S (kp + 1) with the correlation value S (kp) as a reference, and the error width E1. presume. The error width E1 indicates the error range of the correlation value S (k), and the error width E2 indicates the error range of the deviation amount kr of the correlation destination region. As the error width E1, a fixed value obtained by estimating the variation of noise superimposed on the correlation value S (k) may be set, and as the error width E1, a variable value corresponding to the noise factor that fluctuates the error is used. May be. The case where a variable value is used for the error width E1 will be described in the second embodiment. Next, the reliability evaluation unit 4 evaluates the reliability related to the shift amount kr of the correlation destination region based on the error width E2.

図2は、相関値S(k)の誤差幅E1と、相関先候補のズレ量kr(またズレ量の小数部分ksp)の誤差幅E2との関係を示す図である。相関値S(kp−1)が相関値S(kp+1)より大きい場合について説明するが、相関値S(kp−1)が相関値S(kp+1)以下の場合についても同様である。同図(a)においては、基準相関値S(kp)と隣接相関値S(kp−1)との間の変化量ΔS1が比較的大きい。横軸を相関先候補のズレ量k、縦軸を相関値S(k)とする座標上で、誤差幅E2を推定する。すなわち、相関値の変化量ΔS1に応じて、縦軸方向の誤差幅E1を横軸方向の誤差幅E2に投影することによって、誤差幅E2が推定される。点P1(kp,S(kp))は、相関先領域に対応し、点P0(kp−1,S(kp−1))は、相関先領域直前の隣接相関先候補に対応する。また、点P2(kp+1,S(kp+1))は、相関先領域直後の隣接相関先候補に対応する。直線l1は、点P0および点P1を通る傾き−m(mは正の数)の直線であり、直線l2は、点P2を通る傾きmの直線である。ここでは、相関値の変化量ΔS1は、相関先候補のズレ量kの1単位の変化量に対応する。したがって、この変化量ΔS1は、直線l1の傾きの大きさmに等しい。また、特に、真値である相関値S(kp−1),S(kp),S(kp+1)が、各値を中心として誤差幅E1の誤差を有することを仮定する。この仮定の下において、誤差幅E1の投影による誤差幅E2は、同図(a)中、点Aと点Bとを結ぶ線分の長さとして算出される。点Aは、直線l1の下側で誤差幅E1を規定する直線と、直線l2の上側で誤差幅E1規定する直線との交点であり、また、点Bは、直線l1の上側で誤差幅E1規定する直線と、直線l2の下側で誤差幅E1を規定する直線との交点である。   FIG. 2 is a diagram showing the relationship between the error width E1 of the correlation value S (k) and the error width E2 of the deviation amount kr of the correlation destination candidate (and the decimal part ksp of the deviation amount). The case where the correlation value S (kp−1) is larger than the correlation value S (kp + 1) will be described, but the same applies to the case where the correlation value S (kp−1) is equal to or smaller than the correlation value S (kp + 1). In FIG. 9A, the change ΔS1 between the reference correlation value S (kp) and the adjacent correlation value S (kp−1) is relatively large. The error width E2 is estimated on the coordinates where the horizontal axis is the deviation amount k of the correlation destination candidate and the vertical axis is the correlation value S (k). That is, the error width E2 is estimated by projecting the error width E1 in the vertical axis direction onto the error width E2 in the horizontal axis direction in accordance with the change amount ΔS1 of the correlation value. Point P1 (kp, S (kp)) corresponds to the correlation destination region, and point P0 (kp-1, S (kp-1)) corresponds to the adjacent correlation destination candidate immediately before the correlation destination region. Further, the point P2 (kp + 1, S (kp + 1)) corresponds to the adjacent correlation destination candidate immediately after the correlation destination region. The straight line l1 is a straight line having a slope −m (m is a positive number) passing through the points P0 and P1, and the straight line l2 is a straight line having a slope m passing through the point P2. Here, the change amount ΔS1 of the correlation value corresponds to a change amount of one unit of the deviation amount k of the correlation destination candidate. Accordingly, the change amount ΔS1 is equal to the magnitude m of the slope of the straight line l1. In particular, it is assumed that correlation values S (kp−1), S (kp), and S (kp + 1), which are true values, have an error with an error width E1 around each value. Under this assumption, the error width E2 resulting from the projection of the error width E1 is calculated as the length of the line segment connecting the point A and the point B in FIG. Point A is an intersection of a straight line defining the error width E1 below the straight line l1 and a straight line defining the error width E1 above the straight line l2, and point B is an error width E1 above the straight line l1. This is the intersection of the prescribed straight line and the straight line defining the error width E1 below the straight line l2.

同図(b)においては、相関値の変化量ΔS1(直線lの傾きの絶対値)が比較的小さい。当然、相関値S(k)の誤差幅E1は、変化量ΔS1の大小に関わらず、同図(a),(b)を比較すると、変化量ΔS1が小さくなるほど、ズレ量krの誤差幅E2が大きくなることがわかる。   In FIG. 5B, the correlation value change amount ΔS1 (the absolute value of the slope of the straight line l) is relatively small. Naturally, the error width E1 of the correlation value S (k) is equal to the magnitude of the change amount ΔS1, and the error width E2 of the shift amount kr becomes smaller as the change amount ΔS1 becomes smaller, comparing FIGS. It turns out that becomes large.

相関値の変化量ΔS1および誤差幅E1に応じて、誤差幅E2が変動することは、つぎに数式5,6に示す小数部分kspの誤差Nkspからも確認される。

Figure 0004817901
Figure 0004817901
It is also confirmed from the error Nksp of the fractional part ksp shown in equations 5 and 6 that the error width E2 fluctuates in accordance with the change amount ΔS1 of the correlation value and the error width E1.
Figure 0004817901
Figure 0004817901

同数式中、N(0,σ)は、平均が0で、標準偏差がσの所定の分布をする変数を表す。このN(0,σ)は、つぎに数式7,8等に示すように、相関値S(k)に重畳するノイズに基づいて定められる。つまり、このN(0,σ)の標準偏差σは、上述のノイズのばらつきを推定した一定値である誤差幅E1に対応する。σには、予め実験やシミュレーション等を通じて確認した値を用いる。さらに、数式5,6によって算出される誤差Nkspのばらつきが、誤差幅E2に対応する。誤差Nkspはその分母の値に応じて変動し、これによって、誤差幅E2も変動する。   In the formula, N (0, σ) represents a variable having a predetermined distribution with an average of 0 and a standard deviation of σ. This N (0, σ) is determined based on the noise superimposed on the correlation value S (k), as shown in Equations 7 and 8 below. That is, the standard deviation σ of N (0, σ) corresponds to the error width E1, which is a constant value obtained by estimating the noise variation. For σ, a value confirmed in advance through experiments or simulations is used. Further, the variation of the error Nksp calculated by the mathematical expressions 5 and 6 corresponds to the error width E2. The error Nksp varies depending on the value of the denominator, and the error width E2 also varies accordingly.

数式5は、隣接相関値S(kp−1)が、隣接相関値S(kp+1)よりも大きい場合(図2の場合)における小数部分kspの誤差Nkspを示す。この場合、誤差Nkspの分母には、相関値の変化量ΔS1が用いられる。この変化量ΔS1は、上述のように、点P0と点P1とを通る直線lの傾きの絶対値を示す。すなわち、変化量ΔS1が大きく、直線lの傾きが急である場合には、誤差幅E2(誤差kspのとりうる値の範囲)はより小さな値をとる。また、変化量ΔS1が小さく、直線lの傾きがなだらかである場合には、誤差幅E2はより大きな値をとる。数式6は、隣接相関値S(kp−1)が、隣接相関値S(kp+1)以下の場合における誤差Nkspを示している。この場合、誤差Nkspの分母には、基準相関値S(kp)と隣接相関値S(kp+1)との間の変化量ΔS2が用いられる。この変化量ΔS2は、点P1と点P2とを通る直線l’の傾きm’を示す。すなわち、変化量ΔS2が大きく、直線l'の傾きm’が急である場合には、誤差幅E2はより小さな値をとる。また、変化量ΔS2が小さく、直線l’の傾きがなだらかである場合には、誤差幅E2はより大きな値をとる。   Equation 5 shows the error Nksp of the fractional part ksp when the adjacent correlation value S (kp-1) is larger than the adjacent correlation value S (kp + 1) (in the case of FIG. 2). In this case, the correlation value variation ΔS1 is used as the denominator of the error Nksp. As described above, the change amount ΔS1 indicates the absolute value of the slope of the straight line 1 passing through the points P0 and P1. That is, when the amount of change ΔS1 is large and the slope of the straight line l is steep, the error width E2 (range of values that the error ksp can take) takes a smaller value. Further, when the change amount ΔS1 is small and the slope of the straight line l is gentle, the error width E2 takes a larger value. Equation 6 shows the error Nksp when the adjacent correlation value S (kp−1) is equal to or smaller than the adjacent correlation value S (kp + 1). In this case, the amount of change ΔS2 between the reference correlation value S (kp) and the adjacent correlation value S (kp + 1) is used as the denominator of the error Nksp. This change amount ΔS2 indicates the slope m ′ of the straight line 1 ′ passing through the points P1 and P2. That is, when the change amount ΔS2 is large and the slope m ′ of the straight line l ′ is steep, the error width E2 takes a smaller value. Further, when the change amount ΔS2 is small and the slope of the straight line l ′ is gentle, the error width E2 takes a larger value.

この数式5に示す誤差Nkspは、つぎに数式7〜10に示すようにして、数式3のズレ量の小数部分kspから導出される。すなわち、数式3において、相関値S(kp−1),S(kp),S(kp+1)に、それぞれ、ノイズ成分n(kp−1),n(kp),n(kp+1)が含まれているとすると、数式7の小数部分ksp'を得る。さらに、ノイズ成分n(kp−1),n(kp),n(kp+1)が、標準偏差σのガウス分布をすると仮定すると、小数部分ksp’は、数式8に示す通りとなる。ここで、N(0,σ)およびN’(0,σ)は、平均値が0で標準偏差がσの所定の分布をする変数である。

Figure 0004817901
Figure 0004817901
The error Nksp shown in Equation 5 is derived from the decimal part ksp of the deviation amount of Equation 3 as shown in Equations 7-10. That is, in Equation 3, the correlation values S (kp-1), S (kp), and S (kp + 1) include noise components n (kp-1), n (kp), and n (kp + 1), respectively. If so, the decimal part ksp ′ of Equation 7 is obtained. Further, assuming that the noise components n (kp−1), n (kp), and n (kp + 1) have a Gaussian distribution with a standard deviation σ, the decimal part ksp ′ is as shown in Equation 8. Here, N (0, σ) and N ′ (0, σ) are variables having a predetermined distribution with an average value of 0 and a standard deviation of σ.
Figure 0004817901
Figure 0004817901

数式8において数式9を仮定すると、数式10の小数部分ksp’が得られ、この数式10の右辺第2式の第2項から、数式5が導かれる。同様に、数式4から、数式11のksp’を経て、数式6が導かれる。

Figure 0004817901
Figure 0004817901
Figure 0004817901
Assuming Equation 9 in Equation 8, the fractional part ksp ′ of Equation 10 is obtained, and Equation 5 is derived from the second term of the second equation on the right side of Equation 10. Similarly, Expression 6 is derived from Expression 4 through ksp ′ of Expression 11.
Figure 0004817901
Figure 0004817901
Figure 0004817901

図1に示す信頼性評価部4においては、さらに、誤差幅E2に基づいて、相関先領域のズレ量krに関する信頼性が評価される。すなわち、例えば、この信頼性の程度を表す信頼値Rを、つぎの数式12,13によって表す。

Figure 0004817901
Figure 0004817901
In the reliability evaluation unit 4 shown in FIG. 1, the reliability of the correlation amount kr of the correlation destination region is further evaluated based on the error width E2. That is, for example, the reliability value R representing the degree of reliability is expressed by the following equations 12 and 13.
Figure 0004817901
Figure 0004817901

この信頼値Rは、相関先領域のズレ量kr(=kp+ksp)の信頼性が高いほど小さな値をとる。その小数部分kspにおける誤差幅E2が比較的大きい場合、すなわち、誤差Nkspを示す上記数式5,6において、分子が一定値の場合、分母が小さくなるほど、信頼値Rはより大きな値をとり、その結果として信頼性は低く評価される。信頼値Rは、相関先領域のズレ量krに対応付けられて、記憶部5に出力される。記憶部5上のこれらのデータkr,Rに基づいて、各種処理が行われる。なお、ここでは、信頼値Rの分母を、相関値の変化量ΔS1およびΔS2のうちのより大きな値のものとしている。つまり、隣接相関値S(kp−1)とS(kp+1)とのうちのより大きな値と、基準相関値S(kp)との間の変化量を信頼値Rの分母としている。これに代えて、信頼値Rの分母に、変化量ΔS1およびΔS2の平均値を用いることが可能である。また、この信頼値Rの分母に、相関値の変化量ΔS1およびΔS2のうちのより小さな値を用いることが可能である。つまり、隣接相関値S(kp−1)とS(kp+1)とのうちのより値の小さなものと、基準相関値S(kp)との間の変化量を信頼値Rの分母とすることが可能である。さらに、信頼値Rは、相関値の誤差幅E1(ノイズのばらつきの推定値σ)を、変化量ΔS1およびΔS2に関するこれらの値で割ったものに、所定の定数をかけたものとしてもよい。すなわち、信頼値Rは、誤差幅E2の定数倍とし、また、誤差幅E2と同一とすることができる。   The reliability value R takes a smaller value as the reliability of the shift amount kr (= kp + ksp) of the correlation destination region is higher. When the error width E2 in the fractional part ksp is relatively large, that is, in the above Equations 5 and 6 indicating the error Nksp, when the numerator is a constant value, the smaller the denominator, the larger the confidence value R becomes. As a result, the reliability is evaluated low. The confidence value R is output to the storage unit 5 in association with the deviation amount kr of the correlation destination region. Various processes are performed based on the data kr and R on the storage unit 5. Here, the denominator of the confidence value R is set to a larger value of the correlation value changes ΔS1 and ΔS2. That is, the amount of change between the larger value of the adjacent correlation values S (kp−1) and S (kp + 1) and the reference correlation value S (kp) is used as the denominator of the confidence value R. Instead, the average value of the change amounts ΔS1 and ΔS2 can be used as the denominator of the reliability value R. Further, a smaller value of the correlation value changes ΔS1 and ΔS2 can be used as the denominator of the reliability value R. That is, the change amount between the smaller value of the adjacent correlation values S (kp−1) and S (kp + 1) and the reference correlation value S (kp) may be used as the denominator of the confidence value R. Is possible. Further, the reliability value R may be a value obtained by dividing the error width E1 of the correlation value (estimated value σ of noise variation) by these values relating to the changes ΔS1 and ΔS2 and a predetermined constant. That is, the reliability value R can be a constant multiple of the error width E2 and can be the same as the error width E2.

このように、本実施形態では、基準相関値と隣接相関値との間の変化量に応じて、相関先領域の位置に関する信頼性が評価されることになる。この信頼性に関する評価結果を用いて、相関先領域の位置を入力とした各種処理を行うことによって、それらの処理の精度向上を図ることができる。特に、ノイズの影響を考慮して、上記相関値の変化量と信頼性との間の関係を設定することによって、耐ノイズ性に優れた各種処理が可能となる。   Thus, in this embodiment, the reliability regarding the position of a correlation destination area | region is evaluated according to the variation | change_quantity between a reference | standard correlation value and an adjacent correlation value. By using the evaluation results regarding the reliability and performing various processes with the position of the correlation destination region as an input, the accuracy of those processes can be improved. In particular, by considering the influence of noise and setting the relationship between the amount of change in the correlation value and the reliability, various processes with excellent noise resistance can be performed.

(第2の実施形態)
上述した相関評価システムに続き、より具体的な装置について詳述する。以下の実施形態は、ステレオ式監視装置への適用例である。図3は、本ステレオ式監視装置のブロック構成図である。上述の相関評価システムと同一の機能を有する部分には同一の符号を付している。本監視装置は、例えば車両に搭載されており、一対のフレーム画像(フレーム領域)を処理して、自車両前方の道路や監視対象物(先行車や歩行者等)を監視する。この監視装置には、特に、センサ6、アンプ7、データ補正部8、温度センサ9および代表値算出部10が設けられている。信頼性評価部4は、アンプ7、データ補正部8、温度センサ9および代表値算出部10におけるノイズに関するデータに基づいて、誤差幅E1を設定する。つまり、ここでは、誤差幅E1に影響するノイズの要因が、主に、アンプ7におけるゲインG、データ補正部8における補正の度合いγ、周囲の温度T、および、相関評価の対象とする輝度D自体であること前提とする。これらのノイズ要因の状態に応じて、信頼性評価部4は、誤差幅E1を可変に設定する。
(Second Embodiment)
A more specific apparatus will be described in detail following the correlation evaluation system described above. The following embodiment is an example applied to a stereo monitoring device. FIG. 3 is a block diagram of the stereo monitoring device. Parts having the same functions as those of the above-described correlation evaluation system are denoted by the same reference numerals. The monitoring apparatus is mounted on a vehicle, for example, and processes a pair of frame images (frame regions) to monitor a road ahead of the host vehicle and a monitoring object (such as a preceding vehicle or a pedestrian). In particular, the monitoring device includes a sensor 6, an amplifier 7, a data correction unit 8, a temperature sensor 9, and a representative value calculation unit 10. The reliability evaluation unit 4 sets the error width E1 based on data regarding noise in the amplifier 7, the data correction unit 8, the temperature sensor 9, and the representative value calculation unit 10. That is, here, the noise factors affecting the error width E1 are mainly the gain G in the amplifier 7, the correction degree γ in the data correction unit 8, the ambient temperature T, and the luminance D to be subjected to correlation evaluation. Assume that it is itself. Depending on the state of these noise factors, the reliability evaluation unit 4 variably sets the error width E1.

センサ6は、CCDやCMOSセンサといったイメージセンサを内蔵した一対のカメラで構成される。このセンサ6は、一対のフレーム画像(基準画像および比較画像)を出力し、フレーム画像は、それぞれ、所定サイズのデータ群により構成される。アンプ7は、センサ6から出力されたフレーム画像をゲインGによって調整する。このゲインG(現在の設定値)は、信頼性評価部4に出力されている。データ補正部8は、補正前の画像データに対して所定の補正を施した上で、補正後の画像データを出力する。例えば、入力データと出力データとの関係が記述されたテーブルを参照することによって、データ補正が行われる。このデータ補正は、後段での処理に適した適宜の特性に設定すればよく、その特性が非線形的であるか否か、或いは、連続的であるか否かを問わない。この補正前の画像データに対する補正後の画像データの補正の度合いγ(現在の設定値)は、信頼性評価部4に出力されている。マッチング処理部2の前段の各部で処理された後の基準画像データは、所定の輝度階調のデジタルデータとして記憶部1aに格納される。同様に、比較画像データが、記憶部1bに格納される。   The sensor 6 is composed of a pair of cameras incorporating an image sensor such as a CCD or CMOS sensor. The sensor 6 outputs a pair of frame images (a reference image and a comparison image), and each frame image is constituted by a data group having a predetermined size. The amplifier 7 adjusts the frame image output from the sensor 6 with the gain G. This gain G (current set value) is output to the reliability evaluation unit 4. The data correction unit 8 performs predetermined correction on the uncorrected image data and then outputs the corrected image data. For example, data correction is performed by referring to a table in which the relationship between input data and output data is described. This data correction may be set to an appropriate characteristic suitable for processing in the subsequent stage, regardless of whether the characteristic is non-linear or continuous. The degree of correction γ (current set value) of the corrected image data with respect to the uncorrected image data is output to the reliability evaluation unit 4. The reference image data after being processed by each part of the preceding stage of the matching processing unit 2 is stored in the storage unit 1a as digital data of a predetermined luminance gradation. Similarly, the comparison image data is stored in the storage unit 1b.

マッチング処理部2は、基準画像の一部である相関元領域の複数の相関先候補を、比較画像上の探索範囲内に設定する。この相関元領域と探索範囲との輝度特性の相関に関して、第1の実施形態において述べたマッチング処理が行われる。このマッチング処理によって、相関先領域の整数レベルの水平方向のズレ量kpが算出される。ズレ量kpは、ステレオ画像における視差dpに対応し、この整数レベルの視差dpが、相関先領域に対応する相関値S(dp)とともに補間処理部3に出力される。さらに、左側の隣接相関先候補の相関値S(dp−1)と、右側の隣接相関先候補の相関値S(dp+1)とが、補間処理部3に出力される。   The matching processing unit 2 sets a plurality of correlation destination candidates in the correlation source region that is a part of the reference image within the search range on the comparison image. The matching process described in the first embodiment is performed on the correlation between the luminance characteristics of the correlation source area and the search range. By this matching process, the horizontal shift amount kp of the integer level of the correlation destination region is calculated. The shift amount kp corresponds to the parallax dp in the stereo image, and this integer level parallax dp is output to the interpolation processing unit 3 together with the correlation value S (dp) corresponding to the correlation destination region. Further, the correlation value S (dp−1) of the left adjacent correlation destination candidate and the correlation value S (dp + 1) of the right adjacent correlation destination candidate are output to the interpolation processing unit 3.

補間処理部3は、相関先領域の相関値S(dp)と、隣接相関先候補の相関値S(dp−1)およびS(dp+1)とを用いた補間によって、視差について小数レベルの座標値dspを算出する。すなわち、整数レベルの視差dpが、この小数レベル座標値dspによって補間される。この補間された視差d(小数レベルの座標値dr)と、相関先領域の相関値S(dp)と、隣接相関先候補の相関値S(dp−1),S(dp+1)とは、信頼性評価部4に出力される。   The interpolation processing unit 3 performs a coordinate value of a decimal level for the parallax by interpolation using the correlation value S (dp) of the correlation destination region and the correlation values S (dp−1) and S (dp + 1) of the adjacent correlation destination candidates. dsp is calculated. That is, the integer level parallax dp is interpolated by this decimal level coordinate value dsp. The interpolated parallax d (decimal level coordinate value dr), the correlation value S (dp) of the correlation destination region, and the correlation values S (dp−1) and S (dp + 1) of the adjacent correlation destination candidates are reliable. It is output to the sex evaluation unit 4.

温度センサ9は、本ステレオ式監視装置の周囲の温度Tを検出し、その現在の検出値を信頼性評価部4に出力する。例えば、センサ6周辺の温度Tが信頼性評価部4に出力される。   The temperature sensor 9 detects the ambient temperature T of the stereo monitoring device and outputs the current detection value to the reliability evaluation unit 4. For example, the temperature T around the sensor 6 is output to the reliability evaluation unit 4.

代表値算出部10は、相関元領域に含まれる輝度データD1〜Dnに基づいて、当該データを代表する代表値Drepを算出する。例えば、この代表値Drepは、相関元領域内の全部もしくは特定の一部の画素の輝度値の平均値とし、また、これらの輝度値の中間値とすることができる。代表値算出部10は、相関元領域に代えて、探索範囲内に含まれる輝度データに基づいて、代表値を算出することも可能である。   Based on the luminance data D1 to Dn included in the correlation source region, the representative value calculation unit 10 calculates a representative value Drep representing the data. For example, the representative value Drep can be an average value of luminance values of all or a specific part of pixels in the correlation source region, and can be an intermediate value of these luminance values. The representative value calculation unit 10 can calculate a representative value based on luminance data included in the search range instead of the correlation source region.

信頼性評価部4は、基準相関値S(dp)と隣接相関値S(dp−1)またはS(dp+1)との間の変化量と、誤差幅E1とに基づいて、誤差幅E2を推定し、この誤差幅E2に基づいて、視差dの信頼性を評価する。ここでの誤差幅E1は、輝度に乗るノイズ要因の状態に基づいて設定されている。すなわち、信頼性評価部4は、ゲインG、補正の度合いγ、周囲の温度T、および、輝度の代表値Drepを検出して、これらに応じたノイズを推定して誤差幅E1を可変に設定する。例えば、誤差幅E1に基礎値と補正値とを設定して、誤差幅E1を設定することができる。すなわち、ゲインG、補正の度合いγおよび温度Tの組み合わせと、誤差幅E1の基礎値との対応関係を所定の基礎値テーブルに予め記述しておく。また、代表値Drepと、誤差幅E1の補正値との対応関係を所定の補正値テーブルに予め記述しておく。信頼性評価部4は、ゲインG、補正の度合いγおよび周囲の温度Tの現在値に基づいて、基礎値テーブルを参照して基礎値を設定し、また、輝度の代表値Drepに基づいて、補正値テーブルを参照して補正値を設定する。そして、信頼性評価部4は、この基礎値を補正値により補正し、これによって、誤差幅E1を設定する。基礎値テーブルおよび補正値テーブルにおける各対応関係は、予め行った実験やシミュレーション等を通じて適切に設定する。この実験等においては、ノイズのばらつきとノイズ要因の状態値との関係が考慮されつつ、輝度に基づく相関値のノイズ成分が推定されている。当然ながら、基礎値テーブルおよび補正値テーブルの全部または一部に代えて所定の計算式を用い、この計算式によって基礎値を算出することも可能である。また、必ずしも、これら4種のノイズ要因のすべてを考慮する必要はなく、これら以外の要因を含めることも可能である。信頼性評価部4は、視差d、および、信頼性の評価結果を表す信頼値Rを記憶部5に格納する。   The reliability evaluation unit 4 estimates the error width E2 based on the amount of change between the reference correlation value S (dp) and the adjacent correlation value S (dp-1) or S (dp + 1) and the error width E1. Then, the reliability of the parallax d is evaluated based on the error width E2. The error width E1 here is set based on the state of the noise factor riding on the luminance. That is, the reliability evaluation unit 4 detects the gain G, the degree of correction γ, the ambient temperature T, and the representative luminance value Drep, and estimates the noise corresponding to these to set the error width E1 variably. To do. For example, the error width E1 can be set by setting a basic value and a correction value for the error width E1. That is, the correspondence relationship between the combination of the gain G, the correction degree γ, and the temperature T and the basic value of the error width E1 is described in advance in a predetermined basic value table. Further, the correspondence between the representative value Drep and the correction value of the error width E1 is described in advance in a predetermined correction value table. The reliability evaluation unit 4 sets a basic value by referring to the basic value table based on the current value of the gain G, the degree of correction γ, and the ambient temperature T, and based on the representative value Drep of the luminance, The correction value is set with reference to the correction value table. Then, the reliability evaluation unit 4 corrects the basic value with the correction value, thereby setting the error width E1. Each correspondence relationship in the basic value table and the correction value table is appropriately set through experiments or simulations performed in advance. In this experiment or the like, the noise component of the correlation value based on the luminance is estimated while considering the relationship between the noise variation and the state value of the noise factor. Of course, it is also possible to use a predetermined calculation formula instead of all or part of the basic value table and the correction value table and calculate the basic value by this calculation formula. Moreover, it is not always necessary to consider all these four types of noise factors, and it is possible to include other factors. The reliability evaluation unit 4 stores the parallax d and the reliability value R representing the reliability evaluation result in the storage unit 5.

この視差dおよび信頼値Rを用いて、対象物の認識処理をより有効に行うことができる。すなわち、周知の視差dを用いた対象物認識においては、隣接する複数の画素領域の視差dが所定の範囲内にあるか否かに応じて、これらの視差dが同一物体、例えば先行車、歩行者などを示すか否かを判定している。この判定結果に基づいて、視差dがグルーピングされる。従来のグルーピングでは、範囲を画する判定しきい値近傍にある値の視差dが、硬直的に判定される。すなわち、その範囲に属していると判定された視差dは、グルーピングされ、その範囲に属していないと判定された視差dは、グルーピングから外される。これに対して、本実施形態のように、視差dの信頼性を考慮することによって、より柔軟に判定を行うことができる。すなわち、視差dに誤差を考慮したグルーピングが可能となる。特に、ノイズ要因の状態に応じて、相関値Sの誤差幅を可変に設定することによって、より適切な信頼値Rを視差dに対して設定することが可能である。   Using this parallax d and the reliability value R, the object recognition process can be performed more effectively. That is, in the object recognition using the well-known parallax d, the parallax d is the same object, for example, a preceding vehicle, depending on whether or not the parallax d of a plurality of adjacent pixel regions is within a predetermined range. It is determined whether or not a pedestrian or the like is indicated. Based on the determination result, the parallax d is grouped. In the conventional grouping, the parallax d having a value near the determination threshold value that defines the range is determined rigidly. That is, the parallax d determined to belong to the range is grouped, and the parallax d determined to not belong to the range is excluded from the grouping. On the other hand, the determination can be performed more flexibly by considering the reliability of the parallax d as in the present embodiment. In other words, it is possible to perform grouping in consideration of errors in the parallax d. In particular, it is possible to set a more appropriate reliability value R for the parallax d by variably setting the error width of the correlation value S according to the state of the noise factor.

また、視差dおよび信頼値Rを用いて、視差dに重み付けを行うことができる。すなわち、信頼値Rによって、その視差dが有効であるか無効であるかを段階的に判定することが可能である。上述のように、この信頼値Rには、ノイズ要因の状態を適切に反映させることができる。   Also, the parallax d can be weighted using the parallax d and the confidence value R. That is, it is possible to determine stepwise whether the parallax d is valid or invalid based on the reliability value R. As described above, the state of the noise factor can be appropriately reflected in the reliability value R.

このように、本実施形態では、基準相関値と隣接相関値との間の変化量に応じて、視差の信頼性が評価されることになる。この信頼性に関する評価結果を、視差を入力とした各種処理に用いることによって、それらの処理の精度向上を図ることができる。特に、ノイズの影響を考慮して、上記相関値の変化量を設定することによって、耐ノイズ性に優れた各種処理が可能となる。   Thus, in this embodiment, the reliability of parallax is evaluated according to the amount of change between the reference correlation value and the adjacent correlation value. By using this reliability evaluation result for various processes using parallax as input, the accuracy of those processes can be improved. In particular, various processes with excellent noise resistance can be performed by setting the amount of change in the correlation value in consideration of the influence of noise.

なお、本実施形態においては、同時に撮像された一対のフレーム画像(ステレオ画像)への適用例について説明したが、本発明はこれに限定されるものではなく、互いに関連性を有する一対のフレーム領域に対して広く適用可能である。このようなフレーム領域の別の例としては、いわゆるオプティカルフローの処理対象となる、異なる時刻に撮像された一対のフレーム画像が挙げられる。この場合、上述したセンサ6の代わりに単眼カメラが用いられ、フレーム画像が時系列的に出力される。また、このようなフレーム領域は、輝度によって規定される画像(輝度分布)に限定されるものではなく、輝度以外の情報(例えば、遠赤外線カメラより得られる熱分布、レーザレーダやミリ波レーダより得られる反射強度の分布等)によって規定されるものであってもよい。   In the present embodiment, the application example to a pair of frame images (stereo images) captured at the same time has been described. However, the present invention is not limited to this, and a pair of frame regions having relevance to each other. Widely applicable. As another example of such a frame region, there is a pair of frame images captured at different times, which are targets of so-called optical flow processing. In this case, a monocular camera is used instead of the sensor 6 described above, and frame images are output in time series. Further, such a frame region is not limited to an image (luminance distribution) defined by luminance, but information other than luminance (for example, heat distribution obtained from a far-infrared camera, laser radar or millimeter wave radar) It may be defined by the obtained reflection intensity distribution or the like.

また、本発明は、パターンマッチングを行う装置に適用することも可能である。ここでは、相関元領域は所定のテンプレートに対応し、また、相関先領域は、フレーム画像の一部または全部を構成する探索領域中の画素領域に対応する。特に、車線に関する道路標示の画像をテンプレートとして、このテンプレートを移動・拡大させながらフレーム画像上の画素領域との相関を評価して対応画素領域を特定する。特定された対応画素領域の位置に関して、誤差幅を推定し、信頼性を評価することが可能である。   The present invention can also be applied to an apparatus that performs pattern matching. Here, the correlation source region corresponds to a predetermined template, and the correlation destination region corresponds to a pixel region in a search region that constitutes part or all of the frame image. In particular, an image of a road marking relating to a lane is used as a template, and a correlation with a pixel region on a frame image is evaluated while moving and enlarging the template, thereby specifying a corresponding pixel region. With respect to the position of the specified corresponding pixel region, it is possible to estimate the error width and evaluate the reliability.

また、本実施形態では、相関元領域が基準画像中に設定され、また、相関先領域が比較画像中で特定されるものとしたが、これに代えて、一方の画像において、相関元領域と、この相関元領域に隣接する相関先候補との双方を設定し、自己マッチングを行ってもよい。すなわち、例えば、基準画像において、ある領域と、前後に隣接する領域との間の相関値の変化量を算出する。この変化量の大きさ(相関値0付近を中央とする谷形状の深さ)に基づいて、比較画像との間で得られる視差の誤差幅を予測し、さらに、視差の信頼性を設定することが可能である。   In the present embodiment, the correlation source region is set in the reference image, and the correlation destination region is specified in the comparison image. Instead, in one image, the correlation source region and the correlation source region are specified. The self-matching may be performed by setting both the correlation destination candidates adjacent to the correlation source region. That is, for example, in the reference image, the amount of change in the correlation value between a certain area and the areas adjacent to each other is calculated. Based on the magnitude of this change amount (depth of the valley shape centered around the correlation value 0), the error width of the parallax obtained with the comparison image is predicted, and further, the parallax reliability is set. It is possible.

さらに、上述した実施形態では、相関先領域の位置を整数レベルの座標値と、小数レベルの座標値とによって特定しているが、これを一般化すると、前者は離散的な第1の座標値、後者は第1の座標値よりも細かな分解能を有する第2の座標値として表わすことができる。整数レベルの座標値は、離散的な値の代表であるが、これ以外にサブピクセル単位の離散値を用いてもよい。例えば、探索範囲内において0.5画素ずつオフセットさせながら、画素補間にて相関先候補を設定することで、0.5画素毎の離散値として座標値を算出するといった如くである。また、隣接する離散値同士の間隔(オフセット値)が一定である必要は必ずしもない。   Further, in the above-described embodiment, the position of the correlation destination region is specified by the coordinate value of the integer level and the coordinate value of the decimal level. When this is generalized, the former is the discrete first coordinate value. The latter can be expressed as a second coordinate value having a finer resolution than the first coordinate value. The integer level coordinate values are representative of discrete values, but other than these, discrete values in units of subpixels may be used. For example, the coordinate value is calculated as a discrete value for each 0.5 pixel by setting a correlation destination candidate by pixel interpolation while offset by 0.5 pixel within the search range. Further, the interval (offset value) between adjacent discrete values is not necessarily constant.

第1の実施形態の相関評価システムのブロック構成図The block block diagram of the correlation evaluation system of 1st Embodiment 相関値の誤差幅と、相関先領域のズレ量の誤差幅との関係を示す図A diagram showing the relationship between the error width of the correlation value and the error width of the deviation amount of the correlation destination area 第2の実施形態のステレオ式監視装置のブロック構成図Block configuration diagram of the stereo monitoring device of the second embodiment

符号の説明Explanation of symbols

2 マッチング処理部
3 補間処理部
4 信頼性評価部
6 センサ
7 アンプ
8 データ補正部
9 温度センサ
10 代表値算出部
2 Matching Processing Unit 3 Interpolation Processing Unit 4 Reliability Evaluation Unit 6 Sensor 7 Amplifier 8 Data Correction Unit 9 Temperature Sensor 10 Representative Value Calculation Unit

Claims (6)

相関元領域の相関先として特定された相関先領域の位置的な信頼性を評価する相関評価システムにおいて、
前記相関元領域の相関先の候補となる複数の相関先候補を所定の探索範囲内に設定し、前記複数の相関先候補のそれぞれに関して、前記相関元領域との相関の程度を示す相関値を算出し、前記相関先候補毎に算出された相関値同士を比較することによって、前記複数の相関先候補のいずれかを前記相関先領域として特定するとともに、前記相関先領域の位置を離散的な第1の座標値として算出するマッチング処理部と、
前記相関先領域の相関値と、前記相関先領域と位置的に隣接する2以上の前記相関先候補の相関値とを用いた補間によって、前記相関先領域の位置を前記第1の座標値よりも細かな分解能を有する第2の座標値として算出する補間処理部と、
前記相関先領域の相関値と、前記相関先領域と位置的に隣接した前記相関先候補の相関値との間の変化量に応じて、前記相関値の誤差の範囲を示す第1の誤差幅を前記相関先領域の位置を示す座標軸上に投影するとともに、当該投影によって前記座標軸上に生じた第2の誤差幅に基づいて、前記第1の座標値と前記第2の座標値とによって特定される前記相関先領域の位置に関する信頼性を評価する信頼性評価部と
前記マッチング処理部の前段に設けられており、補正前のデータに対して所定の補正を施した上で、補正後のデータを出力するデータ補正部とを有し、
前記信頼性評価部は、前記補正前のデータに対する前記補正後のデータの補正の度合いに応じて、前記第1の誤差幅を可変に設定することを特徴とする相関評価システム。
In the correlation evaluation system that evaluates the positional reliability of the correlation destination area specified as the correlation destination of the correlation source area,
A plurality of correlation destination candidates as correlation destination candidates in the correlation source region are set within a predetermined search range, and a correlation value indicating a degree of correlation with the correlation source region is set for each of the plurality of correlation destination candidates. By calculating and comparing correlation values calculated for each of the correlation destination candidates, one of the plurality of correlation destination candidates is specified as the correlation destination region, and the position of the correlation destination region is discrete A matching processing unit that calculates the first coordinate value;
The position of the correlation destination region is determined from the first coordinate value by interpolation using the correlation value of the correlation destination region and the correlation values of the two or more correlation destination candidates that are positionally adjacent to the correlation destination region. An interpolation processing unit that calculates a second coordinate value having a fine resolution;
The first error width indicating the error range of the correlation value according to the amount of change between the correlation value of the correlation destination region and the correlation value of the correlation destination candidate that is positioned adjacent to the correlation destination region Is projected on the coordinate axis indicating the position of the correlation destination area, and specified by the first coordinate value and the second coordinate value based on the second error width generated on the coordinate axis by the projection. A reliability evaluation unit that evaluates the reliability related to the position of the correlation destination region ,
A data correction unit that is provided in a preceding stage of the matching processing unit, performs a predetermined correction on the data before correction, and outputs data after correction;
The correlation evaluation system, wherein the reliability evaluation unit variably sets the first error width according to a degree of correction of the data after correction with respect to the data before correction .
前記信頼性評価部は、前記相関先領域の相関値と、前記相関先領域と位置的に隣接した前記相関先候補の相関値との間の変化量が小さくなるほど、前記信頼性を低く評価することを特徴とする請求項1に記載された相関評価システム。   The reliability evaluation unit evaluates the reliability lower as the amount of change between the correlation value of the correlation destination region and the correlation value of the correlation destination candidate positioned adjacent to the correlation destination region becomes smaller. The correlation evaluation system according to claim 1, wherein: 前記相関元領域または前記探索範囲内に含まれるデータに基づいて、当該データを代表する代表値を算出する代表値算出部とをさらに有し、
前記信頼性評価部は、前記代表値に応じて、前記第1の誤差幅を可変に設定することを特徴とする請求項1または2に記載された相関評価システム。
A representative value calculating unit that calculates a representative value representing the data based on the data included in the correlation source region or the search range;
The correlation evaluation system according to claim 1 or 2, wherein the reliability evaluation unit variably sets the first error width according to the representative value.
所定サイズのデータ群によって構成され、かつ、前記相関元領域または前記探索範囲として用いられる領域を含むフレーム領域を出力するセンサと、
前記マッチング処理部の前段に設けられており、前記センサから出力されたフレーム領域をゲインによって調整するアンプとをさらに有し、
前記信頼性評価部は、前記ゲインに応じて、前記第1の誤差幅を可変に設定することを特徴とする請求項1から3のいずれかに記載された相関評価システム。
A sensor configured by a data group of a predetermined size and outputting a frame region including a region used as the correlation source region or the search range;
An amplifier that is provided in a preceding stage of the matching processing unit, and that adjusts a frame region output from the sensor by a gain;
The correlation evaluation system according to claim 1, wherein the reliability evaluation unit variably sets the first error width according to the gain.
周囲の温度を検出する温度センサをさらに有し、
前記信頼性評価部は、前記温度に基づいて、前記第1の誤差幅を可変に設定することを特徴とする請求項1から4のいずれかに記載された相関評価システム。
A temperature sensor for detecting the ambient temperature;
The correlation evaluation system according to claim 1, wherein the reliability evaluation unit variably sets the first error width based on the temperature.
相関元領域の相関先として特定された相関先領域の位置的な信頼性を評価する相関評価方法において、
補正前のデータに対して所定の補正を施した上で、補正後のデータを出力するステップと、
前記相関元領域の相関先の候補となる複数の相関先候補を所定の探索範囲内に設定し、前記複数の相関先候補のそれぞれに関して、前記相関元領域との相関の程度を示す相関値を算出するステップと、
前記相関先候補毎に算出された相関値同士を比較することによって、前記複数の相関先候補のいずれかを前記相関先領域として特定するステップと、
前記相関先領域の位置を離散的な第1の座標値として算出するステップと、
前記相関先領域の相関値と、前記相関先領域と位置的に隣接する2以上の前記相関先候補の相関値とを用いた補間によって、前記相関先領域の位置を前記第1の座標値よりも細かな分解能を有する第2の座標値として算出するステップと、
前記相関先領域の相関値と、前記相関先領域と位置的に隣接した前記相関先候補の相関値との間の変化量に応じて、前記相関値の誤差の範囲を示す第1の誤差幅を前記相関先領域の位置を示す座標軸上に投影するとともに、当該投影によって前記座標軸上に生じた第2の誤差幅に基づいて、前記第1の座標値と前記第2の座標値とによって特定される前記相関先領域の位置に関する信頼性を評価するステップとを有し、
前記信頼性を評価するステップでは、前記補正前のデータに対する前記補正後のデータの補正の度合いに応じて、前記第1の誤差幅を可変に設定することを特徴とする相関評価方法。
In the correlation evaluation method for evaluating the positional reliability of the correlation destination area specified as the correlation destination of the correlation source area,
A step of outputting the corrected data after performing a predetermined correction on the data before correction;
A plurality of correlation destination candidates as correlation destination candidates in the correlation source region are set within a predetermined search range, and a correlation value indicating a degree of correlation with the correlation source region is set for each of the plurality of correlation destination candidates. A calculating step;
Identifying any of the plurality of correlation destination candidates as the correlation destination region by comparing correlation values calculated for each of the correlation destination candidates; and
Calculating the position of the correlation destination region as discrete first coordinate values;
The position of the correlation destination region is determined from the first coordinate value by interpolation using the correlation value of the correlation destination region and the correlation values of the two or more correlation destination candidates that are positionally adjacent to the correlation destination region. Calculating a second coordinate value having a fine resolution,
The first error width indicating the error range of the correlation value according to the amount of change between the correlation value of the correlation destination region and the correlation value of the correlation destination candidate that is positioned adjacent to the correlation destination region Is projected on the coordinate axis indicating the position of the correlation destination area, and specified by the first coordinate value and the second coordinate value based on the second error width generated on the coordinate axis by the projection. and a step of the evaluating the reliability of the position of the correlated destination area being,
In the step of evaluating the reliability , the correlation evaluation method is characterized in that the first error width is variably set according to the degree of correction of the corrected data with respect to the uncorrected data .
JP2006073249A 2006-03-16 2006-03-16 Correlation evaluation system and correlation evaluation method Active JP4817901B2 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2006073249A JP4817901B2 (en) 2006-03-16 2006-03-16 Correlation evaluation system and correlation evaluation method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2006073249A JP4817901B2 (en) 2006-03-16 2006-03-16 Correlation evaluation system and correlation evaluation method

Publications (2)

Publication Number Publication Date
JP2007249682A JP2007249682A (en) 2007-09-27
JP4817901B2 true JP4817901B2 (en) 2011-11-16

Family

ID=38593894

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2006073249A Active JP4817901B2 (en) 2006-03-16 2006-03-16 Correlation evaluation system and correlation evaluation method

Country Status (1)

Country Link
JP (1) JP4817901B2 (en)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP5955292B2 (en) * 2013-09-30 2016-07-20 富士重工業株式会社 Filtering device

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP4435525B2 (en) * 2003-09-17 2010-03-17 富士重工業株式会社 Stereo image processing device
JP4394487B2 (en) * 2004-03-05 2010-01-06 富士重工業株式会社 Stereo image processing device

Also Published As

Publication number Publication date
JP2007249682A (en) 2007-09-27

Similar Documents

Publication Publication Date Title
CN110567469B (en) Visual positioning method and device, electronic equipment and system
US9576367B2 (en) Object detection method and device
JP4919036B2 (en) Moving object recognition device
US6985619B1 (en) Distance correcting apparatus of surroundings monitoring system and vanishing point correcting apparatus thereof
JP6167525B2 (en) Distance measuring device and vehicle
US10659762B2 (en) Stereo camera
US20090237491A1 (en) Object Detecting System
JP4554316B2 (en) Stereo image processing device
JP6970577B2 (en) Peripheral monitoring device and peripheral monitoring method
JP2008039491A (en) Stereo image processing apparatus
EP2511875A1 (en) Apparatus and method for refining a value of a similarity measure
JP2009139323A (en) Travel road surface detecting apparatus for vehicle
JP6395429B2 (en) Image processing apparatus, control method thereof, and storage medium
JP2009139324A (en) Travel road surface detecting apparatus for vehicle
US9384576B2 (en) Method and device for computing a change in an image scale of an object
CN114919584A (en) Motor vehicle fixed point target distance measuring method and device and computer readable storage medium
CN111989541A (en) Stereo camera device
JP4817901B2 (en) Correlation evaluation system and correlation evaluation method
CN111260538B (en) Positioning and vehicle-mounted terminal based on long-baseline binocular fisheye camera
JP2009129318A (en) Image processor
JP2014138331A (en) Imaging apparatus and program
JP2001108434A (en) Method and apparatus for measuring distance
JP2023128658A (en) Image processing device and image processing method
WO2020036039A1 (en) Stereo camera device
JP7064400B2 (en) Object detection device

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20081003

RD02 Notification of acceptance of power of attorney

Free format text: JAPANESE INTERMEDIATE CODE: A7422

Effective date: 20101021

RD04 Notification of resignation of power of attorney

Free format text: JAPANESE INTERMEDIATE CODE: A7424

Effective date: 20101119

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20110118

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20110125

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20110311

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20110517

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20110705

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

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

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

Free format text: PAYMENT UNTIL: 20140909

Year of fee payment: 3

R150 Certificate of patent or registration of utility model

Ref document number: 4817901

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150

Free format text: JAPANESE INTERMEDIATE CODE: R150

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

S531 Written request for registration of change of domicile

Free format text: JAPANESE INTERMEDIATE CODE: R313531

R350 Written notification of registration of transfer

Free format text: JAPANESE INTERMEDIATE CODE: R350

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

S533 Written request for registration of change of name

Free format text: JAPANESE INTERMEDIATE CODE: R313533

R350 Written notification of registration of transfer

Free format text: JAPANESE INTERMEDIATE CODE: R350

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250