SDAlgorithmBase.h

Go to the documentation of this file.
00001 //# SDAlgorithmBase.h: Definition for SDAlgorithmBase
00002 //# Copyright (C) 1996,1997,1998,1999,2000,2002
00003 //# Associated Universities, Inc. Washington DC, USA.
00004 //#
00005 //# This library is free software; you can redistribute it and/or modify it
00006 //# under the terms of the GNU Library General Public License as published by
00007 //# the Free Software Foundation; either version 2 of the License, or (at your
00008 //# option) any later version.
00009 //#
00010 //# This library is distributed in the hope that it will be useful, but WITHOUT
00011 //# ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
00012 //# FITNESS FOR A PARTICULAR PURPOSE.  See the GNU Library General Public
00013 //# License for more details.
00014 //#
00015 //# You should have received a copy of the GNU Library General Public License
00016 //# along with this library; if not, write to the Free Software Foundation,
00017 //# Inc., 675 Massachusetts Ave, Cambridge, MA 02139, USA.
00018 //#
00019 //# Correspondence concerning AIPS++ should be adressed as follows:
00020 //#        Internet email: aips2-request@nrao.edu.
00021 //#        Postal address: AIPS++ Project Office
00022 //#                        National Radio Astronomy Observatory
00023 //#                        520 Edgemont Road
00024 //#                        Charlottesville, VA 22903-2475 USA
00025 //#
00026 //#
00027 //# $Id$
00028 
00029 #ifndef SYNTHESIS_SDALGORITHMBASE_H
00030 #define SYNTHESIS_SDALGORITHMBASE_H
00031 
00032 #include <ms/MeasurementSets/MeasurementSet.h>
00033 //#include <synthesis/MeasurementComponents/SkyModel.h>
00034 #include <casa/Arrays/Matrix.h>
00035 #include <images/Images/ImageInterface.h>
00036 #include <images/Images/PagedImage.h>
00037 #include <images/Images/TempImage.h>
00038 #include <casa/Logging/LogMessage.h>
00039 #include <casa/Logging/LogSink.h>
00040 #include <casa/System/PGPlotter.h>
00041 
00042 #include <casa/aips.h>
00043 #include <images/Images/ImageInterface.h>
00044 #include <components/ComponentModels/ComponentList.h>
00045 #include <casa/BasicSL/String.h>
00046 #include <synthesis/TransformMachines/StokesImageUtil.h>
00047 
00048 #include<synthesis/ImagerObjects/SDMaskHandler.h>
00049 #include<synthesis/ImagerObjects/SIImageStore.h>
00050 #include<synthesis/ImagerObjects/SIImageStoreMultiTerm.h>
00051 
00052 namespace casa { //# NAMESPACE CASA - BEGIN
00053 
00054   /* Forware Declaration */
00055   class SIMinorCycleController;
00056 
00057 
00058 class SDAlgorithmBase {
00059 public:
00060 
00061   // Empty constructor
00062   SDAlgorithmBase();
00063  virtual  ~SDAlgorithmBase();
00064 
00065   // Non virtual. Wrapper function implemented only in the base class. 
00066   void deconvolve( SIMinorCycleController& loopController,  
00067                    SHARED_PTR<SIImageStore> &imagestore,
00068                    Int deconvolverid);
00069 
00070   void setRestoringBeam( GaussianBeam restbeam, String usebeam );
00071   //  void setMaskOptions( String maskstring );
00072 
00073   // Base Class contains standard restoration. Overload for more complex behaviour.
00074   virtual void restore( SHARED_PTR<SIImageStore> imagestore );
00075   virtual void pbcor( SHARED_PTR<SIImageStore> imagestore );
00076 
00077   virtual String getAlgorithmName(){return itsAlgorithmName;};
00078 
00079   virtual uInt getNTaylorTerms(){return 1;};
00080 
00081 protected:
00082 
00083   // Pure virtual functions to be implemented by various algorithm deconvolvers.
00084   virtual void takeOneStep( Float loopgain, Int cycleNiter, Float cycleThreshold, 
00085                             Float &peakresidual, Float &modelflux, Int& iterdone )=0;
00086   //  virtual void initializeDeconvolver( Float &peakresidual, Float &modelflux )=0;
00087   virtual void initializeDeconvolver()=0;
00088   virtual void finalizeDeconvolver()=0;
00089 
00090   // Base Class implements the option of single-plane images for the minor cycle.
00091   virtual void queryDesiredShape(Int &nchanchunks, Int& npolchunks, IPosition imshape);
00092 
00093 
00094   // Non virtual. Implemented only in the base class.
00095   Int checkStop( SIMinorCycleController &loopcontrols, Float currentresidual );
00096   Bool findMaxAbs(const Array<Float>& lattice,Float& maxAbs,IPosition& posMaxAbs);
00097   Bool findMaxAbsMask(const Array<Float>& lattice,const Array<Float>& mask,
00098                       Float& maxAbs,IPosition& posMaxAbs);
00099 
00100   // Algorithm name
00101   String itsAlgorithmName;
00102 
00103   SHARED_PTR<SIImageStore> itsImages; //sOriginalImages;
00104 
00105   //    Vector<Slicer> itsDecSlices;
00106   //   SubImage<Float> itsResidual, itsPsf, itsModel, itsImage;
00107   
00108   IPosition itsMaxPos;
00109   Float itsPeakResidual;
00110   Float itsModelFlux;
00111 
00112   //  SDMaskHandler itsMaskHandler;
00113   //Array<Float> itsMatMask;
00114 
00115   GaussianBeam itsRestoringBeam;
00116   String itsUseBeam;
00117   //  String itsMaskString;
00118   //  Bool itsIsMaskLoaded; // Annoying state variable. Remove if possible. 
00119 
00120 };
00121 
00122 } //# NAMESPACE CASA - END
00123 
00124 #endif
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Defines

Generated on 31 Aug 2016 for casa by  doxygen 1.6.1