open_image and load_image fails with twhya image

[Jan-Willem]

Using main branch 20 March:

    img_xds  = open_image({"visibility_normalization":  "twhya_selfcal_5chans_lsrk_robust_0.5_niter_0.sumwt"})
    img_xds  = load_image({"visibility_normalization":  "twhya_selfcal_5chans_lsrk_robust_0.5_niter_0.sumwt"})

fails with ValueError: different number of dimensions on data and dims: 5 vs 3. The image can be opened using carta.

I mage is attached.

twhya_selfcal_5chans_lsrk_robust_0.5_niter_0.sumwt.zip

[dmehring]

open_image works for me using current main, so if it fails for some cases on current main, I'm going to need a bit more detail on how to reproduce that.

In [1]: from xradio.image import open_image, load_image
In [5]: img_xds  = open_image({"visibility_normalization":  "twhya_selfcal_5chans_lsrk_robust_0.5_niter_0.sumwt"})
Successful readonly open of default-locked table twhya_selfcal_5chans_lsrk_robust_0.5_niter_0.sumwt: 1 columns, 1 rows

In [6]: img_xds
Out[6]: 
<xarray.Dataset> Size: 112B
Dimensions:                   (time: 1, frequency: 5, polarization: 1)
Coordinates:
  * time                      (time) float64 8B 5.625e+04
  * frequency                 (frequency) float64 40B 3.727e+11 ... 3.727e+11
  * polarization              (polarization) <U1 4B 'I'
    velocity                  (frequency) float64 40B 1.149e+04 ... 9.531e+03
Data variables:
    VISIBILITY_NORMALIZATION  (time, frequency, polarization) float32 20B dask.array<chunksize=(1, 5, 1), meta=np.ndarray>
Attributes:
    coordinate_system_info:  {'reference_direction': {'attrs': {'frame': 'fk5...
    type:                    image_dataset
    data_groups:             {}

load_image does fail so I will look at that

[Jan-Willem]

Sorry, open_image is working.