lsst.ip.diffim.diffimTools.NbasisEvaluator Class Reference

Inheritance diagram for lsst.ip.diffim.diffimTools.NbasisEvaluator:

Public Member Functions
def	__init__ (self, psfMatchConfig, psfFwhmPixTc, psfFwhmPixTnc)
def	__call__ (self, kernelCellSet, log)

Public Attributes
	psfMatchConfig
	psfFwhmPixTc
	psfFwhmPixTnc

Detailed Description

A functor to evaluate the Bayesian Information Criterion for the number of basis sets
going into the kernel fitting

Definition at line 535 of file diffimTools.py.

Constructor & Destructor Documentation

__init__()

def lsst.ip.diffim.diffimTools.NbasisEvaluator.__init__	(		self,
			psfMatchConfig,
			psfFwhmPixTc,
			psfFwhmPixTnc
	)

Definition at line 539 of file diffimTools.py.

    def __init__(self, psfMatchConfig, psfFwhmPixTc, psfFwhmPixTnc):
        self.psfMatchConfig = psfMatchConfig
        self.psfFwhmPixTc = psfFwhmPixTc
        self.psfFwhmPixTnc = psfFwhmPixTnc
        if not self.psfMatchConfig.kernelBasisSet == "alard-lupton":
            raise RuntimeError("BIC only implemnted for AL (alard lupton) basis")
 

Member Function Documentation

__call__()

def lsst.ip.diffim.diffimTools.NbasisEvaluator.__call__	(		self,
			kernelCellSet,
			log
	)

Definition at line 546 of file diffimTools.py.

    def __call__(self, kernelCellSet, log):
        d1, d2, d3 = self.psfMatchConfig.alardDegGauss
        bicArray = {}
        for d1i in range(1, d1 + 1):
            for d2i in range(1, d2 + 1):
                for d3i in range(1, d3 + 1):
                    dList = [d1i, d2i, d3i]
                    bicConfig = type(self.psfMatchConfig)(self.psfMatchConfig, alardDegGauss=dList)
                    kList = makeKernelBasisList(bicConfig, self.psfFwhmPixTc, self.psfFwhmPixTnc)
                    k = len(kList)
                    visitor = diffimLib.BuildSingleKernelVisitorF(kList,
                                                                  pexConfig.makePropertySet(bicConfig))
                    visitor.setSkipBuilt(False)
                    kernelCellSet.visitCandidates(visitor, bicConfig.nStarPerCell)
 
                    for cell in kernelCellSet.getCellList():
                        for cand in cell.begin(False):  # False = include bad candidates
                            if cand.getStatus() != afwMath.SpatialCellCandidate.GOOD:
                                continue
                            diffIm = cand.getDifferenceImage(diffimLib.KernelCandidateF.RECENT)
                            bbox = cand.getKernel(diffimLib.KernelCandidateF.RECENT).shrinkBBox(
                                diffIm.getBBox(afwImage.LOCAL))
                            diffIm = type(diffIm)(diffIm, bbox, True)
                            chi2 = diffIm.getImage().getArray()**2/diffIm.getVariance().getArray()
                            n = chi2.shape[0]*chi2.shape[1]
                            bic = np.sum(chi2) + k*np.log(n)
                            if cand.getId() not in bicArray:
                                bicArray[cand.getId()] = {}
                            bicArray[cand.getId()][(d1i, d2i, d3i)] = bic
 
        bestConfigs = []
        for candId in bicArray:
            cconfig, cvals = list(bicArray[candId].keys()), list(bicArray[candId].values())
            idx = np.argsort(cvals)
            bestConfig = cconfig[idx[0]]
            bestConfigs.append(bestConfig)
 
        counter = Counter(bestConfigs).most_common(3)
        log.info("B.I.C. prefers basis complexity %s %d times; %s %d times; %s %d times",
                 counter[0][0], counter[0][1],
                 counter[1][0], counter[1][1],
                 counter[2][0], counter[2][1])
        return counter[0][0], counter[1][0], counter[2][0]

Member Data Documentation

psfFwhmPixTc

lsst.ip.diffim.diffimTools.NbasisEvaluator.psfFwhmPixTc

Definition at line 541 of file diffimTools.py.

psfFwhmPixTnc

lsst.ip.diffim.diffimTools.NbasisEvaluator.psfFwhmPixTnc

Definition at line 542 of file diffimTools.py.

psfMatchConfig

lsst.ip.diffim.diffimTools.NbasisEvaluator.psfMatchConfig

Definition at line 540 of file diffimTools.py.

---

The documentation for this class was generated from the following file:

• /j/snowflake/release/lsstsw/stack/lsst-scipipe-0.7.0/Linux64/ip_diffim/22.0.1-24-g1ad7a390+a9625a72a8/python/lsst/ip/diffim/diffimTools.py