我正在尝试根据lat / long中的输入GPS坐标创建连续路径。有没有python库可以做到这一点?最终目标是进行优化,因此我需要连续可微分的路径来确定梯度。
我一直在使用GDAL,但它似乎没有我需要的所有功能。
答案 0 :(得分:2)
这需要scipy.interpolate
包。 Interpolation in the SciPY Cookbook上的此页面显示了如何使用函数splprep
来确定样条曲线结,并使用splev
来绘制平面中噪声数据的结果曲线。
食谱示例(稍微缩短):
from numpy import arange, cos, linspace, pi, sin, random
from scipy.interpolate import splprep, splev
# make ascending spiral in 3-space
t=linspace(0,1.75*2*pi,100)
x = sin(t)
y = cos(t)
z = t
# add noise
x+= random.normal(scale=0.1, size=x.shape)
y+= random.normal(scale=0.1, size=y.shape)
z+= random.normal(scale=0.1, size=z.shape)
# spline parameters
s=3.0 # smoothness parameter
k=2 # spline order
nest=-1 # estimate of number of knots needed (-1 = maximal)
# find the knot points
tckp,u = splprep([x,y,z],s=s,k=k,nest=-1)
# evaluate spline, including interpolated points
xnew,ynew,znew = splev(linspace(0,1,400),tckp)
import pylab
pylab.subplot(1,1,1)
data,=pylab.plot(x,y,'bo-',label='data')
fit,=pylab.plot(xnew,ynew,'r-',label='fit')
pylab.legend()