基于同期光学与微波遥感的茅尾海及其入海口水体悬浮物反演
Inversion of total suspended matter concentration in Maowei Sea and its estuary, Southwest China using contemporaneous optical data and GF SAR data
通讯作者: 李伟(1984-),男,硕士,主要从事区域地质及海岸带地质调查方面的研究。Email:670095930@qq.com。
责任编辑: 陈理
收稿日期: 2021-03-31 修回日期: 2021-07-30
基金资助: |
|
Received: 2021-03-31 Revised: 2021-07-30
作者简介 About authors
丁波(1978-),男,硕士,主要从事遥感地质及海岸带地质方面的研究。Email:
悬浮物是近海海湾及潟湖水质的重要影响因素之一。遥感技术能够准确快速地反演水体的悬浮物浓度,但鲜有同时利用2种不同类型的数据源反演同一研究区的悬浮物浓度。该文基于同一时期(2019年9月)的GF-1C光学影像和GF-3微波数据,采用双波段比值算法和Cloude-Pottier分解算法分别对原始影像进行图像分割和目标提取,并在此基础上,结合同期野外采样数据,利用三次多项式回归算法,开展了茅尾海及其入海口水体悬浮物反演。精度分析显示,GF-1C反演模型相关系数(R2)、均方根误差和平均相对误差分别为0.88,130.25 mg/L和9.65%; 而GF-3反演模型对应结果分别为0.61,230.87 mg/L和15.13%,研究表明,GF-1C光学遥感反演精度(90.35%)要好于GF-3微波遥感反演结果(84.87%),但2种反演结果分布具有较高的相似性和一致性,进一步表明基于2种不同数据源建立的反演模型能够为茅尾海悬浮物反演和海岸带环境监测提供参考。
关键词:
Total suspended matter concentration (TSMC) is one of the important factors influencing water bodies in coastal gulfs and lagoons. The rapid and accurate TSMC inversion can be obtained using remote sensing data. However, it is scarce to conduct TSMC inversion using two different data sources at the same time. This study conducted the inversion of TSMC in Maowei Sea and its estuary based on two data sources. Specifically, this study carried out image segmentation and object extraction using the dual-band ratio algorithm and the Cloude-Pottier target decomposition algorithm, respectively based on GF-1C optical images and GF-3 SAR data of September 2019. Meanwhile, contemporaneous field sample data were utilized. Then, the TSMC inversion was performed using the cubic polynomial regression algorithm. As revealed by the accuracy analysis, the fitting degree (R2), root mean square error, and mean relative percentage error of the GF-1C-based inversion model were 0.88, 130.25 mg/L, and 9.65%, respectively, while those of the GF-3-based inversion model were 0.61, 230.87 mg/L, and 15.13%, respectively. These indicate that the GF-1C-based TSMC inversion had a higher inversion accuracy (90.35%) than the GF-3-based TSMC inversion (84.87%). However, the inversion results of the two models showed highly similar distribution patterns. This further indicates that the inversion models established using two different data sources in this study can serve as references for TSMC inversion of Maowei Sea and its estuary and for the environmental monitoring in coastal zones.
Keywords:
本文引用格式
丁波, 李伟, 胡克.
DING Bo, LI Wei, HU Ke.
0 引言
悬浮物是指悬浮于海水中的固态物质,包括有机成分(如碳水化合物、蛋白质、生物残骸、排泄物和分解物等)和无机组分(如陆源矿物颗粒、悬浮泥沙、黏土以及人类活动产生的颗粒和溶解物等)[1,2]。悬浮物浓度是衡量近海海水清澈度和其他成分的重要参数,尤其是位于海湾、河流入海口、近海等区域,受陆源物质、人类活动和海洋动力等因素的影响[3],较高的悬浮物浓度极大地降低了光线的穿透能力,不仅直接影响了水体生态条件比如浮游生物繁育[4,5]、生蚝养殖、鱼类生长,以及污染物扩散等[6],同时也影响着河口变化及海岸带冲淤过程。悬浮物浓度高低,在光学遥感上表现为光谱反射率、水体清澈度和水色等光学特性[7,8]; 而在微波遥感上反映为后向散射系数、粗糙程度和介电常数等特征。
悬浮物的监测方法主要包括常规实地调查和遥感反演。常规方法是乘船逐点采样-化验成分-分析数据[9],这种方法调查速度慢、周期长,费用高,且只能获得在时间、空间分布上都很离散的少量点数据。遥感反演得益于卫星技术的发展,传感器具有快速、高效且重访周期短等优势,能够快速获取任一区域内水体遥感影像,实时监测其悬浮物的空间分布和动态变化,克服了常规方法的不足,因此受到广泛关注。
光学遥感悬浮物反演理论依据是水体的光谱反射率,不同物质成分的水体拥有不同的光谱反射率,一般来说,纯水的光谱反射率较低,而悬浮物的光谱反射率较高[10,11,12]。过去的20 a里,众多的学者在悬浮物浓度反演方面开展了大量的研究工作并取得了可观的研究成果,反演算法涵盖从单波段[7,13]到多波段[2, 13-17],从低分辨率[2, 14-16, 18-20]到高分辨率[8,12,17,21]。而微波遥感反演悬浮物的理论核心是SAR图像的后向散射系数,一个包含物质微观结构(也包含浑浊水体或悬浮物)表面细节(粗糙度、几何形态、复介电常数等)的基本量 [22,23] 。根据大量学者的研究,近海悬浮物浓度影响了海水的组分、粒径以及海水浑浊度等物理属性,而海水的组分和粒径跟水体的后向散射系数有密切的关系并对其产生着影响,这一特征最终在雷达影像上得到一定的反映[24,25,26,27,28]。
近些年来,学者更多利用多光谱遥感开展各类水体的悬浮物浓度反演,而较少利用微波遥感开展相关研究。本文针对同一目标海域,利用2种同时期获取的不同类型的国产数据源(GF-1C和GF-3)开展悬浮物遥感反演分析,并结合野外采样数据进行对比,讨论并分析了2种反演结果的可靠性和准确性。
1 研究区概况及数据源
1.1 研究区概况
茅尾海地处广西钦州湾,是一个半封闭内海(N21°30'~21° 55',E108°25'~108°45'),由一条狭窄的水道与北部湾相连(图1)。潟湖海岸线长约120 km,东西宽约14 km,南北纵深约17 km,面积约134 km2。平均水深为0.1~5 m,最深处可达29 m。
图1
图1
研究区域与野外采样点分布(2019年GF-1C多光谱波段与全色波段融合影像)
Fig.1
Map of the study area and sampling sites
茅尾海上游汇集大榄江、钦江和茅岭江等3条常年性河流,均发源于上游山区,地质构造复杂,地层主要包括: 下层为下古生界志留系,表层为覆盖较厚的第四系,3条主要河流每年通过径流为茅尾海输入大量陆源泥沙。茅尾海内海常年海面平静,泥沙充沛,饵料丰富,是广西重要的海洋牧场,也是中国著名的生蚝养殖基地,同时兼有对虾、青蟹、鱼类及贝类养殖; 因其内海海域开阔,岛屿众多,也是风景优美的旅游胜地。茅尾海入海口两侧分布着2处重要的基础设施: 东侧为钦州港,包含多个人工岛和防护堤,西侧为防城港核电站,修建了2条短的平行防沙堤和进水通道,以及2条狭长的排水通道。
1.2 采样数据
表1 研究区悬浮物浓度野外采样数据
Tab.1
取样 编号 | 悬浮物浓度/(mg·L-1) | 取样 日期 | 平均潮 高/cm | ||
---|---|---|---|---|---|
0.5 m深度 | 5 m深度 | 10 m深度 | |||
1 | 1 527 | 822 | 1 120 | 2019-09-18 | 324 |
2 | 1 038 | 1 276 | 1 201 | ||
3 | 924 | 1 262 | 1 020 | ||
4 | 1 030 | 1 064 | 989 | ||
5 | 892 | 1 180 | 1 104 | 2019-09-19 | 345 |
6 | 1 124 | 995 | 1 418 | ||
7 | 1 216 | 1 230 | 962 | ||
8 | 1 036 | 1 213 | 1 064 | ||
9 | 1 074 | 1 464 | 948 | ||
10 | 1 248 | 1 394 | 1 174 | ||
11 | 1 090 | 1 177 | 872 | 2019-09-20 | 367 |
12 | 1 120 | 1 226 | 1 186 | ||
13 | 1 149 | 1 374 | 1 281 | ||
14 | 1 352 | 964 | 1 053 | ||
15 | 1 103 | 1 204 | 1 204 |
1.3 遥感数据集
高分专项是我国于2010年5月经由国务院审议并批准实施的高分辨率对地观测系统重大专项之一[30],目前已发射了14个卫星系列(GF-1—GF-14),传感器类型从光学到雷达,轨道从太阳同步到地球同步,拍摄能力从全色、多光谱到高光谱等,构成了一个具有高空间分辨率、高时间分辨率和高光谱分辨率能力的对地观测系统,其中7个卫星系列已投入民用。
GF-1是高分专项的首发光学遥感卫星,包括后续发射的姊妹星GF-1B,C,D,全色影像空间分辨率为2 m; 多光谱影像空间分辨率为8 m,光谱范围包括蓝光波段(0.45~0.52 μm)、绿光波段(0.52~0.59 μm)、红光波段(0.63~0.69 μm)以及近红外波段(0.77~0.89 μm)。GF-3是我国首颗C频段多极化合成孔径雷达(synthetic aperture Radar,SAR)卫星,传感器包括12种成像模式,能够高时效地实现不同应用模式下获取1~500 m空间分辨率、10~650 km幅宽的微波遥感数据。本文所用数据具体参数见表2。
表2 研究区遥感影像详细信息
Tab.2
数据 类型 | 传感 器 | 分辨 率/m | 时相 | 参数 | ||
---|---|---|---|---|---|---|
光学 | GF-1C | 8 | 2019-9-26 11: 44: 59 am | 多光谱波段 | 增益 | 偏差 |
蓝光波段 | 0.029 0 | 0 | ||||
绿光波段 | 0.038 2 | 0 | ||||
红光波段 | 0.042 1 | 0 | ||||
近红外波段 | 0.036 4 | 0 | ||||
SAR | GF-3 | 8 | 2019-9-4 11: 01: 40 am | C波段: (VH,VV) 分辨率: 4.5m×4.8m; 入射角: 42.787 378°~ 47.954 950° |
上述2类数据均可用于悬浮物浓度反演。根据GF-1和GF-3的数据覆盖同步性,最终选取时相为2019年9月份的GF-1多光谱影像和GF-3 SAR影像各1景作为研究区悬浮物反演遥感源数据。同时收集了区内的地形数据及卫星参数信息作为辅助数据,为影像预处理提供基础参考。
2 研究方法
2.1 光学遥感悬浮物水体提取
图2
图2
不同悬浮物浓度水体与其光谱反射率关系
Fig.2
Relationship between TSMC and spectral reflectance
式中: ρR 和ρNIR 分别为GF-1的红光波段和近红外波段反射率; Gain和Bias分别为增益和偏差。
2.2 微波遥感悬浮物水体提取
根据雷达工作原理,自然界大多数地物均可以视为粗糙程度不同的粗糙表面,在雷达上表现为地物的后向散射回波强度,图像上呈现出不同的灰度值,常用后向散射截面积σ或后向散射系数σ0来表达,其公式为:
式中: λ为波长; θ为入射角; Ρ为极化方式; φ为方位角; ε为复介电常数; Γ1为表面粗糙度; Γ2为次表面粗糙度; V为体散射系数。
1)降噪处理(多视与滤波)。斑点噪声是SAR系统的固有缺陷,对SAR图像在实际应用过程中的目标识别、信息提取和分类等产生了严重影响和干扰。地距分辨率是指方位向和距离向相同的地面分辨率。为了得到地距分辨率同时减少SAR图像的噪声干扰,需要对原始SAR图像进行多视和滤波处理。计算公式为:
式中: Rgrr为地距分辨率; r为斜距分辨率。由于距离分辨率随入射角的变化而变化,多视处理尽管提升了数据的辐射分辨率,但降低了SAR数据的空间分辨率。
式中: αi为第i个散射机制类型; Pi为第i个散射机制发生概率; χi为极化相干矩阵的特征值。
基于VH+VV增强后向散射强度图反映出的水体悬浮物浓度微弱信息,再利用已有现状海岸线进行海陆分割,提取目标水体的基础上,采用Hpe-α对目标水体进行空间重分类[44],最终获取了目标海域悬浮物浓度灰度图。
3 结果与分析
3.1 回归模型建立
图3
图3
悬浮物浓度与等效反射率之间的回归关系
Fig.3
Regression relationships of measured TSMC and equivalent remote sensing reflectance
3.2 模型验证
式中: xr,i为第i点悬浮物浓度反演值; xm,i为第i点悬浮物浓度测量值。
基于已建立的评价指标,将2种反演模型分析得到的悬浮物浓度与另外8组采样点实测数据进行精度验证与评估,如表3所示。精度结果表明,GF-1C悬浮物反演精度(均方根误差为130.25 mg/L,平均相对误差为9.65%)明显好于GF-3悬浮物反演结果(均方根误差为230.87 mg/L,平均相对误差高达15.13%)。由于与10号验证点对应的GF-3反演值出现异常,原始数据发现该处存在坏点,故该异常点在对GF-3反演悬浮物精度评价时不参与验证。
表3 反演模型精度验证
Tab.3
反演模型 | 验证点 数/个 | Δrmse/ (mg·L-1) | 相对误差/% | ||
---|---|---|---|---|---|
最大值 | 最小值 | Δmrpe | |||
GF-1C | 8 | 130.25 | 17.70 | 2.15 | 9.65 |
GF-3 | 7 | 230.87 | 34.14 | 5.57 | 15.13 |
3.3 结果分析
3.3.1 研究区悬浮物分布规律
根据已建立的悬浮物反演理论模型,分别对研究区内2种数据源的悬浮物分布进行了反演。图4结果表明,2种方法反演出的结果具有高度的相似性和一致性: 茅尾海内海及入海口两侧悬浮物浓度总体较高,均超过1 000 mg/L,其分布规律呈现由东北向西南依次递减。结合研究区同季节历史资料[45]分析其原因: 在茅尾海内海上游3条主要河流每年向其注入泥沙量高达60万t,而在连接外海通道处分布着大量的天然岛屿和人工养殖蚝排蚝桩,严重阻碍了内海水体与北部湾海域的水体交换; 而在茅尾海入海口,悬浮物浓度分布表现为两侧偏高,中间偏低,原因在于入海口中部海域每年都进行人工航道清淤,但两侧人为活动和船只进出对水体扰动频繁,泥沙沉积过程和海洋动力作用如洋流、潮汐等在该区域的影响微弱。
图4
图4
不同数据源的悬浮物浓度分布
Fig.4
Distributions of retrieved TSMC based on different remote sensing images
3.3.2 反演结果对比
根据反演结果,将悬浮物浓度进一步划分为3个区间: <1 100 mg/L,[1 100,1 200) mg/L,≥1 200 mg/L,分别代表悬浮物浓度低、中、高3个等级,再将实测点分布与悬浮物浓度等级区间进行相对位置关系讨论,如图5所示。
图5
图5
悬浮物浓度反演结果与实测结果对比
Fig.5
Comparisons of retrieved TSMCs and in situ measurements
从精度评价结果看,研究区悬浮物浓度光学遥感反演结果精度(90.35%)明显高于微波遥感反演结果(84.87%),表明光学遥感反演结果要明显好于微波遥感反演结果。
1)从反演效果(图4)看,光学遥感反演结果灰度值过渡均匀,噪声点较少,效果明显好于微波遥感反演结果。
2)从反演结果分区和实测点相对位置(图5)看,光学遥感反演结果高值区和低值区与对应的实测数据完全吻合,仅在悬浮物浓度中值区,存在1处错误; 而微波遥感反演结果各个浓度分区与实测点均存在偏差,其中反演结果中值区和低值区各有1处错误,而在高值区有2处错误。
3)从取样时间和卫星数据获取时间节点看,由于不同传感器对同一区域数据获取时间以及野外取样时间三者无法保证高度一致,光学遥感反演结果精度要高于微波遥感反演结果,这与取样时间和光学遥感影像获取时间间隔(6~8 d)要小于和微波影像获取时间间隔(14~16 d)存在一定的相关性。
4 结论
近海及海湾地区悬浮物浓度的分布对于海岸带综合地质研究和生态环境保护具有重要的意义。利用遥感技术能够快速准确地反演其分布特征,本文利用2种不同的遥感数据源,采用不同的提取算法,并结合同期实地采样数据,准确地反演了茅尾海及其入海口悬浮物的分布。从技术方法看,2种手段均能准确地反演茅尾海悬浮物的分布特征,且反演结果表现出了高度的一致性; 从反演结果看,光学遥感反演结果要优于微波遥感反演结果,但微波遥感反演突破了光学遥感反演诸多的干扰因素,是光学遥感反演的有力技术补充。
另外,需要指出的是,由于野外取样数据有限,反演结果的准确性需要更多野外采样数据进行实验验证。
(责任编辑: 陈理)
参考文献
Using MODIS Terra 250 m imagery to map concentrations of total suspended matter in coastal waters
[J].DOI:10.1016/j.rse.2004.07.012 URL [本文引用: 3]
Influence of suspended particle concentration,composition and size on the variability of inherent optical properties of the Southern North Sea
[J].DOI:10.1016/j.csr.2012.01.007 URL [本文引用: 1]
Effects of spatial and temporal variability of turbidity on phytoplankton blooms
[J].DOI:10.3354/meps254111 URL [本文引用: 1]
The development of a new optical total suspended matter algorithm for the Chesapeake Bay
[J].DOI:10.1016/j.rse.2011.12.018 URL [本文引用: 1]
GEOTRACES:An international study of the global marine biogeochemical cycles of trace elements and their isotopes
[J].DOI:10.1016/j.chemer.2007.02.001 URL [本文引用: 1]
Optical algorithms at satellite wavelengths for total suspended matter in tropical coastal waters
[J].DOI:10.3390/s8074165 URL [本文引用: 2]
Retrieval of total suspended matter concentrations from high resolution WorldView-2 imagery:A case study of inland rivers
[C]//
Suspended solids concentration in wastewater:Influence of sampling conditions
[J].DOI:10.1080/1573062X.2011.630094 URL [本文引用: 1]
The spectral reflectance of water-mineral mixtures at low temperatures
[J].
USGS Spectral Library Version 7[Z].U.S
Detection and delineation of localized flooding from WorldView-2 multispectral data
[J].DOI:10.3390/rs71114853 URL [本文引用: 2]
A SWIR based algorithm to retrieve total suspended matter in extremely turbid waters
[J].DOI:10.1016/j.rse.2015.06.022 URL [本文引用: 2]
Retrieval of total suspended matter from MODIS 250 m imagery in the Bohai Sea of China
[J].DOI:10.1007/s10872-012-0129-5 URL [本文引用: 2]
The use of CBERS(China-Brazil earth resources satellite) to trace the dynamics of total suspended matter at an urbanized coastal area
[J].DOI:10.1590/s1679-87592017135006502 URL
Comparison of remote sensing data,model results and in situ data for total suspended matter(TSM) in the southern Frisian lakes
[J].DOI:10.1016/S0048-9697(00)00679-3 URL [本文引用: 1]
基于GF-1与实测光谱数据鄱阳湖丰水期总悬浮物浓度反演
[J].DOI:10.6046/gtzyyg.2019.01.14 [本文引用: 3]
Retrieval of concentration of total suspended matter from GF-1 satellite and field measured spectral data during flood period in Poyang Lake
[J].DOI:10.6046/gtzyyg.2019.01.14 [本文引用: 3]
A simple empirical band-ratio algorithm to assess suspended particulate matter from remote sensing over coastal and inland waters of Vietnam:Application to the VNREDSat-1/NAOMI sensor
[J].
Retrieval of Chlorophyll a,suspended solids,and colored dissolved organic matter in Tokyo Bay using ASTER data
[J].DOI:10.1016/j.rse.2005.05.016 URL
Modeling and multi-temporal characterization of total suspended matter by the combined use of Sentinel 2-MSI and Landsat 8-OLI Data:The Pertusillo Lake case study(Italy)
[J].DOI:10.3390/rs12010001 URL [本文引用: 2]
Satellite remote sensing of water turbidity
[J].DOI:10.1080/02626668009491950 URL [本文引用: 1]
Tests of backscatter coefficient measurement using broadband pulses
[J].DOI:10.1109/58.238114 URL [本文引用: 1]
On the estimation of backscatter coefficients using single-element focused transducers
[J].DOI:10.1121/1.3557036 URL [本文引用: 1]
InSAR技术及其应用中的若干问题
[J].
In SAR technology and its key problems in applications
[J].
Suspended solids analysis using ERTS-A data
[J].DOI:10.1016/0034-4257(74)90039-X URL [本文引用: 1]
基于实测光谱的杭州湾悬浮物浓度杨反演模式
[J].
Remote sensing inversion mode of suspended particles concentration in Hangzhou Bay based on in situ measurement spectrum
[J].
Deriving total suspended matter concentration from the near-infrared-based inherent optical properties over turbid waters:A case study in Lake Taihu
[J].DOI:10.3390/rs10020333 URL [本文引用: 1]
基于极化SAR的河流有机物污染监测研究
[J].
Monitoring of organic pollutants in river based on polarimetric SAR
[J].
DD2014-03海岸带地质环境调查评价规范
[S].
DD2014-03 Specification for survey and evaluation of coastal geological environment
[S].
USGS digital spectral library splib06a
[Z].
Extraction of spatial and temporal patterns of concentrations of chlorophyll-a and total suspended matter in Poyang Lake using GF-1 satellite data
[J].DOI:10.3390/rs12010001 URL [本文引用: 2]
基于遥感方法反演珠江三角洲西江干流悬浮泥沙分布研究
[J].
Study on inversion model for the suspended sediment concentration in the Pear River delta using remote sensing technology
[J].
A simple speckle smoothing algorithm for synthetic aperture Radar images
[J].
Improved sigma filter for speckle filtering of SAR imagery
[J].DOI:10.1109/TGRS.36 URL [本文引用: 1]
Speckle reduction in multipolarization,multifrequency SAR imagery
[J].DOI:10.1109/36.135815 URL [本文引用: 1]
Scattering-model-based Speckle filtering of polarimetric SAR data
[J].DOI:10.1109/TGRS.2005.859338 URL [本文引用: 1]
Gamma-WMAP:A wavelet based statistical speckle filter for SAR images
[C]//
Noise filtering of SAR interferometric phase based on wavelet transform
[C]//
一种基于小波分析的SAR图像斑点噪声滤波算法
[J].
Speckle filtering method of SAR images based on wavelet analysis
[J].
A review of target decomposition theorems in Radar polarimetry
[J].DOI:10.1109/36.485127 URL [本文引用: 1]
基于Cloude-Pottier目标分解和聚合的层次聚类算法的全极化SAR数据的非监督分类算法研究
[J].
An unsupervised classification for fully polarimetric SAR data using Cloude-Pottier decomposition and agglomerative hierarchical clustering algorithm
[J].
一种利用Cloude-Pottier分解和极化白化滤波的全极化SAR图像分类算法
[J].
An unsupervised Wishart classification for fully polarimetric SAR image based on Cloude-Pottier decomposition and polarimetric whitening filter
[J].
基于Radarsat-2目标极化分解H-α特征空间的海洋悬浮物浓度监测研究
[J].
Monitoring marine suspended particulate matter concentration using full polarimetric SAR data of Radarsat-2
[J].
/
〈 |
|
〉 |
