casa6

// -*- C++ -*-

//# Utils.cc: Implementation of global functions from Utils.h 

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

//# 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 <msvis/MSVis/VisBuffer.h>

#include <casacore/casa/Logging/LogIO.h>

#include <casacore/ms/MeasurementSets/MSColumns.h>

#include <casacore/measures/Measures/MEpoch.h>

#include <casacore/measures/Measures/MeasTable.h>

#include <synthesis/TransformMachines/Utils.h>

#include <synthesis/TransformMachines/StokesImageUtil.h>

#include <casacore/casa/Utilities/Assert.h>

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

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

#include <casacore/lattices/LEL/LatticeExpr.h>

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

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

#include <casacore/casa/Containers/Record.h>

#include <casacore/lattices/Lattices/LatticeIterator.h>

#include <casacore/lattices/Lattices/TiledLineStepper.h> 

#include <casacore/lattices/Lattices/LatticeStepper.h> 

#include <casacore/lattices/LatticeMath/LatticeFFT.h>

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

​

using namespace casacore;

namespace casa{

  //

  //--------------------------------------------------------------------------------------------

  //  

  void storeImg(String fileName,ImageInterface<Complex>& theImg, Bool writeReIm)

  {

    PagedImage<Complex> ctmp(theImg.shape(), theImg.coordinates(), fileName);

    LatticeExpr<Complex> le(theImg);

    ctmp.copyData(le);

    if (writeReIm)

      {

    ostringstream reName,imName;

    reName << "re" << fileName;

    imName << "im" << fileName;

    {