def config_sync():
"""
Synchronise station metadata from inventory/dataless.
"""
import glob
from ..api import connect, get_config, get_stations, update_station,\
preload_instrument_responses
db = connect()
responses = preload_instrument_responses(db)
netsta = []
for id, row in responses.iterrows():
net, sta, loc, chan = row["channel_id"].split(".")
netsta.append("%s.%s"%(net,sta))
responses["netsta"] = netsta
for station in get_stations(db):
id = "%s.%s" % (station.net, station.sta)
coords = responses[responses["netsta"] == id]
lon = float(coords["longitude"].values[0])
lat = float(coords["latitude"].values[0])
update_station(db, station.net, station.sta, lon, lat, 0, "DEG", )
logging.info("Added coordinates (%.5f %.5f) for station %s.%s" %
(lon, lat, station.net, station.sta))
db.close()
def info_stations(db):
"""
Show information about configured stations.
"""
from ..api import get_stations
click.echo('')
click.echo('Stations:')
click.echo(' NET.STA Long. Lat. Alt. Coord. Used?')
s = None
na_sign = '-'
for s in get_stations(db, all=True):
click.echo(' {:>8s} {:>9s} {:>8s} {:>6s} {:3s} {:1s}'.format(
'{}.{}'.format(s.net, s.sta),
'{:.4f}'.format(s.X) if s.X is not None else na_sign,
'{:.4f}'.format(s.Y) if s.Y is not None else na_sign,
'{:.1f}'.format(s.altitude) if s.altitude is not None else na_sign,
s.coordinates or na_sign,
'Y' if s.used else 'N'))
if s is None:
click.echo(' ')
def da_stations_update_loc_chan():
"""EXPERT MODE: Populates the Location & Channel from the Data Availability
table. Warning: rewrites automatically, no confirmation."""
from msnoise.api import connect, get_stations
session = connect()
stations = get_stations(session)
for sta in stations:
data = session.query(DataAvailability). \
filter(text("net=:net")). \
filter(text("sta=:sta")). \
group_by(DataAvailability.net, DataAvailability.sta,
DataAvailability.loc, DataAvailability.chan). \
params(net=sta.net, sta=sta.sta).all()
locids = sorted([d.loc for d in data])
chans = sorted([d.chan for d in data])
print("%s.%s has locids:%s and chans:%s" % (sta.net, sta.sta,
locids, chans))
sta.used_location_codes = ",".join(locids)
sta.used_channel_names = ",".join(chans)
session.commit()