397 def __call__(self, kernelCellSet, log):
398 d1, d2, d3 = self.psfMatchConfig.alardDegGauss
400 for d1i
in range(1, d1+1):
401 for d2i
in range(1, d2+1):
402 for d3i
in range(1, d3+1):
403 dList = [d1i, d2i, d3i]
407 visitor = diffimLib.BuildSingleKernelVisitorF(kList, pexConfig.makePolicy(bicConfig))
408 visitor.setSkipBuilt(
False)
409 kernelCellSet.visitCandidates(visitor, bicConfig.nStarPerCell)
411 for cell
in kernelCellSet.getCellList():
412 for cand
in cell.begin(
False):
413 cand = diffimLib.cast_KernelCandidateF(cand)
414 if cand.getStatus() != afwMath.SpatialCellCandidate.GOOD:
416 diffIm = cand.getDifferenceImage(diffimLib.KernelCandidateF.RECENT)
417 bbox = cand.getKernel(diffimLib.KernelCandidateF.RECENT).shrinkBBox(diffIm.getBBox(afwImage.LOCAL))
418 diffIm = type(diffIm)(diffIm, bbox,
True)
419 chi2 = diffIm.getImage().getArray()**2 / diffIm.getVariance().getArray()
420 n = chi2.shape[0] * chi2.shape[1]
421 bic = np.sum(chi2) + k * np.log(n)
422 if not bicArray.has_key(cand.getId()):
423 bicArray[cand.getId()] = {}
424 bicArray[cand.getId()][(d1i, d2i, d3i)] = bic
427 candIds = bicArray.keys()
428 for candId
in candIds:
429 cconfig, cvals = bicArray[candId].keys(), bicArray[candId].values()
430 idx = np.argsort(cvals)
431 bestConfig = cconfig[idx[0]]
432 bestConfigs.append(bestConfig)
434 counter = Counter(bestConfigs).most_common(3)
435 log.info(
"B.I.C. prefers basis complexity %s %d times; %s %d times; %s %d times" % (counter[0][0], counter[0][1],
436 counter[1][0], counter[1][1],
437 counter[2][0], counter[2][1]))
438 return counter[0][0], counter[1][0], counter[2][0]
