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
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
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)
