(57)【要約】【課題】あらかじめ閾値を設定せずとも、高精度に物体を識別することを可能とする技術が提供されることが望まれる。【解決手段】センサによる計測に基づいて複素数の多次元テンソルを取得する取得部と、前記複素数の多次元テンソルに基づく畳込み処理および次元圧縮それぞれを複数回行ってテンソル要素の特徴を抽出する特徴抽出部と、前記テンソル要素の特徴の次元を逆畳込み処理によって前記多次元テンソルの次元まで拡張して再拡張多次元テンソルを生成する拡張部と、前記再拡張多次元テンソルのテンソル要素ごとにクラスを識別することによって領域分割を行う領域分割部と、を備える、識別装置が提供される。【選択図】図２
2019074478.tif 000002
【特許請求の範囲】
【請求項１】
センサによる計測に基づいて複素数の多次元テンソルを取得する取得部と、
前記複素数の多次元テンソルに基づく畳込み処理および次元圧縮それぞれを複数回行ってテンソル要素の特徴を抽出する特徴抽出部と、
前記テンソル要素の特徴の次元を逆畳込み処理によって前記多次元テンソルの次元まで拡張して再拡張多次元テンソルを生成する拡張部と、
前記再拡張多次元テンソルのテンソル要素ごとにクラスを識別することによって領域分割を行う領域分割部と、
を備える、識別装置。
【請求項２】
前記特徴抽出部は、前記複素数の多次元テンソルの各テンソル要素の絶対値に対して畳込み処理および次元圧縮それぞれを複数回行う、
請求項１に記載の識別装置。
【請求項３】
前記特徴抽出部は、前記複素数の多次元テンソルの各テンソル要素自体に対して畳込み処理および次元圧縮それぞれを複数回行うことによって複素数によるテンソル要素の特徴を抽出し、前記複素数によるテンソル要素の特徴を実数によるテンソル要素の特徴に変換する、
請求項１に記載の識別装置。
【請求項４】
前記領域分割部は、前記テンソル要素ごとのクラスの識別結果に基づいて、物体と偽像とを区別する、
請求項１に記載の識別装置。
【請求項５】
前記取得部は、前記センサによる計測に基づく２次元情報に対応付けられた複素信号を各テンソル要素として有する前記複素数の多次元テンソルを取得する、
請求項１に記載の識別装置。
【請求項６】
前記２次元情報は、前記センサを基準とした物体までの距離と角度との組み合わせ、または、前記センサを基準とした物体までの距離と前記物体の速度との組み合わせによって構成される、
請求項１に記載の識別装置。
【請求項７】
前記特徴抽出部は、学習処理によって得られたパラメータを前記畳込み処理に利用する、
請求項１に記載の識別装置。
【請求項８】
前記特徴抽出部は、高分解能化処理に基づいて生成された教師データを用いた学習処理によって得られたパラメータを前記畳込み処理に利用する、
請求項７に記載の識別装置。
【請求項９】
前記特徴抽出部は、テンソル要素の特徴数をクラスの数に合わせ、
前記拡張部は、前記クラスの数に特徴数が合わせられたテンソル要素の特徴の次元を拡張する、
請求項１に記載の識別装置。
【請求項１０】
センサによる計測に基づいて複素数の多次元テンソルを取得することと、
前記複素数の多次元テンソルに基づく畳込み処理および次元圧縮それぞれを複数回行ってテンソル要素の特徴を抽出することと、
前記テンソル要素の特徴の次元を逆畳込み処理によって前記多次元テンソルの次元まで拡張して再拡張多次元テンソルを生成することと、
前記再拡張多次元テンソルのテンソル要素ごとにクラスを識別することによって領域分割を行うことと、
を含む、識別方法。
【請求項１１】
コンピュータを、
センサによる計測に基づいて複素数の多次元テンソルを取得する取得部と、
前記複素数の多次元テンソルに基づく畳込み処理および次元圧縮それぞれを複数回行ってテンソル要素の特徴を抽出する特徴抽出部と、
前記テンソル要素の特徴の次元を逆畳込み処理によって前記多次元テンソルの次元まで拡張して再拡張多次元テンソルを生成する拡張部と、
前記再拡張多次元テンソルのテンソル要素ごとにクラスを識別することによって領域分割を行う領域分割部と、
を備える識別装置として機能させるためのプログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、識別装置、識別方法およびプログラムに関する。
【背景技術】
【０００２】
近年、レーダによる計測が利用されている。レーダによる計測は、アンテナから送信された電波が物体に当たって反射した場合、その反射波の位相シフトを解析することによって可能となる。例えば、ＦＭ−ＣＷ（Frequency-Modulated continuous wave）レーダでは、ビート信号の周波数はアンテナから物体までの距離に比例するため、信号強度が大きい周波数に対応する距離に物体があると判断することができる。
【０００３】
しかし、物体からの直接の反射波だけではなく、当該物体から他の物体を経由した反射波が受信されてしまう場合も想定される。このような場合には、当該他の物体を経由した反射波に基づく計測によって、実際には物体が存在しない場所に物体があると誤認識されてしまうという偽像（「ゴースト」とも言う。）発生の問題が生じ得る。また、偽像はフーリエ変換のサイドローブによっても発生し得る。
【０００４】
ここで、偽像を除去または抑制するための技術として、様々な技術が開示されている。例えば、第１の間隔で配置された複数の受信アンテナ素子での受信信号の位相差から求めた反射波の到来方向と、第１の間隔とは異なる第２の間隔で配置された複数の受信アンテナ素子での受信信号の位相差から求めた反射波の到来方向とが一致するか否かによって、偽像が発生していないか否かを判定する技術が開示されている（例えば、特許文献１参照）。
【０００５】
また、電波の送信信号と電波の反射波の受信信号とからビート信号を生成し、ビート信号のピーク周波数における振幅値（信号強度）とピーク周波数とに基づいて、物標の存在確率を測定サイクル毎に算出し、算出した存在確率が所定値に達したか否かに応じて、偽像が発生していないか否かを判定する技術が開示されている（例えば、特許文献２参照）。
【先行技術文献】
【特許文献】
【０００６】
【特許文献１】特開２０００−２３０９７４号公報
【特許文献２】特開２００９−２１０３３７号公報
【発明の概要】
【発明が解決しようとする課題】
【０００７】
既に開示されている上記技術などにおいては、ビート信号のピーク周波数の振幅値と閾値との比較結果に基づいて、偽像が発生していないか否かを判定する必要がある。しかし、ビート信号のピーク周波数における振幅値は、電波を反射する物体（反射物体）の材質、電波を反射する部分の面積（反射断面積）、アンテナと物体との距離などによって変化する可能性がある。そのため、あらゆる状況に適用可能な閾値を設定するのが困難である。
【０００８】
そこで、あらかじめ閾値を設定せずとも、高精度に物体を識別することを可能とする技術が提供されることが望まれる。
【課題を解決するための手段】
【０００９】
上記問題を解決するために、本発明のある観点によれば、センサによる計測に基づいて複素数の多次元テンソルを取得する取得部と、前記複素数の多次元テンソルに基づく畳込み処理および次元圧縮それぞれを複数回行ってテンソル要素の特徴を抽出する特徴抽出部と、前記テンソル要素の特徴の次元を逆畳込み処理によって前記多次元テンソルの次元まで拡張して再拡張多次元テンソルを生成する拡張部と、前記再拡張多次元テンソルのテンソル要素ごとにクラスを識別することによって領域分割を行う領域分割部と、を備える、識別装置が提供される。
【００１０】
前記特徴抽出部は、前記複素数の多次元テンソルの各テンソル要素の絶対値に対して畳込み処理および次元圧縮それぞれを複数回行ってよい。
【００１１】
前記特徴抽出部は、前記複素数の多次元テンソルの各テンソル要素自体に対して畳込み処理および次元圧縮それぞれを複数回行うことによって複素数によるテンソル要素の特徴を抽出し、前記複素数によるテンソル要素の特徴を実数によるテンソル要素の特徴に変換してよい。
【００１２】
前記領域分割部は、前記テンソル要素ごとのクラスの識別結果に基づいて、物体と偽像とを区別してよい。
【００１３】
前記取得部は、前記センサによる計測に基づく２次元情報に対応付けられた複素信号を各テンソル要素として有する前記複素数の多次元テンソルを取得してよい。
【００１４】
前記２次元情報は、前記センサを基準とした物体までの距離と角度との組み合わせ、または、前記センサを基準とした物体までの距離と前記物体の速度との組み合わせによって構成されてよい。
【００１５】
前記特徴抽出部は、学習処理によって得られたパラメータを前記畳込み処理に利用してよい。
【００１６】
前記特徴抽出部は、高分解能化処理に基づいて生成された教師データを用いた学習処理によって得られたパラメータを前記畳込み処理に利用してよい。
【００１７】
前記特徴抽出部は、テンソル要素の特徴数をクラスの数に合わせ、前記拡張部は、前記クラスの数に特徴数が合わせられたテンソル要素の特徴の次元を拡張してよい。
【００１８】
また、本発明の別の観点によれば、センサによる計測に基づいて複素数の多次元テンソルを取得することと、前記複素数の多次元テンソルに基づく畳込み処理および次元圧縮それぞれを複数回行ってテンソル要素の特徴を抽出することと、前記テンソル要素の特徴の次元を逆畳込み処理によって前記多次元テンソルの次元まで拡張して再拡張多次元テンソルを生成することと、前記再拡張多次元テンソルのテンソル要素ごとにクラスを識別することによって領域分割を行うことと、を含む、識別方法が提供される。
【００１９】
また、本発明の別の観点によれば、コンピュータを、センサによる計測に基づいて複素数の多次元テンソルを取得する取得部と、前記複素数の多次元テンソルに基づく畳込み処理および次元圧縮それぞれを複数回行ってテンソル要素の特徴を抽出する特徴抽出部と、前記テンソル要素の特徴の次元を逆畳込み処理によって前記多次元テンソルの次元まで拡張して再拡張多次元テンソルを生成する拡張部と、前記再拡張多次元テンソルのテンソル要素ごとにクラスを識別することによって領域分割を行う領域分割部と、を備える識別装置として機能させるためのプログラムが提供される。
【発明の効果】
【００２０】
以上説明したように本発明によれば、あらかじめ閾値を設定せずとも、高精度に物体を識別することを可能とする技術が提供される。
【図面の簡単な説明】
【００２１】
【図１】偽像の発生を説明するための図である。
【図２】第１の実施形態に係る識別装置の構成例を示す図である。
【図３】信号マップの絶対値を画像化した一例を示す図である。
【図４】第１の実施形態に係るクラス識別の処理例を示すフローチャートである。
【図５】学習処理に利用される信号マップの絶対値が画像化された一例を示す図である。
【図６】図５に示した画像に対応する教師データの例を示す図である。
【図７】第２の実施形態に係る識別装置の構成例を示す図である。
【図８】第２の実施形態に係るクラス識別の処理例を示すフローチャートである。
【図９】本発明の実施形態に係る識別装置の例としての情報処理装置のハードウェア構成を示す図である。
【発明を実施するための形態】
【００２２】
以下に添付図面を参照しながら、本発明の好適な実施の形態について詳細に説明する。なお、本明細書及び図面において、実質的に同一の機能構成を有する構成要素については、同一の符号を付することにより重複説明を省略する。
【００２３】
また、本明細書および図面において、実質的に同一の機能構成を有する複数の構成要素を、同一の符号の後に異なる数字を付して区別する場合がある。ただし、実質的に同一の機能構成を有する複数の構成要素等の各々を特に区別する必要がない場合、同一符号のみを付する。また、異なる実施形態の類似する構成要素については、同一の符号の後に異なるアルファベットを付して区別する場合がある。ただし、異なる実施形態の類似する構成要素等の各々を特に区別する必要がない場合、同一符号のみを付する。
【００２４】
（０．概略）
近年、レーダによる計測が利用されている。レーダによる計測は、アンテナから送信された電波が物体に当たって反射した場合、その反射波の位相シフトを解析することによって可能となる。例えば、ＦＭ−ＣＷレーダでは、ビート信号の周波数はアンテナから物体までの距離に比例するため、信号強度が大きい周波数に対応する距離に物体があると判断することができる。しかし、かかる技術によれば、実際には物体がない場所に物体があると誤認識されてしまうという偽像発生の問題が生じ得る。以下、偽像の発生について、図面を参照しながら説明する。
【００２５】
図１は、偽像の発生を説明するための図である。図１を参照すると、実空間に物体９１が存在している。また、実空間には他の物体９０も存在している。ここで、物体９１からの直接の反射波が受信されれば、その反射波の位相シフトを解析することによって、物体９１が存在する位置を正常に認識することが可能である。しかし、図１に示したように、実際には物体９１からの反射波が他の物体９０を経由して受信されてしまう場合も想定される。
【００２６】
このような場合には、他の物体９０を経由した反射波に基づく計測によって、実際には物体９１が存在しない場所（図１に示した例では、偽像９２が存在する場所）に物体９１があると誤認識されてしまうという偽像発生の問題が生じ得る。また、偽像はフーリエ変換のサイドローブによっても発生し得る。
【００２７】
ここで、偽像を除去または抑制するための技術としては、既に説明したように様々な技術が開示されている。既に開示されている上記技術などにおいては、ビート信号のピーク周波数の振幅値と閾値との比較結果に基づいて、偽像が発生していないか否かを判定する必要がある。しかし、ビート信号のピーク周波数における振幅値は、電波を反射する物体（反射物体）の材質、電波を反射する部分の面積（反射断面積）、アンテナと物体との距離などによって変化する可能性がある。そのため、あらゆる状況に適用可能な閾値を設定するのが困難である。
【００２８】
そこで、本明細書においては、閾値をあらかじめ設定せずとも、高精度に物体を識別することを可能とする技術について主に説明する。
【００２９】
（１．第１の実施形態）
まず、第１の実施形態について説明する。
【００３０】
（１．１．構成）
第１の実施形態に係る識別装置１０Ａの構成例について説明する。
【００３１】
図２は、第１の実施形態に係る識別装置１０Ａの構成例を示す図である。図２に示すように、第１の実施形態に係る識別装置１０Ａは、送受信部１１１、取得部１１２、特徴抽出部１１３Ａ、拡張部１１４および領域分割部１１５を備える。ここで、送受信部１１１は、センサによる計測を行い、計測結果を取得部１１２に出力する。そして、取得部１１２は、センサによる計測に基づいて複素数の多次元テンソルを取得する。
【００３２】
特徴抽出部１１３Ａは、複素数の多次元テンソルに基づく畳込み処理および次元圧縮それぞれを複数回行ってテンソル要素の特徴を抽出する。拡張部１１４は、テンソル要素の特徴の次元を逆畳込み処理によって前記多次元テンソルの次元まで拡張して再拡張多次元テンソルを生成する。領域分割部１１５は、再拡張多次元テンソルのテンソル要素ごとにクラスを識別することによって領域分割を行う。
【００３３】
センサの種類は、特に限定されない。以下では、電波を用いたレーダ技術が用いられる場合を主に説明する。かかる場合には、電波を送信するアンテナ（送信アンテナ）と電波を受信するアンテナ（受信アンテナ）とを有するレーダがセンサに該当し得る。しかし、電波を用いたレーダ技術の代わりに、他の技術が用いられてもよい。例えば、電波を用いたレーダ技術の代わりに、音波を用いたソナー技術などが用いられてもよい。音波を用いたソナー技術が用いられる場合、音波受信器がセンサに該当し得る。
【００３４】
以下では、取得部１１２が、センサによる計測に基づく２次元情報に対応付けられた複素信号を各テンソル要素として有する複素数の多次元テンソルを取得する場合を主に説明する。センサによる計測に基づく２次元情報の種類は、特に限定されない。以下では、センサによる計測に基づく２次元情報が、センサから物体までの距離と角度との組み合わせによって構成される場合を主に説明する。しかし、センサによる計測に基づく２次元情報は、センサから物体までの距離と物体の速度との組み合わせによって構成されてもよい。
【００３５】
センサから物体までの距離は、どのようにして求められてもよい。例えば、センサから物体までの距離は、送信アンテナから送信された電波が物体で反射して受信アンテナによって受信されるまでの伝搬遅延から求められてもよい。あるいは、センサから物体までの距離は、周波数が異なる複数の電波が送信された場合、各電波の物体からの反射波によるドップラー信号の位相に基づいて求められてもよい。あるいは、センサから物体までの距離は、送信波の周波数を連続的に変化させ、物体での反射による伝搬遅延に基づき生じる周波数差から求められてもよい。
【００３６】
センサから物体までの角度（方位）は、どのようにして求められてもよい。センサから物体までの角度は、複数のアンテナ間の位相差またはパワー差から求められてもよい。あるいは、センサから物体までの角度は、指向性を有した電波（ビーム）を形成し、形成したビームを電子的または機械的に動かしながら送信した場合に、物体から反射されて検出されるビームの方向に基づいて検出されてもよい。
【００３７】
特に、第１の実施形態においては、特徴抽出部１１３Ａが、複素数の多次元テンソルの各テンソル要素の絶対値に対して畳込み処理および次元圧縮それぞれを複数回行う場合を説明する。各テンソル要素の絶対値は、センサによる計測に基づく２次元情報（例えば、センサによる計測に基づく距離および角度）に対応付けられた複素信号の振幅値（信号強度）に該当し得る。
【００３８】
以上、第１の実施形態に係る識別装置１０Ａの構成例について説明した。
【００３９】
（１．２．動作）
続いて、第１の実施形態に係る識別装置１０Ａの動作例について説明する。まず、送受信部１１１は、送信アンテナによって電波（送信信号）を送信し、送信アンテナによって送信した電波の反射波（受信信号）を受信アンテナによって受信する。送信アンテナの数は、１つであってもよいが、複数であるのがよい。また、受信アンテナの数は、１つであってもよいが、複数であるのがよい。
【００４０】
複数の送信アンテナおよび複数の受信アンテナがＭＩＭＯ（Multiple Input Multiple Output）レーダに用いられた場合、仮想的なアレイアンテナが構成され、角度の検知分解能が向上され得ることが知られている。ＭＩＭＯレーダの仮想的なアレイアンテナの構成については、以下の参考文献１を参照されたい。
【００４１】
参考文献１：菅原  博樹（日本無線株式会社）「電波環境適応レーダーの研究開発」，電波資源拡大のための研究開発，第７回成果発表会，２０１４年
【００４２】
取得部１１２は、送受信部１１１によって送信された電波（送信信号）と送受信部１１１によって受信された反射波（受信信号）とをミキシング処理することによって、ビート信号ｂｒ，ｔを得る。ここで、ｔは、サンプリング時刻を示し、ｒは、仮想アレイの並び順の番号を示す。そして、取得部１１２は、得られたビート信号ｂｒ，ｔを、下記の（数式１）および（数式２）によって、距離ｄ−角度θの２次元の（Ｎ，Ｍ）行列に変換する。
【００４３】
【数１】
2019074478.tif 000003
【００４４】
ここで、Ｎは、距離方向の信号サンプル数を示し、Ｍは、角度方向の信号サンプル数を示し、Ｋは、仮想的な受信アレイアンテナの数を示す。
【００４５】
以上のように、取得部１１２は、ビート信号ｂｒ，ｔを距離ｄ−角度θの２次元平面に分解することによって２次元の（Ｎ，Ｍ）行列を得る。２次元の（Ｎ，Ｍ）行列は、複素数を要素として有している（すなわち、２次元の（Ｎ，Ｍ）行列は、複素信号を各テンソル要素として有する多次元テンソルの一例に相当し得る）。以下では、この２次元の（Ｎ，Ｍ）行列を「信号マップ」と呼ぶこともある。
【００４６】
ここで、信号マップＢｔ（ｄ，θ）の各要素（複素数）の絶対値は振幅値（信号強度）を示している。図３は、信号マップｂｔ（ｄ，θ）の絶対値を画像化した一例を示す図である。図３を参照すると、画像７０−１が示されている。画像７０−１の縦軸は、レーダ中心方向からの角度に相当し、画像７０−１の横軸は、レーダからの距離に相当する。また、画像７０−１において、物体からの反射波の強度が大きい領域ほど、輝度が高くなっている（白っぽくなっている）。一方、物体からの反射波の強度が小さい領域ほど、輝度が低くなっている（黒っぽくなっている）。
【００４７】
より具体的に、画像７０−１において、物体からの反射波の強度の最も大きい領域７１−１は、実際に物体が存在する領域である（すなわち、物体からの直接の反射波の受信によって輝度が最も高くなっている）。一方、画像７０−１において、物体からの反射波の強度が大きい他の領域７２−１〜７２−３は、偽像が発生した領域である（すなわち、物体から他の物体を経由した反射波の受信によって輝度が高くなっている）。
【００４８】
画像７０−１から把握されるように、実際に物体が存在する領域７１−１の信号強度は大きい値を示すこととなる。一方、画像７０−１から把握されるように、実際には物体が存在しない領域７２−１〜７２−３の信号強度も、偽像の発生によって大きめの値を示してしまうこととなる。
【００４９】
続いて、信号マップが有する要素（座標）ごとにクラスが識別される。クラスの種類は特に限定されない。以下では、クラスの種類として、物体に対応するクラスと背景に対応するクラスを主に想定する。背景は、対象とする物体が存在しない領域であり、偽像が発生している領域も含み得る。かかる場合には、領域分割部１１５は、クラスの識別結果に基づいて、物体と偽像とを区別することができる。また、以下においても説明するように、物体の種類は、１つであってもよいし、複数であってもよい。物体の種類が複数である場合、複数の種類それぞれの物体に対応するクラスが存在してよい。
【００５０】
以下、クラス識別の処理例について詳細に説明する。なお、以下で説明するクラス識別の処理例は、以下の参考文献２に記載されたFully convolutional network（ＦＣＮ）をもとにした手法である。すなわち、参考文献２に記載されたFully convolutional network（ＦＣＮ）は、画像のピクセル単位にクラスを識別する手法であり、本実施形態に係るクラス識別の処理例は、この画像のピクセル単位にクラスを識別する手法と同様の考え方を信号マップに適用したものである。
【００５１】
参考文献２：Long, Jonathan,
Evan Shelhamer, and Trevor Darrell. "Fully convolutional networks for
semantic segmentation." Proceedings of the IEEE Conference on Computer
Vision and Pattern Recognition. 2015.
【００５２】
図４は、第１の実施形態に係るクラス識別の処理例を示すフローチャートである。まず、特徴抽出部１１３Ａは、取得部１１２によって得られた信号マップＢｔ（ｄ，θ）（ただし、１≦ｄ≦Ｎ，１≦θ≦Ｍ）を取得する（Ｐ００１）。続いて、特徴抽出部１１３Ａは、信号マップＢｔ（ｄ，θ）の絶対値を算出することによって、絶対値｜Ｂｔ（ｄ，θ）｜を取得する（Ｐ００２）。｜Ｂｔ（ｄ，θ）｜は、２次元行列に該当する。そして、特徴抽出部１１３Ａは、２次元行列｜Ｂｔ（ｄ，θ）｜に対して畳込みフィルタの適用による畳込み処理（Ｐ００３）とｄ軸およびθ軸の各方向の次元圧縮（Ｐ００４）とを複数回行う。
【００５３】
次元圧縮（Ｐ００４）においては、ＭａｘＰｏｏｌｉｎｇなどのプーリングが行われてもよい。また、次元圧縮は、畳込み処理のストライド幅を２以上の整数とすることによって行われてもよい。参考文献２に記載された例では、ｄ軸およびθ軸の各方向のサイズが１／２となるように次元圧縮が行われているが、次元圧縮の度合いは限定されない。
【００５４】
また、図４に示した例では、畳込み処理（Ｐ００３）と次元圧縮（Ｐ００４）とが交互にそれぞれ複数回行われている。しかし、畳込み処理および次元圧縮それぞれが複数回行われればよく、複数回の畳込み処理が連続的に行われてもよい。参考文献２に記載された例では、ｄ軸およびθ軸の各方向のサイズが１／３２となるように、サイズを１／２にする次元圧縮が５回行われているが、畳込み処理および次元圧縮それぞれの回数も限定されない。
【００５５】
上記したようにして、畳込み処理および次元圧縮それぞれが複数回行われることによって、複素数（テンソル要素）の特徴が抽出される（ここでは、複素数（テンソル要素）の特徴が実数として抽出される）。以下では、信号マップに基づいて抽出された特徴を「特徴マップ」とも言う。ここで、特徴マップのテンソルの次元を、（ｌ，ｎ，ｍ）と表すことにする。ｌは、特徴マップの次元数（テンソル要素の特徴数）を示し、ｎは、特徴マップの距離方向の次元数を示し、ｍは、特徴マップの角度方向の次元数を示す。
【００５６】
このとき、識別可能なクラスの数がｎｃｌであり、特徴マップの次元数ｌとクラス数ｎｃｌとが異なる場合が存在し得る。かかる場合には、特徴抽出部１１３Ａは、特徴マップの次元数を変更する（Ｐ００５）。具体的に、特徴抽出部１１３Ａは、特徴マップの次元数ｌをクラス数ｎｃｌに合わせることによって、特徴マップの次元数を変更する。これによって、特徴マップのテンソルの次元（ｌ，ｎ，ｍ）が（ｎｃｌ，ｎ，ｍ）に変更される。
【００５７】
このような特徴マップの次元数の変更は、１回の畳込み処理によって実現されてもよい。あるいは、特徴マップの次元数の変更は、出力次元を（ｋ，ｎ，ｍ）（ただし、ｍ≠ｎｃｌ）とする畳込み処理が複数回行われた後に、出力次元を（ｎｃｌ，ｎ，ｍ）とする畳込み処理が行われることによって実現されてもよい。
【００５８】
続いて、拡張部１１４は、特徴マップの距離ｄおよび角度θそれぞれの軸方向の次元を復元させる（Ｐ００６）。具体的に、拡張部１１４は、テンソルの次元が（ｎｃｌ，ｎ，ｍ）である特徴マップに基づいて、特徴マップの距離方向および角度方向の次元（テンソル要素の特徴の次元）を逆畳込み処理（Ｄｅｃｏｎｖｏｌｕｔｉｏｎ）によって、Ｐ００１において取得された信号マップ（多次元テンソル）の次元まで拡張する。これによって、次元再拡張後の特徴マップ（再拡張多次元テンソル）が生成され、特徴マップのテンソルの次元は、（ｎｃｌ，ｎ，ｍ）から（ｎｃｌ，Ｎ，Ｍ）に変更される。
【００５９】
続いて、領域分割部１１５は、テンソルの次元が（ｎｃｌ，Ｎ，Ｍ）に変更された特徴マップのテンソル要素ごとにクラスを識別することによって領域分割を行う。まず、領域分割部１１５は、テンソルの次元が（ｎｃｌ，Ｎ，Ｍ）に変更された特徴マップの各テンソル要素をＳｏｆｔｍａｘ関数などによって０以上１以下の値に変換する。これによって、特徴マップの各次元に対応して各（ｄ，θ）における識別確率が算出される（Ｐ００７）。なお、Ｓｏｆｔｍａｘ関数の詳細については、以下の参考文献３を参照されたい。
【００６０】
参考文献３：斎藤  康毅，"ゼロから作るDeep Learning-Pythonで学ぶディープラーニングの理論と実装"，O'Reill出版，pp66-71
【００６１】
続いて、領域分割部１１５は、特徴マップの各次元に対応して存在する各（ｄ，θ）における識別確率に基づいて、（ｄ，θ）ごとに識別確率が最大となる特徴マップの次元を検出する。特徴マップの各次元に対しては、あらかじめクラスが対応付けられている。そこで、領域分割部１１５は、（ｄ，θ）ごとに識別確率が最大となる特徴マップの次元を検出すと、検出した特徴マップの次元に対応するクラスを（ｄ，θ）ごとに識別し、識別結果として出力する（Ｐ００８）。
【００６２】
上記したように、クラスの種類は特に限定されない。ここでは、クラスの種類として、物体に対応するクラスと背景に対応するクラスを主に想定する。背景は、対象とする物体が存在しない領域であり、偽像が発生している領域も含み得る。また、物体の種類は、１つであってもよいし、複数であってもよい。物体の種類が複数である場合、複数の種類それぞれの物体に対応するクラスが存在してよい。
【００６３】
例えば、クラスの種類として、物体に対応するクラスと背景に対応するクラスとを想定する場合には、物体に対応するクラス１つと背景に対応するクラス１つとの合計２つがクラス数ｎｃｌとなる。さらに、物体の種類として、人、車、および、自転車の３種類を想定する場合には、物体に対応するクラス３つと背景に対応するクラス１つとの合計４つが、クラス数ｎｃｌとなる。
【００６４】
以上、信号マップの各座標（ｄ，θ）（ただし、１≦ｄ≦Ｎ，１≦θ≦Ｍ）におけるクラス識別の処理例について詳細に説明した。このような推論時におけるクラス識別がなされる前には、畳込みフィルタのパラメータ（例えば、重み、バイアスなど）などがあらかじめ学習される必要がある。すなわち、特徴抽出部１１３Ａは、学習処理によって得られたパラメータを畳込み処理に利用する。このような学習処理には、信号マップの座標に対してクラスの対応付け（ラベル付け）がなされた教師データが必要となる。
【００６５】
図５は、学習処理に利用される信号マップの絶対値が画像化された一例を示す図である。図５を参照すると、画像７０−２が示されている。図３に示した画像７０−１と同様に、画像７０−２の縦軸は、レーダ中心方向からの角度に相当し、画像７０−２の横軸は、レーダからの距離に相当する。また、画像７０−２において、物体からの反射波の強度が大きい領域ほど、輝度が高くなっている。一方、物体からの反射波の強度が小さい領域ほど、輝度が低くなっている。
【００６６】
より具体的に、画像７０−２において、物体からの反射波の強度の最も大きい領域７１−２は、実際に物体が存在する領域である。また、画像７０−２において、物体からの反射波の強度が次に大きい領域７１−３も、実際に物体が存在する領域である。一方、画像７０−２において、物体からの反射波の強度が次に大きい領域７２−４は、偽像が発生した領域である。
【００６７】
図６は、図５に示した画像７０−２に対応する教師データの例を示す図である。図６を参照すると、教師データ８０が示されている。図５を参照しながら説明したように、領域７１−２と領域７１−３とは、実際に物体が存在する領域である。したがって、領域７１−２および領域７１−３それぞれの座標には、物体に対応するクラスが対応付けられるのがよい。図６に示した教師データ８０を参照すると、領域８１−１および領域８１−２それぞれに対して、物体に対応するクラスのラベルが付されている。
【００６８】
一方、図５を参照しながら説明したように、領域７２−４は、偽像が発生した領域である。したがって、領域７２−４の座標には、物体に対応するクラスが対応付けられなくてよい。図６に示した教師データ８０を参照すると、偽像が発生した領域７２−４（図５）に対応する領域に対しては、物体に対応するクラスのラベルが付されていない。
【００６９】
このように、信号マップには偽像が発生した領域が存在し得る。そのため、教師データの生成は、人手を用いたとしても困難な場合もあり得る。そこで、教師データの生成には、以下の参考文献４に記載された手法が適用されてもよい。これによって、精度の高い教師データが生成され得る。
【００７０】
参考文献４：特願２０１７−０６１４２０の出願明細書
【００７１】
さらに、教師データが生成されるに際しては、信号マップ生成時に下記の参考文献５に記載されているような角度方向の高分解能化処理を適用しておいてもよい。そうすれば、高分解能な教師データが生成され得る。このとき、特徴抽出部１１３Ａは、高分解能化処理に基づいて生成された教師データを用いた学習処理によって得られたパラメータを畳込み処理に利用してよい。これによって、より高精度なクラス識別が可能となる。
【００７２】
参考文献５：Y.Wakamatsu, H.
Yamada, Y.Yamaguchi, "MIMO Doppler radar using Khatri-Rao product virtual
array for human location estimation", in Electromagnetics of the 2014 IEEE
International Workshop on, 2014.
【００７３】
なお、高分解能化処理が適用された場合には、高分解能化処理に相当する分の演算量が増大してしまう。しかし、高分解能な教師データを用いた学習処理によって得られたパラメータを推論時に利用するようにすれば、推論時において高分解能化処理は適用されなくとも、クラスの識別結果を高分解能化することが可能となる。
【００７４】
上記した学習処理は、教師データと上記したＰ００１〜Ｐ００８が実行された結果として出力される識別結果（座標ごとに識別されたクラス）とに基づいて行われる。例えば、学習処理は、ＳｏｆｔｍａｘＣｒｏｓｓＥｎｔｒｏｐｙを用いて、教師データと識別結果との誤差を逆伝播することによって行われてよい。ＳｏｆｔｍａｘＣｒｏｓｓＥｎｔｒｏｐｙの詳細は、上記した参考文献３を参照されたい。
【００７５】
（１．３．効果）
第１の実施形態によれば、信号マップの振幅値（信号強度）に対する畳込み処理によって、高精度にクラス識別が可能となる。推論時における畳込み処理に利用される畳込みフィルタには、信号マップの信号強度の分布パターンの学習処理によって得られたパラメータが利用され得る。これによって、信号マップの振幅値（信号強度）の閾値があらかじめ設定されなくても、高精度にクラス識別が可能となる。したがって、クラスを識別する機能を有する識別器（識別装置）の設計が容易となる。
【００７６】
以上、第１の実施形態について説明した。
【００７７】
（２．第２の実施形態）
続いて、第２の実施形態について説明する。
【００７８】
（２．１．構成）
まず、第２の実施形態に係る識別装置１０Ｂの構成例について説明する。
【００７９】
図７は、第２の実施形態に係る識別装置１０Ｂの構成例を示す図である。図２および図７に示すように、第２の実施形態に係る識別装置１０Ｂは、第１の実施形態に係る識別装置１０Ａと比較して、特徴抽出部１１３Ａの代わりに、特徴抽出部１１３Ｂを有する点において異なっている。したがって、第２の実施形態に係る特徴抽出部１１３Ｂについて主に説明し、他の構成についての詳細な説明は省略する。
【００８０】
第２の実施形態においても、特徴抽出部１１３Ｂは、複素数の多次元テンソルに基づく畳込み処理および次元圧縮それぞれを複数回行ってテンソル要素の特徴を抽出する。特に、第２の実施形態においては、特徴抽出部１１３Ｂが、複素数の多次元テンソルの各テンソル要素自体に対して畳込み処理および次元圧縮それぞれを複数回行うことによって複素数によるテンソル要素の特徴を抽出し、複素数によるテンソル要素の特徴を実数によるテンソル要素の特徴に変換する場合を説明する。
【００８１】
具体的に、第１の実施形態においては、複素数の信号マップの絶対値に対する畳込み処理によってクラス識別を行う場合を説明した。信号マップの絶対値は、信号マップの振幅値（信号強度）に相当する。しかし、移動物体からの反射波がレーダによって受信される場合には、反射波の位相が変化することが想定される。そこで、信号マップの位相には、各座標（ｄ，θ）（ただし、１≦ｄ≦Ｎ，１≦θ≦Ｍ）の間を埋める角度または距離に関する情報が含まれていると考えられる。
【００８２】
したがって、信号マップの絶対値だけではなく、信号マップの位相も考慮したほうが物体の存否をより高精度に識別可能になることが想定される。そこで、第２の実施形態においては、信号マップの絶対値だけではなく、信号マップの位相も考慮してクラス識別を行うことによって、さらに高精度に物体を識別することを可能とする技術を説明する。
【００８３】
以上、第２の実施形態に係る識別装置１０Ｂの構成例について説明した。
【００８４】
（２．２．動作）
続いて、第２の実施形態に係る識別装置１０Ｂの動作例について説明する。図８は、第２の実施形態に係るクラス識別の処理例を示すフローチャートである。第１の実施形態と同様にして、特徴抽出部１１３Ｂは、取得部１１２によって得られた複素数の信号マップＢｔ（ｄ，θ）（ただし、１≦ｄ≦Ｎ，１≦θ≦Ｍ）を取得する（Ｐ１０１）。続いて、特徴抽出部１１３Ｂは、複素数の信号マップＢｔ（ｄ，θ）に対して複素畳込みフィルタの適用による畳込み処理（Ｐ１０２）とｄ軸およびθ軸の各方向の次元圧縮（Ｐ１０３）とを複数回行う。
【００８５】
次元圧縮（Ｐ１０３）は、畳込み処理のストライド幅を２以上の整数とすることによって行われてもよい。参考文献２に記載された例では、ｄ軸およびθ軸の各方向のサイズが１／２となるように次元圧縮が行われているが、次元圧縮の度合いは限定されない。
【００８６】
また、図８に示した例では、畳込み処理（Ｐ１０２）と次元圧縮（Ｐ１０３）とが交互にそれぞれ複数回行われている。しかし、畳込み処理および次元圧縮それぞれが複数回行われればよく、複数回の畳込み処理が連続的に行われてもよい。畳込み処理および次元圧縮それぞれの回数も限定されない。
【００８７】
上記したようにして、畳込み処理および次元圧縮それぞれが複数回行われることによって、複素数（テンソル要素）の特徴が抽出される。ここでは、複素数（テンソル要素）の特徴が複素数として抽出される。すなわち、複素数からなる特徴マップ（以下、「複素数特徴マップ」とも言う。）が抽出される。特徴抽出部１１３Ｂは、複素数特徴マップを実数からなる特徴マップ（以下、「実数特徴マップ」とも言う。）に変換する（Ｐ１０４）。
【００８８】
例えば、特徴抽出部１１３Ｂは、複素数特徴マップの実部と虚部とを分けることによって、複素数特徴マップを実数特徴マップに変換する。例えば、複素数特徴マップのテンソルの次元が（ｋ，ｄａ，θａ）である場合、複素数特徴マップの実部と虚部とを分けられることによって、テンソルの次元が（２ｋ，ｄａ，θａ）である実数特徴マップが得られる。なお、複素数特徴マップから実数特徴マップへの変換の処理としては、下記の参考文献５に記載された手法を利用することができる。
【００８９】
参考文献５：特願２０１７−０４４８４６の出願明細書
【００９０】
続いて、図８に示したＰ１０５〜Ｐ１１０が実行される。しかし、図８に示したＰ１０５〜Ｐ１１０は、図４に示したＰ００３〜Ｐ００８と同じである。したがって、図８に示したＰ１０５〜Ｐ１１０についての詳細な説明は省略する。
【００９１】
以上、信号マップの各座標（ｄ，θ）（ただし、１≦ｄ≦Ｎ，１≦θ≦Ｍ）におけるクラス識別の処理例について詳細に説明した。第２の実施形態においても、第１の実施形態と同様にして、推論時におけるクラス識別がなされる前には、畳込みフィルタのパラメータ（例えば、重み、バイアスなど）などがあらかじめ学習される必要がある。
【００９２】
第２の実施形態に係る学習処理は、教師データと上記したＰ１０１〜Ｐ１１０が実行された結果として出力される識別結果（座標ごとに識別されたクラス）とに基づいて行われる。このとき、第２の実施形態においても、第１の実施形態と同様に、ＳｏｆｔｍａｘＣｒｏｓｓＥｎｔｒｏｐｙを用いて、教師データと識別結果との誤差を逆伝播することによって、学習処理が行われてよい。その他、第２の実施形態に係る学習処理（例えば、教師データの生成など）は、第１の実施形態と同様に行われてよい。
【００９３】
（２．３．効果）
第２の実施形態によれば、複素数の信号マップＢｔ（ｄ，θ）それ自体に対する複素畳込み処理によって、複素数特徴マップを得た後、複素数特徴マップを実数特徴マップに変換し（複素数特徴マップを複素数から実数への変換層に通す）、実数特徴マップに基づいてクラス識別を行う。これによって、信号マップの絶対値だけではなく、信号マップの位相も考慮されてクラス識別がされるため、物体の存否をより高精度に識別可能となり、さらに高精度に物体を識別することが可能となる。
【００９４】
推論時における複素畳込み処理（Ｐ１０２）に利用される複素畳込みフィルタには、信号マップの位相の分布パターンの学習処理によって得られたパラメータが利用され得る。これによって、学習処理において、信号マップの振幅値（信号強度）の誤差だけではなく、信号マップの位相の誤差も、複素畳込みフィルタのパラメータに反映され、さらに高精度に物体を識別することが可能となる。また、振幅雑音に対して堅牢となることが期待される。
【００９５】
以上、第２の実施形態について説明した。
【００９６】
（３．第３の実施形態）
続いて、第３の実施形態について説明する。ここで、第１の実施形態および第２の実施形態においては、１つのサンプリング時刻ｔにおける信号マップＢｔ（ｄ，θ）（ただし、１≦ｄ≦Ｎ，１≦θ≦Ｍ）が用いられる例を主に説明した（Ｐ００１，Ｐ１０１）。すなわち、第１の実施形態および第２の実施形態においては、信号マップのテンソルの次元は、（１，Ｎ，Ｍ）であった。
【００９７】
しかし、１つのサンプリング時刻ｔにおける信号マップが用いられる場合には、信号マップに時間的な変化が含まれないため、信号マップを時系列な情報として扱うことができない。信号マップを時系列な情報として扱うことができない場合、移動物体と静止物体との区別を行うのが困難となり得る。そのため、移動物体を物体として認識し、静止物体を背景として認識したい場合などには（移動物体のみを物体として認識したい場合などには）、信号マップから静止物体をあらかじめ除去しておく必要が生じる。
【００９８】
そこで、第３の実施形態においては、識別対象の時刻ｔの前後数時間分の信号マップを用いる。これによって、信号マップに時間的な変化が含まれるようになるため（信号マップに含まれる時間的な変化も学習処理に利用されるようになるため）、移動物体と静止物体との区別を行うのが可能となる。そのため、移動物体を物体として認識し、静止物体を背景として認識したい場合であっても、信号マップから静止物体をあらかじめ除去する必要がなくなる。
【００９９】
具体的に、第３の実施形態においては、第１の実施形態および第２の実施形態において利用される信号マップ（具体的には、Ｐ００１およびＰ１０１において取得される信号マップ）を、下記の（数式３）に示すような信号マップＢｔ+ｋ（ｄ，θ）に変更すればよい。なお、第１の実施形態および第２の実施形態における信号マップＢｔ（ｄ，θ）が、（数式３）に示す信号マップＢｔ+ｋ（ｄ，θ）に変更された場合、信号マップのテンソルの次元は、（１，Ｎ，Ｍ）から（２ｍ+１，Ｎ，Ｍ）に変更される。
【０１００】
【数２】
2019074478.tif 000004
【０１０１】
以上、第３の実施形態について説明した。
【０１０２】
（４．ハードウェア構成例）
続いて、本発明の実施形態に係る識別装置１０のハードウェア構成例について説明する。以下では、本発明の実施形態に係る識別装置１０のハードウェア構成例として、情報処理装置９００のハードウェア構成例について説明する。なお、以下に説明する情報処理装置９００のハードウェア構成例は、識別装置１０のハードウェア構成の一例に過ぎない。したがって、識別装置１０のハードウェア構成は、以下に説明する情報処理装置９００のハードウェア構成から不要な構成が削除されてもよいし、新たな構成が追加されてもよい。
【０１０３】
図９は、本発明の実施形態に係る識別装置１０の例としての情報処理装置９００のハードウェア構成を示す図である。情報処理装置９００は、ＣＰＵ（Ｃｅｎｔｒａｌ  Ｐｒｏｃｅｓｓｉｎｇ  Ｕｎｉｔ）９０１と、ＲＯＭ（Ｒｅａｄ  Ｏｎｌｙ  Ｍｅｍｏｒｙ）９０２と、ＲＡＭ（Ｒａｎｄｏｍ  Ａｃｃｅｓｓ  Ｍｅｍｏｒｙ）９０３と、ホストバス９０４と、ブリッジ９０５と、外部バス９０６と、インタフェース９０７と、入力装置９０８と、出力装置９０９と、ストレージ装置９１０と、通信装置９１１と、を備える。
【０１０４】
ＣＰＵ９０１は、演算処理装置および制御装置として機能し、各種プログラムに従って情報処理装置９００内の動作全般を制御する。また、ＣＰＵ９０１は、マイクロプロセッサであってもよい。ＲＯＭ９０２は、ＣＰＵ９０１が使用するプログラムや演算パラメータ等を記憶する。ＲＡＭ９０３は、ＣＰＵ９０１の実行において使用するプログラムや、その実行において適宜変化するパラメータ等を一時記憶する。これらはＣＰＵバス等から構成されるホストバス９０４により相互に接続されている。
【０１０５】
ホストバス９０４は、ブリッジ９０５を介して、ＰＣＩ（Ｐｅｒｉｐｈｅｒａｌ  Ｃｏｍｐｏｎｅｎｔ  Ｉｎｔｅｒｃｏｎｎｅｃｔ／Ｉｎｔｅｒｆａｃｅ）バス等の外部バス９０６に接続されている。なお、必ずしもホストバス９０４、ブリッジ９０５および外部バス９０６を分離構成する必要はなく、１つのバスにこれらの機能を実装してもよい。
【０１０６】
入力装置９０８は、マウス、キーボード、タッチパネル、ボタン、マイクロフォン、スイッチおよびレバー等ユーザが情報を入力するための入力手段と、ユーザによる入力に基づいて入力信号を生成し、ＣＰＵ９０１に出力する入力制御回路等から構成されている。情報処理装置９００を操作するユーザは、この入力装置９０８を操作することにより、情報処理装置９００に対して各種のデータを入力したり処理動作を指示したりすることができる。
【０１０７】
出力装置９０９は、例えば、ＣＲＴ（Ｃａｔｈｏｄｅ  Ｒａｙ  Ｔｕｂｅ）ディスプレイ装置、液晶ディスプレイ（ＬＣＤ）装置、ＯＬＥＤ（Ｏｒｇａｎｉｃ  Ｌｉｇｈｔ  Ｅｍｉｔｔｉｎｇ  Ｄｉｏｄｅ）装置、ランプ等の表示装置およびスピーカ等の音声出力装置を含む。
【０１０８】
ストレージ装置９１０は、データ格納用の装置である。ストレージ装置９１０は、記憶媒体、記憶媒体にデータを記録する記録装置、記憶媒体からデータを読み出す読出し装置および記憶媒体に記録されたデータを削除する削除装置等を含んでもよい。ストレージ装置９１０は、例えば、ＨＤＤ（Ｈａｒｄ  Ｄｉｓｋ  Ｄｒｉｖｅ）で構成される。このストレージ装置９１０は、ハードディスクを駆動し、ＣＰＵ９０１が実行するプログラムや各種データを格納する。
【０１０９】
通信装置９１１は、例えば、ネットワークに接続するための通信デバイス等で構成された通信インタフェースである。また、通信装置９１１は、無線通信または有線通信のどちらに対応してもよい。
【０１１０】
以上、本発明の実施形態に係る識別装置１０のハードウェア構成例について説明した。
【０１１１】
（５．まとめ）
以上に説明したように、本実施形態によれば、センサによる計測に基づいて複素数の多次元テンソルを取得する取得部と、前記複素数の多次元テンソルに基づく畳込み処理および次元圧縮それぞれを複数回行ってテンソル要素の特徴を抽出する特徴抽出部と、前記テンソル要素の特徴の次元を逆畳込み処理によって前記多次元テンソルの次元まで拡張して再拡張多次元テンソルを生成する拡張部と、前記再拡張多次元テンソルのテンソル要素ごとにクラスを識別することによって領域分割を行う領域分割部と、を備える、識別装置が提供される。
【０１１２】
かかる構成によれば、あらかじめ閾値を設定せずとも、高精度に物体を識別することが可能となる。
【０１１３】
以上、添付図面を参照しながら本発明の好適な実施形態について詳細に説明したが、本発明はかかる例に限定されない。本発明の属する技術の分野における通常の知識を有する者であれば、特許請求の範囲に記載された技術的思想の範疇内において、各種の変更例または修正例に想到し得ることは明らかであり、これらについても、当然に本発明の技術的範囲に属するものと了解される。
【符号の説明】
【０１１４】
１０（１０Ａ，１０Ｂ）  識別装置
１１１  送受信部
１１２  取得部
１１３Ａ，１１３Ｂ  特徴抽出部
１１４  拡張部
１１５  領域分割部
【図１】
2019074478.tif 000005
【図２】
2019074478.tif 000006
【図３】
2019074478.tif 000007
【図４】
2019074478.tif 000008
【図５】
2019074478.tif 000009
【図６】
2019074478.tif 000010
【図７】
2019074478.tif 000011
【図８】
2019074478.tif 000012
【図９】
2019074478.tif 000013