Пример #1
0
    def __init__(self, id, 
                 x=None, y=None, z=None, 
                 covmat=None, 
                 index=None,
                 textTable=None):
        """x, y, z ... coordinates
        covmat ... covariance matrix
        index  ... (xi, yi, zi) indexes of rows in covariance matrix
        textTable ... TextTable instance
        """
        
        if not textTable:
            textTable = coor_var_table()

        super(PointCartCovMat, self).__init__(id, x, y, z, textTable=textTable)
            
        if covmat == None:
            self.covmat = CovMat(3, 2)
        else:
            self.covmat = covmat
        
        i = 0
        if index == None:
            if x == None: self.xi = None
            else: self.xi = i; i += 1
            if y == None: self.yi = None
            else: self.yi = i; i += 1
            if z == None: self.zi = None
            else: self.zi = i
            #self.xi, self.yi, self.zi = 0, 1, 2
        else:
            self.index = index
Пример #2
0
    def __init__(self, covmat=None, textTable=None, duplicateId=DUPLICATE_ID.error, sort=False):
        '''
        duplicateId ... what to do with duplicit point
        textTable ... format of table for text output
        sort   ... sort output by id?
        '''

        if textTable == None:
            textTable = coor_var_table()
        super(PointListCovMat, self).__init__(textTable=textTable, 
                                              duplicateId=duplicateId,
                                              sort=sort)

        if covmat == None:
            self.covmat = CovMat(dim=0, band=-1)
        else:
            self.covmat = covmat

        self._numCoord = 0 # covmat row index of last coordinate
Пример #3
0
    print c1
    print c1.var
    print c1.stdev

    c1.textTable = TextTable([("id","%12s"), ("x","%11.3f"), ("y","%11.3f")], type="plain")
    print c1

    # output tables
    from gizela.data.point_text_table import *
    c1.textTable = coor_table()
    print c1
    
    c1.textTable = coor_stdev_table()
    print c1
    
    c1.textTable = coor_var_table()
    print c1

    c1.textTable = coor_cov_table()
    print c1

    print c1.make_gama_xml()

    # error ellipse
    print c1.errEll
    print c1.stdevz
    
    # graph
    from gizela.pyplot.FigureLayoutErrEll import FigureLayoutErrEll
    fig = FigureLayoutErrEll()
    c1.plot_(fig)