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)
