地中レーダ装置及び画像信号処理方法
【課題】地中内に埋設された物体の検知精度を向上した信頼性の高いものでありながら、その物体を特定するための演算式を簡素化することができる地中レーダ装置地中レーダ装置及び画像信号処理方法を提供する。
【解決手段】地表面を基準とする水平面内方向とこの地表面から地中に向かう深さ方向とで３次元の行列データの各要素の数値を取得する３次元データ取得ステップと、地表面を基準とする水平面内の移動平均を求める移動平均算出ステップと、前記３次元データ取得ステップと前記移動平均算出ステップとに基づいて３次元画像を作成する画像化ステップと、隣り合った画像の相互相関からそのピーク鋭度を求めるピーク鋭度算出ステップと、そのピーク鋭度が最大の条件から地雷画像を特定する地雷画像特定ステップと、その地雷画像の最大値をとる点に基づいて地雷座標を特定する地雷座標特定ステップとを備えている。
【解決手段】地表面を基準とする水平面内方向とこの地表面から地中に向かう深さ方向とで３次元の行列データの各要素の数値を取得する３次元データ取得ステップと、地表面を基準とする水平面内の移動平均を求める移動平均算出ステップと、前記３次元データ取得ステップと前記移動平均算出ステップとに基づいて３次元画像を作成する画像化ステップと、隣り合った画像の相互相関からそのピーク鋭度を求めるピーク鋭度算出ステップと、そのピーク鋭度が最大の条件から地雷画像を特定する地雷画像特定ステップと、その地雷画像の最大値をとる点に基づいて地雷座標を特定する地雷座標特定ステップとを備えている。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、地中内の埋設物を３次元の立体画像として可視化する地中レーダ装置及び画像信号処理方法に関する。
【背景技術】
【０００２】
【特許文献１】特開平０９−０３３１９４号公報
【特許文献２】特開２００４−１９８１９５号公報
【０００３】
従来から、対戦車地雷や対人地雷に代表される地雷は、世界中の紛争地域で使用され、現在においても、無数の地雷が埋設されたままとなっており、紛争後におけるこれらの除去には困難を極めている。
【０００４】
紛争後の地雷除去、所謂、人道的地雷除去においては、地雷原となっている土地での生活を前提としていることから、その地雷原に埋設された地雷のすべてを除去しなければならないが、これらの地雷によっては、その探知及び除去が困難となっているのが実情である。
【０００５】
即ち、地雷の種類として、例えば、地雷探知器の磁場に反応して起爆するものや数回の圧力付加によって起爆する種類のもの等があり、また、金属・非金属によって構成されているもの等がある。
【０００６】
これにより、磁場に反応して起爆してしまう地雷の場合、金属探知器を近付けただけで起爆する可能性がある。また、数回の圧力付加によって起爆する地雷の場合、地雷処理車が通過したときには起爆せず、これに続いて通過する地雷処理車以外の車両等が通過したときに起爆してしまう可能性がある。
【０００７】
しかも、地中内には、既に起爆してしまった地雷の破片や、土中の水分等の反射体、或いは、金属成分を多く含む土質といったように、レーダの種類等によっては誤認識してしまう虞が多分にあった。
【０００８】
一方、このような地雷を探知する装置として、特許文献１に開示のような金属探知機を利用した装置や、特許文献２に開示のように、電磁波を地中に伝播させ、その地中からの反射波をセンサにより検出すると共に、その検出結果に基づいて画像化する装置が知られている。
【発明の開示】
【発明が解決しようとする課題】
【０００９】
ところで、上記の如く構成された地中レーダ装置及び画像信号処理方法にあっては、上述した地雷の種類や構成材料等の様々な要因から、特許文献１に開示のような金属探知機では、その誤認識のほうが多く、撤去作業に時間を要するうえ、埋設位置（深さを含む）の詳細な位置を特定するには熟練を要するといった問題も生じていた。
【００１０】
また、特許文献２に開示の装置の場合、埋設物の３次元位置が視覚的に判り易いように、地中の状態を水平断面視及び垂直断面視で表現しているが、埋設物自体を３次元の立体画像として表示するものではないため、埋設物の形状を容易に把握することが困難で、地雷であるか否かの特定が困難であるという問題が生じていた。
【００１１】
さらに、このような３次元の画像化を可能とする場合、一般的な画像解析の手法を適用した場合、埋設物の有無を特定するには信頼性が低いという問題が生じていた。
【００１２】
尚、埋設物の一つとして、地雷等の特殊で且つある程度の形状が予め特定できる埋設物の場合、予めパターン化された３次元の立体画像と検出された３次元の物体画像パターンとを比較することで地雷であるか否かを特定するパターンマッチング方式なども考えられているが、これらのパターンマッチング化は、例えば、上述した金属成分等を含む土中では、地雷の周囲や地雷に付着した金属成分を含む土までも埋設物として認識してパターン化してしまい、両者間でマッチングしない虞があるなど、信頼性が低いという問題が生じていた。
【００１３】
一方、地中レーダ（Ｇｒｏｕｎｄ ｐｅｎｅｔｒａｔｉｎｇ Ｒａｄａｒ：ＧＰＲ）は、アンテナから電波を地中に放射し、地中の埋設物や地層境界面などから反射を受けた電波を受信することによって地中を可視化する装置である。
【００１４】
この地中レーダによって受信したＧＰＲデータは、地表面を基準とするＸ方向とＹ方向の水平面内の２次元行列と、地表面から地中に向かうＺ方向（深さ方向）とを含む３次元行列で構成されており、このＧＰＲデータに基づいて、図１に示すように、３次元の立体画像を算出することができる。
【００１５】
尚、図１においては、地中の埋設物としての地雷Ｊの存在並びに位置を明確に確認することができるが、多くの場合、地中に含まれる砂礫等の不均質によって不要反射となるクラッタＣが発生し、単純な３次元画像では、地雷ＪとクラッタＣとの区別をするのは困難である。
【００１６】
ところで、上述したＧＰＲデータとしての各深度毎の水平画像用データには、確実に反射体の情報が含まれていることから、地雷の位置を明確に判断することは可能である。しかしながら、これが３次元の膨大なＧＰＲデータとなると、クラッタＣを排除しながら対象物を見いだすことは非常に難しい。
【００１７】
そこで、このような３次元ＧＰＲデータ中に含まれる埋設物からの反射波を効率的に検知し、埋設物の位置を正確且つ簡素に検出することができれば、地雷除去作業現場での迅速な作業並びに装置の小型化を実現することができる。
【００１８】
本発明は、上記問題を解決するため、地中内に埋設された物体の検知精度を向上し得て、信頼性を向上させることができる地中レーダ装置及び画像信号処理方法を提供することを目的とする。
【課題を解決するための手段】
【００１９】
その目的を達成するため、本発明の画像信号処理方法は、地表面を基準とする水平面内方向とこの地表面から地中に向かう深さ方向とで３次元の行列データの各要素の数値を取得する３次元データ取得ステップと、地表面を基準とする水平面内の移動平均を求める移動平均算出ステップと、前記３次元データ取得ステップと前記移動平均算出ステップとに基づいて３次元画像を作成する画像化ステップと、隣り合った画像の相互相関からそのピーク鋭度を求めるピーク鋭度算出ステップと、そのピーク鋭度が最大の条件から埋設物の有無を特定する埋設物特定ステップとを備えていることを特徴とする。
【００２０】
即ち、本発明の画像信号処理方法によれば、１次元或いは２次元の走査線に沿って取得された２次元（１次元位置−深度）或いは３次元（２次元位置−深度）のデータ空間において、特定の深度方向に隣接する範囲における空間相関係数を計算し、空間相関係数が大きい領域では地中レーダデータが連続したパターンを有することから、目標物からの反射波を捉えて目標物の位置を検知すると共に、その目標物からの反射波の連続性から目標物（埋設物）の検知率を高めることができる。
【発明の効果】
【００２１】
本発明の地中レーダ装置及び画像信号処理方法によれば、埋設物の有無に関する信頼性を向上させることができる。
【発明を実施するための最良の形態】
【００２２】
次に、本発明の地中レーダ装置及び画像信号処理方法に係る実施の形態を図面に基づいて説明する。
【００２３】
尚、本発明の地中レーダ装置は、ＧＰＲレーダとこのＧＰＲレーダを搭載した地雷除去車両等の作業車と、ＧＰＲレーダからの検知信号（ＧＰＲデータ）を処理するために作業車に搭載された演算部としてのパーソナルコンピュータとを備えている。尚、このパーソナルコンピュータは作業車に搭載された無線通信装置を介して遠隔地でＧＰＲレーダからの検知信号（ＧＰＲデータ）を受信・処理するようにしても良い。
【００２４】
以下、本発明における地中レーダ装置のパーソナルコンピュータ上での処理方法を説明する。
【００２５】
（地雷画像抽出アルゴリズム）
〔２次元画像化〕
ＧＰＲ測定データは、Ｘ方向，Ｙ方向，Ｚ方向の３次元行列で構成されている。その行列のＸ方向とＹ方向とで構成する平面を画像化することで、Ｚ方向（深さ方向）における水平面画像が得られる。地雷などの埋設物画像は、深度方向に連続して現れ、いったん消えてもまた色が反転して現れる。これは地雷からの反射波の値が正から負に移り変わることに起因する。これに対し、クラッタの場合は、細かく不規則に変化する。
【００２６】
地雷が存在するＸ−Ｙ座標点の値の、深さＺにおけるＧＰＲ波形変化を図２に示す。図２において、深さ０．２ｍ付近の２本の二点鎖線で挟まれている部分が地雷からの反射波の値である。
【００２７】
〔移動平均画像〕
図３（Ａ），（Ｂ）は、深度を変えた場合の任意の深度毎のＧＰＲ水平面画像の例を示し、図３（Ａ）は地雷の場合、図３（Ｂ）はクラッタの場合を示す。尚、実際の画像ではそのパワーが強いほど赤くなるサーモセンサーのようなカラー処理された画像となっている。
【００２８】
図３（Ａ）に示すように、地雷からの反射は、値の大きさは深度に対して連続的に変化するのに対し、図３（Ｂ）に示すように、クラッタからの反射は不規則に変化する。波形のパワーの移動平均をとることで地雷からの反射波は安定した変化となり、クラッタとの識別が容易になる。今、図４に示すように、水平画像の深さ方向に対する深度変化の平均幅をＷ、ずらし幅をＤｚとしてＸ−Ｙ平面の移動平均を計算する。例として、図５において、Ｗ＝５０（４ｃｍ），Ｄｚ＝２５（２ｃｍ）としたときに得られる画像の一部を示す。また、この図５（Ｃ）に示した画像は地雷画像である。
【００２９】
図５（Ｃ）において、地雷は移動平均をとることによってクラッタと明確に区別されていることが分かる。また、連続して現れていることが確認できる。また、図５（Ｃ）より、地雷画像において最大値をとる座標点が地雷位置であることが分かる。これらの特徴から、移動平均画像において隣り合った画像同士の類似度が高ければ、それらが地雷画像であると特定することができる。さらに、その最大値をとる座標点が地雷位置である。
【００３０】
〔空間相関法を用いた類似度の評価〕
画像同士の類似度を評価するために空間的な相互相関関数を利用する。まずＧＰＲ信号を数１とする。
【００３１】
【数１】
この数１のそれぞれを、2次元フーリエ変換することにより、数２に示すように、フーリエスペクトルが得られる。
【００３２】
【数２】
この数２より、数３に示すクロス・パワースペクトルを求めることができる。
【００３３】
【数３】
さらに、この数３を逆フーリエ変換することにより、数４に示す相互相関関数を得ることができる。
【００３４】
【数４】
【００３５】
隣り合う移動平均画像の相互相関を求めた例を図６及び図７に示す。図６（Ａ），（Ｂ）に示したクラッタ同士に基づく相互相関画像（図６（Ｃ））は中心になだらかなピークを持つが、図７（Ａ），（Ｂ）に示した地雷画像同士に基づく相互相関画像（図７（Ｃ））は中心に鋭いピークを持つことがわかる。つまり２つの画像の類似度が高いほど、この相互相関画像のピークは鋭くなる。
【００３６】
こうした類似度を定量化するために、図８，図７（Ｃ）に示すように、Ｘ，Ｙ方向についてピークの半値幅δｘ，δｙを求め、そのピークの鋭度＝δｘ・δｙとして評価する。その値の深度毎の変化をグラフ化した例を図９に示す。
【００３７】
この図９において、ピーク鋭度は深さ０．２ｍ付近で最小値をとっている。このピーク鋭度が小さい場合、連続する反射波の連続性が強い、つまり反射体が存在する可能性が高い。この相関が得られた移動平均画像が地雷画像であり、その画像の最大値をとる座標が地雷位置である。
【００３８】
〔地雷位置特定のアルゴリズムの構築〕
地雷位置を特定するアルゴリズムをまとめると、３次元行列データの各要素のパワーをとる・Ｘ−Ｙ平面の移動平均を求めてそれぞれを画像化する・隣り合った画像の相互相関を求めてそのピーク鋭度を求める・ピーク鋭度が最大の条件から地雷画像を特定する・地雷画像において最大値をとる点を求める・地雷座標を特定する。
【００３９】
〔適用例〕
ここに示したアルゴリズムにおいて、平均をとる幅Ｗと、それをずらす幅Ｄｚの選択は非常に重要である。地雷画像の特徴を抽出するためには、適切な平均幅Ｗとずらし幅Ｄｚを設定しなければならない。次に発明したアルゴリズムを適用した具体例を挙げる。
【００４０】
ここで扱うＧＰＲデータは、乾いた地質において測定されたものである。また、ターゲットはＴＹＰＥ７２型模擬地雷（直径７．６ｃｍ、厚さ４ｃｍ）である。図１０，図１１に示すように、地雷位置を特定することができる。この際、平均幅Ｗ＝３０（２．４ｃｍ），ずらし幅Ｄｚ＝１５（１．２ｃｍ）である。
【００４１】
ところで、以上の説明では埋設物検知として地雷検知を例として使用したが、特に地雷に限定されるものではない。また、地雷のように小さく形が定常ではなく、土壌の地下水分布のような非定常形状の場合についてもＧＰＲ信号の深度方向への連続性がある場合、例えば、地中石油系汚染物質を対象とするなど、異常（異形）物体の検知も可能である。
【図面の簡単な説明】
【００４２】
【図１】本発明の一実施形態を示し、３次元立体画像の模式図である。
【図２】本発明の一実施形態を示し、地雷位置のＸ−Ｙ座標点の値の変化のグラフ図である。
【図３】本発明の一実施形態を示し、（Ａ）は地雷の場合の水平面画像の深さによる変化の説明図、（Ｂ）はクラッタの場合の水平面画像の深さによる変化の説明図である。
【図４】本発明の一実施形態を示し、深さ方向における移動平均の概念の説明図である。
【図５】本発明の一実施形態を示し、（Ａ）〜（Ｅ）は移動平均画像の説明図である。
【図６】本発明の一実施形態を示し、（Ａ）はクラッタの場合の移動平均画像の相互相関の一方となる画像の説明図、（Ｂ）はクラッタの場合の移動平均画像の相互相関の他方となる画像の説明図、（Ｃ）はクラッタの場合の相互相関像の説明図である。
【図７】本発明の一実施形態を示し、（Ａ）は地雷の場合の移動平均画像の相互相関の一方となる画像の説明図、（Ｂ）は地雷の場合の移動平均画像の相互相関の他方となる画像の説明図、（Ｃ）は地雷の場合の相互相関像の説明図である。
【図８】本発明の一実施形態を示し、ピーク断面のグラフ図である。
【図９】本発明の一実施形態を示し、ピーク鋭度−深さの関係のグラフ図である。
【図１０】本発明の一実施形態を示し、地雷特定の説明図である。
【図１１】本発明の一実施形態を示し、地雷特定のためのピーク鋭度−深さの関係のグラフ図である。
【技術分野】
【０００１】
本発明は、地中内の埋設物を３次元の立体画像として可視化する地中レーダ装置及び画像信号処理方法に関する。
【背景技術】
【０００２】
【特許文献１】特開平０９−０３３１９４号公報
【特許文献２】特開２００４−１９８１９５号公報
【０００３】
従来から、対戦車地雷や対人地雷に代表される地雷は、世界中の紛争地域で使用され、現在においても、無数の地雷が埋設されたままとなっており、紛争後におけるこれらの除去には困難を極めている。
【０００４】
紛争後の地雷除去、所謂、人道的地雷除去においては、地雷原となっている土地での生活を前提としていることから、その地雷原に埋設された地雷のすべてを除去しなければならないが、これらの地雷によっては、その探知及び除去が困難となっているのが実情である。
【０００５】
即ち、地雷の種類として、例えば、地雷探知器の磁場に反応して起爆するものや数回の圧力付加によって起爆する種類のもの等があり、また、金属・非金属によって構成されているもの等がある。
【０００６】
これにより、磁場に反応して起爆してしまう地雷の場合、金属探知器を近付けただけで起爆する可能性がある。また、数回の圧力付加によって起爆する地雷の場合、地雷処理車が通過したときには起爆せず、これに続いて通過する地雷処理車以外の車両等が通過したときに起爆してしまう可能性がある。
【０００７】
しかも、地中内には、既に起爆してしまった地雷の破片や、土中の水分等の反射体、或いは、金属成分を多く含む土質といったように、レーダの種類等によっては誤認識してしまう虞が多分にあった。
【０００８】
一方、このような地雷を探知する装置として、特許文献１に開示のような金属探知機を利用した装置や、特許文献２に開示のように、電磁波を地中に伝播させ、その地中からの反射波をセンサにより検出すると共に、その検出結果に基づいて画像化する装置が知られている。
【発明の開示】
【発明が解決しようとする課題】
【０００９】
ところで、上記の如く構成された地中レーダ装置及び画像信号処理方法にあっては、上述した地雷の種類や構成材料等の様々な要因から、特許文献１に開示のような金属探知機では、その誤認識のほうが多く、撤去作業に時間を要するうえ、埋設位置（深さを含む）の詳細な位置を特定するには熟練を要するといった問題も生じていた。
【００１０】
また、特許文献２に開示の装置の場合、埋設物の３次元位置が視覚的に判り易いように、地中の状態を水平断面視及び垂直断面視で表現しているが、埋設物自体を３次元の立体画像として表示するものではないため、埋設物の形状を容易に把握することが困難で、地雷であるか否かの特定が困難であるという問題が生じていた。
【００１１】
さらに、このような３次元の画像化を可能とする場合、一般的な画像解析の手法を適用した場合、埋設物の有無を特定するには信頼性が低いという問題が生じていた。
【００１２】
尚、埋設物の一つとして、地雷等の特殊で且つある程度の形状が予め特定できる埋設物の場合、予めパターン化された３次元の立体画像と検出された３次元の物体画像パターンとを比較することで地雷であるか否かを特定するパターンマッチング方式なども考えられているが、これらのパターンマッチング化は、例えば、上述した金属成分等を含む土中では、地雷の周囲や地雷に付着した金属成分を含む土までも埋設物として認識してパターン化してしまい、両者間でマッチングしない虞があるなど、信頼性が低いという問題が生じていた。
【００１３】
一方、地中レーダ（Ｇｒｏｕｎｄ ｐｅｎｅｔｒａｔｉｎｇ Ｒａｄａｒ：ＧＰＲ）は、アンテナから電波を地中に放射し、地中の埋設物や地層境界面などから反射を受けた電波を受信することによって地中を可視化する装置である。
【００１４】
この地中レーダによって受信したＧＰＲデータは、地表面を基準とするＸ方向とＹ方向の水平面内の２次元行列と、地表面から地中に向かうＺ方向（深さ方向）とを含む３次元行列で構成されており、このＧＰＲデータに基づいて、図１に示すように、３次元の立体画像を算出することができる。
【００１５】
尚、図１においては、地中の埋設物としての地雷Ｊの存在並びに位置を明確に確認することができるが、多くの場合、地中に含まれる砂礫等の不均質によって不要反射となるクラッタＣが発生し、単純な３次元画像では、地雷ＪとクラッタＣとの区別をするのは困難である。
【００１６】
ところで、上述したＧＰＲデータとしての各深度毎の水平画像用データには、確実に反射体の情報が含まれていることから、地雷の位置を明確に判断することは可能である。しかしながら、これが３次元の膨大なＧＰＲデータとなると、クラッタＣを排除しながら対象物を見いだすことは非常に難しい。
【００１７】
そこで、このような３次元ＧＰＲデータ中に含まれる埋設物からの反射波を効率的に検知し、埋設物の位置を正確且つ簡素に検出することができれば、地雷除去作業現場での迅速な作業並びに装置の小型化を実現することができる。
【００１８】
本発明は、上記問題を解決するため、地中内に埋設された物体の検知精度を向上し得て、信頼性を向上させることができる地中レーダ装置及び画像信号処理方法を提供することを目的とする。
【課題を解決するための手段】
【００１９】
その目的を達成するため、本発明の画像信号処理方法は、地表面を基準とする水平面内方向とこの地表面から地中に向かう深さ方向とで３次元の行列データの各要素の数値を取得する３次元データ取得ステップと、地表面を基準とする水平面内の移動平均を求める移動平均算出ステップと、前記３次元データ取得ステップと前記移動平均算出ステップとに基づいて３次元画像を作成する画像化ステップと、隣り合った画像の相互相関からそのピーク鋭度を求めるピーク鋭度算出ステップと、そのピーク鋭度が最大の条件から埋設物の有無を特定する埋設物特定ステップとを備えていることを特徴とする。
【００２０】
即ち、本発明の画像信号処理方法によれば、１次元或いは２次元の走査線に沿って取得された２次元（１次元位置−深度）或いは３次元（２次元位置−深度）のデータ空間において、特定の深度方向に隣接する範囲における空間相関係数を計算し、空間相関係数が大きい領域では地中レーダデータが連続したパターンを有することから、目標物からの反射波を捉えて目標物の位置を検知すると共に、その目標物からの反射波の連続性から目標物（埋設物）の検知率を高めることができる。
【発明の効果】
【００２１】
本発明の地中レーダ装置及び画像信号処理方法によれば、埋設物の有無に関する信頼性を向上させることができる。
【発明を実施するための最良の形態】
【００２２】
次に、本発明の地中レーダ装置及び画像信号処理方法に係る実施の形態を図面に基づいて説明する。
【００２３】
尚、本発明の地中レーダ装置は、ＧＰＲレーダとこのＧＰＲレーダを搭載した地雷除去車両等の作業車と、ＧＰＲレーダからの検知信号（ＧＰＲデータ）を処理するために作業車に搭載された演算部としてのパーソナルコンピュータとを備えている。尚、このパーソナルコンピュータは作業車に搭載された無線通信装置を介して遠隔地でＧＰＲレーダからの検知信号（ＧＰＲデータ）を受信・処理するようにしても良い。
【００２４】
以下、本発明における地中レーダ装置のパーソナルコンピュータ上での処理方法を説明する。
【００２５】
（地雷画像抽出アルゴリズム）
〔２次元画像化〕
ＧＰＲ測定データは、Ｘ方向，Ｙ方向，Ｚ方向の３次元行列で構成されている。その行列のＸ方向とＹ方向とで構成する平面を画像化することで、Ｚ方向（深さ方向）における水平面画像が得られる。地雷などの埋設物画像は、深度方向に連続して現れ、いったん消えてもまた色が反転して現れる。これは地雷からの反射波の値が正から負に移り変わることに起因する。これに対し、クラッタの場合は、細かく不規則に変化する。
【００２６】
地雷が存在するＸ−Ｙ座標点の値の、深さＺにおけるＧＰＲ波形変化を図２に示す。図２において、深さ０．２ｍ付近の２本の二点鎖線で挟まれている部分が地雷からの反射波の値である。
【００２７】
〔移動平均画像〕
図３（Ａ），（Ｂ）は、深度を変えた場合の任意の深度毎のＧＰＲ水平面画像の例を示し、図３（Ａ）は地雷の場合、図３（Ｂ）はクラッタの場合を示す。尚、実際の画像ではそのパワーが強いほど赤くなるサーモセンサーのようなカラー処理された画像となっている。
【００２８】
図３（Ａ）に示すように、地雷からの反射は、値の大きさは深度に対して連続的に変化するのに対し、図３（Ｂ）に示すように、クラッタからの反射は不規則に変化する。波形のパワーの移動平均をとることで地雷からの反射波は安定した変化となり、クラッタとの識別が容易になる。今、図４に示すように、水平画像の深さ方向に対する深度変化の平均幅をＷ、ずらし幅をＤｚとしてＸ−Ｙ平面の移動平均を計算する。例として、図５において、Ｗ＝５０（４ｃｍ），Ｄｚ＝２５（２ｃｍ）としたときに得られる画像の一部を示す。また、この図５（Ｃ）に示した画像は地雷画像である。
【００２９】
図５（Ｃ）において、地雷は移動平均をとることによってクラッタと明確に区別されていることが分かる。また、連続して現れていることが確認できる。また、図５（Ｃ）より、地雷画像において最大値をとる座標点が地雷位置であることが分かる。これらの特徴から、移動平均画像において隣り合った画像同士の類似度が高ければ、それらが地雷画像であると特定することができる。さらに、その最大値をとる座標点が地雷位置である。
【００３０】
〔空間相関法を用いた類似度の評価〕
画像同士の類似度を評価するために空間的な相互相関関数を利用する。まずＧＰＲ信号を数１とする。
【００３１】
【数１】
この数１のそれぞれを、2次元フーリエ変換することにより、数２に示すように、フーリエスペクトルが得られる。
【００３２】
【数２】
この数２より、数３に示すクロス・パワースペクトルを求めることができる。
【００３３】
【数３】
さらに、この数３を逆フーリエ変換することにより、数４に示す相互相関関数を得ることができる。
【００３４】
【数４】
【００３５】
隣り合う移動平均画像の相互相関を求めた例を図６及び図７に示す。図６（Ａ），（Ｂ）に示したクラッタ同士に基づく相互相関画像（図６（Ｃ））は中心になだらかなピークを持つが、図７（Ａ），（Ｂ）に示した地雷画像同士に基づく相互相関画像（図７（Ｃ））は中心に鋭いピークを持つことがわかる。つまり２つの画像の類似度が高いほど、この相互相関画像のピークは鋭くなる。
【００３６】
こうした類似度を定量化するために、図８，図７（Ｃ）に示すように、Ｘ，Ｙ方向についてピークの半値幅δｘ，δｙを求め、そのピークの鋭度＝δｘ・δｙとして評価する。その値の深度毎の変化をグラフ化した例を図９に示す。
【００３７】
この図９において、ピーク鋭度は深さ０．２ｍ付近で最小値をとっている。このピーク鋭度が小さい場合、連続する反射波の連続性が強い、つまり反射体が存在する可能性が高い。この相関が得られた移動平均画像が地雷画像であり、その画像の最大値をとる座標が地雷位置である。
【００３８】
〔地雷位置特定のアルゴリズムの構築〕
地雷位置を特定するアルゴリズムをまとめると、３次元行列データの各要素のパワーをとる・Ｘ−Ｙ平面の移動平均を求めてそれぞれを画像化する・隣り合った画像の相互相関を求めてそのピーク鋭度を求める・ピーク鋭度が最大の条件から地雷画像を特定する・地雷画像において最大値をとる点を求める・地雷座標を特定する。
【００３９】
〔適用例〕
ここに示したアルゴリズムにおいて、平均をとる幅Ｗと、それをずらす幅Ｄｚの選択は非常に重要である。地雷画像の特徴を抽出するためには、適切な平均幅Ｗとずらし幅Ｄｚを設定しなければならない。次に発明したアルゴリズムを適用した具体例を挙げる。
【００４０】
ここで扱うＧＰＲデータは、乾いた地質において測定されたものである。また、ターゲットはＴＹＰＥ７２型模擬地雷（直径７．６ｃｍ、厚さ４ｃｍ）である。図１０，図１１に示すように、地雷位置を特定することができる。この際、平均幅Ｗ＝３０（２．４ｃｍ），ずらし幅Ｄｚ＝１５（１．２ｃｍ）である。
【００４１】
ところで、以上の説明では埋設物検知として地雷検知を例として使用したが、特に地雷に限定されるものではない。また、地雷のように小さく形が定常ではなく、土壌の地下水分布のような非定常形状の場合についてもＧＰＲ信号の深度方向への連続性がある場合、例えば、地中石油系汚染物質を対象とするなど、異常（異形）物体の検知も可能である。
【図面の簡単な説明】
【００４２】
【図１】本発明の一実施形態を示し、３次元立体画像の模式図である。
【図２】本発明の一実施形態を示し、地雷位置のＸ−Ｙ座標点の値の変化のグラフ図である。
【図３】本発明の一実施形態を示し、（Ａ）は地雷の場合の水平面画像の深さによる変化の説明図、（Ｂ）はクラッタの場合の水平面画像の深さによる変化の説明図である。
【図４】本発明の一実施形態を示し、深さ方向における移動平均の概念の説明図である。
【図５】本発明の一実施形態を示し、（Ａ）〜（Ｅ）は移動平均画像の説明図である。
【図６】本発明の一実施形態を示し、（Ａ）はクラッタの場合の移動平均画像の相互相関の一方となる画像の説明図、（Ｂ）はクラッタの場合の移動平均画像の相互相関の他方となる画像の説明図、（Ｃ）はクラッタの場合の相互相関像の説明図である。
【図７】本発明の一実施形態を示し、（Ａ）は地雷の場合の移動平均画像の相互相関の一方となる画像の説明図、（Ｂ）は地雷の場合の移動平均画像の相互相関の他方となる画像の説明図、（Ｃ）は地雷の場合の相互相関像の説明図である。
【図８】本発明の一実施形態を示し、ピーク断面のグラフ図である。
【図９】本発明の一実施形態を示し、ピーク鋭度−深さの関係のグラフ図である。
【図１０】本発明の一実施形態を示し、地雷特定の説明図である。
【図１１】本発明の一実施形態を示し、地雷特定のためのピーク鋭度−深さの関係のグラフ図である。
【特許請求の範囲】
【請求項１】
地表面に向かってレーダ波を照射すると共に地中内で反射した反射波を受信するＧＰＲレーダと、該ＧＰＲレーダで受信した反射波に基づいて任意の深度毎の水平画像データを求めると共にその求めた水平画像データの相互関係からピーク鋭度を求める演算部とを備えていることを特徴とする地中レーダ装置。
【請求項２】
地表面を基準とする水平面内方向とこの地表面から地中に向かう深さ方向とで３次元の行列データの各要素の数値を取得する３次元データ取得ステップと、地表面を基準とする水平面内の移動平均を求める移動平均算出ステップと、前記３次元データ取得ステップと前記移動平均算出ステップとに基づいて３次元画像を作成する画像化ステップと、隣り合った画像の相互相関からそのピーク鋭度を求めるピーク鋭度算出ステップと、そのピーク鋭度が最大の条件から埋設物の有無を特定する埋設物特定ステップとを備えていることを特徴とする画像信号処理方法。
【請求項３】
地表面を基準とする水平面内方向とこの地表面から地中に向かう深さ方向とで３次元の行列データの各要素の数値を取得する３次元データ取得ステップと、地表面を基準とする水平面内の移動平均を求める移動平均算出ステップと、前記３次元データ取得ステップと前記移動平均算出ステップとに基づいて３次元画像を作成する画像化ステップと、隣り合った画像の相互相関からそのピーク鋭度を求めるピーク鋭度算出ステップと、そのピーク鋭度が最大の条件から埋設物の有無を特定する埋設物特定ステップと、その埋設物画像の最大値をとる点に基づいて埋設物座標を特定する埋設物座標特定ステップとを備えていることを特徴とする画像信号処理方法。
【請求項１】
地表面に向かってレーダ波を照射すると共に地中内で反射した反射波を受信するＧＰＲレーダと、該ＧＰＲレーダで受信した反射波に基づいて任意の深度毎の水平画像データを求めると共にその求めた水平画像データの相互関係からピーク鋭度を求める演算部とを備えていることを特徴とする地中レーダ装置。
【請求項２】
地表面を基準とする水平面内方向とこの地表面から地中に向かう深さ方向とで３次元の行列データの各要素の数値を取得する３次元データ取得ステップと、地表面を基準とする水平面内の移動平均を求める移動平均算出ステップと、前記３次元データ取得ステップと前記移動平均算出ステップとに基づいて３次元画像を作成する画像化ステップと、隣り合った画像の相互相関からそのピーク鋭度を求めるピーク鋭度算出ステップと、そのピーク鋭度が最大の条件から埋設物の有無を特定する埋設物特定ステップとを備えていることを特徴とする画像信号処理方法。
【請求項３】
地表面を基準とする水平面内方向とこの地表面から地中に向かう深さ方向とで３次元の行列データの各要素の数値を取得する３次元データ取得ステップと、地表面を基準とする水平面内の移動平均を求める移動平均算出ステップと、前記３次元データ取得ステップと前記移動平均算出ステップとに基づいて３次元画像を作成する画像化ステップと、隣り合った画像の相互相関からそのピーク鋭度を求めるピーク鋭度算出ステップと、そのピーク鋭度が最大の条件から埋設物の有無を特定する埋設物特定ステップと、その埋設物画像の最大値をとる点に基づいて埋設物座標を特定する埋設物座標特定ステップとを備えていることを特徴とする画像信号処理方法。
【図１】
【図２】
【図３】
【図４】
【図５】
【図６】
【図７】
【図８】
【図９】
【図１０】
【図１１】
【図２】
【図３】
【図４】
【図５】
【図６】
【図７】
【図８】
【図９】
【図１０】
【図１１】
【公開番号】特開２００７−２８５７８１（Ｐ２００７−２８５７８１Ａ）
【公開日】平成１９年１１月１日（２００７．１１．１）
【国際特許分類】
【出願番号】特願２００６−１１１４３６（Ｐ２００６−１１１４３６）
【出願日】平成１８年４月１４日（２００６．４．１４）
【新規性喪失の例外の表示】特許法第３０条第１項適用申請有り 平成１７年１０月１６日〜１８日 社団法人物理探査学会主催の「第１１３回（平成１７年度秋季）学術講演会」において文書をもって発表
【出願人】（５０４１５７０２４）国立大学法人東北大学 (2,297)
【Ｆターム（参考）】
【公開日】平成１９年１１月１日（２００７．１１．１）
【国際特許分類】
【出願日】平成１８年４月１４日（２００６．４．１４）
【新規性喪失の例外の表示】特許法第３０条第１項適用申請有り 平成１７年１０月１６日〜１８日 社団法人物理探査学会主催の「第１１３回（平成１７年度秋季）学術講演会」において文書をもって発表
【出願人】（５０４１５７０２４）国立大学法人東北大学 (2,297)
【Ｆターム（参考）】
[ Back to top ]