"宇宙線ミュー粒子を検出する第1の検出器と、前記宇宙線ミュー粒子を検出する第2の検出器と、前記第1の検出器および前記第2の検出器の両方が同じタイミングで前記宇宙線ミュー粒子を検出した場合に、該宇宙線ミュー粒子が該第1の検出器と該第2の検出器とを結ぶ方向を進行経路とするものであるとして計数し、進行経路ごとの前記宇宙線ミュー粒子の数を求めて出力する宇宙線ミュー粒子計数手段と、前記宇宙線ミュー粒子計数手段の出力である進行経路ごとの宇宙線ミュー粒子の数に基づいて、前記宇宙線ミュー粒子の位置分布および角度分布を求める宇宙線ミュー粒子分布取得手段と、前記分布取得手段で求めた宇宙線ミュー粒子の位置分布および角度分布に基づいて三次元のトモグラフィ解析を行うトモグラフィ解析手段と、を備え、前記第1の検出器と前記第2の検出器とを通過する前記宇宙線ミュー粒子のうち、少なくとも前記第1の検出器の中心と前記第2の検出器の中心とを結ぶ中心軸と平行な経路を進行経路とするものと、前記中心軸と所定の入射角を有する経路を進行経路とするものとを区別し、該区別した進行経路ごとの発生割合に基づいた重み係数を求める重み係数算出手段をさらに備え、前記宇宙線ミュー粒子計数手段が、前記重み係数算出手段によって求めた重み係数を用いて、前記第1の検出器と前記第2の検出器とを通過する進行経路ごとの前記宇宙線ミュー粒子の数に重み付けを行い、該重み付けの結果を出力する、ことを特徴とするミュー粒子を利用した三次元地盤探査システム。前記重み係数算出手段が、前記第1の検出器および前記第2の検出器の径および、前記第1の検出器と前記第2の検出器との距離を入力としたモンテカルロシミュレーションによって、前記第1の検出器と前記第2の検出器とを通過する前記宇宙線ミュー粒子の入射角依存性を求め、これによって前記重み係数を求める、ことを特徴とする請求項1に記載のミュー粒子を利用した三次元地盤探査システム。ところで、都市部におけるインフラストラクチャの老朽化による空洞を調べるには、ノイズに影響されず、地表部を占有しないで行うことが重要である。しかし、従来の調査方法を都市部で行うには、予測不可能な交通振動や工場の振動、あるいは送電線からの電気ノイズの影響により十分な成果をあげているとは言い難い。また、さらに地表部の道路などで行うには、その道路を占有するため交通規制が必須であり、交通渋滞などの発生が懸念される。特許文献1に記載の発明においても、予測不可能な交通振動や工場の振動の影響を受けるため、精度の高い探知をするのは困難であるという問題があった。特許文献2に記載のシンチレーション検出器を用いた地盤探査システムによれば、精度の高い探知が可能であるが、所定距離だけ離間して配置した2つのシンチレータを同時に通過した宇宙線ミュー粒子を計数するものであるため、2つのシンチレータの中心同士を結ぶ中心軸の方向に平行なものだけではなく、例えば各シンチレータの端部同士を対角に結ぶ方向など、中心軸に対して入射角を有する方向を進行経路とする宇宙線ミュー粒子も計数してしまう。こうなると、宇宙線ミュー粒子の進行経路を正確に取り扱うことができず、地盤探査システムの精度をより高めることが困難であるという問題があった。本発明は上記の点にかんがみてなされたもので、簡単な構成で、宇宙線ミュー粒子の進行経路を三次元で正確に検出し、地盤探査の精度をより高めることができるミュー粒子を利用した三次元地盤探査システムを提供することを目的とする。上記課題を解決するために、本発明は、宇宙線ミュー粒子を検出する第1の検出器と、前記宇宙線ミュー粒子を検出する第2の検出器と、前記第1の検出器および前記第2の検出器の両方が同じタイミングで前記宇宙線ミュー粒子を検出した場合に、該宇宙線ミュー粒子が該第1の検出器と該第2の検出器とを結ぶ方向を進行経路とするものであるとして計数し、進行経路ごとの前記宇宙線ミュー粒子の数を求めて出力する宇宙線ミュー粒子計数手段と、前記宇宙線ミュー粒子計数手段の出力である進行経路ごとの宇宙線ミュー粒子の数に基づいて、前記宇宙線ミュー粒子の位置分布および角度分布を求める宇宙線ミュー粒子分布取得手段と、前記分布取得手段で求めた宇宙線ミュー粒子の位置分布および角度分布に基づいて三次元のトモグラフィ解析を行うトモグラフィ解析手段と、を備え、前記第1の検出器と前記第2の検出器とを通過する前記宇宙線ミュー粒子のうち、少なくとも前記第1の検出器の中心と前記第2の検出器の中心とを結ぶ中心軸と平行な経路を進行経路とするものと、前記中心軸と所定の入射角を有する経路を進行経路とするものとを区別し、該区別した進行経路ごとの発生割合に基づいた重み係数を求める重み係数算出手段をさらに備え、前記宇宙線ミュー粒子計数手段が、前記重み係数算出手段によって求めた重み係数を用いて、前記第1の検出器と前記第2の検出器とを通過する進行経路ごとの前記宇宙線ミュー粒子の数に重み付けを行い、該重み付けの結果を出力する、ことを特徴とする。また本発明は、前記重み係数算出手段が、前記第1の検出器および前記第2の検出器の径および、前記第1の検出器と前記第2の検出器との距離を入力としたモンテカルロシミュレーションによって、前記第1の検出器と前記第2の検出器とを通過する前記宇宙線ミュー粒子の入射角依存性を求め、これによって前記重み係数を求める、ことを特徴とする。本発明によれば、簡単な構成で、宇宙線ミュー粒子の進行経路を三次元で正確に検出し、地盤探査の精度をより高めることができるミュー粒子を利用した三次元地盤探査システムを提供することができる。以下、本発明を実施するための形態について図面を参照して説明する。図1は本発明によるミュー粒子を利用した三次元地盤探査システムの一実施形態の構成を示す概略図である。図1を参照して分かるように、本実施形態の地盤探査システム1は、台車6に搭載された6つのシンチレータ2a、2b、2c、2d、2eおよび2fを有して構成され、この6つのシンチレータ2a~2fによって宇宙線ミュー粒子を検出する。6つのシンチレータ2a~2fのうちの1つであるシンチレータ2aは親シンチレータとも呼び、他のシンチレータ2b、2c、2d、2eおよび2fは子シンチレータとも呼ぶ。シンチレータ2b、2c、2d、2eおよび2fは、シンチレータ2aを中心とした放射状位置に配置している。中心に配置するシンチレータ2aは、例えば球状のプラスチックシンチレータであり、周囲に配置するシンチレータ2b~2fは、例えば円板円柱形状のプラスチックシンチレータであり、それぞれの円形の面がシンチレータ2aに向くように配置している。また、球体であるシンチレータ2aの直径と、円板円柱形状であるシンチレータ2b~2fの円形の面の直径とは、ほぼ等しい大きさにしている。本実施の形態では、シンチレータ2aの直径および、シンチレータ2b~2fの円形の面の直径は、それぞれ24.5cmである。シンチレータ2aの中心位置には、軸8を設けており、この軸8を中心として放射状に配置したフレーム7の先端にシンチレータ2b~2fを、配置している。フレーム7は軸8において回動可能に軸支されている。本実施の形態では、シンチレータ2aとシンチレータ2b~2fとの距離は1.4mであり、軸8を中心にフレーム7を回動し、シンチレータ2aに対するシンチレータ2b~2fの周方向位置を変えた場合も、このシンチレータ2aとシンチレータ2b~2fとの距離は一定に保たれる。シンチレータ2b~2fは、シンチレータ2aを中心として、例えば互いに15°の角度間隔で配置している。フレーム7を軸支する軸8は、台車6に搭載されて固定されており、台車6の例えば軸8近傍には、台車6の位置を検出する位置検出器4および、フレーム7の回動角度を検出する角度検出器5が搭載されている。フレーム7は、図示しないモータ等の駆動手段によって軸8を中心にして回動させられる。このフレーム7の回動角度は角度検出器5で検出され、図示しない制御手段では、この検出結果に基づきモータの角度位置のフィードバック制御を行う。すなわち、制御手段は、フレーム7をモータによって回動させ、フレーム7を所定角度位置で停止させることができる。角度検出器5としてはモータのエンコーダやジャイロ等を用いることができる。また、台車6は、図示しないモータ等の移動手段によって移動可能に構成されている。台車6の位置は位置検出器4によって検出される。位置検出器4としてはたとえばGPSを用いることができる。シンチレータ2a~2fは宇宙線ミュー粒子の通過に応じて発光し、この発光は光電子増倍管3a、3b、3c、3d、3eおよび3fのそれぞれによって検出され、その検出結果はマルチディスクリミネータ11によって分析される。すなわち、マルチディスクリミネータ11は、発光の検出結果のうち所定のエネルギーよりも高いものを分別して、宇宙線ミュー粒子が通過したときのエネルギーを判別する。マルチディスクリミネータ11は安定化電源10から電源供給される。マルチ同時計数回路12は、マルチディスクリミネータ11による判別結果に基づき、シンチレータ2aと、シンチレータ2b~2fのいずれかと、を通過した宇宙線ミュー粒子を計数する。宇宙線ミュー粒子は直進するため、フレーム7を所定角度位置で停止させ、シンチレータ2aと、シンチレータ2b~2fのいずれかと、を通過した宇宙線ミュー粒子を計数することにより、シンチレータ2aと、シンチレータ2b~2fのいずれかと、を結んだ方向に直進する宇宙線ミュー粒子の数を調べることができ、フレーム7の角度位置を変えてその都度宇宙線ミュー粒子を計数することにより、宇宙線ミュー粒子の角度分布を求めることができる。本実施の形態によれば、シンチレータ2aと2b、2aと2c、2aと2d、2aと2e、2aと2fのそれぞれの方向に直進する宇宙線ミュー粒子を1度に計数することができ、計測にかかる時間を短縮することができる。その後、フレーム7の角度位置を変える際には、上述のように子シンチレータ同士の角度間隔を15°にした場合、例えば0°5°10°と変化させることによって、すべてにおいて重複しない角度方向での宇宙線ミュー粒子を計数することができる。また、台車6の位置を変えてその都度宇宙線ミュー粒子を計数することにより、宇宙線ミュー粒子の位置分布を求めることができる。パソコン14では、収集したデータを保管するとともに、このような宇宙線ミュー粒子の各種分布を求めることができる。地盤探査システム1の台車6は、探査対象の地盤の下方のトンネル坑道等に配置し、探査を行う。安定化電源10、マルチディスクリミネータ11、マルチ同時計数回路12、パソコン14のそれぞれは、台車6に搭載されて台車6とともにトンネル坑道等に配置してもよいし、台車6とは別に、トンネル坑道等の外たとえば地上に配置してもよい。図2、図3は、図1に示した地盤探査システム1の構成の一例を示すブロック図である。上述のようにシンチレータ2a~2fとしては、たとえばプラスチックシンチレータを用いることができる。マルチディスクリミネータ11としては、光電子増倍管3a~3fのそれぞれによる検出結果を分析するディスクリミネータ11a~11fを設けるようにしてもよい。また、ディスクリミネータ11a~11fのそれぞれには、信号を増幅する増幅器を内蔵するようにしてもよい。ディスクリミネータ11a~11fに電源供給する安定化電源10は、たとえば高電圧回路および低電圧回路を有して構成される。また、図1に示したマルチ同時計数回路12は、図2に示すように、ディスクリミネータ11aおよびディスクリミネータ11bの両方からほぼ同時に宇宙線ミュー粒子を検出する信号が入力されたときに、検出信号を出力するマルチ同時計数回路12bと、マルチ同時計数回路12bからの検出信号を計数するカウンタ13bと、を有し、さらに、ディスクリミネータ11aおよびディスクリミネータ11cの両方からほぼ同時に宇宙線ミュー粒子を検出する信号が入力されたときに、検出信号を出力するマルチ同時計数回路12cと、マルチ同時計数回路12cからの検出信号を計数するカウンタ13cと、を有し、さらに、ディスクリミネータ11aおよびディスクリミネータ11dの両方からほぼ同時に宇宙線ミュー粒子を検出する信号が入力されたときに、検出信号を出力するマルチ同時計数回路12dと、マルチ同時計数回路12dからの検出信号を計数するカウンタ13dと、を有し、さらに、ディスクリミネータ11aおよびディスクリミネータ11eの両方からほぼ同時に宇宙線ミュー粒子を検出する信号が入力されたときに、検出信号を出力するマルチ同時計数回路12eと、マルチ同時計数回路12eからの検出信号を計数するカウンタ13eと、を有し、さらに、ディスクリミネータ11aおよびディスクリミネータ11fの両方からほぼ同時に宇宙線ミュー粒子を検出する信号が入力されたときに、検出信号を出力するマルチ同時計数回路12fと、マルチ同時計数回路12fからの検出信号を計数するカウンタ13fと、を有する構成であってもよい。図4は、図1に示した地盤探査システム1の原理を示すタイムチャートであり、aはシンチレータ2aによる検出信号を示し、bはシンチレータ2bによる検出信号を示し、cはシンチレータ2cによる検出信号を示し、dはマルチ同時計数回路12bの計数を示し、eはマルチ同時計数回路12cの計数を示す。例えば、図2に示す宇宙線ミュー粒子20がシンチレータ2aおよび2bを通過すると、その旨が、光電子増倍管3aおよび光電子増倍管3b、ならびにディスクリミネータ11aおよびディスクリミネータ11bを介してマルチ同時計数回路12bに入力され、マルチ同時計数回路12bでは、ディスクリミネータ11aからの検出信号とディスクリミネータ11bからの検出信号とがほぼ同時に報告されたことから、ひとつの宇宙線ミュー粒子20がシンチレータ2aおよび2bを通過したとみなし、宇宙線ミュー粒子20の検出信号をカウンタ13bに出力し、カウンタ13bではカウント数を1増加させる。これに対し、マルチ同時計数回路12bは、ディスクリミネータ11aからの検出信号またはディスクリミネータ11bからの検出信号の一方のみが報告された場合には、その宇宙線ミュー粒子20は、シンチレータ2aまたは2bの一方のみを通過した、すなわち、現在のフレーム7の角度位置でのシンチレータ2bの角度位置とは異なる角度のものであるとみなし、カウンタ13bによるカウント数は増加させない。マルチ同時計数回路12cないし12fも同様に動作し、カウンタ13bないし13fにおいて同様にカウントを行う。図5は、本発明による地盤探査システムの図1とは別の実施形態の構成を示す概略図であって、aは正面図であり、bは側面図である。図5において、図1に示した実施形態と同じ構成については同じ参照番号を付して詳しい説明を省略する。図5を参照して分かるように、本実施形態の地盤探査システム101は、全体を収容する枠体107aを有し、6つのシンチレータ2a、2b、2c、2d、2eおよび2fを有して構成され、この6つのシンチレータ2a~2fによって宇宙線ミュー粒子を検出する。シンチレータ2aの中心位置には、軸108aを設けており、この軸108aを中心として放射状に配置したフレーム107bの径方向外側にシンチレータ2b~2fを、配置している。フレーム107bは軸108aにおいて回動可能に軸支されている。フレーム107bを軸支する軸108aは、枠体107aに搭載されて固定されている。枠体107aの例えば右下部には、軸108a中心としてフレーム107bを回動させるモータ等の駆動手段108bを備えている。このフレーム107bの回動角度は角度検出器105で検出され、図示しない制御手段では、この検出結果に基づきモータの角度位置のフィードバック制御を行う。すなわち、制御手段は、フレーム107bを駆動手段108bによって回動させ、フレーム107bを所定角度位置で停止させることができる。角度検出器105としてはモータのエンコーダやジャイロ等を用いることができる。地盤探査システム101は、モータや車輪等の移動手段106を備え、枠体107aはこの移動手段106によって移動可能に構成されている。枠体107aの位置は、枠体107aの例えば左下部に設けた位置検出器104によって検出される。位置検出器104としてはたとえばGPSを用いることができる。なお、上述の実施の形態では、1つの親シンチレータと複数の子シンチレータとを備える構成を開示したが、本発明はこれに限られるものではなく、1つの親シンチレータと1つの子シンチレータとを備える構成であってもよいし、複数例えばm個の親シンチレータと複数例えばn個の子シンチレータとを備える構成であってもよい図6a参照。この構成によれば、測定時間を短縮し、より効率化することができる。また、親シンチレータや子シンチレータは、平面上に備えるものに限られず、たとえば球面上に備えるものであってもよい図6b参照。この構成によれば、例えば測定結果を三次元の表示する場合に、測定時間を短縮し、より効率化することができる。次に、本実施形態の地盤探査システム1の動作について説明する。図7は、シンチレータに入射するミュー粒子の透過経路の入射角について説明する図であって、aは入射角を示すシンチレータの概略側面図であり、bはシンチレータに入射するミュー粒子の透過経路の、入射角ごとの入射割合に基づいた重み係数を示すグラフである。本実施の形態では、2つのシンチレータをほぼ同時に通過したミュー粒子の透過経路を、その2つのシンチレータを結ぶ線としているが、2つのシンチレータを結ぶ線としては、2つの中心同士を結ぶ中心軸と平行な線のほかに、この中心軸と角度入射角を有する線も存在し、この入射角を有する線は、2つのシンチレータの断面積や、2つのシンチレータ同士の距離に依存する図7a参照。また、ミュー粒子の、入射角シンチレータの垂直からの角度に対するカウント数は一定ではなく、入射角依存性がある。すなわち、どの入射角の透過経路のミュー粒子も同じ数ではなく、入射角ごとに入射割合が異なる。そこで、本実施の形態では、2つのシンチレータの断面積や、2つのシンチレータ同士の距離を入力パラメータとして、モンテカルロシミュレーションを実行し入射角依存性を求め、これに基づいて重み係数を算出する。このミュー粒子の透過経路の入射角と重み係数との関係が、図7bに示すグラフである。本実施の形態では、2つのシンチレータをほぼ同時に通過したミュー粒子の総数を、図7bに示す重み係数を用いて、実際の透過経路に配分する。なお、図7では、二次元で表しているが、本実施の形態では、これを三次元で行う。図8は、シンチレータに入射するミュー粒子の透過経路の入射角について三次元で説明する図であって、aは入射角を示す概略斜視図であり、bはシンチレータに入射するミュー粒子の透過経路であって三次元の経路を点で示す親シンチレータの平面図である。また、図9は、シンチレータに入射するミュー粒子の透過経路を三次元で示す概略斜視図である。ミュー粒子のカウント総数から計算される面密度透過距離×密度は検出器への入射角によって計数率が異なる。すなわちカウント数から計算した面密度には入射角による重みが含まれている。本実施の形態では、この重みを、立体的に例えば1550方向の経路を仮定してトモグラフィ計算解析を行う。図10は、図1や図5、図6に示した本発明の実施形態に係る地盤探査システムによる測定範囲の一例を示す概略側断面図である。以下では、図1に示した地盤探査システム1を例に説明する。本実施形態では、地盤探査システム1を図10に示す坑道内の3つの観測点に順次配置し、各観測点において、フレーム7の角度位置を異ならせながら、各角度位置での宇宙線ミュー粒子のシンチレータ2aと、2b~2fのいずれかとの同時通過数をカウントする角度分布測定。宇宙線ミュー粒子は地盤を通過する際に、その数が減少する特性が知られており、地表からの距離が長ければ宇宙線ミュー粒子の数は減る。ところが、その地盤中の通過経路に空洞があると、その分だけ宇宙線ミュー粒子の数が減らなくなる。坑道内で測定位置観測点およびフレーム7の角度位置を変えることで、各位置、各方向での宇宙線ミュー粒子の数を得ることができ、これにより地盤の空洞との状態を解析するトモグラフィ解析に適用できるデータを取得することができる。これは、通常のトモグラフィを90度回転させたものと同様である。すなわち、本実施形態の地表と坑道が、通常のトモグラフィのボーリング孔に相当する。宇宙線ミュー粒子の物質中の透過特性は、物質の密度と透過距離との積および天頂角に依存する。この特性は既知であるので、ここでは詳しい説明は省略する。本実施形態の地盤探査システム1では、パソコン14においてトモグラフィ解析を行う。図11は、地盤中の空洞によるミュー粒子の角度分布依存性を示す図であり、aは地盤の様子を示す側断面図であり、bはaに示したA点の各角度位置での宇宙線ミュー粒子の同時通過数を示すグラフであり、cはaに示したB点の各角度位置での宇宙線ミュー粒子の同時通過数を示すグラフである。図11b、cのグラフは、縦軸が天頂角真上をゼロ度としたときのミュー粒子の到来方向角度であり、横軸がカウント率すなわち宇宙線ミュー粒子の同時通過数の1日の合計である。宇宙線ミュー粒子は、密度と透過距離が増えると透過カウント率が減少する。逆に空洞があると同じ透過距離でも密度が小さいので透過率が増加する。地盤内に空洞がなければ、図11bおよびcに実線で示すような特性を示すわけだが、図11aに示すような空洞があると、宇宙線ミュー粒子がその空洞を通る角度位置では、空洞がない場合よりも同時通過数が増える。すなわち、均質な地盤の場合、角度に応じて同時計数が変化するが、空洞のある場所では、宇宙線ミュー粒子の透過性が良いため、計数率が増加することが確認できる。図11bに示すように、A点では、天頂角25度から55度の範囲で、破線で示すように透過率が増加し、図11cに示すように、B地点では天頂角30度から50度の範囲で、破線で示すように透過率が増加している。このように位置を変えて観測することにより、空洞の位置をある程度まで絞り込むことができる。この図11a、b、cに示す例では、2点のみの測定結果を示しているが、より多くの測定点でデータを収集することにより、より精度の高いトモグラフィ解析が可能となる。図12は、本実施形態の地盤探査システム1のパソコン14で実行されるトモグラフィ解析のフローチャートを示す図である。本実施形態の地盤探査システム1では、台車6やフレーム7を動かして各位置でのデータを取得したならば、まず初期モデルを設定するS-1。ここで、初期モデルとは、事前に調査した地質断面図や土質試験の結果などに基づいて地盤の平均密度を仮定した、地盤の層構造のモデルである。続いて、地盤が、ステップS-1で設定したモデルであるとした場合に、観測されるであろう宇宙線ミュー粒子の理論値を計算するS-2。本実施の形態では、この理論値の計算を三次元で行う。続いて、ステップS-2で計算した理論値と、すでに取得した各位置でのデータとの差残差を計算するS-3。続いて、ステップS-3で計算した残差があらかじめ定めた許容範囲内であればS-4:Yes、その結果を表示してS-6処理を終了する。ステップS-4において、ステップS-3で計算した残差があらかじめ定めた許容範囲内でない場合にはS-4:No、モデルを修正してS-5、ステップS-2に戻り、処理を繰り返す。本実施形態の地盤探査システム1によれば、宇宙線ミュー粒子の位置分布、角度分布測定技術と、トモグラフィの解析技術とを組み合わせることにより、地盤の密度を可視化することが可能となる。図10は、地盤を可視化したものの一例であり、図10の例では地盤の地質密度が均一で空洞がない場合を示している。このような可視化の別の例について図13を参照して説明する。図13aは初期モデルの例を示す側断面図であり、図13bは位置と角度を変えて宇宙線ミュー粒子の計数を行う、すなわち宇宙線ミュー粒子の透過率を測定する様子を示す側断面図であり、図13cはトモグラフィ解析を行った結果である地盤の密度分布を示す側断面図であり、図13dは複数の断面で測定することによる三次元可視化への応用を示す斜視図である。図13aは通常の地質調査によって得られる断面であり、なんらかの理由たとえばトンネルの老朽化により、空洞の発生が懸念される。そこで、図13aを初期モデルとして、上述の宇宙線ミュー粒子の位置分布、角度分布測定およびトモグラフィ解析を行う図13b、図13c参照。また、上述の図12のステップS-2において三次元の理論値の計算をしたり、複数の断面で測定したりすることで図13dに示すように三次元への拡張が可能である。以下、三次元への適用結果について詳しく説明する。図14は、本実施の形態による地盤探査の三次元解析の結果を示す図であり、aは三次元表示したものであり、b、c、dは、aの断面図である。解析の条件は以下のとおりである。地盤の平均密度:20gcm3土被り:10m縦断方向:20m横断方向:10m空洞の大きさ:直径2m空洞の位置:直方体の中心本実施の形態による地盤探査の三次元解析によれば、図14a~dに示すように、空洞の中央部の低密度が再現された。また、本実施の形態の地盤探査システムを用いて三次元トモグラフィデータを取得するため,地下で実験を行った結果を以下に示す。図15は、地下での実験の実験場所を示す図であり、aは平面図であり、bは側断面図である。図15a、bに示すように、この実験場所は、一番深いところに計測通路があり,地盤の厚さは11mである。その上部に通路と下水管がある。通路は、右端で折れ曲がっている。その下で,縦断方向15m、横断方向5m、鉛直方向11mの直方体部分を透過するミュー粒子をカウントした。本実施の形態により取得した三次元トモグラフィデータを用いて解析した結果を図16、図17、図18に示す。図16は、地下での実験の三次元トモグラフィ解析結果を示す図である。図17は、地下での実験の解析結果を示す図であって、aは図15bの断面図に相当する位置の解析結果であり、bはaを下水管の位置で切断して上から見た解析結果であり、cはaを通路の位置で切断して上から見た解析結果である。図18は、地下での実験の三次元解析結果を示す図であって、内部解析のため一部を切断して示しており、aは下水管の側から見た解析結果であり、bは通路の側から見た解析結果である。計測範囲の大きさは縦断方向15m、横断方向5m、鉛直方向11mである。図16、図17および図18を参照して分かるように、本発明によれば、三次元解析を可能とし、地盤の密度分布様子を三次元で可視化することが可能である。以上、本発明を説明したが、本発明は、この説明に限定されるものではなく、本発明の趣旨を逸脱しない範囲で数々の変形および組み合わせが出来ることは勿論である。"