38#include "boost/timer.hpp"
61template <
typename PixelT>
65 unsigned int rows =
img.getHeight();
66 unsigned int cols =
img.getWidth();
67 Eigen::MatrixXd
M = Eigen::MatrixXd::Zero(
rows,
cols);
68 for (
int y = 0;
y !=
img.getHeight(); ++
y) {
71 ptr !=
img.row_end(
y); ++ptr, ++
x) {
83 unsigned int rows =
mask.getHeight();
85 Eigen::MatrixXi
M = Eigen::MatrixXi::Zero(
rows,
cols);
86 for (
int y = 0;
y !=
mask.getHeight(); ++
y) {
89 ptr !=
mask.row_end(
y); ++ptr, ++
x) {
114template <
typename PixelT,
typename BackgroundT>
143 double time = t.elapsed();
144 LOGL_DEBUG(
"TRACE4.ip.diffim.convolveAndSubtract",
145 "Total compute time to convolve and subtract : %.2f s", time);
165template <
typename PixelT,
typename BackgroundT>
196 double time = t.elapsed();
197 LOGL_DEBUG(
"TRACE4.ip.diffim.convolveAndSubtract",
198 "Total compute time to convolve and subtract : %.2f s", time);
220#define p_INSTANTIATE_convolveAndSubtract(TEMPLATE_IMAGE_T, TYPE) \
222 lsst::afw::image::MaskedImage<TYPE> convolveAndSubtract( \
223 lsst::afw::image::TEMPLATE_IMAGE_T<TYPE> const& templateImage, \
224 lsst::afw::image::MaskedImage<TYPE> const& scienceMaskedImage, \
225 lsst::afw::math::Kernel const& convolutionKernel, \
230 afwImage::MaskedImage<TYPE> convolveAndSubtract( \
231 lsst::afw::image::TEMPLATE_IMAGE_T<TYPE> const& templateImage, \
232 lsst::afw::image::MaskedImage<TYPE> const& scienceMaskedImage, \
233 lsst::afw::math::Kernel const& convolutionKernel, \
234 lsst::afw::math::Function2<double> const& backgroundFunction, \
237#define INSTANTIATE_convolveAndSubtract(TYPE) \
238p_INSTANTIATE_convolveAndSubtract(Image, TYPE) \
239p_INSTANTIATE_convolveAndSubtract(MaskedImage, TYPE)
#define LOGL_DEBUG(logger, message...)
#define INSTANTIATE_convolveAndSubtract(TYPE)
x_iterator row_begin(int y) const
void setDoNormalize(bool doNormalize)
Asseses the quality of a candidate given a spatial kernel and background model.
An include file to include the header files for lsst::ip::diffim.
void convolve(OutImageT &convolvedImage, InImageT const &inImage, KernelT const &kernel, ConvolutionControl const &convolutionControl=ConvolutionControl())
Eigen::MatrixXd imageToEigenMatrix(lsst::afw::image::Image< PixelT > const &img)
Turns a 2-d Image into a 2-d Eigen Matrix.
Eigen::MatrixXi maskToEigenMatrix(lsst::afw::image::Mask< lsst::afw::image::MaskPixel > const &mask)
lsst::afw::image::MaskedImage< PixelT > convolveAndSubtract(lsst::afw::image::MaskedImage< PixelT > const &templateImage, lsst::afw::image::MaskedImage< PixelT > const &scienceMaskedImage, lsst::afw::math::Kernel const &convolutionKernel, BackgroundT background, bool invert=true)
Execute fundamental task of convolving template and subtracting it from science image.