def subcont(inp, out, cont_chan=0):
"""
Subtract continuum from line maps
"""
# Clear temporary space for the continuum channel
cont = "/tmp/cont"
clear_path(cont)
# Do the subtraction
MirTask("imsub")(inp=inp, out=cont, region='image(%d)'%(cont_chan+1))
MirTask("maths")(exp='<%s>-<%s>'%(inp, cont), out=out, options='grow')
return
def subcont(inp, out, cont_chan=0):
"""
Subtract continuum from line maps
"""
# Clear temporary space for the continuum channel
cont = "/tmp/cont"
clear_path(cont)
# Do the subtraction
MirTask("imsub")(inp=inp, out=cont, region='image(%d)' % (cont_chan + 1))
MirTask("maths")(exp='<%s>-<%s>' % (inp, cont), out=out, options='grow')
return
def subcont_convol(inp, out, fwhm): """ Subtract continuum and do convolution in one call """ # Clear temporary space for the continuum-subtracted map sub = "/tmp/sub" clear_path(sub) # Subtract continuum subcont(inp, sub) # Convolve with beam convol(sub, out, fwhm) return
def subcont_convol(inp, out, fwhm):
"""
Subtract continuum and do convolution in one call
"""
# Clear temporary space for the continuum-subtracted map
sub = "/tmp/sub"
clear_path(sub)
# Subtract continuum
subcont(inp, sub)
# Convolve with beam
convol(sub, out, fwhm)
return