[go: up one dir, main page]

JP7384295B2 - 画像解析装置および画像解析方法 - Google Patents

画像解析装置および画像解析方法 Download PDF

Info

Publication number
JP7384295B2
JP7384295B2 JP2022553250A JP2022553250A JP7384295B2 JP 7384295 B2 JP7384295 B2 JP 7384295B2 JP 2022553250 A JP2022553250 A JP 2022553250A JP 2022553250 A JP2022553250 A JP 2022553250A JP 7384295 B2 JP7384295 B2 JP 7384295B2
Authority
JP
Japan
Prior art keywords
parameter
complex matrix
evaluation
candidates
matrix
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.)
Active
Application number
JP2022553250A
Other languages
English (en)
Other versions
JPWO2022070243A5 (ja
JPWO2022070243A1 (ja
Inventor
大地 田中
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.)
NEC Corp
Original Assignee
NEC 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 NEC Corp filed Critical NEC Corp
Publication of JPWO2022070243A1 publication Critical patent/JPWO2022070243A1/ja
Publication of JPWO2022070243A5 publication Critical patent/JPWO2022070243A5/ja
Application granted granted Critical
Publication of JP7384295B2 publication Critical patent/JP7384295B2/ja
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

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
    • G01S13/9021—SAR image post-processing 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/9023—SAR image post-processing techniques combined with interferometric 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
    • G06—COMPUTING OR CALCULATING; COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00—Image analysis
    • G06T7/97—Determining parameters from multiple pictures
    • G—PHYSICS
    • G06—COMPUTING OR CALCULATING; COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00—Indexing scheme for image analysis or image enhancement
    • G06T2207/10—Image acquisition modality
    • G06T2207/10028—Range image; Depth image; 3D point clouds

Landscapes

  • Engineering & Computer Science (AREA)
  • Remote Sensing (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Electromagnetism (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Theoretical Computer Science (AREA)
  • Artificial Intelligence (AREA)
  • Radar Systems Or Details Thereof (AREA)

Description

本発明は、合成開口レーダの受信電磁波から生成される干渉画像に基づいて標高や変位を解析するための画像解析装置および画像解析方法に関する。
合成開口レーダ（ＳＡＲ：Synthetic Aperture Radar）技術は、人工衛星や航空機などの飛翔体が移動しながら電磁波を送受信し、大きな開口を持つアンテナと等価な観測画像が得られる技術である。合成開口レーダは、例えば、地表からの反射波を信号処理して、標高や地表変位を解析するために利用される。ＳＡＲ技術が用いられる場合、解析装置は、合成開口レーダによって得られる時系列のＳＡＲ画像（ＳＡＲデータ）を入力とし、入力されたＳＡＲ画像を時系列解析する。
標高や地表変位を解析するための有効な手法として、干渉ＳＡＲ解析がある。干渉ＳＡＲ解析では、違う時期に撮影された複数（例えば、２枚）のＳＡＲ画像を構成する電波信号の位相差が計算される。そして、撮影時期間で生じた飛翔体と地面間の距離の変化が検出される。
特許文献１には、コヒーレンス行列を用いる解析手法が記載されている。コヒーレンスは、Ｎ（Ｎ≧２）枚のＳＡＲ画像のうちの複数のＳＡＲ画像における同じ位置にあたる画素の複素相関で計算される。ＳＡＲ画像のペアを（ｍ，ｎ）とし、コヒーレンス行列の成分をｃｍ，ｎとする。ｍ，ｎは、それぞれ、Ｎ以下の値であり、Ｎ枚のＳＡＲ画像のいずれかを示す。ＳＡＲ画像のペアについて、位相θｍ，ｎ（具体的には、位相差）が算出される。そして、コヒーレンス算出対象の画素を含む所定領域内の複数の画素についてｅｘｐ（－ｊθｍ，ｎ）が平均化された値が、コヒーレンス行列の成分ｃｍ，ｎとなる。
ｃｍ，ｎの偏角∠ｃｍ，ｎは、平均的な位相（具体的には、位相差）に相当する。ｃｍ ，ｎの絶対値｜ｃｍ，ｎ｜から、位相θｍ，ｎの分散の大きさを把握可能である。
コヒーレンス行列は、ノイズが除去された場合の位相を推測可能な情報を含む。また、位相ノイズ量の程度を推測可能な情報（すなわち、分散）を含む。
地表等の変位解析のために、位相θｍ，ｎが変位速度および撮影時刻差に相関することが利用される。例えば、位相差の平均値に基づいて変位が推定される。なお、位相ノイズ量を利用して変位解析の精度を検証することが可能である。したがって、コヒーレンス行列は、変位解析のために使用可能である。
標高解析のために、位相θｍ，ｎが解析対象物の標高および飛翔体間距離（例えば、飛翔体の２つの撮影位置の間の距離）に相関することが利用される。例えば、位相差の平均値に基づいて標高が推定される。なお、位相ノイズ量を利用して標高解析の精度を検証することが可能である。したがって、コヒーレンス行列は、標高解析のために使用可能である。
国際公開第２０１１／００３８３６号
特許文献１に記載された手法によれば、位相差に含まれるノイズを減らすことができる。また、特許文献１に記載された手法を利用して、地表や物体の変位（変位量）および標高を推定することができる。例えば、下記の（１）式で表される評価式を最大化する、標高に関連する値ｈおよび変位速度に関連する値ｖが算出される。以下、ｈを標高といい、ｖを変位速度という。
Figure 0007384295000001
（１）式において、Ｗは、重み行列を示す。Ｗとして、例えば、コヒーレンス行列の絶対値のべき乗が用いられる。「Ｈ」は随伴行列であることを表す。（１）式におけるｓは、例えば（２）式のように表現される。λは、使用電磁波の波長を示す。Ｂperpは、Ｎ枚の画像のそれぞれが取得されたときの軌道間距離を要素とする行ベクトルの要素を転置したベクトルを示す。Ｂtempは、Ｎ枚の画像のそれぞれの撮影時刻を要素とする行ベクトルの要素を転置したベクトルを示す。ＢperpおよびＢtempは、以下のように表される。
Figure 0007384295000002
Ｂperp＝（b1 b2 ・・・ bN）T
Ｂtemp＝（t1 t2 ・・・ tN）T
Σに関して、例えばコヒーレンス行列として、ｅｘｐ（－ｊθｍ，ｎ）やＡｍ・Ａｎ・ｅｘｐ（－ｊθｍ，ｎ）を平均したものが用いられる。ＡｍおよびＡｎについては後述される。
図１４は、標高および変位の推定結果（解析結果）の一例を示す説明図である。なお、変位は、（変位速度×時間）である。図１４において、縦軸は、評価値を示す。評価式が最大化されることによって、最適値（最適な標高および変位）が推定されるので、最適値は、評価式を最大化する値、すなわち、ノイズの影響が低下している、評価値が最も大きい値である。
ところが、観測時の電磁波の反射環境などに起因してノイズの発生状況は様々に変化する。その結果、例えば、図１４に例示されている２番目に評価値が大きい第２のピークが生ずることがある。実は、第２のピークに対応する標高および変位が真の最適値であったり、真の最適値が求められていないという可能性がある。
本発明は、解析結果の信頼性を提示することが可能になる画像解析装置および画像解析方法を提供することを目的とする。
本発明による画像解析装置は、同一地域が記録された複数の画像における全ての画像の組での位相差が反映された複素行列を算出する複素行列算出手段と、位相の変化を説明するパラメータの候補を複数選定するパラメータ候補選定手段と、複素行列と所定の重み行列とを用いて、複数の候補の尤もらしさを評価する候補評価手段と、パラメータの候補を尤もらしさによって重み付けして、パラメータの候補の統計値を算出する統計値算出手段とを含む。
本発明による画像解析方法は、同一地域が記録された複数の画像における全ての画像の組での位相差が反映された複素行列を算出し、位相の変化を説明するパラメータの候補を複数選定し、複素行列と所定の重み行列とを用いて、複数の候補の尤もらしさを評価し、パラメータの候補を尤もらしさによって重み付けして、パラメータの候補の統計値を算出する。
本発明による画像解析プログラムは、コンピュータに、同一地域が記録された複数の画像における全ての画像の組での位相差が反映された複素行列を算出する処理と、位相の変化を説明するパラメータの候補を複数選定する処理と、複素行列と所定の重み行列とを用いて、複数の候補の尤もらしさを評価する処理と、パラメータの候補を尤もらしさによって重み付けして、パラメータの候補の統計値を算出する処理とを実行させる。
本発明によれば、解析結果の信頼性を提示することが可能になる。
第１の実施形態の画像解析装置を示すブロック図である。 第１の実施形態の画像解析装置の動作を示すフローチャートである。 第１の実施形態の画像解析装置の作用を説明するための説明図である。 第２の実施形態の画像解析装置の構成例を示すブロック図である。 第２の実施形態の画像解析装置の動作を示すフローチャートである。 第３の実施形態の画像解析装置における候補評価部の構成例を示すブロック図である。 第３の実施形態の画像解析装置における候補評価部の動作を示すフローチャートである。 第３の実施形態の効果を説明するための説明図である。 第３の実施形態の効果を説明するための説明図である。 第４の実施形態の画像解析装置の構成例を示すブロック図である。 低周波フィルタ部の構成例を示すブロック図である。 第４の実施形態の画像解析装置の動作を示すフローチャートである。 ＣＰＵを有するコンピュータの一例を示すブロック図である。 画像解析装置の主要部を示すブロック図である。 標高および変位の解析結果の一例を示す説明図である。
以下、本発明の実施形態を図面を参照して説明する。
実施形態１．
図１は、第１の実施形態の画像解析装置の構成例を示すブロック図である。図１に示す画像解析装置１０は、複素行列算出部１１０、パラメータ候補選定部１２０、候補評価部１３０、統計値算出部１４０、および乗算部１５０を含む。
ＳＡＲ画像格納部２００には、あらかじめＮ枚のＳＡＲ画像が格納される。重み行列格納部３００には、重み行列があらかじめ格納される。本実施形態では、重み行列格納部３００には、どの画像のペアを重視して位相のフィッティングを行うかを決定する重み行列が保存されている。なお、ＳＡＲ画像格納部２００および重み行列格納部３００は、画像解析装置１０の外部に設置されていてもよいし、画像解析装置１０に含まれていてもよい。
重み行列は、例えば、重み行列算出部（図１において図示せず）によって算出される。重み行列算出部は、ノイズが多い画像（ＳＡＲ画像）ペアｍ，ｎ（ｍ、ｎ：１～Ｎ）間では小さく、ノイズが少ない画像ペアｍ，ｎ間では大きくなるような値を、ｍ行ｎ列の重み行列の要素として計算する。そのような重み行列として、コヒーレンス行列の絶対値の逆行列に負号を付けたものを使用可能である。また、重み行列算出部は、一例として、複素行列算出部１１０の出力を、加算された画素の数で除算した値をコヒーレンス行列の推定値として使用可能である。また、重み行列算出部は、コヒーレンス行列の絶対値そのもののべき乗を重み行列として使用してもよい。なお、重み行列算出部は、画像解析装置１０に含まれていてもよい。
複素行列算出部１１０は、画像番号ｍの画像と画像番号ｎの画像について、画像番号ｎの画素の値の複素共役を計算する。すなわち、複素行列算出部１１０は、同一地域が記録された複数の画像における全ての画像の組での位相差が反映された複素行列を算出する。また、複素行列算出部１１０は、画像番号ｍの画像の画素の各々と画像番号ｎの画素の値の共役複素数と乗算する。さらに、複素行列算出部１１０は、注目画素から所定範囲内にある画素について、乗算の結果を加算することによって、ｍ行ｎ列の要素に持つような複素行列を作成する。複素行列算出部１１０は、行列を作成する処理を、全ての画素の各々を対象として、すなわち全ての画素の各々を注目画素として実行する。なお、複素行列算出部１１０は、注目画素から所定範囲内にある画素について、乗算の結果を加算する場合に、注目画素の値を乗算結果としてもよい。その場合には、複素行列算出部１１０は、実質的に加算処理を行わない。換言すれば、複素行列算出部１１０は、所定範囲の画素を用いた平均化処理を実行しない。
パラメータ候補選定部１２０は、位相変化に関連するパラメータすなわち位相差を説明するパラメータの候補の値を選定する。位相差を説明するパラメータは、例えば、変位速度と標高である。ただし、それらは一例であって、他の類いのパラメータが用いられたり、変位速度および標高とともに他の類いのパラメータが用いられてもよい。パラメータの候補の値は、一の種類のパラメータ（例えば、変位速度または標高）における選定される値（例えば、標高について、３ｍ、４ｍ、５ｍなど）である。
なお、変位速度を例にすると、（撮影時刻差×変位速度）が位相差に相当する。標高を例にすると、（軌道間距離×標高）が位相差に相当する。
候補評価部１３０は、パラメータ候補選定部１２０が選定した全てのパラメータの候補値の尤もらしさを評価する。
統計値算出部１４０は、尤もらしさで重み付けしてパラメータの候補の平均値や分散を計算する。
乗算部１５０は、複素行列算出部１１０の出力である複素行列と重み行列格納部３００に格納されている重み行列とを乗算し、乗算結果を候補評価部１３０に出力する。
次に、図２のフローチャートを参照して本実施形態の画像解析装置１０の動作を説明する。
複素行列算出部１１０は、複素行列を生成する処理を行う（ステップＳ１０１）。すなわち、複素行列算出部１１０は、画像番号ｍの画像と画像番号ｎの画像について、画素の各々を対象として（注目画素として）、画像ｍにおける画素と、画像ｎの共役複素数との積を、所定範囲で加算する。所定範囲は、例えば、注目画素の近傍の１つまたは複数の画素である。その場合には、複素行列算出部１１０は、注目画素の近傍の１つまたは複数の画素に関する積を加算する。また、所定範囲は、あらかじめ定められている所定の基準に従った範囲であってもよい。複素行列算出部１１０は、その場合には、注目画素に対しての基準を満たす１つまたは複数の画素に関する積を加算する。さらに、複素行列算出部１１０は、所定範囲で加算された値をｍ行ｎ列の要素に持つ複素行列を作成する。そして、複素行列算出部１１０は、作成した複素行列を出力する。
複素行列算出部１１０は、ステップＳ１０１の処理を実行することによって、平均化された位相すなわちノイズの影響が低減された位相の推定値を含む行列であって、ノイズの特性が把握可能な行列を算出することになる。なお、複素行列算出部１１０は、ステップＳ１０１の処理で、コヒーレンス行列を算出してもよい。
コヒーレンス行列を算出する場合、複素行列算出部１１０は、例えばＳＡＲ画像格納部２００に格納されているＮ枚のＳＡＲ画像（複素画像：振幅および位相情報を含む。）を対象としてコヒーレンス行列Ｃを算出する。例えば、ＳＡＲ画像のペアを（ｍ，ｎ）とし、コヒーレンス行列Ｃの成分をｃm，nとする。ｍ，ｎは、それぞれ、Ｎ以下の値であり、Ｎ枚のＳＡＲ画像のいずれかを示す。複素行列算出部１１０は、ＳＡＲ画像のペアについて、位相θm，n（具体的には、位相差）を算出する。そして、複素行列算出部１１０は、コヒーレンス算出対象の画素を含む所定領域内の複数の画素についてｅｘｐ（－ｊθm，n）を平均化した値を、コヒーレンス行列Ｃの成分ｃm，nとする。
また、ＳＡＲ画像m における強度をＡｍ、ＳＡＲ画像n における強度をＡｎとして、Ａｍ・Ａｎ・ｅｘｐ（－ｊθｍ，ｎ）を平均化してもよい。また、Ａｍ・Ａｎ・ｅｘｐ（－ｊθｍ，ｎ）の平均として得られた行列に対して、その対角成分の総和をＮで割った値によって、行列の要素それぞれを割ってもよい。また、Ａｍ・Ａｎ・ｅｘｐ（－ｊθｍ，ｎ）の平均として得られた行列に対して、その対角成分の-1/2乗を対角成分として持つ対角行列を左右から乗算してもよい。
なお、このようなコヒーレンス行列の算出方法は一例であり、コヒーレンス行列算出部は、種々の手法でコヒーレンス行列を算出できる。
パラメータ候補選定部１２０は、位相差を説明するパラメータの候補の値を選定する（ステップＳ１０２）。上述したように、パラメータは、変位速度や標高である。例えば、パラメータが変位速度である場合には、パラメータ候補選定部１２０は、変位速度の複数の値を候補の値とする。パラメータが標高である場合には、パラメータ候補選定部１２０は、標高の複数の値を候補の値とする。なお、パラメータ候補選定部１２０は、候補の値として、画像解析装置１０の使用状況に応じて想定しうる値を、あらかじめ決められている範囲内から一定間隔で選定する。また、パラメータ候補選定部１２０は、候補の値を、あらかじめ決められている範囲内からランダムに選定してもよい。
候補評価部１３０は、パラメータ候補選定部１２０が選定した全ての候補の値を、複素行列算出部１１０が作成した複素行列を用いて評価する（ステップＳ１０３）。候補評価部１３０は、ステップＳ１０３の処理で、（複素行列×重み行列）を含む式を評価式として、パラメータの全ての候補の値を評価する。すなわち、候補評価部１３０は、全ての候補の値の尤もらしさを評価する。なお、（複素行列×重み行列）は、乗算部１５０によって算出される。
候補評価部１３０が尤もらしさを評価することは、候補の値と（複素行列×重み行列）を含む評価式の値との合致度を求めることを意味する。なお、候補評価部１３０が、（複素行列×重み行列）を含む式を評価式として用いるということは、例えば上記の（１）式を用いることに相当する。候補評価部１３０は、（１）式を用いる場合、（１）式におけるΣとして、複素行列算出部１１０が作成した複素行列を用いる。
統計値算出部１４０は、候補の値の統計値を算出する（ステップＳ１０４）。統計値算出部１４０は、ステップＳ１０４の処理で、合致度（尤もらしさ）で重み付けしたパラメータ候補の値の平均値（重み付け平均値）と、合致度で重み付けしたパラメータ候補の値の分散（重み付け分散）とを算出する。
図３は、画像解析装置の作用を説明するための説明図である。図３には、標高または変位の解析結果の一例に対して、統計値（分散および平均値）が記入されている。図３において、分散（具体的には、分散の平方根）が矢印に対応し、平均値が中央の黒丸に対応する。
本実施形態では、画像解析装置１０において、候補評価部１３０が算出した合致度を用いて統計値算出部１４０が統計値を算出するので、統計値によって、パラメータ（例えば、標高や変位）の推定値がどの程度信頼できるのかを把握可能になる。
なお、本実施形態において、１種類のパラメータに関する評価を実行できるが、２種類のパラメータに関する評価を一括して実行することもできる。
実施形態２．
図４は、第２の実施形態の画像解析装置の構成例を示すブロック図である。図４に示す画像解析装置２０は、第１の実施形態の画像解析装置１０の構成要素に加えて、パラメータ複素行列算出部１６０を含む。なお、本実施形態における統計値算出部１４１は、第１の実施形態における統計値算出部１４０の機能に加えて、後述する機能も有する。
位相差を説明する２種類以上のパラメータが存在する場合に、パラメータ候補選定部１２０が選定するパラメータの候補の数が膨大になることもある。そのような場合には、複素行列算出部１１０や候補評価部１３０などの演算量が多大になる。例えば、変位速度に関する候補の数が１０００であり、標高に関する候補の数が１０００である場合には、１０００×１０００についての演算が求められる。本実施形態では、画像解析装置２０は、他の種類のパラメータに関する評価処理に、ある種類のパラメータに関する推定結果を用いることによって、演算量を低減する。
パラメータ複素行列算出部１６０は、例えば、下記の（３）式を用いて複素行列Ｌk,m, nを算出する。以下、パラメータ複素行列算出部１６０が算出した複素行列を、パラメー
タ複素行列という。パラメータ複素行列は、パラメータの候補のそれぞれについて、各画像における位相の値を持つような複素行列である。
Ｌk,m,n＝ｅｘｐ（－ｊｐk（Ｂm－Ｂn）） （３）
（３）式において、ｍ，ｎは、画像番号を示す。ｋは、パラメータの候補番号を示す。よって、"k,m,n"は、画像番号ｍのＳＡＲ画像と画像番号ｎのＳＡＲ画像との間でのパラ
メータの候補を示すインデックスである。ｐkは、パラメータの候補の値を示す。Ｂは、
説明変数に相当する。パラメータが変位速度である場合には、Ｂm，Ｂnは、例えば撮影時刻を示す。その場合、ｐkは、変位速度である。パラメータが標高である場合には、Ｂm，Ｂnは、例えば軌道間距離を示す。その場合、ｐkは、標高である。
なお、複素行列Ｌk,m,nは、位相差を説明するパラメータの候補で説明される位相差を
打ち消すような行列に相当する。
また、本実施形態では、統計値算出部１４１は、候補評価部１３０が算出したパラメータの候補ｐkの尤もらしさ（合致度）をｓkとした場合、ｓkの平均値（具体的には、重み
付け平均値）の他に、例えば（４）式で表されるパラメータ複素行列の平均値Ａも算出する。ただし、（４）式において、Σkはkに対する総和を計算する演算子である。
Ａ＝（ΣkＬk,m,nｓk）／（Σkｓk） （４）
また、本実施形態では、乗算部１５０は、複素行列算出部１１０の出力である複素行列と重み行列格納部３００に格納されている重み行列と統計値算出部１４１の出力であるパラメータ複素行列の平均値Ａとを乗算し、乗算結果を候補評価部１３０に出力する。なお、乗算部１５０は、統計値算出部１４１は、一の種類のパラメータに関する評価が行われているときに、他の種類のパラメータのパラメータ複素行列を乗算する。よって、乗算部１５０の出力において、他の種類のパラメータに関する位相差が除外されている。
次に、図５のフローチャートを参照して本実施形態の画像解析装置２０の動作を説明する。なお、複素行列算出部１１０および候補評価部１３０の動作は、第１の実施形態における動作と同じである。
図５に例示される処理では、一の種類のパラメータに関する処理と他の種類のパラメータに関する処理とが交互に実行される。以下、一の種類のパラメータを変位速度ｖとし、他の種類のパラメータを標高ｈとする。そして、まず、標高ｈに関する処理が実行される場合を例にする。
標高ｈに関する処理が実行されるときに、パラメータ複素行列算出部１６０は、変位速度ｖのパラメータ複素行列を算出する（ステップＳ２０１）。統計値算出部１４１は、変位速度ｖに関する尤もらしさｓkの重み付き平均値を算出する（ステップＳ２０２）。重
み付き平均値は、変位速度ｖが尤もらしいときの位相が反映されたパラメータ複素行列に相当する。また、統計値算出部１４１は、変位速度ｖに関するパラメータ複素行列の平均値Ａも算出する。
乗算部１５０は、複素行列算出部１１０の出力である複素行列と重み行列格納部３００に格納されている重み行列と統計値算出部１４１の出力であるパラメータ複素行列の平均値Ａとを乗算する（ステップＳ２０３）。
候補評価部１３０は、第１の実施形態の場合と同様にして、標高ｈの全ての候補の値の尤もらしさを評価する（ステップＳ２０４）。上述したように、パラメータ複素行列は、位相差を説明するパラメータの候補（この場合には、変位速度ｖの候補）で説明される位相差を打ち消すような行列に相当するので、候補評価部１３０は、変位速度ｖの影響が排除された状態で標高ｈを対象とした評価を実施することができる。
終了条件が成立した場合には、ステップＳ２０６に移行する。終了条件が成立していない場合には、ステップＳ２０１に戻る。なお、ステップＳ２０１～Ｓ２０４の処理をループ処理とする。ステップＳ２０１に戻って次のループ処理が実行される場合に、処理対象が切り替えられる。この例では、処理対象が、標高ｈから変位速度ｖに切り替えられる。
変位速度ｖに関する処理が実行されるときに、パラメータ複素行列算出部１６０は、標高ｈのパラメータ複素行列を算出する（ステップＳ２０１）。統計値算出部１４１は、標高ｈに関する尤もらしさｓkの重み付き平均値を算出する（ステップＳ２０２）。重み付
き平均値は、標高ｈが尤もらしいときの位相が反映されたパラメータ複素行列に相当する。また、統計値算出部１４１は、標高ｈに関するパラメータ複素行列の平均値Ａも算出する。
乗算部１５０は、複素行列算出部１１０の出力である複素行列と重み行列格納部３００に格納されている重み行列と統計値算出部１４１の出力であるパラメータ複素行列の平均値Ａとを乗算する（ステップＳ２０３）。
候補評価部１３０は、変位速度ｖの全ての候補の値の尤もらしさを評価する（ステップＳ２０４）。上述したように、パラメータ複素行列は、位相差を説明するパラメータの候補（この場合には、標高ｈの候補）で説明される位相差を打ち消すような行列に相当するので、候補評価部１３０は、標高ｈの影響が排除された状態で変位速度ｖの評価を実施することができる。
終了条件が成立した場合には、ステップＳ２０６に移行する。終了条件が成立していない場合には、ステップＳ２０１に戻る。
終了条件は、例えば、全てのパラメータ（この例では、変位速度ｖおよび標高ｈ）の尤もらしさが、前回のステップＳ２０４の処理で得られた値に対して、あらかじめ定められた所定値未満であることである。すなわち、終了条件は、パラメータの尤もらしさの変動量が小さくなっていることである。換言すれば、終了条件は、パラメータの尤もらしさが最適値に収束したと判定されたことである。また、終了条件として、例えば、ステップＳ２０１～Ｓ２０４のループ処理が所定回実行されるという条件が用いられてもよい。
本実施形態では、一のパラメータに関する評価処理が他のパラメータに関する評価処理とは独立して実行されるので、複数のパラメータに関する評価処理が一括して実行される場合に比べて、処理が効率化される。その結果、画像解析装置２０における演算量が低減する。
また、一つのパラメータに関する評価処理を行っているときに、他のパラメータの候補で説明される位相差を打ち消すようなパラメータ複素行列を使用して、他のパラメータに依存する位相差の影響を受けないようにすることができる。その結果、統計値算出部１４１が算出する平均値および分散が最適値に収束するまでの時間が短くなる。
なお、本実施形態では、２つの種類のパラメータが例示されたが、画像解析装置２０は、３つ以上の種類のパラメータに関して評価処理を実行することができる。また、画像解析装置２０は、複数種類のパラメータに関する評価処理を、並行して同時に実行するように構成されていてもよい。
実施形態３．
図６は、第３の実施形態の画像解析装置における候補評価部の構成例を示すブロック図である。第３の実施形態では、第１の実施形態および第２の実施形態における候補評価部１３０に代えて、候補評価部１３１が設けられる。
図６に示す候補評価部１３１は、変換部１３２、ランダム位相生成部１３３、評価部１３４、評価部１３５、統合部１３６、および比較部１３７を含む。
変換部１３２は、パラメータの候補の値（例えば、変位速度ｖの候補の値、標高ｈの候補の値）に説明変数Ｂを乗算して、パラメータの候補の値を位相に相当する値（以下、位相という。）に変換する。ランダム位相生成部１３３は、全ＳＡＲ画像の各々に対応して、ランダムな位相（ランダム位相）の組を１つ以上生成する。
評価部１３４は、変換部１３２が出力した位相を評価して評価値（第１の評価値）を出力する。評価部１３５は、ランダム位相評価部１３３が出力したランダム位相の組を評価して評価値（第２の評価値）を出力する。統合部１３６は、ランダム位相の組の評価値を統合して統合評価値とする。比較部１３７は、変換部１３２が出力した各位相（パラメータの候補に対応する各位相）の評価値と統合評価値とを比較する。
候補評価部１３１以外の画像解析装置の構成は、第１の実施形態の画像解析装置１０または第２の実施形態の画像解析装置２０の構成と同じである。
次に、図７のフローチャートを参照して本実施形態の画像解析装置における候補評価部１３１および統計値算出部１４１の動作を説明する。なお、複素行列算出部１１０および候補評価部１３０の動作は、第１の実施形態における動作と同じである。
変換部１３２は、パラメータの候補の値（例えば、変位速度ｖの候補の値、標高ｈの候補の値）に説明変数Ｂを乗算して位相に変換する（ステップＳ３０１）。一例として、変換部１３２は、ＢｖやＢｈを算出する。なお、変換部１３２は、Ｂｖ＋Ｂｈの演算を行ってもよい。上述したように、説明変数Ｂは、一例として、パラメータが変位速度である場合には、Ｂm，Ｂn（例えば、撮影時刻）である。他の例として、説明変数Ｂは、パラメータが標高である場合には、Ｂm，Ｂn（例えば、軌道間距離）である。変換部１３２は、位相に相当する値（以下、位相という。）を評価部１３４に出力する
ランダム位相生成部１３３は、全ＳＡＲ画像の各々に対応するランダムな位相（ランダム位相）の組を１つ以上生成する（ステップＳ３０２）。例えば、ＳＡＲ画像の枚数がＮ枚であるとすると、１つの位相の組には、Ｎ個の位相が含まれる。ランダム位相生成部１３３は、生成したランダム位相の組を評価部１３５に出力する。
なお、ランダム位相生成部１３３は、パラメータの候補を選択する範囲としてあらかじめ決められている範囲内で、ランダム位相を生成することが好ましい。
評価部１３４は、変換部１３２が出力した位相を評価し、評価部１３５は、ランダム位相生成部１３３が出力したランダム位相の組を評価する（ステップＳ３０３）。
統合部１３６は、ランダム位相の組の評価値を統合して統合評価値とする（ステップＳ３０４）。統合の方法として、例えば、各評価値の平均値を統合評価値とする方法がある。統合部１３６は、各評価値のうちの最頻値、各評価値のうちの中心値、各評価値のうちの最大値、評価値の分位点などを統合評価値としてもよい。
比較部１３７は、変換部１３２が出力した各位相（パラメータの候補に対応する各位相）の評価値と統合評価値とを比較する。そして、比較部１３７は、統合評価値を超えているパラメータの候補に対応する位相の評価値については、パラメータの候補に対応する位相の評価値を出力する。比較部１３７は、統合評価値以下のパラメータの候補に対応する位相の評価値については、その評価値を統合評価値に置き換える（ステップＳ３０５）。
図８Ａおよび図８Ｂは、本実施形態の効果を説明するための説明図である。なお、図８Ａおよび図８Ｂにおいて、「パラメータの候補の範囲」は、パラメータの候補を選択する範囲としてあらかじめ決められている範囲である。
図８Ａに示すように、ステップＳ３０５の処理が存在しない場合には、実は最適値ではない評価値のピークＡに対応するパラメータの値が尤も確からしいと判定されるおそれがある。
しかし、本実施形態では、図８Ｂに示すように、パラメータの候補に対応する位相の評価値をそのまま出力するのか統合評価値を出力するのかを決めるためのレベルＢ（例えば、統合評価値に相当）があるので、妥当でないパラメータの値が尤も確からしいと判定される可能性が低減する。
実施形態４．
図９は、第４の実施形態の画像解析装置における候補評価部の構成例を示すブロック図である。図９に示す画像解析装置４０は、第２の実施形態の画像解析装置２０の構成要素に加えて、低周波フィルタ部１７０を含む。
大気中の水分に起因して発生する空間的に低周波の位相成分がＳＡＲ画像に含まれる場合に、変位の推定結果が不正確になるおそれがある。そこで、本実施形態では、低周波フィルタ部１７０を用いて、空間的に低周波の位相成分を除去する。なお、本実施形態では、低周波の位相成分が、滑らかに分布する位相差を形成することが利用される。
図１０は、低周波フィルタ部１７０の構成例を示すブロック図である。低周波フィルタ部１７０は、画素選択部１７１、低次元圧縮部１７２、展開部１７３、および最適化部１７４を含む。
画素選択部１７１は、同一地域が記録され、同じ画素位置が同じ地点を表すような位置合わせが施された複数のＳＡＲ画像から、ある１枚のＳＡＲ画像中の画素を複数選択する。そして、画素選択部１７１は、他のＳＡＲ画像における、上記の１枚のＳＡＲ画像中の画素と同一の位置の画素を特定する。画素選択部１７１は、特定した画素の位置を指定する。低次元圧縮部１７２は、画素選択部１７１が指定した位置の画素を低次元空間に写像する。展開部１７３は、画素選択部１７１が指定した位置の画素を低次元空間から元の画素空間に写像する。
最適化部１７４は、上記の１枚のＳＡＲ画像と他のＳＡＲ画像とのペアについてのコヒーレンス行列を重みとして用いる。なお、画素選択部１７１は、１枚のＳＡＲ画像と他のＳＡＲ画像とを指定する。そして、最適化部１７４は、展開部１７３から出力された演算結果と画素選択部１７１が指定した位置の観測画素の画素値（観測画素値）との両方に対して近くなるような位相を持つ画素（複素ベクトルで表される。後述するｘsamp,sに相当する。）を導出する。
なお、最適化部１７４は、空間的に滑らかに分布する位相差を評価可能な評価関数（以下、空間相関評価関数という。）を最適化（例えば、最大化）することによって、位相を持つ複素ベクトルを展開部１７３から出力された演算結果に近づける。また、最適化部１７４は、観測信号評価関数を最適化（例えば、最大化）することによって、複素ベクトルを画素選択部１７１が指定した位置の観測画素の画素値に近づける。また、観測信号評価関数が最適化されるときに、複素ベクトルは、ノイズが比較的小さい領域における観測画素の画素値に近づく。
なお、本実施形態では、演算量を低減するために、低周波フィルタ部１７０は、全画素について評価処理を行うのではなく、全画素からランダムに選択された画素を対象として評価処理を行う。また、低周波フィルタ部１７０は、後述する空間相関評価関数における空間相関行列Ｋを低次元分解（低ランク近似）して最適化処理を実行する。
なお、詳しくは、
画像解析装置は、空間相関行列Ｋの逆行列Ｋ-1を実質的に低次元分解する。画像解析装置は、例えば、Nystrom 近似と呼ばれる手法を用いて低次元分解を行うことができる。その場合、画像解析装置は、（５）式を用いて空間相関行列Ｋの逆行列Ｋ-1を低次元分解する。
Figure 0007384295000003
（５）式において、Ｇは、ｄ×ｄの行列である（ｄ＜Ｎ）。Ｖは、Ｎ×ｄの行列（Ｖ＝ｖ１，ｖ２，・・・，ｖｄ）である。Λは、Ｎ×Ｎの対角行列である。
次に、図１１のフローチャートを参照して本実施形態の画像解析装置４０の動作を説明する。なお、（５）式で表される空間相関行列Ｋの逆行列Ｋ-1が画像解析装置１０にあらかじめ設定される。なお、低次元圧縮部１７２は、空間相関行列Ｋの逆行列Ｋ-1の一部である行列ＧＶTを使用する。また、展開部１７３は、空間相関行列Ｋの逆行列Ｋ-1の一部
である行列Ｖを使用する。よって、行列ＧＶTと行列Ｖとが、画像解析装置４０に設定さ
れてもよい。
第２の実施形態と同様、本実施形態でも、一の種類のパラメータに関する処理と他の種類のパラメータに関する処理とが交互に実行される。以下、一の種類のパラメータを変位速度ｖとし、他の種類のパラメータを標高ｈとする。そして、まず、標高ｈに関する処理が実行される場合を例にする。すなわち、画像解析装置４０は、まず、標高ｈについて処理を実行する。
画素選択部１７１は、Ｍ個（Ｍ＜Ｎ）の乱数を生成する（ステップＳ４０１）。Ｎは、ＳＡＲ画像の総画素数を示す。なお、Ｍは２以上の値である。
なお、本実施形態では、画素選択部１７１は、ランダムに画素を選択するが、ランダムに選択しなくてもよい。例えば、あらかじめ定められている所定のルールに従って、所定の領域内のＮ個の画素からＭ個の画素を選択してもよい。
以下、選択される画素の画素番号をｎｍとする（ｍ：１～Ｍ）。なお、ｎｍは、１～Ｍのうちのいずれかになる。Ｓは、総ＳＡＲ画像数を示す。画像番号をｓとする（ｓ：１～Ｓ）。画素選択部１７１は、生成した複数の乱数のそれぞれを、選択される画素の画素番号とする。すなわち、画素選択部１７１は、ランダムに画素を選択する。
展開部１７３は、前回のｆs にＶnmを乗算することによってＶmeans,sを算出する（ス
テップＳ４０２）。前回のｆs は、直前に実行されたステップＳ４１０の処理の演算結果を意味する。Ｖnmのそれぞれは、行列Ｖから、ｎ１行、ｎ２行、・・・、ｎＭ行を取り出して生成される行ベクトルである。すなわち、展開部１７３は、画素番号ｎ１，ｎ２，・・・，ｎＭ行のそれぞれの画素について、行列Ｖにおける対応する行ベクトルを用いてχmeans,sを算出する。以下、χmeans,sを、大気位相推定値ということがある。
最適化部１７４は、全パラメータ（この例では、標高ｈとへｎに速度ｖ）に関して、画素番号画素番号ｎ１，ｎ２，・・・，ｎＭの画素のパラメータ複素行列の平均値を乗算して、Ｍm（ｍ：１～Ｍ）とする（ステップＳ４０３）。
最適化部１７４は、（６）式を用いて、一時的な値であるｘtempを得る（ステップＳ４０４）。（６）式において、ｓおよびｓ’は、Ｓ枚のＳＡＲ画像から選択される画像の画像番号に相当する。
Figure 0007384295000004
最適化部１７４は、ｘtempを活性化してｘnm,sを得る（ステップＳ４０５）。なお、（６）式では、ｘnmにおける"n"が下付文字で表されているが、ｘnmは、（６）式における"n"が下付文字で表されている変数と同じである。
最適化部１７４は、ステップＳ４０５の処理で、下記の（７）式のように、非線形関数ｇ（ａ）を用いて、ｘtempを非線形変換する。
Figure 0007384295000005
非線形関数ｇ（ａ）として、以下のような関数を使用可能である。
ｇ（ａ）＝Ｉ1（２ａ）／Ｉ0（２ａ） （８）
Figure 0007384295000006
ｇ（ａ）＝ｔａｎｈ（ａ） （１０）
（８）式におけるＩ０およびＩ１は、それぞれ、０次と１次の第一種変形ベッセル関数である。（８）式を用いる場合、ｇ（ａ）の値は、ａの値が大きくなるほど１に近づく。（９）式を用いる場合にも（１０）式を用いる場合にも、ｇ（ａ）の値は、ａの値が大きくなるほど１に近づく。つまり、本実施形態では、ｘnm,sの最大値（絶対値で）は、１に制限される。
なお、（８）～（１０）式に例示した関数に限られず、最適化部１７４は、出力値が、ａの増加に伴って０から特定の正値に漸近するような他の関数を使用可能である。
なお、ステップＳ４０５の処理で得られるｘnmは、大気位相の不確かさが含まれる推定値に相当する。
本実施形態では、乗算部１５０は、複素行列算出部１１０から出力される画素ｎmにお
ける複素行列と、重み行列格納部３００に格納されている重みＷと、他のパラメータ（標高ｈが扱われているパラメータである場合には、変位速度ｖ）に関するパラメータ複素行列の平均値と、（１１）式で表される値をｓ行ｓ’列の要素として持つ複素行列とを乗算する（ステップＳ４０６）。乗算部１５０による乗算処理で（１１）式で表される値をｓ行ｓ’列の要素として持つ複素行列が乗算されることによって、大気中の水分に依存する位相成分を抑圧することが可能になる。
Figure 0007384295000007
最適化部１７４は、第１の実施形態および第２の実施形態の場合と同様にして、扱っているパラメータ（例えば、標高ｈ）の全ての候補の値の尤もらしさを評価する（ステップＳ４０７）。上述したように、パラメータ複素行列は、位相差を説明するパラメータの候補（この場合には、変位速度ｖの候補）で説明される位相差を打ち消すような行列に相当するので、候補評価部１３０は、変位速度ｖの影響が排除された状態で標高ｈを対象とした評価を実施することができる。
終了条件が成立した場合には、ステップＳ４１０に移行する。終了条件が成立していない場合には、ステップＳ４０３に戻る。なお、ステップＳ４０３～Ｓ４０８の処理を第１ループ処理とする。すなわち、ステップＳ４０９では、第１ループ処理の終了条件が成立したか否か判定される。第１ループ処理の終了条件は、一例として、第１ループ処理があらかじめ決められている所定回実行されたことである。
ステップＳ４１０では、低次元圧縮部１７２は、最適化部１７４によって最適化されたｘnm,sを、ＧＶTnmを用いて圧縮する（ステップＳ４１０）。すなわち、低次元圧縮部１
７２は、ｘsamp,sとＧＶTnmとを乗算する。低次元圧縮部１７２は、さらに、乗算結果を
（Ｎ／Ｍ）倍する。乗算結果が（Ｎ／Ｍ）倍された値を、ｆs と表現する。なお、Ｎ個の全画素から選択されたＭ個の画素を対象として第１ループ処理が実行されるので、ｘsamp ,sとＧＶTnmとの乗算結果は、本来の値の（Ｍ／Ｎ）倍になっている。低次元圧縮部１７
２は、乗算結果を（Ｎ／Ｍ）倍して本来の値に戻す。
終了条件が成立した場合には、処理を終了する。終了条件が成立していない場合には、ステップＳ４０１に戻る。なお、ステップＳ４０１～Ｓ４１０４の処理を第２ループ処理とする。ステップＳ４０１に戻って次の第２ループ処理が実行される場合に、処理対象（扱うパラメータ）が切り替えられる。扱われていたパラメータが標高ｈであった場合には、処理対象が標高ｈから変位速度ｖに切り替えられる。扱われていたパラメータが変位速度ｖであった場合には、処理対象が変位速度ｖから標高ｈに切り替えられる。
終了条件は、例えば、全てのパラメータ（この例では、変位速度ｖおよび標高ｈ）の尤もらしさが、前回のステップＳ４０７の処理で得られた値に対して、あらかじめ定められた所定値未満であることである。すなわち、終了条件は、パラメータの尤もらしさの変動量が小さくなっていることである。換言すれば、終了条件は、パラメータの尤もらしさが最適値に収束したと判定されたことである。また、終了条件として、例えば、第２ループ処理が所定回実行されるという条件が用いられてもよい。
本実施形態では、低周波フィルタ部１７０によって、大気位相の不確かさが含まれる演算が実行されるので、大気中の水分に起因して発生する位相成分を除去することが可能になる。
図１２は、ＣＰＵ（Central Processing Unit）を有するコンピュータの一例を示すブ
ロック図である。コンピュータは、上記の各実施形態の画像解析装置に実装される。ＣＰＵ１０００は、記憶装置１００１に格納されたプログラム（ソフトウェア要素：コード）に従って処理を実行することによって、上記の実施形態における各機能を実現する。すなわち、図１，図４，図６，図９，図１０に示された画像解析装置１０，２０，４０における、複素行列算出部１１０、パラメータ候補選定部１２０、候補評価部１３０，１３１、統計値算出部１４０，１４１、乗算部１５０、パラメータ複素行列算出部１６０、低周波フィルタ部１７０、変換部１３２、ランダム位相生成部１３３、評価部１３４、評価部１３５、統合部１３６、比較部１３７、画素選択部１７１、低次元圧縮部１７２、展開部１７３、および最適化部１７４の機能を実現する。なお、ＣＰＵ１０００に代えて、ＧＰＵ（Graphics Processing Unit）、または、ＣＰＵとＧＰＵとの組み合わせを使用することもできる。
記憶装置１００１は、例えば、非一時的なコンピュータ可読媒体（non-transitory computer readable medium ）である。非一時的なコンピュータ可読媒体は、様々なタイプの実体のある記録媒体（tangible storage medium ）を含む。非一時的なコンピュータ可読媒体の具体例として、磁気記録媒体（例えば、ハードディスク）、光磁気記録媒体（例えば、光磁気ディスク）、ＣＤ－ＲＯＭ（Compact Disc-Read Only Memory ）、ＣＤ－Ｒ（Compact Disc-Recordable ）、ＣＤ－Ｒ／Ｗ（Compact Disc-ReWritable ）、半導体メモリ（例えば、マスクＲＯＭ、ＰＲＯＭ（Programmable ROM）、ＥＰＲＯＭ（Erasable PROM ）、フラッシュＲＯＭ）がある。また、記憶装置１００１として、データを書き換え可能な記憶媒体が用いられる場合には、記憶装置１００１は、ＳＡＲ画像格納部２００および重み行列格納部３００としても使用可能である。
また、プログラムは、様々なタイプの一時的なコンピュータ可読媒体（transitory computer readable medium ）に格納されてもよい。一時的なコンピュータ可読媒体には、例えば、有線通信路または無線通信路を介して、すなわち、電気信号、光信号または電磁波を介して、プログラムが供給される。
メモリ１００２は、例えばＲＡＭ（Random Access Memory）で実現され、ＣＰＵ１０００が処理を実行するときに一時的にデータを格納する記憶手段である。メモリ１００２に、記憶装置１００１または一時的なコンピュータ可読媒体が保持するプログラムが転送され、ＣＰＵ１０００がメモリ１００２内のプログラムに基づいて処理を実行するような形態も想定しうる。
図１３は、画像解析装置の主要部を示すブロック図である。図１３に示す画像解析装置１００は、同一地域が記録された複数の画像における全ての画像の組での位相差が反映された複素行列を算出する複素行列算出部（複素行列算出手段）１１（実施形態では、複素行列算出部１１０で実現される。）と、位相の変化を説明するパラメータの候補を複数選定するパラメータ候補選定部（パラメータ候補選定手段）１２（実施形態では、パラメータ候補選定部１２０で実現される。）と、複素行列と所定の重み行列とを用いて、複数の候補の尤もらしさを評価する候補評価部（候補評価手段）１３（実施形態では、候補評価部１３０で実現される。）と、パラメータの候補を尤もらしさによって重み付けして、パラメータの候補の統計値を算出する統計値算出部（統計値算出手段）１４（実施形態では、統計値算出部１４０で実現される。）とを備えている。
上記の実施形態の一部または全部は、以下の付記のようにも記載され得るが、以下に限定されるわけではない。
（付記１）同一地域が記録された複数の画像における全ての画像の組での位相差が反映された複素行列を算出する複素行列算出手段と、
位相の変化を説明するパラメータの候補を複数選定するパラメータ候補選定手段と、
前記複素行列と所定の重み行列とを用いて、前記複数の候補の尤もらしさを評価する候補評価手段と、
前記パラメータの候補を尤もらしさによって重み付けして、前記パラメータの候補の統計値を算出する統計値算出手段と
を備える画像解析装置。
（付記２）前記複素行列と前記重み行列とを乗算し、乗算結果を前記候補評価手段に対して出力する乗算手段（実施形態では、乗算部１５０で実現される。）を備える
付記１の画像解析装置。
（付記３）パラメータの候補のそれぞれについて、各画像における位相の値を有するパラメータ複素行列を出力するようなパラメータ複素行列算出手段（実施形態では、パラメータ複素行列算出部１６０で実現される。）をさらに備え、
前記統計値算出手段は、前記パラメータ複素行列の平均値を算出し、
前記候補評価手段は、前記パラメータ複素行列の平均値も用いて前記複数の候補の尤もらしさを評価する
付記１の画像解析装置。
（付記４）前記複素行列と前記重み行列と前記パラメータ複素行列の平均値とを乗算し、乗算結果を前記候補評価手段に対して出力する乗算手段（実施形態では、乗算器１５０で実現される。）を備える
付記３の画像解析装置。
（付記５）前記パラメータ候補選定手段は、複数のパラメータの各々について候補を選定し、
前記候補評価手段は、一のパラメータの候補の評価と他のパラメータの候補の評価とを交互に実行し、
前記パラメータ複素行列として、評価対象のパラメータとは異なるパラメータについてのパラメータ複素行列を使用する
付記３または付記４の画像解析装置。
（付記６）前記候補評価手段は、
前記パラメータの候補を位相に変換する変換手段（実施形態では、変換部１３２で実現される。）と、
前記変換手段によって得られた位相を評価して第１の評価値を算出する第１の評価手段（実施形態では、評価部１３４で実現される。）と、
画像の各々に対応して、ランダムな位相の組を１つ以上生成するランダム位相生成手段（実施形態では、ランダム位相生成部１３３で実現される。）と、
ランダムな位相の組を評価して第２の評価値を算出する第２の評価手段（実施形態では、評価部１３５で実現される。）と、
前記第２の評価値を統合する統合手段（実施形態では、統合部１３６で実現される。）と、
前記第１の評価値が、前記統合手段が統合した第２の評価値以下である場合に、前記第１の評価値を前記統合手段が統合した第２の評価値に置き換える比較手段（実施形態では、比較部１３７で実現される。）とを含む
付記１から付記５のうちのいずれかの画像解析装置。
（付記７）前記パラメータの候補の統計値は、前記パラメータの候補の平均値および分散を含む
付記１から付記６のうちのいずれかの画像解析装置。
（付記８）同一地域が記録された複数の画像における全ての画像の組での位相差が反映された複素行列を算出し、
位相の変化を説明するパラメータの候補を複数選定し、
前記複素行列と所定の重み行列とを用いて、前記複数の候補の尤もらしさを評価し、
前記パラメータの候補を尤もらしさによって重み付けして、前記パラメータの候補の統計値を算出する
画像解析方法。
（付記９）前記パラメータの候補のそれぞれについて、各画像における位相の値を有するパラメータ複素行列を出力し、
前記統計値を算出するときに、前記パラメータ複素行列の平均値を算出し、
前記パラメータ複素行列の平均値も用いて前記複数の候補の尤もらしさを評価する
付記８の画像解析方法。
（付記１０）複数のパラメータの各々について候補を選定し、
一のパラメータの候補の評価と他のパラメータの候補の評価とを交互に実行し、
前記パラメータ複素行列として、評価対象のパラメータとは異なるパラメータについてのパラメータ複素行列を使用する
付記９の画像解析方法。
（付記１１）前記複数の候補の尤もらしさを評価するときに、
前記パラメータの候補を位相に変換し、
前記位相を評価して第１の評価値を算出し、
画像の各々に対応して、ランダムな位相の組を１つ以上生成し、
ランダムな位相の組を評価して第２の評価値を算出し、
前記第２の評価値を統合し、
前記第１の評価値が、統合された前記第２の評価値以下である場合に、前記第１の評価値を統合された前記第２の評価値に置き換える
付記８から付記１０のうちのいずれかの画像解析方法。
（付記１２）コンピュータに、
同一地域が記録された複数の画像における全ての画像の組での位相差が反映された複素行列を算出する処理と、
位相の変化を説明するパラメータの候補を複数選定する処理と、
前記複素行列と所定の重み行列とを用いて、前記複数の候補の尤もらしさを評価する処理と、
前記パラメータの候補を尤もらしさによって重み付けして、前記パラメータの候補の統計値を算出する処理と
を実行させる画像解析プログラムが格納されたコンピュータ読み取り可能な記録媒体。
（付記１３）コンピュータに、
前記パラメータの候補のそれぞれについて、各画像における位相の値を有するパラメータ複素行列を出力する処理と、
前記統計値を算出するときに、前記パラメータ複素行列の平均値を算出する処理と、
前記パラメータ複素行列の平均値も用いて前記複数の候補の尤もらしさを評価する処理と
を実行させる画像解析プログラムが格納された付記１２のコンピュータ読み取り可能な記録媒体。
（付記１４）コンピュータに、
複数のパラメータの各々について候補を選定する処理を実行させ、
一のパラメータの候補の評価と他のパラメータの候補の評価とを交互に実行させ、
前記パラメータ複素行列として、評価対象のパラメータとは異なるパラメータについてのパラメータ複素行列を使用させる
画像解析プログラムが格納された付記１３のコンピュータ読み取り可能な記録媒体。
（付記１５）コンピュータに、
前記パラメータの候補のそれぞれについて、各画像における位相の値を有するパラメータ複素行列を出力する処理と、
前記統計値を算出するときに、前記パラメータ複素行列の平均値を算出する処理と、
前記パラメータ複素行列の平均値も用いて前記複数の候補の尤もらしさを評価する処理と
を実行させるための画像解析プログラム。
（付記１６）コンピュータに、
前記パラメータの候補のそれぞれについて、各画像における位相の値を有するパラメータ複素行列を出力する処理と、
前記統計値を算出するときに、前記パラメータ複素行列の平均値を算出する処理と、
前記パラメータ複素行列の平均値も用いて前記複数の候補の尤もらしさを評価する処理と
を実行させる付記１５の画像解析プログラム。
（付記１７）コンピュータに、
複数のパラメータの各々について候補を選定する処理を実行させ、
一のパラメータの候補の評価と他のパラメータの候補の評価とを交互に実行させ、
前記パラメータ複素行列として、評価対象のパラメータとは異なるパラメータについてのパラメータ複素行列を使用させる
付記１６の画像解析プログラム。
以上、実施形態を参照して本願発明を説明したが、本願発明は上記の実施形態に限定されるものではない。本願発明の構成や詳細には、本願発明のスコープ内で当業者が理解し得る様々な変更をすることができる。
１１ 複素行列算出部
１２ パラメータ候補選定部
１３ 候補評価部
１４ 統計値算出部
１０，２０，４０ 画像解析装置
１００ 画像解析装置
１１０ 複素行列算出部
１２０ パラメータ候補選定部
１３０，１３１ 候補評価部
１３２ 変換部
１３３ ランダム位相生成部
１３４，１３５ 評価部
１３６ 統合部
１３７ 比較部
１４０，１４１ 統計値算出部
１５０ 乗算部
１６０ パラメータ複素行列算出部
１７０ 低周波フィルタ部
１７１ 画素選択部
１７２ 低次元圧縮部
１７３ 展開部
１７４ 最適化部
２００ ＳＡＲ画像格納部
３００ 重み行列格納部
１０００ ＣＰＵ
１００１ 記憶装置
１００２ メモリ

Claims (10)

  1. 同一地域が記録された複数の画像における全ての画像の組での位相差が反映された複素行列を算出する複素行列算出手段と、
    位相の変化を説明するパラメータの候補を複数選定するパラメータ候補選定手段と、
    前記複素行列と所定の重み行列とを用いて、前記複数の候補の尤もらしさを評価する候補評価手段と、
    前記パラメータの候補を尤もらしさによって重み付けして、前記パラメータの候補の統計値を算出する統計値算出手段と
    を備える画像解析装置。
  2. 前記複素行列と前記重み行列とを乗算し、乗算結果を前記候補評価手段に対して出力する乗算手段を備える
    請求項１に記載の画像解析装置。
  3. 前記パラメータの候補のそれぞれについて、各画像における位相の値を有するパラメータ複素行列を出力するようなパラメータ複素行列算出手段をさらに備え、
    前記統計値算出手段は、前記パラメータ複素行列の平均値を算出し、
    前記候補評価手段は、前記パラメータ複素行列の平均値も用いて前記複数の候補の尤もらしさを評価する
    請求項１に記載の画像解析装置。
  4. 前記複素行列と前記重み行列と前記パラメータ複素行列の平均値とを乗算し、乗算結果を前記候補評価手段に対して出力する乗算手段を備える
    請求項３に記載の画像解析装置。
  5. 前記パラメータ候補選定手段は、複数のパラメータの各々について候補を選定し、
    前記候補評価手段は、一のパラメータの候補の評価と他のパラメータの候補の評価とを交互に実行し、
    前記パラメータ複素行列として、評価対象のパラメータとは異なるパラメータについてのパラメータ複素行列を使用する
    請求項３または請求項４に記載の画像解析装置。
  6. 前記候補評価手段は、
    前記パラメータの候補を位相に変換する変換手段と、
    前記変換手段によって得られた位相を評価して第１の評価値を算出する第１の評価手段と、
    画像の各々に対応して、ランダムな位相の組を１つ以上生成するランダム位相生成手段と、
    前記ランダムな位相の組を評価して第２の評価値を算出する第２の評価手段と、
    前記第２の評価値を統合する統合手段と、
    前記第１の評価値が、前記統合手段が統合した第２の評価値以下である場合に、前記第１の評価値を前記統合手段が統合した第２の評価値に置き換える比較手段とを含む
    請求項１から請求項５のうちのいずれか１項に記載の画像解析装置。
  7. 前記パラメータの候補の統計値は、前記パラメータの候補の平均値および分散を含む
    請求項１から請求項６のうちのいずれか１項に記載の画像解析装置。
  8. 同一地域が記録された複数の画像における全ての画像の組での位相差が反映された複素行列を算出し、
    位相の変化を説明するパラメータの候補を複数選定し、
    前記複素行列と所定の重み行列とを用いて、前記複数の候補の尤もらしさを評価し、
    前記パラメータの候補を尤もらしさによって重み付けして、前記パラメータの候補の統計値を算出する
    画像解析方法。
  9. 前記パラメータの候補のそれぞれについて、各画像における位相の値を有するパラメータ複素行列を出力し、
    前記統計値を算出するときに、前記パラメータ複素行列の平均値を算出し、
    前記パラメータ複素行列の平均値も用いて前記複数の候補の尤もらしさを評価する
    請求項８に記載の画像解析方法。
  10. コンピュータに、
    同一地域が記録された複数の画像における全ての画像の組での位相差が反映された複素行列を算出する処理と、
    位相の変化を説明するパラメータの候補を複数選定する処理と、
    前記複素行列と所定の重み行列とを用いて、前記複数の候補の尤もらしさを評価する処理と、
    前記パラメータの候補を尤もらしさによって重み付けして、前記パラメータの候補の統計値を算出する処理と
    を実行させるための画像解析プログラム。
JP2022553250A 2020-09-29 2020-09-29 画像解析装置および画像解析方法 Active JP7384295B2 (ja)

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
PCT/JP2020/036844 WO2022070243A1 (ja) 2020-09-29 2020-09-29 画像解析装置および画像解析方法

Publications (3)

Publication Number Publication Date
JPWO2022070243A1 JPWO2022070243A1 (ja) 2022-04-07
JPWO2022070243A5 JPWO2022070243A5 (ja) 2023-06-01
JP7384295B2 true JP7384295B2 (ja) 2023-11-21

Family

ID=80951452

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2022553250A Active JP7384295B2 (ja) 2020-09-29 2020-09-29 画像解析装置および画像解析方法

Country Status (3)

Country Link
US (1) US20230341544A1 (ja)
JP (1) JP7384295B2 (ja)
WO (1) WO2022070243A1 (ja)

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2012533051A (ja) 2009-07-08 2012-12-20 テレ−リレヴァメント エウローパ−ティ．エルレ．エ． エッセ．エルレ．エルレ． 同一領域上で取得されたｓａｒ画像から得られるインターフェログラムのフィルタリング処理方法
CN105717480A (zh) 2016-01-25 2016-06-29 西安电子科技大学 基于相位差分的单站无源快速定位方法
JP2017072473A (ja) 2015-10-07 2017-04-13 三菱電機株式会社 画像処理装置および画像処理方法
WO2019215819A1 (ja) 2018-05-08 2019-11-14 日本電気株式会社 合成開口レーダ画像解析システム、合成開口レーダ画像解析方法および合成開口レーダ画像解析プログラム

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2012533051A (ja) 2009-07-08 2012-12-20 テレ−リレヴァメント エウローパ−ティ．エルレ．エ． エッセ．エルレ．エルレ． 同一領域上で取得されたｓａｒ画像から得られるインターフェログラムのフィルタリング処理方法
JP2017072473A (ja) 2015-10-07 2017-04-13 三菱電機株式会社 画像処理装置および画像処理方法
CN105717480A (zh) 2016-01-25 2016-06-29 西安电子科技大学 基于相位差分的单站无源快速定位方法
WO2019215819A1 (ja) 2018-05-08 2019-11-14 日本電気株式会社 合成開口レーダ画像解析システム、合成開口レーダ画像解析方法および合成開口レーダ画像解析プログラム

Also Published As

Publication number Publication date
WO2022070243A1 (ja) 2022-04-07
JPWO2022070243A1 (ja) 2022-04-07
US20230341544A1 (en) 2023-10-26

Similar Documents

Publication Publication Date Title
RU2503938C2 (ru) Способ и устройство частотного анализа данных
US9223021B2 (en) Method and system for motion compensated target detection using acoustical focusing
Shkvarko Estimation of wavefield power distribution in the remotely sensed environment: Bayesian maximum entropy approach
Lazanu et al. Matter bispectrum of large-scale structure: Three-dimensional comparison between theoretical models and numerical simulations
JP7248109B2 (ja) 合成開口レーダの信号処理装置および信号処理方法
US11709254B2 (en) SAR image analysis system, image processing apparatus, image processing method, and image processing program
CN106772273B (zh) 一种基于动态孔径的sar虚假目标干扰抑制方法及系统
JP5334606B2 (ja) レーダ画像信号処理装置
JP7384295B2 (ja) 画像解析装置および画像解析方法
US9147126B2 (en) Post compression detection (PoCoDe)
JP2012127920A (ja) 目標類識別装置
JP7188594B2 (ja) 画像処理装置および画像処理方法
CN101688772A (zh) 使用空间内插和掩蔽区来估计属性值
KR102279563B1 (ko) 레이더 신호 처리 장치 및 레이더 신호 처리 방법
JP7722453B2 (ja) 画像解析装置および画像解析方法
JP7416277B2 (ja) 画像解析装置および画像解析方法
US12399269B2 (en) Image analyzing device and image analyzing method
Xiong Video-based vibration analysis for structural health monitoring in civil engineering
CN117930238A (zh) 基于最小生成树的sar层析参考网生成方法及装置
JP7768334B2 (ja) 画像解析システムおよび画像解析方法
WO2019123786A1 (ja) 合成開口レーダの信号処理装置及び信号処理方法
Liu et al. On the Value of Terrain Classification Using SAR Altimeter Delay Doppler Image
RU2775154C1 (ru) Способ однолучевого измерения высоты и составляющих скорости воздушного судна
WO2021186557A1 (ja) データ処理装置およびデータ処理方法
Toth et al. Spatial spectrum analysis of various digital elevation models

Legal Events

Date Code Title Description
A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20230313

A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20230313

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: 20231010

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20231023

R151 Written notification of patent or utility model registration

Ref document number: 7384295

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R151