#include <casa/OS/Timer.h>
#include <imageanalysis/ImageAnalysis/ImagePolarimetry.h>
#include <casa/Arrays/Array.h>
#include <casa/Arrays/ArrayMath.h>
#include <casa/Arrays/Vector.h>
#include <casa/Arrays/Matrix.h>
#include <casa/Arrays/MaskedArray.h>
#include <casa/Arrays/MaskArrMath.h>
#include <coordinates/Coordinates/CoordinateSystem.h>
#include <coordinates/Coordinates/StokesCoordinate.h>
#include <coordinates/Coordinates/LinearCoordinate.h>
#include <casa/Exceptions/Error.h>
#include <scimath/Functionals/Polynomial.h>
#include <images/Images/ImageInterface.h>
#include <images/Images/SubImage.h>
#include <images/Images/ImageExpr.h>
#include <imageanalysis/ImageAnalysis/ImageFFT.h>
#include <images/Regions/ImageRegion.h>
#include <images/Images/ImageSummary.h>
#include <images/Images/TempImage.h>
#include <lattices/Lattices/Lattice.h>
#include <lattices/LRegions/LCSlicer.h>
#include <lattices/LEL/LatticeExprNode.h>
#include <lattices/LEL/LatticeExpr.h>
#include <lattices/Lattices/TiledLineStepper.h>
#include <lattices/Lattices/LatticeStepper.h>
#include <lattices/Lattices/LatticeIterator.h>
#include <lattices/Lattices/MaskedLatticeIterator.h>
#include <lattices/LatticeMath/LatticeStatistics.h>
#include <lattices/LRegions/LCPagedMask.h>
#include <casa/Logging/LogIO.h>
#include <casa/Logging/LogOrigin.h>
#include <casa/BasicMath/Math.h>
#include <casa/BasicSL/Constants.h>
#include <scimath/Mathematics/NumericTraits.h>
#include <casa/System/ProgressMeter.h>
#include <casa/Quanta/QC.h>
#include <casa/Quanta/MVAngle.h>
#include <casa/Utilities/GenSort.h>
#include <casa/Utilities/Assert.h>
#include <casa/BasicSL/String.h>
#include <casa/sstream.h>
using namespace casacore;
namespace casa {
ImagePolarimetry::ImagePolarimetry (const ImageInterface<Float>& image)
: _image(image.cloneII())
{
_stokes.resize(4);
_stokesStats.resize(4);
_stokes.set(0);
_stokesStats.set(0);
_findStokes();
_createBeamsEqMat();
}
ImagePolarimetry::ImagePolarimetry(const ImagePolarimetry &other) {
operator=(other);
}