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