LSST Applications g04a91732dc+146a938ab0,g07dc498a13+80b84b0d75,g0fba68d861+4c4f3dcb5c,g1409bbee79+80b84b0d75,g1a7e361dbc+80b84b0d75,g1fd858c14a+f6e422e056,g20f46db602+333b6c0f32,g35bb328faa+fcb1d3bbc8,g42c1b31a95+a1301e4c20,g4d2262a081+f1facf12e5,g4d39ba7253+9b833be27e,g4e0f332c67+5d362be553,g53246c7159+fcb1d3bbc8,g60b5630c4e+9b833be27e,g78460c75b0+2f9a1b4bcd,g786e29fd12+cf7ec2a62a,g7b71ed6315+fcb1d3bbc8,g8852436030+790117df0f,g89139ef638+80b84b0d75,g8d6b6b353c+9b833be27e,g9125e01d80+fcb1d3bbc8,g989de1cb63+80b84b0d75,g9f33ca652e+9c6b68d7f3,ga9baa6287d+9b833be27e,gaaedd4e678+80b84b0d75,gabe3b4be73+1e0a283bba,gb1101e3267+9f3571abad,gb58c049af0+f03b321e39,gb90eeb9370+691e4ab549,gc741bbaa4f+2bcd3860df,gcf25f946ba+790117df0f,gd315a588df+5b65d88fe4,gd6cbbdb0b4+c8606af20c,gd9a9a58781+fcb1d3bbc8,gde0f65d7ad+ee6a3faa19,ge278dab8ac+932305ba37,ge82c20c137+76d20ab76d,gee8db133a9+2a6ae0040b,w.2025.10
LSST Data Management Base Package
|
Classes | |
class | MergeMeasurementsConnections |
Variables | |
schema : `lsst.afw.table.Schema`, optional | |
initInputs : `dict`, optional | |
catalogs : `lsst.afw.table.SourceCatalog` | |
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.
lsst.pipe.tasks.mergeMeasurements.initInputs : `dict`, optional |
Definition at line 131 of file mergeMeasurements.py.
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.