20 namespace afwMath = lsst::afw::math;
21 namespace afwImage = lsst::afw::image;
22 namespace pexExcept = lsst::pex::exceptions;
55 template<
typename PixelT>
59 afwMath::CandidateVisitor(),
64 template<
typename PixelT>
69 int ncomp = eigenImages.
size();
77 for (
int i = 0; i < ncomp; i++) {
88 template<
typename PixelT>
92 if (kCandidate == NULL) {
94 "Failed to cast SpatialCellCandidate to KernelCandidate");
96 LOGL_DEBUG(
"TRACE5.ip.diffim.SetPcaImageVisitor.processCandidate",
97 "Processing candidate %d", kCandidate->
getId());
106 _imagePca->addImage(kImage, 1.0);
112 template<
typename PixelT>
137 LOGL_DEBUG(
"TRACE5.ip.diffim.KernelPcaVisitor.subtractMean",
138 "Subtracting mean feature before Pca");
140 _mean = _imagePca->getMean();
143 end = imageList.
end(); ptr !=
end; ++ptr) {
162 template <
typename ImageT>
168 typename Super::ImageList::const_iterator
iter = eImageList.
begin(),
end = eImageList.
end();
176 afwMath::Statistics stats = afwMath::makeStatistics(*eImage, (afwMath::MIN | afwMath::MAX));
181 if (extreme != 0.0) {