self.errScale*self.stdev.get_conf_scale_1d()*point.stdevz,
style="errorZStyle")
if __name__ == "__main__":
try:
file = open("../../example/xml-epoch/epoch.adj.xml")
except Exception, e:
print e
print "try to run make in ../../example/xml-epoch/ directory"
import sys
sys.exit()
from gizela.data.GamaLocalDataAdj import GamaLocalDataAdj
adj = GamaLocalDataAdj()
adj.parse_file(file)
print adj
print adj.pointListAdjCovMat.covmat.make_gama_xml()
print adj.stdev
from gizela.data.Network import Network
from gizela.util.CoordSystemLocal3D import CoordSystemLocal3D
c3d = CoordSystemLocal3D()
net = Network(c3d, adj, useApriori=True)
fig = FigureLayoutErrEll(figScale=0.0001, errScale=2e3)
net.plot_point(fig)
fig.show_()
# graph
# test orientation of axes
sys.exit(1)
#date = [
# datetime.date(2000,1,1),
# datetime.date(2001,1,1),
# datetime.date(2002,1,1),
# datetime.date(2003,1,1),
# datetime.date(2004,1,1),
# datetime.date(2005,1,1)
# ]
epl = EpochList(stdevUseApriori=True)
#epl = EpochList(stdevUseApriori=False)
for f in file:
adj = GamaLocalDataAdj()
adj.parse_file(f)
ep = Epoch(adj)
epl.add_epoch(ep)
print ep.pointListAdj
#print ep.pointListFix
from gizela.stat.DisplacementTestType import DisplacementTestType
#epl.compute_displacement_test()
epl.compute_displacement_test(testType=DisplacementTestType.xy)
print epl
print epl.get_epoch_table()
print epl.pointListFix
from gizela.tran.Tran2D import Tran2D
tran = Tran2D()
tran.set_trn_key(trnKey)
# open files, read data and add them into epochList
from gizela.data.GamaLocalDataAdj import GamaLocalDataAdj
for name in args:
try:
file = open(name)
except IOError, e:
print >>sys.stderr, "Cannot open file '%s'" % name
print e
sys.exit(1)
data = GamaLocalDataAdj()
try:
data.parse_file(file)
except Exception, e:
print >>sys.stderr, "Cannot parse file '%s'" % name
print e
sys.exit(1)
net = Network(coordSystemLocal=localSystem, data=data)
if opt.tranPoint2D is not None:
#print net.pointListAdj
#print net.pointListAdjCovMat
#print net.pointListAdjCovMat.covmat
net.tran_point(tran)
