Zf޺dZddlZddlZddlZddlZddlmZmZddl Z ddl m Z ddl mZmZmZmZddddd Zeej+Z d$d Zd%d efd Zd efdZd efdZd efdZd efdZ d&dZdZdefdZ dZ!dZ"dZ#d'dZ$GddZ%Gdde%Z&Gdde%Z'Gd d!e%Z(Gd"d#e%Z)y)(z1Class / utilities for GNSS download / operations.N)urlopen urlretrieve) coordinate)ptimereadfile time_funcutils1z7http://geodesy.unr.edu/NGLStationPages/DataHoldings.txtzShttp://garner.ucsd.edu/pub/measuresESESES_products/Velocities/ESESES_Velocities.txtz5https://sideshow.jpl.nasa.gov/post/tables/table2.html)UNRESESES JPL-SIDESHOWGENERICc ^|rtnd}|6|dk(r tdtjj t |}tjj |st||||dk(r t|}n2|dk(r t|}n!|dk(r t|}n|dk(r t|}tjdD cgc]} | jc} |d<|d |d d kDxxd zcc<|d t|dd ddj!|j#|d|dk\|d|dkz|d |dk\z|d |dkz} |d|dtj$| d|r[d|j#vrIt'j(|gdd} | |d| k\z} |d|dtj$| d|r[d|j#vrIt'j(|gdd} | |d| kz} |d|dtj$| d|?d|j#vr-| |d|k\z} |d|dtj$| d||d| |d| |d| |d | fScc} w)aSearch available GNSS sites within the geo bounding box for a given GNSS source. Parameters: SNWE - tuple of 4 float, indicating (South, North, West, East) in degrees source - str, program or institution that processed the GNSS data start_date - str, date in YYYYMMDD format end_date - str, date in YYYYMMDD format site_list_file - str, site list file name min_num_solution - int, minimum number of solutions available Returns: site_names - 1D np.ndarray in string, GNSS station names site_lats - 1D np.ndarray in float32, latitude site_lons - 1D np.ndarray in float32, longitude cyNargskwargss L/home/exouser/operations/rsmas_insar/tools/MintPy/src/mintpy/objects/gnss.pyzsearch_gnss..1r z9Site list file must be specified for GENERIC GNSS source!)source print_msgr r r sitelonhzload dz GNSS sites with fields:  latrzkeep sites within SNWE of z: []end_datezkeep sites with end_date >= start_datezkeep sites with start_date <= num_solutionz"keep sites with # of solutions >= )print ValueErrorospathbasenameGNSS_SITE_LIST_URLSisfiledload_site_listread_UNR_site_listread_ESESES_site_listread_JPL_SIDESHOW_site_listread_GENERIC_site_listnparrayupperlenjoinkeyssumrdate_list2vector) SNWEr'r&rsite_list_filemin_num_solutionrvprintsitesridxstart_dtend_dts r search_gnssrE#s U%AF Y XY YWW--.A&.IJN 77>>. )vK">2 8 %n5 > !+N; 9 &~6HHuV}EtdjjlEFE&M %Lu#$+$ U3uV}%a((A#((5::>( #  /qzM!;LDQYZ[P\] ^M8, Orr>c tj|tdd}|dddf|dddfjtj|dddfjtj|dddf|dddf|ddd fjtj d }tj |d Dcgc]"}tjj|d $c}|d <tj |d Dcgc]"}tjj|d $c}|d <|Scc}wcc}w)z;Return names and lon/lat values for UNR GNSS stations. r") rr"r#r$ )dtypeskiprowsusecolsNrr#rRrSrU)rr!rr'r&r(r'z%Y-%m-%dr&) r5loadtxtstrastypefloat32int16r6dtdatetimestrptime)r>fcrAxs rr1r1}s N#C[ \BAaCAaC 3AaC 3AaCAaCAbD2  E((QVWcQd#eABKK$8$8J$G#efE,((QVWaQb#cABKK$8$8J$G#cdE* L$f#cs >'D8'D=ctj|dt}|dddf|dddfjtj|dddfjtjd}|S)z;Return names and lon/lat values for JPL GNSS stations. )rWrVNrr"r#rrr!r5rYrZr[r\r>rarAs rr2r2s_ NRs ;BAaCAaC +AaC + E Lrctj|ddt}|ddddf|ddddfjtj|ddddfjtjd}|S) zDReturn names and lon/lat values for JPL-SIDESHOW GNSS stations. .rr) displacementvelocityzun-supported obs_type: ror$rO FILE_PATHgnss_rHhorzz.0f_z.csv)SiteLonLat DisplacementVelocityU10f8r"z$default GNSS observation file name: r,T)rV delimiternamesz"read GNSS observations from file: rlz calculating GNSS observation ...incidenceAngle azimuthAnglegeo)work_dircoordz)use incidence / azimuth angle from file: z+use incidence / azimuth angle from metadatar N)maxValuer/r )suffix) url_prefix)r'r& gnss_comp horz_az_anglemodelz!write GNSS observations to file: w)%r)r8lowerr+r,dirnamer9r7r/r5 genfromtxtsizer6floatzerosnan enumerateutget_geometry_filer-get_ESESES_url_prefixr progressBarupdateget_gnss_classget_los_velocityisnanappendrsite_lonsite_latcloseopencsvwriterwriterow writerows)$metaobs_type site_namesr'r&rrrrrredor@num_siteobs_indfile_dircsv_file col_names col_typesnum_rowrasite_obs temp_namestemp_obsi site_name data_list geom_filegeom_objrprog_bargnss_objveldis_tsdisfcwrbs$ r get_los_obsrs4 U%AF:H 3 3Y7Nxj5YY>>#~5a1GwwtK01Hww||Hi]&;? G ww~~h ]]894 P'' BGGNN8,H1D3HQ<@A ]]894 Pi() h IaL)Jxx6Hxx6?H )* 5 G 9 *"*:+B"CA"FHQK Gn Oc 12((*:>)J2:%I  H >rww?O?OPY?Z>[\ ]H @ A17(0B*, $$h)L%j1 ^LAy OOAaC1Q3%q !Ia=(IO J.~f-iJOH"33%!#+ 4KCHHSM"&&vbzF1I/EC   hmmX->->@Q@QSVX[\ ]% ^&  28A,?@ (C  %B**R.C LL # MM) $ % 88;AQwZ;< O % %  !  9 <qz1CDEErctdd}|dz }tjj}d}d}||krV|tj|z j d}|j |} t|t|d |SS#||d z k(rtd |d |d z }t|d YxYw) zlGet the url prefix for the ESESES source, which updates regularly. [Poor design of ESESES website]. z#searching for ESESES url_prefix ...z=http://garner.ucsd.edu/pub/measuresESESES_products/Timeseriesz3/CurrentUntarred/Clean_TrendNeuTimeSeries_comb_{:s}rdays%Y%m%dz [YES!]r"zThe ESESES repository z CANNOT be found!z [no]) r)r^rtoday timedeltastrftimeformatrFileNotFoundError)url_fmtrmax_daynum_dayday_strrs rrrgs /0MG DDG GGMMOEGG G 2<A"b"o N NUQqR[[aqa 0 NE N 1Q3u$$&A 1Q3u$$&A 1Q3u$$&A !Q> Os#Ec$tjtjj||dd}|Dcgc]3}tjj |j dd5}}t j|}|Scc}w)NrHz.*.posrIr")globr+r,r9r-splitryy2yyyy)gnss_dirrfnamesryearss r get_pos_yearsrss YYrww||Ha.?@ AF8> ?1RWW  a & &s +A . ?E ? MM% E L @s8B c t|dd}t|dd}||z dz}ggggf\}}} } t|D]_} t|| z} tjj ||dd| dddd} t | \}}}}||z }||z }| |z } | |z } atj|}tj|}tj| } tj| } tjj|d }tjj|d }tj|jtj }d |||k<d |||kD<||||| || |fS) NrrOr"rrIr#rHz.posrrlF)rrangerZr+r,r9rr5r6r^r_r`onesrbool_)rrr'r&year0year1num_yearrrrrryearirdatesiXiYiZidate0date1flags r read_GSI_F3rs|  1Q E 1 Eu}q HR^NE1a 8_EAI X$q59Q-t'DE*51B  R R R HHUOE  A  A  A KK X 6E KK 8 4E 77177"(( +DDD ;$4!D' 11rceZdZdZddefdZddZddZddZddZ d ed efd Z dd Z dd e de fdZ ddZ ddZ ddZy)GNSSaNGNSS parent class for time-series of daily solution. The GNSS class is solely meant to be a parent class. Child classes, defined below, support functions for downloading and parsing GNSS position based on the processing source (e.g., UNR, etc.). Use the `get_gnss_class` method to determine appropriate child class. The parent class/object will assign the following attributes: source - str, GNSS solution source version - str, GNSS solution version url_prefix - str, GNSS data file url prefix The chile class/object will assign the following attributes: file - str, path of the local data file url - str, path of the remote data file site - str, four-digit site code site_lat/lon - float, site latitude/longitude in degree dates - 1D np.ndarray, dt.datetime object date_list - list(str), dates in YYYYMMDD format dis_e/n/u - 1D np.ndarray, displacement in meters std_e,n,u - 1D np.ndarray, displacement STD in meters Nrc||_||_||_||_d|_d|_|r|n/t jjd|j|_ t jj|js5td|jt j|jd|_d|_d|_d|_d|_d|_d|_d|_y)NzGNSS-zcreate directory:)rrversionrurlfiler+r,abspathr7data_direxistsr)mkdirr date_listdis_edis_ndis_ustd_estd_nstd_u)selfrr rrrs r__init__z GNSS.__init__s   $ $,"''//E&,,.IYBZ2[ ww~~dmm, %t}} 5 HHT]] #       rctjj|js|j ||j |j |y)z^Read the lat/lon and displacement data of the station. Download if necessary. rN)r+r,r/r dload_siteget_site_lat_lonread_displacement)r*rrs rrz GNSS.opensE ww~~dii( OOiO 0  3rc |rtnd}|jr|s)tjj |j s{|d|j dd|jd|j d|}|dkDrC t|j|j |d|j |j S|j S#|d|jd ||z |d z}YxxYw) :Download GNSS site data file. Parameters: overwrite - bool, overwrite existing data file total_tries - int, number of tries to download if failed print_msg - bool, verbose print out msg Returns: self.file - str, path to the local data file cyrrrs rrz!GNSS.dload_site.. rrzdownloading site rHz from rJrzsuccessfully downloaded: zerror downloading z on trial no. r") r)rr+r,r/rrrr)r* overwrite total_triesrr@ remain_triess rr.zGNSS.dload_sites$)E 88  )B &tyym6$++d499UV-X Y'L"$))46txxjAB yytyy /z T`H`Gabc A%Ls <5C "C.ctd)z~Get the GNSS site latitude & longitude into: Returns: site_lat/lon - float, site latitude/longitude in degree zAget_site_lat_lon() is NOT implemented. Override with child class.NotImplementedError)r*rs rr/zGNSS.get_site_lat_lons""effrctd)a]Get the GNSS time/displacement(Std) into: Returns: dates - 1D np.ndarray in datetime.datetime object date_list - list(str), date in YYYYMMDD format dis_e/n/u - 1D np.ndarray in float32, displacement in meters std_e/n/u - 1D np.ndarray in float32, displacement STD in meters zBread_displacement() is NOT implemented. Override with child class.r8)r*r'r&rdisplays rr0zGNSS.read_displacement%s""fggrr'r&ctjt|jt}|r.t j |gdd}d||j|k<|r.t j |gdd}d||j|kD<|j||_|j||_|j||_|j||_ |j||_ |j||_ |j||_ |jDcgc]}|jdc}|_ycc}w)zsCrop the time-series given the start/end_date in format YYYYMMDD, and create date_list from dates. rlrrN)r5r r8rboolrr<r$r%r&r'r(r)rr#)r*r'r&rt0t1rbs r_crop_to_date_rangezGNSS._crop_to_date_range/swws4::d3 '' 5a8;B$%Db ! '' 3A6q9B$%Db !ZZ% ZZ% ZZ% ZZ% ZZ% ZZ% ZZ% 9= C1!**X.CCs)E cddlm}|j|j|j ddd\}}t |dz|}|dj |j|jfi||dj |j|jfi||dj |j|jfi||rt d | }|dj|j|jfd |ji||dj|j|jfd |ji||dj|j|jfd |ji|tgd D]\}} ||j| d |j!|j"dd|j$dd|j'|j)||fS)z+Plot the displacement time-series. rNr$r"T)nrowsncolssharexr#)rHcnone) linestylecoloryerr)EastNorthUpz [m]rHz ())matplotlib.pyplotpyplotr$rsubplotsdictscatterrr%r&errorbarr'r(r)r set_ylabelsuptitlerr tight_layoutshow) r* marker_size marker_colorplot_error_barpltfigaxrrlabels rplotz GNSS.plotHs ( ::  IIK,,Qa,=R Q,71 djj$**771 djj$**771 djj$**77 F,?F BqENN4::tzz M Mf M BqENN4::tzz M Mf M BqENN4::tzz M Mf M""9: -HAu qE  wd^ , -  !}Bt{{1oQ78   Bwr inc_angleaz_anglec|j|jtj|||j |}|j|dz|j |dzz|j |dzz|_|j|dzdz|j|dzdzz|j|dzdzzdz|_ |j|jfS)aConvert displacement in ENU to LOS direction. Parameters: inc_angle - float, LOS incidence angle in degree az_angle - float, LOS aziuth angle in degree from the north, defined as positive in clock-wise direction gnss_comp - str, GNSS components used to convert to LOS direction horz_az_angle - float, fault azimuth angle used to convert horizontal to fault-parallel measured from the north with anti-clockwise as positive Returns: dis_los - 1D np.ndarray for displacement in LOS direction std_los - 1D np.ndarray for displacement standard deviation in LOS direction ) los_inc_angle los_az_anglecomprrr"r#r) r$rr%get_unit_vector4component_of_interestrr%r&dis_losr'r(r)std_los)r*r`rarrunit_vecs rdisplacement_enu2loszGNSS.displacement_enu2losks ::  IIK;;#!"'  **x{2**x{23**x{23 !JJ!4q8 JJ!4q89 JJ!4q89>AB ||T\\))rc|j\}}t|trtj|}t ||}|j |||dd\}}td|}tt|ddz |}td|}tt|ddz |}t|||dz|dzf|} tj|fd d i| dd } tj|fd d i| dd } | | fSt|tr=tj|d| } tjt|d} | | fSt!d|)zNGet the Line-of-Sight geometry info in incidence and azimuth angle in degrees.) lookup_filer-rr#LENGTHr"WIDTH)boxr datasetNamer)rrr) dimensionrHEADINGz(input geom_obj is neither str nor dict: )r/ isinstancerZrread_attributer geo2radarmaxminrrQreadrincidence_angleheading2azimuth_anglerr*) r*rrr!ratrrrrbrr`ras rget_los_geometryzGNSS.get_los_geometryss((*S h $))(3Cs9E??3y?A!AFDAqAq ACH $6q$8! d+good.CCA,,t,a/'//$2G2JJSPA  Q !Mc3 =88rc |j||||||dd\} } d} t| dkrd} nR|rP|rNtj|gdd} tj|gdd} | d| dz | | z dz krd} | r| Dcgc]"}tj j |d $}}tj|| }tjtjj|| d |_ n0tj|_ |rtd |j d |j| fScc}w)aConvert the three-component displacement data into LOS velocity. Parameters: geom_obj - dict / str, metadata of InSAR file, or geometry file path start_date - str, YYYYMMDD format end_date - str, YYYYMMDD format ref_site - str, reference GNSS site gnss_comp - str, GNSS components used to convert to LOS direction horz_az_angle - float, fault azimuth angle used to convert horizontal to fault-parallel model - dict, time function model, e.g. {'polynomial': 1, 'periodic': [1.0, 0.5]} Returns: dates - 1D np.ndarray, datetime.datetime object dis - 1D np.ndarray, displacement in meters )r'r&r~rrNr#TFrrrOr)rr"z& Velocity calculation failed for site rH)rr8rr<r^r_rrget_design_matrix4time_funcr5dotlinalgpinvrprr)r)r*rr'r&r~rrrrrrdis2velr>r?rr#As rrzGNSS.get_los_velocitysG$.. !' /  1 s u:?G H'' 5a8;B'' 3A6q9BRy58#R1}4 DIJq--a:JIJ55iuMAFF299>>!#4c:1=DMFFDM? !}MN}}c!!Ks 'E)NIGS14r N)NTFrPTFNNTF)r#kT)enu2losV)NNNrrF)NNNrrNT)__name__ __module__ __qualname____doc__rZr+rr.r/r0r@r_rrjr|rrrrrrrs.S8 46g hDcDSD2 FOX+/"*U"*U"*J#8W[PU29jSWX\1"rrcPeZdZdZddeffd Zd fd Zd deeffdZd dZ xZ S) rzGNSS child class for daily solutions processed by Nevada Geodetic Lab at University of Nevada, Reno (UNR). Website: http://geodesy.unr.edu/NGLStationPages/GlobalStationList rcLt||||d||dk(rBtjj |j |j dd|dd|_nR|dk(r>tjj |j |j dd|_ntd|d |jsd |j|_ tjj |jtjj|j|_ y) Nr rr rrrIGS08rHrIz.tenv3rzUn-supported GNSS versoin: !z,http://geodesy.unr.edu/gps_timeseries/tenv3/) superr+r+r,r9r rrr*rrr-rr*rr rr __class__s rr+zGNSS_UNR.__init__s !   g  T]]tyym1WQKv4VWDI   T]]tyym64JKDI:7)1EF F  LT\\N[DO77<<1A1A$))1LMrc$t|||tjj |j d|j d}tjjtjj|s4tjtjj|dddd|j}tjj ||j d}t|||jS) r2r4rzpic/z.pngT)exist_okz/http://geodesy.unr.edu/tsplots/IGS08/TimeSeriesz5http://geodesy.unr.edu/tsplots/IGS14/IGS14/TimeSeries)rr) rr.r+r,r9r rr!rmakedirsrrr)r*r4r5r plot_filer plot_file_urlrs rr.zGNSS_UNR.dload_site4s Y)D GGLL$tyyk0FG ww~~bggooi89 KK 2T BHM  ,,  ZDII;d1CD  M9-yyrrFcHtjj|js|j |t j |jtdd}|dddfjt\|_ |_ |j|jfS)zGet station lat/lon from the displacement file. Modifies: self.lat/lon - float Returns: self.lat/lon - float r-r"rU)rVrWmax_rowsrr) r+r,r/rr.r5rYbytesr[rrr)r*rdatas rr/zGNSS_UNR.get_site_lat_lonUsvww~~dii( OOiO 0zz$))51rJ'+Ar"uH~'<'.lrrr->reading time and displacement in east/north/vertical directionr")rVrWzXError occurred while reading, probably due to interuptions during previous downloading. z Remove the file and re-download.TrNz%y%b%d)rSrU )r)r+r,r/rr.r5rYrr[rZr6r^r_r`rr\Tr$r%r&r'r(r)r@r_) r*r'r&rr;r@ramsgrbs rr0zGNSS_UNR.read_displacementds$)Eww~~dii( OOiO 0 OP LDIIUQ?FFsKBXX"QPQT(SQr{{33Ax@ST !//077 CEE          X6  IIK  DJJ  DJJ 4 4- LlC 5 5C #J OOdiO @DIIUQ?FFsKBSs9F#'HA HNrNrrr rrrrrZr+r.rr/r0 __classcell__rs@rrrs1 NSN2B ,E5> ,*4rrcPeZdZdZddeffd Zd fd Zd deeffdZd dZ xZ S) ra%GNSS child class for daily solutions processed for the Enhanced Solid Earth Science ESDR System (ESESES) project by JPL and SOPAC. Website: https://cddis.nasa.gov/Data_and_Derived_Products/GNSS/ESESES_products.html http://garner.ucsd.edu/pub/measuresESESES_products/ rct||||d|tjj |j |j jdd|_y)Nr rrHzCleanTrend.neu.Z) rr+r+r,r9r rrrrs rr+zGNSS_ESESES.__init__sV !  GGLL499??3DQ2GGW0XY rc(ddlm}|jst|_tj j |jtj j|j|_ t|-||||jd5}|j|jddd|jjd|_|rtd|jd|jS#1swYPxYw) z!Download GNSS data file. r)ZipFilerrN.Zz... extracted to rH)zipfilerrrr+r,r9r-rrrr. extractallr stripr))r*r4r5rrfzrs rr.zGNSS_ESESES.dload_sites $35DO77<<1A1A$))1LM Y)DTYY $ ) MM$-- ( )IIOOD)  %dii]3 4yy  ) )s DDrFc.tjj|js|j ||jj dr|jdd|_t |j5}|j}|Dcgc]}|jds|c}djd}t|jd|_ |Dcgc]}|jd s|c}djd}t|jd|_ |xj|jd krdnd zc_ ddd|j|jfScc}wcc}w#1swY+xYw) Get station lat/lon based on processing source. Retrieve data from the displacement file. Modifies: self.lat/lon - float Returns: self.lat/lon - float r-rNz # Latituder rz# East Longituderr)r+r,r/rr.endswithr readlines startswithrrrrr)r*rflinesrblat_linelon_lines rr/zGNSS_ESESES.get_site_lat_lonsEww~~dii( OOiO 0 99  d # #2DI $))_ @KKME$)GaALL,FGJPPQUVH!(.."22"67DM$)MaALL9K,LMaPVVW[\H!(.."22"67DM MM$--3"6QC ?M @}}dmm++HN @ @s7F F,F0:F *FFAF  F  Fc |rtnd}tjj|js|j ||jj dr|jdd|_|dtj|jttdd }|jd}t|Dcgc]L}tjt||d fd d tjt||d f zN} }tj | |_|ddd dfj%tj&j(dz \|_|_|_|_|_|_|j7|||r|j9|j"|j,|j*|j.|j2|j0|j4fScc}w)Read GNSS displacement time-series (defined by start/end_date). Parameters: start/end_date - str, date in YYYYMMDD format Returns: dates - 1D np.ndarray of datetime.datetime object dis_e/n/u - 1D np.ndarray of displacement in meters in float32 std_e/n/u - 1D np.ndarray of displacement STD in meters in float32 cyrrrs rrz/GNSS_ESESES.read_displacement..rrr-rNrrrrT)rXr"r#rr$i)r)r+r,r/rr.rr5rYtupler rr^r_rrr6rr[r\rr%r$r&r(r'r)r@r_) r*r'r&rr;r@rar(rrs rr0zGNSS_ESESES.read_displacements$)Eww~~dii( OOiO 0 99  d # #2DI OP ZZ 5q+< =xx{  -/SAqD]Aq1BLLc"QPQT(m4TT//XXe_ !QqS&z((466=          X6  IIK  DJJ  DJJ 4 4'/sAHrrrrrrs@rrrs0 ZS Z".,E5>,>.4rrcDeZdZdZddeffd ZddeeffdZd dZxZ S) rzoGNSS class for daily solutions processed by JPL-SIDESHOW. Website: https://sideshow.jpl.nasa.gov/pub/ rczt||||d|tjj |j |j dd|_|jsd|_tjj |jtjj|j|_ y)Nr rrHz.serieszJhttps://sideshow.jpl.nasa.gov/pub/JPL_GPS_Timeseries/repro2018a/post/point) rr+r+r,r9r rrrr-rrs rr+zGNSS_JPL_SIDESHOW.__init__s !!  GGLL499Q-w0GH jDO77<<1A1A$))1LMrrFctjjtd}t |5}|D]4}|dd|j k(s|dddk(s|j dd\}}6 dddt|_t|_ |dk(r'td |jd d |jd |j|jfS#1swYmxYw) rr NrOrPrSPOSr#T rz, ) r+r,r-r.rrrrrrr))r*rr> site_listlinerrs rr/z"GNSS_JPL_SIDESHOW.get_site_lat_lon's))*=n*MN. ! ;Y! ;!H )Qq U0B)-a):&Hh ; ; h h   Bt}}Q'r$--):; <}}dmm++ ; ;sC CCC"c Ntjj|js|j ||dk(r t dt j|j}|jd}t jt|Dcgc].}tj||ddfjt0c}|_|ddddfjt j j"\|_|_|_|_|_|_|j1|||dk(r|j3|j|j$|j&|j(|j*|j,|j.fScc}w) aRead GNSS displacement time-series (defined by start/end_date) Parameters: start/end_date - str, date in YYYYMMDD format Returns: dates - 1D np.ndarray of datetime.datetime object dis_e/n/u - 1D np.ndarray of displacement in meters in float32 std_e/n/u - 1D np.ndarray of displacement STD in meters in float32 r-TrriNr"rR)r+r,r/rr.r)r5rYrr6r r^r_r[rrr\rr$r%r&r'r(r)r@r_)r*r'r&rr;rn_datars rr0z#GNSS_JPL_SIDESHOW.read_displacementAsCww~~dii( OOiO 0   R Szz$))$AXXeTZm\r{{D23K,>,>s,CD\] AqsF|**2::688          X6 d? IIK  DJJ  DJJ 4 4#]s3F"rrr) rrrrrZr+rr/r0rrs@rrrs+NSN$,E5>,4'4rrcDeZdZdZddeffd ZddeeffdZd dZxZS) raGNSS class for daily solutions of an otherwise-unsupported source. Required local files: 1. GenericList.txt: site list file, with the following 3, 9 or 11 space-separated columns: SITE lat lon [vel_e vel_n vel_u err_e err_n err_u] [start_date end_date] where site is the four-digit, alphanumeric (uppercase) site code; and lat/lon are in decimal degrees. If included, vel should be in units of m/yr; and dates should be in format YYYYMMDD. 2. .txt: data file for each site, with 7 space-separated columns: date dis_e dis_n dis_u std_e std_n std_u where date is in the format or T or :, displacement values are in meters. rct||||d|tjj |j |j dd|_d|_d|_ y)Nr rrHz.txtrt) rr+r+r,r9r rrrrrs rr+zGNSS_GENERIC.__init__s] !  GGLL499Q-t0DE rrFctd}|djj|j}|d||d|fS)z8Get station lat/lon based on processing source. zGenericList.txtrr!r)r4rindexr)r*rrAinds rr/zGNSS_GENERIC.get_site_lat_lonsM''89Fm""$**4995U|C %,s"333rctjj|js|j ||dk(r t dt j|j}t jtj|dddf|_ |ddttddfjt jj \|_|_|_|_|_|_|j/|||r|j1|j|j"|j$|j&|j(|j*|j,fS)rr-TrNrr"rR)r+r,r/rr.r)r5rYr6rr<rrr r[r\rr$r%r&r'r(r)r@r_)r*r'r&rr;ras rr0zGNSS_GENERIC.read_displacementsww~~dii( OOiO 0   R S ZZ "XXe44R1X>? !U51:../66rzzBDD          X6  IIK  DJJ  DJJ 4 4rrrr) rrrrrZr+r/r0rrs@rrrks)&S"4C:4%4rr)NNr N2T)Nr T)r rrNTF)NN)*rrr_r^rr+urllib.requestrrnumpyr5mintpy.objects.coordr mintpy.utilsrrrr rr.rr:rKrErZr0r1r2r3r4rrrrrrrrrrrrrrrrs!7  /+AAOjL  ',,./ TX/3ADHC,c(   s #6[dEJvr6 F# F L"2