請問scipy.interpolate.griddata是否正確地爲經度和緯度浮點座標插值？

Question

可以說，我有3個np.array在Python中的float值：latitudes,longitudes和values。我想插入一些不在latitude和longitude陣列中的座標值。我可以使用scipy.interpolate.griddata，它工作得很好。

但是我最近開始懷疑interpolate是否按我預期的那樣工作：因爲它不知道歐幾里德網格和緯度和經度網格的球形特性之間的差別，所以它可能計算錯誤的插值至少爲linear。 （因爲2座標之間的距離在兩極也許比在赤道兩個座標之間的距離大大不同。）

我確實需要看，可以發現 at the scipy site但沒有發現有關的處理信息的文檔不等距離的轉換。

座標的球面座標和等距座標之間的差異是scipy.interpolate.griddata的問題，還是不相關？

Answer 1

griddata假設歐幾里得距離測量，所以你有通常的球面映射投影問題。

SciPy的0.11.0具有工具球體上的數據進行插值：http://docs.scipy.org/doc/scipy/reference/generated/scipy.interpolate.RectSphereBivariateSpline.html

這適用於直線緯度/經度網格。

Answer 2

您需要的是Slerp (spherical linear interpolation)，這是一種使用quaternion在球體上線性插值的方法。有一個quaternion dtype for numpy作爲一個附加組件實現，它可以讓你有效地在python中實現slerp。 Here也是一個實現slerp的庫。