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 """ ======== Overview ======== This script checks the residual of LiCSBAS13 step for sbovl data to drop the highly noisy interferograms =============== Input & output files =============== Inputs in TS_GEOCml*/ : - info/ - - 13resid.txt Outputs in TS_GEOCml*/ : - info/ - 11bad_ifg.txt : List of updated bad ifgs 11 to rerun step13 discarded from further processing ===== Usage ===== LiCSBAS11_check_unw.py [-t tsadir] """ #%% Change log ''' M. Nergizci, 2025-08-07 ''' #%% Import from LiCSBAS_meta import * import getopt import os import sys import time import shutil import numpy as np import datetime as dt import LiCSBAS_io_lib as io_lib import LiCSBAS_tools_lib as tools_lib import LiCSBAS_plot_lib as plot_lib from scipy import stats class Usage(Exception): """Usage context manager""" def __init__(self, msg): self.msg = msg #%% Helper def days_between(pair): try: d1_str, d2_str = pair.split('_') d1 = datetime.strptime(d1_str, '%Y%m%d') d2 = datetime.strptime(d2_str, '%Y%m%d') return abs((d2 - d1).days) except: return -1 #%% Main def main(argv=None): #%% Check argv if argv == None: argv = sys.argv start = time.time() 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 tsadir = [] #%% Read options # Parse options try: opts, args = getopt.getopt(argv[1:], "ht:", ["help"]) except getopt.GetoptError as err: print(str(err)) print(__doc__) return 1 for o, a in opts: if o in ("-h", "--help"): print(__doc__) return 0 elif o == "-t": tsadir = a if not tsadir or not os.path.isdir(tsadir): print("ERROR: Please provide a valid TS directory using -t", flush=True) return 1 infodir = os.path.join(tsadir, 'info') restxtfile = os.path.join(infodir,'13resid.txt') badifg_file = os.path.join(infodir, '11bad_ifg.txt') if not os.path.isfile(restxtfile): print("ERROR: Residual file {} does not exist!".format(restxtfile), flush=True) return 1 pairs, rms_values = [], [] with open(restxtfile, 'r') as f: for line in f: if line.startswith('#') or not line.strip(): continue parts = line.strip().split() if len(parts) >= 2: pairs.append(parts[0]) rms_values.append(float(parts[1])) pairs = np.array(pairs) rms_values = np.array(rms_values) # Zero residuals zero_mask = rms_values == 0.0 zero_pairs = pairs[zero_mask] # Remove zeros for stats valid_mask = rms_values != 0.0 pairs = pairs[valid_mask] rms_values = rms_values[valid_mask] # Thresholding mean = np.mean(rms_values) std = np.std(rms_values) threshold = mean + 2 * std outlier_indices = np.where(rms_values > threshold)[0] outlier_pairs = pairs[outlier_indices] # Combine all unique problematic pairs bad_set = set() bad_set.update(zero_pairs) bad_set.update(outlier_pairs) # # Sort and write output # bad_list = sorted(bad_set) # print(f" Total bad interferograms: {len(bad_list)}") # if bad_list: # with open(badifg_file, 'w') as f: # for pair in bad_list: # f.write(pair + '\n') # print(f" Written to: {badifg_file}") # else: # print(" No bad interferograms found.") # Combine with existing if present existing_pairs = set() if os.path.exists(badifg_file): with open(badifg_file, 'r') as f: for line in f: if line.strip(): existing_pairs.add(line.strip()) # Merge and sort all unique bad interferograms all_bad_pairs = existing_pairs.union(bad_set) bad_list = sorted(all_bad_pairs) # Write updated file with open(badifg_file, 'w') as f: for pair in bad_list: f.write(pair + '\n') return 0 #%% main if __name__ == "__main__": sys.exit(main())