西藏林芝市土地利用程度变化及其空间异质性分析

引文格式：郭健斌, 张英, 张志伟, 等. 西藏林芝市土地利用程度变化及其空间异质性分析[J]. 云南农业大学学报(自然科学), 2023, 38(3): 511−519. DOI: 10.12101/j.issn.1004-390X(n).202209064
西藏林芝市土地利用程度变化
及其空间异质性分析*
郭健斌1,2， 张　英
1， 张志伟1， 侯　磊1， 曾维莲2,3
(1. 西藏农牧学院 资源与环境学院，西藏 林芝 860000；2. 西藏农牧学院 西藏乡村振兴研究院，西藏 林芝 860000；
3. 西藏农牧学院 植物科学学院，西藏 林芝 860000)
摘要: 【目的】揭示青藏高原高山峡谷地区土地利用程度空间分异特性，为当地国土空间规划提供理论依据。

【方法】基于2010年和2020年2期Landsat 卫星遥感数据获取30 m 网格土地利用数据，采用空间自相关分析法分析林芝市土地利用程度综合指数及其变化的空间分异特征。

【结果】林芝市较强和强土地利用程度网格单元主要分布在低海拔区域和雅鲁藏布江、尼洋河宽河谷流域，弱和较弱土地利用程度网格单元主要分布在高海拔区域。

2010—2020年，林芝市土地利用程度呈现出由中等和较强等级向强等级转变的态势。

林芝市土地利用程度具有较显著的全域空间正相关性，有明显的集聚特征，且集聚的趋势平稳。

局域空间自相关特征为市域东北部和西北部为低—低集聚区，市域低海拔和雅鲁藏布江、尼洋河宽河谷地区为高—高集聚区。

土地利用程度变化(增加)空间上整体集聚，局域特征以尼洋河宽河谷流域、其他人口活动密集地区和低海拔区为高—高集聚，而低—高集聚区分布在高—高集聚区外围。

【结论】林芝市土地利用程度逐渐增强，土地利用程度及其变化都不是随机分布的，土地利用程度变化高的区域彼此相邻，土地利用程度变化低的区域彼此相邻。

关键词: 土地利用程度；空间异质性；空间自相关；林芝市
中图分类号: U412.14 文献标志码: A 文章编号: 1004–390X (2023) 03−0511−09
Study of Land Use Intensity Change and Its Spatial
Heterogeneity of Nyingchi City in Tibet
GUO Jianbin 1,2，ZHANG Ying 1，ZHANG Zhiwei 1，HOU Lei 1，ZENG Weilian 2,3
(1. Resources & Environment College, Tibet Agriculture & Animal Husbandry University, Nyingchi 860000, China;2. Research Academy for Rural Revitalization in Tibet, Tibet Agriculture & Animal Husbandry University, Nyingchi
860000, China; 3. Plant Sciences College, Tibet Agriculture & Animal Husbandry
University, Nyingchi 860000, China)
Abstract: ［Purpose ］To reveal the spatial heterogeneity characteristics of land use intensity of alpine gorge in Tibetan Plateau, providing a theoretical basis for local territorial spatial planning.［Method ］Based on the Landsat satellite images in 2010 and 2020, the data of land use in 30 meters grid were obtained, and the spatial heterogeneity characteristics of land use intensity and its change in
云南农业大学学报（自然科学），2023，38(3)：511−519
Journal of Yunnan Agricultural University (Natural Science)
E-mail: ********************
收稿日期：2022-10-04 修回日期：2023-05-11 网络首发日期：2023-07-07
*基金项目：西藏自治区自然科学基金项目(XZ202001ZR0049G)；西藏自治区高等学校人文社会科学研究
项目(SK2021-46)。

作者简介：郭健斌(1980—)，男，内蒙古清水河人，硕士，副教授，主要从事西藏土地资源可持续利用研
究。

E-mail ：******************.cn
网络首发地址：https:///kcms2/detail/53.1044.S.20230705.1446.002.html
Nyingchi were analyzed by spatial autocorrelation analysis method. ［Results］The grid units of stronger and strongest level of land use intensity were mainly distributed in the low altitude areas and the wide valleys of the Yarlung Tsangpo River and the Nyang River, the grid units of weakest and weaker level of land use intensity were mainly distributed in the high altitude area. From 2010 to 2020, the increased area of strongest level came from the reduced area of middle and stronger level of land use intensity. The land use intensity in Nyingchi had relatively strong positive global spatial auto-correlation with Moran’s Ⅰ, had obvious spatial aggregation law. The graph of local autocorrelation analysis showed that land use intensity low-low typology was in the northeast and northwest in Ny-ingchi, high-high typology was in the low altitude areas and the wide valleys of the Yarlung Tsangpo River and the Nyang River. The land use intensity change also had positive global spatial autocorrela-tion in Nyingchi from 2010 to 2020, high-high typology was in the wide valleys of the Nyang River, other densely inhabited districts and low altitude areas, but the low-high typology distributed around the high-high typology areas. ［Conclusion］The land use intensity gradually strengthens, the land use intensity and its changes in Nyingchi is not disorderly, the regions with high land use change are adjacent to each other, and the regions with low land use change are adjacent to each other.
Keywords: land use intensity; spatial heterogeneity; spatial autocorrelation; Nyingchi City
土地利用变化研究是近年来全球变化研究的热点[1]，主要包括土地覆被变化和土地利用程度变化2个方面。

但当前学术界普遍关注土地覆被变化的空间形态、类型与数量、演化机理和驱动力研究[2-4]，而对土地利用程度的关注度不高。

土地利用程度是一个地区土地开发数量和程度的度量指标，反映了土地自然属性和人类活动对土地利用的综合影响[5]。

探究区域土地利用程度的变化，有利于人类可持续利用土地资源。

国外关于土地利用程度的研究侧重于农用地的投入产出分析[6]；国内有学者选用1个[5]或多个指标[7]对土地利用程度进行量化研究，为了便于区域之间的横向比较，多数学者采用[8]或改进[9]楚玉山等[10]提出的土地利用程度综合指数量化分析方法对中国不同自然地理区域[11-12]和行政区域[8, 13]的土地利用程度进行综合评价。

随着遥感和地理信息系统技术的快速发展，多源高精度影像数据在土地利用程度变化的研究中也得到普遍运用[13-15]。

空间研究中经常使用的方法是空间自相关分析[16-19]。

由于土地利用变化是一类具有高度空间异质性的地理现象，利用空间自相关技术分析其空间聚集规律的同时，还可以将空间数据可视化。

空间自相关性使用全域和局域2种指标，Moran’s Ⅰ指数在判断区域土地利用程度的空间聚集上具有可靠性[20]和不易受偏离正态分布影响的特点[19]，得到广泛应用[11-12, 21]。

在中国，土地利用程度的研究多聚焦于黄土高原[9, 14, 22]、西南喀斯特山区[12]和西北内陆干旱区[23-24]等生态脆弱地区以及中东部人口密集地区[8, 11, 21]。

青藏高原作为中国典型的生态环境脆弱地区和国家重要生态安全屏障，目前仅见三江源地区土地利用程度变化的时空分析[25]和《西藏自治区土地利用》对西藏各县土地利用程度的对比分析[10]。

林芝市位于西藏东南部的高山峡谷地区，垂直落差大，土质疏松，局部暴雨多，滑坡和泥石流等自然灾害频繁。

随着人口和人类社会经济活动的增加，不但会改变林芝市的社会经济发展模式，也会威胁脆弱的生态环境。

因此，客观和准确地认识林芝市土地利用强度及其空间分异对高山峡谷地区的可持续发展具有重要意义。

本研究以藏东南林芝市为研究区域，以2010年和2020年遥感影像为基础，建立土地利用强度及其变化模型，采用Moran’s Ⅰ指数对林芝市土地利用程度进行空间异质性分析，探究区域土地利用程度的内部差异和空间相关性，可为当地国土空间规划和生态环境保护决策提供理论基础。

1 材料与方法
1.1 研究区概况
林芝市是西藏自治区下辖地级市，位于
512云南农业大学学报第 38 卷
E92°09′~98°47′，N26°52′~30°40′ (图1)。

东与昌都市和云南省迪庆藏族自治州毗邻，西与拉萨市和山南市交界，北与那曲市相连，南与缅甸和印度接壤。

东西长646.7 km ，南北宽353.2 km ，平均海拔3 100 m ，总面积11.7万km 2 (包括印度非
法控制区)[26]。

林芝市地处青藏高原南缘，地势自西北向东南倾斜，雅鲁藏布江和怒江等水系切割高原，形成典型的高山峡谷地貌，海拔落差大。

市域气候分为热带北缘和亚热带山地湿润气候区以及高原温带半湿润气候区，立体气候明显[27]。

林芝史称贡布，是藏民族发祥地之一。

和平解放后，1960—1963年曾设林芝专区，随后撤销，1986年林芝地区行政公署恢复成立，2015年3月国务院批复同意撤销林芝地区，设立林芝市。

2020年全市下辖1区6县，包括54个乡(镇)、2个街道办，共504个村(居)，总人口23.89万，是以藏族为主体，汉、门巴和珞巴等35个民族和僜人多民族聚居区。

境内分布有中国最大的原始林区，耕地散布于河谷地带，草地位于林线上下，垂直地带性明显。

2020年全市实现地区生产总值191.34亿元，可比增长7.9%，三次产业的比例为6.4∶38.7∶54.9[28]。

1.2 数据来源
土地利用数据来源于中国科学院计算机网络
信息中心创立的地学大数据平台——地理空间数据云(https:///)的2010年Landsat-7ETM+和2020年Landsat-8OLI 2期各11景卫星影像(空间分辨率为30 m×30 m)，为确保云量小于10%，影像日期选择当年10—12月或次年1—3月。

利用ENVI 5.1对2期遥感数据进行辐射定标、大气校正和裁剪等预处理，对于部分ETM+数据条带缺失的问题，利用ENVI 软件中的条带修复补丁，并基于单文件三角插值的条带填充进行修复。

结合研究区实际情况，依据《土地利用现状分类标准》(GB/T 21010—2017)，选取典型林地、水域、耕地、建设用地、草地和未利用地为训练样本，对研究区典型地物进行最大似然法监督分类。

同时借助Google Earth 高分辨
100
200 km
50国界 national boundary 省界 province boundary 市界 city boundary
N
林芝市在西藏的位置
location of Nyingchi City in Tibet Naqu City
Lhasa City Shannan City
India
Myanmar
Yunnan Province
Changdu City
高程/m elevation
7 261
101
印度缅甸云南省
昌都市山南市拉萨市那曲市
注：林芝市辖区范围来源于西藏自治区自然资源厅公布的全区系列标准地图，审图号：藏S(2022)004号。

Note: The base map of Nyingchi comes from a series of standard maps of the whole area published by the Department of Natural Resources of Tibetan Autonomous Region, the approval number is Tibet S(2022)004.
图 1 研究区位置
Fig. 1 Location of the study area
第 3 期郭健斌，等：西藏林芝市土地利用程度变化及其空间异质性分析513
率影像和部分实地考察数据，每一地类随机选取100个样本进行精度评价，2010年和2020年2期土地利用类型识别精度分别为90%和92%，kappa 系数分别为0.88和0.90，满足研究需求。

1.3 研究方法
1.3.1 土地利用程度指数
土地利用程度指数是以土地开发利用程度来反映人类活动对某一区域土地利用变化的影响水平[14]。

本研究参考庄大方等[29]数量化土地利用程度分析方法，将土地利用程度按照土地自然综合体在社会因素影响下的自然平衡保持状态分为4级，并分级赋予指数，从而给出土地利用程度的定量表达。

其计算公式为
LA A i i C i i i n 式中：指某一区域土地利用综合指数，为
100~400的连续函数，值越大表明该区域土地利用程度越高；为第级土地利用程度的分级指数，参照已有研究[8, 10-11, 29]中的分级标准并结合研究区实际，将未利用地赋值为1，林地、草地和水域赋值为2，耕地赋值为3，建设用地赋值为4；为第级土地利用程度分级面积百分比，即第级土地利用程度分级类型与区域总面积的比值；为土地利用程度分级数，本研究取值为4。

由于全局空间自相关大小对采样方式和采样尺度较敏感[30]，介于研究区域林草地和未利用地占比大，垂直地带性明显，综合考虑评价单元的实效性和工作量，最终采用10 km×10 km 的网格对研究区2期土地利用类型图进行采样，共计获得1 314个采样网格进行土地利用程度指数计算及土地利用程度变化分析；再利用ArcGIS 10.6，在自然断点基础上结合当地实际，依据土地利用程度指数将林芝市土地利用程度划分为强、较强、中等、较弱和弱5级(表
1)。

1.3.2 土地利用程度变化
本研究利用2个时期土地利用程度指数的变
化来衡量某一区域土地利用程度的变化趋势[23]，CLA b −a LA b LA a 式中：为该区域研究期初到期末土地利用程度变化指数；和分别为该区域研究期末和期初土地利用程度指数。

1.3.3 空间自相关分析
空间自相关是指某一变量在一定区域内空间上的相关性，土地利用程度(变化)作为连续空间的度量指标也具有空间相关性。

本研究选用全域和局域Moran ’s Ⅰ指数衡量土地利用程度(变化)的空间聚集程度。

全域Moran’s Ⅰ统计量用于描述土地利用程度(变化)在整个研究区域的空间依赖程度，可反映土地利用程度(变化)是否在空间上聚类，其表达式为
n X i X j i j X W i j W i j
=1i j W i j =0i j I 式中：为样本数，和 分别为区域内第和第处的土地利用程度(变化)；为区域内所有样本
的平均值；为空间权重矩阵，本研究采用基于距离标准的二进制空间权重矩阵，当时，区域与相邻，当时，区域与不相邻。

值的显著性用以下公式检验[31]
：
E (I )v (I )式中：和分别为Moran’s Ⅰ的数学期望和方差。

全域Moran’s Ⅰ指数范围在[−1,1]，在既定的显著性水平下，I >0时，表示土地利用程度(变化)之间存在显著正相关，值越大相关性越强，空间上聚合； I <0时，表示土地利用程度(变化)之间存在显著负相关，绝对值越大相关性越强，空间上离散； I =0时，表示土地利用程度(变化)不存在空间自相关，随机分布。

虽然全域空间自相关可以反映土地利用程度是否在空间集聚，但它不能进一步显示具体的空间集聚特征[19]。

局域空间自相关则可以进一步展
表 1 土地利用程度分级表
Tab. 1 Grades of land use intensity classification
土地利用程度等级grades of land use intensity 土地利用程度指数land use intensity index
弱 weakest ≤120.00
较弱 weaker 120.01~145.00中等 middle 145.01~180.00较强 stronger 180.01~200.00
强 strongest
≥200.01
514云南农业大学学报第 38 卷
示土地利用程度空间集聚或离散的位置和范围，本研究选用的局部空间关联指标(local indicators spatial association，LISA)为局域Moran’s Ⅰ指数。

局域Moran’s Ⅰ指数用于描述土地利用程度在局域空间是否存在空间自相关，能更直观地反映研究区土地利用程度的空间集聚状况，其表达式为
式中变量的含义与全域Moran’s Ⅰ统计量公式相同。

局域Moran’s Ⅰ指数表示土地利用程度(变化)在不同空间位置上的集聚模式，当I >0时，表示以高—高值或低—低值的模式集聚；当I<0时， 表示以高—低值或低—高值的模式集聚。

2 结果与分析
2.1 土地利用结构及利用程度的时空特征
2.1.1 土地利用结构变化
2010—2020年，林芝市的草地、林地、未利用地和水域4类用地稍有减少，其中，林地的占比从2010年的53.49%降至2020年的53.43%，未利用地的占比从2010年的16.34%降至2020年的16.32%，草地的占比从26.20%降至26.17%，水域的占比从2.86%降至2.85%。

2010—2020年，林芝市的耕地占比从2010年的1.10%增加至2020年的1.17%；建设用地明显增加，从2010的25.05 km2增加至2020年的66.57 km2，增长了2.66倍，占比从0.02%增至0.06% (表2)。

2.1.2 土地利用程度时间变化
由表3可知：2010—2020年，弱等级网格单元面积的占比保持在0.51%；较弱等级网格单元面积的占比从6.50%增加至6.55%；中等等级网格单元面积的占比从26.02%降至25.82%；较强等级网格单元面积的占比从 48.48% 降至 45.61%；强等级网格单元面积的占比从 18.48% 增至 21.51%。

整体而言，土地利用程度呈现出由中等和较强等级向强等级转变的态势。

中等和较强等级土地利用程度网格单元面积从85 486.70 km2降至81 959.90 km2；除弱和较弱等级土地利用程度网格单元面积增加1.07和58.09 km2外，其余3 467.60 km2 (占比98.33%)全部转变为强等级土地利用程度。

2.1.3 土地利用程度空间特征
林芝市土地利用程度空间分布(图2)表明：较强和强等级土地利用程度网格单元主要分布在墨脱县、察隅县南部的低海拔区域和人类活动密集的雅鲁藏布江及其支流尼洋河宽河谷流域，弱和较弱等级土地利用程度网格单元主要分布在高海拔的市域东北部、西北部和中部东喜马拉雅南迦巴瓦区域。

2010年林芝市土地利用程度指数平均值为181.27，2020年为181.37。

2010—2020年，强等级土地利用程度网格单元在尼洋河流域和墨脱县有明显增加。

2.2 土地利用程度指数空间自相关分析
2.2.1 全域空间自相关
z
α
由图3可知：2010年和2020年的Moran’s Ⅰ指数均为正数，数值分别为0.588和0.587，大于0、接近1；值分别为31.21和26.32，均大于临界值(2.58)，并且通过了=0.01的显著性检
表 2 2010—2020年林芝市各用地类型面积及其占比Tab. 2 Area and proportion of land use types in
Nyingchi City from 2010 to 2020
类型land use types
20102020
面积/km2
area
比例/%
proportion
面积/km2
area
比例/%
proportion
草地 grassland29 898.9826.2029 867.3926.17耕地 cropland1 257.37 1.101 332.12 1.17建设用地 construction land25.050.0266.570.06林地 forest61 042.8353.4960 979.7553.43水域 water3 258.32 2.863 253.11 2.85未利用地 unused land18 643.2816.3418 626.8916.32
表 3 2010—2020年林芝市土地利用程度
等级面积及其占比
Tab. 3 Area and proportion of land use intensity in
Nyingchi City from 2010 to 2020
土地利用程度等级
grades of land use intensity
20102020
面积/km2
area
比例/%
proportion
面积/km2
area
比例/%
proportion 弱 weakest588.530.51589.600.51较弱 weaker7 456.36 6.507 514.45 6.55中等 middle29 860.5026.0229 622.5025.82较强 stronger55 626.2048.4852 337.4045.61强 strongest21 208.0018.4824 675.6021.51
第 3 期郭健斌，等：西藏林芝市土地利用程度变化及其空间异质性分析515
验，表明林芝市土地利用程度存在较强的正空间相关性，即土地利用程度空间分布上不是随机的，而是表现为高值或低值的集聚。

可见，林芝市土地利用程度的空间联系特征为：土地利用程度高的网格单元相对趋向于与土地利用程度高的网格单元相邻，土地利用程度低的网格单元相对趋向于与土地利用程度低的网格单元相邻。

2020年Moran’s Ⅰ指数较2010年仅减少了0.001，说明2010—2020年土地利用程度的空间分布呈现集聚的趋势基本稳定。

2.2.2 局域空间自相关
由图4可知：低—低值集聚区主要分布在市域东北部和西北部高海拔地区；高—高值集聚区分布在市域低海拔的南部和东南部部分地区、尼洋河流域和雅鲁藏布江宽河谷地区。

2020年尼洋
河流域高—高值集聚区较2010年明显增加。

2.3 土地利用程度变化空间自相关分析
z α2010—2020年，林芝市土地利用程度变化全
域Moran’s Ⅰ指数为0.223，值为22.21，大于临界值，通过了=0.01的显著性检验，表明其土地利用程度变化量存在较强的空间正相关性，即土地利用程度变化剧烈的网格单元相对趋向于与土地利用程度变化剧烈的网格单元相邻，土地利用程度变化微弱的网格单元相对趋向于与土地利用程度变化微弱的网格单元相邻(图5)。

由局域空间自相关LISA 集聚地图(图6)可知：2010—2020年林芝市土地利用程度变化高—高值集聚区集中分布在尼洋河流域的巴宜区的八一镇、百巴镇和更张门巴民族乡以及工布江达县的巴河镇、错高乡和工布江达镇，在巴宜区的鲁
弱 weakest weaker middle stronger
strongest
图 2 林芝市土地利用程度等级
Fig. 2 Grades of land use intensity in Nyingchi City
莫兰指数 Moran’s Ⅰ: 0.588
空间滞后向量s p a c e l a g v e c t o r
标准化后的土地利用程度指数standardized index of land use intensity
−5−4
−3−2−1012345−8
−6−4−202468
莫兰指数 Moran’s Ⅰ: 0.587
图 3 2010年(左)和2020年(右)林芝市土地利用程度散点图
Fig. 3 Scatterplots of land use intensity in Nyingchi City in 2010 (left) and 2020 (right)
516
云南农业大学学报第 38 卷
朗镇、波密县的扎木镇、墨脱县的墨脱镇和格当乡以及察隅县的下察隅镇和上察隅镇也有分布；低—高值集聚区分布在上述高—高值集聚范围的
周边；低—低值集聚区分布在墨脱县的墨脱镇和背崩乡、米林县的丹娘乡和羌纳乡以及巴宜区的米瑞乡、林芝镇和布久乡，在高海拔区域如工布江达县朱拉乡和察隅县古玉乡等也有分布；高—低值在波密县松宗镇和市域边缘有少量分布。

3 讨论
2010—2020年，林芝市人口从19.51万增加
到23.89万，固定资产投资从65.17亿元增加到125.12亿元[28]，人口和资金集中在海拔低、地势
平坦的河谷地带，导致这一区域建设用地和耕地增加，草地、未利用地和林地减少。

建设用地的增长反映了林芝市非农经济活动的快速增长；耕地的增减变化说明林芝市人地关系渐趋紧张，城市规模扩张态势明显。

在全球变暖背景下，青藏高原内流区得益于冰川融雪和降水增加，湖泊水域扩张明显[32-33]；而位于藏东南外流区的林芝市冰川消融[33]，水域面积减少。

杨玉婷等[23]研究表明张掖市甘州区2005年的土地利用程度指数为255.25，而林芝市2010年和2020年的土地利用程度指数分别为181.27和181.37，远低于西北干旱内陆河流域，土地利用
程度低。

林芝市土地利用程度空间分布明显受地形条件和社会经济发展的影响，高海拔地区难利用土地占比高，土地利用程度低；低海拔区域林草地占比高，土地利用程度较高；平坦的河谷地带耕地占比高，村庄和城镇分布广泛，土地利用程度最高。

土地利用程度综合指数能够诠释林芝市土地利用程度的空间分布与自然环境间的关系，但尚不能准确厘定其作用机制，需要进一步
non-significant (859)
—高 high-high (270)
non-significant (867)—高 high-high (267)0100200 km
50
图 4 2010年(左)和2020(右)年林芝市土地利用程度LISA 集聚图
Fig. 4 LISA agglomeration graphs of land use intensity in Nyingchi City in 2010 (left) and 2020 (right)
−12−7−23
8
13
−30
−20−100102030
空间滞后向量s p a c e l a g v e c t o r
标准化后的土地利用程度指数变化standardized index of the change of land use
莫兰指数 Moran’s I : 0.223
图 5 2010—2020年林芝市土地利用程度变化散点图Fig. 5 Scatterplots of land use intensity change in Nyingchi
City from 2010 to 2020
100200 km
50图 6 2010—2020年林芝市土地利用程度
变化LISA 集聚图Fig. 6 LISA agglomeration graph of land use intensity
change in Nyingchi City from 2010 to 2020
第 3 期
郭健斌，等：西藏林芝市土地利用程度变化及其空间异质性分析
517
研究。

2010—2020年，林芝市土地利用程度由中等和较强等级向强等级转变，说明随着当地人口数量和社会经济活动的不断增加，除人类难以涉足的区域外，林芝市土地利用程度都在不断增强。

林芝市域西北部和东北部为高海拔区域，是土地利用程度低—低聚集区域；而南部和东南部为低海拔区域，是土地利用程度高—高集聚区域；以巴宜市区为中心的尼洋河和雅鲁藏布江宽河谷区域是土地利用程度高—高集聚区域。

这些LISA 集聚特征与实际情况较吻合，能够准确地反映林芝市自然地理环境条件和社会经济发展的特点。

2010—2020年，土地利用程度变化(增加)在低海拔、宽河谷和人类活动密集地区呈高—高集聚的特征，其外围呈低—高集聚的过渡特征，说明随着林芝市城镇化水平的不断提高，农村地区土地利用程度在不断地缓释，而城镇用地利用程度不断增强，且有持续向外扩张的趋势。

在可利用土地资源极其稀缺的藏东南高山峡谷地区，未来土地利用应该以挖掘现有建设用地潜力，提高城镇用地的集约性，引导其有序扩张。

4 结论
本研究利用Landsat 影像解译土地利用数
据，计算土地利用程度及其变化，采用空间自相关法分析土地利用程度及其变化的空间特性，结果表明：林芝市人类社会经济活动密集的低海拔河谷平地土地利用程度高；人类社会经济活动稀少的高海拔区域土地利用程度低。

2010—2020年该区域土地利用程度不断增强，低海拔河谷平地的强等级土地利用程度网格单元增加更为明显。

市政府驻地尼洋河流域是土地利用程度变化的高—高值集聚区，即土地利用程度变化的热点区域，今后国土空间规划和生态环境保护应着重关注这一区域的土地集约利用。

[ 参考文献 ]
Global Land Project. Science plan and impleme-
ntation strategy[R]. Stockholm: International Geosphere-Biosphere Program Secretariat, 2005.
[1]匡文慧, 张树文, 杜国明, 等. 2015—2020年中国土地利
用变化的时空遥感制图及时空特征分析[J]. 地理学报,2022, 77(5): 1056. DOI: 10.11821/dlxb202205002.
[2]杨丽萍, 张静, 贡恩军, 等. GEE 联合多源数据的西安
市土地利用时空格局及驱动力分析[J]. 农业工程学报,
[3]2022, 38(2): 279. DOI: 10.11975/j.issn.1002-6819.2022.02.031.
吴婷婷, 刘学录. 甘肃省庄浪县土地利用结构动态变化研究[J]. 云南农业大学学报(自然科学), 2015, 30(1):112. DOI: 10.16211/j.issn.1004-390X(n).2015.01.019.[4]
HUANG X, HUANG X J, LIU M M, et al. Spatial-tem-poral dynamics and driving forces of land development intensity in the western China from 2000 to 2015[J]. Ch-inese Geographical Science, 2020, 30(1): 16. DOI: 10.1007/s11769-020-1095-2.
[5]
ERB K H, HABERL H, JEPSEN M R, et al. A conceptu-al framework for analysing and measuring land-use in-tensity[J]. Current Opinion in Environmental Sustainabil-ity, 2013, 5(5): 464. DOI: 10.1016/j.cosust.2013.07.010.[6]
曹瑞芬, 蔡银莺. 基于AHP 的土地开发利用程度评价及分析: 以武汉市为例[J]. 华中农业大学学报(社会科学版), 2011(1): 65. DOI: 10.13300/ki.hnwkxb.2011.01.020.
[7]
梁明, 聂拼, 陆胤昊, 等. 淮南市土地利用程度变化过程的时空演化特征[J]. 农业工程学报, 2019, 35(22): 99.DOI: 10.11975/j.issn.1002-6819.2019.22.011.
[8]
吴金华, 李纪伟, 梁晶晶. 土地利用程度与效益关系研究: 以延安市为例[J]. 中国土地科学, 2011, 25(8): 54.DOI: 10.13708/11-2640.2011.08.010.
[9]
楚玉山, 刘纪远. 西藏土地利用[M]. 北京: 科学出版社,1992.
[10]黄木易, 何翔, 吴迪, 等. 巢湖流域土地利用程度变化及其空间异质性分析[J]. 土壤, 2015, 47(5): 994. DOI: 10.13758/ki.tr.2015.05.029.
[11]
张颖, 赵宇鸾. 黔桂岩溶山区土地利用程度演变的空间分异特征[J]. 水土保持研究, 2018, 25(1): 287. DOI: 10.13869/ki.rswc.20171225.001.
[12]
于明雪, 孙建国, 杨维涛, 等. 基于贝叶斯层次时空模型的甘肃省土地利用程度演变分析[J]. 地理科学, 2022,42(5): 918. DOI: 10.13249/ki.sgs.2022.05.017.
[13]
吴琳娜, 杨胜天, 刘晓燕, 等. 1976年以来北洛河流域土地利用变化对人类活动程度的响应[J]. 地理学报,2014, 69(1): 54. DOI: 10.11821/dlxb201401005.
[14]
黄华富, 戴文远, 苏木兰. 海岛生态脆弱区土地利用程度空间格局演化: 以福建省海坛岛为例[J]. 福建师范大学学报(自然科学版), 2016, 32(2): 92.
[15]
HUANG W, WANG M Y, CHEN X Y. Analysis on spa-tial pattern evolution of cultivated land in urban area based on spatial autocorrelation analysis: a case study of Luoyang City[M]//SUN X M, ZHANG X R, XIA Z H, et al. Artificial intelligence and security: part 1. Dublin:Springer, 2021.
[16]
WANG H, ZHU Y X, WANG M Y, et, al. Research of spatial pattern for cultivated land quality in Henan Province based on spatial autocorrelation[M]//SUN X M,ZHANG X R, XIA Z H, et al. Artificial intelligence and security: part 1. Dublin: Springer, 2021.
[17]
518云南农业大学学报第 38 卷。