17、MATLAB空间数据处理与分析

MATLAB空间数据处理与分析

1. 网格化示例

1.1 双调和样条插值方法

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

1.2 示例数据

我们使用合成的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值。使用以下代码