""" 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'], \