"それぞれが所定の周期で識別情報を無線で発信する、または、所定の周期で無線で信号を発信する移動体から該信号を受信する、分散して配置される複数の通信機と、前記通信機が発信する識別情報を受信する移動体、または、前記所定の周期で無線で信号を発信する移動体と、所定の期間、所定の領域内の複数の地点で前記移動体が前記通信機から受信した前記識別情報および異なる2つの前記識別情報を受信した地点が同一か否かを表す情報の集合、あるいは、前記移動体が所定の期間、所定の領域内を移動しながら前記所定の周期で発信する信号を受信した前記通信機の識別情報および異なる前記通信機が受信した信号が同一の地点から発信されたか否かを表す情報の集合、である通信機ID集合を取得する収集手段と、前記通信機ID集合に含まれる前記識別情報ごとの受信回数、および、前記2つの識別情報の組ごとの前記識別情報の受信もしくは発信した地点が同一である受信回数、を計数した受信情報を生成する受信情報生成手段と、前記通信機ごとの発信の到達範囲もしくは受信可能な前記移動体の発信範囲を示す通信範囲情報を取得する範囲情報取得手段と、前記受信情報および前記通信範囲情報から、前記複数の通信機に含まれる2つの前記通信機の間の通信機間距離を算出する距離算出手段と、前記複数の通信機のうち、その位置を頂点とする三角形を構成できる3つの通信機それぞれの前記識別情報および設置座標を含む通信機情報を取得する手段と、前記通信機間距離と前記通信機情報から、前記通信機それぞれの設置座標を算出する座標算出手段と、を備えることを特徴とする位置計測システム。前記通信範囲情報は、前記複数の通信機それぞれの前記到達範囲もしくは前記発信範囲を示すエリア半径を含み、前記距離算出手段は、前記通信機の前記到達範囲もしくは前記発信範囲を、前記通信機を中心とする前記エリア半径の円とみなして、前記通信機間距離を算出することを特徴とする請求項1に記載の位置計測システム。前記距離算出手段は、前記2つの通信機それぞれの前記エリア半径、および、前記到達範囲もしくは前記発信範囲の重複面積から前記通信間距離を算出することを特徴とする請求項2に記載の位置計測システム。前記距離算出手段は、前記2つの通信機それぞれの前記到達範囲もしくは前記発信範囲の面積に、前記2つの通信機の識別情報ごとの受信回数に対する前記受信もしくは発信した地点が同一である受信回数の比率をかけることによって、前記重複面積を算出することを特徴とする請求項3に記載の位置計測システム。前記複数の通信機それぞれの前記エリア半径は同一であることを特徴とする請求項2ないし4のいずれか1項に記載の位置計測システム。前記座標算出手段は、前記通信機間距離から余弦定理を用いて、2つの前記通信機に対する他の1つの前記通信機の相対座標を計算することを特徴とする請求項1ないし5のいずれか1項に記載の位置計測システム。前記座標算出手段は、前記通信機のいずれか1つを原点とし、他のいずれか1つと前記原点を通る直線を座標軸の1つとする仮想座標系を設定して、前記通信機間距離から前記通信機の前記仮想座標系における仮想座標を算出し、前記通信機情報と前記仮想座標から前記通信機の座標を算出することを特徴とする請求項6に記載の位置計測システム。それぞれが所定の周期で識別情報を無線で発信する分散して配置される複数の通信機から、所定の期間、所定の領域内の複数の地点で受信した識別情報および異なる2つの識別情報を受信した地点が同一か否かを表す情報の集合、あるいは、所定の周期で無線で信号を発信する移動体が所定の期間、所定の領域内を移動しながら発信する信号を受信した、分散して配置される複数の通信機の識別情報および異なる前記通信機が受信した信号が同一の地点から発信されたか否かを表す情報の集合、である通信機ID集合を取得する収集手段と、前記通信機ID集合に含まれる前記識別情報ごとの受信回数、および、前記2つの識別情報の組ごとの前記識別情報の受信もしくは発信した地点が同一である受信回数、を計数した受信情報を生成する受信情報生成手段と、前記通信機ごとの発信の到達範囲もしくは受信可能な前記移動体の発信範囲を示す通信範囲情報を取得する範囲情報取得手段と、前記受信情報および前記通信範囲情報から、前記複数の通信機に含まれる2つの前記通信機の間の通信機間距離を算出する距離算出手段と、前記複数の通信機のうち、その位置を頂点とする三角形を構成できる3つの通信機それぞれの前記識別情報および設置座標を含む通信機情報を取得する手段と、前記通信機間距離と前記通信機情報から、前記通信機それぞれの設置座標を算出する座標算出手段と、を備えることを特徴とする位置計測装置。それぞれが所定の周期で識別情報を無線で発信する分散して配置される複数の通信機から、所定の期間、所定の領域内の複数の地点で受信した識別情報および異なる2つの識別情報を受信した地点が同一か否かを表す情報の集合、あるいは、所定の周期で無線で信号を発信する移動体が所定の期間、所定の領域内を移動しながら発信する信号を受信した、分散して配置される複数の通信機の識別情報および異なる前記通信機が受信した信号が同一の地点から発信されたか否かを表す情報の集合、である通信機ID集合を取得する収集ステップと、前記通信機ID集合に含まれる前記識別情報ごとの受信回数、および、前記2つの識別情報の組ごとの前記識別情報の受信もしくは発信した地点が同一である受信回数、を計数した受信情報を生成する受信情報生成ステップと、前記通信機ごとの発信の到達範囲もしくは受信可能な前記移動体の発信範囲を示す通信範囲情報を取得する範囲情報取得ステップと、前記受信情報および前記通信範囲情報から、前記複数の通信機に含まれる2つの前記通信機の間の通信機間距離を算出する距離算出ステップと、前記複数の通信機のうち、その位置を頂点とする三角形を構成できる3つの通信機それぞれの前記識別情報および設置座標を含む通信機情報を取得するステップと、前記通信機間距離と前記通信機情報から、前記通信機それぞれの設置座標を算出する座標算出ステップと、を備えることを特徴とする位置計測方法。コンピュータに、それぞれが所定の周期で識別情報を無線で発信する分散して配置される複数の通信機から、所定の期間、所定の領域内の複数の地点で受信した識別情報および異なる2つの識別情報を受信した地点が同一か否かを表す情報の集合、あるいは、所定の周期で無線で信号を発信する移動体が所定の期間、所定の領域内を移動しながら発信する信号を受信した、分散して配置される複数の通信機の識別情報および異なる前記通信機が受信した信号が同一の地点から発信されたか否かを表す情報の集合、である通信機ID集合を取得する収集ステップと、前記通信機ID集合に含まれる前記識別情報ごとの受信回数、および、前記2つの識別情報の組ごとの前記識別情報の受信もしくは発信した地点が同一である受信回数、を計数した受信情報を生成する受信情報生成ステップと、前記通信機ごとの発信の到達範囲もしくは受信可能な前記移動体の発信範囲を示す通信範囲情報を取得する範囲情報取得ステップと、前記受信情報および前記通信範囲情報から、前記複数の通信機に含まれる2つの前記通信機の間の通信機間距離を算出する距離算出ステップと、前記複数の通信機のうち、その位置を頂点とする三角形を構成できる3つの通信機それぞれの前記識別情報および設置座標を含む通信機情報を取得するステップと、前記通信機間距離と前記通信機情報から、前記通信機それぞれの設置座標を算出する座標算出ステップと、を実行させることを特徴とするプログラム。特許文献4に示されるような赤外線センサなどの信号受信の有無しか用いることができないセンサIN/OUT型センサでは、個別の発信機の位置情報をあらかじめ全て登録した上で、受信機のIN/OUTの判定から位置検出を行うため、発信機の設置コストがかかるという、第1の課題がある。特許文献1の位置決定方法では、受信信号の信号強度を距離を表す指標として用いている。また、特許文献2の位置検出システムでは、測距用信号を用いている。非特許文献1に示されるような電波強度型センサにおいては、電波強度を使って発信機からの距離を推定できるため、より高精度な位置計測が可能であるが、この手法はIN/OUT型センサには適用することができないという、第2の課題がある。本発明は、上述のような事情に鑑みてなされたもので、IN/OUT型のセンサを用いる方式において、センサーノードの設置コストを削減することを目的とする。本発明の第1の観点に係る位置計測システムは、それぞれが所定の周期で識別情報を無線で発信する、または、所定の周期で無線で信号を発信する移動体から該信号を受信する、分散して配置される複数の通信機と、前記通信機が発信する識別情報を受信する移動体、または、前記所定の周期で無線で信号を発信する移動体と、所定の期間、所定の領域内の複数の地点で前記移動体が前記通信機から受信した前記識別情報および異なる2つの前記識別情報を受信した地点が同一か否かを表す情報の集合、あるいは、前記移動体が所定の期間、所定の領域内を移動しながら前記所定の周期で発信する信号を受信した前記通信機の識別情報および異なる前記通信機が受信した信号が同一の地点から発信されたか否かを表す情報の集合、である通信機ID集合を取得する収集手段と、前記通信機ID集合に含まれる前記識別情報ごとの受信回数、および、前記2つの識別情報の組ごとの前記識別情報の受信もしくは発信した地点が同一である受信回数、を計数した受信情報を生成する受信情報生成手段と、前記通信機ごとの発信の到達範囲もしくは受信可能な前記移動体の発信範囲を示す通信範囲情報を取得する範囲情報取得手段と、前記受信情報および前記通信範囲情報から、前記複数の通信機に含まれる2つの前記通信機の間の通信機間距離を算出する距離算出手段と、前記複数の通信機のうち、その位置を頂点とする三角形を構成できる3つの通信機それぞれの前記識別情報および設置座標を含む通信機情報を取得する手段と、前記通信機間距離と前記通信機情報から、前記通信機それぞれの設置座標を算出する座標算出手段と、を備えることを特徴とする。本発明の第2の観点に係る位置計測装置は、それぞれが所定の周期で識別情報を無線で発信する分散して配置される複数の通信機から、所定の期間、所定の領域内の複数の地点で受信した識別情報および異なる2つの識別情報を受信した地点が同一か否かを表す情報の集合、あるいは、所定の周期で無線で信号を発信する移動体が所定の期間、所定の領域内を移動しながら発信する信号を受信した、分散して配置される複数の通信機の識別情報および異なる前記通信機が受信した信号が同一の地点から発信されたか否かを表す情報の集合、である通信機ID集合を取得する収集手段と、前記通信機ID集合に含まれる前記識別情報ごとの受信回数、および、前記2つの識別情報の組ごとの前記識別情報の受信もしくは発信した地点が同一である受信回数、を計数した受信情報を生成する受信情報生成手段と、前記通信機ごとの発信の到達範囲もしくは受信可能な前記移動体の発信範囲を示す通信範囲情報を取得する範囲情報取得手段と、前記受信情報および前記通信範囲情報から、前記複数の通信機に含まれる2つの前記通信機の間の通信機間距離を算出する距離算出手段と、前記複数の通信機のうち、その位置を頂点とする三角形を構成できる3つの通信機それぞれの前記識別情報および設置座標を含む通信機情報を取得する手段と、前記通信機間距離と前記通信機情報から、前記通信機それぞれの設置座標を算出する座標算出手段と、を備えることを特徴とする。本発明の第3の観点に係る位置計測方法は、それぞれが所定の周期で識別情報を無線で発信する分散して配置される複数の通信機から、所定の期間、所定の領域内の複数の地点で受信した識別情報および異なる2つの識別情報を受信した地点が同一か否かを表す情報の集合、あるいは、所定の周期で無線で信号を発信する移動体が所定の期間、所定の領域内を移動しながら発信する信号を受信した、分散して配置される複数の通信機の識別情報および異なる前記通信機が受信した信号が同一の地点から発信されたか否かを表す情報の集合、である通信機ID集合を取得する収集ステップと、前記通信機ID集合に含まれる前記識別情報ごとの受信回数、および、前記2つの識別情報の組ごとの前記識別情報の受信もしくは発信した地点が同一である受信回数、を計数した受信情報を生成する受信情報生成ステップと、前記通信機ごとの発信の到達範囲もしくは受信可能な前記移動体の発信範囲を示す通信範囲情報を取得する範囲情報取得ステップと、前記受信情報および前記通信範囲情報から、前記複数の通信機に含まれる2つの前記通信機の間の通信機間距離を算出する距離算出ステップと、前記複数の通信機のうち、その位置を頂点とする三角形を構成できる3つの通信機それぞれの前記識別情報および設置座標を含む通信機情報を取得するステップと、前記通信機間距離と前記通信機情報から、前記通信機それぞれの設置座標を算出する座標算出ステップと、を備えることを特徴とする。本発明の第4の観点に係るプログラムは、コンピュータに、それぞれが所定の周期で識別情報を無線で発信する分散して配置される複数の通信機から、所定の期間、所定の領域内の複数の地点で受信した識別情報および異なる2つの識別情報を受信した地点が同一か否かを表す情報の集合、あるいは、所定の周期で無線で信号を発信する移動体が所定の期間、所定の領域内を移動しながら発信する信号を受信した、分散して配置される複数の通信機の識別情報および異なる前記通信機が受信した信号が同一の地点から発信されたか否かを表す情報の集合、である通信機ID集合を取得する収集ステップと、前記通信機ID集合に含まれる前記識別情報ごとの受信回数、および、前記2つの識別情報の組ごとの前記識別情報の受信もしくは発信した地点が同一である受信回数、を計数した受信情報を生成する受信情報生成ステップと、前記通信機ごとの発信の到達範囲もしくは受信可能な前記移動体の発信範囲を示す通信範囲情報を取得する範囲情報取得ステップと、前記受信情報および前記通信範囲情報から、前記複数の通信機に含まれる2つの前記通信機の間の通信機間距離を算出する距離算出ステップと、前記複数の通信機のうち、その位置を頂点とする三角形を構成できる3つの通信機それぞれの前記識別情報および設置座標を含む通信機情報を取得するステップと、前記通信機間距離と前記通信機情報から、前記通信機それぞれの設置座標を算出する座標算出ステップと、を実行させることを特徴とする。本発明によれば、IN/OUT型のセンサを用いて基点となる3個の発信機の位置情報を入力することによって、全ての発信機の位置情報を推定することができる。本発明の実施の形態の構成を説明する前に、まず実施の形態の原理を説明する。本発明の実施の形態では、例えばRFタグリーダのような、それぞれの識別情報を所定の周期で発信する複数の発信機が、分散して配置される領域を想定する。複数の発信機のうち、発信機の位置を頂点として三角形を構成する3つの発信機の座標通信機情報が与えられる。発信機から発信される識別情報の到達範囲識別情報を受信できる範囲は、発信機を中心とする円であると仮定し、到達範囲の半径が与えられる。複数の発信機が分散して配置される領域に分布する複数の地点で、発信機から発信される識別情報を受信する。複数の地点で受信した識別情報を収集し、識別情報ごとの受信地点の数受信回数と、任意の異なる2つの識別情報の対ごとに、その2つの識別情報を同一の地点で受信した回数同時受信回数を計数する。任意の2つの発信機識別情報について、識別情報ごとの受信回数と同時受信回数から、それぞれの到達範囲の面積と、到達範囲が重複する範囲の面積重複エリア面積の比を算出する。到達範囲の面積と重複エリア面積の比、および到達範囲の半径から2つの発信機の距離発信機間距離を算出する。重複エリア面積が0を超え、かつ、重複エリア面積がどちらの発信機でもその到達範囲の100未満である2つの発信機について距離を算出できる。そして、三角形を構成する3つの発信機の座標と発信機間距離から、三角形の余弦定理を用いて各発信機の座標を計算する。いずれの発信機も他の2つ以上の発信機との発信機間距離が算出できて、発信機をノード、発信機間距離をエッジとするグラフを考えたときに、複数の発信機が1つの単連結のグラフで表されれば、複数の発信機全部の相対座標本発明では仮想座標というを計算できる。仮想座標は、領域に固定された座標に対して、平行移動、領域面内の回転、および、鏡像表と裏の関係の自由度がある。複数の発信機のうち、三角形を構成する3つの発信機の座標が与えられれば、これらの自由度は一意に決められる。発信機の座標を計算できれば、あとは、移動体RFタグが1つの地点で発信機RFタグリーダなどから受信した識別情報、または、移動体が1つの地点で発信した信号を受信した発信機RFタグリーダなどの識別情報、を取得し、識別情報に対応する発信機の座標の重心点を算出することによって、移動体の位置を推定できる。上述の発信機と受信機を入れ替えて、領域内の複数の地点から発信する信号を、分散して配置された複数の受信機で受信し、受信した信号に受信機の識別情報をつけたものを、識別情報ごとに計数して受信回数としてもよい。その場合、同時受信回数は、2つの受信機で同じ信号を受信した回数である。本発明では、発信機と受信機を総称して通信機という。また、発信機が発信する識別情報を受信する装置、または、領域内の複数の地点から信号を発信する装置を移動体という。本発明の実施の形態では、通信機が発信機の場合を主に説明する。なお、実施の形態では、分散して配置される複数の発信機から、その識別情報を受信する装置移動体をID受信機という。以下、本発明の実施の形態について図面を参照して詳細に説明する。なお図中、同一または同等の部分には同一の符号を付す。実施の形態1図1は、本発明の実施の形態1に係る位置計測システムの構成例を示すブロック図である。位置計測システム1は、複数の発信機3、ID受信機4および位置計測装置2から構成される。発信機3は、それぞれ固有の識別情報を所定の周期で発信する。複数の発信機3は、ID受信機4が移動しうる領域を、その識別情報の到達範囲の集合で被覆するように分散して配置される。ID受信機4は、所定の領域内の複数の地点で、発信機3が発信する識別情報発信機IDを受信する。ID受信機4は、受信した識別情報に受信した時刻をつけて記憶する。または、同一の地点で受信した複数の識別情報を、受信した地点ごとに関連づけてまとめて記憶する。ID受信機4は、記憶した識別情報を、位置計測装置2に送信する。ID受信機4は、識別情報を受信するごとに位置計測装置2に送信して、ID受信機4内に識別情報を記憶しておかなくてもよい。いずれの場合でも、位置計測装置2は、受信した時刻または受信した地点ごとの関連づけから、受信した地点ごとに受信した識別情報を記憶する。位置計測装置2は、受信情報生成部11、受信情報データベース受信情報DB12、発信機間距離算出部13、発信機位置算出部15、発信機情報データベース発信機情報DB14、および、位置推定部16を備える。受信情報生成部11は、ID受信機4から、ID受信機4が受信した識別情報を受信取得して、発信機IDごとに受信した地点の数単独受信回数を計数する。また、異なる2つの識別情報の組対ごとに、その2つの識別情報を同一の地点で受信した回数同時受信回数を計数する。図2は、実施の形態1における発信機ID集合ごとの受信回数の例を示す図である。図2の発信機ID集合は、それに含まれる発信機IDを同一の地点で受信したことを表す。受信回数は、発信ID集合を受信した地点の数である。例えば、1行目の発信機ID集合=S1の場合、S1のみを受信した地点の数が18である。また、2行目は、発信機ID集合={S1,S2}であり、S1とS2を同一の地点で受信した地点の数が10であることを示す。図3は、実施の形態1において受信情報から単独受信回数と同時受信回数を計数する例を示す。図3の発信機ID集合と受信回数は、図2と同じである。例えば、発信機IDがS1とS2の場合、S1の単独受信回数は、S1が発信機ID集合に現れる受信回数を計数して、40になる。S2の単独受信回数は、S2が発信機ID集合に現れる受信回数を計数して、44である。そして、S1とS2の同時受信回数は、{S1,S2}が含まれる発信機ID集合の受信回数を計数して、14になる。図4は、実施の形態1における単独受信回数の例を示す図である。図5は、実施の形態1における同時受信回数の例を示す図である。図3で説明したように、受信情報生成部11は、発信機IDごとに、発信機IDが発信機ID集合に現れる受信回数を計数して図4のようなデータを得る。また、全ての発信機IDの対ペアについて、そのペアが発信機ID集合に含まれる受信回数を計数して、図5のようなデータを得る。図1の受信情報生成部11は、単独受信回数と同時受信回数をまとめて受信情報として、受信情報DB12に送る。受信情報DB12は、受信情報を記憶する。発信機間距離算出部13は、発信機3の信号到達エリアを円と仮定し、単独受信回数と同時受信回数の比率からエリア重複部分の面積重複エリア面積を算出する。信号到達エリアの半径は、例えば、発信機情報DB14に記憶されている。すべての発信機3の信号到達エリアの半径が同一の場合は、半径情報を発信機間距離算出部13が保持していてもよい。そして、信号到達エリアの半径および重複エリア面積から、2つの発信機3の間の距離発信機間距離を算出する。図6は、実施の形態1における重複エリア面積の算出方法を示す。2つの発信機d1、d2の信号到達エリアをそれぞれ、半径r1、r2の円と仮定する。それぞれの信号到達エリアの面積をS1、S2、重複エリアの面積をS1_2とする図6a参照。d1の単独受信回数をn1、d2の単独受信回数をn2、d1とd2の同時受信回数をn12とする。S1:S1_2=n1:n12S2:S1_2=n2:n12とみなすと図6b、c参照、重複エリアの面積S1_2は、S1_2=S1n12/n1S1_2=S2n12/n2の2つの関係で表される。2つの重複エリアの面積を等しいと考えて平均をとれば、S1_2=S1n2S2n1n12/n1n2である。特に、r1=r2=rなら、S1=S2=πr2で、S1_2=πr2n12n1n2/n1n2である。これは、単独受信回数の平均と同時受信回数の比を、単独エリア面積と重複エリア面積の比とみなして計算したことになっている。発信機間距離算出部13は、エリア半径と重複エリア面積から発信機間距離を算出する。図7は、実施の形態1における発信機間の距離の算出方法を示す。2つの発信機d1、d2の信号到達範囲をそれぞれ、半径r1、r2の円とする。それぞれの中心発信機3の位置から、信号到達エリアの円の交点を見込む角度をそれぞれ2α、2βとする。重複エリアの面積S1_2はS1_2=αr12-r12sinαcosαβr22-r22sinβcosβただし、r1sinα=r2sinβである。この方程式からα、βを求めれば、発信機間の距離d12はd12=r1cosαr2cosβで計算できる。特に、r1=r2=rならα=βで、重複エリアの面積はS1_2=2αr2-2r2sin2αであり、距離は、d12=2rcosαである。図8は、実施の形態1において使用される発信機間距離情報の例を示す図である。発信機間距離算出部13は、例えば、図8のような発信機IDペアと発信機間距離からなる発信機間距離情報を生成する。ここで発信機間距離情報をDistanceinfo=<PairIDs,Distance>:PairIDs=発信機IDペア、Distance=発信機間距離と定義する。発信機間距離算出部13は、発信機IDペアと発信機間距離からなる発信機間距離情報を発信機情報DB14に記憶する。図1の発信機情報DB14は、発信機ID、設置座標からなる発信機情報を保持する。図9は、実施の形態1において使用される発信機情報の例を示す。発信機情報DB14は、初期状態では、配置している発信機3の内、三角形を構成することができる3点基点発信機の発信機情報を保持している。他の発信機情報に関しては発信機位置算出部15で生成された発信機情報を取得して保存する。ここで発信機情報SenderInfo=<ID,Coord>ID=発信機ID、Coord=設置座標と定義する。図1の発信機位置算出部15は、発信機情報と発信機間距離情報から、発信機3の位置座標を算出する。発信機3の位置座標は、発信機3が配置される領域に固定された座標系で表す。図10は、実施の形態1における座標候補算出方法を示す。発信機位置算出部15は、発信機情報DB14から発信機情報および発信機間距離情報を取得する。そして、座標が既知の2個の発信機A,Bについて、発信機A,Bまでの距離が発信機間距離情報に存在し、自身と発信機A,Bの3点間の距離が三角形成立条件を満たす発信機Cを選択する。既知の発信機A,Bと選択した発信機Cの3点において、発信機A,Bの座標からX軸に対する辺ABの角度θを算出し、さらに余弦定理を用いて3辺の長さから∠BACαを算出する。これらの角度からX軸に対する辺ACの角度を算出し、発信機Aの座標とX軸に対する辺ACの角度から発信機Cの候補点の座標を算出する。発信機Cの候補点C1,C2は2つ存在する図10。ここで、θ=arctanxB-xA/yB-yAcosα=d12d32-d22/2d1d2α=arccosd12d32-d22/2d1d2c1:xc1=xAd3cosθα,yc1=yAd3sinθαc2:xc2=xAd3cosθ-α,yc2=yAd3sinθ-αである。発信機位置算出部15は、以下のようにして2つの候補のうちの1つを決定する。図11は、実施の形態1において座標候補から最終的な推定座標を決定する方法を示す図である。発信機位置算出部15は、発信機Cとの距離DistanceDCが発信機間距離情報に存在し、かつ発信機3の座標が発信機情報に存在する発信機Dに対して、座標情報を用いて候補点C1、C2との距離を算出する。そして、算出した距離とDistanceDCを比較し、DistanceDCに値が近い方を最終的な発信機Cの位置とする。発信機位置算出部15は、発信機Cについて発信機IDと設置座標からなる発信機情報を生成する。以上の作業を発信機情報が生成できなくなるまで繰返し行う。計算式で表せば、以下のように書ける。DC=DistanceDCDC1=√xD-xc12yD-yc12DC2=√xD-xc22yD-yc22|DC1-DC|<|DC2-DC|ならば、C=C1とする。|DC1-DC|<|DC2-DC|でなければ、C=C2とする。発信機Cの最終的な位置の決定方法としては、他に2通りのやり方が考えられる。一つ目は、発信機Cとの距離が発信機間距離情報に存在し、かつ発信機3の座標が発信機情報に存在する全ての発信機3D={d1,d2,...,dm}に対して、座標情報を用いて各発信機3dk∈DとC1、C2との距離Distancedkc1、Distancedkc2を算出する。Distancedkc1、Distancedkc2とDistancedkcとの差Gapdkc1、Gapdkc2を算出する。そして、その積算値を比較し、値が小さい方を最終的な発信機Cの位置として採用する。二つ目は、一つ目の方法において座標を算出する対象の発信機Cとの距離が発信機間距離算出部13で算出した発信機機間距離に含まれず、かつ、発信機Cの座標の2つの候補C1、C2の一方について座標から算出した距離Distancedkcii=1または2が発信機Cと発信機dkのエリア半径の和よりも小さく、2つの候補の他方について座標から算出した距離Distancedkcjj=2または1が発信機Cと発信機dkのエリア半径の和より大きい場合に、座標から算出した距離Distancedkcii=1または2がそれぞれのエリア半径の和よりも小さい候補に対して、座標から算出した距離Distancedkcii=1または2を距離の差の積算値ΣGapdkciに加える。そして、2つの積算値数1を比較し、値が小さい方を最終的な発信機Cの位置として採用する。発信機Cの座標を決定したら、発信機Aまたは発信機Bと発信機Cを、改めて座標が既知の2つの発信機A、Bにして、上述の処理を繰り返す。このようにして、次々に発信機3の位置座標を算出する。初めに座標が与えられる3つの発信機3のうちのいずれの2つの発信機A,Bについても、発信機A,Bまでの距離が発信機間距離情報に存在し、自身と発信機A,Bの3点間の距離が三角形成立条件を満たす発信機Cがない場合には、発信機Aとの距離が発信機間距離情報に存在し、発信機Aを含んで三角形成立条件を満たす2つの発信機3を選択して、発信機Aに対する仮の座標を算出する。この場合、発信機Aを中心とする回転と鏡像の関係の自由度が残る。そこで、次々に発信機3の仮の座標を算出していって、座標が与えられる発信機3の仮の座標を算出したところで、その仮の座標と与えられた座標を一致させることによって、仮の座標から実際の領域に固定された座標系の座標を算出する。発信機位置算出部15は、他の2つの発信機3と発信機間距離情報が決定されている発信機3すべてについて、以上のようにして発信機IDと設置座標からなる発信機情報を生成する。発信機位置算出部15は、発信機情報を発信機DBに記憶する。位置計測装置2は、すべての発信機3の座標を計算したのち、発信機3の座標を用いて、ID受信機4の位置を推定する。位置推定部16は、1つの地点で前記通信機から受信した発信機ID集合または、移動体が1つの地点で発信した信号を受信した通信機の識別情報を、受信情報生成部11から取得する。そして、発信機ID集合に含まれる発信機IDに対応する発信機3の座標を発信機情報DB14から読み出し、それらの重心点を算出する。例えば、発信機ID集合が{S1,S2,S3}でありそれぞれの座標がS1=x1,y1、S2=x2,y2、S3=x3,y3の場合、最終的な位置座標はx1x2x3/3,y1y2y3/3となる。位置計測装置2は、この重心点の座標をID受信機4の推定位置とする。図12は、実施の形態1に係る発信機座標算出の動作の一例を示す流れ図である。ID受信機4を複数の発信機3が分散して配置される領域内で移動させて、複数の地点で発信機IDを受信するステップS11。例えば、ID受信機4を身につけた人物に計測対象区域発信機3が配置される領域を満遍なく一定回数歩き回ってもらって、ID受信機4で発信機IDを受信する。受信情報生成部11は、ID受信機4が受信した発信機ID集合をID受信機4から取得するステップS12。受信情報生成部11は、発信機ID毎の単独受信回数と、2つの発信機IDの対ごとの同時受信回数を計数し、発信機ID集合、受信回数単独受信回数および同時受信回数からなる受信情報を生成し、受信情報DB12へ登録するステップS13。次に、発信機間距離算出部13は、受信情報DB12から受信情報、発信機情報DB14から発信機情報を取得する。発信機間距離算出部13は、発信機3の到達エリア半径を発信機情報DB14から取得するステップS14。発信機3の信号到達エリアを円と仮定し、単独受信回数と同時受信回数の比を単独エリア面積と重複エリア面積の比として利用して重複エリア面積を算出する。そして、エリア半径、重複エリア面積から発信機間距離を算出するステップS15。発信機間距離算出部13は、発信機IDペアと発信機間距離からなる発信機間距離情報を生成して、発信機情報DB14に記憶する。次に、発信機位置算出部15は、発信機情報DB14から3つの発信機3の座標発信機情報および発信機間距離を取得するステップS16。そして、座標が既知の2個の発信機A,Bについて、発信機A,Bまでの距離が発信機間距離情報に存在し、自身と発信機A,Bの3点間の距離が三角形成立条件を満たす発信機Cを選択する。発信機A,B,Cの3点において、前述のとおり三角形の余弦定理を用いて、発信機Cの座標を算出するステップS17。発信機位置算出部15は、他の2つの発信機3と発信機間距離情報が決定されている発信機3すべてについて、発信機IDと設置座標からなる発信機情報を生成して、発信機情報DB14に記憶する。図13は、実施の形態に係る受信機の位置を推定する動作の一例を示す流れ図である。受信情報生成部11は、例えば、ID受信機4が所定の時間に受信した発信機IDを、ID受信機4から受信するのを待機するステップS21、ステップS22;NO。ID受信機4から発信機IDを受信するとステップS22;YES、受信した発信機IDの発信機3の座標の重心を算出するステップS23。そして、算出した重心の座標をID受信機4の位置として出力するステップS24。以上説明したように、本実施の形態の位置計測システム1によれば、ID受信機4で受信する信号の強度を用いずに、基点となる3個の発信機3の位置情報から、全ての発信機3の位置情報座標を算出することを可能とする。そして、算出した発信機3の座標を用いて、ID受信機4の位置を推定することができる。実施の形態2図14は、本発明の実施の形態2に係る位置計測システムの構成例を示すブロック図である。実施の形態2では、発信機間距離から、発信機3の仮想座標を算出し、そののち、三角形を構成する3つの発信機3の座標実座標を用いて、発信機3の仮想座標からその実座標を算出する。実施の形態2の位置計測装置2は、実施の形態1の構成に加えて、発信機位置算出部15に仮想座標算出部17を備える。その他の構成は、実施の形態1と同様である。受信情報生成部11は、実施の形態1と同様に、発信機ID集合をID受信機4から取得し、単独受信回数および同時受信回数を計数して受信情報を生成する。発信機間距離算出部13は、受信情報と発信機3の到達範囲のエリア半径から、重複エリア面積を算出して、発信機間の距離を算出する。そして、発信機間距離情報を発信機情報DB14に記憶する。仮想座標算出部17は、発信機情報DB14から発信機間距離情報を取得する。発信機3のいずれか1つを原点とし、他のいずれか1つと原点を通る直線を座標軸の1つとする仮想座標系を設定する。そして、発信機間距離から仮想座標系における発信機3の仮想座標を算出する。発信機位置算出部15は、発信機情報DB14から三角形を構成することができる3点基点発信機の実座標を含む発信機情報を取得し、発信機3の仮想座標と発信機情報から発信機3の実座標を算出する。図15は、実施の形態2における基点発信機の仮想座標の算出方法を示す図である。仮想座標算出部17は、発信機間距離から三角形を構成することが可能な任意の3点A,B,Cを選択する。発信機Aの座標を0,0、発信機Bの座標をDistanceAB,0とし、余弦定理を用いて3辺の長さから∠BACαを算出する。そして、発信機Cの座標をDistanceACcosα,DistanceACsinαと算出する。発信機Cの仮想座標について候補が2つ算出されるが、そのうち、最終的な仮想座標は、実施の形態1と同様に、仮想座標における既知の発信機3との距離によって決定する。4点目以降の発信機3について、実施の形態1の発信機位置算出部15の座標算出方法と同様に、仮想座標系における座標を算出する。そして算出された仮想座標から発信機ID、仮想座標からなる仮想座標情報を生成する。ここで仮想座標情報VirtualCoordInfo=<ID,VirtualCoord>ID=発信機ID、Coord=仮想座標と定義する。発信機位置算出部15は、仮想座標算出部17から仮想座標情報を、発信機情報DB14から基点発信機I,J,Kの実座標を含む発信機情報を取得する。発信機I,J,Kの仮想座標と実座標を用いて仮想座標から実座標への変換式を生成する。図16は、実施の形態2における実座標の算出方法を示す図である。実座標系xr-yrと仮想座標系xv-yvとの関係を図16のように表す。すなわち、実座標系の原点に対する仮想座標系原点の位置変位をa,b、実座標系に対する仮想座標系の傾き角をβγ=β|π-βとする。仮想座標系から実座標系への変換式は、xr=axvcosβ-yvsinβ、yr=bxvsinβyvcosβaxIvcosβ-yIvsinγ=xIr、bxIvsinβyIvcosγ=yIr、axJvcosβ-yJvsinγ=xJr、bxJvsinβyJvcosγ=yJr、axKvcosβ-yKvsinγ=xKr、bxKvsinβyKvcosγ=yKr、γ=β|π-βを満たすa、b、β、γとなる。この変換式を用いて、仮想座標情報における仮想座標を実座標に変換した値を設置座標とし、発信機ID、設置座標からなる発信機情報を生成する。図17は、実施の形態2に係る発信機座標算出の動作の一例を示す流れ図である。図17のステップS31からステップS35は、図13のステップS11からステップS15と同様である。発信機間距離を算出したステップS35のち、仮想座標算出部17は、前述のように発信機間距離から発信機3の仮想座標を算出するステップS36。すべての発信機3について、仮想座標を算出すると、発信機位置算出部15は、3つの発信機3の実座標を含む発信機情報を発信機情報DB14から取得するステップS37。そして仮想座標と発信機情報から、前述のように仮想座標から実座標への変換式を生成し、発信機3の座標実座標を算出するステップS38。発信機3の座標実座標を算出したのち、実施の形態1と同様にして、受信機の位置を推定することができる。以上説明したように、本実施の形態2の位置計測システム1によれば、最初に座標が与えられる3つの発信機3の相互の発信機間距離が受信情報から直接的に算出されなくても、いずれの発信機3も他の2つ以上の発信機3との発信機間距離が計算できれば、発信機3の座標を決定することができる。次に、具体的な例を用いて本発明の実施の形態の作用を説明する。図18は、実施の形態の具体例における発信機3と観測点の配置を示す図である。図18の例では、13m×13mの計測対象エリアに対して、信号到達エリア半径が1.5mの発信機100個を、1m間隔の格子状の100点10×10の基準点それぞれから、半径25cmの範囲内に1個ずつランダムに配置した。観測点受信機が発信機IDを受信した場所を計測対象エリアにランダムに5000点用意した場合を考える。計測エリアを1m2の格子に区切って全ての格子を通るようにID受信機4を身につけた人物が歩き回って貰い、各格子で1回発信機IDを受信できるとして、約30回計測エリアを歩き回ってもらった。このような場合において、受信情報生成部11は、ID受信機4から図2のような発信機ID集合を取得し、図4および図5のような受信情報を生成した。発信機間距離情報算出部103は、発信機3の信号到達エリアを円と仮定し全発信機3のエリア半径は同一であり、半径情報は保持している、単独受信回数の平均と同時受信回数の比を単独エリア面積と重複エリア面積の比として利用して重複エリア面積DupAreaを算出した。例えば、エリア半径が1m、単独エリア面積をπm2とし、図3においてS1とS2の発信機間距離情報を生成する場合、それぞれの単独受信回数がSingleRevNumS1=181084=40、SingleRevNumS2=10414763=44、同時受信回数がDoubleRevNumS1S2=104=14となり、重複エリア面積はDupAreaS1S2=π14/4044/2となる。発信機間距離算出部13は、エリア半径、重複エリア面積から発信機間距離を算出し図7、図8のような発信機IDペアと発信機間距離からなる発信機間距離情報を生成した。発信機間距離は、π/3=α-sin2αDistanceS1S2=2cosαの式を解きDistanceS1S2=1.438となる。発信機情報DB14には、三角形を構成できるような任意の3個の発信機3について図9のような発信機IDと設置座標からなる発信機情報を登録しておく。発信機位置算出部15は、発信機情報DB14から発信機情報および発信機間距離情報を取得し、前述の方法で座標が既知の2個の発信機A,Bから出発して、三角形を構成できる発信機間距離の発信機3の座標を算出する。発信機A,Bまでの距離が発信機間距離情報に存在し、自身と発信機A,Bの3点間の距離が三角形成立条件を満たす発信機Cの3点において、発信機A,Bの座標からX軸に対する辺ABの角度を算出し、さらに余弦定理を用いて3辺の長さから∠BACαを算出し、これらの角度からX軸に対する辺ACの角度θを算出し、発信機Aの座標とX軸に対する辺ACの角度から発信機Cの候補点C1,C2の座標を算出する図10。次に、発信機Cとの距離DistanceDCが発信機間距離情報に存在し、かつ発信機3の座標が発信機情報に存在する発信機Dに対して、座標情報を用いて候補点C1、C2との距離を算出し、算出した距離とDistanceDCを比較し、DistanceDCに値が近い方を最終的な発信機Cの位置とし図11、発信機Cについて発信機ID、設置座標からなる発信機情報を生成する。以上の作業を発信機情報が生成できなくなるまで繰返し行う。例えば、図8、図9においてS1、S2の座標、S1S2、S1S4、S2S4の距離からS4の座標を算出する場合、∠S2S1S4はα=arccos1.4520.672-1.122/2×1.45×0.67=48°となる。X軸に対するS1S2の角度はθ=90°となる。S4の候補点S41、S42の座標はそれぞれS41=1.021.45cos48°90°、1.361.45sin48°90°=2.11、2.33、S42=1.021.45cos48°-90°、1.361.45sin48°-90°=-0.1、2.33となる。S42、S41とS3の座標からS3S41,S3S42とS3S4の距離は、S3S41=√2.11-2.0122.33-1.502=0.83、S3S42=√-0.1-2.0122.33-1.502=2.27となる。図9のS3S4=0.82と比較し、S3S41の方が値が近いので最終的にS4の座標を2.11、2.33とする。以上のような計算を繰り返して、図18に示すような発信機3の推定位置を得た。図18では、発信機3の設置位置を○印で、発信機3の推測位置を□印で示す。観測点発信機IDの受信位置を小さい点で示す。観測点の間隔を均等にし、さらに密度を大きく観測点を多くすることによって、発信機3の推定位置の精度を向上できる。位置推定部16は、受信情報生成部11から1つの地点で受信した発信機ID集合を取得し、発信機ID集合の全てのIDに対する座標の重心を最終的な位置として出力する。例えば、発信機ID集合が{S1,S2,S3}であり、それぞれの座標が図9の値の場合、最終的な位置座標は1.021.022.01/3,1.362.081.50/3=1.34,1.31となる。本発明の実施の形態2の仮想座標の算出について、具体的な例を説明する。例えば、図8におけるS1S2、S1S4、S2S4の距離からS1、S2、S4の仮想座標を算出する場合、∠S2S1S4はα=arccos1.4520.672-1.122/2×1.45×0.67=48°となる。S1を原点にして、S1とS2を通る直線をx軸とする仮想座標系を設定する。各点の仮想座標は、S1=0,0、S2=0.72,0、S4=1.45cos48°,1.45sin48°である。S1、S2、S4の3点を元に、4点目以降の仮想座標を算出する。すべての発信機3の仮想座標を算出したのち、発信機位置算出部15は、仮想座標算出部17から仮想座標情報を、発信機情報DB14から基点発信機I,J,Kの発信機情報を取得し、発信機I,J,Kの仮想座標と実座標を用いて仮想座標から実座標への変換式を生成する。実施の形態2で説明し仮想座標と実座標の変換式で、例えば、Ir=3,2、Jr=6,4、Kr=4,5、Iv=1,0.5、Jv=4,1、Kv=2.5,3の場合、a1cosβ-0.5sinγ=3、b1sinβ0.5cosγ=2、a4cosβ-1sinγ=6、b4sinβ1cosγ=4、a2.5cosβ-3sinγ=4、b2.5sinβ3cosγ=5の連立方程式を解き、a,b,β,γを算出して、変換式xr=axvcosβ-yvsinβ、yr=bxvsinβyvcosβを生成する。図19は、本発明の実施の形態に係る位置計測装置の物理的な構成例を示すブロック図である。図19は、図1または図14の位置計測装置2の構成を示す。位置計測装置2は、図19に示すように、制御部21、主記憶部22、外部記憶部23、操作部24、表示部25、入出力部26および送受信部27を備える。主記憶部22、外部記憶部23、操作部24、表示部25、入出力部26および送受信部27はいずれも内部バス20を介して制御部21に接続されている。制御部21はCPUCentralProcessingUnit等から構成され、外部記憶部23に記憶されている制御プログラム29に従って、位置計測のための処理を実行する。主記憶部22はRAMRandomAccessMemory等から構成され、外部記憶部23に記憶されている制御プログラム29をロードし、制御部21の作業領域として用いられる。外部記憶部23は、フラッシュメモリ、ハードディスク、DVD-RAMDigitalVersatileDiscRandomAccessMemory、DVD-RWDigitalVersatileDiscReWritable等の不揮発性メモリから構成され、上述の処理を制御部21に行わせるための制御プログラム29を予め記憶し、また、制御部21の指示に従って、この制御プログラム29が記憶するデータを制御部21に供給し、制御部21から供給されたデータを記憶する。操作部24はキーボードおよびマウスまたはタッチパネルなどのポインティングデバイス等と、キーボードおよびポインティングデバイス等を内部バス20に接続するインタフェース装置から構成されている。操作部24を介して、基点となる3つの発信機3の座標を含む発信機情報および発信機3のエリア半径などの入力を受付ける。表示部25は、LCDLiquidCrystalDisplayもしくは有機ELディスプレイ、およびスピーカなどから構成され、発信機3の座標またはID受信機4の推定位置を表示する。入出力部26は、シリアルインタフェースまたはパラレルインタフェースから構成されている。入出力部26に発信機3が接続される。ID受信機が発信する信号を発信機3が受信する場合、発信機3が受信した情報を入出力部26を介して入力し、制御部21に送る。送受信部27は、無線送受信機、無線モデムまたは網終端装置、およびそれらと接続するシリアルインタフェースまたはLANLocalAreaNetworkインタフェースから構成されている。位置計測装置2は、送受信部27を介して、ID受信機4から発信機ID集合を受信する。位置計測装置2の受信情報生成部11、受信情報DB12、発信機間距離算出部13、発信機情報DB14、発信機位置算出部15、位置推定部16などの処理は、制御プログラム29が、制御部21、主記憶部22、外部記憶部23、操作部24、表示部25、入出力部26および送受信部27などを資源として用いて処理することによって実行する。その他、前記のハードウエア構成やフローチャートは一例であり、任意に変更および修正が可能である。制御部21、主記憶部22、外部記憶部23、操作部24、内部バス20などから構成される制御処理を行う中心となる部分は、専用のシステムによらず、通常のコンピュータシステムを用いて実現可能である。たとえば、前記の動作を実行するためのコンピュータプログラムを、コンピュータが読み取り可能な記録媒体フレキシブルディスク、CD-ROM、DVD-ROM等に格納して配布し、当該コンピュータプログラムをコンピュータにインストールすることにより、前記の処理を実行する位置計測装置2を構成してもよい。また、インターネット等の通信ネットワーク上のサーバ装置が有する記憶装置に当該コンピュータプログラムを格納しておき、通常のコンピュータシステムがダウンロード等することで位置計測装置2を構成してもよい。また、位置計測装置2の機能を、OSオペレーティングシステムとアプリケーションプログラムの分担、またはOSとアプリケーションプログラムとの協働により実現する場合などには、アプリケーションプログラム部分のみを記録媒体や記憶装置に格納してもよい。また、搬送波にコンピュータプログラムを重畳し、通信ネットワークを介して配信することも可能である。たとえば、通信ネットワーク上の掲示板に前記コンピュータプログラムを掲示し、ネットワークを介して前記コンピュータプログラムを配信してもよい。そして、このコンピュータプログラムを起動し、OSの制御下で、他のアプリケーションプログラムと同様に実行することにより、前記の処理を実行できるように構成してもよい。上記の実施形態の一部または全部は、以下の付記のようにも記載されうるが、以下には限られない。付記1それぞれが所定の周期で識別情報を無線で発信する、または、所定の周期で無線で信号を発信する移動体から該信号を受信する、分散して配置される複数の通信機と、前記通信機が発信する識別情報を受信する移動体、または、前記所定の周期で無線で信号を発信する移動体と、所定の期間、所定の領域内の複数の地点で前記移動体が前記通信機から受信した前記識別情報および異なる2つの前記識別情報を受信した地点が同一か否かを表す情報の集合、あるいは、前記移動体が所定の期間、所定の領域内を移動しながら前記所定の周期で発信する信号を受信した前記通信機の識別情報および異なる前記通信機が受信した信号が同一の地点から発信されたか否かを表す情報の集合、である通信機ID集合を取得する収集手段と、前記通信機ID集合に含まれる前記識別情報ごとの受信回数、および、前記2つの識別情報の組ごとの前記識別情報の受信もしくは発信した地点が同一である受信回数、を計数した受信情報を生成する受信情報生成手段と、前記通信機ごとの発信の到達範囲もしくは受信可能な前記移動体の発信範囲を示す通信範囲情報を取得する範囲情報取得手段と、前記受信情報および前記通信範囲情報から、前記複数の通信機に含まれる2つの前記通信機の間の通信機間距離を算出する距離算出手段と、前記複数の通信機のうち、その位置を頂点とする三角形を構成できる3つの通信機それぞれの前記識別情報および設置座標を含む通信機情報を取得する手段と、前記通信機間距離と前記通信機情報から、前記通信機それぞれの設置座標を算出する座標算出手段と、を備えることを特徴とする位置計測システム。付記2前記通信範囲情報は、前記複数の通信機それぞれの前記到達範囲もしくは前記発信範囲を示すエリア半径を含み、前記距離算出手段は、前記通信機の前記到達範囲もしくは前記発信範囲を、前記通信機を中心とする前記エリア半径の円とみなして、前記通信機間距離を算出することを特徴とする付記1に記載の位置計測システム。付記3前記距離算出手段は、前記2つの通信機それぞれの前記エリア半径、および、前記到達範囲もしくは前記発信範囲の重複面積から前記通信間距離を算出することを特徴とする付記2に記載の位置計測システム。付記4前記距離算出手段は、前記2つの通信機それぞれの前記到達範囲もしくは前記発信範囲の面積に、前記2つの通信機の識別情報ごとの受信回数に対する前記受信もしくは発信した地点が同一である受信回数の比率をかけることによって、前記重複面積を算出することを特徴とする付記3に記載の位置計測システム。付記5前記複数の通信機それぞれの前記エリア半径は同一であることを特徴とする付記2ないし4のいずれかに記載の位置計測システム。付記6前記座標算出手段は、前記通信機間距離から余弦定理を用いて、2つの前記通信機に対する他の1つの前記通信機の相対座標を計算することを特徴とする付記1ないし5のいずれかに記載の位置計測システム。付記7前記座標算出手段は、前記相対座標の2つの候補の内、前記相対座標を算出する対象の前記通信機との距離が前記距離算出手段で算出した前記通信機間距離に含まれ、かつ、座標が前記座標算出手段で算出された前記通信機との間で、前記距離算出手段で算出した距離と、前記座標および前記相対座標の候補から算出した距離と、の差が小さい前記候補を前記相対座標として算出することを特徴とする付記6に記載の位置計測システム。付記8前記座標算出手段は、前記相対座標の2つの候補の内、前記相対座標を算出する対象の前記通信機との距離が前記距離算出手段で算出した前記通信機間距離に含まれ、かつ、座標が前記座標算出手段で算出された前記通信機全てに対して、前記距離算出手段で算出した距離と、前記座標および前記相対座標の候補から算出した距離との差を、前記相対座標の候補それぞれについて積算した値を比べ、値が小さい前記候補を前記相対座標として算出することを特徴とする付記7に記載の位置計測システム。付記9前記座標算出手段は、前記相対座標を算出する対象の前記通信機との距離が前記距離算出手段で算出した通信機間距離に含まれず、かつ、前記相対座標の2つの候補の一方について座標から算出した距離がそれぞれの前記エリア半径の和よりも小さく、前記2つの候補の他方について座標から算出した距離がそれぞれの前記エリア半径の和より大きい場合に、前記座標から算出した距離がそれぞれの前記エリア半径の和よりも小さい候補に対して、前記座標から算出した距離を前記積算した値に加える、ことを特徴とする付記8に記載の位置計測システム。付記10前記座標算出手段は、前記通信機のいずれか1つを原点とし、他のいずれか1つと前記原点を通る直線を座標軸の1つとする仮想座標系を設定して、前記通信機間距離から前記通信機の前記仮想座標系における仮想座標を算出し、前記通信機情報と前記仮想座標から前記通信機の座標を算出することを特徴とする付記6ないし9のいずれかに記載の位置計測システム。付記11前記移動体が1つの地点で前記通信機から受信した識別情報、または、前記移動体が1つの地点で発信した信号を受信した前記通信機の識別情報、を取得し、該識別情報に対応する前記座標算出手段で算出した前記通信機の座標の重心点を算出する重心算出手段を備えることを特徴とする付記1ないし10のいずれかに記載の位置計測システム。付記12それぞれが所定の周期で識別情報を無線で発信する分散して配置される複数の通信機から、所定の期間、所定の領域内の複数の地点で受信した識別情報および異なる2つの識別情報を受信した地点が同一か否かを表す情報の集合、あるいは、所定の周期で無線で信号を発信する移動体が所定の期間、所定の領域内を移動しながら発信する信号を受信した、分散して配置される複数の通信機の識別情報および異なる前記通信機が受信した信号が同一の地点から発信されたか否かを表す情報の集合、である通信機ID集合を取得する収集手段と、前記通信機ID集合に含まれる前記識別情報ごとの受信回数、および、前記2つの識別情報の組ごとの前記識別情報の受信もしくは発信した地点が同一である受信回数、を計数した受信情報を生成する受信情報生成手段と、前記通信機ごとの発信の到達範囲もしくは受信可能な前記移動体の発信範囲を示す通信範囲情報を取得する範囲情報取得手段と、前記受信情報および前記通信範囲情報から、前記複数の通信機に含まれる2つの前記通信機の間の通信機間距離を算出する距離算出手段と、前記複数の通信機のうち、その位置を頂点とする三角形を構成できる3つの通信機それぞれの前記識別情報および設置座標を含む通信機情報を取得する手段と、前記通信機間距離と前記通信機情報から、前記通信機それぞれの設置座標を算出する座標算出手段と、を備えることを特徴とする位置計測装置。付記13それぞれが所定の周期で識別情報を無線で発信する分散して配置される複数の通信機から、所定の期間、所定の領域内の複数の地点で受信した識別情報および異なる2つの識別情報を受信した地点が同一か否かを表す情報の集合、あるいは、所定の周期で無線で信号を発信する移動体が所定の期間、所定の領域内を移動しながら発信する信号を受信した、分散して配置される複数の通信機の識別情報および異なる前記通信機が受信した信号が同一の地点から発信されたか否かを表す情報の集合、である通信機ID集合を取得する収集ステップと、前記通信機ID集合に含まれる前記識別情報ごとの受信回数、および、前記2つの識別情報の組ごとの前記識別情報の受信もしくは発信した地点が同一である受信回数、を計数した受信情報を生成する受信情報生成ステップと、前記通信機ごとの発信の到達範囲もしくは受信可能な前記移動体の発信範囲を示す通信範囲情報を取得する範囲情報取得ステップと、前記受信情報および前記通信範囲情報から、前記複数の通信機に含まれる2つの前記通信機の間の通信機間距離を算出する距離算出ステップと、前記複数の通信機のうち、その位置を頂点とする三角形を構成できる3つの通信機それぞれの識別情報および設置座標を含む通信機情報を取得するステップと、前記通信機間距離と前記通信機情報から、前記通信機それぞれの設置座標を算出する座標算出ステップと、を備えることを特徴とする位置計測方法。付記14前記通信範囲情報は、前記複数の通信機それぞれの前記到達範囲もしくは前記発信範囲を示すエリア半径を含み、前記距離算出ステップは、前記通信機の前記到達範囲もしくは前記発信範囲を、前記通信機を中心とする前記エリア半径の円とみなして、前記通信機間距離を算出することを特徴とする付記13に記載の位置計測方法。付記15前記距離算出ステップは、前記2つの通信機それぞれの前記エリア半径、および、前記到達範囲もしくは前記発信範囲の重複面積から前記通信間距離を算出することを特徴とする付記14に記載の位置計測方法。付記16前記距離算出ステップは、前記2つの通信機それぞれの前記到達範囲もしくは前記発信範囲の面積に、前記2つの通信機の識別情報ごとの受信回数に対する前記受信もしくは発信した地点が同一である受信回数の比率をかけることによって、前記重複面積を算出することを特徴とする付記15に記載の位置計測方法。付記17前記複数の通信機それぞれの前記エリア半径は同一であることを特徴とする付記14ないし16のいずれかに記載の位置計測方法。付記18前記座標算出ステップは、前記通信機間距離から余弦定理を用いて、2つの前記通信機に対する他の1つの前記通信機の相対座標を計算することを特徴とする付記13ないし17のいずれかに記載の位置計測方法。付記19前記座標算出ステップは、前記相対座標の2つの候補の内、前記相対座標を算出する対象の前記通信機との距離が前記距離算出ステップで算出した前記通信機間距離に含まれ、かつ、座標が前記座標算出ステップで算出された前記通信機との間で、前記距離算出ステップで算出した距離と、前記座標および前記相対座標の候補から算出した距離と、の差が小さい前記候補を前記相対座標として算出することを特徴とする付記18に記載の位置計測方法。付記20前記座標算出ステップは、前記相対座標の2つの候補の内、前記相対座標を算出する対象の前記通信機との距離が前記距離算出ステップで算出した前記通信機間距離に含まれ、かつ、座標が前記座標算出ステップで算出された前記通信機全てに対して、前記距離算出ステップで算出した距離と、前記座標および前記相対座標の候補から算出した距離との差を、前記相対座標の候補それぞれについて積算した値を比べ、値が小さい前記候補を前記相対座標として算出することを特徴とする付記19に記載の位置計測方法。付記21前記座標算出ステップは、前記相対座標を算出する対象の前記通信機との距離が前記距離算出ステップで算出した通信機間距離に含まれず、かつ、前記相対座標の2つの候補の一方について座標から算出した距離がそれぞれの前記エリア半径の和よりも小さく、前記2つの候補の他方について座標から算出した距離がそれぞれの前記エリア半径の和より大きい場合に、前記座標から算出した距離がそれぞれの前記エリア半径の和よりも小さい候補に対して、前記座標から算出した距離を前記積算した値に加える、ことを特徴とする付記20に記載の位置計測方法。付記22前記座標算出ステップは、前記通信機のいずれか1つを原点とし、他のいずれか1つと前記原点を通る直線を座標軸の1つとする仮想座標系を設定して、前記通信機間距離から前記通信機の前記仮想座標系における仮想座標を算出し、前記通信機情報と前記仮想座標から前記通信機の座標を算出することを特徴とする付記18ないし21のいずれかに記載の位置計測方法。付記23前記移動体が1つの地点で前記通信機から受信した識別情報、または、前記移動体が1つの地点で発信した信号を受信した前記通信機の識別情報、を取得し、該識別情報に対応する前記座標算出ステップで算出した前記通信機の座標の重心点を算出する重心算出ステップを備えることを特徴とする付記13ないし22のいずれかに記載の位置計測方法。付記24コンピュータに、それぞれが所定の周期で識別情報を無線で発信する分散して配置される複数の通信機から、所定の期間、所定の領域内の複数の地点で受信した識別情報および異なる2つの識別情報を受信した地点が同一か否かを表す情報の集合、あるいは、所定の周期で無線で信号を発信する移動体が所定の期間、所定の領域内を移動しながら発信する信号を受信した、分散して配置される複数の通信機の識別情報および異なる前記通信機が受信した信号が同一の地点から発信されたか否かを表す情報の集合、である通信機ID集合を取得する収集ステップと、前記通信機ID集合に含まれる前記識別情報ごとの受信回数、および、前記2つの識別情報の組ごとの前記識別情報の受信もしくは発信した地点が同一である受信回数、を計数した受信情報を生成する受信情報生成ステップと、前記通信機ごとの発信の到達範囲もしくは受信可能な前記移動体の発信範囲を示す通信範囲情報を取得する範囲情報取得ステップと、前記受信情報および前記通信範囲情報から、前記複数の通信機に含まれる2つの前記通信機の間の通信機間距離を算出する距離算出ステップと、前記複数の通信機のうち、その位置を頂点とする三角形を構成できる3つの通信機それぞれの前記識別情報および設置座標を含む通信機情報を取得するステップと、前記通信機間距離と前記通信機情報から、前記通信機それぞれの設置座標を算出する座標算出ステップと、を実行させることを特徴とするプログラム。"