xradio 0.0.24__py3-none-any.whl → 0.0.26__py3-none-any.whl

xradio/_utils/zarr/common.py

@@ -1,9 +1,10 @@
 import copy
 import xarray as xr
 import zarr
+import s3fs
 
 
-def _open_dataset(store, xds_isel=None, data_variables=None, load=False):
+def _open_dataset(store, xds_isel=None, data_variables=None, load=False, **kwargs):
     """
 
     Parameters
@@ -25,7 +26,11 @@ def _open_dataset(store, xds_isel=None, data_variables=None, load=False):
 
     import dask
 
-    xds = xr.open_zarr(store)
+    if "s3" in kwargs.keys():
+        mapping = s3fs.S3Map(root=store, s3=kwargs["s3"], check=False)
+        xds = xr.open_zarr(store=mapping)
+    else:
+        xds = xr.open_zarr(store)
 
     if xds_isel is not None:
         xds = xds.isel(xds_isel)

xradio/image/_util/_casacore/xds_from_casacore.py

@@ -106,7 +106,7 @@ def _add_sky_or_aperture(
         unit = casa_image.unit()
     xda.attrs[_image_type] = image_type
     xda.attrs["units"] = unit
-    name = "sky" if has_sph_dims else "aperture"
+    name = "SKY" if has_sph_dims else "APERTURE"
     xda = xda.rename(name)
     xds[xda.name] = xda
     return xds

xradio/image/_util/_casacore/xds_to_casacore.py

@@ -8,7 +8,7 @@ from astropy.coordinates import Angle
 from casacore import tables
 
 from .common import _active_mask, _create_new_image, _object_name, _pointing_center
-from ..common import _compute_sky_reference_pixel, _doppler_types
+from ..common import _aperture_or_sky, _compute_sky_reference_pixel, _doppler_types
 from ...._utils._casacore.tables import open_table_rw
 
 
@@ -168,7 +168,7 @@ def _history_from_xds(xds: xr.Dataset, image: str) -> None:
 
 def _imageinfo_dict_from_xds(xds: xr.Dataset) -> dict:
     ii = {}
-    ap_sky = "sky" if "l" in xds.coords else "aperture"
+    ap_sky = _aperture_or_sky(xds)
     ii["image_type"] = (
         xds[ap_sky].attrs["image_type"] if "image_type" in xds[ap_sky].attrs else ""
     )
@@ -208,12 +208,12 @@ def _imageinfo_dict_from_xds(xds: xr.Dataset) -> dict:
 
 
 def _write_casa_data(xds: xr.Dataset, image_full_path: str) -> None:
-    sky_ap = "sky" if "sky" in xds else "aperture"
+    sky_ap = _aperture_or_sky(xds)
     if xds[sky_ap].shape[0] != 1:
         raise RuntimeError("XDS can only be converted if it has exactly one time plane")
     trans_coords = (
         ("frequency", "polarization", "m", "l")
-        if sky_ap == "sky"
+        if sky_ap == "SKY"
         else ("frequency", "polarization", "v", "u")
     )
     casa_image_shape = xds[sky_ap].isel(time=0).transpose(*trans_coords).shape[::-1]
@@ -312,7 +312,7 @@ def _write_initial_image(
 ) -> None:
     if not maskname:
         maskname = ""
-    for dv in ["sky", "aperture"]:
+    for dv in ["SKY", "APERTURE"]:
         if dv in xds.data_vars:
             value = xds[dv][0, 0, 0, 0, 0].values.item()
             if xds[dv][0, 0, 0, 0, 0].values.dtype == "float32":
@@ -352,7 +352,7 @@ def _write_pixels(
     value: xr.DataArray = None,
 ) -> None:
     flip = False
-    if v == "sky" or v == "aperture":
+    if v == "SKY" or v == "APERTURE":
         filename = image_full_path
     else:
         # mask

xradio/image/_util/_fits/xds_from_fits.py

@@ -685,7 +685,7 @@ def _add_sky_or_aperture(
     unit = helpers["bunit"]
     xda.attrs[_image_type] = image_type
     xda.attrs["units"] = unit
-    name = "sky" if has_sph_dims else "aperture"
+    name = "SKY" if has_sph_dims else "APERTURE"
     xda = xda.rename(name)
     xds[xda.name] = xda
     if helpers["has_mask"]:

xradio/image/_util/_zarr/zarr_low_level.py

@@ -109,19 +109,55 @@ def write_binary_blob_to_disk(arr, file_path, compressor):
     Returns:
     - None
     """
+    import graphviper.utils.logger as logger
     # Encode the NumPy array using the codec
+    logger.debug('1. Before compressor ' + file_path)
     compressed_arr = compressor.encode(np.ascontiguousarray(arr))
 
+    logger.debug('2. Before makedir')
     # Ensure the directory exists before saving the file
     os.makedirs(os.path.dirname(file_path), exist_ok=True)
 
-    # Save the compressed array to disk
-    with open(file_path, "wb") as file:
-        file.write(compressed_arr)
+    arr_len = len(compressed_arr)
+    logger.debug('3. Before write the len is: ' + str(arr_len))
+    #Save the compressed array to disk
+    # with open(file_path, "wb") as file:
+    #     file.write(compressed_arr)
+    
+    logger.debug('4. Using new writer: ' + str(arr_len))
+    write_to_lustre_chunked(file_path, compressed_arr)
+        
+    # /.lustre/aoc/sciops/pford/CHILES/cube_image/uid___A002_Xee7674_X2844_Cube_3.img.zarr/SKY/0.0.110.0.0
+    # 348192501 bytes
+    # 332.0622453689575 M
+    
+    # from io import BufferedWriter
+    # # Calculate buffer size based on compressed_arr size (adjust multiplier)
+    # buffer_size = min(len(compressed_arr), 1024 * 1024 * 4)  # Max 4 MB buffer
+    # with BufferedWriter(open(file_path, "wb"), buffer_size) as f:
+    #     f.write(compressed_arr)
+    #     f.flush()  # Ensure data gets written to disk 
+    
+    
+    logger.debug('4. Write completed')
 
     # print(f"Compressed array saved to {file_path}")
 
 
+def write_to_lustre_chunked(file_path, compressed_arr, chunk_size=1024 * 1024 * 128):  # 128 MiB chunks
+  """
+  Writes compressed data to a Lustre file path with chunking.
+
+  Args:
+      file_path: Path to the file for writing.
+      compressed_arr: Compressed data array to write.
+      chunk_size: Size of each data chunk in bytes (default: 128 MiB).
+  """
+  with open(file_path, "wb") as f:
+    for i in range(0, len(compressed_arr), chunk_size):
+      chunk = compressed_arr[i:i + chunk_size]
+      f.write(chunk)
+
 def read_binary_blob_from_disk(file_path, compressor, dtype=np.float64):
     """
     Read a compressed binary blob from disk and decode it using Blosc.

xradio/image/_util/casacore.py

@@ -30,7 +30,7 @@ from ._casacore.xds_to_casacore import (
     _imageinfo_dict_from_xds,
     _write_casa_data,
 )
-from .common import _get_xds_dim_order, _dask_arrayize_dv
+from .common import _aperture_or_sky, _get_xds_dim_order, _dask_arrayize_dv
 
 warnings.filterwarnings("ignore", category=FutureWarning)
 
@@ -104,7 +104,7 @@ def _xds_to_casa_image(xds: xr.Dataset, imagename: str) -> None:
     image_full_path = os.path.expanduser(imagename)
     _write_casa_data(xds, image_full_path)
     # create coordinates
-    ap_sky = "sky" if "sky" in xds.data_vars else "aperture"
+    ap_sky = _aperture_or_sky(xds)
     coord = _coord_dict_from_xds(xds)
     ii = _imageinfo_dict_from_xds(xds)
     units = xds[ap_sky].attrs["units"] if "units" in xds[ap_sky].attrs else None

xradio/image/_util/common.py

@@ -13,6 +13,10 @@ _doppler_types = ["RADIO", "Z", "RATIO", "BETA", "GAMMA"]
 _image_type = "image_type"
 
 
+def _aperture_or_sky(xds: xr.Dataset) -> str:
+    return "SKY" if "SKY" in xds.data_vars or "l" in xds.coords else "APERTURE"
+
+
 def _get_xds_dim_order(has_sph: bool) -> list:
     dimorder = ["time", "polarization", "frequency"]
     dir_lin = ["l", "m"] if has_sph else ["u", "v"]

xradio/vis/_processing_set.py

@@ -43,7 +43,10 @@ class processing_set(dict):
             "start_frequency": [],
             "end_frequency": [],
             "shape": [],
+            "field_coords": []
         }
+        from astropy.coordinates import SkyCoord
+        import astropy.units as u
         for key, value in self.items():
             summary_data["name"].append(key)
             summary_data["ddi"].append(value.attrs["ddi"])
@@ -65,6 +68,12 @@ class processing_set(dict):
             )
             summary_data["start_frequency"].append(value["frequency"].values[0])
             summary_data["end_frequency"].append(value["frequency"].values[-1])
+
+            ra_dec_rad = value[data_name].attrs["field_info"]['phase_direction']['data']
+            frame = value[data_name].attrs["field_info"]['phase_direction']['attrs']['frame'].lower()
+            coord = SkyCoord(ra=ra_dec_rad[0]*u.rad, dec=ra_dec_rad[1]*u.rad, frame=frame)
+
+            summary_data["field_coords"].append([frame, coord.ra.to_string(unit=u.hour), coord.dec.to_string(unit=u.deg)])
         summary_df = pd.DataFrame(summary_data)
         return summary_df
 

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

@@ -195,10 +195,8 @@ def read_main_table_chunks(
     n_baseline_chunks = chunks[1]
     # loop over time chunks
     for time_chunk in range(0, n_unique_times, n_time_chunks):
-        time_start = (unique_times[time_chunk] - tol,)
-        time_end = (
-            unique_times[min(n_unique_times, time_chunk + n_time_chunks) - 1] + tol
-        )
+        time_start = unique_times[time_chunk] - tol
+        time_end = unique_times[min(n_unique_times, time_chunk + n_time_chunks) - 1] + tol
 
         # chunk time length
         ctlen = min(n_unique_times, time_chunk + n_time_chunks) - time_chunk

xradio/vis/load_processing_set.py

@@ -11,13 +11,13 @@ def load_processing_set(
     sel_parms: dict,
     data_variables: Union[list, None] = None,
     load_sub_datasets: bool = True,
-)->processing_set:
+) -> processing_set:
     """Loads a processing set into memory.
 
     Parameters
     ----------
     ps_store : str
-        String of the path and name of the processing set. For example '/users/user_1/uid___A002_Xf07bba_Xbe5c_target.lsrk.vis.zarr'.
+        String of the path and name of the processing set. For example '/users/user_1/uid___A002_Xf07bba_Xbe5c_target.lsrk.vis.zarr' for a file stored on a local file system, or 's3://viper-test-data/Antennae_North.cal.lsrk.split.vis.zarr/' for a file in AWS object storage.
     sel_parms : dict
         A dictionary where the keys are the names of the ms_xds's and the values are slice_dicts.
         slice_dicts: A dictionary where the keys are the dimension names and the values are slices.
@@ -37,26 +37,69 @@ def load_processing_set(
     Returns
     -------
     processing_set
-        In memory representation of processing set (data is represented by Dask.arrays). 
-    """    
+        In memory representation of processing set (data is represented by Dask.arrays).
+    """
     from xradio._utils.zarr.common import _open_dataset
+    import s3fs
+    from botocore.exceptions import NoCredentialsError
 
+    s3 = None
     ps = processing_set()
     for ms_dir_name, ms_xds_isel in sel_parms.items():
-        xds = _open_dataset(
-            os.path.join(ps_store, ms_dir_name, "MAIN"),
-            ms_xds_isel,
-            data_variables,
-            load=True,
-        )
-
-        if load_sub_datasets:
-            from xradio.vis.read_processing_set import _read_sub_xds
-
-            xds.attrs = {
-                **xds.attrs,
-                **_read_sub_xds(os.path.join(ps_store, ms_dir_name), load=True),
-            }
+
+        # before the _open_dataset call, check if dealing with an S3 bucket URL
+        if ps_store.startswith("s3"):
+            if not ps_store.endswith("/"):
+                # just for consistency, as there is no os.path equivalent in s3fs
+                ps_store = ps_store + "/"
+
+            try:
+                s3 = s3fs.S3FileSystem(anon=False, requester_pays=False)
+
+                main_xds = ps_store + ms_dir_name + "/MAIN"
+                xds = _open_dataset(
+                    main_xds, ms_xds_isel, data_variables, load=True, s3=s3
+                )
+
+                if load_sub_datasets:
+                    from xradio.vis.read_processing_set import _read_sub_xds
+
+                    xds.attrs = {
+                        **xds.attrs,
+                        **_read_sub_xds(
+                            os.path.join(ps_store, ms_dir_name), load=True, s3=s3
+                        ),
+                    }
+
+            except (NoCredentialsError, PermissionError) as e:
+                # only public, read-only buckets will be accessible
+                s3 = s3fs.S3FileSystem(anon=True)
+
+                main_xds = ps_store + ms_dir_name + "/MAIN"
+                xds = _open_dataset(
+                    main_xds, ms_xds_isel, data_variables, load=True, s3=s3
+                )
+
+                if load_sub_datasets:
+                    from xradio.vis.read_processing_set import _read_sub_xds
+
+                    xds.attrs = {
+                        **xds.attrs,
+                        **_read_sub_xds(
+                            os.path.join(ps_store, ms_dir_name), load=True, s3=s3
+                        ),
+                    }
+        else:
+            # fall back to the default case of assuming the files are on local disk
+            main_xds = os.path.join(ps_store, ms_dir_name, "MAIN")
+            xds = _open_dataset(main_xds, ms_xds_isel, data_variables, load=True)
+            if load_sub_datasets:
+                from xradio.vis.read_processing_set import _read_sub_xds
+
+                xds.attrs = {
+                    **xds.attrs,
+                    **_read_sub_xds(os.path.join(ps_store, ms_dir_name), load=True),
+                }
 
         ps[ms_dir_name] = xds
     return ps
@@ -94,7 +137,7 @@ class processing_set_iterator:
             The list of data variables to load into memory for example ['VISIBILITY', 'WEIGHT, 'FLAGS']. By default None which will load all data variables into memory.
         load_sub_datasets : bool, optional
             If true sub-datasets (for example weather_xds, antenna_xds, pointing_xds, ...) will be loaded into memory, by default True.
-        """        
+        """
 
         self.input_data = input_data
         self.input_data_store = input_data_store
@@ -124,4 +167,4 @@ class processing_set_iterator:
         else:
             xds = self.input_data[xds_name]  # In memory
 
-        return xds
+        return xds

xradio/vis/read_processing_set.py

@@ -3,11 +3,13 @@ import xarray as xr
 from ._processing_set import processing_set
 import graphviper.utils.logger as logger
 from xradio._utils.zarr.common import _open_dataset
+import s3fs
+from botocore.exceptions import NoCredentialsError
 
 
 def read_processing_set(
     ps_store: str, intents: list = None, fields: str = None
-)->processing_set:
+) -> processing_set:
     """Creates a lazy representation of a Processing Set (only meta-data is loaded into memory).
 
     Parameters
@@ -18,44 +20,101 @@ def read_processing_set(
         A list of the intents to be read for example ['OBSERVE_TARGET#ON_SOURCE']. The intents in a processing set can be seem by calling processing_set.summary().
         By default None, which will read all intents.
     fields : str, optional
-       The list of field names that will be read, by default None which will read all fields. 
-       
+       The list of field names that will be read, by default None which will read all fields.
+
     Returns
     -------
     processing_set
-        Lazy representation of processing set (data is represented by Dask.arrays). 
-    """    
-    items = os.listdir(ps_store)
+        Lazy representation of processing set (data is represented by Dask.arrays).
+    """
+    s3 = None
+    ps_store_is_s3dir = None
+
+    if os.path.isdir(ps_store):
+        ps_store_is_s3dir = False
+        # default to assuming the data are accessible on local file system
+        items = os.listdir(ps_store)
+
+    elif ps_store.startswith("s3"):
+        # only if not found locally, check if dealing with an S3 bucket URL
+        ps_store_is_s3dir = True
+        if not ps_store.endswith("/"):
+            # just for consistency, as there is no os.path equivalent in s3fs
+            ps_store = ps_store + "/"
+
+        try:
+            # initialize the S3 "file system", first attempting to use pre-configured credentials
+            s3 = s3fs.S3FileSystem(anon=False, requester_pays=False)
+
+            items = [bd.split(sep="/")[-1] for bd in s3.listdir(ps_store, detail=False)]
+
+        except (NoCredentialsError, PermissionError) as e:
+            # only public, read-only buckets will be accessible
+            # we will want to add messaging and error handling here
+            s3 = s3fs.S3FileSystem(anon=True)
+
+            items = [bd.split(sep="/")[-1] for bd in s3.listdir(ps_store, detail=False)]
+
+    else:
+        raise (
+            FileNotFoundError,
+            f"Could not find {ps_store} either locally or in the cloud.",
+        )
+
     ms_xds = xr.Dataset()
     ps = processing_set()
-    data_group = 'base'
+    data_group = "base"
     for ms_dir_name in items:
         if "ddi" in ms_dir_name:
-            xds = _open_dataset(os.path.join(ps_store, ms_dir_name, "MAIN"))
+            if ps_store_is_s3dir:
+                store_path = ps_store + ms_dir_name
+                store_path_main = store_path + "/MAIN"
+            else:
+                store_path_main = os.path.join(ps_store, ms_dir_name, "MAIN")
+                store_path = os.path.split(store_path_main)[0]
+            if s3 is not None:
+                xds = _open_dataset(store_path_main, s3=s3)
+            else:
+                xds = _open_dataset(store_path_main)
+
             if (intents is None) or (xds.attrs["intent"] in intents):
                 data_name = _get_data_name(xds, data_group)
 
                 if (fields is None) or (
                     xds[data_name].attrs["field_info"]["name"] in fields
                 ):
-                    xds.attrs = {
-                        **xds.attrs,
-                        **_read_sub_xds(os.path.join(ps_store, ms_dir_name)),
-                    }
-                    ps[ms_dir_name] = xds
+                    if s3 is not None:
+                        xds.attrs = {
+                            **xds.attrs,
+                            **_read_sub_xds(store_path, s3=s3),
+                        }
+                        ps[ms_dir_name] = xds
+                    else:
+                        xds.attrs = {
+                            **xds.attrs,
+                            **_read_sub_xds(store_path),
+                        }
+                        ps[ms_dir_name] = xds
+
     return ps
 
 
-def _read_sub_xds(ms_store, load=False):
+def _read_sub_xds(ms_store, load=False, **kwargs):
     sub_xds_dict = {}
 
     sub_xds = {
         "antenna_xds": "ANTENNA",
     }
     for sub_xds_key, sub_xds_name in sub_xds.items():
-        sub_xds_dict[sub_xds_key] = _open_dataset(
-            os.path.join(ms_store, sub_xds_name), load=load
-        )
+        if "s3" in kwargs.keys():
+            joined_store = ms_store + "/" + sub_xds_name
+            sub_xds_dict[sub_xds_key] = _open_dataset(
+                joined_store, load=load, s3=kwargs["s3"]
+            )
+        else:
+            sub_xds_dict[sub_xds_key] = _open_dataset(
+                os.path.join(ms_store, sub_xds_name), load=load
+            )
 
     optional_sub_xds = {
         "weather_xds": "WEATHER",
@@ -65,6 +124,12 @@ def _read_sub_xds(ms_store, load=False):
         sub_xds_path = os.path.join(ms_store, sub_xds_name)
         if os.path.isdir(sub_xds_path):
             sub_xds_dict[sub_xds_key] = _open_dataset(sub_xds_path, load=load)
+        elif "s3" in kwargs.keys():
+            joined_store = ms_store + "/" + sub_xds_name
+            if kwargs["s3"].isdir(joined_store):
+                sub_xds_dict[sub_xds_key] = _open_dataset(
+                    joined_store, load=load, s3=kwargs["s3"]
+                )
 
     return sub_xds_dict
 

{xradio-0.0.24.dist-info → xradio-0.0.26.dist-info}/METADATA

@@ -1,7 +1,7 @@
 Metadata-Version: 2.1
 Name: xradio
-Version: 0.0.24
-Summary: Xarray Radio Astronomy Data IO 
+Version: 0.0.26
+Summary: Xarray Radio Astronomy Data IO
 Author-email: Jan-Willem Steeb <jsteeb@nrao.edu>
 License: BSD 3-Clause License
         
@@ -32,7 +32,7 @@ License: BSD 3-Clause License
         OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT OF THE USE
         OF THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE.
         
-Requires-Python: <3.12,>=3.8
+Requires-Python: <3.12,>=3.9
 Description-Content-Type: text/markdown
 License-File: LICENSE.txt
 Requires-Dist: astropy
@@ -47,6 +47,7 @@ Requires-Dist: prettytable
 Requires-Dist: pytest
 Requires-Dist: pytest-cov
 Requires-Dist: pytest-html
+Requires-Dist: s3fs
 Requires-Dist: scipy
 Requires-Dist: tqdm
 Requires-Dist: xarray

{xradio-0.0.24.dist-info → xradio-0.0.26.dist-info}/RECORD

@@ -3,24 +3,24 @@ xradio/_utils/__init__.py,sha256=47DEQpj8HBSa-_TImW-5JCeuQeRkm5NMpJWZG3hSuFU,0
 xradio/_utils/common.py,sha256=bjUZfZZrYTOt0i_XVfmQ2kvgr3egoYPWKGgnr4vKe-Y,46
 xradio/_utils/_casacore/tables.py,sha256=aq6E_4RRAHdTBCwMKrVil1cWhFU2O980DNH9IlRKXLw,1280
 xradio/_utils/zarr/__init__.py,sha256=47DEQpj8HBSa-_TImW-5JCeuQeRkm5NMpJWZG3hSuFU,0
-xradio/_utils/zarr/common.py,sha256=HzZOky2XNuKkpRi6M4ZQCsTpl9ngiPtH9a6Pw63s4EI,3291
+xradio/_utils/zarr/common.py,sha256=w_orsJ7-nh8nXFllkYPwbEYFS5B-KeMVWnhO3elX_PQ,3470
 xradio/image/__init__.py,sha256=HAD0GfopIbhdxOYckyW6S9US_dSWmZrwIl3FHUzZwrE,435
 xradio/image/image.py,sha256=QoJ_BTLoMfeXJzU1yvtidBIhaMmjNA5_-6C3FWJRUeI,15635
 xradio/image/_util/__init__.py,sha256=M9lxD1Gc7kv0ucDEDbjLRuIEuESev-IG8j9EaCKUAkA,77
-xradio/image/_util/casacore.py,sha256=Ozuq_VgseU1pcX230078TzOSsdYdAbwbjHCkwU79m0g,4242
-xradio/image/_util/common.py,sha256=uH2Po1SycsOK_KZqyLWllNPmRtjTSEr0r2YF4oPRZas,8793
+xradio/image/_util/casacore.py,sha256=DmBTHUQ6870N5ARuFnYSfjZSLniJYgsjrsICUlCREYM,4234
+xradio/image/_util/common.py,sha256=hzVcnq4pYhxvkbMAIBzhWKdLN4nrqRGlmFwABxkww7A,8921
 xradio/image/_util/fits.py,sha256=gyGm06fuCKqVGK7uv-ObvQNfFawUDsIOa_nQyklM3Aw,329
 xradio/image/_util/image_factory.py,sha256=6tPzs20FTm2wEshHc1xqtTV7D0TbKxGLUKAVtvOc68I,10506
 xradio/image/_util/zarr.py,sha256=xTjg-KY-T4vuyua8pvuZZjCL-rI_wAsPjPUOYd5PZr4,1512
 xradio/image/_util/_casacore/__init__.py,sha256=OlsiRE40o1jSbBI4khgQQzgfDYbAlOMKIhO4UFlbGhg,41
 xradio/image/_util/_casacore/common.py,sha256=ky999eTCWta8w-uIs-7P7rPhZRLuh9yTuQXAxPvaPm4,1579
-xradio/image/_util/_casacore/xds_from_casacore.py,sha256=Ubt5PCVLFX2t9YpCEH9jS0g0hMI56Owfr5Ibup7-Hg4,42562
-xradio/image/_util/_casacore/xds_to_casacore.py,sha256=wYam-VgsBfSE5i3YS0dmtuA4xwDpKOkl8r2nYtGFQzA,15357
-xradio/image/_util/_fits/xds_from_fits.py,sha256=SIYyl4dAoCdiB4CSv0_wV_hY7jhLZZE2-39Yn32dd5Y,28339
+xradio/image/_util/_casacore/xds_from_casacore.py,sha256=Rht4A32QLAQ7uizwKfOsZ3Z819shvlUbsZfAuTIVerU,42562
+xradio/image/_util/_casacore/xds_to_casacore.py,sha256=P6c-yoOjuVQkm07ApA7FFKfje4aPwV-MsRFKaRaPq9I,15338
+xradio/image/_util/_fits/xds_from_fits.py,sha256=rsV2OMK3g4U_GYKgj68cZXWIlXi0M0oU2M9FtAUeCk8,28339
 xradio/image/_util/_zarr/common.py,sha256=apMX_bF4Hr3pFGjnDFpp36KgmhTYAPBZquNkjBHrsXk,307
 xradio/image/_util/_zarr/xds_from_zarr.py,sha256=hz6lHlpybfr_r8pn_uObDHOFmN5h75F11bkBv8KCuP0,3192
 xradio/image/_util/_zarr/xds_to_zarr.py,sha256=wogXbwX8n3Sl9PHoc3_Y_LBowQsQ-94HZQFZ5NcxUZA,1624
-xradio/image/_util/_zarr/zarr_low_level.py,sha256=DuAD1kVUbHjOxVH9jRagk9GdPA5uH2fIMzA0djxM6Sc,9537
+xradio/image/_util/_zarr/zarr_low_level.py,sha256=5D8Vu8QZD8CpKSPMkxHhnKpdk3vhUStW3kQFHUQUQns,10984
 xradio/schema/__init__.py,sha256=UpejQegOaCLrsbcR4MLudR0RxeE0sN3zxFXM8rzyJPo,444
 xradio/schema/bases.py,sha256=vcW47jZWpJ0mJdni7eFVY7zJoper2sy2VjX8LE3pUqc,150
 xradio/schema/check.py,sha256=3u79hRL3pGF6dvQE0LF21nGdAVnRXWwgnbMmStGBSSA,16310
@@ -28,11 +28,11 @@ xradio/schema/dataclass.py,sha256=vkc2cqLjGV5QN8j70GbBaNfslT_KLWmebsPGeBEuGcs,88
 xradio/schema/metamodel.py,sha256=RHrihyaetinu7_lGTTZ31Rlv-_Db_EgQCXzk56H004o,3476
 xradio/schema/typing.py,sha256=coF3LuKOlCUJGKTUUH81EcjePZ86koOYzm8qzsAw-HU,9983
 xradio/vis/__init__.py,sha256=AV2WG26NzFB1LEEtFaq1ULQKz9VnluEAjg0Qb5Ju7m8,358
-xradio/vis/_processing_set.py,sha256=3cm-WAuBjzX0_XPUjSR1i2i5_PeXiGwW12XioVji2Wk,3666
+xradio/vis/_processing_set.py,sha256=zYCswmWjsRN_8M61IkbMo6zpI-4ABMWpKZkL8f20KGo,4187
 xradio/vis/convert_msv2_to_processing_set.py,sha256=7vTjqtWFEBOySnLVoadceKCA4VLgOig7eCs1dxrdYQA,3966
-xradio/vis/load_processing_set.py,sha256=DsKjizvcJxXNn5BuXH2fye-dzMoB-BDK2f-8jsvK2QU,4776
+xradio/vis/load_processing_set.py,sha256=ldgg0GqFwYuNc2fPkRlc_lW2Kt7naq2wFGloALl0J4Y,6753
 xradio/vis/model.py,sha256=uBjvvhYEY1p-3H3NStrt1ZKMQACLGLo93OiEBvDVId8,17083
-xradio/vis/read_processing_set.py,sha256=ljQVi8njdsAU--V-3QaZ4_EhWXqjnW3mdLr0VE1N4cA,3043
+xradio/vis/read_processing_set.py,sha256=mhBcwrNEVlu5_dMfRaJQHzD-Yc3n2dtVqhNxvaoq09U,5522
 xradio/vis/vis_io.py,sha256=rCSOt4Max37uFzF3_Ck4U4xWzzYcipdbYcxbhBhQ_Qs,5278
 xradio/vis/_vis_utils/__init__.py,sha256=Scu6rKJ2SpO8aG7F-xdTZcYfyWx0viV8gFh8E8ur_gI,93
 xradio/vis/_vis_utils/ms.py,sha256=0uycYCDmeQku16TdPcnZEBJMdfb_i6xPoieYPhPoVIg,5258
@@ -52,7 +52,7 @@ xradio/vis/_vis_utils/_ms/subtables.py,sha256=hkSa3sXUaSRVEv2CNK65Svx5Uyp3opZHaH
 xradio/vis/_vis_utils/_ms/_tables/load.py,sha256=TXrEf7fKfsP0Wp7fB7DLWeOS_Z8A1OLKGWGNSoiYHOM,1665
 xradio/vis/_vis_utils/_ms/_tables/load_main_table.py,sha256=HyIfjaQXfFA30YRUzU2H_ugRnArQ1XOiVGKxC9PbiOo,13497
 xradio/vis/_vis_utils/_ms/_tables/read.py,sha256=lt3JOre-7WyzTprUuZIowi6aEcXc8R367oLZUB0vdhk,21781
-xradio/vis/_vis_utils/_ms/_tables/read_main_table.py,sha256=cb3vkzZUwYG5h2FOgGipF2z7DPDgx_B56rii0dt6KYk,23985
+xradio/vis/_vis_utils/_ms/_tables/read_main_table.py,sha256=9XasWazrGSkQ6C4jSrItIGidaTUzh-bdSRYkYR5j0Dg,23958
 xradio/vis/_vis_utils/_ms/_tables/read_subtables.py,sha256=_3AlxlLJwdu__a8sfG8a7-wz-yTo2S7JyVit17XaxKs,11741
 xradio/vis/_vis_utils/_ms/_tables/table_query.py,sha256=q8EGFf_zIwHcHnvFJOn8hPh8zFZQ3f7BGbXvL3bHad4,555
 xradio/vis/_vis_utils/_ms/_tables/write.py,sha256=Uuk4QGib_QOc8i6g7OYB9E6SxMUDPlXDYXlHqkzPX5Q,8475
@@ -64,8 +64,8 @@ xradio/vis/_vis_utils/_utils/xds_helper.py,sha256=UhtAZV5DyYzVVBkXzwDAOz6TICxotQ
 xradio/vis/_vis_utils/_zarr/encoding.py,sha256=GENIlThV6a9CUCL6gIGlu9c6NR3OFWNos6mpxZjMwDc,536
 xradio/vis/_vis_utils/_zarr/read.py,sha256=ikNGlOdHuZ_cgWpPAZ4iHzeLdU44I0iBLcqSEiM_hCk,7111
 xradio/vis/_vis_utils/_zarr/write.py,sha256=exvrqNVnVKk6LzbDPm_fm142YzX-7lHGqklXTQB9wh0,8864
-xradio-0.0.24.dist-info/LICENSE.txt,sha256=dvACd-5O67yjSZlnEKcWmu3DqwzBtbC922iPv0KOeAw,1516
-xradio-0.0.24.dist-info/METADATA,sha256=yJW0v_5BpEInE5IVWDbELOpAHGGzgT0a0PmrbLnWYZ8,4070
-xradio-0.0.24.dist-info/WHEEL,sha256=GJ7t_kWBFywbagK5eo9IoUwLW6oyOeTKmQ-9iHFVNxQ,92
-xradio-0.0.24.dist-info/top_level.txt,sha256=dQu27fGBZJ2Yk-gW5XeD-dZ76Xa4Xcvk60Vz-dwXp7k,7
-xradio-0.0.24.dist-info/RECORD,,
+xradio-0.0.26.dist-info/LICENSE.txt,sha256=dvACd-5O67yjSZlnEKcWmu3DqwzBtbC922iPv0KOeAw,1516
+xradio-0.0.26.dist-info/METADATA,sha256=Z3bKRq6dhmg8EIvgddHYafA9TSFoWd1PDerkpkkf8Gw,4089
+xradio-0.0.26.dist-info/WHEEL,sha256=GJ7t_kWBFywbagK5eo9IoUwLW6oyOeTKmQ-9iHFVNxQ,92
+xradio-0.0.26.dist-info/top_level.txt,sha256=dQu27fGBZJ2Yk-gW5XeD-dZ76Xa4Xcvk60Vz-dwXp7k,7
+xradio-0.0.26.dist-info/RECORD,,