基于转移矩阵的重庆市北碚区土地利用变化分析 丁茜+付梅臣+王力摘要：获取土地利用/覆被变化定量信息对于土地管理和规划、促进可持续发展有重要意义以重庆市北碚区2000、2005、2010、2015年的landsat TM 影像为基本数据源，采用监督分类的方法获得了北碚区4个时相的土地利用数据，建立土地利用转移矩阵，从年均变化强度、地类动态度、地类转变强度及倾向度3个层次构建指标，深入挖掘土地利用变化信息，分析北碚区15年土地利用的变化速率、转变强度及变化趋势结果表明：2000—2015年重庆市北碚区的土地利用发生明显变化，变化强度呈递增趋势；各类用地在不同时段内变化分异较大，变化相对剧烈的地类主要是建设用地和耕地，耕地减少的强度与城市扩张速度呈正相关；地类综合动态度与地类占研究区面积的比例有关，单一动态度则反映地类自身的变化情况；地类转变强度和倾向度分别挖掘了地类变化的直观和潜在转变信息，二者结合更能全面地分析地类的来源与转向的情况关键词：土地利用变化；转移矩阵；遥感监测；重庆市北碚区： F301.24文献标志码： A：1002-1302（2017）12-0198-07通信作者：付梅臣，教授，主要从事土地利用和复垦研究。
E-mail：fumeichen@土地利用/覆被变化（LUCC）是全球气候变化和可持续发展研究领域关注的核心问题，也是人类活动与自然环境相互作用最直接的表现形式[1-2]研究区域系统内的土地利用结构变化规律，对土地利用调控、合理利用土地、促进区域社会经济发展具有重要意义[3]当前国内外土地利用变化研究主要涉及3个方面：一是土地利用演化规律研究[4-15]，包括数量和空间演变规律；二是引起土地利用变化的驱动机理分析[16-19]；三是基于变化过程构建有关模型模拟预测土地利用演化[3-4，20-24]土地利用动态过程监测是驱动力分析和模拟预测的基础，对深入开展土地利用变化研究具有重要意义土地利用转移矩阵反映了一个区域在2个时点间土地利用变化的数量结构特征与各用地类型变化的方向，因而在土地利用变化和模拟分析中具有重要意义，并得到广泛应用[6]传统的土地利用转移矩阵应用集中在直接分析土地利用类型面积的变化信息，其中又以土地利用类型面积的净变化、总变化以及动态度分析为主[7]，不能全面而深入挖掘土地利用变化的过程潜在信息近年来，重庆市北碚区社会经济快速增长，重大基础设施建设导致土地利用格局变化具有典型性，有重要研究价值。
本研究以重庆市北碚区为例，引入土地利用年均变化强度、地类动态度、地类转变强度及倾向度的分析计算方法，从3个层次对土地利用变化过程進行分析，并判断土地利用变化强度与倾向性的规则，以期为有效获取北碚区土地利用类型动态变化信息以及为土地利用空间分析和预测建模提供有效信息1研究区域与数据1.1研究区概况北碚区位于2944′N～2933N、10612E～10632E，地处重庆市核心区的西北面嘉陵江下游（图1）全区幅员面积75 500 hm2，常住人口77万人，其中城镇人口61万人，辖5个街道、12个镇北碚区位于四川盆地平行岭谷地区，属西南坳褶带，主要地形为山脉和丘陵，其间形成谷地，嘉陵江从西北向东南横流而过，支流20余条，水系丰富该地区属亚热带季风湿润气候，热量丰富，雨量充沛，最高气温44.3 ℃，最低气温-3.1 ℃，年平均气温18.2 ℃土壤以紫色土、水稻土为主，盛产水稻、玉米、红薯、水果等北碚区距重庆市中心24 km，距重庆市江北国际机场 27 km，是重庆市进出川北的咽喉要地全区公路里程 1 470 km，航道里程44.5 km，港口（码头）23个，襄渝铁路横穿东西，嘉陵江黄金水道纵贯南北，交通极为便捷。
按照重庆市五大功能区域战略布局，北碚区属都市功能拓展区北碚区立足都市功能拓展区定位，围绕“五区四带”发展战略部署，先后引进工业项目180余个，GDP能耗降至0.76 t标煤/万元2015年，完成地区生产总值430亿元，规模以上工业总产值841亿元，全社会固定资产投资595亿元，一般公共预算收入27.6亿元，增长12.1%，城乡常住居民人均可支配收入分别达到30 261、14 499元，扎实推进经济社会发展，进入了加速城市化发展的新阶段1.2数据来源及处理通过美国地质调查局网站，下载重庆市北碚区2000—2015年 landsat TM遥感影像，时段是5年，轨道号均为 128/039，分辨率为30 m30 m利用ENVI 5.1软件，对4个时段的影像进行辐射定标、波段合成、图像裁剪和掩膜等处理，根据研究区域实际情况，结合实地考察和Google Earth提供的高分辨率图像，通过选择样点，用监督分类的方法将研究区分成水域、建设用地、林地、耕地、其他用地5种用地类型，从而得到北碚区4期土地利用/覆被类型图（图2）经过混淆精度检验，4个时段影像的分类精度符合研究的精度要求（表1）利用ENVI对不同时相遥感图像的分类结果进行叠加运算，得到2000—2005年、2005—2010年、2010—2015年32方法2.1年均变化强度分析土地利用转变强度分析研究区每个时段的年均变化速式中：Ut为研究区在研究时段内的年均转变强度，%；Sn为研究区总的转入面积，km2；S为研究区域的土地总面积，km2；T为研究时段的年跨度，年。
2.2地类动态变化分析土地利用的显著特点之一就是其空间位置的固定性与特殊性如图3所示，实线部分表示地类i在T1时点的状态，虚线部分表示地类i在T2时点的状态，地类i从研究时点T1到T2的土地利用空间格局变化包含3种类型：（1）未变化部分（Sii），其土地利用类型与空间区位在时间段T1-T2没有发生变化；（2）转入部分（S+i- Sii），表示其他非i地类转变为地类i的部分；（3）转出部分（Si+- Sii），表示地类i在研究时间段转变成其他非i地类的部分endprint地类动态变化主要进行空间格局变化的转入和转出2个部分分析，研究方法分为2种：一种是地类综合动态度模型，该模型分析各地类相较于整个研究区域的变化程度；另一种是地类单一动态度模型，该模型分析不同地类相较于其自身发生变化之前的状态产生的变化强度Gsi=S+i- SiiS1T100%；（2）Lsi=Si+- SiiS1T100%；（3）Gi=S+i- SiiS+i1T100%；（4）Li=Si+- SiiS+i1T100%5）式中：Gsi和Lsi分別为地类i在综合动态度模型下某研究时段内的年均转入速率和年均转出速率；Gi和Li分别为地类i在单一动态度模型下某研究时段内的年均转入速率和年均转出速率；S为研究区的总面积；S+i和Si+分别为地类i在期初和期末的总面积；Sii为地类i没有发生变化的面积；T为研究时段的年跨度。
式（2）、（3）定义了地类综合动态模型的地类年均转入速率和转出速率，其值越大，表示地类相较于研究区变化速率相对越快；式（4）、（5）定义了单一动态模型的地类年均转入速率和转出速率，其值越大，表示地类相较于自身变化速率相对越快两者关系在于，地类单一动态度是地类综合动态度乘上该地类占研究区总面积的比值所得结果，消除了面积对变化强度的影响，对于分析特定时段内地类相较于自身面积变化速率快慢有着重要意义，两者结合可以较为全面地反映与比较各地类的变化程度2.3地类转变强度和倾向度分析地类转变强度研究针对某一特定地类，其转入和转出面积的主要来源以及转换的程度式（6）、（7）分别定义了地类年均转入强度和年均转出强度，其值越小，表示该地类的转变强度越小，且为非主要转入/转出来源；其值越大，表示该地类的转变强度越大，且为主要转入/转出来源[8]Gji=SjiSi+1T100%；（6）Lij=SijS+i1T100%7）式中：Gji为其他地类j转成地类i（i≠j）的转入强度；Lij为地类i转成其他地类j（i≠j）的转出强度；Sij表示某研究时段内地类i转成地类j的面积；Sji表示研究时段内地类j转成地类i的面积；T为研究时段的年跨度。
地类转变倾向度定义了不同土地利用类型之间转换的倾向性强度，该方法消除了土地利用类型面积对转换面积的影响，能够很好地反映地类之间相互转换的优势信息式（8）、（9）分别从地类转出量和转入量的角度，定义了在随机条件下土地利用类型之间转化的理论值计算方法[7]式（10）、（11）定义了地类实际转换与理论转换的差异强度，其值越小，说明该土地利用转换类型的倾向优势越强，相应地类转换系统性越强；其值越大，表明该土地利用转换类型的倾向优势越弱，相应地类转换随机性越强Wij=（Si+-Sii）Si+S-Si+；（8）Vji=（S+i-Sii）S+iS-S+i；（9）Rji=|Sij-Wij|Wij；（10）Qji=|Sji-Vji|Vji11）式中：Wij为地类i转成地类j（i≠j）的理论转出量；Vji为地类j转成地类i（i≠j）的理论转入量强度；Rij为地类i转成地类j（i≠j）的地类转出倾向度；Qji为地类j转成地类i（i≠j）的地类转入倾向度；Sij表示某研究时段内地类i转成地类j的面积；Sji表示研究时段内地类j转成地类i的面积3结果与分析3.1年均变化强度从地类年均变化强度（图4）来看，2000—2005年土地利用变化面积仅占研究区域的14%，年均变化强度为2.81%，变化相对缓慢；2005—2010年年均变化强度为3.50%，土地利用变化相对加快；2010—2015年土地利用变化最快，年均变化强度达5.32%。
从土地利用变化总体态势来看，2000—2015年北碚区整体的年均变化强度达到3.88%，从分时段的土地利用变化来看，土地利用年均变化强度呈逐渐增大趋势，且相较于整体变化而言，前2个时段变化相对缓慢；2010—2015年土地利用年均变化强度变化迅速，这与北碚区近年来都市功能拓展区的定位及“五区四带”发展战略的不断深入有着密不可分的联系3.2地类动态度从地类净变化来看（图5），耕地在3个时段内的动态度均是转出大于转入，说明耕地在研究时段内呈缩减趋势；建设用地在3个时段内的动态度均是转入大于转出，说明建设用地呈不断增长趋势；林地的净变化呈现先基本平稳到缩减再到增长的趋势；水域净变化呈现先减后增再减的趋势，但变化不明显；其他用地呈现先增后减再平稳的趋势从地类动态度来看（图5），耕地的综合转出动态度均最大，综合转入动态度除2010—2015年较小外，其他时间段变化都相对较大，说明15年间耕地面积变化很大，但是单一动态度均小于7%，这与耕地在研究区内所占比例最大有关，说明对于整个研究区域而言，耕地变化十分剧烈，但是就耕地自身面积来看，由于面积基数较大，变化的部分相对较少；其他用地与耕地变化趋势相反，由于其他用地面积在整个研究区的比例最小，因此综合动态度在3个时段内均最小，但是相较于自身面积来看，变化十分剧烈；建设用地的综合动态度和单一动态度呈增大趋势，且转入均大于转出，说明建设用地在3个时段内的变化强度不断增大，并不断扩张；由于林地面积占研究区比例较大，与其他4种地类相比，综合动态度较大，但单一动态度较小；水域的综合动态度和单一动态度均变化不大，说明水域在研究区内变化强度相对平稳，变化趋势不明显。
3.3地类转变强度和倾向度从转入强度和倾向度来看（图6），3个研究时段内林地一直是耕地增加部分的主要来源，且转入倾向度均小于1%，说明林地转为耕地为系统性转换；转入耕地来源其次是建设用地，转入比例有增大趋势，倾向度为3%～4%，变化不大，说明建设用地转为耕地属于随机性转换；水域和其他用地转入耕地的比例很小，且倾向度逐渐减小，说明2类用地转为耕地有随机到系统变化的趋势。