Hide keyboard shortcuts

Hot-keys on this page

r m x p   toggle line displays

j k   next/prev highlighted chunk

0   (zero) top of page

1   (one) first highlighted chunk

1

2

3

4

5

6

7

8

9

10

11

12

13

14

15

16

17

18

19

20

21

22

23

24

25

26

27

28

29

30

31

32

33

34

35

36

37

38

39

40

41

42

43

44

45

46

47

48

49

50

51

52

53

54

55

56

57

58

59

60

61

62

63

64

65

66

67

68

69

70

71

72

73

74

75

76

77

78

79

80

81

82

83

84

85

86

87

88

89

90

91

92

93

94

95

96

97

98

99

100

101

102

103

104

105

106

107

108

109

110

111

112

113

114

115

116

117

118

119

120

121

122

123

124

125

126

127

128

129

130

131

132

133

134

135

136

137

138

139

140

141

142

143

144

145

146

147

148

149

150

151

152

153

154

155

156

157

158

159

160

161

162

163

164

165

166

167

168

169

170

171

172

173

174

175

176

177

178

179

180

181

182

183

184

185

186

187

188

189

190

191

192

193

194

195

196

197

198

199

200

201

202

203

204

205

206

207

208

209

210

211

212

213

214

215

216

217

218

219

220

221

222

223

224

225

226

227

228

229

230

231

232

233

234

import unittest 

 

import numpy as num 

 

import lsst.utils.tests 

import lsst.afw.image as afwImage 

import lsst.afw.math as afwMath 

import lsst.ip.diffim as ipDiffim 

import lsst.log.utils as logUtils 

import lsst.pex.config as pexConfig 

import lsst.pex.exceptions 

 

# Increase the number for more verbose messages 

logUtils.traceSetAt("ip.diffim", 0) 

 

 

class DiffimTestCases(unittest.TestCase): 

 

def setUp(self): 

self.configAL = ipDiffim.ImagePsfMatchTask.ConfigClass() 

self.configAL.kernel.name = "AL" 

self.subconfigAL = self.configAL.kernel.active 

 

self.configDF = ipDiffim.ImagePsfMatchTask.ConfigClass() 

self.configDF.kernel.name = "DF" 

self.subconfigDF = self.configDF.kernel.active 

 

self.policyAL = pexConfig.makePolicy(self.subconfigAL) 

self.policyDF = pexConfig.makePolicy(self.subconfigDF) 

 

self.kSize = self.policyAL.getInt("kernelSize") 

 

def tearDown(self): 

del self.configAL 

del self.policyAL 

del self.configDF 

del self.policyDF 

 

def deltaFunctionTest(self, ks): 

# right shape 

nk = 0 

for rowi in range(self.kSize): 

for colj in range(self.kSize): 

kernel = ks[nk] 

kimage = afwImage.ImageD(kernel.getDimensions()) 

ksum = kernel.computeImage(kimage, False) 

self.assertEqual(ksum, 1.) 

 

for rowk in range(self.kSize): 

for coll in range(self.kSize): 

if (rowi == rowk) and (colj == coll): 

self.assertEqual(kimage[coll, rowk, afwImage.LOCAL], 1.) 

else: 

self.assertEqual(kimage[coll, rowk, afwImage.LOCAL], 0.) 

nk += 1 

 

def testMakeDeltaFunction(self): 

ks = ipDiffim.makeDeltaFunctionBasisList(self.kSize, self.kSize) 

 

# right size 

self.assertEqual(len(ks), self.kSize * self.kSize) 

 

# right shape 

self.deltaFunctionTest(ks) 

 

def alardLuptonTest(self, ks): 

kim = afwImage.ImageD(ks[0].getDimensions()) 

nBasis = len(ks) 

 

# the first one sums to 1; the rest sum to 0 

ks[0].computeImage(kim, False) 

self.assertAlmostEqual(num.sum(num.ravel(kim.getArray())), 1.0) 

 

for k in range(1, nBasis): 

ks[k].computeImage(kim, False) 

self.assertAlmostEqual(num.sum(num.ravel(kim.getArray())), 0.0) 

 

# the images dotted with themselves is 1, except for the first 

for k in range(1, nBasis): 

ks[k].computeImage(kim, False) 

arr = kim.getArray() 

self.assertAlmostEqual(num.sum(arr*arr), 1.0) 

 

def testMakeAlardLupton(self): 

nGauss = self.policyAL.get("alardNGauss") 

sigGauss = self.policyAL.getDoubleArray("alardSigGauss") 

degGauss = self.policyAL.getIntArray("alardDegGauss") 

self.assertEqual(len(sigGauss), nGauss) 

self.assertEqual(len(degGauss), nGauss) 

self.assertEqual(self.kSize % 2, 1) # odd sized 

kHalfWidth = self.kSize // 2 

 

ks = ipDiffim.makeAlardLuptonBasisList(kHalfWidth, nGauss, sigGauss, degGauss) 

 

# right size 

nTot = 0 

for deg in degGauss: 

nTot += (deg + 1) * (deg + 2) / 2 

self.assertEqual(len(ks), nTot) 

 

# right orthogonality 

self.alardLuptonTest(ks) 

 

def testGenerateAlardLupton(self): 

# defaults 

ks = ipDiffim.generateAlardLuptonBasisList(self.subconfigAL) 

self.alardLuptonTest(ks) 

 

# send FWHM 

ks = ipDiffim.generateAlardLuptonBasisList(self.subconfigAL, targetFwhmPix=3.0, referenceFwhmPix=4.0) 

self.alardLuptonTest(ks) 

 

def testMakeKernelBasisList(self): 

ks = ipDiffim.makeKernelBasisList(self.subconfigAL) 

self.alardLuptonTest(ks) 

 

ks = ipDiffim.makeKernelBasisList(self.subconfigDF) 

self.deltaFunctionTest(ks) 

 

def testRenormalize(self): 

# inputs 

gauss1 = afwMath.GaussianFunction2D(2, 2) 

gauss2 = afwMath.GaussianFunction2D(3, 4) 

gauss3 = afwMath.GaussianFunction2D(0.2, 0.25) 

gaussKernel1 = afwMath.AnalyticKernel(self.kSize, self.kSize, gauss1) 

gaussKernel2 = afwMath.AnalyticKernel(self.kSize, self.kSize, gauss2) 

gaussKernel3 = afwMath.AnalyticKernel(self.kSize, self.kSize, gauss3) 

kimage1 = afwImage.ImageD(gaussKernel1.getDimensions()) 

ksum1 = gaussKernel1.computeImage(kimage1, False) 

kimage2 = afwImage.ImageD(gaussKernel2.getDimensions()) 

ksum2 = gaussKernel2.computeImage(kimage2, False) 

kimage3 = afwImage.ImageD(gaussKernel3.getDimensions()) 

ksum3 = gaussKernel3.computeImage(kimage3, False) 

self.assertNotEqual(ksum1, 1.) 

self.assertNotEqual(ksum2, 1.) 

self.assertNotEqual(ksum3, 1.) 

# no constraints on first kernels norm 

self.assertNotEqual(num.sum(num.ravel(kimage2.getArray())**2), 1.) 

self.assertNotEqual(num.sum(num.ravel(kimage3.getArray())**2), 1.) 

basisListIn = [] 

basisListIn.append(gaussKernel1) 

basisListIn.append(gaussKernel2) 

basisListIn.append(gaussKernel3) 

 

# outputs 

basisListOut = ipDiffim.renormalizeKernelList(basisListIn) 

gaussKernel1 = basisListOut[0] 

gaussKernel2 = basisListOut[1] 

gaussKernel3 = basisListOut[2] 

ksum1 = gaussKernel1.computeImage(kimage1, False) 

ksum2 = gaussKernel2.computeImage(kimage2, False) 

ksum3 = gaussKernel3.computeImage(kimage3, False) 

self.assertAlmostEqual(ksum1, 1.) 

self.assertAlmostEqual(ksum2, 0.) 

self.assertAlmostEqual(ksum3, 0.) 

# no constraints on first kernels norm 

self.assertAlmostEqual(num.sum(num.ravel(kimage2.getArray())**2), 1.) 

self.assertAlmostEqual(num.sum(num.ravel(kimage3.getArray())**2), 1.) 

 

def testCentralRegularization(self): 

# stencil of 1 not allowed 

self.policyDF.set("regularizationType", "centralDifference") 

with self.assertRaises(lsst.pex.exceptions.Exception): 

self.policyDF.set("centralRegularizationStencil", 1) 

ipDiffim.makeRegularizationMatrix(self.policyDF) 

 

# stencil of 5 allowed 

self.policyDF.set("centralRegularizationStencil", 5) 

try: 

ipDiffim.makeRegularizationMatrix(self.policyDF) 

except lsst.pex.exceptions.Exception as e: 

self.fail("Should not raise %s: stencil of 5 is allowed."%e) 

 

# stencil of 9 allowed 

self.policyDF.set("centralRegularizationStencil", 9) 

try: 

ipDiffim.makeRegularizationMatrix(self.policyDF) 

except lsst.pex.exceptions.Exception as e: 

self.fail("Should not raise %s: stencil of 9 is allowed"%e) 

 

# border penalty < 0 

self.policyDF.set("regularizationBorderPenalty", -1.0) 

with self.assertRaises(lsst.pex.exceptions.Exception): 

ipDiffim.makeRegularizationMatrix(self.policyDF) 

 

# border penalty > 0 

self.policyDF.set("regularizationBorderPenalty", 0.0) 

try: 

ipDiffim.makeRegularizationMatrix(self.policyDF) 

except lsst.pex.exceptions.Exception as e: 

self.fail("Should not raise %s: Border penalty > 0"%e) 

 

def testForwardRegularization(self): 

self.policyDF.set("regularizationType", "forwardDifference") 

 

# order 1..3 allowed 

self.policyDF.set("forwardRegularizationOrders", 0) 

with self.assertRaises(lsst.pex.exceptions.Exception): 

ipDiffim.makeRegularizationMatrix(self.policyDF) 

 

self.policyDF.set("forwardRegularizationOrders", 1) 

try: 

ipDiffim.makeRegularizationMatrix(self.policyDF) 

except lsst.pex.exceptions.Exception as e: 

self.fail("Should not raise %s: order 1 allowed"%e) 

 

self.policyDF.set("forwardRegularizationOrders", 4) 

with self.assertRaises(lsst.pex.exceptions.Exception): 

ipDiffim.makeRegularizationMatrix(self.policyDF) 

 

self.policyDF.set("forwardRegularizationOrders", 1) 

self.policyDF.add("forwardRegularizationOrders", 2) 

try: 

ipDiffim.makeRegularizationMatrix(self.policyDF) 

except lsst.pex.exceptions.Exception as e: 

self.fail("Should not raise %s: order 1,2 allowed"%e) 

 

def testBadRegularization(self): 

with self.assertRaises(lsst.pex.exceptions.Exception): 

self.policyDF.set("regularizationType", "foo") 

ipDiffim.makeRegularizationMatrix(self.policyDF) 

 

 

class TestMemory(lsst.utils.tests.MemoryTestCase): 

pass 

 

 

def setup_module(module): 

lsst.utils.tests.init() 

 

 

232 ↛ 233line 232 didn't jump to line 233, because the condition on line 232 was never trueif __name__ == "__main__": 

lsst.utils.tests.init() 

unittest.main()