17、MATLAB空间数据处理：网格化、插值方法及点分布统计分析

MATLAB空间数据处理：网格化、插值方法及点分布统计分析

1. 网格化示例

1.1 数据准备

MATLAB 从一开始就提供了由 Sandwell（1987）开发的双调和样条插值方法，该方法特别适合从具有不均匀分布控制点的噪声数据集中生成平滑表面。

以一个合成的 xyz 数据为例，该数据表示一个被正断层位移的假想地层界面与参考面之间的垂直距离。断层的下盘显示大致水平的地层，而上盘的特征是两个大型沉积盆地的发育。xyz 数据不规则分布，因此需要插值到规则网格上。数据存储在名为 normalfault.txt 的文件中，文件内容如下：

4.3229698e+02   7.4641694e+01   9.7283620e-01
4.4610209e+02   7.2198697e+01   6.0655065e-01
4.5190255e+02   7.8713355e+01   1.4741054e+00
4.6617169e+02   8.7182410e+01   2.2842172e+00
4.6524362e+02   9.7361564e+01   1.1295175e-01
4.5526682e+02   1.1454397e+02   1.9007110e+00
4.2930233e+02   7.3175896e+01   3.3647807e+00

第一列和第二列包含坐标 x（在任意空间坐标系中为 420 到 470 之间）和 y（在 70 到 120 之间），而第三列包含垂直 z 值。使用以下代码加载数据：