def __init__(self, parent, display):
logger.debug("Bilding {}".format(self))
super(SingletViewer, self).__init__(parent, display)
self.xy_changed = Signal()
def __init__(self, parent, display):
logger.debug("Bilding {}".format(self))
super(SingletViewer, self).__init__(parent, display)
self.xy_changed = Signal()
class SingletViewer(WxMPLFitsViewer):
"""
Displays a single FITS image at a time.
"""
def __init__(self, parent, display):
logger.debug("Bilding {}".format(self))
super(SingletViewer, self).__init__(parent, display)
self.xy_changed = Signal()
def _refresh_markers(self, cutout):
self._displayables_by_cutout[cutout].clear_markers(self.ds9)
self.mark_apertures(self.current_cutout)
def mark_apertures(self, cutout, pixel=False):
"""
:param cutout: A SourceCutout that is displayed
:type cutout: SourceCutout
:param pixel: Mark based on pixel locations or based on RA/DEC ?
:type pixel: bool
"""
ra = dec = None
if not self.mark_source:
return
try:
if cutout.reading.mpc_observation is not None:
# Try using the mpc observation attached to this image.
ra = cutout.reading.mpc_observation.coordinate.ra.to('degree').value
dec = cutout.reading.mpc_observation.coordinate.dec.to('degree').value
x, y, hdulist_idx = cutout.world2pix(ra, dec, usepv=True)
(ra, dec) = cutout.pix2world(x, y, hdulist_idx, usepv=False)
else:
ra = cutout.reading.sky_coord.ra.to('degree').value
dec = cutout.reading.sky_coord.dec.to('degree').value
x, y, hdulist_idx = cutout.world2pix(ra, dec, usepv=True)
(ra, dec) = cutout.pix2world(x, y, hdulist_idx, usepv=False)
except Exception as ex:
sys.stderr.write("Exception while computing location to mark aperture: {}\n".format(ex))
try:
# radii = (cutout.apcor.aperture, cutout.apcor.sky, cutout.apcor.swidth+cutout.apcor.sky)
radii = [cutout.apcor.aperture, cutout.apcor.sky, cutout.apcor.swidth+cutout.apcor.sky]
except Exception as ex:
sys.stderr.write("Exception while accessing the aperture correction: {}\n".format(ex))
logger.info("Exception: {0}".format(ex))
logger.warning("Failed trying to get apcor radius values, using defaults for marking.")
radii = [5, 25, 30]
radii = [radius * 0.185 * units.arcsec for radius in radii]
self._displayables_by_cutout[cutout].place_annulus(ra, dec, radii, colour='r')
def place_marker(self, cutout, ra, dec, radius, colour, force=False):
x, y, hdulist_index = cutout.world2pix(ra, dec, usepv=True)
ra, dec = cutout.pix2world(x, y, hdulist_index, usepv=False)
self._displayables_by_cutout[cutout].place_marker(ra, dec, radius, colour=colour, force=force)
def register_xy_changed_event_handler(self, handler):
self.xy_changed.connect(handler)
def _create_displayable(self, cutout):
return DisplayableImageSinglet(cutout.hdulist, self.ds9)
def _attach_handlers(self, displayable):
displayable.xy_changed.connect(self.xy_changed.fire)
displayable.focus_released.connect(self.release_focus)
def _detach_handlers(self, displayable):
if displayable is not None:
displayable.xy_changed.disconnect(self.xy_changed.fire)
displayable.focus_released.disconnect(self.release_focus)
class SingletViewer(WxMPLFitsViewer):
"""
Displays a single FITS image at a time.
"""
def __init__(self, parent, display):
logger.debug("Bilding {}".format(self))
super(SingletViewer, self).__init__(parent, display)
self.xy_changed = Signal()
def _refresh_markers(self, cutout):
self._displayables_by_cutout[cutout].clear_markers(self.ds9)
self.mark_apertures(self.current_cutout)
def mark_apertures(self, cutout, pixel=False):
"""
:param cutout: A SourceCutout that is displayed
:type cutout: SourceCutout
:param pixel: Mark based on pixel locations or based on RA/DEC ?
:type pixel: bool
"""
ra = dec = None
if not self.mark_source:
return
try:
if cutout.reading.mpc_observation is not None:
# Try using the mpc observation attached to this image.
ra = cutout.reading.mpc_observation.coordinate.ra.to('degree').value
dec = cutout.reading.mpc_observation.coordinate.dec.to('degree').value
x, y, hdulist_idx = cutout.world2pix(ra, dec, usepv=True)
(ra, dec) = cutout.pix2world(x, y, hdulist_idx, usepv=False)
else:
ra = cutout.reading.sky_coord.ra.to('degree').value
dec = cutout.reading.sky_coord.dec.to('degree').value
x, y, hdulist_idx = cutout.world2pix(ra, dec, usepv=True)
(ra, dec) = cutout.pix2world(x, y, hdulist_idx, usepv=False)
except Exception as ex:
sys.stderr.write("Exception while computing location to mark aperture: {}\n".format(ex))
try:
# radii = (cutout.apcor.aperture, cutout.apcor.sky, cutout.apcor.swidth+cutout.apcor.sky)
radii = [cutout.apcor.aperture, cutout.apcor.sky, cutout.apcor.swidth+cutout.apcor.sky]
except Exception as ex:
sys.stderr.write("Exception while accessing the aperture correction: {}\n".format(ex))
logger.info("Exception: {0}".format(ex))
logger.warning("Failed trying to get apcor radius values, using defaults for marking.")
radii = [5, 25, 30]
radii = [radius * 0.185 * units.arcsec for radius in radii]
self._displayables_by_cutout[cutout].place_annulus(ra, dec, radii, colour='r')
def place_marker(self, cutout, ra, dec, radius, colour, force=False):
x, y, hdulist_index = cutout.world2pix(ra, dec, usepv=True)
ra, dec = cutout.pix2world(x, y, hdulist_index, usepv=False)
self._displayables_by_cutout[cutout].place_marker(ra, dec, radius, colour=colour, force=force)
def register_xy_changed_event_handler(self, handler):
self.xy_changed.connect(handler)
def _create_displayable(self, cutout):
return DisplayableImageSinglet(cutout.hdulist, self.ds9)
def _attach_handlers(self, displayable):
displayable.xy_changed.connect(self.xy_changed.fire)
displayable.focus_released.connect(self.release_focus)
def _detach_handlers(self, displayable):
if displayable is not None:
displayable.xy_changed.disconnect(self.xy_changed.fire)
displayable.focus_released.disconnect(self.release_focus)
