wproject_regression.py

Go to the documentation of this file.

import sys
import os
import string
import time
from locatescript import copydata
from locatescript import locatescript
import inspect

a=inspect.stack()
stacklevel=0
for k in range(len(a)):
    if (string.find(a[k][1], 'ipython console') > 0):
        stacklevel=k
        break
gl=sys._getframe(stacklevel).f_globals

def data():
    ### return the data files that is needed by the regression script
    return ['coma.ms']

def run(fetch=False):

    #####fetch data
    if fetch:
        for f in data( ):
            copydata( f, os.getcwd( ) )
    
    im=gl['casac'].imager()
    me=gl['me']
    time1=time.time()
    im.open('coma.ms')
    
    im.selectvis(spw=0, field=0);
    mydir=me.direction('J2000', '12h30m48', '12d24m0')
    im.defineimage(nx=200, ny=200, cellx='30arcsec', celly='30arcsec',  phasecenter=mydir); 
    im.make('outlier1re');
    im.defineimage(nx=2000, ny=2000, cellx='30arcsec',celly='30arcsec',phasecenter=0, facets=1);
    im.setoptions(ftmachine='wproject',wprojplanes=400, padding=1.2)
    im.make('bla1re')
    
    im.clean(algorithm='wfclark',model=['bla1re', 'outlier1re'], image=['bla1re.restored', 'outlier1re.restored'], niter=10000);

    im.done()
    time2=time.time();
    print 'Time taken for wproject= ', (time2-time1)/60,'mins'
    #resturn the images that will be templated and compared in future runs
    return ['bla1re.restored', 'outlier1re.restored']