24 from functools
import lru_cache
26 import lsst.afw.cameraGeom
as cameraGeom
27 import lsst.geom
as geom
28 import lsst.afw.geom
as afwGeom
29 from lsst.afw.cameraGeom
import Amplifier, Camera, ReadoutCorner
32 __all__ = [
"makeCamera"]
37 """An imaging camera (e.g. the LSST 3Gpix camera)
42 Camera description YAML file.
46 camera : `lsst.afw.cameraGeom.Camera`
50 with open(cameraFile)
as fd:
51 cameraParams = yaml.load(fd, Loader=yaml.CLoader)
53 cameraName = cameraParams[
"name"]
58 plateScale = geom.Angle(cameraParams[
"plateScale"], geom.arcseconds)
59 nativeSys = cameraGeom.CameraSys(cameraParams[
"transforms"].pop(
"nativeSys"))
62 ccdParams = cameraParams[
"CCDs"]
66 for ccdName, ccdValues
in ccdParams.items():
73 """Make a list of detector configs
77 detectorConfig : `list` of `lsst.afw.cameraGeom.DetectorConfig`
78 A list of detector configs.
81 for name, ccd
in ccdParams.items():
82 detectorConfig = cameraGeom.DetectorConfig()
83 detectorConfigs.append(detectorConfig)
85 detectorConfig.name = name
86 detectorConfig.id = ccd[
'id']
87 detectorConfig.serial = ccd[
'serial']
88 detectorConfig.detectorType = ccd[
'detectorType']
89 if 'physicalType' in ccd:
90 detectorConfig.physicalType = ccd[
'physicalType']
92 detectorConfig.bbox_x0, detectorConfig.bbox_y0 = ccd[
'bbox'][0]
93 detectorConfig.bbox_x1, detectorConfig.bbox_y1 = ccd[
'bbox'][1]
94 detectorConfig.pixelSize_x, detectorConfig.pixelSize_y = ccd[
'pixelSize']
95 detectorConfig.transformDict.nativeSys = ccd[
'transformDict'][
'nativeSys']
96 transforms = ccd[
'transformDict'][
'transforms']
97 detectorConfig.transformDict.transforms =
None if transforms ==
'None' else transforms
98 detectorConfig.refpos_x, detectorConfig.refpos_y = ccd[
'refpos']
99 detectorConfig.offset_x, detectorConfig.offset_y = ccd[
'offset']
100 detectorConfig.transposeDetector = ccd[
'transposeDetector']
101 detectorConfig.pitchDeg = ccd[
'pitch']
102 detectorConfig.yawDeg = ccd[
'yaw']
103 detectorConfig.rollDeg = ccd[
'roll']
104 if 'crosstalk' in ccd:
105 detectorConfig.crosstalk = ccd[
'crosstalk']
107 return detectorConfigs
111 """Construct a list of AmplifierBuilder objects
115 amp = list(ccd[
'amplifiers'].values())[0]
118 xRawExtent, yRawExtent = rawBBox.getDimensions()
120 readCorners = {
"LL": ReadoutCorner.LL,
121 "LR": ReadoutCorner.LR,
122 "UL": ReadoutCorner.UL,
123 "UR": ReadoutCorner.UR}
126 for name, amp
in sorted(ccd[
'amplifiers'].items(), key=
lambda x: x[1][
'hdu']):
127 amplifier = Amplifier.Builder()
128 amplifier.setName(name)
131 perAmpData = amp[
'perAmpData']
135 x0, y0 = ix*xRawExtent, iy*yRawExtent
138 xDataExtent, yDataExtent = rawDataBBox.getDimensions()
139 amplifier.setBBox(geom.BoxI(
140 geom.PointI(ix*xDataExtent, iy*yDataExtent), rawDataBBox.getDimensions()))
143 rawBBox.shift(geom.ExtentI(x0, y0))
144 amplifier.setRawBBox(rawBBox)
147 rawDataBBox.shift(geom.ExtentI(x0, y0))
148 amplifier.setRawDataBBox(rawDataBBox)
151 rawSerialOverscanBBox.shift(geom.ExtentI(x0, y0))
152 amplifier.setRawHorizontalOverscanBBox(rawSerialOverscanBBox)
155 rawParallelOverscanBBox.shift(geom.ExtentI(x0, y0))
156 amplifier.setRawVerticalOverscanBBox(rawParallelOverscanBBox)
159 rawSerialPrescanBBox.shift(geom.ExtentI(x0, y0))
160 amplifier.setRawPrescanBBox(rawSerialPrescanBBox)
163 amplifier.setRawXYOffset(geom.Extent2I(ix*xRawExtent, iy*yRawExtent))
165 amplifier.setRawXYOffset(geom.Extent2I(0, 0))
167 amplifier.setReadoutCorner(readCorners[amp[
'readCorner']])
168 amplifier.setGain(amp[
'gain'])
169 amplifier.setReadNoise(amp[
'readNoise'])
170 amplifier.setSaturation(amp[
'saturation'])
171 amplifier.setSuspectLevel(amp.get(
'suspect', np.nan))
174 flipX, flipY = amp.get(
"flipXY")
176 amplifier.setRawFlipX(flipX)
177 amplifier.setRawFlipY(flipY)
179 amplifier.setLinearityCoeffs([float(val)
for val
in amp[
'linearityCoeffs']])
180 amplifier.setLinearityType(amp[
'linearityType'])
181 amplifier.setLinearityThreshold(float(amp[
'linearityThreshold']))
182 amplifier.setLinearityMaximum(float(amp[
'linearityMax']))
183 amplifier.setLinearityUnits(
"DN")
184 amplifierList.append(amplifier)
189 """Backward compatible name.
195 """Given a list [(x0, y0), (xsize, ysize)], probably from a yaml file,
198 (x0, y0), (xsize, ysize) = ylist
199 return geom.BoxI(geom.PointI(x0, y0), geom.ExtentI(xsize, ysize))
203 """Make a dictionary of TransformPoint2ToPoint2s from yaml, mapping from nativeSys
207 nativeSys : `lsst.afw.cameraGeom.CameraSys`
208 transformDict : `dict`
209 A dict specifying parameters of transforms; keys are camera system names.
210 plateScale : `lsst.geom.Angle`
211 The size of a pixel in angular units/mm (e.g. 20 arcsec/mm for LSST)
216 A dict of `lsst.afw.cameraGeom.CameraSys` : `lsst.afw.geom.TransformPoint2ToPoint2`
218 The resulting dict's keys are `~lsst.afw.cameraGeom.CameraSys`,
219 and the values are Transforms *from* NativeSys *to* CameraSys
222 assert nativeSys == cameraGeom.FOCAL_PLANE,
"Cameras with nativeSys != FOCAL_PLANE are not supported."
226 for key, transform
in transformDict.items():
227 transformType = transform[
"transformType"]
228 knownTransformTypes = [
"affine",
"radial"]
229 if transformType
not in knownTransformTypes:
230 raise RuntimeError(
"Saw unknown transform type for %s: %s (known types are: [%s])" % (
231 key, transform[
"transformType"],
", ".join(knownTransformTypes)))
233 if transformType ==
"affine":
234 affine = geom.AffineTransform(np.array(transform[
"linear"]),
235 np.array(transform[
"translation"]))
237 transform = afwGeom.makeTransform(affine)
238 elif transformType ==
"radial":
243 radialCoeffs = np.array(transform[
"coeffs"])
245 radialCoeffs *= plateScale.asRadians()
246 transform = afwGeom.makeRadialTransform(radialCoeffs)
248 raise RuntimeError(
"Impossible condition \"%s\" is not in: [%s])" % (
249 transform[
"transformType"],
", ".join(knownTransformTypes)))
251 resMap[cameraGeom.CameraSys(key)] = transform
257 pupilFactoryClass=cameraGeom.pupil.PupilFactory):
258 """Construct a Camera instance from a dictionary of
259 detector name : `lsst.afw.cameraGeom.amplifier`
264 The name of the camera
265 detectorConfigList : `list`
266 A list of `lsst.afw.cameraGeom.cameraConfig.DetectorConfig`
267 nativeSys : `lsst.afw.cameraGeom.CameraSys`
268 The native transformation type; must be `lsst.afw.cameraGeom.FOCAL_PLANE`
269 transformDict : `dict`
270 A dict of lsst.afw.cameraGeom.CameraSys : `lsst.afw.geom.TransformPoint2ToPoint2`
271 amplifierDict : `dict`
272 A dictionary of detector name :
273 `lsst.afw.cameraGeom.Amplifier.Builder`
274 pupilFactoryClass : `type`, optional
275 Class to attach to camera;
276 `lsst.default afw.cameraGeom.PupilFactory`
280 camera : `lsst.afw.cameraGeom.Camera`
285 Copied from `lsst.afw.cameraGeom.cameraFactory` with permission and encouragement
294 assert nativeSys == cameraGeom.FOCAL_PLANE,
"Cameras with nativeSys != FOCAL_PLANE are not supported."
296 focalPlaneToField = transformDict[cameraGeom.FIELD_ANGLE]
298 cameraBuilder = Camera.Builder(cameraName)
299 cameraBuilder.setPupilFactoryClass(pupilFactoryClass)
302 for toSys, transform
in transformDict.items():
303 cameraBuilder.setTransformFromFocalPlaneTo(toSys, transform)
305 for detectorConfig
in detectorConfigList:
307 cameraGeom.addDetectorBuilderFromConfig(cameraBuilder, detectorConfig,
308 amplifierDict[detectorConfig.name],
315 detectorNativeSys = detectorConfig.transformDict.nativeSys
316 detectorNativeSys = (cameraGeom.PIXELS
if detectorNativeSys
is None else
317 cameraGeom.CameraSysPrefix(detectorNativeSys))
326 assert detectorNativeSys == cameraGeom.PIXELS, \
327 "Detectors with nativeSys != PIXELS are not supported."
328 detectorNativeSys = cameraGeom.CameraSys(detectorNativeSys, detectorConfig.name)
330 return cameraBuilder.finish()