JP2014067269A - Detector - Google Patents
Detector Download PDFInfo
- Publication number
- JP2014067269A JP2014067269A JP2012212836A JP2012212836A JP2014067269A JP 2014067269 A JP2014067269 A JP 2014067269A JP 2012212836 A JP2012212836 A JP 2012212836A JP 2012212836 A JP2012212836 A JP 2012212836A JP 2014067269 A JP2014067269 A JP 2014067269A
- Authority
- JP
- Japan
- Prior art keywords
- posture
- information
- vector
- pedestrian
- estimation
- 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
Links
Landscapes
- Traffic Control Systems (AREA)
- Image Analysis (AREA)
Abstract
Description
本発明は、撮像画像から識別対象を検出する検出装置に関する。 The present invention relates to a detection device that detects an identification target from a captured image.
従来、車両の周辺環境(車両進行方向に存在する移動体の有無や移動体の位置等)に応じた運転支援を実行するシステムについて、各種検討が行われている。
このようなシステムの一例として、車両に搭載された赤外線カメラ等により撮像された画像から歩行者を検出し、この歩行者が自車両にとって危険な対象であるか否か、具体的には、自車両への衝突の虞があるか否かを判断するもの(以下、従来装置という)が知られている(例えば、特許文献1参照)。
Conventionally, various studies have been made on a system that performs driving support in accordance with the surrounding environment of the vehicle (the presence or absence of a moving body in the vehicle traveling direction, the position of the moving body, and the like).
As an example of such a system, a pedestrian is detected from an image captured by an infrared camera or the like mounted on the vehicle, and whether or not the pedestrian is a dangerous object for the own vehicle, specifically, A device that determines whether or not there is a possibility of a collision with a vehicle (hereinafter referred to as a conventional device) is known (see, for example, Patent Document 1).
この従来装置では、車両に搭載された2台のカメラを用いて歩行者を検出し、この歩行者の全体(全身)のシルエットを水平方向にスライスした直線の中点をサンプル点として、このサンプル点を直線で回帰させる処理を画像上にて実行する。そして、この従来装置では、推定された回帰直線の鉛直方向からの傾きを歩行者の体幹の傾きとして、この体幹の傾きが車道の中心に向かって予め定めた値以上傾いていれば、この人物は車道への飛び出しの虞が高いと判断する。 In this conventional apparatus, a pedestrian is detected by using two cameras mounted on a vehicle, and the sample point is a midpoint of a straight line obtained by slicing the silhouette of the entire pedestrian (whole body) in the horizontal direction. The process of regressing the points with a straight line is executed on the image. And in this conventional device, if the inclination of the estimated regression line from the vertical direction is the inclination of the trunk of the pedestrian, the inclination of the trunk is inclined more than a predetermined value toward the center of the roadway, It is determined that this person has a high risk of jumping out onto the roadway.
ところで、体幹の傾きは、歩行者の一連の動きの中で最後に出現することが一般的である。つまり、例えば歩行者が車道に飛び出そうとしているならば、手が動く又は足が動く等、人体を構成する各パーツに何らかの動きが出現した後に、体幹の傾きが出現すると考えられる。 By the way, the inclination of the trunk generally appears last in a series of pedestrian movements. That is, for example, if a pedestrian is about to jump out of the road, it is considered that the inclination of the trunk appears after some movement appears in each part of the human body, such as a hand moving or a foot moving.
従来装置では、歩行者の姿勢(動き)を判断する際に、体幹の傾きが生じる迄の時間を要してしまう。そこで、車道に飛び出そうとする歩行者等、自車両にとって危険な姿勢(動き)を行う移動体を、より迅速に検出する技術が望まれる。 In the conventional apparatus, it takes time until the trunk tilts when determining the posture (movement) of the pedestrian. Therefore, a technique for more quickly detecting a moving body that performs a posture (movement) that is dangerous to the host vehicle, such as a pedestrian who is about to jump out on the roadway, is desired.
本発明は、上記問題点を解決するためになされたものであり、自車両にとって危険な対象となる移動体を、迅速に検出可能な検出装置を提供することを目的とする。 The present invention has been made to solve the above-described problems, and an object of the present invention is to provide a detection device that can quickly detect a moving body that is a dangerous target for the host vehicle.
上記目的を達成するためになされた本発明の検出装置は、車両に搭載され、パーツ位置検出手段が、撮像画像中の識別対象を部分的に表すパーツの位置(パーツ位置)を検出する。また、情報算出手段が、検出したパーツの位置関係の特徴を示すパーツ特徴情報を算出する。そして、姿勢推定手段が、パーツ特徴情報と識別対象のとる姿勢との関係を示す確率モデルを用いて、パーツ特徴情報から識別対象の姿勢を推定する。姿勢情報提供手段は、姿勢推定手段の推定結果に基づく情報、例えば自車両に対して危険な姿勢をとっている、等の様に識別対象についての情報を提供する。 The detection device of the present invention made to achieve the above object is mounted on a vehicle, and the part position detection means detects the position (part position) of a part that partially represents the identification target in the captured image. In addition, the information calculation means calculates part feature information indicating the feature of the positional relationship of the detected parts. Then, the posture estimation means estimates the posture of the identification target from the part feature information using a probability model indicating the relationship between the part feature information and the posture taken by the identification target. The attitude information providing means provides information on the identification target such as information based on the estimation result of the attitude estimating means, for example, taking a dangerous attitude with respect to the host vehicle.
なお、パーツとは識別対象を部分的に表す部位のことをいい、例えば識別対象が歩行者である場合、頭部、右肩部、左肩部、腰部、右足部、左足部等、身体を部分的に表す部位をいう。 Parts refer to parts that partially represent the identification target.For example, when the identification target is a pedestrian, the body part such as the head, right shoulder, left shoulder, waist, right foot, left foot, etc. This refers to the part that is expressed.
ここで、識別対象が何らかの動きを行うとき、識別対象を構成する各パーツに、まず何らかの動きが出現する。つまり、各パーツの位置関係に変化が生じ、しかもその変化が連動して生じると考えられる。 Here, when the identification target performs some movement, first, some movement appears in each part constituting the identification target. That is, it is considered that a change occurs in the positional relationship between the parts, and the change occurs in conjunction with the change.
本発明の検出装置では、この様なパーツを対象として各パーツ間の位置関係の特徴を抽出し、この特徴に基づき識別対象の姿勢(動き)を推定するため、結果として、より迅速に自車両にとって危険な対象となる識別対象を検出することができる。 In the detection apparatus of the present invention, the feature of the positional relationship between each part is extracted for such a part, and the posture (motion) of the identification target is estimated based on this feature. It is possible to detect an identification object that is a dangerous object for the user.
本発明の検出装置における情報算出手段は、請求項2に記載のように、パーツ位置に基づく複数の変量を成分とするベクトルを観測ベクトルとし、観測ベクトルの各成分について多変量解析における主成分分析を行うことにより得られる各主成分を列成分とした行列である予め用意された変換行列の転置行列を用いて、特定対象の観測ベクトルを座標変換してパーツ特徴情報を算出する。 As described in claim 2, the information calculation means in the detection apparatus of the present invention uses, as an observation vector, a vector having a plurality of variables based on part positions as components, and principal component analysis in multivariate analysis for each component of the observation vector Part feature information is calculated by performing coordinate transformation of the observation vector to be identified using a transposed matrix of a transformation matrix prepared in advance, which is a matrix with each principal component obtained as a column component.
より具体的には、例えば請求項3に記載のように、情報算出手段は、多変量解析における主成分分析を用いて、複数の識別対象に対して観測ベクトルの共分散行列を算出し、算出された該行列の固有値のうち大きい順に予め定められた数までの固有値について、対応する固有ベクトルを順に並べてなる行列を、変換行列として用いる。 More specifically, as described in claim 3, for example, the information calculation unit calculates a covariance matrix of observation vectors for a plurality of identification targets by using principal component analysis in multivariate analysis, and calculates For the eigenvalues up to a predetermined number in the descending order of the eigenvalues of the matrix, a matrix in which corresponding eigenvectors are arranged in order is used as a transformation matrix.
これによると、観測ベクトルがより低次元のパーツ特徴情報に変換されるため、このパーツ特徴情報を用いて歩行者の姿勢を推定する処理では、演算量が低減されると共に、ノイズに頑健な推定結果を得ることができる。結果として、本発明の検出装置は、迅速にかつ精度よく、自車両にとって危険な対象となる移動体を検出することができる。 According to this, since the observation vector is converted into lower-dimensional part feature information, the processing for estimating the pedestrian's posture using this part feature information reduces the amount of computation and makes the estimation robust against noise. The result can be obtained. As a result, the detection device of the present invention can quickly and accurately detect a moving object that is dangerous for the host vehicle.
以下、本発明が適用された実施形態について、図面を用いて説明する。
[実施形態]
<全体構成>
図1に示す移動体検出装置1は、自動車に搭載され、車載カメラ2及び車載表示装置3のそれぞれと通信可能に接続されている。車載カメラ2は、自動車の走行方向前方を撮像するカメラ(本実施形態では、単眼カメラ)であり、撮像画像を移動体検出装置1に出力する。
Embodiments to which the present invention is applied will be described below with reference to the drawings.
[Embodiment]
<Overall configuration>
A moving body detection device 1 shown in FIG. 1 is mounted on a vehicle and is connected to each of a vehicle-mounted camera 2 and a vehicle-mounted display device 3 so as to communicate with each other. The in-vehicle camera 2 is a camera (in this embodiment, a monocular camera) that captures an image in front of the vehicle in the traveling direction, and outputs a captured image to the moving object detection device 1.
移動体検出装置1は、車載カメラ2から入力した撮像画像において、識別対象として歩行者等の移動体を検出し、この移動体について姿勢及び動作等の情報を抽出する。この抽出した情報に基づいて、移動体検出装置1は、この移動体が自車両にとって危険な対象(危険対称)であると判断した場合、危険対象と判断された移動体を強調する表示を撮像画像に加えた画像(加工画像)を生成して車載表示装置3に出力する。 The moving body detection apparatus 1 detects a moving body such as a pedestrian as an identification target in the captured image input from the in-vehicle camera 2 and extracts information such as the posture and motion of the moving body. Based on the extracted information, when the moving body detection device 1 determines that the moving body is a dangerous target (danger symmetry) for the host vehicle, the moving body detection device 1 captures a display that highlights the moving body determined to be a dangerous target. An image (processed image) added to the image is generated and output to the in-vehicle display device 3.
車載表示装置3は、各種画像を表示し視覚的に情報を提示するためのディスプレイを少なくとも備え、移動体検出装置1から入力した加工画像を、自動車の運転者が視認できるように表示する。 The in-vehicle display device 3 includes at least a display for displaying various images and visually presenting information, and displays the processed image input from the mobile body detection device 1 so that the driver of the vehicle can visually recognize the image.
移動体検出装置1は、撮像画像から抽出した移動体について、この移動体の全体(全身)の位置と、この移動体を部分的に表す複数種類のパーツの位置とに基づく観測情報を出力する位置検出部10を備える。また、移動体検出装置1は、この観測情報に基づいて、移動体がどの方向に進もうとする姿勢をとっているのか(姿勢情報)、どのような動作をしようとしているのか(動作情報)、及び自車両からどれくらい離れて位置しているのか(距離情報)等の、移動体に関する情報を抽出する情報抽出部20を備える。さらに、移動体検出装置1は、この情報抽出部20にて抽出された移動体に関する情報を運転者に認識させることができる形のデータ(加工画像)として出力する情報提供部30を備える。 The moving body detection apparatus 1 outputs observation information based on the position of the entire moving body (whole body) and the positions of a plurality of types of parts that partially represent the moving body for the moving body extracted from the captured image. A position detection unit 10 is provided. In addition, based on this observation information, the moving body detection device 1 is in which direction the moving body is going to move (posture information), and what kind of operation is being performed (motion information). And an information extracting unit 20 that extracts information on the moving body such as how far away from the host vehicle (distance information) is provided. Furthermore, the mobile body detection device 1 includes an information providing unit 30 that outputs information on the mobile body extracted by the information extraction unit 20 as data (processed image) in a form that allows the driver to recognize the information.
ここで、移動体として、歩行者、自転車及びオートバイク(自動二輪車)等の車両に乗車する人物等、種々のものが考えられるが、本実施形態では、歩行者(自転車及びオートバイク(自動二輪車)等の車両に乗車していない人物)を検出対象としている。 Here, various types of moving objects such as a pedestrian, a person riding a vehicle such as a bicycle and a motorcycle (motorcycle), and the like can be considered. In this embodiment, a pedestrian (bicycle and motorcycle (motorcycle) is used. )) Is a detection target.
<位置検出部>
位置検出部10は、撮像画像の特徴量を算出する特徴量算出部41と、算出された特徴量を用いて撮像画像から歩行者の全体形状(シルエット)を識別する全体識別部42と、算出された特徴量を用いて撮像画像から歩行者を部分的に表す複数のパーツを識別するパーツ識別部43と、識別された歩行者のシルエット及び該歩行者のパーツの位置に基づく観測情報を出力する観測情報生成部50と、を備える。
<Position detector>
The position detection unit 10 includes a feature amount calculation unit 41 that calculates the feature amount of the captured image, an overall identification unit 42 that identifies the overall shape (silhouette) of the pedestrian from the captured image using the calculated feature amount, and a calculation. A part identifying unit 43 for identifying a plurality of parts partially representing a pedestrian from a captured image using the identified feature value, and outputting observation information based on the identified silhouette of the pedestrian and the position of the part of the pedestrian And an observation information generation unit 50.
特徴量算出部41は、撮像画像の画像データから輝度値で表現された画像データを生成し、この輝度値で表現された画像データに対して、画像を構成する画素ごとに勾配強度及び勾配方向を算出し、複数画素からなる領域(ブロック)のHOG(Histograms of Oriented Gradients)特徴量を算出する。HOG特徴量の算出方法は周知であるため、ここでは説明を省略する。 The feature amount calculation unit 41 generates image data expressed by luminance values from the image data of the captured image, and the gradient strength and gradient direction for each pixel constituting the image with respect to the image data expressed by the luminance values. And an HOG (Histograms of Oriented Gradients) feature amount of an area (block) composed of a plurality of pixels is calculated. Since the calculation method of the HOG feature value is well known, the description is omitted here.
全体識別部42では、予め学習等により定められた基準とすべき歩行者の全体形状(基準シルエット。例えば、歩行者の全身に相当する矩形領域(M0ブロック×N0ブロック。M0、N0は自然数)をいう。)について算出されたHOG特徴量が、全体基準特徴量として記憶されている。 In the overall identification unit 42, the overall shape of the pedestrian to be used as a reference predetermined by learning or the like (reference silhouette. For example, a rectangular area corresponding to the whole body of the pedestrian (M0 block × N0 block. M0 and N0 are natural numbers). The HOG feature value calculated for the above-mentioned is stored as an overall reference feature value.
この全体基準特徴量をテンプレートとして、全体識別部42は、撮像画像においてテンプレートの位置を順にずらしながら、各位置におけるテンプレートとの類似度(スコア)を算出する周知のテンプレートマッチングの処理を行う。各位置で算出された類似度は、その位置に歩行者(の基準シルエット)が存在する可能性の高さを示す。 Using this overall reference feature amount as a template, the overall identification unit 42 performs a well-known template matching process for calculating the similarity (score) with the template at each position while sequentially shifting the position of the template in the captured image. The similarity calculated at each position indicates the high possibility that a pedestrian (its reference silhouette) exists at that position.
このスコアが予め定めた閾値を超える場合、全体識別部42は、その位置に歩行者が存在すると判定し、この位置における基準シルエットの中心(M0×N0ブロックの矩形領域の中心)座標を算出し、算出した座標を全体座標(x_r、y_r)として出力する(図2参照)。 When this score exceeds a predetermined threshold, the overall identification unit 42 determines that a pedestrian exists at the position, and calculates the center of the reference silhouette (the center of the rectangular area of the M0 × N0 block) at this position. The calculated coordinates are output as overall coordinates (x_r, y_r) (see FIG. 2).
パーツ識別部43は、検出すべき移動体が歩行者である場合、撮像画像から頭部を識別する頭部識別部44、右側及び左側の肩部を識別する右肩識別部45及び左肩識別部46、腰部を識別する腰部識別部47、右側及び左側の足部を識別する右足識別部48及び左足識別部49を備える。 When the moving body to be detected is a pedestrian, the parts identifying unit 43 includes a head identifying unit 44 that identifies a head from a captured image, a right shoulder identifying unit 45 that identifies right and left shoulders, and a left shoulder identifying unit. 46, a waist identifying part 47 for identifying the waist, a right foot identifying part 48 for identifying the right and left foot parts, and a left foot identifying part 49.
これらのうち、例えば頭部識別部44には、予め学習等により定められた基準とすべき頭部の形状(例えば、頭部に相当する矩形領域(M1ブロック×N1ブロック。M1、N1は自然数。))について算出されたHOG特徴量が、頭部基準特徴量として記憶されている。 Among these, for example, the head identification unit 44 has a head shape to be used as a reference predetermined by learning or the like (for example, a rectangular area corresponding to the head (M1 block × N1 block. M1 and N1 are natural numbers). .)) Is calculated as a head reference feature value.
頭部識別部44は、頭部基準特徴量をテンプレートとして、撮像画像においてテンプレートの位置を順にずらしながら、全体識別部42にて実施した処理と同様のテンプレートマッチングの処理を行い、各位置における類似度(スコア)を算出する。このスコアが予め定めた閾値を超える場合、頭部識別部44は、その位置に頭部が存在すると判定し、基準頭部の領域の中心(M1×N1ブロックの矩形の中心)座標を算出し、算出した座標を頭部座標(x_p1、y_p1)として出力する。 The head identification unit 44 performs template matching processing similar to the processing performed by the overall identification unit 42 while sequentially shifting the position of the template in the captured image using the head reference feature value as a template. The degree (score) is calculated. When this score exceeds a predetermined threshold, the head identification unit 44 determines that the head is present at the position, and calculates the center of the reference head region (the center of the rectangle of the M1 × N1 block) coordinates. The calculated coordinates are output as head coordinates (x_p1, y_p1).
右肩識別部45、左肩識別部46、腰部識別部47、右足識別部48及び左足識別部49は、使用するテンプレートが検出するパーツによって異なる他は、頭部識別部44と同様の処理を行う。そして、各識別部は、右肩、左肩、胴、右足、左足が存在すると判定された撮像画像上の座標を、それぞれ、右肩座標(x_p2、y_p2)、左肩座標(x_p3、y_p3)、胴座標(x_p4、y_p4)、右足座標(x_p5、y_p5)、左足座標(x_p6、y_p6)として出力する。 The right shoulder identifying unit 45, the left shoulder identifying unit 46, the waist identifying unit 47, the right foot identifying unit 48, and the left foot identifying unit 49 perform the same processing as the head identifying unit 44, except that it differs depending on the parts detected by the template used. . Then, each identification unit determines the coordinates on the captured image determined to have the right shoulder, the left shoulder, the torso, the right foot, and the left foot as the right shoulder coordinates (x_p2, y_p2), the left shoulder coordinates (x_p3, y_p3), and the torso, respectively. Output as coordinates (x_p4, y_p4), right foot coordinates (x_p5, y_p5), and left foot coordinates (x_p6, y_p6).
観測情報生成部50は、頭部、両肩、腰部、両足の6つのパーツの相対的な位置関係から、これら6つのパーツを包含する矩形領域に、歩行者が存在する可能性の高さを示す総合類似度を算出する。算出された総合類似度が予め定めた閾値を越える場合、観測情報生成部50は、その矩形領域に歩行者が存在すると判定する。 From the relative positional relationship of the six parts of the head, both shoulders, waist, and both legs, the observation information generation unit 50 determines the high possibility that a pedestrian exists in a rectangular area that includes these six parts. The total similarity shown is calculated. When the calculated total similarity exceeds a predetermined threshold value, the observation information generation unit 50 determines that a pedestrian exists in the rectangular area.
そして、観測情報生成部50は、歩行者が存在すると判定した矩形領域に対応する全体座標及び各パーツの座標の値を順に並べてなるN(Nは自然数)次元の観測ベクトルvを、(1)式に示すように生成し、この観測ベクトルvを観測情報として出力する。但し、(1)式に示す各x及びyの値は、各座標の値を正規化した値とする。つまり、撮像画像から抽出した歩行者画像を同じ大きさで比較できる様に、全体識別部42及びパーツ識別部43から出力された各座標の値を正規化する。なお本実施形態では、観測ベクトルvは14次元のベクトルであるため(1(全体識別部の数)×2+6(パーツ識別部の数)×2=14)、以下の説明ではN=14とする。 Then, the observation information generation unit 50 sets an N (N is a natural number) -dimensional observation vector v in which the overall coordinates corresponding to the rectangular area determined to have a pedestrian and the coordinate values of each part are arranged in order (1). It produces | generates as shown in a type | formula, and this observation vector v is output as observation information. However, the values of x and y shown in the equation (1) are values obtained by normalizing the values of the coordinates. That is, the value of each coordinate output from the whole identification part 42 and the part identification part 43 is normalized so that the pedestrian image extracted from the captured image can be compared with the same size. In the present embodiment, since the observation vector v is a 14-dimensional vector (1 (number of overall identification units) × 2 + 6 (number of part identification units) × 2 = 14), N = 14 in the following description. .
この様に検出すべき歩行者について少なくとも一つのパーツ位置(座標)が検出されなかった場合、観測情報生成部50は、未検出のパーツ位置(座標)に対応する値として予め定められた代替値を当てはめ、観測ベクトルvを生成する。例えば、左足が抽出されなかった場合であれば、観測ベクトルvの左足の座標に相当する成分の値を代替値である「0」とする。 When at least one part position (coordinates) is not detected for a pedestrian to be detected in this way, the observation information generating unit 50 substitutes a predetermined alternative value as a value corresponding to an undetected part position (coordinates). To generate an observation vector v. For example, if the left foot is not extracted, the value of the component corresponding to the left foot coordinate of the observation vector v is set to “0” as an alternative value.
なお、本実施形態では、位置検出部10における一連の処理において、latentSVM及びDPMを用いた手法(詳細は、P.F.Felzenszwalb et al.,"Object Detection with Discriminatively Trained Part Based Models,"IEEE Trunsactions on Pattern Analysis and Machine Intelligence(2009). 参照)を利用し、頭部、右肩、左肩、胴、右足、左足の6つのパーツの相対的な位置関係を適切に保った状態のものを、人として検出する。 In the present embodiment, in the series of processing in the position detection unit 10, a technique using latent SVM and DPM (for details, see PFFelzenszwalb et al., “Object Detection with Discriminatively Trained Part Based Models,” IEEE Trunsactions on Pattern Analysis. and Machine Intelligence (2009).) to detect humans with the relative positional relationship of the six parts of the head, right shoulder, left shoulder, torso, right foot, and left foot properly maintained .
<情報抽出部>
情報抽出部20は、観測ベクトルvに基づいて、検出された歩行者の姿勢情報を抽出する姿勢情報抽出部60と、動作情報を抽出する動作情報抽出部70と、下端位置情報を出力する下端推定部80と、を備える。
<姿勢情報抽出部>
姿勢情報抽出部60は、歩行者の各パーツの位置関係の特徴を示す姿勢特徴ベクトルfを算出する姿勢特徴抽出部61と、姿勢特徴ベクトルfに基づいて移動体の姿勢を推定する姿勢推定部62とを備える。
<Information extraction unit>
The information extraction unit 20 is based on the observation vector v, the posture information extraction unit 60 that extracts the detected pedestrian posture information, the motion information extraction unit 70 that extracts the motion information, and the lower end that outputs the lower end position information. And an estimation unit 80.
<Attitude information extraction unit>
The posture information extraction unit 60 calculates a posture feature vector f indicating a positional relationship feature of each part of the pedestrian, and a posture estimation unit estimates the posture of the moving body based on the posture feature vector f. 62.
姿勢特徴抽出部61は、(2)式に示すように、N次元の観測ベクトルvに基づいてk(kは自然数。k<N。)次元の姿勢特徴ベクトルfを算出する。観測ベクトルvに比べて次元数が低減された姿勢特徴ベクトルfは、パーツ特徴情報として、姿勢推定部62及び下端推定部80に出力される。 The posture feature extraction unit 61 calculates a k (k is a natural number, k <N.)-Dimensional posture feature vector f based on the N-dimensional observation vector v, as shown in Equation (2). The posture feature vector f having a reduced number of dimensions compared to the observation vector v is output to the posture estimation unit 62 and the lower end estimation unit 80 as part feature information.
まず、事前に収集された多くの学習用の撮像画像から、M体(Mは自然数)の歩行者画像について、位置検出部10で算出したものと同様のN次元の観測ベクトルv(i)((1)式参照。i=1、2、・・M)を算出し、(3)式に示すように、これらM体の観測ベクトルv(i)の平均値を平均ベクトルvaとして算出する。学習用の撮像画像は、歩行者が撮像されている画像で、車載カメラ2により撮像された画像であることが、より望ましい。 First, an N-dimensional observation vector v (i) similar to that calculated by the position detection unit 10 for M pedestrian images (M is a natural number) from many learning captured images collected in advance. (1) see .i = 1, 2, calculates a · · M), as shown in (3), calculates these M of the observation vector v the mean value of (i) a mean vector v a . More preferably, the captured image for learning is an image in which a pedestrian is captured and is an image captured by the in-vehicle camera 2.
次に、(4)式に示すように、M体の各歩行者の観測ベクトルv(i)と平均ベクトルvaとの差分(偏差)を成分としてデータ行列Xを生成し、(5)式に示すように、データ行列Xを用いて分散共分散行列Sを算出する。 Next, (4) As shown in equation generates data matrix X difference M bodies each pedestrian observation vector v of (i) and the average vector v a a (deviation) as the component (5) As shown in FIG. 4, a variance covariance matrix S is calculated using the data matrix X.
ここで、i番目に大きい固有値λiに対応する固有ベクトルa(i)を第i主成分として、第1主成分a(1)〜第N主成分a(N)により構成されるN×Nの行列Aを算出する。そして、この行列Aの内、第1主成分a(1)〜第k主成分a(k)により構成されるN×kの部分基底行列Akを(7)式に示すように算出する(ここでは、N=14より、k<14)。 Here, an N × N matrix composed of the first principal component a (1) to the Nth principal component a (N) with the eigenvector a (i) corresponding to the i-th largest eigenvalue λi as the i-th principal component. A is calculated. Then, an N × k partial basis matrix Ak composed of the first principal component a (1) to the k-th principal component a (k) in the matrix A is calculated as shown in the equation (7) (here) Then, from N = 14, k <14).
なお、kの値は、固有値から算出される累積寄与率により決定する。ここでは、累積寄与率80%以上を満足するkの値を6(k=6)とし、第1主成分a(1)〜第6主成分a(6)により構成される14×6の部分基底行列Akを変換行列Akとした。 The value of k is determined by the cumulative contribution rate calculated from the eigenvalue. Here, the value of k that satisfies the cumulative contribution rate of 80% or more is 6 (k = 6), and a 14 × 6 portion composed of the first principal component a (1) to the sixth principal component a (6) The base matrix Ak is defined as the transformation matrix Ak.
つまり、学習用の多数(M体)の歩行者画像において、全身及び各パーツからなる歩行者の姿勢(形状)を、各パーツ位置を成分とするN次元のベクトル(観測ベクトルv)と見なし、主成分分析を行うことにより、歩行者の姿勢(形状)の特徴をもった比較的低次元(k次元)の部分空間が得られる。この部分空間を構成する基底ベクトルが、変換行列Akの主成分a(1)〜a(k)に対応している。 That is, in a large number of learning pedestrian images (M bodies), the posture (shape) of the pedestrian consisting of the whole body and each part is regarded as an N-dimensional vector (observation vector v) having each part position as a component, By performing principal component analysis, a relatively low-dimensional (k-dimensional) subspace having the characteristics of the posture (shape) of a pedestrian can be obtained. The basis vectors constituting this partial space correspond to the principal components a (1) to a (k) of the transformation matrix Ak.
姿勢特徴抽出部61が算出する姿勢特徴ベクトルfは、N次元の観測ベクトルv(移動体検出装置1の車載カメラ2で撮像した歩行者画像から生成した観測ベクトル)を、次元が低減されたk次元の部分空間に写像したものとなっている。そして、この姿勢特徴ベクトルfを用いて姿勢推定部62にて姿勢推定の処理を行うということは、この次元が低減された部分空間の中で姿勢推定の処理が行われることを意味する。 The posture feature vector f calculated by the posture feature extraction unit 61 is an N-dimensional observation vector v (an observation vector generated from a pedestrian image captured by the vehicle-mounted camera 2 of the mobile object detection device 1), and the dimension is reduced. It is mapped to a subspace of dimensions. Then, performing posture estimation processing by the posture estimation unit 62 using the posture feature vector f means that posture estimation processing is performed in a subspace in which this dimension is reduced.
姿勢推定部62は、k次元の姿勢特徴ベクトルfを有する歩行者画像について、歩行者の姿勢yが、姿勢ラベルrに相当する姿勢(推定姿勢)となっている姿勢確率p(y=r|f)(但し、r=1、2、・・・R。Rは自然数。)を、ロジスティック回帰手法を利用して、(8)式により算出する。 The posture estimation unit 62, for a pedestrian image having a k-dimensional posture feature vector f, a posture probability p (y = r |) in which the posture y of the pedestrian is a posture corresponding to the posture label r (estimated posture). f) (where r = 1, 2,... R, R is a natural number) is calculated by the equation (8) using a logistic regression method.
算出された姿勢確率p(y=r|f)は、姿勢情報として、情報提供部30に出力される。
ここで、(8)式における、各パラメータ(θr0、θr T、θq0、θq Tは、事前の実験によって、次の様に算出されたものである。
The calculated posture probability p (y = r | f) is output to the information providing unit 30 as posture information.
Here, the parameters (θ r0 , θ r T , θ q0 , θ q T ) in the equation (8) are calculated as follows by a prior experiment.
まず、事前に収集された多くの学習用の撮像画像(歩行者が撮像されている画像。車載カメラ2による撮像画像であることがより望ましい)から、多数(G体)の歩行者画像について、姿勢特徴抽出部61で算出したものと同様のk次元の姿勢特徴ベクトルf(j)((2)式参照。j=1、2、・・・G。Gは自然数。)を算出する。また、各歩行者画像における歩行者の姿勢を視認によって確認し、対応する姿勢ラベルrを1つ付与しておく。この様に、予め付与された姿勢ラベルrを「正解タグ」という。 First, from a large number of pre-collected captured images for learning (images in which pedestrians are captured, more preferably captured images by the in-vehicle camera 2), The same k-dimensional posture feature vector f (j) as that calculated by the posture feature extraction unit 61 (see equation (2). J = 1, 2,... G, G is a natural number) is calculated. Further, the posture of the pedestrian in each pedestrian image is confirmed by visual recognition, and one corresponding posture label r is given. In this way, the posture label r given in advance is referred to as a “correct tag”.
そして、学習用の撮像画像から算出された姿勢特徴ベクトルf(j)のそれぞれについて、k個の成分f1〜fkを説明変数とし、撮像画像に対応する「正解タグ」を目的変数として、(8)式が成り立つ様に、各パラメータ(スカラーθr0、θq0、k次元のベクトルθr T、θq T)を設定する。 Then, for each of the posture feature vectors f (j) calculated from the captured image for learning, k components f1 to fk are used as explanatory variables, and “correct answer tags” corresponding to the captured images are used as objective variables. Each parameter (scalar θ r0 , θ q0 , k-dimensional vectors θ r T , θ q T ) is set so that the following equation holds.
ここで、姿勢推定部62による歩行者の姿勢の推定結果を図3に示す。同図(a)は、事前に歩行者の撮像画像を複数準備し、各撮像画像から視認により確認された歩行者の姿勢(正解タグ)を白色で示した図である。同図(b)は、同じ撮像画像における歩行者の姿勢yを、姿勢推定部62によって推定した結果を示し図である。同図(b)では、姿勢確率p(y=r|f)の値は、グレースケールで表示されており、色が濃い(黒い)ほど0に近い値を示し、色が薄い(白い)ほど1に近い値を示す。 Here, the estimation result of the posture of the pedestrian by the posture estimation unit 62 is shown in FIG. The figure (a) is the figure which prepared the captured image of a pedestrian in advance, and showed the pedestrian's attitude | position (correct tag) confirmed by visual recognition from each captured image in white. FIG. 6B is a diagram illustrating a result of estimating the pedestrian's posture y in the same captured image by the posture estimation unit 62. In FIG. 6B, the value of the posture probability p (y = r | f) is displayed in gray scale, and the darker the color (black), the closer to 0, and the lighter the color (white). A value close to 1 is shown.
例えば、最初の数枚(フレーム)についてみると、正解タグとして姿勢ラベル1が付与されており(図3(a)参照)、歩行者が「前」向きの姿勢をとっていたことが確認される。これに対して、同じ最初の数フレームを姿勢推定部62により推定した結果では(同図(b)参照)、正解タグである「前」向の姿勢をとっている姿勢確率p(y=1|f)が、他の姿勢をとっている確率に比べて大きな値として算出されている。つまり、姿勢推定部62による推定結果は高い確率で正解タグと一致することが確認されている。 For example, in the first few frames (frames), posture label 1 is given as a correct tag (see FIG. 3 (a)), and it is confirmed that the pedestrian has taken a “front” posture. The On the other hand, in the result of estimating the same first few frames by the posture estimation unit 62 (see FIG. 5B), the posture probability p (y = 1) taking the “forward” posture that is the correct tag. | F) is calculated as a value larger than the probability of taking another posture. That is, it has been confirmed that the estimation result by the posture estimation unit 62 matches the correct tag with a high probability.
この様に、姿勢情報抽出部60は、多次元(N次元)のデータ(観測ベクトルv)から得られた、より低次元(k次元)の特徴的な指標(姿勢特徴ベクトルf)を用いて、姿勢の推定を行うため、姿勢推定の処理において、演算量が低減されると共に、ノイズに頑健な処理結果が得られる。 In this way, the posture information extraction unit 60 uses a lower-dimensional (k-dimensional) characteristic index (posture feature vector f) obtained from multi-dimensional (N-dimensional) data (observation vector v). Since the posture is estimated, the amount of calculation is reduced in the posture estimation processing, and a processing result robust to noise is obtained.
<動作情報抽出部>
動作情報抽出部70は、動作特徴抽出部71と動作推定部72とを備える。動作情報抽出部70は、事前の実験において観測ベクトルvに代えてL(Lは自然数。L>N)次元の動作観測ベクトルvmを観測して主成分分析を行う点、及び、主成分分析を行った結果を利用し、入力画像における動作観測ベクトルvmに基づいて、歩行者が動作ラベルrm(rmは自然数)に対応する動作をとっている確率を算出する点が異なる他は、姿勢情報抽出部60と同様の処理を行う。
<Operation information extraction unit>
The motion information extraction unit 70 includes a motion feature extraction unit 71 and a motion estimation unit 72. Operation information extracting unit 70, instead of the observation vector v in the pre-experiment L (L is a natural number .L> N) by observing the behavior observation vector v m dimensions that performs principal component analysis, and principal component analysis utilizing the results of, based on the operation observation vector v m in the input image, pedestrian operation label r m (r m is a natural number) in addition to different points to calculate the probability of taking an operation corresponding to the The same processing as that of the posture information extraction unit 60 is performed.
つまり、動作特徴抽出部71は姿勢特徴抽出部61に相当し、動作推定部72は姿勢推定部62に相当する。動作特徴抽出部71は、動作観測ベクトルvmを生成し、この動作観測ベクトルvmに基づいて姿勢特徴ベクトルfに相当する動作特徴ベクトルを算出し、動作特徴情報として動作推定部72に出力する。動作推定部72は、動作特徴情報(動作特徴ベクトル)に基づき、撮像画像の歩行者が各動作ラベルrmに対応する動作をとっている確率を算出し、算出結果を動作情報として情報提供部30に出力する。 That is, the motion feature extraction unit 71 corresponds to the posture feature extraction unit 61, and the motion estimation unit 72 corresponds to the posture estimation unit 62. Operation feature extracting unit 71 generates an operation observation vector v m, and calculates the operation characteristic vector corresponding to the orientation feature vector f on the basis of the behavior observation vector v m, and outputs the movement-estimating unit 72 as the operation characteristic information . Motion estimation unit 72, based on the operation characteristic information (operation feature vectors), to calculate the probability that the pedestrian of the captured image is taking an operation corresponding to the operation labeled r m, the information providing unit the calculation result as operation information Output to 30.
動作観測ベクトルvmとは、検出すべき歩行者について、時系列に連続した撮像画像における観測ベクトルv((1)式参照)を、予め定められたフレーム数H(Hは2以上の自然数。ここでは、H=10とした。)だけ、単純に時系列に並べて構成されるL次元(ここでは、L=N×H=14×10=140)のベクトルをいう。 The motion observation vector v m is a predetermined number of frames H (H is a natural number equal to or greater than 2). Here, only H = 10) is an L-dimensional vector (here, L = N × H = 14 × 10 = 140) configured in a time series.
また、動作ラベルrmは3クラスからなり、動作ラベル1(rm=1)は止まる動作を表し、動作ラベル2(rm=2)は歩く動作を表し、動作ラベル3(rm=3)は走る動作を表わすものとする。 The action label r m is composed of three classes, the action label 1 (r m = 1) represents a stop action, the action label 2 (r m = 2) represents a walking action, and the action label 3 (r m = 3). ) Represents a running action.
ここで、動作情報抽出部70による歩行者の動作の推定結果を図4に示す。同図(a)は、事前に同一の歩行者について時系列に撮像した画像を複数準備し、撮像画像から視認によって確認された該歩行者の動作(正解タグ)を白色で示した図である。同図(b)は、同じ撮像画像における歩行者の動作を、動作推定部72によって推定した結果を示した図である。同図(b)では、各動作ラベルrmに対応する動作をとっていると推定される確率は、色が濃い(黒い)ほど0に近く、色が薄い(白い)ほど1に近い値として示されている。 Here, the estimation result of the motion of the pedestrian by the motion information extraction unit 70 is shown in FIG. The figure (a) is the figure which prepared several images imaged in time series about the same pedestrian beforehand, and showed the operation | movement (correct tag) of this pedestrian confirmed by visual recognition from the captured image in white. . FIG. 6B is a diagram showing a result of estimating the motion of the pedestrian in the same captured image by the motion estimation unit 72. In FIG. (B), the probability that is estimated to have taken an operation corresponding to the operation labeled r m is close to 0 as darker (black), as close to 1 as light color (white) It is shown.
例えば、最初の数フレームについてみると、正解タグとして動作ラベル2が付与されており(図4(a)参照)、歩行者が「歩く」動作をしていたことが正解である。これに対して、同じ最初の数フレームを動作推定部72により推定した結果では(同図(b)参照)、正解タグである「歩く」動作を行っている確率が、他の動作を行っている確率に比べて大きな値として算出されている。つまり、
動作推定部72による推定結果は高い確率で正解タグと一致することが確認されている。
For example, in the first few frames, the action label 2 is given as a correct tag (see FIG. 4A), and the correct answer is that the pedestrian was “walking”. On the other hand, according to the result of estimating the same first few frames by the motion estimation unit 72 (see (b) in the figure), the probability of performing the “walking” motion that is the correct tag is determined by performing other motions. It is calculated as a larger value than the probability of being. That means
It has been confirmed that the estimation result by the motion estimation unit 72 matches the correct tag with a high probability.
この様に、動作情報抽出部70は、姿勢情報抽出部60と同様の効果を得ることができる。つまり、多次元(N×H次元)のデータ(動作観測ベクトルvm)から得られた、より低次元の特徴的な指標(動作特徴ベクトル)を用いて動作の推定が行われるため、動作の推定処理において演算量が低減されると共に、ノイズに頑健な処理結果が得られる。 In this way, the motion information extraction unit 70 can obtain the same effects as the posture information extraction unit 60. That is, since motion estimation is performed using a lower-dimensional characteristic index (motion feature vector) obtained from multi-dimensional (N × H dimension) data (motion observation vector v m ), In the estimation process, the amount of calculation is reduced, and a processing result robust to noise is obtained.
<下端位置推定部>
下端推定部80は、姿勢特徴抽出部61から出力されたk次元の姿勢特徴ベクトルfに、予め用意された上述の変換行列Akを作用させて、(9)式に示すように、N次元の推定ベクトルvdを算出する。
<Lower position estimation part>
The lower end estimation unit 80 applies the above-described transformation matrix Ak prepared in advance to the k-dimensional posture feature vector f output from the posture feature extraction unit 61, and as shown in the equation (9), an N-dimensional An estimated vector v d is calculated.
ここで、上述のように、遮蔽等の影響により未検出となったパーツの情報として予め定められた代替値(ここでは0)を当てはめて観測ベクトルvを生成しておくことにより、未検出であったパーツ位置の情報を推定ベクトルvdの成分として得ることができる。 Here, as described above, an observation vector v is generated by applying a predetermined alternative value (0 in this case) as information on parts that have not been detected due to the influence of shielding or the like. Information on the part position that was present can be obtained as a component of the estimated vector v d .
この様な推定ベクトルvdを用いて、下端推定部80は、撮像画像における歩行者の下端位置座標を推定する。下端位置座標は、右足座標及び左足座標よりさらに下方に位置する。 Using such an estimated vector v d , the lower end estimation unit 80 estimates the lower end position coordinates of the pedestrian in the captured image. The lower end position coordinates are located further below the right foot coordinates and the left foot coordinates.
具体的には、下端推定部80は、推定ベクトルvdの成分を説明変数とし、結果である下端位置座標U=(x、y)Tを、線形回帰モデル(下端推定モデル)に基づき(10)式により推定する。
推定された下端位置座標U=(x、y)Tは、下端位置情報として情報提供部30に出力される。 The estimated lower end position coordinate U = (x, y) T is output to the information providing unit 30 as lower end position information.
ここで、(10)式における回帰係数W、及びバイアスベクトルbの算出方法について説明する。これらは、事前の実験により、収集された多くの撮像画像から最小二乗法により決定されたものである。回帰係数W(i)は回帰係数を示す2×Nの行列であり、バイアスベクトルbは2×1の行列である。 Here, a method of calculating the regression coefficient W and the bias vector b in the equation (10) will be described. These are determined by a least-squares method from a large number of captured images collected by a prior experiment. The regression coefficient W (i) is a 2 × N matrix indicating the regression coefficient, and the bias vector b is a 2 × 1 matrix.
回帰係数W、及びバイアスベクトルbを図5に示す様に決定した場合、推定された下端位置座標Uの推定誤差は、図6に示す様に、ほぼ5画素以内に収まることが確認された。
この様に、下端推定部80では、下端位置座標Uを、各パーツの位置関係に基づいて推定するため、下端位置座標Uが精度よく算出される。
When the regression coefficient W and the bias vector b were determined as shown in FIG. 5, it was confirmed that the estimated error of the estimated lower end position coordinate U was within approximately 5 pixels as shown in FIG.
In this way, since the lower end estimation unit 80 estimates the lower end position coordinate U based on the positional relationship of each part, the lower end position coordinate U is accurately calculated.
<情報提供部>
情報提供部30は、姿勢情報、動作情報、及び下端位置情報に基づき、自車両にとって危険な対象であると判断した歩行者(以下、危険対象という)を強調する表示を撮像画像に加えた加工画像を生成する。
<Information provision department>
The information providing unit 30 adds a display that emphasizes a pedestrian (hereinafter referred to as a dangerous target) that is determined to be a dangerous target for the host vehicle to the captured image based on the posture information, the motion information, and the lower end position information. Generate an image.
本実施形態では、情報提供部30は、図7に示す様に、下端位置情報に基づき、各歩行者91〜94の下端位置を示す下端線101〜104を元の撮像画像に合成した加工画像を生成する。 In the present embodiment, as shown in FIG. 7, the information providing unit 30 combines the lower end lines 101 to 104 indicating the lower end positions of the pedestrians 91 to 94 with the original captured image based on the lower end position information. Is generated.
また、情報提供部30は、歩行者が前後方向、左右方向のうちいずれの方向に進もうとする姿勢を取っているかを示すマーク(姿勢表示マーク)を元の画像に合成した加工画像を生成する。 In addition, the information providing unit 30 generates a processed image in which a mark (posture display mark) indicating whether the pedestrian is going to move in the front-rear direction or the left-right direction is combined with the original image. To do.
具体的には、情報提供部30は、姿勢情報に基づき、識別対象が「前」向及び「後」向の姿勢をとっている姿勢確率p(y=r|f)(但しr=1又は2)の合計値を「縦」向の姿勢をとっている確率とし、「左」向及び「右」向の姿勢をとっている姿勢確率p(y=r|f)(但しr=3又は4)の合計値を「横」向の姿勢をとっている確率として、「縦」「横」の姿勢をとる確率の比を縦横比とした楕円の姿勢表示マークを生成し、該姿勢表示マークを元の画像に合成した加工画像を生成する。つまり、図7に示す様に、「縦」向きの歩行者91については縦長の楕円111が姿勢表示マークとして表示され、「横」向きの歩行者92、93及び94に対しては、横長の楕円112〜114が姿勢表示マークとして表示される。 Specifically, based on the posture information, the information providing unit 30 has a posture probability p (y = r | f) (provided that r = 1 or 2) is the probability that the posture in the “vertical” direction is taken, and the posture probability p (y = r | f) in which the posture in the “left” direction and the “right” direction is taken (where r = 3 or 4) The total value of 4) is used as a probability that the posture in the “horizontal” direction is taken, and an elliptical posture display mark is generated with a ratio of the probability of taking the “vertical” and “horizontal” postures, and the posture display mark A processed image is generated by combining the image with the original image. That is, as shown in FIG. 7, a vertically long ellipse 111 is displayed as a posture display mark for the “vertical” pedestrian 91, and a horizontally long pedestrian 92, 93, and 94 is Ellipses 112 to 114 are displayed as posture display marks.
さらにまた、情報提供部30は、下端位置情報と撮像画像中の無限遠点とを利用して、自車両から歩行者までの距離を算出し、算出した距離を元の撮像画像に追加表示する。撮像画像中の無限遠点を利用した対象物までの距離の算出方法は周知であるため(例えば、特許第4063170号参照)、ここでは説明を省略する。 Furthermore, the information providing unit 30 calculates the distance from the host vehicle to the pedestrian using the lower end position information and the infinity point in the captured image, and additionally displays the calculated distance on the original captured image. . Since the calculation method of the distance to the target object using the infinity point in the captured image is well known (for example, refer to Japanese Patent No. 4063170), description thereof is omitted here.
また、情報提供部30は、「止まる」以外の動作、つまり「歩く」又は「走る」動作を行っている歩行者であって、かつ、「縦」方向より自車両にとって危険が生じる虞のある「横」方向に進もうとしている歩行者を危険対象として、動作推定情報に基づき危険対象と判定された歩行者を囲む枠を元の撮像画像に合成した加工画像を生成する。つまり、図7に示す様に、「歩く」動作を行っており、かつその向きが「横」向きである歩行者92〜94について、これらの歩行者92〜94を強調する様に、枠122〜124が追加表示される。 In addition, the information providing unit 30 is a pedestrian performing an operation other than “stop”, that is, a “walking” or “running” operation, and there is a risk of danger to the own vehicle from the “vertical” direction. A processed image is generated by synthesizing a frame surrounding a pedestrian determined to be a dangerous object based on the motion estimation information with the original captured image, with a pedestrian trying to move in the “lateral” direction as a dangerous object. In other words, as shown in FIG. 7, for the pedestrians 92 to 94 that are performing the “walking” operation and the direction is “lateral”, the frame 122 is emphasized so as to emphasize these pedestrians 92 to 94. -124 are additionally displayed.
この様に情報提供部30は、抽出された姿勢情報、動作情報、及び下端位置情報に基づく歩行者に関する種々の情報を、加工画像という見やすい形で提供することにより、運転者に危険対象についての注意を促す。 In this way, the information providing unit 30 provides various information about the pedestrian based on the extracted posture information, motion information, and lower end position information in a form that is easy to see as a processed image, so that the driver can be informed about the danger target. Call attention.
[効果]
以上説明したように、本実施形態の移動体検出装置1は、姿勢情報抽出部60によって、検出すべき歩行者について各パーツの位置情報を成分とする観測ベクトルvを生成し、この観測ベクトルvを各パーツの位置関係の特徴を示す姿勢特徴ベクトルfに変換し、この姿勢特徴ベクトルfに基づき歩行者の姿勢を推定する。これによると、体幹の傾きに基づき姿勢を推定する場合と比べて、自車両に対して危険な姿勢をとる歩行者を迅速に検出することができる。
[effect]
As described above, in the moving body detection apparatus 1 of the present embodiment, the posture information extraction unit 60 generates the observation vector v having the position information of each part as a component for the pedestrian to be detected, and this observation vector v Is converted into a posture feature vector f indicating the positional relationship feature of each part, and the posture of the pedestrian is estimated based on this posture feature vector f. According to this, it is possible to quickly detect a pedestrian who takes a dangerous posture with respect to the host vehicle, as compared with the case where the posture is estimated based on the inclination of the trunk.
また、姿勢特徴ベクトルfは、撮像画像における観測ベクトルvに比べて次元が低減されるため、姿勢特徴ベクトルfを用いた姿勢の推定処理では、演算量が低減される共に、ノイズに頑健な推定結果を得ることができる。結果として、自車両にとって危険な対象となる歩行者を、迅速にかつ精度よく検出することができる。 Further, since the pose feature vector f has a reduced dimension compared to the observed vector v in the captured image, the pose estimation process using the pose feature vector f reduces the amount of computation and is robust against noise. The result can be obtained. As a result, it is possible to quickly and accurately detect a pedestrian that is a dangerous target for the host vehicle.
なお、移動体検出装置1は、姿勢を推定するにあたり、ロジスティック回帰手法を用いたモデルを確率モデルとして採用した。これにより、歩行者がどの姿勢をとっているかという結果(離散変数)を、姿勢特徴ベクトルfから予測することができる。 Note that the moving object detection apparatus 1 employs a model using a logistic regression method as a probability model in estimating the posture. As a result, the result (discrete variable) indicating which posture the pedestrian is taking can be predicted from the posture feature vector f.
また、移動体検出装置1は、動作情報抽出部70によって、上述の観測ベクトルvに基づいて歩行者の姿勢を推定する考え方を時系列的に拡張し、動作観測ベクトルvmに基づいて歩行者の動作を推定する。これによると、姿勢に加えて動作についての情報を利用できるため、自車両にとって危険な対象となる歩行者を、より精度よく検出することができる。 In addition, the moving body detection device 1 uses the motion information extraction unit 70 to extend the concept of estimating the pedestrian's posture based on the observation vector v described above in a time series, and based on the motion observation vector v m. Estimate the behavior of According to this, since information about an operation in addition to a posture can be used, a pedestrian that is a dangerous target for the host vehicle can be detected with higher accuracy.
さらにまた、移動体検出装置1は、下端推定部80によって算出される推定ベクトルvdに基づいて、下端位置座標Uを推定する。これにより、例えば、パーツ識別部43にて観測された右足座標及び左足座標だけに基づいて下端位置座標を推定する場合に比べて、信頼性の高い推定結果が得られるだけでなく、
右足座標又は左足座標が未検出であった場合でも、精度よく下端位置座標Uを推定することができる。結果として、例えばこの下端位置座標Uの情報を用いて、単眼カメラで撮像した撮像画像においても、歩行者までの距離を精度よく推定することができるようになる。
Furthermore, the moving body detection apparatus 1 estimates the lower end position coordinate U based on the estimated vector v d calculated by the lower end estimation unit 80. Thereby, for example, as compared with the case where the lower end position coordinates are estimated based only on the right foot coordinates and the left foot coordinates observed by the parts identifying unit 43, not only a highly reliable estimation result is obtained,
Even when the right foot coordinate or the left foot coordinate is not detected, the lower end position coordinate U can be estimated with high accuracy. As a result, for example, the distance to the pedestrian can be accurately estimated even in the captured image captured by the monocular camera using the information on the lower end position coordinate U.
[特許請求の範囲との対応]
本実施形態の移動体検出装置1が特許請求の範囲の「検出装置」に相当し、車載カメラ2が特許請求の範囲の「撮像手段」に相当し、位置検出部10が特許請求の範囲の「パーツ位置検出手段」に相当する。また、姿勢特徴抽出部61が特許請求の範囲の「情報算出手段」に相当し、姿勢推定部62が特許請求の範囲の「姿勢推定手段」に相当する。さらにまた、動作特徴抽出部71が特許請求の範囲の「動作情報算出手段」に相当し、動作推定部72が特許請求の範囲の「動作推定手段」に相当し、下端推定部80が特許請求の範囲の「下端推定手段」に相当する。また、情報提供部30が特許請求の範囲の「姿勢情報提供手段」、「動作情報提供手段」、「距離情報提供手段」に相当する。
[Correspondence with Claims]
The mobile body detection device 1 of the present embodiment corresponds to the “detection device” in the claims, the in-vehicle camera 2 corresponds to the “imaging means” in the claims, and the position detection unit 10 corresponds to the claims. This corresponds to “part position detection means”. The posture feature extraction unit 61 corresponds to “information calculation unit” in the claims, and the posture estimation unit 62 corresponds to “posture estimation unit” in the claims. Further, the motion feature extraction unit 71 corresponds to “motion information calculation means” in the claims, the motion estimation unit 72 corresponds to “motion estimation means” in the claims, and the lower end estimation unit 80 claims. Corresponds to “lower end estimation means”. The information providing unit 30 corresponds to “attitude information providing means”, “motion information providing means”, and “distance information providing means” in the claims.
[他の実施形態]
以上、本発明の実施形態について説明したが、本発明は、上記実施形態に限定されることなく、種々の形態をとり得ることは言うまでもない。
[Other Embodiments]
As mentioned above, although embodiment of this invention was described, it cannot be overemphasized that this invention can take a various form, without being limited to the said embodiment.
(イ)上記実施形態では、位置検出部10にて各パーツの座標を求める手法として、latentSVM及びDPMを用いた手法を利用したが、各パーツの座標を算出する手法は、これに限るものではなく、他の学習法やモデルによる手法を適用してもよい。 (B) In the above embodiment, the technique using the latent SVM and DPM is used as the technique for obtaining the coordinates of each part in the position detection unit 10, but the technique for calculating the coordinates of each part is not limited to this. Alternatively, other learning methods and models may be applied.
(ロ)上記実施形態では、情報抽出部20にて姿勢特徴ベクトルf及び動作特徴ベクトルを算出するにあたり、共分散行列の対角化によって得られた主成分を用いた主成分分析を行ったが、共分散行列の代わりに自己相関行列を用いて主成分を得るようにしてもよい。 (B) In the above embodiment, the principal component analysis using the principal component obtained by diagonalizing the covariance matrix is performed when the information extraction unit 20 calculates the posture feature vector f and the motion feature vector. The principal component may be obtained using an autocorrelation matrix instead of the covariance matrix.
(ハ)上記実施形態では、姿勢推定部62にて、姿勢ラベルrのクラス数K=4としたが、クラス数Kは任意に定めてよい。つまり、斜め右、斜め左等、さらに細かく分けてクラス数Kを増加させても良い。また、前後を合わせて縦方向、左右を合わせて横方向として、クラス数をK=2の様に低減してもよい。動作ラベルrmのクラス数についても同様である。 (C) In the above-described embodiment, the posture estimation unit 62 sets the number of classes K of the posture label r = 4, but the number of classes K may be arbitrarily determined. In other words, the number of classes K may be increased in more detail, such as diagonally right and diagonally left. Alternatively, the number of classes may be reduced as K = 2, with the front and rear being the vertical direction and the left and right being the horizontal direction. The same applies to the number of classes operation label r m.
(二)上記実施形態では識別対象を歩行者としていたが、識別対象は、自転車や、オートバイ、三輪車、車椅子等の車両に乗車する人物であってもよい。但し、自動車のように、乗員が車両の内部に乗車するものは除外される。 (2) In the above embodiment, the identification target is a pedestrian, but the identification target may be a person who rides a vehicle such as a bicycle, a motorcycle, a tricycle, or a wheelchair. However, a vehicle in which an occupant gets inside the vehicle, such as an automobile, is excluded.
(ホ)上記実施形態では、歩行者を表すパーツとして、頭、両肩、腰、両足を設定したが、パーツとして選択する部位はこれらに限られるわけではない。パーツは、移動体に何らかの動きが出現するときに、相対的な位置変化が現れやすい部位を含むものであればよい。 (E) In the above embodiment, the head, both shoulders, waist, and both feet are set as the parts representing the pedestrian, but the parts to be selected as the parts are not limited to these. The part only needs to include a part where a relative position change is likely to appear when some movement appears in the moving body.
(ヘ)上記実施系形態で説明した移動体検出装置1の構成要素は、ハードウェアで実現してもよく、ソフトウェアで実現してもよく、ハードウェアとソフトウェアとの組み合わせで実現してもよい。また、これらの構成要素は機能概念的なものであり、その一部又は全部を、機能的又は物理的に分散又は統合してもよい。 (F) The components of the moving object detection apparatus 1 described in the above embodiment may be realized by hardware, software, or a combination of hardware and software. . These components are functionally conceptual, and some or all of them may be functionally or physically distributed or integrated.
1・・・移動体検出装置 2・・・車載カメラ 3・・・車載表示装置 10・・・位置検出部 20・・・情報抽出部 30・・・情報提供部 41・・・特徴量算出部 42・・・全体識別器 43・・・パーツ識別器 50・・・観測情報生成部 60・・・姿勢情報抽出部 70・・・動作情報抽出部 80・・・下端推定部 DESCRIPTION OF SYMBOLS 1 ... Moving body detection apparatus 2 ... Car-mounted camera 3 ... Car-mounted display apparatus 10 ... Position detection part 20 ... Information extraction part 30 ... Information provision part 41 ... Feature-value calculation part 42 ... Whole classifier 43 ... Parts classifier 50 ... Observation information generation unit 60 ... Attitude information extraction unit 70 ... Motion information extraction unit 80 ... Lower end estimation unit
Claims (8)
識別対象の画像を撮像する撮像手段(2)と、
前記識別対象のうち一つを特定対象とし、前記識別対象を部分的に表す複数のパーツの位置をパーツ位置として、前記撮像手段が撮像した画像から、前記特定対象のパーツ位置を検出するパーツ位置検出手段(10)と、
前記パーツの位置関係の特徴を示す情報をパーツ特徴情報として、前記特定対象について、前記パーツ特徴情報を算出する情報算出手段(61)と、
前記識別対象の前記パーツ特徴情報と該識別対象の姿勢との関係を示す確率モデルを用いて、前記撮像手段が取得した画像中の前記特定対象について、前記パーツ特徴情報から前記姿勢を推定する姿勢推定手段(62)と、
前記姿勢推定手段の推定結果に基づく情報を提供する姿勢情報提供手段(30)と、
を備えることを特徴とする検出装置。 A detection device mounted on a vehicle,
Imaging means (2) for imaging an image to be identified;
A part position that detects a part position of the specific target from an image captured by the imaging unit, where one of the identification targets is a specific target and positions of a plurality of parts that partially represent the identification target are part positions Detection means (10);
Information indicating the feature of the positional relationship of the parts as part feature information, information calculating means (61) for calculating the part feature information for the specific target;
A posture for estimating the posture from the part feature information for the specific target in the image acquired by the imaging unit using a probability model indicating a relationship between the part feature information of the identification target and the posture of the identification target. Estimating means (62);
Attitude information providing means (30) for providing information based on the estimation result of the attitude estimation means;
A detection apparatus comprising:
前記情報算出手段は、前記特定対象について、予め用意された前記変換行列の転置行列を用いて、前記観測ベクトルを座標変換して前記パーツ特徴情報を算出すること、を特徴とする請求項1に記載の検出装置。 A matrix having a principal component as a column component obtained by performing a principal component analysis in multivariate analysis on each component of the observation vector as an observation vector and a vector having a plurality of variables based on the part position as a transformation matrix As
The information calculation means calculates the part feature information by performing coordinate conversion of the observation vector using a transpose matrix of the conversion matrix prepared in advance for the specific target. The detection device described.
前記多変量解析における主成分分析を用いて、複数の前記識別対象に対して前記観測ベクトルの共分散行列を算出し、算出された該行列の固有値のうち大きい順に予め定められた数までの固有値について、対応する固有ベクトルを順に並べてなる行列を、前記変換行列として用いること、を特徴とする請求項2に記載の検出装置。 The information calculating means includes
The principal component analysis in the multivariate analysis is used to calculate a covariance matrix of the observation vectors for a plurality of the identification targets, and eigenvalues up to a predetermined number in descending order among the calculated eigenvalues of the matrix The detection device according to claim 2, wherein a matrix in which corresponding eigenvectors are arranged in order is used as the transformation matrix.
前記識別対象の取りうる姿勢として予め定めた少なくとも一つの姿勢を推定姿勢とし、前記識別対象が有する前記パーツ特徴情報の成分を説明変数とし、前記識別対象がそれぞれの前記推定姿勢を取りうる確率を目的変数とした、ロジスティック回帰手法を用いたモデルであること、を特徴とする請求項1から3のいずれか一項に記載の検出装置。 The probability model is
Estimated posture is at least one predetermined posture as the posture that can be taken by the identification target, the component of the part feature information that the identification target has as an explanatory variable, and the probability that the identification target can take the respective estimated posture The detection apparatus according to claim 1, wherein the detection apparatus is a model using a logistic regression method as an objective variable.
前記情報算出手段は、未検出の前記パーツの位置に代替値を当てはめて前記観測ベクトルを生成すること、を特徴とする請求項2から4のいずれか一項に記載の検出装置。 In the part position detecting means, for the specific object where the position of at least one part is not detected,
5. The detection apparatus according to claim 2, wherein the information calculation unit generates the observation vector by applying an alternative value to a position of the undetected part. 6.
前記動作特徴情報に基づいて前記特定対象の動作を推定する動作推定手段(72)と、
前記動作推定手段の推定結果に基づく情報を提供する動作情報提供手段(30)と、
を備えることを特徴とする請求項1から5のいずれか一項に記載の検出装置。 Motion information calculation means for calculating, for the specific target, motion characteristic information indicating the characteristics of the positional relationship of the parts in time series obtained by observing the part positions detected by the part position detection means in time series. 71)
Motion estimation means (72) for estimating the motion of the specific target based on the motion feature information;
Motion information providing means (30) for providing information based on the estimation result of the motion estimation means;
The detection apparatus according to claim 1, further comprising:
前記パーツ位置と前記識別対象の下端位置との関係を示す下端推定モデルを用いて、前記特定対象の前記推定ベクトルに基づいて該特定対象の下端位置を推定する下端推定手段(80)と、
前記下端推定手段により推定された下端位置に基づいて、前記特定対象と自車両との距離を算出し、算出された前記距離に基づく情報を提供する距離情報提供手段(30)と、
を備えることを特徴とする請求項1から6のいずれか一項に記載の検出装置。 As an estimated vector, a vector obtained by operating the transformation matrix on the part feature information,
Lower end estimation means (80) for estimating a lower end position of the specific target based on the estimation vector of the specific target using a lower end estimation model indicating a relationship between the part position and the lower end position of the identification target;
Based on the lower end position estimated by the lower end estimating means, a distance information providing means (30) for calculating a distance between the specific object and the host vehicle and providing information based on the calculated distance;
The detection apparatus according to any one of claims 1 to 6, further comprising:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2012212836A JP5983243B2 (en) | 2012-09-26 | 2012-09-26 | Detection device |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2012212836A JP5983243B2 (en) | 2012-09-26 | 2012-09-26 | Detection device |
Publications (2)
Publication Number | Publication Date |
---|---|
JP2014067269A true JP2014067269A (en) | 2014-04-17 |
JP5983243B2 JP5983243B2 (en) | 2016-08-31 |
Family
ID=50743593
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2012212836A Active JP5983243B2 (en) | 2012-09-26 | 2012-09-26 | Detection device |
Country Status (1)
Country | Link |
---|---|
JP (1) | JP5983243B2 (en) |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2017056385A1 (en) * | 2015-09-29 | 2017-04-06 | ソニー株式会社 | Information processing device, information processing method, and program |
WO2017056382A1 (en) * | 2015-09-29 | 2017-04-06 | ソニー株式会社 | Information processing device, information processing method, and program |
JP2019522187A (en) * | 2016-06-02 | 2019-08-08 | ノキア テクノロジーズ オーユー | Apparatus and related methods |
KR20200077982A (en) * | 2018-12-21 | 2020-07-01 | 네이버랩스 주식회사 | Method and system for predicting future motion of object in still image |
KR20200077942A (en) * | 2018-12-21 | 2020-07-01 | 네이버랩스 주식회사 | Method and system for tracking object in video through prediction of future motion of object |
EP3754608A4 (en) * | 2018-08-01 | 2021-08-18 | Tencent Technology (Shenzhen) Company Limited | Target tracking method, computer device, and storage medium |
WO2022246767A1 (en) * | 2021-05-27 | 2022-12-01 | 华为技术有限公司 | Method and device for determining steering intention of target vehicle |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH0546583A (en) * | 1991-08-15 | 1993-02-26 | Nippon Telegr & Teleph Corp <Ntt> | Confirmation device for moving body action |
JP2010165003A (en) * | 2009-01-13 | 2010-07-29 | Toyota Central R&D Labs Inc | Pedestrian recognition support device and pedestrian recognition support method |
JP2011150381A (en) * | 2010-01-19 | 2011-08-04 | Nippon Telegr & Teleph Corp <Ntt> | Expression discriminator creation device, expression discriminator creation method, expression recognition device, expression recognition method, and program therefor |
WO2012077286A1 (en) * | 2010-12-09 | 2012-06-14 | パナソニック株式会社 | Object detection device and object detection method |
JP2012128638A (en) * | 2010-12-15 | 2012-07-05 | Canon Inc | Image processing device, alignment method and program |
-
2012
- 2012-09-26 JP JP2012212836A patent/JP5983243B2/en active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH0546583A (en) * | 1991-08-15 | 1993-02-26 | Nippon Telegr & Teleph Corp <Ntt> | Confirmation device for moving body action |
JP2010165003A (en) * | 2009-01-13 | 2010-07-29 | Toyota Central R&D Labs Inc | Pedestrian recognition support device and pedestrian recognition support method |
JP2011150381A (en) * | 2010-01-19 | 2011-08-04 | Nippon Telegr & Teleph Corp <Ntt> | Expression discriminator creation device, expression discriminator creation method, expression recognition device, expression recognition method, and program therefor |
WO2012077286A1 (en) * | 2010-12-09 | 2012-06-14 | パナソニック株式会社 | Object detection device and object detection method |
JP2012128638A (en) * | 2010-12-15 | 2012-07-05 | Canon Inc | Image processing device, alignment method and program |
Non-Patent Citations (1)
Title |
---|
角野公亮 外2名: "3次元モデルによるパラメトリック固有空間法を用いた多関節物体の姿勢推定の高度化に関する検討", 電子情報通信学会技術研究報告 IE2009−131−IE2009−178 画像工学, vol. 第109巻 第415号, JPN6015052712, 8 February 2010 (2010-02-08), pages 1 - 4, ISSN: 0003226840 * |
Cited By (16)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US10949656B2 (en) | 2015-09-29 | 2021-03-16 | Sony Corporation | Information processing apparatus and information processing method |
WO2017056382A1 (en) * | 2015-09-29 | 2017-04-06 | ソニー株式会社 | Information processing device, information processing method, and program |
CN108028021A (en) * | 2015-09-29 | 2018-05-11 | 索尼公司 | Information processing equipment, information processing method and program |
JPWO2017056382A1 (en) * | 2015-09-29 | 2018-08-09 | ソニー株式会社 | Information processing apparatus, information processing method, and program |
US20180268700A1 (en) * | 2015-09-29 | 2018-09-20 | Sony Corporation | Information processing apparatus, information processing method, and program |
US10692379B2 (en) | 2015-09-29 | 2020-06-23 | Sony Corporation | Image processing device and image processing method |
US11915522B2 (en) | 2015-09-29 | 2024-02-27 | Sony Corporation | Information processing apparatus and information processing method |
WO2017056385A1 (en) * | 2015-09-29 | 2017-04-06 | ソニー株式会社 | Information processing device, information processing method, and program |
JP2019522187A (en) * | 2016-06-02 | 2019-08-08 | ノキア テクノロジーズ オーユー | Apparatus and related methods |
EP3754608A4 (en) * | 2018-08-01 | 2021-08-18 | Tencent Technology (Shenzhen) Company Limited | Target tracking method, computer device, and storage medium |
US11961242B2 (en) | 2018-08-01 | 2024-04-16 | Tencent Technology (Shenzhen) Company Limited | Target tracking method, computer device, and storage medium |
KR102143031B1 (en) * | 2018-12-21 | 2020-08-10 | 네이버랩스 주식회사 | Method and system for predicting future motion of object in still image |
KR102143034B1 (en) * | 2018-12-21 | 2020-08-10 | 네이버랩스 주식회사 | Method and system for tracking object in video through prediction of future motion of object |
KR20200077942A (en) * | 2018-12-21 | 2020-07-01 | 네이버랩스 주식회사 | Method and system for tracking object in video through prediction of future motion of object |
KR20200077982A (en) * | 2018-12-21 | 2020-07-01 | 네이버랩스 주식회사 | Method and system for predicting future motion of object in still image |
WO2022246767A1 (en) * | 2021-05-27 | 2022-12-01 | 华为技术有限公司 | Method and device for determining steering intention of target vehicle |
Also Published As
Publication number | Publication date |
---|---|
JP5983243B2 (en) | 2016-08-31 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP5983243B2 (en) | Detection device | |
EP3539054B1 (en) | Neural network image processing apparatus | |
JP5647155B2 (en) | Body feature detection and human pose estimation using inner distance shape relation | |
JP5848341B2 (en) | Tracking by monocular 3D pose estimation and detection | |
Rezaei et al. | Look at the driver, look at the road: No distraction! no accident! | |
JP4677046B2 (en) | Fast human pose estimation using appearance and motion via multidimensional boost regression | |
US9213892B2 (en) | Real-time bicyclist detection with synthetic training data | |
US9367735B2 (en) | Object identification device | |
JP2016062610A (en) | Feature model creation method and feature model creation device | |
JP6570786B2 (en) | Motion learning device, skill discrimination device, and skill discrimination system | |
Wang et al. | Point linking network for object detection | |
JP2020135747A (en) | Action analysis device and action analysis method | |
Browning et al. | Cortical dynamics of navigation and steering in natural scenes: Motion-based object segmentation, heading, and obstacle avoidance | |
US20220180555A1 (en) | Estimation device, estimation method, and storage medium | |
US20220164577A1 (en) | Object detection method, object detection apparatus, and non-transitory computer-readable storage medium storing computer program | |
WO2017170761A1 (en) | Three-dimensional video data generation device, three-dimensional video data generation program, and three-dimensional video data generation method | |
US20230410356A1 (en) | Detection apparatus, detection method, and non-transitory storage medium | |
López-Fernández et al. | independent gait recognition through morphological descriptions of 3D human reconstructions | |
Drory et al. | A learning-based markerless approach for full-body kinematics estimation in-natura from a single image | |
Jarraya et al. | Biometric-based security system for smart riding clubs | |
CN115131752A (en) | Learning method, learning device, and program recording medium | |
JP7412514B1 (en) | Cabin monitoring method and cabin monitoring system that implements the above cabin monitoring method | |
KR102459293B1 (en) | Method and apparatus for generating mesh model of human or quadrupeds | |
Huang et al. | Real-Time Human Body Motion Capturing System | |
Hema et al. | Gait Based Human Age Classification using Random Forest Classifier |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20150313 |
|
A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20151224 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20160105 |
|
A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20160303 |
|
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: 20160705 |
|
A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20160718 |
|
R151 | Written notification of patent or utility model registration |
Ref document number: 5983243 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R151 |
|
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 |