信号レベル判定装置及び信号レベル判定方法
インパルス性雑音を含む背景雑音が存在する中で信号をロバストに検出するための装置及び方法が開示される。信号の各値は、２次元空間内の直交軸上の２つの座標によって定義される半円上の点にそれぞれマッピングされる。変換された各点の２つの座標からそれぞれの平均が計算され、それぞれの２つの平均を用いて半円上の点の平均方向が計算される。平均方向は信号領域内に逆マッピングされ、検出閾値と比較される。検出閾値は、半円上の点の密度に基づいて設定されてもよい。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、背景雑音及びインパルス性干渉によって劣化した信号のレベルを判定するための方法及び装置に関する。また、限定はしないが、刻々と変化する荒れた海面によって後方散乱した干渉信号が存在する中で、小さな物体から反射される信号を検出するのに特に適している。
【背景技術】
【０００２】
数多くの実用的な用途において、概ねガウス特性を有する雑音（たとえば、熱雑音）とインパルス性の干渉との混合物によって、目的とする信号が劣化する。そのような組み合わせの確率分布は、多くの場合に所謂「ヘビー」テールを示すことになり、非ガウス現象を特徴付けるために種々の統計モデルが開発されている。たとえば、海面から反射されるマイクロ波信号の大きさは多くの場合に、ワイブル分布、対数正規分布、又はＫ分布によって特徴付けられる。
【０００３】
海洋環境において動作するマイクロ波センサは、海面反射と呼ばれることもある海面から反射される不要な信号が存在する中で、潜在的に目的とする種々の小さな物体を確実に検出することが期待される。検出される小さな物体は、ボート及びゴムボート、ブイ、種々の漂積物及び氷山の小さな破片を含む。それらの物体の中には、船舶が安全に航海するのを大きく脅かすものもあれば、救難任務、沿岸監視、自国保安等の目的となるものもある。
【０００４】
非ガウス海面反射は、多くの場合にガウス雑音に関して最適に動作するように設計されている数多くのセンサの検出性能に悪影響を及ぼす可能性がある。したがって、多年にわたって、海面反射がある中で小さな信号を検出するという問題に対して、種々の解決策が提案されてきた。
【０００５】
海面反射がある中で物体を検出することができる実用的な非コヒーレントシステムの代表的な例が、２００７年１０月２３日の米国特許第７，２８６，０７９号：「Method and Apparatus for Detecting Slow-Moving Target in High-Resolution Sea Clutter」において提示されている。
【０００６】
上記の開示によれば、各レンジセルからのレーダーエコー（radar return）を２値データのストリームに変換するために比較器が用いられる。反射された信号が所定の閾値を超える場合には、その信号は２値データの「１」によって表され、そうでない場合には、その信号は２値データの「０」によって表される。この動作が図１において概略的に示されている。そのような「硬」判定の１つの短所は明らかであり、閾値をわずかに下回る信号は全体的な検出判定に影響を及ぼす可能性があるにも関わらず破棄されるため、有用な情報が無視される。
【０００７】
次に、レンジエクステント（range-extent）フィルタが２値信号を処理し、隣接するレンジセル内に現れる「１」を含むクラスタの存在を示す。これらのクラスタは、レンジ内に延在する物体の存在を示しているものと見なされる。図２は、レンジエクステントフィルタリングの動作を示す。しかしながら、そのようなフィルタを効率的にすべき場合には、仮想される物体の広がりが予め分からなければならない。さらに、単一のレンジセルを占有する程度の小さい物体には、レンジエクステントフィルタリングは適用することができない。
【０００８】
最後に、仮想される物体の時間的な持続性が所定の時間間隔にわたって調べられる。適切な持続積分器（persistence integrator）の出力が予め選択された閾値を超える場合には、物体検出が宣言される。この動作が図３に概略的に示されている。残念なことに、時間における積分は、閾値比較又はハードリミッティングの過程において失われる情報を再生することはできない。
【０００９】
上記のシステムの検出性能が、S. D. Blunt、K. Gerlach、J. Heyerによる「HRR Detector for Slow-Moving Targets in Sea Clutter」（IEEE Transaction on Aerospace and Electronic Systems, vol. 43, Jul. 2007, pp. 965-975）においてさらに詳しく解析されている。
【００１０】
この刊行物は関連する理論的背景を検討し、さらには、検出器の実用的な実施態様に繋がる種々の仮定及び単純化も検討している。
【００１１】
たとえば、検出器によって利用される比較器のステップ状の特性は、理論的な解析から得られる最適な非線形特性の唯一の好都合かつ実用的な近似であることを示している。
【００１２】
さらに、図１ａの観測結果を線形に（算術平均と同様に）結合してその結果を閾値と比較しても、インパルス性の雑音を示す背景に対して必要とされるロバストな検出手順は得られないことを明らかにしている。
【発明の概要】
【発明が解決しようとする課題】
【００１３】
それゆえ、劣化した信号のレベルを非線形かつロバストに判定するための方法及び装置を提供することが望ましい。
【００１４】
特に、この方法及び装置は、従来の技法よりも効率的にスパイク性の海面反射がある中で小さな物体を検出できるようにすることが望ましいが、この目的に限定するものではない。
【課題を解決するための手段】
【００１５】
本発明によれば、信号値を処理してそれら信号値のレベルを判定する方法が提供される。この方法は、シグナルプロセッサが、
複数の信号値のそれぞれを、円の少なくとも一部の円弧上の各点にマッピングするプロセスと、
それらの点を定義する座標のうちの少なくとも１つの座標の値を用いて平均値を求めるプロセスと、
信号値が閾値を超えるか否かを判定するプロセスと
を実行することを含む。
ここで、上記の判定は、
平均値を閾値と比較すること、又は
平均値を信号値の領域内の値に逆マッピングし、その逆マッピングされた値を閾値と比較することによって実施される。
【００１６】
本願の発明者は、このようにして平均値を求めると、その平均値は、他のタイプの平均（たとえば、算術平均）よりも、１つ又は複数の大きな信号値（インパルス性の雑音スパイクから生じる大きな値等）によって及ぼされる影響がはるかに小さいために、信号レベルをよりロバストに推定するものになることを見出した。
【００１７】
その結果として、インパルス性の雑音を有する、雑音のある背景における信号の検出を改善することができる。
【００１８】
本発明はさらに、上記の方法を実行するための装置も提供する。
【図面の簡単な説明】
【００１９】
【図１】各レンジセルからのレーダーエコー（radar return）を２値信号のストリームに変換するための比較器の従来技術における使用法を示す図である。
【図２】クラスタの存在を示すために２値データにおいて実行されるレンジエクステントフィルタの従来技術における動作を示す図である。
【図３】持続積分器の従来技術における動作を示す図である。
【図４】本発明の一実施形態による、円弧平均を求めるステップを示す図である。
【図５】本発明の一実施形態に従って構築される円弧平均計算器ＣＭＣの機能ブロック図である。
【図６】本発明の第２の実施形態に従って構築される円弧平均計算器ＣＭＣの機能ブロック図である。
【図７】ａ）は、本発明の第３の実施形態に従って構築される２つのマッピング関数を示す図である。ｂ）は、信号値の角度位置への変換を示す図である。ｃ）は、信号値の単位半円上へのマッピングを示す図である。
【図８】本発明の第３の実施形態に従って構築される修正型円弧平均計算器ＭＣＭＣの機能ブロック図である。
【図９】信号検出のために円弧平均ＣＭ及び円弧密度の両方が用いられるときの判定領域を示す図である。
【図１０】本発明の第４の実施形態に従って構築される円弧平均計算器ＣＭＣを組み込んだ信号検出システムＳＤＴの機能ブロック図である。
【発明を実施するための形態】
【００２０】
インパルス性雑音がある中での信号検出を分析したところ、図１ａにおいて示されるような信号値の指標として従来の算術平均を利用して検出閾値と比較しても、最適な検出性能を与えることができないことがわかった。それゆえ、信号検出にさらに適した信号レベルを求めるための改善されたロバストな方法が必要とされる。
【００２１】
本発明の実施形態によれば、信号の１組の観測値（たとえば、包絡線、大きさ、或いは電力）の各値ｐが、円弧上の位置を定義する直交軸上の２つの座標（ｘ，ｙ）を有する２次元空間内の１つの点にそれぞれ変換される。この円弧は１８０度（他の単位においても等量）以下の角度を有する。この円弧は半円弧であることが好ましい。個々の平均は、変換された点のｘ座標及びｙ座標について個別に計算される。次に、結果として生成される２つの平均を用いて、２次元空間内の変換された点の平均方向θMDを計算する。その後、この平均方向θMDを用いて、観測された信号値を閾値と比較する。第１の技法では、この比較は、逆変換Ｍ-1を適用して平均方向θMDを観測値の信号領域にマッピングし直して、本明細書において円弧平均ｐCMと呼ばれる量を与え、この円弧平均を信号領域において閾値と比較することによって実行される。その代わりとなるが同等の第２の技法では、この比較は、平均方向θMDを半円の２次元空間において閾値と比較することによって実行される。この閾値は、最初からその２次元空間内で定義してもよいし、或いは、観測値の信号領域において定義して、信号値を変換するために用いられるのと同じ変換を用いて半円の２次元空間内に変換してもよい。
【００２２】
代替の実施形態では、信号の１組の観測値ｐの各値が、２次元空間内の円弧上の点にそれぞれ変換される。２つの座標ｘ、ｙのうちの一方（両方ではない）の値から平均が計算される。その後、上記の２つの技法のうちの１つを用いて、観測される信号値を閾値と比較する際に、この平均が用いられる。
【００２３】
第１の実施形態
本発明の第１の実施形態によれば、以下の３つのステップ（図４においても概略的に示されている）を含む手順を適用することによって、円弧平均ｐCMが計算される。
【００２４】
ステップ１
信号の観測値ｐを、マッピング関数Ｍ（ｐ）を用いて単位半円上にマッピングする。その結果、目的とする所定の範囲（ＰＬ，ＰＨ）からの各値ｐkは、角度位置θkで単位半円上に配置される対応する点によって表される。
【００２５】
したがって、信号のＫ個の観測値
【００２６】
【数１】
【００２７】
は、以下の対応する１組のＫ個の角度によって表される。
【００２８】
【数２】
【００２９】
ただし、αは任意の初期角度である。
【００３０】
ステップ２
半円の２次元空間内にあるマッピングされた各点ｐkのｘ座標及びｙ座標を求めるために、角度｛θk｝の正弦及び余弦を計算する。その後、それらの計算値を別々に平均して２つの個々の平均値を求める。
【００３１】
【数３】
【００３２】
その後、以下の式から、１組の角度における平均方向θMDを求める。
【００３３】
【数４】
【００３４】
ステップ３
平均方向θMDに逆マッピングＭ-1（θ）を適用し、観測された信号値ｐの信号領域内に変換し直すことによって、以下に示す１組のＫ個の元になる信号値の円弧平均ｐCMを得る。
【００３５】
【数５】
【００３６】
続いて、円弧平均ｐCMを信号領域において閾値と比較する。（代替的には、方向θMDを半円の２次元空間において閾値と比較する）。
【００３７】
図４は、所与の１組の信号値の円弧平均ｐCMを求める手順を概略的に示している。この事例では、図示される半円は最初の２つの象限に現れるが（それゆえ、α＝０）、この特定の特徴は不可欠ではない。
【００３８】
第１の実施形態とその利点を理解するのを容易にするために、説明に役に立つ簡単な例を以下に与える。
【００３９】
第１のステップは、ｐの測定値を、（０，１８０）度の区間に入る値ｐ*に変換することである。ｐの測定値を以下のように仮定する。
【００４０】
【数６】
【００４１】
線形シフト演算によって各値から２００が減算され、以下の値が得られる。
【００４２】
【数７】
【００４３】
次に、線形スケーリング演算によって各値が２で除算され、以下の値が得られる。
【００４４】
【数８】
【００４５】
このとき、これらの変換値ｐ*は、（０，１８０）の区間に入っている。変換値ｐ*とそれらに対応する角度位置θ（度単位で表されている）とは、同じ数値を有する。
【００４６】
データをシフト及びスケーリングするためのこの方法は、この演算のために用いることができる数多くの方法の一例にすぎないことに留意されたい。
【００４７】
上記の１組の変換値ｐ*の算術平均ｐAMは、以下の通りである。
【００４８】
【数９】
【００４９】
上記の式から明らかなように、単一の相対的に大きな観測値ｐ6*＝１５０（たとえば、干渉する雑音スパイク）が、残りの５個の観測値ｐ1*〜ｐ5*の平均（３０）を大きくシフトする。
【００５０】
この実施形態における処理の第２のステップでは、変換値ｐ1*〜ｐ6*の平均方向θMDが、以下の式から求められる。
【００５１】
【数１０】
【００５２】
それゆえ、円弧領域内における算術平均ｐAMは５０に等しかったが、観測値の平均方向θMDは４１に等しい。平均方向θMDは、単一のより大きな観測値（すなわち、ｐ6*）によって及ぼされる影響がはるかに小さいので、値ｐの平均レベルのよりロバストな推定量である。
【００５３】
この処理の第３のステップにおいて、本実施形態では、平均方向θMD＝４１度を、逆マッピングＭ-1を用いて変換して測定信号ｐと同じ領域に戻す。
【００５４】
それゆえ、変換された値ｐCM＝（４１×２）＋２００＝２８２である。
【００５５】
一方、信号領域内の信号値ｐの算術平均は、［（５×２６０）＋５００］／６＝３００である。
【００５６】
したがって、本実施形態によって生成される円弧平均ｐCMは、従来の算術平均よりも信号領域内の値ｐの平均レベルをよりロバストに示している。これは、円弧平均ｐCMの値のほうが、算術平均よりも誤った信号値ｐ6によって及ぼされる影響が小さいためである。
【００５７】
ｐCMの上記の計算を実行するための装置が図５に示されている。
【００５８】
図５に示される計算器は、以下の機能ユニットを備える。
− シフト及びスケールユニット５０２；
− ２つの非線形変換器ＳＮＬ５０３及びＣＮＬ５０４；
− ２つのタップ付き遅延線ＤＬＳ５０５及びＤＬＣ５０６；
− ２つの平均化回路ＡＶＳ５０７及びＡＶＣ５０８；
− 演算ユニットＡＴＵ５０９。
【００５９】
包絡線、大きさ、或いは電力のような正の入力信号ｐｐがシフト及びスケールユニット５０２に通されると、ユニット５０２は、入力信号の各値を１８０度にわたる範囲内の値のうちの１つにマッピングする。それにより、正規化された信号ＮＰが生成される。
【００６０】
正規化された信号レベルＮＰは、２つの非線形変換器ＳＮＬ５０３及びＣＮＬ５０４に並列に加えられる。変換器の出力ＳＳ及びＣＣは、共通の入力値NPから２つの適切に選択されたマッピング関数を利用することによって求められる。検討している事例では、以下の通りである。
【００６１】
【数１１】
【００６２】
先に述べたように、値ＳＳ及びＣＣは、２次元空間内の値ＮＰのｘ座標及びｙ座標を定義する。
【００６３】
信号ＳＳ及びＣＣは、それぞれのタップ付き遅延線ＤＬＳ５０５及びＤＬＣ５０６に沿って伝搬する：信号ＳＳのＫ個のサンプルを遅延セルＳ１、Ｓ２、．．．、ＳＫの出力において取得することができ、一方、信号ＣＣのＫ個のサンプルを遅延セルＣ１、Ｃ２、．．．、ＣＫの出力において取得することができる。
【００６４】
平均化回路ＡＶＳ５０７は、その出力において、セルＳ１、Ｓ２、．．．、ＳＫから得られる入力値の和に比例する値ＡＳを生成する。同様に、平均化回路ＡＶＣ５０８は、その出力において、セルＣ１、Ｃ２、．．．、ＣＫから得られる入力値の和に比例する値ＡＣを生成する。これら２つの値ＡＳ及びＡＣは、演算ユニットＡＴＵ５０９によって用いられ、円弧平均ＣＭが求められる。
【００６５】
その後、比較ユニット（図５には示されていない）において、計算された円弧平均ＣＭが閾値と比較される。
【００６６】
第２の実施形態
上記の第１の実施形態では、信号値ｐは「シフト及びスケール」タイプの線形演算を利用することによって角度値θに変換された。しかしながら、実用的な用途では、ダイナミックレンジを非線形に調整するために、観測された信号値に対して最初に非線形（たとえば対数）変換を適用した後、そのように変換された値を単位半円にマッピングすることが好都合である場合がある。そのような処理を実行する実施形態が以下に示される。
【００６７】
たとえば、有用な非線形マッピングは以下の形のものである。
【００６８】
【数１２】
【００６９】
ただし、クリップ関数Ｈ（・）は、その引数の最小値と最大値をそれぞれ−π／２とπ／２に制限し、γは、処理される信号のダイナミックレンジをさらに調整するために用いられるスケーリング係数である。
【００７０】
たとえば、観測値ｐの範囲が０．０１から１００までの場合には、γ＝π／４がｐの全ての値を（−π／２，π／２）内に収めることになる。ただし、両端も範囲に含む。０．０１〜１００の範囲外であるｐの新たな値が後に検出された場合に、その新たなｐ値に対してもγの同じ値が用いられるようにするために、クリップ関数Ｈには、ｐの新たなマッピング値を（−π／２，π／２）の範囲に制限することが要求される。実際には、ｐ値が大抵の場合に０．０１〜１００の範囲に入っており、これらの限界値に達するか又は限界値を超えることがめったにない場合、クリップ関数Ｈは通常は用いられない。
【００７１】
以下の理由から、半円上にマッピングする前に対数変換を適用することが、他の非線形変換を用いることよりも好ましい場合がある。
１．対数変換は圧縮性があり、大きい観測値ほど抑制する。
２．除算が単純な減算に変わるので、元になるデータのスケールパラメータがシフトパラメータに変換され、それにより、信号正規化の演算が簡略化される。
【００７２】
図６は、本発明の第２の実施形態に従って構築される円弧平均計算器ＣＭＣの機能ブロック図である。
【００７３】
この計算器は、以下の機能ユニットを備える。
− 対数変換器ＬＧＩ６０１；
− 減算器ＳＢＴ６０２；
− ２つの非線形変換器ＳＮＬ５０３及びＣＮＬ５０４；
− ２つのタップ付き遅延線ＤＬＳ５０５及びＤＬＣ５０６；
− ２つの平均化回路ＡＶＳ５０７及びＡＶＣ５０８；
− 演算ユニットＡＴＵ６０９。
【００７４】
したがって、図５に示される第１の実施形態の計算器と比べると、シフト及びスケールユニット５０２の代わりに、対数変換器ＬＧＩ６０１及び減算器ＳＢＴ６０２が用いられている。他の全ての構成要素は同じであるが、ＡＴＵ６０９は例外であり、ｐに対して実行される対数演算を考慮するために、第１の実施形態と比べて異なる逆マッピングＭ-1を実行する。
【００７５】
包絡線、大きさ、或いは電力のような正の入力信号ＰＰが対数変換器ＬＧＩ６０１に通されると、入力信号ＰＰのレベルの対数指標である信号LＰが生成される。すなわち、
【００７６】
【数１３】
【００７７】
次に、減算器ＳＢＴ６０２において、信号ＬＰから目的とする或る基準レベルＢＧの対数ＢＬが減算されることによって、信号ＬＰが正規化される。
【００７８】
【数１４】
【００７９】
たとえば、信号検出問題では、基準レベルは長期の観測値から得られる背景雑音の平均レベルとすることができる。ＬＰから値ＢＬを減算する操作は、ＬＰの値を（−π／２，π／２）の角度範囲を有する単位半円上の点にマッピングする。これは、係数γによって与えられるスケーリングと等価である。
【００８０】
正規化された信号レベルＮＰは、２つの非線形変換器ＳＮＬ５０３及びＣＮＬ５０４に並列に加えられる。変換器の出力ＳＳ及びＣＣは、共通の入力値NPから２つの適切に選択されたマッピング関数を利用することによって求められる。本実施形態では、以下の通りである。
【００８１】
【数１５】
【００８２】
信号ＳＳ及びＣＣは、それぞれのタップ付き遅延線ＤＬＳ５０５及びＤＬＣ５０６に沿って伝搬する：信号ＳＳのＫ個のサンプルを遅延セルＳ１、Ｓ２、．．．、ＳＫの出力において取得することができ、一方、信号ＣＣのＫ個のサンプルを遅延セルＣ１、Ｃ２、．．．、ＣＫの出力において取得することができる。
【００８３】
平均化回路ＡＶＳ５０７は、その出力において、セルＳ１、Ｓ２、．．．、ＳＫから得られる入力値の和に比例する値ＡＳを生成する。同様に、平均化回路ＡＶＣ５０８は、その出力において、セルＣ１、Ｃ２、．．．、ＣＫから得られる入力値の和に比例する値ＡＣを生成する。これら２つの値ＡＳ及びＡＣは、演算ユニットＡＴＵ６０９によって用いられ、円弧平均ＣＭが求められる。
【００８４】
その後、比較ユニット（図６には示されていない）において、計算された円弧平均が閾値と比較される。
【００８５】
第３の実施形態
上記の実施形態では、信号値ｐを単位半円上の点にマッピングした後、三角演算子を適用して点の位置を定義する半円の２次元空間内の２つの座標を求める。続いて、これらの座標を用いて平均方向θMDを計算し、必要に応じて円弧平均ｐCMを計算する。しかしながら、信号値ｐの半円への初期マッピングを、結果として生成される半円上の点の２つの座標を直接与えるように実行してもよい。したがって、第１及び第２の実施形態における三角演算を実行することによって座標を計算する必要はない。
【００８６】
そのような処理を実行する一実施形態が以下に説明される。
【００８７】
一般的に、信号値ｐの半円へのマッピングは、適切に構築された２つのマッピング関数Ｓ（ｐ）及びＣ（ｐ）を使用することによって実行することができる。マッピングは半円を生成することを要求されるので、それらのマッピング関数は以下の条件を満たさなければならない。
【００８８】
【数１６】
【００８９】
ただし、ｂは定数である。単位半円の場合、ｂ＝１である。
【００９０】
本発明の第３の実施形態によれば、第１のマッピング関数Ｓ（ｐ）は単調に増加する連続関数になる。この関数は、最小信号レベルＰＬの場合に最小値である−１を有し、最大信号レベルＰＨにおいて＋１に等しい最大値に達すると仮定する。そのような関数は無限に存在し、正弦波の適切に選択された一部分は数ある選択肢のうちの１つであるが、別の有用な関数が以下において検討される。
【００９１】
次に、第２のマッピングＣ（ｐ）が、マッピング関数Ｓ（ｐ）から以下のようにして得られる。
【００９２】
【数１７】
【００９３】
それゆえ、２つのマッピング関数Ｓ（ｐ）及びＣ（ｐ）を利用するマッピング演算Ｍ［Ｓ（ｐ），Ｃ（ｐ）］は、１次元のｐ空間から２次元の（Ｓ，Ｃ）空間への非線形変換と見なすことができる。結果として、各信号値ｐは、その点の一方の座標がＳ（ｐ）によって定義され、その点の他方の座標がＣ（ｐ）によって定義されるように、単位半円上の点に直接マッピングされる。その後、Ｓ（ｐ）及びＣ（ｐ）に算術平均を直接適用して、平均方向のためのパラメータを求めることができる。
【００９４】
より詳細には、平均方向θMDは、以下のように計算される。
【００９５】
【数１８】
【００９６】
その後、平均方向θMDを閾値と比較することができる。代替的には、必要であれば平均方向θMDに逆マッピングＭ-1（θ）を適用することによって、信号値ｐの円弧平均ｐCMが求められ、この円弧平均ｐCMが閾値と比較される。逆マッピングＭ-1は、Ｍ-1を定義する数学的な関数を計算して評価することによって、又は数値技法（反復法等）を用いて要求される値を求めることによって、求めることができる。
【００９７】
双曲線関数ｔａｎｈ ｗ及び１／（ｃｏｓｈ ｗ）は、マッピング関数Ｓ（ｐ）及びＣ（ｐ）として有利に利用することができる。
【００９８】
より詳細には、
【００９９】
【数１９】
【０１００】
であるので、マッピング［ｔａｎｈ ｗ，１／（ｃｏｓｈ ｗ）］は、値ｗを表す点を単位半円上に置く。
【０１０１】
正規化するために、双曲線変換ｔａｎｈ ｗ及び１／（ｃｏｓｈ ｗ）を実行する前に、対数変換を使用してｗ＝ｌｎ（ｐ）とすることが好都合であろう。この場合、円弧平均計算器の機能ブロック図は、図６に示される機能ブロック図と類似の構造を有する。しかしながら、この場合、非線形変換器ＳＮＬ５０３及びＣＮＬ５０４は、以下のマッピング演算を実行する。
【０１０２】
【数２０】
【０１０３】
対数変換と双曲線変換とを組み合わせることによって、円弧平均計算器の構造を単純化することもできる。そのような単純化は以下のように行うことができる。
【０１０４】
【数２１】
【０１０５】
であるので、双曲線変換ｔａｎｈ（ｗ）に先立って対数変換ｗ＝ｌｎ（ｐ）が行なわれる場合には、以下の式が成り立つ。
【０１０６】
【数２２】
【０１０７】
図７ａ及び図７ｂに示されるように、２つのマッピング関数の形状は、信号値ｐと、単位半円上の値を表す点の対応する角度位置θとの間の関係を表す。この関係は非線形であり、「ソフト」リミッタのその形は、以下の数学的な表現から生じる。
【０１０８】
【数２３】
【０１０９】
図７ｃは、信号のいくつかの選択された元になる値ｐに対応する角度位置θを有する単位半円を示す。図に示されるように、信号の値が大きいほど（ｐ＞１０）、対応する角度位置θが限界値π／２に近いクラスタを形成する。その一方で、１よりも小さいｐの値は、象限全体（−π／２，０）を占有する角度位置θを生成する。
【０１１０】
観測値｛ｐ1，ｐ2，．．．，ｐK-1，ｐK｝からパラメータＹK及びＸKが、
【０１１１】
【数２４】
【０１１２】
に従って求められると、信号値ｐの円弧平均ｐCMは、以下の式から得られる。
【０１１３】
【数２５】
【０１１４】
図８は、本発明の第３の実施形態に従って構築される修正型円弧平均計算器ＭＣＭＣ８００の機能ブロック図である。この構成では、図６の円弧平均計算器６００において用いられる対数変換器ＬＧＩ６０１及び減算器ＳＢＴ６０２が、可変利得増幅器ＶＧＩ８０１によって置き換えられている。この増幅器は、目的とする基準レベルを示す信号ＢＧに応じてその利得を調整する。
【０１１５】
この構成では、変換器ＳＮＬ８０２及びＣＮＬ８０３の出力ＳＳ及びＣＣは、共通入力ＣＰから、以下のようにして得られる。
【０１１６】
【数２６】
【０１１７】
修正型システムによって実行される他の機能及び演算は、図６のシステムによって実行されるのと同様であり、したがって、ここでは再び説明しない。
【０１１８】
本発明のこの実施形態に従って構築される円弧平均計算器が対象とする応用例のうちの１つは、背景雑音がある中での信号の検出である。以下の例は、この実施形態のそのような応用例を示す。この例では、対数関数及び双曲線関数がマッピングのために利用される。
【０１１９】
単位電力熱雑音と、時折、雑音レベルの１０倍を超えるスパイクが生じるインパルス性干渉とを含む雑音のある背景内において、単位電力のランダム信号が検出されると仮定する。また、例示のために、検出閾値が１．９に設定されていると仮定する。
【０１２０】
信号が無い場合に、背景雑音の６個の観測値が以下のようであると仮定する。
【０１２１】
【数２７】
【０１２２】
５個の等しいサンプルは熱雑音レベルを表すものであり、サンプル番号６はインパルス性干渉によって生成されたものであろう。この場合、算術平均ＰAMは２．５に等しい。それゆえ、値２．５は１．９の検出閾値よりも大きいので、算術平均ＰAMが検出統計値として用いられる場合、誤検出（false alarm）が宣言されるであろう。
【０１２３】
ここで、算術平均ではなく円弧平均ｐCMを用いて背景雑音のレベルを判定すると仮定する。パラメータＹ6及びＸ6は以下のようになる。
【０１２４】
【数２８】
【０１２５】
それゆえ、円弧平均ｐCMは以下のようになる。
【０１２６】
【数２９】
【０１２７】
上記より、この場合、円弧平均ｐCMは検出閾値を超えないので、誤検出は生じない。
【０１２８】
ここで、検出されるべきランダム信号が存在し、この信号＋雑音の事例においては６個の観測値が全て等しく、たとえばｐ1＝ｐ2＝．．．＝ｐ6＝２であると仮定する。この場合、算術平均ｐAMと円弧平均ｐCMはともに等しく、すなわちｐAM＝ｐCM＝２であり、いずれも正しい検出判定に繋がることは明らかである。
【０１２９】
したがって、円弧平均を検出統計値として利用する主な利点は、時折生じるより大きな観測値を減衰させる能力から得られる。一般的に、そのような特性はインパルス性の雑音によって劣化した信号を処理する際に非常に有用である。
【０１３０】
第４の実施形態
本発明のさらなる実施形態によれば、単位半円上の点の円弧密度ＤCCが計算され、それを用いて円弧平均ｐCMと比較される検出閾値が設定される。かつ／又は、円弧密度は、検出閾値と比較される以前の円弧平均の値の調整に用いられる。
【０１３１】
信号レベルｐの値｛ｐk｝を表す点｛θk｝の円弧密度ＤCCは、以下の式から求められる。
【０１３２】
【数３０】
【０１３３】
ただし、
【０１３４】
【数３１】
【０１３５】
であり、Ｓ（ｐk）及びＣ（ｐk）は、マッピングのために用いられる２つの関数である。
【０１３６】
値ＸK及びＹKは、第３の実施形態において信号値ｐの半円上への初期マッピングによって出力されるので、第３の実施形態はこれらの値ＸK及びＹKを計算するのに最も適している。しかしながら、第１及び第２の実施形態の場合のように、ｐを半円上にマッピングした後に、ＸK及びＹKを計算してもよいことは、当業者には明らかであろう。
【０１３７】
点｛θk｝が半円上の密集したクラスタを形成するとき、円弧密度ＤCCの値は１に近くなる。しかしながら、点｛θk｝が、たとえば支配的な雑音に起因して広範囲に分散するとき、円弧密度ＤCCの値は著しく小さくなるであろう。したがって、円弧密度ＤCCの値を円弧平均と共に用いて、信号検出をさらに改善することができる。
【０１３８】
図９は、円弧平均ＣＭを円弧密度ＤＣと適切に組み合わせて検出判定領域を構築するときの検出性能の潜在的な改善を示す。
【０１３９】
検出のために円弧平均ＣＭだけが利用されるとき、円弧密度にかかわらず、円弧平均の観測値が判定閾値ＴＨを超える場合に、信号の存在が宣言される。この場合、判定領域Ｄ１は、ラインＴＨから延在する長方形の形状を有する。
【０１４０】
しかしながら、円弧平均ＣＭが円弧密度ＤＣと共に用いられるとき、結果として生成される判定領域は、領域Ｄ２だけ拡大される。円弧密度値が或る所定の閾値Ｃ１を超えるとき、円弧平均のための判定閾値ＴＨは新たな下側値Ｔ１まで徐々に下げられるであろう。明らかに、面積（Ｄ１＋Ｄ２）は単独の面積Ｄ１よりも大きく、閾値ＴＨが判定境界ＤＢによって置き換えられると、改善された検出性能が達成される。
【０１４１】
図１０は、円弧密度を用いて検出閾値を設定する、本発明の第４の実施形態に従って構築される信号検出器ＳＤＴ１０００の機能ブロック図である。信号検出器ＳＤＴ１０００の目的とする用途は、海面反射がある中での小さな物体の検出である。背景雑音、又は信号＋雑音の観測サンプルＰＰが、円弧平均計算器ＣＭＣ１００１において処理される。背景雑音の平均レベルを示す補助信号ＢＬを利用することによって、サンプルレベルが正規化される。そのような信号は、たとえば、被試験セルに隣接する複数のレンジセルによって生成される観測値を長い時間間隔にわたって取り込んだものを平均することによって得ることができる。
【０１４２】
円弧平均計算器ＣＭＣ１００１は、それぞれ円弧平均及び円弧密度を示している２つの出力ＣＭ及びＤＣを提供する。これらの信号は、判定ブロックＤＥＴ１０１１によって用いられ、観測サンプルＰＰが海面反射だけによって生成されたのか、それとも海面反射内に埋もれた物体によって生成されたのかが判定される。そのために、判定境界を定義する信号ＤＢが、判定ブロックＤＥＴ１０１１の入力ＤＢに加えられる。このブロックの出力信号ＳＤは、海面反射内の信号の存否に関する包括的な判定である。
【０１４３】
上記の実施形態では、円弧密度ＤCCは、円弧平均ｐCMと比較される閾値を調整するために用いられることに留意されたい（図９に示されている）。しかしながら、これは、円弧密度ＤCCに基づいて円弧平均ｐCMを調整し、調整された平均を無変更の閾値と比較することと等価である。したがって、この実施形態では、閾値の代わりにｐCMを調整してもよい。同様に、円弧密度ＤCCに基づいて閾値及び円弧平均ｐCMの両方を調整してもよい。
【０１４４】
変更及び変形
本発明の好ましい実施形態のこれまでの説明は、例示及び説明のために提示されてきたものであり、本発明を網羅的に述べることも、開示されるのと全く同じ形に限定することも意図していない。これまでの説明に鑑みて、数多くの改変、変更、及び変形によって、当業者は本発明を検討される特定の用途に合わせて種々の実施形態において利用できるようになるのは明らかである。
【０１４５】
たとえば、上記の実施形態では、信号値ｐが半円上の点にマッピングされる。しかしながら、半円全体を用いることは不可欠ではなく、代わりに、半円よりも小さな円弧を用いることもできる。より詳細には、信号値ｐは、半円弧によって定められる角度、すなわちπラジアン（又は同じく、たとえば１８０度）以下の角度を有する任意の円弧にマッピングすることができる。
【０１４６】
上記の全ての実施形態において、ＸK（すなわち、半円上の点のｘ座標の平均）及びＹK（すなわち、半円上の点のｙ座標の平均）の両方を求めた後、それら両方の値を用いてθMDを計算することによって、平均方向θMDが計算される。しかしながら、ロバストではないが、ＸK及びＹKの両方を求める必要はなく、平均方向θMDを計算する代わりに、代替の平均ＭALTを以下のように計算することもできる。
【０１４７】
【数３２】
【０１４８】
より詳細には、点ｐが角度範囲０度〜１８０度又は１８０度〜３６０度を含む単位半円にマッピングされるとき、ＸKの計算に用いられる余弦値がこの範囲内で一意の値を有するので、ＭALT＝ＸKが用いられる。すなわち、半円上の点ｐの値とｘ軸との間に１：１のマッピングが存在する。一方、点ｐが角度範囲９０度〜２７０度又は２７０度〜９０度を含む単位半円にマッピングされるとき、ＹKの計算に用いられる正弦値がこの範囲内で一意の値を有するので、ＭALT＝ＹKが用いられる。
【０１４９】
その後、代替の平均ＭALTを閾値と比較することもできるし、或いは必要に応じて信号領域にマッピングし直した後に、閾値と比較することもできる。たとえば、第１及び第２の実施形態では、このようにマッピングし直すことは、逆余弦マッピング（ＭALT＝ＸKの場合）又は逆正弦マッピング（ＭALT＝ＹKの場合）を用いて半円上の点にマッピングし直した後、逆マッピングＭ-1を適用することによって実行することができる。
【０１５０】
上記の第４の実施形態では、円弧平均ｐCMが、円弧密度ＤCCに基づいて設定される閾値と比較される（又は、円弧平均ｐCMが円弧密度ＤCCに基づいて調整される）。しかしながら、代わりに、平均方向θMDが、円弧密度ＤCCに基づいて設定される閾値と比較されてもよい（又は、平均方向θMDが円弧密度ＤCCに基づいて調整されてもよい）。
【０１５１】
当然のことながら、他の変更も可能である。
【技術分野】
【０００１】
本発明は、背景雑音及びインパルス性干渉によって劣化した信号のレベルを判定するための方法及び装置に関する。また、限定はしないが、刻々と変化する荒れた海面によって後方散乱した干渉信号が存在する中で、小さな物体から反射される信号を検出するのに特に適している。
【背景技術】
【０００２】
数多くの実用的な用途において、概ねガウス特性を有する雑音（たとえば、熱雑音）とインパルス性の干渉との混合物によって、目的とする信号が劣化する。そのような組み合わせの確率分布は、多くの場合に所謂「ヘビー」テールを示すことになり、非ガウス現象を特徴付けるために種々の統計モデルが開発されている。たとえば、海面から反射されるマイクロ波信号の大きさは多くの場合に、ワイブル分布、対数正規分布、又はＫ分布によって特徴付けられる。
【０００３】
海洋環境において動作するマイクロ波センサは、海面反射と呼ばれることもある海面から反射される不要な信号が存在する中で、潜在的に目的とする種々の小さな物体を確実に検出することが期待される。検出される小さな物体は、ボート及びゴムボート、ブイ、種々の漂積物及び氷山の小さな破片を含む。それらの物体の中には、船舶が安全に航海するのを大きく脅かすものもあれば、救難任務、沿岸監視、自国保安等の目的となるものもある。
【０００４】
非ガウス海面反射は、多くの場合にガウス雑音に関して最適に動作するように設計されている数多くのセンサの検出性能に悪影響を及ぼす可能性がある。したがって、多年にわたって、海面反射がある中で小さな信号を検出するという問題に対して、種々の解決策が提案されてきた。
【０００５】
海面反射がある中で物体を検出することができる実用的な非コヒーレントシステムの代表的な例が、２００７年１０月２３日の米国特許第７，２８６，０７９号：「Method and Apparatus for Detecting Slow-Moving Target in High-Resolution Sea Clutter」において提示されている。
【０００６】
上記の開示によれば、各レンジセルからのレーダーエコー（radar return）を２値データのストリームに変換するために比較器が用いられる。反射された信号が所定の閾値を超える場合には、その信号は２値データの「１」によって表され、そうでない場合には、その信号は２値データの「０」によって表される。この動作が図１において概略的に示されている。そのような「硬」判定の１つの短所は明らかであり、閾値をわずかに下回る信号は全体的な検出判定に影響を及ぼす可能性があるにも関わらず破棄されるため、有用な情報が無視される。
【０００７】
次に、レンジエクステント（range-extent）フィルタが２値信号を処理し、隣接するレンジセル内に現れる「１」を含むクラスタの存在を示す。これらのクラスタは、レンジ内に延在する物体の存在を示しているものと見なされる。図２は、レンジエクステントフィルタリングの動作を示す。しかしながら、そのようなフィルタを効率的にすべき場合には、仮想される物体の広がりが予め分からなければならない。さらに、単一のレンジセルを占有する程度の小さい物体には、レンジエクステントフィルタリングは適用することができない。
【０００８】
最後に、仮想される物体の時間的な持続性が所定の時間間隔にわたって調べられる。適切な持続積分器（persistence integrator）の出力が予め選択された閾値を超える場合には、物体検出が宣言される。この動作が図３に概略的に示されている。残念なことに、時間における積分は、閾値比較又はハードリミッティングの過程において失われる情報を再生することはできない。
【０００９】
上記のシステムの検出性能が、S. D. Blunt、K. Gerlach、J. Heyerによる「HRR Detector for Slow-Moving Targets in Sea Clutter」（IEEE Transaction on Aerospace and Electronic Systems, vol. 43, Jul. 2007, pp. 965-975）においてさらに詳しく解析されている。
【００１０】
この刊行物は関連する理論的背景を検討し、さらには、検出器の実用的な実施態様に繋がる種々の仮定及び単純化も検討している。
【００１１】
たとえば、検出器によって利用される比較器のステップ状の特性は、理論的な解析から得られる最適な非線形特性の唯一の好都合かつ実用的な近似であることを示している。
【００１２】
さらに、図１ａの観測結果を線形に（算術平均と同様に）結合してその結果を閾値と比較しても、インパルス性の雑音を示す背景に対して必要とされるロバストな検出手順は得られないことを明らかにしている。
【発明の概要】
【発明が解決しようとする課題】
【００１３】
それゆえ、劣化した信号のレベルを非線形かつロバストに判定するための方法及び装置を提供することが望ましい。
【００１４】
特に、この方法及び装置は、従来の技法よりも効率的にスパイク性の海面反射がある中で小さな物体を検出できるようにすることが望ましいが、この目的に限定するものではない。
【課題を解決するための手段】
【００１５】
本発明によれば、信号値を処理してそれら信号値のレベルを判定する方法が提供される。この方法は、シグナルプロセッサが、
複数の信号値のそれぞれを、円の少なくとも一部の円弧上の各点にマッピングするプロセスと、
それらの点を定義する座標のうちの少なくとも１つの座標の値を用いて平均値を求めるプロセスと、
信号値が閾値を超えるか否かを判定するプロセスと
を実行することを含む。
ここで、上記の判定は、
平均値を閾値と比較すること、又は
平均値を信号値の領域内の値に逆マッピングし、その逆マッピングされた値を閾値と比較することによって実施される。
【００１６】
本願の発明者は、このようにして平均値を求めると、その平均値は、他のタイプの平均（たとえば、算術平均）よりも、１つ又は複数の大きな信号値（インパルス性の雑音スパイクから生じる大きな値等）によって及ぼされる影響がはるかに小さいために、信号レベルをよりロバストに推定するものになることを見出した。
【００１７】
その結果として、インパルス性の雑音を有する、雑音のある背景における信号の検出を改善することができる。
【００１８】
本発明はさらに、上記の方法を実行するための装置も提供する。
【図面の簡単な説明】
【００１９】
【図１】各レンジセルからのレーダーエコー（radar return）を２値信号のストリームに変換するための比較器の従来技術における使用法を示す図である。
【図２】クラスタの存在を示すために２値データにおいて実行されるレンジエクステントフィルタの従来技術における動作を示す図である。
【図３】持続積分器の従来技術における動作を示す図である。
【図４】本発明の一実施形態による、円弧平均を求めるステップを示す図である。
【図５】本発明の一実施形態に従って構築される円弧平均計算器ＣＭＣの機能ブロック図である。
【図６】本発明の第２の実施形態に従って構築される円弧平均計算器ＣＭＣの機能ブロック図である。
【図７】ａ）は、本発明の第３の実施形態に従って構築される２つのマッピング関数を示す図である。ｂ）は、信号値の角度位置への変換を示す図である。ｃ）は、信号値の単位半円上へのマッピングを示す図である。
【図８】本発明の第３の実施形態に従って構築される修正型円弧平均計算器ＭＣＭＣの機能ブロック図である。
【図９】信号検出のために円弧平均ＣＭ及び円弧密度の両方が用いられるときの判定領域を示す図である。
【図１０】本発明の第４の実施形態に従って構築される円弧平均計算器ＣＭＣを組み込んだ信号検出システムＳＤＴの機能ブロック図である。
【発明を実施するための形態】
【００２０】
インパルス性雑音がある中での信号検出を分析したところ、図１ａにおいて示されるような信号値の指標として従来の算術平均を利用して検出閾値と比較しても、最適な検出性能を与えることができないことがわかった。それゆえ、信号検出にさらに適した信号レベルを求めるための改善されたロバストな方法が必要とされる。
【００２１】
本発明の実施形態によれば、信号の１組の観測値（たとえば、包絡線、大きさ、或いは電力）の各値ｐが、円弧上の位置を定義する直交軸上の２つの座標（ｘ，ｙ）を有する２次元空間内の１つの点にそれぞれ変換される。この円弧は１８０度（他の単位においても等量）以下の角度を有する。この円弧は半円弧であることが好ましい。個々の平均は、変換された点のｘ座標及びｙ座標について個別に計算される。次に、結果として生成される２つの平均を用いて、２次元空間内の変換された点の平均方向θMDを計算する。その後、この平均方向θMDを用いて、観測された信号値を閾値と比較する。第１の技法では、この比較は、逆変換Ｍ-1を適用して平均方向θMDを観測値の信号領域にマッピングし直して、本明細書において円弧平均ｐCMと呼ばれる量を与え、この円弧平均を信号領域において閾値と比較することによって実行される。その代わりとなるが同等の第２の技法では、この比較は、平均方向θMDを半円の２次元空間において閾値と比較することによって実行される。この閾値は、最初からその２次元空間内で定義してもよいし、或いは、観測値の信号領域において定義して、信号値を変換するために用いられるのと同じ変換を用いて半円の２次元空間内に変換してもよい。
【００２２】
代替の実施形態では、信号の１組の観測値ｐの各値が、２次元空間内の円弧上の点にそれぞれ変換される。２つの座標ｘ、ｙのうちの一方（両方ではない）の値から平均が計算される。その後、上記の２つの技法のうちの１つを用いて、観測される信号値を閾値と比較する際に、この平均が用いられる。
【００２３】
第１の実施形態
本発明の第１の実施形態によれば、以下の３つのステップ（図４においても概略的に示されている）を含む手順を適用することによって、円弧平均ｐCMが計算される。
【００２４】
ステップ１
信号の観測値ｐを、マッピング関数Ｍ（ｐ）を用いて単位半円上にマッピングする。その結果、目的とする所定の範囲（ＰＬ，ＰＨ）からの各値ｐkは、角度位置θkで単位半円上に配置される対応する点によって表される。
【００２５】
したがって、信号のＫ個の観測値
【００２６】
【数１】
【００２７】
は、以下の対応する１組のＫ個の角度によって表される。
【００２８】
【数２】
【００２９】
ただし、αは任意の初期角度である。
【００３０】
ステップ２
半円の２次元空間内にあるマッピングされた各点ｐkのｘ座標及びｙ座標を求めるために、角度｛θk｝の正弦及び余弦を計算する。その後、それらの計算値を別々に平均して２つの個々の平均値を求める。
【００３１】
【数３】
【００３２】
その後、以下の式から、１組の角度における平均方向θMDを求める。
【００３３】
【数４】
【００３４】
ステップ３
平均方向θMDに逆マッピングＭ-1（θ）を適用し、観測された信号値ｐの信号領域内に変換し直すことによって、以下に示す１組のＫ個の元になる信号値の円弧平均ｐCMを得る。
【００３５】
【数５】
【００３６】
続いて、円弧平均ｐCMを信号領域において閾値と比較する。（代替的には、方向θMDを半円の２次元空間において閾値と比較する）。
【００３７】
図４は、所与の１組の信号値の円弧平均ｐCMを求める手順を概略的に示している。この事例では、図示される半円は最初の２つの象限に現れるが（それゆえ、α＝０）、この特定の特徴は不可欠ではない。
【００３８】
第１の実施形態とその利点を理解するのを容易にするために、説明に役に立つ簡単な例を以下に与える。
【００３９】
第１のステップは、ｐの測定値を、（０，１８０）度の区間に入る値ｐ*に変換することである。ｐの測定値を以下のように仮定する。
【００４０】
【数６】
【００４１】
線形シフト演算によって各値から２００が減算され、以下の値が得られる。
【００４２】
【数７】
【００４３】
次に、線形スケーリング演算によって各値が２で除算され、以下の値が得られる。
【００４４】
【数８】
【００４５】
このとき、これらの変換値ｐ*は、（０，１８０）の区間に入っている。変換値ｐ*とそれらに対応する角度位置θ（度単位で表されている）とは、同じ数値を有する。
【００４６】
データをシフト及びスケーリングするためのこの方法は、この演算のために用いることができる数多くの方法の一例にすぎないことに留意されたい。
【００４７】
上記の１組の変換値ｐ*の算術平均ｐAMは、以下の通りである。
【００４８】
【数９】
【００４９】
上記の式から明らかなように、単一の相対的に大きな観測値ｐ6*＝１５０（たとえば、干渉する雑音スパイク）が、残りの５個の観測値ｐ1*〜ｐ5*の平均（３０）を大きくシフトする。
【００５０】
この実施形態における処理の第２のステップでは、変換値ｐ1*〜ｐ6*の平均方向θMDが、以下の式から求められる。
【００５１】
【数１０】
【００５２】
それゆえ、円弧領域内における算術平均ｐAMは５０に等しかったが、観測値の平均方向θMDは４１に等しい。平均方向θMDは、単一のより大きな観測値（すなわち、ｐ6*）によって及ぼされる影響がはるかに小さいので、値ｐの平均レベルのよりロバストな推定量である。
【００５３】
この処理の第３のステップにおいて、本実施形態では、平均方向θMD＝４１度を、逆マッピングＭ-1を用いて変換して測定信号ｐと同じ領域に戻す。
【００５４】
それゆえ、変換された値ｐCM＝（４１×２）＋２００＝２８２である。
【００５５】
一方、信号領域内の信号値ｐの算術平均は、［（５×２６０）＋５００］／６＝３００である。
【００５６】
したがって、本実施形態によって生成される円弧平均ｐCMは、従来の算術平均よりも信号領域内の値ｐの平均レベルをよりロバストに示している。これは、円弧平均ｐCMの値のほうが、算術平均よりも誤った信号値ｐ6によって及ぼされる影響が小さいためである。
【００５７】
ｐCMの上記の計算を実行するための装置が図５に示されている。
【００５８】
図５に示される計算器は、以下の機能ユニットを備える。
− シフト及びスケールユニット５０２；
− ２つの非線形変換器ＳＮＬ５０３及びＣＮＬ５０４；
− ２つのタップ付き遅延線ＤＬＳ５０５及びＤＬＣ５０６；
− ２つの平均化回路ＡＶＳ５０７及びＡＶＣ５０８；
− 演算ユニットＡＴＵ５０９。
【００５９】
包絡線、大きさ、或いは電力のような正の入力信号ｐｐがシフト及びスケールユニット５０２に通されると、ユニット５０２は、入力信号の各値を１８０度にわたる範囲内の値のうちの１つにマッピングする。それにより、正規化された信号ＮＰが生成される。
【００６０】
正規化された信号レベルＮＰは、２つの非線形変換器ＳＮＬ５０３及びＣＮＬ５０４に並列に加えられる。変換器の出力ＳＳ及びＣＣは、共通の入力値NPから２つの適切に選択されたマッピング関数を利用することによって求められる。検討している事例では、以下の通りである。
【００６１】
【数１１】
【００６２】
先に述べたように、値ＳＳ及びＣＣは、２次元空間内の値ＮＰのｘ座標及びｙ座標を定義する。
【００６３】
信号ＳＳ及びＣＣは、それぞれのタップ付き遅延線ＤＬＳ５０５及びＤＬＣ５０６に沿って伝搬する：信号ＳＳのＫ個のサンプルを遅延セルＳ１、Ｓ２、．．．、ＳＫの出力において取得することができ、一方、信号ＣＣのＫ個のサンプルを遅延セルＣ１、Ｃ２、．．．、ＣＫの出力において取得することができる。
【００６４】
平均化回路ＡＶＳ５０７は、その出力において、セルＳ１、Ｓ２、．．．、ＳＫから得られる入力値の和に比例する値ＡＳを生成する。同様に、平均化回路ＡＶＣ５０８は、その出力において、セルＣ１、Ｃ２、．．．、ＣＫから得られる入力値の和に比例する値ＡＣを生成する。これら２つの値ＡＳ及びＡＣは、演算ユニットＡＴＵ５０９によって用いられ、円弧平均ＣＭが求められる。
【００６５】
その後、比較ユニット（図５には示されていない）において、計算された円弧平均ＣＭが閾値と比較される。
【００６６】
第２の実施形態
上記の第１の実施形態では、信号値ｐは「シフト及びスケール」タイプの線形演算を利用することによって角度値θに変換された。しかしながら、実用的な用途では、ダイナミックレンジを非線形に調整するために、観測された信号値に対して最初に非線形（たとえば対数）変換を適用した後、そのように変換された値を単位半円にマッピングすることが好都合である場合がある。そのような処理を実行する実施形態が以下に示される。
【００６７】
たとえば、有用な非線形マッピングは以下の形のものである。
【００６８】
【数１２】
【００６９】
ただし、クリップ関数Ｈ（・）は、その引数の最小値と最大値をそれぞれ−π／２とπ／２に制限し、γは、処理される信号のダイナミックレンジをさらに調整するために用いられるスケーリング係数である。
【００７０】
たとえば、観測値ｐの範囲が０．０１から１００までの場合には、γ＝π／４がｐの全ての値を（−π／２，π／２）内に収めることになる。ただし、両端も範囲に含む。０．０１〜１００の範囲外であるｐの新たな値が後に検出された場合に、その新たなｐ値に対してもγの同じ値が用いられるようにするために、クリップ関数Ｈには、ｐの新たなマッピング値を（−π／２，π／２）の範囲に制限することが要求される。実際には、ｐ値が大抵の場合に０．０１〜１００の範囲に入っており、これらの限界値に達するか又は限界値を超えることがめったにない場合、クリップ関数Ｈは通常は用いられない。
【００７１】
以下の理由から、半円上にマッピングする前に対数変換を適用することが、他の非線形変換を用いることよりも好ましい場合がある。
１．対数変換は圧縮性があり、大きい観測値ほど抑制する。
２．除算が単純な減算に変わるので、元になるデータのスケールパラメータがシフトパラメータに変換され、それにより、信号正規化の演算が簡略化される。
【００７２】
図６は、本発明の第２の実施形態に従って構築される円弧平均計算器ＣＭＣの機能ブロック図である。
【００７３】
この計算器は、以下の機能ユニットを備える。
− 対数変換器ＬＧＩ６０１；
− 減算器ＳＢＴ６０２；
− ２つの非線形変換器ＳＮＬ５０３及びＣＮＬ５０４；
− ２つのタップ付き遅延線ＤＬＳ５０５及びＤＬＣ５０６；
− ２つの平均化回路ＡＶＳ５０７及びＡＶＣ５０８；
− 演算ユニットＡＴＵ６０９。
【００７４】
したがって、図５に示される第１の実施形態の計算器と比べると、シフト及びスケールユニット５０２の代わりに、対数変換器ＬＧＩ６０１及び減算器ＳＢＴ６０２が用いられている。他の全ての構成要素は同じであるが、ＡＴＵ６０９は例外であり、ｐに対して実行される対数演算を考慮するために、第１の実施形態と比べて異なる逆マッピングＭ-1を実行する。
【００７５】
包絡線、大きさ、或いは電力のような正の入力信号ＰＰが対数変換器ＬＧＩ６０１に通されると、入力信号ＰＰのレベルの対数指標である信号LＰが生成される。すなわち、
【００７６】
【数１３】
【００７７】
次に、減算器ＳＢＴ６０２において、信号ＬＰから目的とする或る基準レベルＢＧの対数ＢＬが減算されることによって、信号ＬＰが正規化される。
【００７８】
【数１４】
【００７９】
たとえば、信号検出問題では、基準レベルは長期の観測値から得られる背景雑音の平均レベルとすることができる。ＬＰから値ＢＬを減算する操作は、ＬＰの値を（−π／２，π／２）の角度範囲を有する単位半円上の点にマッピングする。これは、係数γによって与えられるスケーリングと等価である。
【００８０】
正規化された信号レベルＮＰは、２つの非線形変換器ＳＮＬ５０３及びＣＮＬ５０４に並列に加えられる。変換器の出力ＳＳ及びＣＣは、共通の入力値NPから２つの適切に選択されたマッピング関数を利用することによって求められる。本実施形態では、以下の通りである。
【００８１】
【数１５】
【００８２】
信号ＳＳ及びＣＣは、それぞれのタップ付き遅延線ＤＬＳ５０５及びＤＬＣ５０６に沿って伝搬する：信号ＳＳのＫ個のサンプルを遅延セルＳ１、Ｓ２、．．．、ＳＫの出力において取得することができ、一方、信号ＣＣのＫ個のサンプルを遅延セルＣ１、Ｃ２、．．．、ＣＫの出力において取得することができる。
【００８３】
平均化回路ＡＶＳ５０７は、その出力において、セルＳ１、Ｓ２、．．．、ＳＫから得られる入力値の和に比例する値ＡＳを生成する。同様に、平均化回路ＡＶＣ５０８は、その出力において、セルＣ１、Ｃ２、．．．、ＣＫから得られる入力値の和に比例する値ＡＣを生成する。これら２つの値ＡＳ及びＡＣは、演算ユニットＡＴＵ６０９によって用いられ、円弧平均ＣＭが求められる。
【００８４】
その後、比較ユニット（図６には示されていない）において、計算された円弧平均が閾値と比較される。
【００８５】
第３の実施形態
上記の実施形態では、信号値ｐを単位半円上の点にマッピングした後、三角演算子を適用して点の位置を定義する半円の２次元空間内の２つの座標を求める。続いて、これらの座標を用いて平均方向θMDを計算し、必要に応じて円弧平均ｐCMを計算する。しかしながら、信号値ｐの半円への初期マッピングを、結果として生成される半円上の点の２つの座標を直接与えるように実行してもよい。したがって、第１及び第２の実施形態における三角演算を実行することによって座標を計算する必要はない。
【００８６】
そのような処理を実行する一実施形態が以下に説明される。
【００８７】
一般的に、信号値ｐの半円へのマッピングは、適切に構築された２つのマッピング関数Ｓ（ｐ）及びＣ（ｐ）を使用することによって実行することができる。マッピングは半円を生成することを要求されるので、それらのマッピング関数は以下の条件を満たさなければならない。
【００８８】
【数１６】
【００８９】
ただし、ｂは定数である。単位半円の場合、ｂ＝１である。
【００９０】
本発明の第３の実施形態によれば、第１のマッピング関数Ｓ（ｐ）は単調に増加する連続関数になる。この関数は、最小信号レベルＰＬの場合に最小値である−１を有し、最大信号レベルＰＨにおいて＋１に等しい最大値に達すると仮定する。そのような関数は無限に存在し、正弦波の適切に選択された一部分は数ある選択肢のうちの１つであるが、別の有用な関数が以下において検討される。
【００９１】
次に、第２のマッピングＣ（ｐ）が、マッピング関数Ｓ（ｐ）から以下のようにして得られる。
【００９２】
【数１７】
【００９３】
それゆえ、２つのマッピング関数Ｓ（ｐ）及びＣ（ｐ）を利用するマッピング演算Ｍ［Ｓ（ｐ），Ｃ（ｐ）］は、１次元のｐ空間から２次元の（Ｓ，Ｃ）空間への非線形変換と見なすことができる。結果として、各信号値ｐは、その点の一方の座標がＳ（ｐ）によって定義され、その点の他方の座標がＣ（ｐ）によって定義されるように、単位半円上の点に直接マッピングされる。その後、Ｓ（ｐ）及びＣ（ｐ）に算術平均を直接適用して、平均方向のためのパラメータを求めることができる。
【００９４】
より詳細には、平均方向θMDは、以下のように計算される。
【００９５】
【数１８】
【００９６】
その後、平均方向θMDを閾値と比較することができる。代替的には、必要であれば平均方向θMDに逆マッピングＭ-1（θ）を適用することによって、信号値ｐの円弧平均ｐCMが求められ、この円弧平均ｐCMが閾値と比較される。逆マッピングＭ-1は、Ｍ-1を定義する数学的な関数を計算して評価することによって、又は数値技法（反復法等）を用いて要求される値を求めることによって、求めることができる。
【００９７】
双曲線関数ｔａｎｈ ｗ及び１／（ｃｏｓｈ ｗ）は、マッピング関数Ｓ（ｐ）及びＣ（ｐ）として有利に利用することができる。
【００９８】
より詳細には、
【００９９】
【数１９】
【０１００】
であるので、マッピング［ｔａｎｈ ｗ，１／（ｃｏｓｈ ｗ）］は、値ｗを表す点を単位半円上に置く。
【０１０１】
正規化するために、双曲線変換ｔａｎｈ ｗ及び１／（ｃｏｓｈ ｗ）を実行する前に、対数変換を使用してｗ＝ｌｎ（ｐ）とすることが好都合であろう。この場合、円弧平均計算器の機能ブロック図は、図６に示される機能ブロック図と類似の構造を有する。しかしながら、この場合、非線形変換器ＳＮＬ５０３及びＣＮＬ５０４は、以下のマッピング演算を実行する。
【０１０２】
【数２０】
【０１０３】
対数変換と双曲線変換とを組み合わせることによって、円弧平均計算器の構造を単純化することもできる。そのような単純化は以下のように行うことができる。
【０１０４】
【数２１】
【０１０５】
であるので、双曲線変換ｔａｎｈ（ｗ）に先立って対数変換ｗ＝ｌｎ（ｐ）が行なわれる場合には、以下の式が成り立つ。
【０１０６】
【数２２】
【０１０７】
図７ａ及び図７ｂに示されるように、２つのマッピング関数の形状は、信号値ｐと、単位半円上の値を表す点の対応する角度位置θとの間の関係を表す。この関係は非線形であり、「ソフト」リミッタのその形は、以下の数学的な表現から生じる。
【０１０８】
【数２３】
【０１０９】
図７ｃは、信号のいくつかの選択された元になる値ｐに対応する角度位置θを有する単位半円を示す。図に示されるように、信号の値が大きいほど（ｐ＞１０）、対応する角度位置θが限界値π／２に近いクラスタを形成する。その一方で、１よりも小さいｐの値は、象限全体（−π／２，０）を占有する角度位置θを生成する。
【０１１０】
観測値｛ｐ1，ｐ2，．．．，ｐK-1，ｐK｝からパラメータＹK及びＸKが、
【０１１１】
【数２４】
【０１１２】
に従って求められると、信号値ｐの円弧平均ｐCMは、以下の式から得られる。
【０１１３】
【数２５】
【０１１４】
図８は、本発明の第３の実施形態に従って構築される修正型円弧平均計算器ＭＣＭＣ８００の機能ブロック図である。この構成では、図６の円弧平均計算器６００において用いられる対数変換器ＬＧＩ６０１及び減算器ＳＢＴ６０２が、可変利得増幅器ＶＧＩ８０１によって置き換えられている。この増幅器は、目的とする基準レベルを示す信号ＢＧに応じてその利得を調整する。
【０１１５】
この構成では、変換器ＳＮＬ８０２及びＣＮＬ８０３の出力ＳＳ及びＣＣは、共通入力ＣＰから、以下のようにして得られる。
【０１１６】
【数２６】
【０１１７】
修正型システムによって実行される他の機能及び演算は、図６のシステムによって実行されるのと同様であり、したがって、ここでは再び説明しない。
【０１１８】
本発明のこの実施形態に従って構築される円弧平均計算器が対象とする応用例のうちの１つは、背景雑音がある中での信号の検出である。以下の例は、この実施形態のそのような応用例を示す。この例では、対数関数及び双曲線関数がマッピングのために利用される。
【０１１９】
単位電力熱雑音と、時折、雑音レベルの１０倍を超えるスパイクが生じるインパルス性干渉とを含む雑音のある背景内において、単位電力のランダム信号が検出されると仮定する。また、例示のために、検出閾値が１．９に設定されていると仮定する。
【０１２０】
信号が無い場合に、背景雑音の６個の観測値が以下のようであると仮定する。
【０１２１】
【数２７】
【０１２２】
５個の等しいサンプルは熱雑音レベルを表すものであり、サンプル番号６はインパルス性干渉によって生成されたものであろう。この場合、算術平均ＰAMは２．５に等しい。それゆえ、値２．５は１．９の検出閾値よりも大きいので、算術平均ＰAMが検出統計値として用いられる場合、誤検出（false alarm）が宣言されるであろう。
【０１２３】
ここで、算術平均ではなく円弧平均ｐCMを用いて背景雑音のレベルを判定すると仮定する。パラメータＹ6及びＸ6は以下のようになる。
【０１２４】
【数２８】
【０１２５】
それゆえ、円弧平均ｐCMは以下のようになる。
【０１２６】
【数２９】
【０１２７】
上記より、この場合、円弧平均ｐCMは検出閾値を超えないので、誤検出は生じない。
【０１２８】
ここで、検出されるべきランダム信号が存在し、この信号＋雑音の事例においては６個の観測値が全て等しく、たとえばｐ1＝ｐ2＝．．．＝ｐ6＝２であると仮定する。この場合、算術平均ｐAMと円弧平均ｐCMはともに等しく、すなわちｐAM＝ｐCM＝２であり、いずれも正しい検出判定に繋がることは明らかである。
【０１２９】
したがって、円弧平均を検出統計値として利用する主な利点は、時折生じるより大きな観測値を減衰させる能力から得られる。一般的に、そのような特性はインパルス性の雑音によって劣化した信号を処理する際に非常に有用である。
【０１３０】
第４の実施形態
本発明のさらなる実施形態によれば、単位半円上の点の円弧密度ＤCCが計算され、それを用いて円弧平均ｐCMと比較される検出閾値が設定される。かつ／又は、円弧密度は、検出閾値と比較される以前の円弧平均の値の調整に用いられる。
【０１３１】
信号レベルｐの値｛ｐk｝を表す点｛θk｝の円弧密度ＤCCは、以下の式から求められる。
【０１３２】
【数３０】
【０１３３】
ただし、
【０１３４】
【数３１】
【０１３５】
であり、Ｓ（ｐk）及びＣ（ｐk）は、マッピングのために用いられる２つの関数である。
【０１３６】
値ＸK及びＹKは、第３の実施形態において信号値ｐの半円上への初期マッピングによって出力されるので、第３の実施形態はこれらの値ＸK及びＹKを計算するのに最も適している。しかしながら、第１及び第２の実施形態の場合のように、ｐを半円上にマッピングした後に、ＸK及びＹKを計算してもよいことは、当業者には明らかであろう。
【０１３７】
点｛θk｝が半円上の密集したクラスタを形成するとき、円弧密度ＤCCの値は１に近くなる。しかしながら、点｛θk｝が、たとえば支配的な雑音に起因して広範囲に分散するとき、円弧密度ＤCCの値は著しく小さくなるであろう。したがって、円弧密度ＤCCの値を円弧平均と共に用いて、信号検出をさらに改善することができる。
【０１３８】
図９は、円弧平均ＣＭを円弧密度ＤＣと適切に組み合わせて検出判定領域を構築するときの検出性能の潜在的な改善を示す。
【０１３９】
検出のために円弧平均ＣＭだけが利用されるとき、円弧密度にかかわらず、円弧平均の観測値が判定閾値ＴＨを超える場合に、信号の存在が宣言される。この場合、判定領域Ｄ１は、ラインＴＨから延在する長方形の形状を有する。
【０１４０】
しかしながら、円弧平均ＣＭが円弧密度ＤＣと共に用いられるとき、結果として生成される判定領域は、領域Ｄ２だけ拡大される。円弧密度値が或る所定の閾値Ｃ１を超えるとき、円弧平均のための判定閾値ＴＨは新たな下側値Ｔ１まで徐々に下げられるであろう。明らかに、面積（Ｄ１＋Ｄ２）は単独の面積Ｄ１よりも大きく、閾値ＴＨが判定境界ＤＢによって置き換えられると、改善された検出性能が達成される。
【０１４１】
図１０は、円弧密度を用いて検出閾値を設定する、本発明の第４の実施形態に従って構築される信号検出器ＳＤＴ１０００の機能ブロック図である。信号検出器ＳＤＴ１０００の目的とする用途は、海面反射がある中での小さな物体の検出である。背景雑音、又は信号＋雑音の観測サンプルＰＰが、円弧平均計算器ＣＭＣ１００１において処理される。背景雑音の平均レベルを示す補助信号ＢＬを利用することによって、サンプルレベルが正規化される。そのような信号は、たとえば、被試験セルに隣接する複数のレンジセルによって生成される観測値を長い時間間隔にわたって取り込んだものを平均することによって得ることができる。
【０１４２】
円弧平均計算器ＣＭＣ１００１は、それぞれ円弧平均及び円弧密度を示している２つの出力ＣＭ及びＤＣを提供する。これらの信号は、判定ブロックＤＥＴ１０１１によって用いられ、観測サンプルＰＰが海面反射だけによって生成されたのか、それとも海面反射内に埋もれた物体によって生成されたのかが判定される。そのために、判定境界を定義する信号ＤＢが、判定ブロックＤＥＴ１０１１の入力ＤＢに加えられる。このブロックの出力信号ＳＤは、海面反射内の信号の存否に関する包括的な判定である。
【０１４３】
上記の実施形態では、円弧密度ＤCCは、円弧平均ｐCMと比較される閾値を調整するために用いられることに留意されたい（図９に示されている）。しかしながら、これは、円弧密度ＤCCに基づいて円弧平均ｐCMを調整し、調整された平均を無変更の閾値と比較することと等価である。したがって、この実施形態では、閾値の代わりにｐCMを調整してもよい。同様に、円弧密度ＤCCに基づいて閾値及び円弧平均ｐCMの両方を調整してもよい。
【０１４４】
変更及び変形
本発明の好ましい実施形態のこれまでの説明は、例示及び説明のために提示されてきたものであり、本発明を網羅的に述べることも、開示されるのと全く同じ形に限定することも意図していない。これまでの説明に鑑みて、数多くの改変、変更、及び変形によって、当業者は本発明を検討される特定の用途に合わせて種々の実施形態において利用できるようになるのは明らかである。
【０１４５】
たとえば、上記の実施形態では、信号値ｐが半円上の点にマッピングされる。しかしながら、半円全体を用いることは不可欠ではなく、代わりに、半円よりも小さな円弧を用いることもできる。より詳細には、信号値ｐは、半円弧によって定められる角度、すなわちπラジアン（又は同じく、たとえば１８０度）以下の角度を有する任意の円弧にマッピングすることができる。
【０１４６】
上記の全ての実施形態において、ＸK（すなわち、半円上の点のｘ座標の平均）及びＹK（すなわち、半円上の点のｙ座標の平均）の両方を求めた後、それら両方の値を用いてθMDを計算することによって、平均方向θMDが計算される。しかしながら、ロバストではないが、ＸK及びＹKの両方を求める必要はなく、平均方向θMDを計算する代わりに、代替の平均ＭALTを以下のように計算することもできる。
【０１４７】
【数３２】
【０１４８】
より詳細には、点ｐが角度範囲０度〜１８０度又は１８０度〜３６０度を含む単位半円にマッピングされるとき、ＸKの計算に用いられる余弦値がこの範囲内で一意の値を有するので、ＭALT＝ＸKが用いられる。すなわち、半円上の点ｐの値とｘ軸との間に１：１のマッピングが存在する。一方、点ｐが角度範囲９０度〜２７０度又は２７０度〜９０度を含む単位半円にマッピングされるとき、ＹKの計算に用いられる正弦値がこの範囲内で一意の値を有するので、ＭALT＝ＹKが用いられる。
【０１４９】
その後、代替の平均ＭALTを閾値と比較することもできるし、或いは必要に応じて信号領域にマッピングし直した後に、閾値と比較することもできる。たとえば、第１及び第２の実施形態では、このようにマッピングし直すことは、逆余弦マッピング（ＭALT＝ＸKの場合）又は逆正弦マッピング（ＭALT＝ＹKの場合）を用いて半円上の点にマッピングし直した後、逆マッピングＭ-1を適用することによって実行することができる。
【０１５０】
上記の第４の実施形態では、円弧平均ｐCMが、円弧密度ＤCCに基づいて設定される閾値と比較される（又は、円弧平均ｐCMが円弧密度ＤCCに基づいて調整される）。しかしながら、代わりに、平均方向θMDが、円弧密度ＤCCに基づいて設定される閾値と比較されてもよい（又は、平均方向θMDが円弧密度ＤCCに基づいて調整されてもよい）。
【０１５１】
当然のことながら、他の変更も可能である。
【特許請求の範囲】
【請求項１】
信号処理方法であって、
複数の信号値のそれぞれを、直交軸上の第１の座標及び第２の座標によって定義される２次元空間内における、半円弧によって定められる角度以下の角度を有する円弧上の各点にマッピングすることと、
前記円弧上のマッピングされた各点の第１の座標及び第２の座標のうちの少なくとも一方を用いて平均値を計算することと、
前記計算された平均値を用いて前記信号値が閾値を超えるか否かを判定することと
を含む、信号処理方法。
【請求項２】
前記平均値は、
前記マッピングされた各点の第１の座標値の平均を含む第１の平均値を計算することと、
前記マッピングされた各点の第２の座標値の平均を含む第２の平均値を計算することと、
前記第１の平均値及び前記第２の平均値に基づいて、前記２次元空間内の平均方向を計算することと
によって計算される、請求項１に記載の方法。
【請求項３】
前記信号値のそれぞれは、
前記信号値を前記円弧上の点を定義する角度を表す値にマッピングすることと、
前記角度に三角関数を適用することによって前記点の第１の座標及び第２の座標の値を計算することと
によって、前記円弧上の各点にマッピングされる、請求項２に記載の方法。
【請求項４】
前記信号値のそれぞれは、
前記信号値に第１の関数Ｓ（ｐ）を適用して前記円弧上の前記点の第１の座標のための値を計算すると共に、前記信号値に第２の関数Ｃ（ｐ）を適用して前記円弧上の前記点の第２の座標のための値を計算することによって、前記円弧上の各点にマッピングされ、
ただし、前記第１の関数及び前記第２の関数は、Ｓ2（ｐ）＋Ｃ2（ｐ）＝ｂの関係を満たし、ｐは信号値であり、ｂは定数である、請求項２に記載の方法。
【請求項５】
前記第１の関数及び前記第２の関数のそれぞれは双曲線関数を含む、請求項４に記載の方法。
【請求項６】
前記信号値のそれぞれは、
第１の関数Ｓ（ｐ）＝（ｐ2−１）／（ｐ2＋１）を適用することと、
第２の関数Ｃ（ｐ）＝２ｐ／（ｐ2＋１）を適用することと
によって、前記円弧上の各点にマッピングされる、請求項４に記載の方法。
【請求項７】
前記信号値が閾値を超えるか否かを判定する前記ステップは、
前記計算された平均値を閾値と比較すること、及び
前記計算された平均値を前記信号領域内の値に逆マッピングし、該逆マッピングされた値を閾値と比較すること
のうちの１つを含む、請求項１〜６のいずれか一項に記載の方法。
【請求項８】
前記円弧上の前記マッピングされた点の密度の指標が求められ、
前記求められた密度の指標を用いて前記閾値が設定され、かつ／又は、前記求められた密度の指標を用いて前記閾値と比較される前記平均値が調整される、請求項７に記載の方法。
【請求項９】
前記円弧上の前記マッピングされた点の前記密度の指標は、前記第１の平均値及び前記第２の平均値に基づいて
【数１】
に従って求められ、
ただし、ＤCCは前記密度の指標であり、ＸKは前記第１の平均値であり、ＹKは前記第２の平均値である、請求項２〜６のいずれか一項に従属するときの請求項８に記載の方法。
【請求項１０】
前記円弧密度の指標は、前記閾値を設定するために用いられ、
前記密度の指標を密度閾値と比較することと、
前記密度の指標が前記密度閾値を超えるときに、予め設定された値を小さくすることによって、前記密度の指標と前記密度閾値との間の関係に基づいて前記閾値を設定することと
によって、前記閾値を設定する、請求項８又は９に記載の方法。
【請求項１１】
信号処理装置であって、
複数の信号値のそれぞれを、直交軸上の第１の座標及び第２の座標によって定義される２次元空間内における、半円弧によって定められる角度以下の角度を有する円弧上の各点にマッピングするように動作可能な信号値マップ器と、
前記円弧上のマッピングされた各点の第１の座標及び第２の座標のうちの少なくとも一方を用いて平均値を計算するように動作可能な平均値計算器と、
前記計算された平均値を用いて前記信号値が閾値を超えるか否かを判定するように構成される判定ユニットと
を備える、信号処理装置。
【請求項１２】
前記平均値計算器は、
前記マッピングされた各点の第１の座標値の平均を含む第１の平均値を計算するように動作可能な第１の座標平均値計算器と、
前記マッピングされた各点の第２の座標値の平均を含む第２の平均値を計算するように動作可能な第２の座標平均値計算器と、
前記第１の平均値及び前記第２の平均値に基づいて前記２次元空間内の平均方向を計算するように動作可能な平均方向計算器と
を含む、請求項１１に記載の装置。
【請求項１３】
前記信号値マップ器は、
前記信号値を前記円弧上の点を定義する角度を表す値にマッピングすることと、
前記角度に三角関数を適用することによって前記点の第１の座標及び第２の座標の値を計算することと
によって、前記信号値のそれぞれを前記円弧上の各点にマッピングするように構成される、請求項１２に記載の装置。
【請求項１４】
前記信号値マップ器は、
前記信号値に第１の関数Ｓ（ｐ）を適用して前記円弧上の前記点の第１の座標のための値を計算すると共に、前記信号値に第２の関数Ｃ（ｐ）を適用して前記円弧上の前記点の第２の座標のための値を計算することによって、前記信号値のそれぞれを前記円弧上の各点にマッピングするように構成され、
ただし、前記第１の関数及び前記第２の関数は、Ｓ2（ｐ）＋Ｃ2（ｐ）＝ｂの関係を満たし、ｐは信号値であり、ｂは定数である、請求項１２に記載の装置。
【請求項１５】
前記信号値マップ器は、第１の関数及び第２の関数を適用することによって、前記信号値のそれぞれを前記円弧上の各点にマッピングするように構成され、該第１の関数及び該第２の関数のそれぞれは双曲線関数を含む、請求項１４に記載の装置。
【請求項１６】
前記信号値マップ器は、
第１の関数Ｓ（ｐ）＝（ｐ2−１）／（ｐ2＋１）を適用することと、
第２の関数Ｃ（ｐ）＝２ｐ／（ｐ2＋１）を適用することと
によって、前記信号値のそれぞれを前記円弧上の各点にマッピングするように構成される、請求項１４に記載の装置。
【請求項１７】
前記判定ユニットは、
前記計算された平均値を閾値と比較するように動作可能な比較器を有する検出ユニット、及び
前記計算された平均値を前記信号領域内の値に逆マッピングするように動作可能な逆マップ器と、該逆マッピングされた値を閾値と比較するように動作可能な比較器とを有する検出ユニット、
のうちの一方を含む、請求項１１〜１６のいずれか一項に記載の装置。
【請求項１８】
前記円弧上の前記マッピングされた点の密度の指標を計算するように動作可能な密度指標計算器をさらに含み、
前記判定ユニットは、前記求められた密度の指標を用いて前記閾値を設定すること、及び、前記求められた密度の指標を用いて前記閾値と比較される前記平均値を調整することのうちの少なくとも一方を含む処理を実行するように構成される、請求項１７に記載の装置。
【請求項１９】
前記密度指標計算器は、前記円弧上の前記マッピングされた点の前記密度の指標を、前記第１の平均値及び前記第２の平均値に基づいて
【数２】
に従って計算するように構成され、
ただし、ＤCCは前記密度の指標であり、ＸKは前記第１の平均値であり、ＹKは前記第２の平均値である、請求項１２〜１６のいずれか一項に従属するときの請求項１８に記載の装置。
【請求項２０】
前記密度指標計算器は、前記円弧密度の指標を用いて前記閾値を設定するように構成され、
前記密度の指標を密度閾値と比較することと、
前記密度の指標が前記密度閾値を超えるときに、予め設定された値を小さくすることによって、前記密度の指標と前記密度閾値との間の関係に基づいて前記閾値を設定することと
によって、前記閾値を設定する、請求項１８又は１９に記載の装置。
【請求項１】
信号処理方法であって、
複数の信号値のそれぞれを、直交軸上の第１の座標及び第２の座標によって定義される２次元空間内における、半円弧によって定められる角度以下の角度を有する円弧上の各点にマッピングすることと、
前記円弧上のマッピングされた各点の第１の座標及び第２の座標のうちの少なくとも一方を用いて平均値を計算することと、
前記計算された平均値を用いて前記信号値が閾値を超えるか否かを判定することと
を含む、信号処理方法。
【請求項２】
前記平均値は、
前記マッピングされた各点の第１の座標値の平均を含む第１の平均値を計算することと、
前記マッピングされた各点の第２の座標値の平均を含む第２の平均値を計算することと、
前記第１の平均値及び前記第２の平均値に基づいて、前記２次元空間内の平均方向を計算することと
によって計算される、請求項１に記載の方法。
【請求項３】
前記信号値のそれぞれは、
前記信号値を前記円弧上の点を定義する角度を表す値にマッピングすることと、
前記角度に三角関数を適用することによって前記点の第１の座標及び第２の座標の値を計算することと
によって、前記円弧上の各点にマッピングされる、請求項２に記載の方法。
【請求項４】
前記信号値のそれぞれは、
前記信号値に第１の関数Ｓ（ｐ）を適用して前記円弧上の前記点の第１の座標のための値を計算すると共に、前記信号値に第２の関数Ｃ（ｐ）を適用して前記円弧上の前記点の第２の座標のための値を計算することによって、前記円弧上の各点にマッピングされ、
ただし、前記第１の関数及び前記第２の関数は、Ｓ2（ｐ）＋Ｃ2（ｐ）＝ｂの関係を満たし、ｐは信号値であり、ｂは定数である、請求項２に記載の方法。
【請求項５】
前記第１の関数及び前記第２の関数のそれぞれは双曲線関数を含む、請求項４に記載の方法。
【請求項６】
前記信号値のそれぞれは、
第１の関数Ｓ（ｐ）＝（ｐ2−１）／（ｐ2＋１）を適用することと、
第２の関数Ｃ（ｐ）＝２ｐ／（ｐ2＋１）を適用することと
によって、前記円弧上の各点にマッピングされる、請求項４に記載の方法。
【請求項７】
前記信号値が閾値を超えるか否かを判定する前記ステップは、
前記計算された平均値を閾値と比較すること、及び
前記計算された平均値を前記信号領域内の値に逆マッピングし、該逆マッピングされた値を閾値と比較すること
のうちの１つを含む、請求項１〜６のいずれか一項に記載の方法。
【請求項８】
前記円弧上の前記マッピングされた点の密度の指標が求められ、
前記求められた密度の指標を用いて前記閾値が設定され、かつ／又は、前記求められた密度の指標を用いて前記閾値と比較される前記平均値が調整される、請求項７に記載の方法。
【請求項９】
前記円弧上の前記マッピングされた点の前記密度の指標は、前記第１の平均値及び前記第２の平均値に基づいて
【数１】
に従って求められ、
ただし、ＤCCは前記密度の指標であり、ＸKは前記第１の平均値であり、ＹKは前記第２の平均値である、請求項２〜６のいずれか一項に従属するときの請求項８に記載の方法。
【請求項１０】
前記円弧密度の指標は、前記閾値を設定するために用いられ、
前記密度の指標を密度閾値と比較することと、
前記密度の指標が前記密度閾値を超えるときに、予め設定された値を小さくすることによって、前記密度の指標と前記密度閾値との間の関係に基づいて前記閾値を設定することと
によって、前記閾値を設定する、請求項８又は９に記載の方法。
【請求項１１】
信号処理装置であって、
複数の信号値のそれぞれを、直交軸上の第１の座標及び第２の座標によって定義される２次元空間内における、半円弧によって定められる角度以下の角度を有する円弧上の各点にマッピングするように動作可能な信号値マップ器と、
前記円弧上のマッピングされた各点の第１の座標及び第２の座標のうちの少なくとも一方を用いて平均値を計算するように動作可能な平均値計算器と、
前記計算された平均値を用いて前記信号値が閾値を超えるか否かを判定するように構成される判定ユニットと
を備える、信号処理装置。
【請求項１２】
前記平均値計算器は、
前記マッピングされた各点の第１の座標値の平均を含む第１の平均値を計算するように動作可能な第１の座標平均値計算器と、
前記マッピングされた各点の第２の座標値の平均を含む第２の平均値を計算するように動作可能な第２の座標平均値計算器と、
前記第１の平均値及び前記第２の平均値に基づいて前記２次元空間内の平均方向を計算するように動作可能な平均方向計算器と
を含む、請求項１１に記載の装置。
【請求項１３】
前記信号値マップ器は、
前記信号値を前記円弧上の点を定義する角度を表す値にマッピングすることと、
前記角度に三角関数を適用することによって前記点の第１の座標及び第２の座標の値を計算することと
によって、前記信号値のそれぞれを前記円弧上の各点にマッピングするように構成される、請求項１２に記載の装置。
【請求項１４】
前記信号値マップ器は、
前記信号値に第１の関数Ｓ（ｐ）を適用して前記円弧上の前記点の第１の座標のための値を計算すると共に、前記信号値に第２の関数Ｃ（ｐ）を適用して前記円弧上の前記点の第２の座標のための値を計算することによって、前記信号値のそれぞれを前記円弧上の各点にマッピングするように構成され、
ただし、前記第１の関数及び前記第２の関数は、Ｓ2（ｐ）＋Ｃ2（ｐ）＝ｂの関係を満たし、ｐは信号値であり、ｂは定数である、請求項１２に記載の装置。
【請求項１５】
前記信号値マップ器は、第１の関数及び第２の関数を適用することによって、前記信号値のそれぞれを前記円弧上の各点にマッピングするように構成され、該第１の関数及び該第２の関数のそれぞれは双曲線関数を含む、請求項１４に記載の装置。
【請求項１６】
前記信号値マップ器は、
第１の関数Ｓ（ｐ）＝（ｐ2−１）／（ｐ2＋１）を適用することと、
第２の関数Ｃ（ｐ）＝２ｐ／（ｐ2＋１）を適用することと
によって、前記信号値のそれぞれを前記円弧上の各点にマッピングするように構成される、請求項１４に記載の装置。
【請求項１７】
前記判定ユニットは、
前記計算された平均値を閾値と比較するように動作可能な比較器を有する検出ユニット、及び
前記計算された平均値を前記信号領域内の値に逆マッピングするように動作可能な逆マップ器と、該逆マッピングされた値を閾値と比較するように動作可能な比較器とを有する検出ユニット、
のうちの一方を含む、請求項１１〜１６のいずれか一項に記載の装置。
【請求項１８】
前記円弧上の前記マッピングされた点の密度の指標を計算するように動作可能な密度指標計算器をさらに含み、
前記判定ユニットは、前記求められた密度の指標を用いて前記閾値を設定すること、及び、前記求められた密度の指標を用いて前記閾値と比較される前記平均値を調整することのうちの少なくとも一方を含む処理を実行するように構成される、請求項１７に記載の装置。
【請求項１９】
前記密度指標計算器は、前記円弧上の前記マッピングされた点の前記密度の指標を、前記第１の平均値及び前記第２の平均値に基づいて
【数２】
に従って計算するように構成され、
ただし、ＤCCは前記密度の指標であり、ＸKは前記第１の平均値であり、ＹKは前記第２の平均値である、請求項１２〜１６のいずれか一項に従属するときの請求項１８に記載の装置。
【請求項２０】
前記密度指標計算器は、前記円弧密度の指標を用いて前記閾値を設定するように構成され、
前記密度の指標を密度閾値と比較することと、
前記密度の指標が前記密度閾値を超えるときに、予め設定された値を小さくすることによって、前記密度の指標と前記密度閾値との間の関係に基づいて前記閾値を設定することと
によって、前記閾値を設定する、請求項１８又は１９に記載の装置。
【図１】
【図２】
【図３】
【図４】
【図５】
【図６】
【図７】
【図８】
【図９】
【図１０】
【図２】
【図３】
【図４】
【図５】
【図６】
【図７】
【図８】
【図９】
【図１０】
【公表番号】特表２０１１−５２１２４５（Ｐ２０１１−５２１２４５Ａ）
【公表日】平成２３年７月２１日（２０１１．７．２１）
【国際特許分類】
【出願番号】特願２０１１−５０９９５８（Ｐ２０１１−５０９９５８）
【出願日】平成２１年５月１９日（２００９．５．１９）
【国際出願番号】ＰＣＴ／ＥＰ２００９／０５６０６２
【国際公開番号】ＷＯ２００９／１４１３４１
【国際公開日】平成２１年１１月２６日（２００９．１１．２６）
【出願人】（５０１２５３３１６）ミツビシ・エレクトリック・アールアンドディー・センター・ヨーロッパ・ビーヴィ (77)
【氏名又は名称原語表記】ＭＩＴＳＵＢＩＳＨＩ ＥＬＥＣＴＲＩＣ Ｒ＆Ｄ ＣＥＮＴＲＥ ＥＵＲＯＰＥ Ｂ．Ｖ．
【住所又は居所原語表記】２０ Ｆｒｅｄｅｒｉｃｋ Ｓａｎｇｅｒ Ｒｏａｄ， Ｔｈｅ Ｓｕｒｒｅｙ Ｒｅｓｅａｒｃｈ Ｐａｒｋ， Ｇｕｉｌｄｆｏｒｄ， Ｓｕｒｒｅｙ ＧＵ２ ５ＹＤ， Ｇｒｅａｔ Ｂｒｉｔａｉｎ
【Ｆターム（参考）】
【公表日】平成２３年７月２１日（２０１１．７．２１）
【国際特許分類】
【出願日】平成２１年５月１９日（２００９．５．１９）
【国際出願番号】ＰＣＴ／ＥＰ２００９／０５６０６２
【国際公開番号】ＷＯ２００９／１４１３４１
【国際公開日】平成２１年１１月２６日（２００９．１１．２６）
【出願人】（５０１２５３３１６）ミツビシ・エレクトリック・アールアンドディー・センター・ヨーロッパ・ビーヴィ (77)
【氏名又は名称原語表記】ＭＩＴＳＵＢＩＳＨＩ ＥＬＥＣＴＲＩＣ Ｒ＆Ｄ ＣＥＮＴＲＥ ＥＵＲＯＰＥ Ｂ．Ｖ．
【住所又は居所原語表記】２０ Ｆｒｅｄｅｒｉｃｋ Ｓａｎｇｅｒ Ｒｏａｄ， Ｔｈｅ Ｓｕｒｒｅｙ Ｒｅｓｅａｒｃｈ Ｐａｒｋ， Ｇｕｉｌｄｆｏｒｄ， Ｓｕｒｒｅｙ ＧＵ２ ５ＹＤ， Ｇｒｅａｔ Ｂｒｉｔａｉｎ
【Ｆターム（参考）】
[ Back to top ]