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
|
ConfigClass = ForcedPhotCcdConfig RunnerClass = pipeBase.ButlerInitializedTaskRunner _DefaultName = "forcedPhotCcd" dataPrefix = "" def __init__(self, butler=None, refSchema=None, initInputs=None, **kwds): super().__init__(**kwds) if initInputs is not None: refSchema = initInputs['inputSchema'].schema self.makeSubtask("references", butler=butler, schema=refSchema) if refSchema is None: refSchema = self.references.schema self.makeSubtask("measurement", refSchema=refSchema) # It is necessary to get the schema internal to the forced measurement task until such a time # that the schema is not owned by the measurement task, but is passed in by an external caller if self.config.doApCorr: self.makeSubtask("applyApCorr", schema=self.measurement.schema) self.makeSubtask('catalogCalculation', schema=self.measurement.schema) self.outputSchema = lsst.afw.table.SourceCatalog(self.measurement.schema) def runQuantum(self, butlerQC, inputRefs, outputRefs): inputs = butlerQC.get(inputRefs) tract = butlerQC.quantum.dataId['tract'] skyMap = inputs.pop("skyMap") inputs['refWcs'] = skyMap[tract].getWcs() inputs['refCat'] = self.mergeAndFilterReferences(inputs['exposure'], inputs['refCat'], inputs['refWcs']) inputs['measCat'], inputs['exposureId'] = self.generateMeasCat(inputRefs.exposure.dataId, inputs['exposure'], inputs['refCat'], inputs['refWcs'], "visit_detector") self.attachFootprints(inputs['measCat'], inputs['refCat'], inputs['exposure'], inputs['refWcs']) # TODO: apply external calibrations (DM-17062) outputs = self.run(**inputs) butlerQC.put(outputs, outputRefs) def mergeAndFilterReferences(self, exposure, refCats, refWcs):
# Step 1: Determine bounds of the exposure photometry will # be performed on. expWcs = exposure.getWcs() expRegion = exposure.getBBox(lsst.afw.image.PARENT) expBBox = lsst.geom.Box2D(expRegion) expBoxCorners = expBBox.getCorners() expSkyCorners = [expWcs.pixelToSky(corner).getVector() for corner in expBoxCorners] expPolygon = lsst.sphgeom.ConvexPolygon(expSkyCorners) # Step 2: Filter out reference catalog sources that are # not contained within the exposure boundaries, or whose # parents are not within the exposure boundaries. Note # that within a single input refCat, the parents always # appear before the children. mergedRefCat = lsst.afw.table.SourceCatalog(refCats[0].table) for refCat in refCats: containedIds = {0} # zero as a parent ID means "this is a parent" for record in refCat: if expPolygon.contains(record.getCoord().getVector()) and record.getParent() in containedIds: record.setFootprint(record.getFootprint().transform(refWcs, expWcs, expRegion)) mergedRefCat.append(record) containedIds.add(record.getId()) mergedRefCat.sort(lsst.afw.table.SourceTable.getParentKey()) return mergedRefCat def generateMeasCat(self, exposureDataId, exposure, refCat, refWcs, idPackerName):
expId, expBits = exposureDataId.pack(idPackerName, returnMaxBits=True) idFactory = lsst.afw.table.IdFactory.makeSource(expId, 64 - expBits) measCat = self.measurement.generateMeasCat(exposure, refCat, refWcs, idFactory=idFactory) return measCat, expId def runDataRef(self, dataRef, psfCache=None):
Definition at line 450 of file forcedPhotCcd.py.