def setUp(self):
self.wcs = WcsSol()
self.ra = 1.0
self.dec = 2.0
rotSkyPos = 0.0
self.wcs.setObsMetaData(self.ra, self.dec, rotSkyPos)
def __init__(self, camera):
"""Initialize the camera data class.
Parameters
----------
camera : lsst.afw.cameraGeom.camera.camera.Camera
A collection of Detectors that also supports coordinate
transformation. (the default is None.)
"""
self._wcs = WcsSol(camera=camera)
# List of wavefront sensor CCD name
self._wfsCcd = []
# Dictionary of (x, y) coordinates of detector corners and dimensions
# of detection. The dictonary key is the ccd name.
self._corners = dict()
self._dimension = dict()
def __init__(self, camera):
"""Initialize the camera data class.
Parameters
----------
camera : lsst.afw.cameraGeom.camera.camera.Camera
A collection of Detectors that also supports coordinate
transformation. (the default is None.)
"""
self._wcs = WcsSol(camera=camera)
# List of wavefront sensor CCD name
self._wfsCcd = []
# Dictionary of (x, y) coordinates of detector corners and dimensions
# of detection. The dictonary key is the ccd name.
self._corners = dict()
self._dimension = dict()
class CameraData(object):
def __init__(self, camera, centerCcd="R22_S11"):
"""Initialize the camera data class.
Parameters
----------
camera : lsst.afw.cameraGeom.camera.camera.Camera
A collection of Detectors that also supports coordinate
transformation. (the default is None.)
centerCcd : str, optional
Center Ccd on the camera. (The default is "R22_S11").
"""
self._wcs = WcsSol(camera=camera)
# Center detector on camera
self._centerCcd = centerCcd
# List of wavefront sensor CCD name
self._wfsCcd = []
# Dictionary of (x, y) coordinates of detector corners and dimensions
# of detection. The dictonary key is the ccd name.
self._corners = dict()
self._dimension = dict()
def setWfsCcdList(self, wfsCcdList):
"""Set the wavefront sensor (WFS) charge-coupled devide (CCD) list.
Parameters
----------
wfsCcdList : list
Wavefront sensor list (e.g. ["R22_S11", "R22_S10"]).
"""
self._wfsCcd = wfsCcdList
def setWfsCorners(self, wfsCorners):
"""Set the wavefront sensor corner information.
Parameters
----------
wfsCorners : dict
Wavefront corner information. The dictionary key is the CCD name
(e.g. "R22_S11").
"""
self._corners = wfsCorners
def setCcdDims(self, ccdDims):
"""Set the charge-coupled device (CCD) dimenstions.
Parameters
----------
ccdDims : dict
CCD dimensions. The dictionary key is the CCD name (e.g.
"R22_S11").
"""
self._dimension = ccdDims
def getWfsCcdList(self):
"""Get the list of wavefront sensor CCD list.
Returns
-------
list
Wavefront sensor CCD list.
"""
return self._wfsCcd
def getWfsCorner(self, detectorName):
"""Get the wavefront sensor corner.
Parameters
----------
detectorName : str
Detector Name (e.g. "R22_S11").
Returns
-------
tuple
Wavefront sensor corner.
"""
return self._corners[detectorName]
def getCcdDim(self, detectorName):
"""Get the CCD dimension.
Parameters
----------
detectorName : str
Detector name (e.g. "R22_S11").
Returns
-------
tuple
CCD dimension in pixel.
"""
return self._dimension[detectorName]
def _initDetectors(self, detectorType):
"""Initializes the camera detectors.
Parameters
----------
detectorType : lsst.afw.cameraGeom.detector.detector.DetectorType
Detector type.
"""
for detector in self._wcs.getCamera():
if detector.getType() == detectorType:
# Collect the ccd name
detectorName = detector.getName()
self._wfsCcd.append(detectorName)
# Get the detector corners
bbox = detector.getBBox()
xmin = bbox.getMinX()
xmax = bbox.getMaxX()
ymin = bbox.getMinY()
ymax = bbox.getMaxY()
self._corners[detectorName] = (
np.array([xmin, xmin, xmax, xmax]),
np.array([ymin, ymax, ymin, ymax]),
)
# The CCD dimension here is an estimation.
# Based on LCA-13381, there are three types of sensors.
# e2V CCD250: 40.04 mm x 40.96 mm
# STA 4400: 20.00 mm x 40.72 mm
# STA 3800C: 40.00 mm x 40.72 mm
dim1, dim2 = bbox.getDimensions()
self._dimension[detectorName] = (int(dim1), int(dim2))
def setObsMetaData(self, ra, dec, rotSkyPos):
"""Set the observation meta data.
Parameters
----------
ra : float
Pointing ra in degree.
dec : float
Pointing decl in degree.
rotSkyPos : float
The orientation of the telescope in degrees.
"""
self._wcs.setObsMetaData(ra, dec, rotSkyPos, centerCcd=self._centerCcd)
def populatePixelFromRADecl(self, stars):
"""Populates the RAInPixel and DeclInPixel coordinates to the stars.
Parameters
----------
stars : StarData
The stars to populate.
Returns
-------
StarData
The stars with x-, y-pixel data populated.
"""
# Do the shallow copy
populatedStar = copy.copy(stars)
# Populate the pixel data
ra = populatedStar.getRA()
decl = populatedStar.getDecl()
chipName = np.array([populatedStar.getDetector()] * len(ra))
raInPixel, declInPixel = self._wcs.pixelCoordsFromRaDec(
ra, decl, chipName=chipName)
populatedStar.setRaInPixel(raInPixel)
populatedStar.setDeclInPixel(declInPixel)
return populatedStar
def getStarsOnDetector(self, stars, offset):
"""Get the stars on the detector according to the pixel position.
Parameters
----------
stars : StarData
Star information.
offset : float
Offset to the dimension of camera. If the detector dimension is 10
(assume 1-D), the star's position between -offset and 10+offset
will be seem to be on the detector.
Returns
-------
StarData
The stars on the detector.
"""
# Do the shallow copy
starsOnDet = copy.copy(stars)
# Get the index that will keep the data
starsRaInPixel = starsOnDet.getRaInPixel()
starsDeclInPixel = starsOnDet.getDeclInPixel()
ccdDim = self.getCcdDim(starsOnDet.getDetector())
keep = []
for ii in range(len(starsRaInPixel)):
if (-offset <= starsRaInPixel[ii] <= ccdDim[0] + offset
and -offset <= starsDeclInPixel[ii] <= ccdDim[1] + offset):
keep.append(ii)
# Remove the stars that are not on the detector
starsOnDet.setId(self._getKeepItem(starsOnDet.getId(), keep))
starsOnDet.setRA(self._getKeepItem(starsOnDet.getRA(), keep))
starsOnDet.setDecl(self._getKeepItem(starsOnDet.getDecl(), keep))
starsOnDet.setRaInPixel(self._getKeepItem(starsRaInPixel, keep))
starsOnDet.setDeclInPixel(self._getKeepItem(starsDeclInPixel, keep))
for filterType in FilterType:
if filterType != FilterType.REF:
magArray = starsOnDet.getMag(filterType)
if len(magArray) != 0:
starsOnDet.setMag(filterType,
self._getKeepItem(magArray, keep))
return starsOnDet
def _getKeepItem(self, valArray, keep):
"""Get the keep items in array.
Parameters
----------
valArray : list or 1-D numpy.ndarray
Value array.
keep : list[int]
List of keep index.
Returns
-------
list
The array that has only the keep values.
"""
return [valArray[idx] for idx in keep]
def getWavefrontSensor(self):
"""
Get the corners of LSST curvature wavefront sensors in (ra, dec) based
on the wavefront sensor list.
Returns:
[dict] -- (ra, dec) of four corners of each sensor with the name
of sensor as a list. The dictionary key is the sensor name.
"""
return self._getDetectorRaDec(self._wfsCcd)
def _getDetectorRaDec(self, detectorList):
"""Get the (ra, dec) of CCD corners in the detector list.
Parameters
----------
detectorList : list
List of detectors. For example, ["R22_S11", "R22_S01"].
Returns
-------
dict
This method returns a dict of list. The dict is keyed on the name
of the wavefront sensor. The list contains the (RA, Dec)
coordinates of the corners of that sensor (RA, Dec are paired as
tuples). For example, output['R00_SW0'] = [(23.0, -5.0),
(23.1, -5.0), (23.0, -5.1), (23.1, -5.1)] would mean that the
wavefront sensor named 'R00_SW0' has its corners at RA 23,
Dec -5; RA 23.1, Dec -5; RA 23, Dec -5.1; and RA 23.1, Dec -5.1
Coordinates are in degrees.
"""
ra_dec_out = dict()
for detector in detectorList:
coords = self._corners[detector]
xPix = coords[0]
yPix = coords[1]
chipName = np.array([detector] * len(xPix))
ra, dec = self._wcs.raDecFromPixelCoords(xPix, yPix, chipName)
ra_dec_out[detector] = [
(ra[0], dec[0]),
(ra[1], dec[1]),
(ra[2], dec[2]),
(ra[3], dec[3]),
]
return ra_dec_out
class TestWcsSol(unittest.TestCase):
"""Test the ComCam class."""
def setUp(self):
self.wcs = WcsSol()
self.ra = 1.0
self.dec = 2.0
rotSkyPos = 0.0
self.wcs.setObsMetaData(self.ra, self.dec, rotSkyPos)
def testSetAndGetCamera(self):
camera = "FaultCamera"
self.wcs.setCamera(camera)
self.assertEqual(self.wcs.getCamera(), camera)
def testSetObservationMetaData(self):
ra = 10.0
dec = 20.0
rotSkyPos = 30.0
self.wcs.setObsMetaData(ra, dec, rotSkyPos)
self.assertAlmostEqual(self.wcs._obs.pointingRA, ra)
self.assertAlmostEqual(self.wcs._obs.pointingDec, dec)
self.assertAlmostEqual(self.wcs._obs.rotSkyPos, rotSkyPos)
def testRaDecFromPixelCoordsForSingleChip(self):
xPix = 2032
yPix = 2000
chipName = "R:2,2 S:1,1"
raByWcs, decByWcs = self.wcs.raDecFromPixelCoords(xPix, yPix, chipName)
self.assertAlmostEqual(raByWcs, self.ra, places=2)
self.assertAlmostEqual(decByWcs, self.dec, places=2)
def testRaDecFromPixelCoordsForChipArray(self):
xPix = np.array([2032, 2032])
yPix = np.array([2000, 2000])
chipName = np.array(["R:2,2 S:1,1", "R:2,2 S:1,1"])
raByWcs, decByWcs = self.wcs.raDecFromPixelCoords(xPix, yPix, chipName)
self.assertEqual(len(raByWcs), 2)
self.assertEqual(len(decByWcs), 2)
self.assertAlmostEqual(raByWcs[0], self.ra, places=2)
self.assertAlmostEqual(raByWcs[1], self.ra, places=2)
def testPixelCoordsFromRaDecWithoutChipName(self):
xPix, yPix = self.wcs.pixelCoordsFromRaDec(self.ra, self.dec)
self.assertAlmostEqual(xPix, 2032, places=-1)
self.assertAlmostEqual(yPix, 1994, places=-1)
def testPixelCoordsFromRaDecWithChipName(self):
chipName = "R:2,2 S:1,1"
xPix, yPix = self.wcs.pixelCoordsFromRaDec(self.ra, self.dec,
chipName=chipName)
self.assertAlmostEqual(xPix, 2032, places=-1)
self.assertAlmostEqual(yPix, 1994, places=-1)
def testFocalPlaneCoordsFromRaDecWithZeroRot(self):
self.wcs.setObsMetaData(0, 0, 0)
xInMm, yInMm = self.wcs.focalPlaneCoordsFromRaDec(0, 0)
self.assertEqual(xInMm, 0.0)
self.assertEqual(yInMm, 0.0)
# 0.2 arcsec = 10 um => 1 um = 0.02 arcsec => 1 mm = 20 arcsec
# 1 arcsec = 1/3600 degree
xInMm, yInMm = self.wcs.focalPlaneCoordsFromRaDec(20.0/3600, 0)
self.assertAlmostEqual(xInMm, 1.0, places=3)
self.assertAlmostEqual(yInMm, 0.0, places=3)
def testFocalPlaneCoordsFromRaDecWithNonZeroRot(self):
self.wcs.setObsMetaData(0, 0, 45)
# 0.2 arcsec = 10 um => 1 um = 0.02 arcsec => 1 mm = 20 arcsec
# 1 arcsec = 1/3600 degree
xInMm, yInMm = self.wcs.focalPlaneCoordsFromRaDec(20.0/3600, 0)
self.assertAlmostEqual(xInMm, 1/np.sqrt(2), places=3)
self.assertAlmostEqual(yInMm, -1/np.sqrt(2), places=3)
class TestWcsSol(unittest.TestCase):
"""Test the ComCam class."""
def setUp(self):
self.wcs = WcsSol()
self.ra = 1.0
self.dec = 2.0
rotSkyPos = 0.0
self.wcs.setObsMetaData(self.ra, self.dec, rotSkyPos)
def testSetAndGetCamera(self):
camera = "FaultCamera"
self.wcs.setCamera(camera)
self.assertEqual(self.wcs.getCamera(), camera)
def testsetObsMetaData(self):
ra = 10.0
dec = 20.0
rotSkyPos = 30.0
self.wcs.setObsMetaData(ra, dec, rotSkyPos)
skyOrigin = self.wcs.skyWcs.getSkyOrigin()
self.assertAlmostEqual(skyOrigin.getRa().asDegrees(), ra)
self.assertAlmostEqual(skyOrigin.getDec().asDegrees(), dec)
def testFormatCoordListWiIntFloat(self):
ra = 10
dec = 19.0
raOut, decOut = self.wcs._formatCoordList(ra, dec, "ra", "dec")
self.assertEqual(raOut, [ra])
self.assertEqual(decOut, [dec])
def testFormatCoordListWiLists(self):
ra = [10.0, 20.0]
dec = [5.0, 23.0]
raOut, decOut = self.wcs._formatCoordList(ra, dec, "ra", "dec")
self.assertEqual(ra, raOut)
self.assertEqual(dec, decOut)
def testFormatCoordListAssert(self):
ra = [10, 20]
dec = [10.0, 20.0, 30.0]
with self.assertRaises(AssertionError,
msg="Size of ra not same as dec"):
self.wcs._formatCoordList(ra, dec, "ra", "dec")
def testRaDecFromPixelCoordsForSingleChip(self):
xPix = 2032
yPix = 2000
chipName = "R22_S11"
raByWcs, decByWcs = self.wcs.raDecFromPixelCoords(xPix, yPix, chipName)
self.assertAlmostEqual(raByWcs, self.ra, places=2)
self.assertAlmostEqual(decByWcs, self.dec, places=2)
def testRaDecFromPixelCoordsForChipArray(self):
xPix = np.array([2032, 2032])
yPix = np.array([2000, 2000])
chipName = np.array(["R22_S11", "R22_S11"])
raByWcs, decByWcs = self.wcs.raDecFromPixelCoords(xPix, yPix, chipName)
self.assertEqual(len(raByWcs), 2)
self.assertEqual(len(decByWcs), 2)
self.assertAlmostEqual(raByWcs[0], self.ra, places=2)
self.assertAlmostEqual(raByWcs[1], self.ra, places=2)
def testPixelCoordsFromRaDecWithoutChipName(self):
xPix, yPix = self.wcs.pixelCoordsFromRaDec(self.ra, self.dec)
self.assertAlmostEqual(xPix, 2048, places=-1)
self.assertAlmostEqual(yPix, 2000, places=-1)
def testPixelCoordsFromRaDecWithChipName(self):
chipName = "R22_S11"
xPix, yPix = self.wcs.pixelCoordsFromRaDec(self.ra,
self.dec,
chipName=chipName)
self.assertAlmostEqual(xPix, 2048, places=-1)
self.assertAlmostEqual(yPix, 2000, places=-1)
def testPixelCoordsFromRaDecWithWrongType(self):
with self.assertRaises(
ValueError,
msg=
"chipName is an unallowed type. Can be None, string or array of strings.",
):
self.wcs.pixelCoordsFromRaDec(self.ra, self.dec, chipName=20.0)
def testFocalPlaneCoordsFromRaDecWithZeroRot(self):
self.wcs.setObsMetaData(0, 0, 0)
xInMm, yInMm = self.wcs.focalPlaneCoordsFromRaDec(0, 0)
self.assertEqual(xInMm, 0.0)
self.assertEqual(yInMm, 0.0)
# 0.2 arcsec = 10 um => 1 um = 0.02 arcsec => 1 mm = 20 arcsec
# 1 arcsec = 1/3600 degree
xInMm, yInMm = self.wcs.focalPlaneCoordsFromRaDec(20.0 / 3600, 0)
self.assertAlmostEqual(xInMm, 0.0, places=3)
self.assertAlmostEqual(yInMm, 1.0, places=3)
def testFocalPlaneCoordsFromRaDecWithNonZeroRot(self):
self.wcs.setObsMetaData(0, 0, 45)
# 0.2 arcsec = 10 um => 1 um = 0.02 arcsec => 1 mm = 20 arcsec
# 1 arcsec = 1/3600 degree
xInMm, yInMm = self.wcs.focalPlaneCoordsFromRaDec(20.0 / 3600, 0)
self.assertAlmostEqual(xInMm, -1 / np.sqrt(2), places=3)
self.assertAlmostEqual(yInMm, 1 / np.sqrt(2), places=3)
class CameraData(object):
def __init__(self, camera):
"""Initialize the camera data class.
Parameters
----------
camera : lsst.afw.cameraGeom.camera.camera.Camera
A collection of Detectors that also supports coordinate
transformation. (the default is None.)
"""
self._wcs = WcsSol(camera=camera)
# List of wavefront sensor CCD name
self._wfsCcd = []
# Dictionary of (x, y) coordinates of detector corners and dimensions
# of detection. The dictonary key is the ccd name.
self._corners = dict()
self._dimension = dict()
def setWfsCcdList(self, wfsCcdList):
"""Set the wavefront sensor (WFS) charge-coupled devide (CCD) list.
Parameters
----------
wfsCcdList : list
Wavefront sensor list (e.g. ["R:2,2 S:1,1", "R:2,2 S:1,0"]).
"""
self._wfsCcd = wfsCcdList
def setWfsCorners(self, wfsCorners):
"""Set the wavefront sensor corner information.
Parameters
----------
wfsCorners : dict
Wavefront corner information. The dictionary key is the CCD name
(e.g. "R:2,2 S:1,1").
"""
self._corners = wfsCorners
def setCcdDims(self, ccdDims):
"""Set the charge-coupled device (CCD) dimenstions.
Parameters
----------
ccdDims : dict
CCD dimensions. The dictionary key is the CCD name (e.g.
"R:2,2 S:1,1").
"""
self._dimension = ccdDims
def getWfsCcdList(self):
"""Get the list of wavefront sensor CCD list.
Returns
-------
list
Wavefront sensor CCD list.
"""
return self._wfsCcd
def getWfsCorner(self, detectorName):
"""Get the wavefront sensor corner.
Parameters
----------
detectorName : str
Detector Name (e.g. "R:2,2 S:1,1").
Returns
-------
tuple
Wavefront sensor corner.
"""
return self._corners[detectorName]
def getCcdDim(self, detectorName):
"""Get the CCD dimension.
Parameters
----------
detectorName : str
Detector name (e.g. "R:2,2 S:1,1").
Returns
-------
tuple
CCD dimension in pixel.
"""
return self._dimension[detectorName]
def _initDetectors(self, detectorType):
"""Initializes the camera detectors.
Parameters
----------
detectorType : lsst.afw.cameraGeom.detector.detector.DetectorType
Detector type.
"""
for detector in self._wcs.getCamera():
if (detector.getType() == detectorType):
# Collect the ccd name
detectorName = detector.getName()
self._wfsCcd.append(detectorName)
# Get the detector corners
bbox = detector.getBBox()
xmin = bbox.getMinX()
xmax = bbox.getMaxX()
ymin = bbox.getMinY()
ymax = bbox.getMaxY()
self._corners[detectorName] = \
(np.array([xmin, xmin, xmax, xmax]),
np.array([ymin, ymax, ymin, ymax]))
# The CCD dimension here is an estimation.
# Based on LCA-13381, there are three types of sensors.
# e2V CCD250: 40.04 mm x 40.96 mm
# STA 4400: 20.00 mm x 40.72 mm
# STA 3800C: 40.00 mm x 40.72 mm
dim1, dim2 = bbox.getDimensions()
self._dimension[detectorName] = (int(dim1), int(dim2))
def setObsMetaData(self, ra, dec, rotSkyPos, mjd=59580.0):
"""Set the observation meta data.
Parameters
----------
ra : float
Pointing ra in degree.
dec : float
Pointing decl in degree.
rotSkyPos : float
The orientation of the telescope in degrees.
mjd : float
Camera MJD. (the default is 59580.0.)
"""
self._wcs.setObsMetaData(ra, dec, rotSkyPos, mjd=mjd)
def populatePixelFromRADecl(self, stars):
"""Populates the RAInPixel and DeclInPixel coordinates to the stars.
Parameters
----------
stars : StarData
The stars to populate.
Returns
-------
StarData
The stars with x-, y-pixel data populated.
"""
# Do the shallow copy
populatedStar = copy.copy(stars)
# Populate the pixel data
ra = populatedStar.getRA()
decl = populatedStar.getDecl()
chipName = np.array([populatedStar.getDetector()] * len(ra))
raInPixel, declInPixel = self._wcs.pixelCoordsFromRaDec(
ra, decl, chipName=chipName, epoch=2000.0, includeDistortion=True)
populatedStar.setRaInPixel(raInPixel)
populatedStar.setDeclInPixel(declInPixel)
return populatedStar
def getStarsOnDetector(self, stars, offset):
"""Get the stars on the detector according to the pixel position.
Parameters
----------
stars : StarData
Star information.
offset : float
Offset to the dimension of camera. If the detector dimension is 10
(assume 1-D), the star's position between -offset and 10+offset
will be seem to be on the detector.
Returns
-------
StarData
The stars on the detector.
"""
# Do the shallow copy
starsOnDet = copy.copy(stars)
# Get the index that will keep the data
starsRaInPixel = starsOnDet.getRaInPixel()
starsDeclInPixel = starsOnDet.getDeclInPixel()
ccdDim = self.getCcdDim(starsOnDet.getDetector())
keep = []
for ii in range(len(starsRaInPixel)):
if (-offset <= starsRaInPixel[ii] <= ccdDim[0] + offset and
-offset <= starsDeclInPixel[ii] <= ccdDim[1] + offset):
keep.append(ii)
# Remove the stars that are not on the detector
starsOnDet.setId(self._getKeepItem(starsOnDet.getId(), keep))
starsOnDet.setRA(self._getKeepItem(starsOnDet.getRA(), keep))
starsOnDet.setDecl(self._getKeepItem(starsOnDet.getDecl(), keep))
starsOnDet.setRaInPixel(self._getKeepItem(starsRaInPixel, keep))
starsOnDet.setDeclInPixel(self._getKeepItem(starsDeclInPixel, keep))
for filterType in FilterType:
if (filterType != FilterType.REF):
magArray = starsOnDet.getMag(filterType)
if (len(magArray) != 0):
starsOnDet.setMag(filterType, self._getKeepItem(magArray,
keep))
return starsOnDet
def _getKeepItem(self, valArray, keep):
"""Get the keep items in array.
Parameters
----------
valArray : list or 1-D numpy.ndarray
Value array.
keep : list[int]
List of keep index.
Returns
-------
list
The array that has only the keep values.
"""
return [valArray[idx] for idx in keep]
def getWavefrontSensor(self):
"""
Get the corners of LSST curvature wavefront sensors in (ra, dec) based
on the wavefront sensor list.
Returns:
[dict] -- (ra, dec) of four corners of each sensor with the name
of sensor as a list. The dictionary key is the sensor name.
"""
return self._getDetectorRaDec(self._wfsCcd)
def _getDetectorRaDec(self, detectorList):
"""Get the (ra, dec) of CCD corners in the detector list.
Parameters
----------
detectorList : list
List of detectors. For example, ["R:2,2 S:1,1", "R:2,2 S:0,1"].
Returns
-------
dict
This method returns a dict of list. The dict is keyed on the name
of the wavefront sensor. The list contains the (RA, Dec)
coordinates of the corners of that sensor (RA, Dec are paired as
tuples). For example, output['R:0,0 S:2,2B'] = [(23.0, -5.0),
(23.1, -5.0), (23.0, -5.1), (23.1, -5.1)] would mean that the
wavefront sensor named 'R:0,0 S:2,2B' has its corners at RA 23,
Dec -5; RA 23.1, Dec -5; RA 23, Dec -5.1; and RA 23.1, Dec -5.1
Coordinates are in degrees.
"""
ra_dec_out = dict()
for detector in detectorList:
coords = self._corners[detector]
xPix = coords[0]
yPix = coords[1]
chipName = np.array([detector] * len(xPix))
ra, dec = self._wcs.raDecFromPixelCoords(
xPix, yPix, chipName, epoch=2000.0, includeDistortion=True)
ra_dec_out[detector] = [(ra[0], dec[0]), (ra[1], dec[1]),
(ra[2], dec[2]), (ra[3], dec[3])]
return ra_dec_out