casa6

//# tSubImage.cc: Test program for class SubImage

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

​

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

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

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

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

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

​

#include <imageanalysis/ImageAnalysis/PixelValueManipulator.h>

#include <imageanalysis/ImageAnalysis/SubImageFactory.h>

​

using namespace std;

​

using namespace casacore;

namespace casa {

​

​

std::shared_ptr<ComponentListImage> ImageFactory::createComponentListImage(

    const String& outfile, const Record& cl, const Vector<Int>& shape,

    const Record& csys, Bool overwrite, Bool log, Bool cache

) {

    _checkOutfile(outfile, overwrite);

    ComponentList mycl;

    String err;

    ThrowIf(! mycl.fromRecord(err, cl), err);

    CoordinateSystem mycsys;

    std::unique_ptr<CoordinateSystem> csysPtr;

    if (csys.empty()) {

        mycsys = CoordinateUtil::makeCoordinateSystem(shape, False);

        _centerRefPix(mycsys, shape);

    }

    else {

        csysPtr.reset(

            _makeCoordinateSystem(csys, shape)

        );

        mycsys = *csysPtr;