Coverage for tests/test_PixelFlags.py: 22%

62 statements  

« prev     ^ index     » next       coverage.py v7.3.2, created at 2023-12-12 12:41 +0000

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 lsst.geom 

25import lsst.utils.tests 

26import lsst.meas.base.tests 

27 

28 

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

30 

31 def setUp(self): 

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

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

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

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

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

37 

38 def tearDown(self): 

39 del self.center 

40 del self.bbox 

41 del self.dataset 

42 

43 def testNoFlags(self): 

44 task = self.makeSingleFrameMeasurementTask("base_PixelFlags") 

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

46 task.run(catalog, exposure) 

47 record = catalog[0] 

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

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

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

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

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

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

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

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

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

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

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

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

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

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

62 

63 def testSomeFlags(self): 

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.assertTrue(record.get("base_PixelFlags_flag_cr")) 

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

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

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

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

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

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

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

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

89 

90 

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

92 pass 

93 

94 

95def setup_module(module): 

96 lsst.utils.tests.init() 

97 

98 

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

100 lsst.utils.tests.init() 

101 unittest.main()