|
珊瑚礁是全球气候变化的重要显示器，是热带海洋最突出、最具代表性的生态系统，具有惊人的生物多样性和极高的初级生产力，对于维持着丰富的生物多样性，维护海洋生态环境至关重要。当珊瑚面临巨大的生存压力，珊瑚礁就会发生退化现象，我们称之为珊瑚礁白化，白化后的珊瑚礁可能恢复也可能死亡，白化严重的珊瑚将无法恢复，对珊瑚礁生态系统造成很大的破坏。我国南海珊瑚目前处于总体退化状态，基于卫星遥感可以有效地实现珊瑚礁区的底质变化状态监测，具有重要的研究价值。 , J* ~- t4 c, ?' Q9 o. ]/ e
右侧为赵述岛，供图/朱卫东 早在上世纪七十年代，有学者就提出了利用遥感手段对珊瑚礁进行研究的观点。21世纪以来，珊瑚礁遥感影像分类的相关探究逐渐受到学者的重视。基于遥感影像的珊瑚礁信息提取的主要方法有现场光谱数据样本法、决策树模型法、最大似然法、面向对象法、支持向量机法。现场光谱数据样本法以先验知识为基础，选取样本光谱进行训练，并分析地物光谱和融合信息，但该方法对现场调查要求较高。决策树模型法基于图像光谱带的多个光谱决策规则建立模型，完成了地貌识别和底质监督分类，但决策树的构建过程复杂且该方法的普适性不强，对于深度较深区域识别精度不高，近几年，由于最大似然法的简便性，常用于结合航空摄影、航拍视频和数字测深模型进行分类，有效的提升了分类精度，但最大似然法最大的缺点是需要预先识别分类特征。 面向对象法对于地貌单元和底质分类存在很大的优势，以不同地貌单元的最优分割尺度、光谱参数、形状参数来分割影像并合并成不同地貌单元，分类过程比较繁琐，考验分类者的先验知识。支持向量机方法分类具有良好的分类性能，特别在小样本、非线性及高维特征空间具有较好的泛化能力，但在样本选择上也存在很大的局限性。 针对我国西沙群岛中赵述岛礁盘，本文利用珊瑚礁决策树分类模型和最大似然分类模型，并基于具有中分辨率的Landsat-8和高分辨率的Worldview-2影像，开展珊瑚礁信息识别和演变研究。 ' E, `; m+ V9 N/ y
研究区域与数据源 ⒈研究区概况 本文选取的研究区为赵述岛及其所属礁盘。赵述岛位于七连屿的西北部，是七连屿中的一座珊瑚岛礁(图1)，处于112°15′00″—112°18′00″E、16°57′00″—17°00′00″N。
8 L$ B, U R+ Q2 j8 M 图1 赵述岛位置 赵述岛礁体由浅水造礁珊瑚构成，是典型的珊瑚礁地貌体，具有最完好的珊瑚礁系统。礁盘内部水较深，边缘礁体水较浅，水质为一类水体，其礁盘主要包含礁前带、礁核带、礁后带3带。①礁前坡是礁前带的主要地貌类型之一，可见大量海滩岩；②礁核带主要以礁坪地貌为主，由珊瑚礁岩构成，在该区域鹿角珊瑚繁茂生长；③礁后带多由珊瑚砂组成，主要存在灰沙岛和泻湖坡地貌，泻湖坡有分枝状珊瑚生长。礁盘以外是暗滩，主要由砾石构成。 ⒉数据源 研究使用的遥感数据有两种，一种是Landsat-8多光谱数据，数据采集时间为2014年6月25日、2015年7月30日、2016年6月14日、2017年8月4日和2018年3月16日，另一种是WorldView-2多光谱数据，采集时间为2017年2月18日，本文使用的Landsat-8和WorldView-2的多光谱数据波段参数如表1所示。 表1 Landsat-8和WorldView-2的波段参数 # Q; F9 I4 o3 r
8 e$ @: T5 A) M
9 G, B& e; {2 Z. m; G) b数据处理方法 ⒈技术路线 本文利用决策树分类模型和最大似然分类模型进行珊瑚礁信息提取，其基本思想是：利用水深校正和最佳波段组合法提高遥感影像的分类精度，并结合决策树分类模型的优势提升最大似然法的影像总体分类精度，最后分析近5年赵述岛的地貌和浅海的底质变化。技术流程如图2所示。 / D: `" ^% `. f1 R* ?
图2 技术路线 ⒉遥感影像预处理 ⑴辐射校正 辐射定标可以消除由辐射误差引起的影像畸变，公式如下： Lλ＝(ρλ×Eλ×sinθ)/(πd2)，⑴ Lλ是波段λ的辐射亮度值，ρλ表示大气表观反射率，Eλ表示太阳辐射照度，θ表示太阳高度角，d表示日地距离。 图3为辐射校正前后的植被的像元光谱。图3a表示辐射定标前的光谱曲线，纵坐标表示的是DN值，DN值本身就是一些无量纲的数字。图3b表示辐射定标后的光谱曲线，纵坐标代表辐亮度，辐亮度值主要集中在0～10μW·cm–2·sr–1·nm–1之间，符合辐射定标的标准。
+ m) k( L- W9 b0 q0 ]7 [0 w+ K图3 Landsat-8辐射校正前后光谱曲线 ⑵大气校正 大气校正可以获得物体表面的真实光谱信息，本文采用FLAASH大气校正法，应用MODTRAN4+辐射传输模型进行大气校正。图4为大气校正前后Landsat-8遥感图像中赵述岛上植被像元的光谱，图4a表示Landsat-8大气校正前光谱特征曲线，纵坐标表示的是辐亮度。图4b表示Landsat-8大气校正后光谱特征曲线，其纵坐标表示的是不同波长的植被的反射率值，在0～1之间，与地表真实波谱曲线吻合，即表现为绿波段的反射、红波段的吸收和近红外波段的强反射验证了大气校正的正确性。 ! f8 u6 ?% f9 O+ f
图4 Landsat-8大气校正前后光谱特征曲线 ⒊决策树分类 本文采用决策树算法对赵述岛浅海底质进行监督分类，实现对植被、裸地、海水、珊瑚礁和云的识别。利用归一化差异水体指数(MNDWI)、归一化植被指数(NDVI)、近红外波段(Band5)、蓝光波段(Band2)、及热红外波段(Band10)，构建决策树分类模型(图5)。选用特征数据的主要原因是不同波长之间的穿透能力存在光谱差异。 # M( Y- ~6 z8 T J) n7 X% N
图5 决策树分类的技术路线图 如图6所示，用a、b、c、d和e，分别表示2014年、2015年、2016年、2017年和2018年决策树分类结果。从图中可以看出：①2015年和2014年的分类结果表明4类底质变化较小。②2016年与2015年相比，裸地面积变大，主要是由于人为活动的加剧和围海造田工程。③2017年和2016年相比，总体面积变化不大。④2018年相比于2017年，珊瑚礁面积减少，主要原因与气候变化和海平面上升、人为活动等有关。 ! z1 b: y: {0 y1 Z) M5 ^+ d9 s
图6 决策树分类结果 如图6所示，用a、b、c、d和e，分别表示2014年、2015年、2016年、2017年和2018年决策树分类结果。从图中可以看出：①2015年和2014年的分类结果表明4类底质变化较小。②2016年与2015年相比，裸地面积变大，主要是由于人为活动的加剧和围海造田工程。③2017年和2016年相比，总体面积变化不大。④2018年相比于2017年，珊瑚礁面积减少，主要原因与气候变化和海平面上升、人为活动等有关。 ⒋类型转移矩阵 转移矩阵可以反映某一区域某时间段始末各类型面积之间相互转化的动态过程信息。本文使用转移矩阵分析赵述岛不同地貌之间的面积转移变化，其数学形式为： 5 Q0 F7 G3 z; Y, j+ f% m; [+ g
式中，S表示面积；n代表土地利用的类型数；i、j分别代表研究期初和研究期末的土地利用类型。 利用类型转移矩阵，计算可得表2，由表可知2014年和2018年两年之间的类型转移面积变化： ①2014年海水的面积为1975.59hm2，2018年面积变为2065.93hm2；2014年裸地面积为20.3hm2，2018年面积变为31.3hm2；2014年珊瑚礁面积为187.07hm2，但到2018年面积下降为85.6hm2。②2014年到2018年，有35.89hm2海水转化为珊瑚礁，7.79hm2海水转化为裸地，0.12hm2裸地转化为海水。③结合表2和图7，可以清晰的看出，2018年与2014年相比，赵述岛西北部有35.89hm2海水转化为珊瑚礁，西南部有134.01hm2珊瑚礁转化为海水，珊瑚礁面积急剧下降。 表2 2014—2018年面积转移矩阵 3 Z- u/ }7 [+ {4 T
图7 2014和2018年面积转移图 ⒌水深校正 海水是影响珊瑚礁信息识别的一个重要因素，主要原因是由于光谱会随水深的变化而变化，在水下同样的物质可能具有不同的反射率，这对水下信息的探索造成了较大的干扰，因此需要进行水深校正。针对浅水中底质反射率与水体的衰减系数，Lyzenga提出了校正水深的半经验型方法。 水深校正算法是基于底部反射率和水深指数函数的近似线性函数，利用辐射值的自然对数函数，使其相对于深度的衰减效应线性化，从而建立变换辐射： Xi＝ln(L(λ)i)，(3) 式中，L(λ)i为影像大气校正后波段i的反射率。Xi是在波段i的变换辐射。图8为水深校正前后的影像，图8a表示水深校正前Landsat-8影像，图8b表示水深校正后Landsat-8影像。
) R# V1 e) `3 U7 m图8 水深校正前后Landsat-8影像 ⒍最佳波段组合 为得到精度更高的图像解译结果，选择最大信息量的波段组合较为关键。本次研究使用美国查维茨提出的最佳指数法(OIF)确定波段的组合，该方法计算简单，便于操作，可依据图像数据的标准差和波段间的方差计算得出，公式如下： Y: I. u8 Y, A
式中，Si为第i个波段的标准差；Rij为i、j这2个波段的相关系数。OIF指数越大，则组合波段所包含的信息量就越大，并且波段间的相关性越小，根据计算可得信息量最大的最佳波段组合为：YAYBYC。YA是3和4波段运算后产生的新波段；YB是2和4波段运算后产生的新波段；YC是1和2波段运算后产生的新波段。如图9所示，图9a表示的是YA波段组合影像，图9b表示的是YB波段组合影像，图9c表示的是YC波段组合影像，图9d表示YA波段、YB波段和YC波段重新组合后的新波段。
( U7 R* z! I3 o9 c6 l图9 波段运算 ⒎最大似然法分类 针对赵述岛的珊瑚礁地貌类型，建立相应的分类体系，分为灰沙岛、泻湖坡、礁坪、礁前坡和暗滩5类。分类体系表3所示，可以清晰的区分不同的物质类型。 表3 赵述岛分类样本构建 6 B1 r" ?% Q0 f2 a0 n2 q" O6 q& `
+ v8 E# Y% }+ e' z" U
& J+ a- @, u! q# U: @: E3 p
结果分析与精度评价 ⒈分类结果与分析 基于2017年的Landsat-8数据，分别利用最大似然法和水深校正—波段组合—最大似然法2种方法进行监督分类，可得图10，a、b、c分别表示基于最大似然法的地貌分类结果、基于水深校正—波段组合—最大似然法的地貌分类结果和基于水深校正—波段组合—最大似然法—决策树分类法的地貌分类结果。由图可知：b相比于a，泻湖坡区域面积减少，灰沙岛区域面积也减少，存在灰沙岛的分类精度不高的问题。为了提升b的分类精度，利用决策树分类法对灰沙岛区域进行纠正，可得c，将不同方法的分类结果进行比较(图11)，可以明显看出：基于水深校正—波段组合—最大似然法—决策树分类法的地貌分类效果更好。 ( w% k- N' |: \
图10 地貌分类 1 i+ X- E6 f: T8 P: k/ C) Y
图11 基于三种分类方法的地貌面积统计 ⒉精度评价 本文使用2017年2月18日拍摄的空间分辨率为2m的WorldView-2多光谱影像作为验证数据，对分类结果进行精度评价。 利用混淆矩阵对使用最大似然法监督分类的影像进行精度评价，检验精度见表4，结果表明：①基于最大似然法的赵述岛分类结果总体精度达到92.30%，Kappa系数达到0.905。②灰沙岛的漏分率最高，达到22.35%，制图精度最低，为77.65%；礁坪的制图精度最高，达到98.98%，错分率最高，达到15.11%；泻湖坡的错分率和漏分率都很高，分别为11.07%和20.64%，制图精度只有79.36%。③礁前坡、海面、暗滩三种物质的分类结果良好，漏分率低，制图精度高。 表4 Landsat-8影像的精度验证(最大似然法)
( {2 p- z; n- m将水深校正、波段组合及最大似然法监督分类后影像进行精度验证，检验精度见表5，结果表明：①影像分类后的总体精度由原来的92.30%提升到94.81%，Kappa系数也由原来的0.905提升到0.906。②经过水深校正和波段组合后的影像的错分率和漏分率降低，制图精度明显提升，只有灰沙岛和礁前坡的制图精度降低。 表5 Landsat-8影像的精度验证(水深校正—波段组合—最大似然法) " |* t, ^- A% m# B: Z/ [( }
利用决策树分类法结合水深校正—波段组合—最大似然法，提升灰沙岛的分类精度，对于珊瑚礁系统的信息识别可以得到较好的效果(表6)。结果表明：①灰沙岛的漏分率由50%下降到0，制图精度由50%提升到100%。②影像分类后的总体分类精度由原来的94.81提升到95.07，Kappa系数也由原来的0.906提升到0.939。 表6 Landsat-8最大似然法影像的精度验证(水深校正—波段组合—最大似然法—决策树分类)
# G- T" O0 P" h: Z/ G+ _% x( L
# W+ F! w. f" v& ?2 @2 O
' M" M5 ]: p% c" M. P" W( f* c珊瑚礁时空变化分析 ⒈赵述岛地貌类型时空变化 ⑴赵述岛地貌类型时间变化 基于2014年到2018年5年的Landsat-8影像，进行监督分类，可以获得历年赵述岛的珊瑚礁识别图(图12)，a、b、c、d、e分别表示2014年、2015年、2016年、2017年和2018年珊瑚礁识别结果，从图中可以看出赵述岛的5类地貌变化趋势较小。 对赵述岛的地貌类型面积进行统计(图13)，可以清晰的看出：①灰沙岛的面积呈现增长趋势，这与实际情况不符，由于人类活动的加剧，灰沙岛附近人工建筑在不断扩张，导致灰沙岛的面积持续增长。②2014到2018年间，泻湖坡的面积围绕200hm2面积横线上下浮动，但总体趋势是下降趋势。③5a间，礁坪面积稳定在630hm2，2018年相比于2014年，增长了23.6hm2，总体呈增长趋势。④礁前坡的面积变化呈现下降趋势，在2017年出现最低点，2014年泻湖坡的面积大于2018年其面积，说明泻湖坡面积在退化。⑤海面的面积变化主要取决于其他五类物质的面积变化，从图中还可以看出与暗滩的变化存在互补趋势。⑥暗滩面积波动变化幅度小，面积最高点为2014年的402.56hm2，最低点为2015年的328.81hm2，2016年及其以后，面积略有增加，直到2018年，面积达到380.56hm2，但相比于2014年，暗滩面积在减少。2015年出现最低点，与海水对砾石的冲积和堆积有关，同时也与水深和潮汐有关，当浅海水深上升时，暗滩的识别会更加困难。 3 Q7 U- F* u! \) @
图12 历年珊瑚礁识别结果 ' l9 B( Y: V+ Y
图13 历年底质变化图 ⑵赵述岛地貌类型空间变化 根据2014年和2018年的珊瑚礁分类结果，计算面积转移矩阵，如表7所示，从表中可以看出： ①345.17hm2的暗滩面积未发生变化，但44.55hm2的暗滩转化为礁前坡；海面有361.7hm2面积未转移，16.32hm2海面转移为礁前坡；仅有的20.19hm2的灰沙岛到2018年依然是灰沙岛；36.92hm2礁坪转化为泻湖坡，但仍有556.42hm2礁坪未发生改变；礁前坡转化为暗滩和海面的面积分别为32.41hm2和63.66hm2，400.4hm2礁前坡面积依然是礁前坡；泻湖坡有70.68hm2转移为礁坪，130.9hm2未发生变化。 ②从总体来看，暗滩、礁前坡和泻湖坡面积减少，海面、灰沙岛和礁坪的面积增多；暗滩、礁前坡和泻湖坡面积减少，东北季风、热带气旋温度上升的综合影响是驱动因素；海面面积增多的主要原因与全球气候变暖导致的海平面上升有关；礁坪主要由白化后的珊瑚砂石构成，礁坪的面积增加表明珊瑚白化趋势明显，主要与航道和码头的人为开发活动相关；近年来赵述岛增设人工建筑，导致了灰沙岛面积的增加。 将面积转化矩阵以空间面积变化的形式表现，如图14所示，从图中可以明显的看出赵述岛的面积转移情况，礁坪和礁前坡和泻湖坡的面积转移明显，但仍然符合聚类分布规律。将具体的转移面积进行统计，如表8所示，未发生面积转移的达到1814.42hm2。
0 `" U. G" u4 F; Q0 v" c图14 2014年和2018年地貌变化图 表8 2014—2018地貌面积变化(精确到小数点后两位) ⒉赵述岛珊瑚礁时空演变 结合决策树模型对珊瑚礁的底质分类和最大似然法对珊瑚礁的地貌分类，可得基于地貌分类的珊瑚礁变化(图15)，a和b分别表示2014年和2015年基于地貌分类的珊瑚礁分类，结果表明：①2014年，珊瑚礁主要分布在泻湖坡地貌区，2018年，珊瑚礁的分布范围在扩散，泻湖坡和礁坪都有珊瑚礁分布。②从总体上来讲，珊瑚礁的面积在减少，直观表现出这5年间，我国南海赵述岛地区珊瑚礁退化严重，为我国珊瑚礁保护敲响了警钟。 2 g2 ?7 U- f4 V5 n7 _" l! O
图15 基于地貌分类的珊瑚礁分类
( d8 K* v- ?" u' x2 V6 n9 w6 ], | a5 m2 V1 l$ _9 n
结论 本文利用2014—2018年的Landsat-8影像和2017年的worldview-2影像，以西沙群岛海域赵述岛为研究对象，建立了珊瑚礁决策树底质分类模型，取得的结论如下： ⑴利用最大似然法的监督分类方法对影像地貌进行监督分类，总体分类精度为92.30%，Kappa系数为0.905。⑵使用水深校正算法和波段组合法对影像进行处理能够明显提高浅海地貌分类的识别精度，总体精度由原来的92.30%提升到94.81%，再结合决策树分类结果对灰沙岛分类结果进行改正，总体精度由94.81%提升到95.07%，对于珊瑚礁系统的信息识别可以得到较好的效果。⑶通过对2014和2018分类后影像进行转移矩阵计算，得到近5年赵述岛礁盘珊瑚礁地貌类型的变化情况，暗滩、礁前坡和泻湖坡面积减少，海面、灰沙岛和礁坪的面积增多，该结果同时表明基于多时相遥感影像的变化检测方法的可行性。⑷通过地貌分类基础上的珊瑚礁演变分析，表明珊瑚礁的扩散和退化现象明显。 | 
