"雷の観測地点からみた雷から放射される信号の到来方向を示す仰角および方位角、当該信号を受信した時刻、前記観測地点の位置情報、および送信元を識別する識別情報を送信する複数の雷放射観測装置とネットワークを介して通信可能に接続される雷発生位置標定装置であって、前記雷放射観測装置から、前記仰角の時系列、前記方位角、前記時刻、前記識別情報、および前記位置情報を含む観測データを受信する通信部と、前記雷放射観測装置から受信した前記識別情報の異なる前記仰角の時系列の間で、当該仰角の時系列の時刻の間隔および仰角の値が類似するか否かを判定する類似判定部と、前記類似判定部によって類似すると判定した複数の前記識別情報の異なる前記観測データを用いて、前記雷の発生位置を標定する雷位置標定部とを備えることを特徴とする雷発生位置標定装置。前記雷位置標定部は、前記類似判定部によって類似すると判定された複数の前記識別情報の異なる前記観測データにおいて、類似すると判定された前記仰角と前記方位角とで構成される方向ベクトルの中から2組の方向ベクトルを選択しては、それらの距離が最小になるときの中点を求め、算出した前記中点の位置の平均値を算出し、当該平均値の位置を前記雷の発生位置として標定することを特徴とする請求項1に記載の雷発生位置標定装置。前記雷発生位置標定装置は、雲の動きを少なくとも含む気象情報データに基づいて、前記雷位置標定部によって標定した前記雷の発生位置の移動を予測し、雷発生地域を推定する雷予測部を備えることを特徴とする請求項1または請求項2に記載の雷発生位置標定装置。前記雷放射観測装置は、所定の半径を有する円状の領域の中心に設置され、前記円状の領域が、前記雷の発生位置の観測対象のエリアを含むように配置していることを特徴とする請求項1ないし請求項3のいずれか一項に記載の雷発生位置標定装置。雷の観測地点からみた雷から放射される信号の到来方向を示す仰角および方位角、当該信号を受信した時刻、GPSGlobalPositioningSystem受信機から取得した前記観測地点の位置情報、および送信元を識別する識別情報を雷発生位置標定装置に送信し、移動体に搭載された雷放射観測装置と、ネットワークを介して通信可能に接続される前記雷発生位置標定装置であって、前記雷放射観測装置から、前記仰角の時系列、前記方位角、前記時刻、前記識別情報、および前記位置情報を含む観測データを受信する通信部と、前記雷放射観測装置から受信した前記識別情報の異なる前記仰角の時系列の間で、当該仰角の時系列の時刻の間隔および仰角の値が類似するか否かを判定する類似判定部と、前記類似判定部によって類似すると判定した複数の前記識別情報の異なる前記観測データを用いて、前記雷の発生位置を標定する雷位置標定部とを備えることを特徴とする雷発生位置標定装置。雷の観測地点からみた雷から放射される信号の到来方向を示す仰角および方位角、当該信号を受信した時刻、前記観測地点の位置情報、および送信元を識別する識別情報を送信する複数の雷放射観測装置とネットワークを介して通信可能に接続される雷発生位置標定装置で用いられる雷発生地域推定方法あって、前記雷発生位置標定装置は、前記雷放射観測装置から、前記仰角の時系列、前記方位角、前記時刻、前記識別情報、および前記位置情報を含む観測データを受信するステップと、前記雷放射観測装置から受信した前記識別情報の異なる前記仰角の時系列の間で、当該仰角の時系列の時刻の間隔および仰角の値が類似するか否かを判定する類似判定ステップと、前記類似判定ステップによって類似すると判定した複数の前記識別情報の異なる前記観測データを用いて、前記雷の発生位置を標定する雷位置標定ステップと雲の動きを少なくとも含む気象情報データに基づいて、前記雷位置標定ステップによって標定した雷の発生位置の移動を予測し、雷発生地域を推定する雷予測ステップとを実行することを特徴とする雷発生地域推定方法。しかしながら、いずれの先行技術も雷発生位置とその観測位置の距離が離れるとその観測が困難になる。落雷を光学的に観測する場合、通常は雷雨の存在により観測位置までの距離がのびるにつれその閃光が見えにくくなり正確に雷発生位置への方向ベクトルが得にくくなる。また、電磁波は大気中での減衰が周波数の高周波化にともない著しくなる。落雷リーダの強度にもよるが、そのVHF波帯域の電磁波の到達距離は最大で10km程度と考えられている。従って、画像や映像観測、VHF波帯域電磁波計測いずれの場合においても、数十km以上の広範囲にわたるエリアの雷発生位置を観測するためには何らかの手段が必要となる。そこで、本発明は、雷の活動状況、特にその発生位置を広範囲のエリアにわたって、精度良く測定または推定する技術を提供することを課題とする。前記課題を解決するために、雷の閃光、電磁波あるいはその他の雷の発生に伴う信号例えば雷鳴等も含むの発生位置を推定する本発明の雷発生位置測定システムは、複数の雷放射観測装置上記雷の発生に伴う信号の物理量を受信し、その雷の発生位置への方向ベクトルを出力する計測装置と少なくとも1つの雷発生位置標定装置とで構成する。雷放射観測装置は、雷の観測対象のエリアに複数配置され、雷から放射される信号を受信素子を介して受信して、受信地点からみた信号の到来方向を示す方位角および仰角、当信号を受信した時刻、受信地点の位置情報、および送信元を識別する識別情報を前記雷発生位置標定装置に送信する。また、前記雷発生位置標定装置は、前記雷放射観測装置から、前記仰角の時系列、前記方位角、前記時刻、前記識別情報、および前記位置情報を含む観測データを受信し、前記雷放射観測装置から受信した前記識別情報の異なる前記仰角の時系列の間で、当該仰角の時系列の時刻の間隔および仰角の値が類似するか否かを判定し、類似すると判定した複数の前記識別情報の異なる前記観測データ受信時刻、方位角および仰角の時系列を用いて、雷の発生位置を標定する。本発明によれば、雷の広範囲な観測対象エリアにおいて、雷の発生位置を精度良く標定または推定する技術を提供することができる。次に、本発明を実施するための形態以降、本実施形態と称す。について、適宜図面を参照しながら詳細に説明する。雷発生位置測定システムの構成本実施形態における雷発生位置測定システムの構成例について、図1を用いて説明する。図1に示すように、本実施形態における雷発生位置測定システム1は、雷発生位置標定装置10と、雷発生位置標定装置10とネットワーク11を介して通信可能に接続される雷放射観測装置110110a,110b,110cとで構成される。雷放射観測装置110は、雷の観測対象のエリアに複数配置される。雷放射観測装置110は、360度方位カメラ等を用いた雷の閃光の画像や映像観測手段、無指向性マイク等を用いた雷鳴の検出手段、またはアンテナ等を用いた雷から放出される電磁波の受信手段等を備えることによって、観測地点からみた雷の発生位置への方向ベクトル方位角αおよび仰角βを求めることができる。そして、雷放射観測装置110は、方位角αおよび仰角βの時系列、その時系列の時刻、および観測地点の位置情報を関連付けて記録する。なお、雷放射観測装置110は、予めその設置位置が分かっている場合には位置情報を記憶しておく。または、雷放射観測装置110は、GPSGlobalPositioningSystem受信機不図示を備えて、位置情報を取得しても良い。また、雷放射観測装置110は、他の雷放射観測装置110と同期した時計を備えて時刻を取得する、または、GPS受信機を備えて同期した時刻を取得する。ネットワーク11は、通信方式を限定されることなく、有線であっても無線であっても構わない。雷発生位置標定装置10は、複数の雷放射観測装置110からネットワーク11を介して受信した雷発生位置の方向に関する情報を用いて、雷発生位置を標定する。なお、雷発生位置標定装置10の詳細については後記する。また、図1では、雷放射観測装置110は、3台記載しているが、3台に限られることなく、2台または4台以上であっても構わない。図1において、破線で示す円状の領域13は、その中心に設定される雷放射観測装置110を雷の観測対象のエリアに配置するときに、その円状の領域13が当該エリアをカバーするために用いられる。この円状の領域13の詳細については、後記する。雷発生位置標定装置次に、雷発生位置標定装置10の機能例について、図2を用いて説明する。図2に示すように、雷発生位置標定装置10は、処理部21、通信部22、および記憶部23を備える。通信部22は、雷放射観測装置110から雷発生位置の方向に関する情報を受信し、受信した情報を処理部21へ引き渡す。なお、受信した情報は、少なくとも、図3の受信情報301に示すように、送信元の雷放射観測装置110を識別する識別情報、その雷放射観測装置110の設置されている位置情報例えば、緯度経度等、仰角βおよび方位角αの時系列、および当該時系列の時刻である。また、処理部21は、受信情報301を記憶部23へ記憶する。なお、例えば、ネットワーク11がIPInternetProtocol網の場合には、図3に示す識別情報は送信元のIPアドレスであっても良い。また、送信元のIPアドレスから雷放射観測装置110の位置が取得できる場合には、図3に示す位置情報は不要である。逆に、位置情報が識別情報の代わりになりうるものであれば、識別情報は不要である。図2へ戻り、処理部21は、図示しないCPUCentralProcessingUnitおよびメインメモリによって構成され、記憶部23に記憶されているアプリケーションプログラムをメインメモリに展開して、後記する類似判定部24、雷位置標定部25、および雷予測部26を具現化する。類似判定部24は、識別情報の異なる仰角βの時系列の間で、当該仰角βの時系列の時刻の間隔および仰角βの値が類似するか否か、すなわち類似性を判定する。例えば、図4は、方向ベクトルの仰角βを時系列に則って模式的にプロット□印で表示したものである。ここで、図4aは、仰角βの時系列が類似すると判定する場合を示している。すなわち、図4a1に示す雷放射観測装置110aでの仰角βの時系列は、時刻を所定の時間範囲内で若干ずらすと、図4a2に示す雷放射観測装置110cでの仰角βの時系列とほぼ一致し、仰角βの値の時系列の傾向例えば、仰角βの値を時系列に沿って結んだ線の形状も類似していることが分かる。なお、類似性の判定には、種々の統計処理を用いることが可能であり、特に限定されるものではないが、例えば、2つの仰角βの時系列について相関係数を求め、所定の値以上であるときに類似していると判定しても良い。なお、類似しているということは、図4a1および図4a2に示すプロット図4中の□印が、同じ雷から放射された信号を測定したこと、すなわち、雷発生位置が同じであることを意味している。ちなみに、図4a2の時刻t6付近のプロットに対応する測定点が、図4a1には存在していない。これは、雷発生位置が雷放射観測装置110cにやや近く、雷放射観測装置110aから遠くにあって、雷放射観測装置110aによって観測されなかったケースを表している。当該プロットは、統計処理を行う過程で外れ値として扱い、類似性の判定から除外することが可能である。また、図4bは、仰角βの時系列が類似しないと判定する場合を示している。すなわち、図4b1に示す雷放射観測装置110aでの仰角βの時系列の傾向例えば、仰角βの値を時系列に沿って結んだ線の形状は、時刻をずらしても、図4b2に示す他観測地点での仰角βの時系列の傾向と類似しているとは言えない。仮に、仰角βの値の時系列の傾向例えば、仰角βの値を時系列に沿って結んだ線の形状が一致していても、所定の時間範囲を超えて時刻がずれている場合には、別の雷と判定できる。図2に戻って、雷位置標定部25は、類似判定部24によって類似すると判定した識別情報の異なる2つ以上の仰角βの時系列と、当該仰角βの時系列とともに取得された方位角αの時系列および位置情報を用いて、雷発生位置を標定する。雷発生位置を標定する際には、類似すると判定された複数の方向ベクトルの中から2組の方向ベクトルを抽出し、それらの距離点線の長さが最小になるように演算し、その中点仮の雷発生位置を求める図5参照。さらに他の2組の方向ベクトルについてもそれらの距離点線の長さが最小になるように演算し、中点を算出する。そして、算出した中点の位置例えば、緯度経度の平均値を算出することによって、当該平均値の位置を雷発生位置として決定することができる。この方法の特徴は、複数の方向ベクトルから一括して雷の発生位置を求めるのではなく、2組ずつ総当りで仮の雷発生位置を求めることで、極端に推定位置が異なる結果をもたらす方向ベクトルを統計処理を用いて外れ値として排除することができる点である。このような外れ値となる方向ベクトルの判定方法として、例えば先に述べた2組の方向ベクトルから得られた最小距離を用いる方法がある。外れ値となる方向ベクトルが含まれた2組の方向ベクトルから得られる最小距離は、外れ値とならない2組の方向ベクトルから得られる最小距離に比較して、非常に大きな値になるため、外れ値となる方向ベクトルを容易に判定することができる。雷予測部26は、雲の動きを少なくとも含む気象情報データ風速、風向、気温、気圧等を外部のサーバ不図示等から取得して、雷位置標定部25によって標定した雷発生位置の移動を予測し、雷発生地域を推定する演算を実行する。記憶部23は、前記したアプリケーションプログラム、雷放射観測装置110から受信した情報例えば、受信情報301、および処理部21の処理結果等を記憶している。そして、記憶部23は、雷発生位置標定装置10に内蔵されていても、外付けされていても構わない。雷発生位置標定装置の処理フロー例次に、雷発生位置標定装置10の処理フロー例について、図6を用いて説明する適宜、図2参照。ステップS601では、処理部21の類似判定部24は、雷放射観測装置110から通信部22を介して受信した、仰角βおよび方位角αの時系列、その時系列の時刻、識別情報、および位置情報を観測データとして取得する。ステップS602では、処理部21の類似判定部24は、識別情報の異なる仰角βの時系列の間で、その時系列の時刻の間隔および仰角βの値が類似するか否かを判定する。ステップS603では、処理部21の雷位置標定部25は、類似判定部24によって類似すると判定した2つ以上の仰角βの時系列と、当該仰角βの時系列とともに取得された方位角αの時系列および位置情報とを用いて、雷発生位置を標定する。ここで、ステップS603によって標定された雷発生位置□印の時系列の一例を、雷の観測対象のエリア60とともに模式的に3次元表示すると、図7に示すようになる。図7にプロットした多数の雷発生位置は、図4a1のプロットを表したものである。なお、図7では、水平面の幅方向の軸は緯度、奥行き方向の軸は経度、高さ方向の軸は高度を表している。雷放射観測装置の配置次に、雷放射観測装置110の配置について、図8を用いて説明する。図8に示すように、雷の観測対象のエリア60が設定されたときに、破線で示す円状の領域13が、観測対象のエリア60をカバーするように配置される。ここで、カバーするとは、建造物や山等の陰やトンネル内等の例外的な場所もすべてをカバーするということではなく、観測対象としての範囲を含むようにするものである。なお、雷放射観測装置110は、円状の領域13の中心に設定される。なお、図8は、雷放射観測装置110を、地上ビルの屋上を含むの固定の観測地点に設置した場合を示している。さらに、雷放射観測装置110は、GPS受信機116を備えている場合には位置情報を随時取得できるので、車両や飛行機や船等の移動体に搭載されても構わない。移動体に雷放射観測装置110を搭載する場合には、雷放射観測装置110は、移動体の進行方向に応じて方位角αを補正し、高度に応じて仰角βを地上の数値に補正して、雷発生位置標定装置10に送信する。なお、前記したように、大多数の雷放電は、雷雲内および雷雲間で発生しており、特に、飛行機に雷放射観測装置110を搭載して雷から放射される電磁波を測定することは、雷発生位置を標定する際の精度向上のために、有効である。また、海上であれば、船に雷放射観測装置110を搭載して雷から放射される電磁波を測定することは、陸に近づいてくる雷雲の移動を的確に測定できるので、雷発生地域を推定する際の精度向上のために、有効である。変形例前記した実施形態では、雷発生位置標定装置10が、位置情報を雷放射観測装置110から取得する場合について説明したが、雷放射観測装置110が地上の固定の観測地点に設置されている場合には、位置情報を受信しなくても構わない。以下に、変形例として、雷発生位置標定装置10が、図9に示すように、受信情報901として、識別情報、仰角βおよび方位角αの時系列、および当該時系列の時刻を雷放射観測装置110から受信する場合について、前記した実施形態と異なる点に焦点を当てて説明する適宜、図1,2参照。なお、変形例における雷発生位置測定システムおよび雷発生位置標定装置は、前記した雷発生位置測定システム1および雷発生位置標定装置10と同様の構成であるため、同じ符号を用い、それらの図示は省略する。雷放射観測装置110が地上の固定の観測地点に設置されている場合には、雷放射観測装置110の位置情報が予め既知となっている。すなわち、雷発生位置標定装置10の記憶部23は、受信情報901と、識別情報と位置情報とを関連付けた装置位置情報1001図10参照とを記憶している。そして、雷発生位置標定装置10における処理フローは、図11に示すようになる。ステップS1101では、処理部21の類似判定部24は、通信部22を介して雷放射観測装置110から受信した、仰角βおよび方位角αの時系列、時系列の時刻、および識別情報を観測データとして取得する。ステップS1102では、処理部21の類似判定部24は、識別情報の異なる仰角βの時系列の間で、その時系列の時刻の間隔および仰角βの値が類似するか否かを判定する。ステップS1103では、処理部21の雷位置標定部25は、装置位置情報1001を参照して、類似判定部24によって類似すると判定した2つ以上の仰角βの時系列の識別情報に関連付けられている位置情報を求めるとともに、類似すると判定した2つ以上の仰角βの時系列と、当該仰角βの時系列とともに取得された方位角αの時系列と、識別情報に関連付けられている位置情報とを用いて、雷発生位置を標定する。以上、本実施形態および変形例における雷発生位置測定システム1は、各観測地点雷放射観測装置110を設定している地点において測定した仰角βの時系列の間で類似性を判定し、類似すると判定した仰角β、方位角αの時系列、および位置情報を用いることによって、雷発生位置を精度良く標定することができる。また、雷放射観測装置110を、所定の半径を有する円状の領域13の中心に設定し、その円状の領域13が観測対象のエリア60をカバーするように配置する。このことは、雷雲内および雷雲間の雷放電については、雷発生位置を標定する際の精度向上のために有効である。また、雷放射観測装置110を移動体に搭載することによって、観測対象のエリア60における雷発生位置の推定や雷雲の移動予測の精度向上を図ることができる。なお、雷発生位置標定装置10は、本実施形態に記載した雷発生位置標定装置10の機能と、変形例に記載した雷発生位置標定装置10の機能とを統合して、移動体の観測地点から位置情報を受信する場合および固定の観測地点から位置情報を受信しない場合の双方に対応可能な構成としても構わない。"