JPH10268042A - 補償回路、補償方法およびレーダ装置 - Google Patents

補償回路、補償方法およびレーダ装置

Info

Publication number
JPH10268042A
JPH10268042A JP9075073A JP7507397A JPH10268042A JP H10268042 A JPH10268042 A JP H10268042A JP 9075073 A JP9075073 A JP 9075073A JP 7507397 A JP7507397 A JP 7507397A JP H10268042 A JPH10268042 A JP H10268042A
Authority
JP
Japan
Prior art keywords
range
target
original image
correction
image
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
JP9075073A
Other languages
English (en)
Other versions
JP3642149B2 (ja
Inventor
Kazuhiko Yamamoto
山本  和彦
Masafumi Iwamoto
雅史 岩本
Takahiko Fujisaka
貴彦 藤坂
Ryusaburo Usui
隆三郎 臼井
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Mitsubishi Electric Corp
Original Assignee
Mitsubishi Electric Corp
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Mitsubishi Electric Corp filed Critical Mitsubishi Electric Corp
Priority to JP07507397A priority Critical patent/JP3642149B2/ja
Publication of JPH10268042A publication Critical patent/JPH10268042A/ja
Application granted granted Critical
Publication of JP3642149B2 publication Critical patent/JP3642149B2/ja
Anticipated expiration legal-status Critical
Expired - Fee Related legal-status Critical Current

Links

Landscapes

  • Radar Systems Or Details Thereof (AREA)

Abstract

(57)【要約】 【課題】 目標の移動に伴い目標のレンジ又は反射波の
ドップラー周波数が変化する場合であっても追尾誤差を
押えてレンジ追尾又はドップラー追尾を行なうこと。 【解決手段】このレーダ装置は、目標に高周波パルスを
照射し、受信したエコーをパルス圧縮して距離方向に高
分解能なデータ列(レンジプロフィール)を得る処理を
複数回繰り返し、結果として得られるレンジプロフィー
ルの距離追尾を行い、その結果を用いて同じ反射点から
の受信信号を同じレンジに並べるレンジ補償を行う。レ
ンジ追尾手段501では、レンジプロフィールのヒスト
リをヒット−レンジの二軸を持つ二次元画像とみなし、
これをシフト量推定回路115に入力する。シフト量推
定回路115では、入力された画像内に存在する同じ傾
きを有する複数の線分の傾きを推定して、その結果を出
力する。レンジ追尾手段501では、シフト量推定回路
115より出力された傾きを目標までの距離の時間変化
として出力する。

Description

【発明の詳細な説明】
【0001】
【発明の属する技術分野】この発明は、目標とレーダ装
置との間の相対位置関係の変化に基づく反射信号の変
化、及び送信周波数の変化に基づく反射信号の変化を利
用することにより、高い分解能を得ることのできるレー
ダ装置に関するものである。
【0002】
【従来の技術】従来のこの種のレーダ装置は、一般に高
分解能レーダ装置と呼ばれ、例えば、Donald R. Wehne
r,■■High Resolution Radar■■,Artech House,INC.1
987,pp273-339に記載されているISAR(Inverse Syntheti
c Aperture Radar)、及び、特開平7−92257号に
記載のものがある。
【0003】図17は、上記文献に従って構成したレー
ダ装置のブロック構成図である。図において、1は送信
機、2は送受切換器、3は送受信アンテナ、4は受信
機、5はレンジ圧縮手段、6はレンジ補償回路、7は位
相補償回路、8はクロスレンジ圧縮手段、9はモニタ・
テレビ(以下、モニタTVと呼ぶ)である。
【0004】図18は、図17のレーダ装置のブロック
構成図におけるレンジ補償回路6の内容を記した構成図
である。図において、501はレンジ追尾手段、502
は振幅最大レンジビン検出回路、503は平滑化回路、
504はレンジ補償量算出手段、505レンジ補償手段
である。6は図17のそれと同一である。
【0005】図19は、図17のレーダ装置のブロック
構成図における位相補償回路7の内容を記した構成図で
ある。図において、506は注目レンジビン決定手段、
507は区分周波数分析手段、508はドップラー追尾
手段、509は振幅最大周波数検出回路、510は位相
補償量算出手段、511は位相補償手段である。7は図
17のそれと、503は図18のそれと、それぞれ同一
である。
【0006】図20は、回転運動を行なう目標を観測す
るジオメトリである。図において、520は目標、52
1はレーダ装置である。図21は、図20のジオメトリ
で観測した結果得られたISAR画像の一例である。図
において、522は目標のISAR画像である。
【0007】図22は、並進運動を行なう目標を観測す
るジオメトリである。図において、520、521は、
図20のそれと、それぞれ同一である。図23は、レン
ジプロフィールのヒストリの最大振幅検出結果の一例で
ある。図24は、レンジ補償処理を施した後のレンジプ
ロフィールにおける最大振幅検出結果のヒストリの一例
である。
【0008】図25は、区分周波数分析手段507の処
理内容を説明する図である。図において、530は、注
目レンジビンの受信信号列、531は、区分周波数分布
のヒストリである。図26は、区分周波数分布のヒスト
リの最大振幅検出結果の一例である。
【0009】次に、図17、18、19、20、21、
22、23、24、25、26を用いて、上記従来のレ
ーダ装置の動作原理について説明する。まず、ISAR
の画像再生の原理について説明する。図20に示す通
り、x−y平面の原点に設置されたレーダ装置521
で、レンジr0の点oを通り紙面に垂直な軸を中心に角
速度ωで半時計周りで回転する目標520を観測するジ
オメトリを考える。まず、送信機1では、時間とともに
周波数が変化する信号(チャープ信号)に変調された高
周波パルスを発生し、送受切換器2を介して送受信アン
テナ3に供給する。送受信アンテナ3は、供給された高
周波パルスを送信する。
【0010】送信した高周波パルスは目標520で反射
され、この反射された信号(エコー)は、送受信アンテ
ナ3に入り、送受切換器2を介して受信機4で復調され
る。受信機4で復調された信号は、送信信号の瞬時周波
数に対しレーダ装置521と目標520の間の電波伝搬
の往復に要する時間分遅延したものであるから、レンジ
圧縮手段5において、送信信号s(t)を用いて受信信
号r(t)にマッチドフィルターをかけること、すなわ
ち、式(1)に示すように、送信信号s(t)の共役信
号s*(t)と受信信号r(t)とのコンボリューショ
ンを求めることにより、遅延に相当した時間にインパル
スv(t)(以下では、レンジプロフィールと呼ぶ。)
を得ることができる。このことにより、レンジ分解能が
向上する。
【0011】
【数1】
【0012】レーダが送信を繰り返すごとに上記レンジ
圧縮された信号が得られるから、レーダの送信ごと、即
ちヒットごとにレンジ圧縮された信号をまとめることに
より、レンジrとヒットhを軸とする二次元複素信号v
(h、r)(以下では、レンジプロフィールのヒストリ
と呼ぶ。)が得られる。尚、レンジrと時間tの間に
は、r=(Δr/Δt)・tなる関係がある。ここで、
Δtはサンプリング間隔(=1/B、Bは送信帯域)、
Δrはレンジ分解能(=C/2B、Cは光速)である。
【0013】目標520が図22に示すような運動、す
なわち並進運動を行なう場合には、レーダ装置521か
ら目標520までの距離の変化の影響を補償することに
より、並進運動を行なう目標520を、等価的に、図2
0に示すような回転運動を行なう目標520とみなすこ
とができる。この点については、後で述べる。レーダ5
21と目標520の間の距離の変化を補償するレンジ補
償回路6、位相補償回路7の動作は後述することにし、
以下では、目標520が回転運動を行なうもの、もしく
は、なんらかの方法で、レーダ装置521から目標52
0までの距離の変化の影響を補償されたものとする。
【0014】クロスレンジ圧縮手段8では、レンジプロ
フィールのヒストリv(h、r)を、式(2)に従っ
て、各レンジごとに、ヒット方向にフーリエ変換するこ
とにより、レンジ、クロスレンジの両方について圧縮さ
れた複素信号u(c、r)を得る。ここに、cはクロス
レンジ(方位)方向を示す。
【0015】
【数2】
【0016】上式で、hnumはヒット数、rnumは
レンジビン数である。この処理は、クロスレンジ方向の
分解能を改善する効果がある。以下、この原理を説明す
る。
【0017】図20に示す運動を行なう目標上のある部
位(例えば点a)で反射した信号のドップラー周波数f
dは次の式で表される。
【0018】
【数3】
【0019】ここに、λは送信波長、Lは反射を生じた
部位の回転半径、θは観測の基準となるLOS(Lin
e Of Sight)を基準としたときの目標物の角
度を示している。目標520が図22に示すような並進
運動を行なう場合には、上述の角度θはt=t0〜t2
に変わるにつれて変化することになる。このことから、
レーダ装置521から目標520までの距離の変化の影
響を補償すれば、並進運動を行なう目標520を、等価
的に角度θが変化する回転運動を行なう目標520とみ
なすことができることがわかる。また、式(3)によ
り、同じ角度θ上の点では、回転軸からの距離Lに比例
して、目標上のそれぞれの部位からの反射信号のドップ
ラー周波数が変わるといつ特性を有する。
【0020】レンジ圧縮手段5とクロスレンジ圧縮手段
8により、レンジ方向、クロスレンジ方向の両方につい
て分解能が向上した複素信号u(c、r)は、その絶対
値が目標のレーダ反射断面積に対応するから、モニタT
V9上のレンジr、クロスレンジcの二次元平面にu
(c、r)の絶対値またはその二乗に応じた輝度で表示
を行なうことにより、図21に示すような、レンジとク
ロスレンジの両方について高分解能化された目標のIS
AR画像522を表示することができる。ここで、画像
522上で、例えばa点は、目標上で、レーダに近い位
置にあるので、レンジが小さく、かつ、回転運動により
レーダから遠ざかる運動をしているのでドップラー周波
数が小さくなっている。
【0021】次に、上記説明で省略した、目標の並進運
動の影響を補償する処理について説明する。観測時間内
の、t0、t1、t2という時間において、並進運動に
より、図22のように位置が変化する目標520のIS
AR画像を生成する場合、各ヒットごとに得られたレン
ジプロフィールのヒストリv(h,r)をそのままヒッ
ト方向にフーリエ変換するだけでは、目標上の各点(例
えば点a)が観測時間中にレンジ方向に移動するため、
レンジ、クロスレンジ方向にきちんと圧縮されず、結果
として画像にぼけが生じてしまうのは、式(2)におい
て、各レンジごとにフーリエ変換を行なうというその処
理内容より明らかである。従って、ぼけのない鮮明な画
像を得るためには、目標上の各々の点を観測時間中、同
一レンジビン内に固定するための補償処理を必要とす
る。この処理をレンジ補償処理と呼ぶ。
【0022】レンジ補償処理を行なうレンジ補償処理回
路6の処理内容について説明する。目標が図22に示す
運動を行う場合について考える。ここで、電波の反射を
する点は、図中a、b、cの三点のみとし、このうち、
b点とc点は常に同じレンジビンにあるものとする。ま
ず、レンジ追尾手段501によりレンジ追尾を行なう。
レンジ追尾手段501の中の振幅最大レンジビン検出回
路502では、各ヒットごとに、レンジプロフィールの
振幅が最大となるレンジビンを検出する。その結果の例
を図23に示す。同図で横軸はヒット、縦軸はレンジで
あり、図中太実線で示した部分が、各ヒットのレンジプ
ロフィールで振幅が最大となるレンジビンを示している
ものとする。a、b、cは同一目標上の点であり、実際
は、図中の点線に示されるように、同じ変化率でレンジ
が変化しているはずであるが、見込み角の変化に従う各
点のレーダ断面積の変化や、同一レンジビン内に複数の
反射点が存在する場合の干渉などの影響で、観測時間中
に各点の存在するレンジビンの振幅が変動するため、振
幅最大レンジビンの位置の変化に不連続な部分が発生す
る。
【0023】この振幅最大レンジビンの位置の時間変化
に対して平滑化回路503では、例えば最小二乗法など
を用いて平滑化する。具体的には、平滑化回路503
は、図23における一点鎖線で示すような平滑化直線を
得る。そして、得られた平滑化直線の傾きを示すシフト
量sを算出する。そして、得られたシフト量をレンジ補
償算出手段504に対して通知する。レンジ補償量算出
手段504では、レンジ追尾手段501の平滑化回路5
03で得られたシフト量から、各ヒットにおけるレンジ
補償量sf(h)を式(4)により得る。
【0024】
【数4】
【0025】レンジ補償手段505では、式(4)で得
られたレンジ補償量sf(h)を用いて、各ヒットhに
おけるレンジプロフィールのヒストリv(h、r)をレ
ンジ方向に補償し、レンジ補償後のレンジプロフィール
のヒストリv2(h、r)を得る。レンジ補償後のレン
ジプロフィールのヒストリv2(h、r)で、各ヒット
ごとに、レンジプロフィールの振幅が最大となるレンジ
ビンを検出した結果は、図24に示すようにそれぞれの
点の反射信号が同一レンジビンに並ぶ。
【0026】上記レンジ補償処理により、観測時間中の
各点のレンジビンを超えた距離変化については除去する
ことができたが、レンジビン内の距離変化については除
去できていない。一般に、目標が加速運動、旋回運動を
する時は勿論の事、等速直線運動を行なう場合でも、進
行方向がLOS軸に沿った方向で無い限り、その距離変
化は、線形な成分に加えて、非線形な加速度成分も含
む。このうちの加速度成分の影響で、各点よりの反射信
号のドップラー周波数(クロスレンジ)に広がりが生じ
るため、結果として生成した画像がクロスレンジ方向に
ぼけてしまう。上記加速度成分を除去するための補償処
理が位相補償処理である。
【0027】以下では、この位相補償処理を行なう位相
補償回路7の処理内容について説明する。並進運動に伴
う上記加速度成分は、すべてのレンジビンに対してほぼ
等しく加わるため、ある一つのレンジビンに着目して、
そのレンジビンに加わる加速度成分を推定し、その推定
結果を用いて、すべてのレンジビンの位相補償を行な
う。
【0028】注目レンジビン決定手段506では、レン
ジ補償後のレンジプロフィールv2(h、r)の各レン
ジrにおける平均電力を算出し、その値を最大とするレ
ンジビンを注目レンジビンとして、そのレンジビンの受
信信号列w(h)を出力する。例えば図22のジオメト
リの例では、点b、cを含むレンジビンが注目レンジビ
ンとして選択されたものとする。
【0029】区分周波数分析手段507では、式(5)
に従い、図25に示した注目レンジビンの受信信号列w
(h)530を長さΔhで区分フーリエ変換して周波数
分布のヒストリfs(h’,f)531を求める。
【0030】
【数5】
【0031】ドップラー追尾手段508では、得られた
周波数分布のヒストリfs(h’,f)の追尾を行な
う。まず、振幅最大周波数検出回路509では、各ヒッ
トh’ごとに、周波数分布の振幅が最大となるドップラ
ービンを検出する。その結果の例を図26に示す。同図
で横軸はヒット、縦軸はドップラービンであり、図中太
実線で示した部分が、各ヒットの周波数分布で振幅が最
大となるドップラーを示しているものとする。見込み角
の変化に従うレーダ断面積の変化のため、振幅最大周波
数の位置が変動すること、および、周波数の折り返しの
影響で、その位置の変化に不連続な部分が発生する。
【0032】これに対し、レンジ追尾手段501と同
様、平滑化回路503により平滑化を行なうことにより
図26に一点鎖線で示すような平滑化直線を得る。そし
て、得られた平滑化直線の傾きを示すシフト量sを位相
補償量算出手段510に対して通知する。位相補償量算
出手段510では、式(6)に従い、位相補償量ph
(h)を計算する。
【0033】
【数6】
【0034】位相補償手段511では、位相補償量算出
手段510で得られた位相補償量ph(h)を用いて、
式(7)により、レンジ補償後のレンジプロフィールの
ヒストリv2(h、r)の位相補償を行ない、最終的な
レンジプロフィールのヒストリvL(h、r)を得る。
【0035】
【数7】
【0036】以上の処理を経る事により、並進運動を行
なう目標に関して、目標上の各点のレンジビンを超える
移動、位相の二次の変動を補償することができるため、
並進運動を伴わず、回転運動のみを行なう目標と同様に
鮮明な高分解能画像が得られる。
【0037】
【発明が解決しようとする課題】上記のような従来のレ
ーダ装置では、レーダ断面積の角度特性や干渉の影響の
ために各ヒットにおけるレンジプロフィールの最大点が
急に移動するという現象が起こった場合には、単純に平
滑を行うと図23の一点鎖線に示すように実際の目標の
距離の時間変化を表す直線とは傾きの異なった平滑化直
線を得ることになる。したがって、得られた平滑化直線
によって得られるシフト量を基にレンジ補償を行った場
合には、レンジ補償量に誤差が生じて、再生画像がぼけ
るという問題が有った。
【0038】また、区分周波数分析の結果得られた周波
数分布のヒストリにおいて、各ヒットにおける振幅最大
周波数がレーダ断面積の角度特性や周波数の折り返しの
影響で急に移動するという現象が起こった場合には、単
純に平滑を行うと図26の一点鎖線に示すように実際の
目標上の各点のドップラー周波数の時間変化を表す直線
とは傾きの異なった平滑化直線を得ることになる。した
がって、得られた平滑化直線によって得られるシフト量
を基に位相補償を行った場合には、位相補償量に誤差が
生じて、再生画像がぼけるという問題が有った。
【0039】また、雑音が多い環境では、レンジ追尾に
おける振幅最大レンジビン、ドップラー追尾における、
振幅最大ドップラービンとして、雑音により振幅が大き
くなったセルを選択してしまう場合がある。この場合に
は、雑音部分のセルの影響によって平滑化直線が実際の
目標の距離の時間変化を表す直線又は実際の目標上の各
点のドップラー周波数の時間変化を表す直線とは異なっ
た傾きをもつことになる。このような平滑化直線を用い
てレンジ補償又は位相補償を行うと、当該各補償量に誤
差が生じて、再生画像がぼけるという問題点があった。
【0040】また、目標の加速度運動成分が大きく、レ
ンジプロフィールのヒストリ、もしくは、区分周波数分
布のヒストリにおける振幅最大点のビンの移動が、時間
に対して非線形になる場合に、非線形変動成分まで含め
て追尾を行なうことができないため、この影響で再生画
像がぼけるという問題点が有った。
【0041】また、ドップラー追尾において、区分周波
数分析の点数を少なくすると、ドップラー周波数の分解
能が劣化するため、これに基づく追尾誤差が発生する、
区分周波数分析の点数を多くすると、ドップラーの分解
能は向上するものの、時間の分解能が劣化するため、信
号のチャープ成分によりドップラーが広がって、これに
基づく追尾誤差が発生するという理由で、追尾誤差が発
生し、結果として、再生画像がぼけるという問題点が有
った。
【0042】この発明は、かかる問題点を解決するため
になされたもので、目標の移動に伴い目標のレンジが変
化する場合であってもレンジ追尾誤差を押さえてレンジ
追尾を行なうこと、又は目標の移動に伴い反射波のドッ
プラー周波数が変化する場合であっても追尾誤差を抑え
てドップラー追尾を行うことを共通の目的としている。
【0043】
【課題を解決するための手段】この発明に係る補償回路
は、移動する目標に対して電波を送信し、前記目標から
の反射波を受信して前記目標の画像を得るため、前記目
標の移動に伴う前記目標のレンジの変化を補償する補償
回路であって、前記反射波の受信信号列を基に、移動に
伴う前記目標のレンジの変化を示す原画像を作成する原
画像作成手段と、複数の補正量によって前記原画像を補
正することにより複数の補正画像を作成する補正手段
と、前記補正手段によって作成された複数の補正画像を
評価するための評価値を算出する評価値算出手段と、前
記評価値算出手段によって前記複数の補正画像に対応し
て算出された評価値に基づき前記原画像を補正するのに
適した補正量を決定し、この決定した補正量により前記
目標の移動に伴う前記目標のレンジの変化量を推定する
変化量推定手段と、前記変化量推定手段によって推定し
た前記目標のレンジの変化量に応じて前記受信信号列の
レンジ補償を行うレンジ補償手段とを有するものであ
る。
【0044】この発明に係る補償回路は、移動する目標
に対して電波を送信し、前記目標からの反射波を受信し
て前記目標の画像を得るため、前記目標の移動に伴う前
記反射波の位相の変化を補償する補償回路であって、前
記反射波の受信信号列を基に、前記目標の移動に伴う前
記反射波のドップラー周波数の変化を示す原画像を作成
する原画像作成手段と、前記原画像を複数の補正量によ
って補正することにより複数の補正画像を作成する補正
手段と、前記補正手段によって作成された複数の補正画
像を評価するための評価値を算出する評価値算出手段
と、前記評価値算出手段によって前記複数の補正画像に
対応して算出された評価値に基づき前記原画像を補正す
るのに適した補正量を決定し、この決定した補正量によ
り前記目標の移動に伴う前記反射波のドップラー周波数
の変化量を推定する変化量推定手段と、前記変化量推定
手段によって推定された前記反射波のドップラー周波数
の変化量に応じて、前記受信信号列の位相補償を行う位
相補償手段とを有するものである。
【0045】また、前記評価値算出手段が、前記補正画
像の各レンジ毎に、同一のレンジに存在している複数の
画素の値を積分し、当該積分値の分布を示す指標を評価
値として算出するものである。
【0046】さらに、前記評価値算出手段が、前記補正
画像の各ドップラー周波数毎に、同一のドップラー周波
数に対して存在している複数の画素の値を積分し、当該
積分値の分布を示す指標を評価値として算出するもので
ある。
【0047】この発明に係る補償回路は、移動する目標
に対して電波を送信し、前記目標からの反射波を受信し
て前記目標の画像を得るため、前記目標の移動に伴う前
記目標のレンジの変化を補償する補償回路であって、前
記反射波の受信信号列を基に、移動に伴う前記目標のレ
ンジの変化を示す原画像を作成する原画像作成手段と、
前記原画像作成手段によって作成された原画像を二値化
することによりニ値化画像を作成する画像二値化手段
と、前記二値化画像を複数の補正量によって補正するこ
とにより複数の補正ニ値化画像を作成する補正手段と、
前記補正手段によって作成された複数の補正二値化画像
を射影することにより複数の射影画像を作成し、この複
数の射影画像に対する評価値を算出する評価値算出手段
と、前記評価値算出手段によって前記複数の射影画像に
対応して算出された評価値に基づき前記原画像を補正す
るのに適した補正量を決定し、この決定した補正量によ
り前記目標の移動に伴う前記目標のレンジの変化量を推
定する変化量推定手段と、前記変化量推定手段によって
推定した前記目標のレンジの変化量に応じて、前記受信
信号列のレンジ補償を行うレンジ補償手段とを有するも
のである。
【0048】この発明に係る補償回路は、移動する目標
に対して電波を送信し、前記目標からの反射波を受信し
て前記目標の画像を得るため、前記目標の移動に伴う前
記反射波の位相の変化を補償する補償回路であって、前
記反射波の受信信号列を基に、移動に伴う前記反射波の
ドップラー周波数の変化を示す原画像を作成する原画像
作成手段と、前記原画像作成手段によって作成された原
画像を二値化することによりニ値化画像を作成する画像
二値化手段と、前記二値化画像を複数の補正量によって
補正することにより複数の補正ニ値化画像を作成する補
正手段と、前記複数の補正二値化画像を射影することに
より複数の射影画像を作成し、この複数の射影画像に対
する評価値を算出する評価値算出手段と、前記評価値算
出手段によって前記複数の射影画像に対応して算出され
た評価値に基づいて前記原画像を補正するのに適した補
正量を決定し、この決定した補正量により前記目標の移
動に伴う前記目標のレンジの変化量を推定する変化量推
定手段と、前記変化量推定手段によって推定した前記反
射波のドップラー周波数の変化量に応じて、前記受信信
号列の位相補償を行う位相補償手段とを有するものであ
る。
【0049】また、前記評価値算出手段が、前記射影画
像内の特定の値を有する画素の画素数を評価値として算
出するものである。
【0050】さらに、前記補正手段が、前記複数の補正
量に応じて前記原画像を線形補正するものである。
【0051】さらにまた、前記補正手段が、前記複数の
補正量に応じて前記原画像を非線形補正するものであ
る。
【0052】この発明に係る補償方法は、移動する目標
に対して電波を送信し、前記目標からの反射波を受信し
て前記目標の画像を得るため、前記目標の移動に伴う前
記目標のレンジの変化を補償する補償方法であって、前
記反射波の受信信号列を基に、移動に伴う前記目標のレ
ンジの変化を示す原画像を作成する原画像作成ステップ
と、前記原画像を複数の補正量によって補正することに
より複数の補正画像を作成する補正ステップと、前記補
正ステップにおいて作成された複数の補正画像を評価す
るための評価値を算出する評価値算出ステップと、前記
評価値算出ステップにおいて前記複数の補正画像に対応
して算出された評価値に基づき前記原画像を補正するの
に適した補正量を決定し、この決定した補正量により前
記目標の移動に伴う前記目標のレンジの変化量を推定す
る変化量推定ステップと、前記変化量推定ステップにお
いて推定した前記目標のレンジの変化量に応じて、前記
受信信号列のレンジ補償を行うレンジ補償ステップとを
有するものである。
【0053】この発明に係る補償方法は、移動する目標
に対して電波を送信し、前記目標からの反射波を受信し
て前記目標の画像を得るため、前記目標の移動に伴う前
記反射波の位相の変化を補償する補償方法であって、前
記反射波の受信信号列を基に、前記目標の移動に伴う前
記反射波のドップラー周波数の変化を示す原画像を作成
する原画像作成ステップと、前記原画像を複数の補正量
によって補正することにより複数の補正画像を作成する
補正ステップと、前記補正ステップにおいて作成された
複数の補正画像を評価するための評価値を算出する評価
値算出ステップと、前記評価値算出ステップにおいて前
記複数の補正画像に対応して算出された評価値に基づき
前記原画像を補正するのに適した補正量を決定し、この
決定した補正量により前記目標の移動に伴う前記反射波
のドップラー周波数の変化量を推定する変化量推定ステ
ップと、前記変化量推定ステップにおいて推定した前記
反射波のドップラー周波数の変化量に応じて、前記受信
信号列の位相補償を行う位相補償ステップとを有するも
のである。
【0054】この発明に係るレーダ装置は、移動する目
標に対して電波を送信する送信手段と、前記目標からの
反射波を受信する受信手段と、前記受信手段により受信
された反射波から得られた受信信号列をレンジ圧縮する
レンジ圧縮手段と、前記レンジ圧縮手段においてレンジ
圧縮された受信信号列を基に、前記目標の移動に伴う前
記目標のレンジの変化を示す原画像を作成する原画像作
成手段と、前記原画像を複数の補正量によって補正する
ことにより複数の補正画像を作成する補正手段と、前記
補正手段によって作成された複数の補正画像を評価する
ための評価値を算出する評価値算出手段と、前記評価値
算出手段によって前記複数の補正画像に対応して算出さ
れた評価値に基づき前記原画像を補正するのに適した補
正量を決定し、前記決定した補正量により前記目標の移
動に伴う前記目標のレンジの変化量を推定する変化量推
定手段と、前記変化量推定手段によって推定された前記
目標のレンジの変化量に応じて、前記レンジ圧縮手段に
よってレンジ圧縮された受信信号列のレンジ補償を行う
レンジ補償手段と、前記レンジ補償手段によりレンジ補
償がなされた信号列の中から、注目するレンジビンの信
号列を抽出する信号列抽出手段と、前記信号列抽出手段
において抽出された信号列をクロスレンジ圧縮するクロ
スレンジ圧縮手段とを有するものである。
【0055】この発明に係るレーダ装置は、移動する目
標に対して電波を送信する送信手段と、前記目標からの
反射波を受信する受信手段と、前記受信手段により受信
された反射波から得られた受信信号列をレンジ圧縮する
レンジ圧縮手段と、前記レンジ圧縮手段によりレンジ圧
縮された受信信号列の内、注目するレンジビンの信号列
を抽出する信号列抽出手段と、前記信号列抽出ステップ
により抽出された信号列に基づき、前記目標の移動に伴
う前記反射波のドップラー周波数の変化を示す原画像を
作成する原画像作成手段と、前記原画像を複数の補正量
によって補正することにより複数の補正画像を作成する
補正手段と、前記補正手段によって作成された複数の補
正画像を評価するための評価値を算出する評価値算出手
段と、前記評価値算出手段によって前記複数の補正画像
に対応して算出された評価値に基づいて前記原画像を補
正するのに適した補正量を決定し、この決定した補正量
により前記目標の移動に伴う前記反射波のドップラー周
波数の変化量を推定する変化量推定手段と、前記変化量
推定手段によって推定した前記反射波のドップラー周波
数の変化量に応じて、前記信号列抽出手段により抽出さ
れた信号列の位相補償を行う位相補償手段と、前記位相
補償手段により位相補償された信号列をクロスレンジ圧
縮するクロスレンジ圧縮手段とを有するものである。
【0056】
【発明の実施の形態】
実施の形態1.図1は、レンジ追尾手段の実施の一形態
を示す図である。図1において、111はレンジ圧縮手
段5から出力された信号を入力するピーク画像生成手段
であり、112はピーク画像生成手段111の出力側に
設けられた画像線形補正手段であり、113は画像線形
補正手段112の出力側に設けられた射影セル数算出手
段であり、114は射影セル数算出手段113の出力側
に設けられた推定シフト量算出手段である。
【0057】115は、ピーク画像生成手段111、画
像線形補正手段112、射影セル数算出手段113、推
定シフト量算出手段114を含んで構成される上記シフ
ト量推定回路である。この実施の形態におけるレンジ追
尾手段501は、従来における振幅最大レンジビン検出
回路502、平滑化回路503の代わりにシフト量推定
回路115を備えている点で相違している。このシフト
量推定回路115は、二次元画像中に存在する同じ傾き
を有する直線群の傾きを推定する回路であり、図23で
示したシフト量sを推定する。
【0058】図2は、ピーク画像生成手段111の処理
内容を説明するための原理図である。図において、11
6は原画像、117はピーク画像である。原画像および
ピーク画像については、後で詳述する。
【0059】図3は、画像線形補正手段112の処理内
容を説明するための原理図である。図3において、11
8は補正量0の時の補正後のピーク画像、119は補正
量−4の時の補正後のピーク画像、120は補正量−8
の時の補正後のピーク画像、121は補正量−10の時
の補正後のピーク画像である。
【0060】図4は、射影セル数算出手段113で得ら
れる射影図である。図において、122は、補正量0の
時の射影図、123は、補正量−4の時の射影図、12
4は、補正量−8の時の射影図、125は、補正量−1
0の時の射影図である。
【0061】図5は、射影セル数算出手段113で算出
した、各補正量に対する射影セル数を示す図である。次
に、図1、2、3、4、5、17、18、23を用い
て、本実施の形態の内容を説明する。
【0062】送信機1で高周波パルスを生成し、目標に
照射し、受信信号をレンジ圧縮手段5でレンジ圧縮する
までの処理、及び、レンジ追尾手段501で得られた追
尾結果であるシフト量を用いて、レンジ補償量算出手段
504で、レンジ補償量を算出する部分から、目標のI
SAR画像をモニタTV9に表示するまでの処理は、従
来の技術と同一である。
【0063】以下、レンジ追尾手段501の処理内容に
ついて、順を追って説明する。まず、シフト量推定回路
115の内容について説明する。 <ST01:原画像作成ステップ>ピーク画像生成手段
111は、レンジ圧縮手段5でレンジ圧縮を行うことに
よって得られた画像信号を入力する。レンジ圧縮手段5
によるレンジ圧縮によって得られる画像信号は、複素信
号であるためピーク画像生成手段111は各画素の振幅
値を検出することにより図2(a)に示す原画像116
を作成する。
【0064】ここで、原画像116は、横軸x方向にx
num、縦軸y方向にynumの大きさを持つ二次元画
像であり、画像中の各画素についての振幅を濃淡表示
(白→黒になるに従って振幅が大きくなる)している。
この振幅は、各画素の複素信号の絶対値を表すものであ
る。そして、レンジプロフィールのヒストリを考える場
合には、x軸がレンジ方向の軸に対応し、y軸がヒット
方向の軸に対応する。即ち、原画像g(x,y)は目標
物の移動に伴うレンジ方向の変化の様子を示しているこ
とになる。また、原画像116中に含まれる直線の傾き
は、目標物のレンジの変化量であるシフト量で与えてお
り、このシフト量を推定する必要がある。
【0065】本実施の形態では、以降の処理により、原
画像116に含まれる直線のレンジ方向のシフト量を推
定する、即ち、図2(a)におけるx軸方向のシフト量
を推定する。 <ST02:ピーク画像作成ステップ>ピーク画像生成
手段111は、原画像116のピーク画像117を生成
する。ここで、ピーク画像117とは、各yごとに、原
画像116の振幅を最大とするセル(以下では、ピーク
セルと呼ぶ。)の値を1、それ以外のセルを0とした二
値化画像のことを指す。図2(b)に示したピーク画像
117においては、各yごとに選ばれたピークセルを黒
で示し、ピークセル以外の部分を白で示している。11
7からわかるように、ある一定の傾きを持った複数の直
線が現れることになる。
【0066】<ST03:補正ステップ>次に、画像線
形補正手段112は、ピーク画像生成手段111によっ
て得られたピーク画像117について、図中に存在する
複数の直線の傾きを変えるための補正処理を施す。画像
線形補正手段112では、y=ynum─1セルにおけ
るx方向の補正量をhsとすると、各yにおける補正量
hx(y)を、式(8)で計算する。
【0067】
【数8】
【0068】ここで、式(8)中のAは適当な自然数と
し、通常は1〜2程度の値を用いれば良い。画像線形補
正手段112は、仮の補正量hsを複数設定し、仮の補
正量hsそれぞれを用いた場合のピーク画像117の補
正結果(以下、補正画像と記す)を得る。仮の補正量h
sをそれぞれ、0、―4、−8、─10とした場合の補
正画像を図3の118、119、120、121に示
す。
【0069】<ST04:射影ステップ>各補正量hs
に従い補正した後の各ピーク画像に対して、射影セル数
算出手段113では、まず、射影図を生成する。射影図
とは、補正後の各ピーク画像を、各セルxごとに分析
し、各セルxにおいて、ピークセルが一つでも存在した
場合には1、一つも存在しなかった場合には0とするよ
うな処理を施した結果得られる図をいう。以降の図4で
は、ピークセルが一つでも存在した場合を黒(1に対
応)で示し、ピークセルが一つも存在しなかった場合を
白(0に対応)で示している。
【0070】図4は、各補正量0、─4、─8、─10
に対応して生成された射影図122、123、124、
125を示している。ここで、射影図において、値が1
となるセル(すなわち、ピークセルが一つでも射影され
たセル)の数が少なければ少ないほど、ピーク画像中の
各直線がx軸に垂直になるように補正されていると考え
られる。
【0071】<ST05:評価値算出ステップ>そこ
で、射影セル数算出手段113は、作成した各射影図の
内値が1となったセル数を評価値として算出する。図4
に示した射影図122、123、124、125におい
ては、射影図122はセル数11になり、射影図123
はセル数7になり、射影図124はセル数3になり、射
影図125はセル数5となる。射影セル数算出手段11
3は、算出したセル数および対応する補正量を推定シフ
ト量算出手段114に対して通知する。
【0072】<ST06:変化量推定ステップ>推定シ
フト量算出手段114では、最もセル数の少ない射影図
に対応する補正量に基づいてシフト量を算出する。具体
的には、推定シフト量算出手段114は射影セル数算出
手段113から通知された算出したセル数および対応す
る補正量とから、図5に示すようなグラフを作成する。
図5における縦軸の評価指標とは、射影セル数算出手段
113によって算出されたセル数を示しており、横軸は
対応する補正量を示している。
【0073】推定シフト量算出手段114は、図5に示
したグラフから評価指標を最小とする補正量が−8であ
ることを認識し、この補正量に―の符号をかけた値であ
る8をシフト量として算出する。このようにして得られ
たシフト量が、原画像116中の各直線の傾きの推定結
果であり、こうして目標物のレンジの変化量が推定され
る。推定シフト量算出手段114は、算出したシフト量
をレンジ補償量算出手段504に対して出力する
【0074】<ST07:レンジ補償ステップ>そし
て、レンジ補償量算出手段504は、入力したシフト量
に応じ、式(4)を使ってレンジ補償量を算出する。以
降の処理については、従来技術と同様である。
【0075】以上の処理を用いることにより、レンジ追
尾を行なう際に、目標上の各点の反射信号の振幅値が時
間の経過に伴い変動することによって、振幅最大点がレ
ンジ軸で急に移動した場合であっても、これらの影響に
よって生じるレンジ追尾誤差を抑えてレンジ追尾を行う
ことができ、より鮮明なレーダ画像を得ることができ
る。
【0076】また、推定シフト量算出手段114は、射
影セル数算出手段113から通知された算出したセル数
および対応する補正量とから、図5に示すようなグラフ
を作成し、評価指標である射影されたセル数の変化に基
づいて補正量を得るようにしているので、適正な補正量
を選ぶことができる。
【0077】尚、図5では、射影されたセル数をプロッ
トして各プロットを直線で結ぶことにより適切な補正量
を選んでいるが、評価指標のグラフのピーク近傍の数点
をスプライン補間することにより、評価指標を最小とす
る補正量をより精度良く求めることがもできる。
【0078】また尚、この実施の形態においては画像線
形補正手段112は、仮の補正量hsを複数設定し、仮
の補正量hsそれぞれを用いた場合のピーク画像117
の補正結果に基づいて最適な補正量を得ることとしてい
るが、次にようにして最適な補正量を得るようにしても
よい。まず、粗めに仮の補正量hsを複数設定し、この
補正量それぞれに応じたピーク画像117に基づいて、
設定した補正量の中から適切な補正量hsを得る。次
に、この補正量hsに近い値の補正量を細かく設定し、
この細かく設定した補正量hsについてピーク画像を補
正し、最適な補正量を得る。このように、始めは粗く、
次第に細かく段階的に補正量を設定することによって、
補正量を算出するために必要な処理数を削減することが
できる。
【0079】尚、この実施の形態においては、レンジプ
ロフィールのヒストリをヒットとレンジの二次元画像で
示しているが、一般にヒットは時間に比例するため、時
間とレンジの二次元画像で示すようにすることもでき
る。
【0080】尚、この実施の形態におけるピーク画像作
成手段111の内の振幅検出を行う部分は原画像作成手
段の一例であり、ピーク画像作成手段111のピーク画
像を作成する部分はニ値化画像作成手段の一例である。
また、射影セル数算出手段113は評価値算出手段の一
例であり、推定シフト量算出手段114は変化量推定手
段の一例である。また、この実施の形態では補正画像毎
に射影されたセル数を算出するようにしているが、この
射影されたセル数は補正画像を評価するための評価値の
一例である。したがって、補正画像を評価する評価値と
して他の値を算出するようにしてもよい。
【0081】実施の形態2.図6は、ドップラー追尾手
段の実施の一形態を示す図である。図において、115
は図1に示したシフト量推定回路と同様であり、508
は図19に示したドップラー追尾手段と同様であるの
で、具体的構成については説明を省略する。
【0082】次に図2、6、17、19を用いて本実施
の形態の内容を説明する。送信機1で高周波パルスを生
成し、目標に照射し、受信信号をレンジ圧縮手段5でレ
ンジ圧縮し、これをレンジ補償回路6でレンジ補償する
までの処理は、従来例もしくは実施の形態1と同一であ
る。
【0083】また、レンジ補償後のレンジプロフィール
のヒストリに対して、注目レンジビン決定手段506で
注目レンジビンの受信データ列を抽出し、区分周波数分
析手段507で、区分周波数分布のヒストリを算出する
までの処理と、ドップラー追尾手段508で得られたド
ップラー周波数のシフト量の推定結果から、位相補償量
算出手段510で位相補償量を算出し、これを用いて、
レンジ補償後のレンジプロフィールのヒストリを位相補
償手段511で位相補償し、クロスレンジ圧縮手段8で
ISAR画像を生成しモニタTV9に表示するまでの処
理は、従来例と同一であるので説明は省略する。
【0084】ドップラー追尾手段508には、実施の形
態1で説明したシフト量推定回路115が含まれてい
る。ピーク画像生成手段111には、区分周波数分析手
段507による周波数分析によって得られた周波数分布
のヒストリを示す信号が入力される。この入力された信
号は複素信号であるため、ピーク画像作成手段111が
各画素ごとに振幅を検出することにより図2(a)に示
すような原画像116を作成する。ただし、得られる原
画像116におけるX軸がドップラー周波数であり、y
軸がヒットとなる。即ち、原画像116は反射波のドッ
プラー周波数の変化の様子を示すものである。
【0085】本実施の形態では、原画像116における
x軸をドップラー周波数、y軸をヒットとすることによ
り、実施の形態1でレンジプロフィールのヒストリのシ
フト量を推定したのと同様に、区分周波数分布のヒスト
リのシフト量を推定する。即ち、実施の形態1で説明し
た通り、シフト量推定回路115は、入力した原画像1
16に対して、この画像中に存在する直線群の傾きを推
定する。区分周波数分布のヒストリのシフト量を推定す
る手順については、先の実施の形態と同様であるので説
明は省略する。
【0086】以上の処理を用いることにより、ドップラ
ー追尾を行なう際に、目標上の各点の反射信号の振幅が
時間の経過に伴い変動したり、ドップラー周波数の折り
返しが生じたりした場合であっても、これらの影響によ
って生じるドップラー追尾誤差を抑えてドップラー追尾
を行うことができ、より鮮明なレーダ画像を得ることが
できる。
【0087】また尚、この実施の形態においては、区分
周波数分布ののヒストリをヒットとドップラー周波数の
二次元画像で示しているが、一般にヒットは時間に比例
するため、時間とドップラー周波数の二次元画像で示す
ようにすることもできる。尚、実施の形態1ではレンジ
補償のみについて説明し、実施の形態2では位相補償の
みについて説明したが、実施の形態1におけるレンジ補
償と実施の形態2における位相補償の両方を行うように
レーダを構成することも可能であることはいうまでもな
い。
【0088】実施の形態3.図7は、本実施の形態のシ
フト量推定回路の構成図である。図において、131は
多点ピーク画像生成手段である。112、113、11
4、115は、図1と同一である。図8は、多点ピーク
画像生成手段131で得られた多点ピーク画像である。
図において、132は多点ピーク画像である。ここで、
多点ピークとは、各yに対して選択される点の振幅が最
大振幅からN番目以内に入っている点をいい、このNは
適宜設定される。即ち、多点ピークは、振幅値が大きい
方からN番目以内にあることを示すものである。
【0089】次に、図2、7、8を用いて、本実施の形
態の処理内容を説明する。図7のシフト量推定回路にお
いて、画像線形補正手段112、射影セル数算出手段1
13、推定シフト量算出手段114の処理内容は実施の
形態1、2と同一である。多点ピーク画像生成手段13
1は、入力されたレンジプロフィールのヒストリや、ド
ップラー周波数分布のヒストリより原画像116を生成
し、各yごとに、最大振幅より順にN個のピークセルを
選択する。そして、選択されたピークセルに基づいて図
8に示す多点ピーク画像132を生成する。ただし、図
8の例では、N=2の場合を示しており、多点ピークと
して選択されたセルの値を1、選択されなかったセルの
値を0としたニ値化画像である。
【0090】以降は、多点ピーク画像生成手段131に
よって作成された多点ピーク画像132を用いて、画像
線形補正手段112による処理、射影セル数算出手段1
13による処理、推定シフト量算出手段114による処
理を行う。
【0091】多点ピーク画像132を用いた場合には、
射影セル数算出手段113で得られた評価指標の各補正
量に対する特性が、ピーク画像117を用いた場合に比
べて、画像中の複数の直線がx軸と直交する補正量近傍
でより急峻な変化となる。したがって、評価指標の最小
点をより精度良く求めることができる。すなわち、図1
や図6のシフト量推定回路で本実施の形態の処理を行な
うことにより、実施の形態1、2の効果が得られるのみ
ならず、さらに、上記理由でレンジ追尾、ドップラー追
尾の精度を向上させることができる。
【0092】実施の形態4.図9は、本実施の形態のシ
フト量推定回路の構成図である。図において、141は
画像の各セルの総和算出手段、142はエントロピー算
出手段である。111、112、114、115は図1
と同一である。
【0093】図10は、画像の各セルの総和算出手段1
41の処理内容を示す図である。図において、143は
補正量0の場合のy方向の総和を計算する過程の図、1
44は補正量─4の場合のy方向の総和を計算する過程
の図、145は補正量─8の場合のy方向の総和を計算
する過程の図、146は補正量─10の場合のy方向の
総和を計算する過程の図である。
【0094】図11は推定シフト量算出手段114でシ
フト量を算出するための評価指標の補正量に対する変化
を示す図である。次に、図2、3、9、10、11を用
いて本実施の形態の処理内容を説明する。図9のシフト
量推定回路115において、ピーク画像生成手段111
でレンジプロフィールのヒストリやドップラー周波数分
布のヒストリなどを表す原画像116を作成し、さらに
ピーク画像117を生成する処理、得られたピーク画像
117に対して、想定した各補正量ごとに、画像線形補
正手段112で補正を行なう処理は、実施の形態1.と
同一である。補正結果が、図3における118、11
9、120、121に対応する。
【0095】画像の各セルの総和算出手段141では、
ピーク画像の各補正後の画像118、119、120、
121について、各xごとに、y方向の全セルの電力の
総和の計算を行なう。画像の各セルの総和算出手段14
1によって算出された電力の総和の結果を図10に示
す。このx軸は、レンジ又はドップラー周波数を示して
おり、したがって、この図は各レンジでの総電力又は各
ドップラー周波数での総電力を示していることになる。
図10における143、144、145、146は、各
補正量に対応した各xセル毎の電力(ピーク画像の場合
はピークセル数)の総和の分布t(x)を示している。
【0096】ところで、分布の偏りを評価する指標にエ
ントロピーがある。エントロピーは、画像の情報量を定
量的に示す指標としてしばしば用いられ、ディジタル画
像を、その情報量を損なわずに符号化するために必要な
最も短い符号長を与えることが知られている。その定義
は式(9)で与えられる。
【0097】
【数9】
【0098】ここで、PWRは入力したピーク画像の総
電力、すなわち、ピークセルの数である。エントロピー
はその定義からもわかるように、分布が一様であれば、
エントロピーが大きな値を示し、逆に分布に偏りが有る
とエントロピーが減少する。従って、エントロピーを指
標として、これを最小とする補正量を探索することで、
各点の反射信号列とx軸が直交している場合をみつける
事ができる。そして、このときの補正量に―をかけるこ
とにより、シフト量を得ることができる。
【0099】エントロピー算出手段142は、図10で
得られた各補正量に対するピーク画像の総和の分布t
(x)に対して式(9)によりエントロピーHを算出す
る。各補正量に対するエントロピーを算出した結果を図
11に示す。ここで、横軸は補正量、縦軸は評価指標で
あるエントロピーを示す。エントロピー算出手段142
は、算出結果を推定シフト量算出手段114に対して通
知する。
【0100】推定シフト量算出手段114は、エントロ
ピー算出手段142から通知された算出結果から評価指
標を最小とする補正量が−8であることを認識し、シフ
ト量の推定結果を8を得る。これは、原画像116に与
えたシフト量の値と等しい。そして、推定シフト量算出
手段114は、得られた推定シフト量をレンジ補償量算
出手段504又は位相補償量算出手段510に対して出
力する。
【0101】以上の処理を行なうことにより、雑音が多
い環境で、レンジ追尾における振幅最大レンジビン、も
しくはドップラー追尾における振幅最大ドップラービン
として、雑音により振幅が大きくなったセルを選択して
しまう確率が高くなる結果、これらを用いて平滑化処理
を行なう際の平滑誤差が生じる、すなわち、レンジ追
尾、ドップラ─追尾の精度が劣化してしまうと言う問題
点を除去できる。
【0102】尚、この実施の形態では、式(9)を用い
てエントロピーHを算出するようにしているが、この式
以外を用いてもよい。明細書において、エントロピーと
は、分布の偏りを評価する一つの指標であり、分布の偏
りを評価することができるものであれば式(9)以外の
式を用いて算出してもよい。
【0103】尚、この実施の形態における画像の各セル
の総和算出手段141及びエントロピー算出手段142
は、評価値算出手段の一例である。また、ピーク画像作
成手段111の内の振幅検出を行い原画像を作成する部
分は、原画像作成手段に相当する。
【0104】実施の形態5.図12は、本実施の形態の
シフト量推定回路の構成図である。図において、11
2、114、115、141、142は図9と同一であ
る。次に、図12を用いて、本実施の形態の処理内容を
説明する。実施の形態4では、入力信号に対して、ま
ず、ピーク画像生成手段111によってピーク画像を生
成し、このピーク画像を用いて画像線形補正手段112
で補正処理を施していたが、この実施の形態では振幅検
出手段110によって入力した信号の振幅値を検出する
ことにより二次元の原画像を作成し、画像線形補正手段
112で補正処理する。そして、画像の各セルの総和算
出手段141では、各補正量に対して得られた補正後の
画像について、各セルxにおけるy方向の電力の総和を
計算する。
【0105】得られた各補正量に対する電力分布に対し
て、エントロピー算出手段142で評価指標であるエン
トロピーを算出する処理、各補正量に対して得られた評
価指標から、推定シフト量算出手段114でシフト量を
推定する処理は実施の形態5と同一である。
【0106】すなわち、実施の形態4ではピークセルの
分布に注目して追尾を行なったのに対して、本実施の形
態では、より情報量の多い画像の電力分布に注目して追
尾を行なうので、追尾精度がさらに向上する。尚、この
実施の形態では、式(9)を用いてエントロピーHを算
出するようにしているが、この式以外を用いてもよい。
明細書において、エントロピーとは、分布の偏りを評価
する指標であり、分布の偏りを評価することができるも
のであれば式(9)以外の式を用いて算出してもよい。
【0107】尚、この実施の形態では振幅検出手段11
0によって入力信号の振幅値を検出してから画像線形補
正手段112による補正を行う旨説明したが、画像線形
補正手段112によって入力信号に補正処理を施した後
に振幅検出手段110で振幅検出を行い、さらに画像の
各セルの総和算出を行うように構成することも可能であ
る。
【0108】実施の形態6.図13は、本実施の形態の
レンジ補償回路の構成図である。図において、161は
ピーク画像生成手段111の出力側であり、かつ射影セ
ル数算出手段113の入力側に設けられた画像非線形補
正手段である。162は射影セル数算出手段113の出
力側に設けられた推定シフト量非線形算出手段であり、
163は推定シフト量非線形算出手段162の出力側お
よびレンジ補償手段505の入力側に設けられたレンジ
補償量非線形算出手段である。図において、6、505
は図18と、111、113は図1と同様であるので説
明は省略する。
【0109】図14は、画像非線形補正手段161の処
理内容を示す図である。図において、164は非線形補
正前のピーク画像、165は非線形補正後のピーク画像
を示す。
【0110】次に、図13、14を用いて、本実施の形
態の処理内容を説明する。レンジ補償回路6において、
ピーク画像生成手段111でピーク画像を生成する処理
は、実施の形態1と同一である。画像非線形補正手段1
61には、図14に示すピークセルのx座標がyの変化
に対して非線形に変化するようなピーク画像164が入
力するものとする。このピーク画像164に対して、画
像非線形補正手段161では、非線形な補正を行なう。
すなわち、補正を行なう際に、ピーク画像164に示す
ような、非線形な補正量を仮定する。例えば、二次元の
補正を行なう場合には、補正量sf(y)を次式で定義
する。
【0111】
【数10】
【0112】画像非線形補正手段161は、式(10)
のパラメータa、bを様々に変えて補正後のピーク画像
を生成する。射影セル数算出手段113では、実施の形
態1と同様に、各補正後のピーク画像165に対して、
射影図を作成し、これから、評価指標の値を算出する。
推定シフト量非線形算出手段162では、評価指標を最
小とする式(10)の係数a、bの組(ae、be)を
決定する。
【0113】レンジ補償量非線形算出手段163は、式
(10)に上記(ae、be)を代入して、レンジ補償
量se(y)を決定する。レンジ補償手段505は、以
上で得られたレンジ補償量se(y)を用いて、レンジ
プロフィールのヒストリのレンジ補償を行なう。
【0114】以上の処理を行なうことにより、目標が加
速度を持ち、レンジの移動が非線形になる場合にもレン
ジ追尾が行なえる。また、上記処理は、位相補償回路に
も同様に適用できることは言うまでもない。この実施の
形態において述べた補正方法は非線形補正をよばれるも
のであり、この非線形補正は先の実施の形態で説明した
すべてに適用可能である。
【0115】実施の形態7.従来のレーダ装置において
は、位相補償回路7で区分周波数分析の点数を決定する
場合に次のような問題が生じていた。 (a)区分周波数分析の点数を少なくすると、ドップラー
周波数の分解能が劣化するため、これに基づく追尾誤差
が発生する。 (b)区分周波数分析の点数を多くすると、ドップラーの
分解能は向上するものの、時間の分解能が劣化するた
め、信号のチャープ成分によりドップラーが広がって、
これに基づく追尾誤差が発生する。 この実施の形態では、上記従来の問題を解消するため
に、最初は区分周波数分析の点数を少なくしてドップラ
ー追尾を粗く行ない、注目レンジビンのデータ列の位相
補償を行なう。次に、次第に区分周波数分析の点数を増
やしていきドップラー追尾処理を繰り返していくもので
ある。
【0116】まず、この実施の形態における位相補償回
路7の構成について説明する。図15は、この実施の形
態における位相補償回路7の内部構成を示すブロック図
である。注目レンジビン決定手段506、区分周波数分
析手段507、位相補償量算出手段510、位相補償手
段511は従来例と同様であるので説明は省略する。
【0117】171は注目レンジビン決定手段506の
出力側であり、かつ区分周波分析手段507の入力側に
設けられた区分周波数分析点数更新手段であり、この区
分周波数分布点数更新手段171には区分周波数分析点
数決定用の初期値Nsttが入力される。この実施の形
態における位相補償回路7内には、位相補償手段511
が2つ設けられている点、各回路及び手段相互間の接続
関係が従来の位相補償回路とは異なっている。
【0118】172は位相補償量算出手段510の出力
側であり、かつ位相補償手段511の入力側に設けられ
ている位相補償量蓄積手段である。尚、この実施の形態
におけるシフト量推定回路116は、先の実施の形態に
おいて示したシフト量推定回路116のいずれを用いて
もよい。
【0119】図16は、本実施の形態の処理内容を示す
フローチャートである。注目レンジビン決定手段506
は、入力したレンジ補償後のレンジプロフィールのヒス
トリv2(h、r)の中から、従来技術と同じ方式で、
注目レンジビンのデータ列w(h)を出力する。区分周
波数分析点数更新手段171には、注目レンジビンのデ
ータ列w(h)と、区分周波数分析点数決定用の初期値
Nsttが入力される。この初期値Nsttは予め設定
された値である。
【0120】<ST01>区分周波数分析点数更新手段
171では、区分周波数分析点数決定用数値Nに対して
Nsttを代入するとともに、位相補償量蓄積手段51
0で総位相補償量PT(h)を0にセットする。
【0121】<ST02>区分周波数分析点数更新手段
171は、区分周波数分析点数Mを式(11)により決
定する。ただし、ここで、hnumは総データ長であ
る。
【0122】
【数11】
【0123】<ST03>区分周波数分析手段507
で、注目レンジビンのデータ列w(h)に対してM点の
区分周波数分析を行ない、区分周波数分布のヒストリを
得る。
【0124】<ST04>区分周波数分析手段507に
よって得られた区分周波数分布のヒストリに対して、シ
フト量推定回路115は、実施の形態1〜6と同様にそ
の周波数のシフト量を推定する。そして、得られた推定
シフト量を位相補償量算出手段510に対して通知す
る。位相補償量算出手段510は、シフト量推定回路1
15で得られたシフト量から、位相補償量P(h)を計
算する。
【0125】<ST05>位相補償手段511は、位相
補償量算出手段510によって算出された位相補償量P
(h)によって注目レンジビンのデータ列w(h)を位
相補償する。 <ST06>位相補償量蓄積手段172は、位相補償量
算出手段510によって算出された位相補償量P(h)
を式(12)に代入して、総位相補償量PT(h)に加
算する。
【0126】
【数12】
【0127】<ST07>区分周波数分析点数更新手段
171は、現在の区分周波数分析点数決定用数値Nを調
べ、もし、Nが1ならば、<ST09>へ進む。現在の
区分周波数分析点数決定用数値Nが1でなければ<ST
08>へ進む。
【0128】<ST08>現在の区分周波数分析点数決
定用数値Nが1でない場合には、Nの値を一つ減らし、
<ST02>へ戻り、再度ST03〜ST05の処理を
実行する。
【0129】<ST09>位相補償手段511で、レン
ジ補償後のレンジプロフィールのヒストリを、位相補償
量蓄積手段172に蓄えられた総位相補償量PT(h)
を用いてレンジ補償後のレンジプロフィールのヒストリ
v2(h、r)を位相補償する。
【0130】この実施の形態では、最初は少ない区分周
波数分析の点数でドップラー追尾を粗く行なって注目レ
ンジビンのデータ列の位相補償を行ない、次第に区分周
波数分析の点数を増やして同様の処理を繰り返し、ドッ
プラーの時間変化を逐次的に推定するということによ
り、上述のような問題を克服することができるため、よ
り鮮明な画像が得られるようになる。尚、この実施の形
態では、ST05の後にST06を行っているが、ST
06の後にST05を行うようにしてもよい。また尚、
この実施の形態のST04において位相補償量を算出す
る方法としては、先の実施の形態において説明した方法
以外の方法、例えば従来の方法等を適用することもでき
る。
【0131】
【発明の効果】この発明は、以上に説明したように構成
されているので、以下に記載されるような効果を奏す
る。この発明に係る補償回路は、移動する目標に対して
電波を送信し、前記目標からの反射波を受信して前記目
標の画像を得るため、前記目標の移動に伴う前記目標の
レンジの変化を補償する補償回路であって、前記反射波
の受信信号列を基に、移動に伴う前記目標のレンジの変
化を示す原画像を作成する原画像作成手段と、複数の補
正量によって前記原画像を補正することにより複数の補
正画像を作成する補正手段と、前記補正手段によって作
成された複数の補正画像を評価するための評価値を算出
する評価値算出手段と、前記評価値算出手段によって前
記複数の補正画像に対応して算出された評価値に基づき
前記原画像を補正するのに適した補正量を決定し、この
決定した補正量により前記目標の移動に伴う前記目標の
レンジの変化量を推定する変化量推定手段と、前記変化
量推定手段によって推定した前記目標のレンジの変化量
に応じて前記受信信号列のレンジ補償を行うレンジ補償
手段とを有するため、目標の移動に伴い目標のレンジが
変化する場合であってもレンジ追尾誤差を押さえてレン
ジ追尾を行なうことができる。
【0132】この発明に係る補償回路は、移動する目標
に対して電波を送信し、前記目標からの反射波を受信し
て前記目標の画像を得るため、前記目標の移動に伴う前
記反射波の位相の変化を補償する補償回路であって、前
記反射波の受信信号列を基に、前記目標の移動に伴う前
記反射波のドップラー周波数の変化を示す原画像を作成
する原画像作成手段と、前記原画像を複数の補正量によ
って補正することにより複数の補正画像を作成する補正
手段と、前記補正手段によって作成された複数の補正画
像を評価するための評価値を算出する評価値算出手段
と、前記評価値算出手段によって前記複数の補正画像に
対応して算出された評価値に基づき前記原画像を補正す
るのに適した補正量を決定し、この決定した補正量によ
り前記目標の移動に伴う前記反射波のドップラー周波数
の変化量を推定する変化量推定手段と、前記変化量推定
手段によって推定された前記反射波のドップラー周波数
の変化量に応じて、前記受信信号列の位相補償を行う位
相補償手段とを有するため、目標の移動に伴いドップラ
ー周波数が変化する場合であっても、追尾誤差を押さえ
てドップラー追尾を行なうことができる。
【0133】また、前記評価値算出手段が、前記補正画
像の各レンジ毎に、同一のレンジに存在している複数の
画素の値を積分し、当該積分値の分布を示す指標を評価
値として算出するため、雑音の影響をより低減して精度
良い追尾か可能となる。
【0134】さらに、前記評価値算出手段が、前記補正
画像の各ドップラー周波数毎に、同一のドップラー周波
数に対して存在している複数の画素の値を積分し、当該
積分値の分布を示す指標を評価値として算出するため、
画素の値を用いてドップラー追尾を行なえるので、追尾
精度が向上する。
【0135】この発明に係る補償回路は、移動する目標
に対して電波を送信し、前記目標からの反射波を受信し
て前記目標の画像を得るため、前記目標の移動に伴う前
記目標のレンジの変化を補償する補償回路であって、前
記反射波の受信信号列を基に、移動に伴う前記目標のレ
ンジの変化を示す原画像を作成する原画像作成手段と、
前記原画像作成手段によって作成された原画像を二値化
することによりニ値化画像を作成する画像二値化手段
と、前記二値化画像を複数の補正量によって補正するこ
とにより複数の補正ニ値化画像を作成する補正手段と、
前記補正手段によって作成された複数の補正二値化画像
を射影することにより複数の射影画像を作成し、この複
数の射影画像に対する評価値を算出する評価値算出手段
と、前記評価値算出手段によって前記複数の射影画像に
対応して算出された評価値に基づき前記原画像を補正す
るのに適した補正量を決定し、この決定した補正量によ
り前記目標の移動に伴う前記目標のレンジの変化量を推
定する変化量推定手段と、前記変化量推定手段によって
推定した前記目標のレンジの変化量に応じて、前記受信
信号列のレンジ補償を行うレンジ補償手段とを有するた
め、観測時間中に振幅最大反射点の移動が生じても精度
良いレンジ追尾を行うことができる。
【0136】この発明に係る補償回路は、移動する目標
に対して電波を送信し、前記目標からの反射波を受信し
て前記目標の画像を得るため、前記目標の移動に伴う前
記反射波の位相の変化を補償する補償回路であって、前
記反射波の受信信号列を基に、移動に伴う前記反射波の
ドップラー周波数の変化を示す原画像を作成する原画像
作成手段と、前記原画像作成手段によって作成された原
画像を二値化することによりニ値化画像を作成する画像
二値化手段と、前記二値化画像を複数の補正量によって
補正することにより複数の補正ニ値化画像を作成する補
正手段と、前記複数の補正二値化画像を射影することに
より複数の射影画像を作成し、この複数の射影画像に対
する評価値を算出する評価値算出手段と、前記評価値算
出手段によって前記複数の射影画像に対応して算出され
た評価値に基づいて前記原画像を補正するのに適した補
正量を決定し、この決定した補正量により前記目標の移
動に伴う前記目標のレンジの変化量を推定する変化量推
定手段と、前記変化量推定手段によって推定した前記反
射波のドップラー周波数の変化量に応じて、前記受信信
号列の位相補償を行う位相補償手段とを有するため、観
測時間中に振幅最大反射点の移動が生じても精度良いド
ップラー追尾を行うことができる。
【0137】また、前記評価値算出手段が、前記射影画
像内の特定の値を有する画素の画素数を評価値として算
出するため、評価値を容易に算出することができる。
【0138】さらに、前記補正手段が、前記複数の補正
量に応じて前記原画像を線形補正するため、目標の移動
に伴うレンジ又はドップラー周波数の変化が線形になる
場合に精度よいレンジ追尾やドップラー追尾が行なえ
る。
【0139】さらにまた、前記補正手段が、前記複数の
補正量に応じて前記原画像を非線形補正するため、目標
の移動に伴うレンジ又はドップラー周波数の変化が非線
形になる場合に精度よいレンジ追尾やドップラー追尾が
行なえる。
【0140】この発明に係る補償方法は、移動する目標
に対して電波を送信し、前記目標からの反射波を受信し
て前記目標の画像を得るため、前記目標の移動に伴う前
記目標のレンジの変化を補償する補償方法であって、前
記反射波の受信信号列を基に、移動に伴う前記目標のレ
ンジの変化を示す原画像を作成する原画像作成ステップ
と、前記原画像を複数の補正量によって補正することに
より複数の補正画像を作成する補正ステップと、前記補
正ステップにおいて作成された複数の補正画像を評価す
るための評価値を算出する評価値算出ステップと、前記
評価値算出ステップにおいて前記複数の補正画像に対応
して算出された評価値に基づき前記原画像を補正するの
に適した補正量を決定し、この決定した補正量により前
記目標の移動に伴う前記目標のレンジの変化量を推定す
る変化量推定ステップと、前記変化量推定ステップにお
いて推定した前記目標のレンジの変化量に応じて、前記
受信信号列のレンジ補償を行うレンジ補償ステップとを
有するため、目標の移動に伴い目標のレンジが変化する
場合であってもレンジ追尾誤差を押さえてレンジ追尾を
行なうことができる。
【0141】この発明に係る補償方法は、移動する目標
に対して電波を送信し、前記目標からの反射波を受信し
て前記目標の画像を得るため、前記目標の移動に伴う前
記反射波の位相の変化を補償する補償方法であって、前
記反射波の受信信号列を基に、前記目標の移動に伴う前
記反射波のドップラー周波数の変化を示す原画像を作成
する原画像作成ステップと、前記原画像を複数の補正量
によって補正することにより複数の補正画像を作成する
補正ステップと、前記補正ステップにおいて作成された
複数の補正画像を評価するための評価値を算出する評価
値算出ステップと、前記評価値算出ステップにおいて前
記複数の補正画像に対応して算出された評価値に基づき
前記原画像を補正するのに適した補正量を決定し、この
決定した補正量により前記目標の移動に伴う前記反射波
のドップラー周波数の変化量を推定する変化量推定ステ
ップと、前記変化量推定ステップにおいて推定した前記
反射波のドップラー周波数の変化量に応じて、前記受信
信号列の位相補償を行う位相補償ステップとを有するた
め、目標の移動に伴い反射波のドップラー周波数が変化
する場合であっても追尾誤差を押さえてドップラー追尾
を行なうことができる。
【0142】この発明に係るレーダ装置は、移動する目
標に対して電波を送信する送信手段と、前記目標からの
反射波を受信する受信手段と、前記受信手段により受信
された反射波から得られた受信信号列をレンジ圧縮する
レンジ圧縮手段と、前記レンジ圧縮手段においてレンジ
圧縮された受信信号列を基に、前記目標の移動に伴う前
記目標のレンジの変化を示す原画像を作成する原画像作
成手段と、前記原画像を複数の補正量によって補正する
ことにより複数の補正画像を作成する補正手段と、前記
補正手段によって作成された複数の補正画像を評価する
ための評価値を算出する評価値算出手段と、前記評価値
算出手段によって前記複数の補正画像に対応して算出さ
れた評価値に基づき前記原画像を補正するのに適した補
正量を決定し、前記決定した補正量により前記目標の移
動に伴う前記目標のレンジの変化量を推定する変化量推
定手段と、前記変化量推定手段によって推定された前記
目標のレンジの変化量に応じて、前記レンジ圧縮手段に
よってレンジ圧縮された受信信号列のレンジ補償を行う
レンジ補償手段と、前記レンジ補償手段によりレンジ補
償がなされた信号列の中から、注目するレンジビンの信
号列を抽出する信号列抽出手段と、前記信号列抽出手段
において抽出された信号列をクロスレンジ圧縮するクロ
スレンジ圧縮手段とを有するため、目標の移動に伴い目
標のレンジが変化する場合であってもレンジ追尾誤差を
押さえてレンジ追尾を行なうことができる。
【0143】この発明に係るレーダ装置は、移動する目
標に対して電波を送信する送信手段と、前記目標からの
反射波を受信する受信手段と、前記受信手段により受信
された反射波から得られた受信信号列をレンジ圧縮する
レンジ圧縮手段と、前記レンジ圧縮手段によりレンジ圧
縮された受信信号列の内、注目するレンジビンの信号列
を抽出する信号列抽出手段と、前記信号列抽出ステップ
により抽出された信号列に基づき、前記目標の移動に伴
う前記反射波のドップラー周波数の変化を示す原画像を
作成する原画像作成手段と、前記原画像を複数の補正量
によって補正することにより複数の補正画像を作成する
補正手段と、前記補正手段によって作成された複数の補
正画像を評価するための評価値を算出する評価値算出手
段と、前記評価値算出手段によって前記複数の補正画像
に対応して算出された評価値に基づいて前記原画像を補
正するのに適した補正量を決定し、この決定した補正量
により前記目標の移動に伴う前記反射波のドップラー周
波数の変化量を推定する変化量推定手段と、前記変化量
推定手段によって推定した前記反射波のドップラー周波
数の変化量に応じて、前記信号列抽出手段により抽出さ
れた信号列の位相補償を行う位相補償手段と、前記位相
補償手段により位相補償された信号列をクロスレンジ圧
縮するクロスレンジ圧縮手段とを有するため、目標の移
動に伴い反射波のドップラー周波数が変化する場合であ
っても追尾誤差を押さえてドップラー追尾を行なうこと
ができる。
【図面の簡単な説明】
【図1】 この発明の実施の形態1におけるレンジ追尾
手段の構成図である。
【図2】 この発明の実施の形態1におけるピーク画像
生成手段の処理内容を説明するための原理図である。
【図3】 この発明の実施の形態1における画像線形補
正手段の処理内容を説明するための原理図である。
【図4】 この発明の実施の形態1における射影セル数
算出手段で得られる射影図である。
【図5】 この発明の実施の形態1における射影セル数
算出手段で算出した、各補正量に対する射影セル数を示
す図である。
【図6】 この発明の実施の形態2におけるドップラー
追尾手段の構成図である。
【図7】 この発明の実施の形態3におけるシフト量推
定回路の構成図である。
【図8】 この発明の実施の形態3における多点ピーク
画像生成手段で得られた多点ピーク画像である。
【図9】 この発明の実施の形態4におけるシフト量推
定回路の構成図である。
【図10】 この発明の実施の形態4における画像の各
セルの総和算出手段の処理内容を示す図である。
【図11】 この発明の実施の形態4における推定シフ
ト量算出手段でシフト量を算出するための評価指標の、
補正量に対する変化を示す図である。
【図12】 この発明の実施の形態5におけるシフト量
推定回路の構成図である。
【図13】 この発明の実施の形態6におけるレンジ補
償回路の構成図である。
【図14】 この発明の実施の形態6における画像非線
形補正手段の処理内容を示す図である。
【図15】 この発明の実施の形態7における位相補償
回路の構成図である。
【図16】 この発明の実施の形態7の処理内容を示す
フローチャートである。
【図17】 従来のレーダ装置のブロック構成図であ
る。
【図18】 従来のレーダ装置のレンジ補償を実現する
ための構成図である。
【図19】 従来のレーダ装置の位相補償を実現するた
めの構成図である。
【図20】 回転運動を行なう目標を観測するジオメト
リである。
【図21】 ISAR画像の一例である
【図22】 並進運動を行なう目標を観測するジオメト
リである。
【図23】 レンジプロフィールのヒストリの最大振幅
検出結果の一例である。
【図24】 レンジ補償処理を施した後のレンジプロフ
ィールのヒストリにおける最大振幅検出結果の一例であ
る。
【図25】 区分周波数分析手段の処理内容を説明する
図である。
【図26】 区分周波数分布のヒストリの最大振幅検出
結果の一例である。
【符号の説明】
1 送信機、2 送受切換器、3 送受信アンテナ、4
受信機、5 レンジ圧縮手段、6 レンジ補償回路、
7 位相補償回路、8 クロスレンジ圧縮手段、9 モ
ニタ・テレビ、111 ピーク画像生成手段、112
画像線形補正手段、113 射影セル数算出手段、11
4 推定シフト量算出手段、115 シフト量推定回
路、116 原画像、117 ピーク画像、118 補
正量0の時の補正後のピーク画像、119 補正量−4
の時の補正後のピーク画像、120補正量−8の時の補
正後のピーク画像、121 補正量−10の時の補正後
のピーク画像、122 補正量0の時の射影図、123
補正量−4の時の射影図、124 補正量−8の時の
射影図、125 補正量−10の時の射影図、131
多点ピーク画像生成手段、132 多点ピーク画像、1
41 画像の各セルの総和算出手段、142 エントロ
ピー算出手段、143 補正量0の場合のy方向の総和
を計算する過程の図、144 補正量─4の場合のy方
向の総和を計算する過程の図、145 補正量─8の場
合のy方向の総和を計算する過程の図、146 補正量
─10の場合のy方向の総和を計算する過程の図、16
1 画像非線形補正手段、162 推定シフト量非線形
算出手段、163 レンジ補償量非線形算出手段、16
4 非線形補正前のピーク画像、165 非線形補正後
のピーク画像、171 区分周波数分析点数更新手段、
172 位相補償量蓄積手段、501 レンジ追尾手
段、502 振幅最大レンジビン検出回路、503平滑
化回路、504 レンジ補償量算出手段、505 レン
ジ補償手段、506 注目レンジビン決定手段、507
区分周波数分析手段、508 ドップラー追尾手段、
509 振幅最大周波数検出回路、510 位相補償量
算出手段、511 位相補償手段、520 目標、52
1 レーダ装置、522 目標のISAR画像、530
注目レンジビンの受信信号列、531 区分周波数分
布のヒストリ。
───────────────────────────────────────────────────── フロントページの続き (72)発明者 臼井 隆三郎 東京都千代田区丸の内二丁目2番3号 三 菱電機株式会社内

Claims (13)

    【特許請求の範囲】
  1. 【請求項1】 移動する目標に対して電波を送信し、前
    記目標からの反射波を受信して前記目標の画像を得るた
    め、前記目標の移動に伴う前記目標のレンジの変化を補
    償する補償回路であって、 前記反射波の受信信号列を基に、移動に伴う前記目標の
    レンジの変化を示す原画像を作成する原画像作成手段
    と、 複数の補正量によって前記原画像を補正することにより
    複数の補正画像を作成する補正手段と、 前記補正手段によって作成された複数の補正画像を評価
    するための評価値を算出する評価値算出手段と、 前記評価値算出手段によって前記複数の補正画像に対応
    して算出された評価値に基づき前記原画像を補正するの
    に適した補正量を決定し、この決定した補正量により前
    記目標の移動に伴う前記目標のレンジの変化量を推定す
    る変化量推定手段と、 前記変化量推定手段によって推定した前記目標のレンジ
    の変化量に応じて前記受信信号列のレンジ補償を行うレ
    ンジ補償手段とを有することを特徴とする補償回路。
  2. 【請求項2】 移動する目標に対して電波を送信し、前
    記目標からの反射波を受信して前記目標の画像を得るた
    め、前記目標の移動に伴う前記反射波の位相の変化を補
    償する補償回路であって、 前記反射波の受信信号列を基に、前記目標の移動に伴う
    前記反射波のドップラー周波数の変化を示す原画像を作
    成する原画像作成手段と、 前記原画像を複数の補正量によって補正することにより
    複数の補正画像を作成する補正手段と、 前記補正手段によって作成された複数の補正画像を評価
    するための評価値を算出する評価値算出手段と、 前記評価値算出手段によって前記複数の補正画像に対応
    して算出された評価値に基づき前記原画像を補正するの
    に適した補正量を決定し、この決定した補正量により前
    記目標の移動に伴う前記反射波のドップラー周波数の変
    化量を推定する変化量推定手段と、 前記変化量推定手段によって推定された前記反射波のド
    ップラー周波数の変化量に応じて、前記受信信号列の位
    相補償を行う位相補償手段とを有することを特徴とする
    補償回路。
  3. 【請求項3】 前記評価値算出手段は、前記補正画像の
    各レンジ毎に、同一のレンジに存在している複数の画素
    の値を積分し、当該積分値の分布を示す指標を評価値と
    して算出することを特徴とする請求項1記載の補償回
    路。
  4. 【請求項4】 前記評価値算出手段は、前記補正画像の
    各ドップラー周波数毎に、同一のドップラー周波数に対
    して存在している複数の画素の値を積分し、当該積分値
    の分布を示す指標を評価値として算出することを特徴と
    する請求項2記載の補償回路。
  5. 【請求項5】 移動する目標に対して電波を送信し、前
    記目標からの反射波を受信して前記目標の画像を得るた
    め、前記目標の移動に伴う前記目標のレンジの変化を補
    償する補償回路であって、 前記反射波の受信信号列を基に、移動に伴う前記目標の
    レンジの変化を示す原画像を作成する原画像作成手段
    と、 前記原画像作成手段によって作成された原画像を二値化
    することによりニ値化画像を作成する画像二値化手段
    と、 前記二値化画像を複数の補正量によって補正することに
    より複数の補正ニ値化画像を作成する補正手段と、 前記補正手段によって作成された複数の補正二値化画像
    を射影することにより複数の射影画像を作成し、この複
    数の射影画像に対する評価値を算出する評価値算出手段
    と、 前記評価値算出手段によって前記複数の射影画像に対応
    して算出された評価値に基づき前記原画像を補正するの
    に適した補正量を決定し、この決定した補正量により前
    記目標の移動に伴う前記目標のレンジの変化量を推定す
    る変化量推定手段と、 前記変化量推定手段によって推定した前記目標のレンジ
    の変化量に応じて、前記受信信号列のレンジ補償を行う
    レンジ補償手段とを有することを特徴とする補償回路。
  6. 【請求項6】 移動する目標に対して電波を送信し、前
    記目標からの反射波を受信して前記目標の画像を得るた
    め、前記目標の移動に伴う前記反射波の位相の変化を補
    償する補償回路であって、 前記反射波の受信信号列を基に、移動に伴う前記反射波
    のドップラー周波数の変化を示す原画像を作成する原画
    像作成手段と、 前記原画像作成手段によって作成された原画像を二値化
    することによりニ値化画像を作成する画像二値化手段
    と、 前記二値化画像を複数の補正量によって補正することに
    より複数の補正ニ値化画像を作成する補正手段と、 前記複数の補正二値化画像を射影することにより複数の
    射影画像を作成し、この複数の射影画像に対する評価値
    を算出する評価値算出手段と、 前記評価値算出手段によって前記複数の射影画像に対応
    して算出された評価値に基づいて前記原画像を補正する
    のに適した補正量を決定し、この決定した補正量により
    前記目標の移動に伴う前記目標のレンジの変化量を推定
    する変化量推定手段と、 前記変化量推定手段によって推定した前記反射波のドッ
    プラー周波数の変化量に応じて、前記受信信号列の位相
    補償を行う位相補償手段とを有することを特徴とする補
    償回路。
  7. 【請求項7】 前記評価値算出手段は、前記射影画像内
    の特定の値を有する画素の画素数を評価値として算出す
    ることを特徴とする請求項5又は請求項6記載の補償回
    路。
  8. 【請求項8】 前記補正手段は、前記複数の補正量に応
    じて前記原画像を線形補正することを特徴とする請求項
    1乃至請求項7のいずれかに記載の補償回路。
  9. 【請求項9】 前記補正手段は、前記複数の補正量に応
    じて前記原画像を非線形補正することを特徴とする請求
    項1乃至請求項7のいずれかに記載の補償回路。
  10. 【請求項10】 移動する目標に対して電波を送信し、
    前記目標からの反射波を受信して前記目標の画像を得る
    ため、前記目標の移動に伴う前記目標のレンジの変化を
    補償する補償方法であって、 前記反射波の受信信号列を基に、移動に伴う前記目標の
    レンジの変化を示す原画像を作成する原画像作成ステッ
    プと、 前記原画像を複数の補正量によって補正することにより
    複数の補正画像を作成する補正ステップと、 前記補正ステップにおいて作成された複数の補正画像を
    評価するための評価値を算出する評価値算出ステップ
    と、 前記評価値算出ステップにおいて前記複数の補正画像に
    対応して算出された評価値に基づき前記原画像を補正す
    るのに適した補正量を決定し、この決定した補正量によ
    り前記目標の移動に伴う前記目標のレンジの変化量を推
    定する変化量推定ステップと、 前記変化量推定ステップにおいて推定した前記目標のレ
    ンジの変化量に応じて、前記受信信号列のレンジ補償を
    行うレンジ補償ステップとを有することを特徴とする補
    償方法。
  11. 【請求項11】 移動する目標に対して電波を送信し、
    前記目標からの反射波を受信して前記目標の画像を得る
    ため、前記目標の移動に伴う前記反射波の位相の変化を
    補償する補償方法であって、 前記反射波の受信信号列を基に、前記目標の移動に伴う
    前記反射波のドップラー周波数の変化を示す原画像を作
    成する原画像作成ステップと、 前記原画像を複数の補正量によって補正することにより
    複数の補正画像を作成する補正ステップと、 前記補正ステップにおいて作成された複数の補正画像を
    評価するための評価値を算出する評価値算出ステップ
    と、 前記評価値算出ステップにおいて前記複数の補正画像に
    対応して算出された評価値に基づき前記原画像を補正す
    るのに適した補正量を決定し、この決定した補正量によ
    り前記目標の移動に伴う前記反射波のドップラー周波数
    の変化量を推定する変化量推定ステップと、 前記変化量推定ステップにおいて推定した前記反射波の
    ドップラー周波数の変化量に応じて、前記受信信号列の
    位相補償を行う位相補償ステップとを有することを特徴
    とする補償方法。
  12. 【請求項12】 移動する目標に対して電波を送信する
    送信手段と、 前記目標からの反射波を受信する受信手段と、 前記受信手段により受信された反射波から得られた受信
    信号列をレンジ圧縮するレンジ圧縮手段と、 前記レンジ圧縮手段においてレンジ圧縮された受信信号
    列を基に、前記目標の移動に伴う前記目標のレンジの変
    化を示す原画像を作成する原画像作成手段と、 前記原画像を複数の補正量によって補正することにより
    複数の補正画像を作成する補正手段と、 前記補正手段によって作成された複数の補正画像を評価
    するための評価値を算出する評価値算出手段と、 前記評価値算出手段によって前記複数の補正画像に対応
    して算出された評価値に基づき前記原画像を補正するの
    に適した補正量を決定し、前記決定した補正量により前
    記目標の移動に伴う前記目標のレンジの変化量を推定す
    る変化量推定手段と、 前記変化量推定手段によって推定された前記目標のレン
    ジの変化量に応じて、前記レンジ圧縮手段によってレン
    ジ圧縮された受信信号列のレンジ補償を行うレンジ補償
    手段と、 前記レンジ補償手段によりレンジ補償がなされた信号列
    の中から、注目するレンジビンの信号列を抽出する信号
    列抽出手段と、 前記信号列抽出手段において抽出された信号列をクロス
    レンジ圧縮するクロスレンジ圧縮手段とを有するレーダ
    装置。
  13. 【請求項13】 移動する目標に対して電波を送信する
    送信手段と、 前記目標からの反射波を受信する受信手段と、 前記受信手段により受信された反射波から得られた受信
    信号列をレンジ圧縮するレンジ圧縮手段と、 前記レンジ圧縮手段によりレンジ圧縮された受信信号列
    の内、注目するレンジビンの信号列を抽出する信号列抽
    出手段と、 前記信号列抽出ステップにより抽出された信号列に基づ
    き、前記目標の移動に伴う前記反射波のドップラー周波
    数の変化を示す原画像を作成する原画像作成手段と、 前記原画像を複数の補正量によって補正することにより
    複数の補正画像を作成する補正手段と、 前記補正手段によって作成された複数の補正画像を評価
    するための評価値を算出する評価値算出手段と、 前記評価値算出手段によって前記複数の補正画像に対応
    して算出された評価値に基づいて前記原画像を補正する
    のに適した補正量を決定し、この決定した補正量により
    前記目標の移動に伴う前記反射波のドップラー周波数の
    変化量を推定する変化量推定手段と、 前記変化量推定手段によって推定した前記反射波のドッ
    プラー周波数の変化量に応じて、前記信号列抽出手段に
    より抽出された信号列の位相補償を行う位相補償手段
    と、 前記位相補償手段により位相補償された信号列をクロス
    レンジ圧縮するクロスレンジ圧縮手段とを有するレーダ
    装置。
JP07507397A 1997-03-27 1997-03-27 補償回路、補償方法およびレーダ装置 Expired - Fee Related JP3642149B2 (ja)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP07507397A JP3642149B2 (ja) 1997-03-27 1997-03-27 補償回路、補償方法およびレーダ装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP07507397A JP3642149B2 (ja) 1997-03-27 1997-03-27 補償回路、補償方法およびレーダ装置

Publications (2)

Publication Number Publication Date
JPH10268042A true JPH10268042A (ja) 1998-10-09
JP3642149B2 JP3642149B2 (ja) 2005-04-27

Family

ID=13565660

Family Applications (1)

Application Number Title Priority Date Filing Date
JP07507397A Expired - Fee Related JP3642149B2 (ja) 1997-03-27 1997-03-27 補償回路、補償方法およびレーダ装置

Country Status (1)

Country Link
JP (1) JP3642149B2 (ja)

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2000329849A (ja) * 1999-05-17 2000-11-30 Mitsubishi Electric Corp レーダ信号処理器
JP2006343290A (ja) * 2005-06-10 2006-12-21 Mitsubishi Electric Corp 画像レーダ装置
JP2011174866A (ja) * 2010-02-25 2011-09-08 Toshiba Corp 目標検出・追尾装置
JP2013130472A (ja) * 2011-12-21 2013-07-04 Mitsubishi Electric Corp 画像レーダ装置
JP2013178206A (ja) * 2012-02-29 2013-09-09 Tokyo Keiki Inc 目標運動予測装置及び目標運動予測方法

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112255626B (zh) * 2020-10-10 2023-09-19 中国人民解放军63921部队 基于参数化平动补偿的isar反投影方法

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2000329849A (ja) * 1999-05-17 2000-11-30 Mitsubishi Electric Corp レーダ信号処理器
JP2006343290A (ja) * 2005-06-10 2006-12-21 Mitsubishi Electric Corp 画像レーダ装置
JP2011174866A (ja) * 2010-02-25 2011-09-08 Toshiba Corp 目標検出・追尾装置
JP2013130472A (ja) * 2011-12-21 2013-07-04 Mitsubishi Electric Corp 画像レーダ装置
JP2013178206A (ja) * 2012-02-29 2013-09-09 Tokyo Keiki Inc 目標運動予測装置及び目標運動予測方法

Also Published As

Publication number Publication date
JP3642149B2 (ja) 2005-04-27

Similar Documents

Publication Publication Date Title
US6952178B2 (en) Method of detecting moving objects and estimating their velocity and position in SAR images
Zhou et al. Approach for single channel SAR ground moving target imaging and motion parameter estimation
US7777665B1 (en) Comparing range data across the slow-time dimension to correct motion measurement errors beyond the range resolution of a synthetic aperture radar
JP4046422B2 (ja) 運動補償回路及びレーダ装置
CN116299551B (zh) 一种太赫兹sar二维自聚焦成像算法
US9329264B2 (en) SAR image formation
JP3395606B2 (ja) 位相補償回路、位相補償方法およびレーダ装置
US7551119B1 (en) Flight path-driven mitigation of wavefront curvature effects in SAR images
CN109507669A (zh) 基于相参累积的地面微弱运动目标参数估计方法
CN112255626B (zh) 基于参数化平动补偿的isar反投影方法
JP3642149B2 (ja) 補償回路、補償方法およびレーダ装置
JP3360562B2 (ja) 補償回路、補償方法およびレーダ装置
Zhan et al. An adaptive nonlinear phase error estimation and compensation method for terahertz radar imaging system
JP2003344532A (ja) レーダ装置
JP3946860B2 (ja) レーダ装置
JP2001141821A (ja) レーダ信号処理器
JP3949489B2 (ja) レーダ装置
CN113534143B (zh) 基于地球同步轨道合成孔径雷达的高程空变误差补偿方法
JP2008139206A (ja) 高分解能レーダ装置
JP3709701B2 (ja) レーダ信号処理器
CN115825899A (zh) 基于虚拟孔径的分布式通道自校准方法
JP3832139B2 (ja) レーダ信号処理器
JP2000187074A (ja) レーダ信号処理器
CN118537272B (zh) 合成孔径雷达图像的残余距离徙动矫正方法、装置及设备
KR102874424B1 (ko) Sar 영상 처리 장치 및 방법

Legal Events

Date Code Title Description
RD01 Notification of change of attorney

Free format text: JAPANESE INTERMEDIATE CODE: A7421

Effective date: 20040716

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20041101

TRDD Decision of grant or rejection written
A01 Written decision to grant a patent or to grant a registration (utility model)

Free format text: JAPANESE INTERMEDIATE CODE: A01

Effective date: 20050105

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20050118

LAPS Cancellation because of no payment of annual fees