(57)【要約】【課題】受波器アレイを構成する複数の受波器で音波を受信することで生成された受波信号を処理して、方位−レベル表示のための信号を生成する処理において、微弱な信号を検出し易くする。【解決手段】複数の異なる方位についての整相出力（Ｒ）を時間積分（３０）して時間積分出力（Ｘ）を生成し、方位に対する時間積分出力の変化がより急峻となるように、鮮鋭化（４０）を施し、鮮鋭化の結果（Ｙ）を周波数方向に累加（５０）する。鮮鋭化はラプラシアン処理又はアンシャープマスク処理により行われる。【選択図】図２
2014038028.tif 000002
【特許請求の範囲】
【請求項１】
受波器アレイを構成する複数の受波器で音波を受信することで生成された受波信号を処理して、各方位における受波強度を示す信号を生成する信号処理装置において、
前記受波信号を複数の周波数ビンに分割し、各周波数ビンについての複数の異なる方位についての整相出力を得る整相処理部と、
前記整相出力を時間積分して時間積分出力を生成する時間積分部と、
前記時間積分出力を受け、方位に対する前記時間積分出力の変化がより急峻となるように、鮮鋭化を施す鮮鋭化部と、
前記鮮鋭化の結果を周波数方向に累加することで、前記各方位における受波強度に対応する信号を得る周波数累加部と
を備えた信号処理装置。
【請求項２】
前記鮮鋭化部が、
各方位及びそれに近接する方位における、前記時間積分出力に基づいて、鮮鋭化成分を求める鮮鋭化成分生成部と、
各方位における前記時間積分出力と、当該方位における前記鮮鋭化成分を合成する合成部とを有する
ことを特徴とする請求項１に記載の信号処理装置。
【請求項３】
前記時間積分部が、
第１の積分時間による第１の積分処理を行う第１の積分部と、
前記第１の積分時間よりも長い第２の積分時間による第２の積分処理を行う第２の積分部とを有し、
前記第１の積分処理の結果及び前記第２の積分処理の結果が前記時間積分出力として用いられ、
前記鮮鋭化部が、
前記第２の積分処理の結果を用いて前記鮮鋭化成分を求め、
前記第１の積分処理の結果と、前記鮮鋭化成分に所定の係数を掛けたものを加算又は減算することで前記鮮鋭化処理を行う
ことを特徴とする請求項１に記載の信号処理装置。
【請求項４】
前記周波数方向の累加の結果に対してレベル変換を行うレベル変換部と、
前記レベル変換の結果に対して対数変換を行い、さらに前記対数変換の結果を正規化する正規化部と、
前記正規化の結果の全方位にわたる平均値を求め、該平均値に基づいて、前記レベル変換におけるパラメータを調整するパラメータ調整部とをさらに有し、
前記パラメータの調整により、前記正規化の結果におけるコントラストを改善する
請求項１乃至３のいずれかに記載の信号処理装置。
【請求項５】
前記鮮鋭化部が、ラプラシアン処理によって前記鮮鋭化を行うことを特徴とする請求項１乃至４のいずれかに記載の信号処理装置。
【請求項６】
前記鮮鋭化部が、アンシャープマスク処理によって前記鮮鋭化を行うことを特徴とする請求項１乃至４のいずれかに記載の信号処理装置。
【請求項７】
受波器アレイを構成する複数の受波器で音波を受信することで生成された受波信号を処理して、各方位における受波強度を示す信号を生成する信号処理方法において、
前記受波信号を複数の周波数ビンに分割し、各周波数ビンについての複数の異なる方位についての整相出力を時間積分して時間積分出力を生成し、
前記時間積分出力に基づき、方位に対する前記時間積分出力の変化がより急峻となるように、鮮鋭化を施し、
前記鮮鋭化の結果を周波数方向に累加することで、前記各方位における受波強度に対応する信号を得る
ことを特徴とする信号処理方法。
【請求項８】
前記鮮鋭化が、
各方位及びそれに近接する方位における、前記時間積分出力に基づいて、鮮鋭化成分を求める処理と、
各方位における前記時間積分出力と、当該方位における前記鮮鋭化成分を合成する処理とで行われる
ことを特徴とする請求項７に記載の信号処理方法。
【請求項９】
前記時間積分が、
第１の積分時間による第１の積分処理と、
前記第１の積分時間よりも長い第２の積分時間による第２の積分処理とを含み、
前記第１の積分処理の結果及び前記第２の積分処理の結果が前記時間積分出力として用いられ、
前記鮮鋭化が、
前記第２の積分処理の結果を用いて前記鮮鋭化成分を求める処理と、
前記第１の積分処理の結果と、前記鮮鋭化成分に所定の係数を掛けたものを加算又は減算する処理とを含む
ことを特徴とする請求項７に記載の信号処理方法。
【請求項１０】
前記周波数方向の累加の結果に対してレベル変換を行うレベル変換を行い、
前記レベル変換の結果に対して対数変換を行い、さらに前記対数変換の結果を正規化し、
前記正規化の結果の全方位にわたる平均値を求め、該平均値に基づいて、前記レベル変換におけるパラメータを調整し、
前記パラメータの調整により、前記正規化の結果におけるコントラストを改善する
請求項７乃至９のいずれかに記載の信号処理方法。
【請求項１１】
前記鮮鋭化がラプラシアン処理によって行われることを特徴とする請求項７乃至１０のいずれかに記載の信号処理方法。
【請求項１２】
前記鮮鋭化がアンシャープマスク処理によって行われることを特徴とする請求項７乃至１０のいずれかに記載の信号処理方法。
【請求項１３】
請求項７乃至１２のいずれかの信号処理方法をコンピュータに実行させるためのプログラム。
【請求項１４】
請求項１３のプログラムを記録した、コンピュータで読み取り可能な記録媒体。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、ソーナーシステムにおける信号処理装置及び方法、並びにプログラム及び記録媒体に関する。本発明は、特に、ソーナーシステムにおいて整相出力の周波数−方位スペクトルを周波数方向に累加して画面に方位−レベル出力の表示をする際、周波数累加の前段で鮮鋭化処理を行い、従来方法と比べて微弱な信号をより検出し易くした信号処理に関するものである。
【背景技術】
【０００２】
従来技術における方位−レベル表示のための処理では、整相処理によって得られる整相出力（周波数及び方位の関数としての受波信号のレベル）に対して、時間積分処理及び周波数方向の累加を行うことによって、方位毎のレベルを示す表示（方位−レベル表示）のための処理出力を得ていた。
【先行技術文献】
【非特許文献】
【０００３】
【非特許文献１】Ｒｏｂｅｒｔ  Ｅ．  Ｚａｒｎｉｃｈ、“Ａ  Ｆｒｅｓｈ  Ｌｏｏｋ  ａｔ    ‘Ｂｒｏａｄｂａｎｄ’  Ｐａｓｓｉｖｅ  Ｓｏｎａｒ  Ｐｒｏｃｅｓｓｉｎｇ”，  Ａｄａｐｔｉｖｅ  Ｓｅｎｓｏｒ  Ａｒｒａｙ  Ｐｒｏｃｅｓｓｉｎｇ  Ｗｏｒｋｓｈｏｐ  （ＡＳＡＰ‘  ９９），  ＭＩＴ：  ＬＬ，  Ｍａｒｃｈ，  ｐｐ．  ９９−１０４  （１９９９）
【非特許文献２】Ｓｈｅｆｅｎｇ  Ｙａｎ  ａｎｄ  Ｙｕａｎｌｉａｎｇ  Ｍａ，  “Ｈｉｇｈ−ｒｅｓｏｌｕｔｉｏｎ  ｂｒｏａｄｂａｎｄ  ｂｅａｍｆｏｒｍｉｎｇ  ａｎｄ  ｄｅｔｅｃｔｉｏｎ  ｍｅｔｈｏｄｓ  ｗｉｔｈ  ｒｅａｌ  ｄａｔａ”，  Ａｃｏｕｓｔ．  Ｓｃｉ．  ＆  Ｔｅｃｈ，  ２５，  １，  ｐｐ．  ７３−７６  （２００４）
【発明の概要】
【発明が解決しようとする課題】
【０００４】
しかしながら、このような従来技術では、周波数方向の累加によって、微弱な信号が周波数の異なる強い信号に埋もれてしまう場合がある。図１（ａ）及び（ｂ）は、従来技術による処理の結果の概念図である。図１（ａ）は、整相出力を示し、図１（ｂ）は、周波数方向の累加後の出力を示している。図１（ａ）では、整相出力（方位一周波数スペクトル）において、信号Ａ、Ｂ及びＣの３つが存在していることが示されている。信号Ａはレベルが高いため、背景雑音が大きくても、周波数方向の累加の処理後に、図１（ｂ）に示すように、信号の存在を認識することができる。また、信号Ｂはレベルが大きくないが、周囲にレベルの大きな信号が存在していないため、周波数方向の累加の処理後においても存在を認識することができる。一方、信号Ｃは整相処理後の方位−周波数スペクトルでは認識できるが、周波数方向の累加によって信号Ａに埋もれてしまい、周波数方向の累加の処理後には、信号Ｃの存在を認識することができない。
【０００５】
このように、従来技術による処理においては、微弱な信号（図１（ａ）における信号Ｃ）が、周波数方向の累加後に他のレベルの高い信号（図１（ａ）における信号Ａ）または背景雑音に埋もれてしまい、方位−レベル表示において認識できないという問題があった。
【０００６】
そこで、本発明では、方位−レベル表示において周波数累加後でも微弱な信号を認識できるようにすることを目的とする。
【課題を解決するための手段】
【０００７】
本発明の信号処理装置は、
受波器アレイを構成する複数の受波器で音波を受信することで生成された受波信号を処理して、各方位における受波強度を示す信号を生成する信号処理装置において、
前記受波信号を複数の周波数ビンに分割し、各周波数ビンについての複数の異なる方位についての整相出力を得る整相処理部と、
前記整相出力を時間積分して時間積分出力を生成する時間積分部と、
前記時間積分出力を受け、方位に対する前記時間積分出力の変化がより急峻となるように、鮮鋭化を施す鮮鋭化部と、
前記鮮鋭化の結果を周波数方向に累加することで、前記各方位における受波強度に対応する信号を得る周波数累加部と
を備えたことを特徴とする。
【発明の効果】
【０００８】
本発明によれば、方位−レベル表示において周波数累加後でも微弱な信号の認識がし易くなる。
【図面の簡単な説明】
【０００９】
【図１】（ａ）及び（ｂ）は、従来技術における整相出力と、方位−レベル表示のための処理出力との関係を示すグラフである。
【図２】本発明の実施の形態１の信号処理装置を示すブロック図である。
【図３】（ａ）〜（ｃ）は、ラプラシアンを用いた鮮鋭化処理において生成される信号の例を示す図である。
【図４】従来例により生成される表示出力の一例を示す図である。
【図５】実施の形態１により生成される表示出力の一例を示す図である。
【図６】本発明の実施の形態２の信号処理装置を示すブロック図である。
【図７】実施の形態２により生成される表示出力の一例を示す図である。
【図８】本発明の実施の形態３の信号処理装置を示すブロック図である。
【図９】（ａ）及び（ｂ）は、実施の形態３によるレベル調整による、平均値の変化を示す図である。
【図１０】実施の形態３により生成される表示出力の一例を示す図である。
【発明を実施するための形態】
【００１０】
実施の形態１．
図２は、本発明の実施の形態１の信号処理装置を示す。図示の信号処理装置は、整相処理部２０と、時間積分部３０と、鮮鋭化部４０と、周波数累加部５０と、対数変換部６０とを備える。
【００１１】
整相処理部２０は、受波器アレイ１０からの受波レベルを表す信号（受波信号）を受ける。
受波器アレイ１０は、例えば直線状又は円弧状に配列され、それぞれ音波を受信することで受波信号を生成する複数の受波器で構成される。整相処理部２０は、受波器アレイ１０からの受波信号に対して整相処理を行う。
【００１２】
整相処理は、受波信号の位相もしくは遅延時間を調整して、方位毎の受波ビームを形成する処理である。整相処理においては、例えば、受波信号を、高速フーリエ変換（ＦＦＴ）などによって周波数領域の信号に変換し、周波数ビン毎に、各受波器からの信号に対して位相調整を行う。
【００１３】
各時刻ｔにおいて、整相処理の結果として得られる信号（整相出力）Ｒｋ，ｆ（ｔ）、受波器アレイ１０を中心とする各方位ｋ（ｋ＝１〜Ｋ）（Ｋは方位の総数）における、周波数ビンｆ（ｆ＝１〜Ｆ）毎の各時刻における受波強度（レベル）を表すものである。
【００１４】
時間積分部３０は、整相処理部２０からの整相出力Ｒｋ,ｆ（ｔ）を時間積分する。時間積分部３０による時間積分は、整相出力のバラツキを抑え、雑音成分を抑制する役割を持つ。即ち、目標とする信号成分は一般に比較的安定したレベルを有するのに対して雑音成分のレベルは時間によって変動するので、時間積分することで雑音成分を抑制し、信号成分を抽出することができる。
時間積分は、各方位における、各周波数ビンについて行われ、その結果（積分出力）も各方位の周波数ビン毎に得られる。積分時間をＴとする場合、積分出力は、例えば、現在の時刻ｔを基準として、時間Ｔ前から現在までの整相出力の平均を求めることで得られる。即ち、積分出力Ｘｋ,ｆ（ｔ）は、下記の式（１）で表される。
【００１５】
【数１】
2014038028.tif 000003
【００１６】
但し、時刻ｔは、整相出力Ｒｋ,ｆ（ｔ）が１回得られる毎に１ずつ増加する値で表されている。なお以下では、各時刻（ｔ）における値であることを強調する必要がない限り、時刻（ｔ）を省略して単にＸｋ，ｆで表す。他の値についても同様である。
なお，積分出力の求め方としては、式（１）で表された平均の代わりに、総和でもよいし、加重平均を用いてもよい。
【００１７】
鮮鋭化部４０は、時間積分部３０の出力（積分出力）Ｘｋ，ｆに対して鮮鋭化処理を行う。鮮鋭化処理も、各周波数ビンについて行われる。鮮鋭化処理は、方位ｋに対する積分出力Ｘｋ，ｆの変化がより急峻となるように行われるものであり、例えば、各方位について、当該方位についての積分出力と、当該方位に近接する１又は２以上の方位についての積分出力を用いて行われる。鮮鋭化処理の結果Ｙｋ，ｆは、各方位ｋの受波強度を表す。
【００１８】
周波数累加部５０は、鮮鋭化部４０の出力Ｙｋ，ｆに対して周波数方向の累加を行う。即ち、各方位について、それぞれの周波数ビンについての鮮鋭化処理の結果Ｙｋ，ｆの総和Ｚｋを求める。この演算は下記の式（２）で表される。
【００１９】
【数２】
2014038028.tif 000004
【００２０】
このような処理により得られる、周波数累加の結果Ｚｋは、各方位の全周波数帯域における、受波強度に対応するものであるが、時間積分し、かつ鮮鋭化を行った後に出力を累加することで得られた値を有する。
【００２１】
対数変換部６０は、周波数累加部５０の出力（各方位についての累加結果）Ｚｋを対数に変換する。例えば、１０を底とする対数への変換が行われる。この対数変換は下記の式（３）により表される。
【００２２】
【数３】
2014038028.tif 000005
【００２３】
対数変換部６０の出力Ｌｋも、各方位の全周波数帯域における受波強度に対応する。
対数変換部１０４の出力が、本実施の形態の鮮鋭化処理システムの出力となる。
【００２４】
従来技術による処理と比べ、時間積分部３０の出力に対して鮮鋭化部４０による鮮鋭化処理を適用することが本発明の特徴である。
鮮鋭化部４０は、鮮鋭化成分生成部４１と、合成部４２とを有する。
【００２５】
鮮鋭化成分生成部４１は、鮮鋭化成分を、例えばラプラシアンを用いて求める。
図３（ａ）〜（ｃ）は、実施の形態１でのラプラシアンを用いた鮮鋭化処理の方法を示している。
【００２６】
この処理は、ある方位ｋにおける、ある周波数の時間積分出力Ｘｋ，ｆに対して適用される。ここで、ｋは方位番号を表す整数である。
【００２７】
図３（ａ）では、方位ｋを横軸とし、ある周波数の時間積分出力Ｘｋ，ｆを縦軸とする曲線において、ピークが１つある場合を例に示す。時間積分出力Ｘｋ，ｆに対するラプラシアン処理の結果∇２Ｘｋ，ｆは、下記の式（４）で表される。
【００２８】
【数４】
2014038028.tif 000006
【００２９】
この処理によって、方位−レベル表示出力（各方位におけるレベルを表す値）Ｘｋ，ｆに対する、方位方向の２階微分値が得られる。この処理を図３（ａ）に示される値Ｘｋ，ｆに対して適用すると、図３（ｂ）に示される処理結果が得られる。図３（ｂ）の処理結果は、方位に対するレベルの変化を表す曲線（図３（ａ））のピークの位置で極小値、傾斜が変化しない区間で０となる。このラプラシアン処理結果∇２Ｘｋ，ｆに係数αを掛けて、時間積分処理の結果Ｘｋ，ｆから減算することで鮮鋭化出力Ｙｋ，ｆを求める。即ち、鮮鋭化出力Ｙｋ，ｆは下の式（５）で表される。
【００３０】
【数５】
2014038028.tif 000007
【００３１】
ここで、αは鮮鋭化の程度を調整するパラメータであり、αが大きいほど鮮鋭化の程度が強く、α＝０のとき、
Ｙｋ，ｆ＝Ｘｋ，ｆ
となり、鮮鋭化が一切行われないのと同じで、従来技術による処理結果と同じ結果となる。
【００３２】
図３（ａ）に示される時間積分処理の結果に対して、式（５）で表される鮮鋭化処理を行った場合に得られる処理結果が図３（ｃ）に示されている。
時間積分出力Ｘｋ，ｆ（図３（ａ））にピークがある場合、ラプラシアン処理による鮮鋭化成分∇２Ｘｋ，ｆ（図３（ｂ））は通常、上記のピークと同じ方位においてより鋭く下に（ピークとは逆の向きに）凸となる。従って、鮮鋭化成分に係数αを掛けたものを減算することで得られる鮮鋭化処理の結果Ｙｋ，ｆ（図３（ｃ））は、時間積分出力Ｘｋ，ｆ（図３（ａ））よりも、ピークがより鋭く強調されものとなる。
【００３３】
図４及び図５は、シミュレーションの結果を示す。即ち、疑似信号に疑似雑音を重畳した疑似的なアレイの受波器出力に対して整相処理を行い、その整相出力に対する処理の結果を用いたＢＴＲ（Bearing Time Recording）を示す。ＢＴＲ表示では、縦軸が時間、横軸が方位、濃淡がレベルの強さ（濃いほどレベルが高いことを示す）を表す。なお、疑似信号として、方位番号１４４に信号Ｓ１、方位番号１５０に信号Ｓ２を生成した。
【００３４】
図４は、比較のための従来技術による処理の結果を示し、図５が、実施の形態１の処理の結果を示す。従来処理による結果では、方位番号１４４の信号Ｓ１は確認できるが、方位番号１５０の信号Ｓ２は確認しにくい。一方、実施の形態１による結果では、方位番号１４４の信号Ｓ１のみならず、方位番号１５０の信号Ｓ２をも確認することができ、信号が検出し易くなっていることがわかる。
【００３５】
実施の形態２.
図６は、本発明の実施の形態２の信号処理装置を示す。
実施の形態１との違いは、図２の時間積分部３０の代わりに、時間積分部３５を設けたこと、及び図２の鮮鋭化部４０の代わりに鮮鋭化部４５を設けたことである。時間積分部３５は、第１の積分部３６と第２の積分部３７とを有する。
【００３６】
第１の積分部３６及び第２の積分部３７は、図２の時間積分部３０と同じく、整相出力のバラツキを抑え、雑音を除去することを目的とする。第１の積分部３６の積分時間Ｔ１は、図２の時間積分部３０の積分時間と同程度である。第２の積分部３７の積分時間Ｔ２は、第１の積分部３６の積分時間Ｔ１よりも長い。
【００３７】
鮮鋭化部４５は、実施の形態１における鮮鋭化部４０と同様であるが、第２の積分部３７の出力と、第１の積分部３６の出力とを用いて鮮鋭化処理を行うものであり、鮮鋭化成分生成部４６と、合成部４７とを有する。
鮮鋭化成分生成部４６は、実施の形態１における鮮鋭化成分生成部４１と同じであるが、第２の積分部３７の出力に対して鮮鋭化を行う。
合成部４７は、第１の積分部３６の出力と鮮鋭化成分生成部４６の出力を合成する。この合成は、鮮鋭化成分生成部４６の出力に所定の係数αを掛けたものを第１の積分部３６の出力から減算することで行われる。
【００３８】
第１の積分部３６の出力をＸＴ１，ｋ，ｆと表し、第２の積分部３７の出力をＸＴ２，ｋ，ｆと表す。
鮮鋭化成分生成部４６の出力は下式（６）で表される。
【００３９】
【数６】
2014038028.tif 000008
【００４０】
合成部４７では、第１の積分部３６の出力から、式（６）に示される、鮮鋭化成分∇２ＸＴ２，ｋ，ｆに係数αを掛けたものを減算することで鮮鋭化部４５の出力Ｙｋ，ｆを得る。即ち、鮮鋭化部４５の出力Ｙｋ，ｆは、次の式（７）で表される。
【００４１】
【数７】
2014038028.tif 000009
【００４２】
Ｔ１＝Ｔ２の場合は、式（５）で表される実施の形態１と同じ結果となる。
【００４３】
ラプラシアンを用いた鮮鋭化処理では、微分処理を行うため、雑音の影響によって時間積分出力のバラツキが大きい場合、バラツキの影響を強く受ける可能性がある。そこで、第２の積分部３７では、第１の積分部３６よりも長い積分時間を設定し、雑音の影響によるバラツキをさらに抑制した時間積分出力を得て、それに対してラプラシアン処理を行うこととしている。
【００４４】
図７は、実施の形態２を適用した場合の結果を示している。ここで、第２の積分部３７における積分時間Ｔ２を、第１の積分部３６における積分時間Ｔ１の１０倍に設定した場合を示す。図５と比べても、信号周辺の雑音成分が抑制されていることに加え、信号がより鮮鋭になっており、全体的に信号の見えやすさが実施の形態１よりも向上していることが確認できる。
【００４５】
実施の形態３．
図８は、本発明の実施の形態３の信号処理装置を示す。
実施の形態２との違いは、周波数累加部５０の後段に、レベル変換部７０を設け、さらに対数変換部６０の代わりに、対数変換部６５を設け、対数変換部６５の出力を受ける正規化部８０、及び正規化部８０の出力を受けるパラメータ調整部９０を設けたことである。
【００４６】
実施の形態２では、第１及び第２の積分部３６、３７の出力によっては、周波数累加部５０の出力に、
（１）負の値が含まれる場合や、
（２）極めて小さい数値が含まれる場合がある。
例えば、周波数累加部５０の出力に負の値が含まれる場合、対数変換部１０４での変換処理でエラーが生じ、出力が得られない。また、周波数累加部５０の出力に極めて小さい値が含まれている場合、対数変換は可能であるが、対数変換の結果得られる出力の大きさの分布が偏り、その結果として、処理結果のコントラストが悪くなることがある。
【００４７】
そこで、実施の形態３では、周波数累加部５０の出力Ｚｋに対してレベル変換部７０でレベル変換を行う。この処理では、周波数累加部５０の出力Ｚｋに対して、次式（８）で表される処理を行い、レベル変換出力Ｔｋを得る。
【００４８】
【数８】
2014038028.tif 000010
【００４９】
式（８）で、ｍａｘ（Ｚｋ）及びｍｉｎ（Ｚｋ）は周波数累加の出力Ｚｋにおける最大値及び最小値（各時点で、出力Ｚｋ(ｋ＝１〜Ｋ)が実際の取った値のうちの最大値及び最小値）、ａ及びｂはレベル変換処理のパラメータであり、ａ＞０、ｂ≧０の値が設定される。Ｔｋは
ｂ≦Ｔｋ≦ａ＋ｂ
を満たす。即ち、式（８）で表されるレベル変換によって、周波数累加の出力Ｚｋは、ｂからａ＋ｂの間の値（Ｔｋ）に変換される。
このレベル変換処理によって、対数変換前に負の値を取ることはなくなる。
【００５０】
対数変換部６５は、レベル変換出力Ｔｋの対数を求める。例えば、１０を底とする対数への変換が行われる。この対数変換は下記の式（９）により表される。
【００５１】
【数９】
2014038028.tif 000011
【００５２】
正規化部８０は、対数変換部６５の出力の正規化を行って画面表示出力Ｕｋを求める。正規化は下記の式（１０）で表される。
【００５３】
【数１０】
2014038028.tif 000012
ｍａｘ（Ｓｋ）は、対数Ｓｋの最大値、ｍｉｎ（Ｓｋ）は、対数Ｓｋの最小値である。
【００５４】
式（１０）で得られる画面表示出力Ｕｋは、
０≦Ｕｋ≦１
を満たす。この画面表示出力Ｕｋが方位−レベル表示に用いられる。
【００５５】
パラメータ調整部９０では、画面表示出力Ｕｋの平均的レベル及びコントラストを調整するための処理を行う。
パラメータ調整部９０は、平均値計算部９１と、パラメータ決定部９２とを有する。
平均値計算部９１は、画面表示出力Ｕｋの、全方位についての平均値Ｕａｖｇを求める。平均値Ｕａｖｇは下式（１１）で表される。
【００５６】
【数１１】
2014038028.tif 000013
【００５７】
パラメータ決定部９２は、上記の式（１１）で求めた平均値Ｕａｖｇに応じて、式（８）のパラメータｂを次式（１２）で表されるように変化させる。
【００５８】
【数１２】
2014038028.tif 000014
【００５９】
ここで、Ｍはコントラスト調整用パラメータであり、画面表示出力Ｕｋの平均値Ｕａｖｇはこの値Ｍに近づくように変化する。μはパラメータｂを更新する量である。
【００６０】
図９（ａ）及び（ｂ）は、上記の処理による効果を示す。画面表示出力Ｕｋの平均値Ｕａｖｇが設定値Ｍ以上である場合に、式（８）におけるパラメータｂが大きくされ、この結果、レベル変換出力Ｔｋの取る範囲が大きい方へ変化する。このため、対数変換の出力Ｓｋのレベルの分布が変化し、その結果として画面表示出力Ｕｋのコントラストが変化する。この処理を逐次的に繰り返すことによって、画面表示出力Ｕｋの平均値Ｕａｖｇが設定値Ｍに近づき、処理結果のコントラストが改善される。
【００６１】
図１０は、実施の形態３を適用した場合の結果を示している。
図１０は、Ｍを、画面表示出力ＵＫの最大値ｍａｘ（Ｕｋ）＝１の１／５に設定したときの結果である。図５及び図７と比べ、信号Ｓ２がより鮮明になっており、本実施の形態の効果が確認できる。
【００６２】
なお、実施の形態３を実施の形態２の変形例として説明したが、実施の形態１に対して同様の変形を加えても良い。
【００６３】
変形例１．
実施の形態３では、周波数累加の出力に対してレベル変換を行っているが、鮮鋭化処理の出力（周波数累加の前の出力）に対してレベル変換を行っても良いし、単に鮮鋭化処理の出力の負値を０以上の値へ置き換えても良い。鮮鋭化処理の出力に負の値が含まれる場合、周波数累加の処理によって、目標が存在している方位のレベルが小さくなる可能性があるが、レベル変換をすることでそのような事態を避け得るからである。
【００６４】
変形例２．
これまでの説明では、鮮鋭化処理にラプラシアン処理を適用してきたが、鮮鋭化処理には他の方法を適用しても良い。例えば、アンシャープマスク処理を適用しても良い。以下、実施の形態１において、ラプラシアン処理の代わりに、アンシャープマスク処理を行う場合について説明する。
【００６５】
アンシャープマスク処理により生成される鮮鋭化成分は下記の式（１３）で表される。
【００６６】
【数１３】
2014038028.tif 000015
【００６７】
鮮鋭化成分生成部４１では、上記の式（１３）で表される鮮鋭化成分を計算し、合成部４２では、鮮鋭化成分生成部４１で生成された鮮鋭化成分と、時間積分出力Ｘｋ，ｆとを合成して、鮮鋭化出力を求める。この合成においては、鮮鋭化成分に係数αを掛けたものを、時間積分出力Ｘｋ，ｆに加算する。
その結果、次式（１４）で表される鮮鋭化出力Ｙｋ，ｆが得られる。
【００６８】
【数１４】
2014038028.tif 000016
【００６９】
なお、式（１４）における各方位ｋ及びそれに隣接する方位（ｋ−１）、（ｋ＋１）における値を平均する部分、すなわち、
【数１５】
2014038028.tif 000017
の部分は、重み付きの平均に置き換えてもよい。さらに、ｋ番目の方位の処理に使用するデータは、該方位（ｋ番目の方位）に直に隣接する１つずつ(ｋ−１番目とｋ＋１番目)である必要はなく、用途に応じて増やすことも可能である。
【００７０】
以上、本発明を装置として説明したが上記装置で実施される方法も本発明の一部を成す。また、上記の装置の一部又は全部はソフトウェアにより、即ちプログラムされたコンピュータで実現することができ、コンピュータに上記の処理の実行させるためのプログラム、並びに該プログラムを記録したコンピュータで読み取り可能な記録媒体も本発明の一部を成す。
【符号の説明】
【００７１】
１０  受波器アレイ、  ２０  整相処理部、  ３０、３５  時間積分部、  ３６  第１の積分部、  ３７  第２の積分部、  ４０、４５  鮮鋭化部、  ４１、４６  鮮鋭化成分生成部、  ４２、４７  合成部、  ５０  周波数累加部、  ６０、６５  対数変換部、  ７０  レベル変換部、  ８０  正規化部、  ９０  パラメータ調整部、  ９１  平均値計算部、  ９２  パラメータ決定部。
【図１】
2014038028.tif 000018
【図２】
2014038028.tif 000019
【図３】
2014038028.tif 000020
【図４】
2014038028.tif 000021
【図５】
2014038028.tif 000022
【図６】
2014038028.tif 000023
【図７】
2014038028.tif 000024
【図８】
2014038028.tif 000025
【図９】
2014038028.tif 000026
【図１０】
2014038028.tif 000027