留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

姓名
邮箱
手机号码
标题
留言内容
验证码

基于全球测深数据的中国海岸线周边海域数字水深模型融合

阮晓光 占赵杰 闫兆进 谈秋英 郭美静 杨阳

阮晓光,占赵杰,闫兆进,等. 基于全球测深数据的中国海岸线周边海域数字水深模型融合[J]. 海洋学报,2024,46(7):16–28 doi: 10.12284/hyxb2024062
引用本文: 阮晓光,占赵杰,闫兆进,等. 基于全球测深数据的中国海岸线周边海域数字水深模型融合[J]. 海洋学报,2024,46(7):16–28 doi: 10.12284/hyxb2024062
Ruan Xiaoguang,Zhan Zhaojie,Yan Zhaojin, et al. Digital bathymetric model fusion of offshore waters around China’s coastline based on global bathymetry data[J]. Haiyang Xuebao,2024, 46(7):16–28 doi: 10.12284/hyxb2024062
Citation: Ruan Xiaoguang,Zhan Zhaojie,Yan Zhaojin, et al. Digital bathymetric model fusion of offshore waters around China’s coastline based on global bathymetry data[J]. Haiyang Xuebao,2024, 46(7):16–28 doi: 10.12284/hyxb2024062

基于全球测深数据的中国海岸线周边海域数字水深模型融合

doi: 10.12284/hyxb2024062
基金项目: 浙江省基础公益研究计划项目(LZJWY22E090002);国家自然科学基金(42201451);中国博士后科学基金面上项目(2022M723379);南浔青年学者项目(RC2024021062);浙江省社会科学界联合会研究课题(2024N085);江苏省双创博士项目(JSSCBS20221523)。
详细信息
    作者简介:

    阮晓光(1990—),男,河南省周口市人,主要从事海底地形建模与分析研究。E-mail:ruanxg@zjweu.edu.cn

    通讯作者:

    闫兆进(1991—),男,山东省济宁市人,副教授,主要从事地学大数据挖掘与空间分析建模研究。E-mail:yanzhaojin@cumt.edu.cn

  • 中图分类号: P229.1;P28

Digital bathymetric model fusion of offshore waters around China’s coastline based on global bathymetry data

  • 摘要: 数字水深模型(Digital Bathymetric Models,简称“DBMs”),是近海工程建设、资源开发、环境保护等领域的重要基础地理信息数据。现有全球公开DBMs产品如GEBCO(The General Bathymetric Chart of the Oceans)、SRTM(The Shuttle Radar Topography Mission)、ETOPO(Earth Topography)等在不同海域的数据类型、数据来源和产品精度均存在差异。为利用全球测深数据和DBMs产品重建中国近海水深模型,本文提出一种基于水深分区的加权融合重建框架。首先,从5个维度(整体精度、不同水深、航线剖面、地理分区、局部细节)对比分析6种常用DBMs产品的可靠性和适用性;然后,顾及水深和地形特征对研究区进行分割和分区,并选取分区内最优DBMs产品,以最小误差为约束进行最优加权融合;最后,对融合结果进行实测值恢复、平滑滤波等后处理,形成中国海岸线周边近海海域15″分辨率高精度无缝水深模型。结果表明,融合结果相比SRTM30_PLUS、GEBCO_2022、SRTM15_V2.5.5和ETOPO_2022均方根误差降低了27%、14%、14%和13%,地形细节也得到保留,证明了该融合框架的可行性,可为多数据集大规模海底地形的融合重建和及时更新提供参考。
  • 海洋长期以来被视为人类探索世界最耗时的领域之一,全球海洋面积约为3.62亿km2,占世界面积的约71%,拥有不可估量的战略价值与意义[12]。我国陆架边缘海域约300万km2,占全球陆架海的12.5%,大陆岸线绵延1.8万km,同时还拥有1.4万km的岛屿海岸线[3]。水深是反映海底地形起伏的基本要素,不仅对人类认识和利用海洋起着关键作用,也是进行科学研究和制定海洋政策的重要依据[4]。数字水深模型(Digital Bathymetric Models,简称“DBMs”)是一种反映海洋水深变化描述海底地貌特征的三维数字模型,可用于海底地貌特征以及成因、海洋资源、海洋渔业等诸多领域研究,还在近海海底工程勘察设计、海上导航、海底电缆建设等工程领域具有广泛的应用[56]。随着船测水深数据、海洋重力信息的不断细化和完善,多个国家和组织通过数据融合建成了一些全球公开DBMs产品,如 GEBCO(The General Bathymetric Chart of the Oceans)、SRTM(The Shuttle Radar Topography Mission)等[79], 这些数据集为地学相关研究提供了重要支撑。然而,全球测深数据存在来源、格式和精度等多方面差异,构建高精度数字水深模型面临困难。此外,受重力和测高卫星轨道分辨率的局限,短时间内无法提供更高分辨率、更高精度的重力异常和垂直重力梯度,导致DBMs产品更新速度慢、空间分辨率低(最高约500 m),以南海为例重力反演海底地形精度在100 m以上[10]。因此,充分利用多源水深数据和DBMs产品互补优势,及时更新和优化全球或局部海域水深模型尤为重要,已经成为研究热点。

    近年来,在海洋测深数据融合方面,常用方法是经典的“Remove-Restore”法,利用差值栅格恢复局部的高分辨率数据,实现融合模型最小扰动情况下的平滑过渡,增强局部地形细节[1112]。演变而来的“合并–融合”方法融合多(单)波束声呐测深、海图水深与数字水深产品,亦被用于填补数据空洞,修复空洞区高分辨率海底地形,从而保留海底地形特征,已经成为生产全球大尺度高分辨率海底地形的主流手段[13]。海洋数字水深模型产品融合的研究与陆地DEM(Digital elevation model)类似,按照原理可分为频率域融合、稀疏表达融合、空间域融合等[5, 1416]。其中,直接平均或加权平均方法,简单易行、效率高,且融合效果好,但如何确定最优融合权重是避免融合结果出现伪影和边缘效应的关键[14]。以往研究中权重图通常由源DEM数据与参考数据间的高程误差图(Height error maps,HEMs)估算得到,该方法在缺乏可靠先验数据的海域并不可行,且未意识到分析影响权重图最相关地形参数的重要性,这限制了该方法的广泛适用性[14]

    针对多源水深数据和数字水深模型融合方面研究的不足,本文提出一种基于水深分区的加权融合重建框架,顾及水深和地形特征以最小误差为约束进行最优加权融合,实现中国海岸线周边近海海域高精度数字水深模型的及时更新和优化。

    本文以中国海岸线向外延伸150 km缓冲区的近海海域作为研究区域,如图1a所示,纬度范围在11.676°N 至41.047°N之间,经度范围在105.620°E 至130.355°E之间,国界底图来源于自然资源部标准地图服务系统(审图号:GS(2019)1822号)。

    图  1  a. 研究区;b. 测深数据分布
    该图基于自然资源部标准地图服务网站下载的审图号为GS(2019)1822号的标准地图制作,底图无修改
    Figure  1.  a. Study area ; b. distribution of bathymetry data
    Note: The map is produced based on the standard map with review number GS (2019) 1822 downloaded from the website of the Standard Map Service of the Ministry of Natural Resources of the People’s Republic of China, with no modifications to the base map

    船测水深数据主要取自NCEI(National Centers for Environmental Information)的全球海洋航迹地球物理数据库(http://www.ngdc.noaa.gov/mgg/geodas/geodas.html),该数据库收纳了全球范围内不同研究机构的实测水深数据。本文获取1950年至今所有航次采集的船测单(多)波束声呐测深数据,数据分布如图1b所示。鉴于该海域多波束声呐测深数据分布较少,本研究使用的数据主要是单波束声呐测深数据为主。

    海图作为地图学的一个分支,是以海洋为描绘对象的地图,海部要素包括水深、等深线、海底地质等,可用于支撑海上航行保障、资源开发、海洋划界和海洋军事活动等。本文从电子海图(Electronic Navigational Charts,简称“ENCs”)读取的水深值作为海岸线NCEI船测水深缺失的补充。该ENCs数据按照IHO S-57格式进行收集和汇编,由通过国际海道测量组织(IHO)认证、中国人民解放军海军海道测量局(CNHO)授权的数据服务商提供,其平面基准为 WGS-84 坐标系统,大部分深度在100 m以浅,比例尺为 1∶25 万,共计27幅,数据分布如图1b所示。

    选取6种不同分辨率、不同版本的常用DBMs产品作为融合数据源,其空间分辨率、空间范围、空间基准、数据来源等信息见表1,预处理后的模型结果如图2所示。其中,TOPO指Global Seafloor Topography;DTU指Technical University of Denmark;ETOPO指Earth Topography;MSL指Mean Sea Level(平均海平面);SIO指Scripps Institution of Oceanography。

    表  1  DBMs数据集概况
    Table  1.  An overview of DBMs datasets
    数据集 更新时间 组织机构 国家 分辨率 空间范围 水平基准垂直基准 数据来源
    GEBCO_2022 2022 The Nippon
    Foundation-
    GEBCO
    UK-Japan 15″ 179° 59' 52.5"W~0°
    179° 59' 52.5"E;
    89° 59' 52.5''N~
    89° 59' 52.5''S
    WGS84 MSL https://download.gebco.net/
    SRTM30_PLUS 2014 SIO USA 30″ 经度180°~0°~180°;
    90°N~90°S
    WGS84 MSL https://topex.ucsd.edu/pub/srtm30_plus/srtm30/grd/
    SRTM15_V2.5.5 2023 SIO USA 15″ 经度180°~0°~180°;
    90°N~90°S
    WGS84 MSL https://topex.ucsd.edu/pub/srtm15_plus/
    TOPO_25.1 2023 SIO USA 1' 经度180°~0°~180°;
    90°N~90°S
    WGS84 MSL https://topex.ucsd.edu/pub/global_topo_1min/
    DTU10 2010 DTU Space Denmark 1' 经度180°~0°~180°;
    90°N~90°S
    WGS84 MSL https://ftp.space.dtu.dk/pub/DTU10/1_MIN/
    ETOPO_2022 2023 NCEI USA 15″ 经度180°~0°~180°;
    90°N~90°S
    WGS84 MSL https://www.ncei.noaa.gov/products/etopo-global-relief-model/
    下载: 导出CSV 
    | 显示表格
    图  2  研究区DBMs. a. GEBCO_2022; b. SRTM30_PLUS; c. SRTM15_V2.5.5; d. TOPO_25.1; e. DTU10; f. ETOPO_2022
    Figure  2.  DBMs in the study area. a. GEBCO_2022; b. SRTM30_PLUS; c. SRTM15_V2.5.5; d. TOPO_25.1; e. DTU10; f. ETOPO_2022

    其中,GEBCO_2022数据产品是the Nippon Foundation-GEBCO Seabed 2030 计划的第4个版本,基础主要来源于SRTM15_PLUS,最新区域数据和重力模型的使用,使得GEBCO_2022对海底地形的描述更加详细;SRTM15_V2.5.5汇聚了3亿多条新施测的水深航迹线,其15″分辨率下的数据采样率超过11%;SRTM30_PLUS分辨率为30″,其海洋水深部分由CryoSat-2和Jason-1重力模型反演得到,并利用2.98亿个测深数据对重力与地形的比值进行校正;TOPO_25.1模型融合了Cryosat-2、Jason-1和Envisat多源卫星测高数据,由V29全球重力异常模型反演得到,其1′的分辨率影响了地形细节。DTU10模型由ERS-1、GEOSAT等卫星高度计重力异常反演得到,并集成GEBCO早期数据集,因此存在统计偏差,对海底地形的表达不够清晰,在平面上呈不连续的阶梯状分布。ETOPO_2022采用了最新的ICESat-2数据和其他数据源,其在反映重力异常长波信息的大陆架、海盆等区域精度高,但浅滩、暗礁等浅水区存在不确定性。

    几乎所有的DBMs都是由卫星重力信息(重力异常或垂直重力梯度异常)反演得到,在不同的海底区域存在的重力信息长短波分离误差、密度差常数误差等因素,直接导致反演水深精度的不确定性,例如海岸附近和深海区域[512]。所选用的6个DBMs各有优缺点,不同区域的质量也不尽相同,因此有必要对模型进行融合和优化。首先,对所有数据和产品进行预处理,包括格式转换、平面坐标和水深基准统一、分辨率统一、数据清洗等步骤;然后,对DBMs进行基于水深分区的加权融合,包括地形分割和分区、模型对比与优选、分区加权融合、融合模型后处理等;最后,对模型融合结果进行多维评估。技术流程如图3所示。

    图  3  模型融合框架
    Figure  3.  The model fusion framework

    (1)测深数据格式转换

    NCEI全球海洋航迹地球物理原始数据以M77T方式存储,包含ASCII文本格式的测深数据,该格式包括标题记录和数据记录。其中,标题记录提供每个调查对象的元数据,数据记录提供测量机构的ID、测量日期、时区、测点水深以及经纬度坐标。本文利用通用制图工具GMT(Generic Mapping Tools)集成的MGD7/MGD77T 模块,将M77T格式的测深数据合并转换为XYZ ASCII数据格式[17]

    (2)数学基准统一

    数学基础的统一包括平面坐标、水深基准以及分辨率的统一。

    空间基准的统一包括水平坐标基准和垂直基准统一。本文利用Safe Software 软件集成的FME数据格式转换模块(Feature Manipulation Engine),将NCEI船测水深数据的水平基准转换为WGS-84坐标系,同时将数据转换为Esri shapefile文件格式,垂直基准为平均海平面(MSL)。同样,S-57格式的电子海图原始数据也转为Esri shapefile格式。

    此外,6种水深产品数据类型不尽相同,包括netCDF、GeoTiff、Esri ASCII等栅格类型,本文将其统一为GeoTiff格式,并投影为Asia_North_Albers_Equal_Area_Conic等面积投影坐标系。考虑到本研究待编辑和融合的水深产品最高空间分辨率为15"(约450 m),融合前需将6个产品的空间分辨率统一为 15″,利用反距离权重插值法(IDW)将分辨率为1′和30″的数据产品分辨率插值为15″。船测水深点的原始间距最小不足1 m,空间分辨率过高,需要对原始船测水深点进行格网化,保证每个水深栅格(15″×15″)仅对应一个水深值。若一个栅格范围存在多点,则栅格值用多个水深点的平均值代替。经过预处理的6个水深产品如图2所示。

    (3)数据清洗

    为降低粗差等对精度估计结果的影响,本文利用现有应用最广的水深产品之一SRTM30_PLUS数据集水深值作为参考数据,按国际海道测量组织IHO S-44 第5 版有关规定剔除无效值(corr_= null,或corr_≥0 m),并根据pauta准则对所有有效航次数据进行异常值检验和剔除[18]。具体方法参考SRTM30_PLUS数据集编辑和生产过程[9],共剔除8860个水深异常值,预处理后的数据如图1b所示。

    为验证处理后的NCEI船测水深数据可靠性,构建了NCEI船测水深数据与6个DBMs水深值点位之间一一对应的空间相关矩阵(表2),结果证明相关系数均接近于1,反映出几种数据源互相存在较高的总体相关性。进一步,在船测水深数据与海图数据存在重叠的区域进行空间插值建立水深格网,再利用船测水深提取格网水深值进行一致性验证,结果如图4所示,两者拟合较好,相关系数达到0.99,证明处理后的NCEI船测数据具有一定可靠性和一致性。

    表  2  NCEI船测水深数据与DBMs的空间相关矩阵
    Table  2.  The spatial correlation matrix between NCEI ship-measured bathymetry data and DBMs
    船测水深GEBCO_2022SRTM30_PLUSSRTM15_V2.5.5TOPO_25.1DTU10ETOPO_2022
    船测水深10.99900.99860.99900.99800.99200.9990
    GEBCO_20220.999010.99930.99960.99860.99250.9998
    SRTM30_PLUS0.99860.999310.99920.99910.99280.9991
    SRTM15_V2.5.50.99900.99960.999210.99870.99240.9996
    TOPO_25.10.99800.99860.99910.998710.99330.9985
    DTU100.99200.99250.99280.99240.993310.9924
    ETOPO_20220.99900.99980.99910.99960.99850.99241
    下载: 导出CSV 
    | 显示表格
    图  4  NCEI船测水深数据可靠性验证
    Figure  4.  Reliability verification of NCEI ship-measured bathymetry data

    首先,对基于实测水深数据对中国近海DBMs进行可靠性和适用性评估;然后,顾及水深和地形特征对研究区进行分割和分区,并选取最优DBMs产品进行最优加权融合;最后,对融合结果进行实测值恢复、平滑滤波等后处理,得到最终的融合DBM(图3)。

    (1)地形分割和分区

    研究表明,重力水深反演精度与地形特征高度相关,缓坡区明显高于海山区[5]。因此,为顾及不同模型对融合模型的贡献,有必要对研究区进行分割和分区。本文基于水深范围和海底地形特征,通过多尺度分割选定最优分割尺度,该方法从像素层开始,计算不同尺度参数下对象内均质性结果的局部方差,以局部方差变化率最大值所相应的分割尺度为最优尺度,从而保证分割对象内的均质性和对象间的异质性均得到最大化。最终将研究区划定5个地理子区域,空白区域为陆地,如图5所示。

    图  5  分割分区结果
    Figure  5.  The segmentation and partitioning results

    (2)模型对比与优选

    利用处理后的实测数据对中国近海DBMs产品质量进行定量评估,分析各自的可靠性和适用性,作为分区内选择待融合模型的参考。本文采用平均绝对误差(mean absolute error,简称“MAE”)、标准差(standard deviation,简称“SD”)、均方根误差(root mean square error,简称“RMSE”)等3个误差指标,用于定量分析。其中,MAE能表征在没有异常值情况下的水深值精度鲁棒性;SD可用于用来衡量一组水深值自身的离散程度;RMSE则能更准确地反映存在异常值情况下的水深值误差分布。

    $$ {\mathrm{MAE}} = \frac{1}{n}\sum\limits_{i = 1}^n {\left| {{e_i}} \right|}, $$ (1)
    $$ {\mathrm{SD}} = \sqrt {\frac{1}{n}\sum\limits_{i = 1}^n {{{({e_i} - \bar e)}^2}} }, $$ (2)
    $$ {\mathrm{RMSE}} = \sqrt {\frac{1}{n}\sum\limits_{i = 1}^n {{e_i}^2} }, $$ (3)

    式中,ei为原始DBMs的水深值与实测数据水深值作差得到的误差;$ \bar e $为误差均值;n为参与计算的实测水深值点数。

    此外,采用决定系数R2表征DBMs水深值与实测水深值的拟合程度,R²的值越接近1,说明二者的拟合程度越好。

    $$ {R^2} = 1 - \frac{{\displaystyle\sum\limits_{i = 1}^n {{e_i}^2} }}{{\displaystyle\sum\limits_{i = 1}^n {{{({d_i} - \bar d)}^2}} }}, $$ (4)

    式中,di为原始DBMs的水深值;$ \bar d $为原始DBMs的水深平均值。

    (3)分区加权融合

    加权融合是一种简单实用的融合方法,可以抑制原始DBM中的噪声[14]。本文以最小误差(RMSE)为约束条件,确定每个分区内待融合模型的最优融合权重。通过遍历权值分配确定多个水深产品的最优权值,再利用线性加权平均,计算最优融合结果,具体步骤如图6所示。

    图  6  最优融合权重选取步骤
    Figure  6.  The selection steps for the optimal fusion weights
    $${\mathrm{ Bathymetr}}{y_i} = {\mathrm{DB}}{{\mathrm{M}}_1} \times {w_i} + {\mathrm{DB}}{{\mathrm{M}}_2} \times (1 - {w_i}), $$ (5)

    式中,DBM1和DBM2分别为待融合数字水深模型;wi为DBM1对应的融合权重;Bathymetryi为经过加权融合的数字水深模型水深结果;研究区分区权重分配结果见表3

    表  3  最优融合权重分配表
    Table  3.  Distribution table for the optimal fusion weights
    水深/m RMSE(m) 权重分配
    融合前 融合后 GEBCO_2022 ETOPO_2022 SRTM15_V2.5.5
    0~−200 20.324 20.130 0.1 0.9 0
    −200~−500 61.327 61.325 0.1 0.9 0
    −500~−2000 70.394 69.563 0.1 0.9 0
    2000~−4000 90.683 84.767 0.5 0 0.5
    4000~−8000 78.804 78.730 0.2 0.8 0
    下载: 导出CSV 
    | 显示表格

    (4)融合模型后处理

    不同DBMs在水深分割边界处略有差异,导致融合模型接缝处有空白,为了形成质量更好的融合产品,有必要对模型进行异常值剔除、平滑滤波等后处理,因此本文提出一种实测点恢复和自适应滤波的模型后处理方法。由于研究区实测数据主要为离散的单波束水深数据,经典的“Remove-Restore”后处理方法并不适用[9],提出通过实测点恢复和基于邻域统计的自适应平滑滤波方法,对上一步融合后的地形模型进行优化。将实测点转化为规则格网,与融合模型进行镶嵌可以实现融合模型的实测点恢复,再通过异常值检测与滤波,实现空洞填充,减少地形异常值。具体步骤如下:①利用自然邻域法将Esri shapefile的格式实测数据转为15″×15″分辨率的规则格网;②移去其在融合模型中的对应栅格,二者镶嵌恢复融合模型的实测点;③提取每个像素中心周围3×3窗口大小的邻域补丁块,计算邻域像素与中心像素之间的差值,若差值大于3 倍标准差,则将该中心像素判断为离群值(异常值);④利用邻域像素的均值替代异常值,同时实现原始地形产品中的空洞填充,从而达到滤除异常值的目的(图7)。

    图  7  模型融合后处理
    a. 15″×15″水深栅格;b. 镶嵌栅格;c. 邻域统计结果;d. 最终融合结果
    Figure  7.  Post-processing of model fusion
    a. 15″×15″ bathymetric raster; b. mosaiced raster; c. neighbourhood statistics results; d. final fusion result

    (1)整体分析

    引入实测数数据,以6个DBMs减去实测值得到的差异,进行定量评价,结果如图8所示。整体表现最好的GEBCO_2022、SRTM15_V2.5.5、SRTM30_PLUS和ETOPO_2022等DBMs产品,其误差范围主要在−100~100 m之间,占80%以上。从误差分布频率直方图(图8可见),差异总体呈现正态分布,其中DTU10模型以外,其余模型水深误差都趋于正值,说明模型水深比实测值更浅。其中,GEBCO_2022、SRTM15_V2.5.5和 ETOPO_2022的水深误差分布较为集中,R2接近于1,与实测值拟合结果最好。ETOPO_2022的RMSE、SD以及MAE最小,整体精度较高,精度鲁棒性好。6个DBMs模型中有2个SRTM模型,SRTM15_V2.5.5的分辨率为15″,SRTM30_PLUS的分辨率为30″,但是从整体结果来看,两者的精度相差不大,SRTM15_V2.5.5的各项精度指标略好于SRTM30_PLUS。DTU10误差主要分布在−500~500 m之间,占总体差异的85%以上,误差较为离散,精度最差,R2结果最低(0.99)。

    图  8  DBMs与实测水深值的误差分布
    正值表示 DBM 水深比实测值浅,a. GEBCO_2022;b. SRTM30_PLUS;c. SRTM15_V2.5.5;d. TOPO_25.1;e. DTU10;f. ETOPO_2022
    Figure  8.  Distribution of differences bteween DBMs and measured bathymetry values
    A positive value indicates that the DBM bathymetry is shallower than that of the measured values, a. GEBCO_2022; b. SRTM30_PLUS; c. SRTM15_V2.5.5; d. TOPO_25.1; e. DTU10; f. ETOPO_2022

    (2)航线对比分析

    选取船测航线的实测数据对6个模型的地形细节进行分析,选取EW9509与25874两条航线部分进行对比(图9)。

    图  9  航线分析
    a. 航线位置;b. 25874航线剖面;c. EW9509航线剖面(黑点代表实测水深)
    Figure  9.  Route analysis
    a. Route location; b. profile route 25874; c. profile route EW9509 (the black spots represent measured bathymetry)

    对于航线EW9509,GEBCO_2022、SRTM30_PLUS、SRTM15_V2.5.5、TOPO_25.1、DTU10与ETOPO_2022的R2分别为0.97、0.97、0.98、0.93、0.66和0.96。对于航线25874,GEBCO_2022、SRTM30_PLUS、SRTM15_V2.5.5、TOPO_25.1、DTU10与ETOPO_2022的R2分别为0.70、0.76、0.66、0.63、0.58和0.70。通过计算每个航迹点位置DBMs水深值与实测值之差的最大值(MAX_Difference),证明DBMs的绝对差异与水深关系不大,最大差异主要集中在海底地形变化剧烈的海域。明显可见在地形起伏较缓处,除DTU10,5个模型差异不大,剖面线一致性较好。DTU10模型整体曲线平滑,细节丢失严重。

    (1)总体结果分析

    总体融合结果如图10所示。与原始DBMs相比,融合的DBM总体精度得到了提高,差异范围更小,主要集中在−100~100 m,RMSE相比,SRTM30_PLUS、GEBCO_2022、SRTM15_V2.5.5与ETOPO_2022,降低了27%、14%、14%和13%。融合后的模型其R2拟合结果最好,误差分布比原始模型更为集中,更加接近1。

    图  10  融合结果
    a. 融合模型可视化;b. 融合模型整体精度评估
    Figure  10.  Fusion results
    a. Fusion model visualization; b. overall accuracy evaluation of fusion model

    (2)不同水深对比

    不同水深范围融合前后结果对比如图11所示。与原始DBMs相比,融合DBM在各个水深处的RMSE有所降低,0~−200 m、−500~−2000 m、−2000~−4000 m以及−4000~−8000 m分别降低了1.34 m、0.83 m、5.91 m以及0.07 m。其中在−2000 m~−4000 m处效果提升最好,但在−4000 m~−8000 m的深海区域,该融合模型精度虽有提升,但是提升效果有限,也体现出公开水深模型在不同水深区域存在一定不确定性。

    图  11  不同深度融合结果对比
    Figure  11.  Comparison of fusion results at different depths

    (3)不同海域分析

    不同海域(黄海、南海、东海)融合前后结果对比见表4。实验结果表明,不同海域内的精度指标均有提高。黄海的MAE、SD与RMSE等误差指标分别降低4.5%、17%、17%。南海3个误差指标分别降低33%、24%与24%,R2提高了0.05。东海3个误差指标分别降低60%、10.1%、8.3%,R2提高了0.04。以上数据显示各个海域的提升效果不同,南海精度提升效果最明显,东海次之,黄海精度提升效果最低。

    表  4  不同海域融合前后精度对比
    Table  4.  Accuracy comparison before and after fusion of different sea areas
    地理分区 精度/m R2
    最大值 最小值 平均值 MAE SD RMSE
    黄海 融合前 795.88 −66.64 0.18 1.33 9.19 9.22 0.99
    融合后 675.91 −92.17 0.11 1.27 7.59 7.59 0.99
    南海 融合前 949.08 −948.11 0.09 31.34 86.42 86.42 0.93
    融合后 3601.89 1056.83 1.07 20.88 65.20 65.21 0.98
    东海 融合前 1180.18 −949.14 −1.40 26.94 67.02 67.04 0.89
    融合后 1196.48 1522.43 −0.56 24.77 61.46 61.46 0.93
    下载: 导出CSV 
    | 显示表格

    (4)局部地形分析

    选择西沙群岛和琉球群岛局部及其邻近海域,分别代表大型珊瑚礁发育的浅海复杂地形海域,以及地形起伏大的深浅海过渡海域,对融合前后的DBMs进行局部细节对比分析,以及地形剖面分析。局部细节对比结果如图12图13所示,其中地形细节和异常值等特征用黑色虚线框标出。实验结果表明,在西沙群岛海域,融合后的DBM与SRTM15_V2.5.5和GEBCO_2022结果十分接近,水深梯度分布类似,海平面以上的珊瑚礁大小和位置一致(如图12黑色虚线框A),这也验证了SRTM15_V2.5.5和GEBCO_2022的原始实测数据来源相似。在地形起伏较大的琉球群岛海域,ETOPO_2022与GEBCO_2022保留了深海与浅海过渡区域的更多细节,而ETOPO_2022在深海区域整体深度大于其他模型(如图13黑色虚线框B)。整体来看,GEBCO_2022、ETOPO_2022与融合后模型地形细节较丰富;DTU10模型轮廓相对更为圆滑,虽然整体地形与其余模型类似,但水深细节丢失严重,分辨率不足的劣势凸显;TOPO_25.1整体水深梯度与其余DBMs接近,但分辨率过低,使得地形细节模糊,尤其在在水深变化剧烈处出现大量异常值和明显噪声。

    图  12  西沙群岛结果验证
    a. GEBCO_2022;b. SRTM30_PLUS;c. SRTM15_V2.5.5;d. TOPO_25.1;e. DTU10;f. ETOPO_2022;g. 融合模型
    Figure  12.  The result verification of Xisha Islands
    a. GEBCO_2022; b. SRTM30_PLUS; c. SRTM15_V2.5.5; d. TOPO_25.1; e. DTU10; f. ETOPO_2022; g. Fusion model
    图  13  琉球群岛结果验证
    a. GEBCO_2022;b. SRTM30_PLUS;c. SRTM15_V2.5.5;d. TOPO_25.1;e. DTU10;f. ETOPO_2022;g. 融合模型
    Figure  13.  The result verification of Ryukyu Islands
    a. GEBCO_2022; b. SRTM30_PLUS; c. SRTM15_V2.5.5; d. TOPO_25.1; e. DTU10; f. ETOPO_2022; g. Fusion model

    进一步由西向东、由北向南设置1000个剖面点,绘制的西沙群岛及琉球群岛局部地形剖面如图14图15所示。实验结果表明,DBMs整体剖面相差不大,但细节存在差异。其中GEBCO_2022、SRTM15_V2.5.5、ETOPO_2022与融合模型在细节以及整体地形变化上最为接近且地形细节清晰。DTU10剖面过渡圆滑,但地形细节丢失严重。

    图  14  西沙群岛局部剖面
    a. GEBCO_2022;b. SRTM30_PLUS;c. SRTM15_V2.5.5;d. TOPO_25.1;e. DTU10;f. ETOPO_2022;g. 融合模型
    Figure  14.  Local section of Xisha Islands
    a. GEBCO_2022; b. SRTM30_PLUS; c. SRTM15_V2.5.5; d. TOPO_25.1; e. DTU10; f. ETOPO_2022; g. Fusion model
    图  15  琉球群岛局部剖面
    a. GEBCO_2022;b. SRTM30_PLUS;c. SRTM15_V2.5.5;d. TOPO_25.1;e. DTU10;f. ETOPO_2022;g. 融合模型
    Figure  15.  Local section of Ryukyu Islands
    a. GEBCO_2022; b. SRTM30_PLUS; c. SRTM15_V2.5.5; d. TOPO_25.1; e. DTU10; f. ETOPO_2022; g. Fusion model

    本文提出一种基于水深分区的加权融合重建框架,利用全球测深数据和全球公开的6种DBMs模型,实现了中国近海高精度无缝水深模型的更新和优化,得出以下主要结论:

    (1)从多个维度对全球公开的6种DBMs模型进行可靠性和适用性分析,分析结果表明GEBCO_2022、ETOPO_2022和SRTM30_PLUS精度最优,但ETOPO_2022的整体水深大于其余模型。

    (2)构建的中国近海15″分辨率高精度无缝水深模型,其精度指标优于原始DBMs产品,与SRTM30_PLUS、GEBCO_2022、SRTM15_V2.5.5和ETOPO_2022相比,RMSE分别降低了27%、14%、14%和13%,地形细节也得到保留,地形表达更加真实丰富。

    (3)基于水深分区的加权融合方法的可行性得到验证,可用于多数据集大规模海底地形的融合重建和及时更新。

  • 图  1  a. 研究区;b. 测深数据分布

    该图基于自然资源部标准地图服务网站下载的审图号为GS(2019)1822号的标准地图制作,底图无修改

    Fig.  1  a. Study area ; b. distribution of bathymetry data

    Note: The map is produced based on the standard map with review number GS (2019) 1822 downloaded from the website of the Standard Map Service of the Ministry of Natural Resources of the People’s Republic of China, with no modifications to the base map

    图  2  研究区DBMs. a. GEBCO_2022; b. SRTM30_PLUS; c. SRTM15_V2.5.5; d. TOPO_25.1; e. DTU10; f. ETOPO_2022

    Fig.  2  DBMs in the study area. a. GEBCO_2022; b. SRTM30_PLUS; c. SRTM15_V2.5.5; d. TOPO_25.1; e. DTU10; f. ETOPO_2022

    图  3  模型融合框架

    Fig.  3  The model fusion framework

    图  4  NCEI船测水深数据可靠性验证

    Fig.  4  Reliability verification of NCEI ship-measured bathymetry data

    图  5  分割分区结果

    Fig.  5  The segmentation and partitioning results

    图  6  最优融合权重选取步骤

    Fig.  6  The selection steps for the optimal fusion weights

    图  7  模型融合后处理

    a. 15″×15″水深栅格;b. 镶嵌栅格;c. 邻域统计结果;d. 最终融合结果

    Fig.  7  Post-processing of model fusion

    a. 15″×15″ bathymetric raster; b. mosaiced raster; c. neighbourhood statistics results; d. final fusion result

    图  8  DBMs与实测水深值的误差分布

    正值表示 DBM 水深比实测值浅,a. GEBCO_2022;b. SRTM30_PLUS;c. SRTM15_V2.5.5;d. TOPO_25.1;e. DTU10;f. ETOPO_2022

    Fig.  8  Distribution of differences bteween DBMs and measured bathymetry values

    A positive value indicates that the DBM bathymetry is shallower than that of the measured values, a. GEBCO_2022; b. SRTM30_PLUS; c. SRTM15_V2.5.5; d. TOPO_25.1; e. DTU10; f. ETOPO_2022

    图  9  航线分析

    a. 航线位置;b. 25874航线剖面;c. EW9509航线剖面(黑点代表实测水深)

    Fig.  9  Route analysis

    a. Route location; b. profile route 25874; c. profile route EW9509 (the black spots represent measured bathymetry)

    图  10  融合结果

    a. 融合模型可视化;b. 融合模型整体精度评估

    Fig.  10  Fusion results

    a. Fusion model visualization; b. overall accuracy evaluation of fusion model

    图  11  不同深度融合结果对比

    Fig.  11  Comparison of fusion results at different depths

    图  12  西沙群岛结果验证

    a. GEBCO_2022;b. SRTM30_PLUS;c. SRTM15_V2.5.5;d. TOPO_25.1;e. DTU10;f. ETOPO_2022;g. 融合模型

    Fig.  12  The result verification of Xisha Islands

    a. GEBCO_2022; b. SRTM30_PLUS; c. SRTM15_V2.5.5; d. TOPO_25.1; e. DTU10; f. ETOPO_2022; g. Fusion model

    图  13  琉球群岛结果验证

    a. GEBCO_2022;b. SRTM30_PLUS;c. SRTM15_V2.5.5;d. TOPO_25.1;e. DTU10;f. ETOPO_2022;g. 融合模型

    Fig.  13  The result verification of Ryukyu Islands

    a. GEBCO_2022; b. SRTM30_PLUS; c. SRTM15_V2.5.5; d. TOPO_25.1; e. DTU10; f. ETOPO_2022; g. Fusion model

    图  14  西沙群岛局部剖面

    a. GEBCO_2022;b. SRTM30_PLUS;c. SRTM15_V2.5.5;d. TOPO_25.1;e. DTU10;f. ETOPO_2022;g. 融合模型

    Fig.  14  Local section of Xisha Islands

    a. GEBCO_2022; b. SRTM30_PLUS; c. SRTM15_V2.5.5; d. TOPO_25.1; e. DTU10; f. ETOPO_2022; g. Fusion model

    图  15  琉球群岛局部剖面

    a. GEBCO_2022;b. SRTM30_PLUS;c. SRTM15_V2.5.5;d. TOPO_25.1;e. DTU10;f. ETOPO_2022;g. 融合模型

    Fig.  15  Local section of Ryukyu Islands

    a. GEBCO_2022; b. SRTM30_PLUS; c. SRTM15_V2.5.5; d. TOPO_25.1; e. DTU10; f. ETOPO_2022; g. Fusion model

    表  1  DBMs数据集概况

    Tab.  1  An overview of DBMs datasets

    数据集 更新时间 组织机构 国家 分辨率 空间范围 水平基准垂直基准 数据来源
    GEBCO_2022 2022 The Nippon
    Foundation-
    GEBCO
    UK-Japan 15″ 179° 59' 52.5"W~0°
    179° 59' 52.5"E;
    89° 59' 52.5''N~
    89° 59' 52.5''S
    WGS84 MSL https://download.gebco.net/
    SRTM30_PLUS 2014 SIO USA 30″ 经度180°~0°~180°;
    90°N~90°S
    WGS84 MSL https://topex.ucsd.edu/pub/srtm30_plus/srtm30/grd/
    SRTM15_V2.5.5 2023 SIO USA 15″ 经度180°~0°~180°;
    90°N~90°S
    WGS84 MSL https://topex.ucsd.edu/pub/srtm15_plus/
    TOPO_25.1 2023 SIO USA 1' 经度180°~0°~180°;
    90°N~90°S
    WGS84 MSL https://topex.ucsd.edu/pub/global_topo_1min/
    DTU10 2010 DTU Space Denmark 1' 经度180°~0°~180°;
    90°N~90°S
    WGS84 MSL https://ftp.space.dtu.dk/pub/DTU10/1_MIN/
    ETOPO_2022 2023 NCEI USA 15″ 经度180°~0°~180°;
    90°N~90°S
    WGS84 MSL https://www.ncei.noaa.gov/products/etopo-global-relief-model/
    下载: 导出CSV

    表  2  NCEI船测水深数据与DBMs的空间相关矩阵

    Tab.  2  The spatial correlation matrix between NCEI ship-measured bathymetry data and DBMs

    船测水深GEBCO_2022SRTM30_PLUSSRTM15_V2.5.5TOPO_25.1DTU10ETOPO_2022
    船测水深10.99900.99860.99900.99800.99200.9990
    GEBCO_20220.999010.99930.99960.99860.99250.9998
    SRTM30_PLUS0.99860.999310.99920.99910.99280.9991
    SRTM15_V2.5.50.99900.99960.999210.99870.99240.9996
    TOPO_25.10.99800.99860.99910.998710.99330.9985
    DTU100.99200.99250.99280.99240.993310.9924
    ETOPO_20220.99900.99980.99910.99960.99850.99241
    下载: 导出CSV

    表  3  最优融合权重分配表

    Tab.  3  Distribution table for the optimal fusion weights

    水深/m RMSE(m) 权重分配
    融合前 融合后 GEBCO_2022 ETOPO_2022 SRTM15_V2.5.5
    0~−200 20.324 20.130 0.1 0.9 0
    −200~−500 61.327 61.325 0.1 0.9 0
    −500~−2000 70.394 69.563 0.1 0.9 0
    2000~−4000 90.683 84.767 0.5 0 0.5
    4000~−8000 78.804 78.730 0.2 0.8 0
    下载: 导出CSV

    表  4  不同海域融合前后精度对比

    Tab.  4  Accuracy comparison before and after fusion of different sea areas

    地理分区 精度/m R2
    最大值 最小值 平均值 MAE SD RMSE
    黄海 融合前 795.88 −66.64 0.18 1.33 9.19 9.22 0.99
    融合后 675.91 −92.17 0.11 1.27 7.59 7.59 0.99
    南海 融合前 949.08 −948.11 0.09 31.34 86.42 86.42 0.93
    融合后 3601.89 1056.83 1.07 20.88 65.20 65.21 0.98
    东海 融合前 1180.18 −949.14 −1.40 26.94 67.02 67.04 0.89
    融合后 1196.48 1522.43 −0.56 24.77 61.46 61.46 0.93
    下载: 导出CSV
  • [1] Sandwell D T, Müller R D, Smith W H F, et al. New global marine gravity model from CryoSat-2 and Jason-1 reveals buried tectonic structure[J]. Science, 2014, 346(6205): 65−67. doi: 10.1126/science.1258213
    [2] 吴园涛, 段晓男, 沈刚, 等. 强化我国海洋领域国家战略科技力量的思考与建议[J]. 地球科学进展, 2021, 36(4): 413−420. doi: 10.11867/j.issn.1001-8166.2021.039

    Wu Yuantao, Duan Xiaonan, Shen Gang, et al. Thoughts and suggestions on strengthening the national strategic scientific and technological forces in the marine field of China[J]. Advances in Earth Science, 2021, 36(4): 413−420. doi: 10.11867/j.issn.1001-8166.2021.039
    [3] 张继红, 方建光, 唐启升. 中国浅海贝藻养殖对海洋碳循环的贡献[J]. 地球科学进展, 2005, 20(3): 359−365. doi: 10.3321/j.issn:1001-8166.2005.03.014

    Zhang Jihong, Fang Jianguang, Tang Qisheng. The contribution of shellfish and seaweed mariculture in China to the carbon cycle of coastal ecosystem[J]. Advances in Earth Science, 2005, 20(3): 359−365. doi: 10.3321/j.issn:1001-8166.2005.03.014
    [4] 陈义兰, 唐秋华, 刘晓瑜, 等. 多源水深数据融合的近海数字水深模型构建[J]. 海洋科学进展, 2021, 39(3): 461−469. doi: 10.3969/j.issn.1671-6647.2021.03.012

    Chen Yilan, Tang Qiuhua, Liu Xiaoyu, et al. Construction of offshore digital bathymetric model based on multi-source bathymetric data fusion[J]. Advances in Marine Science, 2021, 39(3): 461−469. doi: 10.3969/j.issn.1671-6647.2021.03.012
    [5] Ruan Xiaoguang, Cheng Liang, Chu Sensen, et al. A new digital bathymetric model of the South China Sea based on the subregional fusion of seven global seafloor topography products[J]. Geomorphology, 2020, 370: 107403. doi: 10.1016/j.geomorph.2020.107403
    [6] 董箭, 彭认灿, 张立华, 等. 滚动球变换的数字水深模型多尺度表达[J]. 地球信息科学学报, 2012, 14(6): 704−711. doi: 10.3724/SP.J.1047.2012.00704

    Dong Jian, Peng Rencan, Zhang Lihua, et al. Multi-scale representation of digital depth model based on rolling ball transform[J]. Journal of Geo-information Science, 2012, 14(6): 704−711. doi: 10.3724/SP.J.1047.2012.00704
    [7] 阳凡林, 沈瑞杰, 梅赛, 等. 联合重力异常和重力垂直梯度异常数据反演皇帝山海域海底地形[J]. 海洋学报, 2022, 44(12): 126−135.

    Yang Fanlin, Shen Ruijie, Mei Sai, et al. Inversion of seafloor topography in Emperor Seamount sea area by combined gravity anomalies and vertical gravity gradient anomalies data[J]. Haiyang Xuebao, 2022, 44(12): 126−135.
    [8] Weatherall P, Marks K M, Jakobsson M, et al. A new digital bathymetric model of the world's oceans[J]. Earth and Space Science, 2015, 2(8): 331−345. doi: 10.1002/2015EA000107
    [9] Becker J J, Sandwell D T, Smith W H F, et al. Global bathymetry and elevation data at 30 arc seconds resolution: SRTM30_PLUS[J]. Marine Geodesy, 2009, 32(4): 355−371. doi: 10.1080/01490410903297766
    [10] Fan Diao, Li Shanshan, Feng Jinkai, et al. A new global bathymetry model: STO_IEU2020[J]. Remote Sensing, 2022, 14(22): 5744. doi: 10.3390/rs14225744
    [11] Leon J X, Phinn S R, Hamylton S, et al. Filling the ‘white ribbon’ – a multisource seamless digital elevation model for Lizard Island, northern Great Barrier Reef[J]. International Journal of Remote Sensing, 2013, 34(18): 6337−6354. doi: 10.1080/01431161.2013.800659
    [12] Tozer B, Sandwell D T, Smith W H F, et al. Global bathymetry and topography at 15 arc sec: SRTM15+[J]. Earth and Space Science, 2019, 6(10): 1847−1864. doi: 10.1029/2019ea000658
    [13] 刘洋, 吴自银, 赵荻能, 等. MF多源测深数据融合方法及大洋水深模型构建[J]. 测绘学报, 2019, 48(9): 1171−1181. doi: 10.11947/j.AGCS.2019.20180495

    Liu Yang, Wu Ziyin, Zhao Dineng, et al. The MF method for multi-source bathymetric data fusion and ocean bathymetric model construction[J]. Acta Geodaetica et Cartographica Sinica, 2019, 48(9): 1171−1181. doi: 10.11947/j.AGCS.2019.20180495
    [14] Okolie C J, Smit J L. A systematic review and meta-analysis of digital elevation model (DEM) fusion: pre-processing, methods and applications[J]. ISPRS Journal of Photogrammetry and Remote Sensing, 2022, 188: 1−29. doi: 10.1016/j.isprsjprs.2022.03.016
    [15] 岳林蔚, 沈焕锋, 袁强强, 等. 基于深度置信网络的多源DEM点面融合模型[J]. 武汉大学学报(信息科学版), 2021, 46(7): 1090−1097.

    Yue Linwei, Shen Huanfeng, Yuan Qiangqiang, et al. A multi-source DEM point-surface fusion model based on deep belief network[J]. Geomatics and Information Science of Wuhan University, 2021, 46(7): 1090−1097.
    [16] Guan Liyi, Hu Jun, Pan Hongbo, et al. Fusion of public DEMs based on sparse representation and adaptive regularization variation model[J]. ISPRS Journal of Photogrammetry and Remote Sensing, 2020, 169: 125−134. doi: 10.1016/j.isprsjprs.2020.09.005
    [17] Wessel P. Tools for analyzing intersecting tracks: the x2sys package[J]. Computers & Geosciences, 2010, 36(3): 348−354. doi: 10.1016/j.cageo.2009.05.009
    [18] Tunalioglu N, Ocalan T, Dogan A H. Precise point positioning with GNSS raw measurements from an Android smartphone in marine environment monitoring[J]. Marine Geodesy, 2022, 45(3): 274−294. doi: 10.1080/01490419.2022.2027831
  • 加载中
图(15) / 表(4)
计量
  • 文章访问数:  239
  • HTML全文浏览量:  89
  • PDF下载量:  25
  • 被引次数: 0
出版历程
  • 收稿日期:  2023-07-14
  • 录用日期:  2024-08-12
  • 修回日期:  2024-01-24
  • 网络出版日期:  2024-08-15
  • 刊出日期:  2024-07-01

目录

/

返回文章
返回