LSST Applications  21.0.0+04719a4bac,21.0.0-1-ga51b5d4+f5e6047307,21.0.0-11-g2b59f77+a9c1acf22d,21.0.0-11-ga42c5b2+86977b0b17,21.0.0-12-gf4ce030+76814010d2,21.0.0-13-g1721dae+760e7a6536,21.0.0-13-g3a573fe+768d78a30a,21.0.0-15-g5a7caf0+f21cbc5713,21.0.0-16-g0fb55c1+b60e2d390c,21.0.0-19-g4cded4ca+71a93a33c0,21.0.0-2-g103fe59+bb20972958,21.0.0-2-g45278ab+04719a4bac,21.0.0-2-g5242d73+3ad5d60fb1,21.0.0-2-g7f82c8f+8babb168e8,21.0.0-2-g8f08a60+06509c8b61,21.0.0-2-g8faa9b5+616205b9df,21.0.0-2-ga326454+8babb168e8,21.0.0-2-gde069b7+5e4aea9c2f,21.0.0-2-gecfae73+1d3a86e577,21.0.0-2-gfc62afb+3ad5d60fb1,21.0.0-25-g1d57be3cd+e73869a214,21.0.0-3-g357aad2+ed88757d29,21.0.0-3-g4a4ce7f+3ad5d60fb1,21.0.0-3-g4be5c26+3ad5d60fb1,21.0.0-3-g65f322c+e0b24896a3,21.0.0-3-g7d9da8d+616205b9df,21.0.0-3-ge02ed75+a9c1acf22d,21.0.0-4-g591bb35+a9c1acf22d,21.0.0-4-g65b4814+b60e2d390c,21.0.0-4-gccdca77+0de219a2bc,21.0.0-4-ge8a399c+6c55c39e83,21.0.0-5-gd00fb1e+05fce91b99,21.0.0-6-gc675373+3ad5d60fb1,21.0.0-64-g1122c245+4fb2b8f86e,21.0.0-7-g04766d7+cd19d05db2,21.0.0-7-gdf92d54+04719a4bac,21.0.0-8-g5674e7b+d1bd76f71f,master-gac4afde19b+a9c1acf22d,w.2021.13
LSST Data Management Base Package
transforms.py
Go to the documentation of this file.
1 # This file is part of meas_base.
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 """Measurement transformations.
23 
24 When a measurement plugin is run, it provides raw, uncalibrated outputs such
25 as pixel positions. A transformation may be run as a post-processing step to
26 convert those outputs to calibrated quantities, such as celestial coordinates.
27 
28 At construction, the transformation is passed the configuration and name of
29 the plugin whose outputs it will be transformaing (all fields in the input
30 table produced by that plugin will have their field names prefixed by the
31 plugin name) and a `~lsst.afw.table.SchemaMapper` which holds the schemata for
32 the input and output catalogs and which may be used to directly map fields
33 between the catalogs.
34 
35 When a transformer is called, it is handed a `~lsst.afw.table.SourceCatalog`
36 containing the measurements to be transformed, a `~lsst.afw.table.BaseCatalog`
37 in which to store results, and information about the WCS and calibration of
38 the data. It may be safely assumed that both are contiguous in memory, thus a
39 ``ColumnView`` may be used for efficient processing. If the transformation is
40 not possible, it should be aborted by throwing an exception; if this happens,
41 the caller should assume that the contents of the output catalog are
42 inconsistent.
43 
44 Transformations can be defined in Python or in C++. Python code should inherit
45 from `MeasurementTransform`, following its interface.
46 """
47 
48 from lsst.afw.table import CoordKey
49 from lsst.pex.exceptions import LengthError
50 from . import CentroidResultKey
51 
52 __all__ = ("MeasurementTransform", "NullTransform", "PassThroughTransform", "SimpleCentroidTransform")
53 
54 
56  """Base class for measurement transformations.
57 
58  Parameters
59  ----------
60  config : subclass of `BasePluginConfig`
61  The configuration of the measurement plugin whose outputs are being
62  transformed.
63  name : `str`
64  The name of the measurement plugin whose outputs are being
65  transformed.
66  mapper : `lsst.afw.table.SchemaMapper`
67  Mapping between the input (pre-transformation) and output
68  (transformed) catalogs.
69 
70  Notes
71  -----
72  Create transformations by deriving from this class, implementing
73  `__call__()` and (optionally) augmenting `__init__()`.
74  """
75 
76  def __init__(self, config, name, mapper):
77  self.namename = name
78  self.configconfig = config
79 
80  def __call__(self, inputCatalog, outputCatalog, wcs, photoCalib):
81  raise NotImplementedError()
82 
83  @staticmethod
84  def _checkCatalogSize(cat1, cat2):
85  if len(cat1) != len(cat2):
86  raise LengthError("Catalog size mismatch")
87 
88 
90  """Null transform which transfers no data from input to output.
91 
92  This is intended as the default for measurements for which no other
93  transformation is specified.
94 
95  Parameters
96  ----------
97  config : subclass of `BasePluginConfig`
98  The configuration of the measurement plugin whose outputs are being
99  transformed.
100  name : `str`
101  The name of the measurement plugin whose outputs are being
102  transformed.
103  mapper : `lsst.afw.table.SchemaMapper`
104  Mapping between the input (pre-transformation) and output
105  (transformed) catalogs.
106  """
107 
108  def __call__(self, inputCatalog, outputCatalog, wcs, photoCalib):
109  self._checkCatalogSize_checkCatalogSize(inputCatalog, outputCatalog)
110 
111 
113  """Copy fields from input to output without transformation.
114 
115  Parameters
116  ----------
117  config : subclass of `BasePluginConfig`
118  The configuration of the measurement plugin whose outputs are being
119  transformed.
120  name : `str`
121  The name of the measurement plugin whose outputs are being
122  transformed.
123  mapper : `lsst.afw.table.SchemaMapper`
124  Mapping between the input (pre-transformation) and output
125  (transformed) catalogs.
126  """
127 
128  def __init__(self, config, name, mapper):
129  MeasurementTransform.__init__(self, config, name, mapper)
130  for key, field in mapper.getInputSchema().extract(name + "*").values():
131  mapper.addMapping(key)
132 
133  def __call__(self, inputCatalog, outputCatalog, wcs, photoCalib):
134  self._checkCatalogSize_checkCatalogSize(inputCatalog, outputCatalog)
135 
136 
138  """Transform pixel centroid, without uncertainty, to celestial coordinates.
139 
140  Parameters
141  ----------
142  config : subclass of `BasePluginConfig`
143  The configuration of the measurement plugin whose outputs are being
144  transformed.
145  name : `str`
146  The name of the measurement plugin whose outputs are being
147  transformed.
148  mapper : `lsst.afw.table.SchemaMapper`
149  Mapping between the input (pre-transformation) and output
150  (transformed) catalogs.
151  """
152 
153  def __init__(self, config, name, mapper):
154  MeasurementTransform.__init__(self, config, name, mapper)
155  self.coordKeycoordKey = CoordKey.addFields(mapper.editOutputSchema(), name, "Position from " + name)
156 
157  def __call__(self, inputCatalog, outputCatalog, wcs, photoCalib):
158  self._checkCatalogSize_checkCatalogSize(inputCatalog, outputCatalog)
159  centroidResultKey = CentroidResultKey(inputCatalog.schema[self.namename])
160  for inSrc, outSrc in zip(inputCatalog, outputCatalog):
161  self.coordKeycoordKey.set(outSrc, wcs.pixelToSky(centroidResultKey.get(inSrc).getCentroid()))
def __init__(self, config, name, mapper)
Definition: transforms.py:76
def __call__(self, inputCatalog, outputCatalog, wcs, photoCalib)
Definition: transforms.py:80
def __call__(self, inputCatalog, outputCatalog, wcs, photoCalib)
Definition: transforms.py:108
def __init__(self, config, name, mapper)
Definition: transforms.py:128
def __call__(self, inputCatalog, outputCatalog, wcs, photoCalib)
Definition: transforms.py:133
def __call__(self, inputCatalog, outputCatalog, wcs, photoCalib)
Definition: transforms.py:157
def __init__(self, config, name, mapper)
Definition: transforms.py:153
Reports attempts to exceed implementation-defined length limits for some classes.
Definition: Runtime.h:76
daf::base::PropertySet * set
Definition: fits.cc:912