def __init__(self, id, dx=None, dy=None, dz=None, index=None, covmat=None, textTable=None): """ dx, dy, dz: displacements covmatd: covariance matrix of displacement """ if textTable == None: textTable = displ_text_table() self.displ = PointCartCovMat(id=id, x=dx, y=dy, z=dz, index=index, covmat=covmat, textTable=textTable)
class PointDisplBase(object): """ base class for classes of points with displacements """ def __init__(self, id, dx=None, dy=None, dz=None, index=None, covmat=None, textTable=None): """ dx, dy, dz: displacements covmatd: covariance matrix of displacement """ if textTable == None: textTable = displ_text_table() self.displ = PointCartCovMat(id=id, x=dx, y=dy, z=dz, index=index, covmat=covmat, textTable=textTable) def set_displacement(self, displ): """ displ: PointCartCovMat instance with displacements """ if not isinstance(displ, PointCartCovMat): raise PointDisplBaseError, "PointCartCovMat instance expected" # set id displ.id = self.displ.id # set textTable #displ.textTable = self.displ.textTable self.displ = displ def __str__(self): return self.displ.__str__()
# line line, = figure.gca().plot([x, x], [y - stdev, y + stdev]) line.set_solid_capstyle("butt") figure.set_style(style, line) if __name__ == "__main__": from gizela.pyplot.FigureLayoutBase import FigureLayoutBase fig = FigureLayoutBase() from gizela.data.PointCartCovMat import PointCartCovMat p = PointCartCovMat(id="A", x=0, y=0) p.var = (9, 4) print p.errEll style = {} PlotPoint.plot_point_xy(fig, p.x, p.y, style) PlotPoint.plot_label_xy(fig, p.id, p.x, p.y, style) PlotPoint.plot_error_ellipse_xy(fig, p.x, p.y, p.errEll, style) PlotPoint.plot_vector_xy(fig, [p.x, p.x + 3, p.x + 1.0], [p.y, p.y - 2, p.y - 1.5], style, style) fig.gca().axis([-4,4,-4,4]) fig.set_aspect_equal() #fig.show_() ori = ("en", "wn", "ne", "nw", "es", "ws", "se", "sw")
point.plot_error_ellipse(figure) def set_use_apriori(self, use): """ sets use of apriori/aposteriori standard deviation use: True/False """ for point in self: point.covmat.useApriori=use if __name__ == "__main__": # points c1 = PointCartCovMat(id="A",x=1,y=2,z=3) c1.var = (3,2,3) print c1 c2 = PointCartCovMat(id="B",x=3,y=5,z=6) c2.var = (9,9,9) c3 = PointCartCovMat(id="C",x=2,y=4) c4 = PointCartCovMat(id="D",z=6) # covariance matrix covmat = CovMat(9,8) for i in xrange(9,0,-1): for j in xrange(i): covmat.append_value((j+10.0-i)*1e-5) print covmat.data # point list pd=PointListCovMat()
def _handle_start_element(self, name, attrs): '''start element handler for method all''' if not self._root: if name == "gama-local": self._root = True else: raise GamaLocalObsParserError, "<gama-local> root element not found" self._tagLLast=self._tagLast self._tagLast=self._tagName self._tagName=name if name=="network" or\ name=="parameters" or\ name=="points-observations": # save attributes if attrs: for aname, val in attrs.items(): #import sys #print >>sys.stderr, "Parameter: %s" % aname #parameters if aname=="update-constrained-coordinates": self.data.param[aname] = val elif aname=="epoch": self.data.param[aname] = float(val) elif aname=="tol-abs": self.data.param[aname] = float(val) / Unit.tolabs # standard deviation elif aname=="direction-stdev" or\ aname=="angle-stdev" or\ aname=="zenith-angle-stdev": self.data.param[aname] = \ float(val) / Unit.angleStdev elif aname=="distance-stdev": val = [float(v) / Unit.distStdev \ for v in val.split()] self.data.param[aname] = val elif aname=="sigma-apr": self.data.stdev["apriori"] = float(val) #self.data.stdev.apriori = float(val) elif aname=="sigma-act": self.data.stdev["used"] = val.strip() #if val.strip() == "apriori": # self.data.stdev.set_use_apriori() #elif val.strip() == "aposteriori": # self.data.stdev.set_use_aposteriori() #else: # raise GamaLocalObsParserError, \ # "unknown sigma-act tag: %s" %val elif aname=="conf-pr": self.data.stdev["probability"] = float(val) #self.data.stdev.set_conf_prob(float(val)) elif aname=="axes-xy": self.data.param[aname] = val #self.data.coordSystemLocal.axesOri = val elif aname=="angles": self.data.param[aname] = val #self.data.coordSystemLocal.bearingOri = val else: import sys sys.stderr.write("unhandled parameter %s = \"%s\"\n" %\ (aname, val)) # points and coordinates elif name=="point": x, y, z, adj, fix = None, None, None, None, None if "x" in attrs: x = float(attrs["x"]) if "y" in attrs: y = float(attrs["y"]) if "z" in attrs: z = float(attrs["z"]) if self._coor: # inside tag <coordinates> point = PointCartCovMat(id = attrs['id'], x=x, y=y,\ z=z) self._cluster.append_obs(point) else: point = PointLocalGama(id = attrs['id'], x=x, y=y, z=z) if "adj" in attrs: point.status = gama_coord_status_attr_adj(attrs["adj"]) self._pointListAdj.add_point(point) elif "fix" in attrs: point.status = gama_coord_status_attr_fix(attrs["fix"]) self._pointListFix.add_point(point) # cluster obs and observations elif name=="obs": if "from" in attrs: fromid = attrs["from"] else: fromid=None if "from_dh" in attrs: fromdh = float(attrs["from_dh"]) else: fromdh=None if "orientation" in attrs: import sys print >>sys.stderr, 'warning: attribute "orientation" not accepted' self._cluster = ObsCluster(fromid=fromid, fromdh=fromdh) self._clusterList.append_cluster(self._cluster) elif name=="direction": if "stdev" in attrs: stdev = float(attrs["stdev"] ) / Unit.angleStdev else: stdev=None if "from_dh" in attrs: fromdh = float(attrs["from_dh"]) else: fromdh=None if "to_dh" in attrs: todh = float(attrs["to_dh"]) else: todh=None val = float(attrs["val"]) / Unit.angleVal direction = ObsDirection(toid=attrs["to"], val=val,\ stdev=stdev, fromdh=fromdh, todh=todh) self._cluster.append_obs(direction) elif name=="distance": if "from" in attrs: fromid = attrs["from"] else: fromid=None if "stdev" in attrs: stdev = float(attrs["stdev"]) / Unit.distStdev else: stdev=None if "from_dh" in attrs: fromdh = float(attrs["from_dh"]) else: fromdh=None if "to_dh" in attrs: todh = float(attrs["to_dh"]) else: todh=None distance = ObsDistance(toid=attrs["to"], val=float(attrs["val"]),\ fromid=fromid, stdev=stdev, fromdh=fromdh, todh=todh) self._cluster.append_obs(distance) elif name=="s-distance": if "from" in attrs: fromid = attrs["from"] else: fromid=None if "stdev" in attrs: stdev = float(attrs["stdev"]) / Unit.distStdev else: stdev=None if "from_dh" in attrs: fromdh = float(attrs["from_dh"]) else: fromdh=None if "to_dh" in attrs: todh = float(attrs["to_dh"]) else: todh=None sdistance = ObsSDistance(toid=attrs["to"], val=float(attrs["val"]),\ fromid=fromid, stdev=stdev, fromdh=fromdh, todh=todh) self._cluster.append_obs(sdistance) elif name=="z-angle": if "from" in attrs: fromid = attrs["from"] else: fromid=None if "stdev" in attrs: stdev = float(attrs["stdev"] ) / Unit.angleStdev else: stdev=None if "from_dh" in attrs: fromdh = float(attrs["from_dh"]) else: fromdh=None if "to_dh" in attrs: todh = float(attrs["to_dh"]) else: todh=None val = float(attrs["val"]) / Unit.angleVal zangle = ObsZAngle(toid=attrs["to"], val=val,\ fromid=fromid, stdev=stdev, fromdh=fromdh, todh=todh) self._cluster.append_obs(zangle) # cluster height-differences and observations elif name=="height-differences": self._cluster = ObsClusterHeightDiff() self._clusterList.append_cluster(self._cluster) elif name=="dh": if "stdev" in attrs: stdev = float(attrs["stdev"]) / Unit.distStdev else: stdev=None if "dist" in attrs: dist = float(attrs["dist"]) / Unit.dhdist else: dist=None dh = ObsHeightDiff(fromid=attrs["from"], toid=attrs["to"], val=float(attrs["val"]), stdev=stdev, dist=dist) self._cluster.append_obs(dh) # cluster vector and observations elif name=="vectors": self._cluster = ObsClusterVector() #self.data["clusterList"].append_cluster(cluster) self._clusterList.append_cluster(self._cluster) elif name=="vec": if "from_dh" in attrs: fromdh = float(attrs["from_dh"]) else: fromdh=None if "to_dh" in attrs: todh = float(attrs["to_dh"]) else: todh=None vec = ObsVector(fromid=attrs["from"], toid=attrs["to"], dx=float(attrs["dx"]), dy=float(attrs["dy"]), dz=float(attrs["dz"])) #fromdh=fromdh, todh=todh) self._cluster.append_obs(vec) # cluster coordinates - observations <point /> are handled in "point" elif name=="coordinates": self._coor = True self._cluster = ObsClusterCoor() self._clusterList.append_cluster(self._cluster) # covariance matrix inside cluster elif name=="cov-mat": self._cluster.covmat = CovMat(dim=int(attrs["dim"]), band=int(attrs["band"])) # tags with nothing to do elif name=="gama-local" or \ name=="description": pass else: import sys sys.stderr.write("unhandled start tag %s\n" % name)