本文介绍如何利用 NumPy 的 np.interp 对 Pandas DataFrame 中的离散数据点进行快速、准确的线性插值查询，无需修改原数据，支持边界值定制，适用于科学计算与工程建模场景。
本文介绍如何利用 numpy 的 `np.interp` 对 pandas dataframe 中的离散数据点进行快速、准确的线性插值查询，无需修改原数据，支持边界值定制，适用于科学计算与工程建模场景。
在数据分析中，我们常遇到已知若干离散 (X, Y) 坐标点，需查询任意 X 值（尤其是非整数或未采样点）对应 Y 的线性插值结果。Pandas 自带的 .interpolate() 方法主要用于填充缺失值且会修改原 DataFrame 或返回新结构，不适用于即查即用的插值查询场景；而 np.interp 正是为此类一维线性插值量身打造的轻量级、高性能工具。
以下为完整实践示例：
import pandas as pd
import numpy as np
# 构建示例数据
d = {'X': [1, 2, 3], 'Y': [220, 187, 170]}
df = pd.DataFrame(data=d)
# 查询 X = 1.5 对应的线性插值 Y 值
x_query = 1.5
y_interp = np.interp(x_query, df['X'], df['Y'])
print(y_interp) # 输出：203.5
✅ 关键前提：df['X'] 必须严格单调递增（或递减），否则 np.interp 行为未定义。若原始数据无序，务必先排序：
df_sorted = df.sort_values('X').reset_index(drop=True)
y_interp = np.interp(1.5, df_sorted['X'], df_sorted['Y'])
? 边界处理灵活可控：
- 默认情况下，超出 X 范围的查询值（如 x=0.5 或 x=4）将被“钳位”至最近端点值（即 df['Y'].iloc[0] 或 df['Y'].iloc[-1]）；
- 可通过 left 和 right 参数自定义越界行为：
# X max(df['X']) → 返回 -999 np.interp([0.5, 3.5], df['X'], df['Y'], left=np.nan, right=-999) # array([nan, -999.])
? 批量查询更高效：np.interp 支持数组输入，一次计算多个插值点：
x_queries = [1.2, 1.8, 2.5] y_results = np.interp(x_queries, df['X'], df['Y']) # array([213.4, 198.6, 178.5])
⚠️ 注意事项：
- np.interp 仅支持一维线性插值，不支持高阶（如样条）插值；如需更复杂拟合，请考虑 scipy.interpolate.interp1d；
- 确保 df['X'] 和 df['Y'] 长度一致且无 NaN；如有缺失，需提前清洗；
- 该方法纯数值计算，不依赖 Pandas 索引，因此性能远优于基于布尔索引+手动公式推导的实现。
综上，np.interp 是 Pandas 生态中实现即时、可复用、零副作用线性插值查询的最佳实践方案——简洁、可靠、高效。
