LiCSBASで計算したSBAS変形結果の時系列曲線グラフのエクスポート方法

LiCSBAS (SBAS を使用した大規模 InSAR 時系列解析) は、合成開口レーダー干渉法 (InSAR) データを処理するための時系列解析手法です。地殻変動、地震、火山活動など、地表の地質変化を監視するように設計されています。LiCSBAS の主な機能と概要は次のとおりです。

時系列解析: LiCSBAS は、複数の InSAR 観測瞬間からのデータを処理して、地表変位の長期および短期変化を検出および解析できる時系列解析手法です。
大規模モニタリング: LiCSBAS は大規模エリアのモニタリングに適しており、都市から国、さらには地球規模に至るまで、広範囲の地理的エリアをカバーできます。
累積変位: LiCSBAS は、表面の累積変位を計算して視覚化します。これは、表面変化の累積効果であり、地殻運動や地震の研究に一般的に使用されます。
速度場推定: LiCSBAS は、表面上の速度場を推定し、表面のさまざまな領域での運動の速度と方向を明らかにします。
ノイズフィルタリング: LiCSBAS には、変位推定の精度を向上させるために、InSAR データに対するノイズフィルタリングステップが含まれています。
マスク処理: マスク処理をサポートしています。これにより、無関係な表面領域を除外して、分析対象の領域に焦点を当てることができます。
グラフィカルインターフェイス: 一般に、LiCSBAS ツールは使いやすいグラフィカルユーザーインターフェイスを備えており、ユーザーは対話的に分析領域を選択し、速度フィールド、累積変位マップや時系列などを描画できます。
科学研究と応用: LiCSBAS は、地質学、地球物理学、自然災害モニタリング、および地表の地質学的および環境変化を研究するその他の分野で広く使用されています。

要約すると、LiCSBAS は、InSAR データを処理および分析して、地表の変化や動きを理解するための強力なツールです。これは、科学研究および応用分野、特に地質学的危険の監視や地殻変動の研究で広く使用されています。

変形結果をエクスポートするには、スクリプトLiCSBAS_plot_ts.pyが必要です。このスクリプトは、速度、累積変位、ノイズインジケーターを表示し、変位の時系列図を描画します。対話的に表示画像・領域を変更したり、時系列プロットを描画する点を選択したりできます。右クリックしてドラッグすることで参照領域を変更することもできます。

コマンドの説明:

LiCSBAS_plot_ts.py [-i cum[_filt].h5] [--i2 cum*.h5] [-m yyyymmdd] [-d results_dir]
    [-u U.geo] [-r x1:x2/y1:y2] [--ref_geo lon1/lon2/lat1/lat2] [-p x/y]
    [--p_geo lon/lat] [-c cmap] [--nomask] [--vmin float] [--vmax float]
    [--auto_crange float] [--dmin float] [--dmax float] [--ylen float]
    [--ts_png pngfile]

 -i    Input cum hdf5 file (Default: ./cum_filt.h5 or ./cum.h5)
 --i2  Input 2nd cum hdf5 file
       (Default: cum.h5 if -i cum_filt.h5, otherwise none)
 -m    Refereference (master) date for time series (Default: first date)
 -d    Directory containing noise indices (e.g., mask, coh_avg, etc.)
       (Default: "results" at the same dir as cum[_filt].h5)
 -u    Input U.geo file to show incidence angle (Default: ../GEOCml*/U.geo)
 -r    Initial reference area (Default: same as info/*ref.txt)
       0 for x2/y2 means all. (i.e., 0:0/0:0 means whole area).
 --ref_geo   Initial reference area in geographical coordinates.
 -p    Initial selected point for time series plot (Default: ref point)
 --p_geo     Initial selected point in geographical coordinates.
 -c    Color map for velocity and cumulative displacement.
       See help of LiCSBAS_disp_img.py.
       (Default: SCM.roma_r, reverse of SCM.roma)
 --nomask     Not use mask (Default: use mask)
 --vmin|vmax  Min|max values of color for velocity map (Default: auto)
 --dmin|dmax  Min|max values of color for cumulative displacement map
              (Default: auto)
 --auto_crange  Percentage of color range used for automatic determination
              (Default: 99 %)
 --ylen       Y Length of time series plot in mm (Default: auto)
 --ts_png     Output png file of time series plot (not display interactive viewers)

コマンド翻訳:

-i：输入的累积位移数据的HDF5文件。默认情况下，如果存在cum_filt.h5则使用它，否则使用cum.h5。

--i2：第二个输入的累积位移数据的HDF5文件（可选）。默认情况下，如果使用了-i参数指定了cum_filt.h5，则默认不使用第二个输入文件。

-m：用于时间序列的参考（主）日期，以便在时间序列中进行相对分析。默认情况下使用第一个日期。

-d：包含噪音指标数据（如掩模、相干度等）的目录路径。默认情况下，这些数据应该位于与cum[_filt].h5相同的目录下的results目录中。

-u：输入的U.geo文件路径，用于显示入射角度信息。默认情况下，该文件应该位于../GEOCml*/U.geo路径下。

-r：初始的参考区域的定义。默认情况下，使用info/*ref.txt中定义的区域。-r参数的格式是x1:x2/y1:y2，其中0表示整个区域（例如，0:0/0:0表示整个区域）。

--ref_geo：初始参考区域的地理坐标定义。

-p：初始选择的用于时间序列绘制的点。默认情况下，使用参考点。

--p_geo：初始选择的时间序列绘制点的地理坐标。

-c：用于速度和累积位移地图的颜色映射。可以参考LiCSBAS_disp_img.py的帮助文档以了解可用的选项。默认使用SCM.roma_r颜色映射。

--nomask：不使用掩模。默认情况下，会使用掩模来排除无效数据。

--vmin和--vmax：速度地图的颜色映射的最小和最大值。默认情况下，自动确定。

--dmin和--dmax：累积位移地图的颜色映射的最小和最大值。默认情况下，自动确定。

--auto_crange：自动确定颜色范围的百分比。默认为99％。

--ylen：时间序列绘制的Y轴长度（以毫米为单位）。默认情况下，自动确定。

--ts_png：时间序列绘制的输出PNG文件的路径。如果指定此参数，将生成PNG文件，否则将在交互式查看器中显示。

操作例:
コマンド:LiCSBAS_plot_ts.py -i TS_GEOCml1maskclip/cum_filt.h5 -r 480:489/400:411 --ts_png TS_GEOCml1maskclip

結果の生成:
ここに画像の説明を挿入します
注意:
元の出力結果の解像度は低いため、次のディレクトリファイルを変更する必要があります。

LiCSBAS/bin/LiCSBAS_plot_ts.py

コードの最後に dpi=600 を追加します。数値はエクスポートされる画像の解像度です。
ここに画像の説明を挿入します

LiCSBASで計算したSBAS変形結果の時系列曲線グラフのエクスポート方法

おすすめ