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 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
