LSSTApplications  10.0-2-g4f67435,11.0.rc2+1,11.0.rc2+12,11.0.rc2+3,11.0.rc2+4,11.0.rc2+5,11.0.rc2+6,11.0.rc2+7,11.0.rc2+8
LSSTDataManagementBasePackage
interpImage.py
Go to the documentation of this file.
1 from __future__ import division, absolute_import
2 #
3 # LSST Data Management System
4 # Copyright 2008, 2009, 2010, 2011, 2012 LSST Corporation.
5 #
6 # This product includes software developed by the
7 # LSST Project (http://www.lsst.org/).
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 LSST License Statement and
20 # the GNU General Public License along with this program. If not,
21 # see <http://www.lsstcorp.org/LegalNotices/>.
22 #
23 import lsst.pex.config as pexConfig
24 import lsst.meas.algorithms as measAlg
25 import lsst.pipe.base as pipeBase
26 import lsst.ip.isr as ipIsr
27 
28 __all__ = ["InterpImageConfig", "InterpImageTask"]
29 
30 class InterpImageConfig(pexConfig.Config):
31  """Config for InterpImageTask
32  """
33  modelPsf = measAlg.GaussianPsfFactory.makeField(doc = "Model Psf factory")
34 
35 class InterpImageTask(pipeBase.Task):
36  """Interpolate over bad image pixels
37  """
38  ConfigClass = InterpImageConfig
39 
40  @pipeBase.timeMethod
41  def interpolateOnePlane(self, maskedImage, planeName, fwhmPixels=None):
42  """Interpolate over one mask plane, in place
43 
44  Note that the interpolation code in meas_algorithms currently
45  doesn't use the input PSF (though it's a required argument),
46  so it's not important to set the input PSF parameters exactly.
47 
48  @param[in,out] maskedImage: MaskedImage over which to interpolate over edge pixels
49  @param[in] planeName: mask plane over which to interpolate
50  @param[in] fwhmPixels: FWHM of core star (pixels); if None then the default is used
51  """
52  self.log.info("Interpolate over %s pixels" % (planeName,))
53  psfModel = self.config.modelPsf.apply(fwhm=fwhmPixels)
54 
55  nanDefectList = ipIsr.getDefectListFromMask(maskedImage, planeName, growFootprints=0)
56  measAlg.interpolateOverDefects(maskedImage, psfModel, nanDefectList, 0.0)
void interpolateOverDefects(MaskedImageT &image, lsst::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:2058