Support subspacing to nearest neighbour of provided value
還沒有人認領這個 Issue。
評估
研究方向
首先閱讀 subspace 和 indices 方法,以及容差內容管理器 cf.atol 和 cf.rtol。將所要求的最近鄰行為與連結的 xarray 文件進行比較。完成的標準是為這些方法支援 nearest 模式,並新增 nearest 和基於容差的子空間劃分文件範例。
由索引模型根據 Issue 內容生成。
描述
A feature request for new functionality for the subspace and indices methods. (xarray supports this functionality for subspacing but we don't, however it seems too useful not to want to steal!)
Picture the scene. You want to do a subspace on a specific data value for some coordinate but due to natural complexity of real-life data, it is a float with many digits to define, as it often is. At the moment in cf-python you have to specify that value exactly e.g. f.subspace(grid_latitude=7.480000078678131) to subspace down to the grid_latitude of 7.480000078678131 (assuming you want to use subspacing by metadata, let's assume you do, you could use indexing but that involves knowing/calculating the appropriate index). But it would be nice to be able to also request by some means 'the value nearest to' a specified value, in this case provide only 7.48 for the grid_latitude to avoid having to put in the exact float, just a close-enough approximation.
xarray supports such 'nearest neighbour' lookups, where you can specify a 'nearest' neighbour method as a method keyword, e.g. 7.48 to go to 7.480000078678131 as the nearest neighbour, plus a tolerance on inexact look-up e.g. 7.5 with +/- 0.1 tolerance to catch this same value that way.
I would like us to support directly the nearest neightbour match, and better advertise how to do the inexact subspace in two lines using our tolerance functions in a context manager.
We already have a new-ish 'halo' subspacing approach and are going to add a 'bounding box' query too. It would be good to consider how to include the above 'nearest neighbour' to the new possibilities for more flexible subspacing. I have made some suggestions below to get the conversation started.
Example
Example set up of a subspaces we'd like to provide the above functionality to simplify:
>>> print(f)
Field: relative_humidity (ncvar%UM_m01s16i204_vn405)
----------------------------------------------------
Data : relative_humidity(air_pressure(17), grid_latitude(30), grid_longitude(24)) %
Cell methods : time(1): mean
Dimension coords: time(1) = [1978-12-16 12:00:00] gregorian
: air_pressure(17) = [1000.0000610351562, ..., 10.0] hPa
: grid_latitude(30) = [7.480000078678131, ..., -5.279999852180481] degrees
: grid_longitude(24) = [-5.720003664493561, ..., 4.399996280670166] degrees
Auxiliary coords: latitude(grid_latitude(30), grid_longitude(24)) = [[61.004354306111864, ..., 48.51422609871432]] degrees_north
: longitude(grid_latitude(30), grid_longitude(24)) = [[-13.762685427418687, ..., 4.622216504491947]] degrees_east
Coord references: grid_mapping_name:rotated_latitude_longitude
>>> f1 = f.subspace(air_pressure=1000.0000610351562)
>>> f2 = f.subspace(grid_latitude=7.480000078678131)
Suggestion for API to provide support
Add a new 'mode' with string identifier 'nearest' to do the nearest neighbour case, e.g. for the above example this would work:
f.subspace(air_pressure=1000, mode="nearest")
The latter case is controlled by our tolerance functions cf.atol and cf.rtol but there are no examples in the documentation to advertise how simply one can control the subspacing tolerance with:
with cf.atol(1e-2):
f.subspace(grid_latitude=7.480000078678131)
so I'd also like us to showcase examples of using a context manager like above as a two-line means to easily do inexact subspacing.
- 主要語言
- Python
- 星號
- 151
- 分支
- 23
- 平均合併
- 1 天 2 小時
- 30 天內合併 PR
- 1
環境準備
- 沒有 Dockerfile 或 Docker Compose 檔案
- 沒有 Pull Request 範本
- 閱讀貢獻指南
從這裡開始
- 先讀完整個 Issue,再讀專案的貢獻指南。
- 在 Issue 下留言說明你要接手 —— 這能避免兩個人做同樣的事。
- Fork 儲存庫,在一個分支上完成修改。
- 送出 Pull Request,並在描述裡引用這個 Issue 編號。
NCAS-CMS/cf-python 的其他 Issue
-
installation
難度 4/5 3-5 天 新手友好度 45/100
-
Attempt at unsupported conservative 3D mesh regridding可能重新可做 @sadielbartholomew 於 40 天前認領,目前沒有進行中的 PR。 未關閉bug?
-
bug
難度 3/5 1-2 天 新手友好度 45/100
-
API review (4.0.0) documentation
難度 4/5 3-5 天 新手友好度 48/100
-
enhancement xarray
難度 4/5 3-5 天 新手友好度 48/100
查看 NCAS-CMS/cf-python 的全部 Issue
相似的 Issue
-
難度 2/5 1-3 小時 新手友好度 76/100
-
[Device Support Request] HOBEIAN ZG-204ZL (Tuya PIR motion sensor with illuminance)可能已有人在做 @lenucksi 今天認領。 未關閉
難度 2/5 1-3 小時 新手友好度 78/100
zigpy/zha-device-handlers#5493 ·
維護者通常 3 天內回覆
-
難度 2/5 1-3 小時 新手友好度 78/100
drevendev/BeautifulRepo#42 ·
維護者通常 1 天內回覆
-
enhancement good first issue
難度 2/5 1-3 小時 新手友好度 84/100
writerslogic/holographic-memory#93 ·
維護者通常 1 天內回覆
-
難度 2/5 1-3 小時 新手友好度 66/100
antoine-tena/kiln#47 ·
維護者通常 1 天內回覆