xradio 0.0.27__py3-none-any.whl → 0.0.29__py3-none-any.whl

xradio/vis/_vis_utils/_ms/_tables/create_field_and_source_xds.py

@@ -0,0 +1,710 @@
+import os
+import time
+from typing import Tuple, Union
+
+import numpy as np
+import xarray as xr
+
+from xradio.vis._vis_utils._ms.msv2_to_msv4_meta import (
+    column_description_casacore_to_msv4_measure,
+)
+from xradio.vis._vis_utils._ms.msv4_sub_xdss import interpolate_to_time
+from xradio.vis._vis_utils._ms.subtables import subt_rename_ids
+from xradio.vis._vis_utils._ms._tables.read import (
+    convert_casacore_time_to_mjd,
+    make_taql_where_between_min_max,
+    read_generic_table,
+)
+from xradio.vis._vis_utils._ms._tables.table_query import open_table_ro
+import graphviper.utils.logger as logger
+
+
+def cast_to_str(x):
+    if isinstance(x, list):
+        return x[0]
+    else:
+        return x
+
+
+def create_field_and_source_xds(
+    in_file,
+    field_id,
+    spectral_window_id,
+    field_times,
+    is_single_dish,
+    time_min_max: Tuple[np.float64, np.float64],
+    ephemeris_interp_time: Union[xr.DataArray, None] = None,
+):
+    """
+    Create a field and source xarray dataset (xds) from the given input file, field ID, and spectral window ID.
+
+    Parameters:
+    ----------
+    in_file : str
+        The path to the input file.
+    field_id : int
+        The ID of the field.
+    spectral_window_id : int
+        The ID of the spectral window.
+    time_min_max : Tuple[np.float64, np.float46]
+        Min / max times to constrain loading (usually to the time range relevant to an MSv4)
+    ephemeris_interp_time : Union[xr.DataArray, None]
+        Time axis to interpolate the ephemeris data vars to (usually main MSv4 time)
+
+    Returns:
+    -------
+    field_and_source_xds : xr.Dataset
+        The xarray dataset containing the field and source information.
+    """
+
+    start_time = time.time()
+
+    field_and_source_xds = xr.Dataset()
+
+    field_and_source_xds, ephemeris_path, ephemeris_table_name = (
+        create_field_info_and_check_ephemeris(
+            field_and_source_xds, in_file, field_id, field_times, is_single_dish
+        )
+    )
+    source_id = field_and_source_xds.attrs["source_id"]
+
+    if ephemeris_path is not None:
+        field_and_source_xds = extract_ephemeris_info(
+            field_and_source_xds,
+            ephemeris_path,
+            ephemeris_table_name,
+            is_single_dish,
+            time_min_max,
+            ephemeris_interp_time,
+        )
+        field_and_source_xds.attrs["is_ephemeris"] = True
+        field_and_source_xds = extract_source_info(
+            field_and_source_xds,
+            in_file,
+            True,
+            source_id,
+            spectral_window_id,
+        )
+
+    else:
+        field_and_source_xds = extract_source_info(
+            field_and_source_xds, in_file, False, source_id, spectral_window_id
+        )
+        field_and_source_xds.attrs["is_ephemeris"] = False
+
+    logger.debug(
+        f"create_field_and_source_xds() execution time {time.time() - start_time:0.2f} s"
+    )
+
+    return field_and_source_xds
+
+
+def extract_ephemeris_info(
+    xds,
+    path,
+    table_name,
+    is_single_dish,
+    time_min_max: Tuple[np.float64, np.float64],
+    interp_time: Union[xr.DataArray, None],
+):
+    """
+    Extracts ephemeris information from the given path and table name and adds it to the xarray dataset.
+
+    Parameters:
+    ----------
+    xds : xr.Dataset
+        The xarray dataset to which the ephemeris information will be added.
+    path : str
+        The path to the input file.
+    table_name : str
+        The name of the ephemeris table.
+    time_min_max : Tuple[np.float46, np.float64]
+        Min / max times to constrain loading (usually to the time range relevant to an MSv4)
+    ephemeris_interp_time : Union[xr.DataArray, None]
+        Time axis to interpolate the data vars to (usually main MSv4 time)
+
+    Returns:
+    -------
+    xds : xr.Dataset
+        The xarray dataset with the added ephemeris information.
+    """
+    # The JPL-Horizons ephemris table implmenation in CASA does not follow the standard way of defining measures.
+    # Consequently a lot of hardcoding is needed to extract the information.
+    # https://casadocs.readthedocs.io/en/latest/notebooks/external-data.html
+
+    min_max_mjd = (
+        convert_casacore_time_to_mjd(time_min_max[0]),
+        convert_casacore_time_to_mjd(time_min_max[1]),
+    )
+    taql_time_range = make_taql_where_between_min_max(
+        min_max_mjd, path, table_name, "MJD"
+    )
+    ephemeris_xds = read_generic_table(
+        path, table_name, timecols=["MJD"], taql_where=taql_time_range
+    )
+
+    # print(ephemeris_xds)
+
+    assert len(ephemeris_xds.ephemeris_id) == 1, "Non standard ephemeris table."
+    ephemeris_xds = ephemeris_xds.isel(ephemeris_id=0)
+
+    ephemeris_meta = ephemeris_xds.attrs["other"]["msv2"]["ctds_attrs"]
+    ephemris_column_description = ephemeris_xds.attrs["other"]["msv2"]["ctds_attrs"][
+        "column_descriptions"
+    ]
+
+    assert (
+        ephemeris_meta["obsloc"] == "GEOCENTRIC"
+    ), "Only geocentric observer ephemeris are supported."
+
+    if "posrefsys" in ephemeris_meta:
+        sky_coord_frame = ephemeris_meta["posrefsys"].replace("ICRF/", "")
+    else:
+        sky_coord_frame = "ICRS"  # We will have to just assume this.
+
+    # Find out witch keyword is used for units (UNIT/QuantumUnits)
+    if "UNIT" in ephemris_column_description["RA"]["keywords"]:
+        unit_keyword = "UNIT"
+    else:
+        unit_keyword = "QuantumUnits"
+
+    coords = {
+        "ellipsoid_pos_label": ["lon", "lat", "dist"],
+        "time": ephemeris_xds["time"].data,
+        "sky_pos_label": ["ra", "dec", "dist"],
+    }
+
+    xds["SOURCE_POSITION"] = xr.DataArray(
+        np.column_stack(
+            (
+                ephemeris_xds["ra"].data,
+                ephemeris_xds["dec"].data,
+                ephemeris_xds["rho"].data,
+            )
+        ),
+        dims=["time", "sky_pos_label"],
+    )
+    # Have to use cast_to_str because the ephemris table units are not consistently in a list or a string.
+    sky_coord_units = [
+        cast_to_str(ephemris_column_description["RA"]["keywords"][unit_keyword]),
+        cast_to_str(ephemris_column_description["DEC"]["keywords"][unit_keyword]),
+        cast_to_str(ephemris_column_description["Rho"]["keywords"][unit_keyword]),
+    ]
+    xds["SOURCE_POSITION"].attrs.update(
+        {"type": "sky_coord", "frame": sky_coord_frame, "units": sky_coord_units}
+    )
+
+    xds["SOURCE_RADIAL_VELOCITY"] = xr.DataArray(
+        ephemeris_xds["radvel"].data, dims=["time"]
+    )
+    xds["SOURCE_RADIAL_VELOCITY"].attrs.update(
+        {
+            "type": "quantity",
+            "units": [
+                cast_to_str(
+                    ephemris_column_description["RadVel"]["keywords"][unit_keyword]
+                )
+            ],
+        }
+    )
+
+    observation_position = [
+        ephemeris_meta["GeoLong"],
+        ephemeris_meta["GeoLat"],
+        ephemeris_meta["GeoDist"],
+    ]
+    xds["OBSERVATION_POSITION"] = xr.DataArray(
+        observation_position, dims=["ellipsoid_pos_label"]
+    )
+    xds["OBSERVATION_POSITION"].attrs.update(
+        {
+            "type": "location",
+            "units": ["deg", "deg", "m"],
+            "data": observation_position,
+            "ellipsoid": "WGS84",
+            "origin_object_name": "Earth",
+            "coordinate_system": ephemeris_meta["obsloc"].lower(),
+        }
+    )  # I think the units are ['deg','deg','m'] and 'WGS84'.
+
+    # Add optional data
+    # NORTH_POLE_POSITION_ANGLE
+
+    if "np_ang" in ephemeris_xds.data_vars:
+        xds["NORTH_POLE_POSITION_ANGLE"] = xr.DataArray(
+            ephemeris_xds["np_ang"].data, dims=["time"]
+        )
+        xds["NORTH_POLE_POSITION_ANGLE"].attrs.update(
+            {
+                "type": "quantity",
+                "units": [
+                    cast_to_str(
+                        ephemris_column_description["NP_ang"]["keywords"][unit_keyword]
+                    )
+                ],
+            }
+        )
+
+    if "np_dist" in ephemeris_xds.data_vars:
+        xds["NORTH_POLE_ANGULAR_DISTANCE"] = xr.DataArray(
+            ephemeris_xds["np_dist"].data, dims=["time"]
+        )
+        xds["NORTH_POLE_ANGULAR_DISTANCE"].attrs.update(
+            {
+                "type": "quantity",
+                "units": [
+                    cast_to_str(
+                        ephemris_column_description["NP_dist"]["keywords"][unit_keyword]
+                    )
+                ],
+            }
+        )
+
+    if "disklong" in ephemeris_xds.data_vars:
+        xds["SUB_OBSERVER_POSITION"] = xr.DataArray(
+            np.column_stack(
+                (
+                    ephemeris_xds["disklong"].data,
+                    ephemeris_xds["disklat"].data,
+                    np.zeros(ephemeris_xds["disklong"].shape),
+                )
+            ),
+            dims=["time", "ellipsoid_pos_label"],
+        )
+
+        if "DiskLong" in ephemris_column_description:
+            units_key_lon = "DiskLong"
+            units_key_lat = "DiskLat"
+        else:
+            units_key_lon = "diskLong"
+            units_key_lat = "diskLat"
+
+        xds["SUB_OBSERVER_POSITION"].attrs.update(
+            {
+                "type": "location",
+                "ellipsoid": "NA",
+                "origin_object_name": ephemeris_meta["NAME"],
+                "coordinate_system": "planetodetic",
+                "units": [
+                    cast_to_str(
+                        ephemris_column_description[units_key_lon]["keywords"][
+                            unit_keyword
+                        ]
+                    ),
+                    cast_to_str(
+                        ephemris_column_description[units_key_lat]["keywords"][
+                            unit_keyword
+                        ]
+                    ),
+                    "m",
+                ],
+            }
+        )
+
+    if "si_lon" in ephemeris_xds.data_vars:
+        xds["SUB_SOLAR_POSITION"] = xr.DataArray(
+            np.column_stack(
+                (
+                    ephemeris_xds["si_lon"].data,
+                    ephemeris_xds["si_lat"].data,
+                    ephemeris_xds["r"].data,
+                )
+            ),
+            dims=["time", "ellipsoid_pos_label"],
+        )
+        xds["SUB_SOLAR_POSITION"].attrs.update(
+            {
+                "type": "location",
+                "ellipsoid": "NA",
+                "origin_object_name": "Sun",
+                "coordinate_system": "planetodetic",
+                "units": [
+                    cast_to_str(
+                        ephemris_column_description["SI_lon"]["keywords"][unit_keyword]
+                    ),
+                    cast_to_str(
+                        ephemris_column_description["SI_lat"]["keywords"][unit_keyword]
+                    ),
+                    cast_to_str(
+                        ephemris_column_description["r"]["keywords"][unit_keyword]
+                    ),
+                ],
+            }
+        )
+
+    if "rdot" in ephemeris_xds.data_vars:
+        xds["HELIOCENTRIC_RADIAL_VELOCITY"] = xr.DataArray(
+            ephemeris_xds["rdot"].data, dims=["time"]
+        )
+        xds["HELIOCENTRIC_RADIAL_VELOCITY"].attrs.update(
+            {
+                "type": "quantity",
+                "units": [
+                    cast_to_str(
+                        ephemris_column_description["rdot"]["keywords"][unit_keyword]
+                    )
+                ],
+            }
+        )
+
+    if "phang" in ephemeris_xds.data_vars:
+        xds["OBSERVER_PHASE_ANGLE"] = xr.DataArray(
+            ephemeris_xds["phang"].data, dims=["time"]
+        )
+        xds["OBSERVER_PHASE_ANGLE"].attrs.update(
+            {
+                "type": "quantity",
+                "units": [
+                    cast_to_str(
+                        ephemris_column_description["phang"]["keywords"][unit_keyword]
+                    )
+                ],
+            }
+        )
+
+    xds = xds.assign_coords(coords)
+    xds["time"].attrs.update(
+        {"type": "time", "units": ["s"], "scale": "UTC", "format": "UNIX"}
+    )
+
+    xds = convert_to_si_units(xds)
+    xds = interpolate_to_time(xds, interp_time, "field_and_source_xds")
+
+    if is_single_dish:
+        xds["FIELD_REFERENCE_CENTER"] = xr.DataArray(
+            add_position_offsets(
+                np.append(xds["FIELD_REFERENCE_CENTER"].data, 0),
+                xds["SOURCE_POSITION"].data,
+            ),
+            dims=["time", "sky_pos_label"],
+        )
+    else:
+        xds["FIELD_PHASE_CENTER"] = xr.DataArray(
+            add_position_offsets(
+                np.append(xds["FIELD_PHASE_CENTER"].data, 0),
+                xds["SOURCE_POSITION"].data,
+            ),
+            dims=["time", "sky_pos_label"],
+        )
+
+    xds["FIELD_PHASE_CENTER"].attrs.update(xds["SOURCE_POSITION"].attrs)
+
+    return xds
+
+
+def add_position_offsets(dv_1, dv_2):
+    new_pos = dv_1 + dv_2
+
+    while np.any(new_pos[:, 0] > np.pi) or np.any(new_pos[:, 0] < -np.pi):
+        new_pos[:, 0] = np.where(
+            new_pos[:, 0] > np.pi, new_pos[:, 0] - 2 * np.pi, new_pos[:, 0]
+        )
+        new_pos[:, 0] = np.where(
+            new_pos[:, 0] < -np.pi, new_pos[:, 0] + 2 * np.pi, new_pos[:, 0]
+        )
+
+    while np.any(new_pos[:, 1] > np.pi / 2) or np.any(new_pos[:, 1] < -np.pi / 2):
+        new_pos[:, 1] = np.where(
+            new_pos[:, 1] > np.pi / 2, new_pos[:, 1] - np.pi, new_pos[:, 1]
+        )
+        new_pos[:, 1] = np.where(
+            new_pos[:, 1] < -np.pi / 2, new_pos[:, 1] + np.pi, new_pos[:, 1]
+        )
+
+    return new_pos
+
+
+def convert_to_si_units(xds):
+    for data_var in xds.data_vars:
+        if "units" in xds[data_var].attrs:
+            for u_i, u in enumerate(xds[data_var].attrs["units"]):
+                if u == "km":
+                    xds[data_var][..., u_i] = xds[data_var][..., u_i] * 1e3
+                    xds[data_var].attrs["units"][u_i] = "m"
+                if u == "km/s":
+                    xds[data_var][..., u_i] = xds[data_var][..., u_i] * 1e3
+                    xds[data_var].attrs["units"][u_i] = "m/s"
+                if u == "deg":
+                    xds[data_var][..., u_i] = xds[data_var][..., u_i] * np.pi / 180
+                    xds[data_var].attrs["units"][u_i] = "rad"
+                if u == "Au" or u == "AU":
+                    xds[data_var][..., u_i] = xds[data_var][..., u_i] * 149597870700
+                    xds[data_var].attrs["units"][u_i] = "m"
+                if u == "Au/d" or u == "AU/d":
+                    xds[data_var][..., u_i] = (
+                        xds[data_var][..., u_i] * 149597870700 / 86400
+                    )
+                    xds[data_var].attrs["units"][u_i] = "m/s"
+                if u == "arcsec":
+                    xds[data_var][..., u_i] = xds[data_var][..., u_i] * np.pi / 648000
+                    xds[data_var].attrs["units"][u_i] = "rad"
+    return xds
+
+
+def extract_source_info(xds, path, is_ephemeris, source_id, spectral_window_id):
+    """
+    Extracts source information from the given path and adds it to the xarray dataset.
+
+    Parameters:
+    ----------
+    xds : xr.Dataset
+        The xarray dataset to which the source information will be added.
+    path : str
+        The path to the input file.
+    is_ephemeris : bool
+        Flag indicating if the source is an ephemeris.
+    source_id : int
+        The ID of the source.
+    spectral_window_id : int
+        The ID of the spectral window.
+
+    Returns:
+    -------
+    xds : xr.Dataset
+        The xarray dataset with the added source information.
+    """
+
+    if source_id == -1:
+        logger.warning(
+            f"Source_id is -1. No source information will be included in the field_and_source_xds."
+        )
+        xds.attrs["source_name"] = "None"
+        return xds
+
+    source_xds = read_generic_table(
+        path,
+        "SOURCE",
+        ignore=["SOURCE_MODEL"],  # Trying to read SOURCE_MODEL causes an error.
+        taql_where=f"where SOURCE_ID = {source_id} AND SPECTRAL_WINDOW_ID = {spectral_window_id}",
+    )
+
+    if len(source_xds.data_vars) == 0:  # The source xds is empty.
+        logger.warning(
+            f"SOURCE table empty for source_id {source_id} and spectral_window_id {spectral_window_id}."
+        )
+        xds.attrs["source_name"] = "None"
+        return xds
+
+    assert (
+        len(source_xds.source_id) == 1
+    ), "Can only process source table with a single source_id and spectral_window_id for a given MSv4 partition."
+    assert (
+        len(source_xds.spectral_window_id) == 1
+    ), "Can only process source table with a single source_id and spectral_window_id for a given MSv4 partition."
+    assert (
+        len(source_xds.time) == 1
+    ), "Can only process source table with a single time entry for a source_id and spectral_window_id."
+    source_xds = source_xds.isel(time=0, source_id=0, spectral_window_id=0)
+
+    xds.attrs["source_name"] = str(source_xds["name"].data)
+    xds.attrs["code"] = str(source_xds["code"].data)
+    source_column_description = source_xds.attrs["other"]["msv2"]["ctds_attrs"][
+        "column_descriptions"
+    ]
+
+    if not is_ephemeris:
+        msv4_measure = column_description_casacore_to_msv4_measure(
+            source_column_description["DIRECTION"]
+        )
+        xds["SOURCE_DIRECTION"] = xr.DataArray(
+            source_xds["direction"].data, dims=["sky_dir_label"]
+        )
+        xds["SOURCE_DIRECTION"].attrs.update(msv4_measure)
+
+        # msv4_measure = column_description_casacore_to_msv4_measure(
+        #     source_column_description["PROPER_MOTION"]
+        # )
+        # xds["SOURCE_PROPER_MOTION"] = xr.DataArray(
+        #     source_xds["proper_motion"].data, dims=["sky_dir_label"]
+        # )
+        # xds["SOURCE_PROPER_MOTION"].attrs.update(msv4_measure)
+
+    # ['DIRECTION', 'PROPER_MOTION', 'CALIBRATION_GROUP', 'CODE', 'INTERVAL', 'NAME', 'NUM_LINES', 'SOURCE_ID', 'SPECTRAL_WINDOW_ID', 'TIME', 'POSITION', 'TRANSITION', 'REST_FREQUENCY', 'SYSVEL']
+    if source_xds["num_lines"] > 0:
+        coords = {"line_name": source_xds["transition"].data}
+        xds = xds.assign_coords(coords)
+
+        optional_data_variables = {
+            "rest_frequency": "LINE_REST_FREQUENCY",
+            "sysvel": "LINE_SYSTEMIC_VELOCITY",
+        }
+        for generic_name, msv4_name in optional_data_variables.items():
+            if generic_name in source_xds:
+                msv4_measure = column_description_casacore_to_msv4_measure(
+                    source_column_description[generic_name.upper()]
+                )
+                xds[msv4_name] = xr.DataArray(
+                    source_xds[generic_name].data, dims=["line_name"]
+                )
+                xds[msv4_name].attrs.update(msv4_measure)
+
+    # Need to add doppler info if present. Add check.
+    try:
+        doppler_xds = read_generic_table(
+            path,
+            "DOPPLER",
+        )
+        assert (
+            False
+        ), "Doppler table present. Please open an issue on https://github.com/casangi/xradio/issues so that we can addd support for this."
+    except:
+        pass
+
+    return xds
+
+
+def create_field_info_and_check_ephemeris(
+    field_and_source_xds, in_file, field_id, field_times, is_single_dish
+):
+    """
+    Create field information and check for ephemeris in the FIELD table folder.
+
+    Parameters:
+    ----------
+    field_and_source_xds : xr.Dataset
+        The xarray dataset to which the field and source information will be added.
+    in_file : str
+        The path to the input file.
+    field_id : int
+        The ID of the field.
+
+    Returns:
+    -------
+    field_and_source_xds : xr.Dataset
+        The xarray dataset with the added field and source information.
+    ephemeris_path : str
+        The path to the ephemeris table.
+    ephemeris_table_name : str
+        The name of the ephemeris table.
+    """
+    field_xds = read_generic_table(
+        in_file,
+        "FIELD",
+        rename_ids=subt_rename_ids["FIELD"],
+    )  # .sel(field_id=field_id)
+    # print('1****',field_xds)
+    # field_xds['field_id'] = np.arange(len(field_xds.field_id))
+
+    assert len(field_xds.poly_id) == 1, "Polynomial field positions not supported."
+    field_xds = field_xds.isel(poly_id=0)
+    field_xds = field_xds.sel(field_id=field_id)
+
+    from xradio._utils.array import check_if_consistent
+
+    source_id = check_if_consistent(field_xds.source_id, "source_id")
+
+    # print('source_id', source_id)
+    # print(field_xds)
+    # print('***')
+    # print(field_xds.field_id)
+    # print('***')
+    # print(field_id)
+
+    field_and_source_xds.attrs.update(
+        {
+            "field_name": str(field_xds["name"].data),
+            "field_code": str(field_xds["code"].data),
+            # "field_id": field_id,
+            "source_id": source_id,
+        }
+    )
+
+    ephemeris_table_name = None
+    ephemeris_path = None
+    is_ephemeris = False
+
+    # Need to check if ephemeris_id is present and if epehemeris table is present.
+    if "ephemeris_id" in field_xds:
+        ephemeris_id = check_if_consistent(
+            field_xds.ephemeris_id, "ephemeris_id"
+        )  # int(field_xds["ephemeris_id"].data)
+        if ephemeris_id > -1:
+            files = os.listdir(os.path.join(in_file, "FIELD"))
+            ephemeris_table_name_start = "EPHEM" + str(ephemeris_id)
+
+            ephemeris_name_table_index = [
+                i for i in range(len(files)) if ephemeris_table_name_start in files[i]
+            ]
+            assert len(ephemeris_name_table_index) == 1, (
+                "More than one ephemeris table which starts with "
+                + ephemeris_table_name_start
+            )
+
+            if len(ephemeris_name_table_index) > 0:  # Are there any ephemeris tables.
+                is_ephemeris = True
+                e_index = ephemeris_name_table_index[0]
+                ephemeris_path = os.path.join(in_file, "FIELD")
+                ephemeris_table_name = files[e_index]
+            else:
+                logger.warning(
+                    f"Could not find ephemeris table for field_id {field_id}. Ephemeris information will not be included in the field_and_source_xds."
+                )
+
+    # if is_ephemeris:
+    #     field_data_variables = {
+    #         "delay_dir": "FIELD_DELAY_CENTER_OFFSET",
+    #         "phase_dir": "FIELD_PHASE_CENTER_OFFSET",
+    #         "reference_dir": "FIELD_REFERENCE_CENTER_OFFSET",
+    #     }
+    #     field_measures_type = "sky_coord_offset"
+    #     field_and_source_xds.attrs["field_and_source_xds_type"] = "ephemeris"
+    # else:
+    #     field_data_variables = {
+    #         "delay_dir": "FIELD_DELAY_CENTER",
+    #         "phase_dir": "FIELD_PHASE_CENTER",
+    #         "reference_dir": "FIELD_REFERENCE_CENTER",
+    #     }
+    #     field_measures_type = "sky_coord"
+    #     field_and_source_xds.attrs["field_and_source_xds_type"] = "standard"
+
+    if is_single_dish:
+        field_data_variables = {
+            "reference_dir": "FIELD_REFERENCE_CENTER",
+        }
+    else:
+        field_data_variables = {
+            # "delay_dir": "FIELD_DELAY_CENTER",
+            "phase_dir": "FIELD_PHASE_CENTER",
+            # "reference_dir": "FIELD_REFERENCE_CENTER",
+        }
+
+    field_measures_type = "sky_coord"
+
+    coords = {}
+    coords["sky_dir_label"] = ["ra", "dec"]
+    field_column_description = field_xds.attrs["other"]["msv2"]["ctds_attrs"][
+        "column_descriptions"
+    ]  # Keys are ['DELAY_DIR', 'PHASE_DIR', 'REFERENCE_DIR', 'CODE', 'FLAG_ROW', 'NAME', 'NUM_POLY', 'SOURCE_ID', 'TIME']
+
+    coords = {}
+    coords["sky_dir_label"] = ["ra", "dec"]
+    if field_times is not None:
+        coords["time"] = field_times
+        dims = ["time", "sky_dir_label"]
+    else:
+        dims = ["sky_dir_label"]
+
+    for generic_name, msv4_name in field_data_variables.items():
+
+        if field_xds.get("delaydir_ref") is None:
+            delaydir_ref = None
+        else:
+            delaydir_ref = check_if_consistent(
+                field_xds.get("delaydir_ref"), "delaydir_ref"
+            )
+        msv4_measure = column_description_casacore_to_msv4_measure(
+            field_column_description[generic_name.upper()], ref_code=delaydir_ref
+        )
+
+        field_and_source_xds[msv4_name] = xr.DataArray.from_dict(
+            {
+                "dims": dims,
+                "data": list(field_xds[generic_name].data),
+                "attrs": msv4_measure,
+            }
+        )
+
+        field_and_source_xds[msv4_name].attrs["type"] = field_measures_type
+
+    field_and_source_xds = field_and_source_xds.assign_coords(coords)
+    return field_and_source_xds, ephemeris_path, ephemeris_table_name