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 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