if isinstance(obs, ObsZAngle):
print obs, "is zenit angle"
if isinstance(obs, ObsSDistance):
print obs, "is slope distance"
print "Making correction"
from gizela.corr.ObsCorrSphere import ObsCorrSphere
from gizela.util.Ellipsoid import Ellipsoid
from gizela.data.PointGeodetic import PointGeodetic
from gizela.data.PointCart import PointCart
from gizela.data.PointList import PointList
from gizela.data.ObsClusterList import ObsClusterList
pA = PointCart("AB", x=0, y=0, z=0)
pl = PointList(textTable=pA.textTable)
pl.add_point(pA)
pl.add_point(PointCart("CD", x=1, y=1, z=1))
print pl
ocl = ObsClusterList()
ocl.append_cluster(cl)
corr = ObsCorrSphere(ellipsoid=Ellipsoid(),
centralPointGeo=PointGeodetic(id="c", lat=0, lon=0, height=0),
centralPointLoc=PointCart(id="c", x=0, y=0, z=0),
axesOri="ne",
bearingOri="righ-handed",
pointList=pl,
obsClusterList=ocl)
cl.compute_corr(corr)
print cl
print cl.make_gama_xml(corrected=True)
# for point in self.list: point.plot_point(figure, pointProp)
if __name__ == "__main__":
from gizela.data.PointCart import PointCart
from gizela.data.PointList import PointList
c1 = PointCart(id="C",x=1,y=2,z=3)
c2 = PointCart(id="B",x=4,y=5,z=6)
pl1=PointList()
pl2=PointList()
pl1.add_point(c1)
pl1.add_point(c2)
pl2.add_point(c2)
print pl1
print pl2
el = EpochPointList()
el.add_(pl1)
el.add_(pl2)
print el.index
el.sortOutput = True
print el
for p in el.iter_point("C"):
# for point in self.list: point.plot_point(figure, pointProp)
if __name__ == "__main__":
from gizela.data.PointCart import PointCart
from gizela.data.PointList import PointList
c1 = PointCart(id="C",x=1,y=2,z=3)
c2 = PointCart(id="B",x=4,y=5,z=6)
pl1=PointList()
pl2=PointList()
pl1.add_point(c1)
pl1.add_point(c2)
pl2.add_point(c2)
print pl1
print pl2
el = EpochPointList()
el.add_epoch(pl1)
el.add_epoch(pl2)
print el.index
el.sortOutput = True
print el
for p in el.iter_point("C"):
"""
print "This is ObsCorrSphere class instance"
# corr = ObsCorrSphere(ellipsoid, centralPointGeo, centralPointLoc, pointList)
print "Making correction"
from gizela.util.Ellipsoid import Ellipsoid
from gizela.data.PointGeodetic import PointGeodetic
from gizela.data.PointCart import PointCart
from gizela.data.PointList import PointList
centralPointLoc = PointCart("99", 0, 0, 0)
pointList = PointList(textTable=centralPointLoc.textTable)
pointList.add_point(centralPointLoc)
pointList.add_point(PointCart("999", 1000, 000, 0))
pointList.add_point(PointCart("1", 2000, 0, 0))
pointList.add_point(PointCart("2", 2000, 1000, 0))
pointList.add_point(PointCart("3", 1000, 1000, 0))
pointList.add_point(PointCart("4", 0, 1000, 0))
pointList.add_point(PointCart("5", 0, 0, 0))
pointList.add_point(PointCart("6", 0, -1000, 0))
pointList.add_point(PointCart("7", 1000, -1000, 0))
pointList.add_point(PointCart("8", 2000, -1000, 0))
pointList.add_point(PointCart("11", 2000, 0, 1000))
pointList.add_point(PointCart("12", 2000, 1000, 1000))
pointList.add_point(PointCart("13", 1000, 1000, 1000))
pointList.add_point(PointCart("14", 0, 1000, 1000))
pointList.add_point(PointCart("15", 0, 0, 1000))
else:
if gama_coord_status_dim(status) != len(coord):
print >>sys.stderr, "The number of coordinates (%i) do not agree with status '%s'" % (len(coord), gama_coord_status_string(status))
line = file.readline()
line_num+=1
continue
if len(coord) == 1:
point = PointLocalGama(id=id, z=coord[0], status=status)
elif len(coord) == 2:
point = PointLocalGama(id=id, x=coord[0], y=coord[1],
status=status)
elif len(coord) == 3:
point = PointLocalGama(id=id, x=coord[0], y=coord[1],
z=coord[2], status=status)
pList.add_point(point)
pointsId.append("%s(%s)" % (id,
gama_coord_status_string(status)))
line = file.readline()
line_num+=1
print >>sys.stderr, "Point(s) read: " + " ".join(pointsId)
# read observation files
from gizela.data.GamaLocalDataObs import GamaLocalDataObs
from gizela.data.NetworkObsList import NetworkObsList
from gizela.data.Network import Network
obsList = NetworkObsList(coordSystemLocal=localSystem)
