(57)【要約】【課題】レーダ装置により取得した画像情報に含まれる変動範囲が広い変動の視認性の向上と変動範囲が狭い変動の視認性の向上を両立する。【解決手段】異なる時刻にレーダ装置により取得された地表の時系列データをもとに干渉解析が行われた結果の位相差情報を第１の位相差ノイズ除去部１２が第１の窓サイズを有するノイズ除去フィルタによりノイズ除去を行い、第２の位相差ノイズ除去部１４が上記位相差情報を第１の窓サイズより小さい第２の窓サイズを有するノイズ除去フィルタによりノイズ除去を行う。干渉縞画像生成部１８は、上記第２の窓サイズを有するノイズ除去フィルタによりノイズ除去を行った第２の位相差情報から干渉縞画像を生成する。表示制御部２０は、上記第１の窓サイズを有するノイズ除去フィルタによりノイズ除去を行った第１の位相差情報と上記干渉縞画像との位置情報に基づいて上記第１の位相差情報と上記干渉縞画像とを結合し表示する。【選択図】図１
2020106460.tif 000002
【特許請求の範囲】
【請求項１】
地表変動を表現するための地表変動表現装置であって、
異なる時刻にレーダ装置により取得された観測対象である地表の時系列データをもとに干渉解析が行われた結果の位相差情報を取得する位相差情報取得手段と、
第１の窓サイズを有し、前記位相差情報取得手段が取得した位相差情報からノイズ除去を行う第１の位相差ノイズ除去手段と、
前記第１の窓サイズより小さい第２の窓サイズを有し、前記位相差情報取得手段が取得した位相差情報からノイズ除去を行う第２の位相差ノイズ除去手段と、
前記第２の位相差ノイズ除去手段から得られた第２の位相差情報から干渉縞画像を生成する干渉縞画像生成手段と、
前記第１の位相差ノイズ除去手段から得られた第１の位相差情報と前記干渉縞画像との位置情報に基づき、前記第１の位相差ノイズ除去手段から得られた第１の位相差情報と前記干渉縞画像とを結合して表示する表示制御手段と、
を備える地表変動表現装置。
【請求項２】
前記表示制御手段は、前記第１の位相差ノイズ除去手段から得られた第１の位相差情報を３次元表示し、前記干渉縞画像を３次元表示した第１の位相差ノイズ除去手段から得られた第１の位相差情報に重畳して表示する、請求項１に記載の地表変動表現装置。
【請求項３】
前記第１の位相差ノイズ除去手段から得られた第１の位相差情報を連続化して連続化位相差情報を生成する連続化手段をさらに備え、
前記表示制御手段は、前記第１の位相差ノイズ除去手段から得られた第１の位相差情報の代わりに連続化位相差情報を使用する、請求項１または請求項２に記載の地表変動表現装置。
【請求項４】
地表変動を表現するための地表変動表現プログラムであって、コンピュータを、
異なる時刻にレーダ装置により取得された観測対象である地表の時系列データをもとに干渉解析が行われた結果の位相差情報を取得する位相差情報取得手段、
第１の窓サイズを有し、前記位相差情報取得手段が取得した位相差情報からノイズ除去を行う第１の位相差ノイズ除去手段、
前記第１の窓サイズより小さい第２の窓サイズを有し、前記位相差情報取得手段が取得した位相差情報からノイズ除去を行う第２の位相差ノイズ除去手段、
前記第２の位相差ノイズ除去手段から得られた第２の位相差情報から干渉縞画像を生成する干渉縞画像生成手段、
前記第１の位相差ノイズ除去手段から得られた第１の位相差情報と前記干渉縞画像との位置情報に基づき、前記第１の位相差ノイズ除去手段から得られた第１の位相差情報と前記干渉縞画像とを結合して表示する表示制御手段、
として機能させることを特徴とする地表変動表現プログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、地表変動表現装置及びプログラムに関する。
【背景技術】
【０００２】
従来、地震等により地表面に変動が発生した場合、地表面の変動前に合成開口レーダ（ＳＡＲ）で観測したデータと地表面の変動後に合成開口レーダ（ＳＡＲ）で観測したデータから、干渉解析を行い地表の変動を把握することが行われている。
【０００３】
しかし、ＳＡＲの干渉解析から得られる情報は、地表面の変動前に地表に放射した電波の位相と地表面の変動後に地表に放射した電波の位相との差を算出した位相差情報（干渉縞画像）であり、判読が難しいという問題がある。
【０００４】
一方、特許文献１では、表示を行いたい輪郭データ（座標情報）を抽出し、抽出した輪郭データから、外周点を結ぶ線分データを補完する。補完した輪郭データから偏角を求め、求めた偏角から高速フーリエ変換により、周波数変換を行う。周波数に変換されたデータから、高周波成分を除去する。除去した周波数データから、逆高速フーリエ変換を行い、除去後の変換結果を輪郭データに戻す。戻した輪郭データから、ある一定の条件で間引きを行い、滑らかな輪郭データを表示させる方法が開示されている。
【０００５】
そこで、下記特許文献１に記載された方法により上記干渉縞画像を連続化処理して視認性を向上することが考えられる。
【０００６】
しかし、干渉解析で得られた位相差情報を高速フーリエ変換により周波数データに変換し、この周波数データに対してフィルタ処理を行い、フィルタ処理を行った周波数データを逆高速フーリエ変換により位相データに戻し、戻した位相データを３次元表示した場合、変動範囲が広い地表変動に適切な表現と、変動範囲が狭い地表変動（断層等の小さな変動）に適切な表現の両立が困難であった。
【先行技術文献】
【特許文献】
【０００７】
【特許文献１】特開２０１７−１８７４９８号公報
【発明の概要】
【発明が解決しようとする課題】
【０００８】
本発明の目的の一つは、異なる時刻にレーダ装置により取得された観測対象である地表の時系列データ（位相情報）をもとに干渉解析が行われた結果の位相差情報に含まれる変動範囲が広い地表変動と変動範囲が狭い地表変動の視認性向上を両立させることにある
【課題を解決するための手段】
【０００９】
上記目的を達成するために、本発明は、地表変動を表現するための地表変動表現装置であって、異なる時刻にレーダ装置により取得された観測対象である地表の時系列データ（位相情報）をもとに干渉解析が行われた結果の位相差情報を取得する位相差情報取得手段と、第１の窓サイズを有し、前記位相差情報取得手段が取得した位相差情報からノイズ除去を行う第１の位相差ノイズ除去手段と、前記第１の窓サイズより小さい第２の窓サイズを有し、前記位相差情報取得手段が取得した位相差情報からノイズ除去を行う第２の位相差ノイズ除去手段と、前記第２の位相差ノイズ除去手段から得られた第２の位相差情報から干渉縞画像を生成する干渉縞画像生成手段と、前記第１の位相差ノイズ除去手段から得られた第１の位相差情報と前記干渉縞画像との位置情報に基づき、前記第１の位相差ノイズ除去手段から得られた第１の位相差情報と前記干渉縞画像とを結合して表示する表示制御手段と、を備える地表変動表現装置を提供する。
【００１０】
上記表示制御手段は、前記第１の位相差ノイズ除去手段から得られた第１の位相差情報を３次元表示し、前記干渉縞画像を３次元表示した第１の位相差ノイズ除去手段から得られた第１の位相差情報に重畳して表示するのが好適である。
【００１１】
また、上記地表変動表現装置は、前記第１の位相差ノイズ除去手段から得られた第１の位相差情報を連続化して連続化位相差情報を生成する連続化手段をさらに備え、前記表示制御手段は、前記第１の位相差ノイズ除去手段から得られた第１の位相差情報の代わりに連続化位相差情報を使用するのが好適である。
【００１２】
また、本発明の他の態様は、地表変動を表現するための地表変動表現プログラムであって、コンピュータを、異なる時刻にレーダ装置により取得された観測対象である地表の時系列データをもとに干渉解析が行われた結果の位相差情報を取得する位相差情報取得手段、第１の窓サイズを有し、前記位相差情報取得手段が取得した位相差情報からノイズ除去を行う第１の位相差ノイズ除去手段、前記第１の窓サイズより小さい第２の窓サイズを有し、前記位相差情報取得手段が取得した位相差情報からノイズ除去を行う第２の位相差ノイズ除去手段、前記第２の位相差ノイズ除去手段から得られた第２の位相差情報から干渉縞画像を生成する干渉縞画像生成手段、前記第１の位相差ノイズ除去手段から得られた第１の位相差情報と前記干渉縞画像との位置情報に基づき、前記第１の位相差ノイズ除去手段から得られた第１の位相差情報と前記干渉縞画像とを結合して表示する表示制御手段、として機能させることを特徴とする地表変動表現プログラムを提供する。
【発明の効果】
【００１３】
本発明によれば、異なる時刻にレーダ装置により取得された観測対象である地表の時系列データ（位相情報）をもとに干渉解析が行われた結果の位相差情報に含まれる変動範囲が広い地表変動と変動範囲が狭い地表変動の視認性の向上を両立することができる。
【図面の簡単な説明】
【００１４】
【図１】実施形態にかかる地表変動表現装置の例の機能ブロック図である。
【図２】実施形態にかかる連続化部の動作の説明図である。
【図３】実施形態にかかる地表変動表現装置の動作例のフロー図である。
【図４】実施形態にかかる第１の位相差ノイズ除去部及び第２の位相差ノイズ除去部にてそれぞれノイズ除去を行った第１、第２の位相差情報を画像化した図である。
【図５】実施形態にかかる地表変動表現装置により地表面に現れた断層線を表現した図である。
【図６】実施形態にかかる地表変動表現装置により地表面に現れた河川部付近の断層線を表現した図である。
【発明を実施するための形態】
【００１５】
以下、本発明を実施するための形態（以下、実施形態という）を、図面に従って説明する。
【００１６】
図１には、実施形態にかかる地表変動表現装置１の例の機能ブロック図が示される。図１において、地表変動表現装置１は、位相差情報取得部１０、第１の位相差ノイズ除去部１２、第２の位相差ノイズ除去部１４、連続化部１６、干渉縞画像生成部１８、表示制御部２０、通信部２２、記憶部２４及びＣＰＵ２６を含んで構成されている。なお、ＣＰＵ２６以外にＧＰＵ等のアクセラレーターを用いてもよい。上記地表変動表現装置１は、ＣＰＵ２６、ＲＯＭ、ＲＡＭ、不揮発性メモリ、Ｉ／Ｏ、通信インターフェース等を備え、装置全体の制御及び各種演算を行うコンピュータとして構成されており、上記各機能は、例えばＣＰＵ２６とＣＰＵ２６の処理動作を制御するプログラムとにより実現される。
【００１７】
位相差情報取得部１０は、異なる時刻にレーダ装置により取得された観測対象である地表の時系列データ（位相情報）をもとに干渉解析が行われた結果の位相差情報を取得する。ここで、レーダ装置としては、例えばＳＡＲ（Synthetic Aperture Radar 合成開口レーダ）が挙げられるが、これには限定されない。また、位相情報は、上記レーダ装置により取得された、異なる時刻に地表の同一点から反射された電波の反射波を含んだ情報である。位相差情報取得部１０は、上記ＳＡＲ等により観測した時系列データの位相情報から干渉解析等が行われて得られた位相差情報を、通信部２２を介して取得することができる。取得した位相差情報は、記憶部２４に記憶させる。位相差情報には、合成開口レーダからの視線方向の距離変化の情報が含まれており、これを使用して地表面の高さの変化の情報を抽出することができる。
【００１８】
第１の位相差ノイズ除去部１２は、上記位相差情報取得部１０が取得した位相差情報を記憶部２４から読み出し、第１の窓サイズを有するノイズ除去フィルタによりノイズ除去を行う。ノイズ除去を行った第１の位相差情報は、記憶部２４に記憶させる。ここで、上記窓サイズとは、処理対象である中心点の位相差情報と中心点の周囲を囲む点の位相差情報とを合わせた位相差情報のデータ数である。この窓サイズが大きい場合は、広範囲に地表面が同じ変動をした断層線などの現象をノイズ除去後の位相差情報に残すことができるが、地表面の狭い領域における変動が残り難くなる。一方、窓サイズが小さい場合は、地表面の狭い領域における変動をノイズ除去後の位相差情報に残すことができるが、ノイズが多くなると共に、広範囲に地表面が同じ変動を示した位相差情報が捉え難くなる。
【００１９】
第２の位相差ノイズ除去部１４は、上記位相差情報取得部１０が取得した位相差情報を記憶部２４から読み出し、第１の窓サイズより小さい第２の窓サイズを有するノイズ除去フィルタによりノイズ除去を行う。ノイズ除去を行った第２の位相差情報は、記憶部２４に記憶させる。
【００２０】
以上に述べた第１の位相差ノイズ除去部１２及び第２の位相差ノイズ除去部１４を構成するノイズ除去フィルタとしては、例えばGoldstein & Wernerのフィルタ等を使用できるが、これには限定されない。
【００２１】
連続化部１６は、上記第１の位相差情報を記憶部２４から読み出し、これを連続化して連続化位相差情報を生成する。生成した連続化位相差情報は、記憶部２４に記憶させる。連続化部１６の動作は後述する。
【００２２】
干渉縞画像生成部１８は、上記第２の位相差情報を記憶部２４から読み出し、これから干渉縞画像を生成する。干渉縞画像生成部１８は、記憶部２４に記憶された第２の位相差情報に基づいて干渉縞を生成し、干渉縞画像とする。干渉縞画像生成部１８が生成した干渉縞画像は、記憶部２４に記憶させる。
【００２３】
表示制御部２０は、上記第１の位相差情報と上記干渉縞画像とを記憶部２４から読み出し、これらの位置情報に基づいて上記第１の位相差情報と上記干渉縞画像とを結合し、液晶表示装置その他の適宜な表示装置を制御して表示する。なお、連続化部１６が生成した連続化位相差情報がある場合には、上記第１の位相差情報の代わりに連続化位相差情報を使用する。
【００２４】
通信部２２は、適宜なインターフェースにより構成され、ＳＡＲ衛星から時系列データ（位相情報）を受信して、これから干渉解析を行い位相差情報を抽出する外部のサーバ等から通信部２２が位相差情報を取得して位相差情報取得部１０に渡す構成とする。
【００２５】
記憶部２４は、ハードディスク装置、ソリッドステートドライブ（ＳＳＤ）等の不揮発性メモリで構成され、上記各種情報等、及びＣＰＵ２６の動作プログラム等の、地表変動表現装置１が行う各処理に必要な情報を記憶する。なお、記憶部２４としては、デジタル・バーサタイル・ディスク（ＤＶＤ）、コンパクトディスク（ＣＤ）、光磁気ディスク（ＭＯ）、フレキシブルディスク（ＦＤ）、磁気テープ、電気的消去および書き換え可能な読出し専用メモリ（ＥＥＰＲＯＭ）、フラッシュ・メモリ等を使用してもよい。また、記憶部２４には、主としてＣＰＵ２６の作業領域として機能するランダムアクセスメモリ（ＲＡＭ）、及びＢＩＯＳ等の制御プログラムその他のＣＰＵ２６が使用するデータが格納される読み出し専用メモリ（ＲＯＭ）を含めるのが好適である。
【００２６】
図２(ａ)、（ｂ）には、連続化部１６の動作の説明図が示される。図２(ａ)には、異なる時刻に取得された位相情報から干渉解析を行った位相差の変化の例が示される。図２(ａ)において、縦軸が位相差であり、横軸が地表面上の位置を表す。
【００２７】
一般に、上記異なる時刻に取得された位相情報から干渉解析を行って得た位相差は、使用する電波の半波長（測定用の電波が往復するので）分の地表変動毎に−πからπまで変化し、地表変動が半波長を超える度にその値がπから−πに戻って値の変化が繰り返されるので、図２(ａ)に示されるように、位相差の変化はのこぎり状の波形となる。従って、位相差の変化をこのまま画像化すると、πから−πに戻る部分で色相、輝度、３次元表示したときの高さ等が急激に変化し、視認性が低い画像となる。
【００２８】
そこで、連続化部１６は、のこぎり状の関数値（−πからπの値）を三角関数（例えばｓｉｎ関数）に代入して正弦波に変換することにより、図２(ｂ)に示されるように、のこぎり波形を正弦波波形に連続化する。
【００２９】
図３には、実施形態にかかる地表変動表現装置１の動作例のフロー図が示される。図３において、位相差情報取得部１０が、異なる時刻にレーダ装置により取得された観測対象である地表の時系列データ（位相情報）をもとに干渉解析が行われた結果の位相差情報を取得して記憶部２４に記憶させると（Ｓ１）、第１の位相差ノイズ除去部１２が位相差情報を記憶部２４から読み出し、第１の窓サイズを有するノイズ除去フィルタによりノイズ除去を行う（Ｓ２）。また、第２の位相差ノイズ除去部１４が位相差情報を記憶部２４から読み出し、第１の窓サイズより小さい第２の窓サイズを有するノイズ除去フィルタによりノイズ除去を行う（Ｓ３）。上記Ｓ２及びＳ３でノイズ除去された第１、第２の位相差情報は、記憶部２４に記憶させる。
【００３０】
連続化部１６は、上記第１の位相差情報を記憶部２４から読み出し、これを連続化して連続化位相差情報を生成する（Ｓ４）。生成した連続化位相差情報は、記憶部２４に記憶させる。なお、Ｓ４の処理は省略してもよい。
【００３１】
また、干渉縞画像生成部１８は、上記第２の位相差情報を記憶部２４から読み出し、これから干渉縞画像を生成する（Ｓ５）。
【００３２】
表示制御部２０は、上記第１の位相差情報と上記干渉縞画像とを記憶部２４から読み出し、これらの位置情報に基づいて上記第１の位相差情報と上記干渉縞画像とを結合し（Ｓ６）、液晶表示装置その他の適宜な表示装置を制御して表示する（Ｓ７）。なお、Ｓ４の処理が実行されている場合には、第１の位相差情報の代わりに連続化位相差情報を使用する。
【００３３】
上記Ｓ６の処理において、第１の位相差情報または連続化位相差情報と干渉縞画像とを結合する際には、第１の位相差情報または連続化位相差情報から各点の高さ情報を抽出し、この高さ情報を干渉縞画像の各画素の中心点の高さ情報として設定して３次元画像を生成する。あるいは、上記高さ情報に基づき、第１の位相差情報または連続化位相差情報を３次元画像とし、これにテクスチャとして第２の位相差情報から作成した干渉縞画像を重畳して３次元画像を生成してもよい。これらの処理を実行する場合、第１の位相差情報または連続化位相差情報の有する位置情報と、干渉縞画像の有する位置情報とに基づき、同じ位置情報を有する点と画素の中心点同士で上記高さ情報の設定あるいはテクスチャとしての重畳を行う。第１の位相差情報または連続化位相差情報と干渉縞画像とは、同じ位相差情報から生成するものであり、相互の有する位置情報は一致するので、上述するように、位置情報に基づく結合が可能となる。
【００３４】
なお、上記結合処理に使用する高さ情報として、第１の位相差情報または連続化位相差情報から抽出した高さ情報を使用するのは、第２の位相差情報のように窓サイズが小さなフィルタによりノイズ除去処理を行った位相差情報から抽出した高さ情報を用いると、残ノイズ等に起因する高さ変化の影響を受けて、見にくい画像になるからである。このため、本実施形態では、第２の窓サイズより大きい第１の窓サイズを有するノイズ除去フィルタを使用する第１の位相差ノイズ除去部１２が出力した第１の位相差情報またはこれを連続化部１６により連続化した連続化位相差情報から高さ情報を抽出している。
【００３５】
図４（ａ）、（ｂ）には、第１の位相差ノイズ除去部１２及び第２の位相差ノイズ除去部１４にてそれぞれノイズ除去を行った第１、第２の位相差情報を干渉縞画像生成部１８が干渉縞画像として画像化した図が示される。図４（ａ）が第１の位相差ノイズ除去部１２を使用した場合であり、図４（ｂ）が第２の位相差ノイズ除去部１４を使用した場合である。
【００３６】
図４（ａ）の例は、第１の窓サイズとして、中心の点とその周囲合わせて縦方向に５１２点の位相差情報と横方向に５１２点（５１２×５１２点）の位相差情報を読み込むノイズ除去フィルタを使用してノイズ除去を行い、上記中心点の位相差情報としたものであり、全ての位相差情報（処理対象の地表面の全点）に対して行われる。また、図４（ｂ）の例は、第１の窓サイズより小さい第２の窓サイズとして、中心の点とその周囲合わせて縦方向に３２点の位相差情報と横方向に３２点（３２×３２点）の位相差情報を読み込むノイズ除去フィルタを使用してノイズ除去を行い、上記中心点の位相差情報としたものであり、全ての位相差情報（処理対象の地表面の全点）に対して行われる。なお、窓サイズは、上記値に限定されるものではなく、処理する位相差情報、得たい地表変動の画像等に応じて適宜決定することができる。
【００３７】
また、図４（ａ）、（ｂ）の例で示された図には、縞模様が表現されているが、縞模様を構成する曲線の位置が、上記位相差情報の位相差がπから−πに変化する位置であり、この曲線の間で使用する電波の半波長分の地表変動が発生していることを表している。
【００３８】
窓サイズが大きいノイズ除去フィルタを使用した第１の位相差ノイズ除去部１２による処理結果は、広範囲の領域の地表変動を捉える一方で、比較的範囲が狭い領域の地表変動はつぶれてしまっている。一方、窓サイズが小さいノイズ除去フィルタを使用した第２の位相差ノイズ除去部１４による処理結果は、比較的範囲が狭い領域の地表変動を捉える一方で、広範囲の領域の地表変動が、範囲が狭い領域の地表変動に隠れてしまい、広範囲の領域の地表変動が確認しにくくなっている。また、第２の位相差ノイズ除去部１４による処理結果には、ノイズ（図に現れる不規則な斑点状の部分）も多くなり、地表面の視認性が低下している。
【００３９】
そこで、本実施形態にかかる地表変動表現装置１では、表示制御部２０が、上記第１の位相差情報または上記連続化位相差情報と、上記干渉縞画像とを結合し、表示することにより、地表面の大きな変動を表示するとともに、画像情報に含まれる小さな変動の視認性を向上することができる。
【００４０】
図５（ａ）、（ｂ）、（ｃ）には、実施形態にかかる地表変動表現装置１により地表面に現れた断層線を表現した図が示される。図５（ａ）は、上記第１の位相差情報を干渉縞画像生成部１８が干渉縞画像として画像化した図であり、図５（ｂ）は、干渉縞画像生成部１８が第２の位相差情報から生成した干渉縞画像を表示制御部２０が表示した図であり、図５（ｃ）は、表示制御部２０が、上記第１の位相差情報を連続化部１６が連続化した連続化位相差情報と干渉縞画像とを結合して表示した図である。
【００４１】
なお、図５（ａ）、（ｂ）において、断層線が発生した領域が破線で、図５（ｃ）では実線で囲まれて示されている。
【００４２】
図５（ａ）、（ｂ）の図は、２次元表示であるのに対して、図５（ｃ）の図は、３次元表示となっている。これは、上述したように、図５（ａ）の図から高さ情報を抽出し、図５（ｂ）の干渉縞画像と結合して得た３次元画像となっている。これにより、干渉ＳＡＲ等により観察した地表面を３次元表示することができ、より見やすい表示とすることができる。
【００４３】
また、窓サイズが小さいノイズ除去フィルタを使用することにより、狭い領域の変化を見やすくした干渉縞画像（図５（ｂ））を、図５（ａ）と結合することにより、断層線のような小さな地表面の変化を見やすくでき、断層線等の空間スケールの大きな変動（広範囲の領域の地表変動）に含まれる局所的な変動の視認性を向上できる。
【００４４】
図６（ａ）、（ｂ）、（ｃ）には、実施形態にかかる地表変動表現装置１により地表面に現れた河川部付近の断層線を表現した図が示される。これらの図も、図５（ａ）、（ｂ）、（ｃ）と同様に、図６（ａ）は、上記第１の位相差情報を干渉縞画像生成部１８が干渉縞画像として画像化した図であり、図６（ｂ）は、干渉縞画像生成部１８が第２の位相差情報から生成した干渉縞画像を表示制御部２０が表示した図であり、図６（ｃ）は、表示制御部２０が、上記連続化位相差情報と干渉縞画像とを結合して表示した図である。
【００４５】
なお、図６（ａ）、（ｂ）、（ｃ）において、断層線が発生した領域が破線で囲まれて示されている。
【００４６】
図６（ｃ）の例でも、３次元表示したこと、窓サイズが小さいノイズ除去フィルタを使用することにより狭い領域の変化を見やすくしたことにより、断層線で発生した細かな変動の視認性も向上することができる。
【符号の説明】
【００４７】
１  地表変動表現装置、１０  位相差情報取得部、１２  第１の位相差ノイズ除去部、１４  第２の位相差ノイズ除去部、１６  連続化部、１８  干渉縞画像生成部、２０  表示制御部、２２  通信部、２４  記憶部、２８  ＣＰＵ。
【図１】
2020106460.tif 000003
【図２】
2020106460.tif 000004
【図３】
2020106460.tif 000005
【図４】
2020106460.tif 000006
【図５】
2020106460.tif 000007
【図６】
2020106460.tif 000008