developmentseed / developmentseed/titiler-multidim

Issue Visualizing Daymet Zarr File with Projected and Chunked Dimensions

Open
#67 1 comment 0 reactions 0 assignees View on GitHub
Dominant language
Python
Stars
88
Forks
13
Avg merge
2d 53m
Merged PRs (30d)
10

Description

Hi,

I’m trying to visualize a Daymet Zarr file using titiler-xarray, but I’m encountering issues related to chunked dimensions and projected coordinates. Below is a description of the two datasets I’ve been working with and the specific error I’m encountering.

Dataset 1: NOAA CoastWatch Geo-Polar SST (Successfully Visualized)

I was able to visualize the Zarr dataset from the following URL:

```
zarr: https://ncsa.osn.xsede.org/Pangeo/pangeo-forge/noaa-coastwatch-geopolar-sst-feedstock/noaa-coastwatch-geopolar-sst.zarr
```
Metadata:
```
Size: 3TB
Dimensions: (time: 7134, lat: 3600, lon: 7200)
Coordinates:
* lat (lat) float32 -89.97 -89.93 -89.88 ... 89.93 89.97
* lon (lon) float32 -180.0 -179.9 -179.9 ... 179.9 180.0
* time (time) datetime64[ns] 2002-09-01T12:00:00 ... 2022...
Data variables:
analysed_sst (time, lat, lon) float32 dask.array
analysis_error (time, lat, lon) float32 dask.array
mask (time, lat, lon) float32 dask.array
sea_ice_fraction (time, lat, lon) float32 dask.array
```

Dataset 2: Daymet Zarr (Issue Encountered)

When trying to visualize the Daymet Zarr dataset from the following URL:
```
zarr: https://daymeteuwest.blob.core.windows.net/daymet-zarr/daily/hi.zarr
```
Metadata:
```
Size: 69GB
Dimensions: (time: 14965, y: 584, x: 284, nv: 2)
Coordinates:
lat (y, x) float32 dask.array
lon (y, x) float32 dask.array
* time (time) datetime64[ns] 1980-01-01T12:00:00 ...
* x (x) float32 -5.802e+06 ... -5.519e+06
* y (y) float32 -3.9e+04 -4e+04 ... -6.22e+05
Dimensions without coordinates: nv
Data variables:
dayl (time, y, x) float32 dask.array
prcp (time, y, x) float32 dask.array
tmax (time, y, x) float32 dask.array
tmin (time, y, x) float32 dask.array
```

Issue:
I received the following error when trying to visualize the Daymet dataset:
```
ValueError: conflicting sizes for dimension 'y': length 256 on the data but length 105 on coordinate 'lat'
```

I noticed that the Daymet Zarr dimensions are chunked, and the lat/lon coordinates are projected to x, y dimensions. I suspect this is the reason why titiler-xarray is unable to visualize the dataset. I’ve also encountered similar issues while attempting to visualize other NetCDF files.

Question:

Does titiler-xarray currently support visualizing datasets with chunked and projected dimensions like the Daymet dataset? If not, I would be happy to contribute to implementing this feature. Could you please provide some guidance or suggestions on how I can modify titiler-xarray or my data to work around this issue?

Contributor guide

No contributing guide indexed for this repository

Assessment

This issue has not been assessed yet.

Get new issues in your inbox

A short digest of beginner-friendly GitHub issues.