casa6

//# Copyright (C) 2015

//# 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  General Public License as published by

//# the Free Software Foundation; either version 3 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/MatrixNACleaner.h>

#include <synthesis/MeasurementEquations/ImageNACleaner.h>

#include <images/Images/ImageInterface.h>

#include <images/Images/PagedImage.h>

#include <images/Images/TempImage.h>

#include <coordinates/Coordinates/CoordinateSystem.h>

#include <coordinates/Coordinates/CoordinateUtil.h>

using namespace std;

using namespace casacore;

namespace casa { //# NAMESPACE CASA - BEGIN

​

  ImageNACleaner::ImageNACleaner(): psf_p(nullptr), dirty_p(nullptr), mask_p(nullptr), nPsfChan_p(0), 

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

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

 {

​

​

  }

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

                 ImageInterface<Float>& dirty):  mask_p(nullptr), 

                                      nMaskChan_p(0), nMaskPol_p(0),maxResidual_p(0.0){

    psf_p=CountedPtr<ImageInterface<Float> >(&psf, false); 

    dirty_p=CountedPtr<ImageInterface<Float> >(&dirty, false); 

    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);

    }

  }

​

​

  ImageNACleaner::ImageNACleaner(const ImageNACleaner& other){