海洋白帽的高空间分辨率光学遥感估算分析

赵碧; 丁静; 刘建强; 焦俊男; 唐君; 陆应诚

doi:10.11834/jrs.20222106

水环境与资源遥感监测应用 | 浏览量 : 1429 下载量: 678 CSCD: 0 更多指标

R-PDF
PDF
导出
分享给微信好友或者朋友圈

分享
收藏
专辑

海洋白帽的高空间分辨率光学遥感估算分析
Estimation of oceanic whitecaps using high spatial-resolution optical remote sensing
- 角色： First author 第一作者
- 机构：
  南京大学国际地球系统科学研究所, 南京 210023
- 邮箱： Zhaobi1929@163.com
- 简介：赵碧，研究方向为海洋环境遥感。E-mail： Zhaobi1929@163.com
赵碧
1 ，
- 机构：
  国家卫星海洋应用中心, 北京 100081
  自然资源部空间海洋遥感与应用重点实验室, 北京 100081
丁静
2 3 ，
- 机构：
  国家卫星海洋应用中心, 北京 100081
  自然资源部空间海洋遥感与应用重点实验室, 北京 100081
刘建强
2 3 ，
- 机构：
  南京大学国际地球系统科学研究所, 南京 210023
焦俊男
1 ，
- 机构：
  南京大学国际地球系统科学研究所, 南京 210023
唐君
1 ，
- 角色： Corresponding author 通讯作者
- 机构：
  南京大学国际地球系统科学研究所, 南京 210023
- 邮箱： Luyc@nju.edu.cn
- 简介：陆应诚，研究方向为海洋遥感研究。E-mail： Luyc@nju.edu.cn
陆应诚
1
2023年27卷第1期页码：92-103
纸质出版日期： 2023-01-07 ，
DOI： 10.11834/jrs.20222106

扫描看全文

引用

阅读全文PDF
赵碧，丁静，刘建强，焦俊男，唐君，陆应诚.2023.海洋白帽的高空间分辨率光学遥感估算分析.遥感学报，27（1）： 92-103

Zhao B，Ding J，Liu J Q，Jiao J N，Tang J and Lu Y C. 2023. Estimation of oceanic whitecaps using high spatial-resolution optical remote sensing. National Remote Sensing Bulletin， 27（1）：92-103
赵碧，丁静，刘建强，焦俊男，唐君，陆应诚.2023.海洋白帽的高空间分辨率光学遥感估算分析.遥感学报，27（1）： 92-103 DOI： 10.11834/jrs.20222106.

Zhao B，Ding J，Liu J Q，Jiao J N，Tang J and Lu Y C. 2023. Estimation of oceanic whitecaps using high spatial-resolution optical remote sensing. National Remote Sensing Bulletin， 27（1）：92-103 DOI： 10.11834/jrs.20222106.

海洋白帽的高空间分辨率光学遥感估算分析

赵碧 ; 丁静 ; 刘建强 ; 焦俊男 ; 唐君 ; 陆应诚 ;

PDF下载

摘要

海洋白帽是海面波浪破碎的直观表现，反映了海面风速等海洋动力环境情况，是海气交换的重要通道，也是海洋研究关注的重点目标之一。海洋白帽对入射光具有较强的反射散射作用，随着海洋光学遥感空间分辨率的提高，海洋白帽能被有效识别与区分。在空间分辨率为10 m的哨兵2号卫星（Sentinel-2）多光谱成像仪MSI（Multi-Spectral Instrument）数据的近红外波段（842 nm）中，不同风速下的海洋白帽影像特征清晰；基于海洋白帽对入射光的辐射传输过程与图像特征，本研究利用一种白帽覆盖率的光学遥感估算方法，实现了MSI弱耀光反射图像上的白帽覆盖率估算。研究表明：（1）基于MSI数据估算的海面白帽覆盖率，在无其他海洋环境动力的影响下，能进一步用于海面风速的遥感反演，且具有较高的空间分辨率；（2）其他海洋动力环境要素，能对海面波浪的破碎程度产生调制作用，因此海洋白帽的识别提取也有助于海洋环境动力遥感监测；（3）海面强耀光反射会对海洋白帽遥感带来较大的干扰，导致海面白帽覆盖率的高估，有效剔除耀光反射是提高海洋白帽识别提取精度的重要前提。对中国首个海洋水色业务卫星星座（海洋一号C/D卫星，HY-1C/D）获取的海岸带成像仪CZI（Coastal Zone Imager）数据开展分析，结果表明，在无其他海洋环境动力的调制下，当风速约高于9 m/s时，海洋白帽能为CZI载荷探测并区分，白帽覆盖率随着风速的增加而增加；CZI瑞利校正反射率图像的统计分析还表明，白帽反射率会比背景海水反射率高约5.8%—8.3%；对CZI图像中海洋白帽的深入研究，不仅为大气校正提供数据参考，也有助于光学遥感对海洋环境的深入理解与应用开发。随着海面白帽光学遥感的不断深入，今后将能为海面风场反演、海洋环境动力监测、海气交换评估提供新的技术方法参考。

译

Abstract

Oceanic whitecaps, generated from wind-wave breaking process, are the medium of air-sea exchange and the indicator of sea surface state. Due to the intense reflection and scattering of incident light, whitecaps can be effectively recognized and distinguished in in-situ photos or videos. Whitecap coverage (W), defined as the proportion of ocean surface covered by whitecaps, is an important parameter for the quantification of whitecaps. Of course, oceanic whitecaps can also be discriminated in high spatial resolution optical remote sensing images, such as Sentinel-2 MSI and Landsat-8 OLI data. This can provide a new research direction in marine environment observation, and may be further used for wind speed monitoring. However, how to estimate oceanic whitecap coverage from these optical remote sensing images is still a challenge. In this study, the formula of whitecap coverage is obtained by converting the form of radiative equation related with the constant and image reflectance of Sentinel-2 MSI data under some assumptions. The background signal of seawater and atmosphere was eliminated by data optimization, and the signal of whitecaps can be distinguished using regional filtering method. Then, whitecap coverage can be estimated. The identification and estimation results indicate that whitecap coverage derived from Sentinel-2 MSI images are consistent with previous studies using in situ observations in the order of magnitude, and can invert sea surface wind speed using a statistical model. Coarse spatial resolution wind speed images converted from MSI inversion were validated by ERA5 wind speed products from European Centre for Medium-Range Weather Forecasts (ECMWF). In addition, whitecap coverage can imply the modulation of other marine environmental dynamic factors (e.g., water mass, ocean fronts, ocean eddies and internal waves). Moreover, it should be noted that sunglint reflection is a non-negligible issue for optical remote sensing of oceanic whitecaps whose signal should be effectively eliminated. Examining the MSI-derived results to HY1-C/D CZI, it indicates that whitecaps can be identified in CZI images when wind speed is greater than 9 m/s, and the reflectance difference between whitecaps and background seawater is 5.8%—8.3%. We hope the above results can be used to improve the accuracy of atmospheric correction, and provide new reference for using high spatial resolution optical remote sensing in sea surface wind speed estimation and marine environmental dynamic factors monitoring. This will hopefully expand the research and application fields of ocean color remote sensing.

译

关键词

光学遥感; 海洋白帽; 耀光; HY-1C/D; Sentinel-2; MSI; CZI

译

Keywords

optical remote sensing; oceanic whitecaps; sunglint; HY-1C/D; Sentinel-2; MSI; CZI

译

1　引言

海面在风的作用下产生海浪，随着海表风速增加到某一临界值，海浪破碎并在波峰处产生大量的气泡，从而形成目视特征显著的白色气泡水体，被称为海洋白帽或白冠（Oceanic whitecaps，本文中统称为白帽）（

Raizer等，1994；Melville，1996；Pivaev等，2021）。海洋白帽是海面波浪破碎的直观表现（Callaghan等，2008a；Salisbury等，2014；Paget等，2015），是海气交换的重要媒介和通道（Frouin等，1996和2001；Callaghan等，2008a），在海气通量交换（Asher和Wanninkhof，1998；Goddijn-Murphy等，2016；Monahan等，1983）、海洋能量传递（Phillips，1985；Melville，1996）、海浪发展及波形破碎（Hanson和Phillips，1999；Melville，1996；Myrhaug和Holmedal，2008）等方面具有指示作用，是海洋环境监测研究的重要对象。

译

描述海洋白帽的属性参数主要有白帽反射率、白帽覆盖率等，海洋白帽在光学到微波谱段内的反射特性基本得到阐明（

Koepke，1984；Kokhanovsky，2004；Moore等，2000；Nicolas等，2001），其在可见光、近红外波段的反射率高于背景水体，因此海洋白帽在现场观测资料中易于目视识别和提取（Koepke，1984；Frouin等，1996和2001）。海洋白帽的现场观测分析基于原位观测的照片或视频资料，通过统计一定范围内白帽面积和发生频率等特征，从而揭示该区域环境动力的现象与影响因素（Monahan等，1983；Goddijn-Murphy等，2011；Callaghan等，2008a；Jia和Zhao，2019；Zhang等，2015）。基于现场观测资料能进一步给出白帽覆盖率，即白帽面积占总观测面积的比例（Whitecap coverage或Whitecap fraction）从而可建立起白帽覆盖率与风速之间的统计关系模型（Monahan和Muircheartaigh，1980；Callaghan等，2008b；Lafon等，2004）。不同海域的“白帽覆盖率—风速”模型研究，进一步揭示了风对海面（波浪破碎、海面粗糙度、海气交换等）的作用机制、影响程度及区域分异特点（Monahan等，1983；Goddijn-Murphy等，2011；Callaghan等，2008a；Jia和Zhao，2019；Zhang等，2015）。

译

卫星光学遥感技术开展海洋白帽的观测分析，是厘清不同区域海洋白帽的分布变化、揭示海洋环境动力特征的有效技术支撑（

Koepke，1984；Kokhanovsky，2004；Moore等，2000；Nicolas等，2001）。在海洋白帽遥感研究中，除了云雾的影响，图像空间分辨率、海洋水色差异、海面耀光反射、海洋环境动力等均会产生影响（Monahan 等，1983；Goddijn-Murphy等，2011；Callaghan等，2008a；Myrhaug和Holmedal，2008）。如在粗空间分辨率的MODIS数据中，难以直接从图像中鉴别海洋白帽，目前多利用统计模型来估算白帽反射信号贡献，提高MODIS数据的大气校正精度（Gordon和Wang，1994；Hu等，2000），因此，要从海洋遥感的高空间分辨率数据中（海洋遥感与陆地遥感在空间分辨率的需求与认知上存在差异，不展开讨论）进行白帽的识别提取、厘清海洋白帽的光学遥感特征。海洋白帽能在美国陆地卫星（Landsat 8）30 m空间分辨率的陆地成像仪OLI（Operational Land Imager）图像中清晰识别（Kubryakov等，2021），也能在欧空局哨兵2号卫星（Sentinel-2 A/B）10 m空间分辨率的多光谱成像仪MSI（Multi Spectral Instrument）图像中清晰展现（如图1所示）。

译

（a）海面照片中的白帽

（a） Whitecaps in sea surface photos

下载: 原图 | 高精图 | 低精图

（b） 8.2 m/s风速时的MSI 842 nm白帽影像

（b） Whitecaps in MSI 842 nm image at wind speed of 8.2 m/s

下载: 原图 | 高精图 | 低精图

（c） 11.2 m/s风速时的MSI 842 nm白帽影像

（c） Whitecaps in MSI 842 nm image at wind speed of 11.2 m/s

下载: 原图 | 高精图 | 低精图

图1 海洋白帽的照片与遥感影像

Fig. 1 Oceanic whitecaps in photograph and remotely sensed images

开展高空间分辨率海洋光学遥感数据中海洋白帽的识别提取与定量反演研究，挖掘光学遥感数据在海洋环境动力（如风场、流场等）监测中的应用潜力，能拓展海洋光学遥感的研究范畴。本研究基于Sentinel-2 MSI高分辨率数据，进行海面白帽覆盖率的光学遥感反演，探索白帽覆盖率对海面环境动力背景的响应，分析白帽覆盖率反演算法的适用条件、不确定性和干扰因素；面向中国首个海洋水色业务卫星星座（海洋一号C/D卫星，HY-1C/D）的高空间分辨率光学载荷——海岸带成像仪CZI（Coastal Zone Imager），开展不同风速背景下，CZI图像中白帽的可识别性研究，分析白帽反射信号在CZI各波段中的信号贡献差异，讨论海洋白帽在CZI数据中的存在特点；从而为国产自主海洋光学卫星的精细应用、海洋环境动力与海气交互监测提供参考。

译

2　数据与方法

2.1　数据及预处理

选用Sentinel-2卫星的MSI数据，其10 m空间分辨率的影像上能够清晰识别出海洋白帽。相对应的背景风速数据，选择欧洲中期天气预报中心ECMWF（European Centre for Medium-Range Weather Forecasts）的ERA5再分析风速产品，其空间分辨率为0.25°、时间分辨率为1 h。数据选择及预处理需基于如下几点考虑：（1）先根据背景风速优选MSI数据。白帽的产生需要一定的风速条件，据现场统计，在6 m/s以上的风速会产生海洋白帽，风速大于7—8 m/s时白帽迅速增多（

Monahan和Zietlow，1969），为了保证在MSI影像中检测到白帽，本研究选择8 m/s以上风速条件的影像。（2）对于高空间分辨率数据，太阳耀光反射的影响较粗空间分辨率数据更为复杂（Wu等，2022）。为了尽量避免耀光反射的影响，选择MSI数据的θ_m大于40°的影像。θ_m为海面太阳镜面反射角与传感器观测角度之间的夹角，由传感器与太阳的几何角度决定，是海面耀光反射强弱的直接指针，其大小可表征平均耀光反射的强弱，θ_m越大表示耀光反射越弱（Hu等，2009；Lu等，2016；Wen等，2018），θ_m的计算公式如下：

译

$c o s (θ_{m}) = c o s (θ_{0}) c o s (θ) - s i n (θ_{0}) s i n (θ) c o s (φ)$

(1)

式中，θ₀、θ和φ分别为太阳天顶角、传感器天顶角和相对方位角。（3）数据预处理。利用ACOLITE（版本20190326）软件对MSI L1C级别数据进行瑞利校正获得10 m分辨率的数据（R_rc）。鉴于云会明显干扰白帽覆盖率的估算，因此选择无云或少云的影像，对图像中少量的云进行掩膜，并对掩膜结果进行5个像元的膨胀处理。

译

对中国东黄海2021年MSI数据进行筛选，研究主要使用MSI 842 nm波段生成R_rc反射率产品，可忽略背景水体差异，具体区域和影像时间如图2所示。图2（a）、2（b）、2（c）分别由2景、5景、2景MSI影像拼接而成。由于影像数据的成像时间在2：16—2：41 UTC，因此将风速产品2：00 UTC及3：00 UTC的数据取平均。

译

下载: | 高精图 | 低精图

图2 Sentinel-2卫星MSI近红外波段R_rc（842 nm）影像图

Fig. 2 Sentinel-2 MSI near infrared R_rc（842 nm） images

2.2　海面白帽覆盖率估算

光学卫星数据中，大气顶层总辐亮度L_t可表示为

译

$L_{t} = L_{r} + L_{a} + L_{r a} + T L_{g} + t L_{W C} + t L_{w}$

(2)

式中，L_r、L_a和L_ra分别表示瑞利散射、气溶胶散射、瑞利—气溶胶散射辐亮度，L_WC示白帽反射辐亮度，L_w表示离水辐亮度，L_g表示耀光反射辐亮度，T、t分别为大气直射、散射透过率（

Gordon和Wang，1994）。根据公式：

译

$R = π L / (F_{0} c o s θ_{0})$

(3)

可将辐亮度转换为反射率，式（3）中R为反射率，L为辐亮度，F₀为太阳辐照度，θ₀为太阳高度角。海面反射率又分为受白帽影响及不受影响部分（Koepk，1984），结合式（2），海面反射率可表达为

译

$R_{t} = R_{r} + R_{a} + R_{r a} + T R_{g} + t W R_{W C} + t (1 - W) R_{w}$

(4)

式中，R_r、R_a和R_ra分别表示瑞利散射、气溶胶散射、瑞利—气溶胶散射反射率，R_wc示白帽反射率，R_w表示离水反射率，R_g表示耀光反射率，T、t分别为大气直射、散射透过率，W为白帽覆盖率（Whitecap coverage），表示海面被白帽覆盖的部分占总面积的比例，由式（4）可得

译

$W = (R_{t} - R_{r} - R_{b} - T R_{g}) / (t R_{W C} - t R_{w})$

(5)

式中， $R_{b} = R_{a} + R_{r a} + t R_{w}$ ，代表了来自背景海水与大气的反射率贡献。由于近红外波段（842 nm）中离水辐亮度可以忽略（即 $t R_{w} \approx 0$ ），且忽略耀光反射的影响（图2中数据的θ_m均大于40°，则 $T R_{g} = 0$ ）。因此，针对图2中白帽覆盖率可根据下式（6）进行估算：

$W = (R_{r c} - R_{b}) / (t R_{W C})$

(6)

近红外波段白帽的反射率（R_wc）约为0.55（

Koepke，1984；Frouin 等，1996和2001），海洋白帽的光学特征相对稳定，具有较为明确的光谱反射率特征（Kokhanovsky，2004；Garaba和Dierssen，2018），而图像上白帽反射率变化则是与图像中白帽覆盖率密切相关。近红外波段大气散射透过率（t）为0.7—0.8（Kneizys 等，1981），上述值仅可用于背景海水与大气均匀稳定且不受海面耀光反射干扰的区域。此外，利用式（6）进行白帽覆盖率（W）反演的关键在于准确的背景反射信号R_b的获取。本研究中，假定在4 km×4 km的区域内，海水和大气的光学属性是均匀的，对图像进行400×400（4 km×4 km）的窗口滤波，获得窗口内的最小值作为背景反射率，可进一步计算得出图像的白帽覆盖率（W）。特别需要说明的是，上述公式和参数有着严格的适用条件，本文原理分析中的数据边界清楚，可以适用，但是在未来的实际业务应用中，还需要发展业务化应用算法。

译

3　基于MSI影像的白帽覆盖率估算与分析

3.1　白帽覆盖率

基于MSI 842 nm波段R_rc数据（图2）反演的白帽覆盖率如图3所示，其中云、陆地、辐射沙脊等已经被掩膜处理。白帽覆盖率作为海面波浪破碎状况的直观体现，一方面反映了风速这一主要动力因素，与风速大小呈正相关；另一方面，又受到海洋涡旋、内波、锋面等的调制（

Monahan等，1983；Goddijn-Murphy等，2011；Callaghan等，2008a）。如在图3（a）中，2021年1月13日的白帽覆盖率的变化主要受到海面风速的影响，而在图3（b）和3（c）中，在2021年1月27日和2021年2月21日的白帽覆盖率在局部区域又受到海洋锋面、涡旋等的影响（具体分析见3.3节）。总体而言，图3中平均白帽覆盖率在1%上下（对应风速8—12 m/s），这与白帽的现场观测（基于海面实拍照片或视频等资料）分析结果较为接近（Blanchard，1983）。

译

下载: | 高精图 | 低精图

图3 基于MSI数据的白帽覆盖率估算结果

Fig. 3 MSI images derived whitecap coverage

为了验证白帽覆盖率算法的估算精度，本文针对图2（a）中红框框选出的12个区域开展精度评价。12个500×400像元验证区域的842 nm波段R_rc图像见图4第一和第三列，对这12个区域图像采用人工目视解译的方法，通过选取合适阈值分割出图像中的海洋白帽，并计算出该区域中的白帽覆盖率（图4第二、第四列）。将12个验证区域的算法反演结果与目视解译结果进行对比分析（如图5所示），算法反演的白帽覆盖率与目视解译的白帽覆盖率具有较好的一致性（R²=0.85），说明提取结果的精度可靠验证拟合线低于1∶1线，说明本研究算法在弱耀光条件下存在一定的低估。

译

下载: | 高精图 | 低精图

图4 MSI图像白帽目视解译结果（（a）—（l）为图2（a）从上到下红色方框部分）

Fig. 4 Visual interpretation of oceanic whitecaps from MSI images （（a）—（l） figures correspond to the red boxes of Fig. 2（a） from top to bottom）

图5 目视解译与算法估算白帽覆盖率对比图

Fig. 5 Comparison of visual interpretation and algorithmic estimation of whitecap coverage from MSI images

下载: 原图 | 高精图 | 低精图

3.2　基于白帽覆盖率的海面风速反演

（a）风速10.4 m/s，白帽覆盖率0.65% （b）风速10.4 m/s，白帽覆盖率0.67%

（a） U₁₀ = 10.4 m/s with W = 0.65% （b） U₁₀ = 10.4 m/s with W = 0.67%

下载: 原图 | 高精图 | 低精图

海面白帽是风生波浪破碎的直观表现，在其他海洋环境动力影响较为均匀的区域，白帽覆盖率与海面风速（海面10 m高处的风速，U₁₀）具有密切的统计关系（

Monahan和Muircheartaigh，1980；Anguelova和Webster，2006；Anguelova和Bettenhausen，2019）。基于海面实测数据，构建“白帽覆盖率—风速”统计模型，有利于对区域海洋环境动力等的深入研究。需要指出的是，本研究的重点在于能否反演出海面风速，并与已有风速产品进行对比分析，但基于本区域海面实测结果的“白帽覆盖率—风速”统计模型缺乏，因此选择Monahan和Muircheartaigh（1980）的白帽覆盖率—风速模型来反演海面风速。模型公式如下：

译

$W = 2.95 \times 10^{- 6} {U_{10}}^{3.52}$

(7)

式中，W代表白帽覆盖率，U₁₀为风速，对其进行变换可以给出风速反演模型如下：

译

$U_{10} = {(W / (2.95 \times 10^{- 6}))}^{0.284}$

(8)

利用式（8）结合图3（a）的海面白帽覆盖率提取结果，可以得到该区域风速反演结果（如图6（a），空间分辨率为4 km），图6（b）为该区域同时段的ERA5海面再分析风速产品（时间差<15 min，空间分辨率为0.25°）。对比图6（a）和6（b），整体风速为东北方向风速高，向西南方向风速递减，MSI推导的风速结果因为江苏沿海辐射沙脊的影响会出现一些反演异常值（见图6（a）左下角）。总体上，基于MSI白帽覆盖率反演的海面风速与ERA5海面风速产品的空间变化趋势与风速变化区间基本一致。给出二者的频率累计曲线对比如图6（c），用于进一步评估风速反演的精度。海面风速的频率累计曲线图可以说明，MSI白帽覆盖率反演的海面风速与ERA5海面风速产品变化基本一致，但MSI反演的海面风速略低于ERA5风速产品的值，二者之间的差异在风速为10 m/s左右达到最大，为1.6 m/s，大部分差异在1 m/s以内，这主要是由于算法对MSI图像中白帽覆盖率的低估造成（图5）。

译

3.3　其他环境动力对白帽覆盖率的调制

海洋白帽主要反映的是风生波浪的破碎过程，虽然海面风场是白帽产生的首要因素，但风生波浪也会被涡旋、内波等其他海洋环境动力所调制，从而改变白帽的发生概率（

Thorpe等，1987；Romero，2019）。图7分别反映了海洋锋面、涡旋、内波和船只尾迹对白帽覆盖率的影响：图7左列为MSI R_rc（842 nm）影像（图7（a）与图7（b）左列分别为图2（b）和2（c）中红框区域，图7（c）左列为图8（a）中红框区域）；图7右列为白帽覆盖率。

译

下载: | 高精图 | 低精图

图6 MSI白帽反演风速及其与ERA5风速产品的对比

Fig. 6 Comparison between whitecaps derived wind speed based on MSI images and wind speed product of ERA5

下载: | 高精图 | 低精图

图7 不同海洋环境动力对白帽的影响

Fig. 7 Influence of different marine environmental dynamics on whitecaps

下载: | 高精图 | 低精图

图8 Sentinel-2 MSI影像中白帽覆盖率与耀光变化分析

Fig. 8 Whitecap coverage and sunglint variation in Sentinel-2 MSI images

MSI 842 nm的开阔海洋观测数据中，虽然该数据的观测角度避开了强耀光反射（图2中数据的θ_m均大于40°），但在高空间分辨率光学影像中，海面耀光反射差异实际上仍存在，并对影像总信号有贡献（

Wu等，2022）。如图7（a）所示，海洋锋线（图中红色虚线）两侧呈现明显的水团图像特征差异，图7（a）中白帽覆盖率在不同水团之间也表现出类似的差异，虽然海水离水辐亮度差异可以忽略，但是白帽覆盖率差异和耀光反射差异同时存在；图7（b）和（c）反映了受海洋涡旋与内波的影响，白帽覆盖率会存在一定程度的升高。易于理解的是，相同风速条件下，内波的波峰处较其他位置，波浪更易于破碎而产生白帽（Thorpe等，1987；Kubryakov等，2021）。

译

海洋白帽是风生波浪破碎的产物，这种波浪破碎也受到海面锋面、涡旋、内波等各种其他海洋动力因素的调制，在高空间分辨率的MSI影像上，海洋白帽及其覆盖率所反映的海洋与大气背景动力特征能清晰展现。

译

3.4　耀光反射对白帽覆盖率估算的影响

海面耀光反射差异一直是海洋光学遥感研究关注的重点问题（

Cox和Munk，1954；Jackson和Alpers，2010；Kay等，2009；Harmel等，2018）。基于经典统计模型（Cox 和Munk，1954）估算的耀光反射率，虽可以准确描述粗空间分辨率卫星光学遥感数据（如250 m的MODIS数据）中耀光反射的强弱；但在高空间分辨率的海洋光学遥感数据中（如10 m的MSI数据），耀光反射规律会有所不同（Wu 等，2022），表现出较大的离散性，无法通过统计模型或简单的空间滤波进行剔除（Jackson和Alpers，2010）。本研究白帽覆盖率估算算法中的4 km × 4 km窗口滤波方法，本质是一种平滑处理（主要用于消除海水和大气的影响），不能直接从总体反射率中有效消除耀光反射的影响，因此无法处理强耀光反射的图像。

译

下载: | 高精图 | 低精图

（a） MSI R_rc （842 nm）图像（图像由6景MSI影像拼接而成，图中红色虚线为剖面线，红色方框部分数据对应图7（c）中图像）（b） MSI估算的白帽覆盖率结果（图中红色虚线为剖面线）

（a） MSI R_rc （842 nm） image （mosaic of 6 MSI images， the dotted line is the profile， the red box is correspond to the image in Fig. 7 （c））（b） Whitecap coverage estimated from MSI images （the dotted line is the profile）

为了阐明耀光反射对海面白帽覆盖率估算的影响，选择2021年6月2日的6景MSI数据（数据预处理与图2和图3一致）。如图8（a）中所示，θ_m从左到右逐渐减小（由32°减小到16°），耀光反射不断增强；图8（b）为算法估算的白帽覆盖率，与耀光强度空间变化相似；图8（c）为该区域同时段的ERA5风速产品，在无明显水团、涡旋等海洋动力环境的调制下，风速与白帽覆盖率的空间分异显然不同（图8（b））。此外，在相似风速下（8—12 m/s），图8（b）中强耀光反射下的白帽覆盖率较图3（a）的结果有一个数量级上的差异；对图8（a）与图8（b）做采样分析，以θ_m为横轴（图8（d）），不难发现样点的R_rc（842 nm）反射率与白帽覆盖率（W）的变化趋势一致。

译

海面耀光反射是白帽识别估算的主要干扰，本研究算法只适用于弱耀光反射下（θ_m大于40°）的MSI图像，如何有效消除耀光反射的影响，也是今后亟待关注的重点问题之一。

译

4　CZI数据中海洋白帽的响应

HY-1C/D是中国首个海洋水色业务卫星星座，CZI载荷能提供星下点空间分辨率为50 m、中心波长位于460 nm、560 nm、650 nm和825 nm的4波段海洋光学遥感数据（

沈亚峰等，2020）。本文所用的CZI数据均为L2A级的瑞利校正反射率产品（R_rc），由于像元的混合效应，相对较低风速条件下的白帽难以被CZI图像有效识别，随着风速的增加，海洋白帽逐渐被有效检测出，即CZI图像中检测出海洋白帽的最低风速要大于MSI图像。

译

图9为选取的开阔海洋CZI R_rc（825 nm）图像，图像都为海面弱耀光反射条件，海洋环境动力影响不明显，且涵盖不同风速条件下的影像，该区域同时段的ERA5风速均已经标注在R_rc（825 nm）图像上。对CZI图像上的海洋白帽进行识别并计算其覆盖率，图9（a）和图9（b）中海洋白帽不能被有效识别；在较高风速下（图9（c）、（d）、（e）、（f））海洋白帽能够从CZI图像上检测出。50 m空间分辨率CZI图像中可识别的海面白帽，其所对应的海面风速约9 m/s以上，随着风速的增加，其白帽覆盖率也随之增加（图9（c）、（d）、（e）、（f）），这也反映了CZI数据具有对高风速下的海面白帽覆盖率进行遥感监测与估算的应用潜力。

译

下载: | 高精图 | 低精图

图9 不同风速条件下HY-1C/D CZI影像中的白帽及其目视解译结果

Fig. 9 Whitecaps in HY-1C/D CZI images and corresponding visual interpretation results at different wind speed （U₁₀）

为了进一步分析阐明CZI图像中的海洋白帽反射率与背景海水反射率的差异，本研究对图9中4幅CZI图像（图9（c）、（d）、（e）、（f））的4个波段（460 nm、560 nm、650 nm和825 nm），对其海洋白帽反射率与背景海水反射率分别进行统计分析，结果如图10所示。CZI的4个波段图像中，白帽反射率均明显高于背景海水反射率，反射率差异分别为5.8%、7.8%、8.3%和7.3%。

译

图10 HY-1C/D CZI不同波段的海洋白帽与背景海水反射率差异

Fig. 10 Reflection difference between whitecaps and background seawater in different HY-1C/D CZI bands

下载: 原图 | 高精图 | 低精图

5　结论

海洋白帽是海面风生波浪破碎的直观表现，其区域性特征参数——白帽覆盖率，不仅与海面风速有着密切的统计关系，也会受到其他海洋环境动力的调制。海洋白帽对入射光具有较高的反射散射特征，是海洋光学遥感数据的大气校正一个重要环节。随着海洋光学遥感数据空间分辨率的提高，海洋白帽能被有效的探测，从而开展更为丰富的研究与应用探索。本文针对Sentinel-2 MSI数据，分析海面白帽反射率的信号贡献，利用一种用于弱耀光反射条件下的白帽覆盖率估算方法，开展MSI图像上白帽覆盖率反演。基于MSI图像估算的不同区域海面白帽覆盖率，与海面风速的变化具有良好的空间一致性和统计关系，基于白帽覆盖率估算的海面风速与ERA5再分析海洋风速产品对比验证表明，其相对偏差不大于1.6 m/s。此外，不同水团、锋面、涡旋、内波等其他海洋环境动力对海面白帽的调制作用，能显著改变遥感数据反演的白帽覆盖率；与此同时，MSI图像中耀光反射对海面白帽遥感的干扰与影响也不容忽视。最后，基于HY-1C/D卫星CZI载荷数据开展海面白帽的响应分析，结果表明：海面风速大于9 m/s时，在无其他海洋环境动力的调制下，海面白帽可被50 m空间分辨率的CZI载荷所探测，且白帽覆盖率随着风速的增加而增加；不同CZI波段图像中，白帽反射率与背景海水反射率有5.8%—8.3%的差异。这表明，HY-1C/D卫星CZI载荷可成为海洋白帽深入研究的可用数据源，能用于较高风速下的海面风速反演。总之，随着海洋光学遥感载荷空间分辨率的提升，通过海洋白帽的有效探测，能为海面风场及其他海洋环境动力监测提供新的技术方法参考，也必将在海面风场反演、海气交换、海洋动力环境研究等方面发挥越来越重要的作用。

译

致谢

致谢:感谢美国南佛罗里达大学海洋科学学院胡传民教授、武汉大学测绘遥感信息工程国家重点实验室田礼乔教授、中国科学院烟台海岸带研究所邢前国研究员对本研究的指导；感谢欧空局和欧洲中期天气预报中心提供的MSI和ERA5数据（https：//scihub.copernicus.eu/dhus/#/home，https：//cds.climate.copernicus.eu/cdsapp#！/dataset/reanalysis-era5-single-levels？tab=form）；感谢国家卫星海洋应用中心提供的CZI数据（http：//www.nsoas.org.cn/）。

参考文献（References）

Anguelova M D and Bettenhausen M H. 2019. Whitecap fraction from satellite measurements: Algorithm description[J]. Journal of Geophysical Research: Oceans, 124(3): 1827-1857 [DOI: 10.1029/2018JC014630] [百度学术]

Anguelova M D and Webster F. 2006. Whitecap coverage from satellite measurements: A first step toward modeling the variability of oceanic whitecaps[J]. Journal of Geophysical Research: Oceans, 111(C3) [DOI: 10.1029/2005JC003158] [百度学术]

Asher W E and Wanninkhof R. 1998. The effect of bubble-mediated gas transfer on purposeful dual-gaseous tracer experiments. Journal of Geophysical Research: Oceans, 103(C5): 10555-10560 [DOI: 10.1029/98jc00245] [百度学术]

Blanchard D C. 1983. The production, distribution, and bacterial enrichment of the sea-salt aerosol//Air-Sea Exchange of Gases and Particles. Dordrecht: Springer: 407-454 [DOI: 10.1007/978-94-009-7169-17] [百度学术]

Callaghan A H, Deane G B and Stokes M D. 2008a. Observed physical and environmental causes of scatter in whitecap coverage values in a fetch-limited coastal zone. Journal of Geophysical Research: Oceans, 113(C5): C05022 [DOI: 10.1029/2007jc004453] [百度学术]

Callaghan A, de Leeuw G, Cohen L and O'Dowd C D. 2008b. Relationship of oceanic whitecap coverage to wind speed and wind history. Geophysical Research Letters, 35(23): L23609 [DOI: 10.1029/2008gl036165] [百度学术]

Cox C and Munk W. 1954. Measurement of the roughness of the sea surface from photographs of the sun’s glitter. Journal of the Optical Society of America, 44(11): 838-850 [DOI: 10.1364/josa.44.000838] [百度学术]

Frouin R, Iacobellis S F and Deschamps P Y. 2001. Influence of oceanic whitecaps on the global radiation budget. Geophysical Research Letters, 28(8): 1523-1526 [DOI: 10.1029/2000gl012657] [百度学术]

Frouin R, Schwindling M and Deschamps P Y. 1996. Spectral reflectance of sea foam in the visible and near-infrared: in situ measurements and remote sensing implications. Journal of Geophysical Research: Oceans, 101(C6): 14361-14371 [DOI: 10.1029/96JC00629] [百度学术]

Garaba S P and Dierssen H M. 2018. An airborne remote sensing case study of synthetic hydrocarbon detection using short wave infrared absorption features identified from marine-harvested macro- and microplastics. Remote Sensing of Environment, 205: 224-235 [DOI: 10.1016/j.rse.2017.11.023] [百度学术]

Goddijn-Murphy L, Woolf D K and Callaghan A H. 2011. Parameterizations and algorithms for oceanic whitecap coverage. Journal of Physical Oceanography, 41(4): 742-756 [DOI: 10.1175/2010jpo4533.1] [百度学术]

Goddijn-Murphy L, Woolf D K, Callaghan A H, Nightingale P D and Shutler J D. 2016. A reconciliation of empirical and mechanistic models of the air-sea gas transfer velocity. Journal of Geophysical Research: Oceans, 121(1): 818-835 [DOI: 10.1002/2015jc011096] [百度学术]

Gordon H R and Wang M H. 1994. Influence of oceanic whitecaps on atmospheric correction of ocean-color sensors. Applied Optics, 33(33): 7754-7763 [DOI: 10.1364/ao.33.007754] [百度学术]

Harmel T, Chami M, Tormos T, Reynaud N and Danis P A. 2018. Sunglint correction of the Multi-Spectral Instrument (MSI)-SENTINEL-2 imagery over inland and sea waters from SWIR bands[J]. Remote Sensing of Environment, 204: 308-321 [DOI: 10.1016/j.rse.2017.10.022] [百度学术]

Hanson J L and Phillips O M. 1999. Wind sea growth and dissipation in the open ocean. Journal of Physical Oceanography, 29(8): 1633-1648 [DOI: 10.1175/1520-0485(1999)029<1633:WSGADI>2.0.CO;2] [百度学术]

Hu C M, Carder K L and Muller-Karger F E. 2000. Atmospheric correction of SeaWiFS imagery over turbid coastal waters: a practical method. Remote Sensing of Environment, 74(2): 195-206 [DOI: 10.1016/s0034-4257(00)00080-8] [百度学术]

Hu C M, Li X F, Pichel W G and Muller-Karger F E. 2009. Detection of natural oil slicks in the NW Gulf of Mexico using MODIS imagery. Geophysical Research Letters, 36(1): L01604 [DOI: 10.1029/2008gl036119] [百度学术]

Jackson C R and Alpers W. 2010. The role of the critical angle in brightness reversals on sunglint images of the sea surface[J]. Journal of Geophysical Research: Oceans, 115(C9) [DOI: 10.1029/2009JC006037] [百度学术]

Jia N and Zhao D L. 2019. The influence of wind speed and sea states on whitecap coverage. Journal of Ocean University of China, 18(2): 282-292 [DOI: 10.1007/s11802-019-3808-7] [百度学术]

Kay S, Hedley J D and Lavender S. 2009. Sun glint correction of high and low spatial resolution images of aquatic scenes: A review of methods for visible and near-infrared wavelengths[J]. Remote sensing, 1(4): 697-730 [DOI: 10.3390/rs1040697] [百度学术]

Kneizys F X, Shettle E P and Gallery W O. 1981. Atmospheric transmittance and radiance: the LOWTRAN 5 code//Proceedings Volume 0277, Atmospheric Transmission. Washington, D.C.: SPIE 277: 116-124 [DOI: 10.1117/12.931910] [百度学术]

Koepke P. 1984. Effective reflectance of oceanic whitecaps. Applied Optics, 23(11): 1816-1824 [DOI: 10.1364/AO.23.001816] [百度学术]

Kokhanovsky A A. 2004. Spectral reflectance of whitecaps. Journal of Geophysical Research: Oceans, 109(C5): C05021 [DOI: 10.1029/2003jc002177] [百度学术]

Kubryakov A A, Kudryavtsev V N and Stanichny S V. 2021. Application of Landsat imagery for the investigation of wave breaking. Remote Sensing of Environment, 253: 112144 [DOI: 10.1016/j.rse.2020.112144] [百度学术]

Lafon C, Piazzola J, Forget P, Le Calve O and Despiau S. 2004. Analysis of the variations of the whitecap fraction as measured in a coastal zone. Boundary-Layer Meteorology, 111(2): 339-360 [DOI: 10.1023/B:BOUN.0000016490.83880.63] [百度学术]

Lu Y C, Sun S J, Zhang M W, Murch B and Hu C M. 2016. Refinement of the critical angle calculation for the contrast reversal of oil slicks under sunglint. Journal of Geophysical Research: Oceans, 121(1): 148-161 [DOI: 10.1002/2015JC011001] [百度学术]

Melville W K. 1996. The role of surface-wave breaking in air-sea interaction. Annual Review of Fluid Mechanics, 28(1): 279-321 [DOI: 10.1146/annurev.fl.28.010196.001431] [百度学术]

Monahan E C, Fairall C W, Davidson K L and Boyle P J. 1983. Observed inter-relations between 10 m winds, ocean whitecaps and marine aerosols. Quarterly Journal of the Royal Meteorological Society, 109(460): 379-392 [DOI: 10.1002/qj.49710946010] [百度学术]

Monahan E C and Muircheartaigh I. 1980. Optimal power-law description of oceanic whitecap coverage dependence on wind speed. Journal of Physical Oceanography, 10(12): 2094-2099 [DOI: 10.1175/1520-0485(1980)010<2094:OPLDOO>2.0.CO;2] [百度学术]

Monahan E C and Zietlow C R. 1969. Laboratory comparisons of fresh-water and salt-water whitecaps. Journal of Geophysical Research, 74(28): 6961-6966 [DOI: 10.1029/JC074i028p06961] [百度学术]

Moore K D, Voss K J and Gordon H R. 2000. Spectral reflectance of whitecaps: their contribution to water-leaving radiance. Journal of Geophysical Research: Oceans, 105(C3): 6493-6499 [DOI: 10.1029/1999jc900334] [百度学术]

Myrhaug D and Holmedal L E. 2008. Effects of wave age and air stability on whitecap coverage. Coastal Engineering, 55(12): 959-966 [DOI: 10.1016/j.coastaleng.2008.03.005] [百度学术]

Nicolas J M, Deschamps P Y and Frouin R. 2001. Spectral reflectance of oceanic whitecaps in the visible and near infrared: aircraft measurements over open ocean. Geophysical Research Letters, 28(23): 4445-4448 [DOI: 10.1029/2001gl013556] [百度学术]

Paget A C, Bourassa M A and Anguelova M D. 2015. Comparing in situ and satellite-based parameterizations of oceanic whitecaps. Journal of Geophysical Research: Oceans, 120(4): 2826-2843 [DOI: 10.1002/2014jc010328] [百度学术]

Phillips O M. 1985. Spectral and statistical properties of the equilibrium range in wind-generated gravity waves. Journal of Fluid Mechanics, 156: 505-531 [DOI: 10.1017/s0022112085002221] [百度学术]

Pivaev P D, Kudryavtsev V N, Korinenko A E and Malinovsky V V. 2021. Field observations of breaking of dominant surface waves. Remote Sensing, 13(16): 3321 [DOI: 10.3390/rs13163321] [百度学术]

Raizer V Y, Novikov V M and Bocharova T Y. 1994. The geometrical and fractal properties of visible radiances associated with breaking waves in the ocean. Annales Geophysicae, 12(12): 1229-1233 [DOI: 10.1007/s00585-994-1229-3] [百度学术]

Romero L. 2019. Distribution of surface wave breaking fronts. Geophysical Research Letters, 46(17/18): 10463-10474 [DOI: 10.1029/2019gl083408] [百度学术]

Salisbury D J, Anguelova M D and Brooks I M. 2014. Global distribution and seasonal dependence of satellite-based whitecap fraction. Geophysical Research Letters, 41(5): 1616-1623 [DOI: 10.1002/2014gl059246] [百度学术]

Shen Y F, Liu J Q, Ding J, Jiao J N, Sun S J and Lu Y C. 2020. HY-1C COCTS and CZI observation of marine oil spills in the South China Sea. Journal of Remote Sensing, 24(8): 933-944 [百度学术]

沈亚峰, 刘建强, 丁静, 焦俊男, 孙绍杰, 陆应诚. 2020. 海洋一号C星光学载荷对海面溢油的识别能力分析. 遥感学报, 24(8): 933-944 [DOI: 10.11834/jrs.20209475] [百度学术]

Thorpe S A, Belloul M B and Hall A J. 1987. Internal waves and whitecaps. Nature, 330(6150): 740-742 [DOI: 10.1038/330740a0] [百度学术]

Wen Y S, Wang M Q, Lu Y C, Sun S J, Zhang M W, Mao Z H, Shi J and Liu Y X. 2018. An alternative approach to determine critical angle of contrast reversal and surface roughness of oil slicks under sunglint. International Journal of Digital Earth, 11(9): 972-979 [DOI: 10.1080/17538947.2018.1470687] [百度学术]

Wu X Y, Lu Y C, Jiao J N, Ding J, Fu W X and Qian W X. 2022. Using sea wave simulations to interpret the sunglint reflection variation with different spatial resolutions. IEEE Geoscience and Remote Sensing Letters, 19: 1501304 [DOI: 10.1109/LGRS.2020.3033700] [百度学术]

Zhang H G, Lou X L, Li Y, Shi A Q, Li D L and Fu B. 2015. Whitecap features induced by submarine sand waves in stereo optical imagery. Journal of Geophysical Research: Oceans, 120(9): 6225-6233 [DOI: 10.1002/2015jc010947] [百度学术]

摘要

Abstract

Oceanic whitecaps

generated from wind-wave breaking process

are the medium of air-sea exchange and the indicator of sea surface state. Due to the intense reflection and scattering of incident light

whitecaps can be effectively recognized and distinguished in in-situ photos or videos. Whitecap coverage (

)

defined as the proportion of ocean surface covered by whitecaps

is an important parameter for the quantification of whitecaps. Of course

oceanic whitecaps can also be discriminated in high spatial resolution optical remote sensing images

such as Sentinel-2 MSI and Landsat-8 OLI data. This can provide a new research direction in marine environment observation

and may be further used for wind speed monitoring. However

how to estimate oceanic whitecap coverage from these optical remote sensing images is still a challenge. In this study

the formula of whitecap coverage is obtained by converting the form of radiative equation related with the constant and image reflectance of Sentinel-2 MSI data under some assumptions. The background signal of seawater and atmosphere was eliminated by data optimization

and the signal of whitecaps can be distinguished using regional filtering method. Then

whitecap coverage can be estimated. The identification and estimation results indicate that whitecap coverage derived from Sentinel-2 MSI images are consistent with previous studies using in situ observations in the order of magnitude

and can invert sea surface wind speed using a statistical model. Coarse spatial resolution wind speed images converted from MSI inversion were validated by ERA5 wind speed products from European Centre for Medium-Range Weather Forecasts (ECMWF). In addition

whitecap coverage can imply the modulation of other marine environmental dynamic factors (

e.g

water mass

ocean fronts

ocean eddies and internal waves). Moreover

it should be noted that sunglint reflection is a non-negligible issue for optical remote sensing of oceanic whitecaps whose signal should be effectively eliminated. Examining the MSI-derived results to HY1-C/D CZI

it indicates that whitecaps can be identified in CZI images when wind speed is greater than 9 m/s

and the reflectance difference between whitecaps and background seawater is 5.8%—8.3%. We hope the above results can be used to improve the accuracy of atmospheric correction

and provide new reference for using high spatial resolution optical remote sensing in sea surface wind speed estimation and marine environmental dynamic factors monitoring. This will hopefully expand the research and application fields of ocean color remote sensing.

关键词

光学遥感海洋白帽耀光HY-1C/DSentinel-2MSICZI

Keywords

optical remote sensingoceanic whitecapssunglintHY-1C/DSentinel-2MSICZI

references

Gordon H R and Wang M H. 1994. Influence of oceanic whitecaps on atmospheric correction of ocean-color sensors. Applied Optics, 33(33): 7754-7763 [DOI: 10.1364/ao.33.007754http://dx.doi.org/10.1364/ao.33.007754]

Jia N and Zhao D L. 2019. The influence of wind speed and sea states on whitecap coverage. Journal of Ocean University of China, 18(2): 282-292 [DOI: 10.1007/s11802-019-3808-7http://dx.doi.org/10.1007/s11802-019-3808-7]

Koepke P. 1984. Effective reflectance of oceanic whitecaps. Applied Optics, 23(11): 1816-1824 [DOI: 10.1364/AO.23.001816http://dx.doi.org/10.1364/AO.23.001816]

Kokhanovsky A A. 2004. Spectral reflectance of whitecaps. Journal of Geophysical Research: Oceans, 109(C5): C05021 [DOI: 10.1029/2003jc002177http://dx.doi.org/10.1029/2003jc002177]

Melville W K. 1996. The role of surface-wave breaking in air-sea interaction. Annual Review of Fluid Mechanics, 28(1): 279-321 [DOI: 10.1146/annurev.fl.28.010196.001431http://dx.doi.org/10.1146/annurev.fl.28.010196.001431]

Monahan E C and Zietlow C R. 1969. Laboratory comparisons of fresh-water and salt-water whitecaps. Journal of Geophysical Research, 74(28): 6961-6966 [DOI: 10.1029/JC074i028p06961http://dx.doi.org/10.1029/JC074i028p06961]

Myrhaug D and Holmedal L E. 2008. Effects of wave age and air stability on whitecap coverage. Coastal Engineering, 55(12): 959-966 [DOI: 10.1016/j.coastaleng.2008.03.005http://dx.doi.org/10.1016/j.coastaleng.2008.03.005]

Pivaev P D, Kudryavtsev V N, Korinenko A E and Malinovsky V V. 2021. Field observations of breaking of dominant surface waves. Remote Sensing, 13(16): 3321 [DOI: 10.3390/rs13163321http://dx.doi.org/10.3390/rs13163321]

Romero L. 2019. Distribution of surface wave breaking fronts. Geophysical Research Letters, 46(17/18): 10463-10474 [DOI: 10.1029/2019gl083408http://dx.doi.org/10.1029/2019gl083408]

Shen Y F, Liu J Q, Ding J, Jiao J N, Sun S J and Lu Y C. 2020. HY-1C COCTS and CZI observation of marine oil spills in the South China Sea. Journal of Remote Sensing, 24(8): 933-944

沈亚峰, 刘建强, 丁静, 焦俊男, 孙绍杰, 陆应诚. 2020. 海洋一号C星光学载荷对海面溢油的识别能力分析. 遥感学报, 24(8): 933-944 [DOI: 10.11834/jrs.20209475http://dx.doi.org/10.11834/jrs.20209475]

Thorpe S A, Belloul M B and Hall A J. 1987. Internal waves and whitecaps. Nature, 330(6150): 740-742 [DOI: 10.1038/330740a0http://dx.doi.org/10.1038/330740a0]

文章被引用时，请邮件提醒。

提交

中国近海绿潮生物量的卫星光学遥感估算

Sentinel-1和Sentinel-2协同反演稀疏非光合植被覆盖度

基于哨兵卫星数据的巴丹吉林—腾格里沙漠辐射准不变场选取与评价