def read_gps(self): for sname in self.site_names: site = {} site['name'] = sname gps_obj = GPS(sname, data_dir=self.gps_dir) gps_obj.open(print_msg=False) site['lat'] = gps_obj.site_lat site['lon'] = gps_obj.site_lon (site['gps_datetime'], site['gps_dis'], site['gps_std']) = gps_obj.read_gps_los_displacement(self.geom_file, self.start_date, self.end_date, ref_site=self.ref_site, gps_comp='enu2los')[0:3] site['reference_site'] = self.ref_site self.ds[sname] = site sys.stdout.write('\rreading GPS {}'.format(sname)) sys.stdout.flush() print() return
def ngps_process(gpsname, geometry, gps_dir): """process ngps sites""" gps_obj = GPS(site=gpsname, data_dir=gps_dir) print('process {} site'.format(gpsname)) gps_obj.open() if not geometry: dis_los = ut.enu2los(gps_obj.dis_e, gps_obj.dis_n, gps_obj.dis_u) dates = gps_obj.dates los_vel = dis_velocity(dates, dis_los) else: dates, dis_los = gps_obj.read_gps_los_displacement(geometry) gps_obj.get_gps_los_velocity(geometry) los_vel = gps_obj.velocity e_vel = dis_velocity(dates, gps_obj.dis_e) n_vel = dis_velocity(dates, gps_obj.dis_n) site_inve = np.array([[ gps_obj.site, gps_obj.site_lat, gps_obj.site_lon, los_vel, n_vel, e_vel ]]) return site_inve