R相當於matlab griddata，scatteredInterpolant和/或TriScatteredInterp

Question

我們做了很多全場3D數值模擬（CFD，FEA等）。解決方案需要很長時間才能運行。我們通常從解決方案插入而不是重新運行每個案例。我們還在多個解之間進行插值，這導致更高維插值（如加入時間，所以x，y，z，t，v）。

Matlab在X，Y，Z座標的不規則網格中讀取數據V，使用griddata，scatterdInterpolan和/或TriScatteredInterp從V進行插值的功能非常出色。由於各種各樣的原因，我已經轉向R.這仍然是我無法找到的與R相當的關鍵領域。 'akima'只有x，y，V（不是，x，y，z，V，更不像x，y，z，t，v那樣更高的尺寸）。

我發現的下一個最好的東西是'krigging'。但是krigging的行爲更像是模型擬合和投影，並且在不規則的網格點之間往往表現不佳。所以它不如簡單的直接線性插值那樣穩健。

Matlab有幾十年的griddata。很難相信R在那裏沒有相應的東西。有什麼建議麼？或者至少有一種方法可以使用krigging來有效地產生與直接線性插值相同的結果？

喬納森

Answer 1

你可以通過看包「tripack」做Delaunay三角，讓您在複製scatteredInterpolant）的第一步（開始。