def __init__(self): self.dateTimeList = [] # list of date and time of observations # joined epoch has more than one value self.description = "" #self.description = ConfigParser.SafeConfigParser() # description as config parser #self.description.optionxform = str # to make options case sensitive #self.fileName = "" # name of file with data #self.coordSystemLocal = CoordSystemLocal3D() # implicit 3D local geodetic coordinates system #self.stdev = StandardDeviation(apriori=1.0, useApriori=True) # standard deviation self.stdev = {"apriori": 1.0, "probability": 0.95, "used": "apriori", "aposteriori": None, "df": None} self.param = {"axes-xy": "ne", "angles" : "right-handed", "epoch" : 0.0, "update-constrained-coordinates" : "no", "tol-abs": 1.0, "direction-stdev": None, "distance-stdev": None, "zenith-angle-stdev": None } self.gamaLocal = {} self.coordinatesSummary = {} self.obsSummary = {} self.projEqn = {} # point list of FIXED POINTS self.pointListFix = PointList(textTable=gama_coor_table()) # list of point TO BE ADJUSTED # point list of adjusted/constrained points without covariance matrix # use for points to be adjusted in input XML document self.pointListAdj = PointList(textTable=gama_coor_table()) # list of ADJUSTED POINTS # point list of adjusted/constrained points with covariance matrix # use for adjusted points in output XML document self.pointListAdjCovMat = PointListCovMat(covmat=CovMatApri(), textTable=gama_coor_stdev_table()) # list of OBSERVATIONS self.obsClusterList = ObsClusterList()
except: print >>sys.stderr, "Can not open file %s" % fn sys.exit(1) if numProt > 1: print >>sys.stderr, "Just one type of protocol for single run supported" sys.exit(1) # text table to use textTable=None if opt.printDisplFile: from gizela.stat.displ_test_text_table import point_displ_text_table textTable = point_displ_text_table() elif opt.printCoordFile: from gizela.data.point_text_table import gama_coor_stdev_table textTable = gama_coor_stdev_table() elif opt.printTestFile: from gizela.stat.displ_test_text_table import displ_test_text_table textTable = displ_test_text_table() #print textTable #print opt.test # scale factors if opt.figScale is not None: opt.figScale = compute_scale(opt.figScale) print >>sys.stderr, "Figure scale set: %s" % opt.figScale if opt.displScale2D is not None: opt.displScale2D = compute_scale(opt.displScale2D) print >>sys.stderr, "XY Displacement scale set: %s" % opt.displScale2D if opt.displScale1D is not None: