def plot(self): extra = [] if self.ra is not None and self.dec is not None: extra.append((self.ra, self.dec, "A")) extra.append((self.prec_ra, self.prec_dec, "P")) ra = 0.0 #self.ra dec = 90.0 #self.dec size = 2.0 w = 640 h = 640 pixscale = size / w wcs = Tan(*[float(x) for x in [ ra, dec, 0.5 + (w / 2.), 0.5 + (h / 2.), -pixscale, 0., 0., -pixscale, w, h, ]]) plot = Plotter(wcs) return plot.plot(extra=extra)
# ra, dec = pp.solve()
# extra.append((ra, dec, ""))
#for i in range (20, 50):
# ra, dec = p.compute2(0, i)
# extra.append((ra, dec, ""))
qp = precession()
ra, dec = qp.inv().transform_ra_dec([0, 90])
extra.append((ra, dec, "p"))
print "precession", ra, dec
ra = 0.
dec = 90.
size = 70.0
w = 800
h = 800
pixscale = size / w
wcs = Tan(*[float(x) for x in [
ra, dec, 0.5 + (w / 2.), 0.5 + (h / 2.),
-pixscale, 0., 0., pixscale, w, h,
]])
plot = Plotter(wcs)
cv2.imshow("plot", plot.plot(extra=extra))
print p.zenith()
cv2.waitKey(0)
