#!/usr/bin/env python3 ############################################################ # Program is part of PlotData # # Author: Falk Amelung, Dec 2023 # ############################################################ import os import sys import re import argparse from plotdata import logging_function ############################################################ EXAMPLE = """example: plot_data.py GalapagosSenDT128 plot_data.py GalapagosSenDT128/mintpy plot_data.py GalapagosSenDT128/mintpy/S1_IW12_128_0593_0597_20181005_XXXXXXXX.he5 plot_data.py GalapagosSenDT128/mintpy --plot-type=velocity --subset-lalo=-0.52:-0.28,-91.7:-91.4 --period=20200131-20231231 --gps --seismicity plot_data.py GalapagosSenDT128/mintpy GalapagosSenAT106/mintpy_orig --plot-type=horzvert --subset-lalo=-1.0:-0.75,-91.55:-91.25 --period=20220101-20230831 --vlim -5 5 plot_data.py MaunaLoaSenDT87/mintpy/ plot_data.py MaunaLoaSenDT87 --plot-type ifgram --seismicity --gps plot_data.py MaunaLoaSenDT87 --plot-type shaded_relief --seismicity --gps plot_data.py MaunaLoaSenDT87 --plot-type velocity --seismicity --gps plot_data.py MaunaLoaSenAT124 MaunaLoaSenDT87 --ref-lalo 19.55,-155.45 plot_data.py MaunaLoaSenDT87/mintpy_5_20 --plot-type velocity plot_data.py MaunaLoaSenDT87/mintpy_5_20 MaunaLoaSenAT124/mintpy_5_20 --plot-type velocity --ref-lalo 19.495,-155.555 --period 20181001-20221122 --subset-lalo 19.43:19.5,-155.62:-155.55 --vlim -5 5 plot_data.py MaunaLoaSenDT87/mintpy_5_20 MaunaLoaSenAT124/mintpy_5_20 --plot-type horzvert --ref-lalo 19.495,-155.555 --period 20181001-20221122 --subset-lalo 19.43:19.5,-155.62:-155.55 --vlim -5 5 plot_data.py MaunaLoaSenDT87/mintpy_5_20 --plot-type shaded-relief --gps --period 20181001-20221122 --dem-file $SCRATCHDIR/MaunaLoa/MLtry/data/demGeo.h5 plot_data.py MaunaLoaSenDT87/mintpy_5_20 --plot-type shaded-relief --gps --gps-scale-fac 200 --gps-key-length 1 plot_data.py MaunaLoaSenDT87/mintpy_5_20 --plot-type shaded-relief --gps --seismicity plot_data.py MaunaLoaSenDT87/mintpy_5_20 --plot-type shaded-relief --subset-lalo 19.43:19.5,-155.62:-155.55 --seismicity plot_data.py MaunaLoaSenDT87/mintpy_5_20 MaunaLoaSenAT124/mintpy_5_20 --plot-type velocity --ref-lalo 19.55,-155.45 --period 20220801-20221127 --vlim -20 20 --save-gbis --gps --seismicity --fontsize 14 plot_data.py GalapagosSenDT128/mintpy --plot-type=velocity --subset-lalo=-0.52:-0.28,-91.7:-91.4 --period=20200131-20221231 --gps --seismicity plot_data.py GalapagosSenDT128/mintpy --plot-type=velocity --subset-lalo=-0.52:-0.28,-91.7:-91.4 --period=20200131-20220430 plot_data.py GalapagosSenDT128/mintpy --plot-type=velocity --period=20200131-20220430 plot_data.py GalapagosSenDT128/mintpy --plot-type=velocity --subset-lalo=-0.52:-0.28,-91.7:-91.4 plot_data.py GalapagosSenDT128/mintpy --subset-lalo=-0.86:-0.77:-91.19:-91.07 --ref-lalo=-0.771,-91.19 """ def create_parser(): synopsis = 'Plotting of InSAR, GPS and Seismicity data' epilog = EXAMPLE parser = argparse.ArgumentParser(description=synopsis, epilog=epilog, formatter_class=argparse.RawTextHelpFormatter) parser.add_argument('data_dir', nargs='*', help='Directory(s) with InSAR data.\n') parser.add_argument('--subset-lalo', nargs='?', dest='plot_box', type=str, default=None, help='geographic area plotted') parser.add_argument('--period', dest='period', metavar='YYYYMMDD-YYYYMMDD', default=None, help='time period (Default: full time period)') parser.add_argument('--seismicity', dest='flag_seismicity', action='store_true', default=False, help='flag to add seismicity') parser.add_argument('--gps', dest='flag_gps', action='store_true', default=False, help='flag to add GPS vectors') parser.add_argument('--plot-type', dest='plot_type', default='velocity', help='Type of plot: velocity, horzvert, ifgram, step, shaded_relief (Default: velocity).') parser.add_argument('--dem', dest='dem_file', default=None, help='external DEM file (Default: geo/geo_geometryRadar.h5)') parser.add_argument('--lines', dest='line_file', default=None, help='fault file (Default: None, but plotdata/data/hawaii_lines_new.mat for Hawaii)') parser.add_argument('--gps-scale-fac', dest='gps_scale_fac', default=500, type=int, help='GPS scale factor (Default: 500)') parser.add_argument('--gps-key-length', dest='gps_key_length', default=4, type=int, help='GPS key length (Default: 4)') parser.add_argument('--gps-units', dest='gps_unit', default="cm", help='GPS units (Default: cm)') parser.add_argument('--unit', dest='unit', default="cm", help='InSAR units (Default: cm)') parser.add_argument('--fontsize', dest='font_size', default=12, type=int, help='fontsize for view.py (Default: 12)') parser.add_argument('--ref-lalo', nargs='*', metavar=('LAT,LON or LAT LON'), type=str, default=None, help='reference point (default: existing reference point)') parser.add_argument("--lalo", nargs='*', metavar=('LAT,LON or LAT LON or LAT1,LON1 LAT2,LON2'), type=str, default=None, help="lat/lon coords of pixel for timeseries (default: ?)") parser.add_argument('--mask-thresh', dest='mask_vmin', type=float, default=0.7, help='coherence threshold for masking (Default: 0.7)') parser.add_argument('--vlim', dest='vlim', nargs=2, metavar=('VMIN', 'VMAX'), type=float, help='colorlimit') parser.add_argument("--noreference", dest="show_reference_point", action='store_false', default=True, help="hide reference point (default: False)" ) parser.add_argument('--save-gbis', dest='flag_save_gbis', action='store_true', default=False, help='save GBIS files') parser.add_argument('--style', dest='style', choices={'image', 'scatter'}, default='image', help='Plot data as image or scatter (default: %(default)s).') parser.add_argument('--scatter-size', dest='scatter_marker_size', type=float, metavar='SIZE', default=10, help='Scatter marker size in points**2 (default: %(default)s).') parser.add_argument('--shade-exag ',dest='shade_exag', type=float, default=1, help='Shade exaggeration factor (Default: 0.5)') inps = parser.parse_args() if len(inps.data_dir) < 1 or len(inps.data_dir) > 2: parser.error('USER ERROR: You must provide 1 or 2 directory paths.') if inps.plot_box: inps.plot_box = [float(val) for val in inps.plot_box.replace(':', ',').split(',')] # converts to plot_box=[19.3, 19.6, -155.8, -155.4] if inps.ref_lalo: inps.ref_lalo = parse_lalo(inps.ref_lalo) if inps.lalo: inps.lalo = parse_lalo(inps.lalo) if inps.dem_file and '$' in inps.dem_file: inps.dem_file = os.path.expandvars(inps.dem_file) # Hardwire line file for Hawaii data if ('Hawaii' in inps.data_dir or 'Mauna' in inps.data_dir or 'Kilauea' in inps.data_dir): inps.line_file = os.getenv('RSMASINSAR_HOME') + '/tools/PlotData' + '/data/hawaii_lines_new.mat' return inps def parse_lalo(str_lalo): """Parse the lat/lon input from the command line.""" if ',' in str_lalo[0]: lalo = [[float(coord) for coord in pair.split(',')] for pair in str_lalo] else: lalo = [[float(str_lalo[i]), float(str_lalo[i+1])] for i in range(0, len(str_lalo), 2)] if len(lalo) == 1: # if given as one string containing ',' lalo = lalo[0] return lalo ############################################################ def main(iargs=None): logging_function.log(os.getcwd(), os.path.basename(__file__) + ' ' + ' '.join(sys.argv[1:])) # parse inps = create_parser() # import from plotdata.prepare_and_plot import run_prepare from plotdata.prepare_and_plot import run_plot # Goto $SCRATHDIR and run os.chdir(os.getenv('SCRATCHDIR')) data_dict = run_prepare(inps) run_plot(data_dict, inps) return ############################################################ if __name__ == '__main__': main(iargs=sys.argv)