JP4032125B1 - Mobile object position estimation system and mobile object position estimation method - Google Patents

Mobile object position estimation system and mobile object position estimation method Download PDF

Info

Publication number
JP4032125B1
JP4032125B1 JP2006185198A JP2006185198A JP4032125B1 JP 4032125 B1 JP4032125 B1 JP 4032125B1 JP 2006185198 A JP2006185198 A JP 2006185198A JP 2006185198 A JP2006185198 A JP 2006185198A JP 4032125 B1 JP4032125 B1 JP 4032125B1
Authority
JP
Japan
Prior art keywords
probability
data
particles
position estimation
moving object
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
JP2006185198A
Other languages
Japanese (ja)
Other versions
JP2008014742A (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.)
Japan Advanced Institute of Science and Technology
Original Assignee
Japan Advanced Institute of Science and Technology
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 Japan Advanced Institute of Science and Technology filed Critical Japan Advanced Institute of Science and Technology
Priority to JP2006185198A priority Critical patent/JP4032125B1/en
Application granted granted Critical
Publication of JP4032125B1 publication Critical patent/JP4032125B1/en
Publication of JP2008014742A publication Critical patent/JP2008014742A/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Traffic Control Systems (AREA)
  • Navigation (AREA)
  • Alarm Systems (AREA)

Abstract

【課題】 移動体の位置を高い精度で推定可能な移動体位置推定システム及び移動体位置推定方法を提供する。
【解決手段】 移動体位置計測装置から移動体の位置情報を取得すると、当該位置情報を用いて地図データの座標系に分散するパーティクルをサンプリングするパーティクルフィルタを備え、当該パーティクルフィルタによりサンプリングしたパーティクルから当該移動体の位置を推定する移動体位置推定システム及び移動体位置推定方法であり、当該地図データの座標系上の各領域に当該移動体が存在する確率を示す確率データを有し、当該確率データに従って当該地図データの座標系にパーティクルを分散させることにより、移動体の位置を高い精度で推定する。
【選択図】 図3
PROBLEM TO BE SOLVED: To provide a moving body position estimation system and a moving body position estimation method capable of estimating the position of a moving body with high accuracy.
When position information of a moving object is acquired from a moving object position measuring device, a particle filter that samples particles dispersed in a coordinate system of map data using the position information is provided, and the particles sampled by the particle filter are used. A mobile object position estimation system and a mobile object position estimation method for estimating the position of the mobile object, the probability data indicating the probability that the mobile object exists in each area on the coordinate system of the map data, and the probability The position of the moving body is estimated with high accuracy by dispersing particles in the coordinate system of the map data according to the data.
[Selection] Figure 3

Description

本発明は、移動体位置計測装置から得られる位置情報を使用して、移動体の位置を推定する移動体位置推定システム、及び、移動体位置推定方法に関する。   The present invention relates to a mobile body position estimation system and a mobile body position estimation method for estimating the position of a mobile body using position information obtained from a mobile body position measurement apparatus.

近年、GPSなどの移動体位置計測装置によって、移動体の位置情報を取得することが容易に可能となってきた。しかし、その位置情報は精度が低く、移動体位置計測装置単体では数十メートルから数百メートルの誤差が発生するのが現状である。   In recent years, it has become possible to easily acquire position information of a moving body by a moving body position measuring device such as a GPS. However, the position information has low accuracy, and an error of several tens of meters to several hundreds of meters is generated in the mobile body position measuring device alone.

その不確実性を軽減するために、移動体位置計測装置から得た位置情報に基づいて、コンピュータにより移動体の位置を推定する移動体位置推定システムが用いられている。下記非特許文献1には、移動体位置推定システムにパーティクルフィルタを用いて位置を確率的に扱う技術が提案されている。パーティクルフィルタとは、複数個のパーティクルの重み付き総和によって全体の確率分布を近似する手法であり、位置の推定は次のように行われる。(1)地図データの座標系上のランダムな位置にN個のパーティクルを生成する。(2)パーティクルをランダムウォークさせる。(3)移動体位置計測装置からの位置情報を確率分布として示す。(4)確率分布に従って個々のパーティクル位置の生起確率を求める。(5)生起確率に比例させて、各パーティクルの重みを更新する。(6)各パーティクルの重みと位置とから移動体の位置について正規分布を最尤推定する。(7)正規分布の平均を移動体の位置と推定する。
J.Hightower,B.Borriello.“The Location Stack:A Layerd Model for Location in Ubiquitous Computing”,4th IEEE Workshop on Modile Computing Systems & Applications (WMCSA 2002),pp.22−28.2002.
In order to reduce the uncertainty, a moving body position estimation system is used in which the position of the moving body is estimated by a computer based on position information obtained from the moving body position measuring device. The following Non-Patent Document 1 proposes a technique for probabilistically handling a position using a particle filter in a moving body position estimation system. The particle filter is a method of approximating the entire probability distribution by a weighted sum of a plurality of particles, and the position is estimated as follows. (1) N particles are generated at random positions on the coordinate system of the map data. (2) Randomly walk particles. (3) Position information from the moving body position measuring device is shown as a probability distribution. (4) The occurrence probability of each particle position is obtained according to the probability distribution. (5) The weight of each particle is updated in proportion to the occurrence probability. (6) Maximum likelihood estimation of the normal distribution for the position of the moving object from the weight and position of each particle. (7) The average of the normal distribution is estimated as the position of the moving object.
J. et al. Highwater, B.W. Borriello. “The Location Stack: A Layer Model for Location in Ubiquitous Computing”, 4th IEEE Workshop on Modular Computing Systems & Applications 2 (WMC). 22-28.2002.

しかしながら、上記従来技術は、移動体の位置を移動体位置計測装置からの測定値のみで推定しているため、推定精度は移動体位置計測装置に依存してしまう。現在の移動体位置計測装置による位置測定は、数十メートルから数百メートルの誤差を含んでいるため、位置推定にもその誤差が影響してしまい、精度の良い推定を行うことが困難であった。また、パーティクルから一つの正規分布を求め、その平均を移動体の位置と推定するため、出力される推定位置は一つである。移動体位置計測装置の誤差が大きいときは、一つの推定値に集約すると推定確率の分散が大きくなり、推定精度が低くなってしまう。   However, since the above prior art estimates the position of the moving body only by the measurement value from the moving body position measuring device, the estimation accuracy depends on the moving body position measuring device. The current position measurement by the mobile object position measurement apparatus includes an error of several tens of meters to several hundreds of meters. Therefore, the error also affects the position estimation, and it is difficult to perform accurate estimation. It was. In addition, since one normal distribution is obtained from the particles and the average is estimated as the position of the moving object, one estimated position is output. When the error of the moving body position measuring device is large, the aggregation of the estimation probabilities becomes large and the estimation accuracy becomes low when the estimation is made into one estimation value.

そこで、本発明は、移動体の位置を高い精度で推定可能な移動体位置推定システム、及び、移動体位置推定方法を提供することを目的とする。   Then, an object of this invention is to provide the mobile body position estimation system which can estimate the position of a mobile body with high precision, and a mobile body position estimation method.

本発明の移動体位置推定システムは、移動体位置計測装置から移動体の位置情報を取得すると、当該位置情報を用いて地図データの座標系に分散するパーティクルをサンプリングするパーティクルフィルタを備え、当該パーティクルフィルタによりサンプリングしたパーティクルから当該移動体の位置を推定する移動体位置推定システムにおいて、当該地図データの座標系上の各領域に当該移動体が存在する確率を示す確率データを有し、当該確率データに従って当該地図データの座標系にパーティクルを分散させる手段を備えることを特徴とする。また、本発明の移動体位置推定方法は、コンピュータが、移動体位置計測装置から移動体の位置情報を取得すると、パーティクルフィルタにより当該位置情報を用いて地図データの座標系に分散するパーティクルをサンプリングし、当該サンプリングしたパーティクルから当該移動体の位置を推定する移動体位置推定方法において、当該地図データの座標系上の各領域に当該移動体が存在する確率を示す確率データを参照し、当該確率データに従って当該地図データの座標系にパーティクルを分散させるステップを備えることを特徴とする。   The mobile object position estimation system of the present invention includes a particle filter that samples particles dispersed in a coordinate system of map data using the position information when the position information of the mobile object is acquired from the mobile object position measurement device. In the moving body position estimation system that estimates the position of the moving body from the particles sampled by the filter, the moving body position estimation system has probability data indicating the probability that the moving body exists in each region on the coordinate system of the map data, and the probability data And a means for dispersing particles in the coordinate system of the map data. Further, according to the moving object position estimation method of the present invention, when the computer acquires the position information of the moving object from the moving object position measuring device, the particle filter uses the position information to sample particles dispersed in the coordinate system of the map data. In the moving object position estimation method for estimating the position of the moving object from the sampled particles, the probability data indicating the probability that the moving object exists in each region on the coordinate system of the map data is referred to. Dispersing particles in the coordinate system of the map data according to the data.

ここで、領域とは、建物や道路や河川といった物理的に区別される領域でも良いし、住所といった人為的に定められた領域でも良いし、ユーザが任意に設定した領域でも良い。地図上の各領域に移動体が存在する確率データは、ユーザにより任意に設定可能としてもよいし、予めシステム内に設定されていてもよい。   Here, the area may be a physically distinguished area such as a building, road, or river, an artificially defined area such as an address, or an area arbitrarily set by the user. The probability data that the moving body exists in each area on the map may be arbitrarily set by the user, or may be set in the system in advance.

移動体が存在する確率は、例えば人が携帯する情報端末であれば、道路や建物が高く、河川や山中は低いなど、領域によって異なる。また、移動体の移動パターンや移動範囲によっても異なる場合がある。本発明によれば、地図上の各領域に移動体が存在する確率を示す確率データを参照し、その確率に応じてパーティクルを移動させるため、その確率を反映させた位置の推定が可能となる。   The probability that a moving object exists differs depending on the area, for example, in the case of an information terminal carried by a person, roads and buildings are high and rivers and mountains are low. Moreover, it may differ depending on the movement pattern and movement range of the moving body. According to the present invention, the probability data indicating the probability that a moving body exists in each area on the map is referred to, and the particles are moved according to the probability. Therefore, the position reflecting the probability can be estimated. .

また、前記確率データは、前記各領域がカテゴリ分けされており、当該カテゴリごとに前記確率を設定可能であることが好ましい。   In the probability data, it is preferable that the respective areas are classified into categories, and the probability can be set for each category.

各領域における移動体が存在する確率は、例えば、移動体が人であれば建物や道路が高く河川は低く、また、移動体が船舶であれば河川が高く陸地は低いというように、領域ごとに区別できる。この発明によれば、前記確率データは、各領域をカテゴリ分けし、そのカテゴリごとに確率を設定することが可能となっているため、移動体に応じて各領域の確率を容易に設定することができる。   The probability that a moving body exists in each area is, for example, that if the moving body is a person, buildings and roads are high and the river is low, and if the moving body is a ship, the river is high and the land is low. Can be distinguished. According to the present invention, since the probability data can categorize each region and set the probability for each category, the probability of each region can be easily set according to the moving object. Can do.

また、前記地図データには前記各領域を識別する識別情報が含まれており、前記確率データは、当該識別情報を用いて前記各領域がカテゴリ分けされていることが好ましい。   The map data preferably includes identification information for identifying each area, and the probability data is preferably classified into the areas using the identification information.

地図は道や建物や川などの領域ごとに色分けされたり、マークが付されたりして、各領域を識別する識別情報が付されていることが多い。この発明によれば、確率データは、これらの識別情報を利用してカテゴリが分けられるため、既存の地図データを用いて簡単にカテゴリ分けを行うことができる。   In many cases, the map is color-coded for each area such as a road, a building, or a river, or is given a mark, and identification information for identifying each area is attached. According to the present invention, the probability data can be easily categorized using existing map data because the categories are divided using the identification information.

前記移動体位置推定システム/方法は、 前記サンプリングされたパーティクルをクラスタリングし、各クラスターの平均値を前記移動体の推定位置として出力する手段/ステップを備えることが好ましい。   The moving body position estimation system / method preferably includes means / steps for clustering the sampled particles and outputting an average value of each cluster as an estimated position of the moving body.

パーティクル全体から移動体の位置を推定した場合、移動体位置計測装置からの位置情報の誤差が大きいときには、推定確率の分散が大きくなり、推定精度が低くなる。この発明によれば、パーティクルをクラスタリングし、各クラスターの平均値を推定位置とするため、各クラスター内における推定確率の分散を抑えることができ、推定精度を高めることができる。   When the position of the moving object is estimated from the whole particle, when the error in the position information from the moving object position measuring device is large, the variance of the estimation probability increases and the estimation accuracy decreases. According to the present invention, particles are clustered, and the average value of each cluster is used as the estimated position. Therefore, dispersion of the estimation probability within each cluster can be suppressed, and the estimation accuracy can be increased.

本発明に係る移動体位置推定ステム/方法によれば、道路や建物や河川といった領域ごとに移動体が存在する確率を示す確率データを有し、移動体の位置推定に際しては、その確率データを反映させるため、移動体の位置をより高い精度で推定可能となる。   According to the moving object position estimation stem / method according to the present invention, the moving object position estimation stem / method has probability data indicating the probability that a moving object exists for each region such as a road, a building, and a river. In order to reflect this, the position of the moving body can be estimated with higher accuracy.

確率データの領域を例えば道路や河川などのようにカテゴリに分け、カテゴリごとに確率を設定可能とすれば、推定対象となる移動体に適する確率を容易に設定することができる。確率データを生成する処理は、カテゴリ分けを地図の各領域に付される識別情報を利用して行うことにより、既存の地図情報を用いて簡単にカテゴリ分けを行うことができる。   If the probability data area is divided into categories such as roads and rivers, and the probability can be set for each category, the probability suitable for the mobile object to be estimated can be easily set. The process of generating the probability data can be easily categorized using the existing map information by performing the categorization using the identification information attached to each area of the map.

また、パーティクルをクラスタリングし、各クラスターの平均値を移動体の位置として出力することにより、各クラスター内における推定確率の分散を抑えることができ、推定精度を高めることができる。また、推定位置として複数の候補を得ることができる。たとえば、移動体を捜索する場合、従来技術のように推定位置が唯一に出力されると捜索場所も一つに限られてしまう。本システム/方法によれば、複数の場所を候補として捜索できるため発見できる可能性も高まる。   Also, by clustering particles and outputting the average value of each cluster as the position of the moving body, dispersion of estimation probabilities within each cluster can be suppressed, and estimation accuracy can be increased. In addition, a plurality of candidates can be obtained as estimated positions. For example, when searching for a mobile object, if the estimated position is output only as in the prior art, the search location is limited to one. According to the present system / method, a plurality of places can be searched as candidates, so that the possibility of discovery can be increased.

本発明は、様々な分野に応用可能である。たとえば、GPS機能付き情報携帯端末等の移動体を人に携帯させ、徘徊する高齢者や遭難者を捜索したり、登下校中の児童の位置を把握して安全を確保したりするのに用いることができる。また、移動体を自動車等に搭載すれば、タクシーや配送車等の位置管理などにも応用可能である。   The present invention can be applied to various fields. For example, a mobile object such as an information portable terminal with GPS function can be carried by a person and used to search for elderly people and victims who are hesitant or to ascertain the location of children who are going to and from school. be able to. In addition, if the moving body is mounted on an automobile or the like, it can be applied to the position management of a taxi or a delivery car.

本実施の形態の移動体位置推定システム1は、移動体位置計測装置から得られる位置情報を取得すると、当該位置情報を用いて地図データの座標系に分散するパーティクルをサンプリングするパーティクルフィルタを備え、当該パーティクルフィルタによりサンプリングしたパーティクルから当該移動体の位置を推定するシステムである。   The mobile body position estimation system 1 according to the present embodiment includes a particle filter that samples the particles dispersed in the coordinate system of the map data using the position information when the positional information obtained from the mobile body position measurement device is acquired. This is a system for estimating the position of the moving object from the particles sampled by the particle filter.

図1は、本発明の第1の実施の形態による移動体位置推定システム1の一構成を示すブロック図である。本実施の形態による移動体位置推定システム1は、内部バス11に、通信インタフェース12、CPU13、ROM14、RAM15、ディスプレイ16、キーボード/マウス17、ドライブ18、ハードディスク19を接続させ、アドレス信号、制御信号、データ等を伝送させ、本実施の形態による移動体位置推定システムを実現するコンピュータシステムとしての構成を備えている。   FIG. 1 is a block diagram showing a configuration of a moving object position estimation system 1 according to the first embodiment of the present invention. The mobile body position estimation system 1 according to the present embodiment connects an internal bus 11 with a communication interface 12, a CPU 13, a ROM 14, a RAM 15, a display 16, a keyboard / mouse 17, a drive 18, and a hard disk 19, and an address signal and a control signal. The computer system is configured to transmit data and the like and realize the moving object position estimation system according to the present embodiment.

通信インタフェース12は、移動体位置計測装置との通信やインターネット等の通信網に接続する各種通信機能を有しており、移動体位置計測装置から移動体の位置情報を受信したり、本発明に係る位置推定をコンピュータに実行させるプログラム19aや地図データ19cをダウンロードしたりすることも可能である。CPU13は、ROM14に格納されたOSにより装置全体の制御を行うとともにハードディスク19に格納された各種のアプリケーションプログラム19aに基づいて処理を実行する機能を司る。   The communication interface 12 has various communication functions for communication with a mobile body position measurement device and a communication network such as the Internet, and receives position information of the mobile body from the mobile body position measurement device. It is also possible to download a program 19a or map data 19c that causes a computer to execute such position estimation. The CPU 13 controls the entire apparatus by the OS stored in the ROM 14 and controls the function of executing processing based on various application programs 19 a stored in the hard disk 19.

ROM14は、OS等のように装置全体の制御を行うためのプログラムを格納しており、これらをCPU13に供給する機能を有している。RAM15は、CPU13による各種プログラムの実行時にワークエリアとして利用されるメモリ機能を有している。   The ROM 14 stores a program for controlling the entire apparatus, such as an OS, and has a function of supplying these to the CPU 13. The RAM 15 has a memory function used as a work area when the CPU 13 executes various programs.

ディスプレイ16は、地図データや確率データや推定位置をグラフィカルに表示したりする機能を有している。キーボード/マウス17は、文字、数字、記号等のデータを入力したり、画面上のポイント位置を操作することが可能である。   The display 16 has a function of graphically displaying map data, probability data, and an estimated position. The keyboard / mouse 17 can input data such as characters, numbers, symbols, etc., and can operate the point position on the screen.

ドライブ18は、各種のプログラム、データを記録したCD、DVD等の記録媒体からインストール作業を実行するための駆動ユニットである。コンピュータを移動体位置推定システム1として機能させるプログラムや地図データを記憶媒体からインストールしたりすることも可能である。   The drive 18 is a drive unit for executing installation work from a recording medium such as a CD or DVD in which various programs and data are recorded. It is also possible to install a program or map data that causes a computer to function as the moving body position estimation system 1 from a storage medium.

ハードディスク19は、プログラム19a、メモリ19b、位置推定に用いられる地図データ19c等の各種データを記憶する外部記憶装置である。プログラム19aは、前述した通信インタフェース12、ドライブ18等からインストールされたプログラムを実行形式で記憶したものに相当する。メモリ19bは、各種プログラムの実行結果等のファイルを保存する記憶部である。   The hard disk 19 is an external storage device that stores various data such as a program 19a, a memory 19b, and map data 19c used for position estimation. The program 19a corresponds to a program stored from the communication interface 12, the drive 18 and the like described above in the execution format. The memory 19b is a storage unit that stores files such as execution results of various programs.

地図データ19cは、例えば図2に示すように、基本データD1と確率データD2とを含む。基本データD1は、建物・道路・河川等や、都道府県、市町村、住所などの領域や、位置、経緯度、方角、距離等のデータが含まれている。また、確率データD2は、基本データD1の各領域に移動体が存在する確率を示す情報が含まれている。確率データD2は、建物がカテゴリa、道路がカテゴリb、河川がカテゴリc、その他がカテゴリdというように、領域がカテゴリ分けされており、カテゴリごとに確率が設定されている。この確率はその領域に推定対象となる移動体が存在する確率である。推定対象として人が携帯する情報端末であれば、建物に対応するカテゴリaや道路に対応するカテゴリbの確率は高く、河川に対応するカテゴリcの確率は低く設定されている。ただし、確率データD2は、カテゴリ分けすることなく、各建物や、各道路ごとというように、領域ごとに確率が設定されていても良いし、住所などの人為的に分けられた領域に設定されていても良い。これにより、移動体固有の行動パターンや行動範囲を確率に反映させることができる。また、確率データD2は地図データ19cの一部として一体となっていても良いし、別データとして記憶されていても良い。   The map data 19c includes basic data D1 and probability data D2, for example, as shown in FIG. The basic data D1 includes data such as areas such as buildings, roads, rivers, prefectures, municipalities, and addresses, positions, longitude and latitude, directions, and distances. The probability data D2 includes information indicating the probability that a moving object exists in each area of the basic data D1. In the probability data D2, the areas are classified into categories such that the building is category a, the road is category b, the river is category c, and the others are category d, and the probability is set for each category. This probability is a probability that a mobile object to be estimated exists in the region. In the case of an information terminal carried by a person as an estimation target, the probability of category a corresponding to a building or category b corresponding to a road is high, and the probability of category c corresponding to a river is set low. However, the probability data D2 may be set for each area, such as each building or each road, without being classified into categories, or set in an artificially divided area such as an address. May be. Thereby, the action pattern and action range peculiar to a moving body can be reflected in the probability. The probability data D2 may be integrated as part of the map data 19c, or may be stored as separate data.

確率データD2は、後述するステップS1の処理前の事前処理で生成され、ハードディスク19に記憶されている。また、確率データD2は予めシステム1内に設定されていても良いし、ユーザが任意に設定や変更を可能としても良い。たとえば、地図データ19cの基本データD1をグラフィカルに表示し、各建物や道路等ごとに確率を直接設定可能とし、ユーザが事前設定として確率やカテゴリを設定するようにして確率データD2を生成可能としても良い。また、基本データD1が建物や道路などの領域を識別する識別情報を含む場合(例えば、建物や道路などが色分けされていたり、マークが付されていたりする場合)、移動体位置推定システム1が、その識別情報を利用して領域をカテゴリ分けし、そのカテゴリを確率データD2のカテゴリとしても良い。これにより、人為的なカテゴリ分けが不要となり、既存の地図データを用いて、領域のカテゴリ分けを簡単に行うことができる。   The probability data D <b> 2 is generated by pre-processing before processing in step S <b> 1 described later, and is stored in the hard disk 19. The probability data D2 may be set in the system 1 in advance, or may be arbitrarily set or changed by the user. For example, the basic data D1 of the map data 19c is displayed graphically, the probability can be set directly for each building, road, etc., and the probability data D2 can be generated by setting the probability and category as a preset by the user. Also good. In addition, when the basic data D1 includes identification information for identifying an area such as a building or a road (for example, when a building or a road is color-coded or marked), the mobile object position estimation system 1 is The areas may be classified into categories using the identification information, and the categories may be used as the categories of the probability data D2. This eliminates the need for artificial categorization and makes it possible to easily categorize areas using existing map data.

本発明におけるパーティクルとは、地図データ19cの座標系上の位置を示す指標であり、属性には少なくとも地図データ19cの座標系において位置を特定可能な位置情報を含む。本発明におけるパーティクルフィルタとは、移動体位置計測装置から得られる位置情報を取得すると、地図データ19cの座標系に分散する複数のパーティクルから、上記位置情報に応じてパーティクルをサンプリングする機能を少なくとも備えるものであり、本実施の形態においては少なくとも下記のステップS1からS5を実行する機能を備えるフィルタである。   The particle in the present invention is an index indicating the position of the map data 19c on the coordinate system, and the attribute includes at least position information capable of specifying the position in the coordinate system of the map data 19c. The particle filter in the present invention has at least a function of sampling particles according to the position information from a plurality of particles dispersed in the coordinate system of the map data 19c when the position information obtained from the moving body position measuring device is acquired. In the present embodiment, the filter has a function of executing at least the following steps S1 to S5.

図3は、移動体位置推定システム1を用いた移動体位置推定方法を説明するフローチャートである。図4は、位置推定処理を概念的に説明する説明図である。移動体位置推定システム1は下記のように機能して、本発明の移動体位置推定方法を実現する。   FIG. 3 is a flowchart for explaining a mobile object position estimation method using the mobile object position estimation system 1. FIG. 4 is an explanatory diagram conceptually illustrating the position estimation process. The mobile body position estimation system 1 functions as follows to realize the mobile body position estimation method of the present invention.

(ステップS1)地図データ19cの座標系のランダムな位置に複数のパーティクルを配置する。パーティクルの配置は、各パーティクルの位置情報をランダムに設定することにより行われる。図4(1)に示すように、パーティクル(図中、白円)はランダムな位置に配置される。   (Step S1) A plurality of particles are arranged at random positions in the coordinate system of the map data 19c. The arrangement of the particles is performed by setting the position information of each particle at random. As shown in FIG. 4A, particles (white circles in the figure) are arranged at random positions.

(ステップS2)地図データ19cの確率データD2を参照し、その確率に依存させてすべてのパーティクルを移動(ランダムウォーク)させる。たとえば、図5に示すように、確率データD2は、基本データD1の緯度経度座標系をピクセル単位に変換し、基本データD1でカテゴリa,b,c,dに分けられた各領域にピクセル単位で確率を設定したものとする。そして、パーティクルの移動可能範囲eをパーティクルの位置を中心とした一定半径サークル内と設定しておき、パーティクルをランダムウォークさせるときには、移動可能範囲内eのすべてのピクセルの確率データを取得し、その確率に比例して次の位置を決定し、その位置にパーティクルを移動させる。パーティクルを移動させるタイミングは予め移動体位置推定システム1に設定されており、そのタイミングになると確率データD2に従ってパーティクルをランダムウォークさせる。図4(2)は、確率データに従ってパーティクルをランダムウォークさせたときの状態を概念的に示す説明図である。確率の高い領域に多くのパーティクルが配置される。   (Step S2) With reference to the probability data D2 of the map data 19c, all particles are moved (random walk) depending on the probability. For example, as shown in FIG. 5, the probability data D2 is obtained by converting the latitude / longitude coordinate system of the basic data D1 into pixel units, and in each area divided into categories a, b, c, and d in the basic data D1. The probability is set with. Then, the particle movable range e is set within a constant radius circle centered on the particle position, and when the particles are randomly walked, the probability data of all the pixels in the movable range e is acquired, The next position is determined in proportion to the probability, and the particle is moved to that position. The timing for moving the particles is set in the moving body position estimation system 1 in advance, and the particles are randomly walked according to the probability data D2 at that timing. FIG. 4B is an explanatory diagram conceptually showing a state when particles are randomly walked according to probability data. Many particles are placed in a region with a high probability.

(ステップS3)移動体位置計測装置から移動体の位置情報を取得しない間は、ステップS2に戻り上記ランダムウォークを続ける。移動体位置計測装置から移動体の位置情報を取得した場合は、ステップS4に進む。   (Step S3) While the position information of the moving object is not acquired from the moving object position measuring device, the process returns to Step S2 and continues the random walk. If the position information of the moving object is acquired from the moving object position measuring device, the process proceeds to step S4.

(ステップS4)移動体位置推定システム1は、位置情報を取得すると、その位置情報を確率分布P1(x)で表されるデータに変換する。位置情報から確率分布P1(x)への変換は、予め設定されている数式等の規則に従い変換する。移動体位置計測装置からの位置情報の形式は、その機器の特性に応じて異なるので、機器に応じた確率分布に変換する。このため、移動体位置推定システム1には移動体位置計測装置ごとに予め数式等による変換規則が設定されており、位置情報から確率分布P1(x)への変換はその規則に従い行う。例えば、GPSの場合は緯度経度情報を出力するので、取得した緯度経度情報を平均とした正規分布に変換する。また、PHSは各地に所在する基地局のどの基地局で端末からのデータを受信したかで位置が特定される。このため、移動体からの電波を受信した基地局が電波を受信できる範囲内は、移動体が存在する確率が一様であり、その他の範囲は確率が0である。そこで、PHSから位置情報を取得した場合は、受信した基地局が電波を受信できる範囲内において移動体が存在する確率を0より大の一定値(0より大)とし、その他の範囲は確率を0とする一様分布に変換する。   (Step S4) When the mobile body position estimation system 1 acquires the position information, the mobile body position estimation system 1 converts the position information into data represented by the probability distribution P1 (x). The conversion from the position information to the probability distribution P1 (x) is performed according to a rule such as a mathematical formula set in advance. Since the format of the position information from the moving body position measuring device differs depending on the characteristics of the device, it is converted into a probability distribution according to the device. For this reason, the conversion rule by numerical formula etc. is preset for every moving body position measuring apparatus in the moving body position estimation system 1, and the conversion from position information to probability distribution P1 (x) is performed according to the rule. For example, in the case of GPS, latitude / longitude information is output, so that the latitude / longitude information obtained is converted into a normal distribution with an average. In addition, the position of the PHS is specified depending on which base station of the base stations located in each place has received data from the terminal. For this reason, the probability that the mobile body exists is uniform within the range in which the base station that has received the radio wave from the mobile body can receive the radio wave, and the probability is 0 in the other ranges. Therefore, when the location information is acquired from the PHS, the probability that the mobile body exists within the range in which the received base station can receive radio waves is set to a constant value (greater than 0) greater than 0, and the other ranges have the probability. Convert to a uniform distribution of zero.

(ステップS5)確率分布P(1)に従って、地図上の母集団からパーティクルをサンプリングする。サンプリングの方法は、まず、その確率分布P1(x)に従った確率で1つのパーティクルを地図上の母集団から抽出し、それを母集団に戻してから再び同様に次のパーティクルを抽出する。それを繰り返すことにより、母集団と同じ数のパーティクルを確率分布P1(x)に従って抽出する。図6は、GPSの位置情報を用いたときの確率分布P1(x)のデータと、地図データ19cの座標系上の母集団と、抽出されたパーティクルとを説明する説明図である。取得した位置情報を平均とした正規分布データP1(x)に従って、母集団からパーティクルが抽出される。パーティクルは、正規分布の頂点(平均)に最も近いところが多く抽出され、遠いところほど抽出されない。図4(3)は、この正規分布によりサンプリングされたパーティクルの分散を概念的に説明する説明図である。パーティクルは、取得した位置情報に最も近いところが多く抽出され、遠いところほど抽出されない。   (Step S5) According to the probability distribution P (1), particles are sampled from the population on the map. In the sampling method, first, one particle is extracted from the population on the map with a probability according to the probability distribution P1 (x), and after returning it to the population, the next particle is similarly extracted again. By repeating this, the same number of particles as the population are extracted according to the probability distribution P1 (x). FIG. 6 is an explanatory diagram for explaining the probability distribution P1 (x) data when using GPS position information, the population on the coordinate system of the map data 19c, and the extracted particles. Particles are extracted from the population according to the normal distribution data P1 (x) with the acquired position information as an average. As for the particle, a part closest to the vertex (average) of the normal distribution is extracted, and a part farther away is not extracted. FIG. 4 (3) is an explanatory diagram for conceptually explaining the dispersion of particles sampled by the normal distribution. As for the particle, a part closest to the acquired position information is extracted, and a part farther away is not extracted.

(ステップS6)サンプリングされたパーティクルをクラスタリングする。図4(4)はパーティクルのクラスタリングを概念的に説明する説明図である。サンプリングされたパーティクルは複数のクラスターにクラスタリングされる。ここで、本実施の形態のクラスタリングについて説明する。本実施の形態のクラスタリングには、正規混合分布P(x)を用いたクラスタリング手法を用いる。正規混合分布とは数1に示される確率分布である。数1のp(x)は正規分布を表す。この数1における正規分布の数cと各正規分布の重みαi、および各正規分布p(x)の平均値と分散を求めることがクラスタリング処理にあたる。本実施の形態では、これらの値を総じて、正規混合分布P(x)のパラメータと呼ぶこととする。クラスタリングによって生成されるクラスターは、正規混合分布P(x)を構成する個々の正規分布に対応する。正規混合分布P(x)のパラメータが求まると、数1によって各パーティクルがどのクラスターに所属するかを確率的に得ることが可能となる。さらに、地図上の任意の位置を入力として与えると、その位置に移動体が居る確率を求めることができるようになる。

Figure 0004032125
(Step S6) The sampled particles are clustered. FIG. 4 (4) is an explanatory diagram for conceptually explaining the clustering of particles. Sampled particles are clustered into multiple clusters. Here, clustering according to the present embodiment will be described. For the clustering of the present embodiment, a clustering method using a normal mixture distribution P (x) is used. The normal mixture distribution is a probability distribution represented by Equation 1. P i (x) in Equation 1 represents a normal distribution. The clustering process is to obtain the number c of normal distributions in Equation 1, the weight α i of each normal distribution, and the average value and variance of each normal distribution p i (x). In the present embodiment, these values are collectively called parameters of the normal mixture distribution P (x). Clusters generated by clustering correspond to the individual normal distributions constituting the normal mixture distribution P (x). When the parameters of the normal mixture distribution P (x) are obtained, it is possible to obtain probabilistically which cluster each particle belongs to according to Equation 1. Furthermore, when an arbitrary position on the map is given as an input, the probability that a moving body is present at that position can be obtained.
Figure 0004032125

図7に、正規混合分布P(x)によるクラスタリング例を示す。水平線の上に並んでいる長方形が個々のパーティクルを表している。P(x)は正規混合分布を表す。個々のパーティクルから求められた正規混合分布P(x)は、パーティクルが集まっているところほど、確率が高く(山が高く)なっている。この例では2つのクラスターが生成されている。   FIG. 7 shows an example of clustering using the normal mixture distribution P (x). The rectangles on the horizon represent individual particles. P (x) represents a normal mixture distribution. The normal mixture distribution P (x) obtained from each particle has a higher probability (a mountain is higher) as the particles are gathered. In this example, two clusters are generated.

すなわちステップS6において移動体位置推定システム1は、以下のように動作する。数1は移動体位置推定システム1に予め設定されており、移動体位置推定システム1は、ステップS5にてサンプリングしたパーティクルの情報を用いて、数1の各パラメータ数cと各正規分布の重みαi、及び、正規分布p(x)の平均値と分散を算出する。詳しくは、各正規分布の重みαi、および正規分布p(x)の平均値と分散については、個々のパーティクルが所属するクラスターを隠れ変数としたEMアルゴリズムに基づいて算出する。EMアルゴリズムとは、隠れ変数がある場合の最尤推定方法として広く知られている手法であり、繰り返し計算を行うことで、解が収束することが証明されている。一方、正規混合分布P(x)の数cについては、初期値を1として数2に示される正規混合分布P(x)における各パーティクルの対数尤度の値が減少するまで1つずつ加算して求める。図4(5)は各クラスターの平均値を概念的に説明する説明図である。算出された各正規分布p(x)の平均値が移動体の推定位置の座標(x、y)となる。

Figure 0004032125
That is, in step S6, the moving body position estimation system 1 operates as follows. Equation 1 is set in advance in the moving object position estimation system 1, and the moving object position estimation system 1 uses the information on the particles sampled in step S5, and the number of parameters c in Equation 1 and the weight of each normal distribution. The average value and variance of α i and normal distribution p i (x) are calculated. Specifically, the weight α i of each normal distribution and the average value and variance of the normal distribution p i (x) are calculated based on an EM algorithm using clusters to which individual particles belong as hidden variables. The EM algorithm is a method widely known as a maximum likelihood estimation method in the case where there are hidden variables, and it has been proved that the solution converges by performing repeated calculations. On the other hand, the number c of the normal mixture distribution P (x) is incremented one by one until the value of the log likelihood of each particle in the normal mixture distribution P (x) shown in Equation 2 is decreased with the initial value set to 1. Ask. FIG. 4 (5) is an explanatory view for conceptually explaining the average value of each cluster. The average value of each calculated normal distribution p i (x) is the coordinates (x, y) of the estimated position of the moving object.
Figure 0004032125

さらに、移動体位置推定システム1は、各推定位置(x、y)における移動体が存在する確率P(x)を算出してもよい。この場合、上記算出された各パラメータを数1に代入し、上記算出された各推定位置の座標(x、y)を数1の入力情報(x)とし、P(x)を計算する。   Furthermore, the moving body position estimation system 1 may calculate the probability P (x) that there is a moving body at each estimated position (x, y). In this case, the calculated parameters are substituted into Equation 1, and the coordinates (x, y) of the calculated estimated positions are used as the input information (x) of Equation 1, and P (x) is calculated.

なお、クラスタリング方法は、上記に限られるものではなく、あらかじめ固定された数のクラスターに分類する非階層的クラスタリングでも良いし、互いの距離の近いものから順にクラスターを融合していくことにより、適当な数のグループに分割するC−Means法などの階層的クラスタリングでも良いし、クラスタリング方法は限定されない。   Note that the clustering method is not limited to the above, and may be non-hierarchical clustering that classifies a fixed number of clusters, or by combining clusters in order from the closest distance to each other. Hierarchical clustering such as a C-Means method for dividing into a large number of groups may be used, and the clustering method is not limited.

(ステップS7)ステップS6にて算出された複数の推定位置(x,y)を出力する。ステップ6にて、各推定位置(x,y)に対する移動体が存在する確率P(x)を算出した場合は、各推定位置(x,y)に確率P(x)を表示したり、推定位置(x,y)を確率P(x)の高い順に順位付けして出力してもよい。出力先はディスプレイ16でも良いし、通信インタフェイス12からネットワークを介して接続可能な情報端末(例えば携帯電話等)に出力しても良い。また、算出された複数の推定位置のすべてを出力してもよいし、確率順に上位数個を出力するようにしてもよい。   (Step S7) A plurality of estimated positions (x, y) calculated in step S6 are output. In step 6, when the probability P (x) that the moving body exists for each estimated position (x, y) is calculated, the probability P (x) is displayed at each estimated position (x, y) or estimated. Position (x, y) may be ranked and output in descending order of probability P (x). The output destination may be the display 16 or may be output from the communication interface 12 to an information terminal (for example, a mobile phone) that can be connected via a network. Also, all of the calculated plurality of estimated positions may be output, or the top several may be output in order of probability.

移動体位置推定システム1による上記移動体位置推定方法によれば、道路や建物や河川といった領域ごとに移動体が存在する確率データD2を反映させるため、移動体の位置をより高い精度で推定可能となる。また、パーティクルをクラスタリングして位置推定を行うため、各クラスター内における推定確率の分散を抑えることができ、推定精度を高めることができる。また、推定位置が唯一に出力される従来技術と比較して、本発明は推定位置として複数の候補を得ることができる。たとえば、ユーザが移動体を捜索している場合は、複数の候補地を捜索することが可能となる。   According to the mobile body position estimation method by the mobile body position estimation system 1, since the probability data D2 in which a mobile body exists for each region such as a road, a building, and a river is reflected, the position of the mobile body can be estimated with higher accuracy. It becomes. Further, since the position is estimated by clustering the particles, dispersion of the estimation probability within each cluster can be suppressed, and the estimation accuracy can be increased. Moreover, compared with the prior art in which the estimated position is output only, the present invention can obtain a plurality of candidates as estimated positions. For example, when the user is searching for a moving body, it is possible to search for a plurality of candidate sites.

(第2の実施の形態)
また、例えば移動体位置推定システム1が徘徊老人の探索など、移動体の探索に使用される場合は、出力される複数の推定位置を探索者が巡回する必要がある。そこで、移動体位置推定システム1は、複数の推定位置の巡回経路を出力するようにしても良い。巡回経路の出力は次のように行う。図8は本実施の形態の処理を説明するフローチャートである。この巡回経路の処理(ステップS6−1)は、上記第1の実施の形態の処理のステップS6の終了後に行われ、巡回経路を算出し、推定位置とともに巡回経路の出力を行うことにより実現される。第1の実施の形態と同様の部分については説明を省略する。
(Second Embodiment)
In addition, for example, when the mobile object position estimation system 1 is used for a mobile object search such as a search for an elderly person, it is necessary for the searcher to go around a plurality of estimated positions to be output. Therefore, the mobile object position estimation system 1 may output a round route of a plurality of estimated positions. The cyclic route is output as follows. FIG. 8 is a flowchart for explaining the processing of this embodiment. This cyclic route processing (step S6-1) is performed after step S6 of the processing of the first embodiment is completed, and is realized by calculating the cyclic route and outputting the cyclic route together with the estimated position. The A description of the same parts as those in the first embodiment will be omitted.

移動体位置推定装置1には、巡回経路の算出のステップS6−1よりも前までに巡回のスタート地点が入力される。スタート地点は、巡回者が任意に入力可能としてもよいし、巡回者が携帯する情報端末のGPS等の機能を用いて移動体位置計測装置から得られる位置情報を入力としても良い。移動体位置測定システム1は、スタート地点の位置情報と、上記算出した複数の推定位置やその推定位置に移動体が存在する確率から、巡回者が移動体に到達する可能性が高く、且つ、短時間に巡回できる経路を算出する。   The moving body position estimation apparatus 1 is input with the starting point of the tour before step S6-1 of the tour route calculation. The start point may be arbitrarily input by the traveler, or may be input as position information obtained from the mobile body position measurement device using a function such as GPS of an information terminal carried by the traveler. The mobile body position measurement system 1 is highly likely that a patrol person will reach the mobile body from the position information of the start point and the calculated plurality of estimated positions and the probability that the mobile body exists at the estimated positions, and A route that can be visited in a short time is calculated.

たとえば、移動体位置推定システム1は、次のように巡回経路を算出する。スタート地点から推定位置(すなわち、推定位置を表す正規混合分布P(x)を構成する各正規分布p(x)の頂点)をすべて巡回する経路を探索する。正規混合分布P(x)を構成する正規分布p(x)の数をCとすると、経路はその階乗C!だけ存在する。これらの経路を評価する関数として、下記数3を用いる。

Figure 0004032125
Routeは評価値Vを求めようとする経路における推定位置r(スタート地点を除く各地点)の集合であり、rはその推定位置の1つを示す。T(r)は、推定位置rに到達するまでの1つ前の地点(スタート地点又は推定位置)からの時間を求める関数である。P(r)は、推定位置rの1つ前の地点までの全ての推定位置における移動体が存在する確率の総和である。評価値Vの値が小さいほど、効率的な巡回経路であるといえる。 For example, the moving body position estimation system 1 calculates a traveling route as follows. A route that circulates all the estimated positions from the start point (that is, the vertices of each normal distribution p i (x) constituting the normal mixture distribution P (x) representing the estimated position) is searched. If the number of normal distributions p i (x) constituting the normal mixed distribution P (x) is C, the path is the factorial C! Only exist. The following formula 3 is used as a function for evaluating these paths.
Figure 0004032125
Route is a set of estimated positions r (points excluding the start point) on the route for which the evaluation value V is to be obtained, and r indicates one of the estimated positions. T (r) is a function for obtaining the time from the previous point (start point or estimated position) until the estimated position r is reached. P (r) is the sum of the probabilities that there are moving objects at all estimated positions up to the point immediately before the estimated position r. It can be said that the smaller the evaluation value V is, the more efficient the traveling route is.

図9は、数3の一実施例をグラフで表した図である。図9の横軸は時間、縦軸はその移動体が存在しない確率(数3における1−P(x))、すなわち、探索者が被探索者を発見できない確率を示している。ここでは、正規混合分布を構成する正規分布の数(すなわち推定位置の数)C=3を例に説明する。すべての経路について次の処理を行う。スタート地点sから一つ目の正規分布の頂点(推定位置)p1に達するときに必要な時間と、その到達した頂点(推定位置)p1において移動体を発見できない確率を算出する。さらに、一つ目の頂点(推定位置)p1から次の頂点(推定位置)p2に達するときに必要な時間と、その到達した頂点p2における移動体を発見できない確率を算出する。このように、スタート地点から巡回する経由地点順に、時間と移動体を発見できない確率を算出する。そして、横軸を時間、縦軸を移動体を発見できない確率として、算出された各頂点(推定位置)の時間と移動体を発見できない確率をプロットし、その積分値(図9の斜線部分の面積)の大きさが評価値Vとなる。これにより、移動体を発見できる確率と巡回時間とを考慮した評価が可能となる。巡回経路の出力に際しては、評価値Vが低い経路のみを提示してもよいし、評価値順に順位付けしてすべての巡回経路を出力しても良い。   FIG. 9 is a graph showing an example of Equation (3). In FIG. 9, the horizontal axis indicates time, and the vertical axis indicates the probability that the moving object does not exist (1-P (x) in Equation 3), that is, the probability that the searcher cannot find the search target. Here, the number of normal distributions constituting the normal mixed distribution (that is, the number of estimated positions) C = 3 will be described as an example. The following processing is performed for all routes. The time required to reach the first normal distribution vertex (estimated position) p1 from the start point s and the probability that a moving object cannot be found at the reached vertex (estimated position) p1 are calculated. Further, the time required to reach the next vertex (estimated position) p2 from the first vertex (estimated position) p1 and the probability that a moving object at the reached vertex p2 cannot be found are calculated. In this way, the time and the probability that the moving object cannot be found are calculated in the order of the waypoints that circulate from the start point. Then, with the horizontal axis representing time and the vertical axis representing the probability that a moving body cannot be found, the calculated time of each vertex (estimated position) and the probability that a moving body cannot be found are plotted, and the integrated value (the hatched portion in FIG. 9). The size of (area) is the evaluation value V. This makes it possible to evaluate in consideration of the probability of finding a moving object and the traveling time. When outputting the cyclic route, only the route with the low evaluation value V may be presented, or all the cyclic routes may be output in the order of the evaluation value.

一例として、図10(a)のように、移動体位置推定の結果、a,b,cの三つの推定位置が得られた場合の巡回経路の算出を説明する。スタート地点は座標(0,0)であり、推定位置aは座標(5,0)、bは座標(4,3)、cは座標(0,5)の位置であり、各推定位置における移動体が存在する確率がaは0.5,bは0.1,cは0.4であったとする。すべての巡回経路について上記評価値Vを算出した結果を図10(b)に示す。a−b−cの順で巡回する経路が最も良い評価値であるという結果が得られ、有効な巡回経路としてa−b−cを出力される。   As an example, as shown in FIG. 10A, calculation of a cyclic route when three estimated positions a, b, and c are obtained as a result of mobile object position estimation will be described. The starting point is the coordinates (0, 0), the estimated position a is the coordinates (5, 0), b is the coordinates (4, 3), c is the position of the coordinates (0, 5), and the movement at each estimated position The probability that a body exists is 0.5 for a, 0.1 for b, and 0.4 for c. FIG. 10B shows the result of calculating the evaluation value V for all the traveling routes. The result that the route that goes around in the order of abc is the best evaluation value is obtained, and abc is outputted as an effective tour route.

これにより、探索者は、出力された複数の推定位置を効率的に巡回できる経路の情報を得ることができる。たとえば、徘徊老人の探索においては、効率的に巡回することにより、より迅速に徘徊老人を見つけ出し、安全を確保することができる。   As a result, the searcher can obtain information on a route through which the plurality of output estimated positions can be efficiently circulated. For example, in searching for an elderly person, it is possible to find the elderly person more quickly and ensure safety by patroling efficiently.

なお、本発明は、上記実施の形態に限られるものではなく、本発明から逸脱しない範囲内で変更が可能である。たとえば、確率データを用いて領域ごとに移動体が存在する確率を反映させる技術(ステップS2)は、位置推定にクラスタリングを用いることなく、従来のようにパラメータの最尤推定を行う技術に適用しても良い。これによっても、領域ごとの確率が反映され、従来よりも高精度に位置推定が行われる。この場合、移動体位置推定システム1は以下のように処理を行う。図11は、そのフローチャートである。移動体位置推定システムは、上記ステップ1からステップ5を行った後、サンプリング結果からパラメータ最尤推定を行って確立分布を算出し、得られた確率分布から推定位置を算出し(ステップS6)、推定位置を出力する(ステップS7)。パラメータ最尤推定や、確率分布からの推定位置の算出は、位置計測装置の種類に応じて行えばよい。例えばGPSであればパラメータ最尤推定を行って正規分布を算出し、その正規分布の平均値を推定位置として出力する。また、PHSの場合は、ステップ4において、受信した基地局が電波を受信できる範囲内において移動体が存在する確率を0より大の一定値(0より大)とし、その他の範囲は確率を0とする一様分布データに変換される。ステップ4以外は、PHSであってもGPSであっても同じ処理である。   Note that the present invention is not limited to the above-described embodiment, and can be modified without departing from the present invention. For example, the technique (step S2) that reflects the probability that a moving object exists for each region using probability data is applied to a conventional technique that performs maximum likelihood estimation of parameters without using clustering for position estimation. May be. This also reflects the probability for each region, and position estimation is performed with higher accuracy than in the past. In this case, the moving body position estimation system 1 performs processing as follows. FIG. 11 is a flowchart thereof. The mobile object position estimation system performs step 1 to step 5 above, then performs parameter maximum likelihood estimation from the sampling result to calculate the probability distribution, calculates the estimated position from the obtained probability distribution (step S6), The estimated position is output (step S7). The parameter maximum likelihood estimation and the calculation of the estimated position from the probability distribution may be performed according to the type of the position measuring device. For example, in the case of GPS, parameter maximum likelihood estimation is performed to calculate a normal distribution, and an average value of the normal distribution is output as an estimated position. In the case of PHS, in step 4, the probability that a moving body exists within the range in which the received base station can receive radio waves is set to a constant value greater than 0 (greater than 0), and the probability of other ranges is 0. Is converted into uniform distribution data. Except for step 4, the process is the same for both PHS and GPS.

また、本発明は、様々な分野に応用可能である。たとえば、GPS機能付き情報携帯端末等の移動体を人に携帯させ、徘徊する高齢者や遭難者を捜索したり、登下校中の児童の位置を把握して安全を確保するのに用いることができる。また、移動体を自動車等に搭載すれば、タクシーや配送車等の位置管理などにも応用可能である。   In addition, the present invention can be applied to various fields. For example, a mobile object such as an information portable terminal with GPS function can be carried by a person to search for an elderly person or a distressed person who is hesitant, or used to ascertain the position of a child who is going to and from school to ensure safety. it can. In addition, if the moving body is mounted on an automobile or the like, it can be applied to the position management of a taxi or a delivery car.

本発明の第1の実施の形態による移動体位置推定システムの一構成を示すブロック図である。It is a block diagram which shows one structure of the moving body position estimation system by the 1st Embodiment of this invention. 地図情報と確率データを概念的に説明する説明図である。It is explanatory drawing which illustrates map information and probability data notionally. 上記実施の形態における移動体位置推定システムを用いた移動体位置推定方法を説明するフローチャートである。It is a flowchart explaining the mobile body position estimation method using the mobile body position estimation system in the said embodiment. 上記実施の形態の位置推定処理を概念的に説明する説明図である。It is explanatory drawing which illustrates notionally the position estimation process of the said embodiment. 上記実施の形態において地図情報に対する確率データの設定を説明する説明図である。It is explanatory drawing explaining the setting of the probability data with respect to map information in the said embodiment. GPSの位置情報を用いたときの確率分布データと、地図情報上のパーティクルの母集団と、抽出されたパーティクルとを説明する説明図である。It is explanatory drawing explaining the probability distribution data when using the positional information of GPS, the population of the particle on map information, and the extracted particle. パーティクルフィルタにより抽出されたパーティクルの確率分布を正規混合分布によって表す例を示す説明図である。It is explanatory drawing which shows the example which represents the probability distribution of the particle extracted by the particle filter by normal mixture distribution. 第2の実施の形態の処理を説明するフローチャートである。It is a flowchart explaining the process of 2nd Embodiment. 上記実施の形態の巡回経路の算出に用いる関数を説明する説明図である。It is explanatory drawing explaining the function used for calculation of the cyclic route of the said embodiment. 上記実施の形態の巡回経路の算出の例を示す説明図である。It is explanatory drawing which shows the example of calculation of the cyclic route of the said embodiment. 他の実施の形態の処理を説明するフローチャートである。It is a flowchart explaining the process of other embodiment.

符号の説明Explanation of symbols

1 移動体位置推定システム
19c 地図データ
D1 基本データ
D2 確率データ
1 Mobile object position estimation system 19c Map data D1 Basic data D2 Probability data

Claims (8)

移動体位置計測装置から移動体の位置情報を取得すると、当該位置情報を用いて地図データの座標系に分散するパーティクルをサンプリングするパーティクルフィルタを備え、当該パーティクルフィルタによりサンプリングしたパーティクルから当該移動体の位置を推定する移動体位置推定システムにおいて、
当該地図データの座標系上の各領域に当該移動体が存在する確率を示す確率データを有し、当該確率データに従って当該地図データの座標系にパーティクルを分散させる手段を備えることを特徴とする移動体位置推定システム。
When the position information of the moving object is acquired from the moving object position measuring device, a particle filter that samples particles dispersed in the coordinate system of the map data using the position information is provided, and the moving object is detected from the particles sampled by the particle filter. In a mobile position estimation system for estimating a position,
The movement having probability data indicating the probability that the moving body exists in each area on the coordinate system of the map data, and having means for dispersing particles in the coordinate system of the map data according to the probability data Body position estimation system.
前記確率データは、前記各領域がカテゴリ分けされており、当該カテゴリごとに前記確率を設定可能であることを特徴とする請求項1記載の移動体位置推定システム。   The mobile object position estimation system according to claim 1, wherein the probability data has the respective areas classified into categories, and the probability can be set for each category. 前記確率データは、前記地図データに含まれる各領域を識別する識別情報を用いて、前記各領域がカテゴリ分けされていることを特徴とする請求項2記載の移動体位置推定システム。   3. The moving object position estimation system according to claim 2, wherein the probability data is classified into categories using identification information for identifying each area included in the map data. 前記サンプリングされたパーティクルをクラスタリングし、各クラスターの平均値を前記移動体の推定位置として出力する手段を備えることを特徴とする請求項1乃至請求項3のいずれか1項に記載の移動体位置システム。   The moving body position according to any one of claims 1 to 3, further comprising means for clustering the sampled particles and outputting an average value of each cluster as an estimated position of the moving body. system. コンピュータが、移動体位置計測装置から移動体の位置情報を取得すると、パーティクルフィルタにより当該位置情報を用いて地図データの座標系に分散するパーティクルをサンプリングし、当該サンプリングしたパーティクルから当該移動体の位置を推定する移動体位置推定方法において、
当該地図データの座標系上の各領域に当該移動体が存在する確率を示す確率データを参照し、当該確率データに従って当該地図データの座標系にパーティクルを分散させるステップを備えることを特徴とする移動体位置推定方法。
When the computer acquires the position information of the moving object from the moving object position measuring device, the particle filter samples the particles dispersed in the coordinate system of the map data using the position information, and the position of the moving object is determined from the sampled particles. In the moving object position estimation method for estimating
The movement comprising the step of referring to probability data indicating the probability that the moving object exists in each area on the coordinate system of the map data, and dispersing particles in the coordinate system of the map data according to the probability data Body position estimation method.
前記確率データは、前記各領域がカテゴリ分けされており、当該カテゴリごとに前記確率を設定可能であることを特徴とする請求項5記載の移動体位置推定方法。   6. The method according to claim 5, wherein in the probability data, each region is classified into categories, and the probability can be set for each category. 前記地図データには前記各領域を識別する識別情報が含まれており、前記確率データは、当該識別情報を用いて前記各領域がカテゴリ分けされていることを特徴とする請求項6記載の移動体位置推定方法。   The movement according to claim 6, wherein the map data includes identification information for identifying each area, and the probability data includes the areas classified into categories using the identification information. Body position estimation method. 前記サンプリングされたパーティクルをクラスタリングし、各クラスターの平均値を前記移動体の推定位置として出力するステップを備えることを特徴とする請求項5乃至請求項7のいずれか1項に記載の移動体位置方法。   The moving body position according to any one of claims 5 to 7, further comprising a step of clustering the sampled particles and outputting an average value of each cluster as an estimated position of the moving body. Method.
JP2006185198A 2006-07-05 2006-07-05 Mobile object position estimation system and mobile object position estimation method Active JP4032125B1 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2006185198A JP4032125B1 (en) 2006-07-05 2006-07-05 Mobile object position estimation system and mobile object position estimation method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2006185198A JP4032125B1 (en) 2006-07-05 2006-07-05 Mobile object position estimation system and mobile object position estimation method

Publications (2)

Publication Number Publication Date
JP4032125B1 true JP4032125B1 (en) 2008-01-16
JP2008014742A JP2008014742A (en) 2008-01-24

Family

ID=39071900

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2006185198A Active JP4032125B1 (en) 2006-07-05 2006-07-05 Mobile object position estimation system and mobile object position estimation method

Country Status (1)

Country Link
JP (1) JP4032125B1 (en)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103604426A (en) * 2013-12-02 2014-02-26 苏州大学张家港工业技术研究院 Estimation method and apparatus for poses of mobile robot
CN103827632A (en) * 2012-09-06 2014-05-28 株式会社东芝 Position detection device, position detection method and position detection program

Families Citing this family (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR100855469B1 (en) * 2006-09-13 2008-09-01 삼성전자주식회사 Pose estimation apparatus used in a mobile robot and method thereof
KR100926783B1 (en) * 2008-02-15 2009-11-13 한국과학기술연구원 Method for self-localization of a robot based on object recognition and environment information around the recognized object
JP4991595B2 (en) 2008-02-21 2012-08-01 株式会社東芝 Tracking system using particle filter
KR100951321B1 (en) 2008-02-27 2010-04-05 아주대학교산학협력단 Method of object tracking in 3D space based on particle filter using acoustic sensors
JP5515647B2 (en) * 2009-11-05 2014-06-11 独立行政法人産業技術総合研究所 Positioning device
JP5396311B2 (en) * 2010-03-02 2014-01-22 日本電信電話株式会社 Behavior prediction apparatus, method, and program
US8731829B2 (en) 2010-05-31 2014-05-20 Nec Corporation Flow line detection system, flow line detection method, and flow line detection program
JP5640867B2 (en) 2011-03-31 2014-12-17 富士通株式会社 Information processing apparatus, information processing method, information processing program, and electronic device
US8583400B2 (en) * 2011-05-13 2013-11-12 Google Inc. Indoor localization of mobile devices
JP5978747B2 (en) * 2012-05-09 2016-08-24 富士通株式会社 Positioning system, positioning method, and program

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP4251545B2 (en) * 2003-07-11 2009-04-08 独立行政法人科学技術振興機構 Route planning system for mobile robot

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103827632A (en) * 2012-09-06 2014-05-28 株式会社东芝 Position detection device, position detection method and position detection program
CN103604426A (en) * 2013-12-02 2014-02-26 苏州大学张家港工业技术研究院 Estimation method and apparatus for poses of mobile robot

Also Published As

Publication number Publication date
JP2008014742A (en) 2008-01-24

Similar Documents

Publication Publication Date Title
JP4032125B1 (en) Mobile object position estimation system and mobile object position estimation method
JP4032126B1 (en) Mobile object position estimation system and mobile object position estimation method
EP3241370B1 (en) Analyzing semantic places and related data from a plurality of location data reports
US11562168B2 (en) Clustering for K-anonymity in location trajectory data
JP4936647B2 (en) Method for determining the approximate location of a device from environmental signals
US9355063B2 (en) Parking lot detection using probe data
US10546043B1 (en) Triangulation for K-anonymity in location trajectory data
EP2712488B1 (en) Prediction of indoor level and location using a three stage process
KR101099204B1 (en) Systems for determining the approximate location of a device from ambient signals
RU2356091C2 (en) Application of device approximate location defined on basis of external signals
CN104584094B (en) Location estimation method and system
CN110793531B (en) Road matching method and device and readable storage medium
CN112148821B (en) City mixed occupation space calculation method and system
St Reconstructing trajectories from sparse call detail records
Au Rss-based wlan indoor positioning and tracking system using compressive sensing and its implementation on mobile devices
JP5358534B2 (en) Action path estimation method and program
Keler et al. Is there a relationship between complicated crossings and frequently visited locations? A case study with boro taxis and OSM in NYC
Jaegal et al. Measuring the structural similarity of network time prisms using temporal signatures with graph indices
Arregui et al. Impact of the road network configuration on map‐matching algorithms for FCD in urban environments
CN117407477B (en) Geographic information data evolution recognition processing method, system and storage medium
Ye et al. Development of a Digital Accident Hotspot Map for ADAS Applications Using Geospatial Methods in GIS
Hindmarch et al. Movement model enhanced RSS localisation
EP3944137A1 (en) Positioning method and positioning apparatus
EP4224116A1 (en) Utilizing models to evaluate geolocation estimate quality without independent test data
Shah Radio map construction using Kriging and functional data analysis for indoor localization

Legal Events

Date Code Title Description
A521 Written amendment

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20061207

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20061201

A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20070802

A871 Explanation of circumstances concerning accelerated examination

Free format text: JAPANESE INTERMEDIATE CODE: A871

Effective date: 20070802

A975 Report on accelerated examination

Free format text: JAPANESE INTERMEDIATE CODE: A971005

Effective date: 20070821

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

R150 Certificate of patent or registration of utility model

Free format text: JAPANESE INTERMEDIATE CODE: R150