"""We use the standard kwargs mechanism for arguments. For example::
kernelname = get_parameter(kwargs, "kernel", "2d")
oversampling = get_parameter(kwargs, "oversampling", 8)
padding = get_parameter(kwargs, "padding", 2)
The kwargs may need to be passed down to called functions.
All functions possess an API which is always of the form::
def processing_function(idatastruct1, idatastruct2, ..., *kwargs):
return odatastruct1, odatastruct2,... other
Processing parameters are passed via the standard Python kwargs approach.
Inside a function, the values are retrieved can be accessed directly from the
kwargs dictionary, or if a default is needed a function can be used::
log = get_parameter(kwargs, 'log', None)
vis = get_parameter(kwargs, 'visibility', None)
Function parameters should obey a consistent naming convention:
======= =======
Name Meaning
======= =======
vis Name of Visibility
sc Name of SkyComponent
gt Name of GainTable
conf Name of Configuration
im Name of input image
qa Name of quality assessment
log Name of processing log
======= =======
If a function argument has a better, more descriptive name
e.g. normalised_gt, newphasecentre, use it.
Keyword=value pairs should have descriptive names. The names should
be lower case with underscores to separate words:
==================== ================================== ===========================
Name Meaning Example
==================== ================================== ===========================
loop_gain Clean loop gain 0.1
niter Number of iterations 10000
eps Fractional tolerance 1e-6
threshold Absolute threshold 0.001
fractional_threshold Threshold as fraction of e.g. peak 0.1
G_solution_interval Solution interval for G term 100
phaseonly Do phase-only solutions True
phasecentre Phase centre (usually as SkyCoord) SkyCoord("-1.0d", "37.0d",
frame='icrs',
equinox='J2000')
spectral_mode Visibility processing mode 'mfs' or 'channel'
==================== ================================== ===========================
"""
__all__ = ["rascil_path", "rascil_data_path", "get_parameter"]
import logging
import os
log = logging.getLogger("rascil-logger")
[docs]
def rascil_path(path):
"""Converts a path that might be relative to RASCIL root into an
absolute path::
rascil_data_path('models/SKA1_LOW_beam.fits')
'/Users/timcornwell/Code/rascil/data/models/SKA1_LOW_beam.fits'
:param path:
:return: absolute path
"""
project_root = os.path.dirname(
os.path.dirname(os.path.abspath(__file__)) + "/../../"
)
rascilhome = os.getenv("RASCIL", project_root)
return os.path.join(rascilhome, path)
[docs]
def rascil_data_path(path, check=True):
"""Converts a path that might be relative to the RASCIL data directory into an
absolute path::
rascil_data_path('models/SKA1_LOW_beam.fits')
'/Users/timcornwell/Code/rascil/data/models/SKA1_LOW_beam.fits'
The data path default is rascil_path('data') but may be overriden with the
environment variable RASCIL_DATA.
:param check: Check path exists
:param path:
:return: absolute path
"""
rascil_data_home = os.getenv("RASCIL_DATA", None)
if rascil_data_home is None:
dp = rascil_path("data/")
else:
dp = rascil_data_home
if check:
if not os.path.exists(dp):
raise FileNotFoundError(
"RASCIL data directory {} does not exist".format(dp)
)
dp = os.path.join(dp, path)
return dp
[docs]
def get_parameter(kwargs, key, default=None):
"""Get a specified named value for this (calling) function
The parameter is searched for in kwargs
:param kwargs: Parameter dictionary
:param key: Key e.g. 'loop_gain'
:param default: Default value
:return: result
"""
if kwargs is None:
return default
value = default
if key in kwargs.keys():
value = kwargs[key]
return value