武汉市湖泊景观动态遥感分析(1973 -2013年)*

陈海珍，石铁柱，邬国锋

(1:武汉大学资源与环境科学学院，武汉430079)

(2:深圳大学海岸带地理环境监测国家测绘地理信息局重点实验室＆深圳大学空间信息智能感知与服务深圳市重点实验室＆深圳大学生命科学学院，深圳518060)

城市湖泊是重要的城市水体形态，在抵御洪水、调节径流、改善气候、维持生态平衡等方面都承担着重要的功能.湖泊是武汉市重要的城市名片之一［1］.武汉市地势低洼，雨量充沛，且季节分配不均匀;但湖区自然地理条件决定了武汉市湖泊多为浅水型湖泊［2］，生态环境较脆弱，极易受到外界的干扰与影响.自1973年以来，由于武汉城市圈经济飞速发展、人口显著增长、湖泊的不合理开发和过度利用以及湖泊的自然演化等因素，武汉湖泊景观结构发生了巨大的变化，生态环境遭受严重破坏.这些变化和破坏影响了湖泊调蓄功能，进而影响到武汉周边地区小气候的变化，带来诸多环境和生态问题［3-4］.研究武汉湖泊长期的时空动态变化，对湖泊生态环境改善和湖区经济可持续发展意义重大.

遥感卫星影像为研究长时间序列的湖泊时空变化提供了良好的数据源，地理信息系统(GIS)作为重要工具广泛应用于湖泊和湿地的研究中［5］.武汉市湖泊属于江汉湖区的重要组成部分，国内外对江汉湖区的动态变化研究较多，主要是应用遥感和GIS技术，结合历史数据和土地利用图，研究湖泊数量和面积的长期变化趋势以及政策和人为干扰的影响［6-8］.但国内外对武汉市湖泊的研究还较少，主要是对典型湖泊，如东湖的富营养化、大型植物以及金属污染的研究［9-11］，武汉市湿地景观格局及其演化相关研究［12-13］等.但总体而言，对武汉市湖泊的时空动态变化，尤其是从景观角度出发、比较系统的研究仍较少.

景观格局是指景观的空间格局，是大小、形状、属性不一的斑块在空间上的分布与组合规律［14］.景观格局分析广泛应用于森林生态系统、土地利用演变等方面的研究［15-16］;近几十年来，从景观角度分析溪流和河流水系的演变也表现出极大潜力［17-21］.湖区总体演变体现在湖泊大小、形状和湖群空间格局变化，所以从景观角度出发分析湖泊动态变化也具有很大优势.景观指数是指能高度浓缩景观格局信息，反映其结构组成和空间配置某些方面特征的定量指标［14，22］，是研究景观格局变化的重要工具.1950年以来，国内外学者提出众多景观指数［23］.从不同水平和尺度角度，这些景观指数可分为3大类:斑块水平景观指数、斑块类型水平指数以及景观水平指数［14］.根据其计算公式又可划分为:信息论类型、面积周长比类型、简单统计学指标类型、空间相邻或自相关类型以及分维类型等［14］.由于景观指数间的相关性很高，故选取适宜的景观指数，避免信息冗余是利用景观指数研究景观动态变化的第一步.

结合遥感和GIS技术，从景观角度出发利用景观指数定量分析研究大范围湖区的时空动态变化具有较大优势.本文基于9景遥感影像进行湖泊提取，选取景观指数，旨在从景观角度分析武汉市主要湖泊1973-2013年近40年的动态变化.

1 研究区概况

武汉市(29°58'～31°22'N，113°41'～115°05'E)位于湖北省东部、长江和汉水的交汇处，是湖北省政治、经济、文化中心，有“九省通衢”的美称.武汉地形属于残丘性河湖冲积平原，地势平坦低洼，以平原为主，中部散列东西向残丘;气候属北亚热带季风湿润气候，冬冷夏热，常年雨量充沛，雨热同期，年降水量为1050～1200 mm，降水多集中在6 -8月，约占全年降水量的40%［24］.武汉市共有13个市辖区，包括7个主城区(江汉、江岸、硚口、汉阳、武昌、青山、洪山)以及6个远城区(蔡甸、江夏、黄陂、新洲、东西湖、汉南).由于影像覆盖区域限制，本研究不考虑黄陂和新洲这2个远城区.武汉市的行政区划范围在40年内没有较大变动.主城区中，洪山区在1987年建制确定，郊区1980年后行政区划范围没有变动.本文以1987年以后的武汉市行政区划范围为准.考虑到武汉市行政边界上的大型湖泊(梁子湖、斧头湖)的完整性，研究区扩展到鄂州梁子湖区，咸宁市西北部分地区和嘉鱼市的东部斧头湖区.

2 数据与方法

2.1 Landsat影像

本研究使用1973 -2013年8 -11月9景覆盖研究区的Landsat影像，平均年际间隔约为4年，影像具体信息如表1所示.

基于ENVI软件，对影像进行了如下预处理:(1)采用多项式法进行相对几何精校正，误差控制在1个象元内;(2)采用双线性内插法进行重采样，将空间分辨率为79 m和30 m的MSS和TM，OLI影像分别重采样为50 m;(3)以2013年影像为标准，采用直方图匹配方法，调节其余影像的色调与反差，以降低色调与色差对目视解译的干扰;(4)根据研究区进行掩膜裁剪，获得研究区的影像.

2.2 水体提取和目视解译

遥感影像分类精度势必会对后续的景观指数计算有重大影响［25］.一般来说，简单的分类体系(更少的类别)会有更高的分类精度［25］，如水陆的二元分类，甚至可以达到95%的分类精度［26］.故本文也仅采用水陆二元分类进行湖泊提取.

水体具有独特的光谱反射特征，在近红外及更长波段的吸收率高，因此Landsat MSS、Landsat 8 OLI和Landsat 5 TM的近红外波段可以用于水体的提取;但水体组分的不同会导致水体反射率的变化，同时，浅水水体底部物质也会影响其反射率［27］.由于受到人为和周围环境的影响，城市浅水湖泊的水体成分复杂，仅使用波段阈值法或(非)监督分类的误差较大.故本文首先采用近红外波段设定阈值得到水陆二元影像［26］，进而结合目视解译来提取湖泊水体.由于影像精度的限制以及目视解译的精度限制，本研究中仅提取＞0.1 km2(采样后的40个像素大小)的湖泊.具体方法为:首先利用ENVI软件选取近红外波段设定直方图阈值，进行水体提取，然后以二值化影像为参考影像，利用ArcGIS对影像进行目视解译，尽量消除因水体富营养化等原因而错漏的湖泊，交互提取湖泊矢量图层.

2.3 景观指数选取

本研究中仅分为湖泊和非湖泊两类，所以应选取斑块水平景观指数;同时，景观指数的选取还需考虑指数与景观的生态过程的联系，所以应选取适宜描述湖泊动态变化的指数.经综合考虑，斑块数量、边界长度、斑块面积、面积周长比、斑块形状指数、分维数等方面共15个景观结构指数可用于湖泊景观动态变化的描述(表2).

本文使用Matlab软件进行指标计算.虽然景观指数数目繁多，但同类型指数之间的相关性往往很高，因此采用更多的指数不能表示更多的信息.当变量较多且变量相关性较大的情况下，常常需要通过选取几个不相关的主成分反映原变量的信息量以达到降维目的［22］.主成分分析法是常用的一种简化数据结构、选择变量子集分析方法.本文利用SPSS统计软件进行主成分分析，计算主成分的矩阵负载.矩阵负载可用来表征原变量属性对主成分的负荷量也可叫做权重，反映各变量对各主成分的重要性.

3 结果

3.1 湖泊提取

从1973 -2013年的9幅影像解译结果(图1)可以观察到，武汉湖泊的面积和数量都有明显变化.1973年湖泊总面积为983.29 km2，到2013年仅为647.47 km2，减少近35%.湖泊面积变化经历了3个阶段:1973 -1979年间湖泊面积急剧减小，减小超过300 km2，萎缩近30%;1980s以后，湖泊面积有所回升，到1996年达到755.34 km2;此后又以较低速率有波动地减少(图2).

斑块是构成景观的基本单元，斑块数量的变化也在很大程度上反映景观格局的变化.40年间，研究区中湖泊斑块数量也有较大波动(图3)，经历了3个阶段:(1)1970s至1980s中期(1973 -1984年间)，湖泊斑块数量急剧增加，从1973年的69个增加到1984年的109个;(2)1990s后，湖泊斑块数量持续减少，到1996年剩余89个;(3)此后以每年约1个湖泊斑块数量的速率缓慢增加，2013年湖泊斑块数量又达到109个.在湖泊斑块数量变化的同时，湖泊大小分布也有明显的变化(图3):前30年间，＜1 km2的湖泊斑块数量明显增加，从1973年的35个增加到1990年的61个，此后先减小后缓慢增加;在1970s初期＞50 km2的大型湖泊有4个，不到10年后仅余2个，1996年和2003年大型湖泊恢复到3个，但2006年后又仅余2个;1～10 km2和10～50 km2的湖泊斑块数量都经历了先增加后减少的阶段.

表2 文中使用的15个景观指数及其公式［28］*Tab.2 15 landscape matrices and the formula used in this paper

结合湖泊总面积和湖泊斑块数量变化，计算湖泊总面积变化趋势与湖泊斑块数量变化趋势的相关性，得到Pearson相关性指数为－0.73，在0.05水平上呈显著负相关.

图1 基于Landsat影像的1973 -2013年武汉湖泊提取结果Fig.1 Lakes of Wuhan from 1973 to 2013 extracted from Landsat images

图2 1973 -2013年武汉市湖泊总面积变化Fig.2 The change of total area of lakes in Wuhan City from 1973 to 2013

图3 1973 -2013年武汉市湖泊斑块数量分布Fig.3 Distribution of the number of lakes patches in Wuhan City from 1973 to 2013

表3 主成分分析结果Tab.3 The result of principal component analysis

3.2 景观空间结构指标

用主成分分析法对相关矩阵进行分析，前3个主成分的累积方差可以达到93.072%，达到表达85%以上景观信息的要求.第一主成分反映斑块面积，主要由PSSD和MPS所决定;第二主成分反映斑块面积变化和边界密度，主要由PSCV和ED所决定，结合指标公式可以得知PSCV直接可以由PSSD与MPS算出，故不选用PSCV，故第二主成分主要反映边界信息;第三主成分反映斑块分维数，主要由MPFD所决定(表3).

本研究中选取具有最大相关性的4个景观结构指数来表达3个主成分，分别是平均斑块面积、斑块面积标准差、边界密度和平均斑块分维数.

3.3 景观动态变化

3.3.1 MPS和PSSD动态变化 从 MPS变化来看，武汉市湖泊1970s的变化最为剧烈，MPS在6年间减少近一半;减少趋势直到1980s中期才停止，1990s后有缓慢回升;此后又以较低速率降低.PSSD变化趋势与MPS变化趋势基本一致，湖泊斑块面积大小趋同(图4).

3.3.2 ED变化 从ED变化来看，武汉市湖泊1970s变化最为剧烈，由1973年的近3300 m/km2减少至1979年的约2100 m/km2;1980s至1990s有缓慢上升趋势，至1990s初回升至近2400 m/km2，而后又以较低速率降低，降至2003年的2200 m/km2;之后有波动地缓慢上升.总体而言，1970s骤减近1/3，此后变化幅度明显减小，在2100～2500 m/km2间波动变化(图4).

3.3.3 MPFD变化 1970s至21世纪初MPFD呈下降趋势，具体可分为4个阶段:1970s和1990s减小速率较快，1980s缓慢减小，21世纪后以较快速率增加.总体而言，分维数下降趋势剧烈，湖泊边界变得更为简单和规则，进入21世纪后情况有所好转(图4).

图4 1973 -2013年武汉市湖泊景观平均斑块面积、斑块面积标准差、边界密度和平均斑块分维数变化Fig.4 Changes of mean patch size，standard deviation of patch size，edge density，mean patch fractal dimension number of lake landscape in Wuhan City from 1973 to 2013

4 讨论

4.1 分类精度讨论

遥感影像的分类精度评估对于研究的可信度意义重大.分类精度评估主要采用的方法有3种［29］:结合实地调查、高分辨率影像检验低分辨率影像的分类和以目视解译为真实信息对随机生成点经过目视解译.然而，遥感影像，尤其是历史时间序列的影像，获取研究区实时实地信息极为困难;获取对应时间的高分辨率影像也存在很大困难;通过随机点目视解译判读的方法可以计算相应的误差矩阵，进行精度分析，但这是假设目视解译的正确率为100%，事实上并非如此，而是受解译者的经验背景和判读水平等因素影响.所以时间序列遥感影像的分类精度评估只能退而求其次，用置信水平来评判分类精度［30-33］.水体的光谱特性较为明显，在红外波段吸收率很高，与其他地物光谱特性具有很大区分度.因此，近红外波段广泛用于水体分类和提取［34］，且分类精度较高，普遍高于 85%［26，35-39］，已经达到分类精度要求［40］.此外，本研究在近红外波段直方图阈值法的基础上，结合了目视解译，理论上精度会有相应提高，故有理由相信本研究中湖泊提取的精度达到分析使用的要求.

4.2 景观格局动态变化及原因分析

结合湖泊面积、湖泊斑块数量以及4个景观指数的动态变化可以看出，1973 -2013年这40年间武汉市湖泊变化经历了4个大的阶段:

1)197 0s期间，湖泊总面积、MPS和PSSD急剧减小，湖泊斑块数量(尤其是小型湖泊)剧烈增加，ED和MPFD剧烈减小.湖泊总面积与湖泊斑块数量呈显著负相关(r=－0.73，P＜0.05)，可以推断湖泊的变化主要体现在大湖破碎、小湖数量上升.探究其原因，主要是由于围湖造田.围湖垦殖、垸田发展自古有之，但程度较轻，且随大洪水的发生，农田被淹，垸田复湖.自新中国成立至1970s，在“以粮为纲”政策的指导下，掀起围湖垦殖高潮［8］，湖泊面积和边界都发生巨大变化，湖泊在衰退和消亡的过程中边界逐渐规则，形状日趋简单.

2)198 0s期间，政策导向发生改变，围湖垦殖的现象基本杜绝，但大湖破碎、小湖数量上升的趋势仍然持续到1980s中期.究其原因，主要是经过前10年的巨变，湖泊围垦和大湖破碎大大降低湖泊的自我调节能力，政策改变后，湖泊萎缩和破碎的趋势尚不能迅速回转.1980s中期后，退田还湖政策有较明显效果，湖泊总面积有了较大回升.

3)199 0s期间，中小型湖泊数量减少，ED和MPFD下降幅度较大，1996年湖泊总面积大幅回升;主要是因为当年大涝，浅水型湖泊的水面面积很大程度上受到降水量的影响.但是从湖泊数量、ED和MPFD的变化上可以推断出该时段的湖泊变化主要体现在中型湖泊破碎或萎缩、小型湖泊消亡，湖泊边界更为规整，对湖泊开发的力度仍然很大.这主要是因为随着武汉城市经济发展，人口继续增长，城市建设用地紧张;同时，渔业发展需要较大力度地湖泊开发利用.

4)进入21世纪后，湖泊斑块数量持续缓慢增加，但湖泊总面积呈减少趋势，故又有可能进入一轮大湖萎缩或破碎、小湖增加的阶段.围湖垦殖、填湖造房，加之天然或环境破坏造成的湖泊淤积，使得湖泊水面面积逐渐减小.即使在“退田还湖”政策实施后，多数围垦的湖泊并非退成天然湖泊，而是改造成精养鱼池、藕田或库塘.

4.3 选取指标的合理性分析

迄今为止，上百种景观指数可供使用，大多数指数都有很高的相关性，这带来一个问题就是使用景观指数进行定量分析前需要进行景观指数的选取.目前选取的方法主要有2种［25］:(1)根据经验或具体的研究对象进行选取［41］;(2)通过主成分分析、因子分析或其他多变量化简方法进行选取［22，42-43］.经验选取具有较大主观性，故本文结合经验选取和主成分分析法进行选取.即，首先选择可以用作衡量变化的15个景观指数，然后对这些指数进行主成分分析，最终选取出4个指数:MPS、PSSD、ED和MPFD.根据其计算公式又可以将4个指标可以归类为3种:面积相关、周长相关以及分维数.

湖泊水面面积是湖泊功能实现的最基本要素.湖泊总面积、MPS、PSSD的变化可以较为全面地描述湖泊面积的变化以及湖群中大小湖泊面积分布的变化.斑块数量和ED可以侧面体现湖泊的破碎化程度.天然湖泊的几何形状偏离欧几里得几何形状，分维数较大，若由于人为因素的影响，边界逐渐规则，则分维数会降低.武汉湖泊多为浅水型湖泊，极易由于淤积和人为围垦而改变湖泊水域的边界形状，分维数的变化速率也可以间接体现出对湖泊开发利用的力度.总体而言，这些参数结合斑块数量和湖泊总面积，可以较好地反映湖泊景观动态变化.

4.4 武汉市湖泊与江汉湖群变化的对比分析

武汉市湖泊是江汉湖群的重要组成部分.江汉湖群湖泊研究表明，江汉湖群湖泊与武汉市湖泊变化都是按照大湖破碎成小湖，小湖逐渐走向人工或自然衰退而消亡的轨迹演化，最后消失［44-45］.利用遥感和GIS技术手段，也有学者对江汉湖群湖泊的演化定量研究［8，44，46］，发现从数量上分析，江汉湖群湖泊数量于1970s迅速减少，1980s上升，1990s先呈下降之势后又略有上升.从湖泊面积上分析，1950s--1970s湖泊面积急剧减小，1980s后湖泊面积有波动地略微增加.对比武汉市湖泊的数量和面积变化以及演化趋势可以发现，武汉市湖泊的变化体现了江汉湖群变化的大体趋势，同时由于城市湖泊受到城市发展和人类活动的影响更大且更为直接，武汉市湖泊变化又与江汉湖群变化不完全同步.1950s--1970s，江汉湖群湖泊由于围湖造田，小型湖泊大量消亡，湖泊数量急剧减少，而1970s武汉市湖泊处于大湖破碎的阶段，可见小湖遭到围垦消亡发生时间较之江汉湖群湖泊更早，武汉人口密集发展快，人为干扰更为显著，围垦的进度也更快［47］.1980s开始政策导向改变，江汉湖群的变化和武汉湖泊变化较为同步，可见围垦力度和进度虽然更大，湖泊衰退和消亡现象更为严重，但是武汉作为城市，对政策的贯彻和执行力度更大，时间效率也更高［48］.

4.5 遥感在湖泊景观研究的优势和发展前景

景观空间格局的自然演变和社会人为干扰演变导致的空间格局变化长久以来都受到生态学家、湖泊学家、规划师和管理者的关注［14］.而遥感和GIS技术在很大程度上改变和主导了生态学家对景观生态的研究方式，并将持续有着重要的影响［5，25］.空间格局分析始于1970s，在几十年内发展迅速，从最初简单的手动地图叠加分析到现在结合GIS和影像处理软件可以在几分钟内自动计算上百种景观格局指数.遥感更是将研究者的视野从局部小区域扩展到了大范围、大尺度的景观格局研究［22-23，28］.Landsat等卫星影像广泛使用，已成为景观研究、环境科学研究的重要组成部分，同时高分辨率遥感影像(如IKONOS、Quick-Bird、SPOT等)使大范围、大尺度景观格局分析也可以达到很高的精度［49-50］.高分辨率遥感影像或航空影像的应用对湖泊景观研究势必有重大推进作用.

5 结论

从武汉市湖泊近40年的动态变化可以看出，城市内湖泊，尤其浅水型湖泊，其演变受政策导向和经济驱动影响剧烈.武汉市湖泊在40年中经历巨大的变化:1970s主要受”围湖造田”政策的影响，湖泊衰退和消亡现象严重;1980s政策导向改变，实行‘退田还湖’，湖泊开发有所缓解和改善;1990s后主要受经济驱动影响，湖泊萎缩和消亡有所缓解，但趋势仍未改变.湖群数量和面积的减少，会严重地破坏湖区生态平衡，使水分涵容与调蓄能力大为降低，春旱夏涝频度增加，局部地区灾害性天气加剧，影响农业增产和人民的生命财产安全.必须采取有效措施，加强湖泊的科学管理，注重综合经济效益与生态效益，合理地开发利用，保护好湖泊自然资源.

［1]夏开元.湖泊中的武汉.城建档案，2007，9:16-19.

［2]张 欢.近十年武汉市中心城区湖泊生态环境演变［学位论文］.武汉:华中师范大学，2012.

［3]陈宝玉.对武汉市水资源开发利用的思考.人民长江，2008，39(12):38-39.

［4]黄利民，刘成武.湖北省近代洪涝灾害的特点及其成因分析.华中师范大学学报:自然科学版，2006，40(1):115-118.

［5]Griffith JA.Geographic techniques and recent applications of remote sensing to landscape-water quality studies.Water，Air，and Soil Pollution，2002，138(1/2/3/4):181-197.

［6]Fang J，Rao S，Zhao S.Human-induced long-term changes in the lakes of the Jianghan Plain，Central Yangtze.Frontiers in Ecology and the Environment，2005，3(4):186-192.

［7]Du Y，Xue H，Wu S et al.Lake area changes in the middle Yangtze region of China over the 20th century.Journal of Environmental Management，2011，92(4):1248-1255.

［8]邓宏兵.江汉湖群演化与湖区可持续发展研究［学位论文］.上海:华东师范大学，2004.

［9]Gan Y，Guo Y.Evaluation analysis and remedy strategy for eutrophication in Wuhan Lake Donghu.Resources and Environment in the Yangtze Valley，2003，13(3):277-281.

［10]Gao J，Xiong Z，Zhang J et al.Phosphorus removal from water of eutrophic Lake Donghu by five submerged macrophytes.Desalination，2009，242(1):193-204.

［11]Yang T，Liu Q，Chan L et al.Magnetic investigation of heavy metals contamination in urban topsoils around the East Lake，Wuhan，China.Geophysical Journal International，2007，171(2):603-612.

［12]宁龙梅，王学雷，吴后建.武汉市湿地景观格局变化研究.长江流域资源与环境，2005，14(1):44-49.

［13]宁龙梅，王学雷，胡望斌.利用马尔科夫过程模拟和预测武汉市湿地景观的动态演变.华中师范大学学报:自然科学版，2004，38(2):255-258.

［14]邬建国.景观生态学:第2版.北京:高等教育出版社，2007:260.

［15]Nagendra H，Munroe DK，Southworth J.From pattern to process:landscape fragmentation and the analysis of land use/land cover change.Agriculture，Ecosystems＆ Environment，2004，101(2):111-115.

［16]Andren H.Corvid density and nest predation in relation to forest fragmentation:a landscape perspective.Ecology，1992，73(3):794-804.

［17]Hynes H.The stream and its valley.Verhandlungen der Internationalen Vereinigung fur theoretische und angewandte Limnologie，1975，19:1-15.

［18]Neville HM，Dunham JB，Peacock MM.Landscape attributes and life history variability shape genetic structure of trout populations in a stream network.Landscape Ecology，2006，21(6):901-916.

［19]Wiens JA.Riverine landscapes:taking landscape ecology into the water.Freshwater Biology，2002，47(4):501-515.

［20]Ficetola，GF，Marziali L，Rossaro B et al.Landscape-stream interactions and habitat conservation for amphibians.Ecological Applications，2011，21(4):1272-1282.

［21]Mark B，Michael L.Testing predictions of stream landscape theory for fish assemblages in highly fragmented watersheds.Folia Zoologica，2010，59(3):231-239.

［22]Riitters KH，O'neill RV，Hunsaker CT et al.A factor analysis of landscape pattern and structure metrics.Landscape Ecology，1995，10(1):23-39.

［23]Turner MG.Landscape ecology:the effect of pattern on process.Annual Review of Ecology and Systematics，1989，20:171-197.

［24]王 艳，吴军玲，王恒亮等.武汉近50年来降雨数据的统计分析.湖北工业大学学报，2007，21(6):98-100.

［25]Wulder MA，Franklin SE.Understanding forest disturbance and spatial pattern:remote sensing and GIS approaches.USA:CRC Press，2006:246.

［26]Benson BJ，MacKenzie MD.Effects of sensor spatial resolution on landscape structure parameters.Landscape Ecology，1995，10(2):113-120.

［27]梅安新，彭望琭，秦其明.遥感导论.北京:高等教育出版社，2001:324.

［28]McGarigal K，Marks BJ.Spatial pattern analysis program for quantifying landscape structure.GenTechRepPNW-GTR-351.US Department of Agriculture，Forest Service，Pacific Northwest Research Station，1995.

［29]柳钦火.定量遥感模型，应用及不确定性研究.北京:科学出版社，2010:305.

［30]Zhan Q，Molenaar M，Tempfli K et al.Quality assessment for geo-spatial objects derived from remotely sensed data.International Journal of Remote Sensing，2005，26(14):2953-2974.

［31]Lyons MB，Phinn SR，Roelfsema CM.Long term land cover and seagrass mapping using Landsat and object-based image analysis from 1972 to 2010 in the coastal environment of South East Queensland，Australia.ISPRS Journal of Photogrammetry and Remote Sensing，2012，71:34-46.

［32]Cohen WB，Yang Z，Kennedy R.Detecting trends in forest disturbance and recovery using yearly Landsat time series:2.TimeSync-Tools for calibration and validation.Remote Sensing of Environment，2010，114(12):2911-2924.

［33]Lyons MB，Roelfsema CM，Phinn SR.Towards understanding temporal and spatial dynamics of seagrass landscapes using time-series remote sensing.Estuarine，Coastal and Shelf Science，2013，120:42-53.

［34]Lillesand TM，Kiefer RW，Chipman JW.Remote sensing and image interpretation.UK:John Wiley ＆ Sons Ltd，2004:763.

［35]Wright C，Gallant A.Improved wetland remote sensing in Yellowstone National Park using classification trees to combine TM imagery and ancillary environmental data.Remote Sensing of Environment，2007，107(4):582-605.

［36]Bwangoy JRB，Hansen MC，Roy DP et al.Wetland mapping in the Congo Basin using optical and radar remotely sensed data and derived topographical indices.Remote Sensing of Environment，2010，114(1):73-86.

［37]Kolios S，Stylios CD.Identification of land cover/land use changes in the greater area of the Preveza peninsula in Greece using Landsat satellite data.Applied Geography，2013，40:150-160.

［38]Sexton J，Urban DL，Donhue MJ et al.Long-term land cover dynamics by multi-temporal classification across the Landsat-5 record.Remote Sensing of Environment，2013，128:246-258.

［39]Tulbure MG，Broich M.Spatiotemporal dynamic of surface water bodies using Landsat time-series data from 1999 to 2011.ISPRS Journal of Photogrammetry and Remote Sensing，2013，79:44-52.

［40]Foody GM.Status of land cover classification accuracy assessment.Remote Sensing of Environment，2002，80(1):185-201.

［41]Langford WT，Gergel SE，Dietterich TG et al.Map misclassification can cause large errors in landscape pattern indices:examples from habitat fragmentation.Ecosystems，2006，9(3):474-488.

［42]Cain DH，Riitters K，Orvis K.A multi-scale analysis of landscape statistics.Landscape Ecology，1997，12(4):199-212.

［43]Lausch A，Herzog F.Applicability of landscape metrics for the monitoring of landscape change:issues of scale，resolution and interpretability.Ecological Indicators，2002，2(1):3-15.

［44]邓宏兵，蔡述明，杜 耘等.近50年来江汉湖群水域演化定量研究.长江流域资源与环境，2006，15(2):244-248.

［45]邓宏兵，蔡述明，薛怀平等.江汉湖区土地资源遥感调查与动态监测.华中师范大学学报:自然科学版，2005，39(2):265-268.

［46]林开愚，杨 凯，兰运超等.江汉湖群围垦变化的遥感测定.环境科学学报，1985，5(1):20-29.

［47]陈 锦，武占云，张志超.武汉市湖泊演变与合理开发利用.高等函授学报:自然科学版，2004，16(5):43-47.

［48]邓 帆，李仁东，王海芳.2000～2007年湖北省耕地变化的遥感调查与分析.长江流域资源与环境，2010，19(10):1171-1176.

［49]Malhi Y，Román-Cuesta RM.Analysis of lacunarity and scales of spatial homogeneity in IKONOS images of Amazonian tropical forest canopies.Remote Sensing of Environment，2008，112(5):2074-2087.

［50]Lu D，Hetrick S，Moran E.Land cover classification in a complex urban-rural landscape with QuickBird imagery.Photogrammetric Engineering＆ Remote Sensing，2010，76(10):1159-1168.