コード例 #1
0
ファイル: imagelog.py プロジェクト: jonathansick/MoAstro
 def chip_footprint_polygon(header, hdulist):
     """Create a Mongo-compatible polygon representing the chip footprint
     in equatorial cordinates. The polygon is a length-4 list, populated
     with length-2 lists of RA, Dec vertices.
     """
     wcs = astropy.wcs.WCS(header=header, fobj=hdulist)
     footprint = wcs.calc_footprint(header=header)
     footprint_lst = footprint.tolist()  # cast as a list of floats
     return footprint_lst
コード例 #2
0
ファイル: imagedb.py プロジェクト: jonathansick/skyoffset
    def add_footprint_from_wcs(self, mosaic_name, wcs):
        """Create a footprint from the WCS embedded in a ``astropy.io.fits``
        header for the named (pre-existing) mosaic document.

        The footprint is stored under the field ``footprint`` in the mosaic's
        document.

        In a sense, this is an attempt to persist a WCS instance. Note that
        we only save a subset of the astropy.wcs data; that is, we're built
        around simple WCS with no PV, etc. This could be fixed though...
        
        .. note:: By default astropy.wcs is 1 based (ie, origins of CRPIX are
           1 and not zero; may need to subtract 1 from crpix when used in
           numpy arrays

        Parameters
        ----------
        mosaic_name : str
            Name of the mosaic (the ``_id`` field for the document).
        wcs : :class:`astropy.wcs.WCS` instance
            WCS for the mosaic.
        """
        doc = {}
        doc['naxis'] = (wcs._naxis1, wcs._naxis2)
        doc['crpix'] = tuple(wcs.wcs.crpix)  # (CRPIX1, CRPIX2)
        doc['crval'] = tuple(wcs.wcs.crval)  # (CRVAL1, CRVAL2)
        doc['ctype'] = tuple(wcs.wcs.ctype)
        if wcs.wcs.has_cd():
            cd = []
            for (cdi, cdj) in wcs.wcs.cd:
                cd.append([cdi, cdj])
            doc['cd'] = cd
        try:
            doc['cdelt'] = tuple(wcs.wcs.cdelt)
        except:
            pass
        try:
            doc['crota'] = tuple(wcs.wcs.crota)
        except:
            pass
        # Make footprint polygon, cast to a list
        raDecFootprintArray = wcs.calc_footprint()
        raDecFootprint = []
        for (ra, dec) in raDecFootprintArray:
            raDecFootprint.append([ra, dec])
        doc['radec_poly'] = raDecFootprint
        # Add footprint to mosaic document
        self.set(mosaic_name, "footprint", doc)
コード例 #3
0
ファイル: fits_utils.py プロジェクト: AstrolabeProject/imdtk
def get_image_corners(wcs):
    """
    Return the image corners as a list, given a valid WCS object, extracted from the image.
    Corners are listed clockwise from the lower left corner.
    """
    return list(wcs.calc_footprint())