cmpute_availabilitymatrix throws ValueError (coords is not dict-like)

Open oruhnau opened this issue 1 year ago • 0 comments

Version Checks (indicate both or one)

[X] I have confirmed this bug exists on the lastest release of Atlite.
[ ] I have confirmed this bug exists on the current master branch of Atlite.

Issue Description

I am trying to run the examples/landuse-availability.ipynb

The cell in which the availability matrix is calculated throws the following error:

ValueError Traceback (most recent call last) Cell In[11], line 1 ----> 1 A = cutout.availabilitymatrix(shapes, excluder) 2 A

File c:\Users\olive\miniconda3\envs\atlite-test\Lib\site-packages\atlite\gis.py:761, in compute_availabilitymatrix(cutout, shapes, excluder, nprocesses, disable_progressbar) 759 availability = np.stack(availability)[:, ::-1] # flip axis, see Notes 760 coords = [(shapes.index), ("y", cutout.data.y.data), ("x", cutout.data.x.data)] --> 761 return xr.DataArray(availability, coords=coords)

File c:\Users\olive\miniconda3\envs\atlite-test\Lib\site-packages\xarray\core\dataarray.py:479, in DataArray.init(self, data, coords, dims, name, attrs, indexes, fastpath) 477 data = _check_data_shape(data, coords, dims) 478 data = as_compatible_data(data) --> 479 coords, dims = _infer_coords_and_dims(data.shape, coords, dims) 480 variable = Variable(dims, data, attrs, fastpath=True) 482 if not isinstance(coords, Coordinates):

File c:\Users\olive\miniconda3\envs\atlite-test\Lib\site-packages\xarray\core\dataarray.py:163, in _infer_coords_and_dims(shape, coords, dims) 156 """All the logic for creating a new DataArray""" 158 if ( 159 coords is not None 160 and not utils.is_dict_like(coords) 161 and len(coords) != len(shape) 162 ): --> 163 raise ValueError( 164 f"coords is not dict-like, but it has {len(coords)} items, " 165 f"which does not match the {len(shape)} dimensions of the " 166 "data" 167 ) 169 if isinstance(dims, str): 170 dims = (dims,)

ValueError: coords is not dict-like, but it has 3 items, which does not match the 4 dimensions of the data

Any idea how this could be solved?

And in the meantime: Do you know a combination of package versions for which this issue should not occur?

Reproducible Example

examples/landuse-availability.ipynb

Expected Behavior

No response

Installed Versions

name: atlite-test channels:

conda-forge
defaults dependencies:
affine=2.4.0=pyhd8ed1ab_0
asttokens=2.4.1=pyhd8ed1ab_0
atlite=0.2.14=pyhd8ed1ab_0
attrs=24.2.0=pyh71513ae_0
aws-c-auth=0.7.31=hdb8b6b4_2
aws-c-cal=0.7.4=h891f644_2
aws-c-common=0.9.29=h2466b09_0
aws-c-compression=0.2.19=h891f644_2
aws-c-event-stream=0.4.3=h9d4aeba_4
aws-c-http=0.8.10=h4d3606f_2
aws-c-io=0.14.19=h5c5bb51_1
aws-c-mqtt=0.10.7=h98e454b_2
aws-c-s3=0.6.7=h23cec54_0
aws-c-sdkutils=0.1.19=h891f644_4
aws-checksums=0.1.20=h891f644_1
aws-crt-cpp=0.28.3=h97eb5cb_8
aws-sdk-cpp=1.11.407=h1b410ea_1
blosc=1.21.6=h85f69ea_0
bokeh=3.6.0=pyhd8ed1ab_0
bottleneck=1.4.2=py312h1a27103_0
branca=0.7.2=pyhd8ed1ab_0
brotli=1.1.0=h2466b09_2
brotli-bin=1.1.0=h2466b09_2
brotli-python=1.1.0=py312h275cf98_2
bzip2=1.0.8=h2466b09_7
c-ares=1.34.2=h2466b09_0
ca-certificates=2024.8.30=h56e8100_0
cads-api-client=1.4.4=pyhd8ed1ab_0
cdsapi=0.7.3=pyhd8ed1ab_0
certifi=2024.8.30=pyhd8ed1ab_0
cffi=1.17.1=py312h4389bb4_0
cftime=1.6.4=py312h1a27103_1
charset-normalizer=3.4.0=pyhd8ed1ab_0
click=8.1.7=win_pyh7428d3b_0
click-plugins=1.1.1=py_0
cligj=0.7.2=pyhd8ed1ab_1
cloudpickle=3.1.0=pyhd8ed1ab_1
colorama=0.4.6=pyhd8ed1ab_0
comm=0.2.2=pyhd8ed1ab_0
contourpy=1.3.0=py312hd5eb7cc_2
cpython=3.12.7=py312hd8ed1ab_0
cycler=0.12.1=pyhd8ed1ab_0
cytoolz=1.0.0=py312h4389bb4_1
dask=2024.10.0=pyhd8ed1ab_0
dask-core=2024.10.0=pyhd8ed1ab_0
dask-expr=1.1.16=pyhd8ed1ab_0
debugpy=1.8.7=py312h275cf98_0
decorator=5.1.1=pyhd8ed1ab_0
distributed=2024.10.0=pyhd8ed1ab_0
exceptiongroup=1.2.2=pyhd8ed1ab_0
executing=2.1.0=pyhd8ed1ab_0
folium=0.17.0=pyhd8ed1ab_0
fonttools=4.54.1=py312h4389bb4_0
freetype=2.12.1=hdaf720e_2
freexl=2.0.0=h8276f4a_0
fsspec=2024.9.0=pyhff2d567_0
geopandas=1.0.1=pyhd8ed1ab_1
geopandas-base=1.0.1=pyha770c72_1
geos=3.13.0=h5a68840_0
geotiff=1.7.3=h496ac4d_3
h2=4.1.0=pyhd8ed1ab_0
hdf4=4.2.15=h5557f11_7
hdf5=1.14.4=nompi_hd5d9e70_101
hpack=4.0.0=pyh9f0ad1d_0
hyperframe=6.0.1=pyhd8ed1ab_0
idna=3.10=pyhd8ed1ab_0
importlib-metadata=8.5.0=pyha770c72_0
importlib_metadata=8.5.0=hd8ed1ab_0
iniconfig=2.0.0=pyhd8ed1ab_0
intel-openmp=2024.2.1=h57928b3_1083
ipykernel=6.29.5=pyh4bbf305_0
ipython=8.28.0=pyh7428d3b_0
jedi=0.19.1=pyhd8ed1ab_0
jinja2=3.1.4=pyhd8ed1ab_0
joblib=1.4.2=pyhd8ed1ab_0
jupyter_client=8.6.3=pyhd8ed1ab_0
jupyter_core=5.7.2=pyh5737063_1
kiwisolver=1.4.7=py312hd5eb7cc_0
krb5=1.21.3=hdf4eb48_0
lcms2=2.16=h67d730c_0
lerc=4.0.0=h63175ca_0
libabseil=20240722.0=cxx17_he0c23c2_1
libaec=1.1.3=h63175ca_0
libarchive=3.7.4=haf234dc_0
libarrow=17.0.0=ha019072_22_cpu
libarrow-acero=17.0.0=hac47afa_22_cpu
libarrow-dataset=17.0.0=hac47afa_22_cpu
libarrow-substrait=17.0.0=ha9530af_22_cpu
libblas=3.9.0=18_win64_mkl
libbrotlicommon=1.1.0=h2466b09_2
libbrotlidec=1.1.0=h2466b09_2
libbrotlienc=1.1.0=h2466b09_2
libcblas=3.9.0=18_win64_mkl
libcrc32c=1.1.2=h0e60522_0
libcurl=8.10.1=h1ee3ff0_0
libdeflate=1.22=h2466b09_0
libevent=2.1.12=h3671451_1
libexpat=2.6.3=he0c23c2_0
libffi=3.4.2=h8ffe710_5
libgdal-core=3.9.2=h042995d_7
libgoogle-cloud=2.30.0=ha00044d_0
libgoogle-cloud-storage=2.30.0=he5eb982_0
libgrpc=1.65.5=ha20e22e_0
libhwloc=2.11.1=default_h8125262_1000
libiconv=1.17=hcfcfb64_2
libjpeg-turbo=3.0.0=hcfcfb64_1
libkml=1.3.0=h538826c_1021
liblapack=3.9.0=18_win64_mkl
libnetcdf=4.9.2=nompi_he239ae6_115
libparquet=17.0.0=h59f2d37_22_cpu
libpng=1.6.44=h3ca93ac_0
libprotobuf=5.27.5=hcaed137_2
libre2-11=2024.07.02=h4eb7d71_1
librttopo=1.1.0=hd4c2148_17
libsodium=1.0.20=hc70643c_0
libspatialite=5.1.0=h939089a_11
libsqlite=3.46.1=h2466b09_0
libssh2=1.11.0=h7dfc565_0
libthrift=0.21.0=hbe90ef8_0
libtiff=4.7.0=hfc51747_1
libutf8proc=2.8.0=h82a8f57_0
libwebp-base=1.4.0=hcfcfb64_0
libxcb=1.16=h013a479_1
libxml2=2.12.7=h0f24e4e_4
libzip=1.11.1=h25f2845_0
libzlib=1.3.1=h2466b09_2
locket=1.0.0=pyhd8ed1ab_0
lz4=4.3.3=py312h0608a1d_1
lz4-c=1.9.4=hcfcfb64_0
lzo=2.10=hcfcfb64_1001
m2w64-gcc-libgfortran=5.3.0=6
m2w64-gcc-libs=5.3.0=7
m2w64-gcc-libs-core=5.3.0=7
m2w64-gmp=6.1.0=2
m2w64-libwinpthread-git=5.0.0.4634.697f757=2
mapclassify=2.8.1=pyhd8ed1ab_0
markupsafe=3.0.2=py312h31fea79_0
matplotlib-base=3.9.2=py312h90004f6_1
matplotlib-inline=0.1.7=pyhd8ed1ab_0
minizip=4.0.6=hb638d1e_0
mkl=2022.1.0=h6a75c08_874
msgpack-python=1.1.0=py312hd5eb7cc_0
msys2-conda-epoch=20160418=1
multiurl=0.3.1=pyhd8ed1ab_0
munkres=1.1.4=pyh9f0ad1d_0
nest-asyncio=1.6.0=pyhd8ed1ab_0
netcdf4=1.7.1=nompi_py312h57e6fe7_103
networkx=3.4.1=pyhd8ed1ab_0
numexpr=2.8.7=mkl_py312h73deb7f_3
numpy=1.26.4=py312h8753938_0
openjpeg=2.5.2=h3d672ee_0
openssl=3.3.2=h2466b09_0
orc=2.0.2=h1c5a4bf_1
packaging=24.1=pyhd8ed1ab_0
pandas=2.2.3=py312h72972c8_1
parso=0.8.4=pyhd8ed1ab_0
partd=1.4.2=pyhd8ed1ab_0
pcre2=10.44=h3d7b363_2
pickleshare=0.7.5=py_1003
pillow=10.4.0=py312h381445a_1
pip=24.2=pyh8b19718_1
platformdirs=4.3.6=pyhd8ed1ab_0
pluggy=1.5.0=pyhd8ed1ab_0
progressbar2=4.5.0=pyhd8ed1ab_0
proj=9.5.0=hd9569ee_0
prompt-toolkit=3.0.48=pyha770c72_0
psutil=6.0.0=py312h4389bb4_2
pthread-stubs=0.4=hcd874cb_1001
pthreads-win32=2.9.1=h2466b09_4
pure_eval=0.2.3=pyhd8ed1ab_0
pyarrow=17.0.0=py312h7e22eef_1
pyarrow-core=17.0.0=py312h6a9c419_1_cpu
pycparser=2.22=pyhd8ed1ab_0
pygments=2.18.0=pyhd8ed1ab_0
pyogrio=0.10.0=py312h8705084_0
pyparsing=3.2.0=pyhd8ed1ab_1
pyproj=3.7.0=py312ha24589b_0
pysocks=1.7.1=pyh0701188_6
pytest=8.3.3=pyhd8ed1ab_0
python=3.12.7=hce54a09_0_cpython
python-dateutil=2.9.0=pyhd8ed1ab_0
python-tzdata=2024.2=pyhd8ed1ab_0
python-utils=3.8.2=pyhd8ed1ab_0
python_abi=3.12=5_cp312
pytz=2024.1=pyhd8ed1ab_0
pywin32=307=py312h275cf98_3
pyyaml=6.0.2=py312h4389bb4_1
pyzmq=26.2.0=py312hd7027bb_3
qhull=2020.2=hc790b64_5
rasterio=1.4.1=py312h1701b51_0
re2=2024.07.02=hd3b24a8_1
requests=2.32.3=pyhd8ed1ab_0
scikit-learn=1.5.2=py312h816cc57_1
scipy=1.14.1=py312h1f4e10d_0
setuptools=75.1.0=pyhd8ed1ab_0
shapely=2.0.6=py312h0c580ee_2
six=1.16.0=pyh6c4a22f_0
snappy=1.2.1=h23299a8_0
snuggs=1.4.7=pyhd8ed1ab_1
sortedcontainers=2.4.0=pyhd8ed1ab_0
sqlite=3.46.1=h2466b09_0
stack_data=0.6.2=pyhd8ed1ab_0
tbb=2021.13.0=hc790b64_0
tblib=3.0.0=pyhd8ed1ab_0
threadpoolctl=3.5.0=pyhc1e730c_0
tk=8.6.13=h5226925_1
tomli=2.0.2=pyhd8ed1ab_0
toolz=1.0.0=pyhd8ed1ab_0
tornado=6.4.1=py312h4389bb4_1
tqdm=4.66.5=pyhd8ed1ab_0
traitlets=5.14.3=pyhd8ed1ab_0
typing-extensions=4.12.2=hd8ed1ab_0
typing_extensions=4.12.2=pyha770c72_0
tzdata=2024b=hc8b5060_0
ucrt=10.0.22621.0=h57928b3_1
uriparser=0.9.8=h5a68840_0
urllib3=2.2.3=pyhd8ed1ab_0
vc=14.3=h8a93ad2_22
vc14_runtime=14.40.33810=hcc2c482_22
vs2015_runtime=14.40.33810=h3bf8584_22
wcwidth=0.2.13=pyhd8ed1ab_0
wheel=0.44.0=pyhd8ed1ab_0
win_inet_pton=1.1.0=pyh7428d3b_7
xarray=2024.9.0=pyhd8ed1ab_1
xerces-c=3.2.5=he0c23c2_2
xorg-libxau=1.0.11=hcd874cb_0
xorg-libxdmcp=1.1.3=hcd874cb_0
xyzservices=2024.9.0=pyhd8ed1ab_0
xz=5.2.6=h8d14728_0
yaml=0.2.5=h8ffe710_2
zeromq=4.3.5=ha9f60a1_6
zict=3.0.0=pyhd8ed1ab_0
zipp=3.20.2=pyhd8ed1ab_0
zlib=1.3.1=h2466b09_2
zstandard=0.23.0=py312h7606c53_1
zstd=1.5.6=h0ea2cb4_0 prefix: C:\Users\olive\miniconda3\envs\atlite-test

Oct 20 '24 21:10 oruhnau