编辑-地球科学江南体育网页版堆栈交换江南电子竞技平台

你不是登录。你的编辑将被放置在一个队列,直到同行评议。

我们欢迎编辑,使读者更容易理解和更有价值。因为社区成员审查编辑,请尽量使文章更好比你怎么发现,例如,通过修正语法或添加额外的资源和超链接。

如何寻找最近的非nan值在2 d xarray数据集

我有一个2 d xarray数据集和我有几个nan值的数据集。这些nan值需要充满最近的non-nan值的二维数组。我怎么做呢?

我试图让我的数据符合grid_spec文件读取的地球系统模型。的土地上有矛盾点是在迫使数据和grid_spec文件导致的爆炸模型。

回答

我能够通过以下功能def fill_mismatch (ds_grid_spec ds_data): ds_data_o = ds_data.copy()数= 0.0我的范围(ds_grid_spec.shape [0]): j的范围(ds_grid_spec.shape[1]):如果ds_grid_spec (I, j)。值> 0.0,(不是ds_data (i, j) > 0.0):数=数+ 1.0 fill_miss2d (ds_data_o, i, j)如果ds_grid_spec (i, j)。值= = 0.0,(不是ds_data (i, j)。值= = 0.0):数=数+ 1.0 ds_data_o (i, j) = 0.0返回ds_data_o def fill_miss2d (ds, lat_i lon_j): #拳头循环周期= 1时真的:bl_i, bl_j = lat_i -周期,lon_j——周期tr_i tr_j = lat_i +循环,lon_j +周期#左下角我如果bl_i < 0: bl_i = 0 #左下角j如果bl_j < 0: bl_j = 0 #右上方我如果tr_i > ds.lat。形状[0]1:tr_i = ds.lat。形状[0]1 #如果tr_j > ds.lon右上方j。形状[0]1:tr_j = ds.lon。形状[0]1 lats_i = np.arange (bl_i, tr_i) lons_j = np.arange lats_i (lats_i_idx bl_j tr_j):在lons_j lons_j_idx:如果ds (lats_i_idx, lons_j_idx) > 0.0: ds [lat_i lon_j] = ds (lats_i_idx lons_j_idx]回报周期=周期+ 1 grid_spec(这是最常用的网格文件交换模块化海洋模型妈妈)有一个变量AREA_LND soiltype和文件。数控变量soiltype。这两个变量有0.0网格点归类为海洋。因此我们不需要使用fill_miss2d发生失配和海洋在哪里遇到grid_spec,但是土地是在soiltype。而我们直接填补soiltype网格点为0.0完整的笔记本可以在https://github.com/manmeet3591/python_class/blob/master/xarray_tutorial/grid_soiltype.ipynb上找到

编辑总结

取消

添加一个评论 |

正确的小错误或错误
澄清的意思,而不更改它
添加相关资源或链接
总是尊重作者的意图
不使用编辑与作者回复

创建代码用引号”或式子~篱笆
' ' '
像这样
' ' '
添加语言标识符来突出显示的代码
”“python
def函数(foo):
打印(foo)
' ' '
把段落之间的回报
在结束linebreak添加2个空间
_italic_或* *的* *
报价将>开始的
做链接(尽量使用https)
< https://example.com >
(例子)(https://example.com)
< a href = " https://example.com " > < / >示例
MathJax方程罪\ \ ^ 2θ美元

格式帮助»
回答帮助»

MathJax帮助»