def redo(self): obj = self.canvas.getObjectByTag(self.areatag) if obj.kind != 'rectangle': self.stop() return True try: image = self.fitsimage.get_image() # calculate center of bbox wd = obj.x2 - obj.x1 dw = wd // 2 ht = obj.y2 - obj.y1 dh = ht // 2 ctr_x, ctr_y = obj.x1 + dw, obj.y1 + dh ra_ctr, dec_ctr = image.pixtoradec(ctr_x, ctr_y, format='str') # Calculate RA and DEC for the three points # origination point ra_org, dec_org = image.pixtoradec(obj.x1, obj.y1) # destination point ra_dst, dec_dst = image.pixtoradec(obj.x2, obj.y2) # "heel" point making a right triangle ra_heel, dec_heel = image.pixtoradec(obj.x1, obj.y2) ht_deg = image.deltaStarsRaDecDeg(ra_org, dec_org, ra_heel, dec_heel) wd_deg = image.deltaStarsRaDecDeg(ra_heel, dec_heel, ra_dst, dec_dst) radius_deg = image.deltaStarsRaDecDeg(ra_heel, dec_heel, ra_dst, dec_dst) # width and height are specified in arcmin sgn, deg, mn, sec = wcs.degToDms(wd_deg) wd = deg * 60.0 + float(mn) + sec / 60.0 sgn, deg, mn, sec = wcs.degToDms(ht_deg) ht = deg * 60.0 + float(mn) + sec / 60.0 sgn, deg, mn, sec = wcs.degToDms(radius_deg) radius = deg * 60.0 + float(mn) + sec / 60.0 except Exception, e: self.fv.showStatus('BAD WCS: %s' % str(e)) return True
def get_starsep_RaDecDeg(self, ra1_deg, dec1_deg, ra2_deg, dec2_deg): sep = self.deltaStarsRaDecDeg(ra1_deg, dec1_deg, ra2_deg, dec2_deg) ## self.logger.debug("sep=%.3f ra1=%f dec1=%f ra2=%f dec2=%f" % ( ## sep, ra1_deg, dec1_deg, ra2_deg, dec2_deg)) sgn, deg, mn, sec = wcs.degToDms(sep) if deg != 0: txt = '%02d:%02d:%06.3f' % (deg, mn, sec) else: txt = '%02d:%06.3f' % (mn, sec) return txt
def redo(self): obj = self.canvas.getObjectByTag(self.areatag) if obj.kind != 'rectangle': self.stop() return True try: image = self.fitsimage.get_image() # calculate center of bbox wd = obj.x2 - obj.x1 dw = wd // 2 ht = obj.y2 - obj.y1 dh = ht // 2 ctr_x, ctr_y = obj.x1 + dw, obj.y1 + dh ra_ctr, dec_ctr = image.pixtoradec(ctr_x, ctr_y, format='str') # Calculate RA and DEC for the three points # origination point ra_org, dec_org = image.pixtoradec(obj.x1, obj.y1) # destination point ra_dst, dec_dst = image.pixtoradec(obj.x2, obj.y2) # "heel" point making a right triangle ra_heel, dec_heel = image.pixtoradec(obj.x1, obj.y2) ht_deg = image.deltaStarsRaDecDeg(ra_org, dec_org, ra_heel, dec_heel) wd_deg = image.deltaStarsRaDecDeg(ra_heel, dec_heel, ra_dst, dec_dst) radius_deg = image.deltaStarsRaDecDeg(ra_heel, dec_heel, ra_dst, dec_dst) # width and height are specified in arcmin sgn, deg, mn, sec = wcs.degToDms(wd_deg) wd = deg*60.0 + float(mn) + sec/60.0 sgn, deg, mn, sec = wcs.degToDms(ht_deg) ht = deg*60.0 + float(mn) + sec/60.0 sgn, deg, mn, sec = wcs.degToDms(radius_deg) radius = deg*60.0 + float(mn) + sec/60.0 #wd, ht, radius = wd_deg, ht_deg, radius_deg except Exception, e: self.fv.showStatus('BAD WCS: %s' % str(e)) return True