issues: 498180644
This data as json
| id | node_id | number | title | user | state | locked | assignee | milestone | comments | created_at | updated_at | closed_at | author_association | active_lock_reason | draft | pull_request | body | reactions | performed_via_github_app | state_reason | repo | type |
|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
| 498180644 | MDU6SXNzdWU0OTgxODA2NDQ= | 3342 | Error when using xarray's built in interp function with dask chunking | 20382006 | closed | 0 | 4 | 2019-09-25T10:00:13Z | 2019-09-26T13:56:19Z | 2019-09-26T13:56:19Z | NONE | This is a copied issue posted on stack overflow: https://stackoverflow.com/questions/57572433/error-when-using-xarrays-built-in-interp-function-with-dask I am trying to use xarray's interp function together with the chunk function, which is a dask array feature. Taking the example under 'Example' at http://xarray.pydata.org/en/stable/interpolation.html I am able to use the interp function just fine. The simplified code I use to do this is: ``` at = xr.tutorial.open_dataset('air_temperature').isel(time=0) x = np.linspace(240, 300, 100) z = np.linspace(20, 70, 100) lat = xr.DataArray(z, dims=['z'], coords={'z': z}) lon = xr.DataArray((x[:, np.newaxis]-270)/np.cos(z*np.pi/180)+270, dims=['x', 'z'], coords={'x': x, 'z': z}) dsi = at.interp(lon=lon, lat=lat)` ``` The problem arises when I try to combine this interp function with xarray chunk by modifying the above code in the following way: ``` at = xr.tutorial.open_dataset('air_temperature').isel(time=0) at = at.chunk({'lat':10}) # added chunking inducing a dask array x = np.linspace(240, 300, 100) z = np.linspace(20, 70, 100) lat = xr.DataArray(z, dims=['z'], coords={'z': z}) lon = xr.DataArray((x[:, np.newaxis]-270)/np.cos(z*np.pi/180)+270, dims=['x', 'z'], coords={'x': x, 'z': z}) dsi = at.interp(lon=lon, lat=lat) ``` The code fails at line
ValueError: Can't specify drop_axis and new_axis together I think what is happening is that the interp function is replacing a coordinate and dask cannot deal with this. Is this a bug that needs fixing within xarray or is there something that I have done wrong? Output of
|
{
"url": "https://api.github.com/repos/pydata/xarray/issues/3342/reactions",
"total_count": 0,
"+1": 0,
"-1": 0,
"laugh": 0,
"hooray": 0,
"confused": 0,
"heart": 0,
"rocket": 0,
"eyes": 0
} |
completed | 13221727 | issue |