JP6261839B1 - 合成開口レーダ信号処理装置 - Google Patents
合成開口レーダ信号処理装置 Download PDFInfo
- Publication number
- JP6261839B1 JP6261839B1 JP2017549349A JP2017549349A JP6261839B1 JP 6261839 B1 JP6261839 B1 JP 6261839B1 JP 2017549349 A JP2017549349 A JP 2017549349A JP 2017549349 A JP2017549349 A JP 2017549349A JP 6261839 B1 JP6261839 B1 JP 6261839B1
- Authority
- JP
- Japan
- Prior art keywords
- curvature
- point
- interpolation processing
- processing unit
- synthetic aperture
- 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.)
- Expired - Fee Related
Links
- 238000004364 calculation method Methods 0.000 claims abstract description 35
- 238000000034 method Methods 0.000 abstract description 80
- 230000008569 process Effects 0.000 abstract description 53
- 230000006835 compression Effects 0.000 description 25
- 238000007906 compression Methods 0.000 description 25
- 238000010586 diagram Methods 0.000 description 12
- 238000004422 calculation algorithm Methods 0.000 description 11
- 238000013500 data storage Methods 0.000 description 8
- 230000006870 function Effects 0.000 description 7
- 238000003672 processing method Methods 0.000 description 4
- 239000004065 semiconductor Substances 0.000 description 4
- 230000005540 biological transmission Effects 0.000 description 3
- 238000006243 chemical reaction Methods 0.000 description 3
- 230000003247 decreasing effect Effects 0.000 description 3
- 230000010354 integration Effects 0.000 description 3
- 230000008859 change Effects 0.000 description 2
- 238000007796 conventional method Methods 0.000 description 2
- 238000010606 normalization Methods 0.000 description 2
- 230000009466 transformation Effects 0.000 description 2
- 230000000694 effects Effects 0.000 description 1
- 238000003384 imaging method Methods 0.000 description 1
- 230000001678 irradiating effect Effects 0.000 description 1
- 238000011158 quantitative evaluation Methods 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 238000010187 selection method Methods 0.000 description 1
- 238000004904 shortening Methods 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 230000001131 transforming effect Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
- G01S13/88—Radar or analogous systems specially adapted for specific applications
- G01S13/89—Radar or analogous systems specially adapted for specific applications for mapping or imaging
- G01S13/90—Radar or analogous systems specially adapted for specific applications for mapping or imaging using synthetic aperture techniques, e.g. synthetic aperture radar [SAR] techniques
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
- G01S13/88—Radar or analogous systems specially adapted for specific applications
- G01S13/89—Radar or analogous systems specially adapted for specific applications for mapping or imaging
- G01S13/90—Radar or analogous systems specially adapted for specific applications for mapping or imaging using synthetic aperture techniques, e.g. synthetic aperture radar [SAR] techniques
- G01S13/9021—SAR image post-processing techniques
- G01S13/9027—Pattern recognition for feature extraction
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
- G01S13/88—Radar or analogous systems specially adapted for specific applications
- G01S13/89—Radar or analogous systems specially adapted for specific applications for mapping or imaging
- G01S13/90—Radar or analogous systems specially adapted for specific applications for mapping or imaging using synthetic aperture techniques, e.g. synthetic aperture radar [SAR] techniques
- G01S13/9021—SAR image post-processing techniques
- G01S13/9029—SAR image post-processing techniques specially adapted for moving target detection within a single SAR image or within multiple SAR images taken at the same time
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
- G01S13/88—Radar or analogous systems specially adapted for specific applications
- G01S13/89—Radar or analogous systems specially adapted for specific applications for mapping or imaging
- G01S13/90—Radar or analogous systems specially adapted for specific applications for mapping or imaging using synthetic aperture techniques, e.g. synthetic aperture radar [SAR] techniques
- G01S13/9094—Theoretical aspects
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Radar, Positioning & Navigation (AREA)
- Physics & Mathematics (AREA)
- Electromagnetism (AREA)
- Computer Networks & Wireless Communication (AREA)
- General Physics & Mathematics (AREA)
- Artificial Intelligence (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
本発明に係る合成開口レーダ信号処理装置は、合成開口レーダの受信信号から得られる離散データを補間処理する低精度補間処理部と、低精度補間処理部より高い精度で離散データを補間処理する高精度補間処理部と、補間処理の対象領域における離散データの曲率である第1の曲率に応じて高精度補間処理部または低精度補間処理部のいずれか一方を選択する曲率判定部と、曲率判定部で選択した補間処理の結果を用いて画像を再生する画像再生処理部とを備え、曲率判定部は、対象領域の近傍で所定方向に隣接して並ぶ離散データA点、B点、C点を選択し、A点とB点とを結ぶ線分を(AB+BC):BCに外分する点とC点との間の距離Δf1Rに基づいて第1の曲率を決定する合成開口レーダ信号処理装置である(ただし、ABはA点とB点との間の距離、BCはB点とC点との間の距離。)。これにより、離散データから曲率を計算する際の演算量を抑えることができる。
Description
この発明は、例えば、航空機や人工衛星などの移動プラットフォームに搭載されて、地表や海面などの高分解能画像である合成開口レーダ(SAR:Synthetic Aperture Radar)画像を再生する合成開口レーダ信号処理装置に関するものである。
航空機や人工衛星などの移動プラットフォームに搭載される合成開口レーダ信号処理装置は、アンテナを有するSARセンサを備える。SARセンサは、移動プラットフォームの移動中に電磁波を繰り返し放射しながら、目標に反射された当該電磁波の反射エコーを受信する。合成開口レーダ信号処理装置は、SARセンサの受信信号に対する信号処理を実施することで、2次元のSAR画像を得る。SARセンサの受信信号から2次元のSAR画像を得る処理は画像再生処理と呼ばれ、その画像再生処理の処理方法は画像再生アルゴリズムと呼ばれる。画像再生アルゴリズムの一例として、以下の非特許文献1に開示されているバックプロジェクションアルゴリズムがある。
このバックプロジェクションアルゴリズムは、プラットフォーム搭載のSARセンサと、電磁波が反射される目標点との間の距離を計算し、この距離に相当する信号を抽出して積分することでSAR画像を得る方法である。
非特許文献1には、上記の距離計算と積分処理の演算量を削減するために、極座標形式グリッドを使用して信号を積分処理し、2次元補間処理によってデカルト座標形式にすることで最終的なSAR画像を得る手法に関する記載がある。
2次元補間処理において、演算量が少ない補間手法(最近傍法、線形法など)を使用すると補間精度が低くなる。一方で、補間精度が高い補間手法(多項式補間)を用いると演算量が多くなる。このように、補間精度と演算量にはトレードオフの関係がある。
そこで、補間の対象となる離散データの区間毎に低精度、高精度の補間処理手法を切り替えることによって、補間精度を保ちつつ、かつ補間処理の演算量を抑えられるような処理が実現できると考えられる。補間または近似処理の切り替え手法として、例えば、次に示す特許文献1の手法がある。
特許文献1には、一定の時間毎に逐次的に発生する離散データを要約して情報量を削減し記録する手法について記載されている。この手法は、発生した全ての離散データをRAMやハードウェアなどの保存手段に格納するのではなく、離散データ点群を精度良く近似式で表し、近似式の係数部分を保存手段に格納することにより、情報量を削減する技術である。
L. M. H. Ulander, H. Hellsten, and G. Stenstrom, "Synthetic Aperture Radar Processing Using Fast Factorized Back−Projection," IEEE Transactions on Aerospace and electronic Systems, vol. 39, no. 3, pp.760−776, 2003.
特許文献1には、離散データの曲率の大きさに応じて近似式を切り替える手法について記載されている。しかしながら、特許文献1に記載されている曲率の計算手法は演算量が多く、仮に補間処理の演算量と補間精度の双方にとって最適な補間手法が選択できたとしても、補間の切り替えを判定する際の演算量がボトルネックとなる。その結果、補間処理手法の切り替えによる補間処理時間の短縮効果が小さくなるという課題があった。
この発明は、上記の課題を解決するためになされたもので、離散データの曲率に応じて補間処理を切り替えるSAR画像の再生処理において、曲率計算の演算量を削減することができる合成開口レーダ信号処理装置を得ることを目的とする。
本発明に係る合成開口レーダ信号処理装置は、合成開口レーダの受信信号から得られる離散データを補間処理する低精度補間処理部と、低精度補間処理部より高い精度で離散データを補間処理する高精度補間処理部と、補間処理の対象領域における離散データの曲率である第1の曲率に応じて高精度補間処理部または低精度補間処理部のいずれか一方を選択する曲率判定部と、曲率判定部で選択した補間処理の結果を用いて画像を再生する画像再生処理部とを備え、曲率判定部は、対象領域の近傍で所定方向に隣接して並ぶ離散データA点、B点、C点を選択し、A点とB点とを結ぶ線分を(AB+BC):BCに外分する点とC点との間の距離Δf1Rに基づいて第1の曲率を決定する合成開口レーダ信号処理装置である(ただし、ABはA点とB点との間の距離、BCはB点とC点との間の距離。)。
この発明によれば、記のように構成したことにより、離散データから曲率を計算する際の演算量を抑えることができる。
本発明は、合成開口レーダの受信信号から得られる離散データを補間処理してSAR画像を再生する際に、補間対象領域の離散データの曲率に基づいて補間処理手法を選択する合成開口レーダ信号処理装置に関するものであり、特に、少ない演算量で曲率を計算することを特徴とするものである。以下、この発明をより詳細に説明するために、この発明を実施するための形態について、添付の図面にしたがって説明する。
実施の形態1.
本実施の形態における合成開口レーダ信号処理装置100の構成について説明する。図1は、この発明の実施の形態1による合成開口レーダ信号処理装置100の構成の一例を示す機能構成図である。合成開口レーダ信号処理装置100は、受信データ格納部2から読み出したSARセンサ1の受信データを処理してSAR画像を再生する。また、合成開口レーダ信号処理装置100は、画像再生処理部3、曲率判定部4、低精度補間処理部5、高精度補間処理部6、SAR画像格納部7から構成される。なお、本実施の形態では、合成開口レーダ信号処理装置100の外部にSARセンサ1と受信データ格納部2とを設けた構成について説明するが、これに限定するものではなく、合成開口レーダ信号処理装置100の内部にSARセンサ1と受信データ格納部2とを備える構成であってもよい。また、本実施の形態では、合成開口レーダ信号処理装置100の内部にSAR画像格納部7を設けた構成について説明するが、合成開口レーダ信号処理装置100の外部に備える構成であってもよい。
本実施の形態における合成開口レーダ信号処理装置100の構成について説明する。図1は、この発明の実施の形態1による合成開口レーダ信号処理装置100の構成の一例を示す機能構成図である。合成開口レーダ信号処理装置100は、受信データ格納部2から読み出したSARセンサ1の受信データを処理してSAR画像を再生する。また、合成開口レーダ信号処理装置100は、画像再生処理部3、曲率判定部4、低精度補間処理部5、高精度補間処理部6、SAR画像格納部7から構成される。なお、本実施の形態では、合成開口レーダ信号処理装置100の外部にSARセンサ1と受信データ格納部2とを設けた構成について説明するが、これに限定するものではなく、合成開口レーダ信号処理装置100の内部にSARセンサ1と受信データ格納部2とを備える構成であってもよい。また、本実施の形態では、合成開口レーダ信号処理装置100の内部にSAR画像格納部7を設けた構成について説明するが、合成開口レーダ信号処理装置100の外部に備える構成であってもよい。
SARセンサ1は、例えば、アンテナ、送信機および受信機(いずれも図示せず)などから構成されており、航空機や人工衛星などの移動プラットフォームに搭載される。SARセンサ1は、送信機により生成された高周波パルス信号をアンテナから空間に放射する一方、観測目標によって反射された高周波パルス信号のエコー信号をアンテナで受信する。受信機は、アンテナで受信した受信信号を増幅し、増幅した受信信号の周波数を中間周波数に変換した後、デジタルの受信データに変換して出力する。なお、SARセンサ1は、合成開口レーダとも言う。
受信データ格納部2は、例えば、RAM(Random Access Memory)やハードディスクなどの記録装置であり、SARセンサ1で取得した受信データを格納する。
画像再生処理部3は、例えば、CPU(Central Processing Unit)を実装している半導体集積回路、あるいはワンチップマイコンなどから構成される。画像再生処理部3は、受信データ格納部2に格納される受信データに対して画像再生処理を実施することでSAR画像を再生する。受信データは、複数の目標点からの様々な反射信号が重なりあっているため、圧縮処理することで画像化し、SAR画像を得ることができる。詳細な画像再生処理については後述する。
曲率判定部4は、例えば、CPUを実装している半導体集積回路、あるいはワンチップマイコンなどから構成される。曲率判定部4は、画像再生処理部3より得られるアジマス圧縮後データを用いて曲率を判定し、曲率の判定結果に基づいて、低精度補間処理部5または高精度補間処理部6のいずれか一方を選択する。画像再生処理部3は、曲率判定部4が選択した補間処理部の補間処理結果を用いてSAR画像を再生する。詳細な曲率の判別手法や補間処理の選択手法については後述する。
低精度補間処理部5は、例えば、CPUを実装している半導体集積回路、あるいはワンチップマイコンなどから構成される。曲率判定部4が低精度補間処理部5を選択した場合には、画像再生処理部3より得られるアジマス圧縮後データに対して低精度補間処理を行い、補間処理結果を画像再生処理部3へ渡す。
高精度補間処理部6は、例えば、CPUを実装している半導体集積回路、あるいは、ワンチップマイコンなどから構成される。曲率判定部4が高精度補間処理部6を選択した場合には、画像再生処理部3より得られるアジマス圧縮後データに対して高精度補間処理を行い、補間処理結果を画像再生処理部3へ渡す。
SAR画像格納部7は、例えば、RAMやハードディスクなどの記録装置であり、画像再生処理部3で再生されたSAR画像を格納する。
以上の説明では、合成開口レーダ信号処理装置の構成要素である画像再生処理部3、曲率判定部4、低精度補間処理部5、高精度補間処理部6、およびSAR画像格納部7は、それぞれ専用のハードウェアで実現される構成を想定しているが、これに限定するものではなく、コンピュータで実現する構成でもよい。
図2は、合成開口レーダ信号処理装置100がコンピュータで構成される場合のハードウェア構成図である。合成開口レーダ信号処理装置100がコンピュータで構成される場合、画像再生処理部3、曲率判定部4、低精度補間処理部5、高精度補間処理部6の処理内容を記述した合成開口レーダ信号処理プログラムがコンピュータのメモリ11に格納される。当該コンピュータのCPUなどのプロセッサ12は、メモリ11に格納されている合成開口レーダ信号プログラムを実行することにより、各処理を実現する。SAR画像格納部7は、ハードディスクなどの記録装置13で実現しても良いし、合成開口レーダ信号処理装置100からデータの読み書きが可能な外部の記録装置で実現しても良い。
次に動作について説明する。図3は、本実施の形態における合成開口レーダ信号処理装置100の動作を示すフローチャートであり、SARセンサ1の受信データからSAR画像を生成する一連の処理を示している。
SARセンサ1は、送信機により生成された高周波パルス信号をアンテナから空間に放射する。アンテナから空間に放射された高周波パルス信号が観測目標で反射され、その高周波パルス信号の反射であるエコー信号が戻ってくると、SARセンサ1のアンテナは、当該エコー信号を受信する。また、SARセンサ1の受信機は、アンテナの受信信号を増幅し、増幅した受信信号の周波数を中間周波数に変換した後、中間周波数の受信信号をA/D変換することでデジタルの受信データを生成し、受信データ格納部2へ出力する。
画像再生処理部3は、受信データ格納部2から受信データを読み出して入力する(ステップST21)。画像再生処理部3は、受信データのレンジ方向の分解能を高めるため、受信データをレンジ圧縮する(ステップST22)。受信データのレンジ圧縮処理は、受信データとレンジ参照関数とをそれぞれレンジ方向にフーリエ変換し、乗算することでレンジ方向の高分解能化を実現する処理手法である。なお、SARセンサ1を搭載するプラットフォームの進行方向をアジマス方向、SARセンサ1から発せられる電波ビームの進行方向をレンジ方向という。レンジ方向に圧縮処理をするのがレンジ圧縮、アジマス方向に圧縮処理をするのがアジマス圧縮である。
さらに、画像再生処理部3は、アジマス方向の分解能を高めるため、レンジ圧縮後のデータをアジマス圧縮する(ステップST23)。アジマス圧縮後のデータは複素の信号であり、これを絶対値で表したものがSAR画像となる。アジマス圧縮処理の原理はレンジ圧縮と同様であるが、使用する画像再生アルゴリズム方式によって具体的なアジマス圧縮の処理内容は異なる。例えば、周波数空間でレンジ圧縮後のデータとアジマス参照関数とを乗算することにより、アジマス圧縮処理をするレンジドップラー型アルゴリズムや、周波数領域でチャープ信号のチャープ率に変調をかけることで、スケーリングとフーリエ変換のみで画像化を行うチャープスケーリングアルゴリズムや、2次元の周波数領域に変換したデータに対し、レンジ方向の角周波数を極座標の波数に変換するω−kアルゴリズムや、上述した時間空間で処理をするバックプロジェクション型アルゴリズムなどがある。これらのアルゴリズムは公知の技術であるため、詳細な説明は省略する。
次に、画像再生処理部3は、ステップST24からステップST33までのループ処理において、アジマス圧縮後データに対して補間処理を行い、最終的なSAR画像を生成する。画像再生処理部3が行う補間処理は、非特許文献1で行う2次元補間処理と同様、座標変換を目的とするものである。
補間処理による座標変換について、図を参照して説明する。図4は、2次元補間処理による座標変換を模式的に示す説明図である。
バックプロジェクションアルゴリズムでは、プラットフォームに搭載したSARセンサ1とSARセンサ1が放射する電磁波を反射する目標点との間の距離と等距離の信号データを受信データから抽出し、抽出した信号を目標点の位置で積分処理をすることでSAR画像を得る。
これに対して、ファストバックプロジェクションでは、図4に示すように、極座標形式グリッドを使用して信号を積分処理した後、2次元補間処理によって極座標形式からデカルト座標形式に変換することで最終的なSAR画像を得る。このような手法を取り入れることで、距離計算と積分処理の演算量を削減することができる。
補間処理の対象単位について説明する。ステップST24からステップST33のループ処理で行われる補間処理の対象単位は、アジマス圧縮後のデータ点間を結ぶ線分を辺とする領域とし、この領域を対象領域と呼ぶ。画像再生処理部3は、対象領域ごとに、低精度補間処理部5あるいは高精度補間処理部6のいずれか一方の処理結果を用いて最終的なSAR画像を生成する処理を繰り返し実行する。つまり、対象領域の総数分のループ処理が行われる。補間処理の選択は、曲率判定部4が行う。曲率判定部4は、対象領域ごとにアジマス圧縮後データの曲率を判定し、曲率に応じて補間処理を選択する。
補間処理の選択方法について、図を参照して説明する。図5は、アジマス圧縮後データの補間処理の対象領域と、補間処理の種類分けを示す図である。アジマス圧縮後データは、離散データであり、図のi軸方向とj軸方向とで構成される平面上で定義される。f(i,j)軸は、位置(i,j)におけるアジマス圧縮後データの信号強度を表す。ここで、i軸方向やj軸方向は任意であり、座標系を問わないため、例えば、極座標形式で表しても構わない。アジマス圧縮後データは、図中に丸印で示すように、i軸方向、j軸方向にそれぞれ等間隔に配置される。離散データであるアジマス圧縮後データをつなぐ破線は、アジマス圧縮後データの真の値を示す。
補間処理の対象領域は、隣接するアジマス圧縮後データの4点を結ぶ線分L1,L2,L3,L4を辺とする領域(図中の太い実線で囲まれる部分)、つまり、隣接する4つのアジマス圧縮後データを頂点とする格子領域である。
曲率判定部4は、対象領域ごとにアジマス圧縮後データの曲率を判定し、曲率に応じて補間処理を選択する。曲率を判定する際には、まず、対象領域を構成する辺ごとに曲率を判定し、その結果を総合して対象領域の曲率とする。図に示す例では、一点鎖線の楕円A1で示す範囲に含まれる領域は曲率が大きく、一点鎖線の楕円A2で示す範囲に含まれる領域は曲率が小さい。曲率が大きい領域では、精度の高い補間処理が選択され、曲率が小さい領域では、精度の低い補間処理が選択される。具体的には、対象領域が曲面状の場合、つまり信号強度の変化が大きい場合には、線形補間のような低精度補間処理では正確に曲面を表現することが難しいため、多項式を用いた高精度補間処理が選択される。一方で、対象領域が平坦な場合、つまり信号強度の変化が小さい場合には、線形補間のような低精度補間処理が選択される。よって、曲率判定部4は、曲率が大きい領域では高精度補間処理部6を選択し、曲率が小さい領域では低精度補間処理部5を選択する。
再び、フローチャートの説明に戻る。ステップST24からステップST33までのループ処理は、対象領域ごとに繰り返し行われる。よって、画像再生処理部3は、アジマス圧縮後データで構成される対象領域の中から、補間処理を行っていない対象領域を選択して曲率判定部4へ伝える。
曲率判定部4は、画像再生処理部3で選択した対象領域におけるアジマス圧縮後データの曲率を判定する。まず、曲率判定部4は、対象領域を構成する辺の中から、曲率の判定処理を行っていない辺を1つ選択する(ステップST25)。
次に、曲率判定部4は、選択した辺の両端にある2つの信号位置において、アジマス圧縮後データの信号強度の曲率を算出する(ステップST26)。一方の信号位置における曲率を第1の曲率、他方の信号位置における曲率を第2の曲率とする。曲率の算出手法については後述する。
その後、曲率判定部4は、第1および第2の曲率の絶対値を基準値(閾値とも言う)と比較する(ステップS27)。その結果、少なくとも一方の曲率の絶対値が基準値以上となる場合には、曲率判定部4は、対象領域の曲率が大きいと判定して高精度補間処理部6を選択する(ステップS29)。また、両方の曲率の絶対値が基準値未満となる場合には、曲率判定部4は、対象領域を構成する全ての辺に対してステップST27の判定処理を実施したか否かを確認する(ステップST28)。全ての辺で実施済みの場合、つまり、曲率が基準値以上となる辺が存在しない場合には、曲率判定部4は、対象領域の曲率が小さいと判定して低精度補間処理部5を選択する(ステップST30)。ステップST27の判定処理を実施していない辺がある場合には、ステップST25へ戻る。曲率判定部4は、対象領域を構成する辺の中から、曲率の判定処理を行っていない辺を選択して処理を続ける。
低精度補間処理部5が行う低精度補間処理の例として、線形補間について説明する。既知の4点が与えられたとき、その内部点(x,y)の値zを2次元補間処理で計算する場合を考える。図6は、2次元補間処理の補間座標点を示す説明図である。4点(x1,y1)、(x1,y2)、(x2,y1)、(x2,y2)は既知の点の座標、(x,y)は内部点の座標を示す。式(1)〜(3)に示すように、内部点の値z(x,y)は、x方向の補間を2回、y方向の補間を1回実施することで計算できる。
高精度補間処理部6が行う高精度補間処理の例として、3次スプライン補間について説明する。スプライン補間は、区間毎に多項式を作り、滑らかに離散データを繋ぐ補間方式である。3次スプライン補間の式は以下のとおりである。ただし、添え字のiはi番目の区間であることを示す。なお、この式は未知数の係数が4個あるので、離散データを最低4つ必要とする。
2次元補間の計算は、線形補間の場合と同様に、x方向の補間を2回、y方向の補間を1回実施することで計算できる。
画像再生処理部3は、高精度補間処理部6または低精度補間処理部5のうち、曲率判定部4が選択した補間処理の結果を用いて最終的なSAR画像を生成する(ステップST31)。その後、画像再生処理部3は、全ての対象領域に対して補間処理を実施したか否かを判断する(ステップST32)。全ての対象領域の補間処理が完了した場合には、ループ処理を終了し(ステップST33)、最終的なSAR画像をSAR画像格納部7へ出力する(ステップST34)。補間処理を実施していない対象領域がある場合には、ステップST25へ戻る。画像再生処理部3は、補間処理を行っていない対象領域を選択して処理を繰り返す。
ここで、曲率の判定方法について、図を用いて詳しく説明する。図7は、図5のi軸方向あるいはj軸方向のアジマス圧縮後データを1次元方向に並べたものである。例えば、j軸方向の位置を固定し、i軸方向の位置を変化させて得られる信号を取り出した場合には、横軸xはi軸方向、縦軸f(x)はアジマス圧縮後データの信号強度を示す。また、図7では、アジマス圧縮後データの横軸方向の位置は等間隔とする。
曲率判定部4は、対象領域の近傍で所定方向に隣接して並ぶアジマス圧縮後データを4点用いて曲率を算出する。ここでは、図に示すように、対象領域の近傍でi軸方向に隣接して並ぶデータA点、B点、C点、D点を使用し、B点とC点とを結ぶ辺の曲率を判定する場合について説明する。なお、以降の説明では、例えば、B点とC点との間の区間を区間BCと言う。また、区間ABと区間BCとをあわせて区間ABCと言う。
まず、曲率判定部4は、B点の位置における曲率を算出する。B点の位置における曲率は、A点とB点とを結ぶ線分を2:1に外分するX1点とC点との間の距離Δf1Rを計算することで求められる。同様に、曲率判定部4は、C点の位置における曲率を算出する。C点の位置における曲率は、B点とC点とを結ぶ線分を2:1に外分するX2点とD点との間の距離Δf2Rを計算することで求められる。これらの計算は、2階微分の計算と等価である。距離Δf1Rが第1の曲率、距離Δf2Rが第2の曲率にそれぞれ対応する。
次に、曲率判定部4は、算出した距離Δf1R,Δf2Rを基準値と比較する。距離Δf1R,Δf2Rは、曲率と比例関係にあることから、距離Δf1R,Δf2Rの少なくとも一方が基準値以上であれば、区間BCの曲率は大きいと判定する。また、曲率判定部4は、距離Δf1R,Δf2Rがともに基準値未満であれば、区間BCの曲率は小さいと判定する。図に示す例では、区間AB,BC,CD,DEは曲率が大きいため高精度補間処理部6が選択され、区間EF,FGは曲率が小さいため低精度補間処理部5が選択されることを示している。
なお、曲率判定部4は、対象領域を構成する全ての辺に対して上記の処理を行う。対象領域が図5に太線で示す格子領域の場合、4本の辺L1,L2,L3,L4に対して上記の処理を行う。その結果、少なくとも1つの辺において曲率が大きいと判定された場合には、その対象領域の曲率は大きいと判断する。全ての辺において曲率が小さいと判定された場合には、その対象領域の曲率は小さいと判断する。
また、本実施の形態においては、区間BCの曲率を判定するために、B点とC点の位置における曲率(第1の曲率と第2の曲率)を算出しているが、これに限定するものではない。例えば、いずれか一方で区間BCの曲率を判定してもよい。ただし、区間BCを含む2つの区間、具体的には、区間ABCと区間BCDとにおいて判定結果が異なる場合、両者は対等な関係にあるため、いずれか一方ではなく、両方の判定結果を考慮しなければならない。よって、B点とC点の位置において算出した曲率を用いることで、区間BCの曲率をより精度よく判定することができる。
次に、本実施の形態の曲率計算と従来の曲率計算とを比較する。信号強度をfの関数で表すと、A点の信号強度はf(A)で表される。このとき、本実施の形態の第1および第2の曲率の計算式は、下記のとおりとなる。
一方、従来手法では、図8に示す離散データ3点(A点,B点,C点)の場合、式(7)〜(9)に示すように、ベクトルABとベクトルBCとの規格化計算および内積計算により曲率が求められる。
このように、従来手法では、ベクトル規格化に伴う除算演算や平方根演算を用いるため、演算量が増大する。これに対して、本発明では、加算や減算のみで曲率が算出できるので、演算量を大幅に削減することができる。
次に、曲率の判定に用いる基準値の決定方法について、図を参照して説明する。図9は、本実施の形態の補間処理に伴う補間誤差の説明図である。破線は、真値と一致する理想的な補間点、一点鎖線は、実際の補間処理によって得られる補間点を表す。両者の間に生じるずれ量は、補間処理に伴う誤差(補間誤差)に相当する。この補間誤差と画像分解能とを関係付けることにより、基準値の決定が可能になる。
図10は、一点のターゲットに電波を照射して観測したときに得られる観測信号を用いて画像再生処理したときに得られるSAR画像であり、横軸はレンジ方向、縦軸はアジマス方向を示す。図10(a)は、全て真値の補間点が得られた場合、つまり、補間誤差が生じない場合のSAR画像で、信号の強度分布がアジマス方向とレンジ方向とで直交し、信号強度がピークとなる中央を基準として、上下左右方向に対称なことから、高分解能なSAR画像が再生されることを示している。一方、図10(b)は、補間処理の精度が悪く、補間誤差が大きい場合のSAR画像で、信号の強度分布がアジマス方向とレンジ方向とで直交しておらず、上下左右方向に非対称なことから、焦点がぼけた低分解能なSAR画像が再生されることを示している。
図10(a)において、ターゲット中心を基準として、レンジ方向、あるいはアジマス方向の信号強度分布に注目すると、図11に示すようなsinc関数と一致する。このsinc関数は、受信信号の電波の式から、理論的に導かれるものである。一方、図10(b)は、sinc関数の形状から崩れた信号強度分布となるため、ぼけたSAR画像となる。
上記では、SAR画像の定性的な見方について示したが、SAR画像の定量的な評価の一つとして、画像分解能が用いられる。画像分解能は、画像のきめ細かさを表す指標であり、値が小さいほど対象物間の区別がしやすく鮮明な画像となる。
画像分解能は、一般に信号強度がピークから半分になるまでの信号区間の長さを意味し、レンジ方向、アジマス方向でそれぞれ定義される。信号強度をデシベルで表すと、図11に示すように、ピーク値で0dB、半分の強度で−3dBとなるので、画像分解能を3dB幅とも言う。スラントレンジ方向の画像分解能の理論値δsrとアジマス方向の画像分解能の理論値δazは、以下の式(10)、(11)で決定される。
ここで、Cは光速、Bは電波の送信帯域幅、λは電波の波長、θSAは合成開口角である。合成開口角とは、固定位置に対して、電波ビームの当たり始めから当たり終わりまでのプラットフォームと固定位置との角度のことを意味する。図12にSAR観測の観測幾何を示す。SAR観測では、図のように、プラットフォームの進行方向に対して斜め下方向に電波ビームを照射する。アジマス方向はプラットフォーム進行方向であり、スラントレンジ方向はプラットフォームから電波ビームを照射する方向である。図では、電波ビームの角度を固定するSAR観測を表しており、この場合、電波ビームの角度が合成開口角と等しくなる。
図10は、1点のターゲットに電波を照射した際にできるSAR画像であるが、無数のターゲットの反射を観測して画像再生処理を行うと、一般に航空写真で良くみられるような地形や植生分布を表すSAR画像が得られる。
SAR画像の分解能が、上記の分解能の理論値に一致するとき、理想的な再生画像となることが分かっているので、SAR画像の分解能と補間誤差とを対応づけることにより、曲率判定用の閾値である基準値を決定することができる。つまり、事前に分解能と補間誤差との関係を調べておき、必要とする分解能に相当する補間誤差を基準値として設定すればよい。例えば、分解能が理論値の数%以内に収まる範囲に相当する補間誤差を基準値として定義することができる。
レンジ方向とアジマス方向では、SAR画像の分解能と補間誤差との対応関係が異なるので、レンジ方向とアジマス方向のそれぞれで対応関係を求めることができる。ここでは、SAR画像の分解能と補間誤差との対応関係について、レンジ方向の例を挙げて詳しく説明する。まず、SAR画像の再生処理前の模擬信号データを連続的な関数で表す。関数の式は、例えば、下記の式(12)で表される受信信号波形を使用する。この信号波形では、振幅が1、exp関数の中身が位相に相当する。
ここで、τはレンジ時間、SR(τ)はレンジ時間τにおける信号データの輝度値、f0は送信チャープパルスの中心周波数、τdは点目標までの電波の往復時間、Krは送信チャープパルスのチャープレートを表す。
式(12)で表される関数上でサンプリングした点を理想的な補間処理により得られる結果と仮定する。次に、理想的な補間処理結果に対し、振幅値を減少させる操作をする。振幅値の減少量は、時間の経過とともに変化させる。図13は、模擬信号データの信号波形を示す図である。図13(a)は理想的な補間処理結果と仮定した信号波形、図13(b)および図13(c)は振幅値を減少させた補間処理結果の信号波形を表す。図13(b)では、平均すると振幅値をd1減少させている。図13(c)では、平均すると振幅値をd2(d2>d1)減少させている。Trは送信チャープパルスの幅を表す。
振幅値の減少量が大きいほど、SAR画像の再生処理後に計測される分解能は低下する。つまり、分解能の理論値から大きくずれる。よって、振幅値の操作量の大きさと分解能の理論値からの誤差との関係を求めることで、SAR画像の画素値の真値からのずれ量と、SAR画像の分解能の理論値からの誤差との関係を求めることができる。上記では、レンジ方向に関する対応関係を説明したが、アジマス方向に関しても同様に、アジマス方向の受信信号波形と合成開口時間とを用いて、アジマス分解能と補間誤差との対応関係を求めればよい。具体的には、送信チャープパルス幅Tr分のレンジ方向の受信信号波形の代わりに、合成開口時間分のアジマス方向の受信信号波形を用いる点が異なる。
図14は、図13に示す3つの信号波形を画像再生処理した結果から得られる分解能をプロットした示した図であり、SAR画像の分解能と補間誤差との関係を示す図である。横軸は、理想的な補間処理結果から振幅値を減少させた量、つまり真値からのずれ量に対応する。縦軸は、分解能の理論値からの誤差を示す。この関係から、画像再生において許容できる分解能の理論値からの誤差を任意に決定し、その誤差に対応する真値からのずれ量を曲率判定の基準値として設定する。
以上で明らかなように、実施の形態1の合成開口レーダ信号処理装置100によれば、合成開口レーダの受信信号から得られる離散データを補間処理する低精度補間処理部と、低精度補間処理部より高い精度で離散データを補間処理する高精度補間処理部と、補間処理の対象領域における離散データの曲率である第1の曲率に応じて高精度補間処理部または低精度補間処理部のいずれか一方を選択する曲率判定部と、曲率判定部で選択した補間処理の結果を用いて画像を再生する画像再生処理部とを備え、曲率判定部は、対象領域の近傍で所定方向に隣接して並ぶ離散データA点、B点、C点を選択し、A点とB点とを結ぶ線分を2:1に外分する点とC点との間の距離Δf1Rに基づいて第1の曲率を決定する。このように、離散データであるアジマス圧縮後データの信号強度の曲率を加算や減算のみで算出することができるので、曲率算出の演算量を抑制することができる。
また、実施の形態1の合成開口レーダ信号処理装置100によれば、処理の対象領域は、離散データを結ぶ線分を辺とする領域であり、曲率判定部は、辺ごとに辺の両端の離散データを点B、点Cとして第1の曲率を決定し、辺ごとに決定した第1の曲率に応じて高精度補間処理部または低精度補間処理部のいずれか一方を選択することを特徴とする。これにより、補間処理の判定にかかる演算量を抑えることができる。
さらに、曲率判定部は、第1の曲率の絶対値が閾値を超える辺が存在する場合には高精度補間処理(例えば2次元スプライン補間)を選択し、そうでない場合には低精度補間処理(例えば2次元線形補間)を選択する。これにより、曲率の大きい領域で高精度な補間処理を行い、曲率の小さい領域で低精度な補間処理を行うので、演算量を抑えるとともに、十分な補間精度をもつ最終的なSAR画像を得ることができる。
また、実施の形態1の合成開口レーダ信号処理装置100によれば、曲率判定部は、さらに対象領域の近傍で所定方向に隣接して並ぶ離散データB点、C点、D点を選択し、B点とC点とを結ぶ線分を2:1に外分する点とD点との間の距離Δf2Rに基づいて第2の曲率を決定し、第1の曲率および第2の曲率に応じて高精度補間処理部または低精度補間処理部のいずれか一方を選択する。これにより、補間処理の判定を高精度に行うことができる。
なお、高精度補間処理および低精度補間処理は、上述した以外の補間式を使用してもよい。例えば、高精度補間処理であれば多項式補間、低精度補間処理であれば最近傍法でもよい。演算量に大小関係があれば、補間処理の種類は問わない。
また、曲率判定部4では、高精度補間処理あるいは低精度補間処理の2種類の判定について説明したが、判定基準の数を増やすことにより、3つ以上の補間処理の中から選択する構成としてもよい。
なお、本実施の形態では、2次元補間処理の判定手法を例として挙げたが、1次元や3次元の補間処理の判定を上記で示した判定手法により実施してもよい。
実施の形態2.
本実施の形態における合成開口レーダ装置信号処理装置について説明する。実施の形態1では、アジマス圧縮後データが縦方向あるいは横方向に等間隔に配置される場合について説明した。本実施の形態では、アジマス圧縮後データが非等間隔、つまり、任意の間隔で配置される点が実施の形態1と異なる。実施の形態2では、この異なる点を中心に説明する。
本実施の形態における合成開口レーダ装置信号処理装置について説明する。実施の形態1では、アジマス圧縮後データが縦方向あるいは横方向に等間隔に配置される場合について説明した。本実施の形態では、アジマス圧縮後データが非等間隔、つまり、任意の間隔で配置される点が実施の形態1と異なる。実施の形態2では、この異なる点を中心に説明する。
図15は、非等間隔に配置された離散データの曲率の算出手法を示す図である。横軸xは離散データであるアジマス圧縮後データの配置位置、縦軸f(x)はアジマス圧縮後データの信号強度を示す。ここでは、x方向に非等間隔に並ぶ離散データA’点、B’点、C’点を用いて、B’点の位置における曲率の算出方法を説明する。なお、A’点とB’点との間の距離をA’B’、B’点とC’点との間の距離をB’C’、A’点とB’点との間のx方向の距離をΔXL、B’点とC’点との間のx方向の距離をΔXRとする。
まず、実施の形態1と同様に、曲率判定部4は、3点A’点、B’点、C’点を用いて、A’点とB’点とを結ぶ線分を(A’B’+B’C’):B’C’に外分するX1R’点とC’点との間の距離Δf1Rを計算する。また、B’点とC’点とを結ぶ線分をA’B’:(A’B’+B’C’)に外分するX1L’点とA’点との間の距離Δf1Lを計算する。さらに、曲率判定部4は、距離Δf1Rと距離Δf1Lとの平均を算出し、その結果を用いてB’点の位置における曲率である第1の曲率を判定する。曲率判定部4は、C’点の位置についても同様に、距離Δf2Rと距離Δf2Lを計算し、その平均を用いてC’点の位置における曲率である第2の曲率を判定する。
距離Δf1Rと距離Δf1Lとの平均をとる理由について説明する。離散データが等間隔の場合には、距離Δf1Rと距離Δf1Lとは等しい。そのため、実施の形態1では距離Δf1Rのみ算出して曲率を判定した。一方、離散データが非等間隔の場合には、距離Δf1Rと距離Δf1LとがA’B’とB’C’との比で変わる。例えば、図15のようにB’点がA’点よりもC’点に近い場合、距離Δf1Rが小さくなり、曲率判定で過小評価される。この問題を回避するため、本実施の形態では、距離Δf1Rと距離Δf1Lとの平均をとり、その結果を用いて曲率を判定する。第1の曲率の計算式は、次式のようになる。
なお、実施の形態1のように離散データが等間隔の場合には、距離Δf1Rと距離Δf1Lは等しいため、距離Δf1Rの代わりに距離Δf1Lを使用して曲率判定してもよい。同様に、距離Δf2Rと距離Δf2Lも等しいため、距離Δf2Rの代わりに距離Δf2Lを使用して曲率判定してもよい。
以上で明らかなように、実施の形態2の合成開口レーダ信号処理装置によれば、合成開口レーダの受信信号から得られる離散データを補間処理する低精度補間処理部と、低精度補間処理部より高い精度で離散データを補間処理する高精度補間処理部と、補間処理の対象領域における離散データの曲率である第1の曲率に応じて高精度補間処理部または低精度補間処理部のいずれか一方を選択する曲率判定部と、曲率判定部で選択した補間処理の結果を用いて画像を再生する画像再生処理部とを備え、曲率判定部は、対象領域の近傍で所定方向に隣接して並ぶ離散データA点、B点、C点を選択し、A点とB点とを結ぶ線分を(AB+BC):BCに外分する点とC点との間の距離Δf1Rに基づいて第1の曲率を決定する。このように、離散データであるアジマス圧縮後データが非等間隔の場合でも、従来と比較して少ない演算量で信号強度の曲率を求めることができる。
また、実施の形態1の合成開口レーダ信号処理装置100によれば、曲率判定部は、B点とC点とを結ぶ線分をAB:(AB+BC)に外分する点とA点との間の距離Δf1Lおよび距離Δf1Rに基づいて第1の曲率を決定する。これにより、離散データであるアジマス圧縮後データが非等間隔の場合でも、補間処理の判定を高精度に行うことができる。
本発明はその発明の範囲内において、各実施の形態の自由な組み合わせ、あるいは各実施の形態の任意の構成要素の変形、もしくは各実施の形態において任意の構成要素の省略が可能である。
この発明に係る合成開口レーダ信号処理装置は、地表や海面等のSAR画像を再生する際、画像再生処理での十分な演算精度を確保しながら、演算量を削減する必要性が高いものに適している。
1 SARセンサ(合成開口レーダ)、2 受信データ格納部、3 画像再生処理部、4 曲率判定部、5 低精度補間処理部、6 高精度補間処理部、7 SAR画像格納部、11 メモリ、12 プロセッサ、13 記録装置、100 合成開口レーダ信号処理装置。
Claims (13)
- 合成開口レーダの受信信号から得られる離散データを補間処理する低精度補間処理部と、
前記低精度補間処理部より高い精度で前記離散データを補間処理する高精度補間処理部と、
補間処理の対象領域における前記離散データの曲率である第1の曲率に応じて前記高精度補間処理部または前記低精度補間処理部のいずれか一方を選択する曲率判定部と、
前記曲率判定部で選択した補間処理の結果を用いて画像を再生する画像再生処理部と
を備え、
前記曲率判定部は、前記対象領域の近傍で所定方向に隣接して並ぶ離散データA点、B点、C点を選択し、前記A点と前記B点とを結ぶ線分を(AB+BC):BCに外分する点と前記C点との間の距離Δf1Rに基づいて前記第1の曲率を決定する合成開口レーダ信号処理装置。
(ただし、ABはA点とB点との間の距離、BCはB点とC点との間の距離。) - 前記曲率判定部は、前記B点と前記C点とを結ぶ線分をAB:(AB+BC)に外分する点と前記A点との間の距離Δf1Lおよび前記距離Δf1Rに基づいて前記第1の曲率を決定することを特徴とする請求項1に記載の合成開口レーダ信号処理装置。
- 前記曲率判定部は、距離Δf1Lと前記距離Δf1Rとの平均値を前記第1の曲率とすることを特徴とする請求項2に記載の合成開口レーダ信号処理装置。
- 前記曲率判定部は、前記第1の曲率の絶対値が閾値を超える場合には前記高精度補間処理部を選択し、そうでない場合には前記低精度補間処理部を選択することを特徴とする請求項1から請求項3のいずれか1項に記載の合成開口レーダ信号処理装置。
- 前記曲率判定部は、さらに前記対象領域の近傍で前記所定方向に隣接して並ぶ離散データB点、C点、D点を選択し、前記B点と前記C点とを結ぶ線分を(BC+CD):CDに外分する点と前記D点との間の距離Δf2Rに基づいて第2の曲率を決定し、前記第1の曲率および前記第2の曲率に応じて前記高精度補間処理部または前記低精度補間処理部のいずれか一方を選択することを特徴とする請求項1に記載の合成開口レーダ信号処理装置。
(ただし、CDはC点とD点との間の距離。) - 前記曲率判定部は、前記C点と前記D点とを結ぶ線分をBC:(BC+CD)に外分する点と前記D点との間の距離Δf2Lおよび前記距離Δf2Rに基づいて前記第2の曲率を決定することを特徴とする請求項5に記載の合成開口レーダ信号処理装置。
- 前記曲率判定部は、距離Δf2Lと前記距離Δf2Rとの平均値を前記第2の曲率とすることを特徴とする請求項6に記載の合成開口レーダ信号処理装置。
- 前記曲率判定部は、前記第1の曲率の絶対値および前記第2の曲率の絶対値の少なくとも一方が閾値を超える場合には前記高精度補間処理部を選択し、そうでない場合には前記低精度補間処理部を選択することを特徴とする請求項5から請求項7のいずれか1項に記載の合成開口レーダ信号処理装置。
- 前記対象領域は、前記離散データを結ぶ線分を辺とする領域であり、
前記曲率判定部は、前記辺ごとに前記辺の両端の前記離散データを前記点B、前記点Cとして前記第1の曲率を決定し、前記辺ごとに決定した前記第1の曲率に応じて前記高精度補間処理部または前記低精度補間処理部のいずれか一方を選択することを特徴とする請求項1から請求項4のいずれか1項に記載の合成開口レーダ信号処理装置。 - 前記曲率判定部は、前記第1の曲率の絶対値が閾値を超える辺が存在する場合には前記高精度補間処理部を選択し、そうでない場合には前記低精度補間処理部を選択することを特徴とする請求項9に記載の合成開口レーダ信号処理装置。
- 前記曲率判定部は、前記離散データを結ぶ線分を辺とする領域により前記対象領域を規定し、前記辺ごとに前記辺の両端の前記離散データを前記点B、前記点Cとして前記第1の曲率および前記第2の曲率を決定し、前記辺ごとに決定した前記第1の曲率および前記第2の曲率に応じて前記高精度補間処理部または前記低精度補間処理部のいずれか一方を選択することを特徴とする請求項5から請求項8のいずれか1項に記載の合成開口レーダ信号処理装置。
- 前記曲率判定部は、前記第1の曲率の絶対値および前記第2の曲率の絶対値の少なくとも一方が閾値を超える辺が存在する場合には前記高精度補間処理部を選択し、そうでない場合には前記低精度補間処理部を選択することを特徴とする請求項11に記載の合成開口レーダ信号処理装置。
- 前記閾値は、前記画像の画素値の真値からのずれ量と、前記画像の分解能の理論値からの誤差との関係を予め求め、前記誤差が所定の値となるときの前記ずれ量を前記閾値とすることを特徴とする請求項4、請求項8、請求項10および請求項12のいずれか1項に記載の合成開口レーダ信号処理装置。
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
PCT/JP2016/057290 WO2017154125A1 (ja) | 2016-03-09 | 2016-03-09 | 合成開口レーダ信号処理装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
JP6261839B1 true JP6261839B1 (ja) | 2018-01-17 |
JPWO2017154125A1 JPWO2017154125A1 (ja) | 2018-03-15 |
Family
ID=59789195
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2017549349A Expired - Fee Related JP6261839B1 (ja) | 2016-03-09 | 2016-03-09 | 合成開口レーダ信号処理装置 |
Country Status (5)
Country | Link |
---|---|
US (1) | US10921441B2 (ja) |
EP (1) | EP3418770B1 (ja) |
JP (1) | JP6261839B1 (ja) |
CA (1) | CA3015743C (ja) |
WO (1) | WO2017154125A1 (ja) |
Families Citing this family (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111338364B (zh) * | 2019-11-21 | 2021-09-21 | 浙江大学 | 快速响应的高超声速飞行器轨迹优化高精度控制器 |
EP3913398B1 (en) * | 2020-05-20 | 2022-11-16 | Airbus Defence and Space GmbH | Synthetic aperture radar data reduction for satellites |
JP7188653B1 (ja) * | 2021-02-25 | 2022-12-13 | 三菱電機株式会社 | データ処理装置およびレーダ装置 |
CN114244907B (zh) * | 2021-11-23 | 2024-01-16 | 华为技术有限公司 | 雷达数据的压缩方法和装置 |
CN115564728B (zh) * | 2022-09-30 | 2023-08-11 | 苏州大学 | 一种图像角点检测方法、装置、设备及应用 |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5243349A (en) * | 1981-03-17 | 1993-09-07 | Westinghouse Electric Corp. | High resolution synthetic aperture radar having rectilinear output image format |
JPH07195253A (ja) * | 1993-11-24 | 1995-08-01 | Hitachi Ltd | 加工用camシステムの処理方法 |
JP2001281529A (ja) * | 2000-03-29 | 2001-10-10 | Minolta Co Ltd | デジタルカメラ |
US7532150B1 (en) * | 2008-03-20 | 2009-05-12 | Raytheon Company | Restoration of signal to noise and spatial aperture in squint angles range migration algorithm for SAR |
JP2009128019A (ja) * | 2007-11-20 | 2009-06-11 | Mitsubishi Electric Corp | 合成開口レーダ画像再生装置、合成開口レーダ画像再生方法及び合成開口レーダ画像再生プログラム |
WO2011018943A1 (ja) * | 2009-08-12 | 2011-02-17 | 日本電気株式会社 | データ要約システム、データ要約方法および記録媒体 |
JP2012063152A (ja) * | 2010-09-14 | 2012-03-29 | Mitsubishi Space Software Kk | 速度推定装置及びコンピュータプログラム及び速度推定方法 |
JP2014153161A (ja) * | 2013-02-07 | 2014-08-25 | Mitsubishi Electric Corp | 画像レーダ装置 |
Family Cites Families (17)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5179383A (en) * | 1991-07-15 | 1993-01-12 | Raney R K | Synthetic aperture radar processor to handle large squint with high phase and geometric accuracy |
US5365465A (en) * | 1991-12-26 | 1994-11-15 | Texas Instruments Incorporated | Floating point to logarithm converter |
US5667373A (en) * | 1994-08-05 | 1997-09-16 | Acuson Corporation | Method and apparatus for coherent image formation |
JPH09288166A (ja) * | 1996-04-24 | 1997-11-04 | Mitsubishi Electric Corp | レーダ走査変換装置 |
US7336820B2 (en) * | 2005-02-04 | 2008-02-26 | Kabushiki Kaisha Toshiba | Method and apparatus for rapid shading in a raster image processor |
US9336628B2 (en) * | 2008-05-08 | 2016-05-10 | United Technologies Corporation | Systems and methods involving surface fitting |
US9075129B2 (en) * | 2008-12-10 | 2015-07-07 | The United States Of America As Represented By The Secretary Of The Army | Method and system for forming images by comparing subsets of image data |
JP2010286359A (ja) * | 2009-06-11 | 2010-12-24 | Furuno Electric Co Ltd | 信号処理装置、この信号処理装置を備えるレーダ装置及びその信号処理装置の画素抜け検出方法 |
JP5731800B2 (ja) * | 2010-11-19 | 2015-06-10 | 古野電気株式会社 | 物標探知方法、物標探知プログラム、物標探知装置、およびレーダ装置 |
JP5699621B2 (ja) * | 2011-01-19 | 2015-04-15 | 株式会社リコー | 画像処理装置、画素補間方法およびプログラム |
US9976887B1 (en) * | 2011-06-22 | 2018-05-22 | Daniel T. Mudd | Wider dynamic accuracy range for gas delivery devices |
MY170961A (en) * | 2012-03-07 | 2019-09-20 | Nippon Steel Corp | Misalignment calculation system |
JP6047968B2 (ja) * | 2012-07-17 | 2016-12-21 | セイコーエプソン株式会社 | プロジェクター、及び、プロジェクターにおける発光制御方法 |
US9179039B2 (en) * | 2014-02-12 | 2015-11-03 | Xerox Corporation | Methods and systems for processing low resolution images via error diffusion |
JP6257434B2 (ja) * | 2014-04-24 | 2018-01-10 | 三菱電機株式会社 | 画像レーダ装置 |
US10624612B2 (en) * | 2014-06-05 | 2020-04-21 | Chikayoshi Sumi | Beamforming method, measurement and imaging instruments, and communication instruments |
JP6493196B2 (ja) * | 2015-12-17 | 2019-04-03 | 株式会社デンソー | 制御装置、制御方法 |
-
2016
- 2016-03-09 US US16/071,760 patent/US10921441B2/en active Active
- 2016-03-09 EP EP16893457.8A patent/EP3418770B1/en active Active
- 2016-03-09 JP JP2017549349A patent/JP6261839B1/ja not_active Expired - Fee Related
- 2016-03-09 WO PCT/JP2016/057290 patent/WO2017154125A1/ja active Application Filing
- 2016-03-09 CA CA3015743A patent/CA3015743C/en active Active
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5243349A (en) * | 1981-03-17 | 1993-09-07 | Westinghouse Electric Corp. | High resolution synthetic aperture radar having rectilinear output image format |
JPH07195253A (ja) * | 1993-11-24 | 1995-08-01 | Hitachi Ltd | 加工用camシステムの処理方法 |
JP2001281529A (ja) * | 2000-03-29 | 2001-10-10 | Minolta Co Ltd | デジタルカメラ |
JP2009128019A (ja) * | 2007-11-20 | 2009-06-11 | Mitsubishi Electric Corp | 合成開口レーダ画像再生装置、合成開口レーダ画像再生方法及び合成開口レーダ画像再生プログラム |
US7532150B1 (en) * | 2008-03-20 | 2009-05-12 | Raytheon Company | Restoration of signal to noise and spatial aperture in squint angles range migration algorithm for SAR |
WO2011018943A1 (ja) * | 2009-08-12 | 2011-02-17 | 日本電気株式会社 | データ要約システム、データ要約方法および記録媒体 |
JP2012063152A (ja) * | 2010-09-14 | 2012-03-29 | Mitsubishi Space Software Kk | 速度推定装置及びコンピュータプログラム及び速度推定方法 |
JP2014153161A (ja) * | 2013-02-07 | 2014-08-25 | Mitsubishi Electric Corp | 画像レーダ装置 |
Also Published As
Publication number | Publication date |
---|---|
CA3015743A1 (en) | 2017-09-14 |
EP3418770A1 (en) | 2018-12-26 |
CA3015743C (en) | 2020-05-26 |
US10921441B2 (en) | 2021-02-16 |
US20200166635A1 (en) | 2020-05-28 |
WO2017154125A1 (ja) | 2017-09-14 |
JPWO2017154125A1 (ja) | 2018-03-15 |
EP3418770B1 (en) | 2020-08-05 |
EP3418770A4 (en) | 2019-02-27 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP6261839B1 (ja) | 合成開口レーダ信号処理装置 | |
US7397418B1 (en) | SAR image formation with azimuth interpolation after azimuth transform | |
CN105842694B (zh) | 一种基于ffbp sar成像的自聚焦方法 | |
Gorham et al. | Scene size limits for polar format algorithm | |
JP4749903B2 (ja) | レーダ画像処理装置 | |
JP6249110B1 (ja) | 合成開口レーダ装置 | |
CN104914440B (zh) | 一种融合动态孔径的压缩感知超声波束合成方法 | |
CN115166714B (zh) | 单通道sar运动舰船二维速度估计与重定位方法及装置 | |
Cheng et al. | Processing of airborne video SAR data using the modified back projection algorithm | |
JP2004198275A (ja) | 合成開口レーダ装置および画像再生方法 | |
CN113608218B (zh) | 一种基于后向投影原理的频域干涉相位稀疏重构方法 | |
JP6289672B2 (ja) | 合成開口レーダ信号処理装置及び合成開口レーダ信号処理プログラム | |
CN112415512B (zh) | 基于进退法和黄金分割法的sar运动目标聚焦方法 | |
JP5664869B2 (ja) | 測定装置、測定システム、測定方法、及びプログラム | |
Xingyu et al. | Approach for ISAR imaging of near-field targets based on coordinate conversion and image interpolation | |
CN116794643A (zh) | 一种基于声反射断层扫描成像的逆投影方法 | |
Kim et al. | Synthetic aperture imaging of dispersive targets | |
CN110632606B (zh) | 一种基于空间重采样的二维快速傅里叶变换三维成像方法 | |
CN112666625A (zh) | 用于毫米波安检的快速成像设备及其方法 | |
Wang et al. | An accelerate FFBP algorithm of bistatic forward-looking SAR based on azimuth equidistant coordinates | |
CN109581377B (zh) | 一种点回波成像质量的分辨率评估方法、设备及存储介质 | |
JP6289388B2 (ja) | 画像レーダ装置 | |
WO2024150806A1 (ja) | 信号処理装置、信号処理方法、及び記録媒体 | |
CN112327299B (zh) | 电磁涡旋前视目标三维成像方法 | |
WO2024150807A1 (ja) | 信号処理装置、信号処理方法、及び記録媒体 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A975 | Report on accelerated examination |
Free format text: JAPANESE INTERMEDIATE CODE: A971005 Effective date: 20171013 |
|
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: 20171114 |
|
A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20171212 |
|
R150 | Certificate of patent or registration of utility model |
Ref document number: 6261839 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R150 |
|
LAPS | Cancellation because of no payment of annual fees |