从具有空白空间的卷中插入散乱数据

时间:2015-09-02 15:52:25

标签: python numpy interpolation

我有从网格点生成的3d数据。网格结构足够复杂,缺乏使用griddata的插值。具体而言,存在没有数据点的区域,这些数据点由griddata给出的值不是fill_value。我需要这些空心区域的值为0.0,我将fill_value设置为。

简化版本如下图所示: cylinider intersecting data region

圆柱体占据的区域没有数据点,但是立方体体积的其余部分确实存在。圆柱内的插值会有数据点,但我需要它们为零。

下面是与实际内插数据的xy平面平行的切片,其中黑色椭圆近似于边缘'。红色和蓝色'流血'插值后的空白。填充值0.0可以在左上角看到:

enter image description here

关于如何实现将这些值设置为0.0的目标的任何想法?请注意,'圆筒'并不是一成不变的。

我想到了z层的z层,找到了一个给出圆柱形状的多边形,然后将多边形内的点设置为零。

我还考虑过对音量进行分区,以便圆柱体的一部分最终位于部分的角落(对于每个z层),并希望插值器不会尝试外推到空白区域。

第一个选项似乎更好,但我想知道Python是否提供了一些更好的功能。

编辑:以下是数据集中的一些实际点:

enter image description here

z比例远小于x或y。你可以看到我感兴趣的地区定义得很好。但是,再次,我如何识别它们以便将网格点设置为0.0?

0 个答案:

没有答案