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
« 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/>.
22import unittest
24import lsst.geom
25import lsst.utils.tests
26import lsst.meas.base.tests
29class PixelFlagsTestCase(lsst.meas.base.tests.AlgorithmTestCase, lsst.utils.tests.TestCase):
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)
38 def tearDown(self):
39 del self.center
40 del self.bbox
41 del self.dataset
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"))
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]
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"))
91class TestMemory(lsst.utils.tests.MemoryTestCase):
92 pass
95def setup_module(module):
96 lsst.utils.tests.init()
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()