ホーム > 特許ランキング > 株式会社カナン・ジオリサーチ
知財求人 - 知財ポータルサイト「IP Force」
▶ 株式会社カナン・ジオリサーチの特許一覧 ▶ 株式会社岩根研究所の特許一覧
(19)【発行国】日本国特許庁(JP)
(12)【公報種別】特許公報(B2)
(11)【特許番号】
(24)【登録日】2022-10-20
(45)【発行日】2022-10-28
(54)【発明の名称】含水比マッピング方法及び含水比マッピング装置
(51)【国際特許分類】
G01V 3/12 20060101AFI20221021BHJP
G01S 13/88 20060101ALI20221021BHJP
G01S 7/41 20060101ALI20221021BHJP
【ＦＩ】
G01V3/12 B
G01S13/88 200
G01S7/41
【請求項の数】
5
(21)【出願番号】P 2018118944
(22)【出願日】2018-06-22
(65)【公開番号】P2019219348
(43)【公開日】2019-12-26
【審査請求日】2021-03-17
(73)【特許権者】
【識別番号】516250753
【氏名又は名称】株式会社カナン・ジオリサーチ
(73)【特許権者】
【識別番号】393021141
【氏名又は名称】株式会社岩根研究所
(74)【代理人】
【識別番号】110002354
【氏名又は名称】弁理士法人平和国際特許事務所
(72)【発明者】
【氏名】篠原 潤
(72)【発明者】
【氏名】加藤 祐悟
(72)【発明者】
【氏名】エルヴィス アヌープ シュクラ
(72)【発明者】
【氏名】馬場 務
(72)【発明者】
【氏名】岩根 和郎
【審査官】福田 裕司
(56)【参考文献】
【文献】特開２０１８－０２８５０７（ＪＰ，Ａ）
【文献】特開２００５－２９２０６９（ＪＰ，Ａ）
【文献】特開２００８－１２８９７４（ＪＰ，Ａ）
【文献】特開２００３－３０２４６５（ＪＰ，Ａ）
【文献】特開２００３－１０７１６９（ＪＰ，Ａ）
【文献】特開２０１４－０９８５９８（ＪＰ，Ａ）
(58)【調査した分野】(Int.Cl.，ＤＢ名)
Ｇ０１Ｖ ３／１２
Ｇ０１Ｓ １３／８８
Ｇ０１Ｓ ７／４１
(57)【特許請求の範囲】
【請求項1】
送信アンテナと受信アンテナが同一の中点と異なるオフセット距離を有する、前記オフセット距離方向に沿って前記中点が異なる複数の送受信アンテナ組を備えた地中レーダを用いて地中の含水性状のマッピングを行う含水比マッピング方法であって、前記地中レーダを地表面に沿って移動させながら地中へ向けて電磁波を放射することにより、前記送受信アンテナのオフセット距離方向に沿って前記中点の異なる前記複数の各送受信アンテナ組によって得られる複数の反射波を、前記地中レーダの進行方向に沿った複数の測定点において取得するレーダ反射波取得工程と、
前記複数の各送受信アンテナ組において送信アンテナから受信アンテナへの到達時間が最短の電磁波を直接波として、その直接波の到達時刻を抽出する直接波フィッティング工程と、
前記直接波の到達時刻から、当該直接波の地表面における到達時刻を時刻ゼロとして算出するゼロ時間補正工程と、
前記複数の各送受信アンテナ組において送信アンテナから受信アンテナへ前記直接波の次に到達する電磁波を反射波として、その反射波の到達時刻を抽出する反射波フィッティング工程と、
前記反射波の到達時刻に基づいて、前記地中の比誘電率を算出する比誘電率算出工程と、
前記比誘電率に基づいて、前記地中の含水比を算出する含水比算出工程と、
前記複数の各送受信アンテナ組に対応する複数の前記含水比に基づいて、含水比分布図を生成する含水比分布図生成工程と、を備え、
前記複数の各送受信アンテナ組のうち、前記送受信アンテナのオフセット距離方向に沿って中点の異なる前記送受信アンテナの組で得られる複数の反射波を、前記地中レーダの進行方向に沿った同一ルートにおける前記複数の測定点となる各中点における測定結果として前記含水比を算出し、
前記同一ルートにおける複数の各中点における前記地中の含水比を、当該複数の各中点の地上と対応した位置情報と対応付けることにより、前記送受信アンテナのオフセット距離方向及び前記地中レーダの進行方向に沿った複数の各中点における含水比の分布を示す前記含水比分布図を生成する
ことを特徴とする含水比マッピング方法。
【請求項2】
送信アンテナと受信アンテナが同一の中点と異なるオフセット距離を有する、前記オフセット距離方向に沿って前記中点が異なる複数の送受信アンテナ組を備えた地中レーダと、前記地中レーダを地表面に沿って移動させながら地中へ向けて電磁波を放射することにより、前記送受信アンテナのオフセット距離方向に沿って前記中点の異なる前記複数の各送受信アンテナ組によって得られる複数の反射波を、前記地中レーダの進行方向に沿った複数の測定点において取得するレーダ反射波取得部と、
前記複数の各送受信アンテナ組において送信アンテナから受信アンテナへの到達時間が最短の電磁波を直接波として、その直接波の到達時刻を抽出する直接波フィッティング部と、
前記直接波の到達時刻から、当該直接波の地表面における到達時刻を時刻ゼロとして算出するゼロ時間補正部と、
前記複数の各送受信アンテナ組において送信アンテナから受信アンテナへ前記直接波の次に到達する電磁波を反射波として、その反射波の到達時刻を抽出する反射波フィッティング部と、
前記反射波の到達時刻に基づいて、前記地中の比誘電率を算出する比誘電率算出部と、
前記比誘電率に基づいて、前記地中の含水比を算出する含水比算出部と、
前記複数の各送受信アンテナ組に対応する複数の前記含水比に基づいて、含水比分布図を生成する含水比分布図生成部と、を備え、
前記複数の各送受信アンテナ組のうち、前記送受信アンテナのオフセット距離方向に沿って中点の異なる前記送受信アンテナの組で得られる複数の反射波を、前記地中レーダの進行方向に沿った同一ルートにおける前記複数の測定点となる各中点における測定結果として前記含水比を算出し、
前記同一ルートにおける複数の各中点における前記地中の含水比を、当該複数の各中点の地上と対応した位置情報と対応付けることにより、前記送受信アンテナのオフセット距離方向及び前記地中レーダの進行方向に沿った複数の各中点における含水比の分布を示す前記含水比分布図を生成する
ことを特徴とする含水比マッピング装置。
【請求項3】
前記地中レーダが、路面を走行可能な探査車両の底面の、前輪と後輪の間、又は複数の車輪の車軸間に配置される、複数の送信アンテナ及び受信アンテナを備えたマルチチャンネル地中レーダからなる
ことを特徴とする請求項２記載の含水比マッピング装置。
【請求項4】
路面を走行可能な探査車両が走行する路面上の地上映像を生成する地上映像生成部と、前記含水比分布図生成部で生成される含水比分布図と、前記地上映像生成部で生成される地上映像を、地中及び地上の位置情報が一元化された一元化情報として生成する地中・地上情報一元化処理部と、を備える
ことを特徴とする請求項２又は３記載の含水比マッピング装置。
【請求項5】
前記地上映像生成部が、前記探査車両が走行する路面上の地上三次元映像の画像データから、所定数の特徴点を自動抽出する特徴点抽出部と、
抽出された特徴点について、動画映像の各フレーム画像内で自動追跡してフレーム画像間での対応関係を求める特徴点対応処理部と、
対応関係が求められた特徴点の三次元位置座標を求め、当該三次元位置座標から、各フレーム画像に対応したカメラの三次元位置座標及び三次元回転座標からなるカメラベクトルを求めるカメラベクトル演算部と、
を有するＣＶ演算部を備え、
前記地中・地上情報一元化処理部が、
前記カメラベクトルと、前記含水比分布図が有する位置情報とに基づいて、前記地上三次元映像と前記含水比分布図を一元化する
ことを特徴とする請求項４記載の含水比マッピング装置。
【発明の詳細な説明】
【技術分野】
【0001】
本発明は、地中レーダ探査に代表される、電磁波や地震波，電気比抵抗，超音波等に対する地盤の物理的性質を利用して、地中（地下）の状況を非破壊で間接的に探査・調査する地中探査技術に関し、特に、地中の土壌の含水比を正確かつ効率的に測定するための技術に関する。
【背景技術】
【0002】
近年、道路やトンネル，橋梁，港湾施設などの社会インフラの老朽化に伴い、それらの施設を補修して寿命を延ばし、安全な社会環境を維持することが求められており、そのために道路等の構造物の調査・診断を行うための地中探査技術が極めて重要となっている。
特に、地震や異常気象，老朽化等による地中や構造物の異常が増加しており、例えばコンクリート内部の劣化状況や道路路面下などの空洞調査を正確・迅速に行うことが急務となっている。
このため、地中探査技術の重要性は益々高まっている。
特に、地震や異常気象，老朽化等による地中や構造物の異常が増加しており、例えばコンクリート内部の劣化状況や道路路面下などの空洞調査を正確・迅速に行うことが急務となっている。
このため、地中探査技術の重要性は益々高まっている。
【0003】
地中探査は、例えば舗装道路下における空洞探査の場合は、地中に向けて電磁波を投射して地中からの反射波を受信する地中レーダ等のセンサを搭載した探査車両を用意し、探査車両を一般車両と同様に道路を走行させながらセンサで得られた探査データを記録し、記録した探査データを解析することで、例えば、地中に存在する空洞や亀裂，埋設管などの存在や位置を特定・把握し、また、地中の土壌の含水比や比誘電率等の性状の測定、性状の異なる地層の境界面の検知などが行われる。
とりわけ近年は、地中の埋設管の老朽化等に伴う漏水件数の増加などが見込まれることから、埋設管等の特定とその周囲の含水比の測定などの地中探査技術の重要性が高まっている。
とりわけ近年は、地中の埋設管の老朽化等に伴う漏水件数の増加などが見込まれることから、埋設管等の特定とその周囲の含水比の測定などの地中探査技術の重要性が高まっている。
【0004】
ここで、地中の含水比は、土壌の比誘電率によって求めることができる。
このような地中の含水比の測定を目的として、例えば特許文献１には、地中の比誘電率を測定する方法が提案されている。
この特許文献１で提案されている方法は、土壌に埋設された深さが既知である対象物が存在する地中に向かって、地面に配置された送信アンテナから地中レーダの電磁波を放射し、送信アンテナと一定のアンテナ間隔で地面に配置された受信アンテナで電磁波の反射波を受信し、対象物からの反射波の到達時間を判断して、到達時間の判断の結果と対象物の深さとアンテナ間隔とを使用して土壌の比誘電率を算出するというものである。
このような地中の含水比の測定を目的として、例えば特許文献１には、地中の比誘電率を測定する方法が提案されている。
この特許文献１で提案されている方法は、土壌に埋設された深さが既知である対象物が存在する地中に向かって、地面に配置された送信アンテナから地中レーダの電磁波を放射し、送信アンテナと一定のアンテナ間隔で地面に配置された受信アンテナで電磁波の反射波を受信し、対象物からの反射波の到達時間を判断して、到達時間の判断の結果と対象物の深さとアンテナ間隔とを使用して土壌の比誘電率を算出するというものである。
【先行技術文献】
【特許文献】
【0005】
【文献】特開２０１６－２１１８９７号公報
【発明の概要】
【発明が解決しようとする課題】
【0006】
しかしながら、特許文献１で提案されている比誘電率の測定方法では、地中レーダを構成する送受信アンテナは、地表面の所定箇所に固定的に配置されるものであり、したがって、比誘電率が測定できる範囲も、送受信アンテナが配置された地表面に対応する地中の一定範囲に限られてしまう。このため、広範囲に広がる道路やその下に設置される埋設管などを監視・測定することはできなかった。
また、特許文献１で用いられているような地中レーダを、例えば手押し型や牽引型の移動可能な構成とすることも考えられるが、道路などの広範囲な地中探査を行うことは、非常に手間と時間がかかり、実際の地中探査に適用することは困難乃至不可能であった。
また、特許文献１で用いられているような地中レーダを、例えば手押し型や牽引型の移動可能な構成とすることも考えられるが、道路などの広範囲な地中探査を行うことは、非常に手間と時間がかかり、実際の地中探査に適用することは困難乃至不可能であった。
【0007】
また、地中探査においては、地中レーダ等のセンサを探査車両によって牽引することにより、探査車両を一般車両と同様に道路を走行させながら探査データを記録し、記録した探査データを解析することで、地中に存在する空洞や亀裂，埋設管などの存在や位置の特定、土壌の含水比の測定などの一次調査が行われる。
そして、一次調査で得られた探査データの解析の結果、地中に例えば空洞が生じていると思われる箇所があれば、その周辺領域を詳細に調査し、空洞の位置や大きさ，規模等を正確に把握するための二次調査が行われ、陥没の危険性等が判定・評価される。その後、道路管理者等によって、空洞の補修等の必要な措置が取られることになる。
そして、一次調査で得られた探査データの解析の結果、地中に例えば空洞が生じていると思われる箇所があれば、その周辺領域を詳細に調査し、空洞の位置や大きさ，規模等を正確に把握するための二次調査が行われ、陥没の危険性等が判定・評価される。その後、道路管理者等によって、空洞の補修等の必要な措置が取られることになる。
【0008】
このような地中探査の一次調査で得られる探査データは、探査車両の走行方向における道路下の縦断面や横断面を波形データや、土壌の含水比データとして得るものであるが、データ中に示される地中の位置や対象物（例えば空洞や埋設管など）が、地表面から見てどの位置に存在するかを特定することが必要となる。従来の地中探査技術では、一次調査で得られた探査データに基づいて、二次調査を実施するにあたって、手押し型のセンサ等を用いた現地での位置特定作業や、地上の構造物の形状や道路周辺物等の計測作業等を行う必要があり、そのような位置特定作業等が大きな負担となっていた。
このため、特許文献１で提案されているような測定方法により土壌の比誘電率に基づく含水比が求められたとしても、地上との位置特定が正確に行われなければ、地中探査のデータとして活用することは困難であった。
このため、特許文献１で提案されているような測定方法により土壌の比誘電率に基づく含水比が求められたとしても、地上との位置特定が正確に行われなければ、地中探査のデータとして活用することは困難であった。
【0009】
ここで、地中探査データに対して地上と対応した位置情報を付与するには、地中探査データを取得する探査車両に装備されるＧＰＳ情報を用いることが考えられる。しかしながら、ＧＰＳで得られる位置情報は、数メートル～１０数メートルの誤差がある。
このため、地中探査データと地表面の位置を正確に精度良く対応させるためには、ＧＰＳを用いた位置情報の取得・生成では不十分であり、特に、空洞や埋設管等の位置をセンチメートルの範囲で特定する必要のある地中探査技術では、より正確で信頼性の高い位置情報の取得が望まれていた。
このため、地中探査データと地表面の位置を正確に精度良く対応させるためには、ＧＰＳを用いた位置情報の取得・生成では不十分であり、特に、空洞や埋設管等の位置をセンチメートルの範囲で特定する必要のある地中探査技術では、より正確で信頼性の高い位置情報の取得が望まれていた。
【0010】
本発明は、以上のような従来の技術が有する課題を解決するために提案されたものであり、地中の土壌の含水比を広範囲にわたって正確かつ効率的に測定できるとともに、測定された含水比データに対して、ＧＰＳによることなく、より正確で高精度な位置情報を付与することができ、それによって地中探査データと地表面の情報を一元化された単一の情報として生成・表示させることができる含水比マッピング方法及び含水比マッピング装置の提供を目的とする。
【課題を解決するための手段】
【0011】
上記目的を達成するため、本発明の含水比マッピング方法は、送信アンテナと受信アンテナが同一の中点と異なるオフセット距離を有する、前記オフセット距離方向に沿って前記中点が異なる複数の送受信アンテナ組を備えた地中レーダを用いて地中の含水性状のマッピングを行う含水比マッピング方法であって、前記地中レーダを地表面に沿って移動させながら地中へ向けて電磁波を放射することにより、前記送受信アンテナのオフセット距離方向に沿って前記中点の異なる前記複数の各送受信アンテナ組によって得られる複数の反射波を、前記地中レーダの進行方向に沿った複数の測定点において取得するレーダ反射波取得工程と、前記複数の各送受信アンテナ組において送信アンテナから受信アンテナへの到達時間が最短の電磁波を直接波として、その直接波の到達時刻を抽出する直接波フィッティング工程と、前記直接波の到達時刻から、当該直接波の地表面における到達時刻を時刻ゼロとして算出するゼロ時間補正工程と、前記複数の各送受信アンテナ組において送信アンテナから受信アンテナへ前記直接波の次に到達する電磁波を反射波として、その反射波の到達時刻を抽出する反射波フィッティング工程と、前記反射波の到達時刻に基づいて、前記地中の比誘電率を算出する比誘電率算出工程と、前記比誘電率に基づいて、前記地中の含水比を算出する含水比算出工程と、前記複数の各送受信アンテナ組に対応する複数の前記含水比に基づいて、含水比分布図を生成する含水比分布図生成工程と、を備え、前記複数の各送受信アンテナ組のうち、前記送受信アンテナのオフセット距離方向に沿って中点の異なる前記送受信アンテナの組で得られる複数の反射波を、前記地中レーダの進行方向に沿った同一ルートにおける前記複数の測定点となる各中点における測定結果として前記含水比を算出し、前記同一ルートにおける複数の各中点における前記地中の含水比を、当該複数の各中点の地上と対応した位置情報と対応付けることにより、前記送受信アンテナのオフセット距離方向及び前記地中レーダの進行方向に沿った複数の各中点における含水比の分布を示す前記含水比分布図を生成する構成としてある。
【0012】
また、本発明の含水比マッピング装置は、送信アンテナと受信アンテナが同一の中点と異なるオフセット距離を有する、前記オフセット距離方向に沿って前記中点が異なる複数の送受信アンテナ組を備えた地中レーダと、前記地中レーダを地表面に沿って移動させながら地中へ向けて電磁波を放射することにより、前記送受信アンテナのオフセット距離方向に沿って前記中点の異なる前記複数の各送受信アンテナ組によって得られる複数の反射波を、前記地中レーダの進行方向に沿った複数の測定点において取得するレーダ反射波取得部と、前記複数の各送受信アンテナ組において送信アンテナから受信アンテナへの到達時間が最短の電磁波を直接波として、その直接波の到達時刻を抽出する直接波フィッティング部と、前記直接波の到達時刻から、当該直接波の地表面における到達時刻を時刻ゼロとして算出するゼロ時間補正部と、前記複数の各送受信アンテナ組において送信アンテナから受信アンテナへ前記直接波の次に到達する電磁波を反射波として、その反射波の到達時刻を抽出する反射波フィッティング部と、前記反射波の到達時刻に基づいて、前記地中の比誘電率を算出する比誘電率算出部と、前記比誘電率に基づいて、前記地中の含水比を算出する含水比算出部と、前記複数の各送受信アンテナ組に対応する複数の前記含水比に基づいて、含水比分布図を生成する含水比分布図生成部と、を備え、前記複数の各送受信アンテナ組のうち、前記送受信アンテナのオフセット距離方向に沿って中点の異なる前記送受信アンテナの組で得られる複数の反射波を、前記地中レーダの進行方向に沿った同一ルートにおける前記複数の測定点となる各中点における測定結果として前記含水比を算出し、前記同一ルートにおける複数の各中点における前記地中の含水比を、当該複数の各中点の地上と対応した位置情報と対応付けることにより、前記送受信アンテナのオフセット距離方向及び前記地中レーダの進行方向に沿った複数の各中点における含水比の分布を示す前記含水比分布図を生成する構成としてある。
【発明の効果】
【0013】
本発明によれば、地中の土壌の含水比を広範囲にわたって正確かつ効率的に測定することができる。そして、測定された含水比データに対して、ＧＰＳによることなく、より正確で高精度な位置情報を付与することができる。
これによって、地中探査データと地表面の情報を一元化された単一の情報として生成・表示させることが可能となり、例えば、地中の土壌の含水比の分布を正確に示すマッピング情報を、地表面を示す正確なオルソ画像とともに一体的に表示・出力させることが可能となる。
したがって、本発明によれば、従来にはない正確・高精度な地中探査を実現することができる。
これによって、地中探査データと地表面の情報を一元化された単一の情報として生成・表示させることが可能となり、例えば、地中の土壌の含水比の分布を正確に示すマッピング情報を、地表面を示す正確なオルソ画像とともに一体的に表示・出力させることが可能となる。
したがって、本発明によれば、従来にはない正確・高精度な地中探査を実現することができる。
【図面の簡単な説明】
【0014】
【図1】本発明の一実施形態に係る含水比マッピング方法を実施するための地中探査の全体像を模式的に示す説明図である。
【図2】本発明の一実施形態に係る含水比マッピング装置を構成する、地中探査データを取得・生成する探査車両の動作イメージを模式的に示す説明図であり、（ａ）は探査車両が走行する路面の断面図、（ｂ）は同じく探査車両が走行する路面の平面図である。
【図4】探査車両と、取得される含水比データの関係を模式的に示す説明図であり、（ａ）は探査車両の進行方向に沿って取得される含水比データを、（ｂ）は探査車両の進行方向と交差する車幅方向に沿って取得される含水比データを、（ｃ）は（ａ）及び（ｂ）に示す含水比データが合成されて生成される含水比の境界面（含水比マッピング）データを示している。
【図5】本発明の一実施形態に係る含水比マッピング方法の各工程を示すフローチャートである。
【図6】本発明の一実施形態に係る含水比マッピング装置を構成する送受信アンテナの中点・オフセット距離の関係を模式的に示す図であり、（ａ）は送受信アンテナの断面正面図、（ｂ）は中点と複数の送受信アンテナ組の位置を示す説明図である。
【図7】オフセット距離が異なる複数の送受信アンテナ組によって取得される直接波・反射波の到達時間とオフセット距離の関係を示すグラフである。
【図8】中点の位置を複数異ならせた場合の、中点と複数の送受信アンテナ組の位置を示す説明図である。
【図9】探査車両（地中レーダ）の移動に伴って移動する中点と複数の送受信アンテナ組の位置を示す図であり、（ａ）は探査車両が直進している場合、（ｂ）は探査車両が路面のカーブに沿って移動している場合を示している。
【図10】探査車両（地中レーダ）の移動に伴って移動する中点と複数の送受信アンテナ組の位置を示す図であり、（ａ）は探査車両が同じ路面を複数回異なるルートで移動した場合を、（ｂ）は（ａ）の図から中点のみを代表して表示させた場合を示している。
【図11】本発明の一実施形態に係る含水比マッピング方法（装置）によって取得される含水比分布図（含水比マッピング）を模式的に示す説明図であり、（ａ）は図１０（ｂ）に示す中点のみを表示させた場合、（ｂ）は中点とともに含水比の分布を等高線状に示した場合、（ｃ）は（ｂ）の図から中点を取り除いて等高線のみを示した場合である。
【図12】本発明の一実施形態に係る含水比マッピング装置で取得される、地上３６０度の全周立体映像の撮像範囲を模式的に示す説明図であり、（ａ）は路面・橋梁等における撮像範囲を示す、（ｂ）は道路側方の斜面（法面）の撮像範囲を示す、それぞれ路面の正面断面図である。また、（ｃ）は３６０度全周映像の一例である。
【図13】本発明の一実施形態に係る含水比マッピング方法（装置）における含水比分布図データと地上映像の一元化処理の手順を示すフローチャートである。
【図14】本発明の一実施形態に係る含水比マッピング方法（装置）で取得・生成される含水比分布図データと地上オルソ画像に一元化処理した画像のイメージ図である。
【図15】本発明の一実施形態に係る含水比マッピング装置において地上３６０度立体映像から生成されるカメラベクトル画像を生成するＣＶ（カメラベクトル）演算部の基本構成を示すブロック図である。
【図16】地上３６０度立体映像から得られる変換画像を示す説明図であり、（ａ）は球面画像が貼り付けられる仮想球面を、（ｂ）は仮想球面に貼り付けられた球面画像の一例を、（ｃ）は（ｂ）に示した球面画像をメルカトール図法に従って平面展開した画像を示している。
【図17】本発明の一実施形態に係るＣＶ演算部における具体的なカメラベクトルの検出方法を示す説明図である。
【図18】本発明の一実施形態に係るＣＶ演算部における具体的なカメラベクトルの検出方法を示す説明図である。
【図19】本発明の一実施形態に係るＣＶ演算部における具体的なカメラベクトルの検出方法を示す説明図である。
【図20】本発明の一実施形態に係るＣＶ演算部によるカメラベクトルの検出方法における望ましい特徴点の指定態様を示す説明図である。
【図21】本発明の一実施形態に係るＣＶ演算部により得られる特徴点の三次元座標とカメラベクトルの例を示すグラフである。
【図22】本発明の一実施形態に係るＣＶ演算部により得られる特徴点の三次元座標とカメラベクトルの例を示すグラフである。
【図23】本発明の一実施形態に係るＣＶ演算部により得られる特徴点の三次元座標とカメラベクトルの例を示すグラフである。
【図24】本発明の一実施形態に係るＣＶ演算部において、カメラから特徴点の距離に応じて複数の特徴点を設定し、複数の演算を繰り返し行う場合を示す説明図である。
【図25】本発明の一実施形態に係るＣＶ演算部求められたカメラベクトルの軌跡をビデオ映像中に表示した場合の図である。
【発明を実施するための形態】
【0015】
以下、本発明に係る含水比マッピング方法及び含水比マッピング装置の実施形態について、図面を参照しつつ説明する。
ここで、以下に示す本発明の含水比マッピング方法及び含水比マッピング装置における地中探査データ及び地上三次元映像の生成・出力・一元化処理・カメラベクトル演算等は、プログラム（ソフトウェア）の命令によりコンピュータで実行される処理，手段，機能によって実現される。プログラムは、コンピュータの各構成要素に指令を送り、以下に示す本発明に係る所定の処理や機能等を行わせることができる。
すなわち、本発明における各処理や手段，機能は、プログラムとコンピュータとが協働した具体的手段によって実現される。
ここで、以下に示す本発明の含水比マッピング方法及び含水比マッピング装置における地中探査データ及び地上三次元映像の生成・出力・一元化処理・カメラベクトル演算等は、プログラム（ソフトウェア）の命令によりコンピュータで実行される処理，手段，機能によって実現される。プログラムは、コンピュータの各構成要素に指令を送り、以下に示す本発明に係る所定の処理や機能等を行わせることができる。
すなわち、本発明における各処理や手段，機能は、プログラムとコンピュータとが協働した具体的手段によって実現される。
【0016】
なお、プログラムの全部又は一部は、例えば、磁気ディスク，光ディスク，半導体メモリ，その他任意のコンピュータで読取り可能な記録媒体により提供され、記録媒体から読み出されたプログラムがコンピュータにインストールされて実行される。
また、プログラムは、記録媒体を介さず、通信回線を通じて直接にコンピュータにロードし実行することもできる。また、本発明に係る含水比マッピング装置に備えられる、単一の情報処理装置（例えば１台のパーソナルコンピュータ等）で構成することもでき、複数の情報処理装置（例えば複数台のコンピュータ群等）で構成することもできる。
また、プログラムは、記録媒体を介さず、通信回線を通じて直接にコンピュータにロードし実行することもできる。また、本発明に係る含水比マッピング装置に備えられる、単一の情報処理装置（例えば１台のパーソナルコンピュータ等）で構成することもでき、複数の情報処理装置（例えば複数台のコンピュータ群等）で構成することもできる。
【0017】
［含水比マッピング装置］
図１に、本発明の一実施形態に係る含水比マッピング装置１の構成を示す。
同図に示す本発明の一実施形態に係る含水比マッピング装置１は、探査車両１０と、探査車両１０に搭載される地中レーダ２０及び全方位カメラ３０を備える構成となっている。
そして、図１に示すように、探査車両１０が探査対象となる任意の路面を走行しながら、地中レーダ２０によって路面下の地中探査データを取得・生成するとともに、全方位カメラ３０によって路面上の地上三次元映像を取得・生成し、その後、地中探査データと地上三次元映像とを、一元化された単一の一元化情報（後述する図１４参照）として生成・出力するようになっている。
以下、含水比マッピング装置１の各構成を具体的に説明する。
図１に、本発明の一実施形態に係る含水比マッピング装置１の構成を示す。
同図に示す本発明の一実施形態に係る含水比マッピング装置１は、探査車両１０と、探査車両１０に搭載される地中レーダ２０及び全方位カメラ３０を備える構成となっている。
そして、図１に示すように、探査車両１０が探査対象となる任意の路面を走行しながら、地中レーダ２０によって路面下の地中探査データを取得・生成するとともに、全方位カメラ３０によって路面上の地上三次元映像を取得・生成し、その後、地中探査データと地上三次元映像とを、一元化された単一の一元化情報（後述する図１４参照）として生成・出力するようになっている。
以下、含水比マッピング装置１の各構成を具体的に説明する。
【0018】
［探査車両１０］
探査車両１０は、地中探査の対象・範囲となる道路や地面などの路面を走行可能な車両であり、例えば乗用車やトラック，ライトバン，各種作業車，軽自動車などの自動車で構成される。なお、探査対象・目的となる所望の路面を走行可能である限り、探査車両１０の具体的な構成は特に限定されるものではない。
図２に、本実施形態の探査車両１０の動作イメージを模式的に示す。
同図（ａ）に示すように、探査車両１０には、屋根上面に車載装備としてのＧＰＳ（高精度ＧＮＳＳ）１１が取り付けられ、車両底面に地中レーダ２０のアンテナ２２が配置され、また、屋根後方に全方位カメラ３０が配置されている。
また、探査車両１０の車内には、地中レーダ２０のレーダ本体２１が配置されるとともに、ＣＶ演算部４０，一元化処理部５０が備えられる。なお、レーダ本体２１，ＣＶ演算部４０及び一元化処理部５０は、具体的には所定のソフトウェア（プログラム）が実装されたＰＣなどの情報処理装置で構成される。
探査車両１０は、地中探査の対象・範囲となる道路や地面などの路面を走行可能な車両であり、例えば乗用車やトラック，ライトバン，各種作業車，軽自動車などの自動車で構成される。なお、探査対象・目的となる所望の路面を走行可能である限り、探査車両１０の具体的な構成は特に限定されるものではない。
図２に、本実施形態の探査車両１０の動作イメージを模式的に示す。
同図（ａ）に示すように、探査車両１０には、屋根上面に車載装備としてのＧＰＳ（高精度ＧＮＳＳ）１１が取り付けられ、車両底面に地中レーダ２０のアンテナ２２が配置され、また、屋根後方に全方位カメラ３０が配置されている。
また、探査車両１０の車内には、地中レーダ２０のレーダ本体２１が配置されるとともに、ＣＶ演算部４０，一元化処理部５０が備えられる。なお、レーダ本体２１，ＣＶ演算部４０及び一元化処理部５０は、具体的には所定のソフトウェア（プログラム）が実装されたＰＣなどの情報処理装置で構成される。
【0019】
［地中レーダ２０］
探査車両１０に備えられる地中レーダ２０は、探査車両１０が走行する路面下の地中探査データを取得・生成するための地中探査手段（センサ）であり、図２に示すように、探査車両１０の車内に搭載されるレーダ本体２１と、探査車両１０の底面に配置されるアンテナ２２とで構成される。
本実施形態に係る地中レーダ２０は、地中の物理的境界面で電磁波が反射する現象を利用して地中の土壌の含水比を測定し、また、地中の空洞や埋設管等の存在を探査・検出する物理探査レーダである。
このような地中レーダ２０を備えることによって、本実施形態では、レーダ本体２１の制御によってアンテナ２２から電磁波を送受信させながら探査車両１０を所定速度（例えば最大時速６０ｋｍ）で走行させることで、地中の所定深度（例えば地下２ｍ）までの含水比の分布及び含水比境界面を示す三次元情報を、簡単にマッピングすることが可能となっている（後述する図４参照）。
探査車両１０に備えられる地中レーダ２０は、探査車両１０が走行する路面下の地中探査データを取得・生成するための地中探査手段（センサ）であり、図２に示すように、探査車両１０の車内に搭載されるレーダ本体２１と、探査車両１０の底面に配置されるアンテナ２２とで構成される。
本実施形態に係る地中レーダ２０は、地中の物理的境界面で電磁波が反射する現象を利用して地中の土壌の含水比を測定し、また、地中の空洞や埋設管等の存在を探査・検出する物理探査レーダである。
このような地中レーダ２０を備えることによって、本実施形態では、レーダ本体２１の制御によってアンテナ２２から電磁波を送受信させながら探査車両１０を所定速度（例えば最大時速６０ｋｍ）で走行させることで、地中の所定深度（例えば地下２ｍ）までの含水比の分布及び含水比境界面を示す三次元情報を、簡単にマッピングすることが可能となっている（後述する図４参照）。
【0020】
具体的には、地中レーダ２０は、電磁波を利用した走る地下ＣＴスキャナーとして、ノルウェーの３ｄ－Ｒａｄａｒ社製の最大４１チャンネルの地中レーダアンテナ（マルチチャンネル）を使用し、探査車両１０を走行させながら地中の含水比分布を三次元で測定し可視化（マッピング）することができるようになっている。
ここで、従来技術の手押し型のセンサ等で使用されていた地中レーダは１チャンネル（シングルチャンネル）であった。シングルチャンネルは、電磁波の発信と受信のセンサが一対となっているアンテナ方式である。
ここで、従来技術の手押し型のセンサ等で使用されていた地中レーダは１チャンネル（シングルチャンネル）であった。シングルチャンネルは、電磁波の発信と受信のセンサが一対となっているアンテナ方式である。
【0021】
これに対して、本実施形態の地中レーダ２０で採用するマルチチャンネルは、アンテナ内に複数の送信アンテナ（送信部）と受信アンテナ（受信部）とが備えられており、送信アンテナと受信アンテナの組み合わせを任意の複数組だけ設定できるようになっている。
本実施形態に係る地中レーダ２０は、マルチチャンネルアンテナが、例えば７．５ｃｍピッチで高密度に配置されており、地中の含水比の分布を高精細な三次元情報として取得できるようになっている。
本実施形態に係る地中レーダ２０は、マルチチャンネルアンテナが、例えば７．５ｃｍピッチで高密度に配置されており、地中の含水比の分布を高精細な三次元情報として取得できるようになっている。
【0022】
このようなマルチチャンネル型の地中レーダ２０を用いることによって、送信アンテナと受信アンテナが同一の中点と異なるオフセット距離を有する複数の送受信アンテナ組を設定することができ、ＣＭＰ（Ｃｏｍｍｏｎ Ｍｉｄｐｏｉｎｔ：コモンミッドポイント）法によるアンテナ配置を、単一の地中レーダ２０によって構成・機能させることができるようになる（後述する図６参照）。
また、マルチチャンネル型の地中レーダ２０は、送信アンテナと受信アンテナのオフセット距離を一定（同一）にした複数の送受信アンテナ組を設定することができ、ＣＯＧ（Ｃｏｍｍｏｎ Ｏｆｆｓｅｔ Ｇａｔｈｅｒ）法によるアンテナ配置も設定することができる。これによって、ＣＭＰ法による測定結果を、ＣＯＧ法による測定と組み合わせることができ、より精度の高い地中探査・含水比測定が可能となる。
また、マルチチャンネル型の地中レーダ２０は、送信アンテナと受信アンテナのオフセット距離を一定（同一）にした複数の送受信アンテナ組を設定することができ、ＣＯＧ（Ｃｏｍｍｏｎ Ｏｆｆｓｅｔ Ｇａｔｈｅｒ）法によるアンテナ配置も設定することができる。これによって、ＣＭＰ法による測定結果を、ＣＯＧ法による測定と組み合わせることができ、より精度の高い地中探査・含水比測定が可能となる。
【0023】
そして、本実施形態では、上記のような地中レーダ２０のアンテナ２２が、図２（ａ）に示すように、探査車両１０の底面の前輪と後輪の間、あるいは複数の車輪の車軸間に配置されるようになっている。
このように地中レーダ２０のアンテナ２２を探査車両１０の底面に配置することにより、図２（ｂ）に示すように、探査対象となる路面を複数回（例えば同じ道路を２回）、車幅方向にずらして走行させることで、アンテナ２２の電磁波の照射範囲（照射幅）を超える道路であっても、路面の全幅にわたる地中探査データを取得することができる。図２（ｂ）では、探査車両１０を同一路面で２回走行させることにより、７．５ｃｍ間隔の地中探査データを、例えば道路幅３．０ｍの全範囲で取得できる場合を示している。
このように地中レーダ２０のアンテナ２２を探査車両１０の底面に配置することにより、図２（ｂ）に示すように、探査対象となる路面を複数回（例えば同じ道路を２回）、車幅方向にずらして走行させることで、アンテナ２２の電磁波の照射範囲（照射幅）を超える道路であっても、路面の全幅にわたる地中探査データを取得することができる。図２（ｂ）では、探査車両１０を同一路面で２回走行させることにより、７．５ｃｍ間隔の地中探査データを、例えば道路幅３．０ｍの全範囲で取得できる場合を示している。
【0024】
ここで、地表面にアンテナを配置する形態は、「グランド型」と呼ばれ、従来のアンテナが車両の前方や後方に大きく突き出たり、アンテナを搭載した荷台等をけん引する形態の「エア型」と比較して多くの利点を有する。
エア型は、未舗装地帯の地雷撤去での使用を想定しており、アンテナを地表より約２０～３０ｃｍ浮かせて測定を行うようになっている。このため、車両の前方に大きく突出するように取り付けたり、車両でけん引するなど、車両への装着等は簡単にできるが、アンテナと地表の距離が大きいことから、電磁波が減衰して探査深度が浅くなるという欠点がある。
このように、従来のエア型は、地表から２０～３０ｃｍ浮かせる必要があるアンテナ特性と、アンテナ自体の高さも約２０ｃｍと厚さがあるため、アンテナを直接車両に配置する車載型（グランド型）としては使用できなかった。
エア型は、未舗装地帯の地雷撤去での使用を想定しており、アンテナを地表より約２０～３０ｃｍ浮かせて測定を行うようになっている。このため、車両の前方に大きく突出するように取り付けたり、車両でけん引するなど、車両への装着等は簡単にできるが、アンテナと地表の距離が大きいことから、電磁波が減衰して探査深度が浅くなるという欠点がある。
このように、従来のエア型は、地表から２０～３０ｃｍ浮かせる必要があるアンテナ特性と、アンテナ自体の高さも約２０ｃｍと厚さがあるため、アンテナを直接車両に配置する車載型（グランド型）としては使用できなかった。
【0025】
これに対して、グランド型は、レーダのアンテナを路面から近距離での使用を前提とするもので、具体的には地表より約３～１０ｃｍで測定することができる。
このため、電磁波が減衰しにくく、探査深度はエア型と比較して約１．５～２倍程度深くなる。
また、グランド型は、可能な限りアンテナを地表に近接させることができ、また、車両に直接配置するために、アンテナ高さも１０ｃｍ程度と薄いため、地中探査における測定の安定性及び安全性に優れる。
このため、電磁波が減衰しにくく、探査深度はエア型と比較して約１．５～２倍程度深くなる。
また、グランド型は、可能な限りアンテナを地表に近接させることができ、また、車両に直接配置するために、アンテナ高さも１０ｃｍ程度と薄いため、地中探査における測定の安定性及び安全性に優れる。
【0026】
本実施形態では、上記のようなメリットの大きいグランド型を採用し、地中探査手段となる地中レーダ２０のアンテナ２２を探査車両１０の底面に配置している。
具体的には、上述したノルウェーの３ｄ－Ｒａｄａｒ社製のグランド型アンテナを採用し、地中レーダ２０のアンテナ２２を、探査車両１０の前後輪の車軸間（ホイルベース間）に取り付けるようにしている。
３ｄ－Ｒａｄａｒ社のグランド型アンテナは、従来製品のエア型アンテナや他社のグランド型アンテナに比べて、アンテナ高さが半分以上薄く、自動車等の前後輪車軸間の最も安定した位置に設置することが可能となる。
具体的には、上述したノルウェーの３ｄ－Ｒａｄａｒ社製のグランド型アンテナを採用し、地中レーダ２０のアンテナ２２を、探査車両１０の前後輪の車軸間（ホイルベース間）に取り付けるようにしている。
３ｄ－Ｒａｄａｒ社のグランド型アンテナは、従来製品のエア型アンテナや他社のグランド型アンテナに比べて、アンテナ高さが半分以上薄く、自動車等の前後輪車軸間の最も安定した位置に設置することが可能となる。
【0027】
ここで、地中レーダ２０のアンテナ２２を探査車両１０のホイルベース間に設置することで、以下のようなメリットがある。
・探査車両１０の車両全長を短くすることにより内輪差を小さくできる。
・探査車両１０の車両全長を短くすることにより小回りがきく。
・けん引型のように後方を気にする必要がない。
・けん引型に比べて後進が容易になる。
・後輪後方設置型よりアンテナが路面に対して一定の高さになり、安定したデータが取得できる。
・探査車両１０の車両全長を短くすることにより内輪差を小さくできる。
・探査車両１０の車両全長を短くすることにより小回りがきく。
・けん引型のように後方を気にする必要がない。
・けん引型に比べて後進が容易になる。
・後輪後方設置型よりアンテナが路面に対して一定の高さになり、安定したデータが取得できる。
【0028】
地中レーダ探査にとって、アンテナの高さが一定であることはデータの精度上きわめて重要であり、地中レーダ２０のアンテナ２２を探査車両１０のホイルベース間に配置することのメリットは非常に大きなものとなる。
なお、探査車両１０のホイルベース間に配置されるアンテナ２２の位置は、探査車両１０の前輪と後輪の間、あるいは複数の車輪の間において、車長方向の最適な位置に配置される。すなわち、探査車両１０の車長や前後輪の配置，レーダの出力や特性などに応じて、探査車両１０の車長方向（車両進行方向の前後）の適切・最適な位置にアンテナ２２を配置する。
したがって、アンテナ２２は、探査車両１０の車長方向に移動可能・調整可能に構成されることが好ましい。
なお、探査車両１０のホイルベース間に配置されるアンテナ２２の位置は、探査車両１０の前輪と後輪の間、あるいは複数の車輪の間において、車長方向の最適な位置に配置される。すなわち、探査車両１０の車長や前後輪の配置，レーダの出力や特性などに応じて、探査車両１０の車長方向（車両進行方向の前後）の適切・最適な位置にアンテナ２２を配置する。
したがって、アンテナ２２は、探査車両１０の車長方向に移動可能・調整可能に構成されることが好ましい。
【0029】
また、本実施形態では、上記のような地中レーダ２０のアンテナ２２を探査車両１０のホイルベース間において、車幅方向（車両進行方向の左右）に移動可能・調整可能に構成してある。
具体的には、図３（ａ），（ｂ）に示すように、地中レーダ２０のアンテナ２２を、探査車両１０の車幅方向の左右にスライド移動できるようになっており、車幅方向の所望の位置で固定可能となっている。これによって、探査対象となる路面の大きさ（道路幅）に応じて、アンテナ２２を移動・調整することができ、様々な大きさ（幅）の路面に対応することができる。
具体的には、図３（ａ），（ｂ）に示すように、地中レーダ２０のアンテナ２２を、探査車両１０の車幅方向の左右にスライド移動できるようになっており、車幅方向の所望の位置で固定可能となっている。これによって、探査対象となる路面の大きさ（道路幅）に応じて、アンテナ２２を移動・調整することができ、様々な大きさ（幅）の路面に対応することができる。
【0030】
図３に示す例では、探査幅１．８ｍのアンテナ２２を、探査車両１０の左右側面からそれぞれ１５ｃｍ突出させるように移動させることで、道路幅３．５ｍの路面について、ｂ）では、２回の走行により地中探査データを取得することができる場合を示している。
このように、地中レーダ２０のアンテナ２２を探査車両１０のホイルベース間において車幅方向に移動可能とすることで、アンテナ幅（例えば１．８ｍ）のアンテナ１台を探査車両の左右寄りにスライドさせて複数回（例えば２回）走行することで、アンテナ幅を超える範囲（例えば最大３．５ｍ）の幅でデータを取得することができる（図３参照）。
また、図３に示すように、探査車両１０の側面左右への飛び出しは、片側１５ｃｍ程度とすることで、例えば従来のエア型アンテナで使用されている２．５ｍ幅のアンテナに比べて、路側の障害物や第三者等との接触事故の危険性を抑制乃至回避することできる。
このように、地中レーダ２０のアンテナ２２を探査車両１０のホイルベース間において車幅方向に移動可能とすることで、アンテナ幅（例えば１．８ｍ）のアンテナ１台を探査車両の左右寄りにスライドさせて複数回（例えば２回）走行することで、アンテナ幅を超える範囲（例えば最大３．５ｍ）の幅でデータを取得することができる（図３参照）。
また、図３に示すように、探査車両１０の側面左右への飛び出しは、片側１５ｃｍ程度とすることで、例えば従来のエア型アンテナで使用されている２．５ｍ幅のアンテナに比べて、路側の障害物や第三者等との接触事故の危険性を抑制乃至回避することできる。
【0031】
なお、以上のようなアンテナ２２の車幅方向への移動可能な範囲は、アンテナ幅や探査車両１０の車幅等に応じて任意に設定することができる。
また、アンテナ２２の移動構造も、スライド構造や、ボルト等による固定位置を複数設ける段階構造など、アンテナ２２を車幅方向に移動できる限り、特に限定されるものではない。
また、アンテナ２２は、探査車両１０の底面に露出した剥き出しの状態で配置することもできるが、異物の跳ね返りや衝突、地面や障害物との接触等による破損などに備えて、アンテナ２２の全部又は一部をカバー等で覆うこともできる。例えばアンテナ２２の電磁波の送受信に影響のないプラスチック製のカバーを備えることができる。
また、アンテナ２２の移動構造も、スライド構造や、ボルト等による固定位置を複数設ける段階構造など、アンテナ２２を車幅方向に移動できる限り、特に限定されるものではない。
また、アンテナ２２は、探査車両１０の底面に露出した剥き出しの状態で配置することもできるが、異物の跳ね返りや衝突、地面や障害物との接触等による破損などに備えて、アンテナ２２の全部又は一部をカバー等で覆うこともできる。例えばアンテナ２２の電磁波の送受信に影響のないプラスチック製のカバーを備えることができる。
【0032】
［含水比測定／境界面探査］
以上のような本実施形態に係る地中レーダ２０で取得・生成される地中探査データによる含水比測定と含水比の境界面探査の概要について、図４を参照して説明する。
本実施形態の地中レーダ２０によれば、まず、図４（ａ）に示すように、探査車両１０が道路などの路面に沿って走行しつつ地中レーダ２０によって地中探査データを取得することにより、走行方向と地中深度方向（縦方向）の二次元に広がる含水比の測定とそれに基づく含水比境界線の探査が可能となる。
以上のような本実施形態に係る地中レーダ２０で取得・生成される地中探査データによる含水比測定と含水比の境界面探査の概要について、図４を参照して説明する。
本実施形態の地中レーダ２０によれば、まず、図４（ａ）に示すように、探査車両１０が道路などの路面に沿って走行しつつ地中レーダ２０によって地中探査データを取得することにより、走行方向と地中深度方向（縦方向）の二次元に広がる含水比の測定とそれに基づく含水比境界線の探査が可能となる。
【0033】
また、本実施形態の地中レーダ２０は、図４（ｂ）に示すように、探査車両の車幅方向に複数の送受信アンテナ組が配置されるマルチチャンネルレーダとなっており、地中深度方向と車幅方向（横方向）の二次元に広がる含水比の測定とそれに基づく含水比境界線の探査が可能となる。
そして、このような走行方向と地中深度方向（図４（ａ））、地中深度方向と車幅方向（図４（ｂ））の二つの測定データを合成することにより、図４（ｃ）に示すように、三次元に広がる含水比の分布と含水比境界面（マップ）の探査・生成が可能となる。
そして、このような走行方向と地中深度方向（図４（ａ））、地中深度方向と車幅方向（図４（ｂ））の二つの測定データを合成することにより、図４（ｃ）に示すように、三次元に広がる含水比の分布と含水比境界面（マップ）の探査・生成が可能となる。
【0034】
このように、探査車両１０の走行に伴って移動する地中レーダ２０では、マルチチャンネルのアンテナ２２によって、所定間隔（例えば７．５ｃｍピッチ）で地中の二次元情報が車両進行方向に沿って取得され、含水比の分布及び境界面を示す地中三次元情報として取得・生成される。
また、このような含水比の分布／境界面を示す三次元情報は、画像生成することで、地中の三次元画像情報として出力することができる。
なお、以上のような地中探査データに基づく含水比測定情報や含水比の分布／境界面情報の生成・出力は、地中レーダ２０のレーダ本体２１や、レーダ本体２１に接続されるＰＣ等によって行われる。
また、このような含水比の分布／境界面を示す三次元情報は、画像生成することで、地中の三次元画像情報として出力することができる。
なお、以上のような地中探査データに基づく含水比測定情報や含水比の分布／境界面情報の生成・出力は、地中レーダ２０のレーダ本体２１や、レーダ本体２１に接続されるＰＣ等によって行われる。
【0035】
【0036】
［含水比測定原理］
まず、地中の土壌の含水比を測定する原理について説明する。
土壌中の電磁場伝搬速度ｖは、以下の式１となる。
式１中、土壌の比誘電率εr、真空中の電磁場伝搬速度ｃとする。
［式１］
上記式１を変形すると、以下の式２となる。
この式２より、電磁波速度ｖから比誘電率εrを計算することができる。
［式２］
まず、地中の土壌の含水比を測定する原理について説明する。
土壌中の電磁場伝搬速度ｖは、以下の式１となる。
式１中、土壌の比誘電率εr、真空中の電磁場伝搬速度ｃとする。
［式１］
この式２より、電磁波速度ｖから比誘電率εrを計算することができる。
［式２］
【0037】
上述のとおり、本実施形態の地中レーダ２０は、送信アンテナと受信アンテナが同一の中点と異なるオフセット距離を有する複数の送受信アンテナ組として設定することができる。これにより、複数の送受信アンテナ組の測定は異なる経路となる。
そして、各送受信アンテナ組の異なるオフセット距離ｘについて、地中レーダ２０の反射波の到達時間τを測定することにより（後述する図７参照）、以下の式３に示すように、地表面から地中の物理的境界面（埋設物や含水比境界面等）までの深さｄと比誘電率εr（上記式２）を求めることができる。
［式３］
そして、各送受信アンテナ組の異なるオフセット距離ｘについて、地中レーダ２０の反射波の到達時間τを測定することにより（後述する図７参照）、以下の式３に示すように、地表面から地中の物理的境界面（埋設物や含水比境界面等）までの深さｄと比誘電率εr（上記式２）を求めることができる。
［式３］
【0038】
上述した特許文献１で提案されている測定方法では、埋設物の深さが既知であることを前提としており、埋設物や含水比境界面の深さが不明な土壌については、比誘電率を測定・算出することはできなかった。
本実施形態の含水比測定方法によれば、深度が未知の土壌について、含水比を求めることが可能となる。
以上のような含水比の測定原理に基づいて、以下に具体的な含水比マッピング方法の各工程について説明する。
本実施形態の含水比測定方法によれば、深度が未知の土壌について、含水比を求めることが可能となる。
以上のような含水比の測定原理に基づいて、以下に具体的な含水比マッピング方法の各工程について説明する。
【0039】
［含水比マッピング方法］
図５は、本発明の一実施形態に係る含水比マッピング方法の各工程（ステップ１～ステップ７）を示すフローチャートである。
同図に示す含水比マッピング方法は、地中レーダ２０において設定・構成される送信アンテナと受信アンテナが同一の中点と異なるオフセット距離を有する複数の送受信アンテナ組を用いて地中の含水比を測定し、含水比分布のマッピングを行う方法である。
なお、以下に示す含水比マッピング方法の各工程は、地中レーダ２０のレーダ本体２１や、レーダ本体２１に接続されるＰＣ等によって処理・実行される。
図５は、本発明の一実施形態に係る含水比マッピング方法の各工程（ステップ１～ステップ７）を示すフローチャートである。
同図に示す含水比マッピング方法は、地中レーダ２０において設定・構成される送信アンテナと受信アンテナが同一の中点と異なるオフセット距離を有する複数の送受信アンテナ組を用いて地中の含水比を測定し、含水比分布のマッピングを行う方法である。
なお、以下に示す含水比マッピング方法の各工程は、地中レーダ２０のレーダ本体２１や、レーダ本体２１に接続されるＰＣ等によって処理・実行される。
【0040】
(ステップ１)
まず、図５に示すステップ１により、地中レーダ２０の反射波が取得される。
具体的には、地中レーダ２０を探査車両１０の移動に伴って地表面に沿って移動させながら、アンテナ２２から地中へ向けて電磁波を放射することにより、複数の各送受信アンテナ組によって得られる複数の反射波を取得する（本発明に係るレーダ反射波取得部）。
地中レーダ２０の送信アンテナが、電磁波を地中へ向けて放射すると、電気的性質が地層境界で不連続に変わるため、放射された電磁波は地層境界で反射をする。
地中レーダ２０の受信アンテナでは、地中で反射された電磁波（反射波）が受信されるとともに、送信アンテナからの電磁波が直接波として受信される。
まず、図５に示すステップ１により、地中レーダ２０の反射波が取得される。
具体的には、地中レーダ２０を探査車両１０の移動に伴って地表面に沿って移動させながら、アンテナ２２から地中へ向けて電磁波を放射することにより、複数の各送受信アンテナ組によって得られる複数の反射波を取得する（本発明に係るレーダ反射波取得部）。
地中レーダ２０の送信アンテナが、電磁波を地中へ向けて放射すると、電気的性質が地層境界で不連続に変わるため、放射された電磁波は地層境界で反射をする。
地中レーダ２０の受信アンテナでは、地中で反射された電磁波（反射波）が受信されるとともに、送信アンテナからの電磁波が直接波として受信される。
【0041】
［ＣＭＰ法による送受信アンテナのオフセット配置］
ここで、地中レーダ２０の複数の送受信アンテナ組におけるオフセット距離の設定について、図６を参照しつつ説明する。
図６は、本発明の一実施形態に係る含水比マッピング装置を構成する送受信アンテナの中点・オフセット距離の関係を模式的に示す図であり、（ａ）は送受信アンテナの断面正面図、（ｂ）は中点と複数の送受信アンテナ組の位置を示す説明図である。
本実施形態の含水比マッピング方法では、ＣＭＰ法を用いている。
ＣＭＰ法では、図６（ａ）に示すように、地上の接地面（実際には探査車両１０の底面）に、送信アンテナと受信アンテナが、所定のオフセット距離ｘだけ離間するようにして配置される。同図は、横断方向からの模式図を示しており、地中の層１と層２は異なる比誘電率を持つ一様な層である。
ここで、地中レーダ２０の複数の送受信アンテナ組におけるオフセット距離の設定について、図６を参照しつつ説明する。
図６は、本発明の一実施形態に係る含水比マッピング装置を構成する送受信アンテナの中点・オフセット距離の関係を模式的に示す図であり、（ａ）は送受信アンテナの断面正面図、（ｂ）は中点と複数の送受信アンテナ組の位置を示す説明図である。
本実施形態の含水比マッピング方法では、ＣＭＰ法を用いている。
ＣＭＰ法では、図６（ａ）に示すように、地上の接地面（実際には探査車両１０の底面）に、送信アンテナと受信アンテナが、所定のオフセット距離ｘだけ離間するようにして配置される。同図は、横断方向からの模式図を示しており、地中の層１と層２は異なる比誘電率を持つ一様な層である。
【0042】
このようなＣＭＰ法による送受信アンテナの配置において、送信アンテナから地中へ向かって電磁波が放射されると、地中の比誘電率が不連続な境界面において電磁波が反射され、受信アンテナによって地表面側に反射された電磁波が検出される。これによって、図６（ａ）に示すように、地表から境界面までの深度ｄが算出できるようになる。
図６（ａ）の縦線丸印と横線丸印は、上面（地表面）からみた送信アンテナと受信アンテナの配置位置であり、白抜丸印は送受信アンテナの中点の位置を表現したものである。
すなわち、図６（ａ）において、左から「送信アンテナ・中点・受信アンテナ」を表しており、送信アンテナと受信アンテナの距離（間隔）が「オフセット距離ｘ」として定義されている。
図６（ａ）の縦線丸印と横線丸印は、上面（地表面）からみた送信アンテナと受信アンテナの配置位置であり、白抜丸印は送受信アンテナの中点の位置を表現したものである。
すなわち、図６（ａ）において、左から「送信アンテナ・中点・受信アンテナ」を表しており、送信アンテナと受信アンテナの距離（間隔）が「オフセット距離ｘ」として定義されている。
【0043】
また、図６（ｂ）は、中点位置Ｃを固定し、３組の異なるオフセット距離のアンテナ配置を表している。このように、送信アンテナと受信アンテナが同一の中点と異なるオフセット距離を有する複数の送受信アンテナ組によって行われる測定方法が、ＣＭＰ法と呼ばれるものである。
図６（ｂ）の例では、送受信アンテナの組が、「Ａ１・Ｂ１」の組，「Ａ２・Ｂ２」の組，「Ａ３・Ｂ３」の組の、異なる３つのオフセット距離で測定が行われる場合を示している。
なお、異なるオフセット距離の設定数は、図６に示す３つの場合に限定されず、地中レーダ２０の構造（アンテナ数・チャンネル数）や大きさ，規模などに応じて、適宜任意の数を設定することができ、オフセット距離の設定数が多くなれば、より高精度な地中探査・含水比測定が行えるようになる。
図６（ｂ）の例では、送受信アンテナの組が、「Ａ１・Ｂ１」の組，「Ａ２・Ｂ２」の組，「Ａ３・Ｂ３」の組の、異なる３つのオフセット距離で測定が行われる場合を示している。
なお、異なるオフセット距離の設定数は、図６に示す３つの場合に限定されず、地中レーダ２０の構造（アンテナ数・チャンネル数）や大きさ，規模などに応じて、適宜任意の数を設定することができ、オフセット距離の設定数が多くなれば、より高精度な地中探査・含水比測定が行えるようになる。
【0044】
(ステップ２)
次に、図５に示すステップ２により、直接波フィッティング処理が行われる。
具体的には、上記ステップ１で取得された反射波（電磁波）について、地中レーダ２０の複数の各送受信アンテナ組において、送信アンテナから受信アンテナへの到達時間が最短の電磁波（送信アンテナから受信アンテナに最初に到達する電磁波）を直接波として、その直接波の到達時刻が抽出される（本発明に係る直接波フィッティング部）。
探査車両１０の移動に伴い、地中レーダ２０では、中点を同一とした送信アンテナと受信アンテナの複数組によって、複数の測定が行われる。
各測定において、受信アンテナへの到達時間が一番早い振幅の大きな波を直接波とみなして（図７参照）、直接波の到達時間の抽出が行われる。この直接波の到達時刻の抽出が、複数の各送信アンテナと受信アンテナの測定に対してそれぞれ行われる。
なお、上記手法以外にも、直接波の到達時間は、例えば、ステップ周波数連続波による周波数毎の強度と位相の測定から直接波の到達時間を求めることもできる。
次に、図５に示すステップ２により、直接波フィッティング処理が行われる。
具体的には、上記ステップ１で取得された反射波（電磁波）について、地中レーダ２０の複数の各送受信アンテナ組において、送信アンテナから受信アンテナへの到達時間が最短の電磁波（送信アンテナから受信アンテナに最初に到達する電磁波）を直接波として、その直接波の到達時刻が抽出される（本発明に係る直接波フィッティング部）。
探査車両１０の移動に伴い、地中レーダ２０では、中点を同一とした送信アンテナと受信アンテナの複数組によって、複数の測定が行われる。
各測定において、受信アンテナへの到達時間が一番早い振幅の大きな波を直接波とみなして（図７参照）、直接波の到達時間の抽出が行われる。この直接波の到達時刻の抽出が、複数の各送信アンテナと受信アンテナの測定に対してそれぞれ行われる。
なお、上記手法以外にも、直接波の到達時間は、例えば、ステップ周波数連続波による周波数毎の強度と位相の測定から直接波の到達時間を求めることもできる。
【0045】
(ステップ３)
次に、図５に示すステップ３により、ゼロ時間補正処理が行われる。
具体的には、上記ステップ２で求められた直接波の到達時刻から、当該直接波の地表面における到達時刻を時刻ゼロとする算出・補正が行われる（本発明に係るゼロ時間補正部）。
地中レーダ２０は、探査車両１０の底面に配置されており（図２（ａ）参照）、地表面から例えば約３～１０ｃｍ上方に位置している。
このため、直接波の到達時刻を、地表面の位置における到達時刻が時刻ゼロとなるように校正・補正を行う。
次に、図５に示すステップ３により、ゼロ時間補正処理が行われる。
具体的には、上記ステップ２で求められた直接波の到達時刻から、当該直接波の地表面における到達時刻を時刻ゼロとする算出・補正が行われる（本発明に係るゼロ時間補正部）。
地中レーダ２０は、探査車両１０の底面に配置されており（図２（ａ）参照）、地表面から例えば約３～１０ｃｍ上方に位置している。
このため、直接波の到達時刻を、地表面の位置における到達時刻が時刻ゼロとなるように校正・補正を行う。
【0046】
ここで、直接波の到達時刻のゼロ補正について具体的に説明する。
図７は、オフセット距離が異なる複数の送受信アンテナ組によって取得される直接波・反射波の到達時間とオフセット距離の関係を示すグラフである。
図７のグラフでは、横軸がオフセット距離、縦軸が時間を表している。
同図では、縦軸方向の複数の各点線が、各オフセット距離における測定と対応しており、それぞれのオフセット距離の振幅強度の時間変化を示している。
なお、図７では、説明の煩雑さなくして理解を容易にするため、直接波と反射波のみを示し、多重反射成分などは割愛している。
図７は、オフセット距離が異なる複数の送受信アンテナ組によって取得される直接波・反射波の到達時間とオフセット距離の関係を示すグラフである。
図７のグラフでは、横軸がオフセット距離、縦軸が時間を表している。
同図では、縦軸方向の複数の各点線が、各オフセット距離における測定と対応しており、それぞれのオフセット距離の振幅強度の時間変化を示している。
なお、図７では、説明の煩雑さなくして理解を容易にするため、直接波と反射波のみを示し、多重反射成分などは割愛している。
【0047】
各測定で最初に現れる信号は、送信アンテナから受信アンテナへ直線で空気中を伝搬する直接波である。
直接波の到達時間は、以下の式４で与えられる。各オフセット距離ｘ１，ｘ２・・・のときに、対応する時刻ｔ１，ｔ２・・・が測定結果として得られる。
受信アンテナで受信された信号から、直接波の到達時間を抽出する。抽出された時刻から、地表面の時刻をゼロとする時間を「時刻ｔ０」とする校正を行う。
［式４］
直接波の到達時間は、以下の式４で与えられる。各オフセット距離ｘ１，ｘ２・・・のときに、対応する時刻ｔ１，ｔ２・・・が測定結果として得られる。
受信アンテナで受信された信号から、直接波の到達時間を抽出する。抽出された時刻から、地表面の時刻をゼロとする時間を「時刻ｔ０」とする校正を行う。
［式４］
【0048】
(ステップ４)
次に、ステップ４により、反射波フィッティング処理が行われる。
具体的には、地中レーダ２０の複数の各送受信アンテナ組において、送信アンテナから受信アンテナへ直接波の次に到達する電磁波を反射波として、その反射波の到達時刻が抽出される（本発明に係る反射波フィッティング部）。
直接波の次に来る反射は、地層境界からの反射と考えることができ、その到達時刻を反射波の到達時刻として抽出する。
次に、ステップ４により、反射波フィッティング処理が行われる。
具体的には、地中レーダ２０の複数の各送受信アンテナ組において、送信アンテナから受信アンテナへ直接波の次に到達する電磁波を反射波として、その反射波の到達時刻が抽出される（本発明に係る反射波フィッティング部）。
直接波の次に来る反射は、地層境界からの反射と考えることができ、その到達時刻を反射波の到達時刻として抽出する。
【0049】
反射波の到達時刻は、境界面が平面である場合には、以下の式５となる。
［式５］
この反射波の到達時刻の抽出が、複数の各送信アンテナと受信アンテナの測定に対してそれぞれ行われる。
［式５］
【0050】
また、含水比測定の精度を高めるために、複数の送受信アンテナのオフセット距離を同一としたＣＯＧ法による測定を行い、その結果から層境界の反射波が、含水比境界面によるものか、埋設管などの埋設物によるものかを判断することも可能である。
例えば、到達時刻が抽出された反射波の近傍について、ＣＯＧ法で測定されたデータを使用することにより、対応する中点で得られた反射波が、埋設物による反射やノイズではないことをＣＯＧ法で得られた反射映像から判断することができる。
例えば、到達時刻が抽出された反射波の近傍について、ＣＯＧ法で測定されたデータを使用することにより、対応する中点で得られた反射波が、埋設物による反射やノイズではないことをＣＯＧ法で得られた反射映像から判断することができる。
【0051】
(ステップ５)
次に、ステップ５により、比誘電率・深度計算処理が行われる。
具体的には、上記ステップ４で求められた反射波の到達時刻に基づいて、地中の比誘電率が算出される（本発明に係る比誘電率算出部）。
上述した式１～式３に示したように、反射波到達時間を、比誘電率εｒと深度ｄ（図６（a）参照）の２つのパラメータとして持つ双曲線モデルで回帰解析を行うことにより、比誘電率と深度を決定することができる。
次に、ステップ５により、比誘電率・深度計算処理が行われる。
具体的には、上記ステップ４で求められた反射波の到達時刻に基づいて、地中の比誘電率が算出される（本発明に係る比誘電率算出部）。
上述した式１～式３に示したように、反射波到達時間を、比誘電率εｒと深度ｄ（図６（a）参照）の２つのパラメータとして持つ双曲線モデルで回帰解析を行うことにより、比誘電率と深度を決定することができる。
【0052】
ここで、比誘電率の算出について具体的に説明する。
上述のとおり、反射波の到達時刻は、境界面が平面である場合には、上記式５のようになる。
そして、式５により抽出された反射波の到達時刻から、層１（図６（ａ）参照）の比誘電率εｒと深度ｄを、最小自乗誤差の最適解をLevenberg-Marquardt法により計算することができる。
得られた比誘電率εｒは、Ｔｏｐｐらによる比誘電率と含水比の関係式である以下の式６から、層１の含水比θｗを計算することができる。
［式６］
上述のとおり、反射波の到達時刻は、境界面が平面である場合には、上記式５のようになる。
そして、式５により抽出された反射波の到達時刻から、層１（図６（ａ）参照）の比誘電率εｒと深度ｄを、最小自乗誤差の最適解をLevenberg-Marquardt法により計算することができる。
得られた比誘電率εｒは、Ｔｏｐｐらによる比誘電率と含水比の関係式である以下の式６から、層１の含水比θｗを計算することができる。
［式６］
【0053】
(ステップ６)
さらに、ステップ６により、体積含水率計算処理が行われる。
具体的には、上記ステップ５で求められた比誘電率に基づいて、地中の含水比が算出される（本発明に係る含水比算出部）。
上記ステップ１～５により得られた比誘電率は、地層境界までの含水比と対応している。
したがって、比誘電率に基づいて、当該測定点における地層境界までの地中の含水比（体積含水率）を求めることができる。
また、既知の土質の場合には、対応した比誘電率と含水比関係を適応して含水比を調べることができる。
さらに、ステップ６により、体積含水率計算処理が行われる。
具体的には、上記ステップ５で求められた比誘電率に基づいて、地中の含水比が算出される（本発明に係る含水比算出部）。
上記ステップ１～５により得られた比誘電率は、地層境界までの含水比と対応している。
したがって、比誘電率に基づいて、当該測定点における地層境界までの地中の含水比（体積含水率）を求めることができる。
また、既知の土質の場合には、対応した比誘電率と含水比関係を適応して含水比を調べることができる。
【0054】
(ステップ７)
最後に、ステップ７により、体積含水分布図作成処理が行われる。
具体的には、上記ステップ１～６で求められた複数の各送受信アンテナ組に対応する複数の含水比に基づいて、含水比分布図が生成される（本発明に係る含水比分布図生成部）。
上記ステップ１～６が繰り返し実行されることにより、同一中点を持つ複数の送受信アンテナ組の測定によって（後述する図６，図９参照）、その中点位置における含水比が得られる。
また、送受信アンテナ組の中点の位置を異ならせて（後述する図８参照）、他の中点の場合についても上記ステップ１～６を繰り返し実行する。
そして、探査車両１０が走行されることにより、車両進行方向における複数の測定点（図４（ａ）参照）と、車両幅方向における複数の測定点（図４（ｂ）参照）での測定が行われ、さらに、複数の異なる中点の位置情報を組み合わせることができる。
これによって、複数走行による測定結果を重ね合わせ含水比分布図（図４（ｃ）参照）が得られる。
最後に、ステップ７により、体積含水分布図作成処理が行われる。
具体的には、上記ステップ１～６で求められた複数の各送受信アンテナ組に対応する複数の含水比に基づいて、含水比分布図が生成される（本発明に係る含水比分布図生成部）。
上記ステップ１～６が繰り返し実行されることにより、同一中点を持つ複数の送受信アンテナ組の測定によって（後述する図６，図９参照）、その中点位置における含水比が得られる。
また、送受信アンテナ組の中点の位置を異ならせて（後述する図８参照）、他の中点の場合についても上記ステップ１～６を繰り返し実行する。
そして、探査車両１０が走行されることにより、車両進行方向における複数の測定点（図４（ａ）参照）と、車両幅方向における複数の測定点（図４（ｂ）参照）での測定が行われ、さらに、複数の異なる中点の位置情報を組み合わせることができる。
これによって、複数走行による測定結果を重ね合わせ含水比分布図（図４（ｃ）参照）が得られる。
【0055】
［中点の複数設定］
上述したステップ１～７の含水比マッピング方法においては、説明の便宜上、複数の送受信アンテナ組における中点を一つに固定した場合について説明している。
但し、同じ送受信アンテナ組について、中点の位置を変えることで、異なる位置の含水比を求めることができる。
図８に、中点の位置を複数異ならせた場合の、中点と複数の送受信アンテナ組の位置を示す。
また、図９に、探査車両（地中レーダ）の移動に伴って移動する中点と複数の送受信アンテナ組の位置を示す。
上述したステップ１～７の含水比マッピング方法においては、説明の便宜上、複数の送受信アンテナ組における中点を一つに固定した場合について説明している。
但し、同じ送受信アンテナ組について、中点の位置を変えることで、異なる位置の含水比を求めることができる。
図８に、中点の位置を複数異ならせた場合の、中点と複数の送受信アンテナ組の位置を示す。
また、図９に、探査車両（地中レーダ）の移動に伴って移動する中点と複数の送受信アンテナ組の位置を示す。
【0056】
ここで、送受信アンテナの組についての方向成分としては、送受信アンテナのオフセット距離方向の成分と、それに垂直となる方向の成分とがある。
図８では、オフセット距離の方向に移動させた３組の中点をとる場合について示している。「Ａ１・Ｂ１」の組は「中点Ｃ１」の測定，「Ａ２・Ｂ２」の組は「中点Ｃ２」の測定，「Ａ３・Ｂ３」の組は「中点Ｃ３」の測定、となり、異なる３つの測定点（中点）についての測定が行えるようになる。
図８では、オフセット距離の方向に移動させた３組の中点をとる場合について示している。「Ａ１・Ｂ１」の組は「中点Ｃ１」の測定，「Ａ２・Ｂ２」の組は「中点Ｃ２」の測定，「Ａ３・Ｂ３」の組は「中点Ｃ３」の測定、となり、異なる３つの測定点（中点）についての測定が行えるようになる。
【0057】
また、送受信アンテナのオフセット距離と垂直な方向については、図９（ａ）に示すような、探査車両１０の進行方向に沿って複数の測定点（中点）についての測定が行えるようになる。
さらに、送受信アンテナのオフセット距離の方向は、図９（ａ）に示すように、各中点で常に同一（一直線）である場合に限られず、異なる方向であっても良い。
具体的には、図９（ｂ）に示すように、探査車両１０は道路などの路面のルートや形状に応じてカーブしながら走行することも一般的であり、探査車両１０の底面に配置された地中レーダ２０の送受信アンテナは、そのオフセット距離の方向も、それに応じて変化することになる。
さらに、送受信アンテナのオフセット距離の方向は、図９（ａ）に示すように、各中点で常に同一（一直線）である場合に限られず、異なる方向であっても良い。
具体的には、図９（ｂ）に示すように、探査車両１０は道路などの路面のルートや形状に応じてカーブしながら走行することも一般的であり、探査車両１０の底面に配置された地中レーダ２０の送受信アンテナは、そのオフセット距離の方向も、それに応じて変化することになる。
【0058】
図１０に、探査車両（地中レーダ）の移動に伴って移動する中点と複数の送受信アンテナ組の位置を示す。
図１０（ａ）は、探査車両が同じ路面を複数回（例えば２回）異なるルートで移動した場合を示している。
図１０（ｂ）は、図１０（ａ）の図から中点のみを代表して表示させた場合である。
これらの図からも明らかなように、各中点について含水比を求めることができ、それぞれの中点の位置情報と組み合わせることで、含水比分布を示すマッピング表示が可能となる。
図１０（ａ）は、探査車両が同じ路面を複数回（例えば２回）異なるルートで移動した場合を示している。
図１０（ｂ）は、図１０（ａ）の図から中点のみを代表して表示させた場合である。
これらの図からも明らかなように、各中点について含水比を求めることができ、それぞれの中点の位置情報と組み合わせることで、含水比分布を示すマッピング表示が可能となる。
【0059】
図１１に、本実施形態に係る含水比マッピング方法（装置）によって取得される含水比分布図（含水比マッピング）を模式的に示す。
図１１（ａ）は、図１０（ｂ）に対応する中点のみを表示させた場合であり、図１１（ｂ）はその中点とともに含水比の分布を等高線状に示した場合である。
また、図１１（ｃ）は、図１１（ｂ）から中点を取り除いて等高線のみを示した場合である。なお、図１１（ｂ），（ｃ）中で、等高線の色付き部分は含水比が○○～△△の値の領域、また、色の無い白色の領域は含水比が○○以下の領域、というように、含水比の値に応じた色分け表示をすることができる。
図１１（ａ）は、図１０（ｂ）に対応する中点のみを表示させた場合であり、図１１（ｂ）はその中点とともに含水比の分布を等高線状に示した場合である。
また、図１１（ｃ）は、図１１（ｂ）から中点を取り除いて等高線のみを示した場合である。なお、図１１（ｂ），（ｃ）中で、等高線の色付き部分は含水比が○○～△△の値の領域、また、色の無い白色の領域は含水比が○○以下の領域、というように、含水比の値に応じた色分け表示をすることができる。
【0060】
このように、測定点や含水比の分布を平面図上に示すことができ、また、含水比を等高線の色付き（色分け）表示で示すことができ、含水比の分布状態や広がり、大きさなどを、視覚により一目で把握・認識できるようにすることができる。
そして、このように生成・出力される含水比マッピングを示す地中探査データは、後述する地上三次元映像とともに、一元化された情報として生成・表示等することができるようになる（後述する図１４参照）。
そして、このように生成・出力される含水比マッピングを示す地中探査データは、後述する地上三次元映像とともに、一元化された情報として生成・表示等することができるようになる（後述する図１４参照）。
【0061】
なお、上述した本実施形態では、同日に測定を行う場合は勿論のこと、異なる日時において測定を行い、その異なる日時における含水比マッピングの差分により、含水比の変化を監視・把握することもできる。そのようにすることで、含水比マッピングの経年変化を監視・管理することにより、例えばその土地の水分動態の変化を知ることができるようになる。
【0062】
また、上述の例では、観測対象となる地中が二層（図６（ａ）の層１・層２参照）である場合について説明したが、三層以上の多数の層の場合には、状況がやや複雑になるものの、本実施形態の含水比マッピング方法によって測定することができる。
実際には、各層で異なる角度で電磁波の伝搬方向が変化するが、理解を容易にするため、各層が同じ角度で伝搬すると近似して扱うことととし、また、境界面は水平、各層の比誘電率は一様と想定して、以下に説明する。
複数の境界面がある場合には、複数の反射波が測定される。
この場合、最初に到達する信号は、多重反射を無視すると、地表に近い境界面からの反射に対応する。そこで、各オフセット距離に対して、電磁波反射の到達時間を測定し、各反射波の時刻を抽出して、回帰分析を行う。
これによって、各層の比誘電率と深度が求められ、その比誘電率の値を用いて、比誘電率－含水比関係から、含水比を計算により求めることができる。
実際には、各層で異なる角度で電磁波の伝搬方向が変化するが、理解を容易にするため、各層が同じ角度で伝搬すると近似して扱うことととし、また、境界面は水平、各層の比誘電率は一様と想定して、以下に説明する。
複数の境界面がある場合には、複数の反射波が測定される。
この場合、最初に到達する信号は、多重反射を無視すると、地表に近い境界面からの反射に対応する。そこで、各オフセット距離に対して、電磁波反射の到達時間を測定し、各反射波の時刻を抽出して、回帰分析を行う。
これによって、各層の比誘電率と深度が求められ、その比誘電率の値を用いて、比誘電率－含水比関係から、含水比を計算により求めることができる。
【0063】
［全方位カメラ３０］
全方位カメラ３０は、探査車両１０の例えば屋根上面などに配置され、探査車両１０が走行する路面上の地上三次元映像を生成する地上映像生成手段である。
具体的には、全方位カメラ３０は、走行する探査車両１０の周囲を３６０°撮影できる１台又は複数台のビデオカメラであり、探査車両１０の走行に伴って、上述した地中レーダ２０が取得する地中データに対応する、地上の三次元映像（全周映像）を撮像できるようになっている。
全方位カメラ３０は、探査車両１０の例えば屋根上面などに配置され、探査車両１０が走行する路面上の地上三次元映像を生成する地上映像生成手段である。
具体的には、全方位カメラ３０は、走行する探査車両１０の周囲を３６０°撮影できる１台又は複数台のビデオカメラであり、探査車両１０の走行に伴って、上述した地中レーダ２０が取得する地中データに対応する、地上の三次元映像（全周映像）を撮像できるようになっている。
【0064】
図１２に、本実施形態の全方位カメラ３０で取得・生成される地上３６０度の立体全周映像の一例を示す。
全方位カメラ３０は、周囲３６０°の全周映像を撮像できるため、例えば図１２（ａ）に示すように、路面を走行する探査車両１０を中心とした半径約５ｍの地上の状態・風景が撮像される。同図の例では、全幅員１４．５５ｍ，有効幅員１３．７５ｍの道路や橋梁等で、片側車線の中央を走行する探査車両１０を中心とした車道部３．５ｍ，路肩０．５ｍ，側帯０．５ｍ，歩道部２．０ｍ，分離帯０．７５，地覆０．４ｍの範囲について、３６０°の全周映像が撮像・取得できることが分かる。
また、図１２（ｂ）に示すように、路面を走行する探査車両１０を中心とした半径約５～２５ｍの範囲で地上の状態・風景が撮像され、路面横にある斜面（法面）についても３６０°の全周映像が撮像・取得できることが分かる。
図１２（ｃ）は、全方位カメラ３０で撮像された３６０°の全周映像の出力例である。
全方位カメラ３０は、周囲３６０°の全周映像を撮像できるため、例えば図１２（ａ）に示すように、路面を走行する探査車両１０を中心とした半径約５ｍの地上の状態・風景が撮像される。同図の例では、全幅員１４．５５ｍ，有効幅員１３．７５ｍの道路や橋梁等で、片側車線の中央を走行する探査車両１０を中心とした車道部３．５ｍ，路肩０．５ｍ，側帯０．５ｍ，歩道部２．０ｍ，分離帯０．７５，地覆０．４ｍの範囲について、３６０°の全周映像が撮像・取得できることが分かる。
また、図１２（ｂ）に示すように、路面を走行する探査車両１０を中心とした半径約５～２５ｍの範囲で地上の状態・風景が撮像され、路面横にある斜面（法面）についても３６０°の全周映像が撮像・取得できることが分かる。
図１２（ｃ）は、全方位カメラ３０で撮像された３６０°の全周映像の出力例である。
【0065】
そして、本実施形態では、全方位カメラ３０で取得・生成される３６０°の全周映像に基づいて、所定のカメラベクトル演算を行うことで、現地での測量を行うことなく、また、ＧＰＳの位置情報に依存することなく、探査車両１０に全方位カメラ３０を取り付けて走行（例えば最大速度６０ｋｍ／ｈ程度）するだけで、簡単に高精度な三次元地上マッピングができるようになっている。
一般に、探査車両を構成する自動車等に搭載されているＧＰＳや４方向カメラによる位置情報は、ＧＰＳの受信状況により例えば０．５ｍ～数ｍ程度の誤差が出る。このため、それだけでは任意の異常箇所等の位置特定を行うことは困難である。また、トンネルの中や高架橋の下ではＧＰＳが受信できないので、位置情報は車両に搭載された４方向カメラの映像のみとなり、任意の箇所の位置特定は更に困難となる。
そこで、本実施形態では、全方位カメラ３０で取得された３６０°全周映像に基づいて、ＣＶ演算部４０において所定のカメラベクトル演算を行うことで、ＧＰＳに依ることなく、また、現地での測量・計測等を行うことなく、高精度な三次元地上マッピングを行うようになっている。
一般に、探査車両を構成する自動車等に搭載されているＧＰＳや４方向カメラによる位置情報は、ＧＰＳの受信状況により例えば０．５ｍ～数ｍ程度の誤差が出る。このため、それだけでは任意の異常箇所等の位置特定を行うことは困難である。また、トンネルの中や高架橋の下ではＧＰＳが受信できないので、位置情報は車両に搭載された４方向カメラの映像のみとなり、任意の箇所の位置特定は更に困難となる。
そこで、本実施形態では、全方位カメラ３０で取得された３６０°全周映像に基づいて、ＣＶ演算部４０において所定のカメラベクトル演算を行うことで、ＧＰＳに依ることなく、また、現地での測量・計測等を行うことなく、高精度な三次元地上マッピングを行うようになっている。
【0066】
［ＣＶ演算部４０］
ＣＶ演算部４０は、上述した全方位カメラ３０で撮像された地上三次元映像の画像データから、高精度の位置情報を演算により算出・生成することができるカメラベクトル演算手段である。
具体的には、本実施形態に係るＣＶ演算部４０は、所定数の特徴点を自動抽出する特徴点抽出部４１と、抽出された特徴点について、動画映像の各フレーム画像内で自動追跡してフレーム画像間での対応関係を求める特徴点対応処理部４２と、対応関係が求められた特徴点の三次元位置座標を求め、当該三次元位置座標から、各フレーム画像に対応したカメラの三次元位置座標及び三次元回転座標からなるカメラベクトルを求めるカメラベクトル演算部４３などの手段で構成される（後述する図１５参照）。
カメラベクトル演算は、本出願人に係る特許第４４４６０４１号公報等で開示している技術であり、３６０°全周映像の中に特徴点を自動抽出し、それを隣接する複数フレームにトラッキングし、カメラ移動のベースとトラッキング点で校正する三角形を構成し、そのトラッキングデータを解析することで、特徴点の三次元座標と、カメラ位置と姿勢を持つ全周ＣＶ（カメラベクトル）映像を取得するものである。
ＣＶ演算部４０は、上述した全方位カメラ３０で撮像された地上三次元映像の画像データから、高精度の位置情報を演算により算出・生成することができるカメラベクトル演算手段である。
具体的には、本実施形態に係るＣＶ演算部４０は、所定数の特徴点を自動抽出する特徴点抽出部４１と、抽出された特徴点について、動画映像の各フレーム画像内で自動追跡してフレーム画像間での対応関係を求める特徴点対応処理部４２と、対応関係が求められた特徴点の三次元位置座標を求め、当該三次元位置座標から、各フレーム画像に対応したカメラの三次元位置座標及び三次元回転座標からなるカメラベクトルを求めるカメラベクトル演算部４３などの手段で構成される（後述する図１５参照）。
カメラベクトル演算は、本出願人に係る特許第４４４６０４１号公報等で開示している技術であり、３６０°全周映像の中に特徴点を自動抽出し、それを隣接する複数フレームにトラッキングし、カメラ移動のベースとトラッキング点で校正する三角形を構成し、そのトラッキングデータを解析することで、特徴点の三次元座標と、カメラ位置と姿勢を持つ全周ＣＶ（カメラベクトル）映像を取得するものである。
【0067】
そして、ＣＶ演算で得られるＣＶ（カメラベクトル）値は、カメラ（全方位カメラ３０）の位置と姿勢を６自由度で表現したものである。すなわち、全方位カメラ３０で撮像される３６０°全周映像の全フレームについて、ｘ，ｙ，ｚ，θｘ，θｙ，θｚ，の６自由度のそれぞれの係数を演算で決定するものであり、ＧＰＳや測量等による位置情報に依ることなく、３６０°映像中に相対的な位置情報・三次元座標情報を演算のみによって生成・付与できるものである。
これによって、任意の全周映像について、ＣＶ値を生成・付与することができ、このＣＶ値に基づいて、全周映像と、上述した地中レーダ２０で取得される地中三次元情報とを一元化する三次元情報一元化処理が実行される。
このＣＶ演算部４０におけるＣＶ演算の具体的な詳細内容については、図１５～１７を参照しつつ後述する。
これによって、任意の全周映像について、ＣＶ値を生成・付与することができ、このＣＶ値に基づいて、全周映像と、上述した地中レーダ２０で取得される地中三次元情報とを一元化する三次元情報一元化処理が実行される。
このＣＶ演算部４０におけるＣＶ演算の具体的な詳細内容については、図１５～１７を参照しつつ後述する。
【0068】
［一元化処理部５０］
一元化処理部５０は、地中レーダ２０で生成される地中三次元情報と、全方位カメラ３０で生成される地上三次元映像を、所定のデータ一元化処理により、地中から地上空間まで一元化された三次元情報として生成する三次元情報一元化処理手段である。
具体的には、本実施形態に係る一元化処理部５０は、上記のＣＶ演算部４０で生成されるカメラベクトルと、地中レーダ２０で取得された地中三次元情報が有する位置座標情報に基づいて、地中三次元情報と地上三次元映像を一元化する。
一元化処理部５０は、地中レーダ２０で生成される地中三次元情報と、全方位カメラ３０で生成される地上三次元映像を、所定のデータ一元化処理により、地中から地上空間まで一元化された三次元情報として生成する三次元情報一元化処理手段である。
具体的には、本実施形態に係る一元化処理部５０は、上記のＣＶ演算部４０で生成されるカメラベクトルと、地中レーダ２０で取得された地中三次元情報が有する位置座標情報に基づいて、地中三次元情報と地上三次元映像を一元化する。
【0069】
図１３は、本実施形態の一元化処理部５０による地中探査データと地上映像の一元化処理の手順を示すフローチャートである。
同図に示すように、本実施形態の一元化処理では、まず、地中レーダ２０で取得された地中探査データ（ステップ１）と、全方位カメラ３０で取得された地上の３６０°全周映像（ステップ２）が、一元化処理部５０に入力される。
このとき、地中レーダ２０で取得された地中探査データには、探査車両１０に搭載されたＧＰＳ（高精度ＧＮＳＳ）１１で取得される位置情報や時間情報などに基づく位置情報が含まれている。
全方位カメラ３０の３６０°全周映像には、探査車両１０に搭載されたＧＰＳ（高精度ＧＮＳＳ）１１で取得される位置情報や時間情報などに基づく位置情報が含まれ、また、上述したＣＶ演算部４０において演算処理され（ステップ３）、三次元一座標となるＣＶ値が付与されている。
同図に示すように、本実施形態の一元化処理では、まず、地中レーダ２０で取得された地中探査データ（ステップ１）と、全方位カメラ３０で取得された地上の３６０°全周映像（ステップ２）が、一元化処理部５０に入力される。
このとき、地中レーダ２０で取得された地中探査データには、探査車両１０に搭載されたＧＰＳ（高精度ＧＮＳＳ）１１で取得される位置情報や時間情報などに基づく位置情報が含まれている。
全方位カメラ３０の３６０°全周映像には、探査車両１０に搭載されたＧＰＳ（高精度ＧＮＳＳ）１１で取得される位置情報や時間情報などに基づく位置情報が含まれ、また、上述したＣＶ演算部４０において演算処理され（ステップ３）、三次元一座標となるＣＶ値が付与されている。
【0070】
そして、これらの地中探査データと３６０°全周映像とが、対応する位置情報に基づいて関連付けられ、ＣＶ演算部４０で演算・付与されたＣＶ値が、３６０°全周映像と対応する地中探査データに付与され、同一の位置座標によって紐付けられる（ステップ４）。
これにより、地中探査データは、位置座標としてＣＶ値が与えられた高精度マッピングされた地中三次元情報となる（ステップ５）。
また、３６０°全周映像は、オルソ画像処理されて、ＣＶ値が付与された高精度のオルソ画像からなる地上三次元映像として生成される（ステップ６）。
これにより、地中探査データは、位置座標としてＣＶ値が与えられた高精度マッピングされた地中三次元情報となる（ステップ５）。
また、３６０°全周映像は、オルソ画像処理されて、ＣＶ値が付与された高精度のオルソ画像からなる地上三次元映像として生成される（ステップ６）。
【0071】
ここで、オルソ画像の生成は、公知の画像生成技術を用いることができるが、本実施形態では、ＣＶ値が付与された全周映像に基づくことで、高精度なオルソ画像を生成・出力することができる。具体的には、ＣＶ値を位置情報として、３６０°全周映像を中心投影から真上から見たような歪みがない正射投影画像に変換し、複数の正射画像をつなぎ目が目立たないよう接合（モザイク）して、統合された一枚のオルソ画像を生成し、所望の範囲や大きさに分割されたオルソ画像データを生成・記憶することができる。
その後、対応するＣＶ値に基づいて、地中三次元情報とオルソ画像からなる地上三次元映像が一元化処理され、地中地上一元化データとして生成・出力される（ステップ７）。
その後、対応するＣＶ値に基づいて、地中三次元情報とオルソ画像からなる地上三次元映像が一元化処理され、地中地上一元化データとして生成・出力される（ステップ７）。
【0072】
以上のような地中三次元情報と地上オルソ画像の一元化処理による高精度三次元マッピングにより、地中及び地上データの三次元位置座標を、ＧＰＳの受信状況等に左右されることなく、ＣＶ演算に基づく約５～１５ｃｍ程度の誤差で取得することができる。
また、オルソ画像は、真上から見たような歪のない画像に変換し、位置情報を付与したもので、路面のオルソ画像は三次元座標データが付加された全方位カメラの映像から生成することができる。このようなオルソ画像は、平面地図と同様の精度を持つため、三次元地中レーダ結果を表示するのに適した画像となる。
また、オルソ画像は、真上から見たような歪のない画像に変換し、位置情報を付与したもので、路面のオルソ画像は三次元座標データが付加された全方位カメラの映像から生成することができる。このようなオルソ画像は、平面地図と同様の精度を持つため、三次元地中レーダ結果を表示するのに適した画像となる。
【0073】
図１４に、地中三次元データ及び地中三次元データを地上オルソ画像に一元化処理した画像の一例を示す。
図１４（ａ）及び（ｂ）は、地中レーダ２０により得られた含水比の分布を等高線状に示した情報と、この地下情報に対応する路面等の地上情報を、ＣＶ演算部４０で求められた三次元位置情報から生成した路面のオルソ画像に一元化した出力結果である。これによって、含水比の分布と路面の状況等を、画像上で対比して視認・解析等することができる。
また、このような対比情報により、路面の状態・状況（例えば、ひび割れや陥没，腐食等）の位置と地中の含水比分布の位置が明確となり、次回調査や補修工事後の時系列のデータ変化を視覚的に比較し、データベース化することができる。また、異常箇所のデータベース化は、例えば補修計画等の維持管理に役立つ有用な情報となる。
このような本実施形態に係る一元化処理画像は、例えばＰＣ等の情報処理装置で任意に加工・編集したり、表示装置（ディスプレイ）で表示させることができ、また、プリンタ等の印刷装置で印刷することができる。
図１４（ａ）及び（ｂ）は、地中レーダ２０により得られた含水比の分布を等高線状に示した情報と、この地下情報に対応する路面等の地上情報を、ＣＶ演算部４０で求められた三次元位置情報から生成した路面のオルソ画像に一元化した出力結果である。これによって、含水比の分布と路面の状況等を、画像上で対比して視認・解析等することができる。
また、このような対比情報により、路面の状態・状況（例えば、ひび割れや陥没，腐食等）の位置と地中の含水比分布の位置が明確となり、次回調査や補修工事後の時系列のデータ変化を視覚的に比較し、データベース化することができる。また、異常箇所のデータベース化は、例えば補修計画等の維持管理に役立つ有用な情報となる。
このような本実施形態に係る一元化処理画像は、例えばＰＣ等の情報処理装置で任意に加工・編集したり、表示装置（ディスプレイ）で表示させることができ、また、プリンタ等の印刷装置で印刷することができる。
【0074】
以上のような本実施形態で得られる地中三次元データを地上オルソ画像に一元化した一元化処理画像により、以下のような効果が得られる。
(1)高精度三次元座標が得られるため、二次元地中レーダによる位置特定が不要になる。
(2)位置補正装置により、トンネルや高架橋のＧＰＳが受信できない場所でも高精度の三次元座標を得ることができる。
(3)効率よく地上の高精度三次元座標が得られるので、現地で構造物の形状や道路周辺物などを計測する必要がない。
(4)三次元地中レーダの地下情報と、全方位カメラ３０及びＣＶ演算部４０の地上情報のデータベース化で一元管理ができる。
(5)高精度のオルソ画像による視覚情報から、次回調査や補修工事後との比較といった情報の時系列管理ができる。
(6)オルソ画像は平面地図と同様の精度を持ちながら、視覚的に地図より優れているので、発注者に効果的な成果が提供できる。
(7)調査の効率化、一元化に及び映像から簡単に高精度の三次元座標情報が得られることから、災害時の供給道路の安全確保などの緊急時に迅速な対応ができる。
(1)高精度三次元座標が得られるため、二次元地中レーダによる位置特定が不要になる。
(2)位置補正装置により、トンネルや高架橋のＧＰＳが受信できない場所でも高精度の三次元座標を得ることができる。
(3)効率よく地上の高精度三次元座標が得られるので、現地で構造物の形状や道路周辺物などを計測する必要がない。
(4)三次元地中レーダの地下情報と、全方位カメラ３０及びＣＶ演算部４０の地上情報のデータベース化で一元管理ができる。
(5)高精度のオルソ画像による視覚情報から、次回調査や補修工事後との比較といった情報の時系列管理ができる。
(6)オルソ画像は平面地図と同様の精度を持ちながら、視覚的に地図より優れているので、発注者に効果的な成果が提供できる。
(7)調査の効率化、一元化に及び映像から簡単に高精度の三次元座標情報が得られることから、災害時の供給道路の安全確保などの緊急時に迅速な対応ができる。
【0075】
［ＣＶ演算処理］
次に、上述した本実施形態に係る地中探査装置のＣＶ演算部４０で実行されるＣＶ演算の具体的な内容について、図１５～図２５を参照しつつ説明する。
上記のとおり、ＣＶ演算とはＣＶ値を求める方法の一つであり、ＣＶ演算により求められた結果をＣＶ値，ＣＶデータと呼ぶ。ＣＶという表記は、カメラベクトル：CameraVectorの略記であり、カメラベクトル（ＣＶ）とは計測等のために映像を取得するビデオカメラ等のカメラの三次元位置と３軸回転姿勢を示す値である。
ＣＶ演算は、動画像（ビデオ映像）を取得し、その映像内の特徴点を検出し、それを隣接する複数のフレームに追跡し、カメラ位置と特徴点の追跡軌跡とが作る三角形を画像内に数多く生成し、その三角形を解析することで、カメラの三次元位置とカメラの３軸回転姿勢を求めるものである。
次に、上述した本実施形態に係る地中探査装置のＣＶ演算部４０で実行されるＣＶ演算の具体的な内容について、図１５～図２５を参照しつつ説明する。
上記のとおり、ＣＶ演算とはＣＶ値を求める方法の一つであり、ＣＶ演算により求められた結果をＣＶ値，ＣＶデータと呼ぶ。ＣＶという表記は、カメラベクトル：CameraVectorの略記であり、カメラベクトル（ＣＶ）とは計測等のために映像を取得するビデオカメラ等のカメラの三次元位置と３軸回転姿勢を示す値である。
ＣＶ演算は、動画像（ビデオ映像）を取得し、その映像内の特徴点を検出し、それを隣接する複数のフレームに追跡し、カメラ位置と特徴点の追跡軌跡とが作る三角形を画像内に数多く生成し、その三角形を解析することで、カメラの三次元位置とカメラの３軸回転姿勢を求めるものである。
【0076】
ＣＶ演算では、ＣＶ値を求める過程で、同時に映像内の特徴点（基準点）についても三次元座標が同時に求まることが重要な特性である。
また、動画像から演算で求められるＣＶ値は、動画像の各フレームに対応して、三次元のカメラ位置と三次元のカメラ姿勢とが同時に求まる。しかも、原理的には一台のカメラで、映像と対応してＣＶ値が求められる特性は、ＣＶ演算でしか実現し得ない、優れた特徴である。
但し、ＣＶ値は、例えば、他の方法による計測手段（ＧＰＳとジャイロ、又はＩＭＵ等）より求めることもできる。従って、ＣＶ値の取得はＣＶ演算による場合には限られない。
なお、ＧＰＳとジャイロ、又はＩＭＵ等によりＣＶ値を取得する場合には、動画像の各フレームと、その三次元的カメラ位置と三次元的カメラ姿勢とを同時に取得するためには画像フレームと計測サンプリング時刻を高精度で、かつ完全に同期させる必要がある。
また、動画像から演算で求められるＣＶ値は、動画像の各フレームに対応して、三次元のカメラ位置と三次元のカメラ姿勢とが同時に求まる。しかも、原理的には一台のカメラで、映像と対応してＣＶ値が求められる特性は、ＣＶ演算でしか実現し得ない、優れた特徴である。
但し、ＣＶ値は、例えば、他の方法による計測手段（ＧＰＳとジャイロ、又はＩＭＵ等）より求めることもできる。従って、ＣＶ値の取得はＣＶ演算による場合には限られない。
なお、ＧＰＳとジャイロ、又はＩＭＵ等によりＣＶ値を取得する場合には、動画像の各フレームと、その三次元的カメラ位置と三次元的カメラ姿勢とを同時に取得するためには画像フレームと計測サンプリング時刻を高精度で、かつ完全に同期させる必要がある。
【0077】
動画像から演算で求められるＣＶデータは、加工しない段階では相対値であるが、短区間であれば高精度で三次元位置情報と３軸回転の角度情報を取得できる。
また、ＣＶデータを画像から取得する場合は、取得されたデータは相対値であり、絶対値に変換するには既知の地点でのキャリブレーションが必要であるが、ＣＶ値と同時に取得される特徴点三次元位置座標と姿勢により、又は新たに画像内に指定して取得した特徴点により、画像内の任意の対象物との位置関係を計測することができるという、他の方法では実現困難な優れた特性を備える。
また、画像に対応したＣＶ値が求まるので、画像との相性が良く、画像から取得したＣＶ値である限り、その画像内では矛盾が無いので、画像以外から求めるＣＶ値と比較して、画像から直接にカメラ位置とその３軸回転姿勢を求めることができるＣＶ演算は画像内計測や画像内測量に好適となる。
そして、本発明では、このＣＶ演算により得られたＣＶ値データに基づいて、上述した地上オルソ画像の生成、地中三次元データと地上オルソ画像の一元化処理画像を生成するものである。
また、ＣＶデータを画像から取得する場合は、取得されたデータは相対値であり、絶対値に変換するには既知の地点でのキャリブレーションが必要であるが、ＣＶ値と同時に取得される特徴点三次元位置座標と姿勢により、又は新たに画像内に指定して取得した特徴点により、画像内の任意の対象物との位置関係を計測することができるという、他の方法では実現困難な優れた特性を備える。
また、画像に対応したＣＶ値が求まるので、画像との相性が良く、画像から取得したＣＶ値である限り、その画像内では矛盾が無いので、画像以外から求めるＣＶ値と比較して、画像から直接にカメラ位置とその３軸回転姿勢を求めることができるＣＶ演算は画像内計測や画像内測量に好適となる。
そして、本発明では、このＣＶ演算により得られたＣＶ値データに基づいて、上述した地上オルソ画像の生成、地中三次元データと地上オルソ画像の一元化処理画像を生成するものである。
【0078】
［ＣＶ演算部］
ＣＶ演算部４０は、上述した全方位カメラ３０で撮像されるビデオ映像について所定のＣＶ演算処理を行うことでＣＶ値を求めるようになっており、具体的には、図１５に示すように、特徴点抽出部４１と、特徴点対応処理部４２と、カメラベクトル演算部４３と、誤差最小化部４４と、三次元情報追跡部４５と、高精度カメラベクトル演算部４６とを備えている。
ＣＶ演算部４０は、上述した全方位カメラ３０で撮像されるビデオ映像について所定のＣＶ演算処理を行うことでＣＶ値を求めるようになっており、具体的には、図１５に示すように、特徴点抽出部４１と、特徴点対応処理部４２と、カメラベクトル演算部４３と、誤差最小化部４４と、三次元情報追跡部４５と、高精度カメラベクトル演算部４６とを備えている。
【0079】
まず、ＣＶ演算に使用する映像としては、原理的にはどのような映像でもよいが、画角の限られた映像では視点方向を移動した場合に映像がとぎれてしまうので、実用的には全周映像（図１２（ｃ），図１６参照）か、広角映像を全周映像の一部として扱うことが望ましい。なお、動画映像は連続する静止画と同様であり、複数の静止画として扱うことができる。
また、映像は、一般には予め記録した動画映像を使うことになるが、自動車等の移動体の移動に合わせてリアルタイムに取り込んだ映像を使用することも勿論可能である。
また、映像は、一般には予め記録した動画映像を使うことになるが、自動車等の移動体の移動に合わせてリアルタイムに取り込んだ映像を使用することも勿論可能である。
【0080】
そこで、本実施形態では、ＣＶ演算に使用する映像として、車両等の移動体の３６０度の全周囲を撮影した全周映像（図１２（ｃ），図１６参照）か、又は全周映像に近い広角映像を用いて、その全周映像を視点方向に平面展開することにより、全周映像の一部として扱うことができる。
全周映像の平面展開とは、全周映像を、通常の画像として遠近法的に表現するものである。ここで、「遠近法」と呼称するのは、全周画像のそのものはメルカトール図法や球面投影図法のように、遠近法とは異なる方法で表示されているので（図１６参照）、これを例えば等距離図法のように平面展開表示することで、通常の遠近法映像に変換表示できるからである。
全周映像の平面展開とは、全周映像を、通常の画像として遠近法的に表現するものである。ここで、「遠近法」と呼称するのは、全周画像のそのものはメルカトール図法や球面投影図法のように、遠近法とは異なる方法で表示されているので（図１６参照）、これを例えば等距離図法のように平面展開表示することで、通常の遠近法映像に変換表示できるからである。
【0081】
本実施形態においては、まず、探査車両１０を走行させて、ＣＶ値データを取得する目的で、探査車両１０に固定された全方位カメラ３０で、探査車両１０の移動とともに探査車両１０の周辺の３６０°全周映像（図１２（ｃ），図１６参照）を撮影する。
なお、探査車両１０には、その位置座標を取得する目的で、例えば、絶対座標を取得するＧＰＳ機器単独やＩＭＵ機器を付加したもの等により構成した位置計測機器等を備えることができる。
また、探査車両１０に搭載される全方位カメラ３０としては、広範囲映像を撮影，取得するカメラであればどのような構成であってもよく、例えば、広角レンズや魚眼レンズ付きカメラ、移動カメラ、固定カメラ、複数のカメラを固定したカメラ、３６０度周囲に回転可能なカメラ等がある。本実施形態では、上述したように、探査車両１０に一台又は複数台のカメラが一体的に固定され、探査車両１０の移動に伴って広範囲映像を撮影する全方位カメラ３０を使用している。
なお、探査車両１０には、その位置座標を取得する目的で、例えば、絶対座標を取得するＧＰＳ機器単独やＩＭＵ機器を付加したもの等により構成した位置計測機器等を備えることができる。
また、探査車両１０に搭載される全方位カメラ３０としては、広範囲映像を撮影，取得するカメラであればどのような構成であってもよく、例えば、広角レンズや魚眼レンズ付きカメラ、移動カメラ、固定カメラ、複数のカメラを固定したカメラ、３６０度周囲に回転可能なカメラ等がある。本実施形態では、上述したように、探査車両１０に一台又は複数台のカメラが一体的に固定され、探査車両１０の移動に伴って広範囲映像を撮影する全方位カメラ３０を使用している。
【0082】
そして、以上のような全方位カメラ３０によれば、探査車両１０の屋根上部等に設置されることで、カメラの３６０度全周囲の映像を一台又は複数のカメラで同時に撮影することができ、探査車両１０が走行・移動することで、広範囲映像を動画データとして取得できる。
ここで、全方位カメラ３０は、カメラの全周映像を直接取得できるビデオカメラであるが、カメラの全周囲の半分以上を映像として取得できれば全周映像として使用できる。
また、画角が制限された通常のカメラの場合でも、ＣＶ演算の精度としては低下するが、全周映像の一部分として取り扱うことが可能である。
ここで、全方位カメラ３０は、カメラの全周映像を直接取得できるビデオカメラであるが、カメラの全周囲の半分以上を映像として取得できれば全周映像として使用できる。
また、画角が制限された通常のカメラの場合でも、ＣＶ演算の精度としては低下するが、全周映像の一部分として取り扱うことが可能である。
【0083】
なお、全方位カメラ３０で撮影された広範囲映像は、一枚の画像として、撮影時の画角に一致する仮想球面に貼り付けることができる。
仮想球面に貼り付けられた球面画像データは、仮想球面に貼り付けた状態の球面画像（３６０度画像）データとして保存・出力される。仮想球面は、広範囲映像を取得するカメラ部を中心点とした任意の球面状に設定することができる。
図１６（ａ）は球面画像が貼り付けられる仮想球面の外観イメージであり、同図（ｂ）は仮想球面に貼り付けられた球面画像の一例である。また、同図（ｃ）は、（ｂ）の球面画像をメルカトール図法に従って平面展開した画像例を示す。
仮想球面に貼り付けられた球面画像データは、仮想球面に貼り付けた状態の球面画像（３６０度画像）データとして保存・出力される。仮想球面は、広範囲映像を取得するカメラ部を中心点とした任意の球面状に設定することができる。
図１６（ａ）は球面画像が貼り付けられる仮想球面の外観イメージであり、同図（ｂ）は仮想球面に貼り付けられた球面画像の一例である。また、同図（ｃ）は、（ｂ）の球面画像をメルカトール図法に従って平面展開した画像例を示す。
【0084】
そして、以上のように生成・取得された全周ビデオ映像が、ＣＶ演算部４０に入力されてＣＶ値データが求められる（図１５参照）。
ＣＶ演算部４０では、まず、特徴点抽出部４１が、全方位カメラ３０で撮影されて一時記録された動画像データの中から、十分な数の特徴点（基準点）を自動抽出する。
特徴点対応処理部４２は、自動抽出された特徴点を、各フレーム間で各フレーム画像内において自動的に追跡することで、その対応関係を自動的に求める。
カメラベクトル演算部４３は、対応関係が求められた特徴点の三次元位置座標から各フレーム画像に対応したカメラベクトルを演算で自動的に求める。
誤差最小化部４４は、複数のカメラ位置の重複演算により、各カメラベクトルの解の分布が最小になるように統計処理し、誤差の最小化処理を施したカメラ位置方向を自動的に決定する。
ＣＶ演算部４０では、まず、特徴点抽出部４１が、全方位カメラ３０で撮影されて一時記録された動画像データの中から、十分な数の特徴点（基準点）を自動抽出する。
特徴点対応処理部４２は、自動抽出された特徴点を、各フレーム間で各フレーム画像内において自動的に追跡することで、その対応関係を自動的に求める。
カメラベクトル演算部４３は、対応関係が求められた特徴点の三次元位置座標から各フレーム画像に対応したカメラベクトルを演算で自動的に求める。
誤差最小化部４４は、複数のカメラ位置の重複演算により、各カメラベクトルの解の分布が最小になるように統計処理し、誤差の最小化処理を施したカメラ位置方向を自動的に決定する。
【0085】
三次元情報追跡部４５は、カメラベクトル演算部４３で得られたカメラベクトルを概略のカメラベクトルと位置づけ、その後のプロセスで順次画像の一部として得られる三次元情報に基づいて、複数のフレーム画像に含まれる部分的三次元情報を隣接するフレームの画像に沿って自動追跡を行う。ここで、三次元情報（三次元形状）とは、主に特徴点の三次元分布情報であり、すなわち、三次元の点の集まりであり、この三次元の点の集まりが三次元形状を構成する。
高精度カメラベクトル演算部４６は、三次元情報追跡部４５で得られた追跡データに基づいて、カメラベクトル演算部４３で得られるカメラベクトルより、さらに高精度なカメラベクトルを生成，出力する。
そして、以上のようにして得られたカメラベクトルが、上述した一元化処理部５０に入力され、地上オルソ画像の生成、地中三次元データと地上オルソ画像の一元化処理に利用されることになる。
高精度カメラベクトル演算部４６は、三次元情報追跡部４５で得られた追跡データに基づいて、カメラベクトル演算部４３で得られるカメラベクトルより、さらに高精度なカメラベクトルを生成，出力する。
そして、以上のようにして得られたカメラベクトルが、上述した一元化処理部５０に入力され、地上オルソ画像の生成、地中三次元データと地上オルソ画像の一元化処理に利用されることになる。
【0086】
複数の画像（動画又は連続静止画）の特徴点からカメラベクトルを検出するにはいくつかの方法があるが、図１５に示す本実施形態のＣＶ演算部４０では、画像内に十分に多くの数の特徴点を自動抽出し、それを自動追跡することで、三角測量幾何学、及びエピポーラ幾何学により、カメラの三次元ベクトル及び３軸回転ベクトルを求めるようにしてある。
特徴点を充分に多くとることにより、カメラベクトル情報が重複することで、三地点の座標で構成される単位の三角形を、その十倍以上の多くの特徴点から求めることになり、重複する情報から誤差を最小化させて、より精度の高いカメラベクトルを求めることができる。
特徴点を充分に多くとることにより、カメラベクトル情報が重複することで、三地点の座標で構成される単位の三角形を、その十倍以上の多くの特徴点から求めることになり、重複する情報から誤差を最小化させて、より精度の高いカメラベクトルを求めることができる。
【0087】
カメラベクトルとは、移動するカメラの持つ六個の自由度のベクトルである。
一般に、静止した三次元物体は、位置座標（Ｘ，Ｙ，Ｚ）と、それぞれの座標軸の回転角（θｘ，θｙ，θｚ）の六個の自由度を持つ。従って、カメラベクトルは、カメラの位置座標（Ｘ，Ｙ，Ｚ）とそれぞれの座標軸の回転角（θｘ，θｙ，θｚ）の六個の自由度のベクトルをいう。なお、カメラが移動する場合は、自由度に移動方向も入るが、これは上記の六個の自由度から微分して導き出すことができる。
このように、本実施形態のカメラベクトルの検出とは、カメラは各フレーム毎に六個の自由度の値をとり、各フレーム毎に異なる六個の係数を決定することである。
一般に、静止した三次元物体は、位置座標（Ｘ，Ｙ，Ｚ）と、それぞれの座標軸の回転角（θｘ，θｙ，θｚ）の六個の自由度を持つ。従って、カメラベクトルは、カメラの位置座標（Ｘ，Ｙ，Ｚ）とそれぞれの座標軸の回転角（θｘ，θｙ，θｚ）の六個の自由度のベクトルをいう。なお、カメラが移動する場合は、自由度に移動方向も入るが、これは上記の六個の自由度から微分して導き出すことができる。
このように、本実施形態のカメラベクトルの検出とは、カメラは各フレーム毎に六個の自由度の値をとり、各フレーム毎に異なる六個の係数を決定することである。
【0088】
以下、ＣＶ演算部４０における具体的なカメラベクトルの検出方法について、図１７以下を参照しつつ説明する。
まず、上述した全方位カメラ３０で取得された画像データは、間接に又は直接に、ＣＶ演算部４０の特徴点抽出部４１に入力され、特徴点抽出部４１で、適切にサンプリングされたフレーム画像中に、特徴点となるべき点又は小領域画像が自動抽出され、特徴点対応処理部４２で、複数のフレーム画像間で特徴点の対応関係が自動的に求められる。
具体的には、カメラベクトルの検出の基準となる、十分に必要な数以上の特徴点を求める。画像間の特徴点とその対応関係の一例を、図１７～図１９に示す。図中「＋」が自動抽出された特徴点であり、複数のフレーム画像間で対応関係が自動追跡される（図１９に示す対応点１～４参照）。
ここで、特徴点の抽出は、図２０に示すように、各画像中に充分に多くの特徴点を指定，抽出することが望ましく（図２０の○印参照）、例えば、１００点程度の特徴点を抽出する。
まず、上述した全方位カメラ３０で取得された画像データは、間接に又は直接に、ＣＶ演算部４０の特徴点抽出部４１に入力され、特徴点抽出部４１で、適切にサンプリングされたフレーム画像中に、特徴点となるべき点又は小領域画像が自動抽出され、特徴点対応処理部４２で、複数のフレーム画像間で特徴点の対応関係が自動的に求められる。
具体的には、カメラベクトルの検出の基準となる、十分に必要な数以上の特徴点を求める。画像間の特徴点とその対応関係の一例を、図１７～図１９に示す。図中「＋」が自動抽出された特徴点であり、複数のフレーム画像間で対応関係が自動追跡される（図１９に示す対応点１～４参照）。
ここで、特徴点の抽出は、図２０に示すように、各画像中に充分に多くの特徴点を指定，抽出することが望ましく（図２０の○印参照）、例えば、１００点程度の特徴点を抽出する。
【0089】
続いて、カメラベクトル演算部４３で、抽出された特徴点の三次元座標が演算により求められ、その三次元座標に基づいてカメラベクトルが演算により求められる。具体的には、カメラベクトル演算部４３は、連続する各フレーム間に存在する、十分な数の特徴の位置と、移動するカメラ間の位置ベクトル、カメラの３軸回転ベクトル、各カメラ位置と特徴点をそれぞれ結んだベクトル等、各種三次元ベクトルの相対値を演算により連続的に算出する。
本実施形態では、例えば、３６０度全周画像のエピポーラ幾何からエピポーラ方程式を解くことによりカメラ運動（カメラ位置とカメラ回転）を計算するようになっている。また、これを三角測量として説明すれば、方位と俯角仰角を測量する機器を車載して、対象となる複数の測量地点を、移動する測量機器により同一対象物を計測して取得したデータから、各対象地点の座標と測量機器の移動の軌跡を求める場合と同じである。
本実施形態では、例えば、３６０度全周画像のエピポーラ幾何からエピポーラ方程式を解くことによりカメラ運動（カメラ位置とカメラ回転）を計算するようになっている。また、これを三角測量として説明すれば、方位と俯角仰角を測量する機器を車載して、対象となる複数の測量地点を、移動する測量機器により同一対象物を計測して取得したデータから、各対象地点の座標と測量機器の移動の軌跡を求める場合と同じである。
【0090】
図１９に示す画像１，２は、３６０度全周画像をメルカトール展開した画像であり、緯度φ、経度θとすると、画像１上の点は（θ１，φ１）、画像２上の点は（θ２，φ２）となる。そして、それぞれのカメラでの空間座標は、ｚ１＝（ｃｏｓφ１ｃｏｓθ１，ｃｏｓφ１ｓｉｎθ１，ｓｉｎφ１）、ｚ２＝（ｃｏｓφ２ｃｏｓθ２，ｃｏｓφ２ｓｉｎθ２，ｓｉｎφ２）である。カメラの移動ベクトルをｔ、カメラの回転行列をＲ、とすると、ｚ１［ｔ］×Ｒｚ２＝０がエピポーラ方程式である。
十分な数の特徴点を与えることにより、線形代数演算により最小自乗法による解としてｔ及びＲを計算することができる。この演算を対応する複数フレームに適用し演算する。
十分な数の特徴点を与えることにより、線形代数演算により最小自乗法による解としてｔ及びＲを計算することができる。この演算を対応する複数フレームに適用し演算する。
【0091】
ここで、カメラベクトルの演算に利用する画像としては、３６０度全周画像を用いることが好ましい。
カメラベクトル演算に用いる画像としては、原理的にはどのような画像でも良いが、全周映像は対象物の方向がそのまま緯度経度で表現されていることから、きわめて有利である。図１９に示す３６０度全周画像のような広角レンズによる画像の方が特徴点を数多く選択し、複数フレームに渡って追跡し易くなる。また、狭角レンズによる映像あっても、同一対象物の周りを移動しながら様々な方向から撮影する場合には特徴点が視界から逃げないので有効である。
カメラベクトル演算に用いる画像としては、原理的にはどのような画像でも良いが、全周映像は対象物の方向がそのまま緯度経度で表現されていることから、きわめて有利である。図１９に示す３６０度全周画像のような広角レンズによる画像の方が特徴点を数多く選択し、複数フレームに渡って追跡し易くなる。また、狭角レンズによる映像あっても、同一対象物の周りを移動しながら様々な方向から撮影する場合には特徴点が視界から逃げないので有効である。
【0092】
そこで、本実施形態では、ＣＶ演算に３６０度全周画像を用いており、これによって、特徴点の追跡距離を長くでき、特徴点を十分に多く選択することができ、遠距離、中距離、短距離それぞれに都合の良い特徴点を選択することができるようになる。また、回転ベクトルを補正する場合には、極回転変換処理を加えることで、演算処理も容易に行えるようになる。これらのことから、より精度の高い演算結果が得られるようになる。
なお、図１９は、ＣＶ演算部４０における処理を理解し易くするために、１台又は複数台のカメラで撮影した画像を合成した３６０度全周囲の球面画像を地図図法でいうメルカトール図法で展開したものを示しているが、実際のＣＶ演算では、必ずしもメルカトール図法による展開画像である必要はない。
なお、図１９は、ＣＶ演算部４０における処理を理解し易くするために、１台又は複数台のカメラで撮影した画像を合成した３６０度全周囲の球面画像を地図図法でいうメルカトール図法で展開したものを示しているが、実際のＣＶ演算では、必ずしもメルカトール図法による展開画像である必要はない。
【0093】
次に、誤差最小化部４４では、各フレームに対応する複数のカメラ位置と複数の特徴点の数により、複数通り生じる演算方程式により、各特徴点に基づくベクトルを複数通り演算して求めて、各特徴点の位置及びカメラ位置の分布が最小になるように統計処理をして、最終的なベクトルを求める。例えば、複数フレームのカメラ位置、カメラ回転及び複数の特徴点について、Levenberg-Marquardt法により最小自乗法の最適解を推定し、誤差を収束してカメラ位置、カメラ回転行列、特徴点の座標を求める。
さらに、誤差の分布が大きい特徴点については削除し、他の特徴点に基づいて再演算することで、各特徴点及びカメラ位置での演算の精度を上げるようにする。
このようにして、特徴点の位置とカメラベクトルを精度良く求めることができる。
さらに、誤差の分布が大きい特徴点については削除し、他の特徴点に基づいて再演算することで、各特徴点及びカメラ位置での演算の精度を上げるようにする。
このようにして、特徴点の位置とカメラベクトルを精度良く求めることができる。
【0094】
図２１～図２３に、ＣＶ演算により得られる特徴点の三次元座標とカメラベクトルの例を示す。図２１～図２３は、本実施形態のＣＶ演算によるベクトル検出方法を示す説明図であり、移動するカメラによって取得された複数のフレーム画像によって得られるカメラ及び対象物の相対的な位置関係を示す図である。
図２１では、図１９の画像１，２に示した特徴点１～４の三次元座標と、画像１と画像２の間で移動するカメラベクトル（Ｘ，Ｙ，Ｚ）が示されている。
図２２及び図２３は、充分に多くの特徴点とフレーム画像により得られた特徴点の位置と移動するカメラの位置が示されている。同図中、グラフ中央に直線状に連続する○印がカメラ位置であり、その周囲に位置する○印が特徴点の位置と高さを示している。
図２１では、図１９の画像１，２に示した特徴点１～４の三次元座標と、画像１と画像２の間で移動するカメラベクトル（Ｘ，Ｙ，Ｚ）が示されている。
図２２及び図２３は、充分に多くの特徴点とフレーム画像により得られた特徴点の位置と移動するカメラの位置が示されている。同図中、グラフ中央に直線状に連続する○印がカメラ位置であり、その周囲に位置する○印が特徴点の位置と高さを示している。
【0095】
ここで、ＣＶ演算部４０におけるＣＶ演算は、より高精度な特徴点とカメラ位置の三次元情報を高速に得るために、図２４に示すように、カメラから特徴点の距離に応じて複数の特徴点を設定し、複数の演算を繰り返し行うようにする。
具体的には、ＣＶ演算部４０では、画像内には映像的に特徴がある特徴点を自動検出し、各フレーム画像内に特徴点の対応点を求める際に、カメラベクトル演算に用いるｎ番目とｎ＋ｍ番目の二つのフレーム画像ＦｎとＦｎ＋ｍに着目して単位演算とし、ｎとｍを適切に設定した単位演算を繰り返すことができる。
ｍはフレーム間隔であり、カメラから画像内の特徴点までの距離によって特徴点を複数段に分類し、カメラから特徴点までの距離が遠いほどｍが大きくなるように設定し、カメラから特徴点までの距離が近いほどｍが小さくなるように設定する。このようにするのは、カメラから特徴点までの距離が遠ければ遠いほど、画像間における位置の変化が少ないからである。
具体的には、ＣＶ演算部４０では、画像内には映像的に特徴がある特徴点を自動検出し、各フレーム画像内に特徴点の対応点を求める際に、カメラベクトル演算に用いるｎ番目とｎ＋ｍ番目の二つのフレーム画像ＦｎとＦｎ＋ｍに着目して単位演算とし、ｎとｍを適切に設定した単位演算を繰り返すことができる。
ｍはフレーム間隔であり、カメラから画像内の特徴点までの距離によって特徴点を複数段に分類し、カメラから特徴点までの距離が遠いほどｍが大きくなるように設定し、カメラから特徴点までの距離が近いほどｍが小さくなるように設定する。このようにするのは、カメラから特徴点までの距離が遠ければ遠いほど、画像間における位置の変化が少ないからである。
【0096】
そして、特徴点のｍ値による分類を、十分にオーバーラップさせながら、複数段階のｍを設定し、画像の進行とともにｎが連続的に進行するのにともなって、演算を連続的に進行させる。そして、ｎの進行とｍの各段階で、同一特徴点について複数回重複演算を行う。
このようにして、フレーム画像ＦｎとＦｎ＋ｍに着目した単位演算を行うことにより、ｍ枚毎にサンプリングした各フレーム間（フレーム間は駒落ちしている）では、長時間かけて精密カメラベクトルを演算し、フレーム画像ＦｎとＦｎ＋ｍの間のｍ枚のフレーム（最小単位フレーム）では、短時間処理で行える簡易演算とすることができる。
このようにして、フレーム画像ＦｎとＦｎ＋ｍに着目した単位演算を行うことにより、ｍ枚毎にサンプリングした各フレーム間（フレーム間は駒落ちしている）では、長時間かけて精密カメラベクトルを演算し、フレーム画像ＦｎとＦｎ＋ｍの間のｍ枚のフレーム（最小単位フレーム）では、短時間処理で行える簡易演算とすることができる。
【0097】
ｍ枚毎の精密カメラベクトル演算に誤差がないとすれば、ｍ枚のフレームのカメラベクトルの両端は、高精度演算をしたＦｎとＦｎ＋ｍのカメラベクトルと重なることになる。従って、ＦｎとＦｎ＋ｍの中間のｍ枚の最小単位のフレームについては簡易演算で求め、簡易演算で求めたｍ枚の最小単位フレームのカメラベクトルの両端を、高精度演算で求めたＦｎとＦｎ＋ｍのカメラベクトルに一致するように、ｍ枚の連続したカメラベクトルのスケール調整をすることができる。
このようにして、画像の進行とともにｎが連続的に進行することにより、同一特徴点について複数回演算されて得られる各カメラベクトルの誤差が最小になるようにスケール調整して統合し、最終のカメラベクトルを決定することができる。これにより、誤差のない高精度のカメラベクトルを求めつつ、簡易演算を組み合わせることにより、演算処理を高速化することができるようになる。
このようにして、画像の進行とともにｎが連続的に進行することにより、同一特徴点について複数回演算されて得られる各カメラベクトルの誤差が最小になるようにスケール調整して統合し、最終のカメラベクトルを決定することができる。これにより、誤差のない高精度のカメラベクトルを求めつつ、簡易演算を組み合わせることにより、演算処理を高速化することができるようになる。
【0098】
ここで、簡易演算としては、精度に応じて種々の方法があるが、例えば、(1)高精度演算では１００個以上の多くの特徴点を用いる場合に、簡易演算では最低限の１０個程度の特徴点を用いる方法や、(2)同じ特徴点の数としても、特徴点とカメラ位置を同等に考えれば、そこには無数の三角形が成立し、その数だけの方程式が成立するため、その方程式の数を減らすことで、簡易演算とすることができる。
これによって、各特徴点及びカメラ位置の誤差が最小になるようにスケール調整する形で統合し、距離演算を行い、さらに、誤差の分布が大きい特徴点を削除し、必要に応じて他の特徴点について再演算することで、各特徴点及びカメラ位置での演算の精度を上げることができる。
これによって、各特徴点及びカメラ位置の誤差が最小になるようにスケール調整する形で統合し、距離演算を行い、さらに、誤差の分布が大きい特徴点を削除し、必要に応じて他の特徴点について再演算することで、各特徴点及びカメラ位置での演算の精度を上げることができる。
【0099】
また、このように高速な簡易演算を行うことにより、カメラベクトルのリアルタイム処理が可能となる。
カメラベクトルのリアルタイム処理は、目的の精度をとれる最低のフレーム数と、自動抽出した最低の特徴点数で演算を行い、カメラベクトルの概略値をリアルタイムで求め、表示し、次に、画像が蓄積するにつれて、フレーム数を増加させ、特徴点の数を増加させ、より精度の高いカメラベクトル演算を行い、概略値を精度の高いカメラベクトル値に置き換えて表示することができる。
カメラベクトルのリアルタイム処理は、目的の精度をとれる最低のフレーム数と、自動抽出した最低の特徴点数で演算を行い、カメラベクトルの概略値をリアルタイムで求め、表示し、次に、画像が蓄積するにつれて、フレーム数を増加させ、特徴点の数を増加させ、より精度の高いカメラベクトル演算を行い、概略値を精度の高いカメラベクトル値に置き換えて表示することができる。
【0100】
さらに、本実施形態では、より高精度のカメラベクトルを求めるために、三次元情報（三次元形状）の追跡を行うことができる。
具体的には、まず、三次元情報追跡部４５で、カメラベクトル演算部４３，誤差最小化部４４を経て得られたカメラベクトルを概略のカメラベクトルと位置づけ、その後のプロセスで生成される画像の一部として得られる三次元情報（三次元形状）に基づいて、複数のフレーム画像に含まれる部分的三次元情報を隣接するフレーム間で連続的に追跡して三次元形状の自動追跡を行う。
そして、この三次元情報追跡部４５で得られた三次元情報の追跡結果から、高精度カメラベクトル演算部４６においてより高精度なカメラベクトルが求められる。
具体的には、まず、三次元情報追跡部４５で、カメラベクトル演算部４３，誤差最小化部４４を経て得られたカメラベクトルを概略のカメラベクトルと位置づけ、その後のプロセスで生成される画像の一部として得られる三次元情報（三次元形状）に基づいて、複数のフレーム画像に含まれる部分的三次元情報を隣接するフレーム間で連続的に追跡して三次元形状の自動追跡を行う。
そして、この三次元情報追跡部４５で得られた三次元情報の追跡結果から、高精度カメラベクトル演算部４６においてより高精度なカメラベクトルが求められる。
【0101】
上述した特徴点抽出部４１及び特徴点対応処理部４２では、特徴点を複数のフレーム間画像内に自動追跡するが、特徴点が消失するなどして特徴点の追跡フレーム数に制限が出てくることがある。また、画像は二次元であり、追跡途中で形状が変化するために追跡精度にも一定の限界がある。
そこで、特徴点追跡で得られるカメラベクトルを概略値と位置づけ、その後のプロセスで得られる三次元情報（三次元形状）を各フレーム画像上に追跡して、その軌跡から高精度カメラベクトルを求めることができる。
三次元形状の追跡は、マッチング及び相関の精度を得やすく、三次元形状はフレーム画像によって、その三次元形状も大きさも変化しないので、多くのフレームに亘って追跡が可能であり、そのことでカメラベクトル演算の精度を向上させることができる。これはカメラベクトル演算部４３により概略のカメラベクトルが既知であり、三次元形状が既に分かっているから可能となるものである。
そこで、特徴点追跡で得られるカメラベクトルを概略値と位置づけ、その後のプロセスで得られる三次元情報（三次元形状）を各フレーム画像上に追跡して、その軌跡から高精度カメラベクトルを求めることができる。
三次元形状の追跡は、マッチング及び相関の精度を得やすく、三次元形状はフレーム画像によって、その三次元形状も大きさも変化しないので、多くのフレームに亘って追跡が可能であり、そのことでカメラベクトル演算の精度を向上させることができる。これはカメラベクトル演算部４３により概略のカメラベクトルが既知であり、三次元形状が既に分かっているから可能となるものである。
【0102】
カメラベクトルが概略値の場合、非常に多くのフレームに亘る三次元座標の誤差は、特徴点追跡による各フレームに関係するフレームが少ないので、誤差が累積して長距離では次第に大きな誤差になるが、画像の一部分を切り取ったときの三次元形状の誤差は相対的に少なく、形状の変化と大きさに及ぼす影響はかなり少ないものとなる。このため、三次元形状での比較や追跡は、二次元形状追跡の時よりも極めて有利となる。追跡において、二次元形状での追跡の場合、複数のフレームにおける形状の変化と大きさの変化を避けられないまま追跡することになるので、誤差が大きかったり、対応点が見つからないなどの問題があったが、三次元形状での追跡においては形状の変化が極めて少なく、しかも原理的に大きさの変化もないので、正確な追跡が可能となる。
【0103】
ここで、追跡の対象となる三次元形状データとしては、例えば、特徴点の三次元分布形状や、特徴点の三次元分布形状から求められるポリゴン面等がある。
また、得られた三次元形状を、カメラ位置から二次元画像に変換して、二次元画像として追跡することも可能である。カメラベクトルの概略値が既知であることから、カメラ視点からの二次元画像に投影変換が可能であり、カメラ視点の移動による対象の形状変化にも追従することが可能となる。
また、得られた三次元形状を、カメラ位置から二次元画像に変換して、二次元画像として追跡することも可能である。カメラベクトルの概略値が既知であることから、カメラ視点からの二次元画像に投影変換が可能であり、カメラ視点の移動による対象の形状変化にも追従することが可能となる。
【0104】
以上のようにして求められたカメラベクトルは、全方位カメラ３０で撮影されたビデオ映像中に重ねて表示することができる。
例えば、図２５に示すように、車載カメラからの映像を平面展開して、各フレーム画像内の目的平面上の対応点を自動で探索し、対応点を一致させるように結合して目的平面の結合画像を生成し、同一の座標系に統合して表示する。
さらに、その共通座標系の中にカメラ位置とカメラ方向を次々に検出し、その位置や方向、軌跡をプロットしていくことができる。ＣＶデータは、その三次元位置と３軸回転を示しており、ビデオ映像に重ねて表示することで、ビデオ映像の各フレームでＣＶ値を同時に観察できる。ＣＶデータをビデオ映像に重ねて表示した画像例を図２５に示す。
なお、ビデオ映像内にカメラ位置を正しく表示すると、ＣＶ値が示すビデオ映像内の位置は画像の中心となり、カメラ移動が直線に近い場合は、すべてのフレームのＣＶ値が重なって表示されてしまうので、例えば図２５に示すように、敢えてカメラ位置から真下に１メートルの位置を表示することが適切である。あるいは道路面までの距離を基準として、道路面の高さにＣＶ値を表示するのがより適切である。
例えば、図２５に示すように、車載カメラからの映像を平面展開して、各フレーム画像内の目的平面上の対応点を自動で探索し、対応点を一致させるように結合して目的平面の結合画像を生成し、同一の座標系に統合して表示する。
さらに、その共通座標系の中にカメラ位置とカメラ方向を次々に検出し、その位置や方向、軌跡をプロットしていくことができる。ＣＶデータは、その三次元位置と３軸回転を示しており、ビデオ映像に重ねて表示することで、ビデオ映像の各フレームでＣＶ値を同時に観察できる。ＣＶデータをビデオ映像に重ねて表示した画像例を図２５に示す。
なお、ビデオ映像内にカメラ位置を正しく表示すると、ＣＶ値が示すビデオ映像内の位置は画像の中心となり、カメラ移動が直線に近い場合は、すべてのフレームのＣＶ値が重なって表示されてしまうので、例えば図２５に示すように、敢えてカメラ位置から真下に１メートルの位置を表示することが適切である。あるいは道路面までの距離を基準として、道路面の高さにＣＶ値を表示するのがより適切である。
【0105】
以上説明したように、本発明の含水比マッピング方法及び含水比マッピング装置によれば、以下のような優れた作用効果を奏することができる。
まず、本発明によれば、マルチチャンネル型の地中レーダを使用することにより、広い範囲における含水比測定を容易・迅速かつ効率的に行うことができるようになる。
まず、本発明によれば、マルチチャンネル型の地中レーダを使用することにより、広い範囲における含水比測定を容易・迅速かつ効率的に行うことができるようになる。
【0106】
上述した特許文献１の測定方法を含めて、従来の地中探査技術では、地中探査データの位置精度が悪いため、仮にＣＭＰ法による測定を行ったとしても、中点位置の精度も悪いままであった。
また、手押し型のセンサ等を用いた従来技術では、広範囲の探査には非常に時間がかかり、また、牽引型の地中レーダを用いる場合には、得られる位置情報はＧＰＳによるものだけとなり、位置精度が良くなかった。特に、地中の水分動態を調べるには、ＧＰＳによる数メートルオーダーの誤差がある位置情報では有効な測定は行えず、数センチの誤差範囲での測定・探査は実現不可能であった。
また、従来の手押し型や牽引型の探査方法では、道路で行う場合には通行規制が必要となってしまい、また、規制をしない場合には、十分な位置精度を持たないために、含水比マッピングを実用的な位置精度を持った情報として得ることはできなかった。
また、手押し型のセンサ等を用いた従来技術では、広範囲の探査には非常に時間がかかり、また、牽引型の地中レーダを用いる場合には、得られる位置情報はＧＰＳによるものだけとなり、位置精度が良くなかった。特に、地中の水分動態を調べるには、ＧＰＳによる数メートルオーダーの誤差がある位置情報では有効な測定は行えず、数センチの誤差範囲での測定・探査は実現不可能であった。
また、従来の手押し型や牽引型の探査方法では、道路で行う場合には通行規制が必要となってしまい、また、規制をしない場合には、十分な位置精度を持たないために、含水比マッピングを実用的な位置精度を持った情報として得ることはできなかった。
【0107】
これに対して、本発明では、送信アンテナと受信アンテナが同一の中点と異なるオフセット距離を有する複数の送受信アンテナ組による測定を複数回繰り返すことにより、高精度の位置情報が付与されたＣＭＰ法による地中探査データを得ることができ、また、ＣＯＧ法による測定結果を組み合わせることができる。さらに、異なる日時の測定結果を比較することができる。
これによって、地中の土壌の含水比を広範囲にわたって正確かつ効率的に測定することができ、測定された含水比データに対して、ＧＰＳによることなく、より正確で高精度な位置情報を付与することができる。
これによって、地中の土壌の含水比を広範囲にわたって正確かつ効率的に測定することができ、測定された含水比データに対して、ＧＰＳによることなく、より正確で高精度な位置情報を付与することができる。
【0108】
また、ＣＶ演算に基づく高精度三次元座標が得られるため、上述した従来技術のような手押し型のセンサ等による位置特定工程が不要になり、地中データが一元化されたオルソ画像のみによって空洞・埋設管等の地中対象物を特定することが可能となる。
特に、従来技術では、例えば車両の前後左右を撮影する４方向カメラとＧＰＳによる簡易位置特定に基づく位置特定作業が必須であったことから、本実施形態による効果は非常に大きなものとなる。
また、このようにＣＶ演算により効率良く地上の高精度三次元座標が得られることから、検出位置特定のために、例えば現地での構造物の形状や道路周辺物等を計測する工程なども不要となり、作業効率は格段に向上することになる。
特に、従来技術では、例えば車両の前後左右を撮影する４方向カメラとＧＰＳによる簡易位置特定に基づく位置特定作業が必須であったことから、本実施形態による効果は非常に大きなものとなる。
また、このようにＣＶ演算により効率良く地上の高精度三次元座標が得られることから、検出位置特定のために、例えば現地での構造物の形状や道路周辺物等を計測する工程なども不要となり、作業効率は格段に向上することになる。
【0109】
また、本実施形態によれば、探査対象・範囲となる地下情報と地上情報について、全方位の情報を連続して取得することができ、地下・地上の全方位三次元情報を一元管理することができる。
また、このように一元管理された地下・地上情報のデータベースは、例えば米国ＥＳＲＩ社製「ＡｒｃＧＩＳ」（登録商標）など、ＧＩＳ（地理情報システム：Geographic Information System）関連のソフトウェア（アプリケーション）に対応させることができ、取扱性や汎用性，拡張性等にも優れたデータベース管理が可能となる。
また、このように一元管理された地下・地上情報のデータベースは、例えば米国ＥＳＲＩ社製「ＡｒｃＧＩＳ」（登録商標）など、ＧＩＳ（地理情報システム：Geographic Information System）関連のソフトウェア（アプリケーション）に対応させることができ、取扱性や汎用性，拡張性等にも優れたデータベース管理が可能となる。
【0110】
また、本実施形態で生成・出力される地下・地上の三次元情報が一元化されたオルソ画像は、時系列に従って蓄積・管理することができ、例えば将来の調査点の抽出や補修工事前後の比較と整合といった、地中探査情報の時系列管理を補正作業無しに容易かつ迅速に行うことができる。
また、本実施形態で得られるオルソ画像は、例えば１／５００の平面図（地図）と同じ精度を持ちながら、映像に基づいて視覚的に現地状況等を把握することでき、視覚的に地図等より優れ、単なる平面図（地図）では得られない有用な付加価値のある成果・情報を提供することができる。
さらに、本実施形態の地中探査装置によれば、探査車両を走行させるだけで、探査結果の全方位連続一元化処理及び映像から、簡単に高精度の三次元座標情報が得られることから、空洞等の地中探査だけでなく、例えば災害時の緊急輸送道路等の安全確保等、緊急時にも迅速な対応ができるという優れた効果もある。
また、本実施形態で得られるオルソ画像は、例えば１／５００の平面図（地図）と同じ精度を持ちながら、映像に基づいて視覚的に現地状況等を把握することでき、視覚的に地図等より優れ、単なる平面図（地図）では得られない有用な付加価値のある成果・情報を提供することができる。
さらに、本実施形態の地中探査装置によれば、探査車両を走行させるだけで、探査結果の全方位連続一元化処理及び映像から、簡単に高精度の三次元座標情報が得られることから、空洞等の地中探査だけでなく、例えば災害時の緊急輸送道路等の安全確保等、緊急時にも迅速な対応ができるという優れた効果もある。
【0111】
このように、本発明の一実施形態に係る地中探査装置によれば、地中探査用のセンサとなる地中レーダ２０を有効に保護しつつ、地中探査データとして高精度で正確なデータの取得を可能とし、また、地中及び地上の位置情報として、ＧＰＳでは得られない、より正確で高精度な位置情報を生成して、地中と地上の情報が一元化された高精度な三次元情報を得ることができる。
したがって、本実施形態によれば、地中探査データと地表面の情報を、従来技術のような地中探査データと平面図（地図等）を単に並べて表示させたものとは異なる、地中探査データと地上のオルソ画像を一元化した単一の三次元情報として生成・表示させることができる。これによって、空洞や埋設管等を示す地中探査データを、地上の映像を示す正確なオルソ画像中に一体的に表示させて特定・把握でき、従来技術では得られない正確・高精度な地中探査が可能となる。
したがって、本実施形態によれば、地中探査データと地表面の情報を、従来技術のような地中探査データと平面図（地図等）を単に並べて表示させたものとは異なる、地中探査データと地上のオルソ画像を一元化した単一の三次元情報として生成・表示させることができる。これによって、空洞や埋設管等を示す地中探査データを、地上の映像を示す正確なオルソ画像中に一体的に表示させて特定・把握でき、従来技術では得られない正確・高精度な地中探査が可能となる。
【0112】
なお、本発明の含水比マッピング方法及び含水比マッピング装置について、好ましい実施形態を示して説明したが、本発明は、上述した実施形態に限定されるものではなく、本発明の範囲で種々の変更実施が可能であることは言うまでもない。
例えば、上述した実施形態では、本発明に係る地中探査手段（センサ）として、電磁波を利用したレーダ（地中レーダ２０）を例に取って説明したが、地中探査手段は、地中の三次元情報を取得・生成できる限り、電磁波を用いるレーダに限定されるものではない。具体的には、地中探査手段としては、電磁波を用いる場合の他、地震波や電気比抵抗，超音波，レーダなどを用いるセンサ等であっても良い。
例えば、上述した実施形態では、本発明に係る地中探査手段（センサ）として、電磁波を利用したレーダ（地中レーダ２０）を例に取って説明したが、地中探査手段は、地中の三次元情報を取得・生成できる限り、電磁波を用いるレーダに限定されるものではない。具体的には、地中探査手段としては、電磁波を用いる場合の他、地震波や電気比抵抗，超音波，レーダなどを用いるセンサ等であっても良い。
【産業上の利用可能性】
【0113】
本発明は、地震波や電気比抵抗，電磁波，超音波等に対する地盤の物理的性質を利用して、地中（地下）の状況を非破壊で間接的に探査・調査する、地中レーダ探査等の地中探査技術に好適に利用可能である。
【符号の説明】
【0114】
１ 含水比マッピング装置
１０ 探査車両
１１ ＧＰＳ（高精度ＧＮＳＳ）
２０ 地中レーダ
２１ レーダ本体
２２ アンテナ
３０ 全方位カメラ
４０ ＣＶ演算部
５０ 一元化処理部
１０ 探査車両
１１ ＧＰＳ（高精度ＧＮＳＳ）
２０ 地中レーダ
２１ レーダ本体
２２ アンテナ
３０ 全方位カメラ
４０ ＣＶ演算部
５０ 一元化処理部
【図1】
【図2】
【図3】
【図4】
【図5】
【図6】
【図7】
【図8】
【図9】
【図10】
【図11】
【図12】
【図13】
【図14】
【図15】
【図16】
【図17】
【図18】
【図19】
【図20】
【図21】
【図22】
【図23】
【図24】
【図25】
