基于Landsat TM_OLI影像的哈尔滨市植被覆盖度与景观格局动态变化分析.pdf

第 39 卷 第 5 期2023 年 9 月森 林 工 程FOREST ENGINEERINGVol.39 No.5Sep.,2023doi:10.3969/j.issn.1006-8023.2023.05.008基于 Landsat TM/OLI 影像的哈尔滨市植被覆盖度与景观格局动态变化分析马建明,常晓晴,邢艳秋,王德军(东北林业大学 森林作业与环境研究中心,哈尔滨 150040)摘 要:以哈尔滨市为研究对象,依据 2004、2010、2015、2020 年 Landsat TM/OLI 遥感影像,利用归一化植被指数和像元二分模型方法提取该区域的植被覆盖度,结合植被覆盖度转移矩阵,分析植被覆盖度的时空动态变化。根据植被覆盖度分级结果,利用 Fragstats4.2 软件计算斑块面积、斑块个数、蔓延度与香农多样性指数 4 个景观格局指数,研究哈尔滨市景观格局的时空演变特征。结果表明,哈尔滨市 20042020 年高植被覆盖度和中高植被覆盖度的区域面积累计和分别占总面积的77.81%、72.54%、73.01%、83.31%,整体处于上升趋势。20042020 年高植被覆盖度区域的景观百分比占据主导地位,其次是中高植被覆盖度区域,且大部分植被覆盖度区域斑块数量先增加后减少。从景观层次上,斑块聚合度降低、多样性略增加。探究哈尔滨市植被覆盖度和景观格局的动态变化,对城市规划设计和景观格局优化具有一定的指导和参考意义。关键词:Landsat TM/OLI;归一化植被指数;像元二分模型;景观格局指数;哈尔滨市中图分类号:S718.55 文献标识码:A 文章编号:1006-8023(2023)05-0065-09Analysis of Dynamic Changes of Vegetation Coverage and Landscape Pattern in Harbin Based on Landsat TM/OLI ImagesMA Jianming,CHANG Xiaoqing,XING Yanqiu,WANG Dejun(Research Center of Forest Operations and Environment,Northeast Forestry University,Harbin 150040,China)Abstract:Taking Harbin city as the research object,the paper extracted the vegetation cover based on Landsat TM/OLI remote sensing images in 2004,2010,2015 and 2020,which combined the normalized vegetation index and dimidiate pixel model.Then,we integrated the transfer matrix to analyze the spatial and temporal dynamic changes of vegetation coverage.Based on the results of vege-tation coverage classification,four landscape pattern indices namely patch area,number of patches,contagion index,and Shannons diversity index,were calculated by Fragstats4.2 software to study the spatial and temporal evolution characteristics of the landscape pat-tern in Harbin.The results showed that the sum of areas with high vegetation coverage and medium-high vegetation cover accounted for 77.81%,72.54%,73.01%and 83.31%of the total area in 20042020,respectively,and the vegetation coverage exhibited an o-verall increasing trend in Harbin City.The landscape percentage of high vegetation coverage was the highest from 2004 to 2020,fol-lowed by medium-high vegetation coverage,and the number of patches exhibited a first increasing and then decreasing trend for the most vegetation coverage region.At the landscape level,the aggregation of patches decreased and the landscape diversity increased slightly.This study explores the dynamic changes of vegetation coverage and landscape patterns in Harbin city,which can provide guid-ance and reference for the planning design and landscape pattern optimization of the city.Keywords:Landsat TM/OLI;normalized difference vegetation index;dimidiate pixel model;landscape pattern index;Harbin City收稿日期:2022-11-22基金项目:国家重点研发计划项目(SQ2021YFE010728)。第一作者简介:马建明,硕士研究生。研究方向为定量遥感。E-mail:jianmingma2023 通信作者:邢艳秋,博士,教授。研究方向为定量遥感。E-mail:yanqiuxing 引文格式:马建明,常晓晴,邢艳秋,等.基于 Landsat TM/OLI 影像的哈尔滨市植被覆盖度与景观格局动态变化分析J.森林工程,2023,39(5):65-73.MA J M,CHANG X Q,XING Y Q,et al.Analysis of dynamic changes of vegetation coverage and landscape pattern in Harbin based on Landsat TM/OLI imagesJ.Forest Engineering,2023,39(5):65-73.森 林 工 程第 39 卷0 引言随着全球人口的持续增长,城市用地不断扩张,对陆地生态系统结构和功能以及区域气候产生了重要的影响1。植被作为陆地生态系统最重要的组成要素之一,城市植被可以为居民提供多种环境和社会效益,包括净化空气、降噪、缓解城市热岛效应、减少雨水径流和维护生物多样性,以及美化环境、缓解压力等2。城市植被能够很好地反映城市生态环境状况,是评价城市化发展健康的主要指标。因此,定量监测城市植被时空变化对评价城市生态环境质量、城乡规划、促进城市可持续发展等方面具有重要意义。植被覆盖度(Fractional vegetation cover,FVC)是评价植被生长状况、分布特征和结构变化的重要参数,在城市化研究、植被变化监测和生态环境调查中得到广泛应用3。城市 FVC 的时空动态变化研究可综合反映城市扩张下的植被变化,进而监测城市环境和城市发展。另外,植被景观格局综合反映了植被空间分布及其在环境异质性和干扰状况综合控制下的动态变化特征,已被证明是评估景观模式的有效手段,通常使用景观格局指数定量化表示4。城市快速发展导致建筑用地大面积扩张,侵占原有自然地表植被,带来严重的生态环境问题。应用景观生态学方法研究植被覆盖度景观格局的演变规律,有助于了解城市发展下的植被变化过程,从而揭示植被覆盖情况与自然环境、人为干扰之间的影响因素,为区域发展提供重要依据5。相较于传统的地面测量,遥感技术可以快速大范围监测城市扩张、植被覆盖度和景观格局的变化6。许玉凤等7根据 15 a 的 MODIS NDVI 数据,分析黎平县植被变化,结果表明植被有明显改善;李恒凯等8根据 19902013 年 5 期 Landsat 影像数据,对 23 a 间定南县岭北稀土矿区的植被覆盖情况和景观格局进行了综合分析;苏凯等9利用 20002015 年4 期 MODIS 数据,对东北森林带景观格局的时空变化展开研究,并结合 MCE-CA-Markov 模型,预测东北森林带 2020 年景观格局变化趋势;苏艳琴等10以福建省三明市沙县作为研究区,以 Landsat为数据源,利用线性光谱混合模型分别提取 1996、2004、2014 年的植被覆盖度,并分析其变化趋势;张晓东等11基于生态景观学原理,利用 19892014年 4 期 Landsat 影像数据,分析了宁夏盐池县植被景观格局时空演变特征,并探讨植被与降雨之间的关系;徐凯健等12利用 19752013 年 6 期 Landsat 系列影像数据,以福建省长汀县为研究区,分析植被覆盖度和景观格局在不同时期的变化情况,并探讨人类活动以及政策的影响。但是上述研究主要集中在生态脆弱区和丘陵地貌的植被覆盖度动态变化研究,而结合植被覆盖度,从景观生态学角度对植被景观格局演变研究较少,尤其是人类活动较为集中的城市区域。因此,开展城市区域植被覆盖度及景观格局变化的研究,对掌握城市发展进程和景观生态保护具有重要意义。哈尔滨市位于我国东北部,地处松嫩平原,是我国最北的省会城市,面积最大,但城市化水平远远低于我国的中南部城市。近些年由于城市的经济发展、新城区规划、棚户区改造以及生态恢复政策实施等影响,生态格局发生重大变化13。本研究以 2004、2010、2005、2020 年 4 期 Landsat TM/OLI 遥感影像作为数据源,利用归一化植被指数和像元二分法提取哈尔滨市的植被覆盖度,结合景观生态学理论分析在快速城市化背景下植被覆盖度的时空变化和景观格局的演变,有助于认清生态景观演替和机制,为该区域今后的城市建设和景观格局优化提供科学依据。1 研究区概况及数据源1.1 研究区概况哈尔滨市位于黑龙江省南部,地理范围是1254213010 E,44044640 N,是中国重要的制造业基地。该地区属中温带大陆性季风气候,年平均温度 3.6,冬长夏短,降水集中在 69 月份,年平均降水量 569.1 mm。整体地势东高西低,东以山丘为主,西以平原为主,四周有大量农田包围。土壤种类众多,以黑土为主。本研究区域为哈尔滨市的中心城区,具体范围包括道里区、南岗区、道外区、平房区、香坊区和松北区,总面积约2243.84 km2,如图1 所示。1.2 数据源及预处理本研究采用 Landsat TM/OLI 影像数据,空间分辨率为 30 m,采集时相分别为 2004 年 9 月 21 日(TM)、2010 年 9 月 22 日(TM)、2015 年 9 月 20 日(OLI)和 2020 年 9 月 15 日(OLI)。本研究忽略在较短时间内由植被生长状况不同对植被覆盖度的影响,选取的 4 个影像中植被属于同一物候期,植被覆盖度的时空变化具有可比性14。66第 5 期马建明,等:基于 Landsat TM/OLI 影像的哈尔滨市植被覆盖度与景观格局动态变化分析126200 EN126400 E126200 E12700 E126400 E45400 N4600 N?46200 N45400 N4600 N46200 N12700 Ekm0?10?20?40图 1 研究区位置图(2020 年 9 月 15 日 Landsat OLI影像,假彩色合成)Fig.1 The location of study area(Landsat OLI image on September 15,2020,false color synthesis)为消除遥感影像因大气传输、平台运行状况及传感器系统等因素造成的误差,本研究利用 ENVI5.3对遥感影像数据进行辐射定标、几何校正和大气校正等预处理,并以 150 000 地形图数据为标准,对 4个影像进行几何精校正,校正精度控制在 1 个像素内,最后利用哈尔滨市中心城区行政边界分别裁剪4 期影像数据。2 研究方法2.1 归一化植被指数归一化植被指数(Normalized Difference Vegeta-tion Index,NDVI)是一种简单有效的量化绿色植被的指标,根据植被的光谱特征,将近红外和红波段反射率进行组合,突出表现植被的遥感特征信息,同时抑制土壤反射或大气等干扰因素的影响,能够较好地反映植被生长状态,被广泛用于植被动态变化监测、农作物长势估算以及植被特征识别等方面15,计算见式(1)。NDVI=(NIR-R)/(NIR+R)。(1)式中:NDVI为归一化植被指数;NIR为近红外波段;R为红波段。NDVI的取值范围是-11,其中 NDVI的负值(接近-1 的值)对应于水域,接近于 0 的值通常表示岩石、裸地或雪的贫瘠区域。正值表示有植被覆盖,值越大表明植被覆盖度越高。2.2 植被覆盖度估算模型遥感影像的每个像元可能由多个组分构成,而每个组分对遥感传感器所观测的信息都有贡献,因此可以将遥感信息分解,建立像元分解模型来估算植被覆盖度。像元二分模型(Dimidiate Pixel Model,DPM)因其简便高效成为混合像元模型法中最常用的植被覆盖度估算模型之一,具有明确的物理含义16。假设影像中每一个像元记录的地物光谱信息只由植被和非植被 2 部分组成,所得的光谱信息也只由这 2 个组分因子线性合成,其各自的面积在像元中所占的比率即为各因子的权重,见式(2)17。NDVI=fNDVIveg+(1-f)NDVIsoil。(2)式中:f 为混合像元中植被所占的比例即为该像元的植被覆盖度;NDVI为混合像元的归一化植被指数;NDVIveg和 NDVIsoil分别为全植被覆盖和裸土对应的NDVI。由于纯像元受土壤、植被类型和叶绿素含量等因素的影响,NDVIveg和 NDVIsoil存在很大的不确定性。为此,对4 期 Landsat 影像提取的 NDVI 值进行直方分布图统计分析,以95%的置信度(影像中 NDVI 值对应的像元数量的累计百分比)确定 NDVI 的上下限阈值,分别近似代表裸土和纯植被覆盖区域的NDVI,取值见表1,此方法可以有效减少遥感影像中的噪声等产生的误差,并在 FVC 研究中得到了广泛使用18-20。植被覆盖度的数学表达形式见式(3)。f=(NDVI-NDVIsoil)(NDVIveg-NDVIsoil)。(3)表 1 20042020 年 NDVIveg和 NDVIsoilTab.1 The values of the NDVIveg and NDVIsoilfrom 2004 to 2020植被指数Vegetation index2004 年2010 年2015 年2020 年NDVIveg0.7480.6960.8730.803NDVIsoil-0.009-0.003-0.037-0.5682.3 景观格局评价方法景观指数高度浓缩了景观格局信息,反映了景观的结构组成和空间配置特征21。参考相关研究,结合研究区域的具体特征和研究的实际需求,从斑块面积、斑块个数、蔓延度与多样性指数 4 个方面定量描述研究区植被覆盖度景观格局特征的时空变化,各景观指数具体含义和计算方式见表 2。利用景观格局分析软件 Fragstats 4.2 计算出哈尔滨市植被覆盖度的景观格局指数。76森 林 工 程第 39 卷表 2 景观指数含义及其计算方法Tab.2 Description and calculation of landscape metrics景观指数Landscape metrics描述Description计算方法Calculation斑块面积百分比(PLAND)Percentage of patch area某一斑块类型占整个景观的面积的相对比PLAND=nj=1aijAaij为斑块 ij 的面积;A 为所有景观的总面积斑块个数(NP)Number of patches某一斑块类型的斑块总个数NP=n,n 为斑块个数蔓延度(CONTAG)Contagion index景观里不同斑块类型的团聚程度或延展趋势CONTAG=1+mi=1mk=1(Pi)gikmk=1gik()ln(Pi)gikmk=1gik()2ln(m)gik为 i 类型斑块和 k 类型斑块毗邻的数目香农多样性指数(SHDI)Shannons diversity index景观的丰富程度SHDI=-mi=1(pilnpi)pi为景观类型 i 所占的面积比例;m 为景观类型数目 3 结果与分析3.1 植被覆盖度变化特征本研究首先利用 Landsat TM/OLI 影像提取了哈尔滨市 2004、2010、2015、2020 年 4 期的 NDVI,其空间分布如图 2 所示。随后利用像元二分法得到哈尔滨市 2004、2010、2015、2020 年的植被覆盖度。为进一步了解哈尔滨市植被覆盖度的空间分布情况,参照 Cai 等3的分级方法,并结合研究区的实际情况,将研究区植被覆盖度划分为 5 个等级:无植被覆盖区(0fc0.2)、低植被覆盖区(0.2fc0.4)、中植被覆盖区(0.4fc0.6)、中高植被覆盖区(0.6fc0.8)、高植被覆盖区(0.8fc1),植被覆盖度空间分布如图 3 所示。由图 3 可以看出,各等级植被覆盖度空间分布具有明显的空间差异性。哈尔滨市区以高植被覆盖分布为主,主要分布在城区的北、东和西 3 个片区,这 3 个片区主要是以林业、农业用地为主,9 月份为农作物的成熟期,植被覆盖度较高;从城区内部来看,以松花江为界线,江北城区整体植被覆盖度均高于江南城区,这可能是因为江北城区属于新开发区,规划合理,同时政府也更加注重城区的生态环境,而江南城区属于老城区,植被覆盖度低于江北城区。图例（a）2004（b）2010kmkm0?10?200?10?20NN图例高：0.992?799低：-0.841?065高：0.993?314低：-0.845?77186第 5 期马建明,等:基于 Landsat TM/OLI 影像的哈尔滨市植被覆盖度与景观格局动态变化分析（c）2015（d）2020kmkm0?10?200?10?20NN图例图例高：0.873?440低：-0.948?306高：0.903?922低：-0.911?765图 2 20042020 年哈尔滨市中心城区 NDVI 分布图Fig.2 The maps of NDVI distribution in Harbin center city from 2004 to 2020无植被覆盖区低度植被覆盖区中度植被覆盖区中高度植被覆盖区高度植被覆盖区（a）2004（b）2010（c）2015（d）2020kmkmkmkm0?10?200?10?200?10?200?10?20无植被覆盖区低度植被覆盖区中度植被覆盖区中高度植被覆盖区高度植被覆盖区无植被覆盖区低度植被覆盖区中度植被覆盖区中高度植被覆盖区高度植被覆盖区无植被覆盖区低度植被覆盖区中度植被覆盖区中高度植被覆盖区高度植被覆盖区NNNN图 3 20042020 年哈尔滨市植被覆盖等级图Fig.3 The maps of vegetation coverage grade in Harbin from 2004 to 202096森 林 工 程第 39 卷 由表 3 可知,20042020 年各等级植被覆盖度区域的面积发生了明显变化。其中,高植被覆盖度和无植被覆盖度的区域面积呈现逐渐增长的趋势,而低、中和中高植被覆盖度区域的面积均逐渐减少。整体来看,在 20042020 年研究区高植被覆盖度和中高植被覆盖度的区域面积累计和占总面积的比例分别为 77.81%、72.54%、73.01%、83.31%,说明哈尔滨市整体植被覆盖情况较好,呈增加趋势。相较于 2004、2010 和 2015,2020 年的植被覆盖情况最好,其中高植被覆盖度地区的面积为 1 647.09 km2,占区域面积的67.40%。20042020 年高、无植被覆盖度区域的面积分别增加了 383.50、91.07 km2,而低植被覆盖度、中植被覆盖度和中高植被覆盖度区域的面积累计减少了474.57 km2,占比19.42%,可以发现高植被覆盖度区域的面积要远超过无植被覆盖区域的面积,说明哈尔滨市总体生态环境向好。3.2 不同植被覆盖等级间的动态变化分析为了进一步分析不同等级植被覆盖度区域面积的变化,本研究计算了 20042010 年、20102015 年和 20152020 年的转移矩阵,见表 4。表 3 各等级植被覆盖度面积Tab.3 The area of each vegetation coverage gradekm2年份Year无植被覆盖Non-vegetation coverage低植被覆盖Low vegetation coverage中植被覆盖Medium vegetation coverage中高植被覆盖Medium-high vegetation coverage高植被覆盖High vegetation coverage200448.59197.49296.14638.031 263.592010128.05241.13301.97620.911 151.782015146.37231.88281.26529.331 255.012020139.6639.36228.86388.871 647.09表 4 20042020 年不同等级植被覆盖转移矩阵Tab.4 Transfer matrix of vegetation coverage at different levels from 2004 to 2010km2时期Period无植被覆盖Non-vegetation coverage低植被覆盖Low vegetation coverage中植被覆盖Mediumvegetation coverage中高植被覆盖Medium-highvegetation coverage高植被覆盖High vegetation coverage总计Total20042010无植被覆盖34.6035.1016.8325.4716.04128.05低植被覆盖10.45108.1161.6033.0427.93241.13中植被覆盖2.5740.29107.0594.1257.94301.97中高植被覆盖0.7210.4570.88236.82302.04620.91高植被覆盖0.243.5539.77248.58859.651 151.78总计48.59197.49296.14638.031 263.592 243.8420102015无植被覆盖90.7517.5012.5610.7214.85146.37低植被覆盖27.68122.0439.0222.5820.56231.88中植被覆盖6.5174.86105.9754.0039.92281.26中高植被覆盖2.5622.4198.65200.20205.51529.33高植被覆盖0.554.3345.78333.41870.941 255.01总计128.05241.13301.97620.911 151.782 243.8420152020无植被覆盖99.664.704.9114.1316.26139.66低植被覆盖14.673.754.138.098.7239.36中植被覆盖26.11132.4330.8320.9518.53228.86中高植被覆盖4.0381.42167.9477.2158.26388.87高植被覆盖1.899.5773.44408.951 153.241 647.09总计146.37231.88281.26529.331 255.012 243.84 20042010 年植被覆盖度等级未发生变化的区域达到 1 346.23 km2,占总面积的 59.99%。在发生变化的区域中,高植被覆盖度区域转向中高植被覆盖度区域面积最大,转换面积为 302.04 km2,这可能是由于哈尔滨市的发展与扩展,人民生产生活需求量大,从自然界获取的资源量增加,对生态植被不可避免地产生负面影响。20102015 年,中高植被覆盖度主要转变为高植被覆盖度区域,转换面07第 5 期马建明,等:基于 Landsat TM/OLI 影像的哈尔滨市植被覆盖度与景观格局动态变化分析积为 333.41 km2。虽然还存在高植被覆盖度区域减少的现象,但是整体发展良好。另外,由中植被覆盖度区域转化为中高植被覆盖度区域,面积为98.65 km2,由低植被覆盖区域转向中植被覆盖区域,面积为 74.86 km2,由无植被覆盖度区域转向低植被覆盖区域,面积为 27.68 km2,植被覆盖度均有所增加。20152020 年,中高植被覆盖区域向高植被覆盖度区域转化的面积仍最多,共转化了408.95 km2,虽然存在中高植被覆盖度转向中植被覆盖度的区域,但面积较少,仅为 20.95 km2。其次是低植被覆盖度的区域向中植被覆盖度的区域转化,面积为 132.43 km2,说明哈尔滨市的植被覆盖度呈现上升的趋势。3.3 哈尔滨市景观格局分析利用 Fragstats 4.2 计算得到哈尔滨市的植被覆盖度的景观指数值。由图 4(a)可以看出,20042020 年高植被覆盖度区域占据了景观中的主导地位,其次为中高植被覆盖度。研究区域内的不同时期,无植被覆盖度区域的景观百分比逐渐增加,低植被覆盖度和中植被覆盖度的景观百分比属于先增加后减少。而中高植被覆盖度的景观百分比一直处于减少状态,但依然处于优势地位。结合图 4(b)斑块个数分析,大多数植被覆盖度区域斑块个数呈现先增加后减少的趋势,说明景观的破碎化程度先增加后减少,最后连片性增加,即城区扩建形成连片区域。8060402002004 年2010 年2015 年2020 年无值被覆盖2004 年2010 年2015 年2020 年景观百分比（%）Percentage?of?vegetation?coverage（b）斑块个数（b）Number?of?patches（a）景观百分比（a）Percentage?of?vegetation?coverage斑块个数（103）Number?of?patch50403020100低值被覆盖中值被覆盖中高值被覆盖高值被覆盖无值被覆盖低值被覆盖中值被覆盖中高值被覆盖高值被覆盖图 4 哈尔滨不同时期植被覆盖度景观百分比和斑块个数Fig.4 Percentage of vegetation coverage and number of patches for each period in Harbin city 多样性指数主要反映景观的异质性和丰富程度。由表 5 可知,20042020 年香农多样性指数先增加后保持稳定,说明哈尔滨市的植被斑块多样性上升后保持平稳,各植被覆盖类型异质性逐渐增强,破碎化程度高。另外,蔓延度指数由 2004 年的53.51%下降至 2020 年的 43.98%,表明植被斑块的表 5 研究区不同时期多样性指数和蔓延度指数Tab.5 Diversity index and spread index in different periods in the study area年份Year香农多样性指数Shannons diversity index蔓延度(%)Contagion index20041.2353.5120101.6745.5620151.9645.4220201.9643.98分布逐渐分散和复杂,景观破碎度加剧,团聚程度不高,景观中大斑块的面积被蚕食减小。虽然哈尔滨市区主要以城区建设为主,但近年来随着旧城区改造和新城区建设以及人们对绿化的重视,景观结构组成趋向复杂化。4 讨论与结论城市的扩张必然以牺牲绿地为代价,但随着社会经济发展,保护环境意识的提高以及各种城市绿化政策的实施,城市区域植被覆盖度逐渐增加。本研究基于长时间遥感影像数据进行城市植被覆盖度和植被景观格局分析,结果表明,随着哈尔滨城市的发展,植被覆盖度逐渐增加,说明近年来哈尔滨政府实施各种植被保护政策取得了显著性的成效,“哈尔滨市城市总体规划”“哈尔滨市城市绿地17森 林 工 程第 39 卷条例”等措施的实施有效保护和增加了绿地面积,改善了生态环境,这与 Dao 等22对哈尔滨土地利用变化研究得出绿地面积逐渐增加的结果基本一致,同时李继红等23研究也表明哈尔滨市近十年植被碳储量明显增强。另外,在本研究中选择分辨率较精细的 TM 影像存在着一定的局限性,城市绿地的斑块较小且破碎,其变化不易检测,虽然不影响结论的正确性,但可能会影响结果的精度24,可在未来研究中借助更高分辨率的遥感影像(如高分系列,WorldView-2 影像等)进行深入研究。Song等25对哈尔滨市中心区景观格局特征研究也同样表明该地区生态斑块面积小,碎片性强。此外,本研究分析了哈尔滨市植被覆盖度和景观格局动态变化,但其变化的驱动因素和影响机制并未考虑。在快速推进城市化进程同时,如何增强植被的生态环境效应,维持城市的可持续发展还需进一步研究。本研究利用 Landsat TM/OLI 影像数据,结合归一化植被指数和像元二分模型,对哈尔滨市 20042020 年内的植被覆盖度进行估算。最后,利用植被覆盖度和景观生态学理论分析哈尔滨景观格局演变,结论如下。1)整体而言,哈尔滨市整体植被覆盖度较好。20042020 年高植被覆盖度和中高植被覆盖度的区域面积之和占总 区 域 面 积 分 别 为 77.81%、72.54%、73.01%、83.31%。不同等级植被覆盖的变化趋势不同,在 20042020 年,高植被覆盖度区域和无植被覆盖度区域变化较大,整体呈现增长趋势。而低、中和中高植被覆盖度区域均有下降趋势。从空间分布来看,哈尔滨市主城区由于建筑物比较集中,植被覆盖度较低,而城郊植被覆盖度较高。2)结合转移矩阵方法,分析了 20042020 年哈尔滨市不同植被覆盖度等级间的动态变化。植被覆盖度的转化整体趋于良性发展,但存在部分区域由于新城区建设和人类生活影响造成的植被覆盖度下降,需要引起相关部门的重视。3)利用景观生态学理论,结合植被覆盖度,分析了哈尔滨市植被覆盖度景观格局动态变化。研究表明 20042020 年高植被覆盖度区域占据了优势地位,其次是中高植被覆盖度区域。大部分植被覆盖区域的斑块数量呈现先增加后减少的趋势。总体来看,20042020 年哈尔滨市植被景观多样性逐渐增加,植被斑块的分布逐渐分散和复杂。【参 考 文 献】1 ZITER C D,PEDERSEN E J,KUCHARIK C J,et al.Scale-dependent interactions between tree canopy cover and impervious surfaces reduce daytime urban heat during summerJ.Proceedings of the National Academy of Sci-ences of the United States of America,2019,116(15):7575-7580.2 MITCHELL M G E,DEVISSCHER T.Strong relationships between urbanization,landscape structure,and ecosystem service multifunctionality in urban forest fragments J.Landscape and Urban Planning,2022,228:104548.3 CAI Y F,ZHANG F,DUAN P,et al.Vegetation cover changes in China induced by ecological restoration-protec-tion projects and land-use changes from 2000 to 2020J.CATENA,2022,217:106530.4 SMIRAGLIA D,CECCARELLI T,BAJOCCO S,et al.Un-raveling landscape complexity:land use/land cover changes and landscape pattern dynamics(1954-2008)in contrasting peri-urban and agro-forest regions of northern ItalyJ.En-vironmental Management,2015,56(4):916-932.5 高彦净,谢余初,钱大文,等.甘肃白龙江流域植被覆盖度及景观格局变化J.水土保持研究,2015,22(1):181-187.GAO Y J,XIE Y C,QIAN D W,et al.Dynamic varia-tions of vegetation coverage and landscape pattern in Bail-ongjiang Basin of southern GansuJ.Research of Soil and Water Conservation,2015,22(1):181-187.6 HEROLD M,SCEPAN J,CLARKE K C.The use of remote sensing and landscape metrics to describe structures and changes in urban land usesJ.Environment and Planning A:Economy and Space,2002,34(8):1443-1458.7 许玉凤,陆承红,潘网生,等.基于 MODIS 的贵州省黎平县植被变化J.森林工程,2017,33(6):6-10,14.XU Y F,LU C H,PAN W S,et al.MODIS-Based vege-tation change in Liping County of GuizhouJ.Forest En-gineering,2017,33(6):6-10,14.8 李恒凯,雷军,杨柳.基于 Landsat 影像的离子稀土矿区植被覆盖度提取及景观格局分析J.农业工程学报,2016,32(10):267-276.LI H K,LEI J,YANG L.Extraction of vegetation cover-age and analysis of landscape pattern in rare earth mining area based on Landsat imageJ.Transactions of the Chi-nese Society of Agricultural Engineering,2016,32(10):267-276.9 苏凯,王茵然,