LSST Applications g044012fb7c+6976b5ec80,g04a91732dc+88a5fc122b,g07dc498a13+7e3c5f68a2,g114c6a66ad+09472d7a76,g1409bbee79+7e3c5f68a2,g1a7e361dbc+7e3c5f68a2,g1fd858c14a+3a43eabc0e,g35bb328faa+fcb1d3bbc8,g3bd4b5ce2c+2647bb081c,g4e0f332c67+5d362be553,g53246c7159+fcb1d3bbc8,g5477a8d5ce+b19c77c7ae,g58d0cdf3ff+4a2e102ff8,g60b5630c4e+09472d7a76,g623d845a50+09472d7a76,g6f0c2978f1+fcf1c0bcd6,g71fabbc107+09472d7a76,g75b6c65c88+d0b1dc44cc,g78460c75b0+2f9a1b4bcd,g786e29fd12+cf7ec2a62a,g7b71ed6315+fcb1d3bbc8,g8852436030+349c7e81d4,g89139ef638+7e3c5f68a2,g9125e01d80+fcb1d3bbc8,g95236ca021+f7a31438ed,g989de1cb63+7e3c5f68a2,g9f33ca652e+f17d666fbc,gaaedd4e678+7e3c5f68a2,gabe3b4be73+1e0a283bba,gb1101e3267+f870f33517,gb58c049af0+f03b321e39,gc99c83e5f0+76d20ab76d,gcf25f946ba+349c7e81d4,gd0fa69b896+f3a65fa83c,gd6cbbdb0b4+c8606af20c,gde0f65d7ad+5bd27d919f,ge278dab8ac+932305ba37,gfba249425e+fcb1d3bbc8,w.2025.07
LSST Data Management Base Package
Loading...
Searching...
No Matches
reinterpolate_pixels.py
Go to the documentation of this file.
1# This file is part of meas_algorithms.
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
22__all__ = ("ReinterpolatePixelsConfig", "ReinterpolatePixelsTask")
23
24import itertools
25
26import lsst.afw.detection as afwDetection
27import lsst.afw.math as afwMath
28import lsst.pex.config as pexConfig
29import lsst.pipe.base as pipeBase
30from . import Defect, interpolateOverDefects
31
32
33class ReinterpolatePixelsConfig(pexConfig.Config):
34 """Config for ReinterpolatePixelsTask"""
35
36 maskNameList = pexConfig.ListField[str](
37 doc="Names of mask planes whose image-plane values should be interpolated.",
38 default=["SAT"],
39 )
40
41 fallbackValue = pexConfig.Field[float](
42 doc="Value of last resort for interpolation. If not provided (or set to None), it is assigned the "
43 "clipped mean value of the exposure image.",
44 default=None,
45 optional=True,
46 )
47
48
49class ReinterpolatePixelsTask(pipeBase.Task):
50 """Task to reinterpolate pixels"""
51
52 ConfigClass = ReinterpolatePixelsConfig
53 _DefaultName = "reinterpolatePixels"
54
55 def run(self, exposure):
56 """Reinterpolate pixels over ``exposure``.
57
58 Parameters
59 ----------
60 exposure : `lsst.afw.image.Exposure`
61 Exposure to interpolate over.
62
63 Notes
64 -----
65 ``exposure`` is modified in place and will become the reinterpolated exposure.
66 """
67 mask = exposure.mask
68
69 thresh = afwDetection.Threshold(
70 mask.getPlaneBitMask(self.config.maskNameList), afwDetection.Threshold.BITMASK
71 )
72 fpSet = afwDetection.FootprintSet(mask, thresh)
73 defectList = self._fromFootprintList(fpSet.getFootprints())
74 self._interpolateDefectList(exposure, defectList)
75
76 def _fromFootprintList(self, fpList):
77 """Compute a defect list from a footprint list.
78
79 Parameters
80 ----------
81 fpList : `list` of `lsst.afw.detection.Footprint`
82 Footprint list to process.
83
84 Returns
85 -------
86 defects : `list` of `lsst.meas.algorithms.Defect`
87 List of defects.
88
89 Notes
90 -----
91 By using `itertools.chain.from_iterable()`, the nested iterables are merged into a flat iterable,
92 allowing convenient iteration over all `lsst.meas.algorithms.Defect` objects.
93 """
94 return list(
95 itertools.chain.from_iterable(map(Defect, afwDetection.footprintToBBoxList(fp)) for fp in fpList)
96 )
97
98 def _interpolateDefectList(self, exposure, defectList):
99 """Interpolate over defects specified in a defect list.
100
101 Parameters
102 ----------
103 exposure : `lsst.afw.image.Exposure`
104 Exposure to process.
105 defectList : `list` of `lsst.meas.algorithms.Defect`
106 List of defects to interpolate over.
107
108 Notes
109 -----
110 ``exposure`` is modified in place and will become the reinterpolated exposure.
111 When reinterpolating following the interpolation in ip_isr, be aware that the masks are intentionally
112 left grown as a side-effect of that interpolation.
113 """
114
115 # Although `interpolateOverDefects` requires the
116 # `lsst.afw.detection.Psf` argument to be passed, it does not actually
117 # utilize it. A dummy `lsst.afw.detection.Psf` object must still be
118 # provided to prevent a TypeError.
119 dummyPsf = afwDetection.Psf()
120 if self.config.fallbackValue is None:
121 fallbackValue = afwMath.makeStatistics(exposure.image, afwMath.MEANCLIP).getValue()
122 else:
123 fallbackValue = self.config.fallbackValue
124 if "INTRP" not in exposure.mask.getMaskPlaneDict():
125 exposure.mask.addMaskPlane("INTRP")
126 interpolateOverDefects(exposure.maskedImage, dummyPsf, defectList, fallbackValue, True)
Statistics makeStatistics(lsst::afw::image::Image< Pixel > const &img, lsst::afw::image::Mask< image::MaskPixel > const &msk, int const flags, StatisticsControl const &sctrl=StatisticsControl())
Handle a watered-down front-end to the constructor (no variance)
Definition Statistics.h:361
void interpolateOverDefects(MaskedImageT &image, afw::detection::Psf const &psf, std::vector< Defect::Ptr > &badList, double fallbackValue=0.0, bool useFallbackValueAtEdge=false)
Process a set of known bad pixels in an image.
Definition Interp.cc:2047