如何使用scipy.interpolate实现按点顺序的插值曲线？

你的问题根源在于scipy.interpolate.interp1d的核心限制：它要求输入的x轴数据必须是单调递增或递减的。而你提供的x数组[1,2,3,4,3,2]是先升后降的非单调序列，这会让interp1d无法正确识别点的连接顺序，最终生成了不符合预期的折线。

要实现你想要的沿点顺序的平滑曲线，我们需要用参数化插值的方式：把点的索引（或一个均匀分布的参数）作为自变量，分别对x和y坐标进行插值，这样就能沿着点的连接路径生成平滑曲线了。这里可以用scipy.interpolate.splprep和splev来实现：

修改后的代码 #

import matplotlib.pyplot as plt
import numpy as np
from scipy import interpolate

# 原始数据点
x = np.array([1,2,3,4,3,2])
y = np.array([1,1,1,1,2,2])

# 创建参数t：用0到1之间的均匀值代表每个点的顺序
t = np.linspace(0, 1, len(x))

# 对x和y进行样条插值，s=0表示严格经过所有原始点
tck, u = interpolate.splprep([x, y], s=0)

# 生成更密集的参数点
t_new = np.linspace(0, 1, 1000)

# 用插值后的样条曲线计算新的x和y坐标
x_new, y_new = interpolate.splev(t_new, tck)

# 绘图
plt.plot(x, y, 'o', label='原始点')
plt.plot(x_new, y_new, '-', label='插值曲线')
plt.legend()
plt.show()

代码解释 #

• splprep函数用于对多变量（这里是x和y）进行参数化样条插值，s=0确保插值曲线严格经过所有原始数据点；
• 我们用参数t来代表点的顺序，而不是直接用x作为自变量，这样就避开了x非单调的问题；
• splev函数根据生成的密集参数点，计算出插值后的x和y坐标，最终得到平滑的曲线。

运行这段代码后，你就能得到类似你期望的那种沿点顺序平滑过渡的曲线：从(1,1)依次经过(2,1)、(3,1)、(4,1)，然后平滑转向(3,2)、(2,2)的连续曲线。

内容的提问来源于stack exchange，提问作者user8267234