html_url,issue_url,id,node_id,user,created_at,updated_at,author_association,body,reactions,performed_via_github_app,issue
https://github.com/pydata/xarray/issues/475#issuecomment-355085766,https://api.github.com/repos/pydata/xarray/issues/475,355085766,MDEyOklzc3VlQ29tbWVudDM1NTA4NTc2Ng==,2443309,2018-01-03T18:18:32Z,2018-01-03T18:18:32Z,MEMBER,http://notes.stefanomattia.net/2017/12/12/The-quest-to-find-the-closest-ground-pixel/,"{""total_count"": 0, ""+1"": 0, ""-1"": 0, ""laugh"": 0, ""hooray"": 0, ""confused"": 0, ""heart"": 0, ""rocket"": 0, ""eyes"": 0}",,95114700
https://github.com/pydata/xarray/issues/475#issuecomment-355084829,https://api.github.com/repos/pydata/xarray/issues/475,355084829,MDEyOklzc3VlQ29tbWVudDM1NTA4NDgyOQ==,2443309,2018-01-03T18:14:51Z,2018-01-03T18:14:51Z,MEMBER,"@stefanomattia - I'd be happy to provide guidance and even to contribute to some of the development. Based on your blog post, I think you may be well on your way. ","{""total_count"": 0, ""+1"": 0, ""-1"": 0, ""laugh"": 0, ""hooray"": 0, ""confused"": 0, ""heart"": 0, ""rocket"": 0, ""eyes"": 0}",,95114700
https://github.com/pydata/xarray/issues/475#issuecomment-354703904,https://api.github.com/repos/pydata/xarray/issues/475,354703904,MDEyOklzc3VlQ29tbWVudDM1NDcwMzkwNA==,2443309,2018-01-02T05:04:08Z,2018-01-02T05:04:08Z,MEMBER,ping @stefanomattia who seems to be interested in the KDTreeIndex concepts described in this issue.,"{""total_count"": 0, ""+1"": 0, ""-1"": 0, ""laugh"": 0, ""hooray"": 0, ""confused"": 0, ""heart"": 0, ""rocket"": 0, ""eyes"": 0}",,95114700
https://github.com/pydata/xarray/issues/475#issuecomment-342558818,https://api.github.com/repos/pydata/xarray/issues/475,342558818,MDEyOklzc3VlQ29tbWVudDM0MjU1ODgxOA==,2443309,2017-11-07T17:28:17Z,2017-11-07T17:28:17Z,MEMBER,"@WeatherGod 

Short answer. We don't have a tool that is production ready. 

Longer answer: `xESMF` may be the best prospect in the near term. There are two main issues with its current implementation. 1) Lack of out-of-core abilities / integration with dask, and 2) lack of a test suite. Conceptually, it would be great to leverage the low-level remapping tools of ESMPy so I think this is a nice way to move forward as a community but I think everyone agrees it isn't ready for use in any sort of production environment.

This issue introduces the concept of point-wise indexing using nearest neighbor lookups on ND coordinates. @shoyer has an example implementation [here](https://github.com/pydata/xarray/issues/475#issuecomment-125468579) but it hasn't moved forward in quite a while.","{""total_count"": 0, ""+1"": 0, ""-1"": 0, ""laugh"": 0, ""hooray"": 0, ""confused"": 0, ""heart"": 0, ""rocket"": 0, ""eyes"": 0}",,95114700
https://github.com/pydata/xarray/issues/475#issuecomment-125849716,https://api.github.com/repos/pydata/xarray/issues/475,125849716,MDEyOklzc3VlQ29tbWVudDEyNTg0OTcxNg==,2443309,2015-07-29T05:44:35Z,2015-07-29T05:44:35Z,MEMBER,"Very nice.  This is the sort of API I was hoping for.  It will be a while before I can come back around on this. In the meantime, if someone else wants to take the `sel_points` method on, that is fine by me.
","{""total_count"": 0, ""+1"": 0, ""-1"": 0, ""laugh"": 0, ""hooray"": 0, ""confused"": 0, ""heart"": 0, ""rocket"": 0, ""eyes"": 0}",,95114700
https://github.com/pydata/xarray/issues/475#issuecomment-125334057,https://api.github.com/repos/pydata/xarray/issues/475,125334057,MDEyOklzc3VlQ29tbWVudDEyNTMzNDA1Nw==,2443309,2015-07-27T20:31:03Z,2015-07-27T20:31:03Z,MEMBER,"Now that the `isel_points` method is implemented, I think it makes sense to discuss the `sel_points` method in a bit more detail. The main outstanding question is - do we want to support spherical nearest neighbor mapping.  The use case is when you are searching for the nearest neighbor using longitudes and latitudes.  [This](http://earthpy.org/interpolation_between_grids_with_ckdtree.html) example shows an example of to do this by projecting the coordinates onto a sphere. If we go this route, which is probably the most common use case here, we are committing to the coordinates being latitudes and longitudes. Maybe it is better to use a `method='spherical'` keyword to fall into this path.
","{""total_count"": 0, ""+1"": 0, ""-1"": 0, ""laugh"": 0, ""hooray"": 0, ""confused"": 0, ""heart"": 0, ""rocket"": 0, ""eyes"": 0}",,95114700
https://github.com/pydata/xarray/issues/475#issuecomment-122198334,https://api.github.com/repos/pydata/xarray/issues/475,122198334,MDEyOklzc3VlQ29tbWVudDEyMjE5ODMzNA==,2443309,2015-07-17T06:54:37Z,2015-07-17T06:54:37Z,MEMBER,"Good point on the dask array business.  From the [dask docs](http://dask.pydata.org/en/latest/slicing.html):

> Dask.array supports most of the NumPy slicing syntax.
> ...
> It does not currently support the following:
> 
>    Slicing one dask.array with another x[x > 0]
>    **Slicing with lists in multiple axes x[[1, 2, 3], [3, 2, 1]]**
> 
> Both of these are straightforward to add though. If you have a use case then raise an issue.

So, from browsing the closed dask issues, it seems like dask has similar support for multi-dimension slicing and indexing as xray. This throws a bit of a wrench in my plan for how I was going to implement `isel_points` as I had not fully considered the dask array complexities.

I'll have to put a bit more thought into this. Any suggestions on how to index the dask array without looping through individual points would be great.
","{""total_count"": 0, ""+1"": 0, ""-1"": 0, ""laugh"": 0, ""hooray"": 0, ""confused"": 0, ""heart"": 0, ""rocket"": 0, ""eyes"": 0}",,95114700
https://github.com/pydata/xarray/issues/475#issuecomment-121998086,https://api.github.com/repos/pydata/xarray/issues/475,121998086,MDEyOklzc3VlQ29tbWVudDEyMTk5ODA4Ng==,2443309,2015-07-16T15:45:59Z,2015-07-16T15:45:59Z,MEMBER,"As a first step, I'll volunteer (unless someone else is more keen on doing this work) to put together a pull request for `isel_points`.

After that, we'll want to add the `sel_points` and `kdtree` API, which will depend on `isel_points`.

Later on, I'm also interested in regridding and resampling between grids - let's open another issue for that. Maybe we use `pyresample` for that.
","{""total_count"": 0, ""+1"": 0, ""-1"": 0, ""laugh"": 0, ""hooray"": 0, ""confused"": 0, ""heart"": 0, ""rocket"": 0, ""eyes"": 0}",,95114700
https://github.com/pydata/xarray/issues/475#issuecomment-121777990,https://api.github.com/repos/pydata/xarray/issues/475,121777990,MDEyOklzc3VlQ29tbWVudDEyMTc3Nzk5MA==,2443309,2015-07-15T23:51:14Z,2015-07-15T23:51:45Z,MEMBER,"I like:

``` python
DataArray.isel_points(x=[1, 2, 3], y=[0, 1, 2], dim='points')
```

I also like the nearest-neighbor / resample API of:

``` python
DataArray.sel_points(lon=[-123.25, -140.0, 72.5], lat=[45.0, 72.25, 65.75],
                     dim='points', method='nearest')
```

How do we want to do the nearest-neighbor selection?  The simplest case would be to follow the cKDTree example from #214.  However, when you're using lat/lon coordinates, it is usually best to map these coordinates from the spherical coordinates to a Cartesian coordinates (see [here](http://earthpy.org/interpolation_between_grids_with_ckdtree.html) for a simple example using cKDTree. Is that a road we want to go down here?

Further along that subject, but not directly relate - has anyone used [pyresample](http://pyresample.readthedocs.org/en/latest/index.html).
","{""total_count"": 0, ""+1"": 0, ""-1"": 0, ""laugh"": 0, ""hooray"": 0, ""confused"": 0, ""heart"": 0, ""rocket"": 0, ""eyes"": 0}",,95114700