西藏林芝市土地利用程度变化及其空间异质性分析
Study of Land Use Intensity Change and Its Spatial Heterogeneity of Nyingchi City in Tibet
-
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 研究区概况
林芝市是西藏自治区下辖地级市,位于E92°09′~98°47′,N26°52′~30°40′ (图1)。东与昌都市和云南省迪庆藏族自治州毗邻,西与拉萨市和山南市交界,北与那曲市相连,南与缅甸和印度接壤。东西长646.7 km,南北宽353.2 km,平均海拔3 100 m,总面积11.7万km2 (包括印度非法控制区)[26]。林芝市地处青藏高原南缘,地势自西北向东南倾斜,雅鲁藏布江和怒江等水系切割高原,形成典型的高山峡谷地貌,海拔落差大。市域气候分为热带北缘和亚热带山地湿润气候区以及高原温带半湿润气候区,立体气候明显[27]。
图 1 研究区位置注:林芝市辖区范围来源于西藏自治区自然资源厅公布的全区系列标准地图,审图号:藏S(2022)004号。Figure 1. Location of the study areaNote: 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.林芝史称贡布,是藏民族发祥地之一。和平解放后,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://www.gscloud.cn/)的2010年Landsat-7 ETM+和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个样本进行精度评价,2010年和2020年2期土地利用类型识别精度分别为90%和92%,kappa系数分别为0.88和0.90,满足研究需求。
1.3 研究方法
1.3.1 土地利用程度指数
土地利用程度指数是以土地开发利用程度来反映人类活动对某一区域土地利用变化的影响水平[14]。本研究参考庄大方等[29]数量化土地利用程度分析方法,将土地利用程度按照土地自然综合体在社会因素影响下的自然平衡保持状态分为4级,并分级赋予指数,从而给出土地利用程度的定量表达。其计算公式为
$ \mathrm{L}\mathrm{A}=100\times \left(\sum _{i=1}^{n}{A}_{i}\times {C}_{i}\right) 。 $
This page contains the following errors:
error on line 1 at column 1: Start tag expected, '<' not foundBelow is a rendering of the page up to the first error.
由于全局空间自相关大小对采样方式和采样尺度较敏感[30],介于研究区域林草地和未利用地占比大,垂直地带性明显,综合考虑评价单元的实效性和工作量,最终采用10 km×10 km的网格对研究区2期土地利用类型图进行采样,共计获得1 314个采样网格进行土地利用程度指数计算及土地利用程度变化分析;再利用ArcGIS 10.6,在自然断点基础上结合当地实际,依据土地利用程度指数将林芝市土地利用程度划分为强、较强、中等、较弱和弱5级(表1)。
表 1 土地利用程度分级表Table 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 1.3.2 土地利用程度变化
本研究利用2个时期土地利用程度指数的变化来衡量某一区域土地利用程度的变化趋势[23],其表达式为
$ \mathrm{C}{\mathrm{L}\mathrm{A}}_{b-a}={\mathrm{L}\mathrm{A}}_{b}-{\mathrm{L}\mathrm{A}}_{a} 。 $
This page contains the following errors:
error on line 1 at column 1: Start tag expected, '<' not foundBelow is a rendering of the page up to the first error.
1.3.3 空间自相关分析
空间自相关是指某一变量在一定区域内空间上的相关性,土地利用程度(变化)作为连续空间的度量指标也具有空间相关性。本研究选用全域和局域Moran’s Ⅰ指数衡量土地利用程度(变化)的空间聚集程度。
全域Moran’s Ⅰ统计量用于描述土地利用程度(变化)在整个研究区域的空间依赖程度,可反映土地利用程度(变化)是否在空间上聚类,其表达式为
$ I=\frac{n\displaystyle\sum _{i=1}^{n}\sum _{j=1}^{n}{W}_{ij}\left({X}_{i}-\overline{X}\right)\left({X}_{j}-\overline{X}\right)}{\displaystyle\sum _{i=1}^{n}\sum _{j=1}^{n}{W}_{ij}\sum _{i=1}^{n}{\left({X}_{i}-\overline{X}\right)}^{2}} 。 $
This page contains the following errors:
error on line 1 at column 1: Start tag expected, '<' not foundBelow is a rendering of the page up to the first error.
$ z=\frac{1-E\left(I\right)}{\sqrt{v\left(I\right)}} 。 $
This page contains the following errors:
error on line 1 at column 1: Start tag expected, '<' not foundBelow is a rendering of the page up to the first error.
全域Moran’s Ⅰ指数范围在[−1,1],在既定的显著性水平下,I >0时,表示土地利用程度(变化)之间存在显著正相关,值越大相关性越强,空间上聚合; I<0时,表示土地利用程度(变化)之间存在显著负相关,绝对值越大相关性越强,空间上离散; I=0时,表示土地利用程度(变化)不存在空间自相关,随机分布。
虽然全域空间自相关可以反映土地利用程度是否在空间集聚,但它不能进一步显示具体的空间集聚特征[19]。局域空间自相关则可以进一步展示土地利用程度空间集聚或离散的位置和范围,本研究选用的局部空间关联指标(local indicators spatial association,LISA)为局域Moran’s Ⅰ指数。
局域Moran’s Ⅰ指数用于描述土地利用程度在局域空间是否存在空间自相关,能更直观地反映研究区土地利用程度的空间集聚状况,其表达式为
$ I=\frac{\displaystyle\sum _{j=1,j\ne i}^{n}{W}_{ij}\left({X}_{i}-\overline{X}\right)\left({X}_{j}-\overline{X}\right)}{\displaystyle\sum _{i=1}^{n}{\left({X}_{i}-\overline{X}\right)}^{2}} 。 $
式中变量的含义与全域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 2010—2020年林芝市各用地类型面积及其占比Table 2. Area and proportion of land use types in Nyingchi City from 2010 to 2020类型
land use types2010 2020 面积/km2
area比例/%
proportion面积/km2
area比例/%
proportion草地 grassland 29 898.98 26.20 29 867.39 26.17 耕地 cropland 1 257.37 1.10 1 332.12 1.17 建设用地 construction land 25.05 0.02 66.57 0.06 林地 forest 61 042.83 53.49 60 979.75 53.43 水域 water 3 258.32 2.86 3 253.11 2.85 未利用地 unused land 18 643.28 16.34 18 626.89 16.32 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%)全部转变为强等级土地利用程度。
表 3 2010—2020年林芝市土地利用程度等级面积及其占比Table 3. Area and proportion of land use intensity in Nyingchi City from 2010 to 2020土地利用程度等级
grades of land use intensity2010 2020 面积/km2
area比例/%
proportion面积/km2
area比例/%
proportion弱 weakest 588.53 0.51 589.60 0.51 较弱 weaker 7 456.36 6.50 7 514.45 6.55 中等 middle 29 860.50 26.02 29 622.50 25.82 较强 stronger 55 626.20 48.48 52 337.40 45.61 强 strongest 21 208.00 18.48 24 675.60 21.51 2.1.3 土地利用程度空间特征
林芝市土地利用程度空间分布(图2)表明:较强和强等级土地利用程度网格单元主要分布在墨脱县、察隅县南部的低海拔区域和人类活动密集的雅鲁藏布江及其支流尼洋河宽河谷流域,弱和较弱等级土地利用程度网格单元主要分布在高海拔的市域东北部、西北部和中部东喜马拉雅南迦巴瓦区域。2010年林芝市土地利用程度指数平均值为181.27,2020年为181.37。2010—2020年,强等级土地利用程度网格单元在尼洋河流域和墨脱县有明显增加。
2.2 土地利用程度指数空间自相关分析
2.2.1 全域空间自相关
This page contains the following errors:
error on line 1 at column 1: Start tag expected, '<' not foundBelow is a rendering of the page up to the first error.
2.2.2 局域空间自相关
由图4可知:低—低值集聚区主要分布在市域东北部和西北部高海拔地区;高—高值集聚区分布在市域低海拔的南部和东南部部分地区、尼洋河流域和雅鲁藏布江宽河谷地区。2020年尼洋河流域高—高值集聚区较2010年明显增加。
2.3 土地利用程度变化空间自相关分析
This page contains the following errors:
error on line 1 at column 1: Start tag expected, '<' not foundBelow is a rendering of the page up to the first error.
由局域空间自相关LISA集聚地图(图6)可知:2010—2020年林芝市土地利用程度变化高—高值集聚区集中分布在尼洋河流域的巴宜区的八一镇、百巴镇和更张门巴民族乡以及工布江达县的巴河镇、错高乡和工布江达镇,在巴宜区的鲁朗镇、波密县的扎木镇、墨脱县的墨脱镇和格当乡以及察隅县的下察隅镇和上察隅镇也有分布;低—高值集聚区分布在上述高—高值集聚范围的周边;低—低值集聚区分布在墨脱县的墨脱镇和背崩乡、米林县的丹娘乡和羌纳乡以及巴宜区的米瑞乡、林芝镇和布久乡,在高海拔区域如工布江达县朱拉乡和察隅县古玉乡等也有分布;高—低值在波密县松宗镇和市域边缘有少量分布。
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,远低于西北干旱内陆河流域,土地利用程度低。林芝市土地利用程度空间分布明显受地形条件和社会经济发展的影响,高海拔地区难利用土地占比高,土地利用程度低;低海拔区域林草地占比高,土地利用程度较高;平坦的河谷地带耕地占比高,村庄和城镇分布广泛,土地利用程度最高。土地利用程度综合指数能够诠释林芝市土地利用程度的空间分布与自然环境间的关系,但尚不能准确厘定其作用机制,需要进一步研究。2010—2020年,林芝市土地利用程度由中等和较强等级向强等级转变,说明随着当地人口数量和社会经济活动的不断增加,除人类难以涉足的区域外,林芝市土地利用程度都在不断增强。
林芝市域西北部和东北部为高海拔区域,是土地利用程度低—低聚集区域;而南部和东南部为低海拔区域,是土地利用程度高—高集聚区域;以巴宜市区为中心的尼洋河和雅鲁藏布江宽河谷区域是土地利用程度高—高集聚区域。这些LISA集聚特征与实际情况较吻合,能够准确地反映林芝市自然地理环境条件和社会经济发展的特点。2010—2020年,土地利用程度变化(增加)在低海拔、宽河谷和人类活动密集地区呈高—高集聚的特征,其外围呈低—高集聚的过渡特征,说明随着林芝市城镇化水平的不断提高,农村地区土地利用程度在不断地缓释,而城镇用地利用程度不断增强,且有持续向外扩张的趋势。在可利用土地资源极其稀缺的藏东南高山峡谷地区,未来土地利用应该以挖掘现有建设用地潜力,提高城镇用地的集约性,引导其有序扩张。
4. 结论
本研究利用Landsat影像解译土地利用数据,计算土地利用程度及其变化,采用空间自相关法分析土地利用程度及其变化的空间特性,结果表明:林芝市人类社会经济活动密集的低海拔河谷平地土地利用程度高;人类社会经济活动稀少的高海拔区域土地利用程度低。2010—2020年该区域土地利用程度不断增强,低海拔河谷平地的强等级土地利用程度网格单元增加更为明显。市政府驻地尼洋河流域是土地利用程度变化的高—高值集聚区,即土地利用程度变化的热点区域,今后国土空间规划和生态环境保护应着重关注这一区域的土地集约利用。
-
图 1 研究区位置
注:林芝市辖区范围来源于西藏自治区自然资源厅公布的全区系列标准地图,审图号:藏S(2022)004号。
Figure 1. Location of the study area
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 土地利用程度分级表
Table 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 表 2 2010—2020年林芝市各用地类型面积及其占比
Table 2 Area and proportion of land use types in Nyingchi City from 2010 to 2020
类型
land use types2010 2020 面积/km2
area比例/%
proportion面积/km2
area比例/%
proportion草地 grassland 29 898.98 26.20 29 867.39 26.17 耕地 cropland 1 257.37 1.10 1 332.12 1.17 建设用地 construction land 25.05 0.02 66.57 0.06 林地 forest 61 042.83 53.49 60 979.75 53.43 水域 water 3 258.32 2.86 3 253.11 2.85 未利用地 unused land 18 643.28 16.34 18 626.89 16.32 表 3 2010—2020年林芝市土地利用程度等级面积及其占比
Table 3 Area and proportion of land use intensity in Nyingchi City from 2010 to 2020
土地利用程度等级
grades of land use intensity2010 2020 面积/km2
area比例/%
proportion面积/km2
area比例/%
proportion弱 weakest 588.53 0.51 589.60 0.51 较弱 weaker 7 456.36 6.50 7 514.45 6.55 中等 middle 29 860.50 26.02 29 622.50 25.82 较强 stronger 55 626.20 48.48 52 337.40 45.61 强 strongest 21 208.00 18.48 24 675.60 21.51 -
[1] Global Land Project. Science plan and implementation strategy[R]. Stockholm: International Geosphere-Biosphere Program Secretariat, 2005.
[2] 匡文慧, 张树文, 杜国明, 等. 2015—2020年中国土地利用变化的时空遥感制图及时空特征分析[J]. 地理学报, 2022, 77(5): 1056. DOI: 10.11821/dlxb202205002. [3] 杨丽萍, 张静, 贡恩军, 等. GEE 联合多源数据的西安市土地利用时空格局及驱动力分析[J]. 农业工程学报, 2022, 38(2): 279. DOI: 10.11975/j.issn.1002-6819.2022.02.031. [4] 吴婷婷, 刘学录. 甘肃省庄浪县土地利用结构动态变化研究[J]. 云南农业大学学报(自然科学), 2015, 30(1): 112. DOI: 10.16211/j.issn.1004-390X(n).2015.01.019. [5] HUANG X, HUANG X J, LIU M M, et al. Spatial-temporal dynamics and driving forces of land development intensity in the western China from 2000 to 2015[J]. Chinese Geographical Science, 2020, 30(1): 16. DOI: 10.1007/s11769-020-1095-2.
[6] ERB K H, HABERL H, JEPSEN M R, et al. A conceptual framework for analysing and measuring land-use intensity[J]. Current Opinion in Environmental Sustainability, 2013, 5(5): 464. DOI: 10.1016/j.cosust.2013.07.010.
[7] 曹瑞芬, 蔡银莺. 基于AHP的土地开发利用程度评价及分析: 以武汉市为例[J]. 华中农业大学学报(社会科学版), 2011(1): 65. DOI: 10.13300/j.cnki.hnwkxb.2011.01.020. [8] 梁明, 聂拼, 陆胤昊, 等. 淮南市土地利用程度变化过程的时空演化特征[J]. 农业工程学报, 2019, 35(22): 99. DOI: 10.11975/j.issn.1002-6819.2019.22.011. [9] 吴金华, 李纪伟, 梁晶晶. 土地利用程度与效益关系研究: 以延安市为例[J]. 中国土地科学, 2011, 25(8): 54. DOI: 10.13708/j.cnki.cn11-2640.2011.08.010. [10] 楚玉山, 刘纪远. 西藏土地利用[M]. 北京: 科学出版社, 1992. [11] 黄木易, 何翔, 吴迪, 等. 巢湖流域土地利用程度变化及其空间异质性分析[J]. 土壤, 2015, 47(5): 994. DOI: 10.13758/j.cnki.tr.2015.05.029. [12] 张颖, 赵宇鸾. 黔桂岩溶山区土地利用程度演变的空间分异特征[J]. 水土保持研究, 2018, 25(1): 287. DOI: 10.13869/j.cnki.rswc.20171225.001. [13] 于明雪, 孙建国, 杨维涛, 等. 基于贝叶斯层次时空模型的甘肃省土地利用程度演变分析[J]. 地理科学, 2022, 42(5): 918. DOI: 10.13249/j.cnki.sgs.2022.05.017. [14] 吴琳娜, 杨胜天, 刘晓燕, 等. 1976年以来北洛河流域土地利用变化对人类活动程度的响应[J]. 地理学报, 2014, 69(1): 54. DOI: 10.11821/dlxb201401005. [15] 黄华富, 戴文远, 苏木兰. 海岛生态脆弱区土地利用程度空间格局演化: 以福建省海坛岛为例[J]. 福建师范大学学报(自然科学版), 2016, 32(2): 92. [16] HUANG W, WANG M Y, CHEN X Y. Analysis on spatial 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.
[17] 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.
[18] 罗芳, 潘安, 陈忠升, 等. 四川省土地利用变化对生态系统服务价值的影响研究[J]. 云南农业大学学报(自然科学), 2021, 36(4): 734. DOI: 10.12101/j.issn.1004-390X(n).202010031. [19] 肖建英, 谭术魁. 中国粮食产量省级尺度下的空间分异规律[J]. 中国土地科学, 2013, 27(8): 26. DOI: 10.13708/j.cnki.cn11-2640.2013.08.008. [20] 张松林, 张昆. 全局空间自相关Moran指数和G系数对比研究[J]. 中山大学学报(自然科学版), 2007, 46(4): 93. DOI: 10.3321/j.issn:0529-6579.2007.04.021. [21] 曹冯, 陈松林. 县域土地利用程度及其空间自相关探析: 以福建省德化县为例[J]. 福建师范大学学报(自然科学版), 2014, 30(3): 119. [22] 鱼坤, 陈学渊, 韩晓静, 等. 西北旱区土地利用变化的时空异质性研究: 以陕西铜川耀州区为例[J]. 西北农林科技大学学报(自然科学版), 2022, 50(9): 1. DOI: 10.13207/j.cnki.jnwafu.2022.09.011. [23] 杨玉婷, 石培基, 潘竟虎. 干旱内陆河流域土地利用程度差异分析: 以张掖市甘州区为例[J]. 干旱区资源与环境, 2012, 26(2): 102. DOI: 10.13448/j.cnki.jalre.2012.02.013. [24] 冯异星, 罗格平, 尹昌应, 等. 干旱区内陆河流域土地利用程度变化与生态安全评价: 以新疆玛纳斯河流域为例[J]. 自然资源学报, 2009, 24(11): 1921. DOI: 10.3321/j.issn:1000-3037.2009.11.007. [25] 冯舒芮, 黄琦, 杨京彪, 等. 三江源地区土地利用程度的时空变化分析[J]. 中央民族大学学报(自然科学版), 2020, 29(3): 21. DOI: 10.3969/j.issn.1005-8036.2020.03.004. [26] 林芝市地方志办公室. 林芝年鉴(2019)[M]. 拉萨: 西藏人民出版社, 2020. [27] 林振耀, 吴样定. 青藏高原气候区划[J]. 地理学报, 1981, 36(1): 22. DOI: 10.3321/j.issn:0375-5444.1981.01.003. [28] 林芝市统计局. 2020年林芝市国民经济和社会发展统计公报[EB/OL]. (2021-04-08)[2022-09-28]. http://www.linzhi.gov.cn/tjj/c103534/202105/73039b332f1045a6a50eecc141f64c80.shtml. [29] 庄大方, 刘纪远. 中国土地利用程度的区域分异模型研究[J]. 自然资源学报, 1997, 12(2): 105. DOI: 10.3321/j.issn:1000-3037.1997.02.002. [30] 邵俊明, 周宝同, 徐中强, 等. 基于Moran’s Ⅰ和Kriging插值在不同采样尺度下的土地利用程度分析: 以重庆市石柱县为例[J]. 四川农业大学学报, 2015, 33(2): 215. DOI: 10.16036/j.issn.1000-2650.2015.02.015. [31] 王劲峰, 廖一兰, 刘鑫. 空间数据分析教程[M]. 2版. 北京: 科学出版社, 2019. [32] 鲁安新, 姚檀栋, 王丽红, 等. 青藏高原典型冰川和湖泊变化遥感研究[J]. 冰川冻土, 2005, 27(6): 783. DOI: 10.3969/j.issn.1000-0240.2005.06.001. [33] 龙笛, 李雪莹, 李兴东, 等. 遥感反演2000—2020年青藏高原水储量变化及其驱动机制[J]. 水科学进展, 2022, 33(3): 375. DOI: 10.14042/j.cnki.32.1309.2022.03.003. -
期刊类型引用(5)
1. 魏伟,益西才措,尹力,薄立明. 西藏自治区国土空间格局演变特征及驱动力分析. 长江流域资源与环境. 2025(03): 558-572 . 百度学术
2. 李浩楠,窦蒙召,郭健斌. 林芝市代耕农空间分布的驱动机制研究. 农业与技术. 2024(16): 156-161 . 百度学术
3. 沈盈佳,李灿松,葛旭瑞,黄宇,刘宇斯,钱镜帆. 澜沧江—湄公河跨境流域土地利用变化及生态环境质量演变. 云南农业大学学报(自然科学). 2024(05): 178-190 . 百度学术
4. 宋玉林,杨泽凡,余敖,李浩. 广州增城区县域尺度的森林生物量碳储量估算及空间分布特征. 陆地生态系统与保护学报. 2024(04): 10-19 . 百度学术
5. 李国庆,赵洪升,张浩东. 淮北市林地景观格局特征与优化研究. 绿色科技. 2023(19): 96-100 . 百度学术
其他类型引用(3)