Exemple #1
0
    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()
Exemple #2
0
    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
Exemple #3
0
    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
Exemple #4
0
    

if __name__ == "__main__":

    c = PointLocalGama(id="A", x=10, y=20, z=30, status=GamaCoordStatus.fix_xyz)
    c.varx =     1.0
    c.vary =     2.0
    c.varz =     3.0
    c.covxy =   -.10
    c.covxz =   -.20
    c.covyz =   -.30

    print c
    
    from gizela.data.point_text_table import gama_coor_table
    c.textTable = gama_coor_table()
    print c
    
    from gizela.data.point_text_table import gama_coor_var_table
    c.textTable = gama_coor_var_table()
    print c
    
    from gizela.data.point_text_table import gama_coor_cov_table
    c.textTable = gama_coor_cov_table()
    print c
    
    # addition and subrtaction
    c2 = PointLocalGama(id="B", x=110, y=120, z=130,
                        status=GamaCoordStatus.adj_XYz)
    c2.var = (1,2,3)
    c2.cov = (.1,.2,.3)
Exemple #5
0
 def __init__(self):
     super(GamaLocalDataObs, self).__init__()
     self.pointListAdj.textTable = gama_coor_table()