###"Frequency" definition of image taskname = "clean" vis = "ngc5921_src.split.ms" imagename = "ngc5921_freq" mode = "frequency" alg = "clark" niter = 6000 gain = 0.1 threshold = 8.0 mask = "" cleanbox = [] nchan = 20 start = "1.412849GHz" width = 1 step = "52kHz" imsize = [256, 256] cell = [15.0, 15.0] stokes = "I" field = "0" spw = "" weighting = "briggs" rmode = "norm" robust = 0.5 uvfilter = False uvfilterbmaj = 1.0 uvfilterbmin = 1.0 uvfilterbpa = 0.0 timerange = "" restfreq = "" noise = "0.0Jy" npixels = 0 scales = [0, 3, 10] npercycle = 100 phasecenter = "" async = False #clean(vis="ngc5921_src.split.ms",imagename="ngc5921_freq",mode="frequency",alg="clark",niter=6000,gain=0.1,threshold=8.0,mask="",cleanbox=[],nchan=20,start="1.412849GHz",width=1,step="52kHz",imsize=[256, 256],cell=[15.0, 15.0],stokes="I",field="0",spw="",weighting="briggs",rmode="norm",robust=0.5,uvfilter=False,uvfilterbmaj=1.0,uvfilterbmin=1.0,uvfilterbpa=0.0,timerange="",restfreq="",noise="0.0Jy",npixels=0,scales=[0, 3, 10],npercycle=100,phasecenter="",async=False) ###GALACTIC coordinates taskname = "clean" vis = "ngc5921_src.split.ms" imagename = "ngc5921_gal" mode = "velocity" alg = "hogbom" niter = 6000 gain = 0.1 threshold = 8.0 mask = [''] cleanbox = [] nchan = 38 start = "1595km/s" width = 1 step = "-5.5km/s" imsize = [256, 256] cell = ['15.0arcsec', '15.0arcsec'] stokes = "I" field = "0" spw = "" weighting = "briggs" rmode = "norm" robust = 0.5 uvfilter = False uvfilterbmaj = 1.0 uvfilterbmin = 1.0 uvfilterbpa = 0.0 timerange = "" restfreq = "" noise = "0.0Jy" npixels = 0 scales = [0, 3, 10] npercycle = 100 phasecenter = "GALACTIC 8deg 47.9deg" async = False #clean(vis="ngc5921_src.split.ms",imagename="ngc5921_gal",mode="velocity",alg="hogbom",niter=6000,gain=0.1,threshold=8.0,mask=[''],cleanbox=[],nchan=38,start="1595km/s",width=1,step="-5.5km/s",imsize=[256, 256],cell=['15.0arcsec', '15.0arcsec'],stokes="I",field="0",spw="",weighting="briggs",rmode="norm",robust=0.5,uvfilter=False,uvfilterbmaj=1.0,uvfilterbmin=1.0,uvfilterbpa=0.0,timerange="",restfreq="",noise="0.0Jy",npixels=0,scales=[0, 3, 10],npercycle=100,phasecenter="GALACTIC 8deg 47.9deg",async=False) #####1950 taskname = "clean" vis = "ngc5921_src.split.ms" imagename = "ngc5921_B1950" mode = "velocity" alg = "hogbom" niter = 6000 gain = 0.1 threshold = 8.0 mask = [''] cleanbox = [] nchan = 38 start = "1595km/s" width = 1 step = "-5.5km/s" imsize = [256, 256] cell = ['15.0arcsec', '15.0arcsec'] stokes = "I" field = "0" spw = "" weighting = "briggs" rmode = "norm" robust = 0.5 uvfilter = False uvfilterbmaj = 1.0 uvfilterbmin = 1.0 uvfilterbpa = 0.0 timerange = "" restfreq = "" noise = "0.0Jy" npixels = 0 scales = [0, 3, 10] npercycle = 100 phasecenter = "B1950 15h19m29 5d15m9.73" async = False #clean(vis="ngc5921_src.split.ms",imagename="ngc5921_B1950",mode="velocity",alg="hogbom",niter=6000,gain=0.1,threshold=8.0,mask=[''],cleanbox=[],nchan=38,start="1595km/s",width=1,step="-5.5km/s",imsize=[256, 256],cell=['15.0arcsec', '15.0arcsec'],stokes="I",field="0",spw="",weighting="briggs",rmode="norm",robust=0.5,uvfilter=False,uvfilterbmaj=1.0,uvfilterbmin=1.0,uvfilterbpa=0.0,timerange="",restfreq="",noise="0.0Jy",npixels=0,scales=[0, 3, 10],npercycle=100,phasecenter="B1950 15h19m29 5d15m9.73",async=False) #### try clark mosaic taskname = "mosaic" vis = "orion.ms" imagename = "clark_mosaic" mode = "mfs" alg = "clark" imsize = [300, 300] cell = ['2arcsec', '2arcsec'] phasecenter = "6" stokes = "I" niter = 7000 gain = 0.1 threshold = 0.0 mask = [''] cleanbox = [] nchan = -1 start = 0 width = 1 step = 1 field = "2~10" spw = "0,1" timerange = "" restfreq = "" sdimage = "" modelimage = "" weighting = "natural" mosweight = False rmode = "none" robust = 0.0 ftmachine = "mosaic" cyclefactor = 1.5 cyclespeedup = -1 scaletype = "NONE" minpb = 0.1 sigma = "0.01Jy" targetflux = "1.0Jy" constrainflux = False prior = [''] negcomponent = -1 scales = [0, 3, 10] npercycle = 100 async = False #mosaic(vis="orion.ms",imagename="clark_mosaic",mode="mfs",alg="clark",imsize=[300, 300],cell=['2arcsec', '2arcsec'],phasecenter="6",stokes="I",niter=6000,gain=0.1,threshold=0.0,mask=[''],cleanbox=[],nchan=-1,start=0,width=1,step=1,field="2~10",spw="0,1",timerange="",restfreq="",sdimage="",modelimage="",weighting="natural",mosweight=False,rmode="none",robust=0.0,ftmachine="mosaic",cyclefactor=1.5,cyclespeedup=-1,scaletype="NONE",minpb=0.1,sigma="0.01Jy",targetflux="1.0Jy",constrainflux=False,prior=[''],negcomponent=-1,scales=[0, 3, 10],npercycle=100,async=False) ################################# ##muliscale ################################## taskname = "mosaic" vis = "orion.ms" imagename = "mscale_mosaic" mode = "mfs" alg = "multiscale" imsize = [300, 300] cell = ['2arcsec', '2arcsec'] phasecenter = "6" stokes = "I" niter = 500 gain = 0.3 threshold = 0.0 mask = [''] cleanbox = [] nchan = -1 start = 0 width = 1 step = 1 field = "2~10" spw = "0,1" timerange = "" restfreq = "" sdimage = "" modelimage = "" weighting = "natural" mosweight = False rmode = "none" robust = 0.0 ftmachine = "mosaic" cyclefactor = 1.5 cyclespeedup = -1 scaletype = "NONE" minpb = 0.1 sigma = "0.01Jy" targetflux = "1.0Jy" constrainflux = False prior = [''] negcomponent = -1 scales = [0, 5, 30] npercycle = 100 async = False #mosaic(vis="orion.ms",imagename="mscale_mosaic",mode="mfs",alg="multiscale",imsize=[300, 300],cell=['2arcsec', '2arcsec'],phasecenter="6",stokes="I",niter=500,gain=0.3,threshold=0.0,mask=[''],cleanbox=[],nchan=-1,start=0,width=1,step=1,field="2~10",spw="0,1",timerange="",restfreq="",sdimage="",modelimage="",weighting="natural",mosweight=False,rmode="none",robust=0.0,ftmachine="mosaic",cyclefactor=1.5,cyclespeedup=-1,scaletype="NONE",minpb=0.1,sigma="0.01Jy",targetflux="1.0Jy",constrainflux=False,prior=[''],negcomponent=-1,scales=[0, 5, 30],npercycle=100,async=False) ############################################ ###playing with different scales ########################################### taskname = "mosaic" vis = "orion.ms" imagename = "mscale_mosaic3" mode = "mfs" alg = "multiscale" imsize = [300, 300] cell = ['2arcsec', '2arcsec'] phasecenter = "6" stokes = "I" niter = 2500 gain = 0.3 threshold = 0.0 mask = "" cleanbox = [] nchan = -1 start = 0 width = 1 step = 1 field = "2~10" spw = "0,1" timerange = "" restfreq = "" sdimage = "" modelimage = "" weighting = "natural" mosweight = False rmode = "none" robust = 0.0 ftmachine = "mosaic" cyclefactor = 3 cyclespeedup = -1 scaletype = "NONE" minpb = 0.1 sigma = "0.01Jy" targetflux = "1.0Jy" constrainflux = False prior = [''] negcomponent = -1 scales = [0, 5, 10, 20, 30] npercycle = 100 async = False #mosaic(vis="orion.ms",imagename="mscale_mosaic3",mode="mfs",alg="multiscale",imsize=[300, 300],cell=['2arcsec', '2arcsec'],phasecenter="6",stokes="I",niter=2500,gain=0.3,threshold=0.0,mask="orion.mask6",cleanbox=[],nchan=-1,start=0,width=1,step=1,field="2~10",spw="0,1",timerange="",restfreq="",sdimage="",modelimage="",weighting="natural",mosweight=False,rmode="none",robust=0.0,ftmachine="mosaic",cyclefactor=3,cyclespeedup=-1,scaletype="NONE",minpb=0.1,sigma="0.01Jy",targetflux="1.0Jy",constrainflux=False,prior=[''],negcomponent=-1,scales=[0, 3, 7, 10, 20, 30],npercycle=100,async=False) #Mem ...note prior is SD deconvolved ###lots of trial and error to get to run right #### will provide better feedback and control in beta 1.0 most probably taskname = "mosaic" vis = "orion.ms" imagename = "orion_mem" mode = "mfs" alg = "entropy" imsize = [300, 300] cell = ['2arcsec', '2arcsec'] phasecenter = "6" stokes = "I" niter = 50 gain = 0.1 threshold = 0.0 mask = "" cleanbox = [] nchan = -1 start = 0 width = 1 step = 1 field = "2~10" spw = "0,1" timerange = "" restfreq = "" sdimage = "orion.gbt.im" modelimage = "" weighting = "briggs" mosweight = True rmode = "norm" robust = -1 ftmachine = "mosaic" cyclefactor = 4 cyclespeedup = -1 scaletype = "SAULT" minpb = 0.01 sigma = "4mJy" targetflux = "500Jy" constrainflux = False prior = "orion.sd.model" negcomponent = -1 scales = [0, 3, 10] npercycle = 100 async = False #mosaic(vis="orion.ms",imagename="orion_mem",mode="mfs",alg="entropy",imsize=[300, 300],cell=['2arcsec', '2arcsec'],phasecenter="6",stokes="I",niter=50,gain=0.1,threshold=0.0,mask="",cleanbox=[],nchan=-1,start=0,width=1,step=1,field="2~10",spw="0,1",timerange="",restfreq="",sdimage="orion.gbt.im",modelimage="",weighting="briggs",mosweight=True,rmode="norm",robust=-1,ftmachine="mosaic",cyclefactor=4,cyclespeedup=-1,scaletype="SAULT",minpb=0.01,sigma="4mJy",targetflux="500Jy",constrainflux=False,prior="orion.sd.model",negcomponent=-1,scales=[0, 3, 10],npercycle=100,async=False) ### feather multiscale with taskname = "feather" imagename = "testfeather.im" highres = "mscale_mosaic.image" lowres = "orion.gbt.im" #feather(imagename="testfeather.im",highres="mscale_mosaic.image",lowres="orion.gbt.im") ###check out total flux of the two images ####regrid single dish image to what the interferometric image is taskname = "regridimage" imagename = "orion.gbt.im" template = "mscale_mosaic.image" output = "start_clean" async = False #regridimage(imagename="orion.gbt.im",template="mscale_mosaic.image",output="start_clean",async=False) ###deconvolve regridded image to get a starting taskname = "deconvolve" imagename = "start_clean" model = "orion_sd_start" psf = ['98arcsec', '98arcsec', '0deg'] alg = "multiscale" niter = 700 gain = 0.3 threshold = "0.1Jy" mask = "" scales = [10, 30, 100] sigma = "1.0mJy" targetflux = "1.0Jy" prior = "" async = False #deconvolve(imagename="start_clean",model="orion_sd_start",psf=['98arcsec', '98arcsec', '0deg'],alg="multiscale",niter=700,gain=0.3,threshold="0.1Jy",mask="",scales=[10, 30, 100],sigma="1.0mJy",targetflux="1.0Jy",prior="",async=False) ####Now let's start using this SD model as initial model for msclean taskname = "mosaic" vis = "orion.ms" imagename = "mscale_sd_mosaic" mode = "mfs" alg = "multiscale" imsize = [300, 300] cell = ['2arcsec', '2arcsec'] phasecenter = "6" stokes = "I" niter = 2500 gain = 0.3 threshold = 0.0 mask = "" cleanbox = [] nchan = -1 start = 0 width = 1 step = 1 field = "2~10" spw = "0,1" timerange = "" restfreq = "" sdimage = "" modelimage = "orion_sd_start" weighting = "natural" mosweight = False rmode = "none" robust = 0.0 ftmachine = "mosaic" cyclefactor = 3 cyclespeedup = -1 scaletype = "NONE" minpb = 0.1 sigma = "0.01Jy" targetflux = "1.0Jy" constrainflux = False prior = [''] negcomponent = -1 scales = [0, 5, 10, 20, 30] npercycle = 100 async = False #mosaic(vis="orion.ms",imagename="mscale_sd_mosaic",mode="mfs",alg="multiscale",imsize=[300, 300],cell=['2arcsec', '2arcsec'],phasecenter="6",stokes="I",niter=2500,gain=0.3,threshold=0.0,mask="",cleanbox=[],nchan=-1,start=0,width=1,step=1,field="2~10",spw="0,1",timerange="",restfreq="",sdimage="",modelimage="orion_sd_start",weighting="natural",mosweight=False,rmode="none",robust=0.0,ftmachine="mosaic",cyclefactor=3,cyclespeedup=-1,scaletype="NONE",minpb=0.1,sigma="0.01Jy",targetflux="1.0Jy",constrainflux=False,prior=[''],negcomponent=-1,scales=[0, 5, 10, 20, 30],npercycle=100,async=False) ###Now try with sdimage as start and see the difference...