@type file: File handler
"""
parser = GamaLocalAdjParser(self)
parser.parse_file(file)
if __name__ == "__main__":
try:
file = open("../../example/xml-epoch/epoch0.adj.xml")
except Exception, e:
print e
print "try to run make in ../../example/xml-epoch/ directory"
import sys
sys.exit()
adj = GamaLocalDataAdj()
adj.parse_file(file)
print adj
print adj.pointListAdjCovMat.covmat.make_gama_xml()
# graph
try:
from gizela.pyplot.FigureLayoutErrEll import FigureLayoutErrEll
except:
print "import of gizela.pyplot.FigureLayoutBase failed"
else:
fig = FigureLayoutErrEll(errScale=2e3)
adj.plot_point(fig)
fig.show_()
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)
c1.plot_error_ellipse(fig)
c1.plot_error_z(fig)
fig.set_aspect_equal()
#fig.show_()
# addition and subtraction of points
c2 = PointCartCovMat("Q", x=10, y=20, z=30)
c2.var = (1, 1, 1)
c2.cov = (-1, -1, -1)
add = c1 + c2
sub = c1 - c2
add.textTable = coor_cov_table()