Coverage for tests/test_PixelFlags.py: 14%

117 statements  

« prev     ^ index     » next       coverage.py v7.4.4, created at 2024-03-20 04:26 -0700

1# This file is part of meas_base. 

2# 

3# Developed for the LSST Data Management System. 

4# This product includes software developed by the LSST Project 

5# (https://www.lsst.org). 

6# See the COPYRIGHT file at the top-level directory of this distribution 

7# for details of code ownership. 

8# 

9# This program is free software: you can redistribute it and/or modify 

10# it under the terms of the GNU General Public License as published by 

11# the Free Software Foundation, either version 3 of the License, or 

12# (at your option) any later version. 

13# 

14# This program is distributed in the hope that it will be useful, 

15# but WITHOUT ANY WARRANTY; without even the implied warranty of 

16# MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the 

17# GNU General Public License for more details. 

18# 

19# You should have received a copy of the GNU General Public License 

20# along with this program. If not, see <https://www.gnu.org/licenses/>. 

21 

22import unittest 

23 

24import numpy as np 

25 

26import lsst.geom 

27import lsst.utils.tests 

28import lsst.meas.base.tests 

29 

30 

31class PixelFlagsTestCase(lsst.meas.base.tests.AlgorithmTestCase, lsst.utils.tests.TestCase): 

32 

33 def setUp(self): 

34 self.center = lsst.geom.Point2D(50.1, 49.8) 

35 self.bbox = lsst.geom.Box2I(lsst.geom.Point2I(-20, -30), 

36 lsst.geom.Extent2I(140, 160)) 

37 self.dataset = lsst.meas.base.tests.TestDataset(self.bbox) 

38 

39 def testNoFlags(self): 

40 self.dataset.addSource(100000.0, self.center) 

41 task = self.makeSingleFrameMeasurementTask("base_PixelFlags") 

42 exposure, catalog = self.dataset.realize(10.0, task.schema, randomSeed=0) 

43 task.run(catalog, exposure) 

44 record = catalog[0] 

45 self.assertFalse(record.get("base_PixelFlags_flag")) 

46 self.assertFalse(record.get("base_PixelFlags_flag_edge")) 

47 self.assertFalse(record.get("base_PixelFlags_flag_edgeCenter")) 

48 self.assertFalse(record.get("base_PixelFlags_flag_edgeCenterAll")) 

49 self.assertFalse(record.get("base_PixelFlags_flag_interpolated")) 

50 self.assertFalse(record.get("base_PixelFlags_flag_interpolatedCenter")) 

51 self.assertFalse(record.get("base_PixelFlags_flag_interpolatedCenterAll")) 

52 self.assertFalse(record.get("base_PixelFlags_flag_saturated")) 

53 self.assertFalse(record.get("base_PixelFlags_flag_saturatedCenter")) 

54 self.assertFalse(record.get("base_PixelFlags_flag_saturatedCenterAll")) 

55 self.assertFalse(record.get("base_PixelFlags_flag_cr")) 

56 self.assertFalse(record.get("base_PixelFlags_flag_crCenter")) 

57 self.assertFalse(record.get("base_PixelFlags_flag_crCenterAll")) 

58 self.assertFalse(record.get("base_PixelFlags_flag_bad")) 

59 self.assertFalse(record.get("base_PixelFlags_flag_badCenter")) 

60 self.assertFalse(record.get("base_PixelFlags_flag_badCenterAll")) 

61 

62 def testSomeFlags(self): 

63 self.dataset.addSource(100000.0, self.center) 

64 task = self.makeSingleFrameMeasurementTask("base_PixelFlags") 

65 exposure, catalog = self.dataset.realize(10.0, task.schema, randomSeed=0) 

66 # one cr pixel outside the center 

67 cosmicray = exposure.mask.getPlaneBitMask("CR") 

68 x = round(self.center.x) 

69 y = round(self.center.y) 

70 exposure.mask[x+3, y+4] |= cosmicray 

71 # one interpolated pixel near the center 

72 interpolated = exposure.mask.getPlaneBitMask("INTRP") 

73 exposure.mask[self.center] |= interpolated 

74 # all pixels in the center are bad 

75 bad = exposure.mask.getPlaneBitMask("BAD") 

76 exposure.mask[x-1:x+2, y-1:y+2] |= bad 

77 task.run(catalog, exposure) 

78 record = catalog[0] 

79 

80 self.assertFalse(record.get("base_PixelFlags_flag_edge")) 

81 self.assertFalse(record.get("base_PixelFlags_flag_edgeCenter")) 

82 self.assertFalse(record.get("base_PixelFlags_flag_edgeCenterAll")) 

83 self.assertTrue(record.get("base_PixelFlags_flag_cr")) 

84 self.assertFalse(record.get("base_PixelFlags_flag_crCenter")) 

85 self.assertFalse(record.get("base_PixelFlags_flag_crCenterAll")) 

86 self.assertTrue(record.get("base_PixelFlags_flag_interpolated")) 

87 self.assertTrue(record.get("base_PixelFlags_flag_interpolatedCenter")) 

88 self.assertFalse(record.get("base_PixelFlags_flag_interpolatedCenterAll")) 

89 self.assertTrue(record.get("base_PixelFlags_flag_bad")) 

90 self.assertTrue(record.get("base_PixelFlags_flag_badCenter")) 

91 self.assertTrue(record.get("base_PixelFlags_flag_badCenterAll")) 

92 

93 def testOffimageNonfinite(self): 

94 """Test that flag_offimage and flag_edge are set correctly for 

95 nonfinite sources. 

96 """ 

97 # These three will be explicitly set to have a non-finite centroid; we 

98 # cannot add sources off the image in TestDataset. 

99 self.dataset.addSource(100000, lsst.geom.Point2D(20, 20)) 

100 self.dataset.addSource(100000, lsst.geom.Point2D(20, 100)) 

101 self.dataset.addSource(100000, lsst.geom.Point2D(100, 30)) 

102 task = self.makeSingleFrameMeasurementTask("base_PixelFlags") 

103 exposure, catalog = self.dataset.realize(10.0, task.schema, randomSeed=0) 

104 catalog[0]["slot_Centroid_x"] = np.inf 

105 catalog[1]["slot_Centroid_x"] = -np.inf 

106 catalog[2]["slot_Centroid_y"] = np.nan 

107 task.run(catalog, exposure) 

108 np.testing.assert_array_equal(catalog["base_PixelFlags_flag_offimage"], 

109 np.array([True, True, True])) 

110 np.testing.assert_array_equal(catalog["base_PixelFlags_flag_edge"], 

111 np.array([True, True, True])) 

112 np.testing.assert_array_equal(catalog["base_PixelFlags_flag_edgeCenter"], 

113 np.array([True, True, True])) 

114 np.testing.assert_array_equal(catalog["base_PixelFlags_flag_edgeCenterAll"], 

115 np.array([True, True, True])) 

116 

117 def testOffimage(self): 

118 """Test that sources at the boundary of the image get flag_offimage 

119 and flag_edge[Center[All]] set appropriately. 

120 """ 

121 # These four will be explicitly set to values at the boundary; we 

122 # cannot add sources off the image in TestDataset. 

123 self.dataset.addSource(100000, lsst.geom.Point2D(20, 20)) 

124 self.dataset.addSource(100000, lsst.geom.Point2D(20, 100)) 

125 self.dataset.addSource(100000, lsst.geom.Point2D(40, 100)) 

126 self.dataset.addSource(100000, lsst.geom.Point2D(80, 30)) 

127 

128 task = self.makeSingleFrameMeasurementTask("base_PixelFlags") 

129 exposure, catalog = self.dataset.realize(10.0, task.schema, randomSeed=0) 

130 

131 # Mask 4 pixels at the x-edges 

132 edgeBit = exposure.mask.getPlaneBitMask("EDGE") 

133 exposure.mask.array[:4, :] |= edgeBit 

134 exposure.mask.array[-4:, ] |= edgeBit 

135 # Mask 4 pixels at the y-edges 

136 nodataBit = exposure.mask.getPlaneBitMask("NO_DATA") 

137 exposure.mask.array[:, :4] |= nodataBit 

138 exposure.mask.array[:, -4:] |= nodataBit 

139 

140 # All of these should get edgeCenter and edgeCenterAll. 

141 catalog[0]["slot_Centroid_x"] = -20.5 # on image 

142 catalog[1]["slot_Centroid_x"] = -20.51 # off image 

143 catalog[2]["slot_Centroid_y"] = 129.4 # on image 

144 catalog[3]["slot_Centroid_y"] = 129.50 # off image 

145 

146 task.run(catalog, exposure) 

147 

148 np.testing.assert_array_equal(catalog["base_PixelFlags_flag_offimage"], 

149 np.array([False, True, False, True])) 

150 np.testing.assert_array_equal(catalog["base_PixelFlags_flag_edge"], 

151 np.array([False, True, False, True])) 

152 np.testing.assert_array_equal(catalog["base_PixelFlags_flag_edgeCenter"], 

153 np.array([True, True, True, True])) 

154 np.testing.assert_array_equal(catalog["base_PixelFlags_flag_edgeCenterAll"], 

155 np.array([True, True, True, True])) 

156 

157 def testEdgeFlags(self): 

158 """Test that edge/edgeCenter/edgeCenterAll are set appropriately. 

159 """ 

160 # Half of the Central 3x3 pixels are on the masked EDGE; no offimage or 

161 # edgeCenterAll, but edge and edgeCenter. 

162 self.dataset.addSource(100000, lsst.geom.Point2D(115, 30)) 

163 # Central 3x3 pixels all masked EDGE; no offimage but all edge flags. 

164 self.dataset.addSource(100000, lsst.geom.Point2D(10, 127)) 

165 

166 task = self.makeSingleFrameMeasurementTask("base_PixelFlags") 

167 exposure, catalog = self.dataset.realize(10.0, task.schema, randomSeed=0) 

168 

169 # Mask 4 pixels at the x-edges 

170 edgeBit = exposure.mask.getPlaneBitMask("EDGE") 

171 exposure.mask.array[:4, :] |= edgeBit 

172 exposure.mask.array[-4:, ] |= edgeBit 

173 # Mask 4 pixels at the y-edges 

174 nodataBit = exposure.mask.getPlaneBitMask("NO_DATA") 

175 exposure.mask.array[:, :4] |= nodataBit 

176 exposure.mask.array[:, -4:] |= nodataBit 

177 

178 task.run(catalog, exposure) 

179 

180 np.testing.assert_array_equal(catalog["base_PixelFlags_flag_offimage"], 

181 np.array([False, False])) 

182 np.testing.assert_array_equal(catalog["base_PixelFlags_flag_edge"], 

183 np.array([True, True])) 

184 np.testing.assert_array_equal(catalog["base_PixelFlags_flag_edgeCenter"], 

185 np.array([True, True])) 

186 np.testing.assert_array_equal(catalog["base_PixelFlags_flag_edgeCenterAll"], 

187 np.array([False, True])) 

188 

189 

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

191 pass 

192 

193 

194def setup_module(module): 

195 lsst.utils.tests.init() 

196 

197 

198if __name__ == "__main__": 198 ↛ 199line 198 didn't jump to line 199, because the condition on line 198 was never true

199 lsst.utils.tests.init() 

200 unittest.main()