LSST Applications 27.0.0,g0265f82a02+469cd937ee,g02d81e74bb+21ad69e7e1,g1470d8bcf6+cbe83ee85a,g2079a07aa2+e67c6346a6,g212a7c68fe+04a9158687,g2305ad1205+94392ce272,g295015adf3+81dd352a9d,g2bbee38e9b+469cd937ee,g337abbeb29+469cd937ee,g3939d97d7f+72a9f7b576,g487adcacf7+71499e7cba,g50ff169b8f+5929b3527e,g52b1c1532d+a6fc98d2e7,g591dd9f2cf+df404f777f,g5a732f18d5+be83d3ecdb,g64a986408d+21ad69e7e1,g858d7b2824+21ad69e7e1,g8a8a8dda67+a6fc98d2e7,g99cad8db69+f62e5b0af5,g9ddcbc5298+d4bad12328,ga1e77700b3+9c366c4306,ga8c6da7877+71e4819109,gb0e22166c9+25ba2f69a1,gb6a65358fc+469cd937ee,gbb8dafda3b+69d3c0e320,gc07e1c2157+a98bf949bb,gc120e1dc64+615ec43309,gc28159a63d+469cd937ee,gcf0d15dbbd+72a9f7b576,gdaeeff99f8+a38ce5ea23,ge6526c86ff+3a7c1ac5f1,ge79ae78c31+469cd937ee,gee10cc3b42+a6fc98d2e7,gf1cff7945b+21ad69e7e1,gfbcc870c63+9a11dc8c8f
LSST Data Management Base Package
Loading...
Searching...
No Matches
cameraConfig.py
Go to the documentation of this file.
1# This file is part of afw.
2#
3# Developed for the LSST Data Management System.
4# This product includes software developed by the LSST Project
5# (https://www.lsst.org).
6# See the COPYRIGHT file at the top-level directory of this distribution
7# for details of code ownership.
8#
9# This program is free software: you can redistribute it and/or modify
10# it under the terms of the GNU General Public License as published by
11# the Free Software Foundation, either version 3 of the License, or
12# (at your option) any later version.
13#
14# This program is distributed in the hope that it will be useful,
15# but WITHOUT ANY WARRANTY; without even the implied warranty of
16# MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
17# GNU General Public License for more details.
18#
19# You should have received a copy of the GNU General Public License
20# along with this program. If not, see <https://www.gnu.org/licenses/>.
21
22__all__ = ["CameraConfig", "DetectorConfig"]
23
24import numpy as np
25import lsst.pex.config as pexConfig
26import lsst.geom as geom
27from ._cameraGeom import Orientation
28from ._transformConfig import TransformMapConfig
29
30
31class DetectorConfig(pexConfig.Config):
32 """A configuration that represents (and can be used to construct) a
33 Detector.
34 """
35 transformDict = pexConfig.ConfigField(
36 "Dictionary of camera transforms keyed on the transform type.", TransformMapConfig)
37 name = pexConfig.Field("Name of detector slot", str)
38 id = pexConfig.Field("ID of detector slot", int)
39 bbox_x0 = pexConfig.Field("x0 of pixel bounding box", int)
40 bbox_y0 = pexConfig.Field("y0 of pixel bounding box", int)
41 bbox_x1 = pexConfig.Field("x1 of pixel bounding box", int)
42 bbox_y1 = pexConfig.Field("y1 of pixel bounding box", int)
43 detectorType = pexConfig.Field(
44 "Detector type: SCIENCE=0, FOCUS=1, GUIDER=2, WAVEFRONT=3", int)
45 physicalType = pexConfig.Field(
46 "How this specific detector is constructed; e.g. CCD, E2V, HgCdTe ", str, default="CCD")
47 serial = pexConfig.Field(
48 "Serial string associated with this specific detector", str)
49 offset_x = pexConfig.Field(
50 "x offset from the origin of the camera in mm in the transposed system.", float)
51 offset_y = pexConfig.Field(
52 "y offset from the origin of the camera in mm in the transposed system.", float)
53 offset_z = pexConfig.Field(
54 "z offset from the origin of the camera in mm in the transposed system.", float, default=0.0)
55 refpos_x = pexConfig.Field("x position of the reference point in the detector in pixels "
56 "in transposed coordinates.", float)
57 refpos_y = pexConfig.Field("y position of the reference point in the detector in pixels "
58 "in transposed coordinates.", float)
59 yawDeg = pexConfig.Field("yaw (rotation about z) of the detector in degrees. "
60 "This includes any necessary rotation to go from "
61 "detector coordinates to camera coordinates "
62 "after optional transposition.", float)
63 pitchDeg = pexConfig.Field(
64 "pitch (rotation about y) of the detector in degrees", float)
65 rollDeg = pexConfig.Field(
66 "roll (rotation about x) of the detector in degrees", float)
67 pixelSize_x = pexConfig.Field("Pixel size in the x dimension in mm", float)
68 pixelSize_y = pexConfig.Field("Pixel size in the y dimension in mm", float)
69
70 # Depending on the choice of detector coordinates, the pixel grid may need
71 # to be transposed before rotation to put it in camera coordinates.
72 transposeDetector = pexConfig.Field(
73 "Transpose the pixel grid before orienting in focal plane?", bool)
74
75 crosstalk = pexConfig.ListField(
76 dtype=float,
77 doc=("Flattened crosstalk coefficient matrix; should have nAmps x nAmps entries. "
78 "Once 'reshape'-ed, ``coeffs[i][j]`` is the fraction of the j-th amp present on the i-th amp."),
79 optional=True
80 )
81
82 # Accessors to get "compiled" versions of parameters.
83 def getCrosstalk(self, numAmps):
84 """Return a 2-D numpy array of crosstalk coefficients of the proper shape"""
85 if not self.crosstalk:
86 return None
87
88 if numAmps != int(np.sqrt(len(self.crosstalk))):
89 numAmps = int(np.sqrt(len(self.crosstalk)))
90 try:
91 return np.array(self.crosstalk, dtype=np.float32).reshape((numAmps, numAmps))
92 except Exception as e:
93 raise RuntimeError(f"Cannot reshape 'crosstalk' coefficients to square matrix: {e}")
94
95 @property
96 def bbox(self):
97 """Return the detector bounding box from the separate box endpoint
98 values.
99 """
100 return geom.BoxI(geom.PointI(self.bbox_x0, self.bbox_y0),
101 geom.PointI(self.bbox_x1, self.bbox_y1))
102
103 @property
104 def offset(self):
105 """Return the detector offset as a Point3D from the separate config
106 values.
107 """
108 return geom.Point3D(self.offset_x, self.offset_y, self.offset_z)
109
110 @property
111 def refPos(self):
112 """Return the detector reference position as a Point2D from the
113 separate config values.
114 """
115 return geom.Point2D(self.refpos_x, self.refpos_y)
116
117 @property
118 def orientation(self):
119 """Return the cameraGeom.Orientation() object defined by the
120 configuration values.
121 """
122 return Orientation(self.offsetoffset, self.refPosrefPos,
123 geom.Angle(self.yawDeg, geom.degrees),
124 geom.Angle(self.pitchDeg, geom.degrees),
125 geom.Angle(self.rollDeg, geom.degrees))
126
127 @property
128 def pixelSize(self):
129 """Return the pixel size as an Extent2D from the separate values.
130 """
131 return geom.Extent2D(self.pixelSize_x, self.pixelSize_y)
132
133
134class CameraConfig(pexConfig.Config):
135 """A configuration that represents (and can be used to construct) a Camera.
136 """
137 detectorList = pexConfig.ConfigDictField(
138 "List of detector configs", keytype=int, itemtype=DetectorConfig)
139 transformDict = pexConfig.ConfigField(
140 "Dictionary of camera transforms keyed on the transform type.", TransformMapConfig)
141 name = pexConfig.Field("Name of this camera", str)
142
143 plateScale = pexConfig.Field(
144 "Plate scale of the camera in arcsec/mm", float)
145 # Note that the radial transform will also apply a scaling, so all coefficients should be
146 # scaled by the plate scale in appropriate units
147 radialCoeffs = pexConfig.ListField(
148 "Coefficients for radial distortion", float)
AmpInfoBoxKey bbox
Definition Amplifier.cc:117
table::Point2DKey pixelSize
Definition Detector.cc:166
A class representing an angle.
Definition Angle.h:128
An integer coordinate rectangle.
Definition Box.h:55