基于BP神经网络的高分辨率海底地形跨层生成模型.pdf

DOI:10.11991/yykj.202305007网络出版地址：https:/ BP 神经网络的高分辨率海底地形跨层生成模型王振1,2，张锡亭3，王建华41.北京航空航天大学数学科学学院，北京1022062.大连理工大学工业装备与结构分析国家重点实验室，辽宁大连1160243.大连理工大学数学科学学院，辽宁大连1160244.长沙矿冶研究院有限责任公司，湖南长沙410006摘要：为了满足海底地形的高分辨率需求及解决测量数据的有限性问题，基于多层前馈神经网络（backpropagation,BP）和跨层网格生成策略，建立了兼顾海底区域地形整体特征和局部地形信息的海底地形跨层生成模型，实现对海底地形数据生成填充。以南海海底地形为例，通过误差对比、假设检验以及海底地形云图的图像清晰度对本文模型生成数据进行有效性验证。结果显示所建立的模型在保证与原始数据之间误差小和数据特征相同的前提下完成了对地形云图的图像清晰度的提升，并且结果优于传统克里金插值方法。本文分析结果可为地形数据相关研究提供参考。关键词：高分辨率海底地形；跨层网格；BP 神经网络；克里金插值；Mann-WhitneyU检验；Levene检验；图像清晰度；误差中图分类号：P229.1文献标志码：A文章编号：1009671X(2024)01014309High-resolution submarine topography cross-layer generation model viaBP neural networkWANGZhen1,2,ZHANGXiting3,WANGJianhua41.SchoolofMathematicalScience,BeihangUniversity,Beijing102206,China2.StateKeyLaboratoryofStructureAnalysisforIndustrial,DalianUniversityofTechnology,Dalian116024,China3.SchoolofMathematicalScience,DalianUniversityofTechnology,Dalian116024,China4.ChangshaResearchInstituteofMiningandMetallurgyCo.,Ltd.,Changsha410006,ChinaAbstract:Tomeetthehigh-resolutiondemandofthesubmarinetopographyunderthelimitedamountofmeasurementdataavailable,across-layergenerationmodelisestablishedviamulti-layerbackpropagation(BP)neuralnetworkandcross-layer grid generation strategy,which incorporates both the overall characteristics of the area of submarinetopographyandlocalterraininformation.TakingthetopographyoftheSouthChinaSeaasanexample,thispaperillustratestheeffectivenessbyemployingerrorcomparison,hypothesistestingandtheclarityofsubmarinetopographycloudimages.Theexperimentresultsshowthattheproposedmodelinthispapercanimprovetheclarityoftopographiccloudimageunderthepremisethattheerrorbetweentheoriginaldataandthemodelsgeneratedsubmarinetopographydataissmallandthedatafeaturesbetweenthemarethesame.Inaddition,theresultoftheproposedmodelisbetterthanthetraditionalKriginginterpolationmethod.Theresultsofthispapercanprovidereferencefortherelatedresearchoftopographicdata.Keywords:high-resolutionsubmarinetopographic;cross-layergrid;BPneuralnetwork;Kriginginterpolation;Mann-WhitneyUtest;Levenetest;clarityofimage;error地球上海洋总面积约为 3.6 亿平方千米，占地球表面积的 71%，蕴含巨大资源1。海底地形作为海洋研究的基础地理信息，对海洋资源开采、水下设备研发以及对海底构造运动和海洋环流的研究具有重要意义2。但由于海底环境的特殊性以及采样困难问题，造成海底地形数据网格精度收稿日期：20230511.网络出版日期：20231207.基金项目：国家重点研发计划项目（2022YFC2806802）；辽宁省“兴辽英才计划”项目(XLYC1907014，XLYC1908027)；中央高校基本科研业务费项目(DUT21ZD205，DUT20TD108)；大连市杰青项目(2019RJ05).作者简介：王振，男，教授，博士.张锡亭，女，硕士.通信作者：王振，E-mail：.第51卷第1期应用科技Vol.51No.12024年1月AppliedScienceandTechnologyJan.2024不足，这给海洋事业发展带来挑战。海底地形复杂多变，高分辨率的水下地形的生成对于船舶、水下航行器等海洋工程装备的路径规划和安全行驶十分重要。提高海底地形图像的分辨率面临诸多困难，如何通过对有限测量数据的有效填充生成高分辨率地形是急需解决的问题。国内外研究人员提出很多种数据插值算法。如最近邻域插值和线性插值，但由于最近邻域插值在运算时直接选用邻近的数据进行插值计算，得到的图像的边缘容易是锯齿状的，视觉效果较差。而线性插值数据量小、速度快，相对于最近邻域插值方法来讲，插值效果较好，但是其插值过程中会丧失一些高频成分，使得图像处理后变得有些模糊。为了改善上述情况，基于高阶多项式的多项式内插算法被应用于图像处理中，如 Keys3提出立方卷积插值和Hou 等4提出三次样条插值改变了图像边缘锯齿状以及图像模糊的问题，但是计算复杂、计算量大。于是，Reed 等5提出了基于最小二乘法和线性插值方法的线性样条插值方法以及 Truong 等6提出基于最小二乘法与立方卷积插值函数的立方卷积样条插值算法，经验证比起先前的插值算法，会获得更清晰的图像，但是其计算复杂度相对较大,特别是算法中的计算循环卷积部分。除上述插值方法外，比如 Schafer 等7提出基于贝叶斯方法和最大似然的插值算法，Li 等8提出基于模糊聚类的插值算法，这些数据插值算法在填充数据量小的问题上取得显著的效果，然而面对大规模数据填充任务时，精度急剧下降，这是由于大规模数据存在着丰富的数据信息维度，而这些数据插值算法不能有效体现大规模数据的深度特征9。在地形数据处理方面，南非矿产地理学家Krige 提出了克里金插值算法。克里金插值算法相比于传统线性插值算法具有更高的精度，在地质统计学及图像处理等方面获得了广泛应用10，但它本质上是一种光滑的内插值算法，生成图像较光滑，然而真实的地表具有复杂性，用这种插值算法建立的模型未能体现真实地表最重要的特征11。为了更好地提高图像的质量，近几年许多新方法被提出，比如 Ramponi 等12提出空间可变的插值算法,该方法利用相邻点之间的不对称性引入加偏差的距离，从而有效保护了图像的局部特征，进而提高重构图像的品质。Thurnhofer 等13提出自适应插值算法，根据原始图像的不同区域的不同的结构特性，采用的不同的插值方法。为了更好完成大规模数据填充任务，大热的机器学习方法被用于数据插值中，如 Narang 等14提出 K近邻插值算法；Wang 等15提出基于支持向量回归方法的数据插值方法；Shi 等16提出 PixelShuffle 算法；Tian 等17提出 DUpsampling 算法；张守建等18提出基于多层前馈神经网络（backpropagation,BP）的数据填充方法，应用于我国西南地区地形数据中。伴随着海底开发工程装备的发展，对海底高分辨率地形的需求变得十分迫切，但现有方法大多是基于邻点的经典数据填充格式，很难反应海底地形的地貌特征，而利用有限采样数据获得相邻样本间的内在联系，同时保持整体特征的地形生成方法值得深入研究探索。本文在有限采样数据下，通过 BP 神经网络建立具有跨层生成网格样本能力的海底地形生成模型，对海底地形样本填充，形成高分辨率海底地形图像。1高分辨率海底地形跨层生成模型在保持样本数据整体特征的前提下生成满足高分辨率要求的海底地形，本文建立具有跨层网格生成效果的 BP 神经网络模型对原始网格型海底样本加密填充，从而获得高分辨率海底地形图像。1.1模型建立k(mn)k(33)kk1本文针对如图 1 所示的规则网格地形样本进行加密填充。图 1 中黑色圆形实心点代表测量地形已有样本网格点，记为第 1 层网格；绿色菱形实心点代表目标生成的网格点，记为第 0 层网格。将第 1 层网格每隔一行采样一行，每隔一列采样一列，得到的更稀疏的网格样本记为第 2 层网格，同理可生成更稀疏的高层网格。在第 层网格中取个网格点组成第 层网格的网格单元，第一层网格的网格单元如图 1 中红框所示。在每个第 层网格单元中，相邻的上下左右 4 个等间隔网格点构成的小单元的中间位置为其对应的第层网格。第 1 层地形网格点第 0 层地形网格点第 1 层网格网格单元图1海底地形样本网格点示意由于空间自相关性，同一区域地形数据具有相同的统计特征和随机性18，因此通过大量网格144应用科技第51卷kk1单元样本训练，神经网络模型已将区域海底地形地貌数据随机特征和统计特征融入跨层网格样本间非线性映射关系中。采用这种方法得到的跨层网格生成模型既能表达网格点之间的相邻关系，也能表达区域地形的随机统计特征。利用 3 层BP 神经网络1925构建高分辨率海底地形跨层生成模型，将第 层网格地形样本和第层网格地形样本作为训练集，通过 BP 神经网络获得两者间的非线性映射关系，形成具有跨层网格生成能力的 BP 神经网络模型。Z(1)=W(1)Xk+b(1)(1)A(2)=f(Z(1)(2)Ok1=W(2)A(2)+b(2)(3)W(1)W(2)b(1)b(2)Z(1)A(2)f式中：、分别为输入层到隐藏层以及隐藏层到输出层权重矩阵；、分别为输入层到隐藏层和隐藏层到输出层的阈值矩阵；为输入层的输出值构成的矩阵；为隐藏层的输出值构成的矩阵；为激活函数，用于隐藏层的传递，以达到增强 BP 神经网络模型非线性逼近能力的目的。kXkk1Xk1Xkk1Ok1假设模型训练样本为第 层原始海底地形网格样本及其对应的第层原始海底地形网格样本，样本以输入层隐藏层输出层的方向26，按照式(1)式(3)运算，得到经 BP 神经网络正向传播运算后的第层海底地形网格样本。Ok1Xk1考虑到正向传播结束后样本与样本之间有误差，根据两者的均方误差利用梯度下降法，沿输出层隐藏层输入层的方向，进行误差反向传播27，从而得到更新后的层间连接权重矩阵和阈值矩阵。Xk1=Fkk1(Xk)(4)XkXk1经过正向传播和反向传播之间的反复迭代，直到误差达到容许的范围或者循环达到设定的次数，则训练过程完毕。BP 神经网络会保存一套层间连接权重矩阵和阈值矩阵，形成如式(4)所示的用于反映训练集中样本与样本之间非线性映射关系的跨层网格生成模型。Fkk1pXpFkk1Xp为检验跨层网格生成模型对第 层海底地形网格样本的适用性，将应用于。Xp1=Fkk1(Xp)(5)p1Xp1XpFkk1,pp1按照式(5)，生成第层海底地形网格样本，完成对样本的加密填充。以上过程记为。Fkk1p1通过已有样本数据检验模型对第层数据的生成质量，确定上述训练模型的适用性，F21X0=F21(X1)也部分说明跨层网格样本之间具有相似的数据结构。因此可以利用训练好的模型应用于第1 层网格样本，也是给定的最精细网格，生成更精细的第 0 层网格样本，即，得到最终加密填充后的样本数据。1.2模型检验为验证模型是否可以在不改变样本统计特征的前提下，完成对海底地形样本的加密填充并提高 海 底 地 形 图 像 的 清 晰 度，本 文 使 用 Mann-WhitneyU 检验、Levene 检验和图像清晰度 3 个指标进行比较、验证。1.2.1Mann-WhitneyU 检验Xk1Ok1WXk1Ok1Xk1Ok1mnXk1Ok1ZMann-WhitneyU 检验2829假设原始海底地形样本和模型生成的海底地形样本分别来自除了均值以外完全相同的 2 个总体，由此检验2 组样本的均值是否有显著的差别。令表示样本观察值小于样本观察值的个数，和 分别表示和的样本个数，则检验统计量 为Z=WXk1Ok1mn2mn(mn+1)12Z检验统计量 服从标准正态分布，当求得的P 值大于显著性水平 0.05 时，认为 2 组地形样本均值相等。1.2.2Levene 检验Xk1Ok1Z1.Z2.Xk1Ok1Z1.Z2.Z.Z1.Z2.WLevene检验3031是一种方差齐性检验的方法。利用该检验方法来检验原始海底地形样本和模型生成的海底地形样本的方差是否相同。令和分别为和经转化后的数值，和分别为转化后 2 组地形样本的均值，表示和全部为地形样本的均值。转化公式及检验统计量为Z1=?Xk1Xk1?Z2=?Ok1Ok1?W=(m+n2)(m(Z1Z)2+n(Z2Z)2)mj=1(Z1jZ1)2+nj=1(Z2jZ2)2WF(1,m+n2)检验统计量服从分布，当求得的 P 值大于显著性水平 0.05 时，认为 2 组地形样本之间方差相等。1.2.3图像清晰度MN本文采用基于频率域特征的评价函数32来计算像素矩阵为型矩阵的海底地形云图的图第1期王振，等：基于 BP 神经网络的高分辨率海底地形跨层生成模型145像清晰度，计算过程为：1）将海底地形云图 RGB 图像转化为对应的灰度图像，利用加权平均法计算图像每一点的灰度值。f(x,y)=0.229R+0.587G+0.114B(x,y)R、G、BR、G、B0 255式 中：为 图 像 矩 阵 中 每 一 点 的 坐 标；为三元光色红、绿、蓝强度值，任何一种颜色都可以由这 3 种颜色表示，取值范围都为。2）对转变为灰度值后的海底地形云图图像数值矩阵进行二维离散傅里叶变换，转换到频域。F(,)=M1x=0N1y=0f(x,y)ei2(xM+yN)3）计算海底地形云图图像清晰度数值。P(,)=R(,)2+I(,)2G=1MN=M1=0N1=02+2P(,)R(,)I(,)F(,)G式中：和分别为傅里叶变换的实部和虚部；为图像的清晰度数值，该数值越高，图像越清晰。2南海海底地形生成与分析(33)(66)(1212)利用 Python 搭建海底地形跨层生成模型，完成对原始网格型海底地形样本加密填充。为了确保模型有效，选取南海海底网格型样本作为算例进行训练分析。采用 ETOPO5(https:/www.ngdc.noaa.gov/mgg/global/relief/ETOPO5/TOPO/ETOPO5/)中南海海域网格型海底地形数据作为样本进行训练，具体为东经 105.0042124.9958,北纬21.362524.9958。选取 3 种跨层取样方式：将整体网格型原始海底网格样本按照紧密相连的个、个、个网格形成的网格单元进行分割。在 3 种网格单元中分别以 1 个、2 个、4 个网格上下左右等间隔取样，每个网格单元都分别含有 16 个第 1、2、3 层样本取样位置，9 个相对应的第 0、1、2 层样本取样位置。取样位置举例说明如图 2 所示。第 1 层样本取值位置与其对应的第 0 层样本取值位置第 2 层样本取值位置与其对应的第 1 层样本取值位置图2取样位置图像F21F21,21F32F32,32F32F32,21将跨层取样后的海底地形样本划分 3 种不同的训练集和测试集的组合验证模型效果。它们分别为：1）将 171570 组第 2 层样本和其对应的第 1层样本随机分为 170570 组训练集和 1000 组测试集。利用模型获得训练集中第 2 层样本与第 1 层样本之间的映射关系，将其作用于测试集的第 2 层样本中，生成相应的新第 1 层样本，与这种样本组合相关的模型运行过程用标记。2）将 83580 组第 3 层样本和其对应的第 2 层样本随机分为 82580 组训练集和 1000 组测试集。通过模型将训练集中的第 3 层样本与第 2 层样本之间的映射关系作用于测试集的第 3 层样本中，生成相应的新第 2 层样本，与这种样本组合相关的模型运行过程用标记。3）将 83580 组第 3 层样本和其对应的第 2 层样本作为训练集，随机选取 1000 组第 2 层样本及其对应的第 1 层样本作为测试集。通过模型将训练集中的第 3 层样本与第 2 层样本之间的映射关系作用于测试集的第 2 层样本中，生成相应的新第 1 层样本，与这种样本组合相关的模型运行过程用标记。在算例运行过程中，为使各变量处于同等重要位置，保证模型运行结果的可靠性，对样本标准化后再进入模型训练。2.1海底地形跨层生成F21,21F32,32F32,21本文 BP 神经网络的输入层、隐藏层和输出层的神经元个数，分别为 16、40 和 9，学习率为0.0001，一次训练选取的训练样本组数为 7，所有训练样本最大训练轮次为 400，选取的激活函数为 tanh 激活函数，根据激活函数，本文使用的权重矩阵初始化方法为 Xavier 初始化3334。、共 3 种组合模型运行结果如图 38 所示。1012345高度/103 m22.3722.3622.3522.3422.3322.32纬度(N)/()经度(E)/()125.0122.5120.0117.5115.0112.5110.0107.5105.05 0004 0003 0002 0001 0001 0000F21,21图3第1层原始海底地形样本构成的三维图像146应用科技第51卷1012345高度/103 m22.3722.3622.3522.3422.3322.32纬度(N)/()经度(E)/()125.0122.5120.0117.5115.0112.5110.0107.5105.05 0004 0003 0002 0001 0001 0000F21,21图4模型产生的第1层海底地形样本构成的三维图像5 0004 0003 0002 0001 0001 00001012345高度/103 m22.3622.3822.4022.4222.3422.32纬度(N)/()经度(E)/()125.0122.5120.0117.5115.0112.5110.0107.5105.0F32,32图5第2层原始海底地形样本构成的三维图像1012345高度/103 m22.3622.3822.4022.4222.3422.32纬度(N)/()经度(E)/()125.0122.5120.0117.5115.0112.5110.0107.5105.05 0004 0003 0002 0001 0001 0000F32,32图6模型产生的第2层海底地形样本构成的三维图像1012345高度/103 m22.3422.3522.3622.3722.3322.32纬度(N)/()经度(E)/()125.0122.5120.0117.5115.0112.5110.0107.5105.05 0004 0003 0002 0001 0001 0000F32,21图7第1层原始海底地形样本构成的三维图像1012345高度/103 m22.3422.3522.3622.3722.3322.32纬度(N)/()经度(E)/()125.0122.5120.0117.5115.0112.5110.0107.5105.05 0004 0003 0002 0001 0001 0000F32,21图8模型产生的第1层海底地形样本构成的三维图像图 3、图 5 和图 7 表示 3 种不同组合方式下测试集中原始海底地形样本构成的三维图像；图 4、图 6 和图 8 表示 3 种不同组合方式下由模型生成的海底地形样本构成的三维图像。由同种组合方式的2 个图像对比可知，模型生成的海底地形样本与原始海底地形样本相差不大。为更客观地评价海底地形跨层生成模型结果，本文不仅对三维海底地形图像进行对比，并且进行理论分析。如表 1 所示，将 3 种不同组合的模型运行结果与克里金插值算法对比，以均方误差、相对误差、决定系数、Mann-WhitneyU检验和 Levene 检验为对比指标进行比较。表1不同填充方式的误差及假设检验结果对比填充方式均方误差相对误差/%决定系数Mann-WhitneyU检验Levene检验F21,210.00161.49520.99840.53540.5368F32,320.00221.78250.99780.25680.8283F32,210.00081.16870.99910.62970.9485克里金插值0.08236.89360.91771.113 910300.2200第1期王振，等：基于 BP 神经网络的高分辨率海底地形跨层生成模型147F21,21F32,32F32,21根据表 1 可知，由、和这3 种组合方式对应的高分辨率海底地形跨层生成模型产生的海底地形样本与原始海底地形样本之间均值和方差相同的假设都成立，而由第 2 层原始海底地形样本和克里金插值方法共同产生的第 1 层海底地形样本与第 1 层原始海底地形样本之间不能通过均值相等的检验。对于模型生成的海底地形样本与原始海底地形样本之间的均方误差和相对误差，模型低于克里金插值。在决定系数方面，模型高于克里金插值。因此，该模型可以在不改变数据统计特征的前提下，利用跨层网格生成策略完成对海底地形数据的填充。2.2有效性验证F32,21F21,10由上文可知，本文建立的海底地形跨层生成模型可以有效获取跨层样本间的非线性映射关系，并运用到其他层网格样本生成中。为完成原始 网 格 型 海 底 样 本 的 加 密 填 充，本 文 参 照运行模式，利用运行模式，获得第 0层海底地形样本，完成对原始海底地形样本的加密填充。为体现高分辨率海底地形生成模型对样本的填充效果，随机选取该海底 2 个局部地区，得到 2 个地区的样本填充前后的三维海底地形对比图像，如图 918 所示。同时，利用 Mann-WhitneyU检验和 Levene 检验来检验 2 个地区样本填充前后的统计特征的一致性，并由克里金插值对相同区域的第 1 层原始海底地形网格样本加密填充，以图像清晰度为依据比较两者生成的海底地形云图图像的质量，具体结果如表 2 和表 3所示。1 7501 5001 2501 0007505002500高度/m25.024.824.624.424.224.0纬度(N)/()经度(E)/()114.0113.5113.0112.5112.0111.5111.0110.5110.02000400 600 800 1 0001 2001 4001 600图9局部地区 1 原始样本海底地形三维图像1 7501 5001 2501 0007505002500高度/m25.024.824.624.424.224.0纬度(N)/()经度(E)/()114.0113.5113.0112.5112.0111.5111.0110.5110.02000400 600 800 1 0001 2001 4001 600图10局部地区 1 模型填充后海底地形三维图像114.0113.0112.0111.0110.024.024.224.4纬度(N)/()24.624.825.02000400600800 1 000 1 200 1 400 1 600高度/m经度(E)/()图11局部地区 1 模型原始样本海底地形云图114.0113.0112.0111.0110.024.025.024.224.4纬度(N)/()24.624.82000400600800 1 000 1 200 1 400 1 600高度/m经度(E)/()图12局部地区 1 模型填充样本后海底地形云图经度(E)/()24.024.224.4纬度(N)/()24.624.825.02000400600800 1 000 1 200 1 400 1 600高度/m114.0113.0112.0111.0110.0图13局部地区 1 克里金插值后海底地形云图148应用科技第51卷1 7501 5001 2501 000750500250高度/m25.024.824.624.424.224.0纬度(N)/()经度(E)/()109.0108.5108.0107.5107.0106.5106.0105.5105.0200 400 600 800 1 0001 2001 4001 6001 800图14局部地区 2 原始样本海底地形三维图像2 0001 7501 5001 2501 000750500250025.024.824.6高度/m24.424.224.0经度(E)/()纬度(N)/()109.0108.5108.0107.5107.0106.5106.0105.5105.0200 400 600 800 1 0001 2001 4001 6001 800图15局部地区 2 模型填充后海底地形三维图像109.0108.0107.0106.0105.0经度(E)/()200400600800 1 000 1 200 1 400 1 600 1 80024.024.224.4纬度(N)/()高度/m24.624.825.0图16局部地区 2 原始样本海底地形云图109.0108.0107.0106.0105.024.025.024.224.4纬度(N)/()24.624.8经度(E)/()200400600800 1 000 1 200 1 400 1 600 1 800高度/m图17局部地区 2 模型填充样本后海底地形云图109.0108.0107.0106.0105.0经度(E)/()24.025.024.224.4纬度(N)/()24.624.8200400600800 1 000 1 200 1 400 1 600 1 800高度/m图18局部地区 2 克里金插值后海底地形云图表2本文模型生成地形和克里金插值结果的数据的有效性对比(局部地区 1)填充方式Mann-WhitneyU检验Levene检验填充前图像清晰度填充后图像清晰度本文模型0.15860.27545.825 81066.511 9106克里金插值1.753 510340.00005.825 81063.814 9106表3本文模型生成地形与克里金插值的数据的有效性结果对比(局部地区 2)填充方式Mann-WhitneyU检验Levene检验填充前图像清晰度填充后图像清晰度本文模型0.67020.23715.304 01065.949 9106克里金插值1.539 110302.335 2101355.304 01063.831 1106由图 910 对比以及图 1415 对比可知，本文高分辨率海底地形跨层生成模型可以对地形数据进 行 合 理 的 填 充，并 且 由 图 1113 对 比 以 及图 1618 对比可知，克里金插值方法获得的图像过于光滑，没有体现出海底的地形特征，而本文模型生成的图像在提升清晰度的同时，保证了生成的海底地形图像有效性。由表 2 和表 3 可知，本文构建的高分辨率海底地形跨层生成模型对样本填充效果良好。Mann-WhitneyU检验和Levene检验结果显示，填充前后样本的均值和方差没有显著差异。通过图像清晰度计算结果可知，该模型对样本填充后图像的清晰度高于原始样本图像以及克里金插值填充后样本图像的清晰度。3结论本文基于 BP 神经网络建立了具有跨层网格生成能力的高分辨率海底地形跨层生成模型，通过南海海底地形验证了该模型产生的数据兼顾了原始样本的相邻数据关系和整体数据特征。第1期王振，等：基于 BP 神经网络的高分辨率海底地形跨层生成模型1491）模型利用 BP 神经网络的泛化能力，将高层网格样本与其对应的低层网格样本之间的非线性映射关系迁移至其他层网格样本生成中，并根据均方误差与相对误差验证了该模型具有有效的网格样本跨层生成能力。2）通过假设检验，验证了该模型生成样本可以表征出原始样本的统计特征。3）由云图图像及其清晰度的计算结果表明，相较于传统克里金插值方法，该模型的结果良好。因此，本文模型可以有效获取并迁移跨层网格间的相关关系至其他层网格数据生成中，完成对海底样本的加密填充，生成高分辨率图像用于展现海底地形。机器学习和跨层网格生成方法的结合方式给海底地形的研究提供了一种新的思路，高分辨率海底地形跨层生成模型在更大地形范围以及更多海山、海脊等海底地貌特征方面具有广阔的应用前景。参考文献：郝瑞杰,万晓云,眭晓虹,等.海底地形探测和模型研制现状及精度分析 J.地球与行星物理论评,2022,53(2):172186.1王轲.海底地形变化影响因子定量分析研究以胶州湾为例 D.山东:自然资源部第一海洋研究所,2021.2KEYS R.Cubic convolution interpolation for digitalimage processingJ.IEEE transactions on acoustics,speech,andsignalprocessing,1981,29(6):11531160.3HOU H,ANDREWS H.Cubic splines for imageinterpolationanddigitalfilteringJ.IEEEtransactionsonacoustics,speech,and signal processing,1978,26(6):508517.4REEDIS,ALFREDCY.Optimalsplineinterpolationforimagecompression:U.S.,582256P.19981013.5TRUONG T K,WANG L J,REED I S,HSIEHWS.Image data compression using cubic convolution splineinterpolationJ.IEEEtransactionsonimageprocessing:apublicationoftheIEEEsignalprocessingsociety,2000,9(11):8895.6SCHAFERJL,GRAHAMJW.Missingdata:ourviewofthestateoftheartJ.Psychologicalmethods,2002,7(2):147177.7LIDan,JITENDERD,BILLS,etal.Towardsmissingdata imputation:a study of fuzzy k-means clusteringmethodC/Rough Sets and Current Trends inComputing.Heidelberg:Springer,2004:573579.8卜范玉,陈志奎,张清辰.基于深度学习的不完整大数据填 充 算 法 J.微 电 子 学 与 计 算 机,2014,31(12):173176.9申静,苏天赟,王国宇,等.基于 Kriging 算法的海底地形插值设计与实现 J.海洋科学,2012,36(5):2428.10于飞,陈斐楠,马慧,等.复合分形插值算法在海底地形仿真中的应用 J.计算机仿真,2015,32(1):254258.11RAMPONI G.Warped distance for space-variant linearimage interpolationJ.IEEE transactions on imageprocessing:a publication of the IEEE signal processingsociety,1999,8(5):629639.12THURNHOFER S,LIGHTSTONE M L,MITRA S K.Adaptive interpolation of images with application tointerlaced-to-progressiveconversionC/VisualCom-munications.Cambridge:SPIE,1993:614625.13NARANGSK,GADDEA,SANOUE,etal.Localizediterative methods for interpolation in graph structureddataC/SignalandInformationProcessing.Austin:IEEE,2013:491494.14WANGXian,LIAo,JIANGZhaohui,etal.MissingvalueestimationforDNAmicroarraygeneexpressiondatabysupport vector regression imputation and orthogonalcodingschemeJ.BMCbioinformatics,2006,7(1):32.15SHIWenzhe,CABALLEROJ,HUSZRF,etal.Real-time single image and video super-resolution using anefficient sub-pixel convolutional neural networkC/Computer Vision and Pattern Recognition.Las Vegas:IEEE,2016:18741883.16TIAN Zhi,HE Tong,SHEN Chunhua,et al.Decodersmatterforsemanticsegmentation:data-dependentdecodingenablesflexiblefeatureaggregationC/ComputerVisionandPatternRecognition.LongBeach:IEEE,2019:31213130.17张守建,李建成,王正涛,等.基于 BP 神经网络的SRTM 空缺区域数据填充方法研究 J.大地测量与地球动力学,2008,28(5):4.18YUAN Mei,LI Chunyang.Research on global highereducationqualitybasedonbpneuralnetworkandanalytichierarchyprocessJ.Journalofcomputerandcommunication,2021,9(6):158173.19LIU Yang.Solar power prediction algorithm based onimprovedBPneuralnetworkJ.Computerinformatizationandmechanicalsystem,2023,6(6)