"""
import sys
import re
sys.path.append('../pyapi/')
from georeader import GeoReader
from sqlitewriter import SqLiteWriter
if len(sys.argv) < 3:
print("Usage {} file.geo database.sqlite".format(sys.argv[0]))
exit(0)
fn = sys.argv[1][:-4]
coo_rdr = GeoReader(fname=fn+'.coo')
geo_rdr = GeoReader(fname=fn+'.geo')
coo_wrt = SqLiteWriter(db=sys.argv[2], table='monitoring_coo',
filt=['id', 'east', 'north', 'elev', 'datetime'])
geo_wrt = SqLiteWriter(db=sys.argv[2], table='monitoring_obs',
filt=['id', 'hz', 'v', 'distance',
'crossincline', 'lengthincline', 'datetime'])
buf = coo_rdr.Load()
for b in buf:
coo_wrt.WriteData(b)
buf = geo_rdr.Load()
for b in buf:
geo_wrt.WriteData(b)
print(
"Usage: freestation.py input_file gama_path station_id station_height"
)
sys.exit(-1)
if fname[-4:] not in ['.geo', '.coo', '.dmp', '.csv']:
fname += '.geo'
if len(sys.argv) > 2:
gama_path = sys.argv[2]
else:
gama_path = '/home/siki/GeoEasy/gama-local'
# load observations and coordinates
fn = fname[:-4] # remove extension
ext = fname[-4:]
if ext in ['.geo', '.coo']:
obs = GeoReader(fname=fn + '.geo')
else:
obs = CsvReader(fname=fn + '.dmp')
# load observations
observations = obs.Load()
# load coordinates and add to adjustment
if ext in ['.geo', '.coo']:
coo = GeoReader(fname=fn + '.coo')
else:
coo = CsvReader(fname=fn + '.csv')
n = 0 # number of points
st = False # station found
coords = coo.Load()
f = Freestation(observations, coords, gama_path)
print(f.Adjustment())
from csvreader import CsvReader
logging.getLogger().setLevel(logging.WARNING)
if len(sys.argv) > 1:
ifname = sys.argv[1]
else:
#ifname = 'test.geo'
print("Usage: blindorientation.py input_file totalstation port")
sys.exit(-1)
if ifname[-4:] != '.dmp' and ifname[-4:] != '.geo':
ifname += '.geo'
if ifname[-4:] == '.geo':
g = GeoReader(fname=ifname)
else:
g = CsvReader(fname=ifname)
data = g.Load()
stationtype = '1100'
if len(sys.argv) > 2:
stationtype = sys.argv[2]
port = '/dev/ttyUSB0'
if len(sys.argv) > 3:
port = sys.argv[3]
if re.search('120[0-9]$', stationtype):
from leicatps1200 import LeicaTPS1200
mu = LeicaTPS1200()
elif re.search('110[0-9]$', stationtype):
from leicatcra1100 import LeicaTCRA1100
mu = LeicaTCRA1100()
elif re.search('180[0-9]$', stationtype):
from leicatca1800 import LeicaTCA1800
mu = LeicaTCA1800()
elif not config:
met_addr = None
if len(sys.argv) > 9:
met_par = sys.argv[9]
elif not config:
met_par = None
# load input data set
coo_filt = ['id', 'east', 'north', 'elev']
if ifname[-4:] in ('.geo', '.coo'):
g = GeoReader(fname=ifname[:-4] + '.geo')
f = GeoReader(fname=ifname[:-4] + '.coo', filt=coo_filt)
else:
g = CsvReader(fname=ifname[:-4] + '.dmp')
f = CsvReader(fname=ifname[:-4] + '.csv', filt=coo_filt)
directions = g.Load()
coordinates = f.Load()
# writers
if ofname[-4:] == '.dmp' or ofname[-4:] == '.csv' or ofname == 'stdout':
# dmp/csv file or console output
if ofname[-4:] == '.dmp' or ofname[-4:] == '.csv':
ofname1 = ofname[:-4] + '.dmp'
ofname2 = ofname[:-4] + '.csv'
else:
ofname1 = ofname2 = ofname
dmp_wrt = CsvWriter(angle = 'DMS', dist = '.4f', \
filt = ['station', 'id','hz','v','distance', 'datetime'], \
fname = ofname1, mode = 'a', sep = ';')
coo_wrt = CsvWriter(dist = '.4f', \
filt = ['id', 'east', 'north', 'elev', 'datetime'], \
