casa6

//# Copyright (C) 1997-2010

//# 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:  $

#include <synthesis/MeasurementEquations/MatrixCleaner.h>

#include <synthesis/MeasurementEquations/ImageMSCleaner.h>

#include <images/Images/PagedImage.h>

#include <coordinates/Coordinates/CoordinateSystem.h>

#include <coordinates/Coordinates/CoordinateUtil.h>

​

using namespace casacore;

namespace casa { //# NAMESPACE CASA - BEGIN

​

  ImageMSCleaner::ImageMSCleaner(): psf_p(0), dirty_p(0), mask_p(0), nPsfChan_p(0), 

                    nImChan_p(0), nPsfPol_p(0), nImPol_p(0), chanAxis_p(-1), 

                    polAxis_p(-1), nMaskChan_p(0), nMaskPol_p(0), maskThresh_p(0.9), maxResidual_p(0.0)

 {

​

​

  }

  ImageMSCleaner::ImageMSCleaner(ImageInterface<Float>& psf, 

                 ImageInterface<Float>& dirty): psf_p(&psf), 

                                      dirty_p(&dirty), mask_p(0), 

                                      nMaskChan_p(0), nMaskPol_p(0),

                                maskThresh_p(0.9), maxResidual_p(0.0){

    chanAxis_p=CoordinateUtil::findSpectralAxis(dirty_p->coordinates());

    Vector<Stokes::StokesTypes> whichPols;

    polAxis_p=CoordinateUtil::findStokesAxis(whichPols, dirty_p->coordinates());

    if(chanAxis_p > -1){

      nPsfChan_p=psf_p->shape()(chanAxis_p);

      nImChan_p=dirty_p->shape()(chanAxis_p);

    }

    if(polAxis_p > -1){

      nPsfPol_p=psf_p->shape()(polAxis_p);

      nImPol_p=dirty_p->shape()(polAxis_p);

    }

  }

​

​

  ImageMSCleaner::ImageMSCleaner(const ImageMSCleaner& other){

    operator=(other);