极坐标图中指定角度未知点的最优插值方法

别上来直接在极坐标域拿角度当横轴、半径当纵轴做插值，这是最常见的坑——极坐标下的等角度线性插值，转成直角坐标看曲线会往原点方向畸变，点位越稀疏误差越大，甚至还不如你现在直接连相邻点的直线靠谱。
按你的精度需求，按优先级选方法就行：

首选方案：弦长参数化的笛卡尔域三次样条插值 #

这是兼顾精度、光滑度和实现难度的最优解，0-180°区间刚好对应直角坐标系y≥0的上半平面，没有跨周期多值的问题，实现步骤非常固定：

• 先把所有已知极坐标点(θ_i, r_i)统一转成直角坐标，计算前记得把角度转成弧度制：

  x_i = r_i * cos(θ_i)
  y_i = r_i * sin(θ_i)
  

• 不要直接用角度作为插值的自变量，要以相邻点的累计弦长作为插值参数t：第一个点t值取0，后续每个点的t值等于前一个点的t值加上两点间的欧氏距离，最后把所有t值归一化到[0,1]区间即可，这种参数化方式不会出现曲线自交、过冲的问题。
• 分别对x关于t、y关于t做自然三次样条拟合，得到两个连续的插值函数x(t)、y(t)。
• 给定待求角度θ_target时，先通过小范围迭代找到满足arctan2(y(t), x(t)) = θ_target的t值，代入插值函数算出对应的(x,y)，再转回极坐标就能得到对应的半径值，0-180°区间函数单调，迭代2-3次就能收敛到精度要求。
  这个方法插出来的曲线二阶连续光滑，没有折线的棱角，只要你的原始点位采样密度不是特别低，插值误差比极坐标域线性插值小一个量级以上，我自己在雷达点云插值、极坐标绘图的场景用了很多年，稳定性足够。

分场景备选方案 #

• 如果你的原始点位是实测数据、带测量噪声：别用普通三次样条，容易过拟合噪声，换成带平滑系数的薄板样条，或者转直角坐标后做局部加权回归（LOWESS），平滑系数根据你的噪声水平调整，既能压掉噪声，也比直线连接的精度高。
• 如果你的原始点位是严格等角度间隔采样、且径向没有剧烈突变：可以偷懒在极坐标域用保形分段三次插值（PCHIP），比直接连直线光滑，也不会出现普通三次样条的边缘过冲问题，但注意这个方法只适合等角采样的场景，不等间隔采样时误差会明显升高。
• 绝对避坑：不要在极坐标域直接做全局高次多项式插值，点位数量超过7-8个就会出龙格现象，靠近0°和180°边界的位置半径会飘得非常离谱，效果比连直线还差。

补充提示：如果你后续要把插值范围扩展到0-360°全周，只需要把三次样条的边界条件改成周期边界，保证首尾点的一阶、二阶导数连续就行，你现在用的0-180°区间不需要做这个处理，自然样条的默认边界条件就够用。

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