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)