(57)【要約】      （修正有）【課題】水量に関するレーダー画像内にある対象物を検知するデータ処理アルゴリズムを改善する自動目標識別システムの提供。【解決手段】水中で目標物を検知する方法は、対象領域の複数の画像を捉える工程と、画像からボクセルデータを抽出する工程と、ボクセルデータを処理して対象領域の対象項目を検知する工程とからなる。【選択図】図２
2008292449.tif 000002
【特許請求の範囲】
【請求項１】
光が透過する媒体内、又は該媒体に浮いている対象物を検知する方法であって、
対象領域の複数の画像を捉える工程と、
画像からボクセルデータを抽出する工程と、
ボクセルデータを処理して、対象領域の対象項目を検知し且つ分類する工程とを有する方法。
【請求項２】
更に、ボクセルデータの歪みを戻す工程を有する、請求項１に記載の方法。
【請求項３】
対象領域の複数の画像を捉える工程は、
送信されたビームをレーダー送信器から対象であるボクセルに伝搬する工程と、
反射されたビームを対象であるボクセルからレーダー送信器に伝搬する工程を有する、請求項１に記載の方法。
【請求項４】
更に、ボクセルデータを処理して、副産物が誘発されたシステムを除去する工程を有する、請求項１に記載の方法。
【請求項５】
画像からボクセルデータを抽出する工程は、連続した画像フレームの一群をスーパーフレームに平均化する工程を有する、請求項１に記載の方法。
【請求項６】
各スーパーフレームは、以前及びそれに続く後のスーパーフレームと、複数の画像フレームを共有する、請求項５に記載の方法。
【請求項７】
更に、光が透過する媒体の表面の歪みを戻す工程と、
標準的な列に基づいて適応衰滅モデルを用いて深さに依存した強度均等化を実行する工程を有する、請求項１に記載の方法。
【請求項８】
減衰係数ｋlidarとして知られる水の混濁度の測定は、標準的な列から評価される、請求項７に記載の方法。
【請求項９】
ボクセルデータを処理して、対象領域の対象項目を検知する工程は、
ボクセルデータの正規化された部分行列を、格納されたテンプレートと比較する工程を有する、請求項１に記載の方法。
【請求項１０】
ボクセルデータの正規化された部分行列を、格納されたテンプレートと比較する工程は、
ボクセルデータの各画素と、テンプレートに格納された対応する画素とポイントごとに掛け合わせて合計して、合計を生成する工程と、
合計を正規化する工程と、
合計をスレッショルドと比較する工程を有する、請求項９に記載の方法。
【請求項１１】
部分的な合計の行列は、ボクセルデータの正規化された部分行列と格納されたテンプレートとの比較演算を加速するのに用いられる、請求項１０に記載の方法。
【請求項１２】
スレッショルドは深さに依存する、請求項１０に記載の方法。
【請求項１３】
検査対象の対象項目を分類する工程は、外観、サイズ及びアスペクト比を用いて目標を識別する、請求項１に記載の方法。
【請求項１４】
検査対象の表面項目は、影一致探知が行われる、請求項１に記載の方法。
【請求項１５】
更に、検査対象の対象項目を集める工程を有する、請求項１に記載の方法。
【請求項１６】
対象の表面項目について、ボクセルデータは深さ次元に折り畳まれ、集める工程は方位角のみに行われる、請求項１５に記載の方法。
【請求項１７】
対象領域の複数画像を捉える複数のセンサと、
画像からボクセルデータを抽出し、該ボクセルデータを処理して対象領域の対象項目を検知するプロセッサを具えた装置。
【請求項１８】
プロセッサは、ボクセルデータの歪みを戻す、請求項１７に記載の装置。
【請求項１９】
プロセッサは、続いた画像フレームのグループをスーパーフレームに平均化する、請求項１７に記載の装置。
【請求項２０】
各スーパーフレームは、前の及び次のスーパーフレームと複数の画像フレームを共有する、請求項１９に記載の装置。
【請求項２１】
プロセッサは、対象項目を分類する、請求項１７に記載の装置。
【発明の詳細な説明】
【技術分野】
【０００１】
政府所有であることの陳述
本発明は、契約番号Ｎ６１３３１−００−Ｃ−００２２に基づいて成された。米国政府は契約の下、本発明に対する権利を有する。
【０００２】
本発明は、水中で対象物を検知し分類する方法及び装置に関する。
【背景技術】
【０００３】
発明の背景
水中の機雷(mine)を捜し出し、識別する種々のシステムが提案されてきた。機雷を検知する１つのシステムは、アメリカ合衆国海軍用に開発され構築されたノースロップグラマン空中レーザー機雷検知システム(ＡＬＭＤＳ)である。ＡＬＭＤＳシステムは、一連の画像の形式で情報を捉えるレーザー検知及び分類システム(ＬＡＤＡＲ)である。各画像は、検査対象の水分量を通して略垂直なスライス画像を表す。
【０００４】
ＡＬＭＤＳシステムは、パルス状のファンビーム送信器、４つのストリーク式受像管、及びデータ処理ユニットを含む。各受信器は、ストリーク式受像管画像化装置、増強装置、画像を捉えるカメラ、及び水量に機雷が含まれているかを判断するデータ処置機器から構成される。各受信器は、送信されたファンビームの約１／４から反射された光を記録する。全てのこのハードウエアは、一般にヘリコプターに取り付けられるポッドに含まれる。
【０００５】
本発明は、水量に関するレーダー画像内にある対象物を検知するデータ処理アルゴリズムを改善することを追求する。
【発明の開示】
【０００６】
発明の要約
本発明は、部分的に透明な材料、一般に海水内にて、ユーザが選択したサイズの対象物を検知する方法を提供する。方法は、検査対象の体積についての複数の画像を捉える工程と、画像からボクセルデータ(voxel data)を抽出する工程と、ボクセルデータを処理して水量内、又は水量の材料の表面に浮いている対象物を検知する工程を有する。水量内の材料が海水であるとき、これらの対象物は、「コンタクト」としてアメリカ海軍に知られている。本発明は更に、「コンタクト」を機雷のような対象物又は非機雷のような対象物の何れかとして分類する。
【発明を実施するための最良の形態】
【０００７】
発明の詳細な記載
図１に関して、図１は本発明の実施例に従って、構成され動作する装置(10)の概略図である。装置は、航空機(14)上に取り付けられたレーザーレーダーシステム(12)を含む。レーダーシステムのレーザーは、信号を空間的なファンビームの一時的なパルスの形式で送信し、信号は対象領域(16)に向けられる。特にこの場合、対象領域は、水域を含む。
【０００８】
図２は、本発明の方法の一実施例を実行するのに用いられる装置のブロック図である。システムは、レーザー(18)及び、航空機に取り付けられる１又は２以上のカメラ(20)の形式である単一又は複数台のセンサを含む。レーザーは、送信された光のパルスで対象領域を照射する。送信された信号は、対象領域内の対象物に反射されて、カメラによって受信される戻り信号を生成し、カメラはボクセルの形式で対象領域の画像を生成する。ボクセルデータは、ポッド装備信号処理デバイス(22)にて処理され(また、データプロセッサとも言われる)、対象領域の対象物を検知し区別する。反射データはカメラ(22)によって受信され、生データレコーダ(ＲＤＲ)によって記録される。カメラによって観察された反射に対応するボクセル値は、データプロセッサ内で予備処理され(26)、画像処理され(28)、検知され(30)及び区別される。検知された機雷のような対象物はデータベース(34)及び／又は他の出力デバイス(36)に送られる。本発明はデータベース内に含まれるアルゴリズムを用いて実行される。
【０００９】
ここで用いられる、「生データ」は電子カメラ(20)を出たデータである。生データは、対象領域の量に関する情報を表す。生データは以下のように収集される。ＡＬＭＤＳの例に於いて、ライダーはライダープラットホームの進行方向に直交する６０度のファンビームの形状でパルスを送信する。センサは「押し方式(pushbroom)」センサとして知られている。水量は、１回で１の「スライス」を画像化する。水量を表すデータは、順序付けられた(ordered)一組のスライスである。これらのスライスは、また「フレーム」と呼ばれる。小さなオーバーラップを除いて、各カメラは送信されたファンビームからの反射の特定の角度セクタから反射を記録する。一般に、ファンビームは約６４８フィートを横切り、各カメラはファンビームの約１５度、即ち約１６２フィートから生じる反射を記録する。データプロセッサでは、これらのフレームは、後でスーパーフレーム(図４)として知られる８つのグループに集約される。
【００１０】
ＡＬＭＤＳによって生成された生のデータは、円柱座標内に収集され、その中では列が円錐角に対応し、行が斜距離(slant range)に対応し(飛行時間に対応する)、フレームは進路沿い(along-track)の光景に対応し、即ち円柱座標の用語で、角度シータ、半径ｒ、軸Ｚに夫々対応する。
【００１１】
送信されたパルスは、送信信号と言われ、受信されたパルスは受信信号と言われる。１つの送信されたパルスは、多くの受信パルスに膨らみを与える、なぜなら水自体の量要素は各々が、送信されたパルスの小さいが探知可能な部分をレーダーに戻し反射し、これらは時間増分で解かれるからである。しばしば、水面からの明るいフラッシュが、水中の対象物からの反射と同様に、この「容積測定」のリターンに追加される。送信された信号は、光学系、空気、対象領域の水上の霞、水面、水、を含むチャンネルを通り、次に(もし、あれば)目標に反射する。
【００１２】
ここで用いられるように、「コンタクト」は水、海底、波の何れでもないものの検知を意味する。コンタクトは、例えば、魚、機雷、又は浮いている破片である。「分類」はコンタクトが機雷のような対象物(ＭＬＯ)又は非機雷のような対象物(ＮＭＬＯ)の何れかであることの表示(declaration)である。この用語は、本発明が機雷を見付けるとの当初の動機付け故に用いられる。しかし、カーネルの形を変えることにより(図６及び図７)、他の形を有する物体、例えば魚のような物体を見つけるのに、本発明が用いられる。「識別」とは機雷のタイプを判断することを意味し、種々の現存のアルゴリズムの変数を用いることにより、より詳細に形を決定することが可能である。
【００１３】
検知されたら、検知と地理位置情報(geolocation)のような更なる情報を有するボクセルの生データキューブ及び水の減衰定数は、データキューブとともに、データベース(34)に格納される。生データレコーダ(ＲＤＲ)(24)と言われる他のデータベースは、１又は２以上のカメラから流れる生画像を同時に記録することができる。
【００１４】
水を透過するレーザーレーダーには、係留している機雷に対してＭＬＯ又はＮＭＬＯであるとして、データキューブの識別を難しくしている２つの様相がある。１つ目の様相は、水(及び他の要素)により起こされる吸収と散乱であり、２つ目の態様は目標からの反射された光のレベルが低いことである。減衰が強いと、ぼやけと結びついて、機雷から受信された反射は、このアルゴリズムによって検知されても、画像データ内には肉眼によって識別されない。ぼやけ及び検知の統計特性が既知であれば、統計的評価により、更なる改善が可能である。
【００１５】
図３に示すように、一実施例において、本発明の方法は４つの工程を有する。前処理工程(40)、画像処理工程(422)、検知工程(44)、及び識別工程(46)である。前処理工程は、データ流れを減じ、受信器の副作用物を除去するのに用いられる。画像処理工程は、対象領域を決定し、検知器の性能を低下させるデータ内のバイアスを除去する。検知工程は、修正された尤度比表現を用いてコンタクトを検知する。識別工程は、外観、サイズのような特徴を演算して、該特徴は機雷のような対象物(ＭＬＯs)又は非機雷のような対象物(ＮＭＬＯs)を区別するのに用いられる。
【００１６】
ここでの自動目標認識(ＡＴＲ)は、警告／確認連続検知器と考えられる。検知器は、各塊量につて、即ちこれはコンタクト又は否かを決定する。検知後に、第２段階はコンタクトを機雷状の対象物(ＭＬＯ)又は非機雷状の対象物(ＮＭＬＯ)かを識別する。信号処理は空中レーザー機雷検知システム(ＡＬＭＤＳ)補助タンクにてデータ処理ユニット内のソフトウエアを用いて実行され得る。
【００１７】
各カメラからのデータは、一般に行列ビューに画素値の行列として表示され、データは直角座標システムに表されるとの誤った印象を与える。例えば、ＡＬＭＤＳが陸上のデータ収集に用いられるならば、データは円筒座標システムに収集される。これは対象領域が海面のみの場合も同じである。列は、円錐角に直に対応し、行は原点(センサ)からの半径距離に対応し、フレーム(進路沿いの動き)は、Ｚ軸(航空機の飛行方向)に沿った動きに対応する。これは、空気と水面の境界に於いてスネルの法則による(Snell)光線の屈折を考慮し、光量子の飛行時間が水中の光速度を考慮に入れて解釈されなければならない点を除き、水面下のデータについても当てはまる。
このように、角度／半径画素の行列は、しばしばカルテシアン配列の行列にて示されるが、実際には完全なカルテシアンではない。及びスネルの法則の屈折故に、正確にはまっすぐな円筒座標システムでもない。このように(カルテシアンフレームの行列として)データを示すことは、レーダーによって捉えられた画像をゆがめることである。特に、海洋の場合は、局部的には平均すると平坦であると誰もが知っているが、ＡＬＭＤＳデータでは、平均すると、曲がって表れる。センサが等間隔でなければ、表面のカーブは数学的にロール角(又は、列)について正弦関係となる。データを本来(native)の座標システムからカルテシアンシステムに変更することは、オルト調整(ortho-rectification)と呼ばれる。これは水面下のデータでは行われない、なぜなら対象物である機雷のサイズの歪みは小さく、オルト調整は、コンピュータ的に高価だからである。この種の処理が必要であれば、目標のテンプレートは、全てのデータを処理するのではなく、データの歪みと整合するように修正されるであろう。
【００１８】
前処理工程は、２つの機能を有する。(１)受信器の副作用物(artifact)を除去すること(２)データ流れを減じることである。例えば、カメラ画像の垂直列中に、データではなく、受信器の副作用物しかないときは、この列は左及び右のデータの平均値と置き換えられる。
レーダーデータは、８つの連続したフレームの各グループを平均化することによって、１つのスーパーフレーム(図４)へとグループ化(binned)される。このように、スーパーフレーム内の各画素は、「生」のフレームに同じ行列指数を有する８つの画素の平均である。各スーパーフレームは、画像処理の為に送られ、次に検知器に送られる２次元の対象物である。グループ化は２つのストライドで実行される。言い換えれば、各スーパーフレームは、その構成の中に、前のスーパーフレーム及び次のスーパーフレームにある６つのフレームを共有する。
【００１９】
ライダーシステムが飛行すると、フレームは、秒当たりの所定のフレーム数の速度で収集される。図４の例は、各スーパーフレームは、時間内に連続して収集された８つの生のフレームの平均であることを示している。次のスーパーフレームは、以前のスーパーフレームと６つのフレームを共有する。これは各画素の移動平均を構成するのと近似している。  数学的に、グループ化の効果は、非干渉性積分(non-coherent integration)の効果と近似しており、進路沿いの次元(along-track dimension)にてパルス整合したフィルタを用いた畳み込みに似ている。
【００２０】
ＡＬＭＤＳカメラフレームデータは、０〜５１１まで番号付けられた５１２の列に配列される。ＡＬＭＤＳレシーバーは四分円境界で副作用物を有し、前のフレーム５１１番目の列は、現在のフレームの構造の０番目の列に複写される。四分円境界は１２８、２５６及び３８４番目の列に位置する。四分円境界はデータがない列である。前処理しないと、それらは、検知器の実行に影響する場合がある。
四分円境界列(１２８、２５６、３８４)のデータは、左右の列のデータの平均値と置き替えられる。５１１列を０列にコピーすることによって、５１１列(右端の列)を０列(左端の列)に複製化することが行われる。
異なる副作用物又は異なるデータ速度を有する異なるＬＩＤＡＲが使用された場合、又はデータプロセッサに異なる計算容量(ＦＬＯＰＳで)があった場合、この手順への適切な変更が為されねばならない。
【００２１】
画像処理段階は、海面が水平に見えるという点でオルト調整に似ている工程を有する。
工程はフレームを「ゆがめない」として知られている。ゆがめないために、画像の各列に相当する円錐角を知る必要がある。
従って、円錐角は、画像内の各列が対応するセンサ内の視角と同様に、プラットフォーム(例えばヘリコプター)の地上に対する向き、及びプラットフォームに対するセンサの向きに依る。各列の大凡の円錐角を以下の公式から得てもよい。
αcone(ｊ)＝ｃｏｓ-1(ｃｏｓ(φroll(ｊ))ｃｏｓ(θplatpitch))         (１)
ここで、θplatpitchはプラットフォームのピッチであり、φrollは特定の列に対応したプラットフォームロールとロール角の合計である。
φroll(ｊ)＝φplatroll＋φcolroll(ｊ)                     (２)
【００２２】
センサ座標システムに於いて(航空機の座標システムに対するものとして)、特定の列に対応するロール角ｊは、以下によって与えられる。
φcolroll(ｊ)＝−０.５０７８９１＋ｒａｄｉａｎｓ／ｃｏｌ×
［(５１２−２０.４８)(ＮcamID−１)＋ｊ＋０.５］,                  (３)
ここで、ＮcamIDは図８に示す約束事に従ったカメラ識別(ＩＤ)である。
【００２３】
海面の行変位(天底(nadir)に於ける海面の行インデックス(row index)に対して)は、以下で与えられる。
【数１】
2008292449.tif 000003
ここで、ｈは高度である。
【００２４】
ソフトウエアは、一実施例に於いてカメラは２０.４８列の重なりがあり、及び視野の縁は１４.１、０.３、−０.３及び−１４.１度の回転角だけ離れていることを考慮に入れる。
【００２５】
Δｉ(ｊ)を直接使用することについての問題は、海面の行が天底から絶対変位であることであり、従って、このｊの関数は、天底から特定カメラでの最高位の海面の行までの変位に等しい一定バイアスを持つということである。
【００２６】
外側寄りのカメラは通常は、天底の列を含まない。従って、外側寄りのカメラは、隣接した内側のカメラを横切る海面の行と同程度に、行に大きな一定のバイアスを持っていてもよい。画像をゆがめないようにするために、海面の最上位の行が見つけられる列は０行にて置換され、海洋表面がより低い他の列は、「引き上げられて」、海面の最上位の行と一致する。
最上位の行は、ＬＩＤＡＲに最接近しているライダーの反射を表わす行、即ち図を見るユーザの視点である。行列インデックスの標準的な数学的な約束事が利用されるので、(０、０)は左上に位置し、即ち、より低位の行インデックスは、ＬＩＤＡＲに一層接近している反射を含んでいる。
【００２７】
特定のカメラについての海面の絶対変位の範囲は、ソフトウェア内の徹底的な探索によって見つかる。数学的には、それは以下に表現される。
Δｉextent＝Δｉ(ｊmax)−Δｉ(ｊmin)                          (５)
ここで、
【数２】
2008292449.tif 000004
である。
【００２８】
各列の行インデックスから引き出される実際の変位は、以下によって与えられる。
Δｉ(ｊ)←Δｉ(ｊ)−Δｉ(ｊmin)                                  (８)
【００２９】
このように、海面が最上位の列ｊminは全く動かない。他の全ての画素値は、以下の公式に従って移動する。
λij→λi-Δi(j)，ｊ，              （９）
(ここで、λijは行ｉ及び列ｊでの画像のグレイレベルであり、)ｉの値が増加しながら、反復されねばならない(データを重ね書きしないために)。これにより、画像内で列は上に引き上げられるように見え、海面は「平坦化された」ように見えるだろう。
【００３０】
航空機データプロセッサのメモリハードウエアが小規模の場合、全ての画素が移動するのではない。フレームの底部は、歪まない操作の後、「削除」される。この点から先は、データプロセッサは一層小さなフレームで動作する。
【００３１】
フレームが歪まない操作がされたので、海面を含む行は、歪んでいないフレームの垂直のシルエット(投影)の中で、最大値の行として容易に見つかる。
【数３】
2008292449.tif 000005
【００３２】
均等化は、先ず最初にそのインデックスが１０−５０１の範囲内の列にあって、正常に動作する(well-behaved)列(ＷＢＣ)を見つけることにより行なわれる。画像の縁の近傍の列は、副作用物、即ちノイズデータを有する。正常に動作する列は、４０００を越える値、又は所定の最小値未満の値を決して有さない。第１の制約は、列が不飽和であることを保証し、第２の制約は、列が現実的な表面反射を有する少なくとも１つの画素を有することを保証する。このように、正常に動作する列は飽和することなく、良好なダイナミックレンジを有する。これらの良い特性故に、正常に動作する列はまた、κlidar(後記する)を演算するのにも用いられる。
【００３３】
均等化のモデルは、各列が乗定数を除いて、同じ減衰曲線に追随することである。減衰曲線は指数関数的である必要はない。実際、指数関数曲線ではなくＷＢＣが用いられるのは正確にこの理由からである。即ち、あらゆる画素(海面下)の値のモデルは、以下である。
λij＝Ａjｆ(ｉ)，                         (11)
ここで、ｆ(ｉ)はＷＢＣの減衰曲線で、Ａjは海面でのカーブの値であり、
【数４】
2008292449.tif 000006
である。j'をＷＢＣのインデックスとする。減衰曲線は以下の如く、形成される。
【数５】
2008292449.tif 000007
【００３４】
このように、減衰曲線は、(海面にて)１つの最大値を有する一組の数であり、インデックスｉが増加するにつれ、ゼロに向かって減少する。
【００３５】
このモデルは、ＷＢＣを均等化するのに完璧である。このモデルが、各列ｊについて完璧であったならば、
λij／Ａjｆ(ｉ)                                      (13)
式(13)は目標を含む画素を除き、常に１つである。これはＷＢＣが目標を含まないと仮定する。勿論、これは現実の場合ではない。しかし、十分に現実に近づいている。水中の外乱又は受信器の副作用物は、均等化に関する問題に繋がる。このように、ストリーク管のブロブ、即ちＭＬＯ画素は正確には等しくないだろう。概して、それは検知器の役に立つ。
【００３６】
システム全体に於いて、目標がＷＢＣ内に時々あるという事実により、検知に失敗する可能性が、少し生じる。この問題は、ＷＢＣをテストして、目標がないことを確かめることにより克服される。これを行う多くの方法がある。
目標がＷＢＣの内にあると判れば、他の列はＷＢＣにあるように選ばれなけばならない。
【００３７】
上記のシステムは、Ａj＝Σλijとして修正され、ここで合計は海面下の行にて実行される。海面下の行のグレイレベルの平均を用いることは、１つのグレイレベルを使用するよりも安定している。
【００３８】
検知器は、データ領域(即ち、大凡ＭＬＯのサイズ)の統計値を計算し、該統計値をスレッショルドと比較することにより作動する。統計値がスレッショルドよりも大きければ、目標はデータのその地域に存在することが示される。統計値がスレッショルドよりも大きくなければ、目標は示されない。この手順は、目標の捜索にて、データの全地域に徹底的に適用される。目標が示されても、それはＭＬＯかＮＭＬＯにまだ分類されなければならない。
【００３９】
本発明は、ストリーク管／ＣＣＤカメラ用のデータモデルがポアソン分布(Poisson)として合理的にモデル化されると仮定する。その場合、非均一な強度の目標は、その十分な統計として重み付け合計を有している。
【数６】
2008292449.tif 000008
異なるデータモデルがデータに、より良く当てはまる場合、十分な計算容量があれば、上記の式は置き換えられてもよい。
【００４０】
計算を加速すべく、目標が均一な強度を持つように、検知器を設計する。その後、上記の式は以下の如く、単純化される。
【数７】
2008292449.tif 000009
【００４１】
この合計が、１の振幅の方形パルスで重畳されたものと、数学的に等価であることに注意されたい。換言すれば、「１を掛け、加える」は「加える」と等価である。十分な統計は、単に画素値の合計である。性能評価は、十分な統計及びポアソン見込みモデルから分析的に導かれる。図５は、(二次元)検知器の一次元の例である。
【００４２】
均一な強度を仮定することは、検知器が相関性受信器であることを意味する。受信器は、入力データの標準化されたバージョンを格納された矩形テンプレートと比較する。これは、格納されたテンプレート中の対応する画素と、データの地域の各画素をポイントごとに掛け合わせて合計することにより実行される。
【００４３】
次に、標準化された合計は、スレッショルドと比較される。スレッショルドは、連続する検知の「警戒」部分の偽警報割合(ＦＡＲ)を決定する。検知性能がシステムの要求を下回ることなく、偽警報割合ができるだけ低いように、スレッショルドが設定される。
【００４４】
ＡＬＭＤＳ検知器用の格納されたテンプレートは、１つのスーパーフレーム中の行列次元で数の矩形領域である。テンプレートを多次元空間のベクトルであると考えてもよい。
【００４５】
'１'と'０'を使用するテンプレートについて、過剰に単純化された見方は、図６に示される。図６は、機雷の単純なＮｘＮテンプレートである。影が付けられた画素は、目標('1')に相当し、また、明るい画素は背景('０')に相当する。各ボックスは１つの画素を表わす。
【００４６】
図７は、動作システムのより現実的なテンプレートの例を示す。図７は、矩形フィルタの下方範囲(行)対進路の交差の次元を示す。計算を促進するために１と０からなるが、それは、図６に示すものよりはるかに大きなテンプレートである。テンプレート次元はソフトウェアへの入力パラメータである。
【００４７】
以前に記載されたように、単純な２進法のテンプレート(１と０)を選択することにより、計算が促進され得る。これは、以下の２つの理由による:
(１)  各テンプレート位置にて、実行されなければならない乗算／加算は単なる「加算」であり、(２)  次のテンプレート位置に於ける合計は、図８に一次元で示すように、１のデータ値を加え、他のデータ値を引くことにより見出されるからである。図８は、演算を加速する一次元の例である。最初の加算が終了した後に、それに続く各合計は１つの数を加え他の数を引くことにより見出される。
【００４８】
二次元に於いて、この工程は以下の如く、実行される。図１０のスーパーフレームを考えると、部分和(図１１)のフレームは、全ての画素を加算することにより計算され、それらのインデックスは部分和のフレームの画素のインデックス未満又はインデックスに等しい。全ての加算は、矩形フィルタによる畳込み計算を加速するために前もって為されている。
部分和のフレームが作成される理由は、部分和は図９に示すように、わずか３つの加算しか要求しない単純な方法で組み合わせられるからであり、矩形フィルタの出力は(矩形Ａ＋矩形Ｄ)の和から(矩形Ｂ＋矩形Ｃ)を引いたものに等しい。
【００４９】
図１０及び図１１は、加算された画像を用いて２ｘ２矩形の整合フィルタの計算の実行例を示す。図１０に強調されるスーパーフレームの画素値の合計は、１８である。
この同じ値は、図９の規則に従って、図１１の部分和のフレームからの値を組み合わせることにより得られる。再び、値１８が得られる。
【００５０】
検知器を信号の強さではなく信号の形状に反応するようにするために、検知器へ送られる画像の全ての領域は単位エネルギーに標準化される。領域内のエネルギーは、Ｌ１ノルム法(L1 norm)によって概算され、換言すれば画素値の合計に接近する。
これにより、信号のエネルギーを計算するのに、正方形を合計するコストの高い計算ではなく、部分和の加算画像を再度使用することができる。このように、検知器の出力は実際には、２つの合計の比率である。部分和のフレームは各合計を計算するのに利用される。
図１２は、検知器の出力が外部の矩形の中で画素値の合計で除算された内部の矩形中の画素値の合計の比率であることを示す。
図１２では、小さな矩形(50)の内部和(整合フィルタの値)は、大きな矩形(52)(データ中のエネルギ)内の外部合計で除算される。
【００５１】
スレッショルドが深さに依存するとの意味に於いて、スレッショルドはまた一般化された。これは、以下のように実行される。各深さ(行)に於ける検知器の出力は、その深さに於ける検知器の平均値及び標準偏差を計算するために利用される。検知器の統計ｋは、その列に於いて、(検知器の出力χijから自身の行平均
【数８】
2008292449.tif 000010
を引いたものを)標準偏差σiで除したものとして、定義される。
【数９】
2008292449.tif 000011
【００５２】
それは実際には、スレッショルドと比較される検知器出力の統計ｋである。
標準的なスレッショルドｒbは４.５と同じくらい高いと経験的に決定され、検知の可能性を下げることなく、偽警報割合を著しく減じた。
【００５３】
検知器は、互いの近傍に多重ヒット(スレッショルド交差)を生じ得る。そのような事象は目標の近傍では一般的である。それが起こる場合、最も大きな統計ｋを有するヒットが保持される。各列には１つの目標だけが許されるように、アルゴリズムに規則が組み込まれる。これは、クラスタリング・アルゴリズムと呼ばれる。これは、レーダーが目標を透かして見ることが物理的に不可能で、各列が視野方向角度に正確に対応しているからである。
【００５４】
分類器は、目標を分類するのに外観、サイズ及びアスペクト比を使用する。分類器は、入力として二次元の領域(「ポステージスタンプ」として知られている)を採り、該領域はコンタクトに集中した(コンタクトキューブとして知られる)５０列、１０行の三次元５０フレームから構成されるテンソルである。全てのポステージスタンプは、５０列毎の５０フレームである。ポステージスタンプは、目標の周りの領域の５行を平均化することにより構成される。
【００５５】
外観(texture)は、クリューゲル空間のグレイレベルに依存する行列の簡略版を使用して計算される。サイズとアスペクト比は、レーベンバーグ−マルクワルト・アルゴリズムを使用して、目標の二次元の投影に適合する二次元のガウス強度表面のパラメーターから導き出される。
【００５６】
外観は、目標と偽警報を区別するために利用される主な特徴である。コンタクトが滑らかであるほど、そのコンタクトは一層確かにＭＬＯであろう。一例に於いて、外観は１２１の数からなるポステージスタンプ内の１１の列領域毎に１１のフレームについて測定される。
【００５７】
関連する外観の識別器は、クリューゲル空間のグレイレベルに依存する行列に依存する。
【数１０】
2008292449.tif 000012
ここで、Ｓ(ｉ，ｊ｜ｄ，θ)は、グレイレベルｉがグレイレベルｊに対して向かう回数であり、
【数１１】
2008292449.tif 000013
及び
【数１２】
2008292449.tif 000014
であれば、
【数１３】
2008292449.tif 000015
であり、且つ
【数１４】
2008292449.tif 000016
である。
【００５８】
Ｓ関数は、近傍の画素(距離ｄ内)中のグレイレベルの同時確率関数(joint probability function)に似ている。全てのこれから、１つの外観番号が、全ての方向の平均を採ることにより生成される。
【数１５】
2008292449.tif 000017
【００５９】
関連する外観が、深さに依存するスレッショルドよりも大きければ、コンタクトはＭＬＯであると示される。
【００６０】
海面の目標は、影一致検知(ＣＳＤ)として知られている付加的な処理段階を経る。図１３は、ＣＳＤアルゴリズムの重要点を記載するトップレベルのブロックダイヤグラムである。アルゴリズムは、信号(60)を生で検知し、影を検知し(62)、クラスタリングし(64)、目標を検知し(66)及び外観による識別(68)の工程を含んでいる。
【００６１】
信号を生で検知するには、表面層を画素ユニット内の深さｎ1及びｎ2の間に収容する。所定の方位角ｍについては、信号が最大であるとき、ｎ1≦ｎmax(ｍ)≦ｎ2を深さとする。これは「海面」を定義する。
【００６２】
「海面信号」はＳ(ｍ、ｎmax(ｍ)、κ)で表示され、ここでκはフレームのインデックスである。所定のスレッショルドＴrawについて、｛ｍp，ｐ＝１，２，３…｝を一組の方位角にし、そこではスレッショルドを越えて信号がピークになる。即ち、
Ｓ(ｍp、ｎmax(ｍp)、κ)≧Ｔraw                        (22)
【００６３】
この工程の出力は、方位角が検知されたピークの関数として、ピークが検知された深さである。
ｎmax(ｍp)(ｐ＝１，２，３…)                         (23)
【００６４】
工程の影検知部では、影の存在は、生の検知段階で見つかったピーク方位角｛ｍp，ｐ＝１，２，３…｝で調べられる。所定のｍpについて、影の確認は、方位角画素ｍpの回りに集中した３・(２Δｍ＋１)画素を用いて「影深さ」で行なわれる。
ｎshadow＝ｎmax(ｍp)＋ｎoffset，                    (24)
【００６５】
図１４は、△ｍ＝２の場合の影検知スキームを描く。影の確認工程の第１段階は、ローパスフィルタリングである。これは、単純な５点の三角形フィルタを用いて、信号を畳み込むことにより達成することができる。濾過された画素値は、以下によって与えられる。
【数１６】
2008292449.tif 000018
【００６６】
【数１７】
2008292449.tif 000019
式(26)を図１４に示す影深さｎshadowに於ける３番目、８番目及び１３番目の方位角に集中した３つの異なる方位角セクションに於ける濾過されたパワー合計にする。データは各端部にある、最初及び最後の２つの画素の反射である、２つの予備の画素によって増加する。換言すれば、
Ｓ(−１)＝Ｓ(２)，Ｓ(０)＝Ｓ(１)，Ｓ(16)＝Ｓ(15)及びＳ(17)＝Ｓ(14)   (27)
【００６７】
Ｐcenter＜Ｐleft  及び  Ｐcenter＞Ｐright                (28)
であれば、影は表現される。
【００６８】
この工程は、ｎoffsetから開始してＮ個の異なる深さで、繰り返される。その後、Ｎ個の影からＭ個を検知する基準が課せられて、調査(interrogation)中のピーク信号が影を有するか否かを決定する。
【００６９】
クラスタリング及び一致した影の検知については、影検知工程の出力はデータキューブであり、三次元中の影が検知された位置の組は(ｍs，ｎs，κ)である。
この例において、海面の目標に厳密に対処しているので、このデータキューブは深さ次元で崩壊し、クラスタリングは方位角でのみ行われる。表面下の目標検知が対象の場合、この手順は有効ではない。結果はデータシートである、即ち、方位角及び軌跡に沿う次元内にある影検出位置の二次元配列(ｍs,κ)である。
【００７０】
方位角「クラスタリング・ウィンドウ」のサイズをｍwとする。一致した影を検知するプリカーサーとして、ウィンドウ内のデータは方位角次元で崩壊する。
図１５に図示するように、クラスタ処理は進む。
【００７１】
図１６は、「スライディングウィンドウ」クラスタリング(単一ストライドを具えた)の結果を示す。図１６では、Ｍ番目の行は、Ｍ番目の方位角から(Ｍ＋ｍw−１)番目の方位角に延びるクラスタリングウィンドウ内に生じる検知の投影(projection)を示す。
【００７２】
一実施例に於いて、検知の基準は、連続した５つの影の最小値、又は７つの影の中から６つの最小値(８つの中の７つ、９つの中の８つの様でもよい)の何れかである。
それらの間に２つ以上の連続する影が見当たらない場合、２つの検知は、軌跡に沿う次元で個別であると考えられる。
【００７３】
この基準を満たす影検知のグループについては、方位角及び検知された対象物のフレームは、グループ内の方位角及び検知のフレーム数によって決定される。
【００７４】
差別化はシンチレーション計数を使用して実行され得る。目標信号の一時的な変化は、目標の中心方位角で検討される。一定数以上のシンチレーションが生じる場合、コンタクトはＭＬＯであると考えられない。これを表示するシンチレーション及びシンチレーションのスレッショルド番号は、経験的に決定される。
【００７５】
指数関数的減衰のモデルが仮定される場合、変数κlidarは水からのライダーエコーの減衰係数である。κlidarを評価するアルゴリズムの微分は、容積測定のリターンが水中の繰り返し往復の減衰に比例した強度を持っているとの観察から始まる。
【００７６】
測定されたリターンｆ(ｚ)は、
【数１８】
2008292449.tif 000020
ここでＧは、全てのレシーバーパラメーター及び水反射率を含む定数である。比(30)を生成してもよい。
【数１９】
2008292449.tif 000021
【００７７】
次に
【数２０】
2008292449.tif 000022
及び
【数２１】
2008292449.tif 000023
ここで測定の記述を単純化すべく、ｆk＝ｆ(ｚk)とする。２つの深さｚ1とｚ2が△ｚだけ離れた隣接した行であるならば、
κlidar＝ｌｎ(ｆ1／ｆ2)／(２△ｚ)                      (33)
【００７８】
更に、次のことに注目する。
【数２２】
2008292449.tif 000024
それで、
【数２３】
2008292449.tif 000025
【００７９】
別の記述の観点から、
【数２４】
2008292449.tif 000026
即ち、２つの隣接した行の容積測定のリターンの比は理想的には、定数である。いくつかの比率は平均化されて、より多くの安定した評価を作り出す。操作上、ソフトウェアでは、４行が平均化される、なぜならこれらの行に最も多くの有効数字があるからである。
【数２５】
2008292449.tif 000027
【００８０】
システム有効度インデックス(ＳＥＩ)は、探索された(任務中に計算された)水量及びデータ初期化中の十分な深さ性能を生成すること(ｋ′ｚ′)についての評価されたｋ値(ｋest)に基づいて計算される。
【００８１】
ＳＥＩは次のように計算される:
ＳＥＩ＝(ｋ′ｚ′)／ｋest)
システムは、更に十分な深さ性能を生成するためのデフォルト値を計算する。
【００８２】
そこで、
【数２６】
2008292449.tif 000028
ここで、Ｋdは拡散減衰の測定値であり、Ｋ′dは任務が計画された時の拡散減衰の予測値であり、ｚ′は深さ探知の仕様である。例えば、Ｋ′dは光学海洋学教科書中のテーブル、又は海軍ＭＥＤＡＬシステムから得られる海洋のその部分についての拡散減衰の値であるかも知れない。
【００８３】
ＳＥＩは距離を測定する。ＳＥＩPIDSは、以下のように開発された。特定の深さｚ′に於けるＡＴＲＰｄ及びＦＡＲ性能は、生成物Ｋ′dＺ′によって決定され、ここでＫ′dは予測された拡散減衰である。
【００８４】
拡散減衰の測定値Ｋdが予測値と異なれば、減衰物が一定に保たれる限り、新たな深さｚで同じ性能が達成される。即ち、
Ｋdｚ＝Ｋ′dｚ′                                     (39)
【００８５】
このように、同じＰｄ及びＦＡＲ性能が達成される新たな深さは、
【数２７】
2008292449.tif 000029
である。
【００８６】
このように、システムの有効性は現実的で、絶対的な深さであり、予想される拡散減衰Ｋ′dを用いて、特定の深さｚ′で予想されたものと、同じ性能が得られた。
【００８７】
発明は幾つかの実施例について記載されてきたが、以下の請求の範囲に示した発明の範囲から離れることなく、記載された実施例に種々の変更が為され得ることは当業者ならば明白であろう。
【図面の簡単な説明】
【００８８】
【図１】本発明の実施例に従って構成されたシステムの概略図である。
【図２】本発明の実施例に従って構成されたシステムのブロック図である。
【図３】本発明の方法を示すフロー図である。
【図４】画像フレームが如何にして結合されるかを示す概略図である。
【図５】一次元例のデータに適用される一次元矩形カーネルの一次元概略図である。
【図６】二次元カーネルの概略図である。
【図７】二次元カーネルの概略図である。
【図８】本発明に用いられるカーネルによる畳み込みが如何に更新されるかを示す一次元概略図である。
【図９】加重画像から整合フィルタのカーネル畳み込み値が如何に演算されるかを示す数値の図表である。
【図１０】加重画像から整合フィルタのカーネル畳み込み値が如何に演算されるかを示す数値の図表である。
【図１１】加重画像から整合フィルタのカーネル畳み込み値が如何に演算されるかを示す数値の図表である。
【図１２】ＡＬＭＤＳシステムにて、目標の特定の画像を検知するのに用いられるカーネルの概略図である。
【図１３】本発明の表面目標検知部の方法を示すフロー図である。
【図１４】本発明に用いられる表面目標データの概略図である。
【図１５】多数の画像フレームを横断する表面目標データの例の概略図である。
【図１６】多数の画像フレームを横断する表面目標データの例の概略図である。
【図１７】多数の画像フレームを横断する表面目標データの例の概略図である。
【図１】
2008292449.tif 000030
【図２】
2008292449.tif 000031
【図３】
2008292449.tif 000032
【図４】
2008292449.tif 000033
【図５】
2008292449.tif 000034
【図６】
2008292449.tif 000035
【図７】
2008292449.tif 000036
【図８】
2008292449.tif 000037
【図９】
2008292449.tif 000038
【図１０】
2008292449.tif 000039
【図１１】
2008292449.tif 000040
【図１２】
2008292449.tif 000041
【図１３】
2008292449.tif 000042
【図１４】
2008292449.tif 000043
【図１５】
2008292449.tif 000044
【図１６】
2008292449.tif 000045
【図１７】
2008292449.tif 000046
【外国語明細書】
2008292449000001.pdf