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
KernelFunctions.cc
Go to the documentation of this file.
1 // -*- LSST-C++ -*-
2 
3 /*
4  * LSST Data Management System
5  * Copyright 2008, 2009, 2010 LSST Corporation.
6  *
7  * This product includes software developed by the
8  * LSST Project (http://www.lsst.org/).
9  *
10  * This program is free software: you can redistribute it and/or modify
11  * it under the terms of the GNU General Public License as published by
12  * the Free Software Foundation, either version 3 of the License, or
13  * (at your option) any later version.
14  *
15  * This program is distributed in the hope that it will be useful,
16  * but WITHOUT ANY WARRANTY; without even the implied warranty of
17  * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
18  * GNU General Public License for more details.
19  *
20  * You should have received a copy of the LSST License Statement and
21  * the GNU General Public License along with this program. If not,
22  * see <http://www.lsstcorp.org/LegalNotices/>.
23  */
24 
32 #include <iostream>
33 
34 #include "boost/format.hpp"
35 
37 
45 void
47  lsst::afw::math::Kernel const &kernel,
48  bool doNormalize,
49  double xPos,
50  double yPos,
51  std::string pixelFmt
52 ) {
53  typedef lsst::afw::math::Kernel::Pixel Pixel;
54 
56  double kSum = kernel.computeImage(kImage, doNormalize, xPos, yPos);
57 
58  for (int y = kImage.getHeight() - 1; y >= 0; --y) {
59  for (lsst::afw::image::Image<Pixel>::const_x_iterator ptr = kImage.row_begin(y);
60  ptr != kImage.row_end(y); ++ptr) {
61  std::cout << boost::format(pixelFmt) % *ptr << " ";
62  }
63  std::cout << std::endl;
64  }
65 
66  if (doNormalize && std::abs(static_cast<double>(kSum) - 1.0) > 1.0e-5) {
67  std::cout << boost::format("Warning! Sum of all pixels = %9.5f != 1.0\n") % kSum;
68  }
69  std::cout << std::endl;
70 }
int y
void printKernel(lsst::afw::math::Kernel const &kernel, bool doNormalize, double x=0, double y=0, std::string pixelFmt="%7.3f")
Print the pixel values of a Kernel to std::cout.
_const_view_t::x_iterator const_x_iterator
A const iterator for traversing the pixels in a row.
Definition: Image.h:158
geom::Extent2I const getDimensions() const
Return the Kernel&#39;s dimensions (width, height)
Definition: Kernel.h:226
Kernels are used for convolution with MaskedImages and (eventually) Images.
Definition: Kernel.h:134
Utility functions for kernels.
double computeImage(lsst::afw::image::Image< Pixel > &image, bool doNormalize, double x=0.0, double y=0.0) const
Compute an image (pixellized representation of the kernel) in place.
Definition: Kernel.cc:94