Coverage for python/lsst/afw/detection/_heavyFootprintContinued.py: 41%
25 statements
« prev ^ index » next coverage.py v7.2.7, created at 2023-06-15 02:47 -0700
« prev ^ index » next coverage.py v7.2.7, created at 2023-06-15 02:47 -0700
1# This file is part of afw.
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#
19import numpy as np
21from lsst.utils import TemplateMeta
23from ._detection import HeavyFootprintI, HeavyFootprintF, HeavyFootprintD, HeavyFootprintU
25__all__ = [] # import this module only for its side effects
28class HeavyFootprint(metaclass=TemplateMeta): # noqa: F811
29 def addTo(self, image):
30 """Add this heavy footprint to an image.
32 Parameters
33 ----------
34 image : `lsst.afw.image`
35 """
36 indices = self.spans.indices()
37 image.array[indices[0, :] - image.getY0(),
38 indices[1, :] - image.getX0()] += self.getImageArray()
40 def subtractFrom(self, image):
41 """Subtract this heavy footprint from an image.
43 Parameters
44 ----------
45 image : `lsst.afw.image`
46 """
47 indices = self.spans.indices()
48 image.array[indices[0, :] - image.getY0(),
49 indices[1, :] - image.getX0()] -= self.getImageArray()
51 def extractImage(self, fill=np.nan, bbox=None, imageType=None):
52 """Create a 2D image of a HeavyFootprint
54 Parameters
55 ----------
56 fill: number
57 Number to fill the pixels in the image that are not
58 contained in `heavy`.
59 bbox : `Box2I`
60 Bounding box of the output image.
61 imageType : `type`
62 This should be either a `MaskedImage` or `Image` and describes
63 the type of the output image.
64 If `imageType` is `None` then `Image` will be used.
66 Returns
67 -------
68 image : `lsst.afw.image.MaskedImage` or `lsst.afw.image.Image`
69 An image defined by `bbox` and padded with `fill` that
70 contains the projected flux in `heavy`.
71 """
72 # Prevent circular import
73 from lsst.afw.image import Image
74 if imageType is None:
75 imageType = Image
76 if bbox is None:
77 bbox = self.getBBox()
78 image = imageType(bbox, dtype=self.getImageArray().dtype)
79 image.set(fill)
80 self.insert(image)
81 return image
84HeavyFootprint.register(np.int32, HeavyFootprintI)
85HeavyFootprint.register(np.float32, HeavyFootprintF)
86HeavyFootprint.register(np.float64, HeavyFootprintD)
87HeavyFootprint.register(np.uint16, HeavyFootprintU)