-
湿地被誉为“地球之肾”,为众多生物提供生存环境,在生态环境发展、地下水补给、防洪防灾、侵蚀控制、保护生物多样性等方面都发挥着巨大作用,有蓄水防洪的天然“海绵”之称,是重要的生态系统。滨海湿地的水体信息提取对于海岸线变化及侵蚀监测,海岸带管理,洪水预测,水资源评估等有重要的意义[1-2]。随着中国城市化程度的加深,人们居住环境的改善,农业及基础设施的发展,使得许多天然海滨湿地变为建设用地和人工湿地,生境破碎化严重[3]。到20世纪90年代中期,50.00%的滨海沼泽已经消失,近1 000个天然湖泊消亡。因此,及时、精确地掌握湿地的面积变化等是进行湿地保护的基本条件。近年来,卫星遥感技术在空间和时间的分辨率已达到较高水平,可在水体监测方面发挥重要作用。将遥感与地理信息系统技术相结合,计算各种水体指数[4-5]、机器学习算法分类[6-9]来开展水体自动或半自动的提取研究[10-11]。由于湿地环境复杂,位于水陆交汇处,水体类型多样,异质化程度高,“同谱异物”等现象影响了湿地信息的提取精度[12]。相对简单快速的水体提取方法主要是单波段阈值法和谱间分析法[4, 13-14]。单波段阈值法是利用水体在近红外波段上反射率低,易与其他地物区分的特点,选取一个阈值来分离水体与其他地物。但是使用该方法提取时,会同时提取山体阴影区域,并且阈值会随时间、环境而变化,不易确定。谱间分析法是通过分析水体与其他地物之间的光谱特征,寻找其中的规律,建立数学模型来分离水体与其他地物,效果优于单波段阈值法,应用较为广泛,但受使用条件和精度范围等限制。将严谨的数学理论、多数据源等相结合,可提高水体提取精度[2, 15-16],如PANG等[17]基于小波变换的图像纹理特征的水体遥感信息提取,在边缘信息提取和分割上取得了很好的效果。虽然上述方法在某些特征水域提取方面具有优势,但实验发现,对于湿地中大量存在的零星破碎水域的提取还不能达到满意的精度,直接影响水体提取的总体精度。鉴于此,本研究构建了新的水体指数模型,以滨海湿地水域为研究对象,开展水域特别是破碎度高的水域信息的提取研究,并验证了构建模型的水域提取有效性。
HTML
-
徐涵秋[5]观察分析影像的光谱曲线发现:建筑物的反射光谱范围为近红外波段到中红外波段(MIR),将McFEETERS[4]的归一化水指数(NDWI)做了修改,用中红外波段替换近红外波段,得到的建筑物的指数值明显减少,从而增大了水体和建筑物的差别,提高了水体提取的精确度。改进的归一化差异水体指数的公式为:
式(1)中:IMNDWI为改进的归一化差异水体指数;G为绿波段;MIR为中红外波段。
MNDWI能够较好地区分建筑物与水体,但是对于水体分布破碎的密集人工围垦区还不适用。
-
凌成星等[14]通过分析NDWI和MNDWI的优缺点和各波段光谱特征的基础上,为了增强水体、植被和建筑的差异,将归一化差分植被指数和近红外波段进行求和,其公式为:
式(2)中:ICIWI为混合水体指数;INDVI为归一化植被指数,其公式为:INDVI=(NIR-R)/(NIR+R),R代表红波段,NIR为近红外波段。
-
上述指数多基于单个像元的数学运算,对于异质性高的破碎水域,像元的数值变化跨度较大,很难通过单一阈值确定破碎水域范围,因此该类型水域提取精度较低。本研究借鉴HUANG等[18]提出的联合植被指数(IFZ)对森林区域快速提取的方法,利用Landsat TM遥感数据,改进并构建新型水体指数WZ5(water zscore-band5)拟提高对破碎水域水体提取的精度。在影像中水体训练样本选取时,不是仅局限于单个像元的样本提取,而是选取典型面状水体或破碎水域中多个像元,统计其均值和标准差,对影像基于水体特征进行标准化,构建水体指数WZ5。湿地水域分布广,但多为浅水区域,可见光、近红外波段特征与其他类型区分度低。中红外B5波段(1.55~1.75 μm)位于2个吸收带之间,对土壤和植被的水分含量敏感。经过重复实验发现,B5波段水体区分度最好,相对于其他各波段,中红外波段(B5)的水体与塘基的像元值差异最大。所以研究选用B5构建水体指数WZ5。由于样本均值反映水体像元辐射值的平均水平,标准化后会产生水体像元WZ5值小于0的情况,为保持一致性,利用样本均值的10%进行标准化。
Table 1. Interpret sign of water in Hangzhou Bay Wetland
式(3)中:IWZ5为新型水体指数;PB5表示中红外波段第5波段B5各像元反射率;M表示水体样本的平均值;SDB5表示第5波段B5标准差。
-
为了验证指数WZ5的有效性,对研究区分别计算MNDWI,CIWI和WZ5指数值,进行水体信息提取实验。由于不同指数计算方法不同,导致量纲不同,不适合采用方差等统计量表达差异;另外,此处需要区分不同指数对水体与其他类型的区分程度,即类间差异比较。因此,本研究利用表达各类均值差异的区分度指标对各指数的提取精度进行评价,其计算公式如下:
式(4)中:DD为区分度,反映2种类型之间的差异程度,其值越大,表示两者的可分离程度越高,反之,则相反[16]; ${{{\bar W}_{{\rm{Ix}}}}}$ , ${{{\bar W}_{{\rm{Iy}}}}}$ 分别表示待区分的2种类型的均值。
2.1. 水体提取方法
2.1.1. 改进的归一化差异水体指数(MNDWI)
2.1.2. 混合水体指数(CIWI)
2.1.3. 构建WZ5水体指数
2.2. 水体提取结果的精度评价
-
开阔水域具有面积大,边界清晰,背景单一等特点。分别计算指数MNDWI,CIWI和WZ5,对水体增强并用阈值法提取水体信息(图 2)。
比较分析提取结果可以看出:对于开阔水域,3个指数的提取结果都较好,但WZ5和CIWI提取的结果边界更加清晰(图 2)。破碎水域水体分布零散,面积较小,边界不易区分,如沼泽中的水体,排水渠,混有塘基的围垦区等分布复杂的背景,极大影响了水体提取的精度。沼泽区,各指数的提取结果没有明显差异;然而对于排水渠,WZ5和CIWI能够完整的分离水体与人工筑坝,提取出不规则水域,MNDWI只能稀疏显示轮廓;对于水产池塘内存在的塘基,WZ5能够很好地剔除,而其他指数明显不能区分。
-
图 3可以看出:对于破碎区域水体,能否有效分离水域及人工建筑(如围垦区的塘基)是关键。因此,研究提取4种主要的地物类型(塘基、水体、植被、建筑)各20个样本,计算MNDWI,CIWI和WZ5,并将值域拉伸至相同范围,比较指数值差异(图 3)。
对于提取结果差异较大的排水渠和水产池塘区域,结合图 3各类指数值可以看出:WZ5各要素区分度最明显,而MNDWI和CIWI的指数值都有不同程度的重合。其中,MNDWI对塘基、水体和建筑的区分度最差,这导致了对排水渠、水产池塘识别的错误率高,特别受塘基的影响较大,不能很好地从大量塘基中提取出水体,提取结果像元分布杂乱,很多像元没有很好地识别;CIWI的植被和建筑区分度较差,而水体和塘基基本无重叠区域,但指数值差异较小,可以将排水渠中的水体、人工堤坝、周围植被很好地分离,塘基与水体虽然也可以区分,但是分界线不是很明显,有些非水体部分也同时被提取,不能辨别水体间较小的塘基;4个地物类型间WZ5数值差异较大,图 3的4种水体类型提取效果均优于另外2个指数,特别是排水渠边界提取光滑,水产池塘水体信息提取完整,分布均匀有序,与实际最为符合,较好地剔除了水体间较小的塘基。
-
以破碎水域区的水产池塘为例,随机选取了25个水体样点和附近的25个塘基样点,统计各像素的指数平均值,计算3种指数破碎水域的区分度(表 2)。
项目 IMNDWI ICIWI IWZ5 水体 塘基 水体 塘基 水体 塘基 最小值 38.48 28.67 2.98 7.15 4.31 10.98 最大值 78.00 38.46 7.13 15.93 9.02 18.43 均值 54.13 34.91 5.21 10.71 6.21 15.12 区分度/% 21.59 34.56 41.78 Table 2. Division of different indexes
从表 2可以看出:WZ5的区分度高达41.78%,大于MNDWI(21.59%)和CIWI(34.56%)的区分度,表明通过WZ5增强后能够有效提高破碎水域水体的提取精度。
-
利用外业采样及Google Earth等高空间分辨率影像转绘部分验证样点进行精度评价,在水体和非水体分别选取了100个样点作为验证信息,其中WZ5湿地水体与非水体的总体精度达82.00%,高于58.00%的MNDWI和77.00%的CIWI,WZ5的Kappa系数为0.64,也高于指数MNDWI(0.16)和CIWI(0.54)。具体分类精度见表 3。
指数 MNDWI精度 CIWI精度 WZ5精度 水体 非水体 水体 非水体 水体 非水体 水体 61 45 91 27 74 26 非水体 39 55 8 73 20 80 总体精度/% 58.00 77.00 82.00 Kappa系数 0.16 0.54 0.64 Table 3. Classification evaluation of different indexes
从表 3得出:在破碎度很大的水体、塘基混合区域,WZ5的识别精度最高,而指数MNDWI和CIWI受塘基的影响较大,提取精度不是很高。为直观比较本研究构建水体指数WZ5在破碎水域提取的准确度高于MNDWI和CIWI,在研究区选取塘基、水体典型区域各若干像元,制作二维散点分布图,并设置各指数提取的阈值,比较3种指数提取结果的差异(图 4)。
结果表明:塘基的MNDWI和CIWI范围跨度较大,导致塘基识别率较低,难以从水体中分离出来,而WZ5值相对较集中,塘基提取较为完整;对破碎水域的提取,阈值设置对结果影响较大,WZ5能够提取完整,而MNDWI和CIWI却不能提取完整。