(57)【特許請求の範囲】
【請求項１】
発生した電磁波を地中に伝播させ、地中からの反射波をセンサにより検出し、検出結果を画像化し、画像化した前記反射波の状態に基づいて地中の物体を探知する方法において、
複数種類のセンサの反射波の画像化情報をそれぞれ２値化し、
２値化画像から、前記複数種類のセンサの検出特性に応じた検出物の特徴量をそれぞれ抽出し、
抽出したそれぞれの種類のセンサに対応する特徴量を、対応するそれぞれのニューラルネットワークの入力層に入力すると共に、抽出した複数種類のセンサに対応する特徴量を、対応する別のニューラルネットワークの入力層に入力して、地中の物体を識別し、
少なくともいずれか１つのニューラルネットワークにより対象とする物体であると識別された場合に、探知したものと判断することを特徴とする方法。
【請求項２】
前記複数種類のセンサの１つは、１つの地中レーダ又は複数のアレイ化された地中レーダを備え、前記特徴量は、該地中レーダによる検出物の２値化画像における走査方向サイズ、深度方向サイズ、断面積、断面円形度、及びレーダ後方断面積の少なくとも１つであることを特徴とする請求項１記載の方法。
【請求項３】
前記複数種類のセンサの１つは、１つの磁気センサ又は複数のアレイ化された磁気センサを備え、前記特徴量は、該磁気センサによる検出物の２値化画像における走査方向サイズ、センサアレイ方向サイズ、平面面積、及び平面円形度の少なくとも１つであることを特徴とする請求項１記載の方法。
【請求項４】
前記ニューラルネットワークは、バックプロパゲーション・ニューラルネットワークであることを特徴とする請求項１乃至３の何れかに記載の方法。
【請求項５】
前記地中の物体は、地中に埋設されている地雷であることを特徴とする請求項１乃至４の何れかに記載の方法。
【請求項６】
発生した電磁波を地中に伝播させ、地中からの反射波をセンサにより検出し、検出結果を画像化し、画像化した前記反射波の状態に基づいて地中の物体を探知する装置において、
複数種類のセンサの反射波の画像化情報をそれぞれ２値化する２値化手段と、
該２値化手段による２値化画像から、前記複数種類のセンサの検出特性に応じた検出物の特徴量をそれぞれ抽出する特徴量抽出手段と、
該特徴量抽出手段により抽出されたそれぞれの種類のセンサに対応する特徴量を、対応するそれぞれのニューラルネットワークの入力層に入力すると共に、抽出された複数種類のセンサに対応する特徴量を、対応する別のニューラルネットワークの入力層に入力して、地中の物体を識別する識別手段と、
を備え、
前記識別手段は、少なくともいずれか１つのニューラルネットワークにより対象とする物体であると識別した場合に、探知したものと判断することを特徴とする装置。
【請求項７】
前記複数種類のセンサの１つは、１つの地中レーダ又は複数のアレイ化された地中レーダを備え、前記特徴量は、該地中レーダによる検出物の２値化画像における走査方向サイズ、深度方向サイズ、断面積、断面円形度、及びレーダ後方断面積の少なくとも１つであることを特徴とする請求項６記載の装置。
【請求項８】
前記複数種類のセンサの１つは、１つの磁気センサ又は複数のアレイ化された磁気センサを備え、前記特徴量は、該磁気センサによる検出物の２値化画像における走査方向サイズ、センサアレイ方向サイズ、平面面積、及び平面円形度の少なくとも１つであることを特徴とする請求項６記載の装置。
【請求項９】
前記ニューラルネットワークは、バックプロパゲーション・ニューラルネットワークであることを特徴とする請求項６乃至８の何れかに記載の装置。
【請求項１０】
前記地中の物体は、地中に埋設されている地雷であることを特徴とする請求項６乃至９の何れかに記載の装置。
【請求項１１】
上記請求項６乃至１０の何れかに記載の装置を備えることを特徴とする移動体。
【発明の詳細な説明】
【０００１】
【発明の属する技術分野】
本願発明は、従来の方法では探知に熟練を要していた地中の物体であっても、軽易な訓練で効果的に探知することができ、探知処理を自動化することが可能な地中の物体を探知する方法及びその実施に使用する装置、並びに当該装置を搭載した移動体に関する。
【０００２】
【従来の技術】
対戦車地雷（ATM: Anti-tank Mines）及び対人地雷（APM: Anti-personnel Mines）に代表される地雷は、世界中の紛争地域で使用されている。現在、世界中で１億個を超える数の地雷が埋設されたままとなっており、紛争後におけるこれらの除去は困難を極めている。
【０００３】
紛争後の地雷除去、所謂、人道的地雷除去においては、地雷原となっている土地での生活を前提としていることから、その地雷原に埋設された地雷のすべてを除去しなければならないが、これらの地雷に付与された例えば「アンチ・ハンドリング・デバイス（Anti-handling Devices）」及び「アンチ・ディスターバンス・デバイス（Anti-disturbance Devices）」等の技術により、埋設された地雷の探知及び除去は困難となっているのが実情である。
【０００４】
これらの技術には、例えば、地雷探知器の磁場に反応して起爆するもの、数回の圧力付加によって起爆する種類のもの等があり、前者の種類では、作動中の地雷探知器を近付けただけで起爆する可能性がある。また、後者の種類では、地雷処理車が通過したときには起爆せず、これに続いて通過する地雷処理車以外の車両が通過したときに起爆する可能性がある。
【０００５】
このような実情の中、地雷を除去する方法としては以下の如きものが一般的である。例えば、地表からの圧力付加により起爆する種類の地雷に対しては、ホイールローラの如き態様の地雷処理車で地雷原を走行することにより爆破除去することができる。しかしながら、地雷原は平坦な土地ばかりではなく、地雷処理車が走行できない土地も多い。また、地雷処理車が走行できるような平坦な土地であっても、前述したような複数回の圧力付加によって起爆する種類の地雷にあっては、その信管の起爆までの圧力付加回数の設定が不明である以上、地雷処理車による地雷の完全な除去は望めない。
【０００６】
また、地雷原において外部から爆発を発生させ、そのときの爆風及び地面の振動等により、埋設された地雷の起爆を誘発させることも行なわれている。しかしながら、地雷の信管がその外部からの爆発により確実に作動するかどうかは不明であり、地雷の中には、非作動率が５％〜３０％と言われる自己不活性化機能付きのものもあり、やはり確実な除去は望めない。
【０００７】
そこで、従来から以上の如き除去方法に加え、残った地雷を探知して手作業にて掘り出すことによって除去することに頼らざるを得ないのである。
【０００８】
地雷の探知方法としては、金属探知器を利用する方法が主流であるが、対人地雷の多くは金属以外の材料を利用しているものが多いため、対人地雷を金属探知器のみで探知することは困難である。
【０００９】
別の探知方法としては、地中レーダ，磁気センサ，Ｘ線センサ，及び赤外線センサ等を利用し、発生した電磁波を地中に伝播させ、その反射波を検出して可聴音に変換して人間の聴覚をもって地雷を識別したり、反射波に基づいて地中の状態を画像表示して人間が視覚的に地雷を識別するものがある。この方法によれば、非金属，非磁性材料等からなる地雷であっても探知することができる（例えば、特許文献１乃至３参照）。
【００１０】
【特許文献１】
特許第3309242号明細書
【特許文献２】
特開平8-122279号明細書
【特許文献３】
特許第3015643号明細書
【００１１】
【発明が解決しようとする課題】
しかしながら、地中には、例えば金属製の異物が埋まっており、特に、紛争地帯であった土地には武器，弾薬，及び爆発物等の破片が多量に埋まっていることから、誤認識の割合が高くなるのが現状である。誤認識の割合を下げるために探知感度（反射波の信号強度の閾値）を下げることも考えられるが、そうすると小型の対人地雷を探知することができなくなるという新たな問題が生じる。
【００１２】
また、例えば、土中の水分が多い場合には、土中の水分が反射体となって誤認識を招いたり、土地によっては金属成分（多くは、鉄分）を多く含む土質である場合もあり、また、土質が不均一である場合もあって、このような土地における探知精度は更に低下する。
【００１３】
さらには、センサ類による聴音又は画像から地雷を識別するのは人間の判断に頼っており、その特定には多くの熟練を要することから、地雷探知にかかる需要に十分応えることができない。
【００１４】
しかし、現状においては上記の如き探知方法に頼らざるを得ないことから、たとえ誤認識であったとしても、その確証がない以上その場所をすべて手作業で掘り起こす必要があり、その作業量は甚大である。
【００１５】
本願発明は、上記状況に鑑みて行なわれたものであり、センサによる反射波の画像化情報を２値化し、２値化画像をニューラルネットワークを利用して識別することにより、探知処理を自動化することができ、また、オペレータの熟練を要しない高い探知精度を提供することができる地中の物体を探知する方法及びその実施に使用する装置、並びに当該装置を搭載した移動体を提供することを目的とする。
【００１６】
【課題を解決するための手段】
本願発明は、上記課題を以下のような構成からなる地中の物体を探知する方法及び装置、並びに当該装置を搭載した移動体によって解決することができる。
【００１７】
本願発明に係る地中の物体を探知する方法は、発生した電磁波を地中に伝播させ、地中からの反射波をセンサにより検出し、検出結果を画像化し、画像化した前記反射波の状態に基づいて地中の物体を探知する方法において、前記反射波の画像化情報を２値化し、２値化画像に基づいて、地中の物体をニューラルネットワークを利用して識別することを特徴とする。
【００１８】
また、本願発明に係る地中の物体を探知する装置は、発生した電磁波を地中に伝播させ、地中からの反射波をセンサにより検出し、検出結果を画像化し、画像化した前記反射波の状態に基づいて地中の物体を探知する装置において、前記反射波の画像化情報を２値化する２値化手段と、該２値化手段による２値化画像に基づいて、地中の物体をニューラルネットワークを利用して識別する識別手段とを備えることを特徴とする。
【００１９】
上記発明によれば、検出後に画像化（可視化）可能な電磁波を地中に伝播させ、センサにより、地中からの反射波に基づいて地中の状態を画像化し、まず、前記センサによる画像化情報（通常は、階調画像である）を２値化することによって、コンピュータ等の演算装置による自動処理を可能とする。そして、前記演算装置により、２値化画像をニューラルネットワークを利用して処理し、探知対象となる地中の物体（例えば、地中に埋設された地雷）を識別することにより、単純な閾値による従来の識別方法と比べて、高度な訓練なしに非常に高精度な探知を行なうことが可能となる。
【００２０】
なお、地中に伝播させる電磁波の発生機能は、本願発明に係る装置に備えさせてもよく、また、当該装置とは別体に設ける構成とすることも可能である。
【００２１】
上記センサとしては、地中レーダ、能動型及び受動型の磁気センサ，Ｘ線センサ，及び赤外線センサ等のその検出信号を画像化可能な種々のセンサを採用することが可能である。また、物体の識別精度を向上させるために、これら複数種類のセンサを組み合わせて利用することも可能である。好ましくは、例えば、インパルスセンサを備えた地中レーダ、及び能動型の磁気センサ（所謂、アクティブ磁気センサ）の２つを採用する。
【００２２】
上記ニューラルネットワークとしては、階層型であり、オンライン又はオフラインでの教育が可能なバックプロパーゲーション・ニューラルネットワーク（BPNN：Back-propagation Neural Network）を採用することが望ましく、このようなバックプロパーゲーション・ニューラルネットワークによれば、学習によって、より高精度な識別を行なうことができる。
【００２３】
このニューラルネットワークの入力層には、上記センサの検出結果（出力）の２値化画像から抽出した検出物の一又は複数の任意の特徴量を入力し、出力層としては、例えば、「地雷」，「非地雷」，「金属対人地雷」，「非金属対人地雷」，「金属対戦車地雷」，及び「非金属対戦車地雷」等の探知対象となる物体の種類に応じたカテゴリを設けることができる。
【００２４】
上記入力層に入力する検出物の特徴量としては、例えば、上記センサが地中レーダである場合には、地中レーダの「検出物の走査方向サイズ」，「検出物の深度方向サイズ」，「検出物の断面積」，「検出物の断面円形度」，及び「検出物のレーダ後方断面積（RCS：Radar Cross Section）」等とすることができる。なお、地中レーダとして、複数の地中レーダを採用することも可能であり、これら複数の地中レーダを走査方向に対して略直交するようにアレイ状に配置することも可能である。
【００２５】
また、上記センサが磁気センサである場合には、磁気センサの「検出物の走査方向サイズ」，「検出物の平面面積」，及び「検出物の平面円形度」等とすることができる。また、複数の磁気センサを採用することも可能であり、これら複数の磁気センサを走査方向に対して直交するようにアレイ状に配置することも可能である。このような場合には、上記特徴量として「検出物のセンサアレイ方向サイズ」を含めることが可能である。
【００２６】
また、複数種類のセンサを採用する場合には、例えば、複数種類のセンサによる検出物の特徴量をそれぞれ入力することができる複数のニューラルネットワークを採用し、これら複数のニューラルネットワークにそれぞれ単独処理又は並列処理するように構成することが可能である。このように、複数のニューラルネットワークを採用することにより、より正確な識別を行なうことができる。
【００２７】
また、例えば、複数種類のセンサによる検出物の特徴量を組み合わせて入力することができる複合的なニューラルネットワークを採用することも可能であり、さらにこれに加えて、このような複合的なニューラルネットワークに対して並列処理することが可能な、一又は複数のセンサによる検出物の特徴量を入力することができる別のニューラルネットワークを採用することもできる。
【００２８】
この具体的な例としては、第１、第２、及び第３の３つのニューラルネットワークを採用し、地中レーダによる検出物の特徴量を第１のニューラルネットワークの入力層に入力し、磁気センサによる検出物の特徴量を第２のニューラルネットワークの入力層に入力し、地中レーダ及び磁気センサの両方による検出物の特徴量を第３のニューラルネットワークの入力層に入力する。
【００２９】
このように、複数のニューラルネットワークで並列に処理した結果を利用することにより、何れかのニューラルネットワークで識別できないような検出物に対して、更に正確な識別を行なうことが可能である。
【００３０】
以上の如き本願発明に係る方法を地雷探知に適用することにより、従来から探知することが困難であった対人地雷でさえも、高度な訓練なしに非常に高精度な探知を行なうことが可能となる。
【００３１】
さらに、上記の如き本願発明に係る方法を実施することができる装置を地雷処理車の如き移動体に搭載して利用することも可能であり、また、オペレータが運搬することが可能な小型のものとすることも可能である。
【００３２】
【発明の実施の形態】
以下、本願発明に係る地中の物体を探知する方法を実施するための装置について、地雷探知の用に供される当該装置を搭載した移動体を例として添付の図面を参照しながら具体的に説明する。
【００３３】
図１は、本願発明の実施の形態に係る方法を実施するための装置（例えば、地雷探知装置）の機能ブロック図である。なお、図１に示した本実施の形態に係る装置又は地雷探知装置は単なる例示であり、本願発明においては、このような構成に限定されるものではない。
【００３４】
図１において、本実施の形態に係る地雷探知装置は、インパルスセンサを有する地中レーダ1と、アクティブ磁気センサ2との２種類のセンサを備えている。これら２種類のセンサは、能動型のセンサであり、地中レーダ1は、地中にパルス状の電磁波を放射し、その反射波を対応するセンサ部分で検出するようになっており、また、アクティブ磁気センサ2は、一定の周波数の送信磁界を与え、それにより得られる誘導磁界を対応するセンサ部分で検出するようになっている。また、上記両センサは、それぞれ複数の素センサを有し、該複数の素センサは、走査方向に略直交する方向にアレイ状に配置されている。
【００３５】
地中レーダ1は、地中の状態を深度方向に検出することができる一方、磁気センサ2は、地中の状態をアレイ方向（即ち、センサチャンネル方向）に検出することができ、地中レーダ1及び磁気センサ2による検出信号は、それぞれの信号処理部11，21に与えられる。
【００３６】
図２（ａ）〜（ｃ）は、本実施の形態に係る地雷探知装置の地中レーダによる後述の各処理における画像データの形態を示す模式図であり、これらの図においては、縦軸に深度方向を、横軸に走査方向をそれぞれとってある。また、図３（ａ）〜（ｃ）は、本実施の形態に係る地雷探知装置の磁気センサによる後述の各処理における画像データの形態を示す模式図であり、これらの図においては、縦軸にチャンネル方向を、横軸に走査方向をそれぞれとってある。
【００３７】
地中レーダ1及び磁気センサ2のそれぞれの検出時点での画像データは、図２（ａ）及び図３（ａ）の如くなっており、各図において、各センサによる検出物は、大きい楕円（破線で示す）で示す範囲内にある。なお、図３（ａ）においては、白色の楕円領域（磁気センサ2による検出信号が「負」の部分）と灰色の楕円領域（磁気センサ2による検出信号が「正」の部分）との間の中間位置が磁気センサ2による検出物である。
【００３８】
信号処理部11，21は、一般的なフィルタリング処理機能を有し、各検出信号からノイズ成分を除去する。信号処理部11によりノイズ除去された地中レーダ1の検出信号は、分解能向上処理部12に与えられ、一方、信号処理部21によりノイズ除去された磁気センサ2の検出信号は、目標領域抽出処理部31へ与えられる。
【００３９】
分解能向上処理部12は、合成開口処理等による高分解能化・画像化機能を有し、信号処理部11から与えられたノイズ除去後の検出信号の空間的分布及び時間的な応答を強調して画像化すると共に、探知目標となる地雷の外形形状（一般的には短寸円柱形状）に対応した特有のパターンを生成する。
【００４０】
本実施の形態においては、合成開口処理が実施され、該合成開口処理では、具体的には、走査方向（所謂、アジマス方向）毎に取得した時間軸上の検出信号群（以後、「検出データ」という）を、或る点反射体と地中レーダ1との間の距離（即ち、伝播距離）の変化量に基づいた参照関数との相関を用いてアジマス方向に圧縮（所謂、アジマス圧縮）することにより、空間分解能を向上させる。この合成開口処理後の画像データの例は、図２（ｂ）に示してある。なお、地中レーダ1の検出信号に対して合成開口処理を実施する場合には、上記アジマス圧縮と併せて、レーダ波の帯域幅に応じたレンジ方向での検出データの圧縮（レンジ圧縮又はパルス圧縮）によるレンジ方向分解能向上がなされる。
【００４１】
分解能向上処理部12による処理結果は、上記目標領域抽出処理部31に与えられ、該目標領域抽出処理部31は、信号処理部21から与えられたノイズ除去後の検出信号に基づいて走査方向毎に取得した時間軸上の検出信号群（以後、「検出データ」という）と、分解能向上処理部12から与えられたノイズ除去及び強調処理後の検出データとからそれぞれ生成される画像データに対して、２値化処理後にラベリング処理を施すことによって、検出した地中の物体（以後、「目標物」という）にラベル番号を付与する。
【００４２】
具体的には、目標領域抽出処理部31は、一般的な画像処理機能を有し、地中レーダ1の画像データに対しては、図４のフローチャートに示す如く、上記ラベリング処理（ステップS14）に先立って、メディアンフィルタによる平滑化処理（ステップS11）、２値化処理（ステップS12）、及び収縮／膨張処理（ステップS13）を施し、目標物画像を背景画像から分離・抽出する。
【００４３】
また、目標領域抽出処理部31は、磁気センサ2の画像データに対しては、図５のフローチャートに示す如く、上記ラベリング処理（ステップS25）に先立って、微分処理（ステップS21）、メディアンフィルタによる平滑化処理（ステップS22）、２値化処理（ステップS23）、及び収縮／膨張処理（ステップS24）を施し、目標物画像を背景画像から分離・抽出する。
【００４４】
なお、磁気センサ2の画像データに対しては、ステップS23の２値化処理に先立って、微分処理を施している。これは、図３（ａ）に示した如き検出時点での画像データを２値化した場合、前述した「正」の部分のみを抽出する可能性があり、また、微分処理を施すことによって、目標物の中心部分を抽出することができるからである。この微分処理後の画像データの例は、図３（ｂ）に示してある。
【００４５】
上記ステップS11及びS22における平滑化処理では、より詳しくは、メディアンフィルタを利用し、画像データから、信号強度が大きいスパイク状のノイズを除去するために、次の（１）式で示すように、注目画素の近傍領域内の中間値（median）を出力する。
【００４６】
【数１】
0003642775.tif 000002
【００４７】
また、ステップS12及びS23における２値化処理では、統計的な振幅変動の分散等を用いたデータ選別を行なうため、p-tile法を利用する。p-tile法において、背景画像よりも明るい目標物画像を切り出す場合の２値化閾値ftは、次の（２）式で与えられる。なお、本願発明においては、２値化の方法をp-tile法に限定するものではないことは言うまでもない。
【００４８】
【数２】
0003642775.tif 000003
【００４９】
但し、上記２値化閾値ftのみでは、目標物が含まれていない背景画像までも切り出す虞があるため、画像データの分散値と、予め求めておいた背景画像の分散値の上限とを比較し、分散値が上限よりも大きい場合にだけ、画像データを２値化処理する。なお、p-tile法に用いる２値化閾値ftは、予め実験等により求めておく。また、平均値m及び標準偏差σを求める領域は、目標物の大きさを基準とし、予め実験等により求めておく。
【００５０】
ステップS13及びS24における収縮／膨張処理では、後に続くラベリング処理での負荷を軽減すべく、次の（３）式及び（４）式を用いて収縮／膨張処理を行ない、孤立ノイズを除去する。
【００５１】
【数３】
0003642775.tif 000004
【００５２】
この収縮／膨張処理は、収縮処理を数回実施し、その後、膨張処理を同じ回数実施することが望ましく、その回数は、実験結果等に応じて調整可能なパラメータとしておく。
【００５３】
ステップS14及びS25におけるラベリング処理では、抽出された領域（即ち、値１の画素群）を個々の目標候補に分離する。具体的には、画像データ全体を走査し、図６に示す如く、「値１の注目画素」とその近傍領域内の「値１の画素」との連結性を調べ、連結した「値１の画素」には同一ラベルを付与する一方、連結していない「値１の画素」には新規のラベルを付与する。次に、湾曲パターンでは、同一目標物に対して複数のラベルが付与されるため、ラベルの統合処理を実施する。
【００５４】
より詳しく説明すれば、図６において画素を矩形領域で示し、注目画素が黒色領域で示される画素であるものとする。この場合、当該注目画素に「aij」のラベルを付与し、４連結型のラベリングを採用した場合には、この注目画素「aij」の四方４つの隣接画素を調べる。図６では注目画素「aij」の上側及び左側の画素がそれぞれ当該注目画素に「aij」に連結しているものとすれば、上側の画素に「ai-1 j」、左側の画素に「ai j-1」という「a」のラベルを付与する（図６においてはそれぞれの画素をクロスハッチング領域で示してある）。また、８連結型のラベリングを採用した場合には、上記注目画素「aij」の周囲８つの隣接画素を調べ、上述と同様の方法にて同一のラベルを付与する（図６においてはそれぞれの画素を白色領域で示してある）。
【００５５】
以上の如く、地中レーダ1及び磁気センサ2のそれぞれの画像データについてラベリング処理が完了すると、目標領域抽出処理部31は、ラベリング後の各画像データを目標位置判定処理部32に与える。
【００５６】
目標位置判定処理部32は、地中レーダ1の画像データと磁気センサ2の画像データとの抽出された目標物の位置の整合と、同一目標物であるか否かの判定とを行なう。
【００５７】
地中レーダ1及び磁気センサ2を併用して使用する場合には、これらを水平面上の同一位置に配置することができない場合、また、両者の計測周期が一致していない場合には、地中レーダ1及び磁気センサ2のそれぞれの画像データを位置的に整合させる必要がある。
【００５８】
まず、図７に示す如く、本実施の形態においては、走査方向に向かって、地中レーダ1，磁気センサ2の順で配置してあり、地中レーダ1及び磁気センサ2の間の距離がdであるものとする。
【００５９】
この場合、図８（ａ）に示す如く、地中レーダ1の計測周期（即ち、計測間隔）がΔSrであるとすれば、Kr回目の計測位置は、KrΔSrとなる。一方、図８（ｂ）に示す如く、磁気センサ2の計測周期（即ち、計測間隔）がΔSmであるとすれば、Km回目の計測位置は、KmΔSmとなる。これら各計測位置の整合は次の（５）式を用いて行なう。
【００６０】
【数４】
0003642775.tif 000005
【００６１】
次いで、地中レーダ1では、１つの目標物を複数のチャンネルで検出する可能性があるため、各チャンネルで検出した目標物が同一の目標物であるか否かを、目標物の位置を推定することによって行なう。
【００６２】
例えば、或るチャンネルn1で検出された目標物の位置（ラベリングされた目標物の最も信号強度が高い位置）の信号強度をTn1（n1，yn1，zn1）とし、これに隣接するチャンネルn2で検出された目標物の位置の信号強度をTn2（n2，yn2，zn2）とし、n1とn2との間の距離をdcとする。また、ここでは、yは走査方向の位置又は距離を表わし、zは深度方向の位置又は距離を表わす。
【００６３】
このとき、図９に示す如く、それぞれのチャンネルで検出された目標物の走査方向の距離が判定値ydよりも近いこと（（６）式参照）、それぞれのチャンネルで検出された目標物の深度方向の距離が判定値zdよりも近いこと（（７）式参照）の２つの条件の双方に合致しているものを１つの目標物としてラベリングし直す。なお、判定距離yd，zdは、予め実験等により求めておく。
【００６４】
【数５】
0003642775.tif 000006
【００６５】
上記のような、検出された目標物の同一性の判断においては、考慮するチャンネルの相互距離が大きいと精度が低下する可能性があるため、可及的に隣接するチャンネル同士で考慮することが望ましい。
【００６６】
次に、図１０（ａ）及び（ｂ）に示す如く、同一目標物としてラベリングされた領域について、隣接チャンネル間でチャンネルアレイ方向に直線近似により内挿補間して見掛け上の分解能を向上させた平面画像データを得る。その後、再びp-tile法により目標物画像を抽出し、その面積重心位置xを採用する。また、１つのチャンネルでしか検出しなかった目標物に対しても、その隣接チャンネルのデータと内挿補間し、その後、p-tile法により目標物画像を抽出する。
【００６７】
一方、磁気センサ2では、目標物の位置を、前述の目標領域抽出処理でラベリングされた目標物の領域の面積重心位置とし、その位置をTm（xm，ym）とする。なお、磁気センサ2では、深度方向の情報は得られず、目標物の水平方向位置のみ特定することが可能である。
【００６８】
次いで、地中レーダ1により検出された目標物の位置Tr（xr，yr，zr）と、磁気センサ2により検出された目標物の位置Tm（xm，ym）とが、同一の目標物のものであるか否かを判定する。即ち、両者の位置が平面上で所定範囲内にある場合（次の（８）式参照）に、双方により検出された目標物に対して、同一目標物としてラベリングをし直す。このとき、磁気センサ2により検出された１つの目標物に対して、上記所定範囲内に、地中レーダ1により検出された複数の目標物がある場合には、その信号強度が強い方を採用する。なお、上記所定範囲は、種々の実験により決定する。
【００６９】
【数６】
0003642775.tif 000007
【００７０】
具体的には、図１１（ａ）に示す如く、例えば、地中レーダ1により検出された目標物▲１▼の領域（図中ハッチングで示す）の面積重心位置を中心とした距離Dを半径とする円形範囲（図中破線で示す）を上記所定範囲とし、この所定範囲内に、磁気センサ2により検出された目標物aの領域の面積重心位置があれば、両者の目標物は同一物であると判定する。
【００７１】
一方、図１１（ｂ）に示す如く、例えば、地中レーダ1により検出された目標物▲２▼の領域（図中ハッチングで示す）の面積重心位置を中心とした距離Dを半径とする円形領域（図中破線で示す）を上記所定範囲とし、この所定範囲内に、磁気センサ2により検出された目標物bの領域の面積重心位置がなければ、両者の目標物は同一物でないと判定する。図１１（ｃ）の例では、目標物c及び目標物dは、上記所定範囲内にあるが、例えば、目標物cの信号強度が目標物dの信号強度よりも強い場合には、この信号強度が強い目標物cが地中レーダ1の目標物▲２▼と同一物であると判定し、目標物dは同一物でないと判定する。
【００７２】
また、地中レーダ1及び磁気センサ2の両方で検出された目標に対しては、走査方向の位置及び深さは、地中レーダ1の情報を、アレイ方向の位置は、磁気センサ2の情報をそれぞれ採用する。
【００７３】
なお、目標位置判定処理後の地中レーダ1及び磁気センサ2の２値化画像は、それぞれ図２（ｃ）及び図３（ｃ）に示す如きものである。
【００７４】
以上のような、地中レーダ1により検出された目標物に関する情報、磁気センサ2により検出された目標物に関する情報、並びに地中レーダ1及び磁気センサ2の双方で検出された目標物に関する情報は、目標位置判定処理部32から特徴量抽出処理部33にそれぞれ与えられ、特徴量抽出処理部33は、ラベリングされた各目標物から、その目標物らしさ（つまり、本願発明においては、地雷らしい形状をなしているということ）を表わす特徴量を、２値化画像データから抽出する。本実施の形態においては、地中レーダ1及び磁気センサ2に関する特徴量として、図１２に示す如きものを採用しており、対応するパラメータを割り当ててある。なお、本願発明においては、これらの特徴量に限定されるものではなく、センサの種類，及び後で詳述するニューラルネットワークの構成等に応じて適宜に決定すればよい（図１４参照）。
【００７５】
２値化画像データから上記特徴量を抽出する方法を次に説明する。
【００７６】
まず、地中レーダ1に関する特徴量にあっては、図１３（ａ）に示す如く演算する。
（１）目標物の走査方向サイズ（即ち、「R size x」）は、当該目標物の走査方向（図においてはx方向）の画素数（図においては8画素）である。
（２）目標物の深度方向サイズ（即ち、「R size y」）は、当該目標物の深度方向（図においてはy方向）の画素数を（図においては5画素）である。
（３）目標物の断面積（即ち、「R area」）は、当該目標物の総画素数（図においては28画素）である。
（４）目標物の断面円形度（即ち、「R roundness」）は、当該目標物の周囲長の画素数（図においては26画素）を「R perimeter」としたときに、（R roundness）＝ 4π×（R area）／（R perimeter）2 で演算される値（図においては、（R roundness）＝ 4π×（28）／（26）2 ＝ 0.5205である）である。
（５）目標物のレーダ後方断面積（低域）（即ち、「Rcs Low」）は、当該目標物の低域RCSであり、次の（９）式で演算される。
（６）目標物のレーダ後方断面積（高域）（即ち、「Rcs High」）は、当該目標物の高域RCSであり、次の（１０）式で演算される。
【００７７】
【数７】
0003642775.tif 000008
【００７８】
一方、磁気センサ2に関する特徴量にあっては、図１３（ｂ）に示す如く演算する。
（１）目標物の走査方向サイズ（即ち、「M size x」）は、当該目標物の走査方向（図においてはx方向）の画素数（図においては8画素）である。
（２）目標物のチャンネル方向サイズ（即ち、「M size z」）は、当該目標物のチャンネル方向（図においてはz方向）の画素数を（図においては5画素）である。
（３）目標物の平面面積（即ち、「M area」）は、当該目標物の総画素数（図においては28画素）である。
（４）目標物の平面円形度（即ち、「M roundness」）は、当該目標物の周囲長の画素数（図においては26画素）を「M perimeter」としたときに、（M roundness）＝ 4π×（M area）／（M perimeter）2 で演算される値（図においては、（M roundness）＝ 4π×（28）／（26）2 ＝ 0.5205である）である。
【００７９】
特徴量抽出処理部33は、以上の如くにそれぞれのセンサの２値化画像データから抽出した特徴量、即ちパラメータの値を識別処理部4に与える。
【００８０】
識別処理部4は、学習機能付きの所謂「階層型バックプロパゲーションニューラルネットワーク（BPNN）」機能を備えた演算装置であり、図１４に示す如く、特徴量抽出処理部33から与えられた各目標物に対応する特徴量をそのニューラルネットワークの入力層に入力し、予め設定してある出力カテゴリへの該当度合を出力するように構成されている。
【００８１】
より詳しくは、図１５に示す如く、識別処理部4は、地中レーダ／磁気センサ同時検出用バックプロパゲーションニューラルネットワーク41，地中レーダ検出用バックプロパゲーションニューラルネットワーク42，及び磁気センサ検出用バックプロパゲーションニューラルネットワーク43の３つのニューラルネットワーク機能を備えており、地中レーダ1及び磁気センサ2の両方で同時に検出された目標物に関する特徴量は、地中レーダ／磁気センサ同時検出用バックプロパゲーションニューラルネットワーク41の入力層へ、地中レーダ1のみで単独に検出された目標物に関する特徴量は、地中レーダ検出用バックプロパゲーションニューラルネットワーク42の入力層へ、そして、磁気センサ2のみで単独に検出された目標物に関する特徴量は、磁気センサ検出用バックプロパゲーションニューラルネットワーク43の入力層へそれぞれ入力されるようになっている。
【００８２】
図１５からも判る通り、図１４に示したニューラルネットワークの構成は、図１５に示した地中レーダ／磁気センサ同時検出用バックプロパゲーションニューラルネットワーク41のものであり、他のニューラルネットワークについては、パラメータを各センサに対応するもののみとする等して、全体としては同様の構成とすることができるため、詳細な説明は省略する。
【００８３】
識別処理部4が備える各ニューラルネットワークは、予めオフラインで学習させておく。具体的には、まず、既知の地雷サンプルを、土質及び埋設深度等の埋設条件を変えながら、本願発明に係る地雷探知装置によりこれらの地雷サンプルの探知トライアルを繰り返す。そして、各トライアルにおいて当該地雷探知装置が検出した目標物の中から、実際に地雷を埋設した位置と略一致する目標物を選択し、この目標物に、埋設した地雷の出力カテゴリを割り当て、ニューラルネットワークの誤差逆伝播学習機能を利用し、結合荷重の修正を行う。なお、本実施の形態においては、上記出力カテゴリを、図１４に示す如く「地雷」，「非地雷」，「金属対人地雷」，「非金属対人地雷」，「金属対戦車地雷」，及び「非金属対戦車地雷」等としてあるが、これらに限定するものではない。なお、地雷の埋設位置以外の位置には、すべて「非地雷」を割り当てるようにする。
【００８４】
ニューラルネットワークの学習アルゴリズムは、図１６のフローチャートに示す一般的な手順である。具体的には、識別処理部4は、まず、記憶している結合荷重の初期化をし（ステップS101）、識別処理部4に接続された入力部5から入力される上記埋設条件、及び該当する出力カテゴリ（所謂、教師信号）等を含む学習パターンをセットする（ステップS102）。
【００８５】
次に、識別処理部4は、探知トライアルを開始し、中間層ユニットの出力計算（ステップS103）、出力層ユニットの出力計算（ステップS104）を順次行なう。
【００８６】
続いて、識別処理部4は、入力された出力カテゴリに対する、出力層ユニットの誤差計算（ステップS105）、中間層ユニットの誤差演算（ステップS106）を順次行ない、演算誤差に基づいて、出力層ユニットの出力が、上記入力された出力カテゴリに近づくように、中間層及び出力層間の結合荷重の修正（ステップS107）、入力層及び中間層間の結合荷重の修正（ステップS108）を順次行なう。
【００８７】
そして、識別処理部4は、探知トライアルの条件を変更すべく、セットした学習パターンを更新し（ステップS109）、学習繰り返し回数が、学習パターンのバリエーションに応じた所定の学習繰り返し回数以内か否かを確認する（ステップS110）。識別処理部4は、学習繰り返し回数以内である場合（ステップS110で“YES”）には、更新した別の学習パターンで学習すべくステップS103に戻り、一方、学習繰り返し回数以内でない場合（ステップS110で“NO”）には、全ての学習パターンのトライアルを終了したものと判断して学習を終了する。
【００８８】
以上の如き学習は、識別処理部4が備える地中レーダ／磁気センサ同時検出用バックプロパゲーションニューラルネットワーク41，地中レーダ検出用バックプロパゲーションニューラルネットワーク42，及び磁気センサ検出用バックプロパゲーションニューラルネットワーク43の３つのニューラルネットワークについてそれぞれ実施する。なお、本実施の形態においてはオフライン学習を想定して説明したが、オンラインで学習を実施することも可能である。
【００８９】
このような識別処理部4による識別結果は、図１７（ａ）及び（ｂ）並びに図１８に示すような形態で、該識別処理部4に接続された表示部6に表示させる。
【００９０】
図１７（ａ）及び（ｂ）においては、地中レーダ1のみにより検出された目標物と、磁気センサ2のみにより検出された目標物と、地中レーダ1及び磁気センサ2の両方で検出された目標物とを視覚的に区別して表示するとともに、各目標物の３次元位置が視覚的に判り易いように、地中の状態を水平断面視及び垂直断面視で表現してある。
【００９１】
より詳しくは、図１７（ａ）においては、横軸に走査方向、縦軸にチャンネル方向の水平距離をそれぞれとった矩形領域内に、検出された目標物が抽出後の形状でそれぞれ表示されている。一方、図１７（ｂ）においては、横軸に走査方向、縦軸に深度方向をそれぞれとった矩形領域内に、検出された目標物が抽出後の形状でそれぞれ表示されている。
【００９２】
また、図１７（ａ）及び図１７（ｂ）に示す如く表示された各目標物に対しては、上記ニューラルネットワークを利用して演算されたそれぞれの出力カテゴリに対する該当度合（識別率）を数値で表示してもよく、また、図１８に示す如く棒グラフの形態で表示してもよい。
【００９３】
図１８の例において、「地雷」及び「非地雷」の２つの出力カテゴリでは、地雷の確率の方が非常に大きく、地雷であることが判る。次いで、地雷であればどのような種類のものであるかは、残りの「金属対人地雷」，「非金属対人地雷」，「金属対戦車地雷」，及び「非金属対戦車地雷」の４つの出力カテゴリを見ることにより判る。つまり、この例では、金属対人地雷の確率が他に比べて大きく、金属対人地雷であると判断できる。
【００９４】
本実施の形態においては、図１７（ａ）及び（ｂ）並びに図１８に示すような画面表示をオペレータが見て、目標物の最終的な判定を行うように説明したが、このような判定処理をコンピュータに行なわせる構成とすることも可能である。
【００９５】
以上に説明した本実施の形態に係る地雷探知装置は、図１９に示す如く、適宜の移動体100に搭載することも可能である。移動体100には、地雷探知装置の全体部分を搭載する構成としてもよいし、また、少なくともそのアクティブセンサ部分（地中レーダ1及び磁気センサ2等）のみを搭載し、別の位置に配置された地雷探知装置の残りの部分と例えば無線通信によりデータを授受をなし得るように構成してもよい。
【００９６】
図１９の例では、移動体100は、ホイールローダの如き態様をなしており、その前方部分に昇降自在に動作可能な一対のアーム101を備え、該アームの先端部には、地中レーダ1及び磁気センサ2を収納した箱形の走査部102が水平に支持されている。該走査部102は、地表から所定距離離隔してアーム101に支持されており、その下側面から地雷探知のための電磁波を地中に放射し、地中からの反射波を受信することが可能なように、走査部102の下側の部位に地中レーダ1及び磁気センサ2が配置されている。
【００９７】
また、他の実施の形態に係る地雷探知装置は、図２０に示す如く、可搬型の小型のものとすることができる。この可搬型地雷探知装置200は、全体として棒状をなし、その上側となる一方の端部にオペレータが把持する把持部203，204，205等を備えている。また、他方の端部には、円盤状の走査部202を備えている。該走査部202の内部には、地中レーダ1及び磁気センサ2を収納しており、これらのアクティブセンサ部分が、その下側面から地雷探知のための電磁波を地中に放射し、地中からの反射波を受信することが可能なように配置されている。これらのアクティブセンサ部分は、棒状の装置本体の上部に配置された箱形のコントローラ201に接続され、該コントローラ201は、入力部5及び表示部6等を含む、上記アクティブセンサ部分以外の地雷探知装置の部分を備えている。
【００９８】
【発明の効果】
本願発明に係る地中の物体を探知する方法及びその実施に使用する装置、並びに当該装置を搭載した移動体によれば、センサによる反射波の画像化情報を２値化し、２値化画像をニューラルネットワークを利用して処理することにより、探知処理を自動化することができ、また、オペレータの熟練を要しない高い探知精度を提供することができる。
【図面の簡単な説明】
【図１】  本願発明の実施の形態に係る地雷探知装置の構成を示す機能ブロック図である。
【図２】  本実施の形態に係る地雷探知装置の地中レーダにより検出された地中の状態を示す画像データの模式図であり、（ａ）は地中レーダによる検出時点での元画像、（ｂ）は合成開口処理後の画像、（ｃ）は目標判定処理後の画像である。
【図３】  本実施の形態に係る地雷探知装置の磁気センサにより検出された地中の状態を示す画像データの模式図であり、（ａ）は磁気センサによる検出時点での元画像、（ｂ）は微分処理後の画像、（ｃ）は目標判定処理後の画像である。
【図４】  図１に示した目標領域抽出処理部が地中レーダによる画像データを処理する手順を示すフローチャートである。
【図５】  図１に示した目標領域抽出処理部が磁気センサによる画像データを処理する手順を示すフローチャートである。
【図６】  図１に示した目標領域抽出処理部のラベリング処理を説明するための図である。
【図７】  図１に示した目標位置判定処理部が、検出された複数の目標物の同一性を判定する処理を説明するための図である。
【図８】  図１に示した目標位置判定処理部が、検出された複数の目標物の同一性を判定する処理を説明するための図である。
【図９】  図１に示した目標位置判定処理部が、検出された複数の目標物の同一性を判定する処理を説明するための図である。
【図１０】  図１に示した目標位置判定処理部が、地中レーダの見掛け上の分解能を向上させる内挿補間を説明するための目標物画像の模式図である。
【図１１】  図１に示した目標位置判定処理部が、地中レーダにより検出された目標物と、磁気センサにより検出された目標物との同一性を判定する処理を説明するための図である。
【図１２】  図１に示した特徴量抽出処理部が識別処理部のニューラルネットワークの入力層に与えるパラメータ、即ち地中レーダ及び磁気センサによる目標物の特徴量を示す図表である。
【図１３】  図１に示した特徴量抽出処理部が２値化画像から各特徴量を演算する方法を説明するための図である。
【図１４】  図１に示した識別処理部のニューラルネットワークの構成を示す模式図である。
【図１５】  図１に示した識別処理部のニューラルネットワークの構成を示すブロック図である。
【図１６】  図１に示した識別処理部のニューラルネットワークの学習アルゴリズムを示すフローチャートである。
【図１７】  図１に示した識別処理部が、地雷探知結果を表示部に表示させる形態を示す模式図であり、（ａ）は地中の状態の水平断面視による表示、（ｂ）は地中の状態の垂直断面視による表示をそれぞれ示している。
【図１８】  図１に示した識別処理部が、地雷探知結果を表示部に表示させる形態を示す模式図であり、各目標物の出力カテゴリに対する該当度合（識別率）を示している。
【図１９】  本願発明の実施の形態に係る地雷探知装置を搭載した移動体の実施の形態を示す正面図である。
【図２０】  本願発明の他の実施の形態に係る地雷探知装置を可搬型とした実施の形態を示す斜視図である。
【符号の説明】
1  地中レーダ
2  磁気センサ
4  識別処理部（バックプロパゲーション・ニューラルネットワーク）
5  入力部
6  表示部
31  目標領域抽出処理部
33  特徴量抽出処理部
100  移動体
200  可搬型地雷探知装置
【図１】
0003642775.tif 000009
【図２】
0003642775.tif 000010
【図３】
0003642775.tif 000011
【図４】
0003642775.tif 000012
【図５】
0003642775.tif 000013
【図６】
0003642775.tif 000014
【図７】
0003642775.tif 000015
【図８】
0003642775.tif 000016
【図９】
0003642775.tif 000017
【図１０】
0003642775.tif 000018
【図１１】
0003642775.tif 000019
【図１２】
0003642775.tif 000020
【図１３】
0003642775.tif 000021
【図１４】
0003642775.tif 000022
【図１５】
0003642775.tif 000023
【図１６】
0003642775.tif 000024
【図１７】
0003642775.tif 000025
【図１８】
0003642775.tif 000026
【図１９】
0003642775.tif 000027
【図２０】
0003642775.tif 000028