lsst::coadd::utils Namespace Reference

Namespaces
	coaddDataIdContainer
	version

Functions
template<typename CoaddPixelT , typename WeightPixelT >
lsst::geom::Box2I	addToCoadd (lsst::afw::image::Image< CoaddPixelT > &coadd, lsst::afw::image::Image< WeightPixelT > &weightMap, lsst::afw::image::Image< CoaddPixelT > const &image, WeightPixelT weight)
	add good pixels from an image to a coadd and associated weight map More...
template<typename CoaddPixelT , typename WeightPixelT >
lsst::geom::Box2I	addToCoadd (lsst::afw::image::MaskedImage< CoaddPixelT, lsst::afw::image::MaskPixel, lsst::afw::image::VariancePixel > &coadd, lsst::afw::image::Image< WeightPixelT > &weightMap, lsst::afw::image::MaskedImage< CoaddPixelT, lsst::afw::image::MaskPixel, lsst::afw::image::VariancePixel > const &maskedImage, lsst::afw::image::MaskPixel const badPixelMask, WeightPixelT weight)
	add good pixels from a masked image to a coadd image and associated weight map More...
template<typename ImagePixelT >
int	copyGoodPixels (lsst::afw::image::Image< ImagePixelT > &destImage, lsst::afw::image::Image< ImagePixelT > const &srcImage)
	copy good pixels from one image to another More...
template<typename ImagePixelT >
int	copyGoodPixels (lsst::afw::image::MaskedImage< ImagePixelT, lsst::afw::image::MaskPixel, lsst::afw::image::VariancePixel > &destImage, lsst::afw::image::MaskedImage< ImagePixelT, lsst::afw::image::MaskPixel, lsst::afw::image::VariancePixel > const &srcImage, lsst::afw::image::MaskPixel const badPixelMask)
	copy good pixels from one masked image to another More...
template<typename WeightPixelT >
void	setCoaddEdgeBits (lsst::afw::image::Mask< lsst::afw::image::MaskPixel > &coaddMask, lsst::afw::image::Image< WeightPixelT > const &weightMap)
	set edge bits of coadd mask based on weight map More...
	PYBIND11_MODULE (addToCoadd, mod)
	PYBIND11_MODULE (copyGoodPixels, mod)
	PYBIND11_MODULE (setCoaddEdgeBits, mod)

Function Documentation

addToCoadd() [1/2]

template<typename CoaddPixelT , typename WeightPixelT >

lsst::geom::Box2I lsst::coadd::utils::addToCoadd	(	lsst::afw::image::Image< CoaddPixelT > &	coadd,
		lsst::afw::image::Image< WeightPixelT > &	weightMap,
		lsst::afw::image::Image< CoaddPixelT > const &	image,
		WeightPixelT	weight
	)

add good pixels from an image to a coadd and associated weight map

The images are assumed to be registered to the same wcs and parent origin, thus: coadd[i+coadd.x0, j+coadd.y0] += image[i+image.x0, j+image.y0] weightMap[i+weightMap.x0, j+weightMap.y0] += weight for all good image pixels that overlap a coadd pixel. Good pixels are those that are not NaN (thus they do include +/- inf).

Returns

overlapBBox: overlapping bounding box, relative to parent image (hence xy0 is taken into account)

Exceptions

pexExcept::InvalidParameterError

if coadd and weightMap dimensions or xy0 do not match.

Parameters

[in,out]	coadd	coadd to be modified
[in,out]	weightMap	weight map to be modified; this is the sum of weights of all images contributing each pixel of the coadd
	image	image to add to coadd
	weight	relative weight of this image

Definition at line 127 of file addToCoadd.cc.

  {
    typedef lsst::afw::image::Image<CoaddPixelT> Image;
    return addToCoaddImpl<Image, WeightPixelT, CheckKnownValue>(coadd, weightMap, image, 0x0, weight);
}

addToCoadd() [2/2]

template<typename CoaddPixelT , typename WeightPixelT >

lsst::geom::Box2I lsst::coadd::utils::addToCoadd	(	lsst::afw::image::MaskedImage< CoaddPixelT, lsst::afw::image::MaskPixel, lsst::afw::image::VariancePixel > &	coadd,
		lsst::afw::image::Image< WeightPixelT > &	weightMap,
		lsst::afw::image::MaskedImage< CoaddPixelT, lsst::afw::image::MaskPixel, lsst::afw::image::VariancePixel > const &	maskedImage,
		lsst::afw::image::MaskPixel const	badPixelMask,
		WeightPixelT	weight
	)

add good pixels from a masked image to a coadd image and associated weight map

The images are assumed to be registered to the same wcs and parent origin, thus: coadd[i+coadd.x0, j+coadd.y0] += image[i+image.x0, j+image.y0] weightMap[i+weightMap.x0, j+weightMap.y0] += weight for all good image pixels that overlap a coadd pixel. Good pixels are those for which mask & badPixelMask == 0.

Returns

overlapBBox: overlapping bounding box, relative to parent image (hence xy0 is taken into account)

Exceptions

pexExcept::InvalidParameterError

if coadd and weightMap dimensions or xy0 do not match.

Parameters

[in,out]	coadd	coadd to be modified
[in,out]	weightMap	weight map to be modified; this is the sum of weights of all images contributing each pixel of the coadd
	maskedImage	masked image to add to coadd
	badPixelMask	skip input pixel if input mask & badPixelMask !=0
	weight	relative weight of this image

Definition at line 139 of file addToCoadd.cc.

  {
    typedef lsst::afw::image::MaskedImage<CoaddPixelT> Image;
    return addToCoaddImpl<Image,WeightPixelT,CheckMask>(coadd, weightMap, maskedImage, badPixelMask, weight);
}

copyGoodPixels() [1/2]

template<typename ImagePixelT >

int lsst::coadd::utils::copyGoodPixels	(	lsst::afw::image::Image< ImagePixelT > &	destImage,
		lsst::afw::image::Image< ImagePixelT > const &	srcImage
	)

copy good pixels from one image to another

Good pixels are those that are not NaN (thus they do include +/- inf).

Only the overlapping pixels (relative to the parent) are copied; thus the images do not have to be the same size.

Returns

number of pixels copied

Parameters

[in,out]	destImage	image to be modified
	srcImage	image to copy

Definition at line 116 of file copyGoodPixels.cc.

  {
    typedef lsst::afw::image::Image<ImagePixelT> Image;
    return copyGoodPixelsImpl<Image, CheckKnownValue>(destImage, srcImage, 0x0);
}

copyGoodPixels() [2/2]

template<typename ImagePixelT >

int lsst::coadd::utils::copyGoodPixels	(	lsst::afw::image::MaskedImage< ImagePixelT, lsst::afw::image::MaskPixel, lsst::afw::image::VariancePixel > &	destImage,
		lsst::afw::image::MaskedImage< ImagePixelT, lsst::afw::image::MaskPixel, lsst::afw::image::VariancePixel > const &	srcImage,
		lsst::afw::image::MaskPixel const	badPixelMask
	)

copy good pixels from one masked image to another

Good pixels are those for which mask & badPixelMask == 0.

Only the overlapping pixels (relative to the parent) are copied; thus the images do not have to be the same size.

Returns

number of pixels copied

Parameters

[in,out]	destImage	image to be modified
	srcImage	image to copy
	badPixelMask	skip input pixel if src mask & badPixelMask != 0

Definition at line 126 of file copyGoodPixels.cc.

  {
    typedef lsst::afw::image::MaskedImage<ImagePixelT> Image;
    return copyGoodPixelsImpl<Image, CheckMask>(destImage, srcImage, badPixelMask);
}

PYBIND11_MODULE() [1/3]

lsst::coadd::utils::PYBIND11_MODULE	(	addToCoadd	,
		mod
	)

Definition at line 56 of file addToCoadd.cc.

                                 {
    py::module::import("lsst.afw.geom");
    py::module::import("lsst.afw.image");
 
    declareAddToCoadd<double, double>(mod);
    declareAddToCoadd<double, float>(mod);
    declareAddToCoadd<double, int>(mod);
    declareAddToCoadd<double, std::uint16_t>(mod);
    declareAddToCoadd<float, double>(mod);
    declareAddToCoadd<float, float>(mod);
    declareAddToCoadd<float, int>(mod);
    declareAddToCoadd<float, std::uint16_t>(mod);
}

PYBIND11_MODULE() [2/3]

lsst::coadd::utils::PYBIND11_MODULE	(	copyGoodPixels	,
		mod
	)

Definition at line 52 of file copyGoodPixels.cc.

                                     {
    py::module::import("lsst.afw.image");
 
    declareCopyGoodPixels<double>(mod);
    declareCopyGoodPixels<float>(mod);
    declareCopyGoodPixels<int>(mod);
    declareCopyGoodPixels<std::uint16_t>(mod);
}

PYBIND11_MODULE() [3/3]

lsst::coadd::utils::PYBIND11_MODULE	(	setCoaddEdgeBits	,
		mod
	)

Definition at line 48 of file setCoaddEdgeBits.cc.

                                       {
    py::module::import("lsst.afw.image");
 
    declareSetCoaddEdgeBits<double>(mod);
    declareSetCoaddEdgeBits<float>(mod);
    declareSetCoaddEdgeBits<int>(mod);
    declareSetCoaddEdgeBits<std::uint16_t>(mod);
}

setCoaddEdgeBits()

template<typename WeightPixelT >

void lsst::coadd::utils::setCoaddEdgeBits	(	lsst::afw::image::Mask< lsst::afw::image::MaskPixel > &	coaddMask,
		lsst::afw::image::Image< WeightPixelT > const &	weightMap
	)

set edge bits of coadd mask based on weight map

Set pixels in the image to the edge pixel when the corresponding pixel in the weight map is zero. The edge pixel is image=nan, variance=inf, mask=NO_DATA for masked images and image=nan for plain images.

Exceptions

pexExcept::InvalidParameterError

if the dimensions of coaddMask and weightMap do not match.

Parameters

[in,out]	coaddMask	mask of coadd
	weightMap	weight map

Definition at line 42 of file setCoaddEdgeBits.cc.

  {
    typedef afwImage::Mask<afwImage::MaskPixel>::x_iterator MaskXIter;
    typedef typename afwImage::Image<WeightPixelT>::const_x_iterator WeightMapConstXIter;
 
    if (coaddMask.getDimensions() != weightMap.getDimensions()) {
        throw LSST_EXCEPT(pexExcept::InvalidParameterError,
            (boost::format("coaddMask and weightMap dimensions differ: %dx%d != %dx%d") %
            coaddMask.getWidth() % coaddMask.getHeight() % weightMap.getWidth() % weightMap.getHeight()
            ).str());
    }
 
    afwImage::MaskPixel const edgeMask = afwImage::Mask<afwImage::MaskPixel>::getPlaneBitMask("NO_DATA");
 
    // Set the pixels row by row, to avoid repeated checks for end-of-row
    for (int y = 0, endY = weightMap.getHeight(); y != endY; ++y) {
        WeightMapConstXIter weightMapPtr = weightMap.row_begin(y);
        WeightMapConstXIter const weightMapEndPtr = weightMap.row_end(y);
        MaskXIter coaddMaskPtr = coaddMask.row_begin(y);
        for (; weightMapPtr != weightMapEndPtr; ++weightMapPtr, ++coaddMaskPtr) {
            if (*weightMapPtr == 0) {
                (*coaddMaskPtr) = (*coaddMaskPtr) | edgeMask;
            }
        }
    }
}