海底地形测量图的插值模型_海底地形

海底地形测量图的插值模型

摘要

随着全球经济一体化和信息技术的发展，企业之间的合作日益加强，跨地区甚至跨国合作制造的趋势日益明显。国际上越来越多的制造企业不断地将大量常规业务“外包”（outsourcing）出去给发展中国家，而只保留最为核心的业务（如市场、关键系统设计和系统集成、总装配以及销售）。在这些合作生产的过程中，大量的物资和信息在更为广阔的地域间转移、储存和交换，国际物流活动将日益频繁，港口作为国际物流活动主要的载体，在国际贸易与国际经济合作中愈来愈发挥着极其重要的作用。但是海底的地形是十分复杂的，它不仅分布有巍峨的海底山脉、平缓的海底平原，而且还有许多陡峭的海底深沟。为使轮船进出港口安全，就需要了解港口航道的海底地形。所以我们就要对海底地形有足够的了解，预测出哪些区域是船只的禁入点，避免船只在危险区域搁浅。为了研究此问题，我们利用低潮时测得14组数据，并进行了6个基本假设，在此基础上我们便利用插入与拟合的思想来利用光滑曲面来模拟海底曲面。但考虑到保凸性及光滑性的要求，我们采用双三次样条函数来模拟海底曲面。首先利用测量的14个随机点的深度，以随机点的坐标将待测矩形区域划分成14?14的网格。然后我们通过某种加权平均来逼近未知网格接点上的深度，采用距离的倒数作权重反映出距离越小影响越大。求出所有14×14个网格点上的深度后，调用IMST中的双三次样条子程序，通过插值得到海底曲面，然后再加细网格，划分成50×50的网格，计算这加细网格接点上的深度;最后找出两个危险区域分别在深度为4Ft的

两个点的周围，并借助于Matlab中的绘图程序，绘出海底的二维、三维网格图及等高线图。以不同的颜色将水深小于船只的吃水深度标示出来，作为船只的避免进入区域，并作出水深小于船只吃水深度的海域范围,并绘出等高线。在找几个近似大小的海底拓扑地图，每个随即取14个数据点，把我们的模型应用上去，通过比较实际地图及由模型画出地图比较接近。这说明该模型具有十分有效的实际作用。该模型充分利用了已知点的数据，给出了求未知网格点上深度的近似方法，用保凸性较好的双三次样条拟合了海底曲面，得到了较为满意的结果。但在实际计算中，三次样条可能会导致数值上的不稳定，遇到这种情况，可以用加细网格点的办法来解决，也可以用稳定性较好的B—样条来拟合。

一、问题的提出

随着全球经济一体化和信息技术的发展，企业之间的合作日益加强，跨地区甚至跨国合作制造的趋势日益明显。国际上越来越多的制造企业不断地将大量常规业务“外包”（outsourcing）出去给发展中国家，而只保留最为核心的业务（如市场、关键系统设计和系统集成、总装配以及销售）。在这些合作生产的过程中，大量的物资和信息在更为广阔的地域间转移、储存和交换，国际物流活动将日益频繁，港口作为国际物流活动主要的载体，在国际贸易与国际经济合作中愈来愈发挥着极其重要的作用。

但是海底的地形是十分复杂的，它不仅分布有巍峨的海底山脉、平缓的海底平原，而且还有许多陡峭的海底深沟。为使轮船进出港口安全，就需要了解港口

航道的海底地形。

表1给出了水面直坐标（x，y）处的水深z，这是在低潮测得的。如果船的吃水深度为5Ft，试问在矩形区域(75，200)×(50，50)中行船应避免进入哪些区域？其中水平方向的坐标x, y以Yd(=0.914m)为单位，水深方向以Ft(=30.48cm)为单位。

二、模型假设

海底是光滑的，且无暗礁;

每个给定的数据都影响着其它未知点的深度，且距离越近，影响越大; 任何两个数据点之间深度的变化都影响着其它未知点的深度;

对于给定两个数据的变化对于某一未知点深度的影响，取决于三个距离两个数据点连线与该未知点的垂直距离、该未知点与离它最近的那个数据点之间的距离及两个给定数据点之间的距离;

两个数据点之间深度的变化对某一未知点深度的影响与这两个数据点连线线性传播; 每一个给定的数据点对某一未知点深度的影响它们之间的距离平方成反比。三、分析与建摸

根据假设条件海底是光滑的，无暗礁，因此很自然地想到利用光滑曲面来拟合海底曲面。例如可以用二维拉格朗日(Lagrange)插值或双三次样条函数来逼近。考虑到保凸性及光滑性的要求，我们采用双三次样条函数来拟合.

为了用双三次样条函数插值，必须知道xoy平面内所

有网格上的深度，而所给数据的14个随机点并不构成任何网格。所以第一步生成网格

我们采取最简单的方法是过个14 数据点分别作平行于x，y 轴的直线，划分成不规则14?14的网格。

第二步确定那些未知点数据网格上的深度应该说所有数据对未知网格的深度都有影响，只是越靠近的数据点影响越大。由于我们对海底所知甚少，所以只能通过某种加权平均来逼近未知网格接点上的深度，采用距离的倒数作权重反映出距离越小影响越大。但仅用加权平均来逼近未知点的深度，它不能反映数据点深度的变化趋势。

设Q1，Q2点的深度分别为4Ft和8Ft，G是Q1Q2连线上未知深度的点(见图1)，

Q1Q2?Q2G?4Ft，求G点的深度zg。

根据光滑假设，由点Q2经点Q1到G点的深度应渐渐变浅，因此，未知点G的深度应小于4Ft。

4

Q2

4

G

下面用三种外推公式加以分析

Z(Q1,Q1,G)Zg?

GQ1

2