LiCSBAS计算的SBAS形变结果导出时序曲线图方法

LiCSBAS（Large-scale InSAR Time-Series Analysis using SBAS）是一种用于处理合成孔径雷达干涉测量（InSAR）数据的时间序列分析方法。它旨在监测地球表面的地质变化，如地壳运动、地震、火山活动等。以下是LiCSBAS的主要特点和介绍：

1. 时间序列分析：LiCSBAS是一种时间序列分析方法，它能够处理多个InSAR观测时刻的数据，以检测和分析地表位移的长期和短期变化。

2. 大范围监测：LiCSBAS适用于大范围地区的监测，可以覆盖广泛的地理区域，从城市到国家甚至全球尺度。

3. 累积位移：LiCSBAS计算并可视化地表的累积位移，这是地表变化的积累效应，通常用于研究地壳运动和地震。

4. 速度场估计：LiCSBAS能够估计地表的速度场，揭示了地表不同区域的运动速度和方向。

5. 噪音过滤：LiCSBAS包括对InSAR数据的噪音过滤步骤，以提高位移估计的准确性。

6. 掩模处理：它支持掩模处理，可以排除不相关的地表区域，以聚焦在感兴趣的地区进行分析。

7. 图形化界面：通常，LiCSBAS工具具有用户友好的图形用户界面，允许用户交互式地选择分析区域、绘制速度场、累积位移图和时间序列等。

8. 科学研究和应用：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，数字为导出的图片分辨率大小。