casa6

import os

from taskinit import *

import pdb

ia, dc, tb=gentools(['ia', 'dc', 'tb'])

​

def deconvolve(imagename,model,psf,alg,niter,gain,threshold,mask,scales,sigma,targetflux,prior):

    """ deconvolve: Image based deconvolver. The psf provided is

    deconvolved out of the image provided.

​

        Keyword arguments:

        imagename -- Name of input image to be deconvoled

        model     -- Name of output to store found model, a.k.a clean components

        psf       -- Name of psf image to use e.g  psf='mypsf.image' .

                     But if the psf has 3 parameter, then

                 a gaussian psf is assumed with the values representing

                     the major , minor and position angle  values

                     e.g  psf=['3arcsec', '2.5arcsec', '10deg']

        alg       -- algorithm to use clark, hogbom or multiscale or mem. if multiscale 

                     the parameter scale is used to define the number of scales.

        niter     -- Number of iteration

        gain      -- CLEAN gain parameter; fraction to remove from peak (< 1.0)

        threshold -- deconvolution stopping threshold: if no peak above

                     found above this level 

        mask      -- mask image (same shape as image and psf) to limit region

                     where deconvoltion is to occur

        ------parameters useful for multiscale only

        scales     -- parameter needed for multiscale clean. default value [0,3,10]

        ------parameters useful for mem only

        sigma     -- Estimated noise for image

        targetflux -- Target total flux in image 

        prior     -- Prior image to guide mem

​

    """

    casalog.origin('deconvolve')

    tmppsf=''

    tmpImagename=''

    if(len(psf)==0):

        raise Exception, "****give some psf please****"

    if(len(psf)==1):

        if (type(psf[0])==str and os.path.exists(psf[0])):

            psf=psf[0]

            ia.open(psf)

            psfaxis=len(ia.shape())

            ia.close()

            if(psfaxis<4):

                modPsf=_add_axes(psf)

                if modPsf is False:

                    raise Exception, "****problem with input psf image****"

                else:

                    tmppsf='__decon_tmp_psf'

                    ia.fromimage(tmppsf,modPsf, overwrite=True)

                    ia.close()

                    psf=tmppsf

            ia.open(imagename)

            imnaxis=len(ia.shape())

            ia.close()

            if(imnaxis<4):

                tmpImagename=_add_axes(imagename)

                if tmpImagename is False:

                    raise Exception, "****problem with input diry image****"