casa6

//# MomentClip.cc:

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

//# 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: MomentCalculator.tcc 19940 2007-02-27 05:35:22Z Malte.Marquarding $

//

​

#include <imageanalysis/ImageAnalysis/MomentClip.h>

​

namespace casa {

​

template <class T> MomentClip<T>::MomentClip(

    shared_ptr<casacore::Lattice<T>> pAncilliaryLattice, MomentsBase<T>& iMom,

    casacore::LogIO& os, const casacore::uInt nLatticeOut

) : _ancilliaryLattice(pAncilliaryLattice), iMom_p(iMom), os_p(os) {

    selectMoments_p = this->selectMoments(iMom_p);

    constructorCheck(calcMoments_p, calcMomentsMask_p, selectMoments_p, nLatticeOut);

    casacore::Int momAxis = this->momentAxis(iMom_p);

    if (_ancilliaryLattice) {

        sliceShape_p.resize(_ancilliaryLattice->ndim());

        sliceShape_p = 1;

        sliceShape_p(momAxis) = _ancilliaryLattice->shape()(momAxis);

    }

    //this->yAutoMinMax(yMinAuto_p, yMaxAuto_p, iMom_p);

​

    this->selectRange(range_p, doInclude_p, doExclude_p, iMom_p);

​

    this->costlyMoments(iMom_p, doMedianI_p, doMedianV_p, doAbsDev_p);

​

    // Are we computing coordinate-dependent moments.  If so

    // precompute coordinate vector if moment axis separable

​

    this->setCoordinateSystem (cSys_p, iMom_p);

    this->doCoordCalc(doCoordProfile_p, doCoordRandom_p, iMom_p);

    this->setUpCoords(

        iMom_p, pixelIn_p, worldOut_p, sepWorldCoord_p, os_p,

        integratedScaleFactor_p, cSys_p, doCoordProfile_p,

        doCoordRandom_p

    );

​

    // What is the axis type of the moment axis

​