(57)【要約】【課題】地表のなめらかな変化と断層等の不連続的な変化との区別が容易な地表変動表現装置及び地表変動表現プログラムを提供する。【解決手段】位相差情報取得部１０が、異なる時刻にレーダー装置により取得された観測対象である地表の時系列データに対して干渉解析を行って生成された位相差情報を取得し、画像化部１２は、位相差の値の範囲の両端で色相が一致し、かつ上記範囲の位相差の任意の値に基準色を設定して各基準色の間で位相差の変化に応じてＲＧＢの値を線形に変化させているパレットに基づき、上記位相差情報に含まれる位相差をピクセル値として有する干渉画像の各ピクセルの位相差に、上記パレットにおいて設定されている色相を対応させて着色画像を生成する。【選択図】図１
2020143910.tif 000002
【特許請求の範囲】
【請求項１】
異なる時刻にレーダー装置により取得された観測対象である地表の時系列データに対して干渉解析を行って生成された位相差情報を取得する位相差情報取得手段と、
前記位相差情報に含まれ、一定範囲の値をとる位相差に応じて画像の色相を決定するパレットであって、前記位相差の値の範囲の両端で色相が一致し、かつ前記範囲の位相差の任意の値に基準色を設定して各基準色の間でＲＧＢの値を線形に変化させているパレットを有し、前記位相差情報を画像化する画像化手段と、
を備える地表変動表現装置。
【請求項２】
前記基準色が位相差の４つの値に対して設定される、請求項１に記載の地表変動表現装置。
【請求項３】
前記位相差の値の範囲が−πからπであり、前記基準色がπ、π／２、０、−π／２である、請求項２に記載の地表変動表現装置。
【請求項４】
前記位相差を三角関数により連続化する連続化手段をさらに有し、
前記画像化手段は、前記三角関数の振幅と接線の傾きとの組み合わせまたは位相差に応じて画像の色相を決定するパレットであって、前記三角関数の１周期の始点と終点で色相が一致し、かつ前記三角関数の振幅と接線の傾きとの組み合わせまたは位相差が異なる任意の点に基準色を設定して各基準色の間でＲＧＢの値を線形に変化させているパレットを有し、前記位相差情報を画像化する、請求項１に記載の地表変動表現装置。
【請求項５】
コンピュータを、
異なる時刻にレーダー装置により取得された観測対象である地表の時系列データに対して干渉解析を行って生成された位相差情報を取得する位相差情報取得手段、
前記位相差情報に含まれ、一定範囲の値をとる位相差に応じて画像の色相を決定するパレットであって、前記位相差の値の範囲の両端で色相が一致し、かつ前記範囲の位相差の任意の値に基準色を設定して各基準色の間でＲＧＢの値を線形に変化させているパレットを有し、前記位相差情報を画像化する画像化手段、
として機能させる地表変動表現プログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、地表変動表現装置及び地表変動表現プログラムに関する。
【背景技術】
【０００２】
従来、異なる時刻に合成開口レーダー（ＳＡＲ）で観測したデータに対して干渉解析を行い（干渉ＳＡＲ）、上記異なる時刻における地表の変動を把握することが行われている（下記特許文献１、特許文献２参照）。
【０００３】
このような干渉ＳＡＲは、２枚のＳＡＲ画像のピクセル毎に観測値として与えられた位相情報の差分を求める技術であり、干渉画像の各ピクセルには上記位相の差分（位相差）として−π〜πの値が与えられる。この位相差は、ＳＡＲで使用する電波の半波長（測定用の電波が往復するので）分の地表変動毎に−πからπまで変化し、地表変動が半波長を超える度にその値がπから−πに戻って値の変化が繰り返される。また、上記位相差は、地形に変化がないときには０（ゼロ）付近の値を取る。ここで、干渉画像の各ピクセルに位相差として与えられた−π〜πの値に適宜な色相を対応させるパレットを設定し、このパレットに基づいて上記干渉画像を画像化（可視化）することができる。このような可視化された干渉画像により、地震等の災害により生じた地形変化箇所を位相差の不連続性から特定することができる。
【０００４】
例えば、地形が緩やかに隆起あるいは沈降した場合、干渉ＳＡＲにおいて−π〜πで与えられる位相差を表す上記干渉画像は、その地域全体がグラデーションをなすように連続的な画像となり、あたかも山や谷の等高線のように表示される。一方、断層のような地形上の不連続箇所は位相差も不連続な値となり、上記可視化された干渉画像も不連続的な画像となるため、この不連続箇所を干渉画像から読み取れば、断層の位置を特定することができる。
【０００５】
しかし、上述したように、位相差として−π〜πの値が与えられるため、これを干渉画像として従来のパレット（位相差に色相を対応させるパレット）を使用して画像化すると、地形が緩やかに隆起あるいは沈降して連続的な変化をしている領域であっても、位相差がπから−πに移行する部分で色相も大きく変化し、断層等に基づく不連続箇所との区別が困難になるという問題がある。
【０００６】
そこで、位相差を三角関数に代入し、連続関数化することで−πとπの２つの値の間に生じる不連続性を解決する方法が考えられる。この方法によれば、干渉画像として位相差を画像化するのではなく、三角関数に代入して得られる振幅に色相を対応させるパレットを設定し、このパレットに基づいて位相差を画像化（可視化）するので、上記−πとπの２つの値の間に生じる不連続性を緩和することができる。
【０００７】
しかし、この方法によっても、例えば三角関数としてｃｏｓ（コサイン）関数を使用した場合、位相差が−π／４とπ／４等のように振幅が同じとなる場合があり、上記振幅に色相を対応させるパレットを使用すると、同じ振幅に対応する色相が同じとなるので、両者を識別できない（解像度が低下する）という課題がある。
【先行技術文献】
【特許文献】
【０００８】
【特許文献１】特開２０１７−２０７４５７号公報
【特許文献２】特開２０１８−４０７２８号公報
【発明の概要】
【発明が解決しようとする課題】
【０００９】
本発明の目的は、地表のなめらかな変化と断層等の不連続的な変化との区別が容易な地表変動表現装置及び地表変動表現プログラムを提供することにある。
【課題を解決するための手段】
【００１０】
上記目的を達成するために、本発明の一実施形態は、地表変動表現装置であって、異なる時刻にレーダー装置により取得された観測対象である地表の時系列データに対して干渉解析を行って生成された位相差情報を取得する位相差情報取得手段と、前記位相差情報に含まれ、一定範囲の値をとる位相差に応じて画像の色相を決定するパレットであって、前記位相差の値の範囲の両端で色相が一致し、かつ前記範囲の位相差の任意の値に基準色を設定して各基準色の間でＲＧＢの値を線形に変化させているパレットを有し、前記位相差情報を画像化する画像化手段と、を備えることを特徴とする。
【００１１】
上記基準色は、位相差の４つの値に対して設定されるのが好適である。
【００１２】
また、上記位相差の値の範囲が−πからπであり、前記基準色は、π、π／２、０、−π／２であるのが好適である。
【００１３】
また、上記位相差を三角関数により連続化する連続化手段をさらに有し、前記画像化手段は、前記三角関数の振幅と接線の傾きとの組み合わせまたは位相差に応じて画像の色相を決定するパレットであって、前記三角関数の１周期の始点と終点で色相が一致し、かつ前記三角関数の振幅と接線の傾きとの組み合わせまたは位相差が異なる任意の点に基準色を設定して各基準色の間でＲＧＢの値を線形に変化させているパレットを有し、前記位相差情報を画像化するのが好適である。
【００１４】
本発明の他の実施形態は、地表変動表現プログラムであって、コンピュータを、異なる時刻にレーダー装置により取得された観測対象である地表の時系列データに対して干渉解析を行って生成された位相差情報を取得する位相差情報取得手段、前記位相差情報に含まれ、一定範囲の値をとる位相差に応じて画像の色相を決定するパレットであって、前記位相差の値の範囲の両端で色相が一致し、かつ前記範囲の位相差の任意の値に基準色を設定して各基準色の間でＲＧＢの値を線形に変化させているパレットを有し、前記位相差情報を画像化する画像化手段、として機能させることを特徴とする。
【発明の効果】
【００１５】
本発明によれば、地表のなめらかな変化と断層等の不連続的な変化との区別を容易に行うことができる。
【図面の簡単な説明】
【００１６】
【図１】実施形態にかかる地表変動表現装置の例の機能ブロック図である。
【図２】実施形態にかかる画像化部が有するパレットの例を示す図である。
【図３】実施形態にかかる連続化部の動作の説明図である。
【図４】実施形態にかかる三角関数の振幅と接線の傾きとの組合せまたは位相差に応じて画像の色相を決定するパレットの例を示す図である。
【図５】実施形態にかかる地表変動表現装置により地表面の変動を画像化した例及び比較例の地表面の変動を画像化した例を示す図である。
【図６】実施形態にかかる地表変動表現装置の動作例のフロー図である。
【発明を実施するための形態】
【００１７】
以下、本発明を実施するための形態（以下、実施形態という）を、図面に従って説明する。
【００１８】
図１には、実施形態にかかる地表変動表現装置１００の例の機能ブロック図が示される。図１において、地表変動表現装置１００は、位相差情報取得部１０、画像化部１２、連続化部１４、表示制御部１６、通信部１８、記憶部２０及びＣＰＵ２２を含んで構成されている。なお、ＣＰＵ２２以外にＧＰＵ等のアクセラレーターを用いてもよい。上記地表変動表現装置１００は、ＣＰＵ２２、ＲＯＭ、ＲＡＭ、不揮発性メモリ、Ｉ／Ｏ、通信インターフェース等を備え、装置全体の制御及び各種演算を行うコンピュータとして構成されており、上記各機能は、例えばＣＰＵ２２とＣＰＵ２２の処理動作を制御するプログラムとにより実現される。
【００１９】
位相差情報取得部１０は、異なる時刻にレーダー装置により取得された観測対象である地表の時系列データ（位相を含む情報であり、以後ＳＡＲ画像ということがある）に対して干渉解析を行って生成された位相差情報（以後干渉画像ということがある）を取得する。ここで、レーダー装置としては、例えばＳＡＲ（Synthetic Aperture Radar 合成開口レーダー）が挙げられるが、これには限定されない。また、上記時系列データは、レーダー装置により取得された、異なる時刻に地表の同一点から反射された電波の反射波を含んだ情報である。位相差情報取得部１０は、上記ＳＡＲ等により観測した時系列データに含まれる反射波の位相から干渉解析が行われて得られた位相差を含む情報（位相差情報）を、通信部１８を介して外部のサーバ等から取得することができる。取得した位相差情報（干渉画像）は、記憶部２０に記憶させる。位相差情報には、合成開口レーダーからの視線方向の距離変化の情報が含まれており、これを使用して地表面の高さの変化の情報を抽出することができる。
【００２０】
ここで、上記位相差情報に含まれる位相差は、２枚のＳＡＲ画像の差分により形成された干渉画像における各ピクセルのピクセル値として与えられ、記憶部２０に記憶される。前述したように、上記位相差は、ＳＡＲで使用する電波の半波長分の地表変動毎に−πからπまで変化し、地表変動が半波長を超える度にその値がπから−πに戻って値の変化が繰り返される。
【００２１】
画像化部１２は、上記位相差情報（干渉画像）を記憶部２０から読み出し、この干渉画像を可視化する。このため、画像化部１２は、上記位相差情報に含まれる位相差、すなわち干渉画像における各ピクセルのピクセル値として与えられた一定範囲の値をとる位相差に応じて画像の色相を決定するパレットであって、上記位相差の値の範囲の両端で色相が一致し、かつ上記範囲の位相差の任意の値に基準色を設定して各基準色の間で位相差の変化に応じてＲ（赤）Ｇ（緑）Ｂ（青）の値を線形に変化させているパレットを有しており、このパレットを使用して上記位相差情報を画像化する。なお、各基準色の色相は、互いに異ならせる。これにより、位相差の値の相違を画像の色として認識できるからである。
【００２２】
ここで、上記位相差情報に含まれる位相差の値は、上述した通り地表面の変動に応じて−πからπの範囲で値が繰り返される。従って、上記パレットでは、位相差の値の範囲の両端である−πとπに同じ色相が設定される。これにより、位相差の値が−πからπに移行する部分で色相が大きく変化することを回避できるので、断層等に基づく不連続箇所と位相差の値の−πからπへの移行箇所を区別する必要がなくなる。また、−πからπの間の位相差の任意の値に基準色が設定される。この基準色は、−πからπの間の位相差の４つの値に対して設定されるのが好適であり、例えばπ、π／２、０、−π／２等とすることができる。このような基準色の設定により、断層等の識別をし易い色相構成とすることができる。上記パレットは、予め作成して記憶部２０に記憶しておき、画像化部１２が読み出して使用する構成とする。
【００２３】
画像化部１２は、上記干渉画像及びパレットを記憶部２０から読み出し、干渉画像の各ピクセルのピクセル値として与えられている位相差に、パレットにおいて設定されている色相を対応させて干渉画像を可視化する。以後、可視化した干渉画像を着色画像ということがある。このように干渉画像に色相を対応させて着色画像を生成することが「位相差情報を画像化する」ことの意味である。画像化した位相差情報（着色画像）は、記憶部２０に記憶させる。なお、この着色画像は、画像化部１２が表示制御部１６に指示して適宜な表示装置に表示させる。
【００２４】
連続化部１４は、上記位相差情報を記憶部２０から読み出し、この位相差情報に含まれる位相差を三角関数により連続化する。連続化した位相差情報は、記憶部２０に記憶させる。この場合には、画像化部１２は、記憶部２０から連続化した位相差情報を読み出して着色画像を生成し、記憶部２０に記憶させる。連続化部１４の動作は後述する。
【００２５】
表示制御部１６は、上記着色画像を記憶部２０から読み出し、液晶表示装置その他の適宜な表示装置を制御して表示する。なお、連続化部１４が生成した連続化した位相差情報から画像化部１２が着色画像を生成した場合には、この着色画像を表示する。
【００２６】
通信部１８は、適宜なインターフェースにより構成され、ＳＡＲ衛星から時系列データ（位相を含む情報）を受信して干渉解析を行い位相差情報を抽出する外部のサーバ等から位相差情報を取得し、位相差情報取得部１０に渡す。
【００２７】
記憶部２０は、ハードディスク装置、ソリッドステートドライブ（ＳＳＤ）等の不揮発性メモリで構成され、上記各種情報等、及びＣＰＵ２２の動作プログラム等の、地表変動表現装置１００が行う各処理に必要な情報を記憶する。なお、記憶部２０としては、デジタル・バーサタイル・ディスク（ＤＶＤ）、コンパクトディスク（ＣＤ）、光磁気ディスク（ＭＯ）、フレキシブルディスク（ＦＤ）、磁気テープ、電気的消去および書き換え可能な読出し専用メモリ（ＥＥＰＲＯＭ）、フラッシュ・メモリ等を使用してもよい。また、記憶部２０には、主としてＣＰＵ２２の作業領域として機能するランダムアクセスメモリ（ＲＡＭ）、及びＢＩＯＳ等の制御プログラムその他のＣＰＵ２２が使用するデータが格納される読み出し専用メモリ（ＲＯＭ）を含めるのが好適である。
【００２８】
図２(ａ)、（ｂ）には、画像化部１２が有するパレットの例が示される。図２(ａ)の円の中心角が位相差の値を示しており、中心から図２(ａ)の右方向に水平に伸びる半径に対して反時計回りに０からπまで変化し、時計回りに０から−πまで変化する。すなわち、πと−πとが円周上の同じ点で重なっている。図２（ｂ）は、図２(ａ)に示された位相差の値に対応して色相が設定されている。図２（ｂ）の例では、色相が赤（Ｒ）、緑（Ｇ）、青（Ｂ）の値（０〜２５５）として示されている。なお、図２(ａ)には、図２（ｂ）に例示された、ＲＧＢの値で示された色相がグレースケールで表されている。
【００２９】
また、基準色は４つ設定されており、位相差がπ、π／２、０、−π／２の場合にそれぞれ基準色ａ、基準色ｂ、基準色ｃ、基準色ｄとなっている。なお、位相差が−πの基準色は位相差がπの時（基準色ａ）と同じである。これにより、位相差の値の範囲の両端（πと−π）で色相が一致する。
【００３０】
上記４つの基準色ａ〜ｄの各間（ａｂ間、ｂｃ間、ｃｄ間及びｄａ間）でＲＧＢの値は線形に変化するように設定される。図２（ｂ）の例では、基準色ａ〜ｄが設定された４つの位相差以外の位相差にも色相の数値（ＲＧＢの値）が記載されているが、基準色ａ〜ｄの各間でＲＧＢの値が線形に変化する条件を満たせば、図２（ｂ）の例に限定されない。
【００３１】
なお、パレットの例は、図２(ａ)、（ｂ）に示されたものに限定されず、位相差の値の範囲の両端で色相が一致し、かつ上記範囲の位相差の任意の値に基準色を設定して各基準色の間で位相差の変化に応じてＲＧＢの値を線形に変化させているものであればよい。
【００３２】
図３(ａ)、（ｂ）には、連続化部１４の動作の説明図が示される。図３(ａ)には、異なる時刻に取得された時系列データから干渉解析を行って得た位相差の変化の例が示される。図３(ａ)において、縦軸が位相差であり、横軸が地表面上の位置を表す。なお、地表面上の位置は、ＳＡＲ画像及び干渉画像におけるピクセルの位置に対応している。
【００３３】
一般に、上記異なる時刻に取得された時系列データから干渉解析を行って得た位相差は、上述した通り、地表変動に応じて−πからπまで変化し、地表変動が使用する電波の半波長を超える度にその値がπから−πに戻って値の変化が繰り返される。従って、地表変動が地表面上の位置の移動とともになめらかに大きくなるまたは小さくなるに従い、図３(ａ)に示されるように、位相差の変化がのこぎり状の波形となる。
【００３４】
そこで、連続化部１４は、のこぎり状の関数値（−πからπの値）を三角関数（例えばｓｉｎ関数またはｃｏｓ関数）に代入して正弦波波形に変換することにより、図３(ｂ)に示されるように、のこぎり波形を正弦波波形に連続化する。ここで、図３(ｂ) の縦軸は、使用した三角関数の振幅（無次元数）であり、図３(ｂ) の横軸は、位相差を表す。
【００３５】
この場合、上記画像化部１２が有するパレットは、上記三角関数の振幅と接線の傾きとの組合せまたは位相差に応じて画像の色相を決定するパレットとなっている。すなわち、図３(ｂ)に示される、上記三角関数の１周期の始点と終点で色相が一致しており、かつ三角関数の振幅と接線の傾きとの組み合わせまたは位相差が異なる任意の点に基準色を設定して各基準色の間でＲＧＢの値を線形に変化させている。
【００３６】
ここで、三角関数の振幅と接線の傾きとの組み合わせまたは位相差が異なる点の例としては、図３(ｂ)に示される始点または終点の一方、頂点１、頂点２、中点等が挙げられる。始点または終点の一方とするのは、終点が次の周期の始点だからである。なお、基準色を設定する点としてはこの４点とすることが好適である。
【００３７】
なお、基準色の色相は、互いに異ならせる。これにより、位相差の値の相違を画像の色として認識できるからである。さらに、上記基準色は、三角関数の振幅が同じで接線の傾きが異なる点、例えば始点（終点）と中点の色相が異なるように設定する。これらの点では振幅が同じであり、同じ色相にすると画像の解像度が低下するからである。
【００３８】
画像化部１２は、上記連続化した位相差情報の正弦波波形の振幅と接線の傾きとの組合せまたは位相差に応じ、上記パレットにおいて設定されている色相を対応させて干渉画像を可視化し、着色画像を生成する。
【００３９】
図４（ａ）、（ｂ）、（ｃ）には、三角関数の振幅と接線の傾きとの組合せまたは位相差に応じて画像の色相を決定するパレットの例が示される。図４（ａ）が振幅と接線の傾きとの組合せを使用した場合のパレットであり、図４（ｂ）が位相差を使用した場合のパレットである。また、図４（ｃ）には、振幅と接線の傾きとの組合せまたは位相差に対して設定された色相の例がＲＧＢの値（０〜２５５）として示されている。
【００４０】
三角関数は、その振幅と接線の傾きとの組み合わせですべての点を特定できるので、図４（ａ）及び図４（ｃ）に示されるように、パレットも振幅と接線の傾きとの組み合わせに応じて色相が決定されている。
【００４１】
接線の傾きについては、位相差と三角関数の関係から傾きを決める。例えば、位相差を正弦波に変換した場合、位相差が−π以上−π／２より小さい範囲と、π／２より大きくπ以下の範囲は、位相差の値が増えると振幅の値は減少するので傾きを−１とする。また、位相差が−π／２より大きくπ／２より小さい範囲は、位相差の値が増えると振幅の値は増加するので傾きを１とする。また位相差が−π／２とπ／２については、傾きを０とする。
【００４２】
また、三角関数は、その角度（図３（ｂ）の横軸で示された位相差）ですべての点を特定できるので、図４（ｂ）及び図４（ｃ）に示されるように、パレットを位相差に応じて色相を決定するように構成することもできる。
【００４３】
図５(ａ)、（ｂ）、（ｃ）、（ｄ）、（ｅ）、（ｆ）、（ｇ）には、実施形態にかかる地表変動表現装置１００により地表面の変動を画像化した例（図５(ａ) 、（ｄ）、（ｅ））及び比較例の地表面の変動を画像化した例（図５（ｂ）、（ｃ）、（ｆ）、（ｇ））が示される。
【００４４】
図５(ａ)の例は、位相差の値の範囲の両端である−πとπに同じ色相が設定され、位相差のπ、π／２、０、−π／２に対して基準色が設定されたパレットを使用して干渉画像を可視化した着色画像の例である。この例では、連続化部１４を使用していない。また、図５（ｄ）、（ｅ）は、図５(ａ)の白線で囲まれた領域の拡大図である。
【００４５】
図５(ａ)において、白線で囲まれた領域には断層線が観測されている。この断層線は、図５（ｄ）に破線で示されている。断層線の抽出は、図５（ｄ）の円で示された位置について、その周囲の位相差に比べて位相差の変化が不連続となっていることが把握できるので、上記円で示された位置の周辺に断層線が存在していると判断して行う。実施形態にかかる地表変動表現装置１００では、上記パレットを使用することにより、周期跨ぎ（πから−πに移行する部分）の際に色相が大きく変化せず、４つの区分に分けて色の値を線形に変化させているため位相差の値の相違を画像の色として認識し易くなるので、位相差の変化が不連続となっている箇所を容易に把握できる。
【００４６】
なお、図５（ｅ）は、図５（ｄ）から破線及び円を削除した図であり、図５（ｄ）と比べて理解を容易にするための参考図である。
【００４７】
以上に対して、図５（ｂ）は、従来のパレット（位相差に色相を対応させただけのパレット）を使用して干渉画像を可視化した着色画像の例であり、図５（ｆ）は、図５(ｂ)の白線で囲まれた領域の拡大図である。図５（ｆ）には、参考として断層線が破線で示されている。また、図５（ｃ）は、位相差を三角関数に代入して連続関数化し、その振幅に色相を対応させただけパレットを使用して干渉画像を可視化した着色画像の例であり、図５（ｇ）は、図５(ｃ)の白線で囲まれた領域の拡大図である。図５（ｇ）にも、参考として断層線が破線で示されている。
【００４８】
図５（ｂ）、（ｃ）、（ｆ）、（ｇ）からわかるように、従来の位相差に色相を対応させただけのパレットあるいは連続関数の振幅に色相を対応させただけのパレットを使用して干渉画像を可視化した場合、位相差の変化が不連続となっている箇所の判別がし難く、断層線等の把握が困難である。
【００４９】
図６には、実施形態にかかる地表変動表現装置１００の動作例のフロー図が示される。図６において、位相差情報取得部１０が、異なる時刻にレーダー装置により取得された観測対象である地表の時系列データ（位相を含む情報）に対して干渉解析を行って生成された位相差情報（干渉画像）を取得して記憶部２０に記憶させる（Ｓ１）。
【００５０】
画像化部１２は、位相差の値の範囲の両端で色相が一致し、かつ上記範囲の位相差の任意の値に基準色を設定して各基準色の間で位相差の変化に応じてＲＧＢの値を線形に変化させているパレットを記憶部２０から読み出す（Ｓ２）。
【００５１】
次に、画像化部１２は、上記干渉画像を記憶部２０から読み出し、干渉画像の各ピクセルのピクセル値として与えられている位相差に、上記パレットにおいて設定されている色相を対応させて着色画像を生成する。画像化部１２が生成した着色画像は、記憶部２０に記憶させる（Ｓ３）。
【００５２】
次に、画像化部１２は、上記生成した着色画像を記憶部２０から読み出し、表示制御部１６に指示して適宜な表示装置に表示させる（Ｓ４）。
【００５３】
上述した、図６の各ステップを実行するためのプログラムは、記録媒体に格納することも可能であり、また、そのプログラムを通信手段によって提供しても良い。その場合、例えば、上記説明したプログラムについて、「プログラムを記録したコンピュータ読み取り可能な記録媒体」の発明または「データ信号」の発明として捉えても良い。
【符号の説明】
【００５４】
１０  位相差情報取得部、１２  画像化部、１４  連続化部、１６  表示制御部、１８通信部、２０  記憶部、２２  ＣＰＵ、１００  地表変動表現装置。
【図１】
2020143910.tif 000003
【図２】
2020143910.tif 000004
【図３】
2020143910.tif 000005
【図４】
2020143910.tif 000006
【図５】
2020143910.tif 000007
【図６】
2020143910.tif 000008