#
#
include 'viewer.g'
include 'interactivemask.g'
include 'imager.g'
include 'image.g'
include 'ms.g'
#
#
# MEM
#
imgr:=imager('orion.ms');
imgr.setdata(fieldid=[3:11],spwid=[1:2],mode="none");
imgr.setimage(nx=300,ny=300,cellx='2.0arcsec',celly='2.0arcsec',stokes='I',
	mode="mfs",fieldid=[7],spwid=[1:2]);
imgr.weight(type='briggs',robust=-0.5,mosaic=T,rmode='norm');
imgr.setvp(dovp=T,usedefaultvp=T,dosquint=F);
imgr.setmfcontrol(cyclefactor=3.0);
#
# first try
#
imgr.mem(algorithm='mfentropy',niter=8,sigma='0.01Jy',
	targetflux='100Jy', displayprogress=T,fixed=F,
	complist='',prior='',
        model='orion-vlamem_MZ1.model',
	image='orion-vlamem_MZ1.image',
        residual='orion-vlamem_MZ1.resid');

#
# dv.gui() should be off before doing this:
# 
mymask:=interactivemask(refimage='orion-vlamem_MZ1.image',
	mask='orion-vlamem_MZ1.mask');
mymask.start();
#
# second round:
#
# Now know how to make masks, use the provided mask now:
#
imgr.mem(algorithm='mfentropy',niter=8,sigma='0.01Jy',
	targetflux='100Jy', displayprogress=T,fixed=F,
	complist='',prior='',
        model='orion-vlamem_MZ2.model',
	image='orion-vlamem_MZ2.image',
        residual='orion-vlamem_MZ2.resid',
	mask='orion.mask6');

# third, fourth,.. round:
#
#
imgr.mem(algorithm='mfentropy',niter=35,sigma='0.0075Jy',
	targetflux='110Jy',displayprogress=T,fixed=F,
	complist='',prior='',
        model='orion-vlamem_MZ5.model',
	image='orion-vlamem_MZ5.image',
        residual='orion-vlamem_MZ5.resid',
	mask='orion.mask6');
# 
#
#
# Feathering
#
imgr.setvp(dovp=T,usedefaultvp=T,dosquint=F);
imgr.feather(image='orion.feather_MZ1.image',
	highres='orion-vlamem_MZ5.image', 
	lowres='orion.gbt2.im');
#
imgr.setvp(dovp=T,usedefaultvp=T,dosquint=F);
imgr.feather(image='orion.feather_MZtest.image',
	highres='orion_vlamem.im', 
	lowres='orion.gbt2.im');
#
#
#
##############################################################
#
# The MS Clean method
#
##############################################################
#
#
imgr:=imager('orion.ms');
imgr.setdata(fieldid=[3:11],spwid=[1:2],mode="none");
imgr.setimage(nx=300,ny=300,cellx='2.0arcsec',celly='2.0arcsec',stokes='I',
	mode="mfs",fieldid=[7],spwid=[1:2]);
imgr.weight(type='briggs',robust=-1,mosaic=T,rmode='norm');
imgr.setvp(dovp=T,usedefaultvp=T,dosquint=F);
imgr.setmfcontrol(cyclefactor=3.0,stoplargenegatives=-1,
   	cyclespeedup=500);
#
#
imgr.setscales(scalemethod='uservector',uservector=[0,2,5,15,25]);
#
imgr.makemodelfromsd(sdimage='orion.gbt2.im',
	modelimage='orion-gbt2_MZ1.model',maskimage='');
# 
imgr.clean(algorithm='mfmultiscale',gain=0.15,niter=1000,threshold='0.01Jy',
 	model='orion-gbt2_MZ1.model',
	mask='orion-mask6', 
	image='orion-gbt2-DI_MZ1.image',
	residual='orion-gbt2_MZ1.resid',
	interactive=T,npercycle=500);

#
# more iterations
#
imgr.clean(algorithm='mfmultiscale',gain=0.10,niter=20000,threshold='0.005Jy',
 	model='orion-gbt2_MZ1.model',
	mask='orion-mask6', 
	image='orion-gbt2-DI_MZ3.image',
	residual='orion-gbt2_MZ3.resid',
	interactive=T,npercycle=500);
#
#
#
#