JP2005250994A - Stereo image processor - Google Patents

Stereo image processor Download PDF

Info

Publication number
JP2005250994A
JP2005250994A JP2004062661A JP2004062661A JP2005250994A JP 2005250994 A JP2005250994 A JP 2005250994A JP 2004062661 A JP2004062661 A JP 2004062661A JP 2004062661 A JP2004062661 A JP 2004062661A JP 2005250994 A JP2005250994 A JP 2005250994A
Authority
JP
Japan
Prior art keywords
image
pixel
stereo
resolution
data
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.)
Granted
Application number
JP2004062661A
Other languages
Japanese (ja)
Other versions
JP4394487B2 (en
Inventor
Tsutomu Tanzawa
勉 丹沢
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 Heavy Industries 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 Fuji Heavy Industries Ltd filed Critical Fuji Heavy Industries Ltd
Priority to JP2004062661A priority Critical patent/JP4394487B2/en
Publication of JP2005250994A publication Critical patent/JP2005250994A/en
Application granted granted Critical
Publication of JP4394487B2 publication Critical patent/JP4394487B2/en
Anticipated expiration legal-status Critical
Expired - Fee Related legal-status Critical Current

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T3/00Geometric image transformation in the plane of the image
    • G06T3/40Scaling the whole image or part thereof
    • G06T3/4007Interpolation-based scaling, e.g. bilinear interpolation

Abstract

<P>PROBLEM TO BE SOLVED: To perform stereo matching with high accuracy to achieve enhancement of data reliability and reduction of variation. <P>SOLUTION: The resolution is expanded by performing interpolation between pixels to a pair of stereo images picked up by a stereo camera 10, to which input processing is performed by an image input part 20 and stored in an original image memory 25 by a resolution expansion part 30, and the stereo matching is performed by an SAD operation to the image whose resolution is expanded by a stereo matching part 35. Then, the stereo matching with high accuracy is performed to achieve enhancement of the data reliability and reduction of the variation by finding parallax at a subpixel level with high accuracy by performing interpolation processing by thee points as a point of agreement defined by the minimum value of an SAD operation value and adjacent points before and after the agreement point by an interpolation processing part 40. <P>COPYRIGHT: (C)2005,JPO&NCIPI

Description

本発明は、互いに相関を有する画像対を用いてステレオマッチングを行うステレオ画像処理装置に関する。   The present invention relates to a stereo image processing apparatus that performs stereo matching using image pairs having a correlation with each other.

一般に、画像による3次元計測技術として、ステレオカメラで対象物を異なる位置から撮像した一対の画像の相関を求め、同一物体に対する視差からステレオカメラの取り付け位置や焦点距離等のカメラパラメータを用いて三角測量の原理により距離を求める、いわゆるステレオ法による画像処理が知られている。   In general, as a three-dimensional measurement technique using an image, a correlation between a pair of images obtained by capturing an object from different positions with a stereo camera is obtained, and a triangle is obtained from parallax with respect to the same object using camera parameters such as the stereo camera mounting position and focal length. Image processing based on a so-called stereo method for obtaining a distance by the principle of surveying is known.

このステレオ法による画像処理としては、本出願人は、先に、特開平5−114099号公報において、車両に搭載したステレオカメラで撮像した画像を処理して車外の対象物の3次元位置を測定する技術を提案しており、この技術では、ステレオ撮像した一対の画像に対し、複数の小領域間の相関度を演算して対応位置を探索するステレオマッチングをハードウエア回路によって高速に行い、対応する小領域の画素ズレ量を視差として出力するようにしている。
特開平5−114099号公報
As image processing by this stereo method, the present applicant previously measured the three-dimensional position of an object outside the vehicle by processing an image picked up by a stereo camera mounted on the vehicle in Japanese Patent Laid-Open No. 5-11409. In this technology, stereo matching that searches the corresponding position by calculating the degree of correlation between multiple small areas is performed at high speed using a hardware circuit for a pair of stereo images. The pixel shift amount of the small area to be output is output as parallax.
Japanese Patent Laid-Open No. 5-114099

ところで、ステレオマッチングにより得られた視差は、画像平面座標から実空間座標への変換により実空間における3次元の距離データとして算出され、車両の車外状況認識システム、高度道路交通システム(ITS)、ロボットビジョン、侵入物監視システム、飛行体の対地高度計測システム等の各種の認識処理に適用される。   By the way, the parallax obtained by the stereo matching is calculated as three-dimensional distance data in the real space by converting the image plane coordinates to the real space coordinates, and the vehicle exterior situation recognition system, intelligent road traffic system (ITS), robot It is applied to various recognition processes such as vision, intruder monitoring system, and aircraft altitude measurement system.

このような認識処理を信頼性高く行うには、視差を算出するステレオマッチングにおけるミスマッチングやデータのばらつきを低減し、マッチング精度を高める必要がある。このミスマッチングを抑制するには、視差の算出単位となるマッチング領域を大きくすること、すなわちマッチングのデータ数を増やすことが有効であるが、従来のステレオマッチングは画像の“画素”を単位としているため、マッチング領域を大きくすると、算出される視差の絶対数が不足し、距離データの空間分解能が低下するという不都合がある。   In order to perform such recognition processing with high reliability, it is necessary to reduce mismatching and data variation in stereo matching for calculating parallax and to increase matching accuracy. In order to suppress this mismatching, it is effective to increase the matching area that is a unit for calculating the parallax, that is, to increase the number of matching data. Conventional stereo matching is based on the “pixel” of the image. Therefore, when the matching area is enlarged, there is a disadvantage that the absolute number of calculated parallaxes is insufficient and the spatial resolution of the distance data is lowered.

本発明は上記事情に鑑みてなされたもので、分解能の低下を招くことなく精度の高いステレオマッチングを可能とし、データ信頼性の向上及びばらつき低減を達成することのできるステレオ画像処理装置を提供することを目的としている。   The present invention has been made in view of the above circumstances, and provides a stereo image processing apparatus that enables high-precision stereo matching without causing a reduction in resolution, and that can achieve improved data reliability and reduced variation. The purpose is that.

上記目的を達成するため、本発明によるステレオ画像処理装置は、互いに相関を有する画像対を用いてステレオマッチングを行うステレオ画像処理装置であって、上記画像対の各画素間に、周辺画素のデータを用いて生成した仮想的な画素を挿入し、上記画像対の解像度を拡張する解像度拡張手段と、上記解像度拡張手段で解像度を拡張した画像対に対し、拡張した解像度の分解能で互いの対応位置を特定するステレオマッチング手段とを備えたことを特徴とする。   In order to achieve the above object, a stereo image processing apparatus according to the present invention is a stereo image processing apparatus that performs stereo matching using image pairs having a correlation with each other, and includes data of peripheral pixels between the pixels of the image pair. The virtual pixel generated by using and inserting the virtual pixel to expand the resolution of the image pair, and the image pair expanded by the resolution expanding unit, the corresponding position of each other at the resolution of the expanded resolution And stereo matching means for specifying the above.

また、本発明によるステレオ画像処理装置は、互いに相関を有する画像対を用いてステレオマッチングを行うステレオ画像処理装置であって、上記画像対の各画素のデータを、周辺画素のデータを用いて置き換える画像変換手段と、上記画像変換手段で画素のデータを置き換えた画像対に対し、互いの対応位置を特定するステレオマッチング手段とを備えたことを特徴とする。   A stereo image processing apparatus according to the present invention is a stereo image processing apparatus that performs stereo matching using image pairs having correlation with each other, and replaces data of each pixel of the image pair with data of peripheral pixels. The image conversion means and stereo matching means for specifying the corresponding positions of the image pairs in which the pixel data is replaced by the image conversion means are provided.

また、本発明によるステレオ画像処理装置は、互いに相関を有する画像対を用いてステレオマッチングを行うステレオ画像処理装置であって、上記画像対のうち、少なくとも一方の画像に対して、各画素のデータを周辺画素のデータを用いて置き換える手段と、置き換えられた画素データを用いて、上記画像対における互いの対応位置を特定するステレオマッチング手段とを備えたことを特徴とする。   A stereo image processing apparatus according to the present invention is a stereo image processing apparatus that performs stereo matching using image pairs having correlation with each other, and each pixel data for at least one of the image pairs. Means for replacing the image data using peripheral pixel data, and stereo matching means for specifying the corresponding positions in the image pair using the replaced pixel data.

本発明のステレオ画像処理装置は、分解能の低下を招くことなく精度の高いステレオマッチングを可能とし、データ信頼性の向上及びばらつき低減を達成することができる。   The stereo image processing apparatus of the present invention enables high-precision stereo matching without causing a reduction in resolution, and can achieve improved data reliability and reduced variation.

以下、図面を参照して本発明の実施の形態を説明する。図1〜図9は本発明の実施の第1形態に係り、図1はステレオ画像処理装置の基本構成を示すブロック図、図2はメイン画像とサブ画像とのマッチングを示す説明図、図3はシティブロック距離の分布を示す説明図、図4はピクセル単位のマッチング後の補間を示す説明図、図5はステレオマッチング及びサブピクセル処理のフローチャート、図6はメイン画像の解像度拡張を示す説明図、図7はサブ画像の解像度拡張を示す説明図、図8はサブピクセル単位のマッチング後の補間を示す説明図、図9は視差のばらつきを示すヒストグラムである。   Embodiments of the present invention will be described below with reference to the drawings. 1 to 9 relate to a first embodiment of the present invention, FIG. 1 is a block diagram showing a basic configuration of a stereo image processing apparatus, FIG. 2 is an explanatory diagram showing matching between a main image and a sub image, and FIG. Is an explanatory diagram showing the distribution of city block distances, FIG. 4 is an explanatory diagram showing interpolation after matching in units of pixels, FIG. 5 is a flowchart of stereo matching and sub-pixel processing, and FIG. 6 is an explanatory diagram showing resolution expansion of the main image FIG. 7 is an explanatory diagram showing the resolution expansion of the sub-image, FIG. 8 is an explanatory diagram showing the interpolation after matching in units of sub-pixels, and FIG. 9 is a histogram showing the parallax variation.

図1において、符号1は、被写体をステレオ撮像した一対の画像を処理して3次元の距離情報を求めるステレオ画像処理装置であり、例えば、車両の車外状況認識システム、高度道路交通システム(ITS)、ロボットビジョン、侵入物監視システム、飛行体の対地高度計測システム等の各種システムに適用される。   In FIG. 1, reference numeral 1 denotes a stereo image processing apparatus that obtains three-dimensional distance information by processing a pair of images obtained by stereo imaging of a subject. For example, a vehicle exterior situation recognition system, an intelligent road traffic system (ITS) It is applied to various systems such as robot vision, intruder monitoring system, and flying object ground altitude measurement system.

ステレオ画像処理装置1は、2台1組のカメラから構成されるステレオカメラ10、このステレオカメラ10で撮像した一対の画像に対する入力処理を行う画像入力部20、この画像入力部で処理した画像を元画像としてストアする元画像メモリ25、一対の元画像の解像度を拡張する解像度拡張部30、一対の画像間の相関度を評価して互いの対応位置を決定するステレオマッチング部35、ステレオマッチング部35で決定した対応位置に基づく補間処理により1画素単位以下の分解能を有するサブピクセルレベルの視差を求める補間処理部40、視差に基づく距離情報を元画像の画像座標と対応付けて二次元的に配置した画像形態(距離画像)としてストアする距離画像メモリ45、元画像メモリ25の元画像情報及び距離画像メモリ45の距離情報を用いて各種認識処理を行う認識処理部50を備えて構成される。   The stereo image processing apparatus 1 includes a stereo camera 10 composed of a pair of cameras, an image input unit 20 that performs input processing on a pair of images captured by the stereo camera 10, and images processed by the image input unit. An original image memory 25 that stores the original image, a resolution expanding unit 30 that expands the resolution of the pair of original images, a stereo matching unit 35 that evaluates the degree of correlation between the pair of images and determines the corresponding positions, and a stereo matching unit Interpolation processing unit 40 that obtains sub-pixel level parallax having a resolution of one pixel unit or less by interpolation processing based on the corresponding position determined in 35, and distance information based on the parallax is associated with the image coordinates of the original image in two dimensions. Distance image memory 45 to store as an arranged image form (distance image), original image information in original image memory 25 and distance image memory With 5 distance information configured with a recognition processing unit 50 for performing various recognition processes.

ステレオカメラ10を構成する2台のカメラ10a,10bは、撮像素子としてCCDやCMOSイメージセンサ等を有し、互いに同期が取れ且つシャッタースピード可変のカメラである。一方のカメラ10aは、ステレオ処理の際の基準画像(メイン画像)を撮像するメインカメラ、他方のカメラ10bは、ステレオ処理の際の比較画像(サブ画像)を撮像するサブカメラとして用いられる。本形態においては、2台のカメラ10a,10bは、互いの光軸が平行で水平面上で所定の基線長をもって配置されているものとする。   The two cameras 10a and 10b constituting the stereo camera 10 have a CCD, a CMOS image sensor, or the like as an image sensor, and are cameras that can be synchronized with each other and have a variable shutter speed. One camera 10a is used as a main camera that captures a reference image (main image) during stereo processing, and the other camera 10b is used as a sub camera that captures a comparative image (sub image) during stereo processing. In the present embodiment, it is assumed that the two cameras 10a and 10b are arranged with a predetermined baseline length on a horizontal plane with their optical axes parallel to each other.

画像入力部20は、2系統のアナログ画像信号を処理するための各種回路(例えば、ゲインコントロールアンプ、ルックアップテーブル、シェーディング補正回路等)や、アナログ信号をデジタル信号に変換するA/Dコンバータ等の画像処理の各種機能回路を備え、更に、カメラ10a,10bの機械的な光学位置の僅かなズレを、画像の回転や平行移動等の幾何学的な変換によって補正する画像調整回路を備えて構成される。   The image input unit 20 includes various circuits (for example, a gain control amplifier, a lookup table, a shading correction circuit, etc.) for processing two systems of analog image signals, an A / D converter that converts analog signals into digital signals, and the like. In addition, an image adjustment circuit that corrects slight misalignment of the mechanical optical positions of the cameras 10a and 10b by geometrical transformation such as image rotation and translation is provided. Composed.

すなわち、画像入力部20に入力されたカメラ10a,10bからの撮像画像は、各カメラ10a,10bの感度差やカメラ光学系の歪みに対する補正がなされた後、所定の輝度階調のデジタル画像データに変換され、更に、画像調整によってカメラ10a,10bの機械的な取り付け位置の誤差が補正され、メイン画像とサブ画像とのステレオ画像対におけるエピポーララインの一致が保証される。このような入力処理を経たステレオ画像対は、元画像メモリ25にストアされる。   That is, the captured images from the cameras 10a and 10b input to the image input unit 20 are corrected for sensitivity differences between the cameras 10a and 10b and distortion of the camera optical system, and then digital image data having a predetermined luminance gradation. Furthermore, the error of the mechanical attachment position of the cameras 10a and 10b is corrected by image adjustment, and the coincidence of epipolar lines in the stereo image pair of the main image and the sub image is guaranteed. The stereo image pair that has undergone such input processing is stored in the original image memory 25.

元画像メモリ25にストアされたステレオ画像対からは、メイン画像とサブ画像との間の相関度を評価しながら対応位置を探索する(ステレオマッチング)ことで視差が得られ、この視差から被写体までの距離情報を求めることができる。メイン画像とサブ画像との相関度は、周知の領域探索法を用い、相関度の評価関数として、メイン画像の小領域とサブ画像の小領域との間のピクセル値の差分(絶対値)の総和(SAD;Sum of Absolute Difference)を演算する。ピクセル値としては、一般的に、各画素の輝度値を用いることが多い。   From the stereo image pair stored in the original image memory 25, parallax is obtained by searching for a corresponding position while evaluating the degree of correlation between the main image and the sub image (stereo matching). Distance information can be obtained. The correlation between the main image and the sub-image is obtained by using a well-known region search method, and the difference (absolute value) of pixel values between the small region of the main image and the small region of the sub-image as an evaluation function of the correlation. Sum (Absolute Difference) is calculated. In general, the luminance value of each pixel is often used as the pixel value.

SADによる評価関数の値は、シティブロック距離と称されるものであり、小領域間の相関が高い程(類似している程)、シティブロック距離の値が小さくなり、シティブロック距離が最小値を取る小領域間の水平方向のシフト量により、視差が与えられる。シティブロック距離Cは、画像平面上の位置を、水平方向をi座標、垂直方向をj座標とする直交座標で定義し、互いの相関度を探索する小領域を、i×j(i=0〜n,j=0〜n)の探索ブロックとするとき、以下の(1)式に示すように、メイン画像の探索ブロックmain(i,j)とサブ画像の探索ブロックsub(i,j)とのSAD値を、i軸上を所定のシフト値ずつずらしながら演算することにより得られる。
C=Σ│main(i,j)−sub(i,j)│…(1)
尚、小領域(探索ブロック)間のステレオマッチングについての詳細は、先に本出願人によって提出された特開平5−114099号公報に詳述されている。
The value of the evaluation function by SAD is referred to as the city block distance. The higher the correlation between small regions (the more similar), the smaller the value of the city block distance, and the minimum city block distance. The parallax is given by the amount of shift in the horizontal direction between the subregions. The city block distance C is defined as an orthogonal coordinate in which the horizontal direction is i-coordinate and the vertical direction is j-coordinate, and a small area in which a correlation degree is searched is i × j (i = 0). ˜n, j = 0 to n), the main image search block main (i, j) and the sub image search block sub (i, j) as shown in the following equation (1): Are calculated by shifting the SAD value by a predetermined shift value on the i-axis.
C = Σ | main (i, j) -sub (i, j) | ... (1)
Details of stereo matching between small regions (search blocks) are described in detail in Japanese Patent Laid-Open No. 5-1114099 previously filed by the present applicant.

以上のSAD演算では、例えば4×4画素の探索ブロックを対象として、メイン画像の探索ブロック(以下、「メインブロック」と記載する)に対してサブ画像の探索ブロック(以下、「サブブロック」と記載する)の位置を水平線上を1画素ずつずらしながらSAD演算を行うと、図2に示すように、シティブロック距離が最小値C0となる点がメインブロックとサブブロックとの相関度が最も高い対応位置(一致点)として求められる。   In the above SAD calculation, for example, for a search block of 4 × 4 pixels, a search block of a sub image (hereinafter referred to as “sub block”) with respect to a search block of a main image (hereinafter referred to as “main block”). When the SAD calculation is performed while shifting the position of (described) on the horizontal line by one pixel at a time, as shown in FIG. 2, the point where the city block distance becomes the minimum value C0 has the highest degree of correlation between the main block and the sub block. It is obtained as a corresponding position (matching point).

この一致点におけるメインブロックとサブブロックとの1画素単位のずれ量(水平方向のメインブロックの位置imとサブブロックの位置isとの差)は、1画素単位の分解能を有する視差(ピクセル視差)を与える。このピクセル視差に基づく距離情報は、ステレオ法における三角測量の原理から、被写体までの距離が大きくなるにつれて分解能が低下するため、遠方の被写体までの距離を精度良く求めるためには、1画素単位以下の分解能を有するサブピクセルレベルの視差を求める必要がある。   The amount of shift in units of one pixel between the main block and the sub block at this coincidence point (difference between the position im of the main block in the horizontal direction and the position is of the sub block) is a parallax (pixel parallax) having a resolution of one pixel unit. give. The distance information based on the pixel parallax is less than one pixel unit in order to accurately determine the distance to a distant subject because the resolution decreases as the distance to the subject increases due to the principle of triangulation in the stereo method. It is necessary to obtain a sub-pixel level parallax having a resolution of.

サブピクセルレベルの視差を求める技術としては、例えば本出願人による特開2000−283753号公報に開示されている補間処理の技術を採用することができる。この技術によれば、シティブロック距離の分布が極小点周りに対称形となることを利用し、ピクセル視差を与える一致点の前後のシティブロック距離の差の大小関係から極小点の位置を特定し、直線近似によってサブピクセルレベルの視差を求めることができる。   As a technique for obtaining the sub-pixel level parallax, for example, an interpolation processing technique disclosed in Japanese Patent Application Laid-Open No. 2000-2833753 by the present applicant can be employed. This technique uses the fact that the distribution of the city block distance is symmetrical around the minimum point, and identifies the position of the minimum point from the magnitude relationship of the difference in the city block distance before and after the coincidence point giving the pixel parallax. Sub-pixel level parallax can be obtained by linear approximation.

すなわち、1画素のサイズが無限に小さいと仮定した場合、図3に示すように、画像平面(i,j)におけるシティブロック距離Cの分布は連続的になり、視差検出方向であるi方向(水平方向)の極小点について考えると、極小点周りに対称形で連続的に分布するとみなすことができる。従って、図4に示すように、座標isにおけるシティブロック距離C0の一致点(ピクセル単位でシティブロック距離が最小となる点)と、その前後の座標is−1,is+1でのシティブロック距離C1,C4の隣接点との3点に基づいて、縦軸に対して線対称となる2本の直線L1,L2を算出し、これらの直線の交点を求めることにより、サブピクセル成分を算出することができる。   That is, assuming that the size of one pixel is infinitely small, as shown in FIG. 3, the distribution of the city block distance C in the image plane (i, j) becomes continuous, and the i direction (the parallax detection direction) Considering the minimum points in the horizontal direction, it can be considered that the distribution is symmetrical and continuously distributed around the minimum points. Therefore, as shown in FIG. 4, the coincidence point of the city block distance C0 at the coordinate is (the point where the city block distance is minimum in pixel units) and the city block distance C1, It is possible to calculate the sub-pixel component by calculating two straight lines L1 and L2 that are symmetrical with respect to the vertical axis based on the three points with the adjacent point of C4, and obtaining the intersection of these straight lines. it can.

しかしながら、ピクセル単位の探索によって得られたシティブロック距離C0の一致点は、サブピクセルレベルで見ると、図4に破線で示す近傍のシティブロック距離C2,C3を誤差として含んで検出されてしまうことから、必ずしも極小点の最近傍に位置するとは言えず、ピクセル単位の離散的なデータを用いた3点補間でサブピクセル成分を求める技術では、シティブロック距離の真の極小点に対して近似精度が粗くなることは避けられない。   However, the coincidence point of the city block distance C0 obtained by the pixel unit search is detected at the subpixel level, including the neighboring city block distances C2 and C3 indicated by broken lines in FIG. 4 as errors. Therefore, it cannot be said that it is necessarily located closest to the local minimum point, and the technique for obtaining the subpixel component by three-point interpolation using discrete data in units of pixels approximates the true local minimum point of the city block distance. It is inevitable that the surface becomes rough.

従って、本形態においては、以下に詳述するように、解像度拡張部30でメイン画像及びサブ画像に対してピクセル間の補間を行って解像度を拡張し、この解像度を拡張した画像に対し、ステレオマッチング部35で上述のSAD演算によるステレオマッチングを行うことにより、図4に破線で示すサブピクセルレベルのシティブロック距離C2,C3を含むデータを求める。そして、補間処理部40で、これらのデータの最小値で定義される一致点と、その前後の隣接点との3点による補間処理を行うことにより、サブピクセルレベルの視差を精度良く求める。   Accordingly, in this embodiment, as will be described in detail below, the resolution expansion unit 30 performs inter-pixel interpolation on the main image and the sub image to expand the resolution, and the image with the expanded resolution is subjected to stereo. By performing stereo matching by the above-mentioned SAD calculation in the matching unit 35, data including city block distances C2 and C3 at sub-pixel levels indicated by broken lines in FIG. 4 is obtained. Then, the interpolation processing unit 40 performs an interpolation process with three points of the coincidence point defined by the minimum value of these data and the adjacent points before and after that, thereby obtaining the sub-pixel level parallax with high accuracy.

本形態においては、視差検出を主目的とすることから、解像度拡張処理を視差検出方向である水平方向についてのみ行い、水平方向の各画素間に仮想的な画素を挿入して水平方向の解像度を1/2ピクセル単位のレベルに拡張する。この場合、演算量が増大するが、水平方向と垂直方向との双方について解像度拡張を行っても良く、マッチング精度をより向上させることができる。   In this embodiment, since the main purpose is parallax detection, the resolution expansion process is performed only in the horizontal direction, which is the parallax detection direction, and virtual pixels are inserted between the horizontal pixels to reduce the horizontal resolution. Expand to 1/2 pixel level. In this case, although the amount of calculation increases, resolution expansion may be performed in both the horizontal direction and the vertical direction, and matching accuracy can be further improved.

また、本形態においては、ステレオマッチング部35では、当初、元画像のメインブロック及びサブブロックに対して通常の1画素単位のマッチングを行い、この1画素単位でマッチングしたブロックの解像度を解像度拡張部30で拡張し、この解像を拡張したインブロック及びサブブロックに対し、再度、ステレオマッチング部35でマッチングを行うようにしている。これは、予めピクセルレベルでの一致点を求めておくことで、その後のサブピクセルレベルの処理の演算量及び演算時間を軽減するためであるが、最初に、解像度拡張部30で元画像全体の解像度を拡張しておき、この解像度を拡張した画像のメインブロック及びサブブロックに対してステレオマッチング部35でマッチング処理を行うようにしても良い。   In this embodiment, the stereo matching unit 35 initially performs normal pixel-by-pixel matching on the main block and sub-block of the original image, and the resolution of the block matched in units of one pixel is set as a resolution expansion unit. The stereo matching unit 35 performs the matching again on the in-block and the sub-block which are expanded at 30 and the resolution is expanded. This is to obtain the coincidence point at the pixel level in advance so as to reduce the calculation amount and calculation time of the subsequent sub-pixel level processing. The resolution may be expanded and the matching processing may be performed by the stereo matching unit 35 on the main block and sub-block of the image whose resolution has been expanded.

一方、認識処理部50は、主としてソフトウエア処理を行うプロセッサにより構成されるものであり、元画像メモリ25の元画像データ及び距離画像メモリ45の距離分布データに基づいて各種認識処理を行う。認識処理部50による認識結果は、車両の車外状況認識システム、ITS、ロボットビジョン、侵入物監視システム、飛行体の高度計測システム等の各種制御系で利用される。   On the other hand, the recognition processing unit 50 is configured mainly by a processor that performs software processing, and performs various recognition processes based on the original image data in the original image memory 25 and the distance distribution data in the distance image memory 45. The recognition result by the recognition processing unit 50 is used in various control systems such as a vehicle exterior situation recognition system, ITS, robot vision, intruder monitoring system, and flying object altitude measurement system.

例えば、車両の車外状況認識システムやITSに適用される場合には、3次元の距離情報に基づいて、自車前方の道路形状や先行車等の認識、道路上の車両の走行状況の認識等を行う。また、ロボットビジョンに適用される場合には、3次元の距離情報に基づいて、移動経路や障害物の有無等を認識する処理を行い、侵入物監視システムに適用される場合には、建造物等に侵入する人物や物体の存在と位置を認識する処理を行う。更に、無人ヘリコプター等の移動体に搭載されて高度を計測する高度計測システムに適用される場合には、距離分布情報から地表面を平面と認識し、鉛直高度を求める処理を行う。   For example, when applied to a vehicle exterior situation recognition system or ITS, based on the three-dimensional distance information, recognition of the road shape in front of the host vehicle, the preceding vehicle, etc., recognition of the running condition of the vehicle on the road, etc. I do. In addition, when applied to robot vision, a process for recognizing the movement path and the presence or absence of obstacles is performed based on the three-dimensional distance information. A process for recognizing the presence and position of a person or an object that enters the user is performed. Furthermore, when applied to an altitude measurement system that is mounted on a moving body such as an unmanned helicopter and measures altitude, the ground surface is recognized as a plane from the distance distribution information, and processing for obtaining the vertical altitude is performed.

次に、画像の解像度拡張によるステレオマッチング及びサブピクセル処理の詳細について、図5のフローチャートを用いて説明する。   Next, details of stereo matching and subpixel processing by image resolution expansion will be described with reference to the flowchart of FIG.

図5のフローチャートは、解像度拡張部30,ステレオマッチング部35,補間処理部40において、時系列的に実行される(一部は、並行して実行される)処理の流れを示すものであり、先ず、ステップS1で、画像平面の左上を原点とする直交座標上で、メインブロックmain(i,j)とサブブロックsub(i,j)との1画素単位の通常のステレオマッチングを行う。メイン及びサブの探索ブロックの大きさは、例えば4×4(i=0〜3,j=0〜3)であり、メインブロックに対してサブブロックの位置を水平方向に1画素ずつシフトさせながらSAD演算を行い、シティブロック距離が最小値C0となる一致点を求めてピクセル視差Dを算出する(図2参照)。   The flowchart of FIG. 5 shows the flow of processing executed in time series (partially executed in parallel) in the resolution extension unit 30, the stereo matching unit 35, and the interpolation processing unit 40. First, in step S1, normal stereo matching in units of one pixel is performed between the main block main (i, j) and the sub-block sub (i, j) on the orthogonal coordinates with the upper left of the image plane as the origin. The size of the main and sub search blocks is, for example, 4 × 4 (i = 0 to 3, j = 0 to 3), and the position of the sub block is shifted by one pixel in the horizontal direction with respect to the main block. SAD calculation is performed to find a coincidence point where the city block distance is the minimum value C0, and the pixel parallax D is calculated (see FIG. 2).

次に、ステップS2へ進み、メイン画像及びサブ画像の各々について、ピクセル間補間によるマッチング領域の解像度拡張処理を行う。この解像度拡張処理は、メイン画像については、図6に示すように、4×4のピクセル単位のブロックmain(i,j)を、8×4の1/2ピクセル単位の拡張ブロックmain_ext(i',j)とし、サブ画像については、図7に示すように、メインブロックの4×4のブロックに対して探索のために水平方向に前後1画素ずつ広げた6×4のピクセル単位のブロックsub(i,j)を、12×4の1/2ピクセル単位のブロックsub_ext(i',j)に拡張する。   Next, the process proceeds to step S <b> 2, and the resolution expansion processing of the matching region by inter-pixel interpolation is performed for each of the main image and the sub image. As shown in FIG. 6, this resolution expansion process is performed on the main image with a block main (i, j) of 4 × 4 pixel units and an expanded block main_ext (i ′ of 8 × 4 pixel units). , j), and for the sub-image, as shown in FIG. 7, a block sub of 6 × 4 pixels, in which the 4 × 4 block of the main block is expanded by one pixel before and after in the horizontal direction for searching. (i, j) is expanded to a block sub_ext (i ′, j) of 12 × 4 1/2 pixel units.

図6,図7に斜線で示す拡張ブロックのi’=(2×i+1)番目のピクセルは、元画像の探索ブロックのi番目のピクセルと(i+1)番目のピクセルとの間に挿入される仮想的なピクセルである。この仮想的なピクセルのデータ(補間データ)は、周囲のピクセル情報に基づいて生成され、 本形態においては、以下の(2),(3)式に示すように、i番目のピクセルデータと水平方向に隣接する(i+1)番目のピクセルデータとの平均値をとって線形補間される。尚、このピクセル補間には、簡易的には最近傍法を採用することも必ずしも不可ではないが、演算能力及び演算時間に余裕がある場合には、3次補間法によるピクセル補間を行っても良い。
main_ext(i',j)=(main(i,j)+main(i+1,j))/2…(2)
sub_ext(i',j)=(sub(i,j)+sub(i+1,j))/2 …(3)
但し、i’=2×i+1
The i ′ = (2 × i + 1) -th pixel of the extension block indicated by hatching in FIGS. 6 and 7 is a virtual inserted between the i-th pixel and the (i + 1) -th pixel of the search block of the original image. Pixel. This virtual pixel data (interpolation data) is generated based on the surrounding pixel information. In this embodiment, as shown in the following equations (2) and (3), the i-th pixel data and the horizontal pixel data Linear interpolation is performed by taking the average value of the (i + 1) th pixel data adjacent in the direction. For this pixel interpolation, it is not always impossible to use the nearest neighbor method. However, if there is a margin in calculation capacity and calculation time, pixel interpolation by the cubic interpolation method may be performed. good.
main_ext (i ′, j) = (main (i, j) + main (i + 1, j)) / 2 (2)
sub_ext (i ′, j) = (sub (i, j) + sub (i + 1, j)) / 2 (3)
However, i ′ = 2 × i + 1

また、元画像の探索ブロックのi番目のピクセルに対応する拡張ブロックのi’=(2×i)番目のピクセルデータは、元画像のピクセルデータを使用することも可能であるが、本形態においては、以下の(3)式に示すように、拡張ブロックのi’=(2×i)番目のピクセルデータを、元画像の(i−1),i,(i+1)番目の3点のピクセルデータ、すなわち、自身のピクセルデータと、水平方向前後に隣接するピクセルデータとの加重平均値で置き換える。
main_ext(i',j)=(m1×main(i-1,j)+m2×main(i,j)
+m3×main(i+1,j))/(m1+m2+m3)…(4)
sub_ext(i',j)=(m1×sub(i-1,j)+m2×sub(i,j)
+m3×sub(i+1,j))/(m1+m2+m3) …(5)
但し、i’=2×i
m1,m2,m3;加重平均の重み
(例えば、m1=1,m2=4,m3=1)
Further, the pixel data of the original image can be used as the i ′ = (2 × i) -th pixel data of the extension block corresponding to the i-th pixel of the search block of the original image. As shown in the following equation (3), i ′ = (2 × i) -th pixel data of the extension block is converted into (i−1), i, (i + 1) -th three pixels of the original image. Data is replaced with a weighted average value of its own pixel data and pixel data adjacent in the horizontal direction.
main_ext (i ′, j) = (m1 × main (i−1, j) + m2 × main (i, j)
+ M3 × main (i + 1, j)) / (m1 + m2 + m3) (4)
sub_ext (i ′, j) = (m1 × sub (i−1, j) + m2 × sub (i, j)
+ M3 × sub (i + 1, j)) / (m1 + m2 + m3) (5)
However, i ′ = 2 × i
m1, m2, m3; weighted average weight
(For example, m1 = 1, m2 = 4, m3 = 1)

この元画像のピクセルデータの置き換えは、画像の加重平均処理であり、ミスマッチングを誘発しやすいノイズ成分を低減する。これらのピクセル間補間と加重平均処理とを併用することにより、マッチングの信頼性向上とばらつきの低減を図ることができる。   The replacement of the pixel data of the original image is a weighted average process of the image, and reduces noise components that easily induce mismatching. By using the inter-pixel interpolation and the weighted average process in combination, it is possible to improve the reliability of matching and reduce variations.

以上の解像度拡張処理を施した後は、ステップS3へ進み、8×4ブロックのSAD演算を、拡張前のブロックに対して水平方向に前後1画素ずつ、1/2ピクセル単位で5回のシフト数で実行し、5個のSAD値(シティブロック距離)を得る。そして、この5個のSAD値を互いに比較し、最小SAD値(シティブロック距離の最小値)をとるシフト値を求める。この最小SAD値をとるシフト値は、メイン拡張ブロックとサブ拡張ブロックとの間で1/2ピクセル単位の視差を与えるずれ量であり、ステップS1におけるマッチングでサブピクセル成分が0であるとすると、シフト値は2となる。   After performing the above resolution expansion processing, the process proceeds to step S3, and the 8 × 4 block SAD calculation is shifted five times in units of 1/2 pixel, one pixel before and after the pixel in the horizontal direction with respect to the block before expansion. Run with numbers to get 5 SAD values (city block distance). Then, these five SAD values are compared with each other to obtain a shift value that takes the minimum SAD value (the minimum value of the city block distance). The shift value that takes the minimum SAD value is a shift amount that gives a parallax in units of 1/2 pixel between the main extension block and the sub extension block, and if the sub pixel component is 0 in the matching in step S1, The shift value is 2.

尚、要求される視差検出精度が1/2ピクセル単位で良い場合には、ステップS2の処理で求めたサブピクセル成分とステップS1で求めたピクセル視差Dとから1/2ピクセル単位の視差を求め、処理を終了するようにしても良い。また、ステップS1における元画像の4×4ブロックのSAD演算に代えて、元画像全体の解像度を拡張した画像(例えば、元画像が512×200ピクセルの場合、1024×200ピクセルに画像全体を拡張)に対して4×4ブロックのSAD演算を行っても良く、直接、1/2ピクセル単位での視差を求めることができる。   If the required parallax detection accuracy is ½ pixel unit, the ½ pixel unit parallax is obtained from the sub-pixel component obtained in step S2 and the pixel parallax D obtained in step S1. The processing may be terminated. Further, instead of the SAD calculation of the 4 × 4 block of the original image in step S1, an image obtained by extending the resolution of the entire original image (for example, when the original image is 512 × 200 pixels, the entire image is expanded to 1024 × 200 pixels). ) May be subjected to SAD calculation of 4 × 4 blocks, and the parallax in units of 1/2 pixel can be directly obtained.

次に、ステップS3からステップS4へ進み、拡張ブロックのSAD演算で求めた最小SAD値をとるシフト値と、最小SAD値の前後のSAD値を与えるシフト値との3点補間により、拡張ブロックにおけるサブピクセルレベルの視差(サブピクセル視差)を求める。   Next, the process proceeds from step S3 to step S4, and the three-point interpolation between the shift value that takes the minimum SAD value obtained by the SAD calculation of the extension block and the shift value that gives the SAD values before and after the minimum SAD value, A sub-pixel level parallax (sub-pixel parallax) is obtained.

具体的には、ステップS3で求めた最小SAD値(シティブロック距離)を中心とするSAD値(シティブロック距離)に対し、シティブロック距離の分布が極小点周りに対称であることを考慮し、シティブロック距離の極小値が、最小SAD値の座標を中心として何れの側に存在するかを調べる。   Specifically, considering that the distribution of the city block distance is symmetric around the minimum point with respect to the SAD value (city block distance) centered on the minimum SAD value (city block distance) obtained in step S3, It is examined on which side the minimum value of the city block distance exists with respect to the coordinates of the minimum SAD value.

そして、最小SAD値をとる一致点周辺で、この一致点の後でSAD値が大きく変化するときには、極小点は、最小SAD値の座標位置と、1/2ピクセル前の座標位置との間にあると判断し、最小SAD値の点と次のSAD値の点とを通る傾きmの直線と、1/2ピクセル前の座標位置のSAD値の点を通る傾き−mの直線との交点を極小点の座標として求め、この極小点の座標位置と最小SAD値の座標位置との差に基いてサブピクセル視差Sが算出される。   When the SAD value changes greatly after the coincidence point around the coincidence point having the minimum SAD value, the local minimum point is between the coordinate position of the minimum SAD value and the coordinate position before 1/2 pixel. It is determined that there is an intersection of a straight line having a slope m passing through the point of the minimum SAD value and the next SAD value and a straight line having a slope of −m passing through the point of the SAD value at the coordinate position 1/2 pixel before. The subpixel parallax S is calculated based on the difference between the coordinate position of the minimum point and the coordinate position of the minimum SAD value.

逆に、最小SAD値をとる一致点の前でSAD値が大きく変化するときには、極小点は、最小SAD値の座標位置と、1/2ピクセル後の座標位置との間にあると判断し、最小SAD値の点と前のSAD値の点とを通る傾きmの直線と、1/2ピクセル後の座標位置のSAD値の点を通る傾き−mの直線との交点を極小点の座標として求め、この極小点の座標位置と最小SAD値の座標位置との差に基いてサブピクセル視差Sが算出される。   Conversely, when the SAD value changes greatly before the coincidence point having the minimum SAD value, it is determined that the minimum point is between the coordinate position of the minimum SAD value and the coordinate position after 1/2 pixel, The intersection of the straight line with the slope m passing through the point with the minimum SAD value and the point with the previous SAD value and the straight line with the slope -m passing through the point with the SAD value at the coordinate position after 1/2 pixel is used as the coordinate of the minimum point. The sub-pixel parallax S is calculated based on the difference between the coordinate position of the minimum point and the coordinate position of the minimum SAD value.

例えば、図8に示すように、1/2ピクセル単位の各座標(i's−2),(i's−1),(i's),(i's+1),(i's+2)における5個のSAD値C'1,C'2,C'0,C'3,C'4のうち、最小となるSAD値がC'3であり、(C'0−C'3)<(C'4−C'3)である場合には、SAD値C'3の点と1/2ピクセル後のSAD値C'4の点を通る直線L'1を求め、1/2ピクセル前のSAD値C'0を通って直線L'1と縦軸に対して線対称となる直線L'2との交点を求め、この交点のi'座標と、最小SAD値C'3の座標(i's+1)との差からサブピクセル視差Sが算出される。   For example, as shown in FIG. 8, five SAD values at coordinates (i′s−2), (i′s−1), (i ′s), (i ′s + 1), and (i ′s + 2) in 1/2 pixel units. Among C′1, C′2, C′0, C′3, and C′4, the minimum SAD value is C′3, and (C′0−C′3) <(C′4−C In the case of '3), a straight line L'1 passing through the point of the SAD value C'3 and the point of the SAD value C'4 after 1/2 pixel is obtained, and the SAD value C'0 before 1/2 pixel is obtained. The intersection of the straight line L′ 1 and the straight line L′ 2 that is symmetric with respect to the vertical axis is obtained, and the i ′ coordinate of this intersection and the coordinate (i ′s + 1) of the minimum SAD value C′3 are obtained. A sub-pixel parallax S is calculated from the difference.

図8の例は、先のピクセル単位のマッチングで最小値として検出されたSAD値C0の座標位置isに対応する座標位置i'sでのSAD値C'0に対し、解像度拡張処理後のサブピクセル単位のマッチングにより、1/2ピクセル後の座標位置(i's+1)におけるSAD値C'3が最小値として検出された例を示すものである。この例からも明らかなように、ピクセル単位のマッチング後に行う3点補間では、SAD値C'0に該当する点を中心として、図中の破線の直線で示すように、SAD値C'1に該当する点、SAD値C'4に該当する点の3点補間によって極小点を推定することになる。従って、解像度拡張後のマッチングによって得られたデータを用いて3点補間を行うことにより、ピクセル単位でのデータを用いた3点補間に比較し、真の極小点により近い位置を求めることができ、誤差を小さくして近似精度を大幅に向上することができる。   In the example of FIG. 8, the SAD value C′0 at the coordinate position i ′s corresponding to the coordinate position is of the SAD value C0 detected as the minimum value in the previous pixel-unit matching is sub-pixel unit after resolution expansion processing. This shows an example in which the SAD value C′3 at the coordinate position (i ′s + 1) after 1/2 pixel is detected as the minimum value by the above matching. As is clear from this example, in the three-point interpolation performed after pixel-by-pixel matching, the SAD value C′1 is centered on the point corresponding to the SAD value C′0, as indicated by the dashed straight line in the figure. The minimum point is estimated by three-point interpolation of the corresponding point and the point corresponding to the SAD value C′4. Therefore, by performing three-point interpolation using data obtained by matching after resolution expansion, a position closer to the true minimum point can be obtained compared to three-point interpolation using data in pixel units. The error can be reduced and the approximation accuracy can be greatly improved.

サブピクセル視差Sを算出した後は、ステップS4からステップ5へ進み、ピクセル視差Dとサブピクセル視差Sとから、1/2ピクセル単位以下の分解能の視差(D+S)を得る。このような処理を繰り返して画像全体に渡って視差(D+S)を求め、距離画像メモリ45に格納する。この距離画像メモリ45に格納された視差(D+S)からは、実空間座標への変換により実空間での距離を求めることができる。   After calculating the sub-pixel parallax S, the process proceeds from step S4 to step 5, and a parallax (D + S) having a resolution of 1/2 pixel unit or less is obtained from the pixel parallax D and the sub-pixel parallax S. By repeating such processing, the parallax (D + S) is obtained over the entire image and stored in the distance image memory 45. From the parallax (D + S) stored in the distance image memory 45, the distance in the real space can be obtained by conversion into the real space coordinates.

図9は、以上の処理によって算出される視差のばらつきを示すヒストグラムであり、本形態の解像度拡張によるサブピクセル単位のマッチング後の3点補間によって得られた視差の真値に対するばらつき分布が白抜きのヒストグラム値で示されている。このヒストグラムに示すように、本形態の処理によって得られる視差は、図9中に塗りつぶしのヒストグラム値で示す従来のピクセル単位のマッチング後の三点補間に比較し、大幅にばらつきが減少していることがわかる。これは、本形態の処理では、元画像の解像度を拡張していることから、サブピクセル算出のための3点補間の近似性が向上し、また、マッチングピクセル数の増加によって安定性が向上しているためであり、更には、加重平均処理を加えることによるミスマッチングの防止が有効に作用しているためであると考えられる。   FIG. 9 is a histogram showing the parallax variation calculated by the above processing, and the variation distribution with respect to the true value of the parallax obtained by the three-point interpolation after matching in sub-pixel units by the resolution extension of this embodiment is outlined. The histogram value is shown. As shown in this histogram, the disparity obtained by the processing of this embodiment is greatly reduced in variation compared to the conventional three-point interpolation after matching in units of pixels indicated by the filled histogram values in FIG. I understand that. This is because, in the processing of this embodiment, the resolution of the original image is expanded, so that the approximation of the three-point interpolation for subpixel calculation is improved, and the stability is improved by increasing the number of matching pixels. Furthermore, it is considered that the prevention of mismatching by adding a weighted average process is effective.

以上のように本実施の形態においては、ステレオカメラで撮像した元画像の解像度を、ピクセル間補間と加重平均処理とを併用して拡張し、拡張した画像に対してステレオマッチングを行うことで、マッチング精度を向上して1/2ピクセル単位の視差を精度良く検出することが可能となる。更に、解像度を拡張した画像をステレオマッチングした後、3点補間を行うことにより、1/2ピクセル単位以下の分解能で精度良く視差を求めることができる。   As described above, in the present embodiment, the resolution of the original image captured by the stereo camera is expanded by using inter-pixel interpolation and weighted average processing, and stereo matching is performed on the expanded image. Matching accuracy can be improved and parallax in units of 1/2 pixel can be accurately detected. Furthermore, parallax can be obtained accurately with a resolution of 1/2 pixel unit or less by performing three-point interpolation after stereo matching of an image with an extended resolution.

これにより、ピクセル単位でのマッチング領域を大きくして空間分解能の低下を招くことなく、サブピクセルレベルの視差データをノイズやばらつきを低減して信頼性高く得ることができ、距離分解能を向上した緻密な距離画像を得て、より高度な画像認識及び制御に繋げることが可能となる。   This makes it possible to obtain subpixel-level parallax data with high reliability by reducing noise and variations without increasing the matching area in pixels and reducing spatial resolution. It is possible to obtain an accurate distance image and to lead to more advanced image recognition and control.

次に、本発明の実施の第2形態について説明する。図10〜図12は本発明の実施の第2形態に係り、図10はステレオ画像処理装置の基本構成を示すブロック図、図11はステレオマッチング及びサブピクセル処理のフローチャート、図12は視差のばらつきを示すヒストグラムである。   Next, a second embodiment of the present invention will be described. 10 to 12 relate to a second embodiment of the present invention, FIG. 10 is a block diagram showing a basic configuration of a stereo image processing apparatus, FIG. 11 is a flowchart of stereo matching and subpixel processing, and FIG. It is a histogram which shows.

第2形態は、ステレオカメラで撮像した元画像を加重平均処理した後、ステレオマッチングを行うものである。以下、第1形態と同様の機能を有する構成要素については、第1形態と同じ符号を付し、必要に応じて第1形態の説明を援用するものとする。   In the second form, stereo matching is performed after weighted averaging of an original image captured by a stereo camera. Hereinafter, the same reference numerals as those of the first embodiment are given to the constituent elements having the same functions as those of the first embodiment, and the description of the first embodiment is incorporated as necessary.

すなわち、第2形態のステレオ画像処理装置1Aは、2台1組のカメラから構成されるステレオカメラ10、このステレオカメラ10で撮像した一対の画像に対する入力処理を行う画像入力部20、この画像入力部で処理した画像を元画像としてストアする元画像メモリ25、元画像の各画素のデータを周辺画素のデータを用いて置き換える画像変換部30A、一対の画像間の相関度を評価して互いの対応位置を決定するステレオマッチング部35、ステレオマッチング部35で決定した対応位置に基づく補間処理により1画素単位以下の分解能を有するサブピクセルレベルの視差を求める補間処理部40、視差に基づく距離情報を元画像の画像座標と対応付けて二次元的に配置した画像形態(距離画像)としてストアする距離画像メモリ45、元画像メモリ25の元画像情報及び距離画像メモリ45の距離情報を用いて各種認識処理を行う認識処理部50を備えて構成され、車両の車外状況認識システム、高度道路交通システム(ITS)、ロボットビジョン、侵入物監視システム、飛行体の対地高度計測システム等の各種システムに適用される。   That is, the stereo image processing apparatus 1A according to the second embodiment includes a stereo camera 10 including a set of two cameras, an image input unit 20 that performs input processing on a pair of images captured by the stereo camera 10, and the image input. The original image memory 25 that stores the image processed by the image processing unit as an original image, the image conversion unit 30A that replaces the data of each pixel of the original image with the data of the surrounding pixels, and evaluates the degree of correlation between the pair of images. A stereo matching unit 35 that determines a corresponding position, an interpolation processing unit 40 that obtains a sub-pixel level parallax having a resolution of one pixel unit or less by an interpolation process based on the corresponding position determined by the stereo matching unit 35, and distance information based on the parallax. Distance image memory 4 for storing as an image form (distance image) arranged two-dimensionally in association with the image coordinates of the original image The recognition processing unit 50 is configured to perform various recognition processes using the original image information in the original image memory 25 and the distance information in the distance image memory 45, and includes a vehicle exterior situation recognition system, an intelligent road traffic system (ITS), It is applied to various systems such as robot vision, intruder monitoring system, and flying object ground altitude measurement system.

画像変換部30Aを除く各部の機能は、前述の第1形態と同様であるが、第2形態においては、ステレオカメラ10で撮像した一対の元画像の各画素のデータを、画像変換部30Aで周辺画素のデータを用いて置き換え、画素のデータを置き換えた一対の画像に対して、ステレオマッチング部35でステレオマッチングを行う。   The functions of the units other than the image conversion unit 30A are the same as those in the first embodiment described above. However, in the second embodiment, the data of each pixel of the pair of original images captured by the stereo camera 10 is received by the image conversion unit 30A. Stereo matching is performed by the stereo matching unit 35 on the pair of images in which the pixel data is replaced using the peripheral pixel data.

画像変換部30Aによる画素データの置き換えは、第1形態で説明した画像の加重平均処理であり、例えば、(4),(5)式に準じる以下の(6)式に示すように、元画像のi番目の画素のデータp(i,j)を、自身のデータ及び隣接する周辺の(i−1)番目の画素のデータp(i-1,j)、(i+1)番目の画素のデータp(i+1,j)で加重平均した新たなデータext(i,j)で置き換える。
ext(i,j)=(m1×p(i-1,j)+m2×p(i,j)
+m3×p(i+1,j))/(m1+m2+m3)…(6)
但し、m1,m2,m3;加重平均の重み
(例えば、m1=1,m2=4,m3=1)
The replacement of the pixel data by the image conversion unit 30A is the weighted average processing of the image described in the first embodiment. For example, as shown in the following equation (6) according to the equations (4) and (5), the original image The data p (i, j) of the i-th pixel of the current data, the data p (i−1, j) of the adjacent (i−1) -th pixel, and the data of the (i + 1) -th pixel Replace with new data ext (i, j) weighted averaged by p (i + 1, j).
ext (i, j) = (m1 * p (i-1, j) + m2 * p (i, j)
+ M3 × p (i + 1, j)) / (m1 + m2 + m3) (6)
Where m1, m2, m3: weighted average weight
(For example, m1 = 1, m2 = 4, m3 = 1)

加重平均処理した画像に対してステレオマッチングを行った後の処理は、第1形態と同様であり、補間処理部40における3点補間を経て1画素単位以下の分解能を有するサブピクセルレベルの視差が得られる。   The processing after performing the stereo matching on the weighted average processed image is the same as that in the first embodiment, and the sub-pixel level parallax having a resolution of one pixel unit or less through the three-point interpolation in the interpolation processing unit 40. can get.

以下、画素データの置き換えによるステレオマッチング及びサブピクセル処理について、図11のフローチャートを用いて説明する。   Hereinafter, stereo matching and sub-pixel processing by pixel data replacement will be described with reference to the flowchart of FIG.

図11のフローチャートは、画像変換部30A,ステレオマッチング部35,補間処理部40において、時系列的に実行される(一部は、並行して実行される)処理の流れを示すものであり、先ず、ステップS11で、メイン画像及びサブ画像の双方の元画像に対し、各画素のデータを、自身の画素データと隣接する周辺の画素のデータとで置き換え、上述の加重平均処理を行う。この加重平均処理は、最初に画像全体に渡って行っても良く、或いは、SAD計算を行うマッチング領域毎に行っても良い。   The flowchart of FIG. 11 shows a flow of processing executed in time series (partially executed in parallel) in the image conversion unit 30A, the stereo matching unit 35, and the interpolation processing unit 40. First, in step S11, with respect to the original images of both the main image and the sub image, the data of each pixel is replaced with its own pixel data and data of neighboring pixels adjacent thereto, and the above-described weighted average processing is performed. This weighted average processing may be performed over the entire image first, or may be performed for each matching region where SAD calculation is performed.

次に、ステップS12へ進み、画像平面の左上を原点とする直交座標上で、メインブロックmain(i,j)とサブブロックsub(i,j)との1画素単位の通常のステレオマッチングを行う。メイン及びサブの探索ブロックの大きさは、例えば4×4(i=0〜3,j=0〜3)であり、メインブロックに対してサブブロックの位置を水平方向に1画素ずつシフトさせながらSAD演算を行い、シティブロック距離が最小値C0となる一致点を求めてピクセル単位の視差(ピクセル視差)Dを算出する(図2参照)。   Next, the process proceeds to step S12, and normal stereo matching in units of one pixel is performed between the main block main (i, j) and the sub-block sub (i, j) on the orthogonal coordinates having the origin at the upper left of the image plane. . The size of the main and sub search blocks is, for example, 4 × 4 (i = 0 to 3, j = 0 to 3), and the position of the sub block is shifted by one pixel in the horizontal direction with respect to the main block. A SAD operation is performed to obtain a coincidence point where the city block distance is the minimum value C0, and a parallax (pixel parallax) D in units of pixels is calculated (see FIG. 2).

次に、ステップS12からステップS13へ進み、ステップS12のSAD演算で求めた最小SAD値をとるシフト値と、最小SAD値の前後のSAD値を与えるシフト値との3点補間により、サブピクセルレベルの視差(サブピクセル視差)Sを求める。この3点補間による処理は、第1形態で説明した通りである(図3、図4参照)。   Next, the process proceeds from step S12 to step S13, and a subpixel level is obtained by three-point interpolation between a shift value that takes the minimum SAD value obtained by the SAD calculation in step S12 and a shift value that gives SAD values before and after the minimum SAD value. The parallax (subpixel parallax) S is obtained. The processing by this three-point interpolation is as described in the first embodiment (see FIGS. 3 and 4).

サブピクセル視差Sを算出した後は、ステップS13からステップ14へ進み、ピクセル視差Dとサブピクセル視差Sとから、1画素単位以下の分解能の視差(D+S)を得る。このような処理を繰り返して画像全体に渡って視差(D+S)を求め、距離画像メモリ45に格納する。この距離画像メモリ45に格納された視差(D+S)からは、実空間座標への変換により実空間での距離を求めることができる。   After calculating the sub-pixel parallax S, the process proceeds from step S 13 to step 14, and a parallax (D + S) having a resolution of one pixel unit or less is obtained from the pixel parallax D and the sub-pixel parallax S. By repeating such processing, the parallax (D + S) is obtained over the entire image and stored in the distance image memory 45. From the parallax (D + S) stored in the distance image memory 45, the distance in the real space can be obtained by conversion into the real space coordinates.

以上の加重平均処理を経てステレオマッチングした後の3点補間による視差は、図12に示すように、真値に対して、図中に白抜きのヒストグラム値で示すばらつき分布となっている。このヒストグラムから明らかなように、第2形態においては、第1形態に比較すると、元画像の解像度を拡張していないことから、3点補間の近似性向上やマッチングピクセル数の増加による安定化が期待できない分、ばらつきの改善度は小さいものの、図12中に塗りつぶしのヒストグラム値で示す、元画像を加重平均処理することなくマッチングした後の三点補間によって得られる従来の視差に比較し、良好にばらつきが減少していることが確認できる。   As shown in FIG. 12, the parallax by the three-point interpolation after the stereo matching through the above weighted average processing has a variation distribution indicated by a white histogram value in the figure with respect to the true value. As is apparent from this histogram, in the second form, compared to the first form, the resolution of the original image is not expanded, so that the approximation by the three-point interpolation is improved and the stabilization by the increase in the number of matching pixels is achieved. Although the degree of improvement in variation is small because it cannot be expected, it is better than the conventional parallax obtained by three-point interpolation after matching the original image without performing weighted average processing, as shown by the histogram value of the fill in FIG. It can be confirmed that the variation is reduced.

以上のように第2形態においては、ステレオカメラで撮像した元画像をぼかし処理した後にステレオマッチングを行うことで、ミスマッチングを減らしてマッチング精度及び信頼性を高め、ばらつきの低減を図ることができる。しかも、ピクセル単位のマッチングにおいて、いたずらにマッチング領域を拡大して距離分解能を低下させることもなく、補間処理によりサブピクセルレベルの視差データをノイズやばらつきを低減して信頼性高く得ることができ、より高度な画像認識及び制御に繋げることが可能となる。   As described above, in the second embodiment, stereo matching is performed after blurring the original image captured by the stereo camera, so that mismatching can be reduced, matching accuracy and reliability can be increased, and variation can be reduced. . Moreover, in matching in units of pixels, it is possible to obtain sub-pixel level parallax data with high reliability by reducing noise and variations without interfering with enlarging the matching area and reducing the distance resolution. It becomes possible to connect to more advanced image recognition and control.

尚、上述した実施の形態には限定されず、一方のみの画像に対して解像度を拡張し、例えば1/2ピクセルにてマッチング処理を行うことも可能である。   Note that the present invention is not limited to the above-described embodiment, and it is possible to extend the resolution of only one image and perform matching processing with, for example, 1/2 pixel.

また、一方の画像に対して加重平均処理を行い、他方の画像の解像度を拡張し、マッチングを行うことも可能である。この場合、他方の画像に対して加重平均処理を施しても、施さなくとも良い。   It is also possible to perform weighted average processing on one image, expand the resolution of the other image, and perform matching. In this case, the weighted average process may or may not be performed on the other image.

更に、撮像手段としてカラー対応型のカメラを用いた場合において、単色画像を抜きだし、この画像に対して本発明を適用することも可能である。この場合には、カラー対応型のカメラにおいて課題となるマッチングの精度を向上させることができる。   Furthermore, when a color-compatible camera is used as the imaging means, it is also possible to extract a monochrome image and apply the present invention to this image. In this case, it is possible to improve the accuracy of matching that is a problem in a color-compatible camera.

本発明の実施の第1形態に係り、ステレオ画像処理装置の基本構成を示すブロック図The block diagram which shows the basic composition of a stereo image processing apparatus concerning 1st Embodiment of this invention. 同上、メイン画像とサブ画像とのマッチングを示す説明図Same as above, explanatory diagram showing matching between main image and sub-image 同上、シティブロック距離の分布を示す説明図Same as above, explanatory diagram showing distribution of city block distance 同上、ピクセル単位のマッチング後の補間を示す説明図Same as above, explanatory diagram showing interpolation after matching in pixel units 同上、ステレオマッチング及びサブピクセル処理のフローチャートSame as above, stereo matching and sub-pixel processing flowchart 同上、メイン画像の解像度拡張を示す説明図As above, an explanatory diagram showing the resolution expansion of the main image 同上、サブ画像の解像度拡張を示す説明図Same as above, explanatory diagram showing resolution expansion of sub-image 同上、サブピクセル単位のマッチング後の補間を示す説明図Same as above, explanatory diagram showing interpolation after matching in sub-pixel units 同上、視差のばらつきを示すヒストグラムSame as above, histogram showing disparity variation 本発明の実施の第2形態に係り、ステレオ画像処理装置の基本構成を示すブロック図The block diagram which shows the basic composition of a stereo image processing apparatus concerning 2nd Embodiment of this invention. 同上、ステレオマッチング及びサブピクセル処理のフローチャートSame as above, stereo matching and sub-pixel processing flowchart 同上、視差のばらつきを示すヒストグラムSame as above, histogram showing disparity variation

符号の説明Explanation of symbols

1,1A ステレオ画像処理装置
10 ステレオカメラ
20 画像入力部
25 元画像メモリ
30 解像度拡張部
30A 画像変換部
35 ステレオマッチング部
40 補間処理部
代理人 弁理士 伊 藤 進
DESCRIPTION OF SYMBOLS 1,1A Stereo image processing apparatus 10 Stereo camera 20 Image input part 25 Original image memory 30 Resolution expansion part 30A Image conversion part 35 Stereo matching part 40 Interpolation processing part
Agent Patent Attorney Susumu Ito

Claims (6)

互いに相関を有する画像対を用いてステレオマッチングを行うステレオ画像処理装置であって、
上記画像対の各画素間に、周辺画素のデータを用いて生成した仮想的な画素を挿入し、上記画像対の解像度を拡張する解像度拡張手段と、
上記解像度拡張手段で解像度を拡張した画像対に対し、拡張した解像度の分解能で互いの対応位置を特定するステレオマッチング手段とを備えたことを特徴とするステレオ画像処理装置。
A stereo image processing apparatus for performing stereo matching using image pairs having correlation with each other,
A resolution extending means for inserting a virtual pixel generated using data of peripheral pixels between each pixel of the image pair, and extending the resolution of the image pair;
A stereo image processing apparatus, comprising: stereo matching means for specifying a corresponding position with the resolution of the expanded resolution for the image pair whose resolution is expanded by the resolution expansion means.
上記解像度拡張手段は、
上記仮想的な画素に隣接する画素のデータを、周辺画素のデータを用いて置き換えることを特徴とする請求項1記載のステレオ画像処理装置。
The resolution expansion means is
2. The stereo image processing apparatus according to claim 1, wherein data of a pixel adjacent to the virtual pixel is replaced by using data of a peripheral pixel.
上記ステレオマッチング手段で特定した対応位置の相関度演算値と、この対応位置に隣接する位置の相関度演算値とに基づいて補間処理を行い、拡張した解像度以下の分解能で対応位置を決定する補間手段を備えたことを特徴とする請求項1又は2記載のステレオ画像処理装置。   Interpolation based on the correlation calculation value of the corresponding position specified by the stereo matching means and the correlation calculation value of the position adjacent to the corresponding position, and determining the corresponding position with a resolution lower than the expanded resolution 3. The stereo image processing apparatus according to claim 1, further comprising means. 互いに相関を有する画像対を用いてステレオマッチングを行うステレオ画像処理装置であって、
上記画像対の各画素のデータを、周辺画素のデータを用いて置き換える画像変換手段と、
上記画像変換手段で画素のデータを置き換えた画像対に対し、互いの対応位置を特定するステレオマッチング手段とを備えたことを特徴とするステレオ画像処理装置。
A stereo image processing apparatus for performing stereo matching using image pairs having correlation with each other,
Image conversion means for replacing data of each pixel of the image pair with data of peripheral pixels;
A stereo image processing apparatus, comprising: stereo matching means for specifying a corresponding position with respect to an image pair in which pixel data is replaced by the image conversion means.
上記ステレオマッチング手段で特定した対応位置の相関度演算値と、この対応位置に隣接する位置の相関度演算値とに基づいて補間処理を行い、1画素以下の分解能で対応位置を決定する補間手段を備えたことを特徴とする請求項4記載のステレオ画像処理装置。   Interpolation means for performing an interpolation process based on the correlation value calculation value of the corresponding position specified by the stereo matching means and the correlation value calculation value of the position adjacent to the corresponding position, and determining the corresponding position with a resolution of one pixel or less. The stereo image processing apparatus according to claim 4, further comprising: 互いに相関を有する画像対を用いてステレオマッチングを行うステレオ画像処理装置であって、
上記画像対のうち、少なくとも一方の画像に対して、各画素のデータを周辺画素のデータを用いて置き換える手段と、
置き換えられた画素データを用いて、上記画像対における互いの対応位置を特定するステレオマッチング手段とを備えたことを特徴とするステレオ画像処理装置。
A stereo image processing apparatus for performing stereo matching using image pairs having correlation with each other,
Means for replacing the data of each pixel with the data of surrounding pixels for at least one of the image pairs;
A stereo image processing apparatus, comprising: stereo matching means for specifying mutual corresponding positions in the image pair using the replaced pixel data.
JP2004062661A 2004-03-05 2004-03-05 Stereo image processing device Expired - Fee Related JP4394487B2 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2004062661A JP4394487B2 (en) 2004-03-05 2004-03-05 Stereo image processing device

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2004062661A JP4394487B2 (en) 2004-03-05 2004-03-05 Stereo image processing device

Publications (2)

Publication Number Publication Date
JP2005250994A true JP2005250994A (en) 2005-09-15
JP4394487B2 JP4394487B2 (en) 2010-01-06

Family

ID=35031401

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2004062661A Expired - Fee Related JP4394487B2 (en) 2004-03-05 2004-03-05 Stereo image processing device

Country Status (1)

Country Link
JP (1) JP4394487B2 (en)

Cited By (13)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2007066842A1 (en) * 2005-12-07 2007-06-14 Electronics And Telecommunications Research Institute Apparatus and method for vision processing on network based intelligent service robot system and the system using the same
JP2007249682A (en) * 2006-03-16 2007-09-27 Fuji Heavy Ind Ltd Correlation evaluation system and correlation evaluation method
JP2009146296A (en) * 2007-12-17 2009-07-02 Konica Minolta Holdings Inc Image corresponding point search device, and distance measuring equipment and image motion detector using the same
JP2011185720A (en) * 2010-03-08 2011-09-22 Ricoh Co Ltd Distance obtaining device
JP2011191905A (en) * 2010-03-12 2011-09-29 Hitachi Automotive Systems Ltd Stereo camera device
JP2013500512A (en) * 2009-07-24 2013-01-07 ピルツ ゲーエムベーハー アンド コー.カーゲー Method and apparatus for monitoring a spatial domain
WO2013175700A1 (en) * 2012-05-22 2013-11-28 株式会社ソニー・コンピュータエンタテインメント Information processing device and information processing method
JP2017142829A (en) * 2017-03-28 2017-08-17 エスゼット ディージェイアイ テクノロジー カンパニー リミテッドSz Dji Technology Co.,Ltd Method, device and program for generating depth map
EP3211368A4 (en) * 2014-10-24 2018-06-06 Kyocera Corporation Stereo camera apparatus and vehicle provided with stereo camera apparatus
US10198830B2 (en) 2013-12-26 2019-02-05 Ricoh Company, Ltd. Parallax operation system, information processing apparatus, information processing method, and recording medium
JP2019120591A (en) * 2018-01-05 2019-07-22 パナソニックIpマネジメント株式会社 Parallax value calculation device, parallax value calculation method and program
CN110298872A (en) * 2019-07-03 2019-10-01 云南电网有限责任公司电力科学研究院 A kind of method for registering of ultraviolet light camera and Visible Light Camera array
US10475209B2 (en) 2014-11-04 2019-11-12 SZ DJI Technology Co., Ltd. Camera calibration

Cited By (18)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2007066842A1 (en) * 2005-12-07 2007-06-14 Electronics And Telecommunications Research Institute Apparatus and method for vision processing on network based intelligent service robot system and the system using the same
US8170324B2 (en) 2005-12-07 2012-05-01 Electronics And Telecommunications Research Institute Apparatus and method for vision processing on network based intelligent service robot system and the system using the same
JP2007249682A (en) * 2006-03-16 2007-09-27 Fuji Heavy Ind Ltd Correlation evaluation system and correlation evaluation method
JP2009146296A (en) * 2007-12-17 2009-07-02 Konica Minolta Holdings Inc Image corresponding point search device, and distance measuring equipment and image motion detector using the same
JP2013500512A (en) * 2009-07-24 2013-01-07 ピルツ ゲーエムベーハー アンド コー.カーゲー Method and apparatus for monitoring a spatial domain
US9292924B2 (en) 2009-07-24 2016-03-22 Pilz Gmbh & Co. Kg Method and device for monitoring a spatial region
JP2011185720A (en) * 2010-03-08 2011-09-22 Ricoh Co Ltd Distance obtaining device
JP2011191905A (en) * 2010-03-12 2011-09-29 Hitachi Automotive Systems Ltd Stereo camera device
JP2013242812A (en) * 2012-05-22 2013-12-05 Sony Computer Entertainment Inc Information processing device and information processing method
WO2013175700A1 (en) * 2012-05-22 2013-11-28 株式会社ソニー・コンピュータエンタテインメント Information processing device and information processing method
US10469829B2 (en) 2012-05-22 2019-11-05 Sony Interactive Entertainment Inc. Information processor and information processing method
US10198830B2 (en) 2013-12-26 2019-02-05 Ricoh Company, Ltd. Parallax operation system, information processing apparatus, information processing method, and recording medium
EP3211368A4 (en) * 2014-10-24 2018-06-06 Kyocera Corporation Stereo camera apparatus and vehicle provided with stereo camera apparatus
US10595003B2 (en) 2014-10-24 2020-03-17 Kyocera Corporation Stereo camera apparatus and vehicle comprising the same
US10475209B2 (en) 2014-11-04 2019-11-12 SZ DJI Technology Co., Ltd. Camera calibration
JP2017142829A (en) * 2017-03-28 2017-08-17 エスゼット ディージェイアイ テクノロジー カンパニー リミテッドSz Dji Technology Co.,Ltd Method, device and program for generating depth map
JP2019120591A (en) * 2018-01-05 2019-07-22 パナソニックIpマネジメント株式会社 Parallax value calculation device, parallax value calculation method and program
CN110298872A (en) * 2019-07-03 2019-10-01 云南电网有限责任公司电力科学研究院 A kind of method for registering of ultraviolet light camera and Visible Light Camera array

Also Published As

Publication number Publication date
JP4394487B2 (en) 2010-01-06

Similar Documents

Publication Publication Date Title
CN107077743B (en) System and method for dynamic calibration of an array camera
KR100513055B1 (en) 3D scene model generation apparatus and method through the fusion of disparity map and depth map
KR100776649B1 (en) A depth information-based Stereo/Multi-view Stereo Image Matching Apparatus and Method
US9300946B2 (en) System and method for generating a depth map and fusing images from a camera array
US10275857B2 (en) Depth image enhancement method
US9374571B2 (en) Image processing device, imaging device, and image processing method
JP4554316B2 (en) Stereo image processing device
WO2020007320A1 (en) Method for fusing multi-visual angle images, apparatus, computer device, and storage medium
KR100938195B1 (en) Method for distance estimation and apparatus for the same using a stereo matching
CN113034568B (en) Machine vision depth estimation method, device and system
CN104685513A (en) Feature based high resolution motion estimation from low resolution images captured using an array source
JP4394487B2 (en) Stereo image processing device
KR20120068470A (en) Apparatus for matching stereo image and method thereof
US20130272600A1 (en) Range image pixel matching method
US8867826B2 (en) Disparity estimation for misaligned stereo image pairs
CN107545586B (en) Depth obtaining method and system based on light field polar line plane image local part
US8340399B2 (en) Method for determining a depth map from images, device for determining a depth map
KR101875532B1 (en) Apparatus and method for hierarchical stereo matching
JP3008875B2 (en) Subject extraction device
Garro et al. A novel interpolation scheme for range data with side information
TWI571099B (en) Device and method for depth estimation
US9674503B2 (en) Stereo matching apparatus using image property
Um et al. Three-dimensional scene reconstruction using multiview images and depth camera
US10726528B2 (en) Image processing apparatus and image processing method for image picked up by two cameras
US9380285B2 (en) Stereo image processing method, stereo image processing device and display device

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20061102

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20090707

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20090907

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

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

R150 Certificate of patent or registration of utility model

Ref document number: 4394487

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150

Free format text: JAPANESE INTERMEDIATE CODE: R150

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

Free format text: PAYMENT UNTIL: 20121023

Year of fee payment: 3

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

Free format text: PAYMENT UNTIL: 20121023

Year of fee payment: 3

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

Free format text: PAYMENT UNTIL: 20131023

Year of fee payment: 4

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

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

LAPS Cancellation because of no payment of annual fees