特許 7377494 属性の識別装置および識別方法

(19)【発行国】日本国特許庁(JP)

(12)【公報種別】特許公報(B2)

(11)【特許番号】

(24)【登録日】2023-11-01

(45)【発行日】2023-11-10

(54)【発明の名称】属性の識別装置および識別方法

(51)【国際特許分類】

   G01S 7/41 20060101AFI20231102BHJP

   G01S 13/34 20060101ALI20231102BHJP

   G01S 13/58 20060101ALI20231102BHJP

   G08G 1/01 20060101ALI20231102BHJP

【ＦＩ】

G01S7/41

G01S13/34

G01S13/58 200

G08G1/01 F

【請求項の数】 4

(21)【出願番号】P 2020083417

(22)【出願日】2020-05-11

(65)【公開番号】P2021179323

(43)【公開日】2021-11-18

【審査請求日】2023-03-30

(73)【特許権者】

【識別番号】318006365

【氏名又は名称】ＪＲＣモビリティ株式会社

(73)【特許権者】

【識別番号】506060258

【氏名又は名称】公立大学法人北九州市立大学

(74)【代理人】

【識別番号】100126561

【弁理士】

【氏名又は名称】原嶋 成時郎

(74)【代理人】

【識別番号】100141678

【弁理士】

【氏名又は名称】佐藤 和彦

(72)【発明者】

【氏名】星 将広

(72)【発明者】

【氏名】矢野 邦哲

(72)【発明者】

【氏名】時枝 幸伸

(72)【発明者】

【氏名】京地 清介

(72)【発明者】

【氏名】千々和 憂希

【審査官】藤田 都志行

(56)【参考文献】

【文献】特開２０１９－２０００８３（ＪＰ，Ａ）

【文献】特開２０１０－２５０５０１（ＪＰ，Ａ）

【文献】特開２０１９－１２８６３４（ＪＰ，Ａ）

【文献】米国特許出願公開第２０１９／０３９１２５０（ＵＳ，Ａ１）

【文献】米国特許第９６６１４７０（ＵＳ，Ｂ１）

(58)【調査した分野】(Int.Cl.，ＤＢ名)

Ｇ０１Ｓ ７／００－ ７／５１

Ｇ０１Ｓ １３／００－１３／９５

Ｇ０１Ｓ １７／００－１７／９５

Ｇ０８Ｇ １／０１

Ｇ０６Ｖ １０／２０－１０／３６

Ｇ０６Ｖ ２０／５０－２０／５９

ＪＳＴＰｌｕｓ／ＪＳＴ７５８０（ＪＤｒｅａｍＩＩＩ）

(57)【特許請求の範囲】

【請求項1】

検知対象エリアへと電波を送信するとともに前記検知対象エリアからの電波を受信してレーダスキャンを行うレーダ装置から出力されるレーダデータを周波数解析して周波数スペクトルを出力する周波数解析部と、
前記周波数スペクトルに基づいて前記検知対象エリア内を移動する物標に関するレーダ信号点それぞれの相対速度を計算する速度計算部と、
前記周波数スペクトルのピークを検出して前記レーダ信号点を抽出する信号点抽出部と、
抽出された前記レーダ信号点のクラスタを生成するクラスタ生成部と、
前記クラスタの代表点を特定する代表点特定部と、
前記代表点を含む関心領域を設定する関心領域設定部と、
前記関心領域に含まれる前記レーダ信号点それぞれの前記相対速度の値別に信号強度情報を計算する信号強度特定部と、
所定の時間長さぶんの前記レーダスキャンごとの前記相対速度の値と前記信号強度情報の値との組み合わせデータの集合を時系列に配列した上で、前記組み合わせデータの時系列の配列を時間軸方向に１つずつシフトさせながら、シフト前の、時系列で配列された前記組み合わせデータの集合との対比における相関値を計算する相関計算部と、
前記相関値を用いて前記物標の属性を判別する識別部と、を有する、
ことを特徴とする属性の識別装置。

【請求項2】

前記識別部が、機械学習のアルゴリズムを用いて前記物標の属性を判別する、
ことを特徴とする請求項１に記載の属性の識別装置。

【請求項3】

検知対象エリアへと電波を送信するとともに前記検知対象エリアからの電波を受信してレーダスキャンを行うレーダ装置から出力されるレーダデータを周波数解析して周波数スペクトルを出力する処理と、
前記周波数スペクトルに基づいて前記検知対象エリア内を移動する物標に関するレーダ信号点それぞれの相対速度を計算する処理と、
前記周波数スペクトルのピークを検出して前記レーダ信号点を抽出する処理と、
抽出された前記レーダ信号点のクラスタを生成する処理と、
前記クラスタの代表点を特定する処理と、
前記代表点を含む関心領域を設定する処理と、
前記関心領域に含まれる前記レーダ信号点それぞれの前記相対速度の値別に信号強度情報を計算する処理と、
所定の時間長さぶんの前記レーダスキャンごとの前記相対速度の値と前記信号強度情報の値との組み合わせデータの集合を時系列に配列した上で、前記組み合わせデータの時系列の配列を時間軸方向に１つずつシフトさせながら、シフト前の、時系列で配列された前記組み合わせデータの集合との対比における相関値を計算する処理と、
前記相関値を用いて前記物標の属性を判別する処理と、を有する、
ことを特徴とする属性の識別方法。

【請求項4】

機械学習のアルゴリズムを用いて前記物標の属性を判別する、
ことを特徴とする請求項３に記載の属性の識別方法。

【発明の詳細な説明】

【技術分野】

【0001】

この発明は、属性の識別装置および識別方法に関し、特に、電波を用いて物標としての歩行者の属性を識別する技術に関する。

【背景技術】

【0002】

歩行者の属性を識別する技術として、例えば、２５台のカメラからなる多視点同期歩容撮影システムを構築して、広い年齢層の被験者の歩容を撮影することで歩容データベースを作成し、性別と年齢とのそれぞれに関して識別実験を行った結果に基づいて性別および年齢の識別に適した分類として子供、成人男性、成人女性、および高齢者という４つのクラスを設定し、前記４つのクラスへの識別を行う、技術が知られている（非特許文献１参照）。

【先行技術文献】

【非特許文献】

【0003】

【文献】万波秀年ほか「歩容における性別・年齢の分類と特徴解析」，電子情報通信学会論文誌 Ｄ，社団法人電子情報通信学会，２００９年，Ｖｏｌ．Ｊ９２－Ｄ，Ｎｏ．８，ｐｐ．１３７３－１３８２

【発明の概要】

【発明が解決しようとする課題】

【0004】

しかしながら、個人の属性を識別するために多数のカメラが必要とされるような大掛かりなシステムは、実用的とは言い難い、という問題がある。また、実環境での使用を考慮した場合、画像では特に人同士の重なりなどの背景変化に弱く、部位ごとの周波数特徴を効果的に抽出することは困難である、という問題がある。

【0005】

そこで本発明は、簡便な機器構成によって歩行者の属性の識別を高精度に行うことが可能な属性の識別装置および識別方法を提供することを目的とする。

【課題を解決するための手段】

【0006】

上記課題を解決するために、請求項１に記載の発明は、検知対象エリアへと電波を送信するとともに前記検知対象エリアからの電波を受信してレーダスキャンを行うレーダ装置から出力されるレーダデータを周波数解析して周波数スペクトルを出力する周波数解析部と、前記周波数スペクトルに基づいて前記検知対象エリア内を移動する物標に関するレーダ信号点それぞれの相対速度を計算する速度計算部と、前記周波数スペクトルのピークを検出して前記レーダ信号点を抽出する信号点抽出部と、抽出された前記レーダ信号点のクラスタを生成するクラスタ生成部と、前記クラスタの代表点を特定する代表点特定部と、前記代表点を含む関心領域を設定する関心領域設定部と、前記関心領域に含まれる前記レーダ信号点それぞれの前記相対速度の値別に信号強度情報を計算する信号強度特定部と、所定の時間長さぶんの前記レーダスキャンごとの前記相対速度の値と前記信号強度情報の値との組み合わせデータの集合を時系列に配列した上で、前記組み合わせデータの時系列の配列を時間軸方向に１つずつシフトさせながら、シフト前の、時系列で配列された前記組み合わせデータの集合との対比における相関値を計算する相関計算部と、前記相関値を用いて前記物標の属性を判別する識別部と、を有する、ことを特徴とする属性の識別装置である。

【0007】

請求項２に記載の発明は、請求項１に記載の属性の識別装置において、前記識別部が、機械学習のアルゴリズムを用いて前記物標の属性を判別する、ことを特徴とする。

【0008】

また、請求項３に記載の発明は、検知対象エリアへと電波を送信するとともに前記検知対象エリアからの電波を受信してレーダスキャンを行うレーダ装置から出力されるレーダデータを周波数解析して周波数スペクトルを出力する処理と、前記周波数スペクトルに基づいて前記検知対象エリア内を移動する物標に関するレーダ信号点それぞれの相対速度を計算する処理と、前記周波数スペクトルのピークを検出して前記レーダ信号点を抽出する処理と、抽出された前記レーダ信号点のクラスタを生成する処理と、前記クラスタの代表点を特定する処理と、前記代表点を含む関心領域を設定する処理と、前記関心領域に含まれる前記レーダ信号点それぞれの前記相対速度の値別に信号強度情報を計算する処理と、所定の時間長さぶんの前記レーダスキャンごとの前記相対速度の値と前記信号強度情報の値との組み合わせデータの集合を時系列に配列した上で、前記組み合わせデータの時系列の配列を時間軸方向に１つずつシフトさせながら、シフト前の、時系列で配列された前記組み合わせデータの集合との対比における相関値を計算する処理と、前記相関値を用いて前記物標の属性を判別する処理と、を有する、ことを特徴とする属性の識別方法である。

【0009】

請求項４に記載の発明は、請求項３に記載の属性の識別方法において、機械学習のアルゴリズムを用いて前記物標の属性を判別する、ことを特徴とする。

【発明の効果】

【0010】

請求項１に記載の発明や請求項３に記載の発明によれば、相対速度の値と信号強度情報の値との組み合わせデータの集合の時系列の配列を時間軸方向に１つずつシフトさせながら、シフト前の、時系列で配列された前記組み合わせデータの集合との対比として計算される相関値を用いて物標の属性を判別するようにしているので、物標の属性に応じた特有の性質を的確に把握して物標の属性の識別を高精度に行うことが可能となる。

【0011】

また、レーダ（電波）であれば、画像処理のような複雑な前処理を施す必要がなく、物標それぞれの速度成分を定量的に抽出することができるので、周波数領域の特徴を簡易かつ高精度に求めることが可能となる。

【0012】

請求項２に記載の発明や請求項４に記載の発明によれば、機械学習のアルゴリズムを用いるようにしているので、相対速度の値と信号強度情報の値との組み合わせデータの集合の時系列の配列に関する、時間軸方向における相関値に基づいて物標の属性の識別を適切に行うことが可能となる。

【図面の簡単な説明】

【0013】

【図1】この発明の実施の形態に係る属性の識別装置の概略構成を示す機能ブロック図である。

【図2】図１の属性の識別装置における処理手順であるとともにこの発明の実施の形態に係る属性の識別方法の処理手順を示すフローチャートである。

【図3】速度・強度の時系列データの一例を示すグラフである。

【図4】相関計算部における相関値の計算の仕方を説明する図である。

【図5】クラスタについての相関値の集合のデータの一例を示すグラフである。

【図6】クラスタの相関値の集合における相関値のピーク間の差分を説明する図である。

【図7】クラスタの相関値の集合における相関値の極大値および極小値を説明する図である。

【図8】クラスタの相関値の集合の、高速フーリエ変換処理後のＬ¹ノルムの例を示す図である。

【図9】クラスタの相関値の集合に関する、データのシフト回ごとの相関値の変動について、男性の場合と女性の場合との比較の例を示す図である。

【図10】クラスタの相関値の集合に関する、高速フーリエ変換処理後のＬ¹ノルムについて、男性の場合と女性の場合との比較の例を示す図である。

【図11】相関値に纏わる３つの指標についての直交３軸座標系に、クラスタ毎にプロットした例を示す図である。

【図12】図１１に示すデータについてサポートベクターマシンを用いて判別処理を行って超平面を計算した結果を示す図である。

【図13】物標を手先に限定した場合の、クラスタの相関値の集合に関する、データのシフト回ごとの相関値の変動の例を示す図である。（Ａ）は成人の場合の図である。（Ｂ）は子供の場合の図である。

【図14】物標をももに限定した場合の、クラスタの相関値の集合に関する、データのシフト回ごとの相関値の変動の例を示す図である。（Ａ）は成人の場合の図である。（Ｂ）は子供の場合の図である。

【図15】物標を手先に限定した場合の、クラスタの相関値の集合に関する、高速フーリエ変換処理後のＬ¹ノルムの例を示す図である。（Ａ）は成人の場合の図である。（Ｂ）は子供の場合の図である。

【図16】物標をももに限定した場合の、クラスタの相関値の集合に関する、高速フーリエ変換処理後のＬ¹ノルムの例を示す図である。（Ａ）は成人の場合の図である。（Ｂ）は子供の場合の図である。

【発明を実施するための形態】

【0014】

以下、この発明を図示の実施の形態に基づいて説明する。

【0015】

図１は、この発明の実施の形態に係る属性の識別装置１の概略構成を示す機能ブロック図である。図２は、属性の識別装置１における処理手順であるとともにこの発明の実施の形態に係る属性の識別方法の処理手順を示すフローチャートである。

【0016】

この属性の識別装置１は、検知対象エリア内の移動物体を物標として検知して前記物標の属性を識別する装置であり、主として、信号処理部２と、物標位置検出部３と、物標計算部４と、特徴量計算部５と、識別部６と、を備える。属性の識別装置１は、具体的には、検知対象エリア内を移動している人の属性を識別する。

【0017】

この実施の形態に係る属性の識別装置１は、検知対象エリアへと電波を送信するとともに前記検知対象エリアからの電波を受信してレーダスキャンを行うレーダ装置から出力されるレーダデータを周波数解析して周波数スペクトルを出力する周波数解析部２１と、周波数スペクトルに基づいて検知対象エリア内を移動する物標に関するレーダ信号点それぞれの相対速度を計算する速度計算部２３と、周波数スペクトルのピークを検出してレーダ信号点を抽出する信号点抽出部３１と、抽出されたレーダ信号点のクラスタを生成するクラスタ生成部３２と、クラスタの代表点を特定する代表点特定部３４と、代表点を含む関心領域を設定する関心領域設定部４１と、関心領域に含まれる前記レーダ信号点それぞれの前記相対速度の値別に信号強度情報を計算する信号強度特定部４２と、所定の時間長さぶんのレーダスキャンごとの相対速度の値と信号強度情報の値との組み合わせデータの集合を時系列に配列した上で、前記組み合わせデータの時系列の配列を時間軸方向に１つずつシフトさせながら、シフト前の、時系列で配列された前記組み合わせデータの集合との対比における相関値を計算する相関計算部５１と、相関値を用いて物標の属性を判別する識別部６と、を有する、ようにしている（図１参照）。

【0018】

また、この実施の形態に係る属性の識別方法は、検知対象エリアへと電波を送信するとともに前記検知対象エリアからの電波を受信してレーダスキャンを行うレーダ装置から出力されるレーダデータを周波数解析して周波数スペクトルを出力する処理（ステップＳ１）と、周波数スペクトルに基づいて検知対象エリア内を移動する物標に関するレーダ信号点それぞれの相対速度を計算する処理（ステップＳ３）と、周波数スペクトルのピークを検出してレーダ信号点を抽出する処理（ステップＳ４）と、抽出されたレーダ信号点のクラスタを生成する処理（ステップＳ５）と、クラスタの代表点を特定する処理（ステップＳ６）と、代表点を含む関心領域を設定する処理（ステップＳ７）と、関心領域に含まれる前記レーダ信号点それぞれの前記相対速度の値別に信号強度情報を計算する処理（ステップＳ８）と、所定の時間長さぶんのレーダスキャンごとの相対速度の値と信号強度情報の値との組み合わせデータの集合を時系列に配列した上で、前記組み合わせデータの時系列の配列を時間軸方向に１つずつシフトさせながら、シフト前の、時系列で配列された前記組み合わせデータの集合との対比における相関値を計算する処理（ステップＳ９）と、相関値を用いて物標の属性を判別する処理（ステップＳ１０、Ｓ１１）と、を有する、ようにしている（図２参照）。

【0019】

この発明では、例えば、レーダ方式としてＦＭＣＷ（Ｆrequency Ｍodulated－Ｃontinuous Ｗave の略；周波数変調連続波）方式のレーダスキャンを行うレーダ装置によって取得されて出力されるレーダデータが用いられ得る。ＦＭＣＷ方式では、周波数変調した連続波を送信するとともに物標の表面で反射された反射波を受信し、送信信号と受信信号との間での周波数差に基づいて物標の位置や相対速度を検出する。なお、ＦＭＣＷ方式では、周波数が三角波状または鋸波状に変化する電波（送信信号）が送信される。

【0020】

この発明で用いられ得るレーダデータを出力するレーダ装置としては、例えば、送信部と受信部とを備えて電波を送信するとともに受信する機能を備え、ＦＭＣＷ方式のレーダスキャンを行ってレーダデータを出力する装置が挙げられる。

【0021】

レーダ装置の送信部としては、例えば、所定電圧を生成して出力する電圧発生器、電圧発生器から出力される前記所定電圧に応じた周波数を有する電波（送信信号）を生成して出力する電圧制御発振器、および、電圧制御発振器から出力される電波（送信信号）を検知対象エリアへと送信波として送信する（別言すると、出射する、放射する）送信アンテナなどを含む仕組みが挙げられる。

【0022】

送信部は、例えば７９ＧＨｚ帯、６０ＧＨｚ帯、或いは７６ＧＨｚ帯の周波数を有する電波（尚、「ミリ波」とも呼ばれる）を生成して送信アンテナを介して送信する。なお、この発明では、高周波数帯の電波を利用することが好ましい。

【0023】

レーダ装置の受信部としては、例えば、送信アンテナから送信された電波（即ち、送信波）が検知対象エリア内を移動する物標の表面で反射された電波（即ち、反射波；「ドップラ反射波」とも呼ばれる）を含む電波を受信波として受信する受信アンテナ、送信部から供給される電波（送信信号）と受信アンテナから出力される電波（受信信号）とをミキシングして差分信号を生成して出力するミキサ、および、ミキサから出力される差分信号に対して所定のサンプリング周波数を用いてサンプリング処理（別言すると、アナログ－デジタル変換処理）を施して前記差分信号をデジタルデータに変換してデジタル信号を出力するＡ／Ｄ変換器などを含む仕組みが挙げられる。

【0024】

なお、ＦＭＣＷ方式では、単数もしくは複数のアンテナによってアレーアンテナが形成され、受信アンテナとしてのアレーアンテナによって電波（反射波）を受信する。そして、ミキサおよびＡ／Ｄ変換器は、受信アンテナごとに設けられる。

【0025】

レーダ装置から出力されるレーダデータとしての差分信号は、送信部から供給される電波（送信信号）の周波数成分と受信アンテナから出力される電波（受信信号）の周波数成分との差の周波数成分を有する信号（つまり、ビート周波数を有する信号であり、「ビート信号」とも呼ばれる）である。

【0026】

レーダ装置は、例えば、あくまで一例として挙げると、自動車などの車両に取り付けられて車両周囲の所定範囲を検知対象エリアとしたり、建物などの構造物に取り付けられて構造体周囲や建物敷地内の所定範囲を検知対象エリアとしたりする。

【0027】

信号処理部２は、レーダ装置から出力されるレーダデータ（即ち、レーダ生データ、レーダＲａｗデータ）の周波数解析を行うとともに物標の位置および相対速度を計算する。信号処理部２は、周波数解析部２１、位置計算部２２、および速度計算部２３を含んで構成され、レーダ装置から出力されるレーダデータとしての差分信号に対して周波数解析を施して、検知対象エリア内を移動する物標を検知して前記物標の位置および相対速度を計算する。

【0028】

属性の識別装置１は、レーダ装置から出力されるレーダデータを記録して格納しておく記憶部を備えて前記記憶部からレーダデータを読み込んで用いるようにしてもよく、また、レーダ装置から出力されてサーバなどの外部記憶装置に格納されているレーダデータを読み込んで用いるようにしてもよい。

【0029】

周波数解析部２１は、レーダ装置（具体的には、受信部のＡ／Ｄ変換器）から出力されるレーダデータ（別言すると、波形データ）としての差分信号（別言すると、ビート信号；尚、デジタル信号である）の周波数解析を行う。周波数解析部２１は、具体的には、１回のレーダスキャンあたりのサンプリング時間幅（例えば、５０ミリ秒間）ぶんの波形データ、具体的には差分信号（ビート信号）の振幅（即ち、受信レベル）に対して高速フーリエ変換処理を施して前記差分信号の振幅の周波数分布を示す周波数スペクトルを生成し、スペクトル強度および位相情報を出力する。周波数スペクトルは、差分信号に含まれる各周波数成分の振幅（受信レベル）を示す。

【0030】

高速フーリエ変換処理に際しては、検知対象とする物標の例えば移動速度が考慮されるなどした上で、前記物標の移動によって発生し得る周波数の範囲が設定されるとともに周波数間隔がドップラ周波数として設定される。周波数解析部２１は、つまり、レーダ装置から出力されるレーダデータ（波形データ）についてドップラ成分、すなわち差分信号のドップラ周波数ごとのスペクトル強度と位相とを周波数解析して求める。

【0031】

位置計算部２２は、周波数解析部２１から出力される周波数スペクトルにおける位相情報に基づいて、レーダ装置（具体的には、レーダ装置の受信アンテナ）と検知対象の物標との間の距離に関する距離情報、および、レーダ装置（具体的には、レーダ装置の受信アンテナ）に対する前記物標の方位角度に関する角度情報を生成して出力する。距離情報および角度情報は、ドップラ周波数ごとに求められる。

【0032】

差分信号の位相情報に基づく距離や方位角度の計算の仕法は、公知の手法が存在し、また、この発明では特定の手法には限定されないので、ここでは詳細の説明は省略する。例えば、レーダ装置と物標との間の距離は、送信波の周波数と受信波の周波数との差がレーダ装置と物標との間の距離に比例して増減する、ことを利用する手法などの公知の手法によって計算され得る。また、レーダ装置に対する物標の方位角度は、具体的には例えば、ＭＵＳＩＣ（ＭＵltiple ＳＩgnal Ｃlassification の略）法やＥＳＰＲＩＴ（Ｅstimation of Ｓignal Ｐarameters via Ｒotational Ｉnvariance Ｔechniques の略）法などの公知の測角方式によって計算され得る。

【0033】

速度計算部２３は、周波数解析部２１における周波数解析で用いられたドップラ周波数のそれぞれを用いて、レーダ装置（具体的には、レーダ装置の受信アンテナ）に対する物標の速度である相対速度（具体的には、送信波が物標の表面で反射した時の瞬時の相対速度）を計算して出力する。すなわち、相対速度は、ドップラ周波数ごとに求められる。

【0034】

相対速度の計算の仕法は、公知の手法が存在し、また、この発明では特定の手法には限定されないので、ここでは詳細の説明は省略する。例えば、相対速度は、レーダに対して相対的に移動している物標の表面で反射して受信アンテナによって受信される電波（即ち、受信波）の周波数は、送信アンテナから送信される電波（即ち、送信波）の周波数に対して、前記物標の表面で反射した際に受信波が前記物標の速度による影響を受け、ドップラ効果により、前記物標とレーダ装置との間の相対速度に応じてシフトする、ことを利用する手法などの公知の手法によって計算され得る。

【0035】

信号処理部２は、上記の処理により、レーダ装置から出力されるレーダデータ（即ち、レーダ生データ、レーダＲａｗデータ）に対して信号処理／周波数解析を施して、速度成分ごとの、直交座標系における強度分布を示した平面図であるＰＰＩ（Ｐlan Ｐosition Ｉndicator の略；平面図表示）を生成し、４次元データ（Ｘ，Ｙ，Ｉ，Ｖ）を出力する。なお、ＸおよびＹは、ｘ軸とｙ軸とが相互に直交する２次元直交座標系におけるｘ座標およびｙ座標（単位：ｍ）である。Ｉは信号強度（受信レベル）（単位：ｄＢ）であり、Ｖは相対速度（単位：ｋｍ／ｈ）である。信号処理部２から出力される４次元データ（Ｘ，Ｙ，Ｉ，Ｖ）の相対速度Ｖ［ｋｍ／ｈ］は、必要に応じ、相対速度Ｖの値が例えば０．１［ｋｍ／ｈ］などの所定のピッチで区分された上で相対速度Ｖの実際の値が前記区分に当てはめられて、相対速度のランクとして出力されるようにしてもよい。

【0036】

物標位置検出部３は、信号処理部２における処理結果を用いて、物標としての移動している人を検出するとともに前記移動している人に関する情報を計算する。物標位置検出部３は、信号点抽出部３１、クラスタ生成部３２、追尾処理部３３、および代表点特定部３４を含んで構成され、信号処理部２から出力される、信号処理部２における処理結果である４次元データ（Ｘ，Ｙ，Ｉ，Ｖ）を用いて、検知対象エリア内を移動している人を検出して前記移動している人に関する情報を計算する。

【0037】

信号点抽出部３１は、信号処理部２の周波数解析部２１によって生成される周波数スペクトルを用いて、所定の大きさの振幅を有するピークに対応するドップラ周波数（別言すると、ピーク周波数）を検出して、検知対象エリア内を移動する物標に関するレーダ信号点の周波数として抽出する。具体的には、信号点抽出部３１は、ピークの検出のための閾値を設定し、周波数スペクトルの振幅について、前記閾値以上の振幅をピークと判定して抽出する。

【0038】

ピークの検出・抽出の仕法は、特定の方法・方式に限定されるものではなく、クラッタと物標とを区別して物標からの反射成分と考えられる周波数を検出・抽出することができればどのような方法・方式であってもよい。ピークの検出・抽出の仕法として、具体的には例えばＣＦＡＲ（Ｃonstant Ｆalse Ａlarm Ｒate の略）アルゴリズムが用いられるようにしてもよい。

【0039】

信号点抽出部３１は、また、位置計算部２２から出力されるドップラ周波数ごとの距離情報および角度情報を用いて、空間位置を座標軸とする座標系における、上記において検出されたドップラ周波数（ピーク周波数）に対応するレーダ信号点の分布に関する分布情報を生成して出力する。座標系としては、例えば、レーダ装置の設置位置を基準とした、動径距離と動径方向とについての極座標系、または、左右方向における距離と奥行方向における距離とについての直交座標系が用いられる。

【0040】

クラスタ生成部３２は、信号点抽出部３１から出力されるレーダ信号点の分布情報を用いてクラスタリング処理を行い、空間位置を座標軸とする座標系における、レーダ信号点のクラスタを生成する。

【0041】

クラスタの生成の仕法は、特定の方法・方式に限定されるものではなく、例えば、信号点抽出部３１によって検出・抽出されたレーダ信号点の密度に基づくなどして物標についてのレーダ信号点の集まりをクラスタリングする、言い換えると、相互の位置関係が一定距離以内で近接しているレーダ信号点の集まりを括って抽出することができればどのような方法・方式であってもよい。クラスタの生成の仕法として、具体的には例えばＤＢＳＣＡＮ（Ｄensity－Ｂased Ｓpatial Ｃlustering of Ａpplications with Ｎoise の略）がクラスタリングアルゴリズムとして用いられるようにしてもよい。クラスタ生成部３２によって生成されるクラスタのそれぞれが１つの物標に対応する。

【0042】

クラスタ生成部３２によって生成される各クラスタを区別して一意に特定するための識別符号をＮとする（Ｎは、例えばシリアル番号）。そして、クラスタＮについての、当該のクラスタＮに含まれるレーダ信号点それぞれの相対速度（即ち、信号処理部２の速度計算部２３によって計算される各ドップラ周波数に対応するレーダ信号点の相対速度）の値の集まりのことを「クラスタＮの相対速度の値の集合」と呼ぶ。

【0043】

なお、物標位置検出部３は、検知対象エリアで検知された物標が、延いてはクラスタ生成部３２によって生成されたクラスタが、移動している人に関するものであるか否かを判別・識別した上で、移動している人に関するクラスタに対してのみ識別符号Ｎを付するようにしてもよい。検知された物標が移動している人であるか否かを判別・識別したり、生成されたクラスタが移動している人に関するものであるか否かを判別・識別したりする仕法は、公知の手法が存在し、また、この発明では特定の手法には限定されないので、ここでは詳細の説明は省略する。

【0044】

追尾処理部３３は、過去に検出されている物標を追尾する処理を、クラスタ生成部３２によって生成されたクラスタＮごと（即ち、物標ごと）に行う。つまり、追尾処理部３３は、連続する時点においてクラスタ生成部３２によって生成される、同一の物標に関するクラスタを関連づけることにより、過去のレーダスキャンによって検出されて識別符号Ｎが付されているクラスタの追尾処理を行う。

【0045】

クラスタ／物標の追尾の仕法は、公知の手法が存在し、また、この発明では特定の手法には限定されないので、ここでは詳細の説明は省略する。例えば、クラスタ／物標の追尾は、追尾フィルタとしてα－βフィルタ、カルマンフィルタ、或いはパーティクルフィルタを用いる手法などの公知の手法によって行われ得る。

【0046】

追尾処理部３３は、必要に応じ、クラスタ生成部３２によって生成されるクラスタそれぞれの代表点を決定した上で追尾処理を行うようにしてもよい。例えば、追尾処理部３３は、１つのクラスタとして認識された信号点群について、当該の信号点群に含まれる信号点それぞれの位置の、信号強度（受信レベル）Ｉ［ｄＢ］を重みとして用いた加重平均の位置を計算するなどして、前記１つのクラスタに含まれる複数の信号点を１点に集約した上で追尾処理を行うようにしてもよい。

【0047】

代表点特定部３４は、クラスタ生成部３２によって生成されるとともに追尾処理部３３によって追尾されるクラスタＮごとに、各レーダスキャンにおける、物標としてのクラスタＮの代表点を特定する。代表点特定部３４は、各クラスタＮの例えば重心位置を代表点の位置として特定して、クラスタＮごとに特定した代表点の位置座標Ｘ［ｍ］，Ｙ［ｍ］の組み合わせを出力する。

【0048】

物標計算部４は、物標位置検出部３における処理結果を用いて、特徴量計算部５による次の処理における関心領域を設定するとともに前記関心領域に関する情報を計算する。物標計算部４は、関心領域設定部４１、信号強度特定部４２、およびＶＴデータ生成部４３を含んで構成され、物標位置検出部３から出力される、物標位置検出部３における処理結果を用いて、特徴量計算部５による次の処理における関心領域を設定して前記関心領域に関する情報を計算する。

【0049】

関心領域設定部４１は、クラスタＮごとに、物標位置検出部３の代表点特定部３４から出力される各クラスタＮの代表点の位置を中心とする所定の範囲を関心領域として設定する。関心領域として設定される範囲は、特定の大きさや形状に限定されるものではなく、検知対象の物標の特に動作速度に纏わる特性を適切に把握することが考慮されるなどした上で、適当な大きさや形状に適宜設定される。関心領域は、具体的には例えば、クラスタＮの代表点の位置座標Ｘ［ｍ］，Ｙ［ｍ］を中心とする±１［ｍ］の四角形で囲まれる範囲に設定されることが考えられる。

【0050】

関心領域の大きさや形状は、予め固定的に設定されて関心領域設定部４１に格納されているようにしてもよく、或いは、関心領域設定部４１に対してユーザが適宜設定可能であるようにしてもよい。

【0051】

関心領域設定部４１は、さらに、信号処理部２から出力される４次元データ（Ｘ，Ｙ，Ｉ，Ｖ）の入力を受け、前記４次元データ（Ｘ，Ｙ，Ｉ，Ｖ）のうち、関心領域として設定されている範囲に含まれているデータのみを抜き出す。関心領域設定部４１は、つまり、関心領域として設定されている範囲に含まれている、全ての速度成分についての強度データを出力する。

【0052】

信号強度特定部４２は、関心領域設定部４１から出力される関心領域に関する４次元データ（Ｘ，Ｙ，Ｉ，Ｖ）を速度成分ごとに参照し、速度成分別に、関心領域に関する強度情報を生成する。信号強度特定部４２は、具体的には、クラスタＮの関心領域に関する４次元データ（Ｘ，Ｙ，Ｉ，Ｖ）について、相対速度Ｖ［ｋｍ／ｈ］ごとに、例えば、信号強度（受信レベル）Ｉ［ｄＢ］のうちの最大値Ｉmaxを特定して関心領域に関する強度情報としたり、信号強度Ｉ［ｄＢ］の平均値Ｉavgを計算して関心領域に関する強度情報としたりする。

【0053】

ここで、信号強度特定部４２によって特定されたり計算されたりする所定の特性をもった信号強度Ｉ［ｋｍ／ｈ］のことを「信号強度情報」と呼ぶ。つまり、信号強度特定部４２は、クラスタＮごとの、関心領域に関する、相対速度Ｖの値と信号強度情報の値（例えば、最大値Ｉmax，平均値Ｉavg［ｄＢ］）との組み合わせデータを出力する。

【0054】

ＶＴデータ生成部４３は、信号強度特定部４２から出力される関心領域に関する信号強度情報を用いて、各レーダスキャンの実行時刻に基づく経過時間Ｔ［秒］および相対速度Ｖ［ｋｍ／ｈ］と、これら経過時間Ｔ［秒］と相対速度Ｖ［ｋｍ／ｈ］との組み合わせに対応する信号強度情報（例えば、最大値Ｉmax，平均値Ｉavg［ｄＢ］）と、の組み合わせデータを生成する。

【0055】

信号処理部２、物標位置検出部３、および物標計算部４による処理は、レーダスキャンごとに行われる。すなわち、１回のレーダスキャンが行われるたびに、クラスタＮごとの、関心領域に関する、経過時間Ｔ［秒］と相対速度Ｖ［ｋｍ／ｈ］と信号強度情報（例えば、最大値Ｉmax，平均値Ｉavg［ｄＢ］）との組み合わせデータが作成される。

【0056】

特徴量計算部５は、物標計算部４における処理結果を用いて、物標としての移動している人の属性を識別するための特徴量を計算する。特徴量計算部５は、相関計算部５１、周期計算部５２、振幅計算部５３、および相関変換部５４を含んで構成され、物標計算部４から出力される、前記物標計算部４おける処理結果を用いて、検知対象エリア内を移動している人に関する特徴量を計算する。下記において説明する特徴量計算部５による処理は、その都度断らないが、クラスタＮごとに（言い換えると、クラスタＮ各々を単位として）行われる。

【0057】

相関計算部５１は、物標計算部４のＶＴデータ生成部４３から出力される経過時間Ｔ［秒］と相対速度Ｖ［ｋｍ／ｈ］との組み合わせに対応する信号強度情報（例えば、最大値Ｉmax，平均値Ｉavg［ｄＢ］）のデータを用いて相関分析を行う。

【0058】

相関計算部５１は、具体的には、まず、所定の時間長さぶん（言い換えると、レーダスキャンの所定の回数ぶん）のレーダスキャンごとの相対速度Ｖの値と信号強度情報の値との組み合わせデータの集合を時系列（言い換えると、レーダスキャンの実行順）に配列する。時系列で配列された、レーダスキャンごとの相対速度Ｖの値と信号強度情報の値との組み合わせデータの集合のことを「速度・強度の時系列データ」と呼ぶ。

【0059】

相関計算部５１から出力される速度・強度の時系列データの一例を図３に示す。図３は、各レーダスキャンの実行時刻に基づく経過時間Ｔ［秒］（即ち、１回のレーダスキャンあたりのサンプリング時間幅（例えば、５０ミリ秒間）とレーダスキャンの実行回数とを掛け合わせたもの）を横軸とするとともに相対速度Ｖ［ｋｍ／ｈ］の値を縦軸として、経過時間Ｔ［秒］と相対速度Ｖ［ｋｍ／ｈ］との組み合わせに対応する信号強度（受信レベル）Ｉ［ｄＢ］のうちの最大値Ｉmax［ｄＢ］を明暗で表して表示したものである。図３では、最大信号強度が、弱い場合は暗く表示され、強い場合は明るく表示されている。速度・強度の時系列データは、時系列における、物標に関する瞬時の速度の分布情報であるとも言える。なお、図３および図５～図１６は、２．５秒ぶん（言い換えると、５０回分のレーダスキャン）のデータが用いられて整理されている。

【0060】

ここで、レーダ装置を用いて得られる相対速度の情報は、物標としての移動している人の胴体速度だけでなく、頭、上腕、前腕・手先、大腿、および下腿・足先などの様々な部位の速度成分を同時に取得している。したがって、移動している人の頭、胴体、上腕、前腕・手先、大腿、および下腿・足先などの各部位の速度の情報が内在することになるため、レーダスキャンごとの相対速度Ｖの値と信号強度情報の値（図３に示す例では、最大信号強度の値）との組み合わせデータを時系列に並べることにより、各部位の運動／動作の周期が重畳して現れるようになる。

【0061】

相関計算部５１は、続いて、速度・強度の時系列データについて、レーダスキャンごとの相対速度Ｖの値と信号強度情報の値との組み合わせデータの時系列の配列を時間軸方向に１つずつシフトさせながら内積を計算して、シフト前の、時系列で配列された組み合わせデータの集合との対比における相関値を計算する（図４参照）。組み合わせデータの時系列の配列を時間軸方向に１つシフトさせる際、時系列の最後の時点に配置されている組み合わせデータは時系列の最初の時点へと配置される（図４参照）。相関値の計算は、時系列で配列された組み合わせデータの数と同じだけシフトを繰り返して行われる。なお、図４に示す例では、時系列で配列された組み合わせデータの数が５０であり、シフト実行回が５０回目のデータは、シフト前の、時系列で配列された組み合わせデータの集合と一致する。

【0062】

相関計算部５１は、速度・強度の時系列データを構成する、レーダスキャンごとの相対速度Ｖの値と信号強度情報の値との組み合わせデータを時間軸方向にシフトさせた回（「データのシフト回」とも表記する）ごとの相関値を出力する。相関計算部５１によって計算される、クラスタＮについてのデータのシフト回ごとの相関値の集まりのことを「クラスタＮの相関値の集合」と呼ぶ。

【0063】

相関計算部５１から出力されるクラスタＮの相関値の集合のデータの一例を図５に示す。図５は、データのシフト回を横軸とするとともに相関値を縦軸として、データのシフト回ごとの相関値の変動を折れ線で表して表示したものである。図５では、物標の運動の周期が重なったときに相関値のピークが現れる。

【0064】

周期計算部５２は、相関計算部５１から出力されるクラスタＮの相関値の集合を用いて、クラスタＮの相関値の集合を時系列（即ち、速度・強度の時系列データについての時間軸方向のシフトの実行順）に配列したときの、相関値のピーク間の差分の平均値を計算して出力する。

【0065】

クラスタＮの相関値の集合（即ち、時系列で配列された相対速度Ｖの値と信号強度情報の値との組み合わせデータを時間軸方向に１つずつシフトさせながら、シフト前の、時系列で配列された組み合わせデータの集合との対比における相関値の集まり）における相関値のピーク間の差分Ｄ１，Ｄ２，Ｄ３は、物標の動作の周期性を表す指標である（図６参照）。すなわち、相関値のピーク間の差分の平均値により、物標の運動周期（ここでは、移動している人の歩行に纏わる動作の周期）の遅速を把握することが可能となる。なお、シフト実行回が５０回目のデータは、シフト前の、時系列で配列された組み合わせデータの集合と一致するので相関値が１であり、周期計算部５２による処理においては相関値のピークとしては扱わない。

【0066】

振幅計算部５３は、相関計算部５１から出力されるクラスタＮの相関値の集合を用いて、クラスタＮの相関値の集合を時系列（即ち、速度・強度の時系列データについての時間軸方向のシフトの実行順）に配列したときの、相関値の極大値および極小値の平均と分散とを計算して出力する。

【0067】

クラスタＮの相関値の集合における相関値の極大値Ｅx１，Ｅx２，Ｅx３，Ｅx４および極小値Ｅn１，Ｅn２，Ｅn３，Ｅn４，Ｅn５は、物標の動作の振幅を表す指標である（図７参照）。すなわち、相関値の極大値および極小値の平均や分散により、物標の運動の振幅や変動（ここでは、移動している人の歩行に纏わる動作の振れ幅やばらつき）の大小を把握することが可能となる。なお、シフト実行回が５０回目のデータは、シフト前の、時系列で配列された組み合わせデータの集合と一致するので相関値が１であり、振幅計算部５３による処理においては相関値の極大値としては扱わない。

【0068】

相関変換部５４は、相関計算部５１から出力されるクラスタＮの相関値の集合を用いて、高速フーリエ変換を行い、Ｌ¹ノルムを計算して出力する。相関変換部５４は、具体的には、以下の処理を行う。

【0069】

１）相関計算部５１から出力される、所定の時間長さぶん（言い換えると、レーダスキャンの所定の回数ぶん）の、データのシフト回ごとの相関値の平均（「相関平均値」と呼ぶ）を算出する。
２）データのシフト回ごとの相関値と相関平均値との差（「シフト回ごとの相関差分」と呼ぶ）を算出する。
３）シフト回ごとの相関差分に対して高速フーリエ変換処理を施す。具体的には、シフト回ごとの相関値の大きさ（言い換えると、振幅）に対して高速フーリエ変換処理を施して相関値の振幅の分布を示すパワースペクトルを生成する。
４）高速フーリエ変換処理の結果の絶対値をとり、前記絶対値を２乗する。
５）上記４の結果として算出される、前記絶対値を２乗した値の集合のうちの最大値で、前記絶対値を２乗した値のそれぞれを割り、前記絶対値を２乗した値を０から１までの範囲で正規化する。

【0070】

上記の処理により、クラスタＮの相関値の集合の、高速フーリエ変換処理後のＬ¹ノルムが求められ、相関変換部５４は、高速フーリエ変換処理後のＬ¹ノルムとして、高速フーリエ変換後の周波数と振幅との組み合わせデータを出力する（図８参照）。

【0071】

上記により、特徴量計算部５は、クラスタＮごとの、相関値のピーク間の差分の平均値、相関値の極大値および極小値の平均と分散、さらに高速フーリエ変換後のＬ¹ノルムとしての周波数および振幅の組み合わせデータを出力する。

【0072】

識別部６は、特徴量計算部５から出力されるクラスタＮごとの各種指標の組み合わせデータに基づいて、クラスタＮごとに物標の属性を識別する。識別部６は、具体的には例えば、クラスタＮに対応する、物標としての移動している人の性別や年代を判別する。

【0073】

人の属性の判別について、識別部６は、例えば、機械学習のアルゴリズムの１つであるサポートベクターマシンを用いて判別処理を行うようにしてもよい。

【0074】

この場合は、識別部６は、例えば、クラスタＮの相関値の集合の、相関値のピーク間の差分の平均値ｘ１、相関値の極大値および極小値の平均の値ｘ２、ならびに相関値の極大値および極小値の分散の値ｘ３を変数として含むとともに係数ｆ１、ｆ２、およびｆ３を含む式（下記の数式１参照）を用いて人の属性の判別処理を行う。数式１について、例えば、クラスタＮが男性についてのレーダ信号点群である場合はｙ＝０とされ、クラスタＮが女性についてのレーダ信号点群である場合はｙ＝１とされる。
（数１） ｙ＝ｆ(ｘ１，ｘ２，ｘ３，ｆ１，ｆ２，ｆ３)

【0075】

数式１の係数ｆ１、ｆ２、およびｆ３は、サポートベクターマシンが用いられて予め決定される。係数ｆ１、ｆ２、およびｆ３の値は、男性についてのレーダ信号点群で構成されるクラスタＮの相関値の集合の、相関値のピーク間の差分の平均値ｘ１、相関値の極大値および極小値の平均の値ｘ２、ならびに相関値の極大値および極小値の分散の値ｘ３、ならびにｙの値（即ち、０）の組み合わせ（複数の組み合わせ）と、女性についてのレーダ信号点群で構成されるクラスタＮの相関値の集合の、相関値のピーク間の差分の平均値ｘ１、相関値の極大値および極小値の平均の値ｘ２、ならびに相関値の極大値および極小値の分散の値ｘ３、ならびにｙの値（即ち、１）の組み合わせ（複数の組み合わせ）と、を含む教師データが用意され、この教師データを男性についてのクラスタと女性についてのクラスタとに分離するための平面（「超平面」と呼ばれる）が最大のマージンを有するように決定される。

【0076】

識別部６は、上記によって予め決定された係数ｆ１、ｆ２、およびｆ３を含む数式１に、特徴量計算部５から出力されるクラスタＮの相関値のピーク間の差分の平均値ｘ１、相関値の極大値および極小値の平均の値ｘ２、ならびに相関値の極大値および極小値の分散の値ｘ３を代入して求められるｙの値に基づいて、前記クラスタＮが、男性についてのレーダ信号点群であるのか、或いは、女性についてのレーダ信号点群であるのかを判別する。

【0077】

ここで、識別部６は、クラスタＮの相関値の集合を時系列に配列したときの、相関値のピーク間の差分の平均値ｘ１、相関値の極大値および極小値の平均の値ｘ２、ならびに相関値の極大値および極小値の分散の値ｘ３のうちの少なくとも一つを用いて人の属性の判別処理を行うようにしてもよい。さらに、識別部６は、クラスタＮの相関値の集合の、高速フーリエ変換処理後のＬ¹ノルムのみを用いて人の属性の判別処理を行うようにしてもよく、或いは、上記平均値ｘ１、平均の値ｘ２、および分散の値ｘ３のうちの少なくとも一つと前記Ｌ¹ノルムとを組み合わせて用いて人の属性の判別処理を行うようにしてもよい。識別部６は、また、人の属性として、例えば、成人であるか子供であるかを判別するようにしてもよい。

【0078】

また、識別部６は、人の属性の判別について、機械学習のアルゴリズムの１つであるニューラルネットワーク、具体的には例えばＲＮＮ（Ｒeccurent Ｎeural Ｎetwork の略）であるＬＳＴＭ（Ｌong－Ｓhort Ｔerm Ｍemory の略）アルゴリズムを用いて判別処理を行うようにしてもよい。

【0079】

次に、このような構成の属性の識別装置１および属性の識別方法の動作や作用などについて図２も用いて説明する。

【0080】

まず、信号処理部２の周波数解析部２１が、レーダ装置から出力されるレーダデータ（別言すると、波形データ；具体的には、差分信号の振幅）としてのデジタル信号に対して高速フーリエ変換を用いて周波数解析を行って周波数スペクトルを生成し、差分信号のドップラ周波数ごとのスペクトル強度および位相情報を出力する（ステップＳ１）。

【0081】

次に、位置計算部２２が、周波数解析部２１から出力される位相情報に基づいて、ドップラ周波数ごとに、レーダ装置と検知対象の物標との間の距離を計算するとともにレーダ装置に対する前記物標の方位角度を計算し、計算の結果得られる距離情報および角度情報を出力する（ステップＳ２）。また、速度計算部２３が、周波数解析部２１における周波数解析で用いられたドップラ周波数を用いて、ドップラ周波数ごとに相対速度を計算して出力する（ステップＳ３）。

【0082】

次に、物標位置検出部３の信号点抽出部３１が、信号処理部２の周波数解析部２１によって生成される周波数スペクトルにおけるピークに対応するドップラ周波数（別言すると、ピーク周波数）を検出するとともに、位置計算部２２から出力されるドップラ周波数ごとの距離情報および角度情報を用いて、前記検出されたドップラ周波数（ピーク周波数）に対応するレーダ信号点の分布情報を生成して出力する（ステップＳ４）。また、クラスタ生成部３２が、前記分布情報を用いてクラスタリング処理を行い、レーダ信号点のクラスタを生成する（ステップＳ５）。さらに、代表点特定部３４が、クラスタＮの代表点（例えば、重心位置）を特定する（ステップＳ６）。この際、追尾処理部３３が、連続する時点においてクラスタ生成部３２によって生成される、同一の物標に関するクラスタを関連づけることにより、クラスタＮの追尾処理を行う。

【0083】

次に、物標計算部４の関心領域設定部４１が、クラスタＮの代表点の位置を中心とする所定の範囲を関心領域として設定した上で、関心領域に関する４次元データ（Ｘ，Ｙ，Ｉ，Ｖ）を出力する（ステップＳ７）。続いて、信号強度特定部４２が、関心領域に関する４次元データ（Ｘ，Ｙ，Ｉ，Ｖ）を速度成分ごとに参照し、速度成分別に、関心領域に関する信号強度情報を生成し（ステップＳ８）、また、ＶＴデータ生成部４３が、前記信号強度情報を用いて、各レーダスキャンの実行時刻に基づく経過時間Ｔ［秒］および相対速度Ｖ［ｋｍ／ｈ］と、これら経過時間Ｔ［秒］と相対速度Ｖ［ｋｍ／ｈ］との組み合わせに対応する信号強度情報（例えば、最大値Ｉmax，平均値Ｉavg［ｄＢ］）と、の組み合わせデータを生成する。

【0084】

次に、特徴量計算部５の相関計算部５１が、物標計算部４のＶＴデータ生成部４３から出力される経過時間Ｔ［秒］と相対速度Ｖ［ｋｍ／ｈ］との組み合わせに対応する信号強度情報（例えば、最大値Ｉmax，平均値Ｉavg［ｄＢ］）のデータを用いて、所定の時間長さぶんのレーダスキャンごとの前記データの集合を時系列に配列して速度・強度の時系列データを作成し、作成した速度・強度の時系列データについて、レーダスキャンごとの相対速度Ｖの値と信号強度情報の値との組み合わせデータの時系列の配列を時間軸方向に１つずつシフトさせながら、シフト前の、時系列で配列された組み合わせデータの集合との対比における相関値を計算する（ステップＳ９）。

【0085】

続いて、周期計算部５２が、相関計算部５１から出力されるクラスタＮの相関値の集合を時系列（即ち、速度・強度の時系列データについての時間軸方向のシフトの実行順）に配列したときの、相関値のピーク間の差分の平均値を計算する。また、振幅計算部５３が、相関計算部５１から出力されるクラスタＮの相関値の集合を時系列（即ち、速度・強度の時系列データについての時間軸方向のシフトの実行順）に配列したときの、相関値の極大値および極小値の平均と分散とを計算する。さらに、相関変換部５４が、相関計算部５１から出力されるクラスタＮの相関値の集合に対して高速フーリエ変換を行ってＬ¹ノルムを計算する（ステップＳ１０）。

【0086】

次に、識別部６が、特徴量計算部５から出力される、クラスタＮごとの、相関値のピーク間の差分の平均値、相関値の極大値および極小値の平均と分散、ならびに高速フーリエ変換後のＬ¹ノルムとしての周波数および振幅の値に基づいて、クラスタＮごとに物標としての移動している人の属性を識別する（ステップＳ１１）。なお、具体的には、相関値のピーク間の差分の平均値、相関値の極大値および極小値の平均と分散、ならびに高速フーリエ変換後のＬ¹ノルムのうちの少なくとも一つに基づいて、物標としての移動している人の属性を識別するようにすればよい。

【0087】

ここで、人の属性の違いによる、レーダデータに表れる特性の比較の例を説明する。

【0088】

図９は、物標が男性である場合の、５０回のレーダスキャンによって取得されたレーダデータについてのクラスタＮの相関値の集合と、物標が女性である場合の、５０回のレーダスキャンによって取得されたレーダデータについてのクラスタＮの相関値の集合とのそれぞれに関して、データのシフト回ごとの相関値の変動を折れ線で表して表示して１つのグラフ上に整理したグラフである。図９から、男性と女性とではデータのシフト回ごとの相関値の変動に明らかな違いがあることが看取され、データのシフト回ごとの相関値の変動に基づいて、男性であるか女性であるかを判別することが、延いては人の属性を識別することが可能であることが確認される。

【0089】

図１０は、物標が男性である場合の、５０回のレーダスキャンによって取得されたレーダデータについてのクラスタＮの相関値の集合と、物標が女性である場合の、５０回のレーダスキャンによって取得されたレーダデータについてのクラスタＮの相関値の集合とのそれぞれに関して（尚、レーダデータは図９に関連して用いたものと同じである）、高速フーリエ変換処理後のＬ¹ノルムの値を折れ線で表示して１つのグラフ上に整理したグラフである。図１０から、男性と女性とではＬ¹ノルムの値に明らかな違いがあることが看取され、Ｌ¹ノルムの値に基づいて男性であるか女性であるかを判別することが、延いては人の属性を識別することが可能であることが確認される。

【0090】

図１１は、相関値のピーク間の差分の平均値（図では「周期」と表記）、相関値の極大値および極小値の平均の値、ならびに相関値の極大値および極小値の分散の値についての直交３軸座標系に、レーダデータから検出・抽出されたクラスタについて計算される値に基づいてクラスタごとにプロットした結果である。図１１に示すクラスタのそれぞれについては、男性についてのクラスタであるのか女性についてのクラスタであるのかが分かっている。図１１に示す結果から、男性についてのクラスタと女性についてのクラスタとに分離するための平面が存在すると考えられる。

【0091】

図１２は、図１１に示すデータについてサポートベクターマシンを用いて判別処理を行って超平面を計算した結果である。図１２に示す結果から、相関値のピーク間の差分の平均値、相関値の極大値および極小値の平均の値、ならびに相関値の極大値および極小値の分散の値を適切に組み合わせて用いることにより、機械学習のアルゴリズムを用いて男性であるか女性であるかを判別することが、延いては人の属性を識別することが可能であることが確認される。

【0092】

以上、この発明の実施の形態について説明したが、具体的な構成は、上記の実施の形態に限られるものではなく、この発明の要旨を逸脱しない範囲の設計の変更等があっても、この発明に含まれる。具体的には、この発明の要旨は、所定の時間長さぶんのレーダスキャンごとの相対速度Ｖの値と信号強度情報の値との組み合わせデータの集合を時系列に配列した上で、前記組み合わせデータの時系列の配列を時間軸方向に１つずつシフトさせながら、シフト前の、時系列で配列された前記組み合わせデータの集合との対比として計算される相関値を用いて物標の属性を判別することであり、この要旨を逸脱しない範囲の設計の変更等があっても、この発明に含まれる。例えば、この発明が適用される装置の具体的な構成は上記の実施の形態における属性の識別装置１に限定されるものではなく、具体的な装置構成は種々の態様があり得る。

【0093】

また、上記の実施の形態では特徴量計算部５がクラスタＮの相関値の集合を時系列に配列したときの、相関値のピーク間の差分の平均値、相関値の極大値および極小値の平均、ならびに相関値の極大値および極小値の分散を計算するとともに識別部６が前記指標のうちの少なくとも一つを用いて人の属性の判別処理を行うようにしているが、相関計算部５１によって計算されて出力されるクラスタＮの相関値の集合を用いて計算される、物標としての移動している人の属性を識別するための特徴量としての指標は前記したものには限定されない。例えば、クラスタＮの相関値の集合における相関値の極大値および極小値の中央値やレンジなどが計算されて判別処理に用いられるようにしてもよい。

【0094】

また、代表点特定部３４および関心領域設定部４１が人の特定の部位に限定して関心領域を設定して前記関心領域に関して取得されたレーダデータが用いられるようにしてもよい。この点に関連して、図１３は、物標を手先に限定した場合の、５０回のレーダスキャンによって取得されたレーダデータについてのクラスタＮの相関値の集合に関して、データのシフト回ごとの相関値の変動を折れ線で表して表示したグラフであり、（Ａ）は成人についての結果であり、（Ｂ）は子供についての結果である。これら図１３の（Ａ）と（Ｂ）との比較から、成人と子供とではデータのシフト回ごとの相関値の変動に明らかな違いがあることが看取され、物標が手先に限定されたレーダデータを用いることにより、或いは、物標として手先を含むレーダデータを用いることにより、データのシフト回ごとの相関値の変動に基づいて成人であるか子供であるかを判別することが、延いては人の属性を識別することが可能であることが確認される。また、図１３の（Ａ）と（Ｂ）との整理から、相関値の極大値および極小値の分散が、成人の手先は１．１７６０となり、子供の手先は０．０７７３となる。また、相関値の極大値および極小値の平均が、成人の手先は０．９９８２となり、子供の手先は０．９９９１となる。これらの結果から、成人と子供とでは相関値の極大値および極小値の特に分散が大きく異なるので、物標が手先に限定されたレーダデータを用いることにより、或いは、物標として手先を含むレーダデータを用いることにより、データのシフト回ごとの相関値の極大値および極小値の分散に基づいて成人であるか子供であるかを判別することが、延いては人の属性を識別することが可能であることが確認される。

【0095】

図１４は、物標をももに限定した場合の、５０回のレーダスキャンによって取得されたレーダデータについてのクラスタＮの相関値の集合に関して、データのシフト回ごとの相関値の変動を折れ線で表して表示したグラフであり、（Ａ）は成人についての結果であり、（Ｂ）は子供についての結果である。これら図１４の（Ａ）と（Ｂ）との比較から、成人と子供とではデータのシフト回ごとの相関値の変動に明らかな違いがあることが看取され、物標がももに限定されたレーダデータを用いることにより、或いは、物標としてももを含むレーダデータを用いることにより、データのシフト回ごとの相関値の変動に基づいて成人であるか子供であるかを判別することが、延いては人の属性を識別することが可能であることが確認される。また、図１４の（Ａ）と（Ｂ）との整理から、相関値の極大値および極小値の分散が、成人のももは１．８９２７となり、子供のももは０．２０７８となる。また、相関値の極大値および極小値の平均が、成人のももは０．９９８８となり、子供のももは０．９９８８となる。これらの結果から、成人と子供とでは相関値の極大値および極小値の特に分散が大きく異なるので、物標がももに限定されたレーダデータを用いることにより、或いは、物標としてももを含むレーダデータを用いることにより、データのシフト回ごとの相関値の極大値および極小値の分散に基づいて成人であるか子供であるかを判別することが、延いては人の属性を識別することが可能であることが確認される。

【0096】

図１５は、物標を手先に限定した場合の、５０回のレーダスキャンによって取得されたレーダデータについてのクラスタＮの相関値の集合に関して、高速フーリエ変換処理後のＬ¹ノルム（具体的には、周波数と振幅との組み合わせデータ）を折れ線で表示したグラフであり、（Ａ）は成人についての結果であり、（Ｂ）は子供についての結果である。これら図１５の（Ａ）と（Ｂ）との比較から、成人と子供とではＬ¹ノルムの値に明らかな違いがあることが看取され、物標が手先に限定されたレーダデータを用いることにより、或いは、物標として手先を含むレーダデータを用いることにより、Ｌ¹ノルムの値に基づいて成人であるか子供であるかを判別することが、延いては人の属性を識別することが可能であることが確認される。

【0097】

図１６は、物標をももに限定した場合の、５０回のレーダスキャンによって取得されたレーダデータについてのクラスタＮの相関値の集合に関して、高速フーリエ変換処理後のＬ¹ノルム（具体的には、周波数と振幅との組み合わせデータ）を折れ線で表示したグラフであり、（Ａ）は成人についての結果であり、（Ｂ）は子供についての結果である。これら図１６の（Ａ）と（Ｂ）との比較から、成人と子供とではＬ¹ノルムの値に明らかな違いがあることが看取され、物標がももに限定されたレーダデータを用いることにより、或いは、物標としてももを含むレーダデータを用いることにより、Ｌ¹ノルムの値に基づいて成人であるか子供であるかを判別することが、延いては人の属性を識別することが可能であることが確認される。

【産業上の利用可能性】

【0098】

この発明は、簡便な機器構成によって歩行者の性別や年代などの属性の識別を高精度に行うことができるので、例えば、自動車の運転支援に関する分野や侵入者検知に関する分野において有用な技術として適用され得る。

【符号の説明】

【0099】

１ 属性の識別装置
２ 信号処理部
２１ 周波数解析部
２２ 位置計算部
２３ 速度計算部
３ 物標位置検出部
３１ 信号点抽出部
３２ クラスタ生成部
３３ 追尾処理部
３４ 代表点特定部
４ 物標計算部
４１ 関心領域設定部
４２ 信号強度特定部
４３ ＶＴデータ生成部
５ 特徴量計算部
５１ 相関計算部
５２ 周期計算部
５３ 振幅計算部
５４ 相関変換部
６ 識別部

【図1】

【図2】

【図3】

【図4】

【図5】

【図6】

【図7】

【図8】

【図9】

【図10】

【図11】

【図12】

【図13】

【図14】

【図15】

【図16】