def __init__(self, id,
x=None, y=None, z=None,
status=GamaCoordStatus.unused,
covmat=None, index=None,
dx=None, dy=None, dz=None,
dcovmat=None, dindex=None,
testStat=None,
testPassed=None,
testPValue=None,
testReliability=None,
testType=DisplacementTestType.none,
textTable=None,
epochIndex=None):
if isinstance(id, PointLocalGama):
p = id
id = p.id; x = p.x; y = p.y; z = p.z;
status = p.status; covmat = p.covmat; index = p.index
if textTable == None:
textTable = displ_test_text_table()
PointLocalGama.__init__(self, id=id, x=x, y=y, z=z,
status=status, covmat=covmat, index=index)
PointDisplBase.__init__(self, id=id, dx=dx, dy=dy, dz=dz,
covmat=dcovmat, index=dindex)
TestResult.__init__(self, testStat=testStat, testPassed=testPassed,
testPValue=testPValue,
testReliability=testReliability,
testType=testType, textTable=textTable)
self.epochIndex = epochIndex
def __init__(self,
stdevUseApriori=True,
confProb=0.95,
reliabProb=0.95,
#config=None,
duplicateFixId=DUPLICATE_ID.compare,
textTable=None):
"""
stdevUseApriori: use apriori/aposteriori standard deviation
confProb: confidence probability
reliabProb: reliability probability
config: project configuration file
duplicateFixId: what to do with duplicate fixed points in epoch
textTable: text table instance for epochPointList
"""
self.numEpoch = 0 # the number of epoch added
#self.epochList = []
self.dateTimeList = [] # list of dates of epochs, datetime instances
self.epochPointList = EpochPointList()
# list of coordinates, displacements and test results
# PointLocalGamaDisplTest instances
self.pointListFix = PointList(textTable=gama_coor_table(),
duplicateId=duplicateFixId)
# list of fixed points for all epochs
self.coordSystemLocal = CoordSystemLocal2D()
# definition of local geodetic system
self.stdevList = []
# list of StandardDeviation instances
# handles degrees of freedom from GamaLocalData instance
self._stdevUseApriori = stdevUseApriori
self._confProb = confProb
self.reliabProb = reliabProb
self.testType = None
if textTable is None:
from gizela.stat.displ_test_text_table import displ_test_text_table
self.textTable = displ_test_text_table()
else:
self.textTable = textTable
def __init__(self,
coordSystemLocal,
stdevUseApriori=True,
confProb=0.95,
reliabProb=0.95,
testId=None,
duplicateIdFix=DUPLICATE_ID.compare,
duplicateIdAdj=DUPLICATE_ID.compare,
textTable=None):
"""
stdevUseApriori: use apriori/aposteriori standard deviation
confProb: confidence probability
reliabProb: reliability probability
testId: Ids of point for testing (currently only for output)
duplicateFixId: what to do with duplicate fixed points in epoch
textTable: text table instance for epochPointList
"""
NetworkList.__init__(self,
coordSystemLocal,
stdevUseApriori=stdevUseApriori,
confProb=confProb,
#reliabProb=reliabProb,
#duplicateIdFix=duplicateIdFix,
duplicateIdAdj=duplicateIdAdj)
self.pointListAdj = PointList(textTable=gama_coor_table(),
duplicateId=duplicateIdFix)
# list of fixed points for all epochs
self.epochPointList = EpochPointList(idList=testId)
# list of coordinates, displacements and test results
# PointLocalGamaDisplTest instances
self.stdevList = []
# list of StandardDeviation instances
# handles degrees of freedom from GamaLocalData instance
##self.testType = None
if textTable is None:
from gizela.stat.displ_test_text_table import displ_test_text_table
self._textTable = displ_test_text_table()
else:
self._textTable = textTable
self._reliabProb = reliabProb
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:
opt.displScale1D = compute_scale(opt.displScale1D)
print >>sys.stderr, "Z Displacement scale set: %s" % opt.displScale1D
