def testLambdaEff(self):
f = afwImage.Filter("g")
g_r = 1.2
c = afwImage.Color(g_r)
self.assertEqual(c.getLambdaEff(f),
1000 * g_r) # XXX Not a real implementation!
def testComputeImage2(self):
"""Test the computation of the PSF's image at a point"""
color = afwImage.Color(1.0)
ccdXY = afwGeom.Point2D(0, 0)
kIm = self.psf.computeImage(ccdXY)
self.assertTrue(kIm.getWidth() == self.ksize)
self.assertAlmostEqual(afwMath.makeStatistics(kIm, afwMath.SUM).getValue(), 1.0)
def determinePsf(self,
exposure,
psfCandidateList,
metadata=None,
flagKey=None):
assert flagKey is not None
flag0 = [c.getSource().get(flagKey) for c in psfCandidateList]
xsrc = np.array([c.getSource().getX() for c in psfCandidateList])
ysrc = np.array([c.getSource().getY() for c in psfCandidateList])
# These are (x,y) values where the psfs will be compared
xvec = np.linspace(np.min(xsrc) + 16.32, np.max(xsrc) - 19.77, 5)
yvec = np.linspace(np.min(ysrc) + 20.45, np.max(ysrc) - 23.89, 5)
print '\nComp2PsfDeterminer: starting first psf run (%s)\n' % self.config.psfName1
determiner1 = self.determiner_class1(self.determiner_config1)
(psf1, psfCellSet1) = determiner1.determinePsf(exposure,
psfCandidateList,
metadata, flagKey)
flag1 = [c.getSource().get(flagKey) for c in psfCandidateList]
print '\nComp2PsfDeterminer: starting second psf run (%s)\n' % self.config.psfName2
for c in psfCandidateList:
c.getSource().set(flagKey, False)
determiner2 = self.determiner_class2(self.determiner_config2)
(psf2, psfCellSet2) = determiner2.determinePsf(exposure,
psfCandidateList,
metadata, flagKey)
flag2 = [c.getSource().get(flagKey) for c in psfCandidateList]
print '\nComp2PsfDeterminer: table of comparison values follows'
for x in xvec:
for y in yvec:
print '(x,y) =', (x, y)
self.compare_images(
'doComputeImage',
psf1.doComputeImage(afwGeom.Point2D(x, y),
afwImage.Color()),
psf2.doComputeImage(afwGeom.Point2D(x, y),
afwImage.Color()))
self.compare_images(
'doComputeKernelImage',
psf1.doComputeKernelImage(afwGeom.Point2D(x, y),
afwImage.Color()),
psf2.doComputeKernelImage(afwGeom.Point2D(x, y),
afwImage.Color()))
bad1 = [i for (i, x) in enumerate(flag1) if not x]
bad2 = [i for (i, x) in enumerate(flag2) if not x]
print '\nComp2PsfDeterminer: the following candidates were marked bad by %s: %s' % (
self.config.psfName1, bad1)
print 'Comp2PsfDeterminer: the following candidates were marked bad by %s: %s\n' % (
self.config.psfName2, bad2)
if bad1 != bad2:
print ' !!! WARNING the two sets of bad candidates differ !!!\n'
return (psf2, psfCellSet2)
def testIsIndeterminate(self):
"""Test that a default-constructed Color tests True, but ones with a g-r value test False"""
self.assertTrue(afwImage.Color().isIndeterminate())
self.assertFalse(afwImage.Color(1.2).isIndeterminate())
def testCtor(self):
afwImage.Color()
afwImage.Color(1.2)
def testCtor(self):
c = afwImage.Color()
c = afwImage.Color(1.2)