JP2001083243A - 干渉型合成開口レーダによる地形の３次元情報抽出装置 - Google Patents
干渉型合成開口レーダによる地形の３次元情報抽出装置Info
- Publication number
- JP2001083243A JP2001083243A JP25870399A JP25870399A JP2001083243A JP 2001083243 A JP2001083243 A JP 2001083243A JP 25870399 A JP25870399 A JP 25870399A JP 25870399 A JP25870399 A JP 25870399A JP 2001083243 A JP2001083243 A JP 2001083243A
- Authority
- JP
- Japan
- Prior art keywords
- filter
- interference
- size
- interference image
- noise
- 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
- Image Processing (AREA)
- Radar Systems Or Details Thereof (AREA)
- Image Analysis (AREA)
Abstract
(57)【要約】
【課題】 衛星や航空機などで３次元情報を抽出する干
渉型合成開口レーダにおいて、電波散乱特性等の影響に
より発生する位相ノイズを除去し、３次元情報抽出の精
度向上を図る。 【解決手段】 合成開口レーダで得られた観測条件の異
なる２つの観測データから、SAR画像再生処理によって
得られた２つの複素SAR画像データを互いに干渉させて
干渉画像を得て、上記干渉画像を地形に対応するように
位相アンラップ処理し、上記アンラップ処理後の位相デ
ータを地形の高さに変換し３次元情報を得る干渉型SAR
による３次元情報抽出装置において、上記干渉画像に存
在する位相ノイズをフィルタリング処理により除去する
ことで、ノイズ除去干渉画像を得るフィルタリング処理
手段を備えた。
渉型合成開口レーダにおいて、電波散乱特性等の影響に
より発生する位相ノイズを除去し、３次元情報抽出の精
度向上を図る。 【解決手段】 合成開口レーダで得られた観測条件の異
なる２つの観測データから、SAR画像再生処理によって
得られた２つの複素SAR画像データを互いに干渉させて
干渉画像を得て、上記干渉画像を地形に対応するように
位相アンラップ処理し、上記アンラップ処理後の位相デ
ータを地形の高さに変換し３次元情報を得る干渉型SAR
による３次元情報抽出装置において、上記干渉画像に存
在する位相ノイズをフィルタリング処理により除去する
ことで、ノイズ除去干渉画像を得るフィルタリング処理
手段を備えた。
Description
【発明の詳細な説明】
【０００１】
【発明の属する技術分野】この発明は、主として等高線
地図作成を目的としたシステムにおいて、人工衛星や航
空機等の飛翔体に搭載された合成開口レーダ（Syntheti
c Aperture Radar; SAR）で観測したデータから地形の
３次元情報を抽出する装置に関するものである。
地図作成を目的としたシステムにおいて、人工衛星や航
空機等の飛翔体に搭載された合成開口レーダ（Syntheti
c Aperture Radar; SAR）で観測したデータから地形の
３次元情報を抽出する装置に関するものである。
【０００２】
【従来の技術】地形の３次元情報を得る技術として、光
学センサによる航空測量やステレオ視等がある。しか
し、上記の技術は雲天時／夜間に観測することが不可能
であった。そこで、近年、天候や観測時間に依存しない
干渉型SARを用いた技術が提案されている。
学センサによる航空測量やステレオ視等がある。しか
し、上記の技術は雲天時／夜間に観測することが不可能
であった。そこで、近年、天候や観測時間に依存しない
干渉型SARを用いた技術が提案されている。
【０００３】図１７は飛翔体を人工衛星としたときの干
渉型SARを用いて地形の３次元情報を抽出するシステム
の概念図である。１a、１ｂはSARを搭載した軌道の異な
る人工衛星である。２は上記SARで得られた観測データ
から地形の３次元情報を得るシステムである。
渉型SARを用いて地形の３次元情報を抽出するシステム
の概念図である。１a、１ｂはSARを搭載した軌道の異な
る人工衛星である。２は上記SARで得られた観測データ
から地形の３次元情報を得るシステムである。
【０００４】図１７で示すように、干渉型SARを用いた
技術では、異なる観測条件で得られた２つのSAR画像を
干渉させるが、電波散乱特性や軌道間距離等の影響で２
つのSAR画像の相関が十分に得られないと、干渉画像に
位相ノイズが発生し精度が劣化する。これについて、い
くつかの手段が提案されている。
技術では、異なる観測条件で得られた２つのSAR画像を
干渉させるが、電波散乱特性や軌道間距離等の影響で２
つのSAR画像の相関が十分に得られないと、干渉画像に
位相ノイズが発生し精度が劣化する。これについて、い
くつかの手段が提案されている。
【０００５】図１８、図１９は従来の干渉型SARによる
地形の３次元情報抽出装置である。図１８において、３
は送受信部、４はSAR画像再生処理を行う画像処理部、
５は高速フーリエ変換（FFT）部、６は周波数カットを行
うフィルタ部、７は逆高速フーリエ変換（IFFT）部、８
は干渉画像を得る干渉処理部、９は干渉画像の差分をと
る差分抽出部である。図１８では送受信部３で複数回の
観測により、画像処理部４で少なくとも２組以上のSAR
画像を再生し、干渉処理部８で干渉画像を作成し、干渉
画像の差分を差分抽出部９でとって３次元画像データを
得る。また、フィルタ部６でSAR再生画像を周波数空間
で相関処理を行い、３次元情報の精度向上を図るもので
ある。
地形の３次元情報抽出装置である。図１８において、３
は送受信部、４はSAR画像再生処理を行う画像処理部、
５は高速フーリエ変換（FFT）部、６は周波数カットを行
うフィルタ部、７は逆高速フーリエ変換（IFFT）部、８
は干渉画像を得る干渉処理部、９は干渉画像の差分をと
る差分抽出部である。図１８では送受信部３で複数回の
観測により、画像処理部４で少なくとも２組以上のSAR
画像を再生し、干渉処理部８で干渉画像を作成し、干渉
画像の差分を差分抽出部９でとって３次元画像データを
得る。また、フィルタ部６でSAR再生画像を周波数空間
で相関処理を行い、３次元情報の精度向上を図るもので
ある。
【０００６】次に、図１９について説明する。図１９に
おいて、１０は送受信部、１１はSAR画像再生処理を行
う画像処理部、１２は干渉画像を得る干渉処理部、１３
はマルチルック処理部、１４は位相を地形に対応させる
位相アンラップ処理部、１５は地形に対応した位相を地
形の高さに変換する変換処理部である。図１９では、干
渉画像をマルチルック処理部１３でマルチルック処理す
ることにより位相ノイズを低減し、３次元情報の精度向
上を図るものである。
おいて、１０は送受信部、１１はSAR画像再生処理を行
う画像処理部、１２は干渉画像を得る干渉処理部、１３
はマルチルック処理部、１４は位相を地形に対応させる
位相アンラップ処理部、１５は地形に対応した位相を地
形の高さに変換する変換処理部である。図１９では、干
渉画像をマルチルック処理部１３でマルチルック処理す
ることにより位相ノイズを低減し、３次元情報の精度向
上を図るものである。
【０００７】
【発明が解決しようとする課題】図１８の場合では、SA
R画像を２組以上必要とするため、複数回の観測が必要
となり、SAR画像再生処理の回数も増える。さらに、 FF
T部５でSAR画像を周波数空間で相関処理を行うため、FF
TおよびIFFTの処理を必要とし、処理時間が増す。ま
た、周波数カットを行うと、分解能が低下するため３次
元情報の精度は向上しない。図１９の場合では、干渉画
像にマルチルック処理を行うため、図１８の場合と同様
に、分解能が低下する。
R画像を２組以上必要とするため、複数回の観測が必要
となり、SAR画像再生処理の回数も増える。さらに、 FF
T部５でSAR画像を周波数空間で相関処理を行うため、FF
TおよびIFFTの処理を必要とし、処理時間が増す。ま
た、周波数カットを行うと、分解能が低下するため３次
元情報の精度は向上しない。図１９の場合では、干渉画
像にマルチルック処理を行うため、図１８の場合と同様
に、分解能が低下する。
【０００８】この発明は、上記のような課題を解決する
ためになされたものであり、１組のSAR画像から分解能
を低下させることなく、短時間で高精度な地形の３次元
情報を容易に得ることができる、干渉型SARによる３次
元情報抽出装置を提供することを目的とするものであ
る。
ためになされたものであり、１組のSAR画像から分解能
を低下させることなく、短時間で高精度な地形の３次元
情報を容易に得ることができる、干渉型SARによる３次
元情報抽出装置を提供することを目的とするものであ
る。
【０００９】
【課題を解決するための手段】第１の発明の干渉型SAR
による地形の３次元情報抽出装置は、SARで得られた観
測条件の異なる２つの観測データから、SAR画像再生処
理によって得られた２つの複素SAR画像データを互いに
干渉させて干渉画像を得る手段と、上記干渉画像を地形
に対応するように位相アンラップ処理し、上記アンラッ
プ処理後の位相データを地形の高さに変換し３次元情報
を得る手段とを備えた干渉型SARによる３次元情報抽出
装置において、上記干渉画像に存在する位相ノイズをフ
ィルタリング処理により除去し、ノイズ除去干渉画像を
得るフィルタリング処理手段を備えたものである。
による地形の３次元情報抽出装置は、SARで得られた観
測条件の異なる２つの観測データから、SAR画像再生処
理によって得られた２つの複素SAR画像データを互いに
干渉させて干渉画像を得る手段と、上記干渉画像を地形
に対応するように位相アンラップ処理し、上記アンラッ
プ処理後の位相データを地形の高さに変換し３次元情報
を得る手段とを備えた干渉型SARによる３次元情報抽出
装置において、上記干渉画像に存在する位相ノイズをフ
ィルタリング処理により除去し、ノイズ除去干渉画像を
得るフィルタリング処理手段を備えたものである。
【００１０】また 第２の発明の干渉型SARによる地形の
３次元情報抽出装置は、上記フィルタリング処理手段
に、干渉画像におけるフィルタ領域の位相が不連続か否
かを判定する第１の手段と、上記第１の手段で不連続と
判定された領域について正の値をもつ画素数が多い領域
か負の画素数が多い領域かを判定する第２の手段と、上
記第２の手段において正の値をもつ画素数が多い領域の
平均値を算出する第３の手段とを備えたものである。
３次元情報抽出装置は、上記フィルタリング処理手段
に、干渉画像におけるフィルタ領域の位相が不連続か否
かを判定する第１の手段と、上記第１の手段で不連続と
判定された領域について正の値をもつ画素数が多い領域
か負の画素数が多い領域かを判定する第２の手段と、上
記第２の手段において正の値をもつ画素数が多い領域の
平均値を算出する第３の手段とを備えたものである。
【００１１】第３の発明の干渉型SARによる地形の３次
元情報抽出装置は、上記フィルタリング処理手段に、上
記第１の手段と、上記第２の手段と、上記第２の手段で
負の値をもつ画素数が多い領域の平均値を算出する第４
の手段とを備えたものである。
元情報抽出装置は、上記フィルタリング処理手段に、上
記第１の手段と、上記第２の手段と、上記第２の手段で
負の値をもつ画素数が多い領域の平均値を算出する第４
の手段とを備えたものである。
【００１２】また 第４の発明の干渉型SARによる地形の
３次元情報抽出装置は、上記フィルタリング処理手段
に、上記第１の手段と、上記第１の手段で不連続でない
と判定された領域の平均値を算出する第５の手段とを備
えたものである。
３次元情報抽出装置は、上記フィルタリング処理手段
に、上記第１の手段と、上記第１の手段で不連続でない
と判定された領域の平均値を算出する第５の手段とを備
えたものである。
【００１３】第５の発明の干渉型SARによる地形の３次
元情報抽出装置は、上記フィルタリング処理手段に、フ
ィルタの大きさをユーザが任意に設定する第６の手段を
備えたものである。
元情報抽出装置は、上記フィルタリング処理手段に、フ
ィルタの大きさをユーザが任意に設定する第６の手段を
備えたものである。
【００１４】また 第６の発明の干渉型SARによる３次元
情報抽出装置は、上記フィルタリング処理手段に、干渉
画像とノイズ除去干渉画像を表示し、ユーザがノイズ除
去５態を評価してフィルタの大きさを設定する第７の手
段を備えたものである。
情報抽出装置は、上記フィルタリング処理手段に、干渉
画像とノイズ除去干渉画像を表示し、ユーザがノイズ除
去５態を評価してフィルタの大きさを設定する第７の手
段を備えたものである。
【００１５】第７の発明の干渉型SARによる地形の３次
元情報抽出装置は、上記フィルタリング処理手段に、干
渉画像のノイズ量およびノイズ除去干渉画像のノイズ量
を表示し、ユーザがノイズ除去状態を評価してフィルタ
の大きさを設定する第８の手段を備えたものである。
元情報抽出装置は、上記フィルタリング処理手段に、干
渉画像のノイズ量およびノイズ除去干渉画像のノイズ量
を表示し、ユーザがノイズ除去状態を評価してフィルタ
の大きさを設定する第８の手段を備えたものである。
【００１６】また 第８の発明の干渉型SARによる地形の
３次元情報抽出装置は、上記フィルタリング処理手段
に、過去のフィルタリング処理において使用したフィル
タの大きさと撮像地域を記録しておき、この記録から同
一撮像地域のフィルタリング処理のときに効果の高いフ
ィルタの大きさを設定する第９の手段を備えたものであ
る。
３次元情報抽出装置は、上記フィルタリング処理手段
に、過去のフィルタリング処理において使用したフィル
タの大きさと撮像地域を記録しておき、この記録から同
一撮像地域のフィルタリング処理のときに効果の高いフ
ィルタの大きさを設定する第９の手段を備えたものであ
る。
【００１７】第９の発明の干渉型SARによる地形の３次
元情報抽出装置は、上記フィルタリング処理手段に、過
去のフィルタリング処理において使用したフィルタの大
きさと干渉画像の位相ノイズ量を記録しておき、この記
録から干渉画像の位相ノイズ量に対応したフィルタの大
きさを設定する第１０の手段を備えたものである。
元情報抽出装置は、上記フィルタリング処理手段に、過
去のフィルタリング処理において使用したフィルタの大
きさと干渉画像の位相ノイズ量を記録しておき、この記
録から干渉画像の位相ノイズ量に対応したフィルタの大
きさを設定する第１０の手段を備えたものである。
【００１８】また 第１０発明の干渉型SARによる地形の
３次元情報抽出装置は、上記フィルタリング処理手段
に、過去のフィルタリング処理において使用した干渉画
像のセンサ名とマルチルック数とフィルタの大きさを記
録しておき、この記録から同一撮像状況（ピクセルスペ
ーシングと波長）に対応したフィルタの大きさを設定す
る第１１の手段を備えたものである。
３次元情報抽出装置は、上記フィルタリング処理手段
に、過去のフィルタリング処理において使用した干渉画
像のセンサ名とマルチルック数とフィルタの大きさを記
録しておき、この記録から同一撮像状況（ピクセルスペ
ーシングと波長）に対応したフィルタの大きさを設定す
る第１１の手段を備えたものである。
【００１９】第１１発明の干渉型SARによる地形の３次
元情報抽出装置は、上記フィルタリング処理手段に、フ
ィルタの大きさを任意に設定する第６の手段と、干渉画
像のノイズ量およびノイズ除去干渉画像のノイズ量を表
示し、ユーザがノイズ除去状態を評価してフィルタの大
きさを調整する第８の手段と、過去のフィルタリング処
理において使用したフィルタの大きさと撮像地域を記録
しておき、この記録から同一撮像地域のフィルタリング
処理のときに効果の高いフィルタの大きさを設定する第
９の手段と、過去のフィルタリング処理において使用し
たフィルタの大きさと干渉画像の位相ノイズ量を記録し
ておき、この記録から干渉画像の位相ノイズ量に対応し
たフィルタの大きさを設定する第 １０の手段と、過去
のフィルタリング処理において使用した干渉画像のセン
サ名と画像のマルチルック数とフィルタの大きさを記録
しておき、この記録からに対応したフィルタの大きさを
設定する第１１の手段と、干渉画像に対する位相ノイズ
除去の１回目の試行において、上記第６の手段、第９の
手段、第１０の手段または第１１の手段のいずれかを選
択する第１２の手段と、上記干渉画像に対する位相ノイ
ズ除去の２回目以降の複数回の再試行において、第６の
手段、第７の手段、第８の手段、第９の手段、第１０の
手段または第１１の手段のいずれかを再試行の度に選択
する第１３の手段とを備えたものである。
元情報抽出装置は、上記フィルタリング処理手段に、フ
ィルタの大きさを任意に設定する第６の手段と、干渉画
像のノイズ量およびノイズ除去干渉画像のノイズ量を表
示し、ユーザがノイズ除去状態を評価してフィルタの大
きさを調整する第８の手段と、過去のフィルタリング処
理において使用したフィルタの大きさと撮像地域を記録
しておき、この記録から同一撮像地域のフィルタリング
処理のときに効果の高いフィルタの大きさを設定する第
９の手段と、過去のフィルタリング処理において使用し
たフィルタの大きさと干渉画像の位相ノイズ量を記録し
ておき、この記録から干渉画像の位相ノイズ量に対応し
たフィルタの大きさを設定する第 １０の手段と、過去
のフィルタリング処理において使用した干渉画像のセン
サ名と画像のマルチルック数とフィルタの大きさを記録
しておき、この記録からに対応したフィルタの大きさを
設定する第１１の手段と、干渉画像に対する位相ノイズ
除去の１回目の試行において、上記第６の手段、第９の
手段、第１０の手段または第１１の手段のいずれかを選
択する第１２の手段と、上記干渉画像に対する位相ノイ
ズ除去の２回目以降の複数回の再試行において、第６の
手段、第７の手段、第８の手段、第９の手段、第１０の
手段または第１１の手段のいずれかを再試行の度に選択
する第１３の手段とを備えたものである。
【００２０】
【発明の実施の形態】実施の形態１．この発明の実施の
形態１について図を用いて説明する。 図１はこの発明
の干渉型合成開口レーダによる３次元情報抽出装置の構
成を示す。１６a、１６ｂはSARアンテナの送受信部、１
７は画像処理部、１８は干渉処理部、１９はフィルタリ
ング処理手段、２０は位相アンラップ処理部、２１は変
換処理部、２２は最終的に出力される３次元情報であ
る。
形態１について図を用いて説明する。 図１はこの発明
の干渉型合成開口レーダによる３次元情報抽出装置の構
成を示す。１６a、１６ｂはSARアンテナの送受信部、１
７は画像処理部、１８は干渉処理部、１９はフィルタリ
ング処理手段、２０は位相アンラップ処理部、２１は変
換処理部、２２は最終的に出力される３次元情報であ
る。
【００２１】はじめに、軌道の異なる２つのアンテナま
たは、１つのアンテナで２回の観測を行い、SARアンテ
ナの送受信部１６a、１６ｂで２つのSARデータを得る。
図１では２つのSARアンテナの送受信部を示した。
たは、１つのアンテナで２回の観測を行い、SARアンテ
ナの送受信部１６a、１６ｂで２つのSARデータを得る。
図１では２つのSARアンテナの送受信部を示した。
【００２２】次に、上記２つのSARデータを画像処理部
１７でSAR画像再生処理を行い、２つの複素SAR画像を得
る。さらに、干渉処理部１８で、上記２つの複素SAR画
像を干渉させて、干渉画像を得る。
１７でSAR画像再生処理を行い、２つの複素SAR画像を得
る。さらに、干渉処理部１８で、上記２つの複素SAR画
像を干渉させて、干渉画像を得る。
【００２３】ここで、フィルタリング処理手段１９によ
り、干渉画像に存在する位相ノイズを除去する。図２は
フィルタリング処理手段１９により位相ノイズが除去さ
れる動作を示す。２３は干渉画像、２４はフィルタリン
グ処理手段のフィルタ領域であり、フィルタ領域２４が
干渉画像２３上を１画素ずつ移動して、位相ノイズを除
去する。
り、干渉画像に存在する位相ノイズを除去する。図２は
フィルタリング処理手段１９により位相ノイズが除去さ
れる動作を示す。２３は干渉画像、２４はフィルタリン
グ処理手段のフィルタ領域であり、フィルタ領域２４が
干渉画像２３上を１画素ずつ移動して、位相ノイズを除
去する。
【００２４】図３はフィルタ領域２４が干渉画像２３上
で、位相が−πからπの２πの範囲において、位相が不
連続であり、正の画素が多い領域に位置するときの動作
を示す。ただし、フィルタは５×５画素の矩形で、−π
からπは−１から１で表示している。２５は正の値をも
つ画素の領域であり、２６は負の値をもつ画素の領域で
ある。また、図４は図３の場合のフィルタリング処理手
段の構成図であり、２７はフィルタリング処理手段前の
干渉画像、２８は第１の手段、２９は第２の手段、 ３
０は第３の手段、３１はノイズ除去画像の画素、３２は
ノイズ除去画像である。
で、位相が−πからπの２πの範囲において、位相が不
連続であり、正の画素が多い領域に位置するときの動作
を示す。ただし、フィルタは５×５画素の矩形で、−π
からπは−１から１で表示している。２５は正の値をも
つ画素の領域であり、２６は負の値をもつ画素の領域で
ある。また、図４は図３の場合のフィルタリング処理手
段の構成図であり、２７はフィルタリング処理手段前の
干渉画像、２８は第１の手段、２９は第２の手段、 ３
０は第３の手段、３１はノイズ除去画像の画素、３２は
ノイズ除去画像である。
【００２５】まず、第１の手段２８により上記領域が不
連続か否かを判定する。判定基準は、正の値をもつ画素
値の平均値をp_avg、負の値をもつ画素値の平均値をm_a
vgとすると、 p_avg - m_avg > π （１） のとき、不連続であり、上記式（１）でないとき、不連
続でないとする。図３の場合では、 p_avg = 0.8、m_avg = -0.8 であり、上記式（１）を満たすから不連続と判定され
る。
連続か否かを判定する。判定基準は、正の値をもつ画素
値の平均値をp_avg、負の値をもつ画素値の平均値をm_a
vgとすると、 p_avg - m_avg > π （１） のとき、不連続であり、上記式（１）でないとき、不連
続でないとする。図３の場合では、 p_avg = 0.8、m_avg = -0.8 であり、上記式（１）を満たすから不連続と判定され
る。
【００２６】第１の手段２８により不連続として判断さ
れると、第２の手段２９により正の値をもつ画素数が多
いか、負の値をもつ画素数が多いかを判定する。判定基
準は、正の画素数をp_num、負の画素数をm_numとする
と、 p_num > m_num （２） のとき、正の画素数が多い領域であり、上記式（２）で
ないとき、負の画素数が多い領域とする。図３の場合で
は、 p_num = 15、m_num = 10 であり、上記式（２）を満たすから正の画素数が多いと
判定される。
れると、第２の手段２９により正の値をもつ画素数が多
いか、負の値をもつ画素数が多いかを判定する。判定基
準は、正の画素数をp_num、負の画素数をm_numとする
と、 p_num > m_num （２） のとき、正の画素数が多い領域であり、上記式（２）で
ないとき、負の画素数が多い領域とする。図３の場合で
は、 p_num = 15、m_num = 10 であり、上記式（２）を満たすから正の画素数が多いと
判定される。
【００２７】第２の手段２９により正の画素数が多いと
判定されると、第３の手段３０により正の画素数が多い
フィルタ領域の平均値を算出する。正の画素数が多いフ
ィルタ領域の平均値をpa_avgとすると、 pa_avg = （p_avg + m_avg + 2π）／2 （３） によって算出される。上記平均値pa_avgをノイズ除去干
渉画像の画素値３１とする。
判定されると、第３の手段３０により正の画素数が多い
フィルタ領域の平均値を算出する。正の画素数が多いフ
ィルタ領域の平均値をpa_avgとすると、 pa_avg = （p_avg + m_avg + 2π）／2 （３） によって算出される。上記平均値pa_avgをノイズ除去干
渉画像の画素値３１とする。
【００２８】図５はフィルタ領域２４が干渉画像２３上
で、位相が−πからπの２πの範囲において、位相が不
連続であり、負の画素が多い領域に位置するときの動作
を示す。ただし、フィルタは５×５画素の矩形で、−π
からπは−１から１で表示している。３３は正の値をも
つ画素の領域であり、３４は負の値をもつ画素の領域で
ある。また、図６は図５の場合のフィルタリング処理手
段の構成図であり、３５はフィルタリング処理手段前の
干渉画像、３６は第１の手段、３７は第２の手段、 ３
８は第４の手段、３９はノイズ除去画像の画素、４０は
ノイズ除去画像である。
で、位相が−πからπの２πの範囲において、位相が不
連続であり、負の画素が多い領域に位置するときの動作
を示す。ただし、フィルタは５×５画素の矩形で、−π
からπは−１から１で表示している。３３は正の値をも
つ画素の領域であり、３４は負の値をもつ画素の領域で
ある。また、図６は図５の場合のフィルタリング処理手
段の構成図であり、３５はフィルタリング処理手段前の
干渉画像、３６は第１の手段、３７は第２の手段、 ３
８は第４の手段、３９はノイズ除去画像の画素、４０は
ノイズ除去画像である。
【００２９】図５の場合では、第１の手段３６により上
記式（１）を満たすことから、不連続と判定され、さら
に、第２の手段３７により上記式（２）を満たさないこ
とから、負の画素数が多い領域と判定される。
記式（１）を満たすことから、不連続と判定され、さら
に、第２の手段３７により上記式（２）を満たさないこ
とから、負の画素数が多い領域と判定される。
【００３０】第２の手段３７により負の画素数が多いと
判定されると、第４の手段３８により負の画素数が多い
フィルタ領域の平均値を算出する。負の画素数が多いフ
ィルタ領域の平均値をma_avgとすると、 m_avg = （p_avg - 2π + m_avg）／2 （４） によって算出される。上記平均値ma_avgをノイズ除去干
渉画像の画素値３９とする。
判定されると、第４の手段３８により負の画素数が多い
フィルタ領域の平均値を算出する。負の画素数が多いフ
ィルタ領域の平均値をma_avgとすると、 m_avg = （p_avg - 2π + m_avg）／2 （４） によって算出される。上記平均値ma_avgをノイズ除去干
渉画像の画素値３９とする。
【００３１】図７は２４のフィルタ領域が２３の干渉画
像上で、位相が−πからπの２πの範囲において、位相
が不連続ではない領域に位置するときの動作を示す。た
だし、フィルタは５×５画素の矩形で、−πからπは−
１から１で表示している。４１は正の値をもつ画素の領
域であり、４２は負の値をもつ画素の領域である。ま
た、図８は図７の場合のフィルタリング処理手段の構成
図であり、４３はフィルタリング処理手段前の干渉画
像、４４は第１の手段、４５は第５の手段、４６はノイ
ズ除去画像の画素、４７はノイズ除去画像である。
像上で、位相が−πからπの２πの範囲において、位相
が不連続ではない領域に位置するときの動作を示す。た
だし、フィルタは５×５画素の矩形で、−πからπは−
１から１で表示している。４１は正の値をもつ画素の領
域であり、４２は負の値をもつ画素の領域である。ま
た、図８は図７の場合のフィルタリング処理手段の構成
図であり、４３はフィルタリング処理手段前の干渉画
像、４４は第１の手段、４５は第５の手段、４６はノイ
ズ除去画像の画素、４７はノイズ除去画像である。
【００３２】図７の場合では、第１の手段４４により上
記式（１）を満たさないことから、不連続でないと判定
される。第１の手段４４により不連続でないと判定され
ると、第５の手段４５により不連続でない領域の平均値
を算出する。不連続でない領域の平均値をa_avgとする
と、 a_avg = （p_avg + m_avg）／2 （５） によって算出される。上記平均値a_avgをノイズ除去干
渉画像の画素値４６とする。
記式（１）を満たさないことから、不連続でないと判定
される。第１の手段４４により不連続でないと判定され
ると、第５の手段４５により不連続でない領域の平均値
を算出する。不連続でない領域の平均値をa_avgとする
と、 a_avg = （p_avg + m_avg）／2 （５） によって算出される。上記平均値a_avgをノイズ除去干
渉画像の画素値４６とする。
【００３３】以上、フィルタリング処理手段１９により
位相ノイズが除去される動作を示した。上記フィルタの
動作は図２で示したように、干渉画像上を１画素ずつフ
ィルタ領域を移動させながら行われる。上記フィルタの
動作が干渉画像の全ての領域で行われた後、ノイズ除去
干渉画像が得られ、位相アンラップ処理部２０により、
上記ノイズ除去干渉画像を地形に対応するように位相ア
ンラップ処理を行う。位相アンラップ処理は、干渉画像
に現れる等位相線と実際の等高線との対応関係を求める
ための処理であり、いくつかの方法が文献（R. M. Go
ldstein, H. A.Zebker, C. L. Werner, "SATELLITE RAD
AR INTERFEROMETRY:TWO DIMENSIONL PHASE UNWRAPPIN
G", RADIO SCIENCE, VOL.23, 4, PP713-720, JULY 198
8、D. C.Ghiglia, L. A. Robust, "TWO-DIMENSIONAL
WEIGHTED AND UNWEIGHTED PHASE UNWRAPPING THAT USES
FAST TRANSFORMS AND ITERATIVE METHODS", JOURNAL O
F SOCIETY AMERICA.A, 11, 1, pp107-117, JAN. 1994）
等で示されている。
位相ノイズが除去される動作を示した。上記フィルタの
動作は図２で示したように、干渉画像上を１画素ずつフ
ィルタ領域を移動させながら行われる。上記フィルタの
動作が干渉画像の全ての領域で行われた後、ノイズ除去
干渉画像が得られ、位相アンラップ処理部２０により、
上記ノイズ除去干渉画像を地形に対応するように位相ア
ンラップ処理を行う。位相アンラップ処理は、干渉画像
に現れる等位相線と実際の等高線との対応関係を求める
ための処理であり、いくつかの方法が文献（R. M. Go
ldstein, H. A.Zebker, C. L. Werner, "SATELLITE RAD
AR INTERFEROMETRY:TWO DIMENSIONL PHASE UNWRAPPIN
G", RADIO SCIENCE, VOL.23, 4, PP713-720, JULY 198
8、D. C.Ghiglia, L. A. Robust, "TWO-DIMENSIONAL
WEIGHTED AND UNWEIGHTED PHASE UNWRAPPING THAT USES
FAST TRANSFORMS AND ITERATIVE METHODS", JOURNAL O
F SOCIETY AMERICA.A, 11, 1, pp107-117, JAN. 1994）
等で示されている。
【００３４】最終的に、変換処理部２１により、上記ア
ンラップ処理後の位相データを地形の高さに変換し、地
形の３次元情報２２を得る。
ンラップ処理後の位相データを地形の高さに変換し、地
形の３次元情報２２を得る。
【００３５】実施の形態２．実施の形態１では、フィル
タリング処理手段におけるフィルタの動作を示したが、
第６の手段によってフィルタ領域の大きさをユーザが任
意に設定してから、フィルタの動作を開始することがで
きる。図９はこの実施の形態のフィルタリング処理手段
の構成図を示す。４８はフィルタリング処理手段前の干
渉画像、４９はユーザであり、第６の手段５０でフィル
タ領域の大きさが設定され、フィルタ５１の動作により
ノイズ除去干渉画像５２を得る。上記大きさの設定され
たフィルタの動作５１およびフィルタリング処理手段の
後の処理は実施の形態１と同じである。
タリング処理手段におけるフィルタの動作を示したが、
第６の手段によってフィルタ領域の大きさをユーザが任
意に設定してから、フィルタの動作を開始することがで
きる。図９はこの実施の形態のフィルタリング処理手段
の構成図を示す。４８はフィルタリング処理手段前の干
渉画像、４９はユーザであり、第６の手段５０でフィル
タ領域の大きさが設定され、フィルタ５１の動作により
ノイズ除去干渉画像５２を得る。上記大きさの設定され
たフィルタの動作５１およびフィルタリング処理手段の
後の処理は実施の形態１と同じである。
【００３６】実施の形態３．実施の形態２では、フィル
タ領域の大きさをユーザが任意に設定したが、第７の手
段によりフィルタリング処理手段前の干渉画像とフィル
タリング処理手段後のノイズ除去干渉画像を表示し、ユ
ーザがノイズ除去状態を評価してフィルタの大きさを設
定することもできる。図１０はこの実施の形態のフィル
タリング処理手段の構成図を示す。５３はフィルタリン
グ処理手段前の干渉画像であり、５４は前述または後述
のフィルタ領域の大きさを設定する手段であり、５５
a、５５ｂはフィルタの動作である。第７の手段５６
で、フィルタリング処理手段前の干渉画像５７とフィル
タリング処理手段後のノイズ除去干渉画像５８がそれぞ
れモニタやプリンタ等に出力され、ユーザ５９が評価す
る。ユーザが評価したノイズ除去状態からフィルタの大
きさを設定する。最終的に、ノイズ除去干渉画像６０を
得る。上記大きさの設定されたフィルタの動作５５a、
５５ｂおよびフィルタリング処理手段の後の処理は実施
の形態１と同じである。
タ領域の大きさをユーザが任意に設定したが、第７の手
段によりフィルタリング処理手段前の干渉画像とフィル
タリング処理手段後のノイズ除去干渉画像を表示し、ユ
ーザがノイズ除去状態を評価してフィルタの大きさを設
定することもできる。図１０はこの実施の形態のフィル
タリング処理手段の構成図を示す。５３はフィルタリン
グ処理手段前の干渉画像であり、５４は前述または後述
のフィルタ領域の大きさを設定する手段であり、５５
a、５５ｂはフィルタの動作である。第７の手段５６
で、フィルタリング処理手段前の干渉画像５７とフィル
タリング処理手段後のノイズ除去干渉画像５８がそれぞ
れモニタやプリンタ等に出力され、ユーザ５９が評価す
る。ユーザが評価したノイズ除去状態からフィルタの大
きさを設定する。最終的に、ノイズ除去干渉画像６０を
得る。上記大きさの設定されたフィルタの動作５５a、
５５ｂおよびフィルタリング処理手段の後の処理は実施
の形態１と同じである。
【００３７】実施の形態４．実施の形態３では干渉画像
とノイズ除去干渉画像を表示してユーザがフィルタ領域
の大きさを設定したが、第８の手段によりフィルタリン
グ処理手段前の干渉画像の位相ノイズ量とフィルタリン
グ処理手段後のノイズ除去干渉画像の位相ノイズ量を表
示し、ユーザがノイズ除去状態を評価してフィルタの大
きさを設定することもできる。上記位相ノイズ量の指標
となるものはいくつかあるが、ここではResidueを測定
する方法について説明する。
とノイズ除去干渉画像を表示してユーザがフィルタ領域
の大きさを設定したが、第８の手段によりフィルタリン
グ処理手段前の干渉画像の位相ノイズ量とフィルタリン
グ処理手段後のノイズ除去干渉画像の位相ノイズ量を表
示し、ユーザがノイズ除去状態を評価してフィルタの大
きさを設定することもできる。上記位相ノイズ量の指標
となるものはいくつかあるが、ここではResidueを測定
する方法について説明する。
【００３８】Residueは位相ノイズを示す指標であり、J
PLのZebker等の提案によるもので、上記文献に発表さ
れた公知の概念である。図１１は干渉画像上の２×２の
領域（４画素）であり、−πからπを−１から１で示し
てある。もし、位相ノイズがなければ、隣り合う画素の
位相変化は±π以下であるので、上記４画素に対する位
相の周回積分は０になる。しかし、位相ノイズが存在す
ると、周回積分は＋２または−２になる。これをResidu
eという。図１１（a）はResidueが存在しない場合であ
り、図１１（ｂ）は正のResidueが存在する場合であ
る。
PLのZebker等の提案によるもので、上記文献に発表さ
れた公知の概念である。図１１は干渉画像上の２×２の
領域（４画素）であり、−πからπを−１から１で示し
てある。もし、位相ノイズがなければ、隣り合う画素の
位相変化は±π以下であるので、上記４画素に対する位
相の周回積分は０になる。しかし、位相ノイズが存在す
ると、周回積分は＋２または−２になる。これをResidu
eという。図１１（a）はResidueが存在しない場合であ
り、図１１（ｂ）は正のResidueが存在する場合であ
る。
【００３９】図１２はこの実施の形態のフィルタリング
処理手段の構成図を示す。６１はフィルタリング処理手
段前の干渉画像であり、６２は前述または後述のフィル
タ領域の大きさを設定する手段であり、６３a、６３ｂ
はフィルタの動作である。第８の手段６４で、フィルタ
リング処理手段前の干渉画像の位相ノイズ量６５とフィ
ルタリング処理手段後のノイズ除去干渉画像の位相ノイ
ズ量がそれぞれモニタやプリンタ等に出力され、ユーザ
６６が評価する。ユーザが評価したノイズ除去状態から
フィルタの大きさを設定する。最終的に、ノイズ除去干
渉画像６７を得る。上記大きさの設定されたフィルタの
動作６３a、６３ｂおよびフィルタリング処理手段の後
の処理は実施の形態１と同じである。
処理手段の構成図を示す。６１はフィルタリング処理手
段前の干渉画像であり、６２は前述または後述のフィル
タ領域の大きさを設定する手段であり、６３a、６３ｂ
はフィルタの動作である。第８の手段６４で、フィルタ
リング処理手段前の干渉画像の位相ノイズ量６５とフィ
ルタリング処理手段後のノイズ除去干渉画像の位相ノイ
ズ量がそれぞれモニタやプリンタ等に出力され、ユーザ
６６が評価する。ユーザが評価したノイズ除去状態から
フィルタの大きさを設定する。最終的に、ノイズ除去干
渉画像６７を得る。上記大きさの設定されたフィルタの
動作６３a、６３ｂおよびフィルタリング処理手段の後
の処理は実施の形態１と同じである。
【００４０】実施の形態５．実施の形態４では位相ノイ
ズ量を出力してユーザがフィルタ領域の大きさを設定し
たが、第９の手段により過去のフィルタリング処理にお
いて使用したフィルタの大きさと撮像地域を記録してお
き、この記録から同一撮像地域のフィルタリング処理の
ときに効果の高いフィルタの大きさを設定することもで
きる。図１３はこの実施の形態のフィルタリング処理手
段の構成図を示す。６８はフィルタリング処理手段前の
干渉画像、６９はフィルタの動作、７０は過去のフィル
タリング処理に使用したフィルタの大きさと撮像地域の
記録データである。第９の手段７１によって、フィルタ
リング処理のときに、記録データ７０からフィルタの大
きさと撮像地域を参照して、効果の高いフィルタの大き
さを設定することができる。また、効果の高いフィルタ
の大きさと撮像地域を記録データ７０として記録してお
き、次回参照することができる。最終的に、ノイズ除去
干渉画像７２を得る。上記大きさの設定されたフィルタ
の動作６９およびフィルタリング処理手段の後の処理は
実施の形態１と同じである。
ズ量を出力してユーザがフィルタ領域の大きさを設定し
たが、第９の手段により過去のフィルタリング処理にお
いて使用したフィルタの大きさと撮像地域を記録してお
き、この記録から同一撮像地域のフィルタリング処理の
ときに効果の高いフィルタの大きさを設定することもで
きる。図１３はこの実施の形態のフィルタリング処理手
段の構成図を示す。６８はフィルタリング処理手段前の
干渉画像、６９はフィルタの動作、７０は過去のフィル
タリング処理に使用したフィルタの大きさと撮像地域の
記録データである。第９の手段７１によって、フィルタ
リング処理のときに、記録データ７０からフィルタの大
きさと撮像地域を参照して、効果の高いフィルタの大き
さを設定することができる。また、効果の高いフィルタ
の大きさと撮像地域を記録データ７０として記録してお
き、次回参照することができる。最終的に、ノイズ除去
干渉画像７２を得る。上記大きさの設定されたフィルタ
の動作６９およびフィルタリング処理手段の後の処理は
実施の形態１と同じである。
【００４１】実施の形態６．実施の形態５では位相ノイ
ズ量を出力してユーザがフィルタ領域の大きさを設定し
たが、第１０の手段により過去のフィルタリング処理に
おいて使用したフィルタの大きさと干渉画像の位相ノイ
ズ量を記録しておき、この記録から干渉画像の位相ノイ
ズ量に対応したフィルタの大きさを設定することもでき
る。図１４はこの実施の形態のフィルタリング処理手段
の構成図を示す。７３はフィルタリング処理手段前の干
渉画像、７４はフィルタの動作、７５は過去のフィルタ
リング処理に使用したフィルタの大きさと干渉画像の位
相ノイズ量の記録データである。第１０の手段７６によ
って、フィルタリング処理のときに、記録データ７５か
らフィルタの大きさと干渉画像の位相ノイズ量を参照し
て、効果の高いフィルタの大きさを設定することができ
る。また、効果の高いフィルタの大きさと干渉画像の位
相ノイズ量を記録データ７５として記録しておき、次回
参照することができる。最終的に、７７のノイズ除去干
渉画像を得る。上記大きさの設定されたフィルタの動作
７４およびフィルタリング処理手段の後の処理は実施の
形態１と同じである。
ズ量を出力してユーザがフィルタ領域の大きさを設定し
たが、第１０の手段により過去のフィルタリング処理に
おいて使用したフィルタの大きさと干渉画像の位相ノイ
ズ量を記録しておき、この記録から干渉画像の位相ノイ
ズ量に対応したフィルタの大きさを設定することもでき
る。図１４はこの実施の形態のフィルタリング処理手段
の構成図を示す。７３はフィルタリング処理手段前の干
渉画像、７４はフィルタの動作、７５は過去のフィルタ
リング処理に使用したフィルタの大きさと干渉画像の位
相ノイズ量の記録データである。第１０の手段７６によ
って、フィルタリング処理のときに、記録データ７５か
らフィルタの大きさと干渉画像の位相ノイズ量を参照し
て、効果の高いフィルタの大きさを設定することができ
る。また、効果の高いフィルタの大きさと干渉画像の位
相ノイズ量を記録データ７５として記録しておき、次回
参照することができる。最終的に、７７のノイズ除去干
渉画像を得る。上記大きさの設定されたフィルタの動作
７４およびフィルタリング処理手段の後の処理は実施の
形態１と同じである。
【００４２】実施の形態７．実施の形態６では位相ノイ
ズ量を出力してユーザがフィルタ領域の大きさを設定し
たが、第１１の手段により過去のフィルタリング処理に
おいて使用した干渉画像のセンサ名と画像のマルチルッ
ク数とフィルタの大きさを記録しておき、この記録から
に対応したフィルタの大きさを設定することもできる。
図１５はこの実施の形態のフィルタリング処理手段の構
成図を示す。７８はフィルタリング処理手段前の干渉画
像、７９はフィルタの動作、８０は過去のフィルタリン
グ処理に使用した干渉画像のセンサ名と画像のマルチル
ック数とフィルタの大きさの記録データである。第１１
の手段８１によって、フィルタリング処理のときに、記
録データ８０から干渉画像のセンサ名とマルチルック数
とフィルタの大きさを参照して、効果の高いフィルタの
大きさを設定することができる。また、効果の高いフィ
ルタの大きさと使用した干渉画像のセンサ名とマルチル
ック数を記録データ８０として記録しておき、次回参照
することができる。最終的に、ノイズ除去干渉画像８２
を得る。上記大きさの設定されたフィルタの動作７９お
よびフィルタリング処理手段の後の処理は実施の形態１
と同じである。
ズ量を出力してユーザがフィルタ領域の大きさを設定し
たが、第１１の手段により過去のフィルタリング処理に
おいて使用した干渉画像のセンサ名と画像のマルチルッ
ク数とフィルタの大きさを記録しておき、この記録から
に対応したフィルタの大きさを設定することもできる。
図１５はこの実施の形態のフィルタリング処理手段の構
成図を示す。７８はフィルタリング処理手段前の干渉画
像、７９はフィルタの動作、８０は過去のフィルタリン
グ処理に使用した干渉画像のセンサ名と画像のマルチル
ック数とフィルタの大きさの記録データである。第１１
の手段８１によって、フィルタリング処理のときに、記
録データ８０から干渉画像のセンサ名とマルチルック数
とフィルタの大きさを参照して、効果の高いフィルタの
大きさを設定することができる。また、効果の高いフィ
ルタの大きさと使用した干渉画像のセンサ名とマルチル
ック数を記録データ８０として記録しておき、次回参照
することができる。最終的に、ノイズ除去干渉画像８２
を得る。上記大きさの設定されたフィルタの動作７９お
よびフィルタリング処理手段の後の処理は実施の形態１
と同じである。
【００４３】実施の形態８．実施の形態２から７では、
各フィルタの大きさ設定について述べたが、上記の各手
段を選択して、フィルタの動作を再試行できる。図１６
はこの実施の形態のフィルタリング処理手段の構成図を
示す。８３はフィルタリング処理前の干渉画像である。
フィルタの動作８４を行う前に、第６の手段８５ａまた
は第９の手段８６ａまたは第１０の手段８７ａまたは第
１１の手段８８ａのいずれかを選択する第１２の手段８
９によってフィルタの大きさを設定する。上記で設定さ
れた大きさのフィルタで８４のフィルタ動作を行い、ユ
ーザの判断で、最終的にノイズ除去干渉画像として出力
するか、再度フィルタの大きさを設定してフィルタの動
作を再試行するかを決定する。フィルタの動作を再試行
する場合は、第６の手段８５ｂまたは第７の手段９０ま
たは第８の手段９１または第９の手段８６ｂまたは第１
０の手段８７ｂまたは第１１の手段８８ｂのいずれかを
選択する第１３の手段９２によってフィルタの大きさを
設定する。最終的にユーザの判断で決定されたノイズ除
去干渉画像９３が出力される。上記大きさの設定された
フィルタの動作およびフィルタリング処理手段の後の処
理は実施の形態１と同じである。
各フィルタの大きさ設定について述べたが、上記の各手
段を選択して、フィルタの動作を再試行できる。図１６
はこの実施の形態のフィルタリング処理手段の構成図を
示す。８３はフィルタリング処理前の干渉画像である。
フィルタの動作８４を行う前に、第６の手段８５ａまた
は第９の手段８６ａまたは第１０の手段８７ａまたは第
１１の手段８８ａのいずれかを選択する第１２の手段８
９によってフィルタの大きさを設定する。上記で設定さ
れた大きさのフィルタで８４のフィルタ動作を行い、ユ
ーザの判断で、最終的にノイズ除去干渉画像として出力
するか、再度フィルタの大きさを設定してフィルタの動
作を再試行するかを決定する。フィルタの動作を再試行
する場合は、第６の手段８５ｂまたは第７の手段９０ま
たは第８の手段９１または第９の手段８６ｂまたは第１
０の手段８７ｂまたは第１１の手段８８ｂのいずれかを
選択する第１３の手段９２によってフィルタの大きさを
設定する。最終的にユーザの判断で決定されたノイズ除
去干渉画像９３が出力される。上記大きさの設定された
フィルタの動作およびフィルタリング処理手段の後の処
理は実施の形態１と同じである。
【００４４】
【発明の効果】この発明は、以上のように構成されるた
め、以下に示す効果を奏する。
め、以下に示す効果を奏する。
【００４５】干渉型SARによる３次元情報抽出装置は、
２つのSAR画像の相関が十分に得られず、干渉画像に位
相ノイズが発生しても、フィルタリング処理により干渉
画像に存在する位相ノイズを除去することで１組のSAR
画像から高精度な地形の３次元情報を容易に抽出でき
る。
２つのSAR画像の相関が十分に得られず、干渉画像に位
相ノイズが発生しても、フィルタリング処理により干渉
画像に存在する位相ノイズを除去することで１組のSAR
画像から高精度な地形の３次元情報を容易に抽出でき
る。
【００４６】また、上記フィルタリング処理手段は、干
渉画像におけるフィルタ領域の位相が不連続か否かを判
定する第１の手段と、上記第１の手段で不連続と判定さ
れた領域について正の値をもつ画素数が多い領域か負の
画素数が多い領域かを判定する第２の手段と、上記第２
の手段において正の値をもつ画素数が多い領域の平均値
を算出する第３の手段によって、フィルタ領域が干渉画
像上で、位相が不連続であり、正の画素が多い領域にお
いて、等位相線情報を喪失することなく、位相ノイズを
除去することができる。
渉画像におけるフィルタ領域の位相が不連続か否かを判
定する第１の手段と、上記第１の手段で不連続と判定さ
れた領域について正の値をもつ画素数が多い領域か負の
画素数が多い領域かを判定する第２の手段と、上記第２
の手段において正の値をもつ画素数が多い領域の平均値
を算出する第３の手段によって、フィルタ領域が干渉画
像上で、位相が不連続であり、正の画素が多い領域にお
いて、等位相線情報を喪失することなく、位相ノイズを
除去することができる。
【００４７】また、上記フィルタリング処理手段は、上
記第１の手段と、上記第２の手段と、上記第２の手段で
負の値をもつ画素数が多い領域の平均値を算出する第４
の手段によって、フィルタ領域が干渉画像上で、位相が
不連続であり、負の画素が多い領域において、等位相線
情報を喪失することなく、位相ノイズを除去することが
できる。
記第１の手段と、上記第２の手段と、上記第２の手段で
負の値をもつ画素数が多い領域の平均値を算出する第４
の手段によって、フィルタ領域が干渉画像上で、位相が
不連続であり、負の画素が多い領域において、等位相線
情報を喪失することなく、位相ノイズを除去することが
できる。
【００４８】また、上記フィルタリング処理手段は、上
記第１の手段と、上記第１の手段で不連続でないと判定
された領域の平均値を算出する第５の手段によって、フ
ィルタ領域が干渉画像上で、位相が不連続でない領域に
おいて、等位相線情報を喪失することなく、位相ノイズ
を除去することができる。
記第１の手段と、上記第１の手段で不連続でないと判定
された領域の平均値を算出する第５の手段によって、フ
ィルタ領域が干渉画像上で、位相が不連続でない領域に
おいて、等位相線情報を喪失することなく、位相ノイズ
を除去することができる。
【００４９】また、上記フィルタリング処理手段は、フ
ィルタの大きさをユーザが任意に設定する第６の手段に
よって、容易にフィルタの大きさを変えることができ、
効果的に位相ノイズを除去することができる。
ィルタの大きさをユーザが任意に設定する第６の手段に
よって、容易にフィルタの大きさを変えることができ、
効果的に位相ノイズを除去することができる。
【００５０】また、上記フィルタリング処理手段は、干
渉画像とノイズ除去干渉画像を表示し、ユーザがノイズ
除去状態を評価してフィルタの大きさを設定する第７の
手段によって、効果的に位相ノイズを除去することがで
きる。
渉画像とノイズ除去干渉画像を表示し、ユーザがノイズ
除去状態を評価してフィルタの大きさを設定する第７の
手段によって、効果的に位相ノイズを除去することがで
きる。
【００５１】また、上記フィルタリング処理手段は、干
渉画像のノイズ量およびノイズ除去干渉画像のノイズ量
を表示し、ユーザがノイズ除去状態を評価してフィルタ
の大きさを調整する第８の手段によって、効果的に位相
ノイズを除去することができ、さらに、ノイズ量の指標
を上記Residueにして、Residueの数を削減することは、
位相アンラップ処理（文献）を高速化および単純化
し、高精度な地形の３次元情報を容易に抽出できる。
渉画像のノイズ量およびノイズ除去干渉画像のノイズ量
を表示し、ユーザがノイズ除去状態を評価してフィルタ
の大きさを調整する第８の手段によって、効果的に位相
ノイズを除去することができ、さらに、ノイズ量の指標
を上記Residueにして、Residueの数を削減することは、
位相アンラップ処理（文献）を高速化および単純化
し、高精度な地形の３次元情報を容易に抽出できる。
【００５２】また、上記フィルタリング処理手段は、過
去のフィルタリング処理において使用したフィルタの大
きさと撮像地域を記録しておき、この記録から同一撮像
地域のフィルタリング処理のときに効果の高いフィルタ
の大きさを設定する第９の手段によって、撮像地域に効
果的なフィルタの大きさを容易に設定することができ、
効果的に位相ノイズを除去することができる。
去のフィルタリング処理において使用したフィルタの大
きさと撮像地域を記録しておき、この記録から同一撮像
地域のフィルタリング処理のときに効果の高いフィルタ
の大きさを設定する第９の手段によって、撮像地域に効
果的なフィルタの大きさを容易に設定することができ、
効果的に位相ノイズを除去することができる。
【００５３】また、上記フィルタリング処理手段は、過
去のフィルタリング処理において使用したフィルタの大
きさと干渉画像の位相ノイズ量を記録しておき、この記
録から干渉画像の位相ノイズ量に対応したフィルタの大
きさを設定する第 １０の手段によって、効果的に位相
ノイズを除去することができ、さらに、ノイズ量の指標
を上記Residueにして、Residueの数を削減することは、
位相アンラップ処理を高速化および単純化し、高精度な
地形の３次元情報を容易に抽出できる。
去のフィルタリング処理において使用したフィルタの大
きさと干渉画像の位相ノイズ量を記録しておき、この記
録から干渉画像の位相ノイズ量に対応したフィルタの大
きさを設定する第 １０の手段によって、効果的に位相
ノイズを除去することができ、さらに、ノイズ量の指標
を上記Residueにして、Residueの数を削減することは、
位相アンラップ処理を高速化および単純化し、高精度な
地形の３次元情報を容易に抽出できる。
【００５４】また、上記フィルタリング処理手段は、過
去のフィルタリング処理において使用した干渉画像のセ
ンサ名と画像のマルチルック数とフィルタの大きさを記
録しておき、この記録からに対応したフィルタの大きさ
を設定する第１１の手段によって、センサの観測条件か
ら干渉画像に効果的なフィルタの大きさを容易に設定す
ることができ、効果的に位相ノイズを除去することがで
きる。
去のフィルタリング処理において使用した干渉画像のセ
ンサ名と画像のマルチルック数とフィルタの大きさを記
録しておき、この記録からに対応したフィルタの大きさ
を設定する第１１の手段によって、センサの観測条件か
ら干渉画像に効果的なフィルタの大きさを容易に設定す
ることができ、効果的に位相ノイズを除去することがで
きる。
【００５５】さらに、上記各フィルタの大きさを設定す
る手段を適宜選択してフィルタの動作を再試行すること
により、効果的に位相ノイズを除去することができる。
る手段を適宜選択してフィルタの動作を再試行すること
により、効果的に位相ノイズを除去することができる。
【図１】 この発明の干渉型SARによる３次元情報抽出
装置を示す図である。
装置を示す図である。
【図２】 この発明の実施の形態１で示すフィルタリン
グ処理手段により位相ノイズが除去される様子である。
グ処理手段により位相ノイズが除去される様子である。
【図３】 この発明の実施の形態１で示すフィルタリン
グ処理手段で用いるフィルタの動作を示す図である。
グ処理手段で用いるフィルタの動作を示す図である。
【図４】 この発明の実施の形態１で示すフィルタリン
グ処理手段の構成図である。
グ処理手段の構成図である。
【図５】 この発明の実施の形態１で示すフィルタリン
グ処理手段で用いるフィルタの動作を示す図である。
グ処理手段で用いるフィルタの動作を示す図である。
【図６】 この発明の実施の形態１で示すフィルタリン
グ処理手段の構成図である。
グ処理手段の構成図である。
【図７】 この発明の実施の形態１で示すフィルタリン
グ処理手段で用いるフィルタの動作を示す図である。
グ処理手段で用いるフィルタの動作を示す図である。
【図８】 この発明の実施の形態１で示すフィルタリン
グ処理手段の構成図である。
グ処理手段の構成図である。
【図９】 この発明の実施の形態２で示すフィルタリン
グ処理手段の構成図である。
グ処理手段の構成図である。
【図１０】 この発明の実施の形態３で示すフィルタリ
ング処理手段の構成図である。
ング処理手段の構成図である。
【図１１】 この発明の実施の形態３で示すResidueの
存在例である。
存在例である。
【図１２】 この発明の実施の形態４で示すフィルタリ
ング処理手段の構成図である。
ング処理手段の構成図である。
【図１３】 この発明の実施の形態５で示すフィルタリ
ング処理手段の構成図である。
ング処理手段の構成図である。
【図１４】 この発明の実施の形態６で示すフィルタリ
ング処理手段の構成図である。
ング処理手段の構成図である。
【図１５】 この発明の実施の形態７で示すフィルタリ
ング処理手段の構成図である。
ング処理手段の構成図である。
【図１６】 この発明の実施の形態８で示すフィルタリ
ング処理手段の構成図である。
ング処理手段の構成図である。
【図１７】 飛翔体を人工衛星としたときの干渉型SAR
を用いて地形の３次元情報を抽出するシステムの概念図
である。
を用いて地形の３次元情報を抽出するシステムの概念図
である。
【図１８】 従来の干渉型SARによる３次元情報抽出装
置を示す図である。
置を示す図である。
【図１９】 従来の干渉型SARによる３次元情報抽出装
置を示す図である。
置を示す図である。
１a SAR搭載人工衛星、１ｂ SAR搭載人工衛星、３ 送
受信部、４ 画像処理部、５ FFT部、６ フィルタ
部、７ IFFT部、８ 干渉処理部、９ 差分抽出部、１
０ 送受信部、１１ 画像処理部、１２ 干渉処理部、
１３ マルチルック処理部、１４ 位相アンラップ処理
部、１５ 変換処理部、１６a 送受信部、１６ｂ 送受
信部、１７ 画像処理部、１８ 干渉処理部、１９ フ
ィルタリング処理手段、２０ 位相アンラップ処理部、
２１ 変換処理部、８ 第１の手段、２９ 第２の手
段、３０ 第３の手段、３６ 第１の手段、３７ 第２
の手段、３８ 第４の手段、４４ 第１の手段、４５
第５の手段、４９ ユーザ、５０ 第６の手段、５４
フィルタの大きさを設定する手段、５６ 第７の手段、
５９ ユーザ、６２ フィルタの大きさを設定する手
段、６４ 第８の手段、６６ ユーザ、７１ 第９の手
段、７６ 第１０の手段、８１ 第１１の手段、８５a
第６の手段、８５ｂ 第６の手段、８６a 第９の手
段、８６ｂ 第９の手段、８７a 第１０の手段、８７
ｂ 第１０の手段、８８a 第１１の手段、８８ｂ 第
１１の手段、８９ 第１２の手段、９０ 第７の手段、
９１ 第８の手段、９２ 第１３の手段。
受信部、４ 画像処理部、５ FFT部、６ フィルタ
部、７ IFFT部、８ 干渉処理部、９ 差分抽出部、１
０ 送受信部、１１ 画像処理部、１２ 干渉処理部、
１３ マルチルック処理部、１４ 位相アンラップ処理
部、１５ 変換処理部、１６a 送受信部、１６ｂ 送受
信部、１７ 画像処理部、１８ 干渉処理部、１９ フ
ィルタリング処理手段、２０ 位相アンラップ処理部、
２１ 変換処理部、８ 第１の手段、２９ 第２の手
段、３０ 第３の手段、３６ 第１の手段、３７ 第２
の手段、３８ 第４の手段、４４ 第１の手段、４５
第５の手段、４９ ユーザ、５０ 第６の手段、５４
フィルタの大きさを設定する手段、５６ 第７の手段、
５９ ユーザ、６２ フィルタの大きさを設定する手
段、６４ 第８の手段、６６ ユーザ、７１ 第９の手
段、７６ 第１０の手段、８１ 第１１の手段、８５a
第６の手段、８５ｂ 第６の手段、８６a 第９の手
段、８６ｂ 第９の手段、８７a 第１０の手段、８７
ｂ 第１０の手段、８８a 第１１の手段、８８ｂ 第
１１の手段、８９ 第１２の手段、９０ 第７の手段、
９１ 第８の手段、９２ 第１３の手段。
フロントページの続き (72)発明者 高橋 史子 東京都千代田区丸の内二丁目２番３号 三 菱電機株式会社内 (72)発明者 水野 政冶 東京都千代田区丸の内二丁目２番３号 三 菱電機株式会社内 (72)発明者 和泉 秀幸 東京都千代田区丸の内二丁目２番３号 三 菱電機株式会社内 Ｆターム(参考） 5B057 BA23 CE02 CE06 DA07 DB03 DC09 5J070 AB01 AC02 AD01 AD05 AE07 AF06 AF08 AH02 AH14 AH19 AH31 AH34 AH50 AJ05 AJ13 AK22 AK28 AK40 AL02 BE02 BE04
Claims (11)
- 【請求項１】 合成開口レーダ（Synthetic Aperture R
adar; SAR）で得られた観測条件の異なる２つの観測デ
ータから、SAR画像再生処理によって得られた２つの複
素SAR画像データを互いに干渉させて干渉画像を得る手
段と、上記干渉画像を地形に対応するように位相アンラ
ップ処理し、上記アンラップ処理後の位相データを地形
の高さに変換し３次元情報を得る手段とを備えた干渉型
合成開口レーダによる地形の３次元情報抽出装置におい
て、上記干渉画像に存在する位相ノイズをフィルタリン
グ処理により除去し、ノイズ除去干渉画像を得るフィル
タリング処理手段を有することを特徴とする干渉型合成
開口レーダによる地形の３次元情報抽出装置。 - 【請求項２】 上記フィルタリング処理手段は、干渉画
像におけるフィルタ領域の位相が不連続か否かを判定す
る第１の手段と、上記第１の手段で不連続と判定された
領域について正の値をもつ画素数が多い領域か負の画素
数が多い領域かを判定する第２の手段と、上記第２の手
段において正の値をもつ画素数が多い領域の平均値を算
出する第３の手段とを有することを特徴とする請求項１
記載の干渉型合成開口レーダによる地形の３次元情報抽
出装置。 - 【請求項３】 上記フィルタリング処理手段は、干渉画
像におけるフィルタ領域の位相が不連続か否かを判定す
る第１の手段と、上記第１の手段で不連続と判定された
領域について正の値をもつ画素数が多い領域か負の画素
数が多い領域かを判定する第２の手段と、上記第２の手
段で負の値をもつ画素数が多い領域の平均値を算出する
第４の手段とを有することを特徴とする請求項１記載の
干渉型合成開口レーダによる地形の３次元情報抽出装
置。 - 【請求項４】 上記フィルタリング処理手段は、干渉画
像におけるフィルタ領域の位相が不連続か否かを判定す
る第１の手段と、上記第１の手段で不連続でないと判定
された領域の平均値を算出する第５の手段とを有するこ
とを特徴とする請求項１記載の干渉型合成開口レーダに
よる地形の３次元情報抽出装置。 - 【請求項５】 上記フィルタリング処理手段は、フィル
タの大きさを任意に設定する第６の手段を有することを
特徴とする請求項２〜４いずれか記載の干渉型合成開口
レーダによる地形の３次元情報抽出装置。 - 【請求項６】 上記フィルタリング処理手段は、干渉画
像とノイズ除去干渉画像を表示し、ユーザがノイズ除去
状態を評価してフィルタの大きさを設定する第７の手段
を有することを特徴とする請求項２〜４いずれか記載の
干渉型合成開口レーダによる地形の３次元情報抽出装
置。 - 【請求項７】 上記フィルタリング処理手段は、干渉画
像のノイズ量およびノイズ除去干渉画像のノイズ量を表
示し、ユーザがノイズ除去状態を評価してフィルタの大
きさを調整する第８の手段を有することを特徴とする請
求項２〜４いずれか記載の干渉型合成開口レーダによる
地形の３次元情報抽出装置。 - 【請求項８】 上記フィルタリング処理手段は、過去の
フィルタリング処理において使用したフィルタの大きさ
と撮像地域を記録しておき、この記録から同一撮像地域
のフィルタリング処理のときに効果の高いフィルタの大
きさを設定する第９の手段を有することを特徴とする請
求項２〜４いずれか記載の干渉型合成開口レーダによる
地形の３次元情報抽出装置置。 - 【請求項９】 上記フィルタリング処理手段は、過去の
フィルタリング処理において使用したフィルタの大きさ
と干渉画像の位相ノイズ量を記録しておき、この記録か
ら干渉画像の位相ノイズ量に対応したフィルタの大きさ
を設定する第１０の手段を有することを特徴とする請求
項２〜４いずれか記載の干渉型合成開口レーダによる地
形の３次元情報抽出装置。 - 【請求項１０】 上記フィルタリング処理手段は、過去
のフィルタリング処理において使用した干渉画像のセン
サ名と画像のマルチルック数とフィルタの大きさを記録
しておき、この記録からに対応したフィルタの大きさを
設定する第１１の手段を有することを特徴とする請求項
２〜４いずれか記載の干渉型合成開口レーダによる地形
の３次元情報抽出装置。 - 【請求項１１】 上記フィルタリング処理手段は、フィ
ルタの大きさを任意に設定する第６の手段と、干渉画像
のノイズ量およびノイズ除去干渉画像のノイズ量を表示
し、ユーザがノイズ除去状態を評価してフィルタの大き
さを調整する第８の手段と、過去のフィルタリング処理
において使用したフィルタの大きさと撮像地域を記録し
ておき、この記録から同一撮像地域のフィルタリング処
理のときに効果の高いフィルタの大きさを設定する第９
の手段と、過去のフィルタリング処理において使用した
フィルタの大きさと干渉画像の位相ノイズ量を記録して
おき、この記録から干渉画像の位相ノイズ量に対応した
フィルタの大きさを設定する第 １０の手段と、過去の
フィルタリング処理において使用した干渉画像のセンサ
名と画像のマルチルック数とフィルタの大きさを記録し
ておき、この記録からに対応したフィルタの大きさを設
定する第１１の手段と、干渉画像に対する位相ノイズ除
去の１回目の試行において、上記第６の手段、第９の手
段、第１０の手段または第１１の手段のいずれかを選択
する第１２の手段と、上記干渉画像に対する位相ノイズ
除去の２回目以降の複数回の再試行において、第６の手
段、第７の手段、第８の手段、第９の手段、第１０の手
段または第１１の手段のいずれかを再試行の度に選択す
る第１３の手段とを有することを特徴とする請求項２〜
４いずれか記載の干渉型合成開口レーダによる地形の３
次元情報抽出装置。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP25870399A JP2001083243A (ja) | 1999-09-13 | 1999-09-13 | 干渉型合成開口レーダによる地形の３次元情報抽出装置 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP25870399A JP2001083243A (ja) | 1999-09-13 | 1999-09-13 | 干渉型合成開口レーダによる地形の３次元情報抽出装置 |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| JP2001083243A true JP2001083243A (ja) | 2001-03-30 |
Family
ID=17323938
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP25870399A Pending JP2001083243A (ja) | 1999-09-13 | 1999-09-13 | 干渉型合成開口レーダによる地形の３次元情報抽出装置 |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JP2001083243A (ja) |
Cited By (16)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2003070763A (ja) * | 2001-08-28 | 2003-03-11 | Ge Medical Systems Global Technology Co Llc | 位相矛盾検出方法および装置、位相矛盾解消方法および装置、並びに、磁気共鳴撮影装置 |
| KR100425283B1 (ko) * | 2003-04-18 | 2004-03-31 | (주)충청측량설계공사 | 간섭측정용 합성 개구 레이다의 기하학적 특성을 이용한지형고도 측정방법 |
| JP2004191053A (ja) * | 2002-12-06 | 2004-07-08 | Mitsubishi Electric Corp | 合成開口レーダ装置及び数値標高モデル作成方法 |
| CN101908206A (zh) * | 2010-07-01 | 2010-12-08 | 西北工业大学 | 一种基于mca的sar图像抑噪方法 |
| CN102628942A (zh) * | 2012-04-24 | 2012-08-08 | 中国科学院遥感应用研究所 | 一种雷达影像双视向信息补偿方法 |
| CN102914773A (zh) * | 2012-09-11 | 2013-02-06 | 电子科技大学 | 一种多航过圆周sar三维成像方法 |
| CN103148842A (zh) * | 2013-02-04 | 2013-06-12 | 国家海洋局第二海洋研究所 | 一种基于遥感图像特征的浅海沙波区多波束测深地形重构方法 |
| JP2013195102A (ja) * | 2012-03-16 | 2013-09-30 | Mitsubishi Space Software Kk | レーダ画像処理装置、レーダ画像処理プログラムおよびレーダ画像処理方法 |
| JP2015184854A (ja) * | 2014-03-24 | 2015-10-22 | 株式会社ジオ技術研究所 | 画像処理装置 |
| CN105954748A (zh) * | 2016-04-26 | 2016-09-21 | 中国科学院电子学研究所 | 基于分块策略的TS-InSAR大气相位滤波方法 |
| CN111854596A (zh) * | 2020-08-26 | 2020-10-30 | 中国地质大学（北京） | 基于地基干涉雷达数据的地表形变实时监测方法及系统 |
| WO2020240720A1 (ja) * | 2019-05-29 | 2020-12-03 | 日本電気株式会社 | 合成開口レーダの信号処理装置および信号処理方法 |
| JPWO2021024336A1 (ja) * | 2019-08-05 | 2021-02-11 | ||
| WO2022145111A1 (ja) * | 2020-12-28 | 2022-07-07 | 日本電気株式会社 | 情報処理装置、情報処理方法、及びコンピュータ読み取り可能な記憶媒体 |
| CN116148853A (zh) * | 2023-02-20 | 2023-05-23 | 中国矿业大学 | 一种矿区沉陷干涉相位滤波方法和装置 |
| JP7505597B2 (ja) | 2020-12-28 | 2024-06-25 | 日本電気株式会社 | 情報処理装置、情報処理方法、及びプログラム |
-
1999
- 1999-09-13 JP JP25870399A patent/JP2001083243A/ja active Pending
Cited By (24)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2003070763A (ja) * | 2001-08-28 | 2003-03-11 | Ge Medical Systems Global Technology Co Llc | 位相矛盾検出方法および装置、位相矛盾解消方法および装置、並びに、磁気共鳴撮影装置 |
| JP2004191053A (ja) * | 2002-12-06 | 2004-07-08 | Mitsubishi Electric Corp | 合成開口レーダ装置及び数値標高モデル作成方法 |
| KR100425283B1 (ko) * | 2003-04-18 | 2004-03-31 | (주)충청측량설계공사 | 간섭측정용 합성 개구 레이다의 기하학적 특성을 이용한지형고도 측정방법 |
| CN101908206A (zh) * | 2010-07-01 | 2010-12-08 | 西北工业大学 | 一种基于mca的sar图像抑噪方法 |
| JP2013195102A (ja) * | 2012-03-16 | 2013-09-30 | Mitsubishi Space Software Kk | レーダ画像処理装置、レーダ画像処理プログラムおよびレーダ画像処理方法 |
| CN102628942A (zh) * | 2012-04-24 | 2012-08-08 | 中国科学院遥感应用研究所 | 一种雷达影像双视向信息补偿方法 |
| CN102628942B (zh) * | 2012-04-24 | 2013-10-23 | 中国科学院遥感与数字地球研究所 | 一种雷达影像双视向信息补偿方法 |
| CN102914773A (zh) * | 2012-09-11 | 2013-02-06 | 电子科技大学 | 一种多航过圆周sar三维成像方法 |
| CN103148842A (zh) * | 2013-02-04 | 2013-06-12 | 国家海洋局第二海洋研究所 | 一种基于遥感图像特征的浅海沙波区多波束测深地形重构方法 |
| JP2015184854A (ja) * | 2014-03-24 | 2015-10-22 | 株式会社ジオ技術研究所 | 画像処理装置 |
| CN105954748A (zh) * | 2016-04-26 | 2016-09-21 | 中国科学院电子学研究所 | 基于分块策略的TS-InSAR大气相位滤波方法 |
| WO2020240720A1 (ja) * | 2019-05-29 | 2020-12-03 | 日本電気株式会社 | 合成開口レーダの信号処理装置および信号処理方法 |
| JPWO2020240720A1 (ja) * | 2019-05-29 | 2020-12-03 | ||
| US20220221577A1 (en) * | 2019-05-29 | 2022-07-14 | Nec Corporation | Synthetic aperture radar signal processing device and signal processing method |
| JP7248109B2 (ja) | 2019-05-29 | 2023-03-29 | 日本電気株式会社 | 合成開口レーダの信号処理装置および信号処理方法 |
| US11841422B2 (en) | 2019-05-29 | 2023-12-12 | Nec Corporation | Synthetic aperture radar signal processing device and signal processing method |
| JPWO2021024336A1 (ja) * | 2019-08-05 | 2021-02-11 | ||
| WO2021024336A1 (ja) * | 2019-08-05 | 2021-02-11 | 日本電気株式会社 | 位相アンラップ装置及び位相アンラップ方法 |
| JP7255690B2 (ja) | 2019-08-05 | 2023-04-11 | 日本電気株式会社 | 位相アンラップ装置及び位相アンラップ方法 |
| CN111854596A (zh) * | 2020-08-26 | 2020-10-30 | 中国地质大学（北京） | 基于地基干涉雷达数据的地表形变实时监测方法及系统 |
| WO2022145111A1 (ja) * | 2020-12-28 | 2022-07-07 | 日本電気株式会社 | 情報処理装置、情報処理方法、及びコンピュータ読み取り可能な記憶媒体 |
| JP7505597B2 (ja) | 2020-12-28 | 2024-06-25 | 日本電気株式会社 | 情報処理装置、情報処理方法、及びプログラム |
| CN116148853A (zh) * | 2023-02-20 | 2023-05-23 | 中国矿业大学 | 一种矿区沉陷干涉相位滤波方法和装置 |
| CN116148853B (zh) * | 2023-02-20 | 2024-02-20 | 中国矿业大学 | 一种矿区沉陷干涉相位滤波方法和装置 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| JP2001083243A (ja) | 干渉型合成開口レーダによる地形の３次元情報抽出装置 | |
| AU2019262089B2 (en) | Method and system for use in colourisation of a point cloud | |
| US5043734A (en) | Discrete autofocus for ultra-high resolution synthetic aperture radar | |
| JP7077961B2 (ja) | 画像解析装置、画像解析方法及びプログラム | |
| US20130230206A1 (en) | Foliage penetration based on 4d lidar datasets | |
| TW201111750A (en) | Method and system of processing data gathered using a range sensor | |
| CA2579898A1 (en) | Method for the processing and representing of ground images obtained by synthetic aperture radar systems (sar) | |
| EP3629055B1 (en) | Method and apparatus for acquiring depth image, and electronic device | |
| JPWO2019202628A1 (ja) | 路面検出装置、路面検出装置を利用した画像表示装置、路面検出装置を利用した障害物検知装置、路面検出方法、路面検出方法を利用した画像表示方法、および路面検出方法を利用した障害物検知方法 | |
| JP2590689B2 (ja) | 干渉型合成開口レーダ装置および地形変動観測方法 | |
| CN114763997A (zh) | 车辆采集到的雷达点云数据的处理方法、装置及电子设备 | |
| JP3695103B2 (ja) | レーダ信号処理装置 | |
| EP2544020B1 (fr) | Procédé et dispositif de détection d'une cible masquée par des réflecteurs de forte énergie | |
| CN110907132B (zh) | 一种波向探测方法、系统、设备、介质 | |
| Mancini et al. | Automating nearshore bathymetry extraction from wave motion in satellite optical imagery | |
| US11175399B2 (en) | Information processing device, information processing method, and storage medium | |
| KR102670901B1 (ko) | 해상 부유체 탐지 장치 및 해상 부유체 탐지 방법 | |
| CN113409457A (zh) | 立体图像的三维重构与可视化方法及设备 | |
| JP3292368B2 (ja) | テクスチャマッピング方法及び装置ならびに記録媒体 | |
| EP1697764B1 (fr) | Procede de visualisation des echos recus par un sonar actif utilisant une emission a spectre de raies | |
| JP2020060496A (ja) | 情報処理装置 | |
| Leberl et al. | Radargrammetric image processing with a softcopy stereo workstation | |
| US20230401678A1 (en) | Information processing device, information processing method, and information processing system | |
| Seemann et al. | Local analysis of inhomogeneous sea surfaces in coastal waters using nautical radar image sequences | |
| Gelautz et al. | Automated matching experiments with different kinds of SAR imagery |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| RD01 | Notification of change of attorney |
Free format text: JAPANESE INTERMEDIATE CODE: A7421 Effective date: 20040628 |
|
| A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20050822 |
|
| A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20070424 |
|
| A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20070619 |
|
| A02 | Decision of refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A02 Effective date: 20071016 |