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:
