00001 //# MosaicFTNew.h: Definition for MosaicFTNew 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_TRANSFORM2_MOSAICFTNEW_H 00030 #define SYNTHESIS_TRANSFORM2_MOSAICFTNEW_H 00031 00032 #include <synthesis/TransformMachines2/MosaicFT.h> 00033 00034 namespace casa { //# NAMESPACE CASA - BEGIN 00035 00036 namespace refim{ 00037 /* 00038 An copy MosaicFT except 00039 Looks like it is just to get differentlently normalized images i.e (image*nx*ny) 00040 which implies somewhere some code is just be using FFT in forward (the toFrequency direction) 00041 when it should be the reverse (or vice-versa) 00042 */ 00043 00044 class MosaicFTNew : public MosaicFT { 00045 public: 00046 00047 // Constructor: cachesize is the size of the cache in words 00048 // (e.g. a few million is a good number), tilesize is the 00049 // size of the tile used in gridding (cannot be less than 00050 // 12, 16 works in most cases). 00051 // <group> 00052 MosaicFTNew(SkyJones* sj, MPosition mloc, String stokes, 00053 Long cachesize, Int tilesize=16, 00054 Bool usezero=True, Bool useDoublePrec=False): 00055 MosaicFT(sj,mloc,stokes,cachesize,tilesize,usezero,useDoublePrec){} 00056 00057 refim::FTMachine* cloneFTM(); 00058 00059 00060 // </group> 00061 00062 // Get the final image: do the Fourier transform and 00063 // grid-correct, then optionally normalize by the summed weights 00064 ImageInterface<Complex>& getImage(Matrix<Float>&, Bool normalize=True); 00065 00066 // Get the final weights image 00067 void getWeightImage(ImageInterface<Float>&, Matrix<Float>&); 00068 00069 virtual String name() const {return "MosaicFTNew";}; 00070 00071 protected: 00072 00073 }; 00074 } //refim ends 00075 } //# NAMESPACE CASA - END 00076 00077 #endif