WO2020090540A1 - 無線通信の伝搬損失を利用した混雑度の推定方法、推定装置および推定プログラム - Google Patents
無線通信の伝搬損失を利用した混雑度の推定方法、推定装置および推定プログラム Download PDFInfo
- Publication number
- WO2020090540A1 WO2020090540A1 PCT/JP2019/041290 JP2019041290W WO2020090540A1 WO 2020090540 A1 WO2020090540 A1 WO 2020090540A1 JP 2019041290 W JP2019041290 W JP 2019041290W WO 2020090540 A1 WO2020090540 A1 WO 2020090540A1
- Authority
- WO
- WIPO (PCT)
- Prior art keywords
- propagation loss
- receiver
- transmitter
- processing
- congestion
- Prior art date
Links
Images
Classifications
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04B—TRANSMISSION
- H04B17/00—Monitoring; Testing
- H04B17/20—Monitoring; Testing of receivers
- H04B17/27—Monitoring; Testing of receivers for locating or positioning the transmitter
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S11/00—Systems for determining distance or velocity not using reflection or reradiation
- G01S11/02—Systems for determining distance or velocity not using reflection or reradiation using radio waves
- G01S11/06—Systems for determining distance or velocity not using reflection or reradiation using radio waves using intensity measurements
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
- G01S13/003—Bistatic radar systems; Multistatic radar systems
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/003—Transmission of data between radar, sonar or lidar systems and remote stations
- G01S7/006—Transmission of data between radar, sonar or lidar systems and remote stations using shared front-end circuitry, e.g. antennas
-
- G—PHYSICS
- G08—SIGNALLING
- G08G—TRAFFIC CONTROL SYSTEMS
- G08G1/00—Traffic control systems for road vehicles
- G08G1/01—Detecting movement of traffic to be counted or controlled
- G08G1/0104—Measuring and analyzing of parameters relative to traffic conditions
- G08G1/0108—Measuring and analyzing of parameters relative to traffic conditions based on the source of data
-
- G—PHYSICS
- G08—SIGNALLING
- G08G—TRAFFIC CONTROL SYSTEMS
- G08G1/00—Traffic control systems for road vehicles
- G08G1/01—Detecting movement of traffic to be counted or controlled
- G08G1/0104—Measuring and analyzing of parameters relative to traffic conditions
- G08G1/0108—Measuring and analyzing of parameters relative to traffic conditions based on the source of data
- G08G1/012—Measuring and analyzing of parameters relative to traffic conditions based on the source of data from other sources than vehicle or roadside beacons, e.g. mobile networks
-
- G—PHYSICS
- G08—SIGNALLING
- G08G—TRAFFIC CONTROL SYSTEMS
- G08G1/00—Traffic control systems for road vehicles
- G08G1/01—Detecting movement of traffic to be counted or controlled
- G08G1/0104—Measuring and analyzing of parameters relative to traffic conditions
- G08G1/0125—Traffic data processing
- G08G1/0133—Traffic data processing for classifying traffic situation
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04B—TRANSMISSION
- H04B17/00—Monitoring; Testing
- H04B17/30—Monitoring; Testing of propagation channels
- H04B17/309—Measuring or estimating channel quality parameters
- H04B17/318—Received signal strength
Abstract
受信電力取得処理と、受信電力から前記送信機と前記受信機との間の伝搬損失を計算する伝搬損失計算処理と、所定時間に前記伝搬損失計算処理により計算された複数の伝搬損失の中央値を求め、当該中央値を以降の処理に用いる伝搬損失として出力する中央値処理と、前記中央値処理により出力される伝搬損失の最大値を取得する損失ピーク取得処理と、前記最大値と、アンテナの高さと、前記無線通信の波長と、に少なくとも基づいて、前記送信機または前記受信機と遮蔽物との距離を計算する距離算出処理と、前記距離に基づいて、遮蔽物の個数を時系列に算出し、前記個数を所定時間だけ積分することにより、単位面積当たりの遮蔽物の個数を求め、混雑度を推定する混雑度計算処理とを実行する。
Description
本発明は、無線通信の伝搬損失を取得することで測定環境における人や車などの混雑の度合いを推定する技術に関する。
従来、人や車などの多い混雑した環境において、混雑の度合い(以降、混雑度と称する)を定量的に評価するために、カメラで撮影した画像から人や車などの数を測定することが行われている。
中村光貴,佐々木元晴,北直樹,鬼沢武,"見通し線上の複数人体による人体遮蔽時の伝搬損失特性",総合大会,B-1-37,2016.
ところが、カメラの高所への設置が困難な場所では、カメラの前が人や車などで遮られて奥の様子が確認できないという問題がある。また、プライバシー配慮の点でカメラの設置が困難な場所では、カメラで撮影した画像による確認ができないため、人や車の混雑度を測定することが難しい。
一方、非特許文献1では、6GHz以上の高周波数帯を利用した無線通信システムにおいて、送信機と受信機との間の通信領域が複数の人体によって遮蔽されたときの人体遮蔽損失の推定モデルが開示されているが、人や車の混雑度の測定に応用する技術は検討されていなかった。
本発明は、対象エリアにおける無線通信の伝搬損失を取得することでエリア内の人や車などの混雑度を推定することができる推定方法、推定装置および推定プログラムを提供することを目的とする。
第1の発明は、無線通信を行う送信機と受信機との間の通信領域における遮蔽物の混雑度を推定する推定方法であって、受信電力を取得する受信電力取得処理と、前記受信電力取得処理により取得された受信電力から前記送信機と前記受信機との間の伝搬損失を計算する伝搬損失計算処理と、予め決められた所定時間に前記伝搬損失計算処理により計算された複数の伝搬損失の中央値を求め、当該中央値を以降の処理に用いる伝搬損失として出力する中央値処理と、前記中央値処理により出力される伝搬損失の最大値を取得する損失ピーク取得処理と、前記損失ピーク取得処理により取得された伝搬損失の前記最大値と、前記送信機および前記受信機のアンテナの高さと、前記無線通信の波長と、に少なくとも基づいて、前記送信機または前記受信機と遮蔽物との距離を計算する距離算出処理と、前記距離算出処理により計算された前記距離に基づいて、遮蔽物の個数を時系列に算出し、前記個数を予め決められた所定時間だけ積分することにより、単位面積当たりの遮蔽物の個数を求め、混雑度を推定する混雑度計算処理とを実行することを特徴とする。
第2の発明は、無線通信を行う送信機と受信機との間の通信領域における遮蔽物の混雑度を推定する推定方法であって、受信電力を取得する受信電力取得処理と、前記受信電力取得処理により取得された受信電力から前記送信機と前記受信機との間の伝搬損失を計算する伝搬損失計算処理と、予め決められた所定時間に前記伝搬損失計算処理により計算された複数の伝搬損失の中央値を求め、当該中央値を以降の処理に用いる伝搬損失として出力する中央値処理と、前記中央値処理の出力結果から伝搬損失の分布を求め、予め保持している単位面積あたりの遮蔽物の数に対応する複数の分布データと突合せを行い、一番近い伝搬損失の分布データを選定する分布突合処理と、分布突合処理により選定された伝搬損失の分布データから当該分布データに対応付けられた単位面積あたりの遮蔽物の数を求め、前記送信機と前記受信機との間における遮蔽物の混雑度として出力する混雑度出力処理とを実行することを特徴とする。
第3の発明は、無線通信を行う送信機と受信機との間の通信領域における遮蔽物の混雑度を推定する推定装置において、受信電力を取得する受信電力取得部と、前記受信電力取得部により取得された受信電力から前記送信機と前記受信機との間の伝搬損失を計算する伝搬損失計算部と、予め決められた所定時間に前記伝搬損失計算部により計算された複数の伝搬損失の中央値を求め、当該中央値を以降の処理に用いる伝搬損失として出力する中央値処理部と、前記中央値処理部により出力された伝搬損失の最大値を取得する損失ピーク取得部と、前記損失ピーク取得部により取得された伝搬損失の前記最大値と、前記送信機および前記受信機のアンテナの高さと、前記無線通信の波長と、に少なくとも基づいて、前記送信機または前記受信機と遮蔽物との距離を計算する距離算出部と、前記距離算出部により計算された前記距離に基づいて、遮蔽物の個数を時系列に算出し、前記個数を予め決められた所定時間だけ積分することにより、単位面積当たりの遮蔽物の個数を求め、混雑度を推定する混雑度計算部とを有することを特徴とする。
第4の発明は、無線通信を行う送信機と受信機との間の通信領域における遮蔽物の混雑度を推定する推定装置において、受信電力を取得する受信電力取得部と、前記受信電力取得部により取得された受信電力から前記送信機と前記受信機との間の伝搬損失を計算する伝搬損失計算部と、予め決められた所定時間に前記伝搬損失計算部により計算された複数の伝搬損失の中央値を求め、当該中央値を以降の処理に用いる伝搬損失として出力する中央値処理部と、前記中央値処理部の出力結果から伝搬損失の分布を求め、予め保持している単位面積あたりの遮蔽物の数に対応する複数の分布データと突合せを行い、一番近い伝搬損失の分布データを選定する分布突合部と、分布突合部により選定された伝搬損失の分布データから当該分布データに対応付けられた単位面積あたりの遮蔽物の数を求め、前記送信機と前記受信機との間における遮蔽物の混雑度として出力する混雑度出力部とを有することを特徴とする。
第5の発明は、第3の発明または第4の発明の推定装置で行う処理をコンピュータに実行させるプログラムであることを特徴とする。
本発明に係る混雑度の推定方法、推定装置および推定プログラムは、対象エリアにおける無線通信の伝搬損失を取得することでエリア内の人や車などの混雑度を推定することができる。
以下、図面を参照して本発明に係る混雑度の推定方法、推定装置および推定プログラムの実施形態について説明する。ここで、本発明に係る混雑度の推定方法、推定装置および推定プログラムは、カメラの設置が困難な高所や個人情報の観点などから画像認識による遮蔽物の検出を行えない場所であっても、既存の無線通信システムを利用することで、人や車の混雑度を測定することが可能となる。
図1は、送信機Txと受信機Rxとの間の通信領域が複数の人体によって遮蔽されたときの人体遮蔽損失の推定モデルの一例を示す。図1において、送信機Txから送信される無線信号は、受信機Rxにより受信される。そして、本実施形態では、受信機Rxに接続される推定装置101が受信機Rxの受信電力から伝搬損失を求めて、送信機Txと受信機Rxとの間を対象エリアとする人や車などの遮蔽物の混雑度を推定する。ここで、本実施形態では、遮蔽物の種類が人や車など予め分かっていることを前提とし、様々な遮蔽物が混在している場合は想定していない。従って、本実施形態に係る推定装置101は、例えば、イベント会場での人の混雑度や駐車場などでの車の混雑度の推定などに利用することができる。なお、図1では、推定装置101が独立した装置として受信機Rxに接続されているが、受信機Rxに内蔵されていてもよい。また、推定装置101は、送信機Txと連携して、送信機Txの送信電力や変調方式などを制御するようにして、既知の信号を送信機Txから送信させるようにしてもよい。
図1の例では、送信機Txと受信機Rxとの間における複数の人の往来により伝搬損失が変化し、受信機Rxでの受信電力も変動する。そこで、本実施形態では、受信機Rxでの受信電力を取得し、送信機Txと受信機Rxとの間の伝搬損失を計算する。そして、伝搬損失から遮蔽物までの距離を求め、対象エリアにおける混雑度を計算する。
[第1実施形態]
図2は、第1実施形態に係る推定装置101の一例を示す。推定装置101は、受信電力取得部201、伝搬損失計算部202、中央値処理部203、損失ピーク取得部204、距離算出部205および混雑度計算部206を有する。
図2は、第1実施形態に係る推定装置101の一例を示す。推定装置101は、受信電力取得部201、伝搬損失計算部202、中央値処理部203、損失ピーク取得部204、距離算出部205および混雑度計算部206を有する。
受信電力取得部201は、送信機Txから送信される無線信号を受信して受信電力を取得する。ここで、無線信号は、無変調キャリア信号であってもよいし、変調信号であってもよい。なお、送信機Txのアンテナの高さ、送信電力、受信機Rxのアンテナの高さ、送信機Txと受信機Rxとの間の距離は、設置時に決められており、既知のパラメータとする。
伝搬損失計算部202は、受信機Rxの受信電力に基づいて、送信機Txと受信機Rxとの間の伝搬損失を計算する。ここで、送信機Txと受信機Rxとの間の伝搬損失は、遮蔽物が無い場合の受信機Rxの受信電力P0を計測して保持しておくことにより、遮蔽物がある場合の受信機Rxの受信電力PrとP0との差(P0-Pr)として計算することができる。
中央値処理部203は、予め決められた所定時間において、伝搬損失計算部202により計算された伝搬損失の中央値を求める。これは、フェージングの影響を除去するための処理であり、所定時間の間に取得した複数の伝搬損失に対する中央値を導出し、中央値の伝搬損失を用いて以降の処理を行う。ここで、所定時間は、例えば数秒間など想定されるフェージングの周期以上の時間が望ましい。
損失ピーク取得部204は、中央値処理部203が出力する所定期間の伝搬損失の中から最大値の伝搬損失(損失ピークとも称する)を取得する。これは、遮蔽物が人や車の場合、移動するので、送信機Txと受信機Rxとの間を結ぶ直線上に遮蔽物が来たときの伝搬損失を取得するためである。このように、伝搬損失の最大値を取得することにより、遮蔽物の位置を正確に推定することが可能になる。もし、伝搬損失の最大値を用いない場合、当該遮蔽物による伝搬損失が低く計算されるので、伝搬損失から遮蔽物までの距離を算出するときの誤差の原因となる。
距離算出部205は、損失ピーク取得部204が取得した伝搬損失の最大値に基づいて、遮蔽物までの距離を算出する。なお、遮蔽物までの距離の算出方法については、後で詳しく説明する。
混雑度計算部206は、距離算出部205が算出した遮蔽物までの距離のデータを予め決められた所定時間だけ蓄積して、送信機Txと受信機Rxとの間における人や車などの遮蔽物の混雑度を推定する。なお、遮蔽物の混雑度の求め方については、後で詳しく説明する。
このようにして、本実施形態に係る推定装置101は、送信機Txと受信機Rxとの間の無線通信を利用して、受信機Rxの受信電力から人や車などの遮蔽物による伝搬損失を計算し、送信機Txまたは受信機Rxと遮蔽物との距離を求めて、送信機Txと受信機Rxとの間の混雑度を推定することができる。
ここで、本実施形態に係る推定装置101は、図2に示した各ブロックを有する装置として説明したが、各ブロックが行う処理に対応するプログラムを実行するコンピュータでも実現できる。なお、プログラムは、記録媒体に記録して提供されてもよいし、ネットワークを通して提供されてもよい。
(遮蔽物までの距離の算出方法)
次に、距離算出部205が行う遮蔽物までの距離の算出方法について説明する。ここで、本実施形態では、遮蔽物を方形状のスクリーンにモデル化して計算する。
次に、距離算出部205が行う遮蔽物までの距離の算出方法について説明する。ここで、本実施形態では、遮蔽物を方形状のスクリーンにモデル化して計算する。
図3は、スクリーンの一例を示す。図3において、スクリーン150は、例えば人の場合、高さ1.4m、横幅0.4mのサイズにモデル化する。そして、送信機Txと受信機Rxとの間の遮蔽物の位置に配置されたスクリーン150の上端、右端、左端の3つのエッジを回折端とする回折損失を用いて伝搬損失を計算する。
図4は、図3に示した遮蔽モデルのパラメータの一例を示す。図4において、送信機Txのアンテナの高さがHt、受信機Rxのアンテナの高さがHr、送信機Txと受信機Rxとの間の距離がDtrである。そして、送信機Txから送信された電波は、スクリーン150の回折端(図4では上端エッジを示す)で回折されて受信機Rxで受信される。このとき、送信機Txのアンテナから回折端までの距離をdi1、回折端から受信機Rxのアンテナまでの距離をdi2とする。また、図4において、送信機Txのアンテナと受信機Rxのアンテナとの間の直線がスクリーン150と交差する位置から回折端までの高さを回折端の見通しの高さhiとする。
次に、図4で説明したパラメータを用いて、スクリーン150による遮蔽損失(伝搬損失に対応)と、スクリーン150と送信機Txまたは受信機Rxとの間の距離との関係を求める方法について説明する。
ここで、図3において、j1(ν)を上端エッジでの回折損失、j2(ν)を左端エッジでの回折損失、j3(ν)を右端エッジでの回折損失とし、回折損失の和が遮蔽損失であるとすると、遮蔽物(スクリーン150)による遮蔽損失J(ν)[dB]は、下式により求まる。
ここで、νはナイフエッジ回折モデルなどで一般に用いられている回折パラメータで、送信機Txと受信機Rxとの間の見通し線上にエッジ先端がある場合、ν=0となる。
なお、ji(ν)は、下式により、デシベルから真値に変換される。
また、上端エッジ、左端エッジ、右端エッジのそれぞれにおける回折損失は、下式で求められることが知られている(非特許文献1参照)。
ここで、i=1のときが上端エッジでの回折損失J1(ν)、i=2のときが左端エッジでの回折損失J2(ν)、i=3のときが右端エッジでの回折損失J3(ν)である。
ここで、回折パラメータνiは、式(4)で表される。
なお、λ:無線波の波長、hi:回折端の見通しの高さ(送信機Txと受信機Rxを結ぶ直線上から回折端までの高さ)、di1:送信機Txから回折端までの距離、di2は回折端から受信機Rxまでの距離である。
このようにして、式(1)から式(4)により、スクリーン150の遮蔽損失(J(ν))と、送信機Txまたは受信機Rxからの距離(di1またはdi2)との関係を求めることができる。ここで、スクリーン150の遮蔽損失が送信機Txと受信機Rxとの間の伝搬損失に相当すると考えると、伝搬損失を取得することにより、スクリーン150までの距離(位置)を推定することが可能になる。
次に、送信機Txと受信機Rxとの間の伝搬損失の取得方法について説明する。ここで、送信機Txと受信機Rxとの間の伝搬損失は、遮蔽物が無い場合の受信機Rxでの受信電力P0と遮蔽物がある場合の受信機Rxでの受信電力Prとの差(P0-Pr)として計算することができるが、フェージングなどの影響により、受信機Rxでの受信電力が変動するので、求められる伝搬損失も変動するという問題が生じる。そこで、本実施形態では、フェージングの影響を除去するために、中央値処理部203は、予め決められた所定時間の間に取得された伝搬損失の中央値を選択する処理を行い、以降の各処理では中央値として選択された伝搬損失の値を用いる。
図5は、損失ピーク取得部204の処理の一例を示す。図5において、横軸は時間(t)を示し、縦軸は中央値処理部203により選択された中央値の伝搬損失(dB)を示す。なお、図5の例では、時間に応じて伝搬損失が変化し、点線円の部分(時刻T1)で伝搬損失が最大(ピーク値)となることを示している。つまり、送信機Txと受信機Rxとの間の遮蔽物(図1の例では人)が移動していると仮定した場合、時刻T1までの時間では、遮蔽物による送信機Txと受信機Rxとの間の通信の遮蔽が弱く、そして時刻T1で送信機Txと受信機Rxとを結ぶ直線を完全に遮蔽し、時刻T1を過ぎると再び遮蔽が弱くなり、伝搬損失が小さくなっていく。
ここで、式(1)から式(4)で説明した遮蔽損失(伝搬損失)と、送信機Txまたは受信機Rxからの距離との関係は、遮蔽物が送信機Txと受信機Rxとを結ぶ直線上にあることが前提であり、遮蔽物が送信機Txと受信機Rxとを結ぶ直線上にない場合(少しずれている場合など)では、伝搬損失が小さくなるので、送信機Txまたは受信機Rxから遮蔽物までの距離を正確に推定することが難しくなる。そこで、本実施形態では、損失ピーク取得部204により、伝搬損失が最大値となる時刻を判定して、当該時刻での伝搬損失を用いて遮蔽物と送信機Txまたは受信機Rxとの間の距離を算出する。
なお、上記の説明では、式(1)から式(4)を用いて、伝搬損失から遮蔽物までの距離を算出するものとして説明したが、予め式(1)から式(4)を用いて、スクリーン150の大きさ毎に伝搬損失と遮蔽物までの距離との関係を計算したデータ表を予め作成しておいてもよい。これにより、受信機Rxで受信電力から求めた伝搬損失を中央値処理して最大となる伝搬損失を取得するだけで、予め計算したデータ表を参照して遮蔽物までの距離を推定することができる。
(遮蔽物の混雑度の求め方)
次に、混雑度の求め方について説明する。
次に、混雑度の求め方について説明する。
図6は、時系列に検出される遮蔽物の一例を示す。ここで、図6では、送信機Txまたは受信機Rxから同じ距離となる所定の範囲(図6では、長方形のボックス領域)の遮蔽物の個数を取得する。
図6(a)は、ある時刻に検出された遮蔽物(ここでは「人」)の一例で、検出された位置に一人の「人」が居ることを示す。図6(a)のボックス領域は、例えば、短辺を1m、長辺を10mとする。ここで、送信機Txまたは受信機Rxからの距離が±0.5mの範囲内で検出された遮蔽物を同じ距離の遮蔽物であると考えると、このときのボックス領域の短辺の長さは、1mとなる。また、ボックス領域の長辺の長さは、送信機Txまたは受信機Rxから同じ距離となる範囲に対応し、例えば10mのように予め設定しておく。図6(a)では、短辺1mおよび長辺10mの10m2のボックス領域に1人の「人」が検出されているので、単位面積当たり人数(混雑度)は、1÷10=0.1となる。
図6(b)は、図6(a)と同じ領域内で所定時間内に検出された「人」の一例で、三人の「人」が検出されている。ここで、例えば、図6(a)で検出された人が領域内を等速度で歩いている場合、所定時間内(例えば10秒間)に同じ人物が等間隔で複数回検出される。図6(b)の例では、同じ人物が等間隔で3回検出されている。この場合、10m2のボックス領域に3人の「人」が居るので、単位面積当たり人数(混雑度)は、3÷10=0.3となる。このように、本実施形態では、瞬間的なある時刻の所定空間の人数を数えて混雑度を求めるのではなく、所定時間だけ時間軸方向に積分した人数(同一人物を含む)から混雑度を求めている。
図6(c)は、前後する異なる距離の複数の領域(同図では6つのボックス領域)の各々の領域で検出された遮蔽物の一例を示す。図6(c)では、例えば、図6(b)と同様の方法で求めた各ボックス領域の人数を予め決められた所定時間だけ累積した人数(積分した人数)としてカウントし、当該人数の「人」が複数のボックス領域に居たと見なしている。このとき、混雑度は、下式のように求めることができる。
混雑度=(複数のボックス領域の人数の和)÷(複数のボックス領域の面積)
図6(c)の例では、6つのボックス領域の人数の和は、11人である。また、6つのボックス領域の面積は、各ボックス領域の長辺を10m、短辺を1m、とすると、10×(1×6)=60m2となる。そして、単位面積当たりの人数(混雑度)は、11÷60=約0.18となる。ここで、複数の領域は、「人」が検出されていない連続した領域を含み、混雑度を求める対象領域の全体に相当する。
混雑度=(複数のボックス領域の人数の和)÷(複数のボックス領域の面積)
図6(c)の例では、6つのボックス領域の人数の和は、11人である。また、6つのボックス領域の面積は、各ボックス領域の長辺を10m、短辺を1m、とすると、10×(1×6)=60m2となる。そして、単位面積当たりの人数(混雑度)は、11÷60=約0.18となる。ここで、複数の領域は、「人」が検出されていない連続した領域を含み、混雑度を求める対象領域の全体に相当する。
このようにして、本実施形態に係る推定装置101は、送信機Txと受信機Rxとの間にある遮蔽物の混雑度を推定することができる。
次に、上述した推定装置101における混雑度の推定方法の処理の流れについて説明する。
図7は、混雑度の推定方法のフローチャートを示す。なお、図7の処理は、例えば図2で説明した推定装置101により実行される。
ステップS101において、受信電力取得部201は、送信機Txから送信される無線信号を受信機Rxで受信したときの受信電力を取得する。
ステップS102において、伝搬損失計算部202は、送信機Txと受信機Rxとの間の伝搬損失を計算する。
ステップS103において、中央値処理部203は、予め決められた所定時間において、伝搬損失計算部202により計算された伝搬損失の中央値を求める。
ステップS104において、損失ピーク取得部204は、中央値処理部203が出力する所定期間の伝搬損失の中から最大値の伝搬損失を取得する。
ステップS105において、距離算出部205は、損失ピーク取得部204が取得した伝搬損失の最大値に基づいて、遮蔽物までの距離を算出する。
ステップS106において、混雑度計算部206は、距離算出部205が算出した遮蔽物までの距離のデータを予め決められた所定時間だけ蓄積して遮蔽物の個数を算出し、送信機Txと受信機Rxとの間における人や車などの遮蔽物の混雑度を推定する。
このようにして、本実施形態に係る混雑度の推定方法では、送信機Txと受信機Rxとの間の無線通信を利用して、受信機Rxの受信電力から人や車などの遮蔽物による伝搬損失を計算し、送信機Txまたは受信機Rxと遮蔽物との距離を求めて所定時間の個数を算出し、送信機Txと受信機Rxとの間の混雑度を推定することができる。
[第2実施形態]
次に、第1実施形態の混雑度の推定方法とは異なる方法で混雑度を推定する方法について説明する。
次に、第1実施形態の混雑度の推定方法とは異なる方法で混雑度を推定する方法について説明する。
図8は、第2実施形態に係る推定装置101aの一例を示す。推定装置101aは、受信電力取得部201、伝搬損失計算部202、中央値処理部203、分布突合部211、混雑度出力部212を有する。ここで、図2と同符号のブロックは、図2の推定装置101と同様の動作を行うブロックであり、重複する説明は省略する。
分布突合部211は、単位面積あたりの遮蔽物の数とそのときの伝搬損失との関係を示す複数の分布データを実測またはシミュレーションにより予め取得し、内部のメモリなどに保持しているものとする。そして、分布突合部211は、中央値処理部203が予め決められた所定時間において算出した出力結果から伝搬損失の分布を求め、予め保持している複数の分布データと突合せを行い、一番近い伝搬損失の分布データを選定する。
混雑度出力部212は、分布突合部211により選定された伝搬損失の分布データから当該分布データに対応付けられた単位面積あたりの遮蔽物の数がわかるので、これを送信機Txと受信機Rxとの間における人や車などの遮蔽物の混雑度として出力する。
このようにして、本実施形態に係る推定装置101aは、送信機Txと受信機Rxとの間の無線通信を利用して、受信機Rxの受信電力から人や車などの遮蔽物による伝搬損失を計算して分布データを算出し、予め保持している単位面積あたりの遮蔽物の数とそのときの伝搬損失との関係を示す複数の分布データと突き合わせることにより、送信機Txと受信機Rxとの間の混雑度を推定することができる。
ここで、本実施形態に係る推定装置101aは、図8に示した各ブロックを有する装置として説明したが、各ブロックが行う処理に対応するプログラムを実行するコンピュータでも実現できる。なお、プログラムは、記録媒体に記録して提供されてもよいし、ネットワークを通して提供されてもよい。
図9は、単位面積あたりの遮蔽物の数とそのときの伝搬損失との関係を示す複数の分布データの一例を示す。図9において、横軸は伝搬損失の確率変数[dB]、縦軸は累積確率分布[%]をそれぞれ示す。図9は、遮蔽物が人の場合の例を示し、実線は、単位面積(1m2)あたりの人数が0人、点線は、単位面積(1m2)あたりの人数が0.5人、破線は、単位面積(1m2)あたりの人数が1.0人、のときのそれぞれの分布データを示す。
図9において、中央値処理部203が計算した伝搬損失の測定結果は、○印でプロットされている。図9の例では、測定結果が予め取得され保持されている点線で示した分布データに最も近い特性であることが分かる。従って、混雑度出力部212は、点線で示した分布データに対応づけられた混雑度(0.5人/m2)である、と推定することができる。
このようにして、本実施形態に係る推定装置101aは、送信機Txと受信機Rxとの間の無線通信を利用して、受信機Rxの受信電力から人や車などの遮蔽物による伝搬損失を計算して分布データを算出し、予め取得された複数の分布データと突き合わせることにより、送信機Txと受信機Rxとの間の混雑度を推定することができる。
次に、上述した推定装置101aにおける混雑度の推定方法の処理の流れについて説明する。
図10は、混雑度の推定方法のフローチャートを示す。なお、図10の処理は、例えば図8で説明した推定装置101aにより実行される。なお、ステップS101からステップS103までの処理は、図7で説明した推定装置101のフローチャートと同じなので重複する説明は省略する。
ステップS201において、分布突合部211は、中央値処理部203が予め決められた所定時間において算出した出力結果から伝搬損失の分布を求め、予め保持している混雑度が既知の複数の分布データと突合せを行い、一番近い伝搬損失の分布データを選定する。
ステップS202において、混雑度出力部212は、ステップS201の処理により選定された伝搬損失の分布データから当該分布データに対応付けられた単位面積あたりの遮蔽物の数を判定し、これを送信機Txと受信機Rxとの間における人や車などの遮蔽物の混雑度として出力する。
このようにして、本実施形態に係る混雑度の推定方法では、送信機Txと受信機Rxとの間の無線通信を利用して、受信機Rxの受信電力から人や車などの遮蔽物による伝搬損失を計算して分布データを算出し、予め取得された複数の分布データと突き合わせることにより、送信機Txと受信機Rxとの間の混雑度を推定することができる。
以上、各実施形態で説明したように、本発明に係る混雑度の推定方法、推定装置および推定プログラムは、カメラの設置が困難な高所や個人情報の観点などから画像認識による遮蔽物の検出を行えない場所であっても、既存の無線通信システムを利用することで、人や車などの遮蔽物の混雑度を測定することができる。
推定装置・・・101,101a;201・・・受信電力取得部;202・・・伝搬損失計算部;203・・・中央値処理部;204・・・損失ピーク取得部;205・・・距離算出部;206・・・混雑度計算部;211・・・分布突合部;212・・・混雑度出力部
Claims (5)
- 無線通信を行う送信機と受信機との間の通信領域における遮蔽物の混雑度を推定する推定方法であって、
受信電力を取得する受信電力取得処理と、
前記受信電力取得処理により取得された受信電力から前記送信機と前記受信機との間の伝搬損失を計算する伝搬損失計算処理と、
予め決められた所定時間に前記伝搬損失計算処理により計算された複数の伝搬損失の中央値を求め、当該中央値を以降の処理に用いる伝搬損失として出力する中央値処理と、
前記中央値処理により出力される伝搬損失の最大値を取得する損失ピーク取得処理と、
前記損失ピーク取得処理により取得された伝搬損失の前記最大値と、前記送信機および前記受信機のアンテナの高さと、前記無線通信の波長と、に少なくとも基づいて、前記送信機または前記受信機と遮蔽物との距離を計算する距離算出処理と、
前記距離算出処理により計算された前記距離に基づいて、遮蔽物の個数を時系列に算出し、前記個数を予め決められた所定時間だけ積分することにより、単位面積当たりの遮蔽物の個数を求め、混雑度を推定する混雑度計算処理と
を実行することを特徴とする推定方法。 - 無線通信を行う送信機と受信機との間の通信領域における遮蔽物の混雑度を推定する推定方法であって、
受信電力を取得する受信電力取得処理と、
前記受信電力取得処理により取得された受信電力から前記送信機と前記受信機との間の伝搬損失を計算する伝搬損失計算処理と、
予め決められた所定時間に前記伝搬損失計算処理により計算された複数の伝搬損失の中央値を求め、当該中央値を以降の処理に用いる伝搬損失として出力する中央値処理と、
前記中央値処理の出力結果から伝搬損失の分布を求め、予め保持している単位面積あたりの遮蔽物の数に対応する複数の分布データと突合せを行い、一番近い伝搬損失の分布データを選定する分布突合処理と、
分布突合処理により選定された伝搬損失の分布データから当該分布データに対応付けられた単位面積あたりの遮蔽物の数を求め、前記送信機と前記受信機との間における遮蔽物の混雑度として出力する混雑度出力処理と
を実行することを特徴とする推定方法。 - 無線通信を行う送信機と受信機との間の通信領域における遮蔽物の混雑度を推定する推定装置において、
受信電力を取得する受信電力取得部と、
前記受信電力取得部により取得された受信電力から前記送信機と前記受信機との間の伝搬損失を計算する伝搬損失計算部と、
予め決められた所定時間に前記伝搬損失計算部により計算された複数の伝搬損失の中央値を求め、当該中央値を以降の処理に用いる伝搬損失として出力する中央値処理部と、
前記中央値処理部により出力された伝搬損失の最大値を取得する損失ピーク取得部と、
前記損失ピーク取得部により取得された伝搬損失の前記最大値と、前記送信機および前記受信機のアンテナの高さと、前記無線通信の波長と、に少なくとも基づいて、前記送信機または前記受信機と遮蔽物との距離を計算する距離算出部と、
前記距離算出部により計算された前記距離に基づいて、遮蔽物の個数を時系列に算出し、前記個数を予め決められた所定時間だけ積分することにより、単位面積当たりの遮蔽物の個数を求め、混雑度を推定する混雑度計算部と
を有することを特徴とする推定装置。 - 無線通信を行う送信機と受信機との間の通信領域における遮蔽物の混雑度を推定する推定装置において、
受信電力を取得する受信電力取得部と、
前記受信電力取得部により取得された受信電力から前記送信機と前記受信機との間の伝搬損失を計算する伝搬損失計算部と、
予め決められた所定時間に前記伝搬損失計算部により計算された複数の伝搬損失の中央値を求め、当該中央値を以降の処理に用いる伝搬損失として出力する中央値処理部と、
前記中央値処理部の出力結果から伝搬損失の分布を求め、予め保持している単位面積あたりの遮蔽物の数に対応する複数の分布データと突合せを行い、一番近い伝搬損失の分布データを選定する分布突合部と、
分布突合部により選定された伝搬損失の分布データから当該分布データに対応付けられた単位面積あたりの遮蔽物の数を求め、前記送信機と前記受信機との間における遮蔽物の混雑度として出力する混雑度出力部と
を有することを特徴とする混雑度の推定装置。 - 請求項3または請求項4に記載の推定装置で行う処理をコンピュータに実行させることを特徴とする混雑度の推定プログラム。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US17/289,214 US11909459B2 (en) | 2018-10-29 | 2019-10-21 | Congestion degree estimation method using propagation loss of wireless communication, estimation device and estimation program |
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2018202849A JP7047705B2 (ja) | 2018-10-29 | 2018-10-29 | 無線通信の伝搬損失を利用した混雑度の推定方法、推定装置および推定プログラム |
JP2018-202849 | 2018-10-29 |
Publications (1)
Publication Number | Publication Date |
---|---|
WO2020090540A1 true WO2020090540A1 (ja) | 2020-05-07 |
Family
ID=70464019
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
PCT/JP2019/041290 WO2020090540A1 (ja) | 2018-10-29 | 2019-10-21 | 無線通信の伝搬損失を利用した混雑度の推定方法、推定装置および推定プログラム |
Country Status (3)
Country | Link |
---|---|
US (1) | US11909459B2 (ja) |
JP (1) | JP7047705B2 (ja) |
WO (1) | WO2020090540A1 (ja) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP7300654B2 (ja) | 2020-08-11 | 2023-06-30 | 日本電信電話株式会社 | 電波遮蔽損失算出方法、電波遮蔽損失算出装置、及び電波遮蔽損失算出プログラム |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2012111138A1 (ja) * | 2011-02-18 | 2012-08-23 | 株式会社日立製作所 | 歩行者移動情報検出装置 |
JP2016014575A (ja) * | 2014-07-01 | 2016-01-28 | 日本信号株式会社 | 車両検知装置 |
JP2017116301A (ja) * | 2015-12-21 | 2017-06-29 | 株式会社リコー | 測位システム |
US20180165965A1 (en) * | 2016-12-08 | 2018-06-14 | Robert Bosch Gmbh | Method and device for detecting at least one pedestrian by a vehicle |
Family Cites Families (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP4528946B2 (ja) | 2005-01-06 | 2010-08-25 | 広島市 | 屋内侵入検知システムおよび屋内侵入検知情報発信システム |
JP6631395B2 (ja) | 2016-04-28 | 2020-01-15 | 富士通株式会社 | 位置・出力決定装置、位置・出力決定方法、および位置・出力決定プログラム |
KR101901359B1 (ko) * | 2016-11-21 | 2018-09-21 | 광운대학교 산학협력단 | 이동 노드를 이용한 수동적 위치 추적 시스템 및 방법 |
CN109151968B (zh) * | 2017-06-16 | 2021-02-23 | 华为技术有限公司 | 一种功率确定方法、设备及系统 |
-
2018
- 2018-10-29 JP JP2018202849A patent/JP7047705B2/ja active Active
-
2019
- 2019-10-21 US US17/289,214 patent/US11909459B2/en active Active
- 2019-10-21 WO PCT/JP2019/041290 patent/WO2020090540A1/ja active Application Filing
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2012111138A1 (ja) * | 2011-02-18 | 2012-08-23 | 株式会社日立製作所 | 歩行者移動情報検出装置 |
JP2016014575A (ja) * | 2014-07-01 | 2016-01-28 | 日本信号株式会社 | 車両検知装置 |
JP2017116301A (ja) * | 2015-12-21 | 2017-06-29 | 株式会社リコー | 測位システム |
US20180165965A1 (en) * | 2016-12-08 | 2018-06-14 | Robert Bosch Gmbh | Method and device for detecting at least one pedestrian by a vehicle |
Also Published As
Publication number | Publication date |
---|---|
US20220014285A1 (en) | 2022-01-13 |
JP7047705B2 (ja) | 2022-04-05 |
US11909459B2 (en) | 2024-02-20 |
JP2020071042A (ja) | 2020-05-07 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP2022028703A (ja) | 物体追跡のための方法、装置、サーバ及びシステム | |
Gomez-de-Gabriel et al. | Monitoring harness use in construction with BLE beacons | |
US7203489B2 (en) | Measurement based prediction method for radiation path loss | |
EP2171890B1 (en) | Method for noise floor and interference estimation | |
JP6251930B2 (ja) | 位置推定システム | |
CN103323838B (zh) | 雷达装置 | |
KR101063318B1 (ko) | 이동 단말의 위치 측정 시스템 및 그 방법 | |
US20130214962A1 (en) | Apparatus and method for traffic lane detection | |
WO2020090540A1 (ja) | 無線通信の伝搬損失を利用した混雑度の推定方法、推定装置および推定プログラム | |
DE102019128729A1 (de) | Fahrunterstützungsvorrichtung, Fahrzeug und Fahrunterstützungsverfahren | |
WO2020126050A1 (en) | Coordination of wireless communication unit and radar unit in a wireless communication network | |
JP6755876B2 (ja) | 測距方法及び装置 | |
KR101694521B1 (ko) | 전파지문지도 생성 장치 및 방법 | |
US8063824B2 (en) | Magnitude-squared coherence filtering in a time difference of arrival (TDOA) or group delay interferometric geolocation system | |
KR102155419B1 (ko) | 인공신경망을 기반으로 하는 적응형 거리측정장치 | |
US20230155422A1 (en) | Information processing device, information processing system, data generation method, and non-transitory computer readable medium | |
KR101423711B1 (ko) | 전파의 전파 경로 분석 방법 및 장치 | |
KR101245522B1 (ko) | 무선 측위 방법 및 시스템 | |
JP7279801B2 (ja) | 電波発信源位置推定システム及び電波発信源位置推定方法 | |
US20190242987A1 (en) | Radar device and method of detecting object using radar device | |
Lai et al. | The characterisation of human body influence on indoor 3.5 GHz path loss measurement | |
DE102021103349A1 (de) | Kommunikationsvorrichtung, Informationsverarbeitungsverfahren und Speichermedium | |
Makhoul et al. | Wideband measurement-based vehicle-to-pedestrian channel models | |
CN110888113A (zh) | 车载雷达控制方法、装置、设备和存储介质 | |
JP7443759B2 (ja) | 位置推定装置、位置推定方法、プログラム及び位置推定システム |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
121 | Ep: the epo has been informed by wipo that ep was designated in this application |
Ref document number: 19877937 Country of ref document: EP Kind code of ref document: A1 |
|
NENP | Non-entry into the national phase |
Ref country code: DE |
|
122 | Ep: pct application non-entry in european phase |
Ref document number: 19877937 Country of ref document: EP Kind code of ref document: A1 |