文章目录
- 前言
- 一、分步拆解：从点到面的过程
-
- [1. 理解"栅格"](#1. 理解“栅格”)
- [2. 如何计算未知网格点的值？------样条的介入](#2. 如何计算未知网格点的值？——样条的介入)
- [3. 两种主要的样条栅格插值](#3. 两种主要的样条栅格插值)
- 二、为什么用样条做栅格插值？
- 总结
前言
样条函数法会生成一张整体光滑过渡的曲面。即使两个测点距离很远且值差异大，它也会用一条优雅的弧线连接起来，而不是生硬的斜坡或台阶。就像在测点间拉了一张光滑的弹性布。
一、分步拆解：从点到面的过程
1. 理解"栅格"
所谓"栅格"，就是像围棋棋盘或Excel表格一样，把整个区域划分成规则的行列网格（像素）。每个格子叫一个"像元"，我们的目标就是算出每个格子中心点的高程值。
问题来了：网格点上是没有实测数据的！怎么知道它的值？
2. 如何计算未知网格点的值？------样条的介入
这就是样条函数大显身手的地方。我们不是为整个区域找一个复杂函数，而是用样条的思路"局部平滑"：
- 局部影响：计算某个网格点P的值时，并不是所有测点都同等重要。通常距离P点越近的测点，影响力越大。
- 平滑过渡：样条函数（特别是薄板样条）就像一个最小化整体弯曲能量的弹性薄板。它要满足两个条件：
精确（或近似）通过已知点：膜要被测点顶到正确高度。
整体尽可能平滑：在没有测点的地方，膜自然过渡，不会出现不该有的陡坡或凹陷。
3. 两种主要的样条栅格插值
- 规则化样条（平滑样条）
像一张有点厚度的橡胶膜。
它不要求曲面严格穿过每一个测点，而是允许有一点微小的误差。
核心思想是在"贴合测点"和"保持曲面平滑"之间找一个平衡。
适用：当数据可能有微小误差（比如仪器误差）时，避免曲面为了强行穿过有误差的点而产生不真实的波动。 - 薄板样条（TPS）
像一张极薄的、完美的弹性钢片膜。
它要求曲面必须精确穿过每一个已知测点（完全精确插值）。
在没有测点的地方，它会让曲面以最平滑的方式（数学上叫"最小曲率"）延展过去。
你可以把它想象成：在测点位置把膜焊死在指定高度，然后让膜自己找到能量最低的平滑形状。
适用：数据非常精确，且我们相信测点值绝对可靠的情况（如高程点）。
二、为什么用样条做栅格插值？
- 优点：
生成表面非常美观光滑：适合需要视觉平滑输出的地图（如地质层面、污染物扩散模拟）。
数学基础严谨：有明确的"最小化曲率"物理意义，结果可解释性强。
适用于渐变性数据：比如温度、气压、高程等自然连续变化的 phenomena。 - 缺点/注意事项：
可能产生"过冲"：如果两个很近的测点值差异巨大，样条为了光滑地连接它们，可能会在中间"拱"出一个比两者都高，或"凹"出一个比两者都低的不真实值。就像拉紧一根穿过一高一低两颗钉子的橡皮筋，中间部分会脱离木板。 - 计算量相对较大：因为要解一个方程组来确定整个曲面的形状。
- 不适合突变数据：比如断层线两侧的地层高程突然变化，用追求全局光滑的样条就不合适了（会模糊掉断层）。
总结
样条函数栅格插值，就是利用样条"分段光滑连接"的思想，把离散的、不规则的空间点数据，"编织"或"蒙皮"成一张覆盖整个区域的、连续光滑的栅格曲面图。它像一位追求优雅和平滑的艺术家，在点与点之间描绘出最自然的过渡。