Listing and manipulating MS metadata (vishead)

NRAO Home > CASA > CASA Cookbook and User Reference Manual

2.2.8 Listing and manipulating MS metadata (vishead)

The vishead task is provided to access keyword information in the Measurement Set. The default inputs are:

#  vishead :: List, get, and put metadata in a measurement set
vis            =          ’’   #  Name of input visibility file
mode           =      ’list’   #  options: list, summary, get, put
   listitems   =          []   #  items to list ([] for all)
async          =       False   #

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 measurement set
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’), {})

[next] [prev] [prev-tail] [front] [up]

More information about CASA may be found at the CASA web page

This code is available under the terms of the GNU General Public Lincense