JP6349938B2 - 測定点情報提供装置、変動検出装置、方法およびプログラム - Google Patents

測定点情報提供装置、変動検出装置、方法およびプログラム Download PDF

Info

Publication number
JP6349938B2
JP6349938B2 JP2014097644A JP2014097644A JP6349938B2 JP 6349938 B2 JP6349938 B2 JP 6349938B2 JP 2014097644 A JP2014097644 A JP 2014097644A JP 2014097644 A JP2014097644 A JP 2014097644A JP 6349938 B2 JP6349938 B2 JP 6349938B2
Authority
JP
Japan
Prior art keywords
measurement point
polarization
candidate
input
information
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
JP2014097644A
Other languages
English (en)
Other versions
JP2015215211A (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
Priority to JP2014097644A priority Critical patent/JP6349938B2/ja
Priority to EP15165444.9A priority patent/EP2942638B1/en
Priority to US14/704,097 priority patent/US10338215B2/en
Publication of JP2015215211A publication Critical patent/JP2015215211A/ja
Application granted granted Critical
Publication of JP6349938B2 publication Critical patent/JP6349938B2/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
    • G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/02—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
    • G01S7/024—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00 using polarisation effects
    • G01S7/025—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00 using polarisation effects involving the transmission of linearly polarised waves
    • 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
    • G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/02—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
    • G01S7/41—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00 using analysis of echo signal for target characterisation; Target signature; Target cross-section
    • G01S7/411—Identification of targets based on measurements of radar reflectivity
    • 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/904—SAR modes
    • G01S13/9076—Polarimetric features in SAR

Description

本発明は、変位を測定可能な測定点の情報をユーザに提供する測定点情報提供装置、地表や物体の変動を検出するための変動検出装置、測定点情報提供方法、変動検出方法、測定点情報提供用プログラムおよび変動検出用プログラムに関する。
地表や物体の変動を計算する技術の１つに、ＰＳ−ＩｎＳＡＲ（Permanent/Persistent Scatters Interferometric Synthetic Aperture Radar）と呼ばれる技術がある（例えば、非特許文献１および非特許文献２）。ＰＳ−ＩｎＳＡＲは、ＰＳ（Permanent/Persistent Scatters）特性を有する点に対して合成開口レーダ（Synthetic Aperture Radar：ＳＡＲ）が取得したデータであるＳＡＲデータに、干渉（Interferometry）技術を適用して、地表上やある物体上の点における変位を測定する技術である。電波は、光波と異なり、雲や雨などを通過するため、天候や夜間に関係なく観測できるという特徴を有している。以下、ＰＳ−ＩｎＳＡＲにおいて変位を測定する地表上や物体上の点を、測定点という。なお、地表上や物体上の各点は、合成開口レーダの撮影視野の各分解能セルに対応している。
ここで、ＰＳ特性は、電波の散乱特性が時間経過によって変化しない特性である。植物や波などはＰＳ特性を有さない例であるが、大抵の人工構造物はその多くの場所がＰＳ特性を有している。しかし、ＰＳ−ＩｎＳＡＲにおいて、ある場所を測定点として用いるためには、ＰＳ特性に加えて、ある程度の後方散乱強度（以下、反射強度という）が必要である。これは、反射強度が弱いと雑音に信号が埋もれてしまうためである。このように測定点が限られてしまうことがＰＳ−ＩｎＳＡＲの価値を低下させていた。
例えば、一面平板のような形状の箇所では、後方散乱が弱く、ＰＳ特性を有していても測定点にはなりにくいことが知られている。また、例えば、コーナーリフレクタ（Corner Reflector：ＣＲ）に用いられるような、正方形または直角二等辺三角形の板３枚をお互いが直角に向くようにつなぎ合わせたような形状の箇所（例えば、窓枠の四隅等）では、後方散乱が返ってきやすいことから、測定点に適していることが知られている。
また、本発明に関連する技術として、例えば、特許文献１には、偏波特性の異なる複数の電波が反射して得られる散乱波の少なくとも３種類の偏波成分を測定した結果を測定位置ごとに偏波特性データとして予め記憶装置に記憶しておき、記憶された偏波特性データを用いて各測定位置にある物体が人工物か否かを判定する技術が記載されている。
また、例えば、非特許文献３には、多偏波ＳＡＲ画像データの後方散乱係数とされる４つのデータであって、受信偏波の楕円方位角（ellipse orientation angle）ψｒ、受信偏波の楕円率角（ellipticity angle）χｒ、送信偏波の楕円方位角ψｔおよび送信偏波の楕円率角χｒの４つのデータの値の異なる３つのセット（例えば、ＨＨ偏波、ＨＶ偏波、ＶＶ偏波）を赤、緑、青に対応させて、多偏波ＳＡＲ画像データの疑似カラー合成表示を行うことが記載されている。また、非特許文献３には、多偏波ＳＡＲ画像データの利用例の１つとして、２つの対象物のコントラストを最大にする送受信条件を求めて、これにより、識別したい特定の対象物とその他の対象物との画素値の比を最大にして表示する例が示されている。
特許第５３０５９８５号公報
Alessandro Ferretti, Claudio Prati, and Fabio Rocca, "Nonliner subsidencd Rate Estimation Using Permanent Svatterers in Differntial SAR Interferometry", IEEE TRANSACTION ON GEOSCIENCE AND REMOTE SENSING, VOL.38, No.5, Sep 2000. Alessandro Ferretti, Claudio Prati, and Fabio Rocca, "Permanent Scatteres in SAR Interferometry", IEEE TRANSACTION ON GEOSCIENCE AND REMOTE SENSING, VOL.39, No.1, Jan 2001. 寛渕 哲也，"グランドトルースデータと光学センサ画像データで検証した災害観測手法としての多周波・多偏波ＳＡＲ画像データの有用性"，防災科学技術研究所研究報告 第６３号，２００２年６月．
ところで、実際のＰＳ−ＩｎＳＡＲ処理では、決められた偏波方向のＳＡＲ画像だけを用いて、地表上または物体上の点（以下、地点という）における変位を測定している。具体的には、代表的な送受信偏波の組合せであるＨＨ偏波、ＨＶ偏波、ＶＨ偏波、ＶＶ偏波の中から一つの送受信偏波の組合せを選び、選んだ送受信偏波の組合せによるＳＡＲ画像を、変動を検出したい期間の前後に渡る複数の時刻分取得して、変動検出のための演算処理が行われている。
ここで、Ｈ（Horizontal）偏波は、水平偏波すなわち電界の振動方向である偏波方向が水平面と平行である偏波を表している。また、Ｖ（Vertival）偏波は、垂直偏波すなわち偏波方向が水平面と垂直である偏波を表している。また、送受信偏波の組合せにおける最初のアルファベットは送信の偏波方向を表し、２つ目のアルファベットは受信の偏波方向を表わしている。したがって、ＨＨ偏波は、水平偏波の電波を送信して水平偏波の電波を受信することを表している。また、ＨＶ偏波は、水平偏波の電波を送信して垂直偏波の電波を受信することを表している。また、ＶＨ偏波は、垂直偏波の電波を送信して水平偏波の電波を受信することを表している。また、ＶＶ偏波は、垂直偏波の電波を送信して垂直偏波の電波を受信することを表している。
一般に、ある偏波方向の電磁波がある地点に照射され、その地点にある構造物に反射して戻ってくるとき、その地点にある構造物の形状や材質等によっては、偏波方向が回転して戻ってくる。そのときの偏波の回転角は、その地点にある構造物の形状や材質等によって異なり、事前に予測するのは困難である。
したがって、道路や鉄道、橋といった広い範囲に渡る人工構造物の歪みや剥がれ等による経年劣化を、ＰＳ点における微小な変動の検出により測定することを考えた場合に、１つの送受信偏波の組合せによるＳＡＲ画像を用いるだけでは、測定可能な点が少なく、十分な精度を得ることができないという問題があった。これは、好適な送受信偏波の組合せは場所ごとに異なるが、１つの送受信偏波の組合せによるＳＡＲ画像しか用いない場合には、他の送受信偏波の組合せによるＳＡＲ画像で反射強度が強い点があってもその点を測定点として用いることができないからである。
なお、特許文献１に記載されている技術は、散乱成分分析により人工物と自然物とを区別するために、測定位置ごとに少なくとも３種類の偏波成分を測定した結果である偏波特性データを用いているものである。そこには、偏波特性データを用いて、変位を測定可能な地点をいかに多くまたは見逃さずに検出して、その情報をユーザに提供しようといったことについては何ら考慮されていない。
また、非特許文献３に記載されている技術は、３種類の偏波成分による後方散乱強度の異なり具合を色相によって判別できるようにしよう、また、反射機構の異なる物質の差を画像上で強調して表示しようというものである。そこには、後方散乱係数を利用して、変位を測定可能な地点をいかに多くまたは見逃さずに検出して、その情報をユーザに提供しようといったことについては何ら考慮されていない。
ユーザに多くのまたは漏れのない測定可能な地点の情報を提供できれば、ユーザがその情報を基に測定点を指定することができ、その結果、変位を測定できる測定点を増やすことができる。また、変位を測定できる測定点が増えれば、変動検出の感度を上げることができ、検出精度の向上に寄与できるものと考えられる。また、測定可能な地点の情報と変動の検出結果とを照らし合わせれば、ある地点に対して変位がなかったのか測定不能なのかを判別できるなど、高感度かつより正確な検出結果が得られるものと考えられる。
そこで、本発明は、変動検出の感度を向上させることが可能な測定点情報提供装置、変動検出装置、測定点情報提供方法、変動検出方法、測定点情報提供用プログラムおよび変動検出用プログラムを提供することを目的とする。
本発明による測定点情報提供装置は、特定領域を含む撮影視野内の分解能セルに対応するピクセルの各々に対応づけて、略同一時刻に合成開口レーダが観測した送受信偏波の組合せがＨＨ偏波、ＨＶ偏波、ＶＨ偏波およびＶＶ偏波である基本４偏波組による観測データから生成される反射強度および位相を示す情報を少なくとも保持するＳＡＲ画像組を入力するＳＡＲ画像組入力手段と、入力されたＳＡＲ画像組を用いて、対象ピクセルの各々に対して、反射強度が所定値以上または最も強い偏波組を決定する偏波組決定手段と、偏波組決定手段により決定された偏波組に基づいて、少なくとも対象ピクセルの決定された偏波組を示す情報を含む測定点情報を生成する測定点情報生成手段と、測定点の候補を入力する測定点候補入力手段と、測定点情報に基づいて、入力された測定点の候補が測定点として適するか否かを、当該候補に対応するピクセルの決定された偏波組による反射強度を少なくとも用いて判定する測定点条件判定手段と、測定点条件判定手段による判定結果を出力する判定結果出力手段とを備え、偏波組決定手段は、入力された測定点の候補に対応するピクセルである対象ピクセルの各々に対して、反射強度が所定値以上または最も強い偏波組を決定し、測定点情報生成手段は、偏波組決定手段による偏波組の決定結果に基づき、少なくとも対象ピクセルの決定された偏波組を示す情報と、当該対象ピクセルの偏波組による反射強度を示す情報とを含む測定点情報を生成することを特徴とする。
また、本発明による変動検出装置は、特定領域を含む撮影視野内の分解能セルに対応するピクセルの各々に対応づけて、略同一時刻に合成開口レーダが観測した送受信偏波の組合せがＨＨ偏波、ＨＶ偏波、ＶＨ偏波およびＶＶ偏波である基本４偏波組による観測データから生成される反射強度および位相を示す情報を少なくとも保持するＳＡＲ画像組を２組以上入力するＳＡＲ画像組入力手段と、測定点の候補を入力する測定点候補入力手段と、入力されたＳＡＲ画像組のうちの少なくとも１組を用いて、入力された測定点の候補に対応するピクセルである対象ピクセルの各々に対して、反射強度が所定値以上または最も強い偏波組を決定する偏波組決定手段と、偏波組決定手段により決定された偏波組に基づいて、少なくとも対象ピクセルの決定された偏波組を示す情報と、当該対象ピクセルの偏波組による反射強度を示す情報とを含む測定点情報を生成する測定点情報生成手段と、測定点情報に基づいて、入力された測定の候補が測定点として適するか否かを、当該候補に対応するピクセルの決定された偏波組による反射強度を少なくとも用いて判定する測定点条件判定手段と、測定点条件判定手段による判定結果を出力する判定結果出力手段と、測定点として適すると判定された候補、または出力後に測定点に指定された候補を測定点に用いて、測定点情報と、入力されたＳＡＲ画像組とに基づいて、測定点における変位を測定する変位測定手段と、変位測定手段による測定結果を示す情報を出力する第１の出力手段と、測定点情報、または測定点情報に基づいて生成される、対象ピクセルに対応する地点が測定点に適するか否かもしくはどの程度測定点に適するかを示す情報を出力する第２の出力手段とを備えたことを特徴とする。
また、本発明による変動検出装置は、特定領域を含む撮影視野内の分解能セルに対応するピクセルの各々に対応づけて、略同一時刻に合成開口レーダが観測した送受信偏波の組合せがＨＨ偏波、ＨＶ偏波、ＶＨ偏波およびＶＶ偏波である基本４偏波組による観測データから生成される反射強度および位相を示す情報を少なくとも保持するＳＡＲ画像組を２組以上入力するＳＡＲ画像組入力手段と、測定点の候補を入力する測定点候補入力手段と、入力されたＳＡＲ画像組のうちの少なくとも１組を用いて、入力された測定点の候補に対応するピクセルである対象ピクセルの各々に対して、反射強度が所定値以上または最も強い偏波組を決定する偏波組決定手段と、偏波組決定手段により決定された偏波組に基づいて、少なくとも対象ピクセルの決定された偏波組を示す情報と、当該対象ピクセルの偏波組による反射強度を示す情報とを含む測定点情報を生成する測定点情報生成手段と、測定点情報に基づいて、入力された測定点の候補が測定点として適するか否かを、当該候補に対応するピクセルの決定された偏波組による反射強度を少なくとも用いて判定する測定点条件判定手段と、測定点条件判定手段による判定結果を出力する判定結果出力手段と、測定点条件判定手段によって測定点として適すると判定された候補を測定点に用いて、測定点情報と、入力されたＳＡＲ画像組とに基づいて、測定点における変位を測定する変位測定手段とを備えたことを特徴とする。
また、本発明による測定点情報提供方法は、特定領域を含む撮影視野内の分解能セルに対応するピクセルの各々に対応づけて、略同一時刻に合成開口レーダが観測した送受信偏波の組合せがＨＨ偏波、ＨＶ偏波、ＶＨ偏波およびＶＶ偏波である基本４偏波組による観測データから生成される反射強度および位相を示す情報を少なくとも保持するＳＡＲ画像組を入力し、測定点の候補を入力し、入力されたＳＡＲ画像組のうちの少なくとも１組を用いて、入力された測定点の候補に対応するピクセルである対象ピクセルの各々に対して、反射強度が所定値以上または最も強い偏波組を決定し、偏波組決定手段により決定された偏波組に基づいて、少なくとも対象ピクセルの決定された偏波組を示す情報と、当該対象ピクセルの偏波組による反射強度を示す情報とを含む測定点情報を生成し、測定点情報に基づいて、入力された測定の候補が測定点として適するか否かを、当該候補に対応するピクセルの決定された偏波組による反射強度を少なくとも用いて判定し、判定結果を出力することを特徴とする。
また、本発明による変動検出方法は、特定領域を含む撮影視野内の分解能セルに対応するピクセルの各々に対応づけて、略同一時刻に合成開口レーダが観測した送受信偏波の組合せがＨＨ偏波、ＨＶ偏波、ＶＨ偏波およびＶＶ偏波である基本４偏波組による観測データから生成される反射強度および位相を示す情報を少なくとも保持するＳＡＲ画像組を２組以上入力し、測定点の候補を入力し、入力されたＳＡＲ画像組を用いて、入力された測定点の候補に対応するピクセルである対象ピクセルの各々に対して、反射強度が所定値以上または最も強い偏波組を決定し、決定された偏波組に基づいて、少なくとも対象ピクセルの決定された偏波組を示す情報と、当該対象ピクセルの偏波組による反射強度を示す情報とを含む測定点情報を生成し、測定点情報に基づいて、入力された測定点の候補が測定点として適するか否かを、当該候補に対応するピクセルの偏波組による反射強度を少なくとも用いて判定し、判定結果を出力し、測定点として適すると判定された候補、または出力後に測定点に指定された候補を測定点に用いて、測定点情報と、入力されたＳＡＲ画像組とに基づいて、測定点における変位を測定し、測定点における変位の測定結果を示す情報とともに、測定点情報、または測定点情報に基づいて生成される、対象ピクセルに対応する地点が測定点に適するか否かもしくはどの程度測定点に適するかを示す情報を出力することを特徴とする。
また、本発明による変動検出方法は、特定領域を含む撮影視野内の分解能セルに対応するピクセルの各々に対応づけて、略同一時刻に合成開口レーダが観測した送受信偏波の組合せがＨＨ偏波、ＨＶ偏波、ＶＨ偏波およびＶＶ偏波である基本４偏波組による観測データから生成される反射強度および位相を示す情報を少なくとも保持するＳＡＲ画像組を２組以上入力し、測定点の候補を入力し、入力されたＳＡＲ画像組を用いて、入力された測定点の候補に対応するピクセルである対象ピクセルの各々に対して、反射強度が所定値以上または最も強い偏波組を決定し、決定された偏波組に基づいて、少なくとも対象ピクセルの決定された偏波組を示す情報と、当該対象ピクセルの偏波組による反射強度を示す情報とを含む測定点情報を生成し、測定点情報に基づいて、入力された測定点の候補が測定点として適するか否かを、当該候補に対応するピクセルの偏波組による反射強度を少なくとも用いて判定し、判定により測定点として適さないと判定された場合に、その旨を出力し、判定により測定点として適すると判定された場合に、当該候補を測定点に用いて、測定点情報と、入力されたＳＡＲ画像組とに基づいて、測定点における変位を測定してもよい。
また、本発明による測定点情報提供用プログラムは、コンピュータに、特定領域を含む撮影視野内の分解能セルに対応するピクセルの各々に対応づけて、略同一時刻に合成開口レーダが観測した送受信偏波の組合せがＨＨ偏波、ＨＶ偏波、ＶＨ偏波およびＶＶ偏波である基本４偏波組による観測データから生成される反射強度および位相を示す情報を少なくとも保持するＳＡＲ画像組を入力するＳＡＲ画像組入力処理、測定点の候補を入力する測定点候補入力処理、入力されたＳＡＲ画像組を用いて、入力された測定点の候補に対応するピクセルである対象ピクセルの各々に対して、反射強度が所定値以上または最も強い偏波組を決定する偏波組決定処理、偏波組決定処理で決定された偏波組に基づいて、少なくとも対象ピクセルの決定された偏波組を示す情報と、当該対象ピクセルの偏波組による反射強度を示す情報とを含む測定点情報を生成する測定点情報生成処理、測定点情報に基づいて、入力された測定点の候補が測定点として適するか否かを、当該候補に対応するピクセルの決定された偏波組による反射強度を少なくとも用いて判定する測定点条件判定処理、および測定点条件判定処理による判定結果を出力する判定結果出力処理を実行させることを特徴とする。
また、本発明による変動検出用プログラムは、コンピュータに、特定領域を含む撮影視野内の分解能セルに対応するピクセルの各々に対応づけて、略同一時刻に合成開口レーダが観測した送受信偏波の組合せがＨＨ偏波、ＨＶ偏波、ＶＨ偏波およびＶＶ偏波である基本４偏波組による観測データから生成される反射強度および位相を示す情報を少なくとも保持するＳＡＲ画像組を２組以上入力するＳＡＲ画像組入力処理、測定点の候補を入力する測定点候補入力処理、入力されたＳＡＲ画像組を用いて、入力された測定点の候補に対応するピクセルである対象ピクセルの各々に対して、反射強度が所定値以上または最も強い偏波組を決定する偏波組決定処理、偏波組決定処理で決定された偏波組に基づいて、少なくとも対象ピクセルの決定された偏波組を示す情報と、当該対象ピクセルの偏波組による反射強度を示す情報とを含む測定点情報を生成する測定点情報生成処理、測定点情報に基づいて、入力された測定点の候補が測定点として適するか否かを、当該候補に対応するピクセルの決定された偏波組による反射強度を少なくとも用いて判定する測定点条件判定処理、測定点条件判定処理による判定結果を出力する判定結果出力処理、測定点として適すると判定された候補、または出力後に測定点に指定された候補を測定点に用いて、測定点情報と、入力されたＳＡＲ画像組とに基づいて、測定点における変位を測定する変位測定処理、測定点における変位の測定結果を示す情報を出力する第１の出力処理、および測定点情報、または測定点情報に基づいて生成される、対象ピクセルに対応する地点における地点が測定点に適するか否かもしくはどの程度測定点に適するかを示す情報を出力する第２の出力処理を実行させることを特徴とする。
また、本発明による変動検出用プログラムは、コンピュータに、特定領域を含む撮影視野内の分解能セルに対応するピクセルの各々に対応づけて、略同一時刻に合成開口レーダが観測した送受信偏波の組合せがＨＨ偏波、ＨＶ偏波、ＶＨ偏波およびＶＶ偏波である基本４偏波組による観測データから生成される反射強度および位相を示す情報を少なくとも保持するＳＡＲ画像組を２組以上入力するＳＡＲ画像組入力処理、測定点の候補を入力する測定点候補入力処理、入力されたＳＡＲ画像組を用いて、入力された測定点の候補に対応するピクセルである対象ピクセルの各々に対して、反射強度が所定値以上または最も強い偏波組を決定する偏波組決定処理、偏波組決定処理で決定された偏波組に基づいて、少なくとも対象ピクセルの決定された偏波組を示す情報と、当該対象ピクセルの偏波組による反射強度を示す情報とを含む測定点情報を生成する測定点情報生成処理、測定点情報に基づいて、入力された測定点の候補が測定点として適するか否かを、当該候補に対応するピクセルの決定された偏波組による反射強度を少なくとも用いて判定する測定点条件判定処理、測定点条件判定処理による判定結果を出力する判定結果出力処理、および測定点条件判定処理で測定点として適すると判定された候補を測定点に用いて、測定点情報と、入力されたＳＡＲ画像組とに基づいて、測定点における変位を測定する変位測定処理を実行させてもよい。
本発明によれば、変位を測定可能な地点をより多くまたは見逃さずに検出して、その情報をユーザに提供するので、その情報を測定点の指定や測定結果の分析に利用することによって、変動検出の感度を向上させることができる。
第１の実施形態の観測システムの構成例を示すブロック図である。 第１の実施形態の測定点情報提供装置１の構成例を示すブロック図である。 ＳＡＲ画像組入力手段１０１が入力するＳＡＲ画像組の例を示す説明図である。 偏波を説明するための説明図である。 第１の実施形態の測定点情報提供装置１の動作の一例を示すフローチャートである。 第２の実施形態の測定点情報提供装置１の構成例を示すブロック図である。 第２の実施形態の測定点情報提供装置１の動作の一例を示すフローチャートである。 第３の実施形態の観測システムの構成例を示す構成図である。 変動検出装置５の構成例を示すブロック図である。 変動検出手段１０５の一例を示すブロック図である。 変動検出手段１０５の他の例を示すブロック図である。 第３の実施形態の変動検出装置５の動作の一例を示すフローチャートである。 第３の実施形態の変動検出装置５の動作の他の例を示すフローチャートである。 第４の実施形態の変動検出装置５の構成例を示すブロック図である。 第４の実施形態の変動検出装置５の動作の一例を示すフローチャートである。 変動検出手段１０５の他の例を示すブロック図である。 変動検出手段１０５の他の例を示すブロック図である。 本発明の概要を示すブロック図である。
実施形態１．
以下、本発明の実施形態を図面を参照して説明する。図１は、本発明の第１の実施形態の観測システムの構成例を示すブロック図である。図１に示す観測システムは、測定点情報提供装置１と、人工衛星や航空機等の飛翔体２に搭載された合成開口レーダ３とを備える。また、図１において、符号４は合成開口レーダ３によるＳＡＲ画像の撮影視野（撮影領域）を表している。
本実施形態において、合成開口レーダ３は、略同一場所に対して、略同一時刻に、ＨＨ偏波、ＨＶ偏波、ＶＨ偏波、ＶＶ偏波の４種類の送受信偏波の組合せ（以下、基本４偏波組という）による散乱波を観測できる多偏波観測モードを有する。このような多偏波観測モードは、例えば、合成開口レーダ３が少なくとも１つの送信機と２つの受信機とを備え、送信機がＨ偏波とＶ偏波の電波パルスを交互に送信し、２つの受信機（Ｈ偏波に対応した受信機とＶ偏波に対応した受信機）がそれぞれの電波パルスの後方散乱エコーを同時に受信することにより実現できる。このようにして、基本４偏波組による反射強度および相対位相を測定する。ここで、「略同一場所」および「略同一時刻」の「略」の範囲は、校正機能により観測データが同一場所および同一時刻で観測されたとみなせるデータに校正可能な範囲であれば特に問わないが、実際は上記送信機における電波パルスの１周期分の照射場所のずれおよび照射時刻のずれを指す場合が多い。
本発明では、略同一場所に対して、略同一時刻に観測された基本４偏波組による観測データから生成されるＳＡＲ画像またはそれと同等の情報をＳＡＲ画像組という。
一般に、ＳＡＲ画像と呼ばれる情報には、各ピクセルが対応する地点における、受信電界の振幅または反射強度を示す情報と、受信したパルスの位相を示す情報である位相情報とが含まれる。なお、反射強度の平方根が受信電界の振幅である。以下、本発明では、受信電界の振幅を示す情報も反射強度を示す情報として扱うものとする。したがって、ＳＡＲ画像は、各ピクセルに対応づけられるデータ（画素情報とも呼ばれる）として、当該ピクセルに対応する地点における反射強度および位相を示す情報を少なくとも有する画像であると換言できる。ＳＡＲ画像において、各ピクセルに対応づけられるデータは、ＳＡＲデータと呼ばれる。本発明では、ＳＡＲデータを、反射強度を示す情報と、位相を示す情報とを含む情報と定義する。合成開口レーダ３は、ＳＡＲデータを生成するための信号処理で、照射時刻のずれや場所のずれに対して所定の校正処理を行ってもよい。また、ＳＡＲ画像には、付随する情報として、送信および受信に用いた偏波の情報や、観測時間や、入射角や、軌道の情報といった観測条件が含まれていてもよい。
また、ＳＡＲ画像組は、略同一場所に対して、略同一時刻に観測された基本４偏波組による観測データから生成されるＳＡＲ画像またはそれと同等の情報を有するものであれば、データ構造は特に問わない。例えば、ＳＡＲ画像組は、各々が基本４偏波組のいずれかの送受信偏波の組合せに対応した４枚のＳＡＲ画像であってもよいし、各ピクセルに、基本４偏波組による４つのＳＡＲデータ（例えば、４種の反射強度および位相の情報）が対応づけられた１枚のＳＡＲ画像であってもよい。
飛翔体２には、合成開口レーダ３を含む合成開口レーダシステムが搭載されている。合成開口レーダシステムは、例えば、合成開口レーダ３と、信号処理手段と、合成開口レーダ３による観測の結果得られたＳＡＲ画像組を時系列に沿って記憶する記憶手段と、要求に応じて記憶されているＳＡＲ画像組またはその時系列データを所定のフォーマットで出力するＳＡＲ画像出力手段とを含む。
また、図２は、本実施形態の測定点情報提供装置１の構成例を示すブロック図である。図２に示す測定点情報提供装置１は、ＳＡＲ画像組入力手段１０１と、最適偏波組検出手段１０２と、測定点情報記憶手段１０３と、偏波回転処理手段１０４と、測定点情報出力手段１０７とを備える。
ＳＡＲ画像組入力手段１０１は、ユーザが指定した特定領域をその撮影視野に含む少なくとも１組のＳＡＲ画像組であって、略同一場所に対して取得された観測データから生成されるＳＡＲ画像組を少なくとも１組入力する。以下、特定領域を、関心領域（Region of Interest：ＲＯＩ）と呼ぶ場合がある。
図３は、ＳＡＲ画像組入力手段１０１が入力するＳＡＲ画像組の例を示す説明図である。ＳＡＲ画像組入力手段１０１は、例えば、図３に示される観測時刻の異なるｎ組（例えば、数十組）のＳＡＲ画像組のうちの任意の１組を入力してもよいし、これらから選択した２組以上のＳＡＲ画像組を入力してもよい。
ＳＡＲ画像組入力手段１０１は、例えば、飛翔体２に搭載されている合成開口レーダシステムのＳＡＲ画像出力手段が出力するＳＡＲ画像組またはＳＡＲ画像組の時系列データを直接受信して、ＳＡＲ画像組を入力してもよい。そのような場合、ＳＡＲ画像組入力手段１０１は、有線または無線による通信ネットワークを介して合成開口レーダシステムに接続し、所定のフォーマットに従って所望するＳＡＲ画像組またはＳＡＲ画像組の時系列データを要求し、その応答として送信されるＳＡＲ画像組またはＳＡＲ画像組の時系列データを受信してもよい。なお、ＳＡＲ画像組入力手段１０１は、一般的なファイル入力手段であって、ユーザが所定の方法で合成開口レーダシステムに接続して取得したＳＡＲ画像組やＳＡＲ画像組の時系列データやそれらを所定のフォーマットに変換したものなどを、ユーザ操作に応じて入力してもよい。
最適偏波組検出手段１０２は、ＳＡＲ画像組入力手段１０１によって入力された少なくとも１組のＳＡＲ画像組を用いて、対象ピクセルについて、最適な送受信偏波の組合せである最適偏波組を検出する。本実施形態の最適偏波組検出手段１０２は、後述する偏波回転処理手段１０４が提供する偏波回転処理を利用して、ピクセルごとに、最適偏波組を検出する。なお、最適偏波組の検出に用いるＳＡＲ画像組は、どの時刻のＳＡＲ画像組であってもよい。ここで、対象ピクセルは、例えば、関心領域に対応する画像領域に含まれる各ピクセルであったり、ユーザが指定されたピクセルである。
例えば、最適偏波組検出手段１０２は、入力されたＳＡＲ画像組の中から最適偏波組の検出処理に用いるＳＡＲ画像組を１組決定し、決定した１組のＳＡＲ画像組を用いて、各ピクセルに対して、偏波回転処理により、予め定めておいた最適偏波組の各候補による反射強度を算出する。そして、各候補の反射強度の算出結果に基づいて、当該ピクセルにおける最適偏波組を決定してもよい。
ここで、偏波回転処理は、偏波合成とも呼ばれる処理であり、基本４偏波組によるＳＡＲデータ（特に、反射強度と相対位相）から、任意の送受信偏波の組合せによる反射強度を求める処理である。合成開口レーダによって略同一場所および略同一時刻に観測されたＨＨ偏波、ＨＶ偏波、ＶＨ偏波およびＶＶ偏波による反射強度と相対位相とを用いると、撮影視野４の各分解能セルの完全な散乱行列を求めることができ、この散乱行列から各分解能セルの完全な偏波特性を計算することができる。ここで、分解能セルは、合成開口レーダが後方散乱を観測する地表の単位区画であり、ＳＡＲ画像におけるピクセルに相当する。この分解能セルは、合成開口レーダの地上分解能より定まる。各ピクセルの完全な偏波特性が求まれば、それを基に任意の送受信偏波の組合せによる反射強度を求めることができる。
ＳＡＲデータから得られるこのような情報の例として、ポラリメトリックＳＡＲデータが挙げられる。ポラリメトリックＳＡＲデータは、散乱波の散乱時の変化位相を表現したデータである。
以下、偏波回転処理について簡単に説明する。平面波の偏波の一般的表現は楕円偏波であり、図４に示されるように角度χと角度ψで記述できる。ここで、角度χは楕円の丸みを表し、角度ψは長軸の水平方向からの角度を表している。また、散乱波の偏波ベクトルｑｓは、送信アンテナの偏波ベクトルｑｔと散乱行列Ｓを用いて以下の式（１）のように表される。
Figure 0006349938
散乱行列Ｓは、２×２の複素行列であり、入射偏波ベクトルと散乱偏波ベクトルを対応づけている。また、散乱行列Ｓの要素Ｓｊｉは偏波ｉの入射波に対する偏波ｊの散乱波の電界の大きさと位相の変化を表している。この散乱行列Ｓの各要素は、基本４偏波組のＳＡＲデータに相当する。
ｑｒを受信アンテナの偏波ベクトルとすると、受信電界の複素振幅Ｖと散乱パワー（反射強度）Ｐは以下の式（２）のように表される。なお、左肩のｔは転置行列を表し、右肩の＊は複素共役を表している。
Ｖ＝ｑｒ ｔＳｑｔ，
Ｐ＝ＶＶ＊ ・・・（２）
散乱行列Ｓが求まれば、ｑｔとｑｒで規定される任意の送受信偏波の組合せについてその応答（受信電界の複素振幅Ｖや散乱パワーＰ）を計算できる。
なお、偏波ベクトルの代わりにストークスベクトルを用いて入射波と散乱波との対応関係を求めることも可能である。例えば、入射波のストークスベクトルＧｔと、散乱波のストークスベクトルＧｒとは、以下の式（３）に示すような関係にある。
Ｐ＝Ｇｒ ｔＦＧｔ ・・・（３）
ここで、Ｆは、ストークス行列と呼ばれる４×４の実行列である。また、ストークスベクトルＧｔとＧｒとはそれぞれ、以下の式（４）に示される４つのストークスパラメータＧ０，Ｇ１，Ｇ２，Ｇ３からなる４次元ベクトルである。なお、ａｈとａｖはそれぞれ、電界ベクトルの水平成分Ｅｈの振幅と、垂直成分Ｅｖの振幅である。また、δは、Ｅｈの位相δｈとＥｖの位相δｖの差（δｈ−δｖ）である。なお、偏波ベクトルが複素数表現であるのに対して、ストークスベクトルは実数表現である。詳細は省略するが、ストークス行列Ｆの要素は、散乱行列Ｓの要素から周知の関係に従って求められる。
Ｇ０＝ａｈ ２＋ａｖ ２ ，
Ｇ１＝ａｈ ２−ａｖ ２ ，
Ｇ２＝２ａｈａｖｃｏｓδ ，
Ｇ３＝２ａｈａｖｓｉｎδ ・・・（４）
このような散乱行列Ｓまたはストークス行列Ｆから任意の送受信偏波の組合せに対する散乱強度を求める方法には、例えば、次のような方法がある。
まず、求めたい送受信偏波の組における送受信アンテナのストークスベクトルＧｔ，Ｇｒの各要素を、次の式（５）に示すように楕円偏波のχとψで記述する。
Ｇ０＝１．０ ，
Ｇ１＝ｃｏｓ（２χ）ｃｏｓ（２ψ） ，
Ｇ２＝ｃｏｓ（２χ）ｓｉｎ（２ψ） ，
Ｇ３＝ｓｉｎ（２χ） ・・・（５）
なお、χ＝０°は直線偏波、｜χ｜＝４５°は円偏波、χ＜０は右回り、χ＞０は左回りの偏波を表す。また、ψ＝０°は水平偏波、ψ＝９０°垂直偏波を表す。χの変域は−４５°から４５°まで、ψの変域は０°から１８０°までである。例えば、χ＝０°かつψ＝０°は直線の水平偏波を表し、χ＝３５°かつψ＝９０°は左回り楕円で長軸が水平方向の偏波を表す。
次に、対象とするターゲットエリアのピクセルのストークス行列Ｆを各要素ごとに平均し、平均ストークス行列Ｆ’を求める。
次に、送信アンテナのストークスベクトルＧｔと、受信アンテナのストークスベクトルＧｒと、求めた平均ストークス行列Ｆとを、上述の式（３）に代入（ただし、ＦをＦ’とする）して散乱パワーＰを求める。なお、ストークス行列を平均ストークス行列Ｆ’とするのは、コヒーレントな電波を用いて面状に分布する対象を計測した場合に生じるスペックルノイズを低減するためである。ここで、複数の対象ストークス行列Ｆを加算平均して散乱パワーＰを求めることにより、複数の対象からの散乱パワーＰを加算平均することができる。
なお、上述したような偏波特性および任意の送受信偏波の組による散乱パワーＰを求める方法は、例えば、上述した非特許文献３にも記載されている。
送受信偏波を変化させると、偏波によってターゲットの散乱強度が変化する。最適偏波組検出手段１０２は、このような偏波回転処理を利用して、予め定めておいた最適偏波組の各候補のうち最も強い反射強度が得られた候補の送受信偏波の組合せを、最適偏波組として検出してもよい。
最適偏波組の候補とする送受信偏波の組合せは、特に限定されない。すなわち、送信波および受信波に対して、それぞれ偏波の回転角が０〜１８０度の範囲に含まれる偏波を任意に設定できる。ここで、偏波の回転角は、水平偏波の偏波角度を０度としたときの当該偏波の偏波角度である。なお、偏波の回転角は、図４に示される角度ψに相当し、例えばＨ偏波であれば０度、Ｖ偏波であれば９０度となる。また、回転角の０度と１８０度とは同じ偏波を表わしている。最適偏波組の候補は、例えば、送信および受信の直線偏波を０度（水平）から１８０度まで＋５度ずつ回転させたものの組合せであってもよい。送信偏波の回転角をαｔとし、受信偏波の回転角をαｒとし、それらの組合せを｛αｔ，αｒ｝と表現すると、上記例の場合、最適偏波組の候補＝｛｛｛０，０｝，．．．，｛０，１７５｝｝，｛｛５，０｝，．．．，｛５，１７５｝｝，・・・，｛｛１７５，０｝，．．．，｛１７５，１７５｝｝｝となる。
また、最適偏波組検出手段１０２は、２組以上のＳＡＲ画像組を用いる場合には、各ＳＡＲ画像組を用いて各候補について反射強度を求め、反射強度の平均が最も高い候補や、反射強度の平均が所定値以上で反射強度のばらつきが最も小さい候補を最適偏波組に選んでもよい。
また、最適偏波組検出手段１０２は、全ての候補に対して反射強度を求めた結果、いずれの候補も反射強度が所定値以上とならなかった場合には、当該ピクセルは測定点には適さないとして、最適偏波組なしとしてもよい。
また、例えば、最適偏波組検出手段１０２は、最適偏波組を求める際に、まず変動が生じていないとされる期間中に撮影された２組以上のＳＡＲ画像組を用いて、ＰＳ特性を有するか否かを判定してもよい。最適偏波組検出手段１０２は、例えば、基本４偏波組のいずれかの送受信偏波の組による上記期間中の相対位相の変化幅が所定値以上であった場合には、ＰＳ特性を有していないとして、上記反射強度に基づく判定の如何にかかわらず、当該ピクセルは測定点には適さないとして、最適偏波組なしとしてもよい。
また、最適偏波組検出手段１０２は、最適偏波組の検出時間を短縮したい場合には、例えば、最適偏波組を、予め定めておいた最適偏波組の候補の中で所定の閾値以上の反射強度となる送受信偏波の組合せとしてもよい。そのような場合、最適偏波組検出手段１０２は、例えば、複数設けられた最適偏波組の候補に対して、所定の順番で偏波回転処理により反射強度を求め、最初に反射強度が所定の閾値以上となった候補を、当該ピクセルにおける最適偏波組としてもよい。以下では、最適偏波組における送受信偏波の回転角の組合せを、最適回転角組と称する場合がある。最適偏波組検出手段１０２による最適偏波組の決定は、最適回転角組を決定していると言うことも可能である。
最適偏波組検出手段１０２は、そのようにして得た各ピクセルにおける最適偏波組を示す情報を含む測定点情報を生成し、後述する測定点情報記憶手段１０３に記憶する。最適偏波組を示す情報は、例えば、送信偏波の回転角と、受信偏波の回転角とを含む情報であってもよい。なお、最適偏波組なしが検出された場合、最適偏波の組を示す情報には、その旨を示す情報が登録される。
測定点情報は、ＰＳ−ＩｎＳＡＲにおける測定点に関する情報であって、例えば、ピクセルまたはピクセルが対応する地点の位置情報と、最適偏波組を示す情報と対応づけた情報であってもよい。また、測定点情報は、位置情報と、最適偏波組を示す情報とに加えて、検出の際に求めた最適偏波組による反射強度（反射強度の平均を含む）や、偏波特性を示す情報（散乱行列Ｓやストークス行列Ｆ等）を含んでいてもよい。また、さらに、最適偏波組の検出日時や、検出に用いたＳＡＲ画像組の情報を含んでいてもよい。
測定点情報記憶手段１０３は、測定点情報を記憶する。
偏波回転処理手段１０４は、偏波回転処理を提供する手段であり、少なくとも１組のＳＡＲ画像組を用いて、指定されたピクセルにおける、指定された送受信偏波の組合せによる反射強度を算出する。
測定点情報出力手段１０７は、測定点情報記憶手段１０３に記憶されている測定点情報を出力する。測定点情報出力手段１０７は、例えば、測定点情報記憶手段１０３に記憶されている測定点情報を基に、対象ピクセルに対応する各地点が測定点に適するか否か、またはどの程度測定点に適するかを示す情報を作成し、出力してもよい。
例えば、測定点情報出力手段１０７は、入力されたＳＡＲ画像組と、測定点情報とを基に回転ＳＡＲ画像を生成し、生成した回転ＳＡＲ画像を、対象ピクセルに対応する各地点がどの程度測定点に適するかを示す情報として出力してもよい。
ここで、回転ＳＡＲ画像は、少なくとも対象とされたピクセルに対応づけて、当該ピクセルの最適偏波組によるＳＡＲデータである回転ＳＡＲデータを少なくとも保持する画像である。なお、回転ＳＡＲ画像は、通常のＳＡＲ画像と違って、ピクセルごとに送受信偏波の組合せが異なりうる画像である。各ピクセルが最適偏波組によるＳＡＲデータを保持する回転ＳＡＲ画像を出力することで、出力先で、各ピクセルの最適偏波組による反射強度画像を表示することができる。
なお、測定点情報出力手段１０７は、画像のデータ構造にとらわれずに、少なくとも対象ピクセルについて、当該ピクセルまたは当該ピクセルに対応する地点の位置情報と、当該ピクセルにおける最適偏波組を示す情報と、回転ＳＡＲデータとを含む情報を生成して、出力してもよい。そのような場合には、出力先で、必要な情報を選択して画像形式にして表示してもよい。
また、測定点情報出力手段１０７は、各ピクセルの偏波特性を示す情報（散乱行列Ｓやストークス行列Ｆ等）を含む情報を生成し、出力してもよい。これにより、出力先で任意の時刻における回転ＳＡＲデータを生成できる。
本実施形態において、ＳＡＲ画像組入力手段１０１は、例えば、データ入力装置と、プログラムに従って動作する情報処理装置とによって実現される。また、最適偏波組検出手段１０２および偏波回転処理手段１０４は、例えば、プログラムに従って動作する情報処理装置によって実現される。また、測定点情報記憶手段１０３は、例えば、記憶装置によって実現される。また、測定点情報出力手段１０７は、例えば、データ出力装置と、プログラムに従って動作する情報処理装置とによって実現される。
次に、本実施形態の動作を説明する。図５は、本実施形態の測定点情報提供装置１の動作の一例を示すフローチャートである。図５は、関心領域に対応するＳＡＲ画像の画像領域に含まれる全てのピクセルを対象に、最適偏波組を求めて、その結果得られる測定点情報を提供する測定点情報提供装置１の動作の例が示されている。
本例では、まず、ユーザに、関心領域を決定させる（ステップＳ１０１）。例えば、測定点情報提供装置１は、関心領域を指定するための関心領域指定画面等のユーザインタフェースを有し、該ユーザインタフェースを介して、ユーザに関心領域の情報を入力させてもよい。関心領域の情報は、例えば、中心地点の緯度経度および半径等であってもよい。また、例えば、合成開口レーダ３の撮影視野４をマッピングした地図等を表示して、その地図上で関心領域の範囲等を指定させてもよい。本例では、関心領域に対応するＳＡＲ画像の画像領域に含まれる全てのピクセルを対象ピクセルとする。
次に、ＳＡＲ画像組入力手段１０１は、関心領域をその撮影視野に含むＳＡＲ画像組を１組以上、入力する（ステップＳ１０２）。
次に、最適偏波組検出手段１０２は、入力されたＳＡＲ画像組を用いて偏波回転処理を行うことにより、対象ピクセルの最適偏波組を検出する（ステップＳ１０３）。また、最適偏波組検出手段１０２は、検出された対象ピクセルの最適偏波組を示す情報を含む測定点情報を生成し、測定点情報記憶手段１０３に記憶する。
最後に、測定点情報出力手段１０７は、測定点情報記憶手段１０３に記憶された測定点情報を出力する（ステップＳ１０４）。測定点情報出力手段１０７は、例えば、ユーザからの要求にしたがって、測定点情報を出力してもよい。また、測定点情報出力手段１０７は、測定点情報に代えてまたは測定点情報に加えて、測定点情報に基づいて生成した回転ＳＡＲ画像を出力してもよい。
以上のように、本実施形態によれば、ＰＳ−ＩｎＳＡＲを行う前に、変位を測定可能な測定点群を検出して、その情報をユーザに提供できるので、ユーザは、特定の送受信偏波の組合せによるＳＡＲ画像を見ただけでは見つけられないような測定点を容易に知ることができる。また、例えば、ユーザがその情報を基に測定点を指定したり、ＰＳ−ＩｎＳＡＲを行って得た変位の情報と照らし合わせたりすれば、高感度かつより正確な検出結果が得られる。
なお、上記説明では、測定点情報提供装置１が、偏波回転処理手段１０４を備える例を示しているが、最適偏波組検出手段１０２における最適偏波組の候補を、ＨＨ偏波、ＨＶ偏波、ＶＨ偏波およびＶＶ偏波の基本４偏波組のみとする場合には、偏波回転処理手段１０４を省略してもよい。
そのような場合には、最適偏波組検出手段１０２は、入力された少なくとも１組のＳＡＲ画像組を用いて、ピクセルごとに、基本４偏波組の中から最適偏波組を検出（選択）すればよい。最適偏波組検出手段１０２は、例えば、１つのＳＡＲ画像組を用いて、当該ＳＡＲ画像組に含まれる対象ピクセルの基本４偏波組によるＳＡＲデータの各々により示される反射強度を比較して、基本４偏波組の中から最も反射強度の強い偏波組を選び、選んだ偏波組を最適偏波組としてもよい。なお、複数のＳＡＲ画像組を用いる場合も同様に、ＳＡＲ画像組に含まれる対象ピクセルの基本４偏波組によるＳＡＲデータの各々により示される反射強度の平均が最も高い偏波組や、反射強度の平均が所定値以上で反射強度のばらつきが最も小さい偏波組を選んでもよい。
実施形態２．
次に、本発明の第２の実施形態について説明する。図６は、本発明の第２の実施形態の測定点情報提供装置１の構成例を示すブロック図である。図６に示す測定点情報提供装置１は、図２に示す第１の実施形態の測定点情報提供装置１と比べて、測定点情報出力手段１０７の代わりに、測定点候補入力手段１０８と、測定点条件判定手段１０９と、判定結果出力手段１１０とを備えている点が異なる。
測定点候補入力手段１０８は、関心領域内においてユーザが特に変位を測定したい地点（測定点）の候補を１つ以上入力する。測定点候補入力手段１０８は、例えば、測定点の候補を指定するための測定点候補指定画面等のユーザインタフェースを有し、該ユーザインタフェースを介して、ユーザに測定点の候補の情報を入力させてもよい。測定点の候補の情報は、例えば、地表における位置を示す情報（緯度経度等）であってもよいし、取得したＳＡＲ画像組に含まれるＳＡＲ画像におけるピクセル座標であってもよい。なお、本例の測定点情報提供装置１は、地表における位置とＳＡＲ画像におけるピクセル位置との対応関係を示す情報を保持しているものとする。
本実施形態では、指定された測定点の候補に対応するピクセルを対象ピクセルとする。
測定点条件判定手段１０９は、入力された測定点の候補が、測定点として適するか否かを判定する。測定点条件判定手段１０９は、指定された測定点の候補とされた地点（以下、候補点という）または当該候補点に対応するピクセルに対応づけられている測定点情報を参照し、当該候補点に対応するピクセルである対象ピクセルの最適偏波組によるＳＡＲデータ（反射強度や位相の情報）を基準にして、当該候補点において変位を測定できるか否かを判定する。なお、測定点条件判定手段１０９は、具体的には、測定点情報に基づいて、当該候補点が測定点として好適とされる所定の条件を満たしているか否かを判定すればよい。なお、当該候補点に対応する測定点情報が測定点情報記憶手段１０３に登録されていない場合、測定点条件判定手段１０９は、最新のＳＡＲ画像組を用いて対象ピクセルについて最適偏波組を検出する処理をするよう、最適偏波組検出手段１０２に要求してもよい。
測定点条件判定手段１０９は、例えば、測定点情報によって示される当該候補点における最適偏波組を示す情報に、最適偏波組なしの旨が登録されている場合には、当該候補点は測定点に適さないと判定してもよい。
また、例えば、測定点条件判定手段１０９は、測定点情報によって示される当該候補点における最適偏波組による反射強度が所定の値以上であれば、当該候補点は測定点として適すると判定してもよい。
判定結果出力手段１１０は、測定点条件判定手段１０９による判定結果を出力する。
本実施形態において、測定点候補入力手段１０８は、例えば、データ入力装置と、プログラムに従って動作する情報処理装置とによって実現される。また、測定点条件判定手段１０９は、プログラムに従って動作する情報処理装置によって実現される。また、判定結果出力手段１１０は、例えば、データ出力装置と、プログラムに従って動作する情報処理装置とによって実現される。
次に、本実施形態の動作を説明する。図７は、本実施形態の測定点情報提供装置１の動作の一例を示すフローチャートである。なお、ステップＳ１０１〜Ｓ１０２の処理は、図５に示した例と同じであるため、説明を省略する。
本例では、ＳＡＲ画像組入力手段１０１が関心領域をその撮影視野に含むＳＡＲ画像組を１組以上入力すると（ステップＳ１０２）、測定点候補入力手段１０８が、ユーザに、測定点の候補を１つ以上、指定させる（ステップＳ２０１）。
次に、最適偏波組検出手段１０２は、入力された少なくとも１組のＳＡＲ画像組を用いて偏波回転処理を行うことにより、測定点の候補とされた地点（候補点）に対応する対象ピクセルの最適偏波組を検出する（ステップＳ２０２）。また、最適偏波組検出手段１０２は、検出された対象ピクセルの最適偏波組を示す情報を含む測定点情報を生成し、測定点情報記憶手段１０３に記憶する。
次に、測定点条件判定手段１０９は、測定点情報に基づいて、指定された候補点が、測定点の条件としてあらかじめ定めておいた条件を満たすか否かを判定する（ステップＳ２０３）。ここでは、測定点条件判定手段１０９は、当該候補点に対応する対象ピクセルの最適偏波組によるＳＡＲデータを基準に、当該候補点が測定点に適するか否かを判定する。
最後に、判定結果出力手段１１０が、測定点条件判定手段１０９による判定結果を出力する（ステップＳ２０４）。
以上のように、本実施形態によれば、ＰＳ−ＩｎＳＡＲを行う前に、ユーザが変位を測定したい地点が測定点に適するか否かを容易に知ることができるので、ユーザは、効率よく変動の検出処理を行うことができる。例えば、ユーザは、変動を検出したい個所がある程度定まっている場合に、その個所やその付近の地点を指定して測定点として適するか否かを判定させ、その結果に基づいて測定点を指定してもよい。このとき、当該地点における最適偏波組の情報に基づいて、当該地点に対応するＳＡＲデータとして、変動を検知したい期間の前後の時刻における最適偏波組によるＳＡＲデータに相当する情報を入力してＰＳ−ＩｎＳＡＲを行うことで、高精度に変動を検出できる。なお、変動を検出したい期間の前後の時刻における最適偏波組によるＳＡＲデータに相当する情報は、当該時刻におおける基本４偏波組のＳＡＲデータを用いて、上記の偏波回転処理等を行うことにより生成できる。
実施形態３．
次に、本発明の第３の実施形態について説明する。図８は、本実施形態の観測システムの構成例を示す構成図である。図８に示す観測システムは、図１に示す第１の実施形態と比べて、測定点情報提供装置１に代えて、変動検出装置５を備えている点が異なる。
また、図９は、変動検出装置５の構成例を示すブロック図である。図９に示すように、変動検出装置５は、図２に示した第１の実施形態の測定点情報提供装置１の構成要素（ＳＡＲ画像組入力手段１０１、最適偏波組検出手段１０２、測定点情報記憶手段１０３、偏波回転処理手段１０４および測定点情報出力手段１０７）に加えて、さらに変動検出手段１０５と、検出結果出力手段１０６とを備えている。
本実施形態において、ＳＡＲ画像組入力手段１０１は、少なくとも観測時間の異なる２組以上のＳＡＲ画像組を入力する。ここで入力されるＳＡＲ画像組には、ユーザが指定した、変動を検出したい時間区間の前後の時刻に対応するＳＡＲ画像組が少なくとも含まれる。以下、ユーザが指定した変動を検出したい時間区間の前後の時刻を指定時刻という場合がある。なお、時間区間は１つとは限らない。また、ＳＡＲ画像組入力手段１０１が入力する２組以上のＳＡＲ画像組を、略同一場所に対して、観測時刻の異なる観測データから生成されるＳＡＲ画像組の集合という意味で「ＳＡＲ画像組系列」という場合がある。
ＳＡＲ画像組入力手段１０１は、例えば、図３に示されるような、観測時刻の異なるｎ組（例えば、数十組）のＳＡＲ画像組からなるＳＡＲ画像組系列を入力してもよい。
変動検出手段１０５は、測定点情報記憶手段１０３に記憶されている測定点情報と、入力されたＳＡＲ画像組系列とに基づいて、少なくとも対象ピクセルに対応する地点における変動を検出する。ここで、対象ピクセルは、例えば、関心領域に対応する画像領域に含まれる各ピクセルまたは関心領域内において特に変位を測定したい地点（測定点）として指定されたピクセルである。
検出結果出力手段１０６は、変動検出手段１０５による検出結果を出力する。
なお、最適偏波組検出手段１０２、測定点情報記憶手段１０３、偏波回転処理手段１０４および測定点情報出力手段１０７は、第１の実施形態と同様でよい。
以下、本実施形態の変動検出手段１０５についてさらに詳しく説明する。図１０は、変動検出手段１０５の一例を示すブロック図である。図１０に示す変動検出手段１０５Ａは、回転ＳＡＲデータ生成手段１５１と、変位測定手段１５２Ａとを含む。
回転ＳＡＲデータ生成手段１５１は、測定点情報記憶手段１０３に記憶されている測定点情報によって示される対象ピクセルの最適偏波組と、入力されたＳＡＲ画像組系列に含まれる指定時刻のＳＡＲ画像組とに基づいて、対象ピクセルの各々に対して、少なくとも２つの異なる時刻に対応する回転ＳＡＲデータを生成する。
なお、回転ＳＡＲデータは、上述したように、ある時刻（ここでは指定時刻に含まれる各時刻）における当該ピクセルの最適偏波組によるＳＡＲデータである。
回転ＳＡＲデータ生成手段１５１は、例えば、偏波回転処理手段１０４が提供する偏波回転処理を利用して、対象ピクセルについて、指定時刻（より具体的には、指定時刻に含まれる各時刻）に対応する回転ＳＡＲデータを生成してもよい。以下、回転ＳＡＲデータ生成手段１５１が生成する対象ピクセルの指定時刻に対応する回転ＳＡＲデータを、対象ピクセルについての観測時刻の異なる観測データから生成される回転ＳＡＲデータの集合という意味で、対象ピクセルの「回転ＳＡＲデータ系列」という場合がある。
回転ＳＡＲデータ生成手段１５１は、例えば、対象ピクセルの各々に対して、ＳＡＲ画像組系列に含まれる指定時刻のＳＡＲ画像組を用いて、指定時刻に含まれる時刻ごとに、最適偏波組を指定した偏波回転処理を行うことにより、各対象ピクセルについての回転ＳＡＲデータ系列を生成できる。
なお、回転ＳＡＲデータは、最適偏波組が基本４偏波組のいずれかである場合には、入力されたＳＡＲ画像組に含まれるＳＡＲデータをそのまま利用してもよい。
変位測定手段１５２Ａは、対象ピクセルについての回転ＳＡＲデータ系列を用いて、対象ピクセルに対応する地点の変位を測定する。変位測定手段１５２Ａは、より具体的には、対象ピクセルについての回転ＳＡＲデータ系列をＰＳ−ＩｎＳＡＲ処理に用いて、対象ピクセルに対応する地点における、変動による位相差を算出する。ＰＳ−ＩｎＳＡＲ処理では、比較対象とされた対象ピクセルのＳＡＲデータ間に現れる位相差のうち、高さによる位相差を除去して、変動による位相差を求めることにより、対象ピクセルに対応する地点における変位を測定できる。
例えば、ｎ組のＳＡＲ画像組を含むＳＡＲ画像組系列が入力され、そのｎ組のＳＡＲ画像組が取得された各時刻が指定時刻として指定された場合、変位測定手段１５２Ａには、１つの対象ピクセルにつき、ｔ１〜ｔｎの時刻にそれぞれ対応するｎ個の回転ＳＡＲデータを含む回転ＳＡＲデータ系列が入力される。
そのような場合に、例えば、変位測定手段１５２Ａは、最初の時刻とそれ以降の各時刻間に現れる位相差の最大値や、隣り合った時刻間に現れる位相差の累計値等を、対象ピクセルの指定時刻間における変位を示す情報として出力してもよい。具体的には、時刻ｔ１に対応する回転ＳＡＲデータと、時刻ｔ２〜ｔｎの各々に対応する回転ＳＡＲデータ間に現れる位相差を求め、求めた位相差の最大値を、対象ピクセルの指定時刻間における変位を示す情報として出力してもよい。また、例えば、変位測定手段１５２Ａは、時刻ｔ１に対応する回転ＳＡＲデータと時刻ｔ２に対応する回転ＳＡＲデータ間に現れる位相差，時刻ｔ２に対応する回転ＳＡＲデータと時刻ｔ３に対応する回転ＳＡＲデータ間に現れる位相差，・・・，時刻ｔ（ｎ−１）に対応する回転ＳＡＲデータと時刻ｔｎに対応する回転ＳＡＲデータ間に現れる位相差を求め、求めた位相差の累計値（ただし、２πを上限とする）を、対象ピクセルの指定時刻間における変位を示す情報として出力してもよい。
このような場合、検出結果出力手段１０６は、変位測定手段１５２Ａから出力される各対象ピクセルの指定時刻間における変位を示す情報を、各対象ピクセルに対応する地点における変動の検出結果を示す情報として出力してもよい。
また、検出結果出力手段１０６は、変位測定手段１５２Ａから出力される各対象ピクセルの指定時刻間における変位を示す情報を基に、差分干渉画像を生成して表示してもよい。差分干渉画像は、少なくとも対象ピクセルの各々に対応づけて、当該ピクセルが対応する地点における変動による位相差の情報または該位相差の大きさに応じた画素情報を保持する画像である。これにより、出力先で、差分干渉画像の表示において、各画素を、当該画素に対応する地点における変位量に応じて色分け等できる。検出結果出力手段１０６は、生成する差分干渉画像において、対象ピクセル以外のピクセルについては、測定対象外もしくは変位なしの旨の情報を保持させてもよい。
また、図１１は、変動検出手段１０５の他の例を示すブロック図である。図１１に示す変動検出手段１０５Ｂは、回転ＳＡＲデータ生成手段１５１と、回転ＳＡＲ画像生成手段１５３Ｂと、変位測定手段１５２Ｂとを含む。
回転ＳＡＲデータ生成手段１５１は、図１０に示した変動検出手段１０５Ａの回転ＳＡＲデータ生成手段１５１と同様である。
回転ＳＡＲ画像生成手段１５３Ｂは、生成された各対象ピクセルについての回転ＳＡＲデータ系列を統合して、指定時刻の各時刻に対応する回転ＳＡＲ画像を生成する。
なお、回転ＳＡＲ画像生成手段１５３Ｂは、各回転ＳＡＲ画像において、対象ピクセル以外のピクセルについては、無効の旨や反射強度をゼロとした回転ＳＡＲデータを保持させてもよいし、例えば、ＳＡＲ画像組系列に含まれる当該時刻の基本４偏波組のいずれか１つの送受信偏波の組によるＳＡＲデータを選び、それを保持させてもよい。
変位測定手段１５２Ｂは、生成された複数時刻に対応する回転ＳＡＲ画像を用いて、少なくとも対象ピクセルに対応する地点における変位を測定し、測定結果を表す差分干渉画像を生成する。変位測定手段１５２Ｂは、例えば、複数時刻に対応するＳＡＲ画像を用いたＰＳ−ＩｎＳＡＲ処理と同様の処理を行えばよい。
なお、上述した変位測定手段１５２Ａでは、入力データが、ピクセル単位のデータである回転ＳＡＲデータ系列であり、出力データが、ピクセル単位の測定結果を示す情報であったのに対して、変位測定手段１５２Ｂでは、入力データが、画像単位のデータである回転ＳＡＲ画像であり、出力データが、画像単位の測定結果である差分干渉画像である点が異なる。このように、変位測定手段への入出力データのフォーマットを選べることができる。
本実施形態において、変動検出手段１０５は、例えば、プログラムに従って動作する情報処理装置によって実現される。また、検出結果出力手段１０６は、例えば、データ出力装置と、プログラムに従って動作する情報処理装置とによって実現される。
次に、本実施形態の動作を説明する。図１２は、本実施形態の変動検出装置５の動作の一例を示すフローチャートである。図１２には、関心領域に対応するＳＡＲ画像の画像領域に含まれる全てのピクセルについて変位測定を行って対象物の歪み等を検出する変動検出装置５の動作の例が示されている。本例では、変動検出装置５が、変動検出手段１０５として、図１０に示される変動検出手段１０５Ａを備えている場合を例に用いて説明する。なお、以下では、ステップＳ１０１は、図５および図７に示した例と同じであるため、説明を省略する。
また、本例では、関心領域に対応するＳＡＲ画像の画像領域に含まれる全てのピクセルを対象ピクセルとする。
本例では、ユーザに、関心領域を決定させると（ステップＳ１０１）、次に、ＳＡＲ画像組入力手段１０１が、関心領域を含む領域に対して異なる時刻に取得されたＳＡＲ画像組を含むＳＡＲ画像組系列を入力する（ステップＳ３０１）。
次に、最適偏波組検出手段１０２は、入力されたＳＡＲ画像組系列に含まれる少なくとも１組のＳＡＲ画像組を用いて、対象ピクセルの各々について偏波回転処理を行うことにより、各対象ピクセルの最適偏波組を検出する（ステップＳ３０２）。また、最適偏波組検出手段１０２は、得られた各対象ピクセルの最適偏波組を示す情報を含む測定点情報を生成し、測定点情報記憶手段１０３に記憶する。
次に、測定点情報出力手段１０７は、測定点情報記憶手段１０３に記憶された測定点情報を出力する（ステップＳ１０４）。測定点情報出力手段１０７は、例えば、ユーザからの要求にしたがって、測定点情報を出力してもよい。また、測定点情報出力手段１０７は、測定点情報に代えてまたは測定点情報に加えて、測定点情報に基づいて生成した回転ＳＡＲ画像を出力してもよい。
次に、変動検出手段１０５の回転ＳＡＲデータ生成手段１５１は、測定点情報によって示される各対象ピクセルの最適偏波組と、ＳＡＲ画像組系列に含まれる指定時刻に対応するＳＡＲ画像組とに基づいて、各対象ピクセルの回転ＳＡＲデータ系列を生成する（ステップＳ３０３）。
次に、変動検出手段１０５の変位測定手段１５２Ａは、生成された各対象ピクセルの回転ＳＡＲデータ系列を用いて、対象ピクセルごとにＰＳ−ＩｎＳＡＲ処理を行うことにより、各対象ピクセルに対応する地点における変位を測定する（ステップＳ３０４）。ここで、変位測定手段１５２Ａは、例えば、測定の結果得られた各対象ピクセルの指定時刻間における変位を示す情報を出力する。なお、対象ピクセルによっては反射強度が小さくて位相差の情報が得られないなど測定点の条件を満たさず、測定不能を示す情報が得られる場合がある。
最後に、検出結果出力手段１０６は、ステップＳ３０４で得られた情報を基に差分干渉画像を生成し、生成した差分干渉画像を、関心領域における変動の検出結果を示す情報として出力する（ステップＳ３０５）。
なお、ステップＳ１０４の処理は、ステップＳ３０３以降に行ってもよい。例えば、ステップＳ３０５による差分干渉画像の出力と併せて、測定点情報や回転ＳＡＲ画像を出力してもよい。
また、図１３は、本実施形態の変動検出装置５の動作の他の例を示すフローチャートである。図１３には、関心領域に対応するＳＡＲ画像の画像領域に含まれる全てのピクセルについて変位測定を行って対象物の歪み等を検出する変動検出装置５の動作の例が示されている。本例では、変動検出装置５が、変動検出手段１０５として、図１１に示される変動検出手段１０５Ｂを備えている場合を例に用いて説明する。なお、ステップＳ１０１〜Ｓ３０３までの処理は、図１２に示した例と同じであるため、説明を省略する。
また、本例でも、関心領域に対応するＳＡＲ画像の画像領域に含まれる全てのピクセルを対象ピクセルとする。
本例では、変動検出手段１０５の回転ＳＡＲデータ生成手段１５１が、各対象ピクセルの回転ＳＡＲデータ系列を生成すると（ステップＳ３０３）、回転ＳＡＲ画像生成手段１５３Ｂが、指定時刻に含まれる時刻ごとに、各対象ピクセルの回転ＳＡＲデータ系列に含まれる当該時刻に対応する回転ＳＡＲデータを統合して、指定時刻に対応する回転ＳＡＲ画像を生成する（ステップＳ４０１）。
次に、変動検出手段１０５の変位測定手段１５２Ｂは、生成された指定時刻に対応する回転ＳＡＲ画像を用いて、ＰＳ−ＩｎＳＡＲ処理を行うことにより、差分干渉画像を生成する（ステップＳ４０２）。
最後に、検出結果出力手段１０６は、少なくとも関心領域における変動の検出結果を示す情報として、ステップＳ４０２で得られた差分干渉画像を出力する（ステップＳ４０３）。
以上のように、本実施形態によれば、ユーザは、関心領域に含まれる各ピクセルが測定点に適するか否かを知ることができる上に、関心領域に含まれる全ての測定可能な地点を使って変動を検出できる。このことは、本発明によって変動検出の感度が向上したことに他ならない。また、変動検出結果を差分干渉画像として出力すれば、関心領域内で変動があった地点が即座にわかるなど検出結果の視認性を高めることができる。
また、図１２や図１３に示した方法のように、変動検出結果を差分干渉画像として出力すれば、関心領域内で変動があった地点が即座にわかるなど検出結果の視認性を高めることができる。また、測定点情報に代えてまたは測定点情報に加えて、回転ＳＡＲ画像を出力すれば、ユーザは、差分干渉画像と回転ＳＡＲ画像とを比較することにより、ＰＳ−ＩｎＳＡＲを行って得た検出結果をより正確に分析できるようになる。
また、図１３に示した方法によれば、複数の時刻に対応したＳＡＲ画像を入力して差分干渉画像を出力するような既存のＰＳ−ＩｎＳＡＲ処理を利用することができる。
実施形態４．
次に、本発明の第４の実施形態について説明する。図１４は、本実施形態の変動検出装置５の構成例を示すブロック図である。図１４に示す変動検出装置５は、図６に示した第２の実施形態の測定点情報提供装置１の構成要素（ＳＡＲ画像組入力手段１０１、最適偏波組検出手段１０２、測定点情報記憶手段１０３、偏波回転処理手段１０４、測定点候補入力手段１０８、測定点条件判定手段１０９および判定結果出力手段１１０）に加えて、変動検出手段１０５と、検出結果出力手段１０６とを備えている。
変動検出手段１０５および検出結果出力手段１０６は、基本的には第３の実施形態と同様である。なお、本実施形態における変動検出手段１０５は、測定点条件判定手段１０９によって測定点として適すると判定された候補点を測定点に用いて、当該測定点における変動を検出する。
本実施形態では、測定点が決定される前までは、指定された測定点の候補に対応するピクセルを対象ピクセルとし、測定点が決定された後は、決定された測定点に対応するピクセルを対象ピクセルとする。
本実施形態において、測定点条件判定手段１０９は、ユーザから指定された測定点の候補とされた地点である候補点が、測定点に適するか否かを判定し、測定点として適すると判定した場合に、当該候補点を測定点に設定する。また、測定点条件判定手段１０９は、候補点が測定点に適さないと判定した場合に、判定結果出力手段１１０にその旨を出力させる。なお、測定点条件判定手段１０９は、候補点が測定点に適するか否かにかかわらず、判定結果出力手段１１０に判定結果を出力させて、当該候補点を測定点として用いるか否かをユーザに確認してもよい。
図１５は、本実施形態の変動検出装置５の動作の一例を示すフローチャートである。図１５には、関心領域内のユーザが指定した地点について変位測定を行って対象物の歪み等を検出する変動検出装置５の動作の例が示されている。なお、以下では、変動検出装置５が、変動検出手段１０５として、図１０に示される変動検出手段１０５Ａを備えている場合を例に用いて説明する。なお、以下では、ステップＳ１０１およびＳ３０１の処理は、図１２および図１３に示した例と同じであるため、説明を省略する。
本例では、ＳＡＲ画像組入力手段１０１が関心領域を含む領域に対して異なる時刻に取得されたＳＡＲ画像組を含むＳＡＲ画像組系列を入力すると（ステップＳ３０１）、次に、測定点候補入力手段１０８が、ユーザに、関心領域内において特に変位を測定したい地点（測定点）の候補を１つ以上、指定させる（ステップＳ２０１）。
次に、最適偏波組検出手段１０２は、入力された少なくとも１組のＳＡＲ画像組を用いて偏波回転処理を行うことにより、指定された候補点に対応するピクセルの最適偏波組を検出する（ステップＳ２０２）。また、最適偏波組検出手段１０２は、検出された当該ピクセルの最適偏波組を示す情報を含む測定点情報を生成し、測定点情報記憶手段１０３に記憶する。
次に、測定点条件判定手段１０９は、測定点情報に基づいて、指定された候補点が、測定点の条件としてあらかじめ定めておいた条件を満たすか否かを判定する（ステップＳ２０３）。なお、ステップＳ２０１からＳ２０３までの処理は、図７に示した第２の実施形態における動作例と同様である。
ステップ２０３の判定の結果、指定された候補点が測定点の条件を満たさない場合には、判定結果出力手段１１０を介してユーザにその旨を出力する（ステップＳ５０１のＮｏ，ステップＳ５０２）。なお、この後は、他の測定点の候補を指定させるため、ステップＳ２０１に戻ってもよい。
一方、指定された候補点が測定点の条件を満たす場合には、当該候補点を測定点に設定する（ステップＳ５０３）。以下、本例において、対象ピクセルといった場合は、当該測定点に対応するピクセルを指す。
次に、変動検出手段１０５の回転ＳＡＲデータ生成手段１５１は、測定点情報によって示される対象ピクセルの最適偏波組と、ＳＡＲ画像組系列に含まれる指定時刻に対応するＳＡＲ画像組とに基づいて、対象ピクセルについての回転ＳＡＲデータ系列を生成する（ステップＳ５０４）。
次に、変動検出手段１０５の変位測定手段１５２Ａは、生成された対象ピクセルの回転ＳＡＲデータ系列を用いてＰＳ−ＩｎＳＡＲ処理を行うことにより、対象ピクセルに対応する地点（測定点の候補として指定された地点）における変位を測定する（ステップＳ５０５）。そして、変位測定手段１５２Ａは、例えば、測定の結果得られた対象ピクセルの指定時刻間における変位を示す情報（指定時刻に含まれる時刻間における位相差の最大値や累計値等）を出力する。なお、対象ピクセルによっては最適偏波組であっても反射強度が小さい、またはＰＳ特性を有していないなど測定点の条件を満たさず、測定不能を示す情報が出力される場合がある。
最後に、検出結果出力手段１０６は、指定された地点における変動の検出結果を示す情報として、ステップＳ５０５で得られた情報を出力する（ステップＳ５０６）。検出結果出力手段１０６は、例えば、変位測定手段１５２Ａから得られた対象ピクセルの指定時刻間における変位を示す情報に基づいて、指定された地点における変動の有無を判定して、その結果（変動有無）とともにステップＳ１０６で得られた情報を出力してもよい。
測定点の候補が複数ある場合は、ステップＳ２０２からステップＳ５０６またはＳ５０２までの処理を繰り返し行う。
以上のように、本実施形態によれば、ユーザが指定した地点における変動を検出できる可能性を高めることができる。なぜなら、変動検出装置５が当該地点に対応するピクセルの最適偏波組を求めた上で、複数時刻における当該地点の最適偏波組によるＳＡＲデータである回転ＳＡＲデータ系列を生成して、ＰＳ−ＩｎＳＡＲ処理を行うからである。例えば、１つの送受信偏波の組合せによるＳＡＲデータを用いるのみの場合、当該地点は該送受信偏波の組による反射強度が小さく測定点として処理できない可能性がある。それに対して、本発明によれば、そのような地点であっても、他に反射強度の大きい送受信偏波の組合せがあれば、その送受信偏波の組合せを用いて当該地点における変位を測定できる。このことは、本発明によって変動検出の感度が向上したことに他ならない。
また、さらに本実施形態によれば、最適偏波組によっても当該地点を測定点として処理できない場合には、その旨をユーザに出力することにより、他の付近の地点を代替させるといったことができる。このことは、本発明によって変動検出の感度がさらに向上したことに他ならない。
また、ユーザは、測定点の指定が容易かつより確実にできるので、対象とする物体または地表の変動をより精密に検出できる。
また、このようにして測定点の指定が容易かつより確実になることによって、ビル、道路、橋等の人口構造物の劣化診断や、地盤沈下測定、地殻変動測定等の土地起伏の変化診断などがより精密に行えるという効果が得られる。
なお、上記第３および第４の実施形態でも、最適偏波組の候補を任意に設定できる場合を例に用いて説明したが、例えば、最適偏波組の候補を、ＨＨ偏波、ＨＶ偏波、ＶＨ偏波およびＶＶ偏波の基本４偏波組のみとしてもよい。そのような場合には、変動検出装置５は、偏波回転処理手段１０４を省略できる。
なお、その場合において、変動検出手段１０５は、測定点情報に基づき、必要に応じてＳＡＲ画像組系列に含まれる指定時刻における対象ピクセルの基本４偏波組によるＳＡＲデータの中から該対象ピクセルの最適偏波組によるＳＡＲデータを選択的に用いて、対象ピクセルに対応する地点における変動を検出すればよい。
図１６および図１７は、最適偏波組の候補を、ＨＨ偏波、ＨＶ偏波、ＶＨ偏波およびＶＶ偏波の基本４偏波組のみとした場合の変動検出手段１０５の構成例を示すブロック図である。図１６に示す変動検出手段１０５Ｃは、変位測定手段１５２Ｃを含む。
変位測定手段１５２Ｃは、図１０に例示した変動検出手段１０５Ａの変位測定手段１５２Ａと基本的に同じであるが、対象ピクセルについての回転ＳＡＲデータ系列を用いる代わりに、ＳＡＲ画像組系列に含まれる指定時刻における対象ピクセルの基本４偏波組によるＳＡＲデータを用いる点が異なる。変位測定手段１５２Ｃは、例えば、ＳＡＲ画像組系列に含まれる指定時刻における対象ピクセルの基本４偏波組によるＳＡＲデータが入力されると、測定点情報に基づいて、各時刻につき最適偏波組によるＳＡＲデータを選択して、対象ピクセルに対応する地点における変位を測定する。
なお、変位測定手段１５２Ｃは、測定点情報に基づき、ＳＡＲ画像組系列に含まれる指定時刻における当該ピクセルの最適偏波組によるＳＡＲデータを選択して用いる点以外は、上述した変位測定手段１５２Ａと同様である。
また、図１７に示す変動検出手段１０５Ｄは、回転ＳＡＲ画像生成手段１５３Ｄと、変位測定手段１５２Ｂとを含む。
回転ＳＡＲ画像生成手段１５３Ｄは、図１１に例示した第１の実施形態の変動検出手段１０５の回転ＳＡＲ画像生成手段１５３Ｂと基本的に同じであるが、対象ピクセルについての回転ＳＡＲデータ系列を用いる代わりに、ＳＡＲ画像組系列に含まれる指定時刻における対象ピクセルの基本４偏波組によるＳＡＲデータを用いる点が異なる。回転ＳＡＲ画像生成手段１５３Ｄは、例えば、ＳＡＲ画像組系列に含まれる指定時刻における対象ピクセルの基本４偏波組によるＳＡＲデータが入力されると、時刻ごとに、測定点情報に基づいて、対象ピクセルの各々に対して最適偏波組によるＳＡＲデータを選択し、それらを統合して、指定時刻に対応する回転ＳＡＲ画像を生成する。
なお、回転ＳＡＲ画像生成手段１５３Ｄは、測定点情報に基づき選択したＳＡＲ画像組系列に含まれる指定時刻における当該ピクセルの最適偏波組によるＳＡＲデータを用いる点以外は、上述した回転ＳＡＲ画像生成手段１５３Ｂと同様である。
また、上記各実施形態では、測定点情報を記憶する手段として、測定点情報提供装置１および変動検出装置５が、測定点情報記憶手段１０３を備える例を示したが、測定点情報記憶手段１０３を省略することも可能である。そのような場合には、最適偏波組検出手段１０２から直接、測定点情報出力手段１０７や測定点条件判定手段１０９や変動検出手段１０５に、測定点情報を出力すればよい。
次に、本発明の概要について説明する。図１８は、本発明の概要を示すブロック図である。図１８に示すように、本発明による測定点情報提供装置または変動検出装置は、ＳＡＲ画像組入力手段６０１と、偏波組決定手段６０２と、測定点情報生成手段６０３とを備えている。
ＳＡＲ画像組入力手段６０１（例えば、ＳＡＲ画像組入力手段１０１）は、特定領域（関心領域）を含む撮影視野内の分解能セルに対応するピクセルの各々に対応づけて、略同一時刻に合成開口レーダが観測した送受信偏波の組合せがＨＨ偏波、ＨＶ偏波、ＶＨ偏波およびＶＶ偏波である基本４偏波組による観測データから生成される反射強度および位相を示す情報を少なくとも保持するＳＡＲ画像組を入力する。
偏波組決定手段６０２（例えば、最適偏波組検出手段１０２）は、入力されたＳＡＲ画像組を用いて、対象ピクセルの各々に対して、反射強度が所定値以上または最も強い偏波組を決定する。
測定点情報生成手段６０３（例えば、最適偏波組検出手段１０２）は、偏波組決定手段６０２により決定された偏波組に基づいて、少なくとも対象ピクセルの決定された偏波組を示す情報を含む測定点情報を生成する。
このような特徴的要素を備えているので、変位を測定可能な地点をより多くまたは見逃さずに検出でき、したがって、その情報をユーザに提供したり、測定点の指定や測定結果の分析に利用することにより、変動検出の感度を向上させることができる。
また、図示省略しているが、本発明による測定点情報提供装置は、さらに、測定点情報生成手段６０２が生成した測定点情報、または測定点情報に基づいて生成される、対象ピクセルに対応する地点が測定点に適するか否かもしくはどの程度測定点に適するかを示す情報を出力する測定点情報出力手段（例えば、測定点情報出力手段１０７）を備えていてもよい。
また、測定点情報出力手段は、入力されたＳＡＲ画像組と、測定点情報とに基づいて、対象ピクセルに対応する地点がどの程度測定点に適するかを示す情報として、少なくとも対象ピクセルに対応づけて、当該ピクセルの決定された偏波組による反射強度および位相を示す情報を少なくとも保持する回転ＳＡＲ画像を生成して、出力してもよい。
また、図示省略しているが、本発明による測定点情報提供装置は、さらに、測定点の候補を入力する測定点候補入力手段（例えば、測定点候補入力手段１０８）と、測定点情報に基づいて、入力された測定点の候補が測定点として適するか否かを、当該候補に対応するピクセルの決定された偏波組による反射強度を少なくとも用いて判定する測定点条件判定手段（例えば、測定点条件判定手段１０９）と、測定点条件判定手段による判定結果を出力する判定結果出力手段（例えば、判定結果出力手段１１０）とを備え、偏波組決定手段６０２は、入力された測定点の候補に対応するピクセルである対象ピクセルの各々に対して、反射強度が所定値以上または最も強い偏波組を決定し、測定点情報生成手段６０３は、偏波組決定手段６０２による偏波組の決定結果に基づき、少なくとも対象ピクセルの決定された偏波組を示す情報と、当該対象ピクセルの該偏波組による反射強度を示す情報とを含む測定点情報を生成してもよい。
以上、実施形態及び実施例を参照して本願発明を説明したが、本願発明は上記実施形態に限定されるものではない。本願発明の構成や詳細には、本願発明のスコープ内で当業者が理解し得る様々な変更をすることができる。
また、上記の実施形態の一部または全部は、以下の付記のようにも記載されうるが、以下には限られない。
（付記１）特定領域を含む撮影視野内の分解能セルに対応するピクセルの各々に対応づけて、略同一時刻に合成開口レーダが観測した送受信偏波の組合せがＨＨ偏波、ＨＶ偏波、ＶＨ偏波およびＶＶ偏波である基本４偏波組による観測データから生成される反射強度および位相を示す情報を少なくとも保持するＳＡＲ画像組を２組以上入力するＳＡＲ画像組入力手段（例えば、ＳＡＲ画像組入力手段１０１）と、入力されたＳＡＲ画像組のうちの少なくとも１組を用いて、対象ピクセルの各々に対して、反射強度が所定値以上または最も強い偏波組を決定する偏波組決定手段（例えば、最適偏波組決定手段１０２）と、偏波組決定手段により決定された偏波組に基づいて、少なくとも対象ピクセルの決定された偏波組を示す情報を含む測定点情報を生成する測定点情報生成手段（例えば、最適偏波組決定手段１０２）と、測定点情報と、入力されたＳＡＲ画像組とに基づいて、対象ピクセルに対応する地点における変位を測定する変位測定手段（例えば、変動検出手段１０５、特に変位測定手段１５２Ａ、変位測定手段１５２Ｂ、変位測定手段１５２Ｃ）と、変位測定手段による測定結果を示す情報を出力する第１の出力手段（例えば、検出結果出力手段１０６）と、測定点情報、または測定点情報に基づいて生成される、対象ピクセルに対応する地点が測定点に適するか否かもしくはどの程度測定点に適するかを示す情報を出力する第２の出力手段（例えば、測定点情報出力手段１０７）とを備えたことを特徴とする変動検出装置。
（付記２）関心領域を含む撮影視野内の分解能セルに対応するピクセルの各々に対応づけて、略同一時刻に合成開口レーダが観測した送受信偏波の組合せがＨＨ偏波、ＨＶ偏波、ＶＨ偏波およびＶＶ偏波である基本４偏波組による観測データから生成される反射強度および位相を示す情報を少なくとも保持するＳＡＲ画像組を２組以上入力するＳＡＲ画像組入力手段（例えば、ＳＡＲ画像組入力手段１０１）と、測定点の候補を入力する測定点候補入力手段（例えば、測定点条件判定手段１０８）と、入力されたＳＡＲ画像組のうちの少なくとも１組を用いて、入力された測定点の候補に対応するピクセルである対象ピクセルの各々に対して、反射強度が所定値以上または最も強い偏波組を決定する偏波組決定手段（例えば、最適偏波組決定手段１０２）と、偏波組決定手段により決定された偏波組に基づいて、少なくとも対象ピクセルの決定された偏波組を示す情報と、当該対象ピクセルの決定された偏波組による反射強度を示す情報とを含む測定点情報を生成する測定点情報生成手段（例えば、最適偏波組決定手段１０２）と、測定点情報に基づいて、入力された測定点の候補が測定点として適するか否かを、当該候補に対応するピクセルの決定された偏波組による反射強度を少なくとも用いて判定する測定点条件判定手段（例えば、測定点条件判定手段１０９）と、測定点条件判定手段による判定結果を出力する判定結果出力手段（例えば、判定結果出力手段１１０）と、測定点条件判定手段によって測定点として適すると判定された候補を測定点に用いて、測定点情報と、入力されたＳＡＲ画像組とに基づいて、測定点における変位を測定する変位測定手段（例えば、変動検出手段１０５、特に変位測定手段１５２Ａ、変位測定手段１５２Ｂ、変位測定手段１５２Ｃ）とを備えたことを特徴とする変動検出装置。
（付記３）偏波組決定手段は、対象ピクセルの各々に対して、予め定められている偏波組の候補の中から、反射強度が所定値以上または最も強い偏波組を決定する付記１または付記２に記載の変動検出装置。
（付記４）偏波組の候補が、ＨＨ偏波、ＨＶ偏波、ＶＨ偏波およびＶＶ偏波の４種類である付記１から付記３のいずれかに記載の変動検出装置。
（付記５）偏波組の候補が、偏波回転角が０度〜１８０度のいずれかである送信偏波と、偏波回転角が０度〜１８０度のいずれかである受信偏波との組合せの集合である付記１から付記３のいずれかに記載の変動検出装置。
（付記６）指定された点を画像領域に含む少なくともＳＡＲ画像組を１組以上用いて、当該点における任意の送受信偏波の組合せによる反射強度を算出する偏波回転処理手段と、対象ピクセルの各々に対して、偏波回転処理手段を利用して、入力されたＳＡＲ画像組から、少なくとも２つの時刻に対応する回転ＳＡＲデータであって、それぞれ対応する時刻における当該ピクセルの決定された偏波組による反射強度および位相を示す情報を含む回転ＳＡＲデータを生成する回転ＳＡＲデータ生成手段とを備え、偏波組決定手段は、対象ピクセルの各々に対して、偏波回転処理手段を利用して、入力されたＳＡＲ画像のうちの少なくとも１組のＳＡＲ画像組から、予め定められている偏波組の各候補の反射強度を算出し、算出した各候補の反射強度に基づいて偏波組を決定し、変位測定手段は、対象ピクセルの各々に対して、回転ＳＡＲデータ生成手段によって生成された当該ピクセルの少なくとも２つの時刻に対応する回転ＳＡＲデータに基づいて、当該ピクセルに対応する地点における変位を測定する付記１から付記５のいずれかに記載の変動検出装置。
（付記７）偏波組決定手段によって決定された偏波組と、入力されたＳＡＲ画像組とに基づいて、少なくとも２つの時刻に対応する回転ＳＡＲ画像であって、それぞれが少なくとも対象ピクセルの各々に対応づけて、対応する時刻における当該ピクセルの決定された偏波組による反射強度および位相を示す情報を少なくとも保持する回転ＳＡＲ画像を生成する回転ＳＡＲ画像生成手段を備え、変位測定手段は、回転ＳＡＲ画像生成手段によって生成された少なくとも２つの時刻に対応する回転ＳＡＲ画像を用いて、対象ピクセルの各々に対して、当該ピクセルに対応する地点における変位を測定する付記６に記載の変動検出装置。
（付記８）変位測定手段による変位の測定結果に基づき、少なくとも対象ピクセルの各々に対応づけて、当該ピクセルが対応する地点における変動による位相差を示す情報または該位相差の大きさに応じた画素情報を保持する画像である差分干渉画像を生成する差分干渉画像生成手段を備えた付記１から付記７のいずれかに記載の変動検出装置。
本発明は、測定点の情報取得用途や、変動検出用途や、ＰＳ−ＩｎＳＡＲを利用する用途に限らず、ある点における電波の散乱波における反射強度および位相を示す情報に基づく当該点における位相の変化を利用するものであれば、好適に適用可能である。
１ 測定点情報提供装置
２ 飛翔体
３ 合成開口レーダ
４ 撮影視野
５ 変動検出装置
１０１ ＳＡＲ画像組入力手段
１０２ 最適偏波組検出手段
１０３ 測定点情報記憶手段
１０４ 偏波回転処理手段
１０５、１０５Ａ、１０５Ｂ、１０５Ｃ、１０５Ｄ 変動検出手段
１０６ 検出結果出力手段
１０７ 測定点情報出力手段
１０８ 測定点候補入力手段
１０９ 測定点条件判定手段
１１０ 判定結果出力手段
１５１ 回転ＳＡＲデータ生成手段
１５２Ａ、１５２Ｂ、１５２Ｃ 変位測定手段
１５３Ｂ、１５３Ｄ 回転ＳＡＲ画像生成手段
６０１ ＳＡＲ画像組入力手段
６０２ 偏波組決定手段
６０３ 測定点情報生成手段

Claims (11)

  1. 特定領域を含む撮影視野内の分解能セルに対応するピクセルの各々に対応づけて、略同一時刻に合成開口レーダが観測した送受信偏波の組合せがＨＨ偏波、ＨＶ偏波、ＶＨ偏波およびＶＶ偏波である基本４偏波組による観測データから生成される反射強度および位相を示す情報を少なくとも保持するＳＡＲ画像組を入力するＳＡＲ画像組入力手段と、
    入力されたＳＡＲ画像組を用いて、対象ピクセルの各々に対して、反射強度が所定値以上または最も強い偏波組を決定する偏波組決定手段と、
    前記偏波組決定手段により決定された偏波組に基づいて、少なくとも対象ピクセルの決定された偏波組を示す情報を含む測定点情報を生成する測定点情報生成手段と、
    測定点の候補を入力する測定点候補入力手段と、
    前記測定点情報に基づいて、入力された測定点の候補が測定点として適するか否かを、当該候補に対応するピクセルの決定された偏波組による反射強度を少なくとも用いて判定する測定点条件判定手段と、
    前記測定点条件判定手段による判定結果を出力する判定結果出力手段とを備え、
    前記偏波組決定手段は、入力された測定点の候補に対応するピクセルである対象ピクセルの各々に対して、反射強度が所定値以上または最も強い偏波組を決定し、
    前記測定点情報生成手段は、前記偏波組決定手段による偏波組の決定結果に基づき、少なくとも前記対象ピクセルの決定された偏波組を示す情報と、当該対象ピクセルの前記偏波組による反射強度を示す情報とを含む測定点情報を生成する
    ことを特徴とする測定点情報提供装置。
  2. 前記測定点情報生成手段が生成した測定点情報、または前記測定点情報に基づいて生成される、対象ピクセルに対応する地点が測定点に適するか否かもしくはどの程度測定点に適するかを示す情報を出力する測定点情報出力手段を備えた
    請求項１に記載の測定点情報提供装置。
  3. 前記測定点情報出力手段は、入力されたＳＡＲ画像組と、測定点情報とに基づいて、対象ピクセルに対応する地点がどの程度測定点に適するかを示す情報として、少なくとも対象ピクセルに対応づけて、当該ピクセルの決定された偏波組による反射強度および位相を示す情報を少なくとも保持する回転ＳＡＲ画像を生成して、出力する
    請求項２に記載の測定点情報提供装置。
  4. 特定領域を含む撮影視野内の分解能セルに対応するピクセルの各々に対応づけて、略同一時刻に合成開口レーダが観測した送受信偏波の組合せがＨＨ偏波、ＨＶ偏波、ＶＨ偏波およびＶＶ偏波である基本４偏波組による観測データから生成される反射強度および位相を示す情報を少なくとも保持するＳＡＲ画像組を２組以上入力するＳＡＲ画像組入力手段と、
    測定点の候補を入力する測定点候補入力手段と、
    入力されたＳＡＲ画像組のうちの少なくとも１組を用いて、入力された測定点の候補に対応するピクセルである対象ピクセルの各々に対して、反射強度が所定値以上または最も強い偏波組を決定する偏波組決定手段と、
    前記偏波組決定手段により決定された偏波組に基づいて、少なくとも前記対象ピクセルの決定された偏波組を示す情報と、当該対象ピクセルの前記偏波組による反射強度を示す情報とを含む測定点情報を生成する測定点情報生成手段と、
    前記測定点情報に基づいて、入力された測定の候補が測定点として適するか否かを、当該候補に対応するピクセルの決定された偏波組による反射強度を少なくとも用いて判定する測定点条件判定手段と、
    前記測定点条件判定手段による判定結果を出力する判定結果出力手段と、
    測定点として適すると判定された候補、または前記出力後に測定点に指定された候補を測定点に用いて、前記測定点情報と、入力されたＳＡＲ画像組とに基づいて、前記測定点における変位を測定する変位測定手段と、
    前記変位測定手段による測定結果を示す情報を出力する第１の出力手段と、
    前記測定点情報、または前記測定点情報に基づいて生成される、前記対象ピクセルに対応する地点が測定点に適するか否かもしくはどの程度測定点に適するかを示す情報を出力する第２の出力手段とを備えた
    ことを特徴とする変動検出装置。
  5. 特定領域を含む撮影視野内の分解能セルに対応するピクセルの各々に対応づけて、略同一時刻に合成開口レーダが観測した送受信偏波の組合せがＨＨ偏波、ＨＶ偏波、ＶＨ偏波およびＶＶ偏波である基本４偏波組による観測データから生成される反射強度および位相を示す情報を少なくとも保持するＳＡＲ画像組を２組以上入力するＳＡＲ画像組入力手段と、
    測定点の候補を入力する測定点候補入力手段と、
    入力されたＳＡＲ画像組のうちの少なくとも１組を用いて、入力された測定点の候補に対応するピクセルである対象ピクセルの各々に対して、反射強度が所定値以上または最も強い偏波組を決定する偏波組決定手段と、
    前記偏波組決定手段により決定された偏波組に基づいて、少なくとも前記対象ピクセルの決定された偏波組を示す情報と、当該対象ピクセルの前記偏波組による反射強度を示す情報とを含む測定点情報を生成する測定点情報生成手段と、
    前記測定点情報に基づいて、入力された測定の候補が測定点として適するか否かを、当該候補に対応するピクセルの決定された偏波組による反射強度を少なくとも用いて判定する測定点条件判定手段と、
    前記測定点条件判定手段による判定結果を出力する判定結果出力手段と、
    前記測定点条件判定手段によって測定点として適すると判定された候補を測定点に用いて、前記測定点情報と、入力されたＳＡＲ画像組とに基づいて、前記測定点における変位を測定する変位測定手段とを備えた
    ことを特徴とする変動検出装置。
  6. 特定領域を含む撮影視野内の分解能セルに対応するピクセルの各々に対応づけて、略同一時刻に合成開口レーダが観測した送受信偏波の組合せがＨＨ偏波、ＨＶ偏波、ＶＨ偏波およびＶＶ偏波である基本４偏波組による観測データから生成される反射強度および位相を示す情報を少なくとも保持するＳＡＲ画像組を入力し、
    測定点の候補を入力し、
    入力されたＳＡＲ画像組のうちの少なくとも１組を用いて、入力された測定点の候補に対応するピクセルである対象ピクセルの各々に対して、反射強度が所定値以上または最も強い偏波組を決定し、
    決定された偏波組に基づいて、少なくとも前記対象ピクセルの決定された偏波組を示す情報と、当該対象ピクセルの前記偏波組による反射強度を示す情報とを含む測定点情報を生成し、
    前記測定点情報に基づいて、入力された測定の候補が測定点として適するか否かを、当該候補に対応するピクセルの決定された偏波組による反射強度を少なくとも用いて判定し、
    前記判定結果を出力する
    ことを特徴とする測定点情報提供方法。
  7. 特定領域を含む撮影視野内の分解能セルに対応するピクセルの各々に対応づけて、略同一時刻に合成開口レーダが観測した送受信偏波の組合せがＨＨ偏波、ＨＶ偏波、ＶＨ偏波およびＶＶ偏波である基本４偏波組による観測データから生成される反射強度および位相を示す情報を少なくとも保持するＳＡＲ画像組を２組以上入力し、
    測定点の候補を入力し、
    入力されたＳＡＲ画像組を用いて、入力された測定点の候補に対応するピクセルである対象ピクセルの各々に対して、反射強度が所定値以上または最も強い偏波組を決定し、
    決定された偏波組に基づいて、少なくとも前記対象ピクセルの決定された偏波組を示す情報と、当該対象ピクセルの前記偏波組による反射強度を示す情報とを含む測定点情報を生成し、
    前記測定点情報に基づいて、入力された測定点の候補が測定点として適するか否かを、当該候補に対応するピクセルの前記偏波組による反射強度を少なくとも用いて判定し、
    前記判定結果を出力し、
    測定点として適すると判定された候補、または前記出力後に測定点に指定された候補を測定点に用いて、前記測定点情報と、入力されたＳＡＲ画像組とに基づいて、前記測定点における変位を測定し、
    前記測定点における変位の測定結果を示す情報とともに、前記測定点情報、または前記測定点情報に基づいて生成される、前記対象ピクセルに対応する地点が測定点に適するか否かもしくはどの程度測定点に適するかを示す情報を出力する
    ことを特徴とする変動検出方法。
  8. 特定領域を含む撮影視野内の分解能セルに対応するピクセルの各々に対応づけて、略同一時刻に合成開口レーダが観測した送受信偏波の組合せがＨＨ偏波、ＨＶ偏波、ＶＨ偏波およびＶＶ偏波である基本４偏波組による観測データから生成される反射強度および位相を示す情報を少なくとも保持するＳＡＲ画像組を２組以上入力し、
    測定点の候補を入力し、
    入力されたＳＡＲ画像組を用いて、入力された測定点の候補に対応するピクセルである対象ピクセルの各々に対して、反射強度が所定値以上または最も強い偏波組を決定し、
    決定された偏波組に基づいて、少なくとも前記対象ピクセルの決定された偏波組を示す情報と、当該対象ピクセルの前記偏波組による反射強度を示す情報とを含む測定点情報を生成し、
    前記測定点情報に基づいて、入力された測定点の候補が測定点として適するか否かを、当該候補に対応するピクセルの前記偏波組による反射強度を少なくとも用いて判定し、
    前記判定により測定点として適さないと判定された場合に、その旨を出力し、
    前記判定により測定点として適すると判定された場合に、当該候補を測定点に用いて、前記測定点情報と、入力されたＳＡＲ画像組とに基づいて、前記測定点における変位を測定する
    ことを特徴とする変動検出方法。
  9. コンピュータに、
    特定領域を含む撮影視野内の分解能セルに対応するピクセルの各々に対応づけて、略同一時刻に合成開口レーダが観測した送受信偏波の組合せがＨＨ偏波、ＨＶ偏波、ＶＨ偏波およびＶＶ偏波である基本４偏波組による観測データから生成される反射強度および位相を示す情報を少なくとも保持するＳＡＲ画像組を入力するＳＡＲ画像組入力処理、
    測定点の候補を入力する測定点候補入力処理、
    入力されたＳＡＲ画像組を用いて、入力された測定点の候補に対応するピクセルである対象ピクセルの各々に対して、反射強度が所定値以上または最も強い偏波組を決定する偏波組決定処理、
    前記偏波組決定処理で決定された偏波組に基づいて、少なくとも対象ピクセルの決定された偏波組を示す情報と、当該対象ピクセルの前記偏波組による反射強度を示す情報とを含む測定点情報を生成する測定点情報生成処理、
    前記測定点情報に基づいて、入力された測定点の候補が測定点として適するか否かを、当該候補に対応するピクセルの決定された偏波組による反射強度を少なくとも用いて判定する測定点条件判定処理、および
    前記測定点条件判定処理による判定結果を出力する判定結果出力処理
    を実行させるための測定点情報提供用プログラム。
  10. コンピュータに、
    特定領域を含む撮影視野内の分解能セルに対応するピクセルの各々に対応づけて、略同一時刻に合成開口レーダが観測した送受信偏波の組合せがＨＨ偏波、ＨＶ偏波、ＶＨ偏波およびＶＶ偏波である基本４偏波組による観測データから生成される反射強度および位相を示す情報を少なくとも保持するＳＡＲ画像組を２組以上入力するＳＡＲ画像組入力処理、
    測定点の候補を入力する測定点候補入力処理、
    入力されたＳＡＲ画像組を用いて、入力された測定点の候補に対応するピクセルである対象ピクセルの各々に対して、反射強度が所定値以上または最も強い偏波組を決定する偏波組決定処理、
    前記偏波組決定処理で決定された偏波組に基づいて、少なくとも前記対象ピクセルの決定された偏波組を示す情報と、当該対象ピクセルの前記偏波組による反射強度を示す情報とを含む測定点情報を生成する測定点情報生成処理、
    前記測定点情報に基づいて、入力された測定点の候補が測定点として適するか否かを、当該候補に対応するピクセルの決定された偏波組による反射強度を少なくとも用いて判定する測定点条件判定処理、
    前記測定点条件判定処理による判定結果を出力する判定結果出力処理、
    測定点として適すると判定された候補、または前記出力後に測定点に指定された候補を測定点に用いて、前記測定点情報と、入力されたＳＡＲ画像組とに基づいて、前記測定点における変位を測定する変位測定処理、
    前記測定点における変位の測定結果を示す情報を出力する第１の出力処理、および
    前記測定点情報、または前記測定点情報に基づいて生成される、前記対象ピクセルに対応する地点における地点が測定点に適するか否かもしくはどの程度測定点に適するかを示す情報を出力する第２の出力処理
    を実行させるための変動検出用プログラム。
  11. コンピュータに、
    特定領域を含む撮影視野内の分解能セルに対応するピクセルの各々に対応づけて、略同一時刻に合成開口レーダが観測した送受信偏波の組合せがＨＨ偏波、ＨＶ偏波、ＶＨ偏波およびＶＶ偏波である基本４偏波組による観測データから生成される反射強度および位相を示す情報を少なくとも保持するＳＡＲ画像組を２組以上入力するＳＡＲ画像組入力処理、
    測定点の候補を入力する測定点候補入力処理、
    入力されたＳＡＲ画像組を用いて、入力された測定点の候補に対応するピクセルである対象ピクセルの各々に対して、反射強度が所定値以上または最も強い偏波組を決定する偏波組決定処理、
    前記偏波組決定処理で決定された偏波組に基づいて、少なくとも前記対象ピクセルの決定された偏波組を示す情報と、当該対象ピクセルの前記偏波組による反射強度を示す情報とを含む測定点情報を生成する測定点情報生成処理、
    前記測定点情報に基づいて、入力された測定点の候補が測定点として適するか否かを、当該候補に対応するピクセルの決定された偏波組による反射強度を少なくとも用いて判定する測定点条件判定処理、
    前記測定点条件判定処理による判定結果を出力する判定結果出力処理、および
    前記測定点条件判定処理で測定点として適すると判定された候補を測定点に用いて、前記測定点情報と、入力されたＳＡＲ画像組とに基づいて、前記測定点における変位を測定する変位測定処理
    を実行させるための変動検出用プログラム。
JP2014097644A 2014-05-09 2014-05-09 測定点情報提供装置、変動検出装置、方法およびプログラム Active JP6349938B2 (ja)

Priority Applications (3)

Application Number Priority Date Filing Date Title
JP2014097644A JP6349938B2 (ja) 2014-05-09 2014-05-09 測定点情報提供装置、変動検出装置、方法およびプログラム
EP15165444.9A EP2942638B1 (en) 2014-05-09 2015-04-28 Measuring point information providing device, change detection device, methods thereof, and recording medium
US14/704,097 US10338215B2 (en) 2014-05-09 2015-05-05 Measuring point information providing device, change detection device, methods thereof, and recording medium

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2014097644A JP6349938B2 (ja) 2014-05-09 2014-05-09 測定点情報提供装置、変動検出装置、方法およびプログラム

Publications (2)

Publication Number Publication Date
JP2015215211A JP2015215211A (ja) 2015-12-03
JP6349938B2 true JP6349938B2 (ja) 2018-07-04

Family

ID=53180522

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2014097644A Active JP6349938B2 (ja) 2014-05-09 2014-05-09 測定点情報提供装置、変動検出装置、方法およびプログラム

Country Status (3)

Country Link
US (1) US10338215B2 (ja)
EP (1) EP2942638B1 (ja)
JP (1) JP6349938B2 (ja)

Families Citing this family (15)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US10871561B2 (en) 2015-03-25 2020-12-22 Urthecast Corp. Apparatus and methods for synthetic aperture radar with digital beamforming
CA2990063A1 (en) 2015-06-16 2017-03-16 King Abdulaziz City Of Science And Technology Efficient planar phased array antenna assembly
US10955546B2 (en) 2015-11-25 2021-03-23 Urthecast Corp. Synthetic aperture radar imaging apparatus and methods
CN106355625B (zh) * 2016-08-31 2019-12-24 甘宗平 一种单极化sar假彩色影像生成方法
JP7020418B2 (ja) * 2016-09-20 2022-02-16 日本電気株式会社 情報処理装置、情報処理方法、およびプログラム
CA3064735C (en) 2017-05-23 2022-06-21 Urthecast Corp. Synthetic aperture radar imaging apparatus and methods
US11378682B2 (en) 2017-05-23 2022-07-05 Spacealpha Insights Corp. Synthetic aperture radar imaging apparatus and methods for moving targets
US11460535B2 (en) * 2017-08-24 2022-10-04 Iceye Oy System and method for transmitting information from synthetic aperture radar satellite to client receiver
CA3083033A1 (en) 2017-11-22 2019-11-28 Urthecast Corp. Synthetic aperture radar apparatus and methods
US11460573B2 (en) * 2017-12-18 2022-10-04 Nec Corporation Synthetic aperture radar signal processing device and method
CN108279404B (zh) * 2018-01-22 2021-06-08 西安电子科技大学 一种基于空间谱估计的双通道sar相位误差校正方法
CN108957454B (zh) * 2018-08-01 2020-10-27 中国科学院电子学研究所 机载干涉合成孔径雷达的干涉相位仿真方法
WO2020202259A1 (ja) * 2019-03-29 2020-10-08 日本電気株式会社 合成開口レーダの画像処理装置及び画像処理方法
CN113408547B (zh) * 2021-07-12 2023-05-23 西南交通大学 一种多时相多极化sar滑坡提取方法
CN115015928B (zh) * 2022-05-30 2023-08-11 中国矿业大学 基于总功率同散射机制的高效极化时序InSAR方法

Family Cites Families (21)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5053778A (en) * 1989-08-10 1991-10-01 The United States Of America As Represented By The Administrator Of The National Aeronautics And Space Administration Generation of topographic terrain models utilizing synthetic aperture radar and surface level data
US5313210A (en) * 1993-02-23 1994-05-17 Ball Corporation Polarimetric radar signal mapping process
US5552787A (en) * 1995-10-10 1996-09-03 The United States Of America As Represented By The Secretary Of The Navy Measurement of topography using polarimetric synthetic aperture radar (SAR)
JP2697720B2 (ja) * 1995-12-27 1998-01-14 日本電気株式会社 ポラリメトリ合成開口レーダ装置
GB2354655A (en) * 1999-09-23 2001-03-28 Matra Marconi Space Uk Ltd Mitigation of Faraday rotation in space bourne radar
JP3964218B2 (ja) * 2002-01-29 2007-08-22 三菱電機株式会社 合成開口レーダ装置及び像再生方法
DE10226508A1 (de) * 2002-06-14 2004-01-08 Dornier Gmbh Verfahren zur Detektion sowie Geschwindigkeits- und Positionsschätzung von bewegten Objekten in SAR-Bildern
US6781540B1 (en) * 2003-02-21 2004-08-24 Harris Corporation Radar system having multi-platform, multi-frequency and multi-polarization features and related methods
JP2007078398A (ja) * 2005-09-12 2007-03-29 Toshiba Corp 合成開口レーダ画像表示装置
US7623064B2 (en) * 2005-12-06 2009-11-24 Arthur Robert Calderbank Instantaneous radar polarimetry
US7492303B1 (en) * 2006-05-09 2009-02-17 Personnel Protection Technologies Llc Methods and apparatus for detecting threats using radar
US9019143B2 (en) * 2006-11-30 2015-04-28 Henry K. Obermeyer Spectrometric synthetic aperture radar
JP4946547B2 (ja) * 2007-03-16 2012-06-06 日本電気株式会社 ポラリメトリｓａｒ画像処理方法、ポラリメトリｓａｒ装置および画像処理装置
US8125370B1 (en) * 2007-04-16 2012-02-28 The United States Of America As Represented By The Secretary Of The Navy Polarimetric synthetic aperture radar signature detector
US7746267B2 (en) * 2007-05-08 2010-06-29 The Johns Hopkins University Synthetic aperture radar hybrid-polarity method and architecture for obtaining the stokes parameters of a backscattered field
JP5040549B2 (ja) * 2007-09-20 2012-10-03 日本電気株式会社 合成開口レーダ及びコンパクト・ポラリメトリｓａｒ処理方法、プログラム
CN102144174B (zh) * 2008-07-04 2015-04-15 电视广播有限公司 Sar图像序列中的永久散射体的识别和分析
US8154435B2 (en) * 2008-08-22 2012-04-10 Microsoft Corporation Stability monitoring using synthetic aperture radar
JP5305985B2 (ja) 2009-02-27 2013-10-02 三菱スペース・ソフトウエア株式会社 人工物検出装置及び人工物検出方法及び人工物検出プログラム
US10107904B2 (en) * 2012-09-04 2018-10-23 Fugro N.V. Method and apparatus for mapping and characterizing sea ice from airborne simultaneous dual frequency interferometric synthetic aperture radar (IFSAR) measurements
US9395437B2 (en) * 2013-06-06 2016-07-19 The United States Of America, As Represented By The Secretary Of The Army Moving multi-polarization multi-transmitter/receiver ground penetrating radar system and signal processing for buried target detection

Also Published As

Publication number Publication date
EP2942638B1 (en) 2020-08-05
US10338215B2 (en) 2019-07-02
US20150323665A1 (en) 2015-11-12
JP2015215211A (ja) 2015-12-03
EP2942638A1 (en) 2015-11-11

Similar Documents

Publication Publication Date Title
JP6349937B2 (ja) 変動検出装置、変動検出方法および変動検出用プログラム
JP6349938B2 (ja) 測定点情報提供装置、変動検出装置、方法およびプログラム
US8125370B1 (en) Polarimetric synthetic aperture radar signature detector
US8154435B2 (en) Stability monitoring using synthetic aperture radar
La et al. Data analysis and visualization for the bridge deck inspection and evaluation robotic system
US9261583B2 (en) Systems and methods for calibrating dual polarization radar systems
US9097805B2 (en) Systems and methods for calibrating dual polarization radar systems
CN109116321B (zh) 一种星载干涉成像高度计的相位滤波方法及高度测量方法
Morrison et al. Laboratory measurement of the DInSAR response to spatiotemporal variations in soil moisture
CN103400137B (zh) 一种sar图像的建筑物几何参数提取方法
US7446705B1 (en) Method and apparatus for determining parameters for a parametric expression characterizing the phase of an acquired signal
CN110441770A (zh) 基于多部mimo-sar联合观测的三维形变测量方法
CN106291542A (zh) 一种隧道三维成像方法
KR101784178B1 (ko) 해양변위 관측용 산란계 시스템
CN105241428A (zh) 一种利用高光谱进行水深反演的方法
Ortiz Arteaga et al. Initial investigation of a low-cost automotive lidar system
Sharma et al. Radiometric calibration stability assessment for the RISAT-1 SAR sensor using a deployed point target array at the Desalpar site, Rann of Kutch, India
US8437509B1 (en) System and method for inferring vessel speed from overhead images
Catapano et al. Structural monitoring via microwave tomography-enhanced GPR: the Montagnole test site
CN108710816B (zh) 一种基于概率融合的极化sar图像建筑物检测方法
Rasol et al. Analysis and calibration of ground penetrating radar shielded antennas
Goel et al. Fusion of monostatic/bistatic InSAR stacks for urban area analysis via distributed scatterers
TWI474029B (zh) 應用微波雷達於海岸線及潮間帶地形測量之技術
RU2630412C1 (ru) Способ дистанционного определения скорости морского течения
JP6234707B2 (ja) 水位計測装置、水位計測プログラム、水位計測方法およびレーダ観測方法

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20170406

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20180214

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20180220

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20180418

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

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20180521

R150 Certificate of patent or registration of utility model

Ref document number: 6349938

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150