1、第 45 卷第 1 期2024 年 1 月Vol.45,No.1Jan.2024遥 测 遥 控Journal of Telemetry,Tracking and CommandWebsite:基于SAR子孔径图像的浅海水下地形探测荣飞成1,2,孟俊敏1,纪永刚2(1 自然资源部第一海洋研究所 青岛 266061;2 中国石油大学(华东)海洋与空间信息学院 青岛 266580)摘要:沿海地区水下地形的精确测量是人类开展海洋活动的关键,合成孔径雷达(SAR)为浅海地形的探测提供了一种新的手段,其中TerraSAR-X 的聚束模式以较长的积分时间得到高分辨率的 SAR 数据,从而能较为精确地反演浅海
2、的海底地形。传统 SAR 图像水下地形探测基于波周期不变的假设,这不仅需要已知初始水深求解波周期,还给探测结果带来一定误差。本文提出一种基于子孔径图像的水下地形探测方法,将一景SAR图像分解成多景时间间隔固定的子孔径图像,利用子孔径图像间的时间间隔,求解不断变化的波周期,从而获得更加精确的水下地形。使用在海南蜈支洲岛的 TerraSAR-X 数据验证了此方法的可行性,将此方法反演得到的结果与 GEBCO 数据进行比较,发现两者吻合较好(MAE 为 2.8 m,MRE为23.91%),证明了此方法在浅海反演水深的巨大潜力。关键词:浅海地形探测;子孔径分解;线性频散关系中图分类号:TP75;P71
3、4+.7 文献标志码:A 文章编号:2095-1000(2024)01-0106-10DOI:10.12347/j.ycyk.20231115002引用格式:荣飞成,孟俊敏,纪永刚.基于SAR子孔径图像的浅海水下地形探测 J.遥测遥控,2024,45(1):106115.Shallow Sea Underwater Topography Detection Based on SAR Subaperture ImageRONG Feicheng1,2,MENG Junmin1,JI Yonggang2(1.First Institute of Oceanography,Ministry of N
4、atural Resources,Qingdao,266061,China;2.College of Oceanography and Space Informatics,China University of Petroleum,Qingdao 266580,China)Abstract:Accurate measurement of underwater terrain in coastal areas is the key to human activities in the ocean.Synthetic aperture radar(SAR)provides a new means
5、for the detection of shallow sea terrain.TerraSAR-X bunching mode can obtain high-resolution SAR data with a long integration time,so that the seabed terrain of shallow sea can be inverted more accurately.Traditional SAR image underwater terrain detection is based on the assumption of constant wave
6、period,which not only needs to know the initial water depth to solve the wave period,but also brings some errors to the detection results.In this paper,an underwater terrain detection method based on subaperture images is proposed.A SAR image is decomposed into subaperture images with fixed multi-sc
7、ene time intervals,and the time interval between subaperture images is used to solve the changing wave period,so as to obtain more accurate underwater terrain.The TerraSAR-X data of Wuzhizhou Island in Hainan Province was used to verify the feasibility of this method.The results obtained by inversio
8、n of this method were compared with the GEBCO data,and it was found that the two were in good agreement(the MAE was 2.8m and the MRE was 23.91%).It is proved that this method has great potential for retrieving water depth in shallow sea.Keywords:Underwater topography detection;Subaperture decomposit
9、ion;Linear dispersion relationCitation:RONG Feicheng,MENG Junmin,JI Yonggang.Shallow Sea Underwater Topography Detection Based on SAR Subaperture ImageJ.Journal of Telemetry,Tracking and Command,2024,45(1):106115.基金项目:国家自然科学基金(51839002)通信作者:孟俊敏()收稿日期:2023-11-15;修回日期:2023-12-042024 年 1 月遥 测 遥 控0引言浅海是
10、人类与海洋接触最频繁的区域,也是人类探索大洋的必经之路,因此,浅海地形的探索,对海洋研究、海上通行等具有重要作用。传统上使用声纳进行浅海地形的现场探测方法,探测范围小、花费高,对于船只难以进入的区域,无法开展现场探测。而遥感能够提供大量、大范围的图像数据,有效弥补了传统方法的不足,因此遥感已成为水下地形探测最有效的手段。遥感探测水下地形主要分为光学探测和 SAR遥感探测,其中光学遥感易受云雾等天气环境的影响,而SAR遥感可以全天时、全天候成像,在进行浅海地形反演时具有更大的优势。水下地形在SAR图像上表现为亮度特征的变化,传统基于SAR图像探测海底地形的方法是对水下地形与潮流、表面流与表面微尺
11、度波、表面微尺度波和雷达入射波之间的相互作用分别建模,进而求解得到水下地形。虽然此方法物理意义明显,但计算量大且需要大量实测数据。当海浪传播到浅海时,会受到海底地形的作用而导致波参数变化,此过程可用频散关系来表示。频散关系描述了浅水波波长、波周期与水深之间的关系。因此,可以从SAR海表图像中提取波 特 征 后,带 入 频 散 关 系 即 可 得 到 水 深13。Pereira 等人的研究证明了快速傅里叶变换(FFT)在提取 SAR 图像波参数的可靠性4。由于单景SAR 图像不包含时间的连续信息,因此需要根据外部已知水深数据的区域求解主波周期。而凭借波在传播过程中周期变化较小的特点,假设波在传播
12、过程中周期保持不变,即在探测区域中共用一个周期,就可以使用少量已知水深点反演大范围水深。Leu等人阐述了关于波周期不变的观点,认为海浪在浅海传播时,波浪波长与水深之间的线性频散关系保持守恒5。黄龙宇分别用 GF-3、Sentinel-1、ALOS PALSAR、ENVISAT ASAR 四种不同的SAR卫星图像对大洲岛附近水深进行了反演,得到的平均相对误差在30%以上,最大误差接近40%6。尽管使用不变波周期方法反演水深简单方便,但浅海中的波周期是随波浪的传播而不断变化的。卞小林从频散关系公式出发,探讨了波长、波浪周期和初始水深的敏感性,发现波周期误差对水深的影响最大,即小的波周期变化就能对反
13、演得到的水深产生较大的影响7。因此,精确获取不同位置的波周期是精确反演水深的关键。利用时间间隔较短的两景图像,得到波在传播过程中变化的波速,进而为水深的精确反演提供了可能。多通道光学图像可以满足此要求,现已在SPOT-58,9、Sentinel-21013等光学图像上证明了此方法的可行性。现有的利用SAR图像反演浅海水深的方法基本是以波周期固定不变的假设为前提,而SAR 是以天为单位的重访时间,阻碍了在 SAR 探测海底地形的方法中求解可变的波周期。为此,进一步从 SAR 成像机制出发,应用子孔径分解技术使一幅 SAR 图像分解得到两幅存在一定时间间隔的子视图像,从而将它们类比于光学数据中不同
14、通道的图像作为实验数据,但这对SAR图像的质量提出了一定要求。2014年和 2016年发射的 Sentinel-1A/B 卫星的产品是近年来最常用的一种 SAR 数据,但由于 Sentinel-1 IW 与 EW 模式使用TopSAR技术进行成像,牺牲了成像积分时间以得到多个条带,使得产品成像积分的时间较短,分解后的子 视 图 像 质 量 较 差 而 无 法 被 使 用14。其 前 身ENVISAT ASAR 卫星的IMAGE模式数据虽然成像时间大于Sentinel-1,但其分辨率低,难以从子视图像中观察到间隔时间内波浪的运动。因此,本研究采用 Terrasar-X聚束模式数据,利用其较长的积
15、分时间确保了子视图像的图像质量,高的分辨率使观测小波速波浪的运动成为了可能,低的运行轨道导致小的截断波长,为使用图像谱反演海浪信息创造了可能15。综上所述,利用求解可变波周期来代替原本基于 SAR 图像波浪的海底地形反演算法中波周期不变的假设,可以反演得到精确的海底地形。本文基于TerraSAR-X卫星遥感图像,提出了一种基于SAR子孔径图像波浪的海底地形反演方法,利用子孔径图像引入的时间信息,得到波浪周期,再结合图像提供的波长信息,求解反演研究区域的海底地形。1研究区域与数据介绍1.1研究区域基于波浪的水下地形反演方法适用于海表存在波浪运动且地形变化较为缓慢的浅海。因此本文选取海南岛东海岸的
16、蜈支洲岛附近水域作为研107第 45 卷第 1 期荣飞成等,基于SAR子孔径图像的浅海水下地形探测究区域。蜈支洲岛位于海南省三亚市北部海棠湾内,此区域中心纬度约 18.31,经度约 109.76,周围深度最大约50 m,经常存在由南海传播过来的波浪。图1所示为研究区域处的水深情况,图中黑色方块区域为SAR影像所处位置。1.2数据介绍1.2.1SAR数据子孔径分解会降低图像的方位向分辨率,因此需要高分辨率图像,本文使用TerraSAR-X聚束模式数据用于实验。TerraSAR-X 卫星于2007年6月成功发射,重访周期 11 天,在 X 波段下存在ScanSAR、StripMap、SpotLig
17、ht 和 HighResolution SpotLight、Staring-Spotlight五种不同的成像模式,本文使用的是TerraSAR-X SpotLight模式数据,其积分时间约1.2 s,距离分辨率约为2 m,方位分辨率约为1 m。研究区域的TerraSAR-X图像如图2所示,获取时间为2018年3月18日10427 UTC,极化方式为VV,中心入射角为29.5,从图2(右图)可以看到清晰的海浪条纹。1.2.2水深数据本 研 究 使 用 全 球 海 洋 总 测 深 图(General Bathymetric Chart of the Oceans,GEBCO)的数据作为参考水深数据
18、,对实验反演结果进行对比和分析。GEBCO是一个包括陆地和海洋的连续的全球地形模型,空间分辨率为0.25,其主要基于多波束测量数据,以网格化的形式提供了全球的地形16,可以在官方网站上免费下载(https:/)。1.2.3ERA5数据本研究使用 ERA5(ECMWF Reanalysis v5)数据来查找波浪参数。ERA5 数据是由 ECMWF(欧洲中期天气预报中心)提供的再分析数据集,空间分辨率为0.125,时间分辨率为1小时,查询研究区域内波浪的有效波高。2方法频散关系描述了浅水波波长、波周期及所处水深三者之间的关系,因此,通过从SAR图像中得到波长与波周期即可求得水深。传统方法假设波周期
19、在传播过程中不变,利用已知水深点结合SAR图像中得到的波长即可求解波周期,从而得到反演水深值,虽然简单,但忽视了波周期的变化,产生一定误差,因此本文通过对SAR图像进图1研究区域Fig.1Study area图2蜈支洲岛附近海域TerraSAR-X SAR图像,获取时间为2018年3月18日10 42 07 UTCFig.2TerraSAR-X SAR image of waters near Wuzhizhou Island,captured at 10 42 07 UTC on March 18,20181082024 年 1 月遥 测 遥 控行子孔径分解,利用子孔径图像间的时间信息求解波
20、周期,克服了传统方法中使用波周期不变假设带来误差的缺陷,以反演更为精确的海底地形。该方法首先在预处理后的SAR图像上确定反演区域,将反演区域中每个反演点对应的子图像进行子孔径分解,利用子孔径图像求解交叉谱,通过交叉谱得到此处反演点的波长与波周期,从而求解得到该处的水深;最后根据网格法求解该区域多个反演点水深,从而得到该区域的海底地形。2.1SAR图像预处理与反演区域确定SAR图像的预处理主要分为辐射定标、斜地距转换和重采样三个步骤。该方法只适用于有限水深区(水深在波长的1/20到1/2的区域),通过目视得到该区域波长约为60 m,再根据参考水深找到符合条件的反演区域,将此区域中陆地部分去除后,
21、得到反演区域。图3中的黄色区域即为反演区域。2.2波浪参数估计与处理2.2.1子孔径分解SAR图像通过多普勒处理将合成孔径时间内不同位置、不同方位角的回波信号进行数字合成,以实现高的方位向分辨率。此过程的逆过程便是子孔径分解,将方位向回波信号组成的多普勒频带进行等间隔划分,每个子多普勒频带对应时间内的回波信号积累成像,便可得到固定时间间隔的多景子孔径图像17。不同子孔径图像的孔径中心不同,使得每个子孔径图像到同一地物的方位视角不同,因此两个子孔径图像的后向散射信息存在显著差异。需要注意的是,子孔径分解会降低图像的分辨率18,其随着子孔径图像的数量成比例下降,并且子孔径频带范围不能超过全孔径数据
22、的多普勒带宽范围,否则会在左右两端出现卷绕现象。两幅子孔径图像的时间间隔可利用分割得到的多普勒频带宽和多普勒调频率FM得到:t=fFM(1)FM=-2V2ER0(2)其中,E为雷达波长、V为卫星运动速度、R0为卫星与目标之间的斜距。通过引入时间信息,可以求得变化波周期。假设两幅子孔径图像的时间间隔为Dt,同一位置处波相位变化为D,则该位置对应的波周期为:T=2DtD(3)2.2.2单点波浪参数求解海浪信息可以通过求解两幅子孔径图像的交叉谱得到。SAR图像的交叉谱能在保持与图像谱相似谱型特征的前提下,大大降低斑点噪声水平,同时也解决了波向 180方向模糊的问题19。将由SAR图像子视分解得到的两
23、幅子孔径图像进行二维FFT,即可获得二维波数谱。F(kxky)=-I(xy)e-i(kxx+kyy)dxdy(4)其中,kx和ky是子图像对应两个方向上的波数,它的范围从-Dx到Dx,Dx为这个方向上的像素间隔,波数分辨率为2Dxm,m为子图像在此方向上像素点的数量,x,y是子图像的中心坐标点。对两个子孔径图像的波数谱利用公式(5)即可得到此区域的交叉谱S(k)。S(k)=F1(k)F*2(k)(5)其中,*表示复数的共轭,主波波长、波向(波动传播方向,即去向)20和相位21可分别通过公式(6)(7)(8)求得。L=2k2xmax+k2ymax(6)=arctan()kymaxkxmax(7)
24、(k)=arctan()Im(S(kymaxkxmax)Re(S(kymaxkxmax)(8)其中,L、为子图像的波长和波向,为对应的相图3反演区域Fig.3Inversion region109第 45 卷第 1 期荣飞成等,基于SAR子孔径图像的浅海水下地形探测位差,Re和Im分别代表取复数的实部与虚部,每一个波数对应一个相位差。kymaxkxmax分别为波数空间中的峰值坐标。2.2.3多点波参数求解反演区域内所有点的波参数可以用网格法求解。以反演区域中一点为中心,确定子图像范围,求解中心点波参数后,将中心点沿x方向或y方向移动固定距离,直到遍历结束整个反演区域,即可得到反演区域中每个反演
25、点的波长与波周期。本文的移动距离设置为50个像素,初始子图像大小设置(250250)个像素,放大因子k=1+0.15 t是距海岸线像素个数t的函数,随着深度的增大而不断增大。图4为多点波参数反演过程,其中图4(a)诠释了网格法,红色框为选定的子图像区域,如图4(b)所示。将子图像分解为两幅子孔径图像,如图 4(c),利用公式(4)、公式(5)得到频域中交叉谱图4(d),再对其进行平滑滤波以降低噪声的影响。根据公式(6)和公式(8)从交叉谱中得到主波的波长与相位变化,再将子图像框在原图像上移动到图4(a)中蓝色框,直到遍历完全部反演点;再根据公式(1)和公式(2),计算得到两幅子孔径图像的时间间
26、隔约为0.48 s;最后利用公式(3)将求得的相位转换为周期,即可得到每个反演点波长与相位信息。2.2.4波浪参数的处理由于SAR图像容易受噪声及其他非海浪特征的干扰,使得反演结果存在一定的异常。对波长的影响因素主要包括 FFT 区域设置和滤波器参数设置。当SAR图像存在大量散斑噪声时,反演的波长会太大(尺度可达公里),所以需要设置一个波长上限,这里将波周期上限对应的波长设定为波长的上限,波长的下限由截断波长得到20,21。海浪在 SAR图像 中成像受非线性调制的速度聚束影响,会引起波在方位向产生位移。如果波在方位向位移小于一个波长,则波在图像上的映射仍为线性,此时速度调制对成像的影响较小。当
27、位移大于一个波长时,波在图像上的映射变为非线性,此时会使图像上波的成像发生扭曲,导致SAR图像模糊,即沿方位向传播的低于某个波长的海浪被截断,这个波长就称为截断波长。距离向的截断波长很小,只依赖于距离向分辨率,通常将其作为距离向截断波长22。而对于方位向传播的波来说,截断波长根据经验公式(9)得到。LMin|Azim=RVHs(9)其中,R 为倾斜距离,V 为平台速度,Hs为有效波高。对于方位向与波传播方向夹角为的波来说,对应的截断波长由公式(10)得到23LMin=LMin|Range2(1-cos(2)+LMin|Azim2(1-cos(2+)(10)利用公式(7)与子图像求得波向与方位向
28、之间图4求解多点波参数流程示意图Fig.4Flow diagram of solving multi-point wave parameters1102024 年 1 月遥 测 遥 控的夹角,通过查询ERA5数据可以得到此处区域的有效波高,代入公式(10)得到截断波长,此图像截断波长为7 m。波 周 期 上 限 设 置 为 18 s,对 应 的 波 速 为28.07 m/s,对于浅水波,不太可能以此速度在沿岸区域进行传播,因此这是一个较为宽松的上限10,下限设置为截断波长对应的波周期。由于小于截断波长的波都认为是错误的,因此,将此波长范围内对应的波周期2 s认为是波周期的下限。将异常点设为非数
29、值后,对波参数矩阵进行二维平滑滤波24,实现保持浅海地形变化趋势的同时,又能减缓相邻反演点处受噪声干扰而变化剧烈。最后对波周期和波长进行频散关系滤波。线性频散关系描述了浅水波的波参数与水深之间的联系。波浪从浅海传播时,随着水深的减小,海表面波逐渐受到海底地形的作用而发生形变,导致波长减小,波速降低,根据能量守恒定理,波高会随之增加,从而使得波浪变陡,当波高达到一定高度时,波浪就会发生破碎现象,使海表面粗糙度突然增大25。在这个过程中,波浪的波长和周期与波传播过程中的水深关系被称为频散关系,由公式(11)来描述。w=gktanh(kh)(11)其中,w为波的角频率()w=2T,g=9.8 m/s
30、2为重力加速度,h为水深。当波浪在深水区传播时,波浪的传播不再受到海底地形的影响,此时波浪的波参数变化与水深无关,满足公式(12):w=gk(12)其中,由于频散关系中的反双曲正切函数是一个到达边界不稳定的函数,在边界处剧烈变化,这会造成一定的误差。因此本研究只选取其自变量小于0.98的部分26,去除不符合频散关系反演点的波长与波周期后,利用剩余反演点的波长与波周期进行水深反演。图 5展示了频散关系滤波的过程,图中每一个点代表了此反演区域中的一个反演点,横轴表示该反演点的波长、竖轴表示该反演点的波周期,颜色代表此波长、波周期对应波数所有具备的强度。从图像中可以看到,反演区域中所有反演点的波周期
31、都大于截断波长,因此所有反演得到的波长都是有效的。左边图像为滤波前,右边图像为滤波后,通过频散滤波去除了一些反演得到的波参数不符合频散关系的无效反演点。2.3水深反演与海底地形估计求解得到每一个有效反演点对应的波参数后,便可利用频散关系得到每一个反演点对应的水深。将公式(11)进行变化可得水深求解公式如下:h=1karctanhw2gk(13)由于去除了一些不符合实际和频散关系的反演点,使得到的反演区域中存在一定的缺漏,因此需要进行插值,从而得到该区域下完整的海底地形。3结果与讨论3.1地形反演结果利用 SAR 图像反演得到的水深结果如图 6所示。从图中可以看到,反演结果基本遵循沿岸越远、水深
32、越深的变化规律,且反演得到的水深主图5频散关系滤波Fig.5Dispersion relation filtering111第 45 卷第 1 期荣飞成等,基于SAR子孔径图像的浅海水下地形探测要分布在20 m以下,在2030 m处水深反演点较少。将反演得到剩余水深进行插值即可得到图7中的地形反演图。3.2讨论3.2.1与GEBCO参考水深比较根据GEBCO数据,查询每一个反演点对应的参考水深。由于参考水深分辨率较低,使用克里金插值将经纬度对应的参考水深数据插值到与反演点经纬度相同的分辨率,利用经纬度查询每一个反演点的参考水深,得到参考水深图如图 8所示。根据反演区域的经纬度和当地时间,通过查
33、询 TMD(Tidal Modal Driver)数据27,可以得到该位置该时刻的潮高,得到该处的潮汐分量为-0.1 m,对反演结果进行潮汐校正后与参考水深相比 较,得 到 MRE 为 23.91%、MAE 为 2.80 m,RMSE为3.30 m和反演率为73.32%。通过与图6反演结果对比发现,反演结果基本符合参考水深的分布,但水深较深处反演点缺少较为严重,这主要是频散滤波造成的。通过图5的频散滤波可以看到,水深越大,即越靠近深水区,波长和周期对水深的影响越大,它们轻微的改变就会造成水深较大的偏差。通过频散滤波虽然能去除这些不符合频散关系的较大偏差点,但也造成了大量水深较深处反演点的丢失。
34、同时,图8中黄色区域显示出此处地形变化较大,而此方法求得的波长是连续变化的,因此也会反演得到一些不符合频散关系的异常点。当波从深海传播到浅海时,由于岛屿的存在,在经过岛屿时会进行绕射,靠近岛屿的左侧图像没有波浪,从而造成图8中白色区域反演点的部分缺失且反演效果较差。图9显示为反演区域中各个反演点与参考水深之间的平均相对误差,将误差分为小于 10%、10%20%、20%30%、30%40%和 40%以上五个部分,通过图9可以看到,当反演点所处浅水区时,平均相对误差较大,这是由于此时波浪受到海底图6反演水深Fig.6Inversion bathymetry图7反演水下地形Fig.7Inversio
35、n underwater topography图8反演区域对应的参考水深 Fig.8The reference water depth corresponding to the inversion area图9反演区域的误差分布Fig.9Error distribution in the inversion region1122024 年 1 月遥 测 遥 控地形的约束力越来越大,导致波长、波速降低、波陡增大,波浪越发不稳定,从而求解得到的周期分布较大,产生了较大的平均相对误差。3.2.2与传统方法的比较使用波周期不变的方法对相同反演区域中的反演点进行反演。首先需要在深海区选取四个区域用于求解
36、波周期,如图10所示,每一个区域大小设置为512像素512像素。将截取的图像进行二维 FFT转换到波数域求解波长,结合查询的水深信息,即可得到对应的波周期。图11分别展示了子图像区域时空域的图像和对应的二维频域图像,各自结合水深得到的周期见表1。从表1中可以看到,四景子图像所对应的波浪全部为深水波,此时波浪波参数变化与水深无关,通过公式(12)利用SAR图像中反演的波长即可得到波周期。四景不同的子图像求解得到的波周期是不同的,分别使用它们作为反演区域中固定的波周期,反演相同反演区域下的水深,反演情况见表2。从表2中发现,将来自不同子图像反演得到的波周期分别应用到水深反演区域得到的水深与参考水深
37、相比具有相似的平均相对误差,但在平均绝对误差、均方根误差和反演率上存在一定的差异。取四个波周期的平均值作为反演区域的波周期,进行水深反演,得到最好的反演效果,但仍差于本文方法取得的反演结果,这主要是由于两方面造成的:一是选取子图像估计波周期时,四幅子图像反演得到的波周期存在一定波动,选取子图像的质量会对波周期的反演产生较大的影响,从而对反演结果造成一定的误差;二是由于反演区域内只使用一个波周期,忽略了波周期较小的变化,从而给反演结果带来了一定的误差,这在小波长、周期的区域下尤为明显,这是因为此区域下频散关系更为严格。图10波周期反演区域在子图像中的位置 Fig.10Position of th
38、e wave period inversion region in the subimage图11子图像及对应的二维频谱图Fig.11Subimage and corresponding two-dimensional spectrogram表1不同子图像求解出来的周期Table 1The period of solving different sub-images图像1234波长/m61.5361.5371.9565.29参考水深/m35404750波周期/s6.2856.2826.7946.470113第 45 卷第 1 期荣飞成等,基于SAR子孔径图像的浅海水下地形探测4结束语本文提出了
39、一种基于子孔径分解的可变波周期的浅海地形探测方法。子孔径分解技术可以获得两幅时间间隔合适的子孔径图像,从而求解得到交叉谱,利用交叉谱的峰值信息确定主波波长与相位变化,带入频散关系即可得到每一个反演点处的水深信息。此方法不需要初始水深数据,并且求得可变的波周期来代替原来波周期不变的假设,从而减小了周期与初始水深带来的误差,得到较高的反演精度。使用TerraSAR-X 聚束模式数据在海南蜈支洲岛附近水域进行实验,并将反演结果与GEBCO参考水深数据进行对比,发现反演效果(MAE为 2.8 m,MRE为 23.91%)优于传统方法(MAE为 4.97 m,MRE为 35.91%)。由于此方法只适用于
40、有限深水区,同时此影像下该处波浪波长位于5070 m,整体较小,该波长对应的频散关系较为敏感,较小的波长变化就会造成水深较大的改变,从而使得在海岛周边浅水区域和有限水深与深水区的交界处,受限于频散关系,反演率较低、误差较高。在未来,可以尝试使用更多较大波长的SAR数据来共同反演相同位置的水深,补充丢失反演点的水深,从而得到更高空间分辨率与精度的浅海地形。参考文献1 BIAN X,SHAO Y,WANG S,et al.Shallow water depth retrieval from multitemporal sentinel-1 SAR dataJ.IEEE Journal of Sel
41、ected Topics in Applied Earth Observations and Remote Sensing,2018,11(9):2991-3000.2 SANTOS D,ABREU T,SILVA P A,et al.Estimation of coastal bathymetry using waveletsJ.Journal of Marine Science and Engineering,2020,8(10):772.3 CLOAREC M,ROEBER V,RANCHIN T,et al.Extraction of bathymetric features usin
42、g multiple SAR images produced by Sentinel-1C/Sixth International Confe-rence on Remote Sensing and Geoinformation of the Environment(RSCy2018).SPIE,2018,10773:274-282.4 PEREIRA P,BAPTISTA P,CUNHA T,et al.Estimation of the nearshore bathymetry from high temporal resolution Sentinel-1A C-band SAR dat
43、a-A case studyJ.Remote Sensing of Environment,2019,223:166-178.5 LEU L G,KUO Y Y,LIU C T.Coastal bathymetry from the wave spectrum of SPOT imagesJ.Coastal Enginee-ring Journal,1999,41(1):21-41.6 HUANG L,MENG J,FAN C,et al.Shallow sea topography detection from multi-source SAR satellites:A case study
44、 of dazhou island in chinaJ.Remote Sensing,2022,14(20):5184.7 BIAN X L,SHAO Y,ZHANG C,et al.The feasibility of assessing swell-based bathymetry using SAR imagery from orbiting satellitesJ.ISPRS Journal of Photogrammetry and Remote sensing,2020,168:124-130.8 HOLMAN R,MOREIRA A.Real-time synthetic ape
45、rture radar(SAR)processing with a new subaperture approachJ.IEEE Transactions on Geoscience and Remote sensing,1992,30(4):714-722.9 DE MICHELE M,LEPRINCE S,THIBOT J,et al.Direct measurement of ocean waves velocity field from a single SPOT-5 datasetJ.Remote Sensing of Environment,2012,119:266-271.10
46、BERGSMA E W J,ALMAR R,MAISONGRANDE P.Radon-augmented sentinel-2 satellite imagery to derive wave-patterns and regional bathymetryJ.Remote Sen-sing,2019,11(16):1918.11 ALMAR R,BERGSMA E W J,THOUMYRE G,et al.Global satellite-based coastal bathymetry from wavesJ.Remote Sensing,2021,13(22):4628.12 BABA
47、M W,THOUMYRE G,BERGSMA E W J,et al.Deriving large-scale coastal bathymetry from Sentinel-2 images using an high-performance cluster:A case study covering North Africas coastal zoneJ.Sensors,2021,21(21):7006.13 DALY C,BABA W,BERGSMA E,et al.The new era of regional coastal bathymetry from space:A show
48、case for west africa using optical sentinel-2 imageryJ.Remote Sensing of Environment,2022,278:113084.14 DE ZAN F,GUARNIERI A M.TOPSAR.Terrain observation by progressive scansJ.IEEE Transactions on Geoscience and Remote Sensing,2006,44(9):2352-2360.15 LI X,LEHNER S,ROSENTHAL W.Investigation 表2使用不同子图像
49、求解的周期反演效果 Table 2The effect of periodic inversion using different subimages图像1234平均MRE/%37.336.9538.9936.2035.91MAE/m6.466.214.975.324.97RMSE/m8.317.946.317.16.24反演率/%56.4056.3771.5367.4466.451142024 年 1 月遥 测 遥 控of ocean surface wave refraction using TerraSAR-X dataJ.IEEE Transactions on Geoscience
50、and Remote Sensing,2009,48(2):16 WEATHERALL P,TOZER B,ARNDT J E,et al.The GEBCO_2021 Grid-a continuous terrain model of the global oceans and landJ.2021.17 FERRO-FAMIL L,REIGBER A,POTTIER E,et al.Scene characterization using subaperture polarimetric SAR dataJ.IEEE Transactions on Geoscience and Remo