casa6

//# SpectralCollapser.cc: Implementation of class SpectralCollapser

//# 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/SpectralCollapser.h>

#include <imageanalysis/ImageAnalysis/ImageCollapser.h>

​

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

#include <casacore/casa/OS/Directory.h>

#include <casacore/casa/OS/RegularFile.h>

#include <casacore/casa/OS/SymLink.h>

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

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

#include <imageanalysis/ImageAnalysis/ImageMoments.h>

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

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

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

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

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

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

#include <casacore/lattices/LRegions/LCSlicer.h>

​

#include <imageanalysis/ImageAnalysis/ImageMoments.h>

​

using namespace casacore;

namespace casa {

SpectralCollapser::SpectralCollapser(const SPCIIF image):

        _image(image), _log(new LogIO()), _storePath(""){

    _setUp();

}

​

SpectralCollapser::SpectralCollapser(const SPCIIF image, const String storePath):

        _image(image), _log(new LogIO()), _storePath(storePath){

    _setUp();

}

​

SpectralCollapser::~SpectralCollapser(){delete _log;}

​

Bool SpectralCollapser::collapse(const Vector<Float> &specVals, const Float startVal, const Float endVal,

        const String &unit, const SpectralCollapser::CollapseType &collType, const SpectralCollapser::CollapseError &collError, String &outname, String &msg){

​

    //Bool ok;

    String unit_(unit);

    String outnameData;

    String outnameError;

​

    *_log << LogOrigin("SpectralCollapser", "collapse");

​

    if (specVals.size() < 1){

        msg = String("No spectral values provided!");

        *_log << LogIO::WARN << msg << LogIO::POST;

        return false;

    }

​

    // in the unit, replace a "/" with "_p_"

    if (unit_.find(String("/"))!=String::npos)  {

        String::size_type slashPos=unit_.find(String("/"));

        unit_.replace(slashPos, 1, String("_p_"));

    }

​

    Bool ascending=true;

    if (specVals(specVals.size()-1)<specVals(0))

        ascending=false;

​

    Int startIndex, endIndex;

    if (ascending){

        if (endVal < specVals(0)){

            msg = String("Start value: ") + String::toString(endVal) + String(" is smaller than all spectral values!");