lsst.pipe.tasks.mergeMeasurements Namespace Reference

Classes
class	MergeMeasurementsConnections

Variables
	schema : `lsst.afw.table.Schema`, optional
	initInputs : `dict`, optional
	catalogs : `lsst.afw.table.SourceCatalog`

Variable Documentation

catalogs

lsst.pipe.tasks.mergeMeasurements.catalogs : `lsst.afw.table.SourceCatalog`

_DefaultName = "mergeCoaddMeasurements"
ConfigClass = MergeMeasurementsConfig

inputDataset = "meas"
outputDataset = "ref"

def __init__(self, schema=None, initInputs=None, **kwargs):
    super().__init__(**kwargs)

    if initInputs is not None:
        schema = initInputs['inputSchema'].schema

    if schema is None:
        raise ValueError("No input schema or initInputs['inputSchema'] provided.")

    inputSchema = schema

    self.schemaMapper = afwTable.SchemaMapper(inputSchema, True)
    self.schemaMapper.addMinimalSchema(inputSchema, True)
    self.instFluxKey = inputSchema.find(self.config.snName + "_instFlux").getKey()
    self.instFluxErrKey = inputSchema.find(self.config.snName + "_instFluxErr").getKey()
    self.fluxFlagKey = inputSchema.find(self.config.snName + "_flag").getKey()

    self.flagKeys = {}
    for band in self.config.priorityList:
        outputKey = self.schemaMapper.editOutputSchema().addField(
            "merge_measurement_%s" % band,
            type="Flag",
            doc="Flag field set if the measurements here are from the %s filter" % band
        )
        peakKey = inputSchema.find("merge_peak_%s" % band).key
        footprintKey = inputSchema.find("merge_footprint_%s" % band).key
        self.flagKeys[band] = pipeBase.Struct(peak=peakKey, footprint=footprintKey, output=outputKey)
    self.schema = self.schemaMapper.getOutputSchema()

    self.pseudoFilterKeys = []
    for filt in self.config.pseudoFilterList:
        try:
            self.pseudoFilterKeys.append(self.schema.find("merge_peak_%s" % filt).getKey())
        except Exception as e:
            self.log.warning("merge_peak is not set for pseudo-filter %s: %s", filt, e)

    self.badFlags = {}
    for flag in self.config.flags:
        try:
            self.badFlags[flag] = self.schema.find(flag).getKey()
        except KeyError as exc:
            self.log.warning("Can't find flag %s in schema: %s", flag, exc)
    self.outputSchema = afwTable.SourceCatalog(self.schema)

def runQuantum(self, butlerQC, inputRefs, outputRefs):
    inputs = butlerQC.get(inputRefs)
    dataIds = (ref.dataId for ref in inputRefs.catalogs)
    catalogDict = {dataId['band']: cat for dataId, cat in zip(dataIds, inputs['catalogs'])}
    inputs['catalogs'] = catalogDict
    outputs = self.run(**inputs)
    butlerQC.put(outputs, outputRefs)

def run(self, catalogs):

Definition at line 201 of file mergeMeasurements.py.

initInputs

lsst.pipe.tasks.mergeMeasurements.initInputs : `dict`, optional

Definition at line 131 of file mergeMeasurements.py.

schema

lsst.pipe.tasks.mergeMeasurements.schema : `lsst.afw.table.Schema`, optional

pseudoFilterList = pexConfig.ListField(
    dtype=str,
    default=["sky"],
    doc="Names of filters which may have no associated detection\n"
    "(N.b. should include MergeDetectionsConfig.skyFilterName)"
)
snName = pexConfig.Field(
    dtype=str,
    default="base_PsfFlux",
    doc="Name of flux measurement for calculating the S/N when choosing the reference band."
)
minSN = pexConfig.Field(
    dtype=float,
    default=10.,
    doc="If the S/N from the priority band is below this value (and the S/N "
    "is larger than minSNDiff compared to the priority band), use the band with "
    "the largest S/N as the reference band."
)
minSNDiff = pexConfig.Field(
    dtype=float,
    default=3.,
    doc="If the difference in S/N between another band and the priority band is larger "
    "than this value (and the S/N in the priority band is less than minSN) "
    "use the band with the largest S/N as the reference band"
)
flags = pexConfig.ListField(
    dtype=str,
    doc="Require that these flags, if available, are not set",
    default=["base_PixelFlags_flag_interpolatedCenter", "base_PsfFlux_flag",
             "ext_photometryKron_KronFlux_flag", "modelfit_CModel_flag", ]
)
priorityList = pexConfig.ListField(
    dtype=str,
    default=[],
    doc="Priority-ordered list of filter bands for the merge."
)
coaddName = pexConfig.Field(
    dtype=str,
    default="deep",
    doc="Name of coadd"
)

def validate(self):
    super().validate()
    if len(self.priorityList) == 0:
        raise RuntimeError("No priority list provided")


class MergeMeasurementsTask(pipeBase.PipelineTask):

Definition at line 129 of file mergeMeasurements.py.