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
classification.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 """Definition and registration of classification plugins.
23 """
24 
25 import numpy as np
26 
27 import lsst.pex.config
28 from .catalogCalculation import CatalogCalculationPluginConfig, CatalogCalculationPlugin
29 from .pluginRegistry import register
30 
31 __all__ = (
32  "CatalogCalculationClassificationConfig", "CatalogCalculationClassificationPlugin",
33 )
34 
35 
37  """Configuration for catalog classification plugin.
38  """
39 
40  fluxRatio = lsst.pex.config.Field(dtype=float, default=.925, optional=True,
41  doc="critical ratio of model to psf flux")
42  modelErrFactor = lsst.pex.config.Field(dtype=float, default=0.0, optional=True,
43  doc="correction factor for modelFlux error")
44  psfErrFactor = lsst.pex.config.Field(dtype=float, default=0.0, optional=True,
45  doc="correction factor for psfFlux error")
46 
47 
48 @register("base_ClassificationExtendedness")
50  """Plugin which calculates a binary measure of source extendedness.
51 
52  Extendedness is based on a simple cut of the ratio of the PSF flux to the
53  model flux.
54 
55  Notes
56  -----
57  Because the fluxes on which this algorithm is based on are slot
58  measurements, they can be provided by different algorithms, and the
59  `~CatalogCalculationClassificationConfig.fluxRatio` threshold used by this
60  algorithm should generally be set differently for different algorithms.
61  To do this, plot the difference between the PSF magnitude and the model
62  magnitude vs. the PSF magnitude, and look for where the cloud of galaxies
63  begins.
64  """
65 
66  ConfigClass = CatalogCalculationClassificationConfig
67 
68  @classmethod
70  return cls.DEFAULT_CATALOGCALCULATIONDEFAULT_CATALOGCALCULATION
71 
72  def __init__(self, config, name, schema, metadata):
73  CatalogCalculationPlugin.__init__(self, config, name, schema, metadata)
74  self.keyProbabilitykeyProbability = schema.addField(name + "_value", type="D",
75  doc="Set to 1 for extended sources, 0 for point sources.")
76  self.keyFlagkeyFlag = schema.addField(name + "_flag", type="Flag", doc="Set to 1 for any fatal failure.")
77 
78  def calculate(self, measRecord):
79  modelFlux = measRecord.getModelInstFlux()
80  psfFlux = measRecord.getPsfInstFlux()
81  modelFluxFlag = (measRecord.getModelFluxFlag()
82  if measRecord.table.getModelFluxSlot().isValid()
83  else False)
84  psfFluxFlag = (measRecord.getPsfFluxFlag()
85  if measRecord.table.getPsfFluxSlot().isValid()
86  else False)
87  flux1 = self.configconfig.fluxRatio*modelFlux
88  if self.configconfig.modelErrFactor != 0:
89  flux1 += self.configconfig.modelErrFactor*measRecord.getModelInstFluxErr()
90  flux2 = psfFlux
91  if not self.configconfig.psfErrFactor == 0:
92  flux2 += self.configconfig.psfErrFactor*measRecord.getPsfInstFluxErr()
93 
94  # A generic failure occurs when either FluxFlag is set to True
95  # A generic failure also occurs if either calculated flux value is NaN:
96  # this can occur if the Flux field itself is NaN,
97  # or the ErrFactor != 0 and the FluxErr is NaN
98  if np.isnan(flux1) or np.isnan(flux2) or modelFluxFlag or psfFluxFlag:
99  self.failfailfail(measRecord)
100  else:
101  measRecord.set(self.keyProbabilitykeyProbability, 0.0 if flux1 < flux2 else 1.0)
102 
103  def fail(self, measRecord, error=None):
104  measRecord.set(self.keyFlagkeyFlag, True)
def fail(self, measRecord, error=None)
Definition: pluginsBase.py:137
bool isValid
Definition: fits.cc:399