casa6

//# tSubImage.cc: Test program for class SubImage

//# Copyright (C) 1998,1999,2000,2001,2003

//# Associated Universities, Inc. Washington DC, USA.

//#

//# This program is free software; you can redistribute it and/or modify it

//# under the terms of the GNU General Public License as published by the Free

//# Software Foundation; either version 2 of the License, or (at your option)

//# any later version.

//#

//# This program is distributed in the hope that it will be useful, but WITHOUT

//# ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or

//# FITNESS FOR A PARTICULAR PURPOSE.  See the GNU General Public License for

//# more details.

//#

//# You should have received a copy of the GNU General Public License along

//# with this program; if not, write to the Free Software Foundation, Inc.,

//# 675 Massachusetts Ave, Cambridge, MA 02139, USA.

//#

//# Correspondence concerning AIPS++ should be addressed as follows:

//#        Internet email: aips2-request@nrao.edu.

//#        Postal address: AIPS++ Project Office

//#                        National Radio Astronomy Observatory

//#                        520 Edgemont Road

//#                        Charlottesville, VA 22903-2475 USA

//#

//# $Id: $

​

#include <imageanalysis/ImageAnalysis/ImageMaskedPixelReplacer.h>

​

#include <imageanalysis/ImageAnalysis/SubImageFactory.h>

​

namespace casa {

​

template <class T>

const casacore::String ImageMaskedPixelReplacer<T>::_class = "ImageMaskedPixelReplacer";

​

template <class T>

ImageMaskedPixelReplacer<T>::ImageMaskedPixelReplacer(

    const SPIIT image,

    const casacore::Record *const &region,

    const casacore::String& mask

) : ImageTask<T>(

        image, "", region, "", "", "", mask, "", false

    ),

    _image(image) {

    this->_construct();

}

​

template <class T>

void ImageMaskedPixelReplacer<T>::replace(

    const casacore::String& expr, casacore::Bool updateMask, casacore::Bool verbose

) {

    casacore::LogOrigin lor(_class, __func__);

    *this->_getLog() << lor;

    ThrowIf(

        expr.empty(),

        "You must specify an expression"

    );

​

    // Modify in place by writing to the image passed to the constructor.

    std::shared_ptr<casacore::SubImage<T> > subImage = SubImageFactory<T>::createSubImageRW(

        *_image, *this->_getRegion(), this->_getMask(),

        (verbose ? this->_getLog().get() : 0),