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)