第33卷，第1期国土资源遥感VoL33,No.l 2021年3月REMOTE SENSING FOR LAND&RESOURCES Mae,2021
doi:10.6046/gtzyyy.2020178
引用格式：田雷，傅文学，孙燕武，等.基于TM影像的西伯利亚北方森林覆盖度近30a空间变化研究［J］.国土资源遥感，2021,33（1）：214-220.（Tian L,Fu W X,Sun YW,et ai.Research on spatial change of the boreal forest cover in Siberio over the past 30years based on TM images）J*.Remote Sensing for Land and Resources,2021,33（1）：214-220.） 基于TM影像的西伯利亚北方森林覆
盖度近30a空间变化研究
田雷1>2，傅文学1，孙燕武1，荆林海1，邱玉宝1，李新武1
（1.中国科学院空天信息创新研究院数字地球重点实验室，北京100094；
2.安徽理工大学空间信息与测绘工程学院，淮南232001）
摘要：在全球变暖的背景下，研究西伯利亚北方森林覆盖的长时间尺度空间变化特征不仅对全球气候变化研究和可持续发展有着重要意义，也为进一步研究北方森林变化对气候变化的响应提供支撑。以1985年和2015年2期Landsat TM/OLI为数据源，选取俄罗斯克拉斯诺亚尔斯克边疆区为西伯利亚北方森林典型研究区，采用决策树分类法得到研究区2期森林覆盖分类图，并采用高分二号（GF-2）影像随机选点验证，分类精度达到94.53%。对2期森林覆盖分类图在N51。~69。纬度范围内以2。为区间进行纬度分割并开展空间叠加分析，定量化分析每个纬度带内森林覆盖空间变化信息及其空间变化规律。结果表明：近30a来，西伯利亚北方森林变化显著，森林总体覆盖度由1985年的75.42%增加到2015年的80.53%，增加了5.11百分点。同时，不同纬度带林地面积变化率有较大的差异：N67。~69。纬度带林地面积变化率最高,N63。~65。纬度带次之,N57。~59。纬度带最低。总体上看，森林覆盖度在各个纬度内都出现了不同程度的增加，在N63。~67。纬度带内增加最显著；N57。~63。纬度带森林覆盖度变化相对平稳；而351。~57。纬度带内森林覆盖度增加量随着纬度的降低而降低。 湖北农业科学关键词：TM影像；北方森林；森林覆盖度；西伯利亚；空间变化
中图法分类号：TP79文献标志码：A文章编号：1001-070X（2021）01-0214-07
0引言
全球变暖导致世界气候产生显著变化,1951—2012年间，地球表面的平均温度上升了约0.72P，石化燃
料的大量使用以及大规模毁林所导致的大气中温室气体尤其是二氧化碳（CO2）浓度的不断升高是造成全球气候变化的主因［1］o据《2015年全球森林资源评估报告》显示：全球森林面积为39.99亿hm2，是陆地重要的碳库）2*。作为地球上分布面积最广、结构最为复杂的陆地生态系统，全球森林在2006—2015年吸收了超过1/4的全球碳排放⑶。因此，森林是全球碳循环系统中一个重要的组成部分，能够对大气层温室气体的浓度带来重大影响，对于全球气候调节起着无可替代的作用［4_6］（
北方森林区约占全球森林面积的30%，主要分布于北美和欧亚大陆，其中：60%位于俄罗斯,28%位于加拿大，其余12%分布于中国、芬兰、挪威、瑞典等10多个北半球国家［7-8］o作为巨大的碳库，北方森林在全球碳平衡和应对气候变化中具有重要的地位，同时也是受气候变化最敏感的地区之一，其动态变化是促进全球变化重要的因子，因此在与全球气候系统的耦合中，存在极强的交互作用。在全球变暖背景下，监测北方森林长时间尺度的空间变化特征，对全球气候变化研究和可持续发展有着实际意义。
森林覆盖变化主要体现在受人为和自然因素导致的林地与非林地之间的变化［9］（遥感技术的快速发展，使得获取大范围地区长时间序列的森林覆盖变化信息成为可能，极大地改善了过去依靠外业耗时耗力的局面。国内外研究学者已在利用多源遥感技术监测森林覆盖变化方面做了大量研究，在不
收稿日期：2020-06-17；修订日期：2020-07-28
基金项目：中国科学院战略先导科技专项项目（A类）子课题“国际三极科学数据共享网络”（编号：XDA19070102）和国家重点研发计划项目“综合自然和人类活动影响的土地利用变化遥感监测研究”（编号：2017YFE0100800）共同资助。
第一作者：田雷（1994-），男，硕士研究生，主要从事遥感图像处理与应用。Emaii:echard_dada@outlook,.om o
通信作者：傅文学（1977-），男，博士，副研究员，主要从事遥感机理与应用研究。Emaii:fu^xfaircas.ac 。
第1期田雷，等：基于TM影像的西伯利亚北方森林覆盖度近30c空间变化研究•215•
同尺度上发展了一系列实用可靠的技术方法，并取得较好的结果。在低空间分辨率数据方面，覃先林等［10*利用MODIS数据，分别采用红光-近红外法、共生矩阵纹理法和基于相似度的变化监测方法，对我国东北林区的森林覆盖变化监测方法进行研究，对比验证表明基于相似度的变化监测方法正确度最高(80%以上)。但由于影像空间分辨率低，结果误差大，较适宜于对大区域的森林覆盖变化进行监测(高空间分辨率数据覆盖范围小，更适用于小尺度精细分析。王荣等〔⑴以IKONOS影像为数据源，利用面向对象多尺度分割算法、Sobel算子边缘检测及骨架线提取等方法，构建纹理线条密度指数(TLDI)对纳版河自然保护区的天然林与人工林进行分类，研究表明：与常用的植被特征指数和纹理特
征指数相比,TLDI指数的离散度更好、分类效果更佳，可以有效地提取森林内部天然林与人工林植被覆盖信息; Wilson等〔12*利用归一化植被指数(nomVized dlfer-encc veyetation index,NDVI)和归一化差异水分指数(normalized dtferencc water index,NDWI)对时间序列Landsat专题图图像进行分类以检测森林采伐水平并比较了两者的多重分类精度；任冲等［⑶以1988#2015年5期夏季Landsat TM/OLI遥感影像为主要数据源，利用随机森林和参数优化支持向量机分类器对土地覆盖类型进行分类，获得了天水市近30a林地动态变化信息；姜洋等〔⑷利用不同时相的Landsat TM/ETM+为数据源，采用面向对象和基于多级决策树的分类方法得到了浙江省2000年、2005年'2010年的森林植被覆被图，并通过专题图分析获取了浙江省森林资源动态变化的空间分布信息。这些研究的结果表明Landsat系列数据可以兼顾大尺度与较精细分辨率的优势，在大范围的森林研究中发挥重要的作用。但是以上研究工作只考虑了森林整体的变化，没有顾及到森林变化随纬度的空间变化规律。
本文选取西伯利亚森林受人类活动影响较少的典型研究区，定量化研究近30a森林覆盖度的总体变化及随纬度的空间变化规律，为分析该区域对全球变化的响应提供参考。
1研究区概况及数据源
1.1研究区概况
本文选取西伯利亚北方森林典型区俄罗斯克拉斯诺亚尔斯克边疆区北方森林区为研究区，地理位置见
图1。研究区位于西伯利亚中部，空间范围在E81。〜1021N51。〜69。之间。该研究区的气候为典型的大陆性气候，气温变化剧烈，冬季漫长，夏季短促。1月份平均温度北部达-36P，南部达-18P；7月份的平均温度北部为10P，南部为20P，年平均降水量约316mm。
图1研究区地理位置示意图
Fig.1Location of the study area
北方森林地理区系主要由云杉(Pices asperata Masi.)、松树(Pinue)、冷杉(Abies fabrs(Mast.) Craii)和落叶松(Lafi gmelini i(Rupc.)Kuzee)等针叶树以及桦木(Betula)、杨树(PopulusL)、柳树(Sal)和桤木(Alnus cremastogyne BurO)等阔叶树组成。尽管北方森林生态系统物种组成相对简单，但北方森林区无论是平均温度还是降水均表现出了与全球一致的变化趋势，并且变化幅度要比全球平均趋势更为剧烈小16*。该研究区人类活动较少，只在N53。〜57。纬度带内存在少量的城市及农田区，因此研究区内森林的变化总体上可以体现自然因素的影响。1.2数据源
黑暗中的幽魂采用的数据为完全覆盖研究区的1985年Landsat5TM与2015年Landsat8OLI2期夏季时相遥感影像，数据来源于USGS(http：//v/)，共149景影像。影像时间主要集中在6—9月，云量小于5%，产品级别L1T。对于当年数据缺失的影像，采用相邻或邻近年份的数据补全，时相差异不超过2a。
2数据处理、分类及评价
2.1数据预处理
影像数据预处理包括辐射定标、大气校正、几何
•216•国土资源遥感2021年
纠正、去云处理、镶嵌以及裁剪等，如图2所示。其中影像去云是通过薄云优化变换(haze optimized transformation,HOT%进行遥感影像薄云识别与去除。主要原理是利用影像中无云区地物的蓝光波段和红光波段的高度相关性确定晴空线，然后计算云区的像元相对于晴空线的偏移距离(HOT%，最后依据HOT的大小实现薄云的自动探测与去除［17*( RST计算公式为：
HOT=P blue S j'一P red C os0-I^cos0,(1%式中：(LUE和(ED分别为TM/OLI影像的蓝光波段和红光波段的反射率；/为晴空线的截距；'为晴空线的倾度。
师OLI影像L陽射定标气校正
顾嵌及裁剪I—去云处理IT几何纠正|
图2TM/OLI影像预处理流程
Fig.2Preprocessing flowchart of TM/OLI scenes 2.2森林覆盖分类
森林覆盖分类采用决策树分类法〔18*，根据影像的不同特征，利用训练空间实体集产生规则和发现规律；其次根据不同取值建立树的分支，在每个分支子集中重复建立下层结点和分支，以树型结构表示
分类或决策集合；最后形成决策树，从而将研究区划分为有林地和无林地2种地类。选用归一化植被指数(normalized difTeTecc vevetation index,ND-VI%、比值植被指数(ratio vevetation index,RVI%以及近红外波段进行有林地和无林地的分类。NDVI作为遥感估算植被信息中最常用的植被指数，能有效地削弱复杂地形对影像信息提取的不利影响，增强植被与其他地类的区分度，有助于提高森林信息提取的精度和可信度，其公式为：
NDVI二(P nir一P red)/(P nir+P red),(2%式中P nir为影像的近红外波段的反射率。
RVI能较好地反映植被生长状况和覆盖度的差异，特别适用于植被生长旺盛、具有高覆盖度的植被监测,<V计算公式为：
<V#P nir/P red°(3%以上2种植被指数可有效区分有林地和无林地，但对于部分植被覆盖度较高的无林地(如农用地)分类效果较差。研究发现不同植物在近红外波段(0.76〜0.90)m%的光谱反射峰值不同〔19*，近红外波段对不同植物叶片内部结构不同引起的反射率差异甚为敏感，因此可用近红外波段来进行有林地和其他植被的区分。根据这一特性，本研究通过样本训练确定有林地和农用地在近红外波段的最佳区分阈值范围，采用决策树算法进行分类，如图3所示。基于样本训练确定定量的决策树分类规则: NDVI值大于0.52且RVI值大于6作为植被对象，否则作为非植被对象；选取上层分类的植被对象的近红外波段值小于3788(反射率经大气校正后扩大了104倍)的作为有林地对象，否则作为非林地
对象。最后，将第一层非植被对象和第二层无林地对象合并作为最终的无林地分类结果，第二层有林地对象作为最终有林地分类结果°
图3TM/OLI影像决策树分类模型
Fig.3TM/OLI scenes classification based on
decision tree classification model
2.3精度评价
为验证分类精度和可靠性，选取了2015年同时期4景高分辨率遥感影像高分二号数据(全谱段空间分辨率为0.81m%进行分类结果验证，在高分二号影像上随机选取均匀分布的987个有代表性的独立检验样本进行精度验证。选用总体分类精度和Kappa系数作为验证指标。
1%总体分类精度。总体分类精度是指被正确分类的像元总数占总像元数的比例，即
N#("N)/N,(4%
"二1
式中：N为分类的总体精度；*为分类的类别数；N 为样本总数；N为第，类被正确分类的样本数目。
2%总体Kappa系数°Kappa系数是由Cohen在1960年提出的用于评价遥感图像分类正确程度和比较图件一致性的指数〔20*，即卩
Kappa#(N-N)/(1-N%,(5%
式中:71,a，…，a*分别为每一类的真实样本个数; L1,L，…,L*分别为预测的每一类样本个数
。
第1期田雷，等：基于TM 影像的西伯利亚北方森林覆盖度近30 a 空间变化研究
• 217 •
3结果与分析
3.1分类结果
对 1985 年、2015 年 Landsat TM/OLI 影像采用 决策树法进行分类，并将分类结果中的图斑进行主
要分析处理，以去除图斑噪声的影响，结果如图4 所示。通过精度评价结果表明本研究所采用的决策 树分类方法总体分类精度为94. 53% , Kappa 系数为 0.870，分类结果可靠，可满足于大区域、复杂地形的
植被/森林覆盖度变化监测。
(a ) 1985 年 (b ) 2015 年
图4 1985年和2015年研究区森林覆盖
Fig. 4 Forest cover in 1985 and 2015
3 - 2结果分析
为了更好地描述研究近30 a 西伯利亚北方森 林覆盖度的空间变化信息，对研究区进行2。纬度带 间隔分割，并分别统计每个纬度带内森林覆盖信息, 统计信息见表1。其中，面积变化率表示1985—
2015年有林地面积增加量与1985年有林地面积之
比；覆盖度变化表示研究区单位面积的林地变化 量，从林地面积变化率和覆盖度的增加量上分别体 现有林地的变化程度。
表1 1985年和2015年森林覆盖信息统计
Tab. 1 Statistics of foress cover information in 1985 and 2015
纬度
总面积/km 2
1985 年
2015 年1985—2015 年
面积/km 2
覆盖度/%
面积/km 2
覆盖度/%
面积增加量 akm 2
面积变 化率/%
覆盖度变 化百分点
N67。〜69。19 596.45
5 838.9329.80
填量词6 731.24
34.35892.31
乐果q715.28 4.55N65。〜67。58 813.1832 255.7754.84
40 192.9968.347 937.23
24.61
13.50
N63。〜65。58 364.0940 083.9968.6845 243.3477.52
5 159.3512.878.84N61。〜63。6
6 364.15
55 911.37
84.2556 952.28
85.82
1 040.92
1.86
1.57N59。〜61。130 507.24106 668.6681.73110 896.3684.974 227.70 3.96 3.24
N57。〜59。155 232.73
139 694.2089.99141 250.6790.99
1 556.47
1.11
1.00N55。〜57。114 456.6779 057.5169.07
88 167.88
77.039 110.3811.527.96
N53。〜55。77 579.5853 261.04
68.6558 131.1974.93
4 870.159.14
6.28N51。〜53。
33 564.27
26 059.5877.6427 744.8082.661 685.21
6.47 5.02总计714 478.36538 831.05
75.42
575 310.77
80.5336 479.72
6.77
5. 11
经统计分析得知研究区内森林1985年的森林 为80. 53% o 经过近30 a 的变化，研究区森林面积 面积为538 831.05 km 2，森林覆盖度为75.42%； 增加了 36 479.72 km 2，森林覆盖度增加了 5.11百
2015年的森林面积为575 310.77 km 2，森林覆盖度 分点，森林面积和覆盖度总体上都呈现增加的趋势
。公益海报设计论文
•218•国土资源遥感2021年
对2期森林覆盖情况进行空间变化叠加分析，得到1985#2015年研究区北方森林覆盖变化的空间分布，如图5所示。
65°N
63°N
61°N
59°N
57°N
55°N
53°N
51°N
67°N
85°E90°E95°E100°E
图5森林覆盖动态变化监测结果Fig.5Monitoring results of dynamic changes of forest coverbatedoncla t nfncatnonretultt
由图5和表1可以看出，近30a研究区在N63。〜69。纬度带内面积变化率较高：其中,N67。〜69。纬度带有林地增加面积为892.31km2，面积变化率为15.28%，森林覆盖度增加了4.55百分点；N65。〜67。纬度带有林地增加面积为7937.23km2，面积变化率为24.61%，森林覆盖度增加了13.50百分点；
N63。〜65。纬度带有林地增加面积为5159.35km2，面积变化率为12.87%，森林覆盖度增加了8.84百分点。在这些纬度带内，森林覆盖度增加最为显著，表明对气候的响应更敏感，受气候变化的影响更强。
在N57。〜63。纬度带，研究区的有林地和无林地空间变化相对平稳，从图5上看，整体上表现为无变化的黄区域，森林覆盖度增加缓慢：不同纬度带增加比例分别为1.57（N61。〜63。）,3.24（N59。〜61。），1.00（N57。〜59。）百分点。表明该纬度带森林受气候变化影响较小，面积和覆盖度都呈现出稳定的趋势。
而在N51。〜57。纬度带内有少量的城市和农业区，人类活动对林地变化产生了一定的影响，但总体上森林存在着增加的趋势。统计发现森林覆盖度增加量随着纬度的降低而降低：N55。〜57。纬度带有林地增加面积为9110.38km2，森林覆盖度增加了7.96百分点,N53。〜55。纬度带有林地增加面积为4870.15km2，森林覆盖度增加了6.28百分点,N51。〜53。纬度带有林地增加面积为1685.21km2，森林覆盖度增加75.02百分点。
图5也展示了有林地与无林地的相互转化，其中绿为无林地转化为有林地，表示林地的增加部分；而红为有林地转化为有林地，表示林地的砍伐或退化；黄为无变化，保持了原有林地的部分；灰为无林地。统计分析表明：研究区1985#2015年无林地转化为有林地为7.94%，有林地转化为无林
地为2.83%，有林地面积净增加5.11%。有林地增加区域主要分布在N53。〜55o,N55。~57O,N63O〜65。和N65。〜67。纬度带；有林地减少区域主要分布在N55。〜57。和N57。〜59。纬度带。
由此可见，研究区内北方森林覆盖总体上呈现增加趋势，但不同纬度带的增加率明显不同。已有研究利用AVHRR数据，分析了CO2大气浓度年增幅40%对北半球高纬度地区植被生长增加的驱动过程［21-22*，与本研究西伯利亚北方森林覆盖度变化规律表现出较强的一致性，表明气候变化是北方森林覆盖增长的一个重要驱动因素。同时不同纬度带内温度和降雨的差异性，可能是西伯利亚北方森林覆盖度变化在不同纬度带表现不同的原因，但对于响应机制以及复杂的耦合过程还需进一步研究。
tv搜查线4结论
利用1985年和2015年2期Landsat遥感数据研究了西伯利亚典型研究区北方森林近30a的总体森林面积空间变化和森林覆盖变化规律。同时对研究区进行2。纬度带间隔分割，定量化分析了不同纬度带内森林面积和覆盖度变化，主要结论如下：1%近30a来，西伯利亚北方森林覆盖度变化显著，总体森林覆盖度由1985年的75.42%增加到2015年的80.53%，森林覆盖度增加了5.11百分点。
2%西伯利亚北方森林在不同纬度带内有林地面积变化出现了较大的差异：在N65。〜67。纬度带内有林地面积变化率最高，变化率达到24.61%，森林覆盖度增加了13.50百分点；在N67。〜69。纬度带
有林地面积变化率次之，变化率为15.28%，森林覆盖度增加了4.55百分点；在N57。〜59。纬度带内有林地面积变化率最低，变化率为1.11%，森林覆盖度增加了1.00百分点。
3%总体上看，西伯利亚北方森林有林地面积在各个纬度带内都出现了不同程度的增加，增加区域
豫公网安备41160202000603
站长QQ:729038198
关于我们
投诉建议