Python中的4D插值

Question

因此，我在x y和z有一個4d力場，每個點都附加了一個Force 。

我對如何插值有點迷茫，盡管我遇到了這個問題。

http://docs.scipy.org/doc/scipy-dev/reference/generated/scipy.ndimage.interpolation.map_coordinates.html

看來我可以制作一個xyz F數組，然后對其進行插值，但我確實不太確定。

任何人都知道如何使用python的任何好的庫？

輸入示例為：

x y z Force  
0 0 0 35  
0 0 1 32  
0 0 2 23  
...  
2 5 9 54  

Answer 1

您描述它的方式（x，y，z），這是3D字段，而不是4D。

插值的目的是在任意點（x，y，z）上都指定一個字段（在本例中為Force ），即使您沒有力也恰好在該點（x，y，z）上。 取而代之的是，您在多個附近的點處具有力，並對它們進行插值以獲得在（x，y，z）點處的有意義的力。 在3D中，您至少需要一個雙線性插值算法，並且對其進行編碼（並使索引保持順序）並非易事。

在Python中，您可以使用Scipy的插值例程 。

Answer 2

因此，如果我了解您的問題，那么您在某些點上施加了力，而在另一些點上卻沒有，並且想提出插值法？ 在NMR中，線性插值的使用一直是最新的。 作為替代方案，您可以嘗試樣條曲線，但是那會花費更多。 但這主要取決於您的力量類型-變化的速度有多快？ 在此基礎上構建插值路由。