[go: up one dir, main page]
More Web Proxy on the site http://driver.im/

JP7377494B2 - Attribute identification device and identification method - Google Patents

Attribute identification device and identification method Download PDF

Info

Publication number
JP7377494B2
JP7377494B2 JP2020083417A JP2020083417A JP7377494B2 JP 7377494 B2 JP7377494 B2 JP 7377494B2 JP 2020083417 A JP2020083417 A JP 2020083417A JP 2020083417 A JP2020083417 A JP 2020083417A JP 7377494 B2 JP7377494 B2 JP 7377494B2
Authority
JP
Japan
Prior art keywords
radar
data
cluster
value
unit
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
JP2020083417A
Other languages
Japanese (ja)
Other versions
JP2021179323A (en
Inventor
将広 星
邦哲 矢野
幸伸 時枝
清介 京地
憂希 千々和
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
THE UNIVERSITY OF KITAKYUSHU
Original Assignee
THE UNIVERSITY OF KITAKYUSHU
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 THE UNIVERSITY OF KITAKYUSHU filed Critical THE UNIVERSITY OF KITAKYUSHU
Priority to JP2020083417A priority Critical patent/JP7377494B2/en
Publication of JP2021179323A publication Critical patent/JP2021179323A/en
Application granted granted Critical
Publication of JP7377494B2 publication Critical patent/JP7377494B2/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Traffic Control Systems (AREA)
  • Radar Systems Or Details Thereof (AREA)

Description

この発明は、属性の識別装置および識別方法に関し、特に、電波を用いて物標としての歩行者の属性を識別する技術に関する。 The present invention relates to an attribute identification device and an identification method, and particularly to a technique for identifying the attributes of a pedestrian as a target object using radio waves.

歩行者の属性を識別する技術として、例えば、25台のカメラからなる多視点同期歩容撮影システムを構築して、広い年齢層の被験者の歩容を撮影することで歩容データベースを作成し、性別と年齢とのそれぞれに関して識別実験を行った結果に基づいて性別および年齢の識別に適した分類として子供、成人男性、成人女性、および高齢者という4つのクラスを設定し、前記4つのクラスへの識別を行う、技術が知られている(非特許文献1参照)。 As a technology for identifying the attributes of pedestrians, for example, we built a multi-view synchronized gait photography system consisting of 25 cameras and created a gait database by photographing the gaits of subjects in a wide range of age groups. Based on the results of a discrimination experiment regarding gender and age, four classes, child, adult male, adult female, and elderly, are set as classifications suitable for gender and age discrimination. There is a known technique for identifying (see Non-Patent Document 1).

万波秀年ほか「歩容における性別・年齢の分類と特徴解析」,電子情報通信学会論文誌 D,社団法人電子情報通信学会,2009年,Vol.J92-D,No.8,pp.1373-1382Hidetoshi Manba et al. “Classification and characteristic analysis of gender and age in gait”, IEICE Transactions D, Institute of Electronics, Information and Communication Engineers, 2009, Vol. J92-D, No. 8, pp. 1373-1382

しかしながら、個人の属性を識別するために多数のカメラが必要とされるような大掛かりなシステムは、実用的とは言い難い、という問題がある。また、実環境での使用を考慮した場合、画像では特に人同士の重なりなどの背景変化に弱く、部位ごとの周波数特徴を効果的に抽出することは困難である、という問題がある。 However, there is a problem in that a large-scale system that requires a large number of cameras to identify personal attributes is hardly practical. Furthermore, when considering use in a real environment, there is a problem in that images are particularly sensitive to background changes such as overlapping people, making it difficult to effectively extract frequency features for each body part.

そこで本発明は、簡便な機器構成によって歩行者の属性の識別を高精度に行うことが可能な属性の識別装置および識別方法を提供することを目的とする。 SUMMARY OF THE INVENTION Therefore, an object of the present invention is to provide an attribute identification device and an identification method that can identify attributes of pedestrians with high accuracy using a simple device configuration.

上記課題を解決するために、請求項1に記載の発明は、検知対象エリアへと電波を送信するとともに前記検知対象エリアからの電波を受信してレーダスキャンを行うレーダ装置から出力されるレーダデータを周波数解析して周波数スペクトルを出力する周波数解析部と、前記周波数スペクトルに基づいて前記検知対象エリア内を移動する物標に関するレーダ信号点それぞれの相対速度を計算する速度計算部と、前記周波数スペクトルのピークを検出して前記レーダ信号点を抽出する信号点抽出部と、抽出された前記レーダ信号点のクラスタを生成するクラスタ生成部と、前記クラスタの代表点を特定する代表点特定部と、前記代表点を含む関心領域を設定する関心領域設定部と、前記関心領域に含まれる前記レーダ信号点それぞれの前記相対速度の値別に信号強度情報を計算する信号強度特定部と、所定の時間長さぶんの前記レーダスキャンごとの前記相対速度の値と前記信号強度情報の値との組み合わせデータの集合を時系列に配列した上で、前記組み合わせデータの時系列の配列を時間軸方向に1つずつシフトさせながら、シフト前の、時系列で配列された前記組み合わせデータの集合との対比における相関値を計算する相関計算部と、前記相関値を用いて前記物標の属性を判別する識別部と、を有する、ことを特徴とする属性の識別装置である。 In order to solve the above problem, the invention according to claim 1 provides radar data output from a radar device that transmits radio waves to a detection target area, receives radio waves from the detection target area, and performs a radar scan. a frequency analysis unit that performs frequency analysis to output a frequency spectrum; a speed calculation unit that calculates the relative speed of each radar signal point regarding a target moving within the detection target area based on the frequency spectrum; a signal point extracting unit that detects a peak of and extracts the radar signal point; a cluster generating unit that generates a cluster of the extracted radar signal points; and a representative point specifying unit that specifies a representative point of the cluster; a region of interest setting unit that sets a region of interest including the representative point; a signal strength specifying unit that calculates signal strength information for each of the relative velocity values of each of the radar signal points included in the region of interest; and a predetermined time length. After arranging a set of combination data of the relative velocity value and the signal strength information value for each of the radar scans in a time series, the time series arrangement of the combination data is arranged once in the time axis direction. a correlation calculation unit that calculates a correlation value in comparison with the set of the combination data arranged in time series before the shift while shifting the data, and an identification unit that determines the attribute of the target using the correlation value. It is an attribute identification device characterized by having the following.

請求項2に記載の発明は、請求項1に記載の属性の識別装置において、前記識別部が、機械学習のアルゴリズムを用いて前記物標の属性を判別する、ことを特徴とする。 The invention according to claim 2 is the attribute identification device according to claim 1, wherein the identification unit determines the attribute of the target object using a machine learning algorithm.

また、請求項3に記載の発明は、検知対象エリアへと電波を送信するとともに前記検知対象エリアからの電波を受信してレーダスキャンを行うレーダ装置から出力されるレーダデータを周波数解析して周波数スペクトルを出力する処理と、前記周波数スペクトルに基づいて前記検知対象エリア内を移動する物標に関するレーダ信号点それぞれの相対速度を計算する処理と、前記周波数スペクトルのピークを検出して前記レーダ信号点を抽出する処理と、抽出された前記レーダ信号点のクラスタを生成する処理と、前記クラスタの代表点を特定する処理と、前記代表点を含む関心領域を設定する処理と、前記関心領域に含まれる前記レーダ信号点それぞれの前記相対速度の値別に信号強度情報を計算する処理と、所定の時間長さぶんの前記レーダスキャンごとの前記相対速度の値と前記信号強度情報の値との組み合わせデータの集合を時系列に配列した上で、前記組み合わせデータの時系列の配列を時間軸方向に1つずつシフトさせながら、シフト前の、時系列で配列された前記組み合わせデータの集合との対比における相関値を計算する処理と、前記相関値を用いて前記物標の属性を判別する処理と、を有する、ことを特徴とする属性の識別方法である。 Further, the invention according to claim 3 provides a frequency analysis method for frequency-analyzing radar data output from a radar device that transmits radio waves to a detection target area, receives radio waves from the detection target area, and performs a radar scan. a process of outputting a spectrum; a process of calculating a relative velocity of each radar signal point regarding a target moving within the detection target area based on the frequency spectrum; and a process of detecting a peak of the frequency spectrum and detecting the radar signal point. a process of generating a cluster of the extracted radar signal points; a process of identifying a representative point of the cluster; a process of setting a region of interest including the representative point; processing for calculating signal strength information for each of the relative velocity values of each of the radar signal points, and combination data of the relative velocity value and the signal strength information value for each radar scan for a predetermined length of time; After arranging the sets of the combined data in chronological order, and then shifting the time-series arrangement of the combined data one by one in the time axis direction, in comparison with the set of the combined data arranged in the chronological order before the shift. This is an attribute identification method characterized by comprising a process of calculating a correlation value, and a process of determining an attribute of the target object using the correlation value.

請求項4に記載の発明は、請求項3に記載の属性の識別方法において、機械学習のアルゴリズムを用いて前記物標の属性を判別する、ことを特徴とする。 The invention according to claim 4 is characterized in that, in the attribute identification method according to claim 3, the attribute of the target object is determined using a machine learning algorithm.

請求項1に記載の発明や請求項3に記載の発明によれば、相対速度の値と信号強度情報の値との組み合わせデータの集合の時系列の配列を時間軸方向に1つずつシフトさせながら、シフト前の、時系列で配列された前記組み合わせデータの集合との対比として計算される相関値を用いて物標の属性を判別するようにしているので、物標の属性に応じた特有の性質を的確に把握して物標の属性の識別を高精度に行うことが可能となる。 According to the invention described in claim 1 or the invention described in claim 3, the time series arrangement of the set of combination data of relative velocity values and signal strength information values is shifted one by one in the time axis direction. However, since the attribute of the target is determined using the correlation value calculated as a comparison with the set of combination data arranged in time series before the shift, the characteristic It becomes possible to accurately grasp the properties of the target object and identify the attributes of the target object with high precision.

また、レーダ(電波)であれば、画像処理のような複雑な前処理を施す必要がなく、物標それぞれの速度成分を定量的に抽出することができるので、周波数領域の特徴を簡易かつ高精度に求めることが可能となる。 In addition, with radar (radio waves), there is no need for complex preprocessing such as image processing, and the velocity components of each target object can be quantitatively extracted. It becomes possible to obtain accuracy.

請求項2に記載の発明や請求項4に記載の発明によれば、機械学習のアルゴリズムを用いるようにしているので、相対速度の値と信号強度情報の値との組み合わせデータの集合の時系列の配列に関する、時間軸方向における相関値に基づいて物標の属性の識別を適切に行うことが可能となる。 According to the invention described in claim 2 or the invention described in claim 4, since a machine learning algorithm is used, the time series of a set of combination data of relative velocity values and signal strength information values is It becomes possible to appropriately identify the attribute of the target object based on the correlation value in the time axis direction regarding the arrangement of the target object.

この発明の実施の形態に係る属性の識別装置の概略構成を示す機能ブロック図である。1 is a functional block diagram showing a schematic configuration of an attribute identification device according to an embodiment of the present invention. FIG. 図1の属性の識別装置における処理手順であるとともにこの発明の実施の形態に係る属性の識別方法の処理手順を示すフローチャートである。2 is a flowchart showing a processing procedure in the attribute identification device of FIG. 1 and a processing procedure of an attribute identification method according to an embodiment of the present invention. 速度・強度の時系列データの一例を示すグラフである。It is a graph showing an example of time series data of speed/intensity. 相関計算部における相関値の計算の仕方を説明する図である。FIG. 3 is a diagram illustrating how a correlation value is calculated in a correlation calculation section. クラスタについての相関値の集合のデータの一例を示すグラフである。7 is a graph showing an example of data of a set of correlation values for clusters. クラスタの相関値の集合における相関値のピーク間の差分を説明する図である。FIG. 3 is a diagram illustrating a difference between peaks of correlation values in a set of correlation values of clusters. クラスタの相関値の集合における相関値の極大値および極小値を説明する図である。FIG. 3 is a diagram illustrating local maximum values and minimum values of correlation values in a set of correlation values of a cluster. クラスタの相関値の集合の、高速フーリエ変換処理後のL1ノルムの例を示す図である。FIG. 7 is a diagram showing an example of the L 1 norm of a set of cluster correlation values after fast Fourier transform processing. クラスタの相関値の集合に関する、データのシフト回ごとの相関値の変動について、男性の場合と女性の場合との比較の例を示す図である。FIG. 7 is a diagram illustrating an example of a comparison between a case of a man and a case of a woman regarding fluctuations in a correlation value for each shift of data regarding a set of correlation values of a cluster. クラスタの相関値の集合に関する、高速フーリエ変換処理後のL1ノルムについて、男性の場合と女性の場合との比較の例を示す図である。FIG. 6 is a diagram illustrating an example of a comparison between a man's case and a woman's case regarding the L 1 norm after fast Fourier transform processing regarding a set of correlation values of clusters. 相関値に纏わる3つの指標についての直交3軸座標系に、クラスタ毎にプロットした例を示す図である。FIG. 7 is a diagram showing an example of three indicators related to correlation values plotted for each cluster on an orthogonal three-axis coordinate system. 図11に示すデータについてサポートベクターマシンを用いて判別処理を行って超平面を計算した結果を示す図である。FIG. 12 is a diagram showing the results of calculating a hyperplane by performing discrimination processing on the data shown in FIG. 11 using a support vector machine. 物標を手先に限定した場合の、クラスタの相関値の集合に関する、データのシフト回ごとの相関値の変動の例を示す図である。(A)は成人の場合の図である。(B)は子供の場合の図である。FIG. 7 is a diagram illustrating an example of variation in correlation values for each shift of data regarding a set of correlation values of clusters when the target is limited to the hand; FIG. (A) is a diagram for an adult. (B) is a diagram for a child. 物標をももに限定した場合の、クラスタの相関値の集合に関する、データのシフト回ごとの相関値の変動の例を示す図である。(A)は成人の場合の図である。(B)は子供の場合の図である。FIG. 7 is a diagram illustrating an example of variation in correlation values for each shift of data regarding a set of correlation values of clusters when the target object is limited to thighs. (A) is a diagram for an adult. (B) is a diagram for a child. 物標を手先に限定した場合の、クラスタの相関値の集合に関する、高速フーリエ変換処理後のL1ノルムの例を示す図である。(A)は成人の場合の図である。(B)は子供の場合の図である。FIG. 7 is a diagram illustrating an example of the L 1 norm after fast Fourier transform processing regarding a set of cluster correlation values when the target is limited to the hand. (A) is a diagram for an adult. (B) is a diagram for a child. 物標をももに限定した場合の、クラスタの相関値の集合に関する、高速フーリエ変換処理後のL1ノルムの例を示す図である。(A)は成人の場合の図である。(B)は子供の場合の図である。FIG. 7 is a diagram showing an example of the L 1 norm after fast Fourier transform processing regarding a set of cluster correlation values when the target object is limited to the thigh. (A) is a diagram for an adult. (B) is a diagram for a child.

以下、この発明を図示の実施の形態に基づいて説明する。 The present invention will be described below based on the illustrated embodiments.

図1は、この発明の実施の形態に係る属性の識別装置1の概略構成を示す機能ブロック図である。図2は、属性の識別装置1における処理手順であるとともにこの発明の実施の形態に係る属性の識別方法の処理手順を示すフローチャートである。 FIG. 1 is a functional block diagram showing a schematic configuration of an attribute identification device 1 according to an embodiment of the present invention. FIG. 2 is a flowchart showing the processing procedure in the attribute identification device 1 and the processing procedure of the attribute identification method according to the embodiment of the present invention.

この属性の識別装置1は、検知対象エリア内の移動物体を物標として検知して前記物標の属性を識別する装置であり、主として、信号処理部2と、物標位置検出部3と、物標計算部4と、特徴量計算部5と、識別部6と、を備える。属性の識別装置1は、具体的には、検知対象エリア内を移動している人の属性を識別する。 This attribute identification device 1 is a device that detects a moving object in a detection target area as a target and identifies the attribute of the target, and mainly includes a signal processing section 2, a target position detection section 3, It includes a target calculation section 4, a feature amount calculation section 5, and an identification section 6. Specifically, the attribute identification device 1 identifies the attributes of a person moving within the detection target area.

この実施の形態に係る属性の識別装置1は、検知対象エリアへと電波を送信するとともに前記検知対象エリアからの電波を受信してレーダスキャンを行うレーダ装置から出力されるレーダデータを周波数解析して周波数スペクトルを出力する周波数解析部21と、周波数スペクトルに基づいて検知対象エリア内を移動する物標に関するレーダ信号点それぞれの相対速度を計算する速度計算部23と、周波数スペクトルのピークを検出してレーダ信号点を抽出する信号点抽出部31と、抽出されたレーダ信号点のクラスタを生成するクラスタ生成部32と、クラスタの代表点を特定する代表点特定部34と、代表点を含む関心領域を設定する関心領域設定部41と、関心領域に含まれる前記レーダ信号点それぞれの前記相対速度の値別に信号強度情報を計算する信号強度特定部42と、所定の時間長さぶんのレーダスキャンごとの相対速度の値と信号強度情報の値との組み合わせデータの集合を時系列に配列した上で、前記組み合わせデータの時系列の配列を時間軸方向に1つずつシフトさせながら、シフト前の、時系列で配列された前記組み合わせデータの集合との対比における相関値を計算する相関計算部51と、相関値を用いて物標の属性を判別する識別部6と、を有する、ようにしている(図1参照)。 The attribute identification device 1 according to this embodiment performs frequency analysis on radar data output from a radar device that transmits radio waves to a detection target area, receives radio waves from the detection target area, and performs a radar scan. a frequency analysis section 21 that outputs a frequency spectrum based on the frequency spectrum, a speed calculation section 23 that calculates the relative speed of each radar signal point regarding a target moving within the detection target area based on the frequency spectrum, and a speed calculation section 23 that detects the peak of the frequency spectrum. a signal point extracting unit 31 that extracts radar signal points by using the radar signal points; a cluster generating unit 32 that generates clusters of the extracted radar signal points; a representative point specifying unit 34 that specifies representative points of clusters; A region of interest setting section 41 that sets a region, a signal strength specifying section 42 that calculates signal strength information for each of the relative velocity values of each of the radar signal points included in the region of interest, and a radar scan for a predetermined length of time. After arranging a set of combination data of relative velocity values and signal strength information values for each time in chronological order, while shifting the time series arrangement of the combination data one by one in the time axis direction, , comprising a correlation calculation unit 51 that calculates a correlation value in comparison with the set of the combination data arranged in time series, and an identification unit 6 that determines the attribute of the target using the correlation value. (See Figure 1).

また、この実施の形態に係る属性の識別方法は、検知対象エリアへと電波を送信するとともに前記検知対象エリアからの電波を受信してレーダスキャンを行うレーダ装置から出力されるレーダデータを周波数解析して周波数スペクトルを出力する処理(ステップS1)と、周波数スペクトルに基づいて検知対象エリア内を移動する物標に関するレーダ信号点それぞれの相対速度を計算する処理(ステップS3)と、周波数スペクトルのピークを検出してレーダ信号点を抽出する処理(ステップS4)と、抽出されたレーダ信号点のクラスタを生成する処理(ステップS5)と、クラスタの代表点を特定する処理(ステップS6)と、代表点を含む関心領域を設定する処理(ステップS7)と、関心領域に含まれる前記レーダ信号点それぞれの前記相対速度の値別に信号強度情報を計算する処理(ステップS8)と、所定の時間長さぶんのレーダスキャンごとの相対速度の値と信号強度情報の値との組み合わせデータの集合を時系列に配列した上で、前記組み合わせデータの時系列の配列を時間軸方向に1つずつシフトさせながら、シフト前の、時系列で配列された前記組み合わせデータの集合との対比における相関値を計算する処理(ステップS9)と、相関値を用いて物標の属性を判別する処理(ステップS10、S11)と、を有する、ようにしている(図2参照)。 Further, the attribute identification method according to this embodiment includes frequency analysis of radar data output from a radar device that transmits radio waves to a detection target area, receives radio waves from the detection target area, and performs a radar scan. a process of outputting a frequency spectrum (step S1), a process of calculating the relative velocity of each radar signal point regarding a target moving within the detection target area based on the frequency spectrum (step S3), and calculating the peak of the frequency spectrum. A process of detecting a radar signal point and extracting a radar signal point (step S4), a process of generating a cluster of the extracted radar signal points (step S5), a process of identifying a representative point of a cluster (step S6), A process of setting a region of interest including points (step S7), a process of calculating signal strength information for each of the relative velocity values of each of the radar signal points included in the region of interest (step S8), and a predetermined time length. After arranging a set of combination data of relative velocity values and signal strength information values for each radar scan in a time series, the time series arrangement of the combination data is shifted one by one in the time axis direction. , a process of calculating a correlation value in comparison with the set of the combination data arranged in time series before shifting (step S9), and a process of determining the attribute of the target using the correlation value (steps S10, S11). ) and (see Figure 2).

この発明では、例えば、レーダ方式としてFMCW(Frequency Modulated-Continuous Wave の略;周波数変調連続波)方式のレーダスキャンを行うレーダ装置によって取得されて出力されるレーダデータが用いられ得る。FMCW方式では、周波数変調した連続波を送信するとともに物標の表面で反射された反射波を受信し、送信信号と受信信号との間での周波数差に基づいて物標の位置や相対速度を検出する。なお、FMCW方式では、周波数が三角波状または鋸波状に変化する電波(送信信号)が送信される。 In the present invention, for example, radar data obtained and output by a radar device that performs radar scanning using an FMCW (Frequency Modulated Continuous Wave) method may be used as a radar method. The FMCW method transmits a frequency-modulated continuous wave and receives the reflected wave reflected from the surface of the target object, and calculates the position and relative speed of the target object based on the frequency difference between the transmitted signal and the received signal. To detect. Note that in the FMCW method, radio waves (transmission signals) whose frequency changes in the shape of a triangular wave or a sawtooth wave are transmitted.

この発明で用いられ得るレーダデータを出力するレーダ装置としては、例えば、送信部と受信部とを備えて電波を送信するとともに受信する機能を備え、FMCW方式のレーダスキャンを行ってレーダデータを出力する装置が挙げられる。 A radar device that outputs radar data that can be used in the present invention is, for example, equipped with a transmitting section and a receiving section and having a function of transmitting and receiving radio waves, and outputting radar data by performing FMCW radar scanning. Examples include devices that do this.

レーダ装置の送信部としては、例えば、所定電圧を生成して出力する電圧発生器、電圧発生器から出力される前記所定電圧に応じた周波数を有する電波(送信信号)を生成して出力する電圧制御発振器、および、電圧制御発振器から出力される電波(送信信号)を検知対象エリアへと送信波として送信する(別言すると、出射する、放射する)送信アンテナなどを含む仕組みが挙げられる。 The transmitter of the radar device includes, for example, a voltage generator that generates and outputs a predetermined voltage, and a voltage that generates and outputs a radio wave (transmission signal) having a frequency according to the predetermined voltage output from the voltage generator. Examples of mechanisms include a controlled oscillator and a transmitting antenna that transmits (in other words, emits or radiates) radio waves (transmission signals) output from the voltage-controlled oscillator to the detection target area as transmission waves.

送信部は、例えば79GHz帯、60GHz帯、或いは76GHz帯の周波数を有する電波(尚、「ミリ波」とも呼ばれる)を生成して送信アンテナを介して送信する。なお、この発明では、高周波数帯の電波を利用することが好ましい。 The transmitter generates radio waves (also referred to as "millimeter waves") having a frequency of, for example, a 79 GHz band, a 60 GHz band, or a 76 GHz band, and transmits the radio waves via a transmitting antenna. Note that in the present invention, it is preferable to use radio waves in a high frequency band.

レーダ装置の受信部としては、例えば、送信アンテナから送信された電波(即ち、送信波)が検知対象エリア内を移動する物標の表面で反射された電波(即ち、反射波;「ドップラ反射波」とも呼ばれる)を含む電波を受信波として受信する受信アンテナ、送信部から供給される電波(送信信号)と受信アンテナから出力される電波(受信信号)とをミキシングして差分信号を生成して出力するミキサ、および、ミキサから出力される差分信号に対して所定のサンプリング周波数を用いてサンプリング処理(別言すると、アナログ-デジタル変換処理)を施して前記差分信号をデジタルデータに変換してデジタル信号を出力するA/D変換器などを含む仕組みが挙げられる。 For example, the receiving section of a radar device receives radio waves (i.e., transmitted waves) transmitted from a transmitting antenna and radio waves (i.e., reflected waves; "Doppler reflected waves") that are reflected on the surface of a target moving within a detection target area. A receiving antenna that receives radio waves (also known as The output mixer and the differential signal output from the mixer undergo sampling processing (in other words, analog-to-digital conversion processing) using a predetermined sampling frequency to convert the differential signal into digital data. One example is a mechanism that includes an A/D converter that outputs a signal.

なお、FMCW方式では、単数もしくは複数のアンテナによってアレーアンテナが形成され、受信アンテナとしてのアレーアンテナによって電波(反射波)を受信する。そして、ミキサおよびA/D変換器は、受信アンテナごとに設けられる。 Note that in the FMCW method, an array antenna is formed by one or more antennas, and radio waves (reflected waves) are received by the array antenna as a receiving antenna. A mixer and an A/D converter are provided for each receiving antenna.

レーダ装置から出力されるレーダデータとしての差分信号は、送信部から供給される電波(送信信号)の周波数成分と受信アンテナから出力される電波(受信信号)の周波数成分との差の周波数成分を有する信号(つまり、ビート周波数を有する信号であり、「ビート信号」とも呼ばれる)である。 A difference signal as radar data output from a radar device is a difference signal between the frequency component of the radio wave (transmission signal) supplied from the transmitter and the frequency component of the radio wave (reception signal) output from the receiving antenna. (that is, a signal having a beat frequency, also called a "beat signal").

レーダ装置は、例えば、あくまで一例として挙げると、自動車などの車両に取り付けられて車両周囲の所定範囲を検知対象エリアとしたり、建物などの構造物に取り付けられて構造体周囲や建物敷地内の所定範囲を検知対象エリアとしたりする。 For example, a radar device may be attached to a vehicle such as a car to detect a predetermined area around the vehicle, or may be attached to a structure such as a building to detect a predetermined area around the structure or within the premises of the building. Set the range as the detection target area.

信号処理部2は、レーダ装置から出力されるレーダデータ(即ち、レーダ生データ、レーダRawデータ)の周波数解析を行うとともに物標の位置および相対速度を計算する。信号処理部2は、周波数解析部21、位置計算部22、および速度計算部23を含んで構成され、レーダ装置から出力されるレーダデータとしての差分信号に対して周波数解析を施して、検知対象エリア内を移動する物標を検知して前記物標の位置および相対速度を計算する。 The signal processing unit 2 performs frequency analysis of radar data (namely, radar raw data, radar raw data) output from the radar device, and calculates the position and relative velocity of the target. The signal processing unit 2 includes a frequency analysis unit 21, a position calculation unit 22, and a velocity calculation unit 23, and performs frequency analysis on the difference signal as radar data output from the radar device to determine the detection target. A target moving within an area is detected and the position and relative velocity of the target are calculated.

属性の識別装置1は、レーダ装置から出力されるレーダデータを記録して格納しておく記憶部を備えて前記記憶部からレーダデータを読み込んで用いるようにしてもよく、また、レーダ装置から出力されてサーバなどの外部記憶装置に格納されているレーダデータを読み込んで用いるようにしてもよい。 The attribute identification device 1 may include a storage section for recording and storing radar data output from the radar device, and read radar data from the storage section for use. Radar data stored in an external storage device such as a server may be read and used.

周波数解析部21は、レーダ装置(具体的には、受信部のA/D変換器)から出力されるレーダデータ(別言すると、波形データ)としての差分信号(別言すると、ビート信号;尚、デジタル信号である)の周波数解析を行う。周波数解析部21は、具体的には、1回のレーダスキャンあたりのサンプリング時間幅(例えば、50ミリ秒間)ぶんの波形データ、具体的には差分信号(ビート信号)の振幅(即ち、受信レベル)に対して高速フーリエ変換処理を施して前記差分信号の振幅の周波数分布を示す周波数スペクトルを生成し、スペクトル強度および位相情報を出力する。周波数スペクトルは、差分信号に含まれる各周波数成分の振幅(受信レベル)を示す。 The frequency analysis unit 21 extracts a difference signal (in other words, a beat signal; in other words, a beat signal; , which is a digital signal). Specifically, the frequency analysis unit 21 analyzes the waveform data for the sampling time width (for example, 50 milliseconds) per one radar scan, specifically the amplitude (i.e., the reception level) of the difference signal (beat signal). ) is subjected to fast Fourier transform processing to generate a frequency spectrum indicating the frequency distribution of the amplitude of the difference signal, and spectral intensity and phase information are output. The frequency spectrum indicates the amplitude (reception level) of each frequency component included in the difference signal.

高速フーリエ変換処理に際しては、検知対象とする物標の例えば移動速度が考慮されるなどした上で、前記物標の移動によって発生し得る周波数の範囲が設定されるとともに周波数間隔がドップラ周波数として設定される。周波数解析部21は、つまり、レーダ装置から出力されるレーダデータ(波形データ)についてドップラ成分、すなわち差分信号のドップラ周波数ごとのスペクトル強度と位相とを周波数解析して求める。 When performing fast Fourier transform processing, for example, the moving speed of the target object to be detected is taken into account, and a range of frequencies that can occur due to the movement of the target object is set, and a frequency interval is set as a Doppler frequency. be done. In other words, the frequency analysis unit 21 performs frequency analysis to obtain the Doppler component, that is, the spectral intensity and phase for each Doppler frequency of the difference signal, for the radar data (waveform data) output from the radar device.

位置計算部22は、周波数解析部21から出力される周波数スペクトルにおける位相情報に基づいて、レーダ装置(具体的には、レーダ装置の受信アンテナ)と検知対象の物標との間の距離に関する距離情報、および、レーダ装置(具体的には、レーダ装置の受信アンテナ)に対する前記物標の方位角度に関する角度情報を生成して出力する。距離情報および角度情報は、ドップラ周波数ごとに求められる。 The position calculation unit 22 calculates the distance between the radar device (specifically, the receiving antenna of the radar device) and the target object to be detected, based on the phase information in the frequency spectrum output from the frequency analysis unit 21. information and angle information regarding the azimuth angle of the target with respect to the radar device (specifically, the receiving antenna of the radar device). Distance information and angle information are obtained for each Doppler frequency.

差分信号の位相情報に基づく距離や方位角度の計算の仕法は、公知の手法が存在し、また、この発明では特定の手法には限定されないので、ここでは詳細の説明は省略する。例えば、レーダ装置と物標との間の距離は、送信波の周波数と受信波の周波数との差がレーダ装置と物標との間の距離に比例して増減する、ことを利用する手法などの公知の手法によって計算され得る。また、レーダ装置に対する物標の方位角度は、具体的には例えば、MUSIC(MUltiple SIgnal Classification の略)法やESPRIT(Estimation of Signal Parameters via Rotational Invariance Techniques の略)法などの公知の測角方式によって計算され得る。 There are known methods for calculating the distance and azimuth angle based on the phase information of the differential signal, and since the present invention is not limited to any particular method, detailed explanation will be omitted here. For example, a method that uses the fact that the distance between the radar device and the target increases or decreases in proportion to the difference between the frequency of the transmitted wave and the frequency of the received wave can be calculated by a known method. In addition, the azimuth angle of the target with respect to the radar device is specifically determined by a known angle measurement method such as the MUSIC (abbreviation for MUltiple Signal Classification) method or the ESPRIT (abbreviation for Estimation of Signal Parameters via Rotational Invariance Techniques) method. can be calculated.

速度計算部23は、周波数解析部21における周波数解析で用いられたドップラ周波数のそれぞれを用いて、レーダ装置(具体的には、レーダ装置の受信アンテナ)に対する物標の速度である相対速度(具体的には、送信波が物標の表面で反射した時の瞬時の相対速度)を計算して出力する。すなわち、相対速度は、ドップラ周波数ごとに求められる。 The velocity calculation unit 23 uses each of the Doppler frequencies used in the frequency analysis in the frequency analysis unit 21 to calculate the relative velocity (specifically, the velocity of the target with respect to the radar device (specifically, the receiving antenna of the radar device)). Specifically, it calculates and outputs the instantaneous relative velocity when the transmitted wave reflects off the surface of the target object. That is, the relative velocity is determined for each Doppler frequency.

相対速度の計算の仕法は、公知の手法が存在し、また、この発明では特定の手法には限定されないので、ここでは詳細の説明は省略する。例えば、相対速度は、レーダに対して相対的に移動している物標の表面で反射して受信アンテナによって受信される電波(即ち、受信波)の周波数は、送信アンテナから送信される電波(即ち、送信波)の周波数に対して、前記物標の表面で反射した際に受信波が前記物標の速度による影響を受け、ドップラ効果により、前記物標とレーダ装置との間の相対速度に応じてシフトする、ことを利用する手法などの公知の手法によって計算され得る。 There are known methods for calculating the relative velocity, and the present invention is not limited to any particular method, so a detailed explanation will be omitted here. For example, the relative speed is determined by the frequency of the radio waves reflected from the surface of a target moving relative to the radar and received by the receiving antenna (i.e., the received waves), and the frequency of the radio waves transmitted from the transmitting antenna (i.e., In other words, when reflected from the surface of the target object, the received wave is influenced by the speed of the target object relative to the frequency of the transmitted wave), and due to the Doppler effect, the relative speed between the target object and the radar device is reduced. It can be calculated by a known method such as a method that utilizes a shift according to .

信号処理部2は、上記の処理により、レーダ装置から出力されるレーダデータ(即ち、レーダ生データ、レーダRawデータ)に対して信号処理/周波数解析を施して、速度成分ごとの、直交座標系における強度分布を示した平面図であるPPI(Plan Position Indicator の略;平面図表示)を生成し、4次元データ(X,Y,I,V)を出力する。なお、XおよびYは、x軸とy軸とが相互に直交する2次元直交座標系におけるx座標およびy座標(単位:m)である。Iは信号強度(受信レベル)(単位:dB)であり、Vは相対速度(単位:km/h)である。信号処理部2から出力される4次元データ(X,Y,I,V)の相対速度V[km/h]は、必要に応じ、相対速度Vの値が例えば0.1[km/h]などの所定のピッチで区分された上で相対速度Vの実際の値が前記区分に当てはめられて、相対速度のランクとして出力されるようにしてもよい。 Through the above processing, the signal processing unit 2 performs signal processing/frequency analysis on the radar data (i.e., radar raw data, radar raw data) output from the radar device, and generates a rectangular coordinate system for each velocity component. A PPI (abbreviation for Plan Position Indicator; plan view display), which is a plan view showing the intensity distribution in the area, is generated, and four-dimensional data (X, Y, I, V) is output. Note that X and Y are the x and y coordinates (unit: m) in a two-dimensional orthogonal coordinate system in which the x and y axes are orthogonal to each other. I is the signal strength (reception level) (unit: dB), and V is the relative speed (unit: km/h). The relative velocity V [km/h] of the four-dimensional data (X, Y, I, V) outputted from the signal processing unit 2 may be changed, if necessary, to a value such that the value of the relative velocity V is 0.1 [km/h], for example. The actual value of the relative velocity V may be applied to the division by a predetermined pitch such as , and then output as a relative velocity rank.

物標位置検出部3は、信号処理部2における処理結果を用いて、物標としての移動している人を検出するとともに前記移動している人に関する情報を計算する。物標位置検出部3は、信号点抽出部31、クラスタ生成部32、追尾処理部33、および代表点特定部34を含んで構成され、信号処理部2から出力される、信号処理部2における処理結果である4次元データ(X,Y,I,V)を用いて、検知対象エリア内を移動している人を検出して前記移動している人に関する情報を計算する。 The target position detection unit 3 uses the processing results in the signal processing unit 2 to detect a moving person as a target and calculates information regarding the moving person. The target position detection section 3 includes a signal point extraction section 31, a cluster generation section 32, a tracking processing section 33, and a representative point specification section 34. Using the four-dimensional data (X, Y, I, V) that is the processing result, a person moving within the detection target area is detected and information regarding the moving person is calculated.

信号点抽出部31は、信号処理部2の周波数解析部21によって生成される周波数スペクトルを用いて、所定の大きさの振幅を有するピークに対応するドップラ周波数(別言すると、ピーク周波数)を検出して、検知対象エリア内を移動する物標に関するレーダ信号点の周波数として抽出する。具体的には、信号点抽出部31は、ピークの検出のための閾値を設定し、周波数スペクトルの振幅について、前記閾値以上の振幅をピークと判定して抽出する。 The signal point extraction unit 31 uses the frequency spectrum generated by the frequency analysis unit 21 of the signal processing unit 2 to detect a Doppler frequency (in other words, a peak frequency) corresponding to a peak having a predetermined amplitude. Then, it is extracted as the frequency of a radar signal point related to a target moving within the detection target area. Specifically, the signal point extracting unit 31 sets a threshold for detecting a peak, and extracts the amplitude of the frequency spectrum by determining that an amplitude equal to or greater than the threshold is a peak.

ピークの検出・抽出の仕法は、特定の方法・方式に限定されるものではなく、クラッタと物標とを区別して物標からの反射成分と考えられる周波数を検出・抽出することができればどのような方法・方式であってもよい。ピークの検出・抽出の仕法として、具体的には例えばCFAR(Constant False Alarm Rate の略)アルゴリズムが用いられるようにしてもよい。 The method of detecting and extracting peaks is not limited to a specific method or method, but it is possible to distinguish between clutter and a target object and detect and extract the frequency that is considered to be the reflected component from the target object. Any method or system may be used. Specifically, for example, a CFAR (abbreviation for Constant False Alarm Rate) algorithm may be used as a method for detecting and extracting peaks.

信号点抽出部31は、また、位置計算部22から出力されるドップラ周波数ごとの距離情報および角度情報を用いて、空間位置を座標軸とする座標系における、上記において検出されたドップラ周波数(ピーク周波数)に対応するレーダ信号点の分布に関する分布情報を生成して出力する。座標系としては、例えば、レーダ装置の設置位置を基準とした、動径距離と動径方向とについての極座標系、または、左右方向における距離と奥行方向における距離とについての直交座標系が用いられる。 The signal point extraction unit 31 also uses the distance information and angle information for each Doppler frequency output from the position calculation unit 22 to calculate the Doppler frequency (peak frequency) detected above in the coordinate system with the spatial position as the coordinate axis. ) generates and outputs distribution information regarding the distribution of radar signal points corresponding to As the coordinate system, for example, a polar coordinate system for the radial distance and the radial direction, or a rectangular coordinate system for the distance in the left-right direction and the distance in the depth direction is used, based on the installation position of the radar device. .

クラスタ生成部32は、信号点抽出部31から出力されるレーダ信号点の分布情報を用いてクラスタリング処理を行い、空間位置を座標軸とする座標系における、レーダ信号点のクラスタを生成する。 The cluster generation unit 32 performs clustering processing using the distribution information of the radar signal points output from the signal point extraction unit 31, and generates clusters of radar signal points in a coordinate system having the spatial position as the coordinate axis.

クラスタの生成の仕法は、特定の方法・方式に限定されるものではなく、例えば、信号点抽出部31によって検出・抽出されたレーダ信号点の密度に基づくなどして物標についてのレーダ信号点の集まりをクラスタリングする、言い換えると、相互の位置関係が一定距離以内で近接しているレーダ信号点の集まりを括って抽出することができればどのような方法・方式であってもよい。クラスタの生成の仕法として、具体的には例えばDBSCAN(Density-Based Spatial Clustering of Applications with Noise の略)がクラスタリングアルゴリズムとして用いられるようにしてもよい。クラスタ生成部32によって生成されるクラスタのそれぞれが1つの物標に対応する。 The method of cluster generation is not limited to a specific method or method. For example, the radar signal points about the target are generated based on the density of radar signal points detected and extracted by the signal point extraction unit 31. Any method or system may be used as long as it can cluster a collection of radar signal points, in other words, collectively extract a collection of radar signal points that are close to each other within a certain distance relative to each other. Specifically, as a method for generating clusters, for example, DBSCAN (abbreviation for Density-Based Spatial Clustering of Applications with Noise) may be used as a clustering algorithm. Each cluster generated by the cluster generation unit 32 corresponds to one target.

クラスタ生成部32によって生成される各クラスタを区別して一意に特定するための識別符号をNとする(Nは、例えばシリアル番号)。そして、クラスタNについての、当該のクラスタNに含まれるレーダ信号点それぞれの相対速度(即ち、信号処理部2の速度計算部23によって計算される各ドップラ周波数に対応するレーダ信号点の相対速度)の値の集まりのことを「クラスタNの相対速度の値の集合」と呼ぶ。 Let N be an identification code for distinguishing and uniquely identifying each cluster generated by the cluster generation unit 32 (N is, for example, a serial number). Then, with respect to the cluster N, the relative speed of each radar signal point included in the cluster N (that is, the relative speed of the radar signal point corresponding to each Doppler frequency calculated by the speed calculation unit 23 of the signal processing unit 2) The collection of values of is called "the collection of relative velocity values of cluster N".

なお、物標位置検出部3は、検知対象エリアで検知された物標が、延いてはクラスタ生成部32によって生成されたクラスタが、移動している人に関するものであるか否かを判別・識別した上で、移動している人に関するクラスタに対してのみ識別符号Nを付するようにしてもよい。検知された物標が移動している人であるか否かを判別・識別したり、生成されたクラスタが移動している人に関するものであるか否かを判別・識別したりする仕法は、公知の手法が存在し、また、この発明では特定の手法には限定されないので、ここでは詳細の説明は省略する。 Note that the target position detection unit 3 determines whether or not the target detected in the detection target area, and by extension the cluster generated by the cluster generation unit 32, is related to a moving person. After identification, the identification code N may be attached only to clusters related to moving people. The method of determining and identifying whether a detected target is a moving person, and determining and identifying whether a generated cluster is related to a moving person is as follows. Since there are known methods and the present invention is not limited to any particular method, detailed explanation will be omitted here.

追尾処理部33は、過去に検出されている物標を追尾する処理を、クラスタ生成部32によって生成されたクラスタNごと(即ち、物標ごと)に行う。つまり、追尾処理部33は、連続する時点においてクラスタ生成部32によって生成される、同一の物標に関するクラスタを関連づけることにより、過去のレーダスキャンによって検出されて識別符号Nが付されているクラスタの追尾処理を行う。 The tracking processing unit 33 performs a process of tracking targets detected in the past for each cluster N generated by the cluster generation unit 32 (that is, for each target). In other words, the tracking processing unit 33 associates clusters related to the same target object generated by the cluster generation unit 32 at successive points in time, thereby identifying the clusters detected by past radar scans and assigned the identification code N. Perform tracking processing.

クラスタ/物標の追尾の仕法は、公知の手法が存在し、また、この発明では特定の手法には限定されないので、ここでは詳細の説明は省略する。例えば、クラスタ/物標の追尾は、追尾フィルタとしてα-βフィルタ、カルマンフィルタ、或いはパーティクルフィルタを用いる手法などの公知の手法によって行われ得る。 There are known cluster/target tracking methods, and the present invention is not limited to any particular method, so a detailed explanation will be omitted here. For example, cluster/target object tracking may be performed by a known method such as a method using an α-β filter, a Kalman filter, or a particle filter as a tracking filter.

追尾処理部33は、必要に応じ、クラスタ生成部32によって生成されるクラスタそれぞれの代表点を決定した上で追尾処理を行うようにしてもよい。例えば、追尾処理部33は、1つのクラスタとして認識された信号点群について、当該の信号点群に含まれる信号点それぞれの位置の、信号強度(受信レベル)I[dB]を重みとして用いた加重平均の位置を計算するなどして、前記1つのクラスタに含まれる複数の信号点を1点に集約した上で追尾処理を行うようにしてもよい。 The tracking processing unit 33 may perform the tracking process after determining representative points of each cluster generated by the cluster generation unit 32, if necessary. For example, for a group of signal points recognized as one cluster, the tracking processing unit 33 uses the signal strength (reception level) I [dB] of the position of each signal point included in the group of signal points as a weight. The tracking process may be performed after a plurality of signal points included in the one cluster are aggregated into one point by calculating a weighted average position.

代表点特定部34は、クラスタ生成部32によって生成されるとともに追尾処理部33によって追尾されるクラスタNごとに、各レーダスキャンにおける、物標としてのクラスタNの代表点を特定する。代表点特定部34は、各クラスタNの例えば重心位置を代表点の位置として特定して、クラスタNごとに特定した代表点の位置座標X[m],Y[m]の組み合わせを出力する。 For each cluster N generated by the cluster generation unit 32 and tracked by the tracking processing unit 33, the representative point identification unit 34 identifies a representative point of the cluster N as a target in each radar scan. The representative point specifying unit 34 specifies, for example, the position of the center of gravity of each cluster N as the position of the representative point, and outputs a combination of position coordinates X[m] and Y[m] of the specified representative point for each cluster N.

物標計算部4は、物標位置検出部3における処理結果を用いて、特徴量計算部5による次の処理における関心領域を設定するとともに前記関心領域に関する情報を計算する。物標計算部4は、関心領域設定部41、信号強度特定部42、およびVTデータ生成部43を含んで構成され、物標位置検出部3から出力される、物標位置検出部3における処理結果を用いて、特徴量計算部5による次の処理における関心領域を設定して前記関心領域に関する情報を計算する。 The target calculation unit 4 uses the processing results in the target position detection unit 3 to set a region of interest for the next processing by the feature amount calculation unit 5 and calculates information regarding the region of interest. The target calculation unit 4 includes a region of interest setting unit 41, a signal strength identification unit 42, and a VT data generation unit 43, and processes output from the target position detection unit 3 in the target position detection unit 3. Using the results, a region of interest is set for the next process by the feature calculation unit 5, and information regarding the region of interest is calculated.

関心領域設定部41は、クラスタNごとに、物標位置検出部3の代表点特定部34から出力される各クラスタNの代表点の位置を中心とする所定の範囲を関心領域として設定する。関心領域として設定される範囲は、特定の大きさや形状に限定されるものではなく、検知対象の物標の特に動作速度に纏わる特性を適切に把握することが考慮されるなどした上で、適当な大きさや形状に適宜設定される。関心領域は、具体的には例えば、クラスタNの代表点の位置座標X[m],Y[m]を中心とする±1[m]の四角形で囲まれる範囲に設定されることが考えられる。 For each cluster N, the region of interest setting section 41 sets a predetermined range centered on the position of the representative point of each cluster N output from the representative point specifying section 34 of the target position detection section 3 as a region of interest. The range set as the region of interest is not limited to a specific size or shape, but may be set as appropriate after taking into account the proper understanding of the characteristics of the target object to be detected, especially related to the movement speed. The size and shape are set as appropriate. Specifically, the region of interest may be set, for example, to a range surrounded by a rectangle of ±1 [m] centered on the position coordinates X [m], Y [m] of the representative point of cluster N. .

関心領域の大きさや形状は、予め固定的に設定されて関心領域設定部41に格納されているようにしてもよく、或いは、関心領域設定部41に対してユーザが適宜設定可能であるようにしてもよい。 The size and shape of the region of interest may be fixedly set in advance and stored in the region of interest setting section 41, or may be set appropriately by the user in the region of interest setting section 41. It's okay.

関心領域設定部41は、さらに、信号処理部2から出力される4次元データ(X,Y,I,V)の入力を受け、前記4次元データ(X,Y,I,V)のうち、関心領域として設定されている範囲に含まれているデータのみを抜き出す。関心領域設定部41は、つまり、関心領域として設定されている範囲に含まれている、全ての速度成分についての強度データを出力する。 The region of interest setting unit 41 further receives input of four-dimensional data (X, Y, I, V) output from the signal processing unit 2, and selects among the four-dimensional data (X, Y, I, V). Extract only the data included in the range set as the region of interest. In other words, the region of interest setting unit 41 outputs intensity data for all velocity components included in the range set as the region of interest.

信号強度特定部42は、関心領域設定部41から出力される関心領域に関する4次元データ(X,Y,I,V)を速度成分ごとに参照し、速度成分別に、関心領域に関する強度情報を生成する。信号強度特定部42は、具体的には、クラスタNの関心領域に関する4次元データ(X,Y,I,V)について、相対速度V[km/h]ごとに、例えば、信号強度(受信レベル)I[dB]のうちの最大値Imaxを特定して関心領域に関する強度情報としたり、信号強度I[dB]の平均値Iavgを計算して関心領域に関する強度情報としたりする。 The signal strength specifying unit 42 refers to the four-dimensional data (X, Y, I, V) regarding the region of interest outputted from the region of interest setting unit 41 for each velocity component, and generates intensity information regarding the region of interest for each velocity component. do. Specifically, the signal strength specifying unit 42 determines, for example, the signal strength (reception level ) The maximum value Imax of I[dB] is specified and used as intensity information regarding the region of interest, or the average value Iavg of the signal intensity I[dB] is calculated and used as intensity information regarding the region of interest.

ここで、信号強度特定部42によって特定されたり計算されたりする所定の特性をもった信号強度I[km/h]のことを「信号強度情報」と呼ぶ。つまり、信号強度特定部42は、クラスタNごとの、関心領域に関する、相対速度Vの値と信号強度情報の値(例えば、最大値Imax,平均値Iavg[dB])との組み合わせデータを出力する。 Here, the signal strength I [km/h] having a predetermined characteristic that is specified or calculated by the signal strength specifying unit 42 is referred to as "signal strength information." That is, the signal strength specifying unit 42 outputs combination data of the value of relative velocity V and the value of signal strength information (for example, maximum value Imax, average value Iavg [dB]) regarding the region of interest for each cluster N. .

VTデータ生成部43は、信号強度特定部42から出力される関心領域に関する信号強度情報を用いて、各レーダスキャンの実行時刻に基づく経過時間T[秒]および相対速度V[km/h]と、これら経過時間T[秒]と相対速度V[km/h]との組み合わせに対応する信号強度情報(例えば、最大値Imax,平均値Iavg[dB])と、の組み合わせデータを生成する。 The VT data generation unit 43 uses the signal intensity information regarding the region of interest output from the signal intensity identification unit 42 to calculate the elapsed time T [seconds] and relative velocity V [km/h] based on the execution time of each radar scan. , and signal strength information (for example, maximum value Imax, average value Iavg [dB]) corresponding to the combination of elapsed time T [seconds] and relative velocity V [km/h].

信号処理部2、物標位置検出部3、および物標計算部4による処理は、レーダスキャンごとに行われる。すなわち、1回のレーダスキャンが行われるたびに、クラスタNごとの、関心領域に関する、経過時間T[秒]と相対速度V[km/h]と信号強度情報(例えば、最大値Imax,平均値Iavg[dB])との組み合わせデータが作成される。 Processing by the signal processing unit 2, target position detection unit 3, and target calculation unit 4 is performed for each radar scan. That is, each time one radar scan is performed, elapsed time T [seconds], relative velocity V [km/h], and signal strength information (for example, maximum value Imax, average value Iavg [dB]) is created.

特徴量計算部5は、物標計算部4における処理結果を用いて、物標としての移動している人の属性を識別するための特徴量を計算する。特徴量計算部5は、相関計算部51、周期計算部52、振幅計算部53、および相関変換部54を含んで構成され、物標計算部4から出力される、前記物標計算部4おける処理結果を用いて、検知対象エリア内を移動している人に関する特徴量を計算する。下記において説明する特徴量計算部5による処理は、その都度断らないが、クラスタNごとに(言い換えると、クラスタN各々を単位として)行われる。 The feature amount calculation section 5 uses the processing results in the target object calculation section 4 to calculate a feature amount for identifying the attributes of the moving person as the target object. The feature calculation section 5 includes a correlation calculation section 51, a period calculation section 52, an amplitude calculation section 53, and a correlation conversion section 54, Using the processing results, feature quantities related to people moving within the detection target area are calculated. The processing by the feature value calculation unit 5 described below is performed for each cluster N (in other words, each cluster N is used as a unit), although this is not specified in each case.

相関計算部51は、物標計算部4のVTデータ生成部43から出力される経過時間T[秒]と相対速度V[km/h]との組み合わせに対応する信号強度情報(例えば、最大値Imax,平均値Iavg[dB])のデータを用いて相関分析を行う。 The correlation calculation unit 51 calculates signal strength information (for example, maximum value Correlation analysis is performed using data of Imax, average value Iavg [dB]).

相関計算部51は、具体的には、まず、所定の時間長さぶん(言い換えると、レーダスキャンの所定の回数ぶん)のレーダスキャンごとの相対速度Vの値と信号強度情報の値との組み合わせデータの集合を時系列(言い換えると、レーダスキャンの実行順)に配列する。時系列で配列された、レーダスキャンごとの相対速度Vの値と信号強度情報の値との組み合わせデータの集合のことを「速度・強度の時系列データ」と呼ぶ。 Specifically, the correlation calculation unit 51 first calculates a combination of the relative velocity V value and the signal strength information value for each radar scan for a predetermined time length (in other words, a predetermined number of radar scans). Arrange the data set in chronological order (in other words, in the order in which radar scans are performed). A set of combination data of the value of relative velocity V for each radar scan and the value of signal strength information arranged in time series is called "time series data of velocity/strength".

相関計算部51から出力される速度・強度の時系列データの一例を図3に示す。図3は、各レーダスキャンの実行時刻に基づく経過時間T[秒](即ち、1回のレーダスキャンあたりのサンプリング時間幅(例えば、50ミリ秒間)とレーダスキャンの実行回数とを掛け合わせたもの)を横軸とするとともに相対速度V[km/h]の値を縦軸として、経過時間T[秒]と相対速度V[km/h]との組み合わせに対応する信号強度(受信レベル)I[dB]のうちの最大値Imax[dB]を明暗で表して表示したものである。図3では、最大信号強度が、弱い場合は暗く表示され、強い場合は明るく表示されている。速度・強度の時系列データは、時系列における、物標に関する瞬時の速度の分布情報であるとも言える。なお、図3および図5~図16は、2.5秒ぶん(言い換えると、50回分のレーダスキャン)のデータが用いられて整理されている。 FIG. 3 shows an example of speed/intensity time series data output from the correlation calculation unit 51. Figure 3 shows the elapsed time T [seconds] based on the execution time of each radar scan (i.e., the product of the sampling time width (for example, 50 milliseconds) per one radar scan and the number of times the radar scan is executed). ) as the horizontal axis and the value of relative speed V [km/h] as the vertical axis, the signal strength (reception level) I corresponding to the combination of elapsed time T [seconds] and relative speed V [km/h] The maximum value Imax [dB] of [dB] is displayed in brightness and darkness. In FIG. 3, when the maximum signal strength is weak, it is displayed darkly, and when it is strong, it is displayed brightly. The speed/intensity time series data can also be said to be instantaneous speed distribution information regarding the target in a time series. Note that FIG. 3 and FIGS. 5 to 16 are organized using data for 2.5 seconds (in other words, 50 radar scans).

ここで、レーダ装置を用いて得られる相対速度の情報は、物標としての移動している人の胴体速度だけでなく、頭、上腕、前腕・手先、大腿、および下腿・足先などの様々な部位の速度成分を同時に取得している。したがって、移動している人の頭、胴体、上腕、前腕・手先、大腿、および下腿・足先などの各部位の速度の情報が内在することになるため、レーダスキャンごとの相対速度Vの値と信号強度情報の値(図3に示す例では、最大信号強度の値)との組み合わせデータを時系列に並べることにより、各部位の運動/動作の周期が重畳して現れるようになる。 Here, the information on the relative speed obtained using a radar device is not only the speed of the torso of a moving person as a target, but also the speed of the head, upper arm, forearm/hand, thigh, lower leg/feet, etc. The velocity components of different parts are acquired simultaneously. Therefore, since information on the speed of each part of a moving person, such as the head, torso, upper arm, forearm/hand, thigh, and lower leg/feet, is included, the value of relative velocity V for each radar scan is included. By arranging the combination data of the signal strength information and the signal strength information value (in the example shown in FIG. 3, the maximum signal strength value) in chronological order, the motion/motion cycles of each part appear superimposed.

相関計算部51は、続いて、速度・強度の時系列データについて、レーダスキャンごとの相対速度Vの値と信号強度情報の値との組み合わせデータの時系列の配列を時間軸方向に1つずつシフトさせながら内積を計算して、シフト前の、時系列で配列された組み合わせデータの集合との対比における相関値を計算する(図4参照)。組み合わせデータの時系列の配列を時間軸方向に1つシフトさせる際、時系列の最後の時点に配置されている組み合わせデータは時系列の最初の時点へと配置される(図4参照)。相関値の計算は、時系列で配列された組み合わせデータの数と同じだけシフトを繰り返して行われる。なお、図4に示す例では、時系列で配列された組み合わせデータの数が50であり、シフト実行回が50回目のデータは、シフト前の、時系列で配列された組み合わせデータの集合と一致する。 Next, regarding the speed/intensity time series data, the correlation calculation unit 51 arranges the time series of the combination data of the relative speed V value and the signal strength information value for each radar scan one by one in the time axis direction. The inner product is calculated while shifting, and the correlation value in comparison with the set of combined data arranged in time series before the shift is calculated (see FIG. 4). When shifting the time series arrangement of the combined data by one in the time axis direction, the combination data placed at the last point in the time series is placed at the first point in the time series (see FIG. 4). Calculation of the correlation value is performed by repeating shifts as many times as the number of combined data arranged in time series. In the example shown in FIG. 4, the number of combination data arranged in chronological order is 50, and the data of the 50th shift execution matches the set of combination data arranged in chronological order before the shift. do.

相関計算部51は、速度・強度の時系列データを構成する、レーダスキャンごとの相対速度Vの値と信号強度情報の値との組み合わせデータを時間軸方向にシフトさせた回(「データのシフト回」とも表記する)ごとの相関値を出力する。相関計算部51によって計算される、クラスタNについてのデータのシフト回ごとの相関値の集まりのことを「クラスタNの相関値の集合」と呼ぶ。 The correlation calculation unit 51 calculates the times when the combination data of the value of relative velocity V and the value of signal strength information for each radar scan, which constitutes the time series data of velocity/intensity, is shifted in the time axis direction ("data shift"). outputs the correlation value for each A collection of correlation values calculated by the correlation calculation unit 51 for each shift of data for cluster N is referred to as a "set of correlation values of cluster N."

相関計算部51から出力されるクラスタNの相関値の集合のデータの一例を図5に示す。図5は、データのシフト回を横軸とするとともに相関値を縦軸として、データのシフト回ごとの相関値の変動を折れ線で表して表示したものである。図5では、物標の運動の周期が重なったときに相関値のピークが現れる。 FIG. 5 shows an example of data of a set of correlation values of cluster N output from the correlation calculation unit 51. In FIG. 5, the horizontal axis is the number of data shifts, and the vertical axis is the correlation value, and the variation in the correlation value for each data shift is displayed using a polygonal line. In FIG. 5, the peak of the correlation value appears when the motion cycles of the targets overlap.

周期計算部52は、相関計算部51から出力されるクラスタNの相関値の集合を用いて、クラスタNの相関値の集合を時系列(即ち、速度・強度の時系列データについての時間軸方向のシフトの実行順)に配列したときの、相関値のピーク間の差分の平均値を計算して出力する。 The period calculation unit 52 uses the set of correlation values of the cluster N output from the correlation calculation unit 51 to convert the set of correlation values of the cluster N into a time series (i.e., in the time axis direction for time series data of speed and intensity). Calculates and outputs the average value of the difference between the peaks of the correlation values when they are arranged in the order of execution of the shift).

クラスタNの相関値の集合(即ち、時系列で配列された相対速度Vの値と信号強度情報の値との組み合わせデータを時間軸方向に1つずつシフトさせながら、シフト前の、時系列で配列された組み合わせデータの集合との対比における相関値の集まり)における相関値のピーク間の差分D1,D2,D3は、物標の動作の周期性を表す指標である(図6参照)。すなわち、相関値のピーク間の差分の平均値により、物標の運動周期(ここでは、移動している人の歩行に纏わる動作の周期)の遅速を把握することが可能となる。なお、シフト実行回が50回目のデータは、シフト前の、時系列で配列された組み合わせデータの集合と一致するので相関値が1であり、周期計算部52による処理においては相関値のピークとしては扱わない。 While shifting the set of correlation values of cluster N (that is, the combination data of relative velocity V values and signal strength information values arranged in time series) one by one in the time axis direction, Differences D1, D2, and D3 between the peaks of the correlation values in the collection of correlation values compared to the set of arranged combination data are indicators representing the periodicity of the movement of the target object (see FIG. 6). That is, the average value of the difference between the peaks of the correlation values makes it possible to grasp the slowness of the motion cycle of the target (here, the cycle of motion related to the walking of a moving person). Note that the data of the 50th shift execution matches the set of combination data arranged in time series before the shift, so the correlation value is 1, and in the processing by the period calculation unit 52, it is regarded as the peak of the correlation value. is not treated.

振幅計算部53は、相関計算部51から出力されるクラスタNの相関値の集合を用いて、クラスタNの相関値の集合を時系列(即ち、速度・強度の時系列データについての時間軸方向のシフトの実行順)に配列したときの、相関値の極大値および極小値の平均と分散とを計算して出力する。 The amplitude calculation unit 53 uses the set of correlation values of the cluster N output from the correlation calculation unit 51 to convert the set of correlation values of the cluster N into a time series (i.e., in the time axis direction for time series data of speed and intensity). The average and variance of the maximum and minimum correlation values are calculated and output when the correlation values are arranged in the order of execution of shifts).

クラスタNの相関値の集合における相関値の極大値Ex1,Ex2,Ex3,Ex4および極小値En1,En2,En3,En4,En5は、物標の動作の振幅を表す指標である(図7参照)。すなわち、相関値の極大値および極小値の平均や分散により、物標の運動の振幅や変動(ここでは、移動している人の歩行に纏わる動作の振れ幅やばらつき)の大小を把握することが可能となる。なお、シフト実行回が50回目のデータは、シフト前の、時系列で配列された組み合わせデータの集合と一致するので相関値が1であり、振幅計算部53による処理においては相関値の極大値としては扱わない。 The local maximum values Ex1, Ex2, Ex3, Ex4 and the local minimum values En1, En2, En3, En4, En5 of the correlation values in the set of correlation values of cluster N are indicators representing the amplitude of the movement of the target object (see FIG. 7). . In other words, it is possible to understand the amplitude and variation of the movement of a target object (here, the amplitude and dispersion of movement related to the walking of a moving person) by the average and variance of the maximum and minimum correlation values. becomes possible. Note that the data of the 50th shift execution matches the set of combination data arranged in time series before the shift, so the correlation value is 1, and in the processing by the amplitude calculation unit 53, the correlation value is the maximum value. It is not treated as such.

相関変換部54は、相関計算部51から出力されるクラスタNの相関値の集合を用いて、高速フーリエ変換を行い、L1ノルムを計算して出力する。相関変換部54は、具体的には、以下の処理を行う。 The correlation transformer 54 performs fast Fourier transform using the set of correlation values of cluster N output from the correlation calculator 51, calculates and outputs the L 1 norm. Specifically, the correlation conversion unit 54 performs the following processing.

1)相関計算部51から出力される、所定の時間長さぶん(言い換えると、レーダスキャンの所定の回数ぶん)の、データのシフト回ごとの相関値の平均(「相関平均値」と呼ぶ)を算出する。
2)データのシフト回ごとの相関値と相関平均値との差(「シフト回ごとの相関差分」と呼ぶ)を算出する。
3)シフト回ごとの相関差分に対して高速フーリエ変換処理を施す。具体的には、シフト回ごとの相関値の大きさ(言い換えると、振幅)に対して高速フーリエ変換処理を施して相関値の振幅の分布を示すパワースペクトルを生成する。
4)高速フーリエ変換処理の結果の絶対値をとり、前記絶対値を2乗する。
5)上記4の結果として算出される、前記絶対値を2乗した値の集合のうちの最大値で、前記絶対値を2乗した値のそれぞれを割り、前記絶対値を2乗した値を0から1までの範囲で正規化する。
1) The average of the correlation values for each data shift for a predetermined length of time (in other words, a predetermined number of radar scans) output from the correlation calculation unit 51 (referred to as "correlation average value") Calculate.
2) Calculate the difference between the correlation value for each data shift and the correlation average value (referred to as "correlation difference for each shift").
3) Perform fast Fourier transform processing on the correlation difference for each shift. Specifically, fast Fourier transform processing is performed on the magnitude (in other words, amplitude) of the correlation value for each shift to generate a power spectrum indicating the distribution of the amplitude of the correlation value.
4) Take the absolute value of the result of the fast Fourier transform process and square the absolute value.
5) Divide each of the squared absolute values by the maximum value of the set of values obtained by squaring the absolute value calculated as a result of 4 above, and calculate the value obtained by squaring the absolute value. Normalize in the range from 0 to 1.

上記の処理により、クラスタNの相関値の集合の、高速フーリエ変換処理後のL1ノルムが求められ、相関変換部54は、高速フーリエ変換処理後のL1ノルムとして、高速フーリエ変換後の周波数と振幅との組み合わせデータを出力する(図8参照)。 Through the above processing, the L 1 norm after the fast Fourier transform of the set of correlation values of cluster N is obtained, and the correlation transform unit 54 calculates the frequency after the fast Fourier transform as the L 1 norm after the fast Fourier transform. The combination data of and amplitude is output (see FIG. 8).

上記により、特徴量計算部5は、クラスタNごとの、相関値のピーク間の差分の平均値、相関値の極大値および極小値の平均と分散、さらに高速フーリエ変換後のL1ノルムとしての周波数および振幅の組み合わせデータを出力する。 As a result of the above, the feature calculation unit 5 calculates the average value of the difference between the peaks of the correlation values, the average and variance of the maximum and minimum values of the correlation values, and the L 1 norm after fast Fourier transform for each cluster N. Outputs combined frequency and amplitude data.

識別部6は、特徴量計算部5から出力されるクラスタNごとの各種指標の組み合わせデータに基づいて、クラスタNごとに物標の属性を識別する。識別部6は、具体的には例えば、クラスタNに対応する、物標としての移動している人の性別や年代を判別する。 The identification unit 6 identifies the attributes of the target object for each cluster N based on the combination data of various indicators for each cluster N output from the feature calculation unit 5. Specifically, the identification unit 6 determines, for example, the gender and age of a person who is moving as a target and corresponds to cluster N.

人の属性の判別について、識別部6は、例えば、機械学習のアルゴリズムの1つであるサポートベクターマシンを用いて判別処理を行うようにしてもよい。 Regarding discrimination of a person's attributes, the identification unit 6 may perform discrimination processing using, for example, a support vector machine, which is one of machine learning algorithms.

この場合は、識別部6は、例えば、クラスタNの相関値の集合の、相関値のピーク間の差分の平均値x1、相関値の極大値および極小値の平均の値x2、ならびに相関値の極大値および極小値の分散の値x3を変数として含むとともに係数f1、f2、およびf3を含む式(下記の数式1参照)を用いて人の属性の判別処理を行う。数式1について、例えば、クラスタNが男性についてのレーダ信号点群である場合はy=0とされ、クラスタNが女性についてのレーダ信号点群である場合はy=1とされる。
(数1) y=f(x1,x2,x3,f1,f2,f3)
In this case, the identification unit 6 calculates, for example, the average value x1 of the difference between the peaks of the correlation values, the average value x2 of the maximum and minimum values of the correlation values, and the average value x2 of the maximum and minimum correlation values of the set of correlation values of cluster N. A process for determining a person's attributes is performed using an equation (see Equation 1 below) that includes the value x3 of the variance of the local maximum value and the local minimum value as a variable, and includes coefficients f1, f2, and f3. Regarding Equation 1, for example, when cluster N is a radar signal point group for men, y=0, and when cluster N is a radar signal point group for women, y=1.
(Math. 1) y=f(x1, x2, x3, f1, f2, f3)

数式1の係数f1、f2、およびf3は、サポートベクターマシンが用いられて予め決定される。係数f1、f2、およびf3の値は、男性についてのレーダ信号点群で構成されるクラスタNの相関値の集合の、相関値のピーク間の差分の平均値x1、相関値の極大値および極小値の平均の値x2、ならびに相関値の極大値および極小値の分散の値x3、ならびにyの値(即ち、0)の組み合わせ(複数の組み合わせ)と、女性についてのレーダ信号点群で構成されるクラスタNの相関値の集合の、相関値のピーク間の差分の平均値x1、相関値の極大値および極小値の平均の値x2、ならびに相関値の極大値および極小値の分散の値x3、ならびにyの値(即ち、1)の組み合わせ(複数の組み合わせ)と、を含む教師データが用意され、この教師データを男性についてのクラスタと女性についてのクラスタとに分離するための平面(「超平面」と呼ばれる)が最大のマージンを有するように決定される。 Coefficients f1, f2, and f3 in Equation 1 are determined in advance using a support vector machine. The values of coefficients f1, f2, and f3 are the average value x1 of the difference between the peaks of the correlation values, the maximum value and the minimum value of the correlation values of the set of correlation values of cluster N consisting of radar signal points for men. It consists of the average value x2, the variance value x3 of the maximum value and minimum value of the correlation value, and combinations (multiple combinations) of the value of y (i.e. 0), and the radar signal point group for women. The average value of the difference between the peaks of the correlation values x1, the average value of the maximum value and the minimum value of the correlation values x2, and the value of the variance of the maximum value and the minimum value of the correlation values x3 of the set of correlation values of the cluster N , and combinations (multiple combinations) of y values (i.e., 1), are prepared, and a plane (“super is determined to have the largest margin.

識別部6は、上記によって予め決定された係数f1、f2、およびf3を含む数式1に、特徴量計算部5から出力されるクラスタNの相関値のピーク間の差分の平均値x1、相関値の極大値および極小値の平均の値x2、ならびに相関値の極大値および極小値の分散の値x3を代入して求められるyの値に基づいて、前記クラスタNが、男性についてのレーダ信号点群であるのか、或いは、女性についてのレーダ信号点群であるのかを判別する。 The identification unit 6 adds the average value x1 of the difference between the peaks of the correlation values of the cluster N output from the feature quantity calculation unit 5, and the correlation value Based on the value of y obtained by substituting the average value x2 of the local maximum value and local minimum value of , and the value x3 of the variance of the local maximum value and local minimum value of the correlation value, the cluster N It is determined whether it is a group of radar signal points or a group of radar signal points about women.

ここで、識別部6は、クラスタNの相関値の集合を時系列に配列したときの、相関値のピーク間の差分の平均値x1、相関値の極大値および極小値の平均の値x2、ならびに相関値の極大値および極小値の分散の値x3のうちの少なくとも一つを用いて人の属性の判別処理を行うようにしてもよい。さらに、識別部6は、クラスタNの相関値の集合の、高速フーリエ変換処理後のL1ノルムのみを用いて人の属性の判別処理を行うようにしてもよく、或いは、上記平均値x1、平均の値x2、および分散の値x3のうちの少なくとも一つと前記L1ノルムとを組み合わせて用いて人の属性の判別処理を行うようにしてもよい。識別部6は、また、人の属性として、例えば、成人であるか子供であるかを判別するようにしてもよい。 Here, the identification unit 6 determines, when the set of correlation values of cluster N is arranged in time series, the average value x1 of the difference between the peaks of the correlation values, the average value x2 of the local maximum value and the local minimum value of the correlation values, Also, at least one of the maximum value of the correlation value and the value x3 of the variance of the minimum value may be used to perform the process of determining a person's attributes. Furthermore, the identification unit 6 may perform a process of determining a person's attributes using only the L 1 norm of the set of correlation values of cluster N after fast Fourier transformation processing, or the above-mentioned average value x1, At least one of the average value x2 and the variance value x3 and the L1 norm may be used in combination to perform the process of determining a person's attributes. The identification unit 6 may also determine whether the person is an adult or a child, for example, as an attribute of the person.

また、識別部6は、人の属性の判別について、機械学習のアルゴリズムの1つであるニューラルネットワーク、具体的には例えばRNN(Reccurent Neural Network の略)であるLSTM(Long-Short Term Memory の略)アルゴリズムを用いて判別処理を行うようにしてもよい。 In addition, the identification unit 6 uses a neural network, which is one of the machine learning algorithms, to determine human attributes, specifically, for example, RNN (recurrent neural network), LSTM (long-short term memory), etc. ) The discrimination process may be performed using an algorithm.

次に、このような構成の属性の識別装置1および属性の識別方法の動作や作用などについて図2も用いて説明する。 Next, the operation and effects of the attribute identification device 1 and the attribute identification method having such a configuration will be described with reference to FIG. 2 as well.

まず、信号処理部2の周波数解析部21が、レーダ装置から出力されるレーダデータ(別言すると、波形データ;具体的には、差分信号の振幅)としてのデジタル信号に対して高速フーリエ変換を用いて周波数解析を行って周波数スペクトルを生成し、差分信号のドップラ周波数ごとのスペクトル強度および位相情報を出力する(ステップS1)。 First, the frequency analysis unit 21 of the signal processing unit 2 performs fast Fourier transform on a digital signal as radar data (in other words, waveform data; specifically, the amplitude of a difference signal) output from a radar device. A frequency spectrum is generated by performing frequency analysis using the spectral signal, and spectral intensity and phase information for each Doppler frequency of the difference signal are output (step S1).

次に、位置計算部22が、周波数解析部21から出力される位相情報に基づいて、ドップラ周波数ごとに、レーダ装置と検知対象の物標との間の距離を計算するとともにレーダ装置に対する前記物標の方位角度を計算し、計算の結果得られる距離情報および角度情報を出力する(ステップS2)。また、速度計算部23が、周波数解析部21における周波数解析で用いられたドップラ周波数を用いて、ドップラ周波数ごとに相対速度を計算して出力する(ステップS3)。 Next, the position calculation section 22 calculates the distance between the radar device and the target object to be detected for each Doppler frequency based on the phase information output from the frequency analysis section 21, and also calculates the distance between the radar device and the target object to be detected. The azimuth angle of the target is calculated, and the distance information and angle information obtained as a result of the calculation are output (step S2). Furthermore, the velocity calculation unit 23 uses the Doppler frequency used in the frequency analysis in the frequency analysis unit 21 to calculate and output relative velocity for each Doppler frequency (step S3).

次に、物標位置検出部3の信号点抽出部31が、信号処理部2の周波数解析部21によって生成される周波数スペクトルにおけるピークに対応するドップラ周波数(別言すると、ピーク周波数)を検出するとともに、位置計算部22から出力されるドップラ周波数ごとの距離情報および角度情報を用いて、前記検出されたドップラ周波数(ピーク周波数)に対応するレーダ信号点の分布情報を生成して出力する(ステップS4)。また、クラスタ生成部32が、前記分布情報を用いてクラスタリング処理を行い、レーダ信号点のクラスタを生成する(ステップS5)。さらに、代表点特定部34が、クラスタNの代表点(例えば、重心位置)を特定する(ステップS6)。この際、追尾処理部33が、連続する時点においてクラスタ生成部32によって生成される、同一の物標に関するクラスタを関連づけることにより、クラスタNの追尾処理を行う。 Next, the signal point extraction unit 31 of the target position detection unit 3 detects the Doppler frequency (in other words, peak frequency) corresponding to the peak in the frequency spectrum generated by the frequency analysis unit 21 of the signal processing unit 2. At the same time, using the distance information and angle information for each Doppler frequency output from the position calculation unit 22, distribution information of radar signal points corresponding to the detected Doppler frequency (peak frequency) is generated and output (step S4). Further, the cluster generation unit 32 performs clustering processing using the distribution information to generate clusters of radar signal points (step S5). Further, the representative point specifying unit 34 specifies the representative point (for example, the position of the center of gravity) of cluster N (step S6). At this time, the tracking processing unit 33 performs tracking processing for cluster N by associating clusters related to the same target object generated by the cluster generation unit 32 at successive points in time.

次に、物標計算部4の関心領域設定部41が、クラスタNの代表点の位置を中心とする所定の範囲を関心領域として設定した上で、関心領域に関する4次元データ(X,Y,I,V)を出力する(ステップS7)。続いて、信号強度特定部42が、関心領域に関する4次元データ(X,Y,I,V)を速度成分ごとに参照し、速度成分別に、関心領域に関する信号強度情報を生成し(ステップS8)、また、VTデータ生成部43が、前記信号強度情報を用いて、各レーダスキャンの実行時刻に基づく経過時間T[秒]および相対速度V[km/h]と、これら経過時間T[秒]と相対速度V[km/h]との組み合わせに対応する信号強度情報(例えば、最大値Imax,平均値Iavg[dB])と、の組み合わせデータを生成する。 Next, the region of interest setting section 41 of the target object calculation section 4 sets a predetermined range centered on the position of the representative point of the cluster N as a region of interest, and then sets four-dimensional data (X, Y, I, V) are output (step S7). Next, the signal strength identifying unit 42 refers to the four-dimensional data (X, Y, I, V) regarding the region of interest for each velocity component, and generates signal strength information regarding the region of interest for each velocity component (step S8). In addition, the VT data generation unit 43 uses the signal strength information to calculate elapsed time T [seconds] and relative velocity V [km/h] based on the execution time of each radar scan, and these elapsed times T [seconds]. and signal strength information (for example, maximum value Imax, average value Iavg [dB]) corresponding to the combination of and relative speed V [km/h].

次に、特徴量計算部5の相関計算部51が、物標計算部4のVTデータ生成部43から出力される経過時間T[秒]と相対速度V[km/h]との組み合わせに対応する信号強度情報(例えば、最大値Imax,平均値Iavg[dB])のデータを用いて、所定の時間長さぶんのレーダスキャンごとの前記データの集合を時系列に配列して速度・強度の時系列データを作成し、作成した速度・強度の時系列データについて、レーダスキャンごとの相対速度Vの値と信号強度情報の値との組み合わせデータの時系列の配列を時間軸方向に1つずつシフトさせながら、シフト前の、時系列で配列された組み合わせデータの集合との対比における相関値を計算する(ステップS9)。 Next, the correlation calculation unit 51 of the feature value calculation unit 5 corresponds to the combination of the elapsed time T [seconds] and the relative velocity V [km/h] output from the VT data generation unit 43 of the target object calculation unit 4. Using data on signal strength information (for example, maximum value Imax, average value Iavg [dB]), the set of data for each radar scan of a predetermined length of time is arranged in time series to calculate speed and strength information. Create time series data, and for the created speed/intensity time series data, arrange the time series of combination data of relative velocity V value and signal strength information value for each radar scan one by one in the time axis direction. While shifting, a correlation value is calculated in comparison with a set of combination data arranged in time series before shifting (step S9).

続いて、周期計算部52が、相関計算部51から出力されるクラスタNの相関値の集合を時系列(即ち、速度・強度の時系列データについての時間軸方向のシフトの実行順)に配列したときの、相関値のピーク間の差分の平均値を計算する。また、振幅計算部53が、相関計算部51から出力されるクラスタNの相関値の集合を時系列(即ち、速度・強度の時系列データについての時間軸方向のシフトの実行順)に配列したときの、相関値の極大値および極小値の平均と分散とを計算する。さらに、相関変換部54が、相関計算部51から出力されるクラスタNの相関値の集合に対して高速フーリエ変換を行ってL1ノルムを計算する(ステップS10)。 Next, the period calculation unit 52 arranges the set of correlation values of cluster N output from the correlation calculation unit 51 in time series (that is, in the order in which shifts in the time axis direction are performed for the time series data of speed and intensity). Then, calculate the average value of the difference between the peaks of the correlation values. Further, the amplitude calculation unit 53 arranges the set of correlation values of cluster N output from the correlation calculation unit 51 in time series (that is, in the order in which shifts in the time axis direction are performed for the time series data of velocity and intensity). The average and variance of the maximum and minimum correlation values are calculated at the time. Further, the correlation transformation unit 54 performs fast Fourier transformation on the set of correlation values of cluster N output from the correlation calculation unit 51 to calculate the L 1 norm (step S10).

次に、識別部6が、特徴量計算部5から出力される、クラスタNごとの、相関値のピーク間の差分の平均値、相関値の極大値および極小値の平均と分散、ならびに高速フーリエ変換後のL1ノルムとしての周波数および振幅の値に基づいて、クラスタNごとに物標としての移動している人の属性を識別する(ステップS11)。なお、具体的には、相関値のピーク間の差分の平均値、相関値の極大値および極小値の平均と分散、ならびに高速フーリエ変換後のL1ノルムのうちの少なくとも一つに基づいて、物標としての移動している人の属性を識別するようにすればよい。 Next, the identification unit 6 calculates, for each cluster N, the average value of the difference between the peaks of the correlation values, the average and variance of the maximum and minimum values of the correlation values, and the fast Fourier Based on the frequency and amplitude values as L 1 norm after conversion, attributes of a moving person as a target are identified for each cluster N (step S11). Specifically, based on at least one of the average value of the difference between the peaks of the correlation values, the average and variance of the maximum and minimum values of the correlation values, and the L 1 norm after fast Fourier transformation, What is necessary is to identify the attributes of a moving person as a target.

ここで、人の属性の違いによる、レーダデータに表れる特性の比較の例を説明する。 Here, an example of comparison of characteristics appearing in radar data due to differences in human attributes will be explained.

図9は、物標が男性である場合の、50回のレーダスキャンによって取得されたレーダデータについてのクラスタNの相関値の集合と、物標が女性である場合の、50回のレーダスキャンによって取得されたレーダデータについてのクラスタNの相関値の集合とのそれぞれに関して、データのシフト回ごとの相関値の変動を折れ線で表して表示して1つのグラフ上に整理したグラフである。図9から、男性と女性とではデータのシフト回ごとの相関値の変動に明らかな違いがあることが看取され、データのシフト回ごとの相関値の変動に基づいて、男性であるか女性であるかを判別することが、延いては人の属性を識別することが可能であることが確認される。 Figure 9 shows a set of correlation values of cluster N for radar data acquired through 50 radar scans when the target is male, and a set of correlation values for cluster N for radar data acquired through 50 radar scans when the target is female. This is a graph in which fluctuations in correlation values for each shift of data are displayed as a polygonal line with respect to a set of correlation values of cluster N for acquired radar data and arranged on one graph. From Figure 9, it can be seen that there is a clear difference in the fluctuation of the correlation value for each data shift between men and women. It is confirmed that it is possible to identify a person's attributes by determining whether the person is a person.

図10は、物標が男性である場合の、50回のレーダスキャンによって取得されたレーダデータについてのクラスタNの相関値の集合と、物標が女性である場合の、50回のレーダスキャンによって取得されたレーダデータについてのクラスタNの相関値の集合とのそれぞれに関して(尚、レーダデータは図9に関連して用いたものと同じである)、高速フーリエ変換処理後のL1ノルムの値を折れ線で表示して1つのグラフ上に整理したグラフである。図10から、男性と女性とではL1ノルムの値に明らかな違いがあることが看取され、L1ノルムの値に基づいて男性であるか女性であるかを判別することが、延いては人の属性を識別することが可能であることが確認される。 Figure 10 shows a set of correlation values of cluster N for radar data acquired through 50 radar scans when the target is male, and a set of correlation values for cluster N for radar data acquired through 50 radar scans when the target is female. For each set of correlation values of cluster N for the acquired radar data (note that the radar data is the same as that used in connection with FIG. 9), the value of L 1 norm after fast Fourier transform processing This is a graph in which the graphs are displayed as a line and arranged on one graph. From Figure 10, it can be seen that there is a clear difference in the value of L 1 norm between men and women, and it is difficult to distinguish between men and women based on the value of L 1 norm. It is confirmed that it is possible to identify the attributes of a person.

図11は、相関値のピーク間の差分の平均値(図では「周期」と表記)、相関値の極大値および極小値の平均の値、ならびに相関値の極大値および極小値の分散の値についての直交3軸座標系に、レーダデータから検出・抽出されたクラスタについて計算される値に基づいてクラスタごとにプロットした結果である。図11に示すクラスタのそれぞれについては、男性についてのクラスタであるのか女性についてのクラスタであるのかが分かっている。図11に示す結果から、男性についてのクラスタと女性についてのクラスタとに分離するための平面が存在すると考えられる。 Figure 11 shows the average value of the difference between the peaks of correlation values (denoted as "period" in the figure), the average value of the local maximum value and local minimum value of the correlation value, and the value of the variance of the local maximum value and local minimum value of the correlation value. These are the results plotted for each cluster on an orthogonal three-axis coordinate system based on values calculated for clusters detected and extracted from radar data. For each cluster shown in FIG. 11, it is known whether it is a cluster about men or a cluster about women. From the results shown in FIG. 11, it is thought that there is a plane for separating clusters for men and clusters for women.

図12は、図11に示すデータについてサポートベクターマシンを用いて判別処理を行って超平面を計算した結果である。図12に示す結果から、相関値のピーク間の差分の平均値、相関値の極大値および極小値の平均の値、ならびに相関値の極大値および極小値の分散の値を適切に組み合わせて用いることにより、機械学習のアルゴリズムを用いて男性であるか女性であるかを判別することが、延いては人の属性を識別することが可能であることが確認される。 FIG. 12 shows the results of calculating a hyperplane by performing discrimination processing on the data shown in FIG. 11 using a support vector machine. From the results shown in Figure 12, the average value of the difference between the peaks of the correlation values, the average value of the local maximum value and local minimum value of the correlation value, and the value of the variance of the local maximum value and the local minimum value of the correlation value are used in appropriate combinations. This confirms that it is possible to use machine learning algorithms to determine whether a person is male or female, and by extension identify a person's attributes.

以上、この発明の実施の形態について説明したが、具体的な構成は、上記の実施の形態に限られるものではなく、この発明の要旨を逸脱しない範囲の設計の変更等があっても、この発明に含まれる。具体的には、この発明の要旨は、所定の時間長さぶんのレーダスキャンごとの相対速度Vの値と信号強度情報の値との組み合わせデータの集合を時系列に配列した上で、前記組み合わせデータの時系列の配列を時間軸方向に1つずつシフトさせながら、シフト前の、時系列で配列された前記組み合わせデータの集合との対比として計算される相関値を用いて物標の属性を判別することであり、この要旨を逸脱しない範囲の設計の変更等があっても、この発明に含まれる。例えば、この発明が適用される装置の具体的な構成は上記の実施の形態における属性の識別装置1に限定されるものではなく、具体的な装置構成は種々の態様があり得る。 Although the embodiments of this invention have been described above, the specific configuration is not limited to the above embodiments, and even if there are changes in the design within the scope of the gist of this invention, Included in invention. Specifically, the gist of the present invention is to arrange a set of combination data of relative velocity V values and signal strength information values for each radar scan for a predetermined length of time in chronological order, and then While shifting the time series arrangement of data one by one in the time axis direction, the attributes of the target are determined using the correlation value calculated as a comparison with the set of the combined data arranged in time series before the shift. However, even if there are changes in the design that do not depart from the scope of this invention, they are included in the present invention. For example, the specific configuration of the device to which the present invention is applied is not limited to the attribute identification device 1 in the above embodiment, and the specific device configuration may take various forms.

また、上記の実施の形態では特徴量計算部5がクラスタNの相関値の集合を時系列に配列したときの、相関値のピーク間の差分の平均値、相関値の極大値および極小値の平均、ならびに相関値の極大値および極小値の分散を計算するとともに識別部6が前記指標のうちの少なくとも一つを用いて人の属性の判別処理を行うようにしているが、相関計算部51によって計算されて出力されるクラスタNの相関値の集合を用いて計算される、物標としての移動している人の属性を識別するための特徴量としての指標は前記したものには限定されない。例えば、クラスタNの相関値の集合における相関値の極大値および極小値の中央値やレンジなどが計算されて判別処理に用いられるようにしてもよい。 Furthermore, in the above embodiment, when the feature calculation unit 5 arranges the set of correlation values of cluster N in time series, the average value of the difference between the peaks of the correlation values, the local maximum value, and the local minimum value of the correlation values are calculated. In addition to calculating the average and the variance of the maximum and minimum values of the correlation values, the identification unit 6 uses at least one of the indicators to perform a process of determining a person's attributes. The index as a feature amount for identifying the attributes of a moving person as a target, which is calculated using the set of correlation values of cluster N calculated and output by . For example, the median value, range, etc. of the maximum and minimum correlation values in the set of correlation values of cluster N may be calculated and used in the discrimination process.

また、代表点特定部34および関心領域設定部41が人の特定の部位に限定して関心領域を設定して前記関心領域に関して取得されたレーダデータが用いられるようにしてもよい。この点に関連して、図13は、物標を手先に限定した場合の、50回のレーダスキャンによって取得されたレーダデータについてのクラスタNの相関値の集合に関して、データのシフト回ごとの相関値の変動を折れ線で表して表示したグラフであり、(A)は成人についての結果であり、(B)は子供についての結果である。これら図13の(A)と(B)との比較から、成人と子供とではデータのシフト回ごとの相関値の変動に明らかな違いがあることが看取され、物標が手先に限定されたレーダデータを用いることにより、或いは、物標として手先を含むレーダデータを用いることにより、データのシフト回ごとの相関値の変動に基づいて成人であるか子供であるかを判別することが、延いては人の属性を識別することが可能であることが確認される。また、図13の(A)と(B)との整理から、相関値の極大値および極小値の分散が、成人の手先は1.1760となり、子供の手先は0.0773となる。また、相関値の極大値および極小値の平均が、成人の手先は0.9982となり、子供の手先は0.9991となる。これらの結果から、成人と子供とでは相関値の極大値および極小値の特に分散が大きく異なるので、物標が手先に限定されたレーダデータを用いることにより、或いは、物標として手先を含むレーダデータを用いることにより、データのシフト回ごとの相関値の極大値および極小値の分散に基づいて成人であるか子供であるかを判別することが、延いては人の属性を識別することが可能であることが確認される。 Alternatively, the representative point specifying section 34 and the region of interest setting section 41 may set a region of interest limited to a specific part of a person, and the radar data acquired regarding the region of interest may be used. In relation to this point, FIG. 13 shows the correlation for each shift of data regarding the set of correlation values of cluster N for radar data acquired by 50 radar scans when the target is limited to the hand. This is a graph showing the fluctuation of values as a line, in which (A) shows the results for adults, and (B) shows the results for children. From the comparison between (A) and (B) in Figure 13, it can be seen that there is a clear difference in the fluctuation of the correlation value for each data shift between adults and children, and the target object is limited to the hand. By using radar data that includes a hand as a target, it is possible to determine whether the person is an adult or a child based on the fluctuation of the correlation value each time the data is shifted. Furthermore, it is confirmed that it is possible to identify a person's attributes. Also, from the arrangement of FIGS. 13A and 13B, the variance of the maximum and minimum correlation values is 1.1760 for an adult's hand and 0.0773 for a child's hand. Further, the average of the maximum and minimum correlation values is 0.9982 for an adult's hand and 0.9991 for a child's hand. These results show that the variance of the maximum and minimum correlation values is significantly different between adults and children. By using data, it is possible to determine whether a person is an adult or a child based on the variance of the maximum and minimum correlation values for each data shift, and by extension identify the attributes of a person. It is confirmed that it is possible.

図14は、物標をももに限定した場合の、50回のレーダスキャンによって取得されたレーダデータについてのクラスタNの相関値の集合に関して、データのシフト回ごとの相関値の変動を折れ線で表して表示したグラフであり、(A)は成人についての結果であり、(B)は子供についての結果である。これら図14の(A)と(B)との比較から、成人と子供とではデータのシフト回ごとの相関値の変動に明らかな違いがあることが看取され、物標がももに限定されたレーダデータを用いることにより、或いは、物標としてももを含むレーダデータを用いることにより、データのシフト回ごとの相関値の変動に基づいて成人であるか子供であるかを判別することが、延いては人の属性を識別することが可能であることが確認される。また、図14の(A)と(B)との整理から、相関値の極大値および極小値の分散が、成人のももは1.8927となり、子供のももは0.2078となる。また、相関値の極大値および極小値の平均が、成人のももは0.9988となり、子供のももは0.9988となる。これらの結果から、成人と子供とでは相関値の極大値および極小値の特に分散が大きく異なるので、物標がももに限定されたレーダデータを用いることにより、或いは、物標としてももを含むレーダデータを用いることにより、データのシフト回ごとの相関値の極大値および極小値の分散に基づいて成人であるか子供であるかを判別することが、延いては人の属性を識別することが可能であることが確認される。 Figure 14 shows, as a polygonal line, the fluctuation of the correlation value for each data shift regarding the set of correlation values of cluster N for radar data acquired through 50 radar scans when the target object is limited to the thigh. These are graphs shown and displayed, in which (A) is the result for adults and (B) is the result for children. From the comparison between (A) and (B) in Figure 14, it can be seen that there is a clear difference in the fluctuation of the correlation value for each data shift between adults and children, and the target object is limited to the thigh. It is possible to determine whether a person is an adult or a child based on the fluctuation of the correlation value each time the data is shifted, by using the radar data that has been detected, or by using the radar data that includes the thigh as a target object. However, it has been confirmed that it is possible to identify a person's attributes. Furthermore, from the arrangement of FIGS. 14A and 14B, the variance of the maximum and minimum correlation values is 1.8927 for an adult's thigh and 0.2078 for a child's thigh. Further, the average of the maximum and minimum correlation values is 0.9988 for an adult's thigh and 0.9988 for a child's thigh. These results show that the variance of the maximum and minimum correlation values is significantly different between adults and children. By using the included radar data, it is possible to determine whether a person is an adult or a child based on the variance of the maximum and minimum correlation values for each data shift, and by extension identify the attributes of a person. It is confirmed that this is possible.

図15は、物標を手先に限定した場合の、50回のレーダスキャンによって取得されたレーダデータについてのクラスタNの相関値の集合に関して、高速フーリエ変換処理後のL1ノルム(具体的には、周波数と振幅との組み合わせデータ)を折れ線で表示したグラフであり、(A)は成人についての結果であり、(B)は子供についての結果である。これら図15の(A)と(B)との比較から、成人と子供とではL1ノルムの値に明らかな違いがあることが看取され、物標が手先に限定されたレーダデータを用いることにより、或いは、物標として手先を含むレーダデータを用いることにより、L1ノルムの値に基づいて成人であるか子供であるかを判別することが、延いては人の属性を識別することが可能であることが確認される。 FIG. 15 shows the L 1 norm (specifically, , combination data of frequency and amplitude) is displayed as a line, in which (A) shows the results for adults, and (B) shows the results for children. From the comparison between (A) and (B) in Figure 15, it can be seen that there is a clear difference in the value of L1 norm between adults and children, and using radar data where the target is limited to the hand. By doing so, or by using radar data that includes hands as targets, it is possible to determine whether a person is an adult or a child based on the value of the L1 norm, and by extension to identify the attributes of a person. It is confirmed that this is possible.

図16は、物標をももに限定した場合の、50回のレーダスキャンによって取得されたレーダデータについてのクラスタNの相関値の集合に関して、高速フーリエ変換処理後のL1ノルム(具体的には、周波数と振幅との組み合わせデータ)を折れ線で表示したグラフであり、(A)は成人についての結果であり、(B)は子供についての結果である。これら図16の(A)と(B)との比較から、成人と子供とではL1ノルムの値に明らかな違いがあることが看取され、物標がももに限定されたレーダデータを用いることにより、或いは、物標としてももを含むレーダデータを用いることにより、L1ノルムの値に基づいて成人であるか子供であるかを判別することが、延いては人の属性を識別することが可能であることが確認される。 FIG. 16 shows the L 1 norm (specifically, is a graph in which the combination data of frequency and amplitude is displayed as a line. (A) is the result for adults, and (B) is the result for children. From the comparison between (A) and (B) in Figure 16, it can be seen that there is a clear difference in the value of the L1 norm between adults and children, and it can be seen that there is a clear difference in the value of the L1 norm between adults and children. By using radar data that includes targets and thighs, it is possible to determine whether a person is an adult or a child based on the value of the L1 norm, and by extension to identify the attributes of a person. It is confirmed that it is possible to do so.

この発明は、簡便な機器構成によって歩行者の性別や年代などの属性の識別を高精度に行うことができるので、例えば、自動車の運転支援に関する分野や侵入者検知に関する分野において有用な技術として適用され得る。 This invention can identify attributes such as gender and age of pedestrians with high precision using a simple device configuration, so it can be applied as a useful technology in, for example, the field of automobile driving support and the field of intruder detection. can be done.

1 属性の識別装置
2 信号処理部
21 周波数解析部
22 位置計算部
23 速度計算部
3 物標位置検出部
31 信号点抽出部
32 クラスタ生成部
33 追尾処理部
34 代表点特定部
4 物標計算部
41 関心領域設定部
42 信号強度特定部
43 VTデータ生成部
5 特徴量計算部
51 相関計算部
52 周期計算部
53 振幅計算部
54 相関変換部
6 識別部
1 Attribute identification device 2 Signal processing unit 21 Frequency analysis unit 22 Position calculation unit 23 Velocity calculation unit 3 Target position detection unit 31 Signal point extraction unit 32 Cluster generation unit 33 Tracking processing unit 34 Representative point identification unit 4 Target calculation unit 41 Region of Interest Setting Unit 42 Signal Strength Specification Unit 43 VT Data Generation Unit 5 Feature Calculation Unit 51 Correlation Calculation Unit 52 Period Calculation Unit 53 Amplitude Calculation Unit 54 Correlation Transformation Unit 6 Identification Unit

Claims (4)

検知対象エリアへと電波を送信するとともに前記検知対象エリアからの電波を受信してレーダスキャンを行うレーダ装置から出力されるレーダデータを周波数解析して周波数スペクトルを出力する周波数解析部と、
前記周波数スペクトルに基づいて前記検知対象エリア内を移動する物標に関するレーダ信号点それぞれの相対速度を計算する速度計算部と、
前記周波数スペクトルのピークを検出して前記レーダ信号点を抽出する信号点抽出部と、
抽出された前記レーダ信号点のクラスタを生成するクラスタ生成部と、
前記クラスタの代表点を特定する代表点特定部と、
前記代表点を含む関心領域を設定する関心領域設定部と、
前記関心領域に含まれる前記レーダ信号点それぞれの前記相対速度の値別に信号強度情報を計算する信号強度特定部と、
所定の時間長さぶんの前記レーダスキャンごとの前記相対速度の値と前記信号強度情報の値との組み合わせデータの集合を時系列に配列した上で、前記組み合わせデータの時系列の配列を時間軸方向に1つずつシフトさせながら、シフト前の、時系列で配列された前記組み合わせデータの集合との対比における相関値を計算する相関計算部と、
前記相関値を用いて前記物標の属性を判別する識別部と、を有する、
ことを特徴とする属性の識別装置。
a frequency analysis unit that performs frequency analysis on radar data output from a radar device that transmits radio waves to a detection target area, receives radio waves from the detection target area, and performs radar scanning, and outputs a frequency spectrum;
a speed calculation unit that calculates the relative speed of each radar signal point regarding a target moving within the detection target area based on the frequency spectrum;
a signal point extraction unit that detects the peak of the frequency spectrum and extracts the radar signal point;
a cluster generation unit that generates a cluster of the extracted radar signal points;
a representative point identifying unit that identifies a representative point of the cluster;
a region of interest setting unit that sets a region of interest including the representative point;
a signal strength identification unit that calculates signal strength information for each of the relative velocity values of each of the radar signal points included in the region of interest;
After arranging a set of combination data of the relative velocity value and the signal strength information value for each radar scan for a predetermined time length in a time series, the time series arrangement of the combination data is arranged on a time axis. a correlation calculation unit that calculates a correlation value in comparison with the set of the combination data arranged in time series before shifting while shifting the data one by one in the direction;
an identification unit that determines an attribute of the target using the correlation value;
An attribute identification device characterized by:
前記識別部が、機械学習のアルゴリズムを用いて前記物標の属性を判別する、
ことを特徴とする請求項1に記載の属性の識別装置。
the identification unit determines attributes of the target using a machine learning algorithm;
The attribute identification device according to claim 1, characterized in that:
検知対象エリアへと電波を送信するとともに前記検知対象エリアからの電波を受信してレーダスキャンを行うレーダ装置から出力されるレーダデータを周波数解析して周波数スペクトルを出力する処理と、
前記周波数スペクトルに基づいて前記検知対象エリア内を移動する物標に関するレーダ信号点それぞれの相対速度を計算する処理と、
前記周波数スペクトルのピークを検出して前記レーダ信号点を抽出する処理と、
抽出された前記レーダ信号点のクラスタを生成する処理と、
前記クラスタの代表点を特定する処理と、
前記代表点を含む関心領域を設定する処理と、
前記関心領域に含まれる前記レーダ信号点それぞれの前記相対速度の値別に信号強度情報を計算する処理と、
所定の時間長さぶんの前記レーダスキャンごとの前記相対速度の値と前記信号強度情報の値との組み合わせデータの集合を時系列に配列した上で、前記組み合わせデータの時系列の配列を時間軸方向に1つずつシフトさせながら、シフト前の、時系列で配列された前記組み合わせデータの集合との対比における相関値を計算する処理と、
前記相関値を用いて前記物標の属性を判別する処理と、を有する、
ことを特徴とする属性の識別方法。
A process of frequency-analyzing radar data output from a radar device that transmits radio waves to a detection target area, receives radio waves from the detection target area, and performs a radar scan, and outputs a frequency spectrum;
a process of calculating the relative speed of each radar signal point regarding a target moving within the detection target area based on the frequency spectrum;
a process of detecting the peak of the frequency spectrum and extracting the radar signal point;
a process of generating a cluster of the extracted radar signal points;
a process of identifying representative points of the cluster;
a process of setting a region of interest including the representative point;
a process of calculating signal strength information for each of the relative velocity values of each of the radar signal points included in the region of interest;
After arranging a set of combination data of the relative velocity value and the signal strength information value for each radar scan for a predetermined time length in a time series, the time series arrangement of the combination data is arranged on a time axis. a process of calculating a correlation value in comparison with the set of the combination data arranged in time series before the shift while shifting the data one by one in the direction;
a process of determining an attribute of the target using the correlation value;
A method for identifying attributes characterized by:
機械学習のアルゴリズムを用いて前記物標の属性を判別する、
ことを特徴とする請求項3に記載の属性の識別方法。
determining attributes of the target using a machine learning algorithm;
The attribute identification method according to claim 3, characterized in that:
JP2020083417A 2020-05-11 2020-05-11 Attribute identification device and identification method Active JP7377494B2 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2020083417A JP7377494B2 (en) 2020-05-11 2020-05-11 Attribute identification device and identification method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2020083417A JP7377494B2 (en) 2020-05-11 2020-05-11 Attribute identification device and identification method

Publications (2)

Publication Number Publication Date
JP2021179323A JP2021179323A (en) 2021-11-18
JP7377494B2 true JP7377494B2 (en) 2023-11-10

Family

ID=78511240

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2020083417A Active JP7377494B2 (en) 2020-05-11 2020-05-11 Attribute identification device and identification method

Country Status (1)

Country Link
JP (1) JP7377494B2 (en)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR102598250B1 (en) * 2023-07-11 2023-11-06 메타빌드주식회사 Target detection method using radar signal, target detection system and computer program for the same

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2010250501A (en) 2009-04-14 2010-11-04 Hitachi Automotive Systems Ltd External environment recognition device for vehicle and vehicle system using the same
US9661470B1 (en) 2015-10-05 2017-05-23 Google Inc. Methods and systems for locating an actor within an environment
JP2019128634A (en) 2018-01-22 2019-08-01 株式会社日立製作所 Optimization apparatus and optimization method for hyper parameter
JP2019200083A (en) 2018-05-15 2019-11-21 株式会社デンソーテン Radar device and signal processing method
US20190391250A1 (en) 2018-06-26 2019-12-26 Zoox, Inc. Radar clustering and velocity disambiguation

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2010250501A (en) 2009-04-14 2010-11-04 Hitachi Automotive Systems Ltd External environment recognition device for vehicle and vehicle system using the same
US9661470B1 (en) 2015-10-05 2017-05-23 Google Inc. Methods and systems for locating an actor within an environment
JP2019128634A (en) 2018-01-22 2019-08-01 株式会社日立製作所 Optimization apparatus and optimization method for hyper parameter
JP2019200083A (en) 2018-05-15 2019-11-21 株式会社デンソーテン Radar device and signal processing method
US20190391250A1 (en) 2018-06-26 2019-12-26 Zoox, Inc. Radar clustering and velocity disambiguation

Also Published As

Publication number Publication date
JP2021179323A (en) 2021-11-18

Similar Documents

Publication Publication Date Title
Jin et al. Multiple patients behavior detection in real-time using mmWave radar and deep CNNs
Mobasseri et al. A time-frequency classifier for human gait recognition
Erol et al. Fall motion detection using combined range and Doppler features
Kim et al. Radar-based human activity recognition combining range–time–Doppler maps and range-distributed-convolutional neural networks
Wang et al. Millimeter-wave frequency modulated continuous wave radar-based soft fall detection using pattern contour-confined Doppler-time maps
Qi et al. Quantitative evaluation of channel micro-Doppler capacity for MIMO UWB radar human activity signals based on time–frequency signatures
JP2019152441A (en) Vital sensor
US20220260702A1 (en) Radar-Based Target Tracker
Lim et al. DNN-based human face classification using 61 GHz FMCW radar sensor
Yang et al. Spatial-temporal-circulated GLCM and physiological features for in-vehicle people sensing based on IR-UWB radar
Kittiyanpunya et al. Millimeter-Wave Radar-Based Elderly Fall Detection Fed by One-Dimensional Point Cloud and Doppler
Shankar et al. Radar-based non-intrusive fall motion recognition using deformable convolutional neural network
JP7377494B2 (en) Attribute identification device and identification method
Yang et al. Fall feature enhancement and fusion using the stockwell transform with dual mmWave radars
Hu et al. Radar-Based Fall Detection: A Survey
JP2024514618A (en) DNN-based human face classification
Alkasimi et al. Geolocation tracking for human identification and activity recognition using radar deep transfer learning
Sona et al. A measurement approach based on micro-Doppler maps for human motion analysis and detection
Raeis et al. InARMS: Individual activity recognition of multiple subjects with FMCW radar
Darlis et al. Autonomous Human and Animal Classification Using Synthetic 2D Tensor Data Based on Dual-Receiver mmWave Radar System
CN114767074B (en) Vital sign measuring method, equipment and storage medium
Zheng et al. Real-Time Fall Recognition Using a Lightweight Convolution Neural Network Based on Millimeter-Wave Radar
Fusco et al. RadarSleepNet: Sleep Pose Classification via PointNet++ and 5D Radar Point Clouds
Nocera et al. Walking pattern identification of FMCW radar data based on a combined CNN and bi-LSTM approach
Zhong et al. Human motion recognition in small sample scenarios based on GaN and CNN models

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20230330

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

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20231020

R150 Certificate of patent or registration of utility model

Ref document number: 7377494

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150