散点场重建

Wendland 紧支撑核的散点场重建方法 Wendland Functions（这篇笔记对 Wendland 核函数的数学构造有更完整的推导） 从散点到连续场 Wendland 并不是“力场重建”专用算法。它更像是一种通用工具：当你手里有一批离散采样点，想把它们变成一个连续场时，Wendland C2 紧支撑核可以作为一个很好的重建核。 问题可以抽象成这样：平面上有一组采样点 $$ (x_i, y_i, v_i), \quad i = 1,2,\dots,N $$其中 $(x_i,y_i)$ 是采样位置，$v_i$ 是这个位置的观测值。我们想得到一个连续函数 $f(x,y)$，让它在没有传感器的位置也能给出合理估计。压力传感器数据只是这个问题的一个例子。换成温度站点、高程采样点、点云属性，甚至仿真里的 collocation points，本质都是“散点样本 → 连续场”。 这里有一个容易混淆的点：这件事到底是“重建”，还是只是“平滑”？我的理解是，平滑更像后期修图，它不关心背后有没有真实函数，只是让数据看起来更顺；重建则先假设背后确实存在一个连续函数，而采样点只是对它的离散观测。对于压力传感器来说，我们认为背后存在一个连续压力分布 $P(x,y)$，传感器只是采到了若干个 $P(x_i,y_i)$，所以目标是从这些离散值恢复连续压力分布的近似，而不只是让图变得柔和。 这也是为什么我更愿意把这个问题叫“散点场重建”。它不是某个物理量专属的问题，而是一类数据重建问题。 为什么会选到 Wendland 从散点构造连续场的方法很多。最朴素的是最近邻：每个位置直接拿最近采样点的值。这很诚实，但画出来会像马赛克。线性插值或三角插值更自然一些，不过它往往需要明确的拓扑结构，边界处理也麻烦。点高斯叠加看起来很直觉：每个采样点盖一个光斑，所有光斑相加。但它有一个致命问题——如果不归一化，采样点密集的地方会天然更亮。亮度来自采样密度，而不是来自真实场值。 KNN-IDW 往前走了一步：找附近 K 个点，按距离倒数加权平均。它已经是“插值”的思路了，但邻居集合会随着像素位置切换，梯度大的地方容易出现细小斑点。局部平面拟合能表达局部趋势，不过不同邻域拟合出的平面一切换，就容易出现片状接缝。全局 RBF 插值理论上更完整，但需要解线性方程组，计算成本会比较高。 所以这个场景真正需要的不是“最花哨”的方法，而是同时满足几个工程约束的方法。第一，它要归一化，否则采样点密度会污染颜色。第二，它要局部计算，否则输出图像一大，每个像素都遍历所有采样点就不现实。第三，它要视觉连续，不能有块状边界、凸包、片状拼缝。Wendland C2 恰好把这几件事放在一起：它是归一化加权平均；它是紧支撑核，超出半径就严格为零；它在支撑边界处 C² 连续，不容易产生拼缝；同时它还是正定 RBF，作为插值或近似核有数学基础。 也就是说，Wendland 不是因为“它专门适合压力场”才被选中，而是因为当前任务具备这些条件：散点采样、需要连续场、采样密度可能不均匀、要求实时、要求视觉上没有明显伪影。 从旧的高斯光斑到 Wendland 核 旧的点高斯方案很容易理解：在每个采样点位置盖一个高斯光斑，然后把所有光斑相加。 $$ \text{field}{[x][y]} = \sum_i \text{value}_i \times \exp\!\left(-\frac{d_i^2}{2\sigma^2}\right) $$这里 $\text{value}_i$ 是第 $i$ 个采样点的值，$d_i$ 是像素到这个采样点的距离，$\sigma$ 控制光斑扩散范围。这个公式最大的问题就在求和符号上。两个采样点靠得近，光斑重叠区域就会叠得更亮。结果是：采样点密度越高的地方越亮，即便真实场值并没有更高。 ...