casa6

//# SepImageConvolver.cc:  separable convolution of an image

//# Copyright (C) 1995,1996,1997,1998,1999,2000,2001,2002,2003

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

//#

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

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

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

//# option) any later version.

//#

//# This library 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 Library General Public

//# License for more details.

//#

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

//# along with this library; 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: SepImageConvolver.tcc 20615 2009-06-09 02:16:01Z Malte.Marquarding $

//   

​

#include <imageanalysis/ImageAnalysis/SepImageConvolver.h>

​

#include <casacore/casa/aips.h>

​

#include <casacore/casa/Arrays/Vector.h>

#include <casacore/casa/Arrays/ArrayMath.h>

#include <casacore/casa/Containers/Block.h>

#include <casacore/casa/Exceptions/Error.h>

#include <casacore/casa/Logging/LogIO.h>

#include <casacore/scimath/Mathematics/Convolver.h>

#include <casacore/casa/Quanta/UnitMap.h>

#include <casacore/casa/Quanta/Quantum.h>

#include <casacore/casa/BasicSL/String.h>

​

#include <casacore/coordinates/Coordinates/CoordinateSystem.h>

#include <casacore/images/Images/PagedImage.h>

#include <casacore/images/Regions/ImageRegion.h>

#include <casacore/images/Images/SubImage.h>

#include <casacore/images/Images/TempImage.h>

#include <casacore/lattices/Lattices/LatticeIterator.h>

#include <casacore/lattices/Lattices/LatticeUtilities.h>

#include <casacore/lattices/Lattices/TiledLineStepper.h>

#include <casacore/scimath/Mathematics/VectorKernel.h>

#include <casacore/casa/System/ProgressMeter.h>

​

#include <sstream>

#include <iomanip>

​

​

namespace casa { //# NAMESPACE CASA - BEGIN

​

template <class T> 

SepImageConvolver<T>::SepImageConvolver(

    const casacore::ImageInterface<T>& image, casacore::LogIO &os, casacore::Bool showProgress

)

: itsImagePtr(image.cloneII()),

  itsOs(os),

  itsAxes(0),

  itsVectorKernels(0),

  itsShowProgress(showProgress) {}

​

template <class T>

SepImageConvolver<T>::SepImageConvolver(const SepImageConvolver<T> &other)

: itsImagePtr(0) {

   operator=(other);

}

​

​

template <class T> 

SepImageConvolver<T>::~SepImageConvolver () {

   delete itsImagePtr;

   itsImagePtr = 0;

   const casacore::uInt n = itsVectorKernels.nelements();

   for (casacore::uInt i=0; i<n; i++) {

      delete itsVectorKernels[i];

      itsVectorKernels[i] = 0;

   }

}

​

template <class T>

SepImageConvolver<T> &SepImageConvolver<T>::operator=(

    const SepImageConvolver<T> &other

) {