F90从非恒定分辨率3D网格到恒定分辨率网格的插值

问题描述

我在分辨率不同的3D笛卡尔网格上有数据。网格仍然是结构化的，因为x分辨率只是x的函数，与y，z相同。我需要将这些数据插值到大小大致相同的恒定分辨率网格上。典型大小为256 ^ 3、512 ^ 3。

由于噪声数据，三线性插值法不会将其剪切掉，但会留下更高的阶数（出于准确性或平滑度，我还没有决定哪种方法最好）。我知道该软件 Mongo Playground，但这说明它可以在常规网格上运行。如果“常规”表示恒定的分辨率，那么这可能无济于事。但是，我可能误解了“常规”的含义。

有人知道可以用于我的问题的开源fortran软件吗？即从一个大的3d网格（非恒定分辨率）到一个新的3D网格（恒定分辨率）进行插值，这不仅仅是三线性的。

谢谢

解决方法

术语“常规”网格的使用频率不如“统一网格”（不是网格）和“结构化网格”（实际上是）。 bspline-fortran的规范似乎表明插值适用于结构化，直线形，拉伸网格，因此我想这是您情况下“常规”网格的定义。好消息是，根据该定义，您的网格是“常规的”。

对于您的情况，重要的问题是您不均匀的网格是否是“平滑变换函数”的结果。如果是这样，您的网格可以被视为“变换空间”中的统一网格，并且可以使用可用的3D插值。

我的建议是尝试一下，找到一些评估准确性的方法。