JP5182676B2 - 合成開口処理システム及び合成開口処理方法 - Google Patents

Description

本発明は、レーダーやソーナーなどの合成開口処理システム及び合成開口処理方法に関するものである。

ソーナーやレーダーは遠隔目標を広範囲に調査できる反面、遠距離の探査ではソーナーやレーダーを搭載するプラットフォームの進行方向（アジマス方向）の分解能が低下するという欠点をもつ。この欠点を解消するために、扱いやすい短い受信アレイを用いて、受信アレイを仮想的に延長することでアジマス分解能を向上させる合成開口処理が知られている。

合成開口処理では、プラットフォームが一定速度で移動しつつ、一定間隔でパルスが地海面に対して繰り返し放射され、複数の反射波が受信される。そしてこれらの複数の反射波を処理することにより、レンジ方向（パルスが放射された地海面上に沿った方向）及びアジマス方向の分解能を向上させたソーナー画像又はレーダー画像が得られる。

合成開口処理の手法としては、１次元アルゴリズムと２次元アルゴリズムが知られている（例えば非特許文献１を参照）。１次元アルゴリズムでは、まず、従来のパルス圧縮技術に基づいて、反射波が時刻方向にパルス圧縮される（これをレンジ圧縮という）。次いで、レンジカーバチャと呼ばれるプラットフォームの移動の影響による反射波の記録時刻のずれが修正される。そして、複数の反射波の波形振幅が加算されて、レンジ方向及びアジマス方向の分解能を向上させたソーナー画像又はレーダー画像として出力される。

２次元アルゴリズムでは、受信した複数の反射波と、反射波の受信可能範囲の複数の所定位置で反射した場合の理想的な受信波に基づき作成された参照信号との２次元相関演算が行われ、受信可能範囲内の複数の所定位置ごとに相関値が算出される。これらの相関値データに基づいて、レンジ方向及びアジマス方向の分解能を向上させたソーナー画像又はレーダー画像が作成される。
大内和夫、"リモートセンシングのための合成開口レーダの基礎"、東京電機大学出版局、２００４

合成開口処理を行いつつレーダーやソーナーの実用性を保つためには、リアルタイムでレーダー画像やソーナー画像を得られることが望ましく、このため、合成開口処理による計算機の演算負荷を低減させたいとの要望がある。

しかしながら、従来の１次元アルゴリズムでは、アジマス方向の処理は単純な総和演算なので演算負荷は小さいが、レンジカーバチャ修正の処理のために反射波のデータ補間が必要であり、条件によっては精度誤差が大きくなるという問題がある。

また、２次元アルゴリズムでは、反射波をそのまま用いて、参照信号と２次元相関演算を行うので精度良く合成開口処理が行えるが、膨大な数の参照信号との相関演算が必要なため、演算負荷が大きくなるという問題がある。

そこで、本発明は、上記の問題点を解決し、合成開口処理を高精度に行うことができ、かつ、演算負荷を低減させることができる合成開口処理システム及び合成開口処理方法を提供することを目的とする。

上記課題を解決するために、本発明の合成開口処理システムは、周波数が時間変化する周波数変調波を生成して、複数の計測位置から照射領域に照射する送信手段と、送信手段により照射された複数の周波数変調波による複数の反射波を受信する受信手段と、受信手段により受信された複数の反射波のそれぞれにレンジ圧縮を行い、複数のシンク関数からなる受信データを生成するレンジ圧縮処理手段と、周波数変調波が前記照射領域内の複数の所定位置で反射した場合の理想的な受信波をレンジ圧縮した複数のシンク関数からなる複数のモデルデータと、レンジ圧縮処理手段により生成された受信データとに基づいて、モデルデータのそれぞれと受信データとの相関度合いを示す相関値データを算出する相関値算出手段と、相関値算出手段により算出された相関値データを出力する出力手段と、を備えている。

同様に、上記課題を解決するために、本発明の合成開口処理方法は、周波数が時間変化する周波数変調波を生成して、複数の計測位置から照射領域に照射する送信ステップと、送信ステップにおいて照射された複数の周波数変調波による複数の反射波を受信する受信ステップと、受信ステップにおいて受信された複数の反射波のそれぞれにレンジ圧縮を行い、複数のシンク関数からなる受信データを生成するレンジ圧縮処理ステップと、周波数変調波が照射領域内の複数の所定位置で反射した場合の理想的な受信波をレンジ圧縮した複数のシンク関数からなる複数のモデルデータと、レンジ圧縮処理ステップにおいて生成された受信データとに基づいて、モデルデータのそれぞれと前記受信データとの相関度合いを示す相関値データを算出する相関値算出ステップと、相関値算出ステップにおいて算出された相関値データを出力する出力ステップと、を備えている。

これらの構成により、まずレンジ圧縮を行って、反射波に比べてデータ数が少ないシンク関数からなる受信データを用いて相関値データを算出するので、演算負荷を低減することができる。また、従来の２次元アルゴリズムと同様に、受信データとモデルデータとの相関値データを算出するので、従来の１次元アルゴリズムのレンジカーバチャ修正が不要となり、精度誤差が大きくなるのを回避できる。なお、相関値データとは、周波数変調波の照射領域内の複数の所定位置において、周波数変調波が実際に反射したか否かを示す信頼度の指標であり、例えばあるモデルデータとサンプルデータとの相関値が大きいということは、このモデルデータの対応する所定位置に実際に物体が存在し、周波数変調波がこの位置で反射した可能性が高いことを表すものである。

また、本発明の合成開口処理システムは、送信手段により生成された周波数変調波に応じて、相関値算出手段で用いられるモデルデータを照射領域内の複数の所定位置ごとに算出するモデルデータ算出手段をさらに備えるのが好適である。

この構成により、実際に使用する周波数変調波に応じてモデルデータが生成されるので、合成開口処理を確実に実施できる。

また、本発明の合成開口処理システムにおいて、モデルデータ算出手段が、モデルデータを構成するシンク関数ごとに、ピークを挟んだ所定の間隔の２点のデータを抽出し、これらの抽出されたデータセットをモデルデータとして算出し、相関値算出手段が、モデルデータ算出手段により算出された複数のモデルデータのそれぞれについて、このモデルデータを構成するシンク関数ごとの２点のデータの位置情報に応じて、受信データを構成するシンク関数の中からそれぞれ２点のデータをサンプルデータとして選択する選択手段と、選択手段により選択された受信データのサンプルデータと、複数のモデルデータのデータセットとの間でそれぞれ相互相関演算を行い、相関値データを算出する相互相関演算手段と、をさらに備えるのが好適である。

この構成により、受信データの中からサンプルデータを抽出して相互相関演算が行われるので、相関値データを算出するためのデータ数がさらに少なくなり、演算負荷をさらに低減することができる。

また、本発明の合成開口処理システムにおいて、モデルデータ算出手段が、モデルデータを構成するシンク関数ごとに、ピークを挟んだ所定の間隔の２点のデータを抽出し、これらの抽出されたデータセットをモデルデータとして算出し、相関値算出手段が、モデルデータ算出手段により算出された複数のモデルデータのそれぞれについて、このモデルデータを構成するシンク関数ごとの２点のデータの位置情報に応じて、受信データを構成するシンク関数の中からそれぞれ２点のデータをサンプルデータとして選択する選択手段と、選択手段により選択された受信データのサンプルデータの総和演算を行い、相関値データを算出する総和演算手段と、をさらに備えるのが好適である。

この構成により、受信データの中からサンプルデータを抽出して相関値データが算出されるので、相関値データを算出するためのデータ数がさらに少なくなる。また、サンプルデータを抽出した後にこれらのデータの総和演算を行うだけで相関値データが算出されるので、演算負荷をさらに低減できる。

本発明の合成開口処理システム及び合成開口処理方法によると、合成開口処理を高精度に行うことができ、かつ、演算負荷を低減させることができる。

合成開口ソーナーの概略を示す図である。 本発明の第１の実施形態に係る合成開口処理システムの構成を示す図である。 本実施形態で照射されるチャープ波の一例を示す図である。 本実施形態で受信される反射波の一例を示す図である。 レンジ圧縮処理部により作成された受信データの一例を示す図である。 モデルデータ算出部により算出される複数のモデルデータの一例を示す図である。 サンプルデータ選択部によりモデルデータごとに受信データから抽出されるサンプルデータの一例を示す図である。 相互相関演算部により算出された相関値データの一例を示す図である。 従来のソーナーによるソーナー画像（ａ）と、画像出力部において出力されたソーナー画像（ｂ）の一例を示す図である。 本実施形態の合成開口処理システムによる処理を示すフローチャートである。 モデルデータ算出部においてシンク関数状の波形から選択される２点と、この波形のパルスピークとの関係を示す図である。 本発明の第２の実施形態に係る合成開口処理システムの構成を示す図である。 本発明の第３の実施形態に係る合成開口処理システムの構成を示す図である。

符号の説明

１０、２０、３０…合成開口処理システム、１１…計測条件入力部、１２…信号送信部、１３…アンテナ、１４…信号受信部、１５…レンジ圧縮処理部、１６…モデルデータ算出部、１７…サンプルデータ選択部、１８、２８…相互相関演算部、１９…画像出力部、２６…モデルデータ算出部、３８…総和演算部。

以下、図面を参照して本発明に係る合成開口処理システム及び合成開口処理方法の好適な実施形態について説明する。なお、図面の説明においては、同一要素には同一符号を付し、重複する説明を省略する。

（合成開口ソーナーの概要）
まず、図１を参照して、本発明に係る合成開口システムが適用される合成開口ソーナー（Ｓｙｎｔｈｅｔｉｃ Ａｐｅｒｔｕｒｅ Ｓｏｎｅｒ：ＳＡＳ）の概要について説明する。

従来のソーナー１は、船や潜水艦などのプラットフォーム２に搭載され、プラットフォーム２が一方向に進行しながら、海底面などの照射領域３に対してパルス音波を照射する。そして、何秒後にどれくらいの強度の反射波が戻ってくるかを計測することで、ソーナー画像を生成し、パルス音波放射方向（レンジ方向）にある物体を検出する。照射するパルス音波には、主に周波数変調波の一種であるチャープパルス（チャープ波）が用いられる。これにより、計測した反射波に周知のパルス圧縮技術を適用して、ソーナー画像のレンジ方向の分解能を向上させることができる。

合成開口ソーナー４は、プラットフォーム２の進行方向（アジマス方向）の複数の位置において、従来のソーナー１で受信された複数の反射波を合成処理することによって、ソーナー１の受信アレイをアジマス方向へ仮想的に延長させたものである。つまり、合成開口ソーナー４は、図１に示すように、従来のソーナー１に比べて大きなアンテナを備えているのと同等であり、従来のソーナー１によるアジマス方向の複数の位置における照射領域３ａ〜３ｃの共通部分５について、アジマス方向の分解能を向上させることができる。

ここで、本発明に係る合成開口システムとは、図１に示す合成開口ソーナー４を指すものとする。また、以下の説明で用いる「アジマス方向」とは、図１に示すように、プラットフォーム２の進行方向を指すものとする。「レンジ方向」とは、アジマス方向と直交し、プラットフォーム２からパルス音波を照射する方向を指すものとする。

（第１実施形態）
次に、本発明に係る合成開口システムの第１の実施形態について説明する。図２に本発明の第1の実施形態に係る合成開口処理システムの構成を示す。図２に示すように、本実施形態の合成開口処理システム１０は、計測条件入力部１１、信号送信部（送信手段）１２、アンテナ（送信手段、受信手段）１３、信号受信部（受信手段）１４、レンジ圧縮処理部（レンジ圧縮処理手段）１５、モデルデータ算出部（モデルデータ算出手段）１６、サンプルデータ選択部（相関値算出手段、選択手段）１７、相互相関演算部（相関値算出手段、相互相関演算手段）１８及び画像出力部（出力手段）１９を含んで構成されている。

合成開口処理システム１０のアンテナ１３以外の構成要素は、物理的には、当該アンテナ１３に接続されたコンピュータシステムとして構成されている。コンピュータシステムは、ＣＰＵ（Ｃｅｎｔｒａｌ Ｐｒｏｃｅｓｓｉｎｇ Ｕｎｉｔ）、主記憶装置であるＲＡＭ（Ｒａｎｄｏｍ Ａｃｃｅｓｓ Ｍｅｍｏｒｙ）およびＲＯＭ(Ｒｅａｄ Ｏｎｌｙ Ｍｅｍｏｒｙ)、ハードディスク等の補助記憶装置、入力デバイスであるキーボード及びマウス等の入力装置、ディスプレイ等の出力装置などを含む。図２に示す合成開口処理システム１０の計測条件入力部１１、信号送信部１２、信号受信部１４、レンジ圧縮処理部１５、モデルデータ算出部１６、サンプルデータ選択部１７、相互相関演算部１８及び画像出力部１９は、ＣＰＵ、ＲＡＭ等のハードウェア上に所定のコンピュータソフトウェアを読み込ませることにより、ＣＰＵの制御のもとで通信モジュール、入力装置、出力装置を動作させるとともに、ＲＡＭや補助記憶装置におけるデータの読み出し及び書き込みを行うことで実現される。以下、図２に示す機能ブロックに基づいて、各機能ブロックを説明する。

計測条件入力部１１は、各種計測条件を入力する部分である。計測条件は、例えばオペレータにより入力される形式でもよいし、記憶装置などに予め記憶されているデータを読み出す形式でもよい。計測条件とは、具体的には、プラットフォームの進行速度、チャープ波を出力する間隔、チャープ波の速度、合成開口処理を行うアジマス方向のデータ数、受信サンプリング周波数、送信チャープ率などである。これらの計測条件のデータは、信号送信部１２、信号受信部１４及びモデルデータ算出部１６に送信される。

信号送信部１２は、周波数が時間変化する周波数変調波（チャープ波）を生成して、複数の計測位置から照射領域に照射する部分である。信号送信部１２は、計測条件入力部１１から受信した各種計測条件のうち、例えばチャープ波の速度などの条件に従って、チャープ波を作成する。そして、プラットフォームの進行速度やチャープ波を出力する間隔などの条件に従って、アジマス方向の複数の位置でアンテナ１３から照射させる。

本実施形態では、図３に示すように、アジマス方向に沿った等間隔の５箇所の位置ｎ１〜ｎ５を１つのセットとして合成開口処理を行う。図３では、アジマス方向の５箇所の位置ｎ１〜ｎ５からレンジ方向（図３では右方向）へそれぞれチャープ波ｃ１〜ｃ５が照射されていることが示されている。図１を参照して説明したように、合成開口処理では、アジマス方向の複数の位置における各照射領域３ａ〜３ｃの共通部分５が検出対象となるものであるので、本実施形態の場合、図３に示す５箇所の位置ｎ１〜ｎ５を１セットとして合成開口処理を行うことにより、これらの位置ｎ１〜ｎ５の中心である位置ｎ３からのレンジ方向の領域Ａが検知対象となる。

信号受信部１４は、信号送信部１２により照射された複数のチャープ波による複数の反射波を受信する部分である。信号受信部１４は、計測条件入力部１１から受信した各種計測条件のうち例えば受信サンプリング周波数などの条件に従って、アンテナ１３で受信された反射波にフィルタ処理、増幅処理、デジタル処理などを施して、これらをレンジ圧縮処理部１５に送信する。

図４に、図３で例示したチャープ波ｃ１〜ｃ５各々に対する反射波ｒ１〜ｒ５の一例を示す。図４では、アジマス方向の位置ｎ３からのレンジ方向に物体（海底面などのものを含む）があった場合に、チャープ波ｃ１〜ｃ５がこの物体で反射されて、信号受信部１４により受信された反射波ｒ１〜ｒ５が示されている。信号受信部１４は、反射波を示すデータを所定の座標系上のデータとして取得する。具体的には、反射波ｒ１〜ｒ５は、レンジ方向及びアジマス方向の位置座標と、信号強度とを含むデータからなるデータ列となる。

データ列のレンジ方向の間隔は、計測条件入力部１１で設定された受信サンプリング周波数などの条件により決められる。各反射波ｒ１〜ｒ５のアジマス方向の位置座標は、チャープ波ｃ１〜ｃ５の照射位置ｎ１〜ｎ５に対応している。反射波ｒ１〜ｒ５のレンジ方向の位置座標は、チャープ波ｃ１〜ｃ５の送信から反射波ｒ１〜ｒ５の受信までの時間差に応じて、信号受信部１４により決められる。時間差が小さく計測位置から物体までの距離が短いほど照射位置に近づくこととなる。図４の例では、計測位置から物体までの距離が最短となる位置ｎ３における反射波ｒ３が最も照射位置寄りに配置され、位置ｎ２及びｎ４の反射波ｒ２及びｒ４が反射波ｒ３より照射位置から遠く配置され、位置ｎ１及びｎ５の反射波ｒ１及びｒ５がさらに照射位置から遠くなるよう配置されている。

レンジ圧縮処理部１５は、信号受信部１４により受信された複数の反射波ｒ１〜ｒ５のそれぞれにレンジ圧縮を行い、複数のシンク関数ｓ１〜ｓ５からなる受信データを生成する部分である。具体的には、レンジ圧縮処理部１５は、自己相関演算などの従来のパルス圧縮手法を用いて、反射波ｒ１〜ｒ５をレンジ方向に沿って圧縮してデータ数を少なくしたシンク関数状の受信データに変換する。図５に、図４で例示した反射波をレンジ圧縮した後の受信データの波形ｓ１〜ｓ５を示す。図５に示すように、受信データは、図４の各反射波ｒ１〜ｒ５の時間差を反映したシンク関数状の波形ｓ１〜ｓ５となる。これらのシンク関数状の波形ｓ１〜ｓ５も、反射波ｒ１〜ｒ５と同様に、レンジ方向及びアジマス方向の位置座標と、信号強度とを含むデータからなるデータ列である。生成された受信データは、サンプルデータ選択部１７に送信される。

モデルデータ算出部１６は、チャープ信号が照射領域内の複数の所定位置で反射した場合の理想的な受信波をレンジ圧縮した複数のシンク関数からなるモデルデータを、照射領域内の複数の所定位置ごとに算出する部分である。ここで、「照射領域内の複数の所定位置」とは、合成開口処理による検出対象の位置であり、本実施形態の場合、図６に示すように、合成開口処理を行う１セットのアジマス方向位置ｎ１〜ｎ５のうち中心の位置ｎ３からのレンジ方向の検出対象領域Ａを５個に区分した領域ａ〜ｅのことを指す。すなわち、モデルデータは、領域ａ〜ｅのそれぞれに応じて算出される。また、「理想的な受信波」とは、物体で反射した後にも何ら散乱及び屈折していないものと仮定されたチャープ波のことを指す。これらの計測位置や受信波は、計測条件入力部１１から受信した各種計測条件に従って設定される。

モデルデータの算出において、まず、各領域ａ〜ｅのそれぞれにおいて、当該領域のレンジ方向の位置座標に応じて求められる距離だけ計測位置ｎ３からレンジ方向に離れた場所に物体があると仮定して、位置ｎ１〜ｎ５からチャープ波ｃ１〜ｃ５が照射されてこの物体に反射してもとの位置に戻るまでの時間差が算出される。この時間差に基づいてレンジ方向の位置座標が決められ、反転したチャープ波ｃ１〜ｃ５が、この位置座標に反射波として配置される。これらの反射波は、具体的には、計測条件入力部１１から受信した各種計測条件のうち例えば受信サンプリング周波数などの条件に従って、信号受信部１４により取得される反射波ｒ１〜ｒ５と同様の処理が施されたデータ列である。データ列はレンジ方向及びアジマス方向の位置座標と、信号強度とを含むデータからなり、データ列のレンジ方向の間隔は信号受信部１４により取得される反射波ｒ１〜ｒ５のものと同一である。

次いで、これらの反射波がレンジ圧縮処理により５個のシンク関数状の波形に変換される。これらのシンク関数状の波形も、レンジ圧縮処理部１５により作成される波形ｓ１〜ｓ５と同様に、レンジ方向及びアジマス方向の位置座標と、信号強度とを含むデータからなるデータ列である。

続いて、これらの５個のシンク関数状の波形をなすデータ列のそれぞれから２点ずつ抽出され、合計１０点のデータからなるモデルデータが作成される。シンク関数は図５や図６に示したように単一のパルスピークをもっており、各シンク関数状のデータ列から、このパルスピークを挟みパルス幅と同じ間隔をとる２点が抽出される。具体的には、例えば、データ列の中から、信号強度の和が最大となる隣り合う２点が抽出される。なお、本実施形態では、シンク関数のパルスピークを挟んでパルス幅と同じ間隔をとる２点が得られるように、予め計測条件入力部１１において受信サンプリング周波数や送信チャープ率が調整されている。

このように、アジマス方向の位置ｎ３におけるレンジ方向の５つの領域ａ〜ｅのそれぞれに対応するモデルデータ（ａ）〜（ｅ）が算出される。図６に、モデルデータ（ａ）〜（ｅ）の一例を示す。図６に示すように、モデルデータ（ａ）〜（ｅ）は、対応する領域ａ〜ｅのレンジ方向の位置に応じて、レンジ方向の配置が異なる。算出されたモデルデータ（ａ）〜（ｅ）は、サンプルデータ選択部１７及び相互相関演算部１８に送信される。

サンプルデータ選択部１７は、モデルデータ算出部１６により算出された複数のモデルデータ（ａ）〜（ｅ）のそれぞれについて、このモデルデータを構成するシンク関数ごとの２点のデータの位置情報に応じて、受信データを構成するシンク関数ｓ１〜ｓ５の中からそれぞれ２点のデータをサンプルデータとして選択する部分である。サンプルデータ選択部１７は、レンジ圧縮処理部１５で算出された受信データの中から、各モデルデータ（ａ）〜（ｅ）のレンジ方向のデータ配置に応じて、１０点のデータをそれぞれ選択する。

図７に、モデルデータ（ａ）〜（ｅ）に応じて選択されるサンプルデータの一例を示す。図７において点線で示され、符号ａ〜ｅを付された領域は、アジマス方向の位置ｎ１〜ｎ５における、図６に示したモデルデータ（ａ）〜（ｅ）のレンジ方向の位置を示している。受信データの各シンク関数ｓ１〜ｓ５のそれぞれにおいて、領域ａ〜ｅに含まれる２点のデータが、それぞれモデルデータ（ａ）〜（ｅ）に対応して選択される。そして、アジマス方向の位置ｎ１〜ｎ５で選択されたデータが領域ａ〜ｅごとに１０点ずつ取りまとめられ、１つのモデルデータにつき１０点のデータがサンプルデータとして選択される。選択されたサンプルデータは、相互相関演算部１８に送信される。

相互相関演算部１８は、サンプルデータ選択部１７により選択された受信データのサンプルデータと、前記モデルデータ算出部１６により算出された複数のモデルデータのデータセットとの間でそれぞれ相互相関演算を行い、相関値データを算出する部分である。

具体的には、まず、複数のモデルデータ（ａ）〜（ｅ）から１つを選択し、選択したモデルデータの各データについて、例えば位置ｎ１の原点寄りのデータを１番目、他方のデータを2番目、位置ｎ５の原点寄りのデータを９番目、他方のデータを１０番目としてデータ列を作成する。同様に、サンプルデータも位置ｎ１の原点寄りのデータを１番目、他方のデータを2番目、位置ｎ５の原点寄りのデータを９番目、他方のデータを１０番目としてデータ列を作成する。そして、モデルデータ及びサンプルデータのデータ列を用いて相互相関演算を行う。すなわち、データ列の順番が同じデータ同士で、各データの信号強度を掛け合わせ、掛け合わせた値を総和演算して、相関値を算出する。このような処理を全てのモデルデータ（ａ）〜（ｅ）について行い、本実施形態では、図８に示すように図６の領域ａ〜ｅに対応する５個の相関値データを算出する。

図８に、図６に示したアジマス位置ｎ３のレンジ方向に区分された領域ａ〜ｅに関する相関値データの一例を示す。図８の符号ａ〜ｅは、それぞれが図６の領域ａ〜ｅにおける相関値であることを表している。これらの相関値ａ〜ｅは、モデルデータとサンプルデータとの相関の度合いを示し、両者の各データの信号強度のパターンが近いほど大きな値となる。ある１つのモデルデータとサンプルデータとの相関値が大きいということは、このモデルデータの対応する領域に実際に物体が存在する可能性が高いことを表す。図８の例では、領域ｃの相関値が最も大きく、このことは領域ｃに物体が存在する可能性が最も高いことを示している。

画像出力部１９は、相互相関演算部１８により算出された相関値データに基づき、ソーナー画像を出力する部分である。具体的には、これまで説明したようにアジマス方向の任意の位置において、その位置の前後２箇所までの受信データを利用して、当該位置のレンジ方向の相関値データを算出する処理を、アジマス方向に１点ずつ移動しながら繰り返して行い、アジマス方向の各位置の相関値データを取得する。そして、これらの相関値データをアジマス方向に沿って並べて、各相関値をアジマス方向及びレンジ方向に関して二次元的に配列する。これらの配列された相関値の大きさに応じて、色が割り当てられ、例えば白黒階調の場合、相関値が大きいほど白に近くなるように色が指定される。画像出力部１９は、このように指定された色の情報に基づいて、ソーナー画像を作成する。

図９に、従来のソーナー１によるソーナー画像（ａ）と、本実施形態の手法により出力されたソーナー画像（ｂ）の一例を示す。図９において、ソーナー画像（ａ）、（ｂ）のそれぞれの横軸はアジマス方向の位置座標を示し、縦軸はレンジ方向の位置座標を示す。図９（ｂ）では、相互相関演算部１９により算出された図８に示すような相関値データは縦方向に向いて配置されており、複数の相関値データが横軸のアジマス方向に沿って並べられて二次元的に配列されている。本実施形態の手法により出力されたソーナー画像（ｂ）は、従来のソーナー１によるソーナー画像（ａ）に比べて、シャープネスやコントラストが改善されており、ソーナー画像の分解能が向上している。

次に、図１０のフローチャートを用いて、本実施形態の合成開口処理システム１０の処理について説明する。

計測条件入力部１１により、各種計測条件を入力され（Ｓ１０１）、入力された計測条件のデータが信号送信部１２及びモデルデータ算出部１６に送信される。

計測条件に従って、信号送信部１２により、チャープ波が作成され、所定のタイミングでアンテナ１３から照射される（Ｓ１０２：送信ステップ）。同様に計測条件に従って、モデルデータ算出部１６により、チャープ信号が照射領域内の複数の所定位置で反射した場合のモデルデータが、照射領域内の複数の所定位置ごとに算出される（Ｓ１０３）。

次に、信号受信部１４により、照射されたチャープ波による複数の反射波が受信され（Ｓ１０４：受信ステップ）、レンジ圧縮処理部１５により、複数の反射波のそれぞれにレンジ圧縮が行われ、複数のシンク関数からなる受信データが生成される（Ｓ１０５：レンジ圧縮処理ステップ）。

そして、サンプルデータ選択部１７により、モデルデータの中から１つのモデルデータが選択され（Ｓ１０６）、この選択されたモデルデータに応じて、受信データを構成するシンク関数の中からそれぞれ２点のデータが選択され、受信データの中から合計で１０点がサンプルデータとして選択される（Ｓ１０７）。

次に、相互相関演算部１８により、ステップＳ１０６で選択されたモデルデータと、ステップＳ１０７で選択されたサンプルデータとの間で相互相関演算が行われ、相関値データが算出される（Ｓ１０８：相関値算出ステップ）。

続いて、サンプルデータ選択部１７により、ステップＳ１０３で算出された全てのモデルデータが既に選択済みか否かが確認される（Ｓ１０９）。まだ未選択のモデルデータがある場合には、ステップＳ１０６に戻り、ステップＳ１０６〜Ｓ１０８の処理を繰り返す。全てのモデルデータが選択済みである場合には、ステップＳ１１０に移行する。

ステップＳ１０９において全てのモデルデータが選択済みであると判定された場合には、画像出力部１９により、相関値データに基づくソーナー画像が出力される（Ｓ１１０：出力ステップ）。

なお、ステップＳ１０６〜Ｓ１０９の処理については、１つのモデルデータを選択して逐次相関値データを算出する代わりに、全てのモデルデータに関するサンプルデータを選択した上で、相互相関演算をまとめて行う形式でもよい。また、モデルデータ算出ステップ（Ｓ１０３）は、計測条件入力ステップ（Ｓ１０１）の後であり、モデルデータ選択ステップ（Ｓ１０６）の前であれば、任意の順番で実行してもよい。

以上説明したように、本実施形態に係る合成開口処理システム１０では、まずレンジ圧縮処理部１５がレンジ圧縮を行って、反射波に比べてデータ数が少ないシンク関数からなる受信データを用いて、相互相関演算部１８が相関値データを算出するので、演算負荷が低減できる。また、従来の２次元アルゴリズムと同様に、相互相関演算部１８が受信データとモデルデータとの相関値データを算出するので、従来の１次元アルゴリズムのレンジカーバチャ修正が不要となり、精度誤差が大きくなるのを回避できる。

また、本実施形態に係る合成開口処理システム１０では、モデルデータ算出部１６が、レンジ圧縮処理部１５により生成されるシンク関数状の波形からなる受信データと同様の波形を作成し、これらのシンク関数状の波形のパルスピークを挟みパルス幅と同じ間隔をとる２点を抽出して、モデルデータを作成する。図１１に示すように、この２点Ｐ１、Ｐ２は、「２点Ｐ１、Ｐ２がパルスピークを挟むこと」及び「２点の間隔Ｌ１がパルス幅（Ｌ／２）と同一であること」という条件を満たせば、点Ｐ１の信号強度Ｌ２と点Ｐ２の信号強度Ｌ３の和が常にパルスピークの信号強度Ｌ４とほぼ同一となる（Ｌ２＋Ｌ３≒Ｌ４）という特徴をもつ。受信データは、サンプリング周波数との兼ね合いにより、シンク関数状の波形のパルスピークの位置のデータがデータ列に含まれない状況が考えられる。相関関係を精度良く把握するために、従来はオーバーサンプリングや周波数領域での補間処理を行ってパルスピークの位置のデータを作成した上で相関演算を行っていた。本実施形態に係る合成開口処理システム１０は、パルスピークを挟みパルス幅と同じ間隔をとる２点を用いてサンプルデータ（受信データ）との相関演算を行えば、パルスピークの位置のデータを利用した場合と実質的に同等の相関関係の傾向を獲得することができる。

このように、受信データやモデルデータに含まれるシンク関数のパルスピークの位置を厳密に一致させなくても両者の相関関係を表すことができ、従来のパルスピーク位置のデータを得るための処理が不要となるので、演算負荷をさらに低減することができる。さらに、実際に使用するチャープ波の情報に応じてモデルデータが生成されるので、合成開口処理を確実に実施できる。

（第２実施形態）
次に、本発明に係る合成開口システムの第２の実施形態について説明する。図１２に本発明の第２の実施形態に係る合成開口処理システムの構成を示す。図１２に示すように、本実施形態の合成開口処理システム２０は、図２の合成開口処理システム１０とは異なり、レンジ圧縮処理部１５により生成された受信データからサンプルデータを抽出せずに、受信データの全てのデータを用いてモデルデータとの相互相関演算を行う。

合成開口処理システム２０における、計測条件入力部１１、信号送信部１２、アンテナ１３、信号受信部１４、レンジ圧縮処理部１５及び画像出力部１９の各機能ブロックは、図２に示した合成開口処理システム１０のものと同一である。

モデルデータ算出部２６は、図６に示したように、アジマス方向の位置ｎ３からのレンジ方向の検出対象領域Ａを５個に区分した領域ａ〜ｅのそれぞれにおいて、５個のシンク関数状の波形を作成するところまでは、図２のモデルデータ算出部１６と同様の処理を行う。モデルデータ算出部２６は、その後シンク関数状の波形のそれぞれからパルスピークを挟んだ２点のデータを選択せずに、５個のシンク関数状の波形をなすデータ列を取りまとめて、１つのモデルデータとして算出する。したがって、モデルデータ算出部２６により算出されるモデルデータのデータ数は、図２のモデルデータ算出部１６によるもの（１０点）よりも多くなる。算出されたモデルデータは、相互相関演算部２８に送信される。

相互相関演算部２８は、レンジ圧縮処理部１５により生成された受信データと、モデルデータ算出部２６により算出された複数のモデルデータとの間でそれぞれ相互相関演算を行い、相関値データを算出する。

本実施形態に係る合成開口処理システム２０によれば、信号送信部１２がチャープ波を生成して、複数の計測位置から照射領域に照射すると、信号受信部１４が、照射された複数のチャープ波による複数の反射波を受信する。レンジ圧縮処理部１５が、受信された複数の反射波のそれぞれにレンジ圧縮を行い、複数のシンク関数からなる受信データを生成する。また、モデルデータ算出部２６は、チャープ波が照射領域内の複数の所定位置で反射した場合の理想的な受信波をレンジ圧縮した複数のシンク関数からなる複数のモデルデータを算出する。相互相関演算部２８は、レンジ圧縮処理部１５により生成された受信データと、複数のモデルデータのデータセットとの間でそれぞれ相互相関演算を行い、相関値データを算出する。画像出力部１９は、相互相関演算部１８により算出された相関値データを出力する。

これらの構成により、まずレンジ圧縮を行って、反射波に比べてデータ数が少ないシンク関数からなる受信データを用いて相関値データを算出するので、演算負荷を低減することができる。また、従来の２次元アルゴリズムと同様に、受信データとモデルデータとの相関値データを算出するので、従来の１次元アルゴリズムのレンジカーバチャ修正が不要となり、精度誤差が大きくなるのを回避できる。また、実際に使用するチャープ波の情報に応じてモデルデータが生成されるので、合成開口処理を確実に実施できる。

（第３実施形態）
次に、本発明に係る合成開口システムの第３の実施形態について説明する。図１３に本発明の第３の実施形態に係る合成開口処理システムの構成を示す。図１３に示すように、本実施形態の合成開口処理システム３０は、図２の合成開口処理システム１０とは異なり、サンプルデータ選択部１７により選択されたサンプルデータに、モデルデータとの相互相関演算を行わず、単なる総和演算を行って相関値データを算出する。

合成開口処理システム３０における、計測条件入力部１１、信号送信部１２、アンテナ１３、信号受信部１４、レンジ圧縮処理部１５、モデルデータ算出部１６及び画像出力部１９の各機能ブロックは、図２に示した合成開口処理システム１０のものと同一である。

総和演算部３８は、サンプルデータ選択部１７により選択された受信データのサンプルデータの総和演算を行い、相関値データを算出する。具体的には、サンプルデータの各データの信号強度を全て足し合わせて、相関値データとして算出する。

本実施形態に係る合成開口処理システム３０によれば、信号送信部１２がチャープ波を生成して、複数の計測位置から照射領域に照射すると、信号受信部１４が、照射された複数のチャープ波による複数の反射波を受信する。レンジ圧縮処理部１５が、受信された複数の反射波のそれぞれにレンジ圧縮を行い、複数のシンク関数からなる受信データを生成する。また、モデルデータ算出部１６は、チャープ波が照射領域内の複数の所定位置で反射した場合の理想的な受信波をレンジ圧縮した複数のシンク関数からなる複数のモデルデータを算出する。また、モデルデータ算出部１６は、モデルデータを構成するシンク関数ごとに、ピークを挟んだ所定の間隔の２点のデータを抽出し、これらの抽出されたデータセットをモデルデータとして算出する。そして、サンプルデータ選択部１７は、モデルデータ算出部１６により算出された複数のモデルデータのそれぞれについて、このモデルデータを構成するシンク関数ごとの２点のデータの位置情報に応じて、受信データを構成するシンク関数の中からそれぞれ２点のデータをサンプルデータとして選択する。総和演算部３８は、サンプルデータ選択部１７により選択された受信データのサンプルデータの総和演算を行い、相関値データを算出する。画像出力部１９は、相互相関演算部１８により算出された相関値データを出力する。

これらの構成により、まずレンジ圧縮を行って、反射波に比べてデータ数が少ないシンク関数からなる受信データを用いて相関値データを算出するので、演算負荷を低減することができる。また、従来の２次元アルゴリズムと同様に、受信データとモデルデータとの相関値データを算出するので、従来の１次元アルゴリズムのレンジカーバチャ修正が不要となり、精度誤差が大きくなるのを回避できる。また、受信データの中からサンプルデータを抽出して相互相関演算が行われるので、相関値データを算出するためのデータ数がさらに少なくなり、演算負荷をさらに低減することができる。また、サンプルデータを抽出した後にこれらのデータの総和演算を行うだけで相関値データが算出されるので、演算負荷をさらに低減できる。さらに、実際に使用するチャープ波の情報に応じてモデルデータが生成されるので、合成開口処理を確実に実施できる。

なお、本発明は上記の実施形態に限定されることなく種々の変形が可能である。例えば、合成開口処理を行うアジマス方向の位置の数は、上記実施形態で示した５箇所以外の数でもよい。また、上記の実施形態ではモデルデータ算出部１６、２６が、計測条件入力部１１により入力された各種計測条件に応じてモデルデータを作成していたが、予め様々な計測条件に応じたモデルデータを記憶装置などに用意しておき、計測条件に応じて適宜選択する形式でもよい。

また、上記の実施形態では、合成開口ソーナーについて説明したが、本発明は、合成開口レーダーにも適用可能である。さらに、本発明は、合成開口処理を行わない従来のパルス圧縮処理にも適用可能である。

Claims (5)

1. 周波数が時間変化する周波数変調波を生成して、複数の計測位置から照射領域に照射する送信手段と、
   前記送信手段により照射された複数の周波数変調波による複数の反射波を受信する受信手段と、
   前記受信手段により受信された複数の反射波のそれぞれにレンジ圧縮を行い、複数のシンク関数からなる受信データを生成するレンジ圧縮処理手段と、
   前記周波数変調波が前記照射領域内の複数の所定位置で反射した場合の理想的な受信波をレンジ圧縮した複数のシンク関数からなる複数のモデルデータと、前記レンジ圧縮処理手段により生成された前記受信データとに基づいて、前記モデルデータのそれぞれと前記受信データとの相関度合いを示す相関値データを算出する相関値算出手段と、
   前記相関値算出手段により算出された相関値データを出力する出力手段と、
   を備える合成開口処理システム。
2. 前記送信手段により生成された前記周波数変調波に応じて、前記相関値算出手段で用いられる前記モデルデータを前記照射領域内の複数の所定位置ごとに算出するモデルデータ算出手段をさらに備える、請求項１に記載の合成開口処理システム。
3. 前記モデルデータ算出手段が、前記モデルデータを構成するシンク関数ごとに、ピークを挟んだ所定の間隔の２点のデータを抽出し、これらの抽出されたデータセットを前記モデルデータとして算出し、
   前記相関値算出手段が、
   前記モデルデータ算出手段により算出された複数のモデルデータのそれぞれについて、該モデルデータを構成するシンク関数ごとの２点のデータの位置情報に応じて、前記受信データを構成するシンク関数の中からそれぞれ２点のデータをサンプルデータとして選択する選択手段と、
   前記選択手段により選択された前記受信データのサンプルデータと、前記複数のモデルデータのデータセットとの間でそれぞれ相互相関演算を行い、相関値データを算出する相互相関演算手段と、
   をさらに備える、請求項１または２に記載の合成開口処理システム。
4. 前記モデルデータ算出手段が、前記モデルデータを構成するシンク関数ごとに、ピークを挟んだ所定の間隔の２点のデータを抽出し、これらの抽出されたデータセットを前記モデルデータとして算出し、
   前記相関値算出手段が、
   前記モデルデータ算出手段により算出された複数のモデルデータのそれぞれについて、該モデルデータを構成するシンク関数ごとの２点のデータの位置情報に応じて、前記受信データを構成するシンク関数の中からそれぞれ２点のデータをサンプルデータとして選択する選択手段と、
   前記選択手段により選択された前記受信データのサンプルデータの総和演算を行い、相関値データを算出する総和演算手段と、
   をさらに備える、請求項１または２記載の合成開口処理システム。
5. 周波数が時間変化する周波数変調波を生成して、複数の計測位置から照射領域に照射する送信ステップと、
   前記送信ステップにおいて照射された複数の周波数変調波による複数の反射波を受信する受信ステップと、
   前記受信ステップにおいて受信された複数の反射波のそれぞれにレンジ圧縮を行い、複数のシンク関数からなる受信データを生成するレンジ圧縮処理ステップと、
   前記周波数変調波が前記照射領域内の複数の所定位置で反射した場合の理想的な受信波をレンジ圧縮した複数のシンク関数からなる複数のモデルデータと、前記レンジ圧縮処理ステップにおいて生成された前記受信データとに基づいて、前記モデルデータのそれぞれと前記受信データとの相関度合いを示す相関値データを算出する相関値算出ステップと、
   前記相関値算出ステップにおいて算出された相関値データを出力する出力ステップと、
   を備える合成開口処理方法。