JPH09281231A - 信号処理装置 - Google Patents
信号処理装置Info
- Publication number
- JPH09281231A JPH09281231A JP8092710A JP9271096A JPH09281231A JP H09281231 A JPH09281231 A JP H09281231A JP 8092710 A JP8092710 A JP 8092710A JP 9271096 A JP9271096 A JP 9271096A JP H09281231 A JPH09281231 A JP H09281231A
- Authority
- JP
- Japan
- Prior art keywords
- unit
- image
- data
- correlation value
- registration
- 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.)
- Pending
Links
Landscapes
- Radar Systems Or Details Thereof (AREA)
Abstract
(57)【要約】
【課題】 画像データの中で標高の変化のほとんどない
領域について後段の位置合わせを行わないようにするこ
とにより、少ない演算量で、より精度の高い位置合わせ
を行う信号処理装置が必要である。 【解決手段】 通常の信号処理装置にレジストレーショ
ンを行いやすいサイズに画像を分割する画像分割部９、
分割画像についてそれぞれ振幅の分散を演算する分散演
算部１０、分散の小さな領域についてはレジストレーシ
ョンを行わないと判定する処理領域判定部１１、処理す
べき領域より画像を切り出す第２の画像切り出し部１２
を付加することにより構成される。
領域について後段の位置合わせを行わないようにするこ
とにより、少ない演算量で、より精度の高い位置合わせ
を行う信号処理装置が必要である。 【解決手段】 通常の信号処理装置にレジストレーショ
ンを行いやすいサイズに画像を分割する画像分割部９、
分割画像についてそれぞれ振幅の分散を演算する分散演
算部１０、分散の小さな領域についてはレジストレーシ
ョンを行わないと判定する処理領域判定部１１、処理す
べき領域より画像を切り出す第２の画像切り出し部１２
を付加することにより構成される。
Description
【０００１】
【発明の属する技術分野】この発明は、複数の合成開口
レーダ画像を干渉させて標高データを取得する信号処理
装置に関する。
レーダ画像を干渉させて標高データを取得する信号処理
装置に関する。
【０００２】
【従来の技術】図６は、例えばＥＮｈｏｌｌａｎｄ著Ｒ
ｅｐｏｒｔ ｏｎ Ｉｎｔｅｒｆｅｒｏｍｅｔｒｉｃ
ＳＡＲ ａｎｄ ｉｔｓ ａｐｐｌｉｃａｔｉｏｎｓ
（ＴＲ−９３／２３５／１．０）に示された従来のイン
タフェロメトリ処理装置の構成ブロック図である。図６
において１は画像データ、２は第２の画像データ、３は
第１の画像切り出し部、４は第１のレジストレーション
部、５は干渉処理部、６は平面位相除去部、７はマルチ
ルック処理部、８はアンラッピング処理部、９はデータ
変換部である。
ｅｐｏｒｔ ｏｎ Ｉｎｔｅｒｆｅｒｏｍｅｔｒｉｃ
ＳＡＲ ａｎｄ ｉｔｓ ａｐｐｌｉｃａｔｉｏｎｓ
（ＴＲ−９３／２３５／１．０）に示された従来のイン
タフェロメトリ処理装置の構成ブロック図である。図６
において１は画像データ、２は第２の画像データ、３は
第１の画像切り出し部、４は第１のレジストレーション
部、５は干渉処理部、６は平面位相除去部、７はマルチ
ルック処理部、８はアンラッピング処理部、９はデータ
変換部である。
【０００３】次に動作について説明する。インタフェロ
メトリでは２個のアンテナで受信されるデータの位相差
を検知することにより、標高データを取得する。標高デ
ータＺは検出された位相φを利用して数１で表わすこと
ができる。また、このときのジオメトリを図７に示す。
メトリでは２個のアンテナで受信されるデータの位相差
を検知することにより、標高データを取得する。標高デ
ータＺは検出された位相φを利用して数１で表わすこと
ができる。また、このときのジオメトリを図７に示す。
【０００４】
【数１】
【０００５】数１においてｄはセンサ１，２それぞれと
観測点の距離の差、λは波長、αはセンサ２から見たセ
ンサ１と水平軸のなす角、θはセンサ１から見た観測点
の仰角、ｒはセンサ１と観測点の距離、Ｂはセンサ１と
センサ２の距離、ｈはセンサ１の高度である。数１より
わかるようにインタフェロメトリにて標高データを求め
る際には様々なパラメータが影響している。特にセンサ
の位置は一意的に定まらず、精度を向上させるためには
２つの画像間の相関を挙げることが重要であり、いくつ
かの処理が行われている。従来のインタフェロメトリ処
理では、時間的にまたは空間的にわずかにずれた画像デ
ータ１，２について第１の画像切り出し部３にて後段の
処理が可能なサイズに切り出す。そして第１のレジスト
レーション部４では切り出した２つの画像データについ
て位置合わせを行い、干渉処理部５にて位相差を求め、
干渉縞を得る。ここで位相差は数２によって求められ
る。
観測点の距離の差、λは波長、αはセンサ２から見たセ
ンサ１と水平軸のなす角、θはセンサ１から見た観測点
の仰角、ｒはセンサ１と観測点の距離、Ｂはセンサ１と
センサ２の距離、ｈはセンサ１の高度である。数１より
わかるようにインタフェロメトリにて標高データを求め
る際には様々なパラメータが影響している。特にセンサ
の位置は一意的に定まらず、精度を向上させるためには
２つの画像間の相関を挙げることが重要であり、いくつ
かの処理が行われている。従来のインタフェロメトリ処
理では、時間的にまたは空間的にわずかにずれた画像デ
ータ１，２について第１の画像切り出し部３にて後段の
処理が可能なサイズに切り出す。そして第１のレジスト
レーション部４では切り出した２つの画像データについ
て位置合わせを行い、干渉処理部５にて位相差を求め、
干渉縞を得る。ここで位相差は数２によって求められ
る。
【０００６】
【数２】
【０００７】数２においてｆ３は位相差データ、ｆ１は
画像データ１の複素データ、ｆ２*が画像データ２の複
素共役である。この干渉縞は観測地域における地球平面
の変動による位相差を含んでいるので平面位相除去部６
にてこの平面位相を除去する。そして位相ノイズを低減
するためにマルチルック処理７にて積分処理を行い、ア
ンラッピング処理部８にて位相積算を行う。これは数２
で求められた位相差は０から２πの範囲となっており、
標高データを求めるにはこの位相差を連続的に求める必
要がある。例えばレンジ方向に０，０．５π，π，１．
７π，０．１π・・・と変化していれば、１．７πの次
は２．１πとなる。最後にデータ変換部９にて位相差デ
ータを標高データに変換する。
画像データ１の複素データ、ｆ２*が画像データ２の複
素共役である。この干渉縞は観測地域における地球平面
の変動による位相差を含んでいるので平面位相除去部６
にてこの平面位相を除去する。そして位相ノイズを低減
するためにマルチルック処理７にて積分処理を行い、ア
ンラッピング処理部８にて位相積算を行う。これは数２
で求められた位相差は０から２πの範囲となっており、
標高データを求めるにはこの位相差を連続的に求める必
要がある。例えばレンジ方向に０，０．５π，π，１．
７π，０．１π・・・と変化していれば、１．７πの次
は２．１πとなる。最後にデータ変換部９にて位相差デ
ータを標高データに変換する。
【０００８】
【発明が解決しようとする課題】従来のこの種の信号処
理装置は、以上のように構成されているので、切り出し
画像データに海域のような標高の変化がほとんどない部
分を含む場合、２つの画像の位置合わせが困難であると
いう問題点があった。
理装置は、以上のように構成されているので、切り出し
画像データに海域のような標高の変化がほとんどない部
分を含む場合、２つの画像の位置合わせが困難であると
いう問題点があった。
【０００９】また、２つの画像の位置合わせを行うとき
に画像全体について位置合わせを調べていたため、大幅
な処理時間がかかるという問題点があった。
に画像全体について位置合わせを調べていたため、大幅
な処理時間がかかるという問題点があった。
【００１０】また、２つの画像の位置合わせを行うとき
にレンジ、アジマス、レンジ伸縮率についてそれぞれ位
置合わせを調べていたため、大幅な処理時間がかかると
いう問題点があった。
にレンジ、アジマス、レンジ伸縮率についてそれぞれ位
置合わせを調べていたため、大幅な処理時間がかかると
いう問題点があった。
【００１１】また、２つの画像の位置合わせを行うとき
に相関値により位置合わせを決定していたため、相関値
が悪い２つの画像の位置合わせが困難であるという問題
点があった。
に相関値により位置合わせを決定していたため、相関値
が悪い２つの画像の位置合わせが困難であるという問題
点があった。
【００１２】さらに、目標信号対ノイズ比（以下Ｓ／Ｎ
比と示す）が悪い画像では位相雑音が多いため、位相差
データがぼやけるという問題点があった。
比と示す）が悪い画像では位相雑音が多いため、位相差
データがぼやけるという問題点があった。
【００１３】この発明は上記のような課題を解消するた
めになされたもので、インタフェロメトリ処理におい
て、精度の高い標高データの取得、そして演算量の減
少、さらには一連の処理を自動的に行う信号処理装置を
得ることを目的とする。
めになされたもので、インタフェロメトリ処理におい
て、精度の高い標高データの取得、そして演算量の減
少、さらには一連の処理を自動的に行う信号処理装置を
得ることを目的とする。
【００１４】
【課題を解決するための手段】第１の発明による信号処
理装置は、通常の第１のレジストレーション部、干渉処
理部、平面位相除去部、マルチルック処理部、アンラッ
ピング部、データ変換部に画像分割する画像分割部と、
分割した画像データについてそれぞれ分散を求める分散
演算部と、分散値により後段の処理を行うか否か判定す
る処理領域判定部と、処理すべき領域の中から画像を切
り出す第２の画像切り出し部とを設けたものである。
理装置は、通常の第１のレジストレーション部、干渉処
理部、平面位相除去部、マルチルック処理部、アンラッ
ピング部、データ変換部に画像分割する画像分割部と、
分割した画像データについてそれぞれ分散を求める分散
演算部と、分散値により後段の処理を行うか否か判定す
る処理領域判定部と、処理すべき領域の中から画像を切
り出す第２の画像切り出し部とを設けたものである。
【００１５】また、第２の発明による信号処理装置は、
通常の第１の画像切り出し部、干渉処理部、平面位相除
去部、マルチルック処理部、アンラッピング処理部、デ
ータ変換部に、切り出した画像データの中の数カ所につ
いて２つの画像の位置合わせを行う第２のレジストレー
ション部とを設けたものである。
通常の第１の画像切り出し部、干渉処理部、平面位相除
去部、マルチルック処理部、アンラッピング処理部、デ
ータ変換部に、切り出した画像データの中の数カ所につ
いて２つの画像の位置合わせを行う第２のレジストレー
ション部とを設けたものである。
【００１６】また、第３の発明による信号処理装置は、
通常の第１の画像切り出し部、干渉処理部に、２つの画
像の位置合わせを行う粗レジストレーション部と、相関
係数を選択する係数選択部と相関値を計算する第１の相
関値演算部、相関値について比較を行う相関値比較部と
相関値から位置合わせを決定するオフセット値決定部と
で構成されるレジストレーション部を設けたものであ
る。
通常の第１の画像切り出し部、干渉処理部に、２つの画
像の位置合わせを行う粗レジストレーション部と、相関
係数を選択する係数選択部と相関値を計算する第１の相
関値演算部、相関値について比較を行う相関値比較部と
相関値から位置合わせを決定するオフセット値決定部と
で構成されるレジストレーション部を設けたものであ
る。
【００１７】また、第４の発明による信号処理装置は、
通常の第１の画像切り出し部、干渉処理部に、２つの画
像の位置合わせを行う粗レジストレーション部と、２つ
の画像データの相関値を算出する第２の相関値演算部と
相関値の値から相関値の最大値を検出する相関値最大値
検出部、相関値の最大値から位置合わせを決定するオフ
セット値決定部とで構成されるレジストレーション部を
設けたものである。
通常の第１の画像切り出し部、干渉処理部に、２つの画
像の位置合わせを行う粗レジストレーション部と、２つ
の画像データの相関値を算出する第２の相関値演算部と
相関値の値から相関値の最大値を検出する相関値最大値
検出部、相関値の最大値から位置合わせを決定するオフ
セット値決定部とで構成されるレジストレーション部を
設けたものである。
【００１８】また、第５の発明による信号処理装置は、
通常の第１の画像切り出し部、第１のレジストレーショ
ン部、干渉処理部、平面位相除去部、マルチルック処理
部、アンラッピング処理部データ変換部にノイズを除去
するフィルタ処理部とを設けたものである。
通常の第１の画像切り出し部、第１のレジストレーショ
ン部、干渉処理部、平面位相除去部、マルチルック処理
部、アンラッピング処理部データ変換部にノイズを除去
するフィルタ処理部とを設けたものである。
【００１９】
実施の形態１．図１はこの発明の実施の形態１を示す構
成ブロック図であり、図において図６の従来例と同一構
成の第１の画像データ１、第２の画像データ２、第１の
レジストレーション部４、干渉処理部５、平面位相除去
部６、マルチルック処理部７、アンラッピング処理部
８、データ変換部９については既に説明してあるのでこ
こでは説明を省略する。１０は画像分割部、１１は分散
演算部、１２は処理領域判定部、１３は第２の画像切り
出し部である。
成ブロック図であり、図において図６の従来例と同一構
成の第１の画像データ１、第２の画像データ２、第１の
レジストレーション部４、干渉処理部５、平面位相除去
部６、マルチルック処理部７、アンラッピング処理部
８、データ変換部９については既に説明してあるのでこ
こでは説明を省略する。１０は画像分割部、１１は分散
演算部、１２は処理領域判定部、１３は第２の画像切り
出し部である。
【００２０】次に動作について説明する。まず、画像分
割部１０では第１の画像データ１について後段の処理が
行いやすいサイズ（例えばレンジ方向１０２４点、アジ
マス方向１０２４点）に画像をｍ個に分割する。そして
分散演算部１１にてそれぞれの分割画像データの振幅分
散値を数３により求める。
割部１０では第１の画像データ１について後段の処理が
行いやすいサイズ（例えばレンジ方向１０２４点、アジ
マス方向１０２４点）に画像をｍ個に分割する。そして
分散演算部１１にてそれぞれの分割画像データの振幅分
散値を数３により求める。
【００２１】
【数３】
【００２２】数３においてＶL は分割画像L （L ＝１，
２，３，・・・ｍ）の振幅分散値、ＭL は画像複素デー
タの振幅値、｜ｆi,j ｜はレンジ座標ｉ、アジマス座標
ｊの画像複素データの振幅値である。一般に海面などの
標高の変化の小さな部分では、振幅分散は小さい。そこ
で処理領域判定部１２にて後段のレジストレーションを
行う領域を数４により判定する。
２，３，・・・ｍ）の振幅分散値、ＭL は画像複素デー
タの振幅値、｜ｆi,j ｜はレンジ座標ｉ、アジマス座標
ｊの画像複素データの振幅値である。一般に海面などの
標高の変化の小さな部分では、振幅分散は小さい。そこ
で処理領域判定部１２にて後段のレジストレーションを
行う領域を数４により判定する。
【００２３】
【数４】
【００２４】数４においてｋはスレッショルドである。
そしてスレッショルドｋを越す領域のみ後段の処理を行
うと判定する。そして処理が可能な領域から画像切り出
し部１３にて画像を切り出す。また、第２の画像データ
についても同様に切り出す。これらの切り出した２つの
画像データについて第１のレジストレーション部４にて
位置合わせを行う。後段の干渉処理部５、平面位相除去
部６、マルチルック処理部７、アンラッピング処理部
８、データ変換部９については従来装置と同一のもので
ある。
そしてスレッショルドｋを越す領域のみ後段の処理を行
うと判定する。そして処理が可能な領域から画像切り出
し部１３にて画像を切り出す。また、第２の画像データ
についても同様に切り出す。これらの切り出した２つの
画像データについて第１のレジストレーション部４にて
位置合わせを行う。後段の干渉処理部５、平面位相除去
部６、マルチルック処理部７、アンラッピング処理部
８、データ変換部９については従来装置と同一のもので
ある。
【００２５】実施の形態２．図２はこの発明の実施の形
態２を示す構成ブロック図であり、図において第１の画
像切り出し部３、干渉処理部５、平面位相除去部６、マ
ルチルック処理部７、アンラッピング処理部８、データ
変換部９は従来装置と同一のものである。１４は第２の
レジストレーション部である。
態２を示す構成ブロック図であり、図において第１の画
像切り出し部３、干渉処理部５、平面位相除去部６、マ
ルチルック処理部７、アンラッピング処理部８、データ
変換部９は従来装置と同一のものである。１４は第２の
レジストレーション部である。
【００２６】次に動作について説明する。第２のレジス
トレーション部１４では、切り出した画像データの数カ
所を使用してレジストレーションを行う。第１のレジス
トレーション部４では、通常、粗レジストレーション及
び精レジストレーションを行うがこの演算量を比較する
と次のようになる。例えば粗レジストレーションではＳ
／Ｎ比（信号対ノズルの比率）によりピクセル単位での
レジストレーションを行い、レンジ方向、アジマス方向
のずれ量であるオフセット値を求め、精レジストレーシ
ョンでは相関値によりサブピクセル単位（０．１ピクセ
ル）で、レンジ方向、アジマス方向のずれ量及びレンジ
伸縮率であるオフセット値を求めるとする。演算量の比
は数５のように表される。
トレーション部１４では、切り出した画像データの数カ
所を使用してレジストレーションを行う。第１のレジス
トレーション部４では、通常、粗レジストレーション及
び精レジストレーションを行うがこの演算量を比較する
と次のようになる。例えば粗レジストレーションではＳ
／Ｎ比（信号対ノズルの比率）によりピクセル単位での
レジストレーションを行い、レンジ方向、アジマス方向
のずれ量であるオフセット値を求め、精レジストレーシ
ョンでは相関値によりサブピクセル単位（０．１ピクセ
ル）で、レンジ方向、アジマス方向のずれ量及びレンジ
伸縮率であるオフセット値を求めるとする。演算量の比
は数５のように表される。
【００２７】
【数５】
【００２８】数５でＰは演算量の比、Ｍ，Ｎは切り出し
画像データのサイズでそれぞれレンジ、アジマスのサイ
ズ、Ｍ’，Ｎ’はレジストレーション用画像データのサ
イズでそれぞれレンジ、アジマスのサイズ、ｎはレジス
トレーション用画像データの個数である。例えばＭ＝１
０２４，Ｎ＝１０２４，Ｍ’＝８，Ｎ’＝８，ｎ＝５で
あればＰ≒１／３２７６であり、演算量がかなり減少で
きると言える。さらに切り出し画像データがＬ倍になっ
てもレジストレーション用画像データの個数ｎを増やす
だけなので演算量はＬ倍より小さくなる。
画像データのサイズでそれぞれレンジ、アジマスのサイ
ズ、Ｍ’，Ｎ’はレジストレーション用画像データのサ
イズでそれぞれレンジ、アジマスのサイズ、ｎはレジス
トレーション用画像データの個数である。例えばＭ＝１
０２４，Ｎ＝１０２４，Ｍ’＝８，Ｎ’＝８，ｎ＝５で
あればＰ≒１／３２７６であり、演算量がかなり減少で
きると言える。さらに切り出し画像データがＬ倍になっ
てもレジストレーション用画像データの個数ｎを増やす
だけなので演算量はＬ倍より小さくなる。
【００２９】実施の形態３．図３はこの発明の実施の形
態３を示す構成ブロック図であり、図において第１の画
像切り出し部３、干渉処理部５は従来装置と同一のもの
である。１５は粗レジストレーション部、１６は係数選
択部、１７は第１の相関値演算部、１８は相関値比較
部、１９はオフセット値決定部である。
態３を示す構成ブロック図であり、図において第１の画
像切り出し部３、干渉処理部５は従来装置と同一のもの
である。１５は粗レジストレーション部、１６は係数選
択部、１７は第１の相関値演算部、１８は相関値比較
部、１９はオフセット値決定部である。
【００３０】次に動作について説明する。粗レジストレ
ーション部１５では、おおまかな（例えば、ピクセル単
位）でのレジストレーションを行い、レンジ方向、アジ
マス方向のずれ量であるオフセット値を求める。係数選
択部１６ではレンジ、アジマス、レンジ伸縮の３つのオ
フセット値を１つの関数におく。例えば関数として数６
を選び、これを満たす係数を選択する。
ーション部１５では、おおまかな（例えば、ピクセル単
位）でのレジストレーションを行い、レンジ方向、アジ
マス方向のずれ量であるオフセット値を求める。係数選
択部１６ではレンジ、アジマス、レンジ伸縮の３つのオ
フセット値を１つの関数におく。例えば関数として数６
を選び、これを満たす係数を選択する。
【００３１】
【数６】
【００３２】数６においてＯ（ｊ）は、アジマスｊでの
オフセット値、Ｏr ，ＯAzはそれぞれレンジ方向、アジ
マス方向オフセット値、ａ〜ｄは係数である。そして第
１の相関値演算部１７にて相関値を演算し、相関値比較
部１８にて相関値を比較する。相関値が小さい場合は、
係数選択部１６に戻り、新たな係数を選択する。そして
相関値が極大となったところで、第２のオフセット値決
定部１９にてオフセット値を決定する。係数選択部１６
では、これによりレンジ、アジマス、レンジ伸縮の３点
についてまとめて検討することができる。
オフセット値、Ｏr ，ＯAzはそれぞれレンジ方向、アジ
マス方向オフセット値、ａ〜ｄは係数である。そして第
１の相関値演算部１７にて相関値を演算し、相関値比較
部１８にて相関値を比較する。相関値が小さい場合は、
係数選択部１６に戻り、新たな係数を選択する。そして
相関値が極大となったところで、第２のオフセット値決
定部１９にてオフセット値を決定する。係数選択部１６
では、これによりレンジ、アジマス、レンジ伸縮の３点
についてまとめて検討することができる。
【００３３】実施の形態４．図４はこの発明の実施の形
態４を示す構成ブロック図であり、図においては第１の
画像切り出し部３、干渉処理部５は従来装置と同一のも
のである。１５，１９〜２１はレジストレーション部４
の詳細を示しており、１５はピクセル単位でレジストレ
ーションを行う粗レジストレーション部、２０は相関値
を演算する第２の相関値演算部、２１は相関値最大検出
部、１９はオフセット値決定部である。
態４を示す構成ブロック図であり、図においては第１の
画像切り出し部３、干渉処理部５は従来装置と同一のも
のである。１５，１９〜２１はレジストレーション部４
の詳細を示しており、１５はピクセル単位でレジストレ
ーションを行う粗レジストレーション部、２０は相関値
を演算する第２の相関値演算部、２１は相関値最大検出
部、１９はオフセット値決定部である。
【００３４】次に動作について説明する。粗レジストレ
ーション部１５は実施例３にて説明したのでここでは省
略する。まず、粗レジストレーション部１５にて大まか
なオフセット値が求められた後、第２の相関値演算部２
０にて相関値が演算される。相関値εは数７に従って計
算する。
ーション部１５は実施例３にて説明したのでここでは省
略する。まず、粗レジストレーション部１５にて大まか
なオフセット値が求められた後、第２の相関値演算部２
０にて相関値が演算される。相関値εは数７に従って計
算する。
【００３５】
【数７】
【００３６】数７においてｆ１i,j はレンジｉ、アジマ
スｊの第１の画像データ１、ｆ２i, j はレンジｉ、アジ
マスｊの第２の画像データ２を表し、それぞれの画像デ
ータのレンジ、アジマス方向のデータの変化量の相関を
求めている。そして、相関値最大検出部２１において相
関値εの最大値を求める。これは、同じ位置にあればデ
ータの変化量が同じであり、このときの相関値は最大に
なることを意味している。オフセット値決定部１９にて
相関値εが最大となる位置をオフセット値と決定する。
スｊの第１の画像データ１、ｆ２i, j はレンジｉ、アジ
マスｊの第２の画像データ２を表し、それぞれの画像デ
ータのレンジ、アジマス方向のデータの変化量の相関を
求めている。そして、相関値最大検出部２１において相
関値εの最大値を求める。これは、同じ位置にあればデ
ータの変化量が同じであり、このときの相関値は最大に
なることを意味している。オフセット値決定部１９にて
相関値εが最大となる位置をオフセット値と決定する。
【００３７】実施の形態５．図５はこの発明の実施の形
態５を示す構成ブロック図であり、図において第１の画
像切り出し部３、第１のレジストレーション部４、干渉
処理部５、平面位相除去部６、マルチルック処理部７、
アンラッピング処理部８、データ変換部９は従来装置と
同一のものである。２２はフィルタ処理部である。
態５を示す構成ブロック図であり、図において第１の画
像切り出し部３、第１のレジストレーション部４、干渉
処理部５、平面位相除去部６、マルチルック処理部７、
アンラッピング処理部８、データ変換部９は従来装置と
同一のものである。２２はフィルタ処理部である。
【００３８】次に動作について説明する。フィルタ処理
部２２ではアンラッピング処理部８の出力データについ
て相関値があるスレッショルドより低い箇所について周
辺サブピクセルの位相データの平均値と置き換える。こ
れにより、位相ノイズを低減できる。また、このフィル
タ処理部２２はアンラッピング処理部８の後段ではな
く、マルチルック処理部７の後段にも再構成可能とす
る。これはマルチルック処理部７はノイズ低減のために
行っているものであるが、これだけではまだ位相ノイズ
が大きいことがあり、位相ノイズが大きいとアンラッピ
ング処理がうまくできないことがあるので、この場合、
フィルタ処理部２２をマルチルック処理部７の後段にて
行う。
部２２ではアンラッピング処理部８の出力データについ
て相関値があるスレッショルドより低い箇所について周
辺サブピクセルの位相データの平均値と置き換える。こ
れにより、位相ノイズを低減できる。また、このフィル
タ処理部２２はアンラッピング処理部８の後段ではな
く、マルチルック処理部７の後段にも再構成可能とす
る。これはマルチルック処理部７はノイズ低減のために
行っているものであるが、これだけではまだ位相ノイズ
が大きいことがあり、位相ノイズが大きいとアンラッピ
ング処理がうまくできないことがあるので、この場合、
フィルタ処理部２２をマルチルック処理部７の後段にて
行う。
【００３９】
【発明の効果】第１の発明によればレジストレーション
を行う前に処理領域の判定を行うので、標高の変化の小
さい部分ではレジストレーションを行わないため、曖昧
なレジストレーションを避けることができる。
を行う前に処理領域の判定を行うので、標高の変化の小
さい部分ではレジストレーションを行わないため、曖昧
なレジストレーションを避けることができる。
【００４０】また、第２の発明によれば切り出し画像デ
ータの数カ所にてのみレジストレーションを行うため、
レジストレーションの演算量を低減することができる。
ータの数カ所にてのみレジストレーションを行うため、
レジストレーションの演算量を低減することができる。
【００４１】また、第３の発明によれば第１のレジスト
レーション部にてオフセット値を１つの関数としている
ので演算量が大幅に減少することができる。
レーション部にてオフセット値を１つの関数としている
ので演算量が大幅に減少することができる。
【００４２】また、第４の発明によれば第２の相関値を
演算するため、位相ノイズの比較的大きな画像に対し、
精度の高いレジストレーションが可能である。
演算するため、位相ノイズの比較的大きな画像に対し、
精度の高いレジストレーションが可能である。
【００４３】また、第５の発明によればフィルタ処理部
を付加することにより、位相ノイズを低減し、精度の高
い標高データを得ることができ、また、アンラッピング
処理が正しく行えるノイズレベルまで押さえることがで
きる。
を付加することにより、位相ノイズを低減し、精度の高
い標高データを得ることができ、また、アンラッピング
処理が正しく行えるノイズレベルまで押さえることがで
きる。
【図１】 この発明による信号処理装置の実施の形態１
を示す図である。
を示す図である。
【図２】 この発明による信号処理装置の実施の形態２
を示す図である。
を示す図である。
【図３】 この発明による信号処理装置の実施の形態３
を示す図である。
を示す図である。
【図４】 この発明による信号処理装置の実施の形態４
を示す図である。
を示す図である。
【図５】 この発明による信号処理装置の実施の形態５
を示す図である。
を示す図である。
【図６】 従来の信号処理装置の図である。
【図７】 インタフェロメトリのジオメトリを示す図で
ある。
ある。
１ 第１の画像データ、２ 第２の画像データ、３ 第
１の画像切り出し部、４ 第１のレジストレーション
部、５ 干渉処理部、６ 平面位相除去部、７マルチル
ック処理部、８ アンラッピング処理部、９ データ変
換部、１０ 画像分割部、１１ 分散演算部、１２ 処
理領域判定部、１３ 第２の画像切り出し部、１４ 第
２のレジストレーション部、１５ 粗レジストレーショ
ン部、１６ 係数選択部、１７ 第１の相関値演算部、
１８ 相関値比較部、１９ オフセット値決定部、２０
第２の相関値演算部、２１ 相関値最大検出部、２２
フィルタ処理部。
１の画像切り出し部、４ 第１のレジストレーション
部、５ 干渉処理部、６ 平面位相除去部、７マルチル
ック処理部、８ アンラッピング処理部、９ データ変
換部、１０ 画像分割部、１１ 分散演算部、１２ 処
理領域判定部、１３ 第２の画像切り出し部、１４ 第
２のレジストレーション部、１５ 粗レジストレーショ
ン部、１６ 係数選択部、１７ 第１の相関値演算部、
１８ 相関値比較部、１９ オフセット値決定部、２０
第２の相関値演算部、２１ 相関値最大検出部、２２
フィルタ処理部。
Claims (5)
- 【請求項１】 近接した２つの軌道で得られた２組のう
ちの１つの画像データから画像を分割する画像分割部、
この画像分割部の出力データである分割されたそれぞれ
の画像データについて分散を演算する分散演算部、この
分散演算部の出力データについて画像処理を行う領域を
判定する処理領域判定部、この処理領域判定部の出力と
画像データから画像を切り出す画像切り出し部、２つの
画像切り出し部の出力データより位置合わせを行う第１
のレジストレーション部、この第１のレジストレーショ
ン部の出力であるオフセット値に合わせて２つの画像デ
ータから位相差を求める干渉処理部、この干渉処理部の
出力である位相差から平面位相を除去する平面位相除去
部、この平面位相除去部の出力データを用いてマルチル
ック処理を行うマルチルック処理部、このマルチルック
処理部の出力データを積算処理するアンラッピング処理
部、このアンラッピング処理部の出力データである位相
差データを標高データに変換するデータ変換部から構成
されることを特徴とする信号処理装置。 - 【請求項２】 近接した２つの軌道で得られた２組の画
像データについてそれぞれ画像を切り出す画像切り出し
部、この画像切り出し部の出力データより位置合わせを
行う第１のレジストレーション部を付加したことを特徴
とする請求項１記載の信号処理装置。 - 【請求項３】 ２つの画像データの位置合わせを行う際
にピクセル単位でレジストレーションを行う粗レジスト
レーション部、この粗レジストレーション部の出力デー
タからレンジ方向オフセット値、アジマス方向オフセッ
ト値、レンジ収縮率を１つの関数として表し係数を選択
する係数選択部、この係数選択部の出力である係数を使
用してレジストレーションを行い、相関値を演算する第
１の相関値演算部、この第１の相関値演算部の出力であ
る相関値の最小を求める相関値比較部、この相関値比較
部の出力データより位置合わせを行うオフセット値を決
定するオフセット値決定部を付加したことを特徴とする
請求項１記載の信号処理装置。 - 【請求項４】 粗レジストレーション部の出力データと
２組の画像切り出し部の出力データにより相関値を演算
する第２の相関値演算部、この第２の相関値演算部の出
力である相関値が最大を求める相関値最大値検出部、こ
の相関値最大値検出部の出力データよりオフセット値を
求めるオフセット値決定部を付加したことを特徴とする
請求項１記載の信号処理装置。 - 【請求項５】 アンラッピング処理部の出力データの雑
音を除去するフィルタ処理部を付加したことを特徴とし
た請求項１記載の信号処理装置。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP8092710A JPH09281231A (ja) | 1996-04-15 | 1996-04-15 | 信号処理装置 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP8092710A JPH09281231A (ja) | 1996-04-15 | 1996-04-15 | 信号処理装置 |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| JPH09281231A true JPH09281231A (ja) | 1997-10-31 |
Family
ID=14062033
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP8092710A Pending JPH09281231A (ja) | 1996-04-15 | 1996-04-15 | 信号処理装置 |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JPH09281231A (ja) |
Cited By (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2004191053A (ja) * | 2002-12-06 | 2004-07-08 | Mitsubishi Electric Corp | 合成開口レーダ装置及び数値標高モデル作成方法 |
| JP2010048651A (ja) * | 2008-08-21 | 2010-03-04 | Mitsubishi Electric Corp | レーダ画像処理装置 |
| JP2015141124A (ja) * | 2014-01-29 | 2015-08-03 | 三菱電機株式会社 | 合成開口レーダ信号処理装置及び合成開口レーダ信号処理方法 |
| CN108802727A (zh) * | 2018-04-13 | 2018-11-13 | 长沙理工大学 | 一种顾及流变参数的时序InSAR公路形变监测模型及解算方法 |
| WO2019106850A1 (ja) * | 2017-12-01 | 2019-06-06 | 日本電気株式会社 | Ｓａｒ画像解析システム、画像処理装置、画像処理方法および画像処理プログラム |
-
1996
- 1996-04-15 JP JP8092710A patent/JPH09281231A/ja active Pending
Cited By (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2004191053A (ja) * | 2002-12-06 | 2004-07-08 | Mitsubishi Electric Corp | 合成開口レーダ装置及び数値標高モデル作成方法 |
| JP2010048651A (ja) * | 2008-08-21 | 2010-03-04 | Mitsubishi Electric Corp | レーダ画像処理装置 |
| JP2015141124A (ja) * | 2014-01-29 | 2015-08-03 | 三菱電機株式会社 | 合成開口レーダ信号処理装置及び合成開口レーダ信号処理方法 |
| WO2019106850A1 (ja) * | 2017-12-01 | 2019-06-06 | 日本電気株式会社 | Ｓａｒ画像解析システム、画像処理装置、画像処理方法および画像処理プログラム |
| JPWO2019106850A1 (ja) * | 2017-12-01 | 2020-12-03 | 日本電気株式会社 | Ｓａｒ画像解析システム、画像処理装置、画像処理方法および画像処理プログラム |
| US11709254B2 (en) | 2017-12-01 | 2023-07-25 | Nec Corporation | SAR image analysis system, image processing apparatus, image processing method, and image processing program |
| CN108802727A (zh) * | 2018-04-13 | 2018-11-13 | 长沙理工大学 | 一种顾及流变参数的时序InSAR公路形变监测模型及解算方法 |
| CN108802727B (zh) * | 2018-04-13 | 2021-01-15 | 长沙理工大学 | 一种顾及流变参数的时序InSAR公路形变监测模型及解算方法 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Bamler | Doppler frequency estimation and the Cramer-Rao bound | |
| Koch | Directional analysis of SAR images aiming at wind direction | |
| US11041941B2 (en) | Method and device for calibrating a radar object detection system | |
| US7961147B1 (en) | Long baseline phase interferometer ambiguity resolution using frequency differences | |
| US5923279A (en) | Method of correcting an object-dependent spectral shift in radar interferograms | |
| US6677885B1 (en) | Method for mitigating atmospheric propagation error in multiple pass interferometric synthetic aperture radar | |
| CN107918127A (zh) | 一种基于车载InSAR的道路边坡形变检测系统及方法 | |
| CN102713503A (zh) | 干涉测距仪中的散斑减轻方法和对应的测距仪 | |
| JP6292772B2 (ja) | レーダ画像処理装置及びレーダ画像処理方法 | |
| JP2017072473A (ja) | 画像処理装置および画像処理方法 | |
| CN117949131B (zh) | 一种拉索全场模态分析和索力识别方法及系统 | |
| US7116265B2 (en) | Recognition algorithm for the unknown target rejection based on shape statistics obtained from orthogonal distance function | |
| Lombardini et al. | Application of array processing techniques to multibaseline InSAR for layover solution | |
| JP2006003302A (ja) | レーダ画像処理装置 | |
| EP3012577B1 (en) | Method for measuring a deformation of a surface of an object, computer program, and computer readable recording medium | |
| US20010049709A1 (en) | Fringe analysis method using fourier transform | |
| US20010019635A1 (en) | Fringe analysis method using fourier transform | |
| CN106526594A (zh) | 一种ati‑sar海面单视复图像配准的方法 | |
| JP3743177B2 (ja) | レーダ信号処理装置 | |
| JP4827330B2 (ja) | レーダ波浪測定方法及び装置 | |
| Sahraoui et al. | Radar interferometry with Sarscape software | |
| JP2000162316A (ja) | レーダ信号処理装置 | |
| JP6278910B2 (ja) | レーダ画像処理装置及びレーダ画像処理方法 | |
| CN117784130A (zh) | 基于动目标检测和去除的InSAR干涉测高方法及装置 | |
| CN115239979B (zh) | 基于结构方向相似度的异构图像匹配方法 |