casa6

import os

import numpy

import pylab as pl

​

from casatasks.private.casa_transition import is_CASA6

if is_CASA6:

    from casatasks import casalog

    from casatools import quanta

    from casatools import image

else:

    from taskinit import casalog

    from taskinit import qatool as quanta

    from taskinit import iatool as image

​

​

qa = quanta()

​

​

def plotprofilemap(imagename=None, figfile=None, overwrite=None, transparent=None,

                   pol=None, spectralaxis=None, restfreq=None, plotrange=None, title=None,

                   linecolor=None, linestyle=None, linewidth=None,

                   separatepanel=None, plotmasked=None, maskedcolor=None,

                   showaxislabel=None, showtick=None, showticklabel=None,

                   figsize=None, numpanels=None):

    casalog.origin('plotprofilemap')

​

    if len(figfile) > 0 and os.path.exists(figfile) and overwrite == False:

        raise RuntimeError('overwrite is False and output file exists: \'%s\''%(figfile))

​

    image = SpectralImage(imagename)

    npol = len(image.stokes)

    if pol < 0 or pol > npol - 1:

        raise RuntimeError('pol {pol} is out of range (Stokes axis {stokes})'.format(pol=pol,stokes=image.stokes))

​

    parsed_size = parse_figsize(figsize)

    nx, ny = parse_numpanels(numpanels)

    if (not isinstance(restfreq, str)) or len(restfreq) == 0:

        restfreq = None

    plot_profile_map(image, figfile, pol, spectralaxis, restfreq, title, linecolor, linestyle, linewidth,

                     separatepanel, plotmasked, maskedcolor,

                     showaxislabel, showtick, showticklabel, parsed_size,

                     nx, ny, transparent, plotrange)

​

​

NoData = -32767.0

NoDataThreshold = NoData + 10000.0

LightSpeedQuantity = qa.constants('c')

LightSpeed = qa.convert(LightSpeedQuantity, 'km/s')['value'] # speed of light in km/s

DPIDetail = 130

​

dsyb = '$^\circ$'

hsyb = ':'

msyb = ':'

​

​

def Deg2HMS(x, arrowance):

    # Transform degree to HHMMSS.sss format

    xx = x % 360 + arrowance

    h = int(xx / 15)

    m = int((xx % 15) * 4)

    s = ((xx % 15) * 4 - m) * 60.0

    return (h, m, s)

​