63 def run(self, catalog, filterName=None, epoch=None):
64 """Load reference objects and match to them.
65
66 Parameters
67 ----------
68 catalog : `lsst.afw.table.SourceCatalog`
69 Catalog to match.
70 filterName : `str`
71 Name of filter loading fluxes.
72 epoch : `astropy.time.Time` or `None`
73 Epoch to which to correct proper motion and parallax, or `None` to
74 not apply such corrections.
75
76 Returns
77 -------
78 result : `lsst.pipe.base.Struct`
79 Result struct with components:
80
81 ``matches``
82 Matched sources with associated reference
83 (`lsst.afw.table.SourceMatchVector`).
84 ``matchMeta``
85 Match metadata (`lsst.meas.astrom.MatchMetadata`).
86 """
87 if self.refObjLoader is None:
88 raise RuntimeError("Running matcher task with no refObjLoader set in __ini__ or setRefObjLoader")
89 circle = self.calculateCircle(catalog)
90 matchMeta = self.refObjLoader.getMetadataCircle(circle.center, circle.radius, filterName, epoch=epoch)
91 emptyResult = Struct(matches=[], matchMeta=matchMeta)
92 sourceSelection = self.sourceSelection.run(catalog)
93 if len(sourceSelection.sourceCat) == 0:
94 self.log.warning("No objects selected from %d objects in source catalog", len(catalog))
95 return emptyResult
96 refData = self.refObjLoader.loadSkyCircle(circle.center, circle.radius, filterName, epoch=epoch)
97 refCat = refData.refCat
98 refSelection = self.referenceSelection.run(refCat)
99 if len(refSelection.sourceCat) == 0:
100 self.log.warning("No objects selected from %d objects in reference catalog", len(refCat))
101 return emptyResult
103 self.config.matchRadius*arcseconds)
104 self.log.info("Matched %d from %d/%d input and %d/%d reference sources",
105 len(matches), len(sourceSelection.sourceCat), len(catalog),
106 len(refSelection.sourceCat), len(refCat))
107 return Struct(matches=matches, matchMeta=matchMeta, refCat=refCat, sourceSelection=sourceSelection,
108 refSelection=refSelection)
109
std::vector< Match< typename Cat1::Record, typename Cat2::Record > > matchRaDec(Cat1 const &cat1, Cat2 const &cat2, lsst::geom::Angle radius, MatchControl const &mc=MatchControl())
Compute all tuples (s1,s2,d) where s1 belings to cat1, s2 belongs to cat2 and d, the distance between...