def _handle_end_element_all(self, name): '''end element handler for method all''' if name=="description": self.data.description=self._characterData.strip() elif name=="id": self._id=self._characterData.strip() elif name=="x": self._coordX=float(self._characterData) self._coordType += name #if self._tagLLast == "fixed": # self._coordType |= GamaCoordStatus.fix_x #else: # "adjusted" # self._coordType |= GamaCoordStatus.adj_x elif name=="X": self._coordX=float(self._characterData) self._coordType += name #if self._tagLLast == "fixed": # self._coordType |= GamaCoordStatus.fix_x #else: # "adjusted" # self._coordType |= GamaCoordStatus.adj_X elif name=="y": self._coordY=float(self._characterData) self._coordType += name #if self._tagLLast == "fixed": # self._coordType |= GamaCoordStatus.fix_y #else: # "adjusted" # self._coordType |= GamaCoordStatus.adj_y elif name=="Y": self._coordY=float(self._characterData) self._coordType += name #if self._tagLLast == "fixed": # self._coordType |= GamaCoordStatus.fix_y #else: # "adjusted" # self._coordType |= GamaCoordStatus.adj_Y elif name=="z": self._coordZ=float(self._characterData) self._coordType += name #if self._tagLLast == "fixed": # self._coordType |= GamaCoordStatus.fix_z #else: # "adjusted" # self._coordType |= GamaCoordStatus.adj_z elif name=="Z": self._coordZ=float(self._characterData) self._coordType += name #if self._tagLLast == "fixed": # self._coordType |= GamaCoordStatus.fix_z #else: # "adjusted" # self._coordType |= GamaCoordStatus.adj_Z elif name=="point": # add point to PointDict point=PointLocalGama(id=self._id, x=self._coordX, y=self._coordY, z=self._coordZ) if self._tagLast=="adjusted": point.status = gama_coord_status_attr_adj(self._coordType) # adds point with covmat index = [None, None, None] # rows in covmat if self._coordX != None: index[0] = self._covmatrow self._covmatrow += 1 if self._coordY != None: index[1] = self._covmatrow self._covmatrow += 1 if self._coordZ != None: index[2] = self._covmatrow self._covmatrow += 1 point.index = index self.data.pointListAdjCovMat.add_point(point) #elif self._tagLast=="approximate": # # adds point without covmat # point.texttable = text_table_coor() # self.data["pointApp"].add_point(point) # elif self._tagLast=="fixed": point.status = gama_coord_status_attr_fix(self._coordType) # adds point without covmat self.data.pointListFix.add_point(point) # erase data self._coordX, self._coordY, self._coordZ = None, None, None self._coordType, self._coordNum = "", 0 # orientation shifts #elif name=="approx": pass #self._orientationApprox=float(self._characterData) #elif name=="adj": pass #self._orientationAdj=float(self._characterData) #elif name=="orientation": pass # add orientation shift as point with status ori # approx "approximate" # adj "adjusted" #self.data["approximate"].add_ori_with_index(self._id, self._orientationApprox) #self.data["adjusted" ].add_ori_with_index(self._id, self._orientationAdj) # stdev elif name=="apriori": self.data.pointListAdjCovMat.covmat.apriori = float(self._characterData) elif name=="aposteriori": self.data.pointListAdjCovMat.covmat.aposteriori = float(self._characterData) elif name=="used": self._characterData = self._characterData.strip() if self._characterData == "apriori": self.data.pointListAdjCovMat.covmat.useApriori = True elif self._characterData == "aposteriori": self.data.pointListAdjCovMat.covmat.useApriori = False else: print "*%s*" % self._characterData raise GamaLocalAdjParserError, "unknown data in tag <used>" elif name=="probability": self.data.stdev["probability"] = float(self._characterData) elif name=="degrees-of-freedom": self.data.stdev["df"] = int(self._characterData) # project equations elif name=="equations" or\ name=="unknowns" or\ name=="defect": self.data.projEqn[name] = int(self._characterData) elif name=="sum-of-squares": self.data.projEqn[name] = float(self._characterData) # observation summary elif name=="distances" or\ name=="directions" or\ name=="angles" or\ name=="xyz-coords" or\ name=="h-diffs" or\ name=="z-angles" or\ name=="s-dists" or\ name=="vectors": self.data.obsSummary[name] = int(self._characterData) # dim and band of covariance matrix elif name=="dim": self.data.pointListAdjCovMat.covmat.dim = int(self._characterData) elif name=="band": self.data.pointListAdjCovMat.covmat.band = int(self._characterData) # covariance matrix float numbers elif name=="flt": self.data.pointListAdjCovMat.covmat.append_value(\ float(self._characterData)/(Unit.distStdev**2)) # cut of orientation shifts in covariance matrix elif name=="cov-mat": self.data.pointListAdjCovMat.covmat.slice(self.data.pointListAdjCovMat.get_num_coord()) #for p in self.data.pointListAdjCovMat: # print p.id, ">", p.index # print p.id, "vc", p.var, p.cov self._characterData="" self._tagName=self._tagLast self._tagLast=self._tagLLast
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)