o tj|tdddddt}|dddf}|ddddftjj\}}|t|d d 8}td d |ddd ftD} td d |dddftD} ||dk||dk||dk||dk} |rt |gdd} | | | k9} |rt |gdd} | | | k9} || || || fS)aHUsed for search available GPS sites within the geo bounding box from CGPS database Parameters: SNWE : tuple of 4 float, indicating (South, North, West, East) in degrees start_date : string in YYYYMMDD format end_date : string in YYYYMMDD format site_list_file : string. min_num_solution : int, minimum number of solutions available Returns: site_names : 1D np.array of string for GPS station names site_lats : 1D np.array for lat site_lons : 1D np.array for lon ,)rrGBKdtype delimiterskiprowsusecolsencodingNrrgv@cSg|]}t|dqSz%Y%m%ddtstrptime.0irH/home/exouser/operations/rsmas_insar/tools/MimtPy/mimtpy/objects/cgps.py &zsearch_gps..r cSrrrrrrrr'rr r) nploadtxtstrastypefloat32Troundarrayptimedate_list2vector)SNWE start_dateend_datesite_list_fileZ pring_msgtxt_data site_names site_lons site_latst_startt_endidxt0t1rrr search_gpss$ "&&    r6c@s\eZdZdZdddZdddZddd Zdd d ZddedefddZ dde fddZ dS)CGPSaxGPS class for CGPS velocity Example: import matplotlib.pyplot as plt from mimtpy.objects.cgps import CGPS from mintpy.utils import utils as ut gps_obj = GPS(site='GV05', data_dir='~/insarlab/GPS') gps_obj.open() vel_los = ut.enu2los(gps_obj.vel_e, gps_obj.vel_n, gps_obj.vel_u) /data/lxr/insarlab/GPS/CGPS/cCs||_||_d|_d|_dS)Nz6GNSS data products of China earthquake adiministrationz'/data/lxr/insarlab/GPS/CGPS/GPS_vel.txt)sitedata_dirsourcefile)selfr9r:rrr__init__Es z CGPS.__init__TcCs&tj|js td|j|ddS)N3GPS velocity dataset cannot be found. Please check! print_msg)ospathisfiler< Exception read_velocity)r=rArrropenLsz CGPS.opencCsv|rtdtj|jstdtj|jtddddd t}|ddd f}||j k}t |ddd f||_ t |ddd f||_ t |ddd f|d |_t |dddf|d |_t |dddf|d |_t |dddf|d |_|dddf| tjd d |_|dddf| tjd d |_|j |j |j|j|j|j|j|jfS)zread velocity datasetzRead velocity datasetr?rr) rrrrr r r  r r Nrrrr drHr r rIrJ)printrBrCrDr<rErr r!r"r9floatsite_latsite_lonvel_estd_evel_nstd_nfloat64vel_ustd_u)r=rAdatar.Z id_positionrrrrFQs& &&  zCGPS.read_velocityFc Cs|j}|j}t|trQt|}t||d}|j|||ddd\}}|||d|df} tj|d| |ddd} tj|d | |ddd} t | } | | fSt|t rst j |d|d } t |d } | d d krot | } | | fStd|)N) lookup_filer@rrrincidenceAngle) datasetNameboxrA)rr azimuthAngle) dimensionrAHEADINGf@gF@z)input geom_obj is neight str nor dict: {})rNrO isinstancer!rread_attributer geo2radarreadutazimuth2heading_angledictincidence_anglerM ValueErrorformat) r=geom_objrAlatlonatrcoordyxr[ inc_angleaz_angle head_anglerrrget_los_geometryjs$        zCGPS.get_los_geometryenu2losrqrscCs|tjd9}|tjd9}t|t|dt|t|t|g}|}|dvr1n|dvr:d|d<n|dvrGd|d<d|d <ntd t||j|d|j|d |j |d|_ |j |dd|j |d d|j |ddd |_|j |jfS) a)Convert displacement in ENU to LOS direction Parameters: inc_angle : float, local incidence angle in degree head_angle : float, satellite orbit heading direction in degree from the north, defined as positive in clock-wise direction gps_comp : string, GPS components used to convert to LOS direction Returns: dis_los : 1D np.array for displacement in LOS direction std_los : 1D np.array for displacement standard deviation in LOS direction r_ru)en2loshz2losgr)u2losup2losrrzUn-known input gps components:g?)rpisincoslowerrhr!rPrRrUvel_losrQrSrVstd_los)r=rqrsgps_compunit_vecrrrdisplacement_enu2loss8       zCGPS.displacement_enu2losrcCs"||\}}|j|||ddS)aRead GPS velocity in LOS direction Parameters: geom_obj : dict / str, metadata of InSAR file, or geometry file path gps_comp : string, GPS components used to convert to LOS direction Retruns: vel : 1D np.array of displacement in meters std : 1D np.array of displacement uncertainty in meters )rN)rtr)r=rjrrArqrsrrrread_gps_los_velocityszCGPS.read_gps_los_velocityN)r8)T)Frw)ruF) __name__ __module__ __qualname____doc__r>rGrFrtrMrr!rrrrrr78s   &r7) rBcodecsrrnumpyrmintpy.objects.coordr mintpy.utilsrZ CGPS_datasetr6r7rrrrs    #