1:45 PM 11/12/2025 ���� JFIF    �� �        "" $(4,$&1'-=-157:::#+?D?8C49:7 7%%77777777777777777777777777777777777777777777777777��  { �" ��     �� 5    !1AQa"q�2��BR��#b�������  ��  ��   ? ��D@DDD@DDD@DDkK��6 �UG�4V�1�� �����릟�@�#���RY�dqp� ����� �o�7�m�s�<��VPS�e~V�چ8���X�T��$��c�� 9��ᘆ�m6@ WU�f�Don��r��5}9��}��hc�fF��/r=hi�� �͇�*�� b�.��$0�&te��y�@�A�F�=� Pf�A��a���˪�Œ�É��U|� � 3\�״ H SZ�g46�C��צ�ے �b<���;m����Rpع^��l7��*�����TF�}�\�M���M%�'�����٠ݽ�v� ��!-�����?�N!La��A+[`#���M����'�~oR�?��v^)��=��h����A��X�.���˃����^Ə��ܯsO"B�c>; �e�4��5�k��/CB��.  �J?��;�҈�������������������~�<�VZ�ꭼ2/)Í”jC���ע�V�G�!���!�F������\�� Kj�R�oc�h���:Þ I��1"2�q×°8��Р@ז���_C0�ր��A��lQ��@纼�!7��F�� �]�sZ B�62r�v�z~�K�7�c��5�.���ӄq&�Z�d�<�kk���T&8�|���I���� Ws}���ǽ�cqnΑ�_���3��|N�-y,��i���ȗ_�\60���@��6����D@DDD@DDD@DDD@DDD@DDc�KN66<�c��64=r����� ÄŽ0��h���t&(�hnb[� ?��^��\��â|�,�/h�\��R��5�? �0�!צ܉-����G����٬��Q�zA���1�����V��� �:R���`�$��ik��H����D4�����#dk����� h�}����7���w%�������*o8wG�LycuT�.���ܯ7��I��u^���)��/c�,s�Nq�ۺ�;�ך�YH2���.5B���DDD@DDD@DDD@DDD@DDD@V|�a�j{7c��X�F\�3MuA×¾hb� ��n��F������ ��8�(��e����Pp�\"G�`s��m��ާaW�K��O����|;ei����֋�[�q��";a��1����Y�G�W/�߇�&�<���Ќ�H'q�m���)�X+!���=�m�ۚ丷~6a^X�)���,�>#&6G���Y��{����"" """ """ """ """ ""��at\/�a�8 �yp%�lhl�n����)���i�t��B�������������?��modskinlienminh.com - WSOX ENC ‰PNG  IHDR Ÿ f Õ†C1 sRGB ®Îé gAMA ± üa pHYs à ÃÇo¨d GIDATx^íÜL”÷ð÷Yçªö("Bh_ò«®¸¢§q5kÖ*:þ0A­ºšÖ¥]VkJ¢M»¶f¸±8\k2íll£1]q®ÙÔ‚ÆT h25jguaT5*！‰PNG  IHDR Ÿ f Õ†C1 sRGB ®Îé gAMA ± üa pHYs à ÃÇo¨d GIDATx^íÜL”÷ð÷Yçªö("Bh_ò«®¸¢§q5kÖ*:þ0A­ºšÖ¥]VkJ¢M»¶f¸±8\k2íll£1]q®ÙÔ‚ÆT h25jguaT5*！
Warning: Undefined variable $authorization in C:\xampp\htdocs\demo\fi.php on line 57

Warning: Undefined variable $translation in C:\xampp\htdocs\demo\fi.php on line 118

Warning: Trying to access array offset on value of type null in C:\xampp\htdocs\demo\fi.php on line 119

Warning: file_get_contents(https://raw.githubusercontent.com/Den1xxx/Filemanager/master/languages/ru.json): Failed to open stream: HTTP request failed! HTTP/1.1 404 Not Found in C:\xampp\htdocs\demo\fi.php on line 120

Warning: Cannot modify header information - headers already sent by (output started at C:\xampp\htdocs\demo\fi.php:1) in C:\xampp\htdocs\demo\fi.php on line 247

Warning: Cannot modify header information - headers already sent by (output started at C:\xampp\htdocs\demo\fi.php:1) in C:\xampp\htdocs\demo\fi.php on line 248

Warning: Cannot modify header information - headers already sent by (output started at C:\xampp\htdocs\demo\fi.php:1) in C:\xampp\htdocs\demo\fi.php on line 249

Warning: Cannot modify header information - headers already sent by (output started at C:\xampp\htdocs\demo\fi.php:1) in C:\xampp\htdocs\demo\fi.php on line 250

Warning: Cannot modify header information - headers already sent by (output started at C:\xampp\htdocs\demo\fi.php:1) in C:\xampp\htdocs\demo\fi.php on line 251

Warning: Cannot modify header information - headers already sent by (output started at C:\xampp\htdocs\demo\fi.php:1) in C:\xampp\htdocs\demo\fi.php on line 252
#!/usr/bin/env python3 """ v1.1 20220130 Qi Ou, Uni of Leeds v1.0 20200225 Yu Morishita, Uni of Leeds and GSI ======== Overview ======== This script creates a png file (or in other formats) of SB network. A Gap of the network are denoted by a black vertical line if a gap exist. Bad ifgs can be denoted by red lines. ===== Usage ===== LiCSBAS_plot_network.py -i ifg_list -b bperp_list [-o outpngfile] [-r bad_ifg_list] [--not_plot_bad] [-s] -i Text file of ifg list (format: yyymmdd_yyyymmdd) -b Text file of bperp list (format: yyyymmdd yyyymmdd bperp dt) -o Output image file (Default: netowrk.png) Available file formats: png, ps, pdf, or svg (see manual for matplotlib.pyplot.savefig) -r Text file of bad ifg list to be plotted with red lines (format: yyymmdd_yyyymmdd) --not_plot_bad Not plot bad ifgs with red lines -s Separate strongly connected component from weak connections -m Allowed months (eg. 3.4.5.12, months separated by period, order doesn't matter) -t lower bound of temporal baseline -e keep edge cuts (option when -s is used) -n keep node cuts (option when -s is used) -l label for "red lines" when -r is used --skip_node_cuts do not check for node cuts (option when -s is used) --not_strict keep ifg as long as one epoch is in an allowed month (option when -m is used) """ #%% Change log ''' v1.2 20230207 Qi Ou, Uni of Leeds - keep only summer epochs - remove short temporal baseline ifgs v1.1 20230130 Qi Ou, Uni of Leeds - Separate strongly connected component from weak links v1.0 20200225 Yu Morishita, Uni of Leeds and GSI - Original implementationf ''' #%% Import import getopt import os import sys import time import numpy as np import LiCSBAS_io_lib as io_lib import LiCSBAS_tools_lib as tools_lib import LiCSBAS_plot_lib as plot_lib class Usage(Exception): """Usage context manager""" def __init__(self, msg): self.msg = msg #%% Main def main(argv=None): #%% Check argv if argv == None: argv = sys.argv start = time.time() ver=1.1; date=20220209; author="Q. Ou" print("\n{} ver{} {} {}".format(os.path.basename(argv[0]), ver, date, author), flush=True) print("{} {}".format(os.path.basename(argv[0]), ' '.join(argv[1:])), flush=True) #%% Set default ifgfile = [] bperpfile = [] pngfile = 'network.png' bad_ifgfile = [] plot_bad_flag = True strong_connected = False suffix = "" strict = True months = False thresh = False remove_edge_cuts = True remove_node_cuts = True skip_node_cuts = False label_name = None #%% Read options try: try: opts, args = getopt.getopt(argv[1:], "hi:b:o:r:sm:t:enl:", ["help", "not_plot_bad", "not_strict", "skip_node_cuts"]) except getopt.error as msg: raise Usage(msg) for o, a in opts: if o == '-h' or o == '--help': print(__doc__) return 0 elif o == '-i': ifgfile = a elif o == '-b': bperpfile = a elif o == '-o': pngfile = a elif o == '-r': bad_ifgfile = a elif o == '--not_plot_bad': plot_bad_flag = False elif o == '-s': strong_connected = True elif o == '-m': months = a elif o == '-t': thresh = float(a) elif o == '--not_strict': strict = False elif o == '-e': remove_edge_cuts = False elif o == '-n': remove_node_cuts = False elif o == '--skip_node_cuts': skip_node_cuts = True elif o == '-l': label_name = a if not ifgfile: raise Usage('No ifg list given, -i is not optional!') elif not os.path.exists(ifgfile): raise Usage('No {} exists!'.format(ifgfile)) elif not bperpfile: raise Usage('No bperp list given, -b is not optional!') elif not os.path.exists(bperpfile): raise Usage('No {} exists!'.format(bperpfile)) except Usage as err: print("\nERROR:", file=sys.stderr, end='') print(" "+str(err.msg), file=sys.stderr) print("\nFor help, use -h or --help.\n", file=sys.stderr) return 2 #%% Read info ifgdates = io_lib.read_ifg_list(ifgfile) basename = os.path.basename(ifgfile).split('.')[0] if bad_ifgfile: bad_ifgdates = io_lib.read_ifg_list(bad_ifgfile) # ifgdates = list(set(ifgdates)-set(bad_ifgdates)) else: bad_ifgdates = [] if thresh: dt = tools_lib.calc_temporal_baseline(ifgdates) shortifg = [ifg for ifg, t in zip(ifgdates, dt) if t <= thresh] ifgdates = list(set(ifgdates)-set(shortifg)) suffix = suffix + "_dt_gt_{}".format(int(thresh)) # export list with open("{}_dt_le_{}.txt".format(basename, int(thresh)), 'w') as f: for i in shortifg: print('{}'.format(i), file=f) with open("{}_dt_gt_{}.txt".format(basename, int(thresh)), 'w') as f: for i in ifgdates: print('{}'.format(i), file=f) if months: ifgdates_allowed_months = tools_lib.select_ifgs_by_months(ifgdates, allowed_month=months, strict=strict) ifgdates_other_months = list(set(ifgdates) - set(ifgdates_allowed_months)) ifgdates = ifgdates_allowed_months suffix = suffix + "_months{}".format(months) # export list if strict: with open("{}{}_either_out.txt".format(basename, suffix), 'w') as f: for i in ifgdates_other_months: print('{}'.format(i), file=f) with open("{}{}_both_in.txt".format(basename, suffix), 'w') as f: for i in ifgdates: print('{}'.format(i), file=f) else: with open("{}{}_both_out.txt".format(basename, suffix), 'w') as f: for i in ifgdates_other_months: print('{}'.format(i), file=f) with open("{}{}_either_in.txt".format(basename, suffix), 'w') as f: for i in ifgdates: print('{}'.format(i), file=f) # extract bperp after modifying ifgdates imdates = tools_lib.ifgdates2imdates(ifgdates) bperp = io_lib.read_bperp_file(bperpfile, imdates) if strong_connected: strong_links, weak_links, edge_cuts, node_cuts = tools_lib.separate_strong_and_weak_links( ifgdates, "{}_stats.txt".format(basename), remove_edge_cuts=remove_edge_cuts, remove_node_cuts=remove_node_cuts, skip_node_cuts=skip_node_cuts) #pngfile = "{}{}_strongly_connected_network.png".format(basename, suffix) plot_lib.plot_strong_weak_cuts_network(ifgdates, bperp, weak_links, edge_cuts, node_cuts, pngfile, plot_weak=True) # export weak links with open("{}{}_weak_links.txt".format(basename, suffix), 'w') as f: for i in weak_links: print('{}'.format(i), file=f) # export strong links with open("{}{}_strong_links.txt".format(basename, suffix), 'w') as f: for i in strong_links: print('{}'.format(i), file=f) # export edge cuts print("{} ifgs are edge cuts".format(len(edge_cuts))) with open("{}_edge_cuts.txt".format(basename), 'w') as f: for i in edge_cuts: print('{}'.format(i), file=f) print('{}'.format(i)) # export edge cuts print("{} epochs are node cuts".format(len(node_cuts))) with open("{}_node_cuts.txt".format(basename), 'w') as f: for i in node_cuts: print('{}'.format(i), file=f) print('{}'.format(i)) else: #%% Plot image #pngfile = "{}{}_network.png".format(basename, suffix) plot_lib.plot_network(ifgdates, bperp, bad_ifgdates, pngfile, plot_bad_flag, label_name=label_name) #%% Finish elapsed_time = time.time()-start hour = int(elapsed_time/3600) minite = int(np.mod((elapsed_time/60),60)) sec = int(np.mod(elapsed_time,60)) print("\nElapsed time: {0:02}h {1:02}m {2:02}s".format(hour,minite,sec)) print('\n{} Successfully finished!!\n'.format(os.path.basename(argv[0]))) print('Output: {}\n'.format(pngfile), flush=True) #%% main if __name__ == "__main__": sys.exit(main())