示例#1
0
 def __init__(self, parent, display):
     logger.debug("Bilding {}".format(self))
     super(SingletViewer, self).__init__(parent, display)
     self.xy_changed = Signal()
示例#2
0
 def __init__(self, parent, display):
     logger.debug("Bilding {}".format(self))
     super(SingletViewer, self).__init__(parent, display)
     self.xy_changed = Signal()
示例#3
0
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)
示例#4
0
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)