MS Metadata List/Change
The vishead task is provided to access keyword information in the MeasurementSet. The default inputs are:
# vishead :: List, get, and put metadata in a MeasurementSet
vis = '' # Name of input visibility file
mode = 'list' # options: list, summary, get, put
listitems = [] # items to list ([] for all)
The mode = ’summary’ option just gives the same output as listobs.
For mode = ’list’ the options are: ’telescope’, ’observer’, ’project’, ’field’, ’freq_group_name’, ’spw_name’, ’schedule’, ’schedule_type’, ’release_date’.
CASA <29>: vishead('ngc5921.demo.ms',mode='list',listitems=[])
Out[29]:
{'cal_grp': (array([-1, -1, -1], dtype=int32), {}),
'field': (array(['1331+30500002_0', '1445+09900002_0', 'N5921_2'],
dtype='|S16'),
{}),
'fld_code': (array(['C', 'A', ''],
dtype='|S2'), {}),
'freq_group_name': (array(['none'],
dtype='|S5'), {}),
'log': ({'r1': False}, {}),
'observer': (array(['TEST'],
dtype='|S5'), {}),
'project': (array([''],
dtype='|S1'), {}),
'ptcs': ({'r1': array([[[-2.74392758]],
[[ 0.53248521]]]),
'r2': array([[[-2.42044692]],
[[ 0.17412604]]]),
'r3': array([[[-2.26020138]],
[[ 0.08843002]]])},
{'MEASINFO': {'Ref': 'J2000', 'type': 'direction'},
'QuantumUnits': array(['rad', 'rad'],
dtype='|S4')}),
'release_date': (array([ 4.30444800e+09]),
{'MEASINFO': {'Ref': 'TAI', 'type': 'epoch'},
'QuantumUnits': array(['s'],
dtype='|S2')}),
'schedule': ({'r1': False}, {}),
'schedule_type': (array([''],
dtype='|S1'), {}),
'source_name': (array(['1331+30500002_0', '1445+09900002_0', 'N5921_2'],
dtype='|S16'),
{}),
'spw_name': (array(['none'],
dtype='|S5'), {}),
'telescope': (array(['VLA'],
dtype='|S4'), {})}
You can use mode=’get’ to retrieve the values of specific keywords, and likewise mode=’put’ to change them. The inputs are:
mode = 'get' # options: list, summary, get, put
hdkey = '' # keyword to get/put
hdindex = '' # keyword index to get/put, counting from zero. ==>all
and
# vishead :: List, summary, get, and put metadata in a MeasurementSet
mode = 'put' # options: list, summary, get, put
hdkey = '' # keyword to get/put
hdindex = '' # keyword index to get/put, counting from zero. ==>all
hdvalue = '' # value of hdkey
For example, a common operation is to change the Telescope name (e.g. if it is unrecognized), e.g.
CASA <36>: vishead('ngc5921.demo.ms',mode='get',hdkey='telescope')
Out[36]:
(array(['VLA'],
dtype='|S4'), {})
CASA <37>: vishead('ngc5921.demo.ms',mode='put',hdkey='telescope',hdvalue='JVLA')
CASA <38>: vishead('ngc5921.demo.ms',mode='get',hdkey='telescope')
Out[38]:
(array(['JVLA'],
dtype='|S5'), {})