def parse_central_point(cpoint): """ parses central point latitude, longitude and height in degrees and meters numbers are separated by commas cpoint: string: latitude,longitude,height (height is optional) return: CoordSystemLocal instance """ cpoint = cpoint.split(",") if len(cpoint) < 2: raise ParseCentralPointError, "latitude and longitude not set: %s"%\ cpoint if len(cpoint) > 3: raise ParseCentralPointError, "too namy numbers set: %s" % cpoint from gizela.util.Converter import Converter lat = Converter.deg2rad_(float(cpoint[0])) lon = Converter.deg2rad_(float(cpoint[1])) if len(cpoint) is 3: height = float(cpoint[2]) else: height = 0.0 #csl = CoordSystemLocal( return lat, lon, height
obs2.parse_file(file2) # set central point if options.central_point is not None: cpoint = options.central_point.split(",") if len(cpoint) is not 2: print >>sys.stderr, "Cannot read central point: %s" % \ options.central_point sys.exit(1) print >>sys.stderr, "Local system: Central point: lat=%s, lon=%s deg" %\ (cpoint[0], cpoint[1]) from gizela.util.Converter import Converter lat = Converter.deg2rad_(float(cpoint[0])) lon = Converter.deg2rad_(float(cpoint[1])) obs1.centralPointGeo.lat=lat obs1.centralPointGeo.lon=lon obs2.centralPointGeo.lat=lat obs2.centralPointGeo.lon=lon # transform vector to local system if options.vecl and options.central_point: obs1.tran_vec_local_ne() obs1.set_axes_ori("ne") #print obs1.make_gama_xml() obs2.tran_vec_local_ne() obs2.set_axes_ori("ne")